MixtureThermal.hpp

Go to the documentation of this file.

 
#ifndef __MIXTURETHERMAL_HEADER__
#define __MIXTURETHERMAL_HEADER__
 
#include <cmath>
 
// OpenCAEPoro header files
#include "Mixture.hpp"
#include "OCPTable.hpp"
 
class MixtureThermal : public Mixture
{
public:
    MixtureThermal() = default;
    void Allocate()
    {
        Ni.resize(numCom);
        phaseExist.resize(numPhase);
        S.resize(numPhase);
        vj.resize(numPhase);
        nj.resize(numPhase);
        xij.resize(numPhase * numCom);
        rho.resize(numPhase);
        xi.resize(numPhase);
        mu.resize(numPhase);
        vfi.resize(numCom);
        rhoP.resize(numPhase);
        rhoT.resize(numPhase);
        rhox.resize(numPhase * numCom);
        xiP.resize(numPhase);
        xiT.resize(numPhase);
        xix.resize(numPhase * numCom);
        muP.resize(numPhase);
        muT.resize(numPhase);
        mux.resize(numPhase * numCom);
        dXsdXp.resize((numCom + 1) * numPhase * (numCom + 2));
        Ufi.resize(numCom);
        H.resize(numPhase);
        HT.resize(numPhase);
        Hx.resize(numPhase * numCom);
    }
    void SetupOptionalFeatures(OptionalFeatures& optFeatures,
                               const OCP_USI&    numBulk) override{};
    // usless in Thermal model
    OCP_DBL GetErrorPEC() override
    {
        OCP_ABORT("Should not be used in Thermal mode!");
        return 0;
    }
    void OutMixtureIters() const override{};
};
 
class MixtureThermal_K01 : public MixtureThermal
{
public:
    MixtureThermal_K01() = default;
    MixtureThermal_K01(const ParamReservoir& param, const USI& tarId);
    void Flash(const OCP_DBL& Pin, const OCP_DBL& Tin, const OCP_DBL* Niin) override;
    void InitFlashIMPEC(const OCP_DBL& Pin,
                        const OCP_DBL& Pbbin,
                        const OCP_DBL& Tin,
                        const OCP_DBL* Sjin,
                        const OCP_DBL& Vpore,
                        const OCP_DBL* Ziin,
                        const OCP_USI& bId) override;
    void InitFlashFIM(const OCP_DBL& Pin,
                      const OCP_DBL& Pbbin,
                      const OCP_DBL& Tin,
                      const OCP_DBL* Sjin,
                      const OCP_DBL& Vpore,
                      const OCP_DBL* Ziin,
                      const OCP_USI& bId) override;
    void InitFlashFIMn(const OCP_DBL& Pin,
                       const OCP_DBL& Pbbin,
                       const OCP_DBL& Tin,
                       const OCP_DBL* Sjin,
                       const OCP_DBL& Vpore,
                       const OCP_DBL* Ziin,
                       const OCP_USI& bId) override{};
    void FlashIMPEC(const OCP_DBL& Pin,
                    const OCP_DBL& Tin,
                    const OCP_DBL* Niin,
                    const USI&     lastNP,
                    const OCP_DBL* xijin,
                    const OCP_USI& bId) override;
    void FlashFIM(const OCP_DBL& Pin,
                  const OCP_DBL& Tin,
                  const OCP_DBL* Niin,
                  const OCP_DBL* Sjin,
                  const USI&     lastNP,
                  const OCP_DBL* xijin,
                  const OCP_USI& bId) override;
    void FlashFIMn(const OCP_DBL& Pin,
                   const OCP_DBL& Tin,
                   const OCP_DBL* Niin,
                   const OCP_DBL* Sjin,
                   const OCP_DBL* xijin,
                   const OCP_DBL* njin,
                   const USI*     phaseExistin,
                   const USI&     lastNP,
                   const OCP_USI& bId) override{};
    OCP_DBL XiPhase(const OCP_DBL& Pin,
                    const OCP_DBL& Tin,
                    const OCP_DBL* Ziin,
                    const USI&     tarPhase) override;
 
    OCP_DBL RhoPhase(const OCP_DBL& Pin,
                     const OCP_DBL& Pbb,
                     const OCP_DBL& Tin,
                     const OCP_DBL* Ziin,
                     const USI&     tarPhase) override;
 
    // for Well
    void CalProdWeight(const OCP_DBL&         Pin,
                       const OCP_DBL&         Tin,
                       const OCP_DBL*         Niin,
                       const vector<OCP_DBL>& prodPhase,
                       vector<OCP_DBL>&       prodWeight) override;
 
    void CalProdRate(const OCP_DBL&   Pin,
                     const OCP_DBL&   Tin,
                     const OCP_DBL*   Niin,
                     vector<OCP_DBL>& prodRate) override;
 
    void    SetupWellOpt(WellOpt&                  wellopt,
                         const vector<SolventINJ>& sols,
                         const OCP_DBL&            Psurf,
                         const OCP_DBL&            Tsurf) override;
    OCP_DBL CalInjWellEnthalpy(const OCP_DBL& Tin, const OCP_DBL* Ziin) override;
 
protected:
    void CalEnthalpy();
 
protected:
    OCP_DBL Pref{PRESSURE_STD};    
    OCP_DBL Tref{TEMPERATURE_STD}; 
 
    vector<OCP_DBL> MWc;    
    vector<OCP_DBL> MWp;    
    vector<OCP_DBL> Tcrit;  
    vector<OCP_DBL> xi_ref; 
    vector<OCP_DBL> cp;     
    vector<OCP_DBL> ct1;    
    vector<OCP_DBL> ct2;    
    vector<OCP_DBL> cpt; 
    OCP_BOOL liquid_based{OCP_TRUE}; 
    OCP_BOOL gas_based{OCP_FALSE};   
    OCP_BOOL simple_hvap{OCP_FALSE}; 
    vector<OCP_DBL> cpl1;            
    vector<OCP_DBL> cpl2;            
    vector<OCP_DBL> cpl3;            
    vector<OCP_DBL> cpl4;            
    vector<OCP_DBL> cpg1;            
    vector<OCP_DBL> cpg2;            
    vector<OCP_DBL> cpg3;            
    vector<OCP_DBL> cpg4;            
    vector<OCP_DBL>
        hvapr; 
    vector<OCP_DBL> hvr; 
    vector<OCP_DBL> ev;  
    vector<OCP_DBL>
        avisc; 
    vector<OCP_DBL>
        bvisc; 
    vector<OCP_DBL> avg; 
    vector<OCP_DBL> bvg; 
    OCP_BOOL        useViscTab; 
    OCPTable        visc;       
    vector<OCP_DBL> data;  
    vector<OCP_DBL> cdata; 
};
 
#endif /* end if __MIXTURETHERMAL_HEADER__ */
 
/*----------------------------------------------------------------------------*/
/*  Brief Change History of This File                                         */
/*----------------------------------------------------------------------------*/
/*  Author              Date             Actions                              */
/*----------------------------------------------------------------------------*/
/*  Shizhe Li           NOV/10/2022      Create file                          */
/*----------------------------------------------------------------------------*/