36 constexpr int numFluidPhases = VolumeVariables::numFluidPhases();
39 Scalar lambdaFluid = 0.0;
40 for (
int phaseIdx = 0; phaseIdx < numFluidPhases; ++phaseIdx)
41 lambdaFluid += volVars.fluidThermalConductivity(phaseIdx)*volVars.saturation(phaseIdx);
43 const Scalar lambdaSolid = volVars.solidThermalConductivity();
44 const Scalar porosity = volVars.porosity();
46 return lambdaSolid*(1-porosity) + lambdaFluid*porosity;