Efficient bound preserving and asymptotic preserving semi-implicit schemes for the fast reaction-diffusion system

TL;DR

This paper develops efficient semi-implicit numerical schemes for fast reaction-diffusion systems that preserve bounds and accurately capture interface dynamics as reaction rates become extremely large.

Contribution

It introduces first- and second-order semi-implicit schemes that are positivity-preserving, stable, and effective for simulating singular limit behaviors in fast reaction-diffusion systems.

Findings

Schemes accurately capture interface propagation for small epsilon.

Numerical tests confirm stability, accuracy, and bound preservation.

Effective simulation of heat transfer and phase change processes.

Abstract

We consider a special type of fast reaction-diffusion systems in which the coefficients of the reaction terms of the two substances are much larger than those of the diffusion terms while the diffusive motion to the substrate is negligible. Specifically speaking, the rate constants of the reaction terms are $O(1/\epsilon)$ while the diffusion coefficients are $O(1)$ where the parameter $\epsilon$ is small. When the rate constants of the reaction terms become highly large, i.e. $\epsilon$ tends to 0, the singular limit behavior of such a fast reaction-diffusion system is inscribed by the Stefan problem with latent heat, which brings great challenges in numerical simulations. In this paper, we adopt a semi-implicit scheme, which is first-order accurate in time and can accurately approximate the interface propagation even when the reaction becomes extremely fast, that is to say, the parameter $\epsilon$ is sufficiently small. The scheme satisfies the positivity, bound preserving properties and has $L^2$ stability and the linearized stability results of the system. For better performance on numerical simulations, we then construct a semi-implicit Runge-Kutta scheme which is second-order accurate in time. Numerous numerical tests are carried out to demonstrate the properties, such as the order of accuracy, positivity and bound preserving, the capturing of the sharp interface with various $\epsilon$ and to simulate the dynamics of the substances and the substrate, and to explore the heat transfer process, such as solid melting or liquid solidification in two dimensions.