【CFD理论】壁面函数粗糙度参数定义与应用
导读：讨论粗糙度参数到底怎么进入模型？
粗糙管实验
粗糙壁函数的许多经验关系，背后都可以追溯到经典粗糙管实验。
实验思路并不复杂：在管壁上布置不同尺寸的均匀砂粒，然后测量湍流管流中的速度剖面和摩擦阻力。随着砂粒尺寸增大，近壁速度剖面变得更“宽”，为了推动同样的流量，需要更大的压力梯度。
这类实验给 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 则把这些粗糙度信息转化为对数律修正。
