Maxwell特征值问题及Cahn-Hilliard方程的两网格算法研究及应用

来源 :湘潭大学 | 被引量 : 2次 | 上传用户:raymond20082002
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
Maxwell特征值问题的快速算法研究一直是科学和工程计算领域的重要研究课题.本文内容共分为两部分,第一部分比较系统的研究了Maxwell特征值问题的两网格算法,针对当前流行的多种形式的变分问题,提出了相应的两网格算法,证明了两网格算法的误差估计,设计了细网格上的快速算法,最后通过多个数值算例验证了理论的正确性和快速算法的高效性.第二部分针对Cahn-Hilliard方程提出了两网格算法.本文具体内容如下:首先,针对Maxwell特征值问题的棱元离散系统,设计了两网格算法.通过在粗网格求解一个原问题的特征值问题,在细网格上解一个不定方程,再由雷利熵求得近似特征值.在误差分析中,通过比较两网格方法的解与细网格上的有限元解,得到一个误差方程,并将误差分解成相互正交的两部分,其一位于特征空间,另一部分位于特征空间的正交补,我们只需要估计位于特征空间补的部分的误差并通过分析该部分的误差,利用相关引理,证得其误差估计.对于细网格上的解法器,采用了Pminres方法求解.我们采用一个正定的Maxwell方程系统作为其预条件子,该正定的系统采用HX求解.数值实验验证了理论分析的结果,并且测试了求解器的效果.接下来,我们针对混合元的离散系统,提出了相应的两网格方法.通过引入一个Lagrange乘子来处理散度为零条件,我们将要求解一个鞍点型离散系统,该系统相对庞大,但不会出现非物理解.类似地,我们在粗网格求解一个特征值问题,在细网格解一个不定方程.利用同样的证明技巧,证得其误差估计.细网格上的不定方程的求解,通过Pgmres求解.其预条件行为相当于求解一个不含低阶项的Maxwell边值问题,对于该不含低阶项的Maxwell边值问题,我们通过基于DMG磨光的多重网格算法求解,且不需要精确求解,只需要1个Vcycle.然后,我们针对面元离散系统,提出了相应的两网格方法.该离散系统所得的解将严格满足散度为零.对于细网格上的不定方程仍然采用Pgmres,其预条件行为类似于求解一个Stokes方程和一个棱元的质量矩阵.对于Stokes方程的求解,我们采用基于DGS磨光的多重网格算法[70],且不需要精确求解,只需要1-3个Vcycle.最后是本文的第二部分内容,关于Cahn-Hilliard方程的两网格算法研究.通过在粗网格求解一个混合变分问题,在细网格上求解两个第二边值条件的Poisson方程得到数值解.Poisson方程的求解采用多重网格法,数值实验证明了该算法的高效性.
其他文献
提要上脸皮肤下垂是指平视时上睑皮肤因松弛而卷垂下来,压住睫毛、使睫毛向下遮挡视线的一种疾病.严重患者上睑下垂的皮肤直接遮挡瞳孔,影响患者的学习、工作和生活.因为下垂
期刊
随着改革开放的深入,越来越多的中国公司不再满足只在中国做大做强,而是不断寻求对中国以外市场的拓展,越来越多的国内施工公司纷纷走出国门,从事国际工程承包。但与国际上成熟的
光老化表现出特征性的皮肤粗糙晦暗,点状色素沉着,毛细血管扩张和粗深皱纹,是面部皮肤哀老的主要原因。临床研究和大量临床实践已证实强脉冲光(IPL)可以改善面部光老化的各种表
目的分析吡拉西坦与尼莫地平联合使用对脑梗死后血管性痴呆患者的治疗效果以及对其认知功能的影响。方法选取200例脑梗死后血管性痴呆患者为研究对象,分成两组,各100例。对照
<正>12月28日,中国深远海海洋工程装备技术产业联盟(简称中国海工联盟)在北京成立。国家工业和信息化部副部长辛国斌、装备司司长李东,中远海运集团董事总经理万敏,中国工程
目的观察曲美他嗪治疗冠状动脉肌桥的效果.方法选择因不同程度胸闷、心悸、气短、心前区疼痛等心脏不适症状在我院心内科住院、经冠脉造影检查确诊为心肌桥的患者76 例,随机
烟幕不仅对目视形成干扰,且能对高技术的光电武器系统进行干扰,已成为对付高科技武器系统的"软胄甲"。而烟幕的构成离不开发烟剂。从常见发烟剂的种类、发烟剂的化学组成和形
<正>一、砭石简介“砭”《说文解字》注:“以石刺病也”,用于治疗的石头称为砭(bian,一声)石。人类进化的标志是学会了制造工具,而砭石正是人类最早制造的医疗工具。砭石疗法
会议
<正>突破性面部除皱技术,采用不开刀、不缝线、不填充、免拉皮一次性除皱,效果可延缓衰老5—10年,此特色技术针法操作简单,安全,效果显著,不留疤痕,是目前最为理想除皱方法。
会议
<正>2008年全国建筑模板脚手架行业年会今天在湖北省襄樊市召开了,首先让我代表中国建筑金属结构协会建筑模板脚手架委员会向来至自国各地的各位代表、各位领导表示热烈的欢