Safety

Correct use of Computational Fluid Dynamics (CFD) tools is essential in order to have confidence in the results. A comprehensive set of Best Practice Guidelines (BPG) in numerical simulations for Fuel Cells and Hydrogen applications has been one of the main outputs of the SUSANA project. These BPG focus on the practical needs of engineers in consultancies and industry undertaking CFD simulations or evaluating CFD simulation results in support of hazard/risk assessments of hydrogen facilities as well as on the needs of regulatory authorities. This contribution presents a summary of the BPG document. All crucial aspects of numerical simulations are addressed such as selection of the physical models domain design meshing boundary conditions and selection of numerical parameters. BPG cover all hydrogen safety relative phenomena i.e. release and dispersion ignition jet fire deflagration and detonation. A series of CFD benchmarking exercises are also presented serving as examples of appropriate modelling strategies.

The recent growth of the net of hydrogen fuelling stations increases the demands to transport compressed hydrogen on road by battery vehicles or tube-trailers both in composite pressure vessels. As a transport regulation the ADR is applicable in Europe and adjoined regions and is used for national transport in the EU. This regulation provides requirements based on the behaviour of each individual pressure vessel regardless of the pressure of the transported hydrogen and relevant consequences resulting from generally possible worst case scenarios such as sudden rupture. In 2012 the BAM (German Federal Institute for Materials Research and Testing) introduced consequence-dependent requirements and established them in national transport requirements concerning the “UN service life checks” etc. to consider the transported volume and pressure of gases. This results in a requirement that becomes more restrictive as the product of pressure and volume increases. In the studies presented here the safety measures for hydrogen road transport are identified and reviewed through a number of safety measures from countries including Japan the USA and China. Subsequently the failure consequences of using trailer vehicles the related risk and the chance are evaluated. A benefit-related risk criterion is suggested to add to regulations and to be defined as a safety goal in standards for hydrogen transport vehicles and for mounted pressure vessels. Finally an idea is given for generating probabilistic safety data and for highly efficient evaluation without a significant increase of effort.

Passive auto-catalytic recombiners (PARs) have the potential to be used in the future for the removal of accidentally released hydrogen inside confined areas. PARs could be operated both as stand-alone or backup safety devices e.g. in case of active ventilation failure.

Recently computational fluid dynamics (CFD) simulations have been performed in order to demonstrate the principal performance of a PAR during a postulated hydrogen release inside a car garage. This fundamental study has now been extended towards a variation of several boundary conditions including PAR location hydrogen release scenario and active venting operation. The goal of this enhanced study is to investigate the sensitivity of the PAR operational behaviour for changing boundary conditions and to support the identification of a suitable PAR positioning strategy. For the simulation of PAR operation the in-house code REKO-DIREKT has been implemented in the CFD code ANSYS-CFX 15.

In a first step the vertical position of the PAR and the thermal boundary conditions of the garage walls have been modified. In a subsequent step different hydrogen release modes have been simulated which result either in a hydrogen-rich layer underneath the ceiling or in a homogeneous hydrogen distribution inside the garage. Furthermore the interaction of active venting and PAR operation has been investigated.

As a result of this parameter study the optimum PAR location was identified to be close underneath the garage ceiling. In case of active venting failure the PAR efficiently reduces the flammable gas volume (hydrogen concentration > 4 vol.%) for both stratified and homogeneous distribution. However the simulations indicate that the simultaneous operation of active venting and PAR may in some cases reduce the overall efficiency of hydrogen removal. Consequently a well-matched arrangement of both safety systems is required in order to optimize the overall efficiency. The presented CFD-based approach is an appropriate tool to support the assessment of the efficiency of PAR application for plant design and safety considerations with regard to the use of hydrogen in confined areas.

Commercial use of hydrogen on-board fuel cell vehicles necessitates the compression of hydrogen gas up to 700 bar raising unique safety challenges. Potential hazards to be addressed include jet fires from high-pressure hydrogen on-board storage. Previous studies investigated effects of jet fires that occur when pressure relief devices (PRDs) on hydrogen fuel cell vehicles activate. This investigation examines plane jets’ axis switching and flame length accounting for compressibility effects and turbulent combustion near the point of release. Comparison with experimental data and previous plane jet simulation results reveal that combustion process does not affect flow dynamics in compressible region of jet flow. Furthermore a theoretical design of a variable aperture pressure relief device is examined which would enable the blow-down time to be minimized while reducing deterministic separation distances is examined using Computational Fluid Dynamics (CFD) techniques. Design recommendations are suggested for a novel PRD design.

A methodology for explosion and fire risk analyses in enclosed rooms is presented. The objectives of this analysis are to accurately predict the risks associated with hydrogen leaks in maritime applications and to use the approach to provide decision support regarding design and risk-prevention and risk mitigating measures. The methodology uses CFD tools and simpler consequence models for ventilation dispersion and explosion scenarios as well as updated frequency for leaks and ignition. Risk is then efficiently calculated with a Monte Carlo routine capturing the transient behavior of the leak. This makes it possible to efficiently obtain effects of sensitivities and design options maintaining safety and reducing costs.

The California Fuel Cell Partnership (CaFCP) is currently working with key stakeholders like the US Department of Energy (DOE) and National Fire Protection Association (NFPA) to develop a national template for educating and training first responders about hydrogen fuel cell-powered vehicles (FCV) and hydrogen fuelling infrastructure. Currently there are several existing programs that either have some related FCV/hydrogen material or have plans to incorporate this in the future. To create a robust national emergency responder (ER) program the strongest elements from these existing programs are considered for incorporation into the template. Working with the key stakeholders the national template will be evaluated on a regular basis to ensure accurate and up to date information and resources and effective teaching techniques for the emergency response community. This paper describes the evaluation process discusses elements of the template and reports on the steps and progress to implementation; all in the effort to effectively support the emergency response community as hydrogen infrastructure develops and FCVs are leased or sold.

Risk assessment of hazardous material transport by road is critical in considering the spatial features of the transport route. However previous studies that focused on hydrogen transport were unable to reflect the spatial features in their risk assessments. Hence this study aims to assess the risk of hydrogen transport by road considering the spatial features of the transport route based on a geographic information system (GIS). This risk assessment method is conducted through a case study in Yokohama which is an advanced city for hydrogen economy in Japan. In our assessment the risk determined by multiplying the frequency of accidents with the consequence was estimated by road segments that constitute the entire transport route. The effects of the road structure and traffic volumes were reflected in the estimation of the frequency and consequence for each road segment. All estimations of frequency consequence and risk were conducted on a GIS compiled with the information regarding the road network and population. In the case study in Yokohama the route for the transport of compressed hydrogen was virtually set from the near-term perspectives. Based on the case study results the risks of the target transport route were assessed at an acceptable level under the previous risk criteria. The results indicated that the risks fluctuated according to the road segments. This implies that the spatial features of the transport route significantly affect the corresponding risks. This finding corroborates the importance of considering spatial features in the risk assessment of hydrogen transport by road. Furthermore the discussion of this importance leads to the capability of introducing hydrogen energy careers with high transport efficiency and transport routing to avoid high risk road segments as risk countermeasures.

This paper describes a CFD study of a scenario involving the vertical downward release of hydrogen from a thermally-activated pressure relief device (TPRD) under a fuel cell car. The volumetric source model is applied to simulate hydrogen release dynamics during the tank blowdown process. Simulations are conducted for both unignited and ignited releases from onboard storage at 35 MPa and 70 MPa with TPRD orifice 4.2 mm. Results show that after TPRD opening the hazards associated with the release of hydrogen lasts less than two minutes and the most hazardous timeframe occurs within ten seconds of the initiation of the release. The deterministic separation distances for unignited releases are longer than those for ignited releases indicating that the separation distances are dominated by delayed ignition events rather than immediate ignition events. The deterministic separation distances for both unignited and ignited hydrogen downward releases under the car are significantly shorter than those of free jets. To ensure the safety of people a deterministic separation distance of at least 10 m for 35 MPa releases is required. This distance should be increased to 12 m for the 70 MPa release case. To ensure that the concentration of hydrogen is always less than 4% at the location of the air intake of buildings the deterministic separation distance should be at least 11 m for 35 MPa releases and 13 m for 70 MPa releases.

The U.S. Department of Energy supports the implementation of hydrogen fuel cell technologies by providing hydrogen safety and emergency response training to first responders. A collaboration was formed to develop and deliver a one-day course that uses a mobile fuel cell vehicle (FCV) burn prop designed and built by Kidde Fire Trainers. This paper describes the development of the training curriculum including the design and operation of the FCV prop; describes the successful delivery of this course to over 300 participants at three training centers in California; and discusses feedback and observations received on the course. Photographs and video clips of the training sessions will be presented.

This paper summarizes the current results of the theoretical and experimental activity carried out by the Italian Working Group on the fire prevention safety issues in the field of the hydrogen transport in pipelines. From the theoretical point of view a draft document has been produced beginning from the regulations in force on the natural gas pipelines; these have been reviewed corrected and integrated with the instructions suitable to the use of hydrogen. From the experimental point of view an apparatus has been designed and installed at the University of Pisa; this apparatus has allowed the simulation of hydrogen releases from a pipeline with and without ignition of hydrogen-air mixture. The experimental data have helped the completion of the above-mentioned draft document with the instructions about the safety distances. The document has been improved for example pipelines above ground (not buried) are allowed due to the knowledge acquired by means of the experimental campaign. The safety distances related to this kind of piping has been chosen on the base of risk analysis. The work on the text contents is concluded and the document is currently under discussion with the Italian stakeholders involved in the hydrogen applications.

In order to determine appropriate value for threshold stress intensity factor for hydrogen-assisted cracking (KIH) constant-displacement and rising-load tests were conducted in high-pressure hydrogen gas for JIS-SCM435 low alloy steel (Cr-Mo steel) used as stationary storage buffer of a hydrogen refuelling station with 0.2% proof strength and ultimate tensile strength equal to 772 MPa and 948 MPa respectively. Thresholds for crack arrest under constant displacement and for crack initiation under rising load were identified. The crack arrest threshold under constant displacement was 44.3 MPa m1/2 to 44.5 MPa m1/2 when small-scale yielding and plane-strain criteria were satisfied and the crack initiation threshold under rising load was 33.1 MPa m1/2 to 41.1 MPa m1/2 in 115 MPa hydrogen gas. The crack arrest threshold was roughly equivalent to the crack initiation threshold although the crack initiation threshold showed slightly more conservative values. It was considered that both test methods could be suitable to determine appropriate value for KIH for this material.

Quantitative Risk Assessment (QRA) can help to establish a set of design and operational requirements in hydrogen codes and standards that will ensure safe operation of hydrogen facilities. By analyzing a complete set of possible accidents in a QRA the risk drivers for these facilities can be identified. Accident prevention and mitigation features can then be analyzed to determine which are the most effective in addressing these risk drivers and thus reduce the risk from possible accidents. Accident prevention features/methods such as proper material selection and preventative maintenance are included in the design and operation of facilities. Accident mitigation features are included to reduce or terminate the potential consequences from unintended releases of hydrogen. Mitigation features can be either passive or active in nature. Passive features do not require any component to function in order to prevent or mitigate a hydrogen release. Examples of passive mitigation features include the use of separation distances barriers and flow limiting orifices. Active mitigation features initiate when specific conditions occur during an accident in order to terminate an accident or reduce its consequences. Examples of active mitigation features include detection and isolation systems fire suppression systems and purging systems. A concept being pursued by the National Fire Protection Association (NFPA) hydrogen standard development is to take credit for prevention and mitigation features as a means to reduce separation distances at hydrogen facilities. By utilizing other mitigation features the risk from accidents can be decreased and risk-informed separation distances can be reduced. This paper presents some preliminary QRA results where the risk reduction potential for several active and passive mitigation features was evaluated. These measures include automatic leak detection and isolation systems the use of flow limiting orifices and the use of barriers. Reducing the number of risk-significant components in a system was also evaluated as an accident prevention method. In addition the potential reduction in separation distances if such measures were incorporated at a facility was also determined.

