87 const FluidState& fluidState)
89 ExtensiveQuantities extQuants;
90 extQuants.updateBoundary(context, bfIdx, timeIdx, fluidState);
91 const auto& insideIntQuants = context.intensiveQuantities(bfIdx, timeIdx);
92 unsigned focusDofIdx = context.focusDofIndex();
93 unsigned interiorDofIdx = context.interiorScvIndex(bfIdx, timeIdx);
98 (*this) = Evaluation(0.0);
99 for (
unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
101 if (fluidState.pressure(phaseIdx) > insideIntQuants.fluidState().pressure(phaseIdx)) {
102 if (focusDofIdx == interiorDofIdx)
103 density = fluidState.density(phaseIdx);
105 density = Opm::getValue(fluidState.density(phaseIdx));
107 else if (focusDofIdx == interiorDofIdx)
108 density = insideIntQuants.fluidState().density(phaseIdx);
110 density = Opm::getValue(insideIntQuants.fluidState().density(phaseIdx));
112 for (
unsigned compIdx = 0; compIdx < numComponents; ++compIdx) {
114 if (fluidState.pressure(phaseIdx) > insideIntQuants.fluidState().pressure(phaseIdx)) {
115 if (focusDofIdx == interiorDofIdx)
116 molarity = fluidState.molarity(phaseIdx, compIdx);
118 molarity = Opm::getValue(fluidState.molarity(phaseIdx, compIdx));
120 else if (focusDofIdx == interiorDofIdx)
121 molarity = insideIntQuants.fluidState().molarity(phaseIdx, compIdx);
123 molarity = Opm::getValue(insideIntQuants.fluidState().molarity(phaseIdx, compIdx));
127 (*this)[conti0EqIdx + compIdx] += extQuants.volumeFlux(phaseIdx)*molarity;
131 Evaluation specificEnthalpy;
132 if (fluidState.pressure(phaseIdx) > insideIntQuants.fluidState().pressure(phaseIdx)) {
133 if (focusDofIdx == interiorDofIdx)
134 specificEnthalpy = fluidState.enthalpy(phaseIdx);
136 specificEnthalpy = Opm::getValue(fluidState.enthalpy(phaseIdx));
138 else if (focusDofIdx == interiorDofIdx)
139 specificEnthalpy = insideIntQuants.fluidState().enthalpy(phaseIdx);
141 specificEnthalpy = Opm::getValue(insideIntQuants.fluidState().enthalpy(phaseIdx));
143 Evaluation enthalpyRate = density*extQuants.volumeFlux(phaseIdx)*specificEnthalpy;
144 EnergyModule::addToEnthalpyRate(*
this, enthalpyRate);
149 EnergyModule::addToEnthalpyRate(*
this, EnergyModule::thermalConductionRate(extQuants));
152 for (
unsigned i = 0; i < numEq; ++i) {
153 Opm::Valgrind::CheckDefined((*
this)[i]);