A semi-adaptive finite difference method for simulating two-sided fractional convection–diffusion quenching problems

This paper investigates quenching solutions of an one-dimensional, two-sided Riemann–Liouville fractional order convection–diffusion problem. Fractional order spatial derivatives are discretized using weighted averaging approximations in conjunction with standard and shifted Grünwald formulas. The advective term is handled utilizing a straightforward Euler formula, resulting in a semi-discretized system of nonlinear ordinary differential equations. The conservativeness of the proposed scheme is rigorously proved and validated through simulation experiments. The study is further advanced to a fully discretized, semi-adaptive finite difference method. Detailed analysis is implemented for the monotonicity, positivity and stability of the scheme. Investigations are carried out to assess the potential impacts of the fractional order on quenching location, quenching time, and critical length. The computational results are thoroughly discussed and analyzed, providing a more comprehensive understanding of the quenching phenomena modeled through two-sided fractional order convection–diffusion problems.