Homogeneous gas phase reaction rate data pertaining to catalytically stabilized combustion systems at different flow velocities
Description
The homogeneous gas phase reaction rate data pertaining to catalytically stabilized combustion systems at different flow velocities are obtained by performing numerical simulations and using fluid mechanics. A segregated solution solver with an under-relaxation method is used to solve the conservation equations. The segregated solver first solves the momentum equations, then solves the continuity equation, and updates the pressure and mass flow rate. The energy and species equations are subsequently solved and convergence is checked. The latter is monitored through both the values of the residuals of the conservation equations and the difference between subsequent iterations of the solution. The boundary conditions are defined as follows. At the inlet, a fixed, flat velocity profile is used. This boundary condition fixes the convective component of the flux of species and energy, but the diffusive component depends on the gradient of the computed temperature or species fields. Symmetry boundary conditions are applied at the centerline between the two plates. At the exit, a fixed pressure is specified and far-field conditions are imposed for the rest of the variables. At the interface between the wall and the fluid, no-slip boundary is employed. The heat flux at the fluid-wall interface is computed using Fourier's law and continuity in temperature and heat flux links the fluid and solid phases. All internal heat transfer between the fluid and the wall is calculated by accounting explicitly for the convective and conductive heat transport in the model within the fluid and within the wall. The wall thermal conductivity and exterior convective heat loss coefficient are taken as independent parameters to understand how important thermal management is. The fluid density is calculated using the ideal gas law. The fluid viscosity, specific heat, and thermal conductivity are calculated from a mass fraction weighted average of species properties, and the specific heat of chemical species is calculated using a piecewise polynomial fit of temperature. Contributor: Junjie Chen, E-mail: koncjj@gmail.com, Department of Energy and Power Engineering, School of Mechanical and Power Engineering, Henan Polytechnic University, 2000 Century Avenue, Jiaozuo, Henan, 454000, P.R. China
Files
Steps to reproduce
The momentum, species, and energy equations are discretized using a first-order upwind approximation. The pressure was discretized using a standard method. The pressure-velocity coupling is discretized using the SIMPLE method. When parallel processing is used, the message passing interface is used to transmit information between nodes. In order to achieve convergence as well as compute extinction points, natural parameter continuation is implemented. Convergence is determined from the residuals of the conservation equations as well as the difference between subsequent iterations of the solution.