Stability of operator splitting methods for systems with indefinite operators: Advection-diffusion-reaction systems
This brief paper presents an A-stability result for operator splitting type time integration methods applied to advection-diffusion-reaction equations with possibly indefinite source terms. These results extend our earlier work on diffusion-reaction systems [D.L. Ropp, J.N. Shadid, Stability of operator splitting methods for systems with indefinite operators: reaction-diffusion systems, J. Comput. Phys. 203 (2) (2005) 449-466]. The A-stability result presents sufficient conditions that control 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. Numerical results are presented that verify second-order convergence for the operator splitting methods and demonstrate control of instabilities on a chemotaxis problem by use of an L-stable diffusion integrator. (C) 2009 Elsevier Inc. All rights reserved.
operator splitting methods time integration advection-diffusion-reaction pdes linear stability numerical methods