70 using FluidSystem =
typename VolumeVariables::FluidSystem;
72 const Scalar sw = volVars.saturation(FluidSystem::wPhaseIdx);
73 const Scalar sn = volVars.saturation(FluidSystem::nPhaseIdx);
74 const Scalar lambdaW = volVars.fluidThermalConductivity(FluidSystem::wPhaseIdx);
75 const Scalar lambdaN = volVars.fluidThermalConductivity(FluidSystem::nPhaseIdx);
76 const Scalar lambdaG = volVars.fluidThermalConductivity(FluidSystem::gPhaseIdx);
77 const Scalar lambdaSolid = volVars.solidThermalConductivity();
78 const Scalar porosity = volVars.porosity();
100 const Scalar lambdaG,
101 const Scalar lambdaSolid,
102 const Scalar porosity)
107 const Scalar satW = max<Scalar>(0.0, sw);
108 const Scalar satN = max<Scalar>(0.0, sn);
111 const Scalar lSw = pow(lambdaSolid, (1.0 - porosity)) * pow(lambdaW, porosity);
112 const Scalar lSn = pow(lambdaSolid, (1.0 - porosity)) * pow(lambdaN, porosity);
113 const Scalar lSg = pow(lambdaSolid, (1.0 - porosity)) * pow(lambdaG, porosity);
114 const Scalar lambdaEff = lSg + sqrt(satW) * (lSw - lSg) + sqrt(satN) * (lSn -lSg);
static Scalar effectiveThermalConductivity(const Scalar sw, const Scalar sn, const Scalar lambdaW, const Scalar lambdaN, const Scalar lambdaG, const Scalar lambdaSolid, const Scalar porosity)
effective thermal conductivity after Somerton (1974)
Definition thermalconductivitysomerton3p.hh:96