非线性方程间断有限元方法的误差估计和保界格式

来源 :中国科学技术大学 | 被引量 : 1次 | 上传用户:xiaowei_0315
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本论文主要研究有界区域中非线性偏微分方程的间断有限元方法。我们首先证明了 Allen-Cahn方程和Cahn-Hilliard方程局部间断有限元方法的能量稳定性和最优误差估计。其次,基于Karush-Kuhn-Tucker(KKT)限制器,我们通过拉格朗日乘子分别构造了反应欧拉方程和非线性退化抛物方程高阶保界的隐式时间离散的间断有限元和局部间断有限元格式。论文的第一部分,我们研究了 Allen-Cahn方程二阶和三阶半隐式谱延迟校正(SDC)时间离散的局部间断有限元格式的能量稳定性和最优误差估计。由于SDC方法是基于一阶凸分裂格式,因此时间离散方法对非线性项的隐式处理会导致每个时间层的方程组都是非线性的,增加了理论分析的难度。对于结合二阶和三阶SDC方法的局部间断有限元离散格式,我们利用有限维空间中的不动点定理证明了数值解的存在唯一性。同时半隐式SDC格式中所涉及的迭代和积分也增加了理论分析的难度。与不包括最左端点的龙格库塔型半隐式格式相比,这里的SDC格式将最左端点作为正交节点。这使得SDC格式测试函数的选取更加复杂,能量方程的构建更加困难。我们提供了两种不同的方法来克服非线性项带来的困难。通过仔细选择测试函数,在时间步长τ仅需要一个正的上限并且与网格大小h无关的意义下,我们得到了二阶和三阶数值格式的能量稳定性和最优误差估计。数值算例验证了我们理论结果的正确性。论文的第二部分,我们主要研究了具有浓度相关迁移率的Cahn-Hilliard方程的一个无条件稳定的局部间断有限元格式的误差分析。我们使用的时间离散是基于不变能量正交化(IEQ)方法,因此我们的全离散格式在每个时间步都是一个线性代数系统。这里误差估计的主要困难是在局部间断有限元格式中缺少对单元边界上某些跳跃项的控制。我们需要特殊处理Cahn-Hilliard方程的初始条件和非恒定迁移率项。对于初始条件的误差估计,我们用一个等价光滑的全局Lipschitz连续函数代替非线性项。这种技巧仅用于初值问题。对于非恒定迁移率项的分析,我们充分利用半隐式时间离散方法的优势,通过数学归纳法得到某些数值变量在L∞-范数下的有界性。我们得到了全离散格式的最优误差估计并给出了数值算例验证此结论。论文的第三部分,我们构造了反应欧拉方程高阶保界的隐式时间离散的间断有限元格式。在反应问题中,由于流体动力时间尺度与反应时间尺度存在较大差异,所以数值计算中的时间步长往往会受到很大的限制。此外,反应问题中的密度和压强都是非负的,质量分数应该在0到1之间。这里我们用分步法分别处理对流问题和反应问题。关于反应欧拉方程,我们主要有三个贡献。首先,我们采用高阶对角隐式龙格库塔(DIRK)方法进行时间离散。与显式时间离散方法相比,隐式方法大大增加了数值计算中具有刚性源项方程的时间步长。其次,在KKT系统的基础上,我们利用拉格朗日乘子将隐式时间离散的数值离散格式与保界约束条件相结合,从而保持数值解的上界0和下界1。最后,由于刚性源项,我们将Harten的子单元分辨技术(SR)推广到反应问题隐式时间离散的间断有限元方法中。数值结果表明,保界DIRK间断有限元格式对于光滑解是高阶精度的,对不连续刚性问题的数值模拟在相对粗的网格中是相当有效的。论文的第四部分,我们针对非线性退化抛物方程提出了一个熵耗散的高阶DIRK局部间断有限元格式。对于非线性抛物方程的一些问题,目前已证明当时间趋于无穷大时,瞬态解会收敛到稳定状态。我们利用简单的交替数值流通量,构造了具有高阶、熵耗散、保稳态和可捕捉长时间行为等优点的隐式DIRK局部间断有限元格式。隐式时间离散方法大大增加了数值格式稳定性所需的时间步长。这里较大的时间步长和简单的交替数值通量极大地简化了数值计算。我们从理论上证明了半离散格式的熵耗散性及一阶全离散格式的熵耗散性和稳态保持性。为了保证数值解的正定性和质量守恒性,我们采用了 KKT限制器,通过拉格朗日乘子将正定不等式约束和质量守恒等式约束与高阶DIRK局部间断有限元格式相耦合。数值结果表明,保正的DIRK局部间断有限元格式具有较高的精度且是高效的。
其他文献
运动是动物适应环境,能够趋利避害,成功地在自然中生存的基本能力。作为动物界中最为广泛的运动形式之一,节律运动是动物执行多样的运动模式的基石。诸如昆虫飞行时翅膀的振动,哺乳动物的呼吸、心跳等都是节律运动的体现。深入理解动物以何种方式产生协调的节律运动,是探索神经控制基本原理的主要路径,可以帮助我们增加对神经系统的认识。动物的节律运动涉及神经至肌肉的驱动力发生,感觉信息反馈以及高层神经环路对运动系统的
参量下转换纠缠光源在线性光学系统中有着不可替代的作用,其可以被用来级联实现大规模的光量子线路。提高纠缠源品质一直是一个亟待解决的问题。之前的参量下转换光源收集效率低的主要原因是由于光子对的频率关联无法解除,所以需要使用窄带滤波片来解除频率所致。本文主要是通过优化参量下转换晶体的参数,找到一个无频率关联在收集效率和亮度都达到最优的参数,为了能够将其应用到实际通信中,我们将参量光波长控制在光纤损耗最小
在标准模型Higgs粒子确认发现之后,当前的高能物理界的发展目标主要集中在新物理的寻找。而新物理的寻找主要有两种途径,1.标准模型的精确检验,即从理论和实验上给出更加精确的结果,以期望找到当前的理论与实验之间的细微偏差;2.给出当前各种新物理模型的存活参数空间,在限制条件允许的情况下给出新物理可能出现的现象学,为实验寻找新物理特征信号提供理论依据与指导。本论文课题就是围绕这两条途径,分别进行了三规
获得1997年诺贝尔奖的激光冷却技术打开了研究冷原子物理的大门。近二十年来冷原子物理与技术快速发展,2001年,2005年,2012年的诺贝尔奖即是冷原子相关技术的发展应用。冷原子可以实现量子存储用于量子通信和量子计算;超冷原子则可以用于量子模拟去研究凝聚态体系的演化现象。磁场会引起原子的能级劈裂,它和激光以及射频,是调控冷原子的重要手段。磁场通过在线圈中通电流产生,实现各种磁场波形需要模拟可控的
实践是加深学生对知识的理解,提高学生能力的唯一途径。但是很多教师认为实践只能运用于带有实验的学科,其实这样的观念是错误的。在其他科目的教学中开展实践同样能发挥作用,使学生更好地成长与发展。根据科学研究表明,在习作单元的教学中开展言语实践,可以加强学生的口语表达能力,可以让学生积累一定的语言、情感、形式,让学生的语文实践能力得到有效的提高,所以教师要将其重视起来。
期刊
带状流是具有一维时空结构的流体速度场。漂移波是非均匀媒质中由压强梯度自由能驱动的高模数的波。人们普遍接受二者是所谓的漂移波-带状流系统的两个基本成分。特别的,托卡马克中的带状流是指由带状的静电势涨落的径向变化产生的EXB速度;轴对称托卡马克中的漂移波是与某个有理面相关的在径向的定域的驻波,在极向的传播是行波。尽管带状流在介观尺度,漂移波中存在两种不同的尺度:微观尺度的线性本征模和介观尺度的漂移波波
最近二十年间,拓扑理论在凝聚态物理方面发挥着越来越重要的作用,拓扑型激发在各种不同的凝聚态系统中都有体现,这些新奇的拓扑量子态极大地丰富了凝聚态物理学。陆续发现和报道了反常量子霍尔效应、拓扑绝缘体、外尔半金属等奇异的物理现象和新材料。此外,还有一种新型的拓扑态,其中包含的准粒子在交换操作下表现出非阿贝尔的统计性质,被认为有希望在未来实现量子计算。总而言之,这些与拓扑序和拓扑态相关的新材料和新现象将
本文主要研究非K?hler流形上的Higgs丛和平坦丛.目的是希望将Corlette-Donaldson-Hitchin-Simpson所建立的陈数消失半稳定(多稳定)Higgs丛和(半单)平坦丛两个模空间之间的对应关系推广到一类非K?hler流形上.在本文中我们总是假设Hermitian底流形(X,ω)满足?(?)ωn-1=0和?(?)ωn-2=0.在本文的第一部分,我们研究非K?hler流形上
第一性原理计算在物理化学数值计算方法中占有重要的地位,近些年来其计算精度与计算速度都取得了很大的进展,也因此在更多领域被广泛应用。在这其中,杂化泛函计算是重要的一类计算,对很多重要的物理性质能得出比较准确的结果。但相比于传统密度泛函理论,杂化泛函的计算量非常大,计算耗时远高于传统密度泛函理论,这制约了其在大体系计算下的应用。为了加速杂化泛函计算,实现其大规模计算,我们在单位元分解(Resoluti