Domain-Dependent Stability Analysis of a Reaction–Diffusion Model on Compact Circular Geometries

Abstract

In this work an activator-depleted reaction–diffusion system is investigated on polar coordinates with the aim of exploring the relationship and the corresponding influence of domain size on the types of possible diffusion-driven instabilities. Quantitative relationships are found in the form of necessary conditions on the area of a disk-shape domain with respect to the diffusion and reaction rates for certain types of diffusion-driven instabilities to occur. Robust analytical methods are applied to find explicit expressions for the eigenvalues and eigenfunctions of the diffusion operator on a disk-shape domain with homogenous Neumann boundary conditions in polar coordinates. Spectral methods are applied using Chebyshev nonperiodic grid for the radial variable and Fourier periodic grid for the angular variable to verify the nodal lines and eigen-surfaces subject to the proposed analytical findings. The full classification of the parameter space in light of the bifurcation analysis is obtained and numerically verified by finding the solutions of the partitioning curves inducing such a classification. Spatiotemporal periodic behavior is demonstrated in the numerical solutions of the system for a proposed choice of parameters and a rigorous proof of the existence of infinitely many such points in the parameter plane is presented under a restriction on the area of the domain, with a lower bound in terms of reaction–diffusion rates.