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

added: MultisegmentWellAssemble::assembleTrivialEq

Browse Source 
extracted from MultisegmentWellEval::assembleICDPressureEq
This commit is contained in:
Arne Morten Kvarving 2022-11-18 12:09:43 +01:00
parent 05a4ca85a7
commit 32dce644d3
3 changed files with 17 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
opm/simulators/wells/MultisegmentWellAssemble.cpp
View File

@ -200,6 +200,16 @@ assemblePressureEq(const int seg,
} }
} }
template<class FluidSystem, class Indices, class Scalar>
void MultisegmentWellAssemble<FluidSystem,Indices,Scalar>::
assembleTrivialEq(const int seg,
const Scalar value,
Equations& eqns) const
{
eqns.resWell_[seg][SPres] = value;
eqns.duneD_[seg][seg][SPres][WQTotal] = 1.;
}
#define INSTANCE(...) \ #define INSTANCE(...) \
template class MultisegmentWellAssemble<BlackOilFluidSystem<double,BlackOilDefaultIndexTraits>,__VA_ARGS__,double>; template class MultisegmentWellAssemble<BlackOilFluidSystem<double,BlackOilDefaultIndexTraits>,__VA_ARGS__,double>;

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

@ -98,6 +98,11 @@ public:
bool wfrac = has_wfrac_variable, bool wfrac = has_wfrac_variable,
bool gfrac = has_gfrac_variable) const; bool gfrac = has_gfrac_variable) const;
//! \brief Assembles a trivial equation.
void assembleTrivialEq(const int seg,
const Scalar value,
Equations& eqns) const;
private: private:
const WellInterfaceIndices<FluidSystem,Indices,Scalar>& well_; //!< Reference to well const WellInterfaceIndices<FluidSystem,Indices,Scalar>& well_; //!< Reference to well
}; };

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

@ -1404,8 +1404,8 @@ assembleICDPressureEq(const int seg,
if (const auto& segment = this->segmentSet()[seg]; if (const auto& segment = this->segmentSet()[seg];
(segment.segmentType() == Segment::SegmentType::VALVE) && (segment.segmentType() == Segment::SegmentType::VALVE) &&
(segment.valve().status() == Opm::ICDStatus::SHUT) ) { // we use a zero rate equation to handle SHUT valve (segment.valve().status() == Opm::ICDStatus::SHUT) ) { // we use a zero rate equation to handle SHUT valve
linSys_.resWell_[seg][SPres] = this->primary_variables_evaluation_[seg][WQTotal].value(); MultisegmentWellAssemble<FluidSystem,Indices,Scalar>(baseif_).
linSys_.duneD_[seg][seg][SPres][WQTotal] = 1.; assembleTrivialEq(seg, this->primary_variables_evaluation_[seg][WQTotal].value(), linSys_);
auto& ws = well_state.well(baseif_.indexOfWell()); auto& ws = well_state.well(baseif_.indexOfWell());
ws.segments.pressure_drop_friction[seg] = 0.; ws.segments.pressure_drop_friction[seg] = 0.;