你这个本身就是写的源代码，你自己改一下不就是了

你的tank水位是不是在增加

我也遇到一模一样的情况，请问解决了吗

openFoam VOF 如何监测进出口的流量和压差

https://gitlab.com/hfchen20/foamMooring

这个你们可以尝试一下，或许我这次讲课回来可以把这个打包到虚拟机里面

b768ef65-777f-4eff-827d-64e7647225e3-image.png ca8f10bf-0a79-46a0-868e-3508b1d64b78-image.png两相流模型为VOF，流场按照浸入边界法对流速场进行了插值，颗粒内部流速和颗粒速度一样为0，但是还是有初始颗粒内部虚拟流体流出，各位大佬有什么解决办法么:135:

其实也不能按照csdn的来。他那个把mpi换了，但是openfoam如果用mpich的话，需要重新编译openfoam。不过楼主那个问题相同算例直接计算可以得到收敛结果，但并行运算时算到大概一半时间时报错，并行计算有时候是这样，单核能跑，多核报错，不太好debug

@htdddwg 感谢回复，我去看看

@李东岳 好的，感谢东岳老师

@zju_klh
尝试在controlDict中增加一行看看，添加lib，其中XXX替换为编译模型库的名字

libs ("libXXX.so");

感谢大家，应该是跑通了，主要是没有使用wedge边界的原因。
用ICEM二维旋转小角度生成三维网格很容易出现报错，使用of7的makeAxialMesh工具实现楔形体网格很方便。
7758ff58-19b6-4f68-8562-ba11e109c72d-46250a307d5087f20b5092c195baecd.png

2楼回帖没看懂。不过你1楼的回帖，是把方程处理成$...=(-\phi)+(\phi)$的形式，$(-\phi)$隐性离散增加对角占优

可以！没问题！这个也挺有意思，可以深入交流！

各位大佬，本人想用rigidbodymotion计算双浮体铰接，铰的位置打算设置在两个浮体中间，但是目前的设置只能铰在两个浮体重心位置，我的浮体1重心位置在全局坐标系中坐标为（0，0，-0.0126），浮体2重心位置在全局坐标系中坐标为（0.5，0，-0.0126），打算把铰设置在在全局坐标系中（0.25，0，-0.0126）的位置，我该如何实现，以下是我的设置1.jpg 2.png 3.png

我先用一维气液激波管试试看看会不会有显著的震荡，谢谢李老师详细的解答~

在OpenFOAM11之后，跨音速求解器rhoCentralFoam（或者说求解模块shockFluid）的求解变量变成了原始变量而非守恒变量。我还没有详细的看算法，粗略的看是OpenFOAM基金会自己做的数值改动。

这知识啊，无穷无尽。

新版：https://github.com/OpenFOAM/OpenFOAM-11/tree/master/applications/modules/shockFluid

HLL类：http://dyfluid.com/hll.html

CFD: 可压 + 密度基：http://dyfluid.com/density.html

@尚善若水 收到，感谢！

@biaowai 不会，跟那个没关系

@李东岳 清华最近有工作“2024-工热学报-基于物理信息神经网络的气膜冷却湍流模型反演学习”，用的Si-LU激活函数:）

原来是这样！谢谢！