107 static void solve(FluidState& fluidState,
108 typename FluidSystem::template ParameterCache<typename FluidState::Scalar>& paramCache,
109 unsigned refPhaseIdx,
115 paramCache.updatePhase(fluidState, refPhaseIdx);
116 fluidState.setDensity(refPhaseIdx,
117 FluidSystem::density(fluidState,
122 fluidState.setEnthalpy(refPhaseIdx,
123 FluidSystem::enthalpy(fluidState,
128 fluidState.setViscosity(refPhaseIdx,
129 FluidSystem::viscosity(fluidState,
134 for (
unsigned compIdx = 0; compIdx < numComponents; ++compIdx) {
135 fluidState.setFugacityCoefficient(refPhaseIdx,
137 FluidSystem::fugacityCoefficient(fluidState,
144 for (
unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
145 if (phaseIdx == refPhaseIdx)
148 ComponentVector fugVec;
149 for (
unsigned compIdx = 0; compIdx < numComponents; ++compIdx) {
150 const auto& fug = fluidState.fugacity(refPhaseIdx, compIdx);
151 fugVec[compIdx] = decay<Evaluation>(fug);
157 fluidState.setViscosity(phaseIdx,
158 FluidSystem::viscosity(fluidState,
163 fluidState.setEnthalpy(phaseIdx,
164 FluidSystem::enthalpy(fluidState,