From 70c9905fec86887f077a481111a9d2a507bb3acf Mon Sep 17 00:00:00 2001
From: tkurz <theresa.kurz@iws.uni-stuttgart.de>
Date: Fri, 27 Mar 2020 12:11:51 +0100
Subject: [PATCH] [brookscorey] fix wrong implementation of dkrn_dswe
---
dumux/material/fluidmatrixinteractions/2p/brookscorey.hh | 6 ++----
1 file changed, 2 insertions(+), 4 deletions(-)
diff --git a/dumux/material/fluidmatrixinteractions/2p/brookscorey.hh b/dumux/material/fluidmatrixinteractions/2p/brookscorey.hh
index 22574dd3f8..404b080c60 100644
--- a/dumux/material/fluidmatrixinteractions/2p/brookscorey.hh
+++ b/dumux/material/fluidmatrixinteractions/2p/brookscorey.hh
@@ -269,10 +269,8 @@ public:
swe = min(max(swe, 0.0), 1.0); // the equation below is only defined for 0.0 <= sw <= 1.0
- return 2.0*(swe - 1)*(1 + pow(swe, 2.0/params.lambda())*(1.0/params.lambda() + 1.0/2
- - swe*(1.0/params.lambda() + 1.0/2)
- )
- );
+ return 2.0*(swe - 1)*(1 + (0.5 + 1/params.lambda())*pow(swe, 2.0/params.lambda())
+ - (3.0/2.0 + 1/params.lambda())*pow(swe, 2.0/params.lambda()+1));
}
};
--
GitLab