The surface finite element method for pattern formation on evolving biological surfaces

In this article we propose models and a numerical method for pattern formation on evolving curved surfaces. We formulate reaction-diffusion equations on evolving surfaces using the material transport formula, surface gradients and diffusive conservation laws. The evolution of the surface is defined by a material surface velocity. The numerical method is based on the evolving surface finite element method. The key idea is based on the approximation of Γ by a triangulated surface Γh consisting of a union of triangles with vertices on Γ. A finite element space of functions is then defined by taking the continuous functions on Γh which are linear affine on each simplex of the polygonal surface. To demonstrate the capability, flexibility, versatility and generality of our methodology we present results for uniform isotropic growth as well as anisotropic growth of the evolution surfaces and growth coupled to the solution of the reaction-diffusion system. The surface finite element method provides a robust numerical method for solving partial differential systems on continuously evolving domains and surfaces with numerous applications in developmental biology, tumour growth and cell movement and deformation.