56 using FluidSystem =
typename VolumeVariables::FluidSystem;
57 static_assert(FluidSystem::numPhases == 2,
"ThermalConductivitySomertonTwoP only works for two-phase fluid systems!");
58 static_assert((FluidSystem::isGas(0) && !FluidSystem::isGas(1)) || (!FluidSystem::isGas(0) && FluidSystem::isGas(1)),
59 "ThermalConductivitySomertonTwoP only works if one phase is gaseous and one is liquid!");
61 constexpr int liquidPhaseIdx = FluidSystem::isGas(0) ? 1 : 0;
62 constexpr int gasPhaseIdx = FluidSystem::isGas(0) ? 0 : 1;
64 const Scalar satLiquid = volVars.saturation(liquidPhaseIdx);
65 const Scalar lambdaLiquid = volVars.fluidThermalConductivity(liquidPhaseIdx);
66 const Scalar lambdaGas = volVars.fluidThermalConductivity(gasPhaseIdx);
67 const Scalar lambdaSolid = volVars.solidThermalConductivity();
68 const Scalar porosity = volVars.porosity();
70 return effectiveThermalConductivity_(satLiquid, lambdaLiquid, lambdaGas, lambdaSolid, porosity);