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,
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,
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,
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]}
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,
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]
}
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]
}
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]
}
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]
}
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,
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]
}
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,
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,
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
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.]
}
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]
}
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]
}
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,
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]}
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]}
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]
}
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,
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]
}
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
}
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
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
}
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,
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()
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,
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,