pzblackoil2p3d.h

Go to the documentation of this file.

#ifndef PZBLACKOIL2P3D_H
#define PZBLACKOIL2P3D_H

#include "TPZMaterial.h"
#include "pzdiscgal.h"

#ifdef _AUTODIFF

#include "fad.h"


class TPZBlackOil2P3D : public TPZDiscontinuousGalerkin{
                
public:
                
                typedef TFad<4, REAL> BFadREAL;
                
protected:
                
                void Interpolate(std::map<REAL,REAL> &dados, double x, double &y, double &dy);
                void Interpolate(std::map<REAL,REAL> &dados, BFadREAL x, BFadREAL &y);
                
                double fDeltaT;
                
                enum EState { ELastState = 0, ECurrentState = 1 };
                
                static EState gState;
                
                void testedoBo();
                
                void testeKrw();
                
public:
                
                static void SetLastState(){ gState = ELastState; }
                static void SetCurrentState(){ gState = ECurrentState; }
                
                TPZBlackOil2P3D(int id, double deltaT);
                
                ~TPZBlackOil2P3D();
                
                TPZBlackOil2P3D(const TPZBlackOil2P3D &cp);
                
                void SetTimeStep(double timestep){
                                this->fDeltaT = timestep;
                }
                
                virtual TPZMaterial * NewMaterial() override;
                
                virtual int Dimension() const override { return 3; }
                
                virtual int NStateVariables() const override{ return 2; }
                
                virtual void Contribute(TPZMaterialData &data, REAL weight, TPZFMatrix<STATE> &ek, TPZFMatrix<STATE> &ef) override;
                
                virtual void ContributeBC(TPZMaterialData &data, REAL weight, TPZFMatrix<STATE> &ek, TPZFMatrix<STATE> &ef, TPZBndCond &bc) override;
                
                virtual void ContributeInterface(TPZMaterialData &data, TPZMaterialData &dataleft, TPZMaterialData &dataright, REAL weight, TPZFMatrix<STATE> &ek, TPZFMatrix<STATE> &ef) override;
                
                virtual void ContributeBCInterface(TPZMaterialData &data, TPZMaterialData &dataleft, REAL weight, TPZFMatrix<STATE> &ek,TPZFMatrix<STATE> &ef,TPZBndCond &bc) override;
                
                
                enum ESolutionVars { ENone = 0, EWaterPressure = 1, EOilPressure, EWaterSaturation, EOilSaturation, EDarcyVelocity };
                
                virtual int VariableIndex(const std::string &name) override;
                
                virtual int NSolutionVariables(int var) override;
                
                virtual void Solution(TPZVec<STATE> &Sol, TPZFMatrix<STATE> &DSol,
                                                                                                  TPZFMatrix<REAL> &axes, int var, TPZVec<STATE> &Solout) override;
                virtual void FillDataRequirements(TPZMaterialData &data) override {
                                data.SetAllRequirements(true);
                                data.fNeedsNeighborSol = false;
                                data.fNeedsNeighborCenter = false;
                }
                
                virtual void FillDataRequirementsInterface(TPZMaterialData &data) override {
                                data.SetAllRequirements(true);
                                data.fNeedsSol = false;
                }
                
                void Kro(double So, double &Kro, double &dKroSo);
                void Kro(BFadREAL So, BFadREAL &Kro);
                
                void Krw(double So, double &Krw, double &dKrwSo);
                void Krw(BFadREAL So, BFadREAL &Krw);
                
                void Bo(double po, double &Bo, double &dBoDpo);
                void Bo(BFadREAL po, BFadREAL &Bo);
                
                void ViscOleo(double po, double &ViscOleo, double &dViscOleoDpo);
                void ViscOleo(BFadREAL po, BFadREAL &ViscOleo);
                
                void PressaoCapilar(double So, double &pc, double &DpcDSo);
                void PressaoCapilar(BFadREAL So, BFadREAL &pc);
                
                void Porosidade(double po, double &poros, double &dPorosDpo);
                void Porosidade(BFadREAL po, BFadREAL &poros);
                
                double RhoOleoSC();
                
                double RhoAguaSC();
                
                double g();
                
                double Bw();
                
                double ViscAgua();
                
                void K(TPZFMatrix<REAL> &K);
                
};

#endif

#endif