class neoC5(MEoS): """Multiparameter equation of state for neopentane""" name = "neopentane" CASNumber = "463-82-1" formula = "C(CH3)4" synonym = "" _refPropName = "NEOPENTN" _coolPropName = "Neopentane" rhoc = unidades.Density(235.9265106) Tc = unidades.Temperature(433.74) Pc = unidades.Pressure(3196.0, "kPa") M = 72.14878 # g/mol Tt = unidades.Temperature(256.6) Tb = unidades.Temperature(282.65) f_acent = 0.1961 momentoDipolar = unidades.DipoleMoment(0.0, "Debye") id = 9 Fi1 = { "ao_log": [1, 3.], "pow": [0, 1], "ao_pow": [0.8702452614, 1.6071746358], "ao_exp": [14.422, 12.868, 17.247, 12.663], "titao": [710 / Tc, 1725 / Tc, 3280 / Tc, 7787 / Tc] } f = 72.151 / 8.3143 CP1 = { "ao": -0.435375 * f, "an": [0.96766e-2 * f, -0.11533e-4 * f, 0.108006e-7 * f, -0.44851e-11 * f], "pow": [1, 2, 3, 4], "ao_exp": [], "exp": [] } lemmon = { "__type__": "Helmholtz", "__name__": "short Helmholtz equation of state for neopentane 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": 550.0, "Pmax": 200000.0, "rhomax": 8.71, "nr1": [1.1136, -3.1792, 1.1411, -0.10467, 0.11754, 0.00034058], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875], "nr2": [0.29553, -0.074765, -0.31474, -0.099401, -0.039569, 0.023177], "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 neopentane 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": 273.0, "Tmax": 498.0, "Pmax": 20000.0, "rhomax": 8.511, "nr1": [ -0.146552261671e1, 0.199230626557e1, -0.500821886276, 0.119809758161e1, -0.363135896710e1, 0.312770556886e1, -2.37405105853, 0.473735725047, 0.101500881659, 0.184937708516, -0.290527628579e-1, -0.258919377284e-1, 0.748831217999e-1, 0.216569936506e-1, -0.100375687935, 0.234924630013e-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.146552261671e1, -0.199230626557e1, 0.500821886276, -0.834410647812, 0.262918341468e1, -0.188136966583e1 ], "d2": [0, 0, 0, 2, 2, 2], "t2": [3, 4, 5, 3, 4, 5], "c2": [2] * 6, "gamma2": [0.968832] * 6 } eq = lemmon, polt _vapor_Pressure = { "eq": 3, "n": [-0.70262e1, 0.20090e1, -0.19932e1, -0.28503e1, -0.53760], "t": [1.0, 1.5, 2.2, 4.8, 6.2] } _liquid_Density = { "eq": 1, "n": [0.56080e1, -0.13549e2, 0.29912e2, -0.28143e2, 0.89021e1], "t": [0.45, 0.7, 1.0, 1.25, 1.6] } _vapor_Density = { "eq": 2, "n": [-.25177e1, -.63565e1, -.11985e3, .43740e3, -.10749e4, .74007e3], "t": [0.366, 1.14, 4.0, 5.0, 6.0, 6.5] }
class Novec649(MEoS): """Multiparameter equation of state for Novec649""" name = "Novec649" CASNumber = "756-13-8" formula = "C6F12O" synonym = "1,1,1,2,2,4,5,5,5-nonafluoro-4-(trifluromethyl)-3-pentanone" _refPropName = "NOVEC649" _coolPropName = "Novec649" rhoc = unidades.Density(606.805248) Tc = unidades.Temperature(441.81) Pc = unidades.Pressure(1.869, "MPa") M = 316.0444 # g/mol Tt = unidades.Temperature(165) Tb = unidades.Temperature(322.202) f_acent = 0.471 momentoDipolar = unidades.DipoleMoment(0.36, "Debye") Fi1 = { "ao_log": [1, 29.8], "pow": [0, 1], "ao_pow": [-30.6610503233, 6.8305296372], "ao_exp": [29.8], "titao": [1940 / Tc] } mclinden = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for Novec649 of McLinden " "(2015).", "__doi__": { "autor": "McLinden, M.O., Perkins, R.A., Lemmon, E.W., Fortin, " "T.J.", "title": "Thermodynamic Properties of 1,1,1,2,2,4,5,5,5-nonafluoro" "-4-(trifluoromethyl)-3-pentanone: Vapor Pressure, (p, " "ρ, T) Behavior, and Speed of Sound Measurements, and an " "Equation of State", "ref": "J. Chem. Eng. Data 60(12) (2015) 3646-3659", "doi": "10.1021/acs.jced.5b00623" }, "R": 8.3144621, "cp": Fi1, "ref": "IIR", "Tmin": Tt, "Tmax": 500.0, "Pmax": 50000.0, "rhomax": 6.24, "nr1": [ 0.05623648, 2.973616, -6.126970, 3.440240, 1.451737, -2.837857, 0.2077767 ], "d1": [4, 1, 1, 1, 2, 2, 3], "t1": [1, 0.25, 0.793, 1.16, 0.75, 1.09, 0.75], "nr2": [2.168307, -2.124648, -1.296704], "d2": [2, 1, 2], "t2": [1.3, 2.25, 1.9], "c2": [1, 2, 2], "gamma2": [1] * 3, "nr3": [ -1.010569, 2.701505, 0.8167202, -1.814579, 0.2075389, -1.009347, -0.04848043 ], "d3": [1, 1, 2, 2, 3, 3, 1], "t3": [0.88, 1.63, 1.3, 2.0, 1.15, 1.66, 1.5], "alfa3": [0.32, 1.32, 1.35, 1.48, 0.51, 1.30, 5.15], "beta3": [0.12, 0.83, 0.19, 0.95, 0.1, 0.11, 65.0], "gamma3": [1.10, 1.04, 1.15, 0.9, 0.8, 1.2, 1.19], "epsilon3": [1.16, 0.793, 1.13, 0.527, 1.19, 0.83, 0.82] } eq = mclinden, _vapor_Pressure = { "eq": 3, "n": [-8.4411, 2.711, -3.6354, -5.3872, -8.1641], "t": [1.0, 1.5, 2.2, 4.4, 15] } _liquid_Density = { "eq": 1, "n": [1.5545, 1.149, 0.51565], "t": [0.297, 0.7, 4.4] } _vapor_Density = { "eq": 2, "n": [-1.6073, -5.8095, -17.824, -61.012, -151.3], "t": [0.291, 0.82, 2.45, 5.5, 12] }
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 Ne(MEoS): """Multiparameter equation of state for neon""" name = "neon" CASNumber = "7440-01-9" formula = "Ne" synonym = "R-720" rhoc = unidades.Density(481.914888) Tc = unidades.Temperature(44.4918) Pc = unidades.Pressure(2678.6, "kPa") M = 20.179 # g/mol Tt = unidades.Temperature(24.556) Tb = unidades.Temperature(27.104) f_acent = -0.0387 momentoDipolar = unidades.DipoleMoment(0.0, "Debye") id = 107 CP1 = { "ao": 2.5, "an": [], "pow": [], "ao_exp": [], "exp": [], "ao_hyp": [], "hyp": [] } helmholtz1 = { "__type__": "Helmholtz", "__name__": u"Helmholtz equation of state for neon of Katti et al. (1986).", "__doi__": { "autor": "Katti, R.S., Jacobsen, R.T, Stewart, R.B., Jahangiri, M.", "title": "Thermodynamic properties for neon for temperatures from the triple point to 700 K at pressures up to 700 MPa", "ref": "Adv. Cryo. Eng. 31 (1986), 1189-1197", "doi": "10.1007/978-1-4613-2213-9_132" }, "R": 8.31434, "cp": CP1, "ref": { "Tref": 298.15, "Pref": 101.325, "ho": 6179, "so": 146.214 }, "Tmin": Tt, "Tmax": 700.0, "Pmax": 700000.0, "rhomax": 90.56, "Pmin": 43.464, "rhomin": 62.059, "nr1": [ 0.3532653449e1, -0.4513954384e1, -0.1524027959, 0.2188568609e1, -0.744299997e1, 0.7755627402e1, -0.3122553128e1, 0.1014206899e1, -0.5289214086e-1, 0.1566849239, -0.222852705, -0.1410150942e-1, 0.7036229719e-1, -0.5882048367e-1, 0.1571172741e-1, 0.1292202769e-2, 0.7902035603e-3, -0.3794403616e-3 ], "d1": [1, 1, 1, 2, 2, 2, 2, 2, 3, 3, 3, 4, 4, 4, 4, 6, 6, 6], "t1": [ 0.5, 0.75, 3.5, 0.5, 0.75, 1, 1.5, 2.5, 0.25, 0.5, 2.5, 1, 3, 4, 5, 1, 5, 6 ], "nr2": [ 0.4652799333e-1, 0.4524001818e-1, -0.2383421991, 0.629359013e-2, -0.1272313644e-2, -0.175235256e-6, 0.7188419232e-2, -0.5403006914e-1, 0.7578222187e-1, -0.3808588254e-1, 0.6034022431e-2 ], "d2": [1, 2, 2, 2, 2, 2, 4, 8, 8, 8, 8], "t2": [4, 1, 5, 8, 12, 32, 10, 6, 7, 8, 9], "c2": [3, 2, 2, 4, 6, 6, 2, 2, 2, 2, 2], "gamma2": [1] * 11, "nr3": [], "nr4": [] } eq = helmholtz1, _surface = { "sigma": [0.012254, 0.02728, -0.025715], "exp": [1.4136, 1.4517, 1.6567] } _dielectric = { "eq": 3, "Tref": 273.16, "rhoref": 1000., "a0": [], "expt0": [], "expd0": [], "a1": [0.9969], "expt1": [0], "expd1": [1], "a2": [-0.109, 0.0708, -2.88, -1.], "expt2": [0, 1, 0, 1], "expd2": [2, 2, 3, 3] } _melting = { "eq": 1, "Tref": Tt, "Pref": 43.36814, "Tmin": Tt, "Tmax": 700.0, "a1": [1., 4437.], "exp1": [0, 1.33], "a2": [], "exp2": [], "a3": [], "exp3": [] } _sublimation = { "eq": 3, "Tref": Tt, "Pref": 43.464, "Tmin": Tt, "Tmax": Tt, "a1": [], "exp1": [], "a2": [-10.65], "exp2": [1], "a3": [], "exp3": [] } _vapor_Pressure = { "eq": 5, "ao": [-0.55805e1, 0.68795e-1, 0.54840e1, -0.83760e1, 0.34276e1], "exp": [1.0, 1.5, 2.3, 2.8, 3.4] } _liquid_Density = { "eq": 1, "ao": [0.10601e1, 0.12076e3, -0.38553e3, 0.81655e3, -0.89907e3, 0.35466e3], "exp": [0.33, 1.4, 1.7, 2.2, 2.6, 3.0] } _vapor_Density = { "eq": 3, "ao": [-0.23338e1, -0.36834e1, -0.85368e2, 0.22769e3, -0.17290e3], "exp": [0.444, 0.95, 3.5, 4.1, 4.5] } visco0 = { "eq": 0, "method": "_visco0", "__name__": "Rabinovich (1988)", "__doi__": { "autor": "Rabinovich, V.A., Vasserman, A.A., Nedostup, V.I. and Veksler, L.S.", "title": "Thermophysical Properties of Neon, Argon, Krypton, and Xenon", "ref": "Hemisphere Publishing Corp., 1988.", "doi": "" } } _viscosity = visco0, def _visco0(self): # FIXME: Da buenos resultados, pero los resultados difierente en la tercera cifra significativa. a = [ 17.67484, -2.78751, 311498.7, -48826500, 3938774000, -1.654629e11, 2.86561e12 ] Tr = self.T / 0.29944 y = 0.68321 * (a[0] + a[1] * log10(Tr) + a[2] / Tr**2 + a[3] / Tr**3 + a[4] / Tr**4 + a[5] / Tr**5 + a[6] / Tr**6) nt = 266.93 * (self.T * self.M)**0.5 / y om = self.rho / 1673.0 c = [1.03010, -0.99175, 2.47127, -3.11864, 1.57066] b = [0.48148, -1.18732, 2.80277, -5.41058, 7.04779, -3.76608] sigma = 0.000000000305 * ( sum([ci * om**i for i, ci in enumerate(c)]) - sum([bi * om**i for i, bi in enumerate(b)]) * log10(self.T / 122.1)) br = 2.0 / 3.0 * pi * Avogadro * sigma**3 brho = self.rho / self.M * 1000 * br d = [1, 0.27676, 0.014355, 2.6480, -1.9643, 0.89161] nd = sum([di * brho**i for i, di in enumerate(d)]) return unidades.Viscosity(nd * nt / 100, "muPas")
class R404a(MEoSBlend): """Multiparameter equation of state for R404A (44% R125, 4% R134a, 52% R143a)""" name = "R404A" CASNumber = "" formula = "R125+R134a+R143a" synonym = "R404A" _refPropName = "R404A" _coolPropName = "R404A" rhoc = unidades.Density(482.162772) Tc = unidades.Temperature(345.27) Pc = unidades.Pressure(3734.8, "kPa") M = 97.6038 # g/mol Tt = unidades.Temperature(200.0) Tb = unidades.Temperature(226.93) f_acent = 0.293 momentoDipolar = unidades.DipoleMoment(0.0, "Debye") Fi1 = { "ao_log": [1, -1], "pow": [0, 1, -0.3], "ao_pow": [7.00407, 7.98695, -18.8664], "ao_exp": [0.63078, 3.5979, 5.0335], "titao": [413 / Tc, 804 / Tc, 1727 / Tc] } lemmon = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for R-404A of Lemmon (2003)", "__doi__": { "autor": "Lemmon, E.W.", "title": "Pseudo-Pure Fluid Equations of State for the " "Refrigerant Blends R-410A, R-404A, R-507A, and " "R-407C", "ref": "Int. J. Thermophys., 24(4) (2003) 991-1006", "doi": "10.1023/A:1025048800563" }, "R": 8.314472, "cp": Fi1, "ref": "IIR", "Tmin": Tt, "Tmax": 500.0, "Pmax": 50000.0, "rhomax": 14.21, "Tj": 345.270, "Pj": 3.7348, "dew": { "i": [0.1 * 2, 0.972 * 2, 3.8 * 2, 9.0 * 2], "n": [-0.00026863, -6.5757, -4.1802, -7.9102] }, "bubble": { "i": [0.54 * 2, 0.965 * 2, 3.7 * 2, 9.0 * 2], "n": [0.061067, -6.5646, -3.6162, 3.9771] }, "nr1": [ 6.10984, -7.79453, 0.0183377, 0.262270, -0.00351688, 0.0116181, 0.00105992 ], "d1": [1, 1, 1, 2, 2, 4, 6], "t1": [0.67, 0.91, 5.96, 0.7, 6, 0.3, 0.7], "nr2": [ 0.850922, -0.520084, -0.0464225, 0.62119, -0.195505, 0.336159, -0.0376062, -0.00636579, -0.0758262, -0.0221041, 0.0310441, 0.0132798, 0.0689437, -0.0507525, 0.0161382 ], "d2": [1, 1, 1, 2, 2, 3, 4, 7, 2, 3, 4, 4, 2, 3, 5], "t2": [1.7, 3.3, 7, 2.05, 4.3, 2.7, 1.8, 1.25, 12, 6, 8.7, 11.6, 13, 17, 16], "c2": [1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3], "gamma2": [1] * 15 } eq = lemmon, _surface = {"sigma": [0.06868, -0.04576], "exp": [1.3142, 2.3084]} _liquid_Density = { "eq": 1, "n": [ -0.011777381218307114, 2.4984547682678704, 5367.932468376006, -3732982.2139762687, 33474569.015986394, -80262239.93679684 ], "t": [0.066, 0.377, 8.99, 13.98, 15.982, 17.954] } _vapor_Density = { "eq": 3, "n": [ -5.211849919504925, 6.798109425092868, -13.038945243933242, 9.654798311425768, -320.35292894913766, 80023.04183470446 ], "t": [0.485, 0.989, 1.431, 2.495, 7.812, 15.329] } thermo0 = { "__name__": "Geller (2001)", "__doi__": { "autor": "Geller, V.Z., Nemzer, B.V., Cheremnykh, U.V.", "title": "Thermal Conductivity of the Refrigerant " "Mixtures R404A, R407C, R410A and R507A", "ref": "Int. J. Termophysics 22(4) (2001) 1035-1043", "doi": "10.1023/a_1010691504352" }, "eq": 1, "Toref": 1, "koref": 1e-3, "no": [-8.624, 7.36e-2], "to": [0, 1], "rhoref_res": 1, "kref_res": 1e-3, "nr": [3.222e-2, 2.569e-5, -2.693e-8, 2.007e-11], "tr": [0, 0, 0, 0], "dr": [1, 2, 3, 4], "critical": 0 } _thermal = thermo0,
class MDM(MEoS): """Multiparamenter equation of state for octamethyltrisiloxane""" name = "octamethyltrisiloxane" CASNumber = "107-51-7" formula = "C8H24O2Si3" synonym = "MDM" _refPropName = "MDM" _coolPropName = "MDM" rhoc = unidades.Density(256.73940949935815) Tc = unidades.Temperature(564.09) Pc = unidades.Pressure(1415.0, "kPa") M = 236.531 # g/mol Tt = unidades.Temperature(187.2) Tb = unidades.Temperature(425.66) f_acent = 0.529 momentoDipolar = unidades.DipoleMoment(1.079, "Debye") # id = 1893 f = 8.314472 CP1 = {"ao": 275.1/f, "an": [], "pow": [], "ao_exp": [], "exp": [], "ao_sinh": [612.9/f], "sinh": [1829.6], "ao_cosh": [413.0/f], "cosh": [802.6]} colonna = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for MDM of Colonna (2006).", "__doi__": {"autor": "Colonna, P., Nannan, N.R., 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": 3.94, "nr1": [1.19735372, -2.40380622, 0.3256564, -0.19971259, 0.11206277, 0.15893999e-3], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875], "nr2": [0.51234323, -0.20660361e-1, -0.38978114, -0.1186931, -0.37203537e-1, 0.18359984e-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.85589e1, 0.20278e1, -0.28501e1, -0.64397e1, -0.85460e1], "t": [1.0, 1.5, 2.3, 4.0, 13.0]} _liquid_Density = { "eq": 1, "n": [0.54145, -0.27650e-1, 0.41558e1, -0.19104e1, 0.67606], "t": [0.12, 0.36, 0.6, 0.8, 2.5]} _vapor_Density = { "eq": 2, "n": [-0.16483e1, -0.71410e1, -0.23088e2, -0.70554e2, 0.19938e1, -0.20193e3], "t": [0.296, 0.905, 2.8, 5.9, 12.0, 13.0]}
class Air(MEoSBlend): """Multiparameter equation of state for Air as pseudocomponent""" name = "air" CASNumber = "1" formula = "N2+Ar+O2" synonym = "R-729" rhoc = unidades.Density(342.60456) Tc = unidades.Temperature(132.6306) Pc = unidades.Pressure(3786.0, "kPa") M = 28.96546 # g/mol Tt = unidades.Temperature(59.75) Tb = unidades.Temperature(78.903) f_acent = 0.0335 momentoDipolar = unidades.DipoleMoment(0.0, "Debye") id = 475 Fi1 = { "ao_log": [1, 2.490888032], "pow": [-3, -2, -1, 0, 1, 1.5], "ao_pow": [ 0.6057194e-7, -0.210274769e-4, -0.158860716e-3, -13.841928076, 17.275266575, -0.19536342e-3 ], "ao_exp": [0.791309509, 0.212236768], "titao": [25.36365, 16.90741], "ao_exp2": [-0.197938904], "titao2": [87.31279], "sum2": [2. / 3] } helmholtz1 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for air of Lemmon et al. (2000)", "__doi__": { "autor": "Lemmon, E.W., Jacobsen, R.T, Penoncello, S.G., and Friend, D.G.", "title": "Thermodynamic Properties of Air and Mixtures of Nitrogen, Argon, and Oxygen From 60 to 2000 K at Pressures to 2000 MPa", "ref": "J. Phys. Chem. Ref. Data 29, 331 (2000)", "doi": "10.1063/1.1285884" }, "__test__": # Table A1, Pag 363 """ >>> print "%0.6f %0.5f" % (Air._bubbleP(59.75).MPa, Air._dewP(59.75).MPa) 0.005265 0.00243 >>> print "%0.5f %0.5f" % (Air._bubbleP(70).MPa, Air._dewP(70).MPa) 0.03191 0.01943 >>> print "%0.5f %0.5f" % (Air._bubbleP(80).MPa, Air._dewP(80).MPa) 0.11462 0.08232 >>> print "%0.5f %0.5f" % (Air._bubbleP(100).MPa, Air._dewP(100).MPa) 0.66313 0.56742 >>> print "%0.5f %0.5f" % (Air._bubbleP(120).MPa, Air._dewP(120).MPa) 2.15573 2.00674 >>> print "%0.5f %0.5f" % (Air._bubbleP(130).MPa, Air._dewP(130).MPa) 3.42947 3.30835 """ # Table A2, Pag 366 """ >>> st=Air(T=100, P=101325) >>> print "%0.0f %0.5g %0.5g %0.5g %0.5g %0.2f %0.2f %0.1f" % (\ st.T, st.rhoM, st.uM.kJkmol, st.hM.kJkmol, st.sM.kJkmolK, st.cvM.kJkmolK, st.cpM.kJkmolK, st.w) 100 0.12449 2028.2 2784.1 166.61 21.09 30.13 198.2 >>> st=Air(T=2000, P=101325) >>> print "%0.0f %0.4g %0.5g %0.5g %0.5g %0.2f %0.2f %0.1f" % (\ st.T, st.rhoM, st.uM.kJkmol, st.hM.kJkmol, st.sM.kJkmolK, st.cvM.kJkmolK, st.cpM.kJkmolK, st.w) 2000 0.006092 48610 65242 259.62 27.90 36.21 863.5 >>> st=Air(T=500, P=2e5) >>> print "%0.0f %0.5g %0.5g %0.5g %0.5g %0.2f %0.2f %0.1f" % (\ st.T, st.rhoM, st.uM.kJkmol, st.hM.kJkmol, st.sM.kJkmolK, st.cvM.kJkmolK, st.cpM.kJkmolK, st.w) 500 0.048077 10418 14578 208.2 21.51 29.84 446.6 >>> st=Air(T=300, P=5e5) >>> print "%0.0f %0.5g %0.5g %0.5g %0.5g %0.2f %0.2f %0.1f" % (\ st.T, st.rhoM, st.uM.kJkmol, st.hM.kJkmol, st.sM.kJkmolK, st.cvM.kJkmolK, st.cpM.kJkmolK, st.w) 300 0.20075 6179.8 8670.5 185.5 20.82 29.33 347.8 >>> st=Air(T=130, P=1e6) >>> print "%0.0f %0.5g %0.5g %0.5g %0.5g %0.2f %0.2f %0.1f" % (\ st.T, st.rhoM, st.uM.kJkmol, st.hM.kJkmol, st.sM.kJkmolK, st.cvM.kJkmolK, st.cpM.kJkmolK, st.w) 130 1.0295 2461.1 3432.5 153.79 22.058 34.69 216.8 >>> st=Air(T=70, P=5e6) >>> print "%0.0f %0.5g %0.5g %0.5g %0.5g %0.2f %0.2f %0.1f" % (\ st.T, st.rhoM, st.uM.kJkmol, st.hM.kJkmol, st.sM.kJkmolK, st.cvM.kJkmolK, st.cpM.kJkmolK, st.w) 70 31.895 -4198 -4041.2 78.907 32.17 54.57 974.6 >>> st=Air(T=2000, P=1e7) >>> print "%0.0f %0.5g %0.5g %0.5g %0.5g %0.2f %0.2f %0.1f" % (\ st.T, st.rhoM, st.uM.kJkmol, st.hM.kJkmol, st.sM.kJkmolK, st.cvM.kJkmolK, st.cpM.kJkmolK, st.w) 2000 0.59094 48600 65522 221.44 27.93 36.25 878.6 >>> st=Air(T=130, P=5e7) >>> print "%0.0f %0.5g %0.5g %0.5g %0.5g %0.2f %0.2f %0.1f" % (\ st.T, st.rhoM, st.uM.kJkmol, st.hM.kJkmol, st.sM.kJkmolK, st.cvM.kJkmolK, st.cpM.kJkmolK, st.w) 130 27.946 -1831.3 -42.096 104.84 27.68 48.19 878.8 >>> st=Air(T=100, P=1e8) >>> print "%0.0f %0.5g %0.5g %0.5g %0.5g %0.2f %0.2f %0.1f" % (\ st.T, st.rhoM, st.uM.kJkmol, st.hM.kJkmol, st.sM.kJkmolK, st.cvM.kJkmolK, st.cpM.kJkmolK, st.w) 100 33.161 -3403.9 -388.34 87.644 31.98 48.22 1192.4 >>> st=Air(T=1000, P=1e9) >>> print "%0.0f %0.5g %0.5g %0.5g %0.5g %0.2f %0.2f %0.1f" % (\ st.T, st.rhoM, st.uM.kJkmol, st.hM.kJkmol, st.sM.kJkmolK, st.cvM.kJkmolK, st.cpM.kJkmolK, st.w) 1000 30.791 21944 54421 156.83 29.07 36.77 1966.3 """, "R": 8.314472, "cp": Fi1, "ref": { "Tref": 298.15, "Pref": 101.325, "ho": 8649.34, "so": 194. }, "M": 28.9586, "Tc": 132.6312, "rhoc": 10.4477, "Tmin": Tt, "Tmax": 2000., "Pmax": 2000000.0, "rhomax": 53.73, "Pmin": 5.2646, "rhomin": 33.067, "Tj": 132.6312, "Pj": 3.78502, "dew": { "i": [1, 2, 5, 8], "n": [-0.1567266, -5.539635, 0.7567212, -3.514322] }, "bubble": { "i": [1, 2, 3, 4, 5, 6], "n": [0.2260724, -7.080499, 5.700283, -12.44017, 17.81926, -10.81364] }, "nr1": [ 0.118160747229, 0.713116392079, -0.161824192067e1, 0.714140178971e-1, -0.865421396646e-1, 0.134211176704, 0.112626704218e-1, -0.420533228842e-1, 0.349008431982e-1, 0.164957183186e-3 ], "d1": [1, 1, 1, 2, 3, 3, 4, 4, 4, 6], "t1": [0, 0.33, 1.01, 0, 0, 0.15, 0, 0.2, 0.35, 1.35], "nr2": [ -0.101365037912, -0.173813690970, -0.472103183731e-1, -0.122523554253e-1, -0.146629609713, -0.316055879821e-1, 0.233594806142e-3, 0.148287891978e-1, -0.938782884667e-2 ], "d2": [1, 3, 5, 6, 1, 3, 11, 1, 3], "t2": [1.6, 0.8, 0.95, 1.25, 3.6, 6, 3.25, 3.5, 15], "c2": [1, 1, 1, 1, 2, 2, 2, 3, 3], "gamma2": [1] * 9 } helmholtz2 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for air of Jacobsen et al. (1992)", "__doi__": { "autor": "Jacobsen, R.T, Penoncello, S.G., Beyerlein, S.W., Clarke, W.P., and Lemmon, E.W.", "title": "A Thermodynamic Property Formulation for Air", "ref": "Fluid Phase Equilibria, 79:113-124, 1992.", "doi": "10.1016/0378-3812(92)85124-Q" }, "R": 8.31451, "cp": Fi1, "ref": { "Tref": 298.15, "Pref": 101.325, "ho": 8649.34, "so": 194. }, "M": 28.9586, "Tc": 132.6312, "rhoc": 10.4477, "Tmin": Tt, "Tmax": 870.0, "Pmax": 70000.0, "rhomax": 34.628, "Pmin": 6.2545, "rhomin": 33.073, "Tj": 132.61738, "Pj": 3.78502, "dew": { "i": [1, 2, 10, 11, 13, 14], "n": [ -0.1537763029, -5.544542064, 312.7182733, -895.9553274, 1834.176566, -1321.892808 ] }, "bubble": { "i": [1, 2, 4, 5, 6, 7, 12], "n": [ 0.2095592444, -6.654905539, 22.13718815, -84.14553609, 135.9753732, -83.66895082, 17.97856602 ] }, "nr1": [ 0.206604930965, 0.367099749382, -0.943192015369, 0.382519513142e-2, -0.865385542309e-1, 0.323019987452, 0.608695449299e-2, 0.128352106296e-3, -0.400058181940e-5 ], "d1": [1, 1, 1, 1, 2, 2, 4, 6, 7], "t1": [0, 0.25, 1, 3.5, 0, 0.25, 0.5, 2, 3], "nr2": [ -0.544697915817, -0.526471065792, -0.608529300347, -0.124174277875, -0.595578533411e-2, -0.157523548353, -0.346463251040e-2, 0.837023084176e-2, -0.316701981142e-1, -0.721856676857e-2, 0.276838040645e-3, 0.160877459321e-4, 0.409235806738e-1, 0.652776125216e-3, -0.952903961290e-2, -0.100337820004e-1, 0.701111041628e-2, -0.472754336912e-2, 0.399257638569e-2, 0.968453675994e-2, -0.106826283630e-1, -0.489679885832e-2 ], "d2": [1, 2, 3, 5, 6, 1, 1, 2, 2, 3, 11, 11, 1, 1, 2, 3, 7, 8, 2, 4, 5, 2], "t2": [ 1.5, 1, 1, 1, 2, 3, 8, 0.5, 5.5, 9, 3, 6, 3, 9, 2, 13, 11, 11, 8, 22, 23, 11 ], "c2": [1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 4, 4, 4, 5], "gamma2": [1] * 22 } eq = helmholtz1, helmholtz2 _surface = {"sigma": [0.03046], "exp": [1.28]} _melting = { "eq": 1, "Tref": Tb, "Pref": 5.265, "Tmin": 59.75, "Tmax": 2000.0, "a1": [1, 0.354935e5, -0.354935e5], "exp1": [0, 0.178963e1, 0], "a2": [], "exp2": [], "a3": [], "exp3": [] } _vapor_Density = { "eq": 3, "ao": [-0.20466e1, -0.4752e1, -0.13259e2, -0.47652e2], "exp": [0.41, 1, 2.8, 6.5] } _vapor_Pressure = { "eq": 5, "ao": [-0.1567266, -0.5539635e1, 0.7567212, -0.3514322e1], "exp": [0.5, 1, 2.5, 4] } visco0 = { "eq": 1, "omega": 1, "__name__": "Lemmon (2004)", "__doi__": { "autor": "Lemmon, E.W. and Jacobsen, R.T.", "title": "Viscosity and Thermal Conductivity Equations for Nitrogen, Oxygen, Argon, and Air", "ref": "Int. J. Thermophys., 25:21-69, 2004.", "doi": "10.1023/B:IJOT.0000022327.04529.f3" }, "__test__": """ >>> st=Air(T=100, rhom=0) >>> print "%0.5f" % st.mu.muPas 7.09559 >>> st=Air(T=300, rhom=0) >>> print "%0.4f" % st.mu.muPas 18.5230 >>> st=Air(T=100, rhom=28) >>> print "%0.3f" % st.mu.muPas 107.923 >>> st=Air(T=200, rhom=10) >>> print "%0.4f" % st.mu.muPas 21.1392 >>> st=Air(T=300, rhom=5) >>> print "%0.4f" % st.mu.muPas 21.3241 >>> st=Air(T=132.64, rhom=10.4) >>> print "%0.4f" % st.mu.muPas 17.7623 """, # Table V, Pag 28 "ek": 103.3, "sigma": 0.36, "Tref": 1, "rhoref": 1. * M, "Tref_res": 132.6312, "rhoref_res": 10.4477 * M, "n_poly": [10.72, 1.122, 0.002019, -8.876, -0.02916], "t_poly": [.2, .05, 2.4, .6, 3.6], "d_poly": [1, 4, 9, 1, 8], "g_poly": [0, 0, 0, 1, 1], "c_poly": [0, 0, 0, 1, 1] } _viscosity = visco0, thermo0 = { "eq": 1, "__name__": "Lemmon (2004)", "__doi__": { "autor": "Lemmon, E.W. and Jacobsen, R.T.", "title": "Viscosity and Thermal Conductivity Equations for Nitrogen, Oxygen, Argon, and Air", "ref": "Int. J. Thermophys., 25:21-69, 2004.", "doi": "10.1023/B:IJOT.0000022327.04529.f3" }, "__test__": """ >>> st=Air(T=100, rhom=0) >>> print "%0.5f" % st.k.mWmK 9.35902 >>> st=Air(T=300, rhom=0) >>> print "%0.4f" % st.k.mWmK 26.3529 >>> st=Air(T=100, rhom=28) >>> print "%0.3f" % st.k.mWmK 119.221 >>> st=Air(T=200, rhom=10) >>> print "%0.4f" % st.k.mWmK 35.3185 >>> st=Air(T=300, rhom=5) >>> print "%0.4f" % st.k.mWmK 32.6062 >>> st=Air(T=132.64, rhom=10.4) >>> print "%0.4f" % st.k.mWmK 75.6231 """, # Table V, Pag 28 "Tref": 132.6312, "kref": 1e-3, "no": [1.308, 1.405, -1.036], "co": [-97, -1.1, -0.3], "Trefb": 132.6312, "rhorefb": 10.4477, "krefb": 1e-3, "nb": [8.743, 14.76, -16.62, 3.793, -6.142, -0.3778], "tb": [0.1, 0, 0.5, 2.7, 0.3, 1.3], "db": [1, 2, 3, 7, 7, 11], "cb": [0, 0, 2, 2, 2, 2], "critical": 3, "gnu": 0.63, "gamma": 1.2415, "R0": 1.01, "Xio": 0.11e-9, "gam0": 0.55e-1, "qd": 0.31e-9, "Tcref": 265.262 } _thermal = thermo0,
class CH4(MEoS): """Multiparameter equation of state for methane""" name = "methane" CASNumber = "74-82-8" formula = "CH4" synonym = "R-50" _refPropName = "METHANE" _coolPropName = "Methane" rhoc = unidades.Density(162.66) Tc = unidades.Temperature(190.564) Pc = unidades.Pressure(4599.2, "kPa") M = 16.0428 # g/mol Tt = unidades.Temperature(90.694) Tb = unidades.Temperature(111.667) f_acent = 0.01142 momentoDipolar = unidades.DipoleMoment(0.0, "Debye") id = 2 _Tr = unidades.Temperature(186.659809) _rhor = unidades.Density(163.413536) _w = 0.010528102 Fi1 = { "ao_log": [1, 3.00160], "pow": [0, 1], "ao_pow": [9.91243972, -6.33270087], "ao_exp": [0.008449, 4.6942, 3.4865, 1.6572, 1.4115], "titao": [648 / Tc, 1957 / Tc, 3895 / Tc, 5705 / Tc, 15080 / Tc], "ao_hyp": [], "hyp": [] } Fi2 = { "R": 8.314510, "ao_log": [1, 3.00088], "pow": [0, 1], "ao_pow": [19.597508817, -83.959667892], "ao_exp": [], "titao": [], "ao_hyp": [0.76315, 0.0046, 8.74432, -4.46921], "hyp": [4.306474465, 0.936220902, 5.577233895, 5.722644361] } Fi3 = { "ao_log": [1, 2.5998324], "pow": [0, -1. / 3, -2. / 3, -1], "ao_pow": [-10.413865, -3.3854083, 1.6900979, -0.3911541], "ao_exp": [4.7206715], "titao": [10.543907], "ao_hyp": [], "hyp": [] } CP4 = { "ao": 0.15438149595e2, "an": [ -0.18044750507e7, 0.77426666393e5, -0.13241658754e4, -0.51479005257e-1, 0.10809172196e-3, -0.65501783437e-7 ], "pow": [-3, -2, -1.001, 1, 2, 3], "ao_exp": [-0.67490056171e1], "exp": [3000], "ao_hyp": [], "hyp": [] } setzmann = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for methane of Setzmann and " "Wagner (1991)", "__doi__": { "autor": "Setzmann, U., Wagner, W.", "title": "A New Equation of State and Tables of " "Thermodynamic Properties for Methane Covering " "the Range from the Melting Line to 625 K at " "Pressures up to 1000 MPa", "ref": "J. Phys. Chem. Ref. Data, 20(6) (1991) 1061-1155", "doi": "10.1063/1.555898" }, "R": 8.31451, "cp": Fi1, "ref": "OTO", "Tmin": Tt, "Tmax": 625.0, "Pmax": 1000000.0, "rhomax": 40.072, "Pmin": 11.696, "rhomin": 28.142, "nr1": [ 0.43679010280e-1, 0.67092361990, -0.17655778590e01, 0.85823302410, -0.12065130520e01, 0.51204672200, -0.40000107910e-3, -0.12478424230e-1, 0.31002697010e-1, 0.17547485220e-2, -0.31719216050e-5, -0.22403468400e-5, 0.29470561560e-6 ], "d1": [1, 1, 1, 2, 2, 2, 2, 3, 4, 4, 8, 9, 10], "t1": [-0.5, 0.5, 1., 0.5, 1., 1.5, 4.5, 0., 1., 3., 1., 3., 3.], "nr2": [ 0.18304879090, 0.15118836790, -0.42893638770, 0.68940024460e-1, -0.14083139960e-1, -0.30630548300e-1, -0.29699067080e-1, -0.19320408310e-1, -0.11057399590, 0.99525489950e-1, 0.85484378250e-2, -0.61505556620e-1, -0.42917924230e-1, -0.18132072900e-1, 0.34459047600e-1, -0.23859194500e-2, -0.11590949390e-1, 0.66416936020e-1, -0.23715495900e-1, -0.39616249050e-1, -0.13872920440e-1, 0.33894895990e-1, -0.29273787530e-2 ], "d2": [1, 1, 1, 2, 4, 5, 6, 1, 2, 3, 4, 4, 3, 5, 5, 8, 2, 3, 4, 4, 4, 5, 6], "t2": [ 0., 1., 2., 0., 0., 2., 2., 5., 5., 5., 2., 4., 12., 8., 10., 10., 10., 14., 12., 18., 22., 18., 14. ], "c2": [1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4], "gamma2": [1] * 23, "nr3": [9.324799946e-5, -6.287171518, 12.71069467, -6.423953466], "d3": [2, 0, 0, 0], "t3": [2., 0., 1., 2.], "alfa3": [20, 40, 40, 40], "beta3": [200, 250, 250, 250], "gamma3": [1.07, 1.11, 1.11, 1.11], "epsilon3": [1] * 4 } MBWR = { "__type__": "MBWR", "__name__": "MBWR equation of state for ethane of Younglove and Ely " "(1987)", "__doi__": { "autor": "Younglove, B.A. and Ely, J.F.", "title": "Thermophysical Properties of Fluids. II. " "Methane, Ethane, Propane, Isobutane, and Normal " "Butane", "ref": "J. Phys. Chem. Ref. Data 16(4) (1987) 577-798", "doi": "10.1063/1.555785" }, "Tmin": 90.68, "Tmax": 600.0, "Pmax": 200000.0, "rhomax": 36.2029, "Pmin": 11.744, "rhomin": 28.147, "R": 8.31434, "cp": CP4, "ref": { "Tref": 298.15, "Pref": 101.325, "ho": 10018, "so": 186.266 }, "b": [ None, 0.9898937956e-4, 0.2199608275, -0.5322788000e1, 0.2021657962e3, -0.2234398926e5, 0.106794028e-3, 0.1457922469e-2, -9.265816666, 0.2915364732e4, 0.2313546209e-5, 0.1387214274e-2, 0.4780467451e-1, 0.1176103833e-3, -0.198209673e-2, -0.2512887756, 0.9748899826e-4, -0.1202192137e-5, 0.4128353939e-3, -0.7215842918e-5, 0.5081738255e4, -0.9198903192e6, -27.32264677, 0.7499024351e6, 0.01114060908, 0.1083955159e2, -0.4490960312e-3, -0.1380337847e2, -0.2371902232e-6, 0.3761652197e-3, -0.2375166954e-8, -0.1237640790e-6, 0.6766926453e-5 ] } GERG = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for methane 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": Fi2, "ref": "OTO", "Tmin": 90.6941, "Tmax": 625.0, "Pmax": 1000000.0, "rhomax": 40.072, "Pmin": 73.476, "rhomin": 29.249, "nr1": [ 0.57335704239162, -0.16760687523730e1, 0.23405291834916, -0.21947376343441, 0.16369201404128e-1, 0.15004406389280e-1 ], "d1": [1, 1, 2, 2, 4, 4], "t1": [0.125, 1.125, 0.375, 1.125, 0.625, 1.5], "nr2": [ 0.98990489492918e-1, 0.58382770929055, -0.7478686756039, 0.30033302857974, 0.20985543806568, -0.18590151133061e-1, -0.15782558339049, 0.12716735220791, -0.32019743894346e-1, -0.68049729364536e-1, 0.24291412853736e-1, 0.51440451639444e-2, -0.019084949733532, 0.55229677241291e-2, -0.44197392976085e-2, 0.040061416708429, -0.33752085907575e-1, -0.25127658213357e-2 ], "d2": [1, 1, 1, 2, 3, 6, 2, 3, 3, 4, 4, 2, 3, 4, 5, 6, 6, 7], "t2": [ 0.625, 2.625, 2.75, 2.125, 2, 1.75, 4.5, 4.75, 5, 4, 4.5, 7.5, 14, 11.5, 26, 28, 30, 16 ], "c2": [1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 3, 3, 3, 6, 6, 6, 6], "gamma2": [1] * 18 } friend = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for methane of Friend (1989)", "__doi__": { "autor": "Friend, D.G., Ely, J.F., Ingham, H.", "title": "Thermophysical Properties of Methane", "ref": "J. Phys. Chem. Ref. Data 18(2) (1989) 583-638", "doi": "10.1063/1.555828" }, "R": 8.31451, "cp": Fi3, "ref": { "Tref": 298.15, "Pref": 101.325, "ho": 10017.7, "so": 186.266 }, "Tt": 90.6854, "Tc": 190.551, "rhoc": 10.139, "M": 16.043, "Tmin": Tt, "Tmax": 620.0, "Pmax": 100000.0, "rhomax": 29.714, "Pmin": 11.694, "rhomin": 28.145, "nr1": [ 0.384436099659, -0.179692598800e1, 0.329444947369, 0.226312728442e-1, 0.759236768798e-1, 0.693758447259e-1, 0.241163263947e-1, 0.107009920854e-1, -0.380933275164e-1, 0.471537561143e-3, 0.556607678805e-3, 0.548759346533e-6, -0.999632699967e-4 ], "d1": [1, 1, 1, 2, 2, 2, 3, 3, 3, 6, 7, 7, 8], "t1": [0, 1.5, 2.5, -0.5, 1.5, 2, 0, 1, 2.5, 0, 2, 5, 2], "nr2": [ -0.128087979280, 0.380198873377e-1, 0.139226650551, -0.874996348859e-1, -0.334894165760e-2, -0.517576297122e-1, 0.252835179116e-1, 0.518703205950e-3, -0.166770594525e-2, -0.607401927389e-3, -0.972915359991e-4, -0.298844010462e-4, -0.130940111124e-1, 0.198175833798e-1, 0.208465762327e-1, -0.358025052631e-1, -0.203486851741, 0.215964755088, -0.429340628249e-2 ], "d2": [1, 1, 2, 2, 3, 3, 5, 6, 7, 8, 10, 2, 3, 3, 4, 4, 5, 5, 5], "t2": [ 5, 6, 3.5, 5.5, 3, 7, 6, 8.5, 4, 6.5, 5.5, 22, 11, 18, 11, 23, 17, 18, 23 ], "c2": [2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 4, 4, 4, 4, 4, 4, 4, 4], "gamma2": [1] * 19 } shortSpan = { "__type__": "Helmholtz", "__name__": "short Helmholtz equation of state for methane 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": Fi1, "ref": "OTO", "M": 16.043, "Tmin": Tt, "Tmax": 750.0, "Pmax": 100000.0, "rhomax": 40.072, "Pmin": 11.661, "rhomin": 28.167, "nr1": [ 0.89269676, -0.25438282e1, 0.64980978, 0.20793471e-1, 0.70189104e-1, 0.23700378e-3 ], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875], "nr2": [ 0.16653334, -0.43855669e-1, -0.1572678, -0.35311675e-1, -0.29570024e-1, 0.14019842e-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 } sun = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for methane 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.31451, "cp": Fi1, "ref": "OTO", "Tmin": Tt, "Tmax": 625.0, "Pmax": 1000000.0, "rhomax": 40.072, "Pmin": 11.696, "rhomin": 28.142, "nr1": [ 1.25595787, 8.48007435e-1, -3.00939285, 5.99544996e-2, 2.57003062e-4, -2.85914246e-2 ], "d1": [1, 1, 1, 3, 7, 2], "t1": [1.5, 0.25, 1.25, 0.25, 0.875, 1.375], "nr2": [ -6.83210861e-2, -3.47523515e-2, 1.04637008e-1, -1.09884198e-2, -0.125124331, -5.53450960e-3, -1.51182884e-2, -2.04800000e-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 = setzmann, MBWR, GERG, friend, shortSpan, sun eq = setzmann, GERG, friend, shortSpan, sun _surface = { "sigma": [0.03825, -0.006024, -0.0007065], "exp": [1.191, 5.422, 0.6161] } _dielectric = { "eq": 3, "Tref": 273.16, "rhoref": 1000., "a0": [], "expt0": [], "expd0": [], "a1": [6.5443, 0.0133], "expt1": [0, 1], "expd1": [1, 1], "a2": [8.4578, 3.7196, -352.97, -100.65], "expt2": [0, 1, 0, 1], "expd2": [2, 2, 3, 3] } _melting = { "eq": 1, "Tref": Tt, "Pref": 11.696, "Tmin": Tt, "Tmax": 625.0, "a1": [1, 0.247568e5, -0.736602e4, -0.247568e5, 0.736602e4], "exp1": [0, 1.85, 2.1, 0, 0], "a2": [], "exp2": [], "a3": [], "exp3": [] } _sublimation = { "eq": 3, "Tref": Tt, "Pref": 11.696, "Tmin": Tt, "Tmax": Tt, "a1": [], "exp1": [], "a2": [-12.84], "exp2": [1], "a3": [], "exp3": [] } _vapor_Pressure = { "eq": 3, "n": [-6.036219, 1.409353, -0.4945199, -1.443048], "t": [1, 1.5, 2, 4.5] } _liquid_Density = { "eq": 2, "n": [1.9906389, -0.78756197, 0.036976723], "t": [0.354, 0.5, 2.5] } _vapor_Density = { "eq": 2, "n": [-1.880284, -2.8526531, -3.000648, -5.251169, -13.191859, -37.553961], "t": [0.354, 5 / 6, 3 / 2, 5 / 2, 25 / 6, 47 / 6] } 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": 190.564, "no": [2.60536, -18.5247, 23.4216], "to": [0, 0.25, 0.5], "a": [-3.12118e-5, 1.99422e-7, 0], "b": [5.98858e-5, -4.91143e-5, 0], "c": [3.49668e-5, -1.73176e-5, 0], "A": [-8.52992e-10, -3.58009e-10, 0], "B": [1.60099e-8, 8.50221e-10, 0], "C": [-3.55631e-7, 2.80326e-7, 0] } visco1 = { "__name__": "Vogel (2000)", "__doi__": { "autor": "Vogel, E., Wilhelm, J., Küchenmeister, C., " "Jaesche, M.", "title": "High-precision viscosity measurements on methane", "ref": "High Temperatures-High Pressures 32(1) (2000) 73-81", "doi": "10.1068/htwu359" }, "eq": 1, "omega": 1, "ek": 160.78, "sigma": 0.37333, "n_chapman": 0.021357, "collision": [0.215309028, -0.46256942, 0.051313823, 0.030320660, -0.0070047029], "Tref_virial": 159.7, "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": 190.564, "rhoref_res": 10.139 * M, "nr": [ -3.02256904347, 17.6965130175, 3.11150846518, -21.5685107769, 0.672852409238, 10.2387524315, -1.09330775541, -1.20030749419 ], "tr": [0, 1, 0, 1, 0, 1, 0, 1], "dr": [2, 2, 3, 3, 4, 4, 5, 5], "gr": [0, 0, 0, 0, 0, 0, 0, 0], "cr": [0, 0, 0, 0, 0, 0, 0, 0], "CPf": -21.1009923406, "CPg1": 5.10860501398, "CPgi": [], "CPti": [] } visco2 = { "eq": 2, "omega": 2, "__name__": "Younglove (1987)", "__doi__": { "autor": "Younglove, B.A., Ely, J.F.", "title": "Thermophysical Properties of Fluids. II. Methane, " "Ethane, Propane, Isobutane, and Normal Butane", "ref": "J. Phys. Chem. Ref. Data 16(4) (1987) 577-798", "doi": "10.1063/1.555785" }, "ek": 168., "sigma": 0.368, # "n_chapman": 0.1069188/M**0.5, "F": [0.16969859271, -0.13337234608e-1, 1.4, 168], "E": [ -0.1620427429e2, 0.4270589027e3, 0.1402596278e2, -0.3916837745e4, -0.347709909e-1, 0.2136542674e2, 0.1436802482e4 ], "rhoc": 10.15 } visco3 = { "__name__": "Friend (1989)", "__doi__": { "autor": "Friend, D.G., Ely, J.F., Ingham, H.", "title": "Thermophysical Properties of Methane", "ref": "J. Phys. Chem. Ref. Data 18(2) (1989) 583-638", "doi": "10.1063/1.555828" }, "eq": 1, "omega": 2, "Tref": 174., "muref": 10.0, "ek": 174., "sigma": 0.36652, "n_chapman": 10.50 / M**0.5 * 0.36652**2, "Tref_res": 190.551, "rhoref_res": 10.139 * M, "muref_res": 12.149, "nr_num": [ 0.41250137, -0.14390912, 0.10366993, 0.40287464, -0.24903524, -0.12953131, 0.06575776, 0.02566628, -0.03716526 ], "tr_num": [0, 1, 0, 1, 1.5, 0, 2, 0, 1], "dr_num": [1, 1, 2, 2, 2, 3, 3, 4, 4], "gr_num": [0, 0, 0, 0, 0, 0, 0, 0, 0], "cr_num": [0, 0, 0, 0, 0, 0, 0, 0, 0], "nr_den": [1.0, -0.38798341, 0.03533815], "tr_den": [0, 0, 1], "dr_den": [0, 1, 1], "gr_den": [0, 0, 0], "cr_den": [0, 0, 0] } _viscosity = visco0, visco1, visco2, visco3 thermo0 = { "eq": 1, "critical": "CH4", "__name__": "Friend (1989)", "__doi__": { "autor": "Friend, D.G., Ely, J.F., and Ingham, H.", "title": "Thermophysical Properties of Methane", "ref": "J. Phys. Chem. Ref. Data 18, 583 (1989)", "doi": "10.1063/1.555828" }, "Tref": 174., "kref": 1e-3, "no": [1.45885, -0.4377162, 0], "co": [0, -1, -96], "Trefb": 190.551, "rhorefb": 10.139, "krefb": 6.29638e-3, "nb": [ 1.5554612, 1., 2.4149207, 0.55166331, -0.52837734, 0.073809553, 0.24465507, -0.047613626 ], "tb": [0, 0, 0, 0, 0, -1, 0, -1], "db": [2, 0, 1, 3, 4, 4, 5, 5], "cb": [0, -99, 0, 0, 0, 0, 0, 0] } thermo1 = { "eq": 2, "omega": 2, "__name__": "Younglove (1987)", "__doi__": { "autor": "Younglove, B.A. and Ely, J.F.", "title": "Thermophysical Properties of Fluids. II. Methane, Ethane, Propane, Isobutane, and Normal Butane ", "ref": "J. Phys. Chem. Ref. Data 16, 577 (1987)", "doi": "10.1063/1.555785" }, "visco": visco2, "n_chapman": 0.1069188, "G": [0.1346953698e1, -0.3254677753], "E": [ 0.2325800819e-2, -0.2477927999, 0.3880593713e2, -0.1579519146e-6, 0.3717991328e-2, -0.9616989434, -0.3017352774e-1, 0.4298153386 ], "critical": 2, "X": [37.42368, 3.16714, 0.78035, 0.60103], "Z": 6.512707e-10 } _thermal = thermo0, thermo1
class mXylene(MEoS): """Multiparameter equation of state for m-xylene""" name = "m-xylene" CASNumber = "108-38-3" formula = "C8H10" synonym = "1,3-dimethylbenzene" _refPropName = "MXYLENE" _coolPropName = "m-Xylene" rhoc = unidades.Density(282.929725) Tc = unidades.Temperature(616.89) Pc = unidades.Pressure(3534.6, "kPa") M = 106.165 # g/mol Tt = unidades.Temperature(225.3) Tb = unidades.Temperature(412.214) f_acent = 0.326 momentoDipolar = unidades.DipoleMoment(0.3, "Debye") id = 43 Fi1 = { "ao_log": [1, 1.169909], "pow": [0, 1], "ao_pow": [12.652887, -0.45975624], "ao_exp": [4.44312, 2.862794, 24.83298, 16.26077], "titao": [160 / Tc, 190 / Tc, 1333 / Tc, 3496 / Tc] } zhou = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for ethylbenzene of Zhou et " "al. (2012).", "__doi__": { "autor": "Zhou, Y., Lemmon, E.W., Wu, J.", "title": "Thermodynamic Properties of o-Xylene, m-Xylene, " "p-Xylene, and Ethylbenzene", "ref": "J. Phys. Chem. Ref. Data 41, 023103 (2012).", "doi": "10.1063/1.3703506" }, "R": 8.314472, "cp": Fi1, "ref": "OTO", "Tmin": Tt, "Tmax": 700.0, "Pmax": 200000.0, "rhomax": 8.677, "nr1": [ 0.000012791017, 0.041063111, 1.505996, -2.3095875, -0.46969, 0.171031 ], "d1": [8, 4, 1, 1, 2, 3], "t1": [1.0, 0.91, 0.231, 0.772, 1.205, 0.323], "nr2": [-1.001728, -0.3945766, 0.6970578, -0.3002876, -0.024311], "d2": [1, 3, 2, 2, 7], "t2": [2.7, 3.11, 0.768, 4.1, 0.818], "c2": [2, 2, 1, 2, 1], "gamma2": [1] * 5, "nr3": [0.815488, -0.330647, -0.123393, -0.54661], "d3": [1, 1, 3, 3], "t3": [2.0, 2.9, 3.83, 0.5], "alfa3": [1.0244, 1.3788, 0.9806, 6.3563], "beta3": [1.66, 1.9354, 1.0323, 78], "gamma3": [1.1013, 0.6515, 0.4975, 1.26], "epsilon3": [0.713, 0.9169, 0.6897, 0.7245] } eq = zhou, _PR = [0.0465, -21.1362] _surface = {"sigma": [0.0661], "exp": [1.29]} _vapor_Pressure = { "eq": 3, "n": [-7.5635, 1.2857, -3.2346, -1.9018], "t": [1.0, 1.5, 3.1, 5.6] } _liquid_Density = { "eq": 1, "n": [0.43346, 3.8716, -3.0144, 1.619], "t": [0.16, 0.6, 1.0, 1.5] } _vapor_Density = { "eq": 2, "n": [-1.1597, -6.0358, -16.712, -45.482, -98.418], "t": [0.26, 0.78, 2.6, 5.7, 11.7] } visco0 = { "__name__": "Cao (2016)", "__doi__": { "autor": "Cao, F.L., Meng, X.Y., Wu, J.T., Vesovic, V.", "title": "Reference Correlation of the Viscosity of " "meta-Xylene from 273 to 673 K and up to 110 MPa", "ref": "J. Phys. Chem. Ref. Data 45(1) (2016) 013103", "doi": "10.1063/1.4941241" }, "eq": 1, "omega": 3, "collision": [-1.4933, 473.2, -57033], "sigma": 1, "n_chapman": 0.22115 / M**0.5, "Tref_res": 616.89, "rhoref_res": 2.665 * M, "nr": [ -0.26895, 0.320971, -0.0290018, 1.72866e-10, 14.7728, -18.9852, 17.1128 ], "dr": [112 / 15, 112 / 15, 2 / 3 + 3.3, 68 / 3, 34 / 15, 19 / 15, 16 / 15], "tr": [-0.5, -0.2, -0.5, 2.7, -0.5, -1.5, -0.5], "special": "_vir" } def _vir(self, rho, T, fase): # The initial density dependence has a different expresion, without muo # and other normal method calculation so hardcoded here muB = 0 if rho: for i, n in enumerate([13.2814, -10862.4, 1664060]): muB += n / T**i return muB * rho / self.M _viscosity = visco0, thermo0 = { "__name__": "Mylona (2014)", "__doi__": { "autor": "Mylona, S.K., Antoniadis, K.D., Assael, M.J., " "Huber, M.L., Perkins, R.A.", "title": "Reference Correlations of the Thermal " "Conductivity of o-Xylene, m-Xylene, p-Xylene, " "and Moderate Pressures", "ref": "J. Phys. Chem. Ref. Data 43(4) (2014) 043104", "doi": "10.1063/1.4901166" }, "eq": 1, "Toref": 616.89, "koref": 1e-3, "no_num": [0.242107, 13.522, -123.168, 296.882, -107.973, 18.686, -1.29167], "to_num": [0, 1, 2, 3, 4, 5, 6], "no_den": [-0.850118, 3.11646, 0.0001], "to_den": [0, 1, 2], "Tref_res": 616.89, "rhoref_res": 282.9297, "kref_res": 1e-3, "nr": [ -6.79314e1, 2.25778e2, -1.85693e2, 6.19006e1, -7.11664, 5.92537e1, -1.62626e2, 1.33036e2, -4.49051e1, 5.6186 ], "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.235e-9, "gam0": 0.057, "qd": 0.713e-9, "Tcref": 925.3 } _thermal = thermo0,
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, _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 HCl(MEoS): """Multiparameter equation of state for hydrogen chloride""" name = "hydrogen chloride" CASNumber = "7647-01-0 " formula = "HCl" synonym = "" rhoc = unidades.Density(410.97) Tc = unidades.Temperature(324.55) Pc = unidades.Pressure(8263.00, "kPa") M = 36.460939 # g/mol Tt = unidades.Temperature(131.4) Tb = unidades.Temperature(188.199) f_acent = 0.12875 momentoDipolar = unidades.DipoleMoment(1.079, "Debye") id = 104 Fi1 = { "ao_log": [1, 2.5], "pow": [0, 1, -1, -2], "ao_pow": [7.913048, -3.217575, -4.149937e-3, 8.019202e-4], "ao_exp": [1.054392], "titao": [1.241138e1] } helmholtz1 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for hydrogen chloride 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" }, "__test__": """ >>> st=HCl(T=170, rho=0.01) >>> print "%0.0f %0.1f %0.9f %0.9f %0.9f %0.9f %0.9f %0.9f" % ( \ st.T, st.rho, st.P.MPa, st.h.kJkg, st.s.kJkgK, st.cv.kJkgK, st.cp.kJkgK, st.w) 170 0.01 0.000387586 -102.414301628 0.820319743 0.571157223 0.799483403 232.898528911 >>> st=HCl(T=170, rho=1230) >>> print "%0.0f %0.0f %0.9f %0.9f %0.9f %0.9f %0.9f %0.8f" % ( \ st.T, st.rho, st.P.MPa, st.h.kJkg, st.s.kJkgK, st.cv.kJkgK, st.cp.kJkgK, st.w) 170 1230 1.229128198 -561.339398154 -2.888318217 1.149719990 1.553944836 999.438819195 >>> st=HCl(T=280, rho=0.1) >>> print "%0.0f %0.1f %0.9f %0.9f %0.9f %0.9f %0.9f %0.9f" % ( \ st.T, st.rho, st.P.MPa, st.h.kJkg, st.s.kJkgK, st.cv.kJkgK, st.cp.kJkgK, st.w) 280 0.1 0.006381824 -14.914381738 0.580027609 0.571466689 0.800079245 298.836360637 >>> st=HCl(T=280, rho=900) >>> print "%0.0f %0.0f %0.9f %0.9f %0.9f %0.9f %0.9f %0.8f" % ( \ st.T, st.rho, st.P.MPa, st.h.kJkg, st.s.kJkgK, st.cv.kJkgK, st.cp.kJkgK, st.w) 280 900 3.421902004 -371.035989751 -2.044095830 0.961502216 2.150825813 577.782761523 """, # Table 9, Pag 26 "R": 8.314472, "cp": Fi1, "ref": "OTO", "Tmin": 155.0, "Tmax": 330.0, "Pmax": 20000.0, "rhomax": 33.8145, "Pmin": 0.7, "rhomin": 33.8145, "nr1": [ -.40937325, 0.943994574, -.178830477e1, 0.128619044, 0.439018427e-2, 0.130480908e-1, 0.169387782e-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": [ 0.751559060, -.800007427, 0.430935939, 0.454319457e-2, -.152172259, -.436174059e-1, -.970625964e-2, 0.101144098e-1, 0.376991644e-2 ], "d2": [1, 1, 2, 5, 1, 3, 4, 5, 2], "t2": [2.375, 3.0, 2.625, 1.875, 4.5, 5.78, 5.375, 2.75, 14.5], "c2": [1, 1, 1, 1, 2, 2, 2, 2, 3], "gamma2": [1] * 9 } eq = helmholtz1, _vapor_Pressure = { "eq": 6, "ao": [-0.01065138, -6.15979914, 1.55860976, -8.42734117], "exp": [1.0, 2.0, 6.0, 11.0] } _liquid_Density = { "eq": 2, "ao": [1.89232034, 0.83621066, -0.22094602, 4.70971253, -5.34396174], "exp": [1.0, 2.0, 4.0, 11.0, 13.0] } _vapor_Density = { "eq": 4, "ao": [-2.95523223, -8.10448179, -14.78392979, -87.13352586], "exp": [1.29, 4.2, 11.1, 24.0] }
class C1Oleate(MEoS): """Multiparameter equation of state for methyl oleate""" name = "methyl oleate" CASNumber = "112-62-9" formula = "C19H36O2" synonym = "" _refPropName = "MOLEATE" _coolPropName = "MethylOleate" rhoc = unidades.Density(241.000222029) Tc = unidades.Temperature(782.0) Pc = unidades.Pressure(1246.0, "kPa") M = 296.48794 # g/mol Tt = unidades.Temperature(253.47) Tb = unidades.Temperature(627.18) f_acent = 0.91 momentoDipolar = unidades.DipoleMoment(1.63, "Debye") # id = 919 CP1 = { "ao": 0.0, "an": [90.2385], "pow": [0.146118], "ao_exp": [234.797, 335.768, 431.66], "exp": [613.529, 1405.31, 2867.76], "ao_hyp": [], "hyp": [] } huber = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for methyl oleate of Huber " "et al. (2009).", "__doi__": { "autor": "Huber, M.L., Lemmon, E.W., Kazakov, A., Ott, " "L.S., and Bruno, T.J.", "title": "Model for the Thermodynamic Properties of a " "Biodiesel Fuel", "ref": "Energy Fuels, 23 (7) (2009) 3790–3797", "doi": "10.1021/ef900159g" }, "R": 8.314472, "cp": CP1, "ref": "NBP", "Tmin": Tt, "Tmax": 1000.0, "Pmax": 50000.0, "rhomax": 3.05, "Pmin": 0.0000000004, "rhomin": 3.05, "nr1": [0.4596121e-1, 2.2954, -3.554366, -0.2291674, 0.6854534e-1], "d1": [4, 1, 1, 2, 3], "t1": [1, 0.34, 1.14, 1.4, 0.6], "nr2": [-1.535778, -0.7334697, 1.712700, -1.471394, -0.1724678e-1], "d2": [1, 3, 2, 2, 7], "t2": [3.3, 4.1, 1.9, 3.8, 1.3], "c2": [2, 2, 1, 2, 1], "gamma2": [1] * 5, "nr3": [0.2115470e1, -0.7555374, -0.4134269], "d3": [1, 1, 3], "t3": [3.4, 3.8, 4.0], "alfa3": [1.1, 1.6, 1.1], "beta3": [0.9, 0.65, 0.75], "gamma3": [1.14, 0.65, 0.77], "epsilon3": [0.79, 0.9, 0.76], "nr4": [] } eq = huber, _vapor_Pressure = { "eq": 3, "n": [-0.13900e2, 0.16246e2, -0.15568e2, -0.73568e1, -0.48739e1], "t": [1.0, 1.5, 1.93, 4.2, 8.0] } _liquid_Density = { "eq": 1, "n": [-0.19920e2, 0.12230e3, -0.23582e3, 0.21009e3, -0.73435e2], "t": [0.461, 0.6, 0.75, 0.91, 1.05] } _vapor_Density = { "eq": 2, "n": [-13.426, 1.8069e2, -1.1108e3, 1.3265e3, -4.6421e2, -2.1070e2], "t": [0.667, 1.71, 2.2, 2.46, 3.0, 9.7] } thermo0 = { "__name__": "Perkins (2010)", "__doi__": { "autor": "Perkins, R.A., Huber, M.L.", "title": "Measurement and Correlation of the Thermal " "Conductivities of Biodiesel Constituent Fluids: " "Methyl Oleate and Methyl Linoleate", "ref": "Energy Fuels 25(5) (2011) 2383-2388", "doi": "10.1021/ef200417x" }, "eq": 1, "Toref": 782.0, "koref": 1, "no": [-2.7125e-4, 2.59365e-3, 0.0350241, -9.02273e-3], "to": [0, 1, 2, 3], "Tref_res": 782.0, "rhoref_res": 241, "kref_res": 1., "nr": [-0.0410106, 0.0606657, 0.0328443, -0.0498407, -0.00418506, 0.0121752], "tr": [0, 1, 0, 1, 0, 1], "dr": [1, 1, 2, 2, 3, 3], "critical": 3, "gnu": 0.63, "gamma": 1.239, "R0": 1.03, "Xio": 0.194e-9, "gam0": 0.0496, "qd": 8.75e-10, "Tcref": 1173 } _thermal = thermo0,
def _calculo(self): self._initialization() x = self._x() m = refprop.wmol(x)["wmix"] self.M = unidades.Dimensionless(m) crit = refprop.critp(x) self.Pc = unidades.Pressure(crit["pcrit"], "kPa") self.Tc = unidades.Temperature(crit["tcrit"]) self.rhoc = unidades.Density(crit["Dcrit"]*self.M) args = self.args() flash = refprop.flsh(*args) # check if ['q'] in fld if 'q' in flash.keys(): x = flash['q'] elif 'h' in flash.keys(): x = refprop.flsh('ph', flash['p'], flash['h'], flash['x'])['q'] elif 's' in flash.keys(): x = refprop.flsh('ps', flash['p'], flash['s'], flash['x'])['q'] if 0 < x < 1: region = 4 else: region = 1 if x < 0: x = 0 elif x > 1: x = 1 self.x = unidades.Dimensionless(x) self.T = unidades.Temperature(flash["t"]) self.P = unidades.Pressure(flash["p"], "kPa") self.Tr = unidades.Dimensionless(self.T/self.Tc) self.Pr = unidades.Dimensionless(self.P/self.Pc) self.rho = unidades.Density(flash["D"]*self.M) self.v = unidades.SpecificVolume(1./self.rho) self.phase = self.getphase(Tc=self.Tc, Pc=self.Pc, T=self.T, P=self.Pc, x=self.x, region=region) if flash["nc"] == 1: name = refprop.name(flash["nc"]) self.name = name["hname"] self.synonim = name["hn80"] self.CAS = name["hcas"] info = refprop.info(flash["nc"]) self.R = unidades.SpecificHeat(info["Rgas"]/self.M) self.Tt = unidades.Temperature(info["ttrp"]) self.Tb = unidades.Temperature(info["tnbpt"]) self.f_accent = unidades.Dimensionless(info["acf"]) self.momentoDipolar = unidades.DipoleMoment(info["dip"], "Debye") self._doc = {} for htype in ['EOS', 'CP0', 'ETA', 'VSK', 'TCX', 'TKK', 'STN', 'DE ', 'MLT', 'SBL', 'PS ', 'DL ', 'DV ']: self._doc[htype] = refprop.getmod(flash["nc"], htype)["hcite"] else: self.name = "" self.synonim = "" self.CAS = "" rmix = refprop.rmix2(flash["x"]) self.R = unidades.SpecificHeat(rmix["Rgas"]/self.M) self.Tt = unidades.Temperature(None) self.Tb = unidades.Temperature(None) self.f_accent = unidades.Dimensionless(None) self.momentoDipolar = unidades.DipoleMoment(None) self._doc = {} self._cp0(flash) self.Liquido = ThermoRefProp() self.Gas = ThermoRefProp() if self.x == 0.: # liquid phase self.fill(self.Liquido, flash["t"], flash["D"], flash["x"]) self.fill(self, flash["t"], flash["D"], flash["x"]) self.fillNone(self.Gas) elif self.x == 1.: # vapor phase self.fill(self.Gas, flash["t"], flash["D"], flash["x"]) self.fill(self, flash["t"], flash["D"], flash["x"]) self.fillNone(self.Liquido) else: # Two phase self.fillNone(self) self.fill(self.Liquido, flash["t"], flash["Dliq"], flash["xliq"]) self.fill(self.Gas, flash["t"], flash["Dvap"], flash["xvap"]) self.v = unidades.SpecificVolume(x*self.Gas.v+(1-x)*self.Liquido.v) self.rho = unidades.Density(1./self.v) self.u = unidades.Enthalpy(flash["e"]/self.M, "Jg") self.h = unidades.Enthalpy(flash["h"]/self.M, "Jg") self.s = unidades.SpecificHeat(flash["s"]/self.M, "JgK") self.a = unidades.Enthalpy(self.u-self.T*self.s) self.g = unidades.Enthalpy(self.h-self.T*self.s) self.rhoM = unidades.MolarDensity(self.rho/self.M) self.hM = unidades.MolarEnthalpy(self.h*self.M) self.sM = unidades.MolarSpecificHeat(self.s*self.M) self.uM = unidades.MolarEnthalpy(self.u*self.M) self.aM = unidades.MolarEnthalpy(self.a*self.M) self.gM = unidades.MolarEnthalpy(self.g*self.M) if self.x < 1 and self.T <= self.Tc: surten = refprop.surten(flash["t"], flash["Dliq"], flash["Dvap"], flash["xliq"], flash["xvap"]) self.sigma = unidades.Tension(surten["sigma"]) else: self.sigma = unidades.Tension(None) if 0 < self.x < 1: self.Hvap = unidades.Enthalpy(self.Gas.h-self.Liquido.h) self.Svap = unidades.SpecificHeat(self.Gas.s-self.Liquido.s) self.K = [] for x, y in zip(self.Liquido.fraccion, self.Gas.fraccion): self.K.append(unidades.Dimensionless(y/x)) else: self.Hvap = unidades.Enthalpy(None) self.Svap = unidades.SpecificHeat(None) self.K = [unidades.Dimensionless(1)]*flash["nc"] self.invT = unidades.InvTemperature(-1/self.T) # NOT supported on Windows if sys.platform != "win32": excess = refprop.excess(flash["t"], flash["D"], flash["x"]) self.vE = unidades.Volume(excess["vE"]/self.M) self.uE = unidades.Enthalpy(excess["eE"]/self.M, "Jg") self.hE = unidades.Enthalpy(excess["hE"]/self.M, "Jg") self.sE = unidades.SpecificHeat(excess["sE"]/self.M, "JgK") self.aE = unidades.Enthalpy(excess["aE"]/self.M, "Jg") self.gE = unidades.Enthalpy(excess["gE"]/self.M, "Jg") else: self.vE = unidades.Volume(0) self.uE = unidades.Enthalpy(0) self.hE = unidades.Enthalpy(0) self.sE = unidades.SpecificHeat(0) self.aE = unidades.Enthalpy(0) self.gE = unidades.Enthalpy(0) self.csat = [] self.dpdt_sat = [] self.cv2p = [] if self.Tt <= flash["t"] <= self.Tc: for i in range(1, flash["nc"]+1): dat = refprop.dptsatk(i, flash["t"], kph=2) cs = unidades.SpecificHeat(dat["csat"]/self.M, "JgK") self.csat.append(cs) self.dpdt_sat.append( unidades.PressureTemperature(dat["dpdt"], "kPaK")) cv2 = refprop.cv2pk(i, flash["t"], flash["D"]) cv = unidades.SpecificHeat(cv2["cv2p"]/self.M, "JgK") self.cv2p.append(cv)
class H2(MEoS): """Multiparameter equation of state for hydrogen (normal)""" name = "hydrogen" CASNumber = "1333-74-0" formula = "H2" synonym = "R-702" rhoc = unidades.Density(31.26226704) Tc = unidades.Temperature(33.145) Pc = unidades.Pressure(1296.4, "kPa") M = 2.01588 # g/mol Tt = unidades.Temperature(13.957) Tb = unidades.Temperature(20.369) f_acent = -0.219 momentoDipolar = unidades.DipoleMoment(0.0, "Debye") id = 1 Fi1 = { "ao_log": [1, 1.5], "pow": [0, 1], "ao_pow": [-1.4579856475, 1.888076782], "ao_exp": [1.616, -0.4117, -0.792, 0.758, 1.217], "titao": [ 16.0205159149, 22.6580178006, 60.0090511389, 74.9434303817, 206.9392065168 ] } Fi2 = { "ao_log": [1, 1.47906], "pow": [0, 1], "ao_pow": [13.796443393, -175.864487294], "ao_exp": [], "titao": [], "ao_hyp": [0.95806, 0.45444, 1.56039, -1.3756], "hyp": [6.891654113, 9.84763483, 49.76529075, 50.367279301] } CP1 = { "ao": 0.72480209e3, "an": [ 0.12155215e11, -0.36396763e10, 0.43375265e9, -0.23085817e8, -0.38680927e4, 0.88240136e5, -0.78587085e4, -0.18426806e3, 0.21801550e2, -0.13051820e1, 0.21003175e-1, 0.23911604e-2, -0.18240547e-3, 0.56149561e-5, -0.73803310e-7, 0.66357755e-11 ], "pow": [-7, -6, -5, -4, -3, -2, -1.001, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 5], "ao_exp": [], "exp": [], "ao_hyp": [], "hyp": [] } helmholtz1 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for normal hydrogen of Leachman et al. (2009).", "__doi__": { "autor": "Leachman, J.W., Jacobsen, R.T, Penoncello, S.G., Lemmon, E.W.", "title": "Fundamental equations of state for parahydrogen, normal hydrogen, and orthohydrogen", "ref": "J. Phys. Chem. Ref. Data, 38 (2009), 721 – 748", "doi": "10.1063/1.3160306" }, "__test__": """ >>> st=H2(T=13.957, x=0.5) >>> print "%0.6g %0.5g %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.kPa, 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) 13.957 7.3580 77.004 0.12985 −53.926 399.83 −3.0723 29.438 5.1616 6.2433 7.0212 10.564 1269.2 307.14 """, # Table 14, Pag 746 "R": 8.314472, "cp": Fi1, "ref": "NBP", "Tmin": Tt, "Tmax": 1000.0, "Pmax": 2000000.0, "rhomax": 102.0, "Pmin": 7.36, "rhomin": 38.2, "nr1": [-6.93643, 0.01, 2.1101, 4.52059, 0.732564, -1.34086, 0.130985], "d1": [1, 4, 1, 1, 2, 2, 3], "t1": [0.6844, 1., 0.989, 0.489, 0.803, 1.1444, 1.409], "nr2": [-0.777414, 0.351944], "d2": [1, 3], "t2": [1.754, 1.311], "c2": [1, 1], "gamma2": [1] * 2, "nr3": [-0.0211716, 0.0226312, 0.032187, -0.0231752, 0.0557346], "d3": [2, 1, 3, 1, 1], "t3": [4.187, 5.646, 0.791, 7.249, 2.986], "alfa3": [1.685, 0.489, 0.103, 2.506, 1.607], "beta3": [0.171, 0.2245, 0.1304, 0.2785, 0.3967], "gamma3": [0.7164, 1.3444, 1.4517, 0.7204, 1.5445], "epsilon3": [1.506, 0.156, 1.736, 0.67, 1.662], "nr4": [] } MBWR = { "__type__": "MBWR", "__name__": "MBWR equation of state for hydrogen of Younglove (1982)", "__doi__": { "autor": "Younglove, B.A.", "title": "Thermophysical Properties of Fluids. I. Argon, Ethylene, Parahydrogen, Nitrogen, Nitrogen Trifluoride, and Oxygen", "ref": "J. Phys. Chem. Ref. Data, Vol. 11, Suppl. 1, pp. 1-11, 1982.", "doi": "" }, "R": 8.31434, "cp": CP1, "ref": "IIR", "Tmin": Tt, "Tmax": 400.0, "Pmax": 121000.0, "rhomax": 38.148, "Pmin": 7.70, "rhomin": 38.3, "b": [ None, 0.4675528393416e-3, 0.4289274251454e-1, -0.5164085596504, 0.2961790279801e1, -0.3027194968412e2, 0.1908100320379e-4, -0.1339776859288e-2, 0.3056473115421, 0.5161197159532e2, 0.1999981550224e-6, 0.2896367059356e-3, -0.2257803939041e-1, -0.2287392761826e-5, 0.2446261478645e-4, -0.1718181601119e-2, -0.5465142603459e-6, 0.4051941401315e-8, 0.1157595123961e-5, -0.1269162728389e-7, -0.4983023605519e2, -0.1606676092098e3, -0.1926799185310, 0.9319894638928e1, -0.3222596554434e-3, 0.1206839307669e-2, -0.3841588197470e-6, -0.4036157453608e-4, -0.1250868123513e-9, 0.1976107321888e-8, -0.2411883474011e-12, -0.4127551498251e-12, 0.8917972883610e-11 ] } GERG = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for hydrogen 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, 2012, 57 (11), pp 3032-3091", "doi": "10.1021/je300655b" }, "R": 8.314472, "cp": Fi2, "ref": "OTO", "Tmin": Tt, "Tmax": 400.0, "Pmax": 121000.0, "rhomax": 38.148, # "Pmin": 0.61166, "rhomin": 55.497, "nr1": [ 0.53579928451252e1, -0.62050252530595e1, 0.13830241327086, -0.71397954896129e-1, 0.15474053959733e-1 ], "d1": [1, 1, 2, 2, 4], "t1": [0.5, 0.625, 0.384, 0.625, 1.125], "nr2": [ -0.14976806405771, -0.26368723988451e-1, 0.56681303156066e-1, -0.60063958030436e-1, -0.45043942027132, 0.42478840244500, -0.21997640827139e-1, -0.1049952137453e-1, -0.28955902866816e-2 ], "d2": [1, 5, 5, 5, 1, 1, 2, 5, 1], "t2": [2.625, 0, 0.25, 1.375, 4, 4.25, 5, 8, 8], "c2": [1, 1, 1, 1, 2, 2, 3, 3, 5], "gamma2": [1] * 9, "nr3": [], "nr4": [] } helmholtz3 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for hydrogen of Bender (1982).", "__doi__": { "autor": "Bender, E.", "title": "Equation of state of normal hydrogen in the range 18 to 700 K and 1 to 500 bar", "ref": "VDI-Forschungsheft, no. 609, 1982, p. 15-20", "doi": "" }, "R": 8.3143, "cp": CP1, "ref": "IIR", "Tmin": 18.0, "Tmax": 700.0, "Pmax": 50000.0, "rhomax": 38.74, "Pmin": 8.736, "rhomin": 38.7, "nr1": [ 0.133442326203e1, -0.104116843433e1, 0.227202245707, 0.300374270906, -0.463984214813, -0.178010492282e1, 0.100460103605e1, -0.187200622541, 0.980276957749e-2, 0.543224866339e-1, -0.263496312610e-1, 0.315432315759e-1, -0.525788294155e-1, -0.685380627808e-2, 0.344540276656e-1, -0.555747275982e-3 ], "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.133442326203e1, 0.104116843433e1, -0.227202245707, -0.378598758038, 0.249888797892, -0.498847982876e-1 ], "d2": [0, 0, 0, 2, 2, 2], "t2": [3, 4, 5, 3, 4, 5], "c2": [2] * 6, "gamma2": [0.711139834571] * 6, "nr3": [], "nr4": [] } eq = helmholtz1, MBWR, GERG, helmholtz3 _PR = -0.004803 _surface = { "sigma": [-1.4165, 0.746383, 0.675625], "exp": [0.63882, 0.659804, 0.619149] } _dielectric = { "eq": 3, "Tref": 273.16, "rhoref": 1000., "a0": [], "expt0": [], "expd0": [], "a1": [2.0306, 0.0056], "expt1": [0, 1], "expd1": [1, 1], "a2": [0.181, 0.021, -7.4], "expt2": [0, 1, 0], "expd2": [2, 2, 3] } _melting = { "eq": 1, "Tref": Tt, "Pref": 7.3578, "Tmin": Tt, "Tmax": 400.0, "a1": [1], "exp1": [0], "a2": [5626.3, 2717.2], "exp2": [1, 1.83], "a3": [], "exp3": [] } _sublimation = { "eq": 3, "Tref": Tt, "Pref": 7.7, "Tmin": Tt, "Tmax": Tt, "a1": [], "exp1": [], "a2": [-8.065], "exp2": [0.93], "a3": [], "exp3": [] } _vapor_Pressure = { "eq": 5, "ao": [-0.489789e1, 0.988558, 0.349689, 0.499356], "exp": [1.0, 1.5, 2.0, 2.85] } _liquid_Density = { "eq": 1, "ao": [0.15456e2, -0.41720e2, 0.50276e2, -0.27947e2, 0.56718e1], "exp": [0.62, 0.83, 1.05, 1.3, 1.6] } _vapor_Density = { "eq": 3, "ao": [-0.29962e1, -0.16724e2, 0.15819e2, -0.16852e2, 0.34586e2, -0.53754e2], "exp": [0.466, 2, 2.4, 4., 7., 8.] } visco0 = { "eq": 0, "method": "_visco0", "__name__": "Muzny (2013)", "__doi__": { "autor": "Muzny, C.D., Huber, M.L., and Kazakov, A.F.", "title": "Correlation for the Viscosity of Normal Hydrogen Obtained from Symbolic Regression", "ref": "J. Chem. Eng. Data, 2013, 58 (4), pp 969–979", "doi": "10.1021/je301273j" } } def _visco0(self, rho, T, fase): sigma = 0.297 ek = 30.41 # Zero-Density Limit, Eq. 3-4 T_ = T / ek ai = [2.0963e-1, -4.55274e-1, 1.43602e-1, -3.35325e-2, 2.76981e-3] suma = 0 for i, a in enumerate(ai): suma += a * log(T_)**i S = exp(suma) no = 0.021357 * (self.M * T)**0.5 / sigma / S # Excess Contribution, Eq. 5-7 bi = [-0.187, 2.4871, 3.7151, -11.0972, 9.0965, -3.8292, 0.5166] B_ = 0 for i, b in enumerate(bi): B_ += b / T**i B = B_ * sigma**3 n1 = B * no # Simbolic Regression, Eq. 9 rhor = rho / 90.5 Tr = T / self.Tc c = [ 6.43449673, 4.56334068e-2, 2.32797868e-1, 9.5832612e-1, 1.27941189e-1, 3.63576595e-1 ] nc = c[0] * rhor**2 * exp(c[1] * Tr + c[2] / Tr + c[3] * rhor**2 / (c[4] + Tr) + c[5] * rhor**6) return unidades.Viscosity(no + n1 + nc, "muPas") visco1 = { "eq": 0, "method": "_visco1", "__name__": "McCarty (1972)", "__doi__": { "autor": "McCarty, R.D. and Weber, L.A.", "title": "Thermophysical properties of parahydrogen from the freezing liquid line to 5000 R for pressures to 10,000 Psia", "ref": "NBS Technical Note 617", "doi": "" } } def _visco1(self, rho, T, fase): DELV = lambda rho1, T1, rho2, T2: DILV(T1) + EXCESV(rho1, T2) \ - DILV(T2)-EXCESV(rho2, T2) def EXVDIL(rho, T): A = exp(5.7694 + log(rho.gcc) + 0.65e2 * rho.gcc**1.5 - 6e-6 * exp(127.2 * rho.gcc)) B = 10 + 7.2 * ( (rho.gcc / 0.07)**6 - (rho.gcc / 0.07)**1.5) - 17.63 * exp(-58.75 * (rho.gcc / 0.07)**3) return A * exp(B / T) * 0.1 def DILV(T): b = [ -0.1841091042788e2, 0.3185762039455e2, -0.2308233586574e2, 0.9129812714730e1, -0.2163626387630e1, 0.3175128582601, -0.2773173035271e-1, 0.1347359367871e-2, -0.2775671778154e-4 ] suma = 0 for i, b in enumerate(b): suma += b * T**((-3. + i) / 3) return suma * 100 def EXCESV(rho, T): c = [ -0.1324266117873e2, 0.1895048470537e2, 0.2184151514282e2, 0.9771827164811e5, -0.1157010275059e4, 0.1911147702539e3, -0.3186427506942e4, 0.0705565000000 ] R2 = rho.gcc**0.5 * (rho.gcc - c[7]) / c[7] A = c[0] + c[1] * R2 + c[2] * rho.gcc**0.1 + c[3] * R2 / T**2 + c[ 4] * rho.gcc**0.1 / T**1.5 + c[5] / T + c[6] * R2 / T B = c[0] + c[5] / T return 0.1 * (exp(A) - exp(B)) if T > 100: n = DILV(100) + EXVDIL(rho, 100) \ + DELV(rho, T, rho, 100) else: n = DILV(T) + EXVDIL(rho, T) return unidades.Viscosity(n, "muPas") visco2 = { "eq": 4, "omega": 1, "__name__": "Quiñones-Cisneros (2011)", "__doi__": { "autor": "S.E.Quinones-Cisneros, M.L. Huber and U.K. Deiters", "title": "model of 1-march-2011", "ref": "unpublished", "doi": "" }, "Tref": 33.145, "muref": 1.0, "ek": 59.7, "sigma": 0.2827, "n_chapman": 0, "n_ideal": [ 72.46400680522131e-1, -352.3929484813708e-1, 664.5332385860778e-1, -566.74979475607415e-1, 265.66570031561248e-1, -54.81307488054635e-1, 4.595978383724549e-1 ], "t_ideal": [0, 0.5, 0.75, 1, 1.25, 1.5], "n_poly": [1], "t_poly": [0.75], "n_polyden": [1, 1], "t_polyden": [0, 1], "nb": [ 1.0, -0.187, 75.6327, 3435.61, -312078, 7.77929e6, -9.95841e7, 4.08557e8 ], "tb": [0.0157768, 0, -1, -2, -3, -4, -5, -6], "a": [ -0.00002348389676311179e3, 0.00002197232806029717e3, 2.4547322430816313e-3, 3.9791170684039065e-8, 4.581319859008102e-3 ], "b": [ 0.000026869839733943842e3, 0.000027387647542474032e3, 0.000013065230652860072e3, 3.0723581102227345e-7, -0.00007033089468735152e3 ], "c": [0, 0, 0, 0, 0], "A": [ -3.912305916140789e-5, -2.1198288980972056e-6, 4.690087618888682e-6, 1.6938783854559677e-11, 9.39021777998824e-5 ], "B": [ -6.381148168720446e-5, 5.178086941554603e-4, -4.5508093750991845e-5 - 1.3780811004280076e-9 - 3.7679840470735697e-4 ], "C": [0, 0, 0, 0, 0], "D": [4.3699367404316626e-7, 0.0, -1.1321685281996792e-8, 0, 0] } visco3 = { "eq": 1, "omega": 1, "__name__": "Vargaftik (1996)", "__doi__": { "autor": "Vargaftik, N.B., Vinogradov, Y.K. and Yargin, V.S.", "title": "Handbook of Physical Properties of Liquids and Gases", "ref": "Hemisphere Publishing Corporation,New York, NY", "doi": "" }, "ek": 59.7, "sigma": 0.2827, "Tref": 32.938, "rhoref": 1. * M, "n_virial": [ -2.1505e-1, 10.727e-1, -16.935e-1, 0.0, 22.702e-1, 2.2123e-1, 0.34163e-1, -0.043206e-1 ], "t_virial": [-1.5, -1, -0.5, 0, 0.5, 1.5, 2.], "Tref_virial": 32.938, "etaref_virial": 1. * M, "Tref_res": 32.938, "rhoref_res": 15.556 * M, "etaref_res": 1., "n_packed": [], "t_packed": [], "n_poly": [ -9.22703e-1, 6.41602, -5.98018, 2.89715e-1, 2.36429, -2.78870e-1, -1.10595e1, 1.11582e1, 7.18928, -7.76971, -1.21827, 1.47193 ], "t_poly": [0, -1, -2, -3, 0, 0, -1, -2, -1, -2, -1, -2], "d_poly": [1, 1, 1, 1, 2, 3, 3, 3, 4, 4, 5, 5], "g_poly": [0] * 12, "c_poly": [0] * 12, "n_num": [], "t_num": [], "d_num": [], "g_num": [], "c_num": [], "n_den": [], "t_den": [], "d_den": [], "g_den": [], "c_den": [] } _viscosity = visco0, visco1, visco2, visco3 thermo0 = { "eq": 1, "__name__": "Assael (2011)", "__doi__": { "autor": " Assael, M.J., Assael. J.-A.M., Huber, M.L., Perkins, R.A. and Takata, Y.", "title": "Correlation of the Thermal Conductivity of Normal and Parahydrogen from the Triple Point to 1000 K and up to 100 MPa", "ref": "J. Phys. Chem. Ref. Data 40, 033101 (2011)", "doi": "10.1063/1.3606499" }, "__test__": """ >>> st=H2(T=298.15, rho=0) >>> print "%0.5g" % st.k.mWmK 185.67 >>> st=H2(T=298.15, rho=0.80844) >>> print "%0.5g" % st.k.mWmK 186.97 >>> st=H2(T=298.15, rho=14.4813) >>> print "%0.5g" % st.k.mWmK 201.35 >>> st=H2(T=35, rho=0) >>> print "%0.5g" % st.k.mWmK 26.988 >>> st=H2(T=35, rho=30) >>> print "%0.5g" % st.k.mWmK 75.594 >>> st=H2(T=35, rho=30) >>> print "%0.5g" % st.k.mWmK 71.854 >>> st=H2(T=18, rho=0) >>> print "%0.5g" % st.k.mWmK 13.875 >>> st=H2(T=18, rho=75) >>> print "%0.5g" % st.k.mWmK 104.48 """, # Table 4, Pag 8 "Tref": 1.0, "kref": 1e-3, "no": [ -1.24159e7, 5.04056e6, -4.80868e4, 3.26394e2, 9.56218e-2, 1.73488e-4, -3.12802e-8 ], "co": [0, 1, 2, 3, 4, 5, 6], "noden": [5.04305e6, -2.43753e4, 1.51523e2, 1.0], "coden": [0, 1, 2, 3], "Trefb": 33.145, "rhorefb": 15.508, "krefb": 1., "nb": [ .363081e-1, -.207629e-1, .31481e-1, -.143097e-1, .17498e-2, .18337e-2, -.886716e-2, .15826e-1, -.106283e-1, .280673e-2 ], "tb": [0, 0, 0, 0, 0, 1, 1, 1, 1, 1], "db": [1, 2, 3, 4, 5, 1, 2, 3, 4, 5], "cb": [0] * 10, "critical": 3, "gnu": 0.63, "gamma": 1.2415, "R0": 1.01, "Xio": 0.15e-9, "gam0": 0.052, "qd": 0.4e-9, "Tcref": 49.7175 } _thermal = thermo0,
class C1Cyclohexane(MEoS): """Multiparameter equation of state for methylcyclohexane""" name = "methylcyclohexane" CASNumber = "108-87-2" formula = "C6H11-CH3" synonym = "" _refPropName = "C1CC6" _coolPropName = "" rhoc = unidades.Density(267.0660832) Tc = unidades.Temperature(572.2) Pc = unidades.Pressure(3470.0, "kPa") M = 98.18606 # g/mol Tt = unidades.Temperature(146.7) Tb = unidades.Temperature(374.) f_acent = 0.23 momentoDipolar = unidades.DipoleMoment(0.0, "Debye") id = 39 CP1 = {"ao": 2.04122, "an": [0.016417, 0.000185315, -3.14826e-7, 1.65567e-10], "pow": [1, 2, 3, 4], "ao_exp": [], "exp": []} lemmon = { "__type__": "Helmholtz", "__name__": "short Helmholtz equation of state for methylcyclohexane " "of Lemmon (2007).", "__doi__": {"autor": "Lemmon, E.W.", "title": "", "ref": "unpublished equation, 2007", "doi": ""}, # Yoneda, Y., Sato, S., Matsumoto, T. # Density of Methylcyclohexane at Temperatures up to 600K and Pressures # up to 200 MPa # Int. J. Thermophys. 38 (2017) 106 # doi: 10.1007/s10765-017-2241-7 # As report Yoneda et al., its new (p, ρ, T) data disabling this eq as # reference equation "R": 8.314472, "cp": CP1, "ref": "NBP", "Tmin": Tt, "Tmax": 600., "Pmax": 500000.0, "rhomax": 9.13, "nr1": [1.3026, -2.6270, 0.68834, -0.16415, 0.092174, 0.0003842], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.38, 1.2, 2.14, 1.6, 0.3, 0.7], "nr2": [-0.29737, -0.078187, -0.049139, -0.30402, -0.074888], "d2": [1, 2, 5, 1, 4], "t2": [2.7, 3.25, 2.35, 3.7, 4.1], "c2": [1, 1, 1, 2, 2], "gamma2": [1]*17} eq = lemmon, _vapor_Pressure = { "eq": 3, "n": [-0.65871e1, -0.56553e1, 0.68947e1, -0.41281e1, -0.25444e1], "t": [1.0, 1.5, 1.6, 3.2, 10.]} _liquid_Density = { "eq": 1, "n": [0.18273e-1, 0.15215e2, -0.21951e2, 0.94466e1, 0.16781], "t": [0.1, 0.64, 0.8, 1.0, 4.5]} _vapor_Density = { "eq": 2, "n": [-0.52572e1, -0.13417e2, -0.24271e1, -0.54482e2, -0.15791e3], "t": [0.544, 2.3, 2.5, 6.1, 15.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, "rhoc": 267.07, "Toref": 572.2, "koref": 1, "no": [2.89968e-3, -1.80666e-2, 7.27576e-2, -1.29778e-2], "to": [0, 1, 2, 3], "Tref_res": 572.2, "rhoref_res": 2.72*M, "kref_res": 1., "nr": [9.19149e-2, -7.90408e-2, -8.17088e-2, 9.23911e-2, 2.96449e-2, -4.28498e-2, -2.99834e-3, 7.2786e-3], "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": 1.5e-10, "gam0": 0.052, "qd": 6.24e-10, "Tcref": 858.3} _thermal = thermo0,
class D5(MEoS): """Multiparameter equation of state for decamethylcyclopentasiloxane""" name = "decamethylcyclopentasiloxane" CASNumber = "541-02-6" formula = "C10H30O5Si5" synonym = "D5" rhoc = unidades.Density(292.570762680819) Tc = unidades.Temperature(619.23462341) Pc = unidades.Pressure(1160, "kPa") M = 370.7697 # g/mol Tt = unidades.Temperature(226.0) Tb = unidades.Temperature(484.05) f_acent = 0.658 momentoDipolar = unidades.DipoleMoment(1.349, "Debye") # id=1671 CP1 = { "ao": -34.898, "an": [1861.5e-3, -1403.4e-6, 500.0e-9], "pow": [1, 2, 3], "ao_exp": [], "exp": [], "ao_hyp": [], "hyp": [] } helmholtz1 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for hexamethyldisiloxane of Colonna et al. (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" }, "__test__": """ >>> st=D5(T=619.23462341, P=1161460) >>> print "%0.6f" % st.v 0.003418 """, # Table 22, Pag 206 "R": 8.314472, "cp": CP1, "ref": "NBP", "Tmin": 300.0, "Tmax": 673.0, "Pmax": 30000.0, "rhomax": 2.83, "Pmin": 0.000005, "rhomin": 2.83, "nr1": [ 1.40844725, -2.29248044, 0.42851607, -0.73506382, 0.16103808, 0.29643278e-3 ], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875], "nr2": [ 0.82412481, 0.15214274, -0.68495890, -0.55703624e-1, 0.13055391e-1, -0.31853761e-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 = helmholtz1, _vapor_Pressure = { "eq": 5, "ao": [-0.99967e1, 0.70091e1, -0.72265e1, -0.62938e1], "exp": [1.0, 1.5, 1.87, 3.8] } _liquid_Density = { "eq": 1, "ao": [0.303988e3, -0.110342e4, 0.134359e4, -0.705243e3, 0.164540e3], "exp": [0.57, 0.65, 0.73, 0.84, 0.96] } _vapor_Density = { "eq": 3, "ao": [ -0.37577e1, -0.47669e1, -0.24233e2, -0.29872e3, 0.34441e3, -0.32498e3 ], "exp": [0.459, 1.02, 2.6, 6.7, 7.7, 11.0] }
class NH3(MEoS): """Multiparameter equation of state for ammonia""" name = "ammonia" CASNumber = "7664-41-7" formula = "NH3" synonym = "R-717" _refPropName = "AMMONIA" _coolPropName = "Ammonia" rhoc = unidades.Density(225.) Tc = unidades.Temperature(405.40) Pc = unidades.Pressure(11333.0, "kPa") M = 17.03026 # g/mol Tt = unidades.Temperature(195.495) Tb = unidades.Temperature(239.823) f_acent = 0.25601 momentoDipolar = unidades.DipoleMoment(1.470, "Debye") id = 63 Fi1 = { "R": 8.314471, "ao_log": [1, -1], "pow": [0, 1, 1 / 3, -1.5, -1.75], "ao_pow": [-15.81502, 4.255726, 11.47434, -1.296211, 0.5706757] } CP2 = { "ao": 5.111814, "an": [ -0.42966650e2, -0.10243792e-1, 0.38750775e-4, -0.46406097e-7, 0.20268561e-10 ], "pow": [-1, 1, 2, 3, 4] } 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, "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., Baehr, H.D.", "title": "The Thermodynamic Properties of Ammonia", "ref": "VDI-Forsch., Number 596, 1979.", "doi": "" }, "R": 8.31434, "M": 17.03026, "Tc": 405.4, "Pc": 11333, "rhoc": 13.212, "cp": CP2, "ref": "IIR", "Tmin": 195.486, "Tmax": 600., "Pmax": 400000.0, "rhomax": 44.0, "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, "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., 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., "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 _PR = [0.1491, -26.0703] _surface = {"sigma": [0.1028, -0.09453], "exp": [1.211, 5.585]} _melting = { "eq": 1, "__doi__": { "autor": "Haar, L., Gallagher, J.S.", "title": "Thermodynamic Properties of Ammonia", "ref": "J. Phys. Chem. Ref. Data 7(3) (1978) 635-792", "doi": "10.1063/1.555579" }, "Tmin": Tt, "Tmax": 700.0, "Tref": Tt, "Pref": 101325, "a3": [2.5e4], "exp3": [1] } _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 Propylene(MEoS): """Multiparameter equation of state for propylene""" name = "propylene" CASNumber = "115-07-1" formula = "CH2=CH-CH3" synonym = "R-1270" rhoc = unidades.Density(230.08) Tc = unidades.Temperature(364.211) Pc = unidades.Pressure(4555.0, "kPa") M = 42.07974 # g/mol Tt = unidades.Temperature(87.953) Tb = unidades.Temperature(225.531) f_acent = 0.146 momentoDipolar = unidades.DipoleMoment(0.366, "Debye") id = 23 Fi1 = { "ao_log": [1, 3], "pow": [0, 1], "ao_pow": [-5.1823279651, 4.3639902765], "ao_exp": [1.544, 4.013, 8.923, 6.020], "titao": [324 / Tc, 973 / Tc, 1932 / Tc, 4317 / Tc] } Fi2 = { "ao_log": [1, 3.07317535], "pow": [0, 1], "ao_pow": [9.48120502357782, -4.47976952867319], "ao_exp": [1.7018443, 3.61342025, 8.83689058, 6.27183616], "titao": [1.01164134251849, 2.75278088800174, 5.16557061703243, 11.68984352477] } CP1 = { "ao": 0.65591381, "an": [ 0.44359641e-1, -.36650786e-4, 0.16822223e-7, -.32651013e-11, 0.33747826e4 ], "pow": [1, 2, 3, 4, -2], "ao_exp": [-4.7032420], "exp": [615.8], "ao_hyp": [], "hyp": [] } helmholtz1 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for propylene of Lemmon et al. (2013).", "__doi__": { "autor": "Lemmon, E.W., Overhoff, U., McLinden, M.O., Wagner, W.", "title": "A reference equation of state for the thermodynamic properties of propene for temperatures from the melting line to 575 K and pressures up to 1000 MPa", "ref": "to be submitted to J. Phys. Chem. Ref. Data", "doi": "" }, "R": 8.314472, "cp": Fi1, "ref": "IIR", "Tmin": Tt, "Tmax": 575.0, "Pmax": 1000000.0, "rhomax": 23.1, "Pmin": 0.00000075, "rhomin": 18.255, "nr1": [ 0.4341002e-1, 0.1136592e1, -0.8528611, 0.5216669, -0.1382953e1, 0.1214347 ], "d1": [4, 1, 1, 2, 2, 3], "t1": [1.0, 0.205, 0.56, 0.676, 1.0, 0.5], "nr2": [ -0.5984662, -0.1391883e1, -0.1008434e1, 0.1961249, -0.3606930, -0.2407175e-2 ], "d2": [1, 1, 3, 2, 2, 8], "t2": [1.0, 1.94, 2.0, 1.0, 2.66, 0.83], "c2": [1, 2, 2, 1, 2, 1], "gamma2": [1] * 6, "nr3": [ 0.7432121, 0.1475162, -0.2503391e-1, -0.2734409, 0.6378889e-2, 0.1502940e-1, -0.3162971e-1, -0.4107194e-1, -0.1190241e1 ], "d3": [1, 1, 2, 3, 3, 2, 1, 2, 3], "t3": [1.6, 2.5, 3.0, 2.5, 2.72, 4.0, 4.0, 1.0, 4.0], "alfa3": [1.07, 0.66, 1.2, 1.12, 1.47, 1.93, 3.3, 15.4, 6], "beta3": [0.77, 0.83, 0.607, 0.4, 0.66, 0.07, 3.1, 387, 41], "gamma3": [1.21, 1.08, 0.83, 0.56, 1.22, 1.81, 1.54, 1.12, 1.4], "epsilon3": [0.78, 0.82, 1.94, 0.69, 1.96, 1.3, 0.38, 0.91, 0.7] } helmholtz2 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for propylene of Overhoff (2006).", "__doi__": { "autor": "Overhoff, U.", "title": "Development of a new equation of state for the fluid region of propene for temperatures from the melting line to 575 K with pressures to 1000 MPa as well as software for the computation of thermodynamic properties of fluids", "ref": "Ph.D. Dissertation, Ruhr University, Bochum, Germany, 2006.", "doi": "" }, "R": 8.314472, "cp": Fi2, "ref": "IIR", "Tmin": Tt, "Tmax": 575.0, "Pmax": 1000000.0, "rhomax": 23.4, "Pmin": 0.00000074, "rhomin": 18.26, "nr1": [ 0.11167427541961e1, -0.76114879497376, -0.18654354344883e1, 0.41500701892893e-1, 0.10706545719025e-1, 0.17481482892991e-1 ], "d1": [1, 1, 1, 3, 4, 4], "t1": [0.125, 0.625, 1.25, 0, 0.25, 1.25], "nr2": [ 0.56509607629258, 0.99156795771235, -0.16341922173416, -0.37037920319844e-1, -0.80058345775777e-1, 0.17004662808796, 0.81351262137108e-1, -0.23817885171378, 0.12962562859214e-1, 0.22577442976798e2, -0.43611886043491e2, 0.21944325628071e2, -0.66234078215924, -0.22258580712469e1, 0.29538388307646e1, -0.10257185828694e1, 0.20521625234481e-1, -0.36462809205891e-1, 0.17625833164005e-1 ], "d2": [2, 3, 3, 3, 4, 4, 5, 5, 6, 1, 1, 1, 1, 2, 2, 2, 5, 6, 1], "t2": [ 2.25, 1.25, 2.125, 2.75, 0.125, 2, 1.125, 1.5, 1.375, 3.5, 3.75, 4, 5, 3, 3.5, 4.5, 4.75, 3.25, 3, ], "c2": [1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3], "gamma2": [1] * 19, "nr3": [ 0.31819374579431, -0.32648950998998, -0.37684374593786e2, 0.72265437094447e2, -0.34814669335983e2, -0.39854778355193e1, 0.37313453915501 ], "d3": [2, 2, 1, 1, 1, 2, 2], "t3": [3, 4, 2, 3, 4, 1, 1], "alfa3": [10, 10, 11, 11, 11, 25, 30], "beta3": [150, 150, 225, 225, 225, 300, 350], "gamma3": [1.13, 1.13, 1.19, 1.19, 1.19, 1.19, 1.19], "epsilon3": [0.85, 0.85, 1, 1, 1, 1, 1] } helmholtz3 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for propylene of Angus et al. (1980)", "__doi__": { "autor": "Angus, S., Armstrong, B., and de Reuck, K.M.", "title": "International Thermodynamic Tables of the Fluid State-7 Propylene", "ref": "International Union of Pure and Applied Chemistry, Pergamon Press, Oxford, 1980.", "doi": "" }, "R": 8.31434, "cp": CP1, "ref": "IIR", "Tmin": 100.0, "Tmax": 600.0, "Pmax": 200000.0, "rhomax": 9.73, "Pmin": 0.48475e-4, "rhomin": 17.938, "nr1": [ 0.631922681460, 0.102655250604, -0.70798923e-2, 0.18624829, -0.1292611017e1, -0.5410160974e-1, 0.5069017035, -0.10606146125e1, 0.763136083, -0.850733053e-1, 0.438262575, 0.2316495716e-1, 0.25503741325e-1, -0.57327581, -0.1141334722e-1, 0.2502895522, -0.468392547833e-1, 0.325228355714e-2 ], "d1": [0, 0, 0, 1, 1, 1, 2, 2, 2, 3, 3, 4, 4, 4, 4, 5, 6, 7], "t1": [3, 4, 5, 1, 2, 3, 0, 1, 2, 2, 3, 0, 1, 3, -1, 3, 3, 3], "nr2": [ -0.63192268146, -0.102655250604, 0.70798923e-2, -0.63192268146, -0.102655250604, -0.11049992895, -0.31596134073, -0.51327625302e-1, -0.4918627871e-1, -0.17109208434e-1, -0.1492467645e-1, -0.42773021085e-2, -0.8554604217e-3, -0.14257673695e-3 ], "d2": [0, 0, 0, 2, 2, 2, 4, 4, 6, 6, 8, 8, 10, 12], "t2": [3, 4, 5, 3, 4, 5, 3, 4, 3, 4, 3, 4, 4, 4], "c2": [2] * 14, "gamma2": [1] * 14 } eq = helmholtz1, helmholtz2, helmholtz3 _surface = {"sigma": [0.05268], "exp": [1.186]} _melting = { "eq": 1, "Tref": Tt, "Pref": 0.48475e-4, "Tmin": Tt, "Tmax": 2000.0, "a1": [-6593000000, 6593000001], "exp1": [0, 2.821], "a2": [], "exp2": [], "a3": [], "exp3": [] } _vapor_Pressure = { "eq": 5, "ao": [-6.75625, 2.02700, -1.35883, -2.74671, -0.936445], "exp": [1.0, 1.5, 1.9, 4.3, 15.0] } _liquid_Density = { "eq": 1, "ao": [0.405430, 2.02481, 0.304022, 0.179159], "exp": [0.195, 0.47, 2.25, 8.0] } _vapor_Density = { "eq": 3, "ao": [-1.59841, -4.73840, -10.8886, -31.0312, -56.9431, -143.544], "exp": [0.309, 0.853, 2.37, 5.2, 10., 20.] }
class iC4(MEoS): """Multiparameter equation of state for isobutane""" name = "isobutane" CASNumber = "75-28-5" formula = "CH(CH3)3" synonym = "R-600a" rhoc = unidades.Density(225.5) Tc = unidades.Temperature(407.81) Pc = unidades.Pressure(3629.0, "kPa") M = 58.1222 # g/mol Tt = unidades.Temperature(113.73) Tb = unidades.Temperature(261.401) f_acent = 0.184 momentoDipolar = unidades.DipoleMoment(0.132, "Debye") id = 5 _Tr = unidades.Temperature(390.355535) _rhor = unidades.Density(228.302484) _w = 0.178714317 Fi1 = { "ao_log": [1, 3.05956619], "pow": [0, 1], "ao_pow": [11.60865546, -5.29450411], "ao_exp": [4.94641014, 4.09475197, 15.6632824, 9.73918122], "titao": [0.9512779015, 2.3878958853, 4.3469042691, 10.3688586351], "ao_hyp": [], "hyp": [] } Fi2 = { "ao_log": [1, 3.06714], "pow": [0, 1], "ao_pow": [20.413726078, -94.467620036], "ao_exp": [], "titao": [], "ao_hyp": [8.97575, 5.25156, 25.1423, 16.1388], "hyp": [1.074673199, 0.485556021, 4.671261865, 2.19158348] } Fi3 = { "ao_log": [1, 3.059347], "pow": [0, 1], "ao_pow": [-5.404217, 4.91136], "ao_exp": [4.940314, 4.090139, 15.68832, 9.739581], "titao": [0.9508183, 2.383449, 10.38655, 4.347095], "ao_hyp": [], "hyp": [] } CP4 = { "ao": -1.7231723278e1, "an": [ 1.7027919006e7, -4.7269724737e5, 4.7301406581e3, 5.8491344291e-2, 8.9440351886e-6, -1.8274599197e-8 ], "pow": [-3, -2, -1.001, 1, 2, 3], "ao_exp": [-1.9283021962e1], "exp": [3000], "ao_hyp": [], "hyp": [] } CP5 = { "ao": 4.06714, "an": [], "pow": [], "ao_exp": [], "exp": [], "ao_hyp": [0.1724067e7, 0.2059196e6, 0.9124395e8, 0.1289193e8], "hyp": [0.4382700e3, 0.1980180e3, 0.1905020e4, 0.8937650e3] } CP6 = { "ao": 0.397893 / 8.3143 * 58.124, "an": [ 0.412501e-2 / 8.3143 * 58.124, -0.196195e-6 / 8.3143 * 58.124, 0.380185e-8 / 8.3143 * 58.124, -0.523950e-11 / 8.3143 * 58.124 ], "pow": [1, 2, 3, 4], "ao_exp": [], "exp": [], "ao_hyp": [], "hyp": [] } helmholtz1 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for isobutane of Buecker and Wagner (2006)", "__doi__": { "autor": "Bücker, D., Wagner, W.", "title": "Reference Equations of State for the Thermodynamic Properties of Fluid Phase n-Butane and Isobutane", "ref": "J. Phys. Chem. Ref. Data 35, 929 (2006)", "doi": "10.1063/1.1901687" }, "__test__": # Table 46, Pag 996 """ >>> st=iC4(T=113.730, x=0.5) >>> print "%0.6g %0.8f %0.5f %0.7f %0.5g %0.5g %0.3f %0.3f %0.3f %0.3f %0.3f %0.3f %0.2f %0.2f" % (\ 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) 113.73 0.00000002 740.33876 0.0000014 -714.04 -233.34 -3.199 1.028 1.174 0.737 1.689 0.880 1999.81 139.39 >>> st=iC4(T=130, x=0.5) >>> print "%0.6g %0.8f %0.4f %0.6f %0.2f %0.2f %0.3f %0.3f %0.3f %0.3f %0.3f %0.3f %0.2f %0.2f" % (\ 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) 130 0.00000088 724.9759 0.000047 -686.14 -218.43 -2.970 0.628 1.209 0.810 1.741 0.953 1863.41 147.92 >>> st=iC4(T=170, x=0.5) >>> print "%0.6g %0.6f %0.3f %0.5f %0.2f %0.2f %0.3f %0.3f %0.3f %0.3f %0.3f %0.3f %0.2f %0.2f" % (\ 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) 170 0.000274 686.932 0.01127 -613.96 -177.00 -2.486 0.084 1.293 0.975 1.868 1.118 1578.65 166.95 >>> st=iC4(T=180, x=0.5) >>> print "%0.6g %0.6f %0.3f %0.5f %0.2f %0.2f %0.3f %0.3f %0.3f %0.3f %0.3f %0.3f %0.2f %0.2f" % (\ 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) 180 0.000737 677.281 0.02867 -595.12 -165.68 -2.379 0.007 1.315 1.014 1.900 1.158 1514.00 171.29 >>> st=iC4(T=300, x=0.5) >>> print "%0.6g %0.5f %0.2f %0.4f %0.2f %0.2f %0.3f %0.3f %0.3f %0.3f %0.3f %0.3f %0.2f %0.2f" % (\ 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 0.37000 548.32 9.6096 -338.17 -11.29 -1.299 -0.209 1.690 1.578 2.442 1.810 810.25 197.74 >>> st=iC4(T=400, x=0.5) >>> print "%0.6g %0.4f %0.2f %0.2f %0.2f %0.2f %0.3f %0.3f %0.3f %0.3f %0.3f %0.3f %0.2f %0.2f" % (\ 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) 400 3.1856 341.03 118.39 -36.27 81.59 -0.459 -0.164 2.250 2.354 6.349 7.555 184.38 128.90 >>> st=iC4(T=405, x=0.5) >>> print "%0.6g %0.4f %0.2f %0.2f %0.2f %0.2f %0.3f %0.3f %0.3f %0.3f %0.3f %0.3f %0.2f %0.2f" % (\ 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) 405 3.4629 305.04 148.28 -9.33 71.55 -0.394 -0.194 2.362 2.481 13.190 18.178 138.43 119.17 >>> st=iC4(T=406, x=0.5) >>> print "%0.6g %0.4f %0.2f %0.2f %0.2f %0.2f %0.3f %0.3f %0.3f %0.3f %0.3f %0.3f %0.2f %0.2f" % (\ 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) 406 3.5210 293.52 158.36 -2.12 67.07 -0.376 -0.206 2.402 2.521 19.368 27.552 128.12 116.78 >>> st=iC4(T=407, x=0.5) >>> print "%0.6g %0.4f %0.2f %0.2f %0.2f %0.2f %0.3f %0.3f %0.3f %0.3f %0.3f %0.3f %0.2f %0.2f" % (\ 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) 407 3.5801 276.83 173.46 7.30 59.64 -0.354 -0.225 2.465 2.575 42.319 61.185 116.97 113.83 """ # Table 47, Pag 1003 """ >>> st=iC4(T=200, P=1e5) >>> print "%0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.6g" % (\ st.T, st.rho, st.u.kJkg, st.h.kJkg, st.s.kJkgK, st.cv.kJkgK, st.cp.kJkgK, st.w) 200 657.78 -556.49 -556.34 -2.1751 1.3627 1.9676 1389.68 >>> st=iC4(T=400, 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) 400 9.2132 129.91 184.18 0.31057 2.0152 2.1959 236.58 >>> st=iC4(T=400, 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) 400 19.586 122.43 173.49 0.19221 2.0358 2.2711 224.3 >>> st=iC4(T=350, 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) 350 476.78 -209.39 -206.24 -0.89941 1.9081 2.8832 530.58 >>> st=iC4(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 46.082 103.55 146.95 0.042963 2.0943 2.5614 194.12 >>> st=iC4(T=400, P=3e6) >>> 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.405 69.772 101.22 -0.11041 2.2417 4.2367 145.97 >>> st=iC4(T=400, P=4e6) >>> 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 374.74 -59.311 -48.637 -0.4951 2.1788 4.0279 270.3 >>> st=iC4(T=575, P=4e6) >>> 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) 575 54.649 522.22 595.42 0.88371 2.7352 3.0054 268.45 >>> st=iC4(T=400, 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) 400 394.05 -67.804 -55.115 -0.51779 2.1541 3.5415 328.2 >>> st=iC4(T=400, P=1e7) >>> 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 439.11 -89.303 -66.53 -0.57611 2.1246 3.0033 486.23 >>> st=iC4(T=500, P=3.5e7) >>> 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) 500 433.69 145.98 226.69 -0.051328 2.522 3.0667 628.75 """, "R": 8.314472, "cp": Fi1, "ref": "OTO", "Tmin": Tt, "Tmax": 575.0, "Pmax": 35000.0, "rhomax": 12.9, "Pmin": 0.0000219, "rhomin": 12.74, "nr1": [ 0.20686820727966e1, -0.36400098615204e1, 0.51968754427244, 0.17745845870123, -0.12361807851599, 0.45145314010528e-1, 0.30476479965980e-1 ], "d1": [1, 1, 1, 2, 3, 4, 4], "t1": [0.50, 1.00, 1.50, 0.00, 0.50, 0.50, 0.75], "nr2": [ 0.75508387706302, -0.85885381015629, 0.36324009830684e-1, -0.19548799450550e-1, -0.44452392904960e-2, 0.46410763666460e-2, -0.71444097992825e-1, -0.80765060030713e-1, 0.15560460945053, 0.20318752160332e-2, -0.10624883571689, 0.39807690546305e-1, 0.16371431292386e-1, 0.53212200682628e-3, -0.78681561156387e-2, -0.30981191888963e-2 ], "d2": [1, 1, 2, 7, 8, 8, 1, 2, 3, 3, 4, 5, 5, 10, 2, 6], "t2": [ 2.00, 2.50, 2.50, 1.50, 1.00, 1.50, 4.00, 7.00, 3.00, 7.00, 3.00, 1.00, 6.00, 0.00, 6.00, 13.00 ], "c2": [1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3], "gamma2": [1] * 16, "nr3": [-0.42276036810382e-1, -0.53001044558079e-2], "d3": [1, 2], "t3": [2., 0.], "alfa3": [10, 10], "beta3": [150, 200], "gamma3": [1.16, 1.13], "epsilon3": [0.85, 1.] } MBWR = { "__type__": "MBWR", "__name__": "MBWR equation of state for isobutane of Younglove and Ely (1987)", "__doi__": { "autor": "Younglove, B.A. and Ely, J.F.", "title": "Thermophysical Properties of Fluids. II. Methane, Ethane, Propane, Isobutane, and Normal Butane ", "ref": "J. Phys. Chem. Ref. Data 16, 577 (1987)", "doi": "10.1063/1.555785" }, "R": 8.31434, "cp": CP4, "ref": { "Tref": 298.15, "Pref": 101.325, "ho": 17932.6, "so": 295.390 }, "Tmin": 113.55, "Tmax": 600.0, "Pmax": 35000.0, "rhomax": 12.89, "Pmin": 1.948e-5, "rhomin": 12.755, "b": [ None, 0.1307325972e-1, 0.3927802742, -0.3185427394e2, 0.7608825192e4, -0.1753919859e7, -0.2090019755e-2, 0.8959557971e1, -0.6816710130e4, -0.1111271045e7, 0.3248737572e-3, -0.1046526456e1, 0.6536598969e3, 0.3726503734e-1, 0.8553649395e1, 0.2109987236e4, -0.1401267363e1, 0.5213089327e-1, -0.1925026382e2, 0.7640067895, 0.3425854273e7, -0.3373475924e9, 0.1180683444e6, 0.1529683738e10, 0.3323837416e4, 0.6423169487e5, 0.3891706042e2, -0.1494755736e7, -0.1720240173e-1, 0.2894195375e3, 0.2005086329e-2, -0.4448393005, 0.8028488415e2 ] } GERG = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for isobutane 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, 2012, 57 (11), pp 3032-3091", "doi": "10.1021/je300655b" }, "R": 8.314472, "cp": Fi2, "ref": "OTO", "Tmin": Tt, "Tmax": 575.0, "Pmax": 35000.0, "rhomax": 12.9, # "Pmin": 7.36, "rhomin": 38.2, "nr1": [ 0.10429331589100e1, -0.28184272548892e1, 0.86176232397850, -0.10613619452487, 0.98615749302134e-1, 0.23948208682322e-3 ], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875], "nr2": [ 0.30330004856950, -0.41598156135099e-1, -0.29991937470058, -0.80369342764109e-1, -0.29761373251151e-1, 0.13059630303140e-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 } helmholtz3 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for isobutane of Miyamoto and Watanabe (2001)", "__doi__": { "autor": "Miyamoto, H. and Watanabe, K.", "title": "A Thermodynamic Property Model for Fluid-Phase Isobutane", "ref": "Int. J. Thermophys., 23(2):477-499, 2002.", "doi": "10.1023/A:1015161519954" }, "R": 8.314472, "cp": Fi3, "ref": "IIR", "Tmin": 113.56, "Tmax": 573.0, "Pmax": 35000.0, "rhomax": 12.9, "Pmin": 0.000021, "rhomin": 12.738, "nr1": [ 2.892737e-1, -1.342570, -7.976713e-3, 2.025793e-1, -4.241612e-2, 2.617971e-3, 5.068955e-5, -1.144596e-6 ], "d1": [1, 1, 2, 2, 3, 5, 8, 8], "t1": [-0.25, 1.5, -0.75, 0, 1.25, 1.5, 0.5, 2.5], "nr2": [ -1.930153, 1.982609, 2.076533e-3, -4.958752e-3, 1.377372e-3, -1.582662e-1, -4.961892e-2, 9.451030e-4, -3.037276e-2, -1.382675e-2, 8.876254e-5 ], "d2": [3, 3, 8, 5, 6, 1, 5, 7, 2, 3, 15], "t2": [1.5, 1.75, -0.25, 3, 3, 4, 2, -1, 2, 19, 5], "c2": [21, 1, 1, 1, 1, 2, 2, 2, 3, 3, 3], "gamma2": [1] * 11 } helmholtz4 = { "__type__": "Helmholtz", "__name__": "short Helmholtz equation of state for isobutane 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=iC4(T=700, rho=200, eq=4) >>> print "%0.4f %0.3f %0.4f" % (st.cp0.kJkgK, st.P.MPa, st.cp.kJkgK) 3.2392 19.108 3.5575 >>> st2=iC4(T=750, rho=100, eq=4) >>> print "%0.2f %0.5f" % (st2.h.kJkg-st.h.kJkg, st2.s.kJkgK-st.s.kJkgK) 210.32 0.37469 """, # Table III, Pag 46 "R": 8.31451, "cp": CP5, "Tmin": 113.55, "Tmax": 600.0, "Pmax": 100000.0, "rhomax": 12.89, "Pmin": 0.000020860, "rhomin": 12.784, "nr1": [ 0.10429332e1, -0.28184273e1, 0.86176232, -0.10613619, 0.986157490e-1, 0.23948209e-3 ], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875], "nr2": [ 0.30330005, -0.41598156e-1, -0.29991937, -0.80369343e-1, -0.29761373e-1, 0.1305963e-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 } helmholtz5 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for isobutane 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": CP6, "ref": "NBP", "Tmin": 120.0, "Tmax": 498.0, "Pmax": 35000.0, "rhomax": 12.89, "Pmin": 0.46491e-4, "rhomin": 12.649, "nr1": [ -0.958589873652, 0.818846326211, -0.115814967179, 0.345513148715, -0.168751721524e1, 0.936693300209, -0.106644545724e1, 0.980958295776e-1, 0.495941129005, -0.261313404262, 0.485109471188, -0.177275820736, -0.209415485311e-1, 0.788178884079e-1, -0.102751671767, 0.178645875838e-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.958589873652, -0.818846326211, 0.115814967179, 0.537585249054, -0.71942446879, 0.245830118086 ], "d2": [0, 0, 0, 2, 2, 2], "t2": [3, 4, 5, 3, 4, 5], "c2": [2] * 6, "gamma2": [1.0071072] * 6 } helmholtz6 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for isobutane 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, "Tmin": Tt, "Tmax": 620.0, "Pmax": 800000.0, "rhomax": 40., "Pmin": 0.1, "rhomin": 40., "nr1": [ 1.18083775, 9.46903331e-1, -2.90618044, 8.51346220e-2, 2.79868503e-4, -1.68266335e-1 ], "d1": [1, 1, 1, 3, 7, 2], "t1": [1.5, 0.25, 1.25, 0.25, 0.875, 1.375], "nr2": [ -2.01202825e-1, -3.32570120e-2, 2.42967225e-1, -4.20931100e-3, -2.24528572e-1, -1.41307663e-2, -5.93401702e-2, -2.27862942e-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, MBWR, GERG, helmholtz3, helmholtz4, helmholtz5, helmholtz6 _surface = {"sigma": [-0.01639, 0.06121], "exp": [2.102, 1.304]} _dielectric = { "eq": 3, "Tref": 273.16, "rhoref": 1000., "a0": [0.388417], "expt0": [-1.], "expd0": [1.], "a1": [20.534, 0.02], "expt1": [0, 1], "expd1": [1, 1], "a2": [126.25, 52.91, -7501.4, -2672.9], "expt2": [0, 1, 0, 1], "expd2": [2, 2, 2.9, 2.9] } _melting = { "eq": 1, "Tref": Tt, "Pref": 0.000022891, "Tmin": Tt, "Tmax": 575.0, "a1": [-1953637129., 1953637130.], "exp1": [0, 6.12], "a2": [], "exp2": [], "a3": [], "exp3": [] } _vapor_Pressure = { "eq": 5, "ao": [-6.85093103, 1.36543198, -1.32542691, -2.56190994], "exp": [1, 1.5, 2.5, 4.5] } _liquid_Density = { "eq": 2, "ao": [2.04025104, 0.850874089, -0.479052281, 0.348201252], "exp": [1.065, 3, 4, 7] } _vapor_Density = { "eq": 6, "ao": [-2.12933323, -2.93790085, -0.89441086, -3.46343707], "exp": [1.065, 2.5, 9.5, 13] } visco0 = { "eq": 1, "omega": 1, "collision": [0.53583008, -0.45629630, 0.049911282], "__name__": "Vogel (2000)", "__doi__": { "autor": "Vogel, E., Kuechenmeister, C., and Bich, E.", "title": "Viscosity Correlation for Isobutane over Wide Ranges of the Fluid Region", "ref": "Int. J. Thermophys, 21(2):343-356, 2000.", "doi": "10.1023/A:1006623310780" }, "ek": 307.55, "sigma": 0.46445, "Tref": 1., "rhoref": 1. * M, "n_chapman": 0.1628213 / M**0.5, "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_virial": 307.55, "etaref_virial": 0.0603345, "Tref_res": 407.817, "rhoref_res": 3.86 * M, "etaref_res": 1, "n_packed": [0.233859774637e1, 0.235255150838e1], "t_packed": [0, 0.5], "n_poly": [ 0.103511763411e3, -0.312670896234e3, 0.145253750239e3, -0.210649894193e3, 0.386269696509e3, -0.214963015527e3, 0.112580360920e3, -0.223242033154e3, 0.119114788598e3, -0.181909745900e2, 0.360438957232e2, -0.213960184050e2, -0.194037606990e4 ], "t_poly": [0, -1, -2, -0, -1, -2, 0, -1, -2, 0 - 1, -2, 0], "d_poly": [2, 2, 2, 3, 3, 3, 4, 4, 4, 5, 5, 5, 1], "g_poly": [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1], "c_poly": [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], "n_num": [0.194037606990e4], "t_num": [0], "d_num": [1], "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] } visco1 = { "eq": 2, "omega": 2, "__name__": "Younglove (1987)", "__doi__": { "autor": "Younglove, B.A. and Ely, J.F.", "title": "Thermophysical Properties of Fluids. II. Methane, Ethane, Propane, Isobutane, and Normal Butane ", "ref": "J. Phys. Chem. Ref. Data 16, 577 (1987)", "doi": "10.1063/1.555785" }, "ek": 418.0, "sigma": 0.509217, "n_chapman": 0.203525266 / M**0.5, "F": [1.687838652, 0.0, 1.40, 407.85], "E": [ -0.2055498053e2, 0.1357076181e4, 0.1893774336e2, -0.1822277344e5, -0.4599387773e-2, 0.6305247065e2, 0.1282253921e5 ], "rhoc": 3.86 } _viscosity = visco0, visco1 thermo0 = { "eq": 1, "__name__": "Perkins (2002)", "__doi__": { "autor": "Perkins, R.A.", "title": "Measurement and Correlation of the Thermal Conductivity of Isobutane from 114 K to 600 K at Pressures to 70 MPa", "ref": "J. Chem. Eng. Data, 2002, 47 (5), pp 1272–1279", "doi": "10.1021/je010121u" }, "Tref": 407.85, "kref": 1, "no": [-2.37901e-3, 1.06601e-2, 2.15811e-2], "co": [0, 1, 2], "Trefb": 407.85, "rhorefb": 3.86, "krefb": 1, "nb": [ -4.11789e-2, 4.76346e-2, 1.46805e-1, -1.28445e-1, -1.19190e-1, 1.07565e-1, 4.10226e-2, -3.85968e-2, -4.88704e-3, 5.20901e-3 ], "tb": [0, 1] * 5, "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.657661e-9, "Tcref": 611.73 } thermo1 = { "eq": 2, "omega": 2, "__name__": "Younglove (1987)", "__doi__": { "autor": "Younglove, B.A. and Ely, J.F.", "title": "Thermophysical Properties of Fluids. II. Methane, Ethane, Propane, Isobutane, and Normal Butane ", "ref": "J. Phys. Chem. Ref. Data 16, 577 (1987)", "doi": "10.1063/1.555785" }, "visco": visco1, "n_chapman": 2.0352526600e-1, "G": [0.1449797353e1, -0.1685643887], "E": [ 0.4307008989e-2, -0.1509010974e1, 0.4693712392e3, -0.3554280979e-3, 0.1841552874, -0.3892338766e2, -0.9354624917e-1, 0.7114330590e1 ], "critical": 2, "X": [0.0034718, 10.1207, 0.466392, 1.00344], "Z": 9.10218e-10 } _thermal = thermo0, thermo1
class nC4(MEoS): """Multiparameter equation of state for n-butane""" name = "n-butane" CASNumber = "106-97-8" formula = "CH3-(CH2)2-CH3" synonym = "R-600" _refPropName = "BUTANE" _coolPropName = "n-Butane" rhoc = unidades.Density(228.) Tc = unidades.Temperature(425.125) Pc = unidades.Pressure(3796.0, "kPa") M = 58.1222 # g/mol Tt = unidades.Temperature(134.895) Tb = unidades.Temperature(272.660) f_acent = 0.201 momentoDipolar = unidades.DipoleMoment(0.05, "Debye") id = 6 _Tr = unidades.Temperature(406.785141) _rhor = unidades.Density(230.384826) _w = 0.194240287 Fi1 = {"ao_log": [1, 3.24680487], "pow": [0, 1], "ao_pow": [12.54882924, -5.46976878], "ao_exp": [5.54913289, 11.4648996, 7.59987584, 9.66033239], "titao": [0.7748404445, 3.3406025522, 4.9705130961, 9.9755537783]} Fi2 = {"ao_log": [1, 3.33944], "pow": [0, 1], "ao_pow": [20.884143364, -91.638478026], "ao_exp": [], "titao": [], "ao_sinh": [9.44893, 24.4618], "sinh": [468.27/Tc, 1914.1/Tc], "ao_cosh": [6.89406, 14.7824], "cosh": [183.636/Tc, 903.185/Tc]} Fi3 = {"ao_log": [1, 3.240207], "pow": [0, 1], "ao_pow": [-5.404217, 4.91136], "ao_exp": [5.513671, 7.388450, 10.250630, 11.061010], "titao": [327.55988/Tc, 1319.06935/Tc, 4138.63184/Tc, 1864.36783/Tc]} CP4 = {"ao": -1.3491511376e1, "an": [3.8802310194e5, -1.5444296890e5, 2.8455082239e3, 6.6142595353e-2, -2.4307965028e-5, 1.5044248429e-10], "pow": [-3, -2, -1, 1, 2, 3], "ao_exp": [-8.3933423467], "exp": [3000]} CP6 = {"ao": 0.801601/8.3143*58.124, "an": [0.655936e-3/8.3143*58.124, 0.12277e-4/8.3143*58.124, -0.165626e-7/8.3143*58.124, 0.67736e-11/8.3143*58.124], "pow": [1, 2, 3, 4], "ao_exp": [], "exp": []} buecker = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for butane of Buecker and " "Wagner (2006)", "__doi__": {"autor": "Bücker, D., Wagner, W.", "title": "Reference Equations of State for the " "Thermodynamic Properties of Fluid Phase " "n-Butane and Isobutane", "ref": "J. Phys. Chem. Ref. Data 35(2) (2006) 929-1019", "doi": "10.1063/1.1901687"}, "R": 8.314472, "cp": Fi1, "ref": "OTO", "Tmin": Tt, "Tmax": 750., "Pmax": 200000.0, "rhomax": 13.86, "nr1": [0.25536998241635e1, -0.44585951806696e1, 0.82425886369063, 0.11215007011442, -0.35910933680333e-1, 0.16790508518103e-1, 0.32734072508724e-1], "d1": [1, 1, 1, 2, 3, 4, 4], "t1": [0.50, 1.00, 1.50, 0.00, 0.50, 0.50, 0.75], "nr2": [0.95571232982005, -0.10003385753419e1, 0.85581548803855e-1, -0.025147918369616, -0.15202958578918e-2, 0.47060682326420e-2, -0.97845414174006e-1, -0.48317904158760e-1, 0.17841271865468, 0.18173836739334e-1, -0.11399068074953, 0.19329896666669e-1, 0.11575877401010e-2, 0.15253808698116e-3, -0.43688558458471e-1, -0.82403190629989e-2], "d2": [1, 1, 2, 7, 8, 8, 1, 2, 3, 3, 4, 5, 5, 10, 2, 6], "t2": [2.00, 2.50, 2.50, 1.50, 1.00, 1.50, 4.00, 7.00, 3.00, 7.00, 3.00, 1.00, 6.00, 0.00, 6.00, 13.00], "c2": [1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3], "gamma2": [1]*16, "nr3": [-0.28390056949441e-1, 0.14904666224681e-2], "d3": [1, 2], "t3": [2., 0.], "alfa3": [10, 10], "beta3": [150, 200], "gamma3": [1.16, 1.13], "epsilon3": [0.85, 1.]} younglove = { "__type__": "MBWR", "__name__": "MBWR equation of state for butane of Younglove and Ely " "(1987)", "__doi__": {"autor": "Younglove, B.A., Ely, J.F.", "title": "Thermophysical Properties of Fluids. II. " "Methane, Ethane, Propane, Isobutane, and Normal " "Butane", "ref": "J. Phys. Chem. Ref. Data 16(4) (1987) 577-798", "doi": "10.1063/1.555785"}, "R": 8.31434, "M": 58.125, "Tt": 134.86, "Tc": 425.16, "Pc": 3796, "rhoc": 3.92, "cp": CP4, "ref": {"Tref": 300, "Pref": 101.325, "ho": 19208.9, "so": 309.95}, "Tmin": 134.86, "Tmax": 600., "Pmax": 70000.0, "rhomax": 13.2, "b": [None, 0.153740104603e-1, -0.160980034611, -0.979782459010e1, 0.499660674504e3, -0.102115607687e7, 0.236032147756e-2, -0.137475757093e1, -0.907038733865e3, 0.385421748213e6, -0.349453710700e-4, 0.157361122714, 0.102301474068e3, 0.182335737331e-1, -0.404114307787e1, 0.187979855783e1, 0.362088795040, -0.738762248266e-2, -0.218618590563e1, 0.118802729027, 0.706854198713e6, -0.219469885796e9, -0.182454361268e5, 0.206790377277e10, 0.111757550145e3, 0.558779925986e5, -0.159579054026e2, -0.148034214622e7, -0.245206328201, 0.218305259309e3, -0.923990627338e-4, -0.205267776639e1, 0.387639044820e2]} GERG = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for butane 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": Fi2, "ref": "OTO", "Tmin": Tt, "Tmax": 575., "Pmax": 69000.0, "rhomax": 13.2, "nr1": [0.10626277411455e1, -0.28620951828350e1, 0.88738233403777, -0.12570581155345, 0.10286308708106, 0.25358040602654e-3], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875], "nr2": [0.32325200233982, -0.037950761057432, -0.32534802014452, -0.079050969051011, -0.020636720547775, 0.57053809334750e-2], "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} miyamoto = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for butane of Miyamoto and " "Watanabe (2001)", "__doi__": {"autor": "Miyamoto, H., Watanabe, K.", "title": "A Thermodynamic Property Model for Fluid-Phase " "n-Butane", "ref": "Int. J. Thermophys., 22(2) (2001) 459-475", "doi": "10.1023/A:1010722814682"}, "R": 8.314472, "cp": Fi3, "ref": "IIR", "Tmin": 134.87, "Tmax": 589., "Pmax": 69000.0, "rhomax": 13.15, "nr1": [2.952054e-1, -1.32636, -2.031317e-3, 2.240301e-1, -3.635425e-2, 1.905841e-3, 7.409154e-5, -1.401175e-6], "d1": [1, 1, 2, 2, 3, 5, 8, 8], "t1": [-0.25, 1.5, -0.75, 0, 1.25, 1.5, 0.5, 2.5], "nr2": [-2.492172, 2.386920, 1.424009e-3, -9.393388e-3, 2.616590e-3, -1.977323e-1, -3.809534e-2, 1.523948e-3, -2.391345e-2, -9.535229e-3, 3.928384e-5], "d2": [3, 3, 8, 5, 6, 1, 5, 7, 2, 3, 15], "t2": [1.5, 1.75, -0.25, 3, 3, 4, 2, -1, 2, 19, 5], "c2": [1, 1, 1, 1, 1, 2, 2, 2, 3, 3, 3], "gamma2": [1]*11} shortSpan = { "__type__": "Helmholtz", "__name__": "short Helmholtz equation of state for butane 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": Fi2, "ref": "OTO", "M": 58.123, "rhoc": 227.84/58.123, "Tmin": 134.86, "Tmax": 750., "Pmax": 100000.0, "rhomax": 13.20, "nr1": [0.10626277e1, -0.28620952e1, 0.88738233, -0.12570581, 0.10286309, 0.25358041e-3], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875], "nr2": [0.323252, -0.37950761e-1, -0.32534802, -0.79050969e-1, -0.20636721e-1, 0.57053809e-2], "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 butane 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": CP6, "ref": "NBP", "Tmin": 140.0, "Tmax": 589., "Pmax": 30000.0, "rhomax": 12.81, "nr1": [-0.504188295325, 0.541067401063, -0.760421383062e-1, 0.846035653528, -0.191317317203e1, 0.521441860186, -0.783511318207, 0.689697797175e-1, 0.947825461055e-1, -0.141401831669, 0.382675021672, -0.423893176684e-1, 0.677591792029e-1, 0.567943363340e-1, -0.131517698401, 0.221136942526e-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.504188295325, -0.541067401063, 0.760421383062e-1, -0.619109535460e-1, 0.423035373804, -0.390505508895], "d2": [0, 0, 0, 2, 2, 2], "t2": [3, 4, 5, 3, 4, 5], "c2": [2]*6, "gamma2": [1.08974964]*6} sun = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for butane 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.314472, "cp": Fi1, "ref": "OTO", "Tmin": Tt, "Tmax": 620.0, "Pmax": 800000.0, "rhomax": 40., "nr1": [1.18936994, 1.05407451, -3.24964532, 8.25263908e-2, 2.76467405e-4, -8.09869214e-2], "d1": [1, 1, 1, 3, 7, 2], "t1": [1.5, 0.25, 1.25, 0.25, 0.875, 1.375], "nr2": [-9.38097492e-2, 1.46213532e-1, 4.01168502e-1, -1.28716120e-2, -0.275191070, -1.62708971e-2, -7.04082962e-2, -2.32871995e-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 = buecker, younglove, GERG, miyamoto, shortSpan, polt, sun _PR = [-0.1332, -15.8278] _surface = {"sigma": [0.05138], "exp": [1.209]} _dielectric = { "eq": 1, "a": [20.611, 0.020], "b": [66.64, 24.44], "c": [-7461.2, -1983.6], "Au": 15.23, "D": 2} _melting = { "eq": 1, "__doi__": buecker["__doi__"], "Tmin": 134.895, "Tmax": 575.0, "Tref": Tt, "Pref": 0.653, "a0": 1, "a2": [5.585582364e8], "exp2": [2.206]} _vapor_Pressure = { "eq": 3, "n": [-0.71897e1, 0.26122e1, -0.21729e1, -0.27230e1], "t": [1, 1.5, 2., 4.5]} _liquid_Density = { "eq": 1, "n": [0.52341e1, -0.62011e1, 0.36063e1, 0.22137], "t": [0.44, 0.6, 0.76, 5.0]} _vapor_Density = { "eq": 2, "n": [-0.27390e1, -0.57347e1, -0.16408e2, -0.46986e2, -0.10090e3], "t": [0.39, 1.14, 3.0, 6.5, 14.0]} visco0 = {"__name__": "Herrmann (2018)", "__doi__": { "autor": "Herrmann, S., Vogel, E.", "title": "New Formulation for the Viscosity of n-Butane", "ref": "J. Phys. Chem. Ref. Data 47(1) (2018) 013104", "doi": "10.1063/1.5020802"}, "eq": 1, "omega": 0, "special0": "_mu0", "Tref_virial": 425.125, # Special term of virial coefficient, with δ term and μPa·s "muref_virial": 4.89736312734e-1/228*M/1e3, "n_virial": [-1.9572881000e1, 1.98887362343e2, -8.3176420912e2, 1.83218450345e3, -2.26510439059e3, 1.51348864395e3, -4.32819866497e2, 5.19698852489, -3.86579291550e-2], "t_virial": [0, -0.25, -0.5, -0.75, -1, -1.25, -1.5, -2.5, -5.5], "Tref_res": 425.125, "rhoref_res": 228, "nr": [2.3460864383872, 7.8632175809804e-1, 1.5823593499816e1, -9.4670516989296, 1.0511496276340, -1.9355799491084e-2, 1.4895031937816e-4], "tr": [2, 5, 0, 0, 0, 4, 5], "dr": [2, 2, 2.5, 3, 5, 7.5, 10], "nr_gaus": [1.2790911462043, 2.5581822924086e-1], "br_gaus": [30, 5], "er_gaus": [220, 400], "special": "_mur"} def _mu0(self, T): """Special term for zero-density viscosity for Herrmann correlation""" tau = self.Tc/T # Eq 8 no = [4.6147656002208, 4.574318591039e-1, 3.0851104723224e-2] suma = 0 for i, n in enumerate(no): suma += n*log(tau)**i muo = 1.0546549635209e3/tau**0.5/exp(suma) return muo def _mur(self, rho, T, fase): """Special exponential term of residual viscosity for Herrmann correlation""" tau = self.Tc/T delta = rho/self.rhoc mur = tau**0.5/delta**(2/3)*1.2280342363570e-3*(delta**5.7*tau)**2 return mur visco1 = {"__name__": "Vogel (1999)", "__doi__": { "autor": "Vogel, E., Küchenmeister, C., Bich, E.", "title": "Viscosity correlation for n-Butane in the Fluid " "Region", "ref": "High Temp. - High Pressures 31(2) (1999) 173-186", "doi": "10.1068/htrt154"}, "eq": 1, "omega": 1, "ek": 280.51, "sigma": 0.57335, "n_chapman": 0.021357, "collision": [0.17067154, -0.48879666, 0.039038856], "Tref_virial": 280.51, "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": 425.125, "rhoref_res": 3.92*M, "nr": [-54.7737770846, 58.0898623034, 35.2658446259, -39.6682203832, -1.83729542151, -0.833262985358, 1.93837020663], "tr": [0, 1, 0, 1, 0, 0, 1], "dr": [2, 2, 3, 3, 4, 5, 5], "CPf": 188.075903903, "CPg1": 2.30873963359, "CPgi": [0.88101765264], "CPti": [-0.5]} visco2 = {"__name__": "Younglove (1987)", "__doi__": { "autor": "Younglove, B.A., Ely, J.F.", "title": "Thermophysical Properties of Fluids. II. Methane, " "Ethane, Propane, Isobutane, and Normal Butane", "ref": "J. Phys. Chem. Ref. Data 16(4) (1987) 577-798", "doi": "10.1063/1.555785"}, "eq": 2, "omega": 2, "ek": 440., "sigma": 0.503103, "F": [0.1630521851e1, 0.0, 1.40, 425.16], "E": [-0.2724386845e2, 0.8012766611e3, 0.2503978646e2, -0.1309704275e5, -0.8313305258e-1, 0.6636975027e2, 0.9849317662e4], "rhoc": 3.920} visco3 = {"__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": 425.125, "no": [18.3983, -57.1255, 49.3197], "to": [0, 0.25, 0.5], "a": [-1.34111e-5, -8.56588e-5, 0], "b": [1.49860e-4, -1.71134e-4, 0], "c": [3.53018e-7, -1.93040e-5, 0], "A": [-3.63389e-9, -7.73717e-10, 0], "B": [3.70980e-8, 2.07659e-9, 0], "C": [-1.12496e-7, 7.66906e-8, 0]} _viscosity = visco0, visco1, visco2, visco3 thermo0 = {"__name__": "Perkins (2002)", "__doi__": { "autor": "Perkins, R.A, Ramires, M.L.V., Nieto de Castro, " "C.A., Cusco, L.", "title": "Measurement and Correlation of the Thermal " "Conductivity of Butane from 135 K to 600 K at " "Pressures to 70 MPa", "ref": "J. Chem. Eng. Data 47(5) (2002) 1263-1271", "doi": "10.1021/je0101202"}, "eq": 1, "Toref": 425.16, "koref": 1., "no": [1.62676e-3, 9.75703e-4, 2.89887e-2], "to": [0, 1, 2], "Tref_res": 425.16, "rhoref_res": 3.92*M, "kref_res": 1., "nr": [-3.04337e-2, 4.18357e-2, 1.65820e-1, -1.47163e-1, -1.48144e-1, 1.33542e-1, 5.25500e-2, -4.85489e-2, -6.29367e-3, 6.44307e-3], "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.03, "Xio": 0.194e-9, "gam0": 0.0496, "qd": 0.875350e-9, "Tcref": 637.68} thermo1 = {"__name__": "Younglove (1987)", "__doi__": { "autor": "Younglove, B.A., Ely, J.F.", "title": "Thermophysical Properties of Fluids. II. Methane," " Ethane, Propane, Isobutane, and Normal Butane", "ref": "J. Phys. Chem. Ref. Data 16(4) (1987) 577-798", "doi": "10.1063/1.555785"}, "eq": 3, "ek": 440., "G": [0.1530992335e1, -0.2114511021], "E": [0.4024170074e-2, 0.1561435847e1, -0.6004381127e3, -0.7547260841e-3, -0.2069676662e-1, 0.9382534978e2, -0.1711371457, 0.3647724935e2], "critical": 2, "Tc": 425.16, "rhoc": 3.92*58.125, "X": [0.000769608, 13.2533, 0.485554, 1.01021], "Z": 9.10218e-10} _thermal = thermo0, thermo1
class RE347mcc(MEoS): """Multiparameter equation of state for RE347mcc""" name = "methyl-heptafluoropropyl-ether" CASNumber = "375-03-1" formula = "CF3CF2CF2OCH3" synonym = "HFE-7000" _refPropName = "RE347MCC" _coolPropName = "" rhoc = unidades.Density(524.143687088) Tc = unidades.Temperature(437.7) Pc = unidades.Pressure(2476.2, "kPa") M = 200.0548424 # g/mol Tt = unidades.Temperature(250) Tb = unidades.Temperature(307.349) f_acent = 0.411 momentoDipolar = unidades.DipoleMoment(3.13, "Debye") CP1 = {"ao": 13.09, "an": [], "pow": [], "ao_exp": [], "exp": [], "ao_hyp": [13.78, 14.21, 0, 0], "hyp": [2045, 850, 0, 0]} zhou = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for RE347mcc of Zhou (2012)", "__doi__": {"autor": "Zhou, Y., Lemmon, E.W., Mahmoud, A.M.", "title": "Equations of state for RE245cb2, RE347mcc, " "RE245fa2 and R1216", "ref": "Preliminary equation", "doi": ""}, "R": 8.314472, "cp": CP1, "ref": "NBP", "Tmin": Tt, "Tmax": 500.0, "Pmax": 20000.0, "rhomax": 7.662, "Pmin": 6.825, "rhomin": 7.66, "nr1": [0.0330627, 2.606165, -4.902937, 2.228012, 1.494115, -2.420459, 0.160067], "d1": [4, 1, 1, 1, 2, 2, 3], "t1": [1, 0.34, 0.77, 1.02, 0.79, 1.017, 0.634], "nr2": [1.383893, -2.092005, -0.5904708], "d2": [2, 1, 2], "t2": [1.35, 2.25, 2.5], "c2": [1, 2, 2], "gamma2": [1]*6, "nr3": [-0.701794, 2.765425, 0.6860982, -2.208170, 0.1739594, -0.9028007, -0.0213123], "d3": [1, 1, 2, 2, 3, 3, 1], "t3": [2, 1.66, 1.33, 2.0, 1.87, 1.75, 1.05], "alfa3": [0.593, 1.36, 1.73, 1.483, 0.617, 1.596, 9.64], "beta3": [0.0872, 1.176, 1.53, 0.78, 0.088, 1.04, 263.0], "gamma3": [1.06, 1.22, 0.92, 1.08, 1.21, 0.85, 1.12], "epsilon3": [1.12, 0.79, 1.055, 0.5, 0.84, 0.85, 0.91]} eq = zhou, _vapor_Pressure = { "eq": 3, "n": [-8.0456, 2.6285, -2.7498, -5.4277, -4.3693], "t": [1.0, 1.5, 2., 4.25, 12.8]} _liquid_Density = { "eq": 1, "n": [1.5144, 2.3745, -2.6363, 2.0830, 0.50537], "t": [0.29, 0.85, 1.5, 2.2, 9.]} _vapor_Density = { "eq": 2, "n": [-2.0640, -6.4226, -18.982, -58.689, -117.64, -253.93], "t": [0.321, 0.96, 2.75, 5.9, 12., 22.]} 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.411, "mur": 0.0, "k": 0.0}
class Cyclohexane(MEoS): """Multiparameter equation of state for cyclohexane""" name = "cyclohexane" CASNumber = "110-82-7" formula = "cyclo(CH2)6" synonym = "" rhoc = unidades.Density(271.33016352) Tc = unidades.Temperature(553.6) Pc = unidades.Pressure(4080.5, "kPa") M = 84.15948 # g/mol Tt = unidades.Temperature(279.47) Tb = unidades.Temperature(353.865) f_acent = 0.2096 momentoDipolar = unidades.DipoleMoment(0.3, "Debye") id = 38 _Tr = unidades.Temperature(526.231121) _rhor = unidades.Density(274.647526) _w = 0.221837522 Fi1 = {"ao_log": [1, 3], "pow": [0, 1], "ao_pow": [0.9891140602, 1.6359660572], "ao_exp": [0.83775, 16.036, 24.636, 7.1715], "titao": [773/Tc, 941/Tc, 2185/Tc, 4495/Tc], "ao_hyp": [], "hyp": []} CP1 = {"ao": 9.3683272, "an": [-0.56214088e8, 0.15261554e-1, -0.36352468e-5], "pow": [-3, 1, 2], "ao_exp": [.23766589e2], "exp": [2000], "ao_hyp": [], "hyp": []} helmholtz1 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for cyclohexane of Zhou et al. (2014)", "__doi__": {"autor": "Zhou, Y., Jun Liu, J., Penoncello, S.G., Lemmon, E.W.", "title": "An Equation of State for the Thermodynamic Properties of Cyclohexane", "ref": "J. Phys. Chem. Ref. Data 43, 043105 (2014)", "doi": "10.1063/1.4900538"}, "__test__": """ >>> st=Cyclohexane(T=300, rhom=9.4) >>> print "%0.1f %0.1f %0.8g %0.8g %0.8g %0.8g %0.8g %0.8g" % (\ st.T, st.rhoM, st.P.MPa, st.cvM.JmolK, st.cpM.JmolK, st.w.ms, st.hM.Jmol, st.sM.JmolK) 300.0 9.4 24.173705 115.286 154.76956 1383.3878 -8400.0834 -28.889069 >>> st=Cyclohexane(T=500, rhom=6.5) >>> print "%0.1f %0.1f %0.8g %0.8g %0.8g %0.8g %0.8g %0.8g" % (\ st.T, st.rhoM, st.P.MPa, st.cvM.JmolK, st.cpM.JmolK, st.w.ms, st.hM.Jmol, st.sM.JmolK) 500.0 6.5 3.9246630 192.52056 255.57087 434.13064 31070.127 70.891447 >>> st=Cyclohexane(T=500, rhom=0.7) >>> print "%0.1f %0.1f %0.8g %0.8g %0.8g %0.8g %0.8g %0.8g" % (\ st.T, st.rhoM, st.P.MPa, st.cvM.JmolK, st.cpM.JmolK, st.w.ms, st.hM.Jmol, st.sM.JmolK) 500.0 0.7 1.9981172 191.96446 235.52281 155.348 52757.706 122.92657 >>> st=Cyclohexane(T=600, rhom=3.5) >>> print "%0.1f %0.1f %0.8g %0.8g %0.8g %0.8g %0.8g %0.8g" % (\ st.T, st.rhoM, st.P.MPa, st.cvM.JmolK, st.cpM.JmolK, st.w.ms, st.hM.Jmol, st.sM.JmolK) 600.0 3.5 6.8225506 232.79222 388.55185 150.53318 70150.132 143.42323 >>> st=Cyclohexane(T=553.6, rhom=3.3) >>> print "%0.1f %0.1f %0.8g %0.8g %0.8g %0.8g %0.8g %0.8g" % (\ st.T, st.rhoM, st.P.MPa, st.cvM.JmolK, st.cpM.JmolK, st.w.ms, st.hM.Jmol, st.sM.JmolK) 553.6 3.3 4.0805433 224.19555 199224.62 87.913911 58532.604 123.59810 >>> st=Cyclohexane(P=101325, x=0) >>> print "%0.9g %0.8g %0.8g %0.8g %0.8g %0.8g %0.5f %0.5f" % (\ st.T, st.rhoM, st.P.MPa, st.cvM.JmolK, st.cpM.JmolK, st.w.ms, st.hM.Jmol, st.sM.JmolK) 353.864939 8.5487851 0.101325 134.6163 179.07223 994.05862 0.00000 -0.00000 >>> st=Cyclohexane(P=101325, x=1) >>> print "%0.9g %0.8g %0.8g %0.8g %0.8g %0.8g %0.8g %0.8g" % (\ st.T, st.rhoM, st.P.MPa, st.cvM.JmolK, st.cpM.JmolK, st.w.ms, st.hM.Jmol, st.sM.JmolK) 353.864939 0.035779032 0.101325 123.4305 133.35895 186.91349 29991.286 84.753484 """, # Table 5, Pag 17 "R": 8.3144621, "cp": Fi1, "ref": {"Tref": 300, "Pref": 1., "ho": 23949.01, "so": 104.2926004}, "Tmin": 279.86, "Tmax": 700.0, "Pmax": 250000.0, "rhomax": 10.3, "Pmin": 5.2402, "rhomin": 9.403, "nr1": [0.05483581, 1.607734, -2.375928, -0.5137709, 0.1858417], "d1": [4, 1, 1, 2, 3], "t1": [1, 0.37, 0.79, 1.075, 0.37], "nr2": [-0.9007515, -0.5628776, 0.2903717, -0.3279141, -0.03177644], "d2": [1, 3, 2, 2, 7], "t2": [2.4, 2.5, 0.5, 3, 1.06], "c2": [2, 2, 1, 2, 1], "gamma2": [1]*5, "nr3": [0.8668676, -0.1962725, -0.1425992, 0.004197016, 0.1776584, -0.04433903, -0.03861246, 0.07399692, 0.02036006, 0.00272825], "d3": [1, 1, 3, 3, 2, 2, 3, 2, 3, 2], "t3": [1.6, 0.37, 1.33, 2.5, 0.9, 0.5, 0.73, 0.2, 1.5, 1.5], "alfa3": [0.99, 1.43, 0.97, 1.93, 0.92, 1.27, 0.87, 0.82, 1.4, 3], "beta3": [0.38, 4.2, 1.2, 0.9, 1.2, 2.6, 5.3, 4.4, 4.2, 25], "gamma3": [0.65, 0.63, 1.14, 0.09, 0.56, 0.4, 1.01, 0.45, 0.85, 0.86], "epsilon3": [0.73, 0.75, 0.48, 2.32, 0.2, 1.33, 0.68, 1.11, 1.47, 0.99]} helmholtz2 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for cyclohexane of Penoncello et al. (1995)", "__doi__": {"autor": "Penoncello, S.G., Goodwin, A.R.H., and Jacobsen, R.T.", "title": "A Thermodynamic Property Formulation for Cyclohexane", "ref": "Int. J. Thermophys., 16(2):519-531, 1995.", "doi": "10.1007/BF01441918"}, "R": 8.31434, "cp": CP1, "ref": {"Tref": 279.47, "Pref": 101.325, "ho": 33884.8, "so": 96.612}, "Tt": 279.47, "Tc": 553.64, "rhoc": 3.244, "M": 84.1608, "Tmin": 279.47, "Tmax": 700.0, "Pmax": 80000.0, "rhomax": 9.77, "Pmin": 5.2538, "rhomin": 9.4045, "nr1": [0.8425412659, -0.3138388327e1, 0.1679072631e1, -0.153819249, 0.1984911143, -0.144532594, 0.3746346428e-3, 0.1861479616e-3, 0.1745721652e-3], "d1": [1, 1, 1, 2, 3, 3, 7, 6, 6], "t1": [0, 1.5, 2.5, 1.5, 1, 2.5, 2, 0.5, 3], "nr2": [-0.6427428062, 0.2280757615, -0.1868116802e1, -0.1028243711e1, 0.5821457418, -0.255891152, 0.1276844113e-1, -0.5158613166e-2, 0.6334794755e-1, -0.6014686589e-1, 0.4439056828, -0.6264920642, 0.2132589969e1, -0.3620300991e-2, 0.2534453992, 0.1669144715e-1, 0.3985052291e-2], "d2": [1, 1, 2, 3, 3, 5, 8, 10, 3, 4, 1, 1, 2, 2, 4, 4, 8], "t2": [5, 6, 5.5, 3, 7, 6, 6.5, 5.5, 11, 11, 0.5, 1, 4, 4, 1.5, 2, 0.5], "c2": [2, 2, 2, 2, 2, 2, 2, 2, 4, 4, 3, 3, 2, 6, 2, 4, 2], "gamma2": [1]*17} helmholtz3 = { "__type__": "Helmholtz", "__name__": "short Helmholtz equation of state for cyclohexane 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=Cyclohexane(T=700, rho=200, eq=2) >>> print "%0.4f %0.3f %0.4f" % (st.cp0.kJkgK, st.P.MPa, st.cp.kJkgK) 3.0278 9.007 3.5927 >>> st2=Cyclohexane(T=750, rho=100, eq=2) >>> print "%0.2f %0.5f" % (st2.h.kJkg-st.h.kJkg, st2.s.kJkgK-st.s.kJkgK) 206.82 0.31448 """, # Table III, Pag 46 "R": 8.31451, "cp": CP1, "ref": {"Tref": 279.47, "Pref": 101.325, "ho": 33884.8, "so": 96.612}, "Tt": 279.47, "Tc": 553.64, "rhoc": 3.244, "M": 84.1608, "Tmin": Tt, "Tmax": 600.0, "Pmax": 100000.0, "rhomax": 9.77, "Pmin": 5.2428, "rhomin":9.3999, "nr1": [0.10232354e1, -0.29204964e1, 0.10736630e1, -0.19573985, 0.12228111, 0.28943321e-3], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875], "nr2": [0.27231767, -0.4483332e-1, -0.38253334, -0.89835333e-1, -0.24874965e-1, 0.10836132e-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} helmholtz4 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for cyclohexane 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.31434, "cp": CP1, "ref": {"Tref": 279.47, "Pref": 101.325, "ho": 33884.8, "so": 96.612}, "Tmin": Tt, "Tmax": 620.0, "Pmax": 800000.0, "rhomax": 40., "Pmin": 0.1, "rhomin": 40., "nr1": [1.27436292, 1.15372124, -3.86726473, 8.84627298e-2, 2.76478090e-4, 7.26682313e-2], "d1": [1, 1, 1, 3, 7, 2], "t1": [1.5, 0.25, 1.25, 0.25, 0.875, 1.375], "nr2": [7.10849914e-2, 4.46376742e-1, 7.64476190e-1, -4.23520282e-2, -3.96468623e-1, -1.41250071e-2, -1.08371284e-1, -2.50082884e-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, helmholtz4 _surface = {"sigma": [0.06485], "exp": [1.263]} _melting = {"eq": 1, "Tref": 1, "Pref": 700, "Tmin": Tt, "Tmax": 370.0, "a1": [0.1329969885, -374.255624], "exp1": [1.41, 0], "a2": [], "exp2": [], "a3": [], "exp3": []} _vapor_Pressure = { "eq": 5, "ao": [-7.0342, 1.7311, -1.7572, -3.3406], "exp": [1., 1.5, 2.3, 4.6]} _liquid_Density = { "eq": 1, "ao": [5.5081, -14.486, 38.241, -64.589, 57.919, -20.55], "exp": [0.51, 0.94, 1.4, 1.9, 2.4, 3.0]} _vapor_Density = { "eq": 3, "ao": [-3.69006, -41.4239, 220.914, -443.72, 491.49, -296.373], "exp": [0.446, 1.98, 2.75, 3.3, 4.1, 4.8]}
class R13I1(MEoS): """Multiparameter equation of state for trifluoroiodomethane""" name = "trifluoroiodomethane" CASNumber = "2314-97-8" formula = "CF3I" synonym = "R13I1" _refPropName = "CF3I" _coolPropName = "R13I1" rhoc = unidades.Density(868.00061824) Tc = unidades.Temperature(396.44) Pc = unidades.Pressure(3953., "kPa") M = 195.9104 # g/mol Tt = unidades.Temperature(120.) Tb = unidades.Temperature(251.3) f_acent = 0.18 momentoDipolar = unidades.DipoleMoment(0.92, "Debye") CP1 = {"ao": 4., "an": [], "pow": [], "ao_exp": [6.2641], "exp": [694]} lemmon = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for R13I1 of Lemmon and " "Span (2013)", "__doi__": { "autor": "Lemmon, E.W., Span, R.", "title": "Thermodynamic Properties of R-227ea, R-365mfc, " "R-115, and R-13I1", "ref": "J. Chem. Eng. Data, 60(12) (2015) 3745-3758", "doi": "10.1021/acs.jced.5b00684" }, "R": 8.3144621, "cp": CP1, "ref": "IIR", "Tmin": Tt, "Tmax": 420., "Pmax": 20000.0, "rhomax": 14.1, "nr1": [1.12191, -3.08087, 1.11307, -0.184885, 0.110971, 0.000325], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875], "nr2": [0.333357, -0.0288288, -0.371554, -0.0997985, -0.0333205, 0.0207882], "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, _surface = {"sigma": [0.05767], "exp": [1.298]} _vapor_Pressure = { "eq": 3, "n": [-6.8642, 1.7877, -1.0619, -2.1677], "t": [1.0, 1.5, 1.9, 3.8] } _liquid_Density = { "eq": 1, "n": [2.0711, 1.562, -2.599, 1.7177], "t": [0.38, 1.3, 1.9, 2.5] } _vapor_Density = { "eq": 2, "n": [-3.0987, -6.8771, -19.701, -46.86, -100.02], "t": [0.41, 1.33, 3.5, 7.4, 16.0] }
class R141b(MEoS): """Multiparameter equation of state for R141b""" name = "1,1-dichloro-1-fluoroethane" CASNumber = "1717-00-6" formula = "CCl2FCH3" synonym = "R141b" _refPropName = "R141B" _coolPropName = "R141b" rhoc = unidades.Density(458.55946) Tc = unidades.Temperature(477.5) Pc = unidades.Pressure(4212.0, "kPa") M = 116.94962 # g/mol Tt = unidades.Temperature(169.68) Tb = unidades.Temperature(305.20) f_acent = 0.2195 momentoDipolar = unidades.DipoleMoment(2.014, "Debye") # id = 1633 Fi1 = { "ao_log": [1, 3.], "pow": [0, 1], "ao_pow": [-15.5074814985, 9.1871858933], "ao_exp": [6.8978, 7.8157, 3.2039], "titao": [502 / Tc, 1571 / Tc, 4603 / Tc] } lemmon = { "__type__": "Helmholtz", "__name__": "short Helmholtz equation of state for R-141b 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": 500.0, "Pmax": 400000.0, "rhomax": 12.56, "nr1": [1.1469, -3.6799, 1.3469, 0.083329, 0.00025137], "d1": [1, 1, 1, 3, 7], "t1": [0.25, 1.25, 1.5, 0.25, 0.875], "nr2": [ 0.32720, 0.46946, -0.029829, -0.31621, -0.026219, -0.078043, -0.020498 ], "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.1122, -17.5406] _surface = {"sigma": [7.3958e-5, 0.059941], "exp": [0.066331, 1.2214]} _vapor_Pressure = { "eq": 3, "n": [-0.73784e1, 0.52955e1, -0.46639e1, -0.31122e1, -0.18972e1], "t": [1.0, 1.5, 1.7, 4.2, 9.0] } _liquid_Density = { "eq": 1, "n": [0.10443e2, -0.24726e2, 0.27718e2, -0.11220e2, 0.75848], "t": [0.49, 0.68, 0.88, 1.1, 2.9] } _vapor_Density = { "eq": 2, "n": [ -0.31177e1, -0.68872e1, -0.18566e2, -0.40311e2, -0.95472e1, -0.12482e3 ], "t": [0.398, 1.33, 3.3, 6.7, 7.0, 14.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": 370.44, "sigma": 0.5493, "omega": 5, "psi": [0.921345, 4.1091e-2], "psi_d": [0, 1], "fint": [5.21722e-4, 2.92456e-6], "fint_t": [0, 1], "chi": [1.0867, -2.16469e-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 Benzene(MEoS): """Multiparameter equation of state for benzene""" name = "benzene" CASNumber = "71-43-2" formula = "C6H6" synonym = "" _refPropName = "BENZENE" _coolPropName = "Benzene" rhoc = unidades.Density(304.7922436) Tc = unidades.Temperature(562.02) Pc = unidades.Pressure(4894, "kPa") M = 78.11184 # g/mol Tt = unidades.Temperature(278.674) Tb = unidades.Temperature(353.22) f_acent = 0.211 momentoDipolar = unidades.DipoleMoment(0.0, "Debye") id = 40 Fi1 = {"ao_log": [1, 2.94645], "pow": [0, 1], "ao_pow": [-0.6740687105, 2.5560188958], "ao_exp": [7.36374, 18.6490, 4.01834], "titao": [4116/Tc, 1511/Tc, 630/Tc]} CP2 = {"ao": -0.478176/8.3143*78.108, "an": [0.618649e-2/8.3143*78.108, -0.380363e-5/8.3143*78.108, 0.699648e-9/8.3143*78.108, 0.42661e-13/8.3143*78.108], "pow": [1, 2, 3, 4], "ao_exp": [], "exp": []} thol = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for benzene of Thol (2013).", "__doi__": {"autor": "Thol, M., Lemmon, E.W., Span, R.", "title": "Equation of State for Benzene for Temperatures " "from the Melting Line up to 750 K with Pressures" " up to 500 MPa", "ref": "High Temperatures-High Pressures 41 (2012) 81-97", "doi": ""}, "R": 8.314472, "cp": Fi1, "ref": "NBP", "Tmin": Tt, "Tmax": 725, "Pmax": 500000.0, "rhomax": 11.45, "nr1": [0.03513062, 2.229707, -3.100459, -0.5763224, 0.2504179], "d1": [4, 1, 1, 2, 3], "t1": [1, 0.3, 0.744, 1.174, 0.68], "nr2": [-0.7049091, -0.1393433, 0.8319673, -0.3310741, -0.02793578], "d2": [1, 3, 2, 2, 7], "t2": [2.5, 3.67, 1.26, 2.6, 0.95], "c2": [2, 2, 1, 2, 1], "gamma2": [1.]*5, "nr3": [0.7087408, -0.3723906, -0.06267414, -0.86295], "d3": [1, 1, 3, 3], "t3": [1, 2.47, 3.35, 0.75], "alfa3": [1.032, 1.423, 1.071, 14.35], "beta3": [1.867, 1.766, 1.824, 297.5], "gamma3": [1.1180, 0.6392, 0.6536, 1.1640], "epsilon3": [0.7289, 0.9074, 0.7655, 0.8711], "nr4": []} polt = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for benzene 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": CP2, "ref": "NBP", "Tmin": 278.7, "Tmax": 635.0, "Pmax": 78000.0, "rhomax": 11.45, "nr1": [-0.918572178424, 0.155357491575e1, -0.356149241161, 0.817273664265, -0.331303917534e1, 0.335336626528e1, -0.256976312022e1, 0.427304812515, 0.406483484297, -0.329744378187, 0.208907540720, 0.777471199254e-1, -0.202621443063, -0.148580350700e-1, 0.503167715817e-1, 0.293012717053e-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.918572178424, -0.155357491575e1, 0.356149241161, -0.447029533153e-1, 0.957712367542, -0.114688433057e1], "d2": [0, 0, 0, 2, 2, 2], "t2": [3, 4, 5, 3, 4, 5], "c2": [2]*6, "gamma2": [0.95481]*6} sun = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for benzene 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.314472, "cp": Fi1, "ref": "NBP", "Tmin": Tt, "Tmax": 620.0, "Pmax": 800000.0, "rhomax": 40., "nr1": [1.76284970, 1.02610647, -3.74263321, 9.57682041e-2, 2.59179321e-4, -1.03082188e-1], "d1": [1, 1, 1, 3, 7, 2], "t1": [1.5, 0.25, 1.25, 0.25, 0.875, 1.375], "nr2": [1.07359246e-1, -1.12562310e-1, 3.18737987e-1, -3.07549016e-2, -3.25082386e-1, 2.28099159e-2, -7.07431076e-2, -1.96809158e-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 = thol, polt, sun _surface = {"sigma": [0.07298, -0.0007802, -0.0001756], "exp": [1.232, 0.8635, 0.3065]} _vapor_Pressure = { "eq": 3, "n": [-0.71661e1, 0.21551e1, -0.20297e1, -0.40668e1, 0.38092], "t": [1.0, 1.5, 2.2, 4.8, 6.2]} _liquid_Density = { "eq": 1, "n": [0.18160e2, -0.56879e2, 0.87478e2, -0.64365e2, 0.18500e2], "t": [0.534, 0.686, 0.84, 1.0, 1.2]} _vapor_Density = { "eq": 2, "n": [-3.1147, -4.6689, -16.161, -146.5, 518.87, -827.72], "t": [0.419, 1.12, 2.8, 7.3, 10, 12]} visco0 = {"__name__": "Avgeri (2014)", "__doi__": { "autor": "Avgeri, S., Assael, M.J., Huber, M.L., Perkins, " "R.A", "title": "Reference Correlation of the Viscosity of Benzene " "from the Triple Point to 675K and up to 300MPa", "ref": "J. Phys. Chem. Ref. Data 43 (2014) 033103", "doi": "10.1063/1.4892935"}, "eq": 1, "omega": 1, "n_chapman": 0.021357, "ek": 412, "sigma": 0.54, "collision": [0.234018, -0.476136, 0, -0.015269], "Tref_virial": 412, "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], "special": "_mur"} def _mur(self, rho, T, fase): """Special form for residual term of Avgeri viscosity correlation, eq 8 in paper""" Tr = T/self.Tc rhor = rho/self.rhoc c = [-9.98945, 86.0626, 2.74872, 1.11130, -1, -134.133, -352.473, 6.60989, 88.4174] F = c[0]*rhor**2 + c[1]*rhor/(c[2]+c[3]*Tr+c[4]*rhor) + \ (c[5]*rhor+c[6]*rhor**2)/(c[7]+c[8]*rhor**2) mur = rhor**(2/3)*Tr**0.5*F return mur _viscosity = visco0, thermo0 = {"__name__": "Assael (2012)", "__doi__": { "autor": "Assael, M.J., Mihailidou, E., Huber, M.L., " "Perkins, R.A.", "title": "Reference Correlation of the Thermal " "Conductivity of Benzene from the Triple Point to " "725 K and up to 500 MPa", "ref": "J. Phys. Chem. Ref. Data 41(4) (2012) 043102", "doi": "10.1063/1.4755781"}, "eq": 1, "Toref": 562.02, "koref": 1e-3, "no_num": [101.404, -521.44, 868.266], "to_num": [0, 1, 2], "no_den": [1, 9.714, 1.467], "to_den": [0, 1, 2], "Tref_res": 562.02, "rhoref_res": 304.792, "kref_res": 1, "nr": [2.82489e-2, -7.73415e-2, 7.14001e-2, -2.36798e-2, 3.00875e-3, -1.19268e-2, 8.33389e-2, -8.98176e-2, 3.63025e-2, -4.90052e-3], "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": 2.16e-10, "gam0": 0.0569, "qd": 6.2e-10, "Tcref": 843} _thermal = thermo0,
class DMC(MEoS): """Multiparameter equation of state for dimethyl carbonate""" name = "dimethyl carbonate" CASNumber = "616-38-6" formula = "C3H6O3" synonym = "" rhoc = unidades.Density(360.3116) Tc = unidades.Temperature(557.) Pc = unidades.Pressure(4908.8, "kPa") M = 90.0779 # g/mol Tt = unidades.Temperature(277.06) Tb = unidades.Temperature(363.256) f_acent = 0.346 momentoDipolar = unidades.DipoleMoment(0.899, "Debye") # id=1798 Fi1 = { "ao_log": [1, 8.28421], "pow": [0, 1], "ao_pow": [4.9916462, -0.1709449], "ao_exp": [1.48525, 0.822585, 16.2453, 1.15925], "titao": [21 / Tc, 1340 / Tc, 1672 / Tc, 7395 / Tc], "ao_hyp": [], "hyp": [] } helmholtz1 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for DMC of Zhou et al. (2011).", "__doi__": { "autor": "Zhou, Y., Wu, J., and Lemmon, E.W.", "title": "Thermodynamic Properties of Dimethyl Carbonate", "ref": "J. Phys. Chem. Ref. Data, Vol. 40, No. 4 2011", "doi": "10.1063/1.3664084" }, "R": 8.314472, "cp": Fi1, "ref": { "Tref": 298.15, "Pref": 1.0, "ho": 26712.371, "so": 109.66202 }, "Tmin": Tt, "Tmax": 400.0, "Pmax": 60000.0, "rhomax": 12.107, "Pmin": 2.2495, "rhomin": 12.107, "nr1": [ 0.52683187e-3, 1.353396, -2.649283, -0.2785412, 0.1742554, 0.031606252 ], "d1": [5, 1, 1, 2, 3, 4], "t1": [1, 0.227, 1.05, 1.06, 0.5, 0.78], "nr2": [0.399866, 1.178144, -0.0235281, -1.015, -0.7880436, -0.12696], "d2": [1, 2, 7, 1, 2, 3], "t2": [1.3, 1.347, 0.706, 2, 2.5, 4.262], "c2": [1, 1, 1, 2, 2, 2], "gamma2": [1] * 6, "nr3": [1.2198, -0.4883, -0.0033293, -0.0035387, -0.51172, -0.16882], "d3": [1, 1, 2, 2, 3, 3], "t3": [1, 2.124, 0.4, 3.5, 0.5, 2.7], "alfa3": [0.9667, 1.5154, 1.0591, 1.6642, 12.4856, 0.9662], "beta3": [1.24, 0.821, 15.45, 2.21, 437., 0.743], "gamma3": [1.2827, 0.4317, 1.1217, 1.1871, 1.1243, 0.4203], "epsilon3": [0.6734, 0.9239, 0.8636, 1.0507, 0.8482, 0.7522], "nr4": [] } eq = helmholtz1, _vapor_Pressure = { "eq": 5, "ao": [-8.3197, 3.4260, -3.5905, -3.3194], "exp": [1.0, 1.5, 2.3, 4.7] } _liquid_Density = { "eq": 1, "ao": [1.1572, 4.969, -14.451, 27.569, -26.223, 10.526], "exp": [0.27, 0.77, 1.29, 1.85, 2.46, 3.16] } _vapor_Density = { "eq": 3, "ao": [-0.54715, -5.19277, -94.048, 327.21, -676.871, 716.072, -379.799], "exp": [0.197, 0.6, 2.86, 3.65, 4.5, 5.4, 6.4] } visco0 = { "eq": 1, "omega": 3, "__name__": "Zhou (2010)", "__doi__": { "autor": "Zhou, Y., Wu, J., and Lemmon, E.W.", "title": "Equations for the Thermophysical Properties of Dimethyl Carbonate", "ref": "AICHE Proceedings, 2009 Annual Meeting", "doi": "" }, "ek": 442.3, "sigma": 0.510747, "Tref": 557.376, "rhoref": 3.9749 * M, "n_chapman": 0.20555, "n_ideal": [], "t_ideal": [], "Tref_res": 557.376, "rhoref_res": 3.9749 * M, "etaref_res": 1, "n_poly": [5.07808, -0.056734, 0.00832177, 35.459838, 0.0513528], "t_poly": [-0.1, -3.0968, -2.8945, 0.0731, -3.9871], "d_poly": [4, 10, 12, 2, 0], "g_poly": [0, 0, 0, 0, 0, 0], "c_poly": [0, 1, 1, 2, 3] } _viscosity = visco0,
class R1216(MEoS): """Multiparameter equation of state for R1216""" name = "hexafluoropropene" CASNumber = "116-15-4" formula = "C3F6" synonym = "R1216" _refPropName = "R1216" _coolPropName = "" rhoc = unidades.Density(583.40757266496) Tc = unidades.Temperature(358.9) Pc = unidades.Pressure(3149.528, "kPa") M = 150.0225192 # g/mol Tt = unidades.Temperature(117.654) Tb = unidades.Temperature(242.81) f_acent = 0.333 momentoDipolar = unidades.DipoleMoment(1.088, "Debye") id = 669 CP1 = { "ao": 5.878676, "ao_exp": [9.351559, 9.192089, 7.983222], "exp": [561, 1486, 7595] } zhou = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for R1216 of Zhou (2010)", "__doi__": { "autor": "Zhou, Y., Lemmon, E.W.", "title": "Preliminary equation, 2010.", "ref": "", "doi": "" }, "R": 8.314472, "cp": CP1, "ref": "NBP", "Tmin": Tt, "Tmax": 400.0, "Pmax": 12000.0, "rhomax": 12.89, "nr1": [.37582356e-1, .14558246e1, -.2701615e1, -.3357347, .1885495], "d1": [4, 1, 1, 2, 3], "t1": [1.0, 0.3, 1.0, 1.35, 0.4], "nr2": [ -0.1689206, 0.1122147e1, -0.6405048, -0.25931535e-1, 0.42940852, -0.10163408e1, -0.43691328e-1 ], "d2": [3, 2, 2, 7, 1, 1, 1], "t2": [1., 1.68, 2.36, 0.615, 1.32, 2.12, 3.], "c2": [2, 1, 2, 1, 1, 2, 3], "gamma2": [1] * 7, "nr3": [0.12530663e1, -0.54254994, -0.15327764, -0.92102535e-2], "d3": [1, 1, 3, 3], "t3": [0.82, 2.85, 2.83, 1.67], "alfa3": [0.9665, 1.503, 0.97, 5.87], "beta3": [1.24, 0.776, 0.86, 478], "gamma3": [1.284, 0.42, 0.434, 1.074], "epsilon3": [0.67, 0.925, 0.75, 0.73], "nr4": [] } eq = zhou, _vapor_Pressure = { "eq": 3, "n": [-7.9011, 3.1506, -3.0852, -4.2112, -15.438], "t": [1.0, 1.5, 2., 4.5, 19.] } _liquid_Density = { "eq": 1, "n": [1.7159, 2.3953, -5.8035, 10.749, -10.537, 4.7535], "t": [0.31, 0.97, 1.7, 2.4, 3.2, 4.1] } _vapor_Density = { "eq": 2, "n": [-2.4969, -5.8935, -16.846, -55.082, -140.43], "t": [0.353, 1.05, 2.74, 6., 13.3] }
class Cyclohexane(MEoS): """Multiparameter equation of state for cyclohexane""" name = "cyclohexane" CASNumber = "110-82-7" formula = "cyclo(CH2)6" synonym = "" _refPropName = "CYCLOHEX" _coolPropName = "CycloHexane" rhoc = unidades.Density(271.33016352) Tc = unidades.Temperature(553.6) Pc = unidades.Pressure(4080.5, "kPa") M = 84.15948 # g/mol Tt = unidades.Temperature(279.47) Tb = unidades.Temperature(353.865) f_acent = 0.2096 momentoDipolar = unidades.DipoleMoment(0.3, "Debye") id = 38 _Tr = unidades.Temperature(526.231121) _rhor = unidades.Density(274.647526) _w = 0.221837522 Fi1 = {"ao_log": [1, 3], "pow": [0, 1], "ao_pow": [0.9891140602, 1.6359660572], "ao_exp": [0.83775, 16.036, 24.636, 7.1715], "titao": [773/Tc, 941/Tc, 2185/Tc, 4495/Tc]} CP1 = {"ao": 9.368327211, "an": [-56214088, 0.01526155409, -3.6352468e-6], "pow": [-3, 1, 2], "ao_exp": [23.766589], "exp": [2000]} zhou = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for cyclohexane of Zhou et " "al. (2014)", "__doi__": {"autor": "Zhou, Y., Liu, J., Penoncello, S.G., Lemmon, " "E.W.", "title": "An Equation of State for the Thermodynamic " "Properties of Cyclohexane", "ref": "J. Phys. Chem. Ref. Data 43 (2014) 043105", "doi": "10.1063/1.4900538"}, "R": 8.3144621, "cp": Fi1, "ref": "NBP", "Tmin": 279.86, "Tmax": 700.0, "Pmax": 250000.0, "rhomax": 10.3, "nr1": [0.05483581, 1.607734, -2.375928, -0.5137709, 0.1858417], "d1": [4, 1, 1, 2, 3], "t1": [1, 0.37, 0.79, 1.075, 0.37], "nr2": [-0.9007515, -0.5628776, 0.2903717, -0.3279141, -0.03177644], "d2": [1, 3, 2, 2, 7], "t2": [2.4, 2.5, 0.5, 3, 1.06], "c2": [2, 2, 1, 2, 1], "gamma2": [1]*5, "nr3": [0.8668676, -0.1962725, -0.1425992, 0.004197016, 0.1776584, -0.04433903, -0.03861246, 0.07399692, 0.02036006, 0.00272825], "d3": [1, 1, 3, 3, 2, 2, 3, 2, 3, 2], "t3": [1.6, 0.37, 1.33, 2.5, 0.9, 0.5, 0.73, 0.2, 1.5, 1.5], "alfa3": [0.99, 1.43, 0.97, 1.93, 0.92, 1.27, 0.87, 0.82, 1.4, 3], "beta3": [0.38, 4.2, 1.2, 0.9, 1.2, 2.6, 5.3, 4.4, 4.2, 25], "gamma3": [0.65, 0.63, 1.14, 0.09, 0.56, 0.4, 1.01, 0.45, 0.85, 0.86], "epsilon3": [0.73, 0.75, 0.48, 2.32, 0.2, 1.33, 0.68, 1.11, 1.47, .99]} penoncello = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for cyclohexane of " "Penoncello et al. (1995)", "__doi__": {"autor": "Penoncello, S.G., Jacobsen, R.T., Goodwin, " "A.R.H.", "title": "A Thermodynamic Property Formulation for " "Cyclohexane", "ref": "Int. J. Thermophys., 16(2) (1995) 519-531", "doi": "10.1007/BF01441918"}, "R": 8.31434, "cp": CP1, "ref": {"Tref": 279.47, "Pref": 101.325, "ho": 33884.8, "so": 96.612}, "Tt": 279.47, "Tc": 553.64, "rhoc": 3.244, "M": 84.1608, "Tmin": 279.47, "Tmax": 700.0, "Pmax": 80000.0, "rhomax": 9.77, "nr1": [0.8425412659, -0.3138388327e1, 0.1679072631e1, -0.153819249, 0.1984911143, -0.144532594, 0.3746346428e-3, 0.1861479616e-3, 0.1745721652e-3], "d1": [1, 1, 1, 2, 3, 3, 7, 6, 6], "t1": [0, 1.5, 2.5, 1.5, 1, 2.5, 2, 0.5, 3], "nr2": [-0.6427428062, 0.2280757615, -0.1868116802e1, -0.1028243711e1, 0.5821457418, -0.255891152, 0.1276844113e-1, -0.5158613166e-2, 0.6334794755e-1, -0.6014686589e-1, 0.4439056828, -0.6264920642, 0.2132589969e1, -0.3620300991e-2, 0.2534453992, 0.1669144715e-1, 0.3985052291e-2], "d2": [1, 1, 2, 3, 3, 5, 8, 10, 3, 4, 1, 1, 2, 2, 4, 4, 8], "t2": [5, 6, 5.5, 3, 7, 6, 6.5, 5.5, 11, 11, 0.5, 1, 4, 4, 1.5, 2, .5], "c2": [2, 2, 2, 2, 2, 2, 2, 2, 4, 4, 3, 3, 2, 6, 2, 4, 2], "gamma2": [1]*17} shortSpan = { "__type__": "Helmholtz", "__name__": "short Helmholtz equation of state for cyclohexane 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": CP1, "ref": {"Tref": 279.47, "Pref": 101.325, "ho": 33884.8, "so": 96.612}, "Tc": 553.60, "rhoc": 273.02/84.161, "M": 84.161, "Tmin": Tt, "Tmax": 750.0, "Pmax": 100000.0, "rhomax": 9.77, "nr1": [0.10232354e1, -0.29204964e1, 0.10736630e1, -0.19573985, 0.12228111, 0.28943321e-3], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875], "nr2": [0.27231767, -0.4483332e-1, -0.38253334, -0.89835333e-1, -0.24874965e-1, 0.10836132e-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} sun = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for cyclohexane 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.31434, "cp": CP1, "ref": {"Tref": 279.47, "Pref": 101.325, "ho": 33884.8, "so": 96.612}, "Tmin": Tt, "Tmax": 620.0, "Pmax": 800000.0, "rhomax": 40., "nr1": [1.27436292, 1.15372124, -3.86726473, 8.84627298e-2, 2.76478090e-4, 7.26682313e-2], "d1": [1, 1, 1, 3, 7, 2], "t1": [1.5, 0.25, 1.25, 0.25, 0.875, 1.375], "nr2": [7.10849914e-2, 4.46376742e-1, 7.64476190e-1, -4.23520282e-2, -0.396468623, -1.41250071e-2, -1.08371284e-1, -2.50082884e-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 = zhou, penoncello, shortSpan, sun _surface = {"sigma": [0.06485], "exp": [1.263]} _melting = { "eq": 1, "__doi__": {"autor": "Wisotzki, K.D., Wǘrflinger, A.", "title": "PVT Data for Liquid and Solid Cyclohexane, " "Cyclohexanone and Cyclopentanol up to 3000 bar", "ref": "J. Phis. Chem. Solids 43(1) (1982) 13-20", "doi": "10.1016/0022-3697(82)90167-6"}, "Tmin": Tt, "Tmax": 700.0, "Tref": 279.7, "Pref": 1, "a2": [3834e5], "exp2": [1.41]} _vapor_Pressure = { "eq": 3, "n": [-7.0342, 1.7311, -1.7572, -3.3406], "t": [1., 1.5, 2.3, 4.6]} _liquid_Density = { "eq": 1, "n": [5.5081, -14.486, 38.241, -64.589, 57.919, -20.55], "t": [0.51, 0.94, 1.4, 1.9, 2.4, 3.0]} _vapor_Density = { "eq": 2, "n": [-3.69006, -41.4239, 220.914, -443.72, 491.49, -296.373], "t": [0.446, 1.98, 2.75, 3.3, 4.1, 4.8]} visco0 = {"__name__": "Tariq (2014)", "__doi__": { "autor": "Tariq, U., Jusoh, A.R.B., Riesco, N., Vesovic, V.", "title": "Reference Correlation of the Viscosity of " "Cyclohexane from the Triple Point to 700K and up " "to 110 MPa", "ref": "J. Phys. Chem. Ref. Data 43(3) (2014) 033101", "doi": "10.1063/1.4891103"}, "eq": 1, "omega": 3, "collision": [-1.5093, 364.87, -39537], "sigma": 1, "n_chapman": 0.19592/M**0.5, "Tref_res": 553.6, "rhoref_res": 3.224*M, "nr": [335.234, 7.8494803, -687.3976, 362.0868, -10.4793856, 2.5521774, 17.2734993, -5.9372242, -10.6186149, 4.3982781, 2.8894928, -1.3468174, -0.2938491, 0.1487134], "dr": [2.2, 2.5, 2.5, 2.8, 10, 10, 11, 11, 12, 12, 13, 13, 14, 14], "tr": [1, 0, 1, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1], "gr": [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], "cr": [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0], "special": "_vir"} def _vir(self, rho, T, fase): # The initial density dependence has a different expresion, without muo # and other normal method calculation so hardcoded here muB = 0 if rho: for i, n in enumerate([5.09643, -3387.21, 337477]): muB += n/T**i return muB*rho/self.M _viscosity = visco0,
class mXylene(MEoS): """Multiparameter equation of state for m-xylene""" name = "m-xylene" CASNumber = "108-38-3" formula = "C8H10" synonym = "1,3-dimethylbenzene" rhoc = unidades.Density(282.929725) Tc = unidades.Temperature(616.89) Pc = unidades.Pressure(3534.6, "kPa") M = 106.165 # g/mol Tt = unidades.Temperature(225.3) Tb = unidades.Temperature(412.214) f_acent = 0.326 momentoDipolar = unidades.DipoleMoment(0.3, "Debye") id = 43 Fi1 = { "ao_log": [1, 1.169909], "pow": [0, 1], "ao_pow": [12.652887, -0.45975624], "ao_exp": [4.44312, 2.862794, 24.83298, 16.26077], "titao": [160 / Tc, 190 / Tc, 1333 / Tc, 3496 / Tc] } helmholtz1 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for m-xylene of Zhou et al. (2012).", "__doi__": { "autor": "Zhou, Y., Lemmon, E.W., and Wu, J.", "title": "Thermodynamic Properties of o-Xylene, m-Xylene, p-Xylene, and Ethylbenzene", "ref": "J. Phys. Chem. Ref. Data 41, 023103 (2012).", "doi": "10.1063/1.3703506" }, "R": 8.314472, "cp": Fi1, "ref": "OTO", "Tmin": Tt, "Tmax": 700.0, "Pmax": 200000.0, "rhomax": 8.677, "Pmin": 0.003123, "rhomin": 8.677, "nr1": [ 0.000012791017, 0.041063111, 1.505996, -2.3095875, -0.46969, 0.171031 ], "d1": [8, 4, 1, 1, 2, 3], "t1": [1.0, 0.91, 0.231, 0.772, 1.205, 0.323], "nr2": [-1.001728, -0.3945766, 0.6970578, -0.3002876, -0.024311], "d2": [1, 3, 2, 2, 7], "t2": [2.7, 3.11, 0.768, 4.1, 0.818], "c2": [2, 2, 1, 2, 1], "gamma2": [1] * 5, "nr3": [0.815488, -0.330647, -0.123393, -0.54661], "d3": [1, 1, 3, 3], "t3": [2.0, 2.9, 3.83, 0.5], "alfa3": [1.0244, 1.3788, 0.9806, 6.3563], "beta3": [1.66, 1.9354, 1.0323, 78], "gamma3": [1.1013, 0.6515, 0.4975, 1.26], "epsilon3": [0.713, 0.9169, 0.6897, 0.7245] } eq = helmholtz1, _surface = {"sigma": [0.0661], "exp": [1.29]} _vapor_Pressure = { "eq": 5, "ao": [-7.5635, 1.2857, -3.2346, -1.9018], "exp": [1.0, 1.5, 3.1, 5.6] } _liquid_Density = { "eq": 1, "ao": [0.43346, 3.8716, -3.0144, 1.619], "exp": [0.16, 0.6, 1.0, 1.5] } _vapor_Density = { "eq": 3, "ao": [-1.1597, -6.0358, -16.712, -45.482, -98.418], "exp": [0.26, 0.78, 2.6, 5.7, 11.7] }
class D2O(MEoS): """Multiparameter equation of state for heavy water""" name = "heavy water" CASNumber = "7789-20-0" formula = "D2O" synonym = "deuterium oxide" Tc = unidades.Temperature(643.847) rhoc = unidades.Density(358) Pc = unidades.Pressure(21671.0, "kPa") M = 20.027508 # g/mol Tt = unidades.Temperature(276.97) Tb = unidades.Temperature(374.563) f_acent = 0.364 momentoDipolar = unidades.DipoleMoment(1.9, "Debye") Fi1 = { "ao_log": [0.5399322597e-2, 0], "pow": [0, 1, 2, 3, 4, 5], "ao_pow": [ 0.3087155964e2, -.3827264031e2, 0.4424799189, -.1256336874e1, 0.2843343470, -.2401555088e-1 ], "tau*logtau": -.1288399716e2, "tau*logdelta": 0.4415884023e1, "ao_exp": [], "titao": [], "ao_hyp": [], "hyp": [] } CP1 = { "ao": 0.39176485e1, "an": [ -0.31123915e-3, 0.41173363e-5, -0.28943955e-8, 0.63278791e-12, 0.78728740 ], "pow": [1.00, 2.00, 3.00, 4.00, -0.99], "ao_exp": [], "exp": [], "ao_hyp": [], "hyp": [] } helmholtz1 = { "__type__": "Helmholtz", "__name__": u"Helmholtz equation of state for heavy water of Hill et al. (1982).", "__doi__": { "autor": "Hill, P.G., MacMillan, R.D.C., and Lee, V.", "title": "A Fundamental Equation of State for Heavy Water", "ref": "J. Phys. Chem. Ref. Data 11, 1 (1982)", "doi": "10.1063/1.555661" }, "__test__": # Pag 17 of IAPWS 2007 update paper """ >>> st=D2O(T=0.5*D2O.Tc, rho=0.0002*D2O.rhoc) >>> print "%0.6f %0.7f %0.4f" % (st.a*D2O.rhoc/D2O.Pc, st.Pr, st.cv.kJkgK*D2O.rhoc*D2O.Tc/D2O.Pc.kPa) -2.644979 0.0004402 14.2768 >>> st=D2O(T=0.5*D2O.Tc, rho=3.18*D2O.rhoc) >>> print "%0.6f %0.7f %0.4f" % (st.a*D2O.rhoc/D2O.Pc, st.Pr, st.cv.kJkgK*D2O.rhoc*D2O.Tc/D2O.Pc.kPa) -0.217388 4.3549719 41.4463 >>> st=D2O(T=0.75*D2O.Tc, rho=0.0295*D2O.rhoc) >>> print "%0.6f %0.7f %0.4f" % (st.a*D2O.rhoc/D2O.Pc, st.Pr, st.cv.kJkgK*D2O.rhoc*D2O.Tc/D2O.Pc.kPa) -7.272543 0.0870308 20.1586 >>> st=D2O(T=0.75*D2O.Tc, rho=2.83*D2O.rhoc) >>> print "%0.6f %0.7f %0.4f" % (st.a*D2O.rhoc/D2O.Pc, st.Pr, st.cv.kJkgK*D2O.rhoc*D2O.Tc/D2O.Pc.kPa) -4.292707 4.4752958 33.4367 >>> st=D2O(T=D2O.Tc, rho=0.3*D2O.rhoc) >>> print "%0.6f %0.7f %0.4f" % (st.a*D2O.rhoc/D2O.Pc, st.Pr, st.cv.kJkgK*D2O.rhoc*D2O.Tc/D2O.Pc.kPa) -15.163326 0.8014044 30.8587 >>> st=D2O(T=D2O.Tc, rho=1.55*D2O.rhoc) >>> print "%0.6f %0.7f %0.4f" % (st.a*D2O.rhoc/D2O.Pc, st.Pr, st.cv.kJkgK*D2O.rhoc*D2O.Tc/D2O.Pc.kPa) -12.643811 1.0976283 33.0103 >>> st=D2O(T=1.2*D2O.Tc, rho=0.4*D2O.rhoc) >>> print "%0.6f %0.7f %0.4f" % (st.a*D2O.rhoc/D2O.Pc, st.Pr, st.cv.kJkgK*D2O.rhoc*D2O.Tc/D2O.Pc.kPa) -25.471535 1.4990994 23.6594 >>> st=D2O(T=1.2*D2O.Tc, rho=1.61*D2O.rhoc) >>> print "%0.6f %0.7f %0.4f" % (st.a*D2O.rhoc/D2O.Pc, st.Pr, st.cv.kJkgK*D2O.rhoc*D2O.Tc/D2O.Pc.kPa) -21.278164 4.5643798 25.4800 """, "R": 8.3143565, "rhoref": 17.875414 * M, "cp": Fi1, "ref": { "Tref": 276.95, "Pref": 0.660096, "ho": 0.598, "so": 0 }, "Tmin": Tt, "Tmax": 800.0, "Pmax": 100000.0, "rhomax": 65., "Pmin": 0.66103, "rhomin": 55.198, "nr1": [ -0.384820628204e3, 0.108213047259e4, -0.110768260635e4, 0.164668954246e4, -0.137959852228e4, 0.598964185629e3, -0.100451752702e3, 0.419192736351e3, -0.107279987867e4, 0.653852283544e3, -0.984305985655e3, 0.845444459339e3, -0.376799930490e3, 0.644512590492e2, -0.214911115714e3, 0.531113962967e3, -0.135454224420e3, 0.202814416558e3, -0.178293865031e3, 0.818739394970e2, -0.143312594493e2, 0.651202383207e2, -0.171227351208e3, 0.100859921516e2, -0.144684680657e2, 0.128871134847e2, -0.610605957134e1, 0.109663804408e1, -0.115734899702e2, 0.374970075409e2, 0.897967147669, -0.527005883203e1, 0.438084681795e-1, 0.406772082680, -0.965258571044e-2, -0.119044600379e-1 ], "d1": [ 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 5, 5, 6, 6, 7, 7, 8, 8 ], "t1": [ 0, 1, 2, 3, 4, 5, 6, 0, 1, 2, 3, 4, 5, 6, 0, 1, 2, 3, 4, 5, 6, 0, 1, 2, 3, 4, 5, 6, 0, 1, 0, 1, 0, 1, 0, 1 ], "nr2": [ 0.382589102341e3, -0.106406466204e4, 0.105544952919e4, -0.157579942855e4, 0.132703387531e4, -0.579348879870e3, 0.974163902526e2, 0.286799294226e3, -0.127543020847e4, 0.275802674911e4, -0.381284331492e4, 0.293755152012e4, -0.117858249946e4, 0.186261198012e3 ], "c2": [1] * 14, "d2": [1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2], "t2": [0, 1, 2, 3, 4, 5, 6, 0, 1, 2, 3, 4, 5, 6], "gamma2": [1.5394] * 14 } eq = helmholtz1, _surface = {"sigma": [0.238, -0.152082], "exp": [1.25, 2.25]} _vapor_Pressure = { "eq": 5, "ao": [-0.80236e1, 0.23957e1, -0.42639e2, 0.99569e2, -0.62135e2], "exp": [1.0, 1.5, 2.75, 3.0, 3.2] } _liquid_Density = { "eq": 1, "ao": [0.26406e1, 0.97090e1, -0.18058e2, 0.87202e1, -0.74487e1], "exp": [0.3678, 1.9, 2.2, 2.63, 7.3] } _vapor_Density = { "eq": 3, "ao": [-0.37651e1, -0.38673e2, 0.73024e2, -0.13251e3, 0.75235e2, -0.70412e2], "exp": [0.409, 1.766, 2.24, 3.04, 3.42, 6.9] } visco0 = { "eq": 0, "method": "_visco0", "__name__": "IAPWS (2007)", "__doi__": { "autor": "J. Kestin, J. V. Sengers, B. Kamgar‐Parsi and J. M. H. Levelt Sengers", "title": "Thermophysical Properties of Fluid D2O", "ref": "J. Phys. Chem. Ref. Data 13, 601 (1984)", "doi": "10.1063/1.555714" }, "__test__": # Pag 17 of IAPWS 2007 update paper """ >>> st=D2O(T=0.431*D2O.Tc, rho=3.09*D2O.rhoc) >>> print "%0.10f" % (st.mu.muPas/55.2651) 36.9123166244 >>> st=D2O(T=0.431*D2O.Tc, rho=3.23*D2O.rhoc) >>> print "%0.10f" % (st.mu.muPas/55.2651) 34.1531546602 >>> st=D2O(T=0.6*D2O.Tc, rho=2.95*D2O.rhoc) >>> print "%0.10f" % (st.mu.muPas/55.2651) 5.2437249935 >>> st=D2O(T=D2O.Tc, rho=0.7*D2O.rhoc) >>> print "%0.10f" % (st.mu.muPas/55.2651) 0.5528693914 >>> st=D2O(T=0.9*D2O.Tc, rho=0.08*D2O.rhoc) >>> print "%0.10f" % (st.mu.muPas/55.2651) 0.3685472578 >>> st=D2O(T=1.1*D2O.Tc, rho=0.98*D2O.rhoc) >>> print "%0.10f" % (st.mu.muPas/55.2651) 0.7816387903 >>> st=D2O(T=1.2*D2O.Tc, rho=0.8*D2O.rhoc) >>> print "%0.10f" % (st.mu.muPas/55.2651) 0.7651099154 """, } def _visco0(self, rho, T, fase=None): Tr = T / 643.847 rhor = rho / 358.0 no = [1.0, 0.940695, 0.578377, -0.202044] fi0 = Tr**0.5 / sum([n / Tr**i for i, n in enumerate(no)]) Li = [ 0, 1, 2, 3, 4, 5, 0, 1, 2, 3, 0, 1, 2, 5, 0, 1, 2, 3, 0, 1, 3, 5, 0, 1, 5, 3 ] Lj = [ 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 5, 5, 5, 6 ] Lij = [ 0.4864192, -0.2448372, -0.8702035, 0.8716056, -1.051126, 0.3458395, 0.3509007, 1.315436, 1.297752, 1.353448, -0.2847572, -1.037026, -1.287846, -0.02148229, 0.07013759, 0.4660127, 0.2292075, -0.4857462, 0.01641220, -0.02884911, 0.1607171, -0.009603846, -0.01163815, -0.008239587, 0.004559914, -0.003886659 ] array = [ lij * (1. / Tr - 1)**i * (rhor - 1)**j for i, j, lij in zip(Li, Lj, Lij) ] fi1 = exp(rhor * sum(array)) return unidades.Viscosity(55.2651 * fi0 * fi1, "muPas") _viscosity = visco0, thermo0 = { "eq": 0, "method": "_thermo0", "__name__": "IAPWS (1994)", "__doi__": { "autor": "J. Kestin, J. V. Sengers, B. Kamgar‐Parsi and J. M. H. Levelt Sengers", "title": "Thermophysical Properties of Fluid D2O", "ref": "J. Phys. Chem. Ref. Data 13, 601 (1984)", "doi": "10.1063/1.555714" }, "__test__": # Pag 17 of IAPWS 2007 update paper """ >>> st=D2O(T=0.431*D2O.Tc, rho=3.09*D2O.rhoc) >>> print "%0.9f" % (st.k/0.742128) 762.915707396 >>> st=D2O(T=0.431*D2O.Tc, rho=3.23*D2O.rhoc) >>> print "%0.9f" % (st.k/0.742128) 833.912049618 >>> st=D2O(T=0.6*D2O.Tc, rho=2.95*D2O.rhoc) >>> print "%0.9f" % (st.k/0.742128) 861.240794445 >>> st=D2O(T=D2O.Tc, rho=0.7*D2O.rhoc) >>> print "%0.9f" % (st.k/0.742128) 469.015122112 >>> st=D2O(T=0.9*D2O.Tc, rho=0.08*D2O.rhoc) >>> print "%0.9f" % (st.k/0.742128) 74.522283066 >>> st=D2O(T=1.1*D2O.Tc, rho=0.98*D2O.rhoc) >>> print "%0.9f" % (st.k/0.742128) 326.652382218 >>> st=D2O(T=1.2*D2O.Tc, rho=0.8*D2O.rhoc) >>> print "%0.9f" % (st.k/0.742128) 259.605241187 """, } def _thermo0(self, rho, T, fase=None): rhor = rho / 358 Tr = T / 643.847 tau = Tr / (abs(Tr - 1.1) + 1.1) no = [1.0, 37.3223, 22.5485, 13.0465, 0.0, -2.60735] Lo = sum([Li * Tr**i for i, Li in enumerate(no)]) nr = [483.656, -191.039, 73.0358, -7.57467] Lr = -167.31 * (1 - exp(-2.506 * rhor)) + sum( [Li * rhor**(i + 1) for i, Li in enumerate(nr)]) f1 = exp(0.144847 * Tr - 5.64493 * Tr**2) f2 = exp(-2.8 * (rhor - 1)**2) - 0.080738543 * exp(-17.943 * (rhor - 0.125698)**2) f3 = 1 + exp(60 * (tau - 1) + 20) f4 = 1 + exp(100 * (tau - 1) + 15) Lc = 35429.6 * f1 * f2 * (1 + f2**2 * (5e9 * f1**4 / f3 + 3.5 * f2 / f4)) Ll = -741.112 * f1**1.2 * (1 - exp(-(rhor / 2.5)**10)) return unidades.ThermalConductivity(0.742128e-3 * (Lo + Lr + Lc + Ll)) _thermal = thermo0,