Next Article in Journal
Circular Economy and Sustainability-Oriented Innovation: Conceptual Framework and Energy Future Avenue
Next Article in Special Issue
Risk Assessment of a Hydrogen Refueling Station in an Urban Area
Previous Article in Journal
Study of the Characterisation Method of Effective Two-Phase Seepage Flow in the Construction of Gas Storage Reservoirs
Previous Article in Special Issue
Tracer Gas Test and CFD Analysis of Semiconductor Gas Box for Flammable Gas Leakage
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on the Explosion of the Hydrogen Fuel Tank of Fuel Cell Electric Vehicles in Semi-Enclosed Spaces

1
Fire Research Department, Korea Institute of Civil Engineering and Building Technology, Hwaseong-si 18544, Republic of Korea
2
Department of Mechanical and Automotive Engineering, Shinhan University, Uijeongbu-si 11644, Republic of Korea
3
Ju Sung G&B Inc., Seongnam-si 13209, Republic of Korea
*
Author to whom correspondence should be addressed.
Energies 2023, 16(1), 241; https://doi.org/10.3390/en16010241
Submission received: 10 November 2022 / Revised: 20 December 2022 / Accepted: 21 December 2022 / Published: 26 December 2022
(This article belongs to the Special Issue Latest Advances and Prospects of Hydrogen Safety)

Abstract

:
The rise in hydrogen fuel cell electric vehicles (FCEVs) is expected to pose a variety of hazards on the road. Vehicles using hydrogen could cause significant damage, owing to hydrogen vapor cloud explosions, jet fires caused by leakage, or hydrogen tank explosions. This risk is expected to further increase in semi-enclosed spaces, such as underground parking lots and road tunnels. Therefore, it is necessary to study the fire safety of hydrogen vehicles in semi-enclosed spaces. In this study, an experiment on hydrogen tank explosion was performed. In addition, the CFD numerical model was verified using the experimental results, and the damaging effect due to pressure propagation during hydrogen tank explosions in underground parking lots and road tunnels was analyzed using numerical analysis. From the experiment results, the hydrogen tank exploded at about 80 Mpa, a maximum incident pressure is generated 267 kPa at a distance of 1.9 m. As a result of numerical analysis based on the experimental results, the limit distance that can cause serious injury due to the explosion of a hydrogen tank in a road tunnel or underground parking lot was analyzed up to about 20 m from the point of explosion.

1. Introduction

As environmental problems accelerate, carbon neutrality policies established through international cooperation have been developed across all industry sectors worldwide. In particular, the use of conventional fossil-fuel-burning internal combustion engines in the automobile industry is being rapidly replaced with eco-friendly electric alternatives. In the EU, the proportion of eco-friendly vehicles is already as high as 30.7% owing to effective policies [1]. Furthermore, the United States is aggressively attempting to increase the proportion of zero-emission vehicles to 50% by 2030 to promote the eco-friendly vehicle industry [2]. In Korea, policies such as the “Hydrogen Economy Activation Roadmap” have been implemented to stimulate the relevant industries, secure technological process, and increase market penetration [3]. As a result, the number of hydrogen vehicles in Korea has gradually increased by 17 times over the period from 2018 to 2021 [4]. However, the transition to hydrogen fuel cell electric vehicles (FCEVs) is expected to pose different risks on the road. The characteristics of hydrogen mean that vehicles using it as their primary fuel could cause significant damage owing to hydrogen cloud explosions or jet fires caused by tank leakage, or even hydrogen tank explosions in the worst-case scenario. Such risks are further increased in semi-enclosed underground parking garages and road tunnels, potentially requiring technical mitigation measures.
LaFleur presented hydrogen FCEV accident scenarios by analyzing various risk factors, such as hydrogen leakage, jet fire, and fire explosion. According to their risk assessment for each accident scenario, the risk of a general fire or jet fire caused by hydrogen leakage was not high, whereas the explosion of hydrogen gas exhibited the highest risk [5].
Glover used a quantitative risk analysis of a road tunnel scenario to determine that the probability of explosion in the event of a hydrogen vehicle accident was 0.092%. These studies reported that hydrogen was relatively safe in the event of a general fire, but extremely dangerous when an explosion occurs [6]. Since the risk of hydrogen explosion increases in a semi-enclosed space, such as a tunnel or underground parking garage, studies have been conducted to predict and prevent accidents in these conditions.
To date, studies predicting the overpressure during a hydrogen explosion have generally been conducted using experimental and analytical methods.
In terms of experimental methods, Zalosh and Weyandt heated a 72 L Type II hydrogen tank pressurized to 34.3 MPa after removing the thermal pressure relief device (TPRD), detonated it through thermal deterioration, and measured the overpressure [7]. The hydrogen tank exploded after 6 min 27 s of heating with approximately 370 kW of energy, and at the moment of explosion, its internal pressure and temperature were 35.7 MPa and 39 °C, respectively. The maximum overpressures were measured to be 300, 83, 61, and 41 kPa at distances of 1.9, 4.2, and 6.5 m, respectively.
The Southwest Research Institute performed a hydrogen tank explosion experiment using propane to heat an 88 L sport utility vehicle-mounted hydrogen tank pressurized to 34.5 MPa. The tank exploded after 12 min 18 s of heating; the overpressures were measured to be 140, 56, 30, 14, and 5.5 kPa at 1.22, 2.5, 5, 9.75, and 24.38 m behind the vehicle, respectively, and 80 and 69 kPa at 2.5 and 5.0 m from the side of the vehicle, respectively [8]. Notably, when the above two experiments were performed by Weyandt, the overpressure decreased by half when the hydrogen tank mounted on a vehicle exploded compared to when the hydrogen tank was detonated alone. Furthermore, when the tank was mounted on a vehicle, there were significant differences depending on the measurement direction.
The equivalent TNT model and numerical analysis represent typical analytical methods applied to predict overpressure. Yoon applied the TNT equivalent method to evaluate the maximum overpressure and impulse produced by the detonation of explosive materials in solid and gaseous states, derived regression equations by fitting the maximum overpressure–converted distance curve, converted the impulse–converted distance curves obtained by a previous experiment under explosion yield conditions (0.01, 0.05, and 0.1), and compared them using the Kingery–Bulmash polynomial [9].
Kwon and Park investigated the TNT equivalent method to evaluate the energy produced by a gas explosion and reported a conservative prediction of overpressure [10].
Representative studies on the analytical prediction of overpressure in the event of a hydrogen explosion in a tunnel include work by Molkov, Shentsov, and Dery [11,12,13]. These researchers compared the results obtained using the equivalent TNT model and numerical analyses with experimental results, reporting that the blast pressure wave in a tunnel exhibited hemispherical propagation characteristics, similar to those observed in the atmosphere; this propagation wave was transformed into a plane wave when it reached a certain distance. The blast pressure wave also exhibited complex patterns as it transformed into a plane wave. These studies also summarized simulation results to present dimensionless correlations that can be used to predict the overpressure in tunnels.
The experiments reported in most previous studies were typically performed in an open space under conditions different from the specifications of hydrogen FCEVs currently available in the market. Therefore, in this study, an explosion experiment was performed to detonate the hydrogen tank of an actual FCEV and account for the conditions inside a semi-enclosed space. The experiment results were then analyzed to verify a computational fluid dynamics (CFD) numerical model, which was subsequently employed to simulate the damage owing to pressure wave propagation when a hydrogen tank explosion occurs in an underground parking garage or and road tunnel.

2. Hydrogen Tank Explosion Experiment

2.1. Experiment Overview

A hydrogen tank explosion experiment was performed to analyze the risk posed by a hydrogen FCEV explosion in a semi-enclosed space, such as a tunnel or underground parking garage. The results were analyzed to investigate the explosion characteristics inside the space as well as their influences on the surrounding environment, such as the atmosphere and structural components.

2.2. Experimental Method

The experiment was carried out in a cylindrical explosion test space with a diameter of 20 m and a height of 15 m. As shown in Figure 1, a hydrogen tank for an FCEV was installed in the center of the space and detonated by heating under specific temperature conditions according to Global Technical Regulation No. 13 (GTR 2013) for engulfing fire tests [14]. The hydrogen tank was a Type IV tank produced by company I (870 mm (L) × 363 mm (D)) that passed both pressure and heat resistance tests. The tank was fixed above the heating source using a wire and filled to a pressure of 70.7 MPa, which is consistent with the typical hydrogen FCEV filling condition. The TPRD was removed from the tank to ensure an explosion, and the tank was finished with an end plug in its place. Finally, a pressure gauge was installed in the tank to observe the change in its internal pressure during the heating process.
The heating source used in the experiment, shown in Figure 2, employed a liquid-petroleum-gas-fueled Bunsen-type burner to provide favorable temperature control and stabilization to meet the GTR (2013) heating condition [14]. The heating source was designed to engulf the entire container with a maximum flame length of 1650 mm, with the tank installed 100 mm above the heating source. Heating was then performed under the set conditions (800 to 1100 °C) until an explosion occurred. The heating temperature was measured by three thermocouples installed 25 mm away from the bottom of the tank (TC1~3). The average temperature of two points was calculated to confirm that the desired heating conditions were achieved. Accordingly, the temperature data were also measured by installing sheathed thermocouples at three points at 100 mm along the bottom centerline of the hydrogen tank.
A pressure gauge (Sensys PHP 1500 bar) and three thermocouples (K-type) were installed in the space to observe the hydrogen tank and heating conditions during the explosion, and the shock wave generated by the explosion was measured using multiple pressure gauges (incident pressure and reflected pressure) to capture its effects on the environment, considered in this study using a passenger car and two protective walls placed in the space, as shown in Figure 3. The incident pressure was measured using pencil-gauge-type static pressure sensors (PCB PIEZOTRONICS, Model 137B24B) placed 1.2 m above the ground and installed at horizontal distances of 1.9, 4.2, and 6.5 m from the hydrogen tank to the east, west, and north directions as in references [7,8,15]. The reflected pressure (PCB PIEZOTRONICS, Model 111A24) was measured at the center of the front door of a 2000 cc passenger car located 2 m to the south, at a concrete protective wall located 6.7 m to the northeast, and at a concrete protective wall located 7.2 m to the northwest of the hydrogen tank center. The DAQ used for pressure measurement was a DEWESoft SIRIUS series data logger, and it measured 500,000 data per second. Images were captured throughout the heating process using a set of general closed-circuit television cameras installed in the test space to provide visual evidence of the explosion as it occurred.

2.3. Hydrogen Tank Explosion Test Results

After 25 s of heating, the heating condition was satisfied as the temperatures of all thermocouples reached approximately 900 °C. The tank exploded after 735 s of heating at this condition, when the pressure in the tank dropped suddenly as shown in Figure 4. At the time of the explosion, the internal pressure was approximately 80 MPa, or 9.3 Mpa higher than the initial value. Figure 5 shows the whole view of the experiment before, during, and after the explosion of the hydrogen tank. When damage to the tank was examined after the explosion, the curved part on the side opposite the TPRD installation location was cracked, and the container had separated into two primary fragments. The relatively larger fragment (including the TPRD installation location) and smaller fragment had been thrown in different directions by the force of the explosion. The larger fragment struck the protective wall located to the northeast (6.9 m away), while the smaller fragment collided with the front door of the passenger car located 2 m adjacent to the tank. The protective wall was moved approximately 5 cm and the passenger car approximately 30 cm; there was considerable damage to the side of the vehicle, as shown in Figure 6. It appears, however, that this damage was caused by the shock wave as well as the fragment impact; if the exploded tank had been installed in a vehicle, the mounting device and the external vehicle structure may have resulted in different scattering directions and intensities.

2.3.1. Incident Pressure Results

Figure 7 and Table 1 show the incident pressure, maximum incident pressure, and impulse caused by the explosion of the hydrogen tank. The maximum incident pressures measured using the pressure gauges installed at the points east (PEs) of the hydrogen tank were found to be 144, 81, and 34 kPa at PE 1.9, 4.2, and 6.5 m, respectively. The maximum incident pressures measured at the points north (PNs) were 97, 62, and 44 kPa at PN 1.9, 4.2, and 6.5 m, respectively. The maximum incident pressures measured at the points west (PWs) were 267, 111, and 73 kPa at PW 1.9, 4.2, and 6.5 m, respectively. The highest maximum incident pressure occurred at PW 1.9 m; indeed, the measurements in the west direction were generally high owing to the pattern of tank damage, as it was judged that a relatively strong shock wave propagated in the westerly direction when the tank ruptured along the north–south direction. Accordingly, the north direction, in which the largest fragment was thrown, exhibited the lowest pressure.
The impulse is the sum of the pressure values applied during the duration and is a factor used together with the overpressure value when analyzing the risk due to an explosion. In contrast to the incident pressure results, the highest impulse occurred at PE 1.9 m, followed by PW 1.9 m, then PN 1.9 m. This result can be attributed to the propagation time of the pressure wave, which is expected to influence the damage to the human body when analyzed using the incident pressure–impulse diagram.

2.3.2. Reflected Pressure Results

Figure 8a shows the reflected pressures measured at the front door of the passenger car located 2 m south of the hydrogen tank. The pressure rapidly increased to a maximum of 960 kPa after the explosion owing to the generated shock wave, then decreased to a negative pressure (below atmospheric pressure). After approximately 4 ms, the pressure soared again, then an abnormal phenomenon indicating measurement error occurred. This appears to reflect a partial measurement of the tank fragment impact as it struck and damaged the installed reflected pressure sensor.
Figure 8b shows the reflected pressures measured by the sensors installed on the protective concrete walls to analyze the effects of the explosion on structures. Maximum reflected pressures of 94 and 85 kPa occurred at the protective wall 6.9 m to the northeast and 7.2 m to the northwest, respectively. These pressures were higher than the incident pressures measured at a relatively closer distance (PE, PN, and PW 6.5 m), because the total pressure, which includes the dynamic pressure, was measured as the reflected pressure at the walls.

3. CFD Analysis of Hydrogen Tank Explosion

The experiment results were used as basic data to construct a CFD model capable of analyzing the hydrogen FCEV explosion phenomenon under various conditions that cannot be easily performed in actual experiments owing to safety and budget constraints. First, the CFD analysis model was verified against the results of this and previous studies. Then, the CFD model was applied to simulate the effects of a hydrogen tank explosion in an underground parking garage and a tunnel.

3.1. CFD Analysis Model Verification

3.1.1. TNT Equivalent Model

The equivalent TNT model converts the target explosive material into an equivalent mass of TNT to predict the overpressure using empirical equations based on TNT and scaled distance as follows [16,17]:
M TNT = E H TNT
Z e = R M TNT 3
where MTNT is the equivalent mass of TNT, E is the explosion energy (kJ), HTNT is the heat of combustion of TNT (1100 kcal/kg), Ze is the scaled distance at which the TNT equivalent model is applied, R is universal gas constant, and MTNT is molecular mass of TNT. The value of E comprises the mechanical energy EM and chemical energy ECH considering the energy conversion efficiency as follows:
E = α E M + β E CH
where α is mechanical energy coefficient, and β is chemical energy coefficient.
During the explosion of a high-pressure hydrogen tank, combustion occurs only after the mechanical energy developed by the expansion pressure has formed into overpressure and been released. Therefore, this study assumed that the contribution of ECH to the explosion energy was insignificant, and EM was assumed to comprise the entire explosion energy when calculating the equivalent TNT mass in Equation (1); the value of EM was determined using the following equation proposed by Molkov [11]:
E M = ( P H 2 P atm ) ( V H 2 m H 2 · b ) γ 1
where P H 2 is the hydrogen pressure, Patm is the atmospheric pressure, V H 2 is the volume of the hydrogen tank (m3), parameter b is a co-volume constant, and 7.69 × 10−3 m3/kg was applied as a constant to correct the abnormal gas equation when calculating the actual gas volume, and m H 2 is calculated using:
m H 2 = ( P H 2 P H 2 · b + R · T H 2 M H 2 ) V H 2
where R is the gas constant (8.314 J/mol‧K), T H 2 is the internal temperature of the hydrogen tank (K), and M H 2 is the molecular weight of hydrogen molecules (2.016 g/mol).
The overpressure PEX, which spreads in all directions in the form of explosion energy, can then be obtained according to MTNT and Ze using an empirical equation or diagram.

3.1.2. Verification of the Equivalent TNT Model

In this section, the overpressure produced by the explosion of a hydrogen tank in the open atmosphere was obtained using a CFD numerical analysis performed according to, and verified through, comparison with previously obtained experimental data from Weyandt (2005) [7]. In the numerical analysis model, a 0.576 m diameter hydrogen tank (72 L volume) was located in the center of a 124.6 m diameter hemispherical space (modeled assuming 3D-axis-symmetric conditions) for consistency with the Weyandt experiment, as shown in Figure 9a. The hydrogen in the fuel tank was assumed to be an ideal gas, and the pressure and temperature immediately before the explosion were set to 35.8 MPa and 312 K, respectively, as measured during the Weyandt experiment. The hydrogen was set to be released simultaneously with the start of the simulation, which was performed in the unsteady state using time increments of 2 × 10−5 s. the simulation was performed using ANSYS Fluent V12 based code and numerical analysis was performed by applying the Spalart–Allmaras model with 3D-axis-symmetric conditions applied for the analysis model. For the sensitivity analysis of the mesh, a sufficient size was set (basic base lattice size: 0.012 m) considering the size of the source of the hydrogen tank explosion, and an appropriate ratio was applied to generate the mesh, considering that the propagation pattern was radioactive.
Figure 9b shows the simulated and experimental overpressures at distances of 1.9, 4.2, and 6.5 m from the explosion point over time. The errors between the maximum overpressures obtained by the numerical analysis and those obtained in the experiment ranged from 0.2 to 5.0% up to a distance of 4.2 m; at a distance of 6.5 m, the error was 31.2%. This is believed to have been affected by the pressure wave generated by the hydrogen tank’s scattering during the hydrogen explosion, but in the case of simulation, it is analyzed that the expansion of the pressure wave into the analysis area was simply numerical analysis.

