A finite element method for anisotropic crystal growth on surfaces

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Garcke, Harald, Nürnberg, Robert
Format: Preprint
Published: 2024
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866913879501570048
author Garcke, Harald
Nürnberg, Robert
author_facet Garcke, Harald
Nürnberg, Robert
contents Phase transition problems on curved surfaces can lead to a panopticon of fascinating patterns. In this paper we consider finite element approximations of phase field models with a spatially inhomogeneous and anisotropic surface energy density. The problems are either posed in $\mathbb R^3$ or on a two-dimensional hypersurface in $\mathbb R^3$. In the latter case, a fundamental choice regarding the anisotropic energy density has to be made. One possibility is to use a density defined in the ambient space $\mathbb R^3$. However, we propose and advocate for an alternative, where a density is defined on a fixed chosen tangent space, and is then moved along geodesics to the other tangent spaces. Our numerical method can be employed in all of the above situations, where for the problems on hypersurfaces the algorithm uses parametric finite elements. We prove an unconditional stability result for our schemes and present several numerical experiments, including for the modelling of ice crystal growth on a sphere.
format Preprint
id arxiv_https___arxiv_org_abs_2403_14206
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle A finite element method for anisotropic crystal growth on surfaces
Garcke, Harald
Nürnberg, Robert
Numerical Analysis
Materials Science
Computational Physics
Phase transition problems on curved surfaces can lead to a panopticon of fascinating patterns. In this paper we consider finite element approximations of phase field models with a spatially inhomogeneous and anisotropic surface energy density. The problems are either posed in $\mathbb R^3$ or on a two-dimensional hypersurface in $\mathbb R^3$. In the latter case, a fundamental choice regarding the anisotropic energy density has to be made. One possibility is to use a density defined in the ambient space $\mathbb R^3$. However, we propose and advocate for an alternative, where a density is defined on a fixed chosen tangent space, and is then moved along geodesics to the other tangent spaces. Our numerical method can be employed in all of the above situations, where for the problems on hypersurfaces the algorithm uses parametric finite elements. We prove an unconditional stability result for our schemes and present several numerical experiments, including for the modelling of ice crystal growth on a sphere.
title A finite element method for anisotropic crystal growth on surfaces
topic Numerical Analysis
Materials Science
Computational Physics
url https://arxiv.org/abs/2403.14206