computing the density of segment mixture and pressure drop earlier

the density of segmetn mixture will be used in the well flux
calculation.
This commit is contained in:
Kai Bao 2015-10-14 14:59:55 +02:00
parent ba7b3e728c
commit d5f917d1ac

View File

@ -516,6 +516,9 @@ namespace Opm {
return; return;
} }
asImpl().computeSegmentDensities(state);
asImpl().computeSegmentPressuresDelta(state);
std::vector<ADB> mob_perfcells; std::vector<ADB> mob_perfcells;
std::vector<ADB> b_perfcells; std::vector<ADB> b_perfcells;
asImpl().extractWellPerfProperties(mob_perfcells, b_perfcells); asImpl().extractWellPerfProperties(mob_perfcells, b_perfcells);
@ -530,8 +533,6 @@ namespace Opm {
V aliveWells; V aliveWells;
std::vector<ADB> cq_s; std::vector<ADB> cq_s;
asImpl().computeWellFlux(state, mob_perfcells, b_perfcells, aliveWells, cq_s); asImpl().computeWellFlux(state, mob_perfcells, b_perfcells, aliveWells, cq_s);
asImpl().computeSegmentDensities(state, well_state);
asImpl().computeSegmentPressuresDelta(state);
asImpl().updatePerfPhaseRatesAndPressures(cq_s, state, well_state); asImpl().updatePerfPhaseRatesAndPressures(cq_s, state, well_state);
asImpl().addWellFluxEq(cq_s, state); asImpl().addWellFluxEq(cq_s, state);
asImpl().addWellContributionToMassBalanceEq(cq_s, state, well_state); asImpl().addWellContributionToMassBalanceEq(cq_s, state, well_state);