Journal of Applied Mathematics

Volume 2012, Article ID 710172, 12 pages

http://dx.doi.org/10.1155/2012/710172

## Theoretical Analysis and Semianalytical Solutions for a Turbulent Buoyant Hydrogen-Air Jet

^{1}Computational Transport Phenomena Laboratory (CTPL), Physical Sciences and Engineering Division (PSE), King Abdullah University of Science and Technology (KAUST), Thuwal, Jeddah 23955-6900, Saudi Arabia^{2}Department of Mathematics, Aswan Faculty of Science, South Valley University, Aswan 81528, Egypt

Received 1 October 2011; Revised 16 November 2011; Accepted 28 November 2011

Academic Editor: Hiroshi Kanayama

Copyright © 2012 M. F. El-Amin et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### 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.

#### References

- M. F. El-Amin and H. Kanayama, “Boundary layer theory approach to the concentration layer adjacent to a ceiling wall at impinging region of a hydrogen leakage,”
*International Journal of Hydrogen Energy*, vol. 33, no. 21, pp. 6393–6400, 2008. View at Publisher · View at Google Scholar · View at Scopus - M. F. El-Amin, M. Inoue, and H. Kanayama, “Boundary layer theory approach to the concentration layer adjacent to a ceiling wall of a hydrogen leakage: far region,”
*International Journal of Hydrogen Energy*, vol. 33, no. 24, pp. 7642–7647, 2008. View at Publisher · View at Google Scholar · View at Scopus - M. F. El-Amin and H. Kanayama, “Boundary layer theory approach to the concentration layer adjacent to the ceiling wall of a hydrogen leakage: axisymmetric impinging and far regions,”
*International Journal of Hydrogen Energy*, vol. 34, no. 3, pp. 1620–1626, 2009. View at Publisher · View at Google Scholar · View at Scopus - M. F. El-Amin and H. Kanayama, “Integral solutions for selected turbulent quantities of small-scale hydrogen leakage: a non-buoyant jet or momentum-dominated buoyant jet regime,”
*International Journal of Hydrogen Energy*, vol. 34, no. 3, pp. 1607–1612, 2009. View at Publisher · View at Google Scholar · View at Scopus - M. F. El-Amin and H. Kanayama, “Similarity consideration of the buoyant jet resulting from hydrogen leakage,”
*International Journal of Hydrogen Energy*, vol. 34, no. 14, pp. 5803–5809, 2009. View at Publisher · View at Google Scholar · View at Scopus - M. F. El-Amin, “Non-Boussinesq turbulent buoyant jet resulting from hydrogen leakage in air,”
*International Journal of Hydrogen Energy*, vol. 34, no. 18, pp. 7873–7882, 2009. View at Publisher · View at Google Scholar · View at Scopus - M. Sánchez-Sanz, M. Rosales, and A. L. Sánchez, “The hydrogen laminar jet,”
*International Journal of Hydrogen Energy*, vol. 35, no. 8, pp. 3919–3927, 2010. View at Publisher · View at Google Scholar · View at Scopus - W. G. Houf and R. W. Schefer, “Analytical and experimental investigation of small-scale unintended releases of hydrogen,”
*International Journal of Hydrogen Energy*, vol. 33, no. 4, pp. 1435–1444, 2008. View at Publisher · View at Google Scholar · View at Scopus - R. W. Sehefer, W. G. Houf, and T. C. Williams, “Investigation of small-scale unintended releases of hydrogen: momentum-dominated regime,”
*International Journal of Hydrogen Energy*, vol. 33, no. 21, pp. 6373–6384, 2008. View at Publisher · View at Google Scholar · View at Scopus - R. W. Schefer, W. G. Houf, and T. C. Williams, “Investigation of small-scale unintended releases of hydrogen: buoyancy effects,”
*International Journal of Hydrogen Energy*, vol. 33, no. 17, pp. 4702–4712, 2008. View at Publisher · View at Google Scholar · View at Scopus - K. Matsuura, “Effects of the geometrical configuration of a ventilation system on leaking hydrogen dispersion and accumulation,”
*International Journal of Hydrogen Energy*, vol. 34, no. 24, pp. 9869–9878, 2009. View at Publisher · View at Google Scholar · View at Scopus - K. Matsuura, M. Nakano, and J. Ishimoto, “Forced ventilation for sensing-based risk mitigation of leaking hydrogen in a partially open space,”
*International Journal of Hydrogen Energy*, vol. 35, no. 10, pp. 4776–4786, 2010. View at Publisher · View at Google Scholar · View at Scopus - K. Matsuura, M. Nakano, and J. Ishimoto, “The sensing-based adaptive risk mitigation of leaking hydrogen in a partially open space,”
*International Journal of Hydrogen Energy*, vol. 34, no. 20, pp. 8770–8782, 2009. View at Publisher · View at Google Scholar · View at Scopus - K. Matsuura, H. Kanayama, H. Tsukikawa, and M. Inoue, “Numerical simulation of leaking hydrogen dispersion behavior in a partially open space,”
*International Journal of Hydrogen Energy*, vol. 33, no. 1, pp. 240–247, 2008. View at Publisher · View at Google Scholar · View at Scopus - S. Kikukawa, “Consequence analysis and safety verification of hydrogen fueling stations using CFD simulation,”
*International Journal of Hydrogen Energy*, vol. 33, no. 4, pp. 1425–1434, 2008. View at Publisher · View at Google Scholar · View at Scopus - V. Agarant, Z. Cheng, and A. Tchouvelev, “CFD modeling of hydrogen releases and dispersion in hydrogen energy station,” in
*Proceedings of the 15th World Hydrogen Energy Conference*, 2004. - M. R. Swain and M. N. Swain, “Passive ventilation systems for the safe use of hydrogen,”
*International Journal of Hydrogen Energy*, vol. 21, no. 10, pp. 823–835, 1996. View at Publisher · View at Google Scholar · View at Scopus - M. R. Swain, P. Filoso, E. S. Grilliot, and M. N. Swain, “Hydrogen leakage into simple geometric enclosures,”
*International Journal of Hydrogen Energy*, vol. 28, no. 2, pp. 229–248, 2003. View at Publisher · View at Google Scholar · View at Scopus - P. F. Crapper and W. D. Baines, “Non Boussinesq forced plumes,”
*Atmospheric Environment*, vol. 11, pp. 415–420, 1977. View at Publisher · View at Google Scholar · View at Scopus - B. R. Morton, G. Taylor, and J. S. Turner, “Turbulent gravitational convection from maintained and instantaneous sources,”
*Proceedings of the Royal Society A*, vol. 234, pp. 1–23, 1956. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - B. R. Morton, “Forced plumes,”
*Journal of Fluid Mechanics*, vol. 5, pp. 151–163, 1959. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - B. R. Morton and J. Middleton, “Scaling diagram for forced plumes,”
*The Journal of Fluid Mechanics*, vol. 58, pp. 165–176, 1973. View at Publisher · View at Google Scholar - M. A. Delichatsios, “Strong turbulent buoyant plumes. I. Similarity,”
*Combustion Science and Technology*, vol. 24, pp. 191–195, 1981. View at Google Scholar - G. G. Rooney and P. F. Linden, “Similarity considerations for non-Boussinesq plumes in an unstratified environment,”
*The Journal of Fluid Mechanics*, vol. 318, pp. 237–250, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - G. R. Hunt and N. G. Kaye, “Virtual origin correction for lazy turbulent plumes,”
*The Journal of Fluid Mechanics*, vol. 435, pp. 377–396, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - G. R. Hunt and N. B. Kaye, “Lazy plumes,”
*The Journal of Fluid Mechanics*, vol. 533, pp. 329–338, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - P. Carlotti and G. R. Hunt, “Analytical solutions for turbulent non-Boussinesq plumes,”
*The Journal of Fluid Mechanics*, vol. 538, pp. 343–359, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - G. Michaux and O. Vauquelin, “Solutions for turbulent buoyant plumes rising from circular sources,”
*Physics of Fluids*, vol. 20, no. 6, Article ID 066601, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - M. F. El-Amin and S. Sun, “Turbulent buoyant jet of a low-density gas leaks into a high-density ambient: hydrogen leakage in air,” in
*Advanced Topics in Mass Transfer*, M. F. El-Amin, Ed., INTECH, 2011. View at Google Scholar - H. B. Fischer, E. J. List, R. C. Y. Koh, J. Imberger, and N. H. Brooks,
*Mixing in Inland and Coastal Waters*, Academic Press, 1979. - H. J. Hussein, S. P. Capp, and W. K. George, “Velocity measurements in a high-Reynolds-number, momentum-conserving, axisymmetric, turbulent jet,”
*The Journal of Fluid Mechanics*, vol. 258, pp. 31–75, 1994. View at Publisher · View at Google Scholar - A. Shabbir and W. K. George, “Experiments on a round turbulent buoyant plume,”
*The Journal of Fluid Mechanics*, vol. 275, pp. 1–32, 1994. View at Publisher · View at Google Scholar - G. K. Batchelor, “Heat convection and buoyancy effects in fluids,”
*Quarterly Journal of the Royal Meteorological Society*, vol. 80, pp. 339–358, 1954. View at Google Scholar - F. P. Ricou and D. B. Spalding, “Measurements of entrainment by axisymmetrical
turbulent jets,”
*The Journal of Fluid Mechanics*, vol. 8, pp. 21–32, 1961. View at Google Scholar · View at Zentralblatt MATH - B. R. Morton, “Modelling fire plumes,” in
*Proceedings of the 10th International Symposium on Combustion*, pp. 973–982, 1965. - P. H. Thomas, “The size of flames from natural fires,” in
*Proceedings of the 9th Symposium on Combustion*, pp. 844–859, 1963. - F. R. Steward, “Prediction of the height of turbulent diffusion buoyant flames,”
*Combustion Science and Technology*, vol. 2, pp. 203–212, 1970. View at Publisher · View at Google Scholar - A. A. Townsend, “The mechanism of entrainment in free turbulent flows,”
*The Journal of Fluid Mechanics*, vol. 26, pp. 689–715, 1966. View at Publisher · View at Google Scholar