Numerical solution of steady-state porous flow free boundary problems