Preconditioning harmonic unsteady potential flow calculations

Journal article

This paper considers finite element discretizations of the Helmholtz equation and its generalization arising from harmonic acoustic perturbations to a nonuniform steady potential flow. A novel elliptic, positive-definite preconditioner with a multigrid implementation is used to accelerate the iterative convergence of Krylov subspace solvers. Both theory and numerical results show that for a model 1-D Helmholtz test problem, the preconditioner clusters the discrete system's eigenvalues and lowers its condition number to a level independent of grid resolution. For the 2-D Helmholtz equation, grid-independent convergence is achieved using a quasi-minimal residual Krylov solver, significantly outperforming the popular symmetric successive over-relaxation preconditioner. Impressive results are also presented on more complex domains, including an axisymmetric aircraft engine inlet with nonstagnant mean flow and modal boundary conditions. Copyright © 2006 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Authors: Laird, A; Giles, M

• Publication status: Published
• Journal: AIAA JOURNAL
• Volume: 44
• Issue: 11
• Pages: 2654-2662
• Publication date: 2006-11-01
• DOI: 10.2514/1.15243
• EISSN: 1533-385X
• ISSN: 0001-1452
• Language: English