OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Merge pull request #1055 from andlaus/fix_rock_compressibility

Browse Source 
fix incorrect derivative  of rock compressibility w.r.t. pressure
This commit is contained in:
Bård Skaflestad 2016-07-20 11:32:19 +02:00 committed by GitHub
commit bb387e16f9
1 changed files with 2 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
opm/core/props/rock/RockCompressibility.cpp
View File

@ -101,7 +101,8 @@ namespace Opm
if (p_.empty()) { if (p_.empty()) {
// Approximating poro multiplier with a quadratic curve, // Approximating poro multiplier with a quadratic curve,
// we must use its derivative. // we must use its derivative.
return rock_comp_ + 2 * rock_comp_ * rock_comp_ * (pressure - pref_); const double cpnorm = rock_comp_*(pressure - pref_);
return rock_comp_ + cpnorm*rock_comp_;
} else { } else {
return Opm::linearInterpolationDerivative(p_, poromult_, pressure); return Opm::linearInterpolationDerivative(p_, poromult_, pressure);
} }