Application of Adaptive Finite Element Method for Elliptic Partial Differential Equations to the Laplace Beltrami Operator on Graphs

Khamron Mekchay

Bulletin of the Malaysian Mathematical Society, Tome 35 (2012) no. 2A

Résumé

The Laplace Beltrami operator, known as an elliptic operator for functions defined on surfaces, appears in some applications in sciences and engineerings. In this paper we consider the Laplace Beltrami operator $\Delta_\Gamma$ on surfaces $\Gamma$ defined as graphs of $C^2$ functions on a flat domain $\Omega \subset \mathbb{R}^{d-1}$ ($d\geq 2$), $ \Delta_\Gamma u = f \ \text{ on }\ \Gamma, \qquad u = 0\ \text{ on }\ \partial\Gamma. $ Based on some properties of differential geometry, we transformed the Laplace Beltrami operator on the surface $\Gamma$ to get an elliptic operator on the flat domain $\Omega$, $ -\text{div}(\mathbf{A}(\nabla u)^T) = F \ \text{ in }\ \Omega, \qquad u = 0 \ \text{ on }\ \partial\Omega. $ We applied an adaptive finite element method (AFEM) for a general second order linear elliptic partial differential equations developed by K. Mekchay and R. H. Nochetto to solve the transformed problem. The a posteriori error estimates in energy norm and the design or algorithm are derived accordingly for the transformed problem in the sense of the elliptic operator. The discretization and mesh generations of the AFEM algorithm rely on indicators and oscillations which now depend on the data $\mathbf{A}$ and $F$ of the elliptic operator on $\Omega$ and do not involve the geometric property of the surface $\Gamma$. A numerical experiment for the AFEM algorithm is provided to illustrate the theoretical results.

Détail

Classification : 65N12, 65N15, 65N30, 65N50, 65Y20.