Esempio n. 1
0
class C3Cyclohexane(MEoS):
    """Multiparameter equation of state for propylcyclohexane"""
    name = "propylcyclohexane"
    CASNumber = "1678-92-8"
    formula = "C6H11-CH2CH2CH3"
    synonym = ""
    _refPropName = "C2CC6"
    _coolPropName = ""
    rhoc = unidades.Density(260.0527932)
    Tc = unidades.Temperature(630.8)
    Pc = unidades.Pressure(2860.0, "kPa")
    M = 126.23922  # g/mol
    Tt = unidades.Temperature(178.2)
    Tb = unidades.Temperature(429.9)
    f_acent = 0.326
    momentoDipolar = unidades.DipoleMoment(0.0, "Debye")
    id = 184

    CP1 = {
        "an": [9.29427],
        "pow": [0.385871],
        "ao_exp": [1.37051, 106.426, 313.713],
        "exp": [173295, 561.14, 1919.52]
    }

    lemmon = {
        "__type__": "Helmholtz",
        "__name__": "short Helmholtz equation of state for propylcyclohexane "
        "of Lemmon (2007).",
        "__doi__": {
            "autor": "Lemmon, E.W.",
            "title": "",
            "ref": "unpublished equation, 2007",
            "doi": ""
        },
        "R": 8.314472,
        "cp": CP1,
        "ref": "NBP",
        "Tmin": Tt,
        "Tmax": 650.,
        "Pmax": 50000.0,
        "rhomax": 7.03,
        "nr1": [1.01911, -2.59762, 0.675152, -0.230891, 0.120966, 0.000309038],
        "d1": [1, 1, 1, 2, 3, 7],
        "t1": [0.2, 1.2, 1.8, 1.5, 0.3, 0.9],
        "nr2": [0.526461, -0.0188462, -0.549272, -0.139233, 0.121242],
        "d2": [2, 5, 1, 4, 1],
        "t2": [1.4, 2.2, 3.7, 4.2, 2.4],
        "c2": [1, 1, 2, 2, 1],
        "gamma2": [1] * 5
    }

    eq = lemmon,

    _vapor_Pressure = {
        "eq": 3,
        "n": [-0.76296e1, 0.16538e1, -0.28518e1, -0.28205e1, -0.28144e1],
        "t": [1.0, 1.5, 2.7, 4.7, 15.0]
    }
    _liquid_Density = {
        "eq": 1,
        "n": [0.39271e-1, 0.38257e2, -0.65743e2, 0.30332e2, 0.17224],
        "t": [0.1, 0.75, 0.87, 1.0, 5.0]
    }
    _vapor_Density = {
        "eq": 2,
        "n": [-0.64572e1, 0.91228e1, -0.25806e2, -0.59044e2, -0.14709e3],
        "t": [0.6, 1.8, 2.2, 6.0, 14.0]
    }

    thermo0 = {
        "__name__":
        "Perkins (2008)",
        "__doi__": {
            "autor":
            "Perkins, R.A. Hammerschmidt, U., Huber, M.L.",
            "title":
            "Measurement and Correlation of the Thermal "
            "Conductivity of Methylcyclohexane and "
            "Propylcyclohexane from 300 to 600 K at Pressures "
            "to 60 MPa",
            "ref":
            "J. Chem. Eng. Data 53(9) (2008) 2120-2127",
            "doi":
            "10.1021/je800255r"
        },
        "eq":
        1,
        "Toref":
        630.8,
        "koref":
        1,
        "no": [1.07402e-2, -6.09829e-2, 1.38204e-1, -3.81213e-2],
        "to": [0, 1, 2, 3],
        "Tref_res":
        630.8,
        "rhoref_res":
        260.05,
        "kref_res":
        1.,
        "nr": [
            0.116524, -0.102821, -0.113871, 0.126431, 0.0445827, -0.05946,
            -0.00545736, 0.0098936
        ],
        "tr": [0, -1, 0, -1, 0, -1, 0, -1],
        "dr": [1, 1, 2, 2, 3, 3, 4, 4],
        "critical":
        3,
        "gnu":
        0.63,
        "gamma":
        1.2415,
        "R0":
        1.01,
        "Xio":
        0.15e-9,
        "gam0":
        0.052,
        "qd":
        6.24e-10,
        "Tcref":
        958.725
    }

    _thermal = thermo0,
Esempio n. 2
0
class iC8(MEoS):
    """Multiparameter equation of state for isooctane"""
    name = "isooctane"
    CASNumber = "540-84-1"
    formula = "(CH3)2CHCH2C(CH3)3"
    synonym = ""
    rhoc = unidades.Density(242.1644624)
    Tc = unidades.Temperature(544)
    Pc = unidades.Pressure(2572.0, "kPa")
    M = 114.22852  # g/mol
    Tt = unidades.Temperature(165.77)
    Tb = unidades.Temperature(372.358)
    f_acent = 0.303
    momentoDipolar = unidades.DipoleMoment(0.0, "Debye")
    id = 82

    CP1 = {"ao": 10.76,
           "an": [], "pow": [],
           "ao_exp": [15.48, 34.42, 21.42],
           "exp": [775, 1900, 5100],
           "ao_hyp": [], "hyp": []}

    helmholtz1 = {
        "__type__": "Helmholtz",
        "__name__": "Helmholtz equation of state for isooctane of Blackham and Lemmon (2011).",
        "__doi__": {"autor": "Blackham, T.M. and Lemmon, E.W.",
                    "title": "",
                    "ref": "to be published in Int. J. Thermophys., 2011.",
                    "doi": ""},

        "R": 8.314472,
        "cp": CP1,
        "ref": "NBP",

        "Tmin": Tt, "Tmax": 600.0, "Pmax": 1000000.0, "rhomax": 6.97,
        "Pmin": 0.00001796, "rhomin": 6.96,

        "nr1": [0.568901e-1, 0.196155e1, -0.281164e1, -0.815112, 0.326583],
        "d1": [4, 1, 1, 2, 3],
        "t1": [1, 0.3, 0.75, 1.11, 0.55],

        "nr2": [-0.160893e1, -0.454734, 0.108306e1, -0.722876, -0.434052e-1],
        "d2": [1, 3, 2, 2, 7],
        "t2": [2.2, 3.7, 1.53, 2.1, 0.9],
        "c2": [2, 2, 1, 2, 1],
        "gamma2": [1]*6,

        "nr3": [0.196648e1, -0.465082, -0.409398, 0.232131e-1],
        "d3": [1, 1, 3, 3, ],
        "t3": [0.88, 1.1, 2.75, 1.0],
        "alfa3": [0.75, 1.13, 0.87, 4.73],
        "beta3": [0.59, 1.45, 0.5, 10.52],
        "gamma3": [1.44, 0.68, 0.51, 0.8],
        "epsilon3": [0.66, 0.9, 0.54, 0.18],
        "nr4": []}

    eq = helmholtz1,
    _PR = -0.0058658

    _surface = {"sigma": [0.0476182, 5.992036e-17], "exp": [1.1914, 2.1914]}
    _vapor_Pressure = {
        "eq": 5,
        "ao": [-7.7985, 8.1280, -7.3106, -3.9392, -1.6732],
        "exp": [1, 1.5, 1.6, 4.0, 16.0]}
    _liquid_Density = {
        "eq": 1,
        "ao": [1.1535, 1.3709, 0.38804],
        "exp": [0.286, 0.54, 3.3]}
    _vapor_Density = {
        "eq": 3,
        "ao": [-2.5793, -6.4934, -18.631, -54.123, -123.58],
        "exp": [0.366, 1.11, 3.0, 6.4, 14.0]}

    thermo0 = {"eq": 5, "omega": 3,
               "__name__": "Chung (1988)",
               "__doi__": {"autor": "T-H. Chung, Ajlan, M., Lee, L.L. and Starling, K.E.",
                           "title": "Generalized Multiparameter Correlation for Nonpolar and Polar Fluid Transport Properties",
                           "ref": "Ind. Eng. Chem. Res., 1988, 27 (4), pp 671–679",
                           "doi": "10.1021/ie00076a024"},
               "w": 0.3, "mur": 0.0, "k": 0.0}

    _viscosity = thermo0,
    _thermal = thermo0,
Esempio n. 3
0
class R116(MEoS):
    """Multiparameter equation of state for R116"""
    name = "hexafluoroethane"
    CASNumber = "76-16-4"
    formula = "CF3CF3"
    synonym = "R116"
    rhoc = unidades.Density(613.3245)
    Tc = unidades.Temperature(293.03)
    Pc = unidades.Pressure(3048.0, "kPa")
    M = 138.01182  # g/mol
    Tt = unidades.Temperature(173.1)
    Tb = unidades.Temperature(195.06)
    f_acent = 0.2566
    momentoDipolar = unidades.DipoleMoment(0.0, "Debye")
    id = 236

    Fi1 = {"ao_log": [1, 3.],
           "pow": [0, 1],
           "ao_pow": [-10.7088650331, 8.9148979056],
           "ao_exp": [2.4818, 7.0622, 7.9951],
           "titao": [190/Tc, 622/Tc, 1470/Tc]}

    CP1 = {"ao": 4,
           "an": [], "pow": [],
           "ao_exp": [2.4818, 7.0622, 7.9951],
           "exp": [190, 655, 1470],
           "ao_hyp": [], "hyp": []}

    CP2 = {"ao": 27.4009901,
           "an": [-2.6057376855e-6, 9.7501305219e-10, -6559.250418,
                  787904.9649, -34166787.86],
           "pow": [2, 3, -1.001, -2, -3],
           "ao_exp": [], "exp": [],
           "ao_hyp": [], "hyp": []}

    helmholtz1 = {
        "__type__": "Helmholtz",
        "__name__": "short Helmholtz equation of state for R-116 of Lemmon and Span (2006)",
        "__doi__": {"autor": "Lemmon, E.W., Span, R.",
                    "title": "Short Fundamental Equations of State for 20 Industrial Fluids",
                    "ref": "J. Chem. Eng. Data, 2006, 51 (3), pp 785–850",
                    "doi":  "10.1021/je050186n"},
        "__test__": """
            >>> st=R116(T=295, rho=4*138.01182)
            >>> print "%0.0f %0.0f %0.3f %0.3f %0.3f %0.3f %0.3f %0.3f" % (st.T, st.rhoM, st.P.kPa, st.hM.kJkmol, st.sM.kJkmolK, st.cvM.kJkmolK, st.cpM.kJkmolK, st.w)
            295 4 3180.336 34509.528 161.389 120.218 2189.730 73.317
            """, # Table 10, Pag 842

        "R": 8.314472,
        "cp": Fi1,
        "ref": "NBP",

        "Tmin": Tt, "Tmax": 425.0, "Pmax": 50000.0, "rhomax": 12.31,
        "Pmin": 26.1, "rhomin": 12.3,

        "nr1": [1.1632, -2.8123, 0.77202, -0.14331, 0.10227, 0.00024629],
        "d1": [1, 1, 1, 2, 3, 7],
        "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875],

        "nr2": [0.30893, -0.028499, -0.30343, -0.068793, -0.027218, 0.010665],
        "d2": [2, 5, 1, 4, 3, 4],
        "t2": [0.625, 1.75, 3.625, 3.625, 14.5, 12],
        "c2": [1, 1, 2, 2, 3, 3],
        "gamma2": [1]*6}

    helmholtz2 = {
        "__type__": "Helmholtz",
        "__name__": "Helmholtz equation of state for R-116 of Kozlov (1996).",
        "__doi__": {"autor": "Kozlov A.D.",
                    "title": "private communication with Dr. Alexander  D. Kozlov, Director, VNITs SMV Russian Research Center for Standartization Information and Certification of Materials, Nahimovsky prospect, 31, bld. 2 Moscow 117418, Russia. [email protected]",
                    "ref": "",
                    "doi":  ""},

        "R": 8.31451,
        "cp": CP2,

        "Tmin": Tt, "Tmax": 425.0, "Pmax": 50000.0, "rhomax": 12.23,
        "Pmin": 32.09, "rhomin": 12.231,
        "Pmin": 32.09, "rhomin": 12.231,

        "nr1": [2.1775273, -5.5052198, -1.3675742, -8.1284229e-1,
                -4.0207525e-1, 2.5890073, 1.4500537, -1.0445036, 9.8965288e-1,
                -8.6794888e-1, 2.8240917e-1, 4.5154220e-2, -3.0294024e-2,
                -1.7668398e-2, 2.0592774e-3],
        "d1": [1, 1, 1, 1, 2, 2, 2, 3, 3, 4, 5, 6, 6, 7, 8],
        "t1": [0.25, 1, 3, 4, 0.25, 1, 3.5, 1.5, 2.5, 3, 3, 1, 3, 1, 1],

        "nr2": [4.2059839, 2.1500380e-1, -1.6449561e-1, -1.2396086e-1,
                1.5814552e-1, -1.4362345e-1, 1.8637877e-2, 1.6342835e-2],
        "d2": [1, 1, 4, 4, 5, 5, 8, 4],
        "t2": [2, 5, 2, 4, 8, 10, 10, 18],
        "c2": [1, 2, 2, 2, 3, 3, 3, 4],
        "gamma2": [1]*8}

    eq = helmholtz1, helmholtz2

    _surface = {"sigma": [0.047593, -0.0073402], "exp": [1.2666, 1.9892]}
    _vapor_Pressure = {
        "eq": 5,
        "ao": [-0.73997e1, 0.22554e1, -0.23385e1, -0.35244e1, 0.40350],
        "exp": [1.0, 1.5, 2.2, 4.8, 6.2]}
    _liquid_Density = {
        "eq": 1,
        "ao": [0.68490e2, -0.24772e3, 0.35824e3, -0.25290e3, 0.76880e2],
        "exp": [0.64, 0.79, 0.95, 1.14, 1.33]}
    _vapor_Density = {
        "eq": 3,
        "ao": [-0.34135e1, -0.14529e3, 0.23651e3, -0.22276e3, 0.23103e3, -0.17433e3],
        "exp": [0.428, 2.0, 2.24, 3.0, 4.0, 5.0]}

    visco0 = {"eq": 1, "omega": 1,
              "collision": [0.355404, -0.464337, 0.257353e-1],
              "__name__": "Huber (2003)",
              "__doi__": {"autor": "Huber, M.L., Laesecke, A., and Perkins, R.A.",
                          "title": "Model for the Viscosity and Thermal Conductivity of Refrigerants, Including a New Correlation for the Viscosity of R134a",
                          "ref": "Ind. Eng. Chem. Res., 2003, 42 (13), pp 3163–3178",
                          "doi": "10.1021/ie0300880"},

              "ek": 226.16, "sigma": 0.5249,
              "Tref": 1., "rhoref": 1.*M,
              "n_chapman": 0.2509/M**0.5,

              "n_virial": [-0.19572881e2, 0.21973999e3, -0.10153226e4,
                           0.24710125e4, -0.33751717e4, 0.24916597e4,
                           -0.78726086e3, 0.14085455e2, -0.34664158],
              "t_virial": [0, -0.25, -0.5, -0.75, -1, -1.25, -1.5, -2.5, -5.5],
              "Tref_virial": 226.16, "etaref_virial": 0.08709,

              "Tref_res": 513.9, "rhoref_res": 5.991*M, "etaref_res": 1000,
              "n_packed": [3.16369563558749, -0.8901733752064137e-1,
                           0.1000352946668359],
              "t_packed": [0, 1, 2],
              "n_poly": [-0.2069007192080741e-1, 0.3560295489828222e-3,
                         0.2111018162451597e-2, 0.1396014148308975e-1,
                         -0.4564350196734897e-2, -0.3515932745836890e-2,
                         -0.2147633195397038],
              "t_poly": [0, -6, -2, -0.5, 2, 0, 0],
              "d_poly": [1, 2, 2, 2, 2, 3, 0],
              "g_poly": [0, 0, 0, 0, 0, 0, -1],
              "c_poly": [0, 0, 0, 0, 0, 0, 0],
              "n_num": [0.2147633195397038],
              "t_num": [0],
              "d_num": [0],
              "g_num": [0],
              "c_num": [0],
              "n_den": [1, -1],
              "t_den": [0, 0],
              "d_den": [0, 1],
              "g_den": [1, 0],
              "c_den": [0, 0]}

    _viscosity = visco0,

    thermo0 = {"eq": 1,
               "__name__": "Perkins (2000)",
               "__doi__": {"autor": "Perkins, R.A., Laesecke, A., Howley, J., Ramires, M.L.V., Gurova, A.N., and Cusco, L.",
                           "title": "Experimental thermal conductivity values for the IUPAC round-robin sample of 1,1,1,2-tetrafluoroethane (R134a)",
                           "ref": "NIST Interagency/Internal Report (NISTIR) - 6605",
                           "doi": ""},

               "Tref": 1, "kref": 1.05,
               "no": [-1.05248e-2, 8.00982e-5],
               "co": [0, 1],

               "Trefb": 1.0, "rhorefb": 4.444, "krefb": 1.64e-3,
               "nb": [1.836526, 5.126143, -1.436883, 6.261441e-1],
               "tb": [0, 0, 0, 0],
               "db": [1, 2, 3, 4],
               "cb": [0, 0, 0, 0],

               "critical": 3,
               "gnu": 0.63, "gamma": 1.239, "R0": 1.03,
               "Xio": 0.194e-9, "gam0": 0.0496e-1, "qd": 0.5e-9, "Tcref": 439.545}

    _thermal = thermo0,
Esempio n. 4
0
class Butene_1(MEoS):
    """Multiparameter equation of state for 1-butene"""
    name = "butene"
    CASNumber = "106-98-9"
    formula = "CH3-CH2-CH=CH2"
    synonym = ""
    _refPropName = "1BUTENE"
    _coolPropName = "1-Butene"
    rhoc = unidades.Density(237.8907968)
    Tc = unidades.Temperature(419.29)
    Pc = unidades.Pressure(4005.1, "kPa")
    M = 56.10632  # g/mol
    Tt = unidades.Temperature(87.8)
    Tb = unidades.Temperature(266.84)
    f_acent = 0.192
    momentoDipolar = unidades.DipoleMoment(0.339, "Debye")
    id = 24

    Fi1 = {"ao_log": [1, 2.9197],
           "pow": [0, 1],
           "ao_pow": [-0.00101126, 2.3869174],
           "ao_exp": [2.9406, 6.5395, 14.535, 5.8971],
           "titao": [274/Tc, 951/Tc, 2127/Tc, 5752/Tc]}

    lemmon = {
        "__type__": "Helmholtz",
        "__name__": "short Helmholtz equation of state for 1-butene of Lemmon "
                    "and Ihmels (2005)",
        "__doi__": {"autor": "Lemmon, E.W., Ihmels, E.C.",
                    "title": "Thermodynamic properties of the butenes: Part "
                             "II. Short fundamental equations of state",
                    "ref": "Fluid Phase Equilibria 228-229 (2005) 173-187",
                    "doi":  "10.1016/j.fluid.2004.09.004"},

        "R": 8.314472,
        "cp": Fi1,
        "ref": "NBP",

        "Tmin": Tt, "Tmax": 525., "Pmax": 70000.0, "rhomax": 14.59,

        "nr1": [0.78084, -2.8258, 0.99403, 0.017951, 0.088889, 0.00024673],
        "d1": [1, 1, 1, 2, 3, 7],
        "t1": [0.12, 1.3, 1.74, 2.1, 0.28, 0.69],

        "nr2": [0.22846, -0.074009, -0.22913, -0.062334, -0.025385, 0.011040],
        "d2": [2, 5, 1, 4, 3, 4],
        "t2": [0.75, 2., 4.4, 4.7, 15., 14.],
        "c2": [1, 1, 2, 2, 3, 3],
        "gamma2": [1]*6}

    eq = lemmon,
    _PR = [-0.2050, -16.0087]

    _surface = {"sigma": [0.05644], "exp": [1.248]}
    _vapor_Pressure = {
        "eq": 3,
        "n": [-0.71727e1, 0.26360e1, -0.20781e1, -0.28860e1, -0.13041e1],
        "t": [1, 1.5, 2, 4.35, 16.]}
    _liquid_Density = {
        "eq": 1,
        "n": [0.16857e2, -0.46280e2, 0.53727e2, -0.23314e2, 0.18889e1],
        "t": [0.547, 0.73, 0.92, 1.14, 2.1]}
    _vapor_Density = {
        "eq": 2,
        "n": [-3.1106, -6.3103, -19.272, -48.739, -99.898, -190.01],
        "t": [0.415, 1.27, 3.34, 7.0, 14.5, 28.0]}
Esempio n. 5
0
class R218(MEoS):
    """Multiparameter equation of R218"""
    name = "octafluoropropane"
    CASNumber = "76-19-7"
    formula = "CF3CF2CF3"
    synonym = "R218"
    _refPropName = "R218"
    _coolPropName = "R218"
    rhoc = unidades.Density(627.9845622)
    Tc = unidades.Temperature(345.02)
    Pc = unidades.Pressure(2640.0, "kPa")
    M = 188.01933  # g/mol
    Tt = unidades.Temperature(125.45)
    Tb = unidades.Temperature(236.36)
    f_acent = 0.3172
    momentoDipolar = unidades.DipoleMoment(0.14, "Debye")
    id = 671

    Fi1 = {
        "ao_log": [1, 3.],
        "pow": [0, 1],
        "ao_pow": [-15.6587335175, 11.4531412796],
        "ao_exp": [7.2198, 7.2692, 11.599],
        "titao": [326 / Tc, 595 / Tc, 1489 / Tc]
    }

    CP1 = {
        "ao": 4.,
        "an": [],
        "pow": [],
        "ao_exp": [7.2198, 7.2692, 11.599],
        "exp": [326, 595, 1489]
    }

    lemmon = {
        "__type__": "Helmholtz",
        "__name__": "short Helmholtz equation of state for R-218 of Lemmon "
        "and Span (2006)",
        "__doi__": {
            "autor": "Lemmon, E.W., Span, R.",
            "title": "Short Fundamental Equations of State for 20 "
            "Industrial Fluids",
            "ref": "J. Chem. Eng. Data, 2006, 51 (3), pp 785–850",
            "doi": "10.1021/je050186n"
        },
        "R": 8.314472,
        "cp": Fi1,
        "ref": "NBP",
        "Tmin": Tt,
        "Tmax": 440.0,
        "Pmax": 20000.0,
        "rhomax": 10.69,
        "nr1": [1.3270, -3.8433, 0.922, 0.1136, 0.00036195],
        "d1": [1, 1, 1, 3, 7],
        "t1": [0.25, 1.25, 1.5, 0.25, 0.875],
        "nr2": [1.1001, 1.1896, -.025147, -.65923, -.027969, -.1833, -.02163],
        "d2": [1, 2, 5, 1, 1, 4, 2],
        "t2": [2.375, 2., 2.125, 3.5, 6.5, 4.75, 12.5],
        "c2": [1, 1, 1, 2, 2, 2, 3],
        "gamma2": [1] * 7
    }

    eq = lemmon,
    _PR = [-0.3225, -15.1811]

    _surface = {"sigma": [0.04322], "exp": [1.224]}
    _vapor_Pressure = {
        "eq": 3,
        "n": [-0.78419e1, 0.28989e1, -0.33458e1, -0.33196e1, 0.25363],
        "t": [1.0, 1.5, 2.2, 4.8, 6.2]
    }
    _liquid_Density = {
        "eq": 1,
        "n": [-0.61027, 0.57453e1, -0.56835e1, 0.32137e1, 0.55194],
        "t": [0.223, 0.39, 0.56, 0.75, 5.0]
    }
    _vapor_Density = {
        "eq": 2,
        "n": [-4.2658, -6.9496, -.18099e2, -.4921e2, -.55945e2, -.74492e2],
        "t": [0.481, 1.53, 3.2, 6.3, 12.0, 15.0]
    }

    trnECS = {
        "__name__": "Huber (2003)",
        "__doi__": {
            "autor":
            "Huber, M.L., Laesecke, A., Perkins, R.A.",
            "title":
            "Model for the Viscosity and Thermal Conductivity "
            "of Refrigerants, Including a New Correlation for "
            "the Viscosity of R134a",
            "ref":
            "Ind. Eng. Chem. Res., 42(13) (2003) 3163-3178",
            "doi":
            "10.1021/ie0300880"
        },
        "eq": "ecs",
        "ref": C3,
        "visco": "visco1",
        "thermo": "thermo0",
        "ek": 266.35,
        "sigma": 0.58,
        "omega": 5,
        "psi": [1.10225, -5.50442e-3],
        "psi_d": [0, 1],
        "fint": [8.92659e-4, 1.14912e-6],
        "fint_t": [0, 1],
        "chi": [1.2877, -7.58811e-2],
        "chi_d": [0, 1],
        "critical": 3,
        "gnu": 0.63,
        "gamma": 1.239,
        "R0": 1.03,
        "Xio": 0.194e-9,
        "gam0": 0.0496,
        "qd": 1.09043e-9,
        "Tcref": 1.5 * Tc
    }

    _viscosity = trnECS,
    _thermal = trnECS,
