62        using FluidSystem = 
typename VolumeVariables::FluidSystem;
 
   64        const Scalar sw = volVars.saturation(FluidSystem::wPhaseIdx);
 
   65        const Scalar sn = volVars.saturation(FluidSystem::nPhaseIdx);
 
   66        const Scalar lambdaW = volVars.fluidThermalConductivity(FluidSystem::wPhaseIdx);
 
   67        const Scalar lambdaN = volVars.fluidThermalConductivity(FluidSystem::nPhaseIdx);
 
   68        const Scalar lambdaG = volVars.fluidThermalConductivity(FluidSystem::gPhaseIdx);
 
   69        const Scalar lambdaSolid = volVars.solidThermalConductivity();
 
   70        const Scalar porosity = volVars.porosity();
 
 
   93                                               const Scalar lambdaSolid,
 
   94                                               const Scalar porosity)
 
   99        const Scalar satW = max<Scalar>(0.0, sw);
 
  100        const Scalar satN = max<Scalar>(0.0, sn);
 
  103        const Scalar lSw = pow(lambdaSolid, (1.0 - porosity)) * pow(lambdaW, porosity);
 
  104        const Scalar lSn = pow(lambdaSolid, (1.0 - porosity)) * pow(lambdaN, porosity);
 
  105        const Scalar lSg = pow(lambdaSolid, (1.0 - porosity)) * pow(lambdaG, porosity);
 
  106        const Scalar lambdaEff = lSg + sqrt(satW) * (lSw - lSg) + sqrt(satN) * (lSn -lSg);