yPlus在openfoam代码里面的实现

[xpqiu]

@vbcwl
请问你计算的时候，nut 的壁面边界条件设置的是什么？
postProcess 算出来的 yPlus 是调用 nut 的边界条件里面定义的 yPlus 函数算出来的。

[vbcwl]

@xpqiu 感谢前辈的回复，我用的laminar做qDNS，没有用到nut这些，只用到了p p_rgh U T。请问您了解过没有nut文件的情况么？

[xpqiu]

@vbcwl
哦，如果是这样的话，那只能通过你描述的
$$
y^+ = \frac{y u_\tau}{\nu}
$$
的方式来算 $y^+$ 了。postProcess 这种计算方式调用的是湍流模型的东西，适合用湍流模型的场景。

[shiyu]

@Calf-Z-DNS 你好，看到你提的这个问题，请问你解决了吗？
我也想计算整个流动区域近壁面好几层网格的yPlus，而不是仅仅贴壁面的一层网格的yPlus，暂时我看到的wallShearStress和yPlus的计算都是贴壁面一层的值；
或者像东岳老师的示例代码那样直接给定了不变的uTau，而其实2维流动中，沿着流线方向uTau是变化的，所以我觉得比较复杂的是，怎么用某一壁面网格上的wallShearStress（或uTau）确定 对应法线方向上的多层网格yPlus，这个对应关系不知道如何寻找。
希望我表达清楚了，非常感谢！