68 using FluidSystem =
typename VolumeVariables::FluidSystem;
69 static_assert(FluidSystem::numPhases == 2,
"ThermalConductivitySomerton only works for two-phase fluid systems!");
72 const Scalar sw = volVars.saturation(volVars.wettingPhase());
73 const Scalar lambdaW = volVars.fluidThermalConductivity(volVars.wettingPhase());
74 const Scalar lambdaN = volVars.fluidThermalConductivity(1-volVars.wettingPhase());
75 const Scalar lambdaSolid = volVars.solidThermalConductivity();
76 const Scalar porosity = volVars.porosity();
77 const Scalar rhoSolid = volVars.solidDensity();
79 return effectiveThermalConductivity_(sw, lambdaW, lambdaN, lambdaSolid, porosity, rhoSolid);