Skip to content
Snippets Groups Projects
Commit 1d004a29 authored by Benjamin Faigle's avatar Benjamin Faigle
Browse files

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
parent 53810381
No related branches found
No related tags found
No related merge requests found
...@@ -72,7 +72,7 @@ public: ...@@ -72,7 +72,7 @@ public:
* - determination of maximum solubilities (mole fractions) according to phase pressures * - determination of maximum solubilities (mole fractions) according to phase pressures
* - comparison with Z1 to determine phase presence => phase mass fractions * - comparison with Z1 to determine phase presence => phase mass fractions
* - round off fluid properties * - 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 pw Pressure of the wetting phase [Pa]
* \param poro Porosity [-] * \param poro Porosity [-]
* \param temperature Temperature [K] * \param temperature Temperature [K]
...@@ -122,11 +122,13 @@ public: ...@@ -122,11 +122,13 @@ public:
{ {
nu_[nPhaseIdx] = 1; // only nPhase nu_[nPhaseIdx] = 1; // only nPhase
massfrac_[nPhaseIdx][wCompIdx] = Z1; // hence, assign complete mass soluted into nPhase massfrac_[nPhaseIdx][wCompIdx] = Z1; // hence, assign complete mass soluted into nPhase
massfrac_[wPhaseIdx][wCompIdx] = 1.;
} }
else // (Z1 >= Xw1) => no nPhase else // (Z1 >= Xw1) => no nPhase
{ {
nu_[nPhaseIdx] = 0; // no second phase nu_[nPhaseIdx] = 0; // no second phase
massfrac_[wPhaseIdx][wCompIdx] = Z1; massfrac_[wPhaseIdx][wCompIdx] = Z1;
massfrac_[nPhaseIdx][wCompIdx] = 0.;
} }
// complete array of mass fractions // complete array of mass fractions
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment