任意凸四边形网格上各向异性扩散问题的Q1型有限体积元法研究

来源 :中国工程物理研究院 | 被引量 : 1次 | 上传用户:zlq5626
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本博士学位论文主要研究Q1型有限体积元方法(FVEM)在任意凸四边形网格条件下关于各向异性扩散问题的强制性.首先,使用梯形求积公式逼近经典Q-FVEM双线性泛函中线积分,得到一个所谓的修正Q1-FVEM(mQ1-FVEM),并进一步研究这个格式的强制性.基于单元刚度矩阵正定的充分必要条件,在任意凸四边形网格上,得到mQ1-FVEM强制性成立的一个充分条件.这个充分条件包含一些已有的标准网格,例如传统的h1+γ-平行四边形网格和一些梯形网格.更有趣的是这个条件有一个显式表达式并且根据这个表达式对任意扩散张量和任意网格尺寸h>0都很容易判断强制性是否成立,这在实际应用中相当具有吸引力.最后,在不需要传统的h1+γ-平行四边形网格假设条件下,严格证明了mQ1-FVEM格式的H1误差估计.因为在对偶网格的边界上Q1型有限元基函数的梯度一般是有理函数,这导致mQ1-FVEM和Q1-FVEM格式是不相同的,并且mQ1-FVEM格式强制性分析方法并不能直接推广到经典Q1-FVEM格式.因此,需要研究一种新的途径在任意凸四边形网格上研究Q1-FVEM强制性.根据单元双线性泛函(?)的表达式,通过一个技巧将原始的单元刚度矩阵转换成一个新的3 × 3单元矩阵.基于这个新单元矩阵正定的充要条件,得到了经典Q1-FVEM强制性成立的一个充分条件.本文发现这个允分条件包含传统的h1+γ-平行网边形条件.同样地,这个条件也有一个显式表达式并且根椐这个表达式对任意扩散张量和网格尺寸h>0都很容易判断强制性是否成立.最后,在不需要h1+γ-平行四边形网格条件下,经典Q1-FVEM的H1误差估计是平凡的.在凸四边形网格上关于各向异性扩散问题,使用一个特殊的求积公式逼近经典Q1-FVEM格式中的线积分,本文提出了一种新的Q1-FVEM格式(sQ1-FVEM).在拟正则网格条件下且不需要网格尺寸充分小,严格证明了sQ1-FVEM格式的强制性,即sQ1-FVEM格式在拟正则网格条件下是无条件稳定的.基于sQ1-FVEM格式强制性结果,本文也提供了一种新的方式证明经典Q1-FVEM格式在正则网格条件下强制性成立.并且在正则网格条件下,本文严格证明了 sQ1-FVEM格式和经典Q1-FVEM格式的H1误差估计.关于sQ1-FVEM格式的L2误差估计,本文构造的一个反例表明通过Aubin-Nitschc技巧证明sQ1-FVEM格式的L2误差估计时额外的网格限制条件(例如h1+γ平行四边形网格)仍然需要.最后,考虑mQ1-FVEM格式的应用.使用mQ1-FVEM格式处理节点未知量,本文构造了各向异性扩散问题的一个稳定的九点格式(NPS-mQ1)并对其进行了理论分析.鉴于九点格式易编码,它已被广泛应用求解一些辐射流体动力学代码,例如LARED-I和MARED[22,69].当NPS-mQ1格式应用到此类问题时,仅需要增加mQ1-FVEM格式求解节点未知量的程序,这很容易实现.另外,基于mQ1-FVEM格式的理论结果和离散泛函分析[81,94],在弱几何条件下可以得到稳定性分析和H1误差估计.与已有的一些中心型和杂交型格式[51]相比较,本文提出的NPS-mQ1格式能克服所谓的数值热障问题.
其他文献
随着量子计算和量子信息技术的发展,越来越多的人开始关注对微观量子态的操作和调控问题。光子作为电磁场量子化的能量实体,与电磁环境没有直接的相互作用,已经成为量子信息载体的理想候选者。相比于其它有质量的微观粒子,光子在应用上具有速度快、容量大、抗干扰能力强和保密性好等优点。近些年,通过利用光子与物质之间的有效相互作用来实现对光子态的操控已经成为量子物理中的重要发展方向。研究表明,将量子发射器放入到受限
相对传统TiAl合金,高Nb-TiAl中部分Ti原子位置或Al原子位置被Nb原子取代,表现出优异强度、抗氧化及抗蠕变性能,应用范围比传统TiAl合金更宽,有望在中等应力和温度区间取代镍基高温合金。同时,由于Nb原子添加,引起TiAl合金中γ/γ+α2和α2/γ相界发生改变,导致高Nb-TiAl合金片层晶团晶界上形成各种不同析出相。另外,热处理工艺及Nb含量不同,也会引起新相析出,这些析出相影响着高
本文主要研究非齐次不可压Navier-Stokes方程密度补丁的全局正则性问题和两类流体方程的适定性问题.全文共分五章,具体如下:第一章为引言部分.主要综述所研究问题的物理背景和相关函数空间的定义,本论文的主要结果及其创新点.第二章研究不可压缩Navier-Stokes-Fokker-Planck方程的Cauchy问题.首先利用Garlakin逼近方法构造系统的近似解,然后利用Littlewood
本文主要研究几类流体力学方程的适定性与周期解问题。在第一章中,我们简要介绍了流体力学方程的适定性问题与周期解问题的研究现状,并概括了本文的主要研究内容与研究意义。在第二章中,我们收集了一些和本文研究内容相关的一些预备知识。在第三章中,基于高低频分解的方法,通过对液晶方程的未知变量的低频部分应用Fourier变换,对高频部分应用加权能量估计,我们证明了对于带外力的液晶方程当外力的某些加权Sobole
在铀铌合金体系中,U-5.4 wt.%Nb(U-5.4Nb)合金因兼具优异的抗腐蚀性能和特殊的力学行为而备受关注。目前,U-5.4Nb合金主要采用压力加工成型和铸造成型获得目标样品。在这两种成型方法中,金属液的凝固过程对样品的铸态组织有重要影响。因此了解、掌握凝固过程变量与合金铸态组织之间的对应关系有助于准确定位关键影响因素,优化控制参数。基于此,本文利用凝固冷却速率可综合反映温度梯度和凝固速率的
在过去的几十年中,原子分子的碰撞碎裂动力学吸引了大量理论和实验科学家的关注,许多开创性的成果被发表。然而这些工作主要局限于简单原子体系(如氢、氦、碱金属原子等)和简单分子体系(氢气、氧气、氮气、一氧化碳等),对于复杂的多电子原子和包含三个以及三个以上原子的分子体系,由于多体系统薛定谔方程无法精确求解,目前还缺乏可靠的理论来预测碰撞碎裂反应的截面信息以及解释碎裂过程物理机制。为此,我们需要发展新的理
强场物理是研究强激光场与物质相互作用的一门科学。超强超短脉冲激光技术,不仅为人们提供了探索微观物质结构及其运动规律的重要手段,而且还使得人们可以对原子分子中的电子动力学过程进行实时探测和操控。过去几十年,强激光场与原子相互作用产生的一系列新颖的高阶非线性物理现象:比如多光子电离(MPI)、阈上电离(ATI)、高阶阈上电离(HATI)、非次序双电离(NSDI)、高次谐波产生(HHG)等等,大大加深了
本文主要研究的是一类满足齐次Carleson-Sj(?)lin条件的Fourier积分算子的局部光滑估计,其中齐次Carleson-Sj(?)lin条件包含两部分内容,一部分是非退化条件,另一部分是曲率条件.通过建立变系数情形下的平方函数不等式,我们改进了 Mockenhaupt-Seeger-Sogge[41,J.Amer.Math.Soc.6(1):65-130,1993.]中的结果.本文证明
准一维纳米结构、二维层状结构等低维组织的发现,为凝聚态物理研究打开了一个全新的科学研究前沿,并推动了微小尺度下材料科学的基础研究。相比于传统材料体系,低维体系在新一代电子器件设计以及器件小型化方面具有无与伦比的优势。它们的奇异物质,将革新微型半导体器件的设计以及生产,对现代半导体技术的进步具有积极且深远的意义。应变调控指通过施加应变(结构变形)来调制材料的电、声子性质。对于低维材料,应变调控具有重
二维Benjanmin方程(?)是在一维Benjamin模型的基础上带有了垂直于传播方向的扰动.相对于二维KP方程,二维Benjamin方程没有伸缩不变性,没有对频率一致的t-1衰减估计,而该t-1衰减对于建立后续的双线性估计是起关键作用的.这使得其初值问题的适定性很具有挑战性.本文主要研究三阶和五阶二维Benjamin方程的适定性主要结果有1.建立了二维Benjamin方程(?)在能量空间E1中