Abstract
The generation of detailed three dimensional meshes for the simulation of groundwater flow in thin layered domains is crucial to capture important properties of the underlying domains and to reach a satisfying accuracy. At the same time, this level of detail poses high demands both on suitable hardware and numerical solver efficiency. Parallel multigrid methods have been shown to exhibit near optimal weak scalability for massively parallel computations of density driven flow. A fully automated parameterized algorithm for prism based meshing of coarse grids from height data of individual layers is presented. Special structures like pinch outs of individual layers are preserved. The resulting grid is used as a starting point for parallel mesh and hierarchy creation through interweaved projected refinement and redistribution. Efficiency and applicability of the proposed approach are demonstrated for a parallel multigrid based simulation of a realistic sample problem.
Original language | English (US) |
---|---|
Journal | Computing and Visualization in Science |
Volume | 23 |
Issue number | 1-4 |
DOIs | |
State | Published - Apr 20 2020 |
Bibliographical note
KAUST Repository Item: Exported on 2020-10-01Acknowledgements: This work has been supported by the German Ministry of Economics and Technology (BMWi, 02E11476B) and by the DFG Priority Program 1648 Software for Exascale Computing (SPPEXA) in the project Exasolvers (WI 1037/24-2). We thank the HLRS for the opportunity to use Hazel Hen and their kind support. The authors also gratefully acknowledge the Gauss Centre for Supercomputing e.V. (www.gauss-centre.eu) for funding this project by providing computing time through the John von Neumann Institute for Computing (NIC) on the GCS Supercomputer JUQUEEN at Jülich Supercomputing Centre (JSC).