scholarly journals Investigations of non-hydrostatic, stably stratified and rapidly rotating flows

2016 ◽  
Vol 801 ◽  
pp. 430-458 ◽  
Author(s):  
David Nieves ◽  
Ian Grooms ◽  
Keith Julien ◽  
Jeffrey B. Weiss

We present an investigation of rapidly rotating (small Rossby number $Ro\ll 1$) stratified turbulence where the stratification strength is varied from weak (large Froude number $Fr\gg 1$) to strong ($Fr\ll 1$). The investigation is set in the context of a reduced model derived from the Boussinesq equations that retains anisotropic inertia-gravity waves with order-one frequencies and highlights a regime of wave–eddy interactions. Numerical simulations of the reduced model are performed where energy is injected by a stochastic forcing of vertical velocity, which forces wave modes only. The simulations reveal two regimes: characterized by the presence of well-formed, persistent and thin turbulent layers of locally weakened stratification at small Froude numbers, and by the absence of layers at large Froude numbers. Both regimes are characterized by a large-scale barotropic dipole enclosed by small-scale turbulence. When the Reynolds number is not too large, a direct cascade of barotropic kinetic energy is observed, leading to total energy equilibration. We examine net energy exchanges that occur through vortex stretching and vertical buoyancy flux and diagnose the horizontal scales active in these exchanges. We find that the baroclinic motions inject energy directly to the largest scales of the barotropic mode, implying that the large-scale barotropic dipole is not the end result of an inverse cascade within the barotropic mode.

2001 ◽  
Vol 427 ◽  
pp. 205-239 ◽  
Author(s):  
GEORGE F. CARNEVALE ◽  
M. BRISCOLINI ◽  
P. ORLANDI

The buoyancy range, which represents a transition from large-scale wave-dominated motions to small-scale turbulence in the oceans and the atmosphere, is investigated through large-eddy simulations. The model presented here uses a continual forcing based on large-scale standing internal waves and has a spectral truncation in the isotropic inertial range. Evidence is presented for a break in the energy spectra from the anisotropic k−3 buoyancy range to the small-scale k−5/3 isotropic inertial range. Density structures that form during wave breaking and periods of high strain rate are analysed. Elongated vertical structures produced during periods of strong straining motion are found to collapse in the subsequent vertically compressional phase of the strain resulting in a zone or patch of mixed fluid.


2021 ◽  
Vol 2099 (1) ◽  
pp. 012033
Author(s):  
G V Zasko ◽  
P A Perezhogin ◽  
A V Glazunov ◽  
E V Mortikov ◽  
Y M Nechepurenko

Abstract Large-scale inclined organized structures in stably stratified turbulent shear flows were revealed in the numerical simulation and indirectly confirmed by the field measurements in the stable atmospheric boundary layer. Spatial scales and forms of these structures coincide with those of the optimal disturbances of a simplified linear model. In this paper, we clarify the relation between the organized structures and the optimal disturbances, analyzing a time series of turbulent fields obtained by the RANS model with eddy viscosity/diffusivity and stochastic forcing generating the small-scale turbulence.


2021 ◽  
Vol 933 ◽  
Author(s):  
Gregory P. Chini ◽  
Guillaume Michel ◽  
Keith Julien ◽  
Cesar B. Rocha ◽  
Colm-cille P. Caulfield

