Arbitrary Lagrangian–Eulerian finite element method for lipid membranes

An arbitrary Lagrangian–Eulerian finite element method and numerical implementation for curved and deforming lipid membranes is presented here. The membrane surface is endowed with a mesh whose in-plane motion need not depend on the in-plane flow of lipids. Instead, in-plane mesh dynamics can be specified arbitrarily. A new class of mesh motions is introduced, where the mesh velocity satisfies the dynamical equations of a user-specified two-dimensional material. A Lagrange multiplier constrains the out-of-plane membrane and mesh velocities to be equal, such that the mesh and material always overlap. An associated numerical inf–sup instability ensues, and is removed by adapting established techniques in the finite element analysis of fluids. In our implementation, the aforementioned Lagrange multiplier is projected onto a discontinuous space of piecewise linear functions. The new mesh motion is compared to established Lagrangian and Eulerian formulations by investigating a pre-eminent numerical benchmark of biological significance: the pulling of a membrane tether from a flat patch and its subsequent lateral translation.