@Samuel-Tu 用了，确实快，但是有时候分网不成功，而且网格质量很难保证

本来想说试试加深水槽。 静水流深嘛，一看是浅水方程……
如果边界条件设定为带反射的，那就应该不停波动。如果是自由边界，不停波动不正常

在cfd online 里逛了一下，有些收货。
似乎non-orthogonal corrector 是在组装矩阵Ax=b的时候用的，它会把扩散项中的正交部分放入系数矩阵Ａ中，而把非正交部分放入方程右边的源项b中，使得系数矩阵更稀疏且主对角线控制，使得计算更稳定？？cfd online

@cccrrryyy 受教受教，授我以渔，万分感谢~ 帖子链接里开篇就阐明了我的疑惑，也难怪k-epsilon可以在工程领域运用的如此广泛~ 我就先对照着openfoam tutorials和CFDEM案例里的k-epsilon取值作参考，计算后再做调整，再次感谢指点~

@东岳 非常感谢老师:xinxin: ！ :xiexie:
就是这样迂回的话好像有些难，我还要再思考思考，再尝试尝试 :144:

@东岳 好滴

Already solved.

@bestucan Thanks for your attention.

谢谢老师

@东岳 这个当时手一抖就点上了，发现没法删掉，真是尴尬:mihu:

@莫晓柒 编译之后得到的库和求解器程序都在 platforms 下面，所以只能重新安装了。

@东岳 我不算运营人员，只是参与了号主的一些活动，有时候帮帮忙。。。

@upc_ngh
1.是的，$\mathbf{U}$代表的是KinematicParcel中的$\mathbf{U_c}$，也就是连续相的在颗粒所处位置的速度插值，具体的插值方式OF里面有cell，cellpoint那些，$\mathbf{U}_\rm{d}$表示的就是颗粒的速度。
2.是的，系数乘以速度差。

新知识点

@沐风 因为是百度云，所以没有下载:xinlei: 不过我有几个小建议：
计算域是否够大，圆柱和上下前后边界的距离是否够远，去除边界的影响；
snGrad建议用corrected；
我印象中上下边界应该用symmetry，用noslip是有什么原因或者依据吗。

@xpqiu 感谢

@chien 是的谢谢🙏，放弃编译了😹

@小考拉 maxU会突然增大只是发散的迹象，不是发散的原因，所以这种问题很难debug :135:

我也是每次都重新设置，很麻烦，前来学习

min(alpha.water) 都出现负值了。

时间步自己调整的特别小是为了让库朗数降低，结果时间步那么小，库朗数还那么大，可能网格有问题