贴体坐标系下模型加力室的大涡模拟

本文利用贴体网格对带V形槽稳定器模型加力燃烧室素流化学反应流流动进行大涡模拟的研究。采用区域法生成模型加力燃烧室的二维贴体网格,并采用多区域耦合法进行区域之间的数据传递,求解加力室整体流场。采用k方程亚网格尺度模型和亚网格EBU燃烧模型分别估算其亚网格紊流粘性和化学反应速率,用热通量辐射模型估算辐射通量,并用交错网格下SIMPLE算法和混合差分格式求解离散方程,壁面函数处理固壁边界条件。计算结果显示了稳定器后面的回流区气流结构,所得的热态流场模拟结果与实验比较吻合,表明采用贴体网格对模型加力燃烧室进行大涡模拟能真实反映流体流动及燃烧过程。     

Effects of buoyancy on turbulent premixed V-flames by large-eddy simulation

Buoyancy effects on turbulent premixed V-flames are investigated under normal gravity (+g) and reversed gravity (–g). Numerical simulations employ large eddy simulation (LES) with a dynamic model for sub-grid scale stress. With the assumption of fast chemistry combustion, a progress variable c-equation is applied to describe the flame front propagation. The equations are solved using a projection-based fractional step method in two dimensions for low-Mach number flows. Computed LES results of buoyancy effects on flame angle and flame brush thickness are consistent with those obtained from experiments. In both +g and –g conditions, the effects of buoyancy become important with increase in Richardson number (Ri). Buoyancy force tends to close up the flame under +g, but has the opposite effect under –g. Buoyancy force also suppresses flame wrinkling in +g and enhances wrinkling in –g. While there is a lack of experimental data available, computed axial velocity is shown to be significantly affected by buoyancy downstream from the flame holder under moderate Reynolds number.     

Large Eddy Simulations of a piloted lean premix jet flame using finite-rate chemistry

A Large Eddy Simulation (LES) model capable of accurately representing finite-rate chemistry effects in turbulent premixed combustion is presented. The LES computations use finite-rate chemistry and implicit LES combustion modelling to simulate an experimentally well-documented lean-premixed jet flame stabilized by a stoichiometric pilot. The validity of the implicit LES assumption is discussed and criteria are expressed in terms of subgrid scale Damköhler and Karlovitz numbers. Simulation results are compared to experimental data for velocity, temperature and species mass fractions of CH₄, CO and OH. The simulation results highlight the validity and capability of the present approach for the flame and in general the combustion regime examined. A sensitivity analysis to the choice of the finite-rate chemistry mechanism is reported, this analysis indicates that the one and two-step global reaction mechanisms evaluated fail to capture the reaction layer with sufficient accuracy, while a 20-species skeletal mechanism reproduces the experimental observations accurately including the key finite-rate chemistry indicators CO and OH. The LES results are shown to be grid insensitive and that the grid resolution within the bounds examined is far less important compared to the sensitivity of the finite-rate chemistry representation. The results are analyzed in terms of the flame dynamics and it is shown that intense small scale mixing (high Karlovitz number) between the pilot and the jet is an important mechanism for the stabilization of the flame.     

Identification and analysis of instability in non-premixed swirling flames using LES

Large eddy simulations (LES) of turbulent non-premixed swirling flames based on the Sydney swirl burner experiments under different flame characteristics are used to uncover the underlying instability modes responsible for the centre jet precession and large scale recirculation zone. The selected flame series known as SMH flames have a fuel mixture of methane-hydrogen (50:50 by volume). The LES solves the governing equations on a structured Cartesian grid using a finite volume method, with turbulence and combustion modelling based on the localised dynamic Smagorinsky model and the steady laminar flamelet model respectively. The LES results are validated against experimental measurements and overall the LES yields good qualitative and quantitative agreement with the experimental observations. Analysis showed that the LES predicted two types of instability modes near fuel jet region and bluff body stabilised recirculation zone region. The mode I instability defined as cyclic precession of a centre jet is identified using the time periodicity of the centre jet in flames SMH1 and SMH2 and the mode II instability defined as cyclic expansion and collapse of the recirculation zone is identified using the time periodicity of the recirculation zone in flame SMH3. Finally frequency spectra obtained from the LES are found to be in good agreement with the experimentally observed precession frequencies.     

Large eddy simulation of a pulverized coal jet flame ignited by a preheated gas flow

