Structure of nonpremixed reaction zones in numerical isotropic turbulence

Results from direct numerical simulation of low heat release, turbulent nonpremixed reacting flows modeled using single-step reactions with constant and temperature-dependent rate laws are discussed, and compared with laminar predictions. The mixture fraction and its dissipation rate are statistically independent in regions of intense reaction, partially supporting a commonly made assumption in flamelet-based models. In the presence of a finite rate reverse reaction, the reaction zone spans the entire range of mixture fraction. The joint pdf of the reactive scalars evolves to an equilibrium that is dictated by a balance between the reactive and mixing fluxes in composition space. When the temperature-dependent rate law is implemented, strain-induced extinction is observed for a Zel'dovich (modified) number of 10. As the ratio of local flow to chemical time scale is decreased below unity, a larger fraction of the flow field experiences this mode of extinction. The critical turbulent scalar dissipation rate is compared with laminar values and asymptotic predictions.The first two authors express their acknowledgment to the donors of The Petroleum Research Fund, administered by the American Chemical Society for partial support of this work through a type G grant. We also acknowledge the Council for Research and Creative Writing at the University of Colorado for providing partial financial support.     

A Generalised Multiple Mapping Conditioning Approach for Turbulent Combustion

This paper follows the evolution in understanding of the multiple mapping conditioning (MMC) approach for turbulent combustion and reviews different implementations of MMC models. As the MMC name suggests, the original version represents a consistent combination of CMC-type conditional equations (conditional moment closure) and generalised mapping closure. It seems that the strength of the MMC model, and especially that of its stochastic version, lies in a more general (and much more transparent) interpretation. In this new generalised interpretation, we can replace complicated derivations by physical reasoning and the model appears to be a natural extension of modelling approaches developed in recent decades. MMC can be seen as a methodology for enforcing certain known characteristics of turbulence on a conventional mixing model. This is achieved by localising the mixing operation in a reference space. The reference space variables are selected to emulate the properties of a turbulent flow which have a strong effect on reactive quantities. The best and simplest example is an MMC model which has a single reference variable emulating the mixture fraction. In diffusion flames turbulent fluctuations of reacting quantities are strongly correlated with fluctuations of the mixture fraction. By making mixing local in the reference mixture fraction space a CMC-type mixing closure is enforced. In the original interpretation of MMC the reference variables are modelled as Markov processes. Since the reference variables should emulate properties of turbulent flows as realistically as possible the next step, and the basis of generalised MMC, is to remove the Markovian restriction and set reference variables equal to traced Lagrangian quantities within DNS or LES flow fields. Indeed, no Markov value can emulate the mixture fraction better than the mixture fraction itself. (Using a Markov vector process of dimension higher than the number of conditioning variables represents a more economical alternative for producing reference variables in generalised MMC.) The generalised MMC approach effectively incorporates the mixture fraction-based models, the PDF methods and LES/DNS techniques into a single methodology with possibility of blending useful features developed previously for conventional models. The generalised approach to MMC stimulates a more flexible understanding of simulations using sparsely placed Lagrangian particles as tools that may provide accurate joint distributions of reactive scalars at relatively low computational cost. The physical reasoning behind the new interpretation of MMC is supported by example computations for a partially premixed methane/air diffusion flame (Sandia Flame D). The scheme utilises LES for the dynamic field and a sparse-Lagrangian filtered density function method with MMC mixing for the scalar field. Two different particle mixing schemes are tested. Simulations are performed using only 35,000 Lagrangian particles (of these only 10,000 are chemically active) on a single workstation. The relatively low computational cost allows the use of realistic chemical kinetics containing 34 reactive species and 219 reactions. Intended for publication in the special issue of Flow, Turbulence and Combustion arising from the 2nd ECCOMAS Thematic Conference on Computational Combustion held at Delft in mid-2007.     

Composition-Space Behavior of Diffusion-Reaction Systems

Species concentrations in homogeneous well-stirred reacting (WSR) systems with multiple time scales mostly access only a small portion (slow manifold) of the permitted composition space. This offers a mathematically sound basis for reducing complex kinetics. Computational investigation of such slow manifolds in inhomogeneous reaction systems undergoing diffusion in addition to chemical reaction is the subject of this paper. We perform simulations of a one-dimensional diffusion–reaction system from random non-premixed initial conditions to study the accessed composition space for a variety of diffusion to reaction time scale ratios. A detailed hydrogen–oxygen chemical kinetic set (6 species, 14 reactions) is used. Our results show that, to first order, the effect of random advection and diffusion is merely to increase the accessed state-space dimensionality over that of the WSR by the number of relevant mixture fractions. Circumstances under which the above simplification is not accurate are also identified.     

Mixing Intensification on the Edge of a Jet by Means of Longitudinal Vortices

The action of an artificially generated spanwise flow in the form of periodical longitudinal vortices on a plane turbulent mixing layer is investigated. It is shown that the disturbances result in a significant increase in the thickness of the mixing region. For two kinds of spanwise flow, namely, vortices whose centers lie in the plane separating the streams and vortices located above this plane, the dependence of the mixing layer thickness on the vortex amplitude and vertical dimension and on the longitudinal coordinate is found.     

Stochastic Simulation of Scalar Mixing Capturing Unsteadiness and Small-scale Structure Based on Mean-flow Properties

A key limitation of Reynolds-Average Navier-Stokes (RANS) simulation of mixing and reaction in turbulent flows is the lack of resolution of small-scale structure and associated unsteadiness. Various subgrid models formulated in state space have been developed to complement the RANS solution in this regard. We here introduce a physical-space formulation that captures unsteady advective and diffusive processes at all scales of the turbulent flow. The approach is a 3D construction based on the Linear Eddy Model (LEM), involving three orthogonally intersecting arrays of 1D LEM domains, and coupled so as to capture the 3D character of fluid trajectories. To illustrate the model performance of the 3D LEM-based formulation, here termed LEM3D, multi-stream mixing in a turbulent round jet is simulated using measured mean-flow properties as input. Comparison to scalar cross-correlation coefficients and other measured properties of this mixing configuration indicate that the LEM3D approach, in conjunction with flow properties that are provided by steady-state models, is a useful representation of complex turbulent mixing processes that would otherwise be difficult to capture within a steady-state CFD framework.     

DNS study of passive plume interference emitting from two parallel line sources in a turbulent channel flow

Turbulent mixing of dual plumes emitting simultaneously from line sources in a turbulent channel flow has been studied using direct numerical simulation (DNS). Three test cases have been compared to investigate the effects of the source separation on turbulent mixing of the two instantaneous plumes. The dispersion and interference of dual plumes are investigated in both physical and spectral spaces, which include an analysis of statistical moments of the concentration field, cross-correlation between the two instantaneous plumes, pre-multiplied spectra of the velocity and concentration fields, and co-spectrum and coherency spectrum of the dual plumes. As the downstream distance from the line source increases, the plume development associated with a single source emission transitions from a turbulent convective stage to a turbulent diffusive stage. It is observed that a plume released from a ground-level source reaches the turbulent diffusive stage faster than that released from an elevated source. It is also observed that a smaller separation between the two line sources tends to facilitate a more rapid growth in the cross-correlation coefficient of two instantaneous plumes. In the near-source region, the maximum coherency spectrum is produced at lower frequencies indicating that dual-plume mixing is dominated by the external flapping effects of large-scale eddy motions. However, in the far downstream region of the sources, the coherency spectrum in the higher frequency range increases significantly, indicating that the spread of the total plume is larger than all scales of turbulent eddies, such that they all contribute to the in-plume mixing of the dual plumes.     

Impact of Turbulent Flow and Mean Mixture Fraction Results on Mixing Model Behavior in Transported Scalar PDF Simulations of Turbulent Non-premixed Bluff Body Flames

Numerical simulation results are presented for three turbulent jet diffusion flames, stabilized behind a bluff body (Sydney Flames HM1-3). Interaction between turbulence and combustion is modeled with the transported joint-scalar PDF approach. The focus of the study is on the impact of the quality of simulation results in physical space on the behavior of two micro-mixing models in composition space: the Euclidean Minimum Spanning Tree (‘EMST’) model and the modified Curl coalescence dispersion (‘CD’) model. Profiles of conditional means and variances of thermo-chemical quantities, conditioned on the mixture fraction, are discussed in the recirculation region and in the neck zone behind. The impact of the flow and mixing fields in physical space on the mixing model behavior in composition space is strong for the CD model and increases as the turbulence – chemistry interaction becomes stronger. The EMST conditional profiles, on the contrary, are hardly affected.     

Invariant analysis of the Reynolds stress tensor for a nuclear fuel assembly with spacer grid and split type vanes

Invariant analysis of the Reynolds stress tensor anisotropy can give an accurate and deep intuitive understanding of the turbulent structure of a turbulent flow. Lumley's triangle has proven to be a powerful representation of the invariant analysis of the second-order statistics collection provided by the Reynolds stress tensor. In the present work the spectral element code Nek5000 has been used to investigate the turbulent structure of the flow across a pressurized water reactor spacer grid with split type mixing vanes. Wall-resolved large eddy simulation of the flow in a prototypical rod bundle geometry at Re = 14,000 and P/D = 1.32 are performed and validated against particle image velocimetry data. The results are then used to perform an in-depth invariant analysis. The results show a reorganization of the Reynolds stresses components in the downstream region of the spacer grid. The mixing vanes orientation produces a symmetric behavior between sub-channels. The turbulent structure in the fully developed region has the typical behavior of fully-developed channel flow turbulence. When averaging the state across regions of the sub-channels, we observed a transition from disk-like turbulence in the mixing vanes region to rod-like turbulence in the fully developed region.     