Passive auto-catalytic recombiners (PARs) are used as safety devices in the containments of nuclear power plants (NPPs) for the removal of hydrogen that may be generated during specific reactor accident scenarios. In the presented study it was investigated whether a PAR designed for hydrogen removal inside a NPP containment would perform principally inside a typical surrounding of hydrogen or fuel cell applications. For this purpose a hydrogen release scenario inside a garage – based on experiments performed by CEA in the GARAGE facility (France) – has been simulated with and without PAR installation. For modelling the operational behaviour of the PAR the in-house code REKO-DIREKT was implemented in the CFD code ANSYS-CFX. The study was performed in three steps: First a helium release scenario was simulated and validated against experimental data. Second helium was replaced by hydrogen in the simulation. This step served as a reference case for the unmitigated scenario. Finally the numerical garage setup was enhanced with a commercial PAR model. The study shows that the PAR works efficiently by removing hydrogen and promoting mixing inside the garage. The hot exhaust plume promotes the formation of a thermal stratification that pushes the initial hydrogen rich gas downwards and in direction of the PAR inlet. The paper describes the code implementation and simulation results.

Conditions are examined under which mechanical stimuli produced by striking controlled blows can result in sparking and ignition of hydrogen in air mixtures. The investigation principally concerns magnesium thermite reaction as the ignition source and focuses on the conditions and thermomechanical parameters that are involved in determining the probability of ignition. It is concluded that the notion of using the kinetic energy of impact as the main criterion in determining whether an ignition event is likely or not is much less useful than considering the parameters which determine the maximum temperature produced in a mechanical stimuli event. The most influential parameter in determining ignition frequency or probability is the velocity of sliding movement during mechanical stimuli. It is also clear that the kinetic energy of a moving hammer head is of lesser importance than the normal force which is applied during contact. This explains the apparent discrepancy in previous studies between the minimum kinetic energy thought to be necessary to allow thermite sparking and gas ignition to occur with drop weight impacts and glancing blow impacts. In any analysis of the likelihood of mechanical stimuli to cause ignition the maximum surface temperature generated should be determined and considered in relation to the temperatures that would be required to initiate hot surface reactions sufficient to cause sparking and ignition.

In future pressure relief devices (PRDs) should be installed on hydrogen vehicles to prevent a hydrogen container burst in the event of a nearby fire. Weakening of the container at elevated temperature could result in such burst. In this case the role of a PRD is to release some or all of the system fluid in the event of an abnormally high pressure. The paper analyzes the possibility of hydrogen self-ignition at PRD operation and ways of its prevention.

During the early phase of the transient process following a hydrogen leak into the atmosphere a contact surface appears separating hot air from cold hydrogen. Locally the interface is approximately planar. Diffusion occurs potentially leading to ignition. This process was analyzed by Lin˜a´n and Crespo (1976) for Lewis number unity and Lin˜a´n and Williams (1993) for Lewis number less than unity. In addition  to conduction these processes are affected by expansion due to the flow which leads to a temperature drop. If chemistry is very temperature-sensitive then the reaction rate peaks close to the hot region where relatively little fuel is present. Indeed the Arrhenius rate drops rapidly as temperature drops much more so than fuel concentration. However the small fuel concentration present close to the airrich side depends crucially upon the balance between fuel diffusion and heat diffusion hence the fuel Lewis number. For Lewis number unity the fuel concentration present due to diffusion is comparable to the rate of consumption due to chemistry. If the Lewis number is less than unity fuel concentration brought in by diffusion is large compared with temperature-controlled chemistry. For a Lewis number greater than unity diffusion is not strong enough to bring in as much fuel as chemistry would be able to burn and combustion is controlled by fuel diffusion. In the former case combustion occurs faster leading to a localized ignition at a finite time determined by the analysis. As long as the temperature drop due to the expansion associated with the multidimensional nature of the jet does not lower significantly the reaction rate up to that point ignition in the jet takes place. For fuel Lewis number greater than unity first the reaction rate is much lower. Second chemistry does not lead to a defined ignition. Eventually expansion will affect the process and ignition does not take place. In summary it appears that the reason why hydrogen is the only fuel for which jet ignition has been observed is a Lewis number effect coupled with a high speed of sound hence a high initial temperature discontinuity.

Deflagration-to-Detonation Transition Ratio (DDTR) is an important parameter in measuring the hazard of hydrogen detonation at given thermodynamic conditions. It’s among the major tasks to evaluate DDTR in the study of hydrogen safety in a nuclear containment. With CFD tools detailed distribution of thermodynamic parameters at each instant can be simulated with considerable reliability. Then DDTR can be estimated using related CFD output. Forstochastic or epistemic reasons uncertainty always exists in input parameters during computations. This lack of accuracy can finally be reflected in the uncertainty of computation results e.g. DDTR in our consideration. The analysis of the influence of the input uncertainty is therefore a key step to understand the model’s response on the output and possibly to improve the accuracy. The increase of computational power makes it possible to perform statistics-based sensitivity and uncertainty (SU) analysis on CFD simulations. This paper aims at presenting some ideas on the procedure in safety analysis on hydrogen in nuclear containment. A hydrogen recombiner case is constructed and simulated with CFD method. DDTR at each instant is computed using a semi-empirical method. RBD-FAST based SU analysis is performed on the result.

Hydrogen has potential applications that require larger-scale storage use and handling systems than currently are employed in emerging-market fuel cell applications. These potential applications include hydrogen generation and storage systems that would support electrical grid systems. There has been extensive work evaluating regulations codes and standards (RCS) for the emerging fuel cell market such as the infrastructure required to support fuel cell electric vehicles. However there has not been a similar RCS evaluation and development process for these larger systems. This paper presents an evaluation of the existing RCS in the United States for large-scale systems and identifies potential RCS gaps. This analysis considers large-scale hydrogen technologies that are currently being employed in limited use but may be more widely used as large-scale applications expand. The paper also identifies areas of potential safety research that would need to be conducted to fill the RCS gaps. U.S. codes define bulk hydrogen storage systems but do not define large-scale systems. This paper evaluates potential applications to define a large-scale hydrogen system relative to the systems employed in emerging technologies such as hydrogen fuelling stations. These large-scale systems would likely be of similar size to or larger than industrial hydrogen systems.

