On computing the jump condition of the dissipation rate in the two-equation turbulence models for two-phase flow and application to air-water waves

Abstract: Traditional turbulence models are derived for single-phase flow. Extension of the family of two-equation turbulence models for two-phase flow is obtained via scaling the transport equations by the density. In the special case of two-phase flow with a sharp interface, jump conditions exist. Two types of jump conditions are found: (1) jump in the partial differential equation (PDE) physical quantities such as density and viscosity and (2) jump in the turbulence frequency. We first derive and clarify the jump in the equations. The jump in the turbulence frequency is proportional to the kinematic viscosity ratio, which is approximately $10$ in the case of air-water. Then a new field, the inverse turbulence area, is considered to model the turbulence effects instead of the turbulence frequency. For the system of air and water, the effect of the jump of the kinematic viscosity is always greater than the effect arising from the jump of velocity gradient. This approximation leads to the assumption of a continuous inverse turbulence area scale. Validation versus experimental measurements from the literature is then presented to demonstrate the improvement of the model. In particular, the wave breaking phenomenon is simulated in two conditions: spilling and plunging wave breakers. The proposed model shows its ability to predict the turbulence in the surf zone accurately. Finally, it explains the low values of the time-averaged turbulent kinetic energy in the surf zone which is caused by the increase of the turbulence frequency in the air.