一种粗网格快速时域有限差分方法技术

技术编号:9568966 阅读:164 留言:0更新日期:2014-01-16 02:32
本发明专利技术公开了一种粗网格快速时域有限差分方法,属于电磁场数值计算领域,本发明专利技术方法的时间步长Δt与空间网格长度Δz无关,空间网格长度Δx只需小于等于模拟频段最小波长的1/2。本发明专利技术能够同时解决传统时域有限差分方法的两大限制条件,即Courant-Friedrich-Levy时间稳定性条件和空间离散间隔限制条件,能够在降低时间稳定性条件的同时,改善波长对空间网格长度的限制,本发明专利技术方法适用于模拟同时具有精细结构和电大尺寸结构的复杂目标,相比于传统时域有限差分方法,该方法具有计算效率高、计算所需内存少两大优点。

【技术实现步骤摘要】
【专利摘要】本专利技术公开了,属于电磁场数值计算领域,本专利技术方法的时间步长Δt与空间网格长度Δz无关,空间网格长度Δx只需小于等于模拟频段最小波长的1/2。本专利技术能够同时解决传统时域有限差分方法的两大限制条件,即Courant-Friedrich-Levy时间稳定性条件和空间离散间隔限制条件,能够在降低时间稳定性条件的同时,改善波长对空间网格长度的限制,本专利技术方法适用于模拟同时具有精细结构和电大尺寸结构的复杂目标,相比于传统时域有限差分方法,该方法具有计算效率高、计算所需内存少两大优点。【专利说明】
本专利技术属于电磁场数值计算领域,具体涉及。
技术介绍
目前,随着科学技术的发展,电磁波频谱的应用范围已从微波波段扩展至光波波段,器件尺寸也缩小至纳米级范围。越来越多的电磁目标呈现出同时具有精细结构和电大尺寸结构的复杂特征,如太赫兹波段下的大规模集成电路器件、大型复杂载体上的发射天线、核电磁脉冲作用下的电子信息系统以及薄层涂覆电大目标等等。这些同时具有精细结构和电大尺寸结构的复杂目标,对电磁场数值计算方法的计算速度以及计算机的内存提出了新的挑战。众所周知,时域有限差分方法是当今计算电磁学中应用最为广泛的数值模拟方法之一。然而,采用该方法来模拟同时具有精细结构和电大尺寸结构的复杂目标,却面临着很大的困难。首先,该方法需满足Courant-Friedrich-Levy时间稳定性条件:【权利要求】1.,其特征在于,包括以下步骤: 1)对待求电磁目标模型进行空间离散:磁场节点和电场节点的空间排布采用Yee元胞,电场节点Ex、Ey和Ez位于元胞的各个棱上,磁场节点Hx垂直于元胞的yz平面,磁场节点Hy与电场节点Ez的空间位置重合,磁场节点Hz与电场节点Ey的空间位置重合; 2)对待求电磁目标模型进行时间取样:电场分量的时间步取值为η和η+1时刻,磁场分量的时间步取值也为η和η+1时刻; 3)对Maxwell方程中的空间求导项從,/&、cEy/cz、clI Jcz、dll' jcz采用混合时间步法进行时间离散; 4)对得到的&:'/&、OEV'/ex, CU:/οχ , ?//:ηΙ/&求导项采用傅立叶变换求解,其余空间求导项采用二阶中心差分近似; 5)利用公式(I)和(2)求解η+1时刻的电场分邏[和磁场分量 2.根据权利要求1所述的,其特征在于,步骤3)所述的对Maxwell方程中的空间求导项故/&、dE^/cz、cl I Jcz、oil yjcz采用混合时间步法进行时间离散;所述的混合时间步法是指在对场量Ex、Ey、Hx和Hy,在时刻的选取上,包含了未知的η+1时刻。3.根据权利要求1所述的,其特征在于,时间步长 At与空间网格长度ΔΖ无关,时间稳定性条件满足: 【文档编号】G06F17/11GK103514143SQ201310433245【公开日】2014年1月15日 申请日期:2013年9月22日 优先权日:2013年9月22日 【专利技术者】陈娟, 王建国, 张安学, 田春明 申请人:西安交通大学本文档来自技高网...

【技术保护点】
一种粗网格快速时域有限差分方法,其特征在于,包括以下步骤:1)对待求电磁目标模型进行空间离散:磁场节点和电场节点的空间排布采用Yee元胞,电场节点Ex、Ey和Ez位于元胞的各个棱上,磁场节点Hx垂直于元胞的yz平面,磁场节点Hy与电场节点Ez的空间位置重合,磁场节点Hz与电场节点Ey的空间位置重合;2)对待求电磁目标模型进行时间取样:电场分量的时间步取值为n和n+1时刻,磁场分量的时间步取值也为n和n+1时刻;3)对Maxwell方程中的空间求导项采用混合时间步法进行时间离散;4)对得到的求导项采用傅立叶变换求解,其余空间求导项采用二阶中心差分近似;5)利用公式(1)和(2)求解n+1时刻的电场分量和磁场分量6)利用公式(3)求解n+1时刻的电场分量7)利用公式(4)求解n+1时刻的电场分量8)利用公式(5)和(6)求解n+1时刻的磁场分量和Hxn+1(i,j+12,k+12)=Hxn(i,j+12,k+12)+Δt2Δμz[Eyn+1(i,j+12,k+1)-Eyn+1(i,j+12,k)+Eyn(i,j+12,k+1)-Eyn(i,j+12,k)]-ΔtμΔy[Ezn+1(i,j+1,k+12)-Ezn+1(i,j,k+12)]---(5);以上各式中,表示傅立叶变换,表示逆傅立叶变换;9)令n=n+1,重复执行步骤5)~8)直至迭代完成。FDA0000385091690000011.jpg,FDA0000385091690000012.jpg,FDA0000385091690000013.jpg,FDA0000385091690000014.jpg,FDA0000385091690000015.jpg,FDA0000385091690000016.jpg,FDA0000385091690000017.jpg,FDA0000385091690000021.jpg,FDA0000385091690000022.jpg,FDA0000385091690000023.jpg,FDA0000385091690000024.jpg,FDA0000385091690000025.jpg,FDA0000385091690000031.jpg,FDA0000385091690000032.jpg,FDA0000385091690000033.jpg...

【技术特征摘要】

【专利技术属性】
技术研发人员:陈娟王建国张安学田春明
申请(专利权)人:西安交通大学
类型:发明
国别省市:

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

1