pzincnskeps.h

Go to the documentation of this file.

#ifndef PZINCNSKEPS
#define PZINCNSKEPS

#include "TPZMaterial.h"

#include <iostream>
#include <string>

#ifdef _AUTODIFF
#include "fadType.h"
#endif

class TPZBndCond;

class TPZIncNavierStokesKEps : public TPZMaterial {
                
private:
                
    int fDimension;
    
    STATE fMU, fRHO, fCmu, fSigmaK, fSigmaEps, fCepsilon1, fCepsilon2;
    
    TPZVec<STATE> fBodyForce; //fc = {0,0,-fc}
    
    STATE Dot(TPZFMatrix<STATE> &A, TPZFMatrix<STATE> &B);
    
    STATE Dot(TPZVec<STATE> &A, TPZFMatrix<STATE> &B, int BRow);
                
    STATE Dot(TPZVec<STATE> &A, TPZVec<STATE> &B);
                
public:
                
    enum VARIABLES {ENone = -1, EK = 0, EEpsilon, EPressure, EVx, EVy, EVz, EVvector};
                
    TPZIncNavierStokesKEps(int id, int dimension);
    
    void SetParameters(STATE MU, STATE RHO, STATE Cmu, STATE SigmaK, STATE SigmaEps, STATE Cepsilon1, STATE Cepsilon2, TPZVec<STATE> &BodyForce );
    
    void GetParameters(STATE &MU, STATE &RHO, STATE &Cmu, STATE &SigmaK, STATE &SigmaEps, STATE &Cepsilon1, STATE &Cepsilon2, TPZVec<STATE> &BodyForce );    
                
    virtual ~TPZIncNavierStokesKEps();
                
    virtual int Dimension() const override;
                
    virtual int NStateVariables() const override;
                
    virtual void Print(std::ostream &out = std::cout) override;
                
    virtual int NSolutionVariables(int var) override;
                
protected:
                virtual void Solution(TPZVec<STATE> &Sol, TPZFMatrix<STATE> &DSol,
                                                                                                  TPZFMatrix<REAL> &axes, int var, TPZVec<STATE> &Solout) override;
public:
                virtual void Solution(TPZMaterialData &data, int var, TPZVec<STATE> &Solout) override
                {
                                TPZMaterial::Solution(data,var,Solout);
                }
                
                
    virtual void Contribute(TPZMaterialData &data,
                                                                                                                REAL weight,
                                                                                                                TPZFMatrix<STATE> &ek,
                                                                                                                TPZFMatrix<STATE> &ef) override;
                
    virtual void Contribute(TPZMaterialData &data,
                                                                                                                REAL weight,
                                                                                                                TPZFMatrix<STATE> &ef) override;
                
                
                virtual void ContributeBC(TPZMaterialData &data,
                                                                                                                  REAL weight,
                                                                                                                  TPZFMatrix<STATE> &ek,
                                                                                                                  TPZFMatrix<STATE> &ef,
                                                                                                                  TPZBndCond &bc) override;
                
                virtual void ContributeBC(TPZMaterialData &data,
                                                                                                                  REAL weight,
                                                                                                                  TPZFMatrix<STATE> &ef,
                                                                                                                  TPZBndCond &bc) override
                {
        TPZMaterial::ContributeBC(data,weight,ef,bc);
    }
                
                
    virtual void Errors(TPZVec<REAL> &x, TPZVec<STATE> &sol, TPZFMatrix<STATE> &dsol,
                        TPZFMatrix<REAL> &axes, TPZVec<STATE> &flux,
                        TPZVec<STATE> &uexact, TPZFMatrix<STATE> &duexact,
                        TPZVec<REAL> &val) override {
        PZError << __PRETTY_FUNCTION__ << std::endl;
        PZError << "Method not implemented! Error comparison not available. Please, implement it." << std::endl;
    }
    public:
virtual int ClassId() const override;

};


#endif