The HyRAM+ software toolkit provides a basis for conducting quantitative risk assessment and consequence modeling for hydrogen, natural gas, and autogas systems. HyRAM+ is designed to facilitate the use of state-of-the-art models to conduct robust, repeatable assessments of safety, hazards, and risk. HyRAM+ integrates deterministic and probabilistic models for quantifying leak sizes and rates, predicting physical effects, characterizing hazards (thermal effects from jet fires, overpressure effects from delayed ignition), and assessing impacts on people. HyRAM+ is developed at Sandia National Laboratories to support the development and revision of national and international codes and standards, and to provide developed models in a publicly-accessible toolkit usable by all stakeholders. This document provides a description of the methodology and models contained in HyRAM+ version 5.0. The most significant change for HyRAM+ version 5.0 from HyRAM+ version 4.1 is the ability to model blends of different fuels. HyRAM+ was previously only suitable for use with hydrogen, methane, or propane, with users having the ability to use methane as a proxy for natural gas and propane as a proxy for autogas/liquefied petroleum gas. In version 5.0, real natural gas or autogas compositions can be modeled as the fuel, or even blends of natural gas with hydrogen. These blends can be used in the standalone physics models, but not yet in the quantitative risk assessment mode of HyRAM+.
Previous research has provided strong evidence that CO2 and H2O gasification reactions can provide non-negligible contributions to the consumption rates of pulverized coal (pc) char during combustion, particularly in oxy-fuel environments. Fully quantifying the contribution of these gasification reactions has proven to be difficult, due to the dearth of knowledge of gasification rates at the elevated particle temperatures associated with typical pc char combustion processes, as well as the complex interaction of oxidation and gasification reactions. Gasification reactions tend to become more important at higher char particle temperatures (because of their high activation energy) and they tend to reduce pc oxidation due to their endothermicity (i.e. cooling effect). The work reported here attempts to quantify the influence of the gasification reaction of CO2 in a rigorous manner by combining experimental measurements of the particle temperatures and consumption rates of size-classified pc char particles in tailored oxy-fuel environments with simulations from a detailed reacting porous particle model. The results demonstrate that a specific gasification reaction rate relative to the oxidation rate (within an accuracy of approximately +/- 20% of the pre-exponential value), is consistent with the experimentally measured char particle temperatures and burnout rates in oxy-fuel combustion environments. Conversely, the results also show, in agreement with past calculations, that it is extremely difficult to construct a set of kinetics that does not substantially overpredict particle temperature increase in strongly oxygen-enriched N2 environments. This latter result is believed to result from deficiencies in standard oxidation mechanisms that fail to account for falloff in char oxidation rates at high temperatures.
Liquid hydrogen (LH2) used as a fuel onboard a heavy-duty vehicle can result in increased storage capacity and faster refueling relative to compressed gas. However, there are concerns about hydrogen losses from boil-off, potential safety issues, gaps in codes and standards for cryogenic hydrogen fuel, and technical challenges with LH2 systems for widespread transportation applications. A failure modes and effects analysis (FMEA), a safety codes and standards review, and a design review of the onboard liquid hydrogen system for a heavy-duty vehicle identified some of these potential safety issues and gaps in the codes and standards. The FMEA identified some medium and low risk failure points of the conceptual design, and the design review identified how carefully pressure relief needs to be considered for LH2 systems. In addition, a conceptual design for a LH2 refueling station was developed. Rough capital costs for the refueling station design were $\$1 million$ and the layout occupied approximately 13,000 ft2. These results can be used to inform future designs and analyses for LH2 heavy-duty vehicles.
The HyRAM+ software toolkit provides a basis for conducting quantitative risk assessment and consequence modeling for hydrogen, methane, and propane systems. HyRAM+ is designed to facilitate the use of state-of-the-art models to conduct robust, repeatable assessments of safety, hazards, and risk. HyRAM+ integrates deterministic and probabilistic models for quantifying accident scenarios, predicting physical effects, characterizing hazards (thermal effects from jet fires, overpressure effects from delayed ignition), and assessing impacts on people. HyRAM+ is developed at Sandia National Laboratories to support the development and revision of national and international codes and standards, and to provide developed models in a publicly-accessible toolkit usable by all stakeholders. This document provides a description of the methodology and models contained in HyRAM+ version 4.1. The two most significant changes for HyRAM+ version 4.1 from HyRAM+ version 4.0 are direct incorporation of unconfined overpressure into the QRA calculations and modification of the models for cryogenic liquid flow through an orifice. In QRA mode, the user no longer needs to input peak overpressure and impulse values that were calculated separately; rather, the unconfined overpressure is estimated for the given system inputs, leak size, and occupant location. The orifice flow model now solves for the maximum mass flux through the orifice at constant entropy while conserving energy, which does not require a direct speed of sound calculation. This does not affect the mass flow for all-gaseous releases; the method results in the same speed of sound for choked flow. However, this method does result in a higher (and more realistic) mass flow rate for a given leak size for liquid releases than was previously calculated.
The HyRAM+ software toolkit provides a basis for conducting quantitative risk assessment and consequence modeling for hydrogen, methane, and propane infrastructure and transportation systems. HyRAM+ is designed to facilitate the use of state-of-the-art science and engineering models to conduct robust, repeatable assessments of safety, hazards, and risk. HyRAM+ includes generic probabilities for equipment failures, probabilistic models for the impact of heat flux on humans and structures, and experimentally validated first-order models of release and flame physics. HyRAM+ integrates deterministic and probabilistic models for quantifying accident scenarios, predicting physical effects, and characterizing hazards (thermal effects from jet fires, overpressure effects from delayed ignition), and assessing impact on people and structures. HyRAM+ is developed at Sandia National Laboratories to support the development and revision of national and international codes and standards. HyRAM+ is a research software in active development and thus the models and data may change. This report will be updated at appropriate developmental intervals. This document provides a description of the methodology and models contained in HyRAM+ version 4.0. The most significant change for HyRAM+ version 4.0 from HyRAM version 3.1 is the incorporation of other alternative fuels, namely methane (as a proxy for natural gas) and propane into the toolkit. This change necessitated significant changes to the installable graphical user interface as well as changes to the back-end Python models. A second major change is the inclusion of physics models for the overpressure associated with the delayed ignition of an unconfined jet/plume of flammable gas.
There are several different calculation approaches and tools that can be used to evaluate the risk of hydrogen energy applications. A comparative study of Air Liquide’s ALDEA (Air Liquide Dispersion and Explosion Assessment) tools suite and Sandia’s HyRAM (Hydrogen Risk Assessment Models) toolkit has been conducted. The purpose of this study was to understand and evaluate the differences between the two calculation approaches, and identify areas for model improvements. There were several scenarios examined in this effort regarding hydrogen release dynamics. These scenarios include free jet release cases at varying pressures, vessel blowdown, and hydrogen build-up scenarios with and without ventilation. For each scenario, the input and output of the HyRAM calculations are documented, along with a comparison to the ALDEA results. Generally, the results from the two different tools were reasonably aligned. However, there were fundamental differences in evaluation methodology and functional limitations in HyRAM that caused discrepancies in some calculations.
Understanding liquid hydrogen tank fluid dynamics is key for modeling liquid hydrogen systems. The tank is the source for nearly all liquid hydrogen systems. Accurate flow modeling out of the tank is needed to predict flows through downstream components. Tank contains liquid and gas that may not be at equilibrium. Questions to be addressed are: Does heat and mass transfer between liquid and vapor affect the flow rate? Is boiling an important consideration? For what conditions is a pressure relief valve (PRV) sufficient to relieve pressure and when is the burst disc needed?
The Hydrogen Risk Assessment Models (HyRAM) software version 3 uses a real gas equation of state rather than the Abel-Noble equation of state that is used in 2.0 and previous versions. This change enables the use of HyRAM 3 for cryogenic hydrogen flows, whereas the Abel-Noble equation of state is not accurate at low temperatures. HyRAM 3.1 results were compared to experimental data from the literature in order to demonstrate the accuracy of the physics models. HyRAM 3.1 results were also compared to HyRAM 2.0 for high-pressure, non-cryogenic flows to highlight the differences in predictions between the two major versions of HyRAM. Validation data sets are from multiple groups and span the range of HyRAM physics models, including tank blowdown, unignited dispersion jet plume, ignited jet flame, and accumulation and overpressure inside an enclosure. Both versions 2.0 and 3.1 of HyRAM are accurate for predictions of blowdowns, diffusion jets, and diffusion flames of hydrogen at pressures up to 900 bar, and HyRAM 3.1 also shows good agreement with cryogenic hydrogen data. Overall, HyRAM 3.1 improves on the accuracy of the physical models relative to HyRAM 2.0. In most cases, this reduces the conservatism in risk calculations using HyRAM.
The HyRAM software toolkit provides a basis for conducting quantitative risk assessment and consequence modeling for hydrogen infrastructure and transportation systems. HyRAM is designed to facilitate the use of state-of-the-art science and engineering models to conduct robust, repeatable assessments of hydrogen safety, hazards, and risk. HyRAM includes generic probabilities for hydrogen equipment failures, probabilistic models for the impact of heat flux on humans and structures, and experimentally validated first-order models of hydrogen release and flame physics. HyRAM integrates deterministic and probabilistic models for quantifying accident scenarios, predicting physical effects, and characterizing hydrogen hazards (thermal effects from jet res, overpressure effects from deflagrations), and assessing impact on people and structures. HyRAM is developed at Sandia National Laboratories for the U.S. Department of Energy to increase access to technical data about hydrogen safety and to enable the use of that data to support development and revision of national and international codes and standards. HyRAM is a research software in active development and thus the models and data may change. This report will be updated at appropriate developmental intervals. This document provides a description of the methodology and models contained in HyRAM version 3.1. There have been several impactful updates since version 3.0. HyRAM 3.1 contains a correction to use the volume fraction for two-phase speed of sound calculations; this only affects cryogenic releases in which two-phase ow (vapor and liquid) is predicted in the orifice. Other changes include clarifications that inputs for tank pressure should be given in absolute pressure, not gauge pressure. Additionally, the interface now rejects invalid inputs to probability distributions, and the less accurate single-point radiative source model selection was removed from the interface.
Li, Xuefang; Chowdhury, Bikram R.; He, Qian; Christopher, David M.; Hecht, Ethan S.
Previous studies have shown that the two-layer model more accurately predicts hydrogen dispersion than the conventional notional nozzle models without significantly increasing the computational expense. However, the model was only validated for predicting the concentration distribution and has not been adequately validated for predicting the velocity distributions. In the present study, particle imaging velocimetry (PIV) was used to measure the velocity field of an underexpanded hydrogen jet released at 10 bar from a 1.5 mm diameter orifice. The two-layer model was the used to calculate the inlet conditions for a two-dimensional axisymmetric CFD model to simulate the hydrogen jet downstream of the Mach disk. The predicted velocity spreading and centerline decay rates agreed well with the PIV measurements. The predicted concentration distribution was consistent with data from previous planar Rayleigh scattering measurements used to verify the concentration distribution predictions in an earlier study. The jet spreading was also simulated using several widely used notional nozzle models combined with the integral plume model for comparison. These results show that the velocity and concentration distributions are both better predicted by the two-layer model than the notional nozzle models to complement previous studies verifying only the predicted concentration profiles. Thus, this study shows that the two-layer model can accurately predict the jet velocity distributions as well as the concentration distributions as verified earlier. Though more validation studies are needed to improve confidence in the model and increase the range of validity, the present work indicates that the two-layer model is a promising tool for fast, accurate predictions of the flow fields of underexpanded hydrogen jets.
Liquid hydrogen is increasingly being used as a delivery and storage medium for stations that provide compressed gaseous hydrogen for fuel cell electric vehicles. In efforts to provide scientific justification for separation distances for liquid hydrogen infrastructure in fire codes, the dispersion characteristics of cryogenic hydrogen jets (50–64 K) from high aspect ratio nozzles have been measured at 3 and 5 barabs stagnation pressures. These nozzles are more characteristic of unintended leaks, which would be expected to be cracks, rather than conventional round nozzles. Spontaneous Raman scattering was used to measure the concentration and temperature field along the major and minor axes. Within the field of interrogation, the axis-switching phenomena was not observed, but rather a self-similar Gaussian-profile flow regime similar to room temperature or cryogenic hydrogen releases through round nozzles. The concentration decay rate and half-widths for the planar cryogenic jets were found to be nominally equivalent to that of round nozzle cryogenic hydrogen jets indicating a similar flammable envelope. The results from these experiments will be used to validate models for cryogenic hydrogen dispersion that will be used for simulations of alternative scenarios and quantitative risk assessment.
Unintentional leaks at hydrogen fueling stations have the potential to form hydrogen jet flames, which pose a risk to people and infrastructure. The heat flux from these jet flames are often used to develop separation distances between hydrogen components and buildings, lot-lines, etc. The heat flux and visible flame length is well understood for releases from round nozzles, but real unintended leaks would be expected to be from higher aspect-ratio cracks. In this work, we measured the visible flame length and heat-flux characteristics of cryogenic hydrogen flames from high-aspect ratio nozzles. Heat flux measurements from 5 radiometers were used to assess the single-point vs the multi-point methods for interpretation of heat flux sensor data, finding the axial distance of the sensor for a single-point heat flux measurement to be important. We compare the flame length and heat flux data to flames of both cryogenic and compressed hydrogen from round nozzles. The aspect ratio of the release does not affect the flame length or heat flux significantly, for a given mass flow under the range of conditions studied. The engineering correlations presented in this work enable the prediction of flame length and heat flux which can be used to assess risk at hydrogen fueling stations with liquid hydrogen and develop science-based separation distances for these stations.
The Hydrogen Risk Assessment Model Plus (HyRAM+) toolkit combines quantitative risk assessment with simulations of unignited dispersion, ignited turbulent diffusion flames, and indoor accumulation with delayed ignition of fuels. HyRAM+ is differentiated from HyRAM in that it includes models and leak data for other alternate fuels. The models of the physical phenomena need to be validated for each of the fuels in the toolkit. This report shows the validation for propane which is being used as a surrogate for autogas, which is a mixture of propane and butane and used in internal combustion engines in vehicles. For flame length comparisons, five previously published experiments from peer reviewed journals were used to validate our models. The validation looked at flame lengths and flame widths with respect to different leak diameters, mass flow rates, and source pressures. Most of the sources included more than one set of experimental data, which were collected using different methods (CCD cameras, IR visualization etc.). In general, HyRAM+ overpredicts the flame lengths by around 65%. For heat and radiation models, we compared the heat flux and radiation data reported from two different sources to the values calculated by HyRAM+. For higher mass flow rates, the HyRAM+ calculated flame length results gave a better estimate of what is found in the experiments (65% error), but a higher error (85%) is observed between the HyRAM+ calculated lengths and the experimental flame lengthsfor lower mass flows. Some differences can be attributed to outdoor environmental effects (i.e. wind speed) and uncertainties in jet flame shapes. The propane flame trajectory is predicted for a high Reynolds number case with Re = 12,500 and a low Reynolds number case where Re = 2,000. The Re=12,500 case which is momentum dominated matches well with the experimental flame trajectory, but the agreement for the bouancy driven low Reynolds number case is not as good. Dispersion modeling for unignited propane was also analyzed. We compared the mole fraction, mixture fraction, mean velocity, concentration half width, and inverse mass concentration over an axial distance from different credible journals to the values calculated by HyRAM+. The results display good agreement but generally, HyRAM+ predicts a wider profile for mole fraction and mixture fraction experiments. Overall, HyRAM+’s results are reasonable for predicting the flame length, heat flux, flame trajectory, and dispersion for propane and can be used in risk analyses
Design and analysis of practical reactors utilizing solid feedstocks rely on reaction rate parameters that are typically generated in lab-scale reactors. Evaluation of the reaction rate information often relies on assumptions of uniform temperature, velocity, and species distributions in the reactor, in lieu of detailed measurements that provide local information. This assumption might be a source of substantial error, since reactor designs can impose significant inhomogeneities, leading to data misinterpretation. Spatially resolved reactor simulations help understand the key processes within the reactor and support the identification of severe variations of temperature, velocity, and species distributions. In this work, Sandia's pressurized entrained flow reactor is modeled to identify inhomogeneities in the reaction zone. Tracer particles are tracked through the reactor to estimate the residence times and burnout ratio of introduced coal char particles in gasifying environments. The results reveal a complex mixing environment for the cool gas and particles entering the reactor along the centerline and the main high-speed hot gas reactor flow. Furthermore, the computational fluid dynamics (CFD) results show that flow asymmetries are introduced through the use of a horizontal gas pre-heating section that connects to the vertical reactor tube. Computed particle temperatures and residence times in the reactor differ substantially from the idealized plug flow conditions typically evoked in interpreting experimental measurements. Furthermore, experimental measurements and CFD analysis of heat flow through porous refractory insulation suggest that for the investigated conditions (1350 °C, <20 atm), the thermal conductivity of the insulation does not increase substantially with increasing pressure.