import sys sys.path.append("/Users/jonathanxavier/Developer/sim42") from ollin.Administrator.AdmOllin import Ollin from pylab import * RK = Ollin.AddModel("RK", "RK") S1 = Ollin.AddCase("S1") Ollin.Add(["N-BUTANE", "N-PENTANE", "N-HEXANE"], "RK") S1.SetX([1, 1, 1]) S1.T(322) plot_x = range(10, 1000, 10) plot_y0 = [] plot_y1 = [] plot_y2 = [] plot_y3 = [] plot_y4 = [] ##plot_y5 = [] Pr = [] print plot_x for P in plot_x: S1.P(P) Ollin.Solve() yf = S1.Get("Ki") plot_y0.append(yf[0]) plot_y1.append(yf[1]) plot_y2.append(yf[2]) ## plot_y3.append( yf[3] ) ## plot_y4.append( yf[4] )
import sys sys.path.append("/Users/jonathanxavier/Developer/sim42") from ollin.Administrator.AdmOllin import Ollin from pylab import * RK = Ollin.AddModel("RK", "PR") S1 = Ollin.AddCase("S1") Ollin.Add(["ETHANE", "PROPANE", "N-BUTANE", "N-PENTANE", "N-HEXANE"], "RK") S1.SetX([0.05, 0.15, 0.25, 0.20, 0.35]) S1.P(101.325) plot_x = range(240, 325, 2) plot_y0 = [] plot_y1 = [] plot_y2 = [] plot_y3 = [] plot_y4 = [] plot_y5 = [] for T in plot_x: S1.T(T) Ollin.Solve() yf = S1.Get("yf") plot_y0.append(yf[0]) plot_y1.append(yf[1]) plot_y2.append(yf[2]) plot_y3.append(yf[3]) plot_y4.append(yf[4])
import sys sys.path.append("/Users/jonathanxavier/Developer/sim42") from ollin.Administrator.AdmOllin import Ollin from ollin.Example.UOS import Stream, Valve,Heater,Flash PR=Ollin.AddModel("PR","PR") Ollin.Add(["METHANE","ETHANE","PROPANE","ISOBUTANE","N-BUTANE","N-PENTANE","N-HEXANE","N-HEPTANE"],"PR") Ollin.LoadConst() S1 = Stream(PR) S1.X([14520,9070,7260,770,2810,950+1630,1540,3180]) ##S1.FracVap(0.9) S1.T(273.15+40) S1.P(500) S1.Mol(455) V1 = Valve(PR) V1.DP = 11030 H1 = Heater(PR) H1.DP = 1 H1.DH = 1 F1 = Flash(PR) S1.Connect(F1) ## ##V1.Connect(H1) ##H1.Connect(F1)
from ollin.Administrator.AdmOllin import Ollin from numpy.oldnumeric import array RK = Ollin.AddModel("PR", "PR", "ANTOINE") Ollin.Add([ "METHANOL", "HYDROGEN", "CARBON MONOXIDE", "CARBON DIOXIDE", "WATER", "METHANE" ], "PR") S1 = Ollin.AddCase("S1") #Alimentacion al compresor S1.SetX([0, 69, 21, 7, 0, 3]) S1.T(25 + 273.15) S1.P(101.325 * 1) Ollin.Solve("S1") #Caculamos S1 Ollin.Resumen("S1") ##print S1.Get("H") S2 = Ollin.AddCase("S2") #salida del compresor S2.SetX(S1.Get("x")) S2.H(S1.Get("H")) #Considerando que es adiabatico dQ=0 S2.P(101.325 * 20) Ollin.Solve("S2") #Caculamos S2 Ollin.Resumen("S2") #Salida del calentador S3 = Ollin.AddCase("S3") #Alimentacion al ractor S3.SetX(S2.Get("x")) S3.T(400 + 273.15) S3.P(S2.Get("P")) Ollin.Solve("S3") #La reaccion se realiza en fase gas adiabaticamente dQ=0