OilfieldToolsNet/opm-simulators
4
0 Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2026-11-03 12:47:42 -06:00
Code Issues Packages Projects Releases Wiki Activity

richardslocalresidual: cosmetics

Browse Source
This commit is contained in:
Arne Morten Kvarving committed 2025-05-22 09:35:30 +02:00
1 parent 0041f89034
commit 0a4472114c
1 file changed
+12 -10
Download Patch File Download Diff File Expand all files Collapse all files
opm/models/richards/richardslocalresidual.hh
+12 -10
View File
@@ -55,7 +55,7 @@ class RichardsLocalResidual : public GetPropType<TypeTag, Properties::DiscLocalR
enum { contiEqIdx = Indices::contiEqIdx };
enum { liquidPhaseIdx = getPropValue<TypeTag, Properties::LiquidPhaseIndex>() };
enum { numEq = getPropValue<TypeTag, Properties::NumEq>() };
using Toolbox = Opm::MathToolbox<Evaluation>;
using Toolbox = MathToolbox<Evaluation>;
public:
/*!
@@ -71,9 +71,9 @@ public:
// partial time derivative of the wetting phase mass
storage[contiEqIdx] =
Toolbox::template decay<LhsEval>(intQuants.fluidState().density(liquidPhaseIdx))
*Toolbox::template decay<LhsEval>(intQuants.fluidState().saturation(liquidPhaseIdx))
*Toolbox::template decay<LhsEval>(intQuants.porosity());
Toolbox::template decay<LhsEval>(intQuants.fluidState().density(liquidPhaseIdx)) *
Toolbox::template decay<LhsEval>(intQuants.fluidState().saturation(liquidPhaseIdx)) *
Toolbox::template decay<LhsEval>(intQuants.porosity());
}
/*!
@@ -86,18 +86,20 @@ public:
{
const auto& extQuants = elemCtx.extensiveQuantities(scvfIdx, timeIdx);
unsigned focusDofIdx = elemCtx.focusDofIndex();
unsigned upIdx = static_cast<unsigned>(extQuants.upstreamIndex(liquidPhaseIdx));
const unsigned focusDofIdx = elemCtx.focusDofIndex();
const unsigned upIdx = static_cast<unsigned>(extQuants.upstreamIndex(liquidPhaseIdx));
const IntensiveQuantities& up = elemCtx.intensiveQuantities(upIdx, timeIdx);
// compute advective mass flux of the liquid phase. This is slightly hacky
// because it is specific to the element-centered finite volume method.
const Evaluation& rho = up.fluidState().density(liquidPhaseIdx);
if (focusDofIdx == upIdx)
flux[contiEqIdx] = extQuants.volumeFlux(liquidPhaseIdx)*rho;
else
flux[contiEqIdx] = extQuants.volumeFlux(liquidPhaseIdx)*Toolbox::value(rho);
if (focusDofIdx == upIdx) {
flux[contiEqIdx] = extQuants.volumeFlux(liquidPhaseIdx) * rho;
}
else {
flux[contiEqIdx] = extQuants.volumeFlux(liquidPhaseIdx) * Toolbox::value(rho);
}
}
/*!