Use a modified version of Darcy's law which allows for blocking of spurious flows.
#include <dumux/porousmediumflow/2p1c/model.hh>
| Public Types | |
| using | type = TwoPOneCDarcysLaw<TypeTag> | 
| using Dumux::Properties::AdvectionType< TypeTag, TTag::TwoPOneCNI >::type = TwoPOneCDarcysLaw<TypeTag> |