Track new results across pure and applied mathematics.
Search papers, labs, and topics across Lattice.
The authors formulate a symmetric weighted interior penalty polytopal discontinuous Galerkin (SWIP-PolyDG) framework for multi-species reaction-diffusion systems governed by non-additive cubic nonlinearities and anisotropic, heterogeneous diffusion on arbitrary 2D polygonal meshes. Handling cubic interactions over complex geometries typically risks catastrophic numerical instabilities unless strong structural dissipativity is preserved at the discrete level. By enforcing a block-hierarchical reaction structure with dissipative diagonal self-interactions, the authors establish recursive local-in-time stability in $L^2$- and dG-norms alongside optimal a priori error bounds in a combined energy norm.
Cubic reaction instabilities on irregular polygonal meshes can be completely subdued by pairing polytopal interior penalty discontinuous Galerkin schemes with block-hierarchical dissipative reaction structures.
We develop and analyse a symmetric weighted interior penalty polytopal discontinuous Galerkin (SWIP-PolyDG) method for multi-species reaction--diffusion systems with nonlinear reaction terms of up to cubic order. Such terms arise naturally when higher-order interactions are incorporated into population models, allowing non-additive effects among multiple species to influence local growth and conversion mechanisms. The proposed framework accommodates heterogeneous and possibly anisotropic diffusion tensors on general polygonal meshes. To control the nonlinear coupling, we consider a hierarchical block structure in the reaction operator, whereby each population block depends only on its own variables and on those associated with preceding blocks. In addition, the cubic self-interactions within each block are assumed to have a dissipative diagonal structure. Under these hypotheses and proceeding recursively over the hierarchy of population blocks, we derive local-in-time stability estimates in two spatial dimensions for the semi-discrete formulation in both the $L^2(Ω)$-and dG-norms. We further establish an a priori error estimate in a combined $L^2(Ω)$-dG energy norm for sufficiently regular solutions for the semi-discrete formulation. Finally, numerical experiments confirm the predicted convergence behaviour and illustrate the robustness of the method under heterogeneous diffusion.