Discontinuous Galerkin method for multicomponent chemically reacting flows and combustion

This paper presents the development of a discontinuous Galerkin (DG) method for application to chemically reacting flows in subsonic and supersonic regimes under the consideration of variable thermo-viscous-diffusive transport properties, detailed and stiff reaction chemistry, and shock capturing. A hybrid-flux formulation is developed for treatment of the convective fluxes, combining a conservative Riemann-solver and an extended double-flux scheme. A computationally efficient splitting scheme is proposed, in which advection and diffusion operators are solved in the weak form, and the chemically stiff substep is advanced in the strong form using a time-implicit scheme. The discretization of the viscous-diffusive transport terms follows the second form of Bassi and Rebay, and the WENO-based limiter due to Zhong and Shu is extended to multicomponent systems. Boundary conditions are developed for subsonic and supersonic flow conditions, and the algorithm is coupled to thermochemical libraries to account for detailed reaction chemistry and complex transport. The resulting DG method is applied to a series of test cases of increasing physico-chemical complexity. Beginning with one- and two-dimensional multispecies advection and shock–fluid interaction problems, computational efficiency, convergence, and conservation properties are demonstrated. This study is followed by considering a series of detonation and supersonic combustion problems to investigate the convergence-rate and the shock-capturing capability in the presence of one- and multistep reaction chemistry. The DG algorithm is then applied to diffusion-controlled deflagration problems. By examining convergence properties for polynomial order and spatial resolution, and comparing these with second-order finite-volume solutions, it is shown that optimal convergence is achieved and that polynomial refinement provides advantages in better resolving the localized flame structure and complex flow-field features associated with multidimensional and hydrodynamic/thermo-diffusive instabilities in deflagration and detonation systems. Comparisons with standard third- and fifth-order WENO schemes are presented to illustrate the benefit of the DG scheme for application to detonation and multispecies flow/shock-interaction problems.

[1]  S. Penner.  Physics of shock waves and high-temperature hydrodynamic phenomena - Ya.B. Zeldovich and Yu.P. Raizer (translated from the Russian and then edited by Wallace D. Hayes and Ronald F. Probstein); Dover Publications, New York, 2002, 944 pp., $34. , 2003 .

[2]  Hiroshi Terashima,et al.  A front-tracking/ghost-fluid method for fluid interfaces in compressible flows , 2009, J. Comput. Phys..

[3]  Chi-Wang Shu,et al.  High Order Strong Stability Preserving Time Discretizations , 2009, J. Sci. Comput..

[4]  Charles K. Westbrook,et al.  Chemical kinetics of hydrocarbon oxidation in gaseous detonations , 1982 .

[5]  Cheng Wang,et al.  Robust high order discontinuous Galerkin schemes for two-dimensional gaseous detonations , 2012, J. Comput. Phys..

[6]  S. Rebay,et al.  GMRES Discontinuous Galerkin Solution of the Compressible Navier-Stokes Equations , 2000 .

[7]  Robert W. Bilger.  Turbulent diffusion flames , 1989 .

[8]  Chi-Wang Shu,et al.  A simple weighted essentially nonoscillatory limiter for Runge-Kutta discontinuous Galerkin methods , 2013, J. Comput. Phys..

[9]  E. Toro,et al.  Restoration of the contact surface in the HLL-Riemann solver , 1994 .

[10]  James J. Quirk,et al.  On the dynamics of a shock–bubble interaction , 1994, Journal of Fluid Mechanics.

[11]  Ralf Deiterding,et al.  Parallel adaptive simulation of multi-dimensional detonation structures , 2003 .

[12]  Elaine S. Oran,et al.  A Numerical Study of a Two-Dimensional H2-O2-Ar Detonation Using a Detailed Chemical Reaction Model , 1998 .

[13]  Forman A. Williams,et al.  Fundamental Aspects of Combustion , 1963, Nature.

[14]  Amable Liñán,et al.  The asymptotic structure of counterflow diffusion flames for large activation energies , 1974 .

[15]  Gary J. Sharpe,et al.  Transverse waves in numerical simulations of cellular detonations , 2001, Journal of Fluid Mechanics.

[16]  Sanford Gordon,et al.  NASA Glenn Coefficients for Calculating Thermodynamic Properties of Individual Species , 2002 .

[17]  Eric Johnsen,et al.  Implementation of WENO schemes in compressible multicomponent flow problems , 2005, J. Comput. Phys..

[18]  R. Borghi.  Turbulent combustion modelling , 1988 .

[19]  Rémi Abgrall,et al.  Computations of compressible multifluids , 2001 .

[20]  R. Abgrall.  How to Prevent Pressure Oscillations in Multicomponent Flow Calculations , 1996 .

[21]  G.,et al.  Combustion theory and modeling , 2008 .

[22]  K. Mazaheri,et al.  High-resolution numerical simulation of the structure of 2-D gaseous detonations , 2011 .

[23]  F. Williams,et al.  Theory of Combustion in Laminar Flows , 1971 .

[24]  D. Goodwin,et al.  Cantera: An Object-oriented Software Toolkit for Chemical Kinetics, Thermodynamics, and Transport Processes. Version 2.2.0 , 2015 .

[25]  John H. S. Lee.  The Detonation Phenomenon , 2008 .

[26]  Keh-Ming Shyue,et al.  An Efficient Shock-Capturing Algorithm for Compressible Multicomponent Problems , 1998 .

[27]  Eric Johnsen,et al.  Preventing numerical errors generated by interface-capturing schemes in compressible multi-material flows , 2012, J. Comput. Phys..

[28]  Xiaolong Gou,et al.  A dynamic multi-timescale method for combustion modeling with detailed and reduced chemical kinetic mechanisms , 2010 .

[29]  Robert K. Cheng,et al.  Numerical simulation of Lewis number effects on lean premixed turbulent flames , 2007 .

[30]  Elaine S. Oran,et al.  Two-dimensional reactive flow dynamics in cellular detonation waves , 1999 .

[31]  M. Y. Hussaini,et al.  An Analysis of the Discontinuous Galerkin Method for Wave Propagation Problems , 1999 .

[32]  Joseph E. Shepherd,et al.  Direct observations of reaction zone structure in propagating detonations , 2003 .

[33]  W. H. Reed,et al.  Triangular mesh methods for the neutron transport equation , 1973 .

[34]  Chi-Wang Shu,et al.  Runge–Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems , 2001, J. Sci. Comput..

[35]  Elaine S. Oran,et al.  Origins of the deflagration-to-detonation transition in gas-phase combustion , 2007 .

[36]  Spencer J. Sherwin,et al.  Dispersion Analysis of the Continuous and Discontinuous Galerkin Formulations , 2000 .

[37]  J. Haas,et al.  Interaction of weak shock waves with cylindrical and spherical gas inhomogeneities , 1987, Journal of Fluid Mechanics.

[38]  Ronald Fedkiw,et al.  Numerical resolution of pulsating detonation waves , 2000 .

[39]  G. Sivashinsky,et al.  Instabilities, Pattern Formation, and Turbulence in Flames , 1983 .

[40]  P. Mulet,et al.  A flux-split algorithm applied to conservative models for multicomponent compressible flows , 2003 .

[41]  R. Bonazza,et al.  Shock-Bubble Interactions , 2011 .

[42]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[43]  J. Austin,et al.  The role of instability in gaseous detonation , 2003 .

[44]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[45]  Kenneth K. Kuo,et al.  A low-dissipation and time-accurate method for compressible multi-component flow with variable specific heat ratios , 2011, J. Comput. Phys..

[46]  Xiangxiong Zhang,et al.  Positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations with source terms , 2011, J. Comput. Phys..

[47]  Rémi Abgrall,et al.  An adaptive shock-capturing algorithm for solving unsteady reactive flows , 2003 .

[48]  Chung King Law,et al.  The hydrodynamic structure of unstable cellular detonations , 2007, Journal of Fluid Mechanics.

[49]  J. Ryan,et al.  A Runge-Kutta discontinuous Galerkin approach to solve reactive flows: The hyperbolic operator , 2011, J. Comput. Phys..

[50]  F. Williams,et al.  Analytical studies of the dynamics of gaseous detonations , 2012, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[51]  George Em Karniadakis,et al.  The Development of Discontinuous Galerkin Methods , 2000 .

[52]  Matthias Ihme,et al.  Development of Discontinuous Galerkin Method for Detonation and Supersonic Combustion , 2013 .

[53]  Ralf Deiterding,et al.  An adaptive high-order hybrid scheme for compressive, viscous flows with detailed chemistry , 2011, J. Comput. Phys..