Turbulent/non-turbulent interfaces in wakes in stably stratified fluids

2016 ◽  
Vol 797 ◽  
Author(s):  
Tomoaki Watanabe ◽  
James J. Riley ◽  
Stephen M. de Bruyn Kops ◽  
Peter J. Diamessis ◽  
Qi Zhou

We report on a study, employing direct numerical simulations, of the turbulent/non-turbulent interface of a wake in a stably stratified fluid. It is found that thresholds for both enstrophy and potential enstrophy are needed to identify the interface. Using conditional averaging relative to the location of the interface, various quantities of interest are examined. The thickness of the interface is found to scale with the Kolmogorov scale. From an examination of the Ozmidov and Kolmogorov length scales as well as the buoyancy Reynolds number, it is found that the buoyancy Reynolds number decreases and becomes of order 1 near the interface, indicating the suppression of the turbulence there by the stable stratification. Finally the overall rate of loss of energy due to internal wave radiation is found to be comparable to the overall rate of loss due to turbulent kinetic energy dissipation.

2015 ◽  
Vol 767 ◽  
pp. 342-363 ◽  
Author(s):  
J. Bertram

AbstractFollowing Malkus’s (J. Fluid Mech., vol. 1, 1956, pp. 521–539) proposal that turbulent Poiseuille channel flow maximises total viscous dissipation $D$, a variety of variational procedures have been explored involving the maximisation of different flow quantities under different constraints. However, the physical justification for these variational procedures has remained unclear. Here we address more recent claims that mean flow viscous dissipation $D_{m}$ should be maximised on the basis of a statistical stability argument, and that maximising $D_{m}$ yields realistic mean velocity profiles (Malkus, J. Fluid Mech., vol. 489, 2003, pp. 185–198). We clarify the connection between maximising $D_{m}$ and other flow quantities, verify Malkus & Smith’s, (J. Fluid Mech., vol. 208, 1989, pp. 479–507) claim that maximising the ‘efficiency’ yields realistic profiles and show that, in contrast, maximising $D_{m}$ does not yield realistic mean velocity profiles as recently claimed. This leads us to revisit Malkus’s statistical stability argument for maximising $D_{m}$ and to address some of its limitations. We propose an alternative statistical stability argument leading to a principle of minimum kinetic energy for fixed pressure gradient, which suggests a principle of maximum $D$ for fixed Reynolds number under certain conditions. We discuss possible ways to reconcile these conflicting results, focusing on the choice of constraints.


2019 ◽  
Vol 49 (7) ◽  
pp. 1873-1887 ◽  
Author(s):  
Kaiwen Zheng ◽  
Maxim Nikurashin

AbstractRecent microstructure observations in the Southern Ocean report enhanced internal gravity waves and turbulence in the frontal regions of the Antarctic Circumpolar Current extending a kilometer above rough bottom topography. Idealized numerical simulations and linear theory show that geostrophic flows impinging on rough small-scale topography are very effective generators of internal waves and estimate vigorous wave radiation, breaking, and turbulence within a kilometer above bottom. However, both idealized simulations and linear theory assume periodic and spatially uniform topography and tend to overestimate the observed levels of turbulent energy dissipation locally at the generation sites. In this study, we explore the downstream evolution and remote dissipation of internal waves generated by geostrophic flows using a series of numerical, realistic topography simulations and parameters typical of Drake Passage. The results show that significant levels of internal wave kinetic energy and energy dissipation are present downstream of the rough topography, internal wave generation site. About 30%–40% of the energy dissipation occurs locally over the rough topography region, where internal waves are generated. The rest of the energy dissipation takes place remotely and decays downstream of the generation site with an e-folding length scale of up to 20–30 km. The model we use is two-dimensional with enhanced viscosity coefficients, and hence it can result in the underestimation of the remote wave dissipation and its decay length scale. The implications of our results for turbulent energy dissipation observations and mixing parameterizations are discussed.


