Stratified turbulence forced with columnar dipoles: numerical study

2015 ◽  
Vol 769 ◽  
pp. 403-443 ◽  
Author(s):  
Pierre Augier ◽  
Paul Billant ◽  
Jean-Marc Chomaz
Keyword(s):  
Reynolds Number ◽  
Energy Spectrum ◽  
Kinetic Energy ◽  
Numerical Simulations ◽  
Numerical Study ◽  
Stationary Flow ◽  
Length Scale ◽  
Computational Domain ◽  
Stratified Turbulence ◽  
Kinetic Energy Spectrum

This paper builds upon the investigation of Augier et al. (Phys. Fluids, vol. 26 (4), 2014) in which a strongly stratified turbulent-like flow was forced by 12 generators of vertical columnar dipoles. In experiments, measurements start to provide evidence of the existence of a strongly stratified inertial range that has been predicted for large turbulent buoyancy Reynolds numbers $\mathscr{R}_{t}={\it\varepsilon}_{\!K}/({\it\nu}N^{2})$, where ${\it\varepsilon}_{\!K}$ is the mean dissipation rate of kinetic energy, ${\it\nu}$ the viscosity and $N$ the Brunt–Väisälä frequency. However, because of experimental constraints, the buoyancy Reynolds number could not be increased to sufficiently large values so that the inertial strongly stratified turbulent range is only incipient. In order to extend the experimental results toward higher buoyancy Reynolds number, we have performed numerical simulations of forced stratified flows. To reproduce the experimental vortex generators, columnar dipoles are periodically produced in spatial space using impulsive horizontal body force at the peripheries of the computational domain. For moderate buoyancy Reynolds number, these numerical simulations are able to reproduce the results obtained in the experiments, validating this particular forcing. For higher buoyancy Reynolds number, the simulations show that the flow becomes turbulent as observed in Brethouwer et al. (J. Fluid Mech., vol. 585, 2007, pp. 343–368). However, the statistically stationary flow is horizontally inhomogeneous because the dipoles are destabilized quite rapidly after their generation. In order to produce horizontally homogeneous turbulence, high-resolution simulations at high buoyancy Reynolds number have been carried out with a slightly modified forcing in which dipoles are forced at random locations in the computational domain. The unidimensional horizontal spectra of kinetic and potential energies scale like $C_{1}{\it\varepsilon}_{\!K}^{2/3}k_{h}^{-5/3}$ and $C_{2}{\it\varepsilon}_{\!K}^{2/3}k_{h}^{-5/3}({\it\varepsilon}_{\!P}/{\it\varepsilon}_{\!K})$, respectively, with $C_{1}=C_{2}\simeq 0.5$ as obtained by Lindborg (J. Fluid Mech., vol. 550, 2006, pp. 207–242). However, there is a depletion in the horizontal kinetic energy spectrum for scales between the integral length scale and the buoyancy length scale and an anomalous energy excess around the buoyancy length scale probably due to direct transfers from large horizontal scale to small scales resulting from the shear and gravitational instabilities. The horizontal buoyancy flux co-spectrum increases abruptly at the buoyancy scale corroborating the presence of overturnings. Remarkably, the vertical kinetic energy spectrum exhibits a transition at the Ozmidov length scale from a steep spectrum scaling like $N^{2}k_{z}^{-3}$ at large scales to a spectrum scaling like $C_{K}{\it\varepsilon}_{\!K}^{2/3}k_{z}^{-5/3}$, with $C_{K}=1$, the classical Kolmogorov constant.

Sensitivity of stratified turbulence to the buoyancy Reynolds number

2013 ◽  
Vol 725 ◽  
pp. 1-22 ◽  
Author(s):  
P. Bartello ◽  
S. M. Tobias
Keyword(s):  
Reynolds Number ◽  
Energy Spectrum ◽  
Kinetic Energy ◽  
Numerical Simulations ◽  
Froude Number ◽  
Relative Size ◽  
Stratified Flow ◽  
Reynolds Numbers ◽  
Stratified Turbulence ◽  
Integral Scale

