Research on hydrogen leakage and diffusion mechanism in hydrogenation station | Scientific Reports

Scientific Reports volume 14, Article number: 3363 (2024) Cite this article

1654 Accesses

Metrics details

As a clean, efficient and sustainable energy carrier, hydrogen energy has been accepted as one of the main directions of future energy development. In this paper, a hydrogenation station providing compressed hydrogen outside was adopted as the research object. Based on finite element method and virtual nozzle model, the influence of leakage of main equipment in hydrogenation station on the distribution of combustible hydrogen was investigated, including hydrogen storage tank group, tube trailer, compressor chamber and hydrogenator. The results showed that the shape and volume of the combustible hydrogen cloud generated by the leak were influenced by obstacles, hydrogen storage pressure, and wind velocity. The disturbance of external wind and the decrease in hydrogen storage pressure will have a positive impact on the reduction of leaked volume. The diffusion of combustible hydrogen clouds can be exacerbated by complex structure of obstacles, while partition wall in the adopted hydrogenation station model can limit the combustible hydrogen cloud in the process area. These conclusions can provide guidance and reference for the risk prevention measures of hydrogenation station.

Hydrogen energy, as a clean, efficient, and sustainable energy carrier, has garnered widespread recognition as a pivotal focal point for future energy development1,2,3. Promoting the utilization of hydrogen energy is beneficial for reducing carbon emissions, protecting the environment, which has environmental and social benefits. As the infrastructure to provide hydrogen for hydrogen fuel cell vehicles, the hydrogenation station becomes an indispensable part of hydrogen energy utilization. However, due to the characteristics of hydrogen, such as a low minimum ignition energy (about 0.02 MJ)4, a wide range of explosion limit (4%∼75%)5,6,7, the safety of hydrogenation has been always questioned. These factors can easily lead to combustion and explosion accidents involving hydrogen, which made the development of hydrogen energy system restricted8,9,10. Therefore, the analysis and research of the influencing factors related to hydrogen leakage, combustion, and explosions in hydrogenation stations hold significant theoretical and practical importance11. This research is crucial for enhancing safety and expediting the healthy development of the hydrogen energy industry.

Relevant scholars have carried out numerous theoretical and experimental researches on the hydrogenation station to reveal the general rules of hydrogen leakage and explosion. De Stefano et al.12 studied the impact of leakage location, leakage speed, and obstacles on hydrogen leakage accidents in confined spaces through experimental methods. The results show that the impact of leakage speed on the distribution of hydrogen concentration in confined spaces is greater than the impact of leakage location. The leaked hydrogen encountering obstacles will cause a loss of kinetic energy and increase the hydrogen concentration gradient in the upper area of the space. Gye et al.13 proposed a QRA (Quantitative Risk Analysis) method for quantitative risk assessment of high-pressure hydrogen refueling stations in densely populated and highly congested urban areas.

The finite element method is also regarded as the common method for hydrogen leakage problems14,15,16. Numerous experiments and simulations conducted by scholars, focusing on hydrogen accident simulation through computational fluid dynamics software, have provided both valid data and a theoretical foundation for comprehending the mechanism of hydrogen leakage17,18,19. Kuroki et al.20 carried out a fire thermal radiation risk assessment study on hybrid hydrogen station. Qian et al.21 compared the unintended hydrogen release in both momentum-buoyancy dominated regime and momentum dominated regime. The results show that the profile of the flammable gas cloud goes upward in the momentum-buoyancy dominated regime. The flammable gas cloud is more likely to accumulate towards the ground in a momentum-dominated environment. Li et al.22 have conducted a comprehensive analysis of the leakage and diffusion distribution of natural gas and hydrogen mixtures within a closed container through numerical simulations. While their study provides valuable insights into confined environments, it is imperative to extend the investigation to the leakage process in open space. Exploring the dynamics of gas mixture dispersion in unconfined settings will contribute significantly to a more thorough understanding of the overall behavior and potential risks associated with these substances. Liang et al.23 employed FLACS software and a computational fluid dynamics approach to simulate hydrogen storage system leakage and explosions in a renewable energy hydrogen production station. The consequences of accidents were analyzed in their study, including the harmful area, lethal area, farthest harmful distance, and longest lethal distance, in relation to variables such as wind velocity, leakage direction, and wind direction. However, there is a lack of necessary suggestions for the structural renovation and optimization of hydrogen station.

