2022/01/23 by Chengchao Zhao, Ruoyu Yang, Zhao, Chengchao +5
Engineering · Mathematics · #65M06 #65M12 #Advanced Numerical Methods in Computational Mathematics #Differential Equations and Numerical Methods #FOS: Mathematics #Numerical Analysis (math.NA) #Numerical methods for differential equations
paper · pdf · doi:10.48550/arxiv.2201.09322
openalex publication_date 2022/01/23 · openalex created_date 2022/04/03 · openalex updated_date 2026/07/28
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, 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 rk: =τk/τk-1 (k≥ 3) < rmax = 4.8645 and a much mild requirement on the first ratio, i.e., r2>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 tk=T(k/N)γ. In this situation, the convergence of order O(N^-min\2,γα\) is achieved with N and α respectively representing the total mesh points and indicating the regularity of the exact solution. This is, the optical convergence will be achieved by taking γopt=2/α. Numerical examples are provided to demonstrate our theoretical analysis.