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