JournalsifbVol. 10, No. 1pp. 119–138

Eulerian finite element method for parabolic PDEs on implicit surfaces

  • Gerhard Dziuk

    Universität Freiburg, Germany
  • Charles M. Elliott

    University of Warwick, Coventry, UK
Eulerian finite element method for parabolic PDEs on implicit surfaces cover
Download PDF

Abstract

We define an Eulerian level set method for parabolic partial differential equations on a stationary hypersurface Γ\Gamma contained in a domain ΩRn+1\Omega \subset \mathbb R^{n+1}. The method is based on formulating the partial differential equations on all level surfaces of a prescribed function Φ\Phi whose zero level set is Γ\Gamma. Eulerian surface gradients are formulated by using a projection of the gradient in Rn+1\mathbb R^{n+1} onto the level surfaces of Φ\Phi. These Eulerian surface gradients are used to define weak forms of surface elliptic operators and so generate weak formulations of surface elliptic and parabolic equations. The resulting equation is then solved in one dimension higher but can be solved on a mesh which is unaligned to the level sets of Φ\Phi. We consider both second order and fourth order elliptic operators with natural second order splittings. The finite element method is applied to the weak form of the split system of second order equations using piece-wise linear elements on a fixed grid. The computation of the mass and element stiffness matrices are simple and straightforward. Numerical experiments are described which indicate the power of the method. We describe how this framework may be employed in applications.

Cite this article

Gerhard Dziuk, Charles M. Elliott, Eulerian finite element method for parabolic PDEs on implicit surfaces. Interfaces Free Bound. 10 (2008), no. 1, pp. 119–138

DOI 10.4171/IFB/182