97 const FluidState& fluidState)
99 ExtensiveQuantities extQuants;
100 extQuants.updateBoundary(context, bfIdx, timeIdx, fluidState);
101 const auto& insideIntQuants = context.intensiveQuantities(bfIdx, timeIdx);
102 unsigned focusDofIdx = context.focusDofIndex();
103 unsigned interiorDofIdx = context.interiorScvIndex(bfIdx, timeIdx);
109 for (
unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
110 if (!FluidSystem::phaseIsActive(phaseIdx)) {
113 const auto& pBoundary = fluidState.pressure(phaseIdx);
114 const Evaluation& pInside = insideIntQuants.fluidState().pressure(phaseIdx);
119 if (pBoundary < pInside)
121 LocalResidual::template evalPhaseFluxes_<Evaluation>(tmp,
123 insideIntQuants.pvtRegionIndex(),
125 insideIntQuants.fluidState());
126 else if (pBoundary > pInside) {
127 using RhsEval =
typename std::conditional<std::is_same<typename FluidState::Scalar, Evaluation>::value,
128 Evaluation, Scalar>::type;
130 LocalResidual::template evalPhaseFluxes_<RhsEval>(tmp,
132 insideIntQuants.pvtRegionIndex(),
137 for (
unsigned i = 0; i < tmp.size(); ++i)
138 (*
this)[i] += tmp[i];
141 if constexpr (enableEnergy) {
143 Evaluation specificEnthalpy;
144 if (pBoundary > pInside) {
145 if (focusDofIdx == interiorDofIdx) {
146 density = fluidState.density(phaseIdx);
147 specificEnthalpy = fluidState.enthalpy(phaseIdx);
150 density = getValue(fluidState.density(phaseIdx));
151 specificEnthalpy = getValue(fluidState.enthalpy(phaseIdx));
154 else if (focusDofIdx == interiorDofIdx) {
155 density = insideIntQuants.fluidState().density(phaseIdx);
156 specificEnthalpy = insideIntQuants.fluidState().enthalpy(phaseIdx);
159 density = getValue(insideIntQuants.fluidState().density(phaseIdx));
160 specificEnthalpy = getValue(insideIntQuants.fluidState().enthalpy(phaseIdx));
163 Evaluation enthalpyRate = density*extQuants.volumeFlux(phaseIdx)*specificEnthalpy;
168 if constexpr (enableSolvent) {
169 (*this)[Indices::contiSolventEqIdx] = extQuants.solventVolumeFlux();
170 if (blackoilConserveSurfaceVolume)
171 (*this)[Indices::contiSolventEqIdx] *= insideIntQuants.solventInverseFormationVolumeFactor();
173 (*
this)[Indices::contiSolventEqIdx] *= insideIntQuants.solventDensity();
177 if constexpr (enablePolymer) {
178 (*this)[Indices::contiPolymerEqIdx] = extQuants.volumeFlux(FluidSystem::waterPhaseIdx) * insideIntQuants.polymerConcentration();
181 if constexpr (enableMICP) {
182 (*this)[Indices::contiMicrobialEqIdx] = extQuants.volumeFlux(FluidSystem::waterPhaseIdx) * insideIntQuants.microbialConcentration();
183 (*this)[Indices::contiOxygenEqIdx] = extQuants.volumeFlux(FluidSystem::waterPhaseIdx) * insideIntQuants.oxygenConcentration();
184 (*this)[Indices::contiUreaEqIdx] = extQuants.volumeFlux(FluidSystem::waterPhaseIdx) * insideIntQuants.ureaConcentration();
188 LocalResidual::adaptMassConservationQuantities_(*
this, insideIntQuants.pvtRegionIndex());
191 if constexpr (enableEnergy)
195 for (
unsigned i = 0; i < numEq; ++i) {
196 Valgrind::CheckDefined((*
this)[i]);
198 Valgrind::CheckDefined(*
this);