MixtureBO2_OW.cpp

Go to the documentation of this file.

 
#include "MixtureBO.hpp"
 
 // BOMixture_OW
 
BOMixture_OW::BOMixture_OW(const ParamReservoir& rs_param, const USI& i)
{
    BOMixtureInit(rs_param);
 
    PVTW.Setup(rs_param.PVTW_T.data[i]);
    PVDO.Setup(rs_param.PVDO_T.data[i]);
 
    data.resize(5, 0);
    cdata.resize(5, 0);
}
 
void BOMixture_OW::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)
{
    phaseExist[0] = true;
    phaseExist[1] = true;
 
    P    = Pin;
    S[1] = Sjin[1];
    // Water Properties
    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[1]  = data[3];
    xi[1]  = 1 / (bw * CONV1);
    rho[1] = std_RhoW / bw;
 
    Ni[1] = Vpore * S[1] * xi[1];
 
    // Oil Properties
    PVDO.Eval_All(0, P, data, cdata);
    OCP_DBL bo  = data[1];
    OCP_DBL bop = cdata[1];
 
    mu[0]  = data[2];
    xi[0]  = 1 / (CONV1 * bo);
    rho[0] = std_RhoO / bo;
    Ni[0]  = Vpore * (1 - S[1]) * xi[0];
 
    xij[0 * 2 + 0] = 1;
    xij[0 * 2 + 1] = 0;
    xij[1 * 2 + 0] = 0;
    xij[1 * 2 + 1] = 1;
 
    v[0]   = CONV1 * Ni[0] * bo;
    v[1]   = CONV1 * Ni[1] * bw;
    vf     = v[0] + v[1];
    S[0]   = v[0] / vf;
    S[1]   = v[1] / vf;
    vfi[0] = CONV1 * bo;
    vfi[1] = CONV1 * bw;
    vfp    = CONV1 * Ni[0] * bop + CONV1 * Ni[1] * bwp;
}
 
void BOMixture_OW::Flash(const OCP_DBL& Pin, const OCP_DBL& Tin, const OCP_DBL* Niin, const USI& ftype, const USI& lastNP,
    const OCP_DBL* lastKs)
{
    phaseExist[0] = true;
    phaseExist[1] = true;
 
    P     = Pin;
    Ni[0] = Niin[0];
    Ni[1] = Niin[1];
 
    // Water Properties
    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[1]  = data[3];
    xi[1]  = 1 / (bw * CONV1);
    rho[1] = std_RhoW / bw;
 
    // Oil Properties
    PVDO.Eval_All(0, P, data, cdata);
    OCP_DBL bo  = data[1];
    OCP_DBL bop = cdata[1];
 
    mu[0]  = data[2];
    xi[0]  = 1 / (CONV1 * bo);
    rho[0] = std_RhoO / bo;
 
    xij[0 * 2 + 0] = 1;
    xij[0 * 2 + 1] = 0;
    xij[1 * 2 + 0] = 0;
    xij[1 * 2 + 1] = 1;
 
    v[0]   = CONV1 * Ni[0] * bo;
    v[1]   = CONV1 * Ni[1] * bw;
    vf     = v[0] + v[1];
    S[0]   = v[0] / vf;
    S[1]   = v[1] / vf;
    vfi[0] = CONV1 * bo;
    vfi[1] = CONV1 * bw;
    vfp    = CONV1 * Ni[0] * bop + CONV1 * Ni[1] * bwp;
}
 
