# Fast Direct Solution of Electromagnetic Scattering from Left-Handed Materials Coated Target over Wide Angle

Authors: Guo-hua Wang; Ying-bao Geng
Journal: Journal of Electrical and Computer Engineering
 (2016)

---

## Abstract

When solving the electromagnetic scattering problems over wide angle, the traditional method of moments (MoM) needs to repeat the solving process of dense systems of linear equations using the iteration method at each incident angle, which proved to be quite inefficient. To circumvent this problem, a fast numerical method based on block LDLT factorization accelerated by adaptive cross approximation (ACA) algorithm is presented to analyze the electromagnetic scattering of left-handed materials (LHM) coated target. The ACA algorithm is applied to impedance matrix filling and all steps of block LDLT factorization process, which can accelerate the computation process and reduce the memory consumption. The numerical results proved that the proposed method is efficient in calculating monostatic RCS of LHM coated target with many required sampling angles. Compared with the traditional MoM, computation time and memory consumption are reduced effectively.


---

## Body

## 1. Introduction

Left-handed materials (LHM) [1, 2] have electromagnetic properties which cannot be generally found in nature; the electric field direction, magnetic field direction, and the propagation vector conform to the laws of the left hand when electromagnetic wave travels in it. LHM are also called double negative materials (DNG) [3] as the dielectric permittivity and permeability at a certain frequency range can be negative together. The application of LHM has been wide in antenna [4], microwave devices, and other fields [5] due to their unusual physical properties, such as negative refraction, inverse Doppler shift, and reversed Cerenkov radiation. For the purpose of stealth, many conductor target surfaces are also coated by LHM; this necessitates a fast and efficient solution to analyze the electromagnetic scattering characteristics of LHM coated targets [6–9].Researchers have proposed various methods to analyze the scattering problem of LHM coated targets. In literature [7], the radar cross section of metallic spheres covered by DNG based on Mie series solution is studied, but this method is valid only in certain circumstances when the geometric structure of the target is simple. In literature [8], the physical optics (PO) method in combination with the impedance boundary condition (IBC) is used to analyze the near-field electromagnetic scattering characteristics of LHM coated targets. In literature [9], the wide band scattering of perfectly conducting target coated with DNG is computed using finite difference time domain (FDTD) method. Apart from the above methods, this paper mainly discusses the fast method to compute the RCS of LHM coated target over wide angle based on MoM [10] framework and hybrid PEC-dielectric formulation. However, when solving the electromagnetic scattering problems with multiple incident angles, traditional MoM needs to compute the current coefficients as the change of incident angle. Usually iteration method is used to solve the linear system equations which proved to be a quite effective method for bistatic RCS computing, but, for monostatic RCS with many sampling angles, the iteration operations must be repeated when the incident angle changes each time; this is inefficient when the number of sampling angles is large; also there are often convergence issues when dealing with complex structure targets. To circumvent this problem, a new block ACA-LDLT method is proposed in this paper to speed up the filling process of impedance matrix and solve the linear system equations more effectively, which is suitable for computing monostatic RCS with many required sampling angles of complex structure targets.


## 2. Hybrid PEC-Dielectric Formulation for LHM Coated Target

Recently, a method [11] based on thin dielectric sheet (TDS) approximation with explicit perfect electric conductor (PEC) boundary conditions at the interfaces of PEC and dielectrics was provided, which can get accurate numerical solution of EM scattering but is only suitable for electrical coating materials. In literature [12], a generalized thin coating equivalent model was proposed to be suitable for both electrical and magnetic coating materials, but there is no further discussion when the coating medium is LHM. In this section, we mainly discuss the building up of the hybrid PEC-dielectric formulation for LHM coated target.For LHM coated PEC target under the irradiation of ideal plane wave shown in Figure1, S stands for the surface of PEC, V is the volume of coating LHM, and τ is the coating thickness. The total tangential component of electric field is zero following the PEC boundary condition of the electric field; we can obtain(1)Eincr+Escartan=0r∈S;here, the subscript “tan” stands for the tangential component, r stands for the field point, Einc stands for the incident field, and Esca is the scattering field which could be described as(2)Escar=Epecscar+Ediescar.The scattering field is produced by conductors and LHM coating medium together, and the scattering field produced by the electric current on S could be described as(3)Epecscar=jωμ0∫SJSr′+1k02∇∇′·JSr′gr,r′dS′;here, g(r,r′)=e-jkr-r′/4πr-r′ stands for scalar Green’s function in free space, JS is surface current density on S, and μ0 and k0 stand for the permeability and wave number in free space, respectively. The scattering field produced by LHM coating medium in V could be described as(4)Ediescar=jωμ0∫Vχr′gr,r′jωDr′dV′-1k02∫Sn+χ∇gr,r′jωn^′·Dr′dS′+1k02∫Sn-χ∇gr,r′jωn^′·Dr′dS′-∫V∇gr,r′×jωξr′Br′dV′,where χ(r′)=1/εr(r′)-1 and ξ(r′)=1/μr(r′)-1 are the dielectric contrast ratios; n^′ is the unit vector normal to the upper surface. When the coating thickness τ is very small compared to the wavelength, we can use an approximate method named as TDS approximation: (5)Dr′≈n^′∇′·JSr′jωBr′≈-μn^′×JSr′r′∈V.Then, the volume integral can be further approximated to surface integral through conversion dV≈τdS. Therefore, the scattering field produced by dielectric can be described as(6)Ediescar=jωμ01-μr2JSτ+∫Sn^′τχg∇′·JSdS′+∇k02∫Sχg-gτ∇′·JSdS′+1-μr∫S-S0τ∇g×n^′×JSdS′.Here, the last term in (6) is the principal value integral. Finally, we substitute (3) and (6) in (2); in combination with (1), we can obtain the hybrid PEC-dielectric formulation which can be expressed by the only unknown JS:(7)Eincr=-jωμ01-μr2JSτ+∫SgJS+n^′τχ∇′·JSdS′+∇k02∫Sg+χg-χgτ∇′·JSdS′+1-μr∫Sτ∇g×n^′×JSdS′r∈S.Here, εr and μr are permittivity and permeability of coating LHM given by Drude model [13, 14]; gτ is scalar Green’s function defined by(8)gτ=e-jkr-r′+τn^′4πr-r′+τn^′.Figure 1
A PEC target coated by LHM.


## 3. Solutions

We will refer to (8) as a hybrid PEC-dielectric formulation, which can be solved using the MoM or MLFMA with iteration method. But, for monostatic RCS with many required sampling angles, iteration method must be repeated with many right-hand sides (RHS), so this part becomes expensive. The block LDLT factorization method in combination with adaptive cross approximation (ACA) is introduced to this paper which can reduce time and memory storage and is proven to be more efficient in computing the monostatic RCS than the method in [12].


### 3.1. Block LDLT Factorization Method

By discretizing the surface currents using RWG vector basis functions [15, 16], the integral equation (8) can be transformed into a dense, complex linear equation:(9)Z·J=V,where Z is the impedance matrix of dimension N×N, J is induced current density, and V is excitation voltage matrix; N is the unknowns number of closed scattering mesh. When using MoM to solve the equation, the storage complexity and matrix filling time is O(N2) and the direct LU factorization complexity is O(N3); both grow rapidly for electrically large target and bring a heavy burden to computer. To overcome the problem, the unknowns in this work have been grouped into multiple local regions, and the rank of the submatrices decreases with the increase of distance of two blocks; then the compressed algorithm can be used to reduce the storage and operations count. When unknowns are grouped into local spatial regions, (9) can be converted into the following form:(10)Z11Z12⋯Z1MZ21Z22⋯Z2M⋮⋮⋱⋮ZM1ZM2⋯ZMMJ1J2⋮JM=V1V2⋮VM,where Zij (i,j=1,2,3,…,M) represent the submatrix of impedance matrix. By using RWG vector basis functions as the test functions, the impedance matrix in (10) is a complex coefficient symmetric matrix, which can be expressed as block LDLT form:(11)Z11Z12⋯Z1MZ21Z22⋯Z2M⋮⋮⋱⋮ZM1ZM2⋯ZMM=10⋯0L211⋯0⋮⋮⋱⋮LM1LM2⋯1D110⋯00D22⋯0⋮⋮⋱⋮00⋯DMM1L21T⋯LM1T01⋯LM2T⋮⋮⋱⋮00⋯1,where Lij are the lower triangular entities and Dii are the diagonal entities; then making the substitutions Uii=Dii and Uij=DiiLjiT, we can obtain the standard block LU form:(12)Z11Z12⋯Z1MZ21Z22⋯Z2M⋮⋮⋱⋮ZM1ZM2⋯ZMM=10⋯0L211⋯0⋮⋮⋱⋮LM1LM2⋯1U11U12⋯U1M0U22⋯U2M⋮⋮⋱⋮00⋯UMM.The iterative formula ofUij is(13)Uij=Zij-∑p=1i-1UpiTDpp-1Upj.The surface currents can be obtained through the forward iterative step(14)Xi=Vi-∑p=1i-1UpiTDpp-1Xpand back iterative step(15)Ji=Dii-1Xi-∑p=i+1nUipJp.


### 3.2. ACA Acceleration

The impedance matrix in (10) is composed of multiple submatrices, and the off-diagonal blocks are of low rank and can be compressed. In literature [17, 18], Shaeffer points out that the upper triangular matrix U and multiple plane wave excitation RHS voltage matrix V also have low rank characteristics and can be compressed too, so the compression operation can be used for all steps of the solutions including impedance filling, LU factorization, and LU solving. In numerous matrix compression methods, the ACA algorithm in literature [19] is widely used [20–22]. The ACA algorithm operates by scanning a row followed by a column of the matrix at each iteration and progressively builds up a low rank estimate of the matrix based on the rows and columns that have been scanned. The algorithm operates as follows:(1)
Initialize the approximate matrix asA~m×n≈0m×n and the iteration count as j=1; arbitrarily choose a row of the matrix.(2)
Scan thejth row of the matrix.(3)
Find the error of the previous field approximation at thejth row R=Am×n(jthrow)-A~m×n(jthrow) and choose the jth column to be the one containing the maximum element of R.(4)
Assignvj1×n=R/R (jth column).(5)
Scan thejth column and find the error of the previous estimate at this column: C=Am×n(jthcolumn)-A~m×n(jthcolumn). Choose the (j+1)th row to be the one containing the maximum element of C (ensuring that the jth row is not chosen again).(6)
Assignujm×1=C.(7)
Update the field approximationA~m×n=A~m×n+ujm×1vj1×n.(8)
Ifujm×1vj1×n<εA~m×n, stop scanning; else increment j and repeat steps (2) to (8).A low rank matrix can be well approximated by the product of two full rank matrices using ACA method:(16)Am×n≈A~m×n=AUm×rAVr×n,where r is the effective rank of matrix and the memory requirement decreased from m×n entries to r(m+n) entries. We can get the block decomposition expression by using the ACA to the iterative process of Uij: (17)Uij=UUUVi,j=ZUZVij-∑p=1i-1UVTUUTpiDpp-1UUUVpj.By using the ACA to the LU solving solution, the induced current densityJi on S can be obtained by the forward iterative steps(18)Xi=Vi-∑p=1i-1UVTUUTpiDpp-1XUXVpand back iterative steps(19)Ji=Dii-1Xi-∑p=i+1nUUUVipJUJVp.


## 3.1. Block LDLT Factorization Method

By discretizing the surface currents using RWG vector basis functions [15, 16], the integral equation (8) can be transformed into a dense, complex linear equation:(9)Z·J=V,where Z is the impedance matrix of dimension N×N, J is induced current density, and V is excitation voltage matrix; N is the unknowns number of closed scattering mesh. When using MoM to solve the equation, the storage complexity and matrix filling time is O(N2) and the direct LU factorization complexity is O(N3); both grow rapidly for electrically large target and bring a heavy burden to computer. To overcome the problem, the unknowns in this work have been grouped into multiple local regions, and the rank of the submatrices decreases with the increase of distance of two blocks; then the compressed algorithm can be used to reduce the storage and operations count. When unknowns are grouped into local spatial regions, (9) can be converted into the following form:(10)Z11Z12⋯Z1MZ21Z22⋯Z2M⋮⋮⋱⋮ZM1ZM2⋯ZMMJ1J2⋮JM=V1V2⋮VM,where Zij (i,j=1,2,3,…,M) represent the submatrix of impedance matrix. By using RWG vector basis functions as the test functions, the impedance matrix in (10) is a complex coefficient symmetric matrix, which can be expressed as block LDLT form:(11)Z11Z12⋯Z1MZ21Z22⋯Z2M⋮⋮⋱⋮ZM1ZM2⋯ZMM=10⋯0L211⋯0⋮⋮⋱⋮LM1LM2⋯1D110⋯00D22⋯0⋮⋮⋱⋮00⋯DMM1L21T⋯LM1T01⋯LM2T⋮⋮⋱⋮00⋯1,where Lij are the lower triangular entities and Dii are the diagonal entities; then making the substitutions Uii=Dii and Uij=DiiLjiT, we can obtain the standard block LU form:(12)Z11Z12⋯Z1MZ21Z22⋯Z2M⋮⋮⋱⋮ZM1ZM2⋯ZMM=10⋯0L211⋯0⋮⋮⋱⋮LM1LM2⋯1U11U12⋯U1M0U22⋯U2M⋮⋮⋱⋮00⋯UMM.The iterative formula ofUij is(13)Uij=Zij-∑p=1i-1UpiTDpp-1Upj.The surface currents can be obtained through the forward iterative step(14)Xi=Vi-∑p=1i-1UpiTDpp-1Xpand back iterative step(15)Ji=Dii-1Xi-∑p=i+1nUipJp.


## 3.2. ACA Acceleration

The impedance matrix in (10) is composed of multiple submatrices, and the off-diagonal blocks are of low rank and can be compressed. In literature [17, 18], Shaeffer points out that the upper triangular matrix U and multiple plane wave excitation RHS voltage matrix V also have low rank characteristics and can be compressed too, so the compression operation can be used for all steps of the solutions including impedance filling, LU factorization, and LU solving. In numerous matrix compression methods, the ACA algorithm in literature [19] is widely used [20–22]. The ACA algorithm operates by scanning a row followed by a column of the matrix at each iteration and progressively builds up a low rank estimate of the matrix based on the rows and columns that have been scanned. The algorithm operates as follows:(1)
Initialize the approximate matrix asA~m×n≈0m×n and the iteration count as j=1; arbitrarily choose a row of the matrix.(2)
Scan thejth row of the matrix.(3)
Find the error of the previous field approximation at thejth row R=Am×n(jthrow)-A~m×n(jthrow) and choose the jth column to be the one containing the maximum element of R.(4)
Assignvj1×n=R/R (jth column).(5)
Scan thejth column and find the error of the previous estimate at this column: C=Am×n(jthcolumn)-A~m×n(jthcolumn). Choose the (j+1)th row to be the one containing the maximum element of C (ensuring that the jth row is not chosen again).(6)
Assignujm×1=C.(7)
Update the field approximationA~m×n=A~m×n+ujm×1vj1×n.(8)
Ifujm×1vj1×n<εA~m×n, stop scanning; else increment j and repeat steps (2) to (8).A low rank matrix can be well approximated by the product of two full rank matrices using ACA method:(16)Am×n≈A~m×n=AUm×rAVr×n,where r is the effective rank of matrix and the memory requirement decreased from m×n entries to r(m+n) entries. We can get the block decomposition expression by using the ACA to the iterative process of Uij: (17)Uij=UUUVi,j=ZUZVij-∑p=1i-1UVTUUTpiDpp-1UUUVpj.By using the ACA to the LU solving solution, the induced current densityJi on S can be obtained by the forward iterative steps(18)Xi=Vi-∑p=1i-1UVTUUTpiDpp-1XUXVpand back iterative steps(19)Ji=Dii-1Xi-∑p=i+1nUUUVipJUJVp.


## 4. Numerical Results

