Additive Semi-Implicit Runge-Kutta Methods for Computing High-Speed Nonequilibrium Reactive Flows

This paper is concerned with time-stepping numerical methods for computing stiff semi-discrete systems of ordinary differential equations for transient hypersonic flows with thermo-chemical nonequilibrium. The stiffness of the equations is mainly caused by the viscous flux terms across the boundary layers and by the source terms modeling finite-rate thermo-chemical processes. Implicit methods are needed to treat the stiff terms while more efficient explicit methods can still be used for the nonstiff terms in the equations. This paper studies three different semi-implicit Runge?Kutta methods for additively split differential equations in the form ofu? =f(u) +g(u), wherefis treated by explicit Runge?Kutta methods andgis simultaneously treated by three implicit Runge?Kutta methods: a diagonally implicit Runge?Kutta method and two linearized implicit Runge?Kutta methods. The coefficients of up to third-order accurate additive semi-implicit Runge?Kutta methods have been derived such that the methods are both high-order accurate and strongly A-stable for the implicit terms. The results of two numerical tests on the stability and accuracy properties of these methods are also presented in the paper.

[1]  Leonhard Kleiser,et al.  Numerical simulation of transition in wall-bounded shear flows , 1991 .

[2]  N. N. Yanenko,et al.  The Method of Fractional Steps , 1971 .

[3]  E. Hairer,et al.  On the stability of semi-implicit methods for ordinary differential equations , 1982 .

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

[5]  J. Butcher.  Implicit Runge-Kutta processes , 1964 .

[6]  T. A. Zang,et al.  Spectral methods for modeling supersonic chemically reacting flowfields , 1985 .

[7]  T. A. Zang,et al.  Spectral methods for fluid dynamics , 1987 .

[8]  X. Zhong,et al.  Simulation of transient hypersonic flow using the ENO schemes , 1995 .

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

[10]  T. A. Zang,et al.  Numerical experiments on subcritical transition mechanisms , 1985 .

[11]  I. Gladwell,et al.  Computational techniques for ordinary differential equations , 1982 .

[12]  Parviz Moin,et al.  On the numerical solution of time-dependent viscous incompressible fluid flows involving solid boundaries , 1980 .

[13]  P. Rentrop,et al.  Generalized Runge-Kutta methods of order four with stepsize control for stiff ordinary differential equations , 1979 .

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

[15]  H. H. Rosenbrock,et al.  Some general implicit processes for the numerical solution of differential equations , 1963, Comput. J..

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

[17]  G. J. Cooper,et al.  Additive Runge-Kutta methods for stiff ordinary differential equations , 1983 .

[18]  José Carlos Goulart de Siqueira,et al.  Differential Equations , 1919, Nature.

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

[20]  R. May,et al.  The Numerical Solution of Ordinary Differential Equations: Initial Value Problems , 1984 .