Spurious behavior of shock-capturing methods by the fractional step approach: Problems containing stiff source terms and discontinuities

The goal of this paper is to relate numerical dissipations that are inherited in high order shock-capturing schemes with the onset of wrong propagation speed of discontinuities. For pointwise evaluation of the source term, previous studies indicated that the phenomenon of wrong propagation speed of discontinuities is connected with the smearing of the discontinuity caused by the discretization of the advection term. The present study focuses only on solving the reactive system by the fractional step method using the Strang splitting. Studies shows that the degree of wrong propagation speed of discontinuities is highly dependent on the accuracy of the numerical method. The manner in which the smearing of discontinuities is contained by the numerical method and the overall amount of numerical dissipation being employed play major roles. Depending on the numerical method, time step and grid spacing, the numerical simulation may lead to (a) the correct solution (within the truncation error of the scheme), (b) a divergent solution, (c) a wrong propagation speed of discontinuities solution or (d) other spurious solutions that are solutions of the discretized counterparts but are not solutions of the governing equations. The findings might shed some light on the reported difficulties in numerical combustion and problems with stiff nonlinear (homogeneous) source terms and discontinuities in general.

[1]  H. C. Yee.  DESIGNING ADAPTIVE LOW-DISSIPATIVE HIGH ORDER SCHEMES FOR LONG-TIME INTEGRATIONS , .

[2]  Wei Wang,et al.  On spurious numerics in solving reactive equations , 2013 .

[3]  D. Nguyen.  A Fully Conservative Ghost Fluid Method & Stiff Detonation Waves , 2002 .

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

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

[6]  Andrew J. Majda,et al.  Theoretical and numerical structure for unstable one-dimensional detonations , 1991 .

[7]  Björn Sjögreen,et al.  Multiresolution Wavelet Based Adaptive Numerical Dissipation Control for High Order Methods , 2004, J. Sci. Comput..

[8]  Joseph Oliger,et al.  Energy and Maximum Norm Es-timates for Nonlinear Conservation Laws , 1994 .

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

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

[11]  H. C. Yee,et al.  Multiresolution Wavelet Based Adaptive Numerical Dissipation Control for Shock-Turbulence Computations , 2001 .

[12]  Randall J. LeVeque,et al.  One-Dimensional Front Tracking Based on High Resolution Wave Propagation Methods , 1995, SIAM J. Sci. Comput..

[13]  Andrea Lani,et al.  Variable High-Order Multiblock Overlapping Grid Methods for Mixed Steady and Unsteady Multiscale Viscous Flows, Part II: Hypersonic Nonequilibrium Flows , 2013 .

[14]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

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

[16]  H. C. Yee.  BUILDING BLOCKS FOR RELIABLE COMPLEX NONLINEAR NUMERICAL SIMULATIONS , .

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

[18]  H. C. Yee,et al.  Numerical wave propagation in an advection equation with a nonlinear source term , 1992 .

[19]  E. F. Kaasschieter,et al.  Detonation capturing for stiff combustion chemistry , 1998 .

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

[21]  H. C. Yee,et al.  Dynamical Approach Study of Spurious Steady-State Numerical Solutions of Nonlinear Differential Equations, Part III. The Effects of Nonlinear Source Terms in Reaction-Convection Equations , 1996 .

[22]  Weizhu Bao,et al.  The Random Projection Method for Hyperbolic Conservation Laws with Stiff Reaction Terms , 2000 .

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

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

[25]  H. Yeea,et al.  Numerical Dissipation and Wrong Propagation Speed of Discontinuities For Stiff Source Terms , 2011 .

[26]  H. C. Yee,et al.  DYNAMICAL APPROACH STUDY OF SPURIOUS STEADY-STATE NUMERICAL SOLUTIONS OF NONLINEAR DIFFERENTIAL EQUATIONS II. GLOBAL ASYMPTOTIC BEHAVIOR OF TIME DISCRETIZATIONS ∗ , 1995 .

[27]  Matania Ben-Artzi,et al.  The generalized Riemann problem for reactive flows , 1989 .

[28]  H. C. Yee,et al.  Dynamical Approach Study of Spurious Steady-State Numerical Solutions of Nonlinear Differential Equations Part IV. Stability vs. Methods of Discretizing Nonlinear Source Terms in Reaction-Convection Equations , 1996 .

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

[30]  Björn Sjögreen,et al.  Comparative Study on High-Order Positivity-preserving WENO Schemes , 2014 .

[31]  Elaine S. Oran,et al.  Determination of detonation cell size and the role of transverse waves in two-dimensional detonations☆ , 1985 .

[32]  Chi-Wang Shu.  Total-variation-diminishing time discretizations , 1988 .

[33]  Hassan Hassan,et al.  Study of Radiation in Electric Arc Shock Tubes , 2010 .

[34]  Prabhu Ramachandran,et al.  Approximate Riemann solvers for the Godunov SPH (GSPH) , 2014, J. Comput. Phys..

[35]  Patrick Jenny,et al.  Correction of Conservative Euler Solvers for Gas Mixtures , 1997 .

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

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

[38]  Albert Edward Honein,et al.  Numerical aspects of compressible turbulence simulations , 2005 .

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

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

[41]  Rolf Jeltsch,et al.  Error estimators for the position of discontinuities in hyperbolic conservation laws with source terms which are solved using operator splitting , 1999 .

[42]  A. R. Humphries,et al.  Dynamical Systems And Numerical Analysis , 1996 .

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

[44]  Barna L. Bihari,et al.  Multiresolution Schemes for the Reactive Euler Equations , 1999 .

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

[46]  P. Colella,et al.  Theoretical and numerical structure for reacting shock waves , 1986 .

[47]  H. C. Yee,et al.  A class of high resolution explicit and implicit shock-capturing methods , 1989 .

[48]  Richard B. Pember,et al.  Numerical Methods for Hyperbolic Conservation Laws With Stiff Relaxation I. Spurious Solutions , 1993, SIAM J. Appl. Math..

[49]  Randall J. LeVeque,et al.  A study of numerical methods for hyperbolic conservation laws with stiff source terms , 1990 .

[50]  H. C. Yee,et al.  High-Resolution Shock-Capturing Schemes for Inviscid and Viscous Hypersonic Flows , 1990 .

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

[52]  Björn Sjögreen,et al.  On Skew-Symmetric Splitting and Entropy Conservation Schemes for the Euler Equations , 2010 .

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

[54]  H. C. Yee,et al.  Dynamical approach study of spurious steady-state numerical solutions of nonlinear differential equations. I. The dynamics of time discretization and its implications for algorithm development in computational fluid dynamics☆ , 1991 .

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

[56]  Miguel R. Visbal,et al.  On spurious behavior of CFD simulations , 1999 .

[57]  H. C. Yee,et al.  Dynamics of Numerics & Spurious Behaviors in CFD Computations. Revised , 1997 .

[58]  Björn Sjögreen,et al.  Comparative Study of Three High Order Schemes for LES of Temporally Evolving Mixing Layers , 2012 .

[59]  V. Ton,et al.  Improved Shock-Capturing Methods for Multicomponent and Reacting Flows , 1996 .

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

[61]  Olivier Chazot,et al.  Electronic Excitation of Atoms and Molecules for the FIRE II Flight Experiment , 2011 .