class iButene(MEoS): """Multiparameter equation of state for isobutene""" name = "isobutene" CASNumber = "115-11-7" formula = "CH2=C(CH3)2" synonym = "" rhoc = unidades.Density(233.9633544) Tc = unidades.Temperature(418.09) Pc = unidades.Pressure(4009.8, "kPa") M = 56.10632 # g/mol Tt = unidades.Temperature(132.4) Tb = unidades.Temperature(266.15) f_acent = 0.193 momentoDipolar = unidades.DipoleMoment(0.5, "Debye") id = 27 Fi1 = { "ao_log": [1, 3.], "pow": [0, 1], "ao_pow": [-0.12737888, 2.3125128], "ao_exp": [4.8924, 7.832, 7.2867, 8.7293], "titao": [399 / Tc, 1270 / Tc, 2005 / Tc, 4017 / Tc] } helmholtz1 = { "__type__": "Helmholtz", "__name__": "short Helmholtz equation of state for isobutene of Lemmon and Ihmels (2005)", "__doi__": { "autor": "Lemmon, E.W., Ihmels, E.C.", "title": "Thermodynamic properties of the butenes: Part II. Short fundamental equations of state", "ref": "Fluid Phase Equilibria 228 – 229 (2004), 173 – 187.", "doi": "10.1016/j.fluid.2004.09.004" }, "__test__": """ >>> st=iButene(T=350, rho=0) >>> print "%0.0f %0.1f %0.1f %0.5g %0.5g %0.5g %0.5g" % (st.T, st.rhoM, st.P.MPa, st.hM.kJkmol, st.cvM.kJkmolK, st.cpM.kJkmolK, st.w) 350 0.0 0.0 29966 92.121 100.44 237.8 >>> st=iButene(T=350, rho=0.3*56.10632) >>> print "%0.0f %0.1f %.5g %.5g %0.5g %0.5g %0.5g %0.5g" % (st.T, st.rhoM, st.P.MPa, st.hM.kJkmol, st.sM.kJkmolK, st.cvM.kJkmolK, st.cpM.kJkmolK, st.w) 350 0.3 0.75754 28666 88.966 96.794 112.57 211.02 >>> st=iButene(T=350, rho=10*56.10632) >>> print "%0.0f %0.1f %.5g %.5g %0.5g %0.5g %0.5g %0.5g" % (st.T, st.rhoM, st.P.MPa, st.hM.kJkmol, st.sM.kJkmolK, st.cvM.kJkmolK, st.cpM.kJkmolK, st.w) 350 10.0 17.776 11782 32.951 101.72 139.45 838.25 >>> st=iButene(T=440, rho=4*56.10632) >>> print "%0.0f %0.1f %.5g %.5g %0.5g %0.5g %0.5g %0.5g" % (st.T, st.rhoM, st.P.MPa, st.hM.kJkmol, st.sM.kJkmolK, st.cvM.kJkmolK, st.cpM.kJkmolK, st.w) 440 4.0 5.4086 30169 82.345 127.28 407 151.13 """, # Table 9, Pag 186 "R": 8.314472, "cp": Fi1, "ref": "NBP", "Tmin": Tt, "Tmax": 550.0, "Pmax": 50000.0, "rhomax": 13.67, "Pmin": 0.00068, "rhomin": 13.67, "nr1": [0.77111, -2.7971, 1.0118, 0.02073, 0.085086, 0.00021968], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.12, 1.3, 1.74, 2.1, 0.28, 0.69], "nr2": [0.20633, -0.078843, -0.23726, -0.080211, -0.027001, 0.013072], "d2": [2, 5, 1, 4, 3, 4], "t2": [0.75, 2., 4.4, 4.7, 15., 14.], "c2": [1, 1, 2, 2, 3, 3], "gamma2": [1] * 6 } eq = helmholtz1, _surface = {"sigma": [0.0545], "exp": [1.23]} _vapor_Pressure = { "eq": 5, "ao": [-0.68973e1, 0.12475e1, -0.25441e1, -0.29282e1, 0.15778e1], "exp": [1., 1.5, 3.16, 6.2, 7.0] } _liquid_Density = { "eq": 1, "ao": [0.62591e2, -0.20805e3, 0.33243e3, -0.29555e3, 0.11148e3], "exp": [0.65, 0.8, 0.98, 1.16, 1.3] } _vapor_Density = { "eq": 3, "ao": [ -0.31841e1, -0.64014e1, -0.93817e1, -0.11160e2, -0.52298e2, -0.12195e3 ], "exp": [0.431, 1.29, 3.3, 3.54, 7.3, 15.8] }
class C1Linoleate(MEoS): """Multiparameter equation of state for methyl linoleate""" name = "methyl linoleate" CASNumber = "112-63-0" formula = "C19H34O2" synonym = "" rhoc = unidades.Density(238.051213304) Tc = unidades.Temperature(799.0) Pc = unidades.Pressure(1341.0, "kPa") M = 294.47206 # g/mol Tt = unidades.Temperature(238.1) Tb = unidades.Temperature(628.84) f_acent = 0.805 momentoDipolar = unidades.DipoleMoment(1.79, "Debye") id = 39 CP1 = { "ao": 0.0, "an": [190.986], "pow": [0.020213], "ao_exp": [437.371, 287.222, 321.956], "exp": [3052.11, 746.631, 1624.33], "ao_hyp": [], "hyp": [] } helmholtz1 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for methyl linoleate 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, 2009, 23 (7), pp 3790–3797", "doi": "10.1021/ef900159g" }, "R": 8.314472, "cp": CP1, "ref": "NBP", "Tmin": Tt, "Tmax": 1000.0, "Pmax": 50000.0, "rhomax": 3.16, "Pmin": 0.7e-14, "rhomin": 3.16, "nr1": [0.3183187e-1, 0.1927286e1, -0.3685053e1, 0.8449312e-1], "d1": [4, 1, 1, 3], "t1": [1, 0.2, 1.2, 1.0], "nr2": [-0.9766643, -0.4323178, 0.2000470e1, -0.1752030e1, -0.1726895e-1], "d2": [1, 3, 2, 2, 7], "t2": [2.2, 2.5, 1.8, 1.92, 1.47], "c2": [2, 2, 1, 2, 1], "gamma2": [1] * 5, "nr3": [0.2116515e1, -0.7884271, -0.3811699], "d3": [1, 1, 3], "t3": [1.7, 2.3, 2.1], "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 = helmholtz1, _vapor_Pressure = { "eq": 5, "ao": [-0.10946e2, 0.48849e1, -0.46773e1, -0.80201e1, -0.89572e1], "exp": [1.0, 1.5, 2.22, 3.6, 8.0] } _liquid_Density = { "eq": 1, "ao": [0.22705e3, -0.66763e3, 0.72323e3, -0.49244e3, 0.21391e3], "exp": [0.83, 0.98, 1.17, 1.5, 1.7] } _vapor_Density = { "eq": 3, "ao": [-0.85880e1, 0.14766e2, -0.24195e2, -0.37474e3, 0.32689e3, -0.19125e3], "exp": [0.568, 1.08, 1.4, 4.8, 5.0, 9.0] } thermo0 = { "eq": 1, "__name__": "Perkins (2010)", "__doi__": { "autor": "Perkins, R.A. and Huber, M.L.", "title": "Measurement and Correlation of the Thermal Conductivities of Biodiesel Constituent Fluids: Methyl Oleate and Methyl Linoleate", "ref": "Energy Fuels, 2011, 25 (5), pp 2383–2388", "doi": "10.1021/ef200417x" }, "__test__": """ >>> st=C1Linoleate(T=450, P=1e2) >>> print "%0.0f %0.4f %0.6g %0.6g" % (st.T, st.P.MPa, st.rho, st.k.WmK) 450 0.0001 0.00787223 0.0122743 >>> st=C1Linoleate(T=450, P=1e6) >>> print "%0.0f %0.0f %0.6g %0.6g" % (st.T, st.P.MPa, st.rho, st.k.WmK) 450 1 778.176 0.122742 >>> st=C1Linoleate(T=450, P=2e7) >>> print "%0.0f %0.0f %0.6g %0.6g" % (st.T, st.P.MPa, st.rho, st.k.WmK) 450 20 799.16 0.131867 """, # Table 3, Pag 2386 "Tref": 799.0, "kref": 1, "no": [-0.10904200e-3, 0.24054300e-2, 0.40736400e-1, -0.10592800e-1], "co": [0, 1, 2, 3], "Trefb": 799.0, "rhorefb": 0.8084 * M, "krefb": 1., "nb": [ -0.713126e-1, 0.466421e-1, -0.557406e-2, 0.0, 0.0, 0.989415e-1, -0.65785e-1, 0.128922e-1, 0.0, 0.0 ], "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.239, "R0": 1.03, "Xio": 0.194e-9, "gam0": 0.0496, "qd": 8.75e-10, "Tcref": 1198.5 } thermo1 = { "eq": 5, "omega": 3, "__name__": "Chung (1988)", "__doi__": { "autor": "T-H. Chung, Ajlan, M., Lee, L.L. and Starling, K.E.", "title": "Generalized Multiparameter Correlation for Nonpolar and Polar Fluid Transport Properties", "ref": "Ind. Eng. Chem. Res., 1988, 27 (4), pp 671–679", "doi": "10.1021/ie00076a024" }, "w": 0.805, "mur": 0.0, "k": 0.0 } _thermal = thermo0, thermo1
class R13(MEoS): """Multiparameter equation of state for R13""" name = "chlorotrifluoromethane" CASNumber = "75-72-9" formula = "CClF3" synonym = "R13" _refPropName = "R13" _coolPropName = "R13" rhoc = unidades.Density(582.88) Tc = unidades.Temperature(302.0) Pc = unidades.Pressure(3879.0, "kPa") M = 104.459 # g/mol Tt = unidades.Temperature(92.0) Tb = unidades.Temperature(191.67) f_acent = 0.1723 momentoDipolar = unidades.DipoleMoment(0.51, "Debye") id = 215 # Cp⁰ coefficient for Tr term, so divide by Tc CP1 = { "ao": 1.86012334, "an": [8.07314520 / Tc, -1.87713639 / Tc**2, 3.17242858e-2 / Tc**3], "pow": [1, 2, 3], "ao_exp": [], "exp": [] } CP2 = { "ao": 2.4766458, "an": [0.018074269, 2.1945535e-5, -8.5810657e-8, 6.3199171e-11], "pow": [1, 2, 3, 4], "ao_exp": [], "exp": [] } magee = { "__type__": "MBWR", "__name__": "MBWR equation of state for R-13 of Magee et al. (2000)", "__doi__": { "autor": "Magee, J.W., Outcalt, S.L., Ely, J.F.", "title": "Molar Heat Capacity Cv, Vapor Pressure, and " "(p, ρ, T) Measurements from 92 to 350 K at " "Pressures to 35 MPa and a New Equation of State " "for Chlorotrifluoromethane (R13)", "ref": "Int. J. Thermophys., 21(5):1097-1121, 2000.", "doi": "10.1023/A:1026446004383" }, "R": 8.314471, "Tc": 302, "Pc": 3879., "rhoc": 5.58, "M": 104.459, "cp": CP1, "ref": "IIR", "Tmin": Tt, "Tmax": 403.0, "Pmax": 35000.0, "rhomax": 17.85, "b": [ None, 0.427710490378e-2, 0.106603397093e1, -0.383065097813e2, 0.661580211522e4, -0.800160780370e6, -0.406405755462e-2, 0.561380767634e1, -0.247694806929e4, -0.639834580892e5, 0.198818486764e-3, -0.206916891385, 0.749317872337e2, -0.431471653965e-2, 0.181741326553e1, -0.206066849491e2, -0.136681208829, 0.260496240940e-2, 0.287244312242, -0.105459756169e-1, 0.582404815872e6, -0.455721947029e8, 0.114174177352e5, 0.265590236008e6, 0.135249873550e3, 0.128289104267e4, 0.800900540368, -0.703307137789e4, 0.235567665577e-2, 0.131830636112e1, -0.115187941781e-4, 0.564530387616e-2, 0.336242130107 ] } platzer = { "__type__": "Helmholtz", "__name__": "Bender equation of state for R-13 of Platzer (1990)", "__doi__": { "autor": "Platzer, B., Polt, A., Maurer, G.", "title": "Thermophysical Properties of refrigerants", "ref": "Berlin: Springer-Verlag, 1990.", "doi": "" }, "R": 8.31451, "cp": CP2, "ref": "IIR", "Tmin": Tt, "Tmax": 450.0, "Pmax": 50000.0, "rhomax": 17.699806, "nr1": [ -0.628346559920, 0.792797111341, -0.134038992692, 0.761143010172, -0.194465098795e1, 0.940938700406, -0.108107050239e1, 0.117501564976, 0.228305167217, -0.403338888789, 0.375585713420, -0.617543677315e-1, 0.170326226881, 0.536612457231e-1, -0.151603010301, 0.252033265074e-1 ], "d1": [0, 0, 0, 1, 1, 1, 1, 1, 2, 2, 2, 3, 3, 4, 4, 5], "t1": [3, 4, 5, 0, 1, 2, 3, 4, 0, 1, 2, 0, 1, 0, 1, 1], "nr2": [ 0.628346559920, -0.792797111341, 0.134038992692, -0.399863840975e-1, 0.436410910529, -0.448724904991 ], "d2": [0, 0, 0, 2, 2, 2], "t2": [3, 4, 5, 3, 4, 5], "c2": [2] * 6, "gamma2": [0.98230055] * 6 } eq = magee, platzer _surface = {"sigma": [0.05045], "exp": [1.269]} _vapor_Pressure = { "eq": 3, "n": [-0.69311e1, 0.18281e1, -0.21901e1, -0.38177e1, 0.20803e1], "t": [1.0, 1.5, 2.5, 6.0, 8.0] } _liquid_Density = { "eq": 1, "n": [0.95469e1, -0.24017e2, 0.33365e2, -0.26837e2, 0.10638e2], "t": [0.51, 0.72, 0.94, 1.2, 1.4] } _vapor_Density = { "eq": 2, "n": [-0.31949e1, -0.73425e1, -0.21966e2, -0.51459e2, -0.85359e2], "t": [0.414, 1.41, 3.7, 7.7, 15.0] } trnECS = { "__name__": "Huber (2003)", "__doi__": { "autor": "Huber, M.L., Laesecke, A., Perkins, R.A.", "title": "Model for the Viscosity and Thermal Conductivity " "of Refrigerants, Including a New Correlation for " "the Viscosity of R134a", "ref": "Ind. Eng. Chem. Res., 42(13) (2003) 3163-3178", "doi": "10.1021/ie0300880" }, "eq": "ecs", "ref": C3, "visco": "visco1", "thermo": "thermo0", "ek": 204.0, "sigma": 0.4971, "omega": 5, "psi": [0.976177, 1.48047e-2], "psi_d": [0, 1], "fint": [1.07447e-3, 6.42373e-7], "fint_t": [0, 1], "chi": [1.1394, -3.65562e-2], "chi_d": [0, 1], "critical": 3, "gnu": 0.63, "gamma": 1.239, "R0": 1.03, "Xio": 0.194e-9, "gam0": 0.0496, "qd": 3.49636e-10, "Tcref": 453.0 } _viscosity = trnECS, _thermal = trnECS,
class Methanol(MEoS): """Multiparameter equation of state for methanol""" name = "methanol" CASNumber = "67-56-1" formula = "CH3OH" synonym = "" rhoc = unidades.Density(275.5626) Tc = unidades.Temperature(512.6) Pc = unidades.Pressure(8103.5, "kPa") M = 32.04216 # g/mol Tt = unidades.Temperature(175.61) Tb = unidades.Temperature(337.632) f_acent = 0.5625 momentoDipolar = unidades.DipoleMoment(1.7, "Debye") id = 117 CP1 = { "ao": 3.9007912, "an": [], "pow": [], "ao_exp": [ 0.10992677e2, 0.18336830e2, -0.16366004e2, -0.62332348e1, 0.28035363e1, 0.10778099e1, 0.96965697 ], "exp": [ 2115.01542, 1676.18569, 1935.16717, 1504.97016, 4222.83691, 5296.17127, 273.36934 ], "ao_hyp": [], "hyp": [] } CP2 = { "ao": 0.964220 / 8.3143 * 32., "an": [ 0.532325e-4 / 8.3143 * 32., 0.672819e-5 / 8.3143 * 32., -0.768411e-8 / 8.3143 * 32., 0.275220e-11 / 8.3143 * 32. ], "pow": [1, 2, 3, 4], "ao_exp": [], "exp": [], "ao_hyp": [], "hyp": [] } helmholtz1 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for methanol of de Reuck and Craven (1993)", "__doi__": { "autor": "de Reuck, K.M. and Craven, R.J.B.", "title": "Methanol, International Thermodynamic Tables of the Fluid State - 12", "ref": "IUPAC, Blackwell Scientific Publications, London, 1993.", "doi": "" }, "R": 8.31448, "cp": CP1, "ref": "NBP", "Tref": 513.38, "rhoref": 8.78517 * M, "Tmin": Tt, "Tmax": 620.0, "Pmax": 800000.0, "rhomax": 35.57, "Pmin": 0.18629e-3, "rhomin": 28.23, "nr1": [ -0.280062505988e1, 0.125636372418e2, -0.130310563173e2, 0.326593134060e1, -0.411425343805e1, 0.346397741254e1, -0.836443967590e-1, -0.369240098923, 0.313180842152e-2, 0.603201474111, -0.231158593638, 0.106114844945, -0.792228164995e-1, -0.422419150975e-4, 0.758196739214e-2, -0.244617434701e-4, 0.115080328802e-5 ], "d1": [1, 1, 1, 1, 2, 2, 2, 2, 2, 3, 3, 3, 4, 4, 5, 6, 7], "t1": [0, 1, 2, 3, 1, 2, 3, 4, 6, 0, 3, 4, 0, 7, 1, 6, 7], "nr2": [ -0.125099747447e2, 0.270392835391e2, -0.212070717086e2, 0.632799472270e1, 0.143687921636e2, -0.287450766617e2, 0.185397216068e2, -0.388720372879e1, -0.416602487963e1, 0.529665875982e1, 0.509360272812, -0.330257604839e1, -0.311045210826, 0.273460830583, 0.518916583979, -0.227570803104e-2, 0.211658196182e-1, -0.114335123221e-1, 0.249860798459e-2 ], "d2": [1, 1, 1, 1, 2, 2, 2, 2, 3, 4, 5, 5, 5, 5, 6, 9, 6, 6, 4], "t2": [1, 2, 3, 4, 1, 2, 3, 5, 1, 2, 1, 2, 4, 5, 2, 5, 9, 14, 19], "c2": [2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 4, 4, 6], "gamma2": [1.01733510223052] * 16 + [1.03497071023039] * 2 + [1.05291203329783], "nr3": [ -0.819291988442e1, 0.478601004557, -0.444161392885, 0.179621810410, -0.687602278259, 0.240459848295e1, -0.688463987466e1, 0.113992982501e1 ], "d3": [1, 1, 1, 1, 1, 3, 3, 3], "t3": [0] * 8, "alfa3": [ 4.06934040892209, 8.20892015621185, 9.15601592007471, 83.8326275286616, 16.2773616356884, 27.705105527215, 16.2773616356884, 264.95250181898 ], "beta3": [ -3.8940745646517, -3.8940745646517, -3.8940745646517, -3.8940745646517, -3.8940745646517, -23.0649031906293, -23.0649031906293, -23.0649031906293 ], "gamma3": [ 1.54080254509371, 1.54080254509371, 1.54080254509371, 1.54080254509371, 1.54080254509371, 1.08389789427588, 1.08389789427588, 1.08389789427588 ], "epsilon3": [0] * 8, "exp1": [2, 3, 2, 4, 2, 3, 2, 4], "exp2": [1] * 8 } helmholtz2 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for methanol 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.3143, "cp": CP1, "ref": "NBP", "Tmin": Tt, "Tmax": 620.0, "Pmax": 800000.0, "rhomax": 40., "Pmin": 0.1, "rhomin": 40., "nr1": [ -2.4578394, 1.39060027, 8.56114069e-1, -4.20843418e-2, 3.63682442e-5, 7.05598662e-1 ], "d1": [1, 1, 1, 3, 7, 2], "t1": [1.5, 0.25, 1.25, 0.25, 0.875, 1.375], "nr2": [ 3.70573369e-1, 2.46303468, 1.50253790, 7.47553687e-2, -3.06417876e-1, -7.48402758e-1, -1.01432849e-1, 8.06830693e-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 } helmholtz3 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for methanol of Polt et al. (1992)", "__doi__": { "autor": "Polt, A., Platzer, B., and Maurer, G.", "title": "Parameter der thermischen Zustandsgleichung von Bender fuer 14 mehratomige reine Stoffe", "ref": "Chem. Technik 22(1992)6 , 216/224", "doi": "" }, "R": 8.3143, "cp": CP2, "ref": "NBP", "Tmin": 298., "Tmax": 703.0, "Pmax": 63000.0, "rhomax": 26.0625, "Pmin": 16.803, "rhomin": 24.576, "nr1": [ -0.412043979985e1, 0.541210456547e1, -0.974639417666, -0.909437999343, -0.143467597275, 0.557052459597e1, -0.697445416557e1, 0.860535902136, 0.244117735035e1, -0.449073510921e1, 0.223855290012e1, -0.71733653794, 0.876135006507, 0.151777405466, -0.233178058896, 0.140022534721e-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.412043979985e1, -0.541210456547e1, 0.974639417666, -0.4642672133, 0.944015617353, -0.449348200461 ], "d2": [0, 0, 0, 2, 2, 2], "t2": [3, 4, 5, 3, 4, 5], "c2": [2] * 6, "gamma2": [0.591872] * 6 } eq = helmholtz1, helmholtz2, helmholtz3 _surface = { "sigma": [0.22421, -0.21408, 0.083233], "exp": [1.3355, 1.677, 4.4402] } _melting = { "eq": 1, "Tref": Tt, "Pref": 0.187e-3, "Tmin": Tt, "Tmax": 620., "a1": [1], "exp1": [0], "a2": [5.320770e9, 4.524780e9, 3.888861e10], "exp2": [1, 1.5, 4], "a3": [], "exp3": [] } _vapor_Pressure = { "eq": 5, "ao": [-0.87414e1, 0.15035e1, -0.28720e1, -0.51345], "exp": [1., 1.5, 2.5, 5.] } _liquid_Density = { "eq": 1, "ao": [0.60230e-1, 0.18855e2, -0.27626e2, 0.11213e2, 0.69039], "exp": [0.1, 0.65, 0.79, 0.95, 4.4] } _vapor_Density = { "eq": 3, "ao": [ -0.81104, -0.55661e1, -0.79326e3, 0.19234e4, -0.29219e4, 0.29660e4, -0.13210e4 ], "exp": [0.25, 0.6, 3.5, 4.0, 5.0, 6.0, 7.0] } visco0 = { "eq": 0, "method": "_visco0", "__name__": "Xiang (2006)", "__doi__": { "autor": "Xiang, H.W., Huber, M.L. and Laesecke, A.", "title": "A New Reference Correlation for the Viscosity of Methanol", "ref": "J. Phys. Chem. Ref. Data 35, 1597 (2006)", "doi": "10.1063/1.2360605" }, "__test__": # Table 5, Pag 15 """ >>> st=Methanol(T=175.63, x=0.5) >>> print "%0.2f %0.3g %0.3g %0.4g %0.5g %0.4g" % (\ st.T, st.P.MPa, st.Liquido.rho, st.Liquido.mu.muPas, st.Gas.rho, st.Gas.mu.muPas) 175.63 1.86e-7 4.09e-6 0.005822 904.56 12.80 >>> st=Methanol(T=200, x=0.5) >>> print "%0.0f %0.4g %0.3g %0.4g %0.5g %0.4g" % (\ st.T, st.P.MPa, st.Liquido.rho, st.Liquido.mu.muPas, st.Gas.rho, st.Gas.mu.muPas) 200 6.098e-6 1.1754e-4 0.006563 880.28 4.506 >>> st=Methanol(T=250, x=0.5) >>> print "%0.0f %0.4g %0.3g %0.4g %0.5g %0.4g" % (\ st.T, st.P.MPa, st.Liquido.rho, st.Liquido.mu.muPas, st.Gas.rho, st.Gas.mu.muPas) 250 0.0008103 0.012577 0.008112 831.52 1.236 >>> st=Methanol(T=300, x=0.5) >>> print "%0.0f %0.5g %0.3g %0.4g %0.5g %0.4g" % (\ st.T, st.P.MPa, st.Liquido.rho, st.Liquido.mu.muPas, st.Gas.rho, st.Gas.mu.muPas) 300 0.018682 0.24623 0.009678 784.51 0.5291 >>> st=Methanol(T=350, x=0.5) >>> print "%0.0f %0.5g %0.3g %0.4g %0.5g %0.4g" % (\ st.T, st.P.MPa, st.Liquido.rho, st.Liquido.mu.muPas, st.Gas.rho, st.Gas.mu.muPas) 350 0.16172 1.9053 0.01118 735.84 0.2838 >>> st=Methanol(T=400, x=0.5) >>> print "%0.0f %0.5g %0.3g %0.4g %0.5g %0.4g" % (\ st.T, st.P.MPa, st.Liquido.rho, st.Liquido.mu.muPas, st.Gas.rho, st.Gas.mu.muPas) 400 0.77374 8.7343 0.01251 678.59 0.1714 >>> st=Methanol(T=450, x=0.5) >>> print "%0.0f %0.5g %0.3g %0.4g %0.5g %0.4g" % (\ st.T, st.P.MPa, st.Liquido.rho, st.Liquido.mu.muPas, st.Gas.rho, st.Gas.mu.muPas) 450 2.5433 30.831 0.01388 600.49 0.1058 >>> st=Methanol(T=500, x=0.5) >>> print "%0.0f %0.5g %0.3g %0.4g %0.5g %0.4g" % (\ st.T, st.P.MPa, st.Liquido.rho, st.Liquido.mu.muPas, st.Gas.rho, st.Gas.mu.muPas) 500 6.5250 109.88 0.01891 451.53 0.05748 >>> st=Methanol(T=512, x=0.5) >>> print "%0.0f %0.5g %0.3g %0.4g %0.5g %0.4g" % (\ st.T, st.P.MPa, st.Liquido.rho, st.Liquido.mu.muPas, st.Gas.rho, st.Gas.mu.muPas) 512 8.0195 202.99 0.02838 341.17 0.04174 """ # Table 6, Pag 16 """ >>> st=Methanol(T=180, P=1e4) >>> print "%0.2f %0.5g %0.4g" % (st.P.MPa, st.rho, st.mu.muPas) 0.01 900.27 10.44 >>> st=Methanol(T=200, P=1e5) >>> print "%0.2f %0.5g %0.4g" % (st.P.MPa, st.rho, st.mu.muPas) 0.10 880.34 4.510 >>> st=Methanol(T=220, P=4e5) >>> print "%0.2f %0.5g %0.4g" % (st.P.MPa, st.rho, st.mu.muPas) 0.40 860.76 2.460 >>> st=Methanol(T=280, P=1e6) >>> print "%0.2f %0.5g %0.4g" % (st.P.MPa, st.rho, st.mu.muPas) 1.00 804.12 0.7228 >>> st=Methanol(T=300, P=1e4) >>> print "%0.2f %0.5g %0.4g" % (st.P.MPa, st.rho, st.mu.muPas) 0.01 0.12955 0.009696 >>> st=Methanol(T=400, P=1e8) >>> print "%0.2f %0.5g %0.4g" % (st.P.MPa, st.rho, st.mu.muPas) 100 784.82 0.2846 >>> st=Methanol(T=500, P=8e8) >>> print "%0.2f %0.5g %0.4g" % (st.P.MPa, st.rho, st.mu.muPas) 800 954.52 0.3908 >>> st=Methanol(T=600, P=3e6) >>> print "%0.2f %0.5g %0.4g" % (st.P.MPa, st.rho, st.mu.muPas) 3.00 20.717 0.01971 """ } def _visco0(self, rho, T, fase): # FIXME: No sale rhoc = 273. ek = 577.87 sigma0 = 0.3408e-9 delta = 0.4575 sigmac = 0.7193422e-9 a = [ 1.16145, -0.14874, 0.52487, -0.77320, 2.16178, -2.43787, 0.95976e-3, 0.10225, -0.97346, 0.10657, -0.34528, -0.44557, -2.58055 ] b = [ -19.572881, 219.73999, -1015.3226, 2471.0125, -3375.1717, 2491.6597, -787.26086, 14.085455, -0.34664158 ] d = [ -1.181909, 0.5031030, -0.6268461, 0.5169312, -0.2351349, 5.3980235e-2, -4.9069617e-3 ] e = [ 0, 4.018368, -4.239180, 2.245110, -0.5750698, 2.3021026e-2, 2.5696775e-2, -6.8372749e-3, 7.2707189e-4, -2.9255711e-5 ] T_ = self.T / ek OmegaLJ = a[0] * T_**a[1] + a[2] * exp(a[3] * T_) + a[4] * exp( a[5] * T_) OmegaD = a[7] * T_**a[8] + a[9] * exp(a[10] * T_) + a[11] * exp( a[12] * T_) OmegaSM = OmegaLJ * (1 + delta**2 * OmegaD / (1 + a[6] * delta**6)) no = 5. * (self.M / Avogadro * Boltzmann * self.T / pi)**0.5 / (16 * sigma0**2 * OmegaSM) B = (sum([b[i] / T_**(0.25 * i) for i in range(7)]) + b[7] / T_**2.5 + b[8] / T_**5.5) * Avogadro * sigma0**3 C = 1.86222085e-3 * T_**3 * exp( 9.990338 / T_**0.5) * (Avogadro * sigma0**3)**2 ng = 1 + B * rho / self.M * 1000 + C * (rho / self.M * 1000)**2 Tr = self.T / self.Tc rhor = rho / rhoc sigmaHS = sigmac * (sum([d[i] / Tr**i for i in range(7)]) + sum([e[i] * rhor**(i) for i in range(1, 10)])) b = 2 * pi * Avogadro * sigmaHS**3 / 3 Xi = b * rho / self.M * 1000 / 4 g = (1 - 0.5 * Xi) / (1 - Xi)**3 ne = 1. / g + 0.8 * b * rho / self.M * 1000 + 0.761 * g * sigmaHS * b**2 * ( rho / self.M * 1000)**2 f = 1 / (1 + exp(5 * (rhor - 1))) n = no * (f * ng + (1 - f) * ne) return unidades.Viscosity(n) _viscosity = visco0, thermo0 = { "eq": 1, "__name__": "Perkins (2002)", "__doi__": { "autor": "", "title": "", "ref": "unpublished preliminary correlation, NIST, MLH, Aug. 2006", "doi": "" }, "Tref": 1., "kref": 1, "no": [5.7992e-7], "co": [1.7862], "Trefb": 513.38, "rhorefb": 8.78517, "krefb": 1., "nb": [0.405435, -0.293791, -0.289002, 0.226890, 0.579019e-1, -0.399576e-1], "tb": [0, 1] * 3, "db": [1, 1, 2, 2, 3, 3], "cb": [0] * 6, "critical": 3, "gnu": 0.63, "gamma": 1.239, "R0": 1.03, "Xio": 0.194e-9, "gam0": 0.0496, "qd": 0.342e-9, "Tcref": 768.9 } _thermal = thermo0,
class nC9(MEoS): """Multiparameter equation of state for n-nonane""" name = "nonane" CASNumber = "111-84-2" formula = "CH3-(CH2)7-CH3" synonym = "" _refPropName = "NONANE" _coolPropName = "n-Nonane" rhoc = unidades.Density(232.1417) Tc = unidades.Temperature(594.55) Pc = unidades.Pressure(2281.0, "kPa") M = 128.2551 # g/mol Tt = unidades.Temperature(219.7) Tb = unidades.Temperature(423.91) f_acent = 0.4433 momentoDipolar = unidades.DipoleMoment(0.07, "Debye") id = 13 Fi1 = { "ao_log": [1, 16.349], "pow": [0, 1], "ao_pow": [10.7927224829, -8.2418318753], "ao_exp": [24.926, 24.842, 11.188, 17.483], "titao": [1221 / Tc, 2244 / Tc, 5008 / Tc, 11724 / Tc] } Fi2 = { "ao_log": [1, 3.0], "pow": [0, 1], "ao_pow": [16.313913248, -102.160247463], "ao_exp": [], "titao": [], "ao_hyp": [18.0241, 38.1235, 53.3415, 0], "hyp": [0.263819696, 1.370586158, 2.848860483, 0] } lemmon = { "__type__": "Helmholtz", "__name__": "short Helmholtz equation of state for nonane 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": 600.0, "Pmax": 800000.0, "rhomax": 6.06, "Pmin": 0.00044, "rhomin": 6.05, "nr1": [1.1151, -2.7020, 0.83416, -0.38828, 0.1376, 0.00028185], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875], "nr2": [0.62037, 0.015847, -0.61726, -0.15043, -0.012982, 0.0044325], "d2": [2, 5, 1, 4, 3, 4], "t2": [0.625, 1.75, 3.625, 3.625, 14.5, 12.], "c2": [1, 1, 2, 2, 3, 3], "gamma2": [1] * 6 } GERG = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for nonane of Kunz and " "Wagner (2008).", "__doi__": { "autor": "Kunz, O., Wagner, W.", "title": "The GERG-2008 Wide-Range Equation of State for " "Natural Gases and Other Mixtures: An Expansion " "of GERG-2004", "ref": "J. Chem.Eng. Data 57(11) (2012) 3032-3091", "doi": "10.1021/je300655b" }, "R": 8.314472, "cp": Fi2, "ref": "OTO", "Tmin": Tt, "Tmax": 600.0, "Pmax": 800000.0, "rhomax": 6.06, "Pmin": 0.00044, "rhomin": 6.05, "nr1": [0.11151e1, -0.27020e1, 0.83416, -0.38828, 0.13760, 0.28185e-3], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875], "nr2": [0.62037, 0.015847, -0.61726, -0.15043, -0.012982, 0.0044325], "d2": [2, 5, 1, 4, 3, 4], "t2": [0.625, 1.75, 3.625, 3.625, 14.5, 12.], "c2": [1, 1, 2, 2, 3, 3], "gamma2": [1] * 6 } eq = lemmon, GERG _surface = {"sigma": [0.053388], "exp": [1.262]} _dielectric = { "eq": 3, "Tref": 273.16, "rhoref": 1000., "a0": [0.10924], "expt0": [-1.], "expd0": [1.], "a1": [44.53, 0.045], "expt1": [0, 1], "expd1": [1, 1], "a2": [286.27, 529.31, -83471, -90493], "expt2": [0, 1, 0, 1], "expd2": [2, 2, 3, 3] } _vapor_Pressure = { "eq": 3, "n": [-0.84804e1, 0.28640e1, -0.37414e1, -0.57479e1, 0.18799e1], "t": [1.0, 1.5, 2.3, 4.6, 5.0] } _liquid_Density = { "eq": 1, "n": [-0.43785, 0.37240e1, -0.23029e1, 0.18270e1, 0.38664], "t": [0.116, 0.32, 0.54, 0.8, 3.5] } _vapor_Density = { "eq": 2, "n": [-0.33199e1, -0.23900e1, -0.15307e2, -0.51788e2, -0.11133e3], "t": [0.461, 0.666, 2.12, 5.1, 11.0] } visco0 = { "__name__": "Huber (2004)", "__doi__": { "autor": "Huber, M.L., Laesecke, A. Xiang, H.W.", "title": "Viscosity correlations for minor constituent " "fluids in natural gas: n-octane, n-nonane and " "n-decane", "ref": "Fluid Phase Equilibria 224 (2004) 263-270", "doi": "10.1016/j.fluid.2004.07.012" }, "eq": 1, "omega": 1, "ek": 472.127, "sigma": 0.66383, "n_chapman": 0.021357, "collision": [0.340344, -0.466455], "Tref_virial": 472.127, "n_virial": [ -19.572881, 219.73999, -1015.3226, 2471.0125, -3375.1717, 2491.6597, -787.26086, 14.085455, -0.34664158 ], "t_virial": [0, -0.25, -0.5, -0.75, -1, -1.25, -1.5, -2.5, -5.5], "Tref_res": 594.55, "rhoref_res": 232.1417, "muref_res": 1000, "nr": [-0.314367e-1, 0.639384e-2, 0.326258e-1, -0.108922e-1], "tr": [1, 1, 2, 2], "dr": [2, 3, 2, 3], "CPf": 192.935, "CPg1": 2.66987, "CPgi": [1.32137 / 2.66987], "CPti": [-0.5] } _viscosity = visco0, thermo0 = { "__name__": "Huber (2005)", "__doi__": { "autor": "Huber, M.L., Perkins, R.A.", "title": "Thermal conductivity correlations for minor " "constituent fluids in natural gas: n-octane, " "n-nonane and n-decane", "ref": "Fluid Phase Equilibria 227 (2005) 47-55", "doi": "10.1016/j.fluid.2004.10.031" }, "eq": 1, "Toref": 594.55, "koref": 1, "no": [0.878765e-2, -0.413510e-1, 0.104791, -0.320032e-1], "to": [0, 1, 2, 3], "Tref_res": 594.55, "rhoref_res": 232.1417, "kref_res": 1, "nr": [0.490088e-2, 0.996486e-2, -0.807305e-2, 0.557431e-2], "tr": [0, -1, 0, 0], "dr": [1, 1, 2, 3], "critical": 3, "gnu": 0.63, "gamma": 1.239, "R0": 1.03, "Xio": 0.194e-9, "gam0": 0.0496, "qd": 1.043054e-9, "Tcref": 891.825 } _thermal = thermo0,
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_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] } 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 R14(MEoS): """Multiparameter equation of state for R14""" name = "tetrafluoromethane" CASNumber = "75-73-0" formula = "CF4" synonym = "R14" _refPropName = "R14" _coolPropName = "R14" rhoc = unidades.Density(625.66) Tc = unidades.Temperature(227.51) Pc = unidades.Pressure(3750.0, "kPa") M = 88.0046 # g/mol Tt = unidades.Temperature(89.54) Tb = unidades.Temperature(145.10) f_acent = 0.1785 momentoDipolar = unidades.DipoleMoment(0.0, "Debye") id = 218 CP1 = { "ao": 3.9465247, "an": [-.88586725e-2, 0.13939626e-3, -0.30056204e-6, 0.20504001e-9], "pow": [1, 2, 3, 4], "ao_exp": [], "exp": [] } platzer = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for R-14 of Platzer (1990)", "__doi__": { "autor": "Platzer, B., Polt, A., Maurer, G.", "title": "Thermophysical Properties of Refrigerants", "ref": "Berlin: Springer-Verlag, 1990.", "doi": "" }, "R": 8.31451, "cp": CP1, "ref": "NBP", "Tmin": 120.0, "Tmax": 623.0, "Pmax": 51000.0, "rhomax": 20.764, "nr1": [ -.334698748966, .586690904687, -.147068929692, .103999039623e1, -.245792025288e1, .799614557889, -.749498954929, .152177772502, -.293408331764, .717794502866, -.0426467444199, .226562749365, -0.391091694003, -0.257394804936e-1, 0.554844884782e-1, 0.610988261204e-2 ], "d1": [0, 0, 0, 1, 1, 1, 1, 1, 2, 2, 2, 3, 3, 4, 4, 5], "t1": [3, 4, 5, 0, 1, 2, 3, 4, 0, 1, 2, 0, 1, 0, 1, 1], "nr2": [ .334698748966, -.586690904687, .147068929692, -.190315426142, .716157133959, -.703161904626 ], "d2": [0, 0, 0, 2, 2, 2], "t2": [3, 4, 5, 3, 4, 5], "c2": [2] * 6, "gamma2": [0.99832625] * 6 } eq = platzer, _surface = {"sigma": [0.0423], "exp": [1.24]} _vapor_Pressure = { "eq": 3, "n": [-0.61905e1, -0.91398e1, 0.12192e2, -0.47215e1, -0.20439e1], "t": [1.0, 1.5, 1.64, 2.5, 7.3] } _liquid_Density = { "eq": 1, "n": [-0.10612e1, 0.44343e1, -0.38753e1, 0.29825e1, 0.30746], "t": [0.1, 0.24, 0.4, 0.6, 3.9] } _vapor_Density = { "eq": 2, "n": [-.55804e2, .10868e3, -.64257e2, -.11954e4, .36688e4, -.25956e4], "t": [0.713, 0.84, 1.0, 5.8, 6.3, 6.6] } 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": N2, "visco": "visco0", "thermo": "thermo0", "ek": 164.44, "sigma": 0.4543, "omega": 5, "psi": [1.10941, -0.0630268], "psi_d": [0, 1], "fint": [1.19864e-3, 1.90048e-7], "fint_t": [0, 1], "chi": [1.0442], "chi_d": [0], "critical": 3, "gnu": 0.63, "gamma": 1.239, "R0": 1.03, "Xio": 0.194e-9, "gam0": 0.0496, "qd": 2.26566e-10, "Tcref": 1.5 * Tc } _viscosity = trnECS, _thermal = trnECS,
class EthylBenzene(MEoS): """Multiparameter equation of state for ethylbenzene""" name = "ethylbenzene" CASNumber = "100-41-4" formula = "C8H10" synonym = "" _refPropName = "EBENZENE" _coolPropName = "EthylBenzene" rhoc = unidades.Density(291.) Tc = unidades.Temperature(617.12) Pc = unidades.Pressure(3622.4, "kPa") M = 106.165 # g/mol Tt = unidades.Temperature(178.2) Tb = unidades.Temperature(409.314) f_acent = 0.305 momentoDipolar = unidades.DipoleMoment(0.6, "Debye") id = 45 Fi1 = { "ao_log": [1, 4.2557889], "pow": [0, 1], "ao_pow": [5.70409, -0.52414353], "ao_exp": [9.7329909, 11.201832, 25.440749], "titao": [585 / Tc, 4420 / Tc, 1673 / 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": 60000.0, "rhomax": 9.124, "nr1": [ 0.0018109418, -0.076824284, 0.041823789, 1.5059649, -2.4122441, -0.47788846, 0.18814732 ], "d1": [5, 1, 4, 1, 1, 2, 3], "t1": [1, 1, 0.92, 0.27, 0.962, 1.033, 0.513], "nr2": [-1.0657412, -0.20797007, 1.1222031, -0.99300799, -0.027300984], "d2": [1, 3, 2, 2, 7], "t2": [2.31, 3.21, 1.26, 2.29, 1.0], "c2": [2, 2, 1, 2, 1], "gamma2": [1] * 5, "nr3": [1.3757894, -0.44477155, -0.07769742, -2.16719], "d3": [1, 1, 3, 3], "t3": [0.6, 3.6, 2.1, 0.5], "alfa3": [1.178, 1.07, 1.775, 15.45], "beta3": [2.437, 1.488, 4, 418.6], "gamma3": [1.2667, 0.4237, 0.8573, 1.15], "epsilon3": [0.5494, 0.7235, 0.493, 0.8566] } eq = zhou, _vapor_Pressure = { "eq": 3, "n": [-7.8411, 2.5921, -3.502, -2.7613], "t": [1.0, 1.5, 2.5, 5.4] } _liquid_Density = { "eq": 1, "n": [3.5146, -3.7537, 5.476, -3.4724, 1.2141], "t": [0.43, 0.83, 1.3, 1.9, 3.1] } _vapor_Density = { "eq": 2, "n": [-3.2877, -3.6071, -15.878, -53.363, -128.57], "t": [0.42, 0.98, 2.48, 5.9, 13.4] } visco0 = { "__name__": "Meng (2017)", "__doi__": { "autor": "Meng, X.Y., Cao, F.L., Wu, J.T., Vesovic, V.", "title": "Reference Correlation of the Viscosity of " "Ethylbenzene from the Triple Point to 673 K and " "up to 110 MPa", "ref": "J. Phys. Chem. Ref. Data 46(1) (2017) 013101", "doi": "10.1063/1.4973501" }, "eq": 1, "omega": 3, "collision": [-1.4933, 473.2, -57033], "sigma": 1, "n_chapman": 0.22115 / M**0.5, "Tref_res": 617.12, "rhoref_res": 2.741016 * M, "nr": [ -0.0376893, 0.168877, 17.9684, 3.57702e-11, 29.996, -8.00082, -25.7468 ], "dr": [209 / 30, 209 / 30, 29 / 30, 731 / 30, 59 / 30, 29 / 30, 455 / 300], "tr": [-0.5, 0.6, -0.5, 2.9, -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 # Special exponential term for residual viscosity, Eq 5 Ei = [-3.29316e-13, -2.92665e-13, 2.97768e-13, 1.76186e-18] ni = [4.6, 11.1, 5.6, 12.4] ki = [20.8, 10.6, 19.7, 21.9] Tr = T / 617.12 rhor = rho / self.M / 2.741016 # Eq 7 g = 0 for E, n, k in zip(Ei, ni, ki): g += E * rhor**n / Tr**k mur = g * exp(rhor**2) return muB * rho / self.M + mur _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": 617.12, "koref": 1e-3, "no_num": [-1.10708, 10.8026, -28.9015, 41.9227, 20.9133, -4.01492], "to_num": [0, 1, 2, 3, 4, 5], "no_den": [0.259475, -0.343879, 1], "to_den": [0, 1, 2], "Tref_res": 617.12, "rhoref_res": 291, "kref_res": 1e-3, "nr": [ -4.97837e1, 1.06739e2, -6.85137e1, 2.26133e1, -2.79455, 6.63073e1, -1.46279e2, 1.21439e2, -4.62245e1, 6.58554 ], "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.056, "qd": 0.706e-9, "Tcref": 925.7 } _thermal = thermo0,
class CO(MEoS): """Multiparameter equation of state for carbon monoxide""" name = "carbon monoxide" CASNumber = "630-08-0" formula = "CO" synonym = "" _refPropName = "CO" _coolPropName = "CarbonMonoxide" rhoc = unidades.Density(303.909585) Tc = unidades.Temperature(132.86) Pc = unidades.Pressure(3494.0, "kPa") M = 28.0101 # g/mol Tt = unidades.Temperature(68.16) Tb = unidades.Temperature(81.64) f_acent = 0.0497 momentoDipolar = unidades.DipoleMoment(0.1, "Debye") id = 48 Fi1 = { "ao_log": [1, 2.5], "pow": [0, 1, -1.5], "ao_pow": [-3.3728318564, 3.3683460039, -9.111274701235156e-5], "ao_exp": [1.0128], "titao": [3089 / Tc] } Fi2 = { "ao_log": [1, 2.50055], "pow": [0, 1], "ao_pow": [10.813340744, -19.834733959], "ao_exp": [], "titao": [], "ao_hyp": [1.02865, 0.00493, 0, 0], "hyp": [11.6698028, 5.302762306, 0, 0] } CP3 = { "ao": 0.36028218e1, "an": [ -0.20871594e5, 0.89208708e3, -0.14157993e2, -0.34021345e-3, 0.44616091e-6, -0.15154703e-9 ], "pow": [-3, -2, -1.001, 1, 2, 3], "ao_exp": [0.90426143], "exp": [30000], "ao_hyp": [], "hyp": [] } lemmon = { "__type__": "Helmholtz", "__name__": "short Helmholtz equation of state for carbon monoxide 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., "Pmax": 100000.0, "rhomax": 33.84, "Pmin": 15.45, "rhomin": 30.33, "nr1": [0.90554, -2.4515, 0.53149, 0.024173, 0.072156, 0.00018818], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875], "nr2": [0.19405, -0.043268, -0.12778, -0.027896, -0.034154, 0.016329], "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 } mccarty = { "__type__": "MBWR", "__name__": "MBWR equation of state for carbon monoxide of McCarty " "(1989)", "__doi__": { "autor": "McCarty, R.D.", "title": "Correlations for the Thermophysical Properties " "of Carbon Monoxide", "ref": "NIST, Boulder, CO, 1989", "doi": "" }, "R": 8.31434, "cp": CP3, "ref": "NBP", "Tmin": Tt, "Tmax": 1000., "Pmax": 30000.0, "rhomax": 30.25, "Pmin": 15.423, "rhomin": 30.249, "b": [ None, 0.8845582109949e-2, -0.2236741566840, 0.1742275796442e1, -0.2169146998363e3, 0.1721504267082e4, -0.3990514770703e-4, 0.1036880040451, -0.3376308165071e2, 0.2061895161095e5, 0.2993711656350e-5, 0.1856003597097e-2, -0.2114419664527, -0.2436986935194e-5, -0.1858029609177e-2, -0.1734563867767e1, 0.1509970839260e-3, -0.2282721433205e-5, 0.2202780295674e-2, -0.3313357789163e-4, -0.1473412120276e5, -0.3141136651147e6, -0.1451168999234e3, 0.6323441221817e5, -0.2203560539926, -0.2087738308480e2, -0.1508165207553e-2, 0.2740740634030e1, 0.8687687989627e-6, -0.1451419251928e-3, -0.3040346241285e-8, 0.4712050805815e-8, -0.2639772456566e-5 ] } GERG = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for carbon monoxide 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": 500., "Pmax": 100000.0, "rhomax": 33.84, "Pmin": 15.45, "rhomin": 30.33, "nr1": [ 0.92310041400851, -0.248858452058e1, 0.58095213783396, 0.028859164394654, 0.070256257276544, 0.21687043269488e-3 ], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.25, 0.125, 1.5, 1.375, 0.25, 0.875], "nr2": [ 0.13758331015182, -0.51501116343466e-1, -0.14865357483379, -0.03885710088681, -0.029100433948943, 0.14155684466279e-1 ], "d2": [2, 5, 1, 4, 3, 4], "t2": [0.625, 1.75, 3.625, 3.625, 14.5, 12], "c2": [1, 1, 2, 2, 3, 3], "gamma2": [1] * 6 } eq = lemmon, mccarty, GERG _surface = {"sigma": [0.02843], "exp": [1.148]} _melting = { "eq": 1, "Tref": 1, "Pref": 1000, "Tmin": Tt, "Tmax": 1000.0, "a1": [-142.941, 0.0195608], "exp1": [0, 2.10747], "a2": [], "exp2": [], "a3": [], "exp3": [] } _vapor_Pressure = { "eq": 3, "n": [-0.61192e1, 0.10411e1, -0.62162e1, 0.10437e2, -0.76813e1], "t": [1.0, 1.5, 3.9, 4.6, 5.4] } _liquid_Density = { "eq": 1, "n": [0.29570e1, -0.42880e1, 0.87643e1, -0.84001e1, 0.36372e1], "t": [0.398, 0.735, 1.08, 1.5, 1.9] } _vapor_Density = { "eq": 2, "n": [ -0.25439e1, -0.55601e1, -0.85276e1, -0.51163e1, -0.17701e2, -0.29858e2 ], "t": [0.395, 1.21, 3.0, 3.5, 6.0, 8.0] }
class DMC(MEoS): """Multiparameter equation of state for dimethyl carbonate""" name = "dimethyl carbonate" CASNumber = "616-38-6" formula = "C3H6O3" synonym = "" _refPropName = "DMC" _coolPropName = "DimethylCarbonate" 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": [] } zhou = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for DMC of Zhou (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 = zhou, _vapor_Pressure = { "eq": 3, "n": [-8.3197, 3.4260, -3.5905, -3.3194], "t": [1.0, 1.5, 2.3, 4.7] } _liquid_Density = { "eq": 1, "n": [1.1572, 4.969, -14.451, 27.569, -26.223, 10.526], "t": [0.27, 0.77, 1.29, 1.85, 2.46, 3.16] } _vapor_Density = { "eq": 2, "n": [-0.54715, -5.19277, -94.048, 327.21, -676.871, 716.072, -379.799], "t": [0.197, 0.6, 2.86, 3.65, 4.5, 5.4, 6.4] }
class Trans_2_butene(MEoS): """Multiparameter equations of state for trans-butene""" name = "trans-butene" CASNumber = "624-64-6" formula = "CH3-CH=CH-CH3" synonym = "" _refPropName = "T2BUTENE" _coolPropName = "trans-2-Butene" rhoc = unidades.Density(236.37592616) Tc = unidades.Temperature(428.61) Pc = unidades.Pressure(4027.3, "kPa") M = 56.10632 # g/mol Tt = unidades.Temperature(167.6) Tb = unidades.Temperature(274.03) f_acent = 0.21 momentoDipolar = unidades.DipoleMoment(0.0, "Debye") id = 26 Fi1 = {"ao_log": [1, 2.9988], "pow": [0, 1], "ao_pow": [0.5917816, 2.1427758], "ao_exp": [5.3276, 13.29, 9.6745, 0.40087], "titao": [362/Tc, 1603/Tc, 3729/Tc, 4527/Tc]} lemmon = { "__type__": "Helmholtz", "__name__": "short Helmholtz equation of state for 1-butene of Lemmon " "and Ihmels (2005)", "__doi__": {"autor": "Lemmon, E.W., Ihmels, E.C.", "title": "Thermodynamic properties of the butenes: Part " "II. Short fundamental equations of state", "ref": "Fluid Phase Equilibria 228-229 (2005) 173-187", "doi": "10.1016/j.fluid.2004.09.004"}, "R": 8.314472, "cp": Fi1, "ref": "NBP", "Tmin": Tt, "Tmax": 525.0, "Pmax": 50000.0, "rhomax": 13.141, "nr1": [0.81107, -2.8846, 1.0265, 0.016591, 0.086511, 0.00023256], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.12, 1.3, 1.74, 2.1, 0.28, 0.69], "nr2": [0.22654, -0.072182, -0.24849, -0.071374, -0.024737, 0.011843], "d2": [2, 5, 1, 4, 3, 4], "t2": [0.75, 2., 4.4, 4.7, 15., 14.], "c2": [1, 1, 2, 2, 3, 3], "gamma2": [1]*6} eq = lemmon, _PR = [-0.1594, -16.6119] _vapor_Pressure = { "eq": 3, "n": [-7.6226, 7.9421, -6.9631, -6.5517, 3.9584], "t": [1.0, 1.5, 1.65, 4.8, 5.3]} _liquid_Density = { "eq": 1, "n": [12.452, -34.419, 52.257, -42.889, 15.463], "t": [0.52, 0.73, 0.97, 1.24, 1.5]} _vapor_Density = { "eq": 2, "n": [-3.1276, -6.0548, -18.243, -60.842, 135.95, -182.70], "t": [0.412, 1.24, 3.2, 7.0, 10.0, 11.0]}
class R227ea(MEoS): """Multiparameter equation of state for R227ea""" name = "1,1,1,2,3,3,3-heptafluoropropane" CASNumber = "431-89-0" formula = "CF3CHFCF3" synonym = "R227ea" _refPropName = "R227EA" _coolPropName = "R227EA" rhoc = unidades.Density(594.2508657) Tc = unidades.Temperature(374.9) Pc = unidades.Pressure(2925.0, "kPa") M = 170.02886 # g/mol Tt = unidades.Temperature(146.35) Tb = unidades.Temperature(256.81) f_acent = 0.357 momentoDipolar = unidades.DipoleMoment(1.456, "Debye") # id = 1872 Fi1 = { "ao_log": [1, 3], "pow": [0, 1], "ao_pow": [-15.8291124137, 11.0879509962], "ao_exp": [11.43, 12.83], "titao": [403 / Tc, 1428 / Tc] } lemmon = { "__type__": "Helmholtz", "__name__": "short Helmholtz equation of state for R-227ea 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": Fi1, "ref": "IIR", "Tmin": Tt, "Tmax": 475.0, "Pmax": 60000.0, "rhomax": 11.05, "nr1": [2.024341, -2.605930, 0.4957216, -0.8240820, 0.06543703], "d1": [1, 1, 2, 2, 4], "t1": [0.34, 0.77, 0.36, 0.9, 1], "nr2": [-1.02461, .6247065, .2997521, -.353917, -1.232043, -.8824483], "d2": [1, 3, 6, 6, 2, 3], "t2": [2.82, 2.1, 0.9, 1.13, 3.8, 2.75], "c2": [1, 1, 1, 1, 2, 2], "gamma2": [1] * 6, "nr3": [ 0.1349661, -0.2662928, 0.1764733, 0.01536163, -0.004667185, -11.70854, 0.9114512 ], "d3": [1, 2, 1, 1, 4, 2, 1], "t3": [1.5, 2.5, 2.5, 5.4, 4, 1, 3.5], "alfa3": [0.83, 2.19, 2.44, 3.65, 8.88, 8.23, 2.01], "beta3": [1.72, 5.2, 2.31, 1.02, 5.63, 50.9, 1.56], "gamma3": [0.414, 1.051, 1.226, 1.7, 0.904, 1.42, 0.926], "epsilon3": [1.13, 0.71, 1.2, 1.7, 0.546, 0.896, 0.747] } eq = lemmon, _surface = { "sigma": [0.06127, -0.009516, -0.00192], "exp": [1.192, 0.9795, 1.421] } _vapor_Pressure = { "eq": 3, "n": [-7.7961, 2.1366, -2.6023, -5.7444, 2.3982], "t": [1.0, 1.5, 2.2, 4.8, 6.2] } _liquid_Density = { "eq": 1, "n": [0.20032e1, 0.49235, 0.13738, 0.21057, -0.12834], "t": [0.345, 0.74, 1.2, 2.6, 7.2] } _vapor_Density = { "eq": 2, "n": [-.2135e1, -.68425e1, -.21447e2, -.20457e3, .51795e3, -.45908e3], "t": [0.324, 1.03, 3.0, 7.4, 9.0, 10.0] } trnECS = { "__name__": "Huber (2003)", "__doi__": { "autor": "Huber, M.L., Laesecke, A., Perkins, R.A.", "title": "Model for the Viscosity and Thermal Conductivity " "of Refrigerants, Including a New Correlation for " "the Viscosity of R134a", "ref": "Ind. Eng. Chem. Res., 42(13) (2003) 3163-3178", "doi": "10.1021/ie0300880" }, "eq": "ecs", "ref": C3, "visco": "visco1", "thermo": "thermo0", "ek": 289.34, "sigma": 0.5746, "omega": 5, "psi": [0.76758, 0.254482, -5.33748e-2], "psi_d": [0, 1, 2], "fint": [1.42313e-3, 8.31496e-9], "fint_t": [0, 1], "chi": [1.3122, -8.74448e-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 R1234ze(MEoS): """Multiparameter equation of state for R1234ze""" name = "trans-1,3,3,3-tetrafluoropropene" CASNumber = "29118-24-9" formula = "CHF=CHCF3" synonym = "R-1234ze" _refPropName = "R1234ZE" rhoc = unidades.Density(489.238464) Tc = unidades.Temperature(382.513) Pc = unidades.Pressure(3634.9, "kPa") M = 114.0416 # g/mol Tt = unidades.Temperature(168.62) Tb = unidades.Temperature(254.177) f_acent = 0.313 momentoDipolar = unidades.DipoleMoment(1.27, "Debye") Fi1 = { "ao_log": [1, 3], "pow": [0, 1], "ao_pow": [-12.558347537, 8.7912297624], "ao_exp": [9.3575, 10.717], "titao": [513 / Tc, 1972 / Tc] } Fi2 = { "ao_log": [1, 5.8887], "pow": [], "ao_pow": [], "ao_exp": [7.0804, 9.3371, 2.5577], "titao": [620 / Tc, 1570 / Tc, 3953 / Tc] } Fi3 = { "ao_log": [1, 3], "pow": [0, 1], "ao_pow": [-10.8724711, -30.1326538], "ao_exp": [6.07536, 9.95795], "titao": [289 / Tc, 1303 / Tc] } thol = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for R1234ze of Thol (2016)", "__doi__": { "autor": "Thol, M., Lemmon, E.W.", "title": "Equation of State for the Thermodynamic" "Properties of trans-1,3,3,3-Tetrafluoroporpene" "[R-1234ze(E)]", "ref": "Int. J. Thermophys. 37(3) (2016) 28", "doi": "10.1007/s10765-016-2040-6" }, "R": 8.3144621, "cp": Fi1, "ref": "IIR", "Tmin": Tt, "Tmax": 420.0, "Pmax": 20000.0, "rhomax": 13.26, "nr1": [0.03982797, 1.812227, -2.537512, -0.5333254, 0.1677031], "d1": [4, 1, 1, 2, 3], "t1": [1.0, 0.223, 0.755, 1.24, 0.44], "nr2": [-1.323801, -0.6694654, 0.8072718, -0.7740229, -0.01843846], "d2": [1, 3, 2, 2, 7], "t2": [2., 2.2, 1.2, 1.5, 0.9], "c2": [2, 2, 1, 2, 1], "gamma2": [1] * 5, "nr3": [ 1.407916, -0.4237082, -0.2270068, -0.805213, 0.00994318, -0.008798793 ], "d3": [1, 1, 3, 3, 2, 1], "t3": [1.33, 1.75, 2.11, 1.0, 1.5, 1.0], "alfa3": [1.0, 1.61, 1.24, 9.34, 5.78, 3.08], "beta3": [1.21, 1.37, 0.98, 171, 47.4, 15.4], "gamma3": [0.943, 0.642, 0.59, 1.2, 1.33, 0.64], "epsilon3": [0.728, 0.87, 0.855, 0.79, 1.3, 0.71], "nr4": [] } mclinden = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for R1234ze of McLinden " "(2010)", "__doi__": { "autor": "McLinden, M.O., Thol, M., Lemmon, E.W.", "title": "Thermodynamic Properties of trans-1,3,3,3-" "Tetrafluoropropene [R1234ze(E)]: Measurements " "of Density and Vapor Pressure and a " "Comprehensive Equation of State", "ref": "International Refrigeration and Air Conditioning " "Conference at Purdue, July 12-15, 2010.", "doi": "10.0000_docs.lib.purdue.edu_generic-99DA7EA2C877" }, "R": 8.314472, "cp": Fi2, "ref": "IIR", "Tc": 382.52, "rhoc": 4.29, "M": 114.0415928, "Tmin": Tt, "Tmax": 420.0, "Pmax": 20000.0, "rhomax": 13.20, "nr1": [0.055563, 1.66927, -2.53408, -0.475075, 0.190055], "d1": [4, 1, 1, 2, 3], "t1": [1.0, 0.34, 0.91, 1.23, 0.46], "nr2": [-1.25154, -0.742195, 0.537902, -0.741246, -0.0355064], "d2": [1, 3, 2, 2, 7], "t2": [2.26, 2.50, 2.0, 2.24, 0.9], "c2": [2, 2, 1, 2, 1], "gamma2": [1] * 5, "nr3": [1.58506, -0.502086, -0.19136, -0.975576], "d3": [1, 1, 3, 3], "t3": [1.06, 1.79, 3.75, 0.92], "alfa3": [1.02, 1.34, 1.08, 6.41], "beta3": [1.19, 2.29, 1.15, 131.8], "gamma3": [1.14, 0.667, 0.505, 1.22], "epsilon3": [0.711, 0.914, 0.694, 0.731] } akasaka = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for R1234yf of Akasaka " "(2011)", "__doi__": { "autor": "Akasaka, R.", "title": "New Fundamental Equations of State with a Common" " Functional Form for 2,3,3,3-Tetrafluoropropene " "(R-1234yf) and trans-1,3,3,3-Tetrafluoropropene " "(R-1234ze(E))", "ref": "Int. J. Thermophys. 32(6) (2011) 1125-1147", "doi": "10.1007/s10765-011-0992-0" }, "R": 8.314472, "cp": Fi3, "ref": "IIR", "Tmin": 240., "Tmax": 420.0, "Pmax": 15000.0, "rhomax": 13.20, "nr1": [8.5579765, -9.4701332, -0.25013623, 0.1378987, 0.012177113], "d1": [1, 1, 1, 2, 5], "t1": [0.66886, 0.83392, 1.6982, 1.8030, 0.36657], "nr2": [ -0.14227996, 0.10096648, 0.017504319, -0.017627303, -0.014705120, 0.37202269, -0.30138266, -0.092927274, 0.087051177, 0.01811377, -0.016018424, 0.005380986 ], "d2": [1, 3, 5, 7, 1, 2, 2, 3, 4, 2, 3, 5], "t2": [ 3.8666, 1.0194, 0, 1.1655, 8.3101, 6.1459, 8.3495, 6.0422, 7.444, 15.433, 21.543, 15.499 ], "c2": [1, 1, 1, 1, 2, 2, 2, 2, 2, 3, 3, 3], "gamma2": [1] * 12 } eq = thol, mclinden, akasaka _surface = { "__doi__": { "autor": "Tanaka, K., Higashi, Y.", "title": "Surface Tensions of trans-1,3,3,3-Tetrafluoropropene " "and trans-1,3,3,3-Tetrafluoropropene + Difluoromethane " "Mixture", "ref": "J. Chem. Eng. Japan 46(6) (2013) 371-375", "doi": "10.1252/jcej.13we021" }, "sigma": [0.05681], "exp": [1.23] } _vapor_Pressure = { "eq": 3, "n": [-7.5888, 1.9696, -2.0827, -4.1238], "t": [1.0, 1.5, 2.2, 4.6] } _liquid_Density = { "eq": 1, "n": [1.1913, 2.2456, -1.7747, 1.3096], "t": [0.27, 0.7, 1.25, 1.9] } _vapor_Density = { "eq": 2, "n": [-1.0308, -5.0422, -11.5, -37.499, -77.945], "t": [0.24, 0.72, 2.1, 4.8, 9.5] } thermo0 = { "__name__": "Perkins (2011)", "__doi__": { "autor": "Perkins, R.A., Huber, M.L.", "title": "Measurement and Correlation of the Thermal " "Conductivity of 2,3,3,3-Tetrafluoroprop-1-ene " "(R1234yf) and trans-1,3,3,3-Tetrafluoropropene " "(R1234ze(E))", "ref": "J. Chem. Eng. Data 56(12) (2011) 4868-4874", "doi": "10.1021/je200811n" }, "eq": 1, "Pc": 3.6363e6, "Toref": 382.52, "koref": 1, "no": [-0.0103589, 0.0308929, 0.000230348], "to": [0, 1, 2], "Tref_res": 382.52, "rhoref_res": 489.24, "kref_res": 1., "nr": [ -0.0428296, 0.0927099, -0.0702107, 0.0249708, -0.00301838, 0.0434288, -0.0605844, 0.0440187, -0.0155082, 0.0021019 ], "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.03, "Xio": 0.194e-9, "gam0": 0.0496, "qd": 5.835e-10, "Tcref": 573.78 } _thermal = thermo0,
class pH2(MEoS): """Multiparameter equation of state for hydrogen (para)""" name = "parahydrogen" CASNumber = "1333-74-0p" formula = "H2" synonym = "R-702p" rhoc = unidades.Density(31.32274344) Tc = unidades.Temperature(32.938) Pc = unidades.Pressure(1285.8, "kPa") M = 2.01588 # g/mol Tt = unidades.Temperature(13.8033) Tb = unidades.Temperature(20.271) f_acent = -0.219 momentoipolar = unidades.DipoleMoment(0.0, "Debye") id = 1 Fi1 = { "ao_log": [1, 1.5], "pow": [0, 1], "ao_pow": [-1.4485891134, 1.884521239], "ao_exp": [4.30256, 13.0289, -47.7365, 50.0013, -18.6261, 0.993973, 0.536078], "titao": [ 15.1496751472, 25.0925982148, 29.4735563787, 35.4059141417, 40.724998482, 163.7925799988, 309.2173173842 ] } CP1 = { "ao": 2.4995169, "an": [ -0.11125185e-2, 0.27491461e-3, -0.10005269e-4, 0.22695404e-8, -0.21031029e-12 ], "pow": [1, 1.5, 2, 3, 4], "ao_exp": [0.12353388e2, -0.17777676e2, 0.64309174e1, 0.73347521e1], "exp": [598, 778, 1101, 6207], "ao_hyp": [], "hyp": [] } helmholtz1 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for parahydrogen of Leachman et al. (2007)", "__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=pH2(T=13.8033, 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.8033 7.0410 76.977 0.12555 −53.741 396.31 −3.0840 29.521 5.1313 6.2265 6.9241 10.534 1263.1 305.65 """, # Table 13, Pag 745 "R": 8.314472, "cp": Fi1, "ref": "NBP", "Tmin": Tt, "Tmax": 1000.0, "Pmax": 2000000.0, "rhomax": 104.0, "Pmin": 7.041, "rhomin": 38.185, "nr1": [-7.33375, 0.01, 2.60375, 4.66279, 0.682390, -1.47078, 0.135801], "d1": [1, 4, 1, 1, 2, 2, 3], "t1": [0.6855, 1, 1, 0.489, 0.774, 1.133, 1.386], "nr2": [-1.05327, 0.328239], "d2": [1, 3], "t2": [1.619, 1.162], "c2": [1, 1], "gamma2": [1] * 2, "nr3": [-0.0577833, 0.0449743, 0.0703464, -0.0401766, 0.119510], "d3": [2, 1, 3, 1, 1], "t3": [3.96, 5.276, 0.99, 6.791, 3.19], "alfa3": [1.7437, 0.5516, 0.0634, 2.1341, 1.777], "beta3": [0.194, 0.2019, 0.0301, 0.2383, 0.3253], "gamma3": [0.8048, 1.5248, 0.6648, 0.6832, 1.493], "epsilon3": [1.5487, 0.1785, 1.28, 0.6319, 1.7104], "nr4": [] } MBWR = { "__type__": "MBWR", "__name__": "MBWR equation of state for parahydrogen 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": 44.0, "Pmin": 7.042, "rhomin": 38.21, "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 ] } eq = helmholtz1, MBWR _surface = {"sigma": [0.005314], "exp": [1.06]} _dielectric = { "eq": 3, "Tref": 273.16, "rhoref": 1000., "a0": [], "expt0": [], "expd0": [], "a1": [2.0297, 0.0069], "expt1": [0, 1], "expd1": [1, 1], "a2": [0.181, 0.021, -7.4], "expt2": [0, 1, 0], "expd2": [2, 2, 3] } _melting = {"Tmin": Tt, "Tmax": 400.0} _sublimation = { "eq": 2, "Tref": 1, "Pref": 0.133332237, "Tmin": Tt, "Tmax": Tt, "a1": [4.009857354, -90.77568949], "exp1": [0, -1], "a2": [], "exp2": [], "a3": [2.48983094], "exp3": [1] } _vapor_Pressure = { "eq": 5, "ao": [-0.487767e1, 0.103359e1, 0.826680, -0.129412], "exp": [1.0, 1.5, 2.65, 7.4] } _liquid_Density = { "eq": 1, "ao": [-0.13509, 0.40739e1, -0.53985e1, 0.55230e1, -0.23643e1], "exp": [0.15, 0.44, 0.7, 0.99, 1.31] } _vapor_Density = { "eq": 3, "ao": [-0.57545e1, 0.38153e1, -0.12293e2, 0.15095e2, -0.17295e2, -0.34190e2], "exp": [0.53, 0.7, 1.7, 2.4, 3.3, 10] } visco0 = { "eq": 0, "method": "_visco0", "__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": "Natl. Bur. Stand., Tech. Note 617, 1972.", "doi": "" } } _viscosity = visco0, def _visco0(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.1099981128e2, 0.1895876508e2, -0.3813005056e3, 0.5950473265e2, 0.1099399458e1, 0.8987269839e1, 0.1231422148e4, 0.311 ] 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") @classmethod def _Melting_Pressure(cls, T=None): Tref = 1 Pref = 1000 Tita = T / Tref # a = [-0.265289115e2, 0.248578596, -0.21272389e2, 0.125746643] # expo = [0, 0.1764739e1, 0, 0.1955e1] if T > 22: suma = -0.265289115e2 + 0.248578596 * Tita**0.1764739e1 else: suma = -0.21272389e2 + 0.125746643 * Tita**0.1955e1 return unidades.Pressure(suma * Pref, "kPa") 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=pH2(T=298.15, rho=0) >>> print "%0.5g" % st.k.mWmK 192.38 >>> st=pH2(T=298.15, rho=0.80844) >>> print "%0.5g" % st.k.mWmK 192.81 >>> st=pH2(T=298.15, rho=14.4813) >>> print "%0.5g" % st.k.mWmK 207.85 >>> st=pH2(T=35, rho=0) >>> print "%0.5g" % st.k.mWmK 27.222 >>> st=pH2(T=35, rho=30) >>> print "%0.5g" % st.k.mWmK 70.335 >>> st=pH2(T=35, rho=30) >>> print "%0.5g" % st.k.mWmK 68.611 >>> st=pH2(T=18, rho=0) >>> print "%0.5g" % st.k.mWmK 13.643 >>> st=pH2(T=18, rho=75) >>> print "%0.5g" % st.k.mWmK 100.52 """, # Table 4, Pag 8 "Tref": 1.0, "kref": 1e-3, "no": [ -1.24500e3, 9.41806e3, -3.05098e2, 6.88449, -5.58871e-2, 2.79243e-4, -4.06944e-7, 3.42309e-10 ], "co": [0, 1, 2, 3, 4, 5, 6, 7], "noden": [ 1.42304e4, -5.88749e2, 1.45983e1, -1.34830e-1, 6.19047e-4, -9.21777e-7, 7.83099e-10 ], "coden": [0, 1, 2, 3, 4, 5, 6], "Trefb": 32.938, "rhorefb": 15.538, "krefb": 1., "nb": [ 0.265975e-1, -0.133826e-2, 0.130219e-1, -0.567678e-2, -0.92338e-4, -0.121727e-2, 0.366663e-2, 0.388715e-2, -0.921055e-2, 0.400723e-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.5e-9, "Tcref": 49.407 } thermo1 = { "eq": 0, "method": "_thermo1", "__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": "Natl. Bur. Stand., Tech. Note 617, 1972.", "doi": "" } } def _thermo1(self, rho, T, fase): # TODO: return unidades.ThermalConductivity(0) _thermal = thermo0, thermo1
class oXylene(MEoS): """Multiparameter equation of state for o-xylene """ name = "o-xylene " CASNumber = "95-47-6" formula = "C8H10" synonym = "1,2-dimethylbenzene" rhoc = unidades.Density(285) Tc = unidades.Temperature(630.259) Pc = unidades.Pressure(3737.5, "kPa") M = 106.165 # g/mol Tt = unidades.Temperature(247.985) Tb = unidades.Temperature(417.521) f_acent = 0.312 momentoDipolar = unidades.DipoleMoment(0.63, "Debye") id = 42 Fi1 = { "ao_log": [1, 2.748798], "pow": [0, 1], "ao_pow": [10.137376, -0.91282993], "ao_exp": [4.754892, 6.915052, 15.84813, 10.93886], "titao": [225 / Tc, 627 / Tc, 1726 / Tc, 4941 / Tc] } helmholtz1 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for o-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": 70000.0, "rhomax": 8.648, "Pmin": 0.0228, "rhomin": 8.647, "nr1": [ 0.0036765156, -0.13918171, 0.014104203, 1.5398899, -2.3600925, -0.44359159, 0.19596977 ], "d1": [5, 1, 4, 1, 1, 2, 3], "t1": [1.0, 0.6, 0.91, 0.3, 0.895, 1.167, 0.435], "nr2": [-1.0909408, -0.21890801, 1.1179223, -0.93563815, -0.018102996], "d2": [1, 3, 2, 2, 7], "t2": [2.766, 3.8, 1.31, 3.0, 0.77], "c2": [2, 2, 1, 2, 1], "gamma2": [1] * 5, "nr3": [1.4172368, -0.57134695, -0.081944041, -40.682878], "d3": [1, 1, 3, 3], "t3": [1.41, 4.8, 1.856, 2.0], "alfa3": [1.1723, 1.095, 1.6166, 20.4], "beta3": [2.442, 1.342, 3.0, 450.0], "gamma3": [1.2655, 0.3959, 0.7789, 1.162], "epsilon3": [0.552, 0.728, 0.498, 0.894] } eq = helmholtz1, _vapor_Pressure = { "eq": 5, "ao": [-7.2834, -1.5813, 7.6516, -7.9953, -2.2277], "exp": [1.0, 1.5, 1.9, 2.4, 6.0] } _liquid_Density = { "eq": 1, "ao": [0.9743, 16.511, -52.934, 87.962, -71.719, 22.569], "exp": [0.3, 0.96, 1.4, 1.9, 2.4, 3.0] } _vapor_Density = { "eq": 3, "ao": [-1.29038, -33.3428, 142.046, -292.211, 293.950, -159.504, -88.2170], "exp": [0.32, 1.14, 1.7, 2.2, 2.8, 3.5, 9.8] }
class Ethanol(MEoS): """Multiparameter equation of state for ethanol""" name = "ethanol" CASNumber = "64-17-5" formula = "C2H6O" synonym = "" _refPropName = "ETHANOL" _coolPropName = "Ethanol" rhoc = unidades.Density(273.1858492) Tc = unidades.Temperature(514.71) Pc = unidades.Pressure(6268., "kPa") M = 46.06844 # g/mol Tt = unidades.Temperature(159) Tb = unidades.Temperature(351.57) f_acent = 0.646 momentoDipolar = unidades.DipoleMoment(1.6909, "Debye") id = 134 Fi1 = {"ao_log": [1, 3.43069], "pow": [0, 1], # Use custom to get the IIR referenve value # "ao_pow": [-12.757491395, 9.3938494264], "ao_pow": [-12.7531, 9.3909425], "ao_exp": [2.14326, 5.09206, 6.60138, 5.70777], "titao": [420.4/Tc, 1334/Tc, 1958/Tc, 4420/Tc]} CP1 = {"ao": 6.4112, "an": [], "pow": [], "ao_exp": [1.95988750679, 7.60084166080, 3.89583440622, 4.23238091363], "exp": [694, 1549, 2911, 4659], "ao_hyp": [], "hyp": []} schroeder = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for ethanol of Schroeder " "(2011).", "__doi__": { "autor": "Schroeder, J.A.; Penoncello, S.G.; Schroeder, J.S.", "title": "A Fundamental Equation of State for Ethanol", "ref": "J. Phys. Chem. Ref. Data 43(4) (2014) 043102", "doi": "10.1063/1.4895394"}, # The paper report a diferent value for R, but the program verification # table work with this ancient value "R": 8.314472, "cp": Fi1, "ref": "IIR", "Tmin": 159.0, "Tmax": 650.0, "Pmax": 280000.0, "rhomax": 19.74, "Pmin": 0.00000088, "rhomin": 19.731, "nr1": [0.58200796e-1, 0.94391227, -0.80941908, 0.55359038, -0.14269032e1, 0.13448717], "d1": [4, 1, 1, 2, 2, 3], "t1": [1, 1.04, 2.72, 1.174, 1.329, 0.195], "nr2": [0.42671978, -0.11700261e1, -0.92405872, 0.34891808, -0.9132772, 0.22629481e-1, -0.15513423, 0.21055146, -0.2199769, -0.65857238e-2], "d2": [1, 1, 1, 3, 3, 2, 2, 6, 6, 8], "t2": [2.43, 1.274, 4.16, 3.3, 4.177, 2.5, 0.81, 2.02, 1.606, 0.86], "c2": [1, 1, 2, 1, 2, 1, 2, 1, 1, 1], "gamma2": [1]*10, "nr3": [.75564749, .1069411, -.69533844e-1, -.24947395, .27177891e-1, -0.9053953e-3, -0.12310953, -0.8977971e-1, -0.39512601], "d3": [1, 1, 2, 3, 3, 2, 2, 2, 1], "t3": [2.5, 3.72, 1.19, 3.25, 3, 2, 2, 1, 1], "alfa3": [1.075, .463, .876, 1.108, .741, 4.032, 2.453, 2.3, 3.143], "beta3": [1.207, .0895, .581, .947, 2.356, 27.01, 4.542, 1.287, 3.09], "gamma3": [1.194, 1.986, 1.583, 0.756, 0.495, 1.002, 1.077, 1.493, 1.542], "epsilon3": [.779, .805, 1.869, .694, 1.312, 2.054, .441, .793, .313], "nr4": []} dillon = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for ethanol of Dillon and " "Penoncello (2004)", "__doi__": {"autor": "Dillon, H.E., Penoncello, S.G.", "title": "A Fundamental Equation for Calculation of the " "Thermodynamic Properties of Ethanol", "ref": "Int. J. Thermophys., 25(2) (2004) 321-335", "doi": "10.1023/B:IJOT.0000028470.49774.14"}, "R": 8.314472, "cp": CP1, "ref": {"Tref": 273.15, "Pref": 1., "ho": 45800, "so": 180}, "Tc": 513.9, "rhoc": 5.991, "Tmin": 250.0, "Tmax": 650.0, "Pmax": 280000.0, "rhomax": 19.4, "Pmin": 0.00000088, "rhomin": 19.4, "nr1": [0.114008942201e2, -0.395227128302e2, 0.413063408370e2, -0.188892923721e2, 0.472310314140e1, -0.778322827052e-2, 0.171707850032, -0.153758307602e1, 0.142405508571e1, 0.132732097050, -0.114231649761, 0.327686088736e-5, 0.495699527725e-3, -0.701090149558e-4, -0.225019381648e-5], "d1": [1, 1, 1, 1, 1, 1, 2, 2, 2, 3, 3, 6, 7, 8, 8], "t1": [-0.5, 0, 0.5, 1.5, 2, 5, -0.5, 1, 2, 0, 2.5, 6, 2, 2, 4], "nr2": [-0.255406026981, -0.632036870646e-1, -0.314882729522e-1, 0.256187828185e-1, -0.308694499382e-1, 0.722046283076e-2, 0.299286406225e-2, 0.972795913095e-3], "d2": [1, 3, 3, 6, 7, 8, 2, 7], "t2": [5, 3, 7, 5.5, 4, 1, 22, 23], "c2": [2, 2, 2, 2, 2, 2, 4, 4], "gamma2": [1]*8} sun = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for ethanol 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.314472, "cp": CP1, "ref": {"name": "CUSTOM", "Tref": 273.15, "Pref": 1., "ho": 45800, "so": 180}, "Tmin": Tt, "Tmax": 650.0, "Pmax": 280000.0, "rhomax": 19.6, "Pmin": 0.00000064, "rhomin": 19.55, "nr1": [-2.95455387, 1.95055493, -1.31612955, -1.47547651e-2, 1.39251945e-4, 5.04178939e-1], "d1": [1, 1, 1, 3, 7, 2], "t1": [1.5, 0.25, 1.25, 0.25, 0.875, 1.375], "nr2": [2.52310166e-1, 1.97074652, 8.73146115e-1, 4.27767205e-2, 9.68966545e-2, -8.39632113e-1, -7.71828521e-2, 1.63430744e-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 = schroeder, dillon, sun _PR = 0.0043733 _surface = {"sigma": [0.05], "exp": [0.952]} _vapor_Pressure = { "eq": 3, "n": [-0.91043e1, 0.47263e1, -0.97145e1, 0.41536e1, -0.20739e1], "t": [1, 1.5, 2.0, 2.55, 4.0]} _liquid_Density = { "eq": 1, "n": [0.11632e2, -0.21866e3, 0.82694e3, -0.13512e4, 0.10517e4, -0.31809e3], "t": [0.66, 1.5, 1.9, 2.3, 2.7, 3.1]} _vapor_Density = { "eq": 2, "n": [0.22543e1, -0.24734e2, 0.48993e2, -0.41689e2, -0.45104e2, -0.10732e3], "t": [0.18, 0.44, 0.68, 0.95, 4.0, 10.0]} visco0 = {"__name__": "Kiselev (2005)", "__doi__": { "autor": "Kiselev, S. B., Ely, J. F., Abdulagatov, I. M., " "Huber, M. L.", "title": "Generalized SAFT-DFT/DMT Model for the " "Thermodynamic, Interfacial, and Transport " "Properties of Associating Fluids: Application for " "n-Alkanols", "ref": "Ind. Eng. Chem. Res. 44(17) (2005) 6916-6927", "doi": "10.1021/ie050010e"}, "eq": 1, "omega": 0, "ek": 362.6, "sigma": 0.453, "no": [-1.03116, 3.48379e-2, -6.50264e-6], "to": [0, 1, 2], "Tref_virial": 362.6, "n_virial": [-19.572881, 219.73999, -1015.3226, 2471.0125, -3375.1717, 2491.6597, -787.26086, 14.085455, -0.34664158], "t_virial": [0, -0.25, -0.5, -0.75, -1, -1.25, -1.5, -2.5, -5.5], "Tref_res": 513.9, "rhoref_res": 5.991*M, "muref_res": 1e3, "nr": [0.131194057, -0.0805700894, -0.382240694, 0.153811778, 0.0, -0.110578307], "tr": [0, 0, 1, 1, 2, 2], "dr": [2, 3, 2, 3, 2, 3], "CPf": 23.7222995e3, "CPg1": -3.38264465, "CPgi": [12.7568864/-3.38264465], "CPti": [-0.5]} _viscosity = visco0, thermo0 = {"__name__": "Assael (2013)", "__doi__": { "autor": "Assael, M.J., Sykioti, E.A., Huber, M.L., " "Perkins, R.A.", "title": "Reference Correlation of the Thermal Conductivity" " of Ethanol from the Triple Point to 600 K and " "up to 245 MPa", "ref": "J. Phys. Chem. Ref. Data 42(2) (2013) 023102", "doi": "10.1063/1.4797368"}, "eq": 1, "Toref": 514.71, "koref": 1e-3, "no_num": [-2.09575, 19.9045, -53.964, 82.1223, -1.98864, -0.495513], "to_num": [0, 1, 2, 3, 4, 5], "no_den": [0.17223, -0.078273, 1.0], "to_den": [0, 1, 2], "Tref_res": 514.71, "rhoref_res": 273.186, "kref_res": 1., "nr": [.267222E-01, .148279, -.130429, .346232E-01, -.244293E-02, .0, .177166E-01, -.893088E-01, .684664E-01, -.145702E-01, .809189E-03, .0], "tr": [0, 0, 0, 0, 0, 0, -1, -1, -1, -1, -1, -1], "dr": [1, 2, 3, 4, 5, 6, 1, 2, 3, 4, 5, 6], "critical": 3, "gnu": 0.63, "gamma": 1.239, "R0": 1.02, "Xio": 0.164296e-9, "gam0": 0.05885, "qd": 0.53e-9, "Tcref": 770.85} thermo1 = {"__name__": "Kiselev (2005)", "__doi__": { "autor": "Kiselev, S. B., Ely, J. F., Abdulagatov, I. M., " "Huber, M. L.", "title": "Generalized SAFT-DFT/DMT Model for the " "Thermodynamic, Interfacial, and Transport " "Properties of Associating Fluids: Application for" " n-Alkanols", "ref": "Ind. Eng. Chem. Res. 44(17) (2005) 6916-6927", "doi": "10.1021/ie050010e"}, "eq": 1, "Toref": 1, "koref": 1, "no_num": [-10.109e-3], "to_num": [0.6475], "no_den": [1.0, -7332, -2.68e5], "to_den": [0, -1, -2], "Tref_res": 514.45, "rhoref_res": 5.988*M, "kref_res": 1., "nr": [1.06917458e-1, -5.95897870e-2, -8.65012441e-2, 6.14073818e-2, 2.12220237e-2, -1.00317135e-2], "tr": [0, -1, 0, -1, 0, -1], "dr": [1, 1, 2, 2, 3, 3], # TODO: Add critical crossover model from paper "critical": 0} _thermal = thermo0, thermo1
class Ne(MEoS): """Multiparameter equation of state for neon""" name = "neon" CASNumber = "7440-01-9" formula = "Ne" synonym = "R-720" _refPropName = "NEON" _coolPropName = "Neon" 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} katti = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for neon of Katti (1986)", "__doi__": { "autor": "Katti, R.S., Jacobsen, R.T, Stewart, R.B., " "Jahangiri, M.", "title": "Thermodynamic Properties of 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, "nr1": [ 0.3532653449e1, -0.4513954384e1, -0.1524027959, 0.2188568609e1, -7.44299997, 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.04652799333, 0.04524001818, -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 = katti, _surface = { "sigma": [0.012254, 0.02728, -0.025715], "exp": [1.4136, 1.4517, 1.6567] } _dielectric = { "eq": 1, "a": [0.9969, 0], "b": [-0.109, 0.0708], "c": [-2.88, -1.0], "Au": 0, "D": 2 } _melting = { "eq": 2, "__doi__": { "autor": "Santamaría-Pérez, D., Mukherjee, G.D., Schwager, B., " "Boehler, R.", "title": "High-pressure melting curve of helium and neon: " "Deviations from corresponding states theory", "ref": "Physical Review B 81 (2010) 214101", "doi": "10.1103/PhysRevB.81.214101" }, "Tmin": 24.4, "Tmax": 700.0, "Tref": 24.4, "Pref": 101325, "a2": [0.17e9], "exp2": [1 / 0.77] } _vapor_Pressure = { "eq": 3, "n": [-0.55805e1, 0.68795e-1, 0.54840e1, -0.83760e1, 0.34276e1], "t": [1, 1.5, 2.3, 2.8, 3.4] } _liquid_Density = { "eq": 1, "n": [1.0601, 120.76, -385.53, 816.55, -899.07, 354.66], "t": [0.33, 1.4, 1.7, 2.2, 2.6, 3.0] } _vapor_Density = { "eq": 2, "n": [-0.23338e1, -0.36834e1, -0.85368e2, 0.22769e3, -0.17290e3], "t": [0.444, 0.95, 3.5, 4.1, 4.5] } visco0 = { "__name__": "Rabinovich (1988)", "__doi__": { "autor": "Rabinovich, V.A., Vasserman, A.A., Nedostup, V.I.," " Veksler, L.S.", "title": "Thermophysical Properties of Neon, Argon, " "Krypton, and Xenon", "ref": "Hemisphere Publishing Corp., 1988.", "doi": "" }, "eq": 0, "method": "_visco0" } _viscosity = visco0, def _visco0(self, rho, T, fase=None): a = [ 17.67484, -2.78751, 311498.7, -48826500, 3938774000, -1.654629e11, 2.86561e12 ] Tr = 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 * (T * self.M)**0.5 / y om = 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] sum1 = sum([ci * om**i for i, ci in enumerate(c)]) sum2 = sum([bi * om**i for i, bi in enumerate(b)]) sigma = 3.05e-10 * (sum1 - sum2 * log10(T / 122.1)) br = 2.0 / 3.0 * pi * Avogadro * sigma**3 brho = 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 SO2(MEoS): """Multiparameter equation of state for sulfur dioxide""" name = "sulfur dioxide" CASNumber = "7446-09-5" formula = "SO2" synonym = "R-764" _refPropName = "SO2" _coolPropName = "SulfurDioxide" rhoc = unidades.Density(525.002841) Tc = unidades.Temperature(430.64) Pc = unidades.Pressure(7886.6, "kPa") M = 64.0638 # g/mol Tt = unidades.Temperature(197.7) Tb = unidades.Temperature(263.13) f_acent = 0.256 momentoDipolar = unidades.DipoleMoment(1.6, "Debye") id = 51 Fi1 = { "ao_log": [1, 3.], "pow": [0, 1, -1], "ao_pow": [-4.5414235721, 4.4732289572, -0.0159272204], "ao_exp": [1.0875, 1.916], "titao": [783 / Tc, 1864 / Tc] } Fi2 = { "ao_log": [1, 3.], "pow": [0, 1, -1], "ao_pow": [-4.5328346436, 4.4777967379, -0.01560057996], "ao_exp": [1.062, 1.9401], "titao": [775 / Tc, 1851 / Tc] } CP2 = { "ao": 0.4021066 / 8.3143 * 64.066, "an": [ 0.87348570e-3 / 8.3143 * 64.066, -0.45968820e-6 / 8.3143 * 64.066, -0.13328400e-11 / 8.3143 * 64.066, 0.23785000e-13 / 8.3143 * 64.066 ], "pow": [1, 2, 3, 4], "ao_exp": [], "exp": [] } gao = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for sulfur dioxide of Gao " "(2016).", "__doi__": { "autor": "Gao, K., Wu, J., Zhang, P., Lemmon, E.W.", "title": "A Helmholtz Energy Equation of State for Sulfur Dioxide", "ref": "J. Chem. Eng. Data, 61(6) (2016) 2859-2872", "doi": "10.1021/acs.jced.6b00195" }, "R": 8.3144621, "cp": Fi1, "ref": "NBP", "rhoc": 8.078, "Tmin": Tt, "Tmax": 525.0, "Pmax": 35000.0, "rhomax": 25.40, "nr1": [0.01744413, 1.814878, -2.246338, -0.4602906, 0.1097049], "d1": [4, 1, 1, 2, 3], "t1": [1, 0.45, 0.9994, 1, 0.45], "nr2": [-0.9485769, -0.8751541, 0.4228777, -0.4174962, -0.002903451], "d2": [1, 3, 2, 2, 7], "t2": [2.907, 2.992, 0.87, 3.302, 1.002], "c2": [2, 2, 1, 2, 1], "gamma2": [1] * 5, "nr3": [1.64041, -0.4103535, -0.08316597, -0.2728126, -0.1075782, -0.4348434], "d3": [1, 1, 3, 2, 2, 1], "t3": [1.15, 0.997, 1.36, 2.086, 0.855, 0.785], "alfa3": [1.061, 0.945, 1.741, 1.139, 1.644, 0.647], "beta3": [0.967, 2.538, 2.758, 1.062, 1.039, 0.41], "gamma3": [1.276, 0.738, 0.71, 0.997, 1.35, 0.919], "epsilon3": [0.832, 0.69, 0.35, 0.961, 0.981, 0.333] } lemmon = { "__type__": "Helmholtz", "__name__": "short Helmholtz equation of state for sulfur dioxide 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": Fi2, "ref": "NBP", "Tmin": Tt, "Tmax": 525.0, "Pmax": 35000.0, "rhomax": 25.30, "nr1": [0.93061, -1.9528, -0.17467, 0.061524, 0.00017711], "d1": [ 1, 1, 1, 3, 7, ], "t1": [0.25, 1.25, 1.5, 0.25, 0.875], "nr2": [ 0.21615, 0.51353, 0.010419, -0.25286, -0.054720, -0.059856, -0.016523 ], "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 } polt = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for sulfur dioxide of Polt " "(1987).", "__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": Tt, "Tmax": 523.0, "Pmax": 32000.0, "rhomax": 22.91, "nr1": [ 0.789407019882, -0.170449580056e1, 0.115984637964e1, -0.576307837294, 0.249237283833e1, -0.518115678632e1, 0.320766081899e1, -0.123636065893e1, 0.144419600938e-1, -.15380705504, 0.386324300525, 0.292550313202, -0.372445361392, -0.636924333910e-1, 0.986166451596e-1, -0.216993783055e-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.789407019882, 0.170449580056e1, -0.115984637964e1, -0.480876182378, 0.164910076886e1, -0.133861069604e1 ], "d2": [0, 0, 0, 2, 2, 2], "t2": [3, 4, 5, 3, 4, 5], "c2": [2] * 6, "gamma2": [1] * 6 } eq = gao, lemmon, polt _PR = [-0.068, -18.12] _surface = { "sigma": [0.0803, 0.0139, -0.0114], "exp": [0.928, 1.57, 0.364] } _vapor_Pressure = { "eq": 3, "n": [-7.303, 1.9794, -2.078, -3.5446, 0.51776], "t": [1, 1.5, 2.2, 4.7, 6] } _liquid_Density = { "eq": 1, "n": [7.2296, -16.928, 29.832, -27.901, 11.085], "t": [0.54, 0.88, 1.23, 1.6, 2] } _vapor_Density = { "eq": 2, "n": [-7.487, 10.118, -13.608, -25.408, -42.04, -38.668], "t": [0.545, 0.85, 1.2, 3.7, 7.5, 10] }
class MD2M(MEoS): """Multiparameter equation of state for decamethyltetrasiloxane""" name = "decamethyltetrasiloxane" CASNumber = "141-62-8" formula = "C10H30Si4O3" synonym = "MD2M" rhoc = unidades.Density(284.1716396703609) Tc = unidades.Temperature(599.40) Pc = unidades.Pressure(1227.0, "kPa") M = 310.685 # g/mol Tt = unidades.Temperature(205.2) Tb = unidades.Temperature(467.51) f_acent = 0.668 momentoDipolar = unidades.DipoleMoment(1.12, "Debye") id = 39 # id = 1837 CP1 = {"ao": 331.9, "an": [], "pow": [], "ao_exp": [], "exp": [], "ao_hyp": [329620742.8, 0, 2556558319.0, 0], "hyp": [795.1, 0, 1813.8, 0]} helmholtz1 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for MD2M of Colonna et al. (2006).", "__doi__": {"autor": "Colonna, P., Nannan, N.R., and Guardone, A.", "title": "Multiparameter equations of state for siloxanes: [(CH3)3-Si-O1/2]2-[O-Si-(CH3)2]i=1,…,3, and [O-Si-(CH3)2]6", "ref": "Fluid Phase Equilibria 263:115-130, 2008", "doi": "10.1016/j.fluid.2007.10.001"}, "R": 8.314472, "cp": CP1, "ref": "NBP", "Tmin": Tt, "Tmax": 673.0, "Pmax": 30000.0, "rhomax": 3.033, "Pmin": 0.0000005, "rhomin": 3.032, "nr1": [1.33840331, -2.62939393, 0.4398383, -0.53496715, 0.1818844, 0.40774609e-3], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875], "nr2": [1.13444506, 0.5774631e-1, -0.5917498, -0.11020225, -0.34942635e-1, 0.7646298e-2], "d2": [2, 5, 1, 4, 3, 4], "t2": [0.625, 1.75, 3.625, 3.625, 14.5, 12.0], "c2": [1, 1, 2, 2, 3, 3], "gamma2": [1]*6} eq = helmholtz1, _vapor_Pressure = { "eq": 5, "ao": [-0.10029e2, 0.44434e1, -0.36753e1, -0.68925e1, -0.32211e1], "exp": [1.0, 1.5, 2.06, 3.5, 10.0]} _liquid_Density = { "eq": 1, "ao": [0.12608e2, -0.32120e2, 0.33559e2, -0.11695e2, 0.76192], "exp": [0.48, 0.64, 0.8, 1.0, 2.6]} _vapor_Density = { "eq": 3, "ao": [-0.24684e1, -0.71262e1, -0.27601e2, -0.49458e2, -0.24106e2, -0.19370e3], "exp": [0.376, 0.94, 2.9, 5.9, 6.2, 13.0]}
class DEE(MEoS): """Multiparameter equation of state for diethyl ether""" name = "diethyl ether" CASNumber = "60-29-7" formula = "C4H10O" synonym = "" _refPropName = "DEE" _coolPropName = "DiethylEther" rhoc = unidades.Density(264) Tc = unidades.Temperature(466.7) Pc = unidades.Pressure(3720.238, "kPa") M = 74.1216 # g/mol Tt = unidades.Temperature(156.92) Tb = unidades.Temperature(307.604) f_acent = 0.282 momentoDipolar = unidades.DipoleMoment(1.151, "Debye") id = 162 Fi1 = {"ao_log": [1, 3.36281], "pow": [0, 1, -1, -2, -3], "ao_pow": [17.099494, -6.160844, -8.943822, 0.54621, -0.016604], "ao_exp": [], "titao": []} thol = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for diethyl ether of Thol " "et al. (2013)", "__doi__": {"autor": "Thol, M., Piazza, L., and Span, R.", "title": "A New Functional Form for Equations of State " "for Some Weakly Associating Fluids", "ref": "Int. J. Thermophys., 35(5):783-811, 2014.", "doi": "10.1007/s10765-014-1633-1"}, "R": 8.314472, "cp": Fi1, "ref": "OTO", "Tmin": 270.0, "Tmax": 500.0, "Pmax": 40000.0, "rhomax": 10.6851, "Pmin": 0.000001, "rhomin": 10.6851, "nr1": [0.376700499, -0.116630334, -0.73801498, -0.2725701, -0.04979231, 0.172267029, 0.0044161891], "d1": [1, 1, 1, 2, 3, 3, 5], "t1": [-0.75, -0.25, 1.25, 0.75, -1.0, -0.375, 1.25], "nr2": [-1.53951612, 1.15606052, -0.0184504019, -0.101800599, -0.403598704, 0.00213055571, -0.154741976, 0.0120950552, -0.0143106371], "d2": [1, 1, 2, 5, 1, 3, 4, 5, 2], "t2": [2.375, 3.0, 2.625, 1.875, 4.5, 5.75, 5.375, 2.75, 14.5], "c2": [1, 1, 1, 1, 2, 2, 2, 2, 3], "gamma2": [1]*9} eq = thol, _vapor_Pressure = { "eq": 3, "n": [-7.3059, 1.1734, 0.7142, -4.3219], "t": [1.0, 1.5, 2.2, 3.0]} _liquid_Density = { "eq": 1, "n": [0.3275, 3.1842, -2.1407, 1.4376], "t": [0.12, 0.55, 1.0, 1.4]} _vapor_Density = { "eq": 2, "n": [-0.35858, -16.843, 32.476, -33.444, -48.036], "t": [0.06, 0.87, 1.3, 1.7, 5.3]}
class NH3(MEoS): """Multiparameter equation of state for ammonia""" name = "ammonia" CASNumber = "7664-41-7" formula = "NH3" synonym = "R-717" 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 = {"ao_log": [1, -1], "pow": [0, 1, 1./3, -1.5, -1.75], "ao_pow": [-15.81502, 4.255726, 11.47434, -1.296211, 0.5706757], "ao_exp": [], "titao": [], "ao_hyp": [], "hyp": []} CP2 = {"ao": 5.111814, "an": [-0.42966650e2, -0.10243792e-1, 0.38750775e-4, -0.46406097e-7, 0.20268561e-10], "pow": [-1.001, 1, 2, 3, 4], "ao_exp": [], "exp": [], "ao_hyp": [], "hyp": []} helmholtz1 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for ammonia of Baehr and Tillner-Roth (1993)", "__doi__": {"autor": "Baehr, H.D. and 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": ""}, "__test__": # Table, Pag 42 """ >>> st=NH3(T=-77.65+273.15, x=0.5) >>> print "%0.2f %0.5f %0.5g %0.3g %0.5g %0.5g %0.5g %0.4g %0.5g %0.5g" % (\ st.T.C, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, \ st.Hvap.kJkg, st.Gas.h.kJkg, st.Liquido.s.kJkgK, st.Svap.kJkgK, st.Gas.s.kJkgK) -77.65 0.00609 732.9 0.0641 -143.14 1484.4 1341.2 -0.4715 7.5928 7.1213 >>> st=NH3(T=-50+273.15, x=0.5) >>> print "%0.0f %0.5f %0.5g %0.4g %0.4g %0.5g %0.5g %0.3g %0.5g %0.5g" % (\ st.T.C, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, \ st.Hvap.kJkg, st.Gas.h.kJkg, st.Liquido.s.kJkgK, st.Svap.kJkgK, st.Gas.s.kJkgK) -50 0.04084 702.09 0.3806 -24.73 1415.9 1391.2 0.0945 6.3451 6.4396 >>> st=NH3(T=-25+273.15, x=0.5) >>> print "%0.0f %0.5f %0.5g %0.5g %0.4g %0.5g %0.5g %0.4g %0.5g %0.5g" % (\ st.T.C, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, \ st.Hvap.kJkg, st.Gas.h.kJkg, st.Liquido.s.kJkgK, st.Svap.kJkgK, st.Gas.s.kJkgK) -25 0.15147 671.53 1.2959 86.01 1344.6 1430.7 0.5641 5.4187 5.9827 >>> st=NH3(T=273.15, x=0.5) >>> print "%0.0f %0.5f %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\ st.T.C, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, \ st.Hvap.kJkg, st.Gas.h.kJkg, st.Liquido.s.kJkgK, st.Svap.kJkgK, st.Gas.s.kJkgK) 0 0.42938 638.57 3.4567 200 1262.2 1462.2 1 4.621 5.621 >>> st=NH3(T=25+273.15, x=0.5) >>> print "%0.0f %0.5f %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\ st.T.C, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, \ st.Hvap.kJkg, st.Gas.h.kJkg, st.Liquido.s.kJkgK, st.Svap.kJkgK, st.Gas.s.kJkgK) 25 1.00324 602.76 7.8069 317.67 1165.8 1483.4 1.4089 3.91 5.3188 >>> st=NH3(T=50+273.15, x=0.5) >>> print "%0.0f %0.5f %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\ st.T.C, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, \ st.Hvap.kJkg, st.Gas.h.kJkg, st.Liquido.s.kJkgK, st.Svap.kJkgK, st.Gas.s.kJkgK) 50 2.03403 562.86 15.785 440.62 1050.5 1491.1 1.799 3.2507 5.0497 >>> st=NH3(T=75+273.15, x=0.5) >>> print "%0.0f %0.5f %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\ st.T.C, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, \ st.Hvap.kJkg, st.Gas.h.kJkg, st.Liquido.s.kJkgK, st.Svap.kJkgK, st.Gas.s.kJkgK) 75 3.71045 516.23 29.923 572.38 907.35 1479.7 2.1823 2.6062 4.7885 >>> st=NH3(T=100+273.15, x=0.5) >>> print "%0.0f %0.5f %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\ st.T.C, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, \ st.Hvap.kJkg, st.Gas.h.kJkg, st.Liquido.s.kJkgK, st.Svap.kJkgK, st.Gas.s.kJkgK) 100 6.25527 456.63 56.117 721 715.63 1436.6 2.5797 1.9178 4.4975 >>> st=NH3(T=125+273.15, x=0.5) >>> print "%0.0f %0.5f %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.4g %0.5g" % (\ st.T.C, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, \ st.Hvap.kJkg, st.Gas.h.kJkg, st.Liquido.s.kJkgK, st.Svap.kJkgK, st.Gas.s.kJkgK) 125 9.97022 357.8 120.73 919.68 389.44 1309.1 3.0702 0.9781 4.0483 >>> st=NH3(T=132+273.15, x=0.5) >>> print "%0.0f %0.5f %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.4g %0.5g" % (\ st.T.C, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, \ st.Hvap.kJkg, st.Gas.h.kJkg, st.Liquido.s.kJkgK, st.Svap.kJkgK, st.Gas.s.kJkgK) 132 11.28976 262.7 193.88 1063 108.59 1171.6 3.416 0.268 3.684 """ # Table, Pag 51 """ >>> st=NH3(T=-75+273.15, P=1e4) >>> print "%0.0f %0.2f %0.5g %0.4g" % (st.T.C, st.rho, st.h.kJkg, st.s.kJkgK) -75 730.10 -131.97 -0.4148 >>> st=NH3(T=-30+273.15, P=2e4) >>> print "%0.0f %0.4f %0.5g %0.5g" % (st.T.C, st.rho, st.h.kJkg, st.s.kJkgK) -30 0.1693 1436.9 6.9812 >>> st=NH3(T=273.15, P=3e4) >>> print "%0.0f %0.4f %0.5g %0.5g" % (st.T.C, st.rho, st.h.kJkg, st.s.kJkgK) 0 0.2260 1498.3 7.0226 >>> st=NH3(T=175+273.15, P=5e4) >>> print "%0.0f %0.4f %0.5g %0.5g" % (st.T.C, st.rho, st.h.kJkg, st.s.kJkgK) 175 0.2288 1884.8 7.8612 >>> st=NH3(T=-50+273.15, P=1e5) >>> print "%0.0f %0.2f %0.4g %0.3g" % (st.T.C, st.rho, st.h.kJkg, st.s.kJkgK) -50 702.11 -24.67 0.0944 >>> st=NH3(T=-30+273.15, P=1e5) >>> print "%0.0f %0.4f %0.5g %0.5g" % (st.T.C, st.rho, st.h.kJkg, st.s.kJkgK) -30 0.8641 1426.1 6.1607 >>> st=NH3(T=273.15, P=2e5) >>> print "%0.0f %0.4f %0.5g %0.5g" % (st.T.C, st.rho, st.h.kJkg, st.s.kJkgK) 0 1.5468 1483.8 6.0557 >>> st=NH3(T=125+273.15, P=3e5) >>> print "%0.0f %0.5g %0.5g %0.5g" % (st.T.C, st.rho, st.h.kJkg, st.s.kJkgK) 125 1.5597 1762.6 6.7009 >>> st=NH3(T=-50+273.15, P=5e5) >>> print "%0.0f %0.5g %0.4g %0.3g" % (st.T.C, st.rho, st.h.kJkg, st.s.kJkgK) -50 702.28 -24.32 0.0934 >>> st=NH3(T=273.15, P=1e6) >>> print "%0.0f %0.5g %0.5g %0.4g" % (st.T.C, st.rho, st.h.kJkg, st.s.kJkgK) 0 638.97 200.37 0.9981 >>> st=NH3(T=273.15, P=2e6) >>> print "%0.0f %0.5g %0.5g %0.4g" % (st.T.C, st.rho, st.h.kJkg, st.s.kJkgK) 0 639.67 201.01 0.9947 >>> st=NH3(T=-50+273.15, P=3e6) >>> print "%0.0f %0.5g %0.4g %0.3g" % (st.T.C, st.rho, st.h.kJkg, st.s.kJkgK) -50 703.33 -22.08 0.0875 """, "R": 8.314471, "cp": Fi1, "ref": "IIR", "Tmin": Tt, "Tmax": 700., "Pmax": 1000000.0, "rhomax": 52.915, "Pmin": 6.09, "rhomin": 43.035, "nr1": [-0.1858814e01, 0.4554431e-1, 0.7238548, 0.1229470e-1, 0.2141882e-10], "d1": [1, 2, 1, 4, 15], "t1": [1.5, -0.5, 0.5, 1., 3.], "nr2": [-0.1430020e-1, 0.3441324, -0.2873571, 0.2352589e-4, -0.3497111e-1, 0.1831117e-2, 0.2397852e-1, -0.4085375e-1, 0.2379275, -0.3548972e-1, -0.1823729, 0.2281556e-1, -0.6663444e-2, -0.8847486e-2, 0.2272635e-2, -0.5588655e-3], "d2": [3, 3, 1, 8, 2, 8, 1, 1, 2, 3, 2, 4, 3, 1, 2, 4], "t2": [0, 3, 4, 4, 5, 5, 3, 6, 8, 8, 10, 10, 5, 7.5, 15, 30], "c2": [1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3], "gamma2": [1]*16} helmholtz2 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for ammonia of Ahrendts and Baehr (1979)", "__doi__": {"autor": "Ahrendts, J. and Baehr, H.D.", "title": "The Thermodynamic Properties of Ammonia", "ref": "VDI-Forsch., Number 596, 1979.", "doi": ""}, "R": 8.31434, "cp": CP2, "ref": "IIR", "Tmin": 195.486, "Tmax": 600., "Pmax": 400000.0, "rhomax": 44.0, "Pmin": 6.0339, "rhomin": 43.137, "nr1": [0.911447599671, -0.382129415537e1, 0.147730246416e1, 0.580205129871e-1, -0.574413226616e-3, 0.153018094697, -0.256626062036, 0.445448838055, -0.1533210545, 0.527996725202e-1, -0.484726581121e-1, 0.246579503330e-2, -0.107999941003e-3, -0.215298673010e-4, -0.306938893790e-4, 0.839163613582e-5, 0.814833533876e-6, -0.314753664228e-7], "d1": [1, 1, 1, 1, 1, 2, 2, 3, 4, 5, 5, 7, 9, 9, 10, 11, 12, 14], "t1": [1, 2, 3, 6, 9, 0, 4, 2, 1, 1, 2, 3, 3, 5, 1, 1, 5, 5], "nr2": [0.642978802435, -0.139510669941e1, 0.956135683432, -0.272787386366, -0.189305337334e1, 0.479043603913e1, -0.245945016980e1, -0.121107723958e1, 0.500552271170e1, -0.615476024667e1, 0.210772481535e1, 0.298003513465, -0.152506723279, 0.115565883925e-2, -0.911244657201e-3, 0.100587210000e-1, -0.120983155888e-1, 0.382694351151e-2], "d2": [1, 1, 1, 1, 2, 2, 2, 3, 3, 3, 3, 4, 4, 0, 0, 0, 0, 0], "t2": [2, 5, 6, 7, 5, 6, 7, 3, 4, 5, 6, 6, 7, 1, 2, 0, 1, 2], "c2": [2]*18, "gamma2": [0.86065403]*13+[506.2670781840292]*2+[50626.70781840292]*3} helmholtz3 = { "__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):111-162, 2003.", "doi": "10.1023/A:1022362231796"}, "__test__": """ >>> st=NH3(T=700, rho=200, eq=2) >>> print "%0.4f %0.3f %0.4f" % (st.cp0.kJkgK, st.P.MPa, st.cp.kJkgK) 2.4759 136.271 4.1916 >>> st2=NH3(T=750, rho=100, eq=2) >>> print "%0.2f %0.5f" % (st2.h.kJkg-st.h.kJkg, st2.s.kJkgK-st.s.kJkgK) 776.69 2.07036 """, # Table III, Pag 117 "R": 8.314471, "cp": Fi1, "ref": "IIR", "Tmin": Tt, "Tmax": 600., "Pmax": 100000.0, "rhomax": 52.915, "Pmin": 6.0531, "rhomin": 43.158, "nr1": [0.7302272, -0.11879116e1, -0.68319136, 0.40028683e-1, 0.90801215e-4], "d1": [1, 1, 1, 3, 7], "t1": [3, 4, 5, 0, 1, 2, 3, 4, 0, 1, 2, 0, 1, 0, 1, 1], "nr2": [-0.56216175e-1, 0.44935601, 0.29897121e-1, -0.18181684, -0.9841666e-1, -0.55083744e-1, -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} helmholtz4 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for ammonia of Sun and Ely (2004)", "__doi__": {"autor": "Sun, L. and Ely, J.F.", "title": "Universal equation of state for engineering application: Algorithm and application to non-polar and polar fluids", "ref": "Fluid Phase Equilib., 222-223:107-118, 2004.", "doi": "10.1016/j.fluid.2004.06.028"}, "R": 8.3143, "cp": Fi1, "ref": "IIR", "Tmin": Tt, "Tmax": 620.0, "Pmax": 800000.0, "rhomax": 40., "Pmin": 0.1, "rhomin": 40., "nr1": [3.29159441e-1, 8.48237019e-1, -2.30706412, 4.08625188e-2, 6.79597481e-5, 4.99412149e-2], "d1": [1, 1, 1, 3, 7, 2], "t1": [1.5, 0.25, 1.25, 0.25, 0.875, 1.375], "nr2": [1.23624654e-1, -3.02129187e-1, 3.31747586e-1, -2.97121254e-3, -1.30202073e-1, -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 = helmholtz1, helmholtz2, helmholtz3, helmholtz4 _melting = {"eq": 1, "Tref": Tt, "Pref": 1000, "Tmin": Tt, "Tmax": 700.0, "a1": [], "exp1": [], "a2": [], "exp2": [], "a3": [0.2533125e4], "exp3": [1]} _surface = {"sigma": [0.1028, -0.09453], "exp": [1.211, 5.585]} _vapor_Pressure = { "eq": 5, "ao": [-0.70993e1, -0.24330e1, 0.87591e1, -0.64091e1, -0.21185e1], "exp": [1., 1.5, 1.7, 1.95, 4.2]} _liquid_Density = { "eq": 1, "ao": [0.34488e2, -0.12849e3, 0.17382e3, -0.10699e3, 0.30339e2], "exp": [0.58, 0.75, 0.9, 1.1, 1.3]} _vapor_Density = { "eq": 3, "ao": [-0.38435, -0.40846e1, -0.66634e1, -0.31881e2, 0.21306e3, -0.24648e3], "exp": [0.218, 0.55, 1.5, 3.7, 5.5, 5.8]} visco0 = {"eq": 1, "omega": 1, "collision": [4.99318220, -0.61122364, 0.0, 0.18535124, -0.11160946], "__name__": "Fenghour (1995)", "__doi__": {"autor": "Fenghour, A., Wakeham, W.A., Vesovic, V., Watson, J.T.R., Millat, J., and Vogel, E.", "title": "The viscosity of ammonia", "ref": "J. Phys. Chem. Ref. Data 24, 1649 (1995)", "doi": "10.1063/1.555961"}, "__test__": # Appendix II, pag 1664 """ >>> st=NH3(T=200, P=1e5) >>> print "%0.2f" % st.mu.muPas 507.47 >>> st=NH3(T=290, P=1e6) >>> print "%0.2f" % st.mu.muPas 142.93 >>> st=NH3(T=250, P=1e7) >>> print "%0.2f" % st.mu.muPas 233.81 >>> st=NH3(T=300, P=1e5) >>> print "%0.2f" % st.mu.muPas 10.16 >>> st=NH3(T=350, P=1.8e7) >>> print "%0.2f" % st.mu.muPas 91.36 >>> st=NH3(T=400, P=5e7) >>> print "%0.2f" % st.mu.muPas 77.29 >>> st=NH3(T=490, P=1e6) >>> print "%0.2f" % st.mu.muPas 17.49 >>> st=NH3(T=550, P=1e5) >>> print "%0.2f" % st.mu.muPas 19.79 >>> st=NH3(T=680, P=5e7) >>> print "%0.2f" % st.mu.muPas 31.90 """ # Appendix III, pag 1667 """ >>> st=NH3(T=196, x=0.5) >>> print "%0.0f %0.4f %0.4f %0.2f %0.4f %0.2f " % ( \ st.T, st.P.MPa, st.Gas.rhoM, st.Gas.mu.muPas, st.Liquido.rhoM, st.Liquido.mu.muPas) 196 0.0063 0.0039 6.85 43.0041 553.31 >>> st=NH3(T=240, x=0.5) >>> print "%0.0f %0.4f %0.4f %0.2f %0.4f %0.2f " % ( \ st.T, st.P.MPa, st.Gas.rhoM, st.Gas.mu.muPas, st.Liquido.rhoM, st.Liquido.mu.muPas) 240 0.1022 0.0527 8.06 40.0318 254.85 >>> st=NH3(T=280, x=0.5) >>> print "%0.0f %0.4f %0.4f %0.2f %0.4f %0.2f " % ( \ st.T, st.P.MPa, st.Gas.rhoM, st.Gas.mu.muPas, st.Liquido.rhoM, st.Liquido.mu.muPas) 280 0.5509 0.2573 9.27 36.9389 158.12 >>> st=NH3(T=300, x=0.5) >>> print "%0.0f %0.4f %0.4f %0.2f %0.4f %0.2f " % ( \ st.T, st.P.MPa, st.Gas.rhoM, st.Gas.mu.muPas, st.Liquido.rhoM, st.Liquido.mu.muPas) 300 1.0617 0.4845 9.89 35.2298 129.33 >>> st=NH3(T=340, x=0.5) >>> print "%0.0f %0.4f %0.4f %0.2f %0.4f %0.2f " % ( \ st.T, st.P.MPa, st.Gas.rhoM, st.Gas.mu.muPas, st.Liquido.rhoM, st.Liquido.mu.muPas) 340 3.0803 1.4325 11.33 31.2641 88.55 >>> st=NH3(T=360, x=0.5) >>> print "%0.0f %0.4f %0.4f %0.2f %0.4f %0.2f " % ( \ st.T, st.P.MPa, st.Gas.rhoM, st.Gas.mu.muPas, st.Liquido.rhoM, st.Liquido.mu.muPas) 360 4.7929 2.3598 12.35 28.7879 65.49 >>> st=NH3(T=380, x=0.5) >>> print "%0.0f %0.4f %0.4f %0.2f %0.4f %0.2f " % ( \ st.T, st.P.MPa, st.Gas.rhoM, st.Gas.mu.muPas, st.Liquido.rhoM, st.Liquido.mu.muPas) 380 7.1403 3.9558 14.02 25.6059 58.31 >>> st=NH3(T=398, x=0.5) >>> print "%0.0f %0.4f %0.4f %0.2f %0.4f %0.2f " % ( \ st.T, st.P.MPa, st.Gas.rhoM, st.Gas.mu.muPas, st.Liquido.rhoM, st.Liquido.mu.muPas) 398 9.9436 7.0447 17.67 21.0667 43.95 >>> st=NH3(T=402, x=0.5) >>> print "%0.0f %0.4f %0.4f %0.2f %0.4f %0.2f " % ( \ st.T, st.P.MPa, st.Gas.rhoM, st.Gas.mu.muPas, st.Liquido.rhoM, st.Liquido.mu.muPas) 402 10.6777 8.5479 19.69 19.0642 39.20 """ , "ek": 386., "sigma": 0.2957, "Tref": 1., "rhoref": 1.*M, "n_chapman": 8.8135503/M**0.5, "n_virial": [-0.17999496e1, 0.46692621e2, -0.53460794e3, 0.33604074e4, -0.13019164e5, 0.33414230e5, -0.58711743e5, 0.71426686e5, -0.59834012e5, 0.33652741e5, -0.1202735e5, 0.24348205e4, -0.20807957e3], "t_virial": [0, -0.5, -1, -1.5, -2, -2.5, -3, -3.5, -4, -4.5, -5, -5.5, -6], "Tref_virial": 386., "etaref_virial": 0.015570557, "Tref_res": 386., "rhoref_res": 1.*M, "etaref_res": 1, "n_poly": [2.19664285e-1, -0.83651107e-1, 0.17366936e-2, -0.64250359e-2, 1.67668649e-4, -1.49710093e-4, 0.77012274e-4], "t_poly": [-2, -4, -0, -1, -2, -3, -4], "d_poly": [2, 2, 3, 3, 4, 4, 4], "g_poly": [0, 0, 0, 0, 0, 0, 0], "c_poly": [0, 0, 0, 0, 0, 0, 0]} _viscosity = visco0, thermo0 = {"eq": 1, "critical": "NH3", "__name__": "Tufeu (1984)", "__doi__": {"autor": "Tufeu, R., Ivanov, D.Y., Garrabos, Y., and 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:422-427, 1984", "doi": "10.1002/bbpc.19840880421"}, "Tref": 1., "kref": 1., "no": [0.3589e-1, -0.1750e-3, 0.4551e-6, 0.1685e-9, -0.4828e-12], "co": [0, 1, 2, 3, 4], "Trefb": 1., "rhorefb": 0.05871901, "krefb": 1., "nb": [0.16207e-3, 0.12038e-5, -0.23139e-8, 0.32749e-11], "tb": [0, 0, 0, 0], "db": [1, 2, 3, 4], "cb": [0, 0, 0, 0]} _thermal = thermo0,
class nC5(MEoS): """Multiparameter equation of state for n-pentane""" name = "pentane" CASNumber = "109-66-0" formula = "CH3-(CH2)3-CH3" synonym = "R-601" _refPropName = "PENTANE" _coolPropName = "n-Pentane" rhoc = unidades.Density(232.) Tc = unidades.Temperature(469.7) Pc = unidades.Pressure(3370.0, "kPa") M = 72.14878 # g/mol Tt = unidades.Temperature(143.47) Tb = unidades.Temperature(309.21) f_acent = 0.251 momentoDipolar = unidades.DipoleMoment(0.07, "Debye") id = 8 _Tr = unidades.Temperature(449.271155) _rhor = unidades.Density(233.873368) _w = 0.247058753 Fi1 = {"ao_log": [1, 3.0], "ao_sinh": [8.95043, 33.4032], "sinh": [0.380391739, 3.777411113], "ao_cosh": [21.836], "cosh": [1.789520971]} CP0 = {"ao": 4, "ao_sinh": [8.95043, 33.4032], "sinh": [178.67, 1774.25], "ao_cosh": [21.836], "cosh": [840.538]} CP1 = {"ao": 10.288132, "an": [-0.2695377e-1, 0.20951065e-3, -0.27910773e-6, 0.12266269e-9], "pow": [1, 2, 3, 4]} f = 4.184/8.3159524 CP2 = {"ao": 22.5012*f, "ao_sinh": [2.057417e8*f], "sinh": [1.71958e3], "ao_cosh": [2.972927e7*f], "cosh": [8.02069e2]} CP3 = {"ao": 4, "ao_exp": [9.751560716, 22.71445741, 11.65392685], "exp": [404.8796661, 1785.491483, 4504.430788]} shortSpan = { "__type__": "Helmholtz", "__name__": "short Helmholtz equation of state for pentane of Span " "and Wagner (2003)", "__doi__": {"autor": "Span, R., Wagner, W.", "title": "Equations of state for technical applications. " "II. Results for nonpolar fluids.", "ref": "Int. J. Thermophys. 24 (1) (2003) 41-109", "doi": "10.1023/A:1022310214958"}, "R": 8.31451, "cp": CP0, "ref": "OTO", "M": 72.15, "Tc": 469.7, "rhoc": 232/72.15, "Tmin": Tt, "Tmax": 750.0, "Pmax": 100000.0, "rhomax": 11.2, "nr1": [0.10968643e1, -0.29988888e1, 0.99516887, -0.16170709, 0.11334460, 0.26760595e-3], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875], "nr2": [0.40979882, -0.40876423e-1, -0.38169482, -0.10931957, -0.32073223e-1, 0.16877016e-1], "d2": [2, 5, 1, 4, 3, 4], "t2": [0.625, 1.75, 3.625, 3.625, 14.5, 12.], "c2": [1, 1, 2, 2, 3, 3], "gamma2": [1]*6} GERG = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for pentane of Kunz and " "Wagner (2004).", "__doi__": {"autor": "Kunz, O., Wagner, W.", "title": "The GERG-2008 Wide-Range Equation of State for " "Natural Gases and Other Mixtures: An Expansion " "of GERG-2004", "ref": "J. Chem.Eng. Data 57(11) (2012) 3032-3091", "doi": "10.1021/je300655b"}, "R": 8.314472, "cp": Fi1, "ref": "OTO", "Tmin": 143.47, "Tmax": 600.0, "Pmax": 100000.0, "rhomax": 10.57, "nr1": [0.10968643098001e1, -0.29988888298061e1, 0.99516886799212, -0.16170708558539, 0.11334460072775, 0.26760595150748e-3], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875], "nr2": [0.40979881986931, -0.40876423083075e-1, -0.38169482469447, -0.10931956843993, -0.32073223327990e-1, 0.16877016216975e-1], "d2": [2, 5, 1, 4, 3, 4], "t2": [0.625, 1.75, 3.625, 3.625, 14.5, 12.], "c2": [1, 1, 2, 2, 3, 3], "gamma2": [1]*6} polt = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for pentane of Polt (1992)", "__doi__": {"autor": "Polt, A., Platzer, B., Maurer, G.", "title": "Parameter der thermischen Zustandsgleichung von " "Bender fuer 14 mehratomige reine Stoffe", "ref": "Chem. Technik 22(1992)6 , 216/224", "doi": ""}, "R": 8.3143, "cp": CP1, "ref": "NBP", "Tmin": 238.0, "Tmax": 573.0, "Pmax": 30000.0, "rhomax": 9.410819, "nr1": [-0.117648900900e1, 0.163499095773e1, -0.366669005817, 0.724947274043, -0.221919300269e1, 0.188671490348e1, -0.195774652096e1, 0.308440851184, 0.437424419722, -0.625853472351, .382868807091, -0.119467393955, .218631441082, 0.485668874195e-1, -0.132198161379, 0.213549844850e-1], "d1": [0, 0, 0, 1, 1, 1, 1, 1, 2, 2, 2, 3, 3, 4, 4, 5], "t1": [3, 4, 5, 0, 1, 2, 3, 4, 0, 1, 2, 0, 1, 0, 1, 1], "nr2": [0.117648900900e1, -0.163499095773e1, 0.366669005817, -0.363660829618e-2, 0.633672105685, -0.705792643982], "d2": [0, 0, 0, 2, 2, 2], "t2": [3, 4, 5, 3, 4, 5], "c2": [2]*6, "gamma2": [0.968832]*6} starling = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for pentane of Starling " "(1973)", "__doi__": {"autor": "Starling, K.E.", "title": "Fluid Thermodynamic Properties for Light " "Petroleum Systems", "ref": "Gulf Publishing Company, 1973.", "doi": ""}, "R": 8.3159524, "cp": CP2, "ref": "NBP", "Tmin": 177.0, "Tmax": 589.0, "Pmax": 55000.0, "rhomax": 10.2534, "nr1": [0.175873733594e1, 0.485604047435, -0.111896446456e1, -0.685918143315, 0.368714111378e-1, -0.167498784887e-2, 0.327765295239, -0.352742092747, -0.999487301826e-1, 0.781999120830e-2, 0.221577806386e-2], "d1": [0, 1, 1, 1, 1, 1, 2, 2, 2, 5, 5], "t1": [3, 0, 1, 3, 4, 5, 0, 1, 2, 1, 2], "nr2": [-0.175873733594e1, -0.411653507564], "d2": [0, 2], "t2": [3, 3], "c2": [2]*2, "gamma2": [0.46812392]*2} sun = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for pentane of Sun and Ely " "(2004)", "__doi__": {"autor": "Sun, L., Ely, J.F.", "title": "Universal equation of state for engineering " "application: Algorithm and application to " "non-polar and polar fluids", "ref": "Fluid Phase Equilib., 222-223 (2004) 107-118", "doi": "10.1016/j.fluid.2004.06.028"}, "R": 8.31451, "cp": Fi1, "ref": "OTO", "Tmin": Tt, "Tmax": 620.0, "Pmax": 800000.0, "rhomax": 40., "nr1": [2.20261753, 1.07797592, -3.82130221, 1.06627357e-1, 3.07513215e-4, -2.84309667e-1], "d1": [1, 1, 1, 3, 7, 2], "t1": [1.5, 0.25, 1.25, 0.25, 0.875, 1.375], "nr2": [-7.28441220e-2, -4.60943732e-1, 8.39360011e-2, -1.50650444e-2, -0.203771872, -7.90244277e-3, -5.68993564e-2, -2.99387974e-2], "d2": [1, 1, 2, 5, 1, 1, 4, 2], "t2": [0, 2.375, 2., 2.125, 3.5, 6.5, 4.75, 12.5], "c2": [1, 1, 1, 1, 2, 2, 2, 3], "gamma2": [1]*8} ratanapisit = { "__type__": "MBWR", "__name__": "MBWR equation of state for pentane of Ratanapisit (1999)", "__doi__": {"autor": "Ratanapisit, J., Ely, J.F.", "title": "Application of New, Modified BWR Equations of " "State to the Corresponding-States Prediction of " "Natural Gas Properties", "ref": "Int. J. Thermophys., 20(6) (1999) 1721-1735", "doi": "10.1023/A:1022610013596"}, "R": 8.31451, "Tc": 469.65, "Pc": 3364.56, "rhoc": 3.2155, "cp": CP3, "ref": "OTO", "Tmin": Tt, "Tmax": 600.0, "Pmax": 70000.0, "rhomax": 11.2, "b": [None, -7.41533782499e-2, 7.54044021950, -1.93328401588e2, 3.39428034054e4, -5.12571561595e6, 1.51195406963e-3, -7.12225059892, 4.12664185793e3, 8.40258305443e5, -4.68416651753e-4, 3.03565637672, -1.42146321204e3, -1.10170659283e-1, -9.80664356304, 1.10979804446e3, 2.98029604130, -1.41484307201e-1, -3.39208006239e1, 2.08782048763, 5.38055429992e5, -6.40401885304e8, -1.19676622034e5, 1.71973349582e10, -3.06383363882e3, 1.43168348944e6, 1.41452433419e1, -2.52955687564e7, -3.85316416299, 2.65416349789e3, 4.76643876980e-3, -8.37595968663, -1.35160880503e3]} eq = shortSpan, GERG, polt, starling, sun, ratanapisit _PR = [-0.0752, -17.5] _surface = {"sigma": [0.08015, 0.004384, -0.03437], "exp": [1.408, 1.031, 1.818]} _dielectric = { "eq": 1, "a": [25.39, 0.025], "b": [78.39, 54.15], "c": [-12480, -4800.0], "Au": 29.84, "D": 2} _melting = { "eq": 2, "__doi__": { "autor": "Reeves, L.E., Scott, G.J., Babb, S.E. Jr.", "title": "Melting Curves of Pressure-Transmitting fluids", "ref": "Fluid Phase Equilib., 222-223 (2004) 107-118", "doi": "10.1063/1.1725068"}, "Tmin": Tt, "Tmax": 2000.0, "Tref": Tt, "Pref": 0.076321, "a2": [6600e5], "exp2": [1.649]} _vapor_Pressure = { "eq": 3, "n": [-0.73918e1, 0.31102e1, -0.22415e1, -0.31585e1, -0.90451], "t": [1., 1.5, 1.74, 3.75, 8.0]} _liquid_Density = { "eq": 1, "n": [0.10178e1, 0.42703, 0.11334e1, 0.41518, -0.47950e-1], "t": [0.27, 0.44, 0.6, 4.0, 5.0]} _vapor_Density = { "eq": 2, "n": [-2.9389, -6.2784, -19.941, -16.709, -36.543, -127.99], "t": [0.4, 1.18, 3.2, 6.6, 7.0, 15.0]} visco0 = {"__name__": u"Quiñones-Cisneros (2006)", "__doi__": { "autor": "Quiñones-Cisneros, S.E., Deiters, U.K.", "title": "Generalization of the Friction Theory for " "Viscosity Modeling", "ref": "J. Phys. Chem. B, 110(25) (2006) 12820-12834", "doi": "10.1021/jp0618577"}, "eq": 4, "omega": 0, "Toref": 469.7, "no": [17.6805, -55.6942, 48.7177], "to": [0, 0.25, 0.5], "a": [1.08193e-5, -4.71699e-5, 0.0], "b": [1.21502e-4, -9.84766e-5, 0.0], "c": [5.08307e-5, -1.07e-5, 0.0], "A": [-2.10025e-10, -1.56583e-9, 0.0], "B": [1.98521e-8, 2.05972e-9, 0.0], "C": [-1.18487e-7, 1.69571e-7, 0.0]} _viscosity = visco0, thermo0 = {"__name__": "Vassiliou (2015)", "__doi__": { "autor": "Vassiliou, C.-M., Assael, M.J., Huber, M.L., " "Perkins, R.A.", "title": "Reference Correlation of the Thermal Conductivity" " of Cyclopentane, iso-pentane, and n-Pentane", "ref": "J. Phys. Chem. Ref. Data 44(3) (2015) 033102", "doi": "10.1063/1.4927095"}, "eq": 1, "Toref": 469.7, "koref": 1e-3, "no_num": [-3.96685, 35.3805, 5.11554, -108.585, 179.573, 39.2128], "to_num": [0, 1, 2, 3, 4, 5], "no_den": [2.71636, -5.76265, 6.77885, -0.59135, 1], "to_den": [0, 1, 2, 3, 4], "Tref_res": 469.7, "rhoref_res": 232, "kref_res": 1e-3, "nr": [7.76054e-1, 1.17655e2, -1.33101e2, 5.34026e1, -6.8793, 7.97696, -7.85888e1, 9.16089e1, -3.70431e1, 5.0962], "tr": [0, 0, 0, 0, 0, -1, -1, -1, -1, -1], "dr": [1, 2, 3, 4, 5, 1, 2, 3, 4, 5], "critical": 3, "gnu": 0.63, "gamma": 1.239, "R0": 1.02, "Xio": 0.227e-9, "gam0": 0.058, "qd": 0.668e-9, "Tcref": 704.55} _thermal = thermo0,
class H2O(MEoS): """Multiparameter equation of state for water (including IAPWS95) # >>> water=H2O(T=300, rho=996.5560) # >>> print("%0.10f %0.8f %0.5f %0.9f" % (water.P.MPa, water.cv.kJkgK, water.w, water.s.kJkgK)) # 0.0992418350 4.13018112 1501.51914 0.393062643 # # >>> water=H2O(T=500, rho=0.435) # >>> print("%0.10f %0.8f %0.5f %0.9f" % (water.P.MPa, water.cv.kJkgK, water.w, water.s.kJkgK)) # 0.0999679423 1.50817541 548.31425 7.944882714 # # >>> water=H2O(T=900., P=700e6) # >>> print("%0.4f %0.8f %0.5f %0.8f" % (water.rho, water.cv.kJkgK, water.w, water.s.kJkgK)) # 870.7690 2.66422350 2019.33608 4.17223802 # # >>> water=H2O(T=300., P=0.1e6) # >>> print("%0.2f %0.5f %0.2f %0.2f %0.5f %0.4f %0.1f %0.6f" % (water.T, water.P.MPa, water.rho, water.h.kJkg, water.s.kJkgK, water.cp.kJkgK, water.w, water.virialB)) # 300.00 0.10000 996.56 112.65 0.39306 4.1806 1501.5 -0.066682 # # >>> water=H2O(T=500., P=0.1e6) # >>> print("%0.2f %0.5f %0.5f %0.1f %0.4f %0.4f %0.2f %0.7f" % (water.T, water.P.MPa, water.rho, water.h.kJkg, water.s.kJkgK, water.cp.kJkgK, water.w, water.virialB)) # 500.00 0.10000 0.43514 2928.6 7.9447 1.9813 548.31 -0.0094137 # # >>> water=H2O(T=450., x=0.5) # >>> print("%0.2f %0.5f %0.4f %0.1f %0.4f %0.6f" % (water.T, water.P.MPa, water.rho, water.h.kJkg, water.s.kJkgK, water.virialB)) # 450.00 0.93220 9.5723 1761.8 4.3589 -0.013028 # # >>> water=H2O(P=1.5e6, rho=1000.) # >>> print("%0.2f %0.4f %0.1f %0.3f %0.5f %0.4f %0.1f %0.6f" % (water.T, water.P.MPa, water.rho, water.h.kJkg, water.s.kJkgK, water.cp.kJkgK, water.w, water.virialB)) # 286.44 1.5000 1000.0 57.253 0.19931 4.1855 1462.1 -0.085566 # # >>> water=H2O(h=3000e3, s=8e3) # >>> print("%0.2f %0.5f %0.5f %0.1f %0.4f %0.4f %0.2f %0.7f" % (water.T, water.P.MPa, water.rho, water.h.kJkg, water.s.kJkgK, water.cp.kJkgK, water.w, water.virialB)) # 536.24 0.11970 0.48547 3000.0 8.0000 1.9984 567.04 -0.0076606 # # >>> water=H2O(h=150e3, s=0.4e3) # >>> print("%0.2f %0.5f %0.2f %0.2f %0.5f %0.4f %0.1f %0.6f" % (water.T, water.P.MPa, water.rho, water.h.kJkg, water.s.kJkgK, water.cp.kJkgK, water.w, water.virialB)) # 301.27 35.50549 1011.48 150.00 0.40000 4.0932 1564.1 -0.065238 # # >>> water=H2O(T=450., rho=300) # >>> print("%0.2f %0.5f %0.2f %0.2f %0.4f %0.6f %0.6f" % (water.T, water.P.MPa, water.rho, water.h.kJkg, water.s.kJkgK, water.x, water.virialB)) # 450.00 0.93220 300.00 770.82 2.1568 0.010693 -0.013028 # # >>> water=H2O(rho=300., P=0.1e6) # >>> print("%0.2f %0.5f %0.2f %0.2f %0.4f %0.7f %0.6f" % (water.T, water.P.MPa, water.rho, water.h.kJkg, water.s.kJkgK, water.x, water.virialB)) # 372.76 0.10000 300.00 420.56 1.3110 0.0013528 -0.025144 # # >>> water=H2O(h=1500e3, P=0.1e6) # >>> print("%0.2f %0.5f %0.4f %0.1f %0.4f %0.5f %0.6f" % (water.T, water.P.MPa, water.rho, water.h.kJkg, water.s.kJkgK, water.x, water.virialB)) # 372.76 0.10000 1.2303 1500.0 4.2068 0.47952 -0.025144 # # >>> water=H2O(s=5e3, P=3.5e6) # >>> print("%0.2f %0.4f %0.3f %0.1f %0.4f %0.5f %0.7f" % (water.T, water.P.MPa, water.rho, water.h.kJkg, water.s.kJkgK, water.x, water.virialB)) # 515.71 3.5000 25.912 2222.8 5.0000 0.66921 -0.0085877 # # >>> water=H2O(T=500., u=900e3) # >>> print("%0.2f %0.2f %0.2f %0.2f %0.1f %0.4f %0.4f %0.1f %0.7f" % (water.T, water.P.MPa, water.rho, water.u.kJkg, water.h.kJkg, water.s.kJkgK, water.cp.kJkgK, water.w, water.virialB)) # 500.00 108.21 903.62 900.00 1019.8 2.4271 4.1751 1576.0 -0.0094137 # # >>> water=H2O(P=0.3e6, u=1550.e3) # >>> print("%0.2f %0.5f %0.4f %0.1f %0.1f %0.4f %0.5f %0.6f" % (water.T, water.P.MPa, water.rho, water.u.kJkg, water.h.kJkg, water.s.kJkgK, water.x, water.virialB)) # 406.67 0.30000 3.3029 1550.0 1640.8 4.3260 0.49893 -0.018263 # # >>> water=H2O(rho=300, h=1000.e3) # >>> print("%0.2f %0.4f %0.2f %0.2f %0.1f %0.4f %0.6f %0.7f" % (water.T, water.P.MPa, water.rho, water.u.kJkg, water.h.kJkg, water.s.kJkgK, water.x, water.virialB)) # 494.92 2.3991 300.00 992.00 1000.0 2.6315 0.026071 -0.0097064 # # >>> water=H2O(rho=30, s=8.e3) # >>> print("%0.2f %0.3f %0.3f %0.1f %0.1f %0.4f %0.4f %0.2f %0.9f" % (water.T, water.P.MPa, water.rho, water.u.kJkg, water.h.kJkg, water.s.kJkgK, water.cp.kJkgK, water.w, water.virialB)) # 1562.42 21.671 30.000 4628.5 5350.9 8.0000 2.7190 943.53 0.000047165 # # >>> water=H2O(rho=30, s=4.e3) # >>> print("%0.2f %0.4f %0.3f %0.1f %0.1f %0.4f %0.5f %0.7f" % (water.T, water.P.MPa, water.rho, water.u.kJkg, water.h.kJkg, water.s.kJkgK, water.x, water.virialB)) # 495.00 2.4029 30.000 1597.3 1677.4 4.0000 0.39218 -0.0097015 # # >>> water=H2O(rho=300, u=1000.e3) # >>> print("%0.2f %0.4f %0.3f %0.1f %0.1f %0.4f %0.5f %0.7f" % (water.T, water.P.MPa, water.rho, water.u.kJkg, water.h.kJkg, water.s.kJkgK, water.x, water.virialB)) # 496.44 2.4691 300.000 1000.0 1008.2 2.6476 0.02680 -0.0096173 # # >>> water=H2O(s=3.e3, h=1000.e3) # >>> print("%0.2f %0.6f %0.5f %0.2f %0.1f %0.4f %0.5f %0.6f" % (water.T, water.P.MPa, water.rho, water.u.kJkg, water.h.kJkg, water.s.kJkgK, water.x, water.virialB)) # 345.73 0.034850 0.73526 952.60 1000.0 3.0000 0.29920 -0.034124 # # >>> water=H2O(u=995.e3, h=1000.e3) # >>> print("%0.2f %0.4f %0.2f %0.2f %0.1f %0.4f %0.5f %0.6f" % (water.T, water.P.MPa, water.rho, water.u.kJkg, water.h.kJkg, water.s.kJkgK, water.x, water.virialB)) # 501.89 2.7329 546.58 995.00 1000.0 2.6298 0.00866 -0.009308 # # >>> water=H2O(u=1000.e3, s=3.e3) # >>> print("%0.2f %0.6f %0.5f %0.2f %0.1f %0.4f %0.5f %0.6f" % (water.T, water.P.MPa, water.rho, water.u.kJkg, water.h.kJkg, water.s.kJkgK, water.x, water.virialB)) # 371.24 0.094712 1.99072 1000.00 1047.6 3.0000 0.28144 -0.025543 """ name = "water" CASNumber = "7732-18-5" formula = "H2O" synonym = "R-718" Tc = unidades.Temperature(647.096) rhoc = unidades.Density(322.) Pc = unidades.Pressure(22064.0, "kPa") M = 18.015268 # g/mol Tt = unidades.Temperature(273.16) Tb = unidades.Temperature(373.1243) f_acent = 0.3443 momentoDipolar = unidades.DipoleMoment(1.855, "Debye") id = 62 Fi1 = { "ao_log": [1, 3.00632], "pow": [0, 1], "ao_pow": [-8.3204464837497, 6.6832105275932], "ao_exp": [0.012436, 0.97315, 1.2795, 0.96956, 0.24873], "titao": [1.28728967, 3.53734222, 7.74073708, 9.24437796, 27.5075105] } Fi2 = { "ao_log": [1, 3.00392], "pow": [0, 1], "ao_pow": [8.203520690, -11.996306443], "ao_exp": [], "titao": [], "ao_hyp": [0.01059, -0.98763, 3.06904, 0], "hyp": [0.415386589, 1.763895929, 3.874803739, 0] } Fi3 = { "ao_log": [1, 3.00632], "pow": [0, 1], "ao_pow": [-8.318441, 6.681816], "ao_exp": [0.012436, 0.97315, 1.2795, 0.96956, 0.24873], "titao": [1.287202151, 3.537101709, 7.740210774, 9.243749421, 27.5056402] } Fi4 = { "ao_log": [1, 3.00632], "pow": [0, 1], "ao_pow": [-8.3177095, 6.6815049], "ao_exp": [0.012436, 0.97315, 1.2795, 0.96956, 0.24873], "titao": [1.287202151, 3.537101709, 7.740210774, 9.243749421, 27.5056402] } helmholtz1 = { "__type__": "Helmholtz", "__name__": u"Helmholtz equation of state for water of Wagner and Pruß (2002).", "__doi__": { "autor": u"Wagner, W., Pruß, A.", "title": "The IAPWS Formulation 1995 for the Thermodynamic Properties of Ordinary Water Substance for General and Scientific Use", "ref": "J. Phys. Chem. Ref. Data 31, 387 (2002)", "doi": "10.1063/1.1461829" }, "__test__": # Table 6.6, Pag 436 """ >>> wt=H2O() >>> tau=wt.Tc/500 >>> delta=838.025/wt.rhoc >>> print "%0.9g %0.9g %0.9g %0.9g %0.9g %0.9g" % wt._phi0(wt._constants["cp"], tau, delta) 2.04797734 9.04611106 -1.93249185 0.384236747 -0.147637878 0 >>> print "%0.9g %0.9g %0.9g %0.9g %0.9g %0.9g" % wt._phir(tau, delta)[:6] -3.42693206 -5.81403435 -2.23440737 -0.36436665 0.856063701 -1.12176915 >>> tau=wt.Tc/647 >>> delta=358/wt.rhoc >>> print "%0.9g %0.9g %0.9g %0.9g %0.9g %0.9g" % wt._phi0(wt._constants["cp"], tau, delta) -1.56319605 9.80343918 -3.43316334 0.899441341 -0.808994726 0 >>> print "%0.9g %0.9g %0.9g %0.9g %0.9g %0.9g" % wt._phir(tau, delta)[:6] -1.21202657 -3.21722501 -9.96029507 -0.714012024 0.475730696 -1.3321472 """ #Table 13.1, Pag 486 """ >>> st=H2O(T=273.16, x=0.5) >>> print "%0.6g %0.3g %0.6g %0.3g %0.3f %0.6g %0.4f %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\ st.T, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, st.Gas.h.kJkg, \ st.Liquido.s.kJkgK, st.Gas.s.kJkgK, st.Liquido.cv.kJkgK, st.Gas.cv.kJkgK, \ st.Liquido.cp.kJkgK, st.Gas.cp.kJkgK, st.Liquido.w, st.Gas.w) 273.16 0.000612 999.793 0.00485 0.001 2500.92 -0.0000 9.1555 4.2174 1.4184 4.2199 1.8844 1402.3 409 >>> st=H2O(T=300, x=0.5) >>> print "%0.6g %0.4g %0.6g %0.4g %0.3f %0.6g %0.4f %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\ st.T, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, st.Gas.h.kJkg, \ st.Liquido.s.kJkgK, st.Gas.s.kJkgK, st.Liquido.cv.kJkgK, st.Gas.cv.kJkgK, \ st.Liquido.cp.kJkgK, st.Gas.cp.kJkgK, st.Liquido.w, st.Gas.w) 300 0.003537 996.513 0.02559 112.565 2549.85 0.3931 8.5174 4.1305 1.4422 4.1809 1.9141 1501.4 427.89 >>> st=H2O(T=400, x=0.5) >>> print "%0.6g %0.5g %0.6g %0.5g %0.3f %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\ st.T, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, st.Gas.h.kJkg, \ st.Liquido.s.kJkgK, st.Gas.s.kJkgK, st.Liquido.cv.kJkgK, st.Gas.cv.kJkgK, \ st.Liquido.cp.kJkgK, st.Gas.cp.kJkgK, st.Liquido.w, st.Gas.w) 400 0.24577 937.486 1.3694 532.953 2715.7 1.6013 7.0581 3.6324 1.6435 4.2555 2.2183 1509.5 484.67 >>> st=H2O(T=500, x=0.5) >>> print "%0.6g %0.5g %0.6g %0.5g %0.3f %0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\ st.T, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, st.Gas.h.kJkg, \ st.Liquido.s.kJkgK, st.Gas.s.kJkgK, st.Liquido.cv.kJkgK, st.Gas.cv.kJkgK, \ st.Liquido.cp.kJkgK, st.Gas.cp.kJkgK, st.Liquido.w, st.Gas.w) 500 2.6392 831.313 13.199 975.431 2802.48 2.581 6.2351 3.2255 2.2714 4.6635 3.4631 1239.6 504.55 >>> st=H2O(T=600, x=0.5) >>> print "%0.6g %0.5g %0.6g %0.5g %0.6g %0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\ st.T, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, st.Gas.h.kJkg, \ st.Liquido.s.kJkgK, st.Gas.s.kJkgK, st.Liquido.cv.kJkgK, st.Gas.cv.kJkgK, \ st.Liquido.cp.kJkgK, st.Gas.cp.kJkgK, st.Liquido.w, st.Gas.w) 600 12.345 649.411 72.842 1505.36 2677.81 3.519 5.4731 3.0475 3.3271 6.9532 9.1809 749.57 457.33 >>> st=H2O(T=646, x=0.5) >>> print "%0.6g %0.5g %0.5g %0.5g %0.2f %0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\ st.T, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, st.Gas.h.kJkg, \ st.Liquido.s.kJkgK, st.Gas.s.kJkgK, st.Liquido.cv.kJkgK, st.Gas.cv.kJkgK, \ st.Liquido.cp.kJkgK, st.Gas.cp.kJkgK, st.Liquido.w, st.Gas.w) 646 21.775 402.96 243.46 1963.49 2238.06 4.2214 4.6465 4.5943 5.1457 204.58 385.23 297.13 331.61 >>> st=H2O(T=647, x=0.5) >>> print "%0.6g %0.5g %0.5g %0.5g %0.2f %0.6g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g %0.5g" % (\ st.T, st.P.MPa, st.Liquido.rho, st.Gas.rho, st.Liquido.h.kJkg, st.Gas.h.kJkg, \ st.Liquido.s.kJkgK, st.Gas.s.kJkgK, st.Liquido.cv.kJkgK, st.Gas.cv.kJkgK, \ st.Liquido.cp.kJkgK, st.Gas.cp.kJkgK, st.Liquido.w, st.Gas.w) 647 22.038 357.34 286.51 2029.44 2148.56 4.3224 4.5065 6.2344 6.274 3905.2 5334.1 251.19 285.32 """ #Table 13.2, Pag 495 """ >>> st=H2O(T=290, P=50000) >>> print "%0.6g %0.5g %0.5g %0.5g %0.4g %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) 290 998.78 70.725 70.775 0.2513 4.1682 4.1868 1472.2 >>> st=H2O(T=600, P=15000000) >>> print "%0.6g %0.6g %0.6g %0.6g %0.5g %0.5g %0.5g %0.5g" % (\ st.T, st.rho, st.u.kJkg, st.h.kJkg, st.s.kJkgK, st.cv.kJkgK, st.cp.kJkgK, st.w) 600 659.407 1474.91 1497.65 3.4994 3.0282 6.583 787.74 >>> st=H2O(T=640, P=22500000) >>> print "%0.6g %0.6g %0.6g %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) 640 531.385 1744.16 1786.5 3.9443 3.1724 12.142 529.59 >>> st=H2O(T=580, P=40000000) >>> print "%0.6g %0.6g %0.6g %0.6g %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) 580 753.362 1306.43 1359.53 3.2061 2.9966 4.9863 1099.6 """, "R": 8.314371357587, "cp": Fi1, "ref": { "Tref": Tt, "Pref": 0.611655, "ho": 0.611872, "so": 0 }, "Tmin": Tt, "Tmax": 2000., "Pmax": 2000000.0, "rhomax": 73.96, "Pmin": 0.61248, "rhomin": 55.49696, "nr1": [ 0.12533547935523e-1, 0.78957634722828e1, -0.87803203303561e1, 0.31802509345418, -0.26145533859358, -0.78199751687981e-2, 0.88089493102134e-2 ], "d1": [1, 1, 1, 2, 2, 3, 4], "t1": [-0.5, 0.875, 1, 0.5, 0.75, 0.375, 1], "nr2": [ -0.66856572307965, 0.20433810950965, -0.66212605039687e-4, -0.19232721156002, -0.25709043003438, 0.16074868486251, -0.4009282892587e-1, 0.39343422603254e-6, -0.75941377088144e-5, 0.56250979351888e-3, -0.15608652257135e-4, 0.11537996422951e-8, 0.36582165144204e-6, -0.13251180074668e-11, -0.62639586912454e-9, -0.10793600908932, 0.17611491008752e-1, 0.22132295167546, -0.40247669763528, 0.58083399985759, 0.49969146990806e-2, -0.31358700712549e-1, -0.74315929710341, 0.47807329915480, 0.20527940895948e-1, -0.13636435110343, 0.14180634400617e-1, 0.83326504880713e-2, -0.29052336009585e-1, 0.38615085574206e-1, -0.20393486513704e-1, -0.16554050063734e-2, 0.19955571979541e-2, 0.15870308324157e-3, -0.16388568342530e-4, 0.43613615723811e-1, 0.34994005463765e-1, -0.76788197844621e-1, 0.22446277332006e-1, -0.62689710414685e-4, -0.55711118565645e-9, -0.19905718354408, 0.31777497330738, -0.11841182425981 ], "c2": [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 4, 6, 6, 6, 6 ], "d2": [ 1, 1, 1, 2, 2, 3, 4, 4, 5, 7, 9, 10, 11, 13, 15, 1, 2, 2, 2, 3, 4, 4, 4, 5, 6, 6, 7, 9, 9, 9, 9, 9, 10, 10, 12, 3, 4, 4, 5, 14, 3, 6, 6, 6 ], "t2": [ 4, 6, 12, 1, 5, 4, 2, 13, 9, 3, 4, 11, 4, 13, 1, 7, 1, 9, 10, 10, 3, 7, 10, 10, 6, 10, 10, 1, 2, 3, 4, 8, 6, 9, 8, 16, 22, 23, 23, 10, 50, 44, 46, 50 ], "gamma2": [1] * 44, "nr3": [-0.31306260323435e2, 0.31546140237781e2, -0.25213154341695e4], "d3": [3] * 3, "t3": [0, 1, 4], "alfa3": [20] * 3, "beta3": [150, 150, 250], "gamma3": [1.21, 1.21, 1.25], "epsilon3": [1.] * 3, "nr4": [-0.14874640856724, 0.31806110878444], "a4": [3.5, 3.5], "b4": [0.85, 0.95], "B": [0.2, 0.2], "C": [28, 32], "D": [700, 800], "A": [0.32, .32], "beta4": [0.3, 0.3] } GERG = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for water 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": 1350.0, "Pmax": 1000000.0, "rhomax": 73.96, "Pmin": 0.61166, "rhomin": 55.497, "nr1": [ 0.82728408749586, -0.18602220416584e1, -0.11199009613744e1, 0.15635753976056, 0.87375844859025, -0.36674403715731, 0.53987893432436e-1 ], "d1": [1, 1, 1, 2, 2, 3, 4], "t1": [0.5, 1.25, 1.875, 0.125, 1.5, 1, 0.75], "nr2": [ 0.10957690214499e1, 0.53213037828563e-1, 0.13050533930825e-1, -0.41079520434476, 0.14637443344120, -0.55726838623719e-1, -0.11201774143800e-1, -0.66062758068099e-2, 0.46918522004538e-2 ], "c2": [1, 1, 1, 2, 2, 2, 3, 5, 5], "d2": [1, 5, 5, 1, 2, 4, 4, 1, 1], "t2": [1.5, 0.625, 2.625, 5, 4, 4.5, 3, 4, 6], "gamma2": [1] * 9 } helmholtz2 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for water of Saul and Wagner-58 coeff (1989).", "__doi__": { "autor": "Saul, A. and Wagner, W.", "title": "A Fundamental Equation for Water Covering the Range from the Melting Line to 1273 K at Pressures up to 25000 MPa", "ref": "J. Phys. Chem. Ref. Data 18, 1537 (1989)", "doi": "10.1063/1.555836" }, "R": 8.31434, "cp": Fi3, "ref": { "Tref": Tt, "Pref": 611.655, "ho": 0.611872, "so": 0 }, "Tmin": Tt, "Tmax": 1273., "Pmax": 400000.0, "rhomax": 55.49, "Pmin": 0.61166, "rhomin": 55.497, "nr1": [ 0.8216377478, -0.2543894379, -0.08830868648, -0.8903097248e-6, -0.1241333357e-5, 0.2895590286e-8, 0.1403610309e-10, 0.8183943371e-12, -0.2397905287e-12 ], "d1": [1, 1, 2, 5, 8, 11, 11, 13, 13], "t1": [0, 2, 0, 9, 0, 0, 12, 7, 13], "nr2": [ -0.7519743341, -0.4151278588, -0.103051374e1, -0.1648036888e1, -0.4686350251, 0.3560258142, -0.6364658294, 0.2227482363, -0.8954849939e-1, 0.1557686788e-2, 0.1347719088e-2, -0.1301353385e-2, 0.9987368673e-6, 0.2263629476e-3, 0.289330495e-5, 0.1995437169, -0.2707767662e-1, 0.1849068216e-1, -0.4402394357e-2, -0.8546876737e-1, 0.1220538576, -0.2562237041, 0.2555034636, -0.6323203907e-1, 0.3351397575e-4, -0.6152834985e-1, -0.3533048208e-3, 0.3146309259e-1, -0.2261795983e-2, 0.18689702e-3, -0.1384614556e-2, 0.2713160073e-2, -0.4866118539e-2, 0.3751789129e-2, -0.5692669373e-3, -0.5876414555, 0.5687838346, -0.1642158198, 0.5878635885, -0.2844301931, -0.2049198337, -0.4039233716e-2, 0.5459049594e-1, -0.8914260146e-2, 0.4974411254e-2 ], "c2": [1] * 15 + [2] * 20 + [3] * 10, "d2": [ 1, 1, 1, 2, 2, 3, 4, 4, 4, 5, 6, 7, 8, 9, 11, 1, 1, 1, 1, 2, 2, 4, 5, 6, 6, 7, 7, 8, 10, 10, 11, 11, 11, 11, 11, 2, 2, 2, 3, 3, 4, 4, 5, 5, 5 ], "t2": [ 0, 1, 3, 1, 5, 5, 2, 3, 5, 6, 4, 1, 8, 0, 1, 0, 9, 10, 11, 0, 8, 5, 4, 2, 12, 3, 10, 3, 2, 8, 0, 1, 3, 4, 6, 13, 14, 15, 14, 16, 13, 26, 15, 23, 25 ], "gamma2": [1] * 45, "nr5": [-0.709318338e-2, 0.1718796342e-1, -0.1482653038e-1, 0.4517292884e-2], "d5": [1, 2, 3, 4], "t5": [50, 40, 32, 26] } helmholtz3 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for water of Saul and Wagner-38 coeff (1989).", "__doi__": { "autor": "Saul, A. and Wagner, W.", "title": "A Fundamental Equation for Water Covering the Range from the Melting Line to 1273 K at Pressures up to 25000 MPa", "ref": "J. Phys. Chem. Ref. Data 18, 1537 (1989)", "doi": "10.1063/1.555836" }, "R": 8.31434, "cp": Fi4, "ref": { "Tref": Tt, "Pref": 611.655, "ho": 0.611872, "so": 0 }, "Tmin": Tt, "Tmax": 1273., "Pmax": 400000.0, "rhomax": 55.49, "Pmin": 0.61166, "rhomin": 55.497, "nr1": [ 0.2330009013, -0.1402091128e1, 0.1172248041, -0.1850749499, 0.1770110422, 0.5525151794e-1, -0.341325738e-3, 0.8557274367e-3, 0.3716900685e-3, -0.1308871233e-3, 0.3216895199e-4, 0.2785881034e-6 ], "d1": [1, 1, 2, 2, 2, 2, 3, 5, 5, 6, 7, 8], "t1": [0, 2, 0, 1, 2, 3, 5, 0, 1, 3, 2, 5], "nr2": [ -0.352151113, 0.7881914536e-1, -0.151966661e-1, -0.1068458586, -0.2055046288, 0.9146198012, 0.3213343569e-3, -0.1133591391e1, -0.3107520749, 0.1217901527e1, -0.4481710831, 0.5494218772e-1, -0.8665222096e-4, 0.3844084088e-1, 0.9853044884e-2, -0.1767598472e-1, 0.1488549222e-2, -0.3070719069e-2, 0.388080328e-2, -0.2627505215e-2, 0.5258371388e-3, -0.1716396901, 0.7188823624e-1, 0.5881268357e-1, -0.145593888e-1, -0.12161394e-1 ], "c2": [ 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3 ], "d2": [ 1, 1, 1, 2, 3, 3, 3, 4, 5, 5, 6, 7, 7, 8, 8, 9, 11, 11, 11, 11, 11, 2, 2, 3, 3, 5 ], "t2": [ 5, 7, 9, 5, 4, 6, 13, 5, 2, 3, 2, 0, 11, 1, 4, 0, 0, 3, 5, 6, 7, 13, 14, 15, 24, 15 ], "gamma2": [1] * 26 } helmholtz4 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for water 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.314371357587, "cp": Fi1, "ref": { "name": "CUSTOM", "Tref": Tt, "Pref": 611.655, "ho": 0.611872, "so": 0 }, "Tmin": Tt, "Tmax": 620.0, "Pmax": 800000.0, "rhomax": 40., "Pmin": 0.1, "rhomin": 40., "nr1": [ 3.46821920e-1, 5.03423025e-1, -3.51059570e-1, 5.07004866e-2, 1.99939129e-4, -5.69888763e-1 ], "d1": [1, 1, 1, 3, 7, 2], "t1": [1.5, 0.25, 1.25, 0.25, 0.875, 1.375], "nr2": [ -1.96198912e-1, -2.02509554, -1.09353609, 7.25785202e-2, 2.16072642e-1, -1.01542630e-1, 7.46926106e-2, 2.18830463e-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 = helmholtz1, GERG, helmholtz2, helmholtz3, helmholtz4 _PR = 0.0043451 _surface = {"sigma": [-0.1306, 0.2151], "exp": [2.471, 1.233]} _melting = {"Tmin": 251.165, "Tmax": 370.0} _sublimation = {"Tmin": 50.0, "Tmax": Tt} # _sublimation={"eq": 2, "Tref": 1, "Pref": 0.133332237, "a1": [-0.212144006e2, 0.273203819e2, -0.61059813e1], "exp1": [-0.9933333333, 0.206667, 0.703333], "a2": [], "exp2": [], "a3": [], "exp3": []} _vapor_Pressure = { "eq": 6, "ao": [ -7.85951783, 1.84408259, -11.7866497, 22.6807411, -15.9618719, 1.80122502 ], "exp": [2, 3, 6, 7, 8, 15] } _liquid_Density = { "eq": 2, "ao": [ 1.99274064, 1.09965342, -0.510839303, -1.75493479, -45.5170352, -6.74694450e5 ], "exp": [1, 2, 5, 16, 43, 110] } _vapor_Density = { "eq": 4, "ao": [ -2.0315024, -2.6830294, -5.38626492, -17.2991605, -44.7586581, -63.9201063 ], "exp": [1, 2, 4, 9, 18.5, 35.5] } visco0 = { "eq": 0, "method": "_visco0", "__name__": "IAPWS (1997)", "__code__": (_Viscosity, ) } visco1 = { "eq": 4, "omega": 1, "__name__": u"Quiñones-Cisneros (2006)", "__doi__": { "autor": "S.E.Quiñones-Cisneros and U.K. Deiters", "title": "Generalization of the Friction Theory for Viscosity Modeling", "ref": "J. Phys. Chem. B, 2006, 110 (25), pp 12820–12834", "doi": "10.1021/jp0618577" }, "Tref": 647.096, "muref": 1.0, "ek": 809.1, "sigma": 0.2641, "n_chapman": 0, "n_ideal": [151.138, -444.318, 398.262, -81.7008], "t_ideal": [0, 0.25, 0.5, 0.75], "a": [-1.17407105202836e-5, -3.78854818708520e-7, 3.56742875797909e-8], "b": [1.62216397984014e-6, -8.36595322447571e-6, 9.10862531286788e-8], "c": [1.92706925578893e-5, -1.28679815491711e-5, 0.0], "A": [-3.30144899918610e-10, 0.0, 1.02931444103415e-11], "B": [5.03139997945133e-10, 1.82304182380560e-10, 0.0], "C": [8.01449084635477e-10, 5.65613687804585e-9, 1.10163426018591e-10], "D": [0.0, 0.0, 0.0] } _viscosity = visco0, visco1 def _visco0(self, rho, T, fase): """IAPWS, Release on the IAPWS Formulation 2008 for the Viscosity of Ordinary Water Substance (International Association for the Properties of Water and Steam, 2008)""" ref = H2O() ref._ref("OTO") estado = ref._Helmholtz(rho, 1.5 * 647.096) drho = 1 / estado["dpdrho"] * 1e3 return _Viscosity(rho, T, fase=fase, drho=drho) thermo0 = { "eq": 0, "method": "_thermo0", "__name__": "IAPWS (1997)", "__code__": (_ThCond, ) } _thermal = thermo0, def _thermo0(self, rho, T, fase): """IAPWS, Release on the IAPWS Formulation 2011 for the Thermal Conductivity of Ordinary Water Substance""" ref = H2O() ref._ref("OTO") estado = ref._Helmholtz(rho, 1.5 * 647.096) drho = 1 / estado["dpdrho"] * 1e3 return _ThCond(rho, T, fase, drho) def _Dielectric(self, rho, T): return unidades.Dimensionless(_Dielectric(rho, T)) @classmethod def _Melting_Pressure(cls, T): # if 251.165 <= T <= 273.16: # Tita = T/cls.Tt # a = [0.119539337e7, 0.808183159e5, 0.33382686e4] # expo = [3., 0.2575e2, 0.10375e3] # suma = 1 # for ai, expi in zip(a, expo): # suma += ai*(1-Tita**expi) # P1 = suma*611.657 # else: # P1 = None if 251.165 <= T <= 256.164: Tref = 251.165 Pref = 208566. Tita = T / Tref P2 = Pref * (1 - 0.299948 * (1 - Tita**60.)) elif 256.164 < T <= 273.31: Tref = 256.164 Pref = 350100. Tita = T / Tref P2 = Pref * (1 - 1.18721 * (1 - Tita**8.)) elif 273.31 < T <= 355: Tref = 273.31 Pref = 632400. Tita = T / Tref P2 = Pref * (1 - 1.07476 * (1 - Tita**4.6)) elif 355. < T: Tref = 355 Pref = 2216000. Tita = T / Tref P2 = Pref * exp(1.73683 * (1 - 1. / Tita) - 0.544606e-1 * (1 - Tita**5) + 0.806106e-7 * (1 - Tita**22)) return unidades.Pressure(P2, "kPa") @classmethod def _Sublimation_Pressure(cls, T): Pref = 611.657 Tita = T / cls.Tt a = [-0.212144006e2, 0.273203819e2, -0.61059813e1] expo = [0.333333333e-2, 1.20666667, 1.70333333] suma = 0 for ai, expi in zip(a, expo): suma += ai * Tita**expi return unidades.Pressure(exp(suma / Tita) * Pref)
class F2(MEoS): """Multiparameter equation of state for fluorine""" name = "fluorine" CASNumber = "7782-41-4" formula = "F2" synonym = "" _refPropName = "FLUORINE" _coolPropName = "Fluorine" rhoc = unidades.Density(592.864) Tc = unidades.Temperature(144.414) Pc = unidades.Pressure(5172.4, "kPa") M = 37.99681 # g/mol Tt = unidades.Temperature(53.4811) Tb = unidades.Temperature(85.0368) f_acent = 0.0449 momentoDipolar = unidades.DipoleMoment(0.0, "Debye") id = 208 CP1 = { "ao": 3.5011231, "an": [ -0.60936946e-4 / 144.414**4, 0.63196690e-3 / 144.414**3, -0.74069617e-4 / 144.414**-2 ], "pow": [4, 3, -2], "ao_exp": [1.0127670], "exp": [1286.12], "ao_hyp": [], "hyp": [] } CP2 = { "ao": 0.7593432 / 8.3143 * 37.997, "an": [ 0.2883653e-3 / 8.3143 * 37.997, -0.4192916e-5 / 8.3143 * 37.997, 0.2309778e-7 / 8.3143 * 37.997, -0.3291582e-10 / 8.3143 * 37.997 ], "pow": [1, 2, 3, 4], "ao_exp": [], "exp": [], "ao_hyp": [], "hyp": [] } reuck = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for fluorine of de Reuck " "(1990).", "__doi__": { "autor": "de Reuck, K.M.", "title": "International thermodynamic tables of the fluid " "state: Vol. 11 - fluorine", "ref": "Pergamon Press, Oxford, 1990.", "doi": "" }, "R": 8.31448, "cp": CP1, "ref": "NBP", "Tmin": Tt, "Tmax": 300.0, "Pmax": 20000.0, "rhomax": 45.47, "Pmin": 0.23881, "rhomin": 44.917, "nr1": [ 0.151144749736e1, -0.298666288409e1, 0.329644905098e1, -0.298458624201e1, -0.228688966459e1, -0.109492193400e1, 0.304775277572e1, 0.115689564208, -0.116100171627e1, 0.295656394476, 0.711482542928e-1, -0.171363832155e-2, 0.665317955515e-3 ], "d1": [1, 1, 1, 1, 2, 2, 3, 3, 4, 4, 5, 8, 9], "t1": [0, 0.5, 1.5, 2., 0.5, 1., 0.5, 2., 0.5, 1., 0., 0.5, 0], "nr2": [ 0.506026676251e1, -0.629268435440e1, 0.617784808739e1, -0.155366191788e1, -0.287170687343e1, 0.317214480494e1, -0.267969025215e1, 0.271865479252e1, -0.107191065039e1, 0.126597342291e1, -0.706244695489, 0.268707888826, 0.527251190274e-1, 0.544411481926e-1, 0.228949994105e-3, -0.547908264304e-9, -0.964273224950e-1, 0.368084486225e-3 ], "d2": [2, 2, 2, 2, 3, 3, 3, 4, 4, 4, 5, 6, 7, 8, 12, 4, 6, 6], "t2": [1, 3, 4, 5, 1, 4, 5, 1, 3, 5, 4, 4, 1, 1, 5, 30, 20, 25], "c2": [2] * 18, "gamma2": [1.07810258] * 15 + [2.15620515, 3.23430773, 3.23430773] } polt = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for fluorine of Polt (1992).", "__doi__": { "autor": "Polt, A., Platzer, B., and Maurer, G.", "title": "Parameter der thermischen Zustandsgleichung von " "Bender fuer 14 mehratomige reine Stoffe", "ref": "Chem. Technik 22(1992)6 , 216/224", "doi": "" }, "R": 8.3143, "cp": CP2, "ref": "NBP", "Tmin": Tt, "Tmax": 300.0, "Pmax": 25000.0, "rhomax": 45.14, "Pmin": 0.25394, "rhomin": 44.89, "nr1": [ 0.862212325175e-2, 0.162286882091, -0.228707299586e-1, 0.624951179331, -0.158918489879e1, 0.195171867807, -0.438453517535, 0.402200928405e-1, 0.319444405579e-1, 0.161784325978e-1, 0.230132378392, 0.819206229044e-1, -0.173741828076, 0.137942204542e-1, -0.449971813506e-2, 0.756554661780e-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.862212325175e-2, -0.162286882091, 0.228707299586e-1, 0.184612089745, -0.425779777811, 0.825656492996e-1 ], "d2": [0, 0, 0, 2, 2, 2], "t2": [3, 4, 5, 3, 4, 5], "c2": [2] * 6, "gamma2": [0.9225328] * 6 } eq = reuck, polt _surface = {"sigma": [0.03978], "exp": [1.218]} _melting = { "eq": 1, "Tref": Tt, "Pref": 1000, "Tmin": Tt, "Tmax": 300.0, "a1": [.000252, 249.975, -249.9750131], "exp1": [0, 2.1845, 0], "a2": [], "exp2": [], "a3": [], "exp3": [] } _vapor_Pressure = { "eq": 3, "n": [-0.37061e1, -0.81517e2, 0.13743e3, -0.58617e2, -0.13528e1], "t": [1.0, 1.50, 1.61, 1.77, 7.3] } _liquid_Density = { "eq": 1, "n": [0.21286, 0.44011e1, -0.53959e1, 0.41347e1, -0.97544], "t": [0.228, 0.58, 0.908, 1.24, 1.6] } _vapor_Density = { "eq": 2, "n": [-3.6218, -55.363, 122.14, -230.92, -338.61, 432.18], "t": [0.454, 2.3, 2.9, 4.0, 6.0, 5.3] }
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") Fi0 = {"ao_log": [1, 3], "pow": [0, 1, 2, 3, 4, 5], "ao_pow": [-8.6739710041, 6.9611755531], "ao_exp": [0.00863, 0.97454, 2.0646, 0.23528, 0.29555], "titao": [0.4255669437, 2.6093155672, 6.0185106089, 11.3380974051, 29.5101165339], "ao_hyp": [], "hyp": []} 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__": "Helmholtz equation of state for heavy water of Herrig" " (2017).", "__doi__": { "autor": "Herrig", "title": "Preliminary helmholtz equation of state for Heavy Water", "ref": "", "doi": ""}, "R": 8.3144621, "rhoref": 17.77555*M, "Tref": 643.847, "cp": Fi0, "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.0105835, 0.99127253, -1.224122, 1.710643, -2.189443, 0.1145315], "d1": [4.0, 1.0, 1.0, 2.0, 2.0, 3.0], "t1": [1.0, 0.463, 1.29, 1.307, 1.2165, 0.587], "nr2": [-0.89875532, -1.597051, -2.804509, 0.33016885, -3.396526, -0.001881], "c2": [1.0, 2.0, 2.0, 1.0, 2.0, 1.0], "d2": [1.0, 1.0, 3.0, 2.0, 2.0, 8.0], "t2": [2.95, 1.713, 1.929, 0.94, 3.033, 0.765], "gamma2": [1]*6, "nr3": [-0.70355957, -0.20345481, -0.70691398, 2.094255, 3.042546, 0.8010728, 0.213384, 0.32335789, -0.0245055, 0.7380677, -0.21484089], "t3": [1.504, 2.85, 1.96, 0.969, 2.576, 2.79, 3.581, 3.67, 1.7, 1.0, 4.1], "d3": [1.0, 2.0, 3.0, 1.0, 3.0, 1.0, 1.0, 2.0, 2.0, 2.0, 1.0], "beta3": [0.907, 0.48, 1.223, 2.61, 4.283, 1.4, 0.735, 0.24, 1067.0, 13.27, 1.48], "alfa3": [0.982, 1.34, 1.658, 1.6235, 1.4, 2.206, 0.84, 1.535, 11.33, 3.86, 7.56], "epsilon3": [2.272, 1.375, 0.648, 0.8925, 0.145, 0.291, 2.01, 1.08, 0.96, 0.181, 0.529], "gamma3": [2.263, 2.343, 0.929, 1.0, 1.383, 0.968, 1.695, 2.23, 1.07, 1.297, 2.41]} helmholtz2 = { "__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, helmholtz2 _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): mu = _D2O_Viscosity(rho, T) return unidades.Viscosity(mu) _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.mWmK/0.742128)) 762.915707396 >>> st=D2O(T=0.431*D2O.Tc, rho=3.23*D2O.rhoc) >>> print("%0.9f" % (st.k.mWmK/0.742128)) 833.912049618 >>> st=D2O(T=0.6*D2O.Tc, rho=2.95*D2O.rhoc) >>> print("%0.9f" % (st.k.mWmK/0.742128)) 861.240794445 >>> st=D2O(T=D2O.Tc, rho=0.7*D2O.rhoc) >>> print("%0.9f" % (st.k.mWmK/0.742128)) 469.015122112 >>> st=D2O(T=0.9*D2O.Tc, rho=0.08*D2O.rhoc) >>> print("%0.9f" % (st.k.mWmK/0.742128)) 74.522283066 >>> st=D2O(T=1.1*D2O.Tc, rho=0.98*D2O.rhoc) >>> print("%0.9f" % (st.k.mWmK/0.742128)) 326.652382218 >>> st=D2O(T=1.2*D2O.Tc, rho=0.8*D2O.rhoc) >>> print("%0.9f" % (st.k.mWmK/0.742128)) 259.605241187 """, } def _thermo0(self, rho, T, fase): k = _D2O_ThCond(rho, T) return unidades.ThermalConductivity(k) _thermal = thermo0,
class MD4M(MEoS): """Multiparameter equation of state for tetradecamethylhexasiloxane""" name = "tetradecamethylhexasiloxane" CASNumber = "107-52-8" formula = "C14H42O5Si6" synonym = "MD4M" rhoc = unidades.Density(285.6576532213632) Tc = unidades.Temperature(653.2) Pc = unidades.Pressure(877.47, "kPa") M = 458.99328 # g/mol Tt = unidades.Temperature(214.15) Tb = unidades.Temperature(532.723) f_acent = 0.825 momentoDipolar = unidades.DipoleMoment(1.308, "Debye") CP1 = { "ao": -20.071, "an": [2228.5e-3, -1311.4e-6, 286.2e-9], "pow": [1, 2, 3], "ao_exp": [], "exp": [], "ao_hyp": [], "hyp": [] } helmholtz1 = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for MD4M 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=MD4M(T=653.2, P=877470) >>> print "%0.6f" % st.v 0.003501 """, # Table 18, Pag 204 "R": 8.314472, "cp": CP1, "ref": "NBP", "Tmin": 300, "Tmax": 673.0, "Pmax": 30000.0, "rhomax": 2.09, "Pmin": 0.000000001, "rhomin": 2.09, "nr1": [ 1.18492421, -1.87465636, -0.65713510e-1, -0.61812689, 0.19535804, 0.50678740e-3 ], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.25, 1.125, 1.5, 1.375, 0.25, 0.875], "nr2": [ 1.23544082, 0.49462708e-1, -0.73685283, -0.19991438, -0.55118673e-1, 0.28325885e-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.10532e2, 0.33939e1, -0.89744e1, -0.56150e1], "exp": [1.0, 1.5, 2.75, 5.1] } _liquid_Density = { "eq": 1, "ao": [0.10453e1, 0.55476, 0.44536e1, -0.76179e1, 0.46237e1], "exp": [0.235, 0.6, 0.95, 1.35, 1.7] } _vapor_Density = { "eq": 3, "ao": [ -0.10890e1, -0.84374e1, -0.35615e2, -0.73478e3, 0.19915e4, -0.16317e4 ], "exp": [0.231, 0.8, 2.9, 7.7, 9.0, 10.0] }
class nC12(MEoS): """Multiparameter equation of state for n-dodecane""" name = "dodecane" CASNumber = "112-40-3" formula = "CH3-(CH2)10-CH3" synonym = "" _refPropName = "C12" _coolPropName = "n-Dodecane" rhoc = unidades.Density(226.5453372) Tc = unidades.Temperature(658.1) Pc = unidades.Pressure(1817.0, "kPa") M = 170.33484 # g/mol Tt = unidades.Temperature(263.6) Tb = unidades.Temperature(489.3) f_acent = 0.574 momentoDipolar = unidades.DipoleMoment(0.0, "Debye") id = 16 CP1 = {"ao": 23.085, "ao_exp": [37.776, 29.369, 12.461, 7.7733], "exp": [1280, 2399, 5700, 13869]} lemmon = { "__type__": "Helmholtz", "__name__": "short Helmholtz equation of state for dodecane of Lemmon " "(2004).", "__doi__": {"autor": "Lemmon, E.W., Huber, M.L.", "title": "Thermodynamic Properties of n-Dodecane", "ref": "Energy & Fuels, 18(4) (2004) 960-967", "doi": "10.1021_ef0341062"}, "R": 8.314472, "cp": CP1, "ref": "NBP", "Tmin": Tt, "Tmax": 700., "Pmax": 700000.0, "rhomax": 4.53, "nr1": [1.38031, -2.85352, .288897, -0.165993, .0923993, .000282772], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.32, 1.23, 1.5, 1.4, 0.07, 0.8], "nr2": [.956627, .0353076, -0.445008, -0.118911, -0.0366475, .0184223], "d2": [2, 5, 1, 4, 3, 4], "t2": [2.16, 1.1, 4.1, 5.6, 14.5, 12.], "c2": [1, 1, 2, 2, 3, 3], "gamma2": [1]*6} eq = lemmon, _PR = [0.1099, -26.8035] _surface = {"sigma": [0.0154, 0.0480], "exp": [4.18, 1.17]} _vapor_Pressure = { "eq": 3, "n": [-0.94217e1, -0.41890e1, 0.54999e1, -0.67789e1, -0.17161e1], "t": [1.0, 1.5, 1.359, 3.56, 9.2]} _liquid_Density = { "eq": 1, "n": [0.92236, 0.92047, 0.55713e1, -0.92253e1, 0.51763e1], "t": [0.21, 0.49, 1.08, 1.49, 1.9]} _vapor_Density = { "eq": 2, "n": [-1.7859, -7.5436, -22.848, -81.355, 92.283, -217.25], "t": [0.298, 0.91, 2.8, 6., 9., 11.]} visco0 = {"__name__": "Huber (2004)", "__doi__": { "autor": "Huber, M.L., Laesecke, A., Perkins, R.A.", "title": "Transport Properties of n-Dodecane", "ref": "Energy & Fuels 18(4) (2004) 968-975.", "doi": "10.1021/ef034109e"}, "eq": 1, "omega": 1, "ek": 522.592, "sigma": 0.735639, "n_chapman": 0.021357, "collision": [0.382987, -0.561050, 0.0313962], "Tref_virial": 522.592, "n_virial": [-19.572881, 219.73999, -1015.3226, 2471.0125, -3375.1717, 2491.6597, -787.26086, 14.085455, -0.34664158], "t_virial": [0, -0.25, -0.5, -0.75, -1, -1.25, -1.5, -2.5, -5.5], "Tref_res": 658.1, "rhoref_res": 1.33*M, "muref_res": 1000, "nr": [-0.0471703, 0.00827816, 0.0298429, -0.0134156], "tr": [1, 1, 2, 2], "dr": [2, 3, 2, 3], "CPf": 503.109, "CPg1": 2.32661, "CPgi": [2.23089/2.32661], "CPti": [-0.5]} _viscosity = visco0, thermo0 = {"eq": 1, "__name__": "Huber (2004)", "__doi__": { "autor": "Huber, M.L., Laesecke, A., Perkins, R.A.", "title": "Transport Properties of n-Dodecane", "ref": "Energy & Fuels 18(4) (2004) 968-975.", "doi": "10.1021/ef034109e"}, "Toref": 658.1, "koref": 1., "no": [0.436343e-2, -0.264054e-1, 0.922394e-1, -0.291756e-1], "to": [0, 1, 2, 3], "Tref_res": 658.1, "rhoref_res": 1.33*M, "kref_res": 1., "nr": [0.693347e-1, -0.280792e-1, -0.331695e-1, 0.173922e-2, 0.676165e-2, 0.309558e-2], "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": 1.52e-9, "Tcref": 987.15} _thermal = thermo0,
class NF3(MEoS): """Multiparameter equation of state for nitrogen trifluoride""" name = "nitrogen trifluoride" CASNumber = "7783-54-2" formula = "NF3" synonym = "" rhoc = unidades.Density(562.47) Tc = unidades.Temperature(234.0) Pc = unidades.Pressure(4460.7, "kPa") M = 71.019 # g/mol Tt = unidades.Temperature(66.36) Tb = unidades.Temperature(144.138) f_acent = 0.126 momentoDipolar = unidades.DipoleMoment(0.235, "Debye") # id = 951 id = 60 CP1 = { "ao": -7.140693612211, "an": [ 0.7427518245951e6, -0.4389825372134e5, 0.1012629224351e4, 0.5481339146452e-1, -0.7677196006769e-4, 0.4203630864340e-7 ], "pow": [-3, -2, -1.001, 1, 2, 3], "ao_exp": [-0.6328752997967], "exp": [3000], "ao_hyp": [], "hyp": [] } MBWR = { "__type__": "MBWR", "__name__": "MBWR equation of state for nitrogen trifluoride 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.314471, "cp": CP1, "Tmin": 85.0, "Tmax": 500.0, "Pmax": 50000.0, "rhomax": 25.3, "Pmin": 0.000186, "rhomin": 26.32, "b": [ None, 0.1774353868e-1, -0.5409379418, 0.3976634466e1, -0.5209476694e3, -0.3286322888e5, -0.5990517411e-3, 0.9217525601, -0.4848977075e3, -0.4235892691e7, -0.9824248063e-5, 0.5432235989e-1, -0.1462388500e2, -0.3366180440e-2, 0.2801374599, 0.8435288597e1, -0.1324421452e-1, 0.1875604377e-3, 0.2959643991, -0.700997687e-2, 0.4365820912e7, -0.1111397536e8, 0.2411866612e5, 0.3179136276e7, 0.6166849090e2, 0.4260854720e2, 0.1090598789, -0.3340951059e2, 0.8597429644e-4, 0.1240544214e-2, 0.1286224248e-6, -0.8941104276e-6, 0.3353054595e-4 ] } eq = MBWR, _vapor_Pressure = { "eq": 5, "ao": [-0.66672e1, 0.33864e1, -0.28222e1, -0.50602e1, 0.32481e1], "exp": [1.0, 1.5, 1.7, 5.5, 7.0] } _liquid_Density = { "eq": 1, "ao": [0.22080e1, 0.35709e2, -0.92868e2, 0.66666e2, -0.93589e1], "exp": [0.35, 2.4, 2.7, 3.0, 4.0] } _vapor_Density = { "eq": 3, "ao": [-0.3061e1, -0.80541e1, -0.19619e2, -0.13432e2, -0.3276e2, -0.67907e2], "exp": [0.421, 1.48, 3.9, 7.0, 8.0, 15.0] }
class RE245fa2(MEoS): """Multiparameter equation of state for RE245fa2""" name = "2,2,2-trifluoroethyl-difluoromethyl-ether" CASNumber = "1885-48-9" formula = "CHF2OCH2CF3" synonym = "HFE-245fa2" _refPropName = "RE245FA2" _coolPropName = "" rhoc = unidades.Density(515.001169364688) Tc = unidades.Temperature(444.88) Pc = unidades.Pressure(3433., "kPa") M = 150.047336 # g/mol Tt = unidades.Temperature(250) Tb = unidades.Temperature(302.4) f_acent = 0.387 momentoDipolar = unidades.DipoleMoment(1.631, "Debye") # id = 1874 CP1 = { "ao": 5.259865, "an": [], "pow": [], "ao_exp": [], "exp": [], "ao_hyp": [12.12843, 13.25677, 0.521867, 0], "hyp": [486, 1762, 7631, 0] } zhou = { "__type__": "Helmholtz", "__name__": "Helmholtz equation of state for RE245fa2 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": 400000.0, "rhomax": 10.02, "Pmin": 8.272, "rhomin": 10., "nr1": [0.47771378e-1, 0.15745383e1, -0.24763491e1, -0.49414564, 0.19380498], "d1": [4, 1, 1, 2, 3], "t1": [1, 0.32, 0.91, 1.265, 0.4266], "nr2": [-0.97863158, -0.42660297, 0.85352583, -0.53380114, -0.29780036e-1], "d2": [1, 3, 2, 2, 7], "t2": [2.24, 1.64, 1.65, 3.28, 0.855], "c2": [2, 2, 1, 2, 1], "gamma2": [1] * 5, "nr3": [0.97659111, -0.33121365, -0.14122591, -0.15312295e2], "d3": [1, 1, 3, 3], "t3": [1.227, 3.0, 4.3, 2.5], "alfa3": [1.005, 1.515, 1.156, 17.7], "beta3": [2, 3.42, 1.37, 471], "gamma3": [1.084, 0.72, 0.49, 1.152], "epsilon3": [0.723, 0.9488, 0.818, 0.891] } eq = zhou, _vapor_Pressure = { "eq": 3, "n": [-8.9235, 10.527, -23.058, 30.291, -20.913, -26.745], "t": [1, 1.5, 1.9, 2.4, 2.9, 3.2] } _liquid_Density = { "eq": 1, "n": [1.2479, 5.5732, -12.26, 13.964, -6.0384], "t": [0.34, 0.75, 1.2, 1.7, 2.3] } _vapor_Density = { "eq": 2, "n": [-0.667, -5.8238, -26.927, 21.574, -65.645], "t": [0.28, 0.66, 2.6, 3.5, 5.2] }
class Cis_2_butene(MEoS): """Multiparameter equation of state for cis-butene""" name = "cis-butene" CASNumber = "590-18-1" formula = "CH3-CH=CH-CH3" synonym = "" _refPropName = "C2BUTENE" _coolPropName = "cis-2-Butene" rhoc = unidades.Density(238.11522208) Tc = unidades.Temperature(435.75) Pc = unidades.Pressure(4225.5, "kPa") M = 56.10632 # g/mol Tt = unidades.Temperature(134.3) Tb = unidades.Temperature(276.87) f_acent = 0.202 momentoDipolar = unidades.DipoleMoment(0.3, "Debye") id = 25 Fi1 = { "ao_log": [1, 2.9687], "pow": [0, 1], "ao_pow": [0.2591542, 2.4189888], "ao_exp": [3.2375, 7.0437, 11.414, 7.3722], "titao": [248 / Tc, 1183 / Tc, 2092 / Tc, 4397 / Tc] } lemmon = { "__type__": "Helmholtz", "__name__": "short Helmholtz equation of state for 1-butene of Lemmon " "and Ihmels (2005)", "__doi__": { "autor": "Lemmon, E.W., Ihmels, E.C.", "title": "Thermodynamic properties of the butenes: Part " "II. Short fundamental equations of state", "ref": "Fluid Phase Equilibria 228-229 (2005) 173-187", "doi": "10.1016/j.fluid.2004.09.004" }, "R": 8.314472, "cp": Fi1, "ref": "NBP", "Tmin": Tt, "Tmax": 525., "Pmax": 50000.0, "rhomax": 14.09, "Pmin": 0.00026, "rhomin": 14.09, "nr1": [0.77827, -2.8064, 1.003, 0.013762, 0.085514, 0.00021268], "d1": [1, 1, 1, 2, 3, 7], "t1": [0.12, 1.3, 1.74, 2.1, 0.28, 0.69], "nr2": [0.22962, -0.072442, -0.23722, -0.074071, -0.026547, 0.012032], "d2": [2, 5, 1, 4, 3, 4], "t2": [0.75, 2., 4.4, 4.7, 15., 14.], "c2": [1, 1, 2, 2, 3, 3], "gamma2": [1] * 6 } eq = lemmon, _vapor_Pressure = { "eq": 3, "n": [-0.70022e1, 0.13695e1, -0.30509e1, 0.10012, -0.15577e1], "t": [1.0, 1.5, 3.2, 3.46, 6.4] } _liquid_Density = { "eq": 1, "n": [0.46849e1, -0.54614e1, 0.34718e1, 0.50511e1, -0.50389e1], "t": [0.402, 0.54, 0.69, 6.6, 7.0] } _vapor_Density = { "eq": 2, "n": [ -0.28918e1, -0.58582e1, -0.17443e2, -0.24566e2, -0.29413e2, -0.11392e3 ], "t": [0.4098, 1.174, 3.11, 6.1, 7.6, 14.8] }