Charles L. Epstein

Charles L. Epstein, Leslie Greengard, Michael O'Neil

In this paper, we develop a new integral representation for the solution of the time harmonic Maxwell equations in media with piecewise constant dielectric permittivity and magnetic permeability in R^3. This representation leads to a coupled system of Fredholm integral equations of the second kind for four scalar densities supported on the material interface. Like the classical Muller equation, it has no spurious resonances. Unlike the classical approach, however, the representation does not suffer from low frequency breakdown. We illustrate the performance of the method with numerical examples.

Charles L. Epstein, Leslie Greengard

In this paper, we develop a new representation for outgoing solutions to the time harmonic Maxwell equations in unbounded domains in $\bbR^3.$ This representation leads to a Fredholm integral equation of the second kind for solving the problem of scattering from a perfect conductor, which does not suffer from spurious resonances or low frequency breakdown, although it requires the inversion of the scalar surface Laplacian on the domain boundary. In the course of our analysis, we give a new proof of the existence of non-trivial families of time harmonic solutions with vanishing normal components that arise when the boundary of the domain is not simply connected. We refer to these as $k$-Neumann fields, since they generalize, to non-zero wave numbers, the classical harmonic Neumann fields. The existence of $k$-harmonic fields was established earlier by Kress.

Charles L. Epstein, Leslie Greengard, Andreas Klöckner

In a recently developed quadrature method (quadrature by expansion or QBX), it was demonstrated that weakly singular or singular layer potentials can be evaluated rapidly and accurately on surface by making use of local expansions about carefully chosen off-surface points. In this paper, we derive estimates for the rate of convergence of these local expansions, providing the analytic foundation for the QBX method. The estimates may also be of mathematical interest, particularly for microlocal or asymptotic analysis in potential theory.

Charles L. Epstein, Leslie Greengard, Michael O'Neil

The generalized Debye source representation of time-harmonic electromagnetic fields yields well-conditioned second-kind integral equations for a variety of boundary value problems, including the problems of scattering from perfect electric conductors and dielectric bodies. Furthermore, these representations, and resulting integral equations, are fully stable in the static limit as $ω\to 0$ in multiply connected geometries. In this paper, we present the first high-order accurate solver based on this representation for bodies of revolution. The resulting solver uses a Nyström discretization of a one-dimensional generating curve and high-order integral equation methods for applying and inverting surface differentials. The accuracy and speed of the solvers are demonstrated in several numerical examples.

Charles L. Epstein, Zydrunas Gimbutas, Leslie Greengard et al.

A classical problem in electromagnetics concerns the representation of the electric and magnetic fields in the low-frequency or static regime, where topology plays a fundamental role. For multiply connected conductors, at zero frequency the standard boundary conditions on the tangential components of the magnetic field do not uniquely determine the vector potential. We describe a (gauge-invariant) consistency condition that overcomes this non-uniqueness and resolves a longstanding difficulty in inverting the magnetic field integral equation.

Charles L. Epstein, Michael O'Neil

We introduce an arbitrary order, computationally efficient method to smooth corners on curves in the plane, as well as edges and vertices on surfaces in $\mathbb R^3$. The method is local, only modifying the original surface in a neighborhood of the geometric singularity, and preserves desirable features like convexity and symmetry. The smoothness of the final surface is an explicit parameter in the method, and the bandlimit of the smoothed surface is proportional to its smoothness. Several numerical examples are provided in the context of acoustic scattering. In particular, we compare scattered fields from smoothed geometries in two dimensions with those from polygonal domains. We observe that significant reductions in computational cost can be obtained if merely approximate solutions are desired in the near- or far-field. Provided that it is sub-wavelength, the error of the scattered field is proportional to the size of the geometry that is modified.

Charles L. Epstein, Jon Wilkening

We study the classical Kimura diffusion operator defined on the n-simplex, $$L^{Kim}=\sum_{1\leq i,j\leq n+1}x_ix_j\partial_{x_i}\partial_{x_j}$$ We give novel constructions for the basis of eigenpolynomials, and the solution to the inhomogeneous Dirichlet problem, which are well adapted to numerical applications. Our solution of the Dirichlet problem is quite explicit and provides a precise description of the singularities that arise along the boundary.

Charles L. Epstein, Yoichiro Mori, Han Zhou

We study a bulk-surface coupled Laplace system involving an embedded open boundary. The problem is reformulated as an integro-differential equation using boundary integral representations, for which we establish existence and uniqueness of the solution. A Wiener-Hopf technique is employed to study the solution regularity and derive asymptotic expressions for the edge singularity. Building on these results, we develop a finite element method that incorporates the singularity structure and provide a rigorous error analysis. Numerical experiments confirm the theoretical convergence rates.

Alex H. Barnett, Charles L. Epstein, Leslie Greengard et al.

We study the stability properties of explicit marching schemes for second-kind Volterra integral equations that arise when solving boundary value problems for the heat equation by means of potential theory. It is well known that explicit finite difference or finite element schemes for the heat equation are stable only if the time step $Δt$ is of the order $O(Δx^2)$, where $Δx$ is the finest spatial grid spacing. In contrast, for the Dirichlet and Neumann problems on the unit ball in all dimensions $d\ge 1$, we show that the simplest Volterra marching scheme, i.e., the forward Euler scheme, is unconditionally stable. Our proof is based on an explicit spectral radius bound of the marching matrix, leading to an estimate that an $L^2$-norm of the solution to the integral equation is bounded by $c_dT^{d/2}$ times the norm of the right hand side. For the Robin problem on the half space in any dimension, with constant Robin (heat transfer) coefficient $κ$, we exhibit a constant $C$ such that the forward Euler scheme is stable if $Δt < C/κ^2$, independent of any spatial discretization. This relies on new lower bounds on the spectrum of real symmetric Toeplitz matrices defined by convex sequences. Finally, we show that the forward Euler scheme is unconditionally stable for the Dirichlet problem on any smooth convex domain in any dimension, in $L^\infty$-norm.

Charles L. Epstein, Leslie Greengard, Thomas Hagstrom

We give a principled approach for the selection of a boundary integral, retarded potential representation for the solution of scattering problems for the wave equation in an exterior domain.