MATHICSE Technical Report : Biomembrane modeling with Isogeometric Analysis

Andrea Bartezzaghi, Alfio Quarteroni
MATHICSE, 2017
Rapport ou document de travail

Résumé

We consider the numerical approximation of lipid biomembranes, including red blood cells, described through the Canham-Helfrich model, according to which the shape minimizes the bending energy under area and volume constraints. Energy minimization is performed via L2- gradient flow of the Canham-Helfrich energy using two Lagrange multipliers to weakly enforce the constraints. This yields a highly nonlinear, high order, time dependent geometric Partial Differential Equation (PDE). We represent the biomembranes as single-patch NURBS closed surfaces. We discretize the geometric PDEs in space with NURBS-based Isogeometric Analysis and in time with Backward Differentiation Formulas. We tackle the nonlinearity in our formulation through a semi-implicit approach by extrapolating, at each time level, the geometric quantities of interest from previous time steps. We report the numerical results of the approximation of the Canham-Helfrich problem on ellipsoids of different aspect ratio, which lead to the classical biconcave shape of lipid vesicles at equilibrium. We show that this framework permits an accurate approximation of the Canham-Helfrich problem, while being computationally efficient.

Official source

https://infoscience.epfl.ch/record/267616?ln=fr

À propos de ce résultat

Cette page est générée automatiquement et peut contenir des informations qui ne sont pas correctes, complètes, à jour ou pertinentes par rapport à votre recherche. Il en va de même pour toutes les autres pages de ce site. Veillez à vérifier les informations auprès des sources officielles de l'EPFL.

Concepts associés

Chargement

Publications associées

Chargement

Biomembrane modeling with isogeometric analysis

Andrea Bartezzaghi, Alfio Quarteroni

We consider the numerical approximation of lipid biomembranes at equilibrium described by the Canham-Helfrich model, according to which the bending energy is minimized under area and volume constraints. Energy minimization is performed via L-2-gradient flow of the Canham-Helfrich energy using two Lagrange multipliers to weakly enforce the constraints. This yields a highly nonlinear, high order, time dependent geometric Partial Differential Equation (PDE). We represent the biomembranes as single-patch NURBS closed surfaces. We discretize the geometric PDEs in space with NURBS-based Isogeometric Analysis and in time with Backward Differentiation Formulas. We tackle the nonlinearity in our formulation through a semi-implicit approach by extrapolating, at each time level, the geometric quantities of interest from previous time steps. We report the numerical results of the approximation of the Canham-Helfrich problem on ellipsoids of different aspect ratio, which leads to the classical biconcave shape of lipid vesicles at equilibrium. We show that this framework permits an accurate approximation of the Canham-Helfrich problem, while being computationally efficient. (C) 2019 Elsevier B.Y. All rights reserved.

ELSEVIER SCIENCE SA2019

Chargement

Isogeometric Analysis of geometric partial differential equations

Andrea Bartezzaghi, Alfio Quarteroni

We consider the numerical approximation of geometric Partial Differential Equations (PDEs) defined on surfaces in the 3D space. In particular, we focus on the geometric PDEs deriving from the minimization of an energy functional by L2L2-gradient flow. We analyze two energy functionals: the area, which leads to the mean curvature flow, a nonlinear second order PDE, and the Willmore energy, leading to the Willmore flow, a nonlinear fourth order PDE. We consider surfaces represented by single-patch tensor product NURBS and discretize the PDEs by means of NURBS-based Isogeometric Analysis in the framework of the Galerkin method. To approximate the high order geometric PDEs we use high order continuous NURBS basis functions. For the time discretization of the nonlinear geometric PDEs, we use Backward Differentiation Formulas (BDF) with extrapolation of the geometric quantities involved in the weak formulation of the problem; in this manner, we solve a linear problem at each time step. We report numerical results concerning the mean curvature and Willmore flows on different geometries of interest and we show the accuracy and efficiency of the proposed approximation scheme.

Elsevier2016

Chargement

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.

EPFL2017

Chargement

EPFL2017