Abstract:The shallow water flows in coastal areas of seas, rivers and reservoirs are simulated usually by 2-D depth averaged models. However, the needs for fine resolution of the computational grids and large scales of the modeling areas require in practical applications to use the algorithms and hardware of HPC. We present comparison of the computational efficiency of the developed parallel 2-D modeling system COASTOX on CPU based multi-processor systems and GPUs.The hydrodynamic module of COASTOX is based on nonlinear shallow water equations (SWE), which describe currents and long waves, including tsunami, river flood waves and wake waves, generated by big vessels in shallow coastal areas. The special pressure term in momentum equations depending from the form of the draft of the vessel is used for wave generation by moving vessels. The currents in the marine nearshore areas generated by wind waves are described by the including into the SWE the wave-radiation stress terms. Sediment and pollutant transport are described by the 2-D advection-diffusion equations with the sink-source terms describing sedimentation-erosion and water-bottom contaminate exchange.Model equations are solved by finite volume method on rectangular grids or unstructured grids with triangular cells. Solution scheme of SWE is Godunov-type, explicit, conservative, has TVD property. Second order in time and space is achieved by Runge-Kutta predictor-corrector method and using different methods for calculating fluxes at predictor and corrector steps. Transport equations schemes are simple upwind and have first order in time and space.Model parallelized for computations on multi-core CPU systems based on domain decomposition approach with halo boundary structures and message-passing updating. To decompose an unstructured model grid, METIS graph partition library is used. For halo values updating the MPI technology is implemented with using of non-blocking send and receive functions.For computations on GPU the model is parallelized using OpenACC directive-based programming interface. Numerical schemes of the model are implemented in the form of loops for cells, nodes, faces with independent iterations because of scheme explicitness and locality. So, OpenACC directives inserted in model code specify for compiler the loops that may be computed in parallel.The efficiency of the developed parallel algorithms is demonstrated for CPU and GPU computing systems by such applications:<ol><li>Simulation of river flooding of July 2008 extreme flood on Prut river (Ukraine).</li> <li>Modeling of ship waves caused by tanker passage on the San Jacinto river near Barbours Cut Container Terminal (USA) and loads on moored container ship.</li> <li>Simulation of the consequences of the breaks of the dikes constructed on the heavy contaminated floodplain of the Pripyat River upstream Chernobyl Nuclear Power Plant.</li> </ol>For parallel performance testing we use Dell 7920 Workstation with 2 Intel Xeon Gold 6230 20 cores processors and NVIDIA Quadro RTX 5000 GPU. We obtain that multi-core computation up to 17.3 times faster than single core with parallel efficiency 43%. And for big computational grid (about or more than a million nodes) GPU faster than single core in 47.5-79.6 times and faster than workstation in 3-4.6 times.

Performance comparison of the A-grid and C-grid shallow-water models on icosahedral grids

HyGrid: A CPU-GPU Hybrid Convolution-Based Gridding Algorithm in Radio Astronomy.

CICE on a C-grid: new momentum, stress, and transport schemes for CICEv6.5

Improving Parallel Performance of A Finite-Difference Agcm on Modern High-Performance Computers

A new adaptive Cartesian-grid CIP method for computation of violent free-surface flows

Comparison of sea ice kinematics at different resolutions modeled with a grid hierarchy in the Community Earth System Model (version 1.2.1)

High-resolution subgrid models: background, grid generation, and implementation

Design of High-Resolution Grid Systems for the Next Generation Coupled Climate System Models and Earth System Models

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

P-CSI V1.0, an Accelerated Barotropic Solver for the High-Resolution Ocean Model Component in the Community Earth System Model V2.0

Multi-Scale Sea Ice Kinematics Modeling with a Grid Hierarchy in Community Earth System Model (version 1.2.1)

Development of a Hybrid Parallel MCV-based High-Order Global Shallow-Water Model

A Scalable Helmholtz Solver in GRAPES over Large‐scale Multicore Cluster

AGCM-3DLF: Accelerating Atmospheric General Circulation Model via 3D Parallelization and Leap-Format

An Adaptive Multimoment Global Model on A Cubed Sphere

High performance computing of waves, currents and contaminants in rivers and coastal areas of seas on multi-processors systems and GPUs

Accuracy and stability analysis of horizontal discretizations used in unstructured grid ocean models

Achieving high performance and portable parallel GMRES algorithm for compressible flow simulations on unstructured grids

Parallel Implicit Solvers for 2D Numerical Models on Structured Meshes

Semi-Lagrangian shallow water equations solver on the cubed-sphere grid as a prototype of new-generation global atmospheric model