An incompressible flow solver on a GPU/CPU heterogeneous architecture parallel computing platform

A computational fluid dynamics (CFD) solver for a GPU/CPU heterogeneous architecture parallel computing platform is developed to simulate incompressible flows on billion-level grid points. To solve the Poisson equation, the conjugate gradient method is used as a basic solver, and a Chebyshev method in combination with a Jacobi sub-preconditioner is used as a preconditioner. The developed CFD solver shows good performance on parallel efficiency, which exceeds $90\%$ in the weak-scalability test when the number of grid points allocated to each GPU card is greater than ${208}^3$. In the acceleration test, it is found that running a simulation with ${1040}^3$ grid points on 125 GPU cards accelerates by 203.6x over the same number of CPU cores. The developed solver is then tested in the context of a two-dimensional lid-driven cavity flow and three-dimensional Taylor-Green vortex flow. The results are consistent with previous results in the literature.