频域电磁计算的IE-FVFD混合方法

doi:10.12305/j.issn.1001-506X.2025.12.05

系统工程与电子技术 ›› 2025, Vol. 47 ›› Issue (12): 4040-4047.doi: 10.12305/j.issn.1001-506X.2025.12.05

• 电子技术 • 上一篇

频域电磁计算的IE-FVFD混合方法

许勇¹, 陈坚强¹^,²^,*, 刘健¹, 刘文君¹, 于新童¹

1. 中国空气动力研究与发展中心计算空气动力研究所，四川绵阳 621000
2. 空天飞行空气动力科学与技术全国重点实验室，四川绵阳 621000

收稿日期:2024-09-24 修回日期:2024-11-22 出版日期:2025-03-13 发布日期:2025-03-13
通讯作者: 陈坚强
作者简介:许　勇（1971—），男，研究员，博士，主要研究方向为气动物理、计算电磁学
刘　健（1984—），男，副研究员，博士，主要研究方向为计算流体力学、线代方程加速算法
刘文君（1994—），女，工程师，硕士，主要研究方向为电磁特性计算、飞行器隐身设计
于新童（1992—），女，初级工程师，硕士，主要研究方向为计算流体力学软件应用测试
基金资助:
国家数值风洞工程资助课题

IE-FVFD hybrid method for frequency domain electromagnetic computation

Yong XU¹, Jianqiang CHEN¹^,²^,*, Jian LIU¹, Wenjun LIU¹, Xintong YU¹

1. Computational Aerodynamics Institute，China Aerodynamics Research and Development Center，Mianyang 621000，China
2. State Key Laboratory of Aerodynamics，Mianyang 621000，China

Received:2024-09-24 Revised:2024-11-22 Online:2025-03-13 Published:2025-03-13
Contact: Jianqiang CHEN

摘要/Abstract

摘要：

为高效模拟电磁散射问题，提出一种积分方程（integral equation，IE）辐射边界条件结合频域有限体积（finite volume frequency domain，FVFD）的电磁场计算方法。FVFD采用全隐格式时间推进和通量迎风分裂计算，IE算法在网格内电流等效面上积分计算截断边界输入波通量，两者以松耦合方式结合。数值实验表明，IE算法对比常用的Mur算法，截断边界到物面距离可缩短到0.5倍波长，模拟多体散射问题可使用整体网格或独立的分离网格，也能模拟复杂外形目标电磁散射。所提方法在保证计算精度同时，可降低网格规模，有效提升效率。

关键词: 频域有限体积法, 积分方程边界条件, 麦克斯韦方程组, 全隐格式, 电磁散射

Abstract:

To efficiently simulate electromagnetic scattering problems, an electromagnetic field calculation method combining integral equation (IE) radiation boundary conditions and finite volume frequency domain (FVFD) in the frequency domain is proposed. FVFD adopts a fully implicit scheme for time advancement and flux upwind splitting calculation, while the IE algorithm integrates the truncated boundary input wave flux on the current equivalent surface within the grid. The two are combined in a loosely coupled manner. Numerical experiments show that compared with the commonly used Mur algorithm, the IE algorithm can shorten the distance from the truncated boundary to the object surface to 0.5 times the wavelength. Simulating multi-body scattering problems can be done using a global grid or independent separated grids, as well as simulating electromagnetic scattering of complex shaped targets. The proposed method can reduce grid size and effectively improve efficiency while ensuring computational accuracy.

Key words: finite volume frequency domain （FVFD） method, integral equation (IE) boundary condition, Maxwell equations, fully implicit scheme, electromagnetic scattering

中图分类号:

O 441.4

许勇, 陈坚强, 刘健, 刘文君, 于新童. 频域电磁计算的IE-FVFD混合方法[J]. 系统工程与电子技术, 2025, 47(12): 4040-4047.

Yong XU, Jianqiang CHEN, Jian LIU, Wenjun LIU, Xintong YU. IE-FVFD hybrid method for frequency domain electromagnetic computation[J]. Systems Engineering and Electronics, 2025, 47(12): 4040-4047.

图/表 13

图1

图2

图3

图4

图5

图6

图7

图8

图9

图10

表1

图11

图12

参考文献 37

1	ENGQUIST B, MAJDA A. Numerical radiation boundary conditions for unsteady transonic flow[J]. Journal of Computational Physics, 1981, 40 (1): 91- 103. doi: 10.1016/0021-9991(81)90201-1
2	MUR G. Absorbing boundary conditions for the finite-difference approximation of the time-domain electromagnetic-field equations[J]. IEEE Trans. on Electromagnetic Compatibility, 1981, EMC-23(4): 377−382.
3	BAYLISS A, TURKEL A. Far field boundary conditions for compressible flows[J]. Journal of Computational Physics, 1982, 48 (2): 182- 199. doi: 10.1016/0021-9991(82)90046-8
4	SOFIA E, JAN N. Exact non-reflecting boundary conditions revisited: well-posedness and stability[J]. Foundations of Computational Mathematics, 2017, 17, 957- 986. doi: 10.1007/s10208-016-9310-3
5	VIANEY V, JACOB C B, ACOSTA S. High order local farfield expansions absorbing boundary conditions for multiple scattering[J]. Journal of Computational Physics, 2022, 460, 111187. doi: 10.1016/j.jcp.2022.111187
6	BERENGER J P. A perfectly matched layer for the absorption of electromagnetic waves[J]. Journal of Computational Physics, 1994, 114 (1): 185- 200.
7	GEDNEY S D. An anisotropic perfectly matched layer-absorbing medium for the truncation of FDTD lattices[J]. IEEE Trans. on Antennas and Propagation, 1996, 44 (12): 1630- 1639. doi: 10.1109/8.546249
8	SANKARAN K, FUMEAUX C, VAHLDIECK R, et al. Cell-centered finite-volume-based perfectly matched layer for time-domain Maxwell system[J]. IEEE Trans. on Microwave Theory and Techniques, 2006, 54 (3): 1269- 1275. doi: 10.1109/TMTT.2006.869704
9	BALSARA D S, et al. An optimized CPML formulation for high order FVTD schemes for CED[J]. IEEE Journal on Multiscale and Multiphysics Computational Techniques, 2021, 6 (1): 183- 200.
10	WANG S Q, WEI X, ZHOU Y G, et al. High-order conformal perfectly matched layer for the DGTD method[J]. IEEE Trans. on Antennas and Propagation, 2021, 69 (11): 7753- 7760. doi: 10.1109/TAP.2021.3084625
11	ZIOLKOWSKI R W, MADSEN N K, CARPENTER R C, et al. Three-dimensional computer modeling of electromagnetic fields: a global lookback lattice truncation scheme[J]. Journal of Computational Physics, 1983, 50 (3): 360- 408. doi: 10.1016/0021-9991(83)90103-1
12	MOERLOOSE J D, ZUTTER D D. Surface integral representation radiation boundary condition for the FDTD method[J]. IEEE Trans. on Antennas and Propagation, 1993, 41 (7): 890- 896. doi: 10.1109/8.237619
13	JIAO D, LU M Y, MICHIELSSEN E, et al. A fast time-domain finite element-boundary integral method for electromagnetic analysis[J]. IEEE Trans. on Antennas and Propagation, 2001, 49 (10): 1453- 1461. doi: 10.1109/8.954934
14	MAYERGOYZ I D, ANDREI P, HAKIM B. A new time-domain approach to the analysis of scattering problems[J]. IEEE Trans. on Magnetics, 2002, 38 (2): 327- 332. doi: 10.1109/20.996089
15	SARTO M S, SCARLATTI A. Efficient numerical calculation of integral equations boundary conditions for the FDTD method[J]. IEEE Trans. on Magnetics, 2003, 39 (3): 1242- 1245. doi: 10.1109/TMAG.2003.810373
16	TENG Z H. Exact boundary condition for time-dependent wave equation based on boundary integral[J]. Journal of Computational Physics, 2003, 190 (2): 398- 418. doi: 10.1016/S0021-9991(03)00281-X
17	SHANKER B, LU M Y, ERGIN A A, et al. Plane-wave time-domain accelerated radiation boundary kernels for FDTD analysis of 3-D electromagnetic phenomena[J]. IEEE Trans. on Antennas and Propagation, 2005, 53 (11): 3704- 3716. doi: 10.1109/TAP.2005.858590
18	FIRSOV D K, LOVETRI J. FVTD—integral equation hybrid for Maxwell’s equations[J]. International Journal of Numerical Modeling: Electronic Networks, Devices and Fields, 2007, 21 (1/2): 29- 42.
19	LIN Y Q, WIELE D S. Multi-region finite-difference time-domain （MR-FDTD） based on domain-optimal Green’s functions[J]. IEEE Trans. on Antennas and Propagation, 2013, 61 (5): 2655- 2663. doi: 10.1109/TAP.2013.2238597
20	TAKAHASHI T. A fast time-domain boundary element method for three-dimensional electromagnetic scattering problems[J]. Journal of Computational Physics, 2023, 482, 112053. doi: 10.1016/j.jcp.2023.112053
21	金建铭. 电磁场有限元方法[M]. 西安: 西安电子科技大学出版社, 1998.
	JIN J M. Finite element method of electromagnetic field[M]. Xi’an: Xidian University Press, 1998.
22	ZHONG Y, WANG H M, HUANG W F, et al. A hybrid loop-tree FEBI method for low-frequency well logging of 3-D structures in layered media[J]. IEEE Trans. on Geoscience and Remote Sensing, 2021, 60, 2001509.
23	SHANKAR V, HALL W F, MOHAMMADIAN A H. A time-domain differential solver for electromagnetics scattering problems[J]. Proceedings of the IEEE, 1989, 77 (5): 709- 721. doi: 10.1109/5.32061
24	SHANG J S, GAITONDE D. Characteristic-based time-dependent Maxwell equations solvers on a general curvilinear frame[J]. AIAA Journal, 1995, 33 (3): 491- 498. doi: 10.2514/3.12376
25	SHANG J S. Time-domain electromagnetic scattering simulations on multicomputers[J]. Journal of Computational Physics, 1996, 128 (2): 381- 390. doi: 10.1006/jcph.1996.0218
26	CAMBEROS J A. Development of finite-volume time-domain, unstructured-grid CEM code on massively parallel platforms[C]//Proc. of the 31th Plasmadynamics & Lasers Conference, 2000.
27	SANKARAN K, FUMEAUX C, VAHLDIECK R. Hybrid PML-ABC truncation techniques for finite-volume time-domain simulations[C]//Proc. of the Asia-Pacific Microwave Conference, 2006.
28	SANKARAN K, FUMEAUX C, VAHLDIECK R, et al. Split and unsplit finite-volume absorbers: formulation and performance comparison[C]//Proc. of the 36th European Microwave Conference: 2006.
29	BAUMANN D, FUMEAUX C, VAHLDIECK R, et al. Conformal perfectly matched absorber for finite-volume time-domain simulations[C]//Proc. of the Asia-Pacific Sympsoium on Electromagnetic Compatibility & 19th International Zurich Symposium on Electromagnetic Compatibility, 2008: 188−191.
30	KROHNE K, BAUMANN D, FUMEAUX C, et al. Frequency-domain finite-volume simulations[C]// Proc. of the 37th European Microwave Conference, 2007.
31	HUH K S, SHU M, AGARWAL R K. A compact high-order finite-volume time-domain/frequency-domain method for electromagnetic scattering[C]//Proc. of the 30th Aerospace Sciences Meeting & Exhibit Conference, 1992.
32	BONNET P, FERRIERES X, GRANDO J, et al. Frequency-domain finite volume method for electromagnetic scattering[C]//Proc. of the IEEE Antennas and Propagation Serenity International Symposium, 1998: 252−255.
33	PULLIAM T H, ZINGG D W. Fundamental algorithms in computational fluid dynamics[M]. Switzerland: Springer, 2014.
34	BLAZEK J. Computational Fluid Dynamics Principles and Applications[M]. 3nd ed. London: Elsevier, 2015.
35	刘久祥, 肖光亮, 郭琨毅, 等. 薄涂覆目标散射中心的精确建模[J]. 系统工程与电子技术, 2023, 45 (6): 1589- 1596.
	LIU J X, XIAO G L, GUO K Y, et al. Accurate scattering center modeling of thin-coated targets[J]. Systems Engineering and Electronics, 2023, 45 (6): 1589- 1596.
36	GIBSON W C. The method of moments in electromagnetics[M]. New York: Chapman and Hall/CRC, 2008.
37	BLAKE D C, BUTER T A. Overset grid methods applied to a finite-volume time-domain maxwell equation solver [C]//Proc. of the 27th Fluid Dynamics, Plasmadynamics and Lasers Conference, 1996: 17−20.

算法	运行时间
IE-FVFD	984.13
Mur-FVFD	1 210.35
MLFMA	441

[1]	王龙, 汪刘应, 刘顾, 葛超群, 张颖, 王文豪, 胡灵杰, 陈海青. 基于PO-PTD优化方法的厢式车辆电磁散射特性研究[J]. 系统工程与电子技术, 2024, 46(12): 3992-3999.
[2]	管灵, 朱晨曦, 董纯柱, 殷红成. 基于模式/区域分解的无人机集群电磁特性快速分析[J]. 系统工程与电子技术, 2021, 43(10): 2697-2705.
[3]	何思远, 朱凌昆, 刘进, 边志丹, 张云华, 朱国强. 雷达目标散射中心正向自动化建模方法研究与实现[J]. 系统工程与电子技术, 2021, 43(10): 2706-2717.
[4]	何姿, 操龙潜, 丁大志, 樊振宏, 陈如山. 典型植被的高效电磁全波分析方法研究[J]. 系统工程与电子技术, 2021, 43(10): 2725-2732.
[5]	杜延磊, 高帆, 刘涛, 杨健. 基于数值仿真的X波段极化SAR海杂波统计建模与特性分析[J]. 系统工程与电子技术, 2021, 43(10): 2742-2755.
[6]	郭杰, 殷红成, 满良. 无源散射单元电磁散射特性可控方法综述[J]. 系统工程与电子技术, 2019, 41(4): 716-723.
[7]	曹欣远，陈明生，吴先良. 矩量法中引入压缩感知求解三维电磁散射问题[J]. Journal of Systems Engineering and Electronics, 2013, 35(6): 1143-1147.
[8]	殷红成，朱国庆，董纯柱，王超. 基于自适应射线管分裂的多次反射计算方法[J]. Journal of Systems Engineering and Electronics, 2013, 35(4): 700-706.
[9]	田炜，任新成，郭立新. 分形海面及其上方目标电磁散射的混合算法[J]. Journal of Systems Engineering and Electronics, 2012, 34(9): 1775-1780.
[10]	余定峰,何思远,付松,朱国强,殷红成,王超. 基于PO-EEC的各向同性介质薄层涂覆目标电磁散射[J]. Journal of Systems Engineering and Electronics, 2012, 34(12): 2428-2432.
[11]	姬金祖, 刘战合, 黄沛霖, 李洁. 反射张量在射线追踪法中的应用[J]. Journal of Systems Engineering and Electronics, 2011, 33(8): 1685-1689.
[12]	张晓燕，盛新庆. 复合目标电磁散射的高效混合计算方法[J]. Journal of Systems Engineering and Electronics, 2011, 33(4): 717-723.
[13]	陈新蕾, 邓小乔, 牛臻弋, 李茁, 顾长青. 电各向异性介质目标散射分析的快速偶极子法[J]. Journal of Systems Engineering and Electronics, 2011, 33(11): 2372-2376.
[14]	丁建军，陈磊，刘志伟，陈如山. 基于时域弹跳射线法分析电大尺寸目标的散射[J]. Journal of Systems Engineering and Electronics, 2010, 32(9): 1846-1849.
[15]	郭立新，麻军，王蕊，刘晓勇. MPI并行矩量法计算二维粗糙面波束电磁散射[J]. Journal of Systems Engineering and Electronics, 2010, 32(9): 1841-1845.

频域电磁计算的IE-FVFD混合方法

IE-FVFD hybrid method for frequency domain electromagnetic computation

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

图/表 13

参考文献 37

相关文章 15

编辑推荐

Metrics

本文评价