Merge pull request #1162 from totto82/add_solvent_to_pu

Implement solvent model in flow_ebos

opm/core/props/BlackoilPhases.hpp

@@ -38,6 +38,7 @@ namespace Opm
         int num_phases;
         int phase_used[MaxNumPhases];
         int phase_pos[MaxNumPhases];
+        bool has_solvent;
     };
 
     /// Check or assign presence of a formed, free phase.  Limited to

opm/core/props/phaseUsageFromDeck.hpp

@@ -66,6 +66,12 @@ namespace Opm
             OPM_THROW(std::runtime_error, "Cannot handle cases with no OIL, i.e. water-gas systems.");
         }
 
+        // Add solvent info
+        pu.has_solvent = false;
+        if (phase.active(Phase::SOLVENT)) {
+            pu.has_solvent = true;
+        }
+
         return pu;
     }
 
@@ -112,6 +118,12 @@ namespace Opm
             OPM_THROW(std::runtime_error, "Cannot handle cases with no OIL, i.e. water-gas systems.");
         }
 
+        // Add solvent info
+        pu.has_solvent = false;
+        if (phase.active(Phase::SOLVENT)) {
+            pu.has_solvent = true;
+        }
+
         return pu;
     }
 

opm/core/simulator/BlackoilState.cpp

@@ -11,6 +11,7 @@ using namespace Opm;
 const std::string BlackoilState::GASOILRATIO = "GASOILRATIO";
 const std::string BlackoilState::RV = "RV";
 const std::string BlackoilState::SURFACEVOL = "SURFACEVOL";
+const std::string BlackoilState::SSOL = "SSOL";
 
 
 BlackoilState::BlackoilState( size_t num_cells , size_t num_faces , size_t num_phases)
@@ -19,6 +20,7 @@ BlackoilState::BlackoilState( size_t num_cells , size_t num_faces , size_t num_p
     registerCellData( GASOILRATIO , 1 );
     registerCellData( RV, 1 );
     registerCellData( SURFACEVOL, num_phases );
+    registerCellData( SSOL , 1 );
     setBlackoilStateReferencePointers();
 }
 

opm/core/simulator/BlackoilState.hpp

@@ -43,6 +43,7 @@ namespace Opm
         static const std::string GASOILRATIO;
         static const std::string RV;
         static const std::string SURFACEVOL;
+        static const std::string SSOL;
 
         /// Main constructor setting the sizes for the contained data
         /// types.