@prometheus10 好，谢谢

已经获得成功。
icem中创建的part对应着openFOAM中的boundarys
icem中创建的body对应着openfoam中的cellzones

各位前辈好，关于外流场计算中的动压、总压等有些疑问，想和大家请教一下：
对于不可压缩的流体，总压=静压+动压。OF的多相流计算中，初始条件中需要给一个参考压力p_rgh，然后经计算得到总压p。
那么，以一个初始位置在水位以上的结构物（即没有静水压力）为例，当结构物受到波浪冲击时，结构物所受到的流体压力可以认为全部是动压吗？或者说，OF计算出的在自由液面处的压力都是动压？

计算结果在WSL里还是？

请问pimpleFoam求解器是否可以对风机尾流场进行数值模拟？有算例或者文献可以提供参考吗？感谢感谢:xinxin: :xinxin: :xinxin:

如下图左边所示（channel395算例）
y+.png
我初步的想法是对yplus在Y方向上采样取一条直线，输出的数据文件第一列就是y轴坐标，第二列就是Y+，这样完成了转换，但是考虑到这个是平均速度剖面，只取样一条线的Y+可能没有那么强的代表性，需要算yPlusMean好像也不太对，请教各位改怎么处理呢？

@杨英狄 在 如何调用上一欧拉时间步的粒子信息 中说：

@星星星星晴 老师，比如当前时刻为4s，欧拉时间步是1s，我想在时刻为4s的时候调用时刻为3s的粒子属性。
我做的东西粒子半径是变化的，我需要用上一时刻的粒子半径一阶时间导数与当前时刻的粒子半径一阶时间导数进行比较。判断粒子的一些状态。
谢谢老师

或者拉格朗日开始和结束的时候做一下比较就好了

你这个要算多少万网格，算多少个算例。

1）云计算走租赁，几千块够了，不知道你要算多久

2）要么我租给你一个 :chouchou: 感兴趣的话邮件

3）因为我不做1万以下的机器，我做评论不太好，我也不方便评论。

@d-benjamin 你好，我想请教一下，你这种边界条件是如何设置的？

将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和其他机器计算时长的对比有吗