A novel ALE scheme with the internal boundary for true free surface simulation in geodynamic models
We propose a novel scheme for modelling the true free surface within the finite element method (FEM), which integrates the “sticky air” approach into the ALE scheme. This method employs an internal boundary to accurately represent the free surface, referred to as ALE-IB. We implement this scheme for free surface simulations in the geodynamic codes Underworld 2 and Underworld 3.
Classification of methods used for simulating a free surface (indicated by the magenta line). Colored points represent markers for different materials. Methods include: (a) ALE scheme, (b) ALE scheme with the internal boundary (ALE-IB) and the “sticky air” method, and (c) Eulerian scheme with the “sticky air” method. The real interface refers to the actual free surface, while the virtual interface represents the surface obtained from numerical modelling.
Research Tags
Associated Publication
A novel ALE scheme with the internal boundary for true free surface simulation in geodynamic models
Neng Lu, Louis Moresi, Julian Giordani
DOI10.5194/gmd-19-5191-2026
Abstract
The accurate simulation of Earth’s surface is essential for understanding lithospheric and mantle dynamics, especially in processes such as subduction and surface deformation. Traditional top boundary conditions, such as free-slip or no-slip, do not fully capture the complex interactions occurring at the surface. The commonly used “Sticky Air” method, while practical, suffers from several limitations, including increased computational cost and marker fluctuation issues. Additionally, free surface numerical fluctuations, known as the “drunken sailor instability”, are characteristic of all free surface simulations, including true Lagrangian free surface treatments and Arbitrary Lagrangian–Eulerian (ALE) methods. In this study, we propose a novel scheme within the finite element framework that integrates the “Sticky Air” concept into an ALE formulation by employing an internal boundary to simulate a true free surface, referred to as the ALE-IB. This approach effectively addresses the limitations of existing methods, notably by reducing marker fluctuation issues and enhancing numerical stability. Moreover, it maintains a true surface in the computational domain that can be further reshaped by surface processes such as erosion and deposition, and provides a foundational scheme for further coupling framework of tectonic modelling and landscape evolution modelling. We detail the theoretical formulation, implementation strategies, and validation through a series of numerical experiments. The results demonstrate that our method achieves higher accuracy and broader applicability compared to conventional techniques. Ultimately, this framework provides a more realistic and robust tool for geodynamic modelling of the Earth’s free surface.
Compute Tags
None specified.
Software
Software information not available.
Model Setup
Model setup for (a) viscous relaxation of sinusoidal topography, (b) Rayleigh–Taylor instability, (c) delamination, (d) rising sphere, and (e) subduction. Ω indicates different material domains, with Ω0 specifically representing the air domain – used only in Eulerian and ALE-IB schemes. The dashed line marks the free surface, and the stars denote tracer locations used in some experiments.
Source repository:
https://github.com/ModelAtlasofTheEarth/Lu-2026-ALEIB-FreeSurface