On the importance of hysteresis and heterogeneity in the numerical simulation of unsaturated flow