PhasePermeability.cpp

Go to the documentation of this file.

 
#include "PhasePermeability.hpp"
 
// Miscible For Compositional Model
 
void Miscible::InputParam(const Miscstr& misterm)
{
    const USI len = misterm.surTenRef.size();
    if (len > 0) {
        ifUseMiscible = OCP_TRUE;
        surTenRef     = misterm.surTenRef[0];
        surTenPc      = 1;
        if (len > 2) {
            surTenPc = misterm.surTenRef[2] / surTenRef;
        }
        Fkexp = 0.25;
    }
}
 
void Miscible::Setup(const OCP_USI& numBulk)
{
    if (!ifSetup) {
        ifSetup = OCP_TRUE;
 
        surTen.resize(numBulk);
        Fk.resize(numBulk);
        Fp.resize(numBulk);
 
        // Last time step
        lsurTen.resize(numBulk);
    }
}
 
OCP_BOOL Miscible::CalFkFp(const OCP_USI& n, OCP_DBL& fk, OCP_DBL& fp)
{
    if (surTen[n] >= surTenRef || surTen[n] <= TINY) {
        Fk[n] = 1;
        Fp[n] = 1;
        return OCP_FALSE; // InMiscible
    } else {
        Fk[n] = fk = min(1.0, pow(surTen[n] / surTenRef, Fkexp));
        Fp[n] = fp = min(surTenPc, surTen[n] / surTenRef);
        return OCP_TRUE; // Miscible
    }
}
 
// Scale The Water-Oil Capillary Pressure Curves (From SWATINIT)
 
void ScalePcow::Setup(const Grid& myGrid)
{
    if (!ifSetup) {
        ifSetup = OCP_TRUE;
 
        if (myGrid.SwatInit.size() > 0) {
            ifScale = OCP_TRUE;
            scaleVal.resize(myGrid.activeGridNum, 1.0);
            swatInit.resize(myGrid.activeGridNum);
 
            for (OCP_USI bIda = 0; bIda < myGrid.activeGridNum; bIda++) {
                OCP_USI bId    = myGrid.map_Act2All[bIda];
                swatInit[bIda] = myGrid.SwatInit[bId];
            }
        }
    }
}
 
/*----------------------------------------------------------------------------*/
/*  Brief Change History of This File                                         */
/*----------------------------------------------------------------------------*/
/*  Author              Date             Actions                              */
/*----------------------------------------------------------------------------*/
/*  Shizhe Li           Dec/26/2022      Create file                          */
/*----------------------------------------------------------------------------*/