81 const SubControlVolume& scv,
82 const VolumeVariables& volVars)
const
84 NumEqVector storage(0.0);
87 for (
int phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx)
89 for (
int compIdx = 0; compIdx < numComponents; ++compIdx)
91 auto eqIdx = Indices::conti0EqIdx + compIdx;
92 storage[eqIdx] += volVars.porosity()
93 * volVars.saturation(phaseIdx)
94 * volVars.molarDensity(phaseIdx)
95 * volVars.moleFraction(phaseIdx, compIdx);
99 EnergyLocalResidual::fluidPhaseStorage(storage, scv, volVars, phaseIdx);
103 EnergyLocalResidual::solidPhaseStorage(storage, scv, volVars);
120 const Element& element,
121 const FVElementGeometry& fvGeometry,
122 const ElementVolumeVariables& elemVolVars,
123 const SubControlVolumeFace& scvf,
124 const ElementFluxVariablesCache& elemFluxVarsCache)
const
126 FluxVariables fluxVars;
127 fluxVars.init(problem, element, fvGeometry, elemVolVars, scvf, elemFluxVarsCache);
128 static constexpr auto referenceSystemFormulation = FluxVariables::MolecularDiffusionType::referenceSystemFormulation();
131 NumEqVector flux(0.0);
134 for (
int phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx)
136 for (
int compIdx = 0; compIdx < numComponents; ++compIdx)
138 auto upwindTerm = [phaseIdx, compIdx](
const VolumeVariables& volVars)
139 {
return volVars.molarDensity(phaseIdx)*volVars.moleFraction(phaseIdx, compIdx)*volVars.mobility(phaseIdx); };
142 auto eqIdx = Indices::conti0EqIdx + compIdx;
143 flux[eqIdx] += fluxVars.advectiveFlux(phaseIdx, upwindTerm);
147 EnergyLocalResidual::heatConvectionFlux(flux, fluxVars, phaseIdx);
151 EnergyLocalResidual::heatConductionFlux(flux, fluxVars);
154 const auto diffusionFluxesWPhase = fluxVars.molecularDiffusionFlux(wPhaseIdx);
155 Scalar jGW = diffusionFluxesWPhase[gCompIdx];
156 Scalar jNW = diffusionFluxesWPhase[nCompIdx];
157 Scalar jWW = -(jGW+jNW);
162 jGW /= FluidSystem::molarMass(gCompIdx);
163 jNW /= FluidSystem::molarMass(nCompIdx);
164 jWW /= FluidSystem::molarMass(wCompIdx);
167 const auto diffusionFluxesGPhase = fluxVars.molecularDiffusionFlux(gPhaseIdx);
168 Scalar jWG = diffusionFluxesGPhase[wCompIdx];
169 Scalar jNG = diffusionFluxesGPhase[nCompIdx];
170 Scalar jGG = -(jWG+jNG);
175 jWG /= FluidSystem::molarMass(wCompIdx);
176 jNG /= FluidSystem::molarMass(nCompIdx);
177 jGG /= FluidSystem::molarMass(gCompIdx);
181 const Scalar jWN = 0.0;
182 const Scalar jGN = 0.0;
183 const Scalar jNN = 0.0;
185 flux[contiWEqIdx] += jWW+jWG+jWN;
186 flux[contiNEqIdx] += jNW+jNG+jNN;
187 flux[contiGEqIdx] += jGW+jGG+jGN;