求解带有界面和不规则区域的椭圆问题的浸入界面有限元方法

来源 :南京师范大学 | 被引量 : 0次 | 上传用户:waxizhaojing
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
在许多实际应用中,所研究的对象经常由不同物质组成,不同的物质通过界面相互分开.如果对这些实际问题建立微分方程的数学模型,那么在微分方程中,不仅参数是间断的,而且在界面上也需要满足一些界面条件.在本文中,我们考虑使用非匹配网格的浸入界面有限元方法来求解这些问题.对于界面问题,Li et al. (Numer Math 96:61-98,2003)提出的非协调浸入界面有限元方法在文献中已有许多深入的研究.该方法实际上是通过对传统的P1协调元进行修正,使其在界面单元上满足界面条件.非协调浸入界面有限元方法不仅程序实现简单而且其离散解在L2和H1范数下能达到最优精度.然而由于修正的基函数的不连续性,非协调浸入界面有限元方法在L∞范数下不能达到最优的二阶精度.尽管Li et al. (Numer Math 96:61-98,2003)也构造了完全二阶精度的协调浸入界面有限元,但是其构造比较复杂,不易程序实现.在本文的第一部分,我们通过改进原来的非协调浸入界面有限元方法,提出了一种对称相容的浸入界面有限元方法.该方法依旧使用非协调浸入界面有限元空间,因此不仅保留了非协调浸入界面有限元的一些优点,而且对称相容,更重要的是具有二阶精度.这个方法的思想是,针对修正的基函数的不连续性,在双线性型中加入一些修正项来保持相容性和对称性.接下来,我们把这个对称相容浸入界面有限元方法推广到带非其次跳跃界面条件的界面问题中.在本文的第二部分,我们提出两个增广的浸入界面有限元方法来求解界面问题和不规则区域问题.该方法实际上是一种快速迭代法.增广技巧首先由Li (SIAM J. Numer. Anal.35:230-254,1998)提出,并且应用在有限差分方法中.我们简单地把增广技巧应用到有限元框架中.得到第一个增广浸入界面有限元方法.在增广方法中,通过引进一个或多个在界面或边界上的增广变量,使得我们能更加容易的离散原来的微分方程.增广变量应该选取到使得界而或边界条件得以满足.增广方法成功的关键经常依靠一种插值来把增广变量和原来的微分方程耦合在一起.这通过最小二乘插值(系数不定)来完成.奇异值分解被用来求解插值系数.接下来,借助有限元的性质,我们提出了第二个增广浸入界面有限元方法.在这个新的方法中,我们避免了使用最小二乘插值.因此,这个新的增广方法比原来使用最小二乘插值的增广方法效率更高,程序更简单.然后我们把这个方法推广到带有狄氏边界条件的不规则区域问题中.我们提供了许多数值实验来展示这个新方法的精度和效率,包括带有任意界面/不规则区域的例子和系数跳跃很大的例子.数值结果表明GMRES迭代次数和网格尺度无关,而且和系数跳跃基本无关.
其他文献
在过去的几十年里,对于偏微分方程数值解的逼近,人们已经提出了各种各样的数值求解方法,如两网格方法,保结构数值方法等.这篇论文主要讨论了这些方法在某些偏微分方程中的应用.文章首先讨论了两网格有限体积元方法在非线性Sobolev方程中的应用,然后介绍了如何利用保结构方法,如离散变分导数方法及哈密尔顿边界值方法等构造保积分的数值方法.首先,我们针对非线性Sobolev方程提出了一类两网格有限体积元方法.
设X是一个实的无穷维的希尔伯特空间,(·,·)x是内积,||·||x恢是其上的范数.A:D(A)(?)X→X是一个无界自伴算子,它的谱集只含有离散谱σ(A)=σd(A),并且假设Φ满足:可微,并且对x∈Z,存在M>0使得|Φ’(x)y|≤ M‖y‖X,(?)y∈Z.(Φ0)意味着对任意x∈Z,都存在X中的元素▽Φ(x)使得对y∈Z都有Φ’(x)y=(▽Φ(x),y)x.我们考虑下面的算子方程:Ax
本文致力于两类界面问题的有限体积元方法的研究,全文共分为三个部分.第一章首先我们介绍了关于界面问题的一些浸入方法,阐述了发展浸入有限体积元方法的目的。然后我们介绍了一维带界面的双相延迟方程、高阶紧有限体积元方法和Pade型紧有限体积方法。第二章第一部分讨论了带界面的泊松方程的浸入有限体积元方法。通过源项移去技巧,将带非齐次跳跃条件的界面问题转化为带齐次跳跃条件的界面问题,与跳跃条件相关的项被转移到
本文定义紧致动力系统中的诱导拓扑压、诱导测度熵,研究它们的性质.具体的安排如下:在引言中,我们介绍动力系统中诱导拓扑压研究的背景.在第一章,我们介绍本文涉及到的遍历论和拓扑动力系统的预备知识.在第二章,我们定义紧致动力系统中的诱导拓扑压,研究诱导拓扑压与拓扑压的关系.在此基础上,得到诱导拓扑压的变分原理.作为诱导拓扑压的-个应用,指出BS维数是诱导拓扑压的特殊情形.我们还研究诱导拓扑压的平衡测度的
许多偏微分方程能被写成一个多辛哈密顿系统,例如:sine-Gordon方程、非线性薛定谔方程、KdV方程、Camassa-Holm方程、麦克斯韦方程、非线性波动方程等.多辛哈密顿系统有三个局部守恒律,即多辛守恒律,局部能量守恒律和局部动量守恒律.如何构造保其中一个或多个守恒律的数值算法是非常有意义的.多辛守恒律是多辛哈密顿系统的一个重要的几何性质.在过去的一、二十年里,人们发展了大量的保离散多辛守
本文结构如下.第一章主要介绍本文所需要的关于正规族和值分布理论的相关知识.第二章中,我们证明了两个涉及例外函数的全纯函数正规定则,改进了庞学诚、方明亮和Zalcman的一个正规定则,同时也是对刘晓俊、叶亚盛的一个相关结果的补充.第三章中,我们得到了一个涉及分担函数的亚纯函数正规定则,推广了常建明的相关结果.第四章中,我们通过减弱已有正规定则的某些条件,研究函数族在不正规点附近的性质,并在此基础上得
界面问题在材料科学、固体力学以及流体动力学中经常出现。例如,带有不同传导系数的热传导问题、描述不同材料行为的弹性问题以及带有不同粘性系数的两相流问题等。目前,对界面问题数值方法的研究已成为科学计算和工程领域内的研究热点之一。本文主要目的是在非拟合网格(即网格剖分和界面位置无关)下,构造一些新的有限元方法来求解界面问题并对其进行数值分析。首先,基于Nitsche方法和LDG方法的思想,我们提出了一类
在“数值算法应尽可能多地保持原问题的本质特征”的指导原则下,冯康先生首先提出了保结构算法的思想.由于其优良的稳定性和精确的长时间数值表现,目前保结构算法在求解哈密顿常微分方程上已经取得了显著的效果.但是,非线性波的传播,电磁场的演化等实际问题经常涉及无穷维哈密顿系统.对该系统的保结构算法研究还不够完善仍处于起步阶段,有许多基础理论与实际应用问题函待解决.因此,本文致力于进一步研究无穷维哈密顿系统的
为了进一步提升中国“无废城市”建设水平,分析了美国“无废”建设的发展历程、建设纲领、建设措施以及旧金山市“无废城市”建设案例,从管理体系和管理策略等方面总结了美国“无废”建设经验,并将其与中国发展现状相结合,提出中国“无废城市”建设的政策建议,包括废物减量管控、社会参与、制度与市场建设、建设措施评估这4个重点发展方向。
以孙中山为首的革命党人成立民国南京临时政府(以下简称:南京临时政府),推行民主共和,创立新的新闻管理体系,其对宣传资产阶级民主观念,促进民初新闻事业的发展,具有重要的进步意义。南京临时政府新闻管理体制的建立过程中,政府、新闻出版业及民众三方就新闻自由与新闻法制的认识和激烈博弈,对此后南京国民政府新闻出版法律体系的形成产生了重要影响。这一时期新闻管理体制在我国近代法制史上占有重要的历史地位,尽管它存