From 1d004a29c5d05e411434f213f297244ceae9e7b7 Mon Sep 17 00:00:00 2001
From: Benjamin Faigle <benjamin.faigle@posteo.de>
Date: Thu, 25 Nov 2010 15:18:44 +0000
Subject: [PATCH] fixed bug where flash yield wrong result in 1p cases
git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@4780 2fb0f335-1f38-0410-981e-8018bf24f1b0
---
dumux/decoupled/2p2c/dec2p2cfluidstate.hh | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/dumux/decoupled/2p2c/dec2p2cfluidstate.hh b/dumux/decoupled/2p2c/dec2p2cfluidstate.hh
index 1328e174d6..72e7d05eda 100644
--- a/dumux/decoupled/2p2c/dec2p2cfluidstate.hh
+++ b/dumux/decoupled/2p2c/dec2p2cfluidstate.hh
@@ -72,7 +72,7 @@ public:
* - determination of maximum solubilities (mole fractions) according to phase pressures
* - comparison with Z1 to determine phase presence => phase mass fractions
* - round off fluid properties
- * \param Z1 Feed mass fraction [-]
+ * \param Z1 Feed mass fraction: Mass of comp1 per total mass [-]
* \param pw Pressure of the wetting phase [Pa]
* \param poro Porosity [-]
* \param temperature Temperature [K]
@@ -122,11 +122,13 @@ public:
{
nu_[nPhaseIdx] = 1; // only nPhase
massfrac_[nPhaseIdx][wCompIdx] = Z1; // hence, assign complete mass soluted into nPhase
+ massfrac_[wPhaseIdx][wCompIdx] = 1.;
}
else // (Z1 >= Xw1) => no nPhase
{
nu_[nPhaseIdx] = 0; // no second phase
massfrac_[wPhaseIdx][wCompIdx] = Z1;
+ massfrac_[nPhaseIdx][wCompIdx] = 0.;
}
// complete array of mass fractions
--
GitLab