Comparison of numerical and analytical models of unsaturated flow around underground openings

The authors benchmarked the integrated finite-volume code MULTIFLO against the analytical solution of Philip et al. (1989) for cylindrical underground openings. They found that the numerical simulation of dripping threshold flux was consistently higher than the analytical solution for rectangular and unstructured grids. The magnitude of this error was proportional to the alpha parameter characterizing capillarity. This error may bias estimates of seepage into drifts for repository performance towards nonconservatism.