煤燃烧模型的收敛性问题

paul2002

我模仿cfx煤燃烧模型的算例，算了一个锅炉的煤燃烧问题，4800个粒子，但总是不收敛，也不发散，特别是焓方程。
我减小了松弛因子，效果也不明显，不知道有什么解决方法吗？
谢谢。

ggbaby

CFX认为那个就是收敛了，因为不停的有源相耦合进来。

paul2002

不会吧，残差很大！
我设置的URF如下：
>>SOLVER DATA
  >>ROGRAM CONTROL
    MAXIMUM NUMBER OF ITERATIONS 200
    MASS SOURCE TOLERANCE 1.0000E-04
    ITERATIONS OF VELOCITY AND PRESSURE EQUATIONS 5
    ITERATIONS OF TEMPERATURE AND SCALAR EQUATIONS 2
    ITERATIONS OF TURBULENCE EQUATIONS 2
  >>UNDER RELAXATION FACTORS
    U VELOCITY 3.0000E-01
    V VELOCITY 3.0000E-01
    W VELOCITY 3.0000E-01
    TEMPERATURE 1.0000E+00
    ENTHALPY 8.0000E-01
    ALL MASS FRACTIONS 0.6
    ALL COMBUSTION SCALARS 8.0000E-01
  >>COURANT NUMBER FALSE TIMESTEPS
   U VELOCITY 1.0
   V VELOCITY 1.0
   W VELOCITY 1.0
   ENTHALPY 5.0
另外，我看了CFX文档solution.pdf里面一大堆关于收敛，叠代的东西，看得稀里糊涂，不知，ggbaby斑竹能不能总结一下？

paul2002

另外，我还想问一下网格划分是不是对收敛问题影响很大，我的几何模型并不复杂，也没有用局部网格种子。
多谢！！

ggbaby

>>COURANT NUMBER FALSE TIMESTEPS
  U VELOCITY 1.0
  V VELOCITY 1.0
  W VELOCITY 1.0
  ENTHALPY 5.0
把这个去掉。这个在你的那个里面没有作用。

paul2002

刚才算了一下CFX的煤燃烧模型的例子，发现它也是收敛很慢，UVW,MASS的残差基本上在1E0到1E-2之间振荡，焓方程残差更大。那又怎么判断结果的正确性呢？CFX是用质量和能量平衡来判断的，不知道这样是否就可以了。
另外还想向版主请教一个问题：
fo文件里的粒子计算部分总是出现：
PARTICLE NO.   233  DIAMETER =  1.229E-04  MASS =  1.253E-09  MASS FLOW RATE =  8.300E-03
***** INTEG: MORE THAN    2 REPEATED CONVERGENCE FAILURES FROM ONESTP.
几乎所有粒子都这样，不知道为什么会这样，会不会影响结果正确性？
多谢，多谢！！！

ggbaby

网格不行，所以可以没有收敛。呵呵。。