BulkConn.hpp

Go to the documentation of this file.

 
#ifndef __BULKCONN_HEADER__
#define __BULKCONN_HEADER__
 
// Standard header files
#include <vector>
 
// OpenCAEPoro header files
#include "Bulk.hpp"
#include "DenseMat.hpp"
#include "Grid.hpp"
#include "LinearSystem.hpp"
#include "OCPStructure.hpp"
 
using namespace std;
 
//  Note: Bulks are the active grid cells.
class BulkPair
{
    friend class BulkConn;
 
public:
    BulkPair() = default;
 
    BulkPair(const OCP_USI& BId,
             const OCP_USI& EId,
             const USI&     direct,
             const OCP_DBL& AreaB,
             const OCP_DBL& AreaE)
        : bId(BId)
        , eId(EId)
        , direction(direct)
        , areaB(AreaB)
        , areaE(AreaE){};
 
    OCP_USI BId() const { return bId; }   
    OCP_USI EId() const { return eId; }   
    OCP_DBL Area() const { return area; } 
    OCP_DBL AreaB() const { return areaB; }
    OCP_DBL AreaE() const { return areaE; }
    USI     Direction() const { return direction; }
 
protected:
    OCP_USI bId;       
    OCP_USI eId;       
    OCP_DBL area;      
    USI     direction; 
    OCP_DBL areaB;     
    OCP_DBL areaE;     
};
 
//  Note: BulkConn is a core component of reservoir, it contains all properties and
//  operations on connections between bulks (active grids). You can traverse all the
//  connections through an effective iterator. Flow calculations between active
// bulks and matrix assembling with contributions from bulks only are included.
class BulkConn
{
    friend class Reservoir;
    // temp
    friend class MyMetisTest;
    friend class Out4VTK;
    friend class IsoT_FIM;
    friend class IsoT_IMPEC;
    friend class IsoT_AIMc;
    friend class IsoT_FIMn;
    friend class T_FIM;
 
public:
    // General
 
public:
    void PrintConnectionInfo(const Grid& myGrid) const;
    void PrintConnectionInfoCoor(const Grid& myGrid) const;
 
    // General Variables
 
public:
    void SetupIsoT(const Grid& myGrid, const Bulk& myBulk);
    void Setup(const Grid& myGrid);
    void CalAkd(const Bulk& myBulk);
    OCP_USI     GetBulkNum() const { return numBulk; }
    const auto& GetNeighborNum() const { return neighborNum; }
 
protected:
    OCP_USI numBulk; 
    OCP_USI numConn; 
 
    //  Note: The i-th entry stores the i-th bulk's neighbors, which is sorted in an
    //  increasing order.
    vector<vector<OCP_USI>> neighbor;
 
    vector<USI> selfPtr;
 
    vector<USI> neighborNum;
 
    //  Note: In each pair, the index of first bulk is greater than the second. The data
    //  in iteratorConn is generated from neighbor.
    vector<BulkPair> iteratorConn;
 
    // Physical Variables
 
protected:
    //  Note: Upblock is identified by difference of pressure between phases.
    vector<OCP_USI>
        upblock; 
    vector<OCP_DBL>
        upblock_Rho; 
    vector<OCP_DBL>
        upblock_Trans; 
    vector<OCP_DBL> upblock_Velocity; 
    vector<OCP_DBL> Adkt;             
 
    // Last time step
    vector<OCP_USI> lupblock;          
    vector<OCP_DBL> lupblock_Rho;      
    vector<OCP_DBL> lupblock_Trans;    
    vector<OCP_DBL> lupblock_Velocity; 
    vector<OCP_DBL> lAdkt;             
 
    // Derivatives
    vector<OCP_DBL> AdktP; 
    vector<OCP_DBL> AdktT; 
    vector<OCP_DBL>
        AdktS; 
 
    // Last time step
    vector<OCP_DBL> lAdktP; 
    vector<OCP_DBL> lAdktT; 
    vector<OCP_DBL> lAdktS; 
 
public:
    void CalFluxFIMS(const Grid& myGrid, const Bulk& myBulk);
    void CalResFIMS(vector<OCP_DBL>& res, const Bulk& myBulk, const OCP_DBL& dt);
};
 
#endif
 
/*----------------------------------------------------------------------------*/
/*  Brief Change History of This File                                         */
/*----------------------------------------------------------------------------*/
/*  Author              Date             Actions                              */
/*----------------------------------------------------------------------------*/
/*  Shizhe Li           Oct/01/2021      Create file                          */
/*  Chensong Zhang      Oct/15/2021      Format file                          */
/*  Chensong Zhang      Jan/17/2022      Update Doxygen                       */
/*----------------------------------------------------------------------------*/