Continuity does not converge
 

Hi people!

I am currently doing a very simple case on natural convection. Basically, I will have a fluid domain from a piping. The 2 faces will have different temperature (hot and cold face) and the circumferential face will be adiabatic.

I am using k-epsilon model with enhanced wall treatment together with buoyancy effect and curvature correction. My y+ value is set to be 5. Y value is calculated from the max velocity ontained from the natural convection and ANSYS Mesh is used.

I have tried several time but I cannot seem to have the continuity converged. May I know what are some of the possible issue which i can look into?

Thank you for your time.