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)在电磁波被介电障碍物散射问题中的应用(可视为时谐麦克斯韦方程组中Lippmann-Schwinger型体积积分方程的一种简单离散化)的启发,我们分析了具有强奇异核的卷积算子的类似离散化方法。对于包括一维有限希尔伯特变换以及二维和三维中DDA所用麦克斯韦体积积分算子主部在内的一类核函数,我们证明该方法虽然不符合已知的投影方法框架,但可视为无限块Toeplitz矩阵的有限截断方法。该矩阵的符号由一个非绝对收敛的傅里叶级数给出。我们利用Ewald方法得到该符号的指数快速收敛级数表达式,并证明其为有界函数,从而能够描述该矩阵的谱和数值范围。结果表明,该数值范围包含积分算子的数值范围,但在某些情况下严格更大。在这些情况下,离散化方法无法提供谱精确的近似;尽管该方法在谱参数λ的较大范围内是稳定的,但存在某些λ值,使得奇异积分方程适定,而离散化方法却是不稳定的。