I had ran the solver before with different porous values and it got converged. In EWF settings we expect the condensation to happen and it did happen in the initial case. But I tried running the solver with EWF off and it is still diverging. Mesh has been improved and minimum orthogonal quality is now above 0.1.