TPZDiffusionConsLaw.h

Go to the documentation of this file.

#ifndef DIFFUSIONFORCONSLAWHH
#define DIFFUSIONFORCONSLAWHH


#include "pzvec.h"
#include "pzfmatrix.h"

class TPZDiffusionConsLaw {
                
                
public:
                static REAL fGamma;
                
                static REAL fDelta;
                
                static REAL fCFL;
                
                static std::string fArtificialDiffusion;
                
private:
                TPZFMatrix<REAL> fA,fB,fC;
                
                int fDimension;
                
public:
                TPZDiffusionConsLaw();
                TPZDiffusionConsLaw(TPZVec<REAL> U,REAL gamma,int dim,const std::string &diff);
                ~TPZDiffusionConsLaw();
                
                //static void SetDelta(REAL delta){TPZDiffusionConsLaw::fDelta = delta;}
                
                void GradientOfTheFlow(TPZFMatrix<REAL> &DF1,TPZFMatrix<REAL> &DF2,TPZFMatrix<REAL> &DF3);
                
                REAL DeltaOtimo();
                
                REAL Delta();
                
                REAL CFL(int degree);
                
                static void JacobFlux(TPZVec<REAL> U,TPZFMatrix<REAL> &A,TPZFMatrix<REAL> &B,TPZFMatrix<REAL> &C);
                
                void Divergence(TPZVec<REAL> &dphi,TPZFMatrix<REAL> &diverg);
                
                void PointOperator(TPZVec<REAL> &dphi,TPZFMatrix<REAL> &diff_term);
                
                void Tau(TPZFMatrix<REAL> &Tx,TPZFMatrix<REAL> &Ty,TPZFMatrix<REAL> &Tz);
                
                void SUPG(TPZFMatrix<REAL> &Tx,TPZFMatrix<REAL> &Ty,TPZFMatrix<REAL> &Tz);
                
                void LS(TPZFMatrix<REAL> &Tx,TPZFMatrix<REAL> &Ty,TPZFMatrix<REAL> &Tz);
                
                void Bornhaus(TPZFMatrix<REAL> &Tx,TPZFMatrix<REAL> &Ty,TPZFMatrix<REAL> &Tz);
                
                static void Roe_Flux(REAL rho_f, REAL rhou_f, REAL rhov_f, REAL rhow_f, REAL rhoE_f, 
                                                                                                 REAL rho_t, REAL rhou_t, REAL rhov_t,REAL rhow_t, REAL rhoE_t, 
                                                                                                 REAL nx, REAL ny, REAL nz, REAL gam, REAL & flux_rho, REAL & flux_rhou, 
                                                                                                 REAL & flux_rhov,REAL & flux_rhow, REAL & flux_rhoE);
                
                static void Roe_Flux(REAL rho_f, REAL rhou_f, REAL rhov_f, REAL rhoE_f,
                                                                                                 REAL rho_t, REAL rhou_t, REAL rhov_t, REAL rhoE_t,
                                                                                                 REAL nx, REAL ny, REAL gam,
                                                                                                 REAL &flux_rho, REAL &flux_rhou,
                                                                                                 REAL &flux_rhov, REAL &flux_rhoE);
};

#endif