论文部分内容阅读
大涡模拟(LES)能实现对非定常湍流流动的精细模拟,优于雷诺平均(RANS),但是计算量较大。随着计算机计算水平的提高,将LES用于实际工程是一种发展趋势。LES的思想是采用低通滤波器对Navier-Stokes(NS)方程进行过滤,仅对大尺度的涡进行直接模拟,丢失的小尺度信息则采用亚格子模型来描述。根据亚格子尺度的不同特性,已经发展出多种亚格子模型,但是目前的亚格子模型均没有说明,在何种情况下,LES能够直接过渡成直接数值模拟(DNS)进行计算。为了克服这个缺点,鉴于旋流强度能够较好的反映旋涡特性,本文结合涡粘性假设,提出了一种基于旋流强度的亚格子模型。该模型具有当地无旋涡时,湍流粘性系数自动为零的特性。并将该模型用于低Re数和高Re数下不可压的湍流方柱绕流和不可压低湍流度的方腔流的模拟。此外,本文采用投影法离散控制方程,并对空间对流项进行改进了,使其达到了4阶精度,并用于LES计算。 采用基于旋流强度的亚格子模型,对低Re数和高Re数下不可压的湍流方柱绕流模拟时,结果表明: (1)在低Re数(Re∈[2,5×103,10×103])和高Re数(Re∈[1.25×105,3.5×105])下,St数均很明显与Re数不相关,平均阻力系数CD也基本与Re数不相关,平均升力系数CL的平均值为0,与Re数无关,与实验结果一致。但是,当Re∈[1,0×103,2.5×103]时,平均阻力系数CD随着Re的增大而增大。另外,在方柱绕流问题中,没有出现类似圆柱绕流中的阻力危机临界现象。 (2) Re∈[1.25×105,3.5×105]时,归一化的平均速度及速度脉动分布受Re数影响较小,计算结果与实验相符。 (3)vST的平均值与均方平均值均随Re数的增大而增大,它的峰值一般出现在方柱下游涡相互作用比较强的2个区域内。 (4) Komogorov微尺度分布显示,近壁面的湍动能耗散较强。湍流间歇因子(FSI)的等值面显示在旋涡脱落的初始位置分成薄片状的旋涡结构,而在尾迹区则表现为很多相对较小的旋涡结构。 (5) t-z平均流场不沿着y=0.5对称,导致顶部的面Nusselt数t-z平均值和底部的不同。 在方腔自然对流问题中,结果表明: (1)通过壁面附近的温度和速度分布及不同截面上的涡量分布可知,对方腔自然对流问题,在方腔的壁面附近会形成较大的气体环流带,伴随着壁面附近的环流,在壁面环流与方腔的核心区之间形成了一些与壁面环流方向相反的旋涡,在方腔的左上角与右下角将出现较多、较复杂的细小旋涡结构。 (2) Ra=1.58×109时,中平面上(z=0)水平中心线的速度分布与实验值在壁面最里层区域和方腔的核心区内符合的较好,但水平中心线左侧热壁附近的平均速度最大值的计算值大于实验值。垂直中心线的速度与实验值符合得较好。中平面水平线上无量纲温度的计算值略低于实验值,垂直中心线上的分布符合得较好。但是它们都会受到Ra数变化的影响。 (3) Ra=1.58×109时,u和(☉)与实验值符合得较好,而u的计算值则明显被低估了,因而最后也导致热壁附近的湍流剪切应力被低估。随着当Ra数的增大,仅壁面环流的主流速度的脉动获得了加强,而垂直于壁面的速度的脉动变化不明显,温度的脉动也被抑制。 (4) Ra=1.58×109时,各面的Nusselt数与实验值符合得较好,与理论值相差小于8%。顶面的Nusselt数沿x方向逐渐增大,底面的Nusselt数沿x方向逐渐减小。左侧热壁的Nusselt数沿y方向逐渐减小,右侧冷壁的Nusselt数沿y方向逐渐增大。 通过对低Re数和高Re数下不可压的湍流方柱绕流及不可压方腔自然对流问题计算结果与实验值和理论值的比较表明,本文提出的基于旋流强度的亚格子模型及改进的有限差分算法在此类问题的计算中是有效的。