Optimized high-order IMEX-RK schemes for degenerate diffusion-reaction problems with application to travelling waves phenomena
For computational scientists modeling biological and physical systems with travelling waves, this work provides a flexible and robust numerical framework, though it is an incremental improvement over existing IMEX-RK methods.
The paper develops high-order IMEX-RK schemes for degenerate diffusion-reaction problems, achieving accurate capture of sharp propagating fronts in travelling wave phenomena. The methods are validated on a hierarchy of test problems, demonstrating robust performance.
We study a class of IMplicit-EXplicit Runge--Kutta (IMEX-RK) schemes for the numerical approximation of reaction and diffusion-reaction problems arising in a variety of biological and physical applications. Such models may admit travelling wave solutions, with the Fisher--Kolmogorov equation representing a prototypical example. Motivated by this feature, the proposed time integration schemes are designed to accurately capture sharp propagating fronts. We also investigate a less standard use of IMEX-RK methods that circumvents a splitting of reaction terms into linear and nonlinear components, while still requiring the solution of linear systems at each stage. This semi-implicit formulation, referred to as SI-IMEX-RK, enables a targeted treatment of stiffness by isolating its relevant contributions. The time discretization is coupled with a high-order polygonal discontinuous Galerkin method for space discretization, resulting in a flexible and robust framework for the treatment of multiscale dynamics in complex geometries. A comprehensive validation strategy is presented to assess the accuracy and stability properties of the proposed schemes across a hierarchy of increasingly challenging test problems.