Skip to main content Accessibility help

Computation of Two-Phase Biomembranes with Phase DependentMaterial Parameters Using Surface Finite Elements

  • Charles M. Elliott and Björn Stinner


The shapes of vesicles formed by lipid bilayers with phase separation are governed by a bending energy with phase dependent material parameters together with a line energy associated with the phase interfaces. We present a numerical method to approximate solutions to the Euler-Lagrange equations featuring triangulated surfaces, isoparametric quadratic surface finite elements and the phase field approach for the phase separation. Furthermore, the method involves an iterative solution scheme that is based on a relaxation dynamics coupling a geometric evolution equation for the membrane surface with a surface Allen-Cahn equation. Remeshing and grid adaptivity are discussed, and in various simulations the influence of several physical parameters is investigated.


Corresponding author



Hide All
[1]Andelman, D., Kawakatsu, T. and Kawasaki, K., Equilibrium shape of two-component unil-amellar membranes and vesicles, Europhys. Lett., 19 (1992), 5762.
[2]Barrett, J. W., Garcke, H. and Nürnberg, R., On the parametric finite element approximation of evolving hypersurfaces in R3, J. Comput. Phys., 227 (2008), 42814307.
[3]Barrett, J. W., Garcke, H. and Nürnberg, R., Parametric approximation of Willmore flow and related geometric evolution equations, SIAM J. Sci. Comput., 31 (2008), 225253.
[4]Barrett, J. W., Garcke, H. and Nürnberg, R., Parametric approximation of surface clusters driven by isotropic and anisotropic surface energies, Interf. Free Bound., 12 (2010), 187234.
[5]Baumgart, T., Das, S., Webb, W. and Jenkins, J., Membrane elasticity in giant vesicles with fluid phase coexistence, Biophys. J., 89 (2005), 10671084.
[6]Baumgart, T., Hess, S. and Webb, W., Imaging coexisting fluid domains in biomembrane models coupling curvature and line tension, Nature, 425 (2003), 821824.
[7]Blowey, J. F. and Elliott, C. M., The Cahn-Hilliard gradient theory for phase separation with non-smooth free energy part ii: Numerical analysis, Euro. J. Appl. Math., 3 (1993), 147179.
[8]Brenner, S. and Scott, L., The Mathematical Theory of Finite Element Methods, No. 15 in Texts in Applied Mathematics, Springer, Third ed., 2007.
[9]Campelo, F., Allain, J.-M. and Amar, M. B., Periodic lipidic membrane tubes, Europhys. Lett., 77 (2007), 38006.
[10]Canham, P., The minimum energy of bending as a possible explanation of the biconcave shape of the red blood cell, J. Theor. Biol., 26 (1970), 6181.
[11]Chen, L.-Q., Phase-field models for microstructure evolution, Ann. Rev. Mat. Res., 32 (2002), 113140.
[12]Das, S., Jenkins, J. and Baumgart, T., Neck geometry and shape transitions in vesicles with co-existing fluid phases: role of gaussian curvature stiffness vs. spontaneous curvature, Europhys. Lett., 86 (2009), 480031-6.
[13]Davies, T. A., Algorithm 832: Umfpack, an unsymmetric-pattern multifrontal method, ACM Trans. Math. Software, 30 (2004), 196199.
[14]Deckelnick, K., Dziuk, G. and Elliott, C. M., Computation of geometric partial differential equations and mean curvature flow, Acta Numer., 14 (2005), 139232.
[15]Demlow, A., Higher-order finite element methods and pointwise error estimates for elliptic problems on surfaces, SIAM J. Numer. Anal., 47 (2009), 805827.
[16]Du, Q., Liu, C. and Wang, X., A phase field approach in the numerical study of the elastic bending energy for vesicle membranes, J. Comput. Phys., 198 (2004), 450468.
[17]Du, Q., Liu, C. and Wang, X., Simulating the deformation of vesicle membranes under elastic bending energy in three dimensions, J. Comput. Phys., 212 (2006), 757777.
[18]Dziuk, G., Finite elements for the Beltrami operator on arbitrary surfaces, in Partial Differential Equations and Calculus of Variations, Hildebrandt, S. and Leis, R., eds., Vol. 1357 of Lecture Notes in Mathematics, Springer, 1988, 142155.
[19]Dziuk, G., An algorithm for evolutionary surfaces, Numerische Mathematik, 58 (1991), 603611.
[20]Dziuk, G., Computational parametric Willmore flow, Numerische Mathematik, 111 (2008), 5580.
[21]Dziuk, G. and Clarenz, U., Numerical methods for conformally parametrized surfaces, CPDw04-Interphase 2003: Numerical Methods for Free Boundary Problems, workshop at the Isaac Newton Institute, /cpdw04/0415/dziuk/, (2003).
[22]Dziuk, G. and Elliott, C. M., Finite elements on evolving surfaces, IMA J. Numer. Anal., 25 (2007), 385407.
[23]Dziuk, G. and Elliott, C. M., Surface finite elements for parabolic equations, J. Comput. Math., 25 (2007), 385407.
[24]Elliott, C. and Stinner, B., Modeling and computation of two phase geometric biomembranes using surface finite elements, J. Comput. Phys., 229 (2010), 65856612.
[25]Elliott, C. and Stinner, B., A surface phase field model for two-phase biological membranes, SIAM J. Appl. Math., 70 (2010), 29042928.
[26]Elliott, C. M. and Stuart, A. M., The global dynamics of discrete semilinear parabolic equations, SIAM J. Numer. Anal., 30 (1993), 16221663.
[27]Evans, E., Bending resistance and chemically induced moments in membrane bilayers, Biophys. J., 14 (1974), 923931.
[28]Feng, F. and Klug, W., Finite element modeling of liquid bilayer membranes, J. Comput. Phys., 220 (2006), 394408.
[29]Garcke, H., Nestler, B., Stinner, B. and Wendler, F., Allen-Cahn systems with volume constraints, Math. Mod. Meth. Appl. Sci., 18 (2008), 13471381.
[30]Heine, C.-J., Computations of Form and Stability of Rotating Drops with Finite Elements, PhD thesis, Faculty for Mathematics, Informatics and Natural Sciences, University of Aachen, 2003.
[31]Heine, C.-J., Isoparametric finite element approximation of curvature on hypersurfaces, Fakultät für Mathematik und Physik, University of Freiburg, Preprint, 26 (2004).
[32]Helfrich, W., Elastic properties of lipid bilayers: theory and possible experiments, Z. Natur-forschung, C28 (1973), 693703.
[33]Helmers, M., Snapping elastic curves as a one-dimensional analogue of two-component lipid bilayers, to appear in Math. Mod. Meth. Appl. Sci., DOI: 10.1142/S0218202511005234, (2011).
[34]Jülicher, F. and Lipowsky, R., Domain-induced budding of vesicles, Phys. Rev. Lett., 70 (1993), 29642967.
[35]Jülicher, F. and Lipowsky, R., Shape transformations of vesicles with intramembrane domains, Phys. Rev. E, 53 (1996), 26702683.
[36]Kawakatsu, T., Andelman, D., Kawasaki, K. and Taniguchi, T., Phase transitions and shapes of two component membranes and vesicles I: strong segregation limit, J. Phys. II (France), 3 (1993), 971997.
[37]Lowengrub, J. S., Rätz, A. and Voigt, A., Phase-field modeling of the dynamics of multicomponent vesicles: spinodal decomposition, coarsening, budding and fission, Phys. Rev. E, 79 (2009), 0319261-13.
[38]Ma, L. and Klug, W., Viscous regularization and r-adaptive remeshing for finite element analysis of lipid membrane mechanics, J. Comput. Phys., 227 (2008), 58165835.
[39]Schmidt, A. and Siebert, K. G., Design of adaptive finite element software: the finite element toolbox ALBERTA, Vol. 42 of Lecture Notes in Computational Science and Engineering, Springer, 2005.
[40]Seifert, U., Configurations of fluid membranes and vesicles, Adv. Phys., 46 (1997), 1137.
[41]Taniguchi, T., Shape deformations and phase separation dynamics of two-component vesicles, Phys. Rev. Lett., 76 (1996), 44444447.
[42]Taniguchi, T., Kawasaki, K., Andelman, D. and Kawakatsu, T., Phase transitions and shapes of two component membranes and vesicles II: weak segregation limit, J. Phys. II (France), 4 (1994), 13331362.
[43]Veatch, S. L. and Keller, S. L., Separation of liquid phases in giant vesicles of ternary mixtures of phospholipids and cholesterol, Biophys. J., 85 (2003), 30743083.
[44]Wang, X. and Du, Q., Modelling and simulations of multi-component lipid membranes and open membranes via diffuse interface approaches, J. Math. Biol., 56 (2008), 347371.


Related content

Powered by UNSILO

Computation of Two-Phase Biomembranes with Phase DependentMaterial Parameters Using Surface Finite Elements

  • Charles M. Elliott and Björn Stinner


Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed.