三种不可压缩流的稳定化数值方法

来源 :华东师范大学 | 被引量 : 0次 | 上传用户:hutianyi199052
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本论文研究了不可压缩流动问题的三类不同方程的有限元稳定化方法。主要内容分为三个部分,第一部分针对不可压缩流动多区域耦合模型,考虑了与大气-海洋耦合系统具有相同数学结构的Stokes-Stokes耦合的简化形式,即具有动力学交界面条件的两个Stokes耦合问题,其交界面上施加了 Nitsche型边界条件。这类边界条件会造成理论分析中高阶积分项很难被有效控制,这是具有Nitsche型边界条件的Stokes-Stokes耦合模型的一个挑战性工作。为此,我们引入了全局Nitsche稳定项和压力松弛项,进而构造了新的离散格式,以保证模型有限元解的连续性和唯一性。由于考虑了界面上的质量守恒,该稳定项和松弛项在连续情况下不会对原始模型产生影响。我们引入弱全局强制性这一种简洁的思路证明了该稳定化方法的稳定性,同时还得到了其最优误差估计。具有光滑真解的算例和模型算例验证了我们的稳定化方法在Stokes-Stokes耦合问题模拟中的有效性,特别的,我们设计并构造了管道流模型,得到了 Stokes-Stokes耦合模型及其算法与物理实际现象非常吻合的效果。论文的第二部分考虑了非牛顿流体的流动,特别是稳态粘弹性流体流动问题,主要研究了具有Oldroyd-B型本构关系的线性化的Oseen型粘弹性流体流动方程。粘弹性流体流动问题的混合格式导致了三类数值不稳定性:一是应力场计算可能面临的模型方程不稳定性,二是对流占优造成的数值格式不稳定,三是混合有限元空间的选择需要满足LBB条件。本文主要采用了间断Galerkin(DG)方法和Streamline Upwind Petrov-Galerkin(SUPG)方法,重点研究了低阶有限元对选择下的稳定化方法。我们拟采用最低等阶有限元方法求解粘弹性流体流动问题。然而,粘弹性流体流动模型中存在三个未知量:速度、压力和应力张量,因此我们对DG格式选择最低等阶有限元三元组P1-P1-Pdc1,对SUPG格式选择最低等阶有限元三元组P1-P1-P1,进而对粘弹性流体流动模型进行数值逼近。众所周知,这些有限元对不满足LBB条件,因此,对应的数值格式不稳定。但是低阶有限元对具有构造简单,易于并行等良好的性质。为了克服低阶有限元对不满足LBB条件这一困难,我们在有限元格式中加入了对应的压力稳定项,构造了新的稳定化方法,该方法具有无惩罚参数、计算灵活、不需要高阶导数等优点。通过证明全局弱强制性,我们得到了该方法的稳定性,同时分析了其最优误差估计,相关的数值试验很好地验证了算法的稳定性和有效性。论文第三部分主要针对磁流体(MHD)方程组,研究了拟最小二乘稳定化混合有限元方法。该方法不需要混合有限元空间满足LBB条件。对于非线性问题,拟最小二乘有限元法比最小二乘法有更多的优势。首先,因为该稳定化方法中只引入了L2范数,易于编程实现。其次,通过线性化,我们可以得到一个简单的带有对称正定阵的迭代步骤,而且这些简单迭代步骤对于初值的选取具有较大收敛域。此外,该方法具有一致收敛性,不依赖于算法设计及粘性系数。我们研究了拟最小二乘有限元方法关于MHD方程组数值解的存在性,并且得到了先验误差估计,数值实验很好地验证了我们的理论分析。
其他文献
各种环境胁迫均会使内质网(Endoplasmic reticulum,ER)的稳态失调,导致内质网胁迫(ER stress)的发生。胁迫下被诱导激活的未折叠蛋白质反应(Unfolded protein response,UPR)可通过减少内质网中错误折叠的蛋白质积聚来维持细胞存活,而内质网胁迫过于严重时,UPR则促进程序性细胞死亡(Programmed cell death,PCD)的发生。动物细
单光子态在量子信息处理、量子通讯以及量子计算等领域的应用前景非常广泛,因此其已经成为量子光学领域最重要的研究对象之一。单光子可以由不同的介质系统产生,比如NV色心、原子气体、量子点、单原子及单分子等。但是无论是由何种源产生的单光子,都需要测量单光子态的品质,比如线宽,谱线明亮度,单光子纯度,不可分辨率,态密度矩阵等。本文所研究的单光子态是在85Rb冷原子系综中通过四波混频过程产生的宣布式单光子源。
在当前“传统复归”的社会思潮和弘扬优秀传统文化的国家政策背景下,传统文化复兴及其儒学建设,已经成为具有当代中国文化特征的社会行动与文化景观。这主要表现在思想文化层面的分歧与论争,社会实践中的多元主体及其不同实践取向,以及儒学复兴中的“制度性焦虑”和“神圣性危机”。上述特征彼此交错且关联密切,形成了儒学复兴的社会学问题。儒学复兴的社会行动包括观念认知、仪式表达、社会行动、制度安排、神圣建构等一整套形
突触可塑性是学习记忆的重要机制,其主要包括长时程增强(long-term potentiation,LTP)和长时程抑制(long-term depression,LTD)。LTP在时间上分为早期LTP和晚期LTP。早期LTP主要涉及蛋白质的翻译后修饰,而晚期LTP则需要蛋白质的重新合成。已有大量研究表明磷酸化修饰参与了早期LTP的发生,而其他类型翻译后修饰在早期LTP形成中的作用还知之甚少。本课
基因转录的精确调控对于生物体的生长发育至关重要,也与疾病的发生发展密切相关。发生在DNA胞嘧啶上的甲基化修饰,即DNA甲基化,是一种重要的参与基因转录调控的表观遗传修饰。位于基因启动子区域的DNA甲基化修饰常与转录抑制相关,然而关于全基因组范围内启动子区域DNA甲基化如何影响基因转录及其精确的分子机制和生物学意义仍然没有被完全理解。DNMT3家族蛋白的ADD结构域能够识别并结合未甲基化修饰的H3(
贝叶斯非参数统计是一个相对年轻并且发展速度比较快的一个统计领域,它不仅产生了大量的理论成果,而且被广泛的应用到各种实质的领域和方向中。相比于传统的频率学派下的非参数统计,贝叶斯非参数是在贝叶斯框架下,基于它的后验或者预测分布来进行统计推断的。而相比于参数贝叶斯,贝叶斯非参数是建立在分布空间上的先验并且具有很大的支撑空间,而传统的参数贝叶斯只是基于参数空间建立的。尽管贝叶斯非参数在实际应用中面临着很
粗几何里面的一个重要问题之一就是粗Baum-Connes猜想,而研究粗Baum-Connes猜想的一种有效方法就是量化K-理论.量化K-理论比经典的K-理论要更为灵活更方便.我们可以利用量化K-理论的Mayer-Vietories序列来证明Baum-Connes猜想成立.或者通过量化K-理论里面传播(Propgation)来寻找Baum-Connes猜想的障碍.全文共分为四章,具体如下:第一章主要
功能加密(Functional Encryption,FE)是一种涵盖了身份基加密(Identity-Based Encryption,IBE)、属性基加密(Attribute-Based Encryption,ABE)和内积加密(Inner-Product Encryption,IPE)等诸多新型公钥加密原语的通用概念。凭借其对于加密数据共享的细粒度访问控制能力,以及迎合各种新型网络服务(如外包
随着生产材料和生产工艺的进步,产品的寿命越来越长.因此人们转而通过退化试验中的退化数据对产品的寿命进行外推.目前为止大量的退化模型被可靠性研究学者提出和研究.但是我们发现许多显示设备,例如等离子显示屏和有机发光二极管等设备在退化过程中表现出两阶段的退化特性.因此变点问题成为一个非常有意义的研究方向.针对该类型的数据我们考虑使用维纳过程进行建模,同时我们假定维纳过程中的漂移函数为关于时间的两段线性函
蝙蝠(翼手目)作为唯一具有真正飞行能力的哺乳动物,其种类在现存哺乳动物中高达20%。为适应飞行,蝙蝠进化出一系列与飞行相关的形态结构,如一对被宽大翼膜包被且特异性延长的“翅膀”、脚踝处支撑尾膜的距。与鸟类带羽毛的翅膀不同的是,蝙蝠的翼手由延长的前臂和第二到五指支撑着薄而有韧性的翼膜形成。蝙蝠近后端四个指特异性延长,而拇指(前肢第一指)和脚趾依旧是短的且长度和宽度相近。蝙蝠前肢形态的巨大改变即翼手的