Large Eddy Simulation of a Novel Gas-Assisted Coal Combustion Chamber

In this work a recently presented combustion chamber that is specifically designed for the investigation of gas-assisted coal combustion and the validation of models is simulated under reactive conditions for the first time. In the configuration coal combustion is assisted and stabilized by a methane flame. In the course of the investigation, the configuration’s complexity is increased successively. Results of the isothermal single-phase flow are discussed first. Subsequently, reproducibility of the single-phase methane flame by means of the applied modeling approach is evaluated. In a further step, coal particles having the same thermal power as the methane flame are injected into the configuration. Particle histories, the conversion of the coal particles as well as its retroactive effect on the gas phase are investigated. Experimental results based on laser diagnostics are provided for all operating points and used for comparison with numerical results. Gas phase velocity fields for all operating points are available. In order to identify the reaction in the reactive single-phase case planar laser induced fluorescence of the OH-radical (OH-PLIF) was applied. Overall good agreement between numerical and experimental results could be obtained. In the Large Eddy Simulation (LES) a Flamelet Generated Manifold (FGM) based model is utilized. The four-dimensional manifold is spanned by two mixture fractions, a reaction progress variable and the enthalpy on which the gas phase chemistry gets mapped onto. Thereby, the model accounts for both, volatiles reaction and char conversion. Furthermore, finite rate chemistry effects as well as non-adiabatic physics are considered.     

Numerical Investigation of a MILD Combustion Burner: Analysis of Mixing Field, Chemical Kinetics and Turbulence-Chemistry Interaction

A numerical study of a jet-in-hot coflow (JHC) burner emulating Moderate or Intense Low-oxygen Dilution (MILD) combustion conditions was carried out by solving the Reynolds Averaged Navier-Stokes equations in a two-dimensional axisymmetric domain and using the Eddy Dissipation Concept (EDC) for the turbulence-chemistry interaction treatment. A systematic methodology was used to analyze all possible sources of discrepancies observed between experimental and numerical data, trying to shedding light on the suitability of specific models for MILD combustion. In this regard, the deficiencies that may come from turbulence model or kinetic scheme have been shown by comparative study on four variants of the k-ε model (i.e. the standard, modified, realizable and RNG) together with the Reynolds stress model and three kinetic schemes namely KEE-58, DRM-19 and DRM-22. A variation of an EDC parameter (i.e. increasing the constant of the fine structure residence time) was proposed for better consideration of MILD combustion features and to overcome the over-prediction of peak temperature observed at downstream. In such a manner encouraging results were also obtained for the prediction of major combustion products as well as for CO and OH.     

The Numerical Simulation of Diesel Spray Combustion with LES-CMC

A Large Eddy Simulation (LES) approach together with the Conditional Moment Closure (CMC) method have been used for the simulation of spray combustion in engine-like conditions. The strategy consists of coupling an academic CMC code with the commercial CFD software Star-CD?(CD-adapco). Two issues have been investigated: firstly, the applicability of conventional spray models to LES and secondly, LES-CMC for spray combustion. Conventional spray models that were originally developed for use in Reynolds-averaged equations have been assessed for their applicability within the LES framework by conducting non-reacting spray computations. Liquid core penetration, spray spreading angle and vapour phase penetration have been compared to the available experimental data and the agreement between LES and experiments is satisfactory. Several reacting spray calculations have been performed with a range of initial mixture and temperature conditions, which mimic Diesel engine configurations. The computed auto-ignition time and flame lift-off length are in good agreement with the experimental data. Despite the uncertainties associated with the spray models and the chemistry, the results illustrate that the LES-CMC methodology can reproduce well the experimental results.     

Numerical Simulation of Two-Phase Inertial Flow in Heterogeneous Porous Media

In this study, non-Darcy inertial two-phase incompressible and non-stationary flow in heterogeneous porous media is analyzed using numerical simulations. For the purpose, a 3D numerical tool was fully developed using a finite volume formulation, although for clarity, results are presented in 1D and 2D configurations only. Since a formalized theoretical model confirmed by experimental data is still lacking, our study is based on the widely used generalized Darcy–Forchheimer model. First, a validation is performed by comparing numerical results of the saturation front kinetics with a semi-analytical solution inspired from the Buckley–Leverett model extended to take into account inertia. Second, we highlight the importance of inertial terms on the evolution of saturation fronts as a function of a suitable Reynolds number. Saturation fields are shown to have a structure markedly different from the classical case without inertia, especially for heterogeneous media, thereby, emphasizing the necessity of a more complete model than the classical generalized Darcy’s one when inertial effects are not negligible.     

