Next Article in Journal
Gelatin–Siloxane Hybrid Monoliths as Novel Heavy Metal Adsorbents
Next Article in Special Issue
Analysis of Hydrogen Filling of 175 Liter Tank for Large-Sized Hydrogen Vehicle
Previous Article in Journal
Borehole Geophysical Joint Inversion to Fully Evaluate Shaly Sandstone Formations
Previous Article in Special Issue
Analysis of Bluetooth RSSI for Proximity Detection of Ship Passengers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Analysis of the Effects of Ship Motion on Hydrogen Release and Dispersion in an Enclosed Area

1
Department of Refrigeration and Air-Conditioning Engineering, Korea Maritime and Ocean University, 727, Taejong-ro, Yeongdo-gu, Busan 49112, Korea
2
Interdisciplinary Major of Maritime AI Convergence, Korea Maritime and Ocean University, 727, Taejong-ro, Yeongdo-gu, Busan 49112, Korea
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(3), 1259; https://doi.org/10.3390/app12031259
Submission received: 23 December 2021 / Revised: 18 January 2022 / Accepted: 19 January 2022 / Published: 25 January 2022

Abstract

:
Hydrogen is an alternative to conventional heavy marine fuel oil following the initial strategy of the International Maritime Organization (IMO) for reducing greenhouse gas emissions. Although hydrogen energy has many advantages (zero-emission, high efficiency, and low noise), it has considerable fire and explosion risks due to its thermal and chemical characteristics (wide flammable concentration range and low ignition energy). Thus, safety is a key concern related to the use of hydrogen. Whereas most previous studies focused on the terrestrial environment, we aim to analyze the effects of the ship’s motion on hydrogen dispersion (using commercial FLUENT code) in an enclosed area. When compared to the steady state, our results revealed that hydrogen reached specific sensors in 63% and 52% less time depending on vessel motion type and direction. Since ships carry and use a large amount of hydrogen as a power source, the risk of hydrogen leakage from collision or damage necessitates studying the correspondence between leakage, diffusion, and motion characteristics of the ship to position the sensor correctly.

1. Introduction

