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

  • Authors:
  • David L. Ropp;John N. Shadid

  • Affiliations:
  • Northrop Grumman Corporation, Electronic Systems, 7323 Aviation Blvd., Linthicum, MD 21240-2001, United States;Computational Sciences R&D Group, MS 0316, P.O. Box 5800, Sandia National Laboratories, Albuquerque, NM 87185-0316, United States

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2009

Quantified Score

Hi-index 31.46

Visualization

Abstract

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.