Example #1
0
def EVOutVaporComposition(m, t):
    return m.EVOutVPMolarFrac[t, 'Oxygen'] == vapor_oxygen_composition(
        m.ircRebPressure[t], m.ircEVOutTemperature[t])
Example #2
0
def ColumnVaporFraction(m, t, tray):
    return m.vaporLeavingMolarFraction[t, tray,
                                       'Oxygen'] == vapor_oxygen_composition(
                                           m.trayPressure[t, tray],
                                           m.trayTemperature[t, tray])
Example #3
0
def rebVaporComposition(m, t):
    return m.ircWNMolarFraction[t, 'Oxygen'] == vapor_oxygen_composition(
        m.ircRebPressure[t], m.ircRebTemperature[t])
Example #4
0
def FeedVaporComposition(m):
    return m.FeedVaporPhaseMolarFrac['Oxygen'] == vapor_oxygen_composition(
        m.FeedPressure, m.FeedTemperature)
Example #5
0
def FeedVaporComposition(m, t):
    return (m.FeedVaporPhaseMolarFrac[t, 'Oxygen'] - vapor_oxygen_composition(
        m.FeedPressure[t], m.FeedTemperature[t])) * 1e2 == 0
Example #6
0
def ColumnAllTraysVapFrac(m,tray):
	return m.ColumnVapLvMFrac[tray,"Oxygen"] == vapor_oxygen_composition(m.ColumnTrayPressure[tray],m.ColumnTrayTemp[tray])
Example #7
0
def ReboilerHdlpVapFrac(m):
	return m.ReboilerVapLvMFrac["Oxygen"] == vapor_oxygen_composition(m.ColumnTrayPressure[1],m.ReboilerTemp)
Example #8
0
def FeedThermCompVap(m):
    return m.FeedVapMFrac["Oxygen"] == vapor_oxygen_composition(
        m.FeedPressure, m.FeedTemp)
Example #9
0
def FeedSelfVapFrac(m):
	return m.FeedVapMFrac["Oxygen"] == vapor_oxygen_composition(m.FeedPressure,m.FeedTemp)
Example #10
0
def VaporPhaseEnthalpy(temperature, pressure, oxygen, nitrogen):
    return vapor_nitrogen_enthalpy(
        pressure, temperature) * nitrogen + vapor_oxygen_enthalpy(
            pressure, temperature) * oxygen


temp_profile = []
enthalpy_profile = []
vapor_frac = []
for ircRebTemperature in np.arange(-190, -186, 0.1):
    temp_profile.append(ircRebTemperature)
    ircRebPressure = 564 - 400
    total_nitrogen = 0.7009695483157175
    total_oxygen = 0.29903045168428255
    EVOutVPMolarFrac_oxygen = vapor_oxygen_composition(ircRebPressure,
                                                       ircRebTemperature)
    EVOutLPMolarFrac_oxygen = liquid_oxygen_composition(
        ircRebPressure, ircRebTemperature)
    EVOutVaporFraction = (EVOutLPMolarFrac_oxygen - total_oxygen) / (
        EVOutLPMolarFrac_oxygen - EVOutVPMolarFrac_oxygen)
    enthalpy =  LiquidPhaseEnthalpy(ircRebTemperature, ircRebPressure, EVOutLPMolarFrac_oxygen,
                                  1-EVOutLPMolarFrac_oxygen)*(1-EVOutVaporFraction)+\
                VaporPhaseEnthalpy(ircRebTemperature, ircRebPressure, EVOutVPMolarFrac_oxygen,
                                   1-EVOutVPMolarFrac_oxygen)*EVOutVaporFraction
    enthalpy_profile.append(enthalpy)
    vapor_frac.append(EVOutVaporFraction)

fig, axes = plt.subplots(1, 2)
axes[0].plot(temp_profile, enthalpy_profile, c='r')
axes[1].plot(temp_profile, vapor_frac, c='b')
plt.show()
Example #11
0
def ReboilerThermCompVap(m, time):
    return m.ReboilerVapLvMFrac[time, "Oxygen"] == vapor_oxygen_composition(
        m.ColumnTrayPressure[time, 1], m.ReboilerTemp[time])
Example #12
0
def ColumnThermCompVap(m, time, tray):
    return m.ColumnVapLvMFrac[time, tray,
                              "Oxygen"] == vapor_oxygen_composition(
                                  m.ColumnTrayPressure[time, tray],
                                  m.ColumnTrayTemp[time, tray])
Example #13
0
def FeedThermCompVap(m, time):
    return m.FeedVapMFrac[time, "Oxygen"] == vapor_oxygen_composition(
        m.FeedPressure[time], m.FeedTemp[time])
Example #14
0
def EVOutVaporComposition(m):
    return m.EVOutVPMolarFrac['Oxygen'] == vapor_oxygen_composition(
        m.ircRebPressure, m.ircRebTemperature)