Full Wave Modeling of Array Structures Using Generalized Single-Source Tangential Equivalence Principle Algorithm

Shao, Hanru; Dong, Jianfeng; Hu, Jun

doi:https://doi.org/10.1155/2017/4903906

International Journal of Antennas and Propagation

On this page

Abstract Introduction Conclusions Conflicts of Interest Acknowledgments References Copyright Related Articles

Special Issue

Advanced Modeling and Simulation Methods for Multiphysics and Multiscale Problems

View this Special Issue

Research Article | Open Access

Volume 2017 | Article ID 4903906 | https://doi.org/10.1155/2017/4903906

Full Wave Modeling of Array Structures Using Generalized Single-Source Tangential Equivalence Principle Algorithm

Hanru Shao,¹Jianfeng Dong,¹and Jun Hu²

Academic Editor: Yumao Wu

Received28 Apr 2017

Accepted28 Jun 2017

Published06 Aug 2017

Abstract

The generalized single-source tangential equivalence principle algorithm (GSST-EPA) is proposed to solve the electromagnetic scattering of array structures with very small distance or even connected elements. For the traditional EPA scheme, it is difficult to deal with those situations. Based on source reconstruction method (SRM), the new GSST-EPA can be used to solve this problem efficiently. This scheme can also reduce the number of unknowns by half via the extinction theorem and keep the good accuracy by using tangential field projection. Furthermore, the multilevel fast multipole algorithm (MLFMA) is used to accelerate the matrix-vector multiplication in the GSST-EPA. Several numerical results are given to demonstrate the accuracy and efficiency of the proposed method.

1. Introduction

The array structure has a great range of applications in practical engineering, such as antenna arrays, metamaterials, and microwave absorber. As a kind of multiscale modeling, the array structure involving large number of elements and fine components has been the challenge for traditional full wave methods [1]. Therefore, more advanced numerical methods with better efficiency and robustness need to be developed.

In recent years, many efforts have been made to address the above problems. Particularly, the full wave integral equation (IE) solvers have attracted much attention due to high accuracy, relatively small matrices, and ability to be accelerated by fast algorithms. The progress in IE solvers for array structures can be divided into two categories: one is macro basis function methods focusing on reducing the number of unknowns, for instance, the sub-entire-domain (SED) basis function method [2], the characteristic basis function (CBF) method [3], and the synthetic basis function (SBF) approach [4]. These methods are similar to each other except the way of generating macro basis function. The other one is domain decomposition methods based on integral equations, including equivalence principle algorithm (EPA) [5–7], linear embedding via Green’s operator (LEGO) [8], and generalized transition matrix (GTM) [9]. The objectives of these methods are similar by using surface equivalence theorem to transform the unknowns on original objects onto the equivalence surfaces enclosing the objects.

The EPA has been used to solve some multiscale problems. Several efforts have also been made to improve the accuracy and extend the application range. For example, the high-order field point sampling scheme [7] and tangential EPA (T-EPA) [10] have been proposed to improve the accuracy of scattering operator. The efficiency of scattering operator has been improved by hierarchical LU decomposition method [11]. Furthermore, the low-frequency problem [12] and time-domain models [13] have been solved using EPA method. The EPA is also combined with BOR to solve large-scale problems [14].

In this paper, the EPA is implemented to solve large-scale array structures and some important improvements are made to achieve better numerical flexibility and efficiency. Firstly, the single-source EPA is used to reduce the number of unknowns by half via extinction theorem, in which only electric currents or magnetic currents on the equivalence surfaces need to be solved [15]. Secondly, the tangential field projection scheme is applied on equivalence surface to improve the accuracy of single-source EPA, which is called single-source tangential EPA (T-EPA) [16]. Thirdly, the source reconstruction method (SRM) is proposed to solve the case that connected structures intercepted by equivalence surfaces [17]. Using this scheme, the EPA equations are still well-conditioned and only equivalence surface currents are the final unknowns to be solved. Finally, the above three improvements are combined together as a new scheme named generalized single-source tangential EPA (GSST-EPA). The GSST-EPA is more efficient and robust than traditional EPA for solving complicated array problems. To improve the ability of solving large-scale array structures further, the multilevel fast multipole algorithm (MLFMA) [18] is applied to accelerate the translation procedure which is the coupling between nonadjacent equivalence surfaces. Numerical results are shown to validate the accuracy and efficiency of the new method.

2. Equivalence Principle Algorithm

The basic idea of EPA is to transform the original problems into a new equivalent problem with scattering operator and translation operator. The electromagnetic scattering of nonconnected perfect electric conductor (PEC) objects in free space is considered to illustrate EPA. Each object is enclosed by an artificial equivalence surface ES (). The EPA scheme is based on the following two procedures: the scattering solution of objects via an equivalence surface and the couplings between objects via the equivalence surfaces, corresponding to the scattering operator and translation operator separately. Next, these two operators are defined in detail.

2.1. Scattering Operator

The first procedure can be derived from a PEC object characterized by equivalence surface with three steps as shown in Figure 1. The first step is outside-in propagation. The equivalence incident electric and magnetic currents on the equivalence surface can be obtained by using equivalence principle as follows:where is the unit outer normal vector on the ES . The currents and generate the original incident fields inside and zero field outside. Then the electric current on the PEC can be solved using method of moments (MoM) which is the second step asHere, and are the surface integral operators defined aswhere is or ; G₀ is Green’s function in free space. is the wave impedance in free space. P.V. stands for the Cauchy principle value integration. The last step is inside-out propagation. Once the PEC current is solved, the scattering electric and magnetic currents can be computed using equivalence principle as where is identity operator. The currents and generate the original scattering fields outside and zero field inside. Combining these three steps, the scattering of a PEC object via an equivalence surface can be characterized by the scattering operator aswhere

It can be seen that the scattering operator can capture the full physics rigorously. Therefore, the unknowns on the inside object are transferred to the unknowns on the equivalence surface. Because the inside objects always have complex structures and different materials, the number of unknowns on the equivalence surface can be much less than that on the object.

2.2. Translation Operator

The second procedure can be illustrated by the scattering of two PEC objects as shown in Figure 2. The field scattered by and produces new fields on ES . So these fields can be interpreted as the incident currents from the ES to ES and can be expressed by a translation operator , , as where Therefore, the scattered currents on ES can be written as For the generalized scattering of N nonadjacent PEC objects enclosed by N equivalence surfaces separately, the EPA equations can be written as

The matrix form can be obtained by expanding the currents with Rao–Wilton–Glisson (RWG) basis function and using Galerkin’s method to discretize the equation [19]. If the array structures have the same element, then each element can be enclosed by the same equivalence surface which leads to the same scattering operator. Therefore, only one scattering operator needs to be solved and stored which can reduce the solving time significantly.

3. Single-Source Equivalence Principle Algorithm

The single-source EPA is similar to original EPA scheme which is based on the above two procedures. The difference is that the equivalence magnetic currents on ES can be expressed by the equivalence electric currents on ES in the single-source EPA. Therefore, only the electric currents on the equivalence surfaces are the final unknowns, which can reduce the number of unknowns by half compared with original EPA.

According to extinction theorem, the equivalence incident currents and generate zero field outside the ES as shown in Figure 3(a):where ES denotes the outside surface of ES . The integral operators , are defined as (3), and is defined as is equivalence incident or scattering surface current or . Then the magnetic current can be expressed by electric current asSimilarly, the relation between equivalence scattering currents and can be established via extinction theorem as shown in Figure 3(b):where ES denotes the inside surface of ES . It can be seen that there are two types of expressions for the magnetic current and electric current, which are named as E-type and H-type. Therefore, the scattering operator of single-source EPA on ES also has two types and can be written as The definition of scattering operator on ES is similar to the above. The two types of translation operator from ES to ES are defined asFinally, the single-source EPA equations have two types of expressions as follows: Once the electric currents have been solved, the magnetic currents can be obtained using (15) or (16).

(a)

(b)

4. Single-Source Tangential Equivalence Principle Algorithm

In the single-source EPA, the equivalence surface currents are the rotated tangential projections of the fields onto the surface as follows:where is the known field. Therefore both scattering operator and translation operator have identity operator in (17) and (18). Numerical results have shown that the discretization of identity operator using conventional Galerkin’s scheme and common divergence-conforming RWG function leads to large error in the field projection procedure [10]. So the accuracy of single-source EPA needs to be improved. In this paper, a different field projection scheme is adopted to improve the accuracy by using surface integral to represent the tangential fields. As shown in Figures 3(a) and 3(b), the equivalent currents are generated from the tangential field projection on the equivalence surface ES by using Huygens equivalence principle as follows:By submitting (13) and (14) into (21) and (22) and submitting (15) and (16) into (23) and (24), there are four types of expressions for the incident current and scattering current separately, which can be denoted as TEE, TEH, THE, and THH. Numerical results have shown that and are more accurate than and when using RWG function as basis and testing functions [15]. It also has been proved that the THH type is the same as single-source EPA [16]. Therefore, only the TEH is used in this paper.

The single-source tangential EPA (T-EPA) based on the TEH type can be derived by submitting (14) into (21) and (16) into (23) as follows:Then the scattering operator of single-source T-EPA on ES can be written as and the translation operator should be changed as

It can be seen that both of scattering operator and translation operator have the combination of and operators instead of using identity operator. So when using conventional Galerkin’s scheme and RWG function to discretize the scattering and translation operators, this will lead to a more accuracy result [10]. Finally, the equations of single-source T-EPA are the same as (19) except for the expressions of scattering and translation operators.

5. The Single-Source Tangential Equivalence Principle Algorithm with Source Reconstruction Method (SRM)

In the practical engineering, the elements of array structures may be very close to each other. It is a challenge for EPA because the near fields between the equivalence surfaces and objects are hard to calculate with high accurate. The other situation is that the elements are connected with each other, which can be solved with tap basis method. But this method needs to increase the number of unknowns and deteriorate the conditioning of EPA equations. In this paper, the source reconstruction method (SRM) is combined with single-source T-EPA to deal with the above two challenges. To illustrate this method, a scattering problem of a PEC strip is considered. The strip is divided into two parts as shown in Figure 4. The left part PEC 1 is enclosed by ES 1, and the right part PEC 2 is enclosed by ES 2. is the current on PEC 1; is the current on PEC 2. The MoM equations can be generated asThe equivalence scattering currents on ES 1 can be obtained through inside-out propagation:where the operator has been defined in (26). Then substituting (29) into (31) yields is the incident plane wave on PEC 1 which is already known, while needs to be solved. Suppose that the equivalence scattering currents and on ES 2 are already known, then the can be got by using source reconstruction method [17] as follows:where denotes the pseudoinverse. In this paper, the pseudoinverse is computed by truncated SVD method. can be replaced by equivalence incident currents asHence, substituting (33) and (34) into (32), we havewhere and the definition of is the same as (17). Similarly, the equation for the equivalence scattering currents on ES 2 can be written as

Finally, the new EPA equations for connected structures can be obtained by combining (35) and (37):The coupling between PEC 1 and PEC 2 is calculated directly and there is no need to use buffer region. So the final unknowns are only equivalence surface currents by using the source reconstruction method. Equation (38) is simple and well-conditioned, which is similar to original EPA equation (10).

6. Generalized Single-Source Tangential Equivalence Principle Algorithm

By combing the above single-source EPA, tangential field projection, and source reconstruction method, the generalized single-source tangential equivalence principle algorithm (GSST-EPA) can be derived. Considering the scattering of connected objects, each object is enclosed by an equivalence surface as shown in Figure 5. Then the GSST-EPA equation can be written aswhereIt should be pointed out that although the GSST-EPA can be used to solve generalized object, the GSST-EPA is much more suited for array structures because is identical when the array element and equivalence surface are the same as each other. For the periodic array structures, and have the property of translation invariance which can be used to reduce the computing time and memory requirement. Therefore, the GSST-EPA is used to solve array structures in this paper.

To improve the ability of GSST-EPA to solve large-scale array structures, multilevel fast multipole algorithm (MLFMA) is applied in GSST-EPA. This is because the calculating of mutual coupling between two equivalence surfaces is divided into two cases: adjacent case and nonadjacent case; then the MLFMA in GSST-EPA is different from that in original EPA. Firstly, the adjacent case is the near interaction actually which is solved by SRM. So the matrix is calculated directly and stored. Secondly, the nonadjacent case is the far interaction which is solved by translation operator. So the translation matrix is implemented using MLFMA instead of MoM.

7. Numerical Results

In this part, several numerical examples are given to show the accuracy and efficiency of the GSST-EPA.

The first example is to demonstrate the accuracy of GSST-EPA. The scattering of two PEC spheres in free space is investigated. The radius of PEC sphere is 0.3 m. The distance between the centers of two spheres is 1.0 m in -axis. Each sphere is enclosed by equivalence surface separately, which is a rectangle surface with m³ as shown in Figure 6. The excitation is a polarized plane wave propagating into the negative direction at 300 MHz. The curvilinear Rao–Wilton–Glisson (CRWG) function [19] is used to discretize the spheres and equivalence surfaces. This case is solved by three kinds of methods separately: MoM, T-EPA [10], and GSST-EPA. The bistatic Radar Cross Section (RCS) results for the HH polarization are shown in Figure 7. It can be seen that the results agree with each other very well.

In the second example, the scattering of two PEC spheres is still investigated but changing the center distance from 1.0 m to 0.605 m as shown in Figure 8. In this case, it is difficult for the T-EPA scheme [10] to solve this problem accurately, because these spheres are too close with each other as shown in Figure 8(a). Then the equivalence surfaces are very close to the spheres which makes the field hard to calculate with good accurate. However, it is easy for the GSST-EPA to solve this case because the equivalence surfaces can intersect with each other as shown in Figure 8(b). In Figure 9, the bistatic RCS for the HH polarization of these two methods are compared with MoM. Obviously, the T-EPA result is incorrect while the GSST-EPA has good accuracy.

(a)

(b)

The third example shows a sphere array solved by the new scheme. Each element is a dielectric-coated PEC sphere. The radius of PEC sphere is 0.2 m. The thickness of the dielectric coating is 0.1 m with , . The center distance between two spheres is 0.605 m in and directions. Each sphere is enclosed by an equivalence surface with m³ as shown in Figure 10. The excitation is a polarized plane wave propagating into the negative direction at 300 MHz. In the GSST-EPA, each sphere is solved by volume surface integral equation (VSIE), and MLFMA is also used to accelerate the solution [20]. The result of VSIE_MLFMA with block diagonal preconditioner (BDP) is given as the benchmark as shown in Figure 11. The computational results are shown in Table 1. Obviously, GSST-EPA is more efficient than VSIE_MLFMA for solving this array structure. It is necessary to mention that the memory requirement has not been reduced much although the number of unknowns has been reduced by about 60%. Because the equivalence surfaces can intersect with each other, couplings between the adjacent equivalence surfaces are calculated and stored directly. Only the coupling between the nonadjacent equivalence surfaces is calculated using MLFMA. Therefore the size of finest cube in MLFMA should be chosen twice as large as equivalence surface, which leads to the neighbor groups having a lot of unknowns.

