Preferential flow behaviour in unsaturated packed beds and heaps: Incorporating into a CFD model