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

not extracting the mob and b inside computeWellFlux() ms wells

Browse Source
This commit is contained in:
Kai Bao 2015-11-19 09:14:18 +01:00
parent 4b9064a0f8
commit 57d796ca68
1 changed files with 2 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
opm/autodiff/BlackoilMultiSegmentModel_impl.hpp
View File

@ -721,8 +721,8 @@ namespace Opm {
template <class Grid>
void
BlackoilMultiSegmentModel<Grid>::computeWellFlux(const SolutionState& state,
const std::vector<ADB>& /* mob_perfcells */,
const std::vector<ADB>& /* b_perfcells */,
const std::vector<ADB>& mob_perfcells,
const std::vector<ADB>& b_perfcells,
V& aliveWells,
std::vector<ADB>& cq_s) const
{
@ -760,14 +760,6 @@ namespace Opm {
// const std::vector<int>& well_cells = well->wellCells();
const std::vector<int>& well_cells = wops_ms_.well_cells;
// extract mob_perfcells and b_perfcells.
std::vector<ADB> mob_perfcells(np, ADB::null());
std::vector<ADB> b_perfcells(np, ADB::null());
for (int phase = 0; phase < np; ++phase) {
mob_perfcells[phase] = subset(rq_[phase].mob, well_cells);
b_perfcells[phase] = subset(rq_[phase].b, well_cells);
}
// determining in-flow (towards well-bore) or out-flow (towards reservoir)
// for mutli-segmented wells and non-segmented wells, the calculation of the drawdown are different.
const ADB& p_perfcells = subset(state.pressure, well_cells);