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

Move hydro p-loss assembly to end of function

Browse Source
This commit is contained in:
Stein Krogstad
2023-05-16 13:26:58 +02:00
parent a3860e7d0c
commit 51604e5956
1 changed files with 11 additions and 11 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@@ -225,23 +225,13 @@ assembleDefaultPressureEq(const int seg,
EvalWell extra_derivatives;
// we need to handle the pressure difference between the two segments
// we only consider the hydrostatic pressure loss first
// hydrostatic pressure loss is assembled seperately at the end
// TODO: we might be able to add member variables to store these values, then we update well state
// after converged
auto& ws = well_state.well(baseif_.indexOfWell());
auto& segments = ws.segments;
// We approximate the hydrostatic pressure loss over the segment by using the mean of the densities
// at the current segment node and its outlet node. Since density derivatives are organized
// differently than what is required for assemblePressureEq, this part needs to be assembled separately.
const int seg_outlet = this->segmentNumberToIndex(this->segmentSet()[seg].outletSegment());
const auto hydro_pressure_drop_seg = segments_.getHalfHydroPressureLoss(seg, seg);
const auto hydro_pressure_drop_outlet = segments_.getHalfHydroPressureLoss(seg, seg_outlet);
MultisegmentWellAssemble<FluidSystem,Indices,Scalar>(baseif_).
assembleHydroPressureLoss(seg, seg_outlet, hydro_pressure_drop_seg, hydro_pressure_drop_outlet, linSys_);
segments.pressure_drop_hydrostatic[seg] = hydro_pressure_drop_seg.value() + hydro_pressure_drop_outlet.value();
if (this->frictionalPressureLossConsidered()) {
const auto friction_pressure_drop = segments_.getFrictionPressureLoss(seg, false);
if (reverseFlow){
@@ -265,6 +255,16 @@ assembleDefaultPressureEq(const int seg,
if (this->accelerationalPressureLossConsidered()) {
handleAccelerationPressureLoss(seg, well_state);
}
// We approximate the hydrostatic pressure loss over the segment by using the mean of the densities
// at the current segment node and its outlet node. Since density derivatives are organized
// differently than what is required for assemblePressureEq, this part needs to be assembled separately.
const int seg_outlet = this->segmentNumberToIndex(this->segmentSet()[seg].outletSegment());
const auto hydro_pressure_drop_seg = segments_.getHalfHydroPressureLoss(seg, seg);
const auto hydro_pressure_drop_outlet = segments_.getHalfHydroPressureLoss(seg, seg_outlet);
MultisegmentWellAssemble<FluidSystem,Indices,Scalar>(baseif_).
assembleHydroPressureLoss(seg, seg_outlet, hydro_pressure_drop_seg, hydro_pressure_drop_outlet, linSys_);
segments.pressure_drop_hydrostatic[seg] = hydro_pressure_drop_seg.value() + hydro_pressure_drop_outlet.value();
}
template<typename FluidSystem, typename Indices, typename Scalar>