用buoyantPimpleFoam求解器在kEqn亚格子LES计算进口速度0.21m/s的圆管流动，管径在3mm。div（phi，u）用的是Gauss LUST grad（U）格式，有壁面热流。入口处用的是mappedPatch，计算的泰勒尺度在0.0004左右，时间步设置为4e-5，开始计算收敛较好，计算一阵子后总会发散，而且查看云图会发现整根管的流体除入口处全部回流，请问这是什么原因啊，是离散格式的原因，还是偏离实际情况吗，我用稳态buoyantsimpleFoam计算可以很好收敛？
细节图片如下：
