Flame Surface Area Research Articles

Equivalence Ratio-Driven Flame Response of an Industrial Premixed Burner: Experiments and Modeling

Abstract This paper presents an analysis of the unsteady heat release rate response of premixed flames to equivalence ratio perturbations for an industrial premixed swirl-based burner. During this investigation, perfectly and technically premixed flames were acoustically forced via fuel/air mixture flow and air flow modulations, respectively, at the same operating conditions. From the resulting flame transfer functions (FTFs), measured using the multimicrophone method, the equivalence ratio driven FTF was isolated and extracted by removing the velocity driven component, i.e., the measured FTF from the perfectly premixed flame, from the technically premixed FTF with two novel extraction techniques. The results are compared with FTFs obtained directly in a previous experimental campaign where the fuel flow was acoustically forced, the resulting equivalence ratio fluctuations measured via an infrared absorption technique, and the heat release rate response to the forcing was quantified using chemiluminescence measurements. The results from both measurement approaches agreed well highlighting the validity of the techniques. Further, to understand the governing features of the equivalence ratio driven FTF, a physics-based analytical model following the G-equation approach was developed. The contributions from flame surface area, flame speed, and heat of reaction oscillations were modeled to describe the heat release rate dynamics. A limited number of physical parameters in the analytical model were anchored on one test condition, optimized and restricted to values, which were all physically reasonable, and were subsequently used for model predictions at other operating conditions. The FTF model predictions compared well with experimental data across a range of different operating conditions. Finally, the relative contributions from flame surface area, flame speed, and heat of reaction oscillations on the features of the FTFs were identified and explored.

Scaling transition of turbulent flame speed for thermodiffusively unstable flames

This work presents an experimental set of Bunsen flames characterized by a moderate Reynolds number and a variable turbulence intensity. Ten lean hydrogen-enriched methane–air mixtures at three turbulence levels are investigated, ranging from methane–air to hydrogen–air mixtures. Such mixtures are selected to have an almost constant laminar flame speed while inducing the onset of thermal-diffusive (TD) instability by gradually increasing the hydrogen content of the blend. The flames' global consumption speed, stretch factor, and flame surface area are investigated and discussed as functions of the effective Lewis number of the mixture. As the interplay between TD instability and turbulence enhances the overall flame propagation, below a transitional Lewis number, flames are observed to be particularly sensitive to external turbulent forcing. This synergistic interaction is discussed in terms of Karlovitz and Lewis numbers. A parameterization of the turbulent flame speed is thus proposed, based on a functional form depending, concurrently, on both Karlovitz and Lewis numbers. The proposed form is shown to fit the experimental results at different turbulence levels and to capture the flame speed enhancement across the transitional Lewis number.

Three-dimensional dynamics of unstable lean premixed hydrogen-air flames: Intrinsic instabilities and morphological characteristics

The 3D swinging laser sheet technique was employed to study the development and morphological characteristics of premixed hydrogen-air unstable flames in a spherical explosion vessel. Pressure dependencies for laminar flame propagation were sought to exploit the role of the Darrieus-Landau (DL) and Thermal-diffusive (TD) instabilities in the unstable self-accelerating flame regime. A sufficiently low Markstein number, as a consequence of the increased pressure, leads to more cracking and smaller cells over the flame surface. The degree of wrinkling on the flame surface is proportional to the increase in flame burning velocity, a relationship that holds true for low pressures but is not applicable under high pressures. External turbulence can significantly alter the extent of flame surface wrinkling even at low root mean square velocities, producing a more wrinkled flame surface compared to intrinsic cellularity, and distinctly affecting flame dynamics. The increased wrinkling and flame speed due to external turbulence can be attributed to the synergistic effects between thermo-diffusive instabilities and turbulence, resulting in higher fuel consumption rates per flame surface area and the formation of finger-like structures that enhance flame displacement speed in curved segments. The parameters, ϵ, deviation of the Lewis number from a critical value, and ω2, obtained through classical linear stability analysis, display a clear linear relationship with the ratio of the wrinkled surface area observed in planar flames. This study enhances the understanding of hydrogen flame instabilities, which is crucial for preventing explosions in hydrogen storage and utilization, and provides valuable insights into flame dynamics, supporting the design of safer and more efficient hydrogen-fueled engines and turbines.

