OCPFluidMethod.hpp

Go to the documentation of this file.

 
#ifndef __OCPFLUIDMETHOD_HEADER__
#define __OCPFLUIDMETHOD_HEADER__
 
// OpenCAEPoro header files
#include "LinearSystem.hpp"
#include "OCPControl.hpp"
#include "Reservoir.hpp"
#include "UtilTiming.hpp"
#include "UtilOutput.hpp"
 
class OCP_IMPEC
{
public:
    void Setup(Reservoir& rs, LinearSystem& myLS, const OCPControl& ctrl);
 
    void InitReservoir(Reservoir& rs) const;
 
    void Prepare(Reservoir& rs, OCP_DBL& dt);
 
    void SolveLinearSystem(LinearSystem& myLS, Reservoir& rs, OCPControl& ctrl);
 
    bool UpdateProperty(Reservoir& rs, OCPControl& ctrl);
    bool UpdateProperty01(Reservoir& rs, OCPControl& ctrl);
 
    bool FinishNR(const Reservoir& rs);
    bool FinishNR01(Reservoir& rs, OCPControl& ctrl);
 
    void FinishStep(Reservoir& rs, OCPControl& ctrl);
 
};
 
class OCP_FIM
{
public:
    void Setup(Reservoir& rs, LinearSystem& myLS, const OCPControl& ctrl);
 
    void InitReservoir(Reservoir& rs) const;
 
    void Prepare(Reservoir& rs, OCP_DBL& dt);
 
    void AssembleMat(LinearSystem& myLS, const Reservoir& rs, const OCP_DBL& dt) const;
 
    void SolveLinearSystem(LinearSystem& myLS, Reservoir& rs, OCPControl& ctrl) const;
 
    bool UpdateProperty(Reservoir& rs, OCPControl& ctrl);
 
    bool FinishNR(Reservoir& rs, OCPControl& ctrl);
 
    void FinishStep(Reservoir& rs, OCPControl& ctrl) const;
 
 
protected:
    ResFIM resFIM;
};
 
 
class OCP_FIMn : public OCP_FIM
{
public:
 
    void InitReservoir(Reservoir& rs) const;
 
    void AssembleMat(LinearSystem& myLS, const Reservoir& rs, const OCP_DBL& dt) const;
 
    void SolveLinearSystem(LinearSystem& myLS, Reservoir& rs, OCPControl& ctrl) const;
 
    bool UpdateProperty(Reservoir& rs, OCPControl& ctrl);
 
};
 
class OCP_AIMc : public OCP_FIM
{
public:
    void Setup(Reservoir& rs, LinearSystem& myLS, const OCPControl& ctrl);
 
    void InitReservoir(Reservoir& rs) const;
 
    void Prepare(Reservoir& rs, OCP_DBL& dt);
 
    void AssembleMat(LinearSystem& myLS, const Reservoir& rs, const OCP_DBL& dt) const;
 
    void SolveLinearSystem(LinearSystem& myLS, Reservoir& rs, OCPControl& ctrl);
 
    bool UpdateProperty(Reservoir& rs, OCPControl& ctrl);
 
    bool FinishNR(Reservoir& rs, OCPControl& ctrl);
};
 
class OCP_AIMs
{
public:
    void Setup(Reservoir& rs, LinearSystem& myLS, const OCPControl& ctrl);
    void Prepare(Reservoir& rs, OCP_DBL& dt);
    void AssembleMat(LinearSystem& myLS, const Reservoir& rs, const OCP_DBL& dt);
    void SolveLinearSystem(LinearSystem& myLS, Reservoir& rs, OCPControl& ctrl);
    bool UpdateProperty(Reservoir& rs, OCPControl& ctrl);
    bool FinishNR(Reservoir& rs, OCPControl& ctrl);
    void FinishStep(Reservoir& rs, OCPControl& ctrl);
 
private:
    ResFIM resFIM;
 
};
 
class OCP_AIMt
{
public:
    void Setup(Reservoir& rs, LinearSystem& myLS, LinearSystem& myAuxLS, const OCPControl& ctrl);
    void Prepare(Reservoir& rs, OCP_DBL& dt);
    bool UpdateProperty(Reservoir& rs, OCPControl& ctrl, LinearSystem& myAuxLS);
 
private:
    ResFIM resFIM;
 
};
 
 
#endif /* end if __OCPFLUIDMETHOD_HEADER__ */
 
/*----------------------------------------------------------------------------*/
/*  Brief Change History of This File                                         */
/*----------------------------------------------------------------------------*/
/*  Author              Date             Actions                              */
/*----------------------------------------------------------------------------*/
/*  Shizhe Li           Oct/01/2021      Create file                          */
/*  Chensong Zhang      Oct/15/2021      Format file                          */
/*----------------------------------------------------------------------------*/