将ICEM网格转化的时候，在0文件出现了cellToRegion子文件，这个有什么影响吗？

我假设了一个0.6×0.6×0.6的立方体，取其顶部为入口，选用可压缩浮力求解器buoyantPimpleFoam对其进行计算。奇怪的是，方腔内部的压力并没有升高？
我猜测其可能与参考压力和参考位置的设置有关，但是为什么可压缩求解器还需要设置参考压力？求解的结果显示参考压力最终形成了低压区。
参考点设置

PIMPLE { momentumPredictor yes; nOuterCorrectors 1; nCorrectors 2; nNonOrthogonalCorrectors 0; pRefCell 0; pRefValue 1e5; }

压力场
微信截图_20220404173941.png
算例文件，Allrun后输入buoyantPimpleFoam即可运行
链接：https://pan.baidu.com/s/1vv2--60uf43oKf2CKg27DQ
提取码：nk3q
感觉这问题还挺基础的，这是什么问题导致的计算结果错误呢:135:

我们在paraview中进行后处理提取某个截面受力，使用clip把某个面截出来，然后计算calculator计算受力，然后plot data over time, 但是这样如果在计算的时候保存的时间步长间隔比较大算的受力变化很不精确，有没有办法在计算中去计算某个截面的受力。
能不能通过代码把某个截面表示出来，然后使用force function object去计算。
或者是paraview中有没有显示类似boundary文件中网格编号的功能，手动在boundary文件中把截面的网格编号拿出来设置一个新的面，可以直接在function object中计算

@李东岳 李老师想请教一下，如果想获得网格信息，比如面网格坐标，是不是前面必须要有如下代码：

微信截图_20220402174442.png
谢谢李老师！

我自己暂时的理解是，程序中fjfj.phase()/fj.x()和fkfk.phase()/fk.x()指代的是j和k的颗粒数目，也就是右边公式中的Nj和Nk。但是其中的fj.phase和fk.phase是什么含义就有点不太明白了。还有就是聚并速率的量纲是m3/s的，程序里面为什么还要再乘以i颗粒体积fi.x()？这样子的话量纲就不是每秒了。

各位老师好，

我用openfoam4.x - pimpleDyMFoam - sixDoFMotionSolver计算刚体运动，我的几何由body1和body2两部分组成，两个body都是刚体，相连并同时运动，因此也可以看做一个刚体。我目前将其看做一个刚体，因此在dynmaicMeshDict文件中写作：

dynamicFvMesh dynamicMotionSolverFvMesh; motionSolverLibs ( "libsixDoFRigidBodyMotion.so" ); solver sixDoFRigidBodyMotion; sixDoFRigidBodyMotionCoeffs { ... patches (BODY); ... }

计算后能得到刚体上的升阻力，但是我要如何根据现有结果分别得到两个body上的升阻力呢？（如果我理解正确的话，BODY1上的气动力 + BODY2上的气动力 = BODY上的气动力）

还是说只能将其看做两个刚体并重新计算，即dynamicMeshDict中改为patches (BODY1 BODY2);

谢谢！

@冠竹simplefoam、pimplefoam、rhosimplefoam和rhopimplefoam的层流计算是不是也可以称为准直接数值模拟呢？

在v2012版本中，TEqn.H中有rhoCpPhi的定义为：

const surfaceScalarField rhoCpPhi ( "rhoCpPhi", rhoPhi*(Cp1 - Cp2) + phi*rho2*Cp2 );

这个定义是如何得到的? 如果Cp1 = Cp2，那rhoCpPhi = phi*rho2*Cp2，这个和 phi*rho*Cp有差距。

如果直接把中心点的rhoCp插值到面上和phi相乘应该表示为 fvc::interpolate(rho1*alpha1*Cp1 + rho2*alpha2*Cp2)*phi。

@杨英狄 这个也是计算粒子与壁面的作用，所以就是在particleerosion函数的基础上改的

@李东岳 李老师，fluent和其他机器计算时长的对比有吗

@cccrrryyy 让运维换了个of2106版本，速度有所提升，相当于自己电脑的四分之一速度，离期望值还是低点

@李东岳 烦请东岳老师有空看下，期待您的回复。谢谢！

@bestucan 谢谢大佬，已经解决了

看下tutorials里的Allrun和Allmesh操作步骤

@wenke 你好，你这个问题解决了吗？我也遇到了这种问题，现在也没有解决。

@qing 是这样的，openfoam完全可以这样任意折腾。它是基于文件的，只要能读取、识别到相应的文件就可以进行计算。

我思考了下，

matrix relax 虽然具体实现方式不同，但是本质不变。

你说的是对的。一楼我说的应该有问题（我更正一下）。不仅仅field relax，matrix relax，也会影响计算结果。我当时把matrix relax理解成不会影响计算结果了。

matrix relax是使得计算结果更慢的趋向于真实值，但换来一个对角占优矩阵

field relax是使得计算结果更慢的趋向于真实值

两个都使得计算结果更慢的趋向于真实值。

当时看这个公式的时候：$\phi^{n+1}=\phi^n+\beta(\phi^{n+1}-\phi^n)$，以为$\beta$变小时候，$\phi^{n+1}$还是原来那个值。然而并不是:136: :136: 好尴尬

@李东岳 李老师，您好！现在使用foam-ext 4.0,已经编译完成，计算成功一次，后续计算时一直报错（可以划分网格，计算中读取信息时出现问题）。
bed30e22-ed0d-4079-b2df-86b3002c91bf-image.png
7502d3ba-1a3d-4601-bd2a-3de6d249d09f-image.png

@zhouxu 谢谢，我尝试下。

@李东岳 东岳老师，请问这个速度能不能从10-2下降到10-4或者-6？