Yes, I have a denser mesh structure near free surface. But residual of Z-Momentum is diverging. I tried solving with the second order but it is the same. Thank you for your interest. Looking forward for your reply.