Numerical Study of the Action of Convection on the Volume and Length of the Flammable Zone Formed by Hydrogen Emissions from the Vent Masts Installed on an International Ship

: International ships carrying liqueﬁed fuel are strongly recommended to install vent masts to control the pressure of cargo tanks in the event of an emergency. However, the gas emitted from a vent mast may be hazardous for the crew of the ship. In the present study, the volume and length of the ﬂammable zone (FZ) created by the emitted gas above the ship was examined. Various scenarios comprising four parameters, namely, relative wind speed, arrangement of vent masts, combination of emissions among four vent masts, and direction of emission from the vent-mast outlet were considered. The results showed that the convection acts on the volume and length of an FZ. The volume of an FZ increases when there is a reduction in convection reaching the FZ and when strong convection brings hydrogen from a nearby FZ. The length of the FZ is also related to convection. An FZ is elongated if the center of a vortex is located inside the FZ, because this vortex traps hydrogen inside the FZ. The length of an FZ decreases if the center of the vortex is located outside the FZ, as such a vortex brings more fresh air into the FZ.


Introduction
The threat of climate change is a significant global issue, according to the Paris Agreement. Most of the countries bound by this agreement must make significant efforts to reduce greenhouse gas (GHG) emissions as soon as possible [1], since large GHG emissions accelerate global warming. In this context, the International Maritime Organization (IMO) has agreed to limit GHG emissions, with the goal of completely eliminating GHGs produced by international shipping by the end of the century [2].
GHG emissions from the use of fossil fuels are unavoidable since GHG contains not only carbon dioxide but also methane, the principal component of natural gas, which is increasingly being used as a substitute for traditional fossil fuels [3,4]. The IMO's goal includes the use of carbon-free fuels such as hydrogen gas in international ships. As a result, they are working to update previous safety requirements, such as the IGF code [5], so that they can be applied to ships that use hydrogen as a fuel [6].
New Energy and Industrial Technology Development Organization (NEDO) also introduced the supply chain of hydrogen [7]. According to the CO 2 -free Hydrogen Energy Supply-chain Technology Research Association (HyStra), the hydrogen supply chain requires the development of a liquefied hydrogen gas transporter [8]. A study on reliquefying systems evaluated the boil-off gas from liquefied hydrogen cargo tanks in these circumstances [9]. It assumed a natural boil-off rate of approximately 0.3% per day.

Modeling
The ship model has four tanks of 40,000 m 3 each, representing a membrane-type liquefied gas carrier of about 50,000 dead weight tonnage (DWT) [19]. The part of the ship below the design load water line (DLWL) is not considered. The four vent masts are set on the top of each tank. It maintains an adequate distance for safety of the inhabitants. In compliance with the safety design rules in the IGC code, the height of each vent mast is 16 m, which is one-third of the width of the ship. The detailed geometry of the ship is shown in Figure 1. The atmospheric space is adequate, as shown in Figure 2. The properties of gas and air are summarized in Table 1. The pressure (P h ) of the computational domain is estimated from Equation (1), based on the height (h) from the DLWL.
where atmospheric pressure (P atm ) is 101,325 Pa; ρ humid is the density of the atmospheric air as calculated from Equation (2); and g is acceleration due to gravity (9.81 m/s 2 ).
ρ atm = x vapor ρ vapor + x dryair ρ dryair (2) where x vapor and x dryair are the volume fractions of vapor and dry air, respectively.
shown in Figure 1. The atmospheric space is adequate, as shown in Figure 2. The properties of gas and air are summarized in Table 1. The pressure (Ph) of the computational domain is estimated from Equation (1), based on the height (h) from the DLWL.
where atmospheric pressure (Patm) is 101,325 Pa; is the density of the atmospheric air as calculated from Equation (2); and is acceleration due to gravity (9.81 m/s ).
= + (2) where and are the volume fractions of vapor and dry air, respectively.

