---
title: 'LinOSS-IMEX: Implicit-Explicit Time Integration'
url: https://www.emergentmind.com/topics/implicit-explicit-time-integration-linoss-imex
type: topic
---

# LinOSS-IMEX: Implicit-Explicit Time Integration

Implicit-Explicit Time Integration (LinOSS-IMEX) schemes are a class of time discretization strategies for stiff systems of partial differential equations (PDEs) or ordinary differential equations (ODEs), particularly suited to problems exhibiting multiple spatial or temporal scales, such as nonlinear parabolic PDEs in finance, multiscale fluid–structure interactions, and advection-diffusion–reaction equations. LinOSS-IMEX schemes combine operator splitting of linear/nonlinear or stiff/nonstiff terms (the “LinOSS” or linear-operator splitting paradigm) with advanced implicit-explicit integration techniques, generally of Runge–Kutta or general linear method (GLM) type, to maximize both stability and computational efficiency. The methodology is characterized by selective, stage-wise implicit treatment of stiff (often diffusive or acoustic) linear components and explicit, high-order accurate, nonstiff handling of advective, nonlinear, or source terms. 

## 1. Formulation and Operator Splitting Principles

The LinOSS-IMEX framework begins with the identification and decomposition of a semi-discrete system derived from a spatial discretization (e.g., finite volume, finite element, spectral) of a PDE. The canonical structure is
\[
\frac{dU}{dt} = F(U) + S(U)
\]
where \(F(U)\) denotes the non-stiff (explicit) part—commonly advective, nonlinear, or reactive—and \(S(U)\) the stiff (implicit) part, typically linear diffusion or fast wave propagation [2409.01125].

A representative example is the one-dimensional parabolic PDE (with possible nonlinearities) in option pricing:
\[
\frac{\partial u}{\partial t} + a(x,t)\frac{\partial u}{\partial x} + b(x,t)\frac{\partial^2 u}{\partial x^2} + H(x,t,u,\partial_x u) = 0
\]
which is recast in conservative form, discretized in space by a finite-volume scheme (explicit advection, implicit diffusion), leading to the split ODE system above [2409.01125, 2409.01131].

IMEX time integrators then advance \(U\) by treating \(F(U)\) explicitly (Runge–Kutta, multistep, or GLM stages) and \(S(U)\) implicitly (often via diagonally implicit RK, BDF, or related solvers) [1307.8303, 1302.2689]. This splitting is fundamental to achieving enlarged stability regions (CFL relaxation for stiff terms) without the complexity or computational cost of fully implicit solves for the entire system.

## 2. IMEX Time Integration Schemes

The LinOSS-IMEX approach employs time integrators defined by a pair of Butcher tableaus—one explicit for non-stiff operators, one implicit for stiff operators—constructed to reach at least second order globally. For example, the L-stable, second order IMEX-SSP2(2,2,2) scheme is characterized by

**Explicit (ERK) tableau:**
\[
\begin{array}{c|cc}
0 & 0 & 0 \\
1 & 1 & 0 \\
\hline
  & \tfrac12 & \tfrac12
\end{array}
\]

**Implicit (DIRK) tableau:**
\[
\begin{array}{c|cc}
\gamma & \gamma & 0 \\
1-\gamma & 1-2\gamma & \gamma \\
\hline
   & \tfrac12 & \tfrac12
\end{array}
\qquad \gamma = 1-\frac{1}{\sqrt{2}}
\]

with stages given, for each time interval \([t^n, t^{n+1}]\), by a sequence of implicit and explicit updates [2409.01125, 2409.01131]. The IMEX methodology generalizes to higher orders via multi-stage ARS (Ascher-Ruuth-Spiteri) [2202.05734], BPR, or GLM-based schemes [1302.2689].

## 3. Stage Equations, Stability, and CFL Conditions