A multiscale reduced description of turbulent free shear flows in the presence of strong stabilizing density stratification is derived via asymptotic analysis of the Boussinesq equations in the simultaneous limits of small Froude and large Reynolds numbers. The analysis explicitly recognizes the occurrence of dynamics on disparate spatiotemporal scales, yielding simplified partial differential equations governing the coupled evolution of slow large-scale hydrostatic flows and fast small-scale isotropic instabilities and internal waves. The dynamics captured by the coupled reduced equations is illustrated in the context of two-dimensional strongly stratified Kolmogorov flow. A noteworthy feature of the reduced model is that the fluctuations are constrained to satisfy quasilinear (QL) dynamics about the comparably slowly varying large-scale fields. Crucially, this QL reduction is not invoked as an ad hoc closure approximation, but rather is derived in a physically relevant and mathematically consistent distinguished limit. Further analysis of the resulting slow–fast QL system shows how the amplitude of the fast stratified-shear instabilities is slaved to the slowly evolving mean fields to ensure the marginal stability of the latter. Physically, this marginal stability condition appears to be compatible with recent evidence of self-organized criticality in both observations and simulations of stratified turbulence. Algorithmically, the slaving of the fluctuation fields enables numerical simulations to be time-evolved strictly on the slow time scale of the hydrostatic flow. The reduced equations thus provide a solid mathematical foundation for future studies of three-dimensional strongly stratified turbulence in extreme parameter regimes of geophysical relevance and suggest avenues for new sub-grid-scale parametrizations.


2019 ◽  
Vol 4 (12) ◽  
Author(s):  
C. Marchioli ◽  
H. Bhatia ◽  
G. Sardina ◽  
L. Brandt ◽  
A. Soldati

Author(s):  
Scott Martin ◽  
Aleksandar Jemcov ◽  
Björn de Ruijter

Here the premixed Conditional Moment Closure (CMC) method is used to model the recent PIV and Raman turbulent, enclosed reacting methane jet data from DLR Stuttgart [1]. The experimental data has a rectangular test section at atmospheric pressure and temperature with a single inlet jet. A jet velocity of 90 m/s is used with an adiabatic flame temperature of 2,064 K. Contours of major species, temperature and velocities along with velocity rms values are provided. The conditional moment closure model has been shown to provide the capability to model turbulent, premixed methane flames with detailed chemistry and reasonable runtimes [2]. The simplified CMC model used here falls into the class of table lookup turbulent combustion models where the chemical kinetics are solved offline over a range of conditions and stored in a table that is accessed by the CFD code. Most table lookup models are based on the laminar 1-D flamelet equations, which assume the small scale turbulence does not affect the reaction rates, only the large scale turbulence has an effect on the reaction rates. The CMC model is derived from first principles to account for the effects of small scale turbulence on the reaction rates, as well as the effects of the large scale mixing, making it more versatile than other models. This is accomplished by conditioning the scalars with the reaction progress variable. By conditioning the scalars and accounting for the small scale mixing, the effects of turbulent fluctuations of the temperature on the reaction rates are more accurately modeled. The scalar dissipation is used to account for the effects of the small scale mixing on the reaction rates. The original premixed CMC model used a constant value of scalar dissipation, here the scalar dissipation is conditioned by the reaction progress variable. The steady RANS 3-D version of the open source CFD code OpenFOAM is used. Velocity, temperature and species are compared to the experimental data. Once validated, this CFD turbulent combustion model will have great utility for designing lean premixed gas turbine combustors.


2014 ◽  
Vol 955-959 ◽  
pp. 2425-2429 ◽  
Author(s):  
Yun Fei Li ◽  
Jian Guo Yang ◽  
Yan Yan Wang ◽  
Xiao Guo Wang

The purpose of this study is to construct a turbulent aggregation device which has specific performance for fine particle aggregation in flue gas. The device consists of two cylindrical pipes and an array of vanes. The pipes extending fully and normal to the gas stream induce large scale turbulence in the form of vortices, while the vanes downstream a certain distance from the pipes induce small one. The process of turbulent aggregation was numerically simulated by coupling the Eulerian multiphase model and population balance model together with a proposed aggregation kernel function taking the size and inertia of particles into account, and based on data of particles’ size distribution measured from the flue of one power plant. The results show that the large scale turbulence generated by pipes favours the aggregation of smaller particles (smaller than 1μm) notably, while the small scale turbulence benefits the aggregation of bigger particles (larger than 1μm) notably and enhances the uniformity of particle size distribution among different particle groups.


2002 ◽  
Vol 14 (7) ◽  
pp. 2475 ◽  
Author(s):  
L. Danaila ◽  
F. Anselmet ◽  
R. A. Antonia

2012 ◽  
Vol 8 (S294) ◽  
pp. 361-363
Author(s):  
A. V. Getling ◽  
O. S. Mazhorova ◽  
O. V. Shcheritsa

AbstractConvection is simulated numerically based on two-dimensional Boussinesq equations for a fluid layer with a specially chosen stratification such that the convective instability is much stronger in a thin subsurface sublayer than in the remaining part of the layer. The developing convective flow has a small-scale component superposed onto a basic large-scale roll flow.


1990 ◽  
Vol 142 ◽  
pp. 60-61
Author(s):  
Sydney D'Silva ◽  
Arnab Rai Choudhuri

Working under the hypothesis that magnetic flux in the sun is generated at the bottom of the convection zone, Choudhuri and Gilman (1987; Astrophys. J. 316, 788) found that a magnetic flux tube symmetric around the rotation axis, when released at the bottom of the convection zone, gets deflected by the Coriolis force and tends to move parallel to the rotation axis as it rises in the convection zone. As a result, all the flux emerges at rather high latitudes and the flux observed at the typical sunspot latitudes remains unexplained. Choudhuri(1989; Solar Physics, in press) finds that non-axisymmetric perturbations too cannot subdue the Coriolis force. In this paper, we no longer treat the convection zone to be passive as in the previous papers, but we consider the role of turbulence in the convection zone in inhibiting the Coriolis force. The interaction of the flux tubes with the turbulence is treated in a phenomenological way as follows: (1) Large scale turbulence on the scale of giant cells can physically drag the tubes outwards, thus pulling the flux towards lower latitudes by dominating over the Coriolis force. (2) Small scale turbulence of the size of the tubes can exchange angular momentum with the tube, thus suppressing the growth of the Coriolis force and making the tubes emerge at lower latitudes. Numerical simulations show that the giant cells can drag the tubes and make them emerge at lower latitudes only if the velocities within the giant cells are unrealistically large or if the radii of the flux tubes are as small as 10 km. However, small scale turbulence can successfully suppress the growth of the Coriolis force if the tubes have radii smaller than about 300 km which may not be unreasonable. Such flux tubes can then emerge at low latitudes where sunspots are seen.


2019 ◽  
Vol 867 ◽  
pp. 146-194 ◽  
Author(s):  
G. L. Richard ◽  
A. Duran ◽  
B. Fabrèges

We derive a two-dimensional depth-averaged model for coastal waves with both dispersive and dissipative effects. A tensor quantity called enstrophy models the subdepth large-scale turbulence, including its anisotropic character, and is a source of vorticity of the average flow. The small-scale turbulence is modelled through a turbulent-viscosity hypothesis. This fully nonlinear model has equivalent dispersive properties to the Green–Naghdi equations and is treated, both for the optimization of these properties and for the numerical resolution, with the same techniques which are used for the Green–Naghdi system. The model equations are solved with a discontinuous Galerkin discretization based on a decoupling between the hyperbolic and non-hydrostatic parts of the system. The predictions of the model are compared to experimental data in a wide range of physical conditions. Simulations were run in one-dimensional and two-dimensional cases, including run-up and run-down on beaches, non-trivial topographies, wave trains over a bar or propagation around an island or a reef. A very good agreement is reached in every cases, validating the predictive empirical laws for the parameters of the model. These comparisons confirm the efficiency of the present strategy, highlighting the enstrophy as a robust and reliable tool to describe wave breaking even in a two-dimensional context. Compared with existing depth-averaged models, this approach is numerically robust and adds more physical effects without significant increase in numerical complexity.


Sign in / Sign up

Export Citation Format

Share Document