Natural and forced ventilation of buoyant gas released in a full-scale garage: Comparison of model predictions and experimental data

An increase in the number of hydrogen-fueled applications in the marketplace will require a better understanding of the potential for fires and explosion associated with the unintended release of hydrogen within a structure. Predicting the temporally evolving hydrogen concentration in a structure, with unknown release rates, leak sizes and leak locations is a challenging task. A simple analytical model was developed to predict the natural and forced mixing and dispersion of a buoyant gas released in a partially enclosed compartment with vents at multiple levels. The model is based on determining the instantaneous compartment over-pressure that drives the flow through the vents and assumes that the helium released under the automobile mixes fully with the surrounding air. Model predictions were compared with data from a series of experiments conducted to measure the volume fraction of a buoyant gas (at 8 different locations) released under an automobile placed in the center of a full-scale garage (6.8 m × 5.4 m × 2.4 m). Helium was used as a surrogate gas, for safety concerns. The rate of helium released under an automobile was scaled to represent 5 kg of hydrogen released over 4 h. CFD simulations were also performed to confirm the observed physical phenomena. Analytical model predictions for helium volume fraction compared favorably with measured experimental data for natural and forced ventilation. Parametric studies are presented to understand the effect of release rates, vent size and location on the predicted volume fraction in the garage. Results demonstrate the applicability of the model to effectively and rapidly reduce the flammable concentration of hydrogen in a compartment through forced ventilation.     

A study on the scaling features for mixing and deflagration potential of stratified layer of hydrogen due to molecular diffusion

Hydrogen deflagration in confined spaces is an important safety issue. The dispersion of a stratified layer of hydrogen due to molecular diffusion is studied. It represents an important class of problems related to long term behaviour of hydrogen release in confined spaces. Diffusion being a slow process, gives an upper bound on the time taken for the stratified layer to mix with air below. A method, based on four indices, namely, average mole fraction (of hydrogen), non-uniformity index, deflagration volume fraction and deflagration pressure ratio, developed recently by the authors, is used to provide vital temporal information on mixing of the stratified layer with air below and formation of flammable cloud in the enclosure. In the present paper, stratified layers of different thickness are considered and the temporal evolutions of the above indices are plotted against diffusion Fourier number. The results in non-dimensional form provide an upper bound of the time that would be required to form a uniform mixture and to attain a state with respect to deflagration potential for enclosures of different sizes. This estimate is an important input for planning mitigation measures before the accident and for post accident investigations.     

Formation of flammable hydrogen–air clouds from hydrogen leakage

The possibilities of the formation of a flammable cloud over the ground in an open atmosphere from the leakage of hydrogen stored at different temperatures are studied. The dispersion of hydrogen in the stable and unstable atmospheric conditions is determined using the Gaussian dispersion model. The efflux of hydrogen from the storage vessel is considered at velocities between 1 m/s and 1500 m/s, the latter corresponding to the upper limit of velocities arising from the choked flow. The dispersion analysis shows that flammable hydrogen–air clouds would not be formed over the ground under unstable atmospheric conditions for all efflux velocities and leakage areas and for the different temperatures of the hydrogen leak. However, under strongly stable atmospheric conditions, such as those associated with clear sky winter nights with low winds and temperature inversion in the planetary boundary layer, a flammable cloud is seen to be formed. This is particularly true for low temperature hydrogen efflux and very low velocities of the efflux.     

Effect of wind and buoyancy on hydrogen release and dispersion in a compartment with vents at multiple levels

The natural and forced mixing and dispersion of hydrogen released in an accidental manner in a partially enclosed compartment with vents at multiple heights is investigated using theoretical tools. The key to the analysis is determination of the position of the neutral buoyancy plane, where the pressure in the compartment is equal to that of the exterior. Air flows in through vents below the position of neutral buoyancy and exits from vents above it. CFD simulations are conducted to confirm the physical phenomena and to compare with the analytical results. The analytical model is useful in understanding the important physical processes involved during hydrogen release and dispersion in a compartment with vents at multiple levels, with and without a steady wind. Parametric studies are conducted to identify the relative importance of various parameters. Model results indicate that the steady-state hydrogen volume fraction in the compartment is lower when the hydrogen release rate is smaller and the vent cross-sectional area is larger. Results also indicate that the fastest way to reduce flammable levels of hydrogen concentration in a compartment can be accomplished by blowing through the vents.     

An inter-comparison exercise on the capabilities of CFD models to predict the short and long term distribution and mixing of hydrogen in a garage

The paper presents the results of the CFD inter-comparison exercise SBEP-V3, performed within the activity InsHyde, internal project of the HySafe network of excellence, in the framework of evaluating the capability of various CFD tools and modelling approaches in predicting the short and long term mixing and distribution of hydrogen releases in confined spaces. The experiment simulated was INERIS-TEST-6C, performed within the InsHyde project by INERIS, consisting of a 1 g/s vertical hydrogen release for 240 s from an orifice of 20 mm diameter into a rectangular room (garage) of dimensions 3.78 × 7.2 × 2.88 m in width, length and height respectively. Two small openings at the bottom of the front side of the room assured constant pressure conditions. During the test hydrogen concentration time histories were measured at 12 positions in the room, for a period up to 5160 s after the end of release, covering both the release and the subsequent diffusion phases. The benchmark was organized in two phases. The first phase consisted of blind simulations performed prior to the execution of the tests. The second phase consisted of post-calculations performed after the tests were concluded and the experimental results made available. The participation in the benchmark was high: 12 different organizations (2 non-HySafe partners), 10 different CFD codes and 8 different turbulence models. Large variation in predicted results was found in the first phase of the benchmark, between the various modelling approaches. This was attributed mainly to differences in turbulence models and numerical accuracy options (time/space resolution and discretization schemes). During the second phase of the benchmark the variation between predicted results was reduced.     

Modeling the diffusion of flammable hydrogen cloud under different liquid hydrogen leakage conditions in a hydrogen refueling station

Constructing hydrogen refueling stations will be popular for hydrogen energy use in the future, and investigating the diffusion characteristics of hydrogen in a leakage incident is quite significant. The instantaneous evolution of flammable hydrogen clouds arising from liquid hydrogen leakage in a hydrogen refueling station is predicted using Ansys Fluent, and parametric analyses are conducted to reveal the effects of storage pressure, source height, and leakage direction on the distributions of the flammable regions. In addition, the feasibilities of heating the ceiling or the ground of the station after the leakage of liquid hydrogen to accelerate the hydrogen dilution are examined. The results show that the flammable region is stabilized at 90 s, the corresponding flammable hydrogen cloud volume is about 333 m³, and the extensions of downwind and vertical directions reach 10 m and 9.3 m. Storage pressure has a finite effect on the downwind diffusion distance of the flammable cloud. A lower source height tends to format the high-concentration hydrogen cloud near the ground while a higher source height helps separate the flammable clouds from the ground. The upward leakage direction leads to the maximum downwind diffusion distance of about 10.2 m while the downward leakage direction makes the high hydrogen concentration region confined below the ceiling. Just maintaining the ceiling at the initial temperature of 300 K is effective for accelerating the hydrogen dilution in the upward leakage. The maximum hydrogen concentration and the flammable volume can be reduced at rates of 0.35 vol % and 8% for every 50 K increase in heating temperature. For the downward leakage, keeping the ground at the initial temperature just works for the first 40 s in reducing the maximum hydrogen concentration, while increasing the heating temperature receives a gradually declined effect on reducing the flammable volume.     

Dispersion modeling and assessment of natural gas containing hydrogen released from a damaged gas transmission pipeline

This paper performs a simulation and assessment of dispersion of natural gas containing hydrogen released from transmission pipeline using a Computational Fluid Dynamics (CFD) approach. A 3D CFD model is established to evaluate the dispersion behavior of hydrogen-enriched natural gas in the hydrogen-natural gas mixing station. The simulations include a matrix of scenarios for hydrogen doping ratios, gas release rates, wind speeds and wind directions. The development process of flammable gas cloud is predicted, and the dangerous area generated in the hydrogen-natural gas mixing station is assessed. Additionally, the effects of some critical factors on flammable gas dispersion behavior are analyzed. The simulations produce some useful outcomes including the parameters of flammable gas cloud and the dangerous area in the station, which are useful for conducting a prior risk assessment and contingency planning.     

Numerical investigation of the flammable extent of semi-confined hydrogen and methane jets

The effect of surfaces on the extent of high pressure horizontal unignited jets of hydrogen and methane is studied using computer fluid dynamics simulations performed with FLACS Hydrogen. Results for constant flow rate through a 6.35 mm diameter pressure relief Device (PRD) orifice from 100 barg, 250 barg, 400 barg, 550 barg and 700 barg compressed gas systems are presented for both horizontal hydrogen and methane jets. To quantify the effect of a horizontal surface on the jet, the jet exit is positioned at various heights above the ground ranging from 0.1 m to 10 m. Free jet simulations are performed for comparison purposes. Also, for cross-validation purposes, a number of cases for 100 barg releases were simulated using proprietary models developed for hydrogen within commercial CFD software PHOENICS. It is found that the presence of a surface and its proximity to the jet centreline result in a pronounced increase in the extent of the flammable cloud compared to a free jet.     

Numerical simulation of hydrogen leakage diffusion in seaport hydrogen refueling station

Ningbo's seaport hydrogen refueling station was used as the research object. The effects of different leakage angles, wind direction, roof shape, leakage hole diameters, temperature, and humidity on the diffusion of hydrogen leakage were studied by numerical simulation. The influence of leakage angle on hydrogen leakage is mainly reflected in the presence or absence of obstacles. The volume of the flammable hydrogen cloud was reduced by 31.16%, and the volume of the hazardous hydrogen cloud was reduced by 63.22% when there was no obstacle. The wind direction can significantly impact hydrogen leakage, with downwind and sidewind accelerating hydrogen discharge and reducing the risk. At the same time, headwind significantly increases the volume of the flammable hydrogen cloud. Compared with no wind, the volume of the flammable hydrogen cloud increased by 71.73% when headwind, but the volume of the hazardous hydrogen cloud decreased by 24.00%. If hydrogen shows signs of accumulation under the roof, the sloping roof can effectively reduce the hydrogen concentration under the roof and accelerate the hydrogen discharge. When the leakage angle θ = 90°, the sloping roof reduced the volume of the flammable hydrogen cloud by 11.74%. The leakage process was similar for different leak hole diameters in the no wind condition. The inverse of the molar fraction of hydrogen on the jet centerline was linearly related to the dimensionless axial distance of the jet in different cases. Using a least squares fit, the decay rate was obtained as 0.0039. In contrast, temperature and humidity have almost no effect on hydrogen diffusion. Hydrogen tends to accumulate on the lower surface of the roof, near the roof pillars and the hydrogen dispenser. In this paper, a set of hydrogen detector layout schemes was developed, and the alarm success rate was verified to be 83.33%.     

Evaluation of the net water transport through electrolytes in Proton Exchange Membrane Fuel Cell

Electro-osmotic drag and back diffusion are the primary water transport mechanisms in PEMFC (Proton Exchange Membrane Fuel Cell) electrolytes. These two phenomena occur competitively in the membrane, and ultimately determine the net water movement. The chemical compositions of reactants and product (i.e., H₂, O₂, and H₂O) in a porous catalyst layer vary with respect to the electro-chemical reactions and water transport through the membrane. The tendency of the chemical compositions was estimated by analyzing the net water transport coefficient (α), defined as the ratio of the reaction rate to the water transport rate. New criteria were suggested for predicting species mole fractions from the flow direction, and these were validated by CFD analysis. The hydrogen mole fraction had different tendencies to rise or fall based on the flow direction at α = 0.5, while the oxygen mole fraction extreme was located at α = −0.75. Finally, α was shown to influence the membrane conductivity and activation losses, which are the main factors that contribute to fuel cell performance.     

