pzblackoilanalysis.h

Go to the documentation of this file.

#ifndef BLACKOILANALH
#define BLACKOILANALH

#include "pzanalysis.h"
#include "pznonlinanalysis.h"
#include "pzcompel.h"
#include "pzfmatrix.h"
#include "pzvec.h"
#include <iostream>

class TPZCompMesh;
class TPZFStructMatrix;

#ifdef _AUTODIFF

class TPZBlackOilAnalysis : public TPZNonLinearAnalysis {
                
private:
                TPZFMatrix<STATE> fLastState;
                
public:
                TPZBlackOilAnalysis(TPZCompMesh *mesh, double TimeStep, std::ostream &out = std::cout);
                ~TPZBlackOilAnalysis();
                
                virtual void Assemble();
                
                virtual void AssembleResidual();
                
                virtual void Solve();
                
//              virtual void Run(std::ostream &out = std::cout, bool linesearch = true);
    virtual void Run(std::ostream &out,bool linesearch);
    virtual void Run(std::ostream &out = std::cout) {
        TPZBlackOilAnalysis::Run(out,true);
    }
                
                virtual void PostProcess(int resolution){ TPZAnalysis::PostProcess(resolution);}
                
                virtual void PostProcess(int resolution, int dimension);
                
                virtual void PostProcess(TPZVec<REAL> &loc, std::ostream &out = std::cout);
                
                void SetConvergence(int niter, REAL eps, bool ForceAllSteps = true);
                
                void SetSaveFrequency(int SaveFrequency, int resolution);
                
                void SetNewtonConvergence(int niter, REAL eps);
                REAL &TimeStep();
                
                void SetInitialSolution(TPZFMatrix<STATE> & InitialSol);
                void SetInitialSolutionAsZero();
                
protected:
                
                REAL fTimeStep;
                
                REAL fSimulationTime;
                
                int fNIter;
                
                int fCurrentStep;
                
                REAL fSteadyTol;
                
                bool fForceAllSteps;
                
                int fSaveFrequency;
                
                int fDXResolution;
                
                int fNewtonMaxIter;
                
                REAL fNewtonTol;
                
                void SetLastState();
                
                void SetCurrentState();
                
                void SetAllMaterialsDeltaT();
                
                static double PressaoMedia(TPZBlackOilAnalysis &an, int matid);
                static void Vazao(TPZBlackOilAnalysis &an, int matid, double & VazaoAguaSC, double  & VazaoOleoSC, double & VazaoAguaFundo, double  & VazaoOleoFundo);
};

#endif

#endif