Ansys Employee

That is the solution, hence my first post in the thread.   Look at  specifically  You're setting rho-rho0 to be zero, so lose the gravity term for the exterior boundaries, otherwise you need to set a pressure gradient on the exterior boundary to account for (rho-rho0)gh and that gets messy.