Large eddy simulation of premixed hydrogen/methane/air flame propagation in a closed duct

In this paper, large eddy simulation (LES) is performed to investigate the propagation characteristics of premixed hydrogen/methane/air flames in a closed duct. In LES, three stoichiometric hydrogen/methane/air mixtures with hydrogen fractions (volume fractions) of 0, 50% and 100% are used. The numerical results have been verified by comparison with experimental data. All stages of flame propagation that occurred in the experiment are reproduced qualitatively in LES. For fuel/air mixtures with hydrogen fractions of 0 and 50%, only four stages of “tulip” flame formation are observed, but when the hydrogen fraction is 100%, the distorted “tulip” flame appears after flame front inversion. In the acceleration stage, the LES and experimental flame speed and pressure dynamic coincide with each other, except for a hydrogen fraction of 0. After “tulip” flame formation, all LES and experimental flame propagation speeds and pressure dynamics exhibit the same trends for hydrogen fractions of 0 and 100%. However, when the hydrogen fraction is 50%, a slight periodic oscillation appears only in the experiment. In general, the different structures displayed in the flame front during flame propagation can be attributed to the interaction between the flame front, the vortex and the reverse flow formed in the unburned and burned zones.     

Numerical studies of premixed hydrogen/air flames in a small-scale combustion chamber with varied area blockage ratio

The increasing use of hydrogen as a renewable source of energy underlines the need to be able to assess the safety risks involved in the event of an accidental explosion. This paper presents numerical studies for hydrogen/air propagating flames at an equivalence ratio of 0.7 in a laboratory-scale combustion chamber equipped with turbulence generating baffles and a solid square cross section obstruction. The large eddy simulation (LES) modelling technique is used with an in-house computational fluid dynamics (CFD) model for compressible flows to study the flow turbulence and the flame propagation characteristics. The study is carried out using four different baffle arrangements and two different solid obstructions with area blockage ratios of 0.24 and 0.5. Results for the generated peak overpressure and the timing at which it occurs following ignition are considered as the primary safety factors. The time histories of the flame speed and position relative to the ignition source are validated against published experimental data. Good agreement is obtained between numerical results and experimental data which enables further predictions where measurements are limited in the study of vented hydrogen explosions. It was concluded that adding successive baffles and increasing the area blockage ratio escalates the maximum rate at which pressure rises and raises the generated peak explosion overpressure.     

Analysis of flame stabilization mechanism in a hydrogen-fueled reacting wall-jet flame

In this study, the novel conservative representation of chemical explosive mode analysis is augmented to analyze the key flame features in the Burrows-Kurkov flames simulated by both Reynolds-Averaged Navier-Stokes (RANS) and large eddy simulation (LES). Subtle difference are revealed in flame stabilization mechanisms resulting from the difference in modeling and spatial resolution. RANS shows that, ahead of the flame onset location, the composition diffusion and shock wave compression play dominant roles in chemical explosion indicating that the flame is stabilized by the assisted-ignition combustion mode. In contrast, LES shows that the flame is stabilized by the auto-ignition mode since the nonchemical contribution counteracts chemical reaction during the development of ignited flame kernels. For RANS, the radical pool builds up through the unphysical back diffusion near the flame stabilization front, which reveals the limitation of RANS method in the resolution and characterization of the key flame features in Burrows-Kurkov flames.     

Control of methane flame properties by hydrogen fuel addition: Application to power plant combustion chamber

This study presents a numerical investigation of the effects of mixing methane/hydrogen on turbulent combustion processes taking place in a burner similar to that integrated in gas turbine power plants. Thereby, in comparison to the reference case where the burner is fuelled by 100% of methane, the variations of the axial velocity field, temperature field and mass fraction of carbon monoxide field are examined for different percentages of hydrogen fuel injection. The computed results, obtained by using the software Fluent-CFD, are compared and validated against experimental reference data. Results show that the hydrogen addition to the methane has an impact on all physical and chemical parameters of the reactive system.     

Dynamics of premixed hydrogen-air flame propagation in the duct with pellets bed

Hydrogen, as the promising clean alternative energy in the future, is in the spotlight now all over the world. However, its flammable and explosive hazards should be highly considered during its practical application. In this study, the experiments are performed to study premixed hydrogen-air flame propagation in the duct with pellets bed, especially for fuel-rich condition. High-speed schlieren photography is employed to capture flame front development during the experiments. As well as the pressure transducer, is used to track the pressure buildup in the flame propagation process. Different diameters of pellets and different concentrations of gas mixture are considered in this experimental study. The typical evolutions about the tulip flame are similar in all cases, although the tulip flame formation time caused by the laminar flame speed are different. The flame propagation velocity is pretty enhanced in fuel-lean mixture under the effect of large diameter pellets bed, but it is significantly suppressed in fuel-rich conditions. While for the small diameter pellets (d = 3 mm), the suppression effect on flame propagation and pressure is obtained over a wider range of equivalence ratios, especially a better suppression effect is generated near the stoichiometric condition.     

Study on flame dynamics with secondary fuel injection control by large eddy simulation

In this study, flame dynamics of secondary fuel injection control is investigated by large eddy simulation (LES) in a lean premixed combustor to elucidate the experimental two-stage oscillation suppression phenomenon. Without secondary fuel injection, large-amplitude longitudinal oscillations were observed in the experiment. With constant injection of secondary fuel, the oscillation amplitude was reduced and with harmonic feedback injection, it was further reduced. In this process, the flame shape and the dynamic behavior of the flame were changed. To fill the gap of the experimental data and to understand the flame dynamics, large eddy simulations are carried out. The LES results show that the oscillation reduction is attributable to the roles of both the main flame and the secondary flame. The interaction of the main flame and the vortices is reduced when the injection is on; namely, heat fluctuation is reduced. The secondary flame helps the flame base stabilization and directly modulates the heat release in the feedback injection. The phase relations between the pressure, secondary heat release, and injection velocity are also shown. With these effects from the main and secondary flames combined, the oscillations are suppressed by secondary fuel injection. The present study indicates that LES can be a tool to understand flame dynamics with control.     

Effects of ignition location on premixed hydrogen/air flame propagation in a closed combustion tube

The dynamics of premixed hydrogen/air flame ignited at different locations in a finite-size closed tube is experimentally studied. The flame behaves differently in the experiments with different ignition positions. The ignition location exhibits an important impact on the flame behavior. When the flame is ignited at one of the tube ends, the heat losses to the end wall reduce the effective thermal expansion and moderate the flame propagation and acceleration. When the ignition source is at a short distance off one of the ends, the tulip flame dynamics closely agrees with that in the theory. And both the tulip and distorted tulip flames are more pronounced than those in the case with the ignition source placed at one of the ends. Besides, the flame–pressure wave coupling is quite strong and a second distorted tulip flame is generated. When the ignition source is in the tube center, the flame propagates in a much gentler way and the tulip flame can not be formed. The flame oscillations are weaker since the flame–pressure wave interaction is weaker.     

Effect of hydrogen enrichment on swirl/bluff-body lean premixed flame stabilization

This paper reports the mechanism of hydrogen enrichment in stabilizing swirl/bluff-body CH₄/air lean premixed flame. Large Eddy Simulation (LES) coupled with Thickened Flame (TF) model was performed to resolve the turbulent reacting flow. A detailed chemistry was used to describe the oxidization of CH₄/H₂/air mixtures. Particle Image Velocimetry (PIV) and Planar Laser-Induced Fluorescence of OH (OH-PLIF) simultaneous measurements were conducted to obtain the velocity fields and flame structures respectively. The numerical methods were validated by experimental data and showing good agreements. Both the experimental and numerical results show that, the flame brush attachment tends to leave the inner shear layer with increasing hydrogen addition, which will reduce the risk of flame lift-off. The chemical analyses prove that the attachment of CH₄/air flame is inherently weak. On the one hand, the CH₄/air flame is stabilized by the hot products inside the recirculation. On the other hand, the burnt gas suppresses the oxidation of H₂ and CO through H₂ + OH = H + H₂O and CO + OH = CO₂ + H, respectively. Although the proportion of CH₄ decomposition through CH₄ + OH = CH₃ + H₂O will be reduced by hydrogen addition, the path of CH₄ + H = CH₃ + H₂ will be enhanced significantly. Hydrogen addition will not only increase the overall reaction rate, but also change the combustion intensity at the nozzle exit from relatively weak to strong, which is also important for flame stabilization. The robust flame attachment obtained by hydrogen addition can attributed to the enhanced reactions of H₂ + OH = H + H₂O and CH₄ + H = CH₃ + H₂.     

Experimental study of flame propagation across a perforated plate

