【CFD理论】从 ΔB 到壁面剪切：CFD 软件如何处理粗糙壁面

导读从 ΔB 到壁面剪切应力。显式 ΔB 和 E′ 写法本质等价在文献、软件手册或源代码中粗糙度修正不一定总以 ΔB 的形式出现。有些代码会把 ΔB 合并到对数律常数 E 中得到一个修正后的 E′。这不是另一个物理模型而是同一个公式的代数变形。U1κln⁡(Ey)−ΔB U^ \frac{1}{\kappa} \ln(E y^) - \Delta BUκ1​ln(Ey)−ΔB等价写成U1κln⁡(E′y) U^ \frac{1}{\kappa} \ln(E^\prime y^)Uκ1​ln(E′y)其中E′Eexp⁡(κΔB) E^\prime \frac{E}{\exp(\kappa \Delta B)}E′exp(κΔB)E​这个变形的好处是实现简单。求解器仍然可以沿用光滑壁面壁函数的框架只是在粗糙壁面条件下把 E 替换为 E′。当 ΔB 0 时E′ E模型自然退回光滑壁面。近壁面剪切计算有限体积求解器在壁面相邻单元里通常仍按单元中心到壁面的距离计算速度梯度。问题是真实速度剖面是非线性的直接用线性梯度会低估或错误描述壁面剪切。壁面函数的做法是不去逐点解析剖面而是修改近壁等效黏性让线性梯度与修正后的黏性相乘时得到符合壁函数模型的壁面剪切应力。τwρνw(UPyP) \tau_w \rho \nu_w \left(\frac{U_P}{y_P}\right)τw​ρνw​(yP​UP​​)其中νwν_wνw​是近壁等效运动黏度。粗糙度通过 ΔB 或 E′ 改变νwν_wνw​最终改变壁面剪切应力τwτ_wτw​。这意味着粗糙度参数不是直接修改主流区速度而是通过壁面边界通量影响动量方程。压降、摩擦力和换热变化都是这一近壁剪切变化进一步传递到整体流场后的结果。按边界面逐个更新在实际网格中一个单元可能有多个边界面。例如二维角落单元可能同时贴着底壁和侧壁两个面的法向距离、粗糙度设置和局部剪切都可能不同。因此粗糙壁面壁函数通常按边界面逐个处理而不是只对单元给一个统一修正。遍历壁面边界面。计算当前边界面的 y。根据输入的 Ks 计算当前面的 Ks。由 Ks 和 Cs 计算 ΔB或进一步得到 E′。更新当前边界面的近壁等效运动黏度νwν_wνw​。回到动量方程、压力方程和其他输运方程求解并在下一次迭代重复。禁止随意调参在 Fluent、CFX、OpenFOAM 或其他 CFD 工具中粗糙壁面参数最终都要服务于工程。填入 Ks 和 Cs 后真正需要检查的不是云图是否更“像样”而是压降、摩擦因子、壁面剪切、换热系数等目标量是否与实验、经验公式或工程判断一致。先判断粗糙结构是否适合等效。如果粗糙元是可分辨的大尺度几何应优先显式建模。没有更多信息时Cs 可以先取 0.5但应说明这是均匀砂粒类粗糙度的常见起点。Ks 应尽量通过压降、摩擦因子或实验数据标定而不是简单取表面最大高度。粗糙壁函数仍依赖合理的近壁网格和 y 范围。网格策略不匹配时参数再细也可能只是数值补偿。对换热问题粗糙度不仅影响阻力也可能影响近壁湍流和热通量需要同时关注压降与换热收益。