若干流动问题的MAC方法及其稳定性和超收敛性分析

来源 :山东大学 | 被引量 : 0次 | 上传用户:qingfeng44
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
Mark and Cell(MAC)方法是交错网格上的一种有限差分方法,是公认的处理Stokes问题和Navier-Stokes问题中最简单且最高效的方法之一,此观点在文献[35]和[40]中都有提到。MAC方法是在20世纪60年代由Lebedev[48]和Daly等人[21]提出的,直到1992年,才由Xicolaides[60]以及Nicolaides和Wu[61]通过将MAC方法转化为有限体积方法得到其理论分析,但是介于对有限体积方法的数学分析工具也很少,所以当时对MAC方法的理论分析并不成熟;1996年,Girault和Lopez[35]提出利用特定的求积公式将MAC方法转化为混合有限元方法,以此来得到方法的理论分析;1998年,Han和Wu[40]也将MAC方法转化为一种新的混合有限元方法,从而得出方法的理论分析;随后在2008年,Kanschat[63]应用特定的积分将MAC方法等价为一阶局部间断有限元方法;紧接着,Minev[58]受到启发利用矩形网格上的一阶Nedelec元得到MAC方法。需要注意的是以上所提到的文献都只证明MAC方法在均匀网格中的速度(H1范数下)和压力(L2范数下)能达到一阶收敛,但是Nicolaides在文献[60]中指出数值计算结果表明MAC方法是可以达到二阶收敛的。2015年Li和Sun在文献[50]中建立了 Stokes方程的MAC格式,文章中证明了速度的稳定性,并得到了 MAC方法在非均匀网格中得到的速度是二阶收敛的,据我们所知这是第一篇证明MAC方法二阶收敛的文章,但是文章中存在几个疑点并没有得到解决,比如在文章中,速度的稳定是基于压力的稳定,但是压力的稳定性在文中并没有给出。随后在2017年,Rui和Li[71]通过建立离散的Ladyzhenskaya-Babuska-Brezzi(LBB)条件证明了压力的稳定性,他们证明了速度和压力在L2范数下均能达到二阶的收敛结果,还得到了离散H1范数下某些项是二阶收敛的,其他项在均匀网格中也能达到二阶收敛。我们这篇文章就是遵循着文献[71]的主要思想,先建立离散的LBB条件来分析稳定性,再利用速度的辅助函数来分析得到速度和压力的收敛性。在本篇论文中,我们主要讨论MAC方法应用于Darcy-Stokes-Brinkman方程,带阻尼项的Stokes方程,以及耦合的Stokes-Darcy方程,并分析其稳定性和收敛性。文章中主要内容如下:第一章中,我们给出整篇文章中通用的一些符号和定义,包括区域剖分,离散差商,离散内积和范数的定义,方便之后的应用。第二章中,我们讨论了 MAC方法用于Darcy-Stokes-Brinkman问题,并分析其在非均匀网格中的稳定性和收敛性。Darcy-Stokes-Brinkman问题,简称Brinkman问题,其模型会随着扰动参数ε的变化在Darcy流和Stokes流之间转化,所以用于解决Stokes问题的MAC方法自然可用于Brinkman问题。在本章中,我们先建立了对应的离散LBB条件,通过将MAC格式转化为混合有限元的变分格式得到了速度和压力的稳定性。收敛性分析中,我们引入基于速度和离散参数的辅助函数,先得到了 MAC格式得到的近似速度在范数‖· ‖ε下是二阶收敛与于辅助速度函数的,再利用辅助函数的定义得到MAC方法的误差估计结论:1、非均匀网格中,速度和压力在离散L2范数下可达到二阶收敛;2、当ε不趋于零时,速度的x方向分量在x方向的差商在L2范数下是可以达到二阶收敛的,同样地,速度y方向分量在y方向的差商也是二阶收敛的;3、当ε不趋于零时,速度x方向分量在y方向的差商的L2范数在非均匀网格中只能达到一阶收敛,在均匀网格中可到达1.5阶,若再去除边界项,则可达到二阶收敛,速度y方向分量在x方向的差商也可得到类似结论。本章的最后,我们会用数值实验来验证理论结果。第三章中,我们考虑带阻尼项的Stokes方程,应用MAC方法来做问题在非均匀网格中的离散。阻尼项、也称作Forchheimer项,是非线性项,对这一项的处理是本章的关键点,我们会参照文献[36,70]中对Forchheimer项的技巧来处理这一项。为了分析方法的稳定性,我们还是先建立对应的离散的LBB条件,再将MAC格式转化为混合有限元的变分格式并得到速度和压力的稳定性。在误差分析中,同样引入了关于速度的辅助函数,先得到了速度和辅助函数之间的二阶收敛,再利用辅助函数的定义得到MAC方法得到的速度与真实速度之间的误差估计。文章的主要结论如下:非均匀网格中,速度和压力在L2范数下都是二阶收敛的;速度在离散H1范数下,x方向分量在x方向的差商和y方向分量在y方向的差商是二阶收敛的;而x方向分量在y方向的差商和y方向分量在x方向的差商只能达到一阶收敛,在均匀网格中,可以达到1.5阶收敛,若再去除边界项,其均能达到二阶收敛。本章的最后也给出了数值实验,其中我们利用Picard迭代来解决非线性部分,验证了理论分析的正确性。第四章中,我们建立了耦合的Stokes-Dary问题在非均匀交错网格中的MAC格式。据我们所知,MAC方法是用于求解Stokes方程最简单有效的方法之一,而块中心差分方法用于解Darcy方程和Darcy-Forchheimer方程时也可达到二阶收敛(分别参见文献[65,83]),并且MAC方法和块中心差分方法都是交错网格上的差分方法,所以块中心有限差分方法也可以被看作为是用于解决多孔介质流的MAC方法,但是解决耦合问题的关键困难在于:在交界面上,Stokes流和Darcy流有不同的正则性质,并且切向的速度也不连续。本章的主要内容如下:首先我们建立耦合问题在交错网格上的MAC差分格式,其中交界面上我们应用有限体积积分来得到对应的差分格式;然后,为了得到整个区域上速度和压力的稳定性,我们建立了对应的离散LBB条件;为了分析收敛性,我们引入速度的辅助函数,且得到结论:1、在非均匀网格中,Stokes流和Darcy流的速度和压力在L2范数下均可达到二阶收敛,2、对于Stokes流的速度,在非均匀网格中,x方向的分量在x方向的差商和y方向的分量在y方向的差商在L2范数下均能达到二阶收敛,3、而x方向的分量在y方向的差商和y方向的分量在x方向的差商在非均匀网格中只能达到一阶收敛,在均匀网格中能达到1.5阶收敛,若再去除边界项,则能达到二阶收敛;最后,还是会给出数值实验来验证理论分析的可靠性。最后在第五章中,我们对整篇论文做一个总结,并且对未来的工作展开做了简单的规划。
其他文献
刀具切削温度过高导致刀具切削性能下降,引发刀发生热变形,加速刀具磨损,降低加工表面质量和精度,增加加工成本,降低加工效率。刀具前刀面与切屑之间的接触热阻影响刀-屑之间的热传导行为,对切削热在刀具、切屑等区域的分配起到重要作用。TiAlN、CrAlN、TiAlCrN、CrAlSiN等刀具涂层可改变刀-屑接触特征,改变刀-屑接触热阻,进而对刀具切削温度产生影响。本文以TiAlN涂层刀具在金属切削过程中
包合水合物,因其在能源储存和器件应用等方面的巨大潜力,有望改善人类社会的能源危机,因此已成为能源和环境领域的一类重要材料。近年来,人们对于包合水合物堆积的笼状结构所发生的包合现象的本质已经进行了大量的研究。此外,一些包合水合物表现出独特的磁性、强的离子导电性、高于室温的熔点、不同寻常的结构转换等显著的物理化学性质也将导致包合水合物转化为极具应用前景的水分子基的新型功能材料。尽管包合水合物展现了广阔
自从上个世纪50年代以硅基半导体为代表的第一代半导体出现,电子产品变得愈发微型化和集成化。而随着器件尺寸的不断减小,传统的硅基半导体材料的性能即将达到物理和技术的极限。石墨烯的成功剥离使得二维材料进入了人们的视野,短短十几年里,许多新型二维材料被相继发现并合成,例如硅烯、六角氮化硼、过渡金属二硫化物(TMDs)和黑磷等。甚至原本被理论预言无法存在的二维磁性半导体也已经合成,如Cr2Ge2Te6、C
背景再生障碍性贫血(再障aplastic anemia,AA)是由多种原因引起的一种骨髓衰竭综合征。骨髓血细胞生成减少导致粒细胞、网织红细胞、单核细胞以及血小板均减少,临床表现为不同程度的贫血、出血和感染。目前公认,AA是T细胞介导的自身免疫性疾病。AA患者体内存在辅助性T细胞0(helper T lymphocyte 0,Th0)向辅助性T细胞1(helper T lymphocyte 1,Th
MicroRNA(miRNA)是一类内源性、长度约为23个核苷酸的非编码单链RNA分子。从二十世纪九十年代首次被发现,miRNA的研究开辟了分子生物学的新领域,也是后基因组时代里生命科学研究的前沿。随着下一代测序(NGS)和生物实验技术的快速发展,miRNA的异常表达已经不断被证实与神经系统疾病、心血管疾病、恶性肿瘤等人类复杂疾病的发生和发展密切相关。此外,随着人类微生物组计划进入第二阶段,作为人
随着电子产品日益向轻薄和智能化的方向发展,对集成度的要求以及对各种电子元器件的微型化要求也相应地提高,纳米材料具有优异的光、电、磁、力学等性质,成为构建功能器件的理想基元。其中,低维碳纳米结构如石墨烯、碳纳米管等在集成电路方面具有良好的应用前景,成为未来计算处理器件的理想候选材料。但是,低维碳材料作为电子器件材料存在诸多缺点,例如,石墨烯具有零带隙性质、碳纳米管的电子性质随手性变化等,制约了该材料
手性纳米结构是自然界中最重要的结构之一,它与生命起源和基因传递等都有着非常重要的关联。手性现象可以出现在不同的尺度,包括微观纳米尺度到宏观的星系尺度。目前人们可以通过在有机材料中引入手性基团、手性或者非手性分子自组装等方式获得具有长程有序的手性超分子结构。而螺旋结构中的电子涡旋束可以产生明显的轨道角动量(OAM)。由于手性螺旋结构产生的OAM,具有手性螺旋结构的手性超分子材料具有自旋选择效应。一般
量子相干性是量子力学最本质的特征,是完成各类量子信息处理任务的必要资源,它已被广泛应用于量子计算、量子通信、量子热力学、量子精密测量以及量子生物学等众多领域。随着量子信息科学的发展,量化量子相干性已成为一个亟待解决的任务,自2014年首个相干度量框架被提出之后,迅速成为了国内外的热门研究课题。量化量子相干性领域的核心任务之一是寻找有效的相干度量并基于相干度量探索相干性与其他量子特性如量子关联之间的
在过去的几十年中,越来越多的扩散过程被证实不满足Fickian定律,例如在生物细胞中的信号传导,神经细胞中的反常的电扩散,动物的觅食行为,黏弹性和粘塑性流动以及地下水中的溶质的迁移等。而分数阶偏微分方程在描述反常扩散方面有着非常重要的应用,所以近些年来,分数阶微分方程受到人们的广泛关注。分数阶微分方程已被成功应用于各学科和工程领域,例如,流体力学,材料力学,物理学,金融学,化学等许多领域。对于布朗
研究背景高血压肾病是导致终末期肾病的主要因素之一,早期评估肾损害危险因素,对高血压进行综合防治,从而延缓高血压肾病的进展,对阻止终末期肾病的发生具有重要的临床意义。流行病学研究表明,高同型半胱氨酸血症(hyperhomocysteinemia,HHcy)和低叶酸(folic acid,FA)现象在中国高血压人群中普遍存在,而HHcy可协同增加高血压患者动脉粥样硬化、脑卒中等疾病的发生。但关于HHc