53 constexpr int numFluidPhases = VolumeVariables::numFluidPhases();
56 Scalar lambdaFluid = 0.0;
57 for (
int phaseIdx = 0; phaseIdx < numFluidPhases; ++phaseIdx)
58 lambdaFluid += volVars.fluidThermalConductivity(phaseIdx)*volVars.saturation(phaseIdx);
60 const Scalar lambdaSolid = volVars.solidThermalConductivity();
61 const Scalar porosity = volVars.porosity();
63 return lambdaSolid*(1-porosity) + lambdaFluid*porosity;