Solving 3D Viscous Navier-Stokes Equations Using CUDA

Santiago D. Costarelli, Mario A. Storti, Rodrigo R. Paz, Lisandro Dalcin, Sergio Idelsohn


A CUDA implementation of the 3D viscous incompressible Navier-Stokes equations is proposed using as advection operator the BFECC (Back and Forth Error Compensation and Correction) schema. The Poisson problem for pressure is solved with a CG (Conjugated Gradient) preconditioning the system with FFTs (Fast Fourier Transforms). Study cases as Lid-Driven Cavity and Flow Past Circular Cilinder, both 2D and 3D, are solved in order to check accuracy and obtain performance meassurements.

Full Text:


Asociación Argentina de Mecánica Computacional
Güemes 3450
S3000GLN Santa Fe, Argentina
Phone: 54-342-4511594 / 4511595 Int. 1006
Fax: 54-342-4511169
E-mail: amca(at)
ISSN 2591-3522