Esempio n. 6
0
class R114(MEoS):
    """Multiparameter equation of state for R114"""
    name = "1,2-dichloro-1,1,2,2-tetrafluoroethane"
    CASNumber = "76-14-2"
    formula = "CClF2CClF2"
    synonym = "R114"
    rhoc = unidades.Density(579.969)
    Tc = unidades.Temperature(418.83)
    Pc = unidades.Pressure(3257.0, "kPa")
    M = 170.921  # g/mol
    Tt = unidades.Temperature(180.63)
    Tb = unidades.Temperature(276.741)
    f_acent = 0.25253
    momentoDipolar = unidades.DipoleMoment(0.658, "Debye")
    id = 231

    CP1 = {
        "ao": 0.97651380e-1,
        "an": [0.32408610e-2, -0.58953640e-5, 0.67379290e-8, -0.35463640e-11],
        "pow": [1, 2, 3, 4],
        "ao_exp": [],
        "exp": [],
        "ao_hyp": [],
        "hyp": []
    }

    helmholtz1 = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Bender equation of state for R-114 of Platzer et al. (1990).",
        "__doi__": {
            "autor": "Platzer, B., Polt, A., and Maurer, G.",
            "title": "Thermophysical properties of refrigerants",
            "ref": "Berlin:  Springer-Verlag, 1990.",
            "doi": ""
        },
        "R":
        8.31451,
        "cp":
        CP1,
        "ref":
        "NBP",
        "Tmin":
        273.15,
        "Tmax":
        507.0,
        "Pmax":
        21000.0,
        "rhomax":
        8.942,
        "Pmin":
        0.2,
        "rhomin":
        10.4,
        "nr1": [
            -0.340776521414, 0.323001398420, -0.424950537596e-1,
            0.107938879710e1, -0.199243619673e1, -0.155135133506,
            -0.121465790553, -0.165038582393e-1, -0.186915808643,
            0.308074612567, 0.115861416115, 0.276358316589e-1, 0.108043243088,
            0.460683793064e-1, -0.174821616881, 0.317530854287e-1
        ],
        "d1": [0, 0, 0, 1, 1, 1, 1, 1, 2, 2, 2, 3, 3, 4, 4, 5],
        "t1": [3, 4, 5, 0, 1, 2, 3, 4, 0, 1, 2, 0, 1, 0, 1, 1],
        "nr2": [
            0.340776521414, -0.323001398420, 0.424950537596e-1,
            -0.166940100976e1, 0.408693082002e1, -0.241738963889e1
        ],
        "d2": [0, 0, 0, 2, 2, 2],
        "t2": [3, 4, 5, 3, 4, 5],
        "c2": [2] * 6,
        "gamma2": [1.21103865] * 6
    }

    eq = helmholtz1,

    _surface = {"sigma": [0.05239], "exp": [1.258]}
    _vapor_Pressure = {
        "eq": 5,
        "ao": [-0.72195e1, 0.16357e1, -0.14576e1, -0.69580e1, 0.57181e1],
        "exp": [1.0, 1.5, 2.2, 4.8, 6.2]
    }
    _liquid_Density = {
        "eq": 1,
        "ao": [0.43023, 0.22722e2, -0.27118e2, 0.13247e2, -0.90529e1],
        "exp": [0.095, 0.93, 1.1, 2.0, 3.0]
    }
    _vapor_Density = {
        "eq": 3,
        "ao":
        [-0.46609, -0.68355e1, -0.16715e3, 0.15805e5, -0.31859e5, 0.21548e5],
        "exp": [0.09, 0.76, 4.0, 6.5, 7.0, 8.0]
    }
Esempio n. 7
0
class iButene(MEoS):
    """Multiparameter equation of state for isobutene"""
    name = "isobutene"
    CASNumber = "115-11-7"
    formula = "CH2=C(CH3)2"
    synonym = ""
    _refPropName = "IBUTENE"
    _coolPropName = "IsoButene"
    rhoc = unidades.Density(233.9633544)
    Tc = unidades.Temperature(418.09)
    Pc = unidades.Pressure(4009.8, "kPa")
    M = 56.10632  # g/mol
    Tt = unidades.Temperature(132.4)
    Tb = unidades.Temperature(266.15)
    f_acent = 0.193
    momentoDipolar = unidades.DipoleMoment(0.5, "Debye")
    id = 27

    Fi1 = {
        "ao_log": [1, 3.],
        "pow": [0, 1],
        "ao_pow": [-0.12737888, 2.3125128],
        "ao_exp": [4.8924, 7.832, 7.2867, 8.7293],
        "titao": [399 / Tc, 1270 / Tc, 2005 / Tc, 4017 / Tc]
    }

    lemmon = {
        "__type__": "Helmholtz",
        "__name__": "short Helmholtz equation of state for 1-butene of Lemmon "
        "and Ihmels (2005)",
        "__doi__": {
            "autor": "Lemmon, E.W., Ihmels, E.C.",
            "title": "Thermodynamic properties of the butenes: Part "
            "II. Short fundamental equations of state",
            "ref": "Fluid Phase Equilibria 228-229 (2005) 173-187",
            "doi": "10.1016/j.fluid.2004.09.004"
        },
        "R": 8.314472,
        "cp": Fi1,
        "ref": "NBP",
        "Tmin": Tt,
        "Tmax": 550.0,
        "Pmax": 50000.0,
        "rhomax": 13.67,
        "nr1": [0.77111, -2.7971, 1.0118, 0.02073, 0.085086, 0.00021968],
        "d1": [1, 1, 1, 2, 3, 7],
        "t1": [0.12, 1.3, 1.74, 2.1, 0.28, 0.69],
        "nr2": [0.20633, -0.078843, -0.23726, -0.080211, -0.027001, 0.013072],
        "d2": [2, 5, 1, 4, 3, 4],
        "t2": [0.75, 2., 4.4, 4.7, 15., 14.],
        "c2": [1, 1, 2, 2, 3, 3],
        "gamma2": [1] * 6
    }

    eq = lemmon,

    _surface = {"sigma": [0.0545], "exp": [1.23]}
    _vapor_Pressure = {
        "eq": 3,
        "n": [-0.68973e1, 0.12475e1, -0.25441e1, -0.29282e1, 0.15778e1],
        "t": [1., 1.5, 3.16, 6.2, 7.0]
    }
    _liquid_Density = {
        "eq": 1,
        "n": [0.62591e2, -0.20805e3, 0.33243e3, -0.29555e3, 0.11148e3],
        "t": [0.65, 0.8, 0.98, 1.16, 1.3]
    }
    _vapor_Density = {
        "eq": 2,
        "n": [-3.1841, -6.4014, -9.3817, -0.11160e2, -0.52298e2, -0.12195e3],
        "t": [0.431, 1.29, 3.3, 3.54, 7.3, 15.8]
    }
Esempio n. 8
0
class MD3M(MEoS):
    """Multiparameter equation of state for dodecamethylpentasiloxane"""
    name = "dodecamethylpentasiloxane"
    CASNumber = "141-63-9"
    formula = "C12H36Si5O4"
    synonym = "MD3M"
    _refPropName = "MD3M"
    _coolPropName = "MD3M"
    rhoc = unidades.Density(263.9218791237794)
    Tc = unidades.Temperature(628.36)
    Pc = unidades.Pressure(945.0, "kPa")
    M = 384.839  # g/mol
    Tt = unidades.Temperature(192.0)
    Tb = unidades.Temperature(503.03)
    f_acent = 0.722
    momentoDipolar = unidades.DipoleMoment(1.223, "Debye")

    CP1 = {
        "ao": 463.2,
        "an": [],
        "pow": [],
        "ao_exp": [],
        "exp": [],
        "ao_hyp": [609372332.2, 0, 4290277999.0, 0],
        "hyp": [908.5, 0, 2117.1, 0]
    }

    colonna = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for MD3M of Colonna (2006).",
        "__doi__": {
            "autor":
            "Colonna, P., Nannan, N.R., and Guardone, A.",
            "title":
            "Multiparameter equations of state for siloxanes:"
            " [(CH3)3-Si-O1/2]2-[O-Si-(CH3)2]i=1,…,3, and "
            "[O-Si-(CH3)2]6",
            "ref":
            "Fluid Phase Equilibria 263:115-130, 2008",
            "doi":
            "10.1016/j.fluid.2007.10.001"
        },
        "R":
        8.314472,
        "cp":
        CP1,
        "ref":
        "NBP",
        "Tmin":
        Tt,
        "Tmax":
        673.0,
        "Pmax":
        30000.0,
        "rhomax":
        2.54,
        "Pmin":
        0.4e-12,
        "rhomin":
        2.54,
        "nr1": [
            1.20540386, -2.42914797, 0.69016432, -0.69268041, 0.18506046,
            0.31161436e-3
        ],
        "d1": [1, 1, 1, 2, 3, 7],
        "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875],
        "nr2": [
            0.99862519, 0.74229034e-1, -0.80259136, -0.20865337,
            -0.36461791e-1, 0.19174051e-1
        ],
        "d2": [2, 5, 1, 4, 3, 4],
        "t2": [0.625, 1.75, 3.625, 3.625, 14.5, 12.0],
        "c2": [1, 1, 2, 2, 3, 3],
        "gamma2": [1] * 6
    }

    eq = colonna,

    _vapor_Pressure = {
        "eq": 3,
        "n": [-0.92608e1, 0.15861e1, -0.32859e1, -0.75194e1, -0.34883e1],
        "t": [1.0, 1.5, 2.46, 3.7, 10.0]
    }
    _liquid_Density = {
        "eq": 1,
        "n": [0.74156, 0.21723e1, 0.66412e2, -0.17125e3, 0.10848e3],
        "t": [0.22, 0.51, 5.5, 6.0, 6.4]
    }
    _vapor_Density = {
        "eq":
        2,
        "n": [
            -0.19054e1, -0.74526e1, -0.10520e3, 0.24548e3, -0.23783e3,
            -0.21226e3
        ],
        "t": [0.332, 0.88, 3.25, 4.0, 4.6, 12.0]
    }
Esempio n. 9
0
class Trans_2_butene(MEoS):
    """Multiparameter equations of state for trans-butene"""
    name = "trans-butene"
    CASNumber = "624-64-6"
    formula = "CH3-CH=CH-CH3"
    synonym = ""
    _refPropName = "T2BUTENE"
    _coolPropName = "trans-2-Butene"
    rhoc = unidades.Density(236.37592616)
    Tc = unidades.Temperature(428.61)
    Pc = unidades.Pressure(4027.3, "kPa")
    M = 56.10632  # g/mol
    Tt = unidades.Temperature(167.6)
    Tb = unidades.Temperature(274.03)
    f_acent = 0.21
    momentoDipolar = unidades.DipoleMoment(0.0, "Debye")
    id = 26

    Fi1 = {
        "ao_log": [1, 2.9988],
        "pow": [0, 1],
        "ao_pow": [0.5917816, 2.1427758],
        "ao_exp": [5.3276, 13.29, 9.6745, 0.40087],
        "titao": [362 / Tc, 1603 / Tc, 3729 / Tc, 4527 / Tc]
    }

    lemmon = {
        "__type__": "Helmholtz",
        "__name__": "short Helmholtz equation of state for 1-butene of Lemmon "
        "and Ihmels (2005)",
        "__doi__": {
            "autor": "Lemmon, E.W., Ihmels, E.C.",
            "title": "Thermodynamic properties of the butenes: Part "
            "II. Short fundamental equations of state",
            "ref": "Fluid Phase Equilibria 228-229 (2005) 173-187",
            "doi": "10.1016/j.fluid.2004.09.004"
        },
        "R": 8.314472,
        "cp": Fi1,
        "ref": "NBP",
        "Tmin": Tt,
        "Tmax": 525.0,
        "Pmax": 50000.0,
        "rhomax": 13.141,
        "nr1": [0.81107, -2.8846, 1.0265, 0.016591, 0.086511, 0.00023256],
        "d1": [1, 1, 1, 2, 3, 7],
        "t1": [0.12, 1.3, 1.74, 2.1, 0.28, 0.69],
        "nr2": [0.22654, -0.072182, -0.24849, -0.071374, -0.024737, 0.011843],
        "d2": [2, 5, 1, 4, 3, 4],
        "t2": [0.75, 2., 4.4, 4.7, 15., 14.],
        "c2": [1, 1, 2, 2, 3, 3],
        "gamma2": [1] * 6
    }

    eq = lemmon,

    _vapor_Pressure = {
        "eq": 3,
        "n": [-7.6226, 7.9421, -6.9631, -6.5517, 3.9584],
        "t": [1.0, 1.5, 1.65, 4.8, 5.3]
    }
    _liquid_Density = {
        "eq": 1,
        "n": [12.452, -34.419, 52.257, -42.889, 15.463],
        "t": [0.52, 0.73, 0.97, 1.24, 1.5]
    }
    _vapor_Density = {
        "eq": 2,
        "n": [-3.1276, -6.0548, -18.243, -60.842, 135.95, -182.70],
        "t": [0.412, 1.24, 3.2, 7.0, 10.0, 11.0]
    }
Esempio n. 10
0
class nC10(MEoS):
    """Multiparameter equation of state for n-decane"""
    name = "decane"
    CASNumber = "124-18-5"
    formula = "CH3-(CH2)8-CH3"
    synonym = ""
    _refPropName = "DECANE"
    _coolPropName = "n-Decane"
    rhoc = unidades.Density(233.342)
    Tc = unidades.Temperature(617.7)
    Pc = unidades.Pressure(2103.0, "kPa")
    M = 142.28168  # g/mol
    Tt = unidades.Temperature(243.5)
    Tb = unidades.Temperature(447.27)
    f_acent = 0.4884
    momentoDipolar = unidades.DipoleMoment(0.07, "Debye")
    id = 14

    Fi1 = {
        "ao_log": [1, 18.109],
        "pow": [0, 1],
        "ao_pow": [13.9361966549, -10.5265128286],
        "ao_exp": [25.685, 28.233, 12.417, 10.035],
        "titao": [1193 / Tc, 2140 / Tc, 4763 / Tc, 10862 / Tc]
    }

    Fi2 = {
        "ao_log": [1, 3.0],
        "pow": [0, 1],
        "ao_pow": [15.870791919, -108.858547525],
        "ao_sinh": [21.0069, 58.3657],
        "sinh": [0.267034159, 2.833479035],
        "ao_cosh": [43.4931],
        "cosh": [1.353835195]
    }

    lemmon = {
        "__type__": "Helmholtz",
        "__name__": "short Helmholtz equation of state for decane of Lemmon "
        "and Span (2006)",
        "__doi__": {
            "autor": "Lemmon, E.W., Span, R.",
            "title": "Short Fundamental Equations of State for 20 "
            "Industrial Fluids",
            "ref": "J. Chem. Eng. Data, 2006, 51 (3), pp 785–850",
            "doi": "10.1021/je050186n"
        },
        "R": 8.314472,
        "cp": Fi1,
        "ref": "NBP",
        "Tmin": Tt,
        "Tmax": 675.0,
        "Pmax": 800000.0,
        "rhomax": 5.41,
        "nr1": [1.0461, -2.4807, 0.74372, -0.52579, 0.15315, 0.00032865],
        "d1": [1, 1, 1, 2, 3, 7],
        "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875],
        "nr2": [0.84178, 0.055424, -0.73555, -0.18507, -0.020775, 0.012335],
        "d2": [2, 5, 1, 4, 3, 4],
        "t2": [0.625, 1.75, 3.625, 3.625, 14.5, 12.],
        "c2": [1, 1, 2, 2, 3, 3],
        "gamma2": [1] * 6
    }

    GERG = {
        "__type__": "Helmholtz",
        "__name__": "Helmholtz equation of state for propane of Kunz and "
        "Wagner (2008).",
        "__doi__": {
            "autor":
            "Kunz, O., Wagner, W.",
            "title":
            "The GERG-2008 Wide-Range Equation of State for "
            "Natural Gases and Other Mixtures: An Expansion "
            "of GERG-2004",
            "ref":
            "J. Chem.Eng. Data 57(11) (2012) 3032-3091",
            "doi":
            "10.1021/je300655b"
        },
        "R": 8.314472,
        "cp": Fi2,
        "ref": "OTO",
        "Tmin": Tt,
        "Tmax": 675.0,
        "Pmax": 800000.0,
        "rhomax": 5.41,
        "nr1": [0.10461e1, -0.24807e1, 0.74372, -0.52579, 0.15315, 0.32865e-3],
        "d1": [1, 1, 1, 2, 3, 7],
        "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875],
        "nr2": [.84178, .55424e-1, -.73555, -.18507, -.20775e-1, .12335e-1],
        "d2": [2, 5, 1, 4, 3, 4],
        "t2": [0.625, 1.75, 3.625, 3.625, 14.5, 12.],
        "c2": [1, 1, 2, 2, 3, 3],
        "gamma2": [1] * 6
    }

    eq = lemmon, GERG
    _PR = [0.0838, -23.0757]

    _surface = {"sigma": [0.05473], "exp": [1.29]}
    _dielectric = {
        "eq": 1,
        "a": [49.32, 0.050],
        "b": [220.15, -316.3],
        "c": [-88358, 53511],
        "Au": 29.84,
        "D": 2
    }

    _vapor_Pressure = {
        "eq": 3,
        "n": [-0.87738e1, 0.40864e1, -0.40775e1, -0.64910e1, 0.15598e1],
        "t": [1.0, 1.5, 1.93, 4.14, 4.7]
    }
    _liquid_Density = {
        "eq": 1,
        "n": [0.92435e1, -0.16288e2, 0.20445e2, -0.17624e2, 0.73796e1],
        "t": [0.535, 0.74, 1.0, 1.28, 1.57]
    }
    _vapor_Density = {
        "eq": 2,
        "n": [-5.0378, -3.4694, -.15906e2, -0.82894e2, 0.29336e2, -0.10985e3],
        "t": [0.4985, 1.33, 2.43, 5.44, 5.8, 11.0]
    }

    visco0 = {
        "__name__":
        "Huber (2004)",
        "__doi__": {
            "autor":
            "Huber, M.L., Laesecke, A. Xiang, H.W.",
            "title":
            "Viscosity correlations for minor constituent "
            "fluids in natural gas: n-octane, n-nonane and "
            "n-decane",
            "ref":
            "Fluid Phase Equilibria 224 (2004) 263-270",
            "doi":
            "10.1016/j.fluid.2004.07.012"
        },
        "eq":
        1,
        "omega":
        1,
        "ek":
        490.51,
        "sigma":
        0.686,
        "n_chapman":
        0.021357,
        "collision": [0.343267, -0.460514],
        "Tref_virial":
        490.51,
        "n_virial": [
            -19.572881, 219.73999, -1015.3226, 2471.0125, -3375.1717,
            2491.6597, -787.26086, 14.085455, -0.34664158
        ],
        "t_virial": [0, -0.25, -0.5, -0.75, -1, -1.25, -1.5, -2.5, -5.5],
        "Tref_res":
        617.7,
        "rhoref_res":
        1.64 * M,
        "muref_res":
        1000,
        "nr": [-.0402094, 0.0404435, -0.0142063],
        "tr": [1, 2, 2],
        "dr": [2, 2, 3],
        "CPf":
        453.387,
        "CPg1":
        2.55105,
        "CPgi": [1.71465 / 2.55105],
        "CPti": [-0.5]
    }

    _viscosity = visco0,

    thermo0 = {
        "__name__":
        "Huber (2005)",
        "__doi__": {
            "autor":
            "Huber, M.L., Perkins, R.A.",
            "title":
            "Thermal conductivity correlations for minor "
            "constituent fluids in natural gas: n-octane, "
            "n-nonane and n-decane",
            "ref":
            "Fluid Phase Equilibria 227 (2005) 47-55",
            "doi":
            "10.1016/j.fluid.2004.10.031"
        },
        "eq":
        1,
        "Toref":
        617.7,
        "koref":
        1,
        "no": [0.105543e-1, -0.514530e-1, 0.118979, -0.372442e-1],
        "to": [0, 1, 2, 3],
        "Tref_res":
        617.7,
        "rhoref_res":
        1.64 * M,
        "kref_res":
        1,
        "nr": [
            -.294394e-1, .150509e-1, .499245e-1, -.142700e-1, -0.138857e-1,
            0.150828e-2, 0.433326e-2
        ],
        "tr": [0, -1, 0, 0, -1, 0, -1],
        "dr": [1, 1, 2, 3, 3, 4, 4],
        "critical":
        3,
        "gnu":
        0.63,
        "gamma":
        1.239,
        "R0":
        1.03,
        "Xio":
        0.194e-9,
        "gam0":
        0.0496,
        "qd":
        7.086368e-10,
        "Tcref":
        926.55
    }

    _thermal = thermo0,
Esempio n. 11
0
class RC318(MEoS):
    """Multiparameter equation of state for RC318"""
    name = "octafluorocyclobutane "
    CASNumber = "406-58-6"
    formula = "cyclo-C4F8"
    synonym = "RC318"
    rhoc = unidades.Density(619.973)
    Tc = unidades.Temperature(388.38)
    Pc = unidades.Pressure(2777.5, "kPa")
    M = 200.0312  # g/mol
    Tt = unidades.Temperature(233.35)
    Tb = unidades.Temperature(267.175)
    f_acent = 0.3553
    momentoDipolar = unidades.DipoleMoment(0.0, "Debye")
    id = 692

    CP1 = {
        "ao": 0.121,
        "an": [0.2903e-2, -0.25327e-5, 0.77191e-9],
        "pow": [1, 2, 3],
        "ao_exp": [],
        "exp": [],
        "ao_hyp": [],
        "hyp": []
    }

    helmholtz1 = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for R-C318 of Platzer et al. (1990)",
        "__doi__": {
            "autor": "Platzer, B., Polt, A., and Maurer, G.",
            "title": "Thermophysical properties of refrigerants",
            "ref": "Berlin:  Springer-Verlag, 1990.",
            "doi": ""
        },
        "R":
        8.31451,
        "cp":
        CP1,
        "ref":
        "NBP",
        "Tmin":
        Tt,
        "Tmax":
        623.0,
        "Pmax":
        60000.0,
        "rhomax":
        8.6452,
        "Pmin":
        19.461,
        "rhomin":
        8.6452,
        "nr1": [
            -0.104729119796e1, 0.138034128822e1, -0.333625769594,
            0.109415569278e1, -0.268265237178e1, 0.173403063905e1,
            -0.163611246876e1, 0.304834499143, 0.102771551787,
            -0.232367895587e-1, 0.166151957803, -0.250103914479e-1,
            0.935680977639e-1, 0.431929127445e-1, -0.133439845861,
            0.255416632165e-1
        ],
        "d1": [0, 0, 0, 1, 1, 1, 1, 1, 2, 2, 2, 3, 3, 4, 4, 5],
        "t1": [3, 4, 5, 0, 1, 2, 3, 4, 0, 1, 2, 0, 1, 0, 1, 1],
        "nr2": [
            0.104729119796e1, -0.138034128822e1, 0.333625769594,
            -0.510485781618, 0.181840728111e1, -0.138530893970e1
        ],
        "d2": [0, 0, 0, 2, 2, 2],
        "t2": [3, 4, 5, 3, 4, 5],
        "c2": [2] * 6,
        "gamma2": [0.99943992] * 6
    }

    eq = helmholtz1,

    _surface = {"sigma": [0.0507], "exp": [1.25]}
    _vapor_Pressure = {
        "eq": 5,
        "ao": [-0.78467e1, 0.24555e1, -0.30824e1, -0.58263e1, 0.35483e1],
        "exp": [1.0, 1.5, 2.2, 4.8, 6.2]
    }
    _liquid_Density = {
        "eq": 1,
        "ao": [-0.30181, 0.29345e1, -0.13741e1, 0.14650e1, 0.16963],
        "exp": [0.11, 0.32, 0.57, 0.84, 2.9]
    }
    _vapor_Density = {
        "eq": 3,
        "ao": [-0.24491e2, 0.53255e2, -0.38863e2, -0.24938e2, -0.90092e2],
        "exp": [0.61, 0.77, 0.92, 3.3, 7.5]
    }
