ELEFANT: a user-friendly multipurpose geodynamics code
- 1Department of Earth Sciences, University of Utrecht, Budapestlaan 4, 3584 CD, Utrecht, the Netherlands
- 2Centre for Earth Evolution and Dynamics (CEED), Postbox 1048, Blindern 0316 Oslo, Norway
Abstract. A new finite element code for the solution of the Stokes and heat transport equations is presented. It has purposely been designed to address geological flow problems in two and three dimensions at crustal and lithospheric scales. The code relies on the Marker-in-Cell technique and Lagrangian markers are used to track materials in the simulation domain which allows recording of the integrated history of deformation; their (number) density is variable and dynamically adapted. A variety of rheologies has been implemented including nonlinear thermally activated dislocation and diffusion creep and brittle (or plastic) frictional models. The code is built on the Arbitrary Lagrangian Eulerian kinematic description: the computational grid deforms vertically and allows for a true free surface while the computational domain remains of constant width in the horizontal direction. The solution to the large system of algebraic equations resulting from the finite element discretisation and linearisation of the set of coupled partial differential equations to be solved is obtained by means of the efficient parallel direct solver MUMPS whose performance is thoroughly tested, or by means of the WISMP and AGMG iterative solvers. The code accuracy is assessed by means of many geodynamically relevant benchmark experiments which highlight specific features or algorithms, e.g., the implementation of the free surface stabilisation algorithm, the (visco-)plastic rheology implementation, the temperature advection, the capacity of the code to handle large viscosity contrasts. A two-dimensional application to salt tectonics presented as case study illustrates the potential of the code to model large scale high resolution thermo-mechanically coupled free surface flows.
11 citations as recorded by crossref.
- Benchmarking numerical models of brittle thrust wedges S. Buiter et al. 10.1016/j.jsg.2016.03.003
- On the choice of finite element for applications in geodynamics C. Thieulot & W. Bangerth 10.5194/se-13-229-2022
- Nonlinear viscoplasticity in ASPECT: benchmarking and applications to subduction A. Glerum et al. 10.5194/se-9-267-2018
- SEPRAN: A versatile finite-element package for a wide variety of problems in geosciences A. van den Berg et al. 10.1007/s12583-015-0508-0
- Geodynamic diagnostics, scientific visualisation and StagLab 3.0 F. Crameri 10.5194/gmd-11-2541-2018
- A community benchmark for viscoplastic thermal convection in a 2-D square box N. Tosi et al. 10.1002/2015GC005807
- The effect of obliquity on temperature in subduction zones: insights from 3-D numerical modeling A. Plunder et al. 10.5194/se-9-759-2018
- 101 geodynamic modelling: how to design, interpret, and communicate numerical studies of the solid Earth I. van Zelst et al. 10.5194/se-13-583-2022
- Using the level set method in geodynamical modeling of multi-material flows and Earth's free surface B. Hillebrand et al. 10.5194/se-5-1087-2014
- Analytical solution for viscous incompressible Stokes flow in a spherical shell C. Thieulot 10.5194/se-8-1181-2017
- Metamorphic Facies and Deformation Fabrics Diagnostic of Subduction: Insights From 2D Numerical Models A. Regorda et al. 10.1029/2021GC009899