定常不可压缩流的若干稳定化方法研究

来源 :新疆大学 | 被引量 : 4次 | 上传用户:lyun888
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
定常不可压缩流可以近似地看做常数的一种流体,它刻划着一些流体的运动规律,如海洋流动、大气运动以及透平机械内部流动等.特别地,它的研究对人们认识和控制湍流至关重要.描述这种流体的控制方程主要是不可压缩Navier–Stokes方程.由于人们对非线性现象本质认识有限,因而数值模拟就成为一种十分重要的研究手段.但直接数值模拟Navier–Stokes方程有一个很大的困难就是巨大的解题规模与有限的计算资源及算法稳定性之间的矛盾.因此,构造和研究具有良好稳定性和收敛性的高效算法就显得尤为重要.本文基于局部高斯积分公式,研究了定常不可压缩流的若干稳定化算法,主要做了下面的工作:一、基于局部高斯积分公式,提出了求解定常Navier–Stokes方程的一种新亏量校正格式.它主要用一个双线性项来代替原先的人工粘性稳定化项.该双线性项是由二次多项式的精确高斯积分公式和线性多项式的精确高斯积分公式做差得到.不同于通常的亏量校正方法,新格式对参数依赖较小且不依赖于网格.并且通过数值实验发现,用我们提出方法得到的解比通常的亏量校正方法得到的解更为精确一些.二、提出了基于三种校正算法的两水平二次等阶稳定化有限元方法.采用的P2P2元不满足离散的inf-sup条件,故而用基于局部高斯积分公式的压力投影技巧来稳定.得到了三种算法: Stokes校正、Newton校正和Oseen校正.数值理论分析和计算结果发现Stokes和Newton校正对大粘性问题处理较好,而Oseen校正求解小粘性系数问题则有一定的优越性.接着,建立了两水平亏量校正稳定化方法.在处理大雷诺数问题的时候,采用亏量校正方法来求解,但是往往需要很多的计算时间.而两水平方法是一种高效的算法,它能节约大量的CPU时间.双赢的组合对求解Navier–Stokes方程的大雷诺数问题有了很大的帮助.接着,基于P1-P1高斯积分稳定化方法,给出了三种算法: m次亏量步采用Oseen迭代,1次校正步采用Oseen迭代(m-Oseen-1-Oseen); m次亏量步采用Oseen迭代,1次校正步采用Simple迭代(m-Oseen-1-Simple); m次亏量步采用Oseen迭代,1次校正步采用Newton迭代(m-Oseen-1-Newton).三、采用L2投影方法,对基于P1-P1元求解定常Navier–Stokes方程的稳定化非协调有限元方法和稳定化有限体积方法,分别建立了它们的超收敛结果.该方法的主要思想是把数值结果投影到另一个不同的较粗网格空间上.经过后处理,用两种网格尺度的差别来达到超收敛结果.数值例子验证了数值理论的准确性.四、基于用局部高斯积分来做稳定化的P1-P1元,提出了解Stokes特征值问题的一种两水平有限元方法.该方法求得的解和一般的稳定化有限元求得的解具有相同的收敛阶,且我们的方法能够节省大量的计算时间.数值试验验证了理论结果.进而,给出了一些求解Stokes特征值问题的基于最低等阶元的稳定化方法,包括加罚方法、正则化方法、丰富多尺度方法、局部高斯积分方法以及非协调局部高斯积分方法.接着,对这些稳定化混合有限元方法在数值求解Stokes特征值问题时的表现做了比较和分析.最后,指出了相比其他方法,用非协调局部高斯积分方法求解Stokes特征值问题具有较好的稳定性和误差结果.
其他文献
代数图论是图论学科的重要研究领域之一,主要运用代数方法来解决图论问题.代数图论有三个主要分支,分别为图与线性代数,图与群论,图不变量.其中图与线性代数的研究核心是图的谱理论.图谱理论是通过研究与图相关的矩阵(邻接矩阵,拉普拉斯矩阵,无符号拉普拉斯矩阵,正规拉普拉斯矩阵和距离矩阵等)的谱的性质来研究图的组合性质.图谱理论研究兴起于20世纪50-60年代,近20-30年来得到迅猛的发展,它是图论,组合
Grobner 基理论是由 Buchberger,Shirshov 和 Bergman 独立引进的.Buchberger 创建的交换代数的Grobner基理论为解决交换代数中的约化问题提供了非常有效的方法.Bergman把Buchberger的理论推广到结合代数上.在李代数上的类似理论由Shirshov创建.后来,Bokut证明了 Buchberger和Bergman的Grobner基理论其实是S
本文考虑了与一个半有限von Neumann代数对应的非交换Lorentz型空间,包含四部分内容.第一部分介绍了文章的研究背景,非交换积分理论的基础知识以及本文的主要结果.第二部分研究了非交换弱Lp空间,给出了非交换弱Lp空间对偶空间的具体形式.在这一部分中还讨论了非交换弱Lp空间上的紧左(右)乘算子.第三部分研究了非交换加权Lorentz空间的对偶空间,得到了类似于经典的加权Lorentz空间对
当今,计算已成为继理论和实验之后的第三种不可或缺的科学研究方法。并且在许多情况下,由于科学计算不受外部因素和实验器材影响的灵活性,它能够最大程度以最小的代价获得与理论和实验相当的结果。这个使得计算在当今科学领域占有非同一般的统治地位。反应扩散方程是一类重要的偏微分方程,在物理、生物、材料以及社会科学中都有广泛的应用。并且,反应扩散方程有着比其他偏微分方程更加好的特性,例如极值原理、比较原理、不变集
设G是一个简单图,M=M(G)是按照某种规定所定义的与G相联系的图矩阵,把利用M的特征值来刻画图G的组合结构的理论称为图谱理论(M-谱理论).定义det(xI-M)为图G的M-特征多项式,其中I为单位矩阵.M-特征多项式的特征根称为图G的M-特征值,由G的所有M-特征值构成的多重集称为M-谱,简记为SpecM(G).图G的最大M-特征值称为M-谱半径.关于图矩阵M具有相同谱的图称为M-同谱图,与G
本文运用算子理论,空间理论及调和分析的方法,研究了非交换Hardy-Lorentz空间的分解问题,非交换Lorentz空间的极大不等式以及一些关于非交换Lp-空间中τ-可测算子的不等式.本文共分为四章,结构如下:第一章介绍了文章的国内外研究现状,主要结论以及常用的符号与定义.第二章研究了非交换Hardy-Lorentz空间的Szego分解及内外型分解.第三章包括两节.在第一节中我们讨论了Hardy
这篇论文分为两部分,分别介绍了有关图中的哈密顿圈和图的列表线性荫度的一些研究成果。第一部分由三章组成。在第一章引言中,我们给出了图的有关定义及概念并介绍了图的哈密顿圈的研究背景。第二章中我们研究了中间图的补图的哈密顿性。图G的中间图M(G)的顶点集为V (G)∪E(G),两个点x和y相邻当且仅当x和y中至少有一个是G的一条边,并且它们在G中相邻或关联。我们定义图G的中间图M(G)的补图为M(G),
分数阶微积分方程是经典微积分方程自然的数学推广,具有深刻的物理背景和丰富的理论内涵,在物理、生物、化学等多个学科领域具有广泛的应用.对分数阶微积分方程的研究,不仅具有重要的学术研究价值,而且具有广阔的工程应用前景.目前,分数阶微积分方程的研究己成为国际上的一个热点研究课题.本文主要研究分数阶偏微分方程的几种近似算法(包括近似解析算法和数值算法).本文的研究内容主要包括三个部分:第一部分研究同伦分析
本论文主要讨论了与具有反霍尔德类势的薛定谔算子相关的一些问题。本论文研究了四类问题:具有反霍尔德类势的薛定谔算子的Riesz变换分别与加权Lipschitz函数和加权BMO函数构成的交换子的加权Lp有界性;与具有反霍尔德类势的薛定谔算子相关的BLOL空间;一个具有反霍尔德类势的抛物型薛定谔算子的Lp有界性;一类具有反霍尔德类势的一致抛物型算子在加权Lp和Morrey空间上的有界性。全文共分五章。第
图谱主要通过研究与图相关的矩阵(邻接矩阵,拉普拉斯矩阵,无符号拉普拉斯矩阵等)的谱(特征值的多重集)的性质来研究图的其它组合性质.图谱是图论,组合,矩阵论及代数理论的一个交叉领域,是代数图论的一个研究分支,它在图论,物理,量子化学,计算机科学,互联网技术等方面有着广泛的应用.谱确定问题是图谱理论中一个著名且困难的问题,本文就此问题展开了研究.第一章首先介绍了图谱理论的研究背景及应用,然后引入了本文