Accidental release from a hydrogen car tank in a confined space like a tunnel poses safety concerns. This Computational Fluid Dynamics (CFD) study focuses on the first seconds of such a release, which are the most critical. Hydrogen leaks through a Thermal Pressure Relief Device (TPRD), forms a high-speed jet that impinges on the street, spreads horizontally, recirculates under the chassis and fills the area below it in about one second. The “fresh-air entrainment effect” at the back of the car changes the concentrations under the chassis and results in the creation of two “tongues” of hydrogen at the rear corners of the car. Two other tongues are formed near the front sides of the vehicle. In general, after a few seconds, hydrogen starts moving upwards around the car mainly in the form of buoyant blister-like structures. The average hydrogen volume concentrations below the car have a maximum of 71%, which occurs at 2 s. The largest “equivalent stoichiometric flammable gas cloud size Q9” is 20.2 m³ at 2.7 s. Smaller TPRDs result in smaller hydrogen flow rates and smaller buoyant structures that are closer to the car. The investigation of the hydrogen dispersion during the initial stages of the leak and the identification of the physical phenomena that occur can be useful for the design of experiments, for the determination of the TPRD characteristics, for potential safety measures and for understanding the further distribution of the hydrogen cloud in the tunnel.

The risks entailed by an accidental spill of Liquefied Natural Gas (LNG) should be indentified and evaluated, in order to design measures for prevention and mitigation in LNG terminals. For this purpose, simulations are considered a useful tool to study LNG spills and to understand the mechanisms that influence the vapor dispersion. In the present study, the ADREA-HF CFD code is employed to simulate the TEEX1 experiment. The experiment was carried out at the Brayton Fire Training Field, which is affiliated with the Texas A&M University system and involves LNG release and dispersion over water surface in open- obstructed environment. In the simulation the source was modeled as a two-phase jet enabling the prediction of both the vapor dispersion and the liquid pool spreading. The conservation equations for the mixture are solved along with the mass fraction for natural gas. Due to the low prevailing temperatures during the spill ambient humidity condenses and this might affect the vapor dispersion. This effect was examined in this work by solving an additional conservation equation for the water mass fraction. Two different models were tested: the hydrodynamic equilibrium model which assumes kinetic equilibrium between the phases and the non hydrodynamic equilibrium model, in order to assess the effect of slip velocity on the prediction. The slip velocity is defined as the difference between the liquid phase and the vapor phase and is calculated using the algebraic slip model. Constant droplet diameter of three different sizes and a lognormal distribution of the droplet diameter were applied and the results are discussed and compared with the measurements.