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 Computational Hydrodynamics Solver PNAOE-SJTU

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

Research on the single-host parallel computing with the local time step scheme for modeling of hydro-sediment-morphodynamic processes

Parallel Implementations of Numerical Simulation of Wind Flow Around Buildings

Parallel finite element method for 2D free surface seepage

Parallel Computing for Water Quality Equation Based on Jppf

Parallel Implicit Solvers for 2D Numerical Models on Structured Meshes

Naoe-Foam-sjtu SOLVER FOR SHIP FLOWS AND OCEAN ENGINEERING FLOWS

A PARALLEL SOLVING DESIGN FOR UNSTEADY NAVIER-STOKES EQUATION

Development of Naoe-Foam-sjtu Solver Based on OpenFOAM for Marine Hydrodynamics

A Parallel Strategy Applied to the Simplified Thermal-Hydraulic Part of System Program

An efficient parallel algorithm with application to computational fluid dynamics

SPHydro: Promoting smoothed particle hydrodynamics method toward extensive applications in ocean engineering

A Parallel Unstructured Finite-Volume Method for All-Speed Flows

A Parallel JPWSPC Algorithm for Hydrodynamic Simulation of River Network

Massively Parallel Simulation of the Flow Characteristics in A Submersible Axial Flow Pump

A parallel solver framework for fully implicit monolithic fluid-structure interaction

THE IMPLEMENTATION OF A PARALLEL SOFTWARE OF PETROLEUM RESERVOIR SIMULATION AND ITS APPLICATION ON NATIVE HIGH PERFORMANCE COMPUTERS

A parallel solver for fluid structure interaction problems with Lagrange multiplier

A Parallel Direct Numerical Method for Solving Unsteady Flow Issue of Tree-type River Networks