[feature][fluxvariables] allow heatconduction calculation for each phase

dumux/porousmediumflow/fluxvariables.hh

@@ -61,6 +61,7 @@ class PorousMediumFluxVariables : public FluxVariablesBase<TypeTag>
     static constexpr bool enableAdvection = GET_PROP_VALUE(TypeTag, EnableAdvection);
     static constexpr bool enableMolecularDiffusion = GET_PROP_VALUE(TypeTag, EnableMolecularDiffusion);
     static constexpr bool enableEnergyBalance = GET_PROP_VALUE(TypeTag, EnableEnergyBalance);
+    static constexpr bool enableThermalNonEquilibrium = GET_PROP_VALUE(TypeTag, EnableThermalNonEquilibrium);
 
 public:
 
@@ -115,6 +116,7 @@ public:
         }
     }
 
+    template <bool enable = !enableThermalNonEquilibrium, typename std::enable_if_t<enable, int> = 0>
     Scalar heatConductionFlux() const
     {
         if (enableEnergyBalance)
@@ -132,6 +134,19 @@ public:
         }
     }
 
+    template <bool enable = enableThermalNonEquilibrium, typename std::enable_if_t<enable, int> = 0>
+    Scalar heatConductionFlux(const int phaseIdx) const
+    {
+            return HeatConductionType::flux(this->problem(),
+                                            this->element(),
+                                            this->fvGeometry(),
+                                            this->elemVolVars(),
+                                            this->scvFace(),
+                                            phaseIdx,
+                                            this->elemFluxVarsCache());
+
+    }
+
 private:
     //! simple caching if advection flux is used twice with different upwind function
     mutable std::bitset<numPhases> advFluxIsCached_;