skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: High Order Finite Difference Methods with Subcell Resolution for Stiff Multispecies Discontinuity Capturing

Journal Article · · Communications in Computational Physics
 [1];  [2];  [3];  [4];  [5]
  1. Florida International Univ., Miami, FL (United States)
  2. Brown Univ., Providence, RI (United States)
  3. NASA Ames Research Center (ARC), Moffett Field, Mountain View, CA (United States)
  4. Center for Turbulence Research, Stanford, CA (United States)
  5. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
+ Show Author Affiliations

Here, we extend the high order finite-difference method with subcell resolution (SR) in [34] for two-species stiff one-reaction models to multispecies and multireaction inviscid chemical reactive flows, which are significantly more difficult because of the multiple scales generated by different reactions. For reaction problems, when the reaction time scale is very small, the reaction zone scale is also small and the governing equations become very stiff. Wrong propagation speed of discontinuity may occur due to the underresolved numerical solution in both space and time. The present SR method for reactive Euler system is a fractional step method. In the convection step, any high order shock-capturing method can be used. In the reaction step, an ODE solver is applied but with certain computed flow variables in the shock region modified by the Harten subcell resolution idea. Several numerical examples of multispecies and multireaction reactive flows are performed in both one and two dimensions. Studies demonstrate that the SR method can capture the correct propagation speed of discontinuities in very coarse meshes.

Research Organization:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Sponsoring Organization:
USDOE National Nuclear Security Administration (NNSA); National Science Foundation (NSF); National Aeronautics and Space Administration (NASA)
Grant/Contract Number:
AC52-07NA27344; NNX12AJ62A; DMS-1112700
OSTI ID:
1811773
Report Number(s):
LLNL-JRNL-823625; 1035349
Journal Information:
Communications in Computational Physics, Vol. 17, Issue 2; ISSN 1815-2406
Publisher:
Global Science PressCopyright Statement
Country of Publication:
United States
Language:
English

References (29)

A Modified Fractional Step Method for the Accurate Approximation of Detonation Waves journal January 2000
Detonation capturing for stiff combustion chemistry journal September 1998
Random choice methods with applications to reacting gas flow journal November 1977
Efficient implementation of essentially non-oscillatory shock-capturing schemes, II journal July 1989
A modified higher order Godunov’s scheme for stiff source conservative hydrodynamics journal June 2007
The generalized Riemann problem for reactive flows journal March 1989
Error estimators for the position of discontinuities in hyperbolic conservation laws with source terms which are solved using operator splitting journal July 1999
A study of numerical methods for hyperbolic conservation laws with stiff source terms journal January 1990
Positivity-preserving method for high-order conservative schemes solving compressible Euler equations journal June 2013
ADER Schemes for Nonlinear Systems of Stiff Advection–Diffusion–Reaction Equations journal October 2010
Heterogeneous Multiscale Methods for Interface Tracking of Combustion Fronts journal January 2006
The Random Projection Method for Hyperbolic Conservation Laws with Stiff Reaction Terms journal September 2000
Numerical Study of the Mechanisms for Initiation of Reacting Shock Waves journal September 1990
One-Dimensional Front Tracking Based on High Resolution Wave Propagation Methods journal March 1995
The Random Projection Method for Stiff Detonation Capturing journal January 2001
Numerical Methods for Hyperbolic Conservation Laws With Stiff Relaxation I. Spurious Solutions journal October 1993
Spurious behavior of shock-capturing methods by the fractional step approach: Problems containing stiff source terms and discontinuities journal May 2013
High order finite difference methods with subcell resolution for advection equations with stiff source terms journal January 2012
Numerical Wave Propagation in an Advection Equation with a Nonlinear Source Term journal October 1992
Multiresolution Schemes for the Reactive Euler Equations journal September 1999
Numerical solution of under-resolved detonations journal February 2008
On the Construction and Comparison of Difference Schemes journal September 1968
Improved Shock-Capturing Methods for Multicomponent and Reacting Flows journal October 1996
Theoretical and Numerical Structure for Reacting Shock Waves journal October 1986
Theoretical and Numerical Structure for Unstable One-Dimensional Detonations journal April 1991
Efficient Implementation of Weighted ENO Schemes journal June 1996
ENO schemes with subcell resolution journal July 1989
Random choice solution of hyperbolic systems journal December 1976
The Random Projection Method for Stiff Multispecies Detonation Capturing journal May 2002

Similar Records

High order finite difference methods with subcell resolution for advection equations with stiff source terms
Journal Article · Sun Jan 01 00:00:00 EST 2012 · Journal of Computational Physics · OSTI ID:1811773
+1 more
On the application of subcell resolution to conservation laws with stiff source terms
Technical Report · Wed Nov 01 00:00:00 EST 1989 · OSTI ID:1811773
Subcell resolution in simplex stochastic collocation for spatial discontinuities
Journal Article · Tue Oct 15 00:00:00 EDT 2013 · Journal of Computational Physics · OSTI ID:1811773

Related Subjects

97 MATHEMATICS AND COMPUTING
mathematics and computing
65M06
76V05
80A32
stiff reaction term
shock capturing
detonation
WENO
ENO subcell resolution
multispecies
multireactions