International Journal of Antennas and Propagation

Volume 2014 (2014), Article ID 281303, 6 pages

http://dx.doi.org/10.1155/2014/281303

Research Article

The MLFMA Equipped with a Hybrid Tree Structure for the Multiscale EM Scattering

Wei-Bin Kong, Hou-Xing Zhou, Wei-Dong Li, Guang Hua, and Wei Hong

State Key Laboratory of Millimeter Waves, Southeast University, Nanjing 210096, China

Received 21 October 2013; Revised 8 January 2014; Accepted 8 January 2014; Published 19 February 2014

Academic Editor: Ladislau Matekovits

Copyright © 2014 Wei-Bin Kong 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.

Abstract.

We present an efficient strategy for reducing the memory requirement for the near-field matrix in the multilevel fast multipole algorithm (MLFMA) for solving multiscale electromagnetic (EM) scattering problems. A multiscale problem can obviously lower the storage efficiency of the MLFMA for the near-field matrix. This paper focuses on overcoming this shortcoming to a certain extent. A hybrid tree structure for the MLFMA that possesses two kinds of bottom-layer boxes with different edge sizes will be built to significantly reduce the memory requirement for the near-field matrix in the multiscale case compared with the single-tree-structure technique. Several numerical examples are provided to demonstrate the efficiency of the proposed scheme in the multiscale EM scattering.

1. Introduction

The method of moments (MoM) [1] has been widely used due to its superior capability to handle arbitrarily shaped objects in electromagnetic (EM) radiation and scattering problems. The computational complexity of the MoM is for a direct solver like the LU and for an iterative solver like the CG, where is the number of unknowns and is the iteration count.

Because of its higher complexity, the traditional MoM is not suitable for solving electrically large EM problems. The multilevel fast multipole algorithm (MLFMA) [2–4] can reduce the computational complexity of a matrix-vector product from to as well as the storage complexity of the MoM matrix from to . Such an advancement thanks to the application of the addition theorem of the Green’s function and the introduction of the octree structure. To further improve the computational efficiency, the MLFMA can be combined with parallel computing technology [5–11].

In EM simulation, the Rao-Wilton-Glisson (RWG) basis functions [12] are frequently selected to build a MLFMA model for an electrically large object. For this, the surface of the electrically large object is generally discretized with about 10 elements per wavelength. In such a case, the octree whose bottom-layer or finest-layer boxes have edge size of about 0.3 wavelength is more suitable for grouping all basis functions. However, a multiscale problem can obviously lower the storage efficiency of the MLFMA for the near-field matrix, because some local regions of the object surface need to be overmeshed to capture tiny geometry features, leading to a very large near-field matrix if the single octree structure is still used. In order to take care of the fine mesh parts, smaller bottom-layer boxes seem to be considered. However, the bottom-layer boxes with edge size of 0.15 wavelength are not suitable for grouping those RWG basis functions with edge size of about 0.1 wavelength. Besides, the box edge size of 0.15 wavelength is already very close to the subwavelength breakdown of the MLFMA (or the HF-MLFMA) [13–19]. Introducing some assistive technology can alleviate this problem to a certain extent. In [20], Pan et al. combined the interpolative decomposition (ID) with the MLFMA to alleviate the multiscale problem. In [21], Vikram et al. combined the accelerated Cartesian expansion (ACE) with the MLFMA to alleviate the multiscale problem.

In this paper, we focus on overcoming the multiscale problem to a certain extent. A hybrid tree structure for the MLFMA that possesses two kinds of bottom-layer boxes with different edge sizes will be built to significantly reduce the memory requirement for the near-field matrix in the multiscale case compared with the single-tree-structure technique. Compared with the previous methods [20, 21], this method is particularly simple and easy to implement. The idea of this method germinated in our previous work [22]. Several numerical examples are provided to demonstrate the efficiency of the proposed method in the multiscale EM scattering.

This paper is organized as follows. In Section 2, the hybrid tree structure (HTS) is built. In Section 3, the MLFMA equipped with a HTS is presented. In Section 4, numerical examples are provided. Finally, the conclusion is given in Section 5. Besides, in this paper, always denotes the wavelength in the free space.

2. The Hybrid Tree Structure

In the EM simulation using the RWG basis functions, the surface of an electrically large perfectly conducting object is discretized with more than or equal to 10 elements per wavelength. To build a MLFMA model, an octree structure for grouping all basis functions must be constructed. For this, the entire object is enclosed within a large box, called 0-layer box. Then the 0-layer box is divided into its eight children boxes, called 1-layer boxes. Each 1-layer box is divided into its eight children boxes, called 2-layer boxes, until the edge size of the finest-layer boxes is about 0.3 wavelength [2–4]. All basis functions are grouped into the corresponding finest-layer or bottom-layer boxes. Only the nonempty boxes are retained at all levels. In this way, we build a hierarchy of the boxes for the entire object.

The MLFMA splits a matrix-vector product into two parts: near-field interactions and far-field interactions. The near-field interactions between adjacent boxes at the bottom layer form the near-field matrix of the MLFMA. Now, the edge size of the bottom-layer boxes is selected to be . At this time, the number of the RWG basis functions with edge size of about in a bottom-layer box is usually less than 50. However, if some local regions of the object are overmeshed, for example, if the discrete sizes are less than or equal to , then the numbers of the RWG basis functions in some bottom-layer boxes can exceed 100. This is why a multiscale problem can lead to a very large near-field matrix. Thus, it is a meaningful work to build a proper tree structure for the MLFMA in the multiscale case.

As elucidated in [18], in most cases, the numerical error produced by bottom-layer-box size is acceptable. If the box size is below , the MLFMA could suffer from subwavelength breakdown [13–19]. Hence, is an appropriate choice for the bottom-layer-box size when the discrete sizes of some local regions are less than or equal to . According to the different surface discretization, a hybrid of the octree structure and the 27-tree structure will be a better choice for the MLFMA in the multiscale case. In such a hybrid tree structure, there are two kinds of bottom-layer boxes. One has the edge size of wavelength and another has the edge size of wavelength, as shown in Figure 1.

Figure 1: Two kinds of the bottom-layer boxes. For example, each purple or black box is a bottom-layer box in the octree structure, while each blue box is a bottom-layer box in the 27-tree structure.

The formation process of the HTS is closely related to the surface discretization and difficult to be described with a few formulas, but this process is very easy for a computer language to describe. Here, our purpose is just to show how a hybrid tree structure is applied in the MLFMA. Then, the situation in Figure 1 can be used as a simple example. In this case, the formation procedure can be outlined as follows.

(1)The entire object is enclosed within a large box, called 0-layer box.(2)Then the 0-layer box is divided into its eight children boxes, called 1-layer boxes. Each 1-layer box is divided into its eight children boxes, called 2-layer boxes, until the edge size of the bottom-layer boxes is about .(3)Group all basis functions into the corresponding current bottom-layer boxes and record the number of the basis functions in each current bottom-layer box.(4)If one bottom-layer box contains more than basis functions, then mark its parent box. Redivide this parent box marked into children boxes with edge size of .(5)After the process described by for all bottom-layer boxes is finished, some basis functions have temporarily lost the bottom-layer boxes containing them. Regroup these basis functions into the corresponding new bottom-layer boxes whose edge size is .

In this way, we build a hierarchy of boxes for all basis functions. Clearly, the finest layer in the above hybrid tree structure is an imprecise concept. In fact, there are two kinds of the bottom-layer edge sizes, that is, and . Figure 2 shows such a tree structure for the MLFMA.

Figure 2: The hybrid tree structure in the 2D case.

3. The MLFMA Equipped with a HTS

Two boxes at the same layer are called near-neighbours if they have at least one common point. In the following, all boxes mentioned are not empty. For convenience, we also need some notations.

Definition 1. Let “ob” be an observed box; then is a collection that consists of “ob” and all its near-neighbours. Here, the subscript in means “near-neighbours.”

Definition 2. Let “ob” be an observed box; then is a collection that consists of those children of all near-neighbours of ob parent box that are well separated from ob. Here, the subscript in means “next near-neighbours.” Figure 3 shows the members of and related to a bottom-layer box ob in the HTS.

Figure 3: An observed box (ob) at the bottom layer and all its near-neighbours (nb) and all next near-neighbours (nnb).

As usual, the near-field matrix formed by the interactions between the near-neighbours at the bottom layer is directly calculated and stored in advance. The product of the near-field matrix with a vector is always directly calculated.

Assume that the RWG functions and belong to the next near-neighbours and , respectively. The location vectors of the centers of and are denoted by and , respectively. The location vectors of any points on supports of and are denoted by and , respectively. Furthermore, define , , and . Then, the far-matrix elements formed by the interactions between these two next near-neighbours for the combined field integral equation (CFIE) are expressed as

where

where is the combination factor and . In this paper, is selected to be .

The formula (1) can be decomposed into three processes: the aggregation process, the translation process, and the disaggregation process. For the sake of clarity, an example in Figure 4 is considered. Compared with the conventional MLFMA, the MLFMA equipped with a HTS contains more details about the aggregation, translation, and disaggregation at the bottom layer.

Figure 4: The process of aggregation, translation, and disaggregation, where each “” represents the center of a bottom-layer box.

In Figure 4, the solid purple arrow line illustrates the aggregation, translation, and disaggregation processes between two next near-neighbours at the bottom layer of different sizes, whose centers are marked with and , respectively. There is an indirect transmission of information between two basis functions marked with and , respectively. That is, . Two solid black arrow lines illustrate the aggregation, translation, and disaggregation processes between two next near-neighbours at the bottom layer of the same sizes. It should be noted that the higher layer calculation is the same as the conventional MLFMA.

4. Numerical Examples

In this section, the conductor surface is meshed with triangular elements and RWG functions [12] are chosen as both basis and test functions, where the edge sizes of the triangular elements are from to . The singularity of the near-field matrix elements is dealt with through the approach in [23].

4.1. A PEC Sphere with Radius 2.0 λ

A perfectly electric conducting (PEC) sphere with radius size is considered in Figure 5. The upper hemisphere surface is discretized with edge size of about , and the lower hemisphere surface is discretized with edge size of about , producing 47148 unknowns.

Figure 5: Mesh of a PEC sphere with radius .

Figure 6 depicts the bistatic RCS curves obtained by the Mie solution [24], the conventional MLFMA, and the MLFMA equipped with a HTS. The numbers of the bottom-layer boxes in the conventional MLFMA and the MLFMA equipped with a HTS are 784 and 1273, respectively, as shown in Table 1. The memory requirement for both the fast methods is 896.796 MB and 450.876 MB, respectively, and hence the new method leads to 50.28% memory reduction.

Table 1: Comparison between the two MLFMAs in Section 4.1.

MLFMA
MLFMA equipped with a HTS

Number of levels
4
4

Number of bottom-layer boxes
784
1273

Memory for the near-field matrix
896.796 MB
450.876 MB

Figure 6: Bistatic RCS curves of the PEC sphere.

4.2. A PEC Rectangular Plate with Thin Slots

Next, we consider the EM scattering from a PEC rectangular plate with thin slots, as shown in Figure 7. The plate surface is discretized into a lot of triangles, yielding 15153 unknowns. Because there are tiny geometric details, the resulting mesh is a multiscale mesh and also an extremely nonuniform mesh. The longest edge size is about , while the shortest edge size is about .

Figure 7: Mesh of a PEC rectangular plate with thin slots. The dimension of the plate is by and that of each slot is by .

Figure 8 depicts the bistatic RCS curves obtained by the MoM, the conventional MLFMA, and the MLFMA equipped with a HTS. The numbers of the bottom-layer boxes in the conventional MLFMA and in the MLFMA equipped with a HTS are 378 and 473, respectively, as shown in Table 2. The memory requirement for both the fast methods is 104.76 MB and 81.85 MB, respectively, and hence the new method leads to 21.9% memory reduction.

Table 2: Comparison between the two MLFMAs in Section 4.2.

MLFMA
MLFMA equipped with a HTS

Number of levels
6
6

Number of bottom-layer boxes
378
473

Memory for the near-field matrix
104.76 MB
81.85 MB

Figure 8: Bistatic RCS curves of the PEC plate with thin slots.

4.3. A PEC Aircraft-Like Model

The last structure is a PEC aircraft-like model in Figure 9. The width is and the length is . The head part of flyer is discretized with edge size of about , and the remaining part is discretized with edge size of about . The number of unknowns is 54159.

Figure 9: Mesh of the PEC aircraft-like model.

In Figure 10, the bistatic RCS results are plotted, which is obtained by the conventional MLFMA and the MLFMA equipped with a HTS. The numbers of the bottom-layer boxes in the conventional MLFMA and the MLFMA equipped with a HTS are 1793 and 1902, respectively, as shown in Table 3. The memory requirement for both the fast methods is 494.83 MB and 409.97 MB, respectively, and hence the new method leads to 17.15% memory reduction.

Table 3: Comparison between the two MLFMAs in Section 4.3.

MLFMA
MLFMA equipped with a HTS

Number of levels
7
7

Number of the bottom-layer boxes
1793
1902

Memory for the near-field matrix
494.83 MB
409.97 MB

Figure 10: Bistatic RCS curves of the PEC aircraft-like model.

5. Conclusions

In the multiscale case, a single octree structure for the MLFMA can produce a very large near-field matrix. In this paper, in order to overcome this shortcoming to a certain extent, a hybrid tree structure scheme has been established. In this scheme, there are two kinds of bottom-layer boxes with different edge sizes. With the hybrid tree structure, the memory requirement for the near-field matrix can be significantly reduced in the multiscale case compared with the MLFMA equipped with a single octree structure.

It should be pointed out that to what extent the total amount of storage is reduced by the proposed method is closely related to the proportion of the local regions with the tiny geometry features to the whole surface domain. This fact has been revealed by the three numerical examples provided in this paper. Therefore, in the general case, to what exact extent the proposed method can reduce the total amount of storage can not be determined in advance.

Conflict of Interests

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

Acknowledgment

This work was supported by the National Basic Research Program of China (2013CB329002, 2010CB327400, and 2009CB320203).

References

R. F. Harrington, Field Computation by Moment Methods, The MacMillian, New York, NY, USA, 1968.

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.

J.-M. Song, C.-C. Lu, and W. C. Chew, “Multilevel fast multipole algorithm for electromagnetic scattering by large complex objects,” IEEE Transactions on Antennas and Propagation, vol. 45, no. 10, pp. 1488–1493, 1997.

W. C. Chew, J. M. Jin, E. Michielssen, and J. M. Song, Fast and Efficient Algorithms in Computational Electro-Magnetics, Artech House, Boston, Mass, USA, 2001.

S. Velamparambil and W. C. Chew, “Analysis and performance of a distributed memory multilevel fast multipole algorithm,” IEEE Transactions on Antennas and Propagation, vol. 53, no. 8, pp. 2719–2727, 2005.

Ö. Ergül and L. Gürel, “Efficient parallelization of the multilevel fast multipole algorithm for the solution of large-scale scattering problems,” IEEE Transactions on Antennas and Propagation, vol. 56, no. 8, pp. 2335–2345, 2008.

L. Gürel and Ö. Ergül, “Hierarchical parallelization of the multilevel fast multipole algorithm (MLFMA),” Proceedings of the IEEE, vol. 101, no. 2, pp. 332–341, 2013.

X. M. Pan, W. C. Pi, M. L. Yang, Z. Peng, and X. Q. Sheng, “Solving problems with over one billion unknowns by the MLFMA,” IEEE Transaction on Antennas and Propagation, vol. 60, no. 5, pp. 2571–2574, 2012.

V. Melapudi, B. Shanker, S. Seal, and S. Aluru, “A scalable parallel wideband MLFMA for efficient electromagnetic simulations on large scale clusters,” IEEE Transactions on Antennas and Propagation, vol. 59, no. 7, pp. 2565–2577, 2011.

B. Michiels, J. Fostier, I. Bogaert, and D. De Zutter, “Weak scalability analysis of the distributed-memory parallel MLFMA,” IEEE Transaction on Antennas and Propagation, vol. 61, no. 11, pp. 5567–5574, 2013.

J. Guan, Y. Su, and J. M. Jin, “An openMP-CUDA implementation of multilevel fast multipole algorithm for electromagnetic simulation on multi-GPU computing systems,” IEEE Transaction on Antennas and Propagation, vol. 61, no. 7, pp. 3607–3616, 2013.

S. M. Rao, D. R. Wilton, and A. W. Glisson, “Electromagnetic scattering by surfaces of arbitrary shape,” IEEE Transactions on Antennas and Propagation, vol. AP-30, no. 3, pp. 409–418, 1982.

J. S. Zhao and W. C. Chew, “Three-dimensional multilevel fast multipole algorithm from static to electrodynamic,” Microwave and Optical Technology Letters, vol. 26, no. 1, pp. 43–48, 2000.

H. Wallén, S. Järvenpää, and P. Ylä-Oijala, “Broadband multilevel fast multipole algorithm for acoustic scattering problems,” Journal of Computational Acoustics, vol. 14, no. 4, pp. 507–526, 2006.

H.-W. Wallén and J. Sarvas, “Translation procedures for broadband MLFMA,” Progress in Electromagnetics Research, vol. 55, pp. 47–78, 2005.

H.-W. Cheng, W. Y. Crutchfield, Z. Gimbutas et al., “A wideband fast multipole method for the Helmholtz equation in three dimensions,” Journal of Computational Physics, vol. 216, no. 1, pp. 300–325, 2006.

L.-J. Jiang and W. C. Chew, “Low-frequency fast inhomogeneous plane-wave algorithm (LF-FIPWA),” Microwave and Optical Technology Letters, vol. 40, no. 2, pp. 117–122, 2004.

L.-J. Jiang and W. C. Chew, “A mixed-form fast multipole algorithm,” IEEE Transactions on Antennas and Propagation, vol. 53, no. 12, pp. 4145–4156, 2005.

H. Wallén, S. Järvenpää, P. Ylä-Oijala, and J. Sarvas, “Broadband Müller-MLFMA for electromagnetic scattering by dielectric objects,” IEEE Transactions on Antennas and Propagation, vol. 55, no. 5, pp. 1423–1430, 2007.

X.-M. Pan, J.-G. Wei, Z. Peng, and X.-Q. Sheng, “A fast algorithm for multiscale electromagnetic problems using interpolative decomposition and multilevel fast multipole algorithm,” Radio Science, vol. 47, no. 1, Article ID RS1011, 2012.

M. Vikram, H. Huang, B. Shanker, and T. Van, “A novel wideband FMM for fast integral equation solution of multiscale problems in electromagnetics,” IEEE Transactions on Antennas and Propagation, vol. 57, no. 7, pp. 2094–2104, 2009.

W.-B. Kong, H.-X. Zhou, W.-D. Li, G. Hua, and W. Hong, “Using a hybrid tree structure in the MLFMA in the multiresolution case,” in International Conference on Computational Problem-Solving (ICCP '12), pp. 2–8, 2012.

P. Ylä-Oijala and M. Taskinen, “Calculation of CFIE impedance matrix elements with RWG and n × RWG functions,” IEEE Transactions on Antennas and Propagation, vol. 51, no. 8, pp. 1837–1846, 2003.

J. J. Bowman, T. B. A. Senior, and P. L. E. Uslenghi, Electromagnetic Scattering by Simple Shapes, Revised, Hemisphere, New York, NY, USA, 1987.