The paper present the work on PIV-measurements of reactant flow velocity in front of propagating flames in hydrogen-air explosions. The experiments was performed with hydrogen-air mixture at atmospheric pressure and room temperature. The experimental rig was a square channel with 45 × 20 mm2 cross section 300 mm long with a single cylindrical obstacle of blockage ratio 1/3. The equipment used for the PIV-measurements was a Firefly diode laser from Oxford lasers Photron SA-Z high-speed camera and a particle seeder producing 1 μm droplets of water. The gas concentrations used in the experiments was 14 and 17 vol% hydrogen in air. The resulting explosion can be characterized as slow since the maximum flow velocity of the reactants was 13 m/s in the 14% mixture and 23 m/s in the 17% mixture. The maximum flow velocities was measured during the flame-vortex interaction and at two obstacle diameters behind the obstacle. The flame-vortex interaction contributed to the flame acceleration by increasing the overall reaction rate and the flow velocity. The flame area as a function of position is the same for both concentrations as the flame passes the obstacle.

Type IV pressure vessels are commonly used for hydrogen on-board stationary or bulk storages. During their lifetime they can be submitted to mechanical impacts creating damage within the composite structure not necessarily correlated to what is visible from the outside. When an impact is suspected or when a cylinder is periodically inspected it is necessary to determine whether it can safely stay in service or not. The FCH JU project Hypactor aims at creating a large database of impacts characterized by various non destructive testing (NDT) methods in order to provide reliable pass-fail criteria for damaged cylinders. This paper presents some of the tests results investigating short term burst) and long term (cycling) performance of impacted cylinders and the recommendations that can be made for impact testing and NDT criteria calibration.

This paper presents results from experiments on gas explosions in inhomogeneous hydrogen-air mixtures. The experimental channel is 3 m with a cross section of 100 mm by 100 mm and a 0.25 mm ID nozzle for hydrogen release into the channel. The channel is open in one end. Spectral analysis of the pressure in the channel is used to determine dynamic load factors for SDOF structures. The explosion pressures in the channel will fluctuate with several frequencies or modes and a theoretical high DLF is seen when the pressure frequencies and eigen frequencies of the structure matches.

The uncertainty on the laminar flame speed extracted from spherically expanding flames can be minimized by using large flame radius data for the extrapolation to zero stretch-rate. However at large radii the hydrodynamic and thermo-diffusive instabilities induce the formation of a complex cellular flame front and limit the range of usable data. In the present study we have employed the flame stability theory of Matalon to optimize the properties of the initial mixture so that transition to cellularity may occur at a pre-determined large radius. This approach was employed to measure the laminar flame speeds of H2/O2/N2/He mixtures with equivalence ratios from 0.6 to 2.0 at pressures of 50/80/100 kPa and a temperature of 300 K. For all the performed experiments the uncertainty related to the extrapolation to zero stretch-rate (performed with the linear curvature model) was below 2% as shown by the position of the data points in the (Lb/Rf;U Lb/Rf;L) plan where Lb is the burned Markstein length; and Rf;L and Rf;U are the flame radii at the lower and upper bounds of the extrapolation range. Comparison of the predictions of four chemical mechanisms with the present unstretched laminar flame speed data indicated an error below 10% for most conditions. In addition unsteady 1-D simulations performed with A-SURF demonstrated that the flame dynamical response to stretch rate could not be captured by the mechanisms. The present work indicates that although the stability theory of Matalon provides a well defined framework to optimize the mixture properties for improved flame speed measurement the uncertainty of some of the required parameters can result in largely over-estimated critical radius for cellularity onset which compromise the accuracy of the optimization procedure.

PEM water electrolysis offers an efficient and flexible way to produce “green-hydrogen” from renewable (intermittent) energy sources. Most research papers published in the open literature on the subject are addressing performances issues and to date very few information is available concerning the mechanisms of performance degradation and the associated consequences. Results reported in this communication have been used to analyze the failure mechanisms of PEM water electrolysis cells which can ultimately lead to the destruction of the electrolyzer. A two-step process involving firstly the local perforation of the solid polymer electrolyte followed secondly by the catalytic recombination of hydrogen and oxygen stored in the electrolysis compartments has been evidenced. The conditions leading to the onset of such mechanism are discussed and some preventive measures are proposed to avoid accidents.

A prospective global market share of Electric vehicle (EV) Hybrid electric vehicle (HEV) and Hydrogen Fuel Cell Vehicle (HFCV) is expected to grow due to stringent emission regulation and oil depletion. However it is essential to secure protection against high-pressure hydrogen gas and high-voltage in fuel cell vehicles and thus needs to develop a technique for safety assessment of HFCV. In this experiment 8 research institutes including the Korea Automobile Testing and Research Institute Hyundai Motor Company took part in. This project was supported by the Ministry of Land Transportation and Maritime Affairs of the Republic of Korea.

Fuel cell vehicles and some compressed natural gas vehicles are equipped with carbon fiber reinforced plastic (CFRP) composite cylinders. Each of the cylinders has a pressure relief device designed to detect heat and release the internal gas to prevent the cylinder from bursting in a vehicle fire accident. Yet in some accident situations the fire may be extinguished before the pressure relief device is activated leaving the high-pressure fuel gas inside the fire-damaged cylinder. To handle such a cylinder safely after an accident it is necessary that the cylinder keeps a sufficient post-fire strength against its internal gas pressure but in most cases it is difficult to accurately determine cylinder strength at the accident site. One way of solving this problem is to predetermine the post-fire burst strengths of cylinders by experiments. In this study automotive CFRP cylinders having no pressure relief device were exposed to a fire to the verge of bursting; then after the fire was extinguished the residual burst strengths and the overall physical state of the test cylinders were examined. The results indicated that the test cylinders all recorded a residual burst strength at least twice greater than their internal gas pressure for tested cylinders with new cylinder burst to nominal working pressure in the range 2.67–4.92 above the regulated ratio of 2.25.

