Premium
A vector‐parallel scheme for Navier‐Stokes computations at multi‐gigaflop performance rates
Author(s) -
Lorber A. A.,
Carey G. F.
Publication year - 1995
Publication title -
international journal for numerical methods in fluids
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.938
H-Index - 112
eISSN - 1097-0363
pISSN - 0271-2091
DOI - 10.1002/fld.1650210602
Subject(s) - domain decomposition methods , stencil , mathematics , solver , incompressible flow , convergence (economics) , flow (mathematics) , algorithm , mathematical optimization , computational science , finite element method , geometry , physics , thermodynamics , economics , economic growth
A class of vector‐parallel schemes for solution of steady compressible or incompressible viscous flow is developed and performance studies carried out. The algorithms employ an artificial transient treatment that permits rapid integration to a steady state. In the present work a four‐stage explicit Runge‐Kutta scheme employing variable local step size is utilized for the ODE system integration. The RK‐4 scheme is restructured to allow vectorization and enhance concurrency in the calculation for a streamfunction‐vorticity formulation of the flow problem. The parameters of the resulting RK scheme can be selected to accelerate convergence of the RK recursion. Four main procedures are considered which permit vector‐parallel solution: a Jacobi update, a hybrid of the Jacobi and Gauss‐Seidel method, red‐black ordering and domain decomposition. Numerical performance studies are conducted with a representative viscous incompressible flow calculation. Results indicate that a scheme involving domain decomposition with a Gauss‐Seidel type of update for the RK four‐stage scheme is most effective and provides performance in excess of 8 Gflops on the Cray C‐90.