Abstract:This paper presents a parallel implementation of the Discrete Unified Gas Kinetic Scheme (DUGKS) on the GPU system using the CUDA Fortran and CUDA C++ programming languages. Firstly, we conducted an extensive revision of our original CPU-based code, resulting in a threefold decrease in memory usage. This new implementation is also paired with a novel approach to compute cell face flux using trilinear interpolation. It is shown analytically that the interpolation-based approach to flux calculation is more accurate compared to the one used in the original DUGKS. The initial simulation results using this new approach suggest that trilinear interpolation can reduce numerical errors on a coarse mesh. For example, in the case of the decaying Taylor-Green vortex flow at a 128 3 mesh resolution, the relative numerical error in the energy dissipation rate at t⁎=2 , using the spectral simulation result as the benchmark, is approximately 30% lower than that of the original implementation. The improved GPU DUGKS method is applied to laminar and turbulent flows in periodic and wall-bounded boundary configurations. A performance comparison of the GPU implementation is also presented and compared to the previous CPU implementation. A maximum speedup of 7.64x was achieved on a desktop-level GPU compared to a 32-core CPU. The strong scaling test, conducted on an eight-GPU node, demonstrated the efficient utilization of available multiple GPU resources by the code. Program summary Program Title: DUGKS-GPU CPC Library link to program files: https://doi.org/10.17632/yykv5s9g2n.1 Developer's repository link: https://github.com/kairzhan/DUGKS-GPU Code Ocean capsule: https://codeocean.com/capsule/3b1e4f74-cdd3-4781-8923-8514d5923dfb/ Licensing provisions: GNU General Public License 3 Programming language: CUDA C++, CUDA Fortran Supplementary material: Nature of problem: DUGKS-GPU is an advanced numerical code designed to accelerate the simulation of turbulent fluid flows through the application of the kinetic method known as the Discrete Unified Gas Kinetic Scheme (DUGKS). This computational tool comprises a CUDA Fortran version, tailored for a single GPU, as well as a CUDA C++ version developed for multi-GPU platforms. Solution method: The code employs a second-order finite-volume discretization of the discrete velocity Boltzmann equation with the BGK collision model. It can be used to simulate small to medium-scale problems on modern multi-GPU platforms with CUDA architecture.

GPU acceleration of an iterative scheme for gas-kinetic model equations with memory reduction techniques

The Implementation of the Three-Dimensional Unified Gas-Kinetic Wave-Particle Method on Multiple Graphics Processing Units

Implementation of the moving particle semi-implicit method for free-surface flows on GPU clusters.

GPU Implementation of the Discrete Unified Gas Kinetic Scheme for Low-Speed Isothermal Flows

Gpu Speed-Up For The Implicit Navier-Stokes Solver

An adaptive finite element multigrid solver using GPU acceleration

Multiple-GPU accelerated high-order gas-kinetic scheme on three-dimensional unstructured meshes

Efficient Same-Dimensional Implicit Time Advancement Parallel Scheme and Optimization Methods for the Iteration Parameters Using a Graphics-Processing Unit

GPU implementation of algorithm SIMPLE-TS for calculation of unsteady, viscous, compressible and heat-conductive gas flows

DUGKS-GPU: An efficient parallel GPU code for 3D turbulent flow simulations using Discrete Unified Gas Kinetic Scheme

GPU-Accelerated Computations for Supersonic Flow Modeling on Hybrid Grids

GPU Optimization for High-Quality Kinetic Fluid Simulation

A Memory Reduction Compact Gas Kinetic Scheme on 3D Unstructured Meshes

Implicit method for the solution of supersonic and hypersonic 3D flow problems with Lower-Upper Symmetric-Gauss-Seidel preconditioner on multiple graphics processing units

A Multi-GPU Parallel Algorithm in Hypersonic Flow Computations

High Performance GPU Parallel Computing of Gas Dynamic by Direct Simulation Monte Carlo

Development of a Flow Solver with Complex Kinetics on the Graphic Processing Units

Accelerating Unstructured Large Eddy Simulation Solver with GPU

Multidisciplinary simulation acceleration using multiple shared memory graphical processing units

Optimization and acceleration of flow simulations for CFD on CPU/GPU architecture

GPU Acceleration of CFD Algorithm: HSMAC and SIMPLE.