Gradient flow dynamics of two-phase biomembranes: Sharp interface variational formulation and finite element approximation
; ;
The SMAI journal of computational mathematics, Volume 4 (2018) , pp. 151-195.

A finite element method for the evolution of a two-phase membrane in a sharp interface formulation is introduced. The evolution equations are given as an L 2 –gradient flow of an energy involving an elastic bending energy and a line energy. In the two phases Helfrich-type evolution equations are prescribed, and on the interface, an evolving curve on an evolving surface, highly nonlinear boundary conditions have to hold. Here we consider both C 0 – and C 1 –matching conditions for the surface at the interface. A new weak formulation is introduced, allowing for a stable semidiscrete parametric finite element approximation of the governing equations. In addition, we show existence and uniqueness for a fully discrete version of the scheme. Numerical simulations demonstrate that the approach can deal with a multitude of geometries. In particular, the paper shows the first computations based on a sharp interface description, which are not restricted to the axisymmetric case.

Published online: 2018-04-27
DOI: https://doi.org/10.5802/smai-jcm.32
Classification: 35R01,  49Q10,  65M12,  65M60,  82B26,  92C10
Keywords: parametric finite elements, Helfrich energy, spontaneous curvature, multi-phase membrane, line energy, C 0 – and C 1 –matching conditions
@article{SMAI-JCM_2018__4__151_0,
     author = {John W. Barrett and Harald Garcke and Robert N\"urnberg},
     title = {Gradient flow dynamics of two-phase biomembranes: Sharp interface variational formulation and finite element approximation},
     journal = {The SMAI journal of computational mathematics},
     publisher = {Soci\'et\'e de Math\'ematiques Appliqu\'ees et Industrielles},
     volume = {4},
     year = {2018},
     pages = {151-195},
     doi = {10.5802/smai-jcm.32},
     language = {en},
     url={smai-jcm.centre-mersenne.org/item/SMAI-JCM_2018__4__151_0/}
}
Barrett, John W.; Garcke, Harald; Nürnberg, Robert. Gradient flow dynamics of two-phase biomembranes: Sharp interface variational formulation and finite element approximation. The SMAI journal of computational mathematics, Volume 4 (2018) , pp. 151-195. doi : 10.5802/smai-jcm.32. https://smai-jcm.centre-mersenne.org/item/SMAI-JCM_2018__4__151_0/

[1] H. Abels; H. Garcke; L. Müller Local well-posedness for volume-preserving mean curvature and Willmore flows with line tension, Math. Nachr., Volume 289 (2016) no. 2–3, pp. 136-174 | Article

[2] J. W. Barrett; H. Garcke; R. Nürnberg A parametric finite element method for fourth order geometric evolution equations, J. Comput. Phys., Volume 222 (2007) no. 1, pp. 441-462 | Article | MR 2298053

[3] J. W. Barrett; H. Garcke; R. Nürnberg On the parametric finite element approximation of evolving hypersurfaces in 3 , J. Comput. Phys., Volume 227 (2008) no. 9, pp. 4281-4307 | Article | MR 2406538

[4] J. W. Barrett; H. Garcke; R. Nürnberg Parametric Approximation of Surface Clusters driven by Isotropic and Anisotropic Surface Energies, Interfaces Free Bound., Volume 12 (2010) no. 2, pp. 187-234 | Article

[5] J. W. Barrett; H. Garcke; R. Nürnberg The Approximation of Planar Curve Evolutions by Stable Fully Implicit Finite Element Schemes that Equidistribute, Numer. Methods Partial Differential Equations, Volume 27 (2011) no. 1, pp. 1-30 | Article

[6] J. W. Barrett; H. Garcke; R. Nürnberg Elastic flow with junctions: Variational approximation and applications to nonlinear splines, Math. Models Methods Appl. Sci., Volume 22 (2012) no. 11, 1250037 pages | Article

[7] J. W. Barrett; H. Garcke; R. Nürnberg On the Stable Numerical Approximation of Two-Phase Flow with Insoluble Surfactant, M2AN Math. Model. Numer. Anal., Volume 49 (2015) no. 2, pp. 421-458 | Article | MR 3342212

[8] J. W. Barrett; H. Garcke; R. Nürnberg Computational Parametric Willmore Flow with Spontaneous Curvature and Area Difference Elasticity effects, SIAM J. Numer. Anal., Volume 54 (2016) no. 3, pp. 1732-1762 | Article

[9] J. W. Barrett; H. Garcke; R. Nürnberg Finite Element Approximation for the Dynamics of Fluidic Two-Phase Biomembranes, M2AN Math. Model. Numer. Anal., Volume 51 (2017) no. 6, pp. 2319-2366 | Article

[10] J. W. Barrett; H. Garcke; R. Nürnberg Stable variational approximations of boundary value problems for Willmore flow with Gaussian curvature, IMA J. Numer. Anal., Volume 37 (2017) no. 4, pp. 1657-1709 | Article

[11] T. Baumgart; S. Das; W. W. Webb; J.T. Jenkins Membrane Elasticity in Giant Vesicles with Fluid Phase Coexistence, Biophys. J., Volume 89 (2005) no. 2, pp. 1067-1080 | Article

[12] T. Baumgart; Hess; S. T.; W. W. Webb Imaging coexisting fluid domains in biomembrane models coupling curvature and line tension, Nature, Volume 425 (2003) no. 6960, pp. 821-824 | Article

[13] R. Choksi; M. Morandotti; M. Veneroni Global minimizers for axisymmetric multiphase membranes, ESAIM Control Optim. Calc. Var., Volume 19 (2013) no. 4, pp. 1014-1029 | Article | MR 3182678

[14] G. Cox; J. Lowengrub The effect of spontaneous curvature on a two-phase vesicle, Nonlinearity, Volume 28 (2015) no. 3, pp. 773-793 | Article | MR 3319383

[15] S.L. Das; J.T. Jenkins; T. Baumgart Neck geometry and shape transitions in vesicles with co-existing fluid phases: Role of Gaussian curvature stiffness vs. spontaneous curvature, Europhys. Lett., Volume 86 (2009) no. 4, 48003 pages

[16] T. A. Davis Algorithm 832: UMFPACK V4.3—an unsymmetric-pattern multifrontal method, ACM Trans. Math. Software, Volume 30 (2004) no. 2, pp. 196-199 | Article | MR 2075981

[17] K. Deckelnick; G. Dziuk; C. M. Elliott Computation of geometric partial differential equations and mean curvature flow, Acta Numer., Volume 14 (2005), pp. 139-232 | Article | MR 2168343

[18] K. Deckelnick; H.-C. Grunau; M. Röger Minimising a relaxed Willmore functional for graphs subject to boundary conditions, Interfaces Free Bound., Volume 19 (2017) no. 1, pp. 109-140

[19] G. Dziuk An algorithm for evolutionary surfaces, Numer. Math., Volume 58 (1991) no. 6, pp. 603-611 | Article | MR 1083523

[20] G. Dziuk Computational parametric Willmore flow, Numer. Math., Volume 111 (2008) no. 1, pp. 55-80 | Article | MR 2448203

[21] G. Dziuk; C. M. Elliott Finite element methods for surface PDEs, Acta Numer., Volume 22 (2013), pp. 289-396 | MR 3038698

[22] C. M. Elliott; B. Stinner Modeling and computation of two phase geometric biomembranes using surface finite elements, J. Comput. Phys., Volume 229 (2010) no. 18, pp. 6585-6612 | Article | MR 2660322

[23] C. M. Elliott; B. Stinner A surface phase field model for two-phase biological membranes, SIAM J. Appl. Math., Volume 70 (2010) no. 8, pp. 2904-2928 | Article | MR 2735109

[24] C.M. Elliott; B. Stinner Computation of two-phase biomembranes with phase dependent material parameters using surface finite elements, Commun. Comput. Phys., Volume 13 (2013) no. 2, pp. 325-360 | MR 2948020

[25] M. Helmers Snapping elastic curves as a one-dimensional analogue of two-component lipid bilayers, Math. Models Methods Appl. Sci., Volume 21 (2011) no. 5, pp. 1027-1042

[26] M. Helmers Kinks in two-phase lipid bilayer membranes, Calc. Var. Partial Differential Equations, Volume 48 (2013) no. 1-2, pp. 211-242 | Article | MR 3090540

[27] M. Helmers Convergence of an approximation for rotationally symmetric two-phase lipid bilayer membranes, Q. J. Math., Volume 66 (2015) no. 1, pp. 143-170 | Article | MR 3356284

[28] F. Jülicher; R. Lipowsky Domain-induced budding of vesicles, Phys. Rev. Lett., Volume 70 (1993) no. 19, pp. 2964-2967

[29] F. Jülicher; R. Lipowsky Shape transformations of vesicles with intramembrane domains, Phys. Rev. E, Volume 53 (1996) no. 3, pp. 2670-2683 | Article

[30] J. S. Lowengrub; A. Rätz; A. Voigt Phase-field modeling of the dynamics of multicomponent vesicles: Spinodal decomposition, coarsening, budding, and fission, Phys. Rev. E, Volume 79 (2009) no. 3, 0311926 pages | MR 2497179

[31] M. Mercker; A. Marciniak-Czochra Bud-Neck Scaffolding as a Possible Driving Force in ESCRT-Induced Membrane Budding, Biophys. J., Volume 108 (2015) no. 4, pp. 833-843 | Article

[32] M. Mercker; A. Marciniak-Czochra; T. Richter; D. Hartmann Modeling and computing of deformation dynamics of inhomogeneous biological surfaces, SIAM J. Appl. Math., Volume 73 (2013) no. 5, pp. 1768-1792 | Article | MR 3095724

[33] J. C. C. Nitsche Boundary value problems for variational integrals involving surface curvatures, Quart. Appl. Math., Volume 51 (1993) no. 2, pp. 363-387 | MR 1218374

[34] A. Schmidt; K. G. Siebert Design of Adaptive Finite Element Software: The Finite Element Toolbox ALBERTA, Lecture Notes in Computational Science and Engineering, Volume 42, Springer-Verlag, Berlin, 2005 | MR 2127659

[35] S. Schmidt; V. Schulz Shape derivatives for general objective functions and the incompressible Navier–Stokes equations, Control Cybernet., Volume 39 (2010) no. 3, pp. 677-713 | MR 2791367

[36] M. E. Taylor Partial differential equations I. Basic theory, Applied Mathematical Sciences, Volume 115, Springer, 2011

[37] F. Tröltzsch Optimal Control of Partial Differential Equations: Theory, Methods and Applications, Graduate Studies in Mathematics, Volume 112, American Mathematical Society, Providence, RI, 2010, xvi+399 pages | MR 2583281

[38] Z.-C. Tu Challenges in theoretical investigations of configurations of lipid membranes, Chin. Phys. B, Volume 22 (2013) no. 2, 28701 pages | Article

[39] Z. C. Tu; Z. C. Ou-Yang A geometric theory on the elasticity of bio-membranes, J. Phys. A, Volume 37 (2004) no. 47, pp. 11407-11429 | Article | MR 2109079

[40] X. Wang; Q. Du Modelling and simulations of multi-component lipid membranes and open membranes via diffuse interface approaches, J. Math. Biol., Volume 56 (2008) no. 3, pp. 347-371 | Article | MR 2358438

[41] C. Wutz Variationsprobleme für elastische Biomembranen unter Berücksichtigung von Linienenergien (2010) (Diploma thesis)