Next, the scattering of a connected pyramid array is investigated. The dimension of a pyramid is shown in Figure 12(a) with , . Each pyramid is enclosed by the same rectangle equivalence surface with cm³. The distance between the centers of two elements is 7 cm which means that the elements connect with each other. The excitation is an polarized plane wave propagating into the negative direction at 1.0 GHz. The array elements are solved by VIE. The comparison of bistatic RCS for HH polarization is shown in Figure 13. It can be seen that these two methods agree with each other very well. The computational results are shown in Table 2. Compared with VIE_MLFMA, the total solution time has been reduced by more than 90% in the GSST-EPA.

(a)

(b)

Finally, to demonstrate the efficiency and capability of the proposed method, the scattering of pyramid array is solved. The parameters of the array are the same as that in the above example, except for the number of elements. The total number of unknowns of the array elements is 105,900, while the number of unknowns on the equivalence surface is 23,400. The error converging to 0.001 only needs 6 iteration steps. The total memory requirement is 3.2 GB and total solution time is 2,501 s. If the VIE_MLFMA is used to solve the array structure, then the memory requirement is more than 16 GB. The result of bistatic RCS for the HH polarization is shown in Figure 14.

8. Conclusions

In this paper, a new scheme named GSST-EPA has been proposed to seek the scattering solution of array structures with very small distance or even connected array elements. This scheme is based on the combination of single-source EPA, tangential field projection, and source reconstruction method. Compared with traditional EPA, this GSST-EPA can reduce the number of unknowns by half, improve the accuracy, and solve the connected structures. The MLFMA is applied further to improve the ability of GSST-EPA to solve large array problems. Numerical results show that the proposed method can solve the connected structures efficiently.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work is supported by National Natural Science Foundation of China (Grant no. 61501267, no. 61631012, no. 61671257, and no. 61425010), Zhejiang Provincial Natural Science Foundation of China (Grant no. LQ16F010003), and Ningbo University Discipline Project (Grant nos. XKL14D2058 and XYL15008) and is partially sponsored by K. C. Wong Magna Fund in Ningbo University.

References

