Use the model after Millington (1961) for the effective diffusivity.
#include <dumux/freeflow/navierstokes/mass/1pnc/model.hh>