- 1Deutsches Zentrum für Luft- und Raumfahrt, Institut für Physik der Atmosphäre, Oberpfaffenhofen, Germany
- 2University of Cologne, Division of Mathematics, Cologne, Germany
- 3Leibniz Universität Hannover, Institute of Applied Mathematics, Hannover, Germany
- 4Deutsches Zentrum für Luft- und Raumfahrt, Simulations- und Softwaretechnik, Cologne, Germany
- 5Universtität Augsburg, Institute of Mathematics, Augsburg, Germany
The currently available computing power severely limits the spatial resolution in chemistry climate simulations, even on upcoming exascale machines. This is mainly due to the large number of prognostic variables which includes chemical tracers. To further enhance the reliability and accuracy of climate projections, small scales need to be better resolved.
Adaptive methods enable a continuously re-adjusted focus of computational power in time and space. This increases the achievable level of detail considerably, while reducing the time to solution and resource consumption. However, adaptivity requires a sophisticated selection of adaptation criteria, algorithms, memory layouts, and communication patterns to fully utilize modern HPC infrastructures.
Additionally, discontinuous Galerkin methods promise to increase the effective resolution, i.e. by employing high order polynomials, so that prognostic variables are better resolved, even on coarser meshes. Most of the additional computation is done locally, so that the overall algorithm is ideally suited for parallel execution. The typical lack of robustness of higher-order methods can be remedied by utilizing state-of-the-art entropy stable schemes.
In this conference contribution, we will present the setup and the interfaces between MESSy, Trixi.jl and t8code as well as results from a prototypical simulation that showcases the interaction, application, and challenges of dynamic adaptive meshes.
Here, MESSy is a software framework that allows to integrate multiple numerical model components to build regional and global chemistry climate models. t8code is a parallel mesh management library written in C++. Finally, Trixi.jl is a computational fluid dynamics solver, build around a modern Discontinuous Galerkin method, and written in the Julia programming language.
This work was performed within the project ADAPTEX (ADAPtive Earth system modelling with strongly reduced computation time for EXascale-supercomputers), which aims to evaluate the potential benefit of dynamic adaptive meshes in ESM.
How to cite: Hartung, K., Geihe, B., Hergl, C., Holke, J., Jöckel, P., Markert, J., and Schlottke-Lakemper, M.: Towards dynamic adaptive mesh refinement in Earth system models, EGU General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025, EGU25-15654, https://doi.org/10.5194/egusphere-egu25-15654, 2025.
