A semi-adaptive preservative scheme for a fractional quenching convective-diffusion problem

Abstract

A preservative scheme is presented and analyzed for the solution of a quenching type convective-diffusion problem modeled through one-sided Riemann-Liouville space-fractional derivatives. Properly weighted Grünwald formulas are employed for the discretization of the fractional derivative. A forward difference approximation is considered in the approximation of the convective term of the nonlinear equation. Temporal steps are optimized via an asymptotic arc-length monitoring mechanism till the quenching point. Under suitable constraints on spatial-temporal discretization steps, the monotonicity, positivity preservations of the numerical solution and numerical stability of the scheme are proved. Three numerical experiments are designed to demonstrate and simulate key characteristics of the semi-adaptive scheme constructed, including critical length, quenching time and quenching location of the fractional quenching phenomena formulated through the one-sided space-fractional convective-diffusion initial-boundary value problem. Effects of the convective function to quenching are discussed. Numerical estimates of the order of convergence are obtained. Computational results obtained are carefully compared with those acquired from conventional integer order quenching convection-diffusion problems for validating anticipated accuracy. The experiments have demonstrated expected accuracy and feasibility of the new method.