Articles | Volume 19, issue 4
https://doi.org/10.5194/gmd-19-1455-2026
https://doi.org/10.5194/gmd-19-1455-2026
Model evaluation paper
 | 
19 Feb 2026
Model evaluation paper |  | 19 Feb 2026

Highly scalable geodynamic simulations with HyTeG

Ponsuganth Ilangovan, Nils Kohl, and Marcus Mohr

Cited articles

Bangerth, W., Dannberg, J., Fraters, M., Gassmoeller, R., Glerum, A., Heister, T., Myhill, R., and Naliboff, J.: geodynamics/aspect: ASPECT 2.5.0, Zenodo [code], https://doi.org/10.5281/ZENODO.8200213, 2023. a
Bauer, S., Mohr, M., Rüde, U., Weismüller, J., Wittmann, M., and Wohlmuth, B.: A two-scale approach for efficient on-the-fly operator assembly in massively parallel high performance multigrid codes, Appl. Numer. Math., 122, 14–38, https://doi.org/10.1016/j.apnum.2017.07.006, 2017. a
Bauer, S., Bunge, H.-P., Drzisga, D., Ghelichkhan, S., Huber, M., Kohl, N., Mohr, M., Rüde, U., Thönnes, D., and Wohlmuth, B.: TerraNeo–Mantle Convection Beyond a Trillion Degrees of Freedom, Springer International Publishing, 569–610, ISBN 9783030479565, https://doi.org/10.1007/978-3-030-47956-5_19, 2020. a
Baumgardner, J.: Three-dimensional treatment of convective flow in the earth's mantle, J. Stat. Phys., 39, 501–511, 1985. a
Baumgardner, J. R. and Frederickson, P. O.: Icosahedral discretization of the two-sphere, SIAM J. Numer. Anal., 22, 1107–1115, https://doi.org/10.1137/0722066, 1985. a
Download
Short summary
Geophysically accurate model of mantle convection requires a mesh with a width on the order of ~1 km. Traditional codes represent the systems of equation by setting up the associated matrix. However, at the scales we want to operate, even forming this matrix is hardly possible. Thus, we use the matrix-free framework HyTeG to create the geophysical model and verify it with numerical experiments while assessing the scalability of the framework and laying out the difficulties involved.
Share