收藏本站
《电子科技大学》 2010年
收藏 | 手机打开
二维码
手机客户端打开本文

预条件算法及在电磁场数值模拟中的应用

任志刚  
【摘要】:计算数学的应用遍及当前科学界的各个领域。在航空航天、生命科学、资源勘探、材料设计等等方面都发挥着重要的作用。利用现代高性能的计算机,从数学理论出发,建立事物的物理模型,经过求解相应的方程,得出最后的期望的结果。在这一系列的过程中,求解线性代数方程组的计算量常常占了整个计算过程的80%以上。若线性系统的系数矩阵的谱性质又并非很好,则其迭代求解就更加困难。如何能够准确、快速的求解这些大规模的线性系统真正成了求解许多实际问题的当务之急。大规模线性系统求解的理论意义毋庸多言,实际价值更是显而易见。本文基于此,对能够加速迭代法收敛速率的预条件方法进行了大量的研究,并针对一些有计算电磁学背景的问题,针对性的构造了多种高效的预条件方法。 正定线性系统的预条件迭代求解一直是许多学者研究的热点。而M-矩阵正是特殊的一类正定矩阵,根据M-矩阵一些特殊的性质有针对性的构造预条件方法更是令许多学者非常着迷。对有块三对角结构的M-矩阵,利用其块结构,进行块不完全LU分解是非常高效的预条件方法之一。为了保证块不完全分解的效率同时又保证分解因子的稀疏性,我们令分解因子保持类似ILU(k)的非零模式。同时根据M-矩阵的不完全分解方法可以构造其一种正规分裂,从而从理论上也证明了预条件子的有效性。而为了一定程度上提高其预条件子构造的可并行性,我们提出了一种重启的方法,使得预条件子兼顾效率的同时,可以节省更多的构造时间。针对二阶椭圆方程的数值实验则进一步的证明了其效率。 对于用棱边元方法离散麦克斯韦方程得到的不定线性系统,根据离散矩阵的特点,利用其刚度矩阵与质量矩阵,我们给出了一类正定的预条件子。这类预条件子的构造方法十分简单,因其正定,求解也不难。针对这种正定的预条件子我们给出了其预条件之后线性系统的谱分布,并通过算例,进一步显示出其效用。 矢量波动方程在研究物体散射问题时常常被用到,利用棱边元方法离散之后得到的通常是大型、稀疏、不定的复线性系统。在对这类线性系统进行不完全分解预条件方法时常常会遇到预条件效果不好的挑战,如何提高不完全分解的加速效果十分重要。通过对系数矩阵的对角元进行扰动是提高其效率的一种非常有效的手段,结合扰动技术与一种修改的不完全分解方法做为迭代求解的预条件子,对几类模型问题的实验表明结合了扰动技术的不完全分解方法对于迭代求解的加速效果十分明显。 有限元与矩量法是离散麦克斯韦方程的两种重要手段。而混合有限元-矩量法更是充分利用了两种方法的各自优势。针对利用混合法离散不同的问题得到的线性系统,我们提出了几种不同的预条件方法。对于利用混合法离散介质体散射问题时,SOR预条件方法有着不俗的表现。我们也提出了几种系数矩阵的近似矩阵做为预条件子,它们的表现也十分好,同SOR相比它们求解也更加简单,而且效率也十分可观。针对离散天线问题时,根据其系数矩阵2×2块的结构,研究了块不完全分解与两层预条件方法,对比给出了最适合求解的线性系统的形式,最后通过相应的数值算例也证明了预条件技术的高效性。 源于Helmholtz方程在电磁计算、声波传播、地球物理等等领域的广泛应用,其预条件迭代求解方法的研究一直备受关注。算子预条件方法不同于从系数矩阵出发的那些经典的预条件方法,它以模型的算子为根本,从物理意义出发,通过对算子进行修改,得到一个近似算子做为其预条件子。与经典的预条件方法相比,更容易从物理意义上解释其有效性。我们根据其预条件算子,利用代数多重网格方法对Helmholtz方程的迭代求解进行加速,实验表明对有阻尼的Helmholtz方程的求解,应用了提出的预条件方法之后,其迭代步数几乎不随着线性系统的规模增大而增长。
【学位授予单位】:电子科技大学
【学位级别】:博士
【学位授予年份】:2010
【分类号】:O441.4

【共引文献】
中国期刊全文数据库 前10条
1 薛秋芳,畅大为,赵广意;相容次序矩阵AOR迭代的最优参数选取[J];安徽大学学报(自然科学版);2004年05期
2 赵佳;吴先良;程媛;;MLC-FDTD在波导带通滤波器设计中的应用[J];安徽大学学报(自然科学版);2008年02期
3 蔡二平 ,蒋全兴;用时域有限差分法模拟千兆赫兹横电磁波室的馈电接头[J];安全与电磁兼容;2004年06期
4 于大鹏;宋航;;PCB的FDTD-PSPICE混合分析方法研究[J];安全与电磁兼容;2010年05期
5 吴琼;;电磁带隙结构的优化设计[J];安庆师范学院学报(自然科学版);2006年04期
6 刘瑜;廖同庆;;现代计算环境下的电磁模拟教学浅议[J];安庆师范学院学报(自然科学版);2011年03期
7 汤炳书;沈廷根;;二维光子晶体透射特性的FDTD数值研究[J];半导体光电;2006年01期
8 李锦屏;杜丽霞;陈伟;;矩形变形脊波导主模截止波长的计算[J];半导体光电;2008年05期
9 李锦屏;杜丽霞;;脊槽波导主模截止波长研究[J];半导体光电;2011年01期
10 武晓海,殷莉,洪先龙;基于不完全分解预优共轭梯度法的电源和地线网络求解器[J];半导体学报;2000年03期
中国重要会议论文全文数据库 前10条
1 吴琼;;人工神经网络在基片集成波导滤波器设计中的应用[A];安徽省第五届“兴皖之光”青年学术年会论文集(工科卷)[C];2005年
2 董国华;贺汉根;胡德文;;非负矩阵谱半径一个严格不等式及其概率证法[A];第二十七届中国控制会议论文集[C];2008年
3 肖映雄;陈鹏;舒适;;两类网格结构模型的预处理方法[A];中国计算力学大会'2010(CCCM2010)暨第八届南方计算力学学术会议(SCCM8)论文集[C];2010年
4 王春;;一种新型宽带HF单极子天线[A];2005年海峡两岸三地无线科技学术会论文集[C];2005年
5 彭朕;盛新庆;;介质填充波导高次模求解的快速边缘元算法[A];2005年海峡两岸三地无线科技学术会论文集[C];2005年
6 郑宏兴;郁道银;;添加各向异性参数修正ADI-FDTD算法中的数值色散[A];2005年海峡两岸三地无线科技学术会论文集[C];2005年
7 常海峰;徐筱麟;鲍峻松;;电磁脉冲对孔缝的耦合效应及防护方法的研究[A];2007北京地区高校研究生学术交流会通信与信息技术会议论文集(上册)[C];2008年
8 高艳;王海锋;;缺陷一维光子晶体传输特性的研究[A];第七届中国功能材料及其应用学术会议论文集(第2分册)[C];2010年
9 钱若军;;网壳结构的稳定分析策略[A];第二届全国现代结构工程学术研讨会论文集[C];2002年
10 曾繁春;陈加清;王盼;;冲击泄流时接地体附近跨步电压的FDTD分析[A];中国通信学会第六届学术年会论文集(上)[C];2009年
中国博士学位论文全文数据库 前10条
1 刘少刚;基于单一平面电容传感器的木材含水率检测系统研究[D];哈尔滨工程大学;2009年
2 易煦农;人工电磁材料的光学特性及应用研究[D];华中科技大学;2010年
3 鲁晶津;地球电磁三维数值模拟的多重网格方法及其应用研究[D];中国科学技术大学;2010年
4 李学识;基于表面波驱动的等离子体天线及其阵列技术研究[D];华南理工大学;2010年
5 芮锡;旋转结构体电磁问题的高效分析与应用[D];电子科技大学;2010年
6 廖臣;三维电磁粒子模拟并行算法及其应用研究[D];电子科技大学;2010年
7 荆燕飞;线性方程组迭代法与预条件技术及在电磁散射计算中的应用[D];电子科技大学;2010年
8 赖生建;计算电磁学中的径向基无网格法[D];电子科技大学;2010年
9 黄卓红;PDE离散方程组和鞍点问题的预处理方法[D];电子科技大学;2010年
10 郑罡;频域有限差分方法的改进研究[D];电子科技大学;2010年
中国硕士学位论文全文数据库 前10条
1 余波;惯性约束聚变的中子半影成像诊断系统和非线性重建方法研究[D];中国工程物理研究院;2010年
2 王瑞英;H(curl)空间的各向异性分析[D];郑州大学;2010年
3 王娅娜;平面弹性问题协调的矩形混合元[D];郑州大学;2010年
4 刘兴鹏;基于微波同轴谐振腔物料密度测量技术的研究[D];哈尔滨工程大学;2010年
5 王化宇;基于时域有限差分法的微带结构电磁兼容研究[D];哈尔滨工程大学;2010年
6 刘肖峰;基于MOM的天线耦合研究[D];哈尔滨工程大学;2010年
7 张宇飞;加稀疏约束的非负矩阵分解[D];大连理工大学;2010年
8 马欢;对偶内点法求解随机线性二次最优控制问题[D];大连理工大学;2010年
9 牛炜;约束矩阵方程AX=B的迭代解法及其最佳逼近[D];长沙理工大学;2010年
10 田静;约束矩阵方程及迭代解法的预处理技术[D];长沙理工大学;2010年
【相似文献】
中国期刊全文数据库 前10条
1 李晓梅;吴建平;;稀疏线性方程组不完全分解预条件方法[J];计算机工程与科学;2006年08期
2 吴建平;宋君强;李晓梅;;块三对角线性方程组不完全分解预条件的一种一维区域分解并行化方法[J];计算物理;2008年06期
3 吴建平;赵军;宋君强;张卫民;李晓梅;;非静力模式GRAPES的预条件技术研究[J];计算机工程与应用;2011年10期
4 陈泽军;肖宏;;预条件GMRES(m)法迭代求解大规模边界元弹性问题[J];固体力学学报;2006年S1期
5 吴建平;王正华;李晓梅;;二维三温能量方程组的预条件迭代软件包研制——离散所得稀疏线性方程组的求解[J];计算机工程与应用;2007年33期
6 项铁铭,梁昌洪;计算电磁学中稠密线性方程组的迭代求解[J];西安电子科技大学学报;2003年06期
7 吴建平;宋君强;张卫民;李晓梅;;块三对角线性方程组的一类二维区域分解并行不完全分解预条件[J];计算物理;2009年02期
8 武瑞婵;;局部正交化在大规模并行计算中的应用[J];襄樊学院学报;2009年11期
9 谭林;江军;舒适;;一种三角形网格下保对称有限体元方法的预条件技术[J];湘潭大学自然科学学报;2006年01期
10 阙肖峰;聂在平;胡俊;;混合场积分方程结合MLFMA分析导体介质复合目标电磁散射问题[J];电子学报;2007年11期
中国重要会议论文全文数据库 前10条
1 樊振宏;朱剑;平学伟;陈如山;;有限元快速多极子混合方法中的求解技术[A];2007年全国微波毫米波会议论文集(上册)[C];2007年
2 郭美珍;张瑗;;基于FEM格式的二维三温辐射热传导方程组的并行预条件子[A];第13届中国系统仿真技术及其应用学术年会论文集[C];2011年
3 肖映雄;陈鹏;舒适;;两类网格结构模型的预处理方法[A];中国计算力学大会'2010(CCCM2010)暨第八届南方计算力学学术会议(SCCM8)论文集[C];2010年
4 吴建平;李晓梅;;块三对角线性方程组不完全分解预条件的一种并行化方法[A];全国计算物理学会第六届年会和学术交流会论文摘要集[C];2007年
5 安翔;吕志清;;迭代子结构法在三维电大散射问题中的应用[A];2007年全国微波毫米波会议论文集(上册)[C];2007年
6 吴梦喜;余进;;岩土工程有限元快速算法的研究与应用[A];中国力学学会学术大会'2009论文摘要集[C];2009年
7 陈璞;肖梃松;孙树立;袁明武;;预条件共轭梯度法的实现以及一些改进[A];第七届全国结构工程学术会议论文集(第Ⅰ卷)[C];1998年
8 梁峰;钱若军;;预优共轭梯度法在结构有限元分析中的应用[A];第十届空间结构学术会议论文集[C];2002年
9 聂存云;舒适;杭旭登;成娟;;柱坐标系下辐射热传导方程的SFVE格式[A];中国核科学技术进展报告——中国核学会2009年学术年会论文集(第一卷·第6册)[C];2009年
10 郑亮;Michael R.Elgersma;David A.Yuen;张怀;石耀霖;;预条件最小残差法GPU求解器用于含有泡状体的Stokes流体问题[A];中国地球物理学会第二十七届年会论文集[C];2011年
中国重要报纸全文数据库 前10条
1 记者 陈超;日发现风湿性关节炎致病机理[N];科技日报;2006年
2 孙立华;新世纪特种枪(上)[N];科技日报;2001年
3 阿娜;烹调不当 营养打折[N];中国医药报;2005年
4 李大刚;常用生石灰调水好处多[N];中国渔业报;2006年
5 阿娜;什么烹调方式营养流失少[N];大众科技报;2006年
6 杨大胜;浅谈非致命性武器的应用[N];人民公安报;2006年
7 王富生;蔬菜泥炭育苗营养块及其育苗新方法[N];瓜果蔬菜报.农业信息周刊;2008年
8 尤春涛;美国造的“卡拉什尼科夫”枪[N];中国国防报;2002年
9 肖培弘;特种水产品的越冬管理[N];中国渔业报;2004年
10 杨子;当心营养从厨房溜走[N];中国医药报;2005年
中国博士学位论文全文数据库 前10条
1 任志刚;预条件算法及在电磁场数值模拟中的应用[D];电子科技大学;2010年
2 颜溯;基于Calderón技术的计算电磁学积分方程方法研究[D];电子科技大学;2011年
3 李良;大型线性方程组求解技术及在计算电磁学中的应用研究[D];电子科技大学;2009年
4 荆燕飞;线性方程组迭代法与预条件技术及在电磁散射计算中的应用[D];电子科技大学;2010年
5 张勇;线性方程组预条件技术及在二维三温问题中的应用和实现[D];电子科技大学;2008年
6 王转德;迭代矩阵的谱分析[D];电子科技大学;2009年
7 寇继生;JFNK方法的若干改进及其在二维河道水流数值模拟中的应用[D];武汉大学;2007年
8 朱剑;复杂电磁问题的有限元、边界积分及混合算法的快速分析技术[D];南京理工大学;2011年
9 杭旭登;偏微分方程迭代并行解法与网格优化方法[D];中国工程物理研究院;2004年
10 江军;二维三温辐射热传导方程的高效自适应算法研究[D];湘潭大学;2006年
中国硕士学位论文全文数据库 前10条
1 张传文;新预条件下矩阵的收敛性分析及其比较[D];扬州大学;2010年
2 李凯;基于新预条件子作用下向后迭代法的收敛性分析[D];扬州大学;2010年
3 丁昌华;ILUT和最小度算法在大型线性方程组求解中的应用研究[D];电子科技大学;2012年
4 李乐波;大规模稀疏线性方程组的预条件迭代法的研究[D];南昌大学;2011年
5 郑振裥;多水平ILU分解及在电磁计算中的应用研究[D];电子科技大学;2012年
6 邓军;基于边界元法的变电站内工频电场计算方法研究[D];重庆大学;2010年
7 王学忠;H矩阵方程组的预条件迭代法和预条件对角占优性[D];电子科技大学;2007年
8 陈丹丹;求解Helmholtz方程的代数多重网格预条件技术的研究[D];电子科技大学;2012年
9 蔡存朋;结构模态重分析的预条件SRQCG方法[D];吉林大学;2009年
10 陈涌频;复杂目标电磁散射的快速非均匀平面波算法[D];电子科技大学;2006年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62791813
  • 010-62985026