On spurious numerics in solving reactive equations

The objective of this study is to gain a deeper understanding of the behavior of high order shock-capturing schemes for problems with stiff source terms and discontinuities and on corresponding numerical prediction strategies. The studies by Yee et al. (2012) and Wang et al. (2012) focus only on solving the reactive system by the fractional step method using the Strang splitting (Strang 1968). It is a common practice by developers in computational physics and engineering simulations to include a cut off safeguard if densities are outside the permissible range. Here we compare the spurious behavior of the same schemes by solving the fully coupled reactive system without the Strang splitting vs. using the Strang splitting. Comparison between the two procedures and the effects of a cut off safeguard is the focus the present study. The comparison of the performance of these schemes is largely based on the degree to which each method captures the correct location of the reaction front for coarse grids. Here "coarse grids" means standard mesh density requirement for accurate simulation of typical non-reacting flows of similar problem setup. It is remarked that, in order to resolve the sharp reaction front, local refinement beyond standard mesh density is still needed.

[1]  P. Olsson.  Summation by parts, projections, and stability. II , 1995 .

[2]  N. A. Adams,et al.  Positivity-preserving flux limiters for high-order conservative schemes , 2012 .

[3]  Chi-Wang Shu,et al.  Construction of low dissipative high-order well-balanced filter schemes for non-equilibrium flows , 2009, J. Comput. Phys..

[4]  Olsson,et al.  SUMMATION BY PARTS, PROJECTIONS, AND STABILITY. I , 2010 .

[5]  Björn Sjögreen,et al.  Development of low dissipative high order filter schemes for multiscale Navier-Stokes/MHD systems , 2006, J. Comput. Phys..

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

[7]  G. Strang.  On the Construction and Comparison of Difference Schemes , 1968 .

[8]  Neil D. Sandham,et al.  Low-Dissipative High-Order Shock-Capturing Methods Using Characteristic-Based Filters , 1999 .

[9]  Björn Sjögreen,et al.  High Order Filter Methods for Wide Range of Compressible Flow Speeds , 2011 .

[10]  Xiangxiong Zhang,et al.  Positivity-preserving high order finite difference WENO schemes for compressible Euler equations , 2012, J. Comput. Phys..

[11]  Luigi Vigevano,et al.  Numerical solution of under-resolved detonations , 2008, J. Comput. Phys..

[12]  Vincent Guinot,et al.  High-Order Fluxes for Conservative Skew-Symmetric-like Schemes in Structured Meshes , 2000 .

[13]  Björn Sjögreen,et al.  Multiresolution Wavelet Based Adaptive Numerical Dissipation Control for Shock-Turbulence Computations , 2001 .

[14]  Nikolaus A. Adams,et al.  Positivity-preserving method for high-order conservative schemes solving compressible Euler equations , 2013, J. Comput. Phys..

[15]  A. Harten.  ENO schemes with subcell resolution , 1989 .

[16]  H. C. Yee,et al.  Entropy Splitting and Numerical Dissipation , 2000 .

[17]  Björn Sjögreen,et al.  ON TENTH-ORDER CENTRAL SPATIAL SCHEMES , 2007, Proceeding of Fifth International Symposium on Turbulence and Shear Flow Phenomena.

[18]  Chi-Wang Shu,et al.  High order finite difference methods with subcell resolution for advection equations with stiff source terms , 2012, J. Comput. Phys..

[19]  Randall J. LeVeque,et al.  A Modified Fractional Step Method for the Accurate Approximation of Detonation Waves , 2000, SIAM J. Sci. Comput..

[20]  H. C. Yee,et al.  Spurious Behavior of Shock-Capturing Methods: Problems Containing Stiff Source Terms and Discontinuities , 2013 .

[21]  Björn Sjögreen,et al.  Adaptive filtering and limiting in compact high order methods for multiscale gas dynamics and MHD systems , 2008 .

[22]  H. C. Yee,et al.  Extension of Efficient Low Dissipative High Order Schemes for 3-D Curvilinear Moving Grids , 2000 .