Inheritance diagram for MixtureComp:

Public Member Functions
OCP_DBL	GetErrorPEC () override

OCP_ULL	GetSSMSTAiters () override

OCP_ULL	GetNRSTAiters () override

OCP_ULL	GetSSMSPiters () override

OCP_ULL	GetNRSPiters () override

OCP_ULL	GetRRiters () override

OCP_ULL	GetSSMSTAcounts () override

OCP_ULL	GetNRSTAcounts () override

OCP_ULL	GetSSMSPcounts () override

OCP_ULL	GetNRSPcounts () override

OCP_ULL	GetRRcounts () override

	MixtureComp (const ParamReservoir &rs_param, const USI &i)

	MixtureComp (const EoSparam &param, const USI &i)

void	InitFlash (const OCP_DBL &Pin, const OCP_DBL &Pbbin, const OCP_DBL &Tin, const OCP_DBL Sjin, const OCP_DBL &Vpore, const OCP_DBL Ziin) override
	flash calculation with saturation of phases.

void	InitFlashDer (const OCP_DBL &Pin, const OCP_DBL &Pbbin, const OCP_DBL &Tin, const OCP_DBL Sjin, const OCP_DBL &Vpore, const OCP_DBL Ziin) override

void	InitFlashDer_n (const OCP_DBL &Pin, const OCP_DBL &Pbbin, const OCP_DBL &Tin, const OCP_DBL Sjin, const OCP_DBL &Vpore, const OCP_DBL Ziin) override

void	Flash (const OCP_DBL &Pin, const OCP_DBL &Tin, const OCP_DBL Niin, const USI &ftype, const USI &lastNP, const OCP_DBL lastKs) override
	Flash calculation with moles of components.

void	CalFlash (const OCP_DBL &Pin, const OCP_DBL &Tin, const OCP_DBL *Niin)

void	FlashDeriv (const OCP_DBL &Pin, const OCP_DBL &Tin, const OCP_DBL Niin, const USI &ftype, const USI &lastNP, const OCP_DBL lastKs) override
	Flash calculation with moles of components and Calculate the derivative.

void	FlashDeriv_n (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 &ftype, const USI phaseExistin, const USI &lastNP, const OCP_DBL lastKs) override

OCP_DBL	XiPhase (const OCP_DBL &Pin, const OCP_DBL &Tin, const OCP_DBL *Ziin) override

OCP_DBL	RhoPhase (const OCP_DBL &Pin, const OCP_DBL &Tin, const OCP_DBL *Ziin) override
	return mass density of phase.

OCP_DBL	GammaPhaseO (const OCP_DBL &Pin, const OCP_DBL &Pbbin) override
	return gamma of oil phase, gamma equals to mass density times gravity factor.

OCP_DBL	GammaPhaseG (const OCP_DBL &Pin) override
	return gamma of gas phase, gamma equals to mass density times gravity factor.

OCP_DBL	GammaPhaseW (const OCP_DBL &Pin) override
	return gamma of water phase, gamma equals to mass density times gravity factor.

OCP_DBL	GammaPhaseOG (const OCP_DBL &Pin, const OCP_DBL &Tin, const OCP_DBL *Ziin) override

void	setPT (const OCP_DBL &p, const OCP_DBL &t)

void	setZi (const OCP_DBL *Ziin)

void	setZi ()

void	setNi (const OCP_DBL *Niin)

void	CallId ()

USI	GetFtype () override

void	CalSurfaceTension ()

OCP_DBL	GetSurTen () override

void	AllocateEoS ()

void	SolEoS (OCP_DBL &ZjT, const OCP_DBL &AjT, const OCP_DBL &BjT) const

void	CalAiBi ()

void	CalAjBj (OCP_DBL &AjT, OCP_DBL &BjT, const vector< OCP_DBL > &xj) const

void	CalAjBj (OCP_DBL &AjT, OCP_DBL &BjT, const OCP_DBL *xj) const

USI	CubicRoot (const OCP_DBL &a, const OCP_DBL &b, const OCP_DBL &c, const bool &NTflag=false) const
	Result is stored in Ztmp.

void	PrintZtmp ()
	test

void	AllocatePhase ()

void	CalFugPhi (vector< OCP_DBL > &phiT, vector< OCP_DBL > &fugT, const vector< OCP_DBL > &xj)

void	CalFugPhi (OCP_DBL phiT, OCP_DBL fugT, const OCP_DBL *xj)

void	CalFugPhi (OCP_DBL fugT, const OCP_DBL xj)

void	CalFugPhiAll ()

void	CalMW ()

void	CalVfXiRho ()

void	CalSaturation ()

USI	FindMWmax ()

void	x2n ()
	x[j][i] -> n[j][i]

void	PrintX ()

void	AllocateMethod ()

void	PhaseEquilibrium ()

void	CalKwilson ()

bool	PhaseStable ()

bool	StableSSM (const USI &Id)
	strict SSM

bool	StableSSM01 (const USI &Id)
	relaxable SSM

bool	StableNR (const USI &Id)

void	CalFugXSTA ()
	Calculate d ln(Fug) / dx for Y.

void	AssembleJmatSTA ()

bool	CheckSplit ()

void	PhaseSplit ()

void	SplitSSM (const bool &flag)

void	SplitSSM2 (const bool &flag)

void	SplitSSM3 (const bool &flag)

void	RachfordRice2 ()
	Used when NP = 2.

void	RachfordRice2P ()
	Used when NP = 2, improved RachfordRice2.

void	RachfordRice3 ()
	Used when NP > 2.

void	UpdateXRR ()
	Update X according to RR.

void	SplitBFGS ()
	Use BFGS to calculate phase splitting.

void	SplitNR ()
	Use NR to calculate phase splitting.

void	CalResSP ()

void	CalFugNAll (const bool &Znflag=true)

void	PrintFugN ()

void	AssembleJmatSP ()

void	CalPhiNSTA ()
	Calculate d ln phi[i][j] / d n[k][j].

void	AssembleSkipMatSTA ()

OCP_DBL	CalStepNRsp ()

OCP_SIN	GetMinEigenSkip () override

bool	GetFlagSkip () override

void	AllocateOthers ()

void	IdentifyPhase ()

void	CopyPhase ()
	Copy the basic properties from MixtureComp to Mixture.

void	CalViscosity ()

void	CalViscoLBC ()

void	CalViscoHZYT ()

void	CalFugXAll ()

void	CalFugPAll (const bool &Zpflag=true)

void	CalVjpVfpVfx_partial ()

void	CalXiPNX_partial ()

void	CalRhoPX_partial ()

void	CalMuPX_partial ()

void	CalMuPXLBC_partial ()

void	CalXiRhoMuPN_pfullx ()

void	CaldXsdXpAPI04 ()

void	CaldXsdXp04 ()

void	CalRhoPNX_full ()

void	CalXiPNX_full01 ()

void	CalRhoPNX_full01 ()

void	CalMuPX_full01 ()

void	CalMuPXLBC_full01 ()

void	CalVfiVfp_full01 ()

void	AssembleMatVfiVfp_full01 ()

void	AssembleRhsVfiVfp_full01 ()

void	CaldXsdXp01 ()

void	CaldXsdXpAPI01 ()

void	CalXiPNX_full02 ()

void	CalVfiVfp_full02 ()

void	AssembleMatVfiVfp_full02 ()

void	AssembleRhsVfiVfp_full02 ()

void	CaldXsdXpAPI02 ()

void	CaldXsdXpAPI02p ()

void	CalVjpVfpVfn_partial ()

void	CalXiPn_partial ()

void	CalRhoPn_partial ()

void	CalMuPn_partial ()

void	CalMuPnLBC_partial ()

void	CalXiRhoMuPN_pfullxn (const bool &xflag=true)

void	CaldXsdXpAPI03 ()

void	CaldXsdXp03 ()

void	CalVfiVfp_full03 ()

void	CalKeyDerx ()

void	CalKeyDern ()

Public Member Functions inherited from Mixture
void	Allocate ()
	Allocate memory for common variables for basic class.

virtual void	SetPVTW ()

USI	GetType () const
	return type of mixture.

virtual bool	IsEmpty_PVDG () const
	Check whether Table PVDG is empty, it will only be used in black oil model.

void	CheckNi (const OCP_DBL *Ni)

Additional Inherited Members
Protected Attributes inherited from Mixture
USI	mixtureType

USI	numPhase
	num of phases.

USI	numCom
	num of components.

OCP_DBL	P
	pressure when flash calculation.

OCP_DBL	T
	temperature when flash calculation.

vector< OCP_DBL >	Ni
	moles of component: numCom

vector< bool >	phaseExist
	existence of phase: numPhase

vector< OCP_DBL >	S
	saturation of phase: numPhase

vector< OCP_DBL >	rho
	mass density of phase: numPhase

vector< OCP_DBL >	xi
	molar density of phase: numPhase

vector< OCP_DBL >	xij

vector< OCP_DBL >	nj
	mole number of phase j

vector< OCP_DBL >	mu
	viscosity of phase: numPhase

vector< OCP_DBL >	v
	volume of phase: numPhase;

OCP_DBL	vf
	volume of total fluids.

OCP_DBL	Nt
	Total moles of Components.

vector< vector< OCP_DBL > >	vji
	dvj / dNi, used in 2 hydrocarbon phase in EOS; or dvj / dnij

vector< OCP_DBL >	vjp
	dvj / dp, used in 2 hydrocarbon phase in EOS

OCP_DBL	vfp

vector< OCP_DBL >	vfi

vector< OCP_DBL >	muP
	d mu / dP: numPhase

vector< OCP_DBL >	xiP
	d xi / dP: numphase

vector< OCP_DBL >	rhoP
	d rho / dP: numphase

vector< OCP_DBL >	muN
	d mu[j] / d N[i]: numphase * numCom

vector< OCP_DBL >	xiN
	d xi[j] / d N[i]: numphase * numCom

vector< OCP_DBL >	rhoN
	d rho[j] / d N[i]: numphase * numCom

vector< OCP_DBL >	mux
	d mu[j] / d x[i][j]: numphase * numCom

vector< OCP_DBL >	xix
	d xi[j] / d x[i][j]: numphase * numCom

vector< OCP_DBL >	rhox
	d rho[j] / d x[i][j]: numphase * numCom

vector< OCP_DBL >	dXsdXp
	the derivates of second variables wrt. primary variables

vector< USI >	pEnumCom
	see pEnumCom in bulk

vector< OCP_DBL >	res
	residual of a set of equations

OCP_DBL	resPc
	a precalculated value

vector< OCP_DBL >	keyDer
	d (xij*xi/mu) / dP or dNk

Detailed Description

Definition at line 127 of file MixtureComp.hpp.

Member Function Documentation

◆ FlashDeriv_n()

void MixtureComp::FlashDeriv_n	(	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 &	ftype,
		const USI *	phaseExistin,
		const USI &	lastNP,
		const OCP_DBL *	lastKs
	)

overridevirtual

Becareful if NP > 2 (temp)

Implements Mixture.

Definition at line 514 of file MixtureComp.cpp.

 {  
     inputNP = 0;
     for (USI j = 0; j < numPhase; j++) {
         if (phaseExistin[j] == 1)
             inputNP++;
     }
     inputNP--;
  
  
     if (inputNP == 1 || true) {
         ftype = myftype;
         lNP = lastNP;
         if (lNP == 2) {
             Dcopy(NC, &lKs[0], lastKs);
         }
         CalFlash(Pin, Tin, Niin);
     }
     else {
         NP = inputNP;
         P = Pin; T = Tin;
         setNi(Niin); 
         CalAiBi();
         Nh = Dnorm1(NC, &Ni[0]);
         for (USI j = 0; j < NP; j++) {
             phaseExist[j] = phaseExistin[j];
             S[j] = Sjin[j];
             nu[j] = njin[j]; 
             nj[j] = njin[j];
             Dcopy(NC, &x[j][0], &xijin[j * numCom]);
             Dcopy(NC, &xij[j * numCom], &xijin[j * numCom]);
         }
         CalFugPhiAll();
         S[numPhase - 1] = Sjin[numPhase - 1];
         phaseLabel[0] = OIL;   phaseLabel[1] = GAS;
         CalMW();
         CalVfXiRho();
         CalViscosity();
         CopyPhase();
         CalSurfaceTension();
     }
  
  
     CalVjpVfpVfn_partial();
  
     // Water Properties
     USI Wpid = numPhase - 1;
     USI Wcid = numCom - 1;
     phaseExist[Wpid] = true;
     xij[Wpid * numCom + Wcid] = 1.0;
     nj[Wpid] = Ni[Wcid];
     Nt = Nh + Ni[Wcid];
     PVTW.Eval_All(0, P, data, cdata);
     OCP_DBL Pw0 = data[0];
     OCP_DBL bw0 = data[1];
     OCP_DBL cbw = data[2];
     OCP_DBL bw = bw0 * (1 - cbw * (P - Pw0));
     OCP_DBL bwp = -cbw * bw0;
     mu[Wpid] = data[3];
     xi[Wpid] = 1 / (CONV1 * bw);
     rho[Wpid] = std_RhoW / bw;
     muP[Wpid] = cdata[3];
     xiP[Wpid] = -bwp / (bw * bw * CONV1);
     rhoP[Wpid] = CONV1 * xiP[Wpid] * std_RhoW;
     v[Wpid] = CONV1 * Ni[Wcid] * bw;
     vfi[Wcid] = CONV1 * bw;
     const OCP_DBL vwp = CONV1 * Ni[Wcid] * bwp;
     vf += v[Wpid];
     vfp += vwp;
     vji[numPhase - 1][numCom - 1] = vfi[Wcid];
     vjp[numPhase - 1] = vwp;
  
     CalSaturation();
  
     // calculate diff    
     if (inputNP == NP && NP == 2 && false) {
         cout << scientific << setprecision(3);
         cout << "--------- " << NP << " --------- " << inputNP << endl;
         for (USI j = 0; j < NP; j++) {
             USI j1 = phaseLabel[j];
             cout << setw(10) << S[j1] - Sjin[j1] << "   "
                 << setw(10) << (nj[j1] - njin[j1]) / njin[j1] << "   ";
             for (USI i = 0; i < NC; i++) {
                 cout << setw(10) << xij[j1 * numCom + i] - xijin[j1 * numCom + i] << "   ";
             }
             cout << endl;
         }
         cout << S[2] - Sjin[2] << endl;
     }
  
  
     CaldXsdXpAPI03();
     CalXiPn_partial();
     CalRhoPn_partial();
     CalMuPn_partial();
     CalVfiVfp_full03();
  
     // Calculate pEnumCom
     fill(pEnumCom.begin(), pEnumCom.end(), 0.0);
     if (phaseExist[0]) pEnumCom[0] = NC;
     if (phaseExist[1]) pEnumCom[1] = NC;
  
 }

