具有间断系数扩散方程的局部间断Galerkin方法

来源 :中国工程物理研究院 | 被引量 : 4次 | 上传用户:liaonianyou
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
在许多实际问题中,如材料科学和流体力学,会遇到具有间断系数的椭圆和抛物问题.如何准确有效地数值求解这些问题是一个很大的挑战.由于扩散系数在界面处间断,这类问题的解在整个定解区域上的正则性较低.真解较低的整体正则性和界面的不规则性,使得用经典的有限元方法求解很难获得高精度的数值解.  局部间断有限元(Local Discontinuous Galerkin,简记为LDG)方法是Runge-Kutta DG方法的推广,它具有局部守恒、易于处理复杂几何区域、允许解出现间断、形式上高阶精度、容易实现并行化和h-p自适应等特点,适用于求解通量比较重要的物理问题.本文研究了LDG方法在具有间断系数椭圆和抛物方程中的应用,给出了这两类方程的一些理论分析结果和数值试验.  第一章介绍了求解具有间断系数椭圆和抛物问题已有的数值方法及其进展情况,阐述了间断有限元方法的发展历史,给出了利用局部间断有限元方法研究界面问题取得的成果.  在第二章,我们详细地讨论了一维含间断系数抛物问题的局部间断有限元方法,给出了方法的L2稳定性分析,并证明出当采用k次多项式空间逼近时,连续时间LDG格式的收敛阶为O(hk+1/2).文中分别用显式和隐式时间离散求解LDG格式,数值试验验证了方法的有效性和理论分析结果.  第三章研究了二维线性椭圆界面问题的最小耗散局部间断有限元(MD-LDG)方法,间断面为任意光滑封闭曲线.尽管解的整体正则性较低,我们证明出该方法的解和通量在L2模下的收敛阶分别为O(h2|logh|)和O(h|logh|1/2).在数值试验中,采用逐次替换迭代方法求解LDG格式,所得的数值结果与理论分析相符合.  在第四章,考虑二维线性抛物界面问题的最小耗散局部间断有限元方法,间断面仍然为任意封闭光滑曲线.我们证明出该方法是L2稳定的,并且在能量模下数值解的收敛阶为O(h|logh|1/2),最后给出了数值试验.  第五章讨论的是二维含间断系数非线性抛物问题的最小耗散局部间断有限元方法,间断面为直线.我们给出了LDG格式的L2稳定性,并证明出能量模下数值解的收敛阶为O(h),数值算例验证了方法的有效性.
其他文献
学位
学生的心理健康受到了社会各界越来越多的关注,而传统的心理教育模式已经不能满足当前的需求,随着信息技术的发展,心理教育空间也需要不断地拓展,而微博无疑是其中一种很好的
该文首先研究了不确定大系统经由状态反馈的鲁棒H控制问题,此类系统带有不确定非线性关联项,且不确定项是时变有界的.基于Riccati方程工作人员可以设计出使系统分散镇定的状
半群的自动理论开始于20世纪末,近年来许多学者研究半群的自动理论.本文收集整理了该理论的近期成果和基本的研究方法,以利于读者进一步的学习和研究该理论.本文给出了Gr(o)bne
矩阵对策是对策论的基础,一直备受关注,是研究的热点。在矩阵对策中广泛研究的就是当策略集模糊、支付值精确或策略集清晰、支付值模糊,或策略集模糊、支付值也模糊的情况下对策
该中针对现代试井分析中存在的问题,提出了利用生产测井同时测量的关井压力和续流,结合最优化手段,用自动历史拟合方法处理关上期资料试井分析的方法,有效地消除了井筒储存效
该文针对压缩图像的JPEG编码标准的一些不足,例如会产生边缘效应等,提出了基于自适应投影算法的压缩编码系统,其主要由三个部分组成:1.自适应投影算法;2.格型矢量量化方法;3.
数据库知识发现是当今计算机领域研究的热点.该文首先概述了知识发现的种类、目标和方法.介绍了粗造理论应用于知识发现过程的研究现状.评述了运用面向属性的归纳方法进行知