摩擦速度比论文偏小

[Samuel-Tu]

各位CFDer，我最近做了一个高低屋面的流场计算，流线如下图：

计算稳定后（低屋面上的分离点位置不再变化），提取低屋面上的摩擦速度，发现摩擦速度比论文里说的偏小，但分离点的位置（摩擦速度为0的点）对应上了。
论文摩擦速度如下（看u*=0.8的线）：

我算出来的摩擦速度如下：

想不通是什么原因，我的yPlus如下：

使用的realizedKE（和论文一致，但论文使用的是FLUENT，未说明壁面模型），我的k壁面使用的是kqRWallFunction（高低雷诺数都适用），我的epsilon壁面使用的是epsilonWallFunction（高低雷诺数都适用），我的nut壁面使用的是nutkWallFunction（高低雷诺数都适用）
请大家帮忙分析下原因

[Samuel-Tu]

另外，我的摩擦速度是用壁面剪切应力wallShearStress的X分量开平方计算出来的，我的模型是不可压流体，因此没有除以密度，这个计算应该没问题吧。

[我是河滩]

@Samuel-Tu 在 摩擦速度比论文偏小 中说：

kqRWallFunction是零梯度，应该只适合高雷诺数，kLowReWallFunction才是高低雷诺数都适合。epsilonWallFunction应该也只适合高雷诺数，epsilonLowReWallFunction才是高低雷诺数都适用。

[Samuel-Tu]

@我是河滩 kqr确实是我记错了，不过epsilon我印象中已经合并了，只有一个epsilonwallfunction了

[hhh]

请问您最后是怎样计算磨擦速度的呢，我的算例计算磨擦速度也偏小一点。