EoSparam Class Reference

EoSParam contains the params for Compositional Model and functions to read them. More...

#include <ParamReservoir.hpp>

Public Member Functions
void	InitEoSparam ()
	Init Params.
void	InputCOM (ifstream &ifs)
	Input the information of components. More...
void	InputCOMPONENTS (ifstream &ifs, const string &keyword)
	Input the information of hydrocarbon components.
Type_A_r< vector< OCP_DBL > > *	FindPtr (const string &varName)
void	InputCNAMES (ifstream &ifs)
	Input the names of hydrocarbon components.
void	InputLBCCOEF (ifstream &ifs)
	Input LBC coefficients for viscosity calculation.
void	InputBIC (ifstream &ifs)
	Input the Binary interaction of components. More...
void	InputSSMSTA (ifstream &ifs)
	TODO: Add Doxygen.
void	InputNRSTA (ifstream &ifs)
	TODO: Add Doxygen.
void	InputSSMSP (ifstream &ifs)
	TODO: Add Doxygen.
void	InputNRSP (ifstream &ifs)
	TODO: Add Doxygen.
void	InputRR (ifstream &ifs)
	TODO: Add Doxygen.

Public Attributes
USI	NTPVT {1}
	num of EoS region, constant now.
USI	numCom {0}
	num of components, water is excluded.
USI	numPhase {2}
	num of phase, water is excluded, constant now.
vector< vector< string > >	COM
	Components information.
vector< string >	Cname
	Name of hydrocarbon components.
Type_A_r< vector< OCP_DBL > >	Tc
	Critical temperature of hydrocarbon components.
Type_A_r< vector< OCP_DBL > >	Pc
	Critical pressure of hydrocarbon components.
Type_A_r< vector< OCP_DBL > >	Vc
	Critical volume of hydrocarbon components.
Type_A_r< vector< OCP_DBL > >	Zc
	Critical Z-factor of hydrocarbon components.
Type_A_r< vector< OCP_DBL > >	MW
	Molecular Weight of hydrocarbon components.
Type_A_r< vector< OCP_DBL > >	Acf
	Acentric factor of hydrocarbon components.
Type_A_r< vector< OCP_DBL > >	OmegaA
	OMEGA_A of hydrocarbon components.
Type_A_r< vector< OCP_DBL > >	OmegaB
	OMEGA_B of hydrocarbon components.
Type_A_r< vector< OCP_DBL > >	Vshift
	Volume shift of hydrocarbon components.
Type_A_r< vector< OCP_DBL > >	Parachor
	PARACHOR of hydrocarbon components.
Type_A_r< vector< OCP_DBL > >	Vcvis
	Critical volume used for viscosity calculations only.
Type_A_r< vector< OCP_DBL > >	Zcvis
	Critical Z-factor used for viscosity calculations only.
vector< OCP_DBL >	LBCcoef
	LBC coefficients for viscosity calculation.
vector< vector< OCP_DBL > >	BIC
	Binary interaction.
bool	miscible {false}
	Miscible treatment of hydrocarbons, used in compositional Model.
vector< string >	SSMparamSTA
	Params for Solving Phase Spliting with SSM.
vector< string >	NRparamSTA
	Params for Solving Phase Spliting with NR.
vector< string >	SSMparamSP
	Params for Solving Phase Spliting with SSM.
vector< string >	NRparamSP
	Params for Solving Phase Spliting with NR.
vector< string >	RRparam
	Params for Solving Rachford-Rice equations.

Detailed Description

EoSParam contains the params for Compositional Model and functions to read them.

Definition at line 65 of file ParamReservoir.hpp.

Member Function Documentation

FindPtr()

Type_A_r< vector< OCP_DBL > > * EoSparam::FindPtr

(

const string &

varName

)

Find corresponding variable according to the name of variable. It is used for the basic properties of hydrocarbon components such as TCRIT

