55 using FluidSystem =
typename VolumeVariables::FluidSystem;
56 static_assert(FluidSystem::numPhases == 2,
"ThermalConductivitySomerton only works for two-phase fluid systems!");
59 const Scalar sw = volVars.saturation(volVars.wettingPhase());
60 const Scalar lambdaW = volVars.fluidThermalConductivity(volVars.wettingPhase());
61 const Scalar lambdaN = volVars.fluidThermalConductivity(1-volVars.wettingPhase());
62 const Scalar lambdaSolid = volVars.solidThermalConductivity();
63 const Scalar porosity = volVars.porosity();
64 const Scalar rhoSolid = volVars.solidDensity();
66 return effectiveThermalConductivity_(sw, lambdaW, lambdaN, lambdaSolid, porosity, rhoSolid);