二维瞬态蓄热，内有相变层，产生自然对流

你好，我想问的是一个融化凝固，蓄热的模型。具体示意图如下，实际计算过程简化为2D！网格数3.2W。
边界是恒热流加热，然后其他壁面绝热。等加热到一段时间，相变材料会融化，而且会产生自然对流现象。需要加热8小时。
我先说下我的问题，本来时间步长为0.2S时候可以收敛。可以正常的算到八小时。
后来我就把材料的密度，小改了一下，结果速度曲线就一直平行，不收敛。算到900S还是一直平行。（原来算到200秒左右的时候速度残差取消就下降，然后慢慢收敛了）