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

Remove segment pressure drop as separate variable

Browse Source
This commit is contained in:
Joakim Hove 2021-05-30 18:58:17 +02:00
parent 1ea70da7e6
commit 2d67d085ef
3 changed files with 5 additions and 21 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -2701,14 +2701,6 @@ namespace Opm
const UnitSystem& unit_system = ebosSimulator.vanguard().eclState().getDeckUnitSystem(); const UnitSystem& unit_system = ebosSimulator.vanguard().eclState().getDeckUnitSystem();
assemblePressureEq(seg, unit_system, well_state, deferred_logger); assemblePressureEq(seg, unit_system, well_state, deferred_logger);
} }
auto seg_press_drop = well_state.segPressDrop(index_of_well_);
auto seg_press_friction = well_state.segPressDropFriction(index_of_well_);
auto seg_press_static = well_state.segPressDropHydroStatic(index_of_well_);
auto seg_press_accel = well_state.segPressDropAcceleration(index_of_well_);
seg_press_drop[seg] = seg_press_static[seg] +
seg_press_friction[seg] +
seg_press_accel[seg];
} }
} }

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

@ -539,7 +539,6 @@ void WellState::init(const std::vector<double>& cellPressures,
} }
//seg_rates_ = wellRates(); //seg_rates_ = wellRates();
seg_rates_.assign(nw*np, 0); seg_rates_.assign(nw*np, 0);
seg_pressdrop_.assign(nw, 0.);
seg_pressdrop_hydorstatic_.assign(nw, 0.); seg_pressdrop_hydorstatic_.assign(nw, 0.);
seg_pressdrop_friction_.assign(nw, 0.); seg_pressdrop_friction_.assign(nw, 0.);
seg_pressdrop_acceleration_.assign(nw, 0.); seg_pressdrop_acceleration_.assign(nw, 0.);
@ -946,7 +945,6 @@ void WellState::initWellStateMSWell(const std::vector<Well>& wells_ecl,
assert(int(seg_press_.size()) == nseg_); assert(int(seg_press_.size()) == nseg_);
assert(int(seg_rates_.size()) == nseg_ * numPhases() ); assert(int(seg_rates_.size()) == nseg_ * numPhases() );
seg_pressdrop_.assign(nseg_, 0.);
seg_pressdrop_hydorstatic_.assign(nseg_, 0.); seg_pressdrop_hydorstatic_.assign(nseg_, 0.);
seg_pressdrop_friction_.assign(nseg_, 0.); seg_pressdrop_friction_.assign(nseg_, 0.);
seg_pressdrop_acceleration_.assign(nseg_, 0.); seg_pressdrop_acceleration_.assign(nseg_, 0.);
@ -1159,7 +1157,7 @@ WellState::reportSegmentResults(const PhaseUsage& pu,
using Value = data::SegmentPressures::Value; using Value = data::SegmentPressures::Value;
auto& segpress = seg_res.pressures; auto& segpress = seg_res.pressures;
segpress[Value::Pressure] = this->segPress(well_id)[seg_ix]; segpress[Value::Pressure] = this->segPress(well_id)[seg_ix];
segpress[Value::PDrop] = this->segPressDrop(well_id)[seg_ix]; segpress[Value::PDrop] = this->segPressDrop(well_id, seg_ix);
segpress[Value::PDropHydrostatic] = this->segPressDropHydroStatic(well_id)[seg_ix]; segpress[Value::PDropHydrostatic] = this->segPressDropHydroStatic(well_id)[seg_ix];
segpress[Value::PDropFriction] = this->segPressDropFriction(well_id)[seg_ix]; segpress[Value::PDropFriction] = this->segPressDropFriction(well_id)[seg_ix];
segpress[Value::PDropAccel] = this->segPressDropAcceleration(well_id)[seg_ix]; segpress[Value::PDropAccel] = this->segPressDropAcceleration(well_id)[seg_ix];

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

@ -219,16 +219,12 @@ public:
return &seg_press_[top_segment_index]; return &seg_press_[top_segment_index];
} }
double * segPressDrop(std::size_t well_index) double segPressDrop(std::size_t well_index, std::size_t segment_index) const
{ {
const int top_segment_index = this->top_segment_index_[well_index]; const int top_segment_index = this->top_segment_index_[well_index];
return &seg_pressdrop_[top_segment_index]; return this->seg_pressdrop_friction_[top_segment_index + segment_index] +
} this->seg_pressdrop_hydorstatic_[top_segment_index + segment_index] +
this->seg_pressdrop_acceleration_[top_segment_index + segment_index];
const double * segPressDrop(std::size_t well_index) const
{
const int top_segment_index = this->top_segment_index_[well_index];
return &seg_pressdrop_[top_segment_index];
} }
double* segPressDropFriction(std::size_t well_index) double* segPressDropFriction(std::size_t well_index)
@ -537,8 +533,6 @@ private:
int nseg_; // total number of the segments int nseg_; // total number of the segments
// The following data are only recorded for output // The following data are only recorded for output
// pressure drop
std::vector<double> seg_pressdrop_;
// frictional pressure drop // frictional pressure drop
std::vector<double> seg_pressdrop_friction_; std::vector<double> seg_pressdrop_friction_;
// hydrostatic pressure drop // hydrostatic pressure drop