Simulation of mildly unsaturated flow