AbstractIn this article we present direct numerical simulations of stratified flow at resolutions of up to $204{8}^{2} \times 513$, to explore scalings for the dynamics of stably stratified turbulence. Recent work suggests that for strong enough stratification, the vertical integral scale of the turbulence adjusts to yield a vertical Froude number, ${F}_{v} $, of order unity at high enough Reynolds number, whilst the horizontal Froude number, ${F}_{h} $, decreases as stratification is increased. Our numerical simulations are consistent with predictions by Lindborg (J. Fluid Mech., vol. 550, 2006, pp, 207–242), and with numerical simulations at lower resolution, in that the horizontal kinetic energy spectrum follows a Kolmogorov spectrum (after replacing the wavenumber with the horizontal wavenumber) and that the horizontal potential energy spectrum similarly follows the Corrsin–Obukhov spectrum for a passive scalar. Most importantly, we build upon these previous results by thoroughly exploring the dependence of the horizontal spectrum of horizontal kinetic energy on both the stratification and the relative size of the vertical dissipation terms, as quantified by the buoyancy Reynolds number. Our most important result is that variations in the power-law exponent scale entirely with the buoyancy Reynolds number and not with the stratification itself, lending considerable support to the Lindborg (2006) hypothesis that horizontal spectra are independent of stratification at large Reynolds numbers. We further demonstrate that even at the large numerical resolution of this study, the spectrum and hence the dynamics are affected by the buoyancy Reynolds number unless it is larger than $O(10)$, indicating that extreme care must be taken when assessing claims made from previous numerical simulations of stratified flow at low or moderate resolution and extrapolating the results to geophysical or astrophysical Reynolds numbers.

The mixing layer and its coherence examined from the point of view of two-dimensional turbulence

1988 ◽  
Vol 192 ◽  
pp. 511-534 ◽  
Author(s):  
Marcel Lesieur ◽  
Chantal Staquet ◽  
Pascal Le Roy ◽  
Pierre Comte
Keyword(s):  
Energy Spectrum ◽  
Kinetic Energy ◽  
White Noise ◽  
Mixing Layer ◽  
Finite Size ◽  
Spanwise Direction ◽  
Computational Domain ◽  
Two Dimensional ◽  
Infinite Domain ◽  
Kinetic Energy Spectrum

A two-dimensional numerical large-eddy simulation of a temporal mixing layer submitted to a white-noise perturbation is performed. It is shown that the first pairing of vortices having the same sign is responsible for the formation of a continuous spatial longitudinal energy spectrum of slope between k−4 and k−3. After two successive pairings this spectral range extends to more than 1 decade. The vorticity thickness, averaged over several calculations differing by the initial white-noise realization, is shown to grow linearly, and eventually saturates. This saturation is associated with the finite size of the computational domain.We then examine the predictability of the mixing layer, considering the growth of decorrelation between pairs of flows differing slightly at the first roll-up. The inverse cascade of error through the kinetic energy spectrum is displayed. The error rate is shown to grow exponentially, and saturates together with the levelling-off of the vorticity thickness growth. Extrapolation of these results leads to the conclusion that, in an infinite domain, the two fields would become completely decorrelated. It turns out that the two-dimensional mixing layer is an example of flow that is unpredictable and possesses a broadband kinetic energy spectrum, though composed mainly of spatially coherent structures.It is finally shown how this two-dimensional predictability analysis can be associated with the growth of a particular spanwise perturbation developing on a Kelvin-Helmholtz billow: this is done in the framework of a one-mode spectral truncation in the spanwise direction. Within this analogy, the loss of two-dimensional predictability would correspond to a return to three-dimensionality and a loss of coherence. We indicate also how a new coherent structure could then be recreated, using an eddy-viscosity assumption and the linear instability of the mean inflexional shear.

scholarly journals Indications of Stratified Turbulence in a Mechanistic GCM

2013 ◽  
Vol 70 (1) ◽  
pp. 231-247 ◽  
Author(s):  
Sebastian Brune ◽  
Erich Becker
Keyword(s):  
Energy Spectrum ◽  
Kinetic Energy ◽  
General Circulation ◽  
Dynamic Instability ◽  
Circulation Model ◽  
Vertical Resolution ◽  
Stratified Turbulence ◽  
Kinetic Energy Spectrum ◽  
Spectral Flux ◽  
Diffusion Scheme

Abstract The horizontal kinetic energy spectrum and its budget are analyzed on the basis of a general circulation model with simplistic parameterizations of radiative and latent heating. A spectral truncation at total wavenumber 330 is combined with a level spacing of either ~200 m or ~1.5 km from the midtroposphere to the lower stratosphere. The subgrid-scale parameterization consists of a Smagorinsky-type anisotropic diffusion scheme that is scaled by a Richardson criterion for dynamic instability and combined with a stress-tensor-based hyperdiffusion that acts only on the very smallest resolved scales. Simulations with both vertical resolutions show a transition from the synoptic −3 to the mesoscale slope in the upper-tropospheric kinetic energy spectrum. Analysis of the spectral budget indicates that the mesoscale slope can be interpreted as stratified turbulence, as has been proposed in recent works of Lindborg and others, only when a high vertical resolution is applied. In this case, the mesoscale kinetic energy around 300–150 hPa is dominated by the nonrotational flow, and the forward horizontal energy cascade is accompanied by an equally strong forward spectral flux due to adiabatic conversion. This adiabatic conversion mainly results from a vertical potential energy flux that originates in the midtroposphere, where the mesoscale adiabatic conversion is negative. For a conventionally coarse vertical resolution, however, the mesoscale slope in the troposphere is dominated by the rotational flow, and the spectral flux due to adiabatic conversion is not comparable to that associated with horizontal advection.

scholarly journals Scaling analysis and simulation of strongly stratified turbulent flows

2007 ◽  
Vol 585 ◽  
pp. 343-368 ◽  
Author(s):  
G. BRETHOUWER ◽  
P. BILLANT ◽  
E. LINDBORG ◽  
J.-M. CHOMAZ
Keyword(s):  
Reynolds Number ◽  
Numerical Simulations ◽  
Turbulent Flows ◽  
Scaling Laws ◽  
Three Dimensional ◽  
Direct Numerical Simulations ◽  
Small Scale ◽  
Scaling Analysis ◽  
Length Scale ◽  
Stratified Turbulence

Direct numerical simulations of stably and strongly stratified turbulent flows with Reynolds number Re ≫ 1 and horizontal Froude number Fh ≪ 1 are presented. The results are interpreted on the basis of a scaling analysis of the governing equations. The analysis suggests that there are two different strongly stratified regimes according to the parameter $\mathcal{R} \,{=}\, \hbox{\it Re} F^2_h$. When $\mathcal{R} \,{\gg}\, 1$, viscous forces are unimportant and lv scales as lv ∼ U/N (U is a characteristic horizontal velocity and N is the Brunt–Väisälä frequency) so that the dynamics of the flow is inherently three-dimensional but strongly anisotropic. When $\mathcal{R} \,{\ll}\, 1$, vertical viscous shearing is important so that $l_v \,{\sim}\, l_h/\hbox{\it Re}^{1/2}$ (lh is a characteristic horizontal length scale). The parameter $\cal R$ is further shown to be related to the buoyancy Reynolds number and proportional to (lO/η)4/3, where lO is the Ozmidov length scale and η the Kolmogorov length scale. This implies that there are simultaneously two distinct ranges in strongly stratified turbulence when $\mathcal{R} \,{\gg}\, 1$: the scales larger than lO are strongly influenced by the stratification while those between lO and η are weakly affected by stratification. The direct numerical simulations with forced large-scale horizontal two-dimensional motions and uniform stratification cover a wide Re and Fh range and support the main parameter controlling strongly stratified turbulence being $\cal R$. The numerical results are in good agreement with the scaling laws for the vertical length scale. Thin horizontal layers are observed independently of the value of $\cal R$ but they tend to be smooth for $\cal R$< 1, while for $\cal R$ > 1 small-scale three-dimensional turbulent disturbances are increasingly superimposed. The dissipation of kinetic energy is mostly due to vertical shearing for $\cal R$ < 1 but tends to isotropy as $\cal R$ increases above unity. When $\mathcal{R}$ < 1, the horizontal and vertical energy spectra are very steep while, when $\cal R$ > 1, the horizontal spectra of kinetic and potential energy exhibit an approximate k−5/3h-power-law range and a clear forward energy cascade is observed.

Vortex-Induced Vibrations of a Cylinder at Subcritical Reynolds Number

2011 ◽  
Author(s):  
Yoann Jus ◽  
Elisabeth Longatte ◽  
Jean-Camille Chassaing ◽  
Pierre Sagaut
Keyword(s):  
Reynolds Number ◽  
Numerical Simulations ◽  
Numerical Study ◽  
Damping Ratio ◽  
Experimental Studies ◽  
Cross Flow ◽  
Physical Analysis ◽  
Arbitrary Lagrangian Eulerian ◽  
Vortex Induced Vibrations ◽  
Low Mass