1997 ◽  
Vol 342 ◽  
pp. 231-261 ◽  
Author(s):  
FRANK G. JACOBITZ ◽  
SUTANU SARKAR ◽  
CHARLES W. VAN ATTA

Direct numerical simulations (DNS) are performed to investigate the evolution of turbulence in a uniformly sheared and stably stratified flow. The spatial discretization is accomplished by a spectral collocation method, and the solution is advanced in time with a third-order Runge–Kutta scheme. The turbulence evolution is found to depend strongly on at least three parameters: the gradient Richardson number Ri, the initial value of the Taylor microscale Reynolds number Reλ, and the initial value of the shear number SK/<ε. The effect of each parameter is individually studied while the remaining parameters are kept constant. The evolution of the turbulent kinetic energy K is found to follow approximately an exponential law. The shear number SK/<ε, whose effect has not been investigated in previous studies, was found to have a strong non-monotone influence on the turbulence evolution. Larger values of the shear number do not necessarily lead to a larger value of the eventual growth rate of the turbulent kinetic energy. Variation of the Reynolds number Reλ indicated that the turbulence growth rate tends to become insensitive to Reλ at the higher end of the Reλ range studied here. The dependence of the critical Richardson number Ricr, which separates asymptotic growth of the turbulent kinetic energy K from asymptotic decay, on the initial values of the Reynolds number Reλ and the shear number SK/<ε was also obtained. It was found that the critical Richardson number varied over the range 0.04<Ricr<0.17 in our DNS due to its strong dependence on Reynolds and shear numbers.


2020 ◽  
Vol 77 (9) ◽  
pp. 3193-3210
Author(s):  
Tiago Pestana ◽  
Matthias Thalhammer ◽  
Stefan Hickel

Abstract We present direct numerical simulations of inertia–gravity waves breaking in the middle–upper mesosphere. We consider two different altitudes, which correspond to the Reynolds number of 28 647 and 114 591 based on wavelength and buoyancy period. While the former was studied by Remmler et al., it is here repeated at a higher resolution and serves as a baseline for comparison with the high-Reynolds-number case. The simulations are designed based on the study of Fruman et al., and are initialized by superimposing primary and secondary perturbations to the convectively unstable base wave. Transient growth leads to an almost instantaneous wave breaking and secondary bursts of turbulence. We show that this process is characterized by the formation of fine flow structures that are predominantly located in the vicinity of the wave’s least stable point. During the wave breakdown, the energy dissipation rate tends to be an isotropic tensor, whereas it is strongly anisotropic in between the breaking events. We find that the vertical kinetic energy spectra exhibit a clear 5/3 scaling law at instants of intense energy dissipation rate and a cubic power law at calmer periods. The term-by-term energy budget reveals that the pressure term is the most important contributor to the global energy budget, as it couples the vertical and the horizontal kinetic energy. During the breaking events, the local energy transfer is predominantly from the mean to the fluctuating field and the kinetic energy production is in balance with the pseudo kinetic energy dissipation rate.


2017 ◽  
Vol 817 ◽  
pp. 61-79 ◽  
Author(s):  
L. Djenidi ◽  
N. Lefeuvre ◽  
M. Kamruzzaman ◽  
R. A. Antonia

The Reynolds number dependence of the non-dimensional mean turbulent kinetic energy dissipation rate$C_{\unicode[STIX]{x1D716}}=\overline{\unicode[STIX]{x1D716}}L/u^{\prime 3}$(where$\unicode[STIX]{x1D716}$is the mean turbulent kinetic energy dissipation rate,$L$is an integral length scale and$u^{\prime }$is the velocity root-mean-square) is investigated in decaying turbulence. Expressions for$C_{\unicode[STIX]{x1D716}}$in homogeneous isotropic turbulent (HIT), as approximated by grid turbulence, and in local HIT, as on the axis of the far field of a turbulent round jet, are developed from the Navier–Stokes equations within the framework of a scale-by-scale energy budget. The analysis shows that when turbulence decays/evolves in compliance with self-preservation (SP),$C_{\unicode[STIX]{x1D716}}$remains constant for a given flow condition, e.g. a given initial Reynolds number. Measurements in grid turbulence, which does not satisfy SP, and on the axis in the far field of a round jet, which does comply with SP, show that$C_{\unicode[STIX]{x1D716}}$decreases in the former case and remains constant in the latter, thus supporting the theoretical results. Further, while$C_{\unicode[STIX]{x1D716}}$can remain constant during the decay for a given initial Reynolds number, both the theory and measurements show that it decreases towards a constant,$C_{\unicode[STIX]{x1D716},\infty }$, as$Re_{\unicode[STIX]{x1D706}}$increases. This trend, in agreement with existing data, is not inconsistent with the possibility that$C_{\unicode[STIX]{x1D716}}$tends to a universal constant.


1997 ◽  
Vol 343 ◽  
pp. 111-130 ◽  
Author(s):  
DAREK BOGUCKI ◽  
J. ANDRZEJ DOMARADZKI ◽  
P. K. YEUNG

Direct numerical simulations of passive scalars, with Prandtl numbers Pr=3, 5, and 7, advected by turbulence at three low Reynolds numbers were performed. The energy spectra are self-similar under the Kolmogorov scaling and exhibit behaviour consistent with many other investigations: a short inertial range for the highest Reynolds number and the universal exponential form of the spectrum for all Reynolds numbers in the dissipation range. In all cases the passive scalar spectra collapse to a single self-similar curve under the Batchelor scaling and exhibit the k−1 range followed by an exponential fall-off. We attribute the applicability of the Batchelor scaling to our low-Reynolds-number flows to the universality of the energy dissipation spectra. The Batchelor range is observed for wavenumbers in general agreement with experimental observations but smaller than predicted by the classical estimates. The discrepancy is caused by the fact that the velocity scales responsible for the generation of the Batchelor range are in the vicinity of the wavenumber of the maximum energy dissipation, which is one order of magnitude less than the Kolmogorov wavenumber used in the classical theory. Two different functional forms of passive scalar spectra proposed by Batchelor and Kraichnan were fitted to the simulation results and it was found that the Kraichnan model agrees very well with the data while the Batchelor formula displays systematic deviations from the data. Implications of these differences for the experimental procedures to measure the energy and passive scalar dissipation rates in oceanographic flows are discussed.


Atmosphere ◽  
2021 ◽  
Vol 12 (2) ◽  
pp. 157
Author(s):  
Duane Rosenberg ◽  
Annick Pouquet ◽  
Raffaele Marino

We study in this paper the correlation between the buoyancy flux, the efficiency of energy dissipation and the linear and nonlinear components of potential vorticity, PV, a point-wise invariant of the Boussinesq equations, contrasting the three identified regimes of rotating stratified turbulence, namely wave-dominated, wave–eddy interactions and eddy-dominated. After recalling some of the main novel features of these flows compared to homogeneous isotropic turbulence, we specifically analyze three direct numerical simulations in the absence of forcing and performed on grids of 10243 points, one in each of these physical regimes. We focus in particular on the link between the point-wise buoyancy flux and the amount of kinetic energy dissipation and of linear and nonlinear PV. For flows dominated by waves, we find that the highest joint probability is for minimal kinetic energy dissipation (compared to the buoyancy flux), low dissipation efficiency and low nonlinear PV, whereas for flows dominated by nonlinear eddies, the highest correlation between dissipation and buoyancy flux occurs for weak flux and high localized nonlinear PV. We also show that the nonlinear potential vorticity is strongly correlated with high dissipation efficiency in the turbulent regime, corresponding to intermittent events, as observed in the atmosphere and oceans.


Sign in / Sign up

Export Citation Format

Share Document