7,745 research outputs found
Convergence of HX Preconditioner for Maxwell's Equations with Jump Coefficients (ii): The Main Results
This paper is the second one of two serial articles, whose goal is to prove
convergence of HX Preconditioner (proposed by Hiptmair and Xu, 2007) for
Maxwell's equations with jump coefficients. In this paper, based on the
auxiliary results developed in the first paper (Hu, 2017), we establish a new
regular Helmholtz decomposition for edge finite element functions in three
dimensions, which is nearly stable with respect to a weight function. By using
this Helmholtz decomposition, we give an analysis of the convergence of the HX
preconditioner for the case with strongly discontinuous coefficients. We show
that the HX preconditioner possesses fast convergence, which not only is nearly
optimal with respect to the finite element mesh size but also is independent of
the jumps in the coefficients across the interface between two neighboring
subdomains.Comment: with 25 pages, 2 figure
Shifted Laplacian multigrid for the elastic Helmholtz equation
The shifted Laplacian multigrid method is a well known approach for
preconditioning the indefinite linear system arising from the discretization of
the acoustic Helmholtz equation. This equation is used to model wave
propagation in the frequency domain. However, in some cases the acoustic
equation is not sufficient for modeling the physics of the wave propagation,
and one has to consider the elastic Helmholtz equation. Such a case arises in
geophysical seismic imaging applications, where the earth's subsurface is the
elastic medium. The elastic Helmholtz equation is much harder to solve than its
acoustic counterpart, partially because it is three times larger, and partially
because it models more complicated physics. Despite this, there are very few
solvers available for the elastic equation compared to the array of solvers
that are available for the acoustic one. In this work we extend the shifted
Laplacian approach to the elastic Helmholtz equation, by combining the complex
shift idea with approaches for linear elasticity. We demonstrate the efficiency
and properties of our solver using numerical experiments for problems with
heterogeneous media in two and three dimensions
Domain Decomposition preconditioning for high-frequency Helmholtz problems with absorption
In this paper we give new results on domain decomposition preconditioners for
GMRES when computing piecewise-linear finite-element approximations of the
Helmholtz equation , with
absorption parameter . Multigrid approximations of
this equation with are commonly used as preconditioners
for the pure Helmholtz case (). However a rigorous theory for
such (so-called "shifted Laplace") preconditioners, either for the pure
Helmholtz equation, or even the absorptive equation (), is
still missing. We present a new theory for the absorptive equation that
provides rates of convergence for (left- or right-) preconditioned GMRES, via
estimates of the norm and field of values of the preconditioned matrix. This
theory uses a - and -explicit coercivity result for the
underlying sesquilinear form and shows, for example, that if , then classical overlapping additive Schwarz will perform optimally for
the absorptive problem, provided the subdomain and coarse mesh diameters are
carefully chosen. Extensive numerical experiments are given that support the
theoretical results. The theory for the absorptive case gives insight into how
its domain decomposition approximations perform as preconditioners for the pure
Helmholtz case . At the end of the paper we propose a
(scalable) multilevel preconditioner for the pure Helmholtz problem that has an
empirical computation time complexity of about for
solving finite element systems of size , where we have
chosen the mesh diameter to avoid the pollution effect.
Experiments on problems with , i.e. a fixed number of grid points
per wavelength, are also given
Compressive Wave Computation
This paper considers large-scale simulations of wave propagation phenomena.
We argue that it is possible to accurately compute a wavefield by decomposing
it onto a largely incomplete set of eigenfunctions of the Helmholtz operator,
chosen at random, and that this provides a natural way of parallelizing wave
simulations for memory-intensive applications.
This paper shows that L1-Helmholtz recovery makes sense for wave computation,
and identifies a regime in which it is provably effective: the one-dimensional
wave equation with coefficients of small bounded variation. Under suitable
assumptions we show that the number of eigenfunctions needed to evolve a sparse
wavefield defined on N points, accurately with very high probability, is
bounded by C log(N) log(log(N)), where C is related to the desired accuracy and
can be made to grow at a much slower rate than N when the solution is sparse.
The PDE estimates that underlie this result are new to the authors' knowledge
and may be of independent mathematical interest; they include an L1 estimate
for the wave equation, an estimate of extension of eigenfunctions, and a bound
for eigenvalue gaps in Sturm-Liouville problems.
Numerical examples are presented in one spatial dimension and show that as
few as 10 percents of all eigenfunctions can suffice for accurate results.
Finally, we argue that the compressive viewpoint suggests a competitive
parallel algorithm for an adjoint-state inversion method in reflection
seismology.Comment: 45 pages, 4 figure
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
- …