发展型方程的连续时空有限元方法及其数值模拟

来源 :内蒙古大学 | 被引量 : 6次 | 上传用户:ldw521
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
连续时空有限元方法是一种高精度的数值方法.其对空间变量和时间变量统一进行处理,即不仅用有限元去离散空间变量而且用有限元去离散时间变量,因此相比于经典的有限元方法其更容易获得关于时间的高精度且其理论分析不会随着时间变量离散方式的改变而改变,即其理论分析对任意次的近似多项式都是一致成立的.此外,连续时空有限元方法特别适合于求解波动问题,因为其相应的离散格式具有重要的能量保守性.连续时空有限元方法可分为下面两种情形:1.每个时间层对应着相同的空间剖分;2.每个时间层允许对应不同的空间剖分.对于情形1由于每个时间层的空间剖分都相同,故在整个时间区间引入时空投影算子后易于获得时空有限元解在各种范数下的误差估计.对于情形2由于各个时空片的时空网格结构允许改变,所以其特别适合于无结构网格上的自适应计算.此外,对于情形2在理论分析中分别引入了由勒让德(Legendre)点和洛巴托(Lobatto)点确定的拉格朗日插值多项式及相应的高斯积分准则,并在理论分析中充分的利用了插值多项式的基本性质以及高斯积分的高精度特点,这使得理论分析变得自然易懂而且在更深层次上抓住了时间步有限元方法的本质.此外,不论是情形1还是情形2其连续时空格式往往是无条件稳定的并且其理论分析通常不受时空网格限制,即不需要时间步长和空间网格参数满足一定的条件.本文主要从理论分析和数值模拟两个方面在情形1和情形2下研究了与时间相关的偏微分方程的连续时空有限元方法.第一章是绪论,其主要阐述了连续时空有限元法的研究现状以及本文的研究内容与文章结构.此外,还给出一些本文理论分析所需的预备知识.第二章和第三章分别在情形1下研究了 Sobolev方程(不含对流项)和粘弹性波动方程的连续时空元方法.我们首先构造了原问题相应的连续时空元格式并证明了连续时空解的存在、唯一及稳定性,然后通过引入时空投影算子用相对简洁的理论分析在没有时空网格条件限制的情况下给出了节点处的L2和H1范数估计以及全局L2(L2)和L2(H1)范数估计.最后,我们给出一个无结构网格上的二维数值算例确认了格式的有效性和可行性.此外,数值算例还说明与传统的时间方向用Euler或者Crank-Nicolson(CN)格式离散的有限元方法相比,连续时空元方法更容易获得时空的高精度.第四章在情形1下我们研究了波动方程的连续时空元方法.在本章我们提出了一种新的用连续时空元求解波动方程的方法.我们通过引入勒让德多项式及其相应的高斯积分准则得到与原连续时空格式等价的格式,然后基于此格式分析了近似解的存在唯一性.此外,通过引入时空投影算子给出了近似解在节点处的L2及H1范数估计.最后,给出两个数值算例验证格式的有效性和可行性.与已有的方法相比,这里的分析方法更加简洁易懂并且容易被推广到其它的波动问题.同时我们需要指出的是在粘弹性波动方程和波动方程的求解过程中,我们首先通过引入辅助函数v=ut得到与原问题等价的耦合系统,然后基于此耦合系统构造了连续时空有限元格式,通过求解此格式可以同时获得u和v的高精度.第五章和第六章分别在情形2下研究了 Sobolev方程(不含对流项)和粘弹性波动方程的变网格连续时空元方法.我们首先构造了原问题相应的变网格连续时空元格式,其可以看作是第二章和第三章连续时空元格式的一种延扩.然后通过引入由勒让德点确定的拉格朗日插值多项式及相应的高斯积分准则给出数值解的适定性分析;通过引入由洛巴托点确定的拉格朗日插值多形式及相应的高斯积分给出了近似解的L∞(L2)和L∞(H1)范数估计.此外,在第六章我们还证明若每个时间层的网格满足一些合理的假设,则可以消去收敛性结果中的跳跃项,从而可以获得关于时空的最优阶L∞(L2)范数估计.第七章我们在情形2下研究了变系数的对流占优的Sobolev方程.首先证明了数值解的存在唯一性,然后在没有时空网格限制的情况下给出了最优阶L∞(H1)范数估计.最后,我们分别给出了原问题在连续时空有限元格式和时间间断的时空有限元格式下的数值模拟,数值实验验证了分析的正确性并展示在实际计算中连续时空有限元方法比时间间断的时空有限元方法更加有效.
其他文献
第四纪冰期–间冰期旋回中,河流与冰川在中低纬度的高海拔山地中相继占据主导地位,但两种外动力的侵蚀方式和能力都存在差异,可能对地形造成不同的影响。因此,了解冰川和河流
江西缫丝厂始建于一九五八年,多年来,由于管理不善,质量低劣,产品滞销,32年连续亏损.到一九九0年,企业处于停顿状态,仅从1988年至1990年的三年中,亏损额达200余万元,工厂濒临
建立了HPLC法测定注射用重酒石酸长春瑞滨的含量。采用氰基键合硅胶柱,流动相为0.05mol/L磷酸盐缓冲液(含0.8%三乙胺,pH2.8)-乙腈(70:30),检测波长268nm。线性范围为0.1~0.75mg/ml,平均回收率为
以6-甲氧基-8-氨基喹啉为原料,经重氮化后与间苯二酚偶联,合成了有机试剂4-(6-甲氧基-8-喹啉偶氮)-间苯二酚(简称MQAR),用DMF重结晶精制。通过TLC,EA,UV-VIS,IR,^1HNMR等分析手段,对合成产物进行分
目的:通过生物信息学挖掘血管再狭窄的潜在靶点,并构建内皮细胞损伤模型和血管再狭窄大鼠模型,应用苏木提取物进行干预,探讨苏木提取物抑制血管炎症及血管再狭窄的作用机制。
介绿了微贮饲料、青贮饲料及氨化饲料的加工原理及加工方法,并对其特点进行了比较。
各地、市、县委,各行署,各市、县人民政府,省委各部门,省直各单位党委、党组:在贯彻党中央、国务院关于治理经济环境、整顿经济秩序、全面深化改革的指导方针中,我省各地区、
目的 :探讨对进行腹腔镜手术的妇科疾病患者实施预防腹胀护理的效果。方法 :选择2013年10月至2016年10月期间在如东县人民医院接受腹腔镜手术的50例妇科疾病患者,对其临床资料
开发专一性β-内酰胺酶抑制剂,对解决细菌耐药性问题具有重要意义.这项研究需要用到纯度较高的β-内酰胺酶.我们将TEM-1型β-内酰胺酶基因克隆到一个高表达载体中,研究了该酶