时域谱元法高效电磁仿真算法的研究

发布时间：2020-09-11 12:07

　　 随着科学技术水平的不断提高,复杂电磁问题的研究开始面向于多尺度和多物理场问题,如何利用数值方法准确高效地分析这类问题成为计算电磁学的一大挑战。时域谱元法结合了谱方法与时域有限元法的优点,既有一定的建模灵活性,又有谱方法谱精度的特性,本文以时域谱元法为计算平台,着重研究了时域谱元法中的高效电磁仿真方法。论文的主要研究工作和成果归纳如下:本文第一部分研究了时域谱元法的基本原理,给出了分别基于麦克斯韦方程和波动方程的时域谱元法的公式的推导,进行了稳定性分析,然后介绍了时域谱元法的并行技术,给出了高效的并行方案设计,通过算例给出了算法的并行效率。并将并行时域谱元法成功地应用于密闭通信舱UWB信道的精确建模和仿真中,将仿真得到的接收信号与实验得到的结果作对比,进一步验证了该方法的正确性。本文第二部分介绍了基于Newmark-Beta格式的无条件稳定的时域谱元法,在此基础上研究了显隐式混合的时域谱元法,推导了迭代公式,在大网格区域,使用显式的中心差分格式,在精细网格处采用隐式的Newmark-Beta差分,以获得一个无条件稳定的迭代格式,在整个求解区域可以使用一个相对较大的时间步长,而不受最小剖分网格的限制,从而提高求解效率。在求解隐式部分形成的矩阵方程时,采用了两种方法,一种直接采用大型稀疏矩阵方程求解器UMFPACK,另一种是本文提出的一种迭代的Newmark-Beta方法,通过算例将两种方法进行了对比。最后将显隐式混合技术应用于色散媒质的仿真中,进一步验证该方法的计算效率和正确性。本文第三部分研究了不连续伽辽金时域谱元法的技术,基于EB为变量介绍了方程的建立过程。该方法可以在不同区域交界面的两边采用不同尺寸的网格进行离散,对非共形的空间网格进行处理,引入迎风通量,使得网格离散更具灵活性,在不需要精细处理的地方采用大网格离散,在含有精细结构的区域用小网格离散,大大减少了未知量个数。然后提出了一种改进的不连续伽辽金时域有限元方法,将连续与不连续伽辽金方法混合,旨在继承连续有限元法未知量小优点的同时又能保留DG方法块对角的特点,最后将不连续伽辽金时域谱元法和不连续伽辽金时域有限元法混合,这样既可以利用六面体剖分未知量少,又可以灵活地采用四面体网格模拟复杂物体的外形。同时引入了多时间步长技术,加快了问题的求解。本文第四部分研究了非线性微波电路的时域谱元仿真算法,建立了场-路耦合的时域谱元法同步求解技术,将时域谱元法应用于含非线性半导体器件微波电路的高效仿真中,首先通过半导体器件的等效物理模型和等效电路模型分析了微波电路的输出响应,然后引入器件的真实物理模型,求解漂移-扩散方程组,分析了含半导体器件微波电路的电热响应。最后将不连续伽辽金技术引入到半导体物理模型的求解中,采用高低阶基函数混合的方法,加速仿真求解的时间。总之,论文详细研究了基于时域谱元法的高效电磁仿真算法,针对大规模电磁计算问题,提出了高效的并行技术;针对多尺度问题,提出了显隐式混合的方法;针对复杂外形难以剖分的问题,提出了谱元法和DG有限元法的混合方法,针对非线性微波电路问题,提出了场-路耦合的同步求解技术以及基于半导体器件物理模型的电热一体化仿真方法。通过大量数值算例,验证了本文各方法的正确性与高效性。
【学位单位】：南京理工大学
【学位级别】：博士
【学位年份】：2018
【中图分类】：O441;O411
【部分图文】：

