[cleanup][nonequilibrium] cleanup and simplify

dumux/porousmediumflow/nonequilibrium/volumevariables.hh

@@ -270,10 +270,6 @@ class NonEquilibriumVolumeVariablesImplementation< Traits,
     static constexpr auto numEnergyEqFluid = ModelTraits::numEnergyEqFluid();
     static constexpr auto numEnergyEqSolid = ModelTraits::numEnergyEqSolid();
 
-    static constexpr auto phase0Idx = FS::phase0Idx;
-    static constexpr auto phase1Idx = FS::phase1Idx;
-    static constexpr auto sPhaseIdx = FS::numPhases;
-
     using DimLessNum = DimensionlessNumbers<Scalar>;
 
 public:

dumux/porousmediumflow/nonisothermal/volumevariables.hh

@@ -179,7 +179,10 @@ public:
             if (numEnergyEq == 2)
             {
                 const Scalar T = elemSol[scv.localDofIndex()][temperatureIdx];
-                fluidState.setTemperature(T);
+                for(int phaseIdx=0; phaseIdx < FluidSystem::numPhases; ++phaseIdx)
+                {
+                    fluidState.setTemperature(phaseIdx, T);
+                }
             }
             //this is for numEnergyEqFluid > 1
             else