一种使用八叉树的半空间MLFMA区域分解算法

doi:10.19665/j.issn1001-2400.2021.06.018

摘要/Abstract

摘要：

为了在有限资源的情况下快速准确地分析半空间电大目标的电磁散射特性,提出一种使用八叉树结构的并行半空间多层快速多极子区域分解算法。通过利用多层快速多极子算法在空间中形成的八叉树结构,自适应地对未知量进行分组来实现区域的划分,从而避免了在区域之间创建虚拟交界面,降低了模型的处理工作量。为保证区域间电流的连续性,在边界处严格计算1/4阻抗,使结果更加精确。对于半空间环境,引入复镜像源来计算半空间近相互作用,并且引入实镜像源来计算半空间远相互作用。通过与商业软件FEKO整体解对比分析,证明了算法的可靠性和精确性。采用核外算法将转移因子等复用数据存储在硬盘中,并与整体解资源消耗进行对比,证明这种算法能有效地减少内存消耗。最后,给出一个千波长半空间舰船区域分解算例,证明所提出的算法能够在有限的计算资源情况下仿真电大尺寸模型。

关键词: 区域分解, 并行计算, 多层快速多极子, 半空间问题, 八叉树策略, 核外算法

Abstract:

In order to quickly and accurately analyze the electromagnetic scattering of electrical large objects in half space under the condition of limited resources,a parallel half-space multi level fast multipole algorithm (MLFMA) with the domain decomposition method utilizing an octree is proposed.By using the octree structure formed by the MLFMA,the unknowns are grouped adaptively to realize the domain decomposition,thereby avoiding the creation of artificial interface between domains and reducing the workload of model processing.To ensure the current continuity between domains,the 1/4 impedance on the boundary of the domain is rigorously calculated,with the results more accurate.To deal with the half-space environment,the complex image source is introduced to calculate the near interaction,and the real image source is introduced to calculate the far interaction.A comparison of numerical results from the proposed algorithm and the commercial software FEKO is given,which proves the reliability and accuracy of the algorithm.An out-of-core algorithm is used to store data such as translators in the hard disk,which can significantly reduce the memory consumption.Finally,a numerical example including a 1000-wavelength ship model in half space is given,which demonstrates that the proposed algorithm can simulate electrical large objects under the condition of limited resources.

Key words: domain decomposition method, parallel algorithm, multi level fast multipole algorithm, half-space problem, octree strategy, out-of-core algorithm

中图分类号:

TN820

翟畅,林中朝,赵勋旺,张玉. 一种使用八叉树的半空间MLFMA区域分解算法[J]. 西安电子科技大学学报, 2021, 48(6): 144-150.

ZHAI Chang,LIN Zhongchao,ZHAO Xunwang,ZHANG Yu. Algorithm for half-space MLFMA domain decomposition utilizing an octree[J]. Journal of Xidian University, 2021, 48(6): 144-150.

图/表 8

图1

图2

图3

图4

图5

图6

表1

图7

参考文献 16

[1]	张森, 邓维波, 杨松岩. 电大尺寸卫星目标RCS的PO法计算与分析[J]. 微计算机信息, 2009, 25(2-1):194-196.
	ZHANG Sen, DENG Weibo, YANG Songyan. Calculation and Analysis of Electronically Large Size Satellite Target's RCS Using PO Method[J]. Microcomputer Information, 2009, 25(2-1):194-196.
[2]	吴三元, 邱新军, 马来群. 海面RCS与雷达海上目标检测[J]. 火控雷达技术, 2008, 37(2):51-53.
	WU Sanyuan, QIU Xinjun, MA Laiqun. Sea Surface RCS and Radar Sea Target Detection[J]. Fire Control Radar Technology, 2008, 37(2):51-53.
