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.

Unified GPU Parallel Framework Based on Discontinuous Galerkin Method

Parallelized Implementation of the Finite Particle Method for Explicit Dynamics in GPU

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.

Advanced Parallelism of DGTD Method with Local Time Stepping Based on Novel MPI + MPI Unified Parallel Algorithm

GPU-Accelerated Discontinuous Galerkin Methods: 30x Speedup on 345 Billion Unknowns

New strategies for developing GPU accelerated disk‐based discontinuous deformation analysis for large‐scale modeling

Manycore Parallel Computing for a Hybridizable Discontinuous Galerkin Nested Multigrid Method

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

Wavefield Separation Algorithm of Helmholtz Theory-Based Discontinuous Galerkin Method Using Unstructured Meshes on GPU

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

An adaptive finite element multigrid solver using GPU acceleration

A scalable multi-GPU method for semi-implicit fractional-step integration of incompressible Navier-Stokes equations

A GPU-accelerated Cartesian grid method for PDEs on irregular domain

An implementation of MPI and hybrid OpenMP/MPI parallelization strategies for an implicit 3D DDG solver

A discontinuous Galerkin implementation of a domain decomposition method for kinetic-hydrodynamic coupling multiscale problems in gas dynamics and device simulations

Accelerating Dissipative Particle Dynamics Simulations on GPUs: Algorithms, Numerics and Applications

A Pipeline-Based ODE Solving Framework

p-adaptive discontinuous Galerkin method for the shallow water equations on heterogeneous computing architectures

Parallel Computing of Discrete Element Method on GPU

Massive parallelization and performance enhancement of an immersed boundary method based unsteady flow solver