Improvement and performance of parallel multilevel fast multipole algorithm

The method of establishing data structures plays an important role in the efficiency of parallel multilevel fast multipole algorithm（PMLFMA）.Considering the main complements of multilevel fast multipole algorithm（MLFMA） memory,a new parallelization strategy and a modified data octree construction scheme are proposed to further reduce communication in order to improve parallel efficiency.For far interaction,a new scheme called dynamic memory allocation is developed.To analyze the workload balancing performance of a parallel implementation,the original concept of workload balancing factor is introduced and verified by numerical examples.Numerical results show that the above measures improve the parallel efficiency and are suitable for the analysis of electrical large-scale scattering objects.

Efficient analysis of dielectric radomes using multilevel fast multipole algorithm with CRWG basis

A full-wave analysis of the electromagnetic problem of a three-dimensional （3-D） antenna radiating through a 3-D dielectric radome is preserued. The problem is formulated using the Poggio-Miller-Chang-Harrington- Wu（PMCHW） approach for homogeneous dielectric objects and the electric field integral equation for conducting objects. The integral equations are discretized by the method of moment （MoM）, in which the conducting and dielectric surface/interfaces are represented by curvilinear triangular patches and the unknown equivalent electric and magnetic currents are expanded using curvilinear RWG basis functions. The resultant matrix equation is then solved by the multilevel fast multipole algorithm （MLFMA） and fast far-field approximation （FAFFA） is used to further accelerate the computation. The radiation patterns of dipole arrays in the presence of radomes are presented. The numerical results demonstrate the accuracy and versatility of this method.

General and efficient parallel approach of finite element-boundary integral-multilevel fast multipole algorithm

A general and efficient parallel approach is proposed for the first time to parallelize the hybrid finiteelement-boundary-integral-multi-level fast multipole algorithm （FE-BI-MLFMA）. Among many algorithms of FE-BI-MLFMA, the decomposition algorithm （DA） is chosen as a basis for the parallelization of FE-BI-MLFMA because of its distinct numerical characteristics suitable for parallelization. On the basis of the DA, the parallelization of FE-BI-MLFMA is carried out by employing the parallelized multi-frontal method for the matrix from the finiteelement method and the parallelized MLFMA for the matrix from the boundary integral method respectively. The programming and numerical experiments of the proposed parallel approach are carried out in the high perfor- mance computing platform CEMS-Liuhui. Numerical experiments demonstrate that FE-BI-MLFMA is efficiently parallelized and its computational capacity is greatly improved without losing accuracy, efficiency, and generality.

Aerodynamic/stealth design of S-duct inlet based on discrete adjoint method

It is a major challenge for the airframe-inlet design of modern combat aircrafts,as the flow and electromagnetic wave propagation in the inlet of stealth aircraft are very complex.In this study,an aerodynamic/stealth optimization design method for an S-duct inlet is proposed.The upwind scheme is introduced to the aerodynamic adjoint equation to resolve the shock wave and flow separation.The multilevel fast multipole algorithm(MLFMA)is utilized for the stealth adjoint equation.A dorsal S-duct inlet of flying wing layout is optimized to improve the aerodynamic and stealth characteristics.Both the aerodynamic and stealth characteristics of the inlet are effectively improved.Finally,the optimization results are analyzed,and it shows that the main contradiction between aerodynamic characteristics and stealth characteristics is the centerline and crosssectional area.The S-duct is smoothed,and the cross-sectional area is increased to improve the aerodynamic characteristics,while it is completely opposite for the stealth design.The radar cross section(RCS)is reduced by phase cancelation for low frequency conditions.The method is suitable for the aerodynamic/stealth design of the aircraft airframe-inlet system.

Fast 3D EM scattering and radiation solvers based on MLFMA

As the fastest integral equation solver to date, the multilevel fast multipole algorithm （MLFMA） has been applied successfully to solve electromagnetic scattering and radiation from 3D electrically large objects. But for very large-scale problems, the storage and CPU time required in MLFMA are still expensive. Fast 3D electromagnetic scattering and radiation solvers are introduced based on MLFMA. A brief review of MLFMA is first given. Then, four fast methods including higher-order MLFMA （HO-MLFMA）, fast far field approximation combined with adaptive ray propagation MLFMA （FAFFA-ARP-MLFMA）, local MLFMA and parallel MLFMA are introduced. Some typical numerical results demonstrate the efficiency of these fast methods.

Overview: Electromagnetic Scattering from Ocean Surface

Understanding the sea surface scattering process is very important in the development of models to detect the target above or under the surface. In this paper, both the analytical and the numerical methods applied in sea surface scattering are summarized. Some important problems concerned in this field are discussed. For numerical study, edge effect brings artificial nonrealistic scattering and therefore must be suppressed. Different edge treatment methods are compared in this paper. Scattering of breaking wave surface at very low grazing angle always needs more attentions than other scattering problems. Some numerical results show the existence of the special phenomena at very low grazing angle, for example, the ＂sea spikes＂ and the Doppler splitting.

Indoor field strength prediction based on FMM

In this article, the fast multipole method （FMM） is applied to analgze the field strength of an ultra-wideband （UWB） signal. Small coverage of UWB communication systems and high efficiency of the algorithm make it possible to calculate the amplitude of electric field accurately. A homogeneous dielectric body and a multilayered dielectric object are studied. The computational results obtained by applying the FMM agree well with those obtained by applying the method of moments （MoM）.