Evaluation of strain-rate effects in transitional round jets using direct numerical simulation

The coherent flame model uses the strain rate to predict reaction rate per unit flame surface area and some procedure that solves for the dynamics of flame surfaces to predict species distributions. The strainrate formula for the reaction rate is obtained from the analytical solution for a flame in a laminar, plane stagnation point flow. Here, the formula's effectiveness is examined by comparisons with data from a direct numerical simulation (DNS) of a round jetlike flow that undergoes transition to turbulence. Significant differences due to general flow features can be understood qualitatively: Model predictions are good in the braids between vortex rings, which are present in the near field of round jets, as the strain rate is extensional and reaction surfaces are isolated. In several other regions, the strain rate is compressive or flame surfaces are folded close together. There, the predictions are poor as the local flow no longer resembles the model flow. Quantitative comparisons showed some discrepancies. A modified, consistent application of the strain-rate solution did not show significant changes in the prediction of mean reaction rate distributions.

Direct Numerical Simulation Of Three-Dimensional Richtmyer-Meshkov Instability

Direct numerical simulation (DNS) is used to study flow characteristics after interaction of a planar shock with a spherical media interface in each side of which the density is different. This interfacial instability is known as the Richtmyer-Meshkov (R-M) instability. The compressible Navier-Stoke equations are discretized with group velocity control (GVC) modified fourth order accurate compact difference scheme. Three-dimensional numerical simulations are performed for R-M instability installed passing a shock through a spherical interface. Based on numerical results the characteristics of 3D R-M instability are analysed. The evaluation for distortion of the interface, the deformation of the incident shock wave and effects of refraction, reflection and diffraction are presented. The effects of the interfacial instability on produced vorticity and mixing is discussed.

