70,726 research outputs found
Adaptive BEM with optimal convergence rates for the Helmholtz equation
We analyze an adaptive boundary element method for the weakly-singular and
hypersingular integral equations for the 2D and 3D Helmholtz problem. The
proposed adaptive algorithm is steered by a residual error estimator and does
not rely on any a priori information that the underlying meshes are
sufficiently fine. We prove convergence of the error estimator with optimal
algebraic rates, independently of the (coarse) initial mesh. As a technical
contribution, we prove certain local inverse-type estimates for the boundary
integral operators associated with the Helmholtz equation
On the Convergence of the Born Series in Optical Tomography with Diffuse Light
We provide a simple sufficient condition for convergence of Born series in
the forward problem of optical diffusion tomography. The condition does not
depend on the shape or spatial extent of the inhomogeneity but only on its
amplitude.Comment: 23 pages, 7 figures, submitted to Inverse Problem
Solving Dirac equations on a 3D lattice with inverse Hamiltonian and spectral methods
A new method to solve the Dirac equation on a 3D lattice is proposed, in
which the variational collapse problem is avoided by the inverse Hamiltonian
method and the fermion doubling problem is avoided by performing spatial
derivatives in momentum space with the help of the discrete Fourier transform,
i.e., the spectral method. This method is demonstrated in solving the Dirac
equation for a given spherical potential in 3D lattice space. In comparison
with the results obtained by the shooting method, the differences in single
particle energy are smaller than ~MeV, and the densities are almost
identical, which demonstrates the high accuracy of the present method. The
results obtained by applying this method without any modification to solve the
Dirac equations for an axial deformed, non-axial deformed, and octupole
deformed potential are provided and discussed.Comment: 18 pages, 6 figure
Fully discrete finite element data assimilation method for the heat equation
We consider a finite element discretization for the reconstruction of the
final state of the heat equation, when the initial data is unknown, but
additional data is given in a sub domain in the space time. For the
discretization in space we consider standard continuous affine finite element
approximation, and the time derivative is discretized using a backward
differentiation. We regularize the discrete system by adding a penalty of the
-semi-norm of the initial data, scaled with the mesh-parameter. The
analysis of the method uses techniques developed in E. Burman and L. Oksanen,
Data assimilation for the heat equation using stabilized finite element
methods, arXiv, 2016, combining discrete stability of the numerical method with
sharp Carleman estimates for the physical problem, to derive optimal error
estimates for the approximate solution. For the natural space time energy norm,
away from , the convergence is the same as for the classical problem with
known initial data, but contrary to the classical case, we do not obtain faster
convergence for the -norm at the final time
Discretely exact derivatives for hyperbolic PDE-constrained optimization problems discretized by the discontinuous Galerkin method
This paper discusses the computation of derivatives for optimization problems
governed by linear hyperbolic systems of partial differential equations (PDEs)
that are discretized by the discontinuous Galerkin (dG) method. An efficient
and accurate computation of these derivatives is important, for instance, in
inverse problems and optimal control problems. This computation is usually
based on an adjoint PDE system, and the question addressed in this paper is how
the discretization of this adjoint system should relate to the dG
discretization of the hyperbolic state equation. Adjoint-based derivatives can
either be computed before or after discretization; these two options are often
referred to as the optimize-then-discretize and discretize-then-optimize
approaches. We discuss the relation between these two options for dG
discretizations in space and Runge-Kutta time integration. Discretely exact
discretizations for several hyperbolic optimization problems are derived,
including the advection equation, Maxwell's equations and the coupled
elastic-acoustic wave equation. We find that the discrete adjoint equation
inherits a natural dG discretization from the discretization of the state
equation and that the expressions for the discretely exact gradient often have
to take into account contributions from element faces. For the coupled
elastic-acoustic wave equation, the correctness and accuracy of our derivative
expressions are illustrated by comparisons with finite difference gradients.
The results show that a straightforward discretization of the continuous
gradient differs from the discretely exact gradient, and thus is not consistent
with the discretized objective. This inconsistency may cause difficulties in
the convergence of gradient based algorithms for solving optimization problems
A finite element data assimilation method for the wave equation
We design a primal-dual stabilized finite element method for the numerical
approximation of a data assimilation problem subject to the acoustic wave
equation. For the forward problem, piecewise affine, continuous, finite element
functions are used for the approximation in space and backward differentiation
is used in time. Stabilizing terms are added on the discrete level. The design
of these terms is driven by numerical stability and the stability of the
continuous problem, with the objective of minimizing the computational error.
Error estimates are then derived that are optimal with respect to the
approximation properties of the numerical scheme and the stability properties
of the continuous problem. The effects of discretizing the (smooth) domain
boundary and other perturbations in data are included in the analysis.Comment: 23 page
- …