def trayHoldupEnthalpyPartialX(m, t, tray, c):
if c == 'Oxygen':
return liquid_oxygen_enthalpy(m.trayPressure[t, tray],
m.trayTemperature[t, tray])
elif c == 'Nitrogen':
return liquid_nitrogen_enthalpy(m.trayPressure[t, tray],
m.trayTemperature[t, tray])
def ColumnPxTrayMEtlpHldp(m, time, tray, comp):
if comp == 'Oxygen':
return liquid_oxygen_enthalpy(m.ColumnTrayPressure[time, tray],
m.ColumnTrayTemp[time, tray])
elif comp == 'Nitrogen':
return liquid_nitrogen_enthalpy(m.ColumnTrayPressure[time, tray],
m.ColumnTrayTemp[time, tray])
def LiquidPhaseEnthalpy(temperature, pressure, oxygen, nitrogen):
return liquid_nitrogen_enthalpy(
pressure, temperature) * nitrogen + liquid_oxygen_enthalpy(
pressure, temperature) * oxygen