计算流体力学中的数值摄动算法研究-计算数学专业论文.docxVIP

摘 要 计算流体力学CCFD)有着广泛的科学和工程应用价值，各种计算 格式和算法也一直是CFD研究的热点。高智提出的数值摄动算法是计 算流体力学计算格式高精度高分辨率重构的一个通用方法，这是一种 精度高、分辨率高、稳定范围大、计算效率高、实用性强的高性能算 法。 数值摄动算法C简称高算法)通过对方程和重构格式本身的运算， 并引入下游不影响上游的对流运动规律及消除重构格式修正微分方程 的截断误差诸项，最终获得数值摄动重构格式C简称高格式)。很显然， 双由方程、Euler方程组、对流扩散方程和Navier-StokesCNS)方程组的高 格式必将形成一个大的格式家族，许多重要的家族成员如三结点高格 式己经提出，但是，许多完全可以预期的高格式，像多结点高格式远没 有构建出来。 本论文主要分两部分: 1. 本文利用数值摄动算法，把对流扩散方程的常用QUICK格式C粘 性和对流项分别用二阶中心和QUICK格式离散)进行了数值摄动高精 度重构，包括利用离散单元内所有结点的全域重构和分别利用离散单 元内上下游结点的上下游重构，即在重构中引入下游不影响上游的对 流运动规律，得到两类新的更高阶精度的数值摄动重构格式，称为高 的QUICK格式C记作G-QUICK格式)。G-QUICK格式与QUICK格式相比 简单性相当，但精度更高:全域重构G-QUICK格式和QUICK格式均为条 件稳定，上下游重构得到一些绝对稳定的G-QUICK格式。解析分析和数 值算例均证实了G-QUICK格式的优良性能，上下游重构的G-QUICK格式 为在对流扩散方程的QUICK格式中避免使用人工粘性提供了新途径。 2. 本文把一阶迎风有限体积格式的高算法重构格式 高的迎风 有限体积格式CGUVS)与不可压缩流和可压缩流SIMPLE算法相结合， 数值模拟了各种宽深比的亚跨声速腔内流动。利用六阶高迎风有限体 积格式CGUVS)在结构和非结构同位网格上模拟不同宽深比亚跨声 速腔内流动，获得较好的数值结果。例如亚声速方腔流动的数值结果 与Ghia等用多重网格技术获得的Benchmark解吻合得非常好:对不同宽 深比腔内流动，数值模拟了大尺度旋涡、腔角小尺度旋涡的形成、分 I I II 岔等演化情况，表明六阶精度GUVS具有分辨率高、稳定性好的良好性 能。 关键词: 计算流体力学，对流扩散方程，QUICK格式，数值摄动高 精度重构，SIMPLE算法 PAGE PAGE III Abstract Computational fluid dynamics(CFD) have wide application in many science and technology fields, and constructing calculation schemes and algorithm have been the research focus of CFD. The numerical perturbation algorithm presented by Gao Zhi is one of the general method for constructing high-accurate and high resolution CFD schemes, which have excellent properties of high accurate, high resolution, large stable range and high efficiency. The numerical perturbation reconstruction schemes (call them Gao’s schemes for brevity) are obtained by operating the equations and reconstruction schemes, introducing the convection law of downstream doing not affect the upstream and eliminating the truncation error of terms in the modified differential equation of reconstruction scheme. Obviously, Gao’s schemes of hyperbolic equation, Eu- ler equations, convective-diffusion equation and Navier-Stokes(NS) equations will form a scheme family with a large