湍流模型及在CFD中的应用.ppt

湍流模型及其在CFD中的应用 一、湍流及其数学描述 1、湍流流动的特征 流体实验表明，当Reynolds数小于某一临界值时，流动是平滑的，相邻的流体层彼此有序地流动，这种流动称为层流（laminar flow）。当 Reynolds数大于临界值时，会出现一系列复杂的变化，最终导致流动特征的本质变化，流动呈无序的混乱状态。这时，即使是边界条件保持不变，流动也是不稳定的，速度、压力、温度等流动特性都随机变化，这种状态称为湍流（turbulent flow）. 湍流流动的两个例子 湍流涡的特点 观测表明，湍流带有旋转流动结构，这就是湍流涡（turbulent eddies），简称涡（eddy）。 从物理结构上看，可以把湍流看成是由各种不同尺寸的涡叠合而成的流动，这些涡的大小和旋转轴的方向分布是随机的。 大尺度的涡主要是由流动的边界条件所决定，其尺寸可以与流场的大小相比拟，它主要受惯性影响而存在，是引起低频脉动的原因； 小尺度的涡主要是由粘性力所决定的，其尺寸可能只是流场尺度的千分之一量级，是引起高频脉动的原因。 涡的生成与耗散 大尺寸的涡不断地从主流中获得能量，通过涡间相互作用，能量逐渐向小尺寸的涡传递。 最后由于流体粘性的作用，小尺度的涡就不断消失，机械能就耗散为流体的热能。 同时由于边界的作用，扰动及速度梯度的作用，新的涡又不断产生，构成了湍流运动。 湍流的作用 由于湍流的存在，速度脉动量在流线方向的分量和垂直于流线方向的分量之间建立了关联量，它代表着一种横向交换通量，也可以认为是由于湍流流动引起的一种附加剪切应力——影响动量的输运过程。 湍流的存在使传热和传质通量提高。 由于湍流会促进这些基本过程，因此对某些物理现象就会产生强烈的影响，如，脉动过程的消衰、均相化学反应率的增加以及液滴蒸发的强化。 2、湍流的基本方程 无论湍流运动多么复杂，非稳态的连续方程和N-S方程对于湍流的瞬时运动仍然是适用的。在此，考虑不可压流动，使用笛卡尔坐标系，速度矢量在x、y和z方向的分量分别为u、v和w，写出湍流瞬时控制方程如下： 定义时均量 为了考察脉动的影响，目前广泛采用的是时间平均法，即把湍流运动看做由两个流动叠加而成，一是时间平均流动，二是瞬时脉动流动。这样，将脉动分离出来，便于处理和进一步探讨。现在，引入Reynolds平均法，任一变量?的时间平均定义为： 时均量与脉动量的关系 湍流时均流动的控制方程 时均输运方程的统一形式 密度脉动的影响 以上是假设流体密度为常数； 但是在实际流动中，密度可能是变化的。 Bradshaw等指出，细微的密度变动并不对流动造成明显的影响 在此，忽略密度脉动的影响，但考虑平均密度的变化，写出可压湍流平均流动的控制方程如下 注意，为方便起见，除脉动值的时均值外，下式中去掉了表示时均值的上划线符号“—”，如 用φ表示 湍流输运方程组 时均形式的连续方程 张量形式的时均输运方程 二、湍流的数值模拟方法简介 1、三维湍流数值模拟方法的分类 湍流数值模拟方法可以分为直接数值模拟方法和非直接数值模拟方法。 所谓直接数值模拟方法是指求解瞬时湍流控制方程。 非直接数值模拟方法就是不直接计算湍流的脉动特性，而是设法对湍流做某种程度的近似和简化处理，例如前面提到的时均性质的 Reynolds方法就是其中的一种典型方法。 根据依赖所采用的近似和简化方法不同，非直接数值模拟方法分为大涡模拟、统计平均法和Reynolds平均法。 2、直接数值模拟（DNS）简介 直接数值模拟方法就是直接用瞬时的N-S方程对湍流进行计算，其最大的好处是无需对湍流流动做任何简化或近似，理论上能得到相对准确的计算结果。 但是，DNS要求网格划分的非常细，对计算机内存空间及计算速度要求非常高，目前还无法用于真正意义上的工程计算。 3、大涡模拟（LES）简介 由于就目前的计算能力而言，能够采用的计算网格的最小尺度仍然比最小涡的尺度要大许多。因此，目前只能放弃对全尺度范围上涡运动的模拟，而只将比网格尺度大的湍流运动通过N-S方程直接计算出来，对于小尺度涡对大尺度运动的影响则通过建立模型来模拟，从而形成了大涡模拟法（LES）。 LES方法的基本思想 用瞬时的N-S方程直接模拟湍流中的大尺度涡，不直接模拟小尺度的涡，而小涡对大涡的影响通过近似的模型来考虑。 LES方法对计算机内存及CPU速度要求仍然很高，但是低于DNS法。 4、Reynolds平均法（RANS）简介 虽然瞬时N-S方程可以描述湍流，但是N-S方程的非线性使得用解析方法精确描写三维时间相关的全部细节极端困难。从工程应用的观点来看，重要的是湍流所引起的平均流场的变化，是整体效果。 Reynolds平均法的核心不是直接求解瞬时的N-S方程，而是想办法求解时均化的Reynolds方