3.2. Verification of Hydrogen Explosion in a Semi-Enclosed Space

A semi-enclosed space with a floor area of 400 m2 and height of 15 m was set as the analysis model space to verify the accuracy of the hydrogen tank explosion simulation in a semi-enclosed space. For comparison with the results of the experiment, domains of the same size as the experiment were selected and shown in Figure 10. The hydrogen tank was located in the center of the space, where the initial temperature and pressure of the air were set to 300 K and 101,325 Pa (atmospheric pressure), respectively. Immediately before the explosion, the hydrogen tank was filled with hydrogen with an initial temperature and pressure of 395 K and 70 MPa, respectively. The explosion energy was described using the equivalent TNT model derived in Section 3.1.1, and the hydrogen was set to be released simultaneously with the start of the simulation. The simulation was performed in the unsteady state using a time increment of 2 × 10−5 s. The numerical analysis method was performed in the same manner as the method presented in Section 3.1.2.
Figure 11a shows the hydrogen explosion pressure obtained by the numerical analysis for the semi-enclosed space condition. The highest pressure wave can be observed at distance of 1.9 m from the hydrogen tank, where a maximum pressure of 172 kPa was observed, followed by the waves at 4.2 and 6.5 m, where the maximum pressure was 78 kPa and 42 kPa, respectively. Furthermore, the generated pressure distribution contours along the horizontal plane and vertical height of the space at the detonation point are shown in Figure 11b.
Figure 11a compares the results of the hydrogen tank explosion test conducted in the semi-enclosed space in this study with those of the CFD numerical analysis when the incident pressure reached a specific location. The figure indicates that the simulated pressure wave propagating to the east slightly exceeded the experimental values, but similar results were predicted within an average error of 10%. There were larger differences in the results for the pressure wave propagating to the west, where the strongest shock wave was observed (average error: 57%), and to the north (average error: 20%). These errors appear to be the result of the explosion occurring after the internal pressure was increased by approximately 10 MPa owing to heating during the experiment, as well as the fact that pressure waves with different intensities were propagated in each direction owing to the non-uniform tank rupture. Indeed, as discussed in Section 2, two large fragments of the tank were thrown to the north and south of the explosion, with the larger fragment moving approximately 30° east of north. This observation suggests that cracking developed from the western side of the tank to the eastern side upon the release of pressure, causing the tank to be completely separated. As a result, the shock waves to the west and north were generated at an overpressure of approximately 80 MPa, a somewhat different result from that obtained by CFD simulation, which induced shock wave propagation using the pressure difference between the tank and the atmosphere (70 Mpa). In addition, the simulated overpressure exhibited an extremely small interval that converged to zero as its propagation was completed within 0.1 s of its generation. Therefore, a comparison of the maximum overpressure rather than the propagation pattern was determined to be more appropriate for this study.
Table 2 shows the maximum generating pressure between the numerical analysis and the experiment accordingly. If the distance is 6.5 m apart, the error between the numerical analysis and the experiment occurs up to 24%, which is believed to be because of the pressure generated by the scattering of the hydrogen tank container in the experiment, but only the expansion of the pressure wave in the simulation was analyzed.
As such, the maximum overpressure induced by the tank explosion was compared with the CFD analysis results and the results of previous studies as shown in Figure 12; the experimental conditions for each case shown in the figure are described in Table 3. The previous studies considered in this comparison comprise open-space explosion tests of different tank types filled to a pressure of 70 MPa conducted by Molcov et al. (2019) (Ref. Tests A and B) and the explosion test of a tank filled to a pressure of 35 MPa conducted by Molcov et al. (2015) (Ref. West and North) [11,18,19].
For Ref. Tests A and B, the initial tank pressure was the same, but the rupture pressure was approximately 14.5–19.5 MPa higher than the experiment results obtained in this study. In the Ref. West and North experiments, the initial tank pressure was 34.3 MPa and the tank rupture pressure was measured to be 35.7 MPa. At a distance of 1.9 m, the overpressure predicted by the simulation was somewhat underestimated compared to the experimental value for the north (PN) but was similar to the experimental value for the east (PE). At distances of 4.2 and 6.5 m, the pressure decreased following a similar slope in all three directions. When the overall errors were compared, the experimental values for PE were found to be in suitable agreement with the Ref. Test B (Type III tank) results.

3.3. CFD Analysis Model Verification

The CFD-based hydrogen tank explosion simulation employing the previously verified TNT equivalent model was applied to analyze the behavior of a hydrogen tank explosion in an underground parking garage and road tunnel, where actual hydrogen vehicles will travel, and the risk of hydrogen explosion exists.

3.3.1. Hydrogen Tank Explosion in an Underground Parking Garage

(1)
Analysis target and method
The pressure wave and overpressure propagation characteristics following a hydrogen tank explosion were analyzed in a 42 m × 42 m × 3 m (width × length × height) underground parking garage space, with the hydrogen tank placed at 21 m (center) in the width direction and 28 m (2/3 point) in the length direction. The capacity of the hydrogen tank was set to 52, 72, or 156 L. Immediately before the explosion, the hydrogen tank was filled with hydrogen at an assumed initial temperature of 395 K and pressure of 70 MPa. The initial temperature and pressure of the underground parking garage were set to 300 K and 101,325 Pa, respectively. In the analysis model, only half of the underground parking garage space was simulated by applying the symmetrical condition along the centerline in the width direction, as shown in Figure 13a. The fluid flow was assumed to be in the unsteady state and the time increment was set to 2 × 10−5 s for the numerical analysis. The numerical analysis method was performed in the same manner as the method presented in Section 3.1.2.
(2)
Analysis results
Figure 14a shows the maximum overpressure according to height (0.5, 1.0, 1.5, 2.0, and 2.5 m) and distance from the center of the underground parking garage when simulating the 52 L hydrogen tank explosion. The overpressure decreased as the horizontal distance from the tank increased and as the height above the floor increased. The maximum overpressures at 0.5, 1.0, 1.5, 2.0, and 2.5 m were 624, 323, 225, 173, and 88 kPa, respectively. The overpressure reduction rate according to the horizontal distance was 81% on average up to a distance of 3 m from the hydrogen tank and 58% from 3 to 10 m. The significant change in reduction rate at 3 m appears to have been caused by the interfering of the pressure wave with the reflected wave generated when the overpressure, radiated in hemispherical form, reached the ceiling (at a height of 3 m), resulting in a complex wave form.
Figure 14b shows the overpressure according to distance from the center of the parking garage at a height of 1.5 m for each hydrogen tank capacity. The maximum overpressure for the 52, 72, and 156 L tanks was found to be 225, 252, and 517 kPa, respectively, indicating that the maximum overpressure for the 72 and 156 L tanks was approximately 1.12 and 2.3 times higher than for the 52 L tank, respectively. When the tank capacity was 52 L, the overpressure sharply decreased up to a distance of 3 m and then gradually decreased.

3.3.2. Hydrogen Tank Explosion in a Road Tunnel

(1)
Analysis target and method
The pressure wave and overpressure propagation characteristics following a hydrogen tank explosion were analyzed using the verified CFD analysis model according to the cross-sectional area of the tunnel and hydrogen tank capacity. The length of the analyzed tunnel was set to 500 m. The tunnel cross section was horseshoe in shape, with an area of either 54 or 72 m2, modeled assuming symmetry about the central plane of the tunnel, as shown in Figure 15. The hydrogen tank was located in the bottom center of the tunnel and evaluated using capacities of 52, 72, or 156 L. The initial temperature and pressure of the air in the tunnel were set to 300 K and 101,325 Pa, respectively. The hydrogen tank was filled with hydrogen immediately before the explosion; its initial temperature and pressure were set to 395 K and 70 MPa, respectively. Once again, the fluid flow was assumed to be in the unsteady state and the time increment was set to 2 × 10−5 s for the numerical analysis. The grid used for analysis was set to a sufficiently small grid size (basic base grid size: 0.012 m) so that the periphery of the hydrogen fuel tank did not affect the quantitative analysis results, and was set in consideration of the analysis economy according to the relationship with the analysis time.
(2)
Analysis results
Figure 16a shows the overpressures along the tunnel length according to the hydrogen tank capacity at a height of 1.8 m when the tunnel cross-sectional area was 54 m2. Over the first approximately 16 m, the overpressures sharply decreased from their maximum values of 259, 348, and 680 kPa to 15.7, 20, and 47.5 kPa for hydrogen tank capacities of 52, 72, and 156 L, respectively. At distances beyond 16 m, the overpressure reduction rate was approximately −0.048 kPa/m on average regardless of tank capacity. Generally, the overpressures for the 72 and 152 L tanks were, respectively, 1.3 and 2.5 times higher than for 52 L tank.
Figure 16b compares the overpressures alone the tunnel length according to the hydrogen tank capacity when the tunnel cross-sectional area was 72 m2. Over the first approximately 18 m, the overpressures sharply decreased from their maximum values of 298, 404, and 752 kPa to 16, 19.5, and 35.4 kPa for hydrogen tank capacities of 52, 72, and 152 L, respectively. Beyond 18 m, the overpressures gradually decreased following nearly the same slope, corresponding to an average reduction rate of −0.035 kPa/m, where the overpressures for the 72 and 152 L tanks were approximately 1.3 and 2.5 times higher, respectively, than for the 52 L tank.

4. Analyzing the Risk Posed by a Hydrogen Tank Explosion in a Semi-Enclosed Space

Using the CFD-obtained pressure wave propagation characteristics presented in Section 3, the damage caused by a hydrogen tank explosion in an underground parking garage or road tunnel was analyzed according to the limit distances for harm to the human body. Many studies have been conducted on the impact of overpressure on the human body. Notably, Kashkarov defined the limit distance owing to a hydrogen tank explosion in an open space, estimating the damage to the human body by dividing the limit distance into four harm stages according to the degree of exposure, as shown in Table 4 [19]. The lower limit of slight injury was defined at 1.35 kPa from a conservative perspective to ensure safety.
Figure 17 shows the limit distances for each harm stage according to the capacity of a hydrogen tank exploding in an underground parking garage. When the tank capacity was 52 L, the fatality limit distance was 0.8 m, and the serious injury limit distance was 10.4 m. As the tank capacity increased to 72 and 156 L, the fatality limit distance increased to 1.0 and 3.2 m, respectively, and the serious injury limit distance increased to 11.2 and greater than 40 m, respectively. Indeed, the serious injury limit distance significantly increased for a 156 L hydrogen tank.
Figure 18 shows the limit distance results for each harm stage according to the capacity of a hydrogen tank exploding in a road tunnel. When the cross-sectional area of the tunnel was 54 m2, the fatality limit distances were 3.0, 3.5, and 9.5 m and the serious injury limit distances were 15.6, 27.3, and 240.5 m for hydrogen tank capacities of 52, 72, and 156 L, respectively. A similar tendency was also observed when the cross-sectional area was 72 m2, with fatality limit distances of 2.9, 3.5, and 10.1 m and serious injury limit distances of 17.6, 24.5, and 71.4 m for hydrogen tank capacities of 52, 72, and 156 L, respectively.
The results in the underground parking garage and road tunnel all demonstrated that the fatality limit distance increased linearly in proportion to the quantity of hydrogen in the tank, but the serious injury limit distance sharply increased when using 156 L of hydrogen. This result reflects the unique characteristics of the semi-enclosed space; the reflected wave generated when the overpressure wave created by the hydrogen explosion reaches the wall (or ceiling) overlaps with the pressure wave to exhibit a complex propagation pattern. The pressure reduction rate significantly decreases beyond the point where this overlap occurs according to the propagation distance. The serious injury limit distance significantly increases for a 156 L hydrogen tank because a large section of relatively high overpressure affected by the reflected wave is generated in which the pressure exceeds the value that determines the limit distance.
Figure 19 shows the limit distance according to the space condition for the explosion of a 52 L hydrogen tank. For the road tunnel, the fatality limit distance was determined to be approximately 3.0 m and the serious injury limit distance was 18.6–20.5 m regardless of the cross-sectional area. The limit distances for the underground parking garage were quite different, indicating a fatality limit distance of 0.8 m and a serious injury limit distance of 10.4 m. This difference appears to be influenced by the significantly lower ceiling height of the parking garage, but these results also differed considerably from those of previous experiments; therefore, detailed additional research is required.
Kashkarov and several other previous studies estimated damage using incident pressure–impulse diagrams [19]. Based only on the limit distance analysis presented above, nearly all human beings will die when exposed to an incident pressure of 200 kPa because lung hemorrhaging will occur with a 99% probability. However, if the impulse, which captures the concept of exposure time to the pressure waves, is used for analysis, more accurate and flexible damage prediction and engineering designs are possible. Figure 20 therefore examines harm using the incident pressure–impulse diagrams for the experiment and CFD results. In this diagram, the CFD results were compared at distances of 1.9, 4.2, and 6.5 m from the hydrogen tank to match the experimental conditions. All the CFD results for the underground parking garage and road tunnel were determined to be higher than the limit indicating death by lung hemorrhage. Among the experiment data, the lung hemorrhage threshold occurred at a distance of 4.2 m regardless of tank capacity, whereas serious injury was found to occur at 6.5 m. This result indicates that the impulse has a larger impact on the degree of damage than the generated incident pressure.

