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