对某些简单的均匀时均流场,如果湍流脉动是均匀的、各向课件.pptVIP

在标准k-ε模型中， k和ε是两个基本未知量，与之相对应的输运方程为： （20） （21） 其中，Gk是由于平均速度梯度引起的湍动能k的产生项，Gb是由于浮力引起的湍动能k的产生项，YM代表可压湍流中脉动扩张的贡献，C1?、C2?和C3?为经验常数，?k和??分别是与湍动能k和耗散率?对应的Prandtl数，Sk和S?是用户定义的源项。 标准k- ?模型 精品文档 Gk是由于平均速度梯度引起的湍动能k的产生项，由下式计算： （22） Gb是由于浮力引起的湍动能k的产生项，对于不可压流体， Gb=0。对于可压流体，有： （23） 标准k- ?模型中的有关公式 精品文档 Prt是湍动Prandtl数，在该模型中可取Prt=0.85，gi是重力加速度在第i方向的分量，β是热膨胀系数，可由可压流体的状态方程求出，其定义为： （24） YM代表可压湍流中脉动扩张的贡献，对于不可压流体，YM=0。对于可压流体，有： （25） 其中，Mt是湍流Mach数， 标准k- ?模型中的有关公式 精品文档 在标准的k-ε模型中，根据Launder等的推荐值及后来的实验验证，模型常数 的取值为： （26） 对于可压缩流体的流动计算中与浮力相关的系数C3?，当主流方向与重力方向平行时，有C3?=1，当主流方向与重力方向垂直时，有C3?=0。 标准k- ?模型中的系数 精品文档 根据以上分析，当流体为不可压，且不考虑用户自定义源项时， （27） （28） 这时，标准k-ε模型变为： 精品文档 对自由射流有： —平面淹没射流 —圆淹没射流 —其中x为沿流动方向，b为射流宽度。lm与横向距离y或r无关 对充分发展管流有 —其中R为管半径，y为距管壁距离 例子：自由射流、充分发展管流 精品文档 Von Karman给出公式 对浮力流动，如为稳定分层（Ri0），则有 其中 Ri 称为梯度Richardson数，为浮力梯度与速度梯度之比，lm0为无浮力时混合长度。?1?7 例子：浮力流 精品文档 对不稳定分层（Ri 0），则有 其中?2=14 浮力或Ri越大则lm越小或湍流粘性越小，即浮力消弱湍流 湍流Prandtl数?T或Schmidt数?Y，由经验来确定 例子：浮力流 精品文档 在浮力流中，浮力对?T的影响可表达为 —其中?T0为无浮力时的湍流Prandtl数 —上式意味着浮力越大（Ri越大）则?T越大或?T及DT越小 —浮力的增大使湍流导热或湍流扩散减弱的程度比使湍流粘性减弱的程度更厉害 例子：浮力流 精品文档 优点：直观、简单，无须附加湍流特性的微分方程——适用于简单流动，如射流、边界层、管流、喷管流动等。另外，研究历史较长，积累了很多经验。 缺点1：在 处必然是湍流粘性?T为零，或剪力、热流、扩散流均为零——与实际不符 通道内中心线处?T按该理论为零 栅网后方均匀流场中的?T按该理论为零 —实际上，均不为零 混合长度模型的优缺点 精品文档 混合长度模型相当于湍流能量达到局部平衡，即湍流的产生等于湍流的耗散，亦即认为湍流的对流（上游影响）和扩散（断面上的混合）均为零。 缺点2：只有简单流动中才能给出lm的表达式。对复杂流动如拐弯或台阶后方有回流的流动，就很难给出lm的规律。 混合长度模型的优缺点 精品文档 湍流动能方程模型（单方程模型） 为了使Reynolds方程组封闭，对其中的关联项 等表征湍流特性的量继续写输运方程，其中的第一个就是Reynolds应力输运方程 周培源先生在四十年代提出的 推导应力输运方程的出发点是瞬态的N-S方程和时均Reynolds方程 精品文档 推导方法1 写出瞬时速度分量vi及vj的N－S方程 将vj乘以vi的N－S方程与vi 乘以vj的N－S方程相加，得到vi vj的方程 对上述方程进行Reynolds展开，取时平均，得到 的方程 将时均速度 乘以 的Reynolds方程与 乘以 的Reynolds方程相加，得到 的方程 上述两者相减，便得到 的方程 精品文档 推导方法2 将瞬时速度vi的N－S方程与时平均速度 的Reynolds方程相减，得到 的方程 用类似办法得到 方程 将v`j乘以v`i的方程与v`i 乘以v`j的方程相加，再取时平均，便得到 的方程 精品文档 N-S方程 右端浮力项按Boussinesq近似应为??gi，此处用体膨胀系数 精品文档 ?ik为粘性应力张量，由广义牛顿定律给出