Use existing formula for compressibility effect.

This commit is contained in:
Atgeirr Flø Rasmussen
2022-07-04 14:31:26 +02:00
parent 8543e677c1
commit 0013c49aa2

View File

@@ -515,8 +515,7 @@ public:
} else {
x = rockCompressibility*(fluidState_.pressure(gasPhaseIdx) - rockRefPressure);
}
//porosity_ *= 1.0 + x + 0.5*x*x;
porosity_ *= 1.0 + (1 + 0.5*x)*x;
porosity_ *= 1.0 + x + 0.5*x*x;
}
// // deal with water induced rock compaction
// porosity_ *= problem.template rockCompPoroMultiplier<Evaluation>(*this, globalSpaceIdx);