Implementation of a parallel ADI algorithm on a finite volume GPU-based elementary porous media flow computation

Abstract

In this work we present flow simulations in laminar and turbulent regime within a representative elementary volume of a simplified porous media by solving the Navier-Stokes equations and a Low-Re turbulence kÀ model. Numerical solution was achieved with an implementation of the SIMPLE algorithm for pressure velocity coupling of variables, and the solution of the tridiagonal systems of algebraic equations was accomplished by a parallelized ADI scheme based on the Thomas algorithm. Implementation of the numerical solution was done with an in-house C code which combined OMP and CUDA technologies for computations based on CPU and GPU, respectively. Exponential structured grids were employed in the wall vicinity to capture the turbulence behavior. Results indicate that similar profiles for velocity, pressure, turbulent kinetic energy and its dissipation were found. Several CUDA grids were tested and their performances measured over two GPUs: GTX 680 and GTX TITAN. Considerable speedup was achieved by the GPUs over the CPU schemes even without the use of the device shared memory which was not explored due to the nature of the algorithm.

Evaluation of Interpolation Scheme Alternatives and Variation of the Number of Slave Cells in a Cut-Cell Methodology

“…The SIMPLE algorithm of Patankar [25] was used to treat the pressure-velocity coupling. Fully implicit discretized equations were solved using a parallel implementation of the ADI/TDMA algorithm [28] in an in-house C/OMP code. All simulations were conducted until the normalized residues were below 10 −8 for the momentum and continuity equations.…”

Numerical investigation of supersonic transverse jet interaction on CPU/GPU system

“…In this study, GPU parallelization is considered. In the past 10 years, GPU has undergone significant updates [9,10]. Researchers have begun to develop the potential of GPUs in scientific computing.…”