A fully discrete evolving surface finite element method

In this paper we consider a time discrete evolving surface finite element method for the advection and diffusion of a conserved scalar quantity on a moving surface. In earlier papers using a suitable variational formulation in time dependent Sobolev space we proposed and analyzed a finite element method using surface finite elements on evolving triangulated surfaces [IMA J. Numer Anal., 25 (2007), pp. 385--407; Math. Comp., to appear]. Optimal order L2(Γ(t)) and H1(Γ(t)) error bounds were proved for linear finite elements. In this work we prove optimal order error bounds for a backward Euler time discretization

L² -estimates for the evolving surface finite element method

In this paper we consider the evolving surface finite element method for the advection and diffusion of a conserved scalar quantity on a moving surface. In an earlier paper using a suitable variational formulation in time dependent Sobolev space we proposed and analysed a finite element method using surface finite elements on evolving triangulated surfaces. An optimal order H¹ -error bound was proved for linear finite elements. In this work we prove the optimal error bound in L² (Γ(t)) uniformly in time

Finite element methods for surface PDEs

In this article we consider finite element methods for approximating the solution of partial differential equations on surfaces. We focus on surface finite elements on triangulated surfaces, implicit surface methods using level set descriptions of the surface, unfitted finite element methods and diffuse interface methods. In order to formulate the methods we present the necessary geometric analysis and, in the context of evolving surfaces, the necessary transport formulae. A wide variety of equations and applications are covered. Some ideas of the numerical analysis are presented along with illustrative numerical examples

Modelling cell motility and chemotaxis with evolving surface finite elements

We present a mathematical and a computational framework for the modelling of cell motility. The cell membrane is represented by an evolving surface, with the movement of the cell determined by the interaction of various forces that act normal to the surface. We consider external forces such as those that may arise owing to inhomogeneities in the medium and a pressure that constrains the enclosed volume, as well as internal forces that arise from the reaction of the cells' surface to stretching and bending. We also consider a protrusive force associated with a reaction-diffusion system (RDS) posed on the cell membrane, with cell polarization modelled by this surface RDS. The computational method is based on an evolving surface finite-element method. The general method can account for the large deformations that arise in cell motility and allows the simulation of cell migration in three dimensions. We illustrate applications of the proposed modelling framework and numerical method by reporting on numerical simulations of a model for eukaryotic chemotaxis and a model for the persistent movement of keratocytes in two and three space dimensions. Movies of the simulated cells can be obtained from http://homepages.warwick.ac.uk/maskae/CV_Warwick/Chemotaxis.html

Stability and error analysis for a diffuse interface approach to an advection-diffusion equation on a moving surface

In this paper we analyze a fully discrete numerical scheme for solving a parabolic PDE on a moving surface. The method is based on a diffuse interface approach that involves a level set description of the moving surface. Under suitable conditions on the spatial grid size, the time step and the interface width we obtain stability and error bounds with respect to natural norms. Furthermore, we present test calculations that confirm our analysis

Semilinear mixed problems on Hilbert complexes and their numerical approximation

Arnold, Falk, and Winther recently showed [Bull. Amer. Math. Soc. 47 (2010), 281-354] that linear, mixed variational problems, and their numerical approximation by mixed finite element methods, can be studied using the powerful, abstract language of Hilbert complexes. In another recent article [arXiv:1005.4455], we extended the Arnold-Falk-Winther framework by analyzing variational crimes (a la Strang) on Hilbert complexes. In particular, this gave a treatment of finite element exterior calculus on manifolds, generalizing techniques from surface finite element methods and recovering earlier a priori estimates for the Laplace-Beltrami operator on 2- and 3-surfaces, due to Dziuk [Lecture Notes in Math., vol. 1357 (1988), 142-155] and later Demlow [SIAM J. Numer. Anal., 47 (2009), 805-827], as special cases. In the present article, we extend the Hilbert complex framework in a second distinct direction: to the study of semilinear mixed problems. We do this, first, by introducing an operator-theoretic reformulation of the linear mixed problem, so that the semilinear problem can be expressed as an abstract Hammerstein equation. This allows us to obtain, for semilinear problems, a priori solution estimates and error estimates that reduce to the Arnold-Falk-Winther results in the linear case. We also consider the impact of variational crimes, extending the results of our previous article to these semilinear problems. As an immediate application, this new framework allows for mixed finite element methods to be applied to semilinear problems on surfaces.Comment: 22 pages; v2: major revision, particularly sharpening of error estimates in Section

Evolving surface finite element method for the Cahn-Hilliard equation

We use the evolving surface finite element method to solve a Cahn- Hilliard equation on an evolving surface with prescribed velocity. We start by deriving the equation using a conservation law and appropriate transport for- mulae and provide the necessary functional analytic setting. The finite element method relies on evolving an initial triangulation by moving the nodes according to the prescribed velocity. We go on to show a rigorous well-posedness result for the continuous equations by showing convergence, along a subse- quence, of the finite element scheme. We conclude the paper by deriving error estimates and present various numerical examples

Parameter identification problems in the modelling of cell motility

We present a novel parameter identification algorithm for the estimation of parameters in models of cell motility using imaging data of migrating cells. Two alternative formulations of the objective functional that measures the difference between the computed and observed data are proposed and the parameter identification problem is formulated as a minimisation problem of nonlinear least squares type. A Levenberg–Marquardt based optimisation method is applied to the solution of the minimisation problem and the details of the implementation are discussed. A number of numerical experiments are presented which illustrate the robustness of the algorithm to parameter identification in the presence of large deformations and noisy data and parameter identification in three dimensional models of cell motility. An application to experimental data is also presented in which we seek to identify parameters in a model for the monopolar growth of fission yeast cells using experimental imaging data. Our numerical tests allow us to compare the method with the two different formulations of the objective functional and we conclude that the results with both objective functionals seem to agree

High-order spectral/hp element discretisation for reaction-diffusion problems on surfaces: application to cardiac electrophysiology

We present a numerical discretisation of an embedded two-dimensional manifold using high-order continuous Galerkin spectral/hp elements, which provide exponential convergence of the solution with increasing polynomial order, while retaining geometric flexibility in the representation of the domain. Our work is motivated by applications in cardiac electrophysiology where sharp gradients in the solution benefit from the high-order discretisation, while the compu- tational cost of anatomically-realistic models can be reduced through the surface representation. We describe and validate our discretisation and provide a demonstration of its application to modeling electrochemical propagation across a human left atrium

A discrete Laplace-Beltrami operator for simplicial surfaces

We define a discrete Laplace-Beltrami operator for simplicial surfaces. It depends only on the intrinsic geometry of the surface and its edge weights are positive. Our Laplace operator is similar to the well known finite-elements Laplacian (the so called ``cotan formula'') except that it is based on the intrinsic Delaunay triangulation of the simplicial surface. This leads to new definitions of discrete harmonic functions, discrete mean curvature, and discrete minimal surfaces. The definition of the discrete Laplace-Beltrami operator depends on the existence and uniqueness of Delaunay tessellations in piecewise flat surfaces. While the existence is known, we prove the uniqueness. Using Rippa's Theorem we show that, as claimed, Musin's harmonic index provides an optimality criterion for Delaunay triangulations, and this can be used to prove that the edge flipping algorithm terminates also in the setting of piecewise flat surfaces.Comment: 18 pages, 6 vector graphics figures. v2: Section 2 on Delaunay triangulations of piecewise flat surfaces revised and expanded. References added. Some minor changes, typos corrected. v3: fixed inaccuracies in discussion of flip algorithm, corrected attributions, added references, some minor revision to improve expositio