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 OCP_DBL& Area)
        : BId(bId)
        , EId(eId)
        , area(Area){};
 
    OCP_USI GetBId() const { return BId; } 
    OCP_USI GetEId() const { return EId; } 
 
private:
    OCP_USI BId; 
    OCP_USI EId; 
    OCP_DBL area;
};
 
//  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
{
 
public:
    BulkConn() = default; 
 
    // General
 
public:
    void Setup(const Grid& myGrid, const Bulk& myBulk);
 
    void SetupWellBulk_K(Bulk& myBulk) const;
 
    void AllocateMat(LinearSystem& myLS) const;
 
    void SetupMatSparsity(LinearSystem& myLS) const;
 
    void UpdateLastStep();
 
    void Reset();
 
    void CheckDiff() const;
 
    OCP_USI GetBulkNum() const { return numBulk; }
 
    void PrintConnectionInfo(const Grid& myGrid) const;
    void PrintConnectionInfoCoor(const Grid& myGrid) const;
 
private:
    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;
 
    // sum of num of neighbors who has the bigger Index than current bulk
    // used to find the location in iteratorConn, upblock, etc.
    vector<USI> neighborNumGacc;
 
    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;
 
 
 
    //  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_USI> lastUpblock;          
    vector<OCP_DBL> lastUpblock_Rho;      
    vector<OCP_DBL> lastUpblock_Trans;    
    vector<OCP_DBL> lastUpblock_Velocity; 
 
    // IMPEC
 
public:
    void AllocateAuxIMPEC(const USI& np);
 
    void AssembleMatIMPEC(LinearSystem& myLS, const Bulk& myBulk,
                          const OCP_DBL& dt) const;
 
    void CalCFL(const Bulk& myBulk, const OCP_DBL& dt) const;
 
    void CalFluxIMPEC(const Bulk& myBulk);
 
    void MassConserveIMPEC(Bulk& myBulk, const OCP_DBL& dt) const;
 
    // FIM
 
public:
    void AllocateAuxFIM(const USI& np);
 
    void AssembleMat_FIM(LinearSystem& myLS, const Bulk& myBulk,
                         const OCP_DBL& dt) const;
 
    void CalFluxFIM(const Bulk& myBulk);
 
    void CalResFIM(vector<OCP_DBL>& res, const Bulk& myBulk, const OCP_DBL& dt);
 
    void CalFluxFIMS(const Bulk& myBulk);
    void CalResFIMS(vector<OCP_DBL>& res, const Bulk& myBulk, const OCP_DBL& dt);
 
 
    // FIM(new)
 
    void AssembleMat_FIM_new(LinearSystem& myLS, const Bulk& myBulk,
        const OCP_DBL& dt) const;
    void AssembleMat_FIM_newS(LinearSystem& myLS, const Bulk& myBulk,
        const OCP_DBL& dt) const;
    void AssembleMat_FIM_new_n(LinearSystem& myLS, const Bulk& myBulk,
        const OCP_DBL& dt) const;
 
    // AIMs, AIMt
 
public:
    void SetupFIMBulk(Bulk& myBulk, const bool& NRflag = false) const;
    void AddFIMBulk(Bulk& myBulk);
    void SetupFIMBulkBoundAIMs(Bulk& myBulk);
    void AllocateAuxAIMt();
    void SetupMatSparsityAIMt(LinearSystem& myLS, const Bulk& myBulk) const;
    void AssembleMat_AIMt(LinearSystem& myLS, const Bulk& myBulk,
        const OCP_DBL& dt) const;
    void CalResAIMt(vector<OCP_DBL>& res, const Bulk& myBulk, const OCP_DBL& dt);
    void CalResAIMs(vector<OCP_DBL>& res, const Bulk& myBulk, const OCP_DBL& dt);
    void AssembleMat_AIMs(LinearSystem& myLS, vector<OCP_DBL>& res, const Bulk& myBulk,
        const OCP_DBL& dt) const;
 
    void AllocateAuxAIMc(const USI& np);
    void AssembleMat_AIMc(LinearSystem& myLS, const Bulk& myBulk, const OCP_DBL& dt) const;
    void AssembleMat_AIMc01(LinearSystem& myLS, const Bulk& myBulk, const OCP_DBL& dt) const;
    void CalResAIMc(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                       */
/*----------------------------------------------------------------------------*/