This study presents a numerical model for simulating dense interflows. The governing equations are provided and the finite difference
method is used with the k-ε turbulence model. The model is used to simulate a dense interflow established in a deep ambient water,
resulting velocity and excess density profiles. It is observed that velocity decreases in the longitudinal direction due to water entrainment
in the vicinity of the outlet and rarely changes for increased Richardson number. Similarity collapses of velocity and excess density are
obtained, but those of turbulent kinetic energy and dissipation rate are not. A shape factor for the dense interflow is obtained from the
simulated profiles. The value of this shape factor can be used in the layer-averaged modeling of dense interflows. In addition, a
buoyancy-related parameter ( ) for the k-ε model and the volume expansion coefficient () are obtained from the simulated results.