Development and acceleration of unstructured mesh-based cfd solver

¹ Faculty of Power Engineering, Baltic State Technical University 1, 1st Krasnoarmeyskaya Str., St. Petersburg 190005, Russia
² Faculty of Mathematics and Mechanics, St. Petersburg State University Universitetsky Prosp., Old Petergof, St. Petersburg 198504, Russia
³ Faculty of Science, Engineering and Computing, Kingston University Friars Av., Roehampton Vale, London SW15 3DW, U.K.

Abstract

The study was undertaken as part of a larger effort to establish a common computational fluid dynamics (CFD) code for simulation of internal and external flows and involves some basic validation studies. The governing equations are solved with ¦nite volume code on unstructured meshes. The computational procedure involves reconstruction of the solution in each control volume and extrapolation of the unknowns to find the flow variables on the faces of control volume, solution of Riemann problem for each face of the control volume, and evolution of the time step. The nonlinear CFD solver works in an explicit time-marching fashion, based on a three-step Runge-Kutta stepping procedure. Convergence to a steady state is accelerated by the use of geometric technique and by the application of Jacobi preconditioning for high-speed flows, with a separate low Mach number preconditioning method for use with low-speed flows. The CFD code is implemented on graphics processing units (GPUs). Speedup of solution on GPUs with respect to solution on central processing units (CPU) is compared with the use of different meshes and different methods of distribution of input data into blocks. The results obtained provide promising perspective for designing a GPU-based software framework for applications in CFD.