A conservative integration of the pseudo-incompressible equations with implicit turbulence parameterization

Durran's pseudo-incompressible equations are integrated in a mass and momentum conserving way with a new implicit turbulence model. This system is soundproof, which has two major advantages over fully compressible systems: the Courant-Friedrichs-Lewy (CFL) condition for stable time advancement is no longer dictated by the speed of sound and all waves in the model are clearly gravity waves (GW). Thus, the pseudo-incompressible equations are an ideal laboratory model for studying GW generation, propagation, and breaking. Gravity wave breaking creates turbulence that needs to be parameterized. For the first time the adaptive local deconvolution method (ALDM) for implicit large-eddy simulation (ILES) is applied to non- Boussinesq stratified flows. ALDM provides a turbulence model that is fully merged with the discretization of the flux function. In the context of non-Boussinesq stratified flows this poses some new numerical challenges-the solution of which is presented in this text. In numerical test cases the authors show the agreement of the results with the literature (Robert's hot-cold bubble test case), they present the sensitivity to the model's resolution and discretization, and they demonstrate qualitatively the behavior of the implicit turbulence model for a 2D breaking gravity wave packet.