Cartilla de Ejercicios

MAY 2015

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Ingenieria de Procesos

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0

Contenido Parte I

INTRODUCCIÓN A LOS SIMULADORES DE PROCESOS Modelos Termodinamicos, propiedades y componentes Administrador básico de la simulación Propiedades de una Mezcla

SIMULACIÓN DE EQUIPOS DE PROCESO Simulacion de Equipos Estacionarios Simulacion de Equipos Rotatorios Siulacion de Tuberias (Pipe Segmet y Gas Pipe)

SIMULACIÓN DE EQUIPOS DE TRANSFERENCIA DE CALOR Air cooler Cooler/Heater Horno (Furnace) Heat Exchanger

USO DE HERRAMIENTAS LOGICAS DE CALCULO Adjust, Set, Balance, Spreadsheet y Recycle.

SIMULACION DE UNIDADES DE PROCESO Simulacion de ciclos de Refrigeracion Simulacion de un ciclo simple de energia de vapor Simulacion de un ciclo de refrigeracion de dos etapas Simulacion de sistemas de Compresion

PRACTICAS PRE-EVALUATIVAS


INTRODUCCIÓN A LOS SIMULADORES DE PROCESOS

En los últimos años, la simulación de procesos ha llegado a ser una herramienta de apoyo para el diseño y evaluación de distintos procesos relacionados a la industria química, hidrocarburos, mineria, etc. La simulación de procesos juega un papel muy importante en las industrias antes mencionadas, como herramienta adecuada y oportuna para el diseño, caracterización, optimización y monitoreo del funcionamiento de procesos industriales.

SOFTWARES DE SIMULACIONES Son herramientas que nos permiten facilitar el diseño, caracterización, optimización y monitoreo del funcionamiento de procesos industriales. Entre los principales tenemos: 

Aspen Hysys



Pro II



Aspen Plus



Chemcad

Figura Nº1: Esquema del paquete AspenOne – Hysys VENTAJAS DEL USO DE SOFTWARE DE SIMULACIÓN 

La simulacion respresenta un sistema real



Es un proceso relativamente eficiente y flexible.

 

Hace posible analizar y sintetizar una compleja y extensa situación real. Permite el diseñador examinar rápidamente varias configuraciones de planta




Permite la experimentación en condiciones que podrían ser peligrosas en el sistema real.



Reduce el tiempo de diseño de una planta.



Desventajas del uso de software de simulación



Un buen modelo de simulación puede resultar bastante costoso; a menudo el proceso a desarrollar un modelo es largo y complicado para su validación.



Por error se producen diferentes resultados en repetidas corridas en el computador.



Cada modelo de simulación es único, las soluciones e inferencias no son usualmente transferibles a otros problemas.



Se requiere gran cantidad de corridas para encontrar “Soluciones óptimas”.



Aplicación de softwares de simulación



Detección de cuellos de botella en la producción.



Predicción de los efectos de cambio en las condiciones de operación de las variables de la planta.



Optimización de las variables de proceso.



Optimización del proceso cuando cambian las características de los insumos y/o las condiciones económicas del mercado.



Evaluación de alternativas de proceso para reducir el consumo de energía.



Análisis de nuevos procesos para nuevos productos



Transformación de un proceso para desarrollar otras materias primas



Análisis de factibilidad y viabilidad de nuevos procesos.



Optimización para minimizar la producción de desechos y contaminantes.



Entrenamiento de operados e ingenieros de procesos.



Investigación de la factibilidad de la automatización de un proceso.

MAPA DE RUTA

Figura Nº2: Ruta para definir una simulación


SIMULACIÓN EN ESTADO ESTACIONARIO Permite evaluar las diferentes alternativas de configuración del proceso y/o determinar la sensibilidad de las variables claves del proceso con respecto a cambios “factibles” en el proceso.

Figura Nº3: Simulación estacionaria – UPD Crudo

SIMULACIÓN EN ESTADO DINÁMICA La simulación Dinámica puede ayudar a mejor el diseño, optimizar y operar procesos. Las plantas nunca están en estado estacionario. Los disturbios de la carga y del medio ambiente, el mal intercambio de calor y la degradación catalítica continuamente afectan las condiciones de un proceso estable. El comportamiento de los sistemas de procesos es mejor estudiado usando una simulación dinámica. Con la simulación dinámica se puede confirmar que la planta puede producir el producto deseado en una forma que sea seguro y fácil de operar. Definiendo detalladamente las especificaciones de los equipos, en la simulación dinámica, se puede verificar que los equipos funcionen como se espera en una situación real de planta también optimizar el diseño de control sin afectar negativamente seguridad de la planta.


Figura Nº4: Simulación dinámica

MODELOS TERMODINAMICOS, COMPONENTES Y PROPIEDADES

Administrador básico de la simulación Ingreso de componentes Aspen Hysys V8.0 contiene una gran variedad de base de datos de componentes químicos, se iniciara el programa Aspen Hysys V8.0, “File”, “New”, “case” y obtendremos la siguiente ventana.

Figura Nº5: Entorno del administrador de propiedades Luego de iniciar el programa en la parte inferior Aspen Hysys V8.0 muestra un panel de mensajes “Message Panel” en el cual indica los errores durante la simulación; para ingresar los componentes para la simulación se ira a “Components Lists”, “Add” y se obtendrá la siguiente ventana


Figura Nº6: Vista de ingreso de componentes

Figura Nº7: Filtro para seleccionar componentes Para seleccionar componentes se utiliza los filtros1 para el tipo de familia o ingresando el nombre formula del componente en el buscador (Nombre de los componentes en idioma ingles)10

Figura Nº8: Vista de ingreso de componentes puros

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Componente hipotético Un componente hipotético puede ser: 

Componentes puros



Mezclas definidas



Mezclas indefinidas



Solidos

Que no se encuentra en la base de datos de Aspen Hysys V8.0; se tiene varios métodos de estimación para asegurar la representación adecuada del comportamiento de os componentes hipotéticos adicionados. También hay métodos para estimar los coeficientes binarios entre componentes hipotéticos. También sirve para clonar los componentes de la lista disponibles en Aspen Hysys y luego modificar sus especificaciones, estos compuestos se pueden usar en reacciones. Los compuestos hipotéticos no pertenecen a un paquete de fluido específico dentro de la simulación porque en su creación son puestos en el Hypo Group. Debido a que estos componentes no pertenecen a un paquete específico estos pueden ser compartidos por varios paquetes en la misma simulación. La información mínima requerida por Aspen Hysys para estimar las propiedades del componente hipotético se resumen en la tabla siguiente: Punto normal de Menor a 700ªF Mayor a 700ªF No se conoce o es un solido

Información mínima Punto normal de ebullición Punto normal de ebullición Densidad del liquido Densidad del liquido Peso molecular

Adicionando un nuevo componente En este caso se trata de estimar propiedades físicas para componente que no está en la base de datos del programa. Para ilustrar el procedimiento estimaremos las propiedades físicas del componente 1,1,2trimetilciclohexano

Figura Nº10: Vista de ingreso de componentes hipotéticos Nota: se debe desactivar el filtro para seguir la búsqueda si no se va a usar más tratar de buscar un componte que no se encuentre dentro del tipo de familia, est ubicada 1


Se construirá la estructura de la molécula en base a los sub-grupos.

SubCH3 CH2 -CCH

Adicionar 3 4 1 1

Figura Nº11: Entorno de componentes hipotéticos

Figura Nº12: Vista de estructura de un componente hipotético


Figura Nº13: Propiedades calculadas - UNIFAC

Figura Nº14: Definición completa de un componente hipotético


Adicionando un nuevo componente especificando propiedades fisicas En este caso se ingresa propiedades físicas para componente deseado que no está en la base de datos del programa. Para ilustrar el procedimiento estimaremos las propiedades físicas de los siguientes componentes.

Figura Nº15: Grados API de alguanas fracciones de Petroleo (Donal Q. Kern)

F IG. 4. Calores específicos de hidrocarburos liquidos. [Holcomb and Brown,] Factor de caracterizacion

Eng.Chem. 34, 505

Figura Nº16: Calor especifico de Hidrocarburos Liquidos

K


Procedimiento Lo mas importante a la hora de ingresar estos compuestos es apoyarse en la definicion de o API, y en la definicion de la correlacion de Gragoe.

�= 𝑀𝑀

6084 44.29 ∙ 𝛾𝛾 = 1.03 − 𝛾𝛾 𝐴𝐴𝐴𝐴𝐴𝐴 − 5.9

Correlacion de Gragoe

1.- Hacer click en la opcion "Hypothetical"

2.- Hacer click en la opcion "Create and Edit Hypos" 3.- Hacer click en la opcion "New Hypo"

4.- Hacer doble click en la opcion "New Hypo2000" y modificar el nombre


PAQUETE DE FLUIDO " FLUID PACKAGE " Clases de Paaquetes de Fluidos Aspen Hysys utiliza el concepto de paquete de fluido o “Fluid Package” como el contenido de toda la información necesaria para desarrollar cálculos de propiedades físicas y evaporaciones espontaneas de corrientes. El paquete de fluido permite definir toda la información (propiedades, componentes hipotéticos, parámetros de interacción, reacciones, datos tabulados, etc.) dentro de un archivo muy sencillo. Son tres las ventajas de esto, a saber:



Toda información asociada se define en una sola localidad, lo que permite la fácil creación y modificación de la información.



Los paquetes fluidos pueden almacenarse como un archivo completo para usarlos en cualquier simulación.



Pueden usarse varios paquetes fluidos n una misma simulación. Sin embargo, todos los paquetes definidos se encuentran dentro del administrador básico de la simulación.

Administrador del paquete básico de la simulación Este paso es muy importante y no se debe tomar al a ligera, ya que definir la base de la simulación. Si tenemos una buena base, tendremos una buena simulación, pero si introducimos un error desde el principio, este se agravara con el desarrollo de la simulación2.

Figura Nº17: Vista de definición de un paquete de fluidos – PR


DETERMINACION DEL PAQUETE DE FLUIDO " FLUID PACKAGE " Clases de Paquetes de Fluidos

Figura Nº18: Paquetes de Fluido Recomendados


ALGORITMO PARA DETERMINACION DEL PAQUETE DE FLUIDO " FLUID PACKAGE " Clases de Paquetes de Fluidos

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Corrientes de mezclas Clases de corrientes en Aspen Hysys Aspen Hysys utiliza el concepto de corrientes de materia y corrientes de energía. Las corrientes de materia requieren, para su completa definición, del a especificación del flujo y de aquellas variables que permitan la estimación de todas sus propiedades físicas y termodinámicas. Las corrientes de energía se utilizan para representar los requerimientos energéticos en unidades como intercambiadores de calor, bombas, etc. Y se especifican, completamente con solo la cantidad de energía intercambia o transferida en dichas unidades. En Aspen Hysys la corriente de materia se observa, por defecto de color azul, mientras que las corrientes de energía es de color rojo.

Corrientes de Energía El elemento más simple que de un diseñador de proceso debe especificar es una simple corriente homogénea. Las variables que definen a una corriente que contiene C componentes son:

Variables Concentraciones Temperatura Presión Flujo Total de Variables

Cantidad C 1 1 1 C+3

Expresando las concentraciones en fracciones molares, Xi, se cumple una restricción de suma entre ellas, es decir que: 𝑁

∑ 𝑋𝑖 = 1 𝑖=1

Por lo tanto, el número de variables de Diseño, 𝑁𝑖𝑒 , que se reuqiern para especificar completamente una corriente de materia es de la diferencia entre el número de variables y el número de restricciones, es decir:

𝑁𝑖𝑒 = 𝐶 + 2 De acuerdo a la ecuación anterior, se define el estado termodinámico de una corriente de materia al conocerse la composición de una corriente de materia y otras dos propiedades, (fracción de vapor, temperatura, presión, entalpía o entropía) una de las cuales debe ser o la temperatura o la presión.

2

Antes de iniciar la simulación se debe considerar la selección de un sistema de unidades


Figura Nº19 Vista de pantalla de simulación

En la versión 8.0 la paleta de objetos se divide en 6 pestañas; Columns, Dynamics, Refining, Common, Custom, Upstream Para adicionar corrientes, hay dos tipos de corrientes, materia y energía (deben especificarse.

A continuación se realizara el ingreso para un corriente de gas y se aplicara todas las aplicaciones que contienen una corriente de materia, utilizando las propiedades de utilities, gráficas y extensiones.

Figura Nº 20: Paleta de objetos


INTRODUCCIÓN DE EQUIPOS DE PROCESO EQUIPOS ESTACIONARIOS Definiendo un separador

1. Agregar un separador a la simulación haciendo clic en el siguiente icono 2. Hacemos doble clic sobre el icono y nos dirigimos a la pestaña “Design” y luego a la pagina “Connections”

separador. 3. En la casilla “Inlets” agregamos la/s corriente/s que van alim a funcionar como 4. En la casilla “Vapour outlet” especificamos la corriente que a vaentar al vapor de salida en el tope del separador

5. En la casilla “Liquid outlet” seleccionamos la corriente que va a funcionar como liquido de salida en el fondo

eparador.

del s agregar una corriente de energía en la casilla 6. De manera opcional podemos “Energy (Optional)”.

Si todas las corrientes anexadas están correctamente definidas, la barra de estato en la parte inferior mostrara un mensaje de Ok tornándose de color verde. En caso de que los cálculos no puedan realizarse, la barra se tornara amarilla indicando que información es necesaria anexar para poder definir completamente la operación. De manera opcional, podemos ingresar a la pagina “Parameters” y agregar una diferencia de presión en la casilla “Delta P”.

Definiendo un separador trifásico

1. Agregar un separador trifásico haciendo clic en el siguiente icono en la paleta de procesos

2.

Hacemos doble clic sobre el icono y nos dirigimos a la pestaña “Design” y luego a la pagina “Connections”

3. En la casilla “Inlets” agregamos la/s corriente/s que van a alimentar al separador. 4. En la casilla “Vapour” especificamos la corriente que va a funcionar como vapor de salida en el tope del separador 5. En la casilla “Light liquid” seleccionamos la corriente que va a funcionar como liquido liviano a la salida del separador. 6. En la casilla “Heavy liquid” seleccionamos la corriente que va a funcionar como liquido pesado a la salida del separador. Si todas las corrientes anexadas están correctamente definidas, la barra de estado en la parte inferior mostrara un mensaje de Ok tornándose de color verde. En caso de que los cálculos no puedan realizarse, la barra se tornara amarilla indicando que información es necesaria anexar para poder definir completamente la operación. De manera opcional, podemos ingresar a la pagina “Parameters” y agregar una diferencia de presión en la casilla “Delta P”.


EQUIPOS ROTATORIOS Definiendo una bomba.

1. Agregar una bomba haciendo clic en el siguiente icono en la paleta 2. Hacemos doble clic sobre el icono y nos dirigimos a la pestaña “Design” y luego a la pagina “Connections”

3. En la casilla “Inlet” especificamos la corriente de entrada. 4. En la casilla “Outlet” especificamos la corriente de salida. 5. En la casilla “Energy” especificamos una corriente de energía para el duty del equipo. De manera opcional podemos especificar los valores para la caída de presión y el duty en la ventana “Parameters”, pero no es totalmente necesario ya que el valor del duty se calcula automáticamente al especificar una corriente de energía y el valor de la caída de presión se calcula automáticamente si se ha especificados la presión en las corrientes de entrada y salida. Si todos los valores han sido especificados correctamente, la barra de estado en la parte inferior mostrara un mensaje de Ok tornándose de color verde. En caso de que los cálculos no puedan realizarse, la barra se tornara amarilla indicando que información es necesaria anexar para poder definir completamente la operación.

Definiendo un compresor/expansor

1. Agregamos un compresor o un expansor desde la paleta de equipos. 2. Hacemos doble clic sobre el icono y nos dirigimos a la pestaña “Design” y luego a la pagina “Connections”

3. 4. 5. 6. 7.

En la casilla “Inlet” especificamos la corriente de entrada. En la casilla “Outlet” especificamos la corriente de salida. En la casilla “Energy” especificamos una corriente de energía. Hacemos clic en la pagina “Parameters” Especificamos un valor para el duty dentro de la casilla “Duty”

Si todos los valores han sido especificados correctamente, la barra de estado en la parte inferior mostrara un mensaje de Ok tornándose de color verde. En caso de que los cálculos no puedan realizarse, la barra se tornara amarilla indicando que información es necesaria anexar para poder definir completamente la operación.

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 TUBERIAS (PIPE SEGMENT Y GAS PIPE) Definiendo una sección de tubería de gas

1. Agregamos una sección de tubería de gas haciendo clic sobre el siguiente icono 2. Hacemos doble clic sobre el icono y nos dirigimos a la pestaña “Design” y luego a la pagina “Connections”

3. 4. 5. 6. 7. 8. 9. 10.

En la casilla “Inlet” especificamos la corriente de entrada. En la casilla “Outlet” especificamos la corriente de salida. Hacemos clic en la pagina “Parameters” En la casilla “Pressure Drop” especificamos un valor para la caída de presión. Hacemos clic en la pestaña “Rating” en la pagina “Sizing”. Si deseamos agregar múltiples secciones de tubería hacemos clic en el botón “Add” Especificamos todos los valores necesarios para el dimensionamiento. Hacemos clic en la pagina “Heat transfer” y especificamos valores para la temperatura ambiental y el coeficiente global de transferencia de calor

Si todos los valores han sido especificados correctamente, la barra de estado en la parte inferior mostrara un mensaje de Ok tornándose de color verde. En caso de que los cálculos no puedan realizarse, la barra se tornara amarilla indicando que información es necesaria anexar para poder definir completamente la operación.

Definiendo un segmento de tubería

1. Agregamos un segmento de tubería haciendo clic en el siguiente icono 2. Hacemos doble clic sobre el icono y nos dirigimos a la pestaña “Design” y luego a la pagina “Connections” 3. En la casilla “Inlet” especificamos la corriente de entrada. 4. En la casilla “Outlet” especificamos la corriente de salida. 5. En la casilla “Energy” especificamos una corriente de energía. 6. Hacemos clic en la pagina “Parameters” 7. En la casilla “Pipe Flow Correlation” seleccionamos la correlación que deseamos usar para la tubería. 8. Especificamos un valor para la caída de presión en la casilla “Delta P” (el valor es automáticamente calculado si las corrientes de entrada y salida ya poseen especificaciones de presión). 9. Hacemos clic en la pestaña “Rating” en la pagina “Sizing”.

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Tuberías (Pipe Segment) El segmento de tubo se utiliza para simular una amplia variedad de tuberías situaciones que van desde tuberías único o de múltiples fases planta con estimación de transferencia de calor riguroso, para una gran capacidad de bucle problemas de tuberías. Ofrece varias correlaciones de caída de presión:

Modos de cálculo La operación segmento de tubería contiene cuatro modos de cálculo:



Caída de presión



Longitud



Flujo



Diámetro

El modo se asigna de forma automática en función de lo se especifica información. Independientemente del modo que utilice, debe especificar el número de incrementos en el tubo. Los cálculos se realizan en cada incremento, por ejemplo, para determinar la caída de presión, los cálculos de los balances de energía y de masa se realizan en cada incremento, y la presión de salida.

Incremental Balances de materia y energía El algoritmo general consta de tres bucles anidados. El exterior estructura se repite en los incrementos (de presión, longitud o de flujo Modo), el bucle medio resuelve para la temperatura, y el bucle interior resuelve por presión. Los bucles medio e interno implementar un método de la secante para acelerar la convergencia.

La presión y la temperatura se calculan como sigue: 𝑄 = 𝑈 ∗ 𝐴 ∗ ∆𝑇𝐿𝑀 𝑄 = 𝑄𝑖𝑛 − 𝑄𝑜𝑢𝑡 Donde: 

Q

: cantidad de calor transferido



U

: coeficiente global de transferencia de calor



A

: área de transferencia de calor exterior



∆𝑇𝐿𝑀

: Log media de diferencia de T



𝑄𝑖𝑛

: El flujo de calor de la corriente de entrada



𝑄𝑜𝑢𝑡

: El flujo de calor de la corriente de salida

Figura Nº 21: Conexión de corrientes

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Resumen de métodos Los métodos anteriores han sido desarrollados para la predicción de dos fases caídas de presión. Algunos métodos se desarrollaron exclusivamente para el flujo en tuberías horizontales, otros exclusivamente para flujo en tuberías verticales, mientras que algunos se pueden usar para cualquiera. Algunos de los métodos definen un mapa régimen de flujo y pueden aplicar específico correlaciones de caída de presión de acuerdo con el tipo de flujo predicho. Algunos de los métodos de calcular el líquido esperado atraco en el flujo de dos fases mientras que otros asumen una homogénea mezcla. La siguiente tabla resume las características de cada modelo. Información más detallada sobre cada modelo se presenta más adelante en esta sección.

Figura Nº22: Modelos para calculos de mecanica de Fluidos

No Utilizar con cuidado Si No

Flujo vertical Si Si Si Si

Almacenamiento de flujo Si No Si Si

Mapa de flujo Si No Si Si

Si

No

Si

Si

No Si Si Si No No

Si Si Si Si Si Si

Si No Si Si Si No

No No No Si Si No

Si

Si

Si

Si

No

Si

Si

Si

Modelo

Flujo Horizontal

Aziz, Govier & Fogarasi Baxendell & Thomas Begg & Brill Duns & Ros Greogory, Aziz, Mandhane Hagedorn & Brown HTFS Homogeneous HTFS Liquid slip Olgas 2000 Orkisewki Poettman & Carpenter Tacite Hydrodynamic Module Tulsa

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Seleccione una de las siguientes



Actual. el diámetro nominal no se puede especificar



Cédula 40



Cedula 80



Cedula 160

Tipo de material de la tubería Drawn Tube Dibujar tubo Mild Steel Acero dulce Asphalted Iron Acero asfaltado Galvanized Iron Acero galvanizado Cast Iron Hierro fundido Smooth Concrete Hormigón liso Rough Concrete Hormigón rugoso Smooth Steel Acero liso Rough Steel Acero rugoso Smooth Wood Stave Madera lisa Rough Wood Stave Madera rugosa

Rugosidad absoluta, m 0.0000015 0.0000457 0.0001220 0.0001520 0.0002590 0.0003050 0.0030500 0.0009140 0.0091400 0.0001830 0.0009140

Heat loos Si se conoce la exigencia de calor total de la tubería, el balance de energía puede calcularse inmediatamente. Cada incremento se supone tener la misma pérdida de calor. Usted entra a la pérdida de calor de la tubería en el campo de la pérdida de calor. Esta suposición es válida cuando el perfil de temperatura es plana, lo que indica las bajas tasas de transferencia de calor en comparación con los flujos de calor de las corrientes. Este es el más rápido método de solución. Si se especifican tanto las temperaturas de entrada y salida, lineal perfil se asume y HYSYS puede calcular el calor total deber. Este método permite el cálculo rápido cuando la corriente se conoce las condiciones. Seleccione el botón de Pérdida de Calor para ver el calculado deber general de calor.

Overall HTC Si el HTC total (coeficiente de transferencia térmica) y una temperatura ambiente representativa son conocidos, el calor riguroso cálculos de transferencia se llevan a cabo en cada incremento.

Figura Nº22: Perdida de calor


Segment HTC7 Si el coeficiente de transferencia de calor y un ambiente representativo temperatura son conocidos para cada segmento. Se puede especificar la temperatura ambiente y HTC para cada segmento de tubería que ha sido creado en la página acerca. HYSYS realiza calor riguroso transferir cálculos en cada incremento.

Cuando se selecciona el botón de radio Estimación HTC, el Heat Transferir los cambios de página a la ventana de propiedades se muestra en la figura a continuación. Métodos se utilizarán únicamente para tuberías monofásicas que operan a números de Reynolds alto (> 10.000). Los métodos Profes y HTFS deben proporcionar mucho mejor resultados para dos y tres sistemas de fase, y en el flujo laminar región a costa de cierto incremento en el tiempo de cálculo. En general la opción Profes se recomienda para la mayoría de tuberías aplicaciones, ya que tiene en cuenta plenamente el régimen de flujo en la tubería y es razonablemente eficiente en el cálculo. La opción es HTFS más cálculo intensivo, sobre todo en dos aplicaciones en fase donde se requieren cálculos de flash adicionales. Es se recomienda su uso en casos con un alto flujo de calor con alta temperaturas delta entre el contenido de tuberías y ambiente condiciones. Las cinco correlaciones proporcionadas son:



Petukov (1970)

ℎ= 

(𝑓⁄8)𝑅𝑒𝑑 𝑃𝑟 𝑘 ∗ 𝑑 1.07 + 12.7(𝑓⁄8)1⁄2 (𝑃𝑟 2⁄3 − 1)

Dittus and Boelter (1930) 𝑘 ∗ 0.023𝑅𝑒𝑑0.8 𝑃𝑟 𝑛 𝑑 Donde: ℎ=

0.4→𝑝𝑎𝑟𝑎 𝑐𝑎𝑙𝑒𝑛𝑡𝑎𝑚𝑖𝑒𝑛𝑡𝑜

𝑛 =0.3→𝑝𝑎𝑟𝑎 𝑒𝑛𝑓𝑟𝑖𝑎𝑚𝑖𝑒𝑛𝑡𝑜 

Sieder and Tate (1936) Para flujo de 2 fases

7 El HTC general y Estimación HTC se pueden utilizar juntos para definir la información de transferencia de calor para el tubo. Si sólo conoce la temperatura ambiente, puede suministrarla en la sección general de HTC y tienen el valor global de HTC calculado por la sección Estimación HTC. Del mismo modo, es necesario para especificar la temperatura ambiente en la estimación de HTC sección para el segmento de tubería para tener la transferencia de calor suficiente información para resolver.


Para flujo de una sola fase



Profes Implementa los métodos utilizados por el programa Profes Pipe Simulación (antes PLAC). Los métodos se basan en los mapas de flujo Profes para flujo horizontal y vertical, y las correlaciones apropiadas se utilizan para determinar el HTC en cada región del mapa de flujo.



HTFS. Implementa los métodos utilizados por los programas HTFS. Correlaciones independientes se utilizan para la ebullición y condensación de transferencia de calor, y para el flujo horizontal y vertical. Los métodos utilizados están documentadas en el manual HTFS.

Usted puede optar por incluir la resistencia térmica de la tubería en su Cálculos HTC seleccionando la casilla de verificación Incluir pared del tubo. La activación de esta opción requiere que se defina la conductividad térmica para el material de la tubería en la ventana de propiedades detalle de cada segmento de tubería. Los valores por defecto de conductividad térmica se proporcionan para los materiales estándar que se pueden seleccionar en el segmento de tubería.

Figura Nº23: Segmento HTC

Estimate HTC

Figura Nº24: Calculo de transferencia de calor


Conducción Fuera / convección Fuera de convección ya sea aire, agua o tierra pueden ser incluidos mediante la selección de la casilla de verificación Incluir HTC exterior. Para el aire y agua, la velocidad del medio ambiente está por defecto en 1 m/s, y es modificable por el usuario. La transferencia de calor por convección fuera coeficiente de correlación es para el flujo pasado tubos horizontales (JP Holman, 1989):

Si se selecciona el suelo como el medio ambiente, el tipo de suelo, entonces se puede seleccionar. La conductividad térmica de este medio aparece, pero también es modificable escribiendo encima del valor predeterminado. Los tipos de tierra y sus correspondientes conductividades se tabulan a continuación:

Tipo de Tierra

Conductivid ad (W/mK)

Tipo de Tierra

Conductivida d (W/mK)

Dry Peat

Turba seca

0.17

Arcilla congelada

2.50

Wet Peat

Turba húmeda

Frozen Clay

0.54

Gravel

Grava

1.10

Icy Peat

Turba helada

1.89

Sandy Gravel

Grava arenosa

2.50

Arena seca

0.50

Limestone

Caliza

1.30

0.95

Sandy Stone

Piedra arenosa

1.95

2.20

Ice

Hielo

2.20

0.48

Cold Ice Loose Snow

Helada

2.66

Nieve suelta

0.15

Hard Snow

Nieve dura

0.80

Dry Sand Moist Sand Wet Sand Dry Clay Moist Clay Wet Clay

Arena húmeda Arena mojada Arcilla seca Arcilla húmeda Arcilla mojada

0.75 1.40

En HYSYS, el valor del coeficiente de transferencia de calor de los alrededores es basado en la siguiente ecuación resistencia a la transferencia de calor:

Donde: 

Hentorno : coeficiente de calor que rodea



Rentorno : entorno de resistencia a la transferencia de calor



Zb

: la profundidad de la cobertura a la línea central de la tubería.



Ks

: conductividad térmica tubería – material circundante (aire, agua, suelo)



Dot

: diámetro exterior de la tubería, incluyendo el aislamiento


Conducción a través de aislamiento Conducción a través del aislamiento o cualquier otro revestimiento de la tubería también se puede especificar. Varios materiales representativos se proporcionan, con sus respectivas conductividades térmicas. Debe especificar un grosor para este recubrimiento.

Profes Método Ceras La deposición de la cera del aceite a granel sobre la pared de la tubería es supone que sólo ser debido a la transferencia de masa, la dispersión de cizallamiento no es considera que es un factor significativo. La tasa de deposición es descrito por:

Aislamiento/tubería Evacuated Annulus Urethane Foam

Anillo evacuado Espuma de uretano

Conductividad (W/mK)

Asphalt

Asfalto

0.700

0.018

Concrete

Hormigón

1.00

0.080

Concrete Insulated

Hormigón con aislamiento

0.500

0.035

Neoprene

Neopreno

0.250

0.070

PVC Foam

Espuma de PVC

0.040

Manta de fibre

0.030

PVC block

Bloque de PVC

0.150

Bloque de plastico

0.036

PolyStyrene Foam

Espuma de poliestireno

0.027

Bloque de vidrio

Fiberglass Block

Bloque de fibra de vidrio Manta de fibra

Plastic Block

Aislamiento/Tubería

0.005

Glass Block

Fiber Blanket Fiber Blanket – Vap Barr

Conductividad (W/mK)


PRACTICANDO FAMILIARIZACION CON EQUIPOS BASICOS DE PROCESO El objetivo de la siguiente practica es que el estudiante, pueda familializarse con los equipos basico de proceso, con el fin de entender los grados de libertad que gobierna cada uno. Introducir los siguientes datos para dar inicio a la simulacion

DIAGRAMA DE FLUJO


PRACTICA DE CLASE 1 De un tanque de almacenamiento cerrado provisto de un respiradero a la atmosfera se desea bombear agua a 20 C, hacia una torre de absorcion. El nivel de liquido en el tanque se encuentra a 1.5 metros sobre el eje de la bomba, el caudal es de 20 m3/h. La conexion de entrada del agua en el tope de la torre se halla a 7.5 m sobre el nivel del eje de la bomba. La linea de succion consiste de tuberia de acero dulce de 2" (5.08 cm) de diametro nominal, No. de cedula 40 y 10 m de longitud, posee 2 codos estandar y una valvula de compuerta " gate" abierta. La linea de descarga tambien es de 2" No, de cedula 40 y 60 m de longitud, tiene 2 codos estandar, y una valvula de control, y la presion manometrica minima requerida es de 10 psig. Calcular la potencia necesaria de la bomba para realizar un presion de descarga de 50 psia, asumiendo que a lo largo del tramo II habra grandes perdidas de presion, por un mal disheno de caneria.

Paso 1.- Ingresar las condiciones y compocision del fluido


Paso 2.- Insertamos un segmento de Tuberia, especificando un flujo de salida y su corriente de energia, e ingresamos a la ventana del pipe segment haciendo doble "click".

Paso 3.- Ahora nos vamos a la opcion " Parameters" y damos en la opcion " Beggs and Brill (1979)" correlacion que se adecua a la mayoria de patrones de flujo

Paso 4.- Ahora nos vamos pestaña "Rating" y damos en la opcion " Append Segement" lo cual es agregar segemento de tuberia


Paso 5.- Especificamos los datos de tuberia correspondiente de la siguiente manera: 1.- Hacemos click en la opcion "View Segment", el cual nos llevara a una ventana de informacion de la tuberia "Pipe Info" 2.-Escogemos el esquema planteado por la pregunta, en este caso "Schedule 40" 3.-Ahora escpgemos el diametro nominal de 2" (50.80 mm). 4.Proximo a ello damos en la opcion especificar "Specify"

Paso 6.- Escogemos la tuberia Acero Dulce " Mild Steel" y agregamos, los accesorios descritos por el problema, comjuntamente con la longitud de la tuberia y el cambio de elevacion "Elevation Change"

Nota: Por lo tanto al agregar todos los paramteros, notaremos que se tornara de color verde por lo que se ha ejecutado correctamente.


Paso 7.- De la paleta de objetos, agregamos la bomba, ahora agregamos el grado de libertad requerido por la bomba, en este caso se agregara 50 psia a la corriente de salida de la bomba.

Ahora vemos que habra corrido correctamente, tornadose de color verde.

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Paso 8.- Ahora agregamos el segundo segmento de tuberia, que vendria a ser especificado de la siguiente manera.

Paso 9.- Seguidamente nos vamos a la opcion "Heat Transfer" y digitamos 0.00 en perdidas de calor "Heat Loss", y observamos que se tornara de color verde.

Por lo tanto nuestro diagrama de flujo quedara de la siguiente manera: DIAGRAMA DE FLUJO DE PROCESO


PRACTICA DE CLASE 2 Se desea transportar en un gasoducto Gas Nitrogeno a una estacion de compresion, la condiciones aguas arriba son 25 C y 206.2 Kpa con un flujo masico de 100 kg/h, el cambio de elevacion que presenta el terreno hasta la estacion es de 15m. La tuberia esta especificada como SCH 40, de material " Mild Steel" de 2" de diametro nominal. Calcule las condiciones a la que llegas a la estacion de compresion y se debe verificar que la tuberia opera en buenas condiciones (Verificar la velocidad).

GASODUCTO

CONSIDERACIONES • Considere que no hay transferencia de calor al ambiente

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 PRACTICA DE CLASE 3 Se desea descargar agua de un Tanque de almacenamiento de 100 m3, contenido en ella un 85% agua de su volumen total, a una altura de 1 m sobre el eje de la bomba, el cual se desea descargar a una unidad de proceso de requerimiento 30 psia. La distancia del tanque a la bomba es de 100 m de longitud, utilizando una cañería de acero al carbono sin revestimiento con las siguientes características.

Características de la tubería de Transporte Material Longitud Conductividad tubería Rugosidad

Acero al Carbono [m] [W/m-K] [mm]

100 38.66 0,04572

CALCULE:

a) El diámetro nominal de la cañería para un esquema 40 STD, tomado como mínimo 4psi/100 ft. (parámetro recomendado por ingeniería). b) Cuál debe ser la potencia de la bomba, para descargar el agua en 5 horas. c) Cuál sería la potencia de la bomba haciendo un diseño “STAND BY”, utilizando 2 bombas en paralelo. CONSIDERACIONES • • • • •

Considere que no hay transferencia de calor al ambiente Considere los “fittings” necesarios para el diseño. Considere que el nivel del tanque se mantiene constante Tome como referencia la boquilla de salida de flujo del tanque Asumir que el nivel del tanque esta 1 metro por encima del nivel de entrada de sucion de la bomba

DIAGRAMA DE FLUJO DEL PROCESO

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 PRACTICA N0 1 Se Almacena Gas Licuado del Petróleo en un Tanque esférico de 6 m de diámetro externo, el volumen contenido de GLP es el 85% de su volumen total, un Ingeniero de diseño debe hallar la potencia necesaria para descargar el GLP en 2 horas. En mantenimiento se dispone de una bomba cuya eficiencia es del 45%, por malas rutinas de operación, el jefe de la unidad indica que debe implementarse dicha bomba, ya que el area de bombeo cuenta con una sola bomba de eficiencia del 75% de tipo centrifuga y potencia conocida. Los tramos de cañería esquema 80 están especificado de la siguiente manera:

LAS CONDICIONES DE OPERACION Y COMPOCISION DEL FLUIDO SON :

CONSIDERACIONES • • • •

Considere una caida de presion de 4psi/100 ft en los tramos de cañería I y II Considere los “fittings” necesarios para el diseño. (Si Corresponde) Considere que el nivel del tanque se mantiene constante Asuma que hay un cambio de elevacion en los tramos I de 1.5 m" Nivel del liquido sobre el eje de la bomba"

