Based on the addition theorem, the principle of a multilevel ray-propagation fastmultipolealgorithm (RPFMA) and fast far-field approximation (FAFFA) has been demonstrated for three-dimensional (3-D) electromagnetic ...
详细信息
Based on the addition theorem, the principle of a multilevel ray-propagation fastmultipolealgorithm (RPFMA) and fast far-field approximation (FAFFA) has been demonstrated for three-dimensional (3-D) electromagnetic scattering problems. From a rigorous mathematical derivation, the relation among RPFMA, FAFFA, and a conventional multilevel fast multipole algorithm (MLFMA) has been clearly stated. For very large-scale problems, the translation between groups in the conventional MLFMA is expensive because the translator is defined on an Ewald sphere with many sampling (k) over cap directions. When two groups are well separated, the translation can be simplified using RPFMA, where only a few sampling (k) over cap directions are required within a cone zone on the Ewald sphere. When two groups are in the far-field region, the translation can be further simplified by using FAFFA where only a single (k) over cap is involved in the translator along the ray-propagation direction. Combining RPFMA and FAFFA with MLFMA, three algorithms RPFMA-MLFMA, FAFFA-MLFMA, and RPFMA-FAFFA-MLFMA have been developed, which are more efficient than the conventional MLFMA in 3-D electromagnetic scattering and radiation for very large structures. Numerical results are given to verify the efficiency of the algorithms.
In practical applications, frequency-selective surfaces (FSSs) are finite, and sometimes even curved. In this paper, we present a hybrid volume-surface integral-equation approach to analyze the transmission and reflec...
详细信息
In practical applications, frequency-selective surfaces (FSSs) are finite, and sometimes even curved. In this paper, we present a hybrid volume-surface integral-equation approach to analyze the transmission and reflection characteristics of finite and curved FFS structures. The hybrid integral equations are established using the surface- and volume-equivalent principles. This approach has two advantages. One is the capability of modeling arbitrarily shaped FSS structures in detail, the other one allows its to easily apply the multilevelfast multiple algorithm to speed up the solution process. The scattering characteristics and frequency responses of several FSSs are analyzed. The simulation results show that for a finite-sized FSS, reducing the radius of curvature causes amplitude variation, frequency shift, and bandwidth change in the reflection and transmission responses. (c) 2005 Wiley Periodicals, Inc.
The numerical solution of wave scattering from large objects or from a large cluster of scatterers requires excessive computational resources and it becomes necessary to use approximate-but fast-methods such as the fa...
详细信息
The numerical solution of wave scattering from large objects or from a large cluster of scatterers requires excessive computational resources and it becomes necessary to use approximate-but fast-methods such as the fastmultipole method;however, since these methods are only approximate, it is important to have an estimate for the error introduced in such calculations. An analysis of the error for the fastmultipole method is presented and estimates for truncation and numerical integration errors are obtained. The error caused by polynomial interpolation in a multilevel fast multipole algorithm is also analyzed. The total error introduced in a multilevel implementation is also investigated numerically.
作者:
Ohnuki, SChew, WCUniv Illinois
Dept Elect & Comp Engn Ctr Computat Electromagnet Urbana IL 61801 USA Univ Illinois
Dept Elect & Comp Engn Electromagnet Lab Urbana IL 61801 USA
The multilevel fast multipole algorithm is based on the multipole expansion, which has numerical error sources such as truncation of the addition theorem, numerical integration, and interpolation/anterpolation. Of the...
详细信息
The multilevel fast multipole algorithm is based on the multipole expansion, which has numerical error sources such as truncation of the addition theorem, numerical integration, and interpolation/anterpolation. Of these, we focus on the truncation error and discuss its control precisely. The conventional selection rule fails when the buffer size is small compared to the desired numerical accuracy. We propose a new approach and show that the truncation error can be controlled and predicted regardless of the number of buffer sizes.
This study presents design and optimization of compact and efficient nanooptical couplers involving photonic crystals. Nanooptical couplers that have single and double input ports are designed to obtain efficient tran...
详细信息
This study presents design and optimization of compact and efficient nanooptical couplers involving photonic crystals. Nanooptical couplers that have single and double input ports are designed to obtain efficient transmission of electromagnetic waves in desired directions. In addition, these nanooptical couplers are cascaded by adding one after another to realize electromagnetic transmission systems. In the design and optimization of all these nanooptical couplers, the multilevel fast multipole algorithm, which is an efficient full-wave solution method, is used to perform electromagnetic analyses and simulations. A heuristic optimization method based on genetic algorithms is employed to obtain effective designs that provide the highest efficiency values. Two types of optimization strategies are applied using nanorods with a fixed length and using nanorods with varying lengths. This way, photonic crystals consisting of irregular arrays of both identical and nonidentical dielectric elements are designed for the realization of nanooptical couplers. The designs and their numerical results show that it is possible to design and further improve efficient nanooptical couplers with simple and compact geometries based on the principles of photonic crystals. Using relatively simple geometries and a single material, the designed nanooptical couplers are more preferable than the available designs in the literature.
The translation matrix for the multilevel fast multipole algorithm (MLFMA) in an FISC (fast Illinois solver code) is calculated directly, and the complexity is O(N-3/2), where N is the number of unknowns. For a proble...
详细信息
The translation matrix for the multilevel fast multipole algorithm (MLFMA) in an FISC (fast Illinois solver code) is calculated directly, and the complexity is O(N-3/2), where N is the number of unknowns. For a problem with a small electrical size, the CPU time for calculating the translation matrix can be negligible. But for large problems, the calculation time increases significantly. In this paper, we use interpolation to calculate the translation matrix, and the complexity is reduced to O(N). Different interpolation techniques are tested, and it is found that the Lagrange polynomial interpolation with high sampling rates is the best. The saving factor is 10 for the VFY218 at 4 GHz. (C) 2001 John Wiley & Sons, Inc.
In this paper, an updated equivalence principle algorithm (EPA) is presented. Compared with previous work by the authors, the high-order point sampling scheme is used to reduce the high-frequency noise in field projec...
详细信息
In this paper, an updated equivalence principle algorithm (EPA) is presented. Compared with previous work by the authors, the high-order point sampling scheme is used to reduce the high-frequency noise in field projection. An updated tap basis scheme is introduced to simplify the formulation. With these schemes, the equivalence principle algorithm can be used to solve multiscale problems by substituting the oversampled region with a proper equivalence surface. Moreover, EPA can also be accelerated using attached unknown accelerations and multilevel fast multipole algorithm so that large multiscale problems can be solved efficiently.
We present computational analysis, optimization, and design of optical couplers that can be useful to improve the transmission along bended nanowires. After demonstrating the deteriorated energy transmission due to sh...
详细信息
We present computational analysis, optimization, and design of optical couplers that can be useful to improve the transmission along bended nanowires. After demonstrating the deteriorated energy transmission due to sharp bends, which lead to out-of-phase nanowires and diffraction, we use a rigorous simulation environment to design efficient couplers made of spherical particles. For this purpose, an optimization module based on genetic algorithms is combined with the multilevel fast multipole algorithm, leading to a full-wave environment for precise designs of couplers. Numerical examples involving silver nanowires are presented to demonstrate the effectiveness of the optimization mechanism.
We consider the preconditioned iterative solution of large dense linear systems, where the coefficient matrix is a complex valued matrix arising from discretizing the integral equation of electromagnetic scattering. F...
详细信息
We consider the preconditioned iterative solution of large dense linear systems, where the coefficient matrix is a complex valued matrix arising from discretizing the integral equation of electromagnetic scattering. For some scattering structures this matrix can be poorly conditioned. The main purpose of this study is to evaluate the efficiency of a class of incomplete LU (ILU) factorization preconditioners for solving this type of matrices. We solve the electromagnetic wave equations using the BiCG method with an ILU preconditioner in the context of a multilevel fast multipole algorithm (MLFMA). The novelty of this work is that the ILU preconditioner is constructed using the near part block diagonal submatrices generated from the MLFMA. Experimental results show that the ILU preconditioner reduces the number of BiCG iterations substantially, compared to the block diagonal preconditioner. The preconditioned iteration scheme also maintains the computational complexity of the MLFMA, and consequently reduces the total CPU time. (C) 2002 Elsevier Science B.V. All rights reserved.
In this article, numerical result obtained using two different commercially available electromagnetic (EM) codes based on an integral equation method are presented. One of the methods solves the integral equation by t...
详细信息
In this article, numerical result obtained using two different commercially available electromagnetic (EM) codes based on an integral equation method are presented. One of the methods solves the integral equation by the method of moments (MoM) using higher order basis functions (the commercially available parallel software named TIDES), and the other one solves the integral equation by MoM with a multilevel fast multipole algorithm (MLFMA) incorporating the RWG basis functions (the commercially available software FEKO). Numerical results indicate that the integral equation solver with higher order basis functions occupies less than 40% of the RAM and the execution little is significantly less when compared with that of using the RWG basis functions along with MLFMA. (C) 2009 Wiley Periodicals, Inc. Microwave Opt Technol Lett 51: 1219-1225, 2009;Published online in Wiley Inter-Science (***). DOI 10.1002/mop.24271
暂无评论