from pyXSteam.XSteam import XSteam
import matplotlib.pyplot as plt

steamTable = XSteam(XSteam.UNIT_SYSTEM_BARE)  # m/kg/sec/K/MPa/W

# input:
# dm_eth: massflow of ethanol gas from dest column [m^3/h] #10
# h_eth1: enthalpy of ethanol gas from dest column [kJ/kg]
# dm_c: mass flow cooling agent [kg/s]
# h_c1: enthalpy cooling agent [kJ/kg]

# output:
# h_eth2: enthalpy ethanol gas to B180 [kJ/kg]
# h_c2: enthalpy cooling agent [kJ/kg]
# Q_l: amount of energy transferred from ethanol gas to cooling agent

dm_eth = 11.37  # [m^3/h]
dm_h2o = 18  # [kg/h]
T1 = 74.3  # [C]
T2 = 13.0  # [C]
P2 = 0.3  # MPa


def heat_exchange():
    # density of ethanol
    den = 785.3  # kg/m^3

    # specific heat capacity {isobaric gas} should be dynamic
    c_eth = 2.18  # [kJ/kg K]

    # heat transferred to cooling agent
    q_1 = dm_eth * c_eth * (T1 - T2) * den

    # entropy of water at atmospheric pressure
    s2 = steamTable.sV_p(0.101325)

    # entropy after compression given ideal conditions
    s3 = s2

    # enthalpy entering heat exchanger
    h3 = steamTable.h_ps(P2, s3)

    return q_1, h3


print(heat_exchange())