A strategy for large-scale scalar advection in large eddy simulations that use the linear eddy sub-grid mixing model

2018 ◽  
Vol 28 (10) ◽  
pp. 2463-2479 ◽  
Author(s):  
Salman Arshad ◽  
Bo Kong ◽  
Alan Kerstein ◽  
Michael Oevermann

PurposeThe purpose of this numerical work is to present and test a new approach for large-scale scalar advection (splicing) in large eddy simulations (LES) that use the linear eddy sub-grid mixing model (LEM) called the LES-LEM.Design/methodology/approachThe new splicing strategy is based on an ordered flux of spliced LEM segments. The principle is that low-flux segments have less momentum than high-flux segments and, therefore, are displaced less than high-flux segments. This strategy affects the order of both inflowing and outflowing LEM segments of an LES cell. The new splicing approach is implemented in a pressure-based fluid solver and tested by simulation of passive scalar transport in a co-flowing turbulent rectangular jet, instead of combustion simulation, to perform an isolated investigation of splicing. Comparison of the new splicing with a previous splicing approach is also done.FindingsThe simulation results show that the velocity statistics and passive scalar mixing are correctly predicted using the new splicing approach for the LES-LEM. It is argued that modeling of large-scale advection in the LES-LEM via splicing is reasonable, and the new splicing approach potentially captures the physics better than the old approach. The standard LES sub-grid mixing models do not represent turbulent mixing in a proper way because they do not adequately represent molecular diffusion processes and counter gradient effects. Scalar mixing in turbulent flow consists of two different processes, i.e. turbulent mixing that increases the interface between unmixed species and molecular diffusion. It is crucial to model these two processes individually at their respective time scales. The LEM explicitly includes both of these processes and has been used successfully as a sub-grid scalar mixing model (McMurtry et al., 1992; Sone and Menon, 2003). Here, the turbulent mixing capabilities of the LES-LEM with a modified splicing treatment are examined.Originality/valueThe splicing strategy proposed for the LES-LEM is original and has not been investigated before. Also, it is the first LES-LEM implementation using unstructured grids.

2014 ◽  
Vol 136 (9) ◽  
Author(s):  
Malcolm J. Andrews ◽  
David L. Youngs ◽  
Daniel Livescu ◽  
Tie Wei

A time-dependent, incompressible, turbulent mixing problem, referred here to as the “tilted-rig,” is defined, based results from an experiment that involved the introduction of a large-scale overturning motion, with a superposed localized Rayleigh-Taylor (RT) driven mixing. The problem serves to examine the development of RT turbulent mixing while being strained by a large-scale two-dimensional confined motion. Care is taken to define the problem in detail so others might use the definition, and the results, to help develop advanced models of buoyancy driven mixing in complex flows. Aside from a careful definition, the problem has been solved using two different implicit-large-Eddy-simulations (ILES) based codes, and with a direct numerical simulations (DNS) code. Two-dimensional and one-dimensional mix metrics are defined, and then used to examine the development of the mixing region, and the overall evolution of the flow. Comparison of simulations with experiment reveals that large-scale overturning can be well captured in all the simulations, similarly central mix widths, and spike/bubble sidewall penetrations are also in good agreement. A comparison between the different simulation methodologies, ILES and DNS, reveals an overall good agreement between mix metrics such as the amount of molecular mixing. The DNS simulations reveal a dependency on Reynolds number that merits further experimental work.


2021 ◽  
Author(s):  
Gaston Latessa ◽  
Angela Busse ◽  
Manousos Valyrakis

<p>The prediction of particle motion in a fluid flow environment presents several challenges from the quantification of the forces exerted by the fluid onto the solids -normally with fluctuating behaviour due to turbulence- and the definition of the potential particle entrainment from these actions. An accurate description of these phenomena has many practical applications in local scour definition and to the design of protection measures.</p><p>In the present work, the actions of different flow conditions on sediment particles is investigated with the aim to translate these effects into particle entrainment identification through analytical solid dynamic equations.</p><p>Large Eddy Simulations (LES) are an increasingly practical tool that provide an accurate representation of both the mean flow field and the large-scale turbulent fluctuations. For the present case, the forces exerted by the flow are integrated over the surface of a stationary particle in the streamwise (drag) and vertical (lift) directions, together with the torques around the particle’s centre of mass. These forces are validated against experimental data under the same bed and flow conditions.</p><p>The forces are then compared against threshold values, obtained through theoretical equations of simple motions such as rolling without sliding. Thus, the frequency of entrainment is related to the different flow conditions in good agreement with results from experimental sediment entrainment research.</p><p>A thorough monitoring of the velocity flow field on several locations is carried out to determine the relationships between velocity time series at several locations around the particle and the forces acting on its surface. These results a relevant to determine ideal locations for flow investigation both in numerical and physical experiments.</p><p>Through numerical experiments, a large number of flow conditions were simulated obtaining a full set of actions over a fixed particle sitting on a smooth bed. These actions were translated into potential particle entrainment events and validated against experimental data. Future work will present the coupling of these LES models with Discrete Element Method (DEM) models to verify the entrainment phenomena entirely from a numerical perspective.</p>


2021 ◽  
Author(s):  
Gregory Wagner ◽  
Andre Souza ◽  
Adeline Hillier ◽  
Ali Ramadhan ◽  
Raffaele Ferrari

<p>Parameterizations of turbulent mixing in the ocean surface boundary layer (OSBL) are key Earth System Model (ESM) components that modulate the communication of heat and carbon between the atmosphere and ocean interior. OSBL turbulence parameterizations are formulated in terms of unknown free parameters estimated from observational or synthetic data. In this work we describe the development and use of a synthetic dataset called the “LESbrary” generated by a large number of idealized, high-fidelity, limited-area large eddy simulations (LES) of OSBL turbulent mixing. We describe how the LESbrary design leverages a detailed understanding of OSBL conditions derived from observations and large scale models to span the range of realistically diverse physical scenarios. The result is a diverse library of well-characterized “synthetic observations” that can be readily assimilated for the calibration of realistic OSBL parameterizations in isolation from other ESM model components. We apply LESbrary data to calibrate free parameters, develop prior estimates of parameter uncertainty, and evaluate model errors in two OSBL parameterizations for use in predictive ESMs.</p>


2013 ◽  
Vol 13 (12) ◽  
pp. 31891-31932 ◽  
Author(s):  
R. Paoli ◽  
O. Thouron ◽  
J. Escobar ◽  
J. Picot ◽  
D. Cariolle

Abstract. Large-eddy simulations of sub-kilometer-scale turbulence in the upper troposphere lower stratosphere (UTLS) are carried out and analyzed using the mesoscale atmospheric model Méso-NH. Different levels of turbulence are generated using a large-scale stochastic forcing technique that was especially devised to treat atmospheric stratified flows. The study focuses on the analysis of turbulence statistics, including mean quantities and energy spectra, as well as on a detailed description of flow topology. The impact of resolution is also discussed by decreasing the grid spacing to 2 m and increasing the number of grid points to 8×109. Because of atmospheric stratification, turbulence is substantially anisotropic, and large elongated structures form in the horizontal directions, in accordance with theoretical analysis and spectral direct numerical simulations of stably stratified flows. It is also found that the inertial range of horizontal kinetic energy spectrum, generally observed at scales larger than a few kilometers, is prolonged into the sub-kilometric range, down to the Ozmidov scales that obey isotropic Kolmorogov turbulence. The results are in line with observational analysis based on in situ measurements from existing campaigns.


2002 ◽  
Vol 455 ◽  
pp. 195-212 ◽  
Author(s):  
DANIELE CARATI ◽  
MICHAEL M. ROGERS ◽  
ALAN A. WRAY

A statistical ensemble of large-eddy simulations (LES) is run simultaneously for the same flow. The information provided by the different large-scale velocity fields is used in an ensemble-averaged version of the dynamic model. This produces local model parameters that only depend on the statistical properties of the flow. An important property of the ensemble-averaged dynamic procedure is that it does not require any spatial averaging and can thus be used in fully inhomogeneous flows. Also, the ensemble of LES provides statistics of the large-scale velocity that can be used for building new models for the subgrid-scale stress tensor. The ensemble-averaged dynamic procedure has been implemented with various models for three flows: decaying isotropic turbulence, forced isotropic turbulence, and the time-developing plane wake. It is found that the results are almost independent of the number of LES in the statistical ensemble provided that the ensemble contains at least 16 realizations.


Sign in / Sign up

Export Citation Format

Share Document