fix bug in boundary conditions for energy by adding the scaling factor

This commit is contained in:
Tor Harald Sandve 2022-10-21 15:34:05 +02:00
parent 11c9c34dcc
commit c23c64893e

View File

@ -161,7 +161,7 @@ public:
} }
Evaluation enthalpyRate = density*extQuants.volumeFlux(phaseIdx)*specificEnthalpy; Evaluation enthalpyRate = density*extQuants.volumeFlux(phaseIdx)*specificEnthalpy;
EnergyModule::addToEnthalpyRate(*this, enthalpyRate); EnergyModule::addToEnthalpyRate(*this, enthalpyRate*getPropValue<TypeTag, Properties::BlackOilEnergyScalingFactor>());
} }
} }
@ -189,7 +189,7 @@ public:
// heat conduction // heat conduction
if constexpr (enableEnergy) if constexpr (enableEnergy)
EnergyModule::addToEnthalpyRate(*this, extQuants.energyFlux()); EnergyModule::addToEnthalpyRate(*this, extQuants.energyFlux()*getPropValue<TypeTag, Properties::BlackOilEnergyScalingFactor>());
#ifndef NDEBUG #ifndef NDEBUG
for (unsigned i = 0; i < numEq; ++i) { for (unsigned i = 0; i < numEq; ++i) {