88 void setFreeFlow(
const Context& context,
unsigned bfIdx,
unsigned timeIdx,
const FluidState& fluidState)
90 ExtensiveQuantities extQuants;
91 extQuants.updateBoundary(context, bfIdx, timeIdx, fluidState);
92 const auto& insideIntQuants = context.intensiveQuantities(bfIdx, timeIdx);
93 unsigned focusDofIdx = context.focusDofIndex();
94 unsigned interiorDofIdx = context.interiorScvIndex(bfIdx, timeIdx);
99 (*this) = Evaluation(0.0);
100 for (
unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
102 if (fluidState.pressure(phaseIdx) > insideIntQuants.fluidState().pressure(phaseIdx)) {
103 if (focusDofIdx == interiorDofIdx)
104 density = fluidState.density(phaseIdx);
106 density = Opm::getValue(fluidState.density(phaseIdx));
108 else if (focusDofIdx == interiorDofIdx)
109 density = insideIntQuants.fluidState().density(phaseIdx);
111 density = Opm::getValue(insideIntQuants.fluidState().density(phaseIdx));
113 for (
unsigned compIdx = 0; compIdx < numComponents; ++compIdx) {
115 if (fluidState.pressure(phaseIdx) > insideIntQuants.fluidState().pressure(phaseIdx)) {
116 if (focusDofIdx == interiorDofIdx)
117 molarity = fluidState.molarity(phaseIdx, compIdx);
119 molarity = Opm::getValue(fluidState.molarity(phaseIdx, compIdx));
121 else if (focusDofIdx == interiorDofIdx)
122 molarity = insideIntQuants.fluidState().molarity(phaseIdx, compIdx);
124 molarity = Opm::getValue(insideIntQuants.fluidState().molarity(phaseIdx, compIdx));
128 (*this)[conti0EqIdx + compIdx] += extQuants.volumeFlux(phaseIdx)*molarity;
132 Evaluation specificEnthalpy;
133 if (fluidState.pressure(phaseIdx) > insideIntQuants.fluidState().pressure(phaseIdx)) {
134 if (focusDofIdx == interiorDofIdx)
135 specificEnthalpy = fluidState.enthalpy(phaseIdx);
137 specificEnthalpy = Opm::getValue(fluidState.enthalpy(phaseIdx));
139 else if (focusDofIdx == interiorDofIdx)
140 specificEnthalpy = insideIntQuants.fluidState().enthalpy(phaseIdx);
142 specificEnthalpy = Opm::getValue(insideIntQuants.fluidState().enthalpy(phaseIdx));
144 Evaluation enthalpyRate = density*extQuants.volumeFlux(phaseIdx)*specificEnthalpy;
145 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]);