5 results
The intermittency boundary in stratified plane Couette flow
- Enrico Deusebio, C. P. Caulfield, J. R. Taylor
-
- Journal:
- Journal of Fluid Mechanics / Volume 781 / 25 October 2015
- Published online by Cambridge University Press:
- 18 September 2015, pp. 298-329
-
- Article
- Export citation
-
We study stratified turbulence in plane Couette flow using direct numerical simulations. Two external dimensionless parameters control the dynamics, the Reynolds number $\mathit{Re}=Uh/{\it\nu}$ and the bulk Richardson number $\mathit{Ri}=g{\it\alpha}_{V}Th/U^{2}$, where $U$ and $T$ are half the velocity and temperature difference between the two walls respectively, $h$ is the half channel depth, ${\it\nu}$ is the kinematic viscosity and $g{\it\alpha}_{V}$ is the buoyancy parameter. We focus on spatio-temporal intermittency due to stratification and we explore the boundary between fully developed turbulence and intermittent flow in the $\mathit{Re}{-}\mathit{Ri}$ plane. The structures populating the intermittent flow regime show coexistence between laminar and turbulent patches, and we demonstrate that there are qualitative differences between the previously studied low-$\mathit{Re}$ low-$\mathit{Ri}$ intermittent regime and the high-$\mathit{Re}$ high-$\mathit{Ri}$ intermittent regime. At low-$\mathit{Re}$ low-$\mathit{Ri}$, turbulent regions span the entire gap, whereas at high-$\mathit{Re}$ high-$\mathit{Ri}$, turbulence is confined vertically with complex dynamics arising from interacting turbulent layers. Consistent with a previous investigation of Flores & Riley (Boundary-Layer Meteorol., vol. 129 (2), 2010, pp. 241–259), we present evidence suggesting that intermittency in the asymptotic regime of high-$\mathit{Re}$ Couette flows appears for $L^{+}<200$, where $L^{+}=Lu_{{\it\tau}}/{\it\nu}$, with $L$ being the Monin–Obukhov length scale, $L=u_{{\it\tau}}^{3}/C_{{\it\kappa}}q_{w}$, $q_{w}$ the wall heat flux, $C_{{\it\kappa}}$ the von Kármán constant and $u_{{\it\tau}}=\sqrt{{\it\tau}_{w}/{\it\rho}_{0}}$ the friction velocity determined from the wall shear stress ${\it\tau}_{w}$, where ${\it\rho}_{0}$ is the constant background density. We also consider the mixing as quantified by various versions of the flux Richardson number $\mathit{Ri}_{f}$, defined as the ratio of the conversion rate from kinetic to potential energy to the turbulent kinetic energy injection rate due to shear. We investigate how laminar and turbulent regions separately contribute to the overall mixing. Remarkably, we find that although fluctuations are greatly suppressed in the laminar regions, $\mathit{Ri}_{f}$ does not change significantly compared with its value in turbulent regions. As we observe a tight coupling between the mean temperature and velocity fields, we demonstrate that both Monin–Obukhov self-similarity theory (Monin & Obukhov, Contrib. Geophys. Inst. Acad. Sci. USSR, vol. 151, 1954, pp. 163–187) and the explicit algebraic model of Lazeroms et al. (J. Fluid Mech., vol. 723, 2013, pp. 91–125) predict the mean profiles well. We thus use these models to trace out the boundary between fully developed turbulence and intermittency in the $\mathit{Re}{-}\mathit{Ri}$ plane.
A numerical study of the unstratified and stratified Ekman layer
- Enrico Deusebio, G. Brethouwer, P. Schlatter, E. Lindborg
-
- Journal:
- Journal of Fluid Mechanics / Volume 755 / 25 September 2014
- Published online by Cambridge University Press:
- 26 August 2014, pp. 672-704
-
- Article
- Export citation
-
We study the turbulent Ekman layer at moderately high Reynolds number, $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}1600 < \mathit{Re} = \delta _{E}G/\nu < 3000$, using direct numerical simulations (DNS). Here, $\delta _{E} = \sqrt{2\nu /f}$ is the laminar Ekman layer thickness, $G$ the geostrophic wind, $\nu $ the kinematic viscosity and $f$ is the Coriolis parameter. We present results for both neutrally, moderately and strongly stably stratified conditions. For unstratified cases, large-scale roll-like structures extending from the outer region down to the wall are observed. These structures have a clear dominant frequency and could be related to periodic oscillations or instabilities developing near the low-level jet. We discuss the effect of stratification and $\mathit{Re}$ on one-point and two-point statistics. In the strongly stratified Ekman layer we observe stable co-existing large-scale laminar and turbulent patches appearing in the form of inclined bands, similar to other wall-bounded flows. For weaker stratification, continuously sustained turbulence strongly affected by buoyancy is produced. We discuss the scaling of turbulent length scales, height of the Ekman layer, friction velocity, veering angle at the wall and heat flux. The boundary-layer thickness, the friction velocity and the veering angle depend on $Lf/u_\tau $, where $u_\tau $ is the friction velocity and $L$ the Obukhov length scale, whereas the heat fluxes appear to scale with $L^+=L u_\tau /\nu $.
Helicity in the Ekman boundary layer
- Enrico Deusebio, Erik Lindborg
-
- Journal:
- Journal of Fluid Mechanics / Volume 755 / 25 September 2014
- Published online by Cambridge University Press:
- 26 August 2014, pp. 654-671
-
- Article
- Export citation
-
Helicity, which is defined as the scalar product of velocity and vorticity, $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}{\mathcal{H}} = {\boldsymbol {u}} \boldsymbol {\cdot }\boldsymbol{\omega}$, is an inviscidly conserved quantity in a barotropic fluid. Mean helicity is zero in flows that are parity invariant. System rotation breaks parity invariance and has therefore the potential of giving rise to non-zero mean helicity. In this paper we study the helicity dynamics in the incompressible Ekman boundary layer. Evolution equations for the mean field helicity and the mean turbulent helicity are derived and it is shown that pressure flux injects helicity at a rate $ 2 \varOmega G^2 $ over the total depth of the Ekman layer, where $ G $ is the geostrophic wind far from the wall and $ {\boldsymbol{\Omega}} = \varOmega {\boldsymbol {e}}_y $ is the rotation vector and $ {\boldsymbol {e}}_y $ is the wall-normal unit vector. Thus right-handed/left-handed helicity will be injected if $ \varOmega $ is positive/negative. We also show that in the uppermost part of the boundary layer there is a net helicity injection with opposite sign as compared with the totally integrated injection. Isotropic relations for the helicity dissipation and the helicity spectrum are derived and it is shown that it is sufficient to measure two transverse velocity components and use Taylor’s hypothesis in the mean flow direction in order to measure the isotropic helicity spectrum. We compare the theoretical predictions with a direct numerical simulation of an Ekman boundary layer and confirm that there is a preference for right-handed helicity in the lower part of the Ekman layer and left-handed helicity in the uppermost part when $ \varOmega > 0 $. In the logarithmic range, the helicity dissipation conforms to isotropic relations. On the other hand, spectra show significant departures from isotropic conditions, suggesting that the Reynolds number considered in the study is not sufficiently large for isotropy to be valid in a wide range of scales. Our analytical and numerical results strongly suggest that there is a turbulent helicity cascade of right-handed helicity in the logarithmic range of the atmospheric boundary layer when $\varOmega >0$, consistent with recent measurements by Koprov, Koprov, Ponomarev & Chkhetiani (Dokl. Phys., vol. 50, 2005, pp. 419–422). The isotropic relations which are derived may facilitate future measurements of the helicity spectrum in the atmospheric boundary layer as well as in controlled wind tunnel experiments.
Third-order structure functions in rotating and stratified turbulence: a comparison between numerical, analytical and observational results
- Enrico Deusebio, P. Augier, E. Lindborg
-
- Journal:
- Journal of Fluid Mechanics / Volume 755 / 25 September 2014
- Published online by Cambridge University Press:
- 19 August 2014, pp. 294-313
-
- Article
- Export citation
-
First, we review analytical and observational studies on third-order structure functions including velocity and buoyancy increments in rotating and stratified turbulence and discuss how these functions can be used in order to estimate the flux of energy through different scales in a turbulent cascade. In particular, we suggest that the negative third-order velocity–temperature–temperature structure function that was measured by Lindborg & Cho (Phys. Rev. Lett., vol. 85, 2000, p. 5663) using stratospheric aircraft data may be used in order to estimate the downscale flux of available potential energy (APE) through the mesoscales. Then, we calculate third-order structure functions from idealized simulations of forced stratified and rotating turbulence and compare with mesoscale results from the lower stratosphere. In the range of scales with a downscale energy cascade of kinetic energy (KE) and APE we find that the third-order structure functions display a negative linear dependence on separation distance $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}} r $, in agreement with observation and supporting the interpretation of the stratospheric data as evidence of a downscale energy cascade. The spectral flux of APE can be estimated from the relevant third-order structure function. However, while the sign of the spectral flux of KE is correctly predicted by using the longitudinal third-order structure functions, its magnitude is overestimated by a factor of two. We also evaluate the third-order velocity structure functions that are not parity invariant and therefore display a cyclonic–anticyclonic asymmetry. In agreement with the results from the stratosphere, we find that these functions have an approximate $ r^{2} $-dependence, with strong dominance of cyclonic motions.
The route to dissipation in strongly stratified and rotating flows
- Enrico Deusebio, A. Vallgren, E. Lindborg
-
- Journal:
- Journal of Fluid Mechanics / Volume 720 / 10 April 2013
- Published online by Cambridge University Press:
- 27 February 2013, pp. 66-103
-
- Article
-
- You have access Access
- Open access
- Export citation
-
We investigate the route to dissipation in strongly stratified and rotating systems through high-resolution numerical simulations of the Boussinesq equations (BQs) and the primitive equations (PEs) in a triply periodic domain forced at large scales. By applying geostrophic scaling to the BQs and using the same horizontal length scale in defining the Rossby and the Froude numbers, $\mathit{Ro}$ and $\mathit{Fr}$, we show that the PEs can be obtained from the BQs by taking the limit ${\mathit{Fr}}^{2} / {\mathit{Ro}}^{2} \rightarrow 0$. When ${\mathit{Fr}}^{2} / {\mathit{Ro}}^{2} $ is small the difference between the results from the BQ and the PE simulations is shown to be small. For large rotation rates, quasi-geostrophic dynamics are recovered with a forward enstrophy cascade and an inverse energy cascade. As the rotation rate is reduced, a fraction of the energy starts to cascade towards smaller scales, leading to a shallowing of the horizontal spectra from ${ k}_{h}^{- 3} $ to ${ k}_{h}^{- 5/ 3} $ at the small-scale end. The vertical spectra show a similar transition as the horizontal spectra and we find that Charney isotropy is approximately valid also at larger wavenumbers than the transition wavenumber. The high resolutions employed allow us to capture both ranges within the same simulation. At the transition scale, kinetic energy in the rotational and in the horizontally divergent modes attain comparable values. The divergent energy is several orders of magnitude larger than the quasi-geostrophic divergent energy given by the $\Omega $-equation. The amount of energy cascading downscale is mainly controlled by the rotation rate, with a weaker dependence on the stratification. A larger degree of stratification favours a downscale energy cascade. For intermediate degrees of rotation and stratification, a constant energy flux and a constant enstrophy flux coexist within the same range of scales. In this range, the enstrophy flux is a result of triad interactions involving three geostrophic modes, while the energy flux is a result of triad interactions involving at least one ageostrophic mode, with a dominant contribution from interactions involving two ageostrophic and one geostrophic mode. Dividing the ageostrophic motions into two classes depending on the sign of the linear wave frequency, we show that the energy transfer is for the largest part supported by interactions within the same class, ruling out the wave–wave–vortex resonant triad interaction as a mean of the downscale energy transfer. The role of inertia-gravity waves is studied through analyses of time-frequency spectra of single Fourier modes. At large scales, distinct peaks at frequencies predicted for linear waves are observed, whereas at small scales no clear wave activity is observed. Triad interactions show a behaviour which is consistent with turbulent dynamics, with a large exchange of energy in triads with one small and two large comparable wavenumbers. The exchange of energy is mainly between the modes with two comparable wavenumbers.