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.

Parallel Multigrid Computations under NOWs

A Parallel Simulator for Massive Reservoir Models Utilizing Distributed-Memory Parallel Systems

Parallel Algorithms for Moving Boundary Problems by Local Remeshing

Power Grid Analysis and Optimization Using Algebraic Multigrid

A Parallel Function Evaluation Approach for Solution to Large-Scale Equation-Oriented Models.

A Two-Level Parallel Decomposition Approach for Transient Stability Constrained Optimal Power Flow

A Parallel Method With Hybrid Algorithms For Mixed Integer Nonlinear Programming

A Realistic Parallel Computational Model

Parallel Computation of Compressible Flows for Implicit Schemes Without Sub-iteration and Modification of the Implicit Solver

Helmholtz Solving and Performance Optimization in Global/Regional Assimilation and Prediction System

Improved parallel finite element methods for the stationary Navier–Stokes problem

Parallel Algorithms for Dynamics Analysis of Structures and Their Impementations

A Parallel I/O System Based on NOW

Parallel Implicit Solvers for 2D Numerical Models on Structured Meshes

Parallelization of a grid-to-grid routing algorithm based on grids layering

Parallel Smoothers for Matrix-Based Geometric Multigrid Methods on Locally Refined Meshes Using Multicore CPUs and GPUs

A Parallel Simulation Framework Incorporating Machine Learning-Based Hotspot Detection for Accelerated Power Grid Analysis

Heterogeneous Parallel Algorithm Design and Performance Optimization for WENO on the Sunway TaihuLight Supercomputer

An efficient parallel algorithm with application to computational fluid dynamics

Implicit Parallel FEM Analysis of Shallow Water Equations

Analysis of a New Space-Time Parallel Multigrid Algorithm for Parabolic Problems