Numerical Simulation of Reacting Flows in Radiant Porous Burners

2005 ◽  
Author(s):  
Timothy W. Tong ◽  
Mohsen M. M. Abou-Ellail ◽  
Yuan Li ◽  
Karam R. Beshay

The present paper presents, numerical computations for flow, heat transfer and chemical reactions in an axisymmetric inert porous burner. The porous media re-radiate the heat absorbed from the gaseous combustion products by convection and conduction. In the present work, the porous burner species mass fraction source terms are computed from an ‘extended’ reaction mechanism, controlled by chemical kinetics of elementary reactions. The porous burner has mingled zones of porous/nonporous reacting flow, i.e. the porosity is not uniform over the entire domain. Therefore, it has to be included inside the partial derivatives of the transport governing equations. Finite-difference equations are obtained by formal integration over control volumes surrounding each grid node. Up-wind differencing is used to insure that the influence coefficients are always positive to reflect the real effect of neighboring nodes on a typical central node. Finite-difference equations are solved, iteratively, for U, V, p’ (pressure correction), enthalpy and species mass fractions, utilizing a fine grid of (80×60) nodes. The eighty grid nodes in the axial direction are needed to resolve the detailed structure of the thin reaction zone inside the porous media. The radial grid is extended inside the annular solid wall of the porous burner, to compute the wall temperature. The porous burner uses a premixed CH4-air mixture, while its radiating characteristics are computed numerically, using a four-flux radiation model. Sixteen species are included, namely CH4, CH3, CH2, CH, CH2O, CHO, CO, CO2, O2, O, OH, H2, H, H2O, H2O, H2O2, involving 49 chemical reaction equations. It was found that 1000 iterations are sufficient for complete conversion of the computed results with errors less than 0.1%. The computed temperature profiles of the gas and the solid show that, heat is conducted from downstream to the upstream of the reaction zone. Most stable species, such as H2O, CO2, H2, keep increasing inside the reaction zone staying appreciable in the combustion products. However, unstable products, such as HO2, H2O2 and CH3, first increase in the preheating region of the reaction zone, they are then consumed fast in the post-reaction zone of the porous burner. Therefore, it appears that their important function is only to help the chemical reactions continue to their inevitable completion of the more stable combustion products.

Author(s):  
Timothy Tong ◽  
Mohsen Abou-Ellail ◽  
Yuan Li ◽  
Karam R. Beshay

The present paper presents, numerical computations for flow, heat transfer and chemical reactions in an axisymmetric inert porous burner. The porous media re-radiate the heat absorbed from the gaseous combustion products by convection and conduction. In the present work, the porous burner species mass fraction source terms are computed from an ‘extended’ reaction mechanism, controlled by chemical kinetics of elementary reactions. The porous burner has mingled zones of porous/nonporous reacting flow, i.e. the porosity is not uniform over the entire domain. Therefore, it has to be included inside the partial derivatives of the transport governing equations. Finite-difference equations are obtained by formal integration over control volumes surrounding each grid node. Up-wind differencing is used to insure that the influence coefficients are always positive to reflect the real effect of neighboring nodes on a typical central node. Finite-difference equations are solved, iteratively, for U, V, p’ (pressure correction), enthalpy and species mass fractions, utilizing a grid of (60×40) nodes. The sixty grid nodes in the axial direction are needed to resolve the detailed structure of the thin reaction zone inside the porous media. The porous burner uses a premixed CH4-air mixture, while its radiating characteristics are computed numerically, using a four-flux radiation model. Sixteen species are included, namely CH4, CH3, CH2, CH, CH2O, CHO, CO, CO2, O2, O, OH, H2, H, H2O, HO2, H2O2, involving 49 chemical reaction equations. It was found that 900 iterations are sufficient for complete conversion of the computed results with errors less than 0.1%. The computed temperature profiles of the gas and the solid show that, heat is conducted from downstream to the upstream of the reaction zone. Most stable species, such as H2O, CO2, H2, keep increasing inside the reaction zone staying appreciable in the combustion products. However, unstable products, such as HO2, H2O2 and CH3, first increase in the preheating region of the reaction zone, they are then consumed fast in the post-reaction zone of the porous burner. Therefore, it appears that their important function is only to help the chemical reactions continue to their inevitable completion of the more stable combustion products.


2008 ◽  
Author(s):  
Timothy W. Tong ◽  
Mohsen M. Abou-Ellail ◽  
Yuan Li ◽  
Karam R. Beshay

The present paper is concerned with the numerical computation of flow, heat transfer and chemical reactions in porous burners. One of the important features of porous burners is their presumed low levels of nitrogen oxides. In the present work, the computed NOx is compared with similar conventional premixed burners and measured nitrogen oxides in porous burners. In order to accurately compute the nitrogen oxides levels in porous burners, both prompt and thermal NOx mechanisms are included. In the present work, the porous burner species mass fraction source terms are computed from an ‘extended’ reaction mechanism, controlled by chemical kinetics of elementary reactions. The porous burner has mingled zones of porous/nonporous reacting flow, i.e. the porosity is not uniform over the entire domain. Finite-volume equations are obtained by formal integration over control volumes surrounding each grid node. Up-wind differencing is used to ensure that the influence coefficients are always positive to reflect the real effect of neighboring nodes on a typical central node. Finite-difference equations are solved iteratively for velocity components, pressure correction, gas enthalpy, species mass fractions and solid matrix temperature. The grid used to solve the solid energy equation is extended inside the zero-porosity solid annular wall of the burner porous disk. This was found useful for computing the solid wall temperature with high accuracy. A two-dimensional, discrete-ordinate, model is used for the computation of thermal radiation emitted from the solid matrix. The porous burner uses a premixed CH4-air mixture, while its radiating characteristics are studied numerically under equivalence ratio ranging from 0.5 to 0.8. Twenty-one species are included, involving 55 chemical reactions. The computed solid wall temperature profiles are compared with experimental data of similar porous burners. The obtained agreement is fairly good. The present numerical results show that as the equivalent ratio decreases, the reaction zone moves downstream. Moreover, as the flame speed increases, the NOx mole fraction increases. Some reacting species, such as H2O, CO2 and H2 increase steadily inside the reaction zone; they stay appreciable in the combustion products. However, unstable products, such as HO2, H2O2 and CH3, first increase in the preheating region of the reaction zone; they are then consumed in the remaining part of the reaction zone. The numerical results show that most of the formed NOx is composed of nitric oxide. The velocity and temperature profiles were accurately predicted using a grid of 80×80 while the nitrogen oxides were computed accurately utilizing a finer grid of 160×160.


2020 ◽  
Vol 7 (1) ◽  
pp. 48-55 ◽  
Author(s):  
Bolat Duissenbekov ◽  
Abduhalyk Tokmuratov ◽  
Nurlan Zhangabay ◽  
Zhenis Orazbayev ◽  
Baisbay Yerimbetov ◽  
...  

AbstractThe study solves a system of finite difference equations for flexible shallow concrete shells while taking into account the nonlinear deformations. All stiffness properties of the shell are taken as variables, i.e., stiffness surface and through-thickness stiffness. Differential equations under consideration were evaluated in the form of algebraic equations with the finite element method. For a reinforced shell, a system of 98 equations on a 8×8 grid was established, which was next solved with the approximation method from the nonlinear plasticity theory. A test case involved computing a 1×1 shallow shell taking into account the nonlinear properties of concrete. With nonlinear equations for the concrete creep taken as constitutive, equations for the quasi-static shell motion under constant load were derived. The resultant equations were written in a differential form and the problem of solving these differential equations was then reduced to the solving of the Cauchy problem. The numerical solution to this problem allows describing the stress-strain state of the shell at each point of the shell grid within a specified time interval.


Sign in / Sign up

Export Citation Format

Share Document