@小火人 帖子看起来比较老了，我想请教下 对于您这个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() )，谢谢！

再次感谢回复。
这样确实可以取得pressure field并进一步可求grad(p)，对于simpleFoam or pimpleFoam是没问题的。
但，对于boundaryFoam，pressure不是通过pEqn求解，i.e. 没有在mesh上定义pressure field，仅如下这样定义还可以获取gradP的值吗？

dimensionedVector gradP
(
    "gradP",
    dimensionSet(0, 1, -2, 0, 0),
    Zero
);