56 using FluidSystem =
typename VolumeVariables::FluidSystem;
58 const Scalar sw = volVars.saturation(FluidSystem::wPhaseIdx);
59 const Scalar sn = volVars.saturation(FluidSystem::nPhaseIdx);
60 const Scalar lambdaW = volVars.fluidThermalConductivity(FluidSystem::wPhaseIdx);
61 const Scalar lambdaN = volVars.fluidThermalConductivity(FluidSystem::nPhaseIdx);
62 const Scalar lambdaG = volVars.fluidThermalConductivity(FluidSystem::gPhaseIdx);
63 const Scalar lambdaSolid = volVars.solidThermalConductivity();
64 const Scalar porosity = volVars.porosity();
87 const Scalar lambdaSolid,
88 const Scalar porosity)
93 const Scalar satW = max<Scalar>(0.0, sw);
94 const Scalar satN = max<Scalar>(0.0, sn);
97 const Scalar lSw = pow(lambdaSolid, (1.0 - porosity)) * pow(lambdaW, porosity);
98 const Scalar lSn = pow(lambdaSolid, (1.0 - porosity)) * pow(lambdaN, porosity);
99 const Scalar lSg = pow(lambdaSolid, (1.0 - porosity)) * pow(lambdaG, porosity);
100 const Scalar lambdaEff = lSg + sqrt(satW) * (lSw - lSg) + sqrt(satN) * (lSn -lSg);