100 void update(
const ElementContext& elemCtx,
104 ParentType::update(elemCtx, dofIdx, timeIdx);
105 ParentType::checkDefined();
107 typename FluidSystem::template ParameterCache<Evaluation> paramCache;
108 const auto& priVars = elemCtx.primaryVars(dofIdx, timeIdx);
111 Evaluation sumSat = 0;
112 for (
unsigned phaseIdx = 0; phaseIdx < numPhases - 1; ++phaseIdx) {
113 const Evaluation& val = priVars.makeEvaluation(saturation0Idx + phaseIdx, timeIdx);
114 fluidState_.setSaturation(phaseIdx, val);
117 fluidState_.setSaturation(numPhases - 1, 1.0 - sumSat);
118 Opm::Valgrind::CheckDefined(sumSat);
121 EnergyIntensiveQuantities::updateTemperatures_(fluidState_, elemCtx, dofIdx, timeIdx);
124 const auto& problem = elemCtx.problem();
125 const MaterialLawParams& materialParams =
126 problem.materialLawParams(elemCtx, dofIdx, timeIdx);
128 Evaluation capPress[numPhases];
129 MaterialLaw::capillaryPressures(capPress, materialParams, fluidState_);
131 const Evaluation& pressure0 = priVars.makeEvaluation(pressure0Idx, timeIdx);
132 for (
unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx)
133 fluidState_.setPressure(phaseIdx, pressure0 + (capPress[phaseIdx] - capPress[0]));
137 for (
unsigned compIdx = 0; compIdx < numComponents; ++compIdx)
138 fug[compIdx] = priVars.makeEvaluation(fugacity0Idx + compIdx, timeIdx);
141 const auto *hint = elemCtx.thermodynamicHint(dofIdx, timeIdx);
142 for (
unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
145 for (
unsigned compIdx = 0; compIdx < numComponents; ++compIdx) {
147 const Evaluation& moleFracIJ = hint->fluidState().moleFraction(phaseIdx, compIdx);
148 fluidState_.setMoleFraction(phaseIdx, compIdx, moleFracIJ);
152 CompositionFromFugacitiesSolver::guessInitial(fluidState_, phaseIdx, fug);
156 CompositionFromFugacitiesSolver::solve(fluidState_, paramCache, phaseIdx, fug);
160 porosity_ = problem.porosity(elemCtx, dofIdx, timeIdx);
161 Opm::Valgrind::CheckDefined(porosity_);
164 MaterialLaw::relativePermeabilities(relativePermeability_, materialParams, fluidState_);
167 for (
unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
169 const Evaluation& mu = FluidSystem::viscosity(fluidState_, paramCache, phaseIdx);
170 fluidState_.setViscosity(phaseIdx, mu);
172 mobility_[phaseIdx] = relativePermeability_[phaseIdx]/mu;
176 intrinsicPerm_ = problem.intrinsicPermeability(elemCtx, dofIdx, timeIdx);
179 FluxIntensiveQuantities::update_(elemCtx, dofIdx, timeIdx);
182 EnergyIntensiveQuantities::update_(fluidState_, paramCache, elemCtx, dofIdx, timeIdx);
185 DiffusionIntensiveQuantities::update_(fluidState_, paramCache, elemCtx, dofIdx, timeIdx);