Improved MMC-LES to compute the structure of a mixed-mode turbulent flame series

An improved version of the Lagrangian multiple mapping conditioning model coupled with large eddy simulation (MMC-LES) is applied to compute the structure of a range of turbulent mixed-mode flames stabilised on the Sydney piloted burner with compositionally inhomogeneous inlet. The calculations are compared with detailed measurements. While premixed MMC-LES models are yet to be developed, the most commonly deployed version and code for MMC-LES uses a mixture-fraction-based approach and a sparse-Lagrangian notional particle method that is most suitable for non-premixed combustion. The improvements reported here are due to (i) a more rigorous density coupling approach, (ii) a more refined mixing time scale, and (iii) selective particle intensification to improve localness in physical space as well as mixture fraction space. This intensification is selectively applied in upstream regions where mixed-mode combustion dominates and is relaxed further downstream as the jet flames transition to a non-premixed flame structure. The computed results show good agreement with experimental data hence confirming the feasibility of the improved MMC-LES approach in being able to account for mixed-mode combustion as well as for finite-rate chemistry effects. These promising results warrant future extension of the selective particle intensification method to become automatically adaptive.     

Filtered density function simulator on unstructured meshes

A new computational filtered density function (FDF) methodology is developed for large eddy simulation (LES) of turbulent reacting flows. This methodology is based on a Lagrangian Monte Carlo (MC) FDF solver constructed on a domain portrayed by an unstructured mesh. The base filtered transport equations on this mesh are solved by a finite-volume (FV) method. The consistency of the hybrid FV-MC solver and the realizability of the simulated results are demonstrated via LES of a temporally developing mixing layer. The overall performance of the model is appraised by comparison with direct numerical simulation (DNS) data. The algorithmic implementation in the commercial software ANSYS-FLUENT facilitates future FDF-LES of turbulent combustion in complex configurations.     

Finding thermoacoustic limit cycles for a ducted Burke-Schumann flame

This paper examines nonlinear thermoacoustic oscillations of a ducted Burke-Schumann diffusion flame. The nonlinear dynamics of the thermoacoustic system are studied using two distinct approaches. In the first approach, a continuation analysis is performed to find limit cycle amplitudes over a range of operating conditions. The strength of this approach is that one can characterize the coupled system’s nonlinear behaviour over a large parameter space with relative ease. It is not able to give physical insight into that behaviour, however. The second approach uses a Flame Describing Function (FDF) to characterize the flame’s response to harmonic velocity fluctuations over a range of forcing frequencies and forcing amplitudes, from which limit cycle amplitudes can be found. A strength of the FDF approach is that it reveals the physical mechanisms responsible for the behaviour observed. However, the calculation of the FDF is time consuming, and it must be recalculated if the flame’s operating conditions change. With the strengths and shortcomings of the two approaches in mind, this paper advocates combining the two to provide the dynamics over a large parameter space and, furthermore, physical insight into that behaviour at judiciously-chosen points in the parameter space. Further physical insight concerning the flame’s near-linear response at all forcing amplitudes is given by studying the forced flame in the time domain. It is shown that, for this flame model, the limit cycles arise because of the flame’s nonlinear behaviour when it is close to the inlet.     

Large eddy simulations of turbulent flames using the filtered density function model

Large eddy simulations (LES) of the Sandia/Sydney swirl burners (SM1 and SMA1) and the Sandia/Darmstadt piloted jet diffusion flame (Flame D) are performed. These flames are part of the database of turbulent reacting flows widely considered as benchmark test cases for validating turbulent-combustion models. In the simulations presented in this paper, the subgrid scale (SGS) closure model adopted for turbulence-chemistry interactions is based on the transport filtered density function (FDF) model. In the FDF model, the transport equation for the joint probability density function (PDF) of scalars is solved. The main advantage of this model is that the filtered reaction rates can be exactly computed. However, the density field, computed directly from the FDF solver and needed in the hydrodynamic equations, is noisy and causes numerical instability. Two numerical approaches that yield a smooth density field are examined. The two methods are based on transport equations for specific sensible enthalpy (h_s) and RT, where R is the gas constant and T is the temperature. Consistency of the two methods is assessed in a bluff-body configuration using Reynolds averaged Navier-Stokes (RANS) methodology in conjunction with the transported PDF method. It is observed that the h_s method is superior to the RT method. Both methods are used in LES of the SM1 burner. In the near-field region, the h_s method produces better predictions of temperature. However, in the far-field region, both methods show deviation from data. Simulations of the SMA1 burner and Flame D are also presented using the h_s method. Some deficiencies are seen in the predictions of the SMA1 burner that may be related to the simple chemical kinetics model and mixing model used in the simulations. Simulations of Flame D show good agreement with data. These results indicate that, while further improvements to the methodology are needed, the LES/FDF method has the potential to accurately predict complex turbulent flames.     

U-RANS／PDF方法计算钝体后漩涡脱落

本研究发展了U-RANS/PDF混合算法研究湍流和化学反应相互作用对燃烧稳定性的影响,采用有限体积和Monte Carlo相结合的方法在非结构网格中求解相容的U-RANS方程和脉动速度-湍流频率-标量的联合PDF方程.本文对钝体火焰驻定器后冷态流场进行了计算,结果表明此混合算法能够捕捉流场中非稳态的漩涡脱落现象.着重研究了湍流频率模型系数的改变对漩涡脱落频率以及拟序结构在动量输运中的作用的影响.     

Critical Exponents of Fully Directed Flights: A Laplacian Transformation Method

A Laplacian transformation method is presented for d-dimensional fully directed flights (FDFs) with continuous steplengths. The correlation exponents we obtained are in agreement with the simulation study for fully directed Levy flight and the analytical results for FDF on a simple lattice. We believe that all FDFs in isotropic space belong to the same universality class.     

Modeling viscoelastic flow with discrete methods

The hydrodynamics of viscoelastic materials (for example, polymer melts and solutions) presents interesting and complex phenomena, for example, instabilities and turbulent flow at very low Reynolds numbers due to normal stress effects and the existence of a finite stress relaxation time. This present work is motivated by renewed interest in instabilities in polymer flow. The majority of currently used numerical methods discretize a constitutive equation on a grid with finite difference or similar methods. We present work in progress in which we simulate viscoelastic flow with dissipative particle dynamics. The advantage of this approach is that many of the numerical instabilities of conventional methods can be avoided, and that the model gives clear physical insight into the origins of many viscoelastic flow instabilities.     

DSM-LPDF两相湍流模型及旋流两相流动的模拟

本文由流体－颗粒速度的拉氏联合概率密度函数（ＰＤＦ）输运方程出发,用Ｓｉｍｏｎｉｎ建议的Ｌａｎｇｅｖｉｎ模型封闭颗粒所遇到流体瞬时速度的条件期望项,并用Ｍｏｎｔｅ Ｃａｒｌｏ方法直接求解 ＰＤＦ输运方程,将其和求解流体雷诺应力方程模型的有限差分方法结合,建立了雷诺应力－拉氏ＰＤＦ（ＤＳＭ－ＬＰＤＦ,简称ＤＬ）两相湍流模型．用此模型模拟了旋流数为０．４７的突扩旋流气粒两相流动,并与文献中ＰＤＰＡ实验和用类似于单相流动湍流模型封闭方法的时平均统一二阶矩（ＵＳＭ）模型的预报进行了对比．     

A posteriori testing of the flame surface density transport equation for LES

Flame Surface Density (FSD) models for Large Eddy Simulation (LES) are implemented and tested for a canonical configuration and a practical bluff body stabilised burner, comparing common algebraic closures with a transport equation closure in the context of turbulent premixed combustion. The transported method is expected to yield advantages over algebraic closures, as the equilibrium of subgrid production and destruction of FSD is no longer enforced and resolved processes of strain, propagation and curvature are explicitly accounted for. These advantages might have the potential to improve the ability to capture large-scale unsteady flame propagation in situations with combustion instabilities or situations where the flame encounters progressive wrinkling with time. The initial study of a propagating turbulent flame in wind-tunnel turbulence shows that the Algebraic Flame Surface Density (FSDA) method can predict an excessively wrinkled flame under fine grid conditions, potentially increasing the consumption rate of reactants to artificially higher levels. In contrast, the Flame Surface Density Transport (FSDT) closure predicts a smooth flame front and avoids the formation of artificial flame cusps when the grid is refined. Five FSDA models and the FSDT approach are then applied to the LES of the Volvo Rig. The predicted mean velocities are found to be relatively insensitive to the use of the FSDT and FSDA approaches, whereas temperature predictions exhibit appreciable differences for different formulations. The FSDT approach yields very similar temperature predictions to two of the tested FSDA models, quantitatively capturing the mean temperature. Grid refinement is found to improve the FSDT predictions of the mean flame spread. Overall, the paper demonstrates that the apparently complicated FSD transport equation approach can be implemented and applied to realistic, strongly wrinkled flames with good success, and opens up the field for further work to improve the models and the overall FSDT approach.     

Effects of combined dimension reduction and tabulation on the simulations of a turbulent premixed flame using a large-eddy simulation/probability density function method

