---
title: Lax–Wendroff Type Methods for Hyperbolic PDEs
url: https://www.emergentmind.com/topics/lax-wendroff-lw-type-method
type: topic
---

# Lax–Wendroff Type Methods for Hyperbolic PDEs

A Lax–Wendroff (LW) type method refers broadly to a class of time integration techniques for hyperbolic partial differential equations (PDEs) or systems of conservation laws, in which the solution is advanced using single-step, typically single-stage, Taylor expansions in time with the PDE recursively used to convert time derivatives to spatial derivatives. This approach, formalized by Lax and Wendroff, has inspired a wide spectrum of high-order, compact algorithms including Taylor-series, ADER (Arbitrary DERivatives), and modern flux reconstruction/discontinuous Galerkin methods. LW-type methods are characterized by strong space–time coupling, explicit construction of high-order temporal derivatives, and, when cast in appropriate frameworks, offer high efficiency and stability compared to conventional multi-stage Runge–Kutta (RK) temporal integrators.

## 1. Core Principles and Mathematical Structure

LW-type schemes originate from the formal Taylor expansion of the solution in time,
\[
u^{n+1}(x) = u^n(x) + \Delta t \, \partial_t u^n(x) + \frac{\Delta t^2}{2} \partial_{tt} u^n(x) + \ldots,
\]
where, for a typical conservation law,
\[
\partial_t u + \partial_x f(u) = 0,
\]
succussive time derivatives of $u$ are recursively expressed as
\[
\partial_t u = -\partial_x f(u), \qquad
\partial_{tt} u = -\partial_t \partial_x f(u),
\]
with the chain- and product-rule recastings leading to all higher time derivatives replaced by combinations of spatial derivatives acting on $u$ and its images under $f$ (Cauchy–Kovalevskaya procedure). This provides a direct update of the form
\[
u^{n+1} = u^n - \sum_{m=1}^{N+1} \frac{\Delta t^m}{m!} \partial_x [\partial_t^{m-1} f(u^n)] + O(\Delta t^{N+2}),
\]
with order controlled by the number of Taylor terms retained.

Notably, in the context of spectral element, flux reconstruction (FR), and discontinuous Galerkin (DG) methods, the time-averaged "Lax–Wendroff" flux polynomial
\[
F_h^\delta(\xi) = a \sum_{k=0}^N \frac{(-a \Delta t)^k}{(k+1)!} \partial_x^k u_h^n(\xi),
\]
or its quadrature- and matrix-free equivalents, is central to the update [2402.18937, 2207.02954].

## 2. LW-Type Methods in High-Order Element Frameworks

LW-type ideas underpin not only classical finite-difference schemes but are intrinsic to state-of-the-art high-order DG and FR frameworks. In these settings:

- Each element is discretized by a local polynomial basis (e.g., Lagrange through Gauss–Legendre nodes);
- The time-stepped update is written as
  \[
  u_j^{n+1} = u_j^n - \Delta t \, \partial_x F_h(\xi_j),
  \]
  where $F_h(\xi)$ is the reconstructed, (typically globally continuous) time-averaged flux polynomial;
- Continuity is imposed via correction functions acting on the difference between the (possibly discontinuous) polynomial flux and appropriately defined numerical fluxes at element interfaces (e.g., using Radau polynomials);
- For general nonlinear systems, a "Jacobian-free" approach is preferred: higher time derivatives of the flux are approximated via symmetric finite difference stencils, requiring only pointwise flux evaluations [2207.02954, 2506.11719].

The ADER method computes an element-local space–time predictor (solving a small implicit multi-variable system), but for linear problems and with the D2 numerical flux (see below), ADER–DG and Lax–Wendroff–FR become algebraically identical [2402.18937].

## 3. Numerical Fluxes and D2 Dissipation

A critical aspect for the stability and equivalence of LW-type schemes to ADER is the construction of the interface numerical flux. The so-called D2 dissipation upwind flux,
\[
F_{e+\frac12} = \frac12\bigl(F_h^\delta(1^-)+F_h^\delta(0^+)\bigr) - \frac{|a|}2\bigl(U_h^n(0^+)-U_h^n(1^-)\bigr),
\]
operates on the time-averaged solution polynomial $U_h^n$. It ensures that the LWFR update coincides with ADER-DG at all polynomial degrees for linear advection and provides the same Fourier stability limit,
\[
\text{CFL}_{\max} \simeq \frac{1}{2N+1},
\]
permitting maximal explicit time steps [2402.18937].

Deviation from D2 (e.g., using only classical upwind on the instantaneous state) results in significant error divergence. This equivalence under D2 is proved through analytical normal-mode analysis and verified by numerical experiments.

## 4. Stability, Dispersion, and DRP Analysis

The spatio-temporal coupling inherent to LW-type methods provides higher stability and dispersion-relation-preserving (DRP) properties compared to stage-based Runge–Kutta schemes for the same formal order. Fourier analysis indicates:

- For the linear advection model, the LW update is stable for CFL up to unity in first-order finite-difference realizations, and to the scheme-specific limit in DG/FR settings;
- Spectral analysis (GSA/DRP) shows that LW-type compact schemes maintain phase and amplitude accuracy over a larger wavenumber band, yielding lower points-per-wavelength requirements for DNS/LES fidelity compared to method-of-lines RK schemes [2208.09428, 2208.09808, 2601.17264].

In the MSMD (multi-stage multi-derivative) perspective, each additional time derivative (spatially constructed) directly lifts the dispersion and stability properties, allowing for higher CFL, more efficient evolution, and lower dissipation in the resolved spectrum.

## 5. Boundary Treatments: Inverse Lax–Wendroff Procedures

Enforcement of high-order boundary conditions in LW-type schemes is accomplished via the inverse Lax–Wendroff (ILW) procedure. Given Dirichlet or characteristic data at the inflow, ghost point values (needed for wide stencils or high-order interpolations) are constructed by:

- Taylor expansion in the normal direction to the boundary;
- Replacing spatial derivatives by time derivatives (and, recursively, by further spatial derivatives using the PDE) up to the desired formal order;
- In modern extensions, a hybrid approach (least-squares or Hermite polynomial fits) is employed to minimize the number of derivatives directly generated from the PDE, further improving computational efficiency and stability, particularly for high-order and multi-dimensional problems [1209.3103, 2402.10152, 2503.21626];

Characteristic-based variants extend this to systems; for strong shocks or non-standard inflow configurations, WENO interpolants and local indicator-based stencil adaptation ensure robust, non-oscillatory boundary closure [1801.00990].

## 6. Extensions: Admissibility, Limiters, Blending, and Source Terms

LW-type methods, when employed for nonlinear systems such as the compressible Euler or relativistic hydrodynamics equations, must enforce admissibility constraints (e.g., positivity of density and pressure). State-of-the-art realizations achieve this via:

- Flux blending: the convex combination of high-order LWFR and robust low-order (e.g., Lax–Friedrichs, finite volume) fluxes on subcells, controlled by smoothness/discontinuity indicators [2305.10781, 2409.15805, 2402.11926];
- Zhang–Shu scaling limiters: ensuring that the solution at all nodal or quadrature points remains within the physical admissible set [1601.08145, 2402.01442, 2409.15805];
- Special strategies for schemes on curvilinear or adaptive mesh refinement grids, including metric-consistent flux reconstruction and error-controlled embedded time stepping [2402.11926];
- Capability for extension to balance laws with source terms, employing Taylor expansions for both flux and source term time derivatives, and flux decomposition followed by admissibility-preserving blending [2402.01442].

Automatic differentiation (AD) provides an alternative to finite-difference stencils for high-order flux time derivatives, being Jacobian-free, efficient, and naturally compatible with positivity constraints [2506.11719].

## 7. Applications, Comparisons, and Ongoing Developments

LW-type methods and their generalizations are now deployed for a diverse range of applications:

- Direct and large eddy simulation (DNS/LES) of complex flows, where their DRP properties give actionable wavenumber and CFL resolution guidelines [2208.09428, 2208.09808];
- Astrophysical flows, relativistic hydrodynamics, plasma physics, and kinetic equations, including sophisticated admissibility- and subcell-limited LWFR methods [2409.15805];
- Structural dynamics in geophysics, with finite-element implementations utilizing Lax–Wendroff/interpolation time stepping for elastodynamics [2502.07797];

Comparison to traditional RK-DG and method-of-lines schemes consistently reveals that, for a given spatial and temporal accuracy, single-stage LW-type methods achieve equivalent or higher accuracy and efficiency, with superior stability limits and reduced dissipation/dispersion at high frequencies [1512.03664, 2601.17264].

### Table: Core Features of Representative LW-Type Schemes

| Variant        | High-Order, Single-Step | Compact Stencil | DRP/Optimal CFL | Admissibility, Limiting        |
|---------------|------------------------|-----------------|----------------|-------------------------------|
| Classical LW  | Y                      | N               | Moderate       | Limited (ad-hoc)              |
| LW-FR/DG      | Y                      | Y               | High           | Blending, Scaling Limiters    |
| ADER-DG       | Y (via predictor)      | Y               | High           | Blending, Scaling Limiters    |
| AD-LW         | Y (automatic diff.)    | Y               | High           | Natural with positivity fix   |
| MSMD-LW       | Y                      | Y               | Highest        | As above                      |

All these developments emphasize the centrality of LW-type (single-stage, space–time coupled) methods as the backbone for efficient, high-fidelity, robust numerical solvers in scientific computing and computational fluid dynamics [2402.18937, 2207.02954, 1601.08145, 2402.11926, 2409.15805]. The choice among the various implementations—exact Taylor, approximate via finite differences, space–time predictor, AD—depends on application context, implementation priorities (e.g., Jacobian-free requirement), and complexity of the physical system addressed.

For nonlinear systems or situations involving boundaries and shocks, the integration of hybrid ILW-type boundary closures and admissibility-preserving subcell limiters has become standard to ensure stability and fidelity across all regimes of interest.

Source: https://www.emergentmind.com/topics/lax-wendroff-lw-type-method