828 research outputs found
Partially Penalized Immersed Finite Element Methods for Elliptic Interface Problems
This article presents new immersed finite element (IFE) methods for solving
the popular second order elliptic interface problems on structured Cartesian
meshes even if the involved interfaces have nontrivial geometries. These IFE
methods contain extra stabilization terms introduced only at interface edges
for penalizing the discontinuity in IFE functions. With the enhanced stability
due to the added penalty, not only these IFE methods can be proven to have the
optimal convergence rate in the H1-norm provided that the exact solution has
sufficient regularity, but also numerical results indicate that their
convergence rates in both the H1-norm and the L2-norm do not deteriorate when
the mesh becomes finer which is a shortcoming of the classic IFE methods in
some situations. Trace inequalities are established for both linear and
bilinear IFE functions that are not only critical for the error analysis of
these new IFE methods, but also are of a great potential to be useful in error
analysis for other IFE methods
Accurate gradient computations at interfaces using finite element methods
New finite element methods are proposed for elliptic interface problems in
one and two dimensions. The main motivation is not only to get an accurate
solution but also an accurate first order derivative at the interface (from
each side). The key in 1D is to use the idea from \cite{wheeler1974galerkin}.
For 2D interface problems, the idea is to introduce a small tube near the
interface and introduce the gradient as part of unknowns, which is similar to a
mixed finite element method, except only at the interface. Thus the
computational cost is just slightly higher than the standard finite element
method. We present rigorous one dimensional analysis, which show second order
convergence order for both of the solution and the gradient in 1D. For two
dimensional problems, we present numerical results and observe second order
convergence for the solution, and super-convergence for the gradient at the
interface
- …