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

ebos: do the gravity correction the same way as opm-autodiff

Browse Source 
i.e., calculate the pressure which the exterior cell of a face would
exhibit at the depth of the interior one (instead of bringing both to
the depth of the face).

because the average fluid density is used, there should not be a
difference because of this.
This commit is contained in:
Andreas Lauser 2015-11-12 19:43:57 +01:00
parent 1c02969608
commit 4f4e25c25f
1 changed files with 10 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
applications/ebos/eclfluxmodule.hh
View File

@ -185,7 +185,7 @@ protected:
int downstreamIndex_(int phaseIdx) const int downstreamIndex_(int phaseIdx) const
{ {
assert(0 <= phaseIdx && phaseIdx < numPhases); assert(0 <= phaseIdx && phaseIdx < numPhases);
return (Toolbox::value(pressureDifferential_[phaseIdx]) >= 0)?interiorDofIdx_:exteriorDofIdx_; return (pressureDifferential_[phaseIdx] < 0)?exteriorDofIdx_:interiorDofIdx_;
} }
/*! /*!
@ -217,22 +217,21 @@ protected:
Scalar zIn = elemCtx.pos(interiorDofIdx_, timeIdx)[dimWorld - 1]; Scalar zIn = elemCtx.pos(interiorDofIdx_, timeIdx)[dimWorld - 1];
Scalar zEx = elemCtx.pos(exteriorDofIdx_, timeIdx)[dimWorld - 1]; Scalar zEx = elemCtx.pos(exteriorDofIdx_, timeIdx)[dimWorld - 1];
Scalar zFace = scvf.integrationPos()[dimWorld - 1];
Scalar distZIn = zIn - zFace; // the distances from the DOF's depths. (i.e., the additional depth of the
Scalar distZEx = zEx - zFace; // exterior DOF)
Scalar distZ = zIn - zEx;
for (int phaseIdx=0; phaseIdx < numPhases; phaseIdx++) { for (int phaseIdx=0; phaseIdx < numPhases; phaseIdx++) {
// do the gravity correction at the face's integration point // do the gravity correction: compute the hydrostatic pressure for the
// external at the depth of the internal one
const Evaluation& rhoIn = intQuantsIn.fluidState().density(phaseIdx); const Evaluation& rhoIn = intQuantsIn.fluidState().density(phaseIdx);
Scalar rhoEx = Toolbox::value(intQuantsEx.fluidState().density(phaseIdx)); Scalar rhoEx = Toolbox::value(intQuantsEx.fluidState().density(phaseIdx));
Evaluation rhoAvg = (rhoIn + rhoEx)/2; Evaluation rhoAvg = (rhoIn + rhoEx)/2;
Evaluation pressureInterior = intQuantsIn.fluidState().pressure(phaseIdx); const Evaluation& pressureInterior = intQuantsIn.fluidState().pressure(phaseIdx);
Evaluation pressureExterior = Toolbox::value(intQuantsEx.fluidState().pressure(phaseIdx)); Evaluation pressureExterior = Toolbox::value(intQuantsEx.fluidState().pressure(phaseIdx));
pressureExterior += rhoAvg*(distZ*g);
pressureInterior += - rhoAvg*(g*distZIn);
pressureExterior += - rhoAvg*(g*distZEx);
pressureDifferential_[phaseIdx] = pressureExterior - pressureInterior; pressureDifferential_[phaseIdx] = pressureExterior - pressureInterior;
@ -243,10 +242,10 @@ protected:
const auto& up = elemCtx.intensiveQuantities(upstreamIdx, timeIdx); const auto& up = elemCtx.intensiveQuantities(upstreamIdx, timeIdx);
if (upstreamIdx == interiorDofIdx_) if (upstreamIdx == interiorDofIdx_)
volumeFlux_[phaseIdx] = volumeFlux_[phaseIdx] =
pressureDifferential_[phaseIdx]*up.mobility(phaseIdx) * (- trans_/faceArea_); pressureDifferential_[phaseIdx]*up.mobility(phaseIdx)*(-trans_/faceArea_);
else else
volumeFlux_[phaseIdx] = volumeFlux_[phaseIdx] =
pressureDifferential_[phaseIdx]*(Toolbox::value(up.mobility(phaseIdx)) * (- trans_/faceArea_)); pressureDifferential_[phaseIdx]*(Toolbox::value(up.mobility(phaseIdx))*(-trans_/faceArea_));
} }
} }