Diseño y Economía en Ingeniería Química
Diseño y Economía en Ingeniería Química
SÍNTESIS DE PROCESOS E
INGENIERÍA ECONÓMICA
1
INDICE
CAPITULO I: INTRODUCCION
1.1 Introducción al diseño de procesos químicos 10
1.2 El aspecto creativo del diseño de procesos 11
1.3 El método ingenieril 13
1.4 Desarrollo de un caso base 16
CAPITULO II: ANALISIS ECONÓMICO DE PROCESOS
2.1 Inversión y costo de producción 18
2.2 Componentes de la inversión de capital 22
2.3 Modelos simplificados de inversión, costo y utilidad 28
2.4 Evaluación de la rentabilidad de un proceso 38
CAPITULO III: REGLAS HEURISTICAS Y METODOS APROXIMADOS
3.1 Reglas heurísticas 46
3.2 Métodos aproximados de diseño 53
CAPITULO IV: DISEÑO Y OPTIMIZACION DE SUB-SISTEMAS
4.1 Diseño óptimo de recuperación de calor 103
4.2 Recuperación de calor con producción de calor 105
4.3 Diseño de un sistema de recuperación de un soluto 106
4.4 Alternativas de procesos 112
4.5 Alternativas del diagrama de flujo 112
2
7.4 Limitaciones de equilibrio 140
CAPITULO VIII: SISTEMA DE SEPARACION
8.1 Estructura de separación de componentes 142
8.2 Cálculos aproximados para una separación flash 144
8.3 Recuperación de vapor 146
8.4 Recuperación de líquidos 148
8.5 Otros tipos de separación de componentes 153
BIBLIOGRAFIA 154
3
INTRODUCCIÓN
A través de los conceptos básicos de síntesis y diseño de procesos se enfatiza el aspecto
creativo de la actividad del ingeniero químico. En el texto se analizan conceptos
económicos básicos y reglas heurísticas para sistemas simples que permiten simplificar el
diseño de un proceso. Se introducen el concepto de diseño conceptual de los procesos
químicos, con la meta específica de encontrar el diagrama de flujos óptimo para un
determinado proceso, lo que significa seleccionar las unidades de procesos y las
interconexiones entre estas unidades, así como las condiciones óptimas de diseño y de
operación del proceso. Se introducen también los conceptos de seguridad y salud,
minimización de residuos y el tratamiento de efluentes de la industria de procesos.
4
PARTE I: INGENIERÍA ECONÓMICA Y REGLAS HEURÍSTICAS
CAPITULO I. INTRODUCCIÓN
A. SINTESIS (DISEÑO)
Síntesis de un proceso es definir las entradas y salidas del sistema; o sea, la estructura del
proceso que se requiere para llevar a cabo la transformación deseada de los reactivos
(materias primas) en productos.
ENTRADAS SALIDAS
PROCESO
B. ANÁLISIS (SIMULACIÓN)
ENTRADAS SALIDAS
PROCESO
C. OPTIMIZACIÓN
ENTRADAS SALIDAS
PROCESO
5
Tabla 1.1. Tecnología de diseño de equipos de procesos
TEORÍAS
FENOMENO DE OPERACIONES EQUIPOS PRINCIPALES
TRANSPORTE UNITARIAS
Bombas
Compresoras
Eyectores
Ciclones
Transporte de
Flujos de Fluidos Tuberías
Momentum
Ventiladores
Mezcladores
Tanques separadores
Sopladores
Transferencia de Calor por Aislamiento
Conducción Material refractario
Intercambiadores de calor
Condensadores
Rehervidores
Evaporadores
Transferencia de Calor por
Transporte de
Convección Equipo de recuperación de
Energía calor
Cambiadores de calor
enfriados por aire
Deaireadores
Hornos
Transferencia de Calor por
Calderas
Radiación
Quemadores
Destilación Columnas de destilación
Absorción-desorción Columnas de absorción
Adsorción Columnas de adsorción
Transporte de Masa Extracción Extractores
Ultrafiltración Mallas moleculares
Osmosis Inversa
Electrodiálisis
Transporte Reactores químicos
Simultáneo de
Cinética química Reactores biológicos
Momentum, Calor y
Masa
Las ideas nuevas generalmente son generadas de la química menos frecuentemente por
otras ciencias, y son utilizadas como base para el desarrollo de la tecnología.
6
Tabla 1.2. Principales ecuaciones de transporte
Cantidad Calor Masa Momento
Flux q Na rz
Transporte Molecular
T C A VZ
Fuerza Impulsora
Z Z r
Ley Fourier’s Fick’s Newton’s
Conductividad Difusividad Viscosidad
Propiedad
Térmica
KT DA
Transporte Total
Fuerza Impulsora T CA P
P 1
Relaciones q hT T N A k L C A f gc D
L 2 v2
* Fuerza impulsora en términos de presión parcial y fracción molar también es comúnmente usada.
En cualquier compañía o empresa podemos tratar de generar nuevas ideas con diferentes
propósitos:
Producir una nueva materia que se pueda comprar
Convertir un subproducto no deseado en un producto valioso
Crear un material completamente nuevo (fibras sintéticas)
Encontrar una nueva vía para producir un producto existente (nuevo catalizador)
Explorar un nuevo material de construcción (para operación de alta temperatura o alta
presión)
Como una identificación del tremendo éxito de la ingeniería, podemos notar que sobre el
50 % de los productos vendidos por la mayoría de empresas químicas han sido
desarrollados en las dos últimas décadas (1960 – 1980).
A pesar de este éxito, podemos visualizar que muy pocas ideas tanto para mejorar procesos
existentes o para desarrollar nuevos procesos nos conducen a nuevos productos. De hecho
en la etapa de investigación, un proyecto tiene un chance de 1 – 3 % de llegar a la
comercialización. El diagrama a continuación explica las etapas y su correspondiente
porcentaje de éxito.
100 % 40 – 60 % 10 – 25 % 1–3%
7
Por supuesto, es de esperar que el éxito para modificaciones de procesos sea mayor que
para plantas completamente nuevas, pero usualmente la rentabilidad económica será menor
en ese caso.
Debido a que tan pocas ideas, prueban ser exitosas, (lo mismo sucede en cualquier
actividad creativa), en cualquier metodología de diseño, la evaluación económica en uno
de los componentes más importantes. En otras palabras, se necesitan procedimientos
eficientes para poder descartar ideas sin rentabilidad y también para modificar ideas y
hacerlas rentables, es decir para generar alternativas de modo de que cualquier esfuerzo
posterior en el desarrollo en un proceso sea justificable económicamente.
A continuación, el artista agrega color, sombra y los detalles de los diferentes objetos en la
pintura, aunque continúa el trabajo y puede introducir grandes modificaciones si es que
parecen convenientes. En una forma análoga, el ingeniero utiliza diseños y estimación de
costos más detallados y rigurosos para los equipos de mayor costo, mejorando la exactitud
de los balances preliminares de materia y energía, y agrega más detalles en término de los
equipos pequeños y de bajo costo, que son necesarios para la operación del proceso pero
que no tienen un gran impacto en el costo total de la planta, es decir, bombas, equipos de
8
flasheo, etc. Así, por lo tanto, podemos ver que tanto una pintura como el diseño de un
proceso pasan a través de etapas sucesivas de diseño y evaluaciones detalladas y rigurosas.
DISEÑO ARTE
Proceso Recreativo
PROCEDIMIENTO BOSQUEO A
PRELIMINAR LAPIZ
ALTERNATIVAS DE MODIFICACIONE
PROCESO S PRELIMINARES
Otra cualidad del arte es que nunca existe una solución única a un problema, es decir,
existen diferentes formas de pintar una obra de pintura, tal como existen diferentes
caminos en ingeniería para cumplir el mismo objetivo por esencialmente el mismo costo.
Otra similitud entre el diseño en ingeniería y arte, es que se requiere bastante criterio para
decidir cuánto detalle debe incluirse en las diferentes etapas de una pintura, al igual que
sucede en el diseño de procesos.
9
Por lo visto anteriormente, podemos ver que existen un gran número de tipos de diseños y
estimaciones de costo que pueden llevarse a cabo. Estos pueden variar desde estimaciones
muy simples, rápidas pero no muy exactas, hasta cálculos muy detallados y con la
exactitud que deseamos. Una clasificación puede ser la siguiente:
También es posible determinar los costos relativos para obtener los tipos de diseño
nombrados anteriormente.
De esta tabla podemos ver como los costos suben rápidamente a medida que se incluyen
mayor cantidad de detalles en los cálculos. Obviamente, evitar grandes costos en el diseño,
a menos que se justifiquen económicamente.
Para el caso de un proceso nuevo, donde no existen datos sobre costos, pareciera imposible
desarrollar el tipo de diseño llamado de orden de magnitud, sin embargo, un ingeniero con
10
experiencia puede sobreponerse a esta dificultad, efectuando analogías entre el proceso
nuevo y otros existentes para los cuales él tiene algunos datos.
Para el ingeniero que comienza, y que no tiene suficiente experiencia, es comenzar con los
tres primeros tipos de diseño o una combinación de ellos. Podemos utilizar un estudio de
orden de magnitud para simplificar muchos de los cálculos y limitar nuestra atención sólo a
los equipos mayores que efectuamos un diseño preliminar de un proceso.
A pesar de que el enfoque del diseño y evaluación se dirige a nuevos procesos, esta
metodología se puede utilizar cualquiera de otras disciplinas utilizadas en ingeniería tales
como investigación básica y asistencia técnica. En investigación básica, debemos guiar
nuestros esfuerzos hacia el estudio de las variables que tendrán el mayor impacto
económico en el proceso, y esas variables sólo se identifican desarrollando diseños y
evaluaciones aproximadas.
Llamamos a este primer diseño que desarrollamos, el caso base. Representa nuestra mejor
conjetura de las condiciones y configuración óptimas del proceso. Por supuesto, que al no
tener suficiente experiencia, nuestra primera estimación no será muy buena, pero eso no es
de gran importancia, ya que esperamos considerar diferentes alternativas. En otras
palabras, tendremos que sustituir la experiencia de un ingeniero experto por el esfuerzo
adicional en la forma de mayor cantidad de estimaciones y cálculos.
12
CAPITULO II: ANALISIS ECONÓMICO DE PROCESOS
Un diseño aceptable de una planta debe presentar un proceso que sea capaz de operar bajo
condiciones tales que produzcan una cierta ganancia. Puesto que las utilidades o rentas son:
Es esencial que el ingeniero químico tenga claro cuáles son los costos asociados en un
proceso de la industria química.
Se debe invertir capital para los gastos directos de la planta, tales como materia prima,
mano de obra y equipos. Los gastos indirectos también se deben incluir en un análisis
completo del costo total. Como ejemplo de costos indirectos están los salarios del personal
administrativo, los costos de distribución de productos, etc.
Una inversión de capital se requiere para cualquier proceso de la inversión de capital fijo,
para los equipos físicos y edificios auxiliares de la planta, más el capital de trabajo que
debe estar disponible para el pago de salarios, materias primas y otros gastos imprevistos
que requieren dinero al contado.
La figura 2.1 muestra el concepto del flujo de caja de una operación industrial basada en
un sistema soporte que sirve como fuente de capital o como receptor de capital. Las
entradas al receptor de capital pueden ser en forma de préstamos, captación de acciones y
otro tipo de recursos incluyendo el flujo neto de caja que retorna al receptor de capital de
cada proyecto.
La salida de la fuente de capital se efectúa en forma de inversión de capital para cada una
de las operaciones de la industria, dividendos para las accionistas, pago de deudas y otras
inversiones.
El diagrama de la figura 2.1 muestra un rectángulo que implica la inversión total de capital
que es necesario para iniciar una operación industrial en general. Esta inversión total de
capital incluye todos los fondos necesarios para poner en marcha el proyecto. Esto
13
comprende la inversión de capital fijo directo o inversión de capital de manufactura y el
capital de trabajo, junto con la inversión requerida para todos los servicios auxiliares de
no-manufacturas.
Capital de
trabajo
Inversión
Capital fijo
total de
directo
capital
Pago de
prestamos Capital fijo
indirecto
Otras
inversiones Dividendos accionistas
14
La Inversión Total de Capital es la suma de la Inversión de Capital Fijo más el Capital de
Trabajo y el Costo de Puesta en Marcha.
Es el costo requerido para construir el proceso, es la suma de los Costos Directos (CD) y
de los Costos Indirectos (CI)
Es la suma del costo del material y de obra requerida para construir la planta
completamente, es aproximadamente un 17 – 80% del ICF, consta de:
Corresponde a los costos de instalación de los equipos que se muestran en los diagramas de
flujo, en una localización geográfica específica (límites de la planta); representa
aproximadamente entre 50 – 60% del ICF.
Incluye todos los equipos listados en un diagrama de flujo, además de los equipos de
repuesto. Es aproximadamente entre el 20 – 40% del ICF.
Costos de instalación de todos los equipos del diagrama de flujo, incluyendo soportes
estructurales, aislación y pintura. Representa entre el 7,3 – 26% del ICF o 35 – 45% del
Costo de Compra.
d) Costo de Cañerías
Incluyen los costos directamente relacionados con el proceso pero construidos en una
localización separada de los equipos principales del proceso.
15
1) Edificios: Entre 6 – 20% del ICF o 10 – 70% del costo de copra de los equipos.
Son los costos que no están directamente envueltos con la mano de obra y materiales de la
planta de procesos, entre 15 – 30% de ICF.
16
3) Supervisión de construcción: Contabilidad, compras, etc.
4) Casetas de personal y guardias.
5) Seguridad y servicio médico
6) Permisos, licencias especiales.
7) Impuestos, seguros e intereses.
8) Equipos y herramientas de construcción.
9) Supervisión de construcción: Contabilidad, compras, etc.
10) Casetas de personal y guardias.
11) Seguridad y servicio médico
12) Permisos, licencias especiales.
13) Impuestos, seguros e intereses.
Para compensar por eventos no predecibles tales como tormentas, huelgas, cambios de
precios, cambios pequeños en el diseño, errores en las estimulaciones, etc. Entre 5 – 20%
de ICF.
3.1. Modificaciones de proceso: Necesarias para cumplir con las especificaciones del
diseño.
3.2. Mano de obra de partida
3.3. Pérdidas en la producción: Pérdidas en las ganancias mientras se pone en operación
el proceso.
La inversión de capital es la cantidad total de dinero necesaria para suministrar los equipos
de la planta y de las instalaciones de manufactura más la cantidad de dinero requerido
como capital de trabajo para la operación de estas instalaciones.
La mayoría de los datos de costos disponibles para diseño preliminar son válidos para el
período en el cual fueron desarrollados. Debido a que los costos pueden cambiar en el
tiempo debido a condiciones económicas, se debe emplear algún método para actualizar
17
estos costos. Esto se efectúa por medio de los llamados índices de costos. Las correlaciones
de costos se comenzaron a publicar en Estados Unidos a comienzos de la mitad del 1900.
Las correlaciones más conocidas son las Chilton (1949), las de Guthrie (1968). Entre las
más recientes son las de Peters and Timmerhaus (2003) y las de ASPEN. Sin embargo
como se tarda entre 2 a 3 años en construir una planta química, y por lo tanto debemos
poder predecir los costos futuros de un proceso ya que los costos aumentan con el tiempo.
Materias Primas
Productos Terminados
Capital de trabajo Dineros en Efectivos
(CPM) Dineros Recibidos
Dinero a Pagar e Impuesto a Pagar
Existen varios métodos para actualizar costos, pero en general son similares en que debe
utilizar un costo base para un año dado y multiplicar este costo por la razón de índices de
los años considerados.
Uno de los índices más populares es el Marshall and Swift (M&S), el que es publicado en
la revista Chemical Engineering. La tabla siguiente entrega los valores más recientes del
Índice M&S. Otros índices son el Engineering News- Record Index, el Nelson Refinery
Index, el Chemical Engineering Plant Construction Index, y el Materials-and-Labor Cost
Index. Algunos de estos índices incluyen factores separados para mano de obra y
materiales, lo cual a menudo experimentan índices inflacionarios diferentes.
18
DESGLOSAMIENTO DE LA INVERISÓN DE CAPITAL FIJO
ICF = CD + CI
Edificios de procesos
Costos directos Edificios Edificios auxiliares
(CD) Talleres de mantenimiento
Edificios de servicios generales
Costos fuera de
Arreglo de áreas
límites de la planta
(CFLP) Servicios
Facilidades
Servicios de
Equipos no de proceso
planta
Distribución y
Inversión Terreno
de capital
fijo (ICF)
Imprevistos
19
Marshal and Swift Nelson- Chemical
Enginnering News Record
Installed –equipment Farrar Engineering
Construction Index
Año Index Index Index
Todas Industrias de
1913=100 1949=100 1967=100 1946=100 1959=100
industrias Procesos
1992 943.1 957.9 4985 1081 464 1277.3 358.2
1993 964.2 971.4 5210 1130 485 1310.8 359.2
1994 993.4 992.8 5408 1173 504 1349.7 368.1
1995 1027.5 1029.0 5471 1187 409 1392.1 381.1
1996 1039.1 1048.5 5620 1219 523 1418.9 381.7
1997 1056.8 1063.7 5825 1264 542 1449.2 386.5
1998 1061.9 1077.1 5920 1284 551 1477.6 389.5
1999 1068.3 1081.9 6060 1315 564 1497.2 390.6
2000 1089.0 1097.7 6221 1350 579 1542.7 394.1
2001 1093.9 1106.9 6342 1376 591 1579.7 394.3
2002 1102.5 1116.9 6490 1408 604 1599.2 390.4
La inversión total se evalúa en función de una capacidad base y de un tiempo base. Ambos
aspectos requieren de ajustes para tener la estimación de la inversión en el tiempo deseado
y para la capacidad requerida. Primero se considera el efecto del tiempo en la inversión y
luego el efecto de la capacidad.
m
Capacidad 2
Costo2 Costo1
Capacidad 1
Dónde:
I2 = Inversión en dólares de una planta de capacidad 2
I1 = Inversión en dólares de una planta de capacidad 1
m = Potencia (0,6)
C1 y C2 = Capacidades de producción que deben tener las mismas unidades
20
El efecto de la inflación en los costos de inversión de plantas químicas puede cuantificarse
mediante el uso de índices especializados. Para la información de la inversión de un
proceso de un año base a un año deseado simplemente se toma de los índices de los años
respectivos.
Indice _ del _ año2
I t 2 I t1
Indice _ del _ año1
La revista Chemical Engineering publica los índices para actualizar el costo de las plantas
químicas. El cálculo de este índice se basa en cuatro rubros con las siguientes
proporciones.
B. CORRELACIONES DE COSTOS
Existen variedad de procedimientos usados para estimar el costo de los procesos químicos.
Buscar en catálogos de tecnologías
Petrochemical Technology Assessment, Wiley, 1981
Hydrocarbon Processing y Chemical Engineering
Cotización de los procesos
Precios de materias primas
Chemical Marketing Reporte (Editorial Schnell Publishing Company)
El diseño para una nueva inversión en una planta industrial, también se estima mediante
los costos de los equipos más caros a instalarse en una planta industrial y a este tipo de
técnica se denomina Estudio Estimado.
La revista Chemical Engineering a partir de 1926, publica los índices y está basada en 8
tipos de industrias
Cemento 2%
Química 48%
Productos de Arcilla 2%
21
Vidrio 3%
Pintura 5%
Papel 10%
Petróleo 22%
Hile 8%
La revista “Oil and Gas Journal” publica estos costos y sus índices se basan en:
Para un proceso nuevo, publicaciones de correlaciones de costo no existe. Para este caso,
debemos diseñar la planta y el costo de cada equipo usando correlaciones de costo de
equipo.
Ejemplo:
M &S
Purchased _ Cost , US $
101,3 A
0 , 65
Fc
280
Dónde:
Fc = Factor de corrección para materiales, presión, etc.
22
Completa fabricación de tubo y carcasa (costo promedio año 1968)
Ejemplo:
Columna de destilación Costo de compra x 4
Reactores a presión Costo de compra x 4
Bombas e instrumentos Costo de compra x 4
Intercambiadores de calor Costo de compra x 3,5
Compresoras Costo de compra x 2,5
M &S
Installed _ Cost , US $
101,3 A
0 , 65
2,29 Fc
280
CORRELACIONES DE GUTHRIE’S
Costo instalado nos proporciona estimaciones separadas para el tubo requerido, concreto
para el soporte estructural, instrumentación convencional y controladores, instalación de
equipos eléctricos auxiliares necesarios, aislación y pintura.
Los factores de instalación por ejemplo para los intercambiadores de calor son de acero al
carbono, independientes de los factores de corrección para presión, materiales, etc.
23
Muchas compañías incluyen los costos de puesta en marcha como parte de la inversión de
capital. Otras compañías consideran la fracción de costos de puesta en marcha destinada a
las modificaciones de equipo como parte de la inversión de capital, mientras que los
dineros utilizados en mano de obra y materiales necesarios en la partida de la planta se
consideran como gastos operación. La decisión depende de la situación tributaria de la
compañía. Nosotros consideraremos la puesta en marcha como parte de la inversión. Por lo
tanto
CAP 0.15IT
CMP 0.11CF
El capital de trabajo, CAPT, representa los fondos necesarios para operar realmente la
planta, es decir, para pagar las materias primas, salarios, etc. Se intenta reponer el capital
de trabajo cada mes mediante los ingresos producidos. Sin embargo, se de disponer de
dinero antes de comenzar a operar para llenar los estanques y pagar los jornales, por esta
razón el capital de trabajo se considera parte de la inversión total. Una estimación
preliminar del capital de trabajo es tomar el costo de abastecimiento de materias primas por
tres meses. Sin embargo, el análisis inicial de inversión se puede simplificar en gran parte
si se supone que el capital de trabajo está relacionado a la inversión, y por esta razón
usaremos CAPT = 0,15IT.
De la definición de inversión de capital fijo tenemos que es la suma de los costos directos e
indirectos.
ICF CD CI
Los costos directos, incluyen los Costos Dentro de los Límites de la Planta (CDLP) y los
Costos Fuera del Límite de la Planta (CFLP).
CD CDLP CFLP
El costo CFLP, se refiere a la planta de vapor, torres de enfriamiento y otros que son
necesarios para la operación del proceso pero que se construyen en otra área geográfica. Es
una práctica común tener áreas centrales para las torres de enfriamiento, y los equipos de
generación de vapor. De la tabla se puede ver que los costos individuales del CFLP varían
en mayor proporción que los costos CDLP.
Usualmente se puede tener dos tipos de proyectos, el de una planta nueva (Grass Root) y el
de una ampliación. Los costos fuera del límite de la planta son muy diferentes en ambos
casos:
24
1. Plantas Nuevas : CFLP 200 – 400% CDLP
2. Ampliación : CFLP 40 -50% CDLP
Plantas Nuevas:
CFLP 3.00 CDLP
Ampliación de planta:
CFLP 0.45 CDLP
Si suponemos:
CP = 0.05CD y CON = 0.20CD
Entonces en
CI = 0.05 CD + 0.20 CD
CI = 0.25 CD
De modo que:
ICF = CDLP + CFLP + 0.25CD
Los factores que hemos seleccionado para nuestro análisis nos dan una estimación
razonable para los procesos petroquímicos. Nuestra meta es desarrollar un método simple
para el diseño preliminar de procesos, de modo que podemos usar otros factores para
simplificar el cálculo de la inversión de un proyecto para la inversión de un proyecto.
25
IT = 1.3ICF = 1.3 [1.25 (CDLP + CFLP)]
Los costos totales de producción pueden ser calculadas en una de las siguientes bases:
Diaria, por unidad de producto, o en base anual. La base anual es el mejor debido a que el
efecto de la variación por temporada se atenúa, permite un cálculo más rápido de los costos
de operación a menos de la capacidad total, etc.
La mejor fuente de información para utilizar en la estimación del costo total de producción,
son los datos de proceso similares o idénticos al estudiado. La mayoría de las compañías
tiene extensos archivos de sus operaciones, de modo que rápidamente se pueden estimar
los costos de producción. La siguiente tabla presenta un desglosamiento de los costos de
producción.
26
1. Costos Directos de Producción (CDP). Cerca del 60% del CTP
a) Materias Primas (entre 10 – 50 del CTP)
b) Servicios (entre 10 – 20% del CTP)
c) Mantenimiento y Reparación (entre 2 – 10% del ICF)
d) Suministros de Operación (entre 10 – 20% Costos de Mantenimiento y Reparación o
0,5 - 1,6% de ICF)
e) Mano de Obra (entre 10 – 20% del CTP)
f) Supervisión Directa (entre 10 – 25% de la mano de obra)
g) Gastos de Laboratorio (entre 10 – 20 de la mano de obra)
h) Patentes y Derechos (entre 10 – 20% del CTP)
1. Costos Administrativos: (Cerca del 15% del costo de mano de obra, supervisión y
mantenimiento, o 2- 3% CTP). Incluye costos de salarios ejecutivos, pagos legales,
accesorios, oficinas y comunicaciones.
2. Distribución y Venta: (2 – 20% del CTP)
Incluye costos de oficina de ventas, vendedores, avisos, etc.
3. Investigación y Desarrollo: (2 – 5 de unidad de venta o cerca del 5% del CTP).
CTP = CFA + GG
27
Los Gastos Generales (GG) corresponde aproximadamente a un 2.5% de los ingresos por
ventas, IV.
GG = 0,025 IV
El costo de fabricación CFA es la suma del costo directo de producción (CDP), los costos
fijos (CF) y los costos generales de fabricación (CGF).
CDP = CMAP + CS + 0,04 ICF + 0,006 ICF + 1.35 CMO + 0.03 CTP
Los costos fijos, CF, incluyen impuestos, IMP, seguros SEG, arriendos, ARR e intereses,
INT.
28
El interés sobre el capital prestado depende de la política financiera de la compañía, y para
un diseño preliminar vamos a suponer que se usan fondos internos para financiar la
operación, en este caso INT = 0; de forma similar, supondremos que no se arrienda el
terreno, entonces ARR = 0 y luego el:
CF = 0.03 ICF
Cuando combinamos todas estas expresiones obtenemos una expresión para el costo total
de producción.
CTP = [CMAP + CS + 0.046 ICF + 1.35 CMO + 0.03 CTP] + 0.03 ICF + [0.72 CMO +
0.024 ICF]+ 0.025IV
Ton US $
IV Pr oducción x Pr ecio _ Pr oducto
año Ton
29
El costo de mano de obra, CMO, depende principalmente de la complejidad del proceso
puede ser estimado a partir del diagrama de flujo, aunque para hacer una estimación
razonable se necesita tener experiencia. Para procesos pequeños (entre 2 y 4 cuadrillas de
trabajadores), siendo el NT el Número de Turnos requeridos, los costos de mano de obra
por turno son de alrededor de US$ 100000/año (considerando que se opera 24 h/día,
durante 7 días/semana y con alrededor de 4.5 operadores por turno), entonces:
CMO = 100000 NT
CMO = 105 NT
Otro intento de cuantificar la mano de obra es el trabajo de Wessel, Chem. Eng, 59, 7.209
(1952), es que correlacionó la mano de obra en (h-hombre)/(día) (etapas de proceso) versus
la capacidad de la planta (Fig. 5.8 Peters Timmerhaus). La dificultad con este
procedimiento está en estimar el número de etapas del proceso, ya que por ejemplo un
reactor discontinuo puede requerir un operario a tiempo completo, mientras un reactor
continuo puede requerir sólo medio tiempo de un operario.
Si consideramos las definiciones de CMO e ICF en CTP se obtiene para una ampliación de
planta:
Por lo tanto, podemos usar las estimaciones del costo de materias primas, servicios,
entradas y el costo instalado de los equipos de nuestro diseño preliminar para calcular el
costo total de producción en forma muy simple para una ampliación de una planta por
medio de:
La utilidad neta antes de impuestos, UNAI, es igual a los ingresos por ventas, IV, menos el
costo total de producción, CTP.
UNAI = IV – CTP
30
Para el caso de ampliaciones y si CMO = 0.15 CTP
[Link] DEPRECIACIÓN
Se define como “la cantidad de dinero que debemos descontar de nuestras utilidades en
orden de acumular suficientes fondos de modo de reemplazar los equipos o vehículos
cuando estos dejen de servir”. Estos fondos se pueden depositar en un banco y ganar un
cierto interés, pero usualmente se espera obtener una mayor ganancia invirtiendo este
dinero en algún otro proyecto de la compañía.
Los terrenos no sufren depreciación, por lo tanto, la inversión de terreno no debe ser
considerada en al cálculo de la depreciación. En forma similar, si reemplazamos el capital
de trabajo cada mes, tendremos la misma cantidad de capital de trabajo al finalizar el
proyecto, y por lo tanto no sufrirán depreciación. Más aun los equipos tendrán un cierto
“Valor Residual” al final de proyecto (a menudo un 10% del costo de compra, 3% ICF) de
modo que podemos contar con estos fondos al final de proyecto.
31
A pesar de los diferentes procedimientos para calcular la depreciación, en estimación
preliminar no se justificar este cálculo y supondremos que se puede estimar la depreciación
de la siguiente forma para el caso de una ampliación.
Por lo tanto, la depreciación debe ser descontada de la utilidad antes de impuestos, porque
representa un costo por reemplazo de equipo. Si queremos obtener la utilidad después de
impuestos, debemos descontar los impuestos, los cuales para el caso de grandes compañías
es alrededor del 48% ([Link].).
UDI = (1 - T) (UNAI - D)
Usando los valores anteriores para el caso de una ampliación de planta, se obtiene:
UDI = 0.52 [0.96 IV – 1.51 (CMAP + CS) – 0.27 CDLP - 0.181 CDLP]
UDI = 0.50 IV – 0.785 (CMAP + CS) – 0.236 CDLP
32
Ahora estamos en condiciones de calcular el flujo neto de caja (FC), que es la utilidad que
retiene la compañía, después de descontar los impuestos y más el costo de depreciación.
FC = UDI + D
FC = 0.52 (UNAI – D) + D
Las utilidades netas no pueden ser utilizadas como un factor de rentabilidad para decir si
una inversión debe efectuarse. La meta en utilidad de una compañía es maximizar las
entradas sobre el costo que el capital que debe invertirse para generar estas entradas. Si la
meta fuera simplemente maximizar utilidades, cualquier inversión sería aceptable, siempre
que entregara utilidades, no importando cuan bajo sea el retorno o cuán grande el costo.
Ejemplo:
Proyecto 1 Proyecto 2
Inversión 100000 1000000
Utilidad 10000 25000
10000 25000
Fórmula de la Tasa de Retorno TRI x100 TRI x100
1000000 1000000
TRI 10% 2.5%
Debido a que las inversiones normales entregan tasas de retorno del 10% la inversión I2 no
es muy atractiva sin embargo la tasa de un 10% sobre los US$ 100000 hacen de esta
inversión más atractiva que la segunda.
33
Además de la tasa de retorno sobre la inversión original, existen otros métodos de
rentabilidad, tales como el periodo de pago, los métodos de flujo de caja descontando
(TIR, VAN), etc.
UDI
TRI x100
IT
La tasa de retomo puede estar basada, en la utilidad antes o después de impuestos, además
debemos tener presente que el capital de trabajo, así como la porción del costo de puesta en
marcha que es considerado como una inversión para efectos de impuestos deben ser
incluidos en la inversión total.
La tasa de retorno, es una medida muy simple de la rentabilidad, pero no considera el valor
del dinero en el tiempo, de modo que no es realista, además debe estar basado en algún
tipo de operación promedio anual, puesto que no se puede utilizar factores de
mantenimiento que aumentan durante el periodo de vida de la planta. A pesar de lo anterior
la tasa de retorno es utilizada frecuentemente en diseño preliminar.
Otra medida de la rentabilidad que es el periodo de pago, que es el tiempo que se toma
para recuperar los fondos que hemos invertido. Debido a que el capital de trabajo se
recupera cada mes, por lo tanto, no se toma en cuenta en este cálculo. Sin embargo, la
fracción del costo de puesta en marcha que se considera en la inversión debe agregarse a la
inversión de capital fijo en orden para encontrar la cantidad de dinero asociado al proyecto.
Los fondos que se recuperan del proyecto son la utilidad después de impuestos más la
depreciación, lo mal hemos llamado Flujo de Caja, de modo que el periodo de pago es:
ICF CPM
PP(años)
FC
Este criterio también es muy simple de calcular, pero sufre el mismo problema que la tasa
de retorno por lo tanto desearíamos contrastar estos procedimientos simples con un análisis
más riguroso que toma en cuenta el valor del dinero en el tiempo. Entre los métodos más
complejos se encuentran los métodos de Flujo de Caja Descontado, y de los más utilizados
son el Valor Actual Neto (VAN) y la Tasa Interna de retorno (TIR).
“Se puede definir como la más alta tasa de interés que se podría pagar por un préstamo que
financia la inversión, y que dejara indiferente su aceptación, si dicho préstamo se fuera
34
abonado con los fondos provenientes de la inversión a medida que se fuera generando”.
También se define como aquella tasa de descuento que iguala el valor presente actual de
los flujos de caja neto positivos con el valor presente de los flujos netos negativos de caja
(inversión) generados por el proyecto de inversión a través del tiempo.
(+)
Período Pago
VR
FC 0
Costo terreno
ICF
CPM
(-)
0 t, anos n
Puesto que toma 3 años en construirse una planta, los costos directos, CD debe colocarse
sobre este período de tiempo. Al comienzo, se debe pagar por el terreno, contratar un
contratista con su personal, pedir los equipos, preparar el terreno. Y comenzar hacer la
fijación para los equipos. Después de esto se comienza la instalación de los equipos a
medida que se entregan. Así, en un período de tiempo menor a 3 años el costo directo
gastado es alrededor de 10 a 15% del total. Durante el período menor a 2 años y menor a
un año se gasta aproximadamente de un 40 a 35 % del costo directo cada año y
normalmente el último año se gasta el restante 10 a 15%. En cuanto a los gastos de
ingeniería y supervisión, imprevistos y salarios se pueden gastar en forma uniforme
durante el periodo de construcción.
El capital de trabajo y los costos de puesta en marcha se invierten en el tiempo cero, pero
hay que recordar que el capital de trabajo se recupera al final del proyecto. De igual
manera el valor residual de los equipos, VR, se puede hacer efectivo al final del proyecto y
es alrededor del 10% del costo de compra de los equipos, o aproximadamente un 3% de la
inversión de capital fijo.
La mayoría de las plantas nuevas no producen a plena capacidad en los primeros años de
operación y se puede considerar que la producción aumenta desde 60% a 90 % y 95%
durante el primer, segundo y tercer año de operación. Después de este tiempo se supone
que la planta opera a plena capacidad.
En forma similar, la depreciación variará cada año en forma lineal, a menos que se utilice
con ventajas variables y depreciación variable. En este caso las variables de las utilidades
anuales, los impuestos y las utilidades netas variarán año a año.
35
El flujo de caja o utilidad anual neta, es la suma de la utilidad después de impuestos más la
depreciación, y se denomina así, puesto que es la cantidad de dinero que retiene la
compañía cada año. Un programa más claro se obtiene al graficar lo anterior de la siguiente
forma (figura 2).
Con lo visto anteriormente, se puede efectuar un análisis para el cálculo de la tasa interna
de retorno, y esto se obtiene, igualando el valor actual o presente de la inversión, y el valor
presente de los flujos de caja, con las siguientes suposiciones.
2. Los ingresos por ventas son constantes, excepto los tres primeros años.
CAPT
100% 100% VR
95%
90%
60% FC
-3 -2 -1 0
1/4CI
1/4CI
1/4CI 1/4CI
CAPT
CPM
D1 = D2 ... DN
36
b) Valor. Presente de los Flujo de Caja.
El análisis se vuelve muchos más simple, si las inversiones, los flujos de caja y la
depreciación son uniformes. Con flujo de caja constante, podemos usar las fórmulas de
interés y anualidades para sumar las series anteriores.
ICF 1 i 1 1 1 i N
4
4 I i
ICF
4
1 i 1 CAPT CPM i 1 i 0.52IV CTP 0.48D 1 i 1 CAPT VRi
N N
4
CAPT + VR
FC
FC
FC
-3 -2 -1 0
ICF/4 1 2 3 n
CAPT + CPM
Dónde: La inversión total (IT), incluye el capital de trabajo y los costos de puesta en
marcha. Lo anterior es una forma de relacionar las utilidades antes de impuestos con la
inversión de capital. Usando esta definición, podemos obtener una expresión de i, en la que
hay solo inversiones.
ICF
4
1 i 1 CAPT CPM i 1 i 0.52FRC IT 0.48D 1 i 1 CAPT VRi
N N
4
ICF
4
1 i 1 CAPT CPM i 0.48D
FRC
4 CAPT VR i 0.48D
0.52 IT 1 i 1 1 i
N N
0.52 IT 1 i 1
N
Si introducimos las siguientes simplificaciones.
VR = 0.03 ICF
IT = 1.30 ICF
D = 0.10 ICF
FRC
0.251 i 4
0.29i 0.298 1 i 0.225i 0.048
N
0676 1 i 1
N
Basados en simples suposiciones, podemos ver que existe una relación directa entre los
siguientes valores:
1
FRC
3
38
2.4.6 MODELO SIMPLIFICADO DE LA RENTABILIDAD:
Por lo tanto, de los valores estimados de los ingresos por ventas, IV, el costo de las
materias primas, CMP, de los servicios, CS, y del costo de los equipos instalados, CDLP,
se puede calcular el factor de recuperación de capital, FRC, correspondiente al diseño. Si
este valor resulta mayor que 0.333 el proyecto es atractivo y se justifica, si el valor de
FRC, es menor que 0.333 se hace FRC = 1/3 en la expresión anterior y se calcula el ingreso
por ventas.
A partir de estos ingresos, se puede calcular el precio del producto requerido para que el
proceso sea rentable, y se efectúa entonces un análisis de oferta y demanda. Si designamos
por (Pp) al precio del producto anual [kg/ año], se tiene que:
IV = PP (P)
39
IV = 1.57 (CMAP + CS) + 1.10 CDLP
1.57CMAP CS 1.10CDLP
PP
P P
Debemos recordar que nuestra primera estimación del potencial económico, basado en los
precios actuales es:
PE1 = IV – CMAP
Después de que se completa los balances de materia y energía, podemos evaluar los
suministros y podemos revisar el potencial económico.
PE2 = IV – CMAP – CS
n
PEn IV CMAP CS CDLP j
j 1
Si este potencial económico es mucho menor a cero en cualquier punto, debemos generar
otras alternativas del proceso, menos costosa, en vez de continuar con el diseño.
Si el análisis del precio del producto nos indica que: Pp es mayor que el precio del
mercado, entonces se debe hacer un estudio de oferta y demanda del mercado, o cambiar
de alternativa, o simplemente posponer la idea. Caso contrario, si el Cp es menor que el
precio del mercado, entonces se debe construir la planta aun cuando el precio sea mayor,
para obtener mayor ganancia.
IV CTP
Volviendo a la definición de FRC. Si el FRC 1
IT
Este factor puede entenderse como la fracción de la inversión total IT que necesitamos
recuperar cada año y para obtener una cierta tasa interna de retorno. Es decir, si nos
imponemos un costo cada año de FRC (IT), el proceso nos dará una tasa de retorno
determinada por FRC.
40
CAPITULO III: REGLAS HEURÍSTICAS Y MÉTODOS
APROXIMADOS
Para los diseños preliminares solo deseamos utilizar los parámetros más importantes de un
proceso y suprimir todos los detalles que tienen efectos muy pequeños en la economía del
proceso. Para lograr el propósito anterior, a menudo utilizamos argumentos de orden de
magnitud para simplificar las ecuaciones de diseño.
Es interesante notar que las soluciones de los problemas de diseño óptimo de muchos
sistemas son esencialmente independientes de los parámetros y de los coeficientes de
costo.
Área lateral Al D H
.D 2
Área de la base A
4
41
w x t x AT (1)
Dónde:
w = Peso
= densidad del material (plancha)
t = Espesor
AT = Área
2. .D 2
w x t. .D.H
4
.D 2
w . t .D.H (4)
2
.D 2 .H
v v .r 2 .H (5)
4
Despejando la altura: H:
4.v
H (6)
.D 2
4v .D 2 4v .D 2
w .t. .D. w . t. (7)
.D 2 2 D 2
42
Entonces, el diámetro del estanque que minimiza el peso será:
w 4v 4v 4.v
0 2 .D 2 .D D 3 (8)
D D D
4 .D 2 .H
D3 D 3 D 2 .H DH (9)
4
H
Regla Heurística: Para estanques atmosféricos usar: 1
D
Este resultado es bastante razonable para estanques de almacenamiento de petróleo en
refinerías. Las dimensiones óptimas son independientes del volumen del material.
Para almacenar productos a alta presión como por ejemplo gas licuado de petróleo.
a. Estanque compresora
El espesor de la pared debe tomar en cuenta la presión para prevenir la rotura del estanque
(soldadura).
43
Figura 3.3. Balance de fuerzas en un estanque a presión
P.D.L P.D
t t (2)
2.S .E.L 2.S .E
P.R
t (3)
.S .E 0,6.P
Dónde: R = Radio
Empleando un factor de corrección por corrosión de al menos de pulgada, de modo
1 "
8
que la ecuación anterior será:
P.R
t t0 (4)
.S .E 0,6.P
D
Si R , entonces la ecuación (4) llega a ser:
2
1.P.D
t 2
t0 (5)
.S .E 0,6.P
C S US$
Kg
Los cabezales son más gruesos que de los costados, es por ello su costo es el doble de los
costados.
44
C E Costo de la plancha * Densidad (solido ) * Espesor * Área
C E C S . S .t. A (6)
A Área lateral Área semiesferica
A Al ASE (7)
Dónde:
2 .D.L
Al 2. .r.L Al Al .D.L (8)
2
Puesto que la presión que actúa sobre los cabezales de un estanque pueden hacer que se
flecten hacia afuera, los estanques a presión se construyen normalmente con cabezales
elipsoidales 2:1. Como el área de un círculo es .R 2 , mientras que el área de una
hemisfera es 2 .R 2 , entonces podemos suponer que el área superficial de un cabezal
elipsoidal será:
D 2
C C S . S .t. .D.L 21,52
4
C C S . S .t. .D.L 1,5.D
2
(11)
.D.L 1,5 .D
1 P.D
C C S . S .t 0 2 2
(12)
SE 0,6 P
b.2.1. Optimizando el costo con respecto al diámetro para hallar las dimensiones
óptimas del estanque:
El costo depende más del área, que el espesor, por lo tanto t constante. Despreciando el
volumen de los cabezales, y considerando solo el volumen del estanque
45
.D 2 .L 4v
v L (13)
4 .D 2
C C S . S .t. .D.
4v
4v
1,5 D 2 C C S . S .t. 1,5 D 2 (14)
.D 2
.D
Optimizando:
C 4V
D
0
D D
4V
1,5 .D 2 2 3. .D
D
4V 4V
3. .D 4.V 3. .D 3 D 3 (15)
D 2
3
4 D 2 L
D
3
.L 3D 3 D 2 .L 3D L 3 (16)
3 4 D
L
Regla Heurística: Para estanques a presión usar: 3
D
b.2.1.2. Optimizando el costo con respecto al espesor:
De ecuación. (12)
.D.L 1,5 .D 2
1 P.D
C C S . S .t 0 2
SE 0,6.P
L P.D
31
D 4.t 0 SE 0,6.P
L
En general una Regla Heurística es: 2 5
D
Para comprimir n (moles/h) de un gas a una temperatura Ti desde una presión inicial Pi
hasta una presión final Pf , podemos usar una compresora de una etapa o una unidad
multicapas.
46
k 1
MRT1 K f k
P
hp w 1 (1)
550 K 1 Pi
k 1
K Pf k
hp 0,00436QPL 1
K 1 Pi
Dónde:
hp = Potencia
w = Trabajo
M = Masa gas o peso del gas lb / s .
T1 = Temperatura ingreso
Pi = presión ingreso lb 2
pu lg
3
Q = Velocidad flujo ft
min
K = Constante (exponente de gas adiabático)
K 1
K 1
K 1
K P2 K P3 K P4 K
w MRT1 3 (2)
K 1 P1 P2 P3
Para un compresor de tres etapas podemos encontrar las presiones intermedias que
minimizan el trabajo total w .
1 K 1 K 1 1 2 K
w
0 P1 K
P K
P K
P K
(3)
P2
2 3 2
1 K 1 K 1 1 2 K
w
0 P2 K
P3 K P4 K
P3 K
(4)
P3
K 1 K 1
w P2 K P3 K
0 K 1
P2 P2 KK1 P2
P P K
1 2
47
K 1 K 1
w K 1 K K 1
1
1 1
0 P2 P3
K
P K
P2 K 1 2
K
K K
P1
w K K K 1
1 K 1 K 1 1 2 K
P1 P2 P3 K P2 K 0
P2 K
1 K 1 K 1 1 2 K
P 1
K
P 2
K
P 3
K
P 2
K
0 (5)
Luego:
1 K 1 K 1 2 K 1 1 K 1 K 1 2 K 1
P 1
K
P
3
K
P 2
K
P2
K
P 1
K
P 3
K
P2 K K
1 K 1 K 2 2 K 1 K 1 K 1 K
2
P1
K
P
3
K
P 2
K
P 1
K
P 3
K
P 2
K
P22 P1 P3 (6)
1 K 1 K 1 1 2 K
P 2
K
P 3
K
P 4
K
P 3
K
0
P32 P12 .P4 1
3
.P4 P3 P1 .P42 1
3
(10)
1
P P P Pfinal N
RS 2 3 4 RS Rt N RS
1
P1 P2 P3 pinicial
REGLA HEURISTICA
“Las razones de compresión de cada etapa en un compresor deben ser iguales”
48
3.2.1 MÉTODO APROXIMADO DE DISEÑO DE GASEODUCTOS
F2
t
F1
Di
Di D0
a) Balance de la fuerza F1 : Esfuerzo sobre la pared del tubo o caño, propio de la presión
de diseño.
De la figura: D0 Di 2t (3.3)
2.S .t
.D0 .[Link]ño 2 .L.S .t Pdiseño (3.5)
D0
Dónde:
Pdiseño = Presión de diseño del caño, Psia
S = Mínima tensión fluencia (Es la máxima presión que el acero puede sostener
antes de perder su elasticidad), Psia.
t = espesor del conducto, Pulg.
49
D0 = diámetro externo de la cañería, Pulg.
2.S .t
Pdiseño .F .L.J .T (3.6)
D0
Dónde:
F = Factor de diseño
L = Factor de localización de la construcción
J = Factor de juntas o soldaduras tipo normal para la línea de transmisión
T = Factor de corrección de temperaturas
2.S .t
Pdiseño .F .E.T (3.7)
D0
Dónde:
FxL F Factor de diseño
E = Factor de soldadura longitudinal
T = Factor temperatura
2.2. Factor de localización: L=1; porque las especificaciones de los aceros según estándar
API es 5LX.
Tipo cañería E
Sin costura 1
Soldado eléctrico 1
Soldado con arco sumergido 1
Soldado al horno con mango 0,6
2.4. Factor de temperatura: T=1; para temperaturas del gas 120 0 C , si la temperatura es
superior se obtiene de normas API el factor correspondiente.
50
2.5. Factor de tensión fluencia: S; de acuerdo al tipo de caño según API estándares para
5LX.
API 5LX grado x65 65000
API 5LX grado x80 80000
En una línea de transmisión de gas entre el 50-55% del costo del proyecto es el peso del
acero usado en la construcción del gaseoducto. Entonces, es importante estimar la cantidad
de acero requerido.
ws
4
D 2
0
Di2 ..L. ws
4
.L. . D0 Di .D0 Di (3.10)
2t
c) Espesor del caño:
D0 Di 2.t D0 2.t Di (3.11)
Reemplazando (3.11) en (3.10)
ws .L. .2t.2t Di Di .L. .2t.2t 2 Di .L. .4t.t Di
4 4 4
ws .L..t.Di t (3.12)
Dónde:
L=m
t =m
ton
m3
D=m
51
En una línea de transmisión de gas el 50-55% del costo del proyecto es el peso del acero
usado en la construcción del gaseoducto. Entonces, es importante estimar la cantidad de
acero requerido.
wd Vs . (3.14)
De ecuación (3.15)
VS D0 Di D0 Di .L (3.18)
4
VS 2t Di 2t Di .L .2t .2 Di 2t .L
4 4
.4t.L.Di t
Luego: VS (3.19)
4
wd .L..t.Di t
Cducto .L..t.Di t .Cacero
El flujo del gas y líquido en un absorbedor deberán cambiar cuando el soluto es transferido
desde el gas a la corriente líquida.
1) Balance total:
G = Flujo total de gas.
L = Flujo total de líquido.
52
Gs G1 y Flujo total de gas de salida.
Ls L1 x Flujo total de líquido de salida.
G, Ys Le, Xe
N-1
G, Ye Ls, Xs
Gs 1 y Ls 1 x (3.20)
y x
Y ;X
1 y 1 x (3.21)
Para presiones moderadas PT (50 – 100 Psig) la mezcla de gases puede considerarse ideal
y podemos escribir como:
PY P X (3.23)
53
P = Presión de vapor (referencia)
X = Fracción molar.
= coeficiente de actividad.
De ecuación (3.21):
Y
y (3.24)
Y 1
X
x (3.25)
X 1
De ecuación (3.23):
P
y x (3.26)
P
Y P X
y (3.27)
Y 1 P X 1
P X
y
P X 1
P
Si hacemos m = la cual es llamada la Ley de Henry; entonces:
P
y n mxn (3.28)
Para mezclas diluidas el problema es simple y podemos obtener una solución analítica de
las ecuaciones de diseño. El resultado es llamado “ecuación de Kremser” y puede ser
escrito como:
ye mxe A N 1 1
(3.29)
y s mxe A 1
Dónde:
54
L
A (3.30)
mG
y mxe
Ln1 A 1 e
y s mxe
N 1 (3.31)
LnA
ye y s A N 1 A
N 1 . (3.32)
y e mxe A 1
ys
f 1 (3.33)
ye
Entonces:
ye y s
f (3.34)
ye
O también:
ys
1 f (3.35)
ye
y e y s A N 1 A y y s A N 1 A
N 1 e N 1 (3.36)
ye 0 A 1 ye A 1
A N 1 A
f (3.37)
A N 1 1
Dónde:
A = Factor de absorción.
f = Recuperación fraccional.
55
“Cualquiera de estas variables pueden determinarse a partir de la figura si se conocen las
otras dos”.
1.0
A N 1 A
f
os
A N 1 1
N=4
lat
∞p
0.5 N=2
N=1
0.0
1.0 L 2.0
A
mG
Figura 3.6. Gráfica de ecuación de Kremser
a) Si A <1 muchas etapas (no se puede obtener una recuperación completa del soluto)
b) Si A >2 los flujos de líquido aumentan y las columnas son bastantes costosas
c) Si f 0 entonces y s y e , y según ecuación (3.33) tendríamos que la columna es
pequeña y existiría una gran pérdida de soluto.
d) Si f 1 entonces y s 0 , con ello la columna sería grande y la pérdida del soluto
pequeña.
Conclusión:
El diseño óptimo de absorbedores se encuentra en el rango de 1 A 2 , y el valor óptimo
usado es con: A 1.4
De ecuación (3.31):
y mxe
Ln1 A 1 e
y s mxe
N 1 (3.38)
LnA
Hacemos:
56
y mxe
B 1 A 1 e (3.39)
y s mxe
ys y y
1 f s 1 0.99 0.01 e 100 (3.41)
ye ye ys
y
B 1 1.4 1 e B 1 1.4 1100 1 40 (3.42)
ys
y mxe
1 A 1 e 1 40
y s mxe
y mxe ye
Ln1 A 1 e Ln A 1 (3.43)
y s mxe ys
y
Ln A 1 e
ys
N 1 . (3.44)
LnA
57
Por lo tanto: LnA A 1 porque A 1 1.4 1 0.4 (3.46)
y
Ln A 1 e
ys
N 1 (3.47)
A 1
De la figura de Kremser N+1 tiene un alto número de etapas y estas resultan muy caras por
eso hacemos que: N 1 N y la ecuación (3.47) llega a ser:
y
Ln A 1 e
ys
N (3.48)
A 1
y
2.3 log A 1 e
ys
N (3.49)
A 1
y y
2.3 log 0.4 e 2.3 log e
ys 2.3 log 0.4 ys
N N
0.4 0.4 0.4
y
N 2 6 log e
ys
ye
N 2 6 log (3.50)
ys
La ecuación (3.50) es llamada la ecuación simplificada de Douglas
y
Ejemplo: Para e 1000
ys
58
N 3.72 2.97Ln1 f
ye y y
1000 s 0.001; f 1 s f 1 0.01 0.999
ys ye ye
G, Ys Le, Xe
US$/año N US$/año
Pérdida soluto Costo solvente
N-1
US$/año
Costo columna 2
G, Ye Ls, Xs
59
y e y s A N 1 A A N 1 A
N 1 f N 1 (3.51)
ye mxe A 1 A 1
Dónde:
L y
A ; f 1 s (3.52)
mG ye
Compromiso Económico:
Si f 0 columna pequeña, entonces gran pérdida de soluto.
Si f 1 columna grande, pequeña perdida de soluto.
Podemos utilizar el modelo de costo total analizado (CTA), para obtener la fracción de
recuperación óptima:
CTA=FRC x IT +CTP
De otra forma:
CTA=FRC x CAP +COP (3.53)
Dónde:
CTP = Costo total de producción
CAP = Costo capital
COP = Costo de operación
Considerando que se utiliza agua como solvente y su costo es bastante bajo, casi
despreciable, entonces tendremos que:
CAP [Link] [Link] US $
año
N
CAP C p AC (3.54)
E0
Dónde:
N = Número de platos.
E 0 = Eficiencia (%).
C p = Costo de plato US $ 3
plato pie
AC = Área de la columna ( pie 3 )
60
COP [Link]ón Costo. pé[Link] US $ año
COP H C y G ye 1 f s
Y
(3.55)
y y
f 1 s s 1 f
ye ye
Dónde:
H = Horas de trabajo al año.
C y = Costo del soluto
G = Flujo de gas (aire o vapor).
1
año 1
FRC (3.56)
3
Reemplazando (3.54), (3.55) y (3.56) en (3.53):
1 N
CTA C p AC H C y G ye 1 f (3.57)
3 E0
G MG
AC (3.58)
1.53600 m
Dónde:
M G = Peso molecular promedio del gas (aire) Lb/Lb mol.
m = Densidad molar del gas (aire) (Lb mol / pie3)
1.5 = Factor para el área de la columna que proviene de suponer un factor de
inundación de 60 %.
MG
1 3.72 2.97 Ln1 f
G
m
CTA Cp H C y G y e 1 f (3.60)
3 E0 1.53600
61
C p 90US $ 2 (3.61)
plato pie
H 8160
h
340dias [Link].8400350días
año
MG
1 3.72 2.97 Ln1 f
G
m
CTA 90 8150 C y G y e 1 f (3.62)
3 0.4 1.53600
MG
90
m
CTA
3.72 2.97 Ln1 f 8150 C y y e 1 f (3.63)
G 3 1.5 3600 0.4
MG
3.72 2.97 Ln1 f
90
CTA G m
8150 C y y e (3.64)
1 f 3 1.5 3600 0.4 1 f
Arreglando:
MG
3.72 2.97 Ln1 f
90
CTA G m
8150 C y y e
1 f 3 1.5 3600 0.4 1 f
MG
CTA 90
G 8150 C y m 2.97
0
1 f 3 1.5 3600 0.4 1 f
y e
Despejando f logramos:
m 2.97
MG
90
f opt 1 (3.65)
3 1.5 3600 0.4 8150 C y y e
ye 0.015
T y 77 F
C y 9.86US $
mol
18US$ Lb
62
M G 29aire
m 0.0802
90 29
0.0802
2.97
f opt 1 0.999
3 1.5 3600 0.4 8150 9.86 0.015
El costo total de una columna de absorción de platos se puede aproximar con la siguiente
ecuación:
Una estimación aproximada del costo de capital de una columna según la correlación de
Guthrie (costo instalado de una columna de acero al carbono incluyendo un 5% del costo
de compra e instalación de platos) está dado por:
H 0.82 3.30 FP
M &S
CAP K C DC
1.066
(3.69)
280
Dónde:
63
FP 1, cuando, P Pmin
a 3.91 10 3 Psia 2 5.68 10 4 Kpa 1
FP 1.0 ...aP Pmin , cuando, P Pmin Pmin 65Psia450Kpa
K C 107US$
pie1.886
1006US$m 1.886
Altura H f N platos, N
El espaciado de platos en una columna normalmente es de 2.0 ft (0.60 m), y podemos
suponer que se necesitan el 15 % de exceso en la altura y del fondo para que se desprenda
el vapor en el tope y el inventario del líquido en el fondo. De este modo tendremos que:
N
H 2 .3 (3.72)
E0
Existe una variedad de métodos para el cálculo de eficiencia. Una de las más populares es
la correlación de O’connell. La eficiencia global de platos está dada en función del
producto de la viscosidad y volatilidad de alimentación, y si aproximamos una línea recta a
la correlación podemos usar:
E0
0.377
0.209
30 60%[Link]áfico (3.73)
m M L L
L
Dónde:
E0 0.550% (3.74)
Diámetro D f [Link]
64
El flujo molar de vapor o aire = G mol
hora
(3.75)
3
Si el flujo volumétrico de vapor o aire es = Q ft , entonces el área lo hallamos de:
h
Q Q
Q AT V V .ó..... AT (3.76)
AT V
Q
G mol
h
ft 3
(3.77)
m mol 3
h
ft
G
m
0.88 AT (3.78)
V
El diámetro de una columna se selecciona de modo que la velocidad de vapor esté entre 60
y 80 % de la velocidad de inundación. Suponiendo que la inundación es causada por el
arrastre de gotas llevadas por el vapor.
FD E
Gota
mg
Balance de fuerzas sobre una gota de líquido (suponiendo que el arrastre comienza cuando
la gota está estacionaria en la corriente de vapor).
1
R 3 L G g C D G V 2 R 2
4
(3.79)
3 2
65
8 g R L G
R 3 g L G 0.44 G V 2 R 2
4 1
V 2
3 2 3 0.44 G
8 g R L G
V (3.80)
30.44 G
Podemos hacer:
8 g R
K0 (3.81)
30.44
L G
V K0 (3.82)
G
Para sistemas que operan cerca del punto crítico ( G L ), por lo que la ecuación (3.82)
llega a ser:
L L
V K0 V K0 V G K 0 L (3.83)
G L
Puesto que la densidad del líquido es aproximadamente constante para muchos sistemas,
con lo cual la ecuación (3.83) es:
Por lo tanto:
1
pie Lbm
2
F 1.5
3
s pie
pie Lbm
2
F 5400
3
(3.85)
h pie
1
m Kg
2
F 6590 3 (3.86)
h m
1
L
2
66
1
L
2
que a menudo se encuentra en el rango de 0.01 y 0.2 (no varía demasiado) G para
G L
''
platos espaciados entre 2 pies ( 24 ,0.6m ) por lo que una estimación de la velocidad de
inundación será promedio del rango (asumiendo L 58 Lb 3 y 20 dinas ).
ft cm
De ecuación (3.84)
V F G 2
F 1
F V V (3.87)
G
G
M G m (3.88)
G
P
.. [Link] (3.89)
m R T
P
MG (3.90)
G R T
A.4) Cálculo del diámetro para el área de una columna circular se tiene:
DC 2
AT (3.91)
4
V F M G (3.92)
R T
DC2 G m
0.88 AT 0.88
4 1 (3.93)
P 2
F MG
RT
67
1
1 2
4 M R T 2 1
DC G
G 2
(3.94)
0.88 F P
1) Primer caso:
1
F 5400 R 10.73
h pie3
y
Lbmol R
M T 14 1
DC 0.0296 G G 2 (3.95)
P
Dónde:
DC = Pie
T=R
P = Psia
G = Lbmol/h
2) Segundo caso:
1
m Kg 2
KJ
F 6590 3 y R 8.314
h m Kgmol K
M T 1
4 1
DC 0.0252 G G 2 (3.96)
P
Dónde:
DC = m
T=K
P = KPa
G = Kgmol/h
0.533
1
M &S M G R T 2
CAP KC 4
G 121.066 1.15 h N 0.82 3.3 FP (3.97)
280 0.88 F P
68
0.533
1
M &S 2
CAP KC
4
M G R T
1.15 h0.82 3.3 FP G 0.533 N 0.82 (3.98)
280 0.88 F P
Si hacemos:
0.533
1
M &S MG R T 2
K1 C
K
4
1.15 h0.82 3.3 FP (3.99)
280 0.88 F P
CS H y C y M G G ye 1 f (3.101)
Dónde:
Hy = Horas por año de trabajo (8400 h/año)
Cy = Costo del soluto por kg
CW H y CW MW L (3.102)
Dónde:
Cw = Costo del solvente por kg
Si:
L
A L A m G (3.103)
mG
CW H y CW MW A m G (3.104)
Si ahora hacemos que:
K2 H y C y M G (3.105)
K3 H y CW MW (3.106)
CS K 2 G ye 1 f (3.107)
CW K 3 A m G (3.108)
69
Finalmente reemplazamos las ecuaciones (3.108, (3.107) y (3.100) en la ecuación (3.68) se
obtiene la ecuación para determinar el costo total de un absorbedor:
K1 G 0.533 N 0.82 K 2 G ye 1 f K 3 m G A
1
CTA (3.109)
3
Estimar separación
componentes no claves 2
Determinar condición de
alimentación 4 Flash adiabático
71
En una columna, el factor de operación es la PRESION que tiene efectos sobre la
separación de los componentes. Entonces en el diseño se debe tener en cuenta tres
aspectos:
a. Presión de diseño elegido.
b. Tipo de enfriamiento o refrigerante a usar en el condensador.
c. Estabilidad térmica de los componentes de la mezcla a separar.
C. Estabilidad térmica
Si los componentes llegan a su Tc pueden: Descomponerse, polimerizarse, condensarse o
interaccionar entre ellos. Para esto se debe tener en cuenta que la T del rehervidor que es el
valor más alto debe cumplir este requisito y luego seleccionar la presión de operación,
incluso puede ser negativo (vacio).
CRITERIO GENERAL:
1. Escoger la presión de operación entre:
0 – 2860 KPa (28 atm, 45 Psia) Estos dependen de
2. Elegir una temperatura mínima factores económicos
50ºC (material, espesor, etc.)
3. Elegir el refrigerante
Agua
4. Estimar un ΔP de carga de la columna y condensador.
35 KPa (0.34 atm, 5 Psia)
72
ΔP = N* 0.1 Psi/Etapa
- Para columna que operan a P de vacio
ΔP = N* 0.05 Psi/Etapa
5. Tipo de condensador
d1
b1
Escoger Refrigerante
para operar Bajar Presión
condensador parcial a PD
PD =28 atm (420 psia)
- Para presiones que se encuentren 0 – 1480 KPa (15 atm, 215 Psia) un condensador
total (Refrigerante agua)
- Para presiones entre 1480 – 2520 KPa (25 atm, 365 Psia) un condensador parcial
(Refrigerante agua)
- Para presión ≥ 2520 KPa utiliza otro refrigerante.
B Pi 5 A2
log P A o ln A1
C t Pc T A3
73
2.- Ecuación Clapeyron (binario)
H vo 1 1
ln P ln P0
R T0 TD
Dónde:
P 0 Punto de ebullicion normal (1atm)
TD Se aproxima a 50ºC
1751
ln P 6,777
TbF
Dónde:
TbF Punto de ebullición de la alimentación
P = MPa
T=K
4.-Ecuación de Clausius-Clapeyron
N H v ,i ,o 1 1
PD P0 . X i , D exp
i R Toi TD
Dónde:
P0 = Presión de referencia
T0i = Punto ebullición de componente y a P0
TD = Punto de burbuja del destilado a TD
H v ,i ,o = Calor latente de vaporización
Pci T
Ki . exp 5,421 ci Ecuación de Mollepop
P T
D) Punto de Burbuja
C
1 X i .K i .P, T
i 1
74
1 Z
Yi .X i i Sistema binario
Z1 Z2
=función de evaporación
1 K 2 K1 1
Isotérmico, adiabático
Z i K i 1
FO
i K i 1 1
Adiabático
Hˆ f Hˆ C
Hˆ V Hˆ L
Curva de equilibrio
q
X A
R min 1 Línea de operación sección
tope
Línea de alimentación
R X L
Línea de operación sección tope: Y X D R
R 1 R 1 D
q XF
Línea alimentación: Y
q 1 q 1
L B
Línea operación sección fondo: Y X XB
L B L B
75
2
Simbolizando:
1
xB xD
Figura 3.10. Cálculo de número de etapas
K CL
aCL ,CP (3.110)
K CP
X CL , N 1 X CP ,1
S (3.111)
X CL ,1 X CP , N 1
Si consideramos que:
X CL , N 1 X CP ,1
1 Y 1 S 1 (3.112)
X X
CL ,1 CP , N 1
X CL , N 1 X CP ,1
PN min
Volatilidad relativa
(3.113)
X CL ,1 X CP , N 1
X CL , N 1 X CP ,1
log
X
CL ,1 CP , N 1
X
log CL ,CP P
N min Ecuación de Fenske (3.114)
76
AB tope PA0 TD
0
PB TD
PB TB
Entonces:
CL ,CP P . CL CL
(3.116)
CP T CP B
P T . B (3.117)
f CL FZ CL , F d CL DX CL , D bCL BX CL , B (3.118)
f CP FZ CP , F d CP DX CP , D bCP BX CP , B
d b
log CL CP
N min bCL d CP
Ecuación de Fenske (3.119)
log CL
CP P
En los cálculos debemos tener en cuenta que el número mínimo de etapas depende del
grado de separación de los componentes claves y de su volatilidad relativa; sin descuidar
que estas variables son independientes de las condiciones de alimentación.
Las ecuaciones (3.114) y (3.119) son válidas en forma exacta para dos etapas teóricas, y
para una etapa corresponde a la ecuación del equilibrio de vaporización. En la práctica, las
columnas de destilación se diseñan hasta para valores de 150 etapas teóricas. Existen otras
expresiones más simples que calculan el Nmin
TCL TCP
N min (3.120)
3.TCP TCL
Dónde:
Ti Punto de ebullición de claves
Se despeja d o b para hallar los componentes no claves a partir de las ecuaciones (3.114
hasta 3.120)
77
f CL
bi (3.121)
d
1 R CL , R P min
N
bR
d
f CL R CL , R P min
N
di bR (3.122)
d
1 R CL , R P min
N
bR
Finalmente podemos representar como:
d CL f CL bCL
Para sistemas binarios con volatilidades constantes la razón de reflujo mínimo (Rm)
corresponde a la intersección común entre línea de operación de la zona tope, la línea de
equilibrio y la línea de alimentación. A partir del gráfico podemos hallar las siguientes
relaciones matemáticas:
Rm XD
Y .X (3.123)
Rm 1 Rm 1
.X
Y (3.124)
1 1. X
q Z
Y .X F (3.125)
q 1 q 1
Rm
1 XD
1 X D
1 X F
Underwood (3.126)
1 X F
78
1 . X D 1 X D
Rm 1
1 X F
Underwood (3.127)
1 YF
1
Rm (3.128)
1. X F
b.1) Vapor saturado
1
Rm 1 (3.129)
1.YF
7.2 MEZCLAS MULTICOMPONENTES (mezclas tipo 2)
Las ecuaciones anteriores son muy conservadoras y es aplicado para casos generales y para
algunos casos específicos se pueden utilizar las siguientes ecuaciones:
7.2.1. Cuando uno de los claves, no se distribuye la relación será
1 X CL , D X CP , D
Rm . (3.130)
1 X CL , F X CP , F
n CP . X F ,CP
CP 1 CP
1 q (3.131)
o
n CL . X CL , F
CL 1 CL
q 1 (3.132)
.Z
q 1
CL , CL , F
(3.133)
CL ,
Dónde:
q = condición térmica de la alimentación
=componente clave pesado
79
La ecuación (3.131) y (3.132) se resuelve para hallar entre los valores de para los
claves ligeros y pesados.
Cuando solo los componentes claves se distribuyen la ecuación (3.133) se resuelve para
hallar (que debe satisfacer CL ,CP 1 ) y este valor se reemplaza en la siguiente
expresión para hallar Rm
n X D ,CP
Rmin 1
CP 1
CP
CP
(3.134)
.X
O R min 1
CL , CL , D
(3.135)
CL ,
CL CP
(3.136)
2
Infinitas
CTA Costo Total
N soluciones
US$/año
S COP
CAP
Nm
Rm R Rm Ropt R
Figura 3.11. Razón de reflujo óptima
De los gráficos, la experiencia indica que los valores óptimos para el reflujo mínimo se
encuentra entre:
R
1,03 1,3 (3.137)
Rm
Sin embargo:
a) Cuando R crece N decrece el CAP de la columna disminuye
b) La pendiente de la curva de costo total usualmente es de muy baja inclinación cerca
del punto óptimo.
80
Como el valor óptimo varía muy poco se usa como regla heurística la razón de reflujo de
R
operación óptima para una columna de destilación: 1,2
Rm
R
1,2 (3.138)
Rmin
N f Z i , F , X i , D , X i , B , q, R,
Dónde:
Z i , F = Alimentación
q = Factor relacionado con la condición térmica de la alimentación.
R = Reflujo
= Volatilidad (cte.)
N m f X CL , D , X CL , B , (3.139)
Rm f Z CL , F , X CL , D , q, (3.140)
Debe cumplirse que:
N f N m f1 X CL , D , X CL , B , ; Rm f 2 Z CL , F , X CL , D , q, (3.141)
Esto tiene cierta justificación para mezclas binarias, pero para mezclas multicomponentes,
existe incertidumbre de la influencia de los componentes no claves.
81
La correlación empírica más famosa y simple corresponde a la desarrollada por Gilliland
(1940) en forma gráfica.
Para reflujo tipo 1:
1.0 1.0
N Nm N Nm
Y Y
N 1 N 1
0.1
0,0 0.0
0.01 0.1 1.0 0.0 1.0
R Rm R Rm
X X
R 1 R 1
Figura 3.12. Correlación gráfica de Gilliland
Gilliland a partir del estudio de numerosos casos ha notado que cuando la razón R/Rm es
aproximadamente 1.2; se obtiene una relación muy simple entre N Y Nm
ln X D 1 X w
1 X
D X w
2
ln s
N 2 N m 2. (3.142)
ln ln
1 F F d Ci
CL ,CP f
d CP
[Link] LR L LF , CL Lf , CP
1
CL CP F
1 1
Rm Y Rm 1
1X F 1.YF
R = 1,20 Rm
82
Utilizando el concepto de eficiencia y el gráfico de O’conell.
N Teórico
E0
N Re ales
E0
Parte lineal
F
Figura 3.13. Correlación para determinar la eficiencia de columnas
0.4983 0.5
Del gráfico. E0 (3.143)
. F
0.252
. F 1
4
Para valores de entre 1.10-3.0; E0 varía entre 63% y 51%, por lo tanto, para presiones
moderadas podemos usar:
E0 0.5 (3.145)
Entonces una estimación rápida es:
N Teorico N Teorico
N Re ales N R 2N (3.146)
E0 0.5
ln s ln s
N Re al 2 2 N Re al 4 Fenske (3.148)
ln ln
0.206
N R Z CP , F B
2
X CL , B
(3.149)
N S Z CL , F X
CP , D
D
Dónde:
NR = Números de platos arriba del plato de alimentación
NS = Número de platos en la sección de agotamiento
N N0 V V0
CTA CTA0 K1 K2 (3.150)
N0 V0
Dónde:
Ki = coeficiente lineal del costo del número de etapas multiplicado por el costo anualizado
de la columna para el caso base.
1) Según Gilliland
N 2N min (3.153)
2) Según Fenske:
X CL , D X CP , B
ln
X
ln s CL , B CP , D
X
N min (3.154)
ln P
1
2
ln CL CL
CP CP
O
X CL , D X CP , B
ln
ln s CL , B X CP , D
X
N min (3.155)
ln P ln P
84
3) El número de platos reales se puede calcular:
V = (R+1) D (3.158)
R = 1,2 Rm (3.159)
6) Razón de reflujo mínimo según Underwood para una mezcla binaria y alimentación
como líquido saturado es:
1 X CL , D X CP , D
Rmin F. (3.160)
F 1 X CL , F X CP , F
1 X CL , D
Rmin . (3.161)
F 1 X CL , F
Reemplazando (3.160) en (3.159)
1 X CL , D
R 1,2 . (3.162)
F 1 X CL , F
X CL , F
D .F (3.163)
X CL , D
1,2 X CL , D X CL , F
V . 1 .F
F 1 X CL , F X CL , D
1,2 X CL , D X CL , F X CL , F 1,2 X CL , F
V . . F .F .F .F
F 1 X CL , F X CL , D X CL , D F 1 X CL , D
85
1,2 X CL , F
V F (3.164)
F 1 X CL , D
ln s 1,2 X CL , F
CTA A0 A1 A2 .F (3.165)
ln F 1 X CL , D
Dónde:
A0 = CTA0-K1-K2
K
A1 = 4 1
N0
K
A2 2
V0
Dónde:
CC= Costo de la columna
CCPN= Costo del condensador
CR=Costo del hervidor
CCA= Costo anual del agua de enfriamiento
CST= Costo total del vapor
FRC=1/3 año 1 que implica i=15% y vida útil entre 11-15 años
0.533
M & S M G .R.T 2
1
86
h1
Sistema de
Separación
h2
h N
H H0 (3.168)
E0
Dónde:
H 0 h1 h 2
En la práctica:
h1 5 10 ft h 2 5 10 ft
Expresando como porcentaje el espacio adicional (H0 = 15%) para evitar el arrastre de
vapor y para el líquido del fondo, la ecuación (3.168) llega a ser:
N
H 1,15h (3.169)
E0
N 0,5 N Re ales
Si reemplazamos que: E0=0,5 N 0,5N Re ales
E0 0,5
4 ln s
Fenske: NR (3.171)
ln
Dónde:
s = Factor de separación entre componentes claves
= Volatilidad relativa
87
Reemplazando (3.171) en (3.170)
4 ln s
H 1,15h (3.172)
ln
1
4 M G .R.T 2
1 2
1
DC .G 2 (3.174)
0,88. .F P
M .T 14 1
DC 0,0296 G .G 2 (3.175)
P
1
m Kg 2
KJ
b) Segundo caso: F 6590 3 R 8,314
y
h m Kmol.K
M .T 14 1
DC 0,0252 G .G 2 (3.176)
P
Límites de diseño
H H 175 ft 53 m
20 30 (Diámetro menor a 0.46 m difícil de construir)
D D 1,5 ft 0.46m
0.533
.1.15h
4
K CL K C
0,82
Si hacemos: (3.177)
0,88. .F
0, 2665
M &S M .R.T
CC K CL G 3.30 FP .N 0.82 .G 0,533 (3.178)
280 P
88
a.2.1) Diseño del condensador
TD 55°C
TAS=49°C
TAe=35°C
Entonces el flujo de calor necesario para enfriar cuando no hay cambio de fase es.
Q mCP T
Q wC C PA TAS TAC (3.180)
Q h. A.TS Tw
3) Balance de calor o ecuación de diseño para dos fluidos que siguen la transferencia de
calor por conducción según la ecuación de Fourier:
Q U . A.T
Q U C . AC .TC (3.181)
Dónde:
1 1 1
(3.182)
U C h Ai h0
i
A
T TAC
TC AS (3.183)
T TAC
ln D
TD TAS
Dónde:
89
hi y h0 = coeficientes que se puede determinar experimentalmente o mediante análisis
dimensional cuando se usa vapor.
4) Ecuación de diseño
H DV U C . AC .TC (3.184)
H DV
AC (3.185)
U C TC
H DV
H DV wC C PA TAS TAC wC (3.186)
C PA TAS TAC
Dónde:
TD = Punto de burbuja del destilado
UC = Coeficiente global de T.C. para el condensador (2000 kJ/h. m 2 .ºC)
Wc = Flujo de agua (kg/h)
V = Flujo de vapor (kg/h)
CPA = Capacidad calorífica (KJ/kgºC)
H D = Calor sensible (KJ/kg)
M &S
CCON K HE . A 3.30 FC
0 , 65
(3.187)
280
0 , 65
M &S H D
C CON K HE . 3.30 FC .V 0,65 (3.188)
280 U C .TC
Dónde:
KHE= Costo de compra condensador= 475 US $ m 1, 30 101,3 pieUS$
1, 30
H D
CCA C A .H . V (3.190)
C PA TAS TAC
Dónde:
CA = Costo del agua US$/Kg
H = Operación h/año
H w .V U R AR TR (3.192)
H w .V
AR (3.194)
U R H S
H w v
wS H w BTU / Kg (3.195)
H S
M &S
CR K HE . AR 3.30 FC
0 , 65
(3.196)
280
91
0 , 65
M &S H w
CR K HE . 3.30 FC .V 0,65 (3.197)
280 U R .TR
Dónde:
KHE = Costo compra hervidor = 475 US$ m1,30
101,3US$ pie 1,30
BTU
U R .TR =127750 KJ / m 2 h 11250 (límite para el flujo de calor en el hervidor)
pie 2 h
U R BTU ; T º F
h. pie 2 .º F
C ST C S .H .wS (3.198)
H w
C ST C S .H . V (3.199)
H S
Dónde:
CS = Costo de vapor US$/Kg
H = Operación h/año
H S = 933BTU/lb
0, 2665
M & S K CL M G .R.T
Hacemos: C0 .3,30 FP (3.201)
280 3 P
M & S K HE H D
0, 65 0, 65
H w
C1 . . 3,30 FC (3.202)
280 3 U C TC U R TR
92
C A .H D C .H w
C2 H . S (3.203)
C PA .TAS TAe H S
Q U . A.Tm (3.204)
Q
A
UTm (3.205)
Dónde:
Q = Flujo de calor
U = Coeficiente global de transferencia de calor
Tm = Gradiente medio logarítmico de temperatura
A = Área del intercambiador de calor
1.1.1. Determinación de U
[Link]. Determinación de U a partir de los coeficientes individuales de película
1 1 Ai
U . (3.206)
hi h0 A0
a) los coeficientes interno (hi) externo (h0) se pueden determinar en forma experimental
(vapor)
b) mediante análisis dimensional
c) mediante correlaciones empíricas como la de: Dittas-Boelter
0 ,8 1 0,14
hi .D D.G C P . 3
0,023 (3.207)
K K w
Dónde:
K = Conductividad térmica
CP = Calores específicos
93
= Viscosidad de los fluidos (carcasa y tubos)
Si el proceso estudiado es rentable, entonces se efectúa diseños rigurosos del cálculo de los
coeficientes (U).
M &S
C1 .K I . A 3,29 FC
0,65
(3.208)
280
0 , 65
M &S Q
C1 .K I . 3,29 FC (3.209)
280 U .Tm
Dónde:
US $ US $
KI = Costo de compra del intercambiador = 101,3 1,30
475
pie m1,30
FC = FP (para intercambiador de acero al carbón)
94
T1 Enfriamos
Condensamos
Q1
TC
Subenfriado
Q2
TWS Q3
T2
t1
t2
Twe
a) Cálculo de los calores para las tres secciones de la corriente del proceso del
intercambiador de calor.
Q1 F .C PV T1 TC (3.210)
Q2 F .H V (3.211)
Q3 F .C PL TC T2 (3.212)
QT Q1 Q2 Q3 (3.213)
QT wS .CPw t ws t we (3.214)
Entonces:
QT
wC
C Pw t ws t we
(3.215)
Q1 w .C t t Q t t
C Pw ws 1 1 ws 1 (3.216)
QT wC .C Pw t ws t we QT t ws t we
95
Q3 w .C t t Q t t
C Pw 2 we 3 2 we
wC .C Pw t ws t we
(3.217)
QT QT t ws t we
4.2 A partir de la ecuación (3.216) y (3.217) determinamos t1 y t2 con los demás valores
conocidos.
e) El cálculo de las áreas para las tres secciones del intercambiador se calculan por
prueba error a partir de las ecuaciones de diseño.
Q1 U 1 . A1 .
t1 t ws TC t1 (3.218)
T t
ln 1 ws
TC t1
Q2 U 2 . A2 .
t1 t 2 (3.219)
T t
ln C 2
TC t1
TC t 2 T2 t we
Q3 U 3 . A3 .
T t 2
(3.220)
ln C
T2 t we
AT A1 A2 A3 (3.221)
Si la sección de condensación es grande comparado con las demás 2 secciones y está llega
a ser un 80 o 90% del QT, entonces el tamaño del intercambiador se puede aproximar
mediante la siguiente relación.
t ws t we
QT U . A. (3.222)
T1 t ws
ln
T2 t we
Se utiliza los mismos conceptos de diseño, pero que es más simple debido a que la
temperatura del medio calefactor es generalmente vapor de agua que se mantiene constante
a través del intercambiador.
96
T
TC
Punto de
ebullición
Q2 Q3
Sobrecalentado
Q1
Vaporización
Ebullición
97
CAPITULO IV: DISEÑO Y OPTIMIZACIÓN DE SUBSISTEMAS
Existe una gran cantidad de procesos que se desean calentar un flujo de alimentación y al
mismo tiempo deseamos enfriar un producto. Estas dos operaciones (calentar y enfriar)
podemos combinar y ahorrar tanto el vapor de calefacción como el agua de enfriamiento,
instalando intercambiadores de calor llamado PRECALENTADORES y
PREENFRIADORES
ts, Ps
F, CP, t1
t0
RECUPERACIÓN (R)
T0 AR
Calentador Ac t2
TAe
TAs F, CP, TR
Agua Ae Condensador
TF
Reactor (R)
Columna
TF t Ae 10º F TF t Ae 10 90 10 TF 100º F
Q F .C P t 0 t1 (4.1)
98
T0
TF
10°F
tAe=90°F
tAs=220°F
Para líquidos orgánicos-agua cerca del punto de ebullición se debe recuperar entre 75-
80% del calor disponible; entonces para sistemas líquido-líquido la ecuación (4.3) se
puede representar como:
e. Recuperación óptima
t 0 t1 r TR t1
Dónde: r = recuperación
99
Dónde:
C SH = Costo de vapor de calentamiento.
Cws = Costo de agua de enfriamiento.
g. Compromiso económico
t 0 TR CE
AR Entonces CC Cmínimo
CR CSH
Cwc
CTA C AR AR C AC AC C Ae AE C wc wc C SH S
CA T1 TS CA
US US
CA CA
UC C C UC C C
w S w S
TS t AS T1 H S T0 t AS T1 H S
100
ts, Ps
F, CP, t1
t0
RECUPERACIÓN (R)
AR Calentador Ac t2
T0
F, CP, TR
Caldera Ws
Ts, AHs
Reactor
T1
TAe
TAs
Ae Condensador
TF
Columna
De una planta de secado se obtiene un flujo de aire a 25 ºC y 1 atm., con un flujo que
contenido 300 kg./h de a cetona con una concentración de 1,5% molar.
Solución
1. Potencial Económico
101
V, aire?
Refrigerante
Acetona-aire
25°C, 1 atm
Flash
300 Kg/h
CONDENSADOR
L, acetona?
El vapor y líquido que sale del tambor están en equilibrio termodinámico por lo tanto:
Presión parcial de acetona en fase vapor = presión vapor de acetona líquido
B
T C
A ln P 0 (2)
Donde.
T=K B = 2940,46
P = mmHg C = -35,93
A = 16,6513 P = 7,5x 10 5 atm = 0.057 mmHg
Reemplazando:
35,93
2940,46
T T 187 K
16,6513 ln 0,057
102
Una alternativa para condensar la acetona es aumentar la presión de alimentación, hasta
que la presión parcial de la acetona sea igual a la presión vapor a 25 ºC. Entonces de tablas
encontramos la presión de vapor de acetona a 25ºC que es PA = 0,33 atm.
De ecuación (1)
PA
PA PT .YA PT
YA
0,33
PT PT 4400 atm PT 4546 Psia
7,5 x10 5
Refrigerante V, aire?
-100°C -86°C
Acetona-aire
25°C, 1 atm
Flash
300 Kg/h -86°C
CONDENSADOR
-100°C
L, acetona?
Si no deseamos enviar a la atmósfera aire a tan baja temperatura (-86ºC) que causaría
condensación y formación de hielo, podemos instalar un preenfriador y así ganar un
beneficio extra al reducir la carga del sistema de refrigeración.
V, aire?
-86°C -100°C -86°C
PREENFRIADOR
Acetona-aire
25°C, 1 atm t1 t2
Flash
300 Kg/h -86°C
CONDENSADOR
Aire
20°C L, acetona?
-86°C, Aire
103
a. Planificación del sistema de refrigeración
El rango de operación del proceso (-86ºC, -100ºC) hace necesario encontrar un refrigerante
a este nivel tan bajo de temperatura, y a partir de diagramas de presión-temperatura (Perry,
figura: 12-30) solo algunos compuestos están disponibles. Si consideramos la menor
temperatura del sistema de refrigeración y que sea ligeramente superior a la atmosférica
(1.36 atm, 20 psia) de modo que no se permita ingreso de aire al sistema de refrigeración,
solamente Etileno (-101ºC, 17 Psia) y Freón-14 (-101ºC, 70 Psia) son aceptables para esta
operación. Pero, además necesitamos que el sistema de refrigeración entregue el calor
extraído al agua de refrigeración a 33 ºC, lo cual significa 38 ºC para el refrigerante (
T 5º C ). Como ni el Etileno ni el Freón-14 puede lograr esta temperatura, se debe usar
un sistema de refrigeración de dos etapas.
Propano Etileno
T = -41ºC T = -101ºC
K=1 K=1
P = 17 Psia P = 17 Psia
T = 38 ºC T = -35 ºC
K=1 K=1
P = 190 Psia P = 250 Psia
Se puede estimar el flujo de material de las corrientes de proceso, aunque no se puede para
los fluidos de refrigeración, figura 4.7.
104
Reemplazando los valores en la ecuación (4)
V, aire?
-86°C -86°C
CONDENSADOR
Acetona-aire
25°C, 1 atm t2 -86°C
Flash
300 Kg/h
PREENFRIADOR
T=-101°C
Aire P=17 Psia
20°C L, acetona?
T=-35°C
P=250 Psia
EVAPORADOR
T=-41°C
P=17 Psia
T=-38°C
P=190 Psia
EVAPORADOR Agua de
49°C 33°C enfriamiento
Kg Kmol KJ
Q pre 29 339,5 0,997 20 86
Kmol h Kg º C
Entonces el calor total es:
KJ
Q pre 1040492
h
2) Calor de Enfriamiento
El punto de rocío de la mezcla alimentada se estima a partir de la presión de vapor
P 0 PT Y 1,36 atm y se obtiene T = -29ºC, Entonces el calor requerido para enfriar la
alimentación hasta el punto de rocío será:
105
Qenf 300
Kg
.2,12
KJ
339,5
Kmol
.0,997
KJ
.29
Kg
25 29º C
h Kg h Kg º C Kmol
KJ
Qenf 530698
h
3) Calor incremental
El calor incremental es la carga usada para enfriar la corriente de alimentación y condensa
algo de acetona.
Qinc H S H C (7)
Calor Calor
Sensible Latente
Simplificando:
Formando una tabla para las variables y por iteración podemos determinar los calores.
Unidad Q t1 t 2 t m U A
MKJ/h ºC ºC ºC KJ / h.m 2 .º C m 2 x10
Pre-enf 1,04 6 27 13 400 200
106
Ref-Cond 0,26 39 12 23 400 30
Eva-Etileno 0,43 86 6 29 400 40
Eva-Propano 0,67 8 6 7 400 240
5) Condensación de Refrigeración
QR 250848 11820
KJ
QR 262668
h
Cuando la potencialidad del proyecto sea establecida claramente, el proyecto estará listo
para la fase de desarrollo. En este punto, puede construirse una planta piloto o una planta
comercial. Una planta piloto es una pequeña replica a escala de la planta final, mientras
que un desarrollo comercial de planta es usualmente hecho sobre la base de piezas
dispersas de equipo las cuales han sido usadas anteriormente.
107
El PFD describe la ruta del proceso, mostrando los flujos de material y energía entre
aquellas unidades de proceso que constituyen la planta.
Después de haber concebido, definido y asignado el problema, la solución del mismo rara
vez es obvia. Frecuentemente hay numerosas posibilidades y suposiciones potenciales.
Aún con esta incertidumbre, en general es posible determinar un diagrama de flujo, esto
conduce a suposiciones que pueden reducir la complejidad del problema. Al efectuar los
balances de materia y energía, a menudo, se exponen las variables más importantes del
proceso. Es una manera eficiente de familiarizarse con el proceso. Se identifica donde falta
información. Si se efectúa y termina adecuadamente el diagrama de flujo contendrá los
datos necesarios para el diseño de partes individuales de equipo. En general aunque se
modifiquen después, es relativamente fácil corregirlo y repetir los cálculos de diseño.
Cuando los estudiantes no sepan por donde un nuevo proyecto, se les aconseja comenzar
por el diagrama de flujo, con la confianza de que por medio de este ejercicio se evaluará un
método potencialmente fructífero. La mecánica de la preparación del diagrama de flujo se
describe más adelante.
Ejemplo 1:
Cuál es la presión de diseño para una cañería NPS 42 con un espesor de pared de 0,354” y
un factor de localización de 0,90. El caño es del tipo 5LX-x80, y la temperatura del flujo
de gas es menor que 250ºF. Asumir un factor de diseño de F=0,80.
Solución:
El tipo de caño 5LX tiene un factor de junta o soldadura de 1.0 y el grado tiene una mínima
tensión de fluencia de 8000 Psia.
2.S .t
Pd .F .L.J .T
D0
Datos:
Pd =?, psia F = 0.80
S = 8000 Psia L = 0.90
t = 0.354 pulg. J = 1.00
D0 = 42 pulg. T = 1.00
2.(80000)(0,354)
Pdiseño .(0,80).(0,90).(1,00).(1,00) 971 Psia
42
Ejemplo 2:
Que diferencia existe es el costo del acero para un NPS 42, L 900 Km. De línea de
transmisión usando x70 o x80 grado del caño. Asumir que la presión de operación máxima
de la línea es 1440 psia, factor L=0.90 y factor F=0.80, la densidad de ambos aceros x70 y
x80 es 7,801ton m3 . El costo del caño en el lugar es: x70 = US$1000/ton; x80 =
US$1030/ton.
108
Solución
a) para el acero tipo x70:
El espesor del caño para el acero de grado x70 es:
t
P.D0
144042 0,60 pu lg .
2.S .F .L 2700000,800,90
Diámetro interno: Di D0 2t Di 42 2(0,60) 40,8 pu lg
Reemplazando en ecuación (3.12) ws .L..t.( Di t )
t
P.D0
144042 0,525 pu lg .
2.S .F .L 2800000,800,90
Diámetro interno: Di D0 2t Di 42 2(0,525) 40,95 pu lg
Reemplazando en ecuación (12) ws .L..t.( Di t )
Ejemplo 3:
Solución:
Kg 1año 1d Kg carbon
QT 126415000 . x QT 14431
año 365d 24h h
109
Kcal
Poder calorífico carbón importado: PC 7500
Kg
Kg Kcal Kcal
Entonces: QT 14431 x7500 QT 108232500
h Kg h
Kcal
Pasando a utilizar gas natural; entonces el PCGN 9300
h
108232500 Kcal h m 3 gas
Luego: Q Q 11638
9300 Kcal m3 h
Qb Pb T f L 0,539 G 0.461
1
2.618
D (2)
2.450Tb P12 P22
0.539
Dónde:
MCF
Qb G = Gravedad específica gas
h
Pb Psia Tb = ºR
Tf º R P1 = Psia
L = ft P2 = Psia
Datos
Condiciones estándar
110
Psia
Pb 101,325Kpa. 14,7 Psia
6,894733
a) distancia Lima-Huancayo por carretera es 310 Km.
Tf = 520ºR
b) presiones de operación: 147 bar
Presión 40 bar
G = 0,6 (promedio-diseño)
c) Para distribución: 60 bar-40 bar
9
Tb = 15ºC + 273 = 288K. = 518,4ºR P1= Presión cabecera (gaseoducto)
5
P2 = Presión sección (Planta)
m 3 gas ft 3 ft 3 MCF
Qb 10102 x 35,314 3 356742 Qb 357
h m h h
m ft
L 310 Km x 1000 x 3,2808 L 1017048 ft
Km m
14,5038Psia
P1 100bar x 1450 Psia P1 1450 Psia
1bar
14,5038Psia
P2 80bar x 1160 Psia P2 1160 Psia
1bar
1
357 14,7 520 x10170480.539 0,60.461 2.618
D 6
2.450518.4 1450 2 1160 2
0.539
2.S .t
Pdiseño .F .E.T (3)
D0
Reemplazando: Pdiseño
2. 650000,126 .0.72. 1. 1 P 1966 Psia
d
6
111
Esta presión es la máxima que debe resistir el caño durante la prueba hidráulica. Esta
presión debe ser siempre mayor (MAPO).
ws .L..t.Di t (4)
Datos:
L = 310000 m
= 7,801 ton/ m 3
m
t = 0.126 pu lg .0.0254 0,0032m
pu lg
Di D0 2t 6 20.126 5,75 x 0,0254 0,14605
C S 3450
US $
caño de acero x65 puesto en Hyo
ton
Luego:
3450US $
C A 3629ton 12520000US $ (50% del cos to total)
ton
25040000US$
Sabiendo que el costo del caño es de: 14 US$/pulg, entonces para el caso anterior
tendremos que:
C A 14x6x310000 26040000US$
112
PARTE II: DISEÑO CONCEPTUAL DE PROCESOS QUÍMICOS
Para diseñar un nuevo proceso de producción se debe tener en cuenta los siguientes
aspectos ambientales:
Regulaciones ambientales cada vez más estrictas
Para tener éxito se debe buscar nuevas rutas en el diseño de procesos, la investigación y
la educación (diseño de tecnologías ambientalmente compatibles con el diseño,
investigación y educación)
Transferencia de tecnologías limpias de naciones industrializadas a países en desarrollo
La industria química está tomando nuevos rumbos a productos de alto valor y de bajo
volumen. (Química Fina, Tecnología Limpia, Biotecnología, etc.)
El método ingenieril nos indica que se resuelve los problemas de diseño primordialmente
desarrollando soluciones muy simples y se va adicionando mayores detalles. Un buen
sistema es la simplificación del flowsheet para atacar los problemas de diseño.
113
Tolueno + H2 Benceno + Metano
2 Benceno Difenil + H2
En la figura 5.1, se presenta el flowsheet inicial. Como en la corriente del reactor aparecen
H2 y CH4 se puede usar un condensador parcial, para que los compuestos aromáticos se
condensen y sean fácilmente separados por flash.
Calentamiento Compresor
H2, CH4 Purga
Calentamiento
Calentamiento Calentamiento
H2, CH4
Benceno
Tolueno
Estabilizador
Producto
Reciclo
Difenil
114
5.2.5 Modificación de flowsheet
Diagrama de los procesos e integración energética que se adiciona a los flujos y composiciones de cada corriente de proceso, además se incluye
las temperaturas de cada corriente de proceso figura 5.2. Con la integración energética se percibe una máxima recuperación de la energía. Sin
embargo este proceso de integración es complicado por ello se realiza al final del análisis.
Purga
H2, CH4
Columna de Estabilización
Columna de Benceno
Columna de Tolueno
Compresora
Reciclo
Gas
Difenil
CW Alimentación
H2
Alimentación
Tolueno
115
5.2.6 Flowsheet simplificado
Sistema Purga
recuperación vapor
H2, CH4
H2, CH4
Sistema
separación fases
Tolueno
Sistema de
reacción
Benceno
Sistema separación
líquido
Difenil
Figura 5.3. Sistema de separación del proceso (HDA)
Reciclo gas
Purga
H2, CH4
Tolueno Difenil
Sistema
Sistema de separación
reacción
Reciclo de Tolueno
Figura 5.4. Estructura del reciclo del proceso (HDA)
Una vez simplificado el flowsheet o realizado la síntesis de los procesos y operaciones del
proceso de producción, o para realizar el diseño de un muevo proceso se realiza mediante
una metodología sistematizada que nos permite en forma simple identificar los posibles
problemas de contaminación y la producción de residuos peligrosos en cada una de las
etapas de los procesos de producción. En consecuencia, una producción limpia nos indica
“No Generar” residuos o desechos.
116
Ejemplo: A + B Producto
A + Producto Sub Producto
Recuperación Tratamiento
(Costo de recuperación) (Costo de tratamiento)
Recuperar Recuperar
y reciclar y remover
a) Acuosos a Tratamiento
Tratamientos b) Sólidos orgánicos a Incineración
c) Sólidos a Vertederos
117
Reciclo-Gaseoso
Purga o Venteno
Gas
Producto
Alimentación Sistema de
Separación Sub-Producto
Reciclo-Líquidos
Figura 5.5. Estructura general de un diagrama de flujo
Se trabaja con datos disponibles y posteriormente se usan los resultados para justificar y
obtener datos adicionales.
118
Estequiometria de todas las reacciones que tienen lugar
El rango de temperatura y presiones para las reacciones.
Las fases del sistema de reacción.
Información sobre la distribución del producto versus conversión.
Información sobre conversión versus velocidad espacial o tiempo de residencia
Estado del catalizador: Desactivación y regeneración
Numerosos procesos se diseñan para operar a las condiciones de máxima conversión, pero
esta operación normalmente no corresponde a la óptima conversión económica. Entonces
se debe estudiar la relación entre máxima conversión versus conversión óptima.
A B C
Producto deseado Producto no deseado
1
A
P Concentración
B
Máxima
Q
0
Reciclo Tiempo
Moles _ de _ B _ Pr oducidos
Selectividad
Moles _ de _ B _ Convertidos
El costo de materia prima está en rango de 35-85 % del costo total del producto.
119
La conversión óptima es fijado por un cambio entre pérdida de selectividad y el costo
elevado del reactor y con las conversiones altas contra el costo de reciclo y
conversiones.
Para diseñar una planta nueva o en ampliación se realiza de acuerdo a las condiciones del
mercado. Como una primera suposición de la velocidad de producción se considera una
planta agrade; pues la producción a gran escala es menor costosa que las materias primas y
los servicios tienen menor costo; y con esto se obtiene grande beneficios.
El tamaño máximo de una planta es usualmente fijada por el tamaño máximo de uno o más
piezas de los equipos que se pueda construir o transportar.
El mercado fija la pureza del producto. Sin embargo, se puede producir de varias purezas a
diferentes precios. Una alta pureza requiere de grande costos y esto debe de comunicarse al
departamento de marketing.
Se debe tener información sobre todos los tipos de servicio disponibles en la planta o en
plantas anexas a la nueva construcción.
Servicios:
- Suministro de combustible.
- Niveles de presión a vapor.
- Temperaturas de ingresos y salidas del agua de enfriamiento.
- Niveles de refrigeración.
- Potencia eléctrica.
Facilidades para la disposición o tratamiento de desechos.
120
- Pesos moleculares.
- Puntos de ebullición.
- Presiones de vapor.
- Calores específicos.
- Calores de reacción.
- Densidades de líquidos.
- Coeficientes de fugacidad.
- Seguridad.
- Toxicidad.
- Impacto ambiental.
121
5.4.4 Operaciones múltiples en un simple reactor
Producto
Destilador
Alimentación
Reactor y separación de Calor
Catalizador
calor
Catalizador
Alimentación Producto
Calor Separador
Reactor
122
CAPITULO VI: ESTRUCTURA DE ENTRADA-SALIDA DEL
DIAGRAMA
Para decidir la estructura de entrada y salida se dibuja como si fuera una caja “negra” el
proceso total y focalizamos nuestra atención en las materias primas, productos
subproductos. Para decidir en este nivel debemos resolver las siguientes interrogantes.
¿Se debe purificar las alimentaciones?
¿Se debe remover o reciclar un subproducto reversible?
¿Se debe utilizar un reciclo gaseoso y purga?
¿Se deben recuperar ciertos reactivos y reciclarlos?
¿Cuántos flujos de productos se deben obtener?
¿Cuáles son las variables de diseño y los compromisos económicos?
Alimentación Producto
Alimentación Proceso Subproducto
Purga
Alimentación Producto
Alimentación Proceso Subproducto
123
6.1.2 Remover o reciclar un producto reversible
T + H2 B + CH4
2B D + H2
Reciclo:
Remover:
Si removemos el componente D, el T se convertirá en D o subproducto por lo tanto
existirá aumento de materia prima y por lo tanto aumento de costo de Tolueno.
Los reactivos gaseosos tienen menores costos que los líquidos orgánicos.
La refrigeración es una operación más costosa.
Las materias primas como el agua y el aire no se recuperan porque son de bajo costo que
los demás compuestos orgánicos (más bien se adicionan en exceso).
Par conocer el número de flujos de los productos, se ordenan los componentes de acuerdo a
su punto de ebullición (de menor a mayor); luego se asigna un código de destino y la
clasificación de componentes de acuerdo a la siguiente tabla.
124
La definición del problema no siempre es completa para determinar el o los balances de
materia total, las corrientes de flujo y sus costos. Solamente un proceso con una sola
reacción y sin pérdida de reactivos no tiene grados de libertad.
Por lo que, los balances de materia y los flujos y sus costos se desarrollan en función de
variables de diseño, tratando de encontrar los valores óptimos de estas variables.
Se desarrollan expresiones para los balances de materia total en términos de las variables
de diseño sin considerar los flujos de reciclo.
125
Calcular los flujos de salida de los reactivos en función del “exceso” de reactivos
que no se reciclan o recuperan (reciclo y purga o aire o agua)
Calcular los flujos de entrada y salida de las impurezas que entran con los reactivos.
Los balances, flujos y los costos se pueden realizar en una hoja electrónica. Las
condiciones de operación, las variables de diseño y las formulas se deben declarar con la
finalidad de facilitar los cambios y con ellos recalcular los flujos de proceso.
Los “mejores” valores de las variables de diseño dependen del factor económico. Por ello,
se debe calcular el costo de las materias primas y de los productos en función de las
variables de diseño. Para determinar el “Potencial económico” definido en el nivel 2 como:
PE2= Valor del producto + Valor del subproducto – Costo de materias primas
En caso de que el PE2 sea negativo, terminamos el estudio de diseño, y se busca una fuente
de materia prima de menor valor o tratamos de buscar otra reacción química que utiliza
materias primas de menor costo.
126
CAPITULO VII: ESTRUCTURA DE RECICLO EN UN DIAGRAMA
DE FLUJO
127
[Link] Reciclo gaseoso
No se estudia en diseño preliminar por los costos de los equipos son de menor costo
(bombas, tuberías, etc.)
Ejemplo 1: Considerar los siguientes componentes y sus destinos los cuales están
ordenados de acuerdo a su punto ebullición.
Reactor
Alimentación R1
acetona
Reactor
R2
Reciclo acetona
128
Componente Punto de Ebullición ºC Código de Destino
CO -191 Combustible
CH4 -161 Combustible
C2H4 -104 Combustible
Queteno -42 Reactivo inestable-conversión completa
Acetona 50 Reactivos-reciclan a R1-liquido
Ácido acético 118 Reactivo reciclan a R2-liquido
Anhídrido acético 139 Producto primario
Reactor
Acetona R1 Anhídrido
Acético
Reactor Separador
R2
Reciclo acetona
129
Condiciones
Mientras mayor es el exceso mayor es la selectividad. Sin embargo existe mayor costo
de recuperación y reciclo de isobutano.
El exceso óptimo se obtiene de un análisis económico.
CO + Cl2 COCl2
130
CH2CH2 +H2O CH3CH2OH
2 CH3CH2OH (CH3CH2)2 + H2O (1)
Supongamos que deseamos producir 783 mol/h del azeotropo EtOH-H2O, que contenga
85.4 % mol EtOH, la alimentación de etileno contiene 4% CH4 y agua pura.
C2H4, CH4
Etileno 4%
C2H4, CH4
Reactor
EtOH-H2O
H2O pura
Separador Yazeo = 85.4% mol EtOH
ETER
H2O
PEtOH = Yagua*Pagua
PEtOH = 0.854(783) = 669 mol ETOH/h (3)
Separador
FW.P+FW.R FW Fw(1-x) FW.P
Reactor
Fw (1-x) – FW.P
131
1.5 La cantidad de agua que sale con el producto: FW,P
FW FW , P FW , R FW 1 X FW , P
FW FW , P FW , R FW FW X FW , P
FW X FW , R
FW , R
X
FW
Después de haber estimado el flujo del reactivo limitante, se usa una razón molar a la
entrada del reactor para calcular los flujos de reciclo de los componentes.
No existe regla heurística para seleccionar un valor de “x” para el caso de reacciones
complejas.
No existe regla heurística para seleccionar la composición de purga “y” y Razón
Molar (RM).
Para reacciones simples un valor razonable es:
132
Para decidir cómo operara el reactor se determina la carga térmica y la temperatura
adiabática.
Para reacciones simples toda la alimentación fresca del reactivo limitante usualmente se
convierte en el proceso (solo se pierde pequeñas cantidades)
QR =ΔHR*FF (7.1)
(CH3)2CHOH (CH3)2CO+ H2
Una vez que se ha determinado la carga térmica del reactor y el flujo a través del mismo en
función de las variables de diseño se estima el cambio de temperatura adiabática.
QR = [Link].(TE,P-TR,S) (7.2)
Solución:
QR = 1.324*106 BTU/h
Remplazando datos:
133
1.324x106 51.3 22.0 2.1 0.922.0572 TR,S
Luego
Si no se puede aplicar operación adiabática, tal como el ejemplo anterior entonces aplicar
calor directo o enfriamiento. Siempre verificar el área superficial de transferencia de calor
que sea adecuado para el tipo de reactor.
BTU KJ
U 20 400 y T 50F 25C
h pie F
2
h m 2 C
Para una carga de calor de QR = 1*106 BTU/h se debe determinar el área de transferencia.
Q
Q U A T A
U T
1x10 6 BTU
A h A 1000 pie 2 100m 2
20 BTU x50F
h pie 2 F
Esta área es muy pequeña, porque la máxima área de transferencia de calor que es capaz de
transferir en un intercambiador de calor de cabeza flotante en la carcasa está en el rango de
6000 - 8000 pie2 (600 - 800 m2). La carga de calor mínima para remover o aplicar calor a
un reactor utilizando un simple intercambiador está limitado a un rango de 6 - 8*106
BTU/h.
134
7.4 LIMITACIONES DE EQUILIBRIO
Los flujos se determinan en función de las variables de diseño (x, RM, etc.). Si se remplaza
los flujos en la relación de equilibrio veremos que la conversión seleccionada se encuentra
sobre o debajo del valor de equilibrio. Si la conversión de equilibrio excede el resultado no
tiene significado. En muchos casos es necesario determinar el valor exacto de la
conversión de equilibrio x en función de las variables de diseño, lo cual aparece en el
modelo cinético que se utiliza para determinar el tamaño del reactor.
Si uno de los productos puede ser removido mientras la reacción sucede, entonces, una
reacción aparentemente limitado por el equilibrio puede ser forzada a una conversión
completa.
Ejemplos:
1) Proceso de ácido sulfúrico
2) Síntesis de amonio
Para altas temperaturas los reactores tienen pequeños volúmenes; entonces en estas
reacciones la conversión de equilibrio decrece al aumentar la temperatura en el reactor. Es
por ello, en estas reacciones frecuentemente se retiran calor mediante una serie de lechos
135
adiabáticos o intercambiadores de calor, enfriando los gases o mediante bypass de
alimentación en frio para bajar la temperatura entre los lechos.
7.4.4 Diluyentes
Para cambiar la conversión de equilibrio se puede usar la T, P y RM, así mismo en algunos
casos un componente extraño (diluyente) se adiciona para variar la conversión de
equilibrio.
Ejemplo:
Producción de estireno
Etilbenceno Estireno H2
Etilbenceno Benceno + Etileno
Etilbenceno Tolueno + Metano 1100 °F, 20psia
La adición de vapor (o metano) al ingreso del reactor baja la presión parcial del estireno e
H2 y así disminuye la velocidad de la reacción reversible primera. El vapor sirve en parte
como un transportador de calor para suministrar calor endotérmico a la reacción.
136
CAPITULO VIII: SISTEMAS DE SEPARACIÓN
Purga
¿Cuál es la fase?
Alimentación Productos
Sistema
Reactor separación
LÍQUIDO
Alimentación Productos
Sistema
Reactor separación
líquidos
RECICLO LÍQUIDOS
137
a. Separar ambas fases en el reactor (vapor y líquido)
Si el efluente del reactor es una mezcla de dos fases, podemos usar al reactor y luego una
fase de separación o ruptura (o poner una columna flash después del reactor). Los líquidos
se van al sistema de separación de líquidos. En cambio la separación del vapor dependerá
de la temperatura de la corriente, si el reactor opera a temperatura del agua de enfriamiento
la corriente del vapor del reactor a 100 ºF y la fase líquida se va al sistema de separación
líquida.
Líquido
Líquido
Sistema reacción
RECICLO LÍQUIDOS
2.2) Si la temperatura del líquido obtenido de la separación flash es baja y esta contiene
todos los reactantes (más no los productos o subproductos), entonces, reciclamos todo el
flujo al reactor.
2.3) En cambio, si a esta temperatura los líquidos tiene otros componentes entonces se
envía al sistema de recuperación del vapor, para luego enviar al sistema de separación de
líquidos.
138
Purga
Sistema conversión
vapor
VAPOR
Reactor
LÍQUIDO LÍQUIDO
Nota: Si no se produce una fase split (ambas fases) en el reactor o a la salida, entonces
presurizar el sistema del reactor para obtener una fase fácil de separarse. Si no se precisa
entonces, se tiene la probabilidad de usar alta presión y condensación parcial refrigerada.
F=V+L (8.1)
139
FZi Vyi
xi (8.3)
L
y i k i xi (8.4)
FZi Vyi
yi k i (8.5)
L
De ecuación (8.1):
L=F–V (8.6)
FZ i Vyi
yi k i
F V
yi F V k i FZ i k iVyi yi F V k iVyi k i FZ i
Luego:
zi
yi (8.7)
V
F
1V
1
F k
i
O También:
zi
xi (8.8)
k j 1V F
1
Estas expresiones son iguales cuando existe una división o separación perfecta, para
componentes con un valor de ki entre 0.1 y 10.
140
La separación perfecta ignora el equilibrio líquido – vapor. Sin embargo, podemos
sobreponer esas relaciones de equilibrio para el flujo; desde este concepto la composición
líquida en equilibrio con un componente vapor en fracción mol:
fi fi
yi (8.13)
V f i
yi fi
xi (8.14)
ki ki f i
fi fi fi fi
li Lxi li (8.15)
ki f i ki f i
i f i li f i 1 f i
k i f i
Entonces:
i f i 1
i f
(8.16)
ki f i
k j f j fi
j (8.17)
f i
k j fi
l j f j 1 (8.18)
fj
Reciclo Purga
Gas
Si la ubicación del sistema de recuperación se realiza según los items 1º y 2º, entonces se
deben corregir los balances. Porque algunos componentes no considerados pueden
perderse en la corriente de purga o pueden ser recuperados en el reciclo del gas (esto
cambia el tamaño del compresor). Esto se realiza para corregir errores; sin embargo las
estimaciones anteriores nos proporcionan buenas aproximaciones. En este caso siempre se
considera el criterio ingenieril para corregir o no.
Si se cambia los balances, entonces debe cambiarse los flujos de proceso para no
incrementar los errores.
142
b.3) Adsorción
b.4) Procesos de separación por membranas
b.5) Sistemas de reacción
Estrategias
Ejemplo 1.- A un absorbedor de gas, se necesita suministrar solvente y en este caso debe
utilizarse un reciclo entre el sistema de separación líquido
Alimentación Sistema
separación fases
Reactor
Se necesita estimar el tamaño y luego el costo de los sistemas de recuperación que sean los
más baratos.
143
Las decisiones que se debe tomar para el sistema de separación de líquidos son:
1. Cómo las fracciones ligeras deben ser removidos si estos contaminan grandemente el o
los productos?
2. ¿Cuál debe ser el destino de las fracciones ligeras?
3. ¿Reciclaremos componentes que forman azeotropos conjuntamente con los reactantes,
o los separamos los azeótropos?
4. ¿Cuál de las separaciones puede ser hecho por destilación?
5. ¿Qué secuencia de columnas se necesita utilizar?
6. ¿Qué tipo de separaciones efectuamos si la destilación no es posible?
A. Fracciones ligeras
Condensador Columna de
parcial pasteurización
Columna de
destilación
Si las fracciones ligeras tienen poco valor económico se pueden evaluar a través de
respiraderos o aberturas desde el proceso y si además son contaminantes éstas se pueden
144
quemar. Si las fracciones ligeras son inflamables se debe recuperar para utilizarlo como
combustible. Sin embargo, si estas fracciones ligeras son valiosas, éstas se pueden tener en
el proceso o reciclados al sistema de recuperación de vapor generándose una nueva
corriente de reciclo.
B. Azeotropos o reactantes
C. Aplicabilidad de la destilación
Cuando las volatilidades relativas son menores a 1.1 se agrupa los componentes para seguir
una secuencia de destilación.
B
B, C Tarea de
C
diseño para
A
α
A 3.2
Juntos B 1.7
C 1.6
D 1.0 D
E 0.4
D, E
Para separar una mezcla de 3 componentes (no azeotropos) lo podemos realizar de dos
formas: recuperar primero los componentes ligeros o los componentes pesados y luego
separar los dos componentes restantes.
145
B A
A
A A
B B
C C
C C B
Directo Indirecto
Número de Componentes 2 3 4 5 6
Número de secuencias 1 2 5 14 42
La primera regla heurística está sustentada en el hecho de que el material para esta
columna es más costoso que el acero al carbono por los componentes corrosivos; por ello
la destilación es más costosa.
Las columnas al vacío son más costosas que las columnas a presión y también se debe
evitarse el incremento de los costos por limpieza.
Por otra parte, se prefiere remover productos del destilado para evitar la contaminación con
las corrientes de reciclo que contiene materiales pesados, mohos, etc., que se acumulan en
el proceso. Para prevenir la obstrucción del reboiler se purga una corriente de líquido
pequeña.
146
D.2 Secuencias de reglas heurísticas para columnas simples
Las principales reglas heurísticas para una secuencia selectiva de columnas simples son:
1º. Los más abundantes primeros.
2º. Los ligeros primeros.
3º. Alta recuperación en separaciones posteriores.
4º. Separaciones dificultosas posteriores.
5º. Favorecer separaciones equimolar
6º. Separaciones siguientes deberían ser baratos
Las reglas anteriores (1º-6º) tienen varias contradicciones (Ejemplo: Si el componente más
abundante es el pesado, el conflicto es entre el 1º y 2º), en estos casos se deben considerar
las variables de diseño.
Si se consideran dos flowsheets alternativos donde existen reciclos de líquidos a los dos
sistemas se tendrían costos distintos. De esta forma, para una conversión óptima,
usualmente se intercambia entre pérdida de selectividad y costo de reciclo.
Reactor
Calor Enfriamiento
Alimentación
Reactante
Flash
Productos Livianos al
final
Pesados al final
147
Reactor
Calor Enfriamiento
Alimentación
Reactante Flash
Livianos al
final
Destilación
Destilación
Destilación
Productos Pesados al final
En este caso se puede subdividir los problemas de secuencia en dos o más partes para
proceder como el caso de separaciones simples.
A. Extracción
148
Para separar una mezcla B y C en el punto 1, que corresponde a la alimentación (Figura
8.9) y el reciclo se requiere la alimentación de un solvente S en una columna de extracción
(punto 2). Por el punto 3 sale el componente 3 y el solvente S que debe recuperarse en una
columna de destilación obteniendo el producto del componente B (punto 5). La otra
corriente, en donde se recura el 99% o más del componente C se encuentra en el punto 4.
En otra columna de destilación removemos el solvente de la mezcla por el fondo, y por el
tope la mezcla C (+B) más concentrada que al inicio (punto 6). Nuevamente por
destilación normal fijamos la especificación del fondo para obtener el componente puro C
(punto 7) logramos separar en tope una mezcla como la original que es reciclado al punto1.
5
B
3 B+S
B+C
S
Columna
extracción
C(+B) 7 C
1
6
B+C
7 C
1
4
4
C+S(+B)
5
3 2
B S
Figura 8.9. Extracción
B. Destilación extractiva
149
B
S C
B
C B = HNO3
C = H2O
S = H2SO4
C
S
C. Destilación azeotrópica
B+C
azeotropo B+C(azeotropo)
B
B
B = Etanol C
C = Agua
S = Benceno
C
Este azeotropo térmico que sale del tope de la primera columna es una mezcla heterogénea
que se condensa. El extracto rico en benceno se utiliza como reflujo a la primera columna
y el otro extracto se alimenta a la segunda columna. En la segunda columna nuevamente
del tope sale el azeotropo ternario que se recicla al condensador y se recupera por el fondo
una mezcla etanol-agua que debe ingresar a la tercera columna en donde se recupera agua
pura y nuevamente por el tope se recupera el etano puro y el azeotropo agua que será
reciclado a la primera columna.
150
D. Destilación reactiva
B S
B,C = Xileno
B α = 1.03
C S = Organometálico
B,S = Complejo
α = 30
C
C+S
Figura 8.12. Destilación reactiva
De esta forma se separa gradualmente, pero el costo del sodio cumeno aumente. Sin
embargo si el costo de extraente es bajo, esta alternativa puede ser puede ser un
fundamental proceso de separación.
E. Cristalización
Vapor C(+B)
B. C
Cristalización Centrifuga Secador
Cristales de B
Líquido C
Figura 8.13. Cristalización
151
separadores (para las corrientes de purga), reactores, recipientes de vaporización
(divisiones de fase), absorbedores de gas (y/o otros unidades de recuperación de vapor), y
columnas de destilación (y/o otros sistemas de separación de líquidos). De esta forma
ahora se desarrolla el set de balances de materia riguroso.
152
BIBLIOGRAFIA
6. Prausnitz J., Eckert C., Orye R., O’Connelll J. “Computer Calculations for
Multicomponent Vapor-Liquid Equilibria”, Pretince-Hall, Inc. 1967
153