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

default value of option for extra derivatives

Browse Source
This commit is contained in:
Stein Krogstad
2023-10-19 15:10:14 +02:00
parent dd6d195a6b
commit 57c960d7c5
3 changed files with 24 additions and 22 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@@ -234,10 +234,10 @@ assembleDefaultPressureEq(const int seg,
auto& segments = ws.segments;
if (this->frictionalPressureLossConsidered()) {
const auto friction_pressure_drop = segments_.getFrictionPressureLoss(seg, false);
const auto friction_pressure_drop = segments_.getFrictionPressureLoss(seg);
if (reverseFlow){
// call function once again to obtain/assemble remaining derivatives
extra_derivatives = -segments_.getFrictionPressureLoss(seg, true);
extra_derivatives = -segments_.getFrictionPressureLoss(seg, /*extra_reverse_flow_derivatives*/ true);
MultisegmentWellAssemble<FluidSystem,Indices,Scalar>(baseif_).
assemblePressureEqExtraDerivatives(seg, seg_upwind, extra_derivatives, linSys_);
}
@@ -292,21 +292,21 @@ assembleICDPressureEq(const int seg,
EvalWell extra_derivatives;
switch(this->segmentSet()[seg].segmentType()) {
case Segment::SegmentType::SICD :
icd_pressure_drop = segments_.pressureDropSpiralICD(seg, /*extra derivatives*/false);
icd_pressure_drop = segments_.pressureDropSpiralICD(seg);
if (reverseFlow){
extra_derivatives = -segments_.pressureDropSpiralICD(seg, /*extra derivatives*/true);
extra_derivatives = -segments_.pressureDropSpiralICD(seg, /*extra_reverse_flow_derivatives*/ true);
}
break;
case Segment::SegmentType::AICD :
icd_pressure_drop = segments_.pressureDropAutoICD(seg, unit_system, /*extra derivatives*/false);
icd_pressure_drop = segments_.pressureDropAutoICD(seg, unit_system);
if (reverseFlow){
extra_derivatives = -segments_.pressureDropAutoICD(seg, unit_system, /*extra derivatives*/true);
extra_derivatives = -segments_.pressureDropAutoICD(seg, unit_system, /*extra_reverse_flow_derivatives*/ true);
}
break;
case Segment::SegmentType::VALVE :
icd_pressure_drop = segments_.pressureDropValve(seg, /*extra derivatives*/false);
icd_pressure_drop = segments_.pressureDropValve(seg);
if (reverseFlow){
extra_derivatives = -segments_.pressureDropValve(seg, /*extra derivatives*/true);
extra_derivatives = -segments_.pressureDropValve(seg, /*extra_reverse_flow_derivatives*/ true);
}
break;
default: {

21
opm/simulators/wells/MultisegmentWellSegments.cpp
View File

@@ -498,7 +498,8 @@ getSurfaceVolume(const EvalWell& temperature,
template<class FluidSystem, class Indices, class Scalar>
typename MultisegmentWellSegments<FluidSystem,Indices,Scalar>::EvalWell
MultisegmentWellSegments<FluidSystem,Indices,Scalar>::
getFrictionPressureLoss(const int seg, const bool return_extra_derivatives) const
getFrictionPressureLoss(const int seg,
const bool extra_reverse_flow_derivatives /*false*/) const
{
EvalWell mass_rate = mass_rates_[seg];
const int seg_upwind = upwinding_segments_[seg];
@@ -511,7 +512,7 @@ getFrictionPressureLoss(const int seg, const bool return_extra_derivatives) cons
// at segment node while fraction derivatives are given at upwind node.
if (seg != seg_upwind) {
if (!return_extra_derivatives){
if (!extra_reverse_flow_derivatives){
constexpr int WQTotal = Indices::numEq + PrimaryVariables::WQTotal;
constexpr int SPres = Indices::numEq + PrimaryVariables::SPres;
density.setDerivative(WQTotal, 0.0);
@@ -550,7 +551,7 @@ template<class FluidSystem, class Indices, class Scalar>
typename MultisegmentWellSegments<FluidSystem,Indices,Scalar>::EvalWell
MultisegmentWellSegments<FluidSystem,Indices,Scalar>::
pressureDropSpiralICD(const int seg,
const bool return_extra_derivatives) const
const bool extra_reverse_flow_derivatives /*false*/) const
{
const auto& segment_set = well_.wellEcl().getSegments();
const SICD& sicd = segment_set[seg].spiralICD();
@@ -593,14 +594,14 @@ pressureDropSpiralICD(const int seg,
if (seg != seg_upwind) {
constexpr int nvar = FluidSystem::numPhases + 1;
std::vector<bool> zero_mask(nvar, false);
if (!return_extra_derivatives){
if (!extra_reverse_flow_derivatives){
zero_mask[PrimaryVariables::WQTotal] = true;
zero_mask[PrimaryVariables::SPres] = true;
} else {
if (PrimaryVariables::has_water){
if constexpr (PrimaryVariables::has_water){
zero_mask[PrimaryVariables::WFrac] = true;
}
if (PrimaryVariables::has_gas){
if constexpr (PrimaryVariables::has_gas){
zero_mask[PrimaryVariables::GFrac] = true;
}
// mass_rate has no extra derivatives (they are organized as in equations)
@@ -654,7 +655,7 @@ typename MultisegmentWellSegments<FluidSystem,Indices,Scalar>::EvalWell
MultisegmentWellSegments<FluidSystem,Indices,Scalar>::
pressureDropAutoICD(const int seg,
const UnitSystem& unit_system,
const bool return_extra_derivatives) const
const bool extra_reverse_flow_derivatives /*false*/) const
{
const auto& segment_set = well_.wellEcl().getSegments();
const AutoICD& aicd = segment_set[seg].autoICD();
@@ -704,7 +705,7 @@ pressureDropAutoICD(const int seg,
if (seg != seg_upwind) {
constexpr int nvar = FluidSystem::numPhases + 1;
std::vector<bool> zero_mask(nvar, false);
if (!return_extra_derivatives){
if (!extra_reverse_flow_derivatives){
zero_mask[PrimaryVariables::WQTotal] = true;
zero_mask[PrimaryVariables::SPres] = true;
} else {
@@ -762,7 +763,7 @@ template<class FluidSystem, class Indices, class Scalar>
typename MultisegmentWellSegments<FluidSystem,Indices,Scalar>::EvalWell
MultisegmentWellSegments<FluidSystem,Indices,Scalar>::
pressureDropValve(const int seg,
const bool return_extra_derivatives) const
const bool extra_reverse_flow_derivatives /*false*/) const
{
const auto& segment_set = well_.wellEcl().getSegments();
const Valve& valve = segment_set[seg].valve();
@@ -777,7 +778,7 @@ pressureDropValve(const int seg,
// For reference: the pressure equation assumes pressure/flow derivatives are given
// at segment node while fraction derivatives are given at upwind node.
if (seg != seg_upwind) {
if (!return_extra_derivatives){
if (!extra_reverse_flow_derivatives){
constexpr int WQTotal = Indices::numEq + PrimaryVariables::WQTotal;
constexpr int SPres = Indices::numEq + PrimaryVariables::SPres;
density.setDerivative(WQTotal, 0.0);

9
opm/simulators/wells/MultisegmentWellSegments.hpp
View File

@@ -71,20 +71,21 @@ public:
const int pvt_region_index,
const int seg_idx) const;
EvalWell getFrictionPressureLoss(const int seg, const bool return_upwind_derivatives) const;
EvalWell getFrictionPressureLoss(const int seg,
const bool extra_reverse_flow_derivatives = false) const;
// pressure drop for Spiral ICD segment (WSEGSICD)
EvalWell pressureDropSpiralICD(const int seg,
const bool return_extra_derivatives) const;
const bool extra_reverse_flow_derivatives = false) const;
// pressure drop for Autonomous ICD segment (WSEGAICD)
EvalWell pressureDropAutoICD(const int seg,
const UnitSystem& unit_system,
const bool return_extra_derivatives) const;
const bool extra_reverse_flow_derivatives = false) const;
// pressure drop for sub-critical valve (WSEGVALV)
EvalWell pressureDropValve(const int seg,
const bool return_extra_derivatives) const;
const bool extra_reverse_flow_derivatives = false) const;
// pressure loss due to acceleration
EvalWell accelerationPressureLoss(const int seg) const;