Convergence of the porous medium model
 

My question is 2-D model, gas flow field. When simulating, firstly, we use first-order upwind scheme, Once the solution approached a stable condition, we switched a second-order upwind scheme. However, the residuals are still high to 0.01. How can we obtain the converged results for the porous medium model? Thanks in advance.