On the Two-point Correlation of Potential Vorticity in Rotating and Stratified Turbulence

A framework is developed to describe the two-point statistics of potential vorticity in rotating and stratified turbulence as described by the Boussinesq equations. The Karman-Howarth equation for the dynamics of the two-point correlation function of potential vorticity reveals the possibility of inertial-range dynamics in certain regimes in the Rossby, Froude, Prandtl and Reynolds number parameters. For the case of large Rossby and Froude numbers, and for the case of quasi-geostrophic dynamics, a linear scaling law with 2/3 prefactor is derived for the third-order mixed correlation between potential vorticity and velocity, a result that is analogous to the Kolmogorov 4/5-law for the third-order velocity structure function in turbulence theory.Comment: 10 pages, to appear in Journal of Fluid Mechanics (2006

Relating statistics to dynamics in axisymmetric homogeneous turbulence

The structure and the dynamics of homogeneous turbulence are modified by the presence of body forces such that the Coriolis or the buoyancy forces, which may render a wide range of turbulence scales anisotropic. The corresponding statistical characterization of such effects is done in physical space using structure functions, as well as in spectral space with spectra of two-point correlations, providing two complementary viewpoints. In this framework, second-order and third-order structure functions are put in parallel with spectra of two-point second- and third-order velocity correlation functions, using passage relations. Such relations apply in the isotropic case, or for isotropically averaged statistics, which, however, do not reflect the actual more complex structure of anisotropic turbulence submitted to rotation or stratification. This complexity is demonstrated in this paper by orientation-dependent energy and energy transfer spectra produced in both cases by means of a two-point statistical model for axisymmetric turbulence. We show that, to date, the anisotropic formalism used in the spectral transfer statistics is especially well-suited to analyze the refined dynamics of anisotropic homogeneous turbulence, and that it can help in the analysis of isotropically computed third-order structure function statistics often used to characterize anisotropic contexts.Comment: Physica

Energy- and flux-budget (EFB) turbulence closure model for the stably stratified flows. Part I: Steady-state, homogeneous regimes

We propose a new turbulence closure model based on the budget equations for the key second moments: turbulent kinetic and potential energies: TKE and TPE (comprising the turbulent total energy: TTE = TKE + TPE) and vertical turbulent fluxes of momentum and buoyancy (proportional to potential temperature). Besides the concept of TTE, we take into account the non-gradient correction to the traditional buoyancy flux formulation. The proposed model grants the existence of turbulence at any gradient Richardson number, Ri. Instead of its critical value separating - as usually assumed - the turbulent and the laminar regimes, it reveals a transition interval, 0.1< Ri <1, which separates two regimes of essentially different nature but both turbulent: strong turbulence at Ri<<1; and weak turbulence, capable of transporting momentum but much less efficient in transporting heat, at Ri>1. Predictions from this model are consistent with available data from atmospheric and lab experiments, direct numerical simulation (DNS) and large-eddy simulation (LES).Comment: 40 pages, 6 figures, Boundary-layer Meteorology, resubmitted, revised versio

Evidence for Bolgiano-Obukhov scaling in rotating stratified turbulence using high-resolution direct numerical simulations

We report results on rotating stratified turbulence in the absence of forcing, with large-scale isotropic initial conditions, using direct numerical simulations computed on grids of up to 4096^3 points. The Reynolds and Froude numbers are respectively equal to Re=5.4 x 10^4 and Fr=0.0242. The ratio of the Brunt-V\"ais\"al\"a to the inertial wave frequency, N/f, is taken to be equal to 4.95, a choice appropriate to model the dynamics of the southern abyssal ocean at mid latitudes. This gives a global buoyancy Reynolds number R_B=ReFr^2=32, a value sufficient for some isotropy to be recovered in the small scales beyond the Ozmidov scale, but still moderate enough that the intermediate scales where waves are prevalent are well resolved. We concentrate on the large-scale dynamics, for which we find a spectrum compatible with the Bolgiano-Obukhov scaling, and confirm that the Froude number based on a typical vertical length scale is of order unity, with strong gradients in the vertical. Two characteristic scales emerge from this computation, and are identified from sharp variations in the spectral distribution of either total energy or helicity. A spectral break is also observed at a scale at which the partition of energy between the kinetic and potential modes changes abruptly, and beyond which a Kolmogorov-like spectrum recovers. Large slanted layers are ubiquitous in the flow in the velocity and temperature fields, with local overturning events indicated by small Richardson numbers, and a small large-scale enhancement of energy directly attributable to the effect of rotation is also observed.Comment: 19 pages, 9 figures (including compound figures

Statistics of incompressible hydrodynamic turbulence: An alternative approach

Using a recent alternative form of the Kolmogorov-Monin exact relation for fully developed hydrodynamics (HD) turbulence, the incompressible energy cascade rate is computed. Under this current theoretical framework, for three-dimensional (3D) freely decaying homogeneous turbulence, the statistical properties of the fluid velocity (u), vorticity (ω= ×u), and Lamb vector (L=ω×u) are numerically studied. For different spatial resolutions, the numerical results show that can be obtained directly as the simple products of two-point increments of u and L, without the assumption of isotropy. Finally, the results for the largest spatial resolutions show a clear agreement with the cascade rates computed from the classical four-thirds law for isotropic homogeneous HD turbulence.Fil: Andrés, Nahuel. Universidad de Buenos Aires. Facultad de Cs.exactas y Naturales. Departamento de Física. Grupo de Plasmas Astrofisicos; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina. Laboratoire de Physique des Plasmas; FranciaFil: Banerjee, Supratik. Indian Institute Of Technology Kanpur; Indi

Generation of turbulence through frontogenesis in sheared stratified flows

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 K\'arm\'an 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^2\times 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.Comment: 19 pages, 15 figures, several simulations added in this new versio

New results on the model problem of the diffusion of turbulence from a plane source

The problem of the diffusion of turbulence from a plane source is addressed in the context of two-equation eddy-viscosity models and Reynolds-stress-transport models. In the steady state, full analytic solutions are given. At second order, they provide the equilibrium value of the anisotropy level obtained with different combinations of return-to-isotropy and turbulent-diffusion schemes and confirm the results obtained by Straatman et al. [AIAA J. 36, 929 (1998)] in an approximate analysis. In addition, all the characteristics of the turbulence decrease can be determined and it is shown that a special constraint on the value of the modeling constants should hold if turbulence fills the whole surrounding space. In a second step, precise results can be given for the unsteady model problem at the first-order-closure level. The evolution cannot be described with a single set of characteristic scales and one has to distinguish the cases of short and large times. In the short-time regime, the flow is governed by the characteristic scales of turbulence at the source and contamination of the flow proceeds as t^1/2. At large times, the flow is governed by time-dependent characteristic scales that correspond to the solution of the steady problem at the instantaneous location of the front. Contamination of the flow proceeds as a power of time that can be related to the value of the modeling constants. The role of a combination of these constants is emphasized whose value can be specified to produce a solution that matches simultaneously the experimental data for the decrease of turbulent kinetic energy in the steady state and the exponent of the propagation velocity in the transient regime

Anisotropy, inhomogeneity and inertial range scalings in turbulent convection

This paper provides a detailed study of scale-by-scale budgets in turbulent Rayleigh-B\'enard convection and aims at testing the applicability of Kolmogorov (1941) and Bolgiano (1959) theories for this flow. Particular emphasis is laid on anisotropic and inhomogeneous effects: the SO(3) decomposition of structure functions (Arad et al 1999) and a method of description of inhomogeneities proposed by Danaila et al (2001) are used to derive inhomogeneous and anisotropic generalizations of Kolmogorov and Yaglom equations applying to RB convection. The various terms in these equations are computed using data from a DNS of turbulent Boussinesq convection at \rayleigh=10^6 and \prandtl=1 with aspect ratio A=5. The analysis of the isotropic component demonstrates that the shape of the third-order velocity structure function is significantly influenced by buoyancy forcing and large-scale inhomogeneities, while the mixed third-order structure function appearing in Yaglom equation exhibits a clear scaling exponent 1 in a small range of scales. The magnitudes of the various low $\ell$ degree anisotropic components of the equations are also estimated and are shown to be comparable to their isotropic counterparts at moderate to large scales. Finally, a qualitative analysis shows that the influence of buoyancy forcing at scales smaller than the Bolgiano scale is likely to remain important up to \rayleigh=10^9, thus preventing Kolmogorov scalings from showing up in convective flows at lower Rayleigh numbers.Comment: 28 pages, 18 figures, accepted for publication in J. Fluid Mec

Turbulent diffusion in rapidly rotating flows with and without stable stratification

In this work, three different approaches are used for evaluating some Lagrangian properties of homogeneous turbulence containing anisotropy due to the application of a stable stratification and a solid-body rotation. The two external frequencies are the magnitude of the system vorticity, chosen vertical here, and the Brunt–Väisälä frequency, which gives the strength of the vertical stratification. Analytical results are derived using linear theory for the Eulerian velocity correlations (single-point, two-time) in the vertical and the horizontal directions, and Lagrangian ones are assumed to be equivalent, in agreement with an additional Corrsin assumption used by Kaneda (2000). They are compared with results from the kinematic simulation model (KS) by Nicolleau & Vassilicos (2000), which also incorporates the wave–vortex dynamics inherited from linear theory, and directly yields Lagrangian correlations as well as Eulerian ones. Finally, results from direct numerical simulations (DNS) are obtained and compared We address the question of the validity of Corrsin's simplified hypothesis, which states the equivalence between Eulerian and Lagrangian correlations. Vertical correlations are found to follow this postulate, but not the horizontal ones. Consequences for the vertical and horizontal one-particle dispersion are examined. In the analytical model, the squared excursion lengths are calculated by time integrating the Lagrangian (equal to the Eulerian) two-time correlations, according to Taylor's procedure. These quantities are directly computed from fluctuating trajectories by both KS and DNS. In the case of pure rotation, the analytical procedure allows us to relate Brownian asymptotic laws of dispersion in both the horizontal and vertical directions to the angular phase-mixing properties of the inertial waves. If stratification is present, the inertia–gravity wave dynamics, which affects the vertical motion, yields a suppressed vertical diffusivity, but not a suppressed horizontal diffusivity, since part of the horizontal velocity field escapes wavy motion

The spatio-temporal spectrum of turbulent flows

Identification and extraction of vortical structures and of waves in a disorganised flow is a mayor challenge in the study of turbulence. We present a study of the spatio-temporal behavior of turbulent flows in the presence of different restitutive forces. We show how to compute and analyse the spatio-temporal spectrum from data stemming from numerical simulations and from laboratory experiments. Four cases are considered: homogeneous and isotropic turbulence, rotating turbulence, stratified turbulence, and water wave turbulence. For homogeneous and isotropic turbulence, the spectrum allows identification of sweeping by the large scale flow. For rotating and for stratified turbulence, the spectrum allows identification of the waves, precise quantification of the energy in the waves and in the turbulent eddies, and identification of physical mechanisms such as Doppler shift and wave absorption in critical layers. Finally, in water wave turbulence the spectrum shows a transition from gravity-capillary waves to bound waves as the amplitude of the forcing is increased.Comment: Added new references and analysi