IncompressibleFluid.h

Go to the documentation of this file.

/*
 * CoolPropFluid.h
 *
 *  Created on: 20 Dec 2013
 *      Author: jowr
 */
 
#ifndef INCOMPRESSIBLEFLUID_H_
#define INCOMPRESSIBLEFLUID_H_
 
#include "CoolProp/DataStructures.h"
#include "CoolProp/fluids/Helmholtz.h"
#include "CoolProp/numerics/Solvers.h"
 
#include <numeric>
#include <string>
#include <utility>
#include <vector>
#include <map>
#include <cassert>
#include <iterator>
 
#include <Eigen/Core>
#include "CoolProp/numerics/PolyMath.h"
#include "CoolProp/numerics/MatrixMath.h"
 
namespace CoolProp {
 
struct IncompressibleData
{
    enum IncompressibleTypeEnum
    {
        INCOMPRESSIBLE_NOT_SET,
        INCOMPRESSIBLE_POLYNOMIAL,
        INCOMPRESSIBLE_EXPPOLYNOMIAL,
        INCOMPRESSIBLE_EXPONENTIAL,
        INCOMPRESSIBLE_LOGEXPONENTIAL,
        INCOMPRESSIBLE_POLYOFFSET
    };
    IncompressibleTypeEnum type;
    Eigen::MatrixXd coeffs;  //TODO: Can we store the Eigen::Matrix objects more efficiently?
    //std::vector<std::vector<double> > coeffs;
    IncompressibleData()
      : type(INCOMPRESSIBLE_NOT_SET) {
 
        };
};
 
#if !defined(NO_FMTLIB) && FMT_VERSION >= 90000
inline int format_as(IncompressibleData::IncompressibleTypeEnum type) {
    return fmt::underlying(type);
}
#endif
 
 
class IncompressibleFluid
{
 
   protected:
    bool strict;
 
    std::string name;
    std::string description;
    std::string reference;
 
    double Tmin, Tmax;
    double xmin, xmax;
    composition_types xid;
 
    double TminPsat;
    double xbase, Tbase;
 
 
 
    IncompressibleData density;
 
    IncompressibleData specific_heat;
 
    IncompressibleData viscosity;
 
    IncompressibleData conductivity;
 
    IncompressibleData p_sat;
 
    IncompressibleData T_freeze;
 
 
    IncompressibleData mass2input;
 
    IncompressibleData volume2input;
 
    IncompressibleData mole2input;
 
    Polynomial2DFrac poly;
 
    // Forward declaration of the some internal functions
    //double h_u(double T, double p, double x);
    //double u_h(double T, double p, double x);
 
   public:
    IncompressibleFluid()
      : strict(true), Tmin(_HUGE), Tmax(_HUGE), xmin(_HUGE), xmax(_HUGE), xid(IFRAC_UNDEFINED), TminPsat(_HUGE), xbase(_HUGE), Tbase(_HUGE) {
 
        };
    virtual ~IncompressibleFluid() = default;
 
    std::string getName() const {
        return name;
    }
    std::string get_name() const {
        return getName();
    }  // For backwards-compatibility.
    std::string getDescription() const {
        return description;
    }
    std::string getReference() const {
        return reference;
    }
 
    double getTmax() const {
        return Tmax;
    }
    double getTmin() const {
        return Tmin;
    }
    double getxmax() const {
        return xmax;
    }
    double getxmin() const {
        return xmin;
    }
    composition_types getxid() const {
        return xid;
    }
    double getTminPsat() const {
        return TminPsat;
    }
    double getTbase() const {
        return Tbase;
    }
    double getxbase() const {
        return xbase;
    }
 
    void setName(const std::string& name) {
        this->name = name;
    }
    void setDescription(const std::string& description) {
        this->description = description;
    }
    void setReference(const std::string& reference) {
        this->reference = reference;
    }
    void setTmax(double Tmax) {
        this->Tmax = Tmax;
    }
    void setTmin(double Tmin) {
        this->Tmin = Tmin;
    }
    void setxmax(double xmax) {
        this->xmax = xmax;
    }
    void setxmin(double xmin) {
        this->xmin = xmin;
    }
    void setxid(composition_types xid) {
        this->xid = xid;
    }
    void setTminPsat(double TminPsat) {
        this->TminPsat = TminPsat;
    }
    void setTbase(double Tbase) {
        this->Tbase = Tbase;
    }
    void setxbase(double xbase) {
        this->xbase = xbase;
    }
 
    void setDensity(IncompressibleData density) {
        this->density = std::move(density);
    }
    void setSpecificHeat(IncompressibleData specific_heat) {
        this->specific_heat = std::move(specific_heat);
    }
    void setViscosity(IncompressibleData viscosity) {
        this->viscosity = std::move(viscosity);
    }
    void setConductivity(IncompressibleData conductivity) {
        this->conductivity = std::move(conductivity);
    }
    void setPsat(IncompressibleData p_sat) {
        this->p_sat = std::move(p_sat);
    }
    void setTfreeze(IncompressibleData T_freeze) {
        this->T_freeze = std::move(T_freeze);
    }
 
    void setMass2input(IncompressibleData mass2input) {
        this->mass2input = std::move(mass2input);
    }
    void setVolume2input(IncompressibleData volume2input) {
        this->volume2input = std::move(volume2input);
    }
    void setMole2input(IncompressibleData mole2input) {
        this->mole2input = std::move(mole2input);
    }
 
    void validate();
    bool is_pure();
 
   protected:
    double baseExponential(const IncompressibleData& data, double y, double ybase);
    double baseLogexponential(const IncompressibleData& data, double y, double ybase);
    double baseExponentialOffset(IncompressibleData data, double y);
    double basePolyOffset(IncompressibleData data, double y, double z = 0.0);
 
   public:
    /* All functions need T and p as input. Might not
     * be necessary, but gives a clearer structure.
     */
    double rho(double T, double p, double x);
    double c(double T, double p, double x);
    double cp(double T, double p, double x) {
        throw ValueError(format("%s (%d): Please use the c-function instead.", __FILE__, __LINE__));
    }
    double cv(double T, double p, double x) {
        throw ValueError(format("%s (%d): Please use the c-function instead.", __FILE__, __LINE__));
    }
    double s(double T, double p, double x) {
        throw ValueError(format("%s (%d): The internal calculations have changed, use the backend to calculate entropy from the partial derivatives.",
                                __FILE__, __LINE__));
    }
    double u(double T, double p, double x) {
        throw ValueError(
          format("%s (%d): The internal calculations have changed, use the backend to calculate internal energy from enthalpy.", __FILE__, __LINE__));
    }
    double h(double T, double p, double x) {
        throw ValueError(
          format("%s (%d): The internal calculations have changed, use the backend to calculate enthalpy from the partial derivatives.", __FILE__,
                 __LINE__));
    }
    double visc(double T, double p, double x);
    double cond(double T, double p, double x);
    double psat(double T, double x);
    double Tfreeze(double p, double x);
 
    /* Below are direct calculations of the derivatives. Nothing
     * special is going on, we simply use the polynomial class to
     * derive the different functions with respect to temperature.
     */
    //  with respect to temperature at constant pressure and composition
    double drhodTatPx(double T, double p, double x);
    //double dsdTatPx  (double T, double p, double x){return c(T,p,x)/T;};
    //double dhdTatPx  (double T, double p, double x){return c(T,p,x);};
    //  with respect to temperature at constant pressure and composition
    //  integrated in temperature
    double dsdTatPxdT(double T, double p, double x);
    //  with respect to temperature at constant pressure and composition
    //  integrated in temperature
    double dhdTatPxdT(double T, double p, double x);
 
 
    double inputFromMass(double T, double x);
 
    double inputFromVolume(double T, double x);
 
    double inputFromMole(double T, double x);
 
    /* Some functions can be inverted directly, those are listed
     * here. It is also possible to solve for other quantities, but
     * that involves some more sophisticated processing and is not
     * done here, but in the backend, T(h,p) for example.
     */
    double T_rho(double Dmass, double p, double x);
    double T_c(double Cmass, double p, double x);
    double T_s(double Smass, double p, double x) {
        throw NotImplementedError(format("%s (%d): T from entropy is not implemented in the fluid, use the backend.", __FILE__, __LINE__));
    }
    double T_u(double Umass, double p, double x) {
        throw NotImplementedError(format("%s (%d): T from internal energy is not implemented in the fluid, use the backend.", __FILE__, __LINE__));
    }
    double T_h(double Hmass, double p, double x) {
        throw NotImplementedError(format("%s (%d): T from enthalpy is not implemented in the fluid, use the backend.", __FILE__, __LINE__));
    }
    double T_visc(double visc, double p, double x) {
        throw NotImplementedError(format("%s (%d): T from viscosity is not implemented.", __FILE__, __LINE__));
    }
    double T_cond(double cond, double p, double x) {
        throw NotImplementedError(format("%s (%d): T from conductivity is not implemented.", __FILE__, __LINE__));
    }
    double T_psat(double psat, double x) {
        throw NotImplementedError(format("%s (%d): T from psat is not implemented.", __FILE__, __LINE__));
    }
    double x_Tfreeze(double Tfreeze, double p) {
        throw NotImplementedError(format("%s (%d): x from T_freeze is not implemented.", __FILE__, __LINE__));
    }
 
   protected:
    /* Define internal energy and enthalpy as functions of the
     * other properties to provide data in case there are no
     * coefficients.
     */
 
    double h_u(double T, double p, double x) {
        return u(T, p, x) + p / rho(T, p, x);
    };
 
 
    double u_h(double T, double p, double x) {
        return h(T, p, x) - p / rho(T, p, x);
    };
 
    /*
     * Some more functions to provide a single implementation
     * of important routines.
     * We start with the check functions that can validate input
     * in terms of pressure p, temperature T and composition x.
     */
 
    bool checkT(double T, double p, double x);
 
 
    bool checkP(double T, double p, double x);
 
   public:
 
    bool checkX(double x);
 
    bool checkTPX(double T, double p, double x) {
        return (checkT(T, p, x) && checkP(T, p, x) && checkX(x));
    };
};
 
} /* namespace CoolProp */
#endif /* INCOMPRESSIBLEFLUID_H_ */