求解偏微分方程的无稳定子格式弱有限元方法

来源 :吉林大学 | 被引量 : 1次 | 上传用户:xltmzzd
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文主要内容是对无稳定子格式的弱有限元方法的研究,弱有限元方法是由王军平和叶秀在2013年最早提出,以求解二阶椭圆问题的一种数值方法,这种方法的主要思想是在数值格式中借助重新定义的广义弱微分算子,对剖分单元的内部和边界分别引入自由度,即使用分片的,间断的多项式函数空间作为近似空间,空间基函数的定义是在各单元分片进行的,且函数的内部抑和边界ub定义相互独立,因此,弱有限元方法具有有限元空间更大,算法更灵活,数值格式更简洁的优点.同时,为了保证数值解的弱连续性,一般的弱有限元方法需要引入刻画相邻单元间差异的项,通常被称为稳定子,稳定子和弱微分算子(以弱梯度和弱散度为例)的一般格式定义为(▽wv,τ)T=-(v0,▽·τ)T+<vb,τ·n>(?)T,(?)τ ∈[H1(T)]d×d,(▽w·v,q)T=~(v0,▽q)T+<vb,qn)(?)T,(?)q∈H1(T),s(v,w)=(?)<vo-vb,ω0-ωb)(?)τ.不同的问题模型,不同的弱函数空间定义都会影响稳定子s(v,ω)的定义方式,由它的定义可以看出,稳定子的加入会导致数值求解的计算复杂度的增加,因此,我们考虑如何在保证弱有限元方法其他优点的同时,可以去掉稳定子,从而降低算法的计算复杂度.在本文中,我们主要介绍两种无稳定子格式的弱有限元方法.第一种是协调间断有限元方法,目前这种格式已被用于求解二阶椭圆方程和Stokes方程.这种格式的主要思想是对弱微分算子引入Raviart-Thomas(RT)元,同时利用内部函数的平均值代替边界函数,这种格式的近似函数空间的定义与传统弱有限元方法相同,并保证了相同的精度,在格式更简单的情况下,减少了边界自由度,去掉了稳定子,有效降低了计算的复杂度.本文第三章中,将使用矩形剖分上的协调间断有限元格式求解二阶椭圆方程.为了使格式简洁,我们考虑Rd上的有界区域Ω上具有Dirichlet边界的Poisson问题:-△u=f,在 Ω,u=g,在(?)Ω.在使用传统的协调有限元方法和弱有限元方法对上述Poisson问题进行求解时,其数值格式分别为:求uh∈H1(Ω),满足uh|(?)Ω=g,以及uh={u0,ub}∈Vh,满足ub|(?)Ω=Qhg,使得(▽uh,▽v)=(f,v),(?)v∈H01(Ω),(▽wuh,▽wv)+s(uh,v)=(f,v),(?)v={v0,vb} ∈ Vh0,其中Qh={Q0,Qb}是投影算子,▽ω是弱梯度算子,它的定义来源是分部积分方法,Vh是相应的弱有限元函数空间,Vh0是Vh的子空间,其在边界(?)Ω的取值为0,具体定义我们将在第二章和第三章中给出.接下来,我们给出相应的协调间断有限元方法的数值算法.算法1求uh={u0,{u0}} ∈ Vh,且在(?)Ω上满足uh=Ihq,使得(▽ωuh,▽ωv)=(f,v),Vv∈Vh0,其中{u0}表示u0的均值,Ih是拉格朗日插值算子,Vh和Vh0是相应的弱函数空间及其子空间.由以上三种数值格式的对比可以看出,协调间断有限元方法的数值格式既具有协调有限元方法的简洁性,又具有弱有限元方法的灵活性.接下来,我们给出了稳定性分析,更近一步地,我们建立了在离散H1和L2范数下最优收敛阶的误差估计,最后我们进行了数值实验,数值结果验证了理论分析所得结论.本文第四章和第五章主要介绍第二种无稳定子格式的弱有限元方法,这种方法的主要思想是通过提高弱微分算子的次数到一定程度,在满足数值解的弱连续性,保证刚度矩阵的大小和全局稀疏性不变的同时,使用更高次数的弱微分算子也可以显著降低数值计算时的复杂度.目前这种格式已被用于求解二阶椭圆问题,Stokes问题和线弹性问题.仍以二阶椭圆问题为例,求解这一问题的无稳定子弱有限元方法的数值格式如下.算法2求uh={u0,ub}∈Vh,且在(?)Ω上满足uh=Qhg,使得(▽ωuh,▽ωv)=(f,v),(?)v ∈Vh0,其中Qh={Q0,Qb}是L2投影算子.可以看出,这种数值格式也具有限元方法的简洁和弱有限元方法的灵活.对于任意v ∈Vh+H1(Ω),其弱梯度▽ω(v|T)∈[Pj(T)]d(j>k),j的取值将在后续章节中给出.在本文中,我们会继续对这种方法进行研究,并将其用于求解Stokes问题和线弹性问题,给出相应的数值格式.我们对数值格式的适定性进行了证明,给出了收敛性分析,并进一步通过数值实验验证了理论分析.特别地,由于在这种格式下,弱微分算子次数的选取与剖分单元的边数有关,因此在数值实验中,我们不仅给出了常规三角形和矩形剖分的情况,同时也给出了六边形剖分的情形,进一步论证了理论分析.
其他文献
本文致力于研究具时滞种群动力学模型的相关问题,我们分别从有界区域和全空间的角度考虑具时滞反应扩散方程.当模型的空间扩散是退化扩散形式时,具有生态学意义且非局部出生率函数涉及全空间,因此本文主要研究具退化扩散和非局部模型的解的存在性、唯一性和正则性.特别地,我们考虑行波解的稳定性以及解收敛到非平凡稳态解的收敛率.在本文的第一章中,我们介绍种群动力学模型在生态学中的重要意义、模型演化过程以及最近的研究
学位
职业认同是特教教师专业成长的关键因素,也是特教教师成长的心理支撑。通过对吉林省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-相依序列在次线性期望中的极限性质和相关
学位
本文研究两种带有交界面的模型的数值求解,分别是Stokes-Darcy耦合模型和线弹性交界面模型.在Stokes-Darcy模型中,交界面两侧是不同的控制方程,方程在交界面上交换信息,这类问题的研究更侧重于不同方程之间的耦合,因此在对这样的耦合问题设计算法时,既需要考虑数值方法对不同区域上不同方程的求解效果,又要考虑数值方法对整体系统的求解效率.而对于线弹性交界面模型,交界面的两侧是同一种方程.在
学位
随着人民生活水平的提高,伴侣动物如犬、猫等受到越来越多家庭的喜爱,对于宠物粮食的需求也由传统粮食逐渐向能够改善其健康的功能粮方面转变。益生菌是一类能够在动物体内定植的活性微生物,具有调控肠道菌群平衡,维持动物机体健康,促进营养吸收等作用。文章对益生菌的定义、菌株选择标准及其在宠物粮食中的应用进行了综述,分析了不同加工方式可能对益生菌的影响,旨在为科学合理地加工含益生菌的功能粮食,保证其生物活性、提
期刊