Home A structure preserving front tracking finite element method for the Mullins–Sekerka problem
Article
Licensed
Unlicensed Requires Authentication

A structure preserving front tracking finite element method for the Mullins–Sekerka problem

  • Robert Nürnberg EMAIL logo
Published/Copyright: March 23, 2022

Abstract

We introduce and analyse a fully discrete approximation for a mathematical model for the solidification and liquidation of materials of negligible specific heat. The model is a two-sided Mullins–Sekerka problem. The discretization uses finite elements in space and an independent parameterization of the moving free boundary. We prove unconditional stability and exact volume conservation for the introduced scheme. Several numerical simulations, including for nearly crystalline surface energies, demonstrate the practicality and accuracy of the presented numerical method.

JEL Classification: 35K55; 35R35; 65M12; 65M50; 65M60; 74E10; 74E15; 80A22

References

[1] N. D. Alikakos, P. W. Bates, and X. Chen, Convergence of the Cahn–Hilliard equation to the Hele-Shaw model, Arch. Rational Mech. Anal. 128 (1994), 165–205.10.1007/BF00375025Search in Google Scholar

[2] F. Almgren and J. E. Taylor, Flat flow is motion by crystalline curvature for curves with crystalline energies, J. Differential Geom. 42 (1995), 1–22.10.4310/jdg/1214457030Search in Google Scholar

[3] W. Bao and Q. Zhao, A structure-preserving parametric finite element method for surface diffusion, SIAM J. Numer. Anal. 59 (2021), 2775–2799.10.1137/21M1406751Search in Google Scholar

[4] J. W. Barrett, H. Garcke, and R. Nürnberg, Numerical approximation of anisotropic geometric evolution equations in the plane, IMA J. Numer. Anal. 28 (2008), 292–330.10.1093/imanum/drm013Search in Google Scholar

[5] J. W. Barrett, H. Garcke, and R. Nürnberg, On the parametric finite element approximation of evolving hypersurfaces in ℝ3, J. Comput. Phys. 227 (2008), 4281–4307.10.1016/j.jcp.2007.11.023Search in Google Scholar

[6] J. W. Barrett, H. Garcke, and R. Nürnberg, A variational formulation of anisotropic geometric evolution equations in higher dimensions, Numer. Math. 109 (2008), 1–44.10.1007/s00211-007-0135-5Search in Google Scholar

[7] J. W. Barrett, H. Garcke, and R. Nürnberg, On stable parametric finite element methods for the Stefan problem and the Mullins–Sekerka problem with applications to dendritic growth, J. Comput. Phys. 229 (2010), 6270–6299.10.1016/j.jcp.2010.04.039Search in Google Scholar

[8] J. W. Barrett, H. Garcke, and R. Nürnberg, Eliminating spurious velocities with a stable approximation of viscous incompressible two-phase Stokes flow, Comput. Methods Appl. Mech. Engrg. 267 (2013), 511–530.10.1016/j.cma.2013.09.023Search in Google Scholar

[9] J. W. Barrett, H. Garcke, and R. Nürnberg, On the stable discretization of strongly anisotropic phase field models with applications to crystal growth, ZAMM Z. Angew. Math. Mech. 93 (2013), 719–732.10.1002/zamm.201200147Search in Google Scholar

[10] J. W. Barrett, H. Garcke, and R. Nürnberg, Stable phase field approximations of anisotropic solidification, IMA J. Numer. Anal. 34 (2014), 1289–1327.10.1093/imanum/drt044Search in Google Scholar

[11] J. W. Barrett, H. Garcke, and R. Nürnberg, A stable parametric finite element discretization of two-phase Navier–Stokes flow, J. Sci. Comp. 63 (2015), 78–117.10.1007/s10915-014-9885-2Search in Google Scholar

[12] J. W. Barrett, H. Garcke, and R. Nürnberg, On the stable numerical approximation of two-phase flow with insoluble surfactant, ESAIM Math. Model. Numer. Anal. 49 (2015), 421–458.10.1051/m2an/2014039Search in Google Scholar

[13] J. W. Barrett, H. Garcke, and R. Nürnberg, Parametric finite element approximations of curvature driven interface evolutions. In: Handb. Numer. Anal., Vol. 21 (Eds. A. Bonito and R. H. Nochetto), Elsevier, Amsterdam, 2020, pp. 275–423.10.1016/bs.hna.2019.05.002Search in Google Scholar

[14] P. W. Bates, X. Chen, and X. Deng, A numerical scheme for the two phase Mullins–Sekerka problem, Electron. J. Diff. Equ. 1995 (1995), 1–28.Search in Google Scholar

[15] G. Bellettini, M. Novaga, and M. Paolini, Facet-breaking for three-dimensional crystals evolving by mean curvature, Interfaces Free Bound. 1 (1999), 39–55.10.4171/IFB/3Search in Google Scholar

[16] J. W. Cahn and J. E. Hilliard, Free energy of a non-uniform system, I. Interfacial free energy, J. Chem. Phys. 28 (1958), 258–267.10.1063/1.1744102Search in Google Scholar

