ARTICLE DETAIL

资讯详情

深耕网站建设与运营推广的一线实战洞察。

【CFD理论】为什么需要壁面函数

【CFD理论】为什么需要壁面函数

在高雷诺数工程计算中,我们经常不想把壁面附近的黏性底层完整解析出来。原因很直接:如果第一层网格要进入y+<1y^+ < 1y+<1,网格数量、长宽比和收敛难度都会明显增加。

壁面函数的作用,就是允许第一层网格中心离壁面更远一些,同时用经验关系补偿“壁面到第一层单元中心”之间那段没有被网格解析出来的速度变化。

但这里有一个前提:传统壁函数默认壁面是光滑的。只要壁面存在砂粒、锈蚀、涂层颗粒、铸造纹理或沉积物,近壁速度剖面和壁面剪切应力都会改变。粗糙壁面修正要解决的,就是这个问题。

壁面函数在补什么

在有限体积求解器中,壁面相邻单元内部通常按线性变化处理。也就是说,从壁面到第一层单元中心P,数值格式看到的是一条近似直线。

真实湍流近壁区不是这样。壁面附近速度从无滑移条件U = 0起步,随后经历黏性底层、缓冲层和对数区,速度剖面明显是非线性的。

图 1 近壁第一层网格与速度剖面

图 1 想表达的不是“网格画得多精细”,而是一个更基本的事实:求解器手里的第一层单元中心速度UPU_PUP,并不能直接代表壁面附近真实速度梯度。

如果直接用线性梯度计算壁面剪切,应力会偏离真实湍流边界层。

所以壁函数本质上是在回答一个问题:已知第一层单元中心速度UPU_PUP和距离ypy_pyp,怎样给出一个合理的壁面剪切应力?

光滑壁面

对光滑壁面,在常见壁函数区间内,近壁速度通常写成无量纲形式:
U+=1κln⁡(Ey+) U^+ = \frac{1}{\kappa}\ln(E y^+)U+=κ1ln(Ey+)
其中:

  • U+U^+U+是无量纲速度;
  • y+y^+y+是无量纲壁面距离;
  • κ\kappaκ是 von Karman 常数;
  • E是光滑壁面对数律中的经验常数。
    该公式通常适用于第一层网格中心落在对数区的情形,例如30<y+<20030 < y^+ < 20030<y+<200。如果网格靠得很近,进入y+<1y^+ < 1y+<1的解析近壁模型路线,壁函数就不是这套逻辑了。

粗糙壁面

当壁面变粗糙后,对数律需要向下修正。常见写法是在光滑壁面对数律后面减去一个粗糙度修正量:
U+=1κln⁡(Ey+)−ΔB U^+ = \frac{1}{\kappa}\ln(E y^+) - \Delta BU+=κ1ln(Ey+)ΔB
这里的ΔB\Delta BΔB可以理解为粗糙壁面对速度剖面的整体偏移量。光滑壁面时,ΔB=0\Delta B = 0ΔB=0;粗糙度效应越明显,ΔB\Delta BΔB越大。


图 2 粗糙度修正使对数律下移

图 2 只需要看一个趋势:ΔB\Delta BΔB增大时,同一个y+y^+y+位置对应的U+U^+U+会降低,因此对数律曲线整体下移。

这个结论容易让人误解。既然U+U^+U+降低了,是不是说明粗糙壁面附近速度反而变小、阻力也变小?实际不是这样。

关键在于,U+U^+U+不是实际速度UUU,它是用摩擦速度uτu_\tauuτ归一化后的量:
U+=Uuτ U^+ = \frac{U}{u_\tau}U+=uτU
摩擦速度和壁面剪切应力相关:
uτ=τwρ u_\tau = \sqrt{\frac{\tau_w}{\rho}}uτ=ρτw
粗糙度增大后,为了让同一个外部流动满足修正后的壁面关系,壁面剪切应力τw\tau_wτw会增加,摩擦速度uτu_\tauuτ也会增加。也就是说,图上的U+U^+U+下移,并不等于实际壁面阻力下降。

小结

  • 壁函数用于补偿壁面到第一层单元中心之间未解析的非线性速度剖面;
  • 光滑壁面对数律可以写成U+=1/κln⁡(Ey+)U^+ = 1/\kappa \ln(E y^+)U+=1/κln(Ey+)
  • 粗糙壁面通过ΔB\Delta BΔB修正对数律,使U+U^+U+曲线下移,但工程上对应的是壁面剪切应力和压降增加。
返回列表