DNS of laboratory-scale turbulent premixed counterflow flames under elevated gravity conditions

In the present work, direct numerical simulations (DNSs) of laboratory-scale turbulent premixed counterflow flames are performed to understand the effect of elevated gravity on flame structures. In the DNS, a turbulent jet of CH4/N2/O2 mixture is injected in opposition to a stream of combustion product. Three cases were considered, i.e., one case with normal gravity and two cases with elevated gravity in different directions. The DNS results of the case under normal gravity conditions were found to be in good agreement with the experimental data. The turbulent flame speed is the highest in the case with elevated gravity aligned with the mean flame propagating direction, which is due to the increase in flame surface area and stretch factor. The gravity levels have a substantial influence on the flame extinction characteristics. It was observed that the flame front is pushed toward the direction of the gravity. When the elevated gravity antialigns with the mean flame propagating direction, the flame front travels across the stagnation plane toward the combustion product stream, leading to significant reactant dilution and flame extinction. The flame normal vector preferentially aligns with the most compressive strain rate for various gravity conditions. The tangential strain rate is the greatest (smallest) when elevated gravity antialigns (aligns) with the mean flame propagating direction. The impact of elevated gravity on the tangential strain rate at the flame fronts ultimately leads to a displacement in the flame front location relative to the stagnation plane, thereby influencing flame extinction.

Numerical study on flame acceleration and deflagration-to-detonation transition affected by the solid obstacles with different shapes

To investigate which shape of solid obstacles is most optimal in terms of detonation-assisting capability within a closed channel, this study utilizes the Large Eddy Simulation (LES) method, based on the OpenFOAM platform, to conduct a detailed numerical study on the effects of different shaped obstacles (rectangular, rhombus, trapezoidal, circular, triangular) on flame acceleration and detonation initiation processes. The results show that the obstacles arranged in the middle of the combustion chamber can create more passages and local flames. Meanwhile, changes in the shape of obstacles can also produce different vortex structures and recirculation zones. These differences affect the flame surface area and the total combustion heat release rate, thereby affecting the flame acceleration effect. Judging from the results, triangular obstacles have the best flame acceleration effect, followed by rhombus, trapezoid, rectangle, and circular obstacles. In terms of detonation initiation, the obstacles with a slope can induce vortex to detach earlier, which produce a better vortex-flame interaction. Simultaneously, the slope can reflect the shock wave more frequently, creating more favourable conditions for the formation of Mach reflection and Mach stem. Based on this, it is found that trapezoidal and triangular obstacles have better detonation-assisting capability. Additionally, it is found that the stronger the leading shock wave and the shorter the distance between it and the flame front, the better the detonation initiation effect. In general, the types of detonation initiation in this study all belong to the shock detonation transition (SDT). However, the detonation initiation process can be further classified into two categories: (I) Detonation induced by shock wave reflection; (II) Detonation triggered by shock wave focusing.

Numerical study on flame acceleration and deflagration-to-detonation transition: Spatial distribution of solid obstacles

This study conducts a detailed numerical investigation on the spatial distribution of solid obstacles using the large eddy simulation method. It is discovered that although flame acceleration induced by solid obstacles is dominated by factors such as flow field disturbances, vortices and recirculation zones, turbulence, flame surface areas and combustion heat release rates, etc., the characteristics of the leading shock wave are key to detonation initiation. Specifically, the intensity of the leading shock wave, its formation time, and its distance from the flame front significantly affect detonation initiation. Depending on the state of the shock wave, the detonation initiation process may occur through various mechanisms such as shock reflection, shock focusing. Overall, the types of detonation initiation in this study all belong to the shock detonation transition. However, the detonation initiation process can be further classified into two categories: (I) Detonation induced by shock wave reflection; (II) detonation triggered by shock wave focusing. Despite certain disparities in the detonation initiation process, all detonation initiation processes conform to the gradient theory, and the flame evolution processes in all cases consistently follow three stages: the laminar slow-ignition stage; the turbulent deflagration stage; the detonation initiation stage. Furthermore, the study further discerns that, compared to positioning obstacles on the wall, placing obstacles inside the combustion chamber can further augment the detonation-assisting effect. However, excessively sparse or dense spatial distributions of solid obstacles fail to yield the optimal detonation effect. An optimal distribution exists, which triggers the fastest detonation initiation.

Development of a multiple laser-sheet imaging technique for the analysis of three-dimensional turbulent explosion flame structures

The present work describes the development of a multiple laser-sheet imaging technique for the study of turbulent, premixed flame surface structures at high Karlovitz stretch factor values. Experiments were conducted using CH4 air mixtures at 365 K and 0.5 MPa and up to a root mean square turbulence velocity of 1.5 m/s. A high-speed Nd:YAG laser capable of pulsing up to 60 kHz in conjunction with a high-speed camera and a rotating mirror was used to reconstruct time-resolved three-dimensional turbulent flames. This has, for the first time, enabled the direct measurement of wrinkled flame surface areas, along with reaction progress variable and flame brush thickness. These are important parameters for the characterization of turbulent burning rates and provide more insight into the dynamic nature of the flames and their structures. In addition, the current data aid toward direct comparison with results from combustion simulation studies.

A Comparison Between Conventional and MILD Combustion Processes of Turbulent Stratified Mixtures

ABSTRACT The relative enhancements of volume-integrated burning rate and flame surface area under turbulence and heat release rate contributions due to different combustion and flame propagation modes have been compared between conventional and Moderate or Intense Low Oxygen Dilution (MILD) combustion of stratified mixtures (for a global mean equivalence ratio of 0.8) for the first time using three-dimensional Direct Numerical Simulation (DNS) data. The burning rate enhancement under turbulence is found to be comparable to the enhancement of flame surface area for the conventional stratified flames considered here, whereas the burning rate enhancement is found to be considerably greater than flame area augmentation in the MILD cases analyzed here. It has been found that almost all the heat release rate arises due to the lean-premixed mode of combustion for the conventional stratified flame cases considered here. In contrast, the lean premixed mode of combustion accounts for about 50% of the total heat release rate for the MILD cases where rich-premixed and non-premixed modes of combustion contribute significantly to the overall heat release rate. Under the conditions analyzed here, more than 75% of heat release rate arises from the propagation-dominated regions for both MILD combustion and stratified flame cases, but the trend is stronger in the conventional stratified mixture combustion. More than 80% of heat release rate in the globally lean conventional stratified flame cases considered here originates from back-supported flame propagation, whereas both front and back-supported flame propagation modes have comparable contributions to heat release in the MILD combustion cases. The non-premixed mode of combustion has been found to play a more important role in the MILD combustion cases than in the case of conventional stratified flames. This implies that the closures, which work reasonably well for MILD combustion of homogeneous mixtures, may not be sufficient for modeling of MILD combustion of stratified/inhomogeneous mixtures.

Impact of H[formula omitted]-enrichment on the response of a partially premixed CH[formula omitted]–air flame to velocity and equivalence ratio fluctuations

The dynamics of hydrogen-enriched partially premixed methane–air flames are investigated, focusing on the frequency response to velocity and equivalence ratio fluctuations. Large-eddy simulations (LES) with conjugate heat transfer are performed on the partially premixed PRECCINSTA combustor for various methane–hydrogen mixtures, maintaining a constant thermal power and global equivalence ratio of 10 kW and 0.8, respectively. To reduce the computational cost of the LES, a three-step global mechanism is calibrated for each hydrogen-enriched partially premixed flame. The accuracy of the LES is assessed against experimental data, showing a good agreement with the mean velocity field and shape of the swirl-stabilized flame. The multiple-input–single-output structure for the flame transfer function is applied to separate the individual contributions of velocity and equivalence ratio fluctuations to the overall heat release rate oscillation in the partially premixed flames. System identification with broadband forcing is used to compute the various flame transfer functions. Hydrogen enrichment significantly influences the flame response to equivalence ratio fluctuations, explained by the respective impact on the response of the heat of reaction, the flame speed, and the flame surface area to equivalence ratio fluctuations. The interaction between velocity and equivalence ratio fluctuations is further characterized by the mixing transfer function, which describes the equivalence ratio fluctuations reaching the flame. The higher mass diffusivity of hydrogen results in lower equivalence ratio fluctuations reaching the flame as the frequency increases, attenuating the pronounced effect of hydrogen on the flame sensitivity to equivalence ratio fluctuations. This study contributes to the understanding of the impact of hydrogen enrichment on the thermoacoustic response of partially premixed flames.

The stability and morphology of laminar premixed hydrogen/air flames under various gravity conditions

In the present work, direct numerical simulations (DNS) were employed to understand the stability and morphology of laminar lean premixed hydrogen/air flames under various gravity conditions. The combined effect of Darrieus–Landau (DL) instability, thermal-diffusive (TD) instability and Rayleigh–Taylor (RT) instability was studied quantitatively for both the initial linear regime and the long-term nonlinear regime, with a particular focus on the effect of RT instability. In the linear regime, it was found that the RT instability has a negligible effect on the linear growth rate of flames with an equivalence ratio of 0.4, but the effect is significant when the equivalence ratio decreases to 0.3. The performance of various models for the dispersion relation was discussed, including a theoretical model considering the effects of all instabilities, which captures well the trend of the simulated results in the cases with unity Lewis numbers. The hydrogen/air flames with an equivalence ratio of 0.3 was investigated comprehensively in the nonlinear regime. For the cases with unity Lewis numbers, the RT-unstable flame exhibits a similar single-cusp structure as the RT-neutral flame, but with a deeper cusp and larger consumption speed as a result of the destabilizing effect of RT instability. The RT-stable flame gradually returns to a planar flame shape owing to the stabilizing effect of RT instability. For the cases with non-unity Lewis numbers, the effect of domain size on the flame evolution was investigated, showing a domain width of 200δT is enough to obtain a domain-independent result, where δT is the laminar flame thickness. It was found that the cases initialized with multimode perturbations reach the quasi-steady state faster, while the flame statistics are consistent with those of the cases initialized with single mode perturbations. The effect of RT instability was explored by comparing the results from various gravity conditions. The flame consumption speed is the largest in RT-unstable flames, and the smallest in RT-stable flames. The variation in flame consumption speed is due to the difference in the flame surface area, while the local reactivity and flame structure are not influenced by the RT instability. Further fractal analysis reveals that although the RT instability causes significant differences in flame surface area, the characteristic fractal dimension remains constant for flames under different gravity conditions.Novelty and significanceIn the present work, direct numerical simulations were employed to understand the stability and morphology of RT-unstable, RT-neutral and RT-stable laminar lean premixed hydrogen/air flames. For the first time, the effect of DL, TD and RT instabilities was quantitatively examined for both the linear regime and the nonlinear regime. A theoretical model considering various instabilities was compared with the simulated dispersion relations in the linear regime. In the non-linear regime, the flame surface area and corresponding flame consumption speed are influenced by the RT instability, while the flame stretch factor remains unaffected. A fractal analysis was conducted in laminar premixed hydrogen/air flames under the influence of various instabilities. Remarkably, it was discovered that the characteristic fractal dimension remains constant for the cases with and without the RT instability.

Choice of reaction progress variable under preferential diffusion effects in turbulent syngas combustion based on detailed chemistry direct numerical simulations

The combustion of hydrogen and carbon-monoxide mixtures, so-called syngas, plays an increasingly important role in the safety context of non-fossil energy generation, more specifically in the risk management of incidents in process engineering plants for ammonia synthesis and in nuclear power plants. In order to characterize and simulate syngas/air combustion on industrially relevant scales, subgrid modelling is required, which is often based on a reaction progress variable. To understand the influence of different fuel compositions, turbulence intensities and flame topologies on different possible definitions of reaction progress variable, detailed chemistry direct numerical simulations data of premixed, lean hydrogen/air and syngas/air flames has been considered. A reaction progress variable based on normalized molecular oxygen mass fraction has been found not to capture the augmentation of the normalized burning rate per unit flame surface area in comparison to the corresponding 1D unstretched premixed flame due to preferential diffusion effects. By contrast, reaction progress variables based on other individual species, such as hydrogen, can capture the augmentation of the rate of burning well, but exhibit a pronounced sensitivity to preferential diffusion effects, especially in response to flame curvatures. However, a reaction progress variable based on the linear combination of the main products can accurately represent the temperature evolution of the flame for different mixtures, turbulence intensities and varying local flame topology, while effectively capturing the augmentation of burning rate due to preferential diffusion effects. However, its tendency to assume values larger than 1.0 in the regions of super-adiabatic temperatures poses challenges for future modeling approaches, whereas the reaction progress variable based on hydrogen mass fraction remains bound between 0.0 and 1.0 despite showing deviations in comparison to corresponding variations obtained from the unstretched laminar flame depending on flame curvature variations.

