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 test_ThCond(self):
# Table 7, pag 12
self.assertEqual(round(_ThCond(0, 298.15)*1000, 7), 18.4341883)
self.assertEqual(round(_ThCond(998, 298.15)*1000, 6), 607.712868)
self.assertEqual(round(_ThCond(1200, 298.15)*1000, 6), 799.038144)
self.assertEqual(round(_ThCond(0, 873.15)*1000, 7), 79.1034659)
# Table 8, pag 13
fluid = H2O(rho=1, T=647.35)
self.assertEqual(round(fluid.k.mWmK, 7), 51.9298924)
fluid = H2O(rho=122, T=647.35)
self.assertEqual(round(fluid.k.mWmK, 6), 130.922885)
fluid = H2O(rho=222, T=647.35)
self.assertEqual(round(fluid.k.mWmK, 6), 367.787459)
fluid = H2O(rho=272, T=647.35)
self.assertEqual(round(fluid.k.mWmK, 6), 757.959776)
fluid = H2O(rho=322, T=647.35)
self.assertEqual(round(fluid.k.mWmK, 5), 1443.75556)
fluid = H2O(rho=372, T=647.35)
self.assertEqual(round(fluid.k.mWmK, 6), 650.319402)
fluid = H2O(rho=422, T=647.35)
self.assertEqual(round(fluid.k.mWmK, 6), 448.883487)
fluid = H2O(rho=750, T=647.35)
self.assertEqual(round(fluid.k.mWmK, 6), 600.961346)
def test_ThCond(self):
# Table 7, pag 12
self.assertEqual(round(_ThCond(0, 298.15)*1000, 7), 18.4341883)
self.assertEqual(round(_ThCond(998, 298.15)*1000, 6), 607.712868)
self.assertEqual(round(_ThCond(1200, 298.15)*1000, 6), 799.038144)
self.assertEqual(round(_ThCond(0, 873.15)*1000, 7), 79.1034659)
# Table 8, pag 13
fluid = H2O(rho=1, T=647.35)
self.assertEqual(round(fluid.k.mWmK, 7), 51.9298924)
fluid = H2O(rho=122, T=647.35)
self.assertEqual(round(fluid.k.mWmK, 6), 130.922885)
fluid = H2O(rho=222, T=647.35)
self.assertEqual(round(fluid.k.mWmK, 6), 367.787459)
fluid = H2O(rho=272, T=647.35)
self.assertEqual(round(fluid.k.mWmK, 6), 757.959776)
fluid = H2O(rho=322, T=647.35)
self.assertEqual(round(fluid.k.mWmK, 5), 1443.75556)
fluid = H2O(rho=372, T=647.35)
self.assertEqual(round(fluid.k.mWmK, 6), 650.319402)
fluid = H2O(rho=422, T=647.35)
self.assertEqual(round(fluid.k.mWmK, 6), 448.883487)
fluid = H2O(rho=750, T=647.35)
self.assertEqual(round(fluid.k.mWmK, 6), 600.961346)
def _thermo0(self, rho, T, fase):
ref = H2O()
ref._ref(False)
estado = ref._eq(rho, 1.5*self.Tc)
delta = estado["delta"]
fird = estado["fird"]
firdd = estado["firdd"]
dpdrho = self.R*estado["T"]*(1+2*delta*fird+delta**2*firdd)
drho = 1/dpdrho*1e6
# convert values to IAPWS units
# ∂ρ/∂P]τ, [kg/m³·MPa]
# cp, [kJ/kg]
if fase:
fase = copy(fase)
fase.drhodP_T *= 1e6
fase.cp /= 1e3
return _ThCond(rho, T, fase, drho)
def _thermo0(self, rho, T, fase):
ref = H2O()
ref._ref(False)
estado = ref._eq(rho, 1.5*self.Tc)
delta = estado["delta"]
fird = estado["fird"]
firdd = estado["firdd"]
dpdrho = self.R*estado["T"]*(1+2*delta*fird+delta**2*firdd)
drho = 1/dpdrho*1e6
# convert values to IAPWS units
# ∂ρ/∂P]τ, [kg/m³·MPa]
# cp, [kJ/kg]
if fase:
fase = copy(fase)
fase.drhodP_T *= 1e6
fase.cp /= 1e3
return _ThCond(rho, T, fase, drho)
