Research Article  Open Access
A Fast Algorithm for Electromagnetic Scattering from OneDimensional Rough Surface
Abstract
In this paper, the Adaptive Modified Characteristic Basis Function Method (AMCBFM) is proposed for quickly simulating electromagnetic scattering from a onedimensional perfectly electric conductor (PEC) rough surface. Similar to the traditional characteristic basis function method (CBFM), FoldyLax multiple scattering equations are applied in order to construct the characteristic basis functions (CBFs). However, the CBFs of the AMCBFM are different from those of the CBFM. In the AMCBFM, the coefficients of the CBFs are first defined. Then, the coefficients and the CBFs are used to structure the total current, which is used to represent the induced current along the rough surface. Moreover, a current criterion is defined to adaptively halt the order of the CBFs. The validity and efficiency of the AMCBFM are assessed by comparing the numerical results of the AMCBFM with the method of moments (MoM). The AMCBFM can effectively reduce the size of the matrix, and it costs less than half the CPU time used by the MoM. Moreover, by comparing it with the traditional CBFM, the AMCBFM can guarantee the accuracy, reduce the number of iterations, and achieve better convergence performance than the CBFM does. The second order of the CBFs is set in the CBFM. Additionally, the first order of the CBFs of the AMCBFM alone is sufficient for this result.
1. Introduction
Electromagnetic (EM) scattering from rough surface has been widely studied and applied in the research areas of marine communication, target detection, and stealth technology [1–5]. The method of moments (MoM) [6–8] has high accuracy and is widely used to simulate the EM scattering from rough surface. However, the dimensions of the solution matrix of the MoM are , where represents the number of unknowns. When solving this problem using the MoM, one has to simulate a sufficiently long rough surface, which results in many unknowns and timeconsumption. Several numerical methods have been developed, such as the fast multipole method (FMM) [9] and the characteristic basis function method (CBFM) [10–12]. In the FMM, an additional algorithm is used to deal with the kernel function. In addition, the iterative method is applied to solve the matrix equation. However, the FMM is an iterative method, which is limited by the convergence of the matrix equation. In the CBFM, the FoldyLax multipath scattering equations [13] are applied in order to structure the characteristic basis functions (CBFs). The selfinteraction is first considered in order to generate the primary characteristic basis functions (PCBFs). Moreover, the mutual coupling effect of the subcells is considered in order to build the secondary characteristic basis functions (SCBFs). In previous work of our teams [14], the CBFM is used to simulate the EM scattering from a dielectric rough surface. In this work, the orders of the SCBFs were set by us and were determined through several tries. In the previous work of others [15], the following criterion was used to determine the order of the SCBFs:where is the impedance matrix, is the known vector that is excited by the incident wave, represents the induced current that is obtained by the CBFM, and is the L2 norm. From (1), it is readily seen that must be refilled and the matrix vector product is calculated in order to determine the order of the SCBFs. Therefore, determining the order of the SCBFs using (1) will be timeconsuming and require more memory.
In this paper, the Adaptive Modified Characteristic Basis Function Method (AMCBFM) [16] with a new current criterion is proposed. Similar to the traditional CBFM, the FoldyLax multipath scattering equations are also used to structure the CBFs. The PCBFs, such as the CBFM, are first obtained. However, the coefficients of the PCBFs are defined. Then, the primary total current of each cell is obtained by combining the coefficients and the PCBFs. Then, the coefficients are taken with the PCBFs in order to get the first order of the SCBFs. The coefficients of the first order of the SCBFs and the first total current of each cell are obtained. Like the primary total current, the first total current is used to get the second order of the SCBFs. Higher orders of the SCBFs are similarly obtained. Moreover, the new current criterion that is defined in this paper is used to adaptively halt the order of the SCBFs.
The remainder of this paper is organized as follows. In Section 2, the theoretical formulations of the AMCBFM are exhibited in detail. In Section 3, the validity and efficiency of the AMCBFM are assessed and compared with the MoM and the CBFM. Section 4 presents a summary of this paper and further research plans on this topic.
2. Theoretical Formulations
2.1. EM Scattering from a Rough Surface Using the MoM
Assume that there is a tapered plane wave incident upon a PEC rough surface, as shown in Figure 1, where is a Gaussian rough surface with a simulated length . and are the incident angle and the scattered angle, respectively. Space with the relative permittivity and the relative permeability is the incident space. Here, space is assumed to be the free space; i.e., , and . represents the scattered wave in space . The position vector is .
For a Horizontal (H) polarized incident wave, the surface integral equation for this scattering problem is [17]where denotes the interface of the PEC rough surface, is the unit normal vector of the rough surface, is the twodimensional Green’s function in the space , and represents the total wave including the incident wave and the scattered wave in space .
For the Vertical (V) case, the surface integral equation is as follows [17].By using the MoM with the pulse basis functions and the point matching technique, the integral equations (2) and (3) can be discretized into the following matrix equations.The elements of the impedance matrix and in detail arewhere is the wavenumbers in space and . The subscripts and are the field point and source point along the rough surface, respectively. The elements of the right vector and are , and the left vectors and are the unknown current vectors.
2.2. Theoretical Formulations of AMCBFM
Analyzing (4) and (5), the two equations can be rewritten aswhere is the impedance matrix with the dimensions , the dimensions of the right vector are , the dimensions of the unknown current vector are , and is the number of the discrete segments along the rough surface.
In the AMCBFM, the segments of the rough surface are first divided into cells; every cell contains segments, and . Equation (8) can be rewritten as [18]where is the mutual impedance matrix of the cell, and the cell with the dimensions . represents the total current of the cell, which represents the induced current along the cell. indicates the incident wave upon the cell.
(1) The Primary Total Current . For the cell, only the selfinteraction is considered to generate the PCBFs such as the CBFMwhere is the selfinteraction impedance matrix of the cell, and and are the PCBFs and the incident wave vector of the cell, respectively. Using the lower upper (LU) decomposition technique, one can obtain the PCBFs . The lower and the upper triangular matrices of are obtained and saved.
Assuming the total current is only from the PCBFs, it is called the primary total current , and it is calculated aswhere is the coefficient of the PCBFs.
Taking (11) into (9), one can get the following.Then, by multiplying (12) by , where is the conjugate transpose of , one can obtain the following.Using the LU decomposition technique, the PCBFs coefficients and the primary total current are acquired. It is obviously found that the mutual coupling effect of the subcells is neglected, and thus is not accurate. However, the coefficients can be used to get the SCBFs.
(2) The First Total Current . Considering the mutual coupling effect of the subcells when building the SCBFs, for the cell, the first order of the SCBFs is defined as where is the first order SCBFs. In addition, the coefficients are obtained using (13).
The first total current is defined aswhere and are the coefficients of the first order of SCBFs.
Taking (15) into (9), one can get the following.It is readily seen that (16) has 2M unknowns and M equations. In order to solve this problem, (16) is multiplied by and , which results in the following.The CBFs coefficients and the secondary total current are obtained.
(3) The Total Current . Similar to the first order of SCBFs, the order of the SCBFs is defined aswhere is the order of the SCBFs, and is the coefficient of .
The total current is defined as follows.
(4) The Current Criterion. The current criterion that determines whether the higherorder SCBFs are needed is defined as [16]where is the L2 norm. It is readily seen that the order of the SCBFs can be adaptively halted.
(5) The Percentage Error. The percentage error of the AMCBFM or the CBFM is defined aswhere represents the induced current that is obtained using the AMCBFM or the CBFM, denotes the induced current that is calculated using the MoM, and is the L2 norm.
3. Numerical Simulations and Discussions
The validity and efficiency of the AMCBFM are first assessed by comparing its results with the results that were obtained from the CBFM and the MoM. Both the H and V polarized incident waves are considered. In Figure 2, the bistatic scattering coefficients that are obtained using these three methods are shown. The induced currents along the rough surface that are obtained by the AMCBFM and the MoM are given in Figure 3. The frequency of the incident wave is , and the incident angle is . The rms height and the correlation length of the simulated rough surface are and , respectively. The length of the simulated Gaussian rough surface is with a discrete interval of . In other words, there are segments along the rough surface. Thus, the dimension of the MoM matrix is . The rough surface is equally divided into 4 cells. Thus, the dimensions of the CBFM matrix and the AMCBFM matrix are only . It is clearly determined that the dimensions of the dealing matrix have been significantly reduced in the AMCBFM and the CBFM. The orders of the SCBFs of the CBFM are set as . The current criterion of the AMCBFM is and the halted order of the SCBFs is . From Figures 2 and 3, it is readily seen that the AMCBFM and the CBFM agree well with the MoM in both the H and V polarization. However, the halted order of SCBFs of AMCBFM is through the criterion . The orders of the SCBFs of the CBFM are set as .
(a)
(b)
(a)
(b)
Moreover, in Table 1, the CPU times of the AMCBFM, the CBFM, and the MoM are shown. It is obviously found that the AMCBFM costs less than half the CPU time of the MoM. In addition, Table 1 compares the percentage errors of the AMCBFM and the CBFM. It is readily seen that the AMCBFM with the order is more efficient than the CBFM with the order . The results are implemented on a computer with an Intel Core i76700 processor, 4 GB of memory, and the Microsoft Windows 10 operating system.

One may note that only one sample of the rough surface is used in Figures 2 and 3. In order to further test the efficiency of the proposed algorithm, the results obtained with and are shown in Table 2. These results are based on the rough surface with a simulated length of , a discrete interval of , and an incident angle of . The CPU time of the MoM is 12.995 seconds for H polarization and 13.01 seconds in the V case. In the AMCBFM, the rough surface is equally divided into 4 cells. Thus, the dimensions of the AMCBFM matrix are only . The current criteria of AMCBFM are set as , , , , and . From Table 2, it is readily seen that the order of the SCBFs can be adaptively halted in the AMCBFM. Furthermore, the current criterion is sufficient for retaining the accuracy of the induced current. By comparing the CPU times between the AMCBFM and the MoM, it is easily determined that the AMCBFM with costs less than half the CPU time of the MoM.

4. Conclusions
In this paper, the AMCBFM was used to simulate the EM scattering from a Gaussian rough surface. Good agreement between the AMCBFM and the MoM was found. In addition, the AMCBFM costs less than half the CPU time of the MoM. Furthermore, compared with the traditional CBFM, the lower order of the CBFs of the AMCBFM is sufficient for guaranteeing the accuracy. Note that the scattering model is only a 2D model. In the future, this method will be expanded to solve the more realistic 3D scattering problem. In addition, this method can also be applied to calculate the dielectric rough surface with a cylinder located above it. Finally, the AMCBFM can be combined with the compressing sensing technology in order to solve the EM scattering over a frequency band.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (Grant Nos. 61501004 and 61722101), the Natural Science Foundation of Anhui Province (Grant Nos. 1608085QF141 and 1608085MF135), and the Provincial Program of Natural Science of Anhui Higher Education (Grant No. KJ2015A073).
References
 R.W. Xu, L.X. Guo, and R. Wang, “An iterative analytic  Numerical method for scattering from a target buried beneath a rough surface,” Chinese Physics B, vol. 23, no. 11, pp. 1674–1056, 2014. View at: Google Scholar
 S. Yang, J. Li, and L. Guo, “Electromagnetic scattering characteristics of ablation rough surface in plasma sheath,” International Symposium on Antennas, pp. 435–438, 2017. View at: Google Scholar
 J. Ma, L.X. Guo, and A.Q. Wang, “Study of MPI based on parallel MOM on PC clusters for EMbeam scattering by 2D PEC rough surfaces,” Chinese Physics B, vol. 18, no. 8, pp. 3431–3437, 2009. View at: Publisher Site  Google Scholar
 R. Bao, A.Q. Wang, and Z.X. Huang, “Fast simulations of electromagnetic scattering from rough surface over a frequency band using asymptotic waveform evaluation technique: horizontal polarization,” Journal of Electromagnetic Waves and Applications, vol. 32, no. 11, pp. 1379–1388, 2018. View at: Publisher Site  Google Scholar
 R. Wang, L. X. Guo, and A. Q. Wang, “Investigation of electromagnetic scattering interaction between the buried target and the rough surface in different types of soil,” Acta Physica Sinica, vol. 59, no. 5, pp. 3179–3186, 2010. View at: Google Scholar
 S. R. Chai and L. X. Guo, “A new fast algorithm based on compressive sensing for composite electromagnetic back scattering from a 2D ship located on a 1D rough sea surface,” Acta Physica Sinica, vol. 64, no. 6, Article ID 060301, 2015. View at: Google Scholar
 Y. Chen, D. Li, J. Hu, and J. Lee, “A nonconformal surface integral equation for electromagnetic scattering by multiscale conducting objects,” IEEE Journal on Multiscale and Multiphysics Computational Techniques, vol. 3, pp. 225–234, 2018. View at: Publisher Site  Google Scholar
 R. F. Harrington, Field Computation by Moment Methods, WileyIEEE Press, 1993. View at: Publisher Site
 V. Rokhlin, “Rapid solution of integral equations of classical potential theory,” Journal of Computational Physics, vol. 60, no. 2, pp. 187–207, 1985. View at: Publisher Site  Google Scholar  MathSciNet
 V. V. S. Prakash and R. Mittra, “Characteristic basis function method: a new technique for efficient solution of method of moments matrix equations,” Microwave and Optical Technology Letters, vol. 36, no. 2, pp. 95–100, 2003. View at: Publisher Site  Google Scholar
 C.S. Park, Y.R. Jeong, I.P. Hong, and J.G. Yook, “Block size optimization of CBFM for scattering problems,” IEEE Transactions on Antennas and Propagation, vol. 66, no. 10, pp. 5370–5377, 2018. View at: Publisher Site  Google Scholar
 C. Li, Y. Sun, G. Wang, and Z. Wang, “Fast analysis for largescale electromagnetic scattering problems by a hybrid approach,” Journal of Electromagnetic Waves and Applications, vol. 31, no. 8, pp. 808–819, 2017. View at: Publisher Site  Google Scholar
 Y. Sun, C. Chan, R. Mittra, and L. Tsang, “Characteristic basis function method for solving large problems arising in dense medium scattering,” in Proceedings of the IEEE Antennas Propag Soc Int Symp, pp. 1068–1071, Columbus, OH, USA, 2003. View at: Publisher Site  Google Scholar
 J. J. Wang, A. Q. Wang, T. Z. Jiang, and Z. X. Huang, “A fast algorithm for electromagnetic scattering from onedimensional dielectric rough surface,” Progress in Electromagnetic Research Symposium, pp. 362–365, 2018. View at: Publisher Site  Google Scholar
 Y. F. Sun, Y. Zhang, S. J. Xu, and X. Q. Chen, “EM scattering analysis of 2D multiple conducting cylinders using characteristic basis function method,” Chiness Journal of Radio Science, vol. 21, no. 2, pp. 229–232, 2006. View at: Google Scholar
 G. D. Han and C. Q. Gu, “Adaptively modified characteristic basis function method for electromagnetic scattering study,” Journal of Electronics & Information Technology, vol. 30, no. 10, pp. 2364–2368, 2008. View at: Google Scholar
 L. Tsang and J. A. Kong, Scattering of Electromagnetic Waves: Numerical Simulations, John Wiley & Sons, New York, NY, USA, 2001. View at: MathSciNet
 G. Han and C. Gu, “A hybrid QR factorization with dualMGS and adaptively modified characteristic basis function method for electromagnetic scattering analysis,” Microwave and Optical Technology Letters, vol. 49, no. 11, pp. 2879–2883, 2007. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Jing Jing Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.