Text this: Numerical model for NAPL migration in double-porosity subsurface systems