残差很奇怪，各位老师们帮忙分析可能的错误

OpenFOAM
1

IMG_20200601_215057
IMG_20200601_2150574000×3000 1.49 MB

2

边界条件
alpha. gas 一次风 二次风 fixedValue 1
返料口 fixedValue 0
wall 出口 零梯度
alpha. particles . 一次风二次风 fixedValue 0
返料口 fixedValue 1
wall 出口 零梯度上传中：U.gas IMG_20200601_215655.jpg…
U. particles

IMG_20200601_215748
IMG_20200601_2157483000×4000 2.37 MB

p
IMG_20200601_215855
IMG_20200601_2158553000×4000 2.88 MB

p_rgh
IMG_20200601_215921
IMG_20200601_2159213000×4000 3.02 MB

T. gas
IMG_20200601_220012
IMG_20200601_2200123000×4000 2.55 MB

T. particles
IMG_20200601_220044
IMG_20200601_2200443000×4000 2.2 MB

其他的k nut epslion 都是仿照reactingTwoPhaseEulerFoam的bubbleEvaluationReacting算例，没有改动

3

IMG_20200601_220457
IMG_20200601_2204573000×4000 1.94 MB

算例是这样子的流化床锅炉

4

其实这个也算正常，你看看每一步是不是达到最大迭代次数了？残差也没有扩大，看看流场怎么样？

5

只看残差收敛量级马马虎虎能接受，嫌它波动大，可以试试加密网格，调小时间步长，检查一下反应模型，选一个精度稍高一点的湍流模型 :joy:

6

这个结果还是可以看，精度不高，残差量级在 10^{-3} 量级，还算可以。

如果要提高精度，减小残差波动，可以考虑楼上同学的建议，还可以

  1. 在 system/fvSolution 里面增加迭代的最大次数。比如：
    p
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-06;
        relTol          0.05;
        maxIter         2000;  // 改这里，如果没有这一项，可以加。没有这一项的话，默认值是 1000
    }
  1. 网格可以调整下，使其过渡平滑，均匀，不要有很大的长宽比的网格。网格质量也是比较重要的。
7

谢谢老师回复 不好意思这个帖子有点乱 我在另一个新帖子整理的边界和出现的问题了

8

谢谢回复，10e-4就算的很慢了我先去检查网格质量

9

谢谢老师 ，我去试一下