A weighted polygonal discontinuous Galerkin method for hierarchically coupled reaction-diffusion systems with cubic interactions
Abstract
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 L2(Ω)-and dG-norms. We further establish an a priori error estimate in a combined L2(Ω)-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.
Related
- From fleet trajectories to single-cell promises: hierarchically weighted conformal prediction for battery SOH monitoring
- Preconditioning for the time-dependent Stokes equations: stage-weighted Radau IIA and weighted discontinuous Galerkin time-stepping
- A super-convergent quasi-discontinuous Galerkin method for advection-diffusion problems
- A Fully Discrete Local Discontinuous Galerkin Method for Quasilinear Stochastic Convection-Diffusion-Type Equations
- Transport-Matched Penalties for Diffusion Synthetic Acceleration of Polytopic Discontinuous Galerkin Discretisations
BibTeX
@misc{corti2026weighted,
title = {{A weighted polygonal discontinuous Galerkin method for hierarchically coupled reaction-diffusion systems with cubic interactions}},
author = {Mattia Corti},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.34625},
eprint = {2609.34625},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.34625},
}
RIS
TY - UNPB
TI - A weighted polygonal discontinuous Galerkin method for hierarchically coupled reaction-diffusion systems with cubic interactions
AU - Mattia Corti
PY - 2026
DA - 2026/09/28
PB - arXiv
DO - 10.48550/arxiv.2609.34625
UR - https://arxiv.org/abs/2609.34625
AB - 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.
ER -