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.

IMEX schemes for a Parabolic-ODE system of European Options with Liquidity Shocks

Numerical analysis and simulation of European options under liquidity shocks: A coupled semilinear system approach with new IMEX methods

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

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

IMEX-RK finite volume methods for nonlinear 1d parabolic PDEs. Application to option pricing

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

An explicitness-preserving IMEX-split multiderivative method

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

On a class of higher-order length preserving and energy decreasing IMEX schemes for the Landau-Lifshitz equation

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

Implicit-explicit time discretization schemes for a class of semilinear wave equations with nonautonomous dampings

Implicit-explicit schemes for incompressible flow problems with variable viscosity

Implicit-Explicit WENO schemes for the equilibrium dispersive model of chromatography

Parallel-in-time high-order multiderivative IMEX solvers

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

On solutions of a partial integro-differential equation in Bessel potential spaces with applications in option pricing models

Efficient numerical valuation of European options under the two-asset Kou jump-diffusion model

A mixed Monte Carlo and PDE variance reduction method for foreign exchange options under the Heston-CIR model

Implicit-explicit Runge–Kutta methods for pricing financial derivatives in state-dependent regime-switching jump-diffusion models

Fourth order compact scheme for option pricing under Merton and Kou jump-diffusion models

Implicit-explicit Crank-Nicolson scheme for Oseen's equation at high Reynolds number