?t??C?其中C??0.09。
k2? (4.1-12)
2 v2?f湍流模型
尽管上述的那些湍流模型对于流场的模拟结果已经相对比较准确了,但是对于低雷诺数
k??湍流模型,需要引入阻力方程来修正固体壁面涡粘系数的不合理的渐进特性。只要这
种阻力效应的目的是表示壁面的动能阻塞,Durbin(1991)提出了使用湍流应力取代湍动能来定义湍流粘度。因此,他介绍了一种新的v?f湍流模型,而这种模型并不需要任何阻力方程,其中v2表示湍流应力,f表示再分配。
目前流行的v?f湍流模型是是由Lien和Kalitzin(2001)提出的,这个模型建立在高雷诺数的k??模型的基础上,并增加了两个额外的方程,即湍流应力方程和再分配方程。
22Dk??Pk?Dt?xj????v?vt???k????k???x?j?? ??(4.2-1)??????x?j?? (4.2-2) ??D?C?1Pk?C?2????DtT?xjDv2???kf?6v2?Dtk?xj????v?vt??????????v?vt???k????v2?????x? (4.2-3) ?j?2L?f?f22?C1?1??v22?????CT?k?3??Pk5v2 (4.2-4) ?kTk涡粘性系数由式vt?C?v2T计算,而湍流时间T和特征长度L由下列两式定义:
????kv?0.6k??????,T?minmax?,6 2?????????6C?vS????? (4.2-5)
3234?????k32?kv?????,C? L?CLmaxmin?,342????????6C?vS????? (4.2-6)?方程中包括应变,这是Durbin(1995)引入的约束条件,以克服停滞点时产生的高应力,湍流时间和特征长度通过Kolgomorov尺度联系到一起。这个模型中的参数如下:C??0.22,
?k?1,???1.3,C?1?1.4?1?0.05??C??70。
????k2?v????,C1?1.4,C2?0.3,CL?0.23,??C?2?1.9,??对二阶矩模型评价如下:
(1)二阶矩模型基于雷诺应力输运方程,它包括雷诺应力发展过程,注入流线曲率、旋转系统等非局部性效应已经包含在雷诺应力的输运方程之中,因此它能较好地预测复杂湍流。
(2)在雷诺应力再分配项的模拟中,快速项?ij仍然采用局部性假设,这对于强非均匀湍流场的预测将产生明显的偏差。
(3)近壁的雷诺应力耗散具有强各向异性,各向同性的耗散格式有待改进。 (4)雷诺应力的梯度耗散模型基本上是合理的,但应当考虑非各向同性的扩散。 目前,RANS湍流模型是目前预测复杂湍流的最主要的工具,即使有朝一日计算机直接数值模拟复杂湍流成为现实,快速而又准确预测湍流特性的模型仍然是工程师们欢迎的方法。毋庸置疑,湍流统计模型需要不断发展和完善。
然而,我们又看到现有的湍流统计模型攒在致命的缺点:没有一个模式能够对所有的湍流运动给出满意的预测结果。一种常用的模式只能对某一类湍流运动给出满意的预测结果。例如,缓变的切边湍流,或近似平衡的湍流,B-L的代数涡粘模型是做够好的。对于非平衡的切边湍流,例如有分离的湍流边界层,非线性k??模型或SA一方程粘性系数模型可优先考虑。对于复杂几何边界的湍流,在计算条件许可的条件下,二阶矩模型可以考虑。
总之,湍流统计模型发展的目标是明确的,即简单易行,适应面广。