Stable rotational symmetric schemes for nonlinear reaction-diffusion equations
COMPUTERS & MATHEMATICS WITH APPLICATIONS(2022)
摘要
In the simulation of biological pattern forming, it has been observed that the numerical solution is more sensitive to the spatial mesh resolution than the temporal one. Such a higher sensitivity to the spatial resolution is mainly originated from an inaccurate approximation of diffusion differential operators, which might violate the rotational symmetry to be seriously erroneous in low spatial resolutions. Also, it has been known that the second-order Crank-Nicolson time-stepping procedure may introduce spurious oscillations when the initial data or the source term is nonsmooth and the temporal step size is set relatively large. This article studies 9-point finite difference schemes for the diffusion operator to enhance the rotational symmetry, employs the variable-theta method to achieve a nonoscillatory second-order time-stepping procedure, and adopts an effective relaxation linear solver to solve the algebraic systems efficiently. The variable-theta method is proved to satisfy the maximum principle, which guarantees that the time-stepping procedure is unconditionally stable. When the successive over-relaxation method with an optimal relaxation parameter is adopted for the algebraic solver, the iteration converges in 2-4 iterations in most time steps. The overall algorithm is second-order in accuracy and scalable in efficiency. Various examples are given to show the accuracy and efficiency of the proposed algorithm for the numerical solution of the system of nonlinear reaction-diffusion equations.
更多查看译文
关键词
Reaction-diffusion equations,Biological pattern formation,Time-stepping procedure,Variable-theta method,Rotational symmetry,Averaging scheme
AI 理解论文
溯源树
样例
![](https://originalfileserver.aminer.cn/sys/aminer/pubs/mrt_preview.jpeg)
生成溯源树,研究论文发展脉络
Chat Paper
正在生成论文摘要