Stability of operator splitting methods for systems with indefinite operators: Advectiondiffusionreaction systems ?
This brief paper presents an A-stability result for operator splitting type time integration methods applied to advectiondiffusionreaction equations with possibly indefinite source terms. These results extend our earlier work on diffusionreaction systems [D.L. Ropp, J.N. Shadid, Stability of operator splitting methods for systems with indefinite operators: reactiondiffusion systems, J. Comput. Phys. 203 (2) (2005) 449466]. 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.
operator splitting methods time integration advectiondiffusionreaction pdes linear stability numerical methods