Bunsen flame simulation by finite elements on adaptively refined, unstructured triangulations

We investigate laminar Bunsen flames with detailed chemistry and multicomponent transport. The governing equations are discretized by a finite element method on a sequence of adaptively refined, unstructured triangulations. The finite element method is an extension to chemically reacting flows of the streamline diffusion method, including least-squares stabilization of the pressure gradient and the low-Mach continuity equation as well as a shock capturing term designed to control species mass fraction undershoots near flame fronts. Unstructured meshes are adaptively refined based on a posteriori estimates of a user specified functional of the numerical error. These estimates are derived from the dual weighted residual method in the form of element-wise residuals weighted by coefficients depending on the solution of a linearized dual problem that accounts for convective error propagation and multicomponent chemistry couplings. Numerical results are presented to illustrate the efficiency of the proposed methodology and to study the impact of inflow velocity profiles on the structure of several hydrogen–air Bunsen flames.

[1]  V. Giovangigli Multicomponent flow modeling , 1999 .

[2]  M. Smooke,et al.  Adaptive continuation algorithms with application to combustion problems , 1989 .

[3]  de Lph Philip Goey,et al.  A Numerical Study of a Premixed Flame on a Slit Burner , 1995 .

[4]  Alexandre Ern,et al.  Impact of Detailed Multicomponent Transport on Planar and Counterflow Hydrogen/Air and Methane/Air Flames , 1999 .

[5]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[6]  L. Franca,et al.  Stabilized finite element methods. II: The incompressible Navier-Stokes equations , 1992 .

[7]  Rolf Rannacher,et al.  A Feed-Back Approach to Error Control in Finite Element Methods: Basic Analysis and Examples , 1996 .

[8]  P. Hansbo,et al.  Edge stabilization for Galerkin approximations of convection?diffusion?reaction problems , 2004 .

[9]  P. Hansbo,et al.  Mathematical Modelling and Numerical Analysis Edge Stabilization for the Generalized Stokes Problem: a Continuous Interior Penalty Method , 2022 .

[10]  Alexandre Ern,et al.  An adaptive finite element method with crosswind diffusion for low Mach, steady, laminar combustion , 2003 .

[11]  O. C. Zienkiewicz,et al.  A simple error estimator and adaptive procedure for practical engineerng analysis , 1987 .

[12]  Rolf Rannacher,et al.  An optimal control approach to a posteriori error estimation in finite element methods , 2001, Acta Numerica.

[13]  M S Day,et al.  Numerical simulation of laminar reacting flows with complex chemistry , 2000 .

[14]  Mitchell D. Smooke,et al.  Solution of burner-stabilized premixed laminar flames by boundary value methods , 1982 .

[15]  B. Bennett,et al.  Local rectangular refinement with application to axisymmetric laminar flames , 1998 .

[16]  Forman A. Williams,et al.  Effects of molecular diffusion and of thermal expansion on the structure and dynamics of premixed flames in turbulent flows of large scale and low intensity , 1982, Journal of Fluid Mechanics.

[17]  Robert J. Kee,et al.  Toward a comprehensive chemical kinetic mechanism for the oxidation of acetylene: Comparison of model predictions with results from flame and shock tube experiments , 1982 .

[18]  Lutz Tobiska,et al.  Analysis of a streamline diffusion finite element method for the Stokes and Navier-Stokes equations , 1996 .

[19]  Alexandre Ern,et al.  Stabilized Galerkin approximation of convection-diffusion-reaction equations: discrete maximum principle and convergence , 2005, Math. Comput..

[20]  S. Rebay Efficient Unstructured Mesh Generation by Means of Delaunay Triangulation and Bowyer-Watson Algorithm , 1993 .

[21]  de Hc Rick Lange,et al.  Numerical flow modelling in a locally refined grid , 1994 .

[22]  Claes Johnson,et al.  Finite element methods for linear hyperbolic problems , 1984 .

[23]  Claes Johnson Numerical solution of partial differential equations by the finite element method , 1988 .

[24]  A. H. Schatz,et al.  Crosswind Smear and Pointwise Errors in Streamline Diffusion Finite Element Methods , 1987 .

[25]  Rolf Rannacher,et al.  Numerical simulation of laminar flames at low Mach number by adaptive finite elements , 1999 .

[26]  de Hc Rick Lange,et al.  MODELING OF CONFINED AND UNCONFINED LAMINAR PREMIXED FLAMES ON SLIT AND TUBE BURNERS , 1995 .

[27]  Alexandre Ern,et al.  Multicomponent transport algorithms , 1994 .

[28]  Pedro J. Coelho,et al.  Calculation of a Confined Axisymmetric Laminar Diffusion Flame Using a Local Grid Refinement Technique , 1993 .

[29]  N. Darabiha,et al.  Vector Computers and Complex Chemistry Combustion , 1988 .

[30]  Alexandre Ern,et al.  Multicomponent transport impact on turbulent premixed H2/O2 flames , 2002 .

[31]  Alexandre Ern,et al.  Thermal diffusion effects in hydrogen-air and methane-air flames , 1998 .