Motivated by the discrete dipole approximation (DDA) for the scattering of electromagnetic waves by a dielectric obstacle that can be considered as a simple discretization of a Lippmann-Schwinger style volume integral equation for time-harmonic Maxwell equations, we analyze an analogous discretization of convolution operators with strongly singular kernels. For a class of kernel functions that includes the finite Hilbert transformation in 1D and the principal part of the Maxwell volume integral operator used for DDA in dimensions 2 and 3, we show that the method, which does not fit into known frameworks of projection methods, can nevertheless be considered as a finite section method for an infinite block Toeplitz matrix. The symbol of this matrix is given by a Fourier series that does not converge absolutely. We use Ewald's method to obtain an exponentially fast convergent series representation of this symbol and show that it is a bounded function, thereby allowing to describe the spectrum and the numerical range of the matrix. It turns out that this numerical range includes the numerical range of the integral operator, but that it is in some cases strictly larger. In these cases the discretization method does not provide a spectrally correct approximation, and while it is stable for a large range of the spectral parameter $\lambda$, there are values of $\lambda$ for which the singular integral equation is well posed, but the discretization method is unstable.
翻译:本文受到离散偶极子近似(DDA)的启发,该方法用于解决电磁波在介质障碍物中的散射问题。DDA被视作Lippmann-Schwinger体积积分方程的简单离散化,在时谐Maxwell方程中描述了体积积分算子。本文对一类核函数的离散化方法进行了分析,该核函数包括1D中的有限希尔伯特变换和2D和3D中用于DDA的Maxwell体积积分算子的主部分。我们证明了虽然该方法不符合投影方法的已知框架,但仍可以将其视为无限块Toeplitz矩阵的有限剖分法。此矩阵的符号由傅里叶级数给出,该级数不绝对收敛。我们使用Ewald方法以指数速度收敛的级数表示该符号,并证明了它是有界函数,从而得出了矩阵的谱和数值范围。结果表明,此数值范围包含了积分算子的数值范围,但在某些情况下,它会严格大于其数值范围。在这些情况下,离散化方法不能提供谱正确的近似,虽然它在大范围的谱参数λ上是稳定的,但在某些λ值上,奇异积分方程是良定义的,而离散化方法则不稳定。