3,422 research outputs found
Analysis of a high order Trace Finite Element Method for PDEs on level set surfaces
We present a new high order finite element method for the discretization of
partial differential equations on stationary smooth surfaces which are
implicitly described as the zero level of a level set function. The
discretization is based on a trace finite element technique. The higher
discretization accuracy is obtained by using an isoparametric mapping of the
volume mesh, based on the level set function, as introduced in [C. Lehrenfeld,
\emph{High order unfitted finite element methods on level set domains using
isoparametric mappings}, Comp. Meth. Appl. Mech. Engrg. 2016]. The resulting
trace finite element method is easy to implement. We present an error analysis
of this method and derive optimal order -norm error bounds. A
second topic of this paper is a unified analysis of several stabilization
methods for trace finite element methods. Only a stabilization method which is
based on adding an anisotropic diffusion in the volume mesh is able to control
the condition number of the stiffness matrix also for the case of higher order
discretizations. Results of numerical experiments are included which confirm
the theoretical findings on optimal order discretization errors and uniformly
bounded condition numbers.Comment: 28 pages, 5 figures, 1 tabl
Trace Finite Element Methods for PDEs on Surfaces
In this paper we consider a class of unfitted finite element methods for
discretization of partial differential equations on surfaces. In this class of
methods known as the Trace Finite Element Method (TraceFEM), restrictions or
traces of background surface-independent finite element functions are used to
approximate the solution of a PDE on a surface. We treat equations on steady
and time-dependent (evolving) surfaces. Higher order TraceFEM is explained in
detail. We review the error analysis and algebraic properties of the method.
The paper navigates through the known variants of the TraceFEM and the
literature on the subject
An adaptive octree finite element method for PDEs posed on surfaces
The paper develops a finite element method for partial differential equations
posed on hypersurfaces in , . The method uses traces of
bulk finite element functions on a surface embedded in a volumetric domain. The
bulk finite element space is defined on an octree grid which is locally refined
or coarsened depending on error indicators and estimated values of the surface
curvatures. The cartesian structure of the bulk mesh leads to easy and
efficient adaptation process, while the trace finite element method makes
fitting the mesh to the surface unnecessary. The number of degrees of freedom
involved in computations is consistent with the two-dimension nature of surface
PDEs. No parametrization of the surface is required; it can be given implicitly
by a level set function. In practice, a variant of the marching cubes method is
used to recover the surface with the second order accuracy. We prove the
optimal order of accuracy for the trace finite element method in and
surface norms for a problem with smooth solution and quasi-uniform mesh
refinement. Experiments with less regular problems demonstrate optimal
convergence with respect to the number of degrees of freedom, if grid
adaptation is based on an appropriate error indicator. The paper shows results
of numerical experiments for a variety of geometries and problems, including
advection-diffusion equations on surfaces. Analysis and numerical results of
the paper suggest that combination of cartesian adaptive meshes and the
unfitted (trace) finite elements provide simple, efficient, and reliable tool
for numerical treatment of PDEs posed on surfaces
A Trace Finite Element Method for Vector-Laplacians on Surfaces
We consider a vector-Laplace problem posed on a 2D surface embedded in a 3D
domain, which results from the modeling of surface fluids based on exterior
Cartesian differential operators. The main topic of this paper is the
development and analysis of a finite element method for the discretization of
this surface partial differential equation. We apply the trace finite element
technique, in which finite element spaces on a background shape-regular
tetrahedral mesh that is surface-independent are used for discretization. In
order to satisfy the constraint that the solution vector field is tangential to
the surface we introduce a Lagrange multiplier. We show well-posedness of the
resulting saddle point formulation. A discrete variant of this formulation is
introduced which contains suitable stabilization terms and is based on trace
finite element spaces. For this method we derive optimal discretization error
bounds. Furthermore algebraic properties of the resulting discrete saddle point
problem are studied. In particular an optimal Schur complement preconditioner
is proposed. Results of a numerical experiment are included
- …