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

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

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

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

问题：如何在自定义边界条件中获取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进行计算，不知具体如何操作？

谢谢大家！

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

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

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

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

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

请问可以在运行过程中，即 边界运算里面提取y方向上一条线的网格中心信息（速度之类的），谢谢

Hi，东岳大佬，请教下openfoam里可以实现分区网格计算，比如两个分区，并设置interface boundary condition来coupling这两个区域吗？如果可以，麻烦给点思路，谢谢。
我暂时想到的只有fluid-solid interaction，即FSI toolbox in Foam-extend，改写solid domain为fluid-fluid interaction，来实现，不过有些繁琐了。

BR，
Shiyu