Direct numerical simulation of stationary homogeneous stratified sheared turbulence

2012 ◽  
Vol 696 ◽  
pp. 434-467 ◽  
Author(s):  
D. Chung ◽  
G. Matheou
Keyword(s):  
Numerical Simulation ◽  
Direct Numerical Simulation ◽  
Large Scale ◽  
Initial Conditions ◽  
Similarity Theory ◽  
Stationary Flow ◽  
Stratified Turbulence ◽  
Diapycnal Diffusivity ◽  
Turbulence Analysis ◽  
Theory Framework

AbstractUsing direct numerical simulation, we investigate stationary and homogeneous shear-driven turbulence in various stratifications, ranging from neutral to very stable. To attain and maintain a stationary flow, we throttle the mean shear so that the net production stays constant for all times. This results in a flow that is characterized solely by its mean shear and its mean buoyancy gradient, independent of initial conditions. The method of throttling is validated by comparison with experimental spectra in the case of neutral stratification. With increasing stratification comes the emergence of vertically sheared large-scale horizontal motions that preclude a straightforward interpretation of flow statistics. However, once these motions are excluded, simply by subtracting the horizontal average, the underlying flow appears amenable to the standard methods of turbulence analysis. It is shown that a direct acknowledgement of the confining influence of the periodic simulation box can lead to a meaningful physical interpretation of the large scales. Once an appropriate confinement scale is identified, many features, including horizontal spectra, flux–gradient relationships and length scales, of stratified sheared turbulence can be readily understood, both qualitatively and quantitatively, in terms of Monin–Obukhov similarity theory. Finally, the similarity-theory framework is used to interpret the scaling of the vertical diapycnal diffusivity in stratified turbulence.

Direct Numerical Simulation of Partial Fuel Stratification Assisted Lean Premixed Combustion for Assessment of Hybrid G-Equation/Well-Stirred Reactor Model

2021 ◽  
Author(s):  
Chao Xu ◽  
Muhsin Ameen ◽  
Pinaki Pal ◽  
Sibendu Som
Keyword(s):  
Numerical Simulation ◽  
Direct Numerical Simulation ◽  
Initial Conditions ◽  
Operating Conditions ◽  
Premixed Combustion ◽  
Lean Premixed Combustion ◽  
Fuel Stratification ◽  
Stirred Reactor ◽  
Lean Premixed ◽  
Flame Structures

Abstract Partial fuel stratification (PFS) is a promising fuel injection strategy to stabilize lean premixed combustion in spark-ignition (SI) engines. PFS creates a locally stratified mixture by injecting a fraction of the fuel, just before spark timing, into the engine cylinder containing homogeneous lean fuel/air mixture. This locally stratified mixture, when ignited, results in complex flame structure and propagation modes similar to partially premixed flames, and allows for faster and more stable flame propagation than a homogeneous lean mixture. This study focuses on understanding the detailed flame structures associated with PFS-assisted lean premixed combustion. First, a two-dimensional direct numerical simulation (DNS) is performed using detailed fuel chemistry, experimental pressure trace, and realistic initial conditions mapped from a prior engine large-eddy simulation (LES), replicating practical lean SI operating conditions. DNS results suggest that conventional triple flame structures are prevalent during the initial stage of flame kernel growth. Both premixed and non-premixed combustion modes are present with the premixed mode contributing dominantly to the total heat release. Detailed analysis reveals the effects of flame stretch and fuel pyrolysis on the flame displacement speed. Based on the DNS findings, the accuracy of a hybrid G-equation/well-stirred reactor (WSR) combustion model is assessed for PFS-assisted lean operation in the LES context. The G-equation model qualitatively captures the premixed branches of the triple flame, while the WSR model predicts the non-premixed branch of the triple flame. Finally, potential needs for improvements to the hybrid G-equation/WSR modeling approach are discussed.

Coherent structures in direct numerical simulation of turbulent boundary layers at Mach 3

