---
title: Stabilized IMEX Crank–Nicolson Method
url: https://www.emergentmind.com/topics/stabilized-imex-crank-nicolson-method
type: topic
---

# Stabilized IMEX Crank–Nicolson Method

The stabilized implicit–explicit (IMEX) Crank–Nicolson method is a class of time-stepping schemes for parabolic, convection–diffusion, incompressible flow, and phase-field problems, designed to achieve second-order temporal accuracy, unconditional or robust stability, and efficient linear algebra at each step. The essential idea is to split the governing equation into linear (stiff) and nonlinear or advective (non-stiff) parts, treating the former implicitly via Crank–Nicolson averaging, and the latter explicitly, often with extrapolation and stabilization. To maintain stability and structure-preservation under severe time step or spatial mesh constraints, targeted stabilization mechanisms such as continuous interior penalty (CIP) terms, biharmonic modification, or doubly-stabilized splittings are incorporated. These methods have seen extensive analysis and application, notably for high-Péclet or high-Reynolds convection–diffusion, Oseen/Navier–Stokes flows, Cahn–Hilliard, and Allen–Cahn equations [2012.05727][1710.03604][2405.19483][2004.05163][2405.12562][2310.19663].

## 1. Abstract Formulation and Model Problems

The stabilized IMEX Crank–Nicolson method applies to evolution equations of the abstract form
$$
\partial_t u = F(u) = \mathcal{L}(u) + \mathcal{N}(u),
$$
where $\mathcal{L}$ is a linear (possibly stiff, such as diffusive or biharmonic) component and $\mathcal{N}$ is a nonlinear, advective, or otherwise explicitly treated part.

Typical model equations include:

- **Convection–Diffusion**: $\partial_t u + \beta\cdot\nabla u - \mu \Delta u = f$ [2012.05727]
- **Oseen/Navier–Stokes**: $\partial_t u + \beta\cdot\nabla u + \nabla p - \nu \Delta u = f,\quad\nabla\cdot u=0$ [2405.12562]
- **Cahn–Hilliard and Thin-Film**: $\partial_t u = \nabla\cdot[M(u)\nabla(W'(u) - \Delta u)]$ [2405.19483]
- **Allen–Cahn with Mobility**: $\partial_t \phi = -M(\phi) (-\varepsilon^2 \Delta \phi + F'(\phi))$ [2310.19663]

In each case, stabilization is introduced to control explicit terms and improve robustness at large time steps or high Reynolds/Péclet numbers.

## 2. Temporal Discretization Scheme

Stabilized IMEX Crank–Nicolson schemes combine implicit Crank–Nicolson averaging for stiff linear terms with explicit, often extrapolated or Adams–Bashforth-type, evaluation for nonlinear/convective components. Stabilizing corrections of $O(\tau^2)$ (where $\tau$ is the time step) are added to preserve energy dissipation, mass conservation, or maximum principles as required.

Generic two-level IMEX Crank–Nicolson scheme:
\[
(u^{n+1}-u^n, v) + \tau\,\mathcal{C}_{\text{stab}}(\hat{u}^{n+1}, v) + \tau\,\mathcal{A}\bigg(\frac{u^{n+1}+u^n}{2}, v\bigg) = \tau(f^{n+1}, v)
\]
where:
- $\mathcal{A}$ is the symmetric, implicit (e.g., diffusive) part,
- $\mathcal{C}_{\text{stab}}$ is the explicitly stabilized advective/nonlinear term, typically evaluated at an extrapolation $\hat{u}^{n+1}$ (for example, $\hat{u}^{n+1} = \tfrac{3}{2}u^n - \tfrac{1}{2}u^{n-1}$).

For problems with additional nonlinear bulk forces or degeneracies (as in Cahn–Hilliard or Allen–Cahn systems), further linear stabilization terms are introduced, often involving higher-order time differences and Laplacians [1710.03604][2310.19663][2004.05163].

## 3. Stabilization Techniques and Spatial Discretization

### Continuous Interior Penalty (CIP) and Jump Stabilization

For finite element methods, especially in convection–dominated regimes, symmetric stabilization terms such as continuous interior penalty are required for robustness. For example, the stabilized bilinear form for the convection term is:
\[
c_h(u_h, v_h) = \int_\Omega (\beta \cdot \nabla u_h)\,v_h\,dx + \gamma \sum_{F\in\mathcal{F}_{\rm int}} \int_F h_F^2 |\beta\cdot n_F| [\nabla u_h] \cdot [\nabla v_h] ds.
\]
Here, $\gamma>0$ is the stabilization parameter, and $[\nabla u_h]$ is the jump of the gradient across interior faces [2012.05727].

### Biharmonic and Doubly-Stabilized Splittings

For fourth-order problems (e.g., Cahn–Hilliard, thin-film), the splitting $\mathcal{L}(u) = -M_2 \Delta^2 u + M_1 \Delta u$ is adopted and further stabilization is enforced via biharmonic terms (with $M_2$ chosen relative to the maximum mobility) [2405.19483].

In phase-field models, doubly-stabilized IMEX CN schemes insert both first-difference Laplacian and second-difference (leapfrog-type) terms, denoted as $-A\tau \Delta (u^{n+1}-u^n) + B(u^{n+1} - 2u^n + u^{n-1})$, with $A$, $B$ chosen for unconditional stability and maximum bound principle (MBP) preservation [1710.03604][2310.19663].

## 4. CFL Conditions, Stability, and Error Analysis

Stability of the scheme is governed by coupled spatial (mesh) and temporal (time step) constraints, generally expressed through CFL-type numbers related to mesh Reynolds/Péclet:
\[
\begin{aligned}
&\mathrm{Co} = (\|\beta\|_\infty+1)\tau/h, \\
&\mathrm{Co}_{4/3} = \tau (\|\beta\|_\infty/h)^{4/3}.
\end{aligned}
\]
For $p$th-degree finite elements:
- Piecewise affine $(p=1)$: standard hyperbolic CFL suffices.
- Higher degree $(p\ge1)$: a stricter $4/3$-CFL, $\tau \le C h^{4/3}$, is required for stability [2012.05727][2405.12562].

Main stability and error results:
- **Discrete energy dissipation**: For appropriate stabilization parameters, modified discrete energies decrease monotonically for any $\tau>0$, i.e., unconditional stability [1710.03604][2004.05163][2310.19663][2405.19483].
- **Uniform bounds**: Solutions are bounded in natural norms (e.g., $L^2$, $H^1$, $L^\infty$ for MBP).
- **Error order**: $L^2$- and $L^\infty$-error estimates of $O(\tau^2 + h^{p+1/2})$ for convection–diffusion/Oseen; $O(\tau^2 + h^q)$ (with $q$ depending on spatial discretization) for phase-field models [2012.05727][2405.12562][1710.03604][2310.19663].

## 5. Algorithmic Details and Practical Implementation

At each time step, the stabilized IMEX Crank–Nicolson method necessitates solving systems with constant-coefficient, self-adjoint (possibly block) elliptic operators, as nonlinearity is handled explicitly. This enables the use of efficient direct or preconditioned iterative solvers, becoming particularly advantageous for spectral or uniform grid discretizations.

- In mixed formulations (Cahn–Hilliard), the coupled $(\phi^{n+1}, \mu^{n+1/2})$ system reduces to a block $2 \times 2$ system with time-invariant structure, admitting efficient block preconditioners [2004.05163].
- Adaptive time-stepping, based either on interface velocity or energy variation, further enhances efficiency during slow coarsening or equilibrium [2004.05163][2310.19663].

Stabilization parameters may be taken as minimal as possible in practice once empirical stability is observed, even though theoretical bounds may prescribe large values (e.g., $A = L^2/(16\varepsilon^2\gamma)$, $B = L/(2\varepsilon)$ for phase-field models, with $L$ a global Lipschitz constant) [1710.03604][2004.05163].

## 6. Applications, Numerical Experiments, and Extensions

The stabilized IMEX Crank–Nicolson method supports a broad spectrum of computational PDEs:

- **High Reynolds/Péclet flows**: Demonstrated with Oseen and Navier–Stokes equations for vortex dynamics, mixing layers, and wake flows at high Reynolds number, yielding robust error orders (e.g., $O(h^{k+1/2} + \tau^2)$ for $P_k$ elements) and stable energy dissipation [2405.12562].
- **Phase-field dynamics**: For Cahn–Hilliard and Allen–Cahn systems with degenerate mobility or sharp interface limits, the method ensures unconditional MBP, energy decay, and spectral/higher-order spatial accuracy, supporting adaptive time stepping for long-term evolution [1710.03604][2004.05163][2405.19483][2310.19663].
- **Thin-film equations**: Adaptation of the biharmonic stabilization yields energy-stable and mass-conserving solvers for dewetting, coarsening, and other fourth-order gradient flows [2405.19483].

Extensions include fractional step (“projection”) methods in incompressible flow, where the IMEX CN scheme underlies operator-splittings for pressure-velocity coupling without loss in order or structure preservation [2405.12562].

## 7. Theoretical Results and Parameter Selection

Central theorems and guiding results across model problems include:

- **Energy stability**: For stabilization parameters $(A,B)$ or $(M_1,M_2)$ satisfying explicit algebraic bounds (e.g., $M_2 \geq \alpha^* \max_x M(U^n(x))$ with $\alpha^*\simeq 0.97$ for IMEX2), unconditional discrete energy dissipation is guaranteed for any time step [1710.03604][2405.19483][2310.19663].
- **Error estimates**: For smooth solutions and appropriate time/space steps, global error bounds are established as $O(\tau^2 + h^{p+\frac12})$ (convection–diffusion), $O(\tau^2 + h^q)$ (phase-field), with prefactors depending algebraically (not exponentially) on small parameters such as $\varepsilon$ [2012.05727][1710.03604][2004.05163].
- **Maximum principle preservation**: For Allen–Cahn dynamics with variable mobility, the scheme preserves $|\Phi^n|\le1$ unconditionally (for sufficiently large stabilizers), with associated $L^\infty$-error estimates [2310.19663].

Tables summarizing recommended stabilization parameters and their empirical effects are reported in [2405.19483][2310.19663].

---

**References:**
- [2012.05727] E. Burman, J. Guzmán. "Implicit-explicit multistep formulations for finite element discretisations using continuous interior penalty"
- [2405.12562] E. Burman, D. Garg, J. Guzmán. "Implicit-explicit Crank-Nicolson scheme for Oseen's equation at high Reynolds number"
- [1710.03604] X. Wang, J. Yu. "Convergence Analysis of an Unconditionally Energy Stable Linear Crank-Nicolson Scheme for the Cahn-Hilliard Equation"
- [2004.05163] P. Cheng, F. Yang, H. Wu. "An Energy Stable Linear Diffusive Crank-Nicolson Scheme for the Cahn-Hilliard Gradient Flow"
- [2405.19483] M. Orizaga, T. Witelski. "IMEX methods for thin-film equations and Cahn-Hilliard equations with variable mobility"
- [2310.19663] X. Xu, P. Gao, Z. Qiao, T. Tang. "A linear doubly stabilized Crank-Nicolson scheme for the Allen-Cahn equation with a general mobility"

Source: https://www.emergentmind.com/topics/stabilized-imex-crank-nicolson-method