For a \(s\)-stage LinOSS-IMEX-RK scheme:
- Non-stiff (explicit) terms are advanced via forward-explicit (Heun-type, RK2, or higher order) updates;
- Stiff (implicit) terms are advanced by inverting diagonal or lower-triangular systems per stage (often only diffusion or linear fast waves).

Typical update for two-stage IMEX-SSP2(2,2,2):
\[
\begin{aligned}
U^{(1)} &= U^n + \gamma\,\Delta t\,S(U^{(1)}) \\
U^{(2)} &= U^n -\Delta t\,F(U^{(1)}) + (1-2\gamma)\,\Delta t\,S(U^{(1)}) + \gamma\,\Delta t\,S(U^{(2)}) \\
U^{n+1} &= U^n - \frac{\Delta t}{2} \big[F(U^{(1)}) + F(U^{(2)})\big] + \frac{\Delta t}{2} \big[S(U^{(1)}) + S(U^{(2)})\big]
\end{aligned}
\]
which is provably second-order in time [2409.01125, 2409.01131].

LinOSS-IMEX schemes only require the explicit part (advection/convection) to satisfy a hyperbolic-type CFL
\[
\alpha\,\frac{\Delta t}{\Delta x} \le \mathrm{CFL}_{\rm adv}\sim0.5
\]
with \(\alpha\) the max explicit advection wave speed. The parabolic (diffusion) CFL of explicit schemes, \(\Delta t \leq C (\Delta x)^2\), is absent due to the implicit treatment of \(S\) [2409.01125]. The implicit tableau is chosen L-stable to encompass the stiff eigenvalues. This structure extends to 2D and systems with mixed derivatives [2409.01131], as well as to various application domains.

## 4. Extensions, Order Reduction, and High-Order Strategies

LinOSS-IMEX supports:
- High-order in time via multi-stage IMEX-RK schemes (e.g., ARS343 for order 3, KC664 for order 4) [2202.05734].
- High-order in space by replacing FV linear reconstructions with CWENO, DG, or compact-difference stencils [1912.10023]. 

A noted phenomenon is order reduction of high-order schemes in problems involving DAEs or extreme stiffness (large diffusion, small relaxation), notably in higher-order IMEX-RKs in Boussinesq convection and Goldstein–Taylor systems [2202.05734, 1307.8303]. For problems with DAEs, globally stiffly accurate (GSA) or special DAE-aware IMEX schemes are preferred to avoid loss of convergence order [2202.05734].

IMEX general linear methods (IMEX-GLMs), particularly DIMSIM-type with stage order \(q=p\), avoid the mixed-term order reduction of classical IMEX-RK, retaining full convergence order even for stiff–nonstiff partitioning [1302.2689].

## 5. Applications and Numerical Performance

LinOSS-IMEX schemes are broadly deployed in:
- Nonlinear parabolic PDEs in option pricing, enabling accurate solution of price and Greeks with non-regular terminal payoffs, while allowing parabolic time steps proportional to \(\Delta x\) rather than \((\Delta x)^2\) [2409.01125, 2409.01131].
- High-Re, stiff-diffusion or DAE systems, such as pseudo-spectral Boussinesq convection, where third- and fourth-order IMEX-RK schemes (ARS343, KC664) outperform classical CNAB2 integrators in both speed and accuracy at their maximum stable \(\Delta t\) [2202.05734].
- ADR equations, where the IMEX midpoint/RK2 with compact-difference stencils achieves large regions of neutral stability, minimal phase error, and suppression of spurious q-waves [1912.10023].
- Coupled free flow–porous media models with nonlinear interface (Lions condition), exploiting SAV-based IMEX schemes for provable unconditional stability and linear, decoupled algebraic solves [2405.11223].
- Multiscale hyperbolic systems, where LinOSS-IMEX (often with TVD or MOOD mechanisms) enables scale-independent CFL, non-oscillatory solutions, and robust detection/limiting [2501.03994].
- Industrial compressible flow simulations, using cell-wise explicit–implicit hybridization (AION/LinOSS-IMEX) to blend Heun and Crank–Nicolson advances optimally in space [1903.12568].

In all applications, LinOSS-IMEX methods yield substantial gains in allowable time steps, CPU time, and solution regularity under stiff or heterogeneous conditions.

## 6. Implementation and Algorithmic Details

Efficient LinOSS-IMEX integrators involve:
- Assembly of explicit and implicit operator splits per time stage;
- Piecewise linear or higher-order reconstruction, often with non-oscillatory TVD or WENO limiters for spatial derivatives;
- For each time step, solution of two (2nd order) or more (higher order) implicit diffusion (or stiff linear) algebraic systems, typically sparse and efficiently solvable by factorization or iterative methods;
- Weighted combination of explicit/implicit stage outputs to update the solution;
- Implicit enforcement of boundary conditions in the solution of the diffusion or pressure operator;

The methods are compatible with a variety of ODE solvers (Newton/fixed-point for nonlinear diffusion, direct or Krylov solvers for linear systems) and can exploit precomputed matrix factorizations for constant-coefficient implicit operators [2409.01125, 2409.01131, 1307.8303].

For problems with severe scale separation (e.g., low Mach flows), LinOSS-IMEX–MOOD defers to stricter TVD variants in troubled cells, reverting to first-order to guarantee non-oscillatory behavior, while maintaining high order in smooth regions [2501.03994].

## 7. Theoretical Properties and Numerical Evidence

LinOSS-IMEX schemes possess:
- Rigorous second-order (or higher) global accuracy, with error estimates in \(L^1\), \(L^\infty\), and energy-based norms under broad regularity assumptions [2409.01125, 2409.01131, 2405.11223].
- Unconditional (often CFL-free) stability with respect to stiff linear terms, owing to L-stability of the implicit tableau [2409.01125, 2405.11223].
- Uniform convergence and asymptotic preservation in singular-perturbation regimes, e.g., in Goldstein–Taylor diffusion, where the discrete scheme preserves the correct limit as \(\varepsilon \rightarrow 0\) [1307.8303].
- Demonstrated computational superiority over fully explicit or fully implicit schemes in both benchmark and real-world test cases; performance gains up to \(10^4\times\) are reported in large-N, diffusion-dominated settings [2409.01125, 2409.01131, 2202.05734].
- Robustness to non-regular initial or boundary data, and preservation of monotonicity or TVD properties with appropriate limiting [2409.01131, 2405.11223, 2501.03994].

A plausible implication, given the breadth of tested problems, is that LinOSS-IMEX methodology is broadly transferrable to other stiff systems, provided proper operator splitting and tableau selection.

---

**References**  
- “IMEX-RK finite volume methods for nonlinear 1d parabolic PDEs. Application to option pricing” [2409.01125]  
- “Second order finite volume IMEX Runge-Kutta schemes for two dimensional parabolic PDEs in finance” [2409.01131]  
- “An assessment of implicit-explicit time integrators for the pseudo-spectral approximation of Boussinesq thermal convection in an annulus” [2202.05734]  
- “Implicit-explicit-compact methods for advection diffusion reaction equations” [1912.10023]  
- “Asymptotic Preserving time-discretization of optimal control problems for the Goldstein-Taylor model” [1307.8303]  
- “A class of new linear, efficient and high-order implicit-explicit methods for the coupled free flow-porous media system based on nonlinear Lions interface condition” [2405.11223]  
- “TVD-MOOD schemes based on implicit-explicit time integration” [2501.03994]  
- “A coupled implicit-explicit time integration method for compressible unsteady flows” [1903.12568]  
- “Partitioned and implicit-explicit general linear methods for ordinary differential equations” [1302.2689]

Source: https://www.emergentmind.com/topics/implicit-explicit-time-integration-linoss-imex