Abdullah Mujahid

Robert Altmann, Abdullah Mujahid, Benjamin Unger

We study the construction and convergence of semi-explicit and iterative decoupling schemes for an elliptic-parabolic problem using higher-order Runge-Kutta methods. For the semi-explicit schemes, which are constructed using a nearby delay system with $k$ time delays, we establish the convergence of $k$th-order Runge-Kutta methods under a weak coupling condition. We develop the convergence analysis by adapting the Fourier stability and perturbation techniques of [Lubich, Ostermann, Math. Comp., 64(210):601--627, 1995]. The key tool is the generating function framework, in which the Runge-Kutta discretization is encoded through an operator-valued function. Stability estimates are then obtained via Parseval's identity on the unit circle. We further present convergence results for iterative (fixed-stress and undrained-split) higher-order Runge-Kutta schemes. Here, a spectral decomposition of the Schur complement operator is central. Finally, we provide numerical examples to verify the proven convergence results.

Robert Altmann, Abdullah Mujahid, Benjamin Unger

For the iterative decoupling of elliptic-parabolic problems such as poroelasticity, we introduce time discretization schemes up to order $5$ based on the backward differentiation formulae. Its analysis combines techniques known from fixed-point iterations with the convergence analysis of the temporal discretization. As main result, we show that the convergence depends on the interplay between the time step size and the parameters for the contraction of the iterative scheme. Moreover, this connection is quantified explicitly, which allows for balancing the single error components. Several numerical experiments illustrate and validate the theoretical results, including a three-dimensional example from biomechanics.