Rich n-heptane and diesel combustion in porous media

Rich n-heptane and diesel flames in two-layer porous media are experimentally investigated in the context of syngas production. The stable operating points of n-heptane reforming have been determined and the mole fractions of H₂, CO, CO₂ and light hydrocarbons have been measured in the exhaust gas at an equivalence ratio of 2 for different thermal input values. The reformer performance has been assessed also from the point of view of the heat losses and the mixture homogeneity. The pre-vapouriser produces an approximately uniform vapour–air mixture upstream of the flame front. The range of flow rates for stable flames decreased with increasing equivalence ratio. Heat losses were about 10% of the thermal input at high firing rates. A 77.2% of the equilibrium H₂ was achieved at a flame speed of 0.82 m/s. The same reactor with a different porous matrix for the reforming stage demonstrates diesel reforming to syngas with a conversion efficiency of 77.3% for a flame speed of 0.65 m/s.     

A LES-CMC formulation for premixed flames including differential diffusion

A finite volume large eddy simulation–conditional moment closure (LES-CMC) numerical framework for premixed combustion developed in a previous studyhas been extended to account for differential diffusion. The non-unity Lewis number CMC transport equation has an additional convective term in sample space proportional to the conditional diffusion of the progress variable, that in turn accounts for diffusion normal to the flame front and curvature-induced effects. Planar laminar simulations are first performed using a spatially homogeneous non-unity Lewis number CMC formulation and validated against physical-space fully resolved reference solutions. The same CMC formulation is subsequently used to numerically investigate the effects of curvature for laminar flames having different effective Lewis numbers: a lean methane–air flame with Le_eff = 0.99 and a lean hydrogen–air flame with Le_eff = 0.33. Results suggest that curvature does not affect the conditional heat release if the effective Lewis number tends to unity, so that curvature-induced transport may be neglected. Finally, the effect of turbulence on the flame structure is qualitatively analysed using LES-CMC simulations with and without differential diffusion for a turbulent premixed bluff body methane–air flame exhibiting local extinction behaviour. Overall, both the unity and the non-unity computations predict the characteristic M-shaped flame observed experimentally, although some minor differences are identified. The findings suggest that for the high Karlovitz number (from 1 to 10) flame considered, turbulent mixing within the flame weakens the differential transport contribution by reducing the conditional scalar dissipation rate and accordingly the conditional diffusion of the progress variable.     

The development and application of a fast-response DC-coupled flame ionization current sensor

The principles of the design of a fast-response DC-coupled ionisation probe and associated electronic circuit are described for the measurement of the mean and fluctuating components of the ion concentration in the reaction zones of stationary and transient flames. The performance of the ionisation probe is evaluated by bench tests and with measurements of ion current mean and fluctuating levels in laminar and turbulent premixed flames and in the combustion chamber of a spark ignition engine. The results demonstrate the wide frequency bandwidth, good spatial resolution and electronic noise immunity of the ionisation probe, making it suitable for a host of combustion applications.     

Mesoscale strain measurement in deformed crystals: A comparison of X-ray microdiffraction with electron backscatter diffraction

Mapping of residual stresses at the mesoscale is increasingly practical thanks to technological developments in electron backscatter diffraction (EBSD) and X-ray microdiffraction using high brilliance synchrotron sources. An analysis is presented of a Cu single crystal deformed in compression to about 10% macroscopic strain. Local orientation measurements were made on sectioned and polished specimens using EBSD and X-ray microdiffraction. In broad strokes, the results are similar to each other with orientations being observed that are on the order of 5° misoriented from that of the original crystallite. At the fine scale it is apparent that the X-ray technique can distinguish features in the structure that are much finer in detail than those observed using EBSD even though the spatial resolution of EBSD is superior to that of X-ray diffraction by approximately two orders of magnitude. The results are explained by the sensitivity of the EBSD technique to the specimen surface condition. Dislocation dynamics simulations show that there is a relaxation of the dislocation structure near the free surface of the specimen that extends approximately 650 Å into the specimen. The high spatial resolution of the EBSD technique is detrimental in this respect as the information volume extends only 200 Å or so into the specimen. The X-rays probe a volume on the order of 2 µm in diameter, thus measuring the structure that is relatively unaffected by the near-surface relaxation.     

Direct numerical simulations of autoignition in turbulent two-phase flows

Three-dimensional direct numerical simulations (DNS) were carried out to investigate the impact of evaporation of droplets on the autoignition process under decaying turbulence. The droplets were taken as point sources and were tracked in a Lagrangian manner. Three cases with the same initial equivalence ratio but different initial droplet size were simulated and the focus was to examine the influence of the droplet evaporation process on the location of autoignition. It was found that an increase in the initial droplet size results in an increase in the autoignition time, that highest reaction rates always occur at a specific mixture fraction ξ_MR, as in purely gaseous flows, and that changes in the initial droplet size did not affect the value of ξ_MR. The conditional correlation coefficient between scalar dissipation rate and reaction rates was only mildly negative, contrary to the strongly negative values for purely gaseous autoigniting flows, possibly due to the continuous generation of mixture fraction by the droplet evaporation process that randomizes both the mixture fraction and the scalar dissipation fields.     

A-Priori Validation of Scalar Dissipation Rate Models for Turbulent Non-Premixed Flames

The modelling of scalar dissipation rate in conditional methods for large-eddy simulations is investigated based on a priori direct numerical simulation analysis using a dataset representing an igniting non-premixed planar jet flame. The main objective is to provide a comprehensive assessment of models typically used for large-eddy simulations of non-premixed turbulent flames with the Conditional Moment Closure combustion model. The linear relaxation model gives a good estimate of the Favre-filtered scalar dissipation rate throughout the ignition with a value of the related constant close to the one deduced from theoretical arguments. Such value of the constant is one order of magnitude higher than typical values used in Reynolds-averaged approaches. The amplitude mapping closure model provides a satisfactory estimate of the conditionally filtered scalar dissipation rate even in flows characterised by shear driven turbulence and strong density variation.

     

Experimental Investigation of the Effects of Turbulence and Mixing on Autoignition Chemistry

The autoignition of acetylene, released from a finite-sized circular nozzle into a turbulent coflow of hot air confined in a pipe, has been the subject of a recent experimental study to supplement previous work for hydrogen and n-heptane. As with hydrogen and n-heptane, autoignition appears in the form of well-defined localized spots. Quantitative information is presented concerning the effects of turbulence intensity, turbulent lengthscale and injector diameter on the location of autoignition. The effects of these parameters on inhomogeneous autoignition have not been investigated experimentally before. The present study establishes that increasing the bulk velocity increases the autoignition length, as was reported for hydrogen and n-heptane. For the same turbulence intensity, the autoignition length increases as the injector diameter increases and as the turbulent lengthscale decreases. A simultaneous decrease in turbulence intensity and increase in lengthscale causes a reduction in autoignition length. Further, the frequency of appearance of the autoignition spots has also been measured. It is found to increase when autoignition occurs closer to the injector, and also at higher velocities. The observed trends are consistent with expectations arising from the dependence of the mixture fraction and the scalar dissipation rate on the geometrical and flow parameters. The data can be used for the validation of turbulent combustion models.     

H2/air autoignition: The nature and interaction of the developing explosive modes

Appropriate algorithmic tools are employed for the analysis of the explosive modes developing during the autoignition of homogeneous mixtures. The ability of these tools to provide significant physical understanding is demonstrated in the case of the homogeneous ignition of a stoichiometric H₂/air mixture, modelled by two different chemical kinetics mechanisms. It is shown that the ignition process evolves in two stages. The first stage is characterised by the development of two explosive timescales (one fast and one slow), that lead the system away from equilibrium. As the end of the first stage is approached, the two explosive timescales converge, they merge and then they disappear. In the second stage only dissipative timescales develop, which drive the system all the way to equilibrium. It is shown that throughout the first stage the fast explosive timescale is generated by chain reactions. The slow explosive timescale is initially generated by an initiation reaction that produces the radicals required for the start-up of the fast mode, while later on it is generated by reactions that are responsible for the heat released. These findings are validated with sensitivity analysis results for the ignition delay time and are employed in order to clarify the discrepancies in the solution provided by the two different chemical kinetics mechanisms considered.     

Pre-Chamber Ignition Mechanism: Simulations of Transient Autoignition in a Mixing Layer Between Reactants and Partially-Burnt Products

The structure of autoignition in a mixing layer between fully-burnt or partially-burnt combustion products from a methane-air flame at ? = 0.85 and a methane-air mixture of a leaner equivalence ratio has been studied with transient diffusion flamelet calculations. This configuration is relevant to scavenged pre-chamber natural-gas engines, where the turbulent jet ejected from the pre-chamber may be quenched or may be composed of fully-burnt products. The degree of reaction in the jet fluid is described by a progress variable c (c = taking values 0.5, 0.8, and 1.0) and the mixing by a mixture fraction ξ (ξ = 1 in the jet fluid and 0 in the CH₄-air mixture to be ignited). At high scalar dissipation rates, N₀, ignition does not occur and a chemically-frozen steady-state condition emerges at long times. At scalar dissipation rates below a critical value, ignition occurs at a time that increases with N₀. The flame reaches the ξ = 0 boundary at a finite time that decreases with N₀. The results help identify overall timescales of the jet-ignition problem and suggest a methodology by which estimates of ignition times in real engines may be made.