B FINAL.sagews

# PARTE B

y = vector([0.2317,0.22511,0.54316]) #Fraciones para cada compuesto
Psat = vector([12.455,17.297, 75.460]) #Presiones saturadas  para cada compuesto
S = 0.32028 #Corriente de los moles de vapor encontrados en la seccion anterior
z1=0.4
z2=0.3
z3=0.3

print(" ")
print("PARTE b .RECICLO")

r1=y[0]*0.95*S

r2=y[1]*0.95*S

r3=y[2]*0.95*S


m4=r1+z1

m5=r2+z2

m6=r3+z3

mt= m4+m5+m6

z4=m4/mt

z5=m5/mt

z6=m6/mt

Pbur=Psat[0]*z4+Psat[1]*z5+Psat[2]*z6

Procio=1/((z4/Psat[0])+(z5/Psat[1])+(z6/Psat[2]))

Pprom=(Pbur+Procio)/2

k =[Psat[i]/Pprom for i in range(3)]

#Cálculo y expresión de V
V=0
fx=(z4*k[0]/(1+x*(k[0]-1)))+(z5*k[1]/(1+x*(k[1]-1)))+(z6*k[2]/(1+x*(k[2]-1)))-1
V=find_root(fx, 0.05, 0.99)

#Cálculo nuevos y

y1=((z4*k[0])/(1+(V*(k[0]-1))))

y2=((z5*k[1])/(1+(V*(k[1]-1))))

y3=((z6*k[2])/(1+(V*(k[2]-1))))

suma=y1+y2+y3

#Cálculo y expresión de x

x1=y1/k[0]

x2=y2/k[1]

x3=y3/k[2]

suma1=x1+x2+x3

print("Corriente V")
show("V= ", V)

print("Corriente L")
show("L= ", ((1-0.05*S)))

print("Presión de operación")
show("P= ", Pprom)


print("Valores de composición gaseosa segunda pasada, y")

show("y1= ", y1)

show("y2= ", y2)

show("y3= ", y3)

show("Suma= ", suma)


print("Valores de composición líquida segunda pasada, x")

show("x1= ", x1)

show("x2= ", x2)

show("x3= ", x3)

show("Suma= ", suma1)

PARTE b .RECICLO Corriente V

V= $\displaystyle 0.352935533282$

Corriente L

L= $\displaystyle 0.983986000000000$

Presión de operación

P= $\displaystyle 28.1443352651515$

Valores de composición gaseosa segunda pasada, y

y1= $\displaystyle 0.198744732195404$

y2= $\displaystyle 0.200977051932366$

y3= $\displaystyle 0.600278215872231$

Suma= $\displaystyle 1.00000000000000$

Valores de composición líquida segunda pasada, x

x1= $\displaystyle 0.449099829393029$

x2= $\displaystyle 0.327014252771362$

x3= $\displaystyle 0.223885917835609$

Suma= $\displaystyle 1.00000000000000$