Numerical simulation of two coalescing turbulent forced plumes in linearly stratified fluids

A computational fluid dynamic model that can solve the Reynolds-averaged Navier-Stokes equations and the species transport equation is developed to simulate two coalescing turbulent forced plumes, which are released with initial momentum and buoyancy flux into a linearly stable stratified environment. The velocity fields, turbulence structures, and entrainment of two plumes with different source separations and source buoyancy fluxes are analyzed quantitatively, in comparison with a series of physical experiments. An empirical parameterization is proposed to predict the amplification of the maximum rise height of two coalescing forced plumes caused by superposition and mutual entrainment. The maximum values of both turbulent kinetic energy and turbulence dissipation rate decrease monotonically with the increase in source separation of the two turbulent plumes. However, the trajectory of the maximum turbulent viscosity attained in the plume cap region presents two notable enhancements. This variation may be attributed to the turbulence transported from the touching region and the strong mixing around the neutrally buoyant layer between two plumes, while the mixing is caused by the lateral convection and the rebound after overshooting. The plume entrainment coefficient in near vent stems has a positive relationship with the source Richardson number. A transition of flow regimes to plume-like flows would occur when the contribution of initial momentum is important. The entrainment coefficient will decrease in the touching region of two plumes due to mutual entrainment, while the superposition of plumes can lead to distortion of the boundary of plume sectors.