Finite element approximation for the dynamics of asymmetric fluidic biomembranes

We present a parametric finite element approximation of a fluidic membrane, whose evolution is governed by a surface Navier–Stokes equation coupled to bulk Navier–Stokes equations. The elastic properties of the membrane are modelled with the help of curvature energies of Willmore and Helfrich type. Forces stemming from these energies act on the surface fluid, together with a forcing from the bulk fluid. Using ideas from PDE constrained optimization, a weak formulation is derived, which allows for a stable semi-discretization. An important new feature of the present work is that we are able to also deal with spontaneous curvature and an area-difference elasticity contribution in the curvature energy. This allows for the modelling of asymmetric membranes, which compared to the symmetric case lead to quite different shapes. This is demonstrated in the numerical computations presented

Isogeometric Analysis for High Order Geometric Partial Differential Equations with Applications

In this thesis, we consider the numerical approximation of high order geometric Partial Differential Equations (PDEs). We first consider high order PDEs defined on surfaces in the 3D space that are represented by single-patch tensor product NURBS. Then, we spatially discretize the PDEs by means of NURBS-based Isogeometric Analysis (IGA) in the framework of the Galerkin method. With this aim, we consider the construction of periodic NURBS function spaces with high degree of global continuity, even on closed surfaces. As benchmark problems for the proposed discretization, we propose Laplace-Beltrami problems of the fourth and sixth orders, as well as the corresponding eigenvalue problems, and we analyze the impact of the continuity of the basis functions on the accuracy as well as on computational costs. The numerical solution of two high order phase field problems on both open and closed surfaces is also considered: the fourth order Cahn-Hilliard equation and the sixth order crystal equation, both discretized in time with the generalized-alpha method. We then consider the numerical approximation of geometric PDEs, derived, in particular, from the minimization of shape energy functionals by L^2-gradient flows. We analyze the mean curvature and the Willmore gradient flows, leading to second and fourth order PDEs, respectively. These nonlinear geometric PDEs are discretized in time with Backward Differentiation Formulas (BDF), with a semi-implicit formulation based on an extrapolation of the geometry, leading to a linear problem to be solved at each time step. Results about the numerical approximation of the two geometric flows on several geometries are analyzed. Then, we study how the proposed mathematical framework can be employed to numerically approximate the equilibrium shapes of lipid bilayer biomembranes, or vesicles, governed by the Canham-Helfrich curvature model. We propose two numerical schemes for enforcing the conservation of the area and volume of the vesicles, and report results on benchmark problems. Then, the approximation of the equilibrium shapes of biomembranes with different values of reduced volume is presented. Finally, we consider the dynamics of a vesicle, e.g. a red blood cell, immersed in a fluid, e.g. the plasma. In particular, we couple the curvature-driven model for the lipid membrane with the incompressible Navier-Stokes equations governing the fluid. We consider a segregated approach, with a formulation based on the Resistive Immersed Surface method applied to NURBS geometries. After analyzing benchmark fluid simulations with immersed NURBS objects, we report numerical results for the investigation of the dynamics of a vesicle under different flow conditions

An adaptive finite element method for the modeling of the equilibrium of red blood cells

International audienceThis contribution is concerned with a the numerical modeling of an isolated red blood cell (RBC), and more generally of phospholipid membranes. We propose an adaptive Eulerian finite element approximation, based on the level set method, of a shape optimization problem arising in the study of RBC's equilibrium. We simulate the equilibrium shapes that minimize the elastic bending energy under prescribed constraints of fixed volume and surface area. An anisotropic mesh adaptation technique is used in the vicinity of the cell's membrane to enhance the robustness of the method. Efficient time and spatial discretizations are considered and implemented. We address in detail the main features of the proposed method and finally we report several numerical experiments in the two-dimensional and the three-dimensional axisymmetric cases. The effectiveness of the numerical method is further demonstrated through numerical comparisons with semi-analytical solutions provided by a reduced order model

Numerical simulation of the dynamics of fluid membranes

The goal of this project is to explore the mechanics of lipid fluid membranes as found in biological and made-man systems through continuum models and numerical simulations. Traditionally, the focus has been on the equilibrium configurations of vesicles through minimization of the curvature energy subject to constraints. Here, the goal is to describe the time-evolution of out-of-equilibrium vesicle configurations, which are of relevance in biological systems. Towards this goal, an accurate description of the dissipative mechanisms is crucial, in particular the viscous dissipation induced by the 2D flow of lipids on the deforming surface that describes a vesicle [1]. The resulting equations can only be solved analytically in very simple settings. So the problem is solved numerically using a B-Spline description of the membrane vesicle. The dynamics are described using two types of viscosity: the L2 or Willmore viscosity and the inner flow viscosity. By comparing the evolution in time of the two systems, it is stated that the dynamics are clearly different. It is found that one describes better the physics of the system

Modelling fluid deformable surfaces with an emphasis on biological interfaces

This article has been published in a revised form in Journal of fluid mechanics, http://dx.doi.org/10.1017/jfm.2019.341. This version is free to view and download for private research and study only. Not for re-distribution, re-sale or use in derivative works. © 2019Fluid deformable surfaces are ubiquitous in cell and tissue biology, including lipid bilayers, the actomyosin cortex or epithelial cell sheets. These interfaces exhibit a complex interplay between elasticity, low Reynolds number interfacial hydrodynamics, chemistry and geometry, and govern important biological processes such as cellular traffic, division, migration or tissue morphogenesis. To address the modelling challenges posed by this class of problems, in which interfacial phenomena tightly interact with the shape and dynamics of the surface, we develop a general continuum mechanics and computational framework for fluid deformable surfaces. The dual solid–fluid nature of fluid deformable surfaces challenges classical Lagrangian or Eulerian descriptions of deforming bodies. Here, we extend the notion of arbitrarily Lagrangian–Eulerian (ALE) formulations, well-established for bulk media, to deforming surfaces. To systematically develop models for fluid deformable surfaces, which consistently treat all couplings between fields and geometry, we follow a nonlinear Onsager formalism according to which the dynamics minimizes a Rayleighian functional where dissipation, power input and energy release rate compete. Finally, we propose new computational methods, which build on Onsager’s formalism and our ALE formulation, to deal with the resulting stiff system of higher-order partial differential equations. We apply our theoretical and computational methodology to classical models for lipid bilayers and the cell cortex. The methods developed here allow us to formulate/simulate these models in their full three-dimensional generality, accounting for finite curvatures and finite shape changes.Peer ReviewedPostprint (author's final draft