From 2006 to 2009 INERIS alongside with CEA PSA PEUGEOT CITROËN and IRPHE were involved in a project called DRIVE. Its objective was to provide data on the whole reaction chain leading to a hydrogen hazard onboard a vehicle. Out of the three types of leakage identified by the consortium (permeation chronic and accidental) the chronic leakage taking place within the engine was judged to be more problematic since it can feature a high probability of occurrence and a significant release flow rate (up to 100 NL/min). Ignition tests carried out within a real and dummy engine compartment showed that pressure effects due to an explosion will be relatively modest provided that the averaged hydrogen concentration in this area is limited to 10% vol/vol which would correspond to a maximum release flow of 10 NL/min. This maximum concentration could be used as a threshold value for detection or as a target while designing the vehicle. Jet fire experiments were also conducted in the frame of the DRIVE project. It was found that pressure-relief devices (PRDs) might be unsuited to protect humans from the explosion of a tank caused by a bonfire. Other solutions are proposed in this paper.

In order to simulate an accidental hydrogen release from the low pressure pipe system of a hydrogen vehicle a systematic study on the nature of transient hydrogen jets into air and their combustion behaviour was performed at the FZK hydrogen test site HYKA. Horizontal unsteady hydrogen jets with an amount of hydrogen up to 60 STP dm3 and initial pressures of 5 and 16 bar have been investigated. The hydrogen jets were ignited with different ignition times and positions. The experiments provide new experimental data on pressure loads and heat releases resulting from the deflagration of hydrogen-air clouds formed by unsteady turbulent hydrogen jets released into a free environment. It is shown that the maximum pressure loads occur for ignition in a narrow position and time window. The possible hazard potential arising from an ignited free transient hydrogen jet is described.

This paper discusses the results of computational fluid dynamics (CFD) modelling of hydrogen releases and dispersion outdoors during venting of hydrogen storage in real environment and geometry of a hydrogen refuelling or energy station for a given flow rate and dimensions of vent stack. The PHOENICS CFD software package was used to solve the continuity momentum and concentration equations with the appropriate boundary conditions buoyancy model and turbulence models. Also thermal effects resulting from potential ignition of flammable hydrogen clouds were assessed using TNO “Yellow Book” recommended approaches. The obtained results were then applied to determine appropriate clearance distances for venting of hydrogen storage for contribution to code development and station design considerations. CFD modelling of hydrogen concentrations and TNO-based modelling of thermal effects have proven to be reliable effective and relatively inexpensive tools to evaluate the effects of hydrogen releases.

Hydrogen jet flames resulting from ignition of unintended releases can be extensive in length and pose significant radiation and impingement hazards. One possible mitigation strategy to reduce exposure to jet flames is to incorporate barriers around hydrogen storage and delivery equipment. While reducing the extent of unacceptable consequences the walls may introduce other hazards if not properly configured. This paper describes experiments carried out to characterize the effectiveness of different barrier wall configurations at reducing the hazards created by jet fires. The hazards that are evaluated are the generation of overpressure during ignition the thermal radiation produced by the jet flame and the effectiveness of the wall at deflecting the flame.<br/>The tests were conducted against a vertical wall (1-wall configuration) and two “3-wall” configurations that consisted of the same vertical wall with two side walls of the same dimensions angled at 135° and 90°. The hydrogen jet impinged on the center of the central wall in all cases. In terms of reducing the radiation heat flux behind the wall the 1-wall configuration performed best followed by the 3-wall 135° configuration and the 3-wall 90°. The reduced shielding efficiency of the three-wall configurations was probably due to the additional confinement created by the side walls that limited the escape of hot gases to the sides of the wall and forced the hot gases to travel over the top of the wall.<br/>The 3-wall barrier with 135° side walls exhibited the best overall performance. Overpressures produced on the release side of the wall were similar to those produced in the 1-wall configuration. The attenuation of overpressure and impulse behind the wall was comparable to that of the three-wall configuration with 90° side walls. The 3-wall 135° configuration’s ability to shield the back side of the wall from the heat flux emitted from the jet flame was comparable to the 1-wall and better than the 3-wall 90° configuration. The ratio of peak overpressure (from in front of the wall and from behind the wall) showed that the 3-wall 135° configuration and the 3-wall 90° configuration had a similar effectiveness. In terms of the pressure mitigation the 3-wall configurations performed significantly better than the 1-wall configuration

Clarification of questions of safety represents a decisive contribution to the successful introduction of vehicles fuelled by hydrogen. At the moment the safety of hydrogen is being discussed and investigated by various bodies. The primary focus is on fuel-cell vehicles with hydrogen stored in gaseous form. This paper looks at the safety of hydrogen-fuelled vehicles with an internal combustion engine and liquefied hydrogen storage. The safety concept of BMW’s hydrogen vehicles is described and the specific aspects of the propulsion and storage concepts discussed. The main discussion emphasis is on the utilization of boil-off parking of the vehicles in an enclosed space and their crash behaviour. Theoretical safety observations are complemented by the latest experimental and test results. Finally reference is made to the topic-areas in the field of hydrogen safety in which cooperative research work could make a valuable contribution to the future of the hydrogen-powered vehicle.

In the course towards a safe future hydrogen based society one of the tasks to be considered is the investigation of the conditions under which the use or storage of hydrogen systems inside buildings becomes too dangerous to be accepted. One of the relevant scenarios which is expected to have a relatively high risk is a slow (and long lasting) hydrogen release from a vehicle stored in a closed private garage without any forced ventilation i.e. only with natural ventilation. This scenario has been earlier investigated experimentally (by M. Swain) using He (helium) to simulate the hydrogen behavior. In the present work the CFD code ADREA-HF is used to simulate three of the abovementioned experiments using the standard k-  turbulence model. For each case modeled the predicted concentration (by vol.) time series are compared against the experimental at the given sensor locations. In addition the structure of the flow is investigated by presenting the helium concentration field.

In this paper a vertical turbulent round jet of helium was studied numerically using the PHOENICS software package. The flow was assumed to be steady incompressible and turbulent. The jet discharge Froude number was 14000 and the turbulent Schmidt number was 0.7. The incompressible Reynolds average Navier-Stokes equations and helium transport equation expressed in 2-D axisymmetric domain were applied to model the underlying helium release. The k-e RNG turbulence model was used for the calculations of the corresponding turbulent viscosity diffusivity velocity and concentration fields in the domain. The simulation results are compared with the experimental measurements from the earlier published studies on helium jets in non-buoyant jet region (NBJ) intermediate region (I) and buoyant plume region (BP). The numerical results show that the radial profiles of mean velocity and mean concentration are consistent with the empirical data scaled by the effective diameter and density-ratio dependence. The mean velocity and concentration fields along the axis of the jet agree with the decay laws correlated from the previous experiments. The discrepancy between the numerical and experimental data is within 10% proving that the current CFD model for gas release and dispersion is robust accurate and reliable and that the CFD technique can be used as an alternative to the experiments with similar helium jets. The authors believe that the current CFD model is well validated through this study and can be further extended to predict similar hydrogen releases and dispersion if the model is properly applied with hydrogen properties.

Wet dust removal systems used to control dust in the polishing or grinding process of Mg alloy products are frequently associated with potential hydrogen explosion caused by magnesium-water reaction. For purpose of avoiding hydrogen explosion risks we try to use a combination of chitosan (CS) and sodium phosphate (SP) to inhibit the hydrogen evolution reaction between magnesium alloy waste dust and water. The hydrogen evolution curves and chemical kinetics modeling for ten different mixing ratios demonstrate that 0.4% wt CS + 0.1% wt SP yields the best inhibition efficiency with hydrogen generation rate of almost zero. SEM and EDS analyses indicate that this composite inhibitor can create a uniform smooth tight protective film over the surface of the alloy dust particles. FTIR and XRD analysis of the chemical composition of the surface film show that this protective film contains CS and SP chemically adsorbed on the surface of ZK60 but no detectable Mg(OH)2 suggesting that magnesium-water reaction was totally blocked. Our new method offers a thorough solution to hydrogen explosion by inhibiting the hydrogen generation of magnesium alloy waste dust in a wet dust removal system.

Abatement of greenhouse gas emissions and diversification of energy sources will probably lead to an economy based on hydrogen. In order to evaluate safety conditions during transport and distribution experimental data is needed on the detonation of Hydrogen/Natural gas blend mixtures. The aim of this study is to constitute detonation and deflagration to detonation transition (DDT) database of H2/CH4-air mixtures. More precisely the detonability of such mixtures is evaluated by the detonation cell size and the DDT run up distance measurements. Large experimental conditions are investigated (i) various equivalence ratios from 0.6 to 3 (ii) various H2 molar fraction x ( ( )2 2 4x H H CH= + ) from 0.5 to 1 (iii) different initial pressure P0 from 0.2 to 2 bar at fixed ambient temperature T0=293 K. Detonation pressures P velocities D and cell sizes ? were measured in two smooth tubes with different i.d. d (52 and 106 mm). For DDT data minimum DDT run up distances LDDT were determined in the d=52 mm tube containing a 2.8 m long Schelkin spiral with a blockage ratio BR = 0.5 and a pitch equal to the diameter. Measured detonation velocities D are very close to the Chapman Jouguet values (DCJ). Concerning the effect of detonation cell size ? follows a classical U shaped- curve with a minimum close to =1 and concerning the effect of x ? decreases when x increases. The ratio ik L?= obtained from different chemical kinetics (Li being the ZND induction length) is well approximated by the value 40 in the range 0.5 < x < 0.9 and 50 for x 0.9. Minimum DDT run up distance LDDT varies from 0.36 to 1.1m when x varies from 1 to 0.8. The results show that LDDT obeys the linear law LDDT ~ 30-40? previously validated in H2/Air mixtures. Adding Hydrogen in Natural Gas promotes the detonability of the mixtures and for x 0.65 these mixtures are considered more sensitive than common heavy Alkane-Air mixtures.

This paper discusses the preliminary results of the Risk Management subtask efforts within the International Energy Agency (IEA) Hydrogen Implementing Agreement (HIA) Task 19 on Hydrogen Safety to develop uniform harm criteria for use in the Quantitative Risk Assessments (QRAs) of hydrogen facilities. The IEA HIA Task 19 efforts are focused on developing guidelines and criteria for performing QRAs of hydrogen facilities. The performance of QRAs requires that the level of harm that is represented in the risk evaluation be established using deterministic models. The level of harm is a function of the type and level of hazard. The principle hazard associated with hydrogen facilities is uncontrolled accumulation of hydrogen in (semi) confined spaces and consecutive ignition. Another significant hazard is combustion of accidentally released hydrogen gas or liquid which may or may not happen instantaneously. The primary consequences from fire hazards consist of personnel injuries or fatalities or facility and equipment damage due to high air temperatures radiant heat fluxes or direct contact with hydrogen flames. The possible consequences of explosions on humans and structures or equipment include blast wave overpressure effects impact from fragments generated by the explosion the collapse of buildings and the heat effects from subsequent fire balls. A harm criterion is used to translate the consequences of an accident evaluated from deterministic models to a probability of harm to people structures or components. Different methods can be used to establish harm criteria including the use of threshold consequence levels and continuous functions that relate the level of a hazard to a probability of damage. This paper presents a survey of harm criteria that can be utilized in QRAs and makes recommendations on the criteria that should be utilized for hydrogen-related hazards.

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 ¼-scale model of a two car garage. The purpose was to provide reference data sets for testing and validating computational fluid dynamics (CFD) models and to experimentally characterize the effects of a number of variables on the mixing behaviour within an enclosure and the exchange of helium with the 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 fuelled automobile with a full tank. Temporal profiles of helium were measured at seven vertical locations within the enclosure during and following one hour and four hour 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 from three different points within the enclosure. It is shown that the National Institute of Standards and Technology (NIST) CFD code Fire Dynamics Simulator (FDS) provides time resolved predictions for helium concentrations that agree well with the experimental measurements.

