D68 /
HorizontalBuoyantJetFor the numerical simulation of turbulent buoyant jets, normally there are two methods available: integral methods and simulation by CFD. Integral methods are based on the basic laws of conservation of mass, momentum, energy, and species concentration. To integrate the partial differential equations, empirical profile shapes (i.e. Gaussian distribution) for velocity, temperature and concentration are assumed. The entrainment of ambient gas into the turbulent buoyant jet is also needed to be assumed to close the equation system [5-8]. Pantokratoras[9] modified the integral Fan-Brooks model[10] to calculate the horizontal penetration of inclined thermal buoyant water jets, and the modified model predictions are in good agreement with the trajectory measurements. Jirka[11-12] formulated an integral model, named CorJet, for turbulent buoyant jets (round and planar) in unbounded stratified flows, i.e., the pure jet, the pure plume, the pure wake. Based on the Boussinesq approximation, CorJet integral model appears to provide a sound, accurate and reliable representation of complex buoyant jet physics. However integral methods are difficult to extend to more complex turbulent buoyant flows which are essentially three dimensional and often obstacles and walls interact with the flow. It is very difficult to prescribe the profile shapes and to relate the entrainment rate to all of the local parameters. Therefore CFD methods become more popular which do not make assumptions for the profile and entrainment but yield them as a part of the solution. However, if not applying DNS, the turbulence model in the engineering CFD code is one of the most significant factors to affect the simulations of turbulent buoyant jets and plumes. Most of the turbulence models were developed and tested for non-buoyant flows. Turbulence models including buoyant effects [13] were reviewed by Hossain and Rodi[14]. In recent studies models with buoyancy related modification were evaluated, assessing in particular buoyancy effect on the production and dissipation of the turbulent kinetic energy in buoyant plumes[15-17]. The standard models were suspected to seriously underpredict the spreading rate of vertical buoyant jets and to overpredict the entrainment of horizontal flows. Comparing to the vertical buoyant jet, much fewer data and calculations for the study of horizontal turbulent buoyant jets can be found in the open literatures[18-23]. Almost all the models and calculations available are based on Boussinesq approximation, which the environmental and atmospheric engineers are more interested in. What we are concerning is hydrogen, helium, or steam released into air, which involves a large density difference between the jet and ambient, making the Boussinesq approximation invalid. In the table below, integral models of horizontal buoyant round and plane jets are introduced. The system of ordinary differential equations can be solved by 4th order Runga-Kutta method to obtain the horizontal buoyant jet trajectory, the velocity, the density, and the tracer concentration. For small density difference cases, the Boussinesq integral model can obtain reasonable results as the CorJet model shows in the table. For hydrogen safety analysis, the large density difference will violate the Boussinesq approximation.
Boussinesq integral model for horizontal buoyant plane jet. |