- Quenching singularity,
- Dimensional splitting,
- Mesh adaptation,
- Positivity,
- Monotonicity,
- Numerical stability
This paper studies the numerical solution of a two-dimensional quenching type nonlinear reaction-diffusion problem via dimensional splitting. The variable coefficient differential equation considered possesses a nonlinear forcing term, and may lead to strong quenching singularities that have profound multiphysics and engineering applications to the energy industry. Our investigations focus on the construction and preservative properties of a semi-adaptive Peaceman-Rachford procedure for solving the aforementioned problem. In the study, the multidimensional differential equation is decomposed into single dimensional subequations so that the computational efficiency can be effectively raised. Temporal adaptation is implemented through arc-length estimations of the rate-of-change of the numerical solution. The positivity, monotonicity and localized linear stability of the variable step splitting scheme are analyzed. These ensure key preservations of underlying physical features of the computed approximations for applications. Computational experiments are presented to illustrate our results as well as to demonstrate the viability and accuracy of the splitting method for solving singular quenching-combustion problems.
Available at: http://works.bepress.com/julienne-kabre/1/
The authors would like to thank the anonymous referees for their time spent and extremely valuable remarks offered. Their suggestions have significantly improved the quality and presentation of this paper. The second author wishes to acknowledge a partial support from Baylor University through a 2021 Research Leave Award (No. BU-PP714-20).
Last but not least, the authors would also thank the editor for the tremendous amount of encouragement received throughout the preparation of this article.