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

don't upwind densitiy

Browse Source
This commit is contained in:
Stein Krogstad 2023-05-09 23:24:25 +02:00
parent e5ac9f9dec
commit 8ee3493aed
2 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -226,7 +226,7 @@ assembleHydroPressureLoss(const int seg,
MultisegmentWellEquationAccess<Scalar,numWellEq,Indices::numEq> eqns(eqns1);
eqns.residual()[seg][SPres] -= hydro_pressure_drop.value();
for (int pv_idx = 0; pv_idx < numWellEq; ++pv_idx) {
eqns.D()[seg][seg_upwind][SPres][pv_idx] -= hydro_pressure_drop.derivative(pv_idx + Indices::numEq);
eqns.D()[seg][seg][SPres][pv_idx] -= hydro_pressure_drop.derivative(pv_idx + Indices::numEq);
}
}

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

@ -323,8 +323,8 @@ typename MultisegmentWellSegments<FluidSystem,Indices,Scalar>::EvalWell
MultisegmentWellSegments<FluidSystem,Indices,Scalar>::
getHydroPressureLoss(const int seg) const
{
const int seg_upwind = upwinding_segments_[seg];
return densities_[seg_upwind] * well_.gravity() * depth_diffs_[seg];
//const int seg_upwind = upwinding_segments_[seg];
return densities_[seg] * well_.gravity() * depth_diffs_[seg];
}