•

Asumir que la reduccion en la TEE a la entrada de area del bombeo es de 3"/2".


NOTA Las Bombas tienen las mismas caracteristicas lo que implica que el conjunto motor-bomba es el mismo para ambos. Use todas las operaciones logicas necesarias para la resolucion del problema.

DIAGRAMA DE FLUJO DE PROCESO

CALCULE: a) Cual es el diametro nominal de Tuberia usando un SCH 80 b) Calcule la potencia equivalente de las bombas P-101 y P-102, para tener una presion de descarga de 20 Bar al final del Tramo II c) Cual es la Presion de descarga de las bombas d) Cuál es la relacion de flujo (R) optimo para satisfacer las condiciones de operacion de la bombas e) Cual es la velocidad adquirida por el fluido en la TEE reductora, si se toma en cuenta como parametro limite de diseño 3 m/s, velocidad tipica en tuberias de proceso.


SIMULACION DE EQUIPOS DE TRANSFERENCIA DE CALOR


INTRODUCCION La Transferencia de calor es una parte esencial de la mayoria de los procesos quimicos. El tipo de equipo de Transferencia de calor mas comunmente usado es el omnipresente intercambiador de calor de tubo y coraza. La palabra "Intercambiador" realmente se aplica a todo tipo de equipos en los que se intercambia calor pero suele emplearse mas especificamente para denominar equipos en los que se intercambia calor pero suele emplearse mas especificamente para denominar equipos en los se intercambian calor entre dos corrientes de un proceso. Los intercambiadores en los que se calienta o se enfria un fluido de proceso mediante corrientes procedentes una planta de servicios se denominan calentadores y refrigerantes. Si el intercambiador se evapora la corriente del proceso, se le denomina evaporizador, caldera si esta instalado en una columna de destilacion, y evaporador si se usa para concentrar una disolucion.



EDGAR MARTIN JAMANCA ANTONIO

40



42



43



45



46



47



48



49



50

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Separación de fases Instantáneas (Separator & Tank) El separador puede tener múltiples entradas, hay dos boquillas de producto  Vapor  Liquido Cuando está en régimen estacionario se define a continuación el balance de energía: 𝐻𝑓𝑒𝑒𝑑 ± 𝐷𝑢𝑡𝑡𝑦 = 𝐻𝑉𝑎𝑝𝑜𝑟 + 𝐻ℎ𝑒𝑎𝑣𝑦 + 𝐻𝑙𝑖𝑔ℎ𝑡 Donde    

Hfeed Hvapor Hlight Hheavy

: el flujo de calor de la corriente de alimentación. : el flujo de calor de la corriente de producto vapor : el flujo de calor de la corriente de producto liquido ligero : el flujo de calor de la corriente de producto liquido pesado

La caída de presión a través del recipiente se define como: 𝑃 = 𝑃𝑙 = 𝑃𝑓𝑒𝑒𝑑 − ∆𝑃 = 𝑃ℎ𝑒𝑎𝑑 + 𝑃𝑣 Donde:      

P Pv Pl Pfeed ΔP Phead

: presión del separador : presión de la corriente de producto de vapor : presión de la corriente liquida : Presión de la corriente de alimentación. : caída de presión en el separador : presión de la carga estática

El volumen del separador, junto con el punto de ajuste para el nivel/flujo de líquido, define la cantidad de atraco, en el recipiente en cualquier momento viene dado por la siguiente expresión. 𝐴𝑡𝑟𝑎𝑐𝑜 = 𝑣𝑜𝑙𝑢𝑚𝑒𝑛 𝑑𝑒𝑙 𝑠𝑒𝑝𝑎𝑟𝑎𝑑𝑜𝑟 𝑥

𝑃𝑉(%𝐹𝑢𝑙𝑙) 100

Donde  PV(%full)

: nivel de líquido en el recipiente en el tiempo t


51

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Figura Nº63: Ventana de conexiones de un separador y tanque

Figura Nº64: Ventana de WorkBook de un separador y un tanque Sizing Se utiliza esta ventana para caracterizar la orientación del separador o tanque

Figura Nº65: Sizing, geometría del equipo Weir (Vertedero) Un vertedero puede especificarse para el separador de cilindro plano horizontal al seleccionar Enable Weir8

La casilla de verificación permitir vertedero solo está disponible para la opción de forma de cilindro plano 8


52


Figura Nº66: Instalando las posiciones del vertedero La vista inicial atraco propiedad permite especificar la altura del vertedero y la posición. La posición del vertedero es la distancia de la presa es desde el lado de alimentación del vaso. Cuando Aspen Hysys simula, el vertedero tiene dos volúmenes en el interior del separador, llamados ángulo de caída 1 y la cámara 2, pero todavía hay sólo un mayor volumen atraco y moles hasta el solucionador de flujo de presión se refiere. Esto significa que las composiciones y propiedades de las fases en los dos volúmenes son los mismos

Figura Nº67: vertedero y los ángulos de caída Boquillas (Nozzles)

Figura Nº68: Dimensión del separador, diámetro de las boquillas


53

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Pérdida de calor (Heat Loss) Modelo simple 𝑄 = 𝑈𝐴(𝑇𝑓 − 𝑇𝑎𝑚𝑏 )

Figura Nº69: Perdida de calor, modelo simple Modelo detallado El modelo detallado le permite especificar los parámetros más detallados de transferencia de calor

Figura Nº70: Perdida de calor, modelo detalloso

Grifos9 de nivel (Level Taps) Ya que el contenido en un recipiente se puede, distribuidos en diferentes fases, la página de los grifos de nivel le permite controlar el nivel de Líquido y acuosa contenidos que coexisten en la zona especificada en un tanque o un separador

9

La información disponible en la página solo de da en casos dinámicos EDGAR MARTIN JAMANCA ANTONIO

54


 PV  OP

Figura Nº71: Especificación de grifos de nivel : límite r de la sección a ser monitoreados. Se expresa en metros : límite de la salida de la escala de normalización

Carry Over Model

Figura Nº72: Carry Over, Feed Basis


55


Figura Nº73: Resultados de Cover over model

Ejemplo N° Separar la siguiente corriente en una fase liquida y otra vapor


56

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Figura Nº74: Composición de la corriente A una presión de 658.8 kPa y una temperatura de 25°C utilizando un flujo molar de 100 kgmole/h

Figura Nº75: Resultados de la separación

SIMULACIÓN DE UNIDADES DE PROCESO (Transferencia de calor) Air cooler El funcionamiento de la unidad del refrigerador de aire utiliza una mezcla ideal de aire como el calor transferir medio para enfriar (o calor) una corriente de proceso de entrada a un necesaria condición corriente de salida. Uno o más ventiladores circulan el aire a través de haces de tubos para enfriar fluidos de proceso. El flujo de aire Se pueden especificar o derivado de los datos de características del ventilador. El refrigerador de aire puede resolver muchos conjuntos diferentes de especificaciones incluyendo la: 

En general, el coeficiente de transferencia de calor, UA EDGAR MARTIN JAMANCA ANTONIO

57

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 

El flujo de aire total



Temperatura de la corriente de salida

Estado Estacionario Utiliza la misma ecuación básica y el Intercambiador de calor, funcionamiento de la unidad, sin embargo, la operación del refrigerador de aire puede calcular el flujo de aire sobre la base de la información de clasificación del ventilador. Los cálculos de aire más frío se basan en un balance de energía entre las corrientes de aire y de proceso. Para un aire contracorriente Enfriador, el balance de energía se calcula como sigue: 𝑀𝑎𝑖𝑟𝑒 (𝐻𝑜𝑢𝑡 − 𝐻𝑖𝑛 )𝑎𝑖𝑟𝑒 = 𝑀𝑝𝑟𝑜𝑐𝑒𝑠𝑜 (𝐻𝑖𝑛 − 𝐻𝑜𝑢𝑡 )𝑝𝑟𝑜𝑐𝑒𝑠𝑜 Donde:  

𝑀𝑎𝑖𝑟𝑒 : caudal másico de la corriente de aire 𝑀𝑝𝑟𝑜𝑐𝑒𝑠𝑜 :caudal másico de la corriente de proceso



H

:entalpía

El deber del refrigerador de aire, Q, se define en términos del calor global coeficiente de transferencia, el área disponible para intercambio de calor, y el registro de diferencia de temperaturas media: 𝑄 = −𝑈𝐴∆𝑇𝐿𝑀 𝐹𝑡 Donde 

U




A

: área de superficie disponible para la transferencia de calor



∆𝑇𝐿𝑀

: iniciar diferencia media de temperatura



𝐹𝑡

: factor de corrección

El factor de corrección LMTD, Ft., se calcula a partir de la geometría y la configuración del refrigerador de aire.

Transferencia de Calor El refrigerador de aire utiliza las mismas ecuaciones básicas del balance de energía como la operación de la unidad de intercambiador de calor. Los cálculos de aire más frío se basan en un balance de energía entre el proceso de aire y arroyos. Para una corriente transversal del refrigerador de aire, se muestra el balance energético como de la siguiente manera: 𝑀𝑝𝑟𝑜𝑐𝑒𝑠𝑜 (𝐻𝑖𝑛 − 𝐻𝑜𝑢𝑡 )𝑝𝑟𝑜𝑐𝑒𝑠𝑜 − 𝑀𝑎𝑖𝑟𝑒 (𝐻𝑜𝑢𝑡 − 𝐻𝑖𝑛 )𝑎𝑖𝑟𝑒 = 𝜌

𝑑(𝑉𝐻𝑜𝑢𝑡 )𝑝𝑟𝑜𝑐𝑒𝑠𝑜 𝑑𝑡

Donde: 

𝑀𝑎𝑖𝑟𝑒 : caudal másico de la corriente de aire EDGAR MARTIN JAMANCA ANTONIO

58


𝑀𝑝𝑟𝑜𝑐𝑒𝑠𝑜



𝜌

: Densidad



𝐻

: Entalpia



𝑉

: volumen de tubo enfriador de aire

:caudal másico de la corriente de proceso

Caída de presión La caída de presión del refrigerante del aire se puede determinar en una de dos maneras: • Especificar la caída de presión. • Definir una relación de flujo de presión en el enfriador de aire por especificando un valor k. Si se elige la opción de flujo de presión para la caída de presión determinación en el refrigerador de aire, un valor k se utiliza para relacionar la pérdida de presión por fricción y el flujo a través del intercambiador. Este relación es similar a la ecuación de la válvula en general: 𝐹𝑙𝑢𝑗𝑜 = √𝑑𝑒𝑛𝑠𝑖𝑑𝑎𝑑 ∗ 𝑘√𝑃1 − 𝑃2 La ecuación de flujo general utiliza la caída de presión a través de la Intercambiador de calor y sin ninguna contribución de cabeza estática. La cantidad, P1 - P2, se define como la pérdida de presión por fricción que se utiliza para el "tamaño" del enfriador de aire con un valor k.

Ejemplo N° 7 Condiciones de operación, 100 bar_g, Temperatura 191°F y un flujo molar de 100 kgmole/h, se requiere enfriar la corriente en un refrigerador de aire, se quiere enfriar la corriente a una temperatura de 110.5°F, se desprecia la caída de presión, solo utilizar un ventilador, determinar la temperatura de la salida del aire

Figura Nº75: Composición molar de la corriente


59


Figura Nº75: Parámetros de Operación

Figura Nº75: Resultados del gas enfriado

Cooler/Heater10

Las operaciones del refrigerador y el calentador son el calor de un solo lado intercambiador. La corriente de entrada se enfría (o calentado) a la condiciones de salida requeridos, y los absorbe corriente de energía (o proporciona) la diferencia de entalpía entre las dos corrientes.

10

La diferencia entre el enfriador y el calentador es la convención de signos balance energético.


60

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Estas operaciones son útiles cuando usted está interesado sólo en cómo Se requiere mucha energía para enfriar o calentar una corriente de proceso con una utilidad, pero usted no está interesado en las condiciones de la utilidad en sí.

Estado Estacionario La diferencia principal entre un enfriador y un calentador es el signo convención. Se especifica el flujo de la energía absoluta de la utilidad corriente, y luego se aplica HYSYS que el valor de la siguiente manera:

Para una, la entalpía o calor de flujo del refrigerador de la energía corriente se resta de la de la corriente de entrada: 𝐹𝑙𝑢𝑗𝑜 𝑑𝑒 𝐶𝑎𝑙𝑜𝑟𝑖𝑛𝑔𝑟𝑒𝑠𝑜 − 𝐷𝑢𝑡𝑦𝑐𝑜𝑜𝑙𝑒𝑟 = 𝐹𝑙𝑢𝑗𝑜 𝑑𝑒 𝑐𝑎𝑙𝑜𝑟𝑆𝑎𝑙𝑖𝑑𝑎

Para un calentador, se añade el flujo de calor de la corriente de energía: 𝐹𝑙𝑢𝑗𝑜 𝑑𝑒 𝐶𝑎𝑙𝑜𝑟𝑖𝑛𝑔𝑟𝑒𝑠𝑜 + 𝐷𝑢𝑡𝑦ℎ𝑒𝑎𝑡𝑒𝑟 = 𝐹𝑙𝑢𝑗𝑜 𝑑𝑒 𝑐𝑎𝑙𝑜𝑟𝑆𝑎𝑙𝑖𝑑𝑎 Caída de presión La caída de presión del enfriador / calentador se puede determinar en una de dos maneras: 

Especificar la pérdida de carga de forma manual.



Definir una relación de flujo de presión en el enfriador o calentador por especificando un valor k.

Si se elige la opción de flujo de presión para la caída de presión determinación en el enfriador o calentador, un valor k se utiliza para relacionar la pérdida de presión por fricción y el flujo a través del enfriador / calentador.

La relación es similar a la ecuación de la válvula en general: 𝐹𝑙𝑢𝑗𝑜 = √𝑑𝑒𝑛𝑠𝑖𝑑𝑎𝑑 ∗ 𝑘√𝑃1 − 𝑃2

Pérdida de calor Página Valoración información con respecto a la pérdida de calor es relevante sólo en El modo dinámico. La página de la pérdida de calor contiene la pérdida de calor parámetros que caracterizan a la cantidad de calor perdido a través de la pared del vaso. En el grupo de pérdida de calor de modelo, se puede elegir entre una simple o Modelo de pérdida de calor detallada o ninguna pérdida de calor a través del vaso paredes.


61

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Modelo simple El modelo simple le permite especificar la pérdida de calor directamente, o haber la pérdida de calor a partir de los especificados valores: 

Buen valor U



Temperatura ambiente

El área de transferencia de calor, A, y la temperatura del fluido, Tf, son calculado por HYSYS usando la siguiente ecuación: 𝑄 = 𝑈𝐴(𝑇𝑓 − 𝑇𝑎𝑚𝑏 ) Para un portátil, los parámetros disponibles para el modelo simple aparecen en la siguiente figura.

Figura Nº76: Perdida de calor Los parámetros simples de pérdida de calor son los siguientes: 

En general Coeficiente de Transferencia de Calor



Temperatura ambiente



En general Área de Transferencia de Calor



Flujo de calor

El flujo de calor se calcula como sigue: 𝐹𝑙𝑢𝑗𝑜 𝑑𝑒 𝑐𝑎𝑙𝑜𝑟 = 𝑈𝐴(𝑇𝑎𝑚𝑏 − 𝑇) Donde 

U




A

: área de transferencia de calor



Tamb

: Temperatura de ambiente



T

:


62

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 El flujo de calor se define como el calor que fluye en el recipiente. El calor área de transferencia se calcula a partir de la geometría del vaso. La temperatura ambiente, Tamb, y el coeficiente global de transferencia de calor, U, se puede modificar los valores por defecto que se muestran en rojo.

Ejemplo N° Se requiere enfriar una corriente de agua de 90°C a 40°C a 14.7 psi, se requiere calcular la energía que se libera, y de igual manera calentar dicha corriente y elevar a 90°C y calcular la perdida de energía.

Figura Nº77: Esquema de simulación Horno (Furnace)

El calentador dinámico Fired operación (Horno) realiza la energía y balances de materia para modelar un tipo de calentador de fuego directo horno. Este tipo de equipo requiere una gran cantidad de calor de entrada. El calor se genera por la combustión de combustible y transferido a procesar streams. Un diagrama esquemático simplificado de un calentador de fuego directo se ilustra en la siguiente figura.


63


Figura Nº78: Esquema de un horno (Zona radiante, convectiva, zona ecomizador)

En general, un calentador encendido11 se puede dividir en tres zonas: 

Zona radiante



Zona convectiva



Zona del economizador

Las conexiones de corriente allowsmultiple operación calentador encendido en el lado de los tubos en cada zona y economizador opcional, y selecciones de la zona de convección. La operación incorpora una sola modelo de quemador, y una sola entrada y la salida en el gas de combustión de alimentación lado. Las siguientes son algunas de las principales características de la dinámica Funcionamiento del calefactor usado: 

Conexión flexible de fluido de proceso asociado en cada Zona calentador encendido. Por ejemplo, la zona radiante y convectivo zona o economizador. Calentador de Fired Diferente configuraciones pueden ser modelados o personalizar mediante tee, mezclador, y operaciones unitarias intercambiador de calor.



Una opción de indicación de presión - flujo en cada lado y pasar realista modelos fluyen a través de calentador encendido funcionamiento de acuerdo con el gradiente de presión en la totalidad de red de la presión de la planta. Posible inversión del flujo Por lo tanto, las situaciones pueden ser modelados.

Para definir el número de zonas que requiere el calentador encendido, introduzca el número en # External Pases campo en Conexiones la página de la ficha Diseño. 11


64


Un cálculo inclusivo integral calor radiante, convección y conducción de transferencia de calor en la zona radiante nos permitirá predecir la temperatura del fluido del proceso, Despedido Temperatura del gas de temperatura de la pared del calentador y de combustión.



Un modelo dinámico que representa energía y materiales atracos en cada zona. La transferencia de calor en cada zona depende de las propiedades de los gases de escape, tubo y Fired Propiedades de la pared del calentador, propiedades de la superficie de metal, el calor pérdida para el ambiente y la física corriente de proceso propiedades.



Un modelo de combustión que representa mezcla imperfecta de combustible, y permite que la llama de encendido automático o extinguido basado en la disponibilidad de oxígeno en el combustible mezcla de aire

REACCION DE COMBUSTION La reacción de combustión en el modelo del quemador del calentador encendido realiza puro de hidrocarburos (CxHy ) Cálculos de combustión solamente. La medida de la combustión depende de la disponibilidad de oxígeno que normalmente se rige por la relación de aire a combustible. Relación aire-combustible (AF) se define de la siguiente manera

Puede establecer los límites de la combustión, como la máxima AFand la AF mínimo, para controlar la llama del quemador. La llama no se puede encender si el aire calculado a combustible cae por debajo del aire mínimo especificado para combustible. El aire mínima para combustible y el máximo de aire a combustible se puede encontrar en la Página Parámetros de la ficha Diseño. El calor liberado por el proceso de combustión es el producto de velocidad de flujo molar, y el calor de formación de los productos menos el calor de formación de los reactivos a la combustión la temperatura y la presión. En la operación de la unidad de calentador encendido, un conjunto reacción tradicional para las reacciones de combustión no es requerida. Usted puede elegir los componentes combustibles (el hidrocarburos y de hidrógeno) para ser considerado en la combustión reacción. Usted puede ver el mixingefficiency de cada combustible componente en la página de parámetros de la ficha Diseño.

TRANSFERENCIA DE CALOR Los cálculos de transferencia de calor calentador encendido se basan en la energía saldos de cada zona. El lado de la carcasa del calentador encendido contiene cinco atracos: 

tres de la zona radiante EDGAR MARTIN JAMANCA ANTONIO

65


una zona de convección



un atraco zona economizador como se indica anteriormente

Por el lado del tubo, cada flujo de individuo que pasa a través de las zonas respectivas se considera como una sola atraco. Términos de calor importantes que subyacen al modelo de calentador encendido son se ilustra en la siguiente figura. Los términos de calor relacionados con la tuve side se ilustran en la figura siguiente.


66


Tomando zona radiante como un sobre, el siguiente de la energía aplica ecuación de balance

TRANSFERENCIA DE CALOR RADIANTE Para un objeto caliente en una habitación grande, la energía radiante emitida es dado como


67

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 TRANSFERENCIA DE CALOR POR CONVECCIÓN La transferencia de calor por convección participar entre un fluido y una de metal se da en la siguiente

El U realmente varía con el flujo de acuerdo con el siguiente flujo de U relación si se utiliza este método de flujo escamado:

La relación de flujo de masa en el tiempo t para hacer referencia de flujo de masa es también conocida como factor de flujo reducido. El factor de escala caudal mínimo es el valor más bajo, lo que se prevé que la proporción con bajo caudal región. Para la operación del calentador Despedido, el caudal mínimo a escala factor puede ser expresada sólo como un valor positive Por ejemplo, si el factor de escala caudal mínimo es 0,001 (0,1%), cuando se logra esta relación de flujo de masa, los Uusedstays como un valor constante. por lo tanto

CAÍDA DE PRESIÓN La caída de presión a través de cualquier paso en la unidad de calentador encendido la operación se puede determinar en una de dos maneras: 

Especificar la caída de presión - delta P.



Definir una relación de flujo de presión para cada pasada por especificando un valor k

Si se elige la opción de flujo de presión para la caída de presión determinación en el paso calentador encendido, un kvalue se utiliza para relacionar la caída de presión por fricción y el EDGAR MARTIN JAMANCA ANTONIO

68

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 flujo molar, Fthrough la Despedido del calentador. Esta relación es similar a la válvula en general ecuación:

Esta ecuación de flujo de caída generalizada de presión a través de la usesthe Despedido pase calentador sin ninguna contribución carga estática. la cantidad, (P1-P2) se define como la pérdida de presión por fricción que se utiliza para el "tamaño" del flujo. El kvalue se calcula sobre la base de dos criterios: 

Si el flujo de la islarger sistema que el valor en kref (k flujo de referencia), el kvalue permanecen sin cambios. es recomienda que el flujo kreference se toma como 40% de diseño de estado estacionario flowfor mejor flujo de presión la estabilidad en el rango de caudal bajo.



Si el flujo del sistema es menor que el kref , La K valor viene dado por: Donde: Factor = valor está determinado por HYSYS internamente para tener en consideración la relación gota flujo y la presión para las regiones de bajo flujo. El efecto de kref es aumentar la estabilidad mediante el modelado de un más relación lineal entre el flujo y la presión. Este es también más realista con bajo caudal


69

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Heat Exchanger

El intercambiador de calor se realiza de energía y material de doble cara equilibrar los cálculos. El intercambiador de calor es muy flexible, y puede calcular temperaturas, presiones, flujos de calor (incluyendo pérdida de calor y la fuga de calor), los flujos de corriente de material, o UA. En HYSYS, usted puede elegir el Intercambiador de Calor Modelo para su análisis. Sus opciones incluyen un diseño de análisis Punto Final modelo, un ideal (Ft = 1) contracorriente modelo de diseño ponderado, un método de clasificación de estado estacionario, y un método de clasificación dinámica para utilizar en simulaciones dinámicas. El método de clasificación dinámica es disponible como un modelo básico o detallado, y también puede ser utiliza en el modo de estado estacionario por calificar Intercambiador de calor. La unidad operación también permite el uso de terceros Intercambiador de calor métodos de diseño a través de OLE extensibilidad Los cálculos de intercambiadores de calor se basan en balances energéticos para los fluidos caliente y frío Estado estacionario: En las siguientes relaciones generales, el fluido caliente suministra el calor Intercambiador deber para con el fluido frío

𝐵𝑎𝑙𝑎𝑛𝑐𝑒 𝐸𝑟𝑟𝑜𝑟 = (𝑀𝑐𝑜𝑙𝑑 [𝐻𝑜𝑢𝑡 − 𝐻𝑖𝑛 ]𝑐𝑜𝑙𝑑 − 𝑄𝑙𝑒𝑎𝑘 ) − (𝑀ℎ𝑜𝑡 [𝐻𝑖𝑛 − 𝐻𝑜𝑢𝑡 ]ℎ𝑜𝑡 − 𝑄𝑙𝑜𝑠𝑠 ) Donde: 

M

: caudal másico del fluido



H

: Entalpia



Qleak : fuga de calor



Qloss : perdida de calor



Bal. Er. :a especificación del intercambiador de calor que es igual a cero para la maypria de las aplicaciones.



Hot,cold: fluidos calientres y frios



In,out : corriente de entrada y de salida.

El calor total transferido entre los lados del tubo y concha (Intercambiador de calor deber) se puede definir en términos de la global coeficiente de transferencia de calor, el área disponible para intercambio de calor, y el registro de diferencia de temperaturas media

𝑄 = 𝑈𝐴∆𝑇𝐿𝑀 𝐹𝑡

Donde: EDGAR MARTIN JAMANCA ANTONIO

70


U




A

: Superficie del área disponible de transferencia de calor



∆𝑇𝐿𝑀

: Diferencia logarítmica de temperatura (LMTD)



𝐹𝑡

: factor de corrección LMTD

El coeficiente de transferencia de calor y el área de la superficie son a menudo combinados por conveniencia en una sola variable denominada UA. La LMTD y su factor de corrección se definen en la sección Rendimiento

Caída de presión La caída de presión del intercambiador de calor se puede determinar en una de tres maneras: • Especificar la caída de presión. • Calcular la caída de presión basado en el intercambiador de calor la geometría y la configuración. • Definir una relación de flujo de presión en el intercambiador de calor por especificando un valor k. Si se elige la opción de flujo de presión para la caída de presión determinación en el intercambiador de calor, un valor k se utiliza para relacionar la pérdida de presión por fricción y fluir a través del intercambiador. Este relación es similar a la ecuación general de la válvula

𝑓 = √𝑑𝑒𝑛𝑠𝑖𝑑𝑎𝑑 ∗ 𝑘 √𝑃1 − 𝑃2 Esta ecuación de flujo de caída generalizada de presión a través de la usesthe Intercambiador de calor y sin ninguna contribución de cabeza estática. la cantidad, P1 - P2 , Se define como la pérdida de presión por fricción que se utiliza para el "tamaño" del Intercambiador de calor con un valor k

Se muestra la ventana de propiedades, Intercambiador de calor


71


El UpdateButton le permite actualizar el intercambiador de calor cálculo en el modo dinámico. Por ejemplo, si usted hace un cambio de configuración en el intercambiador de calor, haga clic en este botón para restablecer las ecuaciones aroundthe intercambiador de calor antes de ejecutar el cálculo de simulación en el modo dinámico

Los principales supuestos del modelo son las siguientes: 

En general, el coeficiente de transferencia de calor, U es constante.



Los calores específicos de los flujos tanto de carcasa y tubos secundarios son constante.

El modelo de Punto Final trata a las curvas de calor tanto para calor Lados Intercambiador como lineal. Para problemas simples, donde no es cambio de fase y los IPC relativamente constante, esta opción puede ser suficiente para modelar su intercambiador de calor. Para el calor no lineal problemas de flujo, el modelo ponderado se debe utilizar en su lugar. Los siguientes parámetros están disponibles cuando el punto final se selecciona el modelo

El modelo ponderado es un excelente modelo para aplicar a los no-lineal problemas de la curva de calor tales asthe cambio de fase de pura componentes en uno o ambos lados intercambiador de calor. con el Modelo ponderado, las curvas de calentamiento se divide en intervalos, y un balance de energía se realiza a lo largo de cada intervalo. A LMTD y UA se calculan para cada intervalo en la curva de calor, y sumados para calcular el intercambiador general UA. El modelo ponderado sólo está disponible para los contra-corriente intercambiadores, y es esencialmente la energia y balance de materiales modelo. Las configuraciones de geometría que afectan al pie factor de corrección no se tienen en cuenta en este modelo. EDGAR MARTIN JAMANCA ANTONIO

72

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Cuando se selecciona el modelo ponderado, la página de Parámetros aparece como se muestra en la figura siguiente


73

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Simulación de reactores: Reactor de conversión El reactor de conversión es un recipiente en el que se llevan a cabo reacciones de conversión. Sólo se puede adjuntar conjuntos de reacción que contienen las reacciones de conversión. Cada reacción en el conjunto continúa hasta que la conversión especificada se alcanza o reactivo limitante se agota. Ingreso de reacción:

Ventana de resultados de reacción

El grupo Resultados Reactor Resumen muestra los siguientes resultados para un reactor de convergencia: Resultado Posición

% de conversión Componente base

Descripción Muestra la posición actual de la reacción. Para múltiples reacciones de menor rango se producen en primer lugar. Cuando hay múltiples reacciones en conjunto de reacción, HYSYS clasifica automáticamente las reacciones. Una reacción con un menor rango de valor que ocurra primero. Cada grupo de reacciones de igual rango puede tener una conversión global especificado entre 0% y 100% Muestra el porcentaje del componente base de la corriente de alimentación, que ha sido consumido en la reacción. El reactivo para el que basa el cálculo de la conversión.


74

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Rxn Extensión

Muestra el consumo de tasa molar del componente de base en la reacción dividido por su coeficiente estequimetrico aparecido en la reacción.

Equilibrio de reacción Cuando se selecciona el botón de opción de balance de reacción, la opción de balance de reacción proporciona un componente global resumen para el reactor de conversión. Todos los componentes que aparecer en el envase de líquidos se muestran aquí.

Los valores aparecen después de la solución del reactor que ha convergido. La tasa de flujo de entrada total que reacciono total y el total de tasa de flujo de salida para cada componente se proporcionan sobre una base molar. Los valores negativos indican el consumo de un reactivo, mientras los valores positivos indican la aparición de un producto.

Ejemplo N°… Se requiere encontrar la masa necesaria de aire para obtener una reacción completa y consumir todo el metano al 100%.

1𝐶𝐻4 + 2𝑂2 → 1𝐶𝑂2 + 2𝐻2 𝑂


75


Utilizamos Adjuste o realizamos un balance químico a mano y obtenemos que se requiere aprox. 400kg/h de oxígeno para consumir todo el metano.

Ejemplo N°.. Composición de la corriente Gas Combustible Componente Metano (C4) Etano (C2) Propano (C3) i-Butano (iC4) n-Butano (nC4) i-Pentano (iC5) n-Pentano (nC5)

Fracción molar 0.8837 0.1032 0.0015 0.0005 0.0007 0.0002 0.0003 EDGAR MARTIN JAMANCA ANTONIO

76

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Agua (H2O) Oxigeno (O2) Dióxido de Carbono (CO2) Nitrógeno (N2)

0.0000 0.0000 0.0024 0.0075

Exceso de aire O2%: 2% Vol. O2 en base seca.

Reacción de conversión: 100% 1𝐶𝐻4 + 2𝑂2 → 1𝐶𝑂2 + 2𝐻2 𝑂

𝑅𝑥𝑛 − 1

1𝐶2 𝐻6 + 7⁄2 𝑂2 → 2𝐶𝑂2 + 3𝐻2 𝑂

𝑅𝑥𝑛 − 2

1𝐶3 𝐻8 + 5𝑂2 → 3𝐶𝑂2 + 4𝐻2 𝑂

𝑅𝑥𝑛 − 3

1𝐶4 𝐻10 + 13⁄2 𝑂2 → 4𝐶𝑂2 + 5𝐻2 𝑂

𝑅𝑥𝑛 − 4 (𝑖𝐶4 )

1𝐶4 𝐻10 + 13⁄2 𝑂2 → 4𝐶𝑂2 + 5𝐻2 𝑂

𝑅𝑥𝑛 − 5 (𝑛𝐶4 )

1𝐶5 𝐻12 + 8𝑂2 → 5𝐶𝑂2 + 6𝐻2 𝑂

𝑅𝑥𝑛 − 6 (𝑖𝐶5 )

1𝐶5 𝐻12 + 8𝑂2 → 5𝐶𝑂2 + 6𝐻2 𝑂

𝑅𝑥𝑛 − 7 (𝑛𝐶5 )

PDF Final


77


Reactor de Equilibrio

El reactor de equilibrio es un recipiente que los modelos de reacciones de equilibrio. Las corrientes de salida del reactor están en un estado químico y equilibrio físico. El conjunto de reacciones que usted adjunta al reactor de equilibrio puede contener un número ilimitado de reacciones de equilibrio, que son simultáneamente o secuencialmente resuelto. Ni los componentes ni la mezcla de proceso tienen por qué ser ideal, ya que HYSYS puede calcular la química la actividad de cada componente en la mezcla a base de la mezcla y fugacidades componentes puros. También puede examinar la conversión real, la base componente, la constante de equilibrio, y el grado de reacción en el conjunto de selección de la reacción. La conversión, la constante de equilibrio y la extensión están calculados en base en la información de reacción de equilibrio que ya ha proporcionado cuando se creó el conjunto de reacción. Página de detalles La página de detalle consiste primordial para los botones de relación: 

Stoichiometry



Basis



Ln[K]



Table

Estequimetria Cuando se selecciona el botón de relación estequimetria, aparece la información estequimetrico del grupo. El grupo estequimetria le permite ver la formula estequimetrica de la reacción seleccionada en la lista desplegable de reacción.

El error de balance (para la estequimetria de la reacción) y el calor de reacción (calor de reacción a 25°C) se muestran también para la reacción actual. EDGAR MARTIN JAMANCA ANTONIO

78

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Basis Cuando seleccionas el botón Basis, muestra el grupo básico de la reacción.

El grupo base le permite ver o editar (a nivel local) varias informacines para cada reacción en el conjutno de reacción que incluye la: 

La base para los cálculos de equilibrio



Fase en la que se produce la reacción



Enfoque de temperatura de la composición de equilibrio

El intervalo de temperatura para la constante de equilibrio, y el fuente para el calculo de la constante de equilibrio es también se muestra.

Keq Cuando se selecciona el botón de radio Keq, el grupo (Keq) Ln y K aparece la tabla.

El grupo Ln(Keq) muestra la relación Ln(Keq) que puede variar dependiendo de la Ln(K) valor de origen seleccionado para la reacción. Cuando se selecciona el botón de radio Ln(Keq) la ecuación en el Ln(K) grupo fuente, los parámetros de la constante de equilibrio aparece en la ecuación, estos valores se especifican o bien cuando la reacción se crea o se calculan en HYSYS. Si un fijo se presto constante de equilibrio, que se muestra aquí. EDGAR MARTIN JAMANCA ANTONIO

79

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Caulquiera de los parámetros en el grupo Ln(K) esta ecuación se puede modificar en esta pagina, los cambios realizados en los parámetros solo afecte a la reacción seleccionado en el reactor actual. Después de un cambio se ha hecho, usted puede tener y HYSYS devolverá el original valor calculado seleccionando el uso por defecto apropiado. Approach Cuando se selecciona el botón de opción Approach, el grupo fraccional de enfoque y el grupo de enfoque de temperatura aparecen:

Para cada reacción en el conjunto de reacción un efonque fraccional se proporciona la ecuación como una función de la temperatura. Calquiera de los parámetros de la ecuación de aproximación% se pueden modificar en esta pagina. Los cambios realizados en los parámetros solo afectan a la seleccionada reacción en el reactor. Después del cambio actual que se ha hecho usted puede tener que HYSYS devuelve el valor calculado original seleccionador el apropiado uso Defaultcheckbokx. Puede editar una reacción haciendo clic en el ver “Reaction Button”, apareceera la ventana de propiedades para la reacción resaltado. Resultados


80


Cuando se selecciona el botón de opción de balance de reacción, ofrece un resumen global de componentes para el reactor de equilibrio. Todos los componentes que aparecen en la lista de componentes relacionados con el paquete de fluido se muestra aquí. Los valores aparecen después de la solución del reactor ha convergido, la tasa global de flujo de entrasa, la tasa total que reacciono, y la tasa total de flujo de salida para cada componente se proporcionan sobre una base molar. Los valores negativos indican el consumo de un reactivo, mientras que los valores positivos indican la aparición de un producto. Resultado Actual % conversion

Componente base Constante de Equlibrio

