#### Abstract

Semianalytical solutions are developed for turbulent hydrogen-air plume. We derived analytical expressions for plume centerline variables (radius, velocity, and density deficit) in terms of a single universal function, called plume function. By combining the obtained analytical expressions of centerline variables with empirical Gaussian expressions of the mean variables, we obtain semianalytical expressions for mean quantities of hydrogen-air plume (velocity, density deficit, and mass fraction).

#### 1. Introduction

One of the important safety issues of hydrogen energy is the hydrogen leakage into ambient air and the associated risk of fire or explosion. In fact, industry has already produced several prototype products using hydrogen as a fuel. Unfortunately, these products are not yet available for commercial use because of safety concerns related to hydrogen leakage. So studying hydrogen-air behavior is very important in order to estimate expected hazards from leakage as well as to propose recommendations when designing hydrogen-related facilities.

Recently, El-Amin and coauthors [1–6] studied the problem of hydrogen leakage in air. In [1–3], they introduced boundary layer theory approach to model the concentration layer adjacent to a ceiling wall at the impinging and far regions in both planar and axisymmetric cases for small-scale hydrogen leakage. While in [4–6], they studied the turbulent hydrogen-air jet/plume resulted from hydrogen leakage in open air. The laminar hydrogen jet is analyzed by Sánchez–Sanz et al. [7]. Also, experimental measurements for turbulent hydrogen jet have been performed by Schefer and coauthors (e.g., [8–10]). On the other hand, CFD simulations of the problem have been done by many researcher such as Matssura and coauthors [11–14], Kikukawa [15], Agarant et al.[16], and Swain et al. [17, 18].

Hydrogen-air jet is an example of non-Boussinesq plume; since the initial fractional density difference is high. The initial fractional density difference is defined as , where is the initial centerline density (density at the source) and is the ambient density. As an example, the initial fractional density differences for selected binary low-density gases at temperature 15°C are 0.93 for H_{2}-Air, 0.86 for He-Air, 0.43 for CH_{4}-Air, and 0.06 for C_{2}H_{2}-N_{2}. Crapper and Baines [19] suggested that the upper bound of applicability of the Boussinesq approximation is that the initial fractional density difference does not exceed 0.05. In general, one can say that the Boussinesq approximation is valid for small initial fractional density difference, (e.g., El**-**Amin and Kanayama [5]). This is correct only for the case of a plume produced by a positive source of buoyancy, that is, a plume composed of fluid less dense than the ambient. For the cases where this criterion is not met, Boussinesq approximation may not be used and a density equation needs to be incorporated. El-Amin [6] introduced a numerical investigation of a non-Boussinesq, low-density gas jet (hydrogen) leaking into a high-density ambient (air). The integral models of jet fluxes are obtained and transformed into a set of ordinary differential equations of the mean centerline quantities. Therefore, mean quantities are obtained in addition to cross-stream velocity, Reynolds stresses, and turbulent Schmidt number. Furthermore, the normalized jet-feed material density and momentum flux density are correlated.

It is worth mentioning that theoretical developments and analysis of jet/plume theory were studied by a number of authors since 1950s (see, e.g., Morton et al. [20]; Morton [21]; Morton and Middleton [22]; Delichatsios [23]; Rooney and Linden [24]; Hunt and Kaye [25, 26]; Carlotti and Hunt [27]). Recently, Michaux and Vauquelin [28] developed analytical solutions for centerlines quantities of turbulent plumes rising from circular sources of positive buoyancy in a quiescent environment of uniform density for both Boussinesq and non-Boussinesq cases.

In this paper, semianalytical solution and theoretical analysis are developed for round hydrogen jet leaking into air based on Michaux and Vauquelin [28]. It is assumed that the rate of entrainment is a function of the plume centerline velocity and the ratio of the mean plume and ambient densities. Analytical expressions for all plume variables (radius, velocity, and density deficit) in terms of plume function for a given source parameter are derived.

#### 2. Mathematical Analysis and Similarity Solutions

##### 2.1. Governing Equations

Consider a vertical axisymmetric hydrogen-air buoyant jet resulting from a small-scale hydrogen leakage in the air with a finite circular source. Using cylindrical polar coordinates with the -axis vertical, the source is located at . The continuity, momentum, and concentration equations in cylindrical coordinate system (Figure 1) for the steady vertical axisymmetric buoyant free jet can be written as [29]:

where *U* is the mean streamwise velocity, and *V* is the mean cross-stream velocity, and *C* is the hydrogen concentration (mass fraction). The overbar denotes the time-averaged quantities, *u, v* are the components of velocity fluctuations in *z*, *r *directions, respectively, *c* is the concentration fluctuation, and is the mixture density.

On the other hand, from the experimental observations, the equations for the vertical velocity, density deficiency, and mass fraction profiles, assuming that the hydrogen-air mixture behaves as an ideal gas, are as follows (Fisher et al. [30], Hussein et al. [31], Shabbir and George [32], and Schefer et al. [9, 10]):

where and are the mean velocity and mean density at any point of the jet body; and are the centerline velocity and density. is the jet/plume width which increases linearly with *z, * is the momentum spread rate of the jet. is the virtual origin, which is the distance above/below the orifice where the flow appears to originate. The experimentally measured spread rate varies in the range 0.1–0.13. The buoyancy spreading factor expresses the ratio of spreading rates between the velocity and density deficiency profiles. The corresponding streamwise concentration for the axisymmetric hydrogen-air, free jet as detected experimentally by Schefer et al. [10] is given as . In general, the spread rate for the concentration is given in the formula . In the work of Schefer et al. [10], the momentum spread rate for the case of hydrogen jet was estimated as , from which one can find , and . It is well known that , that is, velocity and density spread at different rates.

