关于粘性不可压缩流动问题和Maxwell方程组的数值离散方法研究

来源 :四川师范大学 | 被引量 : 0次 | 上传用户:aptx4866
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
关于粘性不可压缩流动问题的数值离散方法研究一直是计算数学研究的热点.Navier-Stokes方程是粘性不可压缩流体问题的基本方程,而Stokes方程是Navier-Stokes方程的定常形式和线性化,对它的数值离散方法研究具有典型性和普遍性意义.Brinkman方程是描述粘性不可压缩流体在渗透系数快速变化的复杂多孔介质中的流动方程,对它的数值离散方法的研究也非常重要.混合有限元方法是研究粘性不可压缩流动问题的一种常用的数值离散方法.由于传统的混合有限元方法需要有限元空间满足inf-sup条件,这个条件限制了工程上非常好用的低阶元的应用.除此,传统混合有限元方法对网格剖分单元的形状要求比较严格,一般只能是三角形或四边形(n=2)单元.这对实际应用中有限元空间逼近满足稳定性条件和复杂区域边界问题求解上带来困难.粘性不可压缩流动问题的数值解严格满足不可压缩条件对解的稳定性、收敛性具有重要意义,而传统的混合有限元方法很难构造无散的有限元格式.为了克服传统混合法遇到的困难,近年来,对粘性不可压缩流动问题的数值离散方法研究转向于非标准的有限元方法的研究,如间断有限元方法、杂交间断有限元方法、弱Galerkin有限元方法等.这些方法的优点在于网格剖分灵活、容易满足稳定性条件、易于构造满足无散的有限元格式.本文的第一部分工作主要研究两类粘性不可压缩流动问题-Stokes方程和Brinkman方程的弱Galerkin方法,分别构造了全局无散的弱Galerkin有限元离散格式,证明了该离散格式的稳定性,得到了与粘性系数一致的误差估计,并用数值算例进行验证.本文的第二部分工作主要面向Maxwell方程组.Maxwell方程组是电磁学的基本方程组,对Maxwell方程组的数值计算一直是计算电磁学的热点问题.时域有限差分方法(Finite-Difference Time-Domain method,简称为FDTD)是最受欢迎的数值方法之一.但它是条件稳定的,即时间步长和空间步长需要满足Courant-Friedrichs-Lewy(CFL)条件.因此当解决高维问题或实际问题要求空间步长很小时,时域有限差分方法的计算量是巨大的,有时候甚至不可实现.因此各种无条件稳定的差分方法被相继提出,如方向交替时域有限差分方法(Alternating Direction Implicit Finite-Difference Time-Domain Method,简记为ADI-FDTD)、分裂时域有限差分方法(Split Finite-Difference Time-Domain method,简记为S-FDTD)等.另一方面,电磁场在传播过程中的能量变化也是计算电磁学研究的一个重点.因此,对无损介质中的Maxwell模型构建稳定的、能量守恒的FDTD方法是具有实际意义的.我们的工作是研究二维Maxwell方程组的时域有限差分方法.其一是针对二维Maxwell方程组构造了两种不同精度的分裂时域有限差分格式,对格式的稳定性进行分析.其二是针对无损介质中Maxwell方程组一种时间4阶的ADI-FDTD格式的能量进行分析.本文的主要结构如下:在第一章中,介绍粘性不可压缩流动问题及Maxwell方程的主要背景以及相关数值离散方法概述.在第二章中,简要介绍本文用到的Sobolev空间、有限元方法分析过程中重要的公式以及重要的定义.在第三章中,针对Stokes方程提出一种改进的弱Galerkin有限元方法.该方法的速度和压力有限元空间在单元内部选择P6)()/P6-1)()(6)≥1),单元边界上速度函数取值于0(??)连续的函数空间.本文证明了改进的方法同样是稳定的,且保持无散的性质,能得到与粘性系数一致的误差估计.最后数值试验的结果展示改进后的方法比原来的方法具有更少的自由度,执行效率更高.在第四章中,针对定常的Brinkman方程提出一种新的弱Galerkin有限元格式.该格式基于一种新的变分格式,并且弱Galerkin有限元格式的数值解满足不可压缩条件,很好的保持了原方程解的物理性质.文中证明了格式的稳定性,推导了一系列的误差,得到与雷诺数无关的误差估计.最后给出数值算例验证理论结果.在第五章中,首先针对磁导率不为零的二维麦克斯韦方程组进行研究.用算子分裂的方法结合插值公式对其构造了两个不同精度的差分格式TS-FDTDI和TS-FDTDII,并用Fourier方法对格式TS-FDTDII的稳定性进行分析.其次我们对无损介质中的二维Maxwell方程组(TE问题)的一种时间四阶ADI-FDTD方法的能量进行分析,推导了该格式的数值能量等式.对能量等式中的两个扰动项进行分析,得到了该格式的能量是渐近守恒的.最后用数值算例验证了理论分析结果.
其他文献
目的:探讨不同浓度七氟醚(Sevo)预处理对小鼠心肌缺血/再灌注(MI/R)损伤后血浆脂联素(APN)水平的影响。方法:将120只C57野生型小鼠随机分为5组:空白对照组(control)、1%七氟醚预处理组(S1)、2%七氟醚预处理组(S2)、3%七氟醚预处理组(S3)、4%七氟醚预处理组(S4)。使用不同浓度七氟醚处理各组小鼠,建立MI/R模型。分别于术后0h、6h、12h、24h使用ELIS
眼肌型重症肌无力(ocular myasthenia gravis,OMG)是一种累及神经肌肉接头突触后膜的自身免疫性疾病。OMG波动性的临床特征易导致较高的漏诊率和误诊率,严重影响患者的生活质量。因此,检测方法在早期确诊OMG中发挥重要作用。本文围绕临床检查、新斯的明试验、自身抗体检测、电生理检测等方面对OMG检测方法作一综述,以期为OMG的早期确诊提供有益的指导。
针对根据陆态网络区域网联测资料来构建速度场复测周期长、不能连续地观测到地壳区域的形变与应变变化的问题,该文利用GAMIT/GLOBK软件解算2015—2019年云南区域45个连续运行GNSS基准站观测数据,得到云南区域在ITRF2014框架下的速度场,进一步获取云南区域整体旋转无基准的速度场,然后利用克里金插值方法对站点的速度进行插值,再利用最小二乘配置法求取云南区域应变场,实时动态地监测云南地区
重症肌无力(myasthenia gravis,MG)是国际公认的难治性疾病。从20世纪70年代首次出现中西医结合治疗MG的报道以来,迄今已走过近50年历程,在本病的治疗思路、治疗方法、机制研究等方面取得了一定进展,促进了治疗水平的提高,但也存在诸多不足。本文在复习文献基础上对近50年中西医结合治疗MG概况进行回顾与反思,以期对进一步提高本病临床治疗及研究水平有所裨益。1中西医结合治疗MG
期刊
在上海市社区治理和服务创新工作持续推进的过程中,凸显了一些发展中的问题。比如:社区工作者对社区要素的了解、掌握、使用程度参差不齐;公共服务没能做到"精准投放",社区居民的满意度和获得感有待提高;社区动员中青年社区居民参与社区公共事务和公益事业的能力有限;社区工作者的稳定性、专业性和服务效能仍有待提升等。为解决上述问题,上海市静
期刊
本论文主要是通过模理论方法刻画w-凝聚环.第一章回顾了w-理论中的一些基本概念,例如w-模,w-包络,w-Noether环,w-凝聚环,w-平坦模.第二章研究了w-模类和w-平坦模类的盖包性质.首先证明了 w-模类既是盖类又是包类,其次得到了w-平坦模类是盖类,最后给出了 VN-正则环的一些新的刻画.第三章引入并研究了关于全体平坦模类F的w-Mittag-Leffler模类w-F-ML.证明了R是
设R(?)T是整环的w-linked扩张.在第一章,我们引入了 T-模M的wR-平坦维数(wR-fdTM)的概念,进而引入了环T的wR-弱整体维数(wR-w.gl.dim(T))的概念.利用环T的wR-弱整体维数可以来刻画域及PwRMD.具体地,设T是R的w-linked扩张,则T为域当且仅当wR-w.gl.dim(T)=0,以及T为PwRMD当且仅当wR-w.gl.dim(T)≤1.设R是有单位
设R是交换Noether环,且记g(R)={M是完全自反R模},ε(R)={M是有限生成R-模|对每个极大理想m,depth(Mm)≥depth(Rm)}.Noether环R称为G-正则的,如果R在每个极大理想m处的局部化Rm是G-正则局部环,即Rm上每个完全自反模是自由的;环同态φ:R→S称为G-消失的,如果Ext1R(g(R),S)=0且Tor1R(g(R),S)=0.本文给出G-正则环与G-
术后认知功能障碍(Postoperative cognitive dysfunction,POCD)是指麻醉和手术后患者出现精神状态改变、人格的改变及记忆力下降等,严重时可能导致认知及社会行为能力的降低,是老龄患者术后常出现的并发症之一,而且远期与痴呆的发生有密切联系。造成POCD的因素尚未完全了解,但许多研究已证实全身麻醉药物如七氟醚、手术创伤等多种因素都与POCD的发生有关。POCD的发生延长