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.     

Large-Eddy Simulation of a Jet-in-Hot-Coflow Burner Operating in the Oxygen-Diluted Combustion Regime

Large-eddy simulations of moderate and intense low-oxygen dilution (MILD) combustion of a jet-in-hot-coflow (JHC) burner are performed. This burner configuration consists of three streams, providing fuel, oxygen-diluted coflow, and air to the burner. To account for the mixing between the three reactant streams, a three-stream flamelet/progress variable (FPV) formulation is utilized. This model was previously applied to a condition corresponding to the upper range of MILD-combustion, and the objective of this contribution is to further investigate this model in application to highly diluted operating conditions. Comparisons of mean and conditional results show that the model accurately captures effects of increasing oxygen-depletion on the flame-structure and heat-release, and predictions for temperature and species mass fractions are in overall good agreement with experiments.     

Influence of Inlet Dilution of Reactants on Premixed Combustion in a Recuperative Furnace

This paper presents a numerical study by RANS modeling that investigates the effect of external dilution on the premixed combustion occurring in a recuperative furnace. Calculations are performed using the detailed GRI-Mech 3.0 mechanism to ensure the accuracy of the modeling. Results of the in-furnace flow, temperature, and concentrations of OH, O₂, CO₂ and NO _x are provided. It is found that the external dilution with the inert gas CO₂ plays a significant role in establishing the premixed MILD (Moderate or Intense Low-oxygen Dilution) combustion. Externally diluting the reactant mixture not only reduces the initial concentration of O₂ but also ensures a stronger internal dilution by recirculation of more hot combustion products. Importantly, the latter effect is more significant for achieving the MILD regime. There is a critical mass fraction of the diluent CO₂ present, below which MILD combustion cannot occur. While the traditional premixed flame produces much more NO _x than the MILD combustion, the emission of NO _x appears to result most from the thermal-NO route and least from the N₂O route no matter which mode occurs. Moreover, the present simulation demonstrates that the MILD mode occurs over a wider range of initial reactant conditions for premixed combustion than for diffusion combustion.     

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.     

Modeling of Turbulent Natural Gas and Biogas Flames of the Delft Jet-in-Hot-Coflow Burner: Effects of Coflow Temperature,Fuel Temperature and Fuel Composition on the Flame Lift-Off Height

The structure of a turbulent non-premixed flame of a biogas fuel in a hot and diluted coflow mimicking moderate and intense low dilution (MILD) combustion is studied numerically. Biogas fuel is obtained by dilution of Dutch natural gas (DNG) with CO₂. The results of biogas combustion are compared with those of DNG combustion in the Delft Jet-in-Hot-Coflow (DJHC) burner. New experimental measurements of lift-off height and of velocity and temperature statistics have been made to provide a database for evaluating the capability of numerical methods in predicting the flame structure. Compared to the lift-off height of the DNG flame, addition of 30 % carbon dioxide to the fuel increases the lift-off height by less than 15 %. Numerical simulations are conducted by solving the RANS equations using Reynolds stress model (RSM) as turbulence model in combination with EDC (Eddy Dissipation Concept) and transported probability density function (PDF) as turbulence-chemistry interaction models. The DRM19 reduced mechanism is used as chemical kinetics with the EDC model. A tabulated chemistry model based on the Flamelet Generated Manifold (FGM) is adopted in the PDF method. The table describes a non-adiabatic three stream mixing problem between fuel, coflow and ambient air based on igniting counterflow diffusion flamelets. The results show that the EDC/DRM19 and PDF/FGM models predict the experimentally observed decreasing trend of lift-off height with increase of the coflow temperature. Although more detailed chemistry is used with EDC, the temperature fluctuations at the coflow inlet (approximately 100K) cannot be included resulting in a significant overprediction of the flame temperature. Only the PDF modeling results with temperature fluctuations predict the correct mean temperature profiles of the biogas case and compare well with the experimental temperature distributions.     

Flamelet-Modeling of Inverse Rich Diffusion Flames

An experimental and numerical investigation of a confined laminar inverse diffusion flame (IDF) with pure oxygen as oxidizer and carbon dioxide diluted methane as fuel with a global stoichiometry of partial oxidation processes (equivalence ratio of 2.5) is presented. The present burner setup allows studying both the flame and the post-flame zone in a simplified geometry considering typical operating conditions as found in large-scale gasifiers. This partial oxidation flame setup is characterized by very high temperatures close to the stoichiometric oxidation zone due to oxy-fuel combustion, whereas lower temperatures and slow endothermic post-flame conversion reactions with long residence times are found in the fuel rich post-flame region. The scope of this paper is to investigate different modeling approaches suitable for both regimes by comparing the simulation results to detailed experimental data. Planar OH laser-induced fluorescence (OH-LIF) was performed for measuring the hydroxyl radical in the reaction zone and the results are compared to CFD calculations. Based on this comparison, the necessary level of detail of diffusion flux modeling, which includes Soret and Dufour effects, is analyzed and established. Finally, steady and unsteady non-premixed flamelet approaches based on a single mixture fraction are used in order to study their applicability for both the oxidation and post-flame zone. Significantly different time scales are obtained using different flamelet paths. Their influence on the results is investigated in the steady flamelet and the Lagrangian flamelet approach.     

Computational investigation of hydrodynamics,coal combustion and NOx emissions in a tangentially fired pulverized coal boiler at various loads

This work presents a computational investigation of hydrodynamics, coal combustion and NOx emissions in a tangentially fired pulverized coal boiler at different loads (630, 440 and 300 MW; relative loads of 100%, 70% and 48%) to clarify the effect of load change on the furnace processes. A computational fluids dynamics model was established; the flow field, temperature profile, species concentration and NOx emissions were predicted numerically; and the influence of burner tilt angles was evaluated. Simulation results indicate that a decrease in boiler load decreases the gas velocity, attenuates the airflow rotations, and increases the tangent circle size. The high-temperature zone and flame moved toward the side walls. Such behaviors impair air–fuel mixing, heat transfer and steady combustion in the furnace. In terms of species concentrations, a decrease in boiler load increased the O₂ content, decreased the CO content, and decreased the char burnout rates only slightly. A change in boiler load from 630 to 440 and 300 MW increased the NOx emissions from 202 to 234 and 247 mg/m³, respectively. Burner tilt angles are important in coal combustion and NOx emissions. A burner angle of –15° favors heat transfer and low NOx emissions (<185 mg/m³) for the current tangentially fired boiler.     

Characteristics of the reaction zone in a combustor operating at mild combustion

The aerodynamic, chemical and thermal aspects of the mild combustion process have been studied with emphasis on mixing rates, flue gas recirculation and strong shear produced by reactants supplied from discrete jets. Time-averaged and instantaneous structures of turbulent flow were examined by visualization and local measurements within a 5400 W burner operating with methane with an overall equivalence ratio varying from 0.8 to 1.2 and at non-premixed and premixed modes. The results showed that the entrainment of the flue gases into the fresh mixture was very important for the initiation and progress of the reaction, and occurred in two successive mechanisms. Initially, the flue gases were driven with the reverse flow towards the annular exit where, by Biot–Savart induction, they acquired some momentum from the supply streams provided at the center. The resulting mixing process in the close vicinity of the burner was less intermittent and this was evident in relatively lower values of the second order moments of the residence time distribution. Slightly downstream, the second order moments were, however, increased by large-scale turbulent fluctuations and this led to the enhancement of the mixing process and introduced some further intermittency. The latter entrainment mechanism caused the flue gases to partially encapsulate the discrete jets, which resulted in islands of flammable mixture surrounded by the inert gases. Hence, as the instantaneous OH radical visualizations revealed, the reaction was only initiated away from the burner and in disconnected regions where the Rayleigh pictures showed strong temperature gradients. As the distance from the nozzle increased further, the reaction seemed to follow local flow patterns in that it progressed radially outwards with large structures, which resulted in an increased space-averaged temperature. Furthermore, the residence time decreased away from the burner and the flame came close to extinction due to the high stretching rates of the large structures. However, the flue gases entrained up to this point increased the inert content of the fresh mixture with chemical time scales comparable to the time scales of the flow. This allowed the reactants to attain temperatures near to those of the flue gases and to ignite with a small temperature rise, which led to a much lower thermal NO formation. The results also showed that when the equivalence ratio of the non- premixed mixtures was increased, the region where the combustion took place was shifted away from the burner and extended further downstream towards the roof. In the case of premixed combustion, however, the reaction started and terminated earlier and was confined to regions in close proximity to the axis. The emissions of OH radical occurred rather patchily and in relatively high concentrations. Received: 6 June 2000/Accepted: 28 October 2000     

2-methylfuran Oxidation in the Absence and Presence of NO

2-methylfuran (2-MF) has become of interest as biofuel because of its properties and the improvement in its production method, and also because it is an important intermediate in the conversion of 2,5-dimethylfuran. In this research, an experimental and kinetic modelling study of the oxidation of 2-MF in the absence and presence of NO has been performed in an atmospheric pressure laboratory installation. The experiments were performed in a flow reactor and covered the temperature range from 800 to 1400 K, for mixtures from very fuel-rich to very fuel-lean, highly diluted in nitrogen. The inlet 2-MF concentration was 100 ppm. In the experiments in the presence of NO, the inlet NO concentration was 900 ppm. An interpretation of the experimental results was performed through a gas-phase chemical kinetic model. A reasonable agreement between the experimental trends and the modelling data is obtained. The results of the concentration profile of 2-MF as a function of temperature indicate that, both in the absence and in the presence of NO, the onset of 2-MF consumption is shifted to lower temperatures only under fuel-lean and very fuel-lean conditions. Furthermore, under these conditions the presence of NO also shifts the onset of 2-MF consumption to lower temperatures. The effect of the 2-MF presence on the NO reduction varies with the oxygen concentration. It is seen that under very fuel-rich and stoichiometric conditions NO is reduced basically by reburn reactions, while under fuel-lean and very fuel-lean conditions, the NO-NO₂ interconversion appears to be dominant.     

Ultra-low-emission steam boiler constituted of reciprocal flow porous burner

This experimental study examined a low-emission steam boiler in which the filtration combustion technology was employed. This new boiler concept is consisted of a reciprocal flow porous burner, in which a combustion wave propagates along the reactor length. The boiler’s burner is filled up by an inert porous material, which leads to a stable burning of ultra-lean fuel/air mixtures, operating below flammability limits of conventional burners. In reciprocal filtration combustion, the reaction zone travels back and forth along the length of the burner, maintaining a typical trapezoidal temperature distribution favorable to the energy extraction. Embedding heat exchangers into the ends of the porous bed results in an alternative low-emission high-efficiency boiler. The heat re-circulation inside the porous matrix and the low degree of thermal non-equilibrium between the gas and the solid phases result in ultra-low levels of CO and NO_x. Over an equivalence ratio range from 0.20 to 1.0 and a gas flow velocity range from 0.2 to 0.6 m/s, burning the technical methane, the developed prototype has reached efficiencies superior to 90% and NO_x and CO emission levels lower than 1.0 and 0.5 ppm, respectively.     

Large Eddy simulation of turbulent hydrogen-fuelled supersonic combustion in an air cross-flow

The main aim of this article is to provide a theoretical understanding of the physics of supersonic mixing and combustion. Research in advanced air-breathing propulsion systems able to push vehicles well beyond $M=4$ M = 4 is of interest around the world. In a scramjet, the air stream flow captured by the inlet is decelerated but still maintains supersonic conditions. As the residence time is very short $(\sim \!\!\mathrm{1ms})$ ( ～ 1 ms ) , the study of an efficient mixing and combustion is a key issue in the ongoing research on compressible flows. Due to experimental difficulties in measuring complex high-speed unsteady flowfields, the most convenient way to understand unsteady features of supersonic mixing and combustion is to use computational fluid dynamics. This work investigates supersonic combustion physics in the Hyshot II combustion chamber within the Large Eddy simulation framework. The resolution of this turbulent compressible reacting flow requires: (1) highly accurate non-dissipative numerical schemes to properly simulate strong gradients near shock waves and turbulent structures away from these discontinuities; (2) proper modelling of the small subgrid scales for supersonic combustion, including effects from compressibility on mixing and combustion; (3) highly detailed kinetic mechanisms (the Warnatz scheme including 9 species and 38 reactions is adopted) accounting for the formation and recombination of radicals to properly predict flame anchoring. Numerical results reveal the complex topology of the flow under investigation. The importance of baroclinic and dilatational effects on mixing and flame anchoring is evidenced. Moreover, their effects on turbulence-scale generation and the scaling law are analysed.     

Heat Transfer Study of the Hencken Burner Flame

The Hencken burner flame is often used in combustion laser diagnostics as a calibration flame because of its near adiabatic condition. For a fast burning H₂ flame, it can tolerate high flow rate and the flame is indeed near adiabatic; however, for a slow burning CH₄ flame, the flow rate is not always high enough to maintain near adiabatic conditions. The heat transfer of the H₂ and CH₄ Hencken burner flames are studied numerically and experimentally. Three heat loss mechanisms are analyzed: the burner surface radiation, the hot gas radiation, and the convection heat transfer between the main flow and the co-flow. The surface radiation produces negligible temperature drop while the gas radiation and the convection heat loss contribute significant temperature drop. Reducing the co-flow rate can decrease the convection heat loss slightly. The temperature drop caused by the heat loss is inversely proportional to the main flow rate. Increasing the burner size and running the flame premixed mode can increase the flow rate and reduce the temperature deviation from the adiabatic equilibrium value. Based on the heat loss and temperature drop analysis, suggestions are given to maintain the flame at near adiabatic conditions.     

Numerical Simulation of Pulverized Coal MILD Combustion Using a New Heterogeneous Combustion Submodel

The scope of this investigation is the application and analysis of a recently developed submodel (Schulze et al., Oil Gas Science Technol, 2013, doi:10.2516/ogst/2012069) for char particle combustion and gasification. The distinguishing feature of this model is a detailed representation of the diffusion and convection processes as well as the homogeneous reactions in the boundary layer around the char particle. These processes are fully coupled to the heterogeneous particle kinetics. The model was implemented into the CFD code ANSYS-Fluent. The coupled solver is used for simulating the IFRF full scale pulverized coal combustion MILD furnace, for which detailed experimental data are available for model evaluation (Orsino et al., IFRF Doc. No F46/y/3, 2000) The new model yields improved agreement with measured data as compared to the standard modeling approach. This can be directly related to the prediction of the char burnout rate. For further analysis, the mixing field in the IFRF furnace is investigated in detail by introducing four mixture fractions for pyrolysis products, char burn-off gases, primary and secondary air, respectively. The solutions of the respective transport equations are used to define the local stoichiometry both in the gas phase and on the particle surface in such a multi-stream system. The conditions in the particle surrounding gas phase as well as on the particle surface are used to define the regime of particle-gas interaction based on the simulations with the new submodel. It can be shown that for certain conditions the homogeneous reactions in the particle boundary must be accounted for.     

Experimental and modeling investigation of the effect of air preheat on the formation of NOx in an RQL combustor

The Rich-burn/Quick–mix/Lean-burn (RQL) combustor concept has been proposed to minimize the formation of oxides of nitrogen (NO_x) in gas turbine systems. The success of this low-NO_x combustor strategy is dependent upon the links between the formation of NO_x, inlet air preheat temperature, and the mixing of the jet air and fuel-rich streams. Chemical equilibrium and kinetics modeling calculations and experiments were performed to further understand NO_x emissions in an RQL combustor. The results indicate that as the temperature at the inlet to the mixing zone increases (due to preheating and/or operating conditions) the fuel-rich zone equivalence ratio must be increased to achieve minimum NO_x formation in the primary zone of the combustor. The chemical kinetics model illustrates that there is sufficient residence time to produce NO_x at concentrations that agree well with the NO_x measurements. Air preheat was found to have very little effect on mixing, but preheating the air did increase NO_x emissions significantly. By understanding the mechanisms governing NO_x formation and the temperature dependence of key reactions in the RQL combustor, a strategy can be devised to further reduce NO_x emissions using the RQL concept.     

LES of a Multi-burner Annular Gas Turbine Combustor

