fluent发散求助

金鱼化艺123

利用fluent自带的DDPM模型模拟煤粉的弥散过程，计算区域：280mm*1000mm，划分网格数13996，网格质量：Minimum Orthogonal Quality =  0.688090  Maximum Aspect Ratio =  8.72361，喷嘴区域网格加密。导入fluent先稳态计算连续相（气相，初始速度10-50m/s）的流场1000步，作为后续计算的初始值，改用非稳态注入颗粒相，采用相间耦合，考虑湍流效果，颗粒平均直径0.1mm，初始速度2-5m/s，质量流量为0.5kg/s。求解器设置：压力-速度耦合 phase couple simple，时间和空间离散均采用一阶迎风格式，松弛因子均减小为0.2，时间步长为0.0001S，计算发散。。。 这是要闹哪样
提示错误:
cx-xy-plot-data：unequal length argument lists
#Divergence detected in AMG solver：pressure correction-->turning off correction scaling!
#Divergence detected in AMG solver：pressure correction-->increasing relaxation sweeps
#Divergence detected in AMG solver：pressure correction
#Divergence detected in AMG solver：K
#Divergence detected in AMG solver：epsilon
.....
Primitive  Error at node 1： floating point exception
Primitive  Error at node 2： floating point exception
Primitive  Error at node 3： floating point exception
....
Error：floating point exception
Error Object：#f
求好心人，熟悉这方面错误的大牛，做fluent模拟的同胞给出建议。
先谢谢大家了

金鱼化艺123

自己先顶一下，这个问题折磨我半个多月了，由于错误一直进步不下去

xrs333

你确定颗粒的质量流量有那么大么？

金鱼化艺123

xrs333 发表于 2018-3-13 22:02
你确定颗粒的质量流量有那么大么？

确实颗粒的质量流量有点大，我尝试过0.04kg/s，但是还是发散

xrs333

应该是由于颗粒与连续相耦合，颗粒进入一个单元时与连续相交换的动量太大，造成动量方程的迭代不收敛。可以尝试把喷管附近的单元尺寸加粗，而不是加密。压力-速度耦合用PRESTO!。

金鱼化艺123

xrs333 发表于 2018-3-13 22:44
应该是由于颗粒与连续相耦合，颗粒进入一个单元时与连续相交换的动量太大，造成动量方程的迭代不收敛。可以 ...

谢谢，尝试一下。

xrs333

也许还需要减小时间步长。
建议把松弛因子改回为原值。

xrs333

增加颗粒流数目。这样，每一颗粒流的质量流量减小。

金鱼化艺123

金鱼化艺123 发表于 2018-3-14 08:29
谢谢，尝试一下。

将喷嘴周围网格加粗，选择DDPM模型后压力速度耦合没有PRESTO！，只有simple、simplec、couple、PISO，我选择的simple，计算还是发散

金鱼化艺123

xrs333 发表于 2018-3-14 14:55
增加颗粒流数目。这样，每一颗粒流的质量流量减小。

颗粒我采用的面入射类型

xrs333

可以发一些图上来，也许有线索。

金鱼化艺123

xrs333 发表于 2018-3-14 16:44
可以发一些图上来，也许有线索。

DDPM模型

湍流模型的选择

DPM模型的设置

颗粒的初始化入射状况

曳力模型的选择

气相湍动对颗粒的影响

颗粒动力学模型

压力速度耦合

松弛因子

步长

网格

xrs333

喷嘴口处气流速度多大？

xrs333

发散后还能显示计算结果吗？如能，比较一下前后的流场，看看是哪里出问题了。

xrs333

改用双精度也可能有改善。