A study of numerical methods for hyperbolic conservation laws with stiff source terms

The proper modeling of nonequilibrium gas dynamics is required in certain regimes of hypersonic flow. For inviscid flow this gives a system of conservation laws coupled with source terms representing the chemistry. Often a wide range of time scales is present in the problem, leading to numerical difficulties as in stiff systems of ordinary differential equations. Stability can be achieved by using implicit methods, but other numerical difficulties are observed. The behavior of typical numerical methods on a simple advection equation with a parameter-dependent source term was studied. Two approaches to incorporate the source term were utilized: MacCormack type predictor-corrector methods with flux limiters, and splitting methods in which the fluid dynamics and chemistry are handled in separate steps. Various comparisons over a wide range of parameter values were made. In the stiff case where the solution contains discontinuities, incorrect numerical propagation speeds are observed with all of the methods considered. This phenomenon is studied and explained.

[1]  H. C. Yee.  Upwind and Symmetric Shock-Capturing Schemes , 1987 .

[2]  H. C. Yee,et al.  Numerical simulation of shock wave diffraction by TVD schemes , 1987 .

[3]  H. C. Yee,et al.  Semi-implicit and fully implicit shock-capturing methods for hyperbolic conservation laws with stiff source terms , 1987 .

[4]  G. C. Carofano,et al.  Blast Computation Using Harten's Total Variation Diminishing Scheme , 1984 .

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

[6]  Jong-Hun Lee,et al.  Basic governing equations for the flight regimes of aeroassisted orbital transfer vehicles , 1984 .

[7]  T. Teichmann,et al.  Introduction to physical gas dynamics , 1965 .

[8]  Paul Kutler,et al.  Second- and Third-Order Noncentered Difference Schemes for Nonlinear Hyperbolic Equations , 1973 .

[9]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[10]  J. M. Watt.  Numerical Initial Value Problems in Ordinary Differential Equations , 1972 .

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

[12]  P. Lax,et al.  Systems of conservation laws , 1960 .

[13]  C. W. Gear,et al.  Numerical initial value problem~ in ordinary differential eqttations , 1971 .

[14]  J. Boris,et al.  A numerical technique for solving stiff ordinary differential equations associated with the chemical kinetics of reactive-flow problems , 1977 .

[15]  R. Maccormack.  The Effect of Viscosity in Hypervelocity Impact Cratering , 1969 .

[16]  Elaine S. Oran,et al.  Numerical Simulation of Reactive Flow , 1987 .

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

[18]  J. Lambert.  Computational Methods in Ordinary Differential Equations , 1973 .