[17] J. W. Cahn and J. E. Taylor, Surface motion by surface diffusion, Acta Metall. Mater. 42 (1994), 1045–1063.10.1016/0956-7151(94)90123-6Search in Google Scholar

[18] S. Chen, B. Merriman, S. Osher, and P. Smereka, A simple level set method for solving Stefan problems, J. Comput. Phys. 135 (1997), 8–29.10.1006/jcph.1997.5721Search in Google Scholar

[19] X. Chen, The Hele-Shaw problem and area-preserving curve-shortening motions, Arch. Rational Mech. Anal. 123 (1993), 117–151.10.1007/BF00695274Search in Google Scholar

[20] X. Chen, J. Hong, and F. Yi, Existence, uniqueness, and regularity of classical solutions of the Mullins–Sekerka problem, Comm. Partial Diff. Equ. 21 (1996), 1705–1727.10.1080/03605309608821243Search in Google Scholar

[21] P. G. Ciarlet, The Finite Element Method for Elliptic Problems, Studies in Mathematics and its Applications, Vol. 4, North-Holland Publishing Co., Amsterdam, 1978.Search in Google Scholar

[22] T. A. Davis, Algorithm 832: UMFPACK V4.3 – an unsymmetric-pattern multifrontal method, ACM Trans. Math. Software 30 (2004), 196–199.10.1145/992200.992206Search in Google Scholar

[23] K. Deckelnick, G. Dziuk, and C. M. Elliott, Computation of geometric partial differential equations and mean curvature flow, Acta Numer. 14 (2005), 139–232.10.1017/S0962492904000224Search in Google Scholar

[24] K. Deckelnick and R. Nürnberg, A novel finite element approximation of anisotropic curve shortening flow, Interfaces Free Bound. (2023), to appear. See also arXiv:2110.04605, 2021.Search in Google Scholar

[25] G. Dziuk, An algorithm for evolutionary surfaces, Numer. Math. 58 (1991), 603–611.10.1007/BF01385643Search in Google Scholar

[26] J. Escher and G. Simonett, Classical solutions for Hele-Shaw models with surface tension, Adv. Diff. Equ. 2 (1997), 619–642.10.57262/ade/1366741151Search in Google Scholar

[27] X. Feng and A. Prohl, Analysis of a fully discrete finite element method for the phase field model and approximation of its sharp interface limits, Math. Comp. 73 (2004), 541–567.10.1090/S0025-5718-03-01588-6Search in Google Scholar

[28] X. Feng and A. Prohl, Error analysis of a mixed finite element method for the Cahn–Hilliard equation, Numer. Math. 99 (2004), 47–84.10.1007/s00211-004-0546-5Search in Google Scholar

[29] Y. Giga, Surface evolution equations, Monographs in Mathematics, Vol. 99, Birkhäuser, Basel, 2006.Search in Google Scholar

[30] W. Jiang and B. Li, A perimeter-decreasing and area-conserving algorithm for surface diffusion flow of curves, J. Comput. Phys. 443 (2021), 110531.10.1016/j.jcp.2021.110531Search in Google Scholar

[31] U. F. Mayer, Two-sided Mullins–Sekerka flow does not preserve convexity. Electron. J. Differ. Equ. Conf. 1 (1998), 171–179.Search in Google Scholar

[32] U. F. Mayer, A numerical scheme for moving boundary problems that are gradient flows for the area functional, European J. Appl. Math. 11 (2000), 61–80.10.1017/S0956792599003812Search in Google Scholar

[33] W. W. Mullins and R. F. Sekerka, Morphological stability of a particle growing by diffusion or heat flow, J. Appl. Phys. 34 (1963), 323–329.10.1063/1.1702607Search in Google Scholar

[34] A. Schmidt and K. G. Siebert, Design of Adaptive Finite Element Software: The Finite Element Toolbox ALBERTA, Lecture Notes in Computational Science and Engineering, Vol. 42, Springer-Verlag, Berlin, 2005.Search in Google Scholar

[35] J. E. Taylor and J. W. Cahn, Linking anisotropic sharp and diffuse surface motion laws via gradient flows, J. Statist. Phys. 77 (1994), 183–197.10.1007/BF02186838Search in Google Scholar

[36] J. E. Taylor, J. W. Cahn, and C. A. Handwerker, Geometric models of crystal growth, Acta Metall. Mater. 40 (1992), 1443–1474.10.1016/0956-7151(92)90090-2Search in Google Scholar

[37] J. Zhu, X. Chen, and T. Y. Hou, An efficient boundary integral method for the Mullins–Sekerka problem, J. Comput. Phys. 127 (1996), 246–267.10.1006/jcph.1996.0173Search in Google Scholar

Received: 2021-11-30
Revised: 2022-03-09
Accepted: 2022-03-09
Published Online: 2022-03-23
Published in Print: 2023-06-27

© 2023 Walter de Gruyter GmbH, Berlin/Boston

Downloaded on 26.11.2025 from https://www.degruyterbrill.com/document/doi/10.1515/jnma-2021-0131/html
Scroll to top button