Simulation of Challenging Electromagnetic Problems Using a Massively Parallel Finite Element Method Solver

Abstract

This communication presents an efficient massively parallel finite element method solver for the solution of complex and electrically large electromagnetic problems with arbitrary structures. The solver makes use of a domain decomposition algorithm to decompose the original problem into several non-overlapping sub-domains that may be solved independently in parallel through the application of the corresponding transmission conditions on the interfaces of the adjacent sub-domains. A numerical exact mesh truncation algorithm called finite element-iterative integral equation evaluation, accelerated with multilevel fast multipole algorithm, is implemented to meet the highly accurate requirements of today's challenging simulations. What's more, a hybrid message passing interface and an open multi-processing parallel framework are designed to achieve large-scale parallel performance on supercomputers. Through several numerical examples, the accuracy, effectiveness, and scalability of the proposed solver will be demonstrated, achieving more than 60% parallel efficiency on an eight times CPU core scale (from 1280 to 10 240 cores). INDEX TERMS Finite element method (FEM), domain decomposition method (DDM), large-scale parallel computing, finite element-iterative integral equation evaluation (FE-IIEE), ten thousand CPU cores. 20346

Remote inspection of adversary-controlled environments

“…EM simulations of our “room” would require meshes with size to guarantee stability and accuracy of the solutions (including near field effects close to the antennas) 9 . Even for the relatively small-sized room used for our experiment , the attacker would need to solve, within a few seconds, a problem involving trillion of unknowns, which is three orders of magnitude larger than current state-of-the-art large-scale EM simulations run on supercomputers 10 , 11 . While mesh optimization and domain decomposition techniques could be used to speed up calculations, they would be inherently limited by the complexity of the problem.…”

On the Efficiency of the Advanced TWA Approach to the 60-GHz Microstrip Antenna Analysis for 5G Wireless Communication Systems

“…The 2D-TWA can be accelerated by applying the anisotropic mesh technique [39] which calls the 2D-NUFFT to manage precisely the non-uniform discretization. Memory and computational complexity for both the Forward and Backward of the ATWA process is guaranteed at O(B CD;B ) time, where B represents the meshing density applied to the structure to be analyzed, whereas the overall computational complexity of iterative solvers is O(B * ), the same with MOMand FEM-based solvers [40]. Table I recapitulates the iterative process of the ATWA approach, highlighting at each step the computational effort and the total time efficiency.…”

A Rigorous Code Verification Process of the Domain Decomposition Method in a Finite Element Method for Electromagnetics

“…Garcia-Castillo are with the Department of Signal Theory and Communications, University Carlos III of Madrid, 28911, Leganes, Spain (e-mail: aamor@ing.uc3m.es, legcasti@ing.uc3m.es) O. Floch, L.L. Toth, and R. Dyczij-Edlinger are with the Lehrstuhl für Theoretische Elektrotechnik, Universität des Saarlandes, Saarbrücken [18], [19], [20], [21], [22], [23], [24], [25], [26], [27], [28], [29], [30], [31], [32], [33], [34], [35], [36], [37], [38], [39], [40], [41], [42]. All these references can be divided into three groups: i) optimized Schwarz methods, focused on the optimization of the so-called transmission conditions, [16], [25], [30]; ii) cement element methods, which use transmission conditions through the definition of cement variables with physical meaning, [13], [14], [17], [18], [19], [20], [21], [22]; iii) finite element tear and interconnecting (FETI) methods, which use Lagrange multipliers to obtain a reduced system of equations, [15], [23], [26], [27], [24], [32],…”