CFD中文网

@欧阳 我以前 OpenFOAM-8 用的上面发的代码，好像没问题的话

@CT 很遗憾，没有。我现在转去了别的方向，也不再使用该求解器了。如果只是希望做DSMC模拟而不是必须使用OpenFOAM的话，可以尝试其它开源软件，如Sparta等。

我不是很理解。不过如果你不想写laplacian，你可以写div(grad())

是的，sensibleEnthalpy是h，sensibleInternalEnergy是e
不过这俩应该不影响rhoCentralFoam，他能量方程自己硬植入了

@Amadeus
以及此图片下方B.2节中第二行

因此关于这部分的介绍非常多,本节重复介绍

是否应该为“本节不再重复介绍”?

您好，请问我用fluent计算这个算例，为什么整体gas hold up相差一个数量级。入口气体速度设置的0.0213m/s。设置基本和论文里相同。
Snipaste_2025-03-24_20-32-07.png eb5ceff2-0be5-4eb1-acbb-826f3f732768-image.png

snappyHexMeshDict文件里面

locationInMesh (x y z);

这一项里面的点坐标要设置为你需要生成网格的部分内，即这个点的坐标在你stl文件内则生成stl文件形状的网格，若点坐标在stl文件外面，则生成背景网格和stl之间部分的网格.

这个要是openfaom的网格，就好处理了。2亿网格确实太多了，icem我也只画过最多几百万的网格。卡的都不行。你这2亿网格太疯狂了，机器内存应该很大。

你这个细化背景网格可以考虑采用topoSet和refineMesh一起进行细化，不一定需要采用SnappyHexMesh。感觉好像snappyHexMesh在重叠网格中主要是运用在floatingBody哪儿作为提取网格特性和细化的工具，在背景网格中还是主要采用topoSet和refineMesh

这里如果使用SMALL，当接近零的时候，Y_s的值会接近零，滑移速度为0，也是符合原理的。但是我不确定OpenFOAM在计算时是否Y_s会真的趋向于0，请大佬们看看是否有更好的处理方法

感谢李老师！

有大佬帮助查看一下case，调一下收敛吗。有偿感谢

@李东岳 响应老师号召:mihu:

cylinder.tar.xz

sdf asdf

@李东岳 东岳老师，您有什么好的解决思路吗？

@李东岳 感觉这个问题挺有意义的，肯定可以从原理上说清楚为啥发散。给的是一个极大的单块网格平动（入口速度零，零法向压力梯度；出口压力零，零速度梯度），一直调不出来就放弃了，可能雷诺数低到某个程度，就不能靠放大计算域来使边界条件匹配了。

祝李老师上课愉快！

@李东岳 确实

compute_node_loop 不能在HOST节点运行，你117行#if RP_HOST 规定了下面的代码在HOST节点里，但是还要进行compute node 循环，所以会卡住，可能？

哦哦 好的，感谢大佬分享这么多资料！