Untitled
Untitled
Untitled
Director
CLAUDIO ALBERTO MORENO ARIAS
Ingeniero Mecánico
_____________________________ _____________________________
Claudio Alberto Moreno Arias
Firma del director Firma del presidente del jurado
_____________________________
Dany José Cárdenas Romay
Firma del jurado
_____________________________
Felipe Correa Mahecha
Firma del jurado
2
DIRECTIVOS DE LA UNIVERSIDAD
Secretario General
Dr. José Luis Macias Rodríguez
Decano de la Facultad
Dr. Julio César Fuentes Arismendi
Director de Programa
Dra. Nubia Liliana Becerra Ospina
3
DEDICATORIA
Dedicamos este trabajo de grado las personas que han estado con nosotros a lo
largo de este camino; a nuestros padres por brindarnos consejos para salir de los
momentos más complicados y a nuestros familiares que siempre estuvieron ahí para
nosotros.
4
AGRADECIMIENTOS
5
Las Directivas de la Universidad de América, los jurados calificadores y el cuerpo
docente no son responsables por los criterios e ideas expuestas en el presente
documento. Estos corresponden únicamente al autor.
6
CONTENIDO
pág.
RESUMEN 20
INTRODUCCIÓN 21
OBJETIVOS 22
1. GENERALIDADES 23
1.1 Lucarvel 23
1.2 Tilapia 23
1.4.1 Temperatura 28
1.4.3 pH 28
1.4.8 Dureza 29
1.4.9 Alcalinidad 29
7
1.5 Enfermedades comunes en la piscicultura (No infecciosas) 30
2.3.1 Temperatura 48
2.3.2 pH 50
8
2.3.4 Concentración de amonio 52
2.3.6 Biofloc 54
3. MATRIZ DE SELECCIÓN 57
4.1 Funcionamiento 64
4.2 Empaques 64
9
4.3.3 Tamaño del empaque 69
4.3.4 Empaque 70
4.4.1 Presión 73
4.4.2 Temperatura 73
10
4.8 Sistema hidrodinámico 108
CONCLUSIONES 134
BIBLIOGRAFÍA 135
ANEXOS 145
11
LISTA DE FIGURAS
pág.
12
Figura 22. Temperatura ambiente 74
13
Figura 45. Selección del paquete termodinámico 122
14
LISTA DE TABLAS
pág.
15
Tabla 22. Resultados de los balances de materia 115
16
LISTA DE ABREVIATURAS
H Altitud [m]
S Factor de stripping
17
at Área específica de empaque [m-1]
Re Número de Reynolds
We Número de Weber
Fr Número de Froude
18
f(kT/εAB) Función de colisión
FS Factor de seguridad
19
RESUMEN
20
INTRODUCCIÓN
Es por esto por lo que el presente trabajo de grado presenta de manera detallada el
diseño de un equipo para aireación y desgasificación del agua de cultivo; partiendo
desde una revisión bibliográfica de las alternativas existentes, un diagnóstico del agua
de cultivo en la empresa Lucarvel, el dimensionamiento del equipo, la simulación del
proceso y la evaluación de los costos de implementación de la alternativa al proceso.
21
OBJETIVOS
OBJETIVO GENERAL
OBJETIVOS ESPECÍFICOS
22
1. GENERALIDADES
1.1 Lucarvel
1.2 Tilapia
23
Figura 1.
Morfología de la tilapia
Cuando el cultivo se realiza en estanques cerrados con una alta densidad de peces
es muy difícil alimentarlos de forma natural (algas, plancton, invertebrados, etc.) por lo
tanto es necesario complementar la dieta con suplementos alimenticios que contengan
proteínas, lípidos, carbohidratos, vitaminas y minerales.
24
Los lípidos son nutrientes que aportan una elevada cantidad de energía al pez, en la
acuicultura se pueden usar para reemplazar la proteína y disminuir los costos de
alimentación; para una tilapia el porcentaje recomendado varía entre el 7-15% y son los
responsables del transporte de vitaminas [1]
Los carbohidratos no son esenciales en la dieta, pero son los menos costosos,
muchas empresas acuícolas los incluyen en la dieta para disminuir los costos del
proceso. En los peces los carbohidratos son almacenados como glucógeno que puede
utilizarse para satisfacer las demandas energéticas, a pesar de eso, el pez solo puede
extraer 1.6 calorías por cada gramo de carbohidrato, mientras que los mamíferos
extraen 4 calorías en la misma cantidad. [6]
25
de tilapia que fueron valoradas en US $36 millones, representando un aumento del
47% en toneladas y 21% en valor respecto al mismo periodo del año pasado, siendo
Estados Unidos el mayor consumidor con el 96,1% del total de exportaciones de tilapia;
el porcentaje restante se exporta a Perú, Reino Unido y Alemania. [8]
Los altos costos de mantenimiento de los cultivos sumado a las bajas ventas y los
impedimentos de la crisis de salud generaron que muchos productores cerraran sus
negocios al entrar en pérdidas.
26
Es recomendado para productores que no tienen grandes cantidades de dinero para
invertir y/o donde no es fácil conseguir alimentos de buena calidad. Los estanques que
se utilizan pueden llenarse y drenarse a gusto del cultivador y suelen tener un tamaño
de pocos metros cuadrados; la tasa de siembra varía entre 50.000 a 100.000 peces/Ha.
Figura 2.
27
1.4 Variables del proceso
1.4.1 Temperatura
1.4.3 pH
28
1.4.4 Concentración de amonio
El amonio (NH4+) y el amonio no ionizado (NH3) son las sustancias más tóxicas en el
cultivo acuático. El amonio es menos tóxico que el amonio no ionizado, por lo tanto, se
debe asegurar una concentración de amonio no ionizado menor a 1 mg/L. [14]
La toxicidad por los nitratos ocurre cuando la relación nitrato-nitrógeno supera los
400 mg/L. El recambio de agua o la remoción de sólidos del tanque generalmente
reducen la concentración de nitratos, pese a esto, algunos sistemas pueden requerir
procesos de desnitrificación adicionales. [14]
1.4.8 Dureza
1.4.9 Alcalinidad
Esta variable mide la capacidad amortiguadora del pH del agua, esta variable debe
mantener unos valores entre 150 y 250 mg/L. Para lograrlo se añade bicarbonato de
sodio ya que es un compuesto de fácil acceso, económico, completamente soluble y de
29
fácil manejo. Se recomienda utilizar en el cultivo una fuente de agua con alta
alcalinidad para disminuir los costos del cultivo. [14]
30
Los peces pueden sobrevivir a altos niveles de pH, pero no toleran una variación
brusca en sus valores, cuando llueve y los estanques no tienen un control de entrada
de agua apropiado el pH disminuye de manera abrupta provocando una lesión aguda
en las branquias que les provocan disfunción respiratoria y muerte a los peces,
adicional también generan lesiones en la piel, las aletas y la córnea. [11]
Figura 3.
31
La sobresaturación ocurre cuando el agua contiene una mayor cantidad de gases
disueltos de la que puede contener a una presión y temperaturas dadas. A altas
presiones y bajas temperaturas la cantidad de gases que puede almacenar el agua es
mayor, no obstante, estas condiciones en un cultivo de peces pueden ser bastante
perjudiciales para el cultivo. [14]
Este tipo de materiales pueden ser obtenidos a través de los afluentes que llenan el
estanque como los residuos de materia orgánica proveniente de los alimentos, teniendo
en cuenta que la tasa de asimilación del alimento por parte de los peces oscila entre el
70% y 75%. [16]
En los estanques en donde la calidad del agua no es controlada, los desechos de los
peces (amoniaco y nitrito) pueden producir un efecto tóxico en los peces en elevadas
concentraciones.
32
Generalmente no es una situación de alto riesgo, pero en ocasiones se generan
ulceras profundas que conllevan al fallo en el control osmorregulador del pez. Para
mitigar este problema se disminuye al máximo la manipulación mecánica del sistema.
Una torre empacada cuenta con cuerpos sólidos o inertes que reciben el nombre de
empaques; el aire entra por la parte inferior de la torre y se va elevando debido a su
densidad, entrando en contacto con los empaques, por otro lado, el agua del proceso
piscícola entra por la parte superior de la torre y desciende entrando en contacto con el
aire y los empaques. El objetivo de los empaques es generar un constante contacto
entre el agua y el aire, minimizando el grosor de la capa del agua y facilitando la
transferencia de masa para eliminar contaminantes como el amoniaco, el dióxido de
carbono y el sulfuro de hidrógeno. [17]
33
aleatoria. Para que un empaque se considere como bueno debe proporcionar espacios
vacíos y caídas de presión bajas del gas.
Los empaques de forma aleatoria son aquellos que se instalan en la torre sin un
orden en específico, para la elección de estos empaques se tiene en cuenta la
resistencia a la corrosión, el peso y el área superficial. En las torres estructuradas los
empaques se organizan en cuadriculas o en mallas lo que permite una menor caída de
presión, sin embargo, son costosos y no recomendables para una torre de pequeñas
dimensiones.
Figura 4.
Figura 5.
35
Se recomienda el uso de estas torres cuando se requiere oxigenar el agua o cuando
se requiere eliminar un gas de una corriente de agua, con la ventaja de tener un diseño
compacto.
Figura 6.
Torre de bandejas
Las torres de pulverización consisten en una torre con aspersores, estas funcionan
dispersando el agua de cultivo en una torre con aire
Figura 7.
Torre de pulverización
37
Las gotas dispersadas entran en contacto con el aire, se oxigenan y liberan los
gases no deseados, volviendo de nuevo al sistema de cultivo con buena calidad. [18]
Son las torres que menor requerimiento energético tiene y es de fácil instalación;
para llevar a cabo el proceso con una de estas torres es necesario tener en cuenta el
tamaño de la gota, la velocidad, el ángulo de dispersión y la presión de dispersión. Si
no se tienen en cuenta estos aspectos el agua puede pasar por el proceso de ebullición
y salir por la parte superior de la columna.
Actualmente la industria química ofrece soluciones para extraer los gases de una
corriente de agua, sin embargo, estos métodos son bastante específicos, por lo que
para eliminar más de un compuesto de una corriente de agua se requiere una alta
inversión. [21]
Para la eliminación del dióxido de carbono del agua de cultivo puede utilizarse la
absorción química con carbonato de potasio (K2CO3), la cual arroja como resultado la
formación de bicarbonato de potasio (KHCO3) [22]. El bicarbonato de potasio es
utilizado en los cultivos acuícolas para mantener la alcalinidad del agua.
Otras técnicas para la desgasificación del agua incluyen tratamientos con hidracina
(N2H4), aminas o sulfito de sodio (Na2SO3) que cumplen la función de secuestrantes
[23], que como su nombre lo indica su función es secuestrar la molécula del gas
disuelto en el agua; sin embargo, así como elimina el dióxido de carbono disuelto
38
también elimina el oxígeno disuelto en el agua; por lo que este tipo de técnicas puede
ser descartada para el objetivo de este trabajo de grado.
39
2. DIAGNÓSTICO DEL PROCESO DE CULTIVO
Lucarvel maneja un cultivo semi-intensivo de tilapias que cuenta con 3 etapas; cada
una de las etapas se lleva a cabo en un tanque de diferente diámetro. La distribución
de los tanques en la empresa se ve en la Figura 8.
Figura 8.
Nota. Planos de vista superior de la distribución del proceso del cultivo de tilapias
en la empresa Lucarvel.
40
La segunda etapa se desarrolla en tanques de 9 metros de diámetro, en los cuales
alcanza un peso promedio de 140 gramos/pez; en esta etapa se presenta una
sobrevivencia del 98% en un tiempo de 6-8 semanas.
41
orgánica, donde las células microbianas tienen del 2 al 20% y del 30 al 40% de materia
orgánica, conteniendo así el agua entre 50 y 200 micras. [25]
Entre las bacterias que componen el biofloc se encuentran las nitrosomonas y las
nitrobacter. Las bacterias nitrosomonas son las que intervienen en la oxidación del ion
amonio a ion nitrito; este tipo de bacterias se desarrollan en un pH de 6 a 9 y una
temperatura entre 20 y 30°C y usan los electrones de la oxidación del amonio para
producir energía. Las nitrobacter intervienen en la oxidación del ion nitrito al ion nitrato;
su temperatura de crecimiento óptimo para este tipo de bacterias es de 28°C en un
rango de pH entre 7,6 y 7,8. [26]
Figura 9.
Reacciones de la nitrificación
42
proceso de nitrificación, obteniendo iones nitratos contenidos en los flóculos que
posteriormente son consumidos por las tilapias, como se puede observar en la Figura
10.
Figura 10.
Para mantener el agua en los niveles de oxígeno adecuados para el cultivo se utiliza
un soplador de aire conocido como blower. Un blower es una turbina que inyecta aire al
sistema mediante la compresión del aire por medio de rotores; la ventaja de este
equipo es que requiere un mínimo mantenimiento ya que no hay fricción de los
componentes en su interior lo que no genera desgaste, y presenta una alta eficiencia
de oxigenación, sin embargo, no es eficiente desgasificando el agua de cultivo. El
sistema de aireación actualmente permanece encendido las 24 horas del día.
43
Figura 11.
44
Para la empresa es complicado mantener los gastos asociados al alto consumo de
agua y energía, es por esta razón que Lucarvel busca una alternativa de aireación y
desgasificación que disminuya los gastos y sea eficiente para su proceso.
En primer lugar, se realiza una estimación del costo del servicio eléctrico, para ello
se tiene en cuenta el valor por kWh definido por Celsia de Colombia S.A, empresa que
presta el servicio eléctrico en la empresa, además de los consumos de los equipos
involucrados en los procesos y el tiempo de operación en cada estanque.
Ecuación 1.
Costo de energía
𝐶𝑜𝑠𝑡𝑜 𝑒𝑛𝑒𝑟𝑔í𝑎
= 𝐶𝑜𝑛𝑠𝑢𝑚𝑜 𝑒𝑞𝑢𝑖𝑝𝑜 (𝑘𝑊) ∗ 𝑡𝑖𝑒𝑚𝑝𝑜 𝑑𝑒 𝑜𝑝𝑒𝑟𝑎𝑐𝑖ó𝑛 (ℎ) ∗ 30 𝑑í𝑎𝑠
𝐶𝑂𝑃
∗ 𝑇𝑎𝑟𝑖𝑓𝑎 𝑑𝑒 𝑒𝑛𝑒𝑟𝑔í𝑎 ( )
𝑘𝑊ℎ
𝐶𝑂𝑃
𝐶𝑜𝑠𝑡𝑜 𝑒𝑛𝑒𝑟𝑔í𝑎 = 0,75 𝑘𝑊 ∗ 24 ℎ ∗ 30 𝑑í𝑎𝑠 ∗ 529,25
𝑘𝑊ℎ
45
Tabla 1.
Por otro lado, debido a que el recambio de agua se realiza con recursos del lago, no
genera grandes costos en el servicio del acueducto; sin embargo, la empresa tiene
previsto utilizar el lago para expandir su capacidad de producción y realizar un cultivo
semi intensivo con jaulas en el lago. De ser así el recambio de agua se realizaría con
agua proveniente del acueducto, por lo cual se realizará un estimado de costos del
servicio de acueducto para suplir las necesidades del recambio.
Tabla 2.
MUNICIPIO DE PLAMIRA
ACUEDUCTO
Estrato C. Fijo (COP) 0-16 m3 ($/m3) > 16 m3 ($/m3)
Estrato 1 $ 3.533,51 $ 1.173,68 2347,35
Estrato 2 $ 6.643,00 $ 2.206,51 2347,35
Estrato 3 $ 7.067,03 $ 2.347,35 2347,35
46
Tabla 2. (Continuación)
Vale la pena recalcar que el proceso de recambio se realiza con la mitad del
volumen de llenado de cada uno de los estanques; de esta forma se procede a hacer el
cálculo respectivo con las tarifas del servicio de acueducto del año 2021.
Teniendo en cuenta los volúmenes de cada uno de los estanques se realizan los
cálculos respectivos para la estimación de costos del servicio de acueducto, ya que de
seguir utilizando los métodos de aireación y desgasificación actuales y contemplando
los planes futuros de la empresa estos gastos deben ser valorados en el diagnóstico
del proceso.
Ecuación 2.
Costo de agua
47
El modelo de cálculo para determinar el costo de operación de los tanques de 5 m se
puede observar a continuación.
𝐶𝑂𝑃
𝐶𝑜𝑠𝑡𝑜 𝑎𝑔𝑢𝑎 = 9,82 𝑚3 ∗ 5 ∗ 30 𝑑í𝑎𝑠 ∗ 2347
𝑚3
Los resultados de los costos asociados al agua para todos los tanques se muestran
en la Tabla 3.
Tabla 3.
Valor acueducto
$ 2.347
($/m3)
48
2.3.1 Temperatura
Figura 12.
Para conocer el comportamiento de esta variable a lo largo del día se tomaron dos
mediciones diarias en el transcurso del mes de abril del año 2021, una a las ocho de la
mañana y otra a las cuatro de la tarde teniendo en cuenta la jornada laboral del
trabajador.
49
Como se observa en la Figura 13, los resultados de las mediciones arrojaron que en
horas de la tarde la temperatura es mayor posiblemente debido a los efectos de la
radiación solar; mientras que en la mañana las temperaturas son más bajas.
Figura 13.
28
Temperatura (°C)
27
26
25
24
23
22
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
Tiempo (Día)
Mañana Tarde
2.3.2 pH
El kit es de la marca API (Figura 14), consiste en diversos reactivos para medir
distintas propiedades del agua y está diseñado para procesos de acuicultura, la
desventaja de este kit es la cantidad limitada de muestras que se pueden tomar; por lo
tanto, para este parámetro se realizó la medición una vez diaria durante todo el mes de
abril del año 2021.
50
Figura 14.
Figura 15.
pH en el agua de cultivo
pH en el agua de cultivo
8,00
7,50
7,00
6,50
pH
6,00
5,50
5,00
4,50
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
Tiempo (Días)
51
2.3.3 Concentración de oxígeno disuelto
La concentración del oxígeno disuelto fue medida con el mismo instrumento que se
usó para la temperatura, por lo tanto, las mediciones se realizaron la misma cantidad
de veces al día en el mismo horario (Figura 16).
Figura 16.
5
4,5
4
3,5
3
2,5
2
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
Tiempo (Día)
Mañana Tarde
La concentración de amonio fue medida con el mismo kit de medición utilizado para
el pH; se tomaron muestras diarias a lo largo del mes de abril de 2021 (Figura 17).
52
Figura 17
8
Amonio (mg/L)
0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
Tiempo (Días)
Para el control de esta variable la empresa utiliza la técnica biofloc, la cual fue
explicada a principios del presente capítulo.
Las mediciones del dióxido de carbono no eran tenidas en cuenta por la empresa,
por lo tanto, no contaban con un equipo que las realizara. Para determinar el punto de
partida de la concentración del dióxido de carbono en el agua se adquirió un equipo de
titulación de la marca HANNA Instruments.
El equipo fue adquirido por parte de los investigadores y debido a las muestras
limitadas se estableció un periodo de tiempo específico para la toma de muestras.
53
Figura 18.
30,5
30,0
29,5
29,0
28,5
28,0
27,5
27,0
26,5
26,0
25,5
25,0
24,5
24,0
1 2 3 4 5 6 7 8 9 10 11 12
Tiempo (Horas)
2.3.6 Biofloc
Figura 19.
Nota. Biofloc en el
agua de cultivo.
54
Para realizar el muestreo de la técnica biofloc se utiliza una botella plástica en donde
se realiza el proceso de sedimentación de la materia orgánica presente en el agua. Las
muestras se toman diariamente al comienzo de la jornada laboral del trabajador.
Vale la pena aclarar que los valores tenidos en cuenta para hacer la comparación
fueron los obtenidos en el mes de abril de 2021, y se tomaron como rango para tener
en cuenta todos los datos, máximos y mínimos.
Tabla 4.
55
Los parámetros que están cercanos o en los valores de letalidad son la
concentración de oxígeno disuelto, la concentración de amonio y la concentración de
dióxido de carbono; si bien el pH no tiene un óptimo valor, no genera un ambiente letal
para el pez.
56
3. MATRIZ DE SELECCIÓN
Esta matriz se realiza con los criterios en la primera columna y se evalúan en todas
las alternativas de diseño con +1, -1 o 0 dependiendo si el cumplimiento del criterio es
superior, inferior o igual al criterio base. [30]
En el método de las sumas ponderadas todos los criterios deben tener un rango de
valoración numérica; a cada criterio se le asigna un peso que va de mayor a menor
57
dependiendo la importancia del criterio en la selección. Los pesos normalizados deben
sumar 1. [32]
Ecuación 3.
𝑀𝑖 = ∑ 𝑎𝑖𝑗 ∗ 𝑤𝑗
Se disponen los criterios en una jerarquía en donde el más alto es el criterio objetivo
seguido de los otros criterios. Se realiza una comparación por pares de los criterios
para obtener su importancia para después comparar cada criterio con cada alternativa
y obtener la mejor alternativa para el proceso. [33]
58
3.3 Criterios de evaluación
1. Espacio utilizado
2. Consumo energético
3. Mantenimiento
4. Inversión
5. Eficiencia en eliminación de dióxido de carbono
6. Seguridad
7. Tiempo de proceso
8. Instalación y movimiento
9. Equipos necesarios
Vale la pena aclarar que cada uno de los criterios mencionados anteriormente tiene
un peso distinto para la empresa, por lo tanto, con la ayuda de los directivos se
establecieron los pesos para cada uno de los criterios.
En primer lugar, los criterios fueron comparados por pares según la clasificación
presente en la Tabla 5.
Tabla 5.
Valor Descripción
20 Mucho más importante
15 Más importante
10 Igual
5 Menos importante
1 Mucho menos importante
Nota. Valores para la clasificación
de los criterios de evaluación de
alternativas.
59
Con los valores definidos se procede a realizar la comparación de criterios por pares.
Tabla 6.
Número de criterio
1 2 3 4 5 6 7 8 9 SUMA PESO
1 5 15 5 5 10 10 5 5 60 0,085
2 20 15 10 10 10 15 10 10 100 0,142
3 15 5 10 10 10 10 10 5 75 0,106
Número de criterio
4 15 10 10 10 10 15 15 10 95 0,135
5 15 10 15 10 15 15 20 10 110 0,156
6 10 5 10 5 10 10 10 10 70 0,099
7 5 1 5 5 5 10 5 5 41 0,058
8 10 5 10 10 5 10 15 10 75 0,106
9 15 5 10 10 5 10 15 10 80 0,113
TOTAL 706 1,000
Nota. Comparación de criterios por pares para determinar el peso de
los criterios.
Tabla 7.
CRITERIO PESO
Eficiencia de eliminación de CO2 0,156
Consumo energético 0,142
60
Tabla 7 (Continuación)
Inversión 0,135
Equipos necesarios 0,113
Instalación y movimiento 0,106
Mantenimiento 0,106
Seguridad 0,099
Espacio utilizado 0,085
Tiempo de proceso 0,058
Nota. Jerarquía de los criterios de
selección ubicados del más importante al
menos importante.
Una vez se tenga la división de criterios se evalúa cada alternativa con respecto a
cada uno de ellos y se aplica la ecuación del método de sumas ponderadas vista en la
sección 3.1. La alternativa de diseño que mejor valoración obtenga va a ser la escogida
y se procederá a realizar el diseño del equipo.
Tabla 8.
61
Tabla 8 (Continuación)
Tabla 9.
CRITERIO
SUMA
ALTERNATIVA 1 2 3 4 5 6 7 8 9
PONDERADA
Torre empaquetada
4 4 3 3 4 4 3 3 4 3,5949
a contracorriente
62
Tabla 9 (Continuación)
Torre empaquetada
4 4 2 3 4 4 3 3 4 3,4887
de flujo transversal
Torre de bandejas 4 4 4 4 2 4 3 3 4 3,5241
Torre de
4 3 4 3 3 4 4 4 4 3,5680
pulverización
Métodos químicos 2 2 1 3 4 3 3 3 2 2,6034
w 0,085 0,142 0,106 0,135 0,156 0,099 0,058 0,106 0,113
Nota. Matriz de selección por medio del método de sumas ponderadas.
63
4. TORRE EMPACADA A CONTRACORRIENTE
4.1 Funcionamiento
4.2 Empaques
Los empaques son sólidos que rellenan las columnas y tienen como objetivo
aumentar el contacto entre la fase gas y la fase liquida, para obtener una mayor área
de contacto con una mínima resistencia al flujo y de esta forma obtener una mayor
transferencia de masa. [36]
Son aquellos que son arrojados en la torre durante la instalación y por lo tanto
quedan distribuidos de forma aleatoria.
64
Los empaques aleatorios pueden construirse de metal, plástico, cerámica u otros
materiales; todo dependerá del tipo de sustancia que se está manejando y los costos
de los materiales. [37] En la Figura 20 se observan los tipos de empaque al azar.
Figura 20.
65
consisten en placas de metal corrugado o laminas delgadas dispuestas de manera que
obligan a los fluidos a recorrer caminos obstaculizados aumentando el área de
transferencia de masa. [38]
Figura 21.
Empaque estructurado
66
Como segunda medida el empaque debe tener las características deseables del flujo
del gas y el líquido; una propiedad de los empaques es el espacio vacío (ε), entre más
grande permite más flujo. Un empaque debe permitir el paso de grandes cantidades de
flujos a través de la torre, esto genera menor caída de presión y poca acumulación de
fluidos.
Un empaque debe ser químicamente inerte con los fluidos que entran al proceso de
aireación y desgasificación, de no ser así la naturaleza química de las sustancias
involucradas se ve comprometida y afectaría al proceso.
La estructura y resistencia del empaque deben ser elevadas para manejar y soportar
los flujos del proceso y el costo debe ser bajo.
Tabla 10.
Tipo de empaque
Criterio Al azar Estructurado
Costo relativo Bajo Alto
Caída de presión Moderada Baja
Eficiencia de transferencia de masa Alta Muy alta
Capacidad del vapor Alta Alta
Nota. Características de los tipos de empaques.
67
Si bien los dos tipos de empaques cumplen con el mismo objetivo existen
parámetros para determinar cuál de los dos tipos de empaques es el más conveniente.
Los empaques al azar son construidos de diversos materiales; cada material tiene
una aplicación en específico en los procesos que se realizan con empaques (Tabla 11),
así mismo poseen diferentes características que influyen en la transferencia de masa.
Tabla 11.
Material Aplicación
Polipropileno Sistemas acuosos hasta 100 °C
Acero al carbón Sistemas de hidrocarburos no corrosivos
Cerámica Sistemas de absorción de ácidos nítrico y sulfúrico
Acero inoxidable Sistemas con hidrocarburos
PFA Sistemas acuosos corrosivos hasta 150 °C
PVC Sistemas acuosos corrosivos hasta 90°C
Nota. Aplicación de los diferentes materiales en empaques al azar.
Tomado de “HAT International”, “Random Packing”, [En línea].
Disponible: https://cutt.ly/Znz0UB0 [Acceso 24 may, 2021]
68
El sistema tratado en la empresa Lucarvel consiste en un sistema acuoso no
corrosivo manejado a temperatura ambiente, por lo cual, los empaques de polipropileno
son los adecuados para el proceso.
Tabla 12.
69
Como medida general los tamaños de empaque de 25 mm o mayores se utilizan
para un flujo de gas de 0,25 m3/s y los de 50 mm o mayores para un flujo de gas de 1
m3/s. [34]
4.3.4 Empaque
Tabla 13
Empaque Descripción
Es el anillo más antiguo, es el más económico por costo
unitario, pero no tiene alta eficiencia. Está disponible en
Anillos Raschig una amplia variedad de materiales. Es estructuralmente
muy sólido y los tamaños más grandes se apilan a mano.
Produce un empuje lateral en la torre.
70
Tabla 13 (Continuación)
Una vez con los diferentes tipos de empaques definidos se realiza una evaluación
respecto a los usos que se les da en la industria, como se observa en la Tabla 14.
Tabla 14
Empaque Uso
Industrias petrolíferas, química, gas de
Anillos Raschig
hulla y protección medio ambiental.
Columnas de secado, columnas de
absorción, torres de enfriamiento en la
Sillas de Berl
industria química, industria metalúrgica y
la producción de oxígeno.
Destilación, regeneración de solventes y
Sillas Intalox
sistemas de desorción.
71
Tabla 14 (Continuación)
Cuando los anillos Pall son comparados con los anillos Raschig se determina un
incremento de la eficiencia de alrededor 45 % y de hasta por lo menos el 20% cuando
se compara con empaques de tipo silla. [40]
72
Tomando en consideración los datos anteriormente expuestos se seleccionan los
anillos Pall como empaque para el diseño del proceso.
4.4.1 Presión
Ecuación 4.
La empresa se encuentra ubicada a una altitud de 1018 m sobre el nivel del mar; al
aplicar la Ecuación 4 se determina una presión de 89669 Pa que son equivalentes a
0,89 atm.
4.4.2 Temperatura
Los datos de temperatura ambiente fueron obtenidos del atlas interactivo del IDEAM
para la región de Cali, Valle del Cauca, más específicamente de la estación de
medición ubicada en la universidad del Valle (Figura 22).
73
Figura 22.
Temperatura ambiente
24,60
24,40
24,20
24,00
23,80
23,60
1 2 3 4 5 6 7 8 9 10 11 12
Mes
Un cultivo de tilapias tiene condiciones variables, sin embargo, para el diseño del
equipo se tiene en cuenta un proceso isotérmico.
74
4.5.1 Suposiciones de diseño
Para el desarrollo del diseño del equipo se realizan suposiciones que permitan
determinar cómo progresar en los cálculos y el dimensionamiento del equipo. Estas
suposiciones se observan a continuación.
75
Tabla 15
Condiciones iniciales
La relación gas-líquido (G/L) se obtiene de la división del flujo de aire sobre el flujo
de agua.
76
Figura 23.
La parte experimental llevada a cabo por el Freshwater Institute indicó que con una
relación G/L de 1,2 la fase gaseosa aumentó su concentración de dióxido de carbono
77
en 690% y se eliminó 20% menos dióxido de carbono que con una relación G/L de 10,
por lo cual, los valores de G/L altos son favorables para el proceso. [41]
Adicional al estudio realizado por el Freshwater Institute, los autores del libro Water
treatment principles and design, establecen un rango óptimo para el factor de stripping
en las torres empacadas; no olvidemos que esta variable surge de la multiplicación de
la relación G/L con la constante de Henry.
Figura 24
20
Relación G/L
15
10
0
0 2 4 6 8 10 12 14 16 18 20 22 24 26
Factor de Stripping
Como se puede observar en la Figura 23 entre mayor sea la relación G/L es más
beneficioso para el proceso debido a la cantidad de empaque necesario; basándonos
78
en ese aspecto se evalúa la relación G/L de 10; con esa relación G/L el factor de
stripping es de 11,8, es decir, se encuentra en el rango óptimo para el proceso.
Considerando los datos del estudio y del libro se opta por tomar una relación G/L de
10 para el proceso de aireación y desgasificación, con la ventaja de que propiciar altos
volúmenes de aire es económico.
Una eficiencia de remoción alta indica una mayor altura del empaque y por lo tanto
un aumento en la inversión de la empresa, sin embargo, se estableció con
representantes de Lucarvel que el equipo debería disminuir adicionalmente 2 mg/L
previendo cualquier cambio en las condiciones ambientales.
79
De acuerdo con las fuentes productoras de dióxido de carbono se determina que en
horas de la noche la concentración de este parámetro aumenta debido a la fase oscura
de la fotosíntesis y a la continua respiración de la especie cultivada, sin embargo, para
efectos del diseño del equipo la velocidad de generación de dióxido de carbono
permanece constante a lo largo del tiempo.
Ecuación 5.
𝑉𝐶𝑂2 = 𝑄 ∗ (𝐶𝑜 − 𝐶𝑒 )
𝐿 𝑚𝑔 𝑚𝑔
𝑉𝐶𝑂2 = 320 ∗ (30 − 17 )
𝑚𝑖𝑛 𝐿 𝐿
𝑚𝑔
𝑉𝐶𝑂2 = 4160
𝑚𝑖𝑛
El dióxido de carbono es el compuesto que se desea extraer del agua por medio de
un flujo de aire; la constante adimensional de Henry provee información del
componente en la fase gas y en la fase líquida.
De acuerdo con la base de datos del NIST (National Institute of Standards and
Technology) la constante de Henry es 0,034 mol/L atm. [42]
80
Ecuación 6.
𝐻𝑃𝐶
𝐻𝑌𝐶 =
𝑅𝑇
1
𝑚𝑜𝑙
0,034
𝐻𝑌𝐶 = 𝐿 𝑎𝑡𝑚
𝑎𝑡𝑚 𝐿
0,082 ∗ 298,15 𝐾
𝑚𝑜𝑙 𝐾
𝐻𝑌𝐶 = 1,1865
Ecuación 7.
Factor de stripping
𝐺
𝑆= ∗ 𝐻𝑌𝐶
𝐿
𝑆 = 10 ∗ 1,1865
𝑆 = 11,86
81
4.5.8 Características del empaque
Figura 25.
82
Figura 26.
Gráfica de Eckert
La relación de Eckert utiliza la relación gas-líquido másica; para el diseño del equipo
se hizo la suposición de la relación gas-líquido volumétrica, por lo cual es necesario
realizar la conversión.
Ecuación 8.
𝐺𝑚 𝑄𝑎 𝜌𝑔
= ( )∗( )
𝐿𝑚 𝑄 𝜌𝑙
𝑘𝑔
𝐺𝑚 1,18 3
= 10 ∗ ( 𝑚 )
𝐿𝑚 𝑘𝑔
997 3
𝑚
𝐺𝑚
= 0,0119
𝐿𝑚
83
Posteriormente se calcula el eje X para la ubicación en la gráfica de Eckert.
Ecuación 9.
𝐿𝑚 𝜌𝑔
𝑥= ∗√
𝐺𝑚 𝜌𝑙 − 𝜌𝑔
𝑘𝑔
1 1,18
𝑥= ∗√ 𝑚3
0,0119 𝑘𝑔 𝑘𝑔
997 − 1,18 3
𝑚3 𝑚
𝑥 = 2,90
El eje y es una expresión matemática que involucra la carga hidráulica del gas, el
coeficiente de empaque, las densidades de los líquidos y su respectiva viscosidad. Con
el valor obtenido se realiza el despeje de la carga hidráulica del gas.
Ecuación 10.
𝐶𝑓
2
𝐺𝑚 ∗ 3,28 ∗ 𝜇𝑙0,1
𝑦=
𝜌𝑔 ∗ (𝜌𝑙 − 𝜌𝑔 )
Ecuación 11.
Carga hidráulica del gas
2
𝑦 ∗ 𝜌𝑔 ∗ (𝜌𝑙 − 𝜌𝑔 )
𝐺𝑚 =√
𝐶𝑓 0,1
3,28 ∗ 𝜇𝑙
84
𝑘𝑔 𝑘𝑔 𝑘𝑔
0,00177 ∗ 1,18 3 ∗ (997 3 − 1,18 3 )
𝑚 𝑚 𝑚
𝐺𝑚 = √
176 𝑚−1 0,1
3,28 ∗ 0,000891 𝑃𝑎 ∗ 𝑠
𝑘𝑔
𝐺𝑚 = 0,280
𝑚2 ∗ 𝑠
Ecuación 12.
Carga hidráulica del agua
𝐺𝑚
𝐿𝑚 =
𝑄𝑎 𝜌𝑔
𝑄 ∗ 𝜌𝑙
𝑘𝑔
0,280 2∗𝑠
𝐿𝑚 = 𝑚
𝑘𝑔
1,18 3
10 ∗ 𝑚
𝑘𝑔
997 3
𝑚
𝑘𝑔
𝐿𝑚 = 23,586 2
𝑚 ∗𝑠
Con los datos obtenidos es posible calcular el área de la torre, para ello es necesario
el caudal de agua utilizado.
Ecuación 13.
Área de la torre
𝑄 ∗ 𝜌𝑙
𝐴=
𝐿𝑚
𝑚3 𝑘𝑔
0,0053 𝑠 ∗ 997 3
𝐴= 𝑚
𝑘𝑔
23,586 2
𝑚 ∗𝑠
85
𝐴 = 0,225 𝑚2
Con el área calculada se determina el diámetro tomando como base la ecuación del
diámetro de un círculo, pues el área calculada es el área transversal.
Ecuación 14.
Diámetro de la torre
𝐴
𝐷 = 2∗√
𝜋
0,225 𝑚2
𝐷 = 2∗√
𝜋
𝐷 = 0,536 𝑚
Figura 27.
Correlación de Onda
86
La extracción de dióxido de carbono del agua depende de la difusividad de este
compuesto en el agua y en el aire. Para el análisis del sistema en cuanto a la
transferencia de masa es necesario conocer el valor de la constante de global de
transferencia de masa del dióxido de carbono en la fase líquida (K La), el cual puede ser
calculado mediante la correlación de Onda.
Tabla 16.
87
Los datos σc y σ corresponden a la tensión superficial del material del empaque y del
agua respectivamente. Para el polipropileno este valor corresponde a 0,033 kg/s 2 y
para el agua a 25°C es de 0,072 kg/s2.
Ecuación 15.
Correlación de Hayduk-Laudie
13,26 ∗ 10−9
𝐷𝑙 =
𝜇𝑙1,14 ∗ 𝑉𝑏0,589
Tabla 17.
Volumen
Elemento atómico Circunstancia
(cm3/mol)
Aire 29,9
Arsénico 30,5
Bismuto 48,0
Bromino 27,0
Carbon 14,8
Flúor 8,7
3,7 En compuestos orgánicos
Hidrógeno
7,2 Hidrógeno molecular
Mercurio 19,0
15,6
Nitrógeno 10,5 Aminas primarias
12,0 Aminas secundarias
88
Tabla 17(Continuación)
Tabla 18
Volumen molar
Compuesto # Moléculas
(cm3/mol)
Carbono 1 14,8
Oxígeno 2 14,8
Dióxido de carbono 29,6
Ecuación 16.
13,26 ∗ 10−9
𝐷𝑙 =
𝑐𝑚3
0,8911,14 𝑐𝑃 ∗ 29,60,589
𝑚𝑜𝑙
𝐷𝑙 = 2,06 ∗ 10−9
89
Para el cálculo del coeficiente de difusión en la fase gaseosa se utiliza la correlación
de Wilke-Lee.
Ecuación 17.
Correlación de Wilke-Lee
1 1 1 1
(1,084 − 0,249 ∗ √𝑀 + 𝑀 ) ∗ (𝑇 1,5 ) ∗ (√𝑀 + 𝑀 )
𝐴 𝐵 𝐴 𝐵
𝐷𝑔 =
𝑘𝑇
𝑃𝑙 ∗ (𝑟𝐴𝐵 )2 ∗ 𝑓 (𝜀 )
𝐴𝐵
Tabla 19
CO2 A
Aire B
Ecuación 18
Cálculo de rA
1
𝑟𝐴 = 1,18 ∗ 𝑉𝐴3
90
1
𝑟𝐴 = 1,18 ∗ 0,02963
𝑟𝐴 = 0,37 𝑛𝑚
Ecuación 19
Cálculo de εA/k
𝜀𝐴
= 1,21 ∗ 𝑇𝐵,𝐴
𝑘
𝜀𝐴
= 1,21 ∗ 194,54
𝑘
𝜀𝐴
= 235,39
𝑘
El cálculo de la función de colisión se realiza con las ecuaciones 20, 21, 22 y 23.
Ecuación 20.
Cálculo de εAB/k
𝜀𝐴𝐵 𝜀𝐴 𝜀𝐵
=√ ∗
𝑘 𝑘 𝑘
𝜀𝐴𝐵
= √235,39 ∗ 78,6
𝑘
𝜀𝐴𝐵
= 136,02
𝑘
Ecuación 21.
Cálculo de ee
𝑘𝑇
𝑒𝑒 = 𝑙𝑜𝑔10 ( )
𝜀𝐴𝐵
91
298
𝑒𝑒 = 𝑙𝑜𝑔10 ( )
136,02
𝑒𝑒 = 0,34
Ecuación 22.
Cálculo de ξ
𝜉 = −0,2825
Ecuación 23.
Cálculo de f(kT/εAB)
𝑘𝑇
𝑓 ( ) = 10𝜉
𝜀𝐴𝐵
𝑘𝑇
𝑓 ( ) = 10−0,2825
𝜀𝐴𝐵
𝑘𝑇
𝑓( ) = 0,5218
𝜀𝐴𝐵
Ecuación 24.
Cálculo de rAB
𝑟𝐴 + 𝑟𝐵
𝑟𝐴𝐵 =
2
92
0,37 𝑛𝑚 + 0,3711 𝑛𝑚
𝑟𝐴𝐵 =
2
𝑟𝐴𝐵 = 0,368 𝑛
Una vez se tiene la totalidad de los datos se aplica la correlación de Wilke-Lee para
el coeficiente de difusión en la fase gaseosa.
Ecuación 25.
1 1 1 1
(1,084 − 0,249 ∗ √44,01 + 28,96 ) ∗ (2981,5 ) ∗ (√44,01 + 28,96)
𝐷𝑔 =
89669 𝑃𝑎 ∗ (0,368 𝑛𝑚)2 ∗ 0,5218
𝑐𝑚2
𝐷𝑔 = 0,1989
𝑠
−5
𝑚2
𝐷𝑔 = 1,99 ∗ 10
𝑠
Tabla 20.
93
Tabla 20 (Continuación)
Dg 1,99*10-5 m2/s
ρg 1,18 kg/m3
μg 1,85*10-5 Pa s
at 209 m-1
dp 0,025 m
g 9,81 m/s2
σc 0,033 kg/s2
σ 0,072 kg/s2
dp < 0,0508 m
0,8 < Lm < 43 kg/m2 s
0,014 < Gm < 1,7 kg/m2 s
Nota. Restricciones de la correlación de
Onda. Tomado de: MWH Water
Treatment. “Strippng and aereation
towers”, [En línea]. Disponible:
https://cutt.ly/XnWjr8F [Acceso 6 jun,
2021.
Los datos presentados cumplen con las restricciones. El primer paso es calcular los
números adimensionales de Reynolds, Froude y Weber.
Ecuación 26.
Número de Reynolds
𝐿𝑚
𝑅𝑒 =
𝑎 𝑡 ∗ 𝜇𝑙
𝑘𝑔
23,59
𝑅𝑒 = 𝑚2 𝑠
𝑘𝑔
209 𝑚−1 ∗ 0,000891
𝑚𝑠
𝑅𝑒 = 126,66
94
Ecuación 27.
Número de Froude
𝐿2𝑚 ∗ 𝑎𝑡
𝐹𝑟 =
𝜌𝑙 ∗ 𝑔
𝑘𝑔 2
(23,59 ) ∗ 209 𝑚−1
𝑚2 𝑠
𝐹𝑟 =
𝑘𝑔 2 𝑚
(997 3 ) ∗ 9,81 2
𝑚 𝑠
𝐹𝑟 = 0,0119
Ecuación 28.
Número de Weber
𝐿2𝑚
𝑊𝑒 =
𝜌𝑙 ∗ 𝑎𝑡 ∗ 𝜎
𝑘𝑔 2
(23,59 )
𝑚2 𝑠
𝑊𝑒 =
𝑘𝑔 𝑘𝑔
997 3 ∗ 209 𝑚−1 ∗ 0,072 2
𝑚 𝑠
𝑊𝑒 = 0,037
Ecuación 29.
𝜎𝑐 0,75
𝑎𝑤 = 𝑎𝑡 ∗ (1 − exp [−1,45 ∗ ( ) ∗ (𝑅𝑒)0,1 ∗ (𝐹𝑟)−0,05 ∗ (𝑊𝑒)0,2 ])
𝜎
95
𝑎𝑤 = 209 𝑚−1
∗ (1
𝑘𝑔 0,75
0,033 2
− exp [−1,45 ∗ ( 𝑠 ) ∗ (126,66)0,1 ∗ (0,0119)−0,05 ∗ (0,037)0,2 ])
𝑘𝑔
0,072 2
𝑠
𝑎𝑤 = 119,34 𝑚−1
El área mojada del empaque es el área de cada empaque que va a estar en contacto
con el líquido.
Ecuación 30.
2 1
𝐿𝑚 3 𝜇𝑙 −0,5 0,4 𝜌𝑙 −3
𝑘𝑙 = 0,0051 ∗ ( ) ∗( ) ∗ (𝑎𝑡 ∗ 𝑑𝑝 ) ∗ ( )
𝑎 𝑤 ∗ 𝜇𝑙 𝜌𝑙 ∗ 𝐷𝑙 𝜇𝑙 ∗ 𝑔
2
−0,5
𝑘𝑔 3 𝑘𝑔
23,59 2 0,000891 𝑚 𝑠
𝑘𝑙 = 0,0051 ∗ ( 𝑚 𝑠 ) ∗( )
−1 𝑘𝑔 𝑘𝑔 −9 𝑚2
119,34 𝑚 ∗ 0,000891 𝑚 𝑠 997 3 ∗ 2,06 ∗ 10
𝑚 𝑠
1
−
𝑘𝑔 3
997 3
∗ (209 𝑚−1 ∗ 0,025 𝑚)0,4 ∗ ( 𝑚 )
𝑘𝑔 𝑚
0,000891 𝑚 𝑠 ∗ 9,81 2
𝑠
𝑚
𝑘𝑙 = 3,6 ∗ 10−4
𝑠
96
Ecuación 31.
1
0,7
𝐺𝑚 𝜇𝑔 3 −2
𝑘𝑔 = 5,23 ∗ (𝑎𝑡 ∗ 𝐷𝑔 ) ∗ ( ) ∗( ) ∗ (𝑎𝑡 ∗ 𝑑𝑝 )
𝑎 𝑡 ∗ 𝜇𝑔 𝜌𝑔 ∗ 𝐷𝑔
0,7
𝑘𝑔
𝑚 2 0,28 2
𝑘𝑔 = 5,23 ∗ (209 𝑚−1 ∗ 1,99 ∗ 10−5 )∗( 𝑚 𝑠 )
𝑠 𝑘𝑔
209 𝑚−1 ∗ 1,85 ∗ 10−5 𝑚 𝑠
1
𝑘𝑔 3
1,85 ∗ 10−5 𝑚 𝑠
∗( ) ∗ (209 𝑚−1 ∗ 0,025 𝑚)−2
𝑘𝑔 −5 𝑚2
1,18 3 ∗ 1,99 ∗ 10
𝑚 𝑠
𝑚
𝑘𝑔 = 0,015
𝑠
Ecuación 32.
1 1 1
= +
𝐾𝐿𝑎 𝑘𝑙 ∗ 𝑎𝑤 𝑘𝑔 ∗ 𝑎𝑤 ∗ 𝐻𝑌𝐶
1 1 1
= +
𝐾𝐿𝑎 3,6 ∗ 10−4 𝑚 ∗ 119,34 𝑚−1 0,015 𝑚 ∗ 119,34 𝑚−1 ∗ 1,1865
𝑠 𝑠
1
= 23,88 𝑠
𝐾𝐿𝑎
1
𝐾𝐿𝑎 = = 0,042 𝑠 −1
23,88
97
En 1998 Djebbar y Narbaitz realizaron una modificación a la correlación de Onda, en
la cual incluían una recalibración y ajuste en las constantes de cálculo, un número de
Reynolds para el flujo gaseoso y la eficiencia del empaque dentro de la determinación
del área superficial del empaque. Las recalibraciones requerían de una serie de
iteraciones suponiendo la altura de la columna empacada; disminuyendo el error un
20% frente a la correlación de Onda. [45]
Debido a la cantidad de tiempo que tomaba determinar la constante KLa por medio de
la recalibración de la correlación de Onda actualmente se utiliza un factor de seguridad.
Ecuación 33.
Factor de seguridad
𝐾𝐿𝑎 𝑑𝑖𝑠𝑒ñ𝑜
𝐹𝑆 =
𝐾𝐿𝑎 𝑂𝑛𝑑𝑎
Ecuación 34.
𝑍 = 𝐻𝑇𝑈 ∗ 𝑁𝑇𝑈
98
Ecuación 35.
𝑄
𝐻𝑇𝑈 =
𝐴 ∗ 𝐾𝐿𝑎
𝑚3
0,005 𝑠
𝐻𝑇𝑈 =
0,225 𝑚2 ∗ 0,029 𝑠 −1
𝐻𝑇𝑈 = 0,807 𝑚
Figura 28
99
la concentración de dióxido de carbono en el efluente con la concentración de dióxido
de carbono en la alimentación, arrojando un valor de 0,57.
Ecuación 36.
𝐶
𝑆 1 + (𝐶𝑜 ) ∗ (𝑆 − 1)
𝑒
𝑁𝑇𝑈 = ∗ ln [ ]
𝑆−1 𝑆
30
11,86 1 + (17) ∗ (11,86 − 1)
𝑁𝑇𝑈 = ∗ ln [ ]
11,86 − 1 11,86
𝑁𝑇𝑈 = 0,58
Para temas del cálculo de la altura de la torre se utiliza el valor de NTU calculado
debido a que el valor obtenido gráficamente puede verse sujeto a errores de lectura
Ecuación 37.
𝑍 = 0,807 𝑚 ∗ 0,58
𝑍 = 0,47 𝑚
100
La Figura 29 muestra el comportamiento típico del agua de acuicultura con bajo
contenido de oxígeno disuelto y alto contenido de dióxido de carbono disuelto.
Figura 29.
A medida que el aire entra en contacto con el agua los gases disueltos se acercan al
equilibrio con las presiones parciales en la atmósfera; como la molécula del agua está
en contacto con aire atmosférico para lograr el equilibrio se transfiere oxígeno dentro
de la molécula y dióxido de carbono fuera de ella. [41]
101
Figura 30.
Concentración de oxígeno en el
efluente
102
Figura 31
12
11
10
9
8
7
6
5
0 5 10 15 20 25 30 35 40
Temperatura (°C)
Ecuación 38.
Ley de Henry
𝑃𝑖 = 𝑥𝑖 ∗ 𝐻𝑖
103
Como la presión es directamente proporcional a la concentración de saturación,
cualquier cambio sobre la presión afecta de igual forma a la concentración de
saturación; la constante de Henry permanece igual con los cambios de presión.
Al conocer los valores asociados a la saturación del oxígeno a una presión de 1 atm
y 25°C se puede determinar la concentración de saturación a una presión de 0,89 atm
por medio de una relación obtenida de la ley de Henry.
Ecuación 39.
La presión parcial que ejerce el gas sobre el líquido depende de la concentración del
gas en el aire y la presión del sistema.
𝑚𝑔
8,244 𝐿 ∗ 0,89 𝑎𝑡𝑚
𝑥𝑖 (0,89 𝑎𝑡𝑚) =
1 𝑎𝑡𝑚
104
De esta forma se obtiene que la concentración de oxígeno disuelto en el efluente
para el proceso diseñado es de aproximadamente 7,34 mg/L.
Figura 32.
Constituyentes internos
105
4.7.1 Eliminadores de arrastre
Figura 33.
Eliminador de arrastre
Nota. Eliminador de
arrastre. Tomado de MTE,
“Demister Pads”, [En
línea]. Disponible:
https://cutt.ly/0nIGWG1
[Acceso 10 jun, 2021]
106
Se ubican en los fondos de la torre y sobre cada redistribuidor del líquido. Una de
sus funciones esenciales es lograr el mezclado del líquido antes de retirarlo del
proceso.
107
Como se puede observar en la Figura 34 una mala distribución del líquido deja
empaque seco, el cual no es eficiente para la transferencia de masa. Generalmente por
cada 0,1 m2 de sección transversal la columna debe disponer de al menos 5 puntos de
distribución del líquido para un proceso eficiente. [34]
Nota. Dimensiones de
la torre.
108
Los fluidos que participan en el proceso son el agua y el aire, para ello es necesario
contar con una bomba, un ventilador y tuberías para el transporte del agua.
4.8.1 Tuberías
4.8.2 Bomba
Al igual que con las tuberías el proceso de recambio de agua utiliza una bomba
centrifuga de la marca AQUA PAK de 1HP; con el fin de realizar una baja inversión
para el proceso se realiza la evaluación de este tipo de bomba.
Figura 36.
109
va a realizar el recorrido de las tuberías desde la bomba hasta la parte superior del
equipo.
Sin embargo, con las curvas de operación (Figura 37) se establece la cabeza
máxima a la que puede llegar la bomba con un caudal de 320 L/min.
Figura 37.
Así, se determina que los equipos de bombeo que posee la empresa son aptos para
el proceso de aireación y desgasificación propuesto.
4.8.3 Ventilador
110
Figura 38.
Tabla 21.
1,21 HP
Potencia
0,9 kW
Tensión Monofásica 220 V
Corriente máx. con carga 5,8 A
111
Tabla 21 (Continuación)
Con el diseño del equipo establecido se procede a realizar los balances de masa del
procedimiento. Para una mejor comprensión se realiza el diagrama de bloques (Figura
39).
112
Figura 39.
Diagrama de bloques
Nota. Diagrama de
bloques del proceso de
aireación y
desgasificación.
Ecuación 40.
𝐿2 ∗ 𝑥𝐻2 𝑂, 𝐿2 = 𝐿1 ∗ 𝑥𝐻2 𝑂, 𝐿1
113
Ecuación 41.
Ecuación 42.
Ecuación 43.
Los datos conocidos son los caudales de agua y aire (L y V) y las concentraciones
de las corrientes de entrada y salida del agua y la entrada del aire. El balance de
materia se realiza en base másica.
114
Tabla 22
Fracción
Corriente Sustancia
másica
CO2 3,02E-05
L2 O2 2,55E-06
H2O 0,99997
CO2 1,64E-05
O2 7,37E-06
L1
H2O 0,99996
N2 1,24E-05
O2 0,21
V1
N2 0,79
CO2 1,28E-03
O2 0,20
V2
N2 0,75
H2O 0,05
115
5. SIMULACIÓN DEL PROCESO
Tabla 23.
116
Tabla 23 (Continuación)
117
5.2 Selección del modelo termodinámico
Para la selección del modelo termodinámico a utilizar se toma como base el artículo
“Do not gamble with physical properties fot simulations” de Eric Carlson, en el cual se
establecen los criterios para una buena elección del modelo termodinámico.
Figura 40.
Árbol de decisión 1
118
Figura 41.
Árbol de decisión 2
En este tercer árbol de decisión se tienen en cuenta otros criterios, sin embargo, en
los modelos termodinámicos a los que se llega por cualquiera de los caminos está
presente el NRTL y el UNIQUAC.
119
Figura 42.
Árbol de decisión 3
120
5.3.1 Definición de sustancias y paquete termodinámico
Figura 43.
El siguiente paso es definir los compuestos que se ven afectados por la ley de
Henry. Como se determinó en la sección 4.7 la concentración del aire en el efluente es
determinada por medio de esta relación; por lo cual el oxígeno y el nitrógeno se ven
afectados por ella.
Figura 44.
121
En la siguiente sección se añade el paquete termodinámico seleccionado en el
apartado 5.2.
Figura 45.
En primer lugar, se definen las corrientes de entrada al proceso; por parte del aire se
establece el flujo de 3200 L/min que ingresa al blower con una composición de 0,21 de
oxígeno y 0,79 de nitrógeno; por parte del agua ingresa un flujo de 320 L/min a una
bomba con una composición de 0,97 de agua, 3*10^-5 de dióxido de carbono y
2,6*10^-6 de oxígeno.
122
Figura 46.
Figura 47.
En total el sistema contiene 6 corrientes, una antes y después del blower, una antes
y después de la bomba y las dos corrientes de salida, el agua limpia y el aire que
arrastró el contaminante.
Con especificar dos de las cuatro corrientes de entrada es suficiente debido a que
las corrientes mantienen sus propiedades al pasar por la bomba y el blower.
Los componentes internos son los relacionados con los empaques y la geometría de
estos; para ello se tiene en cuenta el tipo de anillo seleccionado, el material, el tamaño
y las medidas de la columna empacada.
123
Figura 48.
Figura 49.
124
5.3.3 Resultados
Figura 50.
Resultados de la simulación
Las corrientes de interés son las relacionadas con el agua, es decir, las corrientes
FROMTK y CLEAN, las cuales indican las corrientes de entrada y salida del proceso
respectivamente.
Por parte del oxígeno disuelto en el agua se obtiene una concentración final de 7,33
mg/L a partir de una concentración inicial de 2,5 mg/L y por parte del dióxido de
carbono se obtiene una concentración final de 16,33 mg/L partiendo de una
concentración aproximada de 30 mg/L.
Es necesario realizar la comparación entre los resultados obtenidos por medio de las
correlaciones y los obtenidos por la simulación.
Ecuación 44.
Porcentaje de error
125
Los resultados del porcentaje de error pueden verse reflejados en la Tabla 24.
Tabla 24.
Porcentajes de error
Teórico Experimental
Sustancia %Error
(mg/L) (mg/L)
Dióxido de carbono 17,00 16,33 3,95%
Oxígeno 7,34 7,33 0,11%
Nota. Porcentajes de error obtenidos en la
simulación.
Con los resultados obtenidos se concluye que el diseño realizado es viable para
aplicarlo en la empresa Lucarvel.
126
Figura 51.
40%
30%
20%
10%
0%
1 2 3 4 5 6 7 8 9 10
Relación G/L volumétrica
Por parte de la oxigenación del agua de cultivo la relación G/L no genera cambios
significativos en la concentración de oxígeno disuelto en el efluente debido a que, como
se mencionó anteriormente, la concentración de saturación de oxígeno en el agua se
logra usando una relación G/L menor a 1. Los resultados de ese análisis paramétrico
se reflejan en la Figura 52.
Figura 52.
7,70
7,50
7,30
7,10
6,90
6,70
6,50
1 2 3 4 5 6 7 8 9 10
Relación G/L volumétrica
127
Del análisis anteriormente realizado se concluye que con los datos de altura y
diámetro establecidos en el capítulo 4 la relación G/L mínima para cumplir con los
objetivos de la empresa es de 10.
128
6. COSTOS DE IMPLEMENTACIÓN
En este capítulo se determinarán los costos generales del equipo para realizar el
proceso de aireación y desgasificación en el cultivo semi intensivo de tilapias en la
empresa Lucarvel.
Tabla 25.
Los proveedores venden como mínimo 1 m3 de empaque, para los Pall Rings de 25
mm de polipropileno el m3 tiene un costo de $2.800.000, sin embargo, cada equipo
129
necesita 0,1 m3 por lo que en la inversión inicial se tiene en cuenta un costo
proporcional a la cantidad necesaria.
Los equipos que fueron cotizados en Alibaba tienen precios aproximados debido a
las dimensiones específicas que tiene el diseño propuesto; por otro lado, la bomba
centrifuga se tiene en cuenta, pero no entra en la inversión inicial porque la empresa ya
cuenta con este equipo en sus instalaciones.
En los costos de operación se deben tener en cuenta los gastos asociados a la parte
operacional y al personal necesario para mantener el proceso
Tabla 26.
BÁSICO PRESTACIONES
Auxilio de Salud
Salario básico Pensión (12%) ARL I (0,512%)
transporte (8,5%)
$ 1.200.000 $ 106.454 $ 102.000 $ 144.000 $ 6.264
OTROS
Intereses de
Prima Cesantías Vacaciones Horas extra
cesantías (1%)
$ 100.000 $ 100.000 $ 1.000 $ 50.000 $ 121.136
TOTAL $ 1.930.854
El salario básico del operario consiste en un salario base de 1.200.000 más las
prestaciones reglamentarias por ley, así mismo el pago de horas extra los días que se
considere necesario.
130
La jornada laboral consta de 8 horas diarias de lunes a sábado y 4 horas los
domingos.
De acuerdo con los datos analizados se precisa mantener la torre encendida las 24
horas del día debido al alto nivel de producción de dióxido de carbono provocado por la
respiración y otros procesos metabólicos de los peces.
El mantener la torre encendida trae otra ventaja al proceso, el movimiento del agua
provocado por la corriente proveniente de la torre es beneficioso para la salud de los
peces.
El agua en los sistemas de producción acuícola debe ser mezclada por el viento y en
sistemas intensivos a través de medios mecánicos. [49] La circulación del agua es una
estrategia eficaz para mezclar las aguas superficiales con altos niveles de oxígeno con
aguas de fondo con bajo contenido de oxígeno [50], para así evitar la saturación de
oxígeno en un solo punto del tanque y distribuirlo uniformemente.
131
Figura 53.
Tabla 27.
132
6.2.3 Disminución de recambios
Tabla 28.
Valor total
Cantidad de Volumen Volumen de recambio Recambios
D (m) (Recambios
tanques (m3) (m3) mensuales
diarios)
5 5 19,63 9,82 5 $ 576.188,50
9 6 63,62 31,81 5 $ 2.239.742,10
12 4 113,1 56,55 5 $ 2.654.457,00
TOTAL MES $ 5.470.387,60
Nota. Costo de los recambios de agua cuando se aplica la torre de aireación y
desgasificación.
133
CONCLUSIONES
134
BIBLIOGRAFÍA
[4] L. Towers, «The Fish Site,» 1 Noviembre 2005. [En línea]. Available:
https://thefishsite.com/articles/tilapia-life-history-and-
biology#:~:text=Tilapia%20are%20shaped%20much%20like,dorsal%20fin%2
0is%20heavily%20spined.. [Último acceso: 12 Febrero 2020].
135
Bancolombia, Bogotá D.C, 2020.
[14] C. Webster, Tilapia: Biology, Culture and Nutrition, Food Product Press,
2006.
136
[19] H. G. Maldonado Quimi , «Diseño de planta de tratamiento de agua potable
para el sitio El Recuerdo por el método de aireación y filtración,» Universidad
técnica de Machala, Machala, 2015.
[27] A. Bernhard, «The nitrogen cycle: Processes, Players, and Human Impact,»
Nature Education, vol. 2, nº 12, p. 9, 2010.
137
TILAPIA,» Agencia Española de Cooperación Internacional, Lima, 2004.
[29] S. Burge, «The Systems Engineering Tool Box,» Burge Hughes Walsh,
London, 2009.
138
LÍQUIDO,» INSTITUTO TECNOLOGICO DE CELAYA, Celaya, 1998.
139
[46] J. Staudinger, «Evaluation of mass transfer correlations for packed column
air stripping of volatile organic contaminants from water supplies,» Virginia
Polytechnic Institute and State University, Blacksburg , 1986.
[49] C. Boyd, «Por qué la mezcla y el movimiento del agua es clave en los
sistemas acuicolas,» Global Aquaculture Alliance, 2018.
[50] F. Kubitza, «La correcta circulación del agua en los estanques, Parte 2,»
Global Aquaculture Alliance, 2016.
140
[55] J. A. Hargreaves, «Biofloc Production System for Aquaculture,» Southern
Regional Aquaculture Center, nº 4503, p. 12, 2013.
141
[64] D. A. Chavarro Rodriguez, «DISEÑO DE UNA PLANTA DE TRATAMIENTO
DE AGUA POTABLE DE 2 LTS/S PARA UNA POBLACIÓN DE 750
HABITANTES,» Universidad Militar Nueva Granada, Bogota, 2014.
142
[73] B. P. Lárraga Jurado, «DISEÑO DEL SISTEMA DE AGUA POTABLE PARA
AUGUSTO VALENCIA, CANTON VICES, PROVINCIA DE LOS RIOS,»
Pontificia universidad catolica del Ecuador, Quito, 2016.
[80] D. Karimi , Eding , A. Aarnink y K. Groot , «The effect of gas to liquid ratio on
carbon dioxide removal and heat loss across a forced ventilated trickling filter,»
Wageningen University, Wageningen, 2020.
[81] J. Thom y W. Byers, «Limitations and Practical use of a Mass transfer model
for predicting air stripper performance,» Environmental Progress, Corvallis.
143
Universidad de Carabobo, Carabobo, 2016.
[87] E. Carlson, «Do not gamble with physical properties for simulations,»
Chemical Engineering Progress, Berkeley, 1996.
144
ANEXOS
145
ANEXO 1
RECOMENDACIONES
146