scholarly journals A Mixed Scheme for Subgrid-Scale Fluxes in Cloud-Resolving Models

2010 ◽  
Vol 67 (11) ◽  
pp. 3692-3705 ◽  
Author(s):  
C.-H. Moeng ◽  
P. P. Sullivan ◽  
M. F. Khairoutdinov ◽  
D. A. Randall

Abstract A large-domain large-eddy simulation of a tropical deep convection system is used as a benchmark to derive and test a mixed subgrid-scale (SGS) scheme for scalar and momentum fluxes in cloud-resolving models (CRMs). The benchmark simulation resolves a broad range of scales ranging from mesoscale organizations, through gravity waves and individual clouds, down to energy-containing turbulent eddies. A spectral analysis shows that the vertical-velocity kinetic energy peaks at scales from hundreds of meters in the lower cloud layer to several kilometers higher up; these scales are typical grid sizes of today’s CRMs. The analysis also shows that a significant portion of the scalar and momentum fluxes in the benchmark simulation are carried by motions smaller than several kilometers (i.e., smaller than a typical grid resolution of CRMs). The broad range of scales of the benchmark simulation is split into two components: filter scale (mimicking CRM resolvable scale) and subfilter scale (mimicking CRM SGS), using filter widths characteristic of a typical CRM grid spacing. The local relationship of the subfilter-scale fluxes to the filter-scale variables is examined. This leads to a mixed SGS scheme to represent the SGS fluxes of scalars and momentum in CRMs. A priori tests show that the mixed SGS scheme yields spatial distributions of subfilter-scale fluxes that correlate much better with those retrieved from the benchmark when compared with an eddy viscosity/diffusivity scheme that is commonly used in today’s CRMs.

1970 ◽  
Vol 30 ◽  
pp. 19-31
Author(s):  
M Ashraf Uddin ◽  
M Matiar Rahman ◽  
M Saiful Islam Mallik

Generation of grid-scale (GS) and subgrid-scale (SGS) velocity fields is performed by direct filtering of DNS (Direct Numerical Simulation) data at a low Reynolds number in homogeneous isotropic turbulence in order to assess the spectral accuracy as well as the performance of filter functions for LES (Large Eddy Simulation). The filtering is performed using three classical filter functions: Gaussian, Tophat and Sharp cutoff filters and in all three cases the results are compared with three different filter widths for LES. Comparing the distributions of GS and SGS velocities, and the decay of turbulence with those from DNS fields through out the whole calculation we have found that among the three filter functions, the performance of Sharp cutoff filter is better than that of the other two filter functions in terms of both spatial spectra and the distribution of velocities. Furthermore, it is shown that the accuracy of the filtering approach does not depend only on the filter functions but also on the filter widths for LES. GANIT J. Bangladesh Math. Soc. (ISSN 1606-3694) 30 (2010) 19-31   DOI: http://dx.doi.org/10.3329/ganit.v30i0.8499


2011 ◽  
Vol 679 ◽  
pp. 156-193 ◽  
Author(s):  
EZGI S. TAŞKINOĞLU ◽  
JOSETTE BELLAN

For flows at supercritical pressure, p, the large-eddy simulation (LES) equations consist of the differential conservation equations coupled with a real-gas equation of state, and the equations utilize transport properties depending on the thermodynamic variables. Compared to previous LES models, the differential equations contain not only the subgrid-scale (SGS) fluxes but also new SGS terms, each denoted as a ‘correction’. These additional terms, typically assumed null for atmospheric pressure flows, stem from filtering the differential governing equations and represent differences, other than contributed by the convection terms, between a filtered term and the same term computed as a function of the filtered flow field. In particular, the energy equation contains a heat-flux correction (q-correction) which is the difference between the filtered divergence of the molecular heat flux and the divergence of the molecular heat flux computed as a function of the filtered flow field. We revisit here a previous a priori study where we only had partial success in modelling the q-correction term and show that success can be achieved using a different modelling approach. This a priori analysis, based on a temporal mixing-layer direct numerical simulation database, shows that the focus in modelling the q-correction should be on reconstructing the primitive variable gradients rather than their coefficients, and proposes the approximate deconvolution model (ADM) as an effective means of flow field reconstruction for LES molecular heat-flux calculation. Furthermore, an a posteriori study is conducted for temporal mixing layers initially containing oxygen (O) in the lower stream and hydrogen (H) or helium (He) in the upper stream to examine the benefit of the new model. Results show that for any LES including SGS-flux models (constant-coefficient gradient or scale-similarity models; dynamic-coefficient Smagorinsky/Yoshizawa or mixed Smagorinsky/Yoshizawa/gradient models), the inclusion of the q-correction in LES leads to the theoretical maximum reduction of the SGS molecular heat-flux difference; the remaining error in modelling this new subgrid term is thus irreducible. The impact of the q-correction model first on the molecular heat flux and then on the mean, fluctuations, second-order correlations and spatial distribution of dependent variables is also demonstrated. Discussions on the utilization of the models in general LES are presented.


Author(s):  
Hao Zhou ◽  
Qijing Feng ◽  
Pengcheng Hao ◽  
Zhiwei He ◽  
Li Li

This paper focuses on large eddy simulation of the Richtmyer–Meshkov instability (RMI) in spherical and cylindrical converging geometries with a Mach number [Formula: see text] based on subgrid-scale (SGS) dissipation similar method (SDSM). Based on the converging RMI problem, we obtain from a priori test and theoretical analysis that the suggested method can provide accurate structural correlation while ensuring the computational stability. Comparing the numerical simulation results with direct numerical simulation (DNS) and existing model in converging RMI problem, we could find that the suggested method overcomes some defects of the existing model, such as the Smagorinsky model cannot predict transition accurately and the helicity model can only predicts the quasi-two-dimensional flow precisely. It provides a beneficial tool for the research of converging RMI.


2014 ◽  
Vol 142 (2) ◽  
pp. 703-715 ◽  
Author(s):  
Chin-Hoh Moeng

Abstract A closure relationship between subgrid-scale (SGS) updraft–downdraft differences and resolvable-scale (RS) variables is proposed and tested for cloud-resolving models (CRMs), based on a data analysis of a large-eddy simulation (LES) of deep convection. The LES flow field is partitioned into CRM-RS and CRM-SGS using a cutoff scale that corresponds to a typical CRM grid resolution. This study first demonstrates the capability of an updraft–downdraft model framework in representing the SGS fluxes of heat, moisture, and momentum over the entire deep convection layer. It then formulates a closure scheme to relate SGS updraft–downdraft differences to horizontal gradients of RS variables. The closure is based on the idea that largest SGS and smallest RS motions are spectrally linked and hence their horizontal fluctuations must be strongly communicated. This relation leads to an SGS scheme that expresses vertical SGS fluxes in terms of horizontal gradients of RS variables, which differs from conventional downgradient eddy diffusivity models. The new scheme is shown to better represent the forward and backscatter energy transfer between CRM-RS and CRM-SGS components than conventional eddy-viscosity models.


Author(s):  
Shiwei Sun ◽  
Bowen Zhou ◽  
Ming Xue ◽  
Kefeng Zhu

AbstractIn numerical simulations of deep convection at kilometer-scale horizontal resolutions, in-cloud subgrid-scale (SGS) turbulence plays an important role in the transport of heat, moisture and other scalars. By coarse-graining a 50 m high-resolution large-eddy simulation (LES) of an idealized supercell storm to kilometer-scale grid spacings ranging from 250 m to 4 km, the SGS fluxes of heat, moisture, cloud and precipitating water contents are diagnosed a priori. The kilometer-scale simulations are shown to be within the “gray zone” as in-cloud SGS turbulent fluxes are comparable in magnitude to the resolved fluxes at 4 km spacing, and do not become negligible until ~500 m spacing. Vertical and horizontal SGS fluxes are of comparable magnitudes, both exhibit non-local characteristics associated with deep convection as opposed to local gradient-diffusion type of turbulent mixing. As such, they are poorly parameterized by eddy-diffusivity-based closures. To improve the SGS representation of turbulent fluxes in deep convective storms, a scale-similarity LES closure is adapted to kilometer-scale simulations. The model exhibits good correlations with LES-diagnosed SGS fluxes, and is capable of representing counter-gradient fluxes. In a posteriori tests, supercell storms simulated with the refined similarity closure model at kilometer-scale resolutions show better agreement with the LES benchmark in terms of SGS fluxes than those with a turbulent-kinetic-energy-based gradient-diffusion scheme. However, it underestimates the strength of updraft, which is suggested to be a consequence of the model effective resolution being lower than the native grid resolution.


2021 ◽  
Vol 1802 (4) ◽  
pp. 042088
Author(s):  
Zhipeng Feng ◽  
Huanhuan Qi ◽  
Xuan Huang ◽  
Shuai Liu ◽  
Jian Liu

Sign in / Sign up

Export Citation Format

Share Document