Melting Problem
 

Hi there ... I want to model the melting process. however i don't know how to take the "natural convection" effect into account? (I have specified the gravity and volumetric expansion coefficient). Program converged after 3-4 iteration and results don't look fine. anybody kind enough to give me a hint! Tnx.