Reynolds number influence on turbulent blocking effects of a rigid plane boundary (2nd report, direct numerical simulation using townsend's eddy forcing)

The Reynolds number influence on turbulent blocking effects by a rigid plane boundary is studied using direct numerical simulation (DNS). A new forcing method using 'simple model eddies' (Townsend 1976) for DNS of stationary homogeneous isotropic turbulence is proposed. A force field is obtained in real space by sprinkling many space-filling 'simple model eddies' whose centers are randomly but uniformly distributed in space and whose axes of rotation are random. The method is applied to a shear-free turbulent boundary layer over a rigid plane boundary and the blocking effects are investigated. The results show that stationary homogeneous isotropic turbulence is generated in real space using the present method. By using different model eddies with different sizes and rotation speeds, we could change the turbulence properties such as the integral and micro scales, the turbulent Reynolds number and the isotropy of turbulence. Turbulence intensities near the wall showed good agreements with the previous measurement and the linear analysis based on a rapid distortion theory (RDT). The splat effect (i.e., turbulence intensities of the components parallel to the boundary are amplified) occurs near the boundary and the viscous effect prohibits the splat effect at the quasi steady state at low Reynolds number.

Reynolds number influence on turbulent blocking effects of a rigid plane boundary (3rd report, direct numerical simulation on wall blocking effects for anisotropic turbulence)

The Reynolds number influence on turbulent blocking effects by a rigid plane boundary is studied using direct numerical simulation (DNS). A new forcing method proposed in the second report using Townsend's "simple model eddies" for DNS was extended to generate axisymmetric anisotropic turbulence. A force field is obtained in real space by sprinkling many space-filling "simple model eddies" whose centers are randomly but uniformly distributed in space. The axes of rotation are controlled in this study to generate axisymmetric anisotropic turbulence. The method is applied to a shear-free turbulent boundary layer over a rigid plane boundary and the blocking effects for anisotropic turbulence are investigated. The results show that stationary axisymmetric anisotropic turbulence is generated using the present method. Turbulence intensities near the wall showed good agreements with the rapid distortion theory (RDT) for small t (t ≪ TL), where TL. is the eddy turnover time. The splat effect (i. e. turbulence intensities of the components parallel to the surface are amplified) occurs near the boundary and the viscous effect attenuates the splat effect at the quasi steady state at low Reynolds number as for Isotropic turbulence. Prandtl's secondary flow of the second kind does not occur for low Reynolds number flows, which qualitatively agrees with previous observetion in a mixing-box.

Effects of turbulent reynolds number on the performance of algebraic flame surface density models for large eddy simulation in the thin reaction zones regime: A direct numerical simulation analysis

A direct numerical simulation (DNS) database of freely propagating statistically planar turbulent premixed flames with a range of different turbulent Reynolds numbers has been used to assess the performance of algebraic flame surface density (FSD) models based on a fractal representation of the flame wrinkling factor. The turbulent Reynolds number Ret has been varied by modifying the Karlovitz number Ka and the Damköhler number Da independently of each other in such a way that the flames remain within the thin reaction zones regime. It has been found that the turbulent Reynolds number and the Karlovitz number both have a significant influence on the fractal dimension, which is found to increase with increasing Ret and Ka before reaching an asymptotic value for large values of Ret and Ka. A parameterisation of the fractal dimension is presented in which the effects of the Reynolds and the Karlovitz numbers are explicitly taken into account. By contrast, the inner cut-off scale normalised by the Zel'dovich flame thickness ηi/δz does not exhibit any significant dependence on Ret for the cases considered here. The performance of several algebraic FSD models has been assessed based on various criteria. Most of the algebraic models show a deterioration in performance with increasing the LES filter width. © 2012 Mohit Katragadda et al.

Direct numerical simulation of turbulent flow and mass transfer around a rotating circular cylinder

Turbulent flow around a rotating circular cylinder has numerous applications including wall shear stress and mass-transfer measurement related to the corrosion studies. It is also of interest in the context of flow over convex surfaces where standard turbulence models perform poorly. The main purpose of this paper is to elucidate the basic turbulence mechanism around a rotating cylinder at low Reynolds numbers to provide a better understanding of flow fundamentals. Direct numerical simulation (DNS) has been performed in a reference frame rotating at constant angular velocity with the cylinder. The governing equations are discretized by using a finite-volume method. As for fully developed channel, pipe, and boundary layer flows, a laminar sublayer, buffer layer, and logarithmic outer region were observed. The level of mean velocity is lower in the buffer and outer regions but the logarithmic region still has a slope equal to the inverse of the von Karman constant. Instantaneous flow visualization revealed that the turbulence length scale typically decreases as the Reynolds number increases. Wavelet analysis provided some insight into the dependence of structural characteristics on wave number. The budget of the turbulent kinetic energy was computed and found to be similar to that in plane channel flow as well as in pipe and zero pressure gradient boundary layer flows. Coriolis effects show as an equivalent production for the azimuthal and radial velocity fluctuations leading to their ratio being lowered relative to similar nonrotating boundary layer flows.

Geometrical properties and turbulent flame speed measurements in stationary premixed V-flames using direct numerical simulation

Three dimensional, fully compressible direct numerical simulations (DNS) of premixed turbulent flames are carried out in a V-flame configuration. The governing equations and the numerical implementation are described in detail, including modifications made to the Navier-Stokes Characteristic Boundary Conditions (NSCBC) to accommodate the steep transverse velocity and composition gradients generated when the flame crosses the boundary. Three cases, at turbulence intensities, u′/sL, of 1, 2, and 6 are considered. The influence of the flame holder on downstream flame properties is assessed through the distributions of the surface-conditioned displacement speed, curvature and tangential strain rates, and compared to data from similarly processed planar flames. The distributions are found to be indistinguishable from planar flames for distances greater than about 17δth downstream of the flame holder, where δth is the laminar flame thermal thickness. Favre mean fields are constructed, and the growth of the mean flame brush is found to be well described by simple Taylor type diffusion. The turbulent flame speed, sT is evaluated from an expression describing the propagation speed of an isosurface of the mean reaction progress variable c̃ in terms of the imbalance between the mean reactive, diffusive, and turbulent fluxes within the flame brush. The results are compared to the consumption speed, sC, calculated from the integral of the mean reaction rate, and to the predictions of a recently developed flame speed model (Kolla et al., Combust Sci Technol 181(3):518-535, 2009). The model predictions are improved in all cases by including the effects of mean molecular diffusion, and the overall agreement is good for the higher turbulence intensity cases once the tangential convective flux of c̃ is taken into account. © 2010 Springer Science+Business Media B.V.

Direct Numerical Simulation of Turbulent Flames

The present paper is a report on progress in the simulation of turbulent flames using the Cray T3D and T3E at the Edinburgh parallel computing centre, using codes developed in Cambridge. Two combustion DNS codes are described, ANGUS and SENGA, which solve incompressible and fully compressible reacting flows respectively. The technical background to combustion DNS is presented, and the resource requirements explained in terms of the physic and chemistry of the problem. Results for flame turbulence interaction studies are presented and discussed in terms of their relevance to modelling. Recent work on the fully compressible problem is highlighted and future directions outlined.

Statistically steady turbulence in thin films: direct numerical simulations with Ekman friction

We present a detailed direct numerical simulation (DNS) of the two-dimensional Navier-Stokes equation with the incompressibility constraint and air-drag-induced Ekman friction; our DNS has been designed to investigate the combined effects of walls and such a friction on turbulence in forced thin films. We concentrate on the forward-cascade regime and show how to extract the isotropic parts of velocity and vorticity structure functions and hence the ratios of multiscaling exponents. We find that velocity structure functions display simple scaling, whereas their vorticity counterparts show multiscaling, and the probability distribution function of the Weiss parameter 3, which distinguishes between regions with centers and saddles, is in quantitative agreement with experiments.

Direct numerical simulation of autoignition in a non-premixed, turbulent medium

In this paper, direct numerical simulation of autoignition in an initially non-premixed medium under isotropic, homogeneous, and decaying turbulence is presented. The pressure-based method developed herein is a spectral implementation of the sequential steps followed in the predictor-corrector type of algorithms; it includes the effects of density fluctuations caused by spatial inhomogeneities ill temperature and species. The velocity and pressure field are solved in the spectral space while the scalars and density field are solved in the physical space. The presented results reveal that the autoignition spots originate and evolve at locations where (1) the composition corresponds to a small range around a specific mixture fraction, and (2) the conditional scaler dissipation rate is low. A careful examination of the data obtained indicates that the autoignition spots originate in the vortex cores, and the hot gases travel outward as combustion progresses. Hence, the applicability of the transient laminar flamelet model for this problem is questioned. The dependence of autoignition characteristics on parameters such as (1) die initial eddy-turnover time and (2) the initial ratio of length scale of scalars to that of velocities are investigated. Certain implications of new results on the conditional moment closure modeling are discussed.