Finite element error analysis for a system coupling surface evolution to diffusion on the surface

We consider a numerical scheme for the approximation of a system that couples the evolution of a two-dimensional hypersurface to a reaction–diffusion equation on the surface. The surfaces are assumed to be graphs and evolve according to forced mean curvature flow. The method uses continuous, piecewise linear finite elements in space and a backward Euler scheme in time. Assuming the existence of a smooth solution, we prove optimal error bounds both in L∞(L2) and in L2(H1). We present several numerical experiments that confirm our theoretical findings and apply the method in order to simulate diffusion induced grain boundary motion