FLUENT壁面函数的选择

FLUENT 壁面函数的选择壁面函数的选择 壁面函数问题 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 不合理 就要调整你的边界层网格 面函数 在划分网格时 把第一个内节点布置在对数律成立的范围内 即配置到旺盛湍流区域 11 5 30 y Contours Contours of Turbulence Wall Yplus 如果 y 的值大于该范围 应该加密该区域网格 重新计算 再查看 y 如果仍不在其 范围 继续加密网格 壁面网格加密可采用自适应网格 Adapt Y Y Options 选项 只选 Refine Type 选 Y 点击 Mark 再点击 Adapt 及完成网格加密 非平衡壁面函数 Non Equilibrium Wall Function 主要应用于以下情况 涉及分离 再附着 冲击等受压力梯度影响的远离平衡的复杂流动 Enhanced Wall Treatment 要求 y 30 60 FLUENT 在 y 12 225 时候采用层流 线性 准则 因此网格不必要太密 因为壁面函数在粘性底层更本不起作用 对数区与完全湍流 的交界点随压力梯度和雷诺数变化 如果雷诺数增加 该点远离壁面 但在边界层里 必 须有几个网格点 壁面函数处理时网格划分采用双层模型时近壁网格要求当采用双层模型 时 网格衡量参数是 y 并非 y 最理想的网格划分是需要第一网格在 y 1 位置 如果稍微大点 比如 4 5 只要位于粘性底层内 都是可以接收的 理想的网格划分需 要在粘性影响的区域内 ReyY 是第一层高度一半和 viscous length scale 的比值 for Tetra mesh Y 是第一层高度 1 3 和 viscous length scale 的比值 要准确求解壁面处的流动 需要很细的网格 用壁面函数就是为了避开这一点采用的近似 处理 壁面函数在很多书和 PAPER 里都提到过 但不同模型和不同的人相差很远 而且 没有完整的步骤 我在编程中用到高雷诺数两方程模型 碰到了壁面函数的问题 1 由初 始的速度 U 按对数律计算 U 2 由 U 计算出 Y 3 判断 Y 11 5 第一内点 P 位于 旺盛湍流区 符合对数律 求 P 点 U K E 以及壁面 W 点的 U K E 4 若 Y 11 5 第 一内点 P 位于粘性支层 按 U Y 计算 这是我的理解 但更详细的细节 我还没弄清 比如 P W 点的 U V K E 的具体计算表达式 以上谈到的是规则域的壁面函数法处理 对于贴体坐标转换的壁面函数法处理起来更复杂 因为与壁面平行的速度才满足对数律 希望 CFD 朋友参加讨论 更希望提供详细的步骤和有关壁面函数法的子程序 我的理解 由于 k e 方程要求高雷诺数 所以壁面第一点应布置在粘性支层外 粘性支层 外一定范围内速度分布呈对数分布 这是流体理论的研究结果 而壁面函数主要处理的 是湍流黏度 k e 处理这些要用到这个粘性支层厚度 和速度和切应力 如果知道厚度 了 就可以根据对数分布求出速度 然后计算其他的 所以壁面函数就是要先求出粘性子 层厚度 Y 然后求 U 不要求 u 和 u u 是为了无量纲用的 用以简化推导和计算 事 实上后边用的都是 u y 知道 u 和 y 就可以干求其他 壁面函数要求第一点布置在湍流旺盛区 就是确定 y 有推导出来的表达式 而对数分 布的成立也是有范围的 所以壁面处网格的划分才是最关键的 一般通过试算搞出来 这就个大概思想 主要的理论你还得看看陶的书 peric 的书在怎么推讲得更细致一些 希 望大家补充 壁面函数 在划分网格时 把第一个内节点布置在对数律成立的范围内 即配置到旺盛湍流区域 11 5 30 y Contours Contours of Turbulence Wall Yplus 如果 y 的值大于该范围 应该加密该区域网格 重新计算 再查看 y 如果仍不在其 范围 继续加密网格 壁面网格加密可采用自适应网格 Adapt Y Y Options 选项 只选 Refine Type 选 Y 点击 Mark 再点击 Adapt 及完成网格加密 非平衡壁面函数 Non Equilibrium Wall Function 主要应用于以下情况 涉及分离 再附着 冲击等受压力梯度影响的远离平衡的复杂流动 Enhanced Wall Treatment 要求 y 4 5 1 基本思想 基本思想 对于湍流充分发展的核心流动区域使用标准的 K epsilon 模型或其改进模型求解 对壁面分 子粘性影响明显的区域 直接用半经验公式将壁面上的物理量与湍流核心区内的求解量联 系起