壁面函数.pdf

FLUENT壁面函数的选择 2011-10-09 10:22:05| 分类： 默认分类|举报|字号 订阅 壁面函数问题 1、 无论是标准 k—ε 模型、RNGk—ε 模型，还是 Realizable k—ε 模型,都是针对 充分发展的湍流才有效的，也就是说，这些模型均是高 Re 数的湍流模型。它们 只能用于求解处于湍流核心区的流动。而壁面函数是对近壁区的半经验描述，是 对某些湍流模型的补充（近壁区对整体流动影响较大和低雷诺数 Re 的情况）， 通过壁面函数法和低 Re 数 k—ε 模型与标准 k—ε 模型和 RNGk—ε 模型配合，成 功解决整个整个管道的流动计算问题。 2、 在壁面区，流动情况变化很大。 解决这个问题目前有两个途径： 一、是不对粘性影响比较明显的区域（粘性底层和过渡层）进行求解，而是用一 组半经验的公式（即壁面函数）将壁面上的物理量与湍流核心区内的相应物理量 联系起来。这就是壁面函数法。在划分网格的时候，不需要在壁面区加密，只需 要把第一个节点布置在对数律成立的区域内，即配置在湍流充分发展区域。 如果要用到壁面函数的话，在 define---modle--viscous 面板里有 near wall treatment 一项。可以选择标准壁面函数、不平衡壁面函数等。 二、是采用低 Re 数的 k—ε 模型来求解粘性底层和过渡层，此时需要在壁面区划 分比较细密的网格，越靠近壁面，网格越细。当局部湍流的 Re 数小于 150 时， 就应该使用低 Re 数的 k—ε 模型。 总结： 相对于低 Re 数的 k—ε 模型，壁面函数法计算效率高，工程实用性强。 但当流动分离过大或近壁面流动处于高压之下时，不是很理想。在划分网格的时 候，需要在壁面的位置设置边界层网格，原因也是如此。 ====================================================== ============================================ 为什么要用壁面函数？？就是因为，k-epsilon 模型中，k 的 boundary condition 已知，在壁面上为零，而 epsilon 的 boundary condition 在壁面 上为一未知的非零量，如此如何来解两方程模型？？？所以，我们就需 要壁面函数来确定至少第一内节点上的值，当然也包括壁面上的值。实 际上就是把 epsilon 方程的 boundary condition 放到了流体内部。至于壁 面函数的应用范围，要看它是如何获得的，简单说，他们都是由于，靠 近壁面，雷诺应力在粘性底层内基本消失，所以，navier-stokes 变为可 解，而求得。所以，凡是应用壁面函数求得的节点，都应设置在粘性底 层(y+5-8)或者至少为线性底层(y+30?具体数值忘记了)，当然你放得 越低，精度越高，但是网格越小。我在 matlab 内自己写的 code,在 y+5-8 内放 10 层，fluent 应该可以更高。放在 fully developed region 是完全错 误的。 -------这短话理解得有问题 ====================================================== ============================================ 为什么要使用壁面函数呢？ 首先，在 CFD 中应用湍流模型并不一定需要使用壁面函数，在粘性支 层中可以对 N-S 方程直接求解。在粘性支层中，速度梯度很大，vorticity 不为零，所以要直接求解，就必须在粘性支层中布置较多节点，一般要 10 层以上，这就是一般的低 Re 数湍流模型。当然这样将占用较多的计 算资源。 而在边界层中，是存在解析解的，如果在粘性支层内不求解三维 N-S 方 程，而用一维数学模型代替，将大大降低计算资源的使用，这就是壁面 函数。一般高 Re 数湍流模型都使用壁面函数。第一层网格节点布置在 粘性支层之外。那么你如何判断你的边界层网格节点布置是否合适呢？ 这就要检查你的 y+，y+就是第一层网格质心到壁面的无量纲距离，与 速度、粘度、剪应力等等都有关系。对于 y+的值，各个学者推荐的范 围是不一样的，但一般在 30-60 之内肯定是没有问题的。也有推荐 10-110 甚至 200 的。y+的值合理，意味着你的第一层边界网格布置比较合理， 如果 y+不合理，就要调整你的边界层网格。 ====================================================== =========