void BOMixture_OW::FlashDeriv(const OCP_DBL& Pin, const OCP_DBL& Tin,
    const OCP_DBL* Niin, const USI& ftype, const USI& lastNP,
    const OCP_DBL* lastKs)
{
    phaseExist[0] = true;
    phaseExist[1] = true;
    fill(dXsdXp.begin(), dXsdXp.end(), 0.0);
    fill(pEnumCom.begin(), pEnumCom.end(), 0);
 
    P     = Pin;
    Ni[0] = Niin[0];
    Ni[1] = Niin[1];
    Nt    = Ni[0] + Ni[1];
 
    // Water Properties
    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[1]  = data[3];
    xi[1]  = 1 / (bw * CONV1);
    rho[1] = std_RhoW / bw;
 
    muP[1]  = cdata[3];
    xiP[1]  = -bwp / (bw * bw * CONV1);
    rhoP[1] = CONV1 * xiP[1] * std_RhoW;
 
    // Oil Properties
    PVDO.Eval_All(0, P, data, cdata);
    OCP_DBL bo  = data[1];
    OCP_DBL bop = cdata[1];
 
    mu[0]  = data[2];
    xi[0]  = 1 / (CONV1 * bo);
    rho[0] = std_RhoO / bo;
 
    muP[0]  = cdata[2];
    xiP[0]  = -xi[0] * bop / bo;
    rhoP[0] = -rho[0] * bop / bo;
 
    xij[0 * 2 + 0] = 1;
    xij[0 * 2 + 1] = 0;
    xij[1 * 2 + 0] = 0;
    xij[1 * 2 + 1] = 1;
 
    v[0]   = CONV1 * Ni[0] * bo;
    v[1]   = CONV1 * Ni[1] * bw;
    vf     = v[0] + v[1];
    S[0]   = v[0] / vf;
    S[1]   = v[1] / vf;
    vfi[0] = CONV1 * bo;
    vfi[1] = CONV1 * bw;
    vfp    = CONV1 * Ni[0] * bop + CONV1 * Ni[1] * bwp;
 
    dXsdXp[0] = (CONV1 * Ni[0] * bop - S[0] * vfp) / vf; // dSo / dP
    dXsdXp[1] = (CONV1 * bo - S[0] * vfi[0]) / vf;       // dSo / dNo
    dXsdXp[2] = -S[0] * vfi[1] / vf;                     // dSo / dNw
 
    dXsdXp[3] = (CONV1 * Ni[1] * bwp - S[1] * vfp) / vf; // dSw / dP
    dXsdXp[4] = -S[1] * vfi[0] / vf;                     // dSw / dNo
    dXsdXp[5] = (CONV1 * bw - S[1] * vfi[1]) / vf;       // dSw / dNw
 
    
}
 
OCP_DBL BOMixture_OW::XiPhase(const OCP_DBL& Pin, const OCP_DBL& Tin, const OCP_DBL* Ziin)
{
    if (Ziin[1] > 1 - TINY) {
        // inj fluid is water
        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 xiw = 1 / (CONV1 * bw);
        return xiw;
    } else {
        OCP_ABORT("Wrong Zi!");
    }
}
 
OCP_DBL BOMixture_OW::RhoPhase(const OCP_DBL& Pin, const OCP_DBL& Tin, const OCP_DBL* Ziin)
{
    if (Ziin[1] > 1 - TINY) {
        // inj fluid is water
 
        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 rhow = std_RhoW / bw;
        return rhow;
    } else {
        OCP_ABORT("Wrong Zi!");
    }
}
 
OCP_DBL BOMixture_OW::GammaPhaseO(const OCP_DBL& Pin, const OCP_DBL& Pbbin)
{
    OCP_DBL bo     = PVDO.Eval(0, Pin, 1);
    OCP_DBL gammaO = std_GammaO / bo;
 
    return gammaO;
}
 
OCP_DBL BOMixture_OW::GammaPhaseW(const OCP_DBL& Pin)
{
 
    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 * (Pin - Pw0)));
 
    return std_GammaW / bw;
}
/*----------------------------------------------------------------------------*/
/*  Brief Change History of This File                                         */
/*----------------------------------------------------------------------------*/
/*  Author              Date             Actions                              */
/*----------------------------------------------------------------------------*/
/*  Shizhe Li           Oct/01/2021      Create file                          */
/*----------------------------------------------------------------------------*/