For efficiently solving large dense complex linear systems that arise in the electric field integral equation (EFIE) formulation of electromagnetic wave scattering problems, the multilevel fast multipole algorithm (ML...
详细信息
For efficiently solving large dense complex linear systems that arise in the electric field integral equation (EFIE) formulation of electromagnetic wave scattering problems, the multilevel fast multipole algorithm (MLFMA) is used to speed up the matrix vector product operations, and the sparse approximate inverse (SAI) preconditioning technique is employed to accelerate the convergence rate of the generalized minimal residual (GMRES) iterative method. We show that the convergence rate can be greatly improved by augmenting to the GMRES method a few eigenvectors associated with the smallest eigenvalues of the preconditioned system. Numerical experiments indicate that this new variant GMRES method is very effective with the MLFMA and can reduce both the iteration number and the computational time significantly. (C) 2008 Wiley Periodicals, Inc.
A novel kind of Buffa-Christiansen (BC) functions is developed for improving the surface combined field integral equation with impedance boundary condition (CFIE-IBC) from three-dimensional imperfectly conducting obje...
详细信息
A novel kind of Buffa-Christiansen (BC) functions is developed for improving the surface combined field integral equation with impedance boundary condition (CFIE-IBC) from three-dimensional imperfectly conducting objects. With difference from the original BC function in wide use, the novel BC function is dimensionless as the RWG function, which makes the L and K components of the resulting matrix element from the CFIE-IBC dimensionally homogeneous. As a result, the eigenvalue condition number of the matrix is very low, which leads to a fast-convergence solution whether the surface impedance is low or high. The CFIE-IBC with the novel BC function is combined with the overlapped domain decomposition method for reducing the memory requirement. In its application to the multilevel fast multipole algorithm, the far-field interactions can be simplified for reducing the computational cost.
In this letter, we consider the solution of large electromagnetics problems of composite structures with coexisting open and closed conductors. By modifying the construction of the combined-field integral equation (CF...
详细信息
In this letter, we consider the solution of large electromagnetics problems of composite structures with coexisting open and closed conductors. By modifying the construction of the combined-field integral equation (CFIE), we demonstrate a significant improvement in the iterative solution of these problems compared to the conventional electric-field formulation.
The validity of the use of continuity condition (CC) is studied when it is combined with the volume-surface integral equation (VSIE) to solve the electromagnetic scattering/radiation of composite conductor-dielectric ...
详细信息
The validity of the use of continuity condition (CC) is studied when it is combined with the volume-surface integral equation (VSIE) to solve the electromagnetic scattering/radiation of composite conductor-dielectric objects. It is found that the explicit enforcement of CC in VSIE may be inaccurate if the involved conductor is infinitesimally thin, but it works well for the composite structures that only contain conductors with closed surfaces. Anew combined field integral equation-volume integral equation with enforced CC (CC-CFIE-VIE) is then proposed, and a convenient way to embed it into the context of the multilevel fast multipole algorithm is also provided. Numerical results show that CC-CFIE-VIE not only reduces the overall memory usage and the CPU time of a single iteration, but also decreases the number of iterations to reach desired accuracy.
In this letter, we consider iterative solutions of the three-dimensional electromagnetic scattering problems formulated by surface integral equations. We show that solutions of the electric-field integral equation (EF...
详细信息
In this letter, we consider iterative solutions of the three-dimensional electromagnetic scattering problems formulated by surface integral equations. We show that solutions of the electric-field integral equation (EFIE) can be improved by employing an iterative least-squares QR (LSQR) algorithm. Compared to many other Krylov subspace methods, LSQR provides faster convergence and it becomes an alternative choice to the time-efficient no-restart generalized minimal residual (GMRES) algorithm that requires large amounts of memory. Improvements obtained with the LSQR algorithm become significant for the solution of large-scale problems involving open surfaces that must be formulated using EFIE, which leads to matrix equations that are usually difficult to solve iteratively, even when the matrix-vector multiplications are accelerated via the multilevel fast multipole algorithm.
We present an approximate diagonalization of the Green's function that is stable at arbitrarily short distances with respect to wavelength. The diagonalization is based on scaled spherical functions and plane wave...
详细信息
We present an approximate diagonalization of the Green's function that is stable at arbitrarily short distances with respect to wavelength. The diagonalization is based on scaled spherical functions and plane waves, where a scaling factor is used to stabilize special functions with small arguments. Optimization of the scaling factor leads to accurate diagonalizations, which can be used to implement the multilevel fast multipole algorithm for low-frequency problems.
Although the computational complexity of the finite-element boundary-integral (FE-BI) method is significantly reduced by the multilevel fast multipole algorithm (MLFMA), this MLFMA-enhanced FE-BI solution experiences ...
详细信息
Although the computational complexity of the finite-element boundary-integral (FE-BI) method is significantly reduced by the multilevel fast multipole algorithm (MLFMA), this MLFMA-enhanced FE-BI solution experiences a very slow convergence for some complex inhomogeneous problems. In this paper, a hybrid algorithm, combining direct methods with iterative methods, is designed to speed up the rate of convergence of this MLFMA-enhanced FE-BI solution. This hybrid algorithm is efficiently implemented with the aid of a newly developed package, SuperLU, of the LU decomposition solver. Numerical experiments are performed for a coated Northrop wing to demonstrate the efficiency of this hybrid algorithm. More importantly, the thorough investigation of the numerical experiments clearly shows the better accuracy, stability, and robustness of this hybrid algorithm over the conventional algorithms.
In this letter, a novel single-source equivalence principle algorithm is presented. In the traditional equivalence principle algorithm (EPA), double sources, which are equivalence electric and magnetic currents, are c...
详细信息
In this letter, a novel single-source equivalence principle algorithm is presented. In the traditional equivalence principle algorithm (EPA), double sources, which are equivalence electric and magnetic currents, are constructed on the equivalence surface to replace the currents of the enclosed scatterers. According to the extinction theorem, electric current is related to magnetic current. Therefore, only electric or magnetic current, called single-source, is the final unknown on the equivalence surface. The accuracy of two kinds of single-source EPA is investigated. Multi-level fastmultipolealgorithm (MLFMA) is used to accelerate the couplings between equivalence surfaces. Numerical results demonstrate the accuracy and efficiency of the proposed method.
Accurate electromagnetic (EM) analysis for conductive media requires to consider the finite conductivity of the media. Although the problems can be formulated by surface integral equations (SIEs) with an approximate s...
详细信息
Accurate electromagnetic (EM) analysis for conductive media requires to consider the finite conductivity of the media. Although the problems can be formulated by surface integral equations (SIEs) with an approximate surface impedance, we treat the conductive media as penetrable objects and use volume integral equations (VIEs) to exactly describe their EM features. The VIEs are solved by a point-matching scheme which does not rely on any basis and testing functions and allows the use of nonconforming meshes. Since the VIEs are well-conditioned in general and the integral kernels are free of material parameters, the scheme is flexible to accommodate a wide range of skin depth. Moreover, when the skin depth is small, we can only discretize the skin domain with a current distribution to reduce the cost. For large conductive media, we also incorporate the scheme with the multilevel fast multipole algorithm (MLFMA) to accelerate the solution. Typical numerical examples are presented to illustrate the scheme and its effectiveness has been validated.
This paper reviews the recent progress of fast integral equation solvers at the Center for Computational Electromagnetics and Electromagnetics laboratory, University of Illinois at Urbana-Champaign. We will demonstrat...
详细信息
This paper reviews the recent progress of fast integral equation solvers at the Center for Computational Electromagnetics and Electromagnetics laboratory, University of Illinois at Urbana-Champaign. We will demonstrate the ability to solve a variety of electromagnetic problems for complex structures and low-frequency structures as well as large scale scattering problems with over 10 million unknowns. (C) 2003 Elsevier Ltd. All rights reserved.
暂无评论