Second-order conditional moment closure modeling of a turbulent CH4/H2/N2 jet diffusion flame

The second-order CMC model for a detailed chemical mechanism is used to model a turbulent CH₄/H₂/N₂ jet diffusion flame. Second-order corrections are made to the three rate limiting steps of methane–air combustion, while first-order closure is employed for all the other steps. Elementary reaction steps have a wide range of timescales with only a few of them slow enough to interact with turbulent mixing. Those steps with relatively large timescales require higher-order correction to represent the effect of fluctuating scalar dissipation rates. Results show improved prediction of conditional mean temperature and mass fractions of OH and NO. Major species are not much influenced by second-order corrections except near the nozzle exit. A parametric study is performed to evaluate the effects of the variance parameter in log-normal scalar dissipation PDF and the constants for the dissipation term in conditional variance and covariance equations.     

Rotational CARS for simultaneous measurements of temperature and concentrations of N2, O2, CO, and CO2 demonstrated in a CO/air diffusion flame

Rotational coherent anti-Stokes Raman spectroscopy (CARS) has over the years demonstrated its strong potential to measure temperature and relative concentrations of major species in combustion. A recent work is the development and experimental validation of a CO₂ model for thermometry, in addition to our previous rotational CARS models for other molecules. In the present work, additional calibration measurements for relative CO₂/N₂ concentrations have been made in the temperature range 294-1246 K in standardized CO₂/N₂ mixtures. Following these calibration measurements, rotational CARS measurements were performed in a laminar CO/air diffusion flame stabilized on a Wolfhard-Parker burner. High-quality spectra were recorded from the fuel-rich region to the surrounding hot air in a lateral cross section of the flame. The spectra were evaluated to obtain simultaneous profiles of temperature and concentrations of all major species; N₂, O₂, CO, and CO₂. The potential for rotational CARS as a multi-species detection technique is discussed in relation to corresponding strategies for vibrational CARS.     

Dimensional scaling of flame propagation in discrete particulate clouds

The critical dimension necessary for a flame to propagate in suspensions of fuel particles in oxidiser is studied analytically and numerically. Two types of models are considered: First, a continuum model, wherein the individual particulate sources are not resolved and the heat release is assumed spatially uniform, is solved via conventional finite difference techniques. Second, a discrete source model, wherein the heat diffusion from individual sources is modelled via superposition of the Green's function of each source, is employed to examine the influence of the random, discrete nature of the media. Heat transfer to cold, isothermal walls and to a layer of inert gas surrounding the reactive medium are considered as the loss mechanisms. Both cylindrical and rectangular (slab) geometries of the reactive medium are considered, and the flame speed is measured as a function of the diameter and thickness of the domains, respectively. In the continuum model with inert gas confinement, a universal scaling of critical diameter to critical thickness near 2:1 is found. In the discrete source model, as the time scale of heat release of the sources is made small compared to the interparticle diffusion time, the geometric scaling between cylinders and slabs exhibits values greater than 2:1. The ability of the flame in the discrete regime to propagate in thinner slabs than predicted by continuum scaling is attributed to the flame being able to exploit local fluctuations in concentration across the slab to sustain propagation. As the heat release time of the sources is increased, the discrete source model reverts back to results consistent with the continuum model. Implications of these results for experiments are discussed.     

Combustion properties of H2/N2/O2/steam mixtures

The present work reports new experimental and numerical results of the combustion properties of hydrogen based mixtures diluted by nitrogen and steam. Spherical expanding flames have been studied in a spherical bomb over a large domain of equivalence ratios, initial temperatures and dilutions at an initial pressure of 100 kPa (T_ini = 296, 363, 413 K; N₂/O₂ = 3.76, 5.67, 9; %Steam = 0, 20, 30). From these experiments, the laminar flame speed $S_L^0$, the Markstein length L’, the activation energy Ea and the Zel'dovich β number have been determined. These parameters were also simulated using COSILAB^® in order to verify the validity of the Mével et al. [1] detailed kinetic mechanism. Other parameters as the laminar flame thickness δ and the effective Lewis number Le_eff were also simulated. These new results aim at providing an extended database that will be very useful in the hydrogen combustion hazard assessment for nuclear reactor power plant new design.     

Effect of CO2 gasification reaction on oxy-combustion of pulverized coal char

For oxy-combustion with flue gas recirculation, as is commonly employed, it is recognized that elevated CO₂ levels affect radiant transport, the heat capacity of the gas, and other gas transport properties. A topic of widespread speculation has concerned the effect of the CO₂ gasification reaction with coal char on the char burning rate. To give clarity to the likely impact of this reaction on the oxy-fuel combustion of pulverized coal char, the Surface Kinetics in Porous Particles (SKIPPY) code was employed for a range of potential CO₂ reaction rates for a high-volatile bituminous coal char particle (130 μm diameter) reacting in several O₂ concentration environments. The effects of boundary layer chemistry are also examined in this analysis. Under oxygen-enriched conditions, boundary layer reactions (converting CO to CO₂, with concomitant heat release) are shown to increase the char particle temperature and burning rate, while decreasing the O₂ concentration at the particle surface. The CO₂ gasification reaction acts to reduce the char particle temperature (because of the reaction endothermicity) and thereby reduces the rate of char oxidation. Interestingly, the presence of the CO₂ gasification reaction increases the char conversion rate for combustion at low O₂ concentrations, but decreases char conversion for combustion at high O₂ concentrations. These calculations give new insight into the complexity of the effects from the CO₂ gasification reaction and should help improve the understanding of experimentally measured oxy-fuel char combustion and burnout trends in the literature.     

An intrinsic kinetics model to predict complex ash effects (ash film,dilution, and vaporization) on pulverized coal char burnout in air (O2/N2) and oxy-fuel (O2/CO2) atmospheres

During coal combustion, char chemical reaction is the slowest step, particularly in the last burnout stage, where the char consists of small amounts of carbon in a predominant ash framework. Existing kinetics models tend to deviate from experimental measurements of late char burnout due to the incomplete treatment of ash effects. Ash can improve pore evolution through vaporization, hinder oxygen transport by forming an ash film, and reduce active carbon sites and available surface per unit volume by penetrating into the char matrix. In this work, a sophisticated kinetics model, focusing on these three ash evolution mechanisms (ash vaporization, ash film, and ash dilution) during pulverized coal (PC) char combustion, is developed by integrating them into a thorough mechanistic picture. Further, a detailed comparison of the three distinct ash effects on PC char conversion during air (O₂/N₂) and oxy-fuel (O₂/CO₂) combustion is performed. For the modeled coal, the mass of ash vaporization is approximate 3 orders less than the mass of ash remaining, which participates in ash dilution and ash film formation, both in O₂/N₂ and O₂/CO₂ atmospheres. The influence of these phenomena on burnout time follows the order: ash dilution > ash film > ash vaporization. The influence of ash vaporization on burnout time is minor, but through interactions with the ash dilution and ash film forming processes it can have an impact at high extents of burnout, particularly in O₂/CO₂ atmospheres. In O₂/N₂ atmospheres the residual ash predominately exists as an ash film, whereas it mainly exists as diluted ash in the char matrix in O₂/CO₂ atmospheres. The residual ash particle is encased by a thick film when the ash film forming fraction is high (low ash dilution fraction). These results provide in-depth insights into the conversion of PC char and further utilization of the residual ash.     

Burning rates of turbulent iso-octane aerosol mixtures in spherical flame explosions

