【CFD理论】从 ΔB 到壁面剪切

导读：从 ΔB 到壁面剪切应力，浅谈CFD 软件如何处理粗糙壁面。

在文献、软件手册或源代码中，粗糙度修正不一定总以 ΔB 的形式出现。有些代码会把 ΔB 合并到对数律常数 E 中，得到一个修正后的 E′。这不是另一个物理模型，而是同一个公式的代数变形。

等价写成：

其中：

这个变形的好处是实现简单。求解器仍然可以沿用光滑壁面壁函数的框架，只是在粗糙壁面条件下把 E 替换为 E′。当 ΔB = 0 时，E′ = E，模型自然退回光滑壁面。

有限体积求解器在壁面相邻单元里，通常仍按单元中心到壁面的距离计算速度梯度。问题是，真实速度剖面是非线性的，直接用线性梯度会低估或错误描述壁面剪切。

壁面函数的做法，是不去逐点解析剖面，而是修改近壁等效黏性，让线性梯度与修正后的黏性相乘时，得到符合壁函数模型的壁面剪切应力。

其中      是近壁等效运动黏度。粗糙度通过 ΔB 或 E′ 改变     ，最终改变壁面剪切应力     。

这意味着粗糙度参数不是直接修改主流区速度，而是通过壁面边界通量影响动量方程。压降、摩擦力和换热变化，都是这一近壁剪切变化进一步传递到整体流场后的结果。

在实际网格中，一个单元可能有多个边界面。例如二维角落单元可能同时贴着底壁和侧壁，两个面的法向距离、粗糙度设置和局部剪切都可能不同。因此，粗糙壁面壁函数通常按边界面逐个处理，而不是只对单元给一个统一修正。

1. 遍历壁面边界面。
2. 计算当前边界面的 y+。
3. 根据输入的 Ks 计算当前面的 Ks+。
4. 由 Ks+ 和 Cs 计算 ΔB，或进一步得到 E′。
5. 更新当前边界面的近壁等效运动黏度       。
6. 回到动量方程、压力方程和其他输运方程求解，并在下一次迭代重复。

在 Fluent、CFX、OpenFOAM 或其他 CFD 工具中，粗糙壁面参数最终都要服务于工程。填入 Ks 和 Cs 后，真正需要检查的不是云图是否更“像样”，而是压降、摩擦因子、壁面剪切、换热系数等目标量是否与实验、经验公式或工程判断一致。

• 先判断粗糙结构是否适合等效。如果粗糙元是可分辨的大尺度几何，应优先显式建模。
• 没有更多信息时，Cs 可以先取 0.5，但应说明这是均匀砂粒类粗糙度的常见起点。
• Ks 应尽量通过压降、摩擦因子或实验数据标定，而不是简单取表面最大高度。
• 粗糙壁函数仍依赖合理的近壁网格和 y+ 范围。网格策略不匹配时，参数再细也可能只是数值补偿。
• 对换热问题，粗糙度不仅影响阻力，也可能影响近壁湍流和热通量，需要同时关注压降与换热收益。