No Arabic abstract
We observe the emergence of strong vertical drafts in direct numerical simulations of the Boussinesq equations in a range of parameters of geophysical interest. These structures, which appear intermittently in space and time, generate turbulence and enhance kinetic and potential energy dissipation, providing an explanation for the observed variability of the local energy dissipation in the ocean and the modulation of its probability distribution function. We show how, due to the extreme drafts, in runs with Froude numbers observable in oceans, roughly $10%$ of the domain flow can account for up to $50%$ of the total volume dissipation, consistently with recent estimates based on oceanic models.
We investigate the large-scale intermittency of vertical velocity and temperature, and the mixing properties of stably stratified turbulent flows using both Lagrangian and Eulerian fields from direct numerical simulations, in a parameter space relevant for the atmosphere and the oceans. Over a range of Froude numbers of geophysical interest ($approx 0.05-0.3$) we observe very large fluctuations of the vertical components of the velocity and the potential temperature, localized in space and time, with a sharp transition leading to non-Gaussian wings of the probability distribution functions. This behavior is captured by a simple model representing the competition between gravity waves on a fast time-scale and nonlinear steepening on a slower time-scale. The existence of a resonant regime characterized by enhanced large-scale intermittency, as understood within the framework of the proposed model, is then linked to the emergence of structures in the velocity and potential temperature fields, localized overturning and mixing. Finally, in the same regime we observe a linear scaling of the mixing efficiency with the Froude number and an increase of its value of roughly one order of magnitude.
The large-scale structures in the ocean and the atmosphere are in geostrophic balance, and a conduit must be found to channel the energy to the small scales where it can be dissipated. In turbulence this takes the form of an energy cascade, whereas one possible mechanism in a balanced flow at large scales is through the formation of fronts, a common occurrence in geophysical dynamics. We show in this paper that an iconic configuration in laboratory and numerical experiments for the study of turbulence, that of the Taylor-Green or von Karman swirling flow, can be suitably adapted to the case of fluids with large aspect ratios, leading to the creation of an imposed large-scale vertical shear. To this effect we use direct numerical simulations of the Boussinesq equations without net rotation and with no small-scale modeling, and with this idealized Taylor-Green set-up. Various grid spacings are used, up to $2048^2times 256$ spatial points. The grids are always isotropic, with box aspect ratios of either $1:4$ or $1:8$. We find that when shear and stratification are comparable, the imposed shear layer resulting from the forcing leads to the formation of multiple fronts and filaments which destabilize and further evolve into a turbulent flow in the bulk, with a sizable amount of dissipation and mixing, and with a cycle of front creation, instability, and development of turbulence. The results depend on the vertical length scales for shear and for stratification, with stronger large-scale gradients being generated when the two length scales are comparable.
In this paper we advance physical background of the energy- and flux-budget turbulence closure based on the budget equations for the turbulent kinetic and potential energies and turbulent fluxes of momentum and buoyancy, and a new relaxation equation for the turbulent dissipation time-scale. The closure is designed for stratified geophysical flows from neutral to very stable and accounts for the Earth rotation. In accordance to modern experimental evidence, the closure implies maintaining of turbulence by the velocity shear at any gradient Richardson number Ri, and distinguishes between the two principally different regimes: strong turbulence at Ri << 1 typical of boundary-layer flows and characterised by the practically constant turbulent Prandtl number; and weak turbulence at Ri > 1 typical of the free atmosphere or deep ocean, where the turbulent Prandtl number asymptotically linearly increases with increasing Ri (which implies very strong suppression of the heat transfer compared to the momentum transfer). For use in different applications, the closure is formulated at different levels of complexity, from the local algebraic model relevant to the steady-state regime of turbulence to a hierarchy of non-local closures including simpler down-gradient models, presented in terms of the eddy-viscosity and eddy-conductivity, and general non-gradient model based on prognostic equations for all basic parameters of turbulence including turbulent fluxes.
We present an investigation of rapidly rotating (small Rossby number $Roll 1$) and stratified turbulence where the stratification strength is varied from weak (large Froude number $Frgg1$) to strong ($Frll1$). The investigation is set in the context of a reduced model derived from the Boussinesq equations that efficiently retains anisotropic inertia-gravity waves with order-one frequencies and highlights a regime of wave-eddy interactions. Numerical simulations of the reduced model are performed where energy is injected by a stochastic forcing of vertical velocity, which forces wave modes only. The simulations reveal two regimes characterized by the presence of well-formed, persistent and thin turbulent layers of locally-weakened stratification at small Froude numbers, and by the absence of layers at large Froude numbers. Both regimes are characterized by a large-scale barotropic dipole enclosed by small-scale turbulence. When the Reynolds number is not too large a direct cascade of barotropic kinetic energy is observed, leading to total energy equilibration. We examine net energy exchanges that occur through vortex stretching and vertical buoyancy flux and diagnose the horizontal scales active in these exchanges. We find that the baroclinic motions inject energy directly to the largest scales of the barotropic mode, implying that the large-scale barotropic dipole is not the end result of an inverse cascade within the barotropic mode.
We discuss a mean-field theory of generation of large-scale vorticity in a rotating density stratified developed turbulence with inhomogeneous kinetic helicity. We show that the large-scale nonuniform flow is produced due to ether a combined action of a density stratified rotating turbulence and uniform kinetic helicity or a combined effect of a rotating incompressible turbulence and inhomogeneous kinetic helicity. These effects result in the formation of a large-scale shear, and in turn its interaction with the small-scale turbulence causes an excitation of the large-scale instability (known as a vorticity dynamo) due to a combined effect of the large-scale shear and Reynolds stress-induced generation of the mean vorticity. The latter is due to the effect of large-scale shear on the Reynolds stress. A fast rotation suppresses this large-scale instability.