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");

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

@尚善若水 收到，感谢！

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

原来是这样！谢谢！

@李东岳 东哥，请问您现在知道原因了吗？，我也发现只有背景网格部分的流场信息被映射过来了，35cf3a02-93e9-41bc-837a-bcfe7f6003f3-8a4a6e0d657fca36881f8e9c71f54f2.png ，在进行运算的时候插值这一项也等于0，说明前景网格和背景网格没有进行插值运算。

@tidedrinker 请问用olaFlow在模拟高频波浪与潮流耦合时，总是遇到结构物水动力相应（x和y向受力）随频率增加而震荡的情况（在高频段的波浪工况下无法得到光滑的波浪频率与结构物受力的趋势线），这类情况是否是网格密度的原因？

@李东岳 在 paraview三维算例查看 中说：

paraview上面有一个功能，叫clip，你clip alpha=0.5的值试一下

谢谢东岳老师，确实能看到内部场了

@yanlei interFsiFoam求解器可以分享一下嘛？你说的弹性边界具体是什么，可以交流一下

这是一个非常好的问题。之前上课的时候也有人问我。OpenFOAM其实是不存在网格编号的。你可以看一下无痛笔记里面OpenFOAM的非结构网格。

@李东岳 好的，谢谢老师指点！

有的也不需要alpha。这个没有明确的定义。你可以乘、也可以不乘。