equationofstatebase.hpp

#pragma once
 
#include <armadillo>
 
#include "composition.hpp"
 
class EquationOfStateBase
{
public:
    virtual ~EquationOfStateBase();
 
    explicit EquationOfStateBase(const arma::vec& composition = Composition::defaultComposition);
 
    virtual arma::vec evaluate(const double pressure, const double temperature) const = 0;
 
    virtual arma::vec evaluate(const double pressure, const double temperature, const arma::vec& composition);
 
    virtual double calculateCompressibility(const double pressure, const double temperature) const = 0;
 
    virtual double calculateStandardDensity() const;
 
    virtual bool setComposition(const arma::vec& composition, const bool force = true);
 
    double getMolarMassOfMixture() const; 
    const arma::vec& getComposition() const; 
 
protected:
    // components               C1,         C2,         C3,         iC4,        nC4,        iC5,        nC5,        C6          N2,         CO2
    arma::vec m_molarMass = {   16.04,      30.07,      44.1,       58.12,      58.12,      72.15,      72.15,      86.18,      28.13,      44.01}; // molar weight of components [g/mol]
    arma::vec::fixed<10> m_composition = arma::vec(10);
 
    double m_molarMassOfMixture;
 
    mutable double m_density = 0;
};