A robust and efficient method for steady state patterns in reaction-diffusion systems

JOURNAL OF COMPUTATIONAL PHYSICS

卷:231

Article

Lo, Wing-Cheong^1,2,3,4  Chen, Long^1,2,3  Wang, Ming^1,5  Nie, Qing^1,2,3 

[1] Univ Calif Irvine, Dept Math, Irvine, CA 92717 USA

[2] Univ Calif Irvine, Ctr Complex Biol Syst, Irvine, CA 92717 USA

[3] Univ Calif Irvine, Ctr Math & Computat Biol, Irvine, CA 92717 USA

[4] Ohio State Univ, Math Biosci Inst, Columbus, OH 43210 USA

[5] Peking Univ, Sch Math Sci, LAMA, Beijing 100871, Peoples R China

关键词: Reaction-diffusion;    Steady state;    Patterns;    Hybrid methods;   

DOI  :  10.1016/j.jcp.2012.04.006

来源: Elsevier

【 摘 要 】

An inhomogeneous steady state pattern of nonlinear reaction-diffusion equations with no-flux boundary conditions is usually computed by solving the corresponding time-dependent reaction-diffusion equations using temporal schemes. Nonlinear solvers (e.g., Newton's method) take less CPU time in direct computation for the steady state; however, their convergence is sensitive to the initial guess, often leading to divergence or convergence to spatially homogeneous solution. Systematically numerical exploration of spatial patterns of reaction-diffusion equations under different parameter regimes requires that the numerical method be efficient and robust to initial condition or initial guess, with better likelihood of convergence to an inhomogeneous pattern. Here, a new approach that combines the advantages of temporal schemes in robustness and Newton's method in fast convergence in solving steady states of reaction-diffusion equations is proposed. In particular, an adaptive implicit Euler with inexact solver (AIIE) method is found to be much more efficient than temporal schemes and more robust in convergence than typical nonlinear solvers (e.g., Newton's method) in finding the inhomogeneous pattern. Application of this new approach to two reaction-diffusion equations in one, two, and three spatial dimensions, along with direct comparisons to several other existing methods, demonstrates that AIIE is a more desirable method for searching inhomogeneous spatial patterns of reaction-diffusion equations in a large parameter space. (C) 2012 Elsevier Inc. All rights reserved.

【 授权许可】

Free