Hello,
Is there a reason why the Brooks-Corey saturation curve is implemented
such that effective saturations > 1 are possible?
The function is implemented as
sw = (pc/pe)**(-lambda) for pc > 0
while for example in Bear & Cheng "Modeling Groundwater Flow and
Contaminant Transport" (2010) it is
sw = (pc/pe)**(-lambda) if pc > pe else 1
This restricts the effective saturation to be at most 1.
Best regards,
Samuel
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
