@Calf-Z-DNS 你好,看到你提的这个问题,请问你解决了吗?
我也想计算整个流动区域近壁面好几层网格的yPlus,而不是仅仅贴壁面的一层网格的yPlus,暂时我看到的wallShearStress和yPlus的计算都是贴壁面一层的值;
或者像东岳老师的示例代码那样直接给定了不变的uTau,而其实2维流动中,沿着流线方向uTau是变化的,所以我觉得比较复杂的是,怎么用某一壁面网格上的wallShearStress(或uTau)确定 对应法线方向上的多层网格yPlus,这个对应关系不知道如何寻找。
希望我表达清楚了,非常感谢!
shiyu 发布的最新帖子
-
RE: yPlus在openfoam代码里面的实现
-
RE: chtMultiRegionSimpleFoam前处理未能生constant下finalAgglom文件
@小火人 帖子看起来比较老了,我想请教下 对于您这个case,周围空气的进出口边界条件如何设置呢?特别关心的是压力边界条件,totalPressure吗?因为我使用totalPressure设置后结果不准确也不稳定,谢谢
-
RE: BoundaryFoam求解近壁面流动(倾斜壁面)
@李东岳 谢谢老师回复。
那我如果想要坐标转换边界条件,比如,已知local坐标系下的值,想转换成global坐标系下,方便给定openfoam的边界值,openfoam里有可以直接利用的tools吗? -
RE: BoundaryFoam求解近壁面流动(倾斜壁面)
@李东岳
李老师,因为近壁面区域的控制方程如下:
where the coordinate y is understood to be the direction normal to the wall (local coordinates)也就是说,这个控制方程是建立于local坐标系下,所以下一步是该怎么转换成global坐标系下的求解?控制方程形如boundaryFoam,但boundaryFoam的应用,找到的都是global坐标与local坐标系一致的情况。
-
BoundaryFoam求解近壁面流动(倾斜壁面)
OpenFOAM中边界条件以及solvers求解 默认是全局坐标,即 cartesian coordinates。如果是这样,那么作为1-D solver的boundaryFoam求解倾斜壁面流动时,应该依然采用global坐标系,而不是 local坐标系(水平于wall 的x1 和 垂直于wall的y1)。
但从物理上分析,应采用local坐标系下简化的边界层方程求解这种近似1-D的流动,所以问题是:怎么转换boundaryFoam的求解坐标呢?以及给定倾斜流动的边界条件时应该按照global指定?或许我思路错了,烦请解答。 -
RE: multiRegion case中采用coupled interface边界上信息来计算turbulent viscosity
@李东岳 感谢李老师的快速回复。
不过我还是不太确定,尝试了以下语句: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边界上信息来计算turbulent viscosity
如题,请问如何在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() ),谢谢!