Preferential Flow-Paths Detection for Heterogeneous Reservoirs Using a New Renormalization Technique

We have devised a renormalization scheme which allows very fast determination of preferential flow-paths and of up-scaled permeabilities of 2D heterogeneous porous media. In the case of 2D log-normal and isotropically distributed permeability-fields, the resulting equivalent permeabilities are very close to the geometric mean, which is in good agreement with a rigorous result of Matheron. It is also found to work well for geostatistically anisotropic media when comparing the resulting equivalent permeabilities with a direct solution of the finite-difference equations. The method works exactly as King's does, although the renormalization scheme was modified to obtain tensorial equivalent permeabilities using periodic boundary conditions for the pressure gradient. To obtain an estimation of the local fluxes, the basic idea is that if at each renormalization iteration all the intermediate renormalized permeabilities are stored in memory, we are able to compute -- ad reversum -- an approximation of the small-scale flux map under a given macroscopic pressure gradient. The method is very rapid as it involves a number of calculations that vary linearly with the number of elementary grid blocks. In this sense, the renormalization algorithm can be viewed as a rapid approximate pressure solver. The exact reference flow-rate map (for the finite-difference algorithm) was computed using a classical linear system inversion. It can be shown that the preferential flow paths are well detected by the approximate method, although errors may occur in the local flow direction.     

The Effect of Diluent on the Sustainability of MILD Combustion in a Cyclonic Burner

The present study investigates the characteristics of MILD/flameless combustion in a cyclonic lab-scale burner. Such a configuration is effective for achieving turbulent mixing in a very short time while allowing for a reasonably long residence time for the development of combustion reactions. These two constraints are mandatory in the case of MILD combustion processes (high inlet temperatures and diluted mixtures). Such operating conditions are achieved through massive heat/mass recirculation towards the fresh incoming mixtures by recycling the exhausted gases, featuring a process where chemical kinetics times are elongated because of the dilution levels. Thus, long residence times are needed to achieve a satisfying reaction progress, and the high inlet temperatures result in fast and efficient mixing between disproportionate flows to avoid the onset of oxidation reactions before achieving diluted conditions. Under these constraints, a lab-scale facility was designed and built. The oxidation processes of C₃H₈/O₂ mixtures highly diluted in N₂ or CO₂ were investigated by varying the external parameters of the system, namely, the inlet temperature (up to 1300 K) and the mixture composition (from lean to rich mixtures). Several combustion regimes were experimentally identified. When the MILD regime was established, the combustion process became homogeneous within the burner without luminous emissions. To investigate the distributed nature of the MILD combustion processes, chemical simulations were performed under the assumption of a well-stirred reactor. For both the diluents, good agreement between the experimental and numerical results was obtained for MILD combustion conditions.     

ALE方法在爆炸数值模拟中的应用

本文以ALE（Arbitrary Lagrange—Euler）理论为基础，结合非线性动力学的相关理论，推导了ALE方法描述的控制方程组。最后采用ALE描述方法针对集团装药在半无限土质中爆炸进行了数值模拟，并将模拟计算结果与现有实验研究成果进行了对比。在模拟给出的最大压力时程曲线和爆腔发展时程曲线的基础上，对空腔的形成和发展规律以及次生波的形成等问题进行了研究。通过对比，证明ALE算法综合了Lagrange和Euler法的优点，能够有效地用于对爆炸过程进行数值模拟。     

Numerical Simulation of Non-premixed Turbulent Combustion Using the Eddy Dissipation Concept and Comparing with the Steady Laminar Flamelet Model

Numerical simulations of the Sandia flame CHNa and the Sydney bluff-body stabilized flame HM1E are reported and the results are compared to available experimental data. The numerical method is based on compressible URANS formulations which were implemented recently in the OpenFOAM toolbox. In this study, the calculations are carried out using the conventional compressible URANS approach and a standard k- ?? turbulence model. The Eddy Dissipation Concept with a detailed chemistry approach is used for the turbulence-chemistry interaction. The syngas (CO/H₂) chemistry diluted by 30 % nitrogen in the Sandia flame CHNa and CH₄/H₂ combustion in the Sydney flame HM1E are described by the full GRI-3.0 mechanism. A robust implicit Runge-Kutta method (RADAU5) is used for integrating stiff ordinary differential equations to calculate the reaction rates. The radiation is treated by the P1-approximation model. Both target flames are predicted with the Steady Laminar Flamelet model using the commercial code ANSYS FLUENT as well. In general, there is good agreement between present simulations and measurements for both flames, which indicates that the proposed numerical method is suitable for this type of combustion, provides acceptable accuracy and is ready for further combustion application development.     

一种新型均匀扩散器流场的数值模拟

异性纤维检出设备中扩散器内部流场的情况,直接影响着落棉率和整个清除系统的工作效率.本文用计算流体力学软件FLUENT对新型均匀扩散器的流场进行数值模拟,通过改进扩散器的结构参数,实现流场的均匀和稳定,进而提高整个系统纤维检出设备的清除效率.之后对计算结果进行分析,提出了两种经过结构优化的改进模型,并分析了其内部流场.     

利用格子Boltzmann方法数值模拟Rayleigh-Benard对流

采用格子Boltzmann方法对较大Rayleigh数范围下的二维Rayleigh-Benard对流进行了模拟研究.引入能量分布函数,利用该能量分布函数与粒子速度分布函数耦合来求解一个热流场,能量分布函数与粒子速度分布函数和Boltzmann方程构成了一个新的双分布格子Boltzmann模型.在考虑密度随温度变化的情况下,进行数值模拟,得到了Rayleigh-Benard对流速度、温度随时间的变化规律、系统的流线和等温线分布及平均Nusselt数与Rayleigh数的之间的关系,与相关文献数据进行了对比,模拟结果非常吻合,证明了改进的双分布格子Boltzmann模型的有效性.     

Numerical Simulations of Flows in a Heterogeneous Porous Medium

Numerical simulations are employed to investigate the fluid flow and pressure loss in a heterogeneous block within a composite porous medium. The mean permeability of the heterogeneous block is seen to affect the overall effective flux significantly. The heterogeneous parameters of the permeability field, such as the correlation length and variance, affect it quite differently. Because of the channelling effects, the effective flux depends strongly on the realization of the permeability for larger correlation length. Under a specific permeability field, higher effective flux results from smaller variances. The influences of the inertial factor are found to be insignificant within the range of practical interests.     

H2O and CO2 Dilution in MILD Combustion of Simple Hydrocarbons

MILD combustion is a very attractive technology because of its intrinsic features for energy production from diluted gas deriving from bio- or thermochemical degradation of biomass. An effective use of such a technology for diluted fuel requires a thorough analysis of ignition and oxidation behavior to highlight the potential effects of the different fuel components on the basis of temperature and diluent/oxygen/fuel mixture composition. In this work, ignition and oxidation of a model gas surrogate for the gaseous fraction of biomass pyrolysis products containing C₁-C₂ species, CO and CO₂ were experimentally and numerically studied over a wide range of temperature and overall composition in the presence of large amounts of CO₂ or H₂O. Experimental results showed that such species significantly alter the evolution of the ignition process in dependence on temperature range and mixture composition. Several kinetic models were tested to simulate experimental results. Significant discrepancies occur, especially in the case of steam dilution. Numerical analyses suggested that such diluents acted mainly as third body species at low temperatures, conditioning both radical production pathways and the relative weight of C₁ oxidation/recombination routes, while strongly interacting with the H₂/O₂ high temperature branching mechanisms at high temperatures. Further analyses are mandatory to improve the predictability of the models and extend the applicability of the chemical schemes to non-standard conditions.     

Numerical Evaluation of Combustion Regimes in a GDI Engine

