Globalized Newton-Krylov-Schwarz algorithms and software for parallel implicit CFD.
Implicit solution methods are important in applications modeled by PDEs with disparate temporal and spatial scales. Because such applications require high resolution with reasonable turnaround, parallelization is essential. The pseudo-transient matrix-free Newton-Krylov-Schwarz ({psi}NKS) algorithmic framework is presented as a widely applicable answer. This article shows that for the classical problem of three-dimensional transonic Euler flow about an M6 wing, {psi}NKS can simultaneously deliver globalized, asymptotically rapid convergence through adaptive pseudo-transient continuation and Newton's method; reasonable parallelizability for an implicit method through deferred synchronization and favorable communication-to-computation scaling in the Krylov linear solver; and high per processor performance through attention to distributed memory and cache locality, especially through the Schwarz preconditioner. Two discouraging features of {psi}NKS methods are their sensitivity to the coding of the underlying PDE discretization and the large number of parameters that must be selected to govern convergence. The authors therefore distill several recommendations from their experience and reading of the literature on various algorithmic components of {psi}NKS, and they describe a freely available MPI-based portable parallel software implementation of the solver employed here.
- Research Organization:
- Argonne National Lab. (ANL), Argonne, IL (United States)
- Sponsoring Organization:
- USDOE Office of Science (SC); USDOD; National Science Foundation (NSF)
- DOE Contract Number:
- DE-AC02-06CH11357
- OSTI ID:
- 942826
- Report Number(s):
- ANL/MCS/JA-34971; TRN: US200923%%623
- Journal Information:
- Int. J. High Perform. Comput. Appl., Vol. 14, Issue 2 ; 2000; ISSN 1078-3482
- Country of Publication:
- United States
- Language:
- ENGLISH
Similar Records
Parallel Newton-Krylov-Schwarz algorithms for the transonic full potential equation
Newton-Krylov-Schwarz algorithms for the 2D full potential equation