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

add PhasePresence, make muOil() function happy.

Browse Source
This commit is contained in:
Liu Ming 2014-09-26 14:05:17 +08:00
parent 8620ce2641
commit 1cb1e8475a
1 changed files with 7 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
opm/polymer/fullyimplicit/utilities.cpp
View File

@ -245,12 +245,18 @@ namespace Opm
src[cell] = 1.0; src[cell] = 1.0;
} }
} }
//Add PhasePresence make muOil() happy.
std::vector<PhasePresence> phaseCondition;
for (int c = 0; c < num_cells; ++c) {
phaseCondition[c].setFreeWater();
phaseCondition[c].setFreeOil();
}
const Selector<double> src_selector(src); const Selector<double> src_selector(src);
const V one = V::Constant(num_cells, 1.0); const V one = V::Constant(num_cells, 1.0);
const V zero = V::Zero(num_cells); const V zero = V::Zero(num_cells);
const std::vector<V> kr = props.relperm(sw, so, zero, cells); const std::vector<V> kr = props.relperm(sw, so, zero, cells);
const V muw = props.muWat(p, cells); const V muw = props.muWat(p, cells);
const V muo = props.muOil(p, zero, cells); const V muo = props.muOil(p, zero, phaseCondition, cells);
const V krw_eff = polymer_props.effectiveRelPerm(c, cmax, kr[0]); const V krw_eff = polymer_props.effectiveRelPerm(c, cmax, kr[0]);
const V inv_muw_eff = polymer_props.effectiveInvWaterVisc(c, muw.data()); const V inv_muw_eff = polymer_props.effectiveInvWaterVisc(c, muw.data());
std::vector<V> mob(np); std::vector<V> mob(np);