CFD中文网

@李东岳 谢谢老师，看到了而且觉得效果可以的了，就准备试试这个方法了

目前研究的是溢洪道的空化计算，用的方法是三相的VOF加Schnerr-Sauer空化模型。
最近老师突然想让我考虑水流的含泥沙这一因素，请问有什么好的方法将泥沙这玩意考虑进模型里面呢？
是直接把泥沙考虑成一个流体相加到VOF里头（这样似乎就考虑不了粒径的影响）还是尝试用DPM和VOF耦合计算呢？（这样会不会计算量太大了，毕竟整体模型就挺大尺度的）。

又或者说可以通过一些边界条件将明渠流动化简为管道流然后再用interPhaseChangeFoam计算？

请问现在需要算一个明渠流动（溢洪道）的空化问题，openfoam里面有合适的求解器吗 。
这个算例子是个三相流动（水、空气、水蒸气）而常用来计算不可压空化问题的interPhaseFoam只能算两相流。
只是希望能够像VOF一样能够捕捉空化的位置就可以而不考虑类似空泡破裂这类的可压缩问题。

@youhaoyu 好的我试一试

想获得计算结果的某一个特定值（比如s = 1）的所有点，有啥比较好的方法呢？

@李东岳 谢谢老师

@李东岳 case.rar 老师，这是算例

@李东岳 老师，一旦抬升水位就回0.1s都算不完
多个水位都尝试过，只有之前尝试的2m左右能算到4s左右。

@李东岳 老师，我用setFields把进口后1m内的区域改为了水，可以算到2.5s然后就会出现同样的问题。

尝试了将阶前的部分区域初始化为水，可以算到4.1s出现同样问题

想请教下老师出现那个发散问题的原因是什么呢？初始的alpha场设置的不当还是初始的湍流参数设置不当导致的呢？

@李东岳 好的老师我再试试

@李东岳 blockMeshDict 老师改后的

@李东岳 case.rar 🙏

@李东岳 老师，我试了试将模型简化了，使得网格质量基本不会有问题了：

但是仍然回出现相同的问题，且数值也基本没怎么变化。。
是不是能确定是k和ω初始值的取值问题？
我也通过kOmega的计算器估算了这俩的值，没有解决问题，卡在这里了。。。

@李东岳 好的谢谢老师，我试试

是k和omega初始值的问题还是网格的问题呢？
因为网格也有部分skewness较大的网格：

用interFoam算一个明渠流动的算例，第一个时间步算完k后停止计算然后报错浮点溢出：

Floating point exception (core dumped)

请问大佬们这是个什么情况呢？

@李东岳 好的谢谢老师

烦求各位大佬讲解一下

在pointwise里画得二维结构网格，导入到openfoam里后，就发生了网格变形，边界也变了（比如斜边变成阶梯状），这要怎样解决呢？也没有找到之前别人帖子里的use VTK Polyhedron这个选项：
https://www.cfd-china.com/topic/5109/snappyhexmesh以及paraview显示问题