2007 ◽  
Vol 594 ◽  
pp. 59-69 ◽  
Author(s):  
MATTHEW J. RINGUETTE ◽  
MINWEI WU ◽  
M. PINO MARTÍN
Keyword(s):  
Numerical Simulation ◽  
Direct Numerical Simulation ◽  
Boundary Layers ◽  
Coherent Structures ◽  
Mass Flux ◽  
Large Scale ◽  
Turbulent Boundary Layers ◽  
Turbulence Structure ◽  
Taylor’S Hypothesis ◽  
Mach 3

We demonstrate that data from direct numerical simulation of turbulent boundary layers at Mach 3 exhibit the same large-scale coherent structures that are found in supersonic and subsonic experiments, namely elongated, low-speed features in the logarithmic region and hairpin vortex packets. Contour plots of the streamwise mass flux show very long low-momentum structures in the logarithmic layer. These low-momentum features carry about one-third of the turbulent kinetic energy. Using Taylor's hypothesis, we find that these structures prevail and meander for very long streamwise distances. Structure lengths on the order of 100 boundary layer thicknesses are observed. Length scales obtained from correlations of the streamwise mass flux severely underpredict the extent of these structures, most likely because of their significant meandering in the spanwise direction. A hairpin-packet-finding algorithm is employed to determine the average packet properties, and we find that the Mach 3 packets are similar to those observed at subsonic conditions. A connection between the wall shear stress and hairpin packets is observed. Visualization of the instantaneous turbulence structure shows that groups of hairpin packets are frequently located above the long low-momentum structures. This finding is consistent with the very large-scale motion model of Kim & Adrian (1999).

scholarly journals Anisotropic Characteristics of Turbulence Dissipation in Swirling Flow: A Direct Numerical Simulation Study

2015 ◽  
Vol 2015 ◽  
pp. 1-9 ◽  
Author(s):  
Xingtuan Yang ◽  
Nan Gui ◽  
Gongnan Xie ◽  
Jie Yan ◽  
Jiyuan Tu ◽  
...  
Keyword(s):  
Numerical Simulation ◽  
Energy Dissipation ◽  
Direct Numerical Simulation ◽  
Large Scale ◽  
Isotropic Turbulence ◽  
Vortex Breakdown ◽  
Swirling Flow ◽  
Turbulent Energy ◽  
Swirling Flows ◽  
Small Scale

This study investigates the anisotropic characteristics of turbulent energy dissipation rate in a rotating jet flow via direct numerical simulation. The turbulent energy dissipation tensor, including its eigenvalues in the swirling flows with different rotating velocities, is analyzed to investigate the anisotropic characteristics of turbulence and dissipation. In addition, the probability density function of the eigenvalues of turbulence dissipation tensor is presented. The isotropic subrange of PDF always exists in swirling flows relevant to small-scale vortex structure. Thus, with remarkable large-scale vortex breakdown, the isotropic subrange of PDF is reduced in strongly swirling flows, and anisotropic energy dissipation is proven to exist in the core region of the vortex breakdown. More specifically, strong anisotropic turbulence dissipation occurs concentratively in the vortex breakdown region, whereas nearly isotropic turbulence dissipation occurs dispersively in the peripheral region of the strong swirling flows.

Identification of Bypass Transition Onset Markers Using Direct Numerical Simulation

2018 ◽  
Vol 140 (11) ◽  
Author(s):  
Shanti Bhushan ◽  
D. Keith Walters ◽  
S. Muthu ◽  
Crystal L. Pasiliao
Keyword(s):  
Numerical Simulation ◽  
Boundary Layer ◽  
Direct Numerical Simulation ◽  
Large Scale ◽  
General Purpose ◽  
Navier Stokes ◽  
Bypass Transition ◽  
Flow Parameters ◽  
Rans Models ◽  
Large Scale Flow

