From cc28be8dfbf53b47bc0f92b3c4581bedd3d2f54e Mon Sep 17 00:00:00 2001
From: Simon Scholz <simon.scholz@iws.uni-stuttgart.de>
Date: Wed, 4 Apr 2018 10:57:24 +0200
Subject: [PATCH] [2p2c] unify 2p2c and 2pnc init saturation guess
The 2p2c and 2pnc priVar-switch had different initial saturation
guesses. It was decided to unify them to the symmetry of the 2pnc
priVarswitch.
---
dumux/porousmediumflow/2p2c/primaryvariableswitch.hh | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/dumux/porousmediumflow/2p2c/primaryvariableswitch.hh b/dumux/porousmediumflow/2p2c/primaryvariableswitch.hh
index 3b9282e919..703da6cb14 100644
--- a/dumux/porousmediumflow/2p2c/primaryvariableswitch.hh
+++ b/dumux/porousmediumflow/2p2c/primaryvariableswitch.hh
@@ -107,9 +107,9 @@ protected:
<< xww + xwn << std::endl;
newPhasePresence = bothPhases;
if (formulation == pnsw)
- priVars[switchIdx] = 0.0;
+ priVars[switchIdx] = 0.0001;
else if (formulation == pwsn)
- priVars[switchIdx] = 1.0;
+ priVars[switchIdx] = 0.9999;
}
}
else if (phasePresence == wPhaseOnly)
@@ -135,9 +135,9 @@ protected:
<< xnw + xnn << std::endl;
newPhasePresence = bothPhases;
if (formulation == pnsw)
- priVars[switchIdx] = 0.999;
+ priVars[switchIdx] = 0.9999;
else if (formulation == pwsn)
- priVars[switchIdx] = 0.001;
+ priVars[switchIdx] = 0.0001;
}
}
else if (phasePresence == bothPhases)
--
GitLab