CFD Solvers
In the present simulation, temperature solvers were not used. Instead, a segregated flow solver was used. The gas density was assumed constant. The Reynolds-averaged Navier-Stokes and k-turbulence models were used to solve the fluid mechanics. Fuel emissions were supposed to be ceaseless enough to form a stable FZ at a steady relative wind speed. Atmospheric air was assumed to consist of dry air and vapor. The gravity force had a buoyancy effect on the hydrogen in the air.     Table 2 lists the variables considered for estimating the minimum required discharge flow rate of hydrogen as defined in the IGC code [10]. The temperature at the relieving conditions was assumed as 53.15 K. The size of the LH2 tanks was estimated to be fit for the ship capacity, as shown in Figure 1. The discharging rate of gas from the ith tank was estimated from Equation (3). The gas factor G and the constant coefficient D were defined as per Equations (4) and (5), respectively. The mass flow rate of gas through each vent mast of the tank was estimated from Equation (6), where m i is the minimum required mass flow rate of gas from the ith tank.

Parameters
The following four parameters were used in the numerical study on the volume and length of the FZ. The relative wind speeds considered were 10, 20, and 30 m/s. These speeds assume that the wind blows on directly against the voyage direction of the ship at 0, 10, and 20 m/s, while the ship goes forward as a typical speed of 10 m/s [19]. A fresh gale corresponds to a wind speed of 20 m/s [27].
Four types of arrangement of vent masts were considered ( Figure 3). The reference arrangement is shown in Figure 3a. As shown in Figure 3c     Emissions from four vent masts can yield several possible combinations. Each scenario consists of a combination of emissions from four vent masts. The volumes of the FZs for these combinations were examined to identify how multiple emissions affect the volume and length of an FZ. The relative volume ratio of an FZ is defined as per Equation (7).
where V fz is the total volume of the FZ, and V fz,1 is the volume of the FZ when hydrogen is emitted from the first vent mast alone. The criteria for the FZ are on the volume fraction of hydrogen. This classification is based on the ClassNK group guideline [20] that introduces the criterion of standard flammability limit (FL) on the basis of hydrogen volume fraction. Finally, three types of vent-mast outlets were designed to determine the effect of direction of emissions. The original design of vent mast outlet emits hydrogen into two directions, namely, upward and downward, as shown in Figure 4a, following the work of Kang [15]. The other two kinds of vent-mast outlet emit hydrogen either only downward or upward, as shown in Figure 4b,c, respectively. In all the cases, the vent-mast pipe diameter is set as 0.6 m.    It has been chosen for grid number, from a grid dependency test, on a volume of the FZ of a representative scenario. This scenario is the reference arrangement with a relative wind speed of 10 m/s. All the vent masts emit hydrogen from the outlet in upward and downward directions. It is less than 1.19% in a differential rate of volume of an FZ and is the ratio of the difference between the volumes obtained from two grid numbers of 26 and 18 million. Therefore, the number of grids for the computational domain in all other scenarios was chosen as 18 million.
The mass flow rate of hydrogen for each vent mast was determined from Equation (6). It is 2.37 kg/s for the first vent mast and 2.65 kg/s for the other three vent masts. These correspond to volume flow rates of 28.16 m 3 /s and 31.54 m 3 /s, respectively, under the conditions listed in Table 2.
The FZ is sliced at 5 m intervals to inspect the FZ more closely at the cross sections. These sections begin from the center of the first vent-mast outlet and follow the drifting direction of hydrogen. The contour lines and vectors are drawn on these sections. These contour lines are distinguished by the ranges of volume fractions of hydrogen, as listed in Table 3. These ranges contain not only a standard volume fraction of the lower flammability limit (LFL), but also nearby fractions because the flame can extend up to the region of the lower volume fraction of 0.5%, and the ignition of a hydrogen jet is restrained at the higher volume fraction range of 8-10% [28]. Hydrogen does not explode if the fraction is above 75% of the volume fraction-this limit is known as the upper flammability limit (UFL). In the present study, the FZ is defined in the range of volume fractions between the LFL and UFL.

