On stabilized integration for time-dependent PDEs
Journal of Computational Physics , Volume 224 - Issue 1 p. 3- 16
An integration method is discussed which has been designed to treat parabolic and hyperbolic terms explicitly and stiff reaction terms implicitly. The method is a special two-step form of the one-step IMEX (Implicit-Explicit) RKC (Runge-Kutta-Chebyshev) method. The special two-step form is introduced with the aim of getting a non-zero imaginary stability boundary which is zero for the one-step method. Having a non-zero imaginary stability boundary allows, for example, the integration of pure advection equations space-discretized with centered schemes, the integration of damped or viscous wave equations, the integration of coupled sound and heat flow equations, etc. For our class of methods it also simplifies the choice of temporal step sizes satisfying the von Neumann stability criterion, by embedding a thin long rectangle inside the stability region. Embedding rectangles or other tractable domains with this purpose is an idea of Wesseling.
|Numerical Integration, Stabilized Explicit Integration, Runge-Kutta-Chebyshev Methods, Reactive Flow Problems, Damped Wave Equations, Coupled Sound and Heat Flow|
|Interpolation (acm G.1.1), Ordinary Differential Equations (acm G.1.7), Partial Differential Equations (acm G.1.8)|
|Stability and convergence of numerical methods (msc 65M12), Method of lines (msc 65M20), Initial value problems (msc 65L05)|
|Journal of Computational Physics|
|Development of a State-of-the Art Surface-Capturing Method for Two-Fluid Flow|
Sommeijer, B.P, & Verwer, J.G. (2007). On stabilized integration for time-dependent PDEs. Journal of Computational Physics, 224(1), 3–16.