关于cfdem计算发散

[chapaofan]

使用cfdempiso算法，之前计算空气的时候能够正常计算，结果还可以，但是把流体改成水之后，就发生了两种粒子一旦接触就会丢失。我觉得应该是发散了，但是发散的话我应该从哪里修改入手？我尝试修改了阻力模型、流体速度，还有耦合的A、B、Bfull模型，都有这个发散问题。一般造成发散的原因都有哪些呢？

[chapaofan]

只要有浮力就会发散，为什么呢

[吴建民]

@chapaofan 您好，请问问题现在解决了吗。我也遇到类似问题

[李东岳]

描述的比较粗糙，需要提供一下更详细的信息

[chapaofan]

@吴建民 我有点忘记了，是解决掉了，明天我给您看一下我之前的解决方法，一般来说主要问题是出在网格上边，网格尺寸需要比粒子尺寸大，一旦网格不对，两种不同大小粒子接触就会出现这类问题，可以尝试把网格尺寸调大一些尝试一下

[吴建民]

@chapaofan 好的，麻烦了

[chapaofan]

@吴建民 而且水的粘度很大，对于入口流速大的模型，在遇到比较小的球体会造成球体速度突然增大到指数级增长，这个时候时间步长就会相对不足了，就会导致发散，可以尝试降低一下液体流速。我也不清楚这种是不是能解决这个问题

[吴建民]

@chapaofan 我目前水流流速是0.2m/s，我猜问题主要在颗粒部分，我的颗粒是0.1mm的，大概有十几万个。OF的时间步是10的-5，LIGGGHTS的是10的-6，算一会了就发散。

[吴建民]

@吴建民 我大幅度细化了一下网格，没有任何作用

[chapaofan]

@吴建民 对对，我算过类似的算例，你可以提取一下这些小球收到拖拽力之后的速度变化量，我记得当时我看的是变化幅度太大，因为你进口0.2m/s你的小球只有2e-4m的直径他们两个相距太大了会导致它速度增长量太大，就会导致发散，目前我摸索的可能解决办法是减小时间步长或者降低流速。而且我建议你先在liggghts里设定一个0.2mm的小颗粒让他计算一下，看能不能计算不发散。

[吴建民]

@chapaofan 嗯嗯好的，感谢。我明天再试试

[chapaofan]

@吴建民 不是细化，你的网格大小要大于你的粒子，要粗糙一点，只有IB算法好像是要求网格大小小于粒子，其他都要求网格大小大于粒子

[吴建民]

@chapaofan 也没那么细，我的网格还是比颗粒大的

[chapaofan]

@吴建民
这个确实太早我忘记原因了，我回去翻了以下原来的，他好像是动能实在是太大，原本kineng在0.几，之后增加到了几万，在碰撞接触后就会导致粒子丢失，全部丢失后就会发散。我之前是把速度从3减到0.2好像解决了，但我之后没遇到过这类问题。我也刚学cfdem不久，希望能帮到你。

[吴建民]

@chapaofan 谢谢你。我昨天说错了，我原本的CFD时间步是1e-4，我昨天夜里调成1e-5让它跑，早上看，算了几万步，库朗数稳在0.1左右，目前是不发散了。之后用服务器再试试看。再次感谢您的回复