CG+FFT Poisson solver for CFD problems

DALCÍN LISANDRO; STORTI M.; PAZ RODRIGO; COSTARELLI, SANTIAGO

Buenos Aires

Congreso; MECOM 2010 - IX Congreso Argentino de Mecánica Computacional. II Congreso Sudamericano de Mecánica Computacional. CILAMCE 2010 - XXXI Congreso Ibero-Latino-Americano de Métodos Computacionales en la Ingeniería; 2010

Asociación Argentina de Mecánica Computacional AMCA

 Processing Units (GPU´s) have received much attention in last years. Compute-intensive algorithms operating on multidimensional arrays that have nearest neighbor dependency and/or exploit data locality can achieve massive speedups. Simulation of problems modeled by time-dependent Partial Differential Equations (PDE´s) by using explicit time-stepping methods on structured grids is an instance of such GPU-friendly algorithms. Solvers for transient incompressible fluid flow cannot be developed in a fully explicit manner due to the incompressibility constraint. Fractional step methods require the solution of a Poisson problem for the pressure field at each time level. This stage is usually the most time-consuming one. This work discuss a Poisson solver for the pressure problem on structured grids and fluid+solid domains. This solver is based on standard Conjugate Gradients (CG) iterations and Fast Fourier Transform (FFT) preconditioning. The implementation targets NVIDIA CUDA architecture and takes advantage of freely available software components like *Thrust*, *Cusp*, and *CUFFT* libraries.