scholarly journals Spectral analysis of the transition to turbulence from a dipole in stratified fluid

2012 ◽  
Vol 713 ◽  
pp. 86-108 ◽  
Author(s):  
Pierre Augier ◽  
Jean-Marc Chomaz ◽  
Paul Billant
Keyword(s):  
Reynolds Number ◽  
Spectral Analysis ◽  
Kinetic Energy ◽  
Numerical Simulations ◽  
Scaling Laws ◽  
Stratified Fluid ◽  
Shear Instability ◽  
Transition To Turbulence ◽  
Small Scale ◽  
Wide Range

AbstractWe investigate the spectral properties of the turbulence generated during the nonlinear evolution of a Lamb–Chaplygin dipole in a stratified fluid for a high Reynolds number $Re= 28\hspace{0.167em} 000$ and a wide range of horizontal Froude number ${F}_{h} \in [0. 0225~0. 135] $ and buoyancy Reynolds number $\mathscr{R}= Re{{F}_{h} }^{2} \in [14~510] $. The numerical simulations use a weak hyperviscosity and are therefore almost direct numerical simulations (DNS). After the nonlinear development of the zigzag instability, both shear and gravitational instabilities develop and lead to a transition to small scales. A spectral analysis shows that this transition is dominated by two kinds of transfer: first, the shear instability induces a direct non-local transfer toward horizontal wavelengths of the order of the buoyancy scale ${L}_{b} = U/ N$, where $U$ is the characteristic horizontal velocity of the dipole and $N$ the Brunt–Väisälä frequency; second, the destabilization of the Kelvin–Helmholtz billows and the gravitational instability lead to small-scale weakly stratified turbulence. The horizontal spectrum of kinetic energy exhibits a ${{\varepsilon }_{K} }^{2/ 3} { k}_{h}^{\ensuremath{-} 5/ 3} $ power law (where ${k}_{h} $ is the horizontal wavenumber and ${\varepsilon }_{K} $ is the dissipation rate of kinetic energy) from ${k}_{b} = 2\lrm{\pi} / {L}_{b} $ to the dissipative scales, with an energy deficit between the integral scale and ${k}_{b} $ and an excess around ${k}_{b} $. The vertical spectrum of kinetic energy can be expressed as $E({k}_{z} )= {C}_{N} {N}^{2} { k}_{z}^{\ensuremath{-} 3} + C{{\varepsilon }_{K} }^{2/ 3} { k}_{z}^{\ensuremath{-} 5/ 3} $ where ${C}_{N} $ and $C$ are two constants of order unity and ${k}_{z} $ is the vertical wavenumber. It is therefore very steep near the buoyancy scale with an ${N}^{2} { k}_{z}^{\ensuremath{-} 3} $ shape and approaches the ${{\varepsilon }_{K} }^{2/ 3} { k}_{z}^{\ensuremath{-} 5/ 3} $ spectrum for ${k}_{z} \gt {k}_{o} $, ${k}_{o} $ being the Ozmidov wavenumber, which is the cross-over between the two scaling laws. A decomposition of the vertical spectra depending on the horizontal wavenumber value shows that the ${N}^{2} { k}_{z}^{\ensuremath{-} 3} $ spectrum is associated with large horizontal scales $\vert {\mathbi{k}}_{h} \vert \lt {k}_{b} $ and the ${{\varepsilon }_{K} }^{2/ 3} { k}_{z}^{\ensuremath{-} 5/ 3} $ spectrum with the scales $\vert {\mathbi{k}}_{h} \vert \gt {k}_{b} $.

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.

scholarly journals Onset of secondary instabilities on the zigzag instability in stratified fluids

2011 ◽  
Vol 682 ◽  
pp. 120-131 ◽  
Author(s):  
PIERRE AUGIER ◽  
PAUL BILLANT
Keyword(s):  
Numerical Simulations ◽  
Gravitational Instability ◽  
Dominant Role ◽  
Vortex Pair ◽  
Secondary Instability ◽  
Shear Instability ◽  
Small Scale ◽  
Scaling Analysis ◽  
Wide Range ◽  
Secondary Instabilities

