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

@小火人 帖子看起来比较老了，我想请教下 对于您这个case，周围空气的进出口边界条件如何设置呢？特别关心的是压力边界条件，totalPressure吗？因为我使用totalPressure设置后结果不准确也不稳定，谢谢

@疏影横斜水清浅
谢谢回复，是的，预计和你问题一样的，我正在重新检查代码，编译器不报错的话，只能一行行自己排查吗？

@疏影横斜水清浅
请问最后找到具体问题了吗？是代码还存在问题吗？我也遇到了一样的问题，wmake libso成功，但无法load，谢谢

@李东岳 谢谢老师回复。
那我如果想要坐标转换边界条件，比如，已知local坐标系下的值，想转换成global坐标系下，方便给定openfoam的边界值，openfoam里有可以直接利用的tools吗？

@李东岳
李老师，因为近壁面区域的控制方程如下：

where the coordinate y is understood to be the direction normal to the wall (local coordinates)

也就是说，这个控制方程是建立于local坐标系下，所以下一步是该怎么转换成global坐标系下的求解？控制方程形如boundaryFoam，但boundaryFoam的应用，找到的都是global坐标与local坐标系一致的情况。

@李东岳 感谢回复
不考虑重力的，只考虑坐标变换

OpenFOAM中边界条件以及solvers求解 默认是全局坐标，即 cartesian coordinates。如果是这样，那么作为1-D solver的boundaryFoam求解倾斜壁面流动时，应该依然采用global坐标系，而不是 local坐标系（水平于wall 的x1 和 垂直于wall的y1）。
但从物理上分析，应采用local坐标系下简化的边界层方程求解这种近似1-D的流动，所以问题是：怎么转换boundaryFoam的求解坐标呢？以及给定倾斜流动的边界条件时应该按照global指定？或许我思路错了，烦请解答。

@李东岳 感谢李老师的快速回复。
不过我还是不太确定，尝试了以下语句：

label patchID = this->mesh_.boundaryMesh().findPatchID("A_to_B");  
fvPatch& patch = this->mesh_.boundaryMesh()[patchID];

编译报错如下：

error: type 'Foam::fvPatch' does not provide a call operator
          refCast<const mappedPatchBase>(patch().patch());

该如何解决呢？网上查了半天没有类似的问题。

如题，请问如何在multiRegion case中，获取coupled interface边界上信息，比如速度的梯度snGrad，来计算湍流模型中的turbulent viscosity？
已知在一些coupled的边界条件中，如turbulentTemperatureCoupledBaffleMixedFvPatchScalarField.C，
我们可以采用

// Get the coupling information from the mappedPatchBase
const mappedPatchBase& mpp =
    refCast<const mappedPatchBase>(patch().patch());
const polyMesh& nbrMesh = mpp.sampleMesh();

但问题是 在湍流模型中计算turbulent viscosity的correctNut( )，也需要调用nbrMesh来获取隔壁信息，该怎么定义patch( )呢？(边界条件中不需要额外定义patch() )，谢谢！