Parallelizing a hybrid finite element-boundary integral method for the analysis of scattering and radiation of electromagnetic waves

This paper presents the practical experience of parallelizing a simulator of general scattering and radiation electromagnetic problems. The simulator stems from an existing sequential simulator in the frequency domain, which is based on a finite element analysis. After the analysis of a test case, two steps were carried out: first, a ''hand-crafted'' code parallelization of a convolution-type operation was developed within the kernel of the simulator. Second, the sequential HSL library, used in the existing simulator, was replaced by the parallel MUMPS (MUltifrontal Massively Parallel sparse direct Solver) library in order to solve the associated linear algebra problem in parallel. Such a library allows for the distribution of the factorized matrix and some of the computational load among the available processors. A test problem and three realistic (in terms of the number of unknowns) cases have been run using the parallelized version of the code, and the results are presented and discussed focusing on the memory usage and achieved speed-up.

[1]  Patrick Amestoy,et al.  Hybrid scheduling for the parallel solution of linear systems , 2006, Parallel Comput..

[2]  Tapan K. Sarkar,et al.  Fully coupled hybrid‐method FEM/high‐frequency technique for the analysis of 3D scattering and radiation problems , 2005 .

[3]  M. Salazar-Palma,et al.  Fully Coupled Multi-Hybrid FEM-PO/PTD-UTD Method for the Analysis of Scattering and Radiation Problems , 2006, 2006 12th Biennial IEEE Conference on Electromagnetic Field Computation.

[4]  M. Salazar-Palma,et al.  A finite element method for the analysis of radiation and scattering of electromagnetic waves on complex environments , 2005 .

[5]  B. McDonald,et al.  Finite-Element Solution of Unbounded Field Problems , 1972 .

[6]  Tapan K. Sarkar,et al.  Iterative and Self-Adaptive Finite-Elements in Electromagnetic Modeling , 1998 .

[7]  J. Pasciak,et al.  Computer solution of large sparse positive definite systems , 1982 .

[8]  O. Zienkiewicz,et al.  The coupling of the finite element method and boundary solution procedures , 1977 .

[9]  A. Peirce.  Computer Methods in Applied Mechanics and Engineering , 2010 .

[10]  Joseph W. H. Liu,et al.  The Multifrontal Method for Sparse Matrix Solution: Theory and Practice , 1992, SIAM Rev..

[11]  John B. Shoven,et al.  I , Edinburgh Medical and Surgical Journal.

[12]  Jianming Jin,et al.  Application of hybrid finite element method to electromagnetic scattering from coated cylinders , 1988 .

[13]  Bruno Després,et al.  A Domain Decomposition Method for the Helmholtz equation and related Optimal Control Problems , 1996 .

[14]  G. Borzì,et al.  Iteratively-improved Robin boundary conditions for the finite element solution of scattering problems in unbounded domains , 1998 .

[15]  C. Balanis.  Advanced Engineering Electromagnetics , 1989 .

[16]  Jürgen Schulze.  Towards a Tighter Coupling of Bottom-Up and Top-Down Sparse Matrix Ordering Methods , 2001 .

[17]  S. Sloan.  An algorithm for profile and wavefront reduction of sparse matrices , 1986 .

[18]  R. Fernandez-Recio,et al.  Fully Coupled Multi-Hybrid FEM-PO/PTD-UTD Method for the Analysis of Radiation Problems , 2007, IEEE Transactions on Magnetics.

[19]  R. Fernandez-Recio,et al.  Fully Coupled Hybrid FEM-UTD Method Using NURBS for the Analysis of Radiation Problems , 2008, IEEE Transactions on Antennas and Propagation.

[20]  Jian-Ming Jin,et al.  On the variational formulation of hybrid finite element-boundary integral techniques for electromagnetic analysis , 2004, IEEE Transactions on Antennas and Propagation.

[21]  Magdalena Salazar-Palma,et al.  Iterative and self-adaptive finite-elements in electromagnetic modeling , 1998 .

[22]  Jin-Fa Lee,et al.  A domain decomposition approach for non-conformal couplings between finite and boundary elements for unbounded electromagnetic problems in R3 , 2007, J. Comput. Phys..

[23]  M. Salazar-Palma,et al.  Second-order Ndlec tetrahedral element for computational electromagnetics , 2000 .

[24]  Patrick Amestoy,et al.  A Fully Asynchronous Multifrontal Solver Using Distributed Dynamic Scheduling , 2001, SIAM J. Matrix Anal. Appl..

[25]  Vipin Kumar,et al.  Analysis of Multilevel Graph Partitioning , 1995, Proceedings of the IEEE/ACM SC95 Conference.

[26]  Jennifer A. Scott.  On ordering elements for a frontal solver , 1999 .

[27]  B. Stupfel,et al.  A domain decomposition method for the vector wave equation , 2000 .

[28]  Patrick R. Amestoy,et al.  Multifrontal parallel distributed symmetric and unsymmetric solvers , 2000 .

[29]  George Karypis,et al.  Introduction to Parallel Computing , 1994 .