zbMATH — the first resource for mathematics

Solving degenerate reaction-diffusion equations via variable step Peaceman-Rachford splitting. (English) Zbl 1061.65086
Summary: The authors study the numerical solution of two-dimensional nonlinear degenerate reaction-diffusion differential equations with singular forcing terms over rectangular domains. The equations considered may generate strong quenching singularities. This investigation focuses on a variable time step Peaceman-Rachford splitting method for the aforementioned problem.
The time adaptation is implemented based on arc-length estimations of the first time derivative of the solution. The two-dimensional problem is split into several one-dimensional problems so that the computational cost is significantly reduced. The monotonicity and localized linear stability of the variable step scheme are investigated.
We give some numerical examples to illustrate our results as well as to demonstrate the viability and efficiency of the method over existing methods for the quenching problem. It is also shown that the numerical solution obtained preserves important properties of the physical solution of the given problem.
Reviewer: Reviewer (Berlin)

65M20 Method of lines for initial value and initial-boundary value problems involving PDEs
65M50 Mesh generation, refinement, and adaptive methods for the numerical solution of initial value and initial-boundary value problems involving PDEs
35K65 Degenerate parabolic equations
35K57 Reaction-diffusion equations
Full Text: DOI