A turbulent lean-premixed propane–air flame stabilised by a triangular cylinder as a flame-holder is simulated to assess the accuracy and computational efficiency of combined dimension reduction and tabulation of chemistry. The computational condition matches the Volvo rig experiments. For the reactive simulation, the Lagrangian Large-Eddy Simulation/Probability Density Function (LES/PDF) formulation is used. A novel two-way coupling approach between LES and PDF is applied to obtain resolved density to reduce its statistical fluctuations. Composition mixing is evaluated by the modified Interaction-by-Exchange with the Mean (IEM) model. A baseline case uses In Situ Adaptive Tabulation (ISAT) to calculate chemical reactions efficiently. Its results demonstrate good agreement with the experimental measurements in turbulence statistics, temperature, and minor species mass fractions. For dimension reduction, 11 and 16 represented species are chosen and a variant of Rate Controlled Constrained Equilibrium (RCCE) is applied in conjunction with ISAT to each case. All the quantities in the comparison are indistinguishable from the baseline results using ISAT only. The combined use of RCCE/ISAT reduces the computational time for chemical reaction by more than 50%. However, for the current turbulent premixed flame, chemical reaction takes only a minor portion of the overall computational cost, in contrast to non-premixed flame simulations using LES/PDF, presumably due to the restricted manifold of purely premixed flame in the composition space. Instead, composition mixing is the major contributor to cost reduction since the mean-drift term, which is computationally expensive, is computed for the reduced representation. Overall, a reduction of more than 15% in the computational cost is obtained.     

Multiscale three-point velocity increment correlation in turbulent flows

The three-point velocity increment correlation function is proposed to represent the multiscale correlations in turbulent flows. The inertial–inertial correlation and the inertial–dissipative correlation are discussed due to their endogenetic properties in turbulence and their roles in large-eddy simulation. The zero-correlation points are then emphasized as equilibrium points between them. The credibility of this theoretical result is numerically verified in both isotropic and anisotropic flows. Results imply the universality of this zero-correlation scaling in different turbulent flows. This work is expected to be a dependable theoretical base for creating multiscale subgrid models in large-eddy simulation.     

Modelling of turbulent lifted jet flames using flamelets: a priori assessment and a posteriori validation

This study focuses on the modelling of turbulent lifted jet flames using flamelets and a presumed Probability Density Function (PDF) approach with interest in both flame lift-off height and flame brush structure. First, flamelet models used to capture contributions from premixed and non-premixed modes of the partially premixed combustion in the lifted jet flame are assessed using a Direct Numerical Simulation (DNS) data for a turbulent lifted hydrogen jet flame. The joint PDFs of mixture fraction Z and progress variable c, including their statistical correlation, are obtained using a copula method, which is also validated using the DNS data. The statistically independent PDFs are found to be generally inadequate to represent the joint PDFs from the DNS data. The effects of Z–c correlation and the contribution from the non-premixed combustion mode on the flame lift-off height are studied systematically by including one effect at a time in the simulations used for a posteriori validation. A simple model including the effects of chemical kinetics and scalar dissipation rate is suggested and used for non-premixed combustion contributions. The results clearly show that both Z–c correlation and non-premixed combustion effects are required in the premixed flamelets approach to get good agreement with the measured flame lift-off heights as a function of jet velocity. The flame brush structure reported in earlier experimental studies is also captured reasonably well for various axial positions. It seems that flame stabilisation is influenced by both premixed and non-premixed combustion modes, and their mutual influences.     

Conditional dissipation of scalars in homogeneous turbulence: Closure for MMC modelling

While the mean and unconditional variance are to be predicted well by any reasonable turbulent combustion model, these are generally not sufficient for the accurate modelling of complex phenomena such as extinction/reignition. An additional criterion has been recently introduced: accurate modelling of the dissipation timescales associated with fluctuations of scalars about their conditional mean (conditional dissipation timescales). Analysis of Direct Numerical Simulation (DNS) results for a passive scalar shows that the conditional dissipation timescale is of the order of the integral timescale and smaller than the unconditional dissipation timescale. A model is proposed: the conditional dissipation timescale is proportional to the integral timescale. This model is used in Multiple Mapping Conditioning (MMC) modelling for a passive scalar case and a reactive scalar case, comparing to DNS results for both. The results show that this model improves the accuracy of MMC predictions so as to match the DNS results more closely using a relatively-coarse spatial resolution compared to other turbulent combustion models.     

A consistent LES/filtered-density function formulation for the simulation of turbulent flames with detailed chemistry