Name of Contour Line
Range of Volume Fraction (%) Color 12.5% of LFL 0.5-1.0 Gray In the zigzag arrangement, shifts the first and third vent masts are shifted each by toward the left side of the ship (Figure 3d). Emissions from four vent masts can yield several possible combinations. Each nario consists of a combination of emissions from four vent masts. The volumes of th for these combinations were examined to identify how multiple emissions affect the ume and length of an FZ. The relative volume ratio of an FZ is defined as per equatio , where Vfz is the total volume of the FZ, and Vfz,1 is the volume of the FZ when hydr is emitted from the first vent mast alone. The criteria for the FZ are on the volume fra of hydrogen. This classification is based on the ClassNK group guideline [21] that i duces the criterion of standard flammability limit (FL) on the basis of hydrogen vo fraction.
Finally, three types of vent-mast outlets were designed to determine the effect rection of emissions. The original design of vent mast outlet emits hydrogen into tw rections, namely, upward and downward, as shown in Figure 4a, following the wo Kang [16]. The other two kinds of vent-mast outlet emit hydrogen either only down or upward as shown in Figure 4b and 4c, respectively. In all the cases, the vent-mast diameter is set as 0.6 m.
It has been chosen for grid number from a grid dependency test on a volume o FZ of a representative scenario. This scenario is the reference arrangement with a rel wind speed of 10 m/s. All the vent masts emit hydrogen from the outlet in upward downward directions. It is less than 1.19% in a differential rate of volume of an FZ a the ratio of the difference between the volumes obtained from two grid numbers of 2 18 million. Therefore, the number of grids for the computational domain in all othe narios was chosen as 18 million.
The mass flow rate of hydrogen for each vent mast was determined from formul It is 2.37 kg/s for the first vent mast and 2.65 kg/s for the other three vent masts. T correspond to volume flow rates of 28.16 m 3 /s and 31.54 m 3 /s, respectively, under the ditions listed in Table 2.
The FZ is sliced at 5-m intervals to inspect the FZ more closely at the cross sect These sections begin from the center of the first vent-mast outlet and follow the dr direction of hydrogen. The contour lines and vectors are drawn on these sections. T contour lines are distinguished by the ranges of volume fractions of hydrogen, as list Table 3. These ranges contain not only a standard volume fraction of the lower flam bility limit (LFL), but also nearby fractions because the flame can extend up to the re of the lower volume fraction of 0.5%, and the ignition of a hydrogen jet is restrained a higher volume fraction range of 8%-10% [29]. Hydrogen does not explode if the fra is above 75% of the volume fraction-this limit is called as the upper flammability (UFL). In present study, the FZ is defined in the range of volume fractions betwee LFL and UFL. In the zigzag arrangement, shifts the first and third vent masts are shifted each b toward the left side of the ship (Figure 3d). Emissions from four vent masts can yield several possible combinations. Eac nario consists of a combination of emissions from four vent masts. The volumes of th for these combinations were examined to identify how multiple emissions affect th ume and length of an FZ. The relative volume ratio of an FZ is defined as per equatio , where Vfz is the total volume of the FZ, and Vfz,1 is the volume of the FZ when hyd is emitted from the first vent mast alone. The criteria for the FZ are on the volume fra of hydrogen. This classification is based on the ClassNK group guideline [21] that duces the criterion of standard flammability limit (FL) on the basis of hydrogen vo fraction.
Finally, three types of vent-mast outlets were designed to determine the effect rection of emissions. The original design of vent mast outlet emits hydrogen into tw rections, namely, upward and downward, as shown in Figure 4a, following the wo Kang [16]. The other two kinds of vent-mast outlet emit hydrogen either only down or upward as shown in Figure 4b and 4c, respectively. In all the cases, the vent-mas diameter is set as 0.6 m.
It has been chosen for grid number from a grid dependency test on a volume FZ of a representative scenario. This scenario is the reference arrangement with a re wind speed of 10 m/s. All the vent masts emit hydrogen from the outlet in upwar downward directions. It is less than 1.19% in a differential rate of volume of an FZ a the ratio of the difference between the volumes obtained from two grid numbers of 2 18 million. Therefore, the number of grids for the computational domain in all othe narios was chosen as 18 million.
The mass flow rate of hydrogen for each vent mast was determined from formu It is 2.37 kg/s for the first vent mast and 2.65 kg/s for the other three vent masts. correspond to volume flow rates of 28.16 m 3 /s and 31.54 m 3 /s, respectively, under the ditions listed in Table 2.
The FZ is sliced at 5-m intervals to inspect the FZ more closely at the cross sec These sections begin from the center of the first vent-mast outlet and follow the dr direction of hydrogen. The contour lines and vectors are drawn on these sections. contour lines are distinguished by the ranges of volume fractions of hydrogen, as lis Table 3. These ranges contain not only a standard volume fraction of the lower fla bility limit (LFL), but also nearby fractions because the flame can extend up to the r of the lower volume fraction of 0.5%, and the ignition of a hydrogen jet is restrained higher volume fraction range of 8%-10% [29]. Hydrogen does not explode if the fra is above 75% of the volume fraction-this limit is called as the upper flammability (UFL). In present study, the FZ is defined in the range of volume fractions betwee LFL and UFL. In the zigzag arrangement, shifts the first and third vent masts are shifted each by toward the left side of the ship (Figure 3d). Emissions from four vent masts can yield several possible combinations. Each nario consists of a combination of emissions from four vent masts. The volumes of th for these combinations were examined to identify how multiple emissions affect the ume and length of an FZ. The relative volume ratio of an FZ is defined as per equatio , where Vfz is the total volume of the FZ, and Vfz,1 is the volume of the FZ when hydr is emitted from the first vent mast alone. The criteria for the FZ are on the volume fra of hydrogen. This classification is based on the ClassNK group guideline [21] that i duces the criterion of standard flammability limit (FL) on the basis of hydrogen vo fraction.
Finally, three types of vent-mast outlets were designed to determine the effect o rection of emissions. The original design of vent mast outlet emits hydrogen into tw rections, namely, upward and downward, as shown in Figure 4a, following the wo Kang [16]. The other two kinds of vent-mast outlet emit hydrogen either only down or upward as shown in Figure 4b and 4c, respectively. In all the cases, the vent-mast diameter is set as 0.6 m.
It has been chosen for grid number from a grid dependency test on a volume o FZ of a representative scenario. This scenario is the reference arrangement with a rel wind speed of 10 m/s. All the vent masts emit hydrogen from the outlet in upward downward directions. It is less than 1.19% in a differential rate of volume of an FZ a the ratio of the difference between the volumes obtained from two grid numbers of 26 18 million. Therefore, the number of grids for the computational domain in all other narios was chosen as 18 million.
The mass flow rate of hydrogen for each vent mast was determined from formul It is 2.37 kg/s for the first vent mast and 2.65 kg/s for the other three vent masts. T correspond to volume flow rates of 28.16 m 3 /s and 31.54 m 3 /s, respectively, under the ditions listed in Table 2.
The FZ is sliced at 5-m intervals to inspect the FZ more closely at the cross sect These sections begin from the center of the first vent-mast outlet and follow the dri direction of hydrogen. The contour lines and vectors are drawn on these sections. T contour lines are distinguished by the ranges of volume fractions of hydrogen, as list Table 3. These ranges contain not only a standard volume fraction of the lower flam bility limit (LFL), but also nearby fractions because the flame can extend up to the re of the lower volume fraction of 0.5%, and the ignition of a hydrogen jet is restrained a higher volume fraction range of 8%-10% [29]. Hydrogen does not explode if the fra is above 75% of the volume fraction-this limit is called as the upper flammability (UFL). In present study, the FZ is defined in the range of volume fractions betwee LFL and UFL.  Figure 5a shows the contour lines on a cross section to denote the FL. Two types of lines are used as the criterion for a cross section. The first type is a vent-mast line shown in Figure 3a; it is used for the cross section shown in Figures 5, 6 and 10a,b. The second type is a vertical line (as marked in Figure 10a) perpendicular to the drifting direction. The vertical line is a criterion for the cross sections shown in Figures 10c,d and 13. Figure 5a shows the contour lines on a cross section to denote the FL. Two types of lines are used as the criterion for a cross section. The first type is a vent-mast line shown in Figure 3a; it is used for the cross section shown in Figures 5, 6 and 10a,b. The second type is a vertical line (as marked in Figure 10a) perpendicular to the drifting direction. The vertical line is a criterion for the cross sections shown in Figures 10c,d and 13.

Relative Wind Speed
In Figure 5, the FZ at each relative wind speed is shown. The minimum distance of the FZ from the steer room was measured. The distances were 8.7, 9.2, and 11.5 m at relative wind speeds of 30, 20, and 10 m/s, respectively. Thus, the distance from the steer room decreases as the relative wind speed increases, as shown in Figure 5e.

Relative Wind Speed
In Figure 5, the FZ at each relative wind speed is shown. The minimum distance of the FZ from the steer room was measured. The distances were 8.7, 9.2, and 11.5 m at relative wind speeds of 30, 20, and 10 m/s, respectively. Thus, the distance from the steer room decreases as the relative wind speed increases, as shown in Figure 5e. Figure 6 shows the contour lines of the FL for each arrangement. The LFL merged when any of the two vent masts were close. Figure 7 compares the total volume of the FZ with four types of arrangements using the plot of the average volume. The fluctuation of the volume of the FZ (Vfz) during iteration is plotted as the error bar. At a relative wind speed of 10 m/s, the reference arrangement records the lowest volumes of 7791 m 3 and 4594 m 3 for the 1-23-4 arrangement at a relative wind speed of 20 m/s. Figure 8 shows the area of the FZ at 5 m spacings in the drifting direction. The length of the FZ is defined from this plot as the length for which the area of the FZ becomes zero. This length is 220 m for the reference arrangement and 240 m for the 1-23-4 arrangement at a relative wind speed of 10 m/s; at 20 m/s, these lengths are 245 and 230 m, respectively. In other words, the length increased for the reference arrangement, while the length decreased for the 1-  Figure 6 shows the contour lines of the FL for each arrangement. The LFL merged when any of the two vent masts were close. Figure 7 compares the total volume of the FZ with four types of arrangements using the plot of the average volume. The fluctuation of the volume of the FZ (V fz ) during iteration is plotted as the error bar. At a relative wind speed of 10 m/s, the reference arrangement records the lowest volumes of 7791 m 3 and 4594 m 3 for the 1-23-4 arrangement at a relative wind speed of 20 m/s. Figure 8 shows the area of the FZ at 5 m spacings in the drifting direction. The length of the FZ is defined from this plot as the length for which the area of the FZ becomes zero. This length is 220 m for the reference arrangement and 240 m for the 1-23-4 arrangement at a relative wind speed of 10 m/s; at 20 m/s, these lengths are 245 and 230 m, respectively. In other words, the length increased for the reference arrangement, while the length decreased for the 1-23-4 arrangement as the relative wind speed increased from 10 m/s to 20 m/s. Hence, the rank of volume changed dramatically for the 1-23-4 arrangement; it shows the highest rank at 10 m/s, but the lowest rank at 20 m/s as in Figure 7.

Arrangement of Vent Masts
23-4 arrangement as the relative wind speed increased from 10 m/s to 20 m/s. Hence, the rank of volume changed dramatically for the 1-23-4 arrangement; it shows the highest rank at 10 m/s, but the lowest rank at 20 m/s as in Figure 7.

Combination of Emission Patterns Selected among the Four Vent Masts
Each of the four vent masts can emit hydrogen. Different emission scenarios are generated by selecting vents for emission from among the four available vent masts. Figure 9 compares the relative volume ratio of the scenarios. This is a ratio of the volumes of the FZs obtained using Equation (7). The denominator is the volume in an emission scenario 23-4 arrangement as the relative wind speed increased from 10 m/s to 20 m/s. Hence, the rank of volume changed dramatically for the 1-23-4 arrangement; it shows the highest rank at 10 m/s, but the lowest rank at 20 m/s as in Figure 7.

Combination of Emission Patterns Selected among the Four Vent Masts
Each of the four vent masts can emit hydrogen. Different emission scenarios are generated by selecting vents for emission from among the four available vent masts. Figure 9 compares the relative volume ratio of the scenarios. This is a ratio of the volumes of the FZs obtained using Equation (7). The denominator is the volume in an emission scenario

Combination of Emission Patterns Selected among the Four Vent Masts
Each of the four vent masts can emit hydrogen. Different emission scenarios are generated by selecting vents for emission from among the four available vent masts. Figure 9 compares the relative volume ratio of the scenarios. This is a ratio of the volumes of the FZs obtained using Equation (7). The denominator is the volume in an emission scenario where hydrogen is emitted from the first vent mast alone. The numerator is the volume for a specific scenario of emissions in which specific vents from among the four are selected. In the scenario that the second vent mast alone emits hydrogen, the differences were not greater than 3% compared with the scenario in which hydrogen is emitted from the first vent mast alone. When a scenario of emission from two vents is considered, the volume increases greatly to 232%, which exceeds 32% by 200%. The volume ratio is significantly increased to over 400% for the scenarios of emissions from all the four vent masts, and it increases to at least 497% or up to 796% as the maximum values for relative wind speeds of 10 and 30 m/s, respectively. Two scenarios were used to identify the reason for these increases. In one scenario, emission occurs from the second vent mast alone, as shown in Figure 10a. In the other one, simultaneous emissions occur from the first and second vent masts, as shown in Figure 10b. The contour lines of the FL are plotted on the cross sections, as shown in Figure 10, to distinguish the lower region from the upper region of the FZ. The convection is shown by the fluid vectors. It is seen that convection acts on the FZ. The convection in the lower region of the FZ is weaker in the scenario of simultaneous emissions from the two vent masts compared with the scenario of emission from the second vent mast alone, as shown in Figure 10c,d. The convection-which dilutes the FZ by providing fresh air-acting on the lower region of the FZ is reduced for the scenario of simultaneous emissions. This reduction increases the volume of the FZ. Meanwhile, the upper region of the FZ rapidly floats with a strong uplifting of convection, as shown in Figure 10d where emissions are simultaneous for the two vent masts. This rapid floating allows the hydrogen to overlap with others, namely, the hydrogen emitted from the second vent mast overlaps with the hydrogen emitted from the first vent mast. This overlap broadens the FZ, as shown in Figure 10d. As a result, the volume of the FZ is 232%.

Direction of Emission from the Vent-Mast Outlet
Three directions of emission are shown in Figure 4. Figure 11 compares the total volume of the FZ for these three directions. A scenario of emission wherein all four vent masts emit hydrogen is considered. The differences in the volume of the FZ for each direction of the vent mast outlet are shown in the figure. The volume is 7809 m 3 for the upward and downward outlets and 7738 m 3 for the downward outlet at a relative wind speed of 10 m/s. In other words, the volume of the FZ is lower in the downward direction than in the upward and downward directions of the outlet at a relative wind speed of 10 m/s. for a specific scenario of emissions in which specific vents from among t lected. In the scenario that the second vent mast alone emits hydrogen, were not greater than 3% compared with the scenario in which hydrogen the first vent mast alone. When a scenario of emission from two vents is volume increases greatly to 232%, which exceeds 32% by 200%. The volu nificantly increased to over 400% for the scenarios of emissions from al masts, and it increases to at least 497% or up to 796% as the maximum val wind speeds of 10 and 30 m/s, respectively. Two scenarios were used to ide for these increases. In one scenario, emission occurs from the second vent shown in Figure 10a. In the other one, simultaneous emissions occur fro second vent masts, as shown in Figure 10b. The contour lines of the FL are cross sections, as shown in Figure 10, to distinguish the lower region from gion of the FZ. The convection is shown by the fluid vectors. It is seen that on the FZ. The convection in the lower region of the FZ is weaker in the s ultaneous emissions from the two vent masts compared with the scena from the second vent mast alone, as shown in Figure 10c,d. The convectionthe FZ by providing fresh air-acting on the lower region of the FZ is r scenario of simultaneous emissions. This reduction increases the volume o while, the upper region of the FZ rapidly floats with a strong uplifting of shown in Figure 10d where emissions are simultaneous for the two vent m floating allows the hydrogen to overlap with others, namely, the hydroge the second vent mast overlaps with the hydrogen emitted from the first v overlap broadens the FZ, as shown in Figure 10d. As a result, the volum 232%.

