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

来源 :电子科技大学 | 被引量 : 2次 | 上传用户:jackfbi
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
计算数学的应用遍及当前科学界的各个领域。在航空航天、生命科学、资源勘探、材料设计等等方面都发挥着重要的作用。利用现代高性能的计算机,从数学理论出发,建立事物的物理模型,经过求解相应的方程,得出最后的期望的结果。在这一系列的过程中,求解线性代数方程组的计算量常常占了整个计算过程的80%以上。若线性系统的系数矩阵的谱性质又并非很好,则其迭代求解就更加困难。如何能够准确、快速的求解这些大规模的线性系统真正成了求解许多实际问题的当务之急。大规模线性系统求解的理论意义毋庸多言,实际价值更是显而易见。本文基于此,对能够加速迭代法收敛速率的预条件方法进行了大量的研究,并针对一些有计算电磁学背景的问题,针对性的构造了多种高效的预条件方法。正定线性系统的预条件迭代求解一直是许多学者研究的热点。而M-矩阵正是特殊的一类正定矩阵,根据M-矩阵一些特殊的性质有针对性的构造预条件方法更是令许多学者非常着迷。对有块三对角结构的M-矩阵,利用其块结构,进行块不完全LU分解是非常高效的预条件方法之一。为了保证块不完全分解的效率同时又保证分解因子的稀疏性,我们令分解因子保持类似ILU(k)的非零模式。同时根据M-矩阵的不完全分解方法可以构造其一种正规分裂,从而从理论上也证明了预条件子的有效性。而为了一定程度上提高其预条件子构造的可并行性,我们提出了一种重启的方法,使得预条件子兼顾效率的同时,可以节省更多的构造时间。针对二阶椭圆方程的数值实验则进一步的证明了其效率。对于用棱边元方法离散麦克斯韦方程得到的不定线性系统,根据离散矩阵的特点,利用其刚度矩阵与质量矩阵,我们给出了一类正定的预条件子。这类预条件子的构造方法十分简单,因其正定,求解也不难。针对这种正定的预条件子我们给出了其预条件之后线性系统的谱分布,并通过算例,进一步显示出其效用。矢量波动方程在研究物体散射问题时常常被用到,利用棱边元方法离散之后得到的通常是大型、稀疏、不定的复线性系统。在对这类线性系统进行不完全分解预条件方法时常常会遇到预条件效果不好的挑战,如何提高不完全分解的加速效果十分重要。通过对系数矩阵的对角元进行扰动是提高其效率的一种非常有效的手段,结合扰动技术与一种修改的不完全分解方法做为迭代求解的预条件子,对几类模型问题的实验表明结合了扰动技术的不完全分解方法对于迭代求解的加速效果十分明显。有限元与矩量法是离散麦克斯韦方程的两种重要手段。而混合有限元-矩量法更是充分利用了两种方法的各自优势。针对利用混合法离散不同的问题得到的线性系统,我们提出了几种不同的预条件方法。对于利用混合法离散介质体散射问题时,SOR预条件方法有着不俗的表现。我们也提出了几种系数矩阵的近似矩阵做为预条件子,它们的表现也十分好,同SOR相比它们求解也更加简单,而且效率也十分可观。针对离散天线问题时,根据其系数矩阵2×2块的结构,研究了块不完全分解与两层预条件方法,对比给出了最适合求解的线性系统的形式,最后通过相应的数值算例也证明了预条件技术的高效性。源于Helmholtz方程在电磁计算、声波传播、地球物理等等领域的广泛应用,其预条件迭代求解方法的研究一直备受关注。算子预条件方法不同于从系数矩阵出发的那些经典的预条件方法,它以模型的算子为根本,从物理意义出发,通过对算子进行修改,得到一个近似算子做为其预条件子。与经典的预条件方法相比,更容易从物理意义上解释其有效性。我们根据其预条件算子,利用代数多重网格方法对Helmholtz方程的迭代求解进行加速,实验表明对有阻尼的Helmholtz方程的求解,应用了提出的预条件方法之后,其迭代步数几乎不随着线性系统的规模增大而增长。
其他文献
ue*M#’#dkB4##8#”专利申请号:00109“7公开号:1278062申请日:00.06.23公开日:00.12.27申请人地址:(100084川C京市海淀区清华园申请人:清华大学发明人:隋森芳文摘:本发明属于生物技
慢性心力衰竭(慢性心衰)是临床极为常见的危重症,以利尿剂、扩血管、强心等为主要治疗手段。近年来临床研究证明血管紧张素转换酶抑制剂和β受体阻滞剂能改善患者的症状并降
ue*M#’#dkB4##8#”专利申请号:00109“7公开号:1278062申请日:00.06.23公开日:00.12.27申请人地址:(100084川C京市海淀区清华园申请人:清华大学发明人:隋森芳文摘:本发明属于生物技
在对大体积混凝土进行施工的过程中,裂缝是影响其质量的主要因素,因此施工单位应该严格控制施工中的温度。本文首先分析了大体积混凝土裂缝产生的主要原因,并对大体积混凝土
电力系统的安全稳定运行是保证可靠供电的基础,暂态稳定问题是电力系统稳定性研究的一个重要方向,研究人员一直致力于发展快速、有效的暂态稳定分析与控制方法,并取得了一定
通过2008~2009年连续2年的灯光诱测,共诱集到11目42科168种昆虫,其中林业害虫7目36科161种,占总虫量的99.31%;益虫6目6科7种,占总虫量的0.69%;诱集到的昆虫益害比为1︰144.19。
【正】贡桑诺尔布,字乐亭,1871年6月26日生于内蒙占卓索图盟喀喇沁右翼旗蒙古贵族家庭。清朝历封头等塔布囊、辅国公、喀喇沁郡王、卓索图盟盟长,民国历任蒙藏事务局(院)总裁
元代无疑是山水画发展史上的一个重要阶段,山水画的成熟不仅仅是元代美术发展中的一个重要过程,而且在文人画本身的发展过程中,元代的美术发展也具有承上启下的意义。在元代
利用混合多相流空穴模型进行了喷嘴内气液两相流动的数值模拟。分析了喷油压力、背压、喷嘴几何特征参数对喷孔内空穴流动、空穴区域分布的影响,以及各种参数对流量系数和燃油
君子在古汉语中既指有地位的统治者,又指有德之人,并且随着时代的变化其所指也相应地变化。古人常将君子比德于玉,君子有礼,敬天爱人;君子博学,义以为质;君子求诸己,自强不息