Efficacy of several large-scale flow parameters as transition onset markers are evaluated using direct numerical simulation (DNS) of boundary layer bypass transition. Preliminary results identify parameters (k2D/ν) and u′/U∞ to be a potentially reliable transition onset marker, and their critical values show less than 15% variation in the range of Re and turbulence intensity (TI). These parameters can be implemented into general-purpose physics-based Reynolds-averaged Navier–Stokes (RANS) models for engineering applications.

scholarly journals Direct numerical simulation and large-eddy simulation of stationary buoyancy-driven turbulence

2009 ◽  
Vol 643 ◽  
pp. 279-308 ◽  
Author(s):  
D. CHUNG ◽  
D. I. PULLIN
Keyword(s):  
Numerical Simulation ◽  
Large Eddy Simulation ◽  
Direct Numerical Simulation ◽  
Density Ratio ◽  
Stationary Flow ◽  
Small Scale ◽  
Scalar Flux ◽  
Eddy Simulation ◽  
Large Eddy ◽  
Active Scalar

We report direct numerical simulation (DNS) and large-eddy simulation (LES) of statistically stationary buoyancy-driven turbulent mixing of an active scalar. We use an adaptation of the fringe-region technique, which continually supplies the flow with unmixed fluids at two opposite faces of a triply periodic domain in the presence of gravity, effectively maintaining an unstably stratified, but statistically stationary flow. We also develop a new method to solve the governing equations, based on the Helmholtz–Hodge decomposition, that guarantees discrete mass conservation regardless of iteration errors. Whilst some statistics were found to be sensitive to the computational box size, we show, from inner-scaled planar spectra, that the small scales exhibit similarity independent of Reynolds number, density ratio and aspect ratio. We also perform LES of the present flow using the stretched-vortex subgrid-scale (SGS) model. The utility of an SGS scalar flux closure for passive scalars is demonstrated in the present active-scalar, stably stratified flow setting. The multi-scale character of the stretched-vortex SGS model is shown to enable extension of some second-order statistics to subgrid scales. Comparisons with DNS velocity spectra and velocity-density cospectra show that both the resolved-scale and SGS-extended components of the LES spectra accurately capture important features of the DNS spectra, including small-scale anisotropy and the shape of the viscous roll-off.

scholarly journals A study of the flow-field evolution and mixing in a planar turbulent jet using direct numerical simulation

2002 ◽  
Vol 450 ◽  
pp. 377-407 ◽  
Author(s):  
S. A. STANLEY ◽  
S. SARKAR ◽  
J. P. MELLADO
Keyword(s):  
Numerical Simulation ◽  
Flow Field ◽  
Direct Numerical Simulation ◽  
Large Scale ◽  
Computational Study ◽  
Practical Interest ◽  
Small Scale ◽  
Mixing Process ◽  
Mean Velocity ◽  
Small Scales

Turbulent plane jets are prototypical free shear flows of practical interest in propulsion, combustion and environmental flows. While considerable experimental research has been performed on planar jets, very few computational studies exist. To the authors' knowledge, this is the first computational study of spatially evolving three-dimensional planar turbulent jets utilizing direct numerical simulation. Jet growth rates as well as the mean velocity, mean scalar and Reynolds stress profiles compare well with experimental data. Coherency spectra, vorticity visualization and autospectra are obtained to identify inferred structures. The development of the initial shear layer instability, as well as the evolution into the jet column mode downstream is captured well.The large- and small-scale anisotropies in the jet are discussed in detail. It is shown that, while the large scales in the flow field adjust slowly to variations in the local mean velocity gradients, the small scales adjust rapidly. Near the centreline of the jet, the small scales of turbulence are more isotropic. The mixing process is studied through analysis of the probability density functions of a passive scalar. Immediately after the rollup of vortical structures in the shear layers, the mixing process is dominated by large-scale engulfing of fluid. However, small-scale mixing dominates further downstream in the turbulent core of the self-similar region of the jet and a change from non-marching to marching PDFs is observed. Near the jet edges, the effects of large-scale engulfing of coflow fluid continue to influence the PDFs and non-marching type behaviour is observed.

Sign in / Sign up

Export Citation Format

Share Document