Experimental studies of aerosol combustion under quiescent and turbulence conditions have been conducted to quantify the differences in the flame structure and burning rates between aerosol and gaseous mixtures. Turbulence was generated by variable speed fans to yield rms turbulence velocities between 0.5 and 4.0 m/s and this was uniform and isotropic. Homogeneously distributed and near monodispersed iso-octane-air aerosol clouds were generated using a thermodynamic condensation method. Spherically expanding flames, following central ignition, at near atmospheric pressures were employed to quantify the flame structure and propagation rate. The effects of the diameter of fine fuel droplets on flame propagation were investigated. It is suggested that the inertia of fuel droplets is an important cause of flame enhancement during early flame development. During later stages, cellular flame instability and the effective, gaseous phase, equivalence ratio becomes important. The latter effect leads has increases the flame speed of rich mixtures, but decreases that of lean ones. Droplet enhancement of burning velocity can be significant at low turbulence but is negligible at high turbulence.     

Relay propagation of partially coherent flattened Gaussian beam in turbulent atmosphere

Relay propagation of partially coherent flattened Gaussian–Schell beam in turbulent atmosphere has been studied. The analytical expresses of average intensity distribution at target are derived. The effects of spatial correlation length of initial and relay beam on the average intensity are analyzed in details. Study shows that the effects of the variation of spatial correlation length of relay beam are much larger than that of initial beam. The effects of spatial correlation length on relay propagation become smaller and smaller with the increase of structure constant. When the spatial correlation length is large and turbulence is strong, the effects of the variation of spatial correlation length on relay propagation are so small that can be neglected.     

Flame stabilization in high pressure LOx/GH2 and GCH4 combustion

Planar laser induced fluorescence (PLIF) of OH is used to examine flame stabilization in high pressure cryogenic flames formed by injecting a central jet of low speed liquid oxygen surrounded by a high speed gaseous stream of hydrogen or methane. In the LO_x/GH₂ experiments injection conditions are transcritical as the chamber pressure is above critical but the temperature is below critical . In the LO_x/GCH₄ experiments the chamber pressure and LO_x injection temperature are below critical , . Hydrogen or methane are injected at room temperature LIF images delineate the flame edge in the injector nearfield. The two flames are stabilized in the vicinity of the liquid oxygen injector lip but the anchor point is found to lie closer to the lip in the LO_x/GH₂ case and its displacement from shot to shot is of a smaller amplitude than that corresponding to the LO_x/GCH₄ flame. Interpretation of these data is based on a previous analysis which indicates that stabilization is essentially controlled by a dimensionless group formed by comparing the lip thickness to the flame edge thickness Ψ = h_s/δ_f. It is found that Ψ slightly exceeds unity in the LO_x/GH₂ case essentially fulfilling the stability condition while Ψ < 1 in the LO_x/GCH₄ case. In this last situation the flame is thicker than the characteristic thickness h_s and it is therefore sensitive to the high speed methane stream. Anchoring is imperfect and the flame edge moves with the turbulent eddies shed from the lip. Global stabilization is achieved dynamically but the reactive layer is not well established and the large amplitude motion of the edge is a symptom of a possible lift-off. Theoretical estimates indicate that LO_x/GCH₄ flame stabilization requires a thicker lip size than the LO_x/GH₂ propellant couple.     

Measurements of three-dimensional mean flame surface area ratio in turbulent premixed Bunsen flames

A data processing scheme with particular emphasis on proper flame contour smoothing is developed and applied to measure the three-dimensional mean flame surface area ratio in turbulent premixed flames. The scheme is based on the two-sheet imaging technique such that the mean flame surface area ratio is an average within a window covering a finite section of the turbulent flame brush. This is in contrast to the crossed-plane tomograph technique which applies only to a line. Two sets of Bunsen flames have been investigated in this work with the turbulent Reynolds number up to 4000 and the Damköhler number ranging from less than unity to close to 10. The results show that three-dimensional effects are substantial. The measured three-dimensional mean flame surface area ratio correlates well with a formula similar to the Zimont model for turbulent burning velocity but with different model constants. Also, the mean flame surface area ratio displays a weak dependency on turbulence intensity but a strong positive dependency on the turbulence integral length scale.     

Evaluation of a flamelet/progress variable approach for pulverized coal combustion in a turbulent mixing layer

A steady flamelet/progress variable (FPV) approach for pulverized coal flames is employed to simulate coal particle burning in a turbulent shear and mixing layer. The configuration consists of a carrier-gas stream of air laden with coal particles that mixes with an oxidizer stream of hot products from lean combustion. Carrier-phase DNS (CP-DNS) are performed, where the turbulent flow field is fully resolved, whereas the coal is represented by Lagrangian point particles. CP-DNS with direct chemistry integration is performed first and provides state-of-the-art validation data for FPV modeling. In a second step the control variables for FPV are extracted from the CP-DNS and used to test if the tabulated manifold can correctly describe the reacting flow (a priorianalysis). Finally a fully coupled a posteriori FPV simulation is performed, where only the FPV control variables are transported, and the chemical state is retrieved from the table and fed back to the flow solver. The a priori results show that the FPV approach is suitable for modeling the complex reacting multiphase flow considered here. The a posteriori data is similarly in good agreement with the reference CP-DNS, although stronger deviations than a priori can be observed. These discrepancies mainly appear in the upper flame (of the present DNS), where premixing and highly unsteady extinction and re-ignition effects play a role, which are difficult to capture by steady non-premixed FPV modeling. However, the present FPV model accurately captures the lower, more stable flame that burns in non-premixed mode.     

Ignition enhancement by addition of NO and NO2 from a N2/O2 plasma torch in a supersonic flow

The effects of NO and NO₂ produced by using a plasma jet (PJ) of a N₂/O₂ mixture on ignition of hydrogen, methane, and ethylene in a supersonic airflow were experimentally and numerically investigated. Numerical analysis of ignition delay time showed that the addition of a small amount of NO or NO₂ drastically reduced ignition delay times of hydrogen and hydrocarbon fuels at a relatively low initial temperature. In particular, NO and NO₂ were more effective than O radicals for ignition of a CH₄/air mixture at 1200 K or lower. These ignition enhancement effects were examined by including the low temperature chemistry. Ignition tests by a N₂/O₂ PJ in a supersonic flow (M = 1.7) for using hydrogen, methane, and ethylene injected downstream of the PJ were conducted. The results showed that the ignitability of the N₂/O₂ PJ is affected by the composition of the feedstock and that pure O₂ is not the optimum condition for downstream fuel injection. This result of ignition tests with downstream fuel injection demonstrated a significant difference in ignition characteristics of the PJ from the ignition tests with upstream fuel injection.     

Wing profile measurements of the Na-doublet lines in C2H2/O2/N2, H2/O2/N2 and H2/O2/Ar frames at 1 atm

Fluorescence-excitation (wing) profiles of the Na-D doublet lines were measured over a wavelength range extending from 0.3 to 200 Å from the line center for the red D₁ and blue D₂ wings and from 0.3 to 3 Å for the red D₂ and the blue D₁ wings, respectively. The line profiles were determined with the aid of a tunable CW dye-laser as a background source by measuring the total fluorescence intensity observed on detuning the laser wavelength. The flames were premixed, laminar, shielded flames at 1 atm, with temperatures ranging from 1860 to 2270 K; N₂ and Ar served as diluent gases. The line core and near-wing profiles (i.e. the region covering 0.3<Δλ<7 Å for the outer wings and 0.3<Δλ<3 Å for the inner ones) in all of the flames studied appeared to have the same frequency dependence, regardless of the nature and concentrations of the gases used. The blue D₂-line profile followed an unexpected (-2.2) law, while the other three profiles obeyed the theoretically expected (-2) law (the dispersion profile function). The line profile in the Δλ range between the impact and quasistatic regions was found to depend on the main perturbers involved. We found that the far blue D₂- and red D₁-wings in the Ar-diluted H₂/O₂ flame obeyed the (-$\text{5}\text{4}$) and (-$\text{3}\text{2}$) laws, respectively, as predicted by the quasi-static theory for the Lennard-Jones interaction. For the N₂-diluted C₂H₂/O₂ and H₂/O₂ flames, we did not find these wing dependences in the Δλ range investigated.     

