# Solving Scattering from Conducting Body Coated by Thin-Layer Material by Hybrid Shell Vector Element with Boundary Integral Method

Authors: Lin Lei; Jun Hu; Hao-Quan Hu
Journal: International Journal of Antennas and Propagation
 (2012)

---

## Abstract

The finite element boundary integral (FEM-BI) method is widely used in the scattering and radiating problems. But for the conducting body coated by thin-layer material, plenty of fine meshes are required to discretize the geometry in the traditional FEM. It requires very expensive storage and CPU time. In this paper, the hybrid shell vector element with the boundary integral method is used to expedite the solution of thin coating problems. The shell vector elements are used to discretize thin-layer material instead of traditional tetrahedral elements. Consequently, the volume integral can be simplified into surface integral. This method reduces the number of unknowns greatly and is also extended into the complicated case of multi-thin-layer coating materials. Several numerical results are presented to prove the accuracy and efficiency of this present method.


---

## Body

## 1. Introduction

The electromagnetic analysis of the conducting body coated by thin-layer material has received much attention because these composite dielectric and conductor structures are used in many applications, such as the evaluation of echo from stealth aircraft and microwave integrated circuits design. Plenty of research work has been done on the efficient solution of scattering from conducting body coated by thin-layer material.For the analysis of conductor structures coated by thin-layer material, there are many methods available. Integral equation methods based on thin dielectric sheet approximation (TDS) [1–7] avoid the volumetric discretization of material region, the computational region is only limited as the surface of conductor. A multilevel-TDS extension is also proposed for solution of conducting body coated by multi-thin-layer materials [8–10]. For lossy thin-layer coating, the impedance boundary condition (IBC) is also developed to simplify the electromagnetic analysis by building the relation between the equivalent surface magnetic current and equivalent surface electric current on the surface of conductor [11]. The application of the IBC into the conducting bodies of revolution coated with thin magnetic materials is also developed successfully [12]. A rigorous moment method solution of composite bodies with thin features can achieve accurate results for thin coating problems if the near singularity for near-self term integral can be tackled carefully, as shown in [13, 14].Besides the integral equation methods, the finite element method (FEM) is also widely used for analysis of composite conducting body and dielectric because of its powerful ability of modeling complex inhomogeneous materials [15]. As well known, the absorption boundary condition (ABC) is required in FEM in order to truncate the computational region. In order to combine together the advantage of FEM and integral equation, the hybrid FEM with the boundary integral method (FEM-BI) is proposed [15–23]. No ABC is needed in the FEM-BI, and the number of unknowns is reduced greatly.Though FEM-BI has a good computational property for composite conducting body and dielectric, it is deficient for the analysis of conductor structures coated by thin-layer material. This is because plenty of fine meshes will be required for modeling thin layers if using traditional elements like tetrahedral elements. To further reduce the total number of unknowns, prism elements are used in FEM when solving thin-coating structures [23]. In order to further enhance the computational efficiency, shell vector elements are developed in FEM, instead of prism elements [24, 25]. The volume integral can be simplified into surface integral by using shell vector elements.In this paper, we combine the shell vector element method with the boundary integral method in order to achieve faster solution of scattering from conducting body coated by thin-layer material. The present method does not only reduce the number of unknowns but is also extended into more complicated cases of multi-thin-layer coating. This paper is organized as follows. The basic principle of the traditional FEM-BI method is reviewed in section II. Next, we show the characteristics of the shell elements, and the hybrid shell vector element method with the boundary integral method is implemented. Finally, several examples are shown to demonstrate the efficiency and accuracy of this present method. The conclusions are also given.


## 2. Traditional FEM-BI Method

Considering the electromagnetic scattering problem shown in Figure1, the structure is illuminated by an incident plane wave Einc or Hinc, the n̂ is the unit normal vector on surfaceS and S is the surface enclosing volume V.Figure 1
The 3D electromagnetic scattering problem by the FEM-BI.TheE field inside volume V satisfies the following equation:(1)∇×(1μr∇×E)-k02εrE=0,
where μr and εr denote the relative permeability and permittivity of the media, respectively, and k0 is the wave number in free space.The boundary conditions on surfaceS are written as(2)n̂×E|S-=n̂×E|S+,n̂×(1μr∇×E)|S-=-jk0n̂×H0|S+,
where H0=η0H and η0 is the wave impedance in free space.The FEM matrix can be constructed based on the variational method; the functionalF(E) is expressed as follows:(3)F(E)=12∫∫∫V[1μr(∇×E)⋅(∇×E)-k02εrE⋅E]dV+jk0η0∫∫Sn̂⋅(E×Hs)dS.The tetrahedral element is most widely applied in traditional FEM and the electric fieldE can be expanded by the tetrahedral element as(4)E=∑j=16EjeNje,
where Eje and Nje are the unknown coefficient and the basis function of the jth edge, respectively.For the surface integral term in (3), the magnetic field on the surface can be expanded by three edges of planar triangle on surface S as(5)Hs=∑j=13HjsNjs,
where Hjs and Njs are the unknown coefficient and the basis function of the jth edge, respectively.Finally, the following FEM matrix equation is yielded from (3):(6)[K]{E}+[B]{H}={0}.In the FEM-BI method, surface integral equation method is applied on the surfaceS. The electric field integral equation (EFIE) and the magnetic field integral equation (MFIE) are used as follows:(7)Einc=L(J̅)-K(M),H̅inc=K(J̅)+L(M).The equivalent surface electric and magnetic currents are defined as(8)J=n̂×H,M=E×n̂,J̅=n̂×H̅.
here H̅inc=η0Hinc. The integral operators L and K are expressed, respectively, as(9)L(X)=jk0∫s′[X(r′)+1k02∇∇′⋅X(r′)]G(r,r′)ds′,(10)K(X)=TY(r)+∫s′X(r′)×∇G(r,r′)ds′.In (10), Y(r)=X(r)×n̂, T=1-Ω/4π. For the smooth surface, Ω=2π. G=e-jk0|r-r′|/4π|r-r′| is Green’s function in free space. The integral terms with bars in (10) denote the principal value integrals. The singularity is excluded.By choosingSis=n̂×Ni as the weighting function, the TE formulation can be derived as(11)[PTE]{ES}+[QTE]{HS}={bTE},
where(12)PijTE=-∫SSis⋅K(Nj×n̂)ds,QijTE=∫SSis⋅L(η0n̂×Nj)ds,biTE=∫SSis⋅Eids.Similarly, the NE formulation can be derived by choosingn̂×Sis as the weighting function as(13)[PNE]{ES}+[QNE]{HS}={bNE}.The TH formulation can be derived by choosingSis as the weighting function as:(14)[PTH]{ES}+[QTH]{HS}={bTH}.The NH formulation can be derived by choosingn̂×Sis as the weighting function as:(15)[PNH]{ES}+[QNH]{HS}={bNH}.
To avoid the interior resonance that occurred in some cases, the combined field integral equation (CFIE), which combines the EFIE and MFIE, can be used,(16)[P]{ES}+[Q]{HS}={b},
where(17)[P]=α[PTE]+β[PNE]+γ[PNH],[Q]=α[QTE]+β[QNE]+γ[QNH],{b}=α{bTE}+β{bNE}+γ{bNH}.The parametersα, β, and γ satisfy: α+β+γ=1.By combining (6) and (16), we get the matrix equation of the FEM-BI method as follows:(18)[KIIKIS0KSIKSSB0PQ]{EIESH̅S}={00b},
where EI, ES, and H̅S are the coefficient of the electric field inside volume V, on surface S, and the magnetic field on surface S, respectively.


## 3. Hybrid Shell Element with Boundary Integral (Hybrid Shell Element-BI)

The shell element is the degenerated prism element. As shown in Figure2, there are six edge vectors along the corresponding edges in the upper triangle and bottom triangle and three normal vectors. A linear function β(ς) was used to describe the variation of the field along the normal direction, and ∇β=-n̂/d. The βNj(j=1,2,3) and β′Nj′(j=1,2,3) are the edge basis functions in the upper and bottom triangle, respectively. The Lj, j=1, 2, 3 is the normal basis function at node j.The structure of the prism element and the shell element.
(a)
The prism element(b)
The shell elementBy the shell element, the electric field is expanded as follows:(19)Ee=∑j=13(EjeβNj+E′jeβ′Nj′)+∑j=13EnjeLjen̂,
where β′=1-β, Nj=(Lj1e∇Lj2e-Lj2e∇Lj1e)lje, Nj′=(Lj1e∇Lj2e-Lj2e∇Lj1e)lje,  Eje is the expansion coefficient of the jth edge vector in the upper triangle, Ej′e is the expansion coefficient of the jth edge vector in the bottom triangle, Enje is the expansion coefficient of the normal vector at node j.Based on (3) and (19), we attain(20)∫∫∫v{∑j=13EnjeLjen̂∇×[βNie+β′Nie+Lien̂]⋅∇×[∑j=13(EjeβNje+Ej′eβ′Nje)+∑j=13EnjeLjen̂]μr-k02εr[βNie+β′Nie+Lien̂]⋅[∑j=13(EjeβNje+E′jeβ′Nje)+∑j=13EnjeLjen̂]}dV=jk0η0∑j=13Hjs∫∫sNis⋅(n̂×Njs)ds.For the shell element model, the volume integral can be simplified as:(21)∫∫∫V(∇×βNie)⋅(∇×βNje)dv=∫-1/21/2ββdς∫∫s(∇×Nie)⋅(∇×Nje)ds-∫-1/21/2βddς∫∫s(∇×Nie)⋅Sjeds-∫-1/21/2βddς∫∫sSie(∇×Nje)ds+∫-1/21/21d2dς∫∫sSie⋅Sjeds,(22)∫∫∫V(∇×Lien̂)⋅(∇×β′Nje)dv=∫-1/21/2β′dς∫∫s(∇Lie×n̂)⋅(∇×Nje)ds+∫-1/21/21ddς∫∫s(∇Lie×n̂)⋅Sjeds,
where Sie=n̂×Nie, Sje=n̂×Nje and d is the thickness of the thin-layer material.Obviously, the integration along the normal direction can be calculated analytically. So the volume integral can be simplified into the surface integral.For the conducting object coated by thin-layer material, the electric field in the bottom surface of shell element must be zero, only the integral in the upper surface is needed.The hybrid shell element-BI can also be extended into the conducting body coated by multi-thin-layered material. For example, for the conducting body coated by three-thin-layer material, the matrix equation is expressed as:(23)[Knn1Knu10Kun1Kuu1+bb2Kbn2Kbu2Knb2Knn2Knu2Kub2Kun2Kuu2+bb3Kbn3Kbu3Knb3Knn3Knu3Kub3Kun3Kuu3B0PQ]{En1Eu1En2Eu2En3Eu3Hs}={000000b},
where En1 and Eu1 are the expansion coefficients of the electric field in the normal and upper surface of the first layer, respectively; En2 and Eu2 are the expansion coefficients of the electric field in the normal and upper surface of the second layer, respectively, En3 and Eu3 are the expansion coefficients of the electric field in the normal and upper surface of the third layer, respectively.


## 4. Numerical Results

To demonstrate the accuracy and efficiency of the present method, some typical numerical results are shown here.The first example is a conducting sphere coated by magnetic material, shown in Figure3. The frequency is 300 MHz. The radius of the conducting sphere is 0.3367 m, and the thickness of the coating layer is 0.001 m. The relative permeability is μr=4 and the relative permittivity is εr=1. In the hybrid shell element with BI, the mesh density of 15 elements per wavelength is enough for achieving reasonable accuracy, but in traditional FEM, at least 32 elements per wavelength are required in order to mesh this very thin layer successfully. Table 1 demonstrates the comparisons between the two methods. As shown in Figure 4, the results by the hybrid shell element with BI agree with the ones of MIE method very well. And, the number of the unknowns by the traditional FEM-BI and the present method is 15167 and 2249, respectively. Obviously, the advantage of the present method over the traditional FEM-BI in reducing the number of unknowns is very remarkable.Table 1
The comparison of the FEM-BI and the hybrid shell element-BI. Bistatic RCS of the conducting sphere coated by magnetic material is calculated.
MethodMesh density (wavelength)Node numberElement typeElement numberUnknownsFEM-BI32 elements3033Tetrahedral elements913215167Hybrid shell element-BI15 elements323Shell elements6422249Figure 3
Conducting sphere coated by magnetic material.Figure 4
The bistatic RCS of the conducting sphere coated by magnetic material by the hybrid shell element-BI.To further investigate the accuracy of the hybrid shell element-BI, Figure6 shows the bistatic RCS of the conducting sphere coated by thin dielectric material, the geometry is shown in Figure 5, k0a=6.28, the thickness of the thin layer is a/100, and the relative permittivity is εr=4, where k0 is the wave number in free space and a is the radius of the conducting sphere. The mesh density used here is only 6 elements per wavelength. As shown in Figure 6, the results by the hybrid shell element -BI agree well with the ones of the MIE method. In this case, the number of unknowns is 10530.Figure 5
Coated metal sphere.Figure 6
The bistatic RCS of the conducting sphere coated by thin dielectric material by the hybrid shell element with BI.Figure8 shows the bistatic RCS of the conducting cube coated by thin material in horizontal polarization. The geometry is shown in Figure 7. The length of the cube is 0.1m, and the thickness of the coated layer is 0.001m. The wavelength is 0.1m.The incident angle of the plane wave is θinc=45°,φinc=0°.Figure 7
Conducting cube coated by thin material.Figure 8
The bistatic RCS of the coated metallic cube.From Figure8, it is shown clearly that the results by hybrid shell element-BI agree well with the ones by the commercial software FEKO. The number of unknowns in the hybrid shell element-BI is 5630, but the one in the traditional FEM-BI is 7760.The solving CPU time by the present method is only one-third of the one by the traditional FEM-BI. The details are summarized in Table2. Obviously, compared with the FEM-BI, the hybrid shell element with BI can save memory and CPU time greatly.Table 2
The comparison of the FEM-BI and the hybrid shell element-BI. Bistatic RCS of the conducting cube coated by thin material (mesh density with 10elements/λ).
MethodNode numberMemoryElement numberUnknownsFEM-BI1474480 Mb44047760Hybrid shell element-BI806250 Mb16085630Figure9 shows a conducting sphere with radius of 0.5λ0 coated by three-thin-layer material, the thickness of each layer is 0.01λ0, μr1=1.0, εr1=2.0, μr2=2.0-j2.0, εr2=2.0-j2.0, μr3=2.0, and εr3=2.0. The number of unknowns is 9591 and the total memory required is 736 Mb, but the number of unknowns in the traditional FEM-BI is 18775 and the total memory required is 2.9 Gb. The bistatic RCS computed by the hybrid shell element-BI is shown in Figure 10, the one by the MIE method is also shown for comparison. The lines marked with dot and square denote the curves of results in HH polarization and VV polarization, respectively. Obviously, they agree well with each other.Figure 9
Conducting sphere coated by three-thin-layer materials.Figure 10
The bistatic RCS of the conducting sphere coated by three-thin-layer material.


## 5. Conclusions

In this paper, the hybrid shell element with BI is developed to achieve efficient solution of scattering from the conducting body coated by thin-layer electric or magnetic material. Compared with the traditional FEM-BI, the present method reduces the number of unknowns greatly. Numerical results show that the hybrid shell element with BI is accurate and effective in solving the scattering problems of the conducting body coated with thin layer, also valid for the case of multi-thin-layer material coating.



---

*Source: 854647-2012-04-10.xml*