A transient velocity difference model called the dynamic slip model has been integrated into the FLASH-4 numerical scheme. The numerical scheme resulting from combining this model with the general purpose RELAP4 network computer code is described. The result is a code that calculates slip in all junctions both vertical and horizontal. Results of computations are compared to U.S. Nuclear Regulatory Commission standard problems 1 and 6.