Efficient iterative methods applied to the solution of transonic flows
- Univ. of Minnesota, Minneapolis, MN (United States)
- Purdue Univ., West Lafayette, IN (United States)
- Wayne State Univ., Detroit, MI (United States)
We investigate the use of an inexact Newton`s method to solve the potential equations in the transonic regime. As a test case, we solve the two-dimensional steady transonic small disturbance equation. Approximate factorization/ADI techniques have traditionally been employed for implicit solutions of this nonlinear equation. Instead, we apply Newton`s method using an exact analytical determination of the Jacobian with preconditioned conjugate gradient-like iterative solvers for solution of the linear systems in each Newton iteration. Two iterative solvers are tested; a block s-step version of the classical Orthomin(k) algorithm called orthogonal s-step Orthomin (OSOmin) and the well-known GIVIRES method. The preconditioner is a vectorizable and parallelizable version of incomplete LU (ILU) factorization. Efficiency of the Newton-Iterative method on vector and parallel computer architectures is the main issue addressed. In vectorized tests on a single processor of the Cray C-90, the performance of Newton-OSOmin is superior to Newton-GMRES and a more traditional monotone AF/ADI method (MAF) for a variety of transonic Mach numbers and mesh sizes. Newton- GIVIRES is superior to MAF for some cases. The parallel performance of the Newton method is also found to be very good on multiple processors of the Cray C-90 and on the massively parallel thinking machine CM-5, where very fast execution rates (up to 9 Gflops) are found for large problems. 38 refs., 14 figs., 7 tabs.
- OSTI ID:
- 263368
- Journal Information:
- Journal of Computational Physics, Vol. 123, Issue 2; Other Information: PBD: Feb 1996
- Country of Publication:
- United States
- Language:
- English
Similar Records
Implementation of iterative methods for large sparse nonsymmetric linear systems on a parallel vector machine
User documentation for PVODE, an ODE solver for parallel computers