FLUENT中文学习资料-建模与数值分析.docVIP

ε= k 3/ 2 lε ? Re y Au Re y Aε 4－13 上面的长度尺度根据参考文献[L29]的方法计算： l u = cl y[1 ? e lε = cl y[1 ? e ] ] 4－14 4－15 ? 如果所有的计算区域都在粘性影响的区域以内（ Re y 200 ） ，耗散率的输运方程并不需 要求解，而是用上面的代数方程来就得。上面长度尺度计算过程中的模型常数采用 Chen and Pater [L29]的结果。 ?3 / 4 cl = kC μ ， Au = 70 ， Aε = 2cl 表 4－1，几种壁面处理方法比较 优点 标准壁面函数 应用比较多，计算量较小， 有较好精度 非平衡壁面函数 双层区模型 缺点 适合高雷诺数流动，对低雷诺数 流动问题，有压力梯度，强体积力 及强三维性问题不适合 考虑了压力梯度，可以计算分离， 对低雷诺数流动问题， 重附及撞击问题 有较强压力梯度，强体积力 及强三维性问题不适合 不依赖壁面法则，对于复杂流动， 要求网格密，因而要求计算机处理 特别是低雷诺数流动很适合 时间长，内存大。 第二节，湍流计算中近壁处理对网格的要求 一个成功的湍流计算离不开好的网格。在许多的湍流中，空间的有效粘性系数不同，是平 均动量和其它标量输运的主要决定因素。因此，如果需要有足够的精度，这就需要保证湍流量 要比较精确求解。由于湍流与平均流动有较强的相互作用，因此求解湍流问题比求解层流时候 更依赖网格。 你可以用后处理面板去画出 y ， y 和 Re y 的值来检查网格是否满足自己的计算要求。需 要指出的是计算出来的 y + ， y * 和 Re y 并不是只与几何参数有关的固定量，它们也和最后的收 敛解解有关系。所以，如果你把网格加密一倍（到壁面的距离减少一半） ，计算得到的 y + 并不 一定是加密前计算的 y + 的一半。 对于近壁网格而言，不同的近壁处理对网格要求也不同。下面对常见的几种近壁处理的网 格要求做个说明。 + * 采用壁面函数时候的近壁网格： 49 PDF 文件以 FinePrint pdfFactory Pro 试用版创建 第一网格到壁面距离要在对数区内。通常计算的距离为 y + （ ≡ ρuτ y / μ 或 y * 。如果网格 在对数区内， y + 和 y * 的值差不多大小。 我们知道， 对数区的 y + 30～60。 FLUENT 在 y + 12.225 时候采用层流（线性）准则，因此网格不必要太密，因为壁面函数在粘性底层更本不起作用。 对数区与完全湍流的交界点随压力梯度和雷诺数变化。如果雷诺数增加，该点远离壁面。但在 边界层里，必须有几个网格点。 壁面函数处理时网格划分 采用双层模型时近壁网格要求 当采用双层模型时，网格衡量参数是 y + ，并非 y * 。最理想的网格划分是需要第一网格在 y + ＝1 位置。如果稍微大点，比如 y + ＝4～5，只要位于粘性底层内，都是可以接收的。理想 的网格划分需要在粘性影响的区域内（ Re y 200 ）至少有十个网格，以便可以计算粘性区域 内的平均速度和湍流量。 采用双层区模型时网格划分 采用 Spalart-Allmaras 模型时的近壁网格要求 该模型属于低雷诺数模型。这就要求网格能满足求解粘性影响区域内的流动，引入了阻尼 函数，用以削弱粘性底层的湍流粘性影响。因此，理想的近壁网格要求和采用双层模型时候的 网格要求一致。 50 PDF 文件以 FinePrint pdfFactory Pro 试用版创建 采用大涡模拟的近壁网格要求 对于大涡模拟，壁面条件采用了壁面法则，因此对近壁网格划分没有太多限制。但是，如 果要得到比较好的结果，最好网格要细，最近网格距离壁面在 y + ＝1 的量级上。 第三节，用 FLUENT 求解湍流问题设定 1， 击活粘性模型面板上的湍流模型（Spalart-Allmaras, k-epsilon, Reynolds Stress or Large Eddy Simulation ）,如果选择 k-epsilon 模型， 将需要继续选择采用标准 k ? ε 模型、重整化群 k ? ε 模型或可实现 k ? ε 模型中的一种。 如果流动问题中包含壁面，选择 k ? ε 或者雷诺应力模型，在粘性模型面板上选择 近壁处理方式。 近壁处理方式包括： 标准壁面函数； 非平衡壁面函数和双层区模型。 Option 选项设置 变量的边界条件设置 全场变量赋初始值 2， 3， 4， 5， 任意选项设置 如果选择 Spalart-Allmaras 模型，会出现如下选项： . Vorticity-based production . Strain/vorticity-based production .