Definition at line 811 of file ParamReservoir.cpp.

{
    Type_A_r<vector<OCP_DBL>>* myPtr = nullptr;
 
    switch (Map_Str2Int(&varName[0], varName.size())) 
    {
    case Map_Str2Int("TCRIT", 5):
        myPtr = &Tc;
        break;
 
    case Map_Str2Int("PCRIT", 5):
        myPtr = &Pc;
        break;
 
    case Map_Str2Int("VCRIT", 5):
        myPtr = &Vc;
        break;
 
    case Map_Str2Int("ZCRIT", 5):
        myPtr = &Zc;
        break;
 
    case Map_Str2Int("MW", 2):
        myPtr = &MW;
        break;
 
    case Map_Str2Int("ACF", 3):
        myPtr = &Acf;
        break;
 
    case Map_Str2Int("OMEGAA", 6):
        myPtr = &OmegaA;
        break;
 
    case Map_Str2Int("OMEGAB", 6):
        myPtr = &OmegaB;
        break;
 
    case Map_Str2Int("SSHIFT", 6):
        myPtr = &Vshift;
        break;
 
    case Map_Str2Int("PARACHOR", 8):
        myPtr = &Parachor;
        break;
 
    case Map_Str2Int("VCRITVIS", 8):
        myPtr = &Vcvis;
        break;
 
    case Map_Str2Int("ZCRITVIS", 8):
        myPtr = &Zcvis;
        break;
 
    }
 
    return myPtr;
}

References Acf, Map_Str2Int(), MW, OmegaA, OmegaB, Parachor, Pc, Tc, Vc, Vcvis, Vshift, Zc, and Zcvis.

InputBIC()

void EoSparam::InputBIC

(

ifstream &

ifs

)

Input the Binary interaction of components.

Input Binary Interaction Coefficients Matrix.

Definition at line 958 of file ParamReservoir.cpp.

{
    OCP_ASSERT((numCom > 0) && (NTPVT > 0), "NCNP hasn't been input!");
  
    BIC.resize(NTPVT);
 
    vector<string> vbuf;
    USI nReg = 0;
    while (true) {
        ReadLine(ifs, vbuf);
        if (vbuf[0] == "/") {
            nReg++;
            if (nReg >= NTPVT)
                break;
            continue;
        }
        for (auto& v : vbuf) {
            if (v != "/") {
                BIC[nReg].push_back(stod(v));
                cout << setw(10) << BIC[nReg].back();
            }                
        }
        cout << endl;
        if (vbuf.back() == "/") {
            nReg++;
            if (nReg >= NTPVT)
                break;
        }
    }
}

References BIC, NTPVT, numCom, OCP_ASSERT, and ReadLine().

InputCOM()

void EoSparam::InputCOM

(

ifstream &

ifs

)

Input the information of components.

TODO: Add Doxygen.

Definition at line 777 of file ParamReservoir.cpp.

{
    OCP_ASSERT(numCom > 0, "Wrong NC!");
    COM.resize(numCom);
    USI len = 9;
 
    vector<string> vbuf;
 
    for (USI c = 0; c < numCom; c++) {
        COM[c].resize(len);
        ReadLine(ifs, vbuf);
        for (USI i = 0; i < len; i++) {
            COM[c][i] = vbuf[i];
        }
    }
    OCP_FUNCNAME;
    cout << "Name    "
         << "Pc                "
         << "Tc           "
         << "Acentric    "
         << "MW               "
         << "Vc            "
         << "OmegaA          "
         << "OmegaB       "
         << "Shift" << endl;
    for (auto& c : COM) {
        for (auto& item : c) {
            cout << item << "\t";
        }
        cout << endl;
    }
    cout << endl;
}

References COM, numCom, OCP_ASSERT, OCP_FUNCNAME, and ReadLine().

---

The documentation for this class was generated from the following files:

• ParamReservoir.hpp
• ParamReservoir.cpp