Prepara tus exámenes y mejora tus resultados gracias a la gran cantidad de recursos disponibles en Docsity
Gana puntos ayudando a otros estudiantes o consíguelos activando un Plan Premium
Prepara tus exámenes
Prepara tus exámenes y mejora tus resultados gracias a la gran cantidad de recursos disponibles en Docsity
Prepara tus exámenes con los documentos que comparten otros estudiantes como tú en Docsity
Encuentra los documentos específicos para los exámenes de tu universidad
Estudia con lecciones y exámenes resueltos basados en los programas académicos de las mejores universidades
Responde a preguntas de exámenes reales y pon a prueba tu preparación
Consigue puntos base para descargar
Gana puntos ayudando a otros estudiantes o consíguelos activando un Plan Premium
Comunidad
Pide ayuda a la comunidad y resuelve tus dudas de estudio
Ebooks gratuitos
Descarga nuestras guías gratuitas sobre técnicas de estudio, métodos para controlar la ansiedad y consejos para la tesis preparadas por los tutores de Docsity
Programa para poder calculara una destilacion flash en lenguaje python
Tipo: Ejercicios
1 / 12
Esta página no es visible en la vista previa
¡No te pierdas las partes importantes!
En oferta
V = V - ((d-1)/fs) i = 0 print("{:^10} {:^1.4f}".format("V=", V)) file.write(' ') file.write('V = % s' % V + '\n') FL = F-V print("{:^50}".format("Concentración despues del flash Isotermico: "))
print("====================================================== =========================================")
file.write('==================================================== ==========================================='+'\n') file.write(' '+'# comp.'+' '+"Comp. Líquida" + ' '+"Flujo Liquido" + ' '+'Comp. Vapor'+' '+'Flujo Vapor'+'\n')
file.write('==================================================== ===========================================' + '\n') Fv2 = V Fl2 = F - V print("{:^10} {:^20} {:^20} {:^20}".format("Comp Liq", "Flujo Liq", "Comp Vap", "Flujo Vap", )) for w in range(0, n): Zdd = z[w] Kdd = K[w] x = ZddF/(FL+KddV) y = Kdd*x xr = round(x, 7) fl = xr * Fv yr = round(y, 7) fv = yr * Fv suma += x suma1 += y valor = w + 1 print("{:^5} {:^1.4f} {:^10} {:^10,.4f} {:^10} {:^1,.4f} {:^10} {:^10,.4f}".format("X{}=".format(w+1), x, "Fx{}=".format(w+1), fl, "Y{}=".format(w+1), y, "Fy{}=".format(w+1), fv)) file.write(' ') file.write('# % s' % valor) file.write(' ')
file.write('x= % s' % xr) file.write(' ') file.write('Fl=%s' % fl) file.write(' ') file.write('y= % s' % yr) file.write(' ') file.write('Fv=%s' % fv + '\n') print("Iteraciones totales= ", w) file.write(' ') file.write('Iteraciones totales= ') file.write('i= % s' % w + '\n')
print("====================================================== =========================================")
file.write('==================================================== ===========================================' + '\n') Vr = round(V / F, 4) print("{:^8} {:^1}".format("El V/F es=", Vr)) file.write(' ') file.write('El V/F es= % s' % Vr + '\n')
print("====================================================== =========================================")
file.write('==================================================== ===========================================' + '\n') print("{:^8} {:^1,.4f}".format("Sumatoria x=", suma)) sumar = float(round(suma, 4)) sumar1 = float(round(suma1, 4)) file.write(' ') file.write('Sumatoria de x= % s' % sumar + '\n') file.write(' ') print("{:^8} {:^1,.4f}".format("Sumatoria y=", suma1)) file.write('Sumatoria de y= % s' % sumar1 + '\n') print('Por lo tanto esta en equilibrio') file.write(' Por lo tanto esta en equilibrio'+ '\n') FV = V FL = F-V
suma = 0. suma1 = 0. i = 0 x = Symbol('x') while d != 1: d = 0 while i <= (n-1): d += (z[i]K[i]F)/(V(K[i]-1)+F) i += 1 i = 0 fs = 0 while i <= (n-1): de = sp.diff((z[i]K[i]F)/(x(K[i]-1)+F)-1, x) fs += de.subs(x, V) i += 1 V = V - ((d-1)/fs) i = 0 print("{:^10} {:^1.4f}".format("V=", V)) file.write(' ') file.write('V = % s' % V + '\n') FL = F-V print("{:^50}".format("Concentración despues del flash Isotermico: "))
print("====================================================== =========================================")
file.write('==================================================== ==========================================='+'\n') file.write(' '+'# comp.'+' '+"Comp. Líquida" + ' '+"Flujo Liquido" + ' '+'Comp. Vapor'+' '+'Flujo Vapor'+'\n')
file.write('==================================================== ===========================================' + '\n') Fv2 = V Fl2 = F - V print("{:^10} {:^20} {:^20} {:^20}".format("Comp Liq", "Flujo Liq", "Comp Vap", "Flujo Vap", )) for w in range(0, n): Zdd = z[w] Kdd = K[w]
x = ZddF/(FL+KddV) y = Kdd*x xr = round(x, 7) fl = xr * Fv yr = round(y, 7) fv = yr * Fv suma += x suma1 += y valor = w + 1 print("{:^5} {:^1.4f} {:^10} {:^10,.4f} {:^10} {:^1,.4f} {:^10} {:^10,.4f}".format("X{}=".format(w+1), x, "Fx{}=".format(w+1), fl, "Y{}=".format(w+1), y, "Fy{}=".format(w+1), fv)) file.write(' ') file.write('# % s' % valor) file.write(' ') file.write('x= % s' % xr) file.write(' ') file.write('Fl=%s' % fl) file.write(' ') file.write('y= % s' % yr) file.write(' ') file.write('Fv=%s' % fv + '\n') print("Iteraciones totales= ", w) file.write(' ') file.write('Iteraciones totales= ') file.write('i= % s' % w + '\n')
print("====================================================== =========================================")
file.write('==================================================== ===========================================' + '\n') Vr = round(V / F, 4) print("{:^8} {:^1}".format("El V/F es=", Vr)) file.write(' ') file.write('El V/F es= % s' % Vr + '\n')
print("====================================================== =========================================")
Resultados de los balances para todos los componentes partiendo los datos de prueba y considerando F = 1000 mol/hr.
Captura de pantallas que muestren la funcionalidad del programa.
Tablas que muestren el # de iteraciones y otras variables necesarias para observar la convergencia del método.