Large eddy simulation (LES) is applied to a pulverized coal jet flame ignited by a preheated gas flow. The simulation results are compared to experimental data obtained for the inlet stoichiometric ratios of 0.14, 0.22, and 0.36. An accurate and computationally inexpensive devolatilization model suitable for combustion simulation in LES is proposed and incorporated into the LES. The numerical results of gas temperature and coal burnout on the centerline show good agreement with the experimental data. Two kinds of lift-off heights are introduced to verify the combustion simulation. One is the height from the primary nozzle exit to the starting point of the growing flame region. The other is the height from the primary nozzle exit to the starting point of the continuous flame region. The calculated results of the two lift-off heights show good agreement with the experimental data. In contrast to LES, the standard k–ε model overestimates the lift-off heights because it calculates time-averaged temperature which does not contain information about local flame structure. The stoichiometric ratio in the gas phase at the starting point of the growing flame region is found to be independent of the inlet stoichiometric ratio in the range from 0.14 to 0.36.     

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.     

LES-CMC simulations of a turbulent bluff-body flame

The large Eddy simulations (LES)-conditional moment closure (CMC) method with detailed chemistry is applied to a bluff-body stabilized flame. Computations of the velocity and mixture fraction fields show good agreement with the experiments. Temperature and major species are well-predicted throughout the flame with the exception of the flow regions in the outer shear layer close to the nozzle where the pure mixing between hot recirculating products and fresh oxidizer cannot be captured. LES-CMC generally improves on results obtained with RANS-CMC and on LES that uses one representative flamelet to model the dependence of reactive species on mixture fraction. Simulated CO mass fractions are generally in good agreement with the experimental data although a 10% overprediction can be found at downstream positions. NO predictions show a distinct improvement over the flamelet approach, however, simulations overpredict NO mass fractions at all downstream locations due to an overprediction of temperature close to the nozzle. The potential of LES-CMC to predict unsteady finite rate effects is demonstrated by the prediction of endothermic—or “flame cooling”—regions close to the neck of the recirculation zone that favours ethylene production via the methane fuel decomposition channel.     

Experimental and numerical study of a conical turbulent partially premixed flame

The structure and dynamics of a turbulent partially premixed methane/air flame in a conical burner were investigated using laser diagnostics and large-eddy simulations (LES). The flame structure inside the cone was characterized in detail using LES based on a two-scalar flamelet model, with the mixture fraction for the mixing field and level-set G-function for the partially premixed flame front propagation. In addition, planar laser induced florescence (PLIF) of CH and chemiluminescence imaging with high speed video were performed through a glass cone. CH and CH₂O PLIF were also used to examine the flame structures above the cone. It is shown that in the entire flame the CH layer remains very thin, whereas the CH₂O layer is rather thick. The flame is stabilized inside the cone a short distance above the nozzle. The stabilization of the flame can be simulated by the triple-flame model but not the flamelet-quenching model. The results show that flame stabilization in the cone is a result of premixed flame front propagation and flow reversal near the wall of the cone which is deemed to be dependent on the cone angle. Flamelet based LES is shown to capture the measured CH structures whereas the predicted CH₂O structure is somewhat thinner than the experiments.     

Large-eddy simulations of a stratified-charge direct-injection spark-ignition engine: Comparison with experiment and analysis of cycle-to-cycle variations

Multiple-cycle large-eddy simulations (LES) have been performed for an optically accessible, single-cylinder, four-stroke-cycle, spray-guided direct-injection spark-ignition (SG-DISI) engine operating in a stratified globally fuel-lean mode. The simulations combine a standard Smagorinsky turbulence model, a stochastic Lagrangian parcel method for liquid fuel injection and fuel spray modeling, a simple energy-deposition spark-ignition model, and a modified thickened flame model for turbulent flame propagation through highly stratified reactant mixtures. Comparisons between simulations and experiments include individual-cycle and ensemble-average pressure and apparent-heat-release-rate traces, individual-cycle and ensemble-average indicated mean effective pressures (IMEP), and instantaneous two-dimensional vapor-equivalence-ratio contours. Although the number of LES cycles is small (35), the results show that the simulations are able to capture the global combustion behavior that is observed in the experiments, including cycle-to-cycle variations. The simulation results are then analyzed further to provide insight into the conditions that lead to misfire versus robust combustion. As has been reported in earlier experimental and LES studies for homogeneous-charge SI engines, local conditions in the vicinity of the spark gap at the time of ignition largely determine the subsequent flame development. However, in contrast to homogeneous-charge engines, no single local or global quantity correlates as strongly with the eventual peak pressure or IMEP for each cycle. Rather, it is the interplay among the early flame kernel, the velocity field that it experiences, and the fuel distribution that it encounters that ultimately determines the fate of each combustion event. Deeper analysis and quantitative statistical comparisons between experiments and simulations will require the simulation of larger numbers of engine cycles.     