Large eddy simulation investigation of the effect of radiative heat transfer on the ignition progress in a model combustor

Radiation significantly influences turbulent combustion. This paper investigates the impact of radiative heat transfer on the ignition process in a multi-staged model combustor. The ignition process is simulated both without and with consideration of radiation. Large Eddy Simulations based on the Euler-Lagrange description of the liquid spray and the Dynamic Thickened Flame model coupled with a skeletal chemical reaction mechanism are employed. The Weighted Sum of Gray Gases method are employed to solve the spectral properties and the Discrete Ordinates method is employed to solve the radiative transfer equation. The results indicate that radiative heat transfer influences the flame structure during the ignition process. In the flame propagation phase, both the total heat transfer and the exit average temperature are higher, with the radiative heat transfer during ignition being less than 15% of the total heat transfer. During the kernel generation phase, radiative heat transfer leads to a 7.24% increase in the maximum temperature of the flame kernel, a 7.74% expansion of the flame surface area, and outward movement of the position of the maximum heat release rate. Additionally, the ignition delay time extends by 14.91%. In the flame propagation phase, the flame transfers heat to the mixed gas and droplets near the flame surface through radiation, promoting flame spread.

Implications of Using Scalar Forcing to Sustain Reactant Mixture Stratification in Direct Numerical Simulations of Turbulent Combustion

A recently proposed scalar forcing scheme that maintains the mixture fraction mean, root-mean-square and probability density function in the unburned gas can lead to a statistically quasi-stationary state in direct numerical simulations of turbulent stratified combustion when combined with velocity forcing. Scalar forcing alongside turbulence forcing leads to greater values of turbulent burning velocity and flame surface area in comparison to unforced simulations for globally fuel-lean mixtures. The sustained unburned gas mixture inhomogeneity changes the percentage shares of back- and front-supported flame elements in comparison to unforced simulations, and this effect is particularly apparent for high turbulence intensities. Scalar forcing does not significantly affect the heat release rates due to different modes of combustion and the micro-mixing rate within the flame characterised by scalar dissipation rate of the reaction progress variable. Thus, scalar forcing has a significant potential for enabling detailed parametric studies as well as providing well-converged time-averaged statistics for stratified-mixture combustion using Direct Numerical Simulations in canonical configurations.

Investigation of differential diffusion in lean, premixed, hydrogen-enriched swirl flames

Hydrogen combustion is a promising alternative to fossil fuel combustion in an effort to reduce our carbon footprint. However, hydrogen combustion is prone to thermodiffusive instabilities largely dependent on differential diffusion, a phenomenon that can lead to higher probabilities of flashback in industrial burners, given hydrogen’s high reactivity and diffusivity. This paper evaluates low-swirl flames of methane and air enriched with hydrogen to highlight the onset of differential diffusion. Testing was conducted in a fully controllable swirl burner, where bulk velocity Uav = 13 m/s and swirl number S = 0.6 were kept constant for each hydrogen–methane blend to isolate increases in flame surface area from increases in turbulence intensity. Furthermore, each fuel blend of hydrogen and methane is evaluated at the same laminar flame speed of SLo = 0.267 m/s to isolate flame stretch effects on the turbulent burning rate. Combined hydroxyl (OH) PLIF and stereoscopic PIV at the National Research Council of Canada were used to analyze the OH fluorescence in a 2D-3C velocity field for each flame condition. High-speed PIV at McGill University was used to resolve local flame phenomena, such as local flame displacement velocity and flame stretch rate. Using these techniques, it can be observed that the flame displaces axially in response to turbulent flame speed while exhibiting increases in flamefront wrinkling. This increased corrugation due to flame stretch is highlighted in the PDFs of local curvature and κSf and is further evidenced by a shift towards positive curvatures (κ> 0) for increasing H2 volume fraction. This trend suggests that there is a strong correlation with increases in turbulent burning rate and positive curvature as a result of differential diffusion, but it is not necessarily a control mechanism of the most forward propagating points proposed by the theory of leading points.