Esempio n. 12
0
class R113(MEoS):
    """Multiparameter equation of state for R113"""
    name = "1,1,2-trichloro-1,2,2-trifluoroethane"
    CASNumber = "76-13-1"
    formula = "CCl2FCClF2"
    synonym = "R113"
    rhoc = unidades.Density(560.)
    Tc = unidades.Temperature(487.21)
    Pc = unidades.Pressure(3392.2, "kPa")
    M = 187.375  # g/mol
    Tt = unidades.Temperature(236.93)
    Tb = unidades.Temperature(320.735)
    f_acent = 0.25253
    momentoDipolar = unidades.DipoleMoment(0.803, "Debye")
    id = 232

    CP1 = {
        "ao": 3.9999966,
        "an": [],
        "pow": [],
        "ao_exp": [12.4464495, 2.72181845, 0.692712415, 3.32248298],
        "exp": [5.1143280e2, 1.60676324e3, 4.20292102e3, 1.60618738e3],
        "ao_hyp": [],
        "hyp": []
    }

    helmholtz1 = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for R-113 of Marx et al. (1992).",
        "__doi__": {
            "autor": "Marx, V., Pruss, A., and Wagner, W.",
            "title":
            "Neue Zustandsgleichungen fuer R 12, R 22, R 11 und R 113. Beschreibung des thermodynamishchen Zustandsverhaltens bei Temperaturen bis 525 K und Druecken bis 200 MPa",
            "ref":
            "Duesseldorf: VDI Verlag, Series 19 (Waermetechnik/Kaeltetechnik), No. 57, 1992.",
            "doi": ""
        },
        "R":
        8.314471,
        "cp":
        CP1,
        "Tmin":
        Tt,
        "Tmax":
        525.0,
        "Pmax":
        200000.0,
        "rhomax":
        9.10,
        "Pmin":
        1.87,
        "rhomin":
        9.099,
        "nr1": [
            0.8432092286, -0.2019185967e1, 0.2920612996, 0.5323107661e-1,
            0.3214971931e-2, 0.4667858574e-4, -0.1227522799e-5
        ],
        "d1": [1, 1, 2, 3, 4, 8, 8],
        "t1": [0.5, 1.5, 1.5, -0.5, 2, 0, 3],
        "nr2": [
            0.8167288718, -0.1340790803e1, 0.4065752705, -0.1534754634,
            -0.2414435149e-1, -0.2113056197e-1, -0.3565436205e-1,
            0.1364654968e-2, -0.1251838755e-1, -0.1385761351e-2,
            0.7206335486e-3
        ],
        "d2": [3, 3, 3, 5, 1, 2, 2, 9, 3, 7, 8],
        "t2": [-0.5, 0, 2, 1.5, 6, 2, 10, 6, 18, 15, 33],
        "c2": [1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 4],
        "gamma2": [1] * 11
    }

    helmholtz2 = {
        "__type__":
        "Helmholtz",
        "__name__":
        "short Helmholtz equation of state for R-113 of Span and Wagner (2003).",
        "__doi__": {
            "autor": "Span, R., Wagner, W.",
            "title":
            "Equations of State for Technical Applications. III. Results for Polar Fluids",
            "ref": "Int. J. Thermophys., 24(1):111-162, 2003.",
            "doi": "10.1023/A:1022362231796"
        },
        "__test__":
        """
            >>> st=R113(T=700, rho=200, eq=1)
            >>> print "%0.4f %0.3f %0.4f" % (st.cp0.kJkgK, st.P.MPa, st.cp.kJkgK)
            0.8055 3.962 4.0344
            >>> st2=R113(T=750, rho=100, eq=1)
            >>> print "%0.2f %0.5f" % (st2.h.kJkg-st.h.kJkg, st2.s.kJkgK-st.s.kJkgK)
            131.00 0.26004
            """,  # Table III, Pag 117
        "R":
        8.31451,
        "cp":
        CP1,
        "Tmin":
        Tt,
        "Tmax":
        600.0,
        "Pmax":
        100000.0,
        "rhomax":
        9.09,
        "Pmin":
        1.869,
        "rhomin":
        9.0893,
        "nr1": [
            0.10519071e1, -0.28724742e1, 0.41983153, 0.87107788e-1,
            0.24105194e-3
        ],
        "d1": [1, 1, 1, 3, 7],
        "t1": [0.25, 1.25, 1.5, 0.25, 0.875],
        "nr2": [
            0.70738262, 0.93513411, -0.96713512e-2, -0.52595315, 0.22691984e-1,
            -0.14556325, -0.2741995e-1
        ],
        "d2": [1, 2, 5, 1, 1, 4, 2],
        "t2": [2.375, 2, 2.125, 3.5, 6.5, 4.75, 12.5],
        "c2": [1, 1, 1, 2, 2, 2, 3],
        "gamma2": [1] * 7
    }

    eq = helmholtz1, helmholtz2

    _surface = {"sigma": [0.0556], "exp": [1.24]}
    _vapor_Pressure = {
        "eq": 5,
        "ao": [-0.73838e1, 0.32594e1, -0.27761e1, -0.37758e1, -0.19921],
        "exp": [1.0, 1.5, 1.8, 4.3, 6.2]
    }
    _liquid_Density = {
        "eq": 1,
        "ao": [0.15785e1, 0.12404e1, -0.66933, 0.49775e1, -0.55253e1],
        "exp": [0.3, 0.7, 2.0, 4.0, 5.0]
    }
    _vapor_Density = {
        "eq":
        3,
        "ao": [
            -0.26225e1, -0.60753e1, -0.15768e2, -0.42361e2, -0.79071e1,
            -0.31966e3
        ],
        "exp": [0.379, 1.13, 2.9, 6.0, 7.0, 15.0]
    }

    visco0 = {
        "eq":
        1,
        "omega":
        1,
        "collision": [0.355404, -0.464337, 0.257353e-1],
        "__name__":
        "Huber (2003)",
        "__doi__": {
            "autor": "Huber, M.L., Laesecke, A., and Perkins, R.A.",
            "title":
            "Model for the Viscosity and Thermal Conductivity of Refrigerants, Including a New Correlation for the Viscosity of R134a",
            "ref": "Ind. Eng. Chem. Res., 2003, 42 (13), pp 3163–3178",
            "doi": "10.1021/ie0300880"
        },
        "ek":
        376.035,
        "sigma":
        0.6019,
        "Tref":
        1.,
        "rhoref":
        1. * M,
        "n_chapman":
        0.2509 / M**0.5,
        "n_virial": [
            -0.19572881e2, 0.21973999e3, -0.10153226e4, 0.24710125e4,
            -0.33751717e4, 0.24916597e4, -0.78726086e3, 0.14085455e2,
            -0.34664158
        ],
        "t_virial": [0, -0.25, -0.5, -0.75, -1, -1.25, -1.5, -2.5, -5.5],
        "Tref_virial":
        376.035,
        "etaref_virial":
        0.13132,
        "Tref_res":
        487.21,
        "rhoref_res":
        2.988659 * M,
        "etaref_res":
        1310,
        "n_packed":
        [3.16369563558749, -0.8901733752064137e-1, 0.1000352946668359],
        "t_packed": [0, 1, 2],
        "n_poly": [
            -0.2069007192080741e-1, 0.3560295489828222e-3,
            0.2111018162451597e-2, 0.1396014148308975e-1,
            -0.4564350196734897e-2, -0.3515932745836890e-2, -0.2147633195397038
        ],
        "t_poly": [0, -6, -2, -0.5, 2, 0, 0],
        "d_poly": [1, 2, 2, 2, 2, 3, 0],
        "g_poly": [0, 0, 0, 0, 0, 0, -1],
        "c_poly": [0, 0, 0, 0, 0, 0, 0],
        "n_num": [0.2147633195397038],
        "t_num": [0],
        "d_num": [0],
        "g_num": [0],
        "c_num": [0],
        "n_den": [1, -1],
        "t_den": [0, 0],
        "d_den": [0, 1],
        "g_den": [1, 0],
        "c_den": [0, 0]
    }

    _viscosity = visco0,

    thermo0 = {
        "eq":
        1,
        "__name__":
        "Perkins (2000)",
        "__doi__": {
            "autor":
            "Perkins, R.A., Laesecke, A., Howley, J., Ramires, M.L.V., Gurova, A.N., and Cusco, L.",
            "title":
            "Experimental thermal conductivity values for the IUPAC round-robin sample of 1,1,1,2-tetrafluoroethane (R134a)",
            "ref": "NIST Interagency/Internal Report (NISTIR) - 6605",
            "doi": ""
        },
        "Tref":
        487.21,
        "kref":
        1.1,
        "no": [-0.460820e-2, 0.168688e-1, 0.488345e-2],
        "co": [0, 1, 2],
        "Trefb":
        487.21,
        "rhorefb":
        2.988659,
        "krefb":
        0.66,
        "nb": [
            -0.729410e-2, 0.110497e-1, 0.416339e-1, -0.289236e-1, -0.311487e-1,
            0.278399e-1, 0.112682e-1, -0.121100e-1, -0.138322e-2, 0.211196e-2
        ],
        "tb": [0, 1, 0, 1, 0, 1, 0, 1, 0, 1],
        "db": [1, 1, 2, 2, 3, 3, 4, 4, 5, 5],
        "cb": [0] * 10,
        "critical":
        3,
        "gnu":
        0.63,
        "gamma":
        1.239,
        "R0":
        1.03,
        "Xio":
        0.194e-9,
        "gam0":
        0.0496,
        "qd":
        0.5e-9,
        "Tcref":
        730.8
    }

    _thermal = thermo0,
Esempio n. 13
0
class R161(MEoS):
    """Multiparameter equation of state for R161"""
    name = "fluoroethane"
    CASNumber = "353-36-6"
    formula = "C2H5F"
    synonym = "R161"
    _refPropName = "R161"
    _coolPropName = "R161"
    rhoc = unidades.Density(302.0010921)
    Tc = unidades.Temperature(375.25)
    Pc = unidades.Pressure(5046.0, "kPa")
    M = 48.0595  # g/mol
    Tt = unidades.Temperature(130.0)
    Tb = unidades.Temperature(235.6)
    f_acent = 0.216
    momentoDipolar = unidades.DipoleMoment(1.9397, "Debye")
    id = 247

    Fi1 = {
        "ao_log": [1, 3],
        "pow": [0, 1],
        "ao_pow": [-6.9170707460, 5.4837900434],
        "ao_exp": [1.08888, 1.80842, 8.72417, 5.67715],
        "titao": [329 / Tc, 742 / Tc, 1644 / Tc, 3922 / Tc]
    }

    Fi2 = {
        "ao_log": [1, 3],
        "pow": [0, 1],
        "ao_pow": [-6.9187, 5.4788],
        "ao_exp": [2.059, 9.253, 6.088],
        "titao": [420 / Tc, 1548 / Tc, 3882 / Tc]
    }

    CP1 = {
        "ao": 3.985,
        "ao_exp": [2.077, 9.265, 6.054],
        "exp": [420, 1548, 3882]
    }

    qi = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for R-161 of Qi (2012)",
        "__doi__": {
            "autor": "Qi, H., Fang, D., Gao, K., Meng, X., Wu, J.",
            "title": "Compressed Liquid Densities and Helmholtz Energy"
            " Equation of State for Fluoroethane (R161)",
            "ref": "Int. J. Thermophys. 37(3) (2016) 55",
            "doi": "10.1007/s10765-016-2061-1"
        },
        "R":
        8.3144621,
        "Tc":
        375.25,
        "rhoc":
        6.2839,
        "cp":
        Fi1,
        "ref":
        "IIR",
        "Tmin":
        Tt,
        "Tmax":
        420.0,
        "Pmax":
        100000.0,
        "rhomax":
        18,
        "nr1": [
            0.005145283, -0.001882274, 1.884722, -3.1819965, -0.24432415,
            0.27792467
        ],
        "d1": [5, 4, 1, 1, 2, 3],
        "t1": [1, 0.68, 0.32, 0.92, 1.23, 0.846],
        "nr2": [
            -0.4414064, -0.402065, 0.24171113, -0.16603585, -0.03440867,
            -0.000099185
        ],
        "d2": [1, 3, 2, 2, 7, 5],
        "t2": [4.208, 3.06, 1.85, 4.28, 1.003, 1.12],
        "c2": [2, 2, 1, 2, 1, 1],
        "gamma2": [1] * 6,
        "nr3": [
            1.0146668, -0.03542609, -0.006038245, -0.025437558, -0.00515678,
            0.006396804
        ],
        "d3": [1, 1, 3, 3, 2, 2],
        "t3": [1.055, 0.8, 4.08, 1.6, 3.85, 0.57],
        "alfa3": [0.96212, 3.2147, 2.6288, 0.8657, 2.3839, 1.7814],
        "beta3": [0.62848, 4.5968, 4.9696, 0.239, 0.788, 7.0874],
        "gamma3": [1.9363, 1.5054, 1.3691, 2.3594, 0.5581, 0.6326],
        "epsilon3": [0.70192, 1.23824, 0.73324, 0.6258, 1.564, 1.4861]
    }

    wu = {
        "__type__": "Helmholtz",
        "__name__": "short Helmholtz equation of state for R-161 of Wu and "
        "Zhou (2012)",
        "__doi__": {
            "autor": "Wu, J., Zhou, Y.",
            "title": "An Equation of State for Fluoroethane (R161)",
            "ref": "Int. J. Thermophys. 33(2) (2012) 220-234",
            "doi": "10.1007/s10765-011-1151-3"
        },
        "R": 8.314472,
        "cp": Fi2,
        "ref": {
            "Tref": 273.15,
            "Pref": 1.,
            "ho": 28559.6,
            "so": 167.205
        },
        "Tmin": Tt,
        "Tmax": 450.0,
        "Pmax": 5000.0,
        "rhomax": 20.0,
        "nr1": [1.511, -2.3, -0.457, 0.1683, 0.04133],
        "d1": [1, 1, 2, 3, 4],
        "t1": [0.37, 0.97, 1.14, 0.744, 1.],
        "nr2": [0.62187, -0.0265, -1.03, -0.285, -0.476],
        "d2": [2, 7, 1, 2, 3],
        "t2": [1.26, 1., 1.8, 3., 2.25],
        "c2": [1, 1, 2, 2, 2],
        "gamma2": [1] * 5,
        "nr3": [0.82, -0.3532, -0.116, -0.0220583, -1.63148],
        "d3": [11, 1, 3, 3, 3],
        "t3": [1, 1.2, 5.3, 1, 4],
        "alfa3": [0.96, 1.35, 1.26, 1.23, 16.8],
        "beta3": [2.7, 5.2, 3.9, 4.7, 413],
        "gamma3": [0.9, 0.69, 0.67, 0.67, 1.15],
        "epsilon3": [0.683, 0.892, 0.785, 1.33, 0.86]
    }

    refprop = {
        "__type__": "Helmholtz",
        "__name__": "short Helmholtz equation of state for R-161 of Lemmon "
        "(2005).",
        "__doi__": {
            "autor": "Lemmon, E.W.",
            "title": "preliminary equation, 2005.",
            "ref": "",
            "doi": ""
        },
        "R": 8.314472,
        "cp": CP1,
        "ref": "NBP",
        "Tmin": Tt,
        "Tmax": 400.0,
        "Pmax": 50000.0,
        "rhomax": 20.0,
        "nr1": [0.75688, -1.4110, -0.63922, 0.055685, 0.00028395],
        "d1": [1, 1, 1, 3, 7],
        "t1": [0.25, 1.25, 1.5, 0.25, 0.875],
        "nr2": [.73357, .67596, .011369, -.56406, -.094362, -.1678, .00034215],
        "d2": [1, 2, 5, 1, 1, 4, 2],
        "t2": [2.375, 2.0, 2.125, 3.5, 6.5, 4.75, 12.5],
        "c2": [1, 1, 1, 2, 2, 2, 2],
        "gamma2": [1] * 7
    }

    eq = qi, wu, refprop

    _surface = {"sigma": [0.05385], "exp": [1.111]}

    _vapor_Pressure = {
        "eq": 3,
        "n": [-8.977955, 25.64713, -80.33162, 89.22478, -34.33593],
        "t": [1.0, 1.5, 1.84, 2.11, 2.47]
    }
    _liquid_Density = {
        "eq": 1,
        "n": [2.11404, 6.25555, -16.12805, 116.00294, -233.6455, 149.54793],
        "t": [0.46, 1.08, 1.71, 3.56, 4.44, 5.51]
    }
    _vapor_Density = {
        "eq": 2,
        "n": [-2.26823, 1.39859, -27.93344, 49.91536, -50.51966],
        "t": [0.3, 0.88, 1.46, 2.1, 2.82]
    }

    visco0 = {
        "__name__":
        "Tsolakidou (2017)",
        "__doi__": {
            "autor": "Tsolakidou, C.M., Assael, M.J., Huber, M.L.,"
            "Perkins, R.A.",
            "title": "Correlations for the Viscosity and Thermal "
            "Conductivity of Ethyl Fluoride (R161)",
            "ref": "J. Phys. Chem. Ref. Data 46(2) (2017) 023103",
            "doi": "10.1063/1.4983027"
        },
        "eq":
        1,
        "omega":
        1,
        "n_chapman":
        0.021357,
        "ek":
        320.39,
        "sigma":
        0.4457,
        "collision": [0.24130, -0.45],
        "Tref_virial":
        320.39,
        "n_virial": [
            -19.572881, 219.73999, -1015.3226, 2471.0125, -3375.1717,
            2491.6597, -787.26086, 14.085455, -0.34664158
        ],
        "t_virial": [0, -0.25, -0.5, -0.75, -1, -1.25, -1.5, -2.5, -5.5],
        "Tref_res":
        Tc,
        "rhoref_res":
        302.001,
        "nr": [-10.28373, 7.65563, 4.842, 0.42223],
        "tr": [-0.5, -0.5, -2.5, -0.5],
        "dr": [2 / 3, 5 / 3, 2 / 3, 14 / 3],
        "nr_num": [64.34983, 64.34983],
        "tr_num": [-1.5, -0.5],
        "dr_num": [2 / 3, 5 / 3],
        # The Eq 8 in paper have a typo, the denominator of last term
        # has a plus but the correlation work with a minus
        "nr_den": [10.99213, -1],
        "tr_den": [-2, -2],
        "dr_den": [0, 2]
    }

    _viscosity = visco0,

    thermo0 = {
        "__name__":
        "Tsolakidou (2017)",
        "__doi__": {
            "autor": "Tsolakidou, C.M., Assael, M.J., Huber, M.L.,"
            "Perkins, R.A.",
            "title": "Correlations for the Viscosity and Thermal "
            "Conductivity of Ethyl Fluoride (R161)",
            "ref": "J. Phys. Chem. Ref. Data 46(2) (2017) 023103",
            "doi": "10.1063/1.4983027"
        },
        "eq":
        1,
        "Toref":
        Tc,
        "koref":
        1e-3,
        "no_num": [7.96804, -12.5874, -26.3743, 16.9894, 127.545, -32.548],
        "to_num": [0, 1, 2, 3, 4, 5],
        "no_den": [5.406, -18.8331, 24.868, -9.14139, 1],
        "to_den": [0, 1, 2, 3, 4],
        "Tref_res":
        Tc,
        "rhoref_res":
        302.001,
        "kref_res":
        1e-3,
        "nr": [
            -8.41553, 7.41456, -39.7744, 44.0586, 106.179, -81.9833, -53.2351,
            37.6052, 8.23094, -4.90293
        ],
        "tr": [0, -1, 0, -1, 0, -1, 0, -1, 0, -1],
        "dr": [1, 1, 2, 2, 3, 3, 4, 4, 5, 5],
        "critical":
        3,
        "gnu":
        0.63,
        "gamma":
        1.239,
        "R0":
        1.02,
        "Xio":
        0.183e-9,
        "gam0":
        0.055,
        "qd":
        3.104e-10,
        "Tcref":
        562.88
    }

    _thermal = thermo0,
Esempio n. 14
0
 def editorPerfil(self):
     dialog=inputTable.InputTableDialog(2, data=self.profile_T, title=QtWidgets.QApplication.translate("pychemqt", "Temperature profile"), horizontalHeader=["x", "T, "+unidades.Temperature(None).text()])
     if dialog.exec_():
         self.profile_T=dialog.data
Esempio n. 15
0
class SF6(MEoS):
    """Multiparameter equation of state for sulfur hexafluoride"""
    name = "sulfur hexafluoride"
    CASNumber = "2551-62-4"
    formula = "SF6"
    synonym = ""
    rhoc = unidades.Density(742.3)
    Tc = unidades.Temperature(318.7232)
    Pc = unidades.Pressure(3754.983, "kPa")
    M = 146.0554192  # g/mol
    Tt = unidades.Temperature(223.555)
    Tb = unidades.Temperature(204.9)
    f_acent = 0.21
    momentoDipolar = unidades.DipoleMoment(0.0, "Debye")
    # id = 953
    id = 1
    _Tr = unidades.Temperature(304.013497)
    _rhor = unidades.Density(747.815849)
    _w = 0.181815238

    Fi1 = {
        "ao_log": [1, 3],
        "pow": [0, 1],
        "ao_pow": [11.638611086, -6.392241811],
        "ao_exp": [3.66118232, 7.87885103, 3.45981679],
        "titao": [1.617282065, 2.747115139, 4.232907175],
        "ao_hyp": [],
        "hyp": []
    }

    CP1 = {
        "ao":
        3.9837756784,
        "an": [],
        "pow": [],
        "ao_exp": [
            2.2181851010, -1.0921337374e1, 3.3102497939, 17.5189671483,
            2.8903523803
        ],
        "exp": [1114.38, 925.64, 499.26, 884.9, 1363.93],
        "ao_hyp": [],
        "hyp": []
    }

    CP2 = {
        "ao": -0.376915e-1 / 8.3143 * 146.05,
        "an": [0.305814e-2 / 8.3143 * 146.05, -0.237654e-5 / 8.3143 * 146.05],
        "pow": [1, 2],
        "ao_exp": [],
        "exp": [],
        "ao_hyp": [],
        "hyp": []
    }

    helmholtz1 = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for sulfur hexafluoride of Guder and Wagner (2009)",
        "__doi__": {
            "autor": "Guder, C. and Wagner, W.",
            "title":
            "A Reference Equation of State for the Thermodynamic Properties of Sulfur Hexafluoride for Temperatures from the Melting Line to 625 K and Pressures up to 150 MPa",
            "ref": "J. Phys. Chem. Ref. Data 38, 33 (2009)",
            "doi": "10.1063/1.3037344"
        },
        "__test__":
        # Table 27, Pag 57
        """
            >>> wt=SF6()
            >>> tau=wt.Tc/350
            >>> delta=436.9770888/wt.rhoc
            >>> print "%0.9g %0.9g %0.9g %0.9g %0.9g %0.9g" % wt._phi0(wt._constants["cp"], tau, delta)
            3.30559888 0.91277072 -14.4662979 1.69871606 -2.88563626 0
            >>> print "%0.9g %0.9g %0.9g %0.9g %0.9g %0.9g" % wt._phir(tau, delta)[:6]
            -0.496581463 -1.37926327 -1.37917096 -0.723171558 0.405086373 -2.09574715
            """

        # Table 28, Pag 71
        """
            >>> st=SF6(T=223.555, x=0.5)
            >>> print "%0.6g %0.6g %0.2f %0.2f %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, st.Gas.h.kJkg, \
                st.Liquido.s.kJkgK, st.Gas.s.kJkgK, st.Liquido.cv.kJkgK, st.Gas.cv.kJkgK, \
                st.Liquido.cp.kJkgK, st.Gas.cp.kJkgK, st.Liquido.w, st.Gas.w)
            223.555 0.231425 1845.03 19.56 -158.14 -47.521 -0.72226 -0.22745 0.52753 0.48331 0.83712 0.56309 552.26 112.84
            >>> st=SF6(T=230, x=0.5)
            >>> print "%0.6g %0.5g %0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, st.Gas.h.kJkg, \
                st.Liquido.s.kJkgK, st.Gas.s.kJkgK, st.Liquido.cv.kJkgK, st.Gas.cv.kJkgK, \
                st.Liquido.cp.kJkgK, st.Gas.cp.kJkgK, st.Liquido.w, st.Gas.w)
            230 0.30133 1812.92 25.143 -152.66 -44.796 -0.69829 -0.2293 0.54069 0.49958 0.8573 0.58507 521.49 112.81
            >>> st=SF6(T=240, x=0.5)
            >>> print "%0.6g %0.6g %0.2f %0.2f %0.2f %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, st.Gas.h.kJkg, \
                st.Liquido.s.kJkgK, st.Gas.s.kJkgK, st.Liquido.cv.kJkgK, st.Gas.cv.kJkgK, \
                st.Liquido.cp.kJkgK, st.Gas.cp.kJkgK, st.Liquido.w, st.Gas.w)
            240 0.440064 1760.80 36.21 -143.91 -40.63 -0.66135 -0.23103 0.56043 0.52495 0.89018 0.62193 474.67 112.27
            >>> st=SF6(T=250, x=0.5)
            >>> print "%0.6g %0.6g %0.2f %0.3f %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, st.Gas.h.kJkg, \
                st.Liquido.s.kJkgK, st.Gas.s.kJkgK, st.Liquido.cv.kJkgK, st.Gas.cv.kJkgK, \
                st.Liquido.cp.kJkgK, st.Gas.cp.kJkgK, st.Liquido.w, st.Gas.w)
            250 0.621964 1705.31 50.867 -134.81 -36.589 -0.62463 -0.23177 0.57964 0.55056 0.92638 0.66356 428.57 111.04
            >>> st=SF6(T=260, x=0.5)
            >>> print "%0.6g %0.6g %0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, st.Gas.h.kJkg, \
                st.Liquido.s.kJkgK, st.Gas.s.kJkgK, st.Liquido.cv.kJkgK, st.Gas.cv.kJkgK, \
                st.Liquido.cp.kJkgK, st.Gas.cp.kJkgK, st.Liquido.w, st.Gas.w)
            260 0.854649 1645.54 70.103 -125.33 -32.747 -0.58803 -0.23193 0.5986 0.57618 0.96814 0.71212 382.8 109.02
            >>> st=SF6(T=270, x=0.5)
            >>> print "%0.6g %0.7g %0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, st.Gas.h.kJkg, \
                st.Liquido.s.kJkgK, st.Gas.s.kJkgK, st.Liquido.cv.kJkgK, st.Gas.cv.kJkgK, \
                st.Liquido.cp.kJkgK, st.Gas.cp.kJkgK, st.Liquido.w, st.Gas.w)
            270 1.146251 1580.17 95.322 -115.43 -29.213 -0.55134 -0.23202 0.61759 0.60262 1.0194 0.77387 336.89 106.11
            >>> st=SF6(T=280, x=0.5)
            >>> print "%0.6g %0.7g %0.2f %0.2f %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, st.Gas.h.kJkg, \
                st.Liquido.s.kJkgK, st.Gas.s.kJkgK, st.Liquido.cv.kJkgK, st.Gas.cv.kJkgK, \
                st.Liquido.cp.kJkgK, st.Gas.cp.kJkgK, st.Liquido.w, st.Gas.w)
            280 1.505551 1507.05 128.67 -105.01 -26.139 -0.51431 -0.23262 0.63716 0.63271 1.0875 0.86231 290.16 102.15
            >>> st=SF6(T=290, x=0.5)
            >>> print "%0.6g %0.7g %0.2f %0.2f %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, st.Gas.h.kJkg, \
                st.Liquido.s.kJkgK, st.Gas.s.kJkgK, st.Liquido.cv.kJkgK, st.Gas.cv.kJkgK, \
                st.Liquido.cp.kJkgK, st.Gas.cp.kJkgK, st.Liquido.w, st.Gas.w)
            290 1.942334 1422.51 173.85 -93.928 -23.766 -0.47647 -0.23453 0.6589 0.66504 1.1897 0.99972 241.62 96.919
            >>> st=SF6(T=300, x=0.5)
            >>> print "%0.6g %0.7g %0.2f %0.2f %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, st.Gas.h.kJkg, \
                st.Liquido.s.kJkgK, st.Gas.s.kJkgK, st.Liquido.cv.kJkgK, st.Gas.cv.kJkgK, \
                st.Liquido.cp.kJkgK, st.Gas.cp.kJkgK, st.Liquido.w, st.Gas.w)
            300 2.468172 1319.00 238.43 -81.864 -22.596 -0.43688 -0.23932 0.68633 0.7029 1.3801 1.2722 189.8 90.145
            >>> st=SF6(T=310, x=0.5)
            >>> print "%0.6g %0.7g %0.2f %0.2f %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, st.Gas.h.kJkg, \
                st.Liquido.s.kJkgK, st.Gas.s.kJkgK, st.Liquido.cv.kJkgK, st.Gas.cv.kJkgK, \
                st.Liquido.cp.kJkgK, st.Gas.cp.kJkgK, st.Liquido.w, st.Gas.w)
            310 3.098401 1174.67 344.14 -67.924 -24.009 -0.39285 -0.25119 0.73084 0.76178 1.9613 2.1947 131.47 81.329
            >>> st=SF6(T=318, x=0.5)
            >>> print "%0.6g %0.7g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, st.Gas.h.kJkg, \
                st.Liquido.s.kJkgK, st.Gas.s.kJkgK, st.Liquido.cv.kJkgK, st.Gas.cv.kJkgK, \
                st.Liquido.cp.kJkgK, st.Gas.cp.kJkgK, st.Liquido.w, st.Gas.w)
            318 3.695565 917.87 569.56 -50.87 -33.648 -0.34041 -0.28626 0.8803 0.93693 16.37 26.348 70.232 69.747
            >>> st=SF6(T=318.7, x=0.5)
            >>> print "%0.6g %0.7g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.6g %0.6g %0.5g %0.5g" % (\
                st.T, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, st.Gas.h.kJkg, \
                st.Liquido.s.kJkgK, st.Gas.s.kJkgK, st.Liquido.cv.kJkgK, st.Gas.cv.kJkgK, \
                st.Liquido.cp.kJkgK, st.Gas.cp.kJkgK, st.Liquido.w, st.Gas.w)
            318.7 3.753053 787.75 696.67 -44.789 -40.391 -0.32152 -0.30773 1.0458 1.0721 1027.68 1273.3 60.578 62.816
            """

        # Table 29, Pag 73
        """
            >>> st=SF6(T=390, P=1e5)
            >>> print "%0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.rho, st.u.kJkg, st.h.kJkg, st.s.kJkgK, st.cv.kJkgK, st.cp.kJkgK, st.w)
            390 4.5244 44.587 66.69 0.19538 0.73041 0.78868 154.14
            >>> st=SF6(T=300, P=5e5)
            >>> print "%0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.rho, st.u.kJkg, st.h.kJkg, st.s.kJkgK, st.cv.kJkgK, st.cp.kJkgK, st.w)
            300 31.041 -18.062 -1.9544 -0.094212 0.61818 0.69206 130.24
            >>> st=SF6(T=225, P=1e6)
            >>> print "%0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.rho, st.u.kJkg, st.h.kJkg, st.s.kJkgK, st.cv.kJkgK, st.cp.kJkgK, st.w)
            225 1841.9 -157.3 -156.76 -0.71798 0.53068 0.83889 550.68
            >>> st=SF6(T=600, P=1.5e6)
            >>> print "%0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.rho, st.u.kJkg, st.h.kJkg, st.s.kJkgK, st.cv.kJkgK, st.cp.kJkgK, st.w)
            600 44.273 213.8 247.68 0.41151 0.87699 0.94036 189.91
            >>> st=SF6(T=400, P=2e6)
            >>> print "%0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.rho, st.u.kJkg, st.h.kJkg, st.s.kJkgK, st.cv.kJkgK, st.cp.kJkgK, st.w)
            400 95.798 46.897 67.774 0.032259 0.74992 0.83896 146.21
            >>> st=SF6(T=305, P=3e6)
            >>> print "%0.6g %0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.rho, st.u.kJkg, st.h.kJkg, st.s.kJkgK, st.cv.kJkgK, st.cp.kJkgK, st.w)
            305 1272.55 -78.218 -75.861 -0.41839 0.6977 1.4717 173.33
            >>> st=SF6(T=300, P=4e6)
            >>> print "%0.6g %0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.rho, st.u.kJkg, st.h.kJkg, st.s.kJkgK, st.cv.kJkgK, st.cp.kJkgK, st.w)
            300 1382.35 -86.689 -83.795 -0.44709 0.6729 1.1916 234.55
            >>> st=SF6(T=625, P=5e6)
            >>> print "%0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.rho, st.u.kJkg, st.h.kJkg, st.s.kJkgK, st.cv.kJkgK, st.cp.kJkgK, st.w)
            625 142.37 231.78 266.9 0.37513 0.89076 0.9678 194.91
            >>> st=SF6(T=260, P=6e6)
            >>> print "%0.6g %0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.rho, st.u.kJkg, st.h.kJkg, st.s.kJkgK, st.cv.kJkgK, st.cp.kJkgK, st.w)
            260 1693.94 -128.63 -125.08 -0.59892 0.59673 0.92285 436.07
            >>> st=SF6(T=300, P=1e7)
            >>> print "%0.6g %0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.rho, st.u.kJkg, st.h.kJkg, st.s.kJkgK, st.cv.kJkgK, st.cp.kJkgK, st.w)
            300 1505.13 -92.926 -86.282 -0.46916 0.66145 1.0186 328.66
            >>> st=SF6(T=240, P=2e7)
            >>> print "%0.6g %0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.rho, st.u.kJkg, st.h.kJkg, st.s.kJkgK, st.cv.kJkgK, st.cp.kJkgK, st.w)
            240 1862.82 -150.39 -139.65 -0.68852 0.56357 0.83019 603.24
            >>> st=SF6(T=400, P=5e7)
            >>> print "%0.6g %0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.rho, st.u.kJkg, st.h.kJkg, st.s.kJkgK, st.cv.kJkgK, st.cp.kJkgK, st.w)
            400 1479.33 -19.467 14.333 -0.25764 0.77463 0.98607 437.03
            >>> st=SF6(T=500, P=1e8)
            >>> print "%0.6g %0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.rho, st.u.kJkg, st.h.kJkg, st.s.kJkgK, st.cv.kJkgK, st.cp.kJkgK, st.w)
            500 1515.39 60.333 126.32 -0.082538 0.85386 1.0148 553.13
            >>> st=SF6(T=625, P=1.5e8)
            >>> print "%0.6g %0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\
                st.T, st.rho, st.u.kJkg, st.h.kJkg, st.s.kJkgK, st.cv.kJkgK, st.cp.kJkgK, st.w)
            625 1505.25 171.72 271.37 0.11704 0.91497 1.0498 633.74
            """,
        "R":
        8.314472,
        "cp":
        Fi1,
        "ref":
        "OTO",
        "Tmin":
        Tt,
        "Tmax":
        625.0,
        "Pmax":
        150000.0,
        "rhomax":
        14.5,
        "Pmin":
        231.429,
        "rhomin":
        12.632,
        "nr1": [
            .54958259132835, -.87905033269396, -.84656969731452,
            .27692381593529, -.49864958372345e01, .48879127058055e01,
            .36917081634281e-1, .37030130305087e-3, .39389132911585e-1,
            .42477413690006e-3
        ],
        "d1": [1, 1, 1, 2, 2, 2, 3, 3, 4, 6],
        "t1": [0.125, 1.25, 1.875, 0.125, 1.5, 1.625, 1.5, 5.625, 0.625, 0.25],
        "nr2": [
            -.24150013863890e-1, .59447650642255e-1, -.38302880142267,
            .32606800951983, -.29955940562031e-1, -.86579186671173e-1,
            .41600684707562e01, -.41398128855814e01, -.55842159922714,
            .56531382776891, .82612463415545e-2, -.10200995338080e-1
        ],
        "d2": [1, 2, 2, 2, 3, 6, 2, 2, 4, 4, 2, 2],
        "t2":
        [6., 0.25, 4.75, 5.375, 5.875, 2., 5.875, 6., 5.625, 5.75, 0., 0.5],
        "c2": [1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 3, 3],
        "gamma2": [1] * 12,
        "nr3": [
            -.21662523861406e-1, .34650943893908e-1, -.28694281385812e-1,
            .84007238998053e-2, -.26969359922498, .90415215646344e01,
            -.37233103557977e01, -.27524670823704e04, .57711861697319e04,
            -.30234003119748e04, .22252778435360e07, -.23056065559032e07,
            .63918852944475e07, -.60792091415592e07
        ],
        "d3": [1, 3, 4, 1, 1, 4, 3, 4, 4, 4, 1, 1, 3, 3],
        "t3": [4, 1, 3, 2, 4, 3, 4, 1, 2, 3, 3, 4, 3, 4],
        "alfa3": [10, 10, 10, 10, 11, 25, 30, 30, 30, 30, 30, 30, 30, 30],
        "beta3":
        [150, 150, 150, 150, 225, 300, 350, 350, 350, 350, 400, 400, 400, 400],
        "gamma3": [
            1.13, 1.13, 1.13, 1.16, 1.19, 1.19, 1.16, 1.16, 1.16, 1.16, 1.22,
            1.22, 1.22, 1.22
        ],
        "epsilon3": [0.85, 0.85, 0.85, 0.85, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
    }

    helmholtz2 = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for sulfur hexafluoride of de Reuck et al. (1991)",
        "__doi__": {
            "autor": "de Reuck, K.M., Craven, R.J.B., and Cole, W.A.",
            "title":
            "Report on the Development of an Equation of State for Sulphur Hexafluoride",
            "ref": "IUPAC Thermodynamic Tables Project Centre, London, 1991.",
            "doi": ""
        },
        "R":
        8.31448,
        "cp":
        CP1,
        "ref":
        "OTO",
        "Tmin":
        Tt,
        "Tmax":
        525.0,
        "Pmax":
        55000.0,
        "rhomax":
        12.7,
        "Pmin":
        224.36,
        "rhomin":
        12.677,
        "nr1": [
            0.26945570453, -0.554046585076, -0.929624636454, 0.505661081063,
            -0.683495847809, 0.579161832426, -0.122636218956,
            -0.260339227668e-1, 0.222201648687e-1, -0.118992341472e-2,
            0.292000609763e-2, -0.243315775571e-2, 0.689778297550e-3
        ],
        "d1": [1, 1, 1, 2, 2, 2, 3, 4, 5, 10, 10, 10, 10],
        "t1": [0, 1.5, 2, 0, 1, 2, 0, 2, 0, 0.5, 1, 1.5, 2],
        "nr2": [
            -0.147585329235e1, 0.275952303526e1, -0.142721418498e1,
            0.598794196648e-1, 0.219991168025e-2, 0.746554473361e-2,
            0.345233637389e-2, -0.253226231963e-1, 0.433906886402e-1,
            -0.249349699078e-1, 0.338560952242e-2, 0.539985899700e-3
        ],
        "d2": [2, 2, 2, 3, 7, 7, 9, 4, 4, 4, 6, 4],
        "t2": [3, 4, 5, 5, 1, 5, 1, 9, 14, 24, 24, 9],
        "c2": [2, 2, 2, 2, 2, 2, 2, 4, 4, 4, 4, 6],
        "gamma2": [1] * 12
    }

    helmholtz3 = {
        "__type__":
        "Helmholtz",
        "__name__":
        "short Helmholtz equation of state for sulfur hexafluoride of Span and Wagner (2003)",
        "__doi__": {
            "autor": "Span, R., Wagner, W.",
            "title":
            "Equations of state for technical applications. II. Results for nonpolar fluids.",
            "ref": "Int. J. Thermophys. 24 (2003), 41 – 109.",
            "doi": "10.1023/A:1022310214958"
        },
        "__test__":
        """
            >>> st=SF6(T=700, rho=200, eq=2)
            >>> print "%0.4f %0.3f %0.4f" % (st.cp0.kJkgK, st.P.MPa, st.cp.kJkgK)
            0.9671 8.094 0.9958
            >>> st2=SF6(T=750, rho=100, eq=2)
            >>> print "%0.2f %0.5f" % (st2.h.kJkg-st.h.kJkg, st2.s.kJkgK-st.s.kJkgK)
            52.80 0.10913
            """,  # Table III, Pag 46
        "R":
        8.31451,
        "cp":
        CP1,
        "ref":
        "OTO",
        "Tmin":
        Tt,
        "Tmax":
        750.0,
        "Pmax":
        100000.0,
        "rhomax":
        12.65,
        "Pmin":
        221.22,
        "rhomin":
        12.645,
        "nr1": [
            0.12279403e1, -0.33035623e1, 0.12094019e1, -0.12316, 0.11044657,
            0.32952153e-3
        ],
        "d1": [1, 1, 1, 2, 3, 7],
        "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875],
        "nr2": [
            0.27017629, -0.62910351e-1, -0.3182889, -0.99557419e-1,
            -0.36909694e-1, 0.19136427e-1
        ],
        "d2": [2, 5, 1, 4, 3, 4],
        "t2": [0.625, 1.75, 3.625, 3.625, 14.5, 12],
        "c2": [1, 1, 2, 2, 3, 3],
        "gamma2": [1] * 7
    }

    helmholtz4 = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for sulfur hexafluoride of Polt et al. (1992)",
        "__doi__": {
            "autor": "Polt, A., Platzer, B., and Maurer, G.",
            "title":
            "Parameter der thermischen Zustandsgleichung von Bender fuer 14 mehratomige reine Stoffe",
            "ref": "Chem. Technik 22(1992)6 , 216/224",
            "doi": ""
        },
        "R":
        8.3143,
        "cp":
        CP2,
        "ref":
        "OTO",
        "Tmin":
        Tt,
        "Tmax":
        523.0,
        "Pmax":
        40000.0,
        "rhomax":
        13.133,
        "Pmin":
        236.73,
        "rhomin":
        12.712,
        "nr1": [
            0.131111896375, -0.792338803106, 0.580899809209, 0.153233600406e1,
            -0.485096079094e1, 0.482411603806e1, -0.311285647219e1,
            0.442141211276, 0.206313183222, -0.372305169645, 0.443536383059,
            -0.476354850910e-1, 0.116313319336, 0.570240883234e-1,
            -0.152963195118, 0.259842094503e-1
        ],
        "d1": [0, 0, 0, 1, 1, 1, 1, 1, 2, 2, 2, 3, 3, 4, 4, 5],
        "t1": [3, 4, 5, 0, 1, 2, 3, 4, 0, 1, 2, 0, 1, 0, 1, 1],
        "nr2": [
            -0.131111896375, 0.792338803106, -0.580899809209, -0.744763581796,
            0.204368923925e1, -0.129335324120e1
        ],
        "d2": [0, 0, 0, 2, 2, 2],
        "t2": [3, 4, 5, 3, 4, 5],
        "c2": [2] * 6,
        "gamma2": [1.32678063] * 6
    }

    eq = helmholtz1, helmholtz2, helmholtz3, helmholtz4

    _surface = {"sigma": [0.0538, -4.064e-5], "exp": [1.271, 0.2116]}
    _melting = {
        "eq": 1,
        "Tref": Tt,
        "Pref": 0.48475e-4,
        "Tmin": Tt,
        "Tmax": 800.0,
        "a1": [1., -30.0468473, 30.0468473, 359.771253, -359.771253],
        "exp1": [0, -20., 0, 3.25, 0],
        "a2": [],
        "exp2": [],
        "a3": [],
        "exp3": []
    }
    _sublimation = {
        "eq": 2,
        "Tref": Tt,
        "Pref": 231.429,
        "Tmin": Tt,
        "Tmax": Tt,
        "a1": [-11.6942141, 11.6942141],
        "exp1": [-1.07, 0],
        "a2": [],
        "exp2": [],
        "a3": [],
        "exp3": []
    }
    _vapor_Pressure = {
        "eq": 6,
        "ao": [-7.09634642, 1.676662, -2.3921599, 5.86078302, -9.02978735],
        "exp": [2.0, 3.0, 5.0, 8.0, 9.0]
    }
    _liquid_Density = {
        "eq": 6,
        "ao": [2.31174688, -1.12912486, -1.439347, 0.282489982],
        "exp": [1.065, 1.5, 4., 5.]
    }
    _vapor_Density = {
        "eq":
        6,
        "ao": [
            23.68063442, 0.513062232, -24.4706238, -4.6715244, -1.7536843,
            -6.65585369
        ],
        "exp": [1.044, 0.5, 1.0, 2.0, 8.0, 17.]
    }
Esempio n. 16
0
class R40(MEoS):
    """Multiparameter equation of state for R40"""
    name = "methyl chloride"
    CASNumber = "74-87-3"
    formula = "CH3Cl"
    synonym = "R40"
    _refPropName = "R40"
    _coolPropName = "R40"
    rhoc = unidades.Density(363.219)
    Tc = unidades.Temperature(416.3)
    Pc = unidades.Pressure(6677.3, "kPa")
    M = 50.48752  # g/mol
    Tt = unidades.Temperature(175.0)
    Tb = unidades.Temperature(249.173)
    f_acent = 0.243
    momentoDipolar = unidades.DipoleMoment(1.871, "Debye")
    id = 115

    Fi1 = {
        "ao_log": [1, 2.92518],
        "pow": [0, 1, -1, -2, -3],
        "ao_pow": [7.499423, -2.997533, -6.0842e-2, -1.1525e-1, 1.0843e-2],
        "ao_exp": [3.764997],
        "titao": [3.7101]
    }

    thol = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for R40 of Thol et al (2013)",
        "__doi__": {
            "autor": "Thol, M., Piazza, L., and Span, R.",
            "title": "A New Functional Form for Equations of State "
            "for Some Weakly Associating Fluids",
            "ref": "Int. J. Thermophys., 35(5):783-811, 2014.",
            "doi": "10.1007/s10765-014-1633-1"
        },
        "R":
        8.314472,
        "cp":
        Fi1,
        "ref":
        "OTO",
        "Tmin":
        230.0,
        "Tmax":
        630.0,
        "Pmax":
        100000.0,
        "rhomax":
        21.78756,
        "Pmin":
        0.9,
        "rhomin":
        21.78756,
        "nr1": [
            .274572058, .104235924e-1, -.125727710e1, 2.25609199e-3,
            -3.31830421e-2, .918440878e-1, .261059608e-2
        ],
        "d1": [1, 1, 1, 2, 3, 3, 5],
        "t1": [-0.75, -0.25, 1.25, 0.75, -1.0, -0.375, 1.25],
        "nr2": [
            -.948880966e-1, -.843634836e-1, .226263660, -.470765940e-1,
            -.196610405, -.204318929e-1, -.692145009e-1, .148974844e-1,
            -.642544485e-2
        ],
        "d2": [1, 1, 2, 5, 1, 3, 4, 5, 2],
        "t2": [2.375, 3.0, 2.625, 1.875, 4.5, 5.75, 5.375, 2.75, 14.5],
        "c2": [1, 1, 1, 1, 2, 2, 2, 2, 3],
        "gamma2": [1] * 9
    }

    eq = thol,

    _vapor_Pressure = {
        "eq": 3,
        "n": [-6.5074, 0.7520, -9.4148, 19.654, -20.190],
        "t": [1.0, 1.5, 4.5, 5.8, 7.1]
    }
    _liquid_Density = {
        "eq": 1,
        "n": [2.1809, 0.9228, -2.4615, 7.9722, -13.023, 9.2840],
        "t": [0.37, 1.16, 2.0, 2.9, 3.9, 5.1]
    }
    _vapor_Density = {
        "eq": 2,
        "n": [-0.9433, -6.8001, -82.752, 202.14, -264.16, 99.135],
        "t": [0.18, 0.9, 3.7, 4.6, 5.6, 6.7]
    }
Esempio n. 17
0
class D4(MEoS):
    """Multiparameter equation of state for octamethylcyclotetrasiloxane"""
    name = "octamethylcyclotetrasiloxane"
    CASNumber = "556-67-2"
    formula = "C8H24O4Si4"
    synonym = "D4"
    _refPropName = "D4"
    _coolPropName = "D4"
    rhoc = unidades.Density(307.0335906736056)
    Tc = unidades.Temperature(586.49127187)
    Pc = unidades.Pressure(1332.0, "kPa")
    M = 296.61576  # g/mol
    Tt = unidades.Temperature(290.25)
    Tb = unidades.Temperature(448.504)
    f_acent = 0.592
    momentoDipolar = unidades.DipoleMoment(1.090, "Debye")
    # id=1430

    Fi1 = {
        "ao_log": [1, 3],
        "pow": [0, 1],

        # The paper parameter are incorrect, using alternate parameter
        # to meet reference state OTO
        # "ao_pow": [71.1636049792958, -21.6743650975623],
        "ao_pow": [44.72889170669655, -4.9687471148991968],
        "ao_exp": [0.292757, 38.2456, 58.975],
        "titao": [40 / Tc, 200 / Tc, 1800 / Tc]
    }

    f = 8.314472
    CP1 = {
        "ao": -18.256 / f,
        "an": [1427.2e-3 / f, -990.20e-6 / f, 300.0e-9 / f],
        "pow": [1, 2, 3],
        "ao_exp": [],
        "exp": []
    }

    thol = {
        "__type__": "Helmholtz",
        "__name__": "Helmholtz equation of state for D4 of Thol (2006).",
        "__doi__": {
            "autor":
            "Thol, M.",
            "title":
            "Empirical Multiparameter Equations of State "
            "Based on Molecular Simulation and Hybrid Data "
            "Sets",
            "ref":
            "PhD thesis, Ruhr-Universität Bochum, 2015.",
            "doi":
            ""
        },
        "R": 8.3144621,
        "Tc": 586.5,
        "rhoc": 1.043,
        "Pc": 1347,
        "cp": Fi1,
        "ref": "OTO",
        "Tmin": Tt,
        "Tmax": 1200.0,
        "Pmax": 520000.0,
        "rhomax": 5.266,
        "nr1": [5.273743e-2, 4.176401, -4.737070, -1.289588, 5.272749e-1],
        "d1": [4, 1, 1, 2, 3],
        "t1": [1, 0.27, 0.51, 0.998, 0.56],
        "nr2": [-2.558391, -0.9726737, 0.7208209, -4.789456e-1, -5.563239e-2],
        "d2": [1, 3, 2, 2, 7],
        "t2": [1.75, 3.09, 0.79, 2.71, 0.998],
        "c2": [2, 2, 1, 2, 1],
        "gamma2": [1] * 5,
        "nr3": [3.766589, 8.786997e-2, -1.267646e-1, -1.004246, -1.641887],
        "d3": [1, 1, 3, 2, 2],
        "t3": [0.93, 3.17, 1.08, 1.41, 0.89],
        "alfa3": [0.861, 1.114, 1.01, 1.11, 1.032],
        "beta3": [0.75, 0.55, 1.0, 0.47, 1.36],
        "gamma3": [1.124, 1.388, 1.148, 1.197, 0.817],
        "epsilon3": [0.926, 1.3, 1.114, 0.996, 0.483],
        "nr4": []
    }

    colonna = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for D4 of Colonna (2006).",
        "__doi__": {
            "autor": "Colonna, P., Nannan, N.R., Guardone, A., "
            "Lemmon, E.W.",
            "title": "Multiparameter Equations of State for Selected "
            "Siloxanes",
            "ref": "Fluid Phase Equilibria, 244:193-211, 2006.",
            "doi": "10.1016/j.fluid.2006.04.015"
        },
        "R":
        8.314472,
        "cp":
        CP1,
        "ref":
        "NBP",
        "Tmin":
        300.0,
        "Tmax":
        673.0,
        "Pmax":
        30000.0,
        "rhomax":
        3.21,
        "nr1": [
            1.05392408, -2.22981918, 0.77573923, -0.6937405, 0.18721557,
            0.42193330e-3
        ],
        "d1": [1, 1, 1, 2, 3, 7],
        "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875],
        "nr2": [
            0.70301835, 0.47851888e-1, -0.8025348, -0.18968872, -0.22211781e-1,
            0.60103354e-2
        ],
        "d2": [2, 5, 1, 4, 3, 4],
        "t2": [0.625, 1.75, 3.625, 3.625, 14.5, 12.0],
        "c2": [1, 1, 2, 2, 3, 3],
        "gamma2": [1] * 6
    }

    eq = thol, colonna

    _vapor_Pressure = {
        "eq": 3,
        "n": [-0.87935e1, 0.27204e1, -0.48174e1, -0.69086e1],
        "t": [1.0, 1.5, 2.2, 4.4]
    }
    _liquid_Density = {
        "eq": 1,
        "n": [0.14563e1, -0.94215, 0.45065e1, -0.27688e1, 0.8745],
        "t": [0.24, 0.5, 0.75, 1.0, 2.0]
    }
    _vapor_Density = {
        "eq":
        2,
        "n": [
            -0.16204e1, -0.57888e1, -0.24291e2, 0.53567e2, -0.12135e3,
            -0.10976e4
        ],
        "t": [0.31, 0.78, 2.5, 4.4, 5.0, 15.0]
    }