Ignition of a sequential combustor: Evidence of flame propagation in the autoignitable mixture

Ignition of the second stage in a lab-scale sequential combustor is investigated experimentally. A fuel mixing section between jet-in-cross-flow injection and the second stage chamber allows the fuel and vitiated, hot cross-flow to partially mix upstream of the main heat release zone. The focus of the present work is on the transient ignition process leading to a stable flame in the second stage. High-speed OH-PLIF as well as OH chemiluminescence imaging is applied to obtain complementary planar and line-of-sight integrated information on the ignition. We find experimental evidence for the co-existence of two regimes dominating the chamber ignition, i.e. autoignition and flame propagation. As the mass flow of the dilution air injected downstream of the first stage is increased (i.e. mixing temperatures in the fuel mixing section are decreased), we transition from an autoignition to a flame propagation dominated regime. Hysteresis in the ignition behavior is observed indicating that the first stage in a sequential combustor may be operated at leaner conditions than required for ignition of the second stage. The time traces of integral heat release obtained simultaneously with a photomultiplier tube show distinct features depending on the dominating regime, which is important for high-pressure testing with limited optical access.     

Slanting propagation of stochastic electromagnetic twist anisotropic Gaussian–Schell model beams in turbulent atmosphere

Based on the extended Huygens–Fresnel diffraction integral and the unified theory of coherence and polarization, the slanting propagation properties of stochastic electromagnetic twist anisotropic Gaussian–Schell model (ETAGSM) beams in turbulent atmosphere have been investigated. Both numerical calculation and physical interpretation are obtained. The influence of the atmospheric turbulence and the source parameters on the polarization distribution, intensity distribution and coherence distribution has been studied in great detail. The investigation reveals that slanting propagation can weaken the influence of atmospheric turbulence on the spectral properties of stochastic ETAGSM beams more obviously. It is also shown that we can modulate the spectral degree of polarization in the output plane especially in the far field by simply controlling the coherence, the spot width or the twisty properties of the source beams. The distribution of the spectral degree of polarization and cross-polarization of the beams in the output plane is also given out.     

Prediction of local extinction and re-ignition effects in non-premixed turbulent combustion using a flamelet/progress variable approach

The flamelet/progress variable approach (FPVA) has been proposed by Pierce and Moin as a model for turbulent non-premixed combustion in large-eddy simulation. The filtered chemical source term in this model appears in unclosed form, and is modeled by a presumed probability density function (PDF) for the joint PDF of the mixture fraction Z and a flamelet parameter λ. While the marginal PDF of Z can be reasonably approximated by a beta distribution, a model for the conditional PDF of the flamelet parameter needs to be developed. Further, the ability of FPVA to predict extinction and re-ignition has also not been assessed. In this paper, we address these aspects of the model using the DNS database of Sripakagorn et al. It is first shown that the steady flamelet assumption in the context of FPVA leads to good predictions even for high levels of local extinction. Three different models for the conditional PDF of the flamelet parameter are tested in an a priori sense. Results obtained using a delta function to model the conditional PDF of λ lead to an overprediction of the mean temperature, even with only moderate extinction levels. It is shown that if the conditional PDF of λ is modeled by a beta distribution conditioned on Z, then FPVA can predict extinction and re-ignition effects, and good agreement between the model and DNS data for the mean temperature is observed.     

Experimental and numerical investigation of hetero-/homogeneous combustion of CO/H2/O2/N2 mixtures over platinum at pressures up to 5 bar

The hetero-/homogeneous combustion of fuel-lean CO/H₂/O₂/N₂ mixtures over platinum is investigated at pressures up to 5 bar, inlet temperatures (T_IN) up to 874 K, and a constant CO:H₂ molar ratio of 2:1. Experiments are performed in an optically accessible channel-flow catalytic reactor and involve planar laser induced fluorescence (LIF) of the OH radical for the assessment of homogeneous (gas-phase) ignition and 1-D Raman measurements of major gas-phase species concentrations over the catalyst boundary layer for the evaluation of the heterogeneous (catalytic) processes. Simulations are carried out with an elliptic 2-D model that includes detailed heterogeneous and homogeneous chemical reaction schemes. The predictions reproduce the Raman-measured catalytic CO and H₂ consumption, and it is further shown that for wall temperatures in the range 975 ? T_w ? 1165 K the heterogeneous pathways of CO and H₂ are largely decoupled. However, for wall temperatures below a limiting value of 710–720 K and for the range of pressures and mixture preheats investigated, CO(s) blockage of the surface inhibits the catalytic conversion of both fuel components. The homogeneous ignition distance is well-reproduced by the model for T_IN > 426 K, but it is modestly overpredicted at lower T_IN. Possible reasons for these modest differences can be the values of third body efficiencies in the gas-phase reaction mechanism. The sensitivity of homogeneous ignition distance on the catalytic reactions is weak, while the H₂/O₂ subset of the CO/H₂/O₂ gaseous reaction mechanism controls the onset of homogeneous ignition. Pure hydrogen hetero-/homogeneous combustion results in flames established very close to the catalytic walls. However, in the presence of CO the gaseous combustion of hydrogen extends well-inside the channel core, thus allowing homogeneous consumption of H₂ at considerably shorter reactor lengths. Finally, implications of the above findings for the design of syngas-based catalytic reactors for power generation systems are discussed.     

Simultaneous imaging of Mie scattering,PAHs laser induced fluorescence and soot laser induced incandescence to a lab-scale turbulent jet pulverized coal flame

In this study, transient soot formation processes in a small-scale jet burner (CRIEPI burner) were investigated by simultaneous measurements of coal particles, polycyclic aromatic hydrocarbons (PAHs) and soot. Pairs of simultaneous measurements of “Mie scattering measurement for coal particles with laser induced fluorescence (LIF) for PAHs” and “LIF for PAHs with laser induced incandescence (LII) for soot” were performed to understand the transitive formation processes of soot formation in pulverized coal flame, whose signals were successfully separated. Findings in the present study are as follows. Coal particles, PAHs and soot were distributed in this order in radial direction from the central axis. Existing regions of coal particles, PAHs and soot were overlapped from the time averaged viewpoint while there were few overlapping areas of coal particles, PAHs and soot from the instantaneous viewpoint. This result indicates that a long time is required for the formation of soot from 2 to 3 rings PAHs through larger PAHs.     

Nitrogen release during reaction of coal char with O2, CO2, and H2O

The chemistry of char-N release and conversion to nitrogen-containing products has been probed by studying its release and reactions with O₂, CO₂, and H₂O. The experiments were performed in a fixed bed flow reactor at pressures of up to 1.0 MPa. The results show that the major nitrogen-containing products observed depend on the reactant gas; with O₂, NO, and N₂ being the major species observed. Char-N reaction with CO₂ produces N₂ with very high selectivity over a broad range of pressures and CO₂ concentrations, and reaction with H₂O gives rise to HCN, NH₃, and N₂. Observed distributions of nitrogen-containing products are little affected by pressure when O₂ and CO₂ are the reactant gases, but increasing pressures in the reaction with H₂O results in the formation of increasing proportions of NH₃. Formation of NH₃ is also promoted by increasing concentrations of H₂O in the feed gas. The results suggest that NO and HCN are primary products when O₂ and H₂O, respectively, are used as the reactant gases, and that the other observed products arise from interactions of these primary products with the char surface.     

Nonsequential double ionization of nonaligned diatomic molecules N2 and O2

Nonsequential double ionization (NSDI) processes of nonaligned diatomic molecules N2 and O2 are studied using the S-matrix theory. Our results show that the NSDI process significantly depends on the molecular symmetry and structure. The ratio of NSDI rate to single ionization rate as a function of the field intensity is obtained. It is found that N2 behaves closely with its companion atom Ar in the ratios over the entire intensity range, while O2 exhibits an obvious suppression effect, which is qualitatively consistent with the experiment.