interFoam求解突然浮点溢出

[火山口玩泥巴]

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

Floating point exception (core dumped)

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

[火山口玩泥巴]

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

[李东岳]

如果网格不好的话，最好多调调网格。网格质量好收敛就容易

[火山口玩泥巴]

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

[火山口玩泥巴]

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

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

[李东岳]

@火山口玩泥巴 如果你的网格少，且可公开，可以把算例处理干净发到这里我给你看看

[火山口玩泥巴]

@李东岳 case.rar 🙏

[李东岳]

网格太大了，算例跑不动，你把网格处理到2万以下我试试。

直接发修改后blockMesh的代码就行

[火山口玩泥巴]

@李东岳 blockMeshDict 老师改后的

[李东岳]

1）你的g方向不对 2）你的inlet跟atmosphere不要连起来，在inlet上面增加一个wall 3）初始化一下alpha场，不要让初始化的alpha拳是空气

然后你再试一下

[火山口玩泥巴]

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

[火山口玩泥巴]

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

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

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

[李东岳]

尝试了将阶前的部分区域初始化为水，

你试一下，不仅仅前部，后面也都是水，然后提升水位

[火山口玩泥巴]

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

[李东岳]

把你算例处理干净，打包发上来

[火山口玩泥巴]

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

[李东岳]

ddd.tar.xz

http://dyfluid.com/openfoam.html#id1

[火山口玩泥巴]

@李东岳 谢谢老师