2025/06/12 by Igor Tominec, Lukas Lundgren, Tominec, Igor +5
Engineering · Mathematics · #Advanced Numerical Methods in Computational Mathematics #Navier-Stokes equation solutions #Stability and Controllability of Differential Equations
paper · pdf · doi:10.48550/arxiv.2506.10447
Accurate simulations of ice sheet dynamics, mantle convection, lava flow, and other highly viscous free-surface flows involve solving the coupled Stokes/free-surface equations. In this paper, we theoretically analyze the stability and conservation properties of the weak form of this system for Newtonian fluids and non-Newtonian fluids, at both the continuous and discrete levels. We perform the fully discrete stability analysis for finite element methods used in space with explicit and implicit Euler time-stepping methods used in time. Motivated by the theory, we propose a stabilization term designed for the explicit Euler discretization, which ensures unconditional time stability and permits conservation of the domain volume. Numerical experiments validate and support our theoretical findings.