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

using addition to update segment pressure

Browse Source 
when bhp is changed. It is much safer and more reasonable.

using multiplication is more likely to result in dramatic values.
This commit is contained in:
Kai Bao 2022-03-09 14:14:57 +01:00
parent ae32b44124
commit ccfd4152d3
1 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
opm/simulators/wells/SegmentState.cpp
View File

@ -60,15 +60,16 @@ std::size_t SegmentState::size() const {
}
void SegmentState::scale_pressure(double bhp) {
void SegmentState::scale_pressure(const double bhp) {
if (this->empty())
throw std::logic_error("Tried to pressure scale empty SegmentState");
auto scale_factor = bhp / this->pressure[0];
const auto pressure_change = bhp - this->pressure[0];
std::transform(this->pressure.begin(),
this->pressure.end(),
this->pressure.begin(),
[scale_factor] (const double& p) { return p*scale_factor;});
[pressure_change] (const double& p) { return p + pressure_change;});
}
const std::vector<int>& SegmentState::segment_number() const {