A Spectral-Based Method for Network-Formation PDEs

Abstract

We propose and study a simple and scalable Fourier-based spectral method for a continuum model of network formation under periodic boundary conditions. The model provides the evolution of the pressure $p$ and the conductivity $m$ over time. The evolution of $p$ is given by an anisotropic Poisson equation, while the equation for $m$ contains three terms corresponding to a diffusion and an activation term of the network -- that depends on the gradient of the pressure -- as well as a relaxation term that acts as a decaying term. This system arises as a formal $L^2$-gradient flow of a non-convex energy functional. Our algorithm combines two ingredients: (i) a splitting method for the equation for $m$, where the activation and relaxation parts are solved analytically, and the diffusion part is solved via Fast Fourier Transform (FFT), and (ii) an FFT combined with the Conjugate Gradient (CG) method applied to the equation for the pressure. This makes the scheme easy to implement compared to implicit schemes and naturally extensible to three dimensions on uniform periodic grids. To showcase the method, we recover the previously documented influence of the activation strength $c$, the diffusion coefficient $D$, and the metabolic exponent $γ$ on the morphology of emergent networks, and report grid convergence results.