Mixing efficiency in large-eddy simulations of stratified turbulence

2018 ◽  
Vol 849 ◽  
pp. 373-394 ◽  
Author(s):  
Sina Khani

The irreversible mixing efficiency is studied using large-eddy simulations (LES) of stratified turbulence, where three different subgrid-scale (SGS) parameterizations are employed. For comparison, direct numerical simulations (DNS) and hyperviscosity simulations are also performed. In the regime of stratified turbulence where $Fr_{v}\sim 1$, the irreversible mixing efficiency $\unicode[STIX]{x1D6FE}_{i}$ in LES scales like $1/(1+2Pr_{t})$, where $Fr_{v}$ and $Pr_{t}$ are the vertical Froude number and turbulent Prandtl number, respectively. Assuming a unit scaling coefficient and $Pr_{t}=1$, $\unicode[STIX]{x1D6FE}_{i}$ goes to a constant value $1/3$, in agreement with DNS results. In addition, our results show that the irreversible mixing efficiency in LES, while consistent with this prediction, depends on SGS parameterizations and the grid spacing $\unicode[STIX]{x1D6E5}$. Overall, the LES approach can reproduce mixing efficiency results similar to those from the DNS approach if $\unicode[STIX]{x1D6E5}\lesssim L_{o}$, where $L_{o}$ is the Ozmidov scale. In this situation, the computational costs of numerical simulations are significantly reduced because LES runs require much smaller computational resources in comparison with expensive DNS runs.

2015 ◽  
Vol 773 ◽  
pp. 327-344 ◽  
Author(s):  
Sina Khani ◽  
Michael L. Waite

The dynamic Smagorinsky model for large eddy simulation (LES) of stratified turbulence is studied in this paper. A maximum grid spacing criterion of ${\it\Delta}/L_{b}<0.24$ is found in order to capture several of the key characteristics of stratified turbulence, where ${\it\Delta}$ is the filter scale and $L_{b}$ is the buoyancy scale. These results show that the dynamic Smagorinsky model needs a grid spacing approximately twice as large as the regular Smagorinsky model to reproduce similar results. This improvement on the regular Smagorinsky eddy viscosity approach increases the accuracy of results at small resolved scales while decreasing the computational costs because it allows larger ${\it\Delta}$. In addition, the eddy dissipation spectra in LES of stratified turbulence present anisotropic features, taking energy out of large horizontal but small vertical scales. This trend is not seen in the non-stratified cases, where the subgrid-scale energy transfer is isotropic. Statistics of the dynamic Smagorinsky coefficient $c_{s}$ are investigated; its distribution is peaked around zero, and its standard deviations decrease slightly with increasing stratification. In line with previous findings for unstratified turbulence, regions of increased shear favour smaller $c_{s}$ values; in stratified turbulence, the spatial distribution of the shear, and hence $c_{s}$, is dominated by a layerwise pancake structure. These results show that the dynamic Smagorinsky model presents a promising approach for LES when isotropic buoyancy-scale resolving grids are employed.


2020 ◽  
Vol 148 (10) ◽  
pp. 4299-4311
Author(s):  
Sina Khani ◽  
Michael L. Waite

AbstractSubgrid-scale (SGS) parameterizations in atmosphere and ocean models are often defined independently in the horizontal and vertical directions because the grid spacing is not the same in these directions (anisotropic grids). In this paper, we introduce a new anisotropic SGS model in large-eddy simulations (LES) of stratified turbulence based on horizontal filtering of the equations of motion. Unlike the common horizontal SGS parameterizations in atmosphere and ocean models, the vertical derivatives of the horizontal SGS fluxes are included in our anisotropic SGS scheme, and therefore the horizontal and vertical SGS dissipation mechanisms are not disconnected in the newly developed model. Our model is tested with two vertical grid spacings and various horizontal resolutions, where the horizontal grid spacing is comparatively larger than that in the vertical. Our anisotropic LES model can successfully reproduce the results of direct numerical simulations, while the computational cost is significantly reduced in the LES. We suggest the new anisotropic SGS model as an alternative to current SGS parameterizations in atmosphere and ocean models, in which the schemes for horizontal and vertical scales are often decoupled. The new SGS scheme may improve the dissipative performance of atmosphere and ocean models without adding any backscatter or other energizing terms at small horizontal scales.


2014 ◽  
Vol 754 ◽  
pp. 75-97 ◽  
Author(s):  
Sina Khani ◽  
Michael L. Waite

AbstractIn this paper large-eddy simulations (LES) of forced stratified turbulence using two common subgrid scale (SGS) models, the Kraichnan and Smagorinsky models, are studied. As found in previous studies using regular and hyper-viscosity, vorticity contours show elongated horizontal motions, which are layered in the vertical direction, along with intermittent Kelvin–Helmholtz (KH) instabilities. Increased stratification causes the layer thickness to collapse towards the dissipation scale, ultimately suppressing these instabilities. The vertical energy spectra are relatively flat out to a local maximum, which varies with the buoyancy frequency $\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}}N$. The horizontal energy spectra depend on the grid spacing $\varDelta $; if the resolution is fine enough, the horizontal spectrum shows an approximately $-5/3$ slope along with a bump at the buoyancy wavenumber $k_b = N/u_{rms}$, where $u_{rms}$ is the root-mean-square (r.m.s.) velocity. Our results show that there is a critical value of the grid spacing $\varDelta $, below which dynamics of stratified turbulence are well-captured in LES. This critical $\varDelta $ depends on the buoyancy scale $L_b$ and varies with different SGS models: the Kraichnan model requires $\varDelta < 0.47 L_b$, while the Smagorinsky model requires $\varDelta < 0.17 L_b$. In other words, the Smagorinsky model is significantly more costly than the Kraichnan approach, as it requires three times the resolution to adequately capture stratified turbulence.


Author(s):  
B. Franzelli ◽  
E. Riber ◽  
B. Cuenot ◽  
M. Ihme

Numerical simulations are regarded as an essential tool for improving the design of combustion systems since they can provide information that is complementary to experiments. However, although numerical simulations have already been successfully applied to the prediction of temperature and species concentration in turbulent flames, the production of soot is far from being conclusive due to the complexity of the processes involved in soot production. In this context, first Large Eddy Simulations (LES) of soot production in turbulent flames are reported in the literature in laboratory-scale configurations, thereby confirming the feasibility of the approach. However numerous modeling and numerical issues have not been completely solved. Moreover, validation of the models through comparisons with measurements in realistic complex flows typical of aero-engines is still rare. This work therefore proposes to evaluate the LES approach for the prediction of soot production in an experimental swirl-stabilized non-premixed ethylene/air aero-engine combustor, for which soot and flame data are available. Two simulations are carried out using a two-equation soot model to compare the performance of a hybrid chemical description (reduced chemistry for the flame structure/tabulated chemistry for soot precursor chemistry) to a classical full tabulation method. Discrepancies of soot concentration between the two LES calculations will be analyzed and the sensitivity to the chemical models will be investigated.


1994 ◽  
Vol 47 (6S) ◽  
pp. S163-S165
Author(s):  
Douglas G. Dommermuth ◽  
Rebecca C. Y. Mui

Direct numerical simulations and large-eddy simulations of turbulent free-surface flows are currently being performed to investigate the roughening of the surface, and the scattering, radiation, and dissipation of waves by turbulence. The numerical simulation of turbulent free-surface flows is briefly reviewed. The numerical, modeling, and hardware issues are discussed.


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.


Sign in / Sign up

Export Citation Format

Share Document