Related papers: Fast Multipole Boundary Element Method for Three D…
This work describes a new version of the Fast Multipole Method for summing pairwise particle interactions that arise from discretizing integral transforms and convolutions on the sphere. The kernel approximations use barycentric Lagrange…
Boundary integral equations and Nystrom discretization provide a powerful tool for the solution of Laplace and Helmholtz boundary value problems. However, often a weakly-singular kernel arises, in which case specialized quadratures that…
We consider the efficient numerical solution of the three-dimensional wave equation with Neumann boundary conditions via time-domain boundary integral equations. A space-time Galerkin method with $C^\infty$-smooth, compactly supported basis…
We introduce TERMS, an open-source Fortran program to simulate near-field and far-field optical properties of clusters of particles. The program solves rigorously the Maxwell equations via the superposition T-matrix method, where incident…
We present algorithms for computing weakly singular and near-singular integrals arising when solving the 3D Helmholtz equation with curved boundary elements. These are based on the computation of the preimage of the singularity in the…
A new method for numerical solving of boundary problem for ordinary differential equations with slowly varying coefficients which is aimed at better representation of solutions in the regions of their rapid oscillations or exponential…
In this paper, we study an inverse scattering problem at fixed energy on three-dimensional asymptotically hyperbolic St{\"a}ckel manifolds having the topology of toric cylinders and satisfying the Robertson condition. On these manifolds the…
In boundary element methods (BEM) in $\mathbb{R}^3$, matrix elements and right hand sides are typically computed via analytical or numerical quadrature of the layer potential multiplied by some function over line, triangle and tetrahedral…
In the present paper we describe a method for solving inverse problems for the Helmholtz equation in radially-symmetric domains given multi-frequency data. Our approach is based on the construction of suitable trace formulas which relate…
Integral-equation-based fast direct solvers for electromagnetic scattering can substantially reduce computational costs, especially in the presence of multiple excitations. We recently proposed a new high-frequency fast direct solver…
We develop and implement a new mathematical and computational framework for designing photonic elements with one or more high-$Q$ scattering resonances. The approach relies on solving for the poles of the scattering matrix, which…
We use the tridiagonal representation approach to solve the radial Schr\"odinger equation for the continuum scattering states of the Kratzer potential. We do the same for a radial power-law potential with inverse-square and inverse-cube…
This paper is concerned with the error estimation of the fast multipole method (FMM) for scattering problems in 2-D. The FMM error is caused by truncating Graf's addition theorem in each step of the algorithm, including two expansions and…
The scattering of electromagnetic waves by an obstacle is analyzed through a set of partial differential equations combining the Maxwell's model with the mechanics of fluids. Solitary type EM waves, having compact support, may easily be…
Consider the scattering of a time-harmonic elastic plane wave by a bi-periodic rigid surface. The displacement of elastic wave motion is modeled by the three-dimensional Navier equation in an open domain above the surface. Based on the…
We propose in this work a fast numerical algorithm for solving the equation of radiative transfer (ERT) in isotropic media. The algorithm has two steps. In the first step, we derive an integral equation for the angularly averaged ERT…
In order to generate initial data for nonlinear relativistic simulations, one needs to solve the Einstein constraints, which can be cast into a coupled set of nonlinear elliptic equations. Here we present an approach for solving these…
Advanced computational tools that describe the interaction of electrons with structured nanophotonic materials are crucial for theoretical predictions, specific design tasks, and the interpretation of experimental results. These tools open…
We propose and analyse a hybrid numerical-asymptotic $hp$ boundary element method for time-harmonic scattering of an incident plane wave by an arbitrary collinear array of sound-soft two-dimensional screens. Our method uses an approximation…
A class of algorithms for the solution of discrete material optimization problems in electromagnetic applications is discussed. The idea behind the algorithm is similar to that of the sequential programming. However, in each major iteration…