Computational fluid dynamics is used to simulate the release of high pressure Hydrogen from a reservoir with an exit of increasing diameter. Abel-Noble real gas equation of state is used to accurately simulate this high pressure release. Parallel processing based on Message Passing Interface for domain decomposition is employed to decrease the solution time. The release exit boundary is increased in time to simulate a scenario when the exit area increases during the release. All nodes and elements are moved accordingly at each time step to maintain the quality of the mesh. Different speeds of increasing diameter are investigated to see the impact on this unsteady flow.

High resolution numerical simulations used to study the mechanism of deflagration-to-detonation transition (DDT). The computations solved two-dimensional time-dependent reactive Navier-Stokes equations including the effects of compressibility molecular diffusion thermal conduction viscosity and detailed chemical kinetics for the reactive species with subsequent chain branching production of radicals and energy release. It is shown that from the beginning the flame accelerates exponentially producing shock waves far ahead. On the next stage the flame acceleration decreases and the shocks are formed close ahead of the flame front. The final stage is the actual transition to detonation. During the second stage a compressed unreacted mixture of increased density enters the flame producing a high pressure pulse which enhances reaction rate and the heat release in the reaction zone with a positive feedback coupling between the pressure pulse and the reaction rate. As a result the peak of the pressure pulse grows exponentially steepens into a strong shock which is coupled with the reaction zone forming the overdriven detonation. This new mechanism of DDT is different from the Zel’dovich’s gradient mechanism. The temperature gradients which appear in the form of hot spots and the like are not suitable to initiate detonation.

The effect of vertical surface on the extent of high pressure unignited jets of both hydrogen and methane is studied using computer fluid dynamics simulations performed with FLACS Hydrogen. Results for constant flow rate through a 6.35 mm round leak orifice from 100 barg 250 barg 400 barg 550 barg and 700 barg compressed gas systems are presented for vertical jets. To quantify the effect of the surface on the jet the jet exit is positioned at various distances from the surface ranging from 0.029 m to 12 m. Free jets simulations are performed for comparison purposes.

Because of air quality problems the problem of CO2 related greenhouse gas emissions and shortage of fossil fuels many vehicles with gaseous fuels (CNG biogas hydrogen etc.) are under research and development. Such vehicles have to prove that as well as their exhaust emissions their overall system emissions (including running loss) remain below certain safety limits before they can be used in practice. This paper presents a cost-effective way of monitoring such system emissions from hydrogen or other gaseous fuel powered vehicles within an air-conditioned chassis dynamometer test cell as commonly used for low ambient emission tests on gasoline vehicles. The only additional equipment needed is a low-concentration sensor for the gas of interest (e.g hydrogen). The method is based on concentration measurements and a dynamic mass balance model. This method is based on the fact that atoms cannot vanish. Applied to a room containing a gas mixture this means that the change of mass of a gaseous matter (called gas G subsequently) inside the chamber is the difference of all mass of G flowing into the chamber and all mass of G flowing out of the chamber. This assumes that no chemical reactions of the gas in mind with other matter take place. By measuring the flow rates and concentrations of ventilation-in flow and ventilation-out flow as well as room concentration the emissions of G of a source i.e. the vehicle to be tested can be calculated. These concentrations need to be measured as functions of time to be able to give values of emissions per time unit. It is shown by a real experiment that very low emissions can be recorded. Additionally error bounds and sensitivities on different parameters such as air exchange ratio are quantified.

The present work is concerned with numerical simulations of large scale hydrogen detonations. Euler equations have been solved along with a single step reaction for the chemistry. Total variation diminishing (TVD) numerical schemes are used for shock capturing. The equations are solved in parallel in a decomposed domain. Predictions were firstly conducted with a small domain to ensure that the reaction scheme has been properly tuned to capture the correct detonation pressure and velocity. On this basis simulations were set up for the detonation tests carried out at the RUT tunnel facilities in Russia. This is one of the standard benchmark test cases selected for HYSAFE [1]. Comparison is made between the predictions and measurements. Reasonably good agreement has been obtained on pressure decay and the propagation speed of detonation. Further simulations were then conducted for a hypothetical hydrogen-air cloud in the open to assess the impulse as well as overpressure. The effects of cloud height width were investigated in the safety context.

Buoyant round vertical jet had been investigated using Large Eddy Simulations at low Mach number. For the purpose of comparison with in-house experimental data in the present work helium has been used as a substitute for hydrogen. The influence of the transient concentration fields on the volume of gas with concentration within flammability limits has been investigated and their evolution and relation with average fields ad been characterized. Transient concentration fields created during initial jet development had been considered. Numerical results have been compared with in-house experiments and data published in the literature.

It is important to understand and quantify the potential risk resulting from accidental environmental exposure of condensed phase hydrogen storage materials under differing environmental exposure scenarios. This paper describes a modelling and experimental study with the aim of predicting consequences of the accidental release of 2LiBH₄+MgH₂ from hydrogen storage systems. The methodology and results developed in this work are directly applicable to any solid hydride material and/or accident scenario using appropriate boundary conditions and empirical data.

The ability to predict hydride behaviour for hypothesized accident scenarios facilitates an assessment of the risk associated with the utilization of a particular hydride. To this end an idealized finite volume model was developed to represent the behaviour of dispersed hydride from a breached system. Semi-empirical thermodynamic calculations and substantiating calorimetric experiments were performed in order to quantify the energy released energy release rates and to quantify the reaction products resulting from water and air exposure of a lithium borohydride and magnesium hydride combination.

The hydrides LiBH₄ and MgH₂ were studied individually in the as-received form and in the 2:1 “destabilized” mixture. Liquid water hydrolysis reactions were performed in a Calvet calorimeter equipped with a mixing cell using neutral water. Water vapor and oxygen gas phase reactivity measurements were performed at varying relative humidities and temperatures by modifying the calorimeter and utilizing a gas circulating flow cell apparatus. The results of these calorimetric measurements were compared with standardized United Nations (UN) based test results for air and water reactivity and used to develop quantitative kinetic expressions for hydrolysis and air oxidation in these systems. Thermodynamic parameters obtained from these tests were then inputted into a computational fluid dynamics model to predict both the hydrogen generation rates and concentrations along with localized temperature distributions. The results of these numerical simulations can be used to predict ignition events and the resultant conclusions will be discussed.

