在高雷诺数工程计算中我们经常不想把壁面附近的黏性底层完整解析出来。原因很直接如果第一层网格要进入y1y^ 1y1网格数量、长宽比和收敛难度都会明显增加。壁面函数的作用就是允许第一层网格中心离壁面更远一些同时用经验关系补偿“壁面到第一层单元中心”之间那段没有被网格解析出来的速度变化。但这里有一个前提传统壁函数默认壁面是光滑的。只要壁面存在砂粒、锈蚀、涂层颗粒、铸造纹理或沉积物近壁速度剖面和壁面剪切应力都会改变。粗糙壁面修正要解决的就是这个问题。壁面函数在补什么在有限体积求解器中壁面相邻单元内部通常按线性变化处理。也就是说从壁面到第一层单元中心P数值格式看到的是一条近似直线。真实湍流近壁区不是这样。壁面附近速度从无滑移条件U 0起步随后经历黏性底层、缓冲层和对数区速度剖面明显是非线性的。图 1 近壁第一层网格与速度剖面图 1 想表达的不是“网格画得多精细”而是一个更基本的事实求解器手里的第一层单元中心速度UPU_PUP​并不能直接代表壁面附近真实速度梯度。如果直接用线性梯度计算壁面剪切应力会偏离真实湍流边界层。所以壁函数本质上是在回答一个问题已知第一层单元中心速度UPU_PUP​和距离ypy_pyp​怎样给出一个合理的壁面剪切应力光滑壁面对光滑壁面在常见壁函数区间内近壁速度通常写成无量纲形式U1κln⁡(Ey) U^ \frac{1}{\kappa}\ln(E y^)Uκ1​ln(Ey)其中UU^U是无量纲速度yy^y是无量纲壁面距离κ\kappaκ是 von Karman 常数E是光滑壁面对数律中的经验常数。该公式通常适用于第一层网格中心落在对数区的情形例如30y20030 y^ 20030y200。如果网格靠得很近进入y1y^ 1y1的解析近壁模型路线壁函数就不是这套逻辑了。粗糙壁面当壁面变粗糙后对数律需要向下修正。常见写法是在光滑壁面对数律后面减去一个粗糙度修正量U1κln⁡(Ey)−ΔB U^ \frac{1}{\kappa}\ln(E y^) - \Delta BUκ1​ln(Ey)−ΔB这里的ΔB\Delta BΔB可以理解为粗糙壁面对速度剖面的整体偏移量。光滑壁面时ΔB0\Delta B 0ΔB0粗糙度效应越明显ΔB\Delta BΔB越大。图 2 粗糙度修正使对数律下移图 2 只需要看一个趋势ΔB\Delta BΔB增大时同一个yy^y位置对应的UU^U会降低因此对数律曲线整体下移。这个结论容易让人误解。既然UU^U降低了是不是说明粗糙壁面附近速度反而变小、阻力也变小实际不是这样。关键在于UU^U不是实际速度UUU它是用摩擦速度uτu_\tauuτ​归一化后的量UUuτ U^ \frac{U}{u_\tau}Uuτ​U​摩擦速度和壁面剪切应力相关uττwρ u_\tau \sqrt{\frac{\tau_w}{\rho}}uτ​ρτw​​​粗糙度增大后为了让同一个外部流动满足修正后的壁面关系壁面剪切应力τw\tau_wτw​会增加摩擦速度uτu_\tauuτ​也会增加。也就是说图上的UU^U下移并不等于实际壁面阻力下降。小结壁函数用于补偿壁面到第一层单元中心之间未解析的非线性速度剖面光滑壁面对数律可以写成U1/κln⁡(Ey)U^ 1/\kappa \ln(E y^)U1/κln(Ey)粗糙壁面通过ΔB\Delta BΔB修正对数律使UU^U曲线下移但工程上对应的是壁面剪切应力和压降增加。