Simulating flows of viscous incompressible fluid on graphics processors using the splitting scheme and multigrid method

The use of general purpose graphics processors for the numerical solution of problems in dynamics of viscous incompressible fluid is discussed. Specific features of the parallel implementation of the splitting scheme (projection method) are considered. The system of difference equations produced by the discretization of Poisson's equations for pressure is solved using the multigrid method. A number of benchmark problems are solved on graphics processors, and approaches to the optimization of program code by using different memory types are discussed. The speedup of computations on graphics processors is compared with the computations on the central processor using grids with different resolution and different decompositions of the initial data into blocks.