The possibility of using a risk based approach for the safe installation and siting of stationary fuel cell systems depends upon the availability of normative data and guidance on potential hazards and the probabilities of their occurrence. Such guidance data is readily available for most common hydrocarbon fuels. For hydrogen however data is still required on the hazards associated with different release scenarios. This data can then be related to the probability of different types of scenarios from historical fault data to allow safety distances to be defined and controlled using different techniques. Some data on releases has started to appear but this data generally relates to hydrogen vehicle refuelling systems that are designed for larger throughput higher pressures and the general use of larger pipe diameters than are likely to be used for small fuel cell systems.

As Europe is gradually moving towards a hydrogen based society it has been acknowledged that adding certain amount of hydrogen as a clean energy carrier to the existing natural gas pipeline will help reduce the CO2 emissions which contribute to the greenhouse effect. On the other hand hydrogen has been demonstrated to be able to change the behaviour of the pipeline steel such as lower toughness and faster crack growth due to hydrogen embrittlement. Therefore it is necessary that the risks associated with the failure of the pipeline carrying mixtures of natural gas and hydrogen be assessed.<br/>The study reported in this paper is part of European NATURALHY project whose aim is to investigate the possibility of using the existing natural gas transmission pipelines to convey natural gas/hydrogen mixtures. According to the EGIG database the most common cause of failure for the existing natural gas pipelines is third party damage which mainly refers to a gouge a dent/gouge combination of known geometry. Among third party damage failures 90% are the result of immediate failure i.e. leakage or rupture of the pipeline and only 10% of them are the result of delayed failure. While its not expected that hydrogen will impact the immediate failure it could increase the vulnerability of the pipe to delayed failure through the initiation or activation of crack like defects.<br/>This paper will present a methodology to predict the probability of increased failures and describe a software tool that has been developed to perform the calculations.

A benchmarking exercise on quantitative risk assessment (QRA) methodologies has been conducted within the project HyQRA under the framework of the European Network of Excellence (NoE) HySafe. The aim of the exercise was basically twofold: (i) to identify the differences and similarities in approaches in a QRA and their results for a hydrogen installation between nine participating partners representing a broad spectrum of background in QRA culture and history and (ii) to identify knowledge gaps in the various steps and parameters underlying the risk quantification. In the first step a reference case was defined: a virtual hydrogen refuelling station (HRS) in virtual surroundings comprising housing school shops and other vulnerable objects. All partners were requested to conduct a QRA according to their usual approach and experience. Basically participants were free to define representative release cases to apply models and frequency assessments according their own methodology and to present risk according to their usual format. To enable inter-comparison a required set of results data was prescribed like distances to specific thermal radiation levels from fires and distances to specific overpressure levels. Moreover complete documentation of assumptions base data and references was to be reported. It was not surprising that a wide range of results was obtained both in the applied approaches as well as in the quantitative outcomes and conclusions. This made it difficult to identify exactly which assumptions and parameters were responsible for the differences in results as the paper will show. A second phase was defined in which the QRA was determined by a more limited number of release cases (scenarios). The partners in the project agreed to assess specific scenarios in order to identify the differences in consequence assessment approaches. The results of this phase provide a better understanding of the influence of modelling assumptions and limitations on the eventual conclusions with regard to risk to on-site people and to the off-site public. This paper presents the results and conclusions of both stages of the exercise.

A buoyant round vertical hydrogen jet is investigated using Large Eddy Simulations at low Mach number (M = 0.3). The influence of the transient concentration fields on the extent of the gas envelope with concentrations within the flammability limits is analyzed and their structure are characterized. The transient flammable region has a complex structure that extends up to 30% beyond the time-averaged flammable volume with high concentration pockets that persist sufficiently long for potential ignition. Safety envelopes devised on the basis of simplified time-averaged simulations would need to include a correction factor that accounts for transient incursions of high flammability concentrations.

A 50-litre standard hydrogen gas cylinder was temporarily placed in a laboratory to supply hydrogen gas to a flame ionization detector (FID) for use in gas chromatography (GC). On 20 January 2015 the safety relief valve on the pressure regulator failed and released about 0.34 kg of hydrogen into the laboratory. The gas cloud did not ignite so there was no injury or damage. The results of a full investigation with a complete course of action and reconstruction are presented that verify the cause of the leakage and estimate the gas concentration of the dispersion and gas cloud. A preliminary simulation of the likely explosion is provided. If the gas cloud had ignited the explosion would most likely have caused significant structural damage to doors windows and possibly the walls.

The effect of irregular polygonal grain size and random grain boundary on hydrogen diffusion in polycrystalline nickel is investigated. Hydrogen diffusion behavior in micropolycrystalline nickel is compared with that in nanopolycrystalline nickel through numerical analysis. The two dimensional computational finite element microstructural and nanostructural analyses are based on Fick's law corresponding to heterogeneous polycrystalline model geometry. The heterogeneous polycrystalline model consists of random irregular polygonal grains. These grains are divided into internal grain and grain boundary regions the size of which is determined from the grain size. The computational analysis results show that hydrogen diffusion in nanostructural irregular polycrystalline nickel is higher in magnitude than the microstructural irregular polycrystalline nickel. However models of voids traps and micro and nano clustered grains are yet to be included.

On October 16-17 2012 the International Association for Hydrogen Safety (HySafe) in cooperation with the Institute for Energy and Transport of the Joint Research Centre of the European Commission (JRC IET Petten) held a two-day workshop dedicated to Hydrogen Safety Research Priorities. The workshop was hosted by Federal Institute for Materials Research and Testing (BAM) in Berlin Germany. The main idea of the Workshop was to bring together stakeholders who can address the existing knowledge gaps in the area of the hydrogen safety including identification and prioritization of such gaps from the standpoint of scientific knowledge both experimental and theoretical including numerical. The experience highlighting these gaps which was obtained during both practical applications (industry) and risk assessment should serve as reference point for further analysis. The program included two sections: knowledge gaps as they are addressed by industry and knowledge gaps and state-of-the-art by research. In the current work the main results of the workshop are summarized and analysed.