Stability of operator splitting methods for systems with indefinite operators: reaction-diffusion systems

In this paper numerical results are reviewed [D.L. Ropp, J.N. Shadid, C.C. Ober, Studies of the accuracy of time integration methods for reaction-diffusion equations, J. Comput. Phys. 194 (2) (2004) 544-574] that demonstrate that common second-order operator-splitting methods can exhibit instabilities when integrating the Brusselator equations out to moderate times of about seven periods of oscillation. These instabilities are manifested as high wave number spatial errors. In this paper, we further analyze this problem, and present a theorem for stability of operator-splitting methods applied to linear reaction-diffusion equations with indefinite reaction terms which controls both low and high wave number instabilities. A corollary shows that if L-stable methods are used for the diffusion term the high wave number instability will be controlled more easily. In the absence of L-stability, an additional time step condition that suppresses the high wave number modes appears to guarantee convergence at the asymptotic order for the operator-splitting method. Numerical results for a model problem confirm this theory, and results for the Brusselator problem agree as well.

[1]  J. Verwer,et al.  Numerical solution of time-dependent advection-diffusion-reaction equations , 2003 .

[2]  J. Shadid,et al.  Studies of the Accuracy of Time Integration Methods for Reaction-Diffusion Equations ∗ , 2005 .

[3]  J. Brandts.  [Review of: W. Hundsdorfer, J.G. Verwer (2003) Numerical Solution of Time-Dependent Advection-Diffusion-Reaction Equations] , 2006 .

[4]  J. C. Jorge,et al.  Contractivity results for alternating direction schemes in Hilbert spaces , 1994 .

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

[6]  E. Hairer,et al.  Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems , 1993 .

[7]  Rüdiger Weiner,et al.  Linearly implicit splitting methods for higher space-dimensional parabolic differential equations , 1998 .

[8]  A. Hindmarsh,et al.  CVODE, a stiff/nonstiff ODE solver in C , 1996 .

[9]  J. Shadid,et al.  Studies on the accuracy of time-integration methods for the radiation-diffusion equations , 2004 .

[10]  N. N. Yanenko.  Application of the Method of Fractional Steps to Boundary Value Problems for Laplace’s and Poisson’s Equations , 1971 .

[11]  G. I. Marchuk,et al.  ON THE THEORY OF THE SPLITTING-UP METHOD , 1971 .

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

[13]  J. Verwer,et al.  Stability of Runge-Kutta Methods for Stiff Nonlinear Differential Equations , 1984 .

[14]  J. Lambert.  Numerical Methods for Ordinary Differential Systems: The Initial Value Problem , 1991 .

[15]  Willem Hundsdorfer,et al.  Stability and convergence of the Peaceman-Rachford ADI method for initial-boundary value problems , 1987 .

[16]  I. Prigogine,et al.  Symmetry Breaking Instabilities in Dissipative Systems. II , 1968 .

[17]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[18]  L. G. Stern,et al.  Fractional step methods applied to a chemotaxis model , 2000, Journal of mathematical biology.