CN107016184A - 一种二维高精度迭代的非磁化等离子体中的实现方法 - Google Patents
一种二维高精度迭代的非磁化等离子体中的实现方法 Download PDFInfo
- Publication number
- CN107016184A CN107016184A CN201710205431.5A CN201710205431A CN107016184A CN 107016184 A CN107016184 A CN 107016184A CN 201710205431 A CN201710205431 A CN 201710205431A CN 107016184 A CN107016184 A CN 107016184A
- Authority
- CN
- China
- Prior art keywords
- field component
- ηmax
- component coefficient
- coefficient
- electric field
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 86
- 238000013316 zoning Methods 0.000 claims abstract description 61
- 230000005684 electric field Effects 0.000 claims abstract description 43
- 230000005672 electromagnetic field Effects 0.000 claims abstract description 35
- 230000000977 initiatory effect Effects 0.000 claims abstract description 4
- 238000004088 simulation Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000009795 derivation Methods 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 6
- 239000011159 matrix material Substances 0.000 description 7
- 230000008569 process Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000000205 computational method Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 108010003272 Hyaluronate lyase Proteins 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000005495 cold plasma Effects 0.000 description 1
- 238000004513 sizing Methods 0.000 description 1
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)
- Hall/Mr Elements (AREA)
- Plasma Technology (AREA)
Abstract
本发明公开了一种二维高精度迭代的非磁化等离子体中的实现方法,包括:输入模型文件;初始化参数和设置参数;添加场源到y方向上的电场分量系数中,设置电场分量系数记为初始场值更新计算整个计算区域的y方向上电场分量系数更新计算整个计算区域的x方向上电场分量系数判断迭代次数k是否达到预设值;更新计算整个计算区域的磁场分量系数更新计算整个计算区域的极化电流密度系数更新计算整个计算区域的电磁场分量系数的辅助变量;更新计算观测点处的电磁场分量;判断拉盖尔多项式的阶数q是否达到预设值。本发明一种二维高精度迭代的非磁化等离子体中的实现方法,计算精度高、计算速度快,且对于低频和凋落波具有很好的吸收效果。
Description
技术领域
本发明属于计算电磁学技术领域,具体涉及一种二维高精度迭代的非磁化等离子体中的实现方法。
背景技术
时域有限差分(Finite-difference time-domain,FDTD)方法因其计算简单、容易实现等优点,被广泛用于色散媒质的电磁波传播的仿真中。但是,它的时间步长受柯西稳定性条件的限制,不能选取的较大,在多尺寸复杂精细结构模型中,FDTD方法计算速度慢,计算效率低。为了消除柯西稳定性条件的限制,人们提出了无条件稳定时域有限差分方法,比如:交替方向隐式(Alternating-Direction-Implicit,ADI)的时域有限差分(ADI-FDTD)方法和基于加权拉盖尔多项式的时域有限差分(Weighted-Laguerre-polynomials Finite-difference time-domain,WLP-FDTD)方法。在这些方法中,ADI-FDTD方法在使用较大的时间步长时会产生很大的色散误差,而WLP-FDTD方法既能消除柯西稳定性条件的限制,又能解决ADI-FDTD方法在使用较大的时间步长时会产生很大的色散误差这个难题,因此WLP-FDTD方法可以高效的求解等离子体中的电磁问题。然而,这种WLP-FDTD方法在求解电磁场过程中,会产生一个大型的稀疏矩阵方程,直接求解此方程会使得计算较复杂,内存消耗较大,于是提出了一种因式分裂的WLP-FDTD方法,该方法在计算速度和计算效率上得到了较大的提高,但是由于该方法是通过添加微扰项,进行因式分裂得来的,计算中会产生分裂误差,为了减小分裂误差、提高计算精度,同时保证较快的计算速度,提出了一种迭代的加权拉盖尔多项式时域有限差分方法。
而由于计算机容量的限制,电磁场的计算只能在有限区域进行。为了能模拟开域电磁波传播过程,必须在计算区域的截断边界处给出吸收边界条件。有人提出了完全匹配层(Perfectly matched layer,PML)吸收边界,后来PML被广泛应用于计算区域的截断,而且被证明是非常有效的,但是研究发现这种传统PML对低频以及凋落波的吸收效果并不理想;使用带有复频率偏移(Complex frequency shift,CFS)因子的PML(CFS-PML)吸收边界可以有效地改善传统PML对低频,凋落波与掠射情况的吸收效果。最近,有人提出了一种使用辅助微分方程的近似完全匹配吸收边界的WLP-FDTD方法,来解色散媒质中的电磁场问题,这种近似完全匹配吸收边界的吸收效果非常差、计算时存在误差,而且此算法计算时间长、内存消耗较大。
发明内容
本发明的目的是提供一种二维高精度迭代的非磁化等离子体中的实现方法,计算速度快、内存消耗小、精度高,且对于低频和凋落波具有很好的吸收效果。
本发明所采用的技术方案是,一种二维高精度迭代的非磁化等离子体中的实现方法,按照以下步骤实施:
步骤1:输入模型文件;
步骤2:初始化参数和设置参数;
步骤3:添加场源到y方向上的电场分量系数中,设置电场分量系数记为初始场值
步骤4:更新计算整个计算区域的y方向上电场分量系数
步骤5:更新计算整个计算区域的x方向上电场分量系数
步骤6:将k+1赋值给k,并判断迭代次数k是否达到预设值,若未达到预设值,则返回步骤4,若达到预设值,则执行步骤7;
步骤7:更新计算整个计算区域的磁场分量系数
步骤8:更新计算整个计算区域的极化电流密度系数
步骤9:更新计算整个计算区域的电磁场分量系数的辅助变量;
步骤10:更新计算观测点处的电磁场分量;
步骤11:将q+1赋值给q,并判断拉盖尔多项式的阶数q是否达到预设值,若未达到预设值,则返回步骤3,若达到预设值,则结束。
本发明的特点还在于:
步骤1具体为:
计算区域大小Nx×Ny,其中Nx为x方向的网格数,Ny为y方向的网格数;空间步长Δη,η=x,y,x为横坐标,y为纵坐标;时间步长Δt;真空中的电导率σ,磁导率μ0,介电常数ε0;等离子体碰撞频率υ与等离子体频率ωp;等离子体在计算区域中的位置;吸收边界层数NPML与相关参数κηmax,αηmax,σηmax;κηmax取整数,κηmax取值范围为[1,60];αηmax取值范围为[0,1);σηmax/σopt取值范围为(0,12],σopt=(m+1)/150πΔh,m取值范围为[1,20],Δη取值范围为λ为源的波长;仿真计算时长Tf;迭代次数k,k≥0且为整数;加权拉盖尔多项式的阶数q,q ≥0且为整数;时间尺度因子s,s取值范围为[109,1013];观测点;场源参数。
步骤2初始化的参数具体为:
将整个计算区域的电磁场分量系数整个计算区域的极化电流密度系数整个计算区域的电磁场分量系数的和整个计算区域的极化电流密度系数的和整个计算区域的辅助变量和其中表示 表示拉盖尔多项式其中全部初始化为零;非磁化等离子体参数C5,C6初始化为C5=0,C6=2;PML系数C1η,C2η,C3,C4初始化为C1η=2/ε0s,C2η=1,C3=e0/μ0,C4=2/(ε0s),式中,η=x,y,x为横坐标,y为纵坐标,ε0是空气中的介电常数,s为时间尺度因子,取值范围为[109,1013],μ0是空气中的磁导率。
设置的参数具体为:
设置带有CFS因子的SC-PML吸收边界的参数ση,κη,αη;具体为:
ση=σηmax|η-η0|m/dm;
κη=1+(κηmax-1)|η-η0|m/dm;
αη=αηmax;
其中,η=x,y,η0为PML层与非PML截面位置,d是PML吸收边界的厚度,κηmax取整数,κηmax取值范围为[1,60];αηmax取值范围为[0,1);σηmax根据σopt来设置,σηmax/σopt取值范围为(0,12],σopt=(m+1)/150πΔη,m取值范围为[1,20],Δη取值范围为λ为源的波长;
设置PML系数C1η,C2η和与等离子参数相关的系数C5,C6;具体为:
C1η=1/(κηαη+ση+0.5κηε0s),C2η=(2αη/ε0s+1);
步骤3中所添加的场源的表达式为:
Jy(t)=(t-t0)/t×exp(-(t-t0)2/t2)
其中,t表示时间,单位为秒;t0,t为场源参数。
步骤4具体为:
步骤4.1:电场分量系数在计算区域的方程为:
其中,k表示迭代次数,i表示x轴方向上第i个计算网格的位置,j表示y轴方向上第j个计算网格的位置,i+1/2表示x轴方向上第i个半网格的位置,j表示y轴方向上第j个半网格的位置;C1x|i表示系数C1x在x轴方向上第i个网格处的值;表示第k+1次迭代x轴方向上第i-1个网格y轴方向上第j个半网格处的电场分量系数的值;
步骤4.2:使用追赶法求解步骤4.1的方程,得到整个计算区域的电场分量系数
步骤5具体为:
步骤5.1:电场分量系数在计算区域的方程为:
步骤5.2:使用追赶法求解步骤5.1的方程,得到整个计算区域的电场分量系数
步骤7具体更新方程为:
步骤8更新公式具体为:
步骤9更新公式具体为:
其中表示 表示电场分量系数的和对y的导数、电场分量系数的和对x的导数、磁场分量系数的和对x或y的求导,表示辅助变量 表示辅助变量的和。
步骤10更新公式具体为:
其中,U表示电磁场分量Ex,Ey,Hz,r表示电磁场分量的位置,Uq表示q阶电磁场分量系数,是q阶加权拉盖尔多项式,是带有时间尺度因子s>0的扩展时间,是q阶拉盖尔多项式。
本发明的有益效果是:
1).本发明一种二维高精度迭代的非磁化等离子体中的实现方法,在直角坐标系下,通过用加权拉盖尔多项式表示电磁场分量,来解时域麦克斯韦方程,使得在更新计算整个计算区域的电磁场分量系数时不涉及到时间步长,只是在最后计算观测点处的电磁场分量时用到时间步长,因此计算过程中时间步长可以取得比柯西稳定性条件限制的时间步长更大;
2).本发明一种二维高精度迭代的非磁化等离子体中的实现方法,在求解非磁化等离子体中的电磁场分量系数时,使用迭代的因式分裂的方案将大型稀疏矩阵方程分裂成两个迭代的三对角矩阵方程,使得它在计算时比WLP-FDTD方法更简单、计算速度更快、内存消耗更少而且比无迭代的因式分裂的WLP-FDTD方法精度更高;
3).本发明一种二维高精度迭代的非磁化等离子体中的实现方法,在设置PML系数时,由于采用了CFS因子,并且通过调整CFS因子中的参数,可以使得该吸收边界对低频与凋落波的吸收更加有效;
4).本发明一种二维高精度迭代的非磁化等离子体中的实现方法,由于采用了复扩展坐标系,使得PML在实现时避免了场的分裂且与媒质无关。
附图说明
图1是本发明非磁化等离子体中的实现方法的流程示意图;
图2是本发明非磁化等离子体中点源辐射的计算模型的示意图;
图3是本发明的方法与传统的FDTD方法、因式分裂的WLP-FDTD方法在观测点处y方向上电场分量时域波形对比图;
图4是因式分裂的WLP-FDTD方法和本发明的方法在观测点处y方向上电场分量误差对比图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明的一种二维高精度迭代的非磁化等离子体中的实现方法,原理为:首先导出二维非磁化等离子体中,电磁场所满足的复扩展坐标系下的麦克斯韦方程,然后使用二维非磁化等离子体中迭代的因式分裂的WLP-FDTD方法推导出整个计算区域的电磁场分量系数和电流密度分量系数的更新方程,接着使用追赶法求得观测点处的电磁场分量系数,最后采用公式(20)求解观测点处的电磁场分量。
在求解二维非磁化等离子体中电磁波传播所满足的更新方程时,首先需要推导出复扩展坐标系下,PML中电磁场满足的麦克斯韦方程;
碰撞冷等离子体色散媒质中,扩展坐标下,麦克斯韦方程组和相关的联立方程为
式中,H是磁场强度;E是电场强度;J是极化电流密度;ε0、μ0分别为真空中的介电常数和磁导率;ωp是等离子体频率;υ是等离子体碰撞频率。为修正后的微分算子,可以写成
sx、sy和sz是坐标扩展变量,可以表示成
Sη=1+ση/jωε0 (5)
加入CFS因子后,可以表示成
Sη=kη+ση/(αη+jωε0) (6)
其中η=(x,y,z),kη,ση和αη为PML的有关的参数。
应用扩展坐标的CFS-PML,仅考虑二维TEz的情况,麦克斯韦方程组和相关的联立方程可化为:
下面我们引入几个辅助变量如下:
将(6)式分别代入式(12)-(15),然后做如下变化于是可得下面几式:
我们可以知道
式中U代表Ex、Ey、Hz,是加权拉盖尔多项式, p≥0;t≥0是p阶拉盖尔多项式。
将(20)、(21)代入(16)-(19),然后应用加勒金的测试过程,可得
式中
s>0是时间尺度因子,q是加权拉盖尔多项式的阶数。
将(20)、(21)代入(7)—(11),再应用加勒金的测试过程得到:
将(34)代入(31),(35)代入(32)后得到
式中
将(36)、(37)和(33)式写成矩阵形式如下
式中
如果让那么(43)式可以写为
添加一个微扰项到上式,于是得到
为了强调迭代的特点,将(48)式重写为
上式可以分裂为下面两式
式中是一个非物理中间量,为了解(50)式,我们选择
式中
将(51)式代入(50)式,化简后得到
将上式扩展得到
将(54)式的第四式代入第二式和第五式,第一式和第四式代入第三式得到
对上式进行中心差分得到
上面五式中,i表示x轴方向上第i个计算网格的位置,j表示y轴方向上第j个计算网格的位置,i+1/2表示x轴方向上第i个半网格的位置,j表示y轴方向上第j个半网格的位置,k表示第k次迭代,C1x|i表示系数C1x在x轴方向上第i个网格处的值;表示第k+1次迭代x轴方向上第i-1个网格y轴方向上第j个半网格处的电场分量系数的值;在整个计算区域上,(56)式和(57)式可以写成三对角矩阵差分方程,与WLP-FDTD方法相比,这种迭代的因式分解的WLP-FDTD方法将大型稀疏矩阵方程的求解转变成两个三对角矩阵方程的求解,于是可以使用追赶法,非常简单的解得整个计算区域电磁场分量系数,最后通过公式(20)解得观测点的电磁场分量,而且这种方法添加了迭代的方案,使得其计算精度比因式分裂的WLP-FDTD要高。
本发明一种二维高精度迭代的非磁化等离子体中的实现方法,流程示意图如图1所示,按照以下步骤实施:
步骤1:输入模型文件,具体为:
计算区域大小Nx×Ny,其中Nx为x方向的网格数,Ny为y方向的网格数;空间步长Δη,η=x,y,x为横坐标,y为纵坐标;时间步长Δt;真空中的电导率σ,磁导率μ0,介电常数ε0;等离子体碰撞频率υ与等离子体频率ωp;等离子体在计算区域中的位置;吸收边界层数NPML与相关参数κηmax,αηmax,σηmax;κηmax取整数,κηmax取值范围为[1,60];αηmax取值范围为[0,1);σηmax/σopt取值范围为(0,12],σopt=(m+1)/150πΔη,m取值范围为[1,20],Δη取值范围为λ为源的波长;仿真计算时长Tf;迭代次数k,k≥0且为整数;加权拉盖尔多项式的阶数q,q≥0且为整数;时间尺度因子s,s取值范围为[109,1013];观测点;场源参数。
步骤2:初始化参数和设置参数,具体为:
将整个计算区域的电磁场分量系数整个计算区域的极化电流密度系数整个计算区域的电磁场分量系数的和整个计算区域的极化电流密度系数的和整个计算区域的辅助变量和其中表示 表示拉盖尔多项式其中全部初始化为零;非磁化等离子体参数C5,C6初始化为C5=0,C6=2;PML系数C1η,C2η,C3,C4初始化为C1η=2/ε0s,C2η=1,C3=ε0/μ0,C4=2/(ε0s),式中,η=x,y,x为横坐标,y为纵坐标,ε0是空气中的介电常数,s为时间尺度因子,取值范围为[109,1013],μ0是空气中的磁导率。
设置的参数具体为:
设置带有CFS因子的SC-PML吸收边界的参数ση,κη,αη;具体为:
ση=σηmax|η-η0|m/dm;
κη=1+(κηmax-1)|η-η0|m/dm;
αη=αηmax;
其中,η=x,y,η0为PML层与非PML截面位置,d是PML吸收边界的厚度,κηmax取整数,κηmax取值范围为[1,60];αηmax取值范围为[0,1);σηmax根据σopt来设置,σηmax/σopt取值范围为(0,12],σopt=(m+1)/150pΔη,m取值范围为[1,20],Δη取值范围为λ为源的波长;
设置PML系数C1η,C2η和与等离子参数相关的系数C5,C6;具体为:
C1η=1/(κηαη+ση+0.5κηε0s),C2η=(2αη/ε0s+1);
步骤3:添加场源到y方向上的电场分量系数中,设置电场分量系数记为初始场值
步骤3中所添加的场源的表达式为:
Jy(t)=(t-t0)/t×exp(-(t-t0)2/t2)
其中,t表示时间,单位为秒;t0,t为场源参数。
步骤4:更新计算整个计算区域的y方向上电场分量系数具体为:
步骤4.1:电场分量系数在计算区域的方程为:
其中,k表示迭代次数,i表示x轴方向上第i个计算网格的位置,j表示y轴方向上第j个计算网格的位置,i+1/2表示x轴方向上第i个半网格的位置,j表示y轴方向上第j个半网格的位置;C1x|i表示系数C1x在x轴方向上第i个网格处的值;表示第k+1次迭代x轴方向上第i-1个网格y轴方向上第j个半网格处的电场分量系数的值。
步骤4.2:使用追赶法求解步骤4.1的方程,得到整个计算区域的电场分量系数
步骤5:更新计算整个计算区域的x方向上电场分量系数具体为:
步骤5.1:电场分量系数在计算区域的方程为:
步骤5.2:使用追赶法求解步骤5.1的方程,得到整个计算区域的电场分量系数
步骤6:将k+1赋值给k,并判断迭代次数k是否达到预设值,若未达到预设值,则返回步骤4,若达到预设值,则执行步骤7;
步骤7:更新计算整个计算区域的磁场分量系数具体更新方程为:
步骤8:更新计算整个计算区域的极化电流密度系数更新公式具体为:
步骤9:更新计算整个计算区域的电磁场分量系数的辅助变量,更新公式具体为:
其中表示 表示电场分量系数的和对y的导数、电场分量系数的和对x的导数、磁场分量系数的和对x或y的求导,表示辅助变量 表示辅助变量的和。
步骤10:更新计算观测点处的电磁场分量,更新公式具体为:
其中,U表示电磁场分量Ex,Ey,Hz,Uq表示q阶电磁场分量系数,是q阶加权拉盖尔多项式,是带有时间尺度因子s>0的扩展时间,是q阶拉盖尔多项式。
步骤11:将q+1赋值给q,并判断拉盖尔多项式的阶数q是否达到预设值,若未达到预设值,则返回步骤3,若达到预设值,则结束。
下面通过实验对本发明的效果进行说明:
实验:非磁化等离子体中点源辐射的计算
以图2所示的等离子体中点源辐射的计算模型为例,按照本发明的方法步骤进行实施,实验中整个计算区域为50×50网格,网格大小为0.15mm×0.15mm,即Δx=Δy=0.15mm,等离子体充满整个计算区域,其参数ωp=1.261×1011rad/s,υ=2×1010rad/s。四个边界采用10层网格的PML吸收边界,计算中所加的源位于网格(25,25),所加场源的表达式如下:
Jy(t)=(t-t0)/t×exp(-(t-t0)2/t2) (61)
其中,t0=0.05ns,t=0.01ns。观测点位于(38,38)网格处。时间步长Δt=0.25ps,加权拉盖尔多项式的阶数q=245,时间扩展因子s=1.256×1012,整个仿真时间为Tf=0.75ns,迭代次数k=2,PML吸收边界参数κηmax=1,σηmax=σopt,αηmax=0。采用本发明方法计算的观测点处的电场分量Ey与采用传统FDFD方法和因式分裂的WLP-FDTD方法计算的结果参见图3。采用因式分裂的WLP-FDTD方法和本发明的方法计算观测点处y方向上电场分量误差结果参见图4。从图3可以看出传统FDTD方法与本发明方法计算结果一致,验证了本发明方法的正确性。更进一步的我们可以从图4清晰地发现本发明方法的计算精度相比于因式分裂的WLP-FDTD方法的计算精度要高很多,尤其是在仿真时间的后期。在运行计算时间上,传统FDTD方法需要4.382秒,因式分裂的WLP-FDTD方法需要0.437秒,本发明专利中的计算方法需要0.904秒,由此可知本发明专利中的计算方法相对于FDTD方法计算速度得到较大的提高,虽然计算速度比因式分裂的WLP-FDTD方法要稍慢点但精度比因式分裂的WLP-FDTD方法要高。
Claims (10)
1.一种二维高精度迭代的非磁化等离子体中的实现方法,其特征在于,按照以下步骤实施:
步骤1:输入模型文件;
步骤2:初始化参数和设置参数;
步骤3:添加场源到y方向上的电场分量系数中,设置电场分量系数记为初始场值
步骤4:更新计算整个计算区域的y方向上电场分量系数
步骤5:更新计算整个计算区域的x方向上电场分量系数
步骤6:将k+1赋值给k,并判断迭代次数k是否达到预设值,若未达到预设值,则返回步骤4,若达到预设值,则执行步骤7;
步骤7:更新计算整个计算区域的磁场分量系数
步骤8:更新计算整个计算区域的极化电流密度系数
步骤9:更新计算整个计算区域的电磁场分量系数的辅助变量;
步骤10:更新计算观测点处的电磁场分量;
步骤11:将q+1赋值给q,并判断拉盖尔多项式的阶数q是否达到预设值,若未达到预设值,则返回步骤3,若达到预设值,则结束。
2.根据权利要求1所述的一种二维高精度迭代的非磁化等离子体中的实现方法,其特征在于,所述步骤1具体为:
计算区域大小Nx×Ny,其中Nx为x方向的网格数,Ny为y方向的网格数;空间步长Δη,η=x,y,x为横坐标,y为纵坐标;时间步长Δt;真空中的电导率σ,磁导率μ0,介电常数ε0;等离子体碰撞频率u与等离子体频率ωp;等离子体在计算区域中的位置;吸收边界层数NPML与相关参数κηmax,αηmax,σηmax;κηmax取整数,κηmax取值范围为[1,60];αηmax取值范围为[0,1);σηmax/σopt取值范围为(0,12],σopt=(m+1)/150πΔη,m取值范围为[1,20],Δη取值范围为λ为源的波长;仿真计算时长Tf;迭代次数k,k≥0且为整数;加权拉盖尔多项式的阶数q,q≥0且为整数;时间尺度因子s,s取值范围为[109,1013];观测点;场源参数。
3.根据权利要求2所述的一种二维高精度迭代的非磁化等离子体中的实现方法,其特征在于,所述步骤2初始化的参数具体为:
将整个计算区域的电磁场分量系数整个计算区域的极化电流密度系数整个计算区域的电磁场分量系数的和整个计算区域的极化电流密度系数的和整个计算区域的辅助变量和其中表示 表示拉盖尔多项式其中全部初始化为零;非磁化等离子体参数C5,C6初始化为C5=0,C6=2;PML系数C1η,C2η,C3,C4初始化为C1η=2/ε0s,C2η=1,C3=ε0/μ0,C4=2(ε0s),式中,η=x,y,x为横坐标,y为纵坐标,ε0是空气中的介电常数,s为时间尺度因子,取值范围为[109,1013],μ0是空气中的磁导率;
设置的参数具体为:
设置带有CFS因子的SC-PML吸收边界的参数ση,κη,αη;具体为:
ση=σηmax|η-η0|m/dm;
κη=1+(κηmax-1)|η-η0|m/dm;
αη=αηmax;
其中,η=x,y,η0为PML层与非PML截面位置,d是PML吸收边界的厚度,κηmax取整数,κηmax取值范围为[1,60];αηmax取值范围为[0,1);σηmax根据σopt来设置,σηmax/σopt取值范围为(0,12],σopt=(m+1)/150πΔη,m取值范围为[1,20],Δη取值范围为λ为源的波长;
设置PML系数C1η,C2η和与等离子参数相关的系数C5,C6;具体为:
C1η=1/(κηαη+ση+0.5κηε0s),C2η=(2αη/ε0s+1);
4.根据权利要求1所述的一种二维高精度迭代的非磁化等离子体中的实现方法,其特征在于,所述步骤3中所添加的场源的表达式为:
Jy(t)=(t-t0)t×exp(-(t-t0)2/t2)
其中,t表示时间,单位为秒;t0,t为场源参数。
5.根据权利要求1所述的一种二维高精度迭代的非磁化等离子体中的实现方法,其特征在于,所述步骤4具体为:
步骤4.1:电场分量系数在计算区域的方程为:
其中,k表示迭代次数,i表示x轴方向上第i个计算网格的位置,j表示y轴方向上第j个计算网格的位置,i+1/2表示x轴方向上第i个半网格的位置,j表示y轴方向上第j个半网格的位置;C1x|i表示系数C1x在x轴方向上第i个网格处的值;表示第k+1次迭代x轴方向上第i-1个网格y轴方向上第j个半网格处的电场分量系数的值;
步骤4.2:使用追赶法求解步骤4.1的方程,得到整个计算区域的电场分量系数
6.根据权利要求1所述的一种二维高精度迭代的非磁化等离子体中的实现方法,其特征在于,所述步骤5具体为:
步骤5.1:电场分量系数在计算区域的方程为:
步骤5.2:使用追赶法求解步骤5.1的方程,得到整个计算区域的电场分量系数
7.根据权利要求1所述的一种二维高精度迭代的非磁化等离子体中的实现方法,其特征在于,所述步骤7具体更新方程为:
8.根据权利要求1所述的一种二维高精度迭代的非磁化等离子体中的实现方法,其特征在于,所述步骤8更新公式具体为:
9.根据权利要求1所述的一种二维高精度迭代的非磁化等离子体中的实现方法,其特征在于,所述步骤9更新公式具体为:
其中表示 表示电场分量系数的和对y的导数、电场分量系数的和对x的导数、磁场分量系数的和对x或y的求导,表示辅助变量 表示辅助变量的和。
10.根据权利要求1所述的一种二维高精度迭代的非磁化等离子体中的实现方法,其特征在于,所述步骤10更新公式具体为:
其中,U表示电磁场分量Ex,Ey,Hz,r表示电磁场分量的位置,Uq表示q阶电磁场分量系数,是q阶加权拉盖尔多项式,是带有时间尺度因子s>0的扩展时间,是q阶拉盖尔多项式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710205431.5A CN107016184B (zh) | 2017-03-31 | 2017-03-31 | 一种二维高精度迭代的非磁化等离子体中的实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710205431.5A CN107016184B (zh) | 2017-03-31 | 2017-03-31 | 一种二维高精度迭代的非磁化等离子体中的实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107016184A true CN107016184A (zh) | 2017-08-04 |
CN107016184B CN107016184B (zh) | 2021-02-12 |
Family
ID=59444885
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710205431.5A Active CN107016184B (zh) | 2017-03-31 | 2017-03-31 | 一种二维高精度迭代的非磁化等离子体中的实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107016184B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111259587A (zh) * | 2020-01-16 | 2020-06-09 | 中国人民解放军陆军工程大学 | 一种三维交替迭代无条件稳定fdtd算法 |
CN111783339A (zh) * | 2020-06-30 | 2020-10-16 | 西安理工大学 | 电磁波在随机色散介质中传播的pce-fdtd方法 |
CN113887104A (zh) * | 2021-10-10 | 2022-01-04 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种基于mhd的电磁波与磁化等离子体作用数值模拟方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101449337A (zh) * | 2006-06-01 | 2009-06-03 | 约翰·塞尔登·奥格莱 | 利用多个相交的正离子存储环的先进核聚变方法及组件 |
CN104750990A (zh) * | 2015-03-30 | 2015-07-01 | 西安理工大学 | 二维等离子体中扩展坐标的完全匹配吸收边界的实现方法 |
CN104794289A (zh) * | 2015-04-23 | 2015-07-22 | 西安理工大学 | 一种扩展直角坐标系下完全匹配吸收边界的实现方法 |
CN104809286A (zh) * | 2015-04-23 | 2015-07-29 | 西安理工大学 | 一种等离子体中扩展坐标的完全匹配吸收边界的实现方法 |
CN104809343A (zh) * | 2015-04-23 | 2015-07-29 | 西安理工大学 | 一种等离子体中使用电流密度卷积完全匹配层的实现方法 |
CN104820660A (zh) * | 2015-04-23 | 2015-08-05 | 西安理工大学 | 一种扩展柱坐标系下完全匹配吸收边界的实现方法 |
US20160097722A1 (en) * | 2014-10-07 | 2016-04-07 | Unist Academy-Industry Research Corporation | Plasma diagnostic method and apparatus using raman scattering |
US20160100476A1 (en) * | 2014-10-07 | 2016-04-07 | Unist Academy-Industry Research Corporation | Method and apparatus for generating high-power terahertz wave |
CN105808504A (zh) * | 2016-03-18 | 2016-07-27 | 西安理工大学 | 一种等离子体中使用辅助微分方程完全匹配层的实现方法 |
CN105825015A (zh) * | 2016-03-18 | 2016-08-03 | 中国人民解放军火箭军工程大学 | 一种用于磁化等离子体的时域有限差分方法 |
-
2017
- 2017-03-31 CN CN201710205431.5A patent/CN107016184B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101449337A (zh) * | 2006-06-01 | 2009-06-03 | 约翰·塞尔登·奥格莱 | 利用多个相交的正离子存储环的先进核聚变方法及组件 |
US20160097722A1 (en) * | 2014-10-07 | 2016-04-07 | Unist Academy-Industry Research Corporation | Plasma diagnostic method and apparatus using raman scattering |
US20160100476A1 (en) * | 2014-10-07 | 2016-04-07 | Unist Academy-Industry Research Corporation | Method and apparatus for generating high-power terahertz wave |
CN104750990A (zh) * | 2015-03-30 | 2015-07-01 | 西安理工大学 | 二维等离子体中扩展坐标的完全匹配吸收边界的实现方法 |
CN104794289A (zh) * | 2015-04-23 | 2015-07-22 | 西安理工大学 | 一种扩展直角坐标系下完全匹配吸收边界的实现方法 |
CN104809286A (zh) * | 2015-04-23 | 2015-07-29 | 西安理工大学 | 一种等离子体中扩展坐标的完全匹配吸收边界的实现方法 |
CN104809343A (zh) * | 2015-04-23 | 2015-07-29 | 西安理工大学 | 一种等离子体中使用电流密度卷积完全匹配层的实现方法 |
CN104820660A (zh) * | 2015-04-23 | 2015-08-05 | 西安理工大学 | 一种扩展柱坐标系下完全匹配吸收边界的实现方法 |
CN105808504A (zh) * | 2016-03-18 | 2016-07-27 | 西安理工大学 | 一种等离子体中使用辅助微分方程完全匹配层的实现方法 |
CN105825015A (zh) * | 2016-03-18 | 2016-08-03 | 中国人民解放军火箭军工程大学 | 一种用于磁化等离子体的时域有限差分方法 |
Non-Patent Citations (4)
Title |
---|
HELLMUT SCHMIDT 等: ""Neutron Emission Characteristics of Pinched Dense Magnetized Plasmas"", 《IEEE TRANSACTIONS ON PLASMA SCIENCE》 * |
YUN FANG等: ""An Effective SC-PML Implementation for Iterative WLP-FDTD Method With Reduced Splitting Error"", 《IEEE MICROWAVE AND WIRELESS COMPONENTS LETTERS》 * |
刘少斌 等: ""各向异性磁化等离子体JEC_FDTD算法"", 《物理学报》 * |
殷雄 等: ""高超声速钝锥等离子体鞘套电波传播特性研究"", 《装备环境工程》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111259587A (zh) * | 2020-01-16 | 2020-06-09 | 中国人民解放军陆军工程大学 | 一种三维交替迭代无条件稳定fdtd算法 |
CN111783339A (zh) * | 2020-06-30 | 2020-10-16 | 西安理工大学 | 电磁波在随机色散介质中传播的pce-fdtd方法 |
CN111783339B (zh) * | 2020-06-30 | 2024-04-16 | 西安理工大学 | 电磁波在随机色散介质中传播的pce-fdtd方法 |
CN113887104A (zh) * | 2021-10-10 | 2022-01-04 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种基于mhd的电磁波与磁化等离子体作用数值模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107016184B (zh) | 2021-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104794289B (zh) | 一种扩展直角坐标系下完全匹配吸收边界的实现方法 | |
CN104809343B (zh) | 一种等离子体中使用电流密度卷积完全匹配层的实现方法 | |
CN105825015B (zh) | 一种用于磁化等离子体的时域有限差分方法 | |
Beghein et al. | A space-time mixed Galerkin marching-on-in-time scheme for the time-domain combined field integral equation | |
CN104750990B (zh) | 二维等离子体中扩展坐标的完全匹配吸收边界的实现方法 | |
CN107016184A (zh) | 一种二维高精度迭代的非磁化等离子体中的实现方法 | |
CN110852025B (zh) | 一种基于超收敛插值逼近的三维电磁慢扩散数值模拟方法 | |
Yue et al. | A multi‐grid method of high accuracy surface modeling and its validation | |
CN103413013A (zh) | 一种锂离子电池组系统内电池单体的原位热仿真分析方法 | |
CN104809286B (zh) | 一种等离子体中扩展坐标的完全匹配吸收边界的实现方法 | |
CN106777472A (zh) | 基于拉盖尔多项式的减少分裂误差的完全匹配层实现方法 | |
Holmes et al. | Evolution of flux-conserving tokamak equilibria with preprogrammed cross sections | |
CN105717547A (zh) | 一种各向异性介质大地电磁无网格数值模拟方法 | |
CN105808504B (zh) | 一种等离子体中使用辅助微分方程完全匹配层的实现方法 | |
CN104636553B (zh) | 微波铁氧体元器件的时域谱元仿真方法 | |
CN104820660B (zh) | 一种扩展柱坐标系下完全匹配吸收边界的实现方法 | |
Miranda-Aranguren et al. | An HLLC Riemann solver for resistive relativistic magnetohydrodynamics | |
CN106021184B (zh) | 一种六角形组件几何堆芯的时空动力学求解方法 | |
CN112307639B (zh) | 一种基于高品质算法的贝伦格完全匹配层仿真方法 | |
CN111339688B (zh) | 基于大数据并行算法求解火箭仿真模型时域方程的方法 | |
CN105512387A (zh) | 一种燃料组件精细功率重构方法 | |
CN106991232A (zh) | 一种一维高精度迭代的磁化等离子体中的实现方法 | |
Zhou et al. | An iterative factored topography-dependent eikonal solver for anisotropic media | |
Capilla et al. | A nodal collocation method for the calculation of the lambda modes of the PL equations | |
Yan et al. | Analysis of the overvoltage cooperative control strategy for the small hydropower distribution network |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |