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

Fix output of relative permeability for sequential model.

Browse Source 
Bug caused assert() failures for debug mode runs.
This commit is contained in:
Atgeirr Flø Rasmussen
2017-04-21 15:20:10 +02:00
parent edea16dd74
commit 90d90fb452
1 changed files with 5 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

9
opm/autodiff/BlackoilTransportModel.hpp
View File

@@ -303,14 +303,15 @@ namespace Opm {
// Compute and store mobilities. // Compute and store mobilities.
const int canonical_phase_idx = canph_[ phase_idx ]; const int canonical_phase_idx = canph_[ phase_idx ];
const ADB& phase_pressure = state.canonical_phase_pressures[canonical_phase_idx]; const ADB& phase_pressure = state.canonical_phase_pressures[canonical_phase_idx];
const ADB mu = asImpl().fluidViscosity(canonical_phase_idx, phase_pressure, state.temperature, state.rs, state.rv, cond); sd_.rq[phase_idx].mu = asImpl().fluidViscosity(canonical_phase_idx, phase_pressure, state.temperature, state.rs, state.rv, cond);
sd_.rq[phase_idx].kr = kr[canonical_phase_idx];
// Note that the pressure-dependent transmissibility multipliers are considered // Note that the pressure-dependent transmissibility multipliers are considered
// part of the mobility here. // part of the mobility here.
sd_.rq[ phase_idx ].mob = tr_mult * kr[phase_idx] / mu; sd_.rq[ phase_idx ].mob = tr_mult * sd_.rq[phase_idx].kr / sd_.rq[phase_idx].mu;
// Compute head differentials. Gravity potential is done using the face average as in eclipse and MRST. // Compute head differentials. Gravity potential is done using the face average as in eclipse and MRST.
const ADB rho = asImpl().fluidDensity(canonical_phase_idx, sd_.rq[phase_idx].b, state.rs, state.rv); sd_.rq[phase_idx].rho = asImpl().fluidDensity(canonical_phase_idx, sd_.rq[phase_idx].b, state.rs, state.rv);
const ADB rhoavg = ops_.caver * rho; const ADB rhoavg = ops_.caver * sd_.rq[phase_idx].rho;
sd_.rq[ phase_idx ].dh = ops_.grad * phase_pressure - rhoavg * gdz; sd_.rq[ phase_idx ].dh = ops_.grad * phase_pressure - rhoavg * gdz;
if (is_first_iter_) { if (is_first_iter_) {