Hydrogen permeation from CGH2 vehicles in garages: CFD dispersion calculations and experimental validation

The time and space evolution of the distribution of hydrogen in confined settings was investigated computationally and experimentally for permeation from typical compressed gaseous hydrogen (CGH2) storage systems for buses or cars. The main goal was to examine whether hydrogen is distributed homogeneously within a garage-like facility or whether stratified conditions are developed, under certain conditions. The nominal hydrogen flow rate considered was 1.087 L/min in a bus facility with a volume of 681 m³. The release was assumed to be directed upwards from a 0.15 m diameter hole located at the middle part of the bus cylinders casing. Ventilation rates up to 0.03 air changes per hour (ACH) were considered. Simulated time periods extended up to 20 days. The CFD simulations performed with the ADREA-HF code showed that fully homogeneous conditions exist for low ventilation rates, while stratified conditions prevail for higher ventilation rates. Regarding flow structure it was found that the vertical concentration profiles can be considered as the superposition of the concentration at the floor (driven by diffusion) plus a concentration difference between floor and ceiling (driven by buoyancy forces). In all cases considered this concentration difference was found to be less than 0.5%. The dispersion experiments were performed in a large scale garage-like enclosure of 40 m³ using helium (GARAGE facility). Comparison between CFD simulations and experiments showed that the predicted concentrations were in good agreement with the experimental data. Finally, simulations were performed using two integral models: the fully homogeneous model and a two-layer model and the results were compared both against CFD and the experimental data.     

Comparisons of the hydrogen-rich syngas compositions from wet rice husk slurry steam reforming reactions using different catalysts

Rice husk slurry is pumped into a packed reactor and the products from the steam reforming reactions using different catalysts are studied. The steam/biomass weight ratio of such a system is between 3.47 and 5.25. The solids, liquid and gaseous products are a mass fraction of 2.8-4.1%, a mass fraction of 92.4-93.0% and a mass fraction of 3.5-4.7%, respectively. The hydrogen concentration in the gaseous product is approximate a volume fraction of 41% using the Al₂O₃ catalyst of a CuO mass fraction of 13%, a volume fraction of 38% using the Al₂O₃ catalyst of a Ni mass fraction of 13%, a volume fraction of 31% using the Al₂O₃ catalyst of a ZnO mass fraction of 13%, and a volume fraction of 20% using the Al₂O₃ catalyst at the reactor temperature of 800 °C. In the reactor temperature range studied (350-800 °C), the hydrogen concentration in the product stream increases monotonically with the increasing of the reactor temperature and the steam/carbon molar ratio. The value of dry gas LHV is between 9.4 MJ m⁻³ and 12 MJ m⁻³ at the reaction temperature of 600-800 °C. Considering the simple catalyst used in current study, the syngas of a hydrogen volume fraction of approximate 40% is obtained by pumping the biomass slurry to carry out the catalytic steam reforming reaction.     

Numerical investigation on the effects of load conditions and hydrogen-air ratio on the combustion processes of a HSDI engine

Owing to high specific energy and low emissions production, hydrogen is a desirable alternative fuel. The combustion and emission performance can be improved by hydrogen addition injected in-cylinder, intake manifold and aspirated with air. However, engine loads and hydrogen-air ration have a significant effect on the performance, combustion and emission of the diesel-hydrogen (high speed direct injection) HSDI engine. In this paper, the CFD method is used to calculate the combustion process of a diesel-hydrogen dual HSDI engine working at constant speed of 4000 rpm, at different hydrogen added from intake port (hydrogen volume fraction of 0%–10%) and five engine loads (equivalent to 20%, 40%, 60%, 80% and 100% of its maximum output power), respectively. The modelling results showed that the in-cylinder pressure and temperature under low engine load were more affected by hydrogen addition. With increasing hydrogen volume fraction, the indicated expansion work and in-cylinder peak pressure increased, and combustion duration decreased due to faster fuel-air mixing and spray flame speed.     

Measurements of sooting tendency in laminar diffusion flames of n-heptane at elevated pressure

This research focuses on the effects of an increasing pressure on the soot formation during combustion of vaporized liquid fuel. Therefore soot formation is measured in a laminar diffusion flame, with n-heptane as fuel, over a range of pressures from 1.0 to 3.0 bar. The soot volume fraction in the diffusion flames has been measured using Laser-Induced Incandescence (LII) calibrated by means of the Line Of Sight Attenuation (LOSA) technique. The values of the calibration factors between LII intensities and soot volume fraction from LOSA are slightly varied for different pressure. The integral soot volume fractions show power law dependence on pressures, being proportional to pⁿ, with n being 3.4 ± 0.3 in the pressure range of 1.0–3.0 bar.     

Effect of opening area on the suppression of self-ignition of high-pressure hydrogen gas leaking in the air by an extension tube

The most influential factor for self-ignition of high-pressure hydrogen is known to be the strength of the shock. Thus, the self-ignition can be suppressed by weakening the shock strength, which is possible by reducing the area where the hydrogen is ejected in this study. To confirm the possibility of this method, experiments were done by controlling the burst pressure of up to 302 bar and the ratio of the opening area. The experimental results showed that the minimum burst pressure of self-ignition is increased exponentially as the opening area is reduced. This confirmed that reducing the opening area under the same burst pressure conditions has an effect on the suppression of self-ignition. However, it was also found that the minimum shock speed that causes self-ignition gradually decreases as the opening area becomes smaller, which results from an increasing in mixing. The CFD simulation results showed that the volume of the flammable region in the tube was increased and the hydrogen-air mixing efficiency also increased when the opening area became smaller. The results suggest that reduction of the opening area can suppress a self-ignition by weakening the shock strength, but it should be noted that an increase in mixing effect also occurs.     

Helium dispersion following release in a 1/4-scale two-car residential garage

A series of experiments are described in which helium was released at a constant rate into a 1.5 m × 1.5 m × 0.75 m enclosure designed as a 1/4-scale model of a two-car garage. The purpose was to provide reference datasets for testing and validating computational fluid dynamics (CFD) models and to experimentally characterize the effects of a number of variables on the mixing behavior within an enclosure and the exchange of helium with the outside surroundings. Helium was used as a surrogate for hydrogen, and the total volume released was scaled as the amount that would be released by a typical hydrogen-fueled automobile with a full tank. Temporal profiles of helium were measured at seven vertical locations within the enclosure during and following 1-h and 4-h releases. Idealized vents in one wall sized to provide air exchange rates typical of actual garages were used. The effects of vent size, number, and location were investigated using three different vent combinations. The dependence on leak location was considered by releasing helium at three different points within the enclosure.     

Investigation of the hazardous area in a liquid hydrogen release with or without fence

Constructing a fence for potential liquid hydrogen leakage is a common protective measure to prevent the unrestricted spread, and critically evaluating the effects of a fence is essential for the safe use of hydrogen energy. Based on NASA's large-scale experiment of liquid hydrogen release, a 3D numerical model considering the existence of fence is proposed to study the dynamic spread behaviors of hydrogen flammable cloud and liquid hydrogen pool. The results show that the liquid hydrogen pool keeps a downwind radius of about 1.75 m most of the time during the release, and the hydrogen flammable cloud has a limit spread distance. Three cases are set to comparatively evaluate the effect of the fence, which are full fence, downwind half fence and no fence. Under the three cases, liquid hydrogen pools all evaporate completely within 5 s. Full fence scenario has the modest spread scale of hydrogen flammable cloud in the atmosphere but the longest duration for 86 s, and the upwind fence causes an easily overlooked hazardous region due to the vortex here, thus constructing a full fence is not an ideal strategy when the spread direction is constant. Just constructing a downwind half fence leads to the greatest spread distance of the flammable hydrogen cloud, which are 33.6 m and 68 m in vertical and downwind directions. No fence case corresponds to the smallest diffusion scale of flammable hydrogen cloud in atmosphere but the longest near-ground spread distance for 43.35 m. Four stages can be divided to summarize the entire life period of the flammable cloud for the cases with and without a fence, while three stages for the case with a downwind half fence.     

