CFD中文网

shiyu

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

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());

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

shiyu

如题，请问如何在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() )，谢谢！

shiyu

再次感谢回复。
这样确实可以取得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
);

shiyu

谢谢李老师的回复 @李东岳。
边界类型没问题，但我不知道-- 怎么通过语句获取这个已经定义在solver里的gradP？
比如：想要在边界条件中，获得速度U，可以通过以下语句来获得：

volVectorField& U =
         const_cast<volVectorField&>(mesh.lookupObject<volVectorField>("U"));

但这是因为U定义在mesh上，对于gradP这种变量的获取，我有些迷茫。

shiyu

问题：如何在自定义边界条件中获取solver里定义的某一个参数？
举例：boundaryFoam中定义了pressure gradient ： gradP如下，

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

因为这个gradP每步更新，如下

gradP += (Ubar - UbarStar)/(1.0/UEqn.A())().weightedAverage(mesh.V());

我需要在自定义边界条件中获取gradP进行计算，不知具体如何操作？

谢谢大家！

shiyu

@cfd_上弦月谢谢回复，有道理，但因为是非均匀网格，即越靠近壁面越密集，我其实想要获取上图中绿色区域内每个网格的高度，（假设长度都是一致的），用于高度方向的积分运算。

shiyu

@xpqiu 您好，请问可以发我一份完整版吗？邮箱 lyushiyu@163.com，谢谢

shiyu

感谢东岳老师的付出，貌似金币还差一点，可能的话愿意有偿补充。
想要一本 Turbulence modelling for CFD，
邮箱 lyushiyu@163.com，
再次感谢！

shiyu

谢谢回复，我进一步解释一下，拿你的图来举例，当我设置橘黄色interface边界（介于zone 1 和zone 2）时，我想获取绿色方框内的网格中心处的物理信息，比如速度，他们的特点是这一系列网格垂直于我想要设置interface面，不知该怎么写程序。

shiyu

你好，我暂时不清楚你的问题，不好意思，不过想提一个类似的问题，我的case也有两个cellZone，想问在设置边界的时候可以提取出某一个cellZone里normal方向的所有网格信息吗？谢谢

shiyu

@shiyu

Latest posts made by shiyu