pzspacetimerichardseq.h

Go to the documentation of this file.

#ifndef PZSPACETIMERICHARDSEQ_H
#define PZSPACETIMERICHARDSEQ_H

#include "TPZMaterial.h"

class TPZCompMesh;

class TPZSpaceTimeRichardsEq : public TPZMaterial
{
    static STATE TCoeff, LCoeff, deltaDerivada;
protected:
                
                STATE fAlpha, fN, fThetaS, fThetaR, fKs;
                
                STATE Se(STATE sol);
                
                STATE Theta(STATE Se);
                
                STATE DKDsol(STATE sol);
                
                STATE DCDsol(STATE sol);
                
public:
                
    TPZSpaceTimeRichardsEq();
    
    TPZSpaceTimeRichardsEq(int id);
                
    TPZSpaceTimeRichardsEq(int matid, STATE Alpha, STATE N, STATE ThetaS, STATE ThetaR, STATE Ks);
                
    ~TPZSpaceTimeRichardsEq();
                
    void Set(STATE Alpha, STATE N, STATE ThetaS, STATE ThetaR, STATE Ks);
                
                virtual int Dimension() const override;
                
                virtual int NStateVariables() const override;
                
                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;
                
                STATE C_Coef(STATE sol);
                
                STATE K_Coef(STATE sol);
                
                void AnalysisOfParameters(STATE sol0, STATE solL, char* filename);
                
                static int main();
                
                static TPZCompMesh * CreateMesh(REAL L, REAL Time, int p, int ndiv);
                
                static void DirichletT0(TPZVec<REAL> &x, TPZVec<STATE> &f);
    public:
virtual int ClassId() const override;

                
};

#endif