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

Added support for pressure-dependent viscosity to SinglePvtConstComp class.

Browse Source 
This is the class that handles PVTW and PVCDO inputs.
This commit is contained in:
Atgeirr Flø Rasmussen 2012-05-18 11:08:51 +02:00
parent bfcc856d4b
commit fac1c24fba
1 changed files with 15 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
opm/core/fluid/blackoil/SinglePvtConstCompr.hpp
View File

@ -50,16 +50,15 @@ namespace Opm
ref_B_ = pvtw[region_number][1]; ref_B_ = pvtw[region_number][1];
comp_ = pvtw[region_number][2]; comp_ = pvtw[region_number][2];
viscosity_ = pvtw[region_number][3]; viscosity_ = pvtw[region_number][3];
if (pvtw[region_number].size() > 4 && pvtw[region_number][4] != 0.0) { visc_comp_ = pvtw[region_number][4];
THROW("SinglePvtConstCompr does not support 'viscosibility'.");
}
} }
SinglePvtConstCompr(double visc) SinglePvtConstCompr(double visc)
: ref_press_(0.0), : ref_press_(0.0),
ref_B_(1.0), ref_B_(1.0),
comp_(0.0), comp_(0.0),
viscosity_(visc) viscosity_(visc),
visc_comp_(0.0)
{ {
} }
@ -68,11 +67,20 @@ namespace Opm
} }
virtual void mu(const int n, virtual void mu(const int n,
const double* /*p*/, const double* p,
const double* /*z*/, const double* /*z*/,
double* output_mu) const double* output_mu) const
{ {
std::fill(output_mu, output_mu + n, viscosity_); if (visc_comp_) {
// #pragma omp parallel for
for (int i = 0; i < n; ++i) {
// Computing a polynomial approximation to the exponential.
double x = -visc_comp_*(p[i] - ref_press_);
output_mu[i] = viscosity_/(1.0 + x + 0.5*x*x);
}
} else {
std::fill(output_mu, output_mu + n, viscosity_);
}
} }
virtual void B(const int n, virtual void B(const int n,
@ -132,6 +140,7 @@ namespace Opm
double ref_B_; double ref_B_;
double comp_; double comp_;
double viscosity_; double viscosity_;
double visc_comp_;
}; };
} }