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

adding getPressureEq() to MultisegmentWell

Browse Source
This commit is contained in:
Kai Bao 2017-09-12 13:16:57 +02:00
parent 2bf82b4262
commit d30c1eb2c5
2 changed files with 40 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
opm/autodiff/MultisegmentWell.hpp
View File

@ -45,6 +45,8 @@ namespace Opm
// TODO: for now, not considering the polymer, solvent and so on to simplify the development process. // TODO: for now, not considering the polymer, solvent and so on to simplify the development process.
// TODO: should I begin with the old primary variable or the new fraction based variable systems? // TODO: should I begin with the old primary variable or the new fraction based variable systems?
// Let us begin with the new one // Let us begin with the new one
// TODO: we need to have order for the primary variables and also the order for the well equations.
// sometimes, they are similar, while sometimes, they can have very different forms.
enum WellVariablePositions { enum WellVariablePositions {
GTotal = 0, GTotal = 0,
WFrac = 1, WFrac = 1,
@ -317,6 +319,8 @@ namespace Opm
std::vector<EvalWell>& mob) const; std::vector<EvalWell>& mob) const;
EvalWell getControlEq() const; EvalWell getControlEq() const;
EvalWell getPressureEq(const int seg) const;
}; };
} }

41
opm/autodiff/MultisegmentWell_impl.hpp
View File

@ -245,7 +245,7 @@ namespace Opm
// considering the contributions from the inlet segments // considering the contributions from the inlet segments
{ {
for (const int inlet : segment_inlets_[seg]) { for (const int inlet : segment_inlets_[seg]) {
for (int comp_idx; comp_idx < num_comp; ++comp_idx) { for (int comp_idx = 0; comp_idx < num_comp; ++comp_idx) {
const EvalWell inlet_rate = getSegmentRate(inlet, comp_idx); const EvalWell inlet_rate = getSegmentRate(inlet, comp_idx);
resWell_[seg][comp_idx] -= inlet_rate.value(); resWell_[seg][comp_idx] -= inlet_rate.value();
for (int pv_idx = 0; pv_idx < numWellEq; ++pv_idx) { for (int pv_idx = 0; pv_idx < numWellEq; ++pv_idx) {
@ -304,10 +304,18 @@ namespace Opm
// the fourth dequation, the pressure drop equation // the fourth dequation, the pressure drop equation
// if it is the top segment, it should be the well control equations // if it is the top segment, it should be the well control equations
// if it is not, it will be the pressure drop equation // if it is not, it will be the pressure drop equation
if (seg != 0) { // not the top segment {
; // TODO: currently, we only handle the hydrostatic pressure difference.
} else { // the top segment // We need to add the friction pressure loss and also the acceleration pressure loss
; // with the acceleration pressure loss, there will be inlets flow rates (maybe alos the oulet flow)
// not sure whether to handle them implicitly or explicitly
// TODO: we can try to handle them explicitly first, if it does not work, we can handle them
// implicitly. Even explicily, we can calculate them without considering the derivative first
const EvalWell control_eq = getControlEq();
resWell_[seg][SPres] = control_eq.value();
for (int pv_idx = 0; pv_idx < numWellEq; ++pv_idx) {
invDuneD_[seg][seg][SPres][pv_idx] = control_eq.derivative(pv_idx + numEq);
}
} }
} }
} }
@ -1311,4 +1319,27 @@ namespace Opm
return control_eq; return control_eq;
} }
template<typename TypeTag>
typename MultisegmentWell<TypeTag>::EvalWell
MultisegmentWell<TypeTag>::
getPressureEq(const int seg) const
{
// for top segment, the well control equation will be used.
if (seg == 0) { // for top segment, the well control equation will be used.
return getControlEq();
}
const EvalWell pressure_equation = getSegmentPressure(seg);
const int outlet_segment_location = numberToLocation(segmentSet()[seg].outletSegment());
const EvalWell outlet_pressure = getSegmentPressure(outlet_segment_location);
pressure_equation -= outlet_pressure;
// we need to handle the pressure difference between the two segments
// we only consider the hydrostatic pressure loss first
return pressure_equation;
}
} }