A hybrid large-Eddy simulation/filtered-density function (LES–FDF) methodology is formulated for simulating variable density turbulent reactive flows. An indirect feedback mechanism coupled with a consistency measure based on redundant density fields contained in the different solvers is used to construct a robust algorithm. Using this novel scheme, a partially premixed methane/air flame is simulated. To describe transport in composition space, a 16-species reduced chemistry mechanism is used along with the interaction-by-exchange with the mean (IEM) model. For the micro-mixing model, typically a constant ratio of scalar to mechanical time-scale is assumed. This parameter can have substantial variations and can strongly influence the combustion process. Here, a dynamic time-scale model is used to prescribe the mixing time-scale, which eliminates the time-scale ratio as a model constant. Two different flame configurations, namely, Sandia flames D and E are studied. Comparison of simulated radial profiles with experimental data show good agreement for both flames. The LES–FDF simulations accurately predict the increased extinction near the inlet and re-ignition further downstream. The conditional mean profiles show good agreement with experimental data for both flames.     

Multiple mapping conditioning of turbulent jet diffusion flames

The multiple mapping conditioning (MMC) approach is applied to two non-piloted CH₄/H₂/N₂ turbulent jet diffusion flames with Reynolds numbers of Re = 15,200 and 22,800. The work presented here examines primarily the suitability of MMC to simulate CH₄/H₂ flames with varying Re numbers. The equations are solved in a prescribed Gaussian reference space with only one stochastic reference variable emulating the fluctuations of mixture fraction. The mixture fraction is considered as the only major species on which the remaining minor species are conditioned. Fluctuations around the conditional means are ignored. It is shown that the statistics of the mapped reference field are an accurate model for the statistics of the physical field for both flames. A transformation of the Gaussian reference space introduced in previous work on MMC is used to express the MMC model in the same form as CMC. The most important advantage of this transformation is that the conditionally averaged scalar dissipation term is in a closed form. The corresponding temperature and reactive species predictions are generally in good agreement with experimental data. The application to real laboratory flames and the assessment of the new conditional scalar dissipation model for the closure of the singly conditioned CMC equation is the major novelty of this paper. The results are therefore primarily examined with respect to changes of the conditionally averaged quantities in mixture fraction space.     

Study of mass consistency LES/FDF techniques for chemically reacting flows

A hybrid large eddy simulation/filtered density function (LES/FDF) approach is used for studying chemically reacting flows with detailed chemistry. In particular, techniques utilised for ensuring a mass consistent coupling between LES and FDF are discussed. The purpose of these techniques is to maintain a correct spatial distribution of the computational particles representing specified amounts of fluid. A particular mass consistency technique due to Y.Z. Zhang and D.C. Haworth (A general mass consistency algorithm for hybrid particle/finite-volume PDF methods, J. Comput. Phys. 194 (2004), pp. 156–193) and their associated algorithms are implemented in a pressure-based computational fluid dynamics code suitable for the simulation of variable density flows, representative of those encountered in actual combustion applications. To assess the effectiveness of the referenced technique for enforcing LES/FDF mass consistency, two- and three-dimensional simulations of a temporal mixing layer using detailed and reduced chemistry mechanisms are carried out. The parametric analysis performed focuses on determining the influence on the level of mass consistency errors of parameters such as the initial number of particles per cell and the initial density ratio of the mixing layers. Particular emphasis is put on the computational burden that represents the use of such a mass consistency technique. The results show the suitability of this type of technique for ensuring the mass consistency required when utilising hybrid LES/FDF approaches. The level of agreement of the computed results with experimental data is also illustrated.     

Stratified turbulent Bunsen flames: flame surface analysis and flame surface density modelling

In this paper it is investigated whether the Flame Surface Density (FSD) model, developed for turbulent premixed combustion, is also applicable to stratified flames. Direct Numerical Simulations (DNS) of turbulent stratified Bunsen flames have been carried out, using the Flamelet Generated Manifold (FGM) reduction method for reaction kinetics. Before examining the suitability of the FSD model, flame surfaces are characterized in terms of thickness, curvature and stratification.

All flames are in the Thin Reaction Zones regime, and the maximum equivalence ratio range covers 0.1?φ?1.3. For all flames, local flame thicknesses correspond very well to those observed in stretchless, steady premixed flamelets. Extracted curvature radii and mixing length scales are significantly larger than the flame thickness, implying that the stratified flames all burn in a premixed mode. The remaining challenge is accounting for the large variation in (subfilter) mass burning rate.

In this contribution, the FSD model is proven to be applicable for Large Eddy Simulations (LES) of stratified flames for the equivalence ratio range 0.1?φ?1.3. Subfilter mass burning rate variations are taken into account by a subfilter Probability Density Function (PDF) for the mixture fraction, on which the mass burning rate directly depends. A priori analysis point out that for small stratifications (0.4?φ?1.0), the replacement of the subfilter PDF (obtained from DNS data) by the corresponding Dirac function is appropriate. Integration of the Dirac function with the mass burning rate m=m(φ), can then adequately model the filtered mass burning rate obtained from filtered DNS data. For a larger stratification (0.1?φ?1.3), and filter widths up to ten flame thicknesses, a β-function for the subfilter PDF yields substantially better predictions than a Dirac function. Finally, inclusion of a simple algebraic model for the FSD resulted only in small additional deviations from DNS data, thereby rendering this approach promising for application in LES.     

A PDF projection method: A pressure algorithm for stand-alone transported PDFs

In this paper, a new formulation of the projection approach is introduced for stand-alone probability density function (PDF) methods. The method is suitable for applications in low-Mach number transient turbulent reacting flows. The method is based on a fractional step method in which first the advection–diffusion–reaction equations are modelled and solved within a particle-based PDF method to predict an intermediate velocity field. Then the mean velocity field is projected onto a space where the continuity for the mean velocity is satisfied. In this approach, a Poisson equation is solved on the Eulerian grid to obtain the mean pressure field. Then the mean pressure is interpolated at the location of each stochastic Lagrangian particle. The formulation of the Poisson equation avoids the time derivatives of the density (due to convection) as well as second-order spatial derivatives. This in turn eliminates the major sources of instability in the presence of stochastic noise that are inherent in particle-based PDF methods. The convergence of the algorithm (in the non-turbulent case) is investigated first by the method of manufactured solutions. Then the algorithm is applied to a one-dimensional turbulent premixed flame in order to assess the accuracy and convergence of the method in the case of turbulent combustion. As a part of this work, we also apply the algorithm to a more realistic flow, namely a transient turbulent reacting jet, in order to assess the performance of the method.     

On mean flow universality of turbulent wall flows. I. High Reynolds number flow analysis

ABSTRACT

The universality and mathematical physical structure of wall-bounded turbulent flows is a topic of discussions over many decades. There is no agreement about questions like what is the physical mean flow structure, how universal is it, and how universal are theoretical concepts for local and global flow variations. These questions are addressed by using latest direct numerical simulation (DNS) data at moderate Reynolds numbers Re and experimental data up to extreme Re. The mean flow structure is explained by analytical models for three canonical wall-bounded turbulent flows (channel flow, pipe flow, and the zero-pressure gradient turbulent boundary layer). Thorough comparisons with DNS and experimental data provide support for the validity of models. Criteria for veritable physics derived from observations are suggested. It is shown that the models presented satisfy these criteria. A probabilistic interpretation of the mean flow structure shows that the physical constraints of equal entropies and equally likely mean velocity values in a region unaffected by boundary effects impose a universal log-law structure. The structure of wall-bounded turbulent flows is much more universal than previously expected. There is no discrepancy between local logarithmic velocity variations and global friction law and bulk velocity variations. Flow effects are limited to the minimum: the difference of having a bounded or unbounded domain, and the variation range of mean velocity values allowed by the geometry.     

An abstraction layer for efficient memory management of tabulated chemistry and flamelet solutions

A large number of methods for simulating reactive flows exist, some of them, for example, directly use detailed chemical kinetics or use precomputed and tabulated flame solutions. Both approaches couple the research fields computational fluid dynamics and chemistry tightly together using either an online or offline approach to solve the chemistry domain. The offline approach usually involves a method of generating databases or so-called Lookup-Tables (LUTs). As these LUTs are extended to not only contain material properties but interactions between chemistry and turbulent flow, the number of parameters and thus dimensions increases. Given a reasonable discretisation, file sizes can increase drastically. The main goal of this work is to provide methods that handle large database files efficiently. A Memory Abstraction Layer (MAL) has been developed that handles requested LUT entries efficiently by splitting the database file into several smaller blocks. It keeps the total memory usage at a minimum using thin allocation methods and compression to minimise filesystem operations. The MAL has been evaluated using three different test cases. The first rather generic one is a sequential reading operation on an LUT to evaluate the runtime behaviour as well as the memory consumption of the MAL. The second test case is a simulation of a non-premixed turbulent flame, the so-called HM1 flame, which is a well-known test case in the turbulent combustion community. The third test case is a simulation of a non-premixed laminar flame as described by McEnally in 1996 and Bennett in 2000. Using the previously developed solver ‘flameletFoam’ in conjunction with the MAL, memory consumption and the performance penalty introduced were studied. The total memory used while running a parallel simulation was reduced significantly while the CPU time overhead associated with the MAL remained low.