Direction of Emission from the Vent-Mast Outlet
Three directions of emission are shown in Figure 4. Figure 11 compares the total volume of the FZ for these three directions. A scenario of emission wherein all four vent masts emit hydrogen is considered. The differences in the volume of the FZ for each direction of the vent mast outlet are shown in the figure. The volume is 7809 m 3 for the upward and downward outlets and 7738 m 3 for the downward outlet at a relative wind speed of 10 m/s. In other words, the volume of the FZ is lower in the downward direction than in the upward and downward directions of the outlet at a relative wind speed of 10 m/s. Additional scenarios of emission were created by using the first vent mast alone to identify the different actions of convection on the FZ for two types of outlet directions. The downward emission is compared with the upward and downward emissions to compare the length of the FZ. The length is defined by the drifting distance where the area of the FZ becomes zero, as mentioned in Section 3.2. The drifting distance starts from the center of the first vent mast outlet. It is 10.3° in the clockwise direction from a parallel line of the DLWL, as shown in Figure 10. The area of the FZ is plotted for every 5 m of the cross sections for each scenario, as shown in Figure 12. The figure shows the differences in the lengths of the FZ. The length is longer for outlets with upward and downward emissions than for the scenario with the downward emission only. The length for upward and downward emissions is 160 m, and that for only downward emission is 150 m.
Further, a quantitative method was used to compare the flammable areas of the two outlet types. The summation of the flammable area is defined as an integrated volume.   Additional scenarios of emission were created by using the first vent mast alone to identify the different actions of convection on the FZ for two types of outlet directions. The downward emission is compared with the upward and downward emissions to compare the length of the FZ. The length is defined by the drifting distance where the area of the FZ becomes zero, as mentioned in Section 3.2. The drifting distance starts from the center of the first vent mast outlet. It is 10.3 • in the clockwise direction from a parallel line of the DLWL, as shown in Figure 10. The area of the FZ is plotted for every 5 m of the cross sections for each scenario, as shown in Figure 12. The figure shows the differences in the lengths of the FZ. The length is longer for outlets with upward and downward emissions than for the scenario with the downward emission only. The length for upward and downward emissions is 160 m, and that for only downward emission is 150 m.  Further, a quantitative method was used to compare the flammable areas of the two outlet types. The summation of the flammable area is defined as an integrated volume. The integrated volume is the approximate integration of the area of the FZ (A fz ) based on the drifting distance (d d ). It is estimated using the trapezoidal rule [29]. The angle of the drifting direction from the DLWL differs between the scenarios of the different outlet directions. The angle is 10.3 • for the outlet for upward and downward emission and 4.9 • for the outlet for the outlet for downward emission. The integrated volume for the outlet for upward and downward emission is 2206 m 3 up to a drifting distance of 100 m and 2339 m 3 for the outlet for downward emission. The volume emitted in the downward direction is more than that emitted in both directions only by 6%. Up to 110 m, the outlet for downwards emission showed an increase of 3% compared to the outlet for emission in both directions. This is a decrease of 3% from the integration up to a distance of 100 m. This decrease is originated from a different location of the center of the vortex inside the FZ. Cross sections of the vertical line at drifting distances of 100 and 110 m for each scenario are shown in Figure 13. This figure shows not only the contour lines of the FL, but also the tangential velocity vectors. The center of the vortex is inside the FZ in the scenario of emission in both directions from the outlet, while it is located outside the FZ in the scenario of downward emission. It is expected that the vortex dissipates the FZ more efficiently when the center of the vortex is located outside, as it brings more fresh air into the FZ.

389
In the present study, four parameters that affect on the volume and length of an FZ 390 were examined. These parameters are the relative wind speed, arrangement of vent masts, 391 combination of the emitting vent masts, and direction of emission from the outlet of vent 392 mast. Various scenarios were built with the four parameters. The volume and length of 393 the FZ for each scenario were examined. The results showed that convection contributes 394

Conclusions
In the present study, four parameters that affect the volume and length of an FZ were examined. These parameters are the relative wind speed, arrangement of vent masts, combination of the emitting vent masts, and direction of emission from the outlet of vent mast. Various scenarios were built with the four parameters. The volume and length of the FZ for each scenario were examined. The results showed that convection contributes to the differences between the scenarios. Convection dissipates the FZ by bringing in fresh air. Therefore, a reduction in convection prevents the dissipation of the FZ, and this may result in a significant increase of volume of the FZ. However, this increase cannot be explained only by the reduction in convection. A strong convection can also broaden the FZ by supplying it with hydrogen from a nearby FZ, and the location of the vortex affects the length of the FZ. The FZ is elongated by the trapping of hydrogen via recirculation inside the FZ. On the other hand, the FZ may be shortened by the vortex providing fresh air into the FZ if the center of the vortex is located outside the FZ.

Conflicts of Interest:
The authors declare no conflict of interest.

3D
Three dimensional T rel Temperature at relieving conditions; The boiled gas is assumed to be heated more 32.9 K at relieving conditions (K) UFL Upper flammability limit