Recently, Deloncle, Billant & Chomaz (J. Fluid Mech., vol. 599, 2008, p. 229) and Waite & Smolarkiewicz (J. Fluid Mech., vol. 606, 2008, p. 239) have performed numerical simulations of the nonlinear evolution of the zigzag instability of a pair of counter-rotating vertical vortices in a stratified fluid. Both studies report the development of a small-scale secondary instability when the vortices are strongly bent if the Reynolds number Re is sufficiently high. However, the two papers are at variance about the nature of this secondary instability: it is a shear instability according to Deloncle et al. (J. Fluid Mech., vol. 599, 2008, p. 229) and a gravitational instability according to Waite & Smolarkiewicz (J. Fluid Mech., vol. 606, 2008, p. 239). They also profoundly disagree about the condition for the onset of the secondary instability: ReF2h > O(1) according to the former or ReFh > 80 according to the latter, where Fh is the horizontal Froude number. In order to understand the origin of these discrepancies, we have carried out direct numerical simulations of the zigzag instability of a Lamb–Chaplygin vortex pair for a wide range of Reynolds and Froude numbers. The threshold for the onset of a secondary instability is found to be ReF2h ≃ 4 for Re ≳ 3000 and ReFh = 80 for Re ≲ 1000 in agreement with both previous studies. We show that the scaling analysis of Deloncle et al. (J. Fluid Mech., vol. 599, 2008, p. 229) can be refined to obtain a universal threshold: (Re − Re0)F2h ≃ 4, with Re0 ≃ 400, which works for all Re. Two different regimes for the secondary instabilities are observed: when (Re − Re0)F2h ≃ 4, only the shear instability develops while when (Re − Re0)F2h ≫ 4, both shear and gravitational instabilities appear almost simultaneously in distinct regions of the vortices. However, the shear instability seems to play a dominant role in the breakdown into small scales in the range of parameters investigated.

scholarly journals Surface kinetic energy transfer in surface quasi-geostrophic flows

2008 ◽  
Vol 604 ◽  
pp. 165-174 ◽  
Author(s):  
XAVIER CAPET ◽  
PATRICE KLEIN ◽  
BACH LIEN HUA ◽  
GUILLAUME LAPEYRE ◽  
JAMES C. MCWILLIAMS
Keyword(s):  
Kinetic Energy ◽  
Surface Velocity ◽  
Small Scale ◽  
Surface Dynamics ◽  
Spectral Flux ◽  
Advection Term ◽  
Wide Range ◽  
Geostrophic Flows ◽  
Scale Range ◽  
Density Variance

The relevance of surface quasi-geostrophic dynamics (SQG) to the upper ocean and the atmospheric tropopause has been recently demonstrated in a wide range of conditions. Within this context, the properties of SQG in terms of kinetic energy (KE) transfers at the surface are revisited and further explored. Two well-known and important properties of SQG characterize the surface dynamics: (i) the identity between surface velocity and density spectra (when appropriately scaled) and (ii) the existence of a forward cascade for surface density variance. Here we show numerically and analytically that (i) and (ii) do not imply a forward cascade of surface KE (through the advection term in the KE budget). On the contrary, advection by the geostrophic flow primarily induces an inverse cascade of surface KE on a large range of scales. This spectral flux is locally compensated by a KE source that is related to surface frontogenesis. The subsequent spectral budget resembles those exhibited by more complex systems (primitive equations or Boussinesq models) and observations, which strengthens the relevance of SQG for the description of ocean/atmosphere dynamics near vertical boundaries. The main weakness of SQG however is in the small-scale range (scales smaller than 20–30 km in the ocean) where it poorly represents the forward KE cascade observed in non-QG numerical simulations.

Effect of Size and Slip Coefficient of a Porous Manifold on Thermal Stratification in Storage Tanks

2009 ◽  
Author(s):  
N. M. Brown ◽  
F. C. Lai
Keyword(s):  
Experimental Data ◽  
Reynolds Number ◽  
Numerical Simulations ◽  
Storage Tank ◽  
Richardson Number ◽  
Thermal Stratification ◽  
Temperature Fields ◽  
Storage Tanks ◽  
Slip Coefficient ◽  
Wide Range

Numerical simulations have been performed to study the effects of size and slip coefficient of a porous manifold on the thermal stratification in a storage tank. The model is used to predict the development of flow and temperature fields during a charging process. Computations have covered a wide range of the Grashof number (1.8 × 105 &lt; Gr &lt; 1.8 × 108) and Reynolds number (10 ≤ Re ≤ 104), or in terms of the Richardson number, 10−2 &lt; Ri &lt; 105. The results obtained compare favorably well with the experimental data. In addition, the present results have confirmed the effectiveness of porous manifold in the promotion of thermal stratification and provide useful information for the design of such system.

Fluid Flow and Lateral Fluid Dispersion in Bounded Granular Beds

2002 ◽  
Author(s):  
Azita Soleymani ◽  
Eveliina Takasuo ◽  
Piroz Zamankhan ◽  
William Polashenski
Keyword(s):  
Reynolds Number ◽  
Porous Media ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Numerical Study ◽  
Fluid Velocity ◽  
Three Dimensional ◽  
Numerical Results ◽  
Transition To Turbulence ◽  
Wide Range

Results are presented from a numerical study examining the flow of a viscous, incompressible fluid through random packing of nonoverlapping spheres at moderate Reynolds numbers (based on pore permeability and interstitial fluid velocity), spanning a wide range of flow conditions for porous media. By using a laminar model including inertial terms and assuming rough walls, numerical solutions of the Navier-Stokes equations in three-dimensional porous packed beds resulted in dimensionless pressure drops in excellent agreement with those reported in a previous study (Fand et al., 1987). This observation suggests that no transition to turbulence could occur in the range of Reynolds number studied. For flows in the Forchheimer regime, numerical results are presented of the lateral dispersivity of solute continuously injected into a three-dimensional bounded granular bed at moderate Peclet numbers. Lateral fluid dispersion coefficients are calculated by comparing the concentration profiles obtained from numerical and analytical methods. Comparing the present numerical results with data available in the literature, no evidence has been found to support the speculations by others for a transition from laminar to turbulent regimes in porous media at a critical Reynolds number.

scholarly journals An Experimental Study of Water Flow in Smooth and Rough Rectangular Micro-Channels

2004 ◽  
Author(s):  
R. Baviere ◽  
F. Ayela ◽  
S. Le Person ◽  
M. Favre-Marinet
Keyword(s):  
Reynolds Number ◽  
Water Flow ◽  
Stokes Equations ◽  
Transition To Turbulence ◽  
Friction Coefficients ◽  
Micro Channels ◽  
Wide Range ◽  
Macro Scale ◽  
Experimental Apparatus ◽  
Poiseuille Number

This paper presents experimental results concerning water flow in smooth and rough rectangular micro-channels. It is part of a work intended to test the classical fluid mechanics laws when the characteristic length scale of inner liquid flows falls below 500μm. The method consists in determining experimental friction coefficients as a function of the Reynolds number. This implies simultaneous measurements of pressure drop and flow rates in microstructures. The two experimental apparatus used in this study enabled us to explore a wide range of length scales (7μm to 300μm) and of Reynolds number (0.01 to 8,000). Classical machining technologies were used to make micro-channels of various heights down to a scale of 100μm. Smaller silicon-Pyrex micro-channels were also made by means of silicon-based micro technologies. In these structures, friction coefficients have been measured locally with Cu-Ni strain gauges. For every height tested, both smooth and rough walls were successively used. When compared to macro-scale correlation the results demonstrate that i) In the smooth case, friction is correctly predicted by the Navier-Stokes equations with the classical kinematic boundary conditions, ii) For 200μm high channels, visualizations show transition to turbulence at Reynolds number of about 3,000. The presence of roughness elements did not significantly influence this result and iii) Roughness considerably increases the friction coefficient in the laminar regime. However, the Poiseuille number remains independent of the Reynolds number.

Hot-wire spatial resolution issues in wall-bounded turbulence

2009 ◽  
Vol 635 ◽  
pp. 103-136 ◽  
Author(s):  
N. HUTCHINS ◽  
T. B. NICKELS ◽  
I. MARUSIC ◽  
M. S. CHONG
Keyword(s):  
Reynolds Number ◽  
Boundary Layers ◽  
Spatial Resolution ◽  
Reynolds Numbers ◽  
Energy Spectra ◽  
Small Scale ◽  
Wire Length ◽  
Hot Wire ◽  
Wide Range ◽  
Near Wall

Careful reassessment of new and pre-existing data shows that recorded scatter in the hot-wire-measured near-wall peak in viscous-scaled streamwise turbulence intensity is due in large part to the simultaneous competing effects of the Reynolds number and viscous-scaled wire length l+. An empirical expression is given to account for these effects. These competing factors can explain much of the disparity in existing literature, in particular explaining how previous studies have incorrectly concluded that the inner-scaled near-wall peak is independent of the Reynolds number. We also investigate the appearance of the so-called outer peak in the broadband streamwise intensity, found by some researchers to occur within the log region of high-Reynolds-number boundary layers. We show that the ‘outer peak’ is consistent with the attenuation of small scales due to large l+. For turbulent boundary layers, in the absence of spatial resolution problems, there is no outer peak up to the Reynolds numbers investigated here (Reτ = 18830). Beyond these Reynolds numbers – and for internal geometries – the existence of such peaks remains open to debate. Fully mapped energy spectra, obtained with a range of l+, are used to demonstrate this phenomenon. We also establish the basis for a ‘maximum flow frequency’, a minimum time scale that the full experimental system must be capable of resolving, in order to ensure that the energetic scales are not attenuated. It is shown that where this criterion is not met (in this instance due to insufficient anemometer/probe response), an outer peak can be reproduced in the streamwise intensity even in the absence of spatial resolution problems. It is also shown that attenuation due to wire length can erode the region of the streamwise energy spectra in which we would normally expect to see kx−1 scaling. In doing so, we are able to rationalize much of the disparity in pre-existing literature over the kx−1 region of self-similarity. Not surprisingly, the attenuated spectra also indicate that Kolmogorov-scaled spectra are subject to substantial errors due to wire spatial resolution issues. These errors persist to wavelengths far beyond those which we might otherwise assume from simple isotropic assumptions of small-scale motions. The effects of hot-wire length-to-diameter ratio (l/d) are also briefly investigated. For the moderate wire Reynolds numbers investigated here, reducing l/d from 200 to 100 has a detrimental effect on measured turbulent fluctuations at a wide range of energetic scales, affecting both the broadband intensity and the energy spectra.

scholarly journals Rates of Dissipation of Turbulent Kinetic Energy in a High Reynolds Number Tidal Channel

2016 ◽  
Vol 33 (4) ◽  
pp. 817-837 ◽  
Author(s):  
Justine M. McMillan ◽  
Alex E. Hay ◽  
Rolf G. Lueck ◽  
Fabian Wolk
Keyword(s):  
Reynolds Number ◽  
Kinetic Energy ◽  
Turbulent Kinetic Energy ◽  
High Speed ◽  
Flow Speed ◽  
Small Scale ◽  
Tidal Channel ◽  
Dissipation Rates ◽  
Channel Separation ◽  
Shear Probe

AbstractThe ability to estimate the rate of dissipation (ε) of turbulent kinetic energy at middepth in a high-speed tidal channel using broadband acoustic Doppler current profilers (ADCPs) is assessed by making comparisons to direct measurements of ε obtained using shear probes mounted on a streamlined underwater buoy. The investigation was carried out in Grand Passage, Nova Scotia, Canada, where the depth-averaged flow speed reached 2 m s−1 and the Reynolds number was 8 × 107. The speed bin–averaged dissipation rates estimated from the ADCP data agree with the shear probe data to within a factor of 2. Both the ADCP and the shear probe measurements indicate a linear dependence of ε on the cube of the flow speed during flood and much lower dissipation rates during ebb. The ebb–flood asymmetry and the small-scale intermittency in ε are also apparent in the lognormal distributions of the shear probe data. Possible sources of bias and error in the ε estimates are investigated, and the most likely causes of the discrepancy between the ADCP and shear probe estimates are the cross-channel separation of the instruments and the high degree of spatial variability that exists in the channel.

scholarly journals Developments in turbulence research: a review based on the 1999 Programme of the Isaac Newton Institute, Cambridge

2001 ◽  
Vol 436 ◽  
pp. 353-391 ◽  
Author(s):  
J. C. R. HUNT ◽  
N. D. SANDHAM ◽  
J. C. VASSILICOS ◽  
B. E. LAUNDER ◽  
P. A. MONKEWITZ ◽  
...  
Keyword(s):  
Reynolds Number ◽  
Large Eddy Simulation ◽  
Numerical Simulations ◽  
Turbulent Flows ◽  
Small Scale ◽  
Isaac Newton ◽  
Eddy Simulation ◽  
Isaac Newton Institute ◽  
Large Eddy ◽  
Eddy Structures

Recent research is making progress in framing more precisely the basic dynamical and statistical questions about turbulence and in answering them. It is helping both to define the likely limits to current methods for modelling industrial and environmental turbulent flows, and to suggest new approaches to overcome these limitations. Our selective review is based on the themes and new results that emerged from more than 300 presentations during the Programme held in 1999 at the Isaac Newton Institute, Cambridge, UK, and on research reported elsewhere. A general conclusion is that, although turbulence is not a universal state of nature, there are certain statistical measures and kinematic features of the small-scale flow field that occur in most turbulent flows, while the large-scale eddy motions have qualitative similarities within particular types of turbulence defined by the mean flow, initial or boundary conditions, and in some cases, the range of Reynolds numbers involved. The forced transition to turbulence of laminar flows caused by strong external disturbances was shown to be highly dependent on their amplitude, location, and the type of flow. Global and elliptical instabilities explain much of the three-dimensional and sudden nature of the transition phenomena. A review of experimental results shows how the structure of turbulence, especially in shear flows, continues to change as the Reynolds number of the turbulence increases well above about 104 in ways that current numerical simulations cannot reproduce. Studies of the dynamics of small eddy structures and their mutual interactions indicate that there is a set of characteristic mechanisms in which vortices develop (vortex stretching, roll-up of instability sheets, formation of vortex tubes) and another set in which they break up (through instabilities and self- destructive interactions). Numerical simulations and theoretical arguments suggest that these often occur sequentially in randomly occurring cycles. The factors that determine the overall spectrum of turbulence were reviewed. For a narrow distribution of eddy scales, the form of the spectrum can be defined by characteristic forms of individual eddies. However, if the distribution covers a wide range of scales (as in elongated eddies in the ‘wall’ layer of turbulent boundary layers), they collectively determine the spectra (as assumed in classical theory). Mathematical analyses of the Navier–Stokes and Euler equations applied to eddy structures lead to certain limits being defined regarding the tendencies of the vorticity field to become infinitely large locally. Approximate solutions for eigen modes and Fourier components reveal striking features of the temporal, near-wall structure such as bursting, and of the very elongated, spatial spectra of sheared inhomogeneous turbulence; but other kinds of eddy concepts are needed in less structured parts of the turbulence. Renormalized perturbation methods can now calculate consistently, and in good agreement with experiment, the evolution of second- and third-order spectra of homogeneous and isotropic turbulence. The fact that these calculations do not explicitly include high-order moments and extreme events, suggests that they may play a minor role in the basic dynamics. New methods of approximate numerical simulations of the larger scales of turbulence or ‘very large eddy simulation’ (VLES) based on using statistical models for the smaller scales (as is common in meteorological modelling) enable some turbulent flows with a non-local and non-equilibrium structure, such as impinging or convective flows, to be calculated more efficiently than by using large eddy simulation (LES), and more accurately than by using ‘engineering’ models for statistics at a single point. Generally it is shown that where the turbulence in a fluid volume is changing rapidly and is very inhomogeneous there are flows where even the most complex ‘engineering’ Reynolds stress transport models are only satisfactory with some special adaptation; this may entail the use of transport equations for the third moments or non-universal modelling methods designed explicitly for particular types of flow. LES methods may also need flow-specific corrections for accurate modelling of different types of very high Reynolds number turbulent flow including those near rigid surfaces.This paper is dedicated to the memory of George Batchelor who was the inspiration of so much research in turbulence and who died on 30th March 2000. These results were presented at the last fluid mechanics seminar in DAMTP Cambridge that he attended in November 1999.

Scale-to-scale energy and enstrophy transport in two-dimensional Rayleigh–Taylor turbulence

2015 ◽  
Vol 786 ◽  
pp. 294-308 ◽  
Author(s):  
Quan Zhou ◽  
Yong-Xiang Huang ◽  
Zhi-Ming Lu ◽  
Yu-Lu Liu ◽  
Rui Ni
Keyword(s):  
Kinetic Energy ◽  
Thermal Energy ◽  
Scaling Laws ◽  
Isotropic Turbulence ◽  
Thermal Dissipation ◽  
Small Scale ◽  
Two Dimensional ◽  
Space Technique ◽  
Small Scales ◽  
The Mean

We apply a recently developed filtering approach, i.e. filter-space technique (FST), to study the scale-to-scale transport of kinetic energy, thermal energy, and enstrophy in two-dimensional (2D) Rayleigh–Taylor (RT) turbulence. Although the scaling laws of the energy cascades in 2D RT systems follow the Bolgiano–Obukhov (BO59) scenario due to buoyancy forces, the kinetic energy is still found to be, on average, dynamically transferred to large scales by an inverse cascade, while both the mean thermal energy and the mean enstrophy move towards small scales by forward cascades. In particular, there is a reasonably extended range over which the transfer rate of thermal energy is scale-independent and equals the corresponding thermal dissipation rate at different times. This range functions similarly to the inertial range for the kinetic energy in the homogeneous and isotropic turbulence. Our results further show that at small scales the fluctuations of the three instantaneous local fluxes are highly asymmetrically distributed and there is a strong correlation between any two fluxes. These small-scale features are signatures of the mixing and dissipation of fluids with steep temperature gradients at the fluid interfaces.

Sign in / Sign up

Export Citation Format

Share Document