Abstract:The recently developed technique of DOC kernels has been a great success in the stability and convergence analysis for BDF2 scheme with variable time steps. However, such an analysis technique seems not directly applicable to problems with initial singularity. In the numerical simulations of solutions with initial singularity, variable time-steps schemes like the graded mesh are always adopted to achieve the optimal convergence, whose first adjacent time-step ratio may become pretty large so that the acquired restriction is not satisfied. In this paper, we revisit the variable time-step implicit-explicit two-step backward differentiation formula (IMEX BDF2) scheme presented in [W. Wang, Y. Chen and H. Fang, \emph{SIAM J. Numer. Anal.}, 57 (2019), pp. 1289-1317] to compute the partial integro-differential equations (PIDEs) with initial singularity. We obtain the sharp error estimate under a mild restriction condition of adjacent time-step ratios $r_{k}: =\tau_{k}/\tau_{k-1} \; (k\geq 3) < r_{\max} = 4.8645 $ and a much mild requirement on the first ratio, i.e., $r_2>0$. This leads to the validation of our analysis of the variable time-step IMEX BDF2 scheme when the initial singularity is dealt by a simple strategy, i.e., the graded mesh $t_k=T(k/N)^{\gamma}$. In this situation, the convergence of order $\mathcal{O}(N^{-\min\{2,\gamma \alpha\}})$ is achieved with $N$ and $\alpha$ respectively representing the total mesh points and indicating the regularity of the exact solution. This is, the optical convergence will be achieved by taking $\gamma_{\text{opt}}=2/\alpha$. Numerical examples are provided to demonstrate our theoretical analysis.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to conduct accurate error estimation for the implicit - explicit two - step backward difference method (IMEX BDF2 scheme) with variable time steps when numerically solving parabolic integro - differential equations (PIDEs) with initial singularity. Specifically, the article focuses on how to ensure the optimal convergence of the numerical method when using a graded mesh in the case of non - smooth initial conditions. ### Background and Problem Description 1. **Research Background**: - Parabolic integro - differential equations (PIDEs) have important applications in option pricing theory in finance. - When the underlying asset follows a jump - diffusion process, the solution of PIDEs may exhibit weak singularity at the initial time $ t = 0 $, that is, the regularity of the solution changes as time increases. 2. **Challenges in Numerical Solution**: - For PIDEs with initial singularity, traditional fixed - time - step methods cannot achieve the optimal convergence rate. - Variable - time - step methods (such as graded meshes) can improve convergence, but require strict limitations on the ratio of adjacent time steps. - In the existing literature, the stability analysis of the BDF2 scheme with variable time steps usually depends on the assumption of smooth solutions and may be invalid when dealing with initial singularity. ### Core Problems of the Paper The main objective of the paper is to re - examine and improve the variable - time - step IMEX BDF2 scheme to achieve accurate error estimation for PIDEs with initial singularity. Specifically: - **Error Estimation**: The paper aims to obtain sharp error estimates under relatively loose time - step - ratio conditions (for example, $ r_k\leq4.8645 $). - **Initial Singularity Handling**: By introducing a graded mesh $ t_k = T(k/N)^\gamma $, ensure that the optimal convergence can still be maintained when dealing with initial singularity. - **Influence of the First - Level Time - Step Ratio**: Analyze the influence of the first - level time - step ratio $ r_2 $ on the error bound and propose a more general analysis framework. ### Main Contributions 1. **Extension of Time - Step - Ratio Conditions**: Extend the limitation of the ratio of adjacent time steps from $ r_k\leq1.91 $ to $ r_k\leq4.8645 $, and there are no additional restrictions on the first - level time - step ratio $ r_2 $. 2. **Optimal Convergence**: Under the graded mesh, when $ \gamma = 2/\alpha $, the optimal convergence of $ O(N^{-2}) $ is achieved. 3. **Relationship between Error and $ r_2 $**: Clarify the polynomial relationship between the error and the increase of the first - level time - step ratio $ r_2 $, which is consistent with the numerical experiment results. Through these improvements, the paper provides a more effective and accurate numerical method for dealing with PIDEs with initial singularity.

Sharp error estimate of variable time-step IMEX BDF2 scheme for parabolic integro-differential equations with initial singularity arising in finance

Do we need decay-preserving error estimate for solving parabolic equations with initial singularity?

An Efficient IMEX Compact Scheme for the Coupled Time Fractional Integro-Differential Equations Arising from Option Pricing with Jumps

Optimal error analysis of a non-uniform IMEX-L1 finite element method for time fractional PDEs and PIDEs

On a new class of BDF and IMEX schemes for parabolic type equations

Unconditionally optimal error estimate of a linearized variable-time-step BDF2 scheme for nonlinear parabolic equations

Convergence and Superconvergence Analysis of a Nonconforming Finite Element Variable-Time-Step BDF2 Implicit Scheme for Linear Reaction-Diffusion Equations

A second order finite volume IMEX Runge-Kutta scheme for two dimensional PDEs in finance

On a Non-Uniform $α$-Robust IMEX-L1 Mixed FEM for Time-Fractional PIDEs

Stability and convergence of BDF2-ADI schemes with variable step sizes for parabolic equation

Enhancing accuracy with an adaptive discretization for the non-local integro-partial differential equations involving initial time singularities

Implicit-Explicit difference schemes for nonlinear fractional differential equations with non-smooth solutions

Second order finite volume IMEX Runge-Kutta schemes for two dimensional parabolic PDEs in finance

Uniform convergence analysis of the BDF2 scheme on Bakhvalov-type meshes for a singularly perturbed Volterra integro-differential equation

A Compact Difference Scheme on Graded Meshes for the Nonlinear Fractional Integro-differential Equation with Non-smooth Solutions

An explicitness-preserving IMEX-split multiderivative method

Convergence analysis of high-order IMEX-BDF schemes for the incompressible Navier–Stokes equations

A fourth-order accurate extrapolation nonlinear difference method for fourth-order nonlinear PIDEs with a weakly singular kernel

Analysis of a higher-order scheme for multi-term time-fractional integro-partial differential equations with multi-term weakly singular kernels

Second-order accurate numerical scheme with graded meshes for the nonlinear partial integrodifferential equation arising from viscoelasticity

Stability and error estimates for the variable step-size BDF2 method for linear and semilinear parabolic equations