Study on the flame radiative heat transfer and near-field radiation heat flux predictive model of vehicle fires in a tunnel

As a major form of heat transfer, flame thermal radiation has an important effect on the ignition of adjacent energy, fire spread and environmental damage. Forecasting the flame radiation heat flux can provide theoretical basis for fire hazard assessment, fire extinguishing and salvation. Meanwhile, the distinctive structure of the arc tunnel can result in significant variations in the local heat and mass transfer processes during tunnel fires compared to those in rectangular tunnels, leading to different fire characteristics parameters. In this study, experiments and theoretical analyses were conducted to investigate the flame length scale and near-field radiation heat flux on the lower surface of an arc-shaped ceiling induced by vehicle fires in a tunnel. The findings reveal that the flame length scale in the transverse direction beneath the arc-shaped ceiling decreases with an increasing burner aspect ratio and increases with the increasing heat release rate. A dimensionless model of the flame length scale under the arced ceiling with various burner aspect ratios was proposed based on the amount of remaining unburned fuel after flame impingement and the buoyancy component along the arc-shaped ceiling. Simultaneously, radiative heat fluxes on the lower surface also increase with rising heat release rate. The flame geometry beneath the arc-shaped ceiling is assumed to be a combination of a cuboid and an arc plate model based on the actual flame length scale. Subsequently, the view factor, radiative fraction, and flame surface area are analyzed to predict the radiative heat fluxes on the lower surface. The predictions of the cuboid and arc plate flame model show good agreement with the current experimental data. This work presents a dimensionless correlation of the flame length scale and calculates the view factor between the target and the arc plate to predict the near-field radiative heat flux more accurately. This approach also facilitates the understanding of flame spreading behavior and thermal hazards impacted by arced ceiling.

A Novel Approach for Accurate Determination of Laminar Burning Velocity Using the Constant Volume Method in a Cylindrical Vessel with Full Optical Access

ABSTRACT The constant volume method (CVM) is a potential method for determining laminar burning velocity (LBV) at high pressures and temperatures. However, CVM using spherical vessels gives no optical access to the flame; in cylindrical vessels, with full optical access, the flame does not remain spherical. In either case, one has to rely on semi-analytical models, which have limitations in accuracy, to obtain burned gas volume (BGV) needed for LBV calculation. This study presents a novel approach for accurate determination of BGV and flame surface area (FSA) of non-spherical flames in cylindrical vessels from the measured end-view flame radius. Through CFD simulation, normalized relationships are developed to relate the measured flame end-view radius to the actual FSA and BGV. It is shown that these relationships are purely geometry-dependent and are independent of vessel size and combustion reactions in the vessel. This facilitates accurate determination of LBV using CVM in an optically accessible cylindrical vessel, eliminating the requirement to estimate burned gas mass fraction using burned gas mass fraction models. Experiments using CVM were performed for a stoichiometric methane-air mixture at 1.0 bar and 300 and 350 K. This new approach is demonstrated by using it to determine the LBV of stoichiometric methane-air mixture accurately. The present methodology is expected to provide a way to calculate accurate LBV without the need for additional thermodynamic models for burned gas mass fraction.

A novel Swiss-roll counterflow micro-combustor: Experimental investigation of flame dynamic characteristics by spectroscopy and RGB image processing methods

In this study, a novel Swiss roll counterflow micro-combustor is developed to improve flame stability, preheating, and heat recovery. The study experimentally investigates the combustion characteristics of methane-oxygen mixtures in a constant-temperature aluminum chamber over a range of equivalence ratios using spectroscopic methods and RGB image processing. The results indicate that this micro-combustion flame is a cool flame with an inverse relationship between variations of flame thickness and flame length. Furthermore, the study explores the behavior of some radicals, including OH*, CH*, C2* , H2O*, and CO2* and their radiation intensity, which are influenced by the equivalence ratio and mass flow rate of fuel to the oxidizer, as well as other factors such as flame surface area and heat transfer from the wall. RGB image processing is also utilized to investigate the changes of CH* and C2* in channels G and B, respectively. The study shows a strong correlation between the visible light intensity of the flame and the R channel. The results demonstrate that for the equivalence ratios Between 1.93 and 2.26, the radiation intensity of the radicals is relatively stable, which makes it applicable for micro-combustors.

Effect of hydrogen blending on ammonia/air explosion characteristics under wide equivalence ratio

In order to evaluate the effect of mixing hydrogen on the explosive characteristics of ammonia, this work explored the evolution of explosion overpressure and flame propagation characteristics of ammonia/hydrogen/air mixture under different hydrogen blending ratios in closed pipelines through experimental research. Experimental results show that as the hydrogen blending ratio increase, the instability of flame increases, and flame surface changes from smooth to wrinkled. Flame propagation undergoes the evolution of jellyfish-shaped/near-spherical, finger-shaped, plane-shaped, irregularly protruding flame and tulip flame. Flame surface area is constantly changing, accompanied by violent flame tip velocity pulsations. The maximum explosion pressure and the maximum explosion pressure rise rate are consistent with the evolution law of adiabatic flame temperature. They first increase and then decrease with increase of the equivalence ratio, and their peak value appears near the equivalence ratio of 1.1. On the contrary, the explosion duration and rapid combustion period first decrease and then increase with increase of equivalence ratio, and decrease monotonically with increase of hydrogen blending ratio. The conclusions obtained in this work will contribute to the structural design of flame arresters and the formulation of explosion venting strategies during the preparation, transportation and utilization of ammonia/hydrogen mixtures.

Effects of dilatation and turbulence on tangential strain rates in premixed hydrogen and iso-octane flames

The tangential strain rate in premixed flames impacts significantly the flame surface area generation and thus the combustion process. Studies on incompressible isotropic turbulence have revealed that the mean tangential strain rate at material and iso-scalar surfaces is positive and exhibits a universal value when normalized by the Kolmogorov time. This is associated with the preferential alignment of the surface normal with the most compressive principal strain rate. The present study investigates such effects in premixed hydrogen and iso-octane flame kernels using direct numerical simulations. It is shown that the normalized mean tangential strain rate of the investigated flames has a very similar value compared with the incompressible flows. However, in the reaction zone, the flame surface normal aligns preferentially with the most extensive principal strain rate. Furthermore, this alignment depends on the reaction progress variable and the Lewis number, while the tangential strain rate remains independent of these parameters. Such counter-intuitive behaviour is systematically investigated by decomposing the effects of dilatation and residual solenoidal turbulence. It is found that the solenoidal turbulence influences significantly the tangential strain rate. A general effect of turbulence on the tangential strain rate is identified, which is consistent with incompressible flows and independent of the Lewis number and the reaction progress variable. This is a remarkable finding indicating that models of the tangential strain rate developed based on incompressible flows apply also to premixed flames with different Lewis numbers, and, for the modelling, only the solenoidal turbulence should be considered.

Opposite effects of flame-generated potential and solenoidal velocity fluctuations on flame surface area in moderately intense turbulence

To explore the influence of dilatational and rotational motions on iso-scalar surface area within a premixed turbulent flame, three-dimensional compressible direct numerical simulation data obtained by Dave et al. (2018) from a lean, complex-chemistry, hydrogen-air flame propagating in intense small-scale turbulence (a ratio of laminar flame thickness to Kolmogorov length scale is about 20) in a box are analyzed using Helmholtz–Hodge decomposition of velocity field into solenoidal and potential components. The obtained results indicate that the flame surface area is created by potential velocity fluctuations generated due to combustion-induced thermal expansion, whereas the rotational motion acts to smooth wrinkles of iso-scalar surfaces within local preheat and reaction zones and, consequently, to reduce the flame surface area. The latter finding challenges the classical concept of flame-generated turbulence and is attributed to flame-generated solenoidal velocity fluctuations. Specifically, the incoming tangential (to the flame) vorticity is suppressed by baroclinic torque, which also generates vorticity in another tangential direction, with the latter (flame-generated) solenoidal velocity fluctuations working to smooth the flame surface. Thus, even under conditions of moderately intense turbulence, flame surface area can primarily be created by potential velocity perturbations caused by combustion-induced thermal expansion.