In this section, we simulated three numerical examples to demonstrate the accuracy and validity of the proposed method. All the results are computed on an Intel® Core™2 Duo PC with 3.40 GHz processor and 16 GB RAM; only one core is used; the ACA iteration error threshold is 1e-3. The percentage error of advised algorithm is defined as(20)V-Z·JV×100%.(A) A Dielectric Coated PEC Sphere. We start by considering the scattering problem from a dielectric coated PEC sphere with radius 1 λ. Incident frequency is 300 MHz from θ=0° (φ=0°). The thickness of the coating material is τ=0.01 λ. The relative permittivity of LHM is εr=9.0 and relative permeability is μr=1.77-j4.06. The geometry is divided into 2220 triangular patches with an average length of λ/10, and this gives rise to 3330 unknowns. In Figure 2, the bistatic radar cross section (RCS) is shown from θ=0° to θ=180°. The dashed lines by our proposed TDS-LDLT-ACA method agree well with the solid lines by the commercial EM software FEKO.Figure 2
LHM coated PEC sphere.(B) LHM Coated PEC Cylinder. The second example is the problem of scattering from LHM coated PEC cylinder with radius of 0.5 λ and height of 2 λ as shown in Figure 3, the plane wave is incident from θ=0° (φ=0°) to θ=180° (φ=0°), the relative permittivity of LHM is εr=-9.0, and relative permeability is μr=-1.77-j4.06. As a comparison, we have examined the RHM coated case with the relative permittivity εr=9.0 and relative permeability μr=1.77-j4.06. The geometry is divided into 1456 triangular patches with an average length of λ/10, and this gives rise to 2184 unknowns. In ACA operation, the geometry is subdivided into 8 spatial regions as shown in Figure 3. The monostatic RCS calculated for VV polarization is shown in Figure 3 along with the reference solution provided by the FEKO.Figure 3
LHM coated PEC cylinder.In Figure4(a), the coating thickness is 0.01 λ; a quite good agreement between the reference FEKO solution and the results of advised method can be observed. Also, we can see that the monostatic RCS value of coated cylinder with RHM is less than that of LHM coated circumstance in most of angle range.Figure 4
VV polarization monostatic RCS of coated cylinder.
(a)
Coating thickness: 0.005λ
(b)
Coating thickness: 0.01λThe monostatic RCS is shown in Figure4(b) when coating thickness τ reduced to 0.005 λ. We can find that the RCS change rule with incident plane wave angle is similar to Figure 4(a). When coating thickness reduced, the RCS value of RHM drops about 3 dB in most of incident angle area, while the RCS values of LHM increase about 2 dB in most of incident angle area. The phenomenon is also consistent with the conclusion in literature [23] that the reflectivity and RCS will vary with the coating thickness.(C) LHM Coated PEC Missile Model. The third example is LHM coated PEC missile model shown in Figure 5, which has a length of 9.91 m, a wingspan of 6.28 m, and a height of 2.19 m. The geometry is divided into 11300 triangular patches with an average length of λ/10, and this gives rise to 16950 unknowns.Figure 5
LHM coated PEC missile model.The thickness of the material isτ=0.01 λ. The relative permittivity of LHM is εr=-1.5-j0.15 and relative permeability is μr=-1.0-j0.1. The missile model is illuminated by an incident plane wave with the incident direction of θ=90° and φ=270°. Bistatic RCS is shown in Figure 6 along with the reference solution provided by FEKO. We observe that the agreement between the FEKO results and those obtained with the TDS-LDLT-ACA method is excellent. There exists difference between the incident angles 110°<θ<160°; however, the difference of scattering results is relatively smaller.Figure 6
Bistatic RCS of LHM coated PEC missile model.The total CPU time and RAM requirements of above three examples are shown in Table1. It can be seen that the proposed method has compressed much more CPU time and random memory demand than conventional MoM.Table 1
Time and RAM requirements.
The targets andthe unknowns
Number of incident plane waves
Method
RAM requirements/MB
Time/s
Coated sphere, 3330
1
FEKO
146
56
TDS-LDLT-ACA
132
46
LHM coated cylinder, 2184
180
FEKO
118
948
TDS-LDLT-ACA
28
385
LHM coated missile model, 16950
1
FEKO
1347
924
TDS-LDLT-ACA
1355
816


## 5. Conclusion

A TDS approach together with explicit PEC boundary conditions has been proposed to handle the EM scattering problems of thin LHM coating target efficiently. The modeling process has been greatly simplified compared with the surface or volume integral equation approach. The ACA method is incorporated into block LDLT factorization algorithm to reduce the filling time of impedance matrix and speed up the factorization and solving processes. The numerical results demonstrate that the backscattering monostatic RCS values in specific angle area can be reduced by selecting the appropriate LHM coating thickness. The method can be used for electromagnetic scattering problems of complex structure LHM coated PEC target and has an important reference value in practical engineering applications of arbitrary shapes targets.



---

*Source: 2108586-2016-05-30.xml*