Increasing concerns regarding climate change attributed to air pollution have led to the establishment of stringent regulations for reducing greenhouse gas (GHG) emissions and air pollutants, which are one of the main causes of air pollution [1]. The shipping industry transports about 90% of world trade. The GHGs and pollutants comprising CO2, CH4, SOx, and NOx are generated by the process of fuel consumption in ships [2,3,4]. Based on a IMO GHG study, marine shipping emitted 1056 million tons of CO2 in 2018, representing 2.89% of the world’s total emission. These emissions are predicted to increase by about 90–130% from the 2008 emissions (921 million tons) by 2050 if no action is initiated [5]. Furthermore, at the 72nd Meeting (2018) of the Marine Environment Protection Committee (MEPC), the International Maritime Organization (IMO) adopted the Initial Strategy for IMO GHG Reduction. In accordance with this strategy, the volume of the total annual GHG emissions from international shipping must reduce by at least 50% by 2050, compared with 2008 [6,7]. Based on these strategies, Det Norske Veritas (DNV) announced that the most innovative approach for reducing carbon dioxide emissions is to change marine fuels to eco-friendly fuels [8]. In Lloyd’s Register, hydrogen, biofuel, ammonia, methanol, and electricity are selected as fuel alternatives to reduce CO2; the advantages and disadvantages of each fuel are compared, and the price change is predicted by year. Thus, hydrogen energy was evaluated to be superior in cost and reduced air pollutants, fine dust, and greenhouse gases [9]. Although hydrogen energy is one of the most promising alternative fuels for ships in the future, it has a low ignition energy (0.02 mJ) and a wide combustibility range (4–75%); therefore, there is a high possibility of fire or explosion accidents in the case of a hydrogen leak [10]. Therefore, leakage management, detection, and ventilation systems are important for increasing safety during fires and explosions while using hydrogen energy. In this regard, the IMO adopted resolution MSC.391(95), the International Code of Safety for Ships using Gases or Low-flashpoint Fuels (IGF Code) such as hydrogen [11]. In accordance with this resolution, the IMO recommends using gas dispersal analysis or a physical smoke test for leak detection equipment best arrangement. Therefore, for hydrogen leakage detection, it is necessary to quantitatively analyze hydrogen dispersion characteristics under various parameter conditions (release pressure, velocity, diameter, etc.) for installing a detector based on these characteristics.
Table 1 summarizes various types of hydrogen release scenarios considered in the experimental or numerical methods, including gaseous releases, fuel cell ships, vehicles, semi-closed, and confined environments. Li et al. [12] numerically studied the release of hydrogen in a hydrogen fuel cell ship and analyzed the effects of different ventilation conditions on hydrogen dispersion behavior. Their results indicated the optimal positions for hydrogen sensors and ventilation. Liu et al. [13] conducted a CFD simulation on high-pressure hydrogen released from storage tanks; furthermore, factors such as wind speed, ambient temperature, leak location, leak hole diameter, and obstructions on hydrogen diffusion were evaluated. Dadashzadeh et al. [14] summarized their results on the dispersion behavior of hydrogen gas from a hydrogen fuel cell car in an enclosed area and provided safety measures for using hydrogen in ventilation conditions. Malakhov et al. [15] compared experimental data with CFD computed results for a hydrogen release in a semi-closed ventilation facility and assessed the efficiency of forced ventilation. A simulation was performed using STAR-CCM+. Lacome et al. [16] performed large-scale hydrogen release experiments to analyze the formation of flammable clouds inside a confined area with an internal volume of 80 m3. Earlier studies related to hydrogen leakage and diffusion focused on terrestrial environments.
However, the safety risks and accident levels are considerably higher than those of hydrogen vehicles on land because many electrical facilities can act as ignition sources for hydrogen leaks in hydrogen ships, and the hydrogen storage tank’s capacity is large. The ship has an airtight structure; therefore, hydrogen cannot be easily discharged outside. The risk is considerably greater at a hydrogen refueling station when an accident such as an explosion occurs [12,17]. Despite these risks, few studies have focused on hydrogen diffusion behavior for hydrogen ships. Previous studies only considered the structural conditions of the ship, whereas environmental factors, such as the motion conditions of the ship, have not been published. We analyzed the effects of the ship’s kinetic characteristics on the hydrogen diffusion behavior and conducted safety analysis using simulations attributed to environmental factors of the ship. We detected the time required by hydrogen to reach each sensor due to the ship’s motion conditions. Furthermore, we indicated the need for selecting an appropriate location for installing a detector by considering kinetic characteristics such as motion types, directions, and periods.
The rest of this manuscript is organized as follows: Section 2 presents the methodology for analyzing hydrogen diffusion behavior while considering the ship’s motion. Section 3 provides definitions of the ship’s motion caused by hydrogen diffusion based on the ship motion and hydrogen diffusion characteristics for scenarios including ship motion type, period, and direction. Finally, Section 4 summarizes the main results of our study and reviews its limitations.

2. Methodology

We conducted our study based on the following procedure to analyze the effect of ship motion on hydrogen dispersion behavior:
(1)
A hydrogen leak test was performed in a steady state, and the effectiveness of the numerical analysis tool was verified by comparing the experimental and numerical analyses results;
(2)
The effect of ship motion on hydrogen diffusion was analyzed by applying the ship motion scenario to the CFD simulation.
The overall methodology of the current study, including the verification of the numerical analysis tool and analysis of hydrogen dispersion behavior applied to ship motion, are shown in Figure 1.

2.1. CFD Modeling

The CFD code FLUENT19.2 was used as a numerical tool to evaluate the effects of the ship’s motion on diffusion in the event of a hydrogen leak. The governing equations employed to numerically simulate the hydrogen release and dispersion include the continuity, momentum, energy, species transport, turbulent kinetic energy, and dissipation rate equations of turbulent kinetic energy [18].
The conservation of mass is provided by:
ρ t + · ρ v = 0 ,
where ρ , t , and v represent the density, time, and overall velocity vector, respectively.
The conservation of momentum is expressed as:
t ρ v + · ρ v v = p + · τ = + ρ g + F ,  
where p , τ = , g , and F represent the pressure, stress tensor, gravitational acceleration, and force vector, respectively.
The conservation of energy is provided by:
t ρ E + · v ρ E + p = · k e f f T j h j J j + τ e f f = · v ) ,
where E , k e f f , T , h j , and J j represent the total energy, effective conductivity, temperature, sensible enthalpy of species j, and diffusion flux of species j, respectively.
The species transport equation is:
t ρ Y i + · ρ v Y i = · J i + R i + S i ,
where Y i , R i , and S i denote the mass fraction of each species, net rate of the product of species i by chemical reaction, and the rate of creation by additions from the dispersed phase plus any user-defined sources, respectively.

2.2. Model Validation

2.2.1. Experimental Description

The experimental space was a rectangular box (1.0 m long × 0.5 m wide × 0.75 m high). A vent (diameter = 0.02 m) was located toward the left side of the ceiling to ensure constant pressure conditions. The internal pressure of the hydrogen gas was 120 × 105 Pa, whereas the pressure value after the regulator leaked to 14 × 105 Pa. Hydrogen was released at a mass flow rate of 0.965 × 10−5 kg/s through a 1/16-inch (0.001587 m) diameter tube from the right side of the enclosure and directed horizontally. The experiment was conducted for 240 s to reach a steady-state hydrogen distribution inside the enclosure. A detailed description of the experiment is shown in Figure 2 and Table 2. Seven hydrogen sensors were employed to measure the hydrogen concentration, as shown in Figure 3.

2.2.2. Grid Independence Test

A tetrahedral grid was generated for a shape on the same scale using ANSYS Meshing for the verification experiment. The grid near the hydrogen leakage area and the natural ventilation opening were locally dense to balance the accuracy of the numerical analysis and computing time efficiency. A grid independence test was performed according to the number of grids to determine the numerical validity of the generated grid. Under the aforementioned conditions, an unsteady analysis was performed on six grids from approximately 370,000 to 920,000 grids based on the number of grids. The hydrogen concentration of the ceiling height at a total calculation time of 60 s was compared to examine the effect of the number of grids. The hydrogen gas concentration at the ceiling height from this comparison was used to predict the main flow and development of the concentration field based on the analysis of 492,406 grids or more. According to our results, a grid independence test was conducted for 492,406 grids, as shown in Figure 4.

2.2.3. Numerical Details

The computational domain has the same size as the actual experiment facility (1.0 m long × 0.5 m wide × 0.75 m high) and includes the hydrogen sensor. Boundary conditions such as the leakage flow rate and pressure conditions were applied in the same manner as the experimental conditions. After the geometric model was established, a tetrahedral mesh was generated with a total of 492,406 elements. For the hydrogen leakage and diffusion analysis, the internal flow was assumed to be a three-dimensional, unsteady, incompressible, and turbulent flow. Furthermore, the SIMPLEC algorithm was employed for calculating the pressure and velocity [19,20]. The overall simulation time was 240 s.

2.2.4. Results

Figure 5 shows a comparison between the simulation and experimental results. According to the obtained results, the S1 sensor closest to the hydrogen leak hole responds 6 s after the leak and maintains a concentration of 0.1% thereafter. The sensor responds sequentially as hydrogen leaked from the experimental facility sends discharge through the ventilation opening. In accordance with the sensor position, sensor S4 responded at 33 s after the leak, whereas S6 responded at 73 s after the leak. Therefore, the response time varies based on the sensor position.
In addition, we confirmed that hydrogen moves along S1 through S4 after the leakage because of the low density of hydrogen’s physical characteristics. By comparing the overall hydrogen concentration results, we confirmed that the experimental and numerical analysis results predicted similar concentrations with a difference of approximately 15%.

3. Effect of Ship Motion on Hydrogen Dispersion

3.1. Ship Motion Coordinate System

A ship in contact with the sea level and moving with six degrees of motion is shown in Figure 6. As indicated by Table 3, the coordinate system applied to the calculation has a positive x-axis in the bow direction, a positive y-axis in the starboard direction, and a positive z-axis in the opposite direction of gravity. Here, dynamic motions with periodicity are the roll, pitch, and yaw.

3.2. Ship Motion Scenario