Flame propagation across a single perforated plate was experimentally studied in a square cross-section channel. Experiments were performed in premixed hydrogen-air mixture with different equivalence ratios and initial pressures, aiming at identifying the parametric influence. High-speed schlieren photography and pressure records were used to capture the flame front and obtain the pressure build-up. Four stages for the flame front crossing the perforated plate were obtained, namely, laminar flame, jet flame, turbulent flame and secondary flame front. Following ignition, a laminar flame was obtained, which was nearly not affected by the confinement. This laminar flame was squeezed to pass through the perforated plate, producing the jet flame with a step change on velocity. Turbulent flame was generated by merging the jets, which facilitated the acceleration of the flame front. Secondary flame front induced by Rayleigh-Taylor instability was clearly observed in the process of the turbulent front moving forward. Both velocity and pressure are enhanced in this stage. Parametric studies suggested that the secondary flame front is more obvious in the stoichiometric mixture with higher initial pressure, and characterized by a faster propagation velocity and a bigger pressure rise.     

Detailed simulations of the transient hydrogen mixing,leakage and flammability in air in simple geometries

During an accidental release, hydrogen disperses very quickly in air due to a relatively high density difference. A comprehensive understanding of the transient behavior of hydrogen mixing and the associated flammability limits in air is essential to support the fire safety and prevention guidelines. In this study, a buoyancy diffusion computational model is developed to simultaneously solve for the complete set of equations governing the unsteady flow of hydrogen. A simple vertical cylinder is considered to investigate the transient behavior of hydrogen mixing, especially at relatively short times, for different release scenarios: (i) the sudden release of hydrogen at the cylinder bottom into air with open, partially open, and closed tops, and (ii) small hydrogen jet leaks at the bottom into a closed geometry. Other cases involving the hydrogen releases/leaks at the cylinder top are also explored to quantify the relative roles of buoyancy and diffusion in the mixing process. The numerical simulations display the spatial and temporal distributions of hydrogen for all the configurations studied. The complex flow patterns demonstrate the fast formation of flammable zones with implications in the safe and efficient use of hydrogen in various applications.     

Large eddy simulation and finite-volume conditional moment closure modelling of a turbulent lifted H2/N2 flame

Large eddy simulations with three-dimensional finite-volume Conditional Moment Closure (CMC) model are performed for a hydrogen/nitrogen lifted flame with detailed chemical mechanism. The emphasis is laid on the influences of mesh resolution and convection scheme of finite-volume CMC equations on predictions of reactive scalars and unsteady flame dynamics. The results show that the lift-off height is underestimated and the reactive scalars are over-predicted with coarser CMC mesh. It is also found that further refinement of the CMC mesh would not considerably improve the results. The time sequences of the most reactive and stoichiometric hydroxyl radical mass fractions indicate that finer CMC mesh can capture more unsteady details than the coarser CMC mesh. Moreover, the coarse CMC mesh has lower conditional scalar dissipation rate, which would promote the earlier auto-ignition of the flame base. Besides, the effects of the convection scheme for the CMC equations (i.e., upwind, central differencing and their blends) on the lifted flame characteristics are also investigated. It is shown that different convection schemes lead to limited differences on the time-averaged temperature, mixture fraction and species mass fractions. Moreover, the root-mean square values of hydrogen and hydroxyl mass fractions show larger deviation from the measurements with hybrid upwind and central differencing scheme, especially around the flame base. Furthermore, the distributions of the numerical fluxes on the CMC faces also show obvious distinctions between the upwind and blending schemes. The budget analysis of the individual CMC terms shows that a sequence of CMC faces has comparable contributions with upwind scheme. However, with the hybrid schemes, the instantaneous flux is dominantly from limited CMC faces. The reactivity of a CMC cell is more easily to be affected by its neighbors when the upwind scheme is used.     

Experimental study and three-dimensional simulation of premixed hydrogen/air flame propagation in a closed duct

An experimental and numerical study of premixed hydrogen/air flame propagation in a closed duct is presented. High-speed schlieren photography is used in the experiment to record the changes in flame shape and location. The pressure transient during the combustion is measured using a pressure transducer. A dynamic thickened flame model is applied to model the premixed combustion in the numerical simulation. The four stages of the flame dynamics observed in the experiment are well reproduced in the numerical simulation. The oscillations of the flame speed and pressure growth, induced by the pressure wave, indicate that the pressure wave plays an important role in the combustion dynamics. The predicted pressure dynamics in the numerical simulation is also in good agreement with that in the experiment. The close correspondence between the numerical simulation and experiment demonstrate that the TF approach is quite reliable for the study of premixed hydrogen/air flame propagation in the closed duct. It is shown that the flame wrinkling is important for the flame dynamics at the later stages.     

