VTPRBackend.h

Go to the documentation of this file.

//
//  VTPRBackend.h
//  CoolProp
//
//  Created by Ian on 7/17/16.
//
//
 
#ifndef VTPRBackend_h
#define VTPRBackend_h
 
#include "CoolProp/detail/tools.h"
#include "CoolProp/DataStructures.h"
#include "GeneralizedCubic.h"
#include "CoolProp/AbstractState.h"
#include "Backends/Helmholtz/HelmholtzEOSMixtureBackend.h"
#include "CoolProp/Exceptions.h"
#include <vector>
#include "CubicBackend.h"
#include "CoolProp/Configuration.h"
#include "UNIFACLibrary.h"
#include "UNIFAC.h"
#include "VTPRCubic.h"
 
namespace CoolProp {
 
class VTPRBackend : public PengRobinsonBackend
{
 
   private:
    std::vector<double> Tc, pc, omega, molemass, m_ii;
    double R;
    std::vector<std::string> m_fluid_names;
 
   public:
    VTPRBackend(const std::vector<std::string>& fluid_identifiers, const std::vector<double>& Tc, const std::vector<double>& pc,
                const std::vector<double>& acentric, double R_u, bool generate_SatL_and_SatV = true) {
        const UNIFACLibrary::UNIFACParameterLibrary& lib = LoadLibrary();
        cubic = std::make_shared<VTPRCubic>(Tc, pc, acentric, R_u, lib);
        setup(fluid_identifiers, generate_SatL_and_SatV);
    };
    VTPRBackend(const std::vector<std::string>& fluid_identifiers, const double R_u = get_config_double(R_U_CODATA),
                bool generate_SatL_and_SatV = true) {
        std::vector<double> Tc, pc, acentric;
        N = fluid_identifiers.size();
        components.resize(N);
        // Extract data from the UNIFAC parameter library
        const UNIFACLibrary::UNIFACParameterLibrary& lib = LoadLibrary();
        for (const auto& fluid_identifier : fluid_identifiers) {
            UNIFACLibrary::Component comp = lib.get_component("name", fluid_identifier);
            Tc.push_back(comp.Tc);              // [K]
            pc.push_back(comp.pc);              // [Pa]
            acentric.push_back(comp.acentric);  // [-]
            molemass.push_back(comp.molemass);  // [kg/mol]
        }
        cubic = std::make_shared<VTPRCubic>(Tc, pc, acentric, R_u, lib);
        setup(fluid_identifiers, generate_SatL_and_SatV);
    };
 
    std::string backend_name() override {
        return get_backend_string(VTPR_BACKEND);
    }
 
    HelmholtzEOSMixtureBackend* get_copy(bool generate_SatL_and_SatV = true) override {
        AbstractCubicBackend* ACB =
          new VTPRBackend(calc_fluid_names(), cubic->get_Tc(), cubic->get_pc(), cubic->get_acentric(), cubic->get_R_u(), generate_SatL_and_SatV);
        ACB->copy_k(this);
        ACB->copy_all_alpha_functions(this);
        return static_cast<HelmholtzEOSMixtureBackend*>(ACB);
    }
    void set_alpha_from_components();
 
    std::vector<std::string> calc_fluid_names() override {
        return m_fluid_names;
    }
 
    void setup(const std::vector<std::string>& names, bool generate_SatL_and_SatV = true);
 
    const UNIFACLibrary::UNIFACParameterLibrary& LoadLibrary();
 
    void set_mole_fractions(const std::vector<double>& z) override {
        mole_fractions = z;
        auto* _cubic = static_cast<VTPRCubic*>(cubic.get());
        _cubic->get_unifaq().set_mole_fractions(z);
    };
 
    CoolPropDbl calc_molar_mass() override;
 
    void set_binary_interaction_double(const std::size_t i, const std::size_t j, const std::string& parameter, const double value) override;
 
    double get_binary_interaction_double(const std::size_t i, const std::size_t j, const std::string& parameter) override;
 
    CoolPropDbl calc_fugacity_coefficient(std::size_t i) override;
 
    void set_Q_k(const size_t sgi, const double value);
};
 
}; /* namespace CoolProp */
 
#endif /* VTPRBackend_h */