对某些简单均匀时均流场,如果湍流脉动是均匀、各向.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-ε模型变为： 采用标准k-ε模型求解流动及换热问题时，控制方程包括连续性方程、动量方程、能量方程、k方程、?方程与式（20）。若不考虑热交换的单纯流场计算问题，则不需要包含能量方程。如考虑传质或者有化学变化的情况，则应再加入组分方程。这些方程都可以表示为下面通用形式： （29） 若用散度符号，上式记为： （30） 标准k- ?模型的控制方程组 方程 φ 扩散系数Γ 源项S 连续 1 0 0 x-动量 u y-动量 v z-动量 w 湍动能 k 耗散率 ε 能量 T S按实际问题而定 标准k- ?模型的控制方程组 在将各变量的控制方程都写成式（30）所示的统一形式后，控制方程的离散化及求解方法可以得到统一，不同变量的区别仅在于扩散系数、广义源项及初值、边界条件这三方面。 标准k- ?模型的解法 1）模型中的有关系数，主要根据一些特殊条件下的试验结果而确定的，在不同的文献讨论不同的问题时，这些值可能有出入。在数值计算的过程中，针对特定的问题，参考相关文献，寻求更合理的取值。 标准k- ?模型的适用性 2）上述k- ?模型，是针对湍流发展非常充分的湍流流动来建立的，是一种针对高Re数的湍流计算模型，而当Re数较低时，例如，在近壁区内的流动，湍流发展并不充分，湍流的脉动影响可能不如分子粘性的影响大，在更贴近壁面的底层内，流动可能处于层流状态。因此，对Re数较低的流动使用上面建立的k-?模型进行计算，就会出现问题。这时，必须采用特殊的处理方式，以解决近壁区内的流动计算及低Re数时的流动问题。使用上面的k-?模型可能就会出现问题。常用解决方法有壁面函数法和低Re数的k-?模型。 标准k- ?模型的适用性 3） 标准k-?模型比零方程模型和一方程模型有了很大改进，但是对于强旋流、弯曲壁面流动或弯曲流线流动时，会产生一定失真。原因是在标准k-?模型中，对于Reynolds应力的各个分量，假定粘度系数?t是相同的，即假定μt是各向同性的标量。而在弯曲流线的情况下，湍流是各向异性的，?t应该是各向异性的张量。 标准k-?模型的适用性 需要注意的是，RNG k-?模型仍针对充分发展的湍流是有效的，是高Re数的湍流计算模型，而对近壁区内的流动及Re数较低的流动，必须使用下面将要介绍的壁面函数法或低Re数的k-?模型来模拟。 此外，需要说明一点，在FLUENT手册中，将RNG k-?模型所引入的反映主流的时均应变率Eij的一项，归入了?方程中的C2?系数中，且表达式多了一个系数Cμ，而不像这里归入C1?系数。这两种处理方式实质上是一样的。 RNG k-?模型 标准k-ε模型对时均应变率特别大的情况，有可能导致负的正压力。为使流动符合湍流的物理定律，需要对正应力进行某种数学约束。为保证这种约束的实现，认为湍流粘度计算式中的系数Cμ不应是常数，而应与应变率联系起来。从而，提出了Realizable k-ε模型。这里， Realizable有“可实现”的意思。 Realizble k-?模型 （34） （35） 其中， （36） Realizble k-?模型形式 式中，μ与 Cμ按下式计算： （