发展型方程的混合有限体积元方法及数值模拟

来源 :内蒙古大学 | 被引量 : 3次 | 上传用户:geshufa
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
混合有限体积元方法是将混合有限元方法和有限体积元方法相结合的一种数值方法,该方法又称为混合控制体积方法,最早是由Russell于1995年在求解一类二阶线性椭圆型问题时提出,随后Cai和Jones等人通过数值算例验证该方法的有效性,Chou和Kwak等人在该方法的理论分析方面做了大量的工作.该方法具有如下优点:采用混合元思想引入辅助变量(如梯度函数,流函数等)将高阶问题转化为低阶问题,降低对有限元空间的光滑度的要求;易于处理复杂区域和边界条件,具有有限体积元方法的格式简单性;方法的定义采用混合变分形式,可以从混合变分形式出发进行理论分析;计算量比混合有限元方法小,收敛阶和混合有限元方法相同;能够保持某些物理量(如质量,动量)的局部守恒性质.Rui和Lu于2005年将扩展混合元方法和有限体积元方法相结合,对一类二阶椭圆问题构造了矩形网格剖分下的扩展混合有限体积元方法,该方法继承了扩展元方法和有限体积元方法的优势,可以同时数值计算三个未知变量.目前关于此类方法的研究主要是基于矩形网格剖分,而在三角网格剖分下的研究还非常少,本文主要是在三角网格剖分下将此类方法求解含对流项Sobolev方程,并给出误差分析和数值模拟.近年来随着解决复杂实际数学物理问题的需要,对数值方法在计算效率上也有了更高的要求.本文结合分裂思想对扩展混合有限体积元方法进行简化,提出了一类新型的分裂扩展混合有限体积元方法.该方法在数值计算时可以先求解两个方程的耦合系统得到两个变量的数值解,再利用这两个数值解求解第三个方程得到第三个变量的数值解,这种办法在很大程度上降低了线性方程组的规模从而大幅度的减少计算时间.本文应用混合有限体积元、扩展混合有限体积元以及分裂扩展混合有限体积元方法从理论分析和数值计算两个方面对几类发展型方程进行研究,由于每一类方程都有不同的特点,从而所构造的数值格式也不相同,而且需要根据每一类方程的特点进行相应的理论分析和数值实验.在第一章中简单介绍一下混合有限元方法和混合有限体积元方法的特点和发展现状.在第二章到第四章中应用混合有限体积元方法研究三类发展型方程.其中在第二章研究了一维正则长波方程的混合有限体积元方法.通过引入一维网格剖分下的迁移算子,构造了半离散、非线性和线性向后Euler全离散混合有限体积元格式,利用椭圆投影和L2正交投影算子给出三种离散格式的最优阶误差估计,最后通过数值算例验证格式的有效性和收敛精度.第三章和第四章应用混合有限体积元方法在三角网格剖分下数值求解一类二维伪双曲型方程和非线性阻尼Sine-Gordon方程.选用最低阶Raviart-Thomas有限元空间和分片常函数空间作为解函数空间,并引入迁移算子γh将最低阶Raviart-Thomas有限元空间映射成试探函数空间,对两类方程分别构造了半离散和关于时间隐式的全离散混合有限体积元格式.通过引入广义混合有限体积投影得到半离散和全离散格式的最优误差估计,最后对两类方程分别给出一些数值结果验证了该方法的可行性.第五章将扩展混合元和混合有限体积元方法相结合构造了一类含对流项Sobolev方程的初边值问题的扩展混合有限体积元方法.该方法引入辅助变量λ=-▽u和σ=-(a▽u+6Vut),将原问题降为一阶微分方程系统,选用最低阶Raviart-Thomas有限元空间作为变量λ和σ的解函数空间,并使用分片常函数空间作为u的解函数空间,利用迁移算子γh在三角网格剖分下构造了半离散和向后Euler全离散的扩展混合有限体积元格式.应用微分方程理论证明了半离散格式解的存在唯一性,利用迁移算子的性质和扩展混合有限体积投影得到半离散和全离散格式的最优阶误差估计.最后给出数值算例验证了方法的可行性和理论分析的正确性.第六章将扩展混合有限体积元方法和分裂思想相结合,引入和第五章一样的辅助变量λ和σ,构造了含对流项Sobolev方程的一种新型分裂扩展混合有限体积元格式.这一格式与扩展混合有限体积元格式的区别在于:扩展混合有限体积元格式需要同时求解三个方程的的耦合系统,从而在数值计算过程中生成的线性方程组的系数矩阵规模比较大,而此格式在数值计算时可以先求解两个方程的耦合系统得到变量λ和σ的数值解,再利用这两个数值解求解第三个方程得到变量u的数值解,这种办法在很大程度上降低了线性方程组的规模从而大幅度的减少计算时间.最后给出一些数值结果来验证该方法的有效性和理论结果的正确性.第七章研究了一类抛物型积分微分方程的初边值问题的分裂扩展混合有限体积元方法.引入辅助变量λ(x,t)=-(?)u(x,t)和σ(x,t)=-(a(x)(?)u(x,t)+∫0t(x,t,τ)(?)u(x,τ)dτ),利用迁移算子γh构造了半离散和向后Euler全离散分裂扩展混合有限体积元格式,其中在全离散格式中利用左矩形数值积分公式离散积分项,利用向后Euler格式离散时间导数项.引入Volterra型扩展混合有限体积投影并利用迁移算子的性质得到了两种格式的最优阶误差估计,最后通过数值算例验证了理论分析结果.
其他文献
目的:评价中药晕平方治疗肝阳上亢型慢性主观性头晕(CSD)的临床疗效。方法:采用随机对照方法,将符合纳入标准的CSD患者108例随机分为治疗组、对照组。两组均提供心理教育及盐
随着(移动)互联网技术的变革以及社交网络的广泛普及,人们每天都会接受来自互联网上各种各样的信息。伴随着互联网技术从1.0时代到2.0时代的迈进,每天又会有大量的网络使用者
该文以上证50样本期内38只股票周收益率为实证研究对象,利用GARCH模型预测周收益率和条件方差,并以此分别量化传统BL模型中的观点收益向量和观点误差矩阵,得到的BL预期收益率
本文从当前高校药品管理中的典型问题出发,围绕着高校实验室教学与科研工作的需求,探讨创新性的药品管理方式,运用标准化的方法提高药品管理工作效率,解决药品管理中的实际问
烟草企业作为我国国民经济的特殊而重要的一员,从组建发展到今天,企业管理不断完善,取得了一定成就。但是,在当前市场化经济和现代化发展的背景下,烟草现代企业发展也受到诸
热化学驱作为超稠油/稠油油藏热力采油后的接替手段仍处于论证阶段。相对渗透率曲线反映了驱油剂和原油在多孔介质中的流动规律,是油藏开发中的重要资料。在稠油热化学驱油过
背景:股骨粗隆间骨折是临床中最常见的髋部骨折之一,好发于老年人,主要与老年人易摔倒及骨质疏松有关。现在骨科医师大多主张进行手术治疗,随着医疗技术发展,骨科医学者们从不同影响因素入手进行研究,力求最大限度的避免并发症和畸形的发生。但国内外在对手术时机的选择上却一直未有统一意见,本研究力求从临床病例入手,探寻手术时机与PFNA术后疗效之间的关系。目的:通过对应用PFNA治疗的老年股骨粗隆间骨折的患者选
信息网络传播权是网络时代背景下发挥著作权作用、实现作品经济和社会价值的保障,也是我国《著作权法》2001年修改新增的权利条款。随着信息传播技术的发展,原有的信息网络传
利用抗菌模型和肿瘤细胞毒模型,对浙江舟山群岛12种海洋动物中的生物活性物质进行了筛选.结果表明,7种动物的乙醇提取物有明显的抗金黄色葡萄球菌作用,其中,泥螺乙醇提取物的
关于从宽的讨论由来已久,如自首从宽、立功从宽等,《刑法修正案(八)》的施行,意味着坦白从宽正式被写入刑法,成为我国重要的刑事法律制度之一。自首、坦白从宽的有关规定也被