No abstract

Results of direct numerical simulations of the transitional processes that characterise the evolution of a breaking internal gravity wave to a fully developed and essentially steady turbulent patch are presented. The stationary lee wave was forced by the imposition of an appropriate bottom boundary shape within a density-stratified domain having a uniform upstream velocity and density gradient, and with the ratio of momentum to thermal (or other) diffusivity defined by Pr = 1. An earlier paper considered the eventual, fully developed turbulent patch arising after the breaking process is complete (Yakovenko et al., J. Fluid Mech., vol. 677, 2011, pp. 103-133); the focus in this paper is on the instabilities in the breaking process itself. The flow is analysed using streamlines, density contours and temporal and spatial spectra, as well as second moments of the velocity and density fluctuations, for a Reynolds number of 4000 based on the height of the bottom topography and the upstream velocity. The computations (on a grid using in excess of 10 9 mesh points) yielded sufficient resolution to capture the fine-scale transition processes as well as the subsequent fully developed turbulence discussed earlier. It is shown that the major instability is of Rayleigh-Taylor type (RTI) with a resulting mixing region depth growing in a manner consistent with more classical RTI studies, despite the much more complicated environment. The resolution was sufficient to capture secondary Kelvin-Helmholtz-type instabilities on the developing RTI structures. Overall evolution towards the fully turbulent state characterised by a significant region of − 5 3 subrange in both velocity and density spectra is very rapid. It is much faster than the long time scale characterising the subsequent evolution of the turbulent patch; this latter time scale is sufficiently large that the turbulent patch can itself be viewed as essentially steady.

Results of direct numerical simulations of the development of a breaking internal gravity wave are presented. The wave was forced by the imposition of an appropriate bottom boundary shape (a two-dimensional cosine hill) within a density-stratified domain having a uniform upstream velocity and density gradient. The focus is on turbulence generation and maintenance within the turbulent patch generated by the wave breaking. Pathlines, density contours, temporal and spatial spectra, and second moments of the velocity and density fluctuations and turbulent kinetic energy balance terms obtained from the data averaged over the span in the mixed zone are all used in the analysis of the flow. Typical Reynolds numbers, based on the vertical scale of the breaking region and the upstream velocity, were around 6000 and the fully resolved computations yielded sufficient resolution to capture the fine-scale transition processes as well as the subsequent fully developed turbulence. It is shown that globally, within the turbulent patch, there is an approximate balance in the production, dissipation and transport processes for turbulent kinetic energy, so that the patch remains quasi-steady over a significant time. Although it is far from being axially homogeneous, with turbulence generation occurring largely near the upstream bottom part of the patch where the mean velocity shear is particularly large, it has features not dissimilar to those of a classical turbulent wake.