densityRelations.h

Go to the documentation of this file.

#ifndef DENSITYRELATIONS_H
#define DENSITYRELATIONS_H


class DensityRelation
{
public:
  double rho,drho;
};
 
class ConstantDensity : public DensityRelation
{
public:
  double rho_0;
  inline ConstantDensity(const double* rwork):
    rho_0(rwork[0])
  {
    rho=rho_0;
    drho=0.0;
  }
  inline void calc(const double& psi)
  {
    rho=rho_0;
    drho=0.0;
  }
};

class LinearDensity : public DensityRelation
{
public:
  double rho_0,psi_0,beta;
  inline LinearDensity(const double* rwork):
    rho_0(rwork[0]),
    psi_0(rwork[1]),
    beta(rwork[2])
  {}
  
  inline void calc(const double& psi)
  {
    rho = rho_0*(1.0+beta*(psi - psi_0));
    drho = beta*rho_0;
  }
};

class ExponentialDensity : public LinearDensity
{
public:
  inline ExponentialDensity(const double* rwork):
    LinearDensity(rwork)
  {}
  inline void calc(const double& psi)
  {
    rho = rho_0*exp(beta*(psi-psi_0));
    drho = beta*rho;
  }
};

class RealGasDensity : public DensityRelation
{
public:
  double psi_0,oneOverZRT;
  inline RealGasDensity(const double *rwork):
    psi_0(rwork[0]),
    oneOverZRT(rwork[1])
  {}
  inline void calc(const double& psi)
  {
    rho = psi*oneOverZRT;
    drho = oneOverZRT;
  }
};

class IdealGasDensity : public DensityRelation
{
public:
  double T,W,R,convFactor,WoverRT,rho0,psi0;
  inline IdealGasDensity(const double *rwork):
    T(rwork[0]),
    W(rwork[1]),
    R(rwork[2]),
    convFactor(rwork[3]),
    rho0(rwork[4]),
    psi0(rwork[5])
    {
      WoverRT = convFactor*W/(R*T);
    }
  inline void calc(const double& psi)
  {
    rho = rho0 + (psi-psi0)*WoverRT;
    drho = WoverRT;
  }
};
 

#endif