idealgas.hpp

#pragma once
 
#include <armadillo>
 
#include "constants.hpp"
#include "equationofstate/equationofstatebase.hpp"
 
class IdealGas : public EquationOfStateBase
{
public:
    explicit IdealGas(const Composition&) {}
 
    virtual arma::vec evaluate(const double, const double) const override
    {
        arma::vec Z = arma::zeros<arma::vec>(6);
        Z(0) = 1.0; // Z factor
        Z(1) = 0.0; // derivatives are zero
        Z(2) = 0.0; // derivatives are zero
        Z(3) = 0.0; // derivatives are zero
        Z(4) = 5/2*constants::gasConstant; // c_p of ideal monoatomic gas
        Z(5) = 3/2*constants::gasConstant; // c_v of ideal monoatomic gas
        return Z;
    }
 
    virtual double calculateCompressibility(const double, const double) const override
    {
        return 1.0;
    }
};