#### Abstract

The computation of scattering from multilayer dielectric bodies is studied by using the combined tangential formulation (CTF) of surface integral solution. A simple and efficient preconditioner is designed for the surface integral solution of multilayer dielectric bodies and validated by numerical experiments. Compared with the traditional near field preconditioner, the proposed preconditioner significantly reduce CPU time and memory requirement. Furthermore, the multilevel fast multipole algorithm (MLFMA) is employed to improve the capability of the solutions. The trick of efficiently implementing MLFMA is presented for multilayer dielectric bodies. Numerical examples are presented to verify the accuracy and efficiency of the approach for computing scattering from multilayer dielectric problems.

#### 1. Introduction

Electromagnetic scattering from multilayer dielectric bodies is an important problem in the computational electromagnetics (CEM) since it is widely used in real-life applications such as Luneburg lens, radomes, where each dielectric is fully included in the outer one. The method of moment (MoM) based on surface integral equations (SIE) has shown to be an efficient solution of this problem, since the volume discretization is avoided in this solution. The formulations of SIE have been well studied for homogeneous bodies in the literatures [1–5]. However, the numerical performance of surface integral solutions is not well investigated for multilayer dielectric bodies.

It is known that the preconditioners are usually required for obtaining an efficient surface integral solution of homogeneous bodies [6–8]. There are many preconditioners. A conventional preconditioner is constructed by directly employing the inverse of near field matrix, which is exactly compatible with the multilevel fast multipole algorithm (MLFMA). But it is very time-consuming and large memory requirement. Another general preconditioner is the block-diagonal preconditioner [9], which is constructed by choosing the block-diagonal part of the near-interaction. However, it does not work well for CTF solution as verified in literature [5]. A more efficient preconditioner is the Schur complement preconditioner [7]. Since the discretized matrix of multilayer dielectric problems contains both the interactions of equivalent electric and magnetic currents on the same layer and those on different layers, it is hard to efficiently extend the Schur complement preconditioner from homogeneous bodies to multilayer bodies. In this letter, the combined tangential formulation (CTF) is employed to formulate scattering from multilayer dielectric bodies due to its high accuracy and efficient memory. Furthermore, CTF is free of internal resonances [4]. A simple preconditioner is designed for the efficient CTF solution of multilayer bodies. Compared with the existing preconditioners [6–8], it is easier to be implemented to the multilayer dielectric bodies. MLFMA [10, 11] is efficiently incorporated into the solution to improve its efficiency. Numerical experiments are presented to investigate the performance of the proposed approach.

#### 2. Theory and Formulation

##### 2.1. Surface Integral-Equation Formulations

Consider electromagnetic scattering from a multilayer dielectric object placed in free space and illuminated by an incident plane wave (). The exterior free space region is denoted as and the interior regions are denoted as with for dielectric layers from outside to inside, with each dielectric being fully included in the outer one, as illustrated in Figure 1. The permittivity and permeability of region are and , respectively. The interface between two regions and is denoted by and denotes the unit normal of pointing toward the interior of . The equivalent electric and magnetic currents on are denoted as and , respectively.

The electric field integral equation in the region can be formulated as [12]The magnetic field integral equation in the region can be formulated aswhere , denote the total field in and is the solid angle which equals for common smooth surface. for and for other cases. The operators and are defined asHere indicates the principal value of the integral, is Green’s function in region , and .

By combining the equations in region with those in region , we obtain the following CTF formulation [5]:where , the operator , and , .

##### 2.2. Discretization of Surface Integral-Equation Formulations

To discretize (4) and (5), each dielectric surface/interface is meshed by planner triangular patches. The surface currents are expanded with Rao-Wilton-Glisson (RWG) basis functions [13]. Applying Galerkin’s testing method, we obtain the following discretized matrix equation: where and and are the unknowns of electric and magnetic currents, respectively. To be more specific, the submatrices of , , , have the following structures:wherewhere denotes the RWG basis and testing functions on and stands for the inner product.

#### 3. Implementation of MLFMA for Multilayer Bodies

To efficiently solve (7), MLFMA is employed to speed up the matrix-vector multiplication in iterative solutions. Usually, each sub-matrix-vector multiplication of (7) is independently speeded up by using MLFMA. Each dielectric layer shares the same cubic clusters, but different truncation numbers in MLFMA are employed. However, if the special structure of (9) is employed, a more efficient MLFMA implementation can be achieved. For example, since the far-field interactions of and in (7) act on the same equivalent electric currents and also Green’s function for and is the same, the aggregation and translation components for their far-field interactions are the same; only the receiving components are different. Thus, the calculation of the far-field interactions of and in each matrix-vector multiple of the iteration can be finished with only once computation of the aggregation and translation and twice disaggregations. The same trick can be applied to and , and , and , and , and , and , and and .

#### 4. Construction of Preconditioner

Preconditioners are usually required for iterative solutions of (7) since the increase of the number of dielectric layers leads to the rapid growth of the condition number of impedance matrix. Employing a preconditioner on the matrix equation of (7) yields A conventional preconditioner is constructed by employing the inverse of the near field matrix, which is called near field preconditioner (NFP). NFP is exactly compatible with MLFMA. However, the construction of NFP is time-consuming and large memory is required. More efficient preconditioners are proposed by making use of the approximate Schur complement [7]. Since (7) is not a partitioned system [14, 15], these preconditioners are difficult to be extended from homogeneous bodies to multilayer dielectric object.

In this letter, a simple preconditioner is designed for the iterative solution of (7). In the construction of NFP, the near field matrix is obtained by maintaining the entries of self-interaction in one box and near-interaction with the neighbor boxes and omitting other entries of the original matrix. The boxes are ones in the lowest level of MLFMA. Since the number of the entries of near-interaction is quite large, NFP is resource-consuming, especially for multilayer dielectric targets since multilayer dielectric bodies have interactions between equivalent currents on different interfaces. To overcome such drawbacks, we employ the interaction distances of entries between source and field as criterion to determine whether the entries are retained or omitted, as shown in where and represent the location of the testing and bases elements. To be more specific, if the interaction distance between source and field is less than , the matrix entries are retained; otherwise the matrix entries are omitted. Here is the free space wavelength. For multilayer bodies, the distance between source and basis functions is important because the interactions between different layers are stronger when the layers are more close to each other. The maintained entries include interactions between different layers if the distance between layers is small. Thus the sparsified scheme can retain the main interactions between source and basis functions. Since only those entries within self-elements and neighbor elements are maintained in this sparsified scheme, the number of the maintained entries is much smaller than that in NFP. We employ the inverse of this sparsified matrix as preconditioner and call it distance sparse preconditioner (DSP).

#### 5. Numerical Results

In this section, we present several examples to demonstrate the accuracy and efficiency of the presented approach. All the simulations are at frequency of 300 MHz. The generalized minimum residual (GMRES) solver is used and a residual error of is set to terminate iterations. MUMPS is used to calculate the inverse of the sparsified matrix in the preconditioner construction. MUMPS, short for “Multifrontal Massively Parallel Solver,” is a direct sparse solver based on the parallel multifrontal method [16].

First, we consider a dielectric multilayer body consisting of 5 concentric spheres with radius of , , , , and for each interface from inside to outside. The relative permittivity of regions , , , , and is 3, 2, 3, 2, and 1, respectively. The body is illuminated by a plane wave propagating in direction with the electric field polarized in direction. The body is meshed with element size and has 16,128 unknowns. Figure 2 presents -polarized bistatic RCS in -plane for this multilayer sphere. It shows that the numerical results agree well with the analytical Mie-series solution. The numerical performance of MLFMA with NFP and DSP and no-preconditioner (No PC) is presented in Table 1. The CPU time and the memory requirement in Table 1 are those for the analysis and factorization step of MUMPS for constructing the inverse matrix of the preconditioners. It can be seen that iteration number is greatly reduced by using preconditioners. The total time required by DSP is only about 1/4 of that required by NFP. Furthermore, the memory required by DSP is less than 1/5 of that required by NFP. In addition, the efficient MLFMA implementation presented in this letter can reduce the iteration CPU time from 2,000 s required by the conventional MLFMA implementation to 1,060 s.

To investigate the numerical performance of the solution for bodies with higher dielectric constant, the same multilayer sphere but with different permittivity of the layer is calculated. The parameters of are listed in Table 2. The element sizes used in the meshing and with different are shown in Table 2. The resulted numbers of unknowns are also presented in Table 2. The total solution time as a function of is plotted in Figure 3. It can be seen that DSP has the least total solution time for all dielectric constants. Furthermore, DSP has more advantage over NFP with higher dielectric constant.

The next three examples are computed to further demonstrate the capability of the presented approach for multilayer dielectric objects with different shapes, compared with the commercial software of FEKO. The first object is a cubic shell having three dielectric layers of from inside to outside. Each of the layers has a thickness of and the side length of the cubic shell is , as shown in Figure 4(a). The average mesh size is and the number of unknowns is 122,544. The bistatic RCS patterns for - and -polarizations are computed. Figures 4(b) and 4(c) present the comparison of bistatic RCS obtained by DSP with the results from FEKO, where the mesh size is also set at and the residual error is . It is shown that the computational values of DSP are in agreement with results from FEKO.

**(a)**

**(b)**

**(c)**

The second example is an ellipsoidal shell with radius along -, -, and -axis being , , and , respectively, as shown in Figure 5(a). It consists of three dielectric layers of from inside to outside. The thickness of the innermost layer is and the others are both . Bistatic RCS obtained using DSP with 45,558 unknowns is compared with the result from FEKO in Figures 5(b) and 5(c) and good agreement is observed.

**(a)**

**(b)**

**(c)**

The last example is a cylindrical cavity as shown in Figure 6(a), which is a concave object. The radius of the cylindrical cavity at the base is and the height is . It consists of three dielectric layers; each has thickness of . The middle layer has and the other two layers have . Bistatic RCS obtained using DSP with 80,250 unknowns is compared with the result from FEKO in Figures 6(b) and 6(c) and good agreement is observed. The computation information of these three examples is listed in Table 3. The presented solution is more efficient than FEKO.

**(a)**

**(b)**

**(c)**

#### 6. Conclusions

The surface integral solution of scattering from multilayer dielectric bodies is studied by using CTF. A simple preconditioner of DSP is presented, and its efficiency is validated by numerical experiments. Compared with the conventional near field preconditioner, the proposed preconditioner of DSP greatly reduce computation time and memory requirement. An efficient MLFMA implementation is given for multilayer dielectric bodies. Numerical results show that the presented MLFMA implementation can save half iteration time compared with the conventional MLFMA implementation. Numerical results show that the presented CTF solution is much more efficient than FEKO for multilayer dielectric bodies.

#### Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This work was supported by the National Basic Research Program (973) under Grants no. 61320602 and no. 61327301, the 111 Project of China under Grant B14010, and the NSFC under Grant no. 61421001.