peng-robinson test: correct calculation of the Rs factor

This commit is contained in:
Andreas Lauser
2014-07-08 12:39:12 +02:00
parent 17dbe13a7f
commit 42d0f895e9

View File

@@ -322,6 +322,9 @@ int main(int argc, char** argv)
surfaceAlpha,
flashFluidState,
/*guessInitial=*/false);
Scalar Rs =
surfaceFluidState.saturation(gasPhaseIdx)
/ surfaceFluidState.saturation(oilPhaseIdx);
std::cout << alpha << " "
<< flashFluidState.pressure(oilPhaseIdx) << " "
<< flashFluidState.saturation(gasPhaseIdx) << " "
@@ -329,7 +332,7 @@ int main(int argc, char** argv)
<< flashFluidState.density(gasPhaseIdx) << " "
<< flashFluidState.averageMolarMass(oilPhaseIdx) << " "
<< flashFluidState.averageMolarMass(gasPhaseIdx) << " "
<< surfaceFluidState.saturation(gasPhaseIdx)*surfaceAlpha << " "
<< Rs << " "
<< rho_gRef/flashFluidState.density(gasPhaseIdx) << " "
<< rho_oRef/flashFluidState.density(oilPhaseIdx) << " "
<< "\n";