Source author record

Johan Helsing

Johan Helsing appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2023arXiv

Inverse problem for a planar conductivity inclusion

This paper concerns the inverse problem of determining a planar conductivity inclusion. Our aim is to analytically recover from the generalized polarization tensors (GPTs), which can be obtained from exterior measurements, a homogeneous inclusion with arbitrary constant conductivity. The primary outcome of recovering a homogeneous inclusion is an inversion formula in terms of the GPTs for conformal mapping coefficients associated with the inclusion. To prove the formula, we establish matrix factorizations for the GPTs.

preprint2022arXiv

Solving integral equations on piecewise smooth boundaries using the RCIP method: a tutorial

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.

preprint2021arXiv

Dirac integral equations for dielectric and plasmonic scattering

A new integral equation formulation is presented for the Maxwell transmission problem in Lipschitz domains. It builds on the Cauchy integral for the Dirac equation, is free from false eigenwavenumbers for a wider range of permittivities than other known formulations, can be used for magnetic materials, is applicable in both two and three dimensions, and does not suffer from any low-frequency breakdown. Numerical results for the two-dimensional version of the formulation, including examples featuring surface plasmon waves, demonstrate competitiveness relative to state-of-the-art integral formulations that are constrained to two dimensions. However, our Dirac integral equation performs equally well in three dimensions, as demonstrated in a companion paper.

preprint2020arXiv

An extended charge-current formulation of the electromagnetic transmission problem

A boundary integral equation formulation is presented for the electromagnetic transmission problem where an incident electromagnetic wave is scattered from a bounded dielectric object. The formulation provides unique solutions for all combinations of wavenumbers in the closed upper half-plane for which Maxwell's equations have a unique solution. This includes the challenging combination of a real positive wavenumber in the outer region and an imaginary wavenumber inside the object. The formulation, or variants thereof, is particularly suitable for numerical field evaluations as confirmed by examples involving both smooth and non-smooth objects.

preprint2016arXiv

An explicit kernel-split panel-based Nyström scheme for integral equations on axially symmetric surfaces

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.

preprint2016arXiv

Classification of spectra of the Neumann-Poincaré operator on planar domains with corners by resonance

We study spectral properties of the Neumann-Poincaré 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öm method which makes it possible to construct high-order convergent discretizations of the Neumann-Poincaré 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.

preprint2014arXiv

Variants of an explicit kernel-split panel-based Nyström discretization scheme for Helmholtz boundary value problems

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.

preprint2013arXiv

A higher-order singularity subtraction technique for the discretization of singular integral operators on curved surfaces

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.

preprint2012arXiv

An accurate boundary value problem solver applied to scattering from cylinders with corners

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.

preprint2012arXiv

On the polarizability and capacitance of the cube

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.

preprint2011arXiv

Spectral super-resolution in metamaterial composites

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.

preprint2011arXiv

The effective conductivity of arrays of squares: large random unit cells and extreme contrast ratios

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.