278 const Scalar faceArea,
279 const ConvectiveMixingModuleParam& info)
281 if (info.active_.empty()) {
289 const auto& liquidPhaseIdx = (FluidSystem::phaseIsActive(FluidSystem::waterPhaseIdx)) ?
290 FluidSystem::waterPhaseIdx :
291 FluidSystem::oilPhaseIdx;
292 const Evaluation
SoMax = 0.0;
295 const auto&
t_in =
intQuantsIn.fluidState().temperature(liquidPhaseIdx);
303 const auto bLiquidSatIn = (FluidSystem::phaseIsActive(FluidSystem::waterPhaseIdx)) ?
307 const auto&
densityLiquidIn = (FluidSystem::phaseIsActive(FluidSystem::waterPhaseIdx)) ?
308 FluidSystem::waterPvt().waterReferenceDensity(
intQuantsIn.pvtRegionIndex()) :
309 FluidSystem::oilPvt().oilReferenceDensity(
intQuantsIn.pvtRegionIndex());
323 const auto bLiquidSatEx = (FluidSystem::phaseIsActive(FluidSystem::waterPhaseIdx)) ?
328 const auto&
densityLiquidEx = (FluidSystem::phaseIsActive(FluidSystem::waterPhaseIdx)) ?
329 FluidSystem::waterPvt().waterReferenceDensity(
intQuantsEx.pvtRegionIndex()) :
330 FluidSystem::oilPvt().oilReferenceDensity(
intQuantsEx.pvtRegionIndex());
354 const auto&
Rsup = (FluidSystem::phaseIsActive(FluidSystem::waterPhaseIdx)) ?
355 up.fluidState().Rsw() :
356 up.fluidState().Rs();
358 const Evaluation&
transMult =
up.rockCompTransMultiplier();
359 const auto&
invB =
up.fluidState().invB(liquidPhaseIdx);
360 const auto&
visc =
up.fluidState().viscosity(liquidPhaseIdx);
366 unsigned activeGasCompIdx = Indices::canonicalToActiveComponentIndex(FluidSystem::gasCompIdx);
372 if constexpr (enableEnergy) {
373 const auto&
h =
up.fluidState().enthalpy(liquidPhaseIdx) * FluidSystem::referenceDensity(FluidSystem::gasPhaseIdx,
up.pvtRegionIndex());