Scenarios were selected using the ship motion type, period, and motion direction as variables for analyzing the effects of ship motion on hydrogen dispersion behavior. Since the motion type of the ship was limited to roll and pitch motion, the motion direction was also divided into clockwise and counterclockwise by considering only the x and y axes. Furthermore, the hydrogen concentration for each sensor location was investigated.
Firstly, only the roll and pitch motions were considered for the ship. Secondly, the ship motion period was calculated for a 15,000 gross tonnage ship according to the following equation provided by IMO’s intact stability code. The rolling period was calculated as 15 s and defined as [21]:
T = 2 C B G M
where T , L , B , d , C , and G M represent the rolling period, length of the ship at the waterline (m), molded breadth of the ship (m), mean molded draught of the ship (m), block coefficient (-), and metacentric height corrected for the free surface effect (m), respectively. Furthermore, C = 0.373 + 0.023 B d 0.043 L 100 .
In addition, simulations were performed by applying 60-s and 120-s cycles to confirm the effects of ship motion periods on the hydrogen concentration value. For the roll and pitch motions, the roll and heel angles were set to move from a minimum of −10° to a maximum of 10° [22]. Finally, the scenario was constructed by classifying ship motions in clockwise and counterclockwise directions. When the motion of the vessel was positive, it was expressed as starboard down and pitching to bow up. When the motion of the vessel was negative, it was expressed as starboard up and pitching down. Figure 7 and Figure 8 show the plane view and corresponding motion period graph, respectively, when the ship motion direction is positive and negative.

3.3. Effect of Ship Motion Types

A numerical analysis was performed by applying roll and pitch motions to the enclosed area; our results revealed the effects of vessel motion type on hydrogen diffusion behavior. Table 4 compares the effects of ship motion type on the hydrogen response time for each sensor S2–S6. The percentages, shown below the time, represent differences in the hydrogen sensor’s response time in the motion compared to the steady state. Compared to the steady state, hydrogen reaches the sensor faster to S4 by approximately 63% (25 s) in pitch motion and S2 by approximately 50% (30 s) in roll motion.
Figure 9 and Figure 10 illustrate the hydrogen concentration values of sensors S4 and S6 located at the center of the ceiling’s leak space for roll and pitch motions with a period of 60 s; they move clockwise on each positive axis. Compared to the roll and pitch motions of S4 and S6 in the steady state, hydrogen reaches the sensor faster by 10 s, 25 s, 34 s, and 64 s, respectively. These results indicate that ship motion type significantly affects hydrogen diffusion.
Compared to the roll motion in Figure 9, the hydrogen concentration increases according to the 60-s cycle after 255 s of hydrogen leakage in the pitch motion. This trend can be explained as follows: After leaking from the hydrogen leak hole at the center of the bow’s side, the hydrogen rises to the ceiling due to the buoyancy attributed to the density difference and forms a ceiling jet flow in the longitudinal direction of the leak space.
Furthermore, it moves between the bow and stern under the effects of pitch motion, which is a longitudinal motion centered on S4; the concentration repeatedly rises and falls according to the cycle. In the case of roll motion, moving between the starboard and port compared with pitch motion, the change in position based on S4 is half, and the concentration diffusion speed is fast; therefore, the effect on the concentration value compared with pitch motion is relatively minor. This phenomenon is confirmed in Figure 10, which shows the hydrogen concentration near sensor S6. For roll motion, the hydrogen concentration rises repeatedly and falls in two cycles (30 s) during one pitch motion cycle (60 s). The hydrogen concentration ranges from 0.0634% at 75 s to 0.0374% at 120 s to 41% after hydrogen leakage in the roll motion, and from 0.0664% at 90 s to 0.0318% at 135 s in the pitch motion with a drop rate of 52%. A downward trend was observed; this tendency is attributed to the rotational force generated in the opposite direction of the leakage before the leaked hydrogen diffuses to the ceiling and forms a concentration layer. The latter occurs because hydrogen moves in the port direction in the case of roll motion.

3.4. Effect of Ship Motion Directions

The numerical analysis was performed by applying roll and pitch motions, which move clockwise and counterclockwise on each axis, to the enclosed area for understanding the effect of the ship’s motion direction on hydrogen diffusion behavior. The hydrogen response time for each sensor, S2–S6, is listed in Table 5 to compare the effects of motion direction. Compared to the steady state, hydrogen reaches the S6 sensor faster by approximately 52% (64 s) during pitch motion and 30% (37 s) during roll motion. Figure 11 shows the central longitudinal cross-section of hydrogen concentration distribution contour based on the time required to reach sensor S6. As shown in Figure 11 and Figure 12, the distribution of hydrogen differs according to the direction of roll and pitch motion.
For pitch motion with a period of 15 s, simulations were performed by dividing motions into clockwise and counterclockwise directions. S2 was close to the leak hole, and sensor S6 was located near the ceiling opening of the enclosed area. Hydrogen concentration values are illustrated in Figure 13 and Figure 14. The leakage port in the center of the stern is located higher than the existing position when moving in the same direction as the hydrogen leakage direction (clockwise). The hydrogen leakage attributed to the updraft reaches the ceiling faster. When S6 was located near the opening, the sensor detected the hydrogen concentration at 75 s for pitching to bow; however, in the opposite case, it was detected at 120 s, and a verifiable difference of approximately 45 s occurred. For S2 located near the hydrogen leak hole, the hydrogen concentration was detected at 64 s when moving in the same direction (clockwise) as the hydrogen leak direction; however, in the opposite case, it was detected at 8 s. For the roll motion shown in Figure 13 and Figure 14, lateral motion is less affected by the hydrogen leak hole located in the longitudinal direction than pitch motion, and there is no difference in time to reach the sensor; however, the motion cycle remains affected. According to the 15-s cycle, the concentration repeatedly rises and falls, and the overall hydrogen concentration increases. These results suggest that the ship’s motion characteristics must be considered in selecting the sensor location because the detection time for each sensor location varies depending on the location of the hydrogen leakage and the direction of motion.

3.5. Effect of Ship Motion Periods

The numerical analysis was performed by applying pitch motion with periods of 15 s, 60 s, and 120 s to the enclosed area; the results confirmed our understanding of the effects of the vessel’s motion direction on hydrogen diffusion behavior. Table 6 summarizes the time required by hydrogen to reach each sensor according to the ship motion periods compared to the steady state. As the exercise cycle extends, hydrogen diffuses in a manner similar to the normal state and the time difference to reach the sensors is minor, except at sensor S6. The reason for this difference in S6 is due to the start direction of the pitch motion. Figure 15 shows the hydrogen concentration of S4 for each of the ship motion periods. As shown in Figure 15, the hydrogen concentration increases with the motion periods.

4. Conclusions

Hydrogen leakage diffusion under various ship motion conditions was simulated numerically in our study using the CFD code FLUENT for the enclosed area. Furthermore, we analyzed the hydrogen concentration distribution within the enclosed area and drew the following conclusions:
(1)
When the hydrogen leakage is located in the lateral direction, it reaches the sensors approximately 63% (15 s) faster during pitch (lateral) motion than in the steady state. On the other hand, in the case of longitudinal motion, it reaches the sensor approximately 50% (30 s) faster. When the direction of rotation is the same based on the hydrogen leak’s location, the effect on the hydrogen concentration distribution and diffusion is significant;
(2)
We confirmed that the difference in time required by hydrogen to reach the sensor depends on the direction of motion and location of the hydrogen leak. When moving in the same direction as the hydrogen leakage, it reaches the sensor 52% (64 s) faster than in the steady state due to updraft formation; however, in the opposite case, the time required to reach the sensor compared to the steady state increases by 50% (20 s);
(3)
By comparing the results of hydrogen concentration, we confirmed that the difference in time elapsed to reach the sensor, dependent on the motion period of the vessel, is smaller compared to the steady-state; the longer the cycle, the weaker the effect of the cycle period on hydrogen diffusion.
In summary, motion conditions must be considered to understand the leakage and diffusion characteristics of the hydrogen vessel and install a hydrogen sensor because hydrogen’s diffusion behavior varies according to the vessel’s motion characteristics. However, it is also necessary to evaluate the effects of hydrogen concentration diffusion by applying scenarios that combine the roll and pitch motions of a ship and various periods and amplitudes. These considerations can help overcome limitations in the future, such as analysis targets, leakage, and ship motion conditions.

Author Contributions

Conceptualization, B.K. and K.-I.H.; methodology, B.K.; software, B.K.; validation B.K.; formal analysis, B.K. and K.-I.H.; investigation, B.K.; data curation, B.K.; writing—original draft preparation, B.K.; writing—review and editing, B.K.; visualization, B.K.; supervision, K.-I.H.; project administration, K.-I.H.; funding acquisition, K.-I.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the ”Development of guidance for prevent of leak and mitigation of consequence in hydrogen ships”, funded by the Ministry of Oceans and Fisheries, Korea; and the BK21 Four program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education of Korea (Center for Creative Leaders in Maritime Convergence).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Asikainen, A.; Pärjälä, E.; Jantunen, M.; Tuomisto, J.T.; Sabel, C.E. Effects of local greenhouse gas abatement strategies on air pollutant emissions and on health in Kuopio, Finland. Climate 2017, 5, 43. [Google Scholar] [CrossRef] [Green Version]
  2. Walker, T.R.; Adebambo, O.; Del, M. Environmental effects of marine transportation. In Word Seas: An Environmental Evaluation, 2nd ed.; Elsevier: Amsterdam, The Netherlands, 2018; Chapter 27. [Google Scholar]
  3. Kim, H.; Koo, K.Y.; Joung, T.H. A study on the necessity of integrated evaluation of alternative marine fuels. J. Int. Marit. Saf. Environ. Aff. Shipp. 2020, 4, 26–31. [Google Scholar] [CrossRef]
  4. Chen, J.; Fei, Y.; Wan, Z. The relationship between the development of global maritime fleets and GHG emission from shipping. J. Environ. Manag. 2019, 242, 31–39. [Google Scholar] [CrossRef] [PubMed]
  5. IMO. Fourth IMO GHG Study 2020; International Maritime Organization (IMO): London, UK, 2020. [Google Scholar]
  6. IMO. Initial IMO Strategy on Reduction of GHG Emissions from Ships, Resolution MEPC.304 (72) Adopted on 13 April 2018; International Maritime Organization: London, UK, 2018; p. 13. [Google Scholar]
  7. Joung, T.H.; Kang, S.G.; Lee, J.K.; Ahn, J. The IMO initial strategy for reducing Greenhouse Gas (GHG) emissions, and its follow-up actions towards 2050. J. Int. Marit. Saf. Environ. Aff. Shipp. 2020, 4, 1–7. [Google Scholar] [CrossRef] [Green Version]
  8. Endresen, Ø.; Eide, M.; Longva, T. Maritime Forecast. To 2050. Energy Transition Outlook 2019; DNV GL: Hamburg, Germany, 2019; Available online: https://eto.dnvgl.com/2019/download (accessed on 1 December 2021).
  9. LR/UMAS. Zero-Emission Vessels 2030: Transition Pathways; UMAS/Lloyd’s Register: London, UK, 2019. [Google Scholar]
  10. Lide, D. CRC Handbook of Chemistry and Physics, 74th ed.; CRC Press: Boca Raton, FL, USA, 1993. [Google Scholar]
  11. IMO. International Code of Safety for Ships Using Gases or Other Low-Flashpoint Fuels (IGF Code), MSC.391(95) Resolution; IMO (International Maritime Organization): London, UK, 2015. [Google Scholar]
  12. Li, F.; Yuan, Y.; Yan, X.; Malekian, R.; Li, Z. A study on a numerical simulation of the leakage and diffusion of hydrogen in a fuel cell ship. Renew. Sustain. Energy Rev. 2018, 97, 177–185. [Google Scholar] [CrossRef] [Green Version]
  13. Liu, Y.L.; Zheng, Y.Y.; Xu, P.; Zhao, Y.Z.; Bie, H.Y.; Chen, H.G.; Dryver, H. Numerical simulation on the diffusion of hydrogen due to high pressured storage tanks failure. J. Loss Prev. Process Ind. 2009, 22, 265–270. [Google Scholar] [CrossRef]
  14. Dadashzadeh, M.; Ahmad, A.; Khan, F. Dispersion modelling and analysis of hydrogen fuel gas released in an enclosed area: A CFD-based approach. Fuel 2016, 184, 192–201. [Google Scholar] [CrossRef]
  15. Malakhov, A.A.; Avdeenkov, A.V.; du Toit, M.H.; Bessarabov, D.G. CFD simulation and experimental study of a hydrogen leak in a semi-closed space with the purpose of risk mitigation. Int. J. Hydrog. Energy 2020, 45, 9231–9240. [Google Scholar] [CrossRef]
  16. Lacome, J.M.; Jamois, D.; Perrette, D.; Proust, C.H. Large-scale hydrogen release in an isothermal confined area. Int. J. Hydrog. Energy 2011, 36, 2302–2312. [Google Scholar] [CrossRef] [Green Version]
  17. Mao, X.; Ying, R.; Yuan, Y.; Li, F.; Shen, B. Simulation and analysis of hydrogen leakage and explosion behaviors in various compartments on a hydrogen fuel cell ship. Int. J. Hydrog. Energy 2021, 46, 6857–6872. [Google Scholar] [CrossRef]
  18. ANSYS Fluent Theory Guide. Release 19.2; ANSYS: Canonsburg, PA, USA, 2019. [Google Scholar]
  19. Hussein, H.; Brennan, S.; Molkov, V. Dispersion of hydrogen release in a naturally ventilated covered car park. Int. J. Hydrog. Energy 2020, 45, 23882–23897. [Google Scholar] [CrossRef]
  20. Tang, X.; Pu, L.; Shao, X.; Lei, G.; Li, Y.; Wang, X. Dispersion behavior and safety study of liquid hydrogen leakage under different application situations. Int. J. Hydrog. Energy 2020, 45, 31278–31288. [Google Scholar] [CrossRef]
  21. IMO. International Code on Intact Stability, MSC.267(85) Resolution; IMO (International Maritime Organization): London, UK, 2008. [Google Scholar]
  22. Biran, A.; Pulido, R.L. Ship Hydrostatics and Stability; Butterworth-Heinemann: Oxford, UK, 2013. [Google Scholar]
Figure 1. Flowchart of overall methodology.
Figure 1. Flowchart of overall methodology.
Applsci 12 01259 g001
Figure 2. Schematic of experimental setup.
Figure 2. Schematic of experimental setup.
Applsci 12 01259 g002
Figure 3. Position of hydrogen sensors: layout and numbering.
Figure 3. Position of hydrogen sensors: layout and numbering.
Applsci 12 01259 g003
Figure 4. Grid independence test (H = 0.735 m, Time = 60 s).
Figure 4. Grid independence test (H = 0.735 m, Time = 60 s).
Applsci 12 01259 g004
Figure 5. Comparisons between predicted hydrogen concentrations by experimental and FLUENT data.
Figure 5. Comparisons between predicted hydrogen concentrations by experimental and FLUENT data.
Applsci 12 01259 g005
Figure 6. Six-degree of freedom ship motions.
Figure 6. Six-degree of freedom ship motions.
Applsci 12 01259 g006
Figure 7. Plane view and ship motion period graph: (a) Pitching to bow up; (b) Rolling to port up; (c) Clockwise.
Figure 7. Plane view and ship motion period graph: (a) Pitching to bow up; (b) Rolling to port up; (c) Clockwise.
Applsci 12 01259 g007
Figure 8. Plane view and ship motion period graph: (a) Pitching to bow down; (b) Rolling to port down; (c) Counterclockwise.
Figure 8. Plane view and ship motion period graph: (a) Pitching to bow down; (b) Rolling to port down; (c) Counterclockwise.
Applsci 12 01259 g008
Figure 9. Hydrogen concentration by roll or pitch motion at S4.
Figure 9. Hydrogen concentration by roll or pitch motion at S4.
Applsci 12 01259 g009
Figure 10. Hydrogen concentration by roll or pitch motion at S6.
Figure 10. Hydrogen concentration by roll or pitch motion at S6.
Applsci 12 01259 g010
Figure 11. Central longitudinal cross-section of hydrogen distribution: (a) Steady state; (b) Pitching to bow up; (c) Pitching to bow down.
Figure 11. Central longitudinal cross-section of hydrogen distribution: (a) Steady state; (b) Pitching to bow up; (c) Pitching to bow down.
Applsci 12 01259 g011
Figure 12. Central transversal cross-section of the hydrogen distribution: (a) Steady state; (b) Rolling to port up; (c) Rolling to port down.
Figure 12. Central transversal cross-section of the hydrogen distribution: (a) Steady state; (b) Rolling to port up; (c) Rolling to port down.
Applsci 12 01259 g012
Figure 13. Hydrogen concentration by roll and pitch motion at S2.
Figure 13. Hydrogen concentration by roll and pitch motion at S2.
Applsci 12 01259 g013
Figure 14. Hydrogen concentration by roll and pitch motion at S6.
Figure 14. Hydrogen concentration by roll and pitch motion at S6.
Applsci 12 01259 g014
Figure 15. Hydrogen concentration by pitch motion at S4: (a) 15 s period; (b) 60 s period; (c) 120 s period.
Figure 15. Hydrogen concentration by pitch motion at S4: (a) 15 s period; (b) 60 s period; (c) 120 s period.
Applsci 12 01259 g015
Table 1. Recent studies related to hydrogen release and diffusion.
Table 1. Recent studies related to hydrogen release and diffusion.
ObjectMethodParametersRef
Initial
Pressure
Release
Flow Rate
Diameter
[mm]
Hydrogen fuel cell shipSimulation/FLUENT10 MPa0.46710[12]
Hydrogen storage tankSimulation/FLUENT42 MPa-0.5, 1, 2, 5[13]
Hydrogen fuel cell vehicleSimulation/FDS-6.7 × 10−4 kg/s-[14]
Semi-closed spaceSimulation/STAR-CCM+2, 5, 20 bar-0.4572, 0.8[15]
Simple geometric spacesExperiment-0.2, 1 kg/s5, 10, 20[16]
Table 2. Experiment conditions.
Table 2. Experiment conditions.
Experimental facilityW × L × H (m)0.5 × 1.0 × 0.75
Volume (m3)0.1596
Hydrogen leakDiameter (m)0.001587 (1/16 inch)
Inlet (kg/s)0.965 × 10−5
Pressure (Pa)1,400,000
Natural ventilation openingDiameter (m)0.01
Table 3. Ship motion.
Table 3. Ship motion.
Translation or
Rotation
Type of
Ship Motion
AxisPositive Sense
TranslationSurgeAlong xForwards
SwayAlong yTo starboard
HeaveAlong zUpwards
RotationRollAbout xPort up
PitchAbout yBow up
YawAbout zStarboard to bow
Table 4. Time for hydrogen to reach each sensor according to the ship motion types.
Table 4. Time for hydrogen to reach each sensor according to the ship motion types.
S2S3S4S5S6
Steady state60 s53 s40 s53 s124 s
60 s period Pitching to bow up45 s
(−23%)
30 s
(−43%)
15 s
(−63%)
30 s
(−43%)
60 s
(−52%)
60 s period Rolling to port up30 s
(−50%)
45 s
(−15%)
30 s
(−25%)
45 s
(−15%)
90 s
(−27%)
Table 5. Time for hydrogen to reach each sensor according to the ship motion directions.
Table 5. Time for hydrogen to reach each sensor according to the ship motion directions.
S2S3S4S5S6
Steady state60 s53 s40 s53 s124 s
Pitch motion
15 s period
Bow up68 s
(+13%)
49 s
(−8%)
23 s
(−43%)
49 s
(−8%)
60 s
(−52%)
Bow down53 s
(−12%)
79 s
(+49%)
60 s
(+50%)
79 s
(+49%)
101 s
(−19%)
Roll motion
15 s period
Port up45 s
(−25%)
41 s
(−23%)
65 s
(+63%)
64 s
(+21%)
87 s
(−30%)
Port down45 s
(−25%)
64 s
(+21%)
65 s
(+63%)
41 s
(−23%)
87 s
(−30%)
Table 6. Time for hydrogen to reach each sensor according to ship motion periods.
Table 6. Time for hydrogen to reach each sensor according to ship motion periods.
S2S3S4S5S6
Steady state60 s53 s40 s53 s124 s
15 s period Pitching to bow up68 s
(+13%)
49 s
(−8%)
23 s
(−43%)
49 s
(−8%)
60 s
(−52%)
60 s period Pitching to bow up45 s
(−23%)
30 s
(−43%)
15 s
(−63%)
30 s
(−43%)
60 s
(−52%)
120 s period Pitching to bow up60 s
(0%)
60 s
(+13%)
40 s
(0%)
60 s
(+13%)
60 s
(−52%)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kim, B.; Hwang, K.-I. Numerical Analysis of the Effects of Ship Motion on Hydrogen Release and Dispersion in an Enclosed Area. Appl. Sci. 2022, 12, 1259. https://doi.org/10.3390/app12031259

AMA Style

Kim B, Hwang K-I. Numerical Analysis of the Effects of Ship Motion on Hydrogen Release and Dispersion in an Enclosed Area. Applied Sciences. 2022; 12(3):1259. https://doi.org/10.3390/app12031259

Chicago/Turabian Style

Kim, Byeol, and Kwang-Il Hwang. 2022. "Numerical Analysis of the Effects of Ship Motion on Hydrogen Release and Dispersion in an Enclosed Area" Applied Sciences 12, no. 3: 1259. https://doi.org/10.3390/app12031259

APA Style

Kim, B., & Hwang, K. -I. (2022). Numerical Analysis of the Effects of Ship Motion on Hydrogen Release and Dispersion in an Enclosed Area. Applied Sciences, 12(3), 1259. https://doi.org/10.3390/app12031259

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