make it possible with only one phase in genneral?

This commit is contained in:
hnil
2019-05-20 20:30:40 +02:00
committed by Atgeirr Flø Rasmussen
parent 9c0520aeeb
commit 8572e1d108

View File

@@ -348,7 +348,11 @@ public:
if(FluidSystem::phaseIsActive(oilPhaseIdx)){ if(FluidSystem::phaseIsActive(oilPhaseIdx)){
x = rockCompressibility*(fluidState_.pressure(oilPhaseIdx) - rockRefPressure); x = rockCompressibility*(fluidState_.pressure(oilPhaseIdx) - rockRefPressure);
}else{ }else{
x = rockCompressibility*(fluidState_.pressure(waterPhaseIdx) - rockRefPressure); if(FluidSystem::phaseIsActive(waterPhaseIdx)){
x = rockCompressibility*(fluidState_.pressure(waterPhaseIdx) - rockRefPressure);
}else{
x = rockCompressibility*(fluidState_.pressure(gasPhaseIdx) - rockRefPressure);
}
} }
porosity_ *= 1.0 + x + 0.5*x*x; porosity_ *= 1.0 + x + 0.5*x*x;
} }