Johan Helsing

Johan Helsing, Anders Karlsson

A high-order accurate, explicit kernel-split, panel-based, Fourier-Nyström 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.

Johan Helsing

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.

Johan Helsing, Anders Holst

The incorporation of analytical kernel information is exploited in the construction of Nyström 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.

Johan Helsing, Anders Karlsson

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öm 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.

Johan Helsing

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.

Johan Helsing, Anders Karlsson

An integral equation-based numerical method for scattering from multi-dielectric cylinders is presented. Electromagnetic fields are represented via layer potentials in terms of surface densities with physical interpretations. The existence of null-field representations then adds superior flexibility to the modeling. Local representations are used for fast field evaluation at points away from their sources. Partially global representations, constructed as to reduce the strength of kernel singularities, are used for near-evaluations. A mix of local- and partially global representations is also used to derive the system of integral equations from which the physical densities are solved. Unique solvability is proven for the special case of scattering from a homogeneous cylinder under rather general conditions. High achievable accuracy is demonstrated for several examples found in the literature.

Johan Helsing, Anders Karlsson

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öm scheme. The resulting solver is used to accurately determine eigenwavenumbers and normalized electric eigenfields in the entire computational domain.