竖直平板间自然对流的湍流频谱特性

根据直接数值模拟的计算结果，对竖直平板间湍流自然对流的脉动动能、速度及温度等物理量的时间序列进行频谱分析．结果表明，流动达到充分发展状态后，小尺度到大尺度的能谱很宽，计算的分辨率足够．从能谱分布可以观察到含能区、惯性子区和耗散区的存在，文中对各区的特性进行分析．由于该流动的强各向异性，惯性子区很窄．并讨论了法向位置对脉动动能的影响以及大尺度结构的特性．     

Simultaneous measurements of reactive scalar and velocity in a planar liquid jet with a second-order chemical reaction

This paper presents a new experimental approach for simultaneous measurements of velocity and concentration in a turbulent liquid flow with a chemical reaction. For the simultaneous measurements, we developed a combined probe consisting of an I-type hot-film probe and an optical fiber probe based on the light absorption spectrometric method. In a turbulent planar liquid jet with a second-order chemical reaction $(\mathrm{A}+\mathrm{B}\rightarrow \mathrm{R})$ , streamwise velocity and concentrations of all reactive species are measured by the combined probe. The turbulent mass fluxes of the reactive species are estimated from the simultaneous measurements. The results show that the influence of the chemical reaction on the turbulent mass flux of the reactant species near the jet exit is different from its influence in other regions, and the turbulent mass flux of the product species has a negative value near the jet exit and a positive value in other regions.     

Experimental investigation of ice slurry flow pressure drop in horizontal tubes

Pressure drop behaviour of ice slurry based on ethanol–water mixture in circular horizontal tubes has been experimentally investigated. The secondary fluid was prepared by mixing ethyl alcohol and water to obtain initial alcohol concentration of 10.3% (initial freezing temperature ?4.4 °C). The pressure drop tests were conducted to cover laminar and slightly turbulent flow with ice mass fraction varying from 0% to 30% depending on test conditions. Results from flow tests reveal much higher pressure drop for higher ice concentrations and higher velocities in comparison to the single phase flow. However for ice concentrations of 15% and higher, certain velocity exists at which ice slurry pressure drop is same or even lower than for single phase flow. It seems that higher ice concentration delay flow pattern transition moment (from laminar to turbulent) toward higher velocities. In addition experimental results for pressure drop were compared to the analytical results, based on Poiseulle and Buckingham–Reiner models for laminar flow, Blasius, Darby and Melson, Dodge and Metzner, Steffe and Tomita for turbulent region and general correlation of Kitanovski which is valid for both flow regimes. For laminar flow and low buoyancy numbers Buckingham–Reiner method gives good agreement with experimental results while for turbulent flow best fit is provided with Dodge–Metzner and Tomita methods.Furthermore, for transport purposes it has been shown that ice mass fraction of 20% offers best ratio of ice slurry transport capability and required pumping power.     

Schlieren measurement of turbulent structure in a diffusion flame

The regular and random mixing structures in a turbulent diffusion flame were investigated using the quantitative, dynamic crossed-beam schlieren method. Evidence was found close to the nozzle relating to the vortexlike structure of eddies surrounding the central fuel jet flow. The observations also make possible resolution of turbulent intensity, scales, convection, and spectra within the diffusion flame without the use of seeding or intrusion of measuring probes. It is found that length scales and other turbulence parameters in the diffusion flame progressively revert to values similar to those expected and observed in scalar passive mixing as the combustion reaction intensity reduces with axial distance from the nozzle system.     

A geometric approach for predicting vertical stationary profiles of weakly inertial advecting-diffusing particles in closed incompressible flows

Mixing of weakly inertial particles in closed flows is often addressed by considering individual particles as passive advecting-diffusing tracers, subjected to an additional settling velocity resulting from body forces (e.g. gravity). We show that the qualitative and quantitative features of the vertical particle distribution (i.e. the horizontal cross-sectional averages of particle concentration) can be predicted from the structure of the flow resulting from the superposition of the stirring field and the settling velocity. The prediction is based upon the observation that the resulting flow can be divided into two nonoverlapping regions, namely trajectories that are confined within the mixing space (recirculation loops), and trajectories that cross the mixing space. The spatial extent of these regions is exploited to define an effective vertical convective velocity entering the one-dimensional lumped model. Model two-dimensional flows possessing different flow patterns are used to illustrate the proposed estimate for effective velocity. A CFD-computed three-dimensional turbulent flow inside a baffled stirred vessel is used as a benchmark test to assess the model performance in typical industrial flows.     

条件矩模型模拟湍流扩散燃烧

对条件矩模型模拟湍流扩散燃烧进行了初步的研究.在条件矩模型中,标量的统计是以混合分数为条件的,条件平均使得非线性化学反应源项可以在一阶条件下被封闭.模拟结果和美国Sandia国家实验室的实验结果对比表明:对温度、主要组分浓度的预报结果是令人满意的,NO浓度的预报在趋势上也符合实验结果.误差分析表明,提高标量耗散率的预报精度和二阶条件矩模型都将有助于推动条件矩模型的发展.     

Flow and mixing characteristics of pulsed confined opposed jets in turbulent flow regime

A numerical study is performed on a two-dimensional confined opposed-jet configuration to gain basic understanding of the flow and mixing characteristics of pulsed turbulent opposed-jet streams. The sinusoidal pulsating flows with different temperature are imposed at opposed-jet inlets, which are mixed with each other in a confined flow channel. The current mathematical model taking the effect of temperature-dependent thermo-physical properties of fluid into account can present a good prediction for opposed-jet streams compared with experimental data. The numerical results indicate that introduction of temperature difference between opposed jet flows can lead to an asymmetric flow field immediately after jet impact, and the sinusoidal flow pulsations can effectively enhance mixing rate of opposed jets. Parameter studies are conducted for optimization of pulsed opposed jets. The effect of Reynolds number and flow pulsation as well as the configuration geometry on the mixing performance are discussed in detail. Examination of the flow and thermal field shows that the mixing rate is highly dependent on the vortex-induced mixing and residence time of jet fluid in the exit channel.     

Layered structure of turbulent plane wall jet

This paper investigates the layered structure of a turbulent plane wall jet at a distance from the nozzle exit. Based on the force balances in the mean momentum equation, the turbulent plane wall jet is divided into three regions: a boundary layer-like region (BLR) adjacent to the wall, a half free jet-like region (HJR) away from the wall, and a plug flow-like region (PFR) in between. In the PFR, the mean streamwise velocity is essentially the maximum velocity, and the simplified mean continuity and mean momentum equations result in a linear variation of the mean wall-normal velocity and Reynolds shear stress. In the HJR, as in a turbulent free jet, a proper scale for the mean wall-normal flow is the mean wall-normal velocity far from the wall and a proper scale for the Reynolds shear stress is the product of the maximum mean streamwise velocity and the velocity scale for the mean wall-normal flow. The BLR region can be divided into four sub-layers, similar to those in a canonical pressure-driven turbulent channel flow or shear-driven turbulent boundary layer flow. Building on the log-law for the mean streamwise velocity in the BLR, a new skin friction law is proposed for a turbulent wall jet. The new prediction agrees well with the correlation of Bradshaw and Gee (1960) over moderate Reynolds numbers, but gives larger skin frictions at higher Reynolds numbers.     

Numerical investigation of turbulent-jet primary breakup using one-dimensional turbulence

Primary breakup to form droplets at liquid surfaces is an important fundamental process to study as it determines the initial properties of the dispersed phase, which affect mixing rates, secondary breakup, droplet collisions, and flow separation within the dispersed flow region. Primary breakup can be regarded as one of the least developed model components for simulating and predicting liquid jet breakup. However, it is of paramount importance in many technical applications, e.g. fuel injection in engines and spray painting. This paper presents a numerical investigation of primary breakup of a turbulent liquid jet in still air at standard conditions using the one-dimensional turbulence (ODT) modeling framework. ODT is a stochastic model that simulates turbulent flow evolution along a notional 1D line of sight by applying instantaneous maps to represent the effect of individual turbulent eddies on property profiles. An important feature of ODT is the resolution of all relevant scales, both temporal and spatial. The restriction to one spatial dimension in ODT permits affordable high resolution of interfacial and single-phase property gradients, which is key to capturing the local behavior of the breakup process and allows simulations at high Reynolds and Weber numbers that are currently not accessible to direct numerical simulations (DNS).This paper summarizes our extensions of the ODT model to simulate geometrically simple jet breakup problems, including representations of Rayleigh wave breakup, turbulent breakup, and shear-driven breakup. Each jet breakup simulation consists of a short temporal channel section to initialize a turbulent velocity profile at the nozzle exit followed by an adjacent jet section. The simulations are carried out for jet exit Reynolds number of 11,500, 23,000, 46,000 and 92,000 while the Weber number is varied within the range 10²–10⁷. We present results on breakup statistics including spatial locations of droplet release, droplet sizes and liquid core length. The results on primary breakup are compared to experimental results and models.