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

@xpqiu 感谢

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

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

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

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

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

尝试解释一下:mihu:

就是0.3图里的涡在0.4不见了，然后1里又出来了？

0.3里刚过柱体，柱体后面压力小，然后左右流体往中间挤，挤不动了，要么往上要么往下，往下是墙，只能往上冒，就是0.3里的涡

0.4里不见了，因为这个柱体下面有个空隙，流体在过空隙的时候很挤，过了空隙就释放压力，下边是墙，是能往上释放，然后空隙里出来的流体就把0.3里的涡往上托，托到液面托没了就成0.4了

然后0.6，为了释放压力往上托涡，但又因为惯性，释放过头了，涡把流体都往上带，下面就出现低压区了，然后左右的流体都往中间流，在这里碰头，可以从0.6看见他们碰头后很纠结，不知道往上走还是往下走

然后是1，往上的那部分变成了0.3的涡尾，往下的那部分，又卷出来个涡。

要是水足够深，应该能同时看见0.3和1的涡共存。

再形象点，就是0.3里的涡是柱体后的低压区造成。1里的涡就像柱体下的狭缝放的“屁”，一个屁把0.3里的涡挤走了……:yes:

计算区域为旋转区域与非旋转区域，颗粒由非旋转区域入射，但是后期处理颗粒轨迹时候发现，颗粒无法通过动静交界面（cyclicAMI）.
在颗粒求解文件中设置交界面为none
请问大家有碰到过这问题吗？有什么好的解决方法吗？

@桎梏 好的，我试一下，谢谢建议！

@bestucan 头文件里说的是如果非保守通量（ non-conservative fluxes，应该是翻译成非保守通量吧？）对算法的预测部分有影响的话是需要correctPhi的，但是什么情况下会影响，我完全没有概念。至于后半说多相流的部分，4.1多相流里不开correctPhi的还是有几个的。所以我完全没提取出什么有效信息。

可以用paraview直接算出自适应网格得出的球体体积，然后跟理论值对比一下，求出两者的相对误差就可以了

@波流力不好意思后来赶ddl耽搁了。后来我测试了下，您的算例运行的不错，液面形态很好。然后参考您的设置更改了我的计算，液面有提升，但是在撞到右侧壁，上升时候还是出现了裂隙。然后在您的算例里改成我的case的几何参数，裂隙现象就没了。。。很奇怪，我再研究看看:143:

@xiaolin leftEdge没有意义，leftObstacleEdge（rightObstacleEdge）是指运动的左（右）侧边界的位置
这个类其实是一个不太好用的类，你自己要改一下。这个建议你先去看一下源码，了解里面的操作流程，然后修改一下就可以自己用了，亲测可以

估计是一个吹出来的波，一个荡出来的波，叠一块了，想把荡出来的波滤掉，只看吹出来的波。

滤波器一般按频率滤波，不是按看见的特征滤波。
你花出来的曲线已经是叠加波了。如果不加吹的算出来的才是你要滤掉的部分。这部分应该很规律，能写成三角函数，那么用calculator 把物理量减去三角函数就行了，三角函数的自变量是x的坐标。

如果就要减去这个曲线图，得拟合出来这个函数，在 paraview 里处理才方便。

其他麻烦的方法倒是一大堆，paraview 里有 python。或者导出，用matlab之类的处理。

@东岳 我在速度文件力设置好内部速度为定值（0.6），压力文件里设置内部压强（1e5）通过一段时间计算后，流速越来越大。像我需要指定初始速速，只需要设置速度文件，压力文件中设置为0就可以，或者将指定的初始速度转化为相应的压力添加在压力文件中，速度文件设置为0呢!1.png 2.png

@xpqiu 感谢老师详细指点，试了下这个方法确实简洁:xiexie:

这是没有编译成功吧，我在foam-extend-4.0+2.4.0上面编译成功过这套程序

计算发散报错如下：
报错.jpg

@刀尔東 你好，我想问一下，在后处理中想对一个量做加减运算，这个在后处理里应该是什么指令呢？比如对速度U的模取反，postProcess -func "-mag (U)"，这样是报错的，不知道应该怎么处理？感谢:146:

@cccrrryyy 好的，感谢解答，我再思考一下:quwan: