Block-structured, equal-workload, multi-grid-nesting interface for the Boussinesq wave model FUNWAVE-TVD (Total Variation Diminishing)

• Published in: Geoscientific Model Development Vol. 15; no. 14; pp. 5441 - 5459
• Main Authors: Choi, Young-Kwang, Shi, Fengyan, Malej, Matt, Smith, Jane M, Kirby, James T, Grilli, Stephan T
• Format: Journal Article
• Language: English
• Published: Katlenburg-Lindau Copernicus GmbH 18.07.2022; Copernicus Publications
• Subjects: Algorithms; Allocations; Analysis; Boundary conditions; Boussinesq approximation; Boussinesq equations; Central processing units; Coastal effects; Communication; CPUs; Data management; Earthquakes; Efficiency; Grid refinement (mathematics); Message passing; Microprocessors; Model accuracy; Modelling; Nesting; Propagation; Simulation; Solvers; Tsunamis; Workload; Workloads; United States > US; California
• ISSN: 1991-9603; 1991-959X; 1991-962X; 1991-9603; 1991-962X
• DOI: 10.5194/gmd-15-5441-2022

We describe the development of a block-structured, equal-CPU-load (central processing unit), multi-grid-nesting interface for the Boussinesq wave model FUNWAVE-TVD (Fully Nonlinear Boussinesq Wave Model with Total Variation Diminishing Solver). The new model framework does not interfere with the core solver, and thus the core program, FUNWAVE-TVD, is still a standalone model used for a single grid. The nesting interface manages the time sequencing and two-way nesting processes between the parent grid and child grid with grid refinement in a hierarchical manner. Workload balance in the MPI-based (message passing interface) parallelization is handled by an equal-load scheme. A strategy of shared array allocation is applied for data management that allows for a large number of nested grids without creating additional memory allocations. Four model tests are conducted to verify the nesting algorithm with assessments of model accuracy and the robustness in the application in modeling transoceanic tsunamis and coastal effects.