Exploratory studies of modeling approaches for hydrogen triple flames

A review of triple flame modeling is first presented, which demonstrates the need for additional work in this area. Building on previous methods described in the literature, a hybrid model that uses a weighted average of one-dimensional premixed and diffusion flamelet reaction rates has been proposed and evaluated for a hydrogen triple flame. Results indicated that some type of progress variable is needed for application of the diffusion flamelet contribution. Weighting the premixed flamelet reaction rate contribution at 100%, it is shown that peak temperatures between the model and a case employing detailed chemistry vary 7.5%, while heat release rate, flame speed, and mass fraction contours agree well.A second model, based on a library of reaction rates built from numerical studies which directly resolve the propagating triple flame has also been tested. Computational time for the baseline case is shown to be reduced by a factor of 3 ½ in comparison to use of detailed chemistry. The role of scalar dissipation rate as a necessary independent variable to the library has also been investigated using simulations with variable mixing layer thicknesses. Overall, it is found that large changes in local mixture fraction gradient cause rather small changes in propagation speed and total heat release rate of the hydrogen triple flame. This implies that such a model may be useful for CFD simulations that do not employ spatial resolution capable of resolving the triple flame itself.     

Turbulent non-premixed hydrogen-air flame structure in the pressure range of 1-10 atm

A numerical study of hydrogen turbulent diffusion flame structure is carried out in the pressure range of 1-10 atm with a special emphasis on mixing. The investigation is conducted under constant volumetric fuel and air flows. Mixing is characterized by mixture fraction, its variance and the scalar dissipation rate. The flow field and the chemistry are coupled by the flamelet assumption. Mixture fraction and its variance are transported by computational fluid dynamic (CFD). Computational predictions are analysed at two radial stations (the first one represent the near-field region and the second one the far-field region). The computational results indicate a deterioration of mixing with pressure rise. As a result, flame reaction zone becomes thicker. In addition, mixing and flame structure sensitivity to pressure are found to be high in the first location. Further analysis revealed that the gas becomes increasingly heavy with pressure rise, which hampered its ability to mix.     

Cathode electrocatalyst selection and deposition for a direct borohydride/hydrogen peroxide fuel cell

Catalyst selection, deposition method and substrate material selection are essential aspects for the design of efficient electrodes for fuel cells. Research is described to identify a potential catalyst for hydrogen peroxide reduction, an effective catalyst deposition method, and supporting material for a direct borohydride/hydrogen peroxide fuel cell. Several conclusions are reached. Using Pourbaix diagrams to guide experimental testing, gold is identified as an effective catalyst which minimizes gas evolution of hydrogen peroxide while providing high power density. Activated carbon cloth which features high surface area and high microporosity is found to be well suited for the supporting material for catalyst deposition. Electrodeposition and plasma sputtering deposition methods are compared to conventional techniques for depositing gold on diffusion layers. Both methods provide much higher power densities than the conventional method. The sputtering method however allows a much lower catalyst loading and well-dispersed deposits of nanoscale particles. Using these techniques, a peak power density of 680 mW cm⁻² is achieved at 60 °C with a direct borohydride/hydrogen peroxide fuel cell which employs palladium as the anode catalyst and gold as the cathode catalyst.