Operator-Based Preconditioning of Stiff Hyperbolic Systems
Authors: Daniel R. Reynolds, Ravi Samtaney, and Carol S. Woodward
Abstract:
We introduce an operator-based scheme for preconditioning stiff components encoun- tered in implicit methods for hyperbolic systems of partial differential equations posed on regular grids. The method is based on a directional splitting of the implicit operator, followed by a char- acteristic decomposition of the resulting directional parts. This approach allows for solution to any number of characteristic components, from the entire system to only the fastest, stiffness-inducing waves. We apply the preconditioning method to stiff hyperbolic systems arising in magnetohydro- dynamics and gas dynamics. We then present numerical results showing that this preconditioning scheme works well on problems where the underlying stiffness results from the interaction of fast transient waves with slowly-evolving dynamics, scales well to large problem sizes and numbers of processors, and allows for additional customization based on the specific problems under study.
__________________________________________________
Submitted to: SIAM Journal of Scientific Computing
__________________________________________________
Download PPPL-4377 (pdf 905KB 20pp)
__________________________________________________