A FFT Preconditioning Technique for the Solution of Incompressible Flow on GPU's

Abstract

The resolution of Computational Fluid Dynamics (CFD) problems on Graphic Processing Units (GPU's) requires of specialized algorithms due to the particular hardware architecture of these devices. Algorithms that fall in the category of cellular automata (CA) are the best fitted, for instance explicit Finite Volume or Finite Element methods. But in the case of incompressible flow it is not possible to develop a pure explicit algorithm, due to the essentially non-local character of the incompressibility condition. In this case the algorithms that are closer to an explicit approach, are segregated algorithms, like the Fractional Step Method. In these algorithms the more time consuming stage is (asymptotically for large problems) the solution of the Poisson's equation for pressure. A common choice for it's solution is the IOP (Iterated Orthogonal Projection) method, which requires a series of solutions on the complete mesh. In this work a variant of the IOP, called Accelerated Global Preconditioning (AGP), is proposed. It is based on using a Preconditioned Conjugate Gradient (which is an accelerated iterative method, in contrast with the stationary scheme used in IOP) for the pressure on the fluid, and preconditioning with the solution on the global domain (fluid and solid). [Slides of talk given at CIMNE (Barcelona, Spain), 2012-04-11]