@李东岳 前处理只能在Open FOAM里面建模和画网格吗

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

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

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

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

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

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

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

snappyHexMeshDict文件里面

locationInMesh (x y z);

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

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

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

感谢李老师！

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

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

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

祝李老师上课愉快！

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

就是在存在Y=1的区域的时候，如果Ur不合理，会导致发散

这个帖子我已经忘了大体咋回事了。不过realizableKE有这个问题，你可以看看再尝试一下

https://cfd-china.com/topic/6911

但模拟结果显示的出口压力与想要的出口压力相差很大

waveTransmissive不能保证完全一致。如果是choke住那种的流动（不懂中文怎么说），直接给固定值。

@黒洛克 当时只是发现有这个问题。没有详细debug。感觉应该某个地方需要数值处理一下

可压缩流动瞬态算法要比稳态更容易收敛

@阿飞 过去4年了，今天我才发现foam-extend 4.1用的是绝热compressible VOF。我还进一步看了一下，OpeNFOAM-12还是用的T方程