Author: Brandén Henrik Holmgren Sverker Sundqvist Per
Publisher: Springer Publishing Company
ISSN: 0885-7474
Source: Journal of Scientific Computing, Vol.30, Iss.1, 2007-01, pp. : 35-60
Disclaimer: Any content in publications that violate the sovereignty, the constitution or regulations of the PRC is not accepted or approved by CNPIEC.
Abstract
We present a new preconditioner for the iterative solution of systems of equations arising from discretizations of systems of first order partial differential equations (PDEs) on structured grids. Such systems occur in many important applications, including compressible fluid flow and electromagnetic wave propagation. The preconditioner is a truncated convolution operator, with a kernel that is a fundamental solution of a difference operator closely related to the original discretization. Analysis of a relevant scalar model problem in two spatial dimensions shows that grid independent convergence is obtained using a simple one-stage iterative method. As an example of a more involved problem, we consider the steady state solution of the non-linear Euler equations in a two-dimensional, non-axisymmetric duct. We present results from numerical experiments, verifying that the preconditioning technique again achieves grid independent convergence, both for an upwind discretization and for a centered second order discretization with fourth order artificial viscosity.
Related content
On “Superstable” discrete systems
Automation and Remote Control, Vol. 68, Iss. 4, 2007-04 ,pp. :