Abstract:Solutions exhibiting weak initial singularities arise in various equations, including diffusion and subdiffusion equations. When employing the well-known L1 scheme to solve subdiffusion equations with weak singularities, numerical simulations reveal that this scheme exhibits varying convergence rates for different choices of model parameters (i.e., domain size, final time $T$, and reaction coefficient $\kappa$). This elusive phenomenon is not unique to the L1 scheme but is also observed in other numerical methods for reaction-diffusion equations such as the backward Euler (IE) scheme, Crank-Nicolson (C-N) scheme, and two-step backward differentiation formula (BDF2) scheme. The existing literature lacks an explanation for the existence of two different convergence regimes, which has puzzled us for a long while and motivated us to study this inconsistency between the standard convergence theory and numerical experiences. In this paper, we provide a general methodology to systematically obtain error estimates that incorporate the exponential decaying feature of the solution. We term this novel error estimate the `decay-preserving error estimate' and apply it to the aforementioned IE, C-N, and BDF2 schemes. Our decay-preserving error estimate consists of a low-order term with an exponential coefficient and a high-order term with an algebraic coefficient, both of which depend on the model parameters. Our estimates reveal that the varying convergence rates are caused by a trade-off between these two components in different model parameter regimes. By considering the model parameters, we capture different states of the convergence rate that traditional error estimates fail to explain. This approach retains more properties of the continuous solution. We validate our analysis with numerical results.

Renormalization Group Method for Singular Perturbation Initial Value Problems with Delays

Renormalization Group Approach to a Class of Singularly Perturbed Delay Differential Equations

Delay-Dependent Robust Stability Of Uncertain Discrete Singular Time-Delay Systems With Nonlinear Disturbance

Delay-dependent Robust Stability of Uncertain Discrete Singular Time-Delay Systems

Robust Stabilization For Singular Systems With Time-Delays And Saturating Controls

Robust Stability and Stabilization of Uncertain Discrete Singular Time-Delay Systems Based on PNP Lyapunov Functional

New Delay-Dependent Criteria for Robust Stability of Uncertain Singular Systems ?

Robust Control For Singular Systems With State Delay And Parameter Uncertainty

Singular LQ Sub-Optimal Control Problem for a Class of Discrete-Time Singular Systems with Time-Delay

Analysis of robust stability for uncertain discrete singular systems with time-varying delay

Renormalization group method for singular perturbation system with additive fractional Gaussian noise

Numerical Analysis for a Singularly Perturbed Parabolic Differential Equation with a Time Delay

Robust Stabilization for Discrete-Time Nonlinear Singular Systems with Mixed Time Delays

Delay-Dependent Robust Exponential Stability Of Uncertain Singular Systems With Time Delays

Delay-dependent Robust Stabilization for Uncertain Singular Systems with Discrete and Distributed Delays

A robust numerical method for a two-parameter singularly perturbed time delay parabolic problem

A variational-perturbation method for solving the time-dependent singularly perturbed reaction-diffusion problems

Improved results on reachable set estimation for singularly perturbed systems with time-varying delay

Stability and Bifurcation in Delay–Differential Equations with Two Delays

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

A uniformly convergent difference method for singularly perturbed parabolic partial differential equations with large delay and integral boundary condition