61 using FluidSystem =
typename VolumeVariables::FluidSystem;
62 static_assert(FluidSystem::numPhases == 2,
"ThermalConductivitySomerton only works for two-phase fluid systems!");
65 const Scalar sw = volVars.saturation(volVars.wettingPhase());
66 const Scalar lambdaW = volVars.fluidThermalConductivity(volVars.wettingPhase());
67 const Scalar lambdaN = volVars.fluidThermalConductivity(1-volVars.wettingPhase());
68 const Scalar lambdaSolid = volVars.solidThermalConductivity();
69 const Scalar porosity = volVars.porosity();
70 const Scalar rhoSolid = volVars.solidDensity();
72 return effectiveThermalConductivity_(sw, lambdaW, lambdaN, lambdaSolid, porosity, rhoSolid);