12,549 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
A posteriori error control for discontinuous Galerkin methods for parabolic problems
We derive energy-norm a posteriori error bounds for an Euler time-stepping
method combined with various spatial discontinuous Galerkin schemes for linear
parabolic problems. For accessibility, we address first the spatially
semidiscrete case, and then move to the fully discrete scheme by introducing
the implicit Euler time-stepping. All results are presented in an abstract
setting and then illustrated with particular applications. This enables the
error bounds to hold for a variety of discontinuous Galerkin methods, provided
that energy-norm a posteriori error bounds for the corresponding elliptic
problem are available. To illustrate the method, we apply it to the interior
penalty discontinuous Galerkin method, which requires the derivation of novel a
posteriori error bounds. For the analysis of the time-dependent problems we use
the elliptic reconstruction technique and we deal with the nonconforming part
of the error by deriving appropriate computable a posteriori bounds for it.Comment: 6 figure
- …