A Novel Broadband Multilevel Fast Multipole Algorithm With Incomplete-Leaf Tree Structures for Multiscale Electromagnetic Problems

An efficient and versatile broadband multilevel fast multipole algorithm (MLFMA), which is capable of handling large multiscale electromagnetic problems with a wide dynamic range of mesh sizes, is presented. By invoking a novel concept of incomplete-leaf tree structures, where only the overcrowded boxes are divided into smaller ones for a given population threshold, versatility of using variable-sized boxes is achieved. Consequently, for geometries containing highly overmeshed local regions, the proposed method is always more efficient than the conventional MLFMA for the same accuracy, while it is always more accurate if the efficiency is comparable. Furthermore, in such a population-based clustering scenario, the error is controllable regardless of the number of levels. Several canonical examples are provided to demonstrate the superior efficiency and accuracy of the proposed algorithm in comparison with the conventional MLFMA

Error Control of Multiple-Precision MLFMA

We introduce and demonstrate a new error control scheme for the computation of far-zone interactions in the multilevel fast multipole algorithm when implemented within a multiple-precision arithmetic framework. The proposed scheme provides the optimum truncation numbers as well as the machine precisions given the desired relative error thresholds and the box sizes for the translation operator at all frequencies. In other words, unlike the previous error control schemes which are valid only for high-frequency problems. the proposed scheme can be used to control the error across both low- and high-frequency problems. Optimum truncation numbers and machine precisions are calculated for a wide range of box sizes and desired relative error thresholds with the proposal error control scheme. The results are compared with the previously available methods and numerical surveys

Error Analysis of MLFMA with Closed-Form Expressions

IEEEThe current state-of-the-art error control of Multilevel Fast Multipole Algorithm (MLFMA) is valid for any given error threshold at any frequency, but it requires a multiple-precision arithmetic framework to be implemented. In this work, we use asymptotic approximations and curve-fitting techniques to derive accurate closed-form expressions for the error control of MLFMA that can be implemented in common fixed-precision computers. Moreover, using the proposed closed-form expressions in conjunction with the state-of-the-art scheme, we report novel design curves for MLFMA that can be used to determine achievable error limits, as well as the minimum box sizes that can be solved with a given desired error threshold for a wide range of machine precision levels

Broadband Multilevel Fast Multipole Algorithm For Large-Scale Problems With Nonuniform Discretizations

We present a broadband implementation of the multilevel fast multipole algorithm (MLFMA) for fast and accurate solutions of multiscale problems involving highly nonuniform discretizations. Incomplete tree structures, which are based on population-based clustering with flexible leaf-level boxes at different levels, are used to handle extremely varying triangulation sizes on the same structures. Superior efficiency and accuracy of the developed implementation, in comparison to the standard and broadband MLFMA solvers employing conventional tree structures, are demonstrated on practical problems

Broadband Analysis of Multiscale Electromagnetic Problems: Novel Incomplete-Leaf MLFMA for Potential Integral Equations

Recently introduced incomplete tree structures for the magnetic-field integral equation are modified and used in conjunction with the mixed-form multilevel fast multipole algorithm (MLFMA) to employ a novel broadband incomplete-leaf MLFMA (IL-MLFMA) to the solution of potential integral equations (PIEs) for scattering/radiation from multiscale open and closed surfaces. This population-based algorithm deploys a nonuniform clustering that enables to use deep levels safely and, when necessary, without compromising the accuracy resulting in an improved efficiency and a significant reduction for the memory requirements (order of magnitudes), while the error is controllable. The superiority of the algorithm is demonstrated in several canonical and real-life multiscale geometries

Efficient analysis of phased arrays of microstrip patches using a hybrid generalized forward backward method/Green's function technique with a DFT based acceleration algorithm

A hybrid method based on the combination of generalized forward backward method (GFBM) and Green's function for the grounded dielectric slab together with the acceleration of the combination via a discrete Fourier transform (DFT) based algorithm is developed for the efficient and accurate analysis of electromagnetic radiation/scattering from electrically large, irregularly contoured two-dimensional arrays consisting of finite number of probe-fed microstrip patches. In this method, unknown current coefficients corresponding to a single patch are first solved by a conventional Galerkin type hybrid method of moments (MoM)/Green's function technique that uses the grounded dielectric slab's Green's function. Because the current distribution on the microstrip patch can be expanded using an arbitrary number of subsectional basis functions, the patch can have any shape. The solution for the array currents is then found through GFBM, where it sweeps the current computation element by element. The computational complexity of this method, which is originally O (N-tot(2)) (N-tot being the total number of unknowns) for each iteration, is reduced to O(N-tot) using a DFT based acceleration algorithm making use of the fact that array elements are identical and the array is periodic. Numerical results in the form of array current distribution are given for various sized arrays of probe-fed microstrip patches with elliptical and/or circular boundaries, and are compared with the conventional MoM results to illustrate the efficiency and accuracy of the method

Wireless Monitoring of a Structural Beam To Be Used for Post-Earthquake Damage Assessment

Wireless monitoring of a standard reinforced concrete beam is shown in a simply supported beam experiment. The passive nested split-ring resonator (NSRR) probes are attached on the reinforcing bars (rebars) within the beam, and an antenna interrogates the probes from outside the beam. The results of the experiment show that the plastic deformation region strain/displacement can be detected by the wireless sensing system. The data collected by the system constitutes an important input for the assessment of the damage that can be observed after earthquakes