笛卡尔坐标系下各向异性完全匹配层截断边界的实现方法技术

技术编号:17780133 阅读:129 留言:0更新日期:2018-04-22 08:51
本发明专利技术涉及笛卡尔坐标系下各向异性完全匹配层截断边界的实现方法,该方法是建立在笛卡尔坐标系下的,三维问题中各向异性完全匹配层的外形为球形,二维问题中各向异性完全匹配层的外形为圆形,可以保证网格大小的一致,不会有不稳定的因素出现,又能截断一些多余的网格,减少计算量,提高计算效率,克服了现有笛卡尔坐标系下的方形各向异性完全匹配层计算效率低、计算复杂的问题,将球形的外形和方形的网格相结合,既稳定效率又高。通过

【技术实现步骤摘要】
笛卡尔坐标系下各向异性完全匹配层截断边界的实现方法
本专利技术属于计算电磁学
,具体涉及一种笛卡尔坐标系下各向异性完全匹配层截断边界的实现方法。
技术介绍
电磁场数值计算方法广泛应用于微波电路、天线设计、目标散射计算和电磁兼容等方面的研究,其中时域有限差分方法是最常用的方法之一。对于计算机编程求解电磁场问题,有限的内存是无法模拟无限大的计算区域的,因此需要在计算区域边缘设置吸收层作为截断边界,使电磁波在边界上被吸收,其性能与电磁场数值计算的精度和效率是紧密相关的。StephenD.Gendney于1996年提出了各向异性完全匹配层,可以实现较好的吸收效果,但是固有的立方体网格使得计算区域为立方体形或者正方形,在角区域和棱区域处理的过程中浪费了大量的计算机资源,并且会造成不必要的反射。专利号为201210177288.0的中国专利公开了一种二维柱坐标系完全吸收边界的实现方法,省略了棱区域的计算,但其中涉及过多繁杂的数学表达式,给计算机编程带来了不小的难度。此外,由于其计算区域的形状是圆柱形的,网格形状也是弯曲的,网格的尺寸与半径成正比,因此随着半径的增大,越外侧的网格尺寸也会越大,从而导致计算精度大幅降低,另外网格的尺寸也有可能超出“Courant稳定性条件”的限制导致算法的不稳定。申请号为201410568490.5的中国专利公开了一种阻抗匹配层的截断边界,该截断边界采用阻抗匹配形式进行截断,该过程中所涉及的参数较多,计算十分复杂,在弯曲处理截断边界时,过多的参数将导致计算机程序变得十分冗余,影响计算效率,另一方面,由于阻抗匹配层的自身缺陷,可适用的范围较窄,不能满足使用需求。然而,目前尚未出现建立在笛卡尔坐标系下弯曲形状的各项异性完全匹配层截断边界的实现方法。
技术实现思路
本专利技术的目的是提供一种笛卡尔坐标系下各向异性完全匹配层截断边界的实现方法,在笛卡尔坐标系下对传统的立方体形或者正方形各向异性完全匹配层截断边界进行弯曲处理,有效的避免了计算那些角区域和棱区域,在吸收向外传播的电磁波的同时,实现大幅减少计算量,提高计算效率。本专利技术所采用的技术方案是一种笛卡尔坐标系下各向异性完全匹配层截断边界的实现方法,该方法的步骤为:1)建立求解对象的模型数据和时域有限差分法的计算空间;向计算机申请内存空间,设Xn,Yn,Zn分别为x,y,z方向上的最小位置,其中Xn=Yn=Zn,Xp,Yp,Zp分别为x,y,z方向上的最大位置,其中Xp=Yp=Zp,总体计算的区域范围为(Xn,Yn,Zn)→(Xp,Yp,Zp),呈立方体形,区域范围的大小为(Xp–Xn)×(Yp–Yn)×(Zp–Zn),设x,y,z方向上的空间步长为Δx,Δy和Δz,且Δx=Δy=Δz,构成立方体网格区域;设球形各向异性完全匹配层的直径为a,所对应的网格数量为a/Δx=Xp–Xn,中心坐标为((Xp–Xn)/2,(Yp–Yn)/2,(Zp–Zn)/2),用球形的边界截断立方体网格区域,使计算区域成为球形,设定时间步长为Δt,并对球形内部区域的网格做inner标记,对球形外部区域的网格做outer标记,以区分不同的区域;设定电磁仿真迭代步数为N,n表示时间步,n的范围为1→N,选择时谐场的点源和电流源分别作为二维和三维问题的激励源,激励源函数用Einc表示;2)构建球形各向异性完全匹配层截断边界;在步骤1)所设置的标记为inner的球形内部区域,设定δ为球形各向异性完全匹配层截断边界的层数,将标记为inner的球形内部区域分为两部分,包括半径从0到a/2-δ*Δx的区域,称为自由空间,做inner-free标记,半径从a/2-δ*Δx到a/2的区域,即为球形各向异性完全匹配层截断边界,做inner-CSUPML标记;在球形各向异性完全匹配层截断边界的最内层,即与球心距离为a/2-δ*Δx的位置,电导率σ=0,在球形各向异性完全匹配层截断边界的最外层,即与球心距离为a/2的位置,电导率σ=σmax,中间均匀过度;σmax表示电导率σ可取的最大值,按照计算得到,其中,εr为相对介电常数,π为圆周率;标记为inner-free的自由空间中,σ是某一实数,具体数值根据内部材料的不同而不同;得到了球形的自由空间和球形各向异性完全匹配层;3)对球形内部区域的电磁场系数进行初始化;用(i,j,k)表示位置坐标,位置坐标可取的范围为步骤1)中所设定的区域范围,计算时,(i,j,k)从(Xn,Yn,Zn)点逐点循环至(Xp,Yp,Zp)点,在计算区域内对电场强度分量、电场强度中间变量、上一时刻电场强度中间变量、磁场强度分量、磁场强度中间变量、上一时刻磁场强度中间变量进行初始化操作,即均置为0,再由时域有限差分方法分别计算电场系数和磁场系数;4)更新球形内部区域的磁场强度中间变量;根据步骤3)中所设定的磁场系数,在标记为inner的球形内部区域,给每一点的上一时刻磁场强度中间变量赋值,数值等于上一时刻(即n时刻)的磁场强度中间变量,赋值完成后,计算下一时刻(即n+1时刻)的磁场强度中间变量;5)更新球形内部区域的磁场强度分量;根据步骤4)中的磁场强度中间变量在n时刻和n+1时刻的数值和步骤3)中所设定的磁场系数,在标记为inner的球形内部区域,计算出每一点在n+1时刻的磁场强度;6)更新球形内部区域的电场强度中间变量;根据步骤3)中所设定的电场系数,在标记为inner的球形内部区域,给每一点的上一时刻电场强度中间变量赋值,数值等于上一时刻(即n-1/2时刻)的电场强度中间变量,赋值完成后,计算下一时刻(即n+1/2时刻)的电场强度中间变量;7)更新球形内部区域的电场强度分量;根据步骤6)中的电场强度中间变量在n-1/2时刻和n+1/2时刻的数值和步骤3)中所设定的电场系数,在标记为inner的球形内部区域,计算出每一点在n+1/2时刻的电场强度;8)更新电场的激励源;按照公式逐步更新激励源的数值,其中,n为时间步;J0为脉冲的幅值;τ为常数;脉冲峰值出现在n=n0时刻;将更新后的激励源数值赋值给球形区域内中心位置的电场强度分量Ez,即完成电场激励源位置及数值的更新;9)每次循环时间步n都会加一,判断迭代次数n是否达到步骤1)所设定的电磁仿真迭代步数N,以判断是否满足更新条件;若未达到迭代步数,则在n加一后返回步骤4)继续循环;若达到迭代步数,即n=N,则记录步骤5)中得到的磁场强度和步骤7)中得到的电场强度,并保存作为最终所得结果。与现有技术相比,本专利技术的有益效果是:本专利技术方法是建立在笛卡尔坐标系下的,三维问题中各向异性完全匹配层的外形为球形,二维问题中各向异性完全匹配层的外形为圆形,这样可以保证网格大小的一致,不会有不稳定的因素出现,又能截断一些多余的网格,减少计算量,提高计算效率,克服了现有笛卡尔坐标系下的方形各向异性完全匹配层计算效率低、计算复杂的问题,将球形的外形和方形的网格相结合,既稳定效率又高。本专利技术各向异性完全匹配层截断边界具有表达简单、更稳定、应用更广泛等优势。本专利技术一种笛卡尔坐标系下各向异性完全匹配层截断边界的实现方法,在二维情况下,较传统的正方形各向异性完全匹配层可节省21.545%以上的内存,利用本申请所述的参数和方法设置边界参数本文档来自技高网
...
笛卡尔坐标系下各向异性完全匹配层截断边界的实现方法

【技术保护点】
一种笛卡尔坐标系下各向异性完全匹配层截断边界的实现方法,该方法的步骤为:1)建立求解对象的模型数据和时域有限差分法的计算空间;向计算机申请内存空间,设Xn,Yn,Zn分别为x,y,z方向上的最小位置,其中Xn=Yn=Zn,Xp,Yp,Zp分别为x,y,z方向上的最大位置,其中Xp=Yp=Zp,总体计算的区域范围为(Xn,Yn,Zn)→(Xp,Yp,Zp),呈立方体形,区域范围的大小为(Xp–Xn)×(Yp–Yn)×(Zp–Zn),设x,y,z方向上的空间步长为Δx,Δy和Δz,且Δx=Δy=Δz,构成立方体网格区域;设球形各向异性完全匹配层的直径为a,所对应的网格数量为a/Δx=Xp–Xn,中心坐标为((Xp–Xn)/2,(Yp–Yn)/2,(Zp–Zn)/2),用球形的边界截断立方体网格区域,使计算区域成为球形,设定时间步长为Δt,并对球形内部区域的网格做inner标记,对球形外部区域的网格做outer标记,以区分不同的区域;设定电磁仿真迭代步数为N,n表示时间步,n的范围为1→N,选择时谐场的点源和电流源分别作为二维和三维问题的激励源,激励源函数用Einc表示;2)构建球形各向异性完全匹配层截断边界;在步骤1)所设置的标记为inner的球形内部区域,设定δ为球形各向异性完全匹配层截断边界的层数,将标记为inner的球形内部区域分为两部分,包括半径从0到a/2-δ*Δx的区域,称为自由空间,做inner‑free标记,半径从a/2-δ*Δx到a/2的区域,即为球形各向异性完全匹配层截断边界,做inner‑CSUPML标记;在球形各向异性完全匹配层截断边界的最内层,即与球心距离为a/2-δ*Δx的位置,电导率σ=0,在球形各向异性完全匹配层截断边界的最外层,即与球心距离为a/2的位置,电导率σ=σmax,中间均匀过度;σmax表示电导率σ可取的最大值,按照...

【技术特征摘要】
1.一种笛卡尔坐标系下各向异性完全匹配层截断边界的实现方法,该方法的步骤为:1)建立求解对象的模型数据和时域有限差分法的计算空间;向计算机申请内存空间,设Xn,Yn,Zn分别为x,y,z方向上的最小位置,其中Xn=Yn=Zn,Xp,Yp,Zp分别为x,y,z方向上的最大位置,其中Xp=Yp=Zp,总体计算的区域范围为(Xn,Yn,Zn)→(Xp,Yp,Zp),呈立方体形,区域范围的大小为(Xp–Xn)×(Yp–Yn)×(Zp–Zn),设x,y,z方向上的空间步长为Δx,Δy和Δz,且Δx=Δy=Δz,构成立方体网格区域;设球形各向异性完全匹配层的直径为a,所对应的网格数量为a/Δx=Xp–Xn,中心坐标为((Xp–Xn)/2,(Yp–Yn)/2,(Zp–Zn)/2),用球形的边界截断立方体网格区域,使计算区域成为球形,设定时间步长为Δt,并对球形内部区域的网格做inner标记,对球形外部区域的网格做outer标记,以区分不同的区域;设定电磁仿真迭代步数为N,n表示时间步,n的范围为1→N,选择时谐场的点源和电流源分别作为二维和三维问题的激励源,激励源函数用Einc表示;2)构建球形各向异性完全匹配层截断边界;在步骤1)所设置的标记为inner的球形内部区域,设定δ为球形各向异性完全匹配层截断边界的层数,将标记为inner的球形内部区域分为两部分,包括半径从0到a/2-δ*Δx的区域,称为自由空间,做inner-free标记,半径从a/2-δ*Δx到a/2的区域,即为球形各向异性完全匹配层截断边界,做inner-CSUPML标记;在球形各向异性完全匹配层截断边界的最内层,即与球心距离为a/2-δ*Δx的位置,电导率σ=0,在球形各向异性完全匹配层截断边界的最外层,即与球心距离为a/2的位置,电导率σ=σmax,中间均匀过度;σmax表示电导率σ可取的最大值,按照计算得到,其中,εr为相对介电常数,π为圆周率;标记为inner-free的自由空间中,σ是某一实数,具体数值根据内部材料的不同而不同;得到了球形的自由空间和球形各向异性完全匹配层;3)对球形内部区域的电磁场系数进行初始化;用(i,j,k)表示位置坐标,位置坐标可取的范围为步骤1)中所设定的区域范围,计算时,(i,j,k)从(Xn,Yn,Zn)点逐点循环至(Xp,Yp,Zp)点...

【专利技术属性】
技术研发人员:郑宏兴王辂张玉贤崔文杰王蒙军李尔平
申请(专利权)人:河北工业大学
类型:发明
国别省市:天津,12

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1