# Enhancing the Efficiency and Accuracy of MLFMA-PO Hybrid Method for Analyzing Electrically Large Objects

Authors: Hao Ding; Yong Zhang; Hai Lin
Journal: International Journal of Antennas and Propagation
 (2015)

---

## Abstract

An enhanced hybrid technique combining the multilevel
fast multipole algorithm (MLFMA) and the physical optics
(PO) is presented for analyzing scattering and radiation by
electrically large structures. Based on the previously proposed
hybrid technique, an extra octree structure is built for the
moment method (MoM) region to accelerate the solving process
and an iteration scheme is proposed to further enhance the
efficiency compared to the conventional hybrid method. What
is more, the modified iterative physical optics (IPO) accelerated
by MLFMA technique is adopted to improve the accuracy by
simulating complex structures with strong reflections. Numerical
results are given to demonstrate the improved accuracy and
higher accuracy of the proposed hybrid technique compared to
the previous MLFMA-PO method.


---

## Body

## 1. Introduction

Nowadays, simulations of large-scale electromagnetic problems using hybrid techniques which combine the accurate numerical methods and asymptotic high-frequency methods have achieved more and more attention in the area of computational electromagnetics (CEM). Among them, the hybrid technique which combines the integral numerical method based on the conventional MoM [1] and the high-frequency method based on physical optics (PO) is very popular and efficient. There are usually two different ways to implement the MoM-PO: (1) considering the impedance matrix of the MoM by coupling the PO contributions [2, 3]; (2) iterating the voltage matrix of the MoM by PO currents [4–6]. The former one is self-consistent while the iterative MoM-PO is more efficient without losing any accuracy. Therefore, the iterative MoM-PO is widely studied and used in the CEM area.Since the accuracy of conventional MoM-PO solution is limited, one apparent idea is to enlarge the area of the MoM region at the cost of losing efficiency. More unknowns in MoM region will not only cause a heavy burden for the local solver but also increase the complexity of the coupling between the different regions. Several algorithms have been proposed to alleviate this problem such as high-order hybrid MoM-PO techniques [7], multilevel fast multipole algorithm (MLFMA) [8], precorrected-FFT (P-FFT) [9], and adaptive integral method (AIM) [10]. Among them, the MLFMA-PO hybrid technique [11, 12] seems to be the most efficient one since the MLFMA is proved to be the most efficient fast algorithm for 3D arbitrary objects. In our previous work, MLFMA is adopted to accelerate not only the solution of the matrix equations in the MoM region, but also the interactions between the MoM region and PO region. By using this scheme, the computational complexity can be reduced to less than O(Nlog⁡(N)).Many correction techniques have been proposed to improve the hybrid MoM-PO method since the accuracy of PO method is not sufficiently enough for the analysis of complex structures. Among them, two main methods are introduced here. The first approach is to add certain correction terms on the wedges [13–15] which take the influences of sharp edges into consideration. The other approach is the iterative PO method [16, 17] which considers the mutual influences in the PO region. The second approach is simple and suitable for structures with strong reflections such as concavities.In this paper, we are aiming at enhancing the efficiency and accuracy of the previous MLFMA-PO technique for simulating the radiation and scattering problems in terms of electrically large structures based on [12]. Since the burden of the matrix-vector-product (MVP) based on the octree for the whole structure is heavy and unnecessary for the MoM region, an extra octree is built for the MoM region to decrease the calculation time of the MVP in the MoM region. Based on the additional octree, an iterative scheme is proposed to further enhance the efficiency and the detailed explanation is given. Moreover, due to the limited accuracy of PO when facing the complex structures with strong reflections, the IPO is developed to improve the accuracy by considering the interactions in the high-frequency domain inspired from [16, 17]. Several numerical results are given to demonstrate that the proposed method equipped with these corrections is more accurate and efficient compared to the conventional MLFMA-PO technique.


## 2. Formulations

Firstly, a brief review on the hybrid MoM-PO is given for completeness. As shown in Figure1, the structure is first divided into two parts as follows: the MoM region and PO region, where JMoM and JPO denote the electric currents in the MoM and PO region, respectively. The coupled electric field integral equation (C-EFIE) for the MoM region (SMoM) can be built as [2] (1)n^×LEJMoM+LEJPO=n^×Einc,where n^ denotes the unit vector in the MoM region surface SMoM and Einc stands for the electric field of the incident plane wave if it exists. The integral operator LE is defined as(2)LEJ=-jkZ∫‍J+1k2∇∇′·JGdr′,where k represents the free space wave number and Z=μ0/ɛ0 is wave impedance. G(r,r′)=e-jkr-r′/(4πr-r′) is free space Green’s function. The electric currents in PO region which are induced by the MoM currents and incident wave (if it exists) can be written as(3)JPO=2δn^×LHJMoM+Hinc,where δ is the shadowing effect coefficient, which is 1 in the lit PO region and 0 in the shadowed PO region. Hinc denotes the magnetic field of the incident wave. The integral operator LH is defined as(4)LHJ=-∫‍J×∇Gdτ′.Figure 1
The structure is divided into MoM region and PO region.We use the Rao-Wilton-Glisson (RWG) basis function [18] and Galerkin’s testing procedure to derive the matrix equation which is shown as(5)(ZMoM+AMoM-POBPO-MoM)IMoM=VMoM-AMoM-POIPO,where IMoM is the current weights in the MoM region, ZMoM stands for the impedance matrix in the MoM region, AMoM-PO and BPO-MoM represent the coupling matrix between the MoM region and the PO region, and the VMoM denotes the excitations for the MoM region. IPO stands for the currents in PO region induced by the incident plane wave if it exists. The element expressions of matrices ZMoM and AMoM-PO are the same as the conventional MoM equations with RWG basis functions. However, the element expression of matrix BPO-MoM is a little different which is written as(6)BPO-MoM,ij=δt^i++t^i-·n^i×∇∫Tj±Gr,r′fjr′dτ′,where t^i+ and t^i- are the two unit vectors in the middle of the ith edge which are perpendicular to the edge lying in the plane of the triangular pairs Ti± in the PO region. Tj± is the triangular pairs of the jth RWG basis function.


### 2.1. Extra Octree Structure

In the previous implementation of MLFMA-PO hybrid technique [12], an octree structure is built based on the whole geometry information. The C-EFIE is constructed on the MoM region and the impedance matrix consists of influences from itself and the PO region which are described by Z, A, and B. Z describes the self-interactions in the MoM region and the currents in the PO region are not needed if the MVP process is performed on Z. In fact, all the MVP processes are accelerated by MLFMA based on the total octree in the previous hybrid technique which is time-consuming and unnecessary. What is more, the number of unknowns in the MoM region is always much smaller than the one in PO region, and the scale of the octree built for the MoM region will be much smaller compared to the whole one. Thus, higher efficiency can be achieved in the self-interactive MVP compared to the previous MLFMA-PO hybrid technique. Even though a little more memory usage is needed for the extra tree structures, transfer matrix, interpolation matrix, and the radiation and receiving patterns in the lowest level, it is negligible compared to the whole octree structure. It should be noted that the near-field matrix for the MoM region is also different from the conventional MLFMA-PO technique because of the different near-tree judgements.


### 2.2. Iteration/Noniteration Analysis

An iterative scheme for the MoM-PO hybrid technique is proposed in [4] which is faster but in the same level of accuracy compared to the conventional MoM-PO method. The formulation is a bit different from (5) which is shown as(7)ZMoMIMoMi+1=VMoM+AMoM-POIPO+AMoM-POBPO-MoMIMoMi,where IMoMi is the current weights in the MoM region in the ith iteration process. An iterative process will be needed to derive the solution of (7) until the currents in the MoM region become stable. The error ξi in the ith step is defined as(8)ξi=IMoMi-IMoMi-1IMoMi-1,where the norm is defined in the 2-norm form and IMoMi stands for the current coefficients in the ith iteration. The threshold value of the relative residual error ξT and the maximum number of iteration Niter are also needed to control the iteration computation.The previous MLFMA-PO technique [12] considers the impedance matrix of the MoM by coupling the PO contributions while the iterative method updates the voltage matrix of the MoM by PO currents until it reaches stability. If only one whole octree is built, the computational complexity of the MVP processes of Z, A, and B is nearly in the same level since the MLFMA operations are all implemented based on the whole octree. Moreover, only one solving process is needed for the conventional MLFMA-PO technique while several solving processes are needed for the iterative MLFMA-PO method. Thus, the previous noniterative method may be a better choice if only one octree exists. However, the MVP process of Z based on the newly created octree structure is much faster than the other two if the extra octree for the MoM region is built. In the iterative MLFMA-PO method, several MVP processes of A and B are replaced by Z which makes the iterative method a more efficient solver. Then, the iterative MLFMA-PO method is suggested if two octrees are built which is also demonstrated in the numerical results.An explanation is also given in the mathematical manner. Assume that the computation time of the MVP scanning the whole octree isα while the value is β for the octree of the MoM region. Apparently, β is much smaller than the α. The number of the solving iterations for the conventional MLFMA-PO is defined as n1 while the average value is n2 (n2≈n1) for the iterative method with m (m is always smaller than 10) iterations between the MoM region and PO region. Thus, the computation time χmlfma-po of the conventional MLFMA-PO method is(9)χmlfma-po=3n1α.Since the value ofβ is negligible compared to α, the computation time for the iterative MLFMA-PO method can be defined as(10)χiter-mlfma-po=mn2β+2mα≈mn1β+2α≈2mα.From the above two equations, the computation time of the iterative MLFMA-PO method based on two octrees can be estimated as(11)χiter-mlfma-po≈2m3n1χmlfma-po.Sincem is always much smaller than n1, the computation time of the iterative MLFMA-PO based on two octrees will be much smaller than the previous MLFMA-PO method [12] using the same hardware.


### 2.3. IPO Process

In the real-world problems, the structures are always complicated and nonsmooth. Thus, the currents derived by the PO method may not be accurate enough. This could make the precision of the total solution limited. In particular for the structures with strong reflections in the PO region, the conventional MLFMA-PO technique will produce totally wrong solutions in some angles. This problem can be alleviated by introducing the iterative process for the PO region (IPO) to improve the accuracy of the PO currents.The basic MoM-IPO principle can be sought in [16, 17] and it is not presented here for simplicity. However, the modified IPO method proposed in this paper is a little different from the conventional IPO. Firstly, the RWG basis functions are used to expand the electric currents to maintain the continuity. Secondly, MLFMA is also adopted to decrease the computational complexity of the mutual interaction in PO region to O(Nlog⁡(N)). The proposed IPO process can be defined as(12)IPOi+1=CPOIPOi+IPO0,where IPOi stands for the ith currents in the PO region, IPO0 represents the initial currents produced by the incident wave and the currents in the MoM region, and the matrix CPO is the coupling matrix in the PO region. The elements in the CPO can be expressed by(13)CPO,ij=0,i=jδ(t^i++t^i-)·n^i×∇∫Tj±‍Gr,r′fjr′dτ′,i≠j,where t^i+ and t^i- are the two unit vectors which are perpendicular to the ith edge lying in the plane of the triangular pairs Ti± in the PO region. δ denotes the shadowing effects coefficient and Tj± is the triangular pairs of the jth RWG basis function. Note that the diagonal elements of the coupling matrix are zero since the self-interactions are not considered. In this paper, the MVP process of CPOIPOi in each iteration is accelerated by MLFMA technique to maintain the computational complexity to be O(Nlog⁡(N)). Compared to the conventional MLFMA for the PEC objects, implementation for IPO is straightforward except for the calculation of near-field impedance matrix and lowest-level receiving patterns. The near-field impedance matrix element can be computed based on (13). However, the lowest-level receiving patterns need some corrections compared to the conventional MLFMA which is shown as(14)Rfimk^=-δk^×t^i+×n^i++t^i-×n^i-e-jk·rim.What is more, the extra octree for the PO region is not suggested to build because the PO region always plays a dominant role in the whole structure. The complexity of MVP based on the octree for the PO region will not be decreased apparently compared to the whole one.


## 2.1. Extra Octree Structure

In the previous implementation of MLFMA-PO hybrid technique [12], an octree structure is built based on the whole geometry information. The C-EFIE is constructed on the MoM region and the impedance matrix consists of influences from itself and the PO region which are described by Z, A, and B. Z describes the self-interactions in the MoM region and the currents in the PO region are not needed if the MVP process is performed on Z. In fact, all the MVP processes are accelerated by MLFMA based on the total octree in the previous hybrid technique which is time-consuming and unnecessary. What is more, the number of unknowns in the MoM region is always much smaller than the one in PO region, and the scale of the octree built for the MoM region will be much smaller compared to the whole one. Thus, higher efficiency can be achieved in the self-interactive MVP compared to the previous MLFMA-PO hybrid technique. Even though a little more memory usage is needed for the extra tree structures, transfer matrix, interpolation matrix, and the radiation and receiving patterns in the lowest level, it is negligible compared to the whole octree structure. It should be noted that the near-field matrix for the MoM region is also different from the conventional MLFMA-PO technique because of the different near-tree judgements.


## 2.2. Iteration/Noniteration Analysis

An iterative scheme for the MoM-PO hybrid technique is proposed in [4] which is faster but in the same level of accuracy compared to the conventional MoM-PO method. The formulation is a bit different from (5) which is shown as(7)ZMoMIMoMi+1=VMoM+AMoM-POIPO+AMoM-POBPO-MoMIMoMi,where IMoMi is the current weights in the MoM region in the ith iteration process. An iterative process will be needed to derive the solution of (7) until the currents in the MoM region become stable. The error ξi in the ith step is defined as(8)ξi=IMoMi-IMoMi-1IMoMi-1,where the norm is defined in the 2-norm form and IMoMi stands for the current coefficients in the ith iteration. The threshold value of the relative residual error ξT and the maximum number of iteration Niter are also needed to control the iteration computation.The previous MLFMA-PO technique [12] considers the impedance matrix of the MoM by coupling the PO contributions while the iterative method updates the voltage matrix of the MoM by PO currents until it reaches stability. If only one whole octree is built, the computational complexity of the MVP processes of Z, A, and B is nearly in the same level since the MLFMA operations are all implemented based on the whole octree. Moreover, only one solving process is needed for the conventional MLFMA-PO technique while several solving processes are needed for the iterative MLFMA-PO method. Thus, the previous noniterative method may be a better choice if only one octree exists. However, the MVP process of Z based on the newly created octree structure is much faster than the other two if the extra octree for the MoM region is built. In the iterative MLFMA-PO method, several MVP processes of A and B are replaced by Z which makes the iterative method a more efficient solver. Then, the iterative MLFMA-PO method is suggested if two octrees are built which is also demonstrated in the numerical results.An explanation is also given in the mathematical manner. Assume that the computation time of the MVP scanning the whole octree isα while the value is β for the octree of the MoM region. Apparently, β is much smaller than the α. The number of the solving iterations for the conventional MLFMA-PO is defined as n1 while the average value is n2 (n2≈n1) for the iterative method with m (m is always smaller than 10) iterations between the MoM region and PO region. Thus, the computation time χmlfma-po of the conventional MLFMA-PO method is(9)χmlfma-po=3n1α.Since the value ofβ is negligible compared to α, the computation time for the iterative MLFMA-PO method can be defined as(10)χiter-mlfma-po=mn2β+2mα≈mn1β+2α≈2mα.From the above two equations, the computation time of the iterative MLFMA-PO method based on two octrees can be estimated as(11)χiter-mlfma-po≈2m3n1χmlfma-po.Sincem is always much smaller than n1, the computation time of the iterative MLFMA-PO based on two octrees will be much smaller than the previous MLFMA-PO method [12] using the same hardware.


## 2.3. IPO Process

In the real-world problems, the structures are always complicated and nonsmooth. Thus, the currents derived by the PO method may not be accurate enough. This could make the precision of the total solution limited. In particular for the structures with strong reflections in the PO region, the conventional MLFMA-PO technique will produce totally wrong solutions in some angles. This problem can be alleviated by introducing the iterative process for the PO region (IPO) to improve the accuracy of the PO currents.The basic MoM-IPO principle can be sought in [16, 17] and it is not presented here for simplicity. However, the modified IPO method proposed in this paper is a little different from the conventional IPO. Firstly, the RWG basis functions are used to expand the electric currents to maintain the continuity. Secondly, MLFMA is also adopted to decrease the computational complexity of the mutual interaction in PO region to O(Nlog⁡(N)). The proposed IPO process can be defined as(12)IPOi+1=CPOIPOi+IPO0,where IPOi stands for the ith currents in the PO region, IPO0 represents the initial currents produced by the incident wave and the currents in the MoM region, and the matrix CPO is the coupling matrix in the PO region. The elements in the CPO can be expressed by(13)CPO,ij=0,i=jδ(t^i++t^i-)·n^i×∇∫Tj±‍Gr,r′fjr′dτ′,i≠j,where t^i+ and t^i- are the two unit vectors which are perpendicular to the ith edge lying in the plane of the triangular pairs Ti± in the PO region. δ denotes the shadowing effects coefficient and Tj± is the triangular pairs of the jth RWG basis function. Note that the diagonal elements of the coupling matrix are zero since the self-interactions are not considered. In this paper, the MVP process of CPOIPOi in each iteration is accelerated by MLFMA technique to maintain the computational complexity to be O(Nlog⁡(N)). Compared to the conventional MLFMA for the PEC objects, implementation for IPO is straightforward except for the calculation of near-field impedance matrix and lowest-level receiving patterns. The near-field impedance matrix element can be computed based on (13). However, the lowest-level receiving patterns need some corrections compared to the conventional MLFMA which is shown as(14)Rfimk^=-δk^×t^i+×n^i++t^i-×n^i-e-jk·rim.What is more, the extra octree for the PO region is not suggested to build because the PO region always plays a dominant role in the whole structure. The complexity of MVP based on the octree for the PO region will not be decreased apparently compared to the whole one.


## 3. Numerical Results

In this section, three canonical examples are examined to demonstrate the high efficiency and the reasonable accuracy of the enhanced MLFMA-PO hybrid technique compared to the previous MLFMA-PO method and conventional MLFMA. For all the experiments, the mesh sizes are about 0.1λ (λ represents the wave length). The biconjugate gradient stabilized method (BiCGSTAB) equipped with the incomplete LU factorization with dual-threshold preconditioner (ILUT) is adopted as the iterative solver for the matrix equations in MoM region and the threshold of the relative residual error is chosen as 0.001. All the computations were carried out on a workstation with a 3.4 GHz Core i7 4930 k CPU and 64 GB of RAM with OpenMP technique and the digits were stored in double precision.


### 3.1. A Monopole Mounted on a Large Plane

Firstly, the efficiency and accuracy of the proposed iterative technique based on two octrees without iterations in the PO region are examined. A monopole mounted on an electrically large plane is the first example, and the model is displayed in Figure2. The work frequency of the monopole is 600 MHz. The monopole and a part of the plane which is near to the monopole are divided into the MoM region and the remaining part is given to the PO region. The numbers of unknowns for the MoM region and PO region are 3779 and 117110, respectively.Figure 2
A monopole mounted on an electrically large plane.The far-field gain patterns in the XOZ cut-plane obtained from conventional MLFMA-PO technique, iterative MLFMA-PO technique, and MLFMA are shown in Figure3. It can be seen that the three results agree well with each other. The results demonstrate that both of the MLFMA-PO techniques can achieve high accuracy compared to the low-frequency numerical methods. Moreover, memory usage and the computation time are presented in Table 1. It is obvious that the consumed memory and calculation time of the MLFMA-PO techniques are much less than the conventional MLFMA. What is more, it is shown that the computation time can be further decreased via the extra octree built for the MoM region while the memory usage only increases from 1.46 GB to 1.50 GB. The noniterative MLFMA-PO method is clearly the most efficient if only one whole octree is constructed, but the iterative MLFMA-PO method will be a better choice if the extra octree is built. Compared to the conventional MLFMA, the iterative MLFMA-PO method equipped with two octrees can achieve a speedup ratio of 10.69.Table 1
Calculation time and memory requirement of different methods for the monopole mounted on a plane.
Methods
Memory  (GB)
Calc. time  (seconds)
MLFMA
6.38
435.2
MLFMA-PO (1 tree)
1.46
103.1
MLFMA-PO (2 trees)
1.50
82.5
Iter-MLFMA-PO (1 tree)
1.46
171.9
Iter-MLFMA-PO (2 trees)
1.50
40.7Figure 3
The far-field gain pattern in XOZ plane of the monopole mounted on the plane.


### 3.2. Wire Antenna Mounted on a Ship Model

Secondly, a wire antenna mounted on a ship model is studied. The ship model is displayed in Figure4. The work frequency of the monopole is 300 MHz. In the hybrid analysis, the monopole, the cube, and the part of the ship near to the cube are assigned to the MoM region. The remaining part is placed into the PO region. The numbers of unknowns in the MoM and PO regions are 10621 and 398242, respectively.Figure 4
Wire antenna mounted on a ship model.The far-field gain patterns in the XOZ cut-plane and YOZ cut-plane obtained from conventional MLFMA-PO technique, iterative MLFMA-PO technique, and MLFMA are shown in Figures5 and 6. We can see that the three curves agree well with each other in the figures and the results derived by the two hybrid schemes are almost the same. Moreover, memory usage and the computation time are presented in Table 2. It is also demonstrated that the consumed memory and calculation time of the hybrid techniques are much less than the conventional MLFMA. The memory usage increases from 6.25 GB to 6.31 GB for the extra octree structure which can be negligible. Compared to the conventional MLFMA, the iterative MLFMA-PO method equipped with two octrees can be nearly 100 times faster than MLFMA while the speedup ratio is only 9.2 for the conventional hybrid method.Table 2
Calculation time and memory requirement of different methods for the monopole mounted on a plane.
Methods
Memory  (GB)
Calc. time  (seconds)
MLFMA
35.68
14772.6
MLFMA-PO (1 tree)
6.25
1605.7
MLFMA-PO (2 trees)
6.31
941.4
Iter-MLFMA-PO (1 tree)
6.25
2376.5
Iter-MLFMA-PO (2 trees)
6.31
155.7Figure 5
The far-field gain pattern in XOZ plane of the wire antenna mounted on ship model.Figure 6
The far-field gain pattern in YOZ plane of the wire antenna mounted on ship model.


### 3.3. A Monopole Installed on Folding Planes

Thirdly, the accuracy of the proposed technique with IPO which is called MLFMA-IPO is studied. A monopole installed on a large two-plates body is analyzed which is shown in Figure7. The work frequency of the monopole is also 600 MHz. The monopole and a small part of one plane which is near to the monopole are divided into the MoM region and the remaining part is given to the PO region. The numbers of unknowns for the MoM region and PO region are 7052 and 233783, respectively. The maximum iteration number and relative error threshold for IPO are set to be 3 and 0.1, respectively.Figure 7
A monopole installed on folding planes.In Figure8, the far-field gain patterns in the XOZ cut-plane obtained from MLFMA-IPO technique, MLFMA-PO technique, and MLFMA are plotted. It demonstrates that the results derived from the MLFMA-IPO agree well with results derived from the MLFMA while the ones derived from the conventional MLFMA-PO are not accurate in some angles (especially around 180 degrees). In Table 3, the number of unknowns, memory usage, and the computation time are presented. It is shown that the consumed memory and calculation time of the MLFMA-IPO and MLFMA-PO are much less than the conventional MLFMA. The proposed MLFMA-IPO method uses a little more memory compared to the MLFMA-PO. However, the MLFMA-IPO method is still much faster than the conventional MLFMA-PO method owing to the extra octree structure and the iteration scheme.Table 3
Calculation time and memory requirement of the two different methods for the monopole mounted on the folding planes.
Methods
Unknowns
Memory (GB)
Calc. time (s)
MoM region
PO region
MLFMA
240835
/
18.58
2142.5
MLFMA-PO
7052
233783
2.88
342.6
MLFMA-IPO
7052
233783
4.29
141.7Figure 8
The far-field gain pattern in XOZ plane of the monopole installed on the folding planes.


## 3.1. A Monopole Mounted on a Large Plane

Firstly, the efficiency and accuracy of the proposed iterative technique based on two octrees without iterations in the PO region are examined. A monopole mounted on an electrically large plane is the first example, and the model is displayed in Figure2. The work frequency of the monopole is 600 MHz. The monopole and a part of the plane which is near to the monopole are divided into the MoM region and the remaining part is given to the PO region. The numbers of unknowns for the MoM region and PO region are 3779 and 117110, respectively.Figure 2
A monopole mounted on an electrically large plane.The far-field gain patterns in the XOZ cut-plane obtained from conventional MLFMA-PO technique, iterative MLFMA-PO technique, and MLFMA are shown in Figure3. It can be seen that the three results agree well with each other. The results demonstrate that both of the MLFMA-PO techniques can achieve high accuracy compared to the low-frequency numerical methods. Moreover, memory usage and the computation time are presented in Table 1. It is obvious that the consumed memory and calculation time of the MLFMA-PO techniques are much less than the conventional MLFMA. What is more, it is shown that the computation time can be further decreased via the extra octree built for the MoM region while the memory usage only increases from 1.46 GB to 1.50 GB. The noniterative MLFMA-PO method is clearly the most efficient if only one whole octree is constructed, but the iterative MLFMA-PO method will be a better choice if the extra octree is built. Compared to the conventional MLFMA, the iterative MLFMA-PO method equipped with two octrees can achieve a speedup ratio of 10.69.Table 1
Calculation time and memory requirement of different methods for the monopole mounted on a plane.
Methods
Memory  (GB)
Calc. time  (seconds)
MLFMA
6.38
435.2
MLFMA-PO (1 tree)
1.46
103.1
MLFMA-PO (2 trees)
1.50
82.5
Iter-MLFMA-PO (1 tree)
1.46
171.9
Iter-MLFMA-PO (2 trees)
1.50
40.7Figure 3
The far-field gain pattern in XOZ plane of the monopole mounted on the plane.


## 3.2. Wire Antenna Mounted on a Ship Model

Secondly, a wire antenna mounted on a ship model is studied. The ship model is displayed in Figure4. The work frequency of the monopole is 300 MHz. In the hybrid analysis, the monopole, the cube, and the part of the ship near to the cube are assigned to the MoM region. The remaining part is placed into the PO region. The numbers of unknowns in the MoM and PO regions are 10621 and 398242, respectively.Figure 4
Wire antenna mounted on a ship model.The far-field gain patterns in the XOZ cut-plane and YOZ cut-plane obtained from conventional MLFMA-PO technique, iterative MLFMA-PO technique, and MLFMA are shown in Figures5 and 6. We can see that the three curves agree well with each other in the figures and the results derived by the two hybrid schemes are almost the same. Moreover, memory usage and the computation time are presented in Table 2. It is also demonstrated that the consumed memory and calculation time of the hybrid techniques are much less than the conventional MLFMA. The memory usage increases from 6.25 GB to 6.31 GB for the extra octree structure which can be negligible. Compared to the conventional MLFMA, the iterative MLFMA-PO method equipped with two octrees can be nearly 100 times faster than MLFMA while the speedup ratio is only 9.2 for the conventional hybrid method.Table 2
Calculation time and memory requirement of different methods for the monopole mounted on a plane.
Methods
Memory  (GB)
Calc. time  (seconds)
MLFMA
35.68
14772.6
MLFMA-PO (1 tree)
6.25
1605.7
MLFMA-PO (2 trees)
6.31
941.4
Iter-MLFMA-PO (1 tree)
6.25
2376.5
Iter-MLFMA-PO (2 trees)
6.31
155.7Figure 5
The far-field gain pattern in XOZ plane of the wire antenna mounted on ship model.Figure 6
The far-field gain pattern in YOZ plane of the wire antenna mounted on ship model.


## 3.3. A Monopole Installed on Folding Planes

Thirdly, the accuracy of the proposed technique with IPO which is called MLFMA-IPO is studied. A monopole installed on a large two-plates body is analyzed which is shown in Figure7. The work frequency of the monopole is also 600 MHz. The monopole and a small part of one plane which is near to the monopole are divided into the MoM region and the remaining part is given to the PO region. The numbers of unknowns for the MoM region and PO region are 7052 and 233783, respectively. The maximum iteration number and relative error threshold for IPO are set to be 3 and 0.1, respectively.Figure 7
A monopole installed on folding planes.In Figure8, the far-field gain patterns in the XOZ cut-plane obtained from MLFMA-IPO technique, MLFMA-PO technique, and MLFMA are plotted. It demonstrates that the results derived from the MLFMA-IPO agree well with results derived from the MLFMA while the ones derived from the conventional MLFMA-PO are not accurate in some angles (especially around 180 degrees). In Table 3, the number of unknowns, memory usage, and the computation time are presented. It is shown that the consumed memory and calculation time of the MLFMA-IPO and MLFMA-PO are much less than the conventional MLFMA. The proposed MLFMA-IPO method uses a little more memory compared to the MLFMA-PO. However, the MLFMA-IPO method is still much faster than the conventional MLFMA-PO method owing to the extra octree structure and the iteration scheme.Table 3
Calculation time and memory requirement of the two different methods for the monopole mounted on the folding planes.
Methods
Unknowns
Memory (GB)
Calc. time (s)
MoM region
PO region
MLFMA
240835
/
18.58
2142.5
MLFMA-PO
7052
233783
2.88
342.6
MLFMA-IPO
7052
233783
4.29
141.7Figure 8
The far-field gain pattern in XOZ plane of the monopole installed on the folding planes.


## 4. Conclusion

Some modifications are developed based on the previous proposed MLFMA-PO technique to improve the accuracy and efficiency of the analysis of the electrically large objects. MLFMA is adopted to accelerate all the existing MVP processes in the newly developed hybrid technique to ensure the low computational complexity. An extra octree for the MoM region and alternative iterative MLFMA-PO are also implemented and analyzed. What is more, a modified IPO method is also implemented in the PO region to improve the accuracy of the structures with strong reflections. Numerical results demonstrate the improved accuracy and efficiency of the proposed method compared to the conventional methods.



---

*Source: 826217-2015-03-23.xml*