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

Use well index for segment hydrostatic pressure drop

Browse Source
This commit is contained in:
Joakim Hove
2021-05-27 08:28:09 +02:00
parent ee1a1c6afc
commit 960de96d75
3 changed files with 10 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@@ -1924,7 +1924,7 @@ namespace Opm
// TODO: we might be able to add member variables to store these values, then we update well state
// after converged
const auto hydro_pressure_drop = getHydroPressureLoss(seg);
well_state.segPressDropHydroStatic()[top_segment_index + seg] = hydro_pressure_drop.value();
well_state.segPressDropHydroStatic(index_of_well_)[seg] = hydro_pressure_drop.value();
pressure_equation -= hydro_pressure_drop;
if (frictionalPressureLossConsidered()) {
@@ -2706,7 +2706,8 @@ namespace Opm
auto seg_press_drop = well_state.segPressDrop(index_of_well_);
auto seg_press_friction = well_state.segPressDropFriction(index_of_well_);
seg_press_drop[seg] = well_state.segPressDropHydroStatic()[seg] +
auto seg_press_static = well_state.segPressDropHydroStatic(index_of_well_);
seg_press_drop[seg] = seg_press_static[seg] +
seg_press_friction[seg] +
well_state.segPressDropAcceleration()[seg];
}

2
opm/simulators/wells/WellState.cpp
View File

@@ -1158,7 +1158,7 @@ WellState::reportSegmentResults(const PhaseUsage& pu,
auto& segpress = seg_res.pressures;
segpress[Value::Pressure] = this->segPress(well_id)[seg_ix];
segpress[Value::PDrop] = this->segPressDrop(well_id)[seg_ix];
segpress[Value::PDropHydrostatic] = this->segPressDropHydroStatic()[seg_dof];
segpress[Value::PDropHydrostatic] = this->segPressDropHydroStatic(well_id)[seg_ix];
segpress[Value::PDropFriction] = this->segPressDropFriction(well_id)[seg_ix];
segpress[Value::PDropAccel] = this->segPressDropAcceleration()[seg_dof];
}

10
opm/simulators/wells/WellState.hpp
View File

@@ -241,14 +241,16 @@ public:
return &seg_pressdrop_friction_[top_segment_index];
}
std::vector<double>& segPressDropHydroStatic()
double* segPressDropHydroStatic(std::size_t well_index)
{
return seg_pressdrop_hydorstatic_;
const int top_segment_index = this->top_segment_index_[well_index];
return &seg_pressdrop_hydorstatic_[top_segment_index];
}
const std::vector<double>& segPressDropHydroStatic() const
const double* segPressDropHydroStatic(std::size_t well_index) const
{
return seg_pressdrop_hydorstatic_;
const int top_segment_index = this->top_segment_index_[well_index];
return &seg_pressdrop_hydorstatic_[top_segment_index];
}
std::vector<double>& segPressDropAcceleration()