弱有限元方法求解具有交界面的偏微分方程

来源 :吉林大学 | 被引量 : 0次 | 上传用户:liongliong554
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文研究两种带有交界面的模型的数值求解,分别是Stokes-Darcy耦合模型和线弹性交界面模型.在Stokes-Darcy模型中,交界面两侧是不同的控制方程,方程在交界面上交换信息,这类问题的研究更侧重于不同方程之间的耦合,因此在对这样的耦合问题设计算法时,既需要考虑数值方法对不同区域上不同方程的求解效果,又要考虑数值方法对整体系统的求解效率.而对于线弹性交界面模型,交界面的两侧是同一种方程.在穿过交界面的时候,方程的系数,位移,应力和应变等信息均可能在交界面上产生间断,这会对数值算法的设计带来一些挑战,比如这样的间断性会使得问题具有比较低的全局正则性.因此,基于一些正则性假设的数值方法就会失效.此外,交界面几何形状的复杂性也需要网格具有较强的灵活性,对于线弹性问题的求解,另一个需要关注的点是闭锁现象,即当材料趋于不可压缩时(方程中的其中一个Lam’e常数趋于无穷),数值解的误差会随着常数的增大而改变.本文将根据这两种模型的特点,以弱有限元方法为主对这两种模型设计数值算法.弱有限元方法是近年发展起来的一种新型数值方法,其主要特点是引进了弱函数作为有限元空间中的函数,并对弱函数定义相应的弱微分算子.通常,弱函数在形式上记为{v0,vb}在每个单元上,v0表示函数在单元内部的值,vb表示函数在单元边界上的值,vb与v0在单元边界上的迹之间可以是无关的.因此,我们在作离散的时候,可以根据具体需要选择不同的v0,vb多项式组合.此外,为了保证数值格式的适定性,还需要在弱有限格式中加入稳定子.由于弱有限元的定义是逐单元进行的,该方法允许使用一般的二维多边形网格和三维多面体网格.可以看到,弱有限元方法在离散元的选取和网格剖分方面都具有较强的灵活性,这对于处理具有交界面的问题是非常便利的.Stokes-Darcy模型经常被用来描述自由流和多孔介质流耦合的情形,具体而言,在自由流区域使用Stokes方程描述流体的流动,在多孔介质区域使用Darcy定律描述多孔介质流,在交界面上还需要三个交界面条件.第一个交界面条件来源于质量守恒,第二个交界面条件描述的是法向应力平衡,第三个交界面条件是Beavers-Joseph-Saffman(BJS)条件.Darcy方程有两种不同的形式,即原始格式和混合格式.我们首先考虑标准的Sto-kes方程耦合原始格式的Darcy方程.对于Stokes方程,弱有限元方法在求解时表现出来以下优势.首先,Stokes方程的弱有限元格式构造简单直接且不需要调整参数来保证格式的稳定性.第二点,可以使用低阶的弱有限元对来逼近速度和压力.第三点,弱有限元方法允许使用一般的多边形网格.对于原始格式的Darcy方程,实际上是一个二阶椭圆的形式.标准有限元对于二阶椭圆方程的求解是高效且成熟的.因此,对于这样的Darcy方程,使用有限元方法进行离散是首选.然而对于Stokes-Darcy耦合问题的求解,并不是简单的将两种离散方式结合在一起就可以得到上面所说的优势,因为离散格式在交界面处需要具有一致性.我们可以充分利用弱有限元方法的灵活性,仔细选择Stokes区域的局部弱有限元匹配Darcy区域整体的连续元.此外,由于弱有限元方法允许使用带有悬点的网格,因此交界面两侧的网格也不必是对齐的.在对数值格式的分析中,为了处理Stokes方程中的应力张量,我们证明了针对弱有限元的Korn不等式,并在此基础上,证明数值格式的适定性.进一步,我们还推导了误差估计,并通过一些数值算例验证了理论分析得到的结果.紧接着,我们又考虑Stokes方程与混合格式的Darcy方程的耦合.混合格式的Da-rcy方程涉及速度和压力两个变量,这使得耦合系统在两个子区域上的未知量一致,因此,混合情形比原始情形在描述交界面条件时更加自然.对Stokes方程稳定的元和对混合格式的Darcy方程稳定的元经常是不同的,尽管弱有限元方法可以同时对Stokes方程和Darcy方程稳定,但是出于自由度的考虑,我们仍然在不同的区域使用不同的元来离散.在我们的工作中,我们仍然使用弱有限元方法来离散Stoke方程,而对于混合形式的Darcy方程,我们使用混合有限元方法进行离散,一些经典的混合元可供选择,例如RT元,BDM元和BDFM元等.对于这个离散系统,一个主要的困难是整个区域上的inf-sup条件,我们在Stokes区域和Darcy区域分别构造了投影算子来解决这一难点.进一步,我们推导了两个区域上速度和压力的误差估计,并给出了一些数值算例来说明我们的数值方法的有效性.本文考虑的另一类问题是线弹性交界面问题.线弹性方程广泛应用于固体力学中,用来描述固体在外力的作用下,位移,应变和应力之间的关系.然而在实际生活生产中,物体往往由多种不同材料组成,这些材料的材料性质和物理性在交界面上会发生改变,这就导致了所谓的线弹性交界面模型.在本文中,我们使用弱有限元方法离散线弹性交界面模型.由于方程的系数,位移,应变和应力在穿过交界面时会产生间断,我们首先在交界面上使用双值的弱函数来离散模型,但这样的格式在理论分析和算法实现上比较繁琐,于是在此基础上,我们将交界面条件带入到弱有限元格式中,得到了在交界面上具有单值的弱有限元格式.进一步,我们推导了位移在能量范数和L~2范数下的最优阶误差估计,估计结果说明了所构造的弱有限元格式是无闭锁的.最后,我们通过一些数值算例验证了数值格式的有效性和无闭锁性.
其他文献
医学科普关系全民健康,通过各类医学知识的科普宣传,可以提升公众健康管理意识,推广“每个人是自己健康第一责任人”的健康理念,助力健康中国建设。全媒体时代,医学科普图书在保证医学健康知识传播的专业性、真实性与科学性方面具有明显优势,但也存在内容同质化、创新性不足,难以实现社会效益与经济效益相统一等问题,制约医学科普图书的高质量发展。选题策划是图书出版的源头,具有前瞻性的优质选题是医学科普图书高质量发展
期刊
教师职业认同是教师履行应尽义务、做好本职工作的基础。它密切联系着教师的职业发展、学生的健康成长以及国家建设人力资源强国战略的实施。在国家全面深化教育体制改革,推进教育现代化、促进学生全面发展的进程中,中共中央办公厅、国务院办公厅于2021年7月24日印发了《关于进一步减轻义务教育阶段学生作业负担和校外培训负担的意见》(简称"双减"政策),“双减”政策的提出,减轻了学生负担,规范课后延时服务的同时,
学位
本文致力于研究具时滞种群动力学模型的相关问题,我们分别从有界区域和全空间的角度考虑具时滞反应扩散方程.当模型的空间扩散是退化扩散形式时,具有生态学意义且非局部出生率函数涉及全空间,因此本文主要研究具退化扩散和非局部模型的解的存在性、唯一性和正则性.特别地,我们考虑行波解的稳定性以及解收敛到非平凡稳态解的收敛率.在本文的第一章中,我们介绍种群动力学模型在生态学中的重要意义、模型演化过程以及最近的研究
学位
职业认同是特教教师专业成长的关键因素,也是特教教师成长的心理支撑。通过对吉林省211名特教教师的职业认同现状进行调查的结果表明,其职业认同处于中上水平。不同性别、不同婚姻状态对教师的职业认同有显著影响。不同年龄、不同教龄的教师,其职业认同程度有显著差别。特教专业高级职称教师较多,职称会对教师的职业认同产生较大的影响。特教学校对师资的需求较大,师资队伍的职业认同会受到培训因素的影响。
期刊
本文主要研究了 3-李代数上相对罗巴算子的形变理论和上同调理论,以及3-李代数上相对罗巴算子在3-李双代数和3-李代数的Manin triples等领域的应用.首先,我们通过3-李代数和它的表示及高阶导出括号,构造了一类李3-代数(一类特殊的李无穷代数),证明了此类李3-代数的Maurer-Cartan元素就是3-李代数上的相对罗巴算子.通过相对罗巴算子进行扭曲运算后,还可以得到一类新的李无穷代数
学位
本文研究了几类整值门限自回归过程的建模方法问题.首先,为了刻画整值时间序列同时受多个外部因素影响的数据特征,我们提出一类系数由协变量驱动的一阶随机系数整值门限自回归模型.讨论了模型的基本性质以及模型参数的估计问题,给出了模型的非线性检验和解释变量是否存在的检验方法,并将其应用于一组犯罪数据的实例分析中.其次,为了更好地刻画整值时间序列具有偏大离差(overdisperson)、偏小离差(under
学位
自组织系统来源于人们对动物集群行动、网络信息交互以及市场物价调节等问题的研究,是对具有多个个体相互作用并达到一致的现象的统一描述.自20世纪60年代以来,自组织系统在经济学、生物学和系统工程等领域的广泛应用,极大地促进了相关理论的发展.人们先后建立了著名的Viseck模型、Olfati-Saber模型、Cucker-Smale(C-S)模型等,并取得了一系列重要的研究成果.本文对两类自组织系统进行
学位
<正>一、前言化学镀铜是一种不用外加电流在绝缘材料表面上沉积铜的方法,是印制板制造不可缺少的一步工艺。化学镀铜反应整个过程都存在着电子得失:Cu2++2e-Cu(1)HCHO+3OH--HCOO-+2H2+2e(2)式(1)表明了铜离子得到电子而被还原成金属铜,式(2)表明甲醛在强碱条件下放出电子而被氧化,所以化学镀铜是一种自身催化性氧化还原反应。其主要反应为:
会议
电扩散,即电荷的扩散和迁移,在物理学、化学和生物学等研究领域中都有非常重要的作用.Poisson-Nernst-Planck(PNP)模型是描述离子在电化学梯度作用下运动的常用模型,它由Poisson方程和Nernst-Planck方程组成,其中Poisson方程描述移动离子和固定电荷产生的静电场,Nernst-Planck方程描述离子在静电场中的运动.PNP模型已被广泛应用于半导体中载流子的输运
学位
概率统计、风险度量和数理经济学等领域中出现的模型不确定问题,经典概率理论中的方法已不再适用,各种非线性期望的出现在一定程度上解决了这个问题,目前已经得到了广泛的研究。在实际情况中我们遇到的随机样本往往不是独立的,于是学者逐渐发现突破对独立性的限制,对相依序列展开研究更具有实际意义。本文从非线性期望中最具代表性的次线性期望出发,研究了相对于独立序列条件更弱的m-相依序列在次线性期望中的极限性质和相关
学位