[3]	TASIC M S, KOLUNDZIJA B M. A PO Driven Iterative Solution of MFIE for Large Scatterers[C]// International Conference on Telecommunication in Modern Satellite.Piscataway,IEEE, 2006:24-27.
[4]	GUO L, FAN T. Octree based backward SBR-PO Method for Electromagnetic Scattering of Electrically Large Target[C]// 2015 IEEE International Conference on Computational Electromagnetics.Piscataway:IEEE, 2015:303-305.
[5]	BALDENSPERGER P, LIU J, JIN J. A Hybrid SBR/FE-BI Technique for Computing the RCS of Electrically Large Objects with Deep Cavities[C]// IEEE Antennas and Propagation Society International Symposium.Piscataway:IEEE, 2001:726-729.
[6]	SU Q, ZUO S, ZHAO X, et al. Parallelized Method with MLFMA and Higher-Order MoM for Analysis of Near-Field of Electrically Large Target in Half-Space[C]// IEEE 4th Asia-Pacific Conference on Antennas and Propagation.Piscataway:IEEE, 2015:593-594.
[7]	WANG J, ZHOU F. Electromagnetic Scattering from 2-D Sea Surface with 3-D Electrically Large Ship by Parallel MLFMA[C]// 2017 IEEE International Geoscience and Remote Sensing Symposium.Piscataway:IEEE, 2017:4922-4925.
[8]	罗晨, 葛勇, 王舒申, 等. 基于矩量法的某轰炸机电磁散射特性研究[J]. 空天防御, 2019(3):66-72.
	LUO Chen, GE Yong, WANG Shushen, et al. Electromagnetic Scattering Characteristics of a Certain Bomber Based on MomM[J]. Aerospace Defense, 2019 (3):66-72.
[9]	满明远, 雷振亚, 谢拥军, 等. 电大目标散射问题的预修正多层快速多极子分析[J]. 西安电子科技大学学报, 2012, 39(2):133-137.
	MAN Mingyuan, LEI Zhenya, XIE Yongjun, et al. Precorrected Multilayer Fast Multipole Analysis for Scattering from Electrically Large Targets[J]. Journal of Xidian University, 2012, 39(2):133-137.
[10]	TALIB M F, RAHMAN A K, ANUAR M S, et al. Analysis of Geometrical Effect in Free Space Optical Communication System Using DDM Technique[C]// 2016 3rd International Conference on Electronic Design.Piscataway:IEEE, 2016:128-131.
[11]	TALIB M F, RAHMAN A K, ANUAR M S, et al. Enhancement of Free Space Optical Communication System Modulation Using DDM Technique[C]// 2016 3rd International Conference on Electronic Design.Piscataway:IEEE, 2016:62-65.
[12]	LI Y, LIN Z, ZHANG Y, et al. Analysis of Large Airborne Waveguide Slotted Array Using DDM[C]// IEEE 4th Asia-Pacific Conference on Antennas and Propagation.Piscataway:IEEE, 2015:597-598.
[13]	PENG Z, HIPTMAIR R, SHAO Y, et al. Domain Decomposition Preconditioning for Surface Integral Equations in Solving Challenging Electromagnetic Scattering Problems[J]. IEEE Transactions on Antennas and Propagation, 2016, 64(1):210-223. doi: 10.1109/TAP.2015.2500908
[14]	苏秦, 顾宗静, 林中朝, 等. 并行区域分解法分析千波长目标散射特性[J]. 西安电子科技大学学报, 2018, 45(6):13-19.
	SU Qin, GU Zongjing, LIN Zhongchao, et al. Analysis of Scattering Characteristics of the One-Thousand-Wavelength Object Using the Parallel IE-NDDM[J]. Journal of Xidian University, 2018, 45(6):13-19.
[15]	张玉. 计算电磁学中的超大规模并行矩量法[M]. 西安: 西安电子科技大学出版社, 2016.
[16]	GENG N. Fast Numerical Techniques in Computational Electromagnetics for Planar-Stratified Media[D]. Karlsruhe:University of Karlsruhe, 2001.