Parallelized multilevel fast multipole algorithm for scattering by objects with anisotropic impedance surfaces

A parallelized multilevel fast multipole algorithm (MLFMA) is presented for simulating electromagnetic scattering from complex targets with anisotropic impedance surfaces. By employing both surface electric and magnetic currents as unknowns and weakly enforcing the anisotropic impedance boundary condition, a combined integral equation is formulated to generate a set of well‐conditioned linear systems to be solved by MLFMA. To further improve the iterative convergence of the linear systems, a parallel sparse approximate inverse preconditioner is constructed from the near‐field interaction of the system matrix. The MLFMA is parallelized to enable computation on a large number of processors for large‐scale problems. Several numerical examples are presented to validate the algorithm and demonstrate its accuracy, scalability, and capability in handling large complex objects with anisotropic impedance surfaces. Copyright © 2014 John Wiley & Sons, Ltd.     

SSOR preconditioned GMRES for the FEM analysis of waveguide discontinuities with anisotropic dielectric

The anisotropic media and active properties of the perfectly matched layer (PML) absorbers significantly deteriorate the finite element method (FEM) system condition and as a result, convergence of the iterative solver is substantially slowed down. To address this issue, the symmetric successive over‐relaxation (SSOR) preconditioning scheme is applied to the generalized minimal residual (GMRES) for solving a large sparse and non‐symmetric system of linear equations resulting from the analysis of ferrite waveguide device by use of edge‐based FEM. Consequently, this preconditioned GMRES (PGMRES) approach can reach convergence ten times faster than GMRES for the typical structures. Copyright © 2004 John Wiley & Sons, Ltd.     

Parallel FETI‐DP algorithm for efficient simulation of large‐scale EM problems

An efficient parallelization of the dual‐primal finite‐element tearing and interconnecting (FETI‐DP) algorithm is presented for large‐scale electromagnetic simulations. As a nonoverlapping domain decomposition method, the FETI‐DP algorithm formulates a global interface problem, whose iterative solution is accelerated with a solution of a global corner problem. To achieve a good load balance for parallel computation, the original computational domain is decomposed into subdomains with similar sizes and shapes. The subdomains are then distributed to processors based on their close proximity to minimize inter‐processor communication. The parallel generalized minimal residual method, enhanced with the iterative classical Gram‐Schmidt orthogonalization scheme to reduce global communication, is adopted to solve the global interface problem with a fast convergence rate. The global corner‐related coarse problem is solved iteratively with a parallel communication‐avoiding biconjugate gradient stabilized method to minimize global communication, and its convergence is accelerated by a diagonal preconditioner constructed from the coarse system matrix. To alleviate neighboring communication overhead, the non‐blocking communication approach is employed in both generalized minimal residual and communication‐avoiding biconjugate gradient stabilized iterative solutions. Three numerical examples are presented to demonstrate the accuracy, scalability, and capability of the proposed parallel FETI‐DP algorithm for electromagnetic modeling of general objects and antenna arrays. Copyright © 2016 John Wiley & Sons, Ltd.     

A non‐conformal FETI‐like domain decomposition approach of FE‐BI‐MLFMA for 3‐D electromagnetic scattering/radiation problems

A non‐conformal finite element tearing and interconnecting‐like (FETI‐like) domain decomposition approach (DDA) of the hybrid finite element–boundary integral–multilevel fast multipole algorithm (FE‐BI‐MLFMA) is presented by integrating a series efficient techniques for computing electromagnetic scattering/radiation problems. The Robin transmission condition is employed to cement the non‐conformal meshes on the interconnected surfaces between the interior and exterior regions and between sub‐domains in the interior region. The FETI‐like technique is applied to reduce the FE‐BI matrix equation. Furthermore, a preconditioner is constructed to accelerate the convergent speed of this non‐conformal FETI‐like DDA. The numerical performance of the presented non‐conformal FETI‐like DDA‐FE‐BI‐MLFMA is studied for scattering/radiation problems. Copyright © 2015 John Wiley & Sons, Ltd.     

The modified multilevel compressed block decomposition algorithms for analyzing the scattering of objects in half space

The convergence rate of iterative methods can vary in an unpredictable way. It is related to the matrix condition number, which is notoriously bad for the electric field integral equation in the large‐scale electromagnetic problems. Therefore, an efficient direct solution—a multilevel compressed block decomposition (MLCBD) algorithm based on the adaptive cross‐approximation algorithm—is applied to overcome this problem; it is very efficient for the monostatic problems. Simulation results of the objects up and below ground in half space demonstrate that the proposed MLCBD method is efficient for analyzing electromagnetic problems. Copyright © 2013 John Wiley & Sons, Ltd.     

On the development of efficient FDTD‐PML formulations for general dispersive media

A novel implementation of the perfectly matched layer (PML) absorbing boundary condition (ABC) to terminate the finite‐difference time‐domain (FDTD) algorithm for general dispersive and negative index materials is presented. The proposed formulation also adopts the complex frequency‐shifted (CFS) approach, involves simple FDTD expressions and avoids complex arithmetic. Several FDTD‐PML simulations with different parameters are conducted for the termination of various dispersive media validating the stability, accuracy and effectiveness of the schemes and indicating the advantage of the CFS‐PML. Copyright © 2008 John Wiley & Sons, Ltd.