References CONV1, CopyPhase(), Dcopy(), Dnorm1(), OCPTable::Eval_All(), GAS, Mixture::mu, Mixture::muP, Mixture::Ni, Mixture::nj, Mixture::Nt, Mixture::numCom, Mixture::numPhase, OIL, Mixture::pEnumCom, Mixture::phaseExist, Mixture::rho, Mixture::rhoP, Mixture::S, Mixture::v, Mixture::vf, Mixture::vfi, Mixture::vfp, Mixture::vji, Mixture::vjp, Mixture::xi, Mixture::xij, and Mixture::xiP.

◆ GammaPhaseOG()

OCP_DBL MixtureComp::GammaPhaseOG	(	const OCP_DBL &	Pin,
		const OCP_DBL &	Tin,
		const OCP_DBL *	Ziin
	)

overridevirtual

return gamma of hydrocarbon mixture, gamma equals to mass density times gravity factor.

Implements Mixture.

Definition at line 704 of file MixtureComp.cpp.

 {
     // assume that only single phase exists here, no matter it's oil or gas
     OCP_DBL rhotmp = RhoPhase(Pin, Tin, Ziin);
     return rhotmp * GRAVITY_FACTOR;
 }

References GRAVITY_FACTOR, and RhoPhase().

◆ XiPhase()

OCP_DBL MixtureComp::XiPhase	(	const OCP_DBL &	Pin,
		const OCP_DBL &	Tin,
		const OCP_DBL *	Ziin
	)

overridevirtual

Return molar density of phase, it's used to calculate the molar density of injection fluids in injection wells.

Implements Mixture.

Definition at line 644 of file MixtureComp.cpp.

 {
     // assume that only single phase exists here
     if (Ziin[numCom - 1] > 1 - 1e-6) {
         // water phase
         PVTW.Eval_All(0, Pin, data, cdata);
         OCP_DBL Pw0   = data[0];
         OCP_DBL bw0   = data[1];
         OCP_DBL cbw   = data[2];
         OCP_DBL bw    = bw0 * (1 - cbw * (P - Pw0));
         OCP_DBL xitmp = 1 / (CONV1 * bw);
         return xitmp;
     } else {
         // hydrocarbon phase
         setPT(Pin, Tin);
         setZi(Ziin);
         NP = 1;
         CalAiBi();
         CalAjBj(Aj[0], Bj[0], zi);
         SolEoS(Zj[0], Aj[0], Bj[0]);
         OCP_DBL vtmp = Zj[0] * GAS_CONSTANT * T / P;
         for (USI i = 0; i < NC; i++) {
             vtmp -= zi[i] * Vshift[i];
         }
         OCP_DBL xitmp = 1 / vtmp;
         return xitmp;
     }
 }

References CONV1, OCPTable::Eval_All(), GAS_CONSTANT, and Mixture::numCom.

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

Public Member Functions

Additional Inherited Members

Detailed Description

Member Function Documentation

◆ FlashDeriv_n()

◆ GammaPhaseOG()

◆ XiPhase()