There is significant interest in the gasoline direct-injection engine due to its potential for improvements in fuel consumption but it still remains an area of active research due to a number of challenges including the effect of cycle-by-cycle variations. The current paper presents the use of a 3D-CFD model using both the RANS and LES turbulence modelling approaches, and a Lagrangian DDM to model an early fuel injection event, to evaluate the regimes of combustion in a gasoline direct-injection engine. The velocity fluctuations were investigated as an average value across the cylinder and in the region between the spark plug electrodes. The velocity fluctuations near the spark plug electrodes were seen to be of lower magnitude than the globally averaged fluctuations but exhibited higher levels of cyclic variation due to the influence of the spark plug electrode and the pent-roof geometry on the in-cylinder flow field. Differences in the predicted flame structure due to differences in the predicted velocity fluctuations between RANS and LES modelling approaches were seen as a consequence of the inherently higher dissipation levels present in the RANS methodology. The increased cyclic variation in velocity fluctuations near the spark plug electrodes in the LES predictions suggested significant variation in the relative strength of the in-cylinder turbulence and that may subsequently result in a thickening of the propagating flame front from cycle-to-cycle in this region. Throughout this paper, the numerical results were validated against published experimental data of the same engine geometry under investigation.     

Subgrid Scale Modeling in Large-Eddy Simulation of Turbulent Combustion Using Premixed Flamelet Chemistry

Large-eddy simulation (LES) of turbulent combustion with premixed flamelets is investigated in this paper. The approach solves the filtered Navier–Stokes equations supplemented with two transport equations, one for the mixture fraction and another for a progress variable. The LES premixed flamelet approach is tested for two flows: a premixed preheated Bunsen flame and a partially premixed diffusion flame (Sandia Flame D). In the first case, we compare the LES with a direct numerical simulation (DNS). Four non-trivial models for the chemical source term are considered for the Bunsen flame: the standard presumed beta-pdf model, and three new propositions (simpler than the beta-pdf model): the filtered flamelet model, the shift-filter model and the shift-inversion model. A priori and a posteriori tests are performed for these subgrid reaction models. In the present preheated Bunsen flame, the filtered flamelet model gives the best results in a priori tests. The LES tests for the Bunsen flame are limited to a case in which the filter width is only slightly larger than the flame thickness. According to the a posteriori tests the three models (beta-pdf, filtered flamelet and shift-inversion) show more or less the same results as the trivial model, in which subgrid reaction effects are ignored, while the shift-filter model leads to worse results. Since LES needs to resolve the large turbulent eddies, the LES filter width is bounded by a maximum. For the present Bunsen flame this means that the filter width should be of the order of the flame thickness or smaller. In this regime, the effects of subgrid reaction and subgrid flame wrinkling turn out to be quite modest. The LES-results of the second case (Sandia Flame D) are compared to experimental data. Satisfactory agreement is obtained for the main species. Comparison is made between different eddy-viscosity models for the subgrid turbulence, and the Smagorinsky eddy-viscosity is found to give worse results than eddy-viscosities that are not dominated by the mean shear. Paper presented on the Eccomas Thematic Conference Computational Combustion 2007, submitted for a special issue of Flow, Turbulence and Combustion.     

基于非定常N-S方程的三维机翼颤振特性研究

本研究采用N-S方程为控制方程,求解三维机翼的非定常气动力,并与颤振方程耦合迭代计算,求解机翼广义坐标响应的时间历程,根据广义坐标的时间历程的衰减、等幅和发散振荡等情况确定机翼颤振临界条件.以某导弹翼面作为算例,计算结果与实验结果和理论分析相吻合.     

Numerical investigation of the effects of turbulence on the ignition process in a turbulent MILD flame

Acta Mechanica Sinica - In this study, we conduct three-dimensional nonlinear large-eddy simulation to investigate the interaction between turbulence and reaction during the initial ignition...     

Numerical Simulation of a Turbulent Thermal

Ascent of a large-scale thermal in a standard atmosphere is calculated with the use of the Reynolds equations and the k model of turbulence, which takes into account temperature inhomogeneity and vorticity of the flow, and the Euler equations. Results of numerical calculations of a flow examined experimentally are presented. Gas-dynamic and turbulent flow parameters obtained in calculations and experiments are compared.     

基于有效参数法的固定床化学链燃烧数值模拟

本文基于有效参数法和变核颗粒模型,对固定床化学链燃烧进行数值模拟．计算并分析不同颗粒参数对二维固定床化学链燃烧结果的影响．在有效参数模型中,颗粒填充床被看作均匀的流体域．本文主要模拟了组分场和温度场,通过改变颗粒直径、孔隙率、内部粒子直径三组参数,取得了不同条件下二维固定床化学链燃烧床层所能达到的最高温度以及反应所用的时间．通过对不同参数模拟结果的分析,得出颗粒内部结构对反应特性的影响规律．