65 using FluidSystem =
typename VolumeVariables::FluidSystem;
66 static_assert(FluidSystem::numPhases == 2,
"ThermalConductivitySomerton only works for two-phase fluid systems!");
67 static_assert((FluidSystem::isGas(0) && !FluidSystem::isGas(1)) || (!FluidSystem::isGas(0) && FluidSystem::isGas(1)),
68 "ThermalConductivitySomerton only works if one phase is gaseous and one is liquid!");
70 constexpr int liquidPhaseIdx = FluidSystem::isGas(0) ? 1 : 0;
71 constexpr int gasPhaseIdx = FluidSystem::isGas(0) ? 0 : 1;
73 const Scalar satLiquid = volVars.saturation(liquidPhaseIdx);
74 const Scalar lambdaLiquid = volVars.fluidThermalConductivity(liquidPhaseIdx);
75 const Scalar lambdaGas = volVars.fluidThermalConductivity(gasPhaseIdx);
76 const Scalar lambdaSolid = volVars.solidThermalConductivity();
77 const Scalar porosity = volVars.porosity();
79 return effectiveThermalConductivity_(satLiquid, lambdaLiquid, lambdaGas, lambdaSolid, porosity);