Analysis of Electromagnetic Scattering From Homogeneous Penetrable Objects by a Strong Skeletonization-Based Fast Direct Solver

被引:5
作者
Jiang, Ming [1 ]
Rong, Zhi [1 ]
Yang, Xiong [1 ]
Lei, Lin [1 ]
Li, Pei [2 ]
Chen, Yongpin [1 ]
Hu, Jun [1 ]
机构
[1] Univ Elect Sci & Technol China, Sch Elect Sci & Engn, Chengdu 611731, Peoples R China
[2] Sci & Technol Elect Informat Control Lab, Chengdu 610036, Peoples R China
基金
中国国家自然科学基金;
关键词
Magnetic separation; Surface impedance; Dielectrics; Skeleton; Current; Matrix decomposition; Integral equations; Electromagnetic (EM) scattering; fast direct solver; recursive skeletonization factorization; strong admissibility; HIERARCHICAL LU FACTORIZATION; INTEGRAL-EQUATION SOLVER; FAST MULTIPOLE ALGORITHM; MATRIX DECOMPOSITION; LINEAR-SYSTEMS; COMPRESSION; COMPUTATIONS; ACCURACY; MOMENTS;
D O I
10.1109/TAP.2022.3161474
中图分类号
TM [电工技术]; TN [电子技术、通信技术];
学科分类号
0808 ; 0809 ;
摘要
This article presents a fast direct solver based on surface integral equation (SIE) to solve electromagnetic (EM) scattering from homogeneous penetrable objects. The proposed method relies on a strong admissibility skeletonization factorization (SASF) algorithm and Poggio-Miller-Chang-Harrington-Wu-Tsai (PMCHWT) formulation. In the SASF scheme, only well-separated groups that satisfy the strong admissibility condition are considered to be compressed and thus relatively fewer skeleton basis functions are selected. It is an effective way to compress the matrix with small ranks. An independent compression technique is developed for dielectric problems involving electric and magnetic currents, in which skeleton basis functions representing electric and magnetic currents are obtained separately. Moreover, a novel strategy based on matrix normalization is proposed to treat the arising "fill-in" blocks when far-field interactions are compressed. Ultimately, the impedance matrix can be cast into products of a series of block unit triangular matrices and a block diagonal matrix. Several numerical results show that the proposed approach is effective and stable as well as provides an accurate solution for scattering problems of homogeneous penetrable objects.
引用
收藏
页码:6883 / 6892
页数:10
相关论文
共 42 条
[1]  
Adams R. J., 2011, IEEE T ANTENN PROPAG, V59, P2404
[2]  
Ambikasaran S, 2013, J SCI COMPUT, V57, P477, DOI 10.1007/s10915-013-9714-z
[3]   AIM: Adaptive integral method for solving large-scale electromagnetic scattering and radiation problems [J].
Bleszynski, E ;
Bleszynski, M ;
Jaroszewicz, T .
RADIO SCIENCE, 1996, 31 (05) :1225-1251
[4]  
Borm S., 2004, Computing and Visualization in Science, V7, P173, DOI 10.1007/s00791-004-0135-2
[5]   Fast Direct Solver for Essentially Convex Scatterers Using Multilevel Non-Uniform Grids [J].
Brick, Yaniv ;
Lomakin, Vitaliy ;
Boag, Amir .
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, 2014, 62 (08) :4314-4324
[6]   An H2-Matrix-Based Integral-Equation Solver of Reduced Complexity and Controlled Accuracy for Solving Electrodynamic Problems [J].
Chai, Wenwen ;
Jiao, Dan .
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, 2009, 57 (10) :3147-3159
[7]   On the compression of low rank matrices [J].
Cheng, H ;
Gimbutas, Z ;
Martinsson, PG ;
Rokhlin, V .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2005, 26 (04) :1389-1404
[8]   THE INVERSE FAST MULTIPOLE METHOD: USING A FAST APPROXIMATE DIRECT SOLVER AS A PRECONDITIONER FOR DENSE LINEAR SYSTEMS [J].
Coulier, Pieter ;
Pouransari, Hadi ;
Darve, Eric .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2017, 39 (03) :A761-A796
[9]   Surface-Volume-Surface EFIE Formulation for Fast Direct Solution of Scattering Problems on General 3-D Composite Metal-Dielectric Objects [J].
Gholami, Reza ;
Okhmatovski, Vladimir .
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, 2020, 68 (07) :5742-5747
[10]  
Greengard L, 2009, ACTA NUMER, V18, P243, DOI 10.1017/S0962492906410011