##### 2.2. Similarity Solutions

Integrating the continuity (2.1) radially givesSince is negligible for , then integrating (2.1) for givesThis implies that

where denotes the inflow velocity at the plume edge which is known as the entrainment velocity. Therefore, we have

This equation indicates that the increase in plume volume flux is supplied by a radial influx from the far field which in turn implies a flow across the plume boundary . Batchelor [33] concluded that a vigorous entrainment of the ambient will be obtained as the density ratio tends to unity, . While as the density ratio tends to zero, , the entrainment falls to zero. Between these two limits, there will be a smooth transition of entrainment pattern. The experiments by Ricou and Spalding [34] suggest that the entrainment velocity may be obtained using the following formulawhere is the entrainment coefficient.

Also, Morton [35] assumed that the rate of entrainment into a strongly buoyant plume is a function of both density ratio and Reynolds stresses which have a magnitude proportional to . Therefore, the local entrainment velocity may be obtained as , which has also been suggested by Thomas [36], Steward [37], and Townsend [38]. Therefore, (2.11) can be written in the form

For calculating the momentum flux, let us integrate (2.2) with respect to , from to , noting that and , we get

Similarly, for concentration flux, integrating (2.3) with respect to from to , noting also that and , we get

This equation may be equivalent to the buoyancy flux equation which can be written in the following form [19]: Substituting (2.4)–(2.6) into (2.16), one obtains Substituting (2.5), (2.6), and (2.17) into (2.13), (2.14), and (2.15) gives

In homogeneous surroundings of density , the density deficit flux, (2.17), is equivalent to

which has the dimension of a buoyancy flux.

Recently, Michaux and Vauquelin [28] developed analytical solutions for centerlines quantities of turbulent plumes. Now, let us introduce a modified radius and a dimensionless density deficit [28]:

Therefore, (2.17)–(2.19) can be rewritten in the following form:

In the previous work of El-Amin and Kanayama [5] and El-Amin [6], they developed the similarity formulation and numerical solutions of the centerline quantities such as velocity and concentration.

In the current work, we follow the work of Michaux and Vauquelin [28] to obtain analytical/semianalytical solutions for centerline plume quantities. Using , and , the constants and ; exponents *m*, *n*, *p* can be determined to obtain similarity solutions for , and as

in these expressions may be replaced by to adapt solutions at near-source region and is virtual origin.

##### 2.3. Plume Function and Source Parameter

Now, let us introduce the plume function as [28]:

At source , , corresponding to the source parameter initially introduced by Morton [21] and defined as [26]:

where *Q*, *M*, and are the initial values of specific mass flux, specific momentum flux, and specific buoyancy flux, respectively, defined as

is the inlet diameter, is velocity at source, and is the difference in density between the receiving fluid and the fluid being discharged. Based on source parameter value Morton and Middleton [22] have categorized plumes with positive buoyancy as simple (pure) plume (), forced plume (), and lazy plume (). Other possibilities (Hunt and Kaye [25]) for , flow is pure jet without buoyancy, and, for , flow is weak fountains (negative buoyancy).

For hydrogen-air plume the source parameter is (of order 10^{−4}), so, based on Morton et al. [20] classification, it is a forced plume.

Equation (2.22) can be rewritten in terms of as follow:

Using plume function , (2.24), and (2.27)–(2.29), we may write

One can deduce that for *, * increases monotonically with height and tends asymptotically toward unity.

Using (2.27) and (2.30), one may get
Integrating this equation subject to the source conditions, one obtains
therefore,
Similarly, we can find that
therefore,
Finally, using (2.24), (2.32), and (2.34), we get
or
Now, is a function of *z* to relate each plume variable to *z*. Substituting (2.32) into (2.30), one gets
where , and is the characteristic length defined from initial plume condition. For the case under consideration, we find .

Integrating (2.38), one obtains

The above integration function has no explicit form, so Michaux and Vauquelin [28] computed and tabulated the integral function for several values of *X*. Unfortunately, they did not provide very small values as we find in this investigation when the source parameter is (of order 10^{−4}). Alternatively, we can write the integration function as

where is the hypergeometric function defined by

But is very small (of order 10^{−4}) for the case of hydrogen-air plume, so two terms of the above series can approximate the function. Thus,

This may lead to

According to this relation, , therefore,

From (2.43) and (2.44), one gets

Again from (2.43), we can determine values of from given by (2.44). It is clear that decreases as increases and increases as *z* increases. Thus, decreases as z increases and the maximum value of is located at , which is equal to . Therefore, using (2.36), one can find , and the maximum velocity is .

Finally, by combining (2.4)–(2.6) with (2.33), (2.36), and (2.37), one can obtain analytical expressions for vertical velocity, density deficiency, and mass fraction for hydrogen-air mixture in terms of the universal variable as follows:

#### 3. Conclusion

This paper introduces the reader to a set of features of hydrogen-air plume, which is very important to assess the potential hazard resulting from hydrogen sources upon leakage into the ambient atmosphere. Throughout this work, we derived profiles of the mean quantities for a turbulent hydrogen-air plume. These mean quantities, such as plume radius, velocity, and density deficit, are expressed in terms of the plume function for a given source parameter. These quantities are determined by integral relations and by analysis using similarity variables. Therefore, mean quantities are expressed solely in terms of the plume function and the source parameter. The plume function is valid for a range of small values of as required for hydrogen-air plume. The hypergeometric function is exploited, and the profiles of the mean quantities are obtained. These results may be generalized and extended in a future work to cover more complex flow types found in a prober accident of hydrogen leaks such as leakage in hydrogen station.