gerg04.hpp

#pragma once
 
#include <armadillo>
#include <string>
 
#include "equationofstate/equationofstatebase.hpp"
 
class GERG04 : public EquationOfStateBase
{
public:
    explicit GERG04(const arma::vec& composition = Composition::defaultComposition);
 
    virtual arma::vec evaluate(const double pressure, const double temperature) const override;
 
    arma::vec evaluateAllProperties(const double pressure, const double temperature) const;
 
    virtual double calculateCompressibility(const double pressure, const double temperature) const override;
 
    virtual bool setComposition(const arma::vec& composition, const bool force = true) override;
 
    double findDensity(const double pressure, const double temperature) const;
 
    double findSpeedOfSound(const double temperature, const double density) const;
 
//    void useBadHeatCapacities(bool useBadHeatCapacities);
 
//    arma::vec getOutput() const;
 
    const arma::uvec& indicesOfNonZeroComponents() const { return m_indices; }
 
private:
    double findDensity(
            const double pressure,
            const double temperature,
            double& density, // output
            double& aroidelta, // output
            double& arijdelta, // output
            double& aroideltadelta, // output
            double& arijdeltadelta // output
            ) const;
 
    void evaluateAlpha_roi_deltas(
            const double tred_temperature,
            const double start_rhored,
            const double start_rhored_pow_minusOne,
            const double start_rhored_pow_minusTwo,
            double& aroidelta, // output
            double& arijdelta, // output
            double& aroideltadelta, // output
            double& arijdeltadelta // output
            ) const;
 
    void setNonZeroComponents();
 
    void calculateCoefficients();
 
    const arma::uword N = 10; 
    arma::vec X; 
 
    arma::vec Tc =   {190.56,         126.19,         304.13,         305.32,         369.83,         425.13,         407.82,         469.7,          460.35,         507.82};
    arma::vec rhoc = {10.14*16.04,    11.18*28.01,    10.62*44.01,    6.87*30.07,     5*44.1,         3.92*58.12,     3.86*58.12,     3.21*72.15,     3.27*72.15,     2.71*86.18};
    //                ch4,            n2,             c02,            c2h6,           c3h8,           nc4h10,         ic4h10,         nc5h12,         ic5h12          nc6h14
 
    arma::uvec m_indices;
 
    arma::uvec m_firstIndices;
 
    arma::uvec m_lastIndices;
 
    static arma::mat betav; 
    static arma::mat betat; 
    static arma::mat gammav; 
    static arma::mat gammat; 
 
    // Kpol coefficients - used in the first sum in alpha^r_{oi} and derivatives
    static arma::mat noipol; 
    static arma::mat doipol; 
    static arma::mat toipol; 
 
    // Kexp coefficient - used in the second sum in alpha^r_{oi} and derivatives
    static arma::mat noiexp; 
    static arma::mat doiexp; 
    static arma::mat coiexp; 
    static arma::mat toiexp; 
 
    static arma::mat Fij; 
 
    static arma::cube nijpol; 
    static arma::cube dijpol; 
    static arma::cube tijpol; 
 
    static arma::cube nijexp; 
    static arma::cube dijexp; 
    static arma::cube tijexp; 
 
    static arma::cube nuijexp; 
    static arma::cube epijexp; 
    static arma::cube beijexp; 
    static arma::cube gaijexp; 
 
    // other coefficients
    static arma::mat noik; 
    static arma::mat voik; 
 
    // constant terms
    static arma::cube nijpol_times_tijpol; 
    static arma::cube nijpol_times_tijpol_times_tijpol_minus_one; 
    static arma::cube nijexp_times_tijexp; 
    static arma::cube nijexp_times_tijexp_times_tijexp_minus_one; 
    static arma::cube nijpol_times_dijpol; 
    static arma::cube nijpol_times_dijpol_times_tijpol; 
    static arma::cube nijpol_times_dijpol_times_dijpol_minus_one; 
 
    double rhored; 
    double tred; 
 
    double Ra; 
};