W. C. Chew and L. J. Jiang, “Overview of large-scale computing: The past, the present, and the future,” Proceedings of the IEEE, vol. 101, no. 2, pp. 227–241, 2013.
View at: Publisher Site | Google Scholar
W. B. Lu, T. J. Cui, Z. G. Qian, X. X. Yin, and W. Hong, “Accurate analysis of large-scale periodic structures using an efficient sub-entire-domain basis function method,” IEEE Transactions on Antennas and Propagation, vol. 52, no. 11, pp. 3078–3085, 2004.
View at: Publisher Site | Google Scholar
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
L. Matekovits, V. A. Laza, and G. Vecchi, “Analysis of large complex structures with the synthetic-functions approach,” IEEE Transactions on Antennas and Propagation, vol. 55, no. 9, pp. 2509–2521, 2007.
View at: Publisher Site | Google Scholar
M.-K. Li, W. C. Chew, and L. J. Jiang, “A domain decomposition scheme based on equivalence theorem,” Microwave and Optical Technology Letters, vol. 48, no. 9, pp. 1853–1857, 2006.
View at: Publisher Site | Google Scholar
M. Li and W. C. Chew, “Wave-field interaction with complex structures using equivalence principle algorithm,” IEEE Transactions on Antennas and Propagation, vol. 55, no. 1, pp. 130–138, 2007.
View at: Publisher Site | Google Scholar
M.-K. Li and W. C. Chew, “Multiscale simulation of complex structures using equivalence principle algorithm with high-order field point sampling scheme,” Institute of Electrical and Electronics Engineers. Transactions on Antennas and Propagation, vol. 56, no. 8, part 1, pp. 2389–2397, 2008.
View at: Publisher Site | Google Scholar | MathSciNet
A. M. Van De Water, B. P. De Hon, M. C. Van Beurden, A. G. Tijhuis, and P. De Maagt, “Linear embedding via Greenâs operators: A modeling technique for finite electromagnetic band-gap structures,” Physical Review E - Statistical, Nonlinear, and Soft Matter Physics, vol. 72, no. 5, Article ID 056704, 2005.
View at: Publisher Site | Google Scholar
S. Xiang, G. Xiao, X. Tian, and J. Mao, “Analysis of large-scale phased antenna array with generalized transition matrix,” IEEE Transactions on Antennas and Propagation, vol. 61, no. 11, pp. 5456–5464, 2013.
View at: Publisher Site | Google Scholar
P. Ylä-Oijala and M. Taskinen, “Electromagnetic scattering by large and complex structures with surface equivalence principle algorithm,” Waves in Random and Complex Media, vol. 19, no. 1, pp. 105–125, 2009.
View at: Publisher Site | Google Scholar
H. Shao, J. Hu, H. Guo, F. Ye, W. Lu, and Z. Nie, “Fast simulation of array structures using T-EPA with hierarchical LU decomposition,” IEEE Antennas and Wireless Propagation Letters, vol. 11, pp. 1560–1563, 2012.
View at: Publisher Site | Google Scholar
Z.-H. Ma, L. J. Jiang, and W. C. Chew, “Loop-tree free augmented equivalence principle algorithm for low-frequency problems,” Microwave and Optical Technology Letters, vol. 55, no. 10, pp. 2475–2479, 2013.
View at: Publisher Site | Google Scholar
Y. Shi, J. Wang, and C.-H. Liang, “A time-domain equivalence principle and its marching-on-in-degree solution,” Microwave and Optical Technology Letters, vol. 56, no. 10, pp. 2415–2422, 2014.
View at: Publisher Site | Google Scholar
T. Su, L. Du, and R. Chen, “Electromagnetic scattering for multiple PEC bodies of revolution using equivalence principle algorithm,” Institute of Electrical and Electronics Engineers. Transactions on Antennas and Propagation, vol. 62, no. 5, pp. 2736–2744, 2014.
View at: Publisher Site | Google Scholar | MathSciNet
H. Shao, J. Hu, and W. C. Chew, “Single-source equivalence principle algorithm for the analysis of complex structures,” IEEE Antennas and Wireless Propagation Letters, vol. 13, pp. 1255–1258, 2014.
View at: Publisher Site | Google Scholar
H. Shao, J. Dong, and J. Hu, “Improving the Accuracy of Single-Source Equivalence Principle Algorithm by Using Tangential Field Projection Scheme,” Electromagnetics, vol. 36, no. 8, pp. 534–545, 2016.
View at: Publisher Site | Google Scholar
H. Shao and J. Hu, “Analysis of connected structures using equivalence principle algorithm with source reconstruction method,” Journal of Electromagnetic Waves and Applications, vol. 30, no. 13, pp. 1740–1754, 2016.
View at: Publisher Site | Google Scholar
J. M. Song and W. C. Chew, “Multilevel fast-multipole algorithm for solving combined field integral equations of electromagnetic scattering,” Microwave and Optical Technology Letters, vol. 10, no. 1, pp. 14–19, 1995.
View at: Publisher Site | Google Scholar
S. Wandzura, “Electric current basis functions for curved surfaces,” Electromagnetics, vol. 12, no. 1, pp. 77–91, 1992.
View at: Publisher Site | Google Scholar
C. C. Lu and W. C. Chew, “Electromagnetic scattering from material coated PEC objects: a hybrid volume and surface integral equation approach,” in Proceedings of the IEEE Antennas and Propagation Society International Symposium, IEEE, pp. 2562–2565, Orlando, FL, USA, July, 1999.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2017 Hanru Shao 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.

PDF Download Citation

Download other formats

Order printed copies

Views

645

Downloads

1088

Citations