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

Authors: Hanru Shao; Jianfeng Dong; Jun Hu
Journal: International Journal of Antennas and Propagation
 (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.


---

## Body

## 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 ofN nonconnected perfect electric conductor (PEC) objects in free space is considered to illustrate EPA. Each object is enclosed by an artificial equivalence surface ES l (l=1,…,N). 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 Figure1. 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:(1)Jli=-n^l×Hinc,Mli=n^l×Einc,where n^l is the unit outer normal vector on the ES l. The currents Jli and Mli 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 as(2)η0LJps=-η0Lplη0Kpl·Jli1η0Mli.Here, L and K are the surface integral operators defined as(3)Lr,r′·Xr′=ik∫sXr′+1k2∇∇′·Xr′·G0r,r′ds′,Kr,r′·Xr′=P.V.∫s∇G0r,r′×Xr′ds′,where X is J or M;G0 is Green’s function in free space. η0 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 (4)Ill00Ill·Jls1η0Mls=n^l×Hsca-n^l×Esca=n^l×KlpLlp×n^l·Jps,where I is identity operator. The currents Jls and Mls 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 Sll as(5)Cls=Ill00Ill-1·n^l×KlpLlp×n^l·η0Lpp-1·-η0Lplη0Kpl·Cli=Sll·Cli,where (6)Cls/i=Jls/i1η0Mls/i.Figure 1
A PEC object characterized by equivalence surface with three steps.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 Figure2. The field scattered by Jks and Mks produces new fields on ES l. So these fields can be interpreted as the incident currents from the ES k to ES l and can be expressed by a translation operator Tlk, l≠k, as (7)Clki=Tlk·Cks,where (8)Tlk=Ill00Ill-1·-n^l×Klk-n^l×Llkη0n^l×Llk-η0n^l×Klk.Therefore, the scattered currents on ES l can be written as (9)Cls=Sll·Cli+Tlk·Cks.For the generalized scattering ofN nonadjacent PEC objects enclosed byN equivalence surfaces separately, the EPA equations can be written as(10)I-S11T12⋯-S11T1N-S22T21I⋯-S22T2N⋮⋮⋱⋮-SNNTN1-SNNTN2⋯I·C1sC2s⋮CNs=S11·C1iS22·C2i⋮SNN·CNi.Figure 2
The interaction of two PEC objects using EPA scheme.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.


## 2.1. Scattering Operator

The first procedure can be derived from a PEC object characterized by equivalence surface with three steps as shown in Figure1. 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:(1)Jli=-n^l×Hinc,Mli=n^l×Einc,where n^l is the unit outer normal vector on the ES l. The currents Jli and Mli 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 as(2)η0LJps=-η0Lplη0Kpl·Jli1η0Mli.Here, L and K are the surface integral operators defined as(3)Lr,r′·Xr′=ik∫sXr′+1k2∇∇′·Xr′·G0r,r′ds′,Kr,r′·Xr′=P.V.∫s∇G0r,r′×Xr′ds′,where X is J or M;G0 is Green’s function in free space. η0 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 (4)Ill00Ill·Jls1η0Mls=n^l×Hsca-n^l×Esca=n^l×KlpLlp×n^l·Jps,where I is identity operator. The currents Jls and Mls 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 Sll as(5)Cls=Ill00Ill-1·n^l×KlpLlp×n^l·η0Lpp-1·-η0Lplη0Kpl·Cli=Sll·Cli,where (6)Cls/i=Jls/i1η0Mls/i.Figure 1
A PEC object characterized by equivalence surface with three steps.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 Figure2. The field scattered by Jks and Mks produces new fields on ES l. So these fields can be interpreted as the incident currents from the ES k to ES l and can be expressed by a translation operator Tlk, l≠k, as (7)Clki=Tlk·Cks,where (8)Tlk=Ill00Ill-1·-n^l×Klk-n^l×Llkη0n^l×Llk-η0n^l×Klk.Therefore, the scattered currents on ES l can be written as (9)Cls=Sll·Cli+Tlk·Cks.For the generalized scattering ofN nonadjacent PEC objects enclosed byN equivalence surfaces separately, the EPA equations can be written as(10)I-S11T12⋯-S11T1N-S22T21I⋯-S22T2N⋮⋮⋱⋮-SNNTN1-SNNTN2⋯I·C1sC2s⋮CNs=S11·C1iS22·C2i⋮SNN·CNi.Figure 2
The interaction of two PEC objects using EPA scheme.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 ESl can be expressed by the equivalence electric currents on ES l 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 currentsJli and Mli generate zero field outside the ES l as shown in Figure 3(a):(11)Ei=η0LllJli+12IllrMli-KllMli=0,on  ESl+,Hi=-12IllrJli+KllJli+1η0LllMli=0,on  ESl+,where ES l+ denotes the outside surface of ES l. The integral operators L, K are defined as (3), and Ir is defined as(12)IrX=n^×X.X is equivalence incident or scattering surface current Ji/s or Mi/s. Then the magnetic current can be expressed by electric current as(13)1η0Mli=-12Illr+Kll-1·LllJli=ZEiJli,on  ESl+,(14)1η0Mli=-Lll-1·-12Illr+KllJli=ZHiJli,on  ESl+.Similarly, the relation between equivalence scattering currents Jls and Mls can be established via extinction theorem as shown in Figure 3(b):(15)1η0Mls=12Illr+Kll-1·LllJls=ZEsJls,on  ESl-,(16)1η0Mls=-Lll-1·12Illr+KllJls=ZHsJls,on  ESl-,where ES l- denotes the inside surface of ES l. It can be seen that there are two types of expressions for the magnetic current and electric current, which are named asE-type andH-type. Therefore, the scattering operator of single-source EPA on ES l also has two types and can be written as (17)Jls=Ill-1·n^l×Klp·η0Lpp-1·η0Lpl-η0Kpl·1ZE/HiJli=S~llE/H·Jli.The definition of scattering operator on ES k is similar to the above. The two types of translation operator from ES k to ES l are defined as(18)T~lkE/H=Ill-1·-n^l×Klk-n^l×Llk·1ZE/Hs.Finally, the single-source EPA equations have two types of expressions as follows: (19)I-S~11E/HT~12E/H⋯-S~11E/HT~1NE/H-S~22E/HT~21E/HI⋯-S~22E/HT~2NE/H⋮⋮⋱⋮-S~NNE/HT~N1E/H-S~NNE/HT~N2E/H⋯I·J1sJ2s⋮JNs=S~11E/H·J1iS~22E/H·J2i⋮S~NNE/H·JNi.Once the electric currents have been solved, the magnetic currents can be obtained using (15) or (16).Figure 3
(a) Outside-in equivalence process. (b) Inside-out equivalence process.
(a)
(b)


## 4. Single-Source Tangential Equivalence Principle Algorithm

In the single-source EPA, the equivalence surface currentsJ are the rotated tangential projections of the fields onto the surface as follows:(20)IJs/i=±n^×Hs/i,where Hs/i 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 l by using Huygens equivalence principle as follows:(21)η0LllJli-12IllrMli-KllMli=Ei,on  ESl-,(22)12IllrJli+KllJli+1η0LllMli=Hi,on  ESl-,(23)η0LllJls+12IllrMls-KllMls=Es,on  ESl+,(24)-12IllrJls+KllJls+1η0LllMls=Hs,ESl+.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 Jlli and scattering current Jlls separately, which can be denoted asTEE,TEH,THE, andTHH. Numerical results have shown that ZHi and ZHs are more accurate than ZEi and ZEs when using RWG function as basis and testing functions [15]. It also has been proved that theTHH type is the same as single-source EPA [16]. Therefore, only theTEH is used in this paper.The single-source tangential EPA (T-EPA) based on theTEH type can be derived by submitting (14) into (21) and (16) into (23) as follows:(25)η0LllJli-12η0Illr·ZHiJli-η0Kll·ZHiJli=ZTEHiJli=Ei,on  ESl-,(26)η0LllJls+12η0Illr·ZHsJls-η0Kll·ZHsJls=ZTEHsJls=Es,on  ESl+.Then the scattering operator of single-source T-EPA on ES l can be written as (27)Jls=ZTEHs-1·η0Llp·η0Lpp-1·η0Lpl-η0Kpl·1ZHi·Jli=S~llTEH·Jliand the translation operator should be changed as(28)T~lkTEH=ZTEHi-1·η0Llk-η0Klk·1ZHs.It can be seen that both of scattering operator and translation operator have the combination ofL and ±1/2Ir+K 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 Figure4. The left part PEC 1 is enclosed by ES 1, and the right part PEC 2 is enclosed by ES 2. J1p is the current on PEC 1; J2p is the current on PEC 2. The MoM equations can be generated as(29)η0Lp1p1J1p+η0Lp1p2J2p=-E1i,(30)η0Lp2p2J2p+η0Lp2p1J1p=-E2i.The equivalence scattering currents J1s on ES 1 can be obtained through inside-out propagation:(31)J1s=ZTEH1s-1·η0Ls1p1·J1p,where the operator ZTEH1s has been defined in (26). Then substituting (29) into (31) yields(32)J1s=ZTEH1s-1·η0Ls1p1·η0Lp1p1-1·-η0Lp1p2·J2p-E1i.E1i is the incident plane wave on PEC 1 which is already known, while J2p needs to be solved. Suppose that the equivalence scattering currents J2s and M2s on ES 2 are already known, then the J2p can be got by using source reconstruction method [17] as follows:(33)J2p=η0Ls2p2+·ZTEH2s·J2s,where (·)+ denotes the pseudoinverse. In this paper, the pseudoinverse is computed by truncated SVD method. E1i can be replaced by equivalence incident currents as(34)E1i=η0Lp1s1-η0Kp1s1·1ZH1i·J1i.Hence, substituting (33) and (34) into (32), we have(35)J1s=T~s1s2SRM·J2s+S~s1s1TEH·J1i,where (36)T~s1s2SRM=ZTEH1s-1·η0Ls1p1·η0Lp1p1-1·-η0Lp1p2·η0Ls2p2+·ZTEH2sand the definition of S~s1s1TEH is the same as (17). Similarly, the equation for the equivalence scattering currents on ES 2 can be written as(37)J2s=T~s2s1SRM·J1s+S~s2s2TEH·J2i.Figure 4
A PEC strip solved using source reconstruction method.Finally, the new EPA equations for connected structures can be obtained by combining (35) and (37):(38)I-T~s1s2SRM-T~s2s1SRMI·J1sJ2s=S~s1s1TEH·J1iS~s2s2TEH·J2i.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 ofN connected objects, each object is enclosed by an equivalence surface as shown in Figure 5. Then the GSST-EPA equation can be written as(39)I-Z12⋯-Z1N-Z21I⋯-Z2N⋮⋮⋱⋮-ZN1-ZN2⋯I·J1sJ2s⋮JNs=S~11TEH·J1iS~22TEH·J2i⋮S~NNTEH·JNi,where(40)Zij=S~iiTEHT~ijTEH,ifiandjare  non-adjacent,T~ijSRM,ifiandjare  adjacent.It 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 S~iiTEH is identical when the array element and equivalence surface are the same as each other. For the periodic array structures, T~ijTEH and T~ijSRM 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.Figure 5
Each object is enclosed by an equivalence surface.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 matrixT~ijSRM is calculated directly and stored. Secondly, the nonadjacent case is the far interaction which is solved by translation operator. So the translation matrix T~ijTEH 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 inx-axis. Each sphere is enclosed by equivalence surface separately, which is a rectangle surface with 0.8×0.8×0.8 m3 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.Figure 6
Two spheres are enclosed by two equivalence surfaces.Figure 7
The bistatic RCS of two PEC spheres with center distance 1.0 m.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 Figure8. 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.Figure 8
Two PEC spheres enclosed by two ES, (a) T-EPA scheme, and (b) GSST-EPA scheme.
(a)
(b)Figure 9
The bistatic RCS of two PEC spheres with center distance 0.605 m.The third example shows a5×5 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 εr=2.2, μr=1.0. The center distance between two spheres is 0.605 m in x^ and y^ directions. Each sphere is enclosed by an equivalence surface with 0.8×0.8×0.8 m3 as shown in Figure 10. The excitation is a polarized plane wave propagating into the negative z^ 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.Table 1
The computational results of two methods.
Method
Unknowns
Memory (MB)
Iteration number (0.001)
Total time (s)
VSIE_MLFMA
47,314
2,220.6
251
4,442.4
GSST-EPA
13,050
1,888.3
22
2,421.3Figure 10
The model of a5×5 sphere array.Figure 11
The bistatic RCS of a5×5 sphere array; HH polarization at 300 MHz.Next, the scattering of a connected pyramid array is investigated. The dimension of a pyramid is shown in Figure12(a) with εr=2.0-j0.04, μr=2.0-j0.04. Each pyramid is enclosed by the same rectangle equivalence surface with 9×9×17 cm3. The distance between the centers of two elements is 7 cm which means that the elements connect with each other. The excitation is an x^ polarized plane wave propagating into the negative z^ 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.Table 2
The comparison of computational results.
Method
Unknowns
Memory (MB)
Iteration number (0.001)
Total time (s)
VIE_MLFMA
27,121
3,285.7
11
8,112.7
GSST-EPA
5,850
422.5
5
575.9Figure 12
The model of pyramid array: (a) element; (b)5×5 array and equivalence surface.
(a)
(b)Figure 13
The bistatic RCS of5×5 pyramid array; HH polarization at 1.0 GHz.Finally, to demonstrate the efficiency and capability of the proposed method, the scattering of10×10 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.Figure 14
The bistatic RCS of10×10 pyramid array; HH polarization at 1.0 GHz.


## 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.



---

*Source: 4903906-2017-08-06.xml*