• Recursively compressed inverse preconditioning (RCIP) is a kernel-independent and purely numerical method for solving Fredholm second kind boundary integral equations in situations where the boundary shape induces a non-smooth behavior in the solution. The method originated in 2008 within a scheme for Laplace's equation in two-dimensional domains with corners. In a series of subsequent papers the method was then refined and extended as to apply to integral equation formulations of a broad range of boundary value problems in physics and engineering. The purpose of the present tutorial is threefold: First, to review the RCIP method in a simple setting. Second, to show how easily the method can be implemented in Matlab. Third, to present new applications.
  • A high-order convergent and robust numerical solver is constructed and used to find complex eigenwavenumbers and electromagnetic eigenfields of dielectric objects with axial symmetry. The solver is based on Fourier--Nystr\"om discretization of combined integral equations for the transmission problem and can be applied to demanding resonance problems at microwave, terahertz, and optical wavelengths. High achievable accuracy, even at very high wavenumbers, makes the solver ideal for benchmarking and for assessing the performance of general purpose commercial software.
  • We consider the perturbation of an electric potential due to an insulating inclusion with corners. This perturbation is known to admit a multipole expansion whose coefficients are linear combinations of generalized polarization tensors. We define new geometric factors of a simple planar domain in terms of a conformal mapping associated with the domain. The geometric factors share properties of the generalized polarization tensors and are the Fourier series coefficients of a kind of generalized external angle of the inclusion boundary. Since the generalized external angle contains the Dirac delta singularity at corner points, we can determine the criterion for the existence of corner points on the inclusion boundary in terms of the geometric factors. We illustrate and validate our results with numerical examples computed to a high degree of precision using integral equation techniques, Nystr\"om discretization, and recursively compressed inverse preconditioning.
  • A particular mix of integral equations and discretization techniques is suggested for the solution of a planar Helmholtz transmission problem with relevance to the study of surface plasmon waves. The transmission problem describes the scattering of a time harmonic transverse magnetic wave from an infinite dielectric cylinder with complex permittivity and sharp edges. Numerical examples illustrate that the resulting scheme is capable of obtaining total magnetic and electric fields to very high accuracy in the entire computational domain.
  • We study the adjoint of the double layer potential associated with the Laplacian (the adjoint of the Neumann-Poincar\'e operator), as a map on the boundary surface $\Gamma$ of a domain in $\mathbb{R}^3$ with conical points. The spectrum of this operator directly reflects the well-posedness of related transmission problems across $\Gamma$. In particular, if the domain is understood as an inclusion with complex permittivity $\epsilon$, embedded in a background medium with unit permittivity, then the polarizability tensor of the domain is well-defined when $(\epsilon+1)/(\epsilon-1)$ belongs to the resolvent set in energy norm. We study surfaces $\Gamma$ that have a finite number of conical points featuring rotational symmetry. On the energy space, we show that the essential spectrum consists of an interval. On $L^2(\Gamma)$, i.e. for square-integrable boundary data, we show that the essential spectrum consists of a countable union of curves, outside of which the Fredholm index can be computed as a winding number with respect to the essential spectrum. We provide explicit formulas, depending on the opening angles of the conical points. We reinforce our study with very precise numerical experiments, computing the energy space spectrum and the spectral measures of the polarizability tensor in two different examples. Our results indicate that the densities of the spectral measures may approach zero extremely rapidly in the continuous part of the energy space spectrum.
  • A high-order accurate, explicit kernel-split, panel-based, Fourier-Nystr\"om discretization scheme is developed for integral equations associated with the Helmholtz equation in axially symmetric domains. Extensive incorporation of analytic information about singular integral kernels and on-the-fly computation of nearly singular quadrature rules allow for very high achievable accuracy, also in the evaluation of fields close to the boundary of the computational domain.
  • We study spectral properties of the Neumann-Poincar\'e operator on planar domains with corners with particular emphasis on existence of continuous spectrum and pure point spectrum. We show that the rate of resonance at continuous spectrum is different from that at eigenvalues, and then derive a method to distinguish continuous spectrum from eigenvalues. We perform computational experiments using the method to see whether continuous spectrum and pure point spectrum appear on domains with corners. For the computations we use a modification of the Nystr\"om method which makes it possible to construct high-order convergent discretizations of the Neumann-Poincar\'e operator on domains with corners. The results of experiments show that all three possible spectra, absolutely continuous spectrum, singularly continuous spectrum, and pure point spectrum, may appear depending on domains. We also prove rigorously two properties of spectrum which are suggested by numerical experiments: symmetry of spectrum (including continuous spectrum), and existence of eigenvalues on rectangles of high aspect ratio.
  • The magnetic field integral equation for axially symmetric cavities with perfectly conducting piecewise smooth surfaces is discretized according to a high-order convergent Fourier--Nystr\"om scheme. The resulting solver is used to accurately determine eigenwavenumbers and normalized electric eigenfields in the entire computational domain.
  • The magnetic field integral equation for axially symmetric cavities with perfectly conducting surfaces is discretized according to a high-order convergent Fourier--Nystr\"om scheme. The resulting solver is used to determine eigenwavenumbers and normalized magnetic eigenfields to very high accuracy in the entire computational domain.
  • The incorporation of analytical kernel information is exploited in the construction of Nystr\"om discretization schemes for integral equations modeling planar Helmholtz boundary value problems. Splittings of kernels and matrices, coarse and fine grids, high-order polynomial interpolation, product integration performed on the fly, and iterative solution are some of the numerical techniques used to seek rapid and stable convergence of computed fields in the entire computational domain.
  • This note is about promoting singularity subtraction as a helpful tool in the discretization of singular integral operators on curved surfaces. Singular and nearly singular kernels are expanded in series whose terms are integrated on parametrically rectangular regions using high-order product integration, thereby reducing the need for spatial adaptivity and precomputed weights. A simple scheme is presented and an application to the interior Dirichlet Laplace problem on some tori gives around ten digit accurate results using only two expansion terms and a modest programming- and computational effort.
  • In this paper we consider the classic problems of scattering of waves from perfectly conducting cylinders with piecewise smooth boundaries. The scattering problems are formulated as integral equations and solved using a Nystr\"om scheme where the corners of the cylinders are efficiently handled by a method referred to as Recursively Compressed Inverse Preconditioning (RCIP). This method has been very successful in treating static problems in non-smooth domains and the present paper shows that it works equally well for the Helmholtz equation. In the numerical examples we specialize to scattering of E- and H-waves from a cylinder with one corner. Even at a size kd=1000, where k is the wavenumber and d the diameter, the scheme produces at least 13 digits of accuracy in the electric and magnetic fields everywhere outside the cylinder.
  • An efficient integral equation based solver is constructed for the electrostatic problem on domains with cuboidal inclusions. It can be used to compute the polarizability of a dielectric cube in a dielectric background medium at virtually every permittivity ratio for which it exists. For example, polarizabilities accurate to between five and ten digits are obtained (as complex limits) for negative permittivity ratios in minutes on a standard workstation. In passing, the capacitance of the unit cube is determined with unprecedented accuracy. With full rigor, we develop a natural mathematical framework suited for the study of the polarizability of Lipschitz domains. Several aspects of polarizabilities and their representing measures are clarified, including limiting behavior both when approaching the support of the measure and when deforming smooth domains into a non-smooth domain. The success of the mathematical theory is achieved through symmetrization arguments for layer potentials.
  • An integral equation based scheme is presented for the fast and accurate computation of effective conductivities of two-component checkerboard-like composites with complicated unit cells at very high contrast ratios. The scheme extends recent work on multi-component checkerboards at medium contrast ratios. General improvement include the simplification of a long-range preconditioner, the use of a banded solver, and a more efficient placement of quadrature points. This, together with a reduction in the number of unknowns, allows for a substantial increase in achievable accuracy as well as in tractable system size. Results, accurate to at least nine digits, are obtained for random checkerboards with over a million squares in the unit cell at contrast ratio 10^6. Furthermore, the scheme is flexible enough to handle complex valued conductivities and, using a homotopy method, purely negative contrast ratios. Examples of the accurate computation of resonant spectra are given.
  • We investigate the optical properties of periodic composites containing metamaterial inclusions in a normal material matrix. We consider the case where these inclusions have sharp corners, and following Hetherington and Thorpe, use analytic results to argue that it is then possible to deduce the shape of the corner (its included angle) by measurements of the absorptance of such composites when the scale size of the inclusions and period cell is much finer than the wavelength. These analytic arguments are supported by highly accurate numerical results for the effective permittivity function of such composites as a function of the permittivity ratio of inclusions to matrix. The results show that this function has a continuous spectral component with limits independent of the area fraction of inclusions, and with the same limits for both square and staggered square arrays.