高等工程流体力学(新).ppt

在一方程模型中的雷诺应力或 项仍用涡粘性模型，不过其中的涡粘性系数要用κ与l来表示 于是 (3-39) (3-40) 其中常数 。 在一方程模型范围内，同样也有很多模型。 3.7 双尺度二阶湍流模型 尽管上述的二阶封闭模型比之于一阶封闭模型大大改进了计算的精度与进行预报的能力，然而当用同一个模型于不同的流动时，时常需要调整经验常数的数值。这是不能令人满意的，我们仍需要继续探讨二阶模型的改进。 最不能令人满意的是对ε方程的模拟。其中对方程右边每一个项的模拟都使人不能放心。陈景仁（1986）提出应考虑对不同的项采取不同的湍流尺度。某些项，例如扩散项，主要贡献来自含能涡，所需的尺度应通过κ和ε来表示，如 而另一些项，如小涡拉伸产生项 与粘性破坏项 则主要由耗散范围内的小尺度涡决定，应采用柯尔莫戈洛夫定义的尺度，有关尺度可用ε和ν来表示，如 仅仅根据这一想法，其他都不变，则修改后的方程为 (3-41) 其中的常数仍要通过实验来确定。得到的值是 当用双尺度的 模型来计算各种自由剪切湍流时，发现比单尺度的 模型有令人满意的改进，现在用同一组常数可以计算二维的与轴对称的射流，以及平面尾流与平面混合层，都得到与实验能很好符合的结果。 3.8 湍流模型评价 应该期望雷诺应力模型比其他较简单的模型有更广泛的适用范围，更高的预报能力和更好的计算精度。但由于模型的复杂性，计算工作量远高于其他模型，特别对复杂的三维湍流流动，至今已有的算例仍不算多，目前还不能做出明确肯定的结论，尤其是从获得的收益与付出的代价相比来看。 对于大多数适用薄剪切层近似的简单剪切流，在雷诺应力各分量中只有一个剪切应力分量对平均流计算有影响。在这种情况下，涡粘性模型与代数应力模型实际上给出相同的剪应力表达式，尽管对正应力分量的表达式有很大区别。因此，就平均流计算而言，两者没有区别。如能将适用于壁面邻近区域的低雷诺数修正考虑进去，则 涡粘性模型可给出相当满意的结果。而在数值计算上 涡粘性模型显然要比代数应力模型简单得多。这就是为什么 涡粘性模型得到最广泛应用的原因。 但是对以下这些较复杂的流动，如用 涡粘性模型，则或者精度下降或者完全失败。例如，在激波后的边界层内，在边界层分离点附近以及有分离区的流动，流动沿流向有显著变化，正应力的梯度要起重要作用，涡粘性模型不能恰当地描写流动的发展；在不对称槽流。弯曲的混合层或壁射流等情形，剪应力等于零的点与 等于零的点不重合，涡粘性模型不可能反映此流动特点；在非圆截面管道或明渠流中，在横截面内由湍流应力场的强各向异性引起的二次流是很重要的，而各向同性的涡粘性模型则完全不能产生二次流。 在这些情形中，比 模型更简单的模型当然都不适用，只有代数应力模型和雷诺应力模型取得了成功。而现有的算例表明雷诺应力模型并未显示出比代数应力模形有显著的改进。 主要变量控制方程的通用形式 质量守恒方程 [单位时间内微元体中流体质量的增加] =[同一时间间隔内流入该微元体的净质量] 不可压条件 动量守恒方程 [微元体内流体动量的增加率]=[作用在微元体上各种力之和] 惯性力=质量力+表面力 牛顿广义切应力公式 X- 方向动量 y- 方向动量 z- 方向动量 X- 方向动量 y- 方向动量 z- 方向动量 再整理，得： 能量守恒方程 通用形式 瞬态 对流 扩散 源项 下面就依次讨论需要模拟的各项的湍流模型。 1 湍流扩散项的模拟