Síntesis y Economía en Ing. Química
Síntesis y Economía en Ing. Química
INGENIERÍA ECONÓMICA
2
EVER INGARUCA ALVAREZ
SÍNTESIS DE PROCESOS E
INGENIERÍA ECONÓMICA
3
PRESENTACIÓN
El área de diseño de procesos actualmente está reflejada por los resultados de los esfuerzos
científicos, soportados por equipos de cómputo cada vez más poderosos. En enfoque de los
conceptos inherentes a la síntesis de procesos es una parte medular de la educación del
ingeniero químico, pues integra los conocimientos de varias disciplinas y proporciona una
visión global del campo de la ingeniería química. El material contenido en este texto
refleja mi experiencia en la enseñanza del curso de síntesis de procesos químicos de la
Universidad Nacional del Centro del Perú.
El propósito de esta obra es brindar una introducción al diseño preliminar de los procesos
químicos y contribuir con la formación profesional de los Ingenieros de Procesos
Industriales
4
Ms. Ever Florencio Ingaruca Alvarez
5
AGRADECIMIENTO
DEDICATORIA
A Dios.
A la memoria de mis padres Pedro y
Flora.
A mi esposa Gloria, por todo su esfuerzo
y dedicación diaria.
A mis hijos Flora, Jones, Katherine y
Vladimir
A mis maestros, colegas y alumnos.
6
INDICE
CAPITULO 1: INTRODUCCION
1.1 Introducción al diseño de procesos químicos 1O
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 11: 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 111: 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
7
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
CAPITULO IX: SINTESIS DE REDES DE INTERCAMBIADORES DE CALOR
9.1 Calor mínimo y requerimiento de enfriamiento 158
9.2 El concepto de delta T mínima 158
9.3 Uso de diagramas de contenido de calor 159
9.4 El método del punto de pliegue 162
CAPITULO X: PREDICCIÓN DE AREAS EN REDES DE INTERCAMBIADORES
DE CALOR
10.1 Construcción de las curvas compuestas 169
10.2 Curvas compuestas balanceadas 171
10.3 La gran curva compuesta 171
10.4 predicción de aéreas de transferencia de calor 172
CAPITULO XI: ANÁLISIS DE REDES DE INTERCAMBIADORES DE CALOR
EXISTENTES
11.1 Diagnóstico de la red 176
11.2 El principio básico de la revisión 176
11.3 Establecimiento de objetivos 182
BIBLIOGRAFIA 184
8
INTRODUCCIÓN
9
PARTE 1: INGENIERÍA ECONÓMICA Y REGLAS HEURÍSTICAS
CAPITULO l. 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
B. ANÁLISIS (SIMULACIÓN)
ENTRADAS SALIDAS
C. OPTIMIZACIÓN
ENTRADAS SALIDAS
10
Tabla 1.1. Tecno ogía de diseño de eqmpos de procesos
TE ORlAS
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 • Equipo de recuperación de
Convección
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 • Reactores biológicos
Cinética química
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.
11
..
T a bl a 12 P.nnc1pa es ecuacwnes d e transporte
Cantidad Calor Masa Momento
Flux q Na 'trz
Transporte Molecular
8T acA avz
Fuerza Impulsora - --
az az 8r
Ley Fourier's Fick's Newton's
Conductividad Difusividad Viscosidad
Propiedad
Térmica
KT DA J1
Transporte Total
Fuerza Impulsora ~T ~CA ~p
M 1
Relaciones q = hr ·~T NA= k¿ ·~CA f=g ·D·-
e L 2· p·v 2
* 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.
12
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
13
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.
Proceso Recreativo
PROCEDIMIENTO BOSQUEOA
PRELIMINAR LAPIZ
1
ALTERNATIVAS DE
l
MODIFICACIONE
PROCESO S PRELIMINARES
DIS~L
ESTIMACION DE
COSTOS RIGUROSOS
ARTISTA1REGAN
COLOR SOMBRA Y
DETALLES
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.
14
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:
..
T a bl a 13 T1pos de d.1seno . ' d e procesos
- en mgemena
No ERROR O COSTO
TIPO DE DISEÑO
EXACTITUD INGENIERIL
Orden de Magnitud (Comparación con otra
1 ±40% BAJO
planta)
2 Estudio Estimado (Basado en equipos caros) ±25% !
Estimación Preliminar (Estimación de
3 ± 12% !
presupuestos)
Estimación Definitiva (Detalles sin
4 ±6% !
confección de planos)
Estimación Detallada (Planos de ingeniería
con todas las especificaciones requeridas
5 ±3% ALTO
incluyendo la búsqueda del lugar de
construcción)
También es posible determinar los costos relativos para obtener los tipos de diseño
nombrados anteriormente.
..
T a bl a 14 e osto de mgemena
. para preparar tipos de d.1seno
-
Costos del Proyecto (MM
0,5 1,0 5,0 10,0 15,0 20,0
US$)
Rango de Exactitud Costo de la evaluación (M US$)
-3% a + 12% 15 25 45 70 90 120
-20% a +40% 0,9 1,7 2,8 4,5 6,3 8,0
-50% a +70% 0,25 0,35 0,60 0,90 1,30 1,70
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
15
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.
17
CAPITULO 11: 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 fisicos 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 retoma 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
18
capital incluye todos los fondos necesarios para poner en marcha el proyecto. Esto
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.
d = Depreciación ...
...
~~ Utilidad = IV - CTP
(antes de depreciacdión)
,r
Capital de
.. trabajo
Inversión
... Capital fijo
total de
directo
capital
Pago de
prestamos
.. Capital fijo
indirecto
Otras
inversiones Dividendos accionistas
19
2.1.1.1 INVERSIÓN TOTAL DE CAPITAL (IT)
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
20
Incluyen los costos directamente relacionados con el proceso pero construidos en una
localización separada de los equipos principales del proceso.
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.
21
1) Construcciones temporales: Construcción, operaciones y mantenimiento de
construcciones temporales, oficinas, camiones, áreas de estacionamiento, líneas
férreas, comunicaciones, cercas, etc.
2) Equipos y herramientas de construcción.
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.
1O) 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.
• 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 lndex, 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.
23
DESGLOSAMIENTO DE LA INVERISÓN DE CAPITAL FIJO
ICF=CD+CI
{
• Costo de compra de equipos
Costos dentro de • Costos de instalación de equipos
límites de la planta • Costo instrumentación y control
(CDLP) • Costo cañerías
• Costo equipos y material eléctrico.
• Edificios de procesos
Costos directos Edificios • Edificios auxiliares
(CD) { ~ ;~,l!e~es ~e tn~~t~nimie~t<
Costos fuera de
Arreglo de áreas
límites de la planta
(CFLP) • Servicios
• Facilidades
Servicios de
• Equipos no de proceso
• n1.,tt·;}mf'1nn "
Inversión Terreno
de capital
fiio (lCF)
• Construcciones temporales
• Equipos y herramientas de
Costos de construcción
Costos
indirectos (1 C) construcción • Supervisión de construcción
• Caseta personal y de guardias
• SP.mlrirhil v ~P.rv1r.1o mP.il1r.o
Costos contratistas
Imprevistos
24
Marshal and Swift Nelson- Chemical
Enginnering News Record
Installed -equipment Farrar Engineering
Construction lndex
Año lndex lndex lndex
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
Ca acidad2
Costo2 = Costo1 • ( 'P .
Capaczdad1 )
Dónde:
h = Inversión en dólares de una planta de capacidad 2
11 = 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
25
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.
Existen varios índices que se publican en forma periódica. A continuación se describen tres
de los índices más comunes.
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%
26
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
eqmpo.
Ejemplo:
Purchased _ Cost,US$ = ( ~:OS)101,3·A 0 •65 XFc)
Dónde:
Fe= Factor de corrección para materiales, presión, etc.
27
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
CI =(Costo_ instalado_ del_ equipo_ de_ acero_ al_ carbono)+ (Costo_ incremental_ de_ materiales)
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.
28
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=O.llCF
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+Cl
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:
29
l. 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
CJ = 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.
30
IT = 1.3/CF = 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.
31
l. Costos Directos de Producción (CDP). Cerca del 60% del CTP
a) Materias Primas (entre 1O- 50 del CTP)
b) Servicios (entre 1O- 20% del CTP)
e) Mantenimiento y Reparación (entre 2-10% del ICF)
d) Suministros de Operación (entre 1O - 20% Costos de Mantenimiento y Reparación o
0,5- 1,6% de ICF)
e) Mano de Obra (entre 1O- 20% del CTP)
f) Supervisión Directa (entre 10-25% de la mano de obra)
g) Gastos de Laboratorio (entre 1O- 20 de la mano de obra)
h) Patentes y Derechos (entre 10-20% del CTP)
l. 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 deiS% del CTP).
CTP= CFA + GG
32
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).
eDP = eMAP +es+ 0,04 JeF + 0,006 JeF + 1.35 eMO + 0.03 eTP
eDP = eMAP + es+ es+ 0,046 JeF + 1.35 eMO + 0.03 eTP
Los costos fijos, CF, incluyen impuestos, IMP, seguros SEG, arriendos, ARRe intereses,
INT.
33
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 = O; de forma similar, supondremos que no se arrienda el
terreno, entonces ARR = Oy luego el:
CF= 0.03ICF
Cuando combinamos todas estas expresiones obtenemos una expresión para el costo total
de producción.
CTP = [CMAP + CS + 0.046ICF + 1.35 CMO + 0.03 CTP] + 0.03 ICF + [0.72 CMO +
0.024 ICF]+ 0.025IV
34
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 Tumos requeridos, los costos de mano de obra
por tumo 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 tumo), 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.
UNA!= 0.98 IV- [1.03 (CMAP + CS)+ 0.186 CDLP + 2.13 CMO}
35
Para el caso de ampliaciones y si CMO = 0.15 CTP
2.3.3.2 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.
.
Tabl a 2 2. V'd 1 a est'1mad a para d'1versos eqmpos
1 a med'd
EQUIPOS AÑOS EQUIPOS AÑOS
Calderas 23 Hornos rotatorios 22
Edificios de ladrillos y aceros 35 Molinos 12
Compresores 20 Mezcladoras 12
Condensadores 17 Motores 14
Enfriadoras 17 Cañerías 15
Trituradores 12 Bombas 20
Secadores 25 Retortas 22
Hornos eléctricos 20 Cribas 12
Evaporadores 17 Depósitos 20
Filtro prensa 17 Espesadores 5
Hornos de gas 8 transformadores 15
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 ..
36
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.
.. 1 a med"1a de fi'b
T abl a 2 3 V"d .
a ncas comp11et as
Tipo de fábrica Años
Acidos 15
Productos alcalinos 22
Tintura de anilina 20
Nitrógeno atmosférico 15
Carburo 15
Productos de alquitrán 21
Productos de goma 14
Productos de vidrio 14
Electroquímica 17
Producto de oxigeno 18
Farmacéuticas 20
Celulosa 17
Refinería 25
Jabón 20
Cemento 20
Aceites vegetales 18
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 del48% (EE.UU.).
Donde tes el impuesto de las utilidades (en EE.UU. t = 48%, en Perú t = 30%)
Usando los valores anteriores para el caso de una ampliación de planta, se obtiene:
UD!= 0.52 {0.96 IV -1.51 (CMAP + CS)- 0.27 CDLP- 0.181 CDLP}
UD!= 0.50 IV- O. 785 (CMAP + CS)- 0.236 CDLP
37
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 1z 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.
38
Además de la tasa de retomo 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.
TRI = UD! xl 00
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 retomo, 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 retomo 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:
Este criterio también es muy simple de calcular, pero sufre el mismo problema que la tasa
de retomo 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 retomo (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
39
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
FC O
Costo terreno
ICF
(-) L-------------~~-----------------------------------+--~
o 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 de que se entregan. Así, en un período de tiempo menor a 3 años el costo directo
gastado es alrededor de 1O 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 1O 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.
40
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 retomo, 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.
t CAPT
100% 100% ~VR
95%
90% A A
60% Fe
-3 -2 -1 o
10%CD 40%CI: 40%CD 10%CD 1 2 3 4 n
~
114C1
1/4C1
l/4C1 1/4C1
CAPT
~, CPM
±
}=0
[a¡(CD+ CI)fl+i)i +CAPT+CPM}
41
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.
CAPT+VR
FC
j~
•
FC
FC
o t
¡
-3 -2 -1
ICF/4 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.
VR = 0.03 ICF
CPM= 0.10ICF
IT= 1.30ICF
D = 0.10ICF
Basados en simples suposiciones, podemos ver que existe una relación directa entre los
siguientes valores:
FRC := _!_
3
43
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)
44
IV= 1.57 (CMAP + CS) + 1.10 CDLP
Pp = 1.57(CMAP+CS) + l.IOCDLP
p p
Debemos recordar que nuestra primera estimación del potencial económico, basado en los
precios actuales es:
Después de que se completa los balances de materia y energía, podemos evaluar los
suministros y podemos revisar el potencial económico.
n
PEn =IV -CMAP-CS- L(CDLP}j
j=l
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.
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 retomo. Es decir si nos
imponemos un costo cada año de FRC (IT), el proceso nos dará una tasa de retomo
determinada por FRC.
45
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.
A1 = n · D · H
1
2
Área de la base A= !02__
4
46
w=px txAr (1)
Dónde:
w =Peso
p =densidad del material (plancha)
t =Espesor
AT =Área
2
w = pxt.[ TC.D.H + 2.1C.D ]
4
1C.D2 .H
v=--- v = 1C.r 2 .H (5)
4
Despejando la altura: H:
(6)
2 2
4v 1C.D
w=.p.t. TC.D.--2 + - -] 4v
1C.D ] (7)
[ w=.p.t.[ D +-2-
1C.D 2
47
Entonces, el diámetro del estanque que minimiza el peso será:
8w 4v 4v 4.v
-=o =--+TC.D => -=TC.D => D =- 3
(8)
an 2
D 2
D 1[
Para almacenar productos a alta presión como por ejemplo gas licuado de petróleo.
a. Estanque compresora
\·
.. - :\_ __
Figura 3.2. Compresora de GLP
El espesor de la pared debe tomar en cuenta la presión para prevenir la rotura del estanque
(soldadura).
48
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
que la ecuación anterior será:
P.R
t= +t0 (4)
.S.E-0,6.P
)-i.P.D
t= +t0 (5)
.S.E-0,6.P
es =us~g
Los cabezales son más gruesos que de los costados, es por ello su costo es el doble de los
costados.
49
e = Costo de la plancha* Densidad (solido)* Espesor* Area
E
Dónde:
2tr.D.L
A1 =2.tr.r.L => A1 =--2- => A 1 = tr.D.L (8)
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 1r.R 2 , mientras que el área de una
hemisfera es 2 · 1r.R 2 , entonces podemos suponer que el área superficial de un cabezal
elipsoidal será:
2
C=Cs·Ps·[t0 + }iP.D ](tr.D.L+(1,5Xtr.D )] (12)
SE-0,6P
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
50
;r.D2 .L
V=--- (13)
4
Optimizando:
.. !:_=3 (16)
D
De ecuación. (12)
2
e= cs·Ps·[to + YzP.D ][1r.D.L + (I,sX;r.D )]
SE-0,6.P
!:_ = 3[1 +
D
P.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 T; desde una presión inicial P¡
hasta una presión final P1 , podemos usar una compresora de una etapa o una unidad
multicapas.
51
h = w = MRT;
'P 550
(__!5_){(Pr
K -1 P¡
Jk;l -1} (1)
Dónde:
hp ~ 0,00436QP,(K~i* ~ -1} r
hp = Potencia
w =Trabajo
M = Masa gas o peso del gas (lb 1s).
1'¡ =Temperatura ingreso
---D---0--D---0
w~MRT,(K~i){(;, f'l +(;,f'l +( ~ f'l -3} (2)
Para un compresor de tres etapas podemos encontrar las presiones intermedias que
minimizan el trabajo total w.
0w-
--0-~
- c~K) -(~)
Pz -~
(K;I) c-:K)
Pz (3)
8P2
0w - - c~K) -(~)
- - 0 - P2 ~ -P4
(K;I) ~ c-:K) (4)
a~
52
8w =O=
BP
1 (K -1)p-~ _(KK-1)p(3K;l) p2-(
(K-1) K 2
1
K; }
1
2 P, K
1
Luego:
:. ~2 = P2.P4 (8)
Al reemplazar la ecuación (8) en (7)
Tendremos:
ó (9)
(lO)
REGLA HEURISTICA
"Las razones de compresión de cada etapa en un compresor deben ser iguales"
53
3.2.1 MÉTODO APROXIMADO DE DISEÑO DE GASEODUCTOS
F, t
D;
Do
~----- L-----~
a) Balance de la fuerza F;: Esfuerzo sobre la pared del tubo o caño, propio de la presión
de diseño.
(3.1)
(3.2)
De la figura: (3.3)
Dónde:
~iseño = Presiónde 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.
D 0 = diámetro externo de la cañería, Pulg.
54
1.2.Presión de diseño considerando los factores de seguridad
p 2.S.t
diseño = --.F.L.J.T (3.6)
Do
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)
Do
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 ( l20°C, si la temperatura es
superior se obtiene de normas API el factor correspondiente.
2.5. Factor de tensión fluencia: S; de acuerdo al tipo de caño según API estándares para
5LX.
55
Especificación Mínima tensión fluencia (Psia)
API 5LX grado x42 42000
API 5LX grado x46 46000
• •
• •
• •
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
1[(Do2-D¡2) .L.p
=4 => ws = : .L.p.(D0 - D; )(D0 + D;) (3.10)
:. ws = ;r.L.p.t.(D; + t) (3.12)
Dónde:
L=m
t=m
p = tofm3
D=m
56
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.
(3.14)
De ecuación (3 .15)
(3.18)
2
V5 = 7r (2t XD; + 2t + D; ).L = "· t .(2D; + 2t }L
4 4
Luego:
V = _n_A_t.L_.("--D_:_;_+-'-t) (3.19)
S 4
wd = tr.L.p.t.(D; + t)
:. educto = tr.L.p.t.(D; +t}eacero
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.
57
Gs = G(1- y) Flujo total de gas de salida.
Ls = L(1- x) Flujo total de líquido de salida.
G, Y,
-······---····-·--- N
N-1
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .¡. 2
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=Poxr (3.23)
58
Po= Presión de vapor (referencia)
X = Fracción molar.
Y = coeficiente de actividad.
De ecuación (3.21):
y
y= y +1 (3.24)
X
x=-- (3.25)
X+l
De ecuación (3.23):
(3.26)
y yPO X
y=-=--- (3.27)
Y+l P X+l
.po
Si hacemos m = L = la cual es llamada la Ley de Henry; entonces:
p
Yn =mx" (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:
Y e- mxe A N+J - }
=--- (3.29)
Ys-mxe A-l
Dónde:
59
A=~ (3.30)
mG
L++(A-it: =::]}
N+ 1 = ---=--------------=--- (3.31)
LnA
Una gráfica para estas expresiones puede ser dada en función de N vs YL - mX. La
ecuación (3 .29) también puede escribirse de la siguiente forma:
AN+l- A
Ye-
___;: _ Ys
_;;__ = -N+l
--1.
- (3.32)
Y e- mxe A
/=1-~ (3.33)
Ye
Entonces:
J =Ye- Ys (3.34)
Ye
O también:
Ys=(1-J) (3.35)
Ye
y -y
Ye- Ys = AN+I - A ~ e s = AN+I -A (3.36)
Ye -0 AN+I -1 Ye AN+I -1
(3.37)
Dónde:
A= Factor de absorción.
f= Recuperación fracciona!.
60
"Cualquiera de estas variables pueden determinarse a partir de la figura si se conocen las
otras dos".
1.0
O.S N=2
N=l
0.0
1.0 A=~ 2.0
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 so luto)
b) Si A >2 los flujos de líquido aumentan y las columnas son bastantes costosas
e) Si f ~O entonces Ys =Ye, 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 = O , con ello la columna sería grande y la pérdida del so luto
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):
Ln{1+(A-1{ Ye =mxe]}
N+ 1 = _ _::____ 1Ys mxe
--=:____--~ (3.38)
LnA
Hacemos:
61
B = 1+(A -1{ Y e -mxel (3.39)
1Ys -mxe
(3.43)
(3.44)
LnA = Ln(1 +e)=> LnA = Lnl + Lne => LnA = Lne =e=> A= 1 +e=> e= A -1 (3.45)
62
Por lo tanto: LnA =A-l porque A -1 = 1.4-1 = 0.4 (3.46)
(3.47)
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:
N= Ln{(A-1{~)}
(3.48)
A-1
N= 2.3!og{(A-1{~ )}
(3.49)
A-1
N=
2.3log{(o.4f
\Ys
~J} :::::>N= 2.3log(0.4) + 2.3log(~J
Ys
0.4 0.4 0.4
:. N+ 26lot: J
= (3.50)
63
N~ -3.72- 2.97Ln(l- J)
G, Ys
US$/año
Costo columna 2
G, Ye ··-1·-·--·--c-=.
Figura 3.7. Absorción de un soluto
64
Ye- Ys ]= AN+I -A=> = AN+I -A (3.51)
[
Y e -mxe AN+I -1 f AN+I -1
Dónde:
A=~;/=1-Ys (3.52)
mG Ye
Compromiso Económico:
Si f ~ Ocolumna 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 bÓo, casi
despreciable, entonces tendremos que:
(3.54)
Dónde:
N= Número de platos.
E 0 = Eficiencia (%).
65
Ac = Área de la columna ( pie 3 )
(3.55)
Dónde:
H = Horas de trabajo al año.
C y= Costo del so luto
G =Flujo de gas (aire o vapor).
(3.57)
(3.58)
Dónde:
M G =Peso molecular promedio del gas (aire) Lb/Lb mol.
Pm =Densidad molar del gas (aire) (Lb mol 1 pie3)
(1.5)= Factor para el área de la columna que proviene de suponer un factor de
inundación de 60 %.
66
E 0 =varía entre 40 - 50 % ~ E0 = 0.40
H = 8160~(340dias)puede.ser.8400(350días)
ano
3
3 72 2 97
0.4
1
CTA=-!_[- · - · Ln( - J)]90[G~j+8150·C
{1.5X3600) Y
·G·
Ye
(1- J) (3.62)
90~M~
CTA= Pm [-3.72-2.97Ln(1-J)]+8150·Cy·Ye(1-f) (3.63)
G 3 x 1.5 x 3600x 0.4
90~Ma/
CTAjG = 1 Pm ¡(- 3.72- 2.97Ln(1- J)]] + 8150 . CY. Y e (3.64)
1- j 3 X 1.5 X 3600 X 0.4 1- j
Arreglando:
90~M~m
CTAjG = 8150. e . y + --~~:::____ ¡(-3.72-2.97Ln(1- J)]l
---'--- -"---------'-----'------'-"
1- j y e 3 X 1.5 X 3600 X 0.4 1- j
aCTA/ 90~M a/ 7
2 ·9 =0
/G =8150·C · - /Pm
-8...,...(1--'_'---!"""") Y Ye 3 x 1.5 x 3600 x 0.4 1- f
Despejando f logramos:
90~Ma/ (2.97)
J. =1 - 1 Pm (3.65)
opt 3 X 1.5 X 3600x 0.4x 8150· cy. Y e
Y e= 0.015
Ty = 77° F
67
Cy = 9.86US7,;wl(t8USYrb)
MG =29(aire)
Pm = 0.0802
90~2%.0802 (2.97)
fa ¡ = 1- = 0.999
'P 3 X 1.5 X 3600 X 0.4 X 8150 · 9.86 · 0.015
El costo total de una columna de absorción de platos se puede aproximar con la siguiente
ecuación:
CTA=CAP+COP
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:
68
Dónde:
Fp = l,cuando,P ~ Pmin a= 3.91x 1o-3Psia2(5.68x 1u-4 Kpé)
Fp = l.O+ .. a(P- Pmin),cuandqP~ Pmin Pmin = 65 Psia (450 Kpa )
Ke = 107US$/ . 1.886 (1 006US$/ 1 886)
/ pze /m·
H = 2.3!!__ (3.72)
Ea
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:
0 377
Ea = · a. 2a9 (30- 60% )del.gráfico (3.73)
(m·~: p, J
Dónde:
Diámetro = D = f(flujo.de.vapor)
69
El flujo molar de vapor o aire= a(moy,;ora) (3.75)
Si el flujo volumétrico de vapor o aire es= Q(flX), entonces el área lo hallamos de:
(3.77)
0.88AT =%m V
(3.78)
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.
Balance de fuerzas sobre una gota de líquido (suponiendo que el arrastre comienza cuando
la gota está estacionaria en la corriente de vapor).
(3.79)
70
Para flujos turbulentos (partículas grandes) el coeficiente de arrastre es constante, en estos
casos CD = 0.44 ; reemplazando esta condición en ecuación (3. 79) y luego simplificando
tendremos:
4
- 1r • R
3 (
·g p - p
) g.-
= 0.44-1 p ·V 2 · ;r · R 2 => -8 · - R (pL - po)
· =V 2
3 L G 2 G 3 0.44 Pa
=>V= 8·g·R~-p
L G (3.80)
3(0.44) Pa
Podemos hacer:
8·g·R
(3.81)
Ko = 3(0.44)
(3.82).
Para sistemas que operan cerca del punto crítico (Pa<«PL ), por lo que la ecuación (3.82)
llega a ser:
(3.83)
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:
F= t.s(pie)( L~m )X
3
s pze
F= L~m )Yz
54oo(pie)( pze 3
(3.85)
h
F ~ 6590(: )( ~ r (3.86)
71
Fair a graficado la ecuación (3.80) con respecto a ~(PG )Yz para encontrar la inundación
G PL
límite en platos perforados y burbujas, para factores de absorción A=1.4 y ha encontrado
que a menudo se encuentra en el rango de 0.01 y 0.2 (no varía demasiado) ~(PG )Yz para
G PL
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 PL ~ 58Llfl3 y a= 2odinay:m ).
De ecuación (3.84)
(3.87)
(3.88)
p
p = - - ..parael.gas (3.89)
m R·T
(3.90)
A.4) Cálculo del diámetro para el área de una columna circular se tiene:
2
tr·D
Ar = e (3.91)
4
~
p )-Yz
V= MG·-- (3.92)
R·T
(3.93)
1) Primer caso:
= 54oo(pie)( L~m3 JX
3
F R = 10.73 Psia · Pie
y
h pze Lbmol·R
(3.95)
Dónde:
De= Pie
T=R
P = Psia
G = Lbmol/h
2) Segundo caso:
y R = 8.314 KJ
Kgmol·K
Dónde:
Dc=m
T=K
P=KPa
G=Kgmol/h
~]0.533
CAP=(M &S)·Kc.
280
4 (M aP·R-T)
[0.88·n·F 2 (a)XJ.o66 ·(1.15·h·N)0.82(3.3+Fp) (3.97)
73
~]0.533
CAP=(M&S)·K
280
. 4
e [ 0.88 · 1r · F
(Ma·R·T)
P
2 (1.15·h)o.s2(3.3+F )·G0.533N0.82
P
(3.98)
Si hacemos:
1/ ]0.533
M a · R · T 12
K1 =( M & S
280
4
)·Kc ·[ 0.88 · 1r · F
( P
J (1.15·h)
o.s2
(3.3+Fp) (3.99)
(3.101)
Dónde:
Hy = Horas por año de trabajo (8400 h!año)
Cy = Costo del soluto por kg
(3.102)
Dónde:
Cw = Costo del solvente por kg
Si:
L
A=--~L=A·m·G (3.103)
m·G
(3.104)
Si ahora hacemos que:
K 2 =Hy ·Cy ·M a (3.105)
K3 = Hy ·Cw ·Mw (3.106)
(3.107)
74
(3.108)
(3.109)
75
b. Si los valores no dan buenos resultados, se estima nuevo valores y se continúa con el
cálculo iterativo.
Especificar separación 1
componentes claves
Alimentación especifica
Estimar separación
componentes no claves 2
........
Determinar condición de
alimentación 4 Flash adiabático
'+----------- ~~~~~:,~,~~,~~- J 6
Ecuación de
Fenske
-~:::=~: ~-~:::I_: - -- =
[ Calcular razón de reflujo
mínimo
---····--- ····--··-· -·······-··· ----
7
Ecuación de
Underwood
Correlación
Calcular números de etapas 8
- - -- I ---- - -
teóricos de Gilliland
76
ETAPA No 3: DETERMINAR LA PRESION DE COLUMNA Y CONDENSADOR
C. Estabilidad térmica
Si los componentes llegan a su Te 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:
l. Escoger la presión de operación entre:
O- 2860 KPa (28 atm, 45 Psia) Estos dependen de
2. Elegir una temperatura mínima factores económicos
sooc (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)
77
La pérdida de carga ~p se puede calcular conociendo N etapas.
Para columnas que operan a P ~ que la atmósfera.
~p =N* 0.1 Psi/Etapa
Para columna que operan a P de vacio
~p =N* 0.05 Psi/Etapa
5. Tipo de condensador
dl
bl
Po <2S atm
Calcular Punto de Ts <Te
Estimar Calcular Punto Burbuja
Rocío del Destilado
Usar Presión de Fondo P, Fondo a presión P,
(Po) aTo= SO °C O.K.
Condensador
Parcial
Escoger Refrigerante
para operar Bajar Presión
condensador parcial a Po
Po =28 atm (420 psia)
Para presiones que se encuentren O - 1480 K.Pa (15 atm, 215 Psia) un condensador
total (Refrigerante agua)
Para presiones entre 1480 - 2520 K.Pa (25 atm, 365 Psia) un condensador parcial
(Refrigerante agua)
Para presión ~ 2520 K.Pa utiliza otro refrigerante.
78
B
logP=A--- o
C+t
Dónde:
P 0 =Punto de ebullicion normal (1 atm)
TD = Se aproxima a 50°C
3.- Ecuación de Redder y Radal
In P = 1751 - 6 777
TbF '
Dónde:
TbF = Punto de ebullición de la alimentación
P=MPa
T=K
4.-Ecuación de Clausius-Clapeyron
~
PD =P0 .~X;,Dexp {Mfvio(
·· - 1 J}
1 --
; R Toi TD
Dónde:
Po= Presión de referencia
T oí= Punto ebullición de componente y a Po
T 0 = Punto de burbuja del destilado a T 0
Mfv,i,o =Calor latente de vaporización
Ecuación de Mollepop
D) Punto de Burbuja
e
1 = LX;.K;.(P,T)
¡;¡
79
Una vez determinado la presión de operación y temperatura óptima se vuelve a elaborar las
nuevas condiciones de alimentación y se realiza los cálculos para comparar si ocurre una
separación flash adiabático.
Sistema binario
Z¡
cp=l-K cp =función de evaporación
2
Isotérmico, adiabático
Adiabático
Curva de equilibrio
Línea de alimentación
80
Simbolizando:
Xs Xo
Figura 3.10. Cálculo de número de etapas
KCL
a cr,cP =K- - (3.110)
CP
Si consideramos que:
X CL,N+I J;?: 1 y [ X CP,I J ;?: 1 (3.112)
[
XcL,I XCP,N+I
a:mm = [ X '
. CL N+l
X CL,I
]( X
CP 1
X CP,N+I
' JVolatilidad relativa
• • •
(3.113)
N . = X CL,I X CP,N+l
Ecuación de F enske (3.114)
mm log[aCL,CP ]p
También se puede escribir como:
N. = logS (3.115)
mm l ogaP
81
Según la relación de Antoine
po(T ) .
(a AB )tope -_ p: (TD)
A D
PJ(T 8 )
(a AB )fondo = p: (TB)
Entonces:
[aCL,cP L= · (aCL)
CP T
(aCL)
CP B
(3.116)
(3.117)
Jo~(~::)(!:;)]
N min = -~lo=-g-(:-a_c_Yc_P-,-)-P--= Ecuación de Fenske (3.119)
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
(3.120)
Dónde:
1'¡ =Punto de ebullición de claves
Se despejad o b para hallar los componentes no claves a partir de las ecuaciones (3.114
hasta 3.120)
82
b.
1
= fCL (3.121)
b ~ CL,R )Nm;n
1+(dRJ{a p
R
d. = J CL
1' (d¡;;J(
R
aCL,R
)Nm;n
p
(3.122)
J(a CL,R )Nm;n
1
1+ ( db R p
R
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:
(3.123)
Y= a.X (3.124)
1+(a -1}X
Y=-q-.x- zF (3.125)
q-1 q-1
Underwood (3.126)
83
Underwood (3.127)
(3.128)
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á
(3.130)
(3.131)
o
...¿., aCL .XCL F
L.,¡--------''- =q - 1 (3.132)
CL=J aCL-e
(3.133)
Dónde:
q = condición térmica de la alimentación
r =componente clave pesado
84
La ecuación (3.131) y (3.132) se resuelve para hallar e entre los valores de a para los
claves ligeros y pesados.
Cuando solo los componentes claves se distribuyen la ecuación (3.133) se resuelve para
hallar e (que debe Satisfacer aCL,CP ) e) 1) y este Valor Se reemplaza en la siguiente
expresión para hallar Rm
n
1a X )
R . + 1 = "" ~ CP D,CP (3.134)
mm ~
cP;J acP-
e
Ü (R ) .
a mm
+1 =""~ (acL,y
(a
t).X-e CL,D
(3.135)
CL,y a
Infinitas
CTA
N soluciones
US$/año
1 i/
/ ¡
''
'' ''
l CAP
Ropt-+R
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.
85
Como el valor óptimo varía muy poco se usa como regla heurística la razón de reflujo de
operación óptima para una columna de destilación:
Rm
=
_B._ 1,2
R
--=1,2 (3.138)
Rmin
Dónde:
z.'· F = Alimentación
q =Factor relacionado con la condición térmica de la alimentación.
R =Reflujo
a =Volatilidad (cte.)
Nm = f(XcL,D•XcL,B•a) (3.139)
(3.140)
Debe cumplirse que:
(3.141)
Esto tiene cierta justificación para mezclas binarias, pero para mezclas multicomponentes,
existe incertidumbre de la influencia de los componentes no claves.
86
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:
l. O
Y=N-Nm Y= N-Nm
N+l N+!
0.1
0,0
0.01 0.1 l. O 0.0 l. O
R-R X= R-Rm
X=--m
R+l R+!
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 N m
1{[Xnf ][(1-XJ/ ]}
N ~ 2N = 2. ln S = 2
m n /(1 - X D) 1X w (3.142)
lna lna
(1 ) { do { )~r dcP }
-cpF F LF,CL -\aCL,CP J Lf,CP
D.Rm =LR =Loo= ( )
aCLI -1
ICP F
1
R = 1 y R = -1
m (a-1)XF m (a -1).YF
R = 1,20Rm
87
Utilizando el concepto de eficiencia y el gráfico de O'conell.
_ NTeórico
E o-
N Reales
Ea
Parte lineal
a ·f.Jp
Figura 3.13. Correlación para determinar la eficiencia de columnas
0.4983 0.5
Del gráfico. Eo = (a.f.lF )o.2s2
.
=(a.f.1F )Y. (3.143)
Para valores de a entre 1.10-3.0; Eo varia entre 63% y 51%, por lo tanto para presiones
moderadas podemos usar:
(3.145)
Entonces una estimación rápida es:
lns) lns
N Real= 2( 2-- ~ NReal =4-- Fenske (3.148)
lna lna
0 206
2
NR
Ns
= {[ZcP,F J[XcL,B J
ZcL,F XcP,D
(B)}
D
. (3.149)
Dónde:
NR =Números de platos arriba del plato de alimentación
Ns =Número de platos en la sección de agotamiento
Dónde:
K;= 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 Ri 2Nmin (3.153)
2) Según Fenske:
(3.154)
o
l[[X J[ CL,D X CP,B ]]
N. =~= n x;;; x;;;; (3.155)
mm 1naP 1nap
89
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:
(3.160)
R . = 1 xCL,D
(3.161)
mm (
aF -1 ) ' XcL,F
Reemplazando (3 .160) en (3 .159)
X
D ~ ____!!:L .F (3.163)
XCLD
V= [ 1,2 . X CL,D
a F - 1 X cL,F
+ 1](xCL,F
X cL,D
.FJ
12 X X X ] 12 X
V = ' ._____!!::!!____ • ____!!:L F + ____!!:L .F = ' .F + _____!!:_L .F
[ aF -1 XcLF XcLv XcLv aF -1 XCLD
' ' ' '
90
V =[ 1,2 + X eL,F ] • F (3.164)
aF -1 XCLn
Dónde:
(3.166)
Dónde:
Ce= Costo de la columna
CepN= Costo del condensador
CR=Costo del hervidor
CeA= Costo anual del agua de enfriamiento
Csr Costo total del vapor
FRC=1/3 año- 1 que implica i=15% y vida útil entre 11-15 años
~]0.533
C =(M
e
&S)K
280
4
e [ 0.88.1r.F
(Mc.R.T)
P
2
Go.s33(1. 15 h)o.s2 No.s2(3 .3 0+F)
R P
(3.167)
91
Sistema de
Separación
h·N
H=--+H0 (3.168)
Eo
Dónde:
En la práctica:
h 1 =5-10fl h 2 =5-10fl
Expresando como porcentaje el espacio adicional (Ho = 15%) para evitar el arrastre de
vapor y para el líquido del fondo, la ecuación (3 .168) llega a ser:
H = 1,15h!!_ (3.169)
Eo
_ 4lns
Fenske: NR - (3.171)
lna
Dónde:
s =Factor de separación entre componentes claves
a = Volatilidad relativa
92
Reemplazando (3.171) en (3.170)
41
H = 1,15h( ns) (3.172)
lna
De=
[ 4
0,88.;rr.F
(Ma.R.T)
P
x]x .GX
2 (3.174)
(3.175)
t
h m ~oU
Límites de diseño
H H (175ft = 53m
( 20-30 (Diámetro menor a 0.46 m dificil de construir)
D D ) 1,5 ft = 0.46m
4 ]0.533
Si hacemos: KeL = Ke [ .(1.15h)o.s2 (3.177)
0,88.;rr.F
0,2665
C = M &S K Ma.R.T (3.30+F }No.s2.Go,533 (3.178)
e ( 28Q ) eL [ p ] p
93
a.2.1) Diseño del condensador
Entonces el flujo de calor necesario para enfriar cuando no hay cambio de fase es.
Q=mCp11T
Q = WcCPA(TAS -TAC) (3.180)
Q = h.A.(T8 - 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 = Uc.Ac.!::.Tc (3.181)
Dónde:
(3.182)
(3.183)
Dónde:
94
h¡ y ho = coeficientes que se puede determinar experimentalmente o mediante análisis
dimensional cuando se usa vapor.
4) Ecuación de diseño
(3.184)
(3.185)
(3.186)
Dónde:
T 0 = 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)
ePA= Capacidad calorífica (KJ/kg°C)
MiD = Calor sensible (KJ/kg)
CcoN = (M&S)KHE.A'
280
o6s(
3.30+Fc ) (3.187)
(3.188)
Dónde:
KHE= Costo de compra condensador= 475(us$/ 1,30 )= 101,3
/m
~~~o
pze'
Fe= Factor de corrección por presión dado por la siguiente tabla
95
a.2.4) costo anual del agua de enfriamiento
(3.190)
Dónde:
CA= Costo del agua US$/Kg
H =Operación h/año
(3.192)
(3.194)
(3.196)
96
e =(M&S)K ·( M-Iw )o,65(3.30+F \V0,65 (3.197)
R 280 HE U R./':.TR e}
Dónde:
KHE =Costo compra hervidor =475IUS$/ 130
\ /m'
)= 101,3US$/.
/ pze'130
(3.198)
(3.199)
Dónde:
Cs = Costo de vapor US$/Kg
H =Operación h/año
M-I S = 933BTU/lb
(M +(M
02665
&S)K (MG.R.T) ' (3 _30+F \No.sz.vo,533 &S)K .
1 280 eL p p} 280 HE
eTA = - o 65 o 65
3 [ M-Iv], (330+F ).vo,65+(M&S)K ( M-Iw), (330+F J.vo,65
U e I':.Te ' e 280 HE U I':.T
R R
' e
+ eA.H[
ePA rAs
(M-I~ TAe )]v + es.H·[M-I
M-Is
w].v
(3.200)
0,2665
Hacemos: e = M &s KCL MG.R.T .(330 +F) (3.201)
o ( 280 ) 3 ( p ] ' p
(3.202)
97
e2 _
-
H •[ eA .MIn )
es .MIw
+ --=------"- J (3.203)
ePA .(TAS - TAe Mi S
Q=U.A.dTm (3.204)
(3.205)
Dónde:
Q = Flujo de calor
U =Coeficiente global de transferencia de calor
dTm = Gradiente medio logarítmico de temperatura
A = Área del intercambiador de calor
1.1.1. Determinación de U
1.1.1.1. Determinación de U a partir de los coeficientes individuales de película
1 1 A.
U=-+-.-' (3.206)
h¡ h0 A0
a) los coeficientes interno (hi) externo (ho) se pueden determinar en forma experimental
(vapor)
b) mediante análisis dimensional
e) mediante correlaciones empíricas como la de: Dittas-Boelter
Dónde:
K = Conductividad térmica
eP = Calores específicos
98
JI =Viscosidad de los fluidos (carcasa y tubos)
Sistema de transferencia
%2oc BT~pie 2.h.° F KJ:n2.h.K
Vapor condensante 1 liquido en ebullición 1500 250 5000
Vapor condensante 1 fluido líquido 900 150 3000
Vapor condensante 1 gas 120 20 400
Liquido 1 liquido 300 50 1000
Liquido 1 gas 120 20 400
Gas 1 gas 60 10 200
Condensado parcial 180 30 600
Si el proceso estudiado es rentable, entonces se efectúa diseños rigurosos del cálculo de los
coeficientes (U).
C1 = (M&S).K .A o'65(3,29 + Fe )
280 1 (3.208)
o 65
C1 =( M&S
280
)x [ 1•
Q
U.I1Tm J
'
(3,29+Fc) (3.209)
Dónde:
K1 = .
eosto de compra de1mtercam b"1ador = 101,3---:-)"30
US$ US$
= 475~
pw, m,
Fe = Fp (para intercambiador de acero al carbón)
99
Enfriamos
/ Condensamos
Te 1 Subenfriado
/
Tws
T,
a) Cálculo de los calores para las tres secciones de la corriente del proceso del
intercambiador de calor.
(3.210)
Q2 = F.M-Iv (3.211)
(3.212)
(3.213)
(3.214)
Entonces:
(3.215)
lOO
Q3 Wc.CPw(tz -twe) Q3 lz -twe
QT = Wc.CPw(tws -twe) ~ QT = tws -twe (3.217)
4.2 A partir de la ecuación (3.216) y (3.217) determinamos t 1 y tz 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.
=U A (t 1 -tws)-(Tc -t1 )
Q1 (3.218)
¡• ¡· [( )]
ln 1'¡ - tws
Te -t1
(3.219)
(3.220)
(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.
(3.222)
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.
101
T
Punto de
ebullición
~-----=---~
Q.~
~~ Sobrecalentado
Vaporización
Ebullición
102
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
t., P,
F, Ce, t, el
to
RECUPERACIÓN (R) 1
~
To AR
Calentador Ac t2
TAe
TAs F, Ce, T.
~
Agua Condensador
,---..,.
A.
------
T, ------
------
Reactor (R
------
Columna
(4.1)
103
Figura 4.2. Perfil de temperaturas de enfriamiento
(4.2)
(4.3)
• 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:
(4.4)
• Para sistema gas-gas (a temperatura ::::; 1000°F) se debe recuperar un 50%, entonces:
1
(t 0 -t1 )=-(TR -t1 ) = 0,50(TR -t1 ) (4.5)
2
e. Recuperación óptima
Dónde: r = recuperación
(4.6)
(4.7)
104
Dónde:
CsH =Costo de vapor de calentamiento.
Cws = Costo de agua de enfriamiento.
g. Compromiso económico
CE
Entonces Ce Cmínimo
CsH
Cwc
105
t., P,
F, e,, t,
RECUPERACIÓN {R)
Calentador Ac t,
F, C,, T.
W,
TSIAH,
Reactor
T,
Ae Condensador
Columna
De una planta de secado se obtiene un flujo de aire a 25 oc y 1 atm., con un flujo que
contenido 300 kg./h de a cetona con una concentración de 1,5% molar.
Solución
l. Potencial Económico
106
V, aire?
Refrigerante
Acetona-aire
2s·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
(1)
T= B -C
A -lnP 0 (2)
Donde.
T=K B = 2940,46
P=mmHg C= -35,93
A= 16,6513 P = 7,5x10-5 atm = 0.057 mmHg
Reemplazando:
T= 2940,46 _ (- 35 93 ) T = 187 K
16,6513 -ln(0,057) '
107
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 P A= 0,33 atm.
De ecuación (1)
~ ~ = 4400 atm ~
0 33
p
T
= 7,5x1o-s
• PT PT = 4546 Psia
Refrigerante V, aire?
-loo· e -s6·c
Acetona-aire
2s·c, 1 atm
Flash
300 Kg/h -s6·c
CONDENSADOR
-loo· e
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?
-s6·c -1oo·c rl-s6·c
Acetona-aire
2s·c, 1 atm t,
PREENFRIADOR
t,
r i -----------
-----------
Flash
300 Kg/h -s6·c -----------
---------
CONDENSADOR
Aire
2o·c L, acetona?
-s6·c, Aire
108
a. Planificación del sistema de refrigeración
El rango de operación del proceso (-86°C, -1 00°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 oc, 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 oc, lo cual significa 38 oc para el refrigerante (
f).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 oc T = -35 oc
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.
109
Reemplazando los valores en la ecuación (4)
-sG·c
CONDENSADOR
Acetona-aire
2s·c, 1 atm t2 -sG·c
Flash
300 Kg/h
PREENFRIADOR
T=-101•c
Aire P=17 Psia
2o·c L, acetona?
Valvula Etileno
T=-3s·c
P=250 Psia
Agua de
33•c enfriamiento
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 = PrY = l,36atm y se obtiene T = -29°C, Entonces el calor requerido para enfriar la
alimentación hasta el punto de rocío será:
110
_
Qenf -
Kg KJ Kmol
300-.2,12- + 339,5--.0,997
KJ Kg ( ( )lo
.29-- 25- -29 J e
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.
Simplificando:
Formando una tabla para las variables y por iteración podemos determinar los calores.
111
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
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.
El PFD describe la ruta del proceso, mostrando los flujos de material y energía entre
aquellas unidades de proceso que constituyen la planta.
112
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 familiarizars~ 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.
Datos:
Pd =?, psia F= 0.80
S = 8000 Psia L = 0.90
t = 0.354 pulg. J= 1.00
Do= 42 pulg. T = 1.00
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 p = 7,801 10fm¡. El costo del caño en el lugar es: x70 = US$1000/ton; x80 =
US$1 030/ton.
Solución
a) para el acero tipo x70:
113
El espesor del caño para el acero de grado x70 es:
t = P.D0 = (1440X42) =O
60 1
-2.-s.____,F.'-----L (2X7ooooXo,8oXo,9o) ' pu g.
Ejemplo 3:
Solución:
114
Kg Kcal Kcal
Entonces: Qr =14431-x7500-- => Q =108232500--
h Kg T h
3
Q = 126415000 Kg xO 001 ton = 126415 toncarbon x700 m gas 1año
año ' Kg año ton carbon
Q = (2 540)
p
J[ (P.zT Lpz )]
(
T,
_b 1 - 2
0 549
• { D
2 618
,
':f!.____j_
\
(1)
b '
b f 0 o.461
(2)
Dónde:
_MCF
Qb - G = Gravedad específica gas
h
0
Pb =Psia Tb= R
T¡ = R0
P1 = Psia
L=ft Pz = Psia
Datos
Condiciones estándar
Psia .
~ = 101,325Kpa. = 14,7 Psza
6,894733
115
a) distancia Lima-Huancayo por carretera es 31 O Km.
Tr= 520°R
b) presiones de operación: 147 bar
Presión 40 bar
G = 0,6 (promedio-diseño)
e) Para distribución: 60 bar-40 bar
3 3 3
Q6 =10102m gas x35314fl 3 = 356742.ft => Q6 =357MCF
h ' m h h
0 539 0 61 618
D = [(357Xt4,7X520x1017048) · (0,6) .4 ]Yz. =
6
"
2 2 0 539
(2.450 Xst8.4 )[1450 -1160 ] -
2.S.t
pdiseño =D .F.E.T (3)
o
Para un caño API 5LX grado x65, F=0.72, de diámetro Do= 6 => t = 0.126pulg.
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).
116
1.4. Cálculo Del Peso Del Acero:
ws = ;r.L.p.t.(D; + t) (4)
Datos:
L = 310000 m
p = 7,801 ton/ m 3
m
t = 0.126pu lg.0.0254-~ = 0,0032m
pulg
D; = D 0 - 2t =6- 2(0.126) =5,75x0,0254 = 0,14605
25040000US$
Sabiendo que el costo del caño es de: 14 US$/pulg, entonces para el caso anterior
tendremos que:
117
PARTE 11: 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.
118
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.
Compresor
Purga
Reactor
Tolueno
Difenil
119
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.
1 1 1 Purga
Hz, CH4
e
·O
o
e
Ql
..
'ü
~
o
eQl
u
:J
o1- .
:e
..,
VI
e
Ql
ce
.
Ql
.
Ql
"U UJ
"U
.
Ql
e "U e
E e E
:J :J
ou E
:J ou
ou
Alimentación
Hz
Alimentación
Tolueno
1
5.2.6 Flowsheet simplificado
Sistema Purga
recuperación vapor
Sistema
separación fases
Tolueno
Sistema de
reacción
Benceno
Sistema separación
líquido 1----•
Difenil
Reciclo gas
Benceno
Difenil
Sistema
separación
Reciclo de Tolueno
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.
1
Ejemplo: A+ B - - • • Producto
A + Producto•ot~~----.
... Sub Producto
Recuperación Tratamiento
(Costo de recuperación) (Costo de tratamiento)
/
Recuperar
~
Recuperar
y reciclar y remover
2
Reciclo-Gaseoso
Purga o Venteno
Gas
Producto
Alimentación Sistema de
Separación Sub-Producto
Recuperación de Reactivos
Reciclo-Líquidos
Se trabaja con datos disponibles y posteriormente se usan los resultados para justificar y
obtener datos adicionales.
3
• 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 e
Producto deseado Producto no deseado
1 '
''
' Concentración
p
B J--··-····-····-··········
Máxima
L•'"',/
., . -----·¡
/ '
/--"~:
// ~ e"
i ~ / _____________ ///
/ :' Q ~-----{
___-<~
Ü L/_//---~------=-~1_~-~---------------------------~¡-_-_-____________~
Reciclo Tiempo
. .d d
Selectzvz a = Moles de B Pr aducidos
- -
Moles de B Convertidos
• El costo de materia prima está en rango de 35-85 %del costo total del producto.
4
• 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.
5
- 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.
6
5.4.4 Operaciones múltiples en un simple reactor
Producto
o
m
V'l
!:!'.
QJ
a.
..,o
Ali mentación ..
Reactor y separación de ...
Ca talizador .. Calor
... calor
Alimentación Producto
Calor Separador
Reactor
7
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
8
6.1.2 Remover o reciclar un producto reversible
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.
9
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.
10
• 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:
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.
11
CAPITULO VII: ESTRUCTURA DE RECICLO EN UN DIAGRAMA
DE FLUJO
12
7.1.2.1 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 Rl
A-------------~/
acetona
Reciclo acetona
13
Componente Punto de Ebullición oc 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
Acido acético 118 Reactivo reciclan a R2-liquido
Anhídrido acético 139 Producto primario
CO, CH 4 , C,H 4
Ácido acético
Acetona Anhídrido
Acético
Separador
Reciclo acetona
14
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.
15
CH2CH2 +H20 -----+ CH3CH20H
2 CH3CH20H -----+ (CH3CH2)2 + H20 (1)
Supongamos que deseamos producir 783 mollh del azeotropo EtOH-H20, que contenga
85.4% mol EtOH, la alimentación de etileno contiene 4% CH4 y agua pura.
C2H., CH.
Etileno 4%
C2H., CH 4
Reactor
EtOH-H20
H20 pura
Separador Y.,.o= 85.4% mol EtOH
ETER
H20
PEtOH = Yagua*Pagua
PEtOH = 0.854(783) = 669 mol ETOH!h (3)
Separador
Fw.P
Reactor
Fw (1-x}- FwP
16
1.5 La cantidad de agua que sale con el producto: Fw,P
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:
17
• Necesita Portador de Calor
Para reacciones simples toda la alimentación fresca del reactivo limitante usualmente se
convierte en el proceso (solo se pierde pequeñas cantidades)
(7.1)
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.
(7.2)
Solución:
QR = 1.324*106 BTU/h
Remplazando datos:
18
1.324x1 0
6
= [(51.3 + 22.0) + (2.1 + 0.9 )K22.0Xsn- 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.
U = 20 BTU = 400 KJ
. 2 · op
h ·pze h ·m 2 · oc
Para una carga de calor de QR = 1*10 6 BTU/h se debe determinar el área de transferencia.
1x10 6 BT%
A~ BT% h ""A~ 1000pie 2 ~100m 2
20 h . 2 opx50°F
· pze ·
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- 800m2 ). 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*10 6
BTU/h.
19
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.
Isopropanol ~ Acetona +H 2
Ejemplos:
1) Proceso de ácido sulfúrico
2) Síntesis de amonio
Síntesis 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
20
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
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.
21
CAPITULO VIII: SISTEMAS DE SEPARACIÓN
Purga
¿Cuál es la fase?
Alimentación Productos
Sistema
separación
LÍQUIDO
Alimentación Productos
Sistema
separación
líquidos
RECICLO LÍQUIDOS
22
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.
i
35°C Sistema de Vapor Sistema de
Separación Fases conversión vapor
-
Líquido
1
Líquido
r-----------------------------------------------------
~
/
~
1:~ ~
~?$
Alimentación LÍQUIDO .1 Sistema separación 1 p d t ..
-~ líquido ro uc os
~
s;;tema reacc;,
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.
23
VAPOR
Alimentación
3s·c Sistema de
Separación Fases
Reactor
LÍQUIDO LÍQUIDO
~___ _RECICLO
_ _ _LÍQUIDOS
_ _ _ _ _---¡ s·1st ema separac1on
·' Productos
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)
Fz = Vy. +Lx.
j l l
(8.2)
De donde:
Lx. l
= Fz j
-Vy. l
Entonces:
24
F2 -Vy¡
X - ---''--
¡- L (8.3)
(8.4)
F2 -Vy¡
Y 1. = k.1 -----''--
L (8.5)
De ecuación (8.1):
L=F-V (8.6)
Luego:
(8.7)
O También:
(8.8)
Estas expresiones son iguales cuando existe una división o separación perfecta, para
componentes con un valor de k; entre 0.1 y 1O.
25
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:
(8.13)
(8.14)
(8.15)
Entonces:
(8.16)
(8.17)
(8.18)
. Gas/mmJTlmnnnm< .
Reciclo ~-------------------¡ :------------------: Purga
Si la ubicación del sistema de recuperación se realiza según los items 1o 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.
27
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
Sistema separación
Productos
líquido
Se necesita estimar el tamaño y luego el costo de los sistemas de recuperación que sean los
más baratos.
28
Las decisiones que se debe tomar para el sistema de separación de líquidos son:
l. 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 economiCo se pueden evaluar a través de
respiraderos o aberturas desde el proceso y si además son contaminantes éstas se pueden
29
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.
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.
30
A
A A
B B
e e
e '----.B
Directo Indirecto
1 Número de Componentes 1
2 3 4 5 1 6
1 Número de secuencias 1
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.
31
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.
so. 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
Livianos al
final
32
Reactor
Calor Enfriamiento
Alimentación
Reactante Flash
Livianos al
final
e------------- e------------- e -------------
•0 •O ·0
"ü ------------- "ü "ü -------------
..!!! 14----+-_J ..!!! ..!!!
Q)~ ·~
·~
.....
0
------------- Q)
o
-------------
0
"'
Q) -------------
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
33
Para separar una mezcla B y e 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 e 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 e (+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 e
(punto 7) logramos separar en tope una mezcla como la original que es reciclado al punto l.
B+C
Columna
extracción
B+C
e
0
C+S{+B)
B
Figura 8.9. Extracción
B. Destilación extractiva
34
S
B =HN03
e= H20
e s = H2so.
S
C. Destilación azeotrópica
B+e
azeotropo --------- B+C( azeotropo)
B =Etanol e
e= Agua
S= Benceno
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 temario 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.
35
D. Destilación reactiva
,----..B S
B,C = Xileno
a= 1.03
B
e S= Organometálico
B,S = Complejo
a= 30
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
Cristales de B
Líquido e
Figura 8.13. Cristalización
36
separadores (para las corrientes de purga), reactores, recipientes de vaporizacwn
(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.
37
CAPITULO IX: SINTESIS DE REDES DE INTERCAMBIADORES DE CALOR
Los altos costos de energía durante la década de los setenta motivaron el desarrollo de los
procesos con un uso eficiente de energía. El ahorro de servicios ofrece un potencial
interesante para mejorar la economía de un proceso. Se ha demostrado que el costo de
servicios en columnas de destilación representa el componente principal en el costo total
anual de esos procesos de separación.
Considérese, por ejemplo, la siguiente situación. Una corriente caliente debe enfriarse de
180 a 120 °F, mientras que una corriente fría debe calentarse de 100 a 180 °F. Estos
requerimientos de energía pueden satisfacerse mediante un enfriador con agua de
enfriamiento y un calentador con vapor, respectivamente. Si se desea un sistema que
reduzca el consumo de servicios, se puede comenzar por diseñar un intercambiador entre
las dos corrientes. Consideremos por un momento la opción de que ambas corrientes se
procesaran en forma completa; en el extremo frío del intercambiador se tendría un
diferencia de temperaturas entre las corrientes de 20 °F, mientras que en el otro extremo se
tendría el límite de calentar la corriente fría hasta 180 °F con otra que se encuentra a la
misma temperatura. El costo de servicios desaparece, pero el diseño de ese intercambiador
implica un área infinita. La solución que deseamos se encuentra entre estos dos extremos.
38
Lo anterior establece que debe existir una diferencia de temperaturas en el extremo crítico
del intercambiar que proporciona el mejor compromiso entre el costo de servicios y el
costo del equipo. Este concepto, el de establecer una diferencia mínima de temperaturas
entre corrientes que intercambian calor (~Tmin), es esencial para el diseño de cualquier red
se intercambiadores de calor. Es muy común tomar, en forma heurística, un valor de ~T min
de 10 a 20 °F.
Según Rudd et al, ha propuesto el uso de diagramas de contenido de calor como ayuda en
la asignación de los intercambios que generan una red de recuperación de energía. Aunque
el enfoque es heurístico y no garantiza por lo tanto una solución óptica, proporciona una
manera bastante clara de entender la lógica para una buena solución del problema.
El método se ilustra con el problema mostrado en la tabla 9.1 y 9.2. Se tienen dos
corrientes calientes y dos corrientes frías. Si se usan solamente servicios para satisfacer las
necesidades térmicas de las corrientes, el costo de la operación del sistema se calcula como
98,116 $/año, del cual el98 por ciento (96,025 $/año) es debido a servicios (el lector puede
fácilmente corroborar estos cálculos). El incentivo existe para sintetizar una red de
intercambiadores de calor que reduzca el costo de servicios, y por lo tanto el costo anual
del proceso que utiliza esas corrientes.
La figura 9.1 muestra el diagrama de contenido de calor de las corrientes. Lo primero que
debe corroborarse son los límites de temperatura que termodinámicamente se puede lograr
mediante integración de energía. Suponiendo una ~ T mínima de 20 °F y teniendo la
corriente más caliente de 480 °F, lo más que puede calentar cualquiera corriente fría es a
460 °F. La corriente C2 necesita entonces un servicio externo de calentamiento para llegar
a su temperatura final de 500 °F. Por otro lado, teniendo la corriente más fría es a 140 °F
no existe restricción termodinámica para enfriar la corriente h 1 hasta 200 °F. Usando la
segunda regla heurística se propone un intercambio entre h 2 y C2; la corriente C2 puede
pasar de 240 hasta 460 °F, mientras que la temperatura de h2 que sale del intercambiador es
de 353 °F
39
Tabla 9.2. Datos de diseño
Presión de vapor saturado 962.5 lb/pulgL
Temperatura del agua de enfriamiento 100 °F
Temperatura máxima de salida del agua 180 °F
Diferencia mínima de temperatura
Intercambiadores de calor 20 °F
Calentadores 25 °F
Enfriadores 20 °F
Coeficientes globales de transferencia de calor
Intercambiadores de calor 150 BTU/1Yft2 ·oF
Calentadores 200 BTU!h·ft2 ·°F
Enfriadores 150 BTU!h·ft2 ·°F
Paro de equipo 380 h/año
Parámetros de costo de intercambiadores de calor Costo= a Ab; a= 3250; b = 0.6
A en ft 2
Factor de depreciación 0.1
Costo de agua de enfriamiento 5x10- 5 $/lb
Costo de vapor 1x1o-j $/lb
h2
500 480
400 hl
320
300
280
200
200
500 500
400
320
300
200 240
C2
140
100
C1
Figura 9.1. Diagrama de contenido de calor
La corriente C 1 intercambia calor con la fuente más caliente disponible para lograr su
temperatura final de 320 °F. Usando h 2 para este fin, se tiene que C 1 debe entrar al
intercambio a 219 °F, lo cual satisface la necesidad de enfriamiento de h2 • Finalmente, se
hace el intercambio entre C 1 y h 1, con el cual C 1 se calienta de 140 a 219 °F, y h 1 sale del
40
intercambiador a 252 °F. El resto de energía disponible en h 1 debe proponerse en agua de
enfriamiento. Lo asignación de intercambios se muestra en la Figura 9.2. Nótese que ~T
mínima debe comprobarse cuidadosamente en cada intercambio.
41
h2--~
480
200
140
C1
Figura 9.3. Red de intercambiadores de calor
Comenzamos por ajustar las temperaturas de las corrientes calientes respetando a cada una
de ellas el valor del ~T mínimo especificado. Suponiendo un ~T mínimo de 20 °F, se
obtienen los ajustes que se muestra en la tabla 9.5. Nótese que las temperaturas de las
corrientes frías permanecen inalteradas. Con referencia a las temperaturas ajustadas, se
establece el orden decreciente que se muestra en la tabla (los valores duplicados no se
toman en cuenta). A partir de las cinco temperaturas modificadas que se han obtenido, se
definen cuatro intervalos de temperatura y para cada uno de ellos se efectúa un balance de
entalpia. Para cualquier intervalo i el balance de entalpía está dado por:
42
Tabla 9.5. Temperaturas modificadas
Temperatura Temperatura
Corriente Orden
original eF) a_justada eF)
1 250 230 T¡
100 80 Ts
2 180 160 T3
100 80 Duplicada
3 110 110 T4
200 200 T2
4 110 110 Duplicada
230 230 Duplicada
Los balances para cada intervalo se muestran en la tabla 9.6. La idea ahora es formar una
cascada de calor que vaya desde T 1 hasta T5; si suponemos inicialmente que no hay una
fuente de calor externa (Qh = Q 1 = 0), podemos evaluar la cantidad de calor que fluye desde
un nivel de temperatura T 1 hasta un nivel Ti+! mediante el siguiente balance:
º.
1+
1 =Q.1 +Mi.1
43
e) Punto de pliegue: 110 °F (temperaturas modificadas)
T1
o ----42=Qn
T2
1.5 - - - - 4 3 . 5
T3
- - - - -36.5 - - - - 5.5
T4
-42 - - - - 0 Punto de
pliegue
T5
- - - - 11.7 - - - - 53.7=Qc
Por lo tanto, con respecto a las temperaturas originales, el punto de pliegue para las
corrientes calientes es de 130 °F, y para las corrientes frías es de 11 O °F
El objeto ahora es el de diseñar una red se intercambiadores de calor que cumpla con los
requerimientos mínimos de servicios que han estimado.
El punto del pliegue divide a la red en dos zonas, y los siguientes criterios son esenciales
para obtener la red que se busca:
a) No transferir calor a través del punto de pliegue
Cualquier cantidad de calor que se transfiera debe satisfacer un balance global de
energía, lo cual implica que el consumo de servicios va a exceder a los requerimientos
mínimos, tanto de calentamiento como de enfriamiento, en una cantidad equivalente a
la que se transfiere a través del punto de pliegue
b) No usar calentamiento abajo del punto de pliegue
Cualquier cantidad que se use en este nivel debe rechazarse eventualmente a una
fuente de enfriamiento externa, alterando así el consumo de servicios con respecto a
los requerimientos mínimos
e) No usar enfriamiento arriba del punto de pliegue
Cualquier cantidad de enfriamiento arriba del punto de pliegue debe compensarse con
una cantidad equivalente de calentamiento
Para aplicar el método conviene usar la representación de rejilla que se muestra en la figura
9.5 para indicar que intercambio de calor entre dos corrientes. El número mínimo de
unidades que se requieren para cada lado de la red puede estimarse mediante:
(9.1)
44
h
h----j
=> 1------C
Nótese que no se incluye el servicio de enfriamiento en esta zona. En la zona de abajo del
punto de pliegue o de enfriamiento, el número mínimo de unidades es de 2. A continuación
procederemos con el diseño de la red.
WCp
230 0.9
( 1100
Figura 9.6. Corrientes del problema de diseño ubicadas con respecto al punto de pliegue
45
intercambiador a contracorriente. Tenemos dos posibilidades para comenzar el diseño de la
red, combinar 1 con 4 y 2 con 3, o combinar 1 con 3 y 2 con 4. Para elegir el intercambio
adecuado se debe cumplir el siguiente criterio
Resulta claro entonces que los intercambios iniciales deben ser entre 1 y 3 y entre 2 y 4.
Haciendo estos intercambios, se satisfacen los requerimientos energéticos de la corriente 2
y de la corriente 3. Intercambiando ahora la corriente 1 con la corriente 4 (este intercambio
es ahora factible debido a que la diferencia de temperatura del lado frío del intercambiador
es mayor que ~T), se obtienes que la corriente 1 satisface sus requerimientos energéticos,
mientras que la corriente 4 se calienta hasta 184 op (su temperatura final debe ser de 230
°F). Como ya no existen corrientes calientes disponibles, ajustamos la temperatura de la
corriente 4 hasta su nivel final mediante el uso de servicios de calentamiento. La cantidad
de calentamiento requerida es de 42x 104 lo cual coincide con la cantidad mínima que se
había estimado al hacer los balances de energía entre los intervalos de temperatura
modificados. La red resultante se muestra en la figura 9.7. El número de unidades que se
requieren en la red es de 6, lo cual satisface también el número mínimo requerido.
100
2 1--18_0_ _-+---t
200
110 3
110 4
Figura 9.7. Red de corrientes del problema de diseño ubicadas con respecto al punto de
pliegue
46
En contraste al problema que acabamos de mostrar, donde dos corrientes calientes
intercambiaron calor con dos corrientes frías por encima del punto de pliegue, existen
casos en los cuales los intercambios de calor entre las corrientes originales del proceso
parecen no garantizar el objetivo de generar una red con el consumo mínimo de servicio el
siguiente ejemplo ilustra este situación
EJEMPL09.1
Se requieren usar dos corrientes calientes y una corriente fría para diseñar una red de
intercambiadores de calor. La siguiente tabla proporciona los daos del problema.
Considerando un valor de ~ Tmín de 10°C, use el método del punto de pliegue para integrar
la máxima cantidad de calor entre las corrientes.
SOLUCIÓN
Desarrollando una cascada de calor al valor proporcionado de ~ Tmín se obtiene que la carga
mínima de calentamiento es de 60 kcallh, la carga mínima de enfriamiento es de 350
kcal/h, y que el punto de pliegue es de 80 °C para las corrientes calientes y 70 °C para las
corrientes frías. El siguiente diagrama proporciona la estructura de las corrientes con
respecto al punto de pliegue.
[8_!1 200 80 i
~¡
80 ¡
~
180
..:
1
80 30
.. 170 70~
47
toda la energía de las corrientes calientes es integrable. Esto implica que
termodinámicamente h 1 y h 2 son integrables en esta sección de la red de manera completa.
200 80
H•~------------~
180
H~~---------------4---4
WCp=7
WCp=3
Notar que si el número de corrientes frías es mayor que el número de corrientes calientes,
no necesariamente se necesita dividir corrientes por que la energía de las corrientes
calientes es integrable en forma completa arriba del punto de pliegue.
Similarmente. Puede establecer que:
48
CAPITULO X: PREDICCIÓN DE ÁREAS EN REDES DE INTERCAMBIADORES
DE CALOR
Los principios del método del punto de pliegue tiene la virtud de lograr, mediante la
cascada de calor, una predicción confiable del consumo mínimo de energía. Los
requerimientos de áreas para los problemas mostrados en el capítulo anterior se obtienen
después de que la red ha sido diseñada. Un complemento importante a los requerimientos
de energía lo constituye la predicción del área de transferencia de calor que tendrá el
diseño de una red. En este capítulo mostramos como establecer ese pronóstico de áreas.
Comenzamos por mostrar algunos diagramas auxiliares basados en las propiedades de las
corrientes.
En seguida, se acumulan los ~H's para cada tipo de corriente para preparar los datos a
graficar. Las corrientes calientes comienzan con una entalpía acumulada de O a la
temperatura más baja, mientras que a las corrientes frías se les agrega la cantidad mínima
de calentamiento a todos los valores de entalpía, lo que ocasiona que la curva fría se
ubique por debajo de la caliente a una distancia vertical igual a ~ T mín en algún punto del
diagrama (el cual representa el punto de pliegue). Los datos a graficar quedan como sigue.
La figura 10.1 muestra las curvas compuestas que se obtienen para este problema; notar
que el acercamiento máximo permisible entre las dos curvas está dado por el punto en el
cual la distancia vertical es exactamente igual a ~ Tmín lo cual identifica el punto de pliegue.
La región comprendida entre ambas curvas representa la máxima cantidad de energía
integrable termodinámicamente. El excedente de la curva compuesta caliente, para el cual
49
no hay una contraparte de la curva compuesta fría, debe de procesarse con servicios
externos de enfriamiento, lo cual presenta la cantidad mínima de enfriamiento.
Similarmente, la parte de las corrientes frías que excede horizontalmente a la terminación
de la curva compuesta caliente indica que no existen corrientes de proceso suficientemente
calientes para integrarse con esta parte de la curva compuesta fría, y por lo tanto debe
tratarse con servicios externos, reflejando así la cantidad mínima de calentamiento.
Cantidad mínima de
calentamiento= 42 Btu/h
200
re
.....
....re;;;¡ Corrientes frías
.....
(1)
a.
E
(1)
f- lOo--
~ Punto de pliegue
Cantidad mínima de enfriamiento = 53.7
0 ~------~----B~tu~/_h,_------,-------,-------,-----~
o lOO 200 300
Entalpía
Figura 10.1. Curvas compuestas
Para interpretar el diagrama, leemos que el punto máximo de acercamiento vertical entre
las dos curvas, equivalente a ~ Tmín (20 °F), representa el punto de pliegue igual a 11 O °F
para las corrientes frías y de 130 °F para las corrientes calientes, la distancia horizontal
entre los extremos de las curvas en 1 parte baja de temperatura representa cantidad mínima
de enfriamiento de 53.7x 104 Btulh, y la distancia horizontal entre los dos extremos de las
50
curvas en la parte alta de temperatura representa la cantidad mínima de calentamiento de
42x104 Btulh.
Las curvas compuestas que muestra la figura 10.1 contienen únicamente las corrientes de
proceso que se analizan para integrar energía entre ellas. El diagrama puede
complementarse al incluir los niveles de temperatura de los servicios extremos disponibles.
Supongamos por ejemplo para el problema en cuestión se dispone de vapor de
calentamiento a 280 °F y de agua de enfriamiento a una temperatura promedio de 60 °F
(que puede representar por ejemplo una temperatura de entrada de 50 °F con una
temperatura de salida de 70 °F). La curva compuesta balanceada que se obtiene se muestra
en la figura 10.2. Los servicios sierran el diagrama completamente, de tal forma que toda la
parte de las corrientes calientes tiene una contraparte fría para ceder calor.
Vapor
200
~
::::¡
...... Corrientes frías
....ra
Q)
a.
E
Q)
1-
Entalpía
Figura 10.2. Curvas compuestas balanceadas
51
.. . tes frí as
T a bl a 10 2 Dat os para grafi1car d e 1as comen
T H Original H Modificada
230 o 42
200 1.5 43.5
160 -36.5 5.5
110 -42 o
80 11.7 53.7
300--~------------------------------------~
200--
ro
....
:::1
.....
ro
....
Q)
c.
E
Q)
1-
0--~-----,------,-----,------,-----,----~
o 10 20 30 40 50 60
H
Figura 10.3. Gran curvas compuesta
Un aspecto adicional al del consumo de servicio para la estimación del costo anual de una
red es el conocimiento de las áreas del equipo necesario. El consumo mínimo de servicio
se puede presidir con precisión mediante el concepto de cascada de calor, representado
numéricamente o a través de algún diagrama de curvas compuestas. El área total de
transferencia o a través de algún diagrama de curvas compuestas. El área total de
transferencia de calor también puede estimarse antes del diseño de la red haciendo uso de
las curvas compuestas. Esta característica reviste especial valor, porque al estimar áreas y
consumo de energía antes del diseño puede tenerse el pronóstico del costo de la red antes
de conocer su estructura.
52
contra corriente. La ecuación de diseño para intercambiadores de calor se aplica de manera
directa
(1 0.1)
Si se supone intercambio verticales, lo cual implica que una sección caliente intercambia
calor directamente con su contra parte fría localizada verticalmente debajo de ella, se
espera un uso eficiente de temperatura, tal que el área del equipo tienda a un valor mínimo.
El único artificio necesario es aplicar esta ecuación de diseño en cada sección k dada por el
cambio de pendiente de cualquiera de las dos curvas compuestas, de tal forma que el área
predicha para la red estaría dada por.
(10.2)
Esta ecuación, conocida también como fórmula Bath, se aplica sol para el caso en que el
coeficiente global de transferencia de calor es el mismo para todos los intercambios. El
valor de coeficiente global depende de los coeficientes de película de los fluidos
involucrados en el proceso de intercambio de calor:
1 1 1
-=--+-- (10.3)
u hcal h¡río
(10.4)
De tal manera que la estimación del área total se obtiene al sumar las áreas de cada
intervalo k. de las suma de las áreas definidas por la sección de intercambio comprendida
entre ambas curvas compuestas se obtiene la predicción del área de intercambio necesaria
entre las corrientes del proceso. Para completar la estimación y obtener el área total debe
agregarse la predicción de las áreas necesarias para calentadores y enfriadores. Para lo cual
es necesario usar las curvas compuestas balaceadas y estimar las áreas para calentadores y
enfriadores de la misma manera.
EJEMPLO
Considere el problema de 4 corrientes que se usó en el capítulo 9 para mostrar el método
del punto de pliegue. Complemente la cascada de calor con la predicción del área requerida
por la red usando la formula Bath. Suponga que los coeficientes de películas son
constantes e iguales para todas las corrientes con un valor de 400 Btu/h re
op.
SOLUCIÓN
53
El primer paso consiste en la construcción de las curvas compuestas las cuales se muestran
en la figura 10.1 para las corrientes de procesos y 10.2 para corrientes y servicio. Tomemos
las curvas compuestas balanceadas para predecir el área total requerida incluyendo
calentadores y enfriadores. La figura 10.4 ilustra la forma en que se dividen los bloques de
las curvas compuestas para estimar las áreas de cada uno de ellos. En total, se tienen 5
divisiones al considerar los cambios de pendientes de algunas de las curvas.
Con estos datos y el uso de los valores de los coeficientes de película, la ecuación 10.4
puede aplicarse para este intervalo, para obtener una estimación del área
A2 = 208.5ft2
Similarmente, los demás intervalos se procesan para obtener los valores mostrados en la
tabla siguiente
....ra
....::::¡ra
....
ClJ
c.
E
~
H
Figura 10.4. División de las curvas compuestas para la estimación de áreas requeridas
El área mínima estimada para una red es de 376.6 ft2 . De este valor, 30.1 ft 2 corresponden
a calentadores, 50.1 ft 2 a enfriadores y 296.4 ft 2 a correspondientes de proceso; el área
54
requerida proceso-proceso está ubicada en su totalidad en la selección de arriba del punto
de pliegue.
En el capítulo 9 se obtuvo en red con el mínimo consumo se servicios para este problema.
El lector puede corroborar que el área total requerida por la red diseñada es de 396.4 ft2 , lo
cual está encima de la predicción hecha en un 5%.
Para obtener un costo estimado de la investigación en equipo antes del diseño de la red, se
puede dividir el área correspondiente a cada sección de la red con respecto al punto de
pliegue por el número mínimo de unidades que se requiere en cada sección, y aplicar la
correspondiente ecuación de costos para cada unidad estimada. Esta predicción de área y
su costo permite un procedimiento de pre optimización de valor de L1Tmin que se debe usar
para el diseño de la red, con el fin de obtener el costo anual mínimo.
Debe de enfatizarse que la formula Bath produce una estimación valida del área mínima
necesaria para una red solo para el caso en que los coeficientes de la película sean iguales
para todas las corrientes. De otra manera, su aplicación para la estimación de áreas
requeridas puede ser aceptable, pero el valor obtenido no representa una cota mínima a los
requerimientos de área de intercambio de calor.
55
CAPITULO XI: ANÁLISIS DE REDES DE INTERCAMBIADORES DE CALOR
EXISTENTES
Existen procesos instalados que operan con alguna red de intercambiadores de calor que
fue diseñada bajo esquemas de conocimiento o de escenarios económicos diferentes a los
actuales. En esos casos, no es raro encontrar que el ahorro de energía que se obtiene
representa una fracción del máximo ahorro que puede obtenerse con un diseño nuevo. En
este capítulo se muestra cómo los principios fundamentales de diseño de redes de
intercambiadores de calor a partir del método del punto de pliegue pueden usarse para
modificar redes existentes.
El primer paso para revisar alguna red existente es cuantificar su eficiencia. Una manera de
hacerlo es comparando el consumo de servicios externos de la red con el mínimo consumo
de servicios; el exceso de energía utilizada se buscara integrar, en forma total o parcial,
mediante cambios estratégicos de la red existente.
El consumo mínimo de energía depende del valor de ~T mín· Para fines de diseño, es posible
realizar una pre-optimización de esta variable. Para fines de revisión de redes existentes,
sin embargo, el problema es menos claro. Pudiera pensarse en tres opciones:
1- Formular un problema de diseño y optimizar ~ T mín· usar ese valor para revisar la red.
2- Examinar la red existente y detectar la menor diferencia de temperaturas que se tenga
en alguno de los intercambiadores de calor; usar ese valor como ~Tmín para la revisión
de la red.
3- Usar algún tipo de pronóstico del efecto de ~Tmín en la red al revisarse. Cabe la
posibilidad de usar algún valor heurístico de ~ T mín de acuerdo a la experiencia y
conocimiento del escenario económico de algún diseñador.
Pudiera pensar en que si la red existente es ineficiente, entonces lo más conveniente sería
diseñar una nueva red y dejar a un lado la red anterior. Este procedimiento, aunque tendría
cada vez más incentivo a medida que los precios de los energéticos aumentan, tiene el gran
inconveniente de no aprovechar la inversión que ya se ha hecho en el equipo instalado,
además de requerir de una inversión fija adicional mayor. Una mejor estrategia consiste en
buscar utilizar el equipo existente en mayor parte, y corregir solamente algún punto crítico
de la red. En términos de los conceptos de diseño que se han presentado, la estrategia
consiste en detectar el punto de pliegue de la red existente y corregir los intercambios que
se llevan a cabo a través de ese punto.
Ejemplo
Considérese la red de intercambiadores mostrada en la figura 11.1, la cual es parte de un
proceso industrial que se encuentra en operación. Los datos energéticos y de coeficientes
de película de cada corriente se dan en la tabla 11.1. Usando un valor ~Tmín de 19 °C, haga
un diagnóstico de la red y proponga alguna red modificada que mejore la eficiencia de la
red existente.
56
159
137 77
1 t--------l 4 ¡ - - - - - - - - 1 C1 ¡---~
267
2 1-----------~--~
80
90
¡-----~2~6--j 4
118 S
17,597
Figura 11.1. Red de intercambiadores de calor existente
Datos de costos:
Costo de un intercambiador = 8,600 + 670 A0.83; (A en m 2)
Costo de combustible= 63,360 $/ (MW) (año)
Solución
El primer punto es desarrollar una cascada de calor para determinar los requerimientos
mínimos de servicios y el punto de pliegue. Al hacer este procedimiento, se encuentran las
siguientes propiedades de la red para el valor establecido de ~Tmín de 19 °C:
El siguiente objetivo es corregir el diseño de la red para lograr la mayor fracción posible de
este ahorro potencial. Una vez que se ha desarrollado la cascada de calor, conviene
redibujar la red con respecto al punto de pliegue para detectar con mayor claridad los
intercambios que violan el principio básico de no intercambiar calor atreves de él; estos
57
intercambios son los candidatos naturales a corregir y es de esperarse que sean los
responsables, en alto grado cuando menos, del exceso de energía que consume la red.
La Figura 11.2 muestra la red existente con respecto al punto de pliegue. Puede verse que
existen cuatro equipos que cruzan este punto, un enfriador (C 2 ) y tres intercambiadores (1,
2 y 3). Estos casos deben corregirse. Al desconectar estos intercambiadores se tiene el
esquema mostrado en la figura 11.3, en donde puede observarse que el intercambiador 4, el
enfriador C 1 y el calentador se conservan en las posiciones que tenían originalmente. Al
igual que para fines de diseño, conviene hacer la revisión de la red en las dos secciones que
se generan con respecto al punto de pliegue. Un cálculo útil es el de los calores que tiene
cada corriente en cada sección de la red; por ejemplo la corriente 2 tiene una cantidad de
calor por encima del punto de pliegue de :
77
~------~4 ~----~(1~------_.
80
90
73
128 118 5
140
Figura 11.2. Red existente con respecto al punto de pliegue
Tabla 11.2. Calores disponibles y requeridos por cada corriente con respecto al punto de
r
p11egue (kW)
Corriente Qarriba Qabaio
1 o 18,737
2 2,203 1,612
3 9,899 3,712
4 o 9,423
5 24,512 4,314
58
Las modificaciones que se pueden proponer no son únicas; lo importante, sin embargo, es
conservar los principios básicos del método del punto de pliegue. A continuación se
describe una posible modificación de la red existente, dividiendo el proceso de revisión en
dos secciones.
77
r-------~C1r--------.
171 90
127
73
118m
140
Figura 11.3. Identificación de intercambiadores mal colocados
l. Dejar el enfriador C2 abajo del punto de pliegue con una carga térmica reducida.
Qc2 = 1,612 kW
_15_9_--j@}-----'8=-=0~
..
2. Prevenir el uso de calentamiento por debajo del punto de pliegue; incorporar por lo
tanto un nuevo intercambiador (A) para la corriente 5 e integrarla con la corriente 3 o
con la corriente l. La corriente 1 tiene suficiente carga térmica. Después del
intercambio, la corriente 5 sale 140 °C. A continuación se muestra esta sección de la
red revisada:
159
140
59
3. Procesar la corriente 3 con la corriente 4 usando el intercambiador disponible. Si
fijamos las temperaturas de entrada y de salida de cada corriente, entonces debe
dividirse la corriente 4. El resultado de este intercambio se muestra en seguida.
159
q =5711
127 26
q =3712
159
127 26
Como resultado de esta revisión, la red modificada por debajo del punto de pliegue queda
como se muestra en la Figura 11.4. La nueva carga térmica de enfriamiento es de 8,683 +
1,612 = 10,295 Kw, lo cual coincide con los requerimientos mínimos establecidos por la
cascada de calor.
l - - - - - - - - - 1 C11-----'-7..:__7_.
8683
~1=59~--+-------+--~c2~----------8_o.
1612
159 90
127
140 118
)-------------------------------15
60
Revisión por encima del punto de pliegue.
Existen dos corrientes calientes y una corriente fría. Se puede dividir la corriente fría y usar
los intercambiadores existentes como base. Ambas corrientes calientes se satisfacen y la
temperatura de salida de la corriente fría después de ambos intercambios resulta ser 202
oc. Lo único que queda es terminar de procesar la corriente fría hasta su temperatura
terminal de 265 °C usando el calentador disponible. La carga térmica de calentamiento se
reduce a 12,354 kW, la mínima cantidad requerida. La Figura 11.5 muestra la parte de la
red modificada.
159
159
140
Red modificada
La red completa que se obtiene después de las modificaciones propuestas consume la
mínima cantidad de servicios, y se muestra en la Figura 11.6. Queda por llevar a cabo el
análisis económico que implica esta propuesta.
115 77
C1
8,683
159 80
C2
1612
343 159 90
3 2 3
127
118 S
61
En base a las áreas adicionales requeridas y al escenario económico proporcionado, el
costo de inversión adicional para este proyecto sería de $0.539 X 106 . Los ahorros extra de
energía que se obtienen equivalen a
Antes de llevar a cabo el ajuste de la red existente, pueden establecerse los objetivos que
deben esperarse mediante el uso de la relación entre áreas y energía. Para fines de diseño,
esa relación tiene la forma típica dada por la Figura 11.7. Cuando el ~ Tmín es pequeño, los
requerimientos mínimos de energía son relativamente bajos a expensas de áreas requeridas
en los intercambiadores de calor relativamente grandes; este compromiso se invierte para
valores de ~ Tmín grandes. El problema es poder definir una curva conceptualmente
equivalente para poder establecer objetivos que guíen la modificación de redes existentes.
~Tmín
L __ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ __ _ . grande
Figura 11.7. Requerimientos mínimos de áreas frente a energía para diseño de redes de
intercambiadores
62
energía determinado. A su vez, estos valores dependen del valor de L1Tmín que se use para
la revisión de la red.
a= ~x ::;;1.0
4r
Energía Ex
Figura 11.8. Ubicación de las características de la red existente
Un procedimiento que puede usarse para pronosticar los resultados del ajuste de la red es
trazar una curva auxiliar con a constante y tomarla como el lugar geométrico del
compromiso esperado entre áreas y energía.
De esta manera, puede hacerse un listado o una nueva curva que relacione los ahorros
potenciales contra la inversión adicional y tomarla como una posible fuente de
establecimiento de L1Tmín· Un parámetro económico que puede usarse como una medida de
respaldo para establecer este valor buscado de L1Tmín es el tiempo de recuperación de
capital. Un inversionista esperaría un compromiso adecuado que le permitiera un ahorro
suficiente de energía sin esperar un largo tiempo para recuperar la inversión hecha. Es de
esperarse que mayores ahorros de energía impliquen mayores niveles de inversiones extra.
Si se deseara una alternativa que tuviera un año de recuperación de capital como
restricción, el punto de la curva que refleja este compromiso proporcionaría el valor de
L1Tmín que debería usarse para el ajuste de la red.
63
BffiLIOGRAFIA
2 ... ·Guthrie, Ke11neth, "Capital Costs Estimating", Chemical Engineering, March, 24,
1969, 142-144
6 ... Prausnitz J., Eckert . . C., Orye R., O'Connem· J~ ''Computer Calculations for
. :rv.t:ulticomponent Vapor-LiquidEquilibria", Pretince-Hall, Inc. 1967
64