class MKS_FunctionTester(unittest.TestCase):
def setUp(self):
self.maxError = 1E-6
self.steamTable = XSteam(XSteam.UNIT_SYSTEM_MKS)
def tearDown(self):
pass
def test_tsat_p(self):
error = self.steamTable.tsat_p(1.0) - 99.60591861
self.assertLess(error, self.maxError, 'tsat_p not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_t_ph(self):
error = self.steamTable.t_ph(1.0, 100.0) - 23.84481908
self.assertLess(error, self.maxError, 't_ph not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_t_ps(self):
error = self.steamTable.t_ps(1.0, 1.0) - 73.70859421
self.assertLess(error, self.maxError, 't_ps not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_t_hs(self):
error = self.steamTable.t_hs(100.0, 0.2) - 13.84933511
self.assertLess(error, self.maxError, 't_hs not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_psat_t(self):
error = self.steamTable.psat_t(100.0) - 1.014179779
self.assertLess(error, self.maxError, 'psat_t not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_p_hs(self):
error = self.steamTable.p_hs(84.0, 0.296) - 2.295498269
self.assertLess(error, self.maxError, 'p_hs not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_hV_p(self):
error = self.steamTable.hV_p(1.0) - 2674.949641
self.assertLess(error, self.maxError, 'hV_p not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_hL_p(self):
error = self.steamTable.hL_p(1.0) - 417.4364858
self.assertLess(error, self.maxError, 'hL_p not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_hV_t(self):
error = self.steamTable.hV_t(100.0) - 2675.572029
self.assertLess(error, self.maxError, 'hV_t not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_hL_t(self):
error = self.steamTable.hL_t(100.0) - 419.099155
self.assertLess(error, self.maxError, 'hL_t not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_h_pt(self):
error = self.steamTable.h_pt(1.0, 20.0) - 84.01181117
self.assertLess(error, self.maxError, 'h_pt not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_h_ps(self):
error = self.steamTable.h_ps(1.0, 1.0) - 308.6107171
self.assertLess(error, self.maxError, 'h_ps not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_h_px(self):
error = self.steamTable.h_px(1.0, 0.5) - 1546.193063
self.assertLess(error, self.maxError, 'h_px not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_h_prho(self):
error = self.steamTable.h_prho(1.0, 2.0) - 1082.773391
self.assertLess(error, self.maxError, 'h_prho not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_h_tx(self):
error = self.steamTable.h_tx(100.0, 0.5) - 1547.33559211
self.assertLess(error, self.maxError, 'h_tx not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_vV_p(self):
error = self.steamTable.vV_p(1.0) - 1.694022523
self.assertLess(error, self.maxError, 'vV_p not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_vL_p(self):
error = self.steamTable.vL_p(1.0) - 0.001043148
self.assertLess(error, self.maxError, 'vL_p not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_vV_t(self):
error = self.steamTable.vV_t(100.0) - 1.671860601
self.assertLess(error, self.maxError, 'vV_t not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_vL_t(self):
error = self.steamTable.vL_t(100.0) - 0.001043455
self.assertLess(error, self.maxError, 'vL_t not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_v_pt(self):
error = self.steamTable.v_pt(1.0, 100.0) - 1.695959407
self.assertLess(error, self.maxError, 'v_pt not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_v_ph(self):
error = self.steamTable.v_ph(1.0, 1000.0) - 0.437925658
self.assertLess(error, self.maxError, 'v_ph not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_v_ps(self):
error = self.steamTable.v_ps(1.0, 5.0) - 1.03463539
self.assertLess(error, self.maxError, 'v_ps not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_rhoV_p(self):
error = self.steamTable.rhoV_p(1.0) - 0.590310924
self.assertLess(error, self.maxError, 'rhoV_p not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_rhoL_p(self):
error = self.steamTable.rhoL_p(1.0) - 958.6368897
self.assertLess(error, self.maxError, 'rhoL_p not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_rhoV_t(self):
error = self.steamTable.rhoV_t(100.0) - 0.598135993
self.assertLess(error, self.maxError, 'rhoV_t not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_rhoL_t(self):
error = self.steamTable.rhoL_t(100.0) - 958.3542773
self.assertLess(error, self.maxError, 'rhoL_t not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_rho_pt(self):
error = self.steamTable.rho_pt(1.0, 100.0) - 0.589636754
self.assertLess(error, self.maxError, 'rho_pt not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_rho_ph(self):
error = self.steamTable.rho_ph(1.0, 1000.0) - 2.283492601
self.assertLess(error, self.maxError, 'rho_ph not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_rho_ps(self):
error = self.steamTable.rho_ps(1.0, 1.0) - 975.6236788
self.assertLess(error, self.maxError, 'rho_ps not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_sV_p(self):
error = self.steamTable.sV_p(0.006117) - 9.155465556
self.assertLess(error, self.maxError, 'sV_p not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_sL_p(self):
error = self.steamTable.sL_p(0.0061171) - 1.8359e-05
self.assertLess(error, self.maxError, 'sL_p not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_sV_t(self):
error = self.steamTable.sV_t(0.0001) - 9.155756716
self.assertLess(error, self.maxError, 'sV_t not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_sL_t(self):
error = self.steamTable.sL_t(100.0) - 1.307014328
self.assertLess(error, self.maxError, 'sL_t not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_s_pt(self):
error = self.steamTable.s_pt(1.0, 20.0) - 0.296482921
self.assertLess(error, self.maxError, 's_pt not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_s_ph(self):
error = self.steamTable.s_ph(1.0, 84.01181117) - 0.296813845
self.assertLess(error, self.maxError, 's_ph not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_uV_p(self):
error = self.steamTable.uV_p(1.0) - 2505.547389
self.assertLess(error, self.maxError, 'uV_p not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_uL_p(self):
error = self.steamTable.uL_p(1.0) - 417.332171
self.assertLess(error, self.maxError, 'uL_p not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_uV_t(self):
error = self.steamTable.uV_t(100.0) - 2506.015308
self.assertLess(error, self.maxError, 'uV_t not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_uL_t(self):
error = self.steamTable.uL_t(100.0) - 418.9933299
self.assertLess(error, self.maxError, 'uL_t not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_u_pt(self):
error = self.steamTable.u_pt(1.0, 100.0) - 2506.171426
self.assertLess(error, self.maxError, 'u_pt not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_u_ph(self):
error = self.steamTable.u_ph(1.0, 1000.0) - 956.2074342
self.assertLess(error, self.maxError, 'u_ph not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_u_ps(self):
error = self.steamTable.u_ps(1.0, 1.0) - 308.5082185
self.assertLess(error, self.maxError, 'u_ps not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_CpV_p(self):
error = self.steamTable.CpV_p(1.0) - 2.075938025
self.assertLess(error, self.maxError, 'cpV_p not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_CpL_p(self):
error = self.steamTable.CpL_p(1.0) - 4.216149431
self.assertLess(error, self.maxError, 'cpL_p not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_CpV_t(self):
error = self.steamTable.CpV_t(100.0) - 2.077491868
self.assertLess(error, self.maxError, 'cpV_t not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_CpL_t(self):
error = self.steamTable.CpL_t(100.0) - 4.216645119
self.assertLess(error, self.maxError, 'cpL_t not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_Cp_pt(self):
error = self.steamTable.Cp_pt(1.0, 100.0) - 2.074108555
self.assertLess(error, self.maxError, 'cp_pt not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_Cp_ph(self):
error = self.steamTable.Cp_ph(1.0, 200.0) - 4.17913573169
self.assertLess(error, self.maxError, 'Cp_ph not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_Cp_ps(self):
error = self.steamTable.Cp_ps(1.0, 1.0) - 4.190607038
self.assertLess(error, self.maxError, 'Cp_ps not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_CvV_p(self):
error = self.steamTable.CvV_p(1.0) - 1.552696979
self.assertLess(error, self.maxError, 'CvV_p not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_CvL_p(self):
error = self.steamTable.CvL_p(1.0) - 3.769699683
self.assertLess(error, self.maxError, 'CvL_p not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_CvV_t(self):
error = self.steamTable.CvV_t(100.0) - 1.553698696
self.assertLess(error, self.maxError, 'CvV_t not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_CvL_t(self):
error = self.steamTable.CvL_t(100.0) - 3.76770022
self.assertLess(error, self.maxError, 'CvL_t not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_Cv_pt(self):
error = self.steamTable.Cv_pt(1.0, 100.0) - 1.551397249
self.assertLess(error, self.maxError, 'Cv_pt not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_Cv_ph(self):
error = self.steamTable.Cv_ph(1.0, 200.0) - 4.035176364
self.assertLess(error, self.maxError, 'Cv_ph not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_Cv_ps(self):
error = self.steamTable.Cv_ps(1.0, 1.0) - 3.902919468
self.assertLess(error, self.maxError, 'Cv_ps not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_wV_p(self):
error = self.steamTable.wV_p(1.0) - 472.0541571
self.assertLess(error, self.maxError, 'wV_p not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_wL_p(self):
error = self.steamTable.wL_p(1.0) - 1545.451948
self.assertLess(error, self.maxError, 'wL_p not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_wV_t(self):
error = self.steamTable.wV_t(100.0) - 472.2559492
self.assertLess(error, self.maxError, 'wV_t not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_wL_t(self):
error = self.steamTable.wL_t(100.0) - 1545.092249
self.assertLess(error, self.maxError, 'wL_t not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_w_pt(self):
error = self.steamTable.w_pt(1.0, 100.0) - 472.3375235
self.assertLess(error, self.maxError, 'w_pt not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_w_ph(self):
error = self.steamTable.w_ph(1.0, 200.0) - 1542.682475
self.assertLess(error, self.maxError, 'w_ph not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_w_ps(self): # TODO: Check values and calculation
error = self.steamTable.w_ps(1.0, 1.0) - 1557.8585
self.assertLess(error, self.maxError, 'w_ps not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_my_pt(self):
error = self.steamTable.my_pt(1.0, 100.0) - 1.22704e-05
self.assertLess(error, self.maxError, 'my_pt not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_my_ph(self):
error = self.steamTable.my_ph(1.0, 100.0) - 0.000914003770302
self.assertLess(error, self.maxError, 'my_ph not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_my_ps(self):
error = self.steamTable.my_ps(1.0, 1.0) - 0.000384222
self.assertLess(error, self.maxError, 'my_ps not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_tcL_p(self):
error = self.steamTable.tcL_p(1.0) - 0.677593822
self.assertLess(error, self.maxError, 'tcL_p not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_tcV_p(self):
error = self.steamTable.tcV_p(1.0) - 0.024753668
self.assertLess(error, self.maxError, 'tcV_p not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_tcL_t(self):
error = self.steamTable.tcL_t(25.0) - 0.607458162
self.assertLess(error, self.maxError, 'tcL_t not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_tcV_t(self):
error = self.steamTable.tcV_t(25.0) - 0.018326723
self.assertLess(error, self.maxError, 'tcV_t not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_tc_pt(self):
error = self.steamTable.tc_pt(1.0, 25.0) - 0.607509806
self.assertLess(error, self.maxError, 'tc_pt not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_tc_ph(self):
error = self.steamTable.tc_ph(1.0, 100.0) - 0.605710062
self.assertLess(error, self.maxError, 'tc_ph not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_tc_hs(self):
error = self.steamTable.tc_hs(100.0, 0.34) - 0.606283124
self.assertLess(error, self.maxError, 'tc_hs not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_st_t(self):
error = self.steamTable.st_t(100.0) - 0.0589118685877
self.assertLess(error, self.maxError, 'st_t not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_st_p(self):
error = self.steamTable.st_p(1.0) - 0.058987784
self.assertLess(error, self.maxError, 'st_p not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_x_ph(self):
error = self.steamTable.x_ph(1.0, 1000.0) - 0.258055424
self.assertLess(error, self.maxError, 'x_ph not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_x_ps(self):
error = self.steamTable.x_ps(1.0, 4.0) - 0.445397961
self.assertLess(error, self.maxError, 'x_ps not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_vx_ph(self):
error = self.steamTable.vx_ph(1.0, 418.0) - 0.288493093
self.assertLess(error, self.maxError, 'vx_ph not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
def test_vx_ps(self):
error = self.steamTable.vx_ps(1.0, 4.0) - 0.999233827
self.assertLess(error, self.maxError, 'vx_ps not passed Error %(error)e allowed: %(max)e' % {'error': error, 'max': self.maxError})
print(f"p2dash: {round(float(p2),1)} bar")
print(f"H2dash: {round(float(h2dash),1)} kJ/kg")
print(f"S2dash: {round(float(s2dash),3)} kJ/kg K")
h3dash = steamTable.hV_p(p2)
s3dash = steamTable.sV_p(p2)
T3dash = T2dash
print('\nPoint 3dash')
print(f"T3dash: {round(float(T3dash),1)} degC")
print(f"H3dash: {round(float(h3dash),1)} kJ/kg")
print(f"S3dash: {round(float(s3dash),3)} kJ/kg K")
p3 = p2
T3 = 540
h3 = steamTable.h_pt(p3, T3)
s3 = steamTable.s_pt(p3, T3)
print('\nPoint 3')
print(f"T3: {round(float(T3),1)} degC")
print(f"p3: {round(float(p3),1)} bar")
print(f"H3: {round(float(h3),1)} kJ/kg")
print(f"S3: {round(float(s3),3)} kJ/kg K")
p4 = p1
s4 = s3
T4 = steamTable.t_ps(p4, s4)
x4 = steamTable.x_ps(p4, s4)
h4 = steamTable.h_px(p4, x4)
print('\nPoint 4')
print(f"T4: {round(float(T4),1)} degC")
print(f"P4: {round(float(p4),1)} bar")
class FLS_FunctionTester(unittest.TestCase):
def setUp(self):
self.maxError = 1E-6
self.steamTable = XSteam(XSteam.UNIT_SYSTEM_FLS)
def tearDown(self):
pass
def test_tsat_p(self):
error = self.steamTable.tsat_p(1.0) - 99.60591861
self.assertLess(
error, self.maxError,
'tsat_p not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_t_ph(self):
error = self.steamTable.t_ph(1.0, 100.0) - 23.84481908
self.assertLess(
error, self.maxError,
't_ph not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_t_ps(self):
error = self.steamTable.t_ps(1.0, 1.0) - 73.70859421
self.assertLess(
error, self.maxError,
't_ps not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_t_hs(self):
error = self.steamTable.t_hs(100.0, 0.2) - 13.84933511
self.assertLess(
error, self.maxError,
't_hs not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_psat_t(self):
error = self.steamTable.psat_t(100.0) - 1.014179779
self.assertLess(
error, self.maxError,
'psat_t not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_p_hs(self):
error = self.steamTable.p_hs(84.0, 0.296) - 2.295498269
self.assertLess(
error, self.maxError,
'p_hs not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_hV_p(self):
error = self.steamTable.hV_p(1.0) - 2674.949641
self.assertLess(
error, self.maxError,
'hV_p not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_hL_p(self):
error = self.steamTable.hL_p(1.0) - 417.4364858
self.assertLess(
error, self.maxError,
'hL_p not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_hV_t(self):
error = self.steamTable.hV_t(100.0) - 2675.572029
self.assertLess(
error, self.maxError,
'hV_t not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_hL_t(self):
error = self.steamTable.hL_t(100.0) - 419.099155
self.assertLess(
error, self.maxError,
'hL_t not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_h_pt(self):
error = self.steamTable.h_pt(1.0, 20.0) - 84.01181117
self.assertLess(
error, self.maxError,
'h_pt not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_h_ps(self):
error = self.steamTable.h_ps(1.0, 1.0) - 308.6107171
self.assertLess(
error, self.maxError,
'h_ps not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_h_px(self):
error = self.steamTable.h_px(1.0, 0.5) - 1546.193063
self.assertLess(
error, self.maxError,
'h_px not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_h_prho(self):
error = self.steamTable.h_prho(1.0, 2.0) - 1082.773391
self.assertLess(
error, self.maxError,
'h_prho not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_h_tx(self):
error = self.steamTable.h_tx(100.0, 0.5) - 1547.33559211
self.assertLess(
error, self.maxError,
'h_tx not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_vV_p(self):
error = self.steamTable.vV_p(1.0) - 1.694022523
self.assertLess(
error, self.maxError,
'vV_p not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_vL_p(self):
error = self.steamTable.vL_p(1.0) - 0.001043148
self.assertLess(
error, self.maxError,
'vL_p not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_vV_t(self):
error = self.steamTable.vV_t(100.0) - 1.671860601
self.assertLess(
error, self.maxError,
'vV_t not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_vL_t(self):
error = self.steamTable.vL_t(100.0) - 0.001043455
self.assertLess(
error, self.maxError,
'vL_t not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_v_pt(self):
error = self.steamTable.v_pt(1.0, 100.0) - 1.695959407
self.assertLess(
error, self.maxError,
'v_pt not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_v_ph(self):
error = self.steamTable.v_ph(1.0, 1000.0) - 0.437925658
self.assertLess(
error, self.maxError,
'v_ph not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_v_ps(self):
error = self.steamTable.v_ps(1.0, 5.0) - 1.03463539
self.assertLess(
error, self.maxError,
'v_ps not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_rhoV_p(self):
error = self.steamTable.rhoV_p(1.0) - 0.590310924
self.assertLess(
error, self.maxError,
'rhoV_p not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_rhoL_p(self):
error = self.steamTable.rhoL_p(1.0) - 958.6368897
self.assertLess(
error, self.maxError,
'rhoL_p not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_rhoV_t(self):
error = self.steamTable.rhoV_t(100.0) - 0.598135993
self.assertLess(
error, self.maxError,
'rhoV_t not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_rhoL_t(self):
error = self.steamTable.rhoL_t(100.0) - 958.3542773
self.assertLess(
error, self.maxError,
'rhoL_t not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_rho_pt(self):
error = self.steamTable.rho_pt(1.0, 100.0) - 0.589636754
self.assertLess(
error, self.maxError,
'rho_pt not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_rho_ph(self):
error = self.steamTable.rho_ph(1.0, 1000.0) - 2.283492601
self.assertLess(
error, self.maxError,
'rho_ph not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_rho_ps(self):
error = self.steamTable.rho_ps(1.0, 1.0) - 975.6236788
self.assertLess(
error, self.maxError,
'rho_ps not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_sV_p(self):
error = self.steamTable.sV_p(0.006117) - 9.155465556
self.assertLess(
error, self.maxError,
'sV_p not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_sL_p(self):
error = self.steamTable.sL_p(0.0061171) - 1.8359e-05
self.assertLess(
error, self.maxError,
'sL_p not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_sV_t(self):
error = self.steamTable.sV_t(0.0001) - 9.155756716
self.assertLess(
error, self.maxError,
'sV_t not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_sL_t(self):
error = self.steamTable.sL_t(100.0) - 1.307014328
self.assertLess(
error, self.maxError,
'sL_t not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_s_pt(self):
error = self.steamTable.s_pt(1.0, 20.0) - 0.296482921
self.assertLess(
error, self.maxError,
's_pt not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_s_ph(self):
error = self.steamTable.s_ph(1.0, 84.01181117) - 0.296813845
self.assertLess(
error, self.maxError,
's_ph not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_uV_p(self):
error = self.steamTable.uV_p(1.0) - 2505.547389
self.assertLess(
error, self.maxError,
'uV_p not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_uL_p(self):
error = self.steamTable.uL_p(1.0) - 417.332171
self.assertLess(
error, self.maxError,
'uL_p not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_uV_t(self):
error = self.steamTable.uV_t(100.0) - 2506.015308
self.assertLess(
error, self.maxError,
'uV_t not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_uL_t(self):
error = self.steamTable.uL_t(100.0) - 418.9933299
self.assertLess(
error, self.maxError,
'uL_t not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_u_pt(self):
error = self.steamTable.u_pt(1.0, 100.0) - 2506.171426
self.assertLess(
error, self.maxError,
'u_pt not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_u_ph(self):
error = self.steamTable.u_ph(1.0, 1000.0) - 956.2074342
self.assertLess(
error, self.maxError,
'u_ph not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_u_ps(self):
error = self.steamTable.u_ps(1.0, 1.0) - 308.5082185
self.assertLess(
error, self.maxError,
'u_ps not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_CpV_p(self):
error = self.steamTable.cpV_p(1.0) - 2.075938025
self.assertLess(
error, self.maxError,
'cpV_p not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_CpL_p(self):
error = self.steamTable.cpL_p(1.0) - 4.216149431
self.assertLess(
error, self.maxError,
'cpL_p not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_CpV_t(self):
error = self.steamTable.cpV_t(100.0) - 2.077491868
self.assertLess(
error, self.maxError,
'cpV_t not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_CpL_t(self):
error = self.steamTable.cpL_t(100.0) - 4.216645119
self.assertLess(
error, self.maxError,
'cpL_t not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_Cp_pt(self):
error = self.steamTable.cp_pt(1.0, 100.0) - 2.074108555
self.assertLess(
error, self.maxError,
'cp_pt not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_Cp_ph(self):
error = self.steamTable.Cp_ph(1.0, 200.0) - 4.17913573169
self.assertLess(
error, self.maxError,
'Cp_ph not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_Cp_ps(self):
error = self.steamTable.Cp_ps(1.0, 1.0) - 4.190607038
self.assertLess(
error, self.maxError,
'Cp_ps not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_CvV_p(self):
error = self.steamTable.CvV_p(1.0) - 1.552696979
self.assertLess(
error, self.maxError,
'CvV_p not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_CvL_p(self):
error = self.steamTable.CvL_p(1.0) - 3.769699683
self.assertLess(
error, self.maxError,
'CvL_p not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_CvV_t(self):
error = self.steamTable.CvV_t(100.0) - 1.553698696
self.assertLess(
error, self.maxError,
'CvV_t not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_CvL_t(self):
error = self.steamTable.CvL_t(100.0) - 3.76770022
self.assertLess(
error, self.maxError,
'CvL_t not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_Cv_pt(self):
error = self.steamTable.Cv_pt(1.0, 100.0) - 1.551397249
self.assertLess(
error, self.maxError,
'Cv_pt not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_Cv_ph(self):
error = self.steamTable.Cv_ph(1.0, 200.0) - 4.035176364
self.assertLess(
error, self.maxError,
'Cv_ph not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_Cv_ps(self):
error = self.steamTable.Cv_ps(1.0, 1.0) - 3.902919468
self.assertLess(
error, self.maxError,
'Cv_ps not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_wV_p(self):
error = self.steamTable.wV_p(1.0) - 472.0541571
self.assertLess(
error, self.maxError,
'wV_p not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_wL_p(self):
error = self.steamTable.wL_p(1.0) - 1545.451948
self.assertLess(
error, self.maxError,
'wL_p not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_wV_t(self):
error = self.steamTable.wV_t(100.0) - 472.2559492
self.assertLess(
error, self.maxError,
'wV_t not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_wL_t(self):
error = self.steamTable.wL_t(100.0) - 1545.092249
self.assertLess(
error, self.maxError,
'wL_t not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_w_pt(self):
error = self.steamTable.w_pt(1.0, 100.0) - 472.3375235
self.assertLess(
error, self.maxError,
'w_pt not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_w_ph(self):
error = self.steamTable.w_ph(1.0, 200.0) - 1542.682475
self.assertLess(
error, self.maxError,
'w_ph not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_w_ps(self):
error = self.steamTable.w_ps(1.0, 1.0) - 1557.8585
self.assertLess(
error, self.maxError,
'w_ps not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_my_pt(self):
error = self.steamTable.my_pt(1.0, 100.0) - 1.22704e-05
self.assertLess(
error, self.maxError,
'my_pt not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_my_ph(self):
error = self.steamTable.my_ph(1.0, 100.0) - 0.000914003770302
self.assertLess(
error, self.maxError,
'my_ph not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_my_ps(self):
error = self.steamTable.my_ps(1.0, 1.0) - 0.000384222
self.assertLess(
error, self.maxError,
'my_ps not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_tcL_p(self):
error = self.steamTable.tcL_p(1.0) - 0.677593822
self.assertLess(
error, self.maxError,
'tcL_p not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_tcV_p(self):
error = self.steamTable.tcV_p(1.0) - 0.024753668
self.assertLess(
error, self.maxError,
'tcV_p not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_tcL_t(self):
error = self.steamTable.tcL_t(25.0) - 0.607458162
self.assertLess(
error, self.maxError,
'tcL_t not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_tcV_t(self):
error = self.steamTable.tcV_t(25.0) - 0.018326723
self.assertLess(
error, self.maxError,
'tcV_t not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_tc_pt(self):
error = self.steamTable.tc_pt(1.0, 25.0) - 0.607509806
self.assertLess(
error, self.maxError,
'tc_pt not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_tc_ph(self):
error = self.steamTable.tc_ph(1.0, 100.0) - 0.605710062
self.assertLess(
error, self.maxError,
'tc_ph not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_tc_hs(self):
error = self.steamTable.tc_hs(100.0, 0.34) - 0.606283124
self.assertLess(
error, self.maxError,
'tc_hs not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_st_t(self):
error = self.steamTable.st_t(100.0) - 0.0589118685877
self.assertLess(
error, self.maxError,
'st_t not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_st_p(self):
error = self.steamTable.st_p(1.0) - 0.058987784
self.assertLess(
error, self.maxError,
'st_p not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_x_ph(self):
error = self.steamTable.x_ph(1.0, 1000.0) - 0.258055424
self.assertLess(
error, self.maxError,
'x_ph not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_x_ps(self):
error = self.steamTable.x_ps(1.0, 4.0) - 0.445397961
self.assertLess(
error, self.maxError,
'x_ps not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_vx_ph(self):
error = self.steamTable.vx_ph(1.0, 418.0) - 0.288493093
self.assertLess(
error, self.maxError,
'vx_ph not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
def test_vx_ps(self):
error = self.steamTable.vx_ps(1.0, 4.0) - 0.999233827
self.assertLess(
error, self.maxError,
'vx_ps not passed Error %(error)e allowed: %(max)e' % {
'error': error,
'max': self.maxError
})
print(f"P2dash: {round(float(p2),1)} bar")
print(f"H2dash: {round(float(h2dash),1)} kJ/kg")
print(f"S2dash: {round(float(s2dash),3)} kJ/kg K")
h3dash = steamTable.hV_p(p2)
s3dash = steamTable.sV_p(p2)
T3dash = T2dash
print('\nPoint 3dash')
print(f"T3dash: {round(float(T3dash),1)} degC")
print(f"H3dash: {round(float(h3dash),1)} kJ/kg")
print(f"S3dash: {round(float(s3dash),3)} kJ/kg K")
p3 = p2
T3 = 540
h3 = steamTable.h_pt(p3, T3)
s3 = steamTable.s_pt(p3, T3)
print('\nPoint 3')
print(f"T3: {round(float(T3),1)} degC")
print(f"P3: {round(float(p3),1)} bar")
print(f"H3: {round(float(h3),1)} kJ/kg")
print(f"S3: {round(float(s3),3)} kJ/kg K")
p4 = p2/8
print(f"Reheat Pressure: {round(float(p4),1)} bar")
s4 = s3
T4 = steamTable.t_ps(p4, s4)
h4 = steamTable.h_pt(p4, T4)
HPturbeff = 0.90 # HP turbine isentropic efficiency can be entered here (typically 0.85 - 0.95)
h4r = h3 - (HPturbeff * h3) + (HPturbeff * h4)
DeltaP_low = (P - atm) * 100000 # pa
DeltaP_hig = (max_press - P) * 100000 # pa
# pump_power = DeltaP * m_dot / rho
pump_power_low = DeltaP_low * m_dot_low / rho / nP
pump_power_low *= 10 ** -3 # kpa
pump_powers_low.append(pump_power_low)
pump_power_hig = DeltaP_hig * m_dot_hig / rho / nP
pump_power_hig *= 10 ** -3 # kpa
pump_powers_hig.append(pump_power_hig)
pump_power = pump_power_low + pump_power_hig
pump_powers.append(pump_power)
s1 = steamTable.s_pt(atm, T)
h1 = steamTable.h_pt(atm, T)
h2 = h1 + pump_power_low / m_dot_low * nP
# h2 = steamTable.h_ps(P, s1)
h4 = h3 + pump_power_hig / m_dot_hig * nP
Q_core_low = m_dot_low * (h3 - h2) * 10 ** -3 # mw
core_heats_low.append(Q_core_low)
Q_core_hig = m_dot_hig * (h6 - h4) * 10 ** -3 #mw
core_heats_hig.append(Q_core_hig)
Q_core = Q_core_hig + Q_core_low
core_heats.append(Q_core)
efficiency = power / Q_core * 10 ** -3
eff.append(efficiency)
def demo_reheat_ms_ofwh():
"""
converted example from Stu Blair, found at https://github.com/stu314159/xsteam/blob/42767648a05c6759ad11aea95256cb24e4fc9499/Examples/Rankine_Reheat_MS_OFWH.m
"""
print("converted example from Stu Blair / https://github.com/stu314159")
print("Purpose: use fminsearch / scipy.optimize.fmin")
print(
"""Problem Description:
(a picture would be better....)
"""
)
steam_table = XSteam(XSteam.UNIT_SYSTEM_FLS)
# %% State point 1 - condenser exit
state_point_1 = dict()
state_point_1["P"] = 1.5 # % psia
state_point_1["h"] = steam_table.hL_p(state_point_1["P"])
state_point_1["s"] = steam_table.sL_p(state_point_1["P"])
# %% State point 1 -> 2, main condensate pump
eta_mcp = 0.84
state_point_2 = dict()
state_point_2["P"] = 164.0 # % psia
state_point_2["s_s"] = state_point_1["s"]
state_point_2["h_s"] = steam_table.h_ps(state_point_2["P"], state_point_2["s_s"])
state_point_2["h"] = (
state_point_1["h"] - (state_point_1["h"] - state_point_2["h_s"]) / eta_mcp
)
# %% State point 3, OFWH exit, saturated liquid
state_point_3 = dict()
state_point_3["P"] = state_point_2["P"]
state_point_3["h"] = steam_table.hL_p(state_point_3["P"])
state_point_3["s"] = steam_table.sL_p(state_point_3["P"])
# %% State point 3 -> 4, main feed pump
eta_mfp = 0.84
state_point_4 = dict()
state_point_4["P"] = 820.0 # % psia
state_point_4["s_s"] = state_point_3["s"]
state_point_4["h_s"] = steam_table.h_ps(state_point_4["P"], state_point_4["s_s"])
state_point_4["h"] = (
state_point_3["h"] - (state_point_3["h"] - state_point_4["h_s"]) / eta_mfp
)
# %% State point 5 - Steam generator exit
state_point_5 = dict()
state_point_5["P"] = state_point_4["P"]
state_point_5["h"] = steam_table.hV_p(state_point_5["P"])
state_point_5["s"] = steam_table.sV_p(state_point_5["P"])
# %% State point 6 - HP Turbine Exhaust
eta_hpt = 0.94
state_point_6 = dict()
state_point_6["P"] = 164.0 # % psia
state_point_6["s_s"] = state_point_5["s"]
state_point_6["h_s"] = steam_table.h_ps(state_point_6["P"], state_point_6["s_s"])
state_point_6["h"] = state_point_5["h"] - eta_hpt * (
state_point_5["h"] - state_point_6["h_s"]
)
state_point_6["x"] = steam_table.x_ph(state_point_6["P"], state_point_6["h"])
# %% State point 7 - Moisture Separator Exit
state_point_7 = dict()
state_point_7["P"] = state_point_6["P"]
state_point_7["h"] = steam_table.hV_p(state_point_7["P"])
state_point_7["s"] = steam_table.sV_p(state_point_7["P"])
# %% State point 8 - Reheater Mid-Pressure Steam exit
state_point_8 = dict()
state_point_8["P"] = state_point_7["P"]
state_point_8["T"] = 490.0 # % degrees F
state_point_8["h"] = steam_table.h_pt(state_point_8["P"], state_point_8["T"])
state_point_8["s"] = steam_table.s_pt(state_point_8["P"], state_point_8["T"])
# %% State point 9 - LP Turbine Exhaust
state_point_9 = dict()
state_point_9["P"] = state_point_1["P"]
eta_lpt = 0.94
state_point_9["s_s"] = state_point_8["s"]
state_point_9["h_s"] = steam_table.h_ps(state_point_9["P"], state_point_9["s_s"])
state_point_9["h"] = state_point_8["h"] - eta_lpt * (
state_point_8["h"] - state_point_9["h_s"]
)
# %% State point 10 - Reheater HP Steam exit
state_point_10 = dict()
state_point_10["P"] = state_point_5["P"]
# % assume steam exits as a saturated liquid.
state_point_10["h"] = steam_table.hL_p(state_point_10["P"])
# %% State point 11 - pressure trap exit to OFWH
state_point_11 = dict()
state_point_11["P"] = state_point_2["P"]
# % assume isenthalpic expansion in the trap.
state_point_11["h"] = state_point_10["h"]
# %% State point 12 - Moisture Separator liquid drain to OFWH
state_point_12 = dict()
state_point_12["P"] = state_point_6["P"]
state_point_12["h"] = steam_table.hL_p(state_point_12["P"])
# %% Heat Balance - find the flow fractions
def RH_heatBalance(f):
return (
f[0] * state_point_5["h"]
+ (1 - f[0]) * state_point_6["x"] * (1 - f[1]) * state_point_7["h"]
) - (
f[0] * state_point_10["h"]
+ (1 - f[0]) * state_point_6["x"] * (1 - f[1]) * state_point_8["h"]
)
def OFWH_heatBalance(f):
return (
(1 - f[0]) * (1 - f[1]) * state_point_6["x"] * state_point_2["h"]
+ f[0] * state_point_11["h"]
+ (1 - f[0]) * state_point_6["x"] * f[1] * state_point_7["h"]
+ (1 - f[0]) * (1 - state_point_6["x"]) * state_point_12["h"]
) - state_point_3["h"]
def totalFunctional(f):
return math.fabs(RH_heatBalance(f)) + math.fabs(OFWH_heatBalance(f))
# % the strategy is to minimize the total functional. The minimum value is
# % when they are both equal to zero.
initialGuess = [0.1, 0.1]
f = fmin(func=totalFunctional, x0=initialGuess)
# % % calculate heat and energy balances
w_mcp = (
(state_point_1["h"] - state_point_2["h"])
* (1 - f[0])
* (1 - f[1])
* state_point_6["x"]
)
w_mfp = state_point_3["h"] - state_point_4["h"]
w_hpt = (state_point_5["h"] - state_point_6["h"]) * (1 - f[0])
w_lpt = (
(state_point_8["h"] - state_point_9["h"])
* (1 - f[0])
* (1 - f[1])
* state_point_6["x"]
)
w_net = w_mcp + w_mfp + w_hpt + w_lpt
q_cond = (
(state_point_1["h"] - state_point_9["h"])
* (1 - f[0])
* (1 - f[1])
* state_point_6["x"]
)
q_sg = state_point_5["h"] - state_point_4["h"]
q_net = q_cond + q_sg
eta_th = w_net / q_sg
# % % report the results:
print(f"Net heat: {q_net:.3f} BTU/lbm Net work: {w_net:.3f} BTU/lbm")
print(f"Thermal efficiency: {eta_th*100:.2f} percent")
s2 = s1
v = 1/steamTable.rhoL_p(p1)
w_p = v*(p2-p1)
print('\nPoint 2')
h2 = h1+w_p
print(f"H2: {round(float(h2),1)} kJ/kg")
T2 = steamTable.t_ph(p2, h2)
print(f"T2: {round(float(T2),1)} degC")
p3 = p2
T3 = 600
h3 = steamTable.h_pt(p3, T3)
s3 = steamTable.s_pt(p3, T3)
print('\nPoint 3')
print(f"T3: {round(float(T3),1)} degC")
print(f"P3: {round(float(p3),1)} bar")
print(f"H3: {round(float(h3),1)} kJ/kg")
print(f"S3: {round(float(s3),3)} kJ/kg K")
p4 = 80
print(f"Reheat 1 Pressure: {round(float(p4),1)} bar")
s4 = s3
T4 = steamTable.t_ps(p4, s4)
h4 = steamTable.h_pt(p4, T4)
print('\nPoint 4')
print(f"T4: {round(float(T4),1)} degC")
print(f"P4: {round(float(p4),1)} bar")
print(f"H4: {round(float(h4),1)} kJ/kg")
