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);