Esempio n. 18
0
class nC5(MEoS):
    """Multiparameter equation of state for n-pentane"""
    name = "pentane"
    CASNumber = "109-66-0"
    formula = "CH3-(CH2)3-CH3"
    synonym = "R-601"
    _refPropName = "PENTANE"
    _coolPropName = "n-Pentane"
    rhoc = unidades.Density(232.)
    Tc = unidades.Temperature(469.7)
    Pc = unidades.Pressure(3370.0, "kPa")
    M = 72.14878  # g/mol
    Tt = unidades.Temperature(143.47)
    Tb = unidades.Temperature(309.21)
    f_acent = 0.251
    momentoDipolar = unidades.DipoleMoment(0.07, "Debye")
    id = 8
    _Tr = unidades.Temperature(449.271155)
    _rhor = unidades.Density(233.873368)
    _w = 0.247058753

    Fi1 = {
        "ao_log": [1, 3.0],
        "pow": [0, 1],
        "ao_pow": [],
        "ao_exp": [],
        "titao": [],
        "ao_sinh": [8.95043, 33.4032],
        "sinh": [0.380391739, 3.777411113],
        "ao_cosh": [21.836],
        "cosh": [1.789520971]
    }

    CP0 = {
        "ao": 4,
        "an": [],
        "pow": [],
        "ao_exp": [],
        "exp": [],
        "ao_sinh": [8.95043, 33.4032],
        "sinh": [178.67, 1774.25],
        "ao_cosh": [21.836],
        "cosh": [840.538]
    }

    CP1 = {
        "ao": 10.288132,
        "an": [-0.2695377e-1, 0.20951065e-3, -0.27910773e-6, 0.12266269e-9],
        "pow": [1, 2, 3, 4],
        "ao_exp": [],
        "exp": []
    }

    f = 4.184 / 8.3159524
    CP2 = {
        "ao": 22.5012 * f,
        "an": [],
        "pow": [],
        "ao_exp": [],
        "exp": [],
        "ao_sinh": [2.057417e8 * f],
        "sinh": [1.71958e3],
        "ao_cosh": [2.972927e7 * f],
        "cosh": [8.02069e2]
    }

    CP3 = {
        "ao": 4,
        "an": [],
        "pow": [],
        "ao_exp": [9.751560716, 22.71445741, 11.65392685],
        "exp": [404.8796661, 1785.491483, 4504.430788]
    }

    shortSpan = {
        "__type__":
        "Helmholtz",
        "__name__":
        "short Helmholtz equation of state for pentane of Span "
        "and Wagner (2003)",
        "__doi__": {
            "autor": "Span, R., Wagner, W.",
            "title": "Equations of state for technical applications. "
            "II. Results for nonpolar fluids.",
            "ref": "Int. J. Thermophys. 24 (1) (2003) 41-109",
            "doi": "10.1023/A:1022310214958"
        },
        "R":
        8.31451,
        "cp":
        CP0,
        "ref":
        "OTO",
        "M":
        72.15,
        "Tc":
        469.7,
        "rhoc":
        232 / 72.15,
        "Tmin":
        Tt,
        "Tmax":
        750.0,
        "Pmax":
        100000.0,
        "rhomax":
        11.2,
        "nr1": [
            0.10968643e1, -0.29988888e1, 0.99516887, -0.16170709, 0.11334460,
            0.26760595e-3
        ],
        "d1": [1, 1, 1, 2, 3, 7],
        "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875],
        "nr2": [
            0.40979882, -0.40876423e-1, -0.38169482, -0.10931957,
            -0.32073223e-1, 0.16877016e-1
        ],
        "d2": [2, 5, 1, 4, 3, 4],
        "t2": [0.625, 1.75, 3.625, 3.625, 14.5, 12.],
        "c2": [1, 1, 2, 2, 3, 3],
        "gamma2": [1] * 6
    }

    GERG = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for pentane of Kunz and "
        "Wagner (2004).",
        "__doi__": {
            "autor":
            "Kunz, O., Wagner, W.",
            "title":
            "The GERG-2008 Wide-Range Equation of State for "
            "Natural Gases and Other Mixtures: An Expansion "
            "of GERG-2004",
            "ref":
            "J. Chem.Eng. Data 57(11) (2012) 3032-3091",
            "doi":
            "10.1021/je300655b"
        },
        "R":
        8.314472,
        "cp":
        Fi1,
        "ref":
        "OTO",
        "Tmin":
        143.47,
        "Tmax":
        600.0,
        "Pmax":
        100000.0,
        "rhomax":
        10.57,
        "nr1": [
            0.10968643098001e1, -0.29988888298061e1, 0.99516886799212,
            -0.16170708558539, 0.11334460072775, 0.26760595150748e-3
        ],
        "d1": [1, 1, 1, 2, 3, 7],
        "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875],
        "nr2": [
            0.40979881986931, -0.40876423083075e-1, -0.38169482469447,
            -0.10931956843993, -0.32073223327990e-1, 0.16877016216975e-1
        ],
        "d2": [2, 5, 1, 4, 3, 4],
        "t2": [0.625, 1.75, 3.625, 3.625, 14.5, 12.],
        "c2": [1, 1, 2, 2, 3, 3],
        "gamma2": [1] * 6
    }

    polt = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for pentane of Polt (1992)",
        "__doi__": {
            "autor": "Polt, A., Platzer, B., Maurer, G.",
            "title": "Parameter der thermischen Zustandsgleichung von "
            "Bender fuer 14 mehratomige reine Stoffe",
            "ref": "Chem. Technik 22(1992)6 , 216/224",
            "doi": ""
        },
        "R":
        8.3143,
        "cp":
        CP1,
        "ref":
        "NBP",
        "Tmin":
        238.0,
        "Tmax":
        573.0,
        "Pmax":
        30000.0,
        "rhomax":
        9.410819,
        "nr1": [
            -0.117648900900e1, 0.163499095773e1, -0.366669005817,
            0.724947274043, -0.221919300269e1, 0.188671490348e1,
            -0.195774652096e1, 0.308440851184, 0.437424419722, -0.625853472351,
            .382868807091, -0.119467393955, .218631441082, 0.485668874195e-1,
            -0.132198161379, 0.213549844850e-1
        ],
        "d1": [0, 0, 0, 1, 1, 1, 1, 1, 2, 2, 2, 3, 3, 4, 4, 5],
        "t1": [3, 4, 5, 0, 1, 2, 3, 4, 0, 1, 2, 0, 1, 0, 1, 1],
        "nr2": [
            0.117648900900e1, -0.163499095773e1, 0.366669005817,
            -0.363660829618e-2, 0.633672105685, -0.705792643982
        ],
        "d2": [0, 0, 0, 2, 2, 2],
        "t2": [3, 4, 5, 3, 4, 5],
        "c2": [2] * 6,
        "gamma2": [0.968832] * 6
    }

    starling = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for pentane of Starling "
        "(1973)",
        "__doi__": {
            "autor": "Starling, K.E.",
            "title": "Fluid Thermodynamic Properties for Light "
            "Petroleum Systems",
            "ref": "Gulf Publishing Company, 1973.",
            "doi": ""
        },
        "R":
        8.3159524,
        "cp":
        CP2,
        "ref":
        "NBP",
        "Tmin":
        177.0,
        "Tmax":
        589.0,
        "Pmax":
        55000.0,
        "rhomax":
        10.2534,
        "nr1": [
            0.175873733594e1, 0.485604047435, -0.111896446456e1,
            -0.685918143315, 0.368714111378e-1, -0.167498784887e-2,
            0.327765295239, -0.352742092747, -0.999487301826e-1,
            0.781999120830e-2, 0.221577806386e-2
        ],
        "d1": [0, 1, 1, 1, 1, 1, 2, 2, 2, 5, 5],
        "t1": [3, 0, 1, 3, 4, 5, 0, 1, 2, 1, 2],
        "nr2": [-0.175873733594e1, -0.411653507564],
        "d2": [0, 2],
        "t2": [3, 3],
        "c2": [2] * 2,
        "gamma2": [0.46812392] * 2
    }

    sun = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for pentane of Sun and Ely "
        "(2004)",
        "__doi__": {
            "autor":
            "Sun, L., Ely, J.F.",
            "title":
            "Universal equation of state for engineering "
            "application: Algorithm and  application to "
            "non-polar and polar fluids",
            "ref":
            "Fluid Phase Equilib., 222-223 (2004) 107-118",
            "doi":
            "10.1016/j.fluid.2004.06.028"
        },
        "R":
        8.31451,
        "cp":
        Fi1,
        "ref":
        "OTO",
        "Tmin":
        Tt,
        "Tmax":
        620.0,
        "Pmax":
        800000.0,
        "rhomax":
        40.,
        "nr1": [
            2.20261753, 1.07797592, -3.82130221, 1.06627357e-1, 3.07513215e-4,
            -2.84309667e-1
        ],
        "d1": [1, 1, 1, 3, 7, 2],
        "t1": [1.5, 0.25, 1.25, 0.25, 0.875, 1.375],
        "nr2": [
            -7.28441220e-2, -4.60943732e-1, 8.39360011e-2, -1.50650444e-2,
            -0.203771872, -7.90244277e-3, -5.68993564e-2, -2.99387974e-2
        ],
        "d2": [1, 1, 2, 5, 1, 1, 4, 2],
        "t2": [0, 2.375, 2., 2.125, 3.5, 6.5, 4.75, 12.5],
        "c2": [1, 1, 1, 1, 2, 2, 2, 3],
        "gamma2": [1] * 8
    }

    ratanapisit = {
        "__type__":
        "MBWR",
        "__name__":
        "MBWR equation of state for pentane of Ratanapisit (1999)",
        "__doi__": {
            "autor":
            "Ratanapisit, J., Ely, J.F.",
            "title":
            "Application of New, Modified BWR Equations of "
            "State to the Corresponding-States Prediction of "
            "Natural Gas Properties",
            "ref":
            "Int. J. Thermophys., 20(6) (1999) 1721-1735",
            "doi":
            "10.1023/A:1022610013596"
        },
        "R":
        8.31451,
        "Tc":
        469.65,
        "Pc":
        3364.56,
        "rhoc":
        3.2155,
        "cp":
        CP3,
        "ref":
        "OTO",
        "Tmin":
        Tt,
        "Tmax":
        600.0,
        "Pmax":
        70000.0,
        "rhomax":
        11.2,
        "b": [
            None, -7.41533782499e-2, 7.54044021950, -1.93328401588e2,
            3.39428034054e4, -5.12571561595e6, 1.51195406963e-3,
            -7.12225059892, 4.12664185793e3, 8.40258305443e5,
            -4.68416651753e-4, 3.03565637672, -1.42146321204e3,
            -1.10170659283e-1, -9.80664356304, 1.10979804446e3, 2.98029604130,
            -1.41484307201e-1, -3.39208006239e1, 2.08782048763,
            5.38055429992e5, -6.40401885304e8, -1.19676622034e5,
            1.71973349582e10, -3.06383363882e3, 1.43168348944e6,
            1.41452433419e1, -2.52955687564e7, -3.85316416299, 2.65416349789e3,
            4.76643876980e-3, -8.37595968663, -1.35160880503e3
        ]
    }

    eq = shortSpan, GERG, polt, starling, sun, ratanapisit
    _PR = [-0.0752, -17.5]

    _surface = {
        "sigma": [0.08015, 0.004384, -0.03437],
        "exp": [1.408, 1.031, 1.818]
    }
    _dielectric = {
        "eq": 1,
        "a": [25.39, 0.025],
        "b": [78.39, 54.15],
        "c": [-12480, -4800.0],
        "Au": 29.84,
        "D": 2
    }

    _melting = {
        "eq": 2,
        "__doi__": {
            "autor": "Reeves, L.E., Scott, G.J., Babb, S.E. Jr.",
            "title": "Melting Curves of Pressure-Transmitting fluids",
            "ref": "Fluid Phase Equilib., 222-223 (2004) 107-118",
            "doi": "10.1063/1.1725068"
        },
        "Tmin": Tt,
        "Tmax": 2000.0,
        "Tref": Tt,
        "Pref": 0.076321,
        "a2": [6600e5],
        "exp2": [1.649]
    }

    _vapor_Pressure = {
        "eq": 3,
        "n": [-0.73918e1, 0.31102e1, -0.22415e1, -0.31585e1, -0.90451],
        "t": [1., 1.5, 1.74, 3.75, 8.0]
    }
    _liquid_Density = {
        "eq": 1,
        "n": [0.10178e1, 0.42703, 0.11334e1, 0.41518, -0.47950e-1],
        "t": [0.27, 0.44, 0.6, 4.0, 5.0]
    }
    _vapor_Density = {
        "eq": 2,
        "n": [-2.9389, -6.2784, -19.941, -16.709, -36.543, -127.99],
        "t": [0.4, 1.18, 3.2, 6.6, 7.0, 15.0]
    }

    visco0 = {
        "__name__": u"Quiñones-Cisneros (2006)",
        "__doi__": {
            "autor": "Quiñones-Cisneros, S.E., Deiters, U.K.",
            "title": "Generalization of the Friction Theory for "
            "Viscosity Modeling",
            "ref": "J. Phys. Chem. B, 110(25) (2006) 12820-12834",
            "doi": "10.1021/jp0618577"
        },
        "eq": 4,
        "omega": 0,
        "Toref": 469.7,
        "no": [17.6805, -55.6942, 48.7177],
        "to": [0, 0.25, 0.5],
        "a": [1.08193e-5, -4.71699e-5, 0.0],
        "b": [1.21502e-4, -9.84766e-5, 0.0],
        "c": [5.08307e-5, -1.07e-5, 0.0],
        "A": [-2.10025e-10, -1.56583e-9, 0.0],
        "B": [1.98521e-8, 2.05972e-9, 0.0],
        "C": [-1.18487e-7, 1.69571e-7, 0.0]
    }

    _viscosity = visco0,

    thermo0 = {
        "__name__":
        "Vassiliou (2015)",
        "__doi__": {
            "autor": "Vassiliou, C.-M., Assael, M.J., Huber, M.L., "
            "Perkins, R.A.",
            "title": "Reference Correlation of the Thermal Conductivity"
            " of Cyclopentane, iso-pentane, and n-Pentane",
            "ref": "J. Phys. Chem. Ref. Data 44(3) (2015) 033102",
            "doi": "10.1063/1.4927095"
        },
        "eq":
        1,
        "Toref":
        469.7,
        "koref":
        1e-3,
        "no_num": [-3.96685, 35.3805, 5.11554, -108.585, 179.573, 39.2128],
        "to_num": [0, 1, 2, 3, 4, 5],
        "no_den": [2.71636, -5.76265, 6.77885, -0.59135, 1],
        "to_den": [0, 1, 2, 3, 4],
        "Tref_res":
        469.7,
        "rhoref_res":
        232,
        "kref_res":
        1e-3,
        "nr": [
            7.76054e-1, 1.17655e2, -1.33101e2, 5.34026e1, -6.8793, 7.97696,
            -7.85888e1, 9.16089e1, -3.70431e1, 5.0962
        ],
        "tr": [0, 0, 0, 0, 0, -1, -1, -1, -1, -1],
        "dr": [1, 2, 3, 4, 5, 1, 2, 3, 4, 5],
        "critical":
        3,
        "gnu":
        0.63,
        "gamma":
        1.239,
        "R0":
        1.02,
        "Xio":
        0.227e-9,
        "gam0":
        0.058,
        "qd":
        0.668e-9,
        "Tcref":
        704.55
    }

    _thermal = thermo0,
Esempio n. 19
0
class R161(MEoS):
    """Multiparameter equation of state for R161"""
    name = "fluoroethane"
    CASNumber = "353-36-6"
    formula = "C2H5F"
    synonym = "R161"
    rhoc = unidades.Density(301.81366)
    Tc = unidades.Temperature(375.25)
    Pc = unidades.Pressure(5010.0, "kPa")
    M = 48.0595  # g/mol
    Tt = unidades.Temperature(130.0)
    Tb = unidades.Temperature(235.6)
    f_acent = 0.216
    momentoDipolar = unidades.DipoleMoment(1.9397, "Debye")
    id = 247

    Fi1 = {"ao_log": [1, 3],
           "pow": [0, 1],
           "ao_pow": [-6.9187, 5.4788],
           "ao_exp": [2.059, 9.253, 6.088],
           "titao": [420/Tc, 1548/Tc, 3882/Tc]}

    CP1 = {"ao": 3.985,
           "an": [], "pow": [],
           "ao_exp": [2.077, 9.265, 6.054], "exp": [420, 1548, 3882],
           "ao_hyp": [], "hyp": []}

    helmholtz1 = {
        "__type__": "Helmholtz",
        "__name__": "short Helmholtz equation of state for R-161 of Wu and Zhou (2012).",
        "__doi__": {"autor": "Wu, J. and Zhou, Y.",
                    "title": "An Equation of State for Fluoroethane (R161)",
                    "ref": "Int. J. Thermophys. 33:220-234, 2012.",
                    "doi": "10.1007/s10765-011-1151-3"},
        "R": 8.314472,
        "cp": Fi1,
        "ref": {"Tref": 273.15, "Pref": 1., "ho": 28559.6, "so": 167.205},

        "Tmin": Tt, "Tmax": 450.0, "Pmax": 5000.0, "rhomax": 20.0,
        "Pmin": 0.005512, "rhomin": 19.91,

        "nr1": [1.511, -2.3 , -0.457, 0.1683, 0.04133],
        "d1": [1, 1, 2, 3, 4],
        "t1": [0.37, 0.97, 1.14, 0.744, 1.],

        "nr2": [0.62187, -0.0265, -1.03, -0.285, -0.476],
        "d2": [2, 7, 1, 2, 3],
        "t2": [1.26, 1., 1.8, 3., 2.25],
        "c2": [1, 1, 2, 2, 2],
        "gamma2": [1]*5,

        "nr3": [0.82, -0.3532, -0.116, -0.0220583, -1.63148],
        "d3": [11, 1, 3, 3, 3],
        "t3": [1, 1.2, 5.3, 1, 4],
        "alfa3": [0.96, 1.35, 1.26, 1.23, 16.8],
        "beta3": [2.7, 5.2, 3.9, 4.7, 413],
        "gamma3": [0.9, 0.69, 0.67, 0.67, 1.15],
        "epsilon3": [0.683, 0.892, 0.785, 1.33, 0.86]}

    helmholtz2 = {
        "__type__": "Helmholtz",
        "__name__": "short Helmholtz equation of state for R-161 of Lemmon (2005).",
        "__doi__": {"autor": "Lemmon, E.W.",
                    "title": "preliminary equation, 2005.",
                    "ref": "",
                    "doi": ""},
        "R": 8.314472,
        "cp": CP1,

        "Tmin": Tt, "Tmax": 400.0, "Pmax": 50000.0, "rhomax": 20.0,
        "Pmin": 0.006, "rhomin": 19.95,

        "nr1": [0.75688, -1.4110, -0.63922, 0.055685, 0.00028395],
        "d1": [1, 1, 1, 3, 7],
        "t1": [0.25, 1.25, 1.5, 0.25, 0.875],

        "nr2": [.73357, .67596, .011369, -.56406, -.094362, -.1678, .00034215],
        "d2": [1, 2, 5, 1, 1, 4, 2],
        "t2": [2.375, 2.0, 2.125, 3.5, 6.5, 4.75, 12.5],
        "c2": [1, 1, 1, 2, 2, 2, 2],
        "gamma2": [1]*7}

    eq = helmholtz1, helmholtz2

    _surface = {"sigma": [0.05385], "exp": [1.111]}
    _vapor_Pressure = {
        "eq": 5,
        "ao": [-0.75224e1, 0.29140e1, -0.30129e1, -0.44497e1, 0.24207e1],
        "exp": [1.0, 1.5, 2.3, 6.0, 7.0]}
    _liquid_Density = {
        "eq": 1,
        "ao": [-.22587e2, .13424e3, -.2671e3, .3389e3, -.31059e3, .13009e3],
        "exp": [0.56, 0.7, 0.9, 1.2, 1.5, 1.7]}
    _vapor_Density = {
        "eq": 3,
        "ao": [-0.62548e1, 0.10499e2, -0.20353e2, -0.36709e2, -0.86781e2],
        "exp": [0.56, 1.3, 1.7, 5.0, 11.0]}
Esempio n. 20
0
class R14(MEoS):
    """Multiparameter equation of state for R14"""
    name = "tetrafluoromethane"
    CASNumber = "75-73-0"
    formula = "CF4"
    synonym = "R14"
    rhoc = unidades.Density(625.66)
    Tc = unidades.Temperature(227.51)
    Pc = unidades.Pressure(3750.0, "kPa")
    M = 88.0046  # g/mol
    Tt = unidades.Temperature(89.54)
    Tb = unidades.Temperature(145.10)
    f_acent = 0.1785
    momentoDipolar = unidades.DipoleMoment(0.0, "Debye")
    id = 218

    CP1 = {"ao": 3.9465247,
           "an": [-.88586725e-2, 0.13939626e-3, -0.30056204e-6, 0.20504001e-9],
           "pow": [1, 2, 3, 4],
           "ao_exp": [], "exp": [],
           "ao_hyp": [], "hyp": []}

    helmholtz1 = {
        "__type__": "Helmholtz",
        "__name__": "Helmholtz equation of state for R-14 of Platzer et al. (1990)",
        "__doi__": {"autor": "Platzer, B., Polt, A., and Maurer, G.",
                    "title": "Thermophysical properties of refrigerants", 
                    "ref": "Berlin:  Springer-Verlag, 1990.",
                    "doi": ""}, 
        "R": 8.31451,
        "cp": CP1,
        "ref": "NBP", 
        
        "Tmin": 120.0, "Tmax": 623.0, "Pmax": 51000.0, "rhomax": 20.764, 
        "Pmin": 0.64144, "rhomin": 20.764, 

        "nr1": [-.334698748966, .586690904687, -.147068929692, .103999039623e1,
                -.245792025288e1, .799614557889, -.749498954929, .152177772502,
                -.293408331764, .717794502866, -.426467444199e-1, .226562749365,
                -0.391091694003, -0.257394804936e-1, 0.554844884782e-1,
                0.610988261204e-2],
        "d1": [0, 0, 0, 1, 1, 1, 1, 1, 2, 2, 2, 3, 3, 4, 4, 5],
        "t1": [3, 4, 5, 0, 1, 2, 3, 4, 0, 1, 2, 0, 1, 0, 1, 1],

        "nr2": [.334698748966, -.586690904687, .147068929692, -.190315426142,
                .716157133959, -.703161904626],
        "d2": [0, 0, 0, 2, 2, 2],
        "t2": [3, 4, 5, 3, 4, 5],
        "c2": [2]*6,
        "gamma2": [0.99832625]*6}

    eq = helmholtz1,

    _surface = {"sigma": [0.0423], "exp": [1.24]}
    _vapor_Pressure = {
        "eq": 5,
        "ao": [-0.61905e1, -0.91398e1, 0.12192e2, -0.47215e1, -0.20439e1],
        "exp": [1.0, 1.5, 1.64, 2.5, 7.3]}
    _liquid_Density = {
        "eq": 1,
        "ao": [-0.10612e1, 0.44343e1, -0.38753e1, 0.29825e1, 0.30746],
        "exp": [0.1, 0.24, 0.4, 0.6, 3.9]}
    _vapor_Density = {
        "eq": 3,
        "ao": [-.55804e2, .10868e3, -.64257e2, -.11954e4, .36688e4, -.25956e4],
        "exp": [0.713, 0.84, 1.0, 5.8, 6.3, 6.6]}
