---
title: Global-in-Time Iterative Decoupled Methods
url: https://www.emergentmind.com/topics/global-in-time-iterative-decoupled-algorithms
type: topic
---

# Global-in-Time Iterative Decoupled Methods

Global-in-time iterative decoupled algorithms are numerical methods for evolutionary and multiphysics problems in which the object of iteration is the solution over an entire time interval, or over large time blocks, while the constituent subphysics, subdomains, or operator blocks are solved separately inside each iteration. In the strict sense, the iteration acts on time trajectories, space–time interface variables, or all-at-once space–time algebraic systems; in a weaker but common sense, the algorithm is local in time yet admits stability or error estimates that are uniform in time. This distinction is explicit in recent work on Biot and MPET systems, nonlinear Stokes–Darcy coupling, and time-parallel parabolic solvers [2409.18391] [2310.15457] [2007.01736] [1802.08126].

## 1. Conceptual scope and defining distinctions

In strict usage, “global-in-time” means that the coupled problem is posed and iterated on the full interval \((0,T)\) at once, rather than advanced sequentially with localized interface or coupling iterations. This is the formulation adopted for the nonlinear Stokes–Darcy interface problem, where the interface unknown \(\lambda\) is defined on \(\Gamma\times(0,T)\) and the interface equation is solved over the whole time interval by a nested Newton–Krylov method [2007.01736]. The same terminology is used for Crank–Nicolson- and BDF-based Biot algorithms in which one outer fixed-stress iteration updates the full temporal history \(\{(\mathbf u^n,\xi^n,p^n)\}_{n=0}^N\), with a forward flow sweep and a mechanics step that is independent across time levels [2409.18391] [2508.04082].

A weaker usage appears in time-stepping iterative decouplings. In the total-pressure MPET algorithm, the splitting is performed at each time step \(t^n\), and the inner iteration converges to the monolithic solution at that time level; nevertheless, the stability and error bounds are uniform in time, and the contraction factor is independent of \(\Delta t\) [2310.15457]. The three-field Biot algorithm based on total pressure has the same character: the iteration is local in time, but its convergence requires no restriction on \(\Delta t\) or \(h\), and the coupled backward Euler scheme satisfies discrete energy estimates that are uniform over the entire interval [2204.08403]. This suggests that the phrase “global-in-time” in the literature refers both to genuine space–time iteration and to time-marching schemes with uniform-in-time control, and those two meanings should not be conflated.

A related but distinct class consists of all-at-once time-parallel solvers. For implicit Euler discretizations of linear parabolic problems, the standard nonsymmetric time-global system \(\mathbf B\mathbf u=\mathbf f\) is reformulated as a symmetric saddle-point system, and the resulting inexact Uzawa iteration acts on all time levels simultaneously [1802.08126]. For block \(\alpha\)-circulant preconditioned all-at-once schemes, the global error iteration is analyzed directly in terms of the stability function of the underlying one-step or multistep integrator [2102.04646].

## 2. Reformulation patterns and operator structure

A recurring structural device is the introduction of a total-pressure variable. In Biot’s three-field formulation,
\[
\xi = -\lambda\,\nabla\cdot\mathbf u + \alpha\,p,
\]
while in the quasi-static MPET system,
\[
\xi = \vec{\alpha}^{\mathsf T}\vec p - \lambda \operatorname{div}\mathbf u.
\]
These substitutions transform the original displacement–pressure formulation into a generalized Stokes problem for \((\mathbf u,\xi)\) coupled to a parabolic problem for the pressure variables [2204.08403] [2310.15457]. The immediate algorithmic consequence is a canonical two-block split: a Stokes-type saddle-point solve for mechanics and a diffusion–reaction solve for transport.

The same separation appears in the Crank–Nicolson and BDF-based global-in-time Biot schemes. In the global-in-time Crank–Nicolson algorithm, the flow step computes the entire sequence \(p_h^{n,i}\) using the previous iterate of the total pressure differences, and the mechanics step then solves \((\mathbf u_h^{n,i},\xi_h^{n,i})\) independently for each \(n\) [2409.18391]. In the BDF-2 full-order global-in-time algorithm, the diffusion subproblem is
\[
a_3 ( \mathbb{D}_{t}^2 p_{h,f}^{n,i} , q_h ) + d ( p_{h,f}^{n,i},q_h )  =  ( g^{n} , q_h ) + c( q_h , \mathbb{D}_{t}^2 \xi_{h,f}^{n,i-1} ),
\]
followed by a generalized Stokes problem for \((\bm u_{h,f}^{n,i},\xi_{h,f}^{n,i})\) at all time levels [2508.04082].

A second major reformulation pattern is interface reduction. In nonlinear Stokes–Darcy coupling, a Lagrange multiplier
\[
\lambda = -n_f\cdot(\nu_f(|D(u_f)|)D(u_f)-p_f I)\cdot n_f = p_p
\]
represents the interface normal stress. This yields space–time Steklov–Poincaré operators \(S_f(\lambda)\) and \(S_p(\lambda)\), and the full coupled evolution is reduced to the interface equation \(S_f(\lambda)+S_p(\lambda)=0\) on \((0,T)\times\Gamma\) [2007.01736]. Closely related interface formulations appear in recent fluid–structure domain decomposition methods with local time stepping, where the coupled system is again posed as a space–time interface problem [2408.02746].

A third pattern is all-at-once operator symmetrization. For implicit Euler discretizations of linear parabolic equations, the nonsymmetric time-global matrix \(\mathbf B=\mathbf K+\mathbf A\) is left-preconditioned by \(\mathbf P=\mathbf A^{-1}\mathbf K+\mathbf I\), producing the symmetric positive definite operator
\[
\mathbf S=\mathbf P^{\top}\mathbf B
= \mathbf{K}^\top \mathbf{A}^{-1}\mathbf{K} + \mathbf{K} + \mathbf{K}^\top + \mathbf{A},
\]
which is then embedded in a symmetric saddle-point system [1802.08126]. This algebraic recasting is the basis for robust global-in-time preconditioning.

## 3. Principal algorithmic families

The contemporary literature contains several distinct but structurally related families.

| Family | Representative structure | Decoupling mechanism |
|---|---|---|
| Total-pressure poromechanics | generalized Stokes block + parabolic block | alternate transport and mechanics solves |
| Space–time interface methods | interface equation on \((0,T)\times\Gamma\) | independent subdomain solves with interface updates |
| All-at-once time-parallel solvers | global space–time linear system | Fourier/DST or circulant diagonalization across time |
| DG correction iterations | DG weak formulation on each time slab | preconditioned Picard or correction equations |

In total-pressure poromechanics, the split may be time-stepping or genuinely global-in-time. The MPET method alternates a symmetric positive definite pressure solve and a generalized Stokes solve at each time level [2310.15457]. The Biot Crank–Nicolson and BDF-2 methods extend this idea to strict global-in-time fixed-stress iterations, with the mechanics stage partially parallel-in-time [2409.18391] [2508.04082]. The unified thermo-/multiple-network poroelasticity model further generalizes the pattern to displacement, two generalized pressures, and total pressure, with a transport block solved over the full interval and a mechanics block decoupled from it at each iteration [2603.29275].

In interface methods, the nonlinear Stokes–Darcy problem is solved by a nested Newton–Krylov procedure in which each Krylov matrix–vector product requires one linearized Stokes and one linearized Darcy solve over the full time interval [2007.01736]. In global fluid–structure domain decomposition, Steklov–Poincaré and Robin formulations again shift the iteration to interface traces on space–time slabs [2408.02746].

In all-at-once time-parallel methods, the emphasis is on linear algebra. The inf-sup-based saddle-point reformulation for parabolic equations uses an inexact Uzawa iteration with a time-parallel preconditioner built from a Discrete Sine Transform [1802.08126]. The block \(\alpha\)-circulant class instead diagonalizes an \(\alpha\)-circulant temporal preconditioner by FFT, so each iteration reduces to independent solves for the time modes [2102.04646].

Finally, DG-based iterative time integration starts from the weak DG formulation on each interval and applies a preconditioned Picard iteration. This yields explicit, implicit, and semi-implicit SDG schemes with arbitrary order, and a multilevel extension intended for coupling with space–time multigrid and parallel-in-time methods [1610.01324].

## 4. Convergence theory, contraction mechanisms, and uniform-in-time bounds

The sharpest results concern contraction in norms built on time derivatives of the coupling variable. For the Crank–Nicolson global-in-time Biot iteration,
\[
\sum_{n=1}^N \|e_{\xi,g}^{n,i}-e_{\xi,g}^{n-1,i}\|_{L^2}^2
\le L_g\sum_{n=1}^N \|e_{\xi,g}^{n,i-1}-e_{\xi,g}^{n-1,i-1}\|_{L^2}^2,
\]
with
\[
L_g = \frac{1}{\big(\frac{c_0\lambda}{\alpha^2}+1\big)\big(C(\beta,\mu)\lambda+1\big)}<1,
\]
and \(L_g\) is independent of \(\Delta t,h,N\) [2409.18391]. The BDF-2 global-in-time full-order model satisfies the analogous contraction
\[
\sum_{n=2}^N \| \mathbb{D}_{t}^2 e_{\xi,f}^{n,i} \|_{L^2}^2
\le K^2 \sum_{n=2}^N \| \mathbb{D}_{t}^2 e_{\xi,f}^{n,i-1} \|_{L^2}^2,
\]
with
\[
K = \frac{1}{\big( \frac{c_0 \lambda }{\alpha^2} + 1\big) \big( C(\tilde{\beta}, \mu)\lambda + 1 \big)} < 1
\]
[2508.04082].

The thermo-/multiple-network poroelasticity scheme obtains a continuous contraction
\[
\|\partial_t e_{\xi}^{i}\|_{L^2(0,T_f; W)} \leq L_{\mathrm{con}}\,
\|\partial_t e_{\xi}^{i-1}\|_{L^2(0,T_f; W)},
\qquad
L_{\mathrm{con}} = \frac{1}{\sqrt{1+2\lambda C_\mu}}<1,
\]
and a fully discrete contraction
\[
\bigg( \Delta t \sum_{n=1}^N \| \bar{\partial}_t e_\xi^{n,i} \|^2 \bigg)^{1/2}
\le L_{\mathrm{dis}}
\bigg( \Delta t \sum_{n=1}^N \| \bar{\partial}_t e_\xi^{n,i-1} \|^2 \bigg)^{1/2},
\]
with
\[
L_{\mathrm{dis}} = \frac{1}{\sqrt{(1+2\lambda C_{\tilde\mu})(1+2\lambda h^2 C_{\mathrm{stab}})}}<1
\]
[2603.29275]. Here the transport stabilizations enter the contraction factor explicitly.

In local-in-time total-pressure splittings, the iteration is still unconditional with respect to \(\Delta t\). For MPET, if \(\delta=\min_i c_i>0\), then
\[
\|e_{\xi}^k\| \le C^* \|e_{\xi}^{k-1}\|,
\qquad
C^* = \frac{\|\vec{\alpha}\|^2/\lambda}{\delta + \|\vec{\alpha}\|^2/\lambda}<1,
\]
while convergence for \(\min_i c_i=0\) is proved by contradiction [2310.15457]. For the three-field Biot iteration, the theorem gives
\[
\|e_\xi^i\|_{L^2(\Omega)} \le C\,\|e_\xi^{i-1}\|_{L^2(\Omega)},
\]
with
\[
C = \frac{(\alpha^2/\lambda)^2}{\big(c_0+\alpha^2/\lambda\big)^2}\le 1,
\]
and convergence holds even for \(c_0=0\) [2204.08403].

In all-at-once parallel-in-time methods, the contraction constants are tied to time-integrator stability rather than PDE parameters. For block \(\alpha\)-circulant iterations, stable one-step methods satisfy
\[
\|(I_t\otimes P)\,\mathbf{err}^{k+1}\|_\infty
\le \frac{\alpha}{1-\alpha}\,\|(I_t\otimes P)\,\mathbf{err}^{k}\|_\infty,
\]
and stable linear multistep methods satisfy
\[
\|(I_t\otimes P)\,\mathbf{err}^{k+1}\|_\infty
\le \frac{c\alpha}{1-c\alpha}\,\|(I_t\otimes P)\,\mathbf{err}^{k}\|_\infty,
\]
with bounds independent of \(\Delta t\) and \(\sigma(A)\) [2102.04646]. For implicit-Euler all-at-once parabolic solvers, spectral bounds for the DST preconditioner yield Uzawa convergence rates independent of the number of time steps, final time, and spatial mesh sizes [1802.08126].

Not all strict global-in-time algorithms have complete proofs. For the nonlinear Stokes–Darcy space–time interface method, there is no full analytic convergence theory for the nonlinear global-in-time interface algorithm; the work instead relies on well-posedness and numerical evidence [2007.01736].

## 5. Parallelism, local time stepping, and model reduction

A principal motivation for strict global-in-time decoupling is parallelism in the time direction. In the global-in-time Crank–Nicolson Biot algorithm, the mechanics solves for \(n=1,\dots,N\) are independent and in parallel, while the flow solve remains a forward sweep; the paper therefore characterizes the method as having a partially parallel-in-time feature [2409.18391]. The same work gives the speedup estimate
\[
\text{Speedup} = \frac{T_p+T_{(\mathbf u,\xi)}}{T_p+2^{-m}T_{(\mathbf u,\xi)}},
\]
when \(2^m\) processors are used for the mechanics stage [2409.18391].

Global-in-time interface methods also support asynchronous temporal resolution. In nonlinear Stokes–Darcy coupling, different time grids \(\mathcal T_f\) and \(\mathcal T_p\) are coupled by \(L^2\) projections \(\Pi_{p,f}\) and \(\Pi_{f,p}\), so nonmatching time grids are handled without introducing a separate mortar time grid [2007.01736]. The fluid–structure domain decomposition methods are explicitly designed for local time stepping in the subsystems [2408.02746].

Reduced-order acceleration has been incorporated directly into global-in-time decouplings. The POD-based BDF-2 Biot algorithm keeps the diffusion step at full order but replaces most generalized Stokes solves by solves in iteration-dependent POD spaces built from snapshot times \(\bm{\Lambda}\) [2508.04082]. Theoretical estimates show convergence to a neighborhood controlled by the ROM error,
\[
\lim_{i \rightarrow \infty}
\sum_{n=2}^N \| \mathbb{D}_{t}^2 e_{\xi,r}^{n,i} \|_{L^2}^2
\leq C_2 \epsilon_{ROM}(\bm{\Lambda}, N_r),
\]
and the heterogeneous injection–production test reports about \(3\) minutes/iteration for the FOM, about \(6\) minutes/iteration for ROM with the full index set, and about \(1\) minute/iteration for a sparse Legendre-based snapshot set, achieving \(\approx 3\times\) speed-up while maintaining high accuracy [2508.04082].

At the linear-algebraic end, the DST-preconditioned all-at-once parabolic solver has theoretical parallel complexity \(O(\log N)\) in the number of time steps, and the cost per iteration is
\[
O\big((\log N+1)c_v + c_m + c_s\big)
\]
[1802.08126].

## 6. Applications, limitations, and recurrent misunderstandings

The application range is broad. Total-pressure MPET decoupling has been used for brain flow simulations [2310.15457]. Global-in-time Biot algorithms have been tested on manufactured solutions, Barry–Mercer, and Mandel’s problem [2409.18391]. The unified thermo-/multiple-network formulation covers both linear thermo-poroelasticity and Barenblatt–Biot dual-network models [2603.29275]. Space–time interface methods address nonlinear Stokes–Darcy coupling with nonconforming time grids [2007.01736]. Global iterative solvers have also been studied for the Kadanoff–Baym equations in nonequilibrium DMFT, where the unknowns are two-time Green’s functions rather than PDE fields [2512.11371].

Several misconceptions recur in this literature. First, not every decoupled iterative algorithm is global-in-time in the strict sense: the MPET and three-field Biot total-pressure schemes iterate per time step, even though their stability and error estimates are uniform in time [2310.15457] [2204.08403]. Second, global-in-time does not automatically mean fully parallel-in-time: in Biot fixed-stress schemes, the flow sweep remains sequential and only the mechanics stage is independent across time levels [2409.18391]. Third, global-in-time formulations do not automatically outperform time-stepping. For the Kadanoff–Baym equations, several stable global solvers exhibit an iteration count that scales roughly linearly with the number of time steps, and with dense matrix–matrix convolutions the total work scales as \(O(N_t^4)\), which is asymptotically worse than standard time-stepping [2512.11371]. The same work identifies a propagating residual front as the mechanism behind that scaling, and shows that a standard forward fixed point iteration does not converge stably at long propagation times [2512.11371].

A final limitation concerns parameter robustness. Total-pressure formulations control locking and enable stable Stokes pairs, but the iteration rate need not be parameter-robust in the strongest sense. In MPET, the contraction factor depends on \(\lambda\), \(\vec\alpha\), and the storage coefficients; the paper explicitly does not claim a fully parameter-robust rate independent of those physical parameters [2310.15457].

Taken together, these developments show that the modern theory of global-in-time iterative decoupled algorithms is organized around a small set of structural ideas: reformulation by total pressure or interface variables, iteration in space–time norms, decoupling into canonical subproblems, and preconditioning or stabilization strong enough to preserve uniform-in-time control. The main unresolved questions concern nonlinear convergence theory in the strict global-in-time setting, the interaction between compression and iteration in genuinely large space–time systems, and the extent to which time-parallel efficiency can be obtained without sacrificing the robust stability properties that have made these methods technically attractive [2007.01736] [2512.11371].

Source: https://www.emergentmind.com/topics/global-in-time-iterative-decoupled-algorithms