The innovation of this paper lies in extending above researches to the leakage process in open spaces. Moreover, a set of optimization suggestions for obstacle and ceiling in the hydrogen stations has been proposed. In this work, the leakage of Hydrogen station in an open space was simulated by finite element method. The distribution of the most unfavorable combustible gas cloud was calculated under no-wind conditions through the leakage of various equipment in the hydrogen station. The effects of different external wind speeds, leakage directions, obstacle structures, and ceiling structures were compared, and the mechanism of leakage and diffusion was analyzed. Finally, preventive measures for hydrogen leakage and explosion accidents in the hydrogen station were proposed. Necessary suggestion has been drawn in the conclusion.

The FLACS (flame acceleration simulation)24,25 is a computational fluid dynamics software adopting the finite volume method. The temperature, fuel concentration, combustion products, overpressure, and other variables in the computational domain can be determined by the semi-implicit method of pressure coupled equations combined with boundary conditions. The governing equations was established based on energy conservation, momentum conservation, mass conservation, and component conservation. The interaction between the flame and overpressure generated by explosion and surrounding environment was considered as following equation:

where, φ is the variable of the conservation equations of mass; momentum, and energy; ρ is the density of air, kg/m3; ∂xj represent differential along the j direction, m; ui is the velocity along the i direction, m/s; Γφ is the diffusion coefficient; Sφ is the source term.

The concentration of gas mixture is determined by the global equivalence ratio θ, where the θact represents the actual mass ratio of fuel to oxygen, kg/kg; the θsto represents the stoichiometric mass ratio of fuel to oxygen, kg/kg; moxygen is the mass of oxygen in the gas diffusion zone, kg. The detailed calculation is shown in Eq. (2).

A series of excitation force induced by fluid fluctuations at leakage location of high-pressure hydrogen storage tanks, which make the fluid parameters for numerical calculations difficult to converge. Therefore, the virtual nozzle model26,27 was used to simplify the study of high-pressure hydrogen storage vessel leakage. The equivalent outlet was substituted for leakage outlet. It should be noted that there is not such a location in the actual leakage process, which is only an assumption made to meet the needs of the model. The calculation principle is shown in Fig. 1, where P1, P2, P3 represent the pressure of interior medium in the hydrogen storage container, leakage location, and the equivalent outlet, Pa; T1, T2, T3 represent the temperature of medium, K; ρ1, ρ2, ρ3 represent the density of medium, kg/m3. The pressure of the jet-flow at the equivalent outlet is the environmental pressure. According to the conservation of mass, the flow rate at the equivalent outlet is equal to the flow rate at the actual leakage location. Other flow parameters can be calculated based on the virtual nozzle model.

Virtual nozzle model.

As depicted in Fig. 1, the mass flow rate at the leakage location Q2 can be expressed as Eq. (3).

where V2 local acoustic velocity of hydrogen gas at the actual leakage location, m/s. The mass flow rate at the equivalent outlet Q3 can be expressed as Eq. (4).

where V3 surrounding environmental local acoustic velocity of hydrogen gas, m/s. According to the continuity equation, the relationship between the parameters at two location is as follows.

Due to the instantaneous completion of the flow process from hydrogen storage container to actual leakage location, the flow can be regarded as adiabatic flow. Therefore, based on the isentropic relationship equation, the fluid parameters at actual leakage location can be obtained.

where P2 is the pressure of actual leakage location, Pa; T2 is the thermodynamic temperature of actual leakage location, K; γ is the specific heat ratio, taken as 1.40. Other fluid parameters were derived from ideal gas state equation.

where M is the molar mass of hydrogen; R represents the universal gas constant, taken as 8.314 J·mol−1·K−1.

According to the proportional relationship between local sound velocity and the square root of thermodynamic temperature, the velocity of leakage location and equivalent outlet can be expressed as follow:

The length relationship between the actual leakage diameter and the equivalent outlet diameter can be obtained by substituting Eqs. (9) and (10) into Eq. (5).

To investigate the impact of different leak locations, leak directions, leakage pressure, and external wind speed on the consequences of accidents, high-pressure hydrogen leakage simulation was carried out. A 1:1 physical model of a hydrogen station was simplified and established in pre-processing module CASD, as shown in Fig. 2. The hydrogen station is generally in the form of a right triangle. The three-dimensional coordinate system of hydrogen station was established centered on the southwest corner. The east direction is the positive direction of the X axis, the north direction is the positive direction of the Y axis, and the vertical direction is the positive direction of the Z axis.

Plane layout of hydrogenation station created by FLACS [version: 9.0; URL: https://www.gexcon.com/support/flacs-cfd/downloads/].

Approximately 90 m from east to west and 60 m from south to north, the west and south sides are occupied by certain public buildings, and the north and south sides by a highway. The station equipment mainly includes 2 hydrogen tube trailers, 2 unloading columns, 2 compressors, 2 fixed gas storage cylinder sets, 4 hydrogenation machines, station buildings, shelters etc.

The pressure of seamless steel gas cylinder during transportation in tube trailer is 20 MPa, while the exhaust pressure is 5 MPa. The volume of gas cylinder is 26 m3. The type of hydrogen refueling machine is dual gun, dual metering, and three-line, with a filling pressure of 35 MPa. The compressor model is PDC-13-7500 (100). The specific parameters are shown in Table 1.

The parameters of the hydrogen storage bottle group are listed in Table 2.

According to the formulas (2), (3), (4), (5), (6), (7), (8), (9), (10) and (11), the boundary conditions at the leakage location of the equipment in the hydrogen station can be found in Table 3. The values of temperature were all set to 15 °C, which is determined by ambient temperature.

To simulate the diffusion of hydrogen leakage under the most unfavorable conditions, the pressure at which hydrogen leakage occurred in the hydrogen storage cylinder group was taken as 45 MPa. FLACS v9.0 can only simulate leakage diffusion at a constant leakage rate during simulation, and cannot fully simulate the leakage process of pressure vessels. Therefore, during the simulation setup process, it is assumed that hydrogen gas leaks at a constant velocity calculated at the highest pressure.

The leakages of combustible gas under different environmental wind speeds and leakage locations for 5 s were shown in Fig. 3. Distribution of combustible gas leaked from hydrogen storage cylinder group in + Z, + X, -Y, and -Z direction under no wind conditions were depicted as Fig. 3a–d respectively. Under external wind conditions of 5 m/s, Fig. 3e depicted the distribution of combustible gas l that escaped from the hydrogen storage cylinder group in the -Z direction. Under external wind conditions of 10 m/s, distribution of combustible gas was shown as Fig. 3f.

Leakage of hydrogen at the hydrogen storage cylinder group (t = 5 s) created by FLACS [version: 9.0; URL: https://www.gexcon.com/support/flacs-cfd/downloads/].

When the leakage was along the + Z direction, the leakage of hydrogen was perpendicular to the ground and upward, with no obstacles or external wind disturbance. It can be found that the combustible hydrogen cloud was in the shape of a vertical and elongated "rugby", as depicted in Fig. 3a. The possibility of an explosion occurring when encountering an ignition source was close to zero since the lack of contact between the leaked gas and any facilities in the station.

When the leakage was along the + X direction, the leakage of hydrogen was squeezed towards the partition wall by high-pressure, as shown in Fig. 3b. The building in the direction of leakage was a partition wall with a simple structure. The hydrogen diffused to all sides as the leaking gas hit the partition wall, accompanied by a loss of kinetic energy during the process. Due to the low density and rapid diffusion of hydrogen, it rapidly diluted upwards, resulting in a smaller volume of combustible gas clouds.

When the leakage was along the -Y direction, the hydrogen leaked from the hydrogen storage bottle group was blocked by another hydrogen storage bottle group, as illustrated in Fig. 3c. After reaching another hydrogen storage bottle group, the high-pressure jet was diffused between the gaps of the hydrogen storage bottles, ultimately forming a large area of irregular combustible hydrogen cloud. However, due to the obstruction of the partition wall, the combustible hydrogen cloud was restricted within the process area.

When the leakage was along the -Z direction, the leakage of hydrogen was perpendicular to the ground and downward, as indicated in Fig. 3d. Under the influence of air buoyancy, hydrogen dispersed across the ground and to nearby facilities or buildings, generating irregular combustible hydrogen cloud. Due to the obstruction of the partition wall, combustible hydrogen clouds were similarly confined within the process area.

It can be found that the influence of hydrogen leakage direction on the shape and volume of combustible hydrogen clouds was mainly determined by the structure and density of obstacles in the leakage direction. If there were no obstacles in the leakage direction or if the structure of the obstacles was regular and simple, the leaked combustible hydrogen cloud had a regular shape and small volume. The probability of explosion when encountering an ignition source is relatively low. On the contrary, if obstacles were densely congested in the leakage direction, the leaked combustible hydrogen cloud had an irregular shape and larger volume, making it more likely to cause an explosion accident.

As Fig. 3e and findicated, the distribution of combustible hydrogen cloud varied from the influence of different velocity of external wind. As the external wind velocity increased, the volume of combustible gas gradually decreased. It can be concluded that the external wind can accelerate the diffusion and dilution of hydrogen in the air. The higher the wind velocity, the stronger the turbulence effect of the atmosphere. The safety of the hydrogen station was guaranteed due to the smaller volume of leaked combustible gas. Even in the presence of external wind, the combustible hydrogen cloud can still be contained within the process area of the hydrogen station, preventing the leaked hydrogen from spreading to other regions.

To further investigated the diffusion of hydrogen leakage under the most unfavorable conditions, the influence of external wind on the other key facilities was neglected. The velocity of external wind was set to zero.

The leakages of tube trailer in different directions for 5 s were indicated in Fig. 4. Distributions of combustible gas leaked from tube trailer in + Z, + X, and − Z direction under no wind conditions were depicted as Fig. 4a–c.

Leakage of hydrogen at tube trailer. (t = 5 s) created by FLACS [version: 9.0; URL: https://www.gexcon.com/support/flacs-cfd/downloads/].

The pressure of hydrogen leakage from the tube trailer was 20 MPa, while from the hydrogen storage tank group is 45 MPa. It can be found that the shapes of combustible hydrogen clouds leaked from tube trailer and hydrogen storage tank group in the same direction were similar. Due to the lower pressure inside the tube trailer, the volume of combustible hydrogen cloud is smaller. Therefore, probability of explosion accidents and damage level caused by explosion accidents were lower.

The physical model of compressor chamber was shown in Fig. 5a. A ventilation vent with a length of 0.3 m was installed above the compressor chamber. In order to simplify the model and computational burden, internal devices such as compressors, pipelines, and valves were ignored. The leakage of compressor chamber in + Z direction for 5 s was indicated in Fig. 5b.

Leakage of hydrogen at compressor chamber. (t = 5 s) created by FLACS [version: 9.0; URL: https://www.gexcon.com/support/flacs-cfd/downloads/].

The compressor chamber was a confined space. The leaked hydrogen was initially restricted by the wall before diffusing into the surrounding environment via the vents. As shown in Fig. 5b, the hydrogen escaping from the vent diffused around the compressor chamber. The combustible hydrogen cloud spread to the areas of circulating pumps, tube trailer, and hydrogen storage bottle group. The shape of combustible hydrogen cloud was irregular. Despite its large volume, the combustible hydrogen cloud was still limited by the partition walls in the process area.

Sectional views of hydrogenator with flat roof shelter and concave roof shelter were depicted in Fig. 6a and b respectively. The leakage direction was + Z, the leakage pressure was 35 MPa, and the leakage time was 5 s.

Sectional view of the shelter in the hydrogenation machine area.

The distributions of hydrogen concentration under different roof shelter at Y = 16.3 m plane were shown in Fig. 7. Areas with hydrogen concentrations over the lower explosive limit (4%) were marked in red. The shapes of red area in the two types of shelter were similar, which both were beneath the shelter and above the hydrogenator. Therefore, the probability of explosion occurring when encountering an ignition source is relatively low.

The distributions of hydrogen concentration under different roof shelter at Y = 16.3 m plane.

The high-pressure hydrogen jet leaking in the + Z direction diffused along the bottom of the hydrogenator shelter. Along the flat roof shelter, hydrogen diffused into the atmosphere via bottom surface of shelter and gradually diluted upwards. Along the concave roof shelter, hydrogen was obstructed by the vertical planes around the shelter, causing it to diffuse to the area between the ground and the shelter. The hydrogen concentration inside the area was around 2%. The concentration was below the lower explosion limit. However, if there was a problem of filling vehicles or sundries storage in the hydrogenator, the hydrogen concentration in the area will be increased by the poor ventilation. In addition, long-term leakage may also cause the local concentration of hydrogen to reach the lower explosion limit. As a result, using a flat roof shelter in a hydrogenator is safer than using a concave roof shelter.

In this paper, the diffusion of hydrogen leakage under the most unfavorable conditions was simulated by finite element analysis method. The main conclusions are drawn as follows:

The distribution shape of leaked hydrogen was mainly determined by the structure and density of obstacles in the leakage direction. The dense blockage of obstacles in the direction of leakage can cause the diffusion of combustible hydrogen clouds to be irregular in shape, larger in volume, and more likely to cause explosion accidents. Optimizing equipment layout, minimizing obstacles in the direction of leakage should be considered to ensure smooth ventilation systems

The device with higher hydrogen storage pressure produces a larger volume of combustible hydrogen gas cloud in the event of a leakage. Safety measures should be categorized based on the hydrogen storage pressure of different devices.

The turbulence effect caused by external wind has a positive effect on the dilution of leaked gas. As the external wind velocity increased, the volume of combustible gas gradually decreased.

Prioritize choosing a flat roof shelter over a concave roof shelter. The combustible gas diffusion performance of flat roof shelter is better than that of concave roof shelter. A large amount of combustible gas will accumulate inside the concave roof shelter, which will cause potential safety hazards to the hydrogenator.

The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Rong, Y. et al. Numerical investigation on the liquid hydrogen leakage and protection strategy. Process 11, 1173. https://doi.org/10.3390/PR11041173 (2023).

Article CAS Google Scholar

Sakamoto, J. et al. Leakage-type-based analysis of accidents involving hydrogen fueling stations in Japan and USA. Int. J. Hydrogen Energy 41, 21564–21570. https://doi.org/10.1016/J.IJHYDENE.2016.08.060 (2016).

Article ADS CAS Google Scholar

Chen, Q. et al. Analysis of unsteady heat transfer in the pre-cooling process of 300 M3 liquid hydrogen tank. Process. 11, 787. https://doi.org/10.3390/PR11030787 (2023).

Article CAS Google Scholar

Li, Y., Wang, Z., Shi, X. & Fan, R. Safety analysis of hydrogen leakage accident with a mobile hydrogen refueling station. Process Saf. Environ. Protect. 171, 619–629. https://doi.org/10.1016/J.PSEP.2023.01.051 (2023).

Article CAS Google Scholar

Cao, X. et al. Effect of obstacle on the H2/CO/air explosion characteristics under lean-fuel conditions. Fuel 319, 123834. https://doi.org/10.1016/J.FUEL.2022.123834 (2022).

Article CAS Google Scholar

Hansen, O. R. Hydrogen infrastructure—Efficient risk assessment and design optimization approach to ensure safe and practical solutions. Process Saf. Environ. Protect. 143, 164–176. https://doi.org/10.1016/J.PSEP.2020.06.028 (2020).

Article CAS Google Scholar

Astbury, G. R. A review of the properties and hazards of some alternative fuels. Process Saf. Environ. Protect. 86, 397–414. https://doi.org/10.1016/J.PSEP.2008.05.001 (2008).

Article CAS Google Scholar

Chang, X., Bai, C. & Zhang, B. The effect of gas jets on the explosion dynamics of hydrogen-air mixtures. Process Saf. Environ. Protect. 162, 384–394. https://doi.org/10.1016/J.PSEP.2022.04.032 (2022).

Article CAS Google Scholar

Machniewski, P. & Molga, E. CFD analysis of large-scale hydrogen detonation and blast wave overpressure in partially confined spaces. Process Saf. Environ. Prot. 158, 537–546. https://doi.org/10.1016/J.PSEP.2021.12.032 (2022).

Article CAS Google Scholar

Rui, S., Wang, C., Guo, S., Jing, R. & Li, Q. Hydrogen-air explosion with concentration gradients in a cubic enclosure. Process Saf. Environ. Protect. 151, 141–150. https://doi.org/10.1016/J.PSEP.2021.05.003 (2021).

Article CAS Google Scholar

Li, H. et al. Safety of hydrogen storage and transportation: An overview on mechanisms, techniques, and challenges. Energy Reports 8, 6258–6269. https://doi.org/10.1016/J.EGYR.2022.04.067 (2022).

Article Google Scholar

De Stefano, M., Rocourt, X., Sochet, I. & Daudey, N. Hydrogen dispersion in a closed environment. Int. J. Hydrogen Energy 44, 9031–9040. https://doi.org/10.1016/J.IJHYDENE.2018.06.099 (2019).

Article ADS Google Scholar

Gye, H. R., Seo, S. K., Bach, Q. V., Ha, D. & Lee, C. J. Quantitative risk assessment of an urban hydrogen refueling station. Int. J. Hydrogen Energy 44, 1288–1298. https://doi.org/10.1016/J.IJHYDENE.2018.11.035 (2019).

Article ADS CAS Google Scholar

Gong, L. et al. Experimental investigation on the dispersion characteristics and concentration distribution of unignited low-temperature hydrogen release. Process Saf. Environ. Protect. 160, 676–682. https://doi.org/10.1016/J.PSEP.2022.02.055 (2022).

Article CAS Google Scholar

Lee, J. et al. CFD modeling on natural and forced ventilation during hydrogen leaks in a pressure regulator process of a residential area. Process Saf. Environ. Prot. 161, 436–446. https://doi.org/10.1016/J.PSEP.2022.03.065 (2022).

Article CAS Google Scholar

Wang, Q. et al. Analytical and numerical predictions of hydrogen gas flow induced by wall and corner leakages in confined space. Int. J. Hydrogen Energy 45, 6848–6862. https://doi.org/10.1016/J.IJHYDENE.2019.12.219 (2020).

Article ADS CAS Google Scholar

Kim, E., Park, J., Cho, J. H. & Moon, I. Simulation of hydrogen leak and explosion for the safety design of hydrogen fueling station in Korea. Int. J. Hydrogen Energy 38, 1737–1743. https://doi.org/10.1016/J.IJHYDENE.2012.08.079 (2013).

Article ADS CAS Google Scholar

Kim, E. et al. Development of Korean hydrogen fueling station codes through risk analysis. Int. J. Hydrogen Energy 36, 13122–13131. https://doi.org/10.1016/J.IJHYDENE.2011.07.053 (2011).

Article ADS CAS Google Scholar

Bauwens, C. R. & Dorofeev, S. B. CFD modeling and consequence analysis of an accidental hydrogen release in a large scale facility. Int. J. Hydrogen Energy 39, 20447–20454. https://doi.org/10.1016/J.IJHYDENE.2014.04.142 (2014).

Article ADS CAS Google Scholar

Kuroki, T., Sakoda, N., Shinzato, K., Monde, M. & Takata, Y. Temperature rise of hydrogen storage cylinders by thermal radiation from fire at hydrogen-gasoline hybrid refueling stations. Int. J. Hydrogen Energy 43, 2531–2539. https://doi.org/10.1016/J.IJHYDENE.2017.12.072 (2018).

Article ADS CAS Google Scholar

Qian, J., Li, X., Gao, Z. & Jin, Z. A Numerical study of unintended hydrogen release in a hydrogen refueling station. Int. J. Hydrogen Energy 45, 20142–20152. https://doi.org/10.1016/J.IJHYDENE.2020.05.063 (2020).

Article ADS CAS Google Scholar

Li, H. et al. Numerical simulation of leakage and diffusion distribution of natural gas and hydrogen mixtures in a closed container. Int. J. Hydrogen Energy 47, 35928–35939. https://doi.org/10.1016/J.IJHYDENE.2022.08.142 (2022).

Article ADS CAS Google Scholar

Liang, Y., Pan, X., Zhang, C., Xie, B. & Liu, S. The simulation and analysis of leakage and explosion at a renewable hydrogen refuelling station. Int. J. Hydrogen Energy 44, 22608–22619. https://doi.org/10.1016/J.IJHYDENE.2019.05.140 (2019).

Article ADS CAS Google Scholar

Lucas, M., Atanga, G., Hisken, H., Mauri, L. & Skjold, T. Simulating vented hydrogen deflagrations: Improved modelling in the CFD Tool FLACS-hydrogen. Int. J. Hydrogen Energy 46, 12464–12473. https://doi.org/10.1016/J.IJHYDENE.2020.09.073 (2021).

Article ADS CAS Google Scholar

Moen, A., Mauri, L. & Narasimhamurthy, V. D. Comparison of K-ε models in gaseous release and dispersion simulations using the CFD code FLACS. Process Saf. Environ. Prot. 130, 306–316. https://doi.org/10.1016/J.PSEP.2019.08.016 (2019).

Article CAS Google Scholar

Yüceil, K. B. & Ötügen, M. V. Scaling parameters for underexpanded supersonic jets. Phys. Fluids 14, 4206–4215. https://doi.org/10.1063/1.1513796 (2002).

Article ADS CAS Google Scholar

Birch, A. D., Brown, D. R., Dodson, M. G. & Swaffield, F. The structure and concentration decay of high pressure jets of natural gas. Combust. Sci. Technol. 36, 249–261. https://doi.org/10.1080/00102208408923739 (2007).

Article Google Scholar

Download references

School of Thermal Engineering, Shandong Jianzhu University, Jinan, 250101, Shandong, China

Zenggang Zhang

Shandong Provincial Architecture Design & Research Institute CO., Ltd., Jinan, 250101, Shandong, China

Mingheng Shang

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

Z.Z. and M.S. wrote the main manuscript text and prepared figures. All authors reviewed the manuscript.

Correspondence to Zenggang Zhang.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

Zhang, Z., Shang, M. Research on hydrogen leakage and diffusion mechanism in hydrogenation station. Sci Rep 14, 3363 (2024). https://doi.org/10.1038/s41598-023-50917-4

Download citation

Received: 11 July 2023

Accepted: 28 December 2023

Published: 09 February 2024

DOI: https://doi.org/10.1038/s41598-023-50917-4

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative