1,841 research outputs found
Near-optimal perfectly matched layers for indefinite Helmholtz problems
A new construction of an absorbing boundary condition for indefinite
Helmholtz problems on unbounded domains is presented. This construction is
based on a near-best uniform rational interpolant of the inverse square root
function on the union of a negative and positive real interval, designed with
the help of a classical result by Zolotarev. Using Krein's interpretation of a
Stieltjes continued fraction, this interpolant can be converted into a
three-term finite difference discretization of a perfectly matched layer (PML)
which converges exponentially fast in the number of grid points. The
convergence rate is asymptotically optimal for both propagative and evanescent
wave modes. Several numerical experiments and illustrations are included.Comment: Accepted for publication in SIAM Review. To appear 201
Refined a posteriori error estimation for classical and pressure-robust Stokes finite element methods
Recent works showed that pressure-robust modifications of mixed finite
element methods for the Stokes equations outperform their standard versions in
many cases. This is achieved by divergence-free reconstruction operators and
results in pressure independent velocity error estimates which are robust with
respect to small viscosities. In this paper we develop a posteriori error
control which reflects this robustness.
The main difficulty lies in the volume contribution of the standard
residual-based approach that includes the -norm of the right-hand side.
However, the velocity is only steered by the divergence-free part of this
source term. An efficient error estimator must approximate this divergence-free
part in a proper manner, otherwise it can be dominated by the pressure error.
To overcome this difficulty a novel approach is suggested that uses arguments
from the stream function and vorticity formulation of the Navier--Stokes
equations. The novel error estimators only take the of the
right-hand side into account and so lead to provably reliable, efficient and
pressure-independent upper bounds in case of a pressure-robust method in
particular in pressure-dominant situations. This is also confirmed by some
numerical examples with the novel pressure-robust modifications of the
Taylor--Hood and mini finite element methods
Div First-Order System LL* (FOSLL*) for Second-Order Elliptic Partial Differential Equations
The first-order system LL* (FOSLL*) approach for general second-order
elliptic partial differential equations was proposed and analyzed in [10], in
order to retain the full efficiency of the L2 norm first-order system
least-squares (FOSLS) ap- proach while exhibiting the generality of the
inverse-norm FOSLS approach. The FOSLL* approach in [10] was applied to the
div-curl system with added slack vari- ables, and hence it is quite
complicated. In this paper, we apply the FOSLL* approach to the div system and
establish its well-posedness. For the corresponding finite ele- ment
approximation, we obtain a quasi-optimal a priori error bound under the same
regularity assumption as the standard Galerkin method, but without the
restriction to sufficiently small mesh size. Unlike the FOSLS approach, the
FOSLL* approach does not have a free a posteriori error estimator, we then
propose an explicit residual error estimator and establish its reliability and
efficiency bound
- …