Abstract:This paper presents the parallelization of two widely used implicit numerical solvers for the solution of partial differential equations on structured meshes, namely, the ADI (Alternating-Direction Implicit) solver for tridiagonal linear systems and the SIP (Strongly Implicit Procedure) solver for the penta-diagonal systems. Both solvers were parallelized using CUDA (Computer Unified Device Architecture) Fortran on GPGPUs (General-Purpose Graphics Processing Units). The parallel ADI solver (P-ADI) is based on the Parallel Cyclic Reduction (PCR) algorithm, while the parallel SIP solver (P-SIP) uses the wave front method (WF) following a diagonal line calculation strategy. To map the solution schemes onto the hierarchical block-threads framework of the CUDA on the GPU, the P-ADI solver adopted two mapping methods, one block thread with iterations (OBM-it) and multi-block threads (MBMs), while the P-SIP solver also used two mappings, one conventional mapping using effective WF lines (WF-e) with matrix coefficients and solution variables defined on original computational mesh, and a newly proposed mapping using all WF mesh (WF-all), on which matrix coefficients and solution variables are defined. Both the P-ADI and the P-SIP have been integrated into a two-dimensional (2D) hydrodynamic model, the CCHE2D (Center of Computational Hydroscience and Engineering) model, developed by the National Center for Computational Hydroscience and Engineering at the University of Mississippi. This study for the first time compared these two parallel solvers and their efficiency using examples and applications in complex geometries, which can provide valuable guidance for future uses of these two parallel implicit solvers in computational fluids dynamics (CFD). Both parallel solvers demonstrated higher efficiency than their serial counterparts on the CPU (Central Processing Unit): 3.73~4.98 speedup ratio for flow simulations, and 2.166~3.648 speedup ratio for sediment transport simulations. In general, the P-ADI solver is faster than but not as stable as the P-SIP solver; and for the P-SIP solver, the newly developed mapping method WF-all significantly improved the conventional mapping method WF-e.

Development of an efficient large time step unsteady solver for incompressible flows using the IDEAL algorithm in OpenFOAM

High-fidelity Simulation of Regular Waves Based on Multi-Moment Finite Volume Formulation and THINC Method

Application of IDEAL algorithm based on the collocated unstructured grid for incompressible flows

Development of a Volume of Fluid Method for Computing Interfacial Incompressible Fluid Flows

Inflow and Outflow Numerical Simulation Using Least-Square Moving Particle Semi-Implicit Method on GPU

A PARALLEL SOLVING DESIGN FOR UNSTEADY NAVIER-STOKES EQUATION

An Efficient Semi-implicit Solver for Direct Numerical Simulation of Compressible Flows at All Speeds

Development of a Chemistry Dynamic Load Balancing Solver with Sparse Analytical Jacobian Approach for Rapid and Accurate Reactive Flow Simulations

A Scalable Fully Implicit Method with Adaptive Time Stepping for Unsteady Compressible Inviscid Flows

An Improved CIP-based Numerical Model for Simulating Free-Surface Flow with Adaptive Mesh

Implementation of IDEAL algorithm based on Delaunay triangular mesh for 2D-compressible flows

Implementation and assessment of a low-dissipative OpenFOAM solver for compressible multi-species flows

Numerical experiments using one high-resolution scheme for unsteady inviscid compressible flows

Development of low-dissipative projection method framework integrating various high-order time integration schemes using OpenFOAM

Accelerating flow simulations in the built environment by using the fast fluid dynamics initializer

Efficient implementation of complex equations of state in a high-order framework

Development of an accurate central finite-difference scheme with a compact stencil for the simulation of unsteady incompressible flows on staggered orthogonal grids

Research on Computational Method of Supersonic Inlet/Isolator Internal Flow

Parallel Implicit Solvers for 2D Numerical Models on Structured Meshes

Nonintrusive proper generalised decomposition for parametrised incompressible flow problems in OpenFOAM

Accurate and rapid reactive flow simulations using dynamic load balancing and sparse analytical Jacobian approach