In this study, Large Eddy Simulation (LES) has been used to predict the flow, mixing and combustion in both a single burner laboratory gas turbine combustor and in an 18 burner annular combustor, having identical cross sections. The LES results for the single burner laboratory combustor are compared with experimental data for a laboratory model of this combustor, and with other LES predictions, with good agreement. An explicit finite volume based LES model, using the mixed subgrid model together with a partially stirred reactor model for the turbulence chemistry interactions, is used. For the annular combustor, with the swirlers parameterized by jet inflow boundary conditions, we have investigated the influence of the a-priori unknown combustor exit impedance, the influence of the swirler characteristics and the fuel type. The combustion chemistry of methane–air and n-decane–air combustion is modeled by a two-step reaction mechanism, whereas NOx is separately modeled with a one-step mechanism. No experimental data exists for the annular combustor, but these results are compared with the single burner LES and experimental results available. The combustor exit impedance, the swirler- and fuel characteristics all seem to influence the combusting flow through the acoustics of the annular combustor. To examine this in greater detail time-series and eigenmodes of the combustor flow fields are analyzed and comparisons are made also with results from conventional thermoacoustic eigenmode analysis, with reasonable agreement. The flow and pressure distributions in the annular combustor are described in some detail and the mechanisms by which the burners interact are outlined.     

A multi-phase theory for the detonation of granular explosives containing an arbitrary number of solid components

Multiphase continuum models are commonly used to predict the shock, combustion and detonation behavior of granular energetic mixtures containing solid reactants and gaseous products. These models often include phase interaction terms that formally satisfy the strong form of the Second Law of Thermodynamics and provide flexibility in distributing dissipation between phases arising from non-equilibrium phenomena. This work presents a thermodynamically compatible constitutive theory for reactive systems containing an arbitrary number of solid components. The theory represents a rigorous extension of the two-phase theory formulated by Bdzil et al., based on the well-studied Baer–Nunziato model. Forms of the gas–solid and solid–solid interphase sources suggested by general reactions of type A → B are considered, where the combustion processes discussed in Bdzil et al. are treated as a special case. The model energetics are augmented by supplemental evolutionary equations that track local changes in phase temperatures due to dissipative and transport processes allowing for the identification of dominant energetic processes. This capability provides a mean to identify system parameters (e.g., metal particle size and mass fraction in metalized energetic mixtures) which optimize performance metrics. Detonation predictions are given for mixtures of granular HMX and aluminum to demonstrate model features and to highlight the effect of aluminum particle self-heating by oxidation on detonation. Predicted spatial profiles for mechanical fields, and the heating contributions from individual dissipative processes, illustrate how aluminum particle size can affect the coupling of oxidative heating to the explosive reaction zone.     

Numerical Simulation of Delft-Jet-in-Hot-Coflow (DJHC) Flames Using the Eddy Dissipation Concept Model for Turbulence–Chemistry Interaction

In this paper, we report results of a numerical investigation of turbulent natural gas combustion for a jet in a coflow of lean combustion products in the Delft-Jet-in-Hot-Coflow (DJHC) burner which emulates MILD (Moderate and Intense Low Oxygen Dilution) combustion behavior. The focus is on assessing the performance of the Eddy Dissipation Concept (EDC) model in combination with two-equation turbulence models and chemical kinetic schemes for about 20 species (Correa mechanism and DRM19 mechanism) by comparing predictions with experimental measurements. We study two different flame conditions corresponding to two different oxygen levels (7.6% and 10.9% by mass) in the hot coflow, and for two jet Reynolds number (Re = 4,100 and Re = 8,800). The mean velocity and turbulent kinetic energy predicted by different turbulence models are in good agreement with data without exhibiting large differences among the model predictions. The realizable k-ε model exhibits better performance in the prediction of entrainment. The EDC combustion model predicts too early ignition leading to a peak in the radial mean temperature profile at too low axial distance. However the model correctly predicts the experimentally observed decreasing trend of lift-off height with jet Reynolds number. A detailed analysis of the mean reaction rate of the EDC model is made and as possible cause for the deviations between model predictions and experiments a low turbulent Reynolds number effect is identified. Using modified EDC model constants prediction of too early ignition can be avoided. The results are weakly sensitive to the sub-model for laminar viscosity and laminar diffusion fluxes.