89 const FluidState& fluidState)
91 ExtensiveQuantities extQuants;
92 extQuants.updateBoundary(context, bfIdx, timeIdx, fluidState);
93 const auto& insideIntQuants = context.intensiveQuantities(bfIdx, timeIdx);
94 unsigned focusDofIdx = context.focusDofIndex();
95 unsigned interiorDofIdx = context.interiorScvIndex(bfIdx, timeIdx);
100 (*this) = Evaluation(0.0);
101 for (
unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
103 if (fluidState.pressure(phaseIdx) > insideIntQuants.fluidState().pressure(phaseIdx)) {
104 if (focusDofIdx == interiorDofIdx)
105 density = fluidState.density(phaseIdx);
107 density = Opm::getValue(fluidState.density(phaseIdx));
109 else if (focusDofIdx == interiorDofIdx)
110 density = insideIntQuants.fluidState().density(phaseIdx);
112 density = Opm::getValue(insideIntQuants.fluidState().density(phaseIdx));
114 for (
unsigned compIdx = 0; compIdx < numComponents; ++compIdx) {
116 if (fluidState.pressure(phaseIdx) > insideIntQuants.fluidState().pressure(phaseIdx)) {
117 if (focusDofIdx == interiorDofIdx)
118 molarity = fluidState.molarity(phaseIdx, compIdx);
120 molarity = Opm::getValue(fluidState.molarity(phaseIdx, compIdx));
122 else if (focusDofIdx == interiorDofIdx)
123 molarity = insideIntQuants.fluidState().molarity(phaseIdx, compIdx);
125 molarity = Opm::getValue(insideIntQuants.fluidState().molarity(phaseIdx, compIdx));
129 (*this)[conti0EqIdx + compIdx] += extQuants.volumeFlux(phaseIdx)*molarity;
133 Evaluation specificEnthalpy;
134 if (fluidState.pressure(phaseIdx) > insideIntQuants.fluidState().pressure(phaseIdx)) {
135 if (focusDofIdx == interiorDofIdx)
136 specificEnthalpy = fluidState.enthalpy(phaseIdx);
138 specificEnthalpy = Opm::getValue(fluidState.enthalpy(phaseIdx));
140 else if (focusDofIdx == interiorDofIdx)
141 specificEnthalpy = insideIntQuants.fluidState().enthalpy(phaseIdx);
143 specificEnthalpy = Opm::getValue(insideIntQuants.fluidState().enthalpy(phaseIdx));
145 Evaluation enthalpyRate = density*extQuants.volumeFlux(phaseIdx)*specificEnthalpy;
146 EnergyModule::addToEnthalpyRate(*
this, enthalpyRate);
151 EnergyModule::addToEnthalpyRate(*
this, EnergyModule::thermalConductionRate(extQuants));
154 for (
unsigned i = 0; i < numEq; ++i) {
155 Opm::Valgrind::CheckDefined((*
this)[i]);