2,726 research outputs found
A fast and well-conditioned spectral method for singular integral equations
We develop a spectral method for solving univariate singular integral
equations over unions of intervals by utilizing Chebyshev and ultraspherical
polynomials to reformulate the equations as almost-banded infinite-dimensional
systems. This is accomplished by utilizing low rank approximations for sparse
representations of the bivariate kernels. The resulting system can be solved in
operations using an adaptive QR factorization, where is
the bandwidth and is the optimal number of unknowns needed to resolve the
true solution. The complexity is reduced to operations by
pre-caching the QR factorization when the same operator is used for multiple
right-hand sides. Stability is proved by showing that the resulting linear
operator can be diagonally preconditioned to be a compact perturbation of the
identity. Applications considered include the Faraday cage, and acoustic
scattering for the Helmholtz and gravity Helmholtz equations, including
spectrally accurate numerical evaluation of the far- and near-field solution.
The Julia software package SingularIntegralEquations.jl implements our method
with a convenient, user-friendly interface
Structure-preserving mesh coupling based on the Buffa-Christiansen complex
The state of the art for mesh coupling at nonconforming interfaces is
presented and reviewed. Mesh coupling is frequently applied to the modeling and
simulation of motion in electromagnetic actuators and machines. The paper
exploits Whitney elements to present the main ideas. Both interpolation- and
projection-based methods are considered. In addition to accuracy and
efficiency, we emphasize the question whether the schemes preserve the
structure of the de Rham complex, which underlies Maxwell's equations. As a new
contribution, a structure-preserving projection method is presented, in which
Lagrange multiplier spaces are chosen from the Buffa-Christiansen complex. Its
performance is compared with a straightforward interpolation based on Whitney
and de Rham maps, and with Galerkin projection.Comment: 17 pages, 7 figures. Some figures are omitted due to a restricted
copyright. Full paper to appear in Mathematics of Computatio
Certified and fast computations with shallow covariance kernels
Many techniques for data science and uncertainty quantification demand
efficient tools to handle Gaussian random fields, which are defined in terms of
their mean functions and covariance operators. Recently, parameterized Gaussian
random fields have gained increased attention, due to their higher degree of
flexibility. However, especially if the random field is parameterized through
its covariance operator, classical random field discretization techniques fail
or become inefficient. In this work we introduce and analyze a new and
certified algorithm for the low-rank approximation of a parameterized family of
covariance operators which represents an extension of the adaptive cross
approximation method for symmetric positive definite matrices. The algorithm
relies on an affine linear expansion of the covariance operator with respect to
the parameters, which needs to be computed in a preprocessing step using, e.g.,
the empirical interpolation method. We discuss and test our new approach for
isotropic covariance kernels, such as Mat\'ern kernels. The numerical results
demonstrate the advantages of our approach in terms of computational time and
confirm that the proposed algorithm provides the basis of a fast sampling
procedure for parameter dependent Gaussian random fields
An overlapped grid method for multigrid, finite volume/difference flow solvers: MaGGiE
The objective is to develop a domain decomposition method via overlapping/embedding the component grids, which is to be used by upwind, multi-grid, finite volume solution algorithms. A computer code, given the name MaGGiE (Multi-Geometry Grid Embedder) is developed to meet this objective. MaGGiE takes independently generated component grids as input, and automatically constructs the composite mesh and interpolation data, which can be used by the finite volume solution methods with or without multigrid convergence acceleration. Six demonstrative examples showing various aspects of the overlap technique are presented and discussed. These cases are used for developing the procedure for overlapping grids of different topologies, and to evaluate the grid connection and interpolation data for finite volume calculations on a composite mesh. Time fluxes are transferred between mesh interfaces using a trilinear interpolation procedure. Conservation losses are minimal at the interfaces using this method. The multi-grid solution algorithm, using the coaser grid connections, improves the convergence time history as compared to the solution on composite mesh without multi-gridding
- …