Evaluation saturatedViscosity(unsigned regionIdx,
const Evaluation& temperature,
const Evaluation& pressure) const
{
if (!enableThermalViscosity())
return isothermalPvt_->saturatedViscosity(regionIdx, temperature, pressure);
// compute the viscosity deviation due to temperature
const auto &muGasvisct = gasvisctCurves_[regionIdx].eval(temperature);
return muGasvisct;
}
Evaluation saturatedViscosity(unsigned regionIdx,
const Evaluation& temperature,
const Evaluation& pressure) const
{
const auto& isothermalMu = isothermalPvt_->saturatedViscosity(regionIdx, temperature, pressure);
if (!enableThermalViscosity())
return isothermalMu;
// compute the viscosity deviation due to temperature
const auto &muOilvisct = oilvisctCurves_[regionIdx].eval(temperature);
return muOilvisct/viscRef_[regionIdx]*isothermalMu;
}
Evaluation viscosity(unsigned regionIdx,
const Evaluation& temperature,
const Evaluation& pressure) const
{
const auto& isothermalMu = isothermalPvt_->viscosity(regionIdx, temperature, pressure);
if (!enableThermalViscosity())
return isothermalMu;
Scalar x = -pvtwViscosibility_[regionIdx]*(viscrefPress_[regionIdx] - pvtwRefPress_[regionIdx]);
Scalar muRef = pvtwViscosity_[regionIdx]/(1.0 + x + 0.5*x*x);
// compute the viscosity deviation due to temperature
const auto& muWatvisct = watvisctCurves_[regionIdx].eval(temperature);
return isothermalMu * muWatvisct/muRef;
}
