removal of incomplete interface for basing mass transfer on chemical...

removal of incomplete interface for basing mass transfer on chemical potentials in stead of mole fractions. 
Future work ...

reviwed by Christoph

git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@12744 2fb0f335-1f38-0410-981e-8018bf24f1b0

dumux/implicit/mpnc/mass/mpncvolumevariablesmasskinetic.hh

@@ -36,6 +36,8 @@ namespace Dumux
 /*!
  * \brief The compositional part of the volume variables if chemical
  *        equilibrium is _not_ assumed
+ *
+ *        The interface for basing mass transfer on chemical potentials was present in revision 12743
  */
 template <class TypeTag>
 class MPNCVolumeVariablesMass<TypeTag, /*bool enableKinetic=*/true>
@@ -130,24 +132,6 @@ public:
                 }
             }
 
-            // Setting the equilibrium chemical potential (in a kinetic model not necessarily the same as the actual chemical potential)
-            for(int smallLoopPhaseIdx=0; smallLoopPhaseIdx<numPhases; ++smallLoopPhaseIdx){
-                for (int compIdx=0; compIdx< numComponents; ++ compIdx){
-                    chemicalPotentialEquil_[smallLoopPhaseIdx][compIdx] = FluidSystem::chemicalPotential(equilFluidState,
-                                                                                                         smallLoopPhaseIdx,
-                                                                                                         compIdx);
-                }
-            }
-
-            // Setting the actual chemical potential (in a kinetic model not necessarily the same as the equilibrium chemical potential)
-            for(int smallLoopPhaseIdx=0; smallLoopPhaseIdx<numPhases; ++smallLoopPhaseIdx){
-                for (int compIdx=0; compIdx< numComponents; ++ compIdx){
-                    chemicalPotential_[smallLoopPhaseIdx][compIdx] = FluidSystem::chemicalPotential(actualFluidState,
-                                                                                                    smallLoopPhaseIdx,
-                                                                                                    compIdx);
-                }
-            }
-
             // compute densities of all phases
             for(int smallLoopPhaseIdx=0; smallLoopPhaseIdx<numPhases; ++smallLoopPhaseIdx){
                 const Scalar rho = FluidSystem::density(actualFluidState, paramCache, smallLoopPhaseIdx);
@@ -170,28 +154,6 @@ public:
         return xEquil_[phaseIdx][compIdx] ;
     }
 
-    /*!
-     * \brief The chemical potential we would have in the case of chemical equilibrium / on the interface
-     *
-     *     \param phaseIdx The index of the fluid phase
-     *     \param compIdx The local index of the component
-     */
-    const Scalar chemicalPotential(const unsigned int phaseIdx, const unsigned int compIdx) const
-    {
-        return chemicalPotential_[phaseIdx][compIdx] ;
-    }
-
-    /*!
-     * \brief The actual chemical potential we currently have. In the case of non-equilibrium this is not
-     *        necessarily equilibrium.
-     *
-     *     \param phaseIdx The index of the fluid phase
-     *     \param compIdx The local index of the component
-     */
-    const Scalar chemicalPotentialEquil(const unsigned int phaseIdx, const unsigned int compIdx) const
-    {
-        return chemicalPotentialEquil_[phaseIdx][compIdx] ;
-    }
 
     /*!
      * \brief If running in valgrind this makes sure that all
@@ -201,15 +163,11 @@ public:
     {
 #if HAVE_VALGRIND && !defined NDEBUG
         Valgrind::CheckDefined(xEquil_);
-        Valgrind::CheckDefined(chemicalPotentialEquil_);
-        Valgrind::CheckDefined(chemicalPotential_);
 #endif
     }
 
 private:
     Scalar xEquil_[numPhases][numComponents];
-    Scalar chemicalPotentialEquil_[numPhases][numComponents];
-    Scalar chemicalPotential_[numPhases][numComponents];
 };
 
 } // end namespace