OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-10 01:55:33 -06:00
Code Issues Packages Projects Releases Wiki Activity

Avoid call to saturationPressure

Browse Source
This commit is contained in:
Stein Krogstad 2023-12-19 09:54:52 +01:00
parent eec2fbaaa1
commit e908c95ca2
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
opm/simulators/wells/WellInterface_impl.hpp
View File

@ -1541,8 +1541,9 @@ namespace Opm
const double connection_pressure = ws.perf_data.pressure[perf];
// viscosity is evaluated at connection pressure
const auto& rv = getValue(intQuants.fluidState().Rv());
const double psat = FluidSystem::gasPvt().saturationPressure(this->pvtRegionIdx(), ws.temperature, rv);
const double mu = connection_pressure < psat ?
const double rv_sat = FluidSystem::gasPvt().saturatedOilVaporizationFactor(this->pvtRegionIdx(), ws.temperature, connection_pressure);
// note that rv here is from grid block with typically p_block > connection_pressure (so we may very well have rv > rv_sat)
const double mu = rv >= rv_sat ?
FluidSystem::gasPvt().saturatedViscosity(this->pvtRegionIdx(), ws.temperature, connection_pressure) :
FluidSystem::gasPvt().viscosity(this->pvtRegionIdx(), ws.temperature, connection_pressure, rv, getValue(intQuants.fluidState().Rvw()));
double rho = FluidSystem::referenceDensity(FluidSystem::gasPhaseIdx, this->pvtRegionIdx());