- Journal Home
- Volume 42 - 2024
- Volume 41 - 2023
- Volume 40 - 2022
- Volume 39 - 2021
- Volume 38 - 2020
- Volume 37 - 2019
- Volume 36 - 2018
- Volume 35 - 2017
- Volume 34 - 2016
- Volume 33 - 2015
- Volume 32 - 2014
- Volume 31 - 2013
- Volume 30 - 2012
- Volume 29 - 2011
- Volume 28 - 2010
- Volume 27 - 2009
- Volume 26 - 2008
- Volume 25 - 2007
- Volume 24 - 2006
- Volume 23 - 2005
- Volume 22 - 2004
- Volume 21 - 2003
- Volume 20 - 2002
- Volume 19 - 2001
- Volume 18 - 2000
- Volume 17 - 1999
- Volume 16 - 1998
- Volume 15 - 1997
- Volume 14 - 1996
- Volume 13 - 1995
- Volume 12 - 1994
- Volume 11 - 1993
- Volume 10 - 1992
- Volume 9 - 1991
- Volume 8 - 1990
- Volume 7 - 1989
- Volume 6 - 1988
- Volume 5 - 1987
- Volume 4 - 1986
- Volume 3 - 1985
- Volume 2 - 1984
- Volume 1 - 1983
Cited by
- BibTex
- RIS
- TXT
In this article we define a surface finite element method (SFEM) for the numerical solution of parabolic partial differential equations on hypersurfaces $\Gamma$ in $\mathbb R^{n+1}$. The key idea is based on the approximation of $\Gamma$ by a polyhedral surface $\Gamma_h$ consisting of a union of simplices (triangles for $n=2$, intervals for $n=1$) with vertices on $\Gamma$. A finite element space of functions is then defined by taking the continuous functions on $\Gamma_h$ which are linear affine on each simplex of the polygonal surface. We use surface gradients to define weak forms of elliptic operators and naturally generate weak formulations of elliptic and parabolic equations on $\Gamma$. Our finite element method is applied to weak forms of the equations. The computation of the mass and element stiffness matrices is simple and straightforward. We give an example of error bounds in the case of semi-discretization in space for a fourth order linear problem. Numerical experiments are described for several linear and nonlinear partial differential equations. In particular, the power of the method is demonstrated by employing it to solve highly nonlinear second and fourth order problems such as surface Allen-Cahn and Cahn-Hilliard equations and surface level set equations for geodesic mean curvature flow.
}, issn = {1991-7139}, doi = {https://doi.org/}, url = {http://global-sci.org/intro/article_detail/jcm/8700.html} }In this article we define a surface finite element method (SFEM) for the numerical solution of parabolic partial differential equations on hypersurfaces $\Gamma$ in $\mathbb R^{n+1}$. The key idea is based on the approximation of $\Gamma$ by a polyhedral surface $\Gamma_h$ consisting of a union of simplices (triangles for $n=2$, intervals for $n=1$) with vertices on $\Gamma$. A finite element space of functions is then defined by taking the continuous functions on $\Gamma_h$ which are linear affine on each simplex of the polygonal surface. We use surface gradients to define weak forms of elliptic operators and naturally generate weak formulations of elliptic and parabolic equations on $\Gamma$. Our finite element method is applied to weak forms of the equations. The computation of the mass and element stiffness matrices is simple and straightforward. We give an example of error bounds in the case of semi-discretization in space for a fourth order linear problem. Numerical experiments are described for several linear and nonlinear partial differential equations. In particular, the power of the method is demonstrated by employing it to solve highly nonlinear second and fourth order problems such as surface Allen-Cahn and Cahn-Hilliard equations and surface level set equations for geodesic mean curvature flow.