时域有限元法由于其对复杂媒质和材料灵活的建模能力，近年来成为计算电一种重要的全波仿真技术，然而由于需要在每个时间步求解一全局的矩阵方程，及到大型稀疏矩阵的求逆运算，所以计算效率不高。为了解决这一问题，国内者提出了多种质量矩阵块对角化技术，近年来比较热门的一种就是时域谱元法，采用正交的基函数，获得的质量矩阵具有对角或者块对角特性，其逆矩阵很容易此可获得待求量显式的求解表达式，极大地减少计算时间。另一方面，谱元法方法的谱精度特性，随着多项式阶数的提高，计算误差呈指数下降，因而相同度条件下，谱元法的剖分尺寸可以放大，使得需求解的未知量个数显著减少。逡逑时域谱元法的基本原理逡逑时域谱元法作为一种特殊的有限元方法，近年来逐渐成为一种求解时域微分要方法，得到了广泛的应用。谱元法采用曲六面体建模，如图２．２．１所示，因而ＴＤ而言，能更好地拟合复杂物体的外形，同时生成的质量矩阵具有块对角的此不需要像有限元一样，每个时间步求解一个大的矩阵方程，更加省时。逡逑

图２．２．１标准参量空间中的二阶ＧＬＬ矢量基函数分布图逡逑数的定义域是每条边－１到１的标准立方体单元，而仿真一的是直六面体或者曲六面体单元建模，此时需要建立一量坐标系中标准的正方体单元映射成ｘｙｚ坐标系单的直六面体单元，单元内任何一点的几何坐标满足关ｘ邋＝邋０．５［ｘ０＋ｘ，＋（＾－ｘ０）＾］逡逑－ｊ邋＝邋０．５［ｊ０＋＾，＋邋（＞；，－ｊ０）７７］ｚ邋＝邋０．５［ｚ０＋ｚ，邋＋邋（ｚ１－ｚ０）＾］逡逑ｆ，？７，４＂，点（尤。，；＾，＾）为直六面体在三个方向上的起见式（２＿２．７）代表的参量映射是一个简单的坐标缩放操作为：逡逑０．５（戋－工０）逡逑Ｊ＝逦0－５Ｕ－＾0）

总共离散出４３２个曲六面体单元，如图２．２．３所示。使用二阶基函数计算，圆柱谐振观察点的时域波形如图２．２．４所示逡逑图２．２．３圆柱谐振腔剖分网格图逡逑１５００邋逦逡逑

【相似文献】

相关期刊论文 前10条

1 杨植震;;用β反散射测量半导体器件底座的镀金层厚度[J];腐蚀与防护;1980年04期

2 晋李华;;国产塑封半导体器件可靠性试验分析[J];信息系统工程;2017年04期

3 王俊生;;半导体器件分析仪任意波形发生电路的设计[J];安徽电子信息职业技术学院学报;2017年03期

4 毕经海;;Arduino创意设计:点亮LED灯[J];少年电脑世界;2017年05期

5 刘哲;;谈常用的电子半导体器件[J];科技创业家;2012年20期

6 ;2010年全国半导体器件技术研讨会[J];微纳电子技术;2010年04期

7 唐莹;李万清;;《半导体器件》课程教学的改革探索[J];重庆电子工程职业学院学报;2010年02期

8 张永华;杨平雄;高婉丽;;“半导体器件原理”课程建设的研究[J];电气电子教学学报;2009年04期

9 任振;;基于有限元方法的半导体器件电-热分析[J];计算机仿真;2009年10期

10 余望;;半导体器件的工艺及测试分析[J];科技资讯;2008年35期

相关会议论文 前10条

1 王桂珍;林东生;杨善潮;郭晓强;李瑞宾;白小燕;刘岩;;半导体器件γ瞬时辐射的脉冲宽度效应研究[A];中国核科学技术进展报告——中国核学会2009年学术年会论文集（第一卷·第7册）[C];2009年

2 许铭真;马金源;谭长华;谭映;王洁;靳磊;;用于半导体器件参数检测的比例差值谱仪[A];第十五届全国化合物半导体材料，微波器件和光电器件学术会议论文集[C];2008年

3 邓永孝;;半导体器件的静电损伤及防护[A];电子产品静电防护技术研讨会论文集[C];1992年

4 戴精武;;微波半导体器件的陶瓷——金属封装[A];第一届陶瓷与金属封接会议论文集[C];1990年

5 陈盘训;周开明;;半导体器件的X射线辐射效应及加固[A];中国工程物理研究院科技年报（1998）[C];1998年

6 赵宇;谢劲松;;半导体器件的通用电迁移失效物理模型[A];第十一届全国可靠性物理学术讨论会论文集[C];2005年

7 母继荣;祁阳;;半导体器件抗辐射涂层研究[A];TFC'07全国薄膜技术学术研讨会论文摘要集[C];2007年

8 陈军全;陈星;;半导体器件和电路温度效应的多物理场协同计算[A];2011年全国微波毫米波会议论文集（下册）[C];2011年

9 胡新;;纳米尺度半导体器件的数值模拟[A];北京力学学会第12届学术年会论文摘要集[C];2006年

10 王礼贤;;微型焊环结构及工艺改进[A];中国机械工程学会锻压学会第六届学术年会论文集[C];1995年

相关重要报纸文章 前10条

1 本报记者 潘锋;界面是发展碳基半导体器件关键[N];科学时报;2011年

2 记者 施铭;云南半导体器件用镍铂靶材制备关键技术取得重大突破[N];云南日报;2018年

3 河南 杨远西 编译;2012年最值得关注的五类半导体器件[N];电子报;2012年

4 记者 刘其丕;新型半导体器件生产线通过验收[N];中国有色金属报;2010年

5 宋跃 王祥明;四川半导体器件出口突破1亿美元[N];中国电子报;2003年

6 左辉 李春丽;证明中国半导体器件技术实力[N];科技日报;2004年

7 邵庚;2000年半导体器件销售额创纪录[N];人民邮电;2001年

8 记者 刘其丕;天津新型半导体器件产业化项目获国家支持[N];中国有色金属报;2008年

9 杨春;我国半导体器件年销售额将达170亿美元[N];电子资讯时报;2008年

10 记者 张平阳　实习生 徐帆;新型半导体器件研发中心在我市成立[N];西安日报;2008年

相关博士学位论文 前9条

1 许浩;时域谱元法高效电磁仿真算法的研究[D];南京理工大学;2018年

2 包军林;半导体器件噪声—可靠性诊断方法研究[D];西安电子科技大学;2005年

3 闫大为;宽禁带Ⅲ族氮化物半导体器件输运与界面特性研究[D];南京大学;2011年

4 任兴荣;半导体器件的电磁损伤效应与机理研究[D];西安电子科技大学;2014年

5 李伟华;半导体器件噪声频域和时域分析的新方法研究[D];西安电子科技大学;2010年

6 赵启凤;双极型半导体器件低剂量率增强效应的1/f噪声表征研究[D];西安电子科技大学;2015年

7 范菊平;典型半导体器件的高功率微波效应研究[D];西安电子科技大学;2014年

8 汤启云;模板诱导的嵌段共聚物和纳米粒子自组装[D];南京大学;2011年

9 张丽;有机/纳米复合物场效应光电探测器研究[D];北京理工大学;2015年

相关硕士学位论文 前10条

1 沈岙;深亚微米半导体器件瞬态电热特性的谱元法分析[D];南京理工大学;2018年

2 姜昱彤;中小型半导体器件企业成本管理研究[D];沈阳工业大学;2016年

3 Adnan Ghaffar;基于MSP430的半导体器件测试设备的设计[D];兰州交通大学;2015年

4 刘影;半导体器件的二维数值模拟[D];电子科技大学;2012年

5 刘波;半导体器件的高功率微波毁伤阈值数值计算研究[D];电子科技大学;2004年

6 江代力;半导体器件的电磁热一体化建模[D];华东交通大学;2012年

7 许娜颖;进化计算在半导体器件仿真中的应用研究[D];河北工业大学;2007年

8 苟文建;半导体器件的二维仿真[D];电子科技大学;2013年

9 马云林;半导体器件交流模拟研究[D];电子科技大学;2014年

10 刘昱兵;半导体器件的1/f噪声检测方法理论与研究[D];吉林大学;2015年

本文编号：2816654