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

Merge pull request #5952 from hnil/consistent_primary_variables

Browse Source 
Consistent primary variables
This commit is contained in:
Atgeirr Flø Rasmussen 2025-02-13 18:24:53 +01:00 committed by GitHub
commit c015df56f9
No known key found for this signature in database
GPG Key ID: B5690EEEBB952194
2 changed files with 5 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -152,6 +152,7 @@ update(const WellState<Scalar>& well_state,
} }
} }
} }
init();
} }
template<class FluidSystem, class Indices> template<class FluidSystem, class Indices>
@ -208,6 +209,7 @@ updateNewton(const BVectorWell& dwells,
if (stop_or_zero_rate_target) { if (stop_or_zero_rate_target) {
value_[0][WQTotal] = 0.; value_[0][WQTotal] = 0.;
} }
init();
} }
template<class FluidSystem, class Indices> template<class FluidSystem, class Indices>

3
opm/simulators/wells/StandardWellPrimaryVariables.cpp
View File

@ -223,6 +223,7 @@ update(const WellState<Scalar>& well_state,
// BHP // BHP
value_[Bhp] = ws.bhp; value_[Bhp] = ws.bhp;
init();
} }
template<class FluidSystem, class Indices> template<class FluidSystem, class Indices>
@ -301,6 +302,7 @@ updateNewton(const BVectorWell& dwells,
// so that bhp constaint can be an active control when needed. // so that bhp constaint can be an active control when needed.
constexpr Scalar bhp_lower_limit = 1. * unit::barsa - 1. * unit::Pascal; constexpr Scalar bhp_lower_limit = 1. * unit::barsa - 1. * unit::Pascal;
value_[Bhp] = std::max(value_[Bhp] - dx1_limited, bhp_lower_limit); value_[Bhp] = std::max(value_[Bhp] - dx1_limited, bhp_lower_limit);
init();
} }
template<class FluidSystem, class Indices> template<class FluidSystem, class Indices>
@ -320,6 +322,7 @@ updateNewtonPolyMW(const BVectorWell& dwells)
value_[pskin_index] -= relaxation_factor * dx_pskin; value_[pskin_index] -= relaxation_factor * dx_pskin;
} }
} }
init();
} }
template<class FluidSystem, class Indices> template<class FluidSystem, class Indices>