Oscar P. Bruno,Agustin G. Fernandez-Lado
Abstract:This paper presents a full-spectrum Green function methodology (which is valid, in particular, at and around Wood-anomaly frequencies) for evaluation of scattering by periodic arrays of cylinders of arbitrary cross section-with application to wire gratings, particle arrays and reflectarrays and, indeed, general arrays of conducting or dielectric bounded obstacles under both TE and TM polarized illumination. The proposed method, which, for definiteness is demonstrated here for arrays of perfectly conducting particles under TE polarization, is based on use of the shifted Green-function method introduced in the recent contribution (Bruno and Delourme, Jour. Computat. Phys. pp. 262--290 (2014)). A certain infinite term arises at Wood anomalies for the cylinder-array problems considered here that is not present in the previous rough-surface case. As shown in this paper, these infinite terms can be treated via an application of ideas related to the Woodbury-Sherman-Morrison formulae. The resulting approach, which is applicable to general arrays of obstacles even at and around Wood-anomaly frequencies, exhibits fast convergence and high accuracies. For example, a few hundreds of milliseconds suffice for the proposed approach to evaluate solutions throughout the resonance region (wavelengths comparable to the period and cylinder sizes) with full single-precision accuracy.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: **How to solve the scattering problem caused by periodically arranged cylinders efficiently and accurately near the Wood abnormal frequency?** Specifically, the author proposes a new Green's function method, which can converge rapidly and provide high - precision solutions near and at the Wood abnormal frequency.
### Problem Background
When dealing with the scattering problems of periodically arranged cylinders (such as gratings, particle arrays, and reflective arrays), the classical quasi - periodic Green's function will fail or converge extremely slowly at the Wood abnormal frequency. The Wood abnormal frequency refers to the situation when the wave vector of the incident wave matches a certain diffraction order of the obstacle array. At this time, the propagation directions of some modes are parallel to the obstacle array, causing the denominator of the classical Green's function to approach zero, thus leading to difficulties in numerical calculations.
### Main Contributions of the Paper
1. **Introduction of the Shifted Green's Function**:
- By introducing the shifted Green's function, the author solves the divergence problem of the classical Green's function near the Wood abnormal frequency. The shifted Green's function achieves rapid convergence by adding several new poles outside the physical propagation domain, and can maintain good convergence and precision even near the Wood abnormal frequency.
2. **Treatment of the Singular Behavior at the Wood Abnormal Frequency**:
- Near the Wood abnormal frequency, although the shifted Green's function exists, the denominators of some of its terms will approach zero. Using the idea of the Woodbury - Sherman - Morrison formula, the author separates these singular terms and deals with these problems by the method of finite - rank operators, thus obtaining an effective solution method that can be used at and near the Wood abnormal frequency.
3. **Efficient Numerical Implementation**:
- The author proposes a numerical method based on the half - space quasi - periodic Green's function, the smoothing window function method, and the high - order singular integral quadrature formula, achieving fast and high - precision solutions. Experimental results show that this method can obtain single - precision solutions within milliseconds, and can maintain high precision even in the resonance region (when the wavelength is comparable to the period and the cylinder size).
### Summary of Mathematical Formulas
- **Shifted Green's Function**:
\[
G_{j,k}(X,Y)=\frac{i}{4}\sum_{\ell = 0}^{j}(- 1)^{\ell}\binom{j}{\ell}H_{0}^{(1)}\left(k\sqrt{X^{2}+(Y + \ell h)^{2}}\right)
\]
where \(H_{0}^{(1)}\) is the first - kind zero - order Hankel function.
- **Quasi - periodic Green's Function**:
\[
G_{q}^{j,k}(X,Y)=\sum_{n\in\mathbb{Z}}e^{-i\alpha_{n}L}G_{j,k}(X + nL,Y)
\]
- **Rayleigh Series Operator**:
\[
I_{n}^{k}[\psi]=\int_{\partial D_{0}}\left(\frac{\partial}{\partial\nu(r')}\left(e^{-i\alpha_{n}(k)x'-i\beta_{n}(k)y'}\right)-i\gamma e^{-i\alpha_{n}(k)x'-i\beta_{n}(k)y'}\right)\psi(r')dS(r')
\]
- **Mixed Integral Equation/Rayleigh Series Operator Formula**:
\[
\frac{1}{2}\psi(r)+\int_{\partial D_{0}}\left(\frac{\partial G_{q}^{j,k}}{\partial\nu(r')}-i\gamma G_{q}^{j,k}\right)(r - r')\psi(r')dS(r')-\frac{i}{2L}\sum_{n\in\mathbb{Z}}\frac{\s