A tensor-based exponential integrator for diffusion--reaction equations in common curvilinear coordinates

In this paper, we study a tensor-based method for the numerical solution of a class of diffusion--reaction equations defined on spatial domains that admit common curvilinear coordinate representations. Typical examples in 2D include disks (polar coordinates), and in 3D balls or cylinders (spherical or cylindrical coordinates) as well as spheres for problems involving the Laplace--Beltrami operator. The proposed approach is based on a carefully chosen finite difference discretization of the Laplace operators that yields matrices with a structured representation as sums of Kronecker products. For the time integration, we introduce a novel split variant of the exponential Euler method that effectively handles the stiffness and avoids the severe time step size restriction of classical explicit methods. By exploiting the peculiar form of the obtained discretized operators and the chosen splitting strategy, we compute the needed action of the $\varphi_1$ matrix function through suitable tensor-matrix products in a $μ$-mode framework. We demonstrate the efficiency the approach on a wide range of physically relevant 2D and 3D examples of coupled diffusion--reaction systems generating Turing patterns with up to $10^6$ degrees of freedom.

Source: arXiv:2604.11558v1 - http://arxiv.org/abs/2604.11558v1 PDF: https://arxiv.org/pdf/2604.11558v1 Original Link: http://arxiv.org/abs/2604.11558v1