pzeuler.h

Go to the documentation of this file.

#ifndef PZEULER_H
#define PZEULER_H

#include <iostream>
#include "TPZMaterial.h"
#include "pzfmatrix.h"
#include "pzvec.h"
#include "pzdiscgal.h"

#include "pzausmflux.h"
#include "pzgradientflux.h"
#include "pzlog.h"
class TPZCompMesh;

class TPZEulerEquation : public TPZDiscontinuousGalerkin{
                
public:
                enum BCType{EFreeSlip = 1};
                enum CALCType{EFlux = 1, EGradient = 2};
                static void SetComputeFlux(){
                                gType = EFlux;
                }
                
                static void SetComputeGradient(){
                                gType = EGradient;
                }
                
#ifdef LinearConvection
                static void SetLinearConvection(TPZCompMesh * cmesh, TPZVec<REAL> &Celerity);
#endif
                
private:
                
#ifdef LinearConvection
                TPZVec<STATE> fCelerity;
#endif
                
                static CALCType gType;
                
                static STATE gGamma;
                
                TPZAUSMFlux fAUSMFlux;
                
                TPZGradientFlux fGradientFlux;
                
                void ComputeEulerFlux(TPZVec<STATE> &sol, TPZFMatrix<STATE> & F);
                
public:
                
                static REAL Gamma(){ return gGamma; }
                
                static void FromPrimitiveToConservative(TPZVec<STATE> &sol,STATE gamma);
                
                static void FromConservativeToPrimitive(TPZVec<STATE> &sol,STATE gamma);
                
                TPZEulerEquation(int nummat, STATE gamma);
                
                ~TPZEulerEquation();
                
                TPZEulerEquation();
                
                TPZEulerEquation(const TPZEulerEquation &cp);
                
                TPZMaterial * NewMaterial() override;
                
                virtual int NStateVariables() const override;
                
                virtual int Dimension() const override;
                
                static REAL Pressure(TPZVec<STATE> &U, double gamma);
                
                STATE cSpeed(TPZVec<STATE> & sol);
                
                STATE uRes(TPZVec<STATE> & sol);
                
                virtual void Print(std::ostream & out) override;
                
                virtual std::string Name() override {return "TPZEulerEquation";}
                
                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 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 ContributeInterface(TPZMaterialData &data, TPZMaterialData &dataleft, TPZMaterialData &dataright, REAL weight, TPZFMatrix<STATE> &ek, TPZFMatrix<STATE> &ef) override;
                
                virtual void ContributeInterface(TPZMaterialData &data, TPZMaterialData &dataleft, TPZMaterialData &dataright, REAL weight, TPZFMatrix<STATE> &ef) override;
                
                virtual void ContributeBC(TPZMaterialData &data,
                                                                                                                  REAL weight,
                                                                                                                  TPZFMatrix<STATE> &ek, TPZFMatrix<STATE> &ef,
                                                                                                                  TPZBndCond &bc) override;
                
                virtual void ContributeBCInterface(TPZMaterialData &data, TPZMaterialData &dataleft,
                                                                                                                                                   REAL weight,
                                                                                                                                                   TPZFMatrix<STATE> &ek,TPZFMatrix<STATE> &ef,
                                                                                                                                                   TPZBndCond &bc) override;
                
                virtual void ContributeBCInterface(TPZMaterialData &data, TPZMaterialData &dataleft,
                                                                                                                                                   REAL weight,
                                                                                                                                                   TPZFMatrix<STATE> &ef,
                                                                                                                                                   TPZBndCond &bc) override;
    public:
virtual int ClassId() const override;

};

#endif