Runge--Kutta Based gCQ
- Runge--Kutta based generalized Convolution Quadrature (gCQ) is a high-order, stable method for approximating convolution integrals, notably handling arbitrary time meshes unlike classical CQ.
- This method excels in solving time-dependent problems with memory and sectorial kernels, such as fractional evolution equations and boundary integral equations, especially with singular initial data.
- gCQ allows for fast O(N log N) implementations with low O(log N) memory usage, making it practical for long-time simulations on complex grids.
Runge–Kutta based generalized Convolution Quadrature (gCQ) is a family of high-order, stable, and computationally efficient methods for approximating convolution integrals involving causal, sectorial kernels—fundamental to the numerical treatment of time-dependent boundary integral equations, fractional evolution equations, and parabolic problems with memory. In contrast to Lubich’s classical convolution quadrature (CQ), which is restricted to uniform time grids and fixed-step multistep or Runge–Kutta methods, gCQ allows for arbitrary time meshes (including nonuniform and graded meshes), optimal accuracy for singular data, and fast memory- and computation-efficient implementation. The gCQ approach is especially potent when based on A-stable Runge–Kutta methods, enabling reliable simulation of problems with highly irregular or singular behavior at initial times.
1. Generalization of Convolution Quadrature and Time Mesh Flexibility
Classical convolution quadrature methods, as developed by Lubich, approximate causal convolutions
by applying a linear multistep or Runge–Kutta time-stepping scheme to an associated operational equation involving the Laplace transform of , typically on a uniform time grid. The gCQ framework extends this paradigm to arbitrary, potentially highly nonuniform meshes: allowing for mesh grading or local time refinement tailored to singularities in the input data or solution.
This generalization is crucial for practical applications where solution regularity is lost at initial times—such as subdiffusion, wave propagation in non-smooth domains, or viscoelastic models—since uniform time stepping leads to significant order reduction. gCQ leverages mesh flexibility to restore optimal convergence.
2. High-Order Runge–Kutta Integration in gCQ
Runge–Kutta based gCQ replaces the discrete convolution weights of classical CQ with those generated by high-order, A-stable Runge–Kutta schemes applied to the ODE
for each parameterizing the Laplace variable, reflecting the Laplace inversion representation of sectorial kernels.
Given a multistage Runge–Kutta method with nodes and coefficients satisfying A-stability, classical order , stage order , and stiff accuracy, the method propagates the solution in time according to the possibly variable time steps , computing at each . The numerical approximation to the target convolution is then
with a function determined by the sectorial Laplace transform of .
This structure allows the use of high-order RK schemes such as Radau IIA or Lobatto IIIC, ensuring superior convergence and stability properties compared to lower-order, single-step CQ (such as backward Euler).
3. Sectorial Problem Classes and Analytical Foundations
The method's strongest theoretical guarantees and performance are realized for so-called sectorial problems—those where the Laplace transform of the kernel is holomorphic for (in practice, for a sector ) and satisfies decay conditions
This encompasses kernels from fractional integrals, subdiffusion, boundary integral operators for parabolic PDEs, and many memory effects.
A central insight is that such admit a real inversion integral
with
enabling accurate and efficient quadrature for the convolution and the design of the gCQ scheme.
4. Convergence, Stability, and Graded Mesh Optimization
gCQ methods based on suitably chosen RK schemes preserve the classical convergence order of CQ (, the order of the underlying scheme), even for arbitrary time grids and for data with algebraic singularities , provided optimal mesh grading is enforced. Theoretical results (Theorems 1 and 2 in the source) state:
- For sufficiently smooth data and sectorial kernels:
- For singular data and a graded mesh with
the scheme nearly attains the optimal order
where is the stage order, the kernel decay, and the strength of singularity in the data.
This refined control is unattainable for uniform-step CQ and is decisive for accurate simulation of memory- or singular-data-driven phenomena.
5. Fast and Memory-Efficient Implementation
A central advantage of gCQ is its amenability to fast, large-scale simulation. The convolution weights can be computed using the "fast and oblivious" algorithm, which decomposes history interactions into "local" (recent) and "history" (distant past) contributions. The latter exploits the real inversion formula for the kernel and makes use of analytical or quadrature-based recursion, yielding the following computational features:
- operations for time steps,
- memory requirements, even on nonuniform and graded grids.
These features make gCQ practical for long-time integration and for problems with significant history-dependence, such as those found in viscoelasticity and time-domain boundary integral equations.
6. Applications and Numerical Examples
Runge–Kutta based gCQ has demonstrable performance across multiple challenging application domains:
- Parabolic and subdiffusive PDEs with singular initial data,
- Time-dependent boundary integral equations (for instance, in time-domain BEM discretizations of wave problems),
- Nonlinear memory-type equations, illustrated by damped nonlinear wave equations and the Westervelt equation,
- Fractional-order problems where uniform CQ order would be significantly degraded.
Numerical experiments confirm that, with sectorial kernels and optimal mesh grading, full theoretical convergence rates are achieved and order reduction is avoided. These experiments are reproducible with the provided MATLAB implementation (github.com/jingguo-math/gCQrK).
7. Summary Table: Method Properties
Feature/Aspect | Runge–Kutta gCQ | Classical CQ (BE/multistep) |
---|---|---|
Time mesh | Arbitrary (graded/adaptive) | Uniform only |
Achievable order (for sectorial) | Up to RK order | Up to scheme order |
Order with singular data | Full with optimal grading | Reduced |
Stability | Optimal (for A-stable RK) | Good (for A-stable multistep) |
Memory usage | (fast, oblivious) | (classic) |
Computational cost | (classic) |
References
The methodology and theoretical results are as given in (2506.21242), with background in Lubich (1988), LoSau16, BanLo19, and closely related literature on RK-based time stepping and convolution quadrature.
Runge–Kutta based generalized Convolution Quadrature provides a high-order, stable, and flexible approach to discretizing convolution and memory terms on arbitrary time grids for sectorial problems, optimizing convergence for both smooth and singular data, and enabling efficient long-time simulation through advanced quadrature and recursion algorithms.