pith. sign in

arxiv: 2605.22485 · v1 · pith:4WWYZOJSnew · submitted 2026-05-21 · 🧮 math.NA · cs.NA

Decoupling Runge-Kutta schemes for elliptic-parabolic problems

Robert Altmann , Abdullah Mujahid , Benjamin Unger This is my paper

Pith reviewed 2026-05-22 03:25 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Runge-Kutta methodsdecoupling schemeselliptic-parabolic problemsconvergence analysissemi-explicit schemesiterative decouplinggenerating functions
0
0 comments X

The pith

Semi-explicit Runge-Kutta schemes converge at full order for elliptic-parabolic problems under weak coupling.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops semi-explicit decoupling schemes for elliptic-parabolic problems by replacing the original system with a nearby delay system that has k time delays. It proves that applying a kth-order Runge-Kutta method to this delayed system yields convergence of order k provided a weak coupling condition is met. The proof uses a generating function to represent the discretization and applies Fourier stability analysis with Parseval's identity. Convergence is also established for iterative decoupling schemes such as fixed-stress and undrained-split methods through spectral analysis of the Schur complement. Numerical examples illustrate the predicted convergence rates.

Core claim

We establish the convergence of kth-order Runge-Kutta methods under a weak coupling condition for semi-explicit schemes constructed using a nearby delay system with k time delays. The analysis adapts Fourier stability and perturbation techniques using the generating-function framework, obtaining stability estimates via Parseval's identity on the unit circle. Convergence results for iterative higher-order Runge-Kutta schemes follow from a spectral decomposition of the Schur complement operator.

What carries the argument

The nearby delay system with k time delays used to construct the semi-explicit decoupling scheme, which permits the application of standard Runge-Kutta methods while maintaining order under the weak coupling condition.

If this is right

  • Decoupling becomes feasible for higher-order time integrators without sacrificing accuracy in weakly coupled elliptic-parabolic systems.
  • The generating function approach yields stability bounds directly from the unit circle integral.
  • Iterative schemes achieve the same order through spectral properties of the Schur complement.
  • Practical computations can use these schemes to avoid solving coupled systems at each stage.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The weak coupling condition likely holds in many poroelastic or fluid flow applications, enabling broader use of these schemes.
  • Similar delay-based constructions could apply to other coupled PDE systems where full coupling is expensive.
  • Testing the method on nonlinear extensions would check if the linear analysis carries over.

Load-bearing premise

The Fourier stability and perturbation techniques can be adapted to the elliptic-parabolic setting using the generating-function framework without further restrictions beyond the weak coupling condition.

What would settle it

Run a simulation with a known exact solution for an elliptic-parabolic problem where the coupling satisfies the weak condition and check if the observed convergence rate equals the theoretical order k for the chosen Runge-Kutta method.

read the original abstract

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.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript constructs semi-explicit decoupling schemes for elliptic-parabolic problems by approximating the system with a nearby delay equation involving k time delays. It proves that kth-order Runge-Kutta methods achieve full convergence order under a weak coupling condition by adapting the generating-function framework, Fourier stability analysis, and perturbation arguments of Lubich and Ostermann (1995), with stability obtained via Parseval's identity on the unit circle. Convergence results are also given for iterative schemes (fixed-stress and undrained-split) that rely on a spectral decomposition of the Schur complement operator. Numerical examples are provided to illustrate the theoretical rates.

Significance. If the central adaptation succeeds, the work supplies a rigorous justification for high-order time-stepping with decoupling in elliptic-parabolic systems, which appear in poroelasticity and related applications. The explicit use of an established 1995 technique together with a verifiable weak coupling condition and supporting numerics constitutes a concrete contribution to the numerical analysis of coupled evolution problems.

major comments (1)
  1. [Convergence analysis for semi-explicit schemes] In the convergence analysis for the semi-explicit schemes (the section adapting Lubich-Ostermann via the operator-valued generating function), the manuscript states that the weak coupling condition suffices to inherit the required resolvent bounds. However, an explicit verification is needed that the perturbation arising from the elliptic operator and the Schur complement remains O(Δt^k) uniformly with respect to the spatial mesh size and any elliptic regularity parameter; without this uniform bound the contour-integral or maximum-modulus argument may lose uniformity and the claimed order may suffer mesh-dependent degradation.
minor comments (2)
  1. [Abstract] The abstract refers to 'k time delays' without immediately relating k to the stage order or the number of stages of the Runge-Kutta method; a brief clarifying sentence would improve readability.
  2. [Numerical examples] In the numerical section, the convergence plots would be clearer if the theoretical slope lines (e.g., O(Δt^k)) were superimposed on the error curves for direct visual comparison.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comment. The positive assessment of the work's significance is appreciated. We address the major comment below and will revise the manuscript to incorporate an explicit verification as suggested.

read point-by-point responses

  1. Referee: In the convergence analysis for the semi-explicit schemes (the section adapting Lubich-Ostermann via the operator-valued generating function), the manuscript states that the weak coupling condition suffices to inherit the required resolvent bounds. However, an explicit verification is needed that the perturbation arising from the elliptic operator and the Schur complement remains O(Δt^k) uniformly with respect to the spatial mesh size and any elliptic regularity parameter; without this uniform bound the contour-integral or maximum-modulus argument may lose uniformity and the claimed order may suffer mesh-dependent degradation.

    Authors: We agree that an explicit verification of the uniform bound is necessary for full rigor. The weak coupling condition together with the spectral decomposition of the Schur complement ensures that the perturbation term is O(Δt^k) with a constant independent of the spatial mesh size h and elliptic regularity parameters. This follows from the operator-valued generating function representation and the application of Parseval's identity, where the elliptic contribution is controlled uniformly by the weak coupling assumption. In the revised manuscript we will insert a dedicated lemma providing this explicit estimate, confirming that the resolvent bounds are inherited uniformly and that the contour-integral argument preserves the full order without mesh-dependent degradation. revision: yes

Circularity Check

0 steps flagged

No circularity: convergence analysis adapts external Lubich-Ostermann techniques to elliptic-parabolic setting under stated weak coupling condition.

full rationale

The derivation encodes the Runge-Kutta discretization in an operator-valued generating function, applies Parseval's identity on the unit circle, and adapts Fourier stability plus perturbation arguments from the 1995 Lubich-Ostermann reference (external authors). The weak coupling condition is introduced to justify the adaptation for the Schur complement arising from elliptic-parabolic coupling, but this does not reduce any central claim to a self-definition, a fitted input renamed as prediction, or a self-citation chain. The spectral decomposition for the iterative schemes is likewise presented as an independent tool. No quoted step equates the target convergence order to its own inputs by construction; the analysis is self-contained against the external benchmark once the weak coupling assumption is granted.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The work rests on standard assumptions from numerical analysis of PDEs and the cited 1995 stability framework; no free parameters or new entities are introduced in the abstract description.

axioms (1)
  • domain assumption Fourier stability and perturbation techniques from Lubich and Ostermann (1995) adapt to the present elliptic-parabolic decoupling setting.
    Central to the convergence proof for semi-explicit schemes as described in the abstract.

pith-pipeline@v0.9.0 · 5706 in / 1271 out tokens · 44088 ms · 2026-05-22T03:25:23.877964+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

24 extracted references · 24 canonical work pages

  1. [1]

    Altmann and M

    R. Altmann and M. Deiml. A novel iterative time integration scheme for linear poroelasticity. Electron. Trans. Numer. Anal. , 60:256--275, 2024

    work page 2024

  2. [2]

    Altmann and M

    R. Altmann and M. Deiml. A second-order iterative time integration scheme for linear poroelasticity. SIAM J. Sci. Comput. , 47(4):B875--B898, 2025

    work page 2025

  3. [3]

    Altmann, R

    R. Altmann, R. Maier, and B. Unger. Semi-explicit discretization schemes for weakly-coupled elliptic-parabolic problems. Math. Comp. , 90:1089--1118, 2021

    work page 2021

  4. [4]

    Altmann, R

    R. Altmann, R. Maier, and B. Unger. Semi-explicit integration of second order for weakly coupled poroelasticity. BIT Numer. Math. , 64:20, 2024

    work page 2024

  5. [5]

    Altmann, A

    R. Altmann, A. Mujahid, and B. Unger. Higher-order iterative decoupling for poroelasticity. Adv. Comput. Math. , 50:11, 2024

    work page 2024

  6. [6]

    Altmann, A

    R. Altmann, A. Mujahid, and B. Unger. Decoupling multistep schemes for elliptic--parabolic problems. SMAI J. Comput. Math. , 2026. accepted for publication

    work page 2026

  7. [7]

    M. A. Biot. General theory of three-dimensional consolidation. J. Appl. Phys , 12(2):155--164, 1941

    work page 1941

  8. [8]

    H. Brezis. Functional analysis, Sobolev spaces and partial differential equations . Universitext. Springer, New York, London, 2011

    work page 2011

  9. [9]

    Ern and S

    A. Ern and S. Meunier. A posteriori error analysis of E uler- G alerkin approximations to coupled elliptic-parabolic problems. ESAIM: Math. Model. Numer. Anal. , 43(2):353--375, 2009

    work page 2009

  10. [10]

    Eliseussen, M

    E. Eliseussen, M. E. Rognes, and T. B. Thompson. A posteriori error estimation and adaptivity for multiple-network poroelasticity. ESAIM : Math. Model. Numer. Anal. , 57(4):1921--1952, 2023

    work page 1921

  11. [11]

    Hairer and G

    E. Hairer and G. Wanner. Solving Ordinary Differential Equations II : Stiff and Differential--Algebraic Problems . Springer Berlin Heidelberg, Berlin, Heidelberg, 1996

    work page 1996

  12. [12]

    Kunkel and V

    P. Kunkel and V. Mehrmann. Differential-Algebraic Equations. Analysis and Numerical Solution . European Mathematical Society, Z \"u rich, 2006

    work page 2006

  13. [13]

    Kim, H.A

    J. Kim, H.A. Tchelepi, and R. Juanes. Stability and convergence of sequential methods for coupled flow and geomechanics: Drained and undrained splits. Comput. Meth. Appl. Mech. Eng. , 200(23):2094--2116, 2011

    work page 2094

  14. [14]

    Kim, H.A

    J. Kim, H.A. Tchelepi, and R. Juanes. Stability and convergence of sequential methods for coupled flow and geomechanics: Fixed-stress and fixed-strain splits. Comput. Meth. Appl. Mech. Eng. , 200(13):1591--1606, 2011

    work page 2011

  15. [15]

    Lubich and A

    C. Lubich and A. Ostermann. Runge--Kutta methods for parabolic equations and convolution quadrature. Math. Comp. , 60(201):105--131, 1993

    work page 1993

  16. [16]

    Lubich and A

    C. Lubich and A. Ostermann. Runge--Kutta approximation of quasi-linear parabolic equations. Math. Comp. , 64(210):601--627, 1995

    work page 1995

  17. [17]

    C. Lubich. Convolution quadrature and discretized operational calculus. I . Numer. Math. , 52(2):129--145, 1988

    work page 1988

  18. [18]

    C. Lubich. Convolution quadrature and discretized operational calculus. II . Numer. Math. , 52(4):413--425, 1988

    work page 1988

  19. [19]

    Mikeli\' c and M

    A. Mikeli\' c and M. F. Wheeler. Convergence of iterative coupling for coupled flow and geomechanics. Comput. Geosci. , 17(3):455--461, 2013

    work page 2013

  20. [20]

    Nevanlinna and F

    O. Nevanlinna and F. Odeh. Multiplier techniques for linear multistep methods. Numer. Funct. Anal. Optim. , 3(4):377--423, 1981

    work page 1981

  21. [21]

    R. E. Showalter. Diffusion in poro-elastic media. J. Math. Anal. Appl. , 251(1):310--340, 2000

    work page 2000

  22. [22]

    J. C. Vardakis, D. Chou, B. J. Tully, C. C. Hung, T. H. Lee, P. H. Tsui, and Y. Ventikos. Investigating cerebral oedema using poroelasticity. Med. Eng. Phys. , 38(1):48--57, 2016

    work page 2016

  23. [23]

    E. Zeidler. Nonlinear functional analysis and its applications. 2A : Linear monotone operators . Springer, New York, Berlin, Heidelberg, 1990

    work page 1990

  24. [24]

    M. D. Zoback. Reservoir geomechanics . Cambridge University Press, Cambridge, 2010

    work page 2010