Ok I will change the viscous model from k-epsilon to laminar.nFor the density definition I used Piecewise-linear and I introduced 2 points depending the temperature (K) (910 kg/m3 at 313 K and 790 kg/m3 at 398 K).nI don't have any natural convection. All faces are adiabatic. The only one boundary condition that I have in my model is a heat flux in a wall inside the deposit.nWhat do you think?n