Abstract:This paper is concerned with numerical methods for a class of nonhomogeneous Neumann problems of time-fractional reaction–diffusion equations with variable coefficients. The solutions of this kind of problems often have weak singularity at the initial time. This makes the existing numerical methods with uniform time mesh often lose accuracy. In this paper, we propose and analyze a high-order compact finite difference method with nonuniform time mesh. The time-fractional derivative is approximated by Alikhanov’s high-order approximation on a class of fitted time meshes. For the spatial variable coefficient differential operator, a new fourth-order boundary discretization is developed under the nonhomogeneous Neumann boundary condition, and then a new fourth-order compact finite difference approximation on a space uniform mesh is obtained. Under the assumption of the weak initial singularity of solution, we prove that for the general case of the variable coefficients, the proposed method is unconditionally stable and the numerical solution converges to the solution of the problem under consideration. The convergence result also gives an optimal error estimate of the numerical solution in the discrete L2-norm, which shows that the method has the spatial fourth-order convergence, while it attains the temporal optimal second-order convergence provided a proper mesh grading parameter is employed. Numerical results that confirm the sharpness of the error analysis are presented.

A Splitting Moving Mesh Method for Reaction-Diffusion Equations of Quenching Type

High Accuracy Non-Equidistant Method for Singular Perturbation Reaction-Diffusion Problem

Numerical Solution of Quenching Problems Using Mesh-Dependent Variable Temporal Steps

A moving mesh method for quenching problem

Numerical quenching for a reaction-diffusion system with coupled singular absorption terms

The Quenching Behavior of a General Singular Electrostatic Micro-Electro-Mechanical-System

Quenching for a Reaction-Diffusion System with Coupled Inner Singular Absorption Terms

Quenching for a reaction-diffusion equation with nonlinear memory

A high-order compact finite difference method on nonuniform time meshes for variable coefficient reaction–subdiffusion problems with a weak initial singularity

Simple Positivity Preserving Nonlinear Finite Volume Scheme for Subdiffusion Equations on General Non-conforming Distorted Meshes

Moving Mesh Method for Problems with Blow-Up on Unbounded Domains

Quenching for a reaction–diffusion system with logarithmic singularity

A High-Order Compact Difference Method on Fitted Meshes for Neumann Problems of Time-Fractional Reaction–diffusion Equations with Variable Coefficients

Non-Simultaneous Quenching for A Slow Diffusion System Coupled at the Boundary

Quenching Rates for Heat Equations with Coupled Singular Nonlinear Boundary Flux

Quenching Behavior of Solutions to Heat Equations with Coupled Boundary Singularities.

Novel numerical methods based on graded, adaptive and uniform meshes for a time-fractional advection-diffusion equation subjected to weakly singular solution

Adaptive temporal mesh equidistribution for studying the singularly perturbed delay-in-time parabolic equations

High-order splitting finite element methods for the subdiffusion equation with limited smoothing property

On the Quenching Behavior of the Mems with Fringing Field

Quenching Behavior of Solutions in Coupled Heat Equations with Singular Multi-Nonlinearities.