OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-28 18:21:00 -06:00
Code Issues Packages Projects Releases Wiki Activity

use diffusivity

Browse Source
This commit is contained in:
Tor Harald Sandve 2021-01-14 15:32:41 +01:00
parent eadcf4f99d
commit 87abc2c36b
2 changed files with 7 additions and 70 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

75
opm/models/blackoil/blackoildiffusionmodule.hh
View File

@ -388,27 +388,13 @@ protected:
*/
void update_(const ElementContext& elemCtx, unsigned faceIdx, unsigned timeIdx)
{
// This uses a simple finite difference approach and is only correct for a Cartesian grids.
// TODO compute diffusive "transmissibilities" by re-using the logic in EclTransmissibility
const auto& stencil = elemCtx.stencil(timeIdx);
const auto& face = stencil.interiorFace(faceIdx);
const auto& extQuants = elemCtx.extensiveQuantities(faceIdx, timeIdx);
const auto& intQuantsInside = elemCtx.intensiveQuantities(extQuants.interiorIndex(), timeIdx);
const auto& intQuantsOutside = elemCtx.intensiveQuantities(extQuants.exteriorIndex(), timeIdx);
// distance between the centers
const auto& insideScv = stencil.subControlVolume(face.interiorIndex());
const auto& outsideScv = stencil.subControlVolume(face.exteriorIndex());
DimVector distVec = outsideScv.globalPos();
distVec -= insideScv.globalPos();
Scalar dist = distVec.two_norm();
// if the following assertation triggers, the center of the
// center of the interior SCV was not inside the element!
assert(dist > 0);
Scalar diffusivity = elemCtx.problem().diffusivity(elemCtx, face.interiorIndex(), face.exteriorIndex());
for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
if (!FluidSystem::phaseIsActive(phaseIdx)) {
continue;
@ -422,7 +408,7 @@ protected:
(intQuantsOutside.fluidState().moleFraction(phaseIdx, compIdx)
-
intQuantsInside.fluidState().moleFraction(phaseIdx, compIdx))
/ dist;
* diffusivity;
Opm::Valgrind::CheckDefined(moleFractionGradientNormal_[phaseIdx][compIdx]);
// use the arithmetic average for the effective
@ -439,59 +425,12 @@ protected:
}
template <class Context, class FluidState>
void updateBoundary_(const Context& context,
unsigned bfIdx,
unsigned timeIdx,
const FluidState& fluidState)
void updateBoundary_(const Context& context OPM_UNUSED,
unsigned bfIdx OPM_UNUSED,
unsigned timeIdx OPM_UNUSED,
const FluidState& fluidState OPM_UNUSED)
{
// This uses a simple finite difference approach and is only correct for a Cartesian grids.
// TODO compute diffusive "transmissibilities" by re-using the logic in EclTransmissibility
const auto& stencil = context.stencil(timeIdx);
const auto& face = stencil.boundaryFace(bfIdx);
const auto& elemCtx = context.elementContext();
unsigned insideScvIdx = face.interiorIndex();
const auto& insideScv = stencil.subControlVolume(insideScvIdx);
const auto& intQuantsInside = elemCtx.intensiveQuantities(insideScvIdx, timeIdx);
const auto& fluidStateInside = intQuantsInside.fluidState();
// distance between the center of the SCV and center of the boundary face
DimVector distVec = face.integrationPos();
distVec -= context.element().geometry().global(insideScv.localGeometry().center());
Scalar dist = distVec.two_norm();
// if the following assertation triggers, the center of the
// center of the interior SCV was not inside the element!
assert(dist > 0);
for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
if (!FluidSystem::phaseIsActive(phaseIdx)) {
continue;
}
// no diffusion in water for blackoil models
if (FluidSystem::waterPhaseIdx == phaseIdx) {
continue;
}
for (unsigned compIdx = 0; compIdx < numComponents; ++compIdx) {
// calculate mole fraction gradient using two-point
// gradients
moleFractionGradientNormal_[phaseIdx][compIdx] =
(fluidState.moleFraction(phaseIdx, compIdx)
-
fluidStateInside.moleFraction(phaseIdx, compIdx))
/ dist;
Opm::Valgrind::CheckDefined(moleFractionGradientNormal_[phaseIdx][compIdx]);
// use effective diffusion coefficients of the interior finite
// volume.
effectiveDiffusionCoefficient_[phaseIdx][compIdx] =
intQuantsInside.effectiveDiffusionCoefficient(phaseIdx, compIdx);
Opm::Valgrind::CheckDefined(effectiveDiffusionCoefficient_[phaseIdx][compIdx]);
}
}
throw std::runtime_error("Not implemented: Diffusion across boundary not implemented for blackoil");
}
public:

2
opm/models/blackoil/blackoilextensivequantities.hh
View File

@ -94,8 +94,6 @@ public:
MultiPhaseParent::updateBoundary(ctx, bfIdx, timeIdx, fluidState);
asImp_().updateEnergyBoundary(ctx, bfIdx, timeIdx, fluidState);
DiffusionExtensiveQuantities::updateBoundary_(ctx, bfIdx, timeIdx, fluidState);
}
protected: