导读讨论粗糙度参数到底怎么进入模型粗糙管实验粗糙壁函数的许多经验关系背后都可以追溯到经典粗糙管实验。实验思路并不复杂在管壁上布置不同尺寸的均匀砂粒然后测量湍流管流中的速度剖面和摩擦阻力。随着砂粒尺寸增大近壁速度剖面变得更“宽”为了推动同样的流量需要更大的压力梯度。这类实验给 CFD 建模提供了一个重要启发不一定需要把每一个砂粒都画进网格只要粗糙元的平均效应可以通过壁面剪切应力体现出来就可以用等效粗糙度参数进入壁函数。粗糙度高度 Ks粗糙度高度 Ks 可以理解为粗糙元向流体中伸出的代表性高度。如果壁面上是尺寸均匀、紧密排列的砂粒那么 Ks 可以近似对应砂粒直径如果表面粗糙度来自加工纹理、腐蚀坑、涂层颗粒或沉积物Ks 就更像一个等效参数。这也是很多工程计算容易出问题的地方。真实表面可能有峰值高度、平均高度、均方根高度等多种统计描述壁面函数需要的是能够反映壁面阻力效应的等效高度。Ks最好通过压降、摩擦因子或已有实验数据来标定而不是简单地取表面最高点。的形式与 y 类似只是把壁面法向距离 y 换成了粗糙度高度 。粗糙度常数 Cs仅有 Ks 还不够因为同样高度的粗糙元形状和排布不同对流动的影响也不同。均匀砂粒、规则沟槽、随机腐蚀点、涂层颗粒给近壁动量交换带来的效果并不一样。粗糙度常数 Cs 的作用就是用一个经验系数补偿粗糙元形状和分布差异。在很多 CFD 软件手册中均匀、紧密排列的砂粒粗糙度常以 Cs 0.5 作为推荐值。对于非均匀、非规则粗糙表面Cs 可能需要更高或需要校准。但从工程稳健性看通常不建议同时随意调整 Ks 和 Cs更常见的做法是先固定 Cs再通过 Ks 标定目标压降或壁面剪切。ΔB 的三段区域ΔB 是粗糙度对对数律的修正量但它并不是随 Ks 从零开始线性增加。根据粗糙元相对黏性底层和近壁尺度的大小通常可以把粗糙壁面分为三个区域。中间过渡区的函数形式往往来自实验数据拟合看起来比较复杂。Ks 很小时粗糙元对主流几乎不可见Ks 进入过渡区后粗糙度开始增强壁面剪切Ks 很大时壁函数本身逐渐接近适用极限。粗糙壁面函数的适用边界粗糙壁函数是一种等效模型它假设粗糙元可以被平均化为壁面边界条件的一部分。如果粗糙元尺寸已经接近近壁第一层网格高度或者已经是明显的离散几何结构例如凸台、沟槽、大颗粒、障碍物那么继续用 Ks 等效就会偏离实际。这时更合理的做法是把这些几何特征显式建模并通过网格解析其绕流和分离而不是把所有影响都压缩进一个粗糙度高度。Ks 决定粗糙元的代表高度Ks 决定它在近壁摩擦尺度下是否“有效”Cs 补偿粗糙元的形状和排布ΔB 则把这些粗糙度信息转化为对数律修正。