ON THE NUMERICAL MODELING OF BUOYANCY- DOMINATED TURBULENT DIFFUSION FLAMES BY USING URANS WITH A k ε TURBULENCE MODEL - PDF

Description
MCS 7 Chia Laguna, Cagliari, Sardinia, Italy, September -5, 2 ON THE NUMERICAL MODELING OF BUOYANCY- DOMINATED TURBULENT DIFFUSION FLAMES BY USING URANS WITH A k ε TURBULENCE MODEL B. Sapa*/**, P. Plion**,

Please download to get full document.

View again

of 12
All materials on our website are shared by users. If you have any questions about copyright issues, please report us to resolve them. We are always happy to assist you.
Information
Category:

Entertainment & Humor

Publish on:

Views: 25 | Pages: 12

Extension: PDF | Download: 0

Share
Transcript
MCS 7 Chia Laguna, Cagliari, Sardinia, Italy, September -5, 2 ON THE NUMERICAL MODELING OF BUOYANCY- DOMINATED TURBULENT DIFFUSION FLAMES BY USING URANS WITH A k ε TURBULENCE MODEL B. Sapa*/**, P. Plion**, L. Gay**, F. Nmira**, H. Y. Wang* * Institut P, Fluides, Thermique, Combustion, CNRS, ENSMA, Université de Poitiers, BP 49 F8696 Futuroscope Chasseneuil Cedex, France **EDF R\&D, 6 quai Watier, 784 Chatou, France Abstract This work concerns the modelling of buoyancy-dominated turbulent diffusion flames through an URANS approach in the software. It is found that URANS successfully predicts the important features of highly oscillating buoyancy-induced flows (recirculation zone, narrowing and broadening of the flame). Also, the mean velocity, temperature and chemical species are correctly reproduced by URANS. While, the standard k ε model should be improved to model the turbulent mixing near the burning wall for a low Froude number flame.. Introduction A large number of studies [-7] have been performed on the characteristics of buoyant diffusion flames above a horizontal burning surface. Becker et al. [2] have demonstrated that, in the natural convection limit, as the Froude number decreases, coherent structures appear in the reactive zone, the flame presenting a pronounced instability due to buoyancy. Besides, coupling between pyrolysis rate and heat transfer becomes rather complex when fire is confined with severe ventilation condition, where the phenomenon such as extinction and reinflamation or ghost flames appear. Thus, a low Froude number flame in free conditions has already been largely studied [3-7] by using a porous burner. The detailed measurements of the flame and flow structures have been obtained with medium-scale fires [4]. For a low Froude number flame, shear-stresses between hot combustion products and fresh air make the flow unstable and amplify oscillations near the fire base (Weckman [6]) due to air entrainment variation and flame flicker, inducing large eddy structures corresponding to hot gases puffs burning. Fresh air entrained by these vortexes feeds the flame with oxygen and cools the smoke influencing natural convection and then air entrainment. Computational Fluid Dynamics are expected to improve understanding and prediction, by its local description of such buoyancy-controlled flow, of the gas composition and thermal exchanges. Several authors [7] have shown the limit of the use of Reynolds Averaged Navier-Stokes (RANS) k ε turbulence model to simulate the highly oscillating buoyancy-induced flow. In order to capture the buoyancy-controlled flow at industrial scale, Large Eddy Simulation (LES) and Unsteady Reynolds Average Navier-Stokes (URANS) method are currently employed. The accuracy of LES depends on filtering scale related to mesh size, and a few millions cells is usually needed for a good description of combustion scales. URANS method is based on Reynolds average equations and the definition of ensemble (or time, under ergodic hypothesis) averaging implies that scales smaller than a characteristic time are modelled and the other are solved through mean flow. The solution is assumed to be an averaged one in which turbulent effects are considered through turbulence models. In general, URANS allows less refined meshes than LES can do and time calculation becomes reasonable for cells width adapted to the modelling. With this main motivation, EDF R&D has developed an URANS approach in a CFD freeware [8] for industry application such as gas turbines, coal furnaces, compressible flows, gas-particles interactions or electric arcs. It is found that a fluctuating, laminar, diffusion flame at the fire base, and a turbulent, intermittent one as it evolves upward along the fire plume are better described by URANS approach. The mean temperature, velocity and main chemical species, provided from URANS, are also compared with the experimental data, and a good agreement is observed. 2. Physical Modelling The simulations were performed using the software [8], in which combustion is modelled by a turbulent diffusion flame model using a presumed probability density function. Chemistry is described by an infinitely fast global reaction at a constant stœchiometric ratio. The starting point of the analysis is the set of three-dimensional elliptic, reacting flow equations that governs the phenomena of interest here. This set consists of the conservation equations of mass, momentum, energy and species. Density is calculated by the perfect gas law. Two equations are solved for turbulent kinetic energy k and its dissipation rate ε [9]. A finite volume technique is used to discretize the partial differential equations. The precise formulation of the differential equations describing the model and the numerical technique can be found elsewhere [8], and will not be repeated here. In the following section, only URANS approach and radiation associated with soot formation, developed by the authors in, are briefly described. 2. Unsteady RANS approach URANS is based on a Favre meaning operation which decomposes a full flow field Φ (x, t) into a mean component Φ (x, t) (ensemble average) and a fluctuation Φ (x, t). Meaning the exact equations of mass and momentum, averaged equations are obtained: d + div u = Γ () dt du + div ( u.u ) = grad p + div τ + Su (2) dt In, thermal expansion is not taken into account in SIMPLEC algorithm for Navier-Stokes system resolution due to neglecting the mass accumulation term in continuity equation. The new procedure, named as URANS, accounts for the thermal expansion in continuity equation through density-velocity coupling resolution. Based on a low Mach algorithm, this approach is developed still from a two-step scheme with a prediction step which considers momentum equation treated with an explicit pressure and uses mass equation for density derivative versus time. du + u [ Γ div( u )] + div ( u.u ) = grad p + div τ + Su (3) dt It gives a predicted velocity u*, and a correction step which considers mass equation () and time-discretizised momentum correction equation (δp = p - p n ): d u * t = grad δp u u = grad δp (4) dt Replacing corrected velocity in divergence of equation (4) by adequate form of equation (), an equation for pressure increment is obtained: t * d * div grad δp = div u div u = Γ div u (5) dt For a density dependant of one scalar (φ), from non-conservative forms of mass and scalar equation: d + div u = Γ (6) dt dφ = div( D grad φ) + S (7) φ dt By multiplying Eq.(6) and making density particular derivative appearing in Eq.(7), the velocity divergence can be expressed as: div u = Γ [ div( D grad φ) + S ] 2 (8) φ φ and is used in correction step of Navier-Stokes resolution : t div grad δp = Γ div( D grad ) S 2 φ + φ * [ ] div u φ (9) Corrected velocity is deduced from predicted one and pressure increment through: * t u u = grad δp () 2.2 Radiative heat transfer In a heavily sooting flame as fire, as the radiation spectrum of soot is continuous, it is possible to assume that the gas behaves as a gray medium. The spectral dependence is then lumped into one absorption coefficient, κ, and the following radiative transfer equation without scattering is solved. r r σ Ω + κ = κ T 4. I I () π The overall absorption coefficient is calculated from a grey gas Modak s regression [] in function of the temperature and the mixture concentration of soot and gas (CO 2 and H 2 O). For this, the soot formation and its oxidation are incorporated into a turbulent flow calculation in two convection-diffusion equations for the soot number density, n, and soot volume fraction, f v, expressed as follows, n ( µ + u i n) t n = ω& n (2) t x j x j Sct x j f t v ( u i + f x j v ) x j µ t Sc t f x v j = ω & fv (3) The local formation rate, ω& n, for the soot number density, and ω& fv, for the soot volume fraction is calculated as, C (36 ) W T X e T / T T n n 4 / 3 2 / 2 β / 2 2 π ox = C N α ω& (4) n α F / 3 N soot f v / 3 / 3 π 2 / 3 ω = C δ 2 / 2 γ v /T C / 2 / 3 2 / 3 ox α + /T γ (36 ) W n f & fv α F / 3 Fe T n f v C T X e T T soot X (5) soot soot where N ( /kmol) is the Avogadro s number, and soot, the soot density, given as 8 kg/m 3. The empirical constants, such as C α ( m 3.kg -2.K -/2.s - ), C β (.3 7 m 3.K /2.s - ), C δ (44 kg), C γ (. m 3.kg -2/3.K -/2.s - ), T α (4.6 4 K) and T γ (.26 4 K) were experimentally determined by Moss et al. [] for a range of fuel. The first term in Eq.(4) represents the rate of particle nucleation and the second the coagulation of soot. The second term in Eq.(5) represents the surface growth of soot which contained a linear dependence on aerosol surface area. The term, ω ox (kg/m 2 s), on the right hand side of Eqs.(4-5) corresponds to soot oxidation which is evaluated from the rate for oxidation of pyrolitic graphite by O 2 [2]. 3. Results and discussion The propane diffusion flames are chosen for comparison between predictions and experiments by using URANS with a turbulence model, k ε. The free boundaries are located sufficiently far away from the fire to minimize a numerical perturbation. The 3D computational domain are 2.4 m x 2.4 m x.5 m in the x, y and z directions, respectively. / 3 Here we concentrate on a comparison of the simulation predictions with the experimental data [4] for a pool diameter, D, of.3 m. The flame is stabilized on a horizontal circular porous burner, providing a heat release rate (HRR) in a range from 6, 23 to 38 kw. Temperatures were obtained by thermocouples and velocities were determined using a two-component laser Doppler velocimetry (LDV) system. Only the experimental data in the symmetrical plane is available, and the numerical results are examined in the plane, x-z. Pool fire Figure. Schematic diagram of the release examined experimentally including the free boundary and the coordinate system. z y Free boundaries Wall x In order to have a true predictive capability, grid refinement studies were performed, by using a series (Mesh -4) of a typical 3D grid points containing, 24(x) x 24(y) x 28(z), 3(x) x 3(y) x 35(z), 9(x) x 9(y) x 6(z) and 3(x) x 3(y) x 9(z) cells, respectively. For a fire plume, the characteristic length scale can be estimated from a given HRR [3], and the HRR considered here is higher than 2 kw, thus the characteristic length is in an order of. m. The extra grid points are added in the fire region where a strongly stratified layer is developed. At the fire source base, the length scale,, is of the order of 3, 2, to.5 cm by using the four grid systems (Mesh -4), respectively. Cell sizes are uniform near the fire source and stretch to about 5 cm far away from the pool-fire region at the free boundary. About 3 days CPU time are needed for a real minute simulation on a 42 processors. Influence of number of grid cells on the predicted temperature and velocity is checked in Figures 2 and 3 for a heat release rate of 23 kw. The changes between the coarse Mesh and fine Mesh 3 in the calculated plots of temperature and velocity above the fire source base are higher than 6%. A further reduction in the grid size (Mesh 4) results in a significant reduction in the time step ( t. s) for satisfying the CFL stability condition. An adequate resolution of the fire plume in large-scale can be achieved with a spatial resolution of about. m. Finally, the Mesh 3 with 9(x) x 9(y) x 6(z), are employed because it gives the best trade-off between accuracy and cost for the present purpose. T (K) 24 Mesh 2 Mesh 2 Mesh 3 Mesh Figure 2. Influence of number of grid cells on the predicted temperature w (m/s) Figure 3. Influence of number of grid cells on the predicted vertical velocity Based on the grid system Mesh 3, effects of soot formation and gas radiation on the mean temperature T, are shown in Fig.4, as compared with the experimental data. It can be seen that in the plume zone, the temperature peak decrease progressively by accounting for the radiation and soot formation. The calculated maximum values are in good agreement with experiments, while the position of the flame is far away the measured one. The use of the URANS coupled with one-step irreversible chemical reaction is not sufficient to improve the prediction of the flame structure, and the lateral spread of the diffusion flame is usually under-predicted. T (K) adiabatic gas radiation soot and gas radiation exp Figure 4. Effects of soot formation and gas radiation on the mean temperature The calculated characteristics of the instantaneous temperature output on the middle plane (x-z, y=.2 m) is plotted in Figs.5 and 6 by using RANS and URANS, respectively. The RANS approach generates a stable diffusion flame, and the flow behaves as a boundary layer type. While three regions of the pool fires are reproduced by URANS approach. At the flame base, immediately above the burner surface, is the persistent flame zone, which appears to have a constant shape and structure. This zone is followed by the intermittent zone that has a fluctuating character due to buoyant instability. In this region, ambient air is sucked into the flame from a downstream location, creating a large recirculating zone above the horizontal burner. The last region is the plume zone that is non reacting and turbulent, giving a good idea of the flow development and the boundary layer broadening. In general, the large scale structure can be completely described when the characteristic length is spanned by roughly ten computational cells. The greater temporal and spatial fidelity can be considered to be extracted by using URANS with Mesh 3 for an application to three-dimensional pool-like fire. Figures 5. Instantaneous temperature output from RANS Figures 6. Instantaneous temperature output from URANS URANS approach can provide a detailed description of the oscillatory behaviour for the buoyancy-driven fire. Determination of the oscillatory behaviour is also important for the spread of fire and hot gases. The simulated time histories of the temperature in the combusting portions of the plume at z =.3 m, for a heat release of 23 kw are depicted in Figure 7. The time traces indicate that a large oscillatory fluctuation in the temperature is captured by the time-resolved calculations. This behaviour could manifest itself as the puffs of the plume, and the amplitude of the oscillation of the temperature is as large as about 5 C. The spectral analysis (cf. Fig.8) from prediction suggest that the buoyancy-driven flow produces oscillatory behaviour in cyclic fashion, with a frequency between 2 and 3 Hz, which is commonly observed in experiments and agrees with Zukoski s correlation [5]. Even if the grid size in Mesh 3 is sufficiently small, an incompressible approach (RANS) or a compressible one (URANS) influences significantly the profiles of the mean temperature, T, and the longitudinal velocity, w, along the axial distance, as compared with the experimental data in Figs.9 and. Profiles of the longitudinal velocity fluctuation, w, and the transversal one, u, obtained from the two approaches RANS and URANS are also compared with the experimental data in Figs. and 2. The velocity fluctuation predicted by the RANS is unrealistically low at the fire base, indicating that this approach can not reproduce the expanded zone of hot gases which coincides with the periodic formation of large coherent structures (cf. Fig.6). The significant over-prediction of the mean temperature and velocity (cf. Figs.9 and ) with RANS is certainly due to underestimation of the resolved buoyant instability. Subsequently, the flame behaves as a laminar type, causing a higher flame 24 T (K) t (s) Figure 7. Time traces of a large oscillatory fluctuation in the temperature f (Hz) Figure 8. Spectral evolution of the oscillatory behaviour in the temperature T (K) URANS RANS Figure 9. Comparison between the measured and predicted mean temperature from RANS and URANS.8.6 URANS RANS Figure. Comparison between the measured and predicted longitudinal velocity fluctuation from RANS and URANS w' ' (m/s) w (m/s) RANS URANS Figure. Comparison between the measured and predicted mean longitudinal velocity from RANS and URANS u' ' (m/s) RANS URANS Figure 2. Comparison between the measured and predicted mean transversal velocity fluctuation from RANS and URANS temperature as compared to the experiment. At the fire region, the underprediction of u and w highlights the serious shortcoming in the RANS approach [7] which fails to model the turbulent mixing there even with a fine Mesh 3. For the buoyancy-controlled fire, the oscillation of the flame from buoyant instability is initially taking place near the edge of the horizontal burner surface, where turbulence is weak due to a low fuel injection velocity. The high levels of the velocity fluctuation (u and w ) due to the presence of large-scale coherent structures are correctly predicted by URANS. The profiles of the mean temperature, T and the longitudinal velocity along the axial distance for the HRRs of 6 and 38 kw are also compared with the experimental data in Figs.3-6. In general terms, the comparison between URANS prediction and experiment is reasonably good, although the longitudinal velocity is clearly underestimated. Globally, the approach of URANS is capable of reproducing the mechanism generating the buoyant instability present in the early development of the flame and the transition to turbulence far away from the fire source (z .5 m) T (K) 8 w (m/s) Figure 3. Comparison between the measured and predicted mean temperature for HRR of 6 kw Figure 4. Comparison between the measured and predicted vertical velocity for HRR of 6 kw T (K) 8 w (m/s) Figure 5. Comparison between the measured and predicted mean temperature for HRR of 38 kw.5.5 Figure 6. Comparison between the measured and predicted vertical velocity for HRR of 38 kw Figures 7-22 show the profiles for the mean mass fraction of the chemical species C 3 H 8, CO 2 and O 2 along the axial distance. The experimental data show that just above the burning wall, a fuel-rich zone develops and sharply disappears as we move away from the pool fire surface. A strong mixing between the fuel vapour and air probably results from a convective type of mixing rather than from a diffusive one, which gives to the reaction zone a premixed rather than diffusive character. Evident in Figs.7, 9 and 2 is the good agreement between the predicted and the measured values of O 2. As shown in Fig.2, the quality of the agreement between the URANS prediction and experiment for fuel mass fraction deteriorates with an increase of the HRR. The URANS simulation suggests that the reacting flowfield is displaced in a stronger flapping manner. As we move away from the burning surface, the fuel and CO 2 concentrations seem significantly over-predicted. This divergence between prediction and experiment is certainly due to the standard k ε model which fails to model correctly the turbulent mixing near the reacting region. Mass fraction C 3 H 8 O 2 Mass fraction CO Figure 7. Comparison between the measured and predicted mean fuel and oxygen concentrations for HRR of 6 kw.5.5 Figure 8. Comparison between the measured and predicted CO 2 concentration for HRR of 6 kw.4.2 Mass fraction.3.2. C 3 H 8 O 2 Mass fraction CO Figure 9. Comparison between the measured and predicted mean fuel and oxygen concentrations for HRR of 23 kw.5.5 Figure 2. Comparison between the measured and predicted CO 2 concentration for HRR of 23 kw .4.2 Mass fraction.3.2. C 3 H 8 O 2 Mass fraction CO Figure 2. Comparison between the measured and predicted mean fuel and oxygen concentrations for HRR of 38 kw.5.5 Figure 22. Comparison between the measured and predicted CO 2 concentration for HRR of 38 kw One of the practical interests of pool-like fire simulation is the prediction of air entrainment from outside of the plume boundary. The entrainment rate is de
Related Search
Similar documents
View more...
We Need Your Support
Thank you for visiting our website and your interest in our free products and services. We are nonprofit website to share and download documents. To the running of this website, we need your help to support us.

Thanks to everyone for your continued support.

No, Thanks