On the Spatiotemporal Pattern Formation in Nonlinear Coupled Reaction–Diffusion Systems

Nonlinear coupled reaction–diffusion (NCRD) systems have played a crucial role in the emergence of spatiotemporal patterns across various scientific and engineering domains. The NCRD systems considered in this study encompass various models, such as linear, Gray–Scott, Brusselator, isothermal chemical, and Schnakenberg, with the aim of capturing the spatiotemporal patterns they generate. These models cover a diverse range of intricate spatiotemporal patterns found in nature, including spots, spot replication, stripes, hexagons, and more. A mixed-type modal discontinuous Galerkin approach is employed for solving one- and two-dimensional NCRD systems. This approach introduces a mathematical formulation to handle the occurrence of second-order derivatives in diffusion terms. For spatial discretization, hierarchical modal basis functions premised on orthogonal scaled Legendre polynomials are used. Moreover, a novel reaction term treatment is proposed for the NCRD systems, demonstrating an intrinsic feature of the new DG scheme and preventing erroneous solutions due to extremely nonlinear reaction terms. The proposed approach reduces the NCRD systems into a framework of ordinary differential equations in time, which are addressed by an explicit third-order TVD Runge–Kutta algorithm. The spatiotemporal patterns generated with the present approach are comparable to those found in the literature. This approach can readily be expanded to handle large multi-dimensional problems that appear as model equations in developed biological and chemical applications.