Hydrogen decompression analysis by multi-stage Tesla valves for hydrogen fuel cell

Hydrogen fuel cell is an ideal power source for electric vehicles. For a hydrogen fuel cell electric vehicle, the hydrogen is reserved in a high pressure level to promote the recharge mileage while relatively low-pressure hydrogen is demanded for proper functioning of the fuel cell stack, so that decompression of hydrogen is needed before hydrogen flowing into the fuel cell. With a reverse flow through Tesla valves, there appears a large pressure drop between the inlet and outlet, which can be used for hydrogen decompression nicely. However, a single-stage Tesla valve cannot meet the pressure drop requirement, so multi-stage Tesla valves are utilized. In this paper, numerical simulations of reversed hydrogen flow through multi-stage Tesla valves are carried out. The stage number of multi-stage Tesla valves and the inlet/outlet pressure ratio are both studied, and the distributions of temperature, pressure, and velocity inside multi-stage Tesla valves are all investigated. Results show that as the stage number increases or the inlet/outlet pressure ratio decreases, the pressure and the velocity inside multi-stage Tesla valves decrease, and the less the stage number, the more possibility for the velocity higher than local acoustic speed. Besides, a power-law relationship between the flow rate, the stage number and pressure ratio is summarized.     

Buoyancy-driven single-sided natural ventilation in buildings with large openings

Full-scale experimental and computational fluid dynamics (CFD) methods were used to investigate buoyancy-driven single-sided natural ventilation with large openings. Detailed airflow characteristics inside and outside of the room and the ventilation rate were measured. The experimental data were used to validate two CFD models: Reynolds averaged Navier-Stokes equation (RANS) modeling and large eddy simulation (LES). LES provides better results than the RANS modeling. With LES, the mechanism of single-sided ventilation was examined by turbulence statistical analysis. It is found that most energy is contained in low-frequency regions, and mean flow fields play an important role.     

Experimental study on spontaneous ignition and flame propagation of high-pressure hydrogen release through tubes

An experimental study was conducted to research the mechanism of spontaneous ignition induced by high-pressure hydrogen release through tubes with a diameter of 10 mm and varying lengths from 0.3 to 3 m. The pressure and light signals inside the tube were collected. The propagation of shock wave inside and outside the tube was also systematically investigated. The development process of the jet flame in the atmosphere was completely recorded, and the multiple Mach disks at the tube exit were observed by using a high-speed camera. The results show that the minimum release pressure, at which the jet flame is formed, is found to be 3.87 MPa with the tube length of 1.7 m. When the tube length was longer than 1.7 m, the critical pressure for forming jet flame increased rapidly. The velocity attenuation of the shock wave is mainly affected by the burst pressure but not sensitive to the tube length, and the flame propagates in the tube at a slower velocity than the shock wave. The compression of the hydrogen-air mixture by the Mach disk causes it to burn more violently after passing through the Mach disk. It is confirmed that the flame at the tube exit is lifted in the atmosphere, then a jet flame initiates behind the second Mach disk.     

DNS of a premixed turbulent V flame and LES of a ducted flame using a FSD-PDF subgrid scale closure with FPI-tabulated chemistry

Two complementary simulations of premixed turbulent flames are discussed. Low Reynolds number two-dimensional direct numerical simulation of a premixed turbulent V flame is first performed, to further analyze the behavior of various flame quantities and to study key ingredients of premixed turbulent combustion modeling. Flame surface density, subgrid-scale variance of progress variables, and unresolved turbulent fluxes are analyzed. These simulations include fully detailed chemistry from a flame-generated tabulation (FPI) and the analysis focuses on the dynamics of the thin flame front. Then, a novel subgrid scale closure for large eddy simulation of premixed turbulent combustion (FSD-PDF) is proposed. It combines the flame surface density (FSD) approach with a presumed probability density function (PDF) of the progress variable that is used in FPI chemistry tabulation. The FSD is useful for introducing in the presumed PDF the influence of the spatially filtered thin reaction zone evolving within the subgrid. This is achieved via the exact relation between the PDF and the FSD. This relation involves the conditional filtered average of the magnitude of the gradient of the progress variable. In the modeling, this conditional filtered mean is approximated from the filtered gradient of the progress variable of the FPI laminar flame. Balance equations providing mean and variance of the progress variable together with the measure of the filtered gradient are used to presume the PDF. A three-dimensional larger Reynolds number flow configuration (ORACLES experiment) is then computed with FSD-PDF and the results are compared with measurements.     

Numerical studies on heat coupling and configuration optimization in an industrial hydrogen production reformer

Steam methane reforming furnaces are the most important devices in the hydrogen production industry. The highly endothermic reaction system requires reaction tubes in the furnace to have a large heat transfer area and to be operated under high temperature and pressure conditions. In order to enhance heat transfer efficiency and protect reaction tubes, the controlling and optimization of the furnace structure have increasingly received more and more research attention. As known from the furnace structure, it is essential to couple the exothermic combustion with the endothermic reforming reactions due to the highly interactive nature of the two processes. Thus, in this paper, the combustion process in the furnace was numerically studied by using computational fluid dynamics (CFD) to model the combustion chamber, coupled with methane steam reforming reaction inside the reaction tubes, defined by a plug flow model. A set of combustion models were compared for the furnace chamber and a plug flow reaction model was employed for reforming reaction tubes, and then a heat coupling process was established. The predicted flue gas temperature distribution showed that the heat transfer in the furnace was not uniform, resulting in hot spots and heat losses on the tube wall. Therefore, structure optimization schemes were proposed. Optimization on arrangements of the tubes and the nozzles promoted the uniform distribution of flue-gas temperature and then improved heat transfer efficiency, thereby enhancing performance of the steam reforming process.     

Pressure dynamics,self-ignition,and flame propagation of hydrogen jet discharged under high pressure

The self-ignition of hydrogen released from a high-pressure tank using extension tubes (2200 mm) with different diameters was studied. The processes of flame transition at a nozzle and jet flame development were characterized using a high-speed camera. The results indicated that the intensity of a shockwave and the Mach number decay faster in a 10-mm-diameter tube than that in a 15-mm-diameter tube. The pressure in a 15-mm-diameter tube was weaker than that in the 10-mm-diameter tube at the initial stage; however, it became higher in the later stage. Spontaneous ignition was more likely to happen in a 15-mm-diameter tube. The formation of a stabilized flame at the tube exit and Mach disk were observed during the transition of the flame to a jet fire. The stabilized flame showed a triangular shape because of the influence of a Prandtl–Meyer flow when a hydrogen jet entered a suddenly expanding environment. The formation and separation of a spherical flame were recorded during jet flame development. Large vortexes were formed in front of the flame because of the Kelvin–Helmholtz instability, which resulted in the separation of the spherical flame. The vortexes stopped rotating until the separated flame disappeared.     

Implementation of a dynamic thickened flame model for large eddy simulations of turbulent premixed combustion

A dynamic version of the thickened flame model for large eddy simulations (TFLES) of turbulent premixed combustion is implemented by coupling two independent codes through MPI (Message Passing Interface) running on massively parallel machines. The combustion code, AVBP, solves the usual filtered Navier–Stokes, mass fractions and energy balance equations on unstructured meshes while a second code takes advantage of the knowledge of resolved flow structures to automatically adjust the combustion model parameter from filtered resolved flow fields. This approach is tested against data from the F1, F2 and F3 pilot stabilised jet-flames experimentally investigated by Chen and coworkers. The model parameters determined dynamically are quite different for these cases, evidencing the decisive advantage of the dynamic procedure while statistical results are in good agreement with experimental data. The extra cost induced by the dynamic model is here about 20% but further optimisation remains possible.     

Numerical study of hydrogen dispersion in a fuel cell vehicle under the effect of ambient wind

In the rescue of hydrogen-fueled vehicle accidents, once accidental leakage occurs and hydrogen enters the cabin, the relatively closed environment of the vehicle is prone to hydrogen accumulation. Excessive hydrogen concentration inside the vehicle cabin may cause suffocation death of injured passengers and rescue crews, or explosion risk. Based on hydrogen fuel cell vehicle (HFCV) with hydrogen storage pressure 70 MPa, four different scenarios (i. with opened sunroof, ii. opened door windows, iii. opened sunroof and door windows and iv. opened sunroof, door windows and rear windshield) under the condition of accidental leakage were simulated using computational fluid dynamics (CFD) tools. The hydrogen concentration inside the vehicle and the distribution of flammable area (>4% hydrogen mole fraction) were analyzed, considering the effect of ambient wind. The results show that in the case of convection between interior and exterior of the vehicle via the sunroof, door windows or rear windshield, the distribution of hydrogen inside the vehicle is strongly affected by the ambient wind speed. In the least risk case, ambient wind can reduce the hydrogen mole fraction in the front of the vehicle to less than 4%, however the rear of the vehicle is always within flammable risk.