49 research outputs found
Harmonic density interpolation methods for high-order evaluation of Laplace layer potentials in 2D and 3D
We present an effective harmonic density interpolation method for the
numerical evaluation of singular and nearly singular Laplace boundary integral
operators and layer potentials in two and three spatial dimensions. The method
relies on the use of Green's third identity and local Taylor-like
interpolations of density functions in terms of harmonic polynomials. The
proposed technique effectively regularizes the singularities present in
boundary integral operators and layer potentials, and recasts the latter in
terms of integrands that are bounded or even more regular, depending on the
order of the density interpolation. The resulting boundary integrals can then
be easily, accurately, and inexpensively evaluated by means of standard
quadrature rules. A variety of numerical examples demonstrate the effectiveness
of the technique when used in conjunction with the classical trapezoidal rule
(to integrate over smooth curves) in two-dimensions, and with a Chebyshev-type
quadrature rule (to integrate over surfaces given as unions of non-overlapping
quadrilateral patches) in three-dimensions
Second-Kind integral solvers for TE and TM problems of diffraction by open-arcs
We present a novel approach for the numerical solution of problems of
diffraction by open arcs in two dimensional space. Our methodology relies on
composition of {\em weighted versions} of the classical integral operators
associated with the Dirichlet and Neumann problems (TE and TM polarizations,
respectively) together with a generalization to the open-arc case of the well
known closed-surface Calder\'on formulae. When used in conjunction with
spectrally accurate discretization rules and Krylov-subspace linear algebra
solvers such as GMRES, the new second-kind TE and TM formulations for open arcs
produce results of high accuracy in small numbers of iterations and short
computing times---for low and high frequencies alike.Comment: 20 page
Weighted integral solvers for elastic scattering by open arcs in two dimensions
We present new methodologies for the numerical solution of problems of elastic scattering by open arcs in two dimensions. The algorithms utilize weighted versions of the classical elastic integral operators associated with Dirichlet and Neumann boundary conditions, where the integral weight accounts for (and regularizes) the singularity of the integralâequation solutions at the openâarc endpoints. Crucially, the method also incorporates a certain âopenâarc elastic CalderĂłn relationâ introduced in this paper, whose validity is demonstrated on the basis of numerical experiments, but whose rigorous mathematical proof is left for future work. (In fact, the aforementioned openâarc elastic CalderĂłn relation generalizes a corresponding elastic CalderĂłn relation for closed surfaces, which is also introduced in this paper, and for which a rigorous proof is included.) Using the openâsurface CalderĂłn relation in conjunction with spectrally accurate quadrature rules and the Krylovâsubspace linear algebra solver GMRES, the proposed overall openâarc elastic solver produces results of high accuracy in small number of iterations, for both low and high frequencies. A variety of numerical examples in this paper demonstrate the accuracy and efficiency of the proposed methodology
A high-order integral solver for scalar problems of diffraction by screens and apertures in three-dimensional space
We present a novel methodology for the numerical solution of problems of diffraction by infinitely thin screens in three-dimensional space. Our approach relies on new integral formulations as well as associated high-order quadrature rules. The new integral formulations involve weighted versions of the classical integral operators related to the thin-screen Dirichlet and Neumann problems as well as a generalization to the open-surface problem of the classical CalderĂłn formulae. The high-order quadrature rules we introduce for these operators, in turn, resolve the multiple Green function and edge singularities (which occur at arbitrarily close distances from each other, and which include weakly singular as well as hypersingular kernels) and thus give rise to super-algebraically fast convergence as the discretization sizes are increased. When used in conjunction with Krylov-subspace linear algebra solvers such as GMRES, the resulting solvers produce results of high accuracy in small numbers of iterations for low and high frequencies alike. We demonstrate our methodology with a variety of numerical results for screen and aperture problems at high frequenciesâincluding simulation of classical experiments such as the diffraction by a circular disc (featuring in particular the famous Poisson spot), evaluation of interference fringes resulting from diffraction across two nearby circular apertures, as well as solution of problems of scattering by more complex geometries consisting of multiple scatterers and cavities
Fast Numerical Methods for Non-local Operators
[no abstract available
Windowed Green function method for wave scattering by periodic arrays of 2D obstacles
This paper introduces a novel boundary integral equation (BIE) method for the
numerical solution of problems of planewave scattering by periodic line arrays
of two-dimensional penetrable obstacles. Our approach is built upon a direct
BIE formulation that leverages the simplicity of the free-space Green function
but in turn entails evaluation of integrals over the unit-cell boundaries. Such
integrals are here treated via the window Green function method. The windowing
approximation together with a finite-rank operator correction -- used to
properly impose the Rayleigh radiation condition -- yield a robust second-kind
BIE that produces super-algebraically convergent solutions throughout the
spectrum, including at the challenging Rayleigh-Wood anomalies. The corrected
windowed BIE can be discretized by means of off-the-shelf Nystr\"om and
boundary element methods, and it leads to linear systems suitable for iterative
linear-algebra solvers as well as standard fast matrix-vector product
algorithms. A variety of numerical examples demonstrate the accuracy and
robustness of the proposed methodolog
A new integral representation for quasiperiodic fields and its application to two-dimensional band structure calculations
In this paper, we consider band-structure calculations governed by the
Helmholtz or Maxwell equations in piecewise homogeneous periodic materials.
Methods based on boundary integral equations are natural in this context, since
they discretize the interface alone and can achieve high order accuracy in
complicated geometries. In order to handle the quasi-periodic conditions which
are imposed on the unit cell, the free-space Green's function is typically
replaced by its quasi-periodic cousin. Unfortunately, the quasi-periodic
Green's function diverges for families of parameter values that correspond to
resonances of the empty unit cell. Here, we bypass this problem by means of a
new integral representation that relies on the free-space Green's function
alone, adding auxiliary layer potentials on the boundary of the unit cell
itself. An important aspect of our method is that by carefully including a few
neighboring images, the densities may be kept smooth and convergence rapid.
This framework results in an integral equation of the second kind, avoids
spurious resonances, and achieves spectral accuracy. Because of our image
structure, inclusions which intersect the unit cell walls may be handled easily
and automatically. Our approach is compatible with fast-multipole acceleration,
generalizes easily to three dimensions, and avoids the complication of
divergent lattice sums.Comment: 25 pages, 6 figures, submitted to J. Comput. Phy