The present work focusses on the numerical study of Vortex-Induced Vibrations (VIV) of an elastically mounted cylinder in a cross flow at moderate Reynolds numbers. Low mass-damping experimental studies show that the dynamic behavior of the cylinder exhibits a three-branch response model, depending on the range of the reduced velocity. However, few numerical simulations deal with accurate computations of the VIV amplitudes at the lock-in upper branch of the bifurcation diagram. In this work, the dynamic response of the cylinder is investigated by means of three-dimensional Large Eddy Simulation (LES). An Arbitrary Lagrangian Eulerian framework is employed to account for fluid solid interface boundary motion and grid deformation. Numerous numerical simulations are performed at a Reynolds number of 3900 for both no damping and low-mass damping ratio and various reduced velocities. A detailed physical analysis is conducted to show how the present methodology is able to capture the different VIV responses.

scholarly journals Waves and vortices in the inverse cascade regime of stratified turbulence with or without rotation

2016 ◽  
Vol 806 ◽  
pp. 165-204 ◽  
Author(s):  
Corentin Herbert ◽  
Raffaele Marino ◽  
Duane Rosenberg ◽  
Annick Pouquet
Keyword(s):  
Reynolds Number ◽  
Numerical Simulations ◽  
Direct Numerical Simulations ◽  
Vertical Shear ◽  
Energy Spectra ◽  
Stratified Turbulence ◽  
Viscous Effects ◽  
Inverse Cascade ◽  
Main Effects ◽  
The Waves

We study the partition of energy between waves and vortices in stratified turbulence, with or without rotation, for a variety of parameters, focusing on the behaviour of the waves and vortices in the inverse cascade of energy towards the large scales. To this end, we use direct numerical simulations in a cubic box at a Reynolds number $Re\approx 1000$, with the ratio between the Brunt–Väisälä frequency $N$ and the inertial frequency $f$ varying from $1/4$ to 20, together with a purely stratified run. The Froude number, measuring the strength of the stratification, varies within the range $0.02\leqslant Fr\leqslant 0.32$. We find that the inverse cascade is dominated by the slow quasi-geostrophic modes. Their energy spectra and fluxes exhibit characteristics of an inverse cascade, even though their energy is not conserved. Surprisingly, the slow vortices still dominate when the ratio $N/f$ increases, also in the stratified case, although less and less so. However, when $N/f$ increases, the inverse cascade of the slow modes becomes weaker and weaker, and it vanishes in the purely stratified case. We discuss how the disappearance of the inverse cascade of energy with increasing $N/f$ can be interpreted in terms of the waves and vortices, and identify the main effects that can explain this transition based on both inviscid invariants arguments and viscous effects due to vertical shear.

scholarly journals Calculating Energy Spectra from Drifters

2016 ◽  
Author(s):  
Joseph H. LaCasce
Keyword(s):  
Energy Spectrum ◽  
Kinetic Energy ◽  
Structure Function ◽  
Two Dimensions ◽  
Longitudinal Structure ◽  
Kinetic Energy Spectrum ◽  
Longitudinal Structure Function ◽  
Large Pair ◽  
Non Local ◽  
Lagrangian Data

The relations between the kinetic energy spectrum and the second order longitudinal structure function in two dimensions are derived, and several examples are considered. The forward conversion (from spectrum to structure function) is illustrated first with idealized power law spectra, representing turbulent inertial ranges. The forward conversion is also applied to the zonal kinetic energy spectrum of Nastrom and Gage (1985) and the result agrees well with the longitudinal structure function of Lindborg (1999). The inverse conversion (from structure function to spectrum) is tested with data from 2D turbulence simulations. When applied to the theoretical structure function (derived from the forward conversion of the spectrum), the result closely resembles the original spectrum, except at the largest wavenumbers. However the inverse conversion is much less successful when applied to the structure function obtained from pairs of particles in the flow. This is because the inverse conversion favors large pair separations, which are typically noisy with particle data. Fitting the structure function to a polynomial improves the result, but not sufficiently to distinguish the correct inertial range dependencies. Furthermore the inversion of non-local spectra is largely unsuccessful. Thus it appears that focusing on structure functions with Lagrangian data is preferable to estimating spectra.

Sign in / Sign up

Export Citation Format

Share Document