Esempio n. 21
0
class MM(MEoS):
    """Multiparameter equation of state for hexamethyldisiloxane"""
    name = "hexamethyldisiloxane"
    CASNumber = "107-46-0"
    formula = "C6H18OSi2"
    synonym = "MM"
    _refPropName = "MM"
    _coolPropName = "MM"
    rhoc = unidades.Density(304.4043888253152)
    Tc = unidades.Temperature(518.69997204)
    Pc = unidades.Pressure(1939.39, "kPa")
    M = 162.3768  # g/mol
    Tt = unidades.Temperature(204.93)
    Tb = unidades.Temperature(373.401)
    f_acent = 0.418
    momentoDipolar = unidades.DipoleMoment(0.801, "Debye")
    id = 1376

    Fi1 = {
        "ao_log": [1, 3],
        "pow": [0, 1],
        "ao_pow": [72.110754, -10.431499],
        "ao_exp": [18.59, 29.58, 19.74, 4.87],
        "titao": [20 / Tc, 1400 / Tc, 3600 / Tc, 6300 / Tc]
    }

    CP1 = {
        "ao": -51.894,
        "an": [741.34e-3, -416e-6, 70e-9],
        "pow": [1, 2, 3],
        "ao_exp": [],
        "exp": [],
        "ao_hyp": [],
        "hyp": []
    }

    thol = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for hexamethyldisiloxane of "
        "Thol (2015).",
        "__doi__": {
            "autor":
            "Thol, M., Dubberke, F.H., Rutkai, G., Windmann, "
            "T., Köster, A., Span, R., Vrabec, J.",
            "title":
            "Fundamental equation of state correlation for "
            "hexamethyldisiloxane based on experimental and "
            "molecular simulation data",
            # TODO: Search reference
            "ref":
            "Fluid Phase Equilibria, ",
            "doi":
            "10.1016/j.fluid.2015.09.047"
        },
        "R":
        8.3144621,
        "cp":
        Fi1,
        "ref":
        "OTO",
        "Tmin":
        220.0,
        "Tmax":
        1200.0,
        "Pmax":
        600000.0,
        "rhomax":
        5.266,
        "M":
        162.3768,
        "Tc":
        518.7,
        "rhoc":
        1.653,
        "Pc":
        1931.1,
        "nr1": [0.5063651e-1, 0.8604724, -0.9179684, -0.1146325, 0.4878559],
        "d1": [4, 1, 1, 2, 3],
        "t1": [1, 0.346, 0.46, 1.01, 0.59],
        "nr2": [-0.2434088, -0.1621326, 0.6239872, -0.2306057, -0.5555096e-1],
        "d2": [1, 3, 2, 2, 7],
        "t2": [2.600, 3.330, 0.750, 2.950, 0.930],
        "c2": [2, 2, 1, 2, 1],
        "gamma2": [1] * 5,
        "nr3": [
            0.9385015, -0.2493508, -0.3308032, -0.1885803, -0.9883865e-1,
            0.1111090, 0.1061928, -0.1452454e-1
        ],
        "d3": [1, 1, 3, 3, 1, 2, 3, 1],
        "t3": [1.33, 1.68, 1.7, 3.08, 5.41, 1.4, 1.1, 5.3],
        "alfa3": [1.0334, 1.544, 1.113, 1.113, 1.11, 7.2, 1.45, 4.73],
        "beta3": [0.4707, 0.32, 0.404, 0.517, 0.432, 7.2, 1.2, 35.8],
        "gamma3": [1.7754, 0.692, 1.242, 0.421, 0.406, 0.163, 0.795, 0.88],
        "epsilon3": [0.8927, 0.5957, 0.559, 1.056, 1.3, 0.106, 0.181, 0.525],
        "nr4": []
    }

    colonna = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for hexamethyldisiloxane of "
        "Colonna (2006).",
        "__doi__": {
            "autor": "Colonna, P., Nannan, N.R., Guardone, A., "
            "Lemmon, E.W.",
            "title": "Multiparameter Equations of State for Selected "
            "Siloxanes",
            "ref": "Fluid Phase Equilibria, 244:193-211, 2006.",
            "doi": "10.1016/j.fluid.2006.04.015"
        },
        "R":
        8.314472,
        "cp":
        CP1,
        "ref":
        "NBP",
        "Tmin":
        273.0,
        "Tmax":
        673.0,
        "Pmax":
        30000.0,
        "rhomax":
        5.21,
        "Pmin":
        0.00269,
        "rhomin":
        5.2,
        "nr1": [
            1.01686012, -2.19713029, 0.75443188, -0.68003426, 0.19082162,
            0.10530133e-2
        ],
        "d1": [1, 1, 1, 2, 3, 7],
        "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875],
        "nr2": [
            0.6284595, 0.30903042e-1, -0.83948727, -0.20262381, -0.35131597e-1,
            0.25902341e-1
        ],
        "d2": [2, 5, 1, 4, 3, 4],
        "t2": [0.625, 1.75, 3.625, 3.625, 14.5, 12.0],
        "c2": [1, 1, 2, 2, 3, 3],
        "gamma2": [1] * 6
    }

    eq = thol, colonna

    _vapor_Pressure = {
        "eq": 3,
        "n": [-0.850230e1, 0.380300e1, -0.341500e1, -0.467900e1, -0.310600e1],
        "t": [1.0, 1.5, 1.98, 3.86, 14.6]
    }
    _liquid_Density = {
        "eq": 1,
        "n": [0.4003e1, -0.6406e1, 0.115e2, -0.1004e2, 0.4e1],
        "t": [0.436, 0.827, 1.24, 1.7, 2.23]
    }
    _vapor_Density = {
        "eq": 2,
        "n": [-0.37421e1, -0.37087e2, 0.7546e2, -0.7167e2, -68.69, -178.4],
        "t": [0.428, 1.79, 2.28, 2.8, 7, 15.4]
    }
Esempio n. 22
0
class RE245cb2(MEoS):
    """Multiparameter equation of state for RE245cb2"""
    name = "methyl-pentafluoroethyl-ether"
    CASNumber = "22410-44-2"
    formula = "CF3CF2OCH3"
    synonym = "HFE-245cb2"
    rhoc = unidades.Density(499.507581544)
    Tc = unidades.Temperature(406.813)
    Pc = unidades.Pressure(2886.4, "kPa")
    M = 150.047336  # g/mol
    Tt = unidades.Temperature(250)
    Tb = unidades.Temperature(278.76)
    f_acent = 0.354
    momentoDipolar = unidades.DipoleMoment(2.785, "Debye")
    id = 671
    # id = 1817

    CP1 = {
        "ao": 10.196438,
        "an": [],
        "pow": [],
        "ao_exp": [],
        "exp": [],
        "ao_hyp": [10.214789, 10.503071, 0.98682562, 0],
        "hyp": [814, 2031, 3040, 0]
    }

    helmholtz1 = {
        "__type__": "Helmholtz",
        "__name__":
        "Helmholtz equation of state for RE245cb2 of Zhou et al. (2010)",
        "__doi__": {
            "autor": "Zhou, Y. and Lemmon, E.W.",
            "title": "preliminary equation, 2010.",
            "ref": "",
            "doi": ""
        },
        "R": 8.314472,
        "cp": CP1,
        "ref": "NBP",
        "Tmin": Tt,
        "Tmax": 500.0,
        "Pmax": 400000.0,
        "rhomax": 10.02,
        "Pmin": 8.272,
        "rhomin": 10.,
        "nr1":
        [0.41453162e-1, 0.15010352e1, -0.23142144e1, -0.471412, 0.17182],
        "d1": [4, 1, 1, 2, 3],
        "t1": [1, 0.25, 0.786, 1.32, 0.338],
        "nr2": [-0.98793, -0.392049, 0.6848583, -0.32413816, -0.2414796e-1],
        "d2": [1, 3, 2, 2, 7],
        "t2": [2.82, 2., 1., 3., 0.766],
        "c2": [2, 2, 1, 2, 1],
        "gamma2": [1] * 5,
        "nr3": [0.82792487, -0.31833343, -0.11929747, -0.65010212],
        "d3": [1, 1, 3, 3],
        "t3": [1.75, 3.5, 3.86, 2.75],
        "alfa3": [1.023, 1.384, 0.998, 6.9],
        "beta3": [1.727, 1.543, 1.075, 88],
        "gamma3": [1.1, 0.64, 0.5, 1.26],
        "epsilon3": [0.713, 0.917, 0.69, 0.743]
    }

    eq = helmholtz1,

    _vapor_Pressure = {
        "eq": 5,
        "ao": [-7.8026, 1.8804, -2.8375, -4.3077],
        "exp": [1, 1.5, 2.5, 5]
    }
    _liquid_Density = {
        "eq": 1,
        "ao": [1.8378, 2.5311, -7.084, 18.678, -30.228, 22.985],
        "exp": [0.32, 1.08, 1.9, 2.8, 3.8, 4.9]
    }
    _vapor_Density = {
        "eq": 3,
        "ao": [-1.5224, -5.7245, -15.972, -50.473, -6.8916],
        "exp": [0.286, 0.82, 2.5, 5.6, 7.3]
    }

    visco0 = {
        "eq": 5,
        "omega": 3,
        "__doi__": {
            "autor": "T-H. Chung, Ajlan, M., Lee, L.L. and Starling, K.E",
            "title":
            "Generalized Multiparameter Correlation for Nonpolar and Polar Fluid Transport Properties",
            "ref": "Ind. Eng. Chem. Res., 1988, 27 (4), pp 671–679",
            "doi": "10.1021/ie00076a024"
        },
        "__name__": "Chung (1988)",
        "w": 0.354,
        "mur": 0.0,
        "k": 0.0
    }

    _viscosity = visco0,
Esempio n. 23
0
class R21(MEoS):
    """Multiparameter equation of state for R21"""
    name = "dichlorofluoromethane"
    CASNumber = "75-43-4"
    formula = "CHCl2F"
    synonym = "R21"
    rhoc = unidades.Density(526.0138)
    Tc = unidades.Temperature(451.48)
    Pc = unidades.Pressure(5181.2, "kPa")
    M = 102.9227  # g/mol
    Tt = unidades.Temperature(142.8)
    Tb = unidades.Temperature(282.01)
    f_acent = 0.2061
    momentoDipolar = unidades.DipoleMoment(1.37, "Debye")
    id = 642

    CP1 = {
        "ao":
        0.2376576 / 8.31451 * 102.92,
        "an": [
            0.12714330e-2 / 8.31451 * 102.92, 0.32413520e-6 / 8.31451 * 102.92,
            -0.24924280e-8 / 8.31451 * 102.92,
            0.17172080e-11 / 8.31451 * 102.92
        ],
        "pow": [1, 2, 3, 4],
        "ao_exp": [],
        "exp": [],
        "ao_hyp": [],
        "hyp": []
    }

    helmholtz1 = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for R-21 of Platzer et al. (1990)",
        "__doi__": {
            "autor": "Platzer, B., Polt, A., and Maurer, G.",
            "title": "Thermophysical properties of refrigerants",
            "ref": "Berlin:  Springer-Verlag, 1990.",
            "doi": ""
        },
        "R":
        8.31451,
        "cp":
        CP1,
        "ref":
        "NBP",
        "Tmin":
        200.0,
        "Tmax":
        473.19,
        "Pmax":
        137900.0,
        "rhomax":
        15.36,
        "Pmin":
        0.6828e-4,
        "rhomin":
        16.519,
        "nr1": [
            -.44386484873e2, .926505600935e1, -.551709104376, .504676623431,
            -.732431415692, -.868403860387, .146234705555, -.280576335053,
            0.864743656093, -0.270767233732e1, 0.330476390706e1,
            -0.210878239171, 0.449531449589, 0.120779813143, -0.277297953777,
            0.305441291172e-1
        ],
        "d1": [0, 0, 0, 1, 1, 1, 1, 1, 2, 2, 2, 3, 3, 4, 4, 5],
        "t1": [3, 4, 5, 0, 1, 2, 3, 4, 0, 1, 2, 0, 1, 0, 1, 1],
        "nr2": [
            0.443864848730e2, -0.926505600935e1, 0.551709104376,
            0.121128809552e1, 0.167119476587, -0.504876793028e-1
        ],
        "d2": [0, 0, 0, 2, 2, 2],
        "t2": [3, 4, 5, 3, 4, 5],
        "c2": [2] * 6,
        "gamma2": [0.07470252] * 6
    }

    eq = helmholtz1,

    _surface = {"sigma": [0.06924], "exp": [1.259]}
    _vapor_Pressure = {
        "eq": 5,
        "ao": [-0.70336e1, 0.15672e1, -0.33932e1, 0.17582e1, -0.86765e1],
        "exp": [1.0, 1.5, 3.0, 7.0, 10.0]
    }
    _liquid_Density = {
        "eq": 1,
        "ao": [0.33546, 0.18208e2, -0.26400e2, 0.10586e2],
        "exp": [0.09, 0.78, 0.92, 1.1]
    }
    _vapor_Density = {
        "eq": 3,
        "ao": [-0.38213, -0.55559e1, -0.15886e2, -0.44766e2, -0.27606e3],
        "exp": [0.09, 0.667, 2.5, 6.0, 15.0]
    }
Esempio n. 24
0
class C4F10(MEoS):
    """Multiparameter equation of state for perfluorobutane"""
    name = "perfluorobutane"
    CASNumber = "355-25-9"
    formula = "C4F10"
    synonym = ""
    rhoc = unidades.Density(599.8356)
    Tc = unidades.Temperature(386.326)
    Pc = unidades.Pressure(2323.4, "kPa")
    M = 238.03  # g/mol
    Tt = unidades.Temperature(145.0)
    Tb = unidades.Temperature(271.061)
    f_acent = 0.371
    momentoDipolar = unidades.DipoleMoment(0.0, "Debye")
    id = 693

    CP1 = {
        "ao": 2.0150709,
        "R": 8.31451,
        "an": [0.96863193e-1, -0.99797537e-4, 0.37348060e-7],
        "pow": [1, 2, 3],
        "ao_exp": [],
        "exp": [],
        "ao_hyp": [],
        "hyp": []
    }

    ecs = {
        "__type__": "ECS",
        "__name__":
        "Thermodynamic Extended Corresponding States model w/ T- and rho-dependent shape factors.",
        "__doc__":
        u"""Huber, M.L. and Ely, J.F., "A predictive extended corresponding states model for pure and mixed refrigerants including an equation of state for R134a," Int. J. Refrigeration, 17:18-31, 1994.""",
        "cp": CP1,
        "ref": R113,
        "eq": "helmholtz1",
        "R": 8.314471,
        "Tmin": 189.,
        "Tmax": 500.,
        "Pmax": 30000.0,
        "rhomax": 7.64,
        #        "Pmin": 0.61166, "rhomin": 55.497,
        "ft": [0.776042865e-2, -0.641975631],
        "ft_add": [],
        "ft_add_exp": [],
        "fd": [],
        "fd_exp": [],
        "ht": [0.278313281e-2, -0.593657910],
        "ht_add": [],
        "ht_add_exp": [],
        "hd": [-0.236093735e-2],
        "hd_exp": [1]
    }

    eq = ecs,

    _surface = {"sigma": [0.04429], "exp": [1.242]}
    _vapor_Pressure = {
        "eq": 5,
        "ao": [-0.72217e1, -0.18886e2, 0.47288e2, -0.29794e2, -0.50457e1],
        "exp": [1.0, 1.5, 1.65, 1.8, 4.8]
    }
    _liquid_Density = {
        "eq": 1,
        "ao": [0.36787e1, -0.20581e1, 0.98945, 0.60478e-1],
        "exp": [0.4, 0.6, 0.8, 2.0]
    }
    _vapor_Density = {
        "eq": 3,
        "ao": [-0.42967e1, -0.10715e2, -0.33457e2, -0.72206e2],
        "exp": [0.486, 1.7, 4.2, 8.0]
    }

    trnECS = {
        "eq": "ecs",
        "__name__": "Extended Corresponding States model",
        "__doc__":
        """Huber, M.L., Laesecke, A., and Perkins, R.A., Model for the viscosity and thermal conductivity of refrigerants, including a new correlation for the viscosity of R134a, Ind.Eng.Chem.Res. 42: 3163-3178 (2003).""",
        "ref": R113,
        "ref_eq": "helmholtz1",
        "eq_visco": "visco0",
        "eq_thermo": "thermo0",
        "f_int": [1.32e-3],
        "psi": [1.0],
        "phi": [1.0],
        "critical": 3,
        "gnu": 0.63,
        "gamma": 1.239,
        "R0": 1.03,
        "Xio": 0.194e-9,
        "gam0": 0.0496,
        "qd": 1.5e-9,
        "Tcref": 579.49
    }
Esempio n. 25
0
class NH3(MEoS):
    """Multiparameter equation of state for ammonia"""
    name = "ammonia"
    CASNumber = "7664-41-7"
    formula = "NH3"
    synonym = "R-717"
    _refPropName = "AMMONIA"
    _coolPropName = "Ammonia"
    rhoc = unidades.Density(225.)
    Tc = unidades.Temperature(405.40)
    Pc = unidades.Pressure(11333.0, "kPa")
    M = 17.03026  # g/mol
    Tt = unidades.Temperature(195.495)
    Tb = unidades.Temperature(239.823)
    f_acent = 0.25601
    momentoDipolar = unidades.DipoleMoment(1.470, "Debye")
    id = 63

    Fi1 = {
        "R": 8.314471,
        "ao_log": [1, -1],
        "pow": [0, 1, 1 / 3, -1.5, -1.75],
        "ao_pow": [-15.81502, 4.255726, 11.47434, -1.296211, 0.5706757],
        "ao_exp": [],
        "titao": [],
        "ao_hyp": [],
        "hyp": []
    }

    CP2 = {
        "ao":
        5.111814,
        "an": [
            -0.42966650e2, -0.10243792e-1, 0.38750775e-4, -0.46406097e-7,
            0.20268561e-10
        ],
        "pow": [-1.001, 1, 2, 3, 4],
        "ao_exp": [],
        "exp": [],
        "ao_hyp": [],
        "hyp": []
    }

    tillner = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for ammonia of Baehr and "
        "Tillner-Roth (1993)",
        "__doi__": {
            "autor":
            "Baehr, H.D., Tillner-Roth, R.",
            "title":
            "Thermodynamic Properties of Environmentally "
            "Acceptable Refrigerants: Equations of State and "
            "Tables for Ammonia, R22, R134a, R152a, and R123",
            "ref":
            "Springer-Verlag, Berlin, 1994.",
            "doi":
            "10.1007/978-3-642-79400-1"
        },
        "R":
        8.314471,
        "cp":
        Fi1,
        "ref":
        "IIR",
        "Tmin":
        Tt,
        "Tmax":
        700.,
        "Pmax":
        1000000.0,
        "rhomax":
        52.915,
        "Pmin":
        6.09,
        "rhomin":
        43.035,
        "nr1":
        [-0.1858814e01, 0.4554431e-1, 0.7238548, 0.1229470e-1, 0.2141882e-10],
        "d1": [1, 2, 1, 4, 15],
        "t1": [1.5, -0.5, 0.5, 1., 3.],
        "nr2": [
            -0.1430020e-1, 0.3441324, -0.2873571, 0.2352589e-4, -0.3497111e-1,
            0.1831117e-2, 0.2397852e-1, -0.4085375e-1, 0.2379275,
            -0.3548972e-1, -0.1823729, 0.2281556e-1, -0.6663444e-2,
            -0.8847486e-2, 0.2272635e-2, -0.5588655e-3
        ],
        "d2": [3, 3, 1, 8, 2, 8, 1, 1, 2, 3, 2, 4, 3, 1, 2, 4],
        "t2": [0, 3, 4, 4, 5, 5, 3, 6, 8, 8, 10, 10, 5, 7.5, 15, 30],
        "c2": [1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3],
        "gamma2": [1] * 16
    }

    ahrendts = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for ammonia of Ahrendts "
        "and Baehr (1979)",
        "__doi__": {
            "autor": "Ahrendts, J. and Baehr, H.D.",
            "title": "The Thermodynamic Properties of Ammonia",
            "ref": "VDI-Forsch., Number 596, 1979.",
            "doi": ""
        },
        "R":
        8.31434,
        "cp":
        CP2,
        "ref":
        "IIR",
        "Tmin":
        195.486,
        "Tmax":
        600.,
        "Pmax":
        400000.0,
        "rhomax":
        44.0,
        "Pmin":
        6.0339,
        "rhomin":
        43.137,
        "nr1": [
            0.911447599671, -0.382129415537e1, 0.147730246416e1,
            0.580205129871e-1, -0.574413226616e-3, 0.153018094697,
            -0.256626062036, 0.445448838055, -0.1533210545, 0.527996725202e-1,
            -0.484726581121e-1, 0.246579503330e-2, -0.107999941003e-3,
            -0.215298673010e-4, -0.306938893790e-4, 0.839163613582e-5,
            0.814833533876e-6, -0.314753664228e-7
        ],
        "d1": [1, 1, 1, 1, 1, 2, 2, 3, 4, 5, 5, 7, 9, 9, 10, 11, 12, 14],
        "t1": [1, 2, 3, 6, 9, 0, 4, 2, 1, 1, 2, 3, 3, 5, 1, 1, 5, 5],
        "nr2": [
            0.642978802435, -0.139510669941e1, 0.956135683432, -0.272787386366,
            -0.189305337334e1, 0.479043603913e1, -0.245945016980e1,
            -0.121107723958e1, 0.500552271170e1, -0.615476024667e1,
            0.210772481535e1, 0.298003513465, -0.152506723279,
            0.115565883925e-2, -0.911244657201e-3, 0.100587210000e-1,
            -0.120983155888e-1, 0.382694351151e-2
        ],
        "d2": [1, 1, 1, 1, 2, 2, 2, 3, 3, 3, 3, 4, 4, 0, 0, 0, 0, 0],
        "t2": [2, 5, 6, 7, 5, 6, 7, 3, 4, 5, 6, 6, 7, 1, 2, 0, 1, 2],
        "c2": [2] * 18,
        "gamma2":
        [0.86065403] * 13 + [506.2670781840292] * 2 + [50626.70781840292] * 3
    }

    shortSpan = {
        "__type__":
        "Helmholtz",
        "__name__":
        "short Helmholtz equation of state for ammonia of Span "
        "and Wagner (2003)",
        "__doi__": {
            "autor": "Span, R., Wagner, W.",
            "title": "Equations of State for Technical Applications. "
            "III. Results for Polar Fluids",
            "ref": "Int. J. Thermophys., 24(1) (2003) 111-162",
            "doi": "10.1023/A:1022362231796"
        },
        "R":
        8.31451,
        "cp":
        Fi1,
        "ref":
        "IIR",
        "Tmin":
        Tt,
        "Tmax":
        600.,
        "Pmax":
        100000.0,
        "rhomax":
        52.915,
        "Pmin":
        6.0531,
        "rhomin":
        43.158,
        "M":
        17.031,
        "rhoc":
        13.211203,
        "nr1": [0.7302272, -1.1879116, -0.68319136, 0.040028683, 9.0801215e-5],
        "d1": [1, 1, 1, 3, 7],
        "t1": [0.25, 1.25, 1.5, 0.25, 0.875],
        "nr2": [
            -0.056216175, 0.44935601, 0.029897121, -0.18181684, -0.09841666,
            -0.055083744, -0.88983219e-2
        ],
        "d2": [1, 2, 5, 1, 1, 4, 2],
        "t2": [2.375, 2, 2.125, 3.5, 6.5, 4.75, 12.5],
        "c2": [1, 1, 1, 2, 2, 2, 3],
        "gamma2": [1] * 7
    }

    sun = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for ammonia of Sun and Ely "
        "(2004)",
        "__doi__": {
            "autor":
            "Sun, L. and Ely, J.F.",
            "title":
            "Universal equation of state for engineering "
            "application: Algorithm and  application to "
            "non-polar and polar fluids",
            "ref":
            "Fluid Phase Equilib., 222-223 (2004) 107-118",
            "doi":
            "10.1016/j.fluid.2004.06.028"
        },
        "R":
        8.3143,
        "cp":
        Fi1,
        "ref":
        "IIR",
        "Tmin":
        Tt,
        "Tmax":
        620.0,
        "Pmax":
        800000.0,
        "rhomax":
        40.,
        "Pmin":
        0.1,
        "rhomin":
        40.,
        "nr1": [
            3.29159441e-1, 8.48237019e-1, -2.30706412, 4.08625188e-2,
            6.79597481e-5, 4.99412149e-2
        ],
        "d1": [1, 1, 1, 3, 7, 2],
        "t1": [1.5, 0.25, 1.25, 0.25, 0.875, 1.375],
        "nr2": [
            1.23624654e-1, -3.02129187e-1, 3.31747586e-1, -2.97121254e-3,
            -0.130202073, -7.45181207e-2, -4.73506171e-2, -9.70095484e-3
        ],
        "d2": [1, 1, 2, 5, 1, 1, 4, 2],
        "t2": [0, 2.375, 2., 2.125, 3.5, 6.5, 4.75, 12.5],
        "c2": [1, 1, 1, 1, 2, 2, 2, 3],
        "gamma2": [1] * 8
    }

    eq = tillner, ahrendts, shortSpan, sun

    _melting = {
        "eq": 1,
        "Tref": Tt,
        "Pref": 1000,
        "Tmin": Tt,
        "Tmax": 700.0,
        "a1": [],
        "exp1": [],
        "a2": [],
        "exp2": [],
        "a3": [0.2533125e4],
        "exp3": [1]
    }
    _surface = {"sigma": [0.1028, -0.09453], "exp": [1.211, 5.585]}
    _vapor_Pressure = {
        "eq": 3,
        "n": [-0.70993e1, -0.24330e1, 0.87591e1, -0.64091e1, -0.21185e1],
        "t": [1., 1.5, 1.7, 1.95, 4.2]
    }
    _liquid_Density = {
        "eq": 1,
        "n": [0.34488e2, -0.12849e3, 0.17382e3, -0.10699e3, 0.30339e2],
        "t": [0.58, 0.75, 0.9, 1.1, 1.3]
    }
    _vapor_Density = {
        "eq": 2,
        "n": [-0.38435, -4.0846, -6.6634, -31.881, 213.06, -246.48],
        "t": [0.218, 0.55, 1.5, 3.7, 5.5, 5.8]
    }

    visco0 = {
        "__name__":
        "Monogenidou (2018)",
        "__doi__": {
            "autor": "Monogenidou, S.A., Assael, M.J., Huber, M.L.",
            "title": "Reference Correlation for the Viscosity of Ammonia"
            " from the Triple Point to 725K and up to 50 MPa",
            "ref": "J. Phys. Chem. Ref. Data 47(2) (2018) 023102",
            "doi": "10.1063/1.5036724"
        },
        "eq":
        1,
        "omega":
        1,
        "ek":
        386,
        "sigma":
        0.2957,
        "M":
        17.03052,
        "n_chapman":
        0.021357,
        "collision": [0.39175, -0.59918, -0.00022, 0.19871, -0.06942],
        "Tref_virial":
        386,
        "n_virial": [
            -19.572881, 219.73999, -1015.3226, 2471.01251, -3375.1717,
            2491.6597, -787.26086, 14.085455, -0.34664158
        ],
        "t_virial": [0, -0.25, -0.5, -0.75, -1, -1.25, -1.5, -2.5, -5.5],
        "Tref_res":
        405.56,
        "rhoref_res":
        233.25,
        "nr": [0.0393308, 16.724735, 1.1975934, 0.0016995, -4.2399794],
        "tr": [-0.5, -0.5, -0.5, 3.5, -1.5],
        "dr": [2 / 3, 5 / 3, 14 / 3, 26 / 3, 8 / 3]
    }

    visco1 = {
        "__name__":
        "Fenghour (1995)",
        "__doi__": {
            "autor": "Fenghour, A., Wakeham, W.A., Vesovic, V., Watson, "
            "J.T.R., Millat, J., Vogel, E.",
            "title": "The Viscosity of Ammonia",
            "ref": "J. Phys. Chem. Ref. Data 24(5) (1995) 1649-1667",
            "doi": "10.1063/1.555961"
        },
        "eq":
        1,
        "omega":
        1,
        "ek":
        386,
        "sigma":
        0.2957,
        # Missing parameter of 100 in Chapman-Enskog term
        "n_chapman":
        0.021357 * 100,
        "collision": [4.9931822, -0.61122364, 0, .18535124, -0.11160946],
        "Tref_virial":
        386,
        "n_virial": [
            -1.7999496, 46.692621, -534.60794, 3360.4074, -13019.164,
            33414.230, -58711.743, 71426.686, -59834.012, 33652.741, -12027.35,
            2434.8205, -208.07957
        ],
        "t_virial":
        [0, -0.5, -1, -1.5, -2, -2.5, -3, -3.5, -4, -4.5, -5, -5.5, -6],
        "Tref_res":
        386.,
        "rhoref_res":
        M,
        "nr": [
            0.219664285, -0.083651107, 0.0017366936, -0.0064250359,
            1.67668649e-4, -1.49710093e-4, 7.7012274e-5
        ],
        "tr": [2, 4, 0, 1, 2, 3, 4],
        "dr": [2, 2, 3, 3, 4, 4, 4]
    }

    _viscosity = visco0, visco1

    thermo0 = {
        "__name__":
        "Monogenidou (2018)",
        "__doi__": {
            "autor":
            "Monogenidou, S.A., Assael, M.J., Huber, M.L.",
            "title":
            "Reference Correlations for the Thermal "
            "Conductivity of Ammonia from the Triple Point to "
            "680 K and Pressures up to 80 MPa",
            "ref":
            "J. Phys. Chem. Ref. Data 47(4) (2018) 043101",
            "doi":
            "10.1063/1.5053087"
        },
        "eq":
        1,
        "Toref":
        405.56,
        "koref":
        1e-3,
        "no_num": [86.9294, -170.5502, 608.0287, -100.9764, 85.1986],
        "to_num": [0, 1, 2, 3, 4],
        "no_den": [4.68994, 9.21307, -1.53637, 1],
        "to_den": [0, 1, 2, 3],

        # The table 2 in paper report values as mW/mK, it's a typo,
        # really is in W/mK
        "Tref_res":
        405.56,
        "rhoref_res":
        233.25,
        "kref_res":
        1,
        "nr": [
            0.103432, -0.112597, 0.233301, -0.112536, 0.141129e-1,
            -0.283976e-1, 0.482520e-1, -0.644124e-1, 0.529376e-2, 0.891203e-2
        ],
        "tr": [0, 0, 0, 0, 0, -1, -1, -1, -1, -1],
        "dr": [1, 2, 3, 4, 5, 1, 2, 3, 4, 5],
        "critical":
        3,
        "gnu":
        0.63,
        "gamma":
        1.239,
        "R0":
        1.02,
        "Xio":
        0.14e-9,
        "gam0":
        0.053,
        "qd":
        0.4e-9,
        "Tcref":
        608.34
    }

    thermo1 = {
        "__name__": "Tufeu (1984)",
        "__doi__": {
            "autor":
            "Tufeu, R., Ivanov, D.Y., Garrabos, Y., Le "
            "Neindre, B.",
            "title":
            "Thermal Conductivity of Ammonia in a Large "
            "Temperature and Pressure Range Including the "
            "Critical Region",
            "ref":
            "Ber. Bunsenges. Phys. Chem. 88 (1984) 422-427",
            "doi":
            "10.1002/bbpc.19840880421"
        },
        "eq": 1,
        "Tref": 1.,
        "kref": 1.,
        "no": [0.03589, -1.750e-4, 4.551e-7, 1.685e-10, -4.828e-13],
        "to": [0, 1, 2, 3, 4],
        "rhoref_res": 1,
        "kref_res": 1.,
        "nr": [0.16207e-3, 0.12038e-5, -0.23139e-8, 0.32749e-11],
        "tr": [0, 0, 0, 0],
        "dr": [1, 2, 3, 4],
        "critical": "_ThCondCritical"
    }

    def _ThCondCritical(self, rho, T, fase):
        # Custom Critical enhancement

        # The paper use a diferent rhoc value to the EoS
        rhoc = 235

        t = abs(T - 405.4) / 405.4
        dPT = 1e5 * (2.18 - 0.12 / exp(17.8 * t))
        nb = 1e-5 * (2.6 + 1.6 * t)

        DL = 1.2*Boltzmann*T**2/6/pi/nb/(1.34e-10/t**0.63*(1+t**0.5)) * \
            dPT**2 * 0.423e-8/t**1.24*(1+t**0.5/0.7)

        # Add correction for entire range of temperature, Eq 10
        DL *= exp(-36 * t**2)

        X = 0.61 * rhoc + 16.5 * log(t)
        if rho > 0.6 * rhoc:
            # Eq 11
            DL *= X**2 / (X**2 + (rho - 0.96 * rhoc)**2)
        else:
            # Eq 14
            DL = X**2 / (X**2 + (0.6 * rhoc - 0.96 * rhoc)**2)
            DL *= rho**2 / (0.6 * rhoc)**2

        return DL

    _thermal = thermo0, thermo1
Esempio n. 26
0
class R236ea(MEoS):
    """Multiparameter equation of state for R236ea"""
    name = "1,1,1,2,3,3-hexafluoropropane"
    CASNumber = "431-63-0"
    formula = "CF3CHFCHF2"
    synonym = "R236ea"
    rhoc = unidades.Density(565.)
    Tc = unidades.Temperature(412.44)
    Pc = unidades.Pressure(3420.0, "kPa")
    M = 152.0384  # g/mol
    Tt = unidades.Temperature(170.0)
    Tb = unidades.Temperature(279.322)
    f_acent = 0.369
    momentoDipolar = unidades.DipoleMoment(1.129, "Debye")
    id = 693

    Fi1 = {
        "ao_log": [1, 2.762],
        "pow": [0, 1],
        "ao_pow": [-14.121424135, 10.2355589225],
        "ao_exp": [0.7762, 10.41, 12.18, 3.332],
        "titao": [144 / Tc, 385 / Tc, 1536 / Tc, 7121 / Tc]
    }

    CP1 = {
        "ao": 5.30694,
        "an": [0.03973, -1.859e-5],
        "pow": [1, 2],
        "ao_exp": [],
        "exp": [],
        "ao_hyp": [],
        "hyp": []
    }

    helmholtz1 = {
        "__type__": "Helmholtz",
        "__name__":
        "Helmholtz equation of state for R236ea of Rui et al. (2013).",
        "__doi__": {
            "autor": "Rui, X., Pan, J., Wang, Y.",
            "title":
            "An Equation of State for Thermodynamic Properties of 1,1,1,2,3,3-Hexafluoropropane (R236ea)",
            "ref": "Fluid Phase Equilibria, Volume 348, 25 June 2013, Page 83",
            "doi": "10.1016/j.fluid.2012.12.026"
        },
        "R": 8.314472,
        "cp": Fi1,
        "ref": {
            "Tref": 273.15,
            "Pref": 1.,
            "ho": 56317.4970978844,
            "so": 282.8465334259
        },
        "Tmin": 240.0,
        "Tmax": 412.0,
        "Pmax": 6000.0,
        "rhomax": 10.5,
        "Pmin": 0.02,
        "rhomin": 11.7,
        "nr1": [0.051074, 2.5584, -2.9180, -0.71485, 0.15534],
        "d1": [4, 1, 1, 2, 3],
        "t1": [1., 0.264, 0.5638, 1.306, 0.2062],
        "nr2": [-1.5894, -0.784, 0.85767, -0.67235, -0.017953],
        "d2": [1, 3, 2, 2, 7],
        "t2": [2.207, 2.283, 1.373, 2.33, 0.6376],
        "c2": [2, 2, 1, 2, 1],
        "gamma2": [1] * 5,
        "nr3": [1.3165, -0.42023, -0.28053, -1.4134, -0.0000062617],
        "d3": [1, 1, 3, 3, 2],
        "t3": [1.08, 1.67, 3.502, 4.357, 0.6945],
        "alfa3": [1.019, 1.341, 1.034, 5.264, 24.44],
        "beta3": [1.3, 2.479, 1.068, 79.85, 49.06],
        "gamma3": [1.13, 0.6691, 0.465, 1.28, 0.8781],
        "epsilon3": [0.7119, 0.9102, 0.678, 0.7091, 1.727]
    }

    ecs = {
        "__type__": "ECS",
        "__name__":
        "Thermodynamic Extended Corresponding States model w/ T- and rho-dependent shape factors.",
        "__doc__":
        u"""Huber, M.L. and Ely, J.F., "A predictive extended corresponding states model for pure and mixed refrigerants including an equation of state for R134a," Int. J. Refrigeration, 17:18-31, 1994.""",
        "cp": CP1,
        "ref": R134a,
        "eq": "helmholtz1",
        # "eq": "MBWR",
        "R": 8.314471,
        "Tmin": 242.0,
        "Tmax": 500.0,
        "Pmax": 60000.0,
        "rhomax": 10.465,
        #            "Pmin": aaaaaaa, "rhomin": aaaaaaa,
        "ft": [-0.67786992, -0.52182651],
        "ft_add": [],
        "ft_add_exp": [],
        "fd": [0.113833347e-1],
        "fd_exp": [1],
        "ht": [0.142369159e1, 0.870214752e-1],
        "ht_add": [0.195298641e-1],
        "ht_add_exp": [1],
        "hd": [],
        "hd_exp": []
    }

    eq = helmholtz1, ecs

    _surface = {"sigma": [0.306974, -0.247277], "exp": [1.12614, 1.09899]}
    _vapor_Pressure = {
        "eq": 5,
        "ao": [-7.9095, 2.3374, -2.6453, -5.7058],
        "exp": [1, 1.5, 2.15, 4.75]
    }
    _liquid_Density = {
        "eq": 1,
        "ao": [1.6074, 1.5021, -1.106, 0.91146],
        "exp": [0.31, 0.75, 1.3, 1.9]
    }
    _vapor_Density = {
        "eq": 3,
        "ao": [-2.7426, -6.2268, -15.109, -49.524, -114.11],
        "exp": [0.376, 1.1, 2.7, 5.5, 11]
    }

    trnECS = {
        "eq": "ecs",
        "__name__": "Extended Corresponding States model",
        "__doc__":
        """Huber, M.L., Laesecke, A., and Perkins, R.A., Model for the viscosity and thermal conductivity of refrigerants, including a new correlation for the viscosity of R134a, Ind.Eng.Chem.Res. 42: 3163-3178 (2003).""",
        "ref": R134a,
        "ref_eq": "helmholtz1",
        "eq_visco": "visco0",
        "eq_thermo": "thermo0",
        "f_int": [1.32e-3],
        "psi": [1.0],
        "phi": [1.0],
        "critical": 3,
        "gnu": 0.63,
        "gamma": 1.239,
        "R0": 1.03,
        "Xio": 0.194e-9,
        "gam0": 0.0496,
        "qd": 1.5e-9,
        "Tcref": 579.49
    }
Esempio n. 27
0
class R13(MEoS):
    """Multiparameter equation of state for R13"""
    name = "chlorotrifluoromethane"
    CASNumber = "75-72-9"
    formula = "CClF3"
    synonym = "R13"
    _refPropName = "R13"
    _coolPropName = "R13"
    rhoc = unidades.Density(582.88)
    Tc = unidades.Temperature(302.0)
    Pc = unidades.Pressure(3879.0, "kPa")
    M = 104.459  # g/mol
    Tt = unidades.Temperature(92.0)
    Tb = unidades.Temperature(191.67)
    f_acent = 0.1723
    momentoDipolar = unidades.DipoleMoment(0.51, "Debye")
    id = 215

    # Cp⁰ coefficient for Tr term, so divide by Tc
    CP1 = {
        "ao": 1.86012334,
        "an": [8.07314520 / Tc, -1.87713639 / Tc**2, 3.17242858e-2 / Tc**3],
        "pow": [1, 2, 3]
    }

    CP2 = {
        "ao": 2.4766458,
        "an": [0.018074269, 2.1945535e-5, -8.5810657e-8, 6.3199171e-11],
        "pow": [1, 2, 3, 4]
    }

    magee = {
        "__type__":
        "MBWR",
        "__name__":
        "MBWR equation of state for R-13 of Magee et al. (2000)",
        "__doi__": {
            "autor":
            "Magee, J.W., Outcalt, S.L., Ely, J.F.",
            "title":
            "Molar Heat Capacity Cv, Vapor Pressure, and "
            "(p, ρ, T) Measurements from 92 to 350 K at "
            "Pressures to 35 MPa and a New Equation of State "
            "for Chlorotrifluoromethane (R13)",
            "ref":
            "Int. J. Thermophys., 21(5):1097-1121, 2000.",
            "doi":
            "10.1023/A:1026446004383"
        },
        "R":
        8.314471,
        "Tc":
        302,
        "Pc":
        3879.,
        "rhoc":
        5.58,
        "M":
        104.459,
        "cp":
        CP1,
        "ref":
        "IIR",
        "Tmin":
        Tt,
        "Tmax":
        403.0,
        "Pmax":
        35000.0,
        "rhomax":
        17.85,
        "b": [
            None, 0.427710490378e-2, 0.106603397093e1, -0.383065097813e2,
            0.661580211522e4, -0.800160780370e6, -0.406405755462e-2,
            0.561380767634e1, -0.247694806929e4, -0.639834580892e5,
            0.198818486764e-3, -0.206916891385, 0.749317872337e2,
            -0.431471653965e-2, 0.181741326553e1, -0.206066849491e2,
            -0.136681208829, 0.260496240940e-2, 0.287244312242,
            -0.105459756169e-1, 0.582404815872e6, -0.455721947029e8,
            0.114174177352e5, 0.265590236008e6, 0.135249873550e3,
            0.128289104267e4, 0.800900540368, -0.703307137789e4,
            0.235567665577e-2, 0.131830636112e1, -0.115187941781e-4,
            0.564530387616e-2, 0.336242130107
        ]
    }

    platzer = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Bender equation of state for R-13 of Platzer (1990)",
        "__doi__": {
            "autor": "Platzer, B., Polt, A., Maurer, G.",
            "title": "Thermophysical Properties of refrigerants",
            "ref": "Berlin:  Springer-Verlag, 1990.",
            "doi": ""
        },
        "R":
        8.31451,
        "cp":
        CP2,
        "ref":
        "IIR",
        "Tmin":
        Tt,
        "Tmax":
        450.0,
        "Pmax":
        50000.0,
        "rhomax":
        17.699806,
        "nr1": [
            -0.628346559920, 0.792797111341, -0.134038992692, 0.761143010172,
            -0.194465098795e1, 0.940938700406, -0.108107050239e1,
            0.117501564976, 0.228305167217, -0.403338888789, 0.375585713420,
            -0.617543677315e-1, 0.170326226881, 0.536612457231e-1,
            -0.151603010301, 0.252033265074e-1
        ],
        "d1": [0, 0, 0, 1, 1, 1, 1, 1, 2, 2, 2, 3, 3, 4, 4, 5],
        "t1": [3, 4, 5, 0, 1, 2, 3, 4, 0, 1, 2, 0, 1, 0, 1, 1],
        "nr2": [
            0.628346559920, -0.792797111341, 0.134038992692,
            -0.399863840975e-1, 0.436410910529, -0.448724904991
        ],
        "d2": [0, 0, 0, 2, 2, 2],
        "t2": [3, 4, 5, 3, 4, 5],
        "c2": [2] * 6,
        "gamma2": [0.98230055] * 6
    }

    eq = magee, platzer
    _PR = [-0.2395, -15.5913]

    _surface = {"sigma": [0.05045], "exp": [1.269]}
    _vapor_Pressure = {
        "eq": 3,
        "n": [-0.69311e1, 0.18281e1, -0.21901e1, -0.38177e1, 0.20803e1],
        "t": [1.0, 1.5, 2.5, 6.0, 8.0]
    }
    _liquid_Density = {
        "eq": 1,
        "n": [0.95469e1, -0.24017e2, 0.33365e2, -0.26837e2, 0.10638e2],
        "t": [0.51, 0.72, 0.94, 1.2, 1.4]
    }
    _vapor_Density = {
        "eq": 2,
        "n": [-0.31949e1, -0.73425e1, -0.21966e2, -0.51459e2, -0.85359e2],
        "t": [0.414, 1.41, 3.7, 7.7, 15.0]
    }

    trnECS = {
        "__name__": "Huber (2003)",
        "__doi__": {
            "autor":
            "Huber, M.L., Laesecke, A., Perkins, R.A.",
            "title":
            "Model for the Viscosity and Thermal Conductivity "
            "of Refrigerants, Including a New Correlation for "
            "the Viscosity of R134a",
            "ref":
            "Ind. Eng. Chem. Res., 42(13) (2003) 3163-3178",
            "doi":
            "10.1021/ie0300880"
        },
        "eq": "ecs",
        "ref": C3,
        "visco": "visco1",
        "thermo": "thermo0",
        "ek": 204.0,
        "sigma": 0.4971,
        "omega": 5,
        "psi": [0.976177, 1.48047e-2],
        "psi_d": [0, 1],
        "fint": [1.07447e-3, 6.42373e-7],
        "fint_t": [0, 1],
        "chi": [1.1394, -3.65562e-2],
        "chi_d": [0, 1],
        "critical": 3,
        "gnu": 0.63,
        "gamma": 1.239,
        "R0": 1.03,
        "Xio": 0.194e-9,
        "gam0": 0.0496,
        "qd": 3.49636e-10,
        "Tcref": 453.0
    }

    _viscosity = trnECS,
    _thermal = trnECS,
Esempio n. 28
0
    def __init__(self, T, P, mezcla):
        self.T = unidades.Temperature(T)
        self.P = unidades.Pressure(P, "atm")
        self.componente = mezcla.componente
        self.zi = mezcla.fraccion
        self.kij = Kij("BWRS")

        Aoi = []
        Boi = []
        Coi = []
        Doi = []
        Eoi = []
        ai = []
        bi = []
        ci = []
        di = []
        alfai = []
        gammai = []
        for compuesto in self.componente:
            Ao_, Bo_, Co_, Do_, Eo_, a_, b_, c_, d_, alfa_, gamma_ = self._lib(
                compuesto)
            Aoi.append(Ao_)
            Boi.append(Bo_)
            Coi.append(Co_)
            Doi.append(Do_)
            Eoi.append(Eo_)
            ai.append(a_)
            bi.append(b_)
            ci.append(c_)
            di.append(d_)
            alfai.append(alfa_)
            gammai.append(gamma_)

        Ao = Co = Do = Eo = Bo = a = b = c = d = alfa = gamma = 0
        for i in range(len(self.componente)):
            Bo += self.zi[i] * Boi[i]
            a += self.zi[i] * ai[i]**(1. / 3)
            b += self.zi[i] * bi[i]**(1. / 3)
            c += self.zi[i] * ci[i]**(1. / 3)
            d += self.zi[i] * di[i]**(1. / 3)
            alfa += self.zi[i] * alfai[i]**(1. / 3)
            gamma += self.zi[i] * gammai[i]**0.5
        a = a**3
        b = b**3
        c = c**3
        d = d**3
        alfa = alfa**3
        gamma = gamma**2

        for i in range(len(self.componente)):
            for j in range(len(self.componente)):
                Ao += self.zi[i] * self.zi[j] * Aoi[i]**0.5 * Aoi[j]**0.5 * (
                    1 - self.kij[i][j])
                Co += self.zi[i] * self.zi[j] * Coi[i]**0.5 * Coi[j]**0.5 * (
                    1 - self.kij[i][j])**3
                Do += self.zi[i] * self.zi[j] * Doi[i]**0.5 * Doi[j]**0.5 * (
                    1 - self.kij[i][j])**4
                Eo += self.zi[i] * self.zi[j] * Eoi[i]**0.5 * Eoi[j]**0.5 * (
                    1 - self.kij[i][j])**5

        self.Aoi = Aoi
        self.Boi = Boi
        self.Coi = Coi
        self.Doi = Doi
        self.Eoi = Eoi
        self.ai = ai
        self.bi = bi
        self.ci = ci
        self.di = di
        self.alfai = alfai
        self.gammai = gammai
        self.Ao = Ao
        self.Co = Co
        self.Do = Do
        self.Eo = Eo
        self.Bo = Bo
        self.a = a
        self.b = b
        self.c = c
        self.d = d
        self.alfa = alfa
        self.gamma = gamma

        Vm = lambda V: self.P.atm - R_atml * self.T / V - (
            Bo * R_atml * self.T - Ao - Co / self.T**2 + Do / self.T**3 - Eo /
            self.T**4) / V**2 - (
                b * R_atml * self.T - a - d / self.T) / V**3 - alfa * (
                    a + d / self.T) / V**6 - c / self.T**2 / V**3 * (
                        1 + gamma / V**2) * exp(-gamma / V**2)

        # Usamos SRK para estimar los volumenes de ambas fases usados como valores iniciales en la iteeración
        srk = cubic.SRK(T, P, mezcla)
        Z_srk = srk.Z
        Vgo = Z_srk[0] * R_atml * T / P
        Vlo = Z_srk[1] * R_atml * T / P
        Vg = fsolve(Vm, Vgo)
        Vl = fsolve(Vm, Vlo)
        self.V = r_[Vg, Vl]  #mol/l
        self.Z = P * self.V / R_atml / T

        self.H_exc = (
            Bo * R_atml * self.T - 2 * Ao - 4 * Co / self.T**2 +
            5 * Do / self.T**3 - 6 * Eo / self.T**4) / self.V + (
                2 * b * R_atml * self.T - 3 * a - 4 * d / self.T
            ) / 2 / self.V**2 + alfa / 5 * (
                6 * a + 7 * d / self.T) / self.V**5 + c / gamma / self.T**2 * (
                    3 - (3 + gamma / self.V**2 / 2 - gamma**2 / self.V**4) *
                    exp(-gamma / self.V**2))

        self.x, self.xi, self.yi, self.Ki = self._Flash()
Esempio n. 29
0
class R32(MEoS):
    """Multiparameter equation of state for R32"""
    name = "difluoromethane"
    CASNumber = "75-10-5"
    formula = "CH2F2"
    synonym = "R32"
    _refPropName = "R32"
    _coolPropName = "R32"
    rhoc = unidades.Density(424.)
    Tc = unidades.Temperature(351.255)
    Pc = unidades.Pressure(5782., "kPa")
    M = 52.024  # g/mol
    Tt = unidades.Temperature(136.34)
    Tb = unidades.Temperature(221.499)
    f_acent = 0.2769
    momentoDipolar = unidades.DipoleMoment(1.978, "Debye")
    id = 645

    Fi1 = {
        "R": 8.314471,
        "ao_log": [1, 3.004486],
        "pow": [0, 1],
        "ao_pow": [-8.258096, 6.353098],
        "ao_exp": [1.160761, 2.645151, 5.794987, 1.129475],
        "titao": [2.2718538, 11.9144210, 5.1415638, 32.7682170]
    }

    Fi2 = {
        "ao_log": [1, 2.999660],
        "pow": [0, 1],
        "ao_pow": [-8.253834, 6.351918],
        "ao_exp": [3.12115, 0.9994221, 2.412721, 3.055435],
        "titao": [4.559777, 2.164788, 1.234687e1, 5.877902]
    }

    CP2 = {
        "ao":
        36.79959 / 8.314471,
        "an": [
            -0.06304821 / 8.314471, 3.757936e-4 / 8.314471,
            -3.219812e-7 / 8.314471
        ],
        "pow": [1, 2, 3],
        "ao_exp": [],
        "exp": []
    }

    # Expression in tau term, dividing by Tc in all terms
    CP3 = {
        "ao": 4.3914,
        "an": [-2.5143 / 351.35, 5.3885 / 351.35**2, -1.6057 / 351.35 * 3],
        "pow": [1, 2, 3],
        "ao_exp": [],
        "exp": []
    }

    tillner = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for R-32 of Tillner-Roth "
        "and Yokozeki (1997)",
        "__doi__": {
            "autor":
            "Tillner-Roth, R., Yokozeki, A.",
            "title":
            "An International Standard Equation of State for "
            "Difluoromethane (R-32) for Temperatures from "
            "the Triple Point at 136.34 K to 435 K at "
            "Pressures up to 70 MPa",
            "ref":
            "J. Phys. Chem. Ref. Data 26(6) (1997) 1273-1328",
            "doi":
            "10.1063/1.556002"
        },
        "R":
        8.314471,
        "cp":
        Fi1,
        "ref":
        "IIR",
        "Tmin":
        Tt,
        "Tmax":
        435.0,
        "Pmax":
        70000.0,
        "rhomax":
        27.4734,
        "nr1": [
            0.1046634e1, -0.5451165, -0.2448595e-2, -0.4877002e-1,
            0.3520158e-1, 0.1622750e-2, 0.2377225e-4, 0.2914900e-1
        ],
        "d1": [1, 2, 5, 1, 1, 3, 8, 4],
        "t1": [0.25, 1., -0.25, -1., 2., 2., 0.75, 0.25],
        "nr2": [
            0.3386203e-2, -0.4202444e-2, 0.4782025e-3, -0.5504323e-2,
            -0.2418396e-1, 0.4209034, -0.4616537, -0.1200513e1, -0.2591550e1,
            -0.1400145e1, 0.8263017
        ],
        "d2": [4, 4, 8, 3, 5, 1, 1, 3, 1, 2, 3],
        "t2": [18., 26., -1., 25., 1.75, 4., 5., 1., 1.5, 1., 0.5],
        "c2": [4, 3, 1, 4, 1, 2, 2, 1, 1, 1, 1],
        "gamma2": [1] * 11
    }

    shortSpan = {
        "__type__":
        "Helmholtz",
        "__name__":
        "short Helmholtz equation of state for R-32 of Span and "
        "Wagner (2003).",
        "__doi__": {
            "autor": "Span, R., Wagner, W.",
            "title": "Equations of State for Technical Applications. "
            "III. Results for Polar Fluids",
            "ref": "Int. J. Thermophys., 24(1) (2003) 111-162",
            "doi": "10.1023/A:1022362231796"
        },
        "R":
        8.31451,
        "cp":
        Fi1,
        "ref":
        "IIR",
        "Tc":
        351.35,
        "rhoc":
        427 / M,
        "Pc":
        5795,
        "Tmin":
        Tt,
        "Tmax":
        600.0,
        "Pmax":
        100000.0,
        "rhomax":
        27.41,
        "nr1": [0.92876414, -2.4673952, 0.40129043, 0.055101049, 1.1559754e-4],
        "d1": [1, 1, 1, 3, 7],
        "t1": [0.25, 1.25, 1.5, 0.25, 0.875],
        "nr2": [
            -0.25209758, 0.42091879, 0.0037071833, -0.10308607, -0.11592089,
            -0.044350855, -0.012788805
        ],
        "d2": [1, 2, 5, 1, 1, 4, 2],
        "t2": [2.375, 2, 2.125, 3.5, 6.5, 4.75, 12.5],
        "c2": [1, 1, 1, 2, 2, 2, 3],
        "gamma2": [1] * 7
    }

    astina = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for R-32 of Astina (2003)",
        "__doi__": {
            "autor":
            "Astina, I.M., Sato, H.",
            "title":
            "A Rational Helmholtz Fundamental Equation of "
            "State for Difluoromethane with an Intermolecular"
            " Potential Background",
            "ref":
            "Int. J. Thermophys. 24(4) (2003) 963-990",
            "doi":
            "10.1023/A:1025096716493"
        },
        "R":
        8.314472,
        "cp":
        Fi2,
        "ref":
        "IIR",
        "Tmin":
        Tt,
        "Tmax":
        450.0,
        "Pmax":
        72000.0,
        "rhomax":
        27.48,
        "nr1": [
            2.118688, -4.531096, 1.442456, 2.053906e-1, -1.311675e-1,
            1.022272e-2
        ],
        "d1": [1, 1, 1, 3, 3, 4],
        "t1": [0.5, 1.125, 1.625, 0.875, 1.5, 1.75],
        "nr2": [
            4.873982e-1, -1.062213, -4.542051e-3, -6.933347e-4, -3.510307e-2,
            -5.606161e-2, 8.849625e-2, -1.850758e-2, 7.878071e-3, -3.384115e-2,
            1.641979e-4, -1.459172e-3
        ],
        "d2": [1, 1, 5, 5, 6, 1, 2, 5, 6, 2, 2, 8],
        "t2": [1.75, 2.75, 0.25, 3.75, 1, 6.5, 2.5, 7.5, 7.5, 11, 16, 13],
        "c2": [1, 1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3],
        "gamma2": [1] * 12
    }

    outcalt = {
        "__type__":
        "MBWR",
        "__name__":
        "MBWR equation of state for R-32 of Outcalt and McLinden "
        "(1995)",
        "__doi__": {
            "autor":
            "Outcalt, S.L., McLinden, M.O.",
            "title":
            "Equations of State for the Thermodynamic "
            "Properties of R32 (Difluoromethane) and R125 "
            "(Pentafluoroethane)",
            "ref":
            "Int. J. Thermophysics 16(1) (1995) 79-89.",
            "doi":
            "10.1007/BF01438959"
        },
        "R":
        8.314471,
        "Tc":
        351.35,
        "Pc":
        5795,
        "rhoc":
        8.2078,
        "cp":
        CP2,
        "ref":
        "IIR",
        "Tmin":
        Tt,
        "Tmax":
        500.0,
        "Pmax":
        60000.0,
        "rhomax":
        27.48,
        "b": [
            None, -0.131275405202e-3, 0.899927934911, -0.281400805178e2,
            0.436091182784e4, -0.837235280004e6, -0.782176408963e-6,
            -0.111226606825e1, 0.539331431878e3, 0.288600276863e6,
            -0.352264609289e-4, 0.189661830119, -0.686549003993e2,
            -0.349007064245e-2, -0.749983559476e-1, -0.321524283063e2,
            0.913057921906e-2, -0.171082181849e-3, 0.503986984347e-1,
            -0.830354867752e-3, -0.245522676708e6, -0.107859056038e8,
            -0.429514279646e4, 0.808724729567e8, -0.125945229993e2,
            -0.105735009761e4, -0.904064745354e-1, -0.183578733048e4,
            -0.169690612464e-3, 0.639250820631e-1, -0.204925767440e-6,
            -0.165629700870e-3, -0.932607493424e-2
        ]
    }

    sun = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for R-32 of Sun and Ely "
        "(2004)",
        "__doi__": {
            "autor":
            "Sun, L., Ely, J.F.",
            "title":
            "Universal equation of state for engineering "
            "application: Algorithm and  application to "
            "non-polar and polar fluids",
            "ref":
            "Fluid Phase Equilib., 222-223 (2004) 107-118",
            "doi":
            "10.1016/j.fluid.2004.06.028"
        },
        "R":
        8.314471,
        "cp":
        Fi1,
        "ref":
        "IIR",
        "Tmin":
        Tt,
        "Tmax":
        620.0,
        "Pmax":
        800000.0,
        "rhomax":
        40.,
        "nr1": [
            2.75866232e-1, 9.26526641e-1, -2.44296579, 5.34289357e-2,
            1.06739638e-4, 3.46487335e-2
        ],
        "d1": [1, 1, 1, 3, 7, 2],
        "t1": [1.5, 0.25, 1.25, 0.25, 0.875, 1.375],
        "nr2": [
            9.07435007e-2, -1.93104843e-1, 5.11370826e-1, 3.09453923e-3,
            -1.53328967e-1, -1.03816916e-1, -3.8066998e-2, -1.16075825e-2
        ],
        "d2": [1, 1, 2, 5, 1, 1, 4, 2],
        "t2": [0, 2.375, 2., 2.125, 3.5, 6.5, 4.75, 12.5],
        "c2": [1, 1, 1, 1, 2, 2, 2, 3],
        "gamma2": [1] * 8
    }

    vasserman = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for R-32 of Vasserman and"
        "Fominsky (2001)",
        "__doi__": {
            "autor": "Vasserman A.A., Fominsky D.V.",
            "title": "Equations of State for the Ozone-Safe "
            "Refrigerants R32 and R125",
            "ref": "Int. J. Thermophysics 22(4) (2001) 1089-1098",
            "doi": "10.1023/a_1010699806169"
        },
        "R":
        0.159821 * M,
        "Tc":
        351.35,
        "rhoc":
        427 / M,
        "cp":
        CP3,
        "ref":
        "IIR",
        "Tmin":
        Tt,
        "Tmax":
        620.0,
        "Pmax":
        800000.0,
        "rhomax":
        40.,
        "nr1": [
            1.183486, -2.430934, -1.472179e-2, -4.506743e-1, 1.721527,
            -1.349166, -6.052212e-1, 9.265910e-1, 8.081905e-2, -1.999587e-1,
            3.655934e-3, 8.217181e-3, -3.230880e-3, 5.778584e-3, -2.536027e-6
        ],
        "d1": [1, 1, 1, 2, 2, 2, 3, 3, 4, 4, 4, 5, 6, 6, 10],
        "t1": [0, 1, 4, 0, 1, 2, 1, 2, 0, 2, 3, 2, 0, 1, 0],
        "nr2": [
            -6.546357e-2, -2.784785e-1, 1.113400, -2.954417, 4.898234,
            -2.354906, -7.709682e-1, 6.502963e-1, 2.168338e-1, -5.499117e-1,
            1.978099e-2, 9.535163e-2, -1.425744e-2, 3.921874e-3
        ],
        "d2": [1, 1, 2, 2, 2, 2, 3, 4, 5, 5, 6, 6, 8, 9],
        "t2": [4, 5, 1, 2, 4, 5, 5, 5, 3, 4, 3, 5, 4, 2],
        "c2": [2] * 14,
        "gamma2": [1] * 14
    }

    eq = tillner, outcalt, shortSpan, astina, vasserman, sun
    _PR = [0.18086, -25.5000]

    _surface = {"sigma": [0.07147], "exp": [1.246]}

    _vapor_Pressure = {
        "eq": 3,
        "n": [-0.74883e1, 0.19697e1, -0.17496e1, -0.40224e1, 0.15209e1],
        "t": [1.0, 1.5, 2.2, 4.8, 6.2]
    }
    _liquid_Density = {
        "eq": 1,
        "n": [0.12584e1, 0.46410e1, -0.54870e1, 0.33115e1, -0.61370],
        "t": [0.27, 0.8, 1.1, 1.5, 1.8]
    }
    _vapor_Density = {
        "eq": 2,
        "n": [-.22002e1, -.5972e1, -.14571e2, -.42598e2, .42686e1, -.73373e2],
        "t": [0.336, 0.98, 2.7, 5.7, 6.5, 11.0]
    }

    trnECS = {
        "__name__": "Huber (2003)",
        "__doi__": {
            "autor":
            "Huber, M.L., Laesecke, A., Perkins, R.A.",
            "title":
            "Model for the Viscosity and Thermal Conductivity "
            "of Refrigerants, Including a New Correlation for "
            "the Viscosity of R134a",
            "ref":
            "Ind. Eng. Chem. Res., 42(13) (2003) 3163-3178",
            "doi":
            "10.1021/ie0300880"
        },
        "eq": "ecs",
        "ref": C3,
        "visco": "visco1",
        "thermo": "thermo0",
        "ek": 289.65,
        "sigma": 0.4098,
        "omega": 5,
        "psi": [0.7954, 5.42658e-2],
        "psi_d": [0, 1],
        "fint": [4.36654e-4, 1.78134e-6],
        "fint_t": [0, 1],
        "chi": [1.2942, -9.24549e-2],
        "chi_d": [0, 1],
        "critical": 3,
        "gnu": 0.63,
        "gamma": 1.239,
        "R0": 1.03,
        "Xio": 0.194e-9,
        "gam0": 0.0496,
        "qd": 5e-10,
        "Tcref": 1.5 * Tc
    }

    _viscosity = trnECS,
    _thermal = trnECS,
Esempio n. 30
0
class Toluene(MEoS):
    """Multiparameter equation of state for toluene
    """
    name = "toluene"
    CASNumber = "108-88-3"
    formula = "C6H5-CH3"
    synonym = ""
    rhoc = unidades.Density(291.98665)
    Tc = unidades.Temperature(591.75)
    Pc = unidades.Pressure(4126.3, "kPa")
    M = 92.13842  # g/mol
    Tt = unidades.Temperature(178.0)
    Tb = unidades.Temperature(383.75)
    f_acent = 0.2657
    momentoDipolar = unidades.DipoleMoment(0.36, "Debye")
    id = 41

    Fi1 = {
        "ao_log": [1, 3.],
        "pow": [0, 1],
        "ao_pow": [3.5241174832, 1.1360823464],
        "ao_exp": [1.6994, 8.0577, 17.059, 8.4567, 8.6423],
        "titao": [190 / Tc, 797 / Tc, 1619 / Tc, 3072 / Tc, 7915 / Tc]
    }

    CP1 = {
        "ao": 4.,
        "an": [],
        "pow": [],
        "ao_exp": [1.6994, 8.0577, 17.059, 8.4567, 8.6423],
        "exp": [190, 797, 1619, 3072, 7915],
        "ao_hyp": [],
        "hyp": []
    }

    CP2 = {
        "ao":
        -0.321892 / 8.3143 * 92.142,
        "an": [
            0.579338e-2 / 8.3143 * 92.142, -0.348446e-5 / 8.3143 * 92.142,
            0.143577e-8 / 8.3143 * 92.142, -0.71935e-12 / 8.3143 * 92.142
        ],
        "pow": [1, 2, 3, 4],
        "ao_exp": [],
        "exp": [],
        "ao_hyp": [],
        "hyp": []
    }

    helmholtz1 = {
        "__type__": "Helmholtz",
        "__name__":
        "short Helmholtz equation of state for toluene of Lemmon and Span (2006)",
        "__doi__": {
            "autor": "Lemmon, E.W., Span, R.",
            "title":
            "Short Fundamental Equations of State for 20 Industrial Fluids",
            "ref": "J. Chem. Eng. Data, 2006, 51 (3), pp 785–850",
            "doi": "10.1021/je050186n"
        },
        "__test__": """
            >>> st=Toluene(T=593, rho=3*92.13842)
            >>> print "%0.0f %0.0f %0.3f %0.3f %0.3f %0.3f %0.3f %0.3f" % (st.T, st.rhoM, st.P.kPa, st.hM.kJkmol, st.sM.kJkmolK, st.cvM.kJkmolK, st.cpM.kJkmolK, st.w)
            593 3 4186.620 52937.550 105.422 214.488 7705.724 89.464
            """,  # Table 10, Pag 842
        "R": 8.314472,
        "cp": Fi1,
        "ref": "NBP",
        "Tmin": Tt,
        "Tmax": 700.0,
        "Pmax": 500000.0,
        "rhomax": 10.581,
        "Pmin": 0.000039,
        "rhomin": 10.58,
        "nr1": [0.96464, -2.7855, 0.86712, -0.18860, 0.11804, 0.00025181],
        "d1": [1, 1, 1, 2, 3, 7],
        "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875],
        "nr2": [0.57196, -0.029287, -0.43351, -0.12540, -0.028207, 0.014076],
        "d2": [2, 5, 1, 4, 3, 4],
        "t2": [0.625, 1.75, 3.625, 3.625, 14.5, 12],
        "c2": [1, 1, 2, 2, 3, 3],
        "gamma2": [1] * 6
    }

    helmholtz2 = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for toluene of Polt et al. (1992)",
        "__doi__": {
            "autor": "Polt, A., Platzer, B., and Maurer, G.",
            "title":
            "Parameter der thermischen Zustandsgleichung von Bender fuer 14 mehratomige reine Stoffe",
            "ref": "Chem. Technik 22(1992)6 , 216/224",
            "doi": ""
        },
        "R":
        8.3143,
        "cp":
        Fi1,
        "ref":
        "NBP",
        "Tmin":
        298.0,
        "Tmax":
        673.0,
        "Pmax":
        25000.0,
        "rhomax":
        9.7242,
        "Pmin":
        3.774,
        "rhomin":
        9.3606,
        "nr1": [
            -0.343905499875, 0.737562743137, -0.158601557810, 0.113243121503e1,
            -0.253681929563e1, 0.104584338973e1, -0.115732119380e1,
            0.176205273278, -0.242942016719, 0.398925293195, 0.193881828889,
            0.199426230143, -0.306598708746, -0.114697533947e-1,
            0.230068676459e-1, 0.658341220591e-2
        ],
        "d1": [0, 0, 0, 1, 1, 1, 1, 1, 2, 2, 2, 3, 3, 4, 4, 5],
        "t1": [3, 4, 5, 0, 1, 2, 3, 4, 0, 1, 2, 0, 1, 0, 1, 1],
        "nr2": [
            0.343905499875, -0.737562743137, 0.15860155781, 0.40707928397,
            -0.68140614165, 0.110425925004
        ],
        "d2": [0, 0, 0, 2, 2, 2],
        "t2": [3, 4, 5, 3, 4, 5],
        "c2": [2] * 6,
        "gamma2": [0.841] * 6
    }

    helmholtz3 = {
        "__type__":
        "Helmholtz",
        "__name__":
        "Helmholtz equation of state for toluene of Sun and Ely (2004)",
        "__doi__": {
            "autor": "Sun, L. and Ely, J.F.",
            "title":
            "Universal equation of state for engineering application: Algorithm and  application to non-polar and polar fluids",
            "ref": "Fluid Phase Equilib., 222-223:107-118, 2004.",
            "doi": "10.1016/j.fluid.2004.06.028"
        },
        "R":
        8.314472,
        "cp":
        Fi1,
        "ref":
        "NBP",
        "Tmin":
        Tt,
        "Tmax":
        620.0,
        "Pmax":
        800000.0,
        "rhomax":
        40.,
        "Pmin":
        0.1,
        "rhomin":
        40.,
        "nr1": [
            1.34060172, 1.01624262, -3.27810202, 9.69209624e-2, 2.61950176e-4,
            -1.58891991e-1
        ],
        "d1": [1, 1, 1, 3, 7, 2],
        "t1": [1.5, 0.25, 1.25, 0.25, 0.875, 1.375],
        "nr2": [
            6.28559812e-2, -8.42364946e-2, 4.49701117e-1, -1.08658876e-2,
            -3.83733669e-1, 2.21127543e-2, -9.54658223e-2, -1.77905259e-2
        ],
        "d2": [1, 1, 2, 5, 1, 1, 4, 2],
        "t2": [0, 2.375, 2., 2.125, 3.5, 6.5, 4.75, 12.5],
        "c2": [1, 1, 1, 1, 2, 2, 2, 3],
        "gamma2": [1] * 8
    }

    eq = helmholtz1, helmholtz2, helmholtz3

    _surface = {"sigma": [0.06897], "exp": [1.291]}
    _vapor_Pressure = {
        "eq": 5,
        "ao": [-7.45201, 2.03681, -1.43777, -3.51652, -1.75818],
        "exp": [1.0, 1.5, 2.13, 4.0, 12.0]
    }
    _liquid_Density = {
        "eq": 1,
        "ao": [14.0531, -32.5072, 35.1091, -16.0694, 2.38699],
        "exp": [0.54, 0.72, 0.93, 1.2, 2.0]
    }
    _vapor_Density = {
        "eq": 3,
        "ao": [-2.97587, -5.34939, -19.1781, -24.0058, -32.4034, -140.645],
        "exp": [0.425, 1.06, 3.0, 6.3, 7.0, 15.0]
    }

    visco0 = {
        "eq":
        1,
        "omega":
        1,
        "__name__":
        "Lemmon (2010",
        "__doc__":
        """Lemmon, E.W. and Laesecke, A., 2010. Unpublished preliminary equation for the viscosity of toluene.""",
        "ek":
        469.90,
        "sigma":
        0.5507,
        "Tref":
        1.,
        "rhoref":
        1. * M,
        "n_chapman":
        9.598876,
        "Tref_res":
        591.75,
        "rhoref_res":
        3.169 * M,
        "n_poly": [
            0.157560701809e2, 0.658234203776e2, -0.909162962259e2,
            -0.806740654754e2, 0.395093273404e1, 0.867277691823e-1,
            -0.928414042924e-2, 0.982264892850e-5, -0.785434913708e-3,
            0.169683455336e-7
        ],
        "t_poly": [
            -0.2843, -2.4238, -2.7667, -3.0019, -3.2869, -6.0789, -6.1564,
            -6.8541, -5.5123, -4.1175
        ],
        "d_poly": [1, 2, 2, 4, 6, 9, 11, 12, 17, 19],
        "g_poly": [0, 0, 1, 1, 1, 1, 1, 0, 1, 0],
        "c_poly": [0, 0, 1, 1, 2, 1, 1, 0, 2, 0]
    }

    _viscosity = visco0,

    thermo0 = {
        "eq":
        1,
        "__name__":
        "Lemmon (2010)",
        "__doc__":
        """Lemmon, E.W. and Laesecke, A., 2010. Unpublished preliminary equation for the thermal conductivity of toluene.""",
        "Tref":
        591.75,
        "kref":
        1e-3,
        "no": [28.96745197, -167.24996945, 180.04690463],
        "co": [1.20532335, 1.58866032, 1.71267964],
        "Trefb":
        591.75,
        "rhorefb":
        3.169,
        "krefb":
        1e-3,
        "nb": [
            -0.318905053658e1, 0.258544682121e2, -0.263059677817e2,
            -0.691196173614, 0.542428651638e-1, -0.326501347819
        ],
        "tb": [-0.53316, -0.27224, -0.09974, -5.53274, -6.84315, -0.39659],
        "db": [4, 3, 5, 7, 8, 3],
        "cb": [0, 0, 1, 2, 2, 2],
        "critical":
        3,
        "gnu":
        0.63,
        "gamma":
        1.2415,
        "R0":
        1.01,
        "Xio":
        0.33442441e-9,
        "gam0":
        0.55e-1,
        "qd":
        0.71763799e-9,
        "Tcref":
        1183.5
    }

    _thermal = thermo0,