@李东岳 嗯嗯，这个应该没问题了，现在解决方案就是reconstructPar就效果很好了
谢谢李老师！

最近涉及到一些化学反应，并不十分理解。在比如reactingTwoPhaseEulerFoam求解器的tutorial，bubbleColumnEvaporatingReacting中，展示了所有物质都在气相状态下的反应，CO^0.93 + H2O^0.24 = CO2^0.69 + H2^1。因为该反应被定义在reactions.gas文件中，反应物和产物的属性也在thermo.gas文件中，所有物质都为气态。

我想问问如果我的反应产物是有液态或者固态参与，比如A(气态) + B(气态) = C(气态) + D(液态)，应该如何设置这个reaction和thermo等文件。有没有相关的教程或者资料可以参考？

@李东岳 :xinxin:

@李东岳 解决问题了 谢谢李老师！

wallshear要怎么在tecplot显示啊

@学流体的小明 在 toposet之后出现空集。。。 中说：

感觉是topoSet选择网格，提取各种面的时候有些问题了，我没用过这方面的功能，不太了解。
你如果只是想把边界面提取出来的话，在blockMesh当中就可以设置吧。

@李东岳 在 toposet之后出现空集。。。 中说：

豁牙子注意一下toposet点位置的精度。

感谢两位大佬提点，豁牙子问题解决了...

@田畔的风 感谢您的回复 可以咨询一下您如果我使用mppic方法时还可以怎么样发生粒子相变吗 我目前想实现水滴入射时发生蒸气相变

@李东岳 谢谢哥，现在OK了

@李东岳 感谢老师回复，我的是课题组服务器，应该是有48核心的。我用GPT教的测试并行c程序是可以成功运行的，但是openfoam每个案例都有并行问题。不知道怎么搞

有关于dsmcfoam使用的相关教程吗？我运行官方的案例是可以的，但是有些文件/代码不知道具体代表什么。并且我想模拟一个稀薄气体的流场，以及稀薄气体吹小的固体颗粒的demo，不知道如何去设置边界条件等，有人知道怎么做吗？谢谢谢谢！！！

@vien锁 您好，看到您的帖子，最近在学习dsmcfoam的使用，可以向你请教相关问题吗？谢谢！

您好，请问您这个问题解决了吗？能分享一下解决办法吗？

有没有前辈使用过ParaView脚本编辑图片，需要实现自定义图例和坐标轴，超过4位数的坐标和图例刻度值需要表示为例如：10,000样式，即三位一个逗号分隔。请问各位有没有实现过？

https://github.com/louisgag/openFoam-Overset-SpeedUp

@学流体的小明 是的，这个算例中确实使用了动网格，使用reconstructParMesh后就成功了

谢谢李老师！我再去仔细看下

viv振幅突然特别大导致发散是什么原因呢？一开始振幅几厘米，后边突然特别大

@zhou 您好，我也是做这些方面的 可以留一个联系方式方便进行交流嘛 我刚开始接触这部分

@zzkluck 天哪，还真是...细节了

@Mina-Lee 在 precice tourorials中的受力问题 中说：

@RolandLent precice中没有计算OF力的能力，力都是在OF里计算的。OF中的力也是通过force函数计算的。
pecice只干三件事:
1、引用of中的力变量，提取流固耦合交界面的力变量；引用位移变量。2、precice内部运算。3、覆盖引用的位移变量。

我阅读了第一位回答者的源代码，力的计算是precice的of适配器实现的，我在适配器的forcebase类里边找到了计算力的方法，这里附上计算压力和黏性力的部分

if (solverType_.compare("incompressible") == 0) { forceField.boundaryFieldRef()[patchID] = surface * pb[patchID] * rhob[patchID]; } else if (solverType_.compare("compressible") == 0) { forceField.boundaryFieldRef()[patchID] = surface * pb[patchID]; } forceField.boundaryFieldRef()[patchID] += surface & devRhoReffb[patchID];

如果是用of的solid求解器就是向您这样的执行方式了

// Check if a force field with the requested name exists. // If yes (e.g., solids4Foam), bind Force_ to that field. // If not (e.g., pimpleFoam without the Forces function object), create it. if (mesh_.foundObject<volVectorField>(nameForce)) { Force_ = &const_cast<volVectorField&>( mesh_.lookupObject<volVectorField>(nameForce)); } else { ForceOwning_.reset(new volVectorField( IOobject( nameForce, mesh_.time().timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), mesh, dimensionedVector( "fdim", dimensionSet(1, 1, -2, 0, 0, 0, 0), Foam::vector::zero))); Force_ = ForceOwning_.get(); }