5. Conclusions

In this study, the characteristics of the pressure wave generated by the explosion of a hydrogen tank in a semi-enclosed space were analyzed using experiments and numerical analyses to examine the risks posed by a hydrogen FCEV explosion in an underground parking garage or road tunnel. In the experiment, a 52 L commercial hydrogen tank for an FCEV was filled with hydrogen to pressure of 70.7 MPa and exploded in an enclosed test space through heating. The resulting pressure wave propagation was measured in each direction at a variety of distances. A CFD numerical analysis of the hydrogen explosion energy was then conducted using the equivalent TNT model. The influence of the pressure wave propagated by a hydrogen explosion on the human body was subsequently analyzed in underground parking garage and road tunnel environments using the same technique. The following conclusions were identified based on the results:
  • When the Type IV hydrogen container was heated to a certain temperature (800–1100 °C) by an engulfing fire, it exploded after 735 s. At the time of the explosion, the pressure inside the container was 80.0 MPa, approximately 9.3 MPa higher than the initial pressure (70.7 MPa).
  • The largest incident pressure occurred in the direction (west) perpendicular to the separation plane of the hydrogen tank during the explosion; a maximum pressure of 267 kPa was measured by the western pressure gauge 1.9 m from the hydrogen tank. The maximum reflected pressure (967 kPa) was measured at the passenger car, which was closest to the hydrogen tank (2 m away). Reflected pressures of 94 and 85 kPa were observed at the protective walls located approximately 7 m from the hydrogen tank.
  • When the hydrogen tank exploded in the simulated underground parking garage and road tunnels, the overpressure wave exhibited hemispherical propagation until reaching the wall (or ceiling), at which time it formed a reflected wave that exhibited a complex pattern as it interacted with the pressure wave; afterward, the pressure reduction rate according to distance from the tank decreased.
  • For the underground parking garage and a given road tunnel cross-sectional area, the overpressure increased approximately 2.3 and 2.5 times, respectively, as the hydrogen tank capacity increased from 52 to 156 L.
  • The limit distances dictating the degree of harm to the human body increased as the hydrogen tank capacity increased. When the tank capacity was 52 L, the limit distance for fatality was approximately 1 m in the underground parking garage and 3.0 m in the road tunnel; the limit distance for serious injury was approximately 10 m in the underground parking garage and 18.65–20.5 m in the road tunnel.

Author Contributions

Conceptualization, J.P. and Y.Y.; methodology, J.P. and Y.Y.; experiment, J.P. and Y.Y.; software, J.R. and H.L.; validation, J.P. and H.L.; formal analysis, J.R.; data curation, J.P. and H.L.; writing—original draft preparation, J.P. and H.L.; writing—review and editing, J.P. and Y.Y. and J.R.; visualization, H.L.; supervision, Y.Y.; project administration, Y.Y.; funding acquisition, Y.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Fire Safety Agency’s ESS Hydrogen Facility Fire Safety Technology Research and Development Project (No. 20011645).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Korea Automotive Technology Institute. 2021 EV Sales Performance and Market Trends; Industry trend; Korea Automotive Technology Institute: Cheonan, Republic of Korea, 2022; p. 87. [Google Scholar]
  2. Korea Automotive Technology Institute. Review of US Eco-Friendly Car Industry Promotion Policy; Korea Automotive Technology Institute: Cheonan, Republic of Korea, 2021; p. 73. [Google Scholar]
  3. Ministry of Trade, Industry and Energy of the Republic of Korea. Hydrogen Economy Activation Roadmap; Ministry of Trade, Industry and Energy of the Republic of Korea: Sejong, Republic of Korea, 2019.
  4. Ministry of Trade, Industry and Energy of the Republic of Korea. 1st Basic Plan for Implementation of Hydrogen Economy; Ministry of Trade, Industry and Energy of the Republic of Korea: Sejong, Republic of Korea, 2021.
  5. LaFleur, C.B.; Bran Anleu, G.A.; Muna, A.B.; Ehrhart, B.D.; Blaylock, M.L.; Houf, W.G. Hydrogen Fuel Cell Electric Vehicle Tunnel Safety Study; SAND2017-11157; Sandia National Laboratories: Albuquerque, NM, USA, 2017; pp. 36–71.
  6. Glover, A.M.; Baird, A.R.; LaFleur, C.B. Hydrogen Fuel Cell Vehicles in Tunnels; SAND2020-4507 R; Sandia National Laboratories: Albuquerque, NM, USA, 2020; pp. 54–59.
  7. Zalosh, R.; Weyandt, N. Hydrogen fuel tank fire exposure burst test. SAE Trans. 2005, 114, 2338–2343. [Google Scholar]
  8. Weyandt, N. Vehicle Bonfire to Induce Catastrophic Failure of a 5,000-PSIG Hydrogen Cylinder Installed on a Typical SUV without Standard Required Safety Devices; Southwest Research Institute: San Antonio, TX, USA, 2006; p. 8. [Google Scholar]
  9. Yoon, Y.K. Evaluation of peak overpressure and impulse induced by explosion. J. Korean Soc. Explos. Blasting Eng. Explos. Blasting 2016, 34, 28–34. [Google Scholar]
  10. Kwon, S.; Park, J.C. A review of TNT equivalent method for evaluating explosion energy due to gas explosion. J. Korean Soc. Explos. Blasting Eng. Explos. Blasting 2015, 33, 1–13. [Google Scholar]
  11. Molkov, V.; Kashkarov, S. Blast wave from a high-pressure gas tank rupture in a fire: Stand-alone and under-vehicle hydrogen tanks. Int. J. Hydrogen Energy 2015, 40, 12581–12603. [Google Scholar] [CrossRef]
  12. Shentsov, V.; Makarov, D.; Molkov, V. Blast wave after hydrogen storage tank rupture in a tunnel fire. In Proceedings of the 8th International Symposium on Tunnel Safety and Security, Boras, Sweden, 14–16 March 2018; pp. 243–254. [Google Scholar]
  13. Molkov, V.; Dery, W. The blast wave decay correlation for hydrogen tank rupture in a tunnel fire. Int. J. Hydrogen Energy 2020, 45, 31289–31302. [Google Scholar] [CrossRef]
  14. United Nations. Global Technical Regulation No. 13; Global Technical Regulation on Hydrogen and Fuel Cell Vehicles; United Nations: Geneva, Switzerland, 2013. [Google Scholar]
  15. Shentsov, V.; Kim, W.; Makarov, D.; Molkov, V. Numerical Simulations of Experimental Fireball and Blast Wave from a High-Pressure Tank Rupture in a Fire. In Proceedings of the Eighth International Seminar on Fire and Explosion Hazards (ISFEH8), Hefei, China, 25 April 2016; p. 10. [Google Scholar]
  16. Nehzat, N. Gas Explosion Modeling for Complex Geometries. Ph.D. Thesis, University of New South Wales, Sydney, Australia, 1998. [Google Scholar]
  17. Kumar, A. Guidelines for Evaluating the Characteristics of Vapor Cloud Explosions, Flash Fires, and BLEVEs; Center for Chemical Process Safety of AIChE: New York, NY, USA, 1994. [Google Scholar]
  18. Molkov, V.; Cirrone, D.; Shentsov, V.; Dery, W.; Kim, W.; Makarov, D. Blast wave and fireball after hydrogen tank rupture in a fire. In Advances in Pulsed and Continuous Detonations; Torus Press: Newtownabbey, UK, 2018; p. 16. [Google Scholar]
  19. Molkov, V.; Cirrone, D.; Shentsov, V.; Dery, W.; Kim, W.; Makarov, D. Dynamics of blast wave and fireball after hydrogen tank rupture in a fire in the open atmosphere. Int. J. Hydrogen Energy 2020, 46, 4644–4665. [Google Scholar] [CrossRef]
Figure 1. Schematic of experimental setup.
Figure 1. Schematic of experimental setup.
Energies 16 00241 g001
Figure 2. Hydrogen tank and heating source: (a) hydrogen tank (type Ⅳ); (b) heating source.
Figure 2. Hydrogen tank and heating source: (a) hydrogen tank (type Ⅳ); (b) heating source.
Energies 16 00241 g002
Figure 3. Sensors installed for the experiment: (a) incident pressure sensor; (b) reflect pressure sensor.
Figure 3. Sensors installed for the experiment: (a) incident pressure sensor; (b) reflect pressure sensor.
Energies 16 00241 g003
Figure 4. Tank internal pressure and heating temperature during the test: (a) internal pressure of hydrogen tank during heating; (b) temperature of heating source during experiment.
Figure 4. Tank internal pressure and heating temperature during the test: (a) internal pressure of hydrogen tank during heating; (b) temperature of heating source during experiment.
Energies 16 00241 g004
Figure 5. Scenes of the hydrogen tank explosion test in each direction: (a) left side; (b) right side; (c) front side.
Figure 5. Scenes of the hydrogen tank explosion test in each direction: (a) left side; (b) right side; (c) front side.
Energies 16 00241 g005
Figure 6. Observed damage to the hydrogen tank and vehicle: (a) exploded hydrogen tank; (b) vehicle damage.
Figure 6. Observed damage to the hydrogen tank and vehicle: (a) exploded hydrogen tank; (b) vehicle damage.
Energies 16 00241 g006
Figure 7. Incident pressures induced by the hydrogen tank explosion: (a) east direction; (b) north direction; (c) west direction.
Figure 7. Incident pressures induced by the hydrogen tank explosion: (a) east direction; (b) north direction; (c) west direction.
Energies 16 00241 g007
Figure 8. Reflected pressures induced by the hydrogen tank explosion: (a) at passenger car; (b) at concrete walls.
Figure 8. Reflected pressures induced by the hydrogen tank explosion: (a) at passenger car; (b) at concrete walls.
Energies 16 00241 g008
Figure 9. Verification of CFD analysis model using the open explosion condition: (a) numerical analysis model schematic; (b) comparison of simulation and experimental [7] results.
Figure 9. Verification of CFD analysis model using the open explosion condition: (a) numerical analysis model schematic; (b) comparison of simulation and experimental [7] results.
Energies 16 00241 g009
Figure 10. Schematic of simulation model with photograph for comparison.
Figure 10. Schematic of simulation model with photograph for comparison.
Energies 16 00241 g010
Figure 11. Simulation results for a hydrogen explosion in a semi-enclosed space: (a) pressure profile; (b) pressure contours—floor plane (left) and vertical plane (right).
Figure 11. Simulation results for a hydrogen explosion in a semi-enclosed space: (a) pressure profile; (b) pressure contours—floor plane (left) and vertical plane (right).
Energies 16 00241 g011
Figure 12. Comparison of maximum overpressures generated by hydrogen tank explosions.
Figure 12. Comparison of maximum overpressures generated by hydrogen tank explosions.
Energies 16 00241 g012
Figure 13. Explosion simulation model in underground parking lot: (a) schematic; (b) mesh generation.
Figure 13. Explosion simulation model in underground parking lot: (a) schematic; (b) mesh generation.
Energies 16 00241 g013
Figure 14. Explosion simulation results at underground parking lot: (a) overpressure according to height (52 L tank); (b) overpressure according to tank volume (1.5 m height).
Figure 14. Explosion simulation results at underground parking lot: (a) overpressure according to height (52 L tank); (b) overpressure according to tank volume (1.5 m height).
Energies 16 00241 g014
Figure 15. Road tunnel simulation model schematic and mesh generation: (a) tunnel area: 54 m2; (b) tunnel area: 72 m2.
Figure 15. Road tunnel simulation model schematic and mesh generation: (a) tunnel area: 54 m2; (b) tunnel area: 72 m2.
Energies 16 00241 g015
Figure 16. Explosion simulation results at road tunnel: (a) tunnel area: 54 m2; (b) tunnel area: 72 m2.
Figure 16. Explosion simulation results at road tunnel: (a) tunnel area: 54 m2; (b) tunnel area: 72 m2.
Energies 16 00241 g016
Figure 17. Limit distances according to tank capacity in the underground parking garage.
Figure 17. Limit distances according to tank capacity in the underground parking garage.
Energies 16 00241 g017
Figure 18. Limit distances in the road tunnel according to cross-sectional area and tank capacity.
Figure 18. Limit distances in the road tunnel according to cross-sectional area and tank capacity.
Energies 16 00241 g018
Figure 19. Limit distances according to space conditions for the explosion of a 52 L hydrogen tank.
Figure 19. Limit distances according to space conditions for the explosion of a 52 L hydrogen tank.
Energies 16 00241 g019
Figure 20. Comparison of overpressure–impulse curves (harm criteria for humans).
Figure 20. Comparison of overpressure–impulse curves (harm criteria for humans).
Energies 16 00241 g020
Table 1. Maximum overpressure and impulse at each measurement point.
Table 1. Maximum overpressure and impulse at each measurement point.
Distance [m]EastWestNorth
Max Overpressure [kPa]Max Impulse [Pa·s]Max Overpressure [kPa]Max Impulse [Pa·s]Max Overpressure [kPa]Max Impulse [Pa·s]
1.91445366.5922673214.082972187.842
4.2811594.4571112019.31621565.239
6.5341640.278731510.218441676.127
Table 2. Comparison of Maximum Overpressure Experiments and Simulations.
Table 2. Comparison of Maximum Overpressure Experiments and Simulations.
Maximum OverpressureRate of Error [%]
Experimental ResultsSimulation Results
Distance 1.9 m1.44 kPa1.72 kPa19
Distance 4.2 m0.81 kPa0.78 kPa4
Distance 6.5 m0.34 kPa0.42 kPa24
Table 3. Capacities, filling pressures, and overpressure conditions of hydrogen tank explosions.
Table 3. Capacities, filling pressures, and overpressure conditions of hydrogen tank explosions.
TestTank Volume [L]
(Mass [kg])
Initial Tank Pressure
[MPa]
Tank Rupture Pressure
[MPa]
Pressure Difference
[MPa]
Ref.
Experiments in this study52.2
(2.0)
708010This study
CFD simulation52.27070-This study
Ref. Test A
(Type IV tank)
357094.5424.54[11,19]
Ref. Test B
(Type III tank)
367099.4729.47[11,19]
Ref. West and North72.4
(1.654)
34.335.71.4[18]
Table 4. Limit distances according to the effects of pressure on the human body.
Table 4. Limit distances according to the effects of pressure on the human body.
Pressure Range [kPa]≥100100–16.516.5–1.35≤1.35
Effect on human bodyFatalitySerious injurySlight injuryNo injury
Type of harm1% (100 kPa)–99% (200 kPa) probability of fatality owing to lung hemorrhage1% (16.5 kPa)–90% (84 kPa) probability of eardrum ruptureIndirect effects (injuries likely from broken glass, structure debris, body translation)-
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Park, J.; Yoo, Y.; Ryu, J.; Lee, H. Study on the Explosion of the Hydrogen Fuel Tank of Fuel Cell Electric Vehicles in Semi-Enclosed Spaces. Energies 2023, 16, 241. https://doi.org/10.3390/en16010241

AMA Style

Park J, Yoo Y, Ryu J, Lee H. Study on the Explosion of the Hydrogen Fuel Tank of Fuel Cell Electric Vehicles in Semi-Enclosed Spaces. Energies. 2023; 16(1):241. https://doi.org/10.3390/en16010241

Chicago/Turabian Style

Park, Jinouk, Yongho Yoo, Jioh Ryu, and Hohyung Lee. 2023. "Study on the Explosion of the Hydrogen Fuel Tank of Fuel Cell Electric Vehicles in Semi-Enclosed Spaces" Energies 16, no. 1: 241. https://doi.org/10.3390/en16010241

APA Style

Park, J., Yoo, Y., Ryu, J., & Lee, H. (2023). Study on the Explosion of the Hydrogen Fuel Tank of Fuel Cell Electric Vehicles in Semi-Enclosed Spaces. Energies, 16(1), 241. https://doi.org/10.3390/en16010241

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop