Direct numerical simulation of a turbulent hydraulic jump: turbulence statistics and air entrainment

Abstract

We present direct numerical simulation (DNS) of a stationary turbulent hydraulic jump with inflow Froude number of 2, Weber number of 1820 and density ratio of 831, consistent with ambient water–air systems, all based on the inlet height and inlet velocity. A non-dissipative geometric volume of fluid (VOF) method is used to track the detailed interactions between turbulent flow structures and the nonlinear interface dynamics. Level set equations are also solved concurrent with VOF in order to calculate the interface curvature and surface tension forces. The mesh resolution is set to resolve a wide range of interfacial scales including the Hinze scale. Calculations are compared against experimental data of void fraction and interfacial scales indicating, reasonable agreement despite a Reynolds number mismatch. Multiple calculations are performed confirming weak sensitivity of low-order statistics and void fraction on the Reynolds number. The presented results provide, for the first time, a comprehensive quantitative data for a wide range of phenomena in a turbulent breaking wave using DNS. These include mean velocity fields, Reynolds stresses, turbulence production and dissipation, velocity spectra and air entrainment data. In addition, we present the energy budget as a function of streamwise location by keeping track of various energy exchange processes in the wake of the jump. The kinetic energy is mostly transferred to pressure work, potential energy and dissipation while surface energy plays a less significant role. Our results indicate that the rate associated with various energy exchange processes peak at different streamwise locations, with exchange to pressure work flux peaking first, followed by potential energy flux and then dissipation. The energy exchange process spans a streamwise length of order ${\sim}10$ jump heights. Furthermore, we report statistics associated with bubble transport downstream of the jump. The bubble formation is found to have a periodic nature. Meaning that the bubbles are generated in patches with a specific frequency associated with the roll-up frequency of the roller at the toe of the jump, with its footprint apparent in the velocity energy spectrum. Our study also provides the ensemble-averaged statistics of the flow which we present in this paper. These results are useful for the development and validation of reduced-order models such as dissipation models in wave dynamics simulations, Reynolds-averaged Navier–Stokes models and air entrainment models.

A two-phase mixing layer between parallel gas and liquid streams: multiphase turbulence statistics and influence of interfacial instability

“…Similar TKE equations have also been shown by Vallet et al (2001) and Mortazavi et al (2016). The terms on the right hand side are advection, pressure di↵usion, turbulent di↵usion, viscous di↵usion, dissipation, production, and surface-tension induced di↵usion, respectively.…”

Numerical Simulation of Hydraulic Jumps. Part 2: Recent Results and Future Outlook

“…The function f can be solved with different methods, and different schemes have been implemented to improve the accuracy of the free surface reconstruction. More complicated free surface tracking schemes are also possible; see for instance the methods used in the study of Mortazavi et al [19], Fuster et al [20], or the extensive review of Scardovelli and Zaleski [21]. It must be noted that the free surface is a movable interface where both density and viscosity change abruptly, being prone to significant numerical diffusion.…”

“…Mortazavi et al [19] recently presented the first Direct Numerical Simulation (DNS) of a stationary turbulent hydraulic jump with F 1 = 2 and a low Reynolds number (11,000) with no boundary layer development at the inlet. A non-dissipative geometric VOF method was used to track the strongly-distorted interface, while level-set equations were also solved in order to calculate the interface curvature for the surface tension forces' computation.…”

Wake behind a three-dimensional dry transom stern. Part 1. Flow structure and large-scale air entrainment

“…It is tempting to compare the average interface structure on the ship centerline of a dry transom stern to low Froude number two-dimensional hydraulic jumps (e.g. Mortazavi et al 2016;Lin et al 2012)) and two-dimensional sterns (Maki et al 2008;Rodriguez-Rodriguez et al 2011). However, the average interface of a two-dimensional hydraulic jump increases sharply over a distance twice the fluid initial depth h 1 and then asymptotically increases to an analytic value of ( 1 + q 1 + 8F r 2 h1 )/2 (Mortazavi et al 2016).…”