Computationally-efficient and accurate particle PDF simulations of turbulent combustion using coupled pre-partitioned adaptive chemistry and tabulation

LES/PDF methods are known to provide accurate results for challenging turbulent combustion configurations with strong turbulence-chemistry interactions. These methods are generally applicable as they do not make any assumptions on the topology of the underlying flame structure. However, this added generality comes at an increased computational cost. To mitigate this added cost, the majority of the LES/PDF computations performed to date utilize reduced mechanisms. We recently presented a coupled pre-partitioned adaptive chemistry (PPAC) and tabulation (ISAT) methodology (Newale et al., Comb. Th. Mod., 2019), which retains the fidelity of the detailed mechanism, while keeping the computational cost affordable. This methodology was tested in a partially-stirred reactor configuration. In this work, we describe the developments required for a holistic integration of PPAC-ISAT with a LES/PDF framework. We examine the performance of this coupled methodology in two LES/PDF configurations of Sandia flame D. A smaller simulation domain is initially utilized to characterize the efficiency and accuracy of standalone PPAC and coupled PPAC-ISAT in detail. Then, the performance of PPAC-ISAT is examined in a full-scale LES/PDF simulation. We show that the coupled PPAC-ISAT LES/PDF captures the resolved mean and RMS profiles of temperature and major species mass fractions to within 2% and OH to within 5%, with a reduction in the average simulation wall clock time per time step of 39% over an ISAT implementation using the detailed mechanism.     

DME/Oxygen wall-stabilized premixed cool flame

Wall-stabilized cool flames have been studied through numerical analysis and a series of experiments. One- and two-dimensional numerical simulations were performed to estimate the characteristics of the wall-stabilized cool flames, such as flammability and temperature/species distributions. Based on the computational results, the ignition condition of the cool flame at a fixed wall temperature has been identified with the strain rate between the van't Hoff point and the cool flame extinction point. In experiments with an impinging jet burner and a heated plate, the spontaneous ignition of the cool flame on the heated wall has been successfully established under the conditions predicted by the present numerical simulations. Spatial distributions of the HCHO concentration and flame temperature were measured through formaldehyde Planar Laser-Induced Fluorescence (HCHO-PLIF) and thermocouple measurements, respectively. It is found that the measurement data show a reasonable accordance with the simulation results with reduced low-temperature reactivity.     

Flamelet structure in turbulent premixed swirling oxy-combustion of methane

Two key flame macrostructures in swirling flows have been observed in experiments of oxy-combustion (as well as air-combustion); as the equivalence ratio is raised, the flame moves from being stabilized on just the inner shear layer (Flame III) to getting stabilized on both the inner and outer shear layers (Flame IV). We report results of an LES investigation of two different inlet oxy-fuel mixtures, in a turbulent swirling flow at $\text{Re}=20,000,$ that capture these two macrostructures. Previous work on the effects of heat loss have mostly focused on its impact on macro-scale observations. In this paper, we examine how heat loss impacts the flame microstructures as well for these two macrostructures. For both flames, the flamelet structure, as represented by a scatter plot of the normalized fuel concentration against the normalized temperature, depends on whether the combustor walls are adiabatic or non-adiabatic. For the adiabatic case, the flamelets of both macrostructures behave like strained flames. When wall heat transfer is included, Flame III microstructure is more bimodal. Since this flame extends farther downstream and part of it propagates along the walls, heat transfer has a greater impact on it’s microstructure. These results show that heat loss impacts not just the macro properties of the flame such as its shape or interactions with the wall, but also fundamentally changes its internal structure. Scatter plots of the turbulent flames are constructed and compared to different 1D laminar flame profiles (e.g., strained or with heat loss), and comparisons suggest the important role of the wall thermal boundary conditions in the accurate simulations of combustion dynamics and interpretations of experimental data, including data reduction and scaling.     

Predictions of Conjugate Heat Transfer in Turbulent Channel Flow Using Advanced Wall-Modeled Large Eddy Simulation Techniques

In this paper, advanced wall-modeled large eddy simulation (LES) techniques are used to predict conjugate heat transfer processes in turbulent channel flow. Thereby, the thermal energy transfer process involves an interaction of conduction within a solid body and convection from the solid surface by fluid motion. The approaches comprise a two-layer RANS–LES approach (zonal LES), a hybrid RANS–LES representative, the so-called improved delayed detached eddy simulation method (IDDES) and a non-equilibrium wall function model (WFLES), respectively. The results obtained are evaluated in comparison with direct numerical simulation (DNS) data and wall-resolved LES including thermal cases of large Reynolds numbers where DNS data are not available in the literature. It turns out that zonal LES, IDDES and WFLES are able to predict heat and fluid flow statistics along with wall shear stresses and Nusselt numbers accurately and that are physically consistent. Furthermore, it is found that IDDES, WFLES and zonal LES exhibit significantly lower computational costs than wall-resolved LES. Since IDDES and especially zonal LES require considerable extra work to generate numerical grids, this study indicates in particular that WFLES offers a promising near-wall modeling strategy for LES of conjugated heat transfer problems. Finally, an entropy generation analysis using the various models showed that the viscous entropy production is zero inside the solid region, peaks at the solid–fluid interface and decreases rapidly with increasing wall distance within the fluid region. Except inside the solid region, where steep temperature gradients lead to high (thermal) entropy generation rates, a similar behavior is monitored for the entropy generation by heat transfer process.     

Large Eddy Simulation of a turbulent lifted flame using multi-modal manifold-based models: Feasibility and interpretability

A general model for multi-modal turbulent combustion is achievable with two-dimensional manifold equations that use the mixture fraction and a generalized progress variable as coordinates. Information about the underlying mode of combustion is encoded in three scalar dissipation rates that appear as parameters in the two-dimensional equations. In this work, Large Eddy Simulation (LES) of a multi-modal turbulent lifted hydrogen jet flame in a vitiated coflow is performed using this new turbulent combustion model, leveraging both convolution-on-the-fly and In-Situ Adaptive Tabulation for computational tractability. The simulation predicts a lifted flame consistent with observations from past experiments. The feasibility of such a model implemented in LES is examined, and the cost per timestep is found to be comparable to conventional one-dimensional manifold-based models describing one asymptotic mode of combustion. Additionally, the model provides clear interpretability, allowing for combustion mode analysis to be performed with ease by evaluating the scalar dissipation rates and generalized progress variable source term. This analysis is used to show that the flame is stabilized by autoignition and has a trailing nonpremixed flame. Furthermore, transport of progress variable from the most reactive mixture fraction towards richer mixtures at the centerline is found to be important.     

Partially premixed flamelet in LES of acetone spray flames

An effective partially premixed flamelet model for large eddy simulation (LES) of turbulent spray combustion is formulated. Different flame regimes are identified with a flame index defined by budget terms in a 2-D multi-phase flamelet formulation, and the application in LES of partially pre-vaporized spray flames shows a favorable agreement with experiments. Simulations demonstrate that, compared to the conventional single-regime flamelets, the present partially premixed flamelet formulation shows its ability in capturing the subgrid regime transitions, yielding a well prediction of peak gas temperature and the downstream flame spreading. A propagating premixed flame front is found coupled with a trailing diffusion burning through the spray evaporation, and the spray effect on regime discrimination is manifested with transport budget analysis. A two-phase regime indicator is then proposed, by which the evaporation-dictated regime is properly described. Its intended use will rely on both gas and spray flamelet structures.     

Computational Modeling of Boundary Layer Flashback in a Swirling Stratified Flame Using a LES-Based Non-Adiabatic Tabulated Chemistry Approach

When operating under lean fuel–air conditions, flame flashback is an operational safety issue in stationary gas turbines. In particular, with the increased use of hydrogen, the propagation of the flame through the boundary layers into the mixing section becomes feasible. Typically, these mixing regions are not designed to hold a high-temperature flame and can lead to catastrophic failure of the gas turbine. Flame flashback along the boundary layers is a competition between chemical reactions in a turbulent flow, where fuel and air are incompletely mixed, and heat loss to the wall that promotes flame quenching. The focus of this work is to develop a comprehensive simulation approach to model boundary layer flashback, accounting for fuel–air stratification and wall heat loss. A large eddy simulation (LES) based framework is used, along with a tabulation-based combustion model. Different approaches to tabulation and the effect of wall heat loss are studied. An experimental flashback configuration is used to understand the predictive accuracy of the models. It is shown that diffusion-flame-based tabulation methods are better suited due to the flashback occurring in relatively low-strain and lean fuel–air mixtures. Further, the flashback is promoted by the formation of features such as flame tongues, which induce negative velocity separated boundary layer flow that promotes upstream flame motion. The wall heat loss alters the strength of these separated flows, which in turn affects the flashback propensity. Comparisons with experimental data for both non-reacting cases that quantify fuel–air mixing and reacting flashback cases are used to demonstrate predictive accuracy.     

Large eddy simulation of piloted pulverised coal combustion using extended flamelet/progress variable model

An extended flamelet/progress variable (EFPV) model for simulating pulverised coal combustion (PCC) in the context of large eddy simulation (LES) is proposed, in which devolatilisation, char surface reaction and radiation are all taken into account. The pulverised coal particles are tracked in the Lagrangian framework with various sub-models and the sub-grid scale (SGS) effects of turbulent velocity and scalar fluctuations on the coal particles are modelled by the velocity-scalar joint filtered density function (VSJFDF) model. The presented model is then evaluated by LES of an experimental piloted coal jet flame and comparing the numerical results with the experimental data and the results from the eddy break up (EBU) model. Detailed quantitative comparisons are carried out. It is found that the proposed model performs much better than the EBU model on radial velocity and species concentrations predictions. Comparing against the adiabatic counterpart, we find that the predicted temperature is evidently lowered and agrees well with the experimental data if the conditional sampling method is adopted.     

旋流和无旋突扩流动的LES和RANS模拟

本文用smagorinsky-Lilly亚网格尺度湍流模型对旋流突扩流动(s=0.53)和无旋突扩流动(s=0)进行了大涡模拟(LES模拟),同时分别用压力应变项为IPCM和IPCM+Wall模型的雷诺应力方程模型进行了RANS模拟,和LES的统计结果对比。LES的统计结果与雷诺应力模型的模拟结果及实验对照表明,LES结果与实验结果的吻合比雷诺应力模型的好,说明所用的亚网格尺度湍流模型对旋流流动是适用的,LES结果是可信的。LES的瞬态结果揭示出在旋流作用下,流场中存在复杂的旋涡脱落现象。大涡结构极易破碎成小涡,而在无旋突扩流动的情况下,由于剪切的作用更强,大涡结构的尺寸和范围比旋流流动的要大。     

Large-eddy simulation of MILD combustion using partially stirred reactor approach

Subgrid-scale (SGS) parameterization and method for calculating filtered reaction rate are critical components of an accurate large-eddy simulation (LES) of turbulent flames. In this study, we integrate gradient-type structural SGS models with a partially stirred reactor approach by using detailed chemical kinetics to simulate a turbulent methane/hydrogen jet flame under moderate or intense low-oxygen dilution (MILD) conditions. The study examines two oxygen dilution levels. The framework is assessed through a systematic and comprehensive comparison of temperature, and mass fractions of major and minor species with experimental data and other reference simulation results. Overall, the statistics of the combustion field show excellent agreement with measurements at different axial locations, and a significant improvement compared to some previous simulations. It suggests that the proposed nonlinear LES framework is able to accurately model MILD combustion with reasonable computational cost.     

An assessment of large eddy simulations of premixed flames propagating past repeated obstacles

This paper presents an assessment of Large Eddy Simulations (LES) in calculating the structure of turbulent premixed flames propagating past solid obstacles. One objective of the present study is to evaluate the LES simulations and identify the drawbacks in accounting the chemical reaction rate. Another objective is to analyse the flame structure and to calculate flame speed, generated overpressure at different time intervals following ignition of a stoichiometric propane/air mixture. The combustion chamber has built-in repeated solid obstructions to enhance the turbulence level and hence increase the flame propagating speed. Various numerical tests have also been carried out to determine the regimes of combustion at different stages of the flame propagation. These have been identified from the calculated results for the flow and flame characteristic parameters. It is found that the flame lies within the ‘thin reaction zone’ regime which supports the use of the laminar flamelet approach for modelling turbulent premixed flames. A submodel to calculate the model coefficient in the algebraic flame surface density model is implemented and examined. It is found that the LES predictions are slightly improved owing to the calculation of model coefficient by using submodel. Results are presented and discussed in this paper are for the flame structure, position, speed, generated pressure and the regimes of combustion during all stages of flame propagation from ignition to venting. The calculated results are validated against available experimental data.