绝热过程这个实现起来非常简单，等温过程我看也类似，过几天可以试一下
方程我都写在这里了，你可以自己尝试植入一下
http://dyfluid.com/comadiabatic.html

@xue 这个链接里使用的系泊模型是准静态系泊模型，我测试了浮式平台的自由衰减，精度不理想

03a93b4e-34b4-4d36-8d29-72da1776f74d-image.png 大佬们在用wave2foam求解波浪与浮体作用时，alpha water出现负值及大于1的情况是为什么呢，时间步长0.001，如果是网格问题的话该如何进行修改呢

按理说颗粒应力的作用是阻止颗粒体积分数过大，但是就目前的结果看添加了颗粒应力，颗粒反而更加聚集，体积分数更大

才看见这个帖子。看起来是颗粒应力并没有起作用。不过半年过去了，估计已经解决了。

想提取截面上随时间变化的流场数据，然后放到另外一个流场中作为入流条件，但没找到合适的方法，请问大家有什么建议吗

我自己一般在本地服务器做计算，超算上传下载比较麻烦。
内存1TB算1000万网格有点浪费吧，李老师卖的服务器，内存256G，我也能算2700万的网格，8000万网格的也不会爆。大网格情况下，CPU限制了计算速度，内存是确保程序能运行下去的。

还有1000万网格时需要开到64核以上计算吗？

可以往上开，计算速度应该会更快。可以先研究一下加速比再做决定。

waveProperties.input文件如下所示：
80357645591c3eb0d80586e29884d5e.png

打算在OpenFOAM里提取计算域内一条线上的所有网格。目前想到了几个方法：

通过topoSet来做，做一个box，恰好包围一个狭长的矩形，然后创建zone，这个zone可以包含所有线上的网格。难点就是这个矩形的点不太好找，要固定的特别精细，稍微大一点，网格就多了。

用paraview挑出来某条线上的网格ID，简单，就是麻烦。网格一多要疯。

给一系列线上的点，通过nearestCell判断网格ID。但是寻找Cell比较慢。并且这线上的一系列点也需要给一堆。

哪位大佬有什么更好的方式，能把一条线上的网格ID都给提取出来么。

@EricLiu 是纯气体

太高端了 知识盲区了老铁

瞬态的话网格比较细的话，时间步长应该足够小

这个问题是在编译的时候，postProcessingWaves没有被编译，导致编译postProcessWvaes2Foam中无法正确使用该函数。

@李东岳 嗯嗯，这个应该没问题了，现在解决方案就是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了