OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2026-08-13 14:35:18 -05:00
Code Issues Packages Projects Releases Wiki Activity

Merge pull request #1869 from GitPaean/adding_comment_getsegmentrateupwinding

Browse Source 
adding some comments about getSegmentRateUpwinding()
This commit is contained in:
Atgeirr Flø Rasmussen
2019-05-31 20:37:16 +02:00
committed by GitHub
parent 8796935f7c 8fd1a8b7da
commit 2d0cd0ecfe
1 changed files with 15 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files
opm/simulators/wells/MultisegmentWell_impl.hpp
+15 -4
View File
@@ -1391,7 +1391,14 @@ namespace Opm
const int comp_idx) const
{
const int seg_upwind = upwinding_segments_[seg];
return primary_variables_evaluation_[seg][GTotal] * volumeFractionScaled(seg_upwind, comp_idx);
// the result will contain the derivative with resepct to GTotal in segment seg,
// and the derivatives with respect to WFrac GFrac in segment seg_upwind.
// the derivative with respect to SPres should be zero.
const EvalWell segment_rate = primary_variables_evaluation_[seg][GTotal] * volumeFractionScaled(seg_upwind, comp_idx);
assert(segment_rate.derivative(SPres + numEq) == 0.);
return segment_rate;
}
@@ -1984,6 +1991,8 @@ namespace Opm
const EvalWell segment_rate = getSegmentRateUpwinding(seg, comp_idx);
const int seg_upwind = upwinding_segments_[seg];
// segment_rate contains the derivatives with respect to GTotal in seg,
// and WFrac and GFrac in seg_upwind
resWell_[seg][comp_idx] -= segment_rate.value();
duneD_[seg][seg][comp_idx][GTotal] -= segment_rate.derivative(GTotal + numEq);
duneD_[seg][seg_upwind][comp_idx][WFrac] -= segment_rate.derivative(WFrac + numEq);
@@ -1998,11 +2007,13 @@ namespace Opm
for (int comp_idx = 0; comp_idx < num_components_; ++comp_idx) {
const EvalWell inlet_rate = getSegmentRateUpwinding(inlet, comp_idx);
const int seg_upwind = upwinding_segments_[inlet];
const int inlet_upwind = upwinding_segments_[inlet];
// inlet_rate contains the derivatives with respect to GTotal in inlet,
// and WFrac and GFrac in inlet_upwind
resWell_[seg][comp_idx] += inlet_rate.value();
duneD_[seg][inlet][comp_idx][GTotal] += inlet_rate.derivative(GTotal + numEq);
duneD_[seg][seg_upwind][comp_idx][WFrac] += inlet_rate.derivative(WFrac + numEq);
duneD_[seg][seg_upwind][comp_idx][GFrac] += inlet_rate.derivative(GFrac + numEq);
duneD_[seg][inlet_upwind][comp_idx][WFrac] += inlet_rate.derivative(WFrac + numEq);
duneD_[seg][inlet_upwind][comp_idx][GFrac] += inlet_rate.derivative(GFrac + numEq);
// pressure derivative should be zero
}
}