CN107016174A - 一种应用于时域有限差分法的透明激励源的实现方法 - Google Patents

一种应用于时域有限差分法的透明激励源的实现方法 Download PDF

Info

Publication number
CN107016174A
CN107016174A CN201710176708.6A CN201710176708A CN107016174A CN 107016174 A CN107016174 A CN 107016174A CN 201710176708 A CN201710176708 A CN 201710176708A CN 107016174 A CN107016174 A CN 107016174A
Authority
CN
China
Prior art keywords
excitation source
field
grid
updating
transparent
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710176708.6A
Other languages
English (en)
Other versions
CN107016174B (zh
Inventor
金晓林
谷晓梁
黄桃
李斌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201710176708.6A priority Critical patent/CN107016174B/zh
Publication of CN107016174A publication Critical patent/CN107016174A/zh
Application granted granted Critical
Publication of CN107016174B publication Critical patent/CN107016174B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Measuring Magnetic Variables (AREA)
  • Complex Calculations (AREA)

Abstract

本发明属于计算电磁学技术领域，具体为一种应用于时域有限差分法的透明激励源的实现方法。本发明采用Yee网格，通过在激励源网格处加入两个幅值相反的电场激励源，减去反馈电场，使之等效为硬源且透明；并在输入及输出端口加完美匹配层吸收边界截断波导结构，适用于一维、二维及三维结构。本发明只需对激励源面上的电磁场分量进行额外的计算，相对于其它方法需要的辅助计算量是微乎其微的，对计算效率的影响可以忽略不计，且具有透明的特性，对于波导类计算模拟问题避免了辅助计算所花费的时间，高效得到入射场和散射场，对提高模拟计算效率具有很好的效果。

Description

一种应用于时域有限差分法的透明激励源的实现方法
技术领域
本发明属于计算电磁学技术领域。在采用时域有限差分法来模拟计算波导类结构的电磁波特性问题时，涉及一种加载在波导端口的透明激励源，具体为一种应用于时域有限差分法的透明激励源的实现方法。
背景技术
在使用时域有限差分法进行模拟计算的过程中，需要加入激励源来激发起特定频段和幅值的电磁波，进而得到模拟计算结构的电磁波传输特性及电磁场分布。
使用非透明激励源激励电磁波，在电磁波传播过程中，当遇到非均匀结构或边界等情况下，会产生散射。散射场会传播至激励源处，非透明激励源相当于导电反射屏，会将散射场反射回计算区域，而这部分反射场是非透明激励源数值实现过程中产生的非物理效应。所以在模拟计算波导类结构的问题时，为了消除或避免上述非透明激励源产生的非物理场反射，需要在波导端口加载一种透明的激励源，这样计算区域内部计算的散射场才能通过激励源网格而不反射回计算区域。
目前，应用于时域有限差分法的激励源主要有以下三种。设波导结构沿z向分布，波导端口位于x-y平面。
第一种激励源为硬源，就是直接在激励源网格处设置电场幅值为我们想要激励起的电场值，具体公式表示为
其中，下标T表示x或y方向，n表示第n个时间步长，f表示所加激励源方程，ksrc表示激励源网格z向的编号。这种方法的缺点是，由于直接设置激励源网格处的电场值，激励源等效为一个理想的导电反射屏，所有散射场都会被反射。克服反射的方法是在散射场传播到激励源网格处之前，将激励源去掉，这种方法对激励源作用的时间长度有非常大的限制。
第二种激励源为透明电流激励源，此种激励源是将激励项看成是有源麦克斯韦方程中的电流项，具体公式表示为
其中，J表示加载的激励电流，表示在T方向的离散方程，具体表达为
这种方法的缺点是无法设置激励源网格处的电场为我们设置的场值。为了计算散射问题，需要进行两次计算，首先通过一次辅助计算，计算一个与实际波导结构端口相同尺寸的无限长均匀波导得到激励源网格处的电场值，即入射场，然后对实际波导结构进行模拟计算得到激励源网格处的总场，减去入射场就得到了散射场。这样通过两次计算的方法会大大降低计算效率。
第三种激励源为透明场激励源，这种方法的主要思想是在麦克斯韦方程中减去反馈场，使其等效为硬源，具体表达式为
其中，In(i,j)需要通过一次辅助计算得到。应用硬源激励的方法，激励源函数采用克罗内克函数δ[n]，计算一个与实际波导结构端口相同尺寸的无限长均匀波导得到激励源网格处的响应，计算公式为
这种方法的缺点也是采用了一次辅助计算，降低了计算效率。
发明内容
针对上述存在问题或不足，为解决两次计算降低计算效率的技术问题，本发明提供了一种应用于时域有限差分法的透明激励源的实现方法。
本发明采用Yee网格，设波导结构沿z向分布，波导端口位于x-y平面。通过在激励源网格处加入两个幅值相反的电场激励源，减去反馈电场，使之等效为硬源且透明；并在输入及输出端口加完美匹配层(PML)吸收边界截断波导结构，如图2所示，图2给出了一个波导类结构问题求解的平面示意图。
本发明适用于一维、二维及三维结构。分布在激励源网格处的电磁场分量具有三套值，分别以下标P、N和R来区别，即(EPx,EPy,HPz)、(ENx,ENy,HNz)和(ERx,ERy,HRz)。P表示与正源相关，N表示与负源相关，R表示与散射场相关。下标PT表示Px或Py，对于NT、RT同理。适用于一维和二维结构时，电磁场分量及其对应公式相应退化。
适用于三维时的具体方案如下：
步骤1、所有分布在激励源网格处的电磁场分量具有三套值，分别以下标P、N和R来区别，即(EPx,EPy,HPz)、(ENx,ENy,HNz)和(ERx,ERy,HRz)。P表示与正源相关，N表示与负源相关，R表示与散射场相关。下标PT表示Px或Py，对于NT、RT同理。
步骤2、
根据激励源方程fT(i,j)来计算表示纵向磁场的辅助项fH(i,j)，由于fT(i,j)为已知项，所有时刻的fH(i,j)项都可以求得
计算所有时刻的反馈电场RNT(i,j)、RPT(i,j)
步骤3、
更新激励源网格点处的电场
其中，表示在求解时，当涉及激励源网格上的磁场时使用正源对应的磁场值，对于同理，具体公式为
散射电场可得
然后应用标准离散麦克斯韦方程更新其它网格点的电场，同时更新完美匹配层吸收边界电场。
步骤4、
更新激励源网格点处的磁场
其中，表示在求解时，当涉及激励源网格上的电场时使用正源对应的电场值，对于同理，具体公式为
接着更新与激励源网格点相邻的横向磁场
其中，
散射磁场可得
然后按照标准离散麦克斯韦方程更新其它网格点的磁场，同时更新完美匹配层吸收边界磁场。
步骤5、
重复3、4步骤，从n＝0迭代直至达到n＝nt，nt为所设置的模拟步数。
适用于二维时的具体方案如下，如图3所示：
步骤1、
所有分布在激励源网格处的电磁场分量具有三套值，即(EPy,HPz)、(ENy,HNz)和(ERy,HRz)。
步骤2、
根据激励源方程fy(i)来计算表示纵向磁场的辅助项fH(i)
计算反馈电场RNy(i)、RPy(i)
步骤3、更新激励源网格点处的电场
其中，
计算散射电场
然后应用标准离散麦克斯韦方程更新其它网格点的电场，同时更新完美匹配层吸收边界电场。
步骤4、更新激励源网格点处的磁场
更新与激励源网格点相邻的横向磁场
计算散射磁场可得
然后按照标准离散麦克斯韦方程更新其它网格点的磁场，同时更新完美匹配层吸收边界磁场。
步骤5、
重复3、4步骤，从n＝0迭代直至达到n＝nt，nt为所设置的模拟步数。
适用于一维时的具体方案如下，如图4所示：
步骤1、
所有分布在激励源网格处的电磁场分量具有三套值，即EPy、ENy和ERy。
步骤2、
更新激励源网格点处的电场
计算散射电场
然后应用标准离散麦克斯韦方程更新其它网格点的电场，同时更新完美匹配层吸收边界电场。
步骤3、
更新与激励源网格点相邻的横向磁场
然后按照标准离散麦克斯韦方程更新其它网格点的磁场，同时更新完美匹配层吸收边界磁场。
步骤4、重复2、3步骤，从n＝0迭代直至达到n＝nt，nt为所设置的模拟步数。
本发明通过在激励源网格处加载一对正负激励源，并减去反馈电场，等效为硬源且透明，从而可以轻易设置入射场并得到散射场。实际上，反馈电场项由横向磁场和纵向磁场共同求得，由坡应廷矢量可以知道与激励源平面相邻的两个方向上的横向磁场幅值相同，使得激励源网格点处的横向磁场沿纵向中心差分为零，而表示纵向磁场的辅助项fH(i,j)通过已知的激励源fT(i,j)得到，从而避免了辅助计算。
如果没有散射场，那么激励源网格处的电场有
联立(7)式和(47)式得
那么会有
对于同理可以得到相同的结论。
如果让激励源为零，即
那么可以得到
此方程即为标准的时域有限差分法的方程，对于同理可以得到相同的结论。根据电磁波的叠加性原理可知，对于散射场来说，此激励源是透明的，任何其它的场都可以穿透激励源网格而不发生反射。
本发明只需对激励源面上的电磁场分量进行额外的计算，相对于其它方法需要的辅助计算量是微乎其微的，对计算效率的影响可以忽略不计。
综上所述，本发明具有透明的特性，对于波导类计算模拟问题避免了辅助计算所花费的时间，高效得到入射场和散射场，对提高模拟计算效率具有很好的效果。
附图说明
图1为三维Yee网格电磁场分量位置分布及编号示意图；
图2为三维透明激励源x-z平面俯视示意图；
图3为二维透明激励源示意图；
图4为一维透明激励源示意图；
图5为两个波导连接散射测试实例示意图；
图6为波导连接散射实例测试点总场EPy和散射电场ERy的幅值结果图；
图7为波导连接散射实例测试点激励源方程fy和实际激励起入射电场Ey的幅值比较图。
具体实施方式
下面结合附图和实施例对本发明作进一步详细说明。
实例中测试了两个二维均匀波导在连接处发生散射的模型，如图5所示。第一个波导宽度为a＝22.86mm，第二个波导宽度为a/2，模拟参数为：空间步长Δx＝Δz＝1.143mm，Δt＝0.0019ns，模拟时间步数nt＝900。激励采用了TE10模式的余弦调制高斯脉冲，中心频率f0＝16GHz，激励源加载在第一个波导输入端口上，所加激励源方程为
其中，n0＝450，nd＝125。
步骤1、
所有分布在激励源网格处的电磁场分量具有三套值，即EPy、HPz，ENy、HNz和ERy、HRz。
步骤2、根据激励源方程fy(i)来计算表示纵向磁场的辅助项fH(i)
计算反馈电场RNy(i)、RPy(i)
步骤3、更新激励源网格点处的电场
计算散射电场
然后应用标准离散麦克斯韦方程更新其余网格点的电场，同时更新完美匹配层吸收边界电场。
步骤4、更新激励源网格点处的磁场
更新与激励源网格点相邻的横向磁场
计算散射磁场可得
然后按照标准离散麦克斯韦方程更新其余网格点的磁场，同时更新完美匹配层吸收边界磁场。
步骤5、
重复3、4步骤，从n＝0迭代直至达到n＝900。
通过一次计算，我们就得到了散射场，图6显示了在激励源平面上一个测试点的总场EPy和散射电场ERy的幅值。通过总场EPy与散射电场ERy相减，就可以得到在激励源网格测试点实际激励起的入射场，如图7。图7给出了激励源方程与实际激励起的入射场在测试点场幅值的比较，可以看出两者完全吻合。
综上可见，本发明的透明激励源具有透明的特性，只需要一次模拟计算就可以同时得到总场、入射场及散射场，避免了辅助计算所花费的时间，大大提高了计算效率。

Claims (3)

1.一种应用于时域有限差分法的透明激励源的实现方法，具体包括以下步骤：
步骤1、采用Yee网格，设波导结构沿z向分布，波导端口位于x-y平面，所有分布在激励源网格处的电磁场分量具有三套值，分别以下标P、N和R来区别，即EPx、EPy、HPz，ENx、ENy、HNz和ERx、ERy、HRz；P表示与正源相关，N表示与负源相关，R表示与散射场相关，下标PT表示Px或Py，对于NT、RT同理；
步骤2、
根据激励源方程fT(i,j)来计算表示纵向磁场的辅助项fH(i,j)，由于fT(i,j)为已知项，所有时刻的fH(i,j)项都可以求得
计算所有时刻的反馈电场RNT(i,j)、RPT(i,j)
步骤3、
更新激励源网格点处的电场
其中，表示在求解时，当涉及激励源网格上的磁场时使用正源对应的磁场值，对于同理，具体公式为
散射电场可得
然后应用标准离散麦克斯韦方程更新其余网格点的电场，同时更新完美匹配层吸收边界电场；
步骤4、
更新激励源网格点处的磁场
其中，表示在求解时，当涉及激励源网格上的电场时使用正源对应的电场值，对于同理，具体公式为
接着更新与激励源网格点相邻的横向磁场
其中，
散射磁场可得
然后按照标准离散麦克斯韦方程更新其余网格点的磁场，同时更新完美匹配层吸收边界磁场；
步骤5、重复3、4步骤，从n＝0迭代直至达到n＝nt，nt为所设置的模拟步数。
2.如权利要求1所述应用于时域有限差分法的透明激励源的实现方法，其特征在于：适用于二维时，分布在激励源网格处的三套电磁场分量分别对应为EPy、HPz，ENy、HNz和ERy、HRz；电磁场分量对应公式相应退化。
3.如权利要求1所述应用于时域有限差分法的透明激励源的实现方法，其特征在于：适用于一维时，分布在激励源网格处的三套电磁场分量对应为EPy、ENy和ERy；电磁场分量对应公式相应退化。
CN201710176708.6A 2017-03-23 2017-03-23 一种应用于时域有限差分法的透明激励源的实现方法 Active CN107016174B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710176708.6A CN107016174B (zh) 2017-03-23 2017-03-23 一种应用于时域有限差分法的透明激励源的实现方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710176708.6A CN107016174B (zh) 2017-03-23 2017-03-23 一种应用于时域有限差分法的透明激励源的实现方法

Publications (2)

Publication Number Publication Date
CN107016174A true CN107016174A (zh) 2017-08-04
CN107016174B CN107016174B (zh) 2020-03-27

Family

ID=59440245

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710176708.6A Active CN107016174B (zh) 2017-03-23 2017-03-23 一种应用于时域有限差分法的透明激励源的实现方法

Country Status (1)

Country Link
CN (1) CN107016174B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109492284A (zh) * 2018-10-30 2019-03-19 电子科技大学 一种波导端口共形卷积完美匹配层吸收边界算法
CN111209249A (zh) * 2020-01-10 2020-05-29 中山大学 一种时域有限差分法硬件加速器架构及其实现方法
CN114996991A (zh) * 2022-05-07 2022-09-02 华中科技大学 一种二维结构散射场的分析方法及系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1144022A (zh) * 1994-03-18 1997-02-26 布朗大学研究基金会 含有强散射增益介质以提供类似激光器作用的光源
CN101819235A (zh) * 2010-01-07 2010-09-01 南京大学 基于时域有限差分分析方法的短波电台电磁脉冲测试电路
CN103310069A (zh) * 2013-06-25 2013-09-18 西安电子科技大学 面向时域有限差分电磁计算的载体网格划分方法
CN103605633A (zh) * 2013-09-22 2014-02-26 西安交通大学 一种粗网格大时间步时域有限差分方法
US20140365188A1 (en) * 2013-06-06 2014-12-11 Acacia Communications Inc. Sparse finite-difference time domain simulation
CN105069247A (zh) * 2015-08-21 2015-11-18 西安空间无线电技术研究所 一种基于时域有限差分的微波部件无源互调数值分析方法
CN105825015A (zh) * 2016-03-18 2016-08-03 中国人民解放军火箭军工程大学 一种用于磁化等离子体的时域有限差分方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1144022A (zh) * 1994-03-18 1997-02-26 布朗大学研究基金会 含有强散射增益介质以提供类似激光器作用的光源
CN101819235A (zh) * 2010-01-07 2010-09-01 南京大学 基于时域有限差分分析方法的短波电台电磁脉冲测试电路
US20140365188A1 (en) * 2013-06-06 2014-12-11 Acacia Communications Inc. Sparse finite-difference time domain simulation
CN103310069A (zh) * 2013-06-25 2013-09-18 西安电子科技大学 面向时域有限差分电磁计算的载体网格划分方法
CN103605633A (zh) * 2013-09-22 2014-02-26 西安交通大学 一种粗网格大时间步时域有限差分方法
CN105069247A (zh) * 2015-08-21 2015-11-18 西安空间无线电技术研究所 一种基于时域有限差分的微波部件无源互调数值分析方法
CN105825015A (zh) * 2016-03-18 2016-08-03 中国人民解放军火箭军工程大学 一种用于磁化等离子体的时域有限差分方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
R.M.JOSEPH 等: "FDTD Maxwell"s equations models for nonlinear electrodynamic and optics", 《IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》 *
张建萍: "二维多区域伪谱时域算法的研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109492284A (zh) * 2018-10-30 2019-03-19 电子科技大学 一种波导端口共形卷积完美匹配层吸收边界算法
CN109492284B (zh) * 2018-10-30 2022-05-03 电子科技大学 一种波导端口共形卷积完美匹配层吸收边界算法
CN111209249A (zh) * 2020-01-10 2020-05-29 中山大学 一种时域有限差分法硬件加速器架构及其实现方法
CN111209249B (zh) * 2020-01-10 2021-11-02 中山大学 一种时域有限差分法硬件加速器系统及其实现方法
CN114996991A (zh) * 2022-05-07 2022-09-02 华中科技大学 一种二维结构散射场的分析方法及系统
CN114996991B (zh) * 2022-05-07 2024-06-21 华中科技大学 一种二维结构散射场的分析方法及系统

Also Published As

Publication number Publication date
CN107016174B (zh) 2020-03-27

Similar Documents

Publication Publication Date Title
Taflove et al. Computational electromagnetics: the finite-difference time-domain method
CN107016174B (zh) 一种应用于时域有限差分法的透明激励源的实现方法
Sirenko et al. Incorporation of exact boundary conditions into a discontinuous Galerkin finite element method for accurately solving 2D time-dependent Maxwell equations
Blanc et al. Wave simulation in 2D heterogeneous transversely isotropic porous media with fractional attenuation: a Cartesian grid approach
Riano On persistence properties in weighted spaces for solutions of the fractional Korteweg–de Vries equation
Mukhin Analysis of the dispersion characteristics in the vicinity of the passband boundaries of the slow-wave structures that represent chained cavities
Iemma et al. On the use of the analogue transformation acoustics in aeroacoustics
CN113987792A (zh) 一种fdtd算法中实现模式源精确输入的方法
Nikolova et al. $ S $-parameter sensitivities for electromagnetic optimization based on volume field solutions
Yevick et al. Impedance-matched absorbers for finite-difference parabolic equation algorithms
Nikkhah et al. Inverse-designed metastructures together with reconfigurable couplers to compute forward scattering
Garcia et al. A new hybrid DGTD/FDTD method in 2-D
Pereira et al. Coupled plate energy models at mid-and high-frequency vibrations
Wu et al. Analysis of resonance in jet screech with large-eddy simulations
Hayashi et al. Combination of beam propagation method and mode expansion propagation method for bidirectional optical beam propagation analysis
Cicuttin et al. Numerical investigation of a 3D hybrid high-order method for the indefinite time-harmonic Maxwell problem
Kondrat'ev et al. Excitation and propagation of electromagnetic waves in nonuniform density ducts
Hagness et al. Finite-difference time-domain methods
Malevsky et al. Source decomposition as a diakoptic boundary condition in FDTD with reflecting external regions
Hsiao et al. An algorithm for beam propagation method in matrix form
Doc et al. Bremmer series for the multimodal sound propagation in inhomogeneous waveguides
Bleszynski et al. Nonsingular Laplacian representation of the asymptotic part of the layered medium Green function in the mixed potential formulation
Mayfield et al. AN ASYMPTOTIC GREEN’S FUNCTION METHOD FOR VECTOR WAVE EQUATIONS
Curt et al. Enhanced functionality for hardware-based FDTD accelerators
Zekios et al. A finite element eigenanalysis of electrically large THz and photonic structures

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Jin Xiaolin

Inventor after: Gu Xiaoliang

Inventor after: Huang Tao

Inventor after: Yang Zhonghai

Inventor after: Li Bin

Inventor before: Jin Xiaolin

Inventor before: Gu Xiaoliang

Inventor before: Huang Tao

Inventor before: Li Bin

SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant