论文部分内容阅读
Navier-Stokes方程是流体力学的基本方程,它在气象学、海洋学、生命科学、地球物理等涉及流体运动的各学科中有着广泛的应用。此外,由物理、化学等学科提出的模拟燃烧、爆炸、湍流、多相流、复杂流体及弱可压流体也均基于N-S方程。因此,研究N-S方程的数值模拟,对各种流体现象的科学探知和相关实践有着重要的意义和影响。
由于Navier-Stokes方程为一个非线性的鞍点问题,在高Reynolds数流情形下,标准Galerkin有限元方法会导致伪数值振荡。迎风有限元方法,因其在算法结构上表征了流体“上游”决定“下游”的流动性态,使得它能够有效地消除对流占优-扩散奇异摄动问题的非物理振荡。鉴于N-S方程高Reynolds流具有对流占优的性质,考虑N-S方程的迎风有限元方法是自然而适当的选择。
然而,迎风方法的数值模拟一般只具有一阶精度,为了提高计算精度并且保持迎风方法计算的稳定性,一种可行的途径是研究迎风方法的自适应算法,藉以自动识别解函数的局部奇异区域,对其采用更高精度的逼近,以达到整体求解精度的提高。
后验误差估计,即利用数值解对真实误差作出的估计,能够根据误差估计结果判断出解函数的局部奇异区域,因此后验误差估计是进行自适应计算的基础。本文在第二章和第三章分别对定常和非定常N-S方程的后验误差估计进行了研究,第四章的数值算例验证了估计器的有效性和构造的自适应算法的高效性。
本文的主要结果如下:
(1)对于定常不可压N-S方程,提出了一种基于对偶剖分的迎风有限元格式,证明了此种迎风格式解的存在唯一性。在一种新的度量方式之下,给出了真实误差的后验误差估计。通过求解每一个形函数支集上的局部边值问题,解决了抽象误差估计器的计算问题。
(2)对于非定常不可压N-S方程的全离散迎风有限元格式进行了研究,讨论了该格式的唯一可解性和稳定性估计,借助于时间半离散格式,分别对时间方向和空间方向做出了后验误差估计,估计器给出了数值解与真解之误差的整体上界和局部下界。
(3)对定常不可压N-S方程的迎风有限元之后验误差估计及自适应算法,进行了题例众多的数值实验与测试性研究。首先对具有精确解的题例作了数值模拟,考察了定常N-S方程迎风方法的后验误差估计器的实用效率,验证了定常迎风格式所具有的一阶精度性能。其次构造了定常N-S方程的自适应算法,通过在规则区域上采用保持弱锐剖分的网格加密策略和不规则区域上采用自动生成自适应网格的方法,对方腔流动和圆柱绕流模型问题分别采用迎风方法、迎风自适应方法进行了计算,对计算的结果作了比较性的研究。此外,对于圆柱绕流问题中的兴趣量如拉力系数和升力系数的计算,也给出了后验的计算方法。