Descripcion Muestra el porcentaje de compoente de base en la alimenteacion corriente(s) que ha sido consumido en la reacción. La conversión real se calcula como el porcentaje del componente base que se consume en la reacción. 𝑁𝐴 − 𝑁𝐴𝑜𝑢𝑡 𝑋 = 𝑖𝑛 𝑥100% 𝑁𝐴𝑖𝑛 Donde: X : Conversion real % 𝑁𝐴𝑖𝑛 : Componente del caudal de base en el reactor 𝑁𝐴𝑜𝑢𝑡 : Caudal del componente base (mismas condiciones que el entrada de índice) del reactor. El reactivo al que se aplica la conversión La constante de equilibrio se calcula a la temperatura del reactor por la ecuacuin siguiente: 𝐵 𝐿𝑛𝐾 = 𝐴 + + 𝐶 𝐿𝑛𝑇 + 𝐷𝑇 𝑇 Donde: T : Temperatura del reactor, K A,B,C,D : parámetros de la ecuación Los cuatros parámetros en la ecuación anterior se calcula como el HYSYS si no se especifican durante la instalacon de la reacción de equilibrio. Los cuatro parámetros para cada ecuación de equilibrio se enumeran en la pagina de Rxn Ln(K)


81

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Rxn Extns.

Lista el consumo de tasa molar del componete de base en la reacción divido por su coeficiente estequimetrico aparecido en la pagina de reacción.

Reacciones de equilibrio incluidas en HYSYS

Ejemplo

Reactor CSTR,

Reactor PFR

CSTR El CSTR es un recipiente en el que Kinetic, heterogénea catalítica, y reacciones Puntúa simples se pueden realizar. La conversión en el reactor depende de la rateexpression de las reacciones asociado con el tipo de reacción. Se supone que la corriente de entrada para ser perfectamente (y de forma instantánea) mezclado con el material ya en el reactor, de modo que la composición de la corriente de salida es idéntica a la de los contenidos del reactor. Teniendo en cuenta el reactor


82

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 volumen, una velocidad consistente expressionfor cada reacción y la estequiometría de la reacción, el CSTR calcula la la conversión de cada componente de entrar en el reactor. En la ficha reacciones, se puede seleccionar una reacción fijada para la peración. También puede ver los resultados del reactor resuelto incluyendo la conversión real del componente de base. la conversión real se calcula como el porcentaje de la base componente que se consume en la reacción.

𝑋=

𝑁𝐴𝑖𝑛 − 𝑁𝐴𝑜𝑢𝑡 𝑥 100% 𝑁𝐴𝑖𝑛

Donde  

X 𝑁𝐴𝑖𝑛

: conversión real : componente caudal de base en el reactor



𝑁𝐴𝑜𝑢𝑡

: caudal componente base (misma base que la tasa de entrada) del reactor.

Reactor de equilibrio El reactor de Equilibrio es un recipiente que los modelos de equilibrio reacciones. Las corrientes de salida Ofthe del reactor están en un estado de química y equilibrio físico. El conjunto de reacción que usted adjuntar a la Equilibrium reactor puede contener un número ilimitado número de reacciones de equilibrio, que son simultáneamente o secuencialmente resuelto. Ni los componentes ni la mezcla proceso tiene por qué ser ideal, ya HYSYS puede calcular la química la actividad de cada componente en la mezcla a base de la mezcla y fugacidades componentes puros. También puede examinar la conversión real, la base componente, la constante de equilibrio, y el grado de reacción para cada reacción en el conjunto selectedreaction. La conversión, la constante de equilibrio y la extensión están calculados en base en la información de reacción de equilibrio que ya ha proporcionado cuando se creó el conjunto de reacción

Las reacciones de los reactores de Gibbs El reactor de Gibbs calcula las composiciones que salen de tal manera que la fase y químicas equilibrios de las corrientes de salida son alcanzado. Sin embargo, el Gibbs reactor no necesita para hacer utilizar de una estequiometría de la reacción especificada para calcular la salida composición de la corriente. La condición de que la energía libre de Gibbs de el sistema de reacción está en un mínimo en el equilibrio se utiliza para calcular la composición de la mezcla producto. Al igual que con la Equilibrio reactor, ni los componentes puros ni la reacción mezcla se supone que se comportan muy bien. EDGAR MARTIN JAMANCA ANTONIO

83


Reactor de flujo de piston (PFR) El PFR (Plug Flow Reactor o Reactor tubular) generalmente consiste en un banco de tubos cilíndricos o tubos. El campo de flujo es modela como flujo de pistón, lo que implica que el flujo es radial isotrópica (sin masa o energygradients). Esto también implica que la mezcla axial es insignificante.

Como los reactivos fluyen de la longitud del reactor, que están consumido continuamente, por lo tanto, hay una variación axial en concentración. Dado que la velocidad de reacción es una función de la concentración, la velocidad de reacción también varía axialmente (a excepción de orden cero reacciones). Para obtener la solución para el PFR (perfiles axiales de composiciones, temperatura, y así sucesivamente), el reactor se divide en varios subvolúmenes. Dentro de cada volumen secundario, la velocidad de reacción es considera que es espacialmente uniforme. Un balance molar se realiza en cada subvolumen j EDGAR MARTIN JAMANCA ANTONIO

84


𝐹𝑗0 − 𝐹𝑗 + ∫ 𝑟𝑗 𝑑𝑉 =

𝑑𝑁𝑗 𝑑𝑡

Debido a que la velocidad de reacción se considera espacialmente uniforme en cada uno subvolumen, el tercer término se reduce a RJV. En el estado estacionario, la lado derecho de este equilibrio equalszero, y la ecuación se reduce a 𝐹𝑗 = 𝐹𝑗0 + 𝑟𝐽 𝑉


85


SIMULACIÓN DE COLUMNAS DE DESTILACIÓN Columnas de destilación corta Uno de los procedimientos más usualmente utilizados para obtener estimativos simplificados de numero de etapas teóricas requeridas en una separación por destilación es el propuesto por Fenske, Underwood y Gililand. Correlación de Gililand 𝑁 − 𝑁𝑚 𝑅 − 𝑅𝑚 0.5688 = 0.75 ∗ [1 − ( ) ] 𝑁+1 𝑅+1 Donde: 

Nm

: número de etapas a reflujo total



Rm

: la relación de reflujo mínimo



R

: relación de reflujo de operación

Ecuación de Fenske para calcular el número mínimo de platos

𝑁𝑚 =

𝑋 𝑋 𝑙𝑛 [(𝑋𝐷,𝐿𝐾 ) ( 𝑋𝑊,𝐻𝐾 )] 𝐷,𝐻𝐾

𝑊,𝐿𝐾

𝑙𝑛𝛼𝐿𝐾/𝐻𝐾

Donde: 

D

: destilado



W

: Fondos



LK

: componente como clave liviano



HK

: componente como clave pesado



𝛼

: volatilidad

Ecuación de Underwood para calcular la relación de reflujo mínimo Para mezclas multicomponentes con volatilidad relativa constante, la relación de reflujo mínimo en una columna de destilación puede calcularse mediante dos ecuaciones: 𝑛

∑ 𝑖=1

𝛼𝑖 𝑋𝑖,𝐹 =1−𝑞 𝛼𝑖 − 𝜃

Donde: EDGAR MARTIN JAMANCA ANTONIO

86


n

: número de componentes



q

: calidad del alimento

La segunda ecuación se utiliza el valor 𝜃, calculado con la ecuación nº ….., para estimar la relación de reflujo mínimo. 𝑛

∑ 𝑖=1

𝛼𝑖 𝑋𝑖,𝐷 = 1 + 𝑅𝑚 𝛼𝑖 − 𝜃

Relación de reflujo12 La experiencia ha demostrado que el valor óptimo de la relación de reflujo se encuentra en un intervalo dado por:

1.03 < 𝑅⁄𝑅 < 1.3 𝑚

Ejemplo N° Se quiere destilar la siguiente composición a las condiciones siguientes T 200ºF, presión 100 psia y 1300 lbmole/h Composición Etano Propano i-Butano n-Butano i-Pentano n-Pentano n-Hexano

Fracción Molar 0.0148 0.7315 0.0681 0.1462 0.0173 0.0150 0.0071

Se determina las corrientes de entrada y los productos, para calcular el número de platos y el plato de la alimentación.

Una regla de diseño sugiere que se asigne, para una separación deseada, una relación de reflujo de 1.2 veces la mínima 12


87


Se determina los productos ligeros y pesados en el destilado y en los fondos

El reflujo externo se calcula mediante la multiplicación de 0.997*1.3, y de inmediato la simulación calculara el número de platos, plato de alimentación y las demás variables de operación.


88

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Columna de destilación simplificada

Se presente 3 tipos de equipos: 

Destilación simple (absorbedor)



Destilación con reboiler



Destilación con condensador.

Para la convergencia, del siguiente equipo, se debe determinar el número de platos de la columna, el número de platos en el cual deberá ingresar la corriente de alimentación, configurar las corrientes de productos y de energía, las condiciones operaciones, y determinar los grados de libertad en la simulación para que ello converja.

Página de configuración del reboiler


89


Ejemplo: Se requiere destilar un mezcla de 100 kgmol/h de agua y 50 kmolg/h de etanol a condiciones normales, determinar la composición del destilado y de los fondos, la columna es de 15 platos y la alimentación ingresa en el plato 8.

Comportamiento de la composición en los platos

PFD del proceso de destilación EDGAR MARTIN JAMANCA ANTONIO

90

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Ambiente de la columna de destilación

Ventana de configuración del reboiler

Ventana de configuración del condensador total.


91

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Dentro del ambiente de la configuración de la columna de destilación, se encuentra una paleta de herramientas especialmente para esta operación, cuenta con reboiler, condensadores total, columna absorbedor, condensador parcial, condensador reflux, Cooler, Heater, bombas, válvulas, separadores, intercambiadores, mixer, tee.

Ejercicio N° El propano y propileno son muy difíciles de separa uno de otro, ya que son componentes con puntos de ebullición cercanos. No obstante, la destilación a presión elevada es una tecnologia

Columnas de Absorbedor El absorbedor tiene cuatro corrientes de contorno y así requiere cuatro especificaciones de flujo de presión. Una especificación de presión siempre se requiere para la corriente de producto líquido que sale de la parte inferior de la columna. Una segunda especificación de presión debe añadirse a la del producto de vapor de la columna, con las dos corrientes de alimentación que tiene especificaciones de flujo

Laa única operación de la unidad contenida en el absorbedor es la sección de la bandeja, y las únicas corrientes son de vapor de cabeza y los productos líquidos de fondo.

En conclusión la absorción de gases es una operación en la cual una mezcla de gases se pone en contacti con un liquido, a fin e disolver de forma selectiva uno o mas componente en el gas y


92

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 obtener un solución de estos en el liquidi. En HYSYS están disponibles columnas de absorción de gases.

ejemplo El CO2 es absorbido en carbonato de propileno (propylenearbonato), la corriente de gas de entrada es 20% mol de CO” y 80% mol de metano, este fluye a razón de 2m3/s y la columna funciona en 60°F y 60.1 atm; el flujo de solventee de la entrada es 2000 kmol/h.


93

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 MODULO VI: CASOS DE ESTUDIO Caracterización de crudos de refinerías Objetivos 

Obtener el rendimiento y determinar las propiedades más importantes de cada uno de los cortes establecidos para cualquier tipo de Crudo o mezclas de diferentes tipos de Crudo.



Estimar las principales condiciones de operación en la Columna de Destilación de interés, para luego completar el balance de materia y energía correspondiente

Bases de simulación 

Modelo Termodinámico: Peng-Robinson



Método de Corte de crudo: Autocortes



Método de Extrapolación: Mínimos cuadrados



Especificación de Cortes Corte Gases de Refinería Gasolina (NL) Solvente3 (NP) Kerosene Diesel

PFE (°C) 32 150 210 300 385

Columna de fraccionamiento de crudo


94


Componentes principales  Carga a la columna  Zona Flash  Platos (Números de platos)  Pumparound (reflujo externos)  Reflujo de tope  Sistema de condensación de tope  Stripper (Despojadores)  Vapor despojante Bases de simulación Corriente de ingreso a la columna  

Carga Requerimiento de vapor (columna y despojadores)

Especificaciones en la columna      

Numero de platos Plato de alimentación Corrientes principales de extracción (Tope y fondo) Presión a la salida del sistema de condensación Perdida de presión en el sistema de condensación Presión en el fondo de la columna. EDGAR MARTIN JAMANCA ANTONIO

95

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Especificaciones mínimas de equipos complementarios Pumparound  Plato de extracción  Plato de retorno Stripper de Vapor  Plato de extracción  Plato de retorno Convergenvia y monitoreo Los parámetros de convergencia y monitoreo se definirán desde el monitor hysys.          

Extracción de los productos de destilación Flujos de pumparound Variación de temperatura en el pumparound Overflash GAP entre cortes Temperatura de cortes Temperatura de los platos correspondientes a los productos obtenidos Temperatura de las corrientes de extracción y retorno correspondiente a los pamparound Calidad de los cortes (PIE, 10%, 50%, 90%, PFE, Flasf point) Duty del sistema de condensación y de los pamparound

Datos de partida Carga de la UDP Tipo de Crudo Crudo A Carga (BPD) 20, 000 T (°F) 630 P (psig) 23 UDP Presion cond. (psig) 8 Variacion de presion 6 (psi) Pesion fondo (psig) 19 Strippers N° Platos 6


96


Equipo Flujo (lb/h) T (°F) P(psig)

UDP 800 350 100

Vapor despojante Stripper – NP Stripper – Kerosene 330 100 350 350 100 100

Stripper – Diesel 100 350 100

Datos de convergencia Gas Refinería (GPM) Naphta Ligera (GPM) Naphta Pesada (GPM) Kerosene (GPM) Diesel (GPM) Reflujo PA (GPM) ∆T PA °F) Flash NP (°C) PIE NP (°C) 50% NP (°C) 90% NP (°C) PFE NP (°C)

40 149 177 199 210

T (°F) plato 27

350

0.5 60 40 52 66 170 200


97

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 T (°F) plato 20 T (°F) plato 12

400 600

PFE Kerosene (°C) 90% Diesel (°C) PFE Diesel (°C)

300 357 385

Crudo A: °API@60°F Viscosidad T(°F) cSt T(°F) cSt

: 23.6 Cinematica 100 59.87 122 33.06

TBP %V destilado 0.98 6.92 13.14 16.63 20.55 28.18 32.41 41.69 46.09 52.93 61.95 73.93

°F 82 200 300 350 400 500 550 650 700 800 900 1050

Comp. Ligeros Methane Ethane Propane i-butane n-butane i-pentane n-pentane

%V 0 0 0.24 0.16 0.58 0.59 0.69

BSW Agua

0


98

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Eficiencia de Hornos o calderas

Objetivos  Construir un paquete de reacciones en el entorno básico del simulador  Utilizar un reactor de combustión, divisor de corrientes entre otros equipos y la herramienta Spreadsheet, así como el operador lógico “Adjust”.  Calcular la eficiencia de un horno que utilice gas combustible. Ecuaciones de cálculo de eficiencia:

𝐸= 𝑥𝑟 =

𝐶𝑎𝑙𝑜𝑟 𝑡𝑜𝑡𝑎𝑙 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 − 𝑝𝑒𝑟𝑑𝑖𝑑𝑎 𝑡𝑜𝑡𝑎𝑙𝑒𝑠 𝑥 100 𝐶𝑎𝑙𝑜𝑟 𝑡𝑜𝑡𝑎𝑙 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎

𝐶𝑎𝑙𝑜𝑟 𝑡𝑜𝑡𝑎𝑙 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 − 𝑝𝑒𝑟𝑑𝑖𝑑𝑎 𝑒𝑛 𝑙𝑎 𝑧𝑜𝑛𝑎 𝑟𝑎𝑑𝑖𝑎𝑛𝑡𝑒 𝑥 100 𝐶𝑎𝑙𝑜𝑟 𝑡𝑜𝑡𝑎𝑙 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 𝐸=

(𝐿𝐻𝑉 + 𝐻𝑎 + 𝐻𝑓) − (𝑄𝑠, 𝑟 + 𝑄𝑠, 𝑐) − 𝑄𝑟 𝑥 100 (𝐿𝐻𝑉 + 𝐻𝑎 + 𝐻𝑓)

𝑥𝑟 =

(𝐿𝐻𝑉 + 𝐻𝑎 + 𝐻𝑓) − (𝑄𝑠, 𝑟) − 𝑄𝑟 𝑥 100 (𝐿𝐻𝑉 + 𝐻𝑎 + 𝐻𝑓) 𝑥𝑐 = 𝐸 − 𝑥𝑟

Esquema del horno


99


Condiciones de proceso Variables de proceso     

T, entrada gas combustibles T, entrada de aire T, salida productos del combustible %Vol. O2 en exceso, base seca T, zona de choque (BWT)

77 °F 77 °F 808°F 3.0 1440°F

Composición del gas combustible Metano 0.8832 Etano 0.1052 Propano 0.0014 N2 0.0078 CO2 0.0024 Consideraciones 2.5% de pérdidas de calor al ambiente por radiación Combustión completa de gas natural Solución Dentro de la simulación hay que incluir al Agua y oxigeno que aparecerán durante el proceso mediante las reacciones que se darán, se utilizara el modelo termodinámico PR. Ingreso de los componentes

Ingresando las reacciones – reacción de combustión (100%) 𝟏𝑪𝑯𝟒 + 𝟐𝑶𝟐 → 𝟏𝑪𝑶𝟐 + 𝟐𝑯𝟐 𝑶


100


Se observa con la cantidad de 10000lbmol/h la reacción no converge y se obtiene 18.66% de Oxigeno en la corriente de combustión, y también se observa que no se encuentra en base seca, ya que se encuentra 2.07% de agua. Para ello se utilizara la herramienta lógica “Adjust” para obtener la cantidad de aire suficiente para que en la corriente de producto de combustión se obtenga 3%vol. De oxígeno.


101


Usando la herramienta Adjust se calcula un flujo másico de aire de 3.485x104 lb/h

Para verificar si el producto de combustión se encuentra en base seca se utilizara un Splitter para verificar la composición, esta herramienta nos permite hacer un análisis de transferencia de masa.


102


Se especifica en la corriente base seca 0 Agua

Para converger el Splitter, se especifica que la corriente calcular la temperatura igual (sólo dos productos) e igualar todas las presiones corrientes

Se observa la composición de oxigeno de 3.58% en la corriente de base seca EDGAR MARTIN JAMANCA ANTONIO

103


La temperatura de la corriente de producto de combustión es 3294.2°F, esto se requiere enviar hasta 808°F, para ello utilizaremos un cooler.

Para enfriar la corriente se requiere 29.67 MMBTU/h.

La corriente que se enfrio a 808°F, se tiene que enfriar a 77°F para ello se necesita otro cooler para bajar la temperatura.

La corriente enfriada, se separara mediante un splitter para separar el agua de la corriente, la corriente de agua que se separe tiene que estar a una temperatura de 77°F


104


Siendo Q-101, viene hacer calor de perdida

A la corriente de agua se requiere calcular el calor latente y sensible para realizar el cálculo de la eficiencia del horno


105


Utilizando un Spreadsheet para el cálculo de la eficiencia de calor, esta herramienta nos permitirá extraer los datos de la simulación.


106


La eficiencia del horno es 78.33%

Sistema multiefecto de evaporadores

Realizar un sistema de evaporación triple efecto paralelo, que permita concentrar el jugo en 60°brix, para ello calcular el vapor vivo suficiente para tener la concentración final en °brix. Datos de Partida Jugo mezclado Flujo másico

: 40 TN

Brix

: 15°

Presión

: 14.7 psia

Temperatura : 30°C

El vapor vivo viene de la caldera a una presión de 100 psia y 167°C, pasando por una válvula y reduciendo su presión a 25psia, el vapor se introduce a un sistema de precalentamiento para el jugo mezclado de 48 a 105 y la diferencia ingresa al primer tacho


107



108

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 Planta de producción de etanol

El objetivo del presente trabajo es presentar un programa de simulación para obtener etanol a partir de un sistema de purificación, haciendo el uso de la tecnología de HYSYS, con el cual se puedo dar solución del sistema plantado. El plato óptimo de salida de la corriente “Fusel”, es el plato 20, siendo el componente que contiene una mayor concentración de 1 – Propanol, que es el componente principal de extracción de tal corriente. Se determinó un perfil de flujos molares y másicos de corriente de condensado y rehervidor, en función a la relación de reflujo en un rango de 2300 a 10000. De manera análoga también se determinaron perfiles de calor a agregar al rehervidor y calor a retirar del condensador. Esto con objetivo de en un futuro analizar el reflujo óptimo, teniendo en consideración el aspecto económico.

Solución: Modelado Considerar la columna de absorción de gas mostrada en la Figura Nº 01. Los componentes que entran al fondo de la columna en la corriente de alimentación de gas son absorbidos por la corriente de líquido, de tal manera que la corriente de gas producto (saliendo por el tope de la columna) es más "puro". Las columnas de absorción a menudo contienen "platos" con una corriente de líquido fluyendo a través de los platos; estos platos son a menudo modelados como etapas de equilibrio.


109

SIMULACIÓN DE PROCESOS APLICADO A LA INDUSTRIA USANDO ASPEN HYSYS V8.0 DEFINICIÓN DE VARIABLES Usamos las siguientes definiciones de las variables: 

L = moles de liquido / tiempo



V = moles de vapor / tiempo



M = moles de liquido / etapa



W = moles de vapor / etapa



xij = fracción molar componente i en liquido (etapa j)



yij = fracción molar componente i en gas (etapa j)

ETAPA DE EQUILIBRIO El concepto de una etapa de equilibrio es importante para el desarrollo de un modelo dinámico de la columna de absorción. Una etapa de equilibrio es representada esquemáticamente

1. Balances de materia para cada componente (c ecuaciones para cada etapa): 𝑀𝑖𝑗 = 𝐿𝑗−1 𝑥𝑖,𝑗−1 + 𝑉𝑗+1 𝑦𝑖,𝑗+1 − 𝐿𝑗 𝑥𝑖,𝑗 − 𝑉𝑗 𝑦𝑖,𝑗 = 0 … (1) 2. E- relaciones de equilibrio entre fases para cada componente (c ecuaciones por etapa):

𝐸𝑖𝑗 = 𝑦𝑖𝑗 − 𝐾𝑖𝑗 𝑥𝑖𝑗 … (2)

3. S - Sumatorios de las fracciones molares (una para cada etapa): 𝐶

(𝑆𝑦 )𝑗 = ∑ 𝑦𝑖𝑗 − 1.0 = 0 … (3) 𝑖=1 𝐶

(𝑆𝑥 )𝑗 = ∑ 𝑥𝑖𝑗 − 1.0 = 0 … (4) 𝑖=1


110


4. H- Balance de energía (uno para cada etapa):

𝐻𝑗 = 𝐿𝑗−1 𝐻𝐿𝑗−1 + 𝑉𝑗+1 𝐻𝑉𝑗+1 − 𝐿𝑗 𝐻𝐿𝑗 − 𝑉𝑗 𝐻𝑉𝑗 = 0 … (5)

Frecuentemente, los sistemas de ecuaciones que resultan al aplicar los métodos componente a componente suelen dar lugar a una matriz de coeficientes en forma de matriz tridiagonal, y se resuelven mediante un método de eliminación progresiva que recibe el nombre de algoritmo de Thomas.

A continuación se plantea el sistema de ecuaciones MESH. NOMENCLATURA SISTEMA DE ECUACIONES MESH


111


FLUJO CALOR RETIRADO DEL CONDENSADOR Y CALOR AGREGADO EN EL REHERVIDOR VS RELACIÓN DE REFLUJO

FLUJO MÁSICO DEL CONDENSADOR Y REHERVIDOR VS RELACIÓN DE REFLUJO


112


FIG URA N º 12 – FLUJO MOLAR DEL CONDENSADOR Y REHERVIDOR VS RELACIÓN DE REFLUJO

FLUJO MÁSICO DE 1-PROPANOL EN FUSEL VS. PLATO DE EXTRACCIÓN DE FUSEL


113


Atentamente: Ing. Edgar Martin Jamanca Antonio UNJFSC – Ingenieria Química

https://www.facebook.com/jamancaantonio https://www.youtube.com/channel/UCkchegYyLut6GKaxEONmW5Q


114

Cartilla de Ejercicios

Recommend Documents