fluid_model.hpp

Go to the documentation of this file.

#pragma once

#include "../../Common/include/mpi_structure.hpp"

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <string>
#include <cmath>

#define LEN_COMPONENTS 32

#include "stdio.h"
#include "math.h"

using namespace std;

#include "../include/transport_model.hpp"
#include "../../Common/include/config_structure.hpp"


class CFluidModel {
protected:
su2double      StaticEnergy,      
       Entropy,          
       Density,          
       Pressure,         
       SoundSpeed2,       
       Temperature,      
       dPdrho_e,         
       dPde_rho,         
       dTdrho_e,         
       dTde_rho,         
             Cp,                    
       Mu,          
         dmudrho_T,       
         dmudT_rho,        
         Kt,          
         dktdrho_T,       
         dktdT_rho;        
CViscosityModel *LaminarViscosity;            
CConductivityModel *ThermalConductivity;    
public:

    CFluidModel(void);

    virtual ~CFluidModel(void);

    su2double GetPressure ();

    su2double GetTemperature ();

    su2double GetEntropy ();

    su2double GetStaticEnergy ();

    su2double GetDensity ();

    su2double GetSoundSpeed ();

    su2double GetSoundSpeed2 ();

    su2double GetCp ();

    su2double GetLaminarViscosity ();

    su2double GetThermalConductivity ();

    su2double GetdPdrho_e ();

    su2double GetdPde_rho ();

    su2double GetdTdrho_e ();

    su2double GetdTde_rho ();

    su2double Getdmudrho_T ();

    su2double GetdmudT_rho ();

    su2double Getdktdrho_T ();

    su2double GetdktdT_rho ();

    void SetLaminarViscosityModel (CConfig *config);

    void SetThermalConductivityModel (CConfig *config);

    virtual void SetTDState_rhoe (su2double rho, su2double e );

    virtual void SetTDState_PT (su2double P, su2double T );

    virtual void SetTDState_Prho (su2double P, su2double rho );

    virtual void SetEnergy_Prho (su2double P, su2double rho );

    virtual void SetTDState_hs (su2double h, su2double s );


    virtual void SetTDState_rhoT (su2double rho, su2double T );


    virtual void SetTDState_Ps (su2double P, su2double s );

};


class CIdealGas : public CFluidModel {

protected:
  su2double Gamma,             
          Gamma_Minus_One,       
          Gas_Constant;        
public:

    CIdealGas(void);

    CIdealGas(su2double gamma, su2double R);


    virtual ~CIdealGas(void);

    void SetTDState_rhoe (su2double rho, su2double e );

    void SetTDState_PT (su2double P, su2double T );

    void SetTDState_Prho (su2double P, su2double rho );

    void SetEnergy_Prho (su2double P, su2double rho );

    void SetTDState_hs (su2double h, su2double s );


    void SetTDState_rhoT (su2double rho, su2double T );

    void SetTDState_Ps (su2double P, su2double s );
};


class CVanDerWaalsGas : public CIdealGas {

protected:
  su2double
      a, b, Zed;             
public:

    CVanDerWaalsGas(void);

    CVanDerWaalsGas(su2double gamma, su2double R, su2double Pstar, su2double Tstar);


    virtual ~CVanDerWaalsGas(void);

    void SetTDState_rhoe (su2double rho, su2double e );

    void SetTDState_PT (su2double P, su2double T );

    void SetTDState_Prho (su2double P, su2double rho );

    void SetEnergy_Prho (su2double P, su2double rho );

    void SetTDState_hs (su2double h, su2double s );


    void SetTDState_rhoT (su2double rho, su2double T );

    void SetTDState_Ps (su2double P, su2double s );

};


class CPengRobinson : public CIdealGas {

protected:
  su2double  a,             
        b,             
        k,             
        Zed,             
        TstarCrit;        
private:

      su2double  alpha2 (su2double T);


    su2double  T_v_h (su2double v, su2double h);
    su2double T_P_rho(su2double P, su2double rho);



public:

    CPengRobinson(void);

    CPengRobinson(su2double gamma, su2double R, su2double Pstar, su2double Tstar, su2double w);

    virtual ~CPengRobinson(void);

    void SetTDState_rhoe (su2double rho, su2double e );

    void SetTDState_PT (su2double P, su2double T );

    void SetTDState_Prho (su2double P, su2double rho );

    void SetEnergy_Prho (su2double P, su2double rho );
    void SetTDState_hs (su2double h, su2double s );

    void SetTDState_rhoT (su2double rho, su2double T );

    void SetTDState_Ps (su2double P, su2double s );

};



#include "fluid_model.inl"