Sobolev型方程的局部间断Galerkin方法

来源 :南京大学 | 被引量 : 0次 | 上传用户:fairycx
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
Sobolev型方程在数学、物理及工程等领域中应用广泛,相关的数值方法研究一直颇受重视。此类方程包含时空混合导数,使得相应的格式设计和理论分析变得较为困难。本文拟建立Sobolev型方程的局部间断有限元(LDG)方法,并从数值理论分析和数值实验两个角度阐释算法的有效性。该方法的主要亮点是辅助变量时空转换关系的充分应用,使得相应的时间推进简便易行。本文共分五章。第一章回顾Sobolev型方程及LDG方法的研究概况,最后一章是总结和展望。第二章到第四章是文章的主体部分,具体内容如下:第二章采用广义数值通量,构造二维线性Sobolev方程的LDG格式,并配合三阶全变差不增的显式Runge-Kutta(TVDRK3)时间离散,给出半离散与全离散格式的稳定性与最优误差估计。该方法巧妙运用辅助变量的时空转换关系,将高阶的发展型方程转化为一阶的发展型方程组。该做法具有两个数值优势:一是数值解的时间演化可以解耦为自由度的直接推进;二是数值通量的定义变得清晰和自然。本章给出了两种初值定义方式,并指出它们的差异:当原始变量和辅助变量在初始时刻满足DG离散关系时,LDG格式关于μ-范数是强稳定的;否则相应格式仅具有普通意义的稳定性。本章借助二维广义Gauss-Radau(GGR)投影及其修正技术,协调了不同变量在单元边界处数值通量的耦合关系,进而获得与实际相符的最佳估计。着重指出,误差估计中的界定常数与高阶项系数的倒数无关,换言之,本章结果同样适用于对流占优的Sobolev型方程。第三章将LDG方法应用于MBL方程。该方程为一类典型的非线性Sobolev方程,含有非凸对流项且为对流占优的,相应问题的解往往呈现出大梯度变化的瞬变层。本章旨在解决此类瞬变层问题数值模拟的准确性。实验表明,当方程粘性系数ε不是很小时,LDG方法可以准确捕捉到行波解的瞬变层。但在ε逐步变小的过程中,瞬变层附近的数值解将变得不再理想,出现剧烈的数值振荡。为此,我们引入三种非线性限制器:TVB限制器、WENO限制器和Moe限制器来改善数值效果。实验结果显示,选取适当参数的Moe限制器更为有效,它在减少数值振荡的同时保证了数值给出的上溢出平台高度与理论行波解相吻合。为捕捉瞬变层结构,采用全局加密将在光滑区域引起不必要的计算。为此,本章末还考虑一种自适应移动网格方法,使其在瞬变层附近进行局部加密。数值实验表明:在同样数值精度的要求下,移动网格方法有效较少了空间自由度,提升了计算效率。第四章将前两章的LDG方法推广至其它两类方程。一类是线性化的耦合BBM系统,另一类是含有五阶时空混合导数的KDV方程,它们具有或可转化为Sobolev型方程的形式。第一节我们结合TVDRK3时间离散方法,构造耦合BBM系统的半离散与全离散LDG格式。利用广义交替型通量关于DG空间离散的反对称性,可以证明,相应格式关于L~2-范数是强稳定的。同时,借助一维GGR投影技术,我们获得了L~2-范数下的最优误差估计。第二节研究LDG方法求解五阶KDV方程的理论分析。其中,格式设计的核心是引入适当的辅助变量和时间求导,将五阶方程转化为Sobolev方程组,进而给出方程组的半离散LDG格式。若数值初值满足DG空间离散关系,该格式具有与前一节类似的稳定性。在此基础上,利用一维GGR投影及其修正投影,我们证得了次最优的L~2-范数收敛阶。
其他文献
疾病是机体在一定病因作用下自稳调节絮乱而发生的异常生命活动过程。从分子生物学角度来看,疾病通常是由于蛋白质质和量的改变导致的细胞功能絮乱,最终引起的人体某个或某些器官的异常工作。在分子生物学中,与疾病相关的生物表征、关联性问题等研究在疾病分析、诊断和疾病治疗等方面都发挥了重要作用。基于分子生物学实验的方法能够得到较为准确且置信度高的实验结果,但实验成本高、效率较低且无法验证所有可能的样本。而基于计
背景:T细胞过度活化和克隆扩增是系统性红斑狼疮(SLE)发病进程中常见的免疫失衡现象,在自身免疫反应中扮演着重要角色。过度活化T细胞衰老可能是SLE发病进程中抑制免疫反应的一个重要机制。人脐带间充质干细胞(hUC-MSCs)移植治疗SLE取得较好疗效,但其对T细胞衰老的作用及机制尚不清楚。目的:本课题旨在研究hUC-MSCs在治疗SLE中对CD4+T细胞衰老的具体作用,深入探讨hUC-MSCs对C
在特定的有机半导体中,单线激发态的能量大约是三线激发态的两倍,单个单线态激子可以转换成一对三线态激子,也即单线态裂分过程。单线态裂分实现了光电转换量子效率的倍增,原则上能将单节太阳能电池的光电转换能量效率的理论极限从33%提高到45%。单线态裂分已被成功引入原型器件结构中实现了特定波段光电转换的效率提升。为进一步实现实用化,设计和合成合适的单线态裂分材料,我们亟需理解单线态裂分的物理机制。一般认为
新文化史学是批判与质疑新史学而产生的史学潮流,它注重文化的解读与阐释。新文化史学的研究理路为教育史研究开辟了新的路径并提供了诸多启示。教育史学研究的转向路径为:研究对象下移到普通民众;研究领域拓展得更为广域丰富;研究方法注重跨学科采借;表达方式转为生动形象的叙事;史料来源变为多元化搜集。这些路径转向也是教育史研究摆脱自身困境的应对之措。
锕系元素的成键和反应性研究对未来核能的发展、核废料的循环利用及理解f区元素的成键本质至关重要。贫铀(238U)是一种稳定的锕系金属核素。因此,化学家常以金属铀配合物为模型,来探究锕系元素的基本化学性质。进入21世纪以来,金属有机铀化学实现了快速发展,一系列铀化合物在催化、小分子活化、单分子磁体等方面展现出独特性能。众所周知,配体对金属有机化学的发展至关重要。人们发现基于不同配体的金属铀配合物具有不
肠外营养(Parenteral nutrition,PN)是因各种疾病或病理因素导致的胃肠道功能障碍病人的首要营养支持治疗方式,挽救了数以万计病人的生命。PN长期实施可导致糖代谢紊乱,主要表现是血糖稳定下降,高血糖和低血糖发生率均显著增加,严重影响病人的预后。传统观念认为PN时高血糖和低血糖分别是由于机械性输注含糖物质和胰岛素不当使用、PN突然停用所致。然而,PN糖代谢紊乱缺乏进一步的机制探索。肠
全球性同期造山带是连接超大陆中各陆块的重要纽带,其造山带信息可经剥蚀-沉积系统保存于这些陆块的陆源沉积记录中。从源汇关系新视角,研究华南地块自形成后最早一套沉积记录——南华系-奥陶系,不仅可丰富地质资料、澄清古地理认识争议,而且可通过辨析造山带信息,有助于深入理解华南地块在同期全球构造体系中的位置。基于野外广泛调查及获得的来自华南地块南华系至奥陶系中部新素材,本文进行了锆石U-Pb-Hf同位素及几
北极海冰是全球气候系统的重要组成部分,海冰体积直接与能量平衡有关,是气候变化的敏感指示器。由于海冰厚度数据限制,缺乏对长时间序列海冰体积变化分析。本文利用多源海冰厚度数据,通过数据评估、比较、回归分析和融合获得1979-2018年海冰厚度,利用厚度和面积数据计算近40年体积的时空变化特征。影响海冰变化的因素包括动力学(变形和输出)和热力学(消融和生长),哪个因素对体积变化的影响较大,及不同时期影响
研究和开发具有各种功能的新型材料是解决能源危机、环境问题和信息需求的关键所在。晶体结构作为凝聚态物理学中最基本的一个概念,是对晶体材料的数学描述,它决定了晶体材料几乎所有的物理化学性质,是理论上对材料物性进行预测和性能调控方案探测的基础。随着计算机软硬件能力的提升和凝聚态物理理论的发展,特别是大型高性能计算集群和基于密度泛函理论的第一性原理计算方法的出现,人们已经开始有能力根据实验结果解析出材料的
α,β-不饱和醛选择性加氢制备的α,β-不饱和醇作为一种重要的化工原料和医药中间体具有广泛的应用和研究价值。α,β-不饱和醛由于其内部存在的共轭双键C=O和C=C,其选择性加氢一直以来是研究的热点。目前对于α,β-不饱和醛的研究中主要有均相催化和多相催化两种体系。均相催化虽然有较好的不饱和醇选择性,但是催化剂和溶液之间的分离很困难,三废比较多,工艺条件较为苛刻,不符合现代绿色催化的理念。多相催化多