Diseño y Construcción de Un Banco de Pruebas de Un Sistema de Refrigeración
Diseño y Construcción de Un Banco de Pruebas de Un Sistema de Refrigeración
Benavides Gallegos, Joel Sebastian; Chávez Guerrero, Bryan Fernando; Reyes Utreras, José Daniel y
Certificación
sistema de refrigeración” fue realizado por los señores Benavides Gallegos, Joel Sebastian, Chávez
Guerrero, Bryan Fernando, Reyes Utreras, José Daniel y Yépez Rosero, Gandhy David, el mismo que
cumple con los requisitos legales, teóricos, científicos, técnicos y metodológicos establecidos por la
Universidad de las Fuerzas Armadas ESPE, además fue revisado y analizada en su totalidad por la
herramienta de prevención y/o verificación de similitud de contenidos; razón por la cual me permito
Firma:
__________________________________
C. C. 110357467 - 7
4
RESPONSABILIDAD DE AUTORÍA
Nosotros, Benavides Gallegos Joel Sebastian, Chávez Guerrero Bryan Fernando, Reyes Utreras José
Daniel y Yépez Rosero Gandhy David, con cédulas de identidad n° 172624987 - 1, 172210687 - 7,
172606302 - 5, 171854559 - 1, declaramos que el contenido, ideas y criterios del trabajo de titulación:
requisitos legales, teóricos, científicos, técnicos, y metodológicos establecidos por la Universidad de las
Fuerzas Armadas ESPE, respetando los derechos intelectuales de terceros y referenciando las citas
bibliográficas.
Firmas:
__________________________ __________________________
Benavides Gallegos, Joel Sebastian Chávez Guerrero, Bryan Fernando
C.C. 172624987 - 1 C.C. 172210687 - 7
__________________________ __________________________
Reyes Utreras, José Daniel Yépez Rosero, Gandhy David
C.C. 172606302 - 5 C.C. 171854559 - 1
5
AUTORIZACIÓN DE PUBLICACIÓN
Nosotros, Benavides Gallegos Joel Sebastian, Chávez Guerrero Bryan Fernando, Reyes Utreras José
Daniel y Yépez Rosero Gandhy David, con cédulas de identidad n° 172624987 - 1, 172210687 - 7,
Firmas:
__________________________ __________________________
Benavides Gallegos, Joel Sebastian Chávez Guerrero, Bryan Fernando
C.C. 172624987 - 1 C.C. 172210687 - 7
__________________________ __________________________
Reyes Utreras, José Daniel Yépez Rosero, Gandhy David
C.C. 172606302 - 5 C.C. 171854559 - 1
6
Dedicatoria
Dedico este trabajo y todo mi esfuerzo durante mi vida académica a los principales precursores
de mis logros, mis padres Anita y Javier, quienes con su ejemplo de dedicación y sacrificio plantaron en
mi conciencia ideales de superación y valores que me han permitido llegar a cumplir mis anhelos, son
quienes me brindaron todo su apoyo y amor en todo momento con lo que saco fuerzas cada día para
seguir adelante.
A mi hermano Juan Daniel, por ser la luz que alegra nuestro hogar, ser mi cómplice en todo
momento y ser por quien me supero cada día para poder brindarle el mejor ejemplo de vida.
A mis amigos, mis fieles compañeros que con cada aventura junto a ellos he logrado percibir su
A mis padres, Luis y Paulina, quienes han sido el vínculo directo para lograr el desarrollo de este
trabajo, que con sus palabras de aliento y apoyo incondicional supieron guiar y abrir el camino de la
educación durante estos 25 años, decirles que la semilla que sembraron finalmente ha empezado a dar
frutos, a mi hermana que con todo el apoyo material y sentimental me guio en este camino, y sin duda a
mis abuelitos quienes han estado detrás de mi camino por todo este tiempo.
El presente trabajo lo quiero dedicar a Dios por estar conmigo desde siempre y amarme por
como soy. A mi madre Soraya, quien es la persona que más amo y valoro en este mundo, pues ha estado
conmigo brindándome su apoyo, sus consejos y demostrándome cuanto me ama con cada uno de sus
actos. A mi hermano Christopher, por ser mi mayor motor para seguir adelante, mi mano derecha y mi
persona de confianza y por quien yo daría la vida. A mi padre Franklin, por su apoyo incondicional.
A mi angelito del cielo, mi abuelita Marianita que, aunque su partida me rompió el corazón, ella
siempre estuvo orgullosa de mí y sabía que cada uno de mis logros alcanzados serán por ella y para ella.
A mis tíos Graciela, José y Enma; a mis primos Steev, Jheremy, Adriana y Luana, por siempre estar para
mí en las situaciones buenas y malas. A mis amigos Evelin, Ricardo, Andrea, Jhon, Joha y Erika, que han
Este trabajo está dedicado a mis padres, Hugo y Dorys, quienes me han sabido apoyar
hermano Kevin y a mi familia, quienes siempre han creído en mí y en la frase que alguna vez un niño
A mis abuelos, “Papá Yayo” y “Mamá Anita”, quienes con su saber de vida me han ayudado a
entender que, a pesar de las adversidades siempre existirá, por mínimo que sea, un momento para
- Paul Valéry
Agradecimientos
En primer lugar, quiero agradecer a mis padres y hermano por ser la luz que guía mi vida, son el
principal apoyo para cada uno de mis logros y la razón de seguir luchando para conseguirlos, espero la
vida me premie con todo el tiempo necesario para poder recompensarles todo lo que por mí han hecho.
Quiero también agradecer a mi familia porque en cada momento y recuerdo de mi vida siempre
han estado presente, siempre me han recibido en su hogar con ese cariño y aprecio que espero nunca
me falte.
Por último, quiero agradecer a las personas por las que este proyecto fue posible, a mi tutor,
Ing. Luis Carrión PhD, por brindarnos su apoyo y conocimiento a lo largo de todo el proceso, a mis
compañeros Bryan, Daniel y Gandhy por todo el esfuerzo que pusieron para poder culminarlo y a las
empresas Acero de los Andes y Termo-Sistemas quienes nos apoyaron con su experticia y talleres para
En primer lugar, quiero agradecer a Dios por darme la capacidad y voluntad para obtener el
anhelo más deseado, a mis padres agradecerles por sus conocimientos, consejos y la provisión necesaria
para culminar todo este trabajo. A mis familiares quienes brindaron su apoyo y conocimiento indudable
en cualquier momento.
Agradecer también al Ingeniero Luis Carrión que con el aporte de sus conocimientos se logró
cumplir el proyecto y como no agradecer a mis amigos, quienes estuvieron pendientes de mis logros y
saben el esfuerzo realizado en las aulas para finalmente decir: ¡Lo logramos!
Primero que nada, agradezco a Dios por permitirme llegar hasta esta etapa de mi vida, a pesar
de existir situaciones complicadas durante mi trayecto académico, éstas me han ayudado a crecer como
persona.
Carrera de Ingeniería Mecánica que me han posibilitado desarrollar mi etapa universitaria para
A mi tutor, Ing. Luis Carrión PhD por apoyarme en todo el proceso de integración curricular, por
A la Dra. Lucy Mora por brindarme todo el apoyo psicológico en las etapas más difíciles de mi
carrera.
Quiero agradecer a toda mi familia que, a lo largo de toda mi vida, han demostrado un apoyo
incondicional, han inculcado consejos y experiencias, y han sido siempre mi motivo de superación para
poder llegar a este día. A Vicente (QEPD) y Cecilia, mi tío y tía que me mostraron con cariño que un libro
A mis grandes amigos que la universidad me ayudo a encontrar, aquellos que han demostrado
que ser un amigo, es brindarte una sonrisa a pesar de la distancia en el tiempo. Por último, pero no
menos importante, quiero agradecer al Ing. Luis Carrión PhD, quien con su conocimiento, experiencia y
Índice de Contenidos
Dedicatoria .................................................................................................................................................... 6
Agradecimientos ........................................................................................................................................... 8
Nomenclatura ............................................................................................................................................. 28
Subíndices ................................................................................................................................................... 30
Resumen ..................................................................................................................................................... 31
Abstract ....................................................................................................................................................... 32
Antecedentes .......................................................................................................................................... 33
Objetivos ................................................................................................................................................. 37
Alcance .................................................................................................................................................... 37
Justificación ......................................................................................................................................... 38
Importancia ......................................................................................................................................... 39
Entalpia ............................................................................................................................................... 42
Energía Interna.................................................................................................................................... 42
Energía de Flujo................................................................................................................................... 43
Ebullición ............................................................................................................................................. 60
Condensación...................................................................................................................................... 63
Sistemas de Refrigeración....................................................................................................................... 75
Ciclos de Refrigeración........................................................................................................................ 75
Compresor........................................................................................................................................... 78
Condensador ....................................................................................................................................... 82
Refrigerantes........................................................................................................................................... 84
Mezclas ............................................................................................................................................... 86
Conceptos Básicos............................................................................................................................... 91
Capítulo III. Diseño y Construcción del Banco de Pruebas del Sistema De Refrigeración ........................ 127
Capítulo IV. Diseño e Implementación del Sistema de Control y Adquisición de Datos .......................... 219
Referencias................................................................................................................................................ 285
Índice de Tablas
Tabla 9 Parámetros para el cálculo de constantes para el control P, PI, PID. ......................................... 119
Tabla 11 Parámetros para selección de tubería en función de la capacidad frigorífica .......................... 128
Tabla 12 Valores admitidos de perdida de carga equivalente a temperatura de vaporización .............. 128
Tabla 15 Velocidad del refrigerante dentro del sistema para diferentes diámetros de tubería ............. 130
Tabla 18 Evaluación de los refrigerantes según criterio medio ambiental ............................................. 135
Tabla 19 Cálculos de la energía requerida para enfriar una diferencia de temperatura. ....................... 137
Tabla 24 Parámetros del agua y refrigerante para diseño del evaporador ............................................. 150
Tabla 26 Propiedades termo físicas del refrigerante estado líquido ....................................................... 155
Tabla 29 Propiedades Termo físicas del agua (unidades inglesas) .......................................................... 158
Tabla 30 Parámetros para el aceite y refrigerante en el diseño del condensador .................................. 161
Tabla 32 Propiedades termo físicas del refrigerante estado líquido ....................................................... 166
Tabla 36 Dimensiones finales de diseño para los serpentines sumergidos en coraza ............................ 171
Tabla 39 Datos para el cálculo de perdidas mayores en el sistema de alimentación ............................. 173
Tabla 40 Datos para el cálculo de perdidas menores en el sistema de alimentación ............................. 173
Tabla 52 Conexión ESP32 – WROOM – 32 a MCP6002 (U11:A, U12:A, U12;B, U13:A, U13:B) ............... 227
Tabla 53 Conexión ESP32 – WROOM – 32 a PC817 (U2, U14, U15, U16) ............................................... 228
Tabla 56 Conexión ESP32 – WROOM – 32 a MAX31855JASA (U6, U7, U8, U9) ...................................... 231
Tabla 66 Toma de datos para intercambiadores de calor de tubo concéntrico en flujo paralelo .......... 261
Tabla 67 Análisis de eficiencia para intercambiador de calor de tubo concéntrico ................................ 263
Tabla 68 Toma de datos de la temperatura de ingreso del agua a temperatura setting de 19,5 °C ...... 265
Tabla 69 Toma de datos de la temperatura de ingreso del agua a temperatura setting de 17,0 °C ...... 267
Tabla 70 Toma de datos de la temperatura de ingreso del agua a temperatura setting de 14,5 °C ...... 269
Tabla 71 Toma de datos de la temperatura de ingreso del agua a temperatura setting de 12,0 °C ...... 271
20
Tabla 73 Honorarios para profesionales de la Universidad de las Fuerzas Armadas “ESPE” .................. 273
Índice de Figuras
Figura 2 Interfaz HMI y esquematización gráfica del equipo Túnel de Aire Multipropósito ..................... 36
Figura 3 Trabajo de flujo por fluido en movimiento dentro de un volumen de control ........................... 43
Figura 4 Razón de incremento de energía en flujo másico dentro de un volumen de control ................. 47
Figura 15 Factor de fricción para flujo laminar totalmente desarrollado en tuberías. ............................. 67
Figura 20 Ciclo ideal por compresión de vapor, diagramas T-s y P-h ........................................................ 76
Figura 36 Plano compuesto por un eje real σ y un eje imaginario 𝑗𝜔 ..................................................... 100
Figura 37 Representación de los componentes del SoC SnapDragon 801 de Qualcom .......................... 102
Figura 60 Comparativa, modos de control proporcionales, integrales y proporcionales + integrales .... 125
Figura 62 Potencia Requerida vs Flujo Másico, para determinada variación de temperatura ............... 138
Figura 67 Esquema y distribución geométrica del evaporador serpentín coraza. .................................. 151
Figura 70 Esquema y distribución geométrica del condensador serpentín coraza ................................. 162
Figura 73 Curvas de rendimiento, bomba centrifuga SHIMGE – Modelo CPm ....................................... 177
Figura 74 Curvas de rendimiento, bomba periférica EBARA – Modelo PRA ........................................... 178
Figura 104 Línea de succión de agua desde evaporador hacia la bomba SHIMGE ................................. 215
Figura 105 Línea de descarga de agua desde bomba SHIMGE hacia intercambiador de calor ............... 216
Figura 106 Bomba periférica EBARA modelo PRA 0.50 M ....................................................................... 217
Figura 107 Línea de succión de aceite desde condensador hacia la bomba EBARA ............................... 217
Figura 108 Línea de descarga de aceite desde bomba EBARA hacia intercambiador de calor ............... 218
Figura 112 Placa con implementación de componentes del circuito de control .................................... 223
Figura 114 Conexión ESP32 – WROOM – 32 a MCP6002 (U12:A, U12;B, U13:A, U13:B)........................ 227
Figura 116 Conexión ESP32 – WROOM – 32 a PC817 (U2, U14, U15, U16) ............................................ 229
Figura 119 Conexión ESP32 – WROOM – 32 a MAX31855JASA (U6, U7, U8, U9) ................................... 232
26
Figura 129 Interfaz del programa Nextion Editor HMI ............................................................................ 244
Figura 132 Programación de botones de actuadores y salida de texto en HMI ...................................... 246
Figura 136 Páginas para parámetros de entrada flujo paralelo y contraflujo ......................................... 248
Figura 142 Conexión física de los sensores a la placa de control ............................................................ 254
Figura 143 Constantes proporcional, integral y derivativo para los controles PID ................................. 255
27
Figura 144 Ingreso de temperatura de refrigerante y aceite para el control PID ................................... 255
Figura 147 Temperatura de entrada y salida del Agua vs. Tiempo ......................................................... 262
Figura 148 Temperatura de entrada y salida del Aceite vs. tiempo ........................................................ 262
Figura 149 Eficiencia del intercambiador de calor de tubos concéntricos .............................................. 264
Figura 150 Comparativa de valores obtenidos para una temperatura seteada de 19,5 °C .................... 266
Figura 151 Comparativa de valores obtenidos para una temperatura seteada de 17 °C ....................... 268
Figura 152 Comparativa de valores obtenidos para una temperatura seteada de 14,5 °C .................... 270
Figura 153 Comparativa de valores obtenidos para una temperatura seteada de 12 °C ....................... 272
28
Glosario de Términos
Nomenclatura
𝑄̇ Flujo de calor 𝐽
̇
𝑄𝑚𝑎𝑥 Flujo de calor máximo (Potencia máxima) 𝑊
𝑄𝐶 Calor cedido al entorno exterior 𝑊
𝑅 Constante del gas 𝑘𝑃𝑎 𝑚3⁄𝑘𝑔 𝐾
𝑅𝑒 Resistencia Térmica °𝐶/𝑊
𝑅𝑒𝑚𝑎𝑥 Número de Reynolds máximo Adimensional
𝑇 Temperatura absoluta del gas 𝐾
𝑇𝑝 Temperatura promedio °𝐶
𝑈 Coeficiente global de transferencia de calor 𝑊/𝑚2 𝐾
𝑢 Energía interna 𝑘𝐽/𝑘𝑔
V Velocidad de la corriente libre, velocidad del fluido 𝑚/𝑠
𝑉̇ Flujo Volumétrico 𝑚3 /𝑠
𝑉𝑎 Volumen aspirado 𝑚3
𝑉𝑏 Volumen de barrido 𝑚3
𝑉𝑜 Voltaje 𝑉
𝑉𝑝 Volumen desplazado 𝑚3 ⁄𝑟𝑒𝑣
𝑣´ Volumen específico del gas 𝑚3 ⁄𝑘𝑔
𝜈 Viscosidad cinemática del fluido m2 /s
𝑥 Calidad del vapor Adimensional
Δ𝐸𝐶 Cambio de energía cinética del sistema 𝐽
𝛥𝐸𝑠𝑖𝑠𝑡𝑒𝑚𝑎 Cambio neto de la energía total del sistema 𝐽
Δℎ Diferencia de entalpias 𝑘𝐽/𝑘𝑔
Δ𝑃 Cambio de energía potencial del sistema 𝐽
∆𝑇 Diferencia de temperatura media logarítmica °𝐶
Δ𝑈 Cambio de energía interna del sistema 𝐽
30
Símbolos Griegos
α Constante Adimensional
𝛽 Coeficiente de expansión térmica, constante 1⁄°𝐶
𝛾 Constante Adimensional
𝜀 Constante Adimensional
𝜃 Temperatura
𝜆 Constante Adimensional
𝜇 Viscosidad dinámica del fluido kg/𝑚 s
𝜌 Densidad 𝑘𝑔⁄𝑚3
𝜏 Constante de tiempo
Subíndices
1 Entrada
2 Salida
𝑎 Agua
𝑎𝑐 Aceite
𝑐 Condensador
𝑐𝑜 Frío
𝑐𝑢 Cobre
𝑒 Evaporador
𝑓 Fluido
ℎ Caliente
𝑖 Referido al área interior
𝑙 Liquido
𝑜 Referido al área exterior
𝑟 Refrigerante
𝑅𝑒 Resistencia eléctrica
𝑠 Superficie
𝑠𝑎𝑡 Saturación
∞ Entorno exterior
31
Resumen
Energía y Mecánica cuenta con dos bancos de pruebas para intercambiadores de calor Agua – Aceite,
equipos que son utilizados para la instrucción académica, siendo importantes en el proceso de
calor Agua – Aceite, busca contribuir en el mejoramiento de dicha formación estudiantil permitiendo el
estudio didáctico de un nuevo equipo en su totalidad. Por lo tanto, el proyecto presenta los análisis
modelamientos dinámicos del sistema, útiles para la implementación de los sistemas de control,
logrando el correcto funcionamiento del equipo. Para llevar a cabo el banco de pruebas, se construyo: la
bancada mediante el diseño y simulación de estructuras metálicas que permita ser el responsable de
soportar el peso del equipo en su totalidad, tanques de almacenamiento aislados térmicamente y con
las propiedades necesarias para evitar la corrosión debido a los fluidos de trabajo, sistema de
refrigeración, sistemas de bombeo y sistemas de control, donde la disposición de los componentes que
conforman el banco de pruebas, buscan que el estudiante cuente con un entorno educativo.
Abstract
The University of the Armed Forces - ESPE, specifically the Department of Energy and Mechanics
Sciences, has two test benches for Water-Oil heat exchangers, equipment that is used for academic
instruction, being important in the learning process and training of students of Mechanical Engineering
and Mechatronics. The present project, design and construction of a test bench for a cooling system for
Water-Oil heat exchangers, seeks to contribute to the improvement of said student training, allowing
the didactic study of a new equipment in its entirety. Therefore, the project presents the mathematical
analyzes carried out for the design and selection of components that intervene in the refrigeration
system to be implemented, using theory taught throughout the Mechanical Engineering career such as
Thermodynamics, Heat Transfer and Fluid Mechanics. In addition, the dynamic modeling of the system is
exposed, useful for the implementation of control systems, achieving the correct operation of the
equipment. To carry out the test bench, the following was built: the bench through the design and
simulation of metallic structures that allow it to be responsible for supporting the weight of the
equipment in its entirety, thermally insulated storage tanks and with the necessary properties to avoid
corrosion due to working fluids, cooling system, pumping systems and control systems, where the layout
of the components that make up the test bench, seek to provide the student with an educational
environment.
Capítulo I. Generalidades
Antecedentes
El uso de los equipos de experimentación que se dispone en los laboratorios del Departamento
Existen equipos que han sido adquiridos con el fin de sustentar el conocimiento teórico junto
con la práctica, pero a pesar de que estos equipos reciben el mantenimiento para seguir funcionando
necesitan adaptarse a los nuevos cambios tecnológicos para que los estudiantes estén familiarizados
con el entorno industrial y profesional, cumpliendo con la exigencia para mejorar el proceso educativo.
procesos termo-físicos más utilizados en la industria para cumplir con diferentes propósitos dentro de
los diferentes procesos requeridos con el fin de desarrollar actividades de producción industrial, es por
ello por lo que el estudio de intercambiadores de calor y el análisis para mejorar su eficiencia se
desarrolla frecuentemente.
a la obtención del Título de Ingeniero Mecánico denominado “Ingeniería conceptual, básica y de detalle
de un banco de pruebas para intercambiadores de calor tipo tubo concéntrico y multi - tubos para flujo
paralelo y contraflujo utilizando como fluidos de trabajo aceites API y agua para el laboratorio de
conversión de la energía del DECEM.” Bajo la responsabilidad de los ingenieros Jorge P. Mejía y Juan F.
Pérez en el año 2011, el cual busca ser un medio experimental para la comprobación de valores
Calor impartida a los estudiantes de Ingeniería Mecánica (Mejía & Pérez, 2011).
34
Figura 1
intercambiadores de calor tipo tubo concéntrico y multi - tubos para flujo paralelo y contraflujo
utilizando como fluidos de trabajo aceites API y agua para el laboratorio de conversión de la energía del
DECEM (p. 113), por Mejía & Pérez, 2011, Repositorio digital ESPE.
35
cual se pueda realizar un estudio de los intercambiadores de calor de tipo tubo concéntrico y multi -
Universidad de las Fuerzas Armadas “ESPE”, además de ser capaz de trabajar con múltiples aceites API y
Otro trabajo realizado hace no más de un año por: Ambato Guasumba Pablo Xavier y Bravo Ríos
Juan Fernando, con el tema: “Recuperación energética, puesta a punto, automatización y adquisición de
datos del equipo “Túnel de Aire Multipropósito”, del Laboratorio de Conversión de Energía”. Realizaron
HMI, permitiendo que el equipo esté disponible para realizar las prácticas en laboratorio con una
interfaz programada en LabVIEW para controlar la unidad de refrigeración, finalmente se realizaron las
respectivas pruebas de funcionamiento para la puesta a punto del equipo (Ambato & Bravo, 2021).
A continuación, en Figura 2, se muestra la interfaz gráfica del proyecto realizado por los
ingenieros Ambato y Bravo en la cual se puede observar una esquematización gráfica del Túnel de Aire
Multipropósito.
36
Figura 2
Nota. Adaptado de Recuperación energética, puesta a punto, automatización y adquisición de datos del
equipo “Túnel de Aire Multipropósito”, del Laboratorio de Conversión de Energía (p. 106), por Ambato &
El Departamento de Ciencias de la Energía y Mecánica cuenta con dos bancos de pruebas para
intercambiadores de calor Agua – Aceite, los cuales son utilizados para la formación académica de los
independientes para lograr las condiciones de temperatura objetivo, desaprovechando el calor que
transfieren las unidades condensadoras al ambiente. Además, los equipos carecen de un entorno
amigable que permita mejorar la formación estudiantil, tanto en la disposición de sus componentes en
la estructura del equipo como en la adquisición de datos brindada. Ante este problema se presenta la
intercambiadores de calor Agua – Aceite, que permita el estudio didáctico del equipo.
37
Objetivos
Objetivo General
Objetivos Específicos
adecuado.
• Diseñar el sistema de control en base a las señales de entrada y de salida del sistema de
refrigeración.
Alcance
componentes principales. Una vez culminada la investigación bibliográfica del Sistema de Refrigeración,
Se seleccionarán, con el previo estudio bibliográfico, todos los componentes que conforman el
sistema de refrigeración. Se analizará la dinámica de sistemas del proceso fluido – térmico, los
resultados del mismo indicaran y descartaran los posibles sistemas de control más adecuados.
38
El banco de pruebas estará disponible para un flujo másico para los fluidos de trabajo entre 0.1 y
0.8 GPM en el agua, mientras que de 1 a 11 GPM para el aceite. En el caso de las temperaturas se
trabajo donde podrá manipular el banco de pruebas según sus necesidades. Finalmente, los datos
obtenidos a partir de las condiciones dispuestas serán obtenidos mediante termocuplas, dispuestas en
cada entrada y salida del intercambiador de calor, las cuales se procederán a mostrar a través de una
interfaz HMI.
Justificación e Importancia
Justificación
Los intercambiadores de calor tienen una amplía aplicación en la industria como, por ejemplo;
en la industria alimentaria se usa con el fin de esterilizar alimentos elevando la temperatura para
Como se mencionó existen diversas aplicaciones y usos de estos equipos en donde se necesita
elevar o disminuir la temperatura de un fluido de trabajo, para el buen funcionamiento de los procesos y
Los intercambiadores de calor en los procesos industriales tienen gran demanda productiva, por
lo que, dentro de sus procesos, en donde el fluido de trabajo deba ser calentado o refrigerado en
grandes cantidades, requieren de intercambiadores de calor eficientes para que exista la mayor
transferencia de calor.
39
relación con el área de Energía y Termo fluidos, analizando la eficiencia de los intercambiadoras de calor
formación estudiantil permitiendo el estudio didáctico de un nuevo equipo en su totalidad. Por tal
motivo, se busca en el proyecto presentar los análisis matemáticos realizados para el diseño y selección
Mecánica de Fluidos. Además, de exponer los modelamientos dinámicos del sistema, útiles para la
La finalidad de este proyecto es diseñar y construir un banco de pruebas que funcione para
realizar ensayos en intercambiadores de calor y poder analizar los resultados obtenidos, con el fin de
Importancia
importante para estudiar el comportamiento térmico de equipos de transferencia de calor, en donde los
parámetros se controlan bajo ciertas condiciones, lo que permite analizar las variables y elementos del
sistema.
Resulta claro que al disponer de un banco de pruebas que permita controlar temperaturas y
flujos de entrada y salida mediante un sistema de control y adquisición de datos, ayuda a realizar el
intercambiadores de calor y en base a los resultados desarrollar equipos de transferencia de calor más
eficientes.
Por otro lado, la implementación del banco de pruebas ayudara al aprendizaje de los
estudiantes de la carrera de ingeniería mecánica durante las prácticas que se realizan en el área de
Diseño térmico y Transferencia de calor, permitiendo familiarizarse con el proceso fluido – térmico y con
La termodinámica es la ciencia que estudia la energía y los cambios que esta puede generar en
la materia.
La primera ley de la termodinámica establece que la energía no se crea ni se destruye, sino que,
solo se transforma. Esto indica que, cualquier cantidad de energía debe justificarse dentro de la
interacción total de un sistema, cumpliendo con el principio de conservación de la energía. Por tal
Balance de Energía
Çengel & Boles (2012) señala que el principio de conservación de la energía esta expresado
como: “el cambio neto (aumento o disminución) de la energía total del sistema durante un proceso es
igual a la diferencia entre la energía total que entra y la energía total que sale del sistema durante el
cualquier tipo de sistema que experimenta un determinado proceso. Además, como se observa en
ecuación 2, el cambio neto también se expresa como el compendio de los cambios de energía interna,
cinética y potencial.
𝛥𝐸𝑠𝑖𝑠𝑡𝑒𝑚𝑎 = Δ𝑈 + Δ𝐸𝐶 + Δ𝑃 ( 2 )
42
𝑃𝑣´ = 𝑅𝑇 ; 𝑃 = 𝜌𝑅𝑇 ( 3 )
Entalpia
combinación de las sustancias con las que se trabajan. En este caso es significativo el conocer que la
entalpia es la combinación de dos propiedades: Energía Interna y la energía de flujo, representada como
la siguiente ecuación:
ℎ = 𝑢 + 𝑃𝑣´ ( 4 )
Energía Interna
La energía interna de un sistema hace relación a una parte de la energía total, desde el punto de
vista microscópica que abarca a la estructura y grado de actividad molecular, que no depende de los
marcos de referencia externos o macroscópicos como son las energías cinética y potencial. De este
modo, el conjunto de todas las formas de energía microscópicas se las denomina como energía interna y
interna está en función de la temperatura y no del volumen específico y la presión. (Çengel & Boles,
2012, p. 180)
𝑢 = 𝑢(𝑇) ( 5 )
43
Energía de Flujo
Cuando se analiza un volumen de control, existe un flujo de masa que pasa a través de sus
fronteras, para esto se requiere un trabajo que permita ingresar o sacarlo del volumen de control.
En Figura 3 se tiene un volumen de control, por el que se va a ingresar un fluido que pasa a
Figura 3
Nota. Adaptado de Termodinámica (p. 226), por Y. Cengel y M, Boles, 2012, Mc Graw Hill.
Considerando que el fluido tiene un Volumen, V, cuando este ingrese al volumen de control
mediante el empuje de un embolo imaginario que ejerce una presión P a través del área transversal del
embolo (Çengel & Boles, 2012), se obtiene una fuerza que aplica el émbolo al fluido, definido como:
𝐹 = 𝑃𝐴 ( 6 )
44
En efecto, para que el fluido sea empujado al volumen de control, la fuerza que actúa en el
émbolo debe desplazarse a lo largo de la distancia L. Quedando expresado el trabajo aplicado al fluido o
𝑊𝑓𝑙𝑢𝑗𝑜 = 𝑃𝑉 ( 7 )
Ahora bien, mediante la definición de entalpía y la ecuación que rige el estado de un gas ideal,
ecuación 4, se tiene:
ℎ = 𝑢 + 𝑅𝑇 ( 8 )
El cambio de energía interna y entalpia para gases ideales cuando pasa de un estado 1 a 2
𝑑𝑢 = 𝐶𝑣 (𝑇)𝑑𝑇 ( 9 )
2
∆𝑢 = 𝑢2 − 𝑢1 = ∫ 𝐶𝑣 (𝑇)𝑑𝑇 ( 10 )
1
𝑑ℎ = 𝐶𝑝 (𝑇)𝑑𝑇 ( 11 )
2
∆ℎ = ℎ2 − ℎ1 = ∫ 𝐶𝑝 (𝑇)𝑑𝑇 ( 12 )
1
integración, por lo que estos valores de energía interna y entalpia ya están tabulados para diversos
gases a pequeños intervalos de temperatura. De acuerdo con Shapiro la variación de los calores
45
específicos en intervalos pequeños de temperatura se pueden aproximar como lineal (Moran & Shapiro,
2011).
𝑢2 − 𝑢1 = 𝐶𝑣 ,𝑝𝑟𝑜𝑚 (𝑇2 − 𝑇1 ) ( 13 )
ℎ2 − ℎ1 = 𝐶𝑝 ,𝑝𝑟𝑜𝑚 (𝑇2 − 𝑇1 ) ( 14 )
Para gases ideales, la relación que existe entre los calores específicos a volumen y presión
𝐶𝑝
𝐾= ( 15 )
𝐶𝑣
Proceso isentrópico
decir una sustancia mantendrá el mismo valor de entropía al inicio y al final del proceso que se llevará a
∆𝑠 = 0 ó 𝑠2 = 𝑠1 ( 16 )
𝑘−1
𝑇2 𝑃2 𝑘
( 17 )
( ) =( )
𝑇1 𝑠=𝑐𝑜𝑛𝑠𝑡 𝑃1
46
aproximado de suponer que las propiedades de un fluido en un ducto son unidimensionales, donde las
propiedades solo varían en la dirección del flujo y por lo tanto en la sección transversal del tubo serán
𝑚̇ = 𝜌𝑉𝐴𝑐 ( 18 )
𝑚̇
𝑉̇ = 𝑉𝐴𝑐 = ( 19 )
𝜌
Para un flujo estacionario se tiene que las variables de flujo másico tanto de entrada como de
salida son iguales, por lo tanto, la energía que gane o pierda el fluido será igual a la razón de incremento
de energía obtenido por dicho flujo másico dentro de un volumen de control como se muestra en
𝑚̇𝑒𝑛𝑡𝑟𝑎𝑑𝑎 = 𝑚̇𝑠𝑎𝑙𝑖𝑑𝑎 = 𝑚̇ ( 20 )
𝑄̇ = 𝑚̇Δℎ = 𝑚̇𝐶𝑝 Δ𝑇 ( 21 )
47
Figura 4
Nota. Adaptado de Transferencias de calor y masa. Fundamentos y Aplicaciones (p. 13), por Y. Cengel,
tanto cinética como potencial, serán nulas, debido a la ausencia de cambio en velocidad y elevación. Por
lo que, el cambio de energía del sistema se definiría como se muestra en ecuación 22.
𝛥𝐸𝑠𝑖𝑠𝑡𝑒𝑚𝑎 = Δ𝑈 ( 22 )
Sabiendo que el estado estacionario que presenta un sistema se define como la ausencia de
cambios con respecto al tiempo. La energía del sistema no cambiara debido a que el estado del sistema
𝛥𝐸𝑠𝑖𝑠𝑡𝑒𝑚𝑎 = 0 ( 23 )
𝑑𝐸𝑠𝑖𝑠𝑡𝑒𝑚𝑎
=0 ( 24 )
𝑑𝑡
48
define que el balance de energía para sistemas en estado estacionario esta dado por la ecuación 25.
𝐸̇𝑒𝑛𝑡𝑟𝑎𝑑𝑎 = 𝐸̇𝑠𝑎𝑙𝑖𝑑𝑎 ( 26 )
En ecuación 26, se puede observar el balance energía para un sistema estacionario. Misma que
Figura 5
Nota. Adaptado de Transferencias de calor y masa. Fundamentos y Aplicaciones (p. 11), por Y. Cengel,
𝜕𝑇
𝑞 = −𝑘𝐴 ( 27 )
𝜕𝑛
En cuanto a la transición de estado estacionario, ocurre siempre que la temperatura interna del
cuerpo en todos los puntos, incluidas las superficies, no dependan del tiempo. Si la temperatura cambia
con el tiempo, el cuerpo transfiere o almacena energía (Pitts & Sissom, 1998). La tasa de ganancias
acumulada será:
𝜕𝑇
𝑞𝑎𝑐𝑢𝑚𝑢𝑙𝑎𝑑𝑎 = 𝑚𝑐𝑝 ( 28 )
𝜕𝑥
Se produce cuando un cuerpo solido está en contacto con un fluido en movimiento de diferente
temperatura a la del cuerpo, el fluido transportará energía por convección por el fluido (Pitts & Sissom,
𝑞 = ℎ𝐴(𝑇𝑠 − 𝑇∞ ) ( 29 )
Cuando existe propagación de ondas electromagnéticas, se puede transferir calor, dichas ondas
electromagnéticas ubicadas en el vacío completo, así como en cualquier medio, viene dado por la ley de
Stefan-Boltzmann.
50
𝑞 = 𝜎𝐴𝑇 4 ( 30 )
10-8 W/m2K4.
Debido a la emisividad de los cuerpos, los cuales no son ideales, es decir que no tienen emisividad
𝑞 = 𝜀𝜎𝐴𝑇 4 ( 31 )
Números Adimensionales
puede correlacionar variables en una sola expresión matemática para caracterizar determinadas
La ventaja del análisis dimensional recae en que se puede estudiar el comportamiento del
proceso. En el presente documento se puntualizará los números adimensionales que tienen relación con
Antes de mencionar cada uno de los números adimensionales es necesario recordar ciertas
definiciones.
• Velocidad de corriente libre: velocidad promedio del fluido en dirección del flujo.
• Longitud característica: longitud característica de interés (por ejemplo, longitud de una placa
plana o al diámetro).
• Viscosidad: Medida de la resistencia interna a la deformación. (Çengel & Ghajar, 2011, p. 382)
• Densidad del fluido: Cantidad de masa por unidad de volumen que ocupa determinado fluido.
51
década de 1880, quien publico una serie de artículos que describen el flujo en tuberías, llegó a
determinar que el régimen de flujo de un fluido depende de la razón que existe entre las fuerzas de
En Figura 6 observar los regímenes de flujo sobre una placa plana. En donde inicialmente el
flujo es laminar, a medida que aumenta la distancia desde el inicio de borde de ataque de la placa se
producen perturbaciones, en donde empieza a ocurrir un flujo de transición, para finalmente volverse
Figura 6
Nota. Adaptado de Fundamentals of Heat and Mass Transfer (p. 389), por F. P. Incropera & D. P. DeWitt
Adicionalmente también se puede caracterizar a un flujo que circula por un tubo horizontal
Figura 7
Nota. Adaptado de Procesos de Transferencia de Calor (p. 48), por D. Q. Kern, 1999, Mc Graw Hill.
• Flujo turbulento: Cuando en un fluido predominan las fuerzas de inercia sobre las fuerzas
observar que el fluido se encuentra en una fase caótica y desordenada por el incremento de su
velocidad.
• Flujo laminar: Cuando en un fluido predominan las fuerzas viscosas sobre las de inercia se dice
que el flujo se encuentra en régimen laminar. En Figura 7 se observa que el movimiento del
líquido es ordenado y las líneas de corriente fluyen paralelamente a lo largo del eje del tubo.
• Fujo transitorio: En esta región el flujo varía entre laminar y turbulento, antes de ingresar a la
La caracterización del flujo de un líquido dentro de una tubería está en función de la velocidad,
densidad, viscosidad y diámetro del tubo por donde circula el fluido. El número de Reynolds se lo
𝐹𝑢𝑒𝑟𝑧𝑎𝑠 𝑑𝑒 𝑖𝑛𝑒𝑟𝑐𝑖𝑎
𝑅𝑒 = ( 32 )
𝐹𝑢𝑒𝑟𝑧𝑎𝑠 𝑣𝑖𝑠𝑐𝑜𝑠𝑎𝑠
𝑉 𝐿𝑐 𝜌 𝑉 𝐿𝑐 4𝑚̇ 𝐺𝑓 𝐷
𝑅𝑒 = = = = ( 33 )
𝑣 𝜇 𝜋𝜇𝐷 𝜇
53
De acuerdo con Cengel, los números de Reynolds para determinar en qué régimen se encuentra
Tabla 1
Sin embargo, en muchos de los casos en donde el flujo opera en régimen de flujo turbulentos
Número de Prandtl. Cuando un fluido fluye sobre una superficie, se desarrolla una capa límite
de velocidad, pero cuando este fluye a una temperatura diferente a la temperatura de la superficie, se
desarrolla una capa límite térmica como se muestra en Figura 8 y es aquí en donde las partículas del
Figura 8
Nota. Adaptado de Transferencia de Calor y Masa (p. 383), por Y. Cengel y A. Ghajar, 2011, Mc Graw Hill.
momento y energía por difusión en las capas límite de velocidad y térmica. (Incropera & DeWitt, 1999)
Es decir, el espesor relativo entre las dos capas de velocidad y térmica, mediante la siguiente
expresión adimensional.
𝐷𝑖𝑓𝑢𝑠𝑖𝑣𝑖𝑑𝑎𝑑 𝑑𝑒 𝑚𝑜𝑚𝑒𝑛𝑡𝑜
𝑃𝑟 = ( 34 )
𝐷𝑖𝑓𝑢𝑠𝑖𝑣𝑖𝑑𝑎𝑑 𝑡é𝑟𝑚𝑖𝑐𝑎
𝐶𝑝 𝜇
𝑃𝑟 = ( 35 )
𝑘
En Tabla 2 se menciona el rango de los números de Prandtl para metales, aceites y gases.
Tabla 2
Números de Prandtl
0.7-1 Gases
Es decir que para valores muy pequeños menores que uno, el calor se difunde con mayor
físicas del fluido con la longitud característica por donde circula. Este número adimensional se lo expresa
de la siguiente manera:
ℎ 𝐿𝑐
𝑁𝑢 = ( 36 )
𝑘
adimensional de Nusselt, según el rango para el número de Prandtl y Reynolds, es decir, que la relación
dependerá de saber en qué estado se encuentra el fluido, laminar, transición o turbulento y también de
en estudios experimentales por la dificultad que se tiene para encontrar de manera teórica este valor.
estas varían según el sentido en el que circula el fluido de mayor temperatura al de menor temperatura
de un punto a otro. En Figura 9 se muestra la distribución del intercambio térmico cuando se considera
Figura 9
Transferencia de Calor (p. 498), por F. Kreith y R. Manglik y M. Bohn, 2012. Cengage Learning
∆𝑇𝑎 − ∆𝑇𝑏
∆𝑇 = ( 37 )
∆𝑇
ln (∆𝑇𝑎 )
𝑏
En un intercambiador de calor interactúan dos fluidos que generalmente están divididos por una
pared que suele ser un tubo según el tipo de intercambiador. La transferencia de calor siempre se
produce de un foco caliente hacia un foco frío, por lo que el calor se ve obligado a pasar por los
diferentes procesos existentes como: convección, conducción y radiación. Las resistencias asociadas a
estos procesos son útiles a evaluar cuando se realiza un diseño de un intercambiador. Al realizar un
diseño térmico es importante analizar los coeficientes de transferencia de calor que intervienen en todo
Figura 10
encuentran a la misma temperatura. Adaptado de Transferencia de Calor y Masa (p. 155), por Y. Cengel
y A. Ghajar, 2011.
En Figura 10 se muestran las resistencias que intervienen en el flujo de calor cuando se tienen
dos fluidos separados por la pared de un tubo, expresándose al coeficiente global de transferencia
como:
1
𝑈𝐴𝑜 = ( 38 )
𝑅𝑐𝑜𝑛 ,𝑅 + 𝑅𝑐𝑖𝑙 + 𝑅𝑐𝑜𝑛𝑣 ,∞
Donde:
1
𝑈= 𝑟
ln ( 𝑜 )
1 𝑟𝑖 1
𝐴𝑜 ( + + )
2𝜋𝑟𝑖 𝐿ℎ𝑅 2𝜋𝐿𝑘𝑐𝑢 2𝜋𝑟𝑜 𝐿ℎ∞
1
𝑈= 𝑟
1 ln ( 𝑟𝑜 ) 1
𝑖
2𝜋𝑟𝑜 𝐿 ( + + )
2𝜋𝑟𝑖 𝐿ℎ𝑅 2𝜋𝐿𝑘𝑐𝑢 2𝜋𝑟𝑜 𝐿ℎ∞
1
𝑈= ( 39 )
𝑑𝑜 do ln (𝑟𝑜 /𝑟𝑖 ) 1
+ +
𝑑𝑖 ℎ𝑅 𝑘𝑐𝑢 ℎ∞
58
En el apartado anterior se analizó los fenómenos o mecanismos relacionados con el flujo interno
en tubos, ampliando los conceptos en relación con los fenómenos que suceden cuando se estudia un
fluido bifásico. En esta sección se busca conocer cómo se produce la transferencia de calor sobre un
banco de tubos o superficies externas con el fin de calcular el coeficiente de convección exterior.
que diferentes problemas de ingeniería se relacionan con flujos alrededor de cilindros. Recae aquí el
estudio en intercambiadores de calor de coraza y tubos con el fin de encontrar correlaciones que
En aplicaciones con intercambiadores de calor, se recomienda evitar que el fluido que circula
por la coraza sea laminar, ya que a diferencia con el flujo turbulento se ve obligado a necesitar mayor
área de transferencia de calor para intercambiar el mismo calor entre los fluidos.
La convección natural sobre superficies no está restringida por las corrientes de convección
cuando se encuentran al interior de un tubo, por lo que en el transcurso que el fluido se calienta puede
elevarse hasta zonas donde hay fluido frio permitiendo aumentar la transferencia de calor por
convección.
tiene la pared caliente del tubo a una temperatura considerablemente alta y a medida que circula el
fluido frio al exterior del tubo se observa que en cada corriente del fluido externo tiende a elevarse la
temperatura.
59
Figura 11
Nota. Adaptado de Procesos de Transferencia de Calor (p. 257), por D. Kern, 1999. Mc Graw Hill.
convección libre analizados sobre un solo tubo o alambres, por lo que no es suficiente para usarlos en
aplicaciones industriales.
determinar el coeficiente de convección cuando se tiene flujo laminar al exterior de banco de tubos sin
encontrar un error considerable en el cálculo. Los autores desarrollaron un nomograma que permite
0.25
𝑘𝑓3 𝜌𝑓2 𝐶𝑝𝑓 𝛽 ∆𝑡
ℎ𝑜 = 116 (( ) ( )) ( 40 )
𝜇𝑓′ 𝑑𝑜
Donde 𝜇𝑓′ esta en centipoises, y ∆𝑡, se calcula mediante la diferencia entre la temperatura de la
superficie y la del fluido frío. Para el cálculo y uso de la ecuación se debe trabajar en unidades inglesas.
60
En los sistemas de refrigeración el refrigerante está sometido a cambios de fase durante el ciclo
líquido a gas y de vapor a líquido, se los denomina como fluidos que trabajan en flujo bifásico o en
cambio de fase. Esto ocurre cuando trabaja el refrigerante en el evaporador produciéndose la ebullición
como la condensación.
Analizar estos dos fenómenos se vuelve más complicado que cuando se estudia convección en
una sola fase, ya que a parte de las variables generales del proceso de convección se aumentan las
Las propiedades como la densidad, conductividad térmica, viscosidad o geometría del sistema
son suficientes para describir el proceso de convección en una sola fase. Sin embargo, al analizar un
cambio de fase, características como la presión, el calor latente de vaporización juegan un papel
importante en el cálculo del coeficiente de convección, además de otras propiedades que pueden
de condensación y ebullición tienen un valor más alto comparados con los procesos de una sola fase.
Ebullición
cambia de estado líquido a gaseoso. Para entender cómo se produce el cambio de fase es importante
conocer los mecanismos físicos que se relacionan en este proceso, sin antes recalcar que se abordaran
superficie es mayor a la temperatura de saturación del fluido que cambiara de fase, por lo que se
empiezan a formar las primeras burbujas en la capa limite térmica. (Kreith et al., 2012)
estacionario y por medio de una superficie se transfiere calor al líquido, hasta llegar a la evaporación,
ebullición de estanque, depende de que tan rápido se forman las burbujas sobre la superficie por lo que
es complicado de predecir.
Ebullición en flujo. Este mecanismo se produce cuando el fluido tiende a moverse por ayuda de
una fuente externa como una bomba o con la ayuda de un compresor que empuja al fluido mientras se
estanque. En este caso puede existir dos tipos de ebullición en flujo, interno o externo, dependiendo si
el fluido que cambia de fase se encuentra al interior de una tubería o sobre una superficie caliente.
Ebullición de Flujo Interno. El análisis en flujo interno se torna compleja ya que el vapor no se
escapa a una superficie, sino que tanto el líquido como vapor se mueven juntos.
Régimen de Flujo. En Figura 12 se observa las distintas etapas que se producen en la ebullición
en flujo dentro de un tubo y la variación del coeficiente de convección a lo largo de los regímenes de
flujo. En primer lugar, el calor que se transfiere al liquido al ingresar por el tubo es mediante convección
62
forzada, inmediatamente al llegar a la temperatura de saturación se generan las primeras burbujas que
crecen y son arrastradas por la corriente principal denominándose la zona de flujo con burbujas. Estas
burbujas empiezan a unirse para formar una sola masa de vapor conocido como régimen de flujo en
grupo o tapón, para después continuar al régimen de flujo anular , aquí se forma una película de líquido
continúan por el régimen de transición en donde se empiezan a notar manchas totalmente secas en la
superficie interna del tubo hasta llegar a vapor saturado en donde el coeficiente convección disminuye
Figura 12
Nota. Adaptado de Fundamentals of Heat and Mass Transfer (p. 671), por F. P. Incropera & D. P. DeWitt
Condensación
Este proceso se produce cuando el vapor se encuentra en contacto con una superficie a una
Condensación por Gotas. Este mecanismo a diferencia del anterior tiende a formar gotas sobre
la superficie, esto se produce cuando se recubre a la superficie con una sustancia que impide que se
una película de líquido que por acción de la gravedad cae hacia abajo, que aumenta al condesarse mayor
cantidad de vapor.
Figura 13
Modos de condensación
Nota. a) Condensación de película, b) Condensación por gotas. Adaptado de Fundamentals of Heat and
Mass Transfer (p. 671), por F. P. Incropera & D. P. DeWitt & T. L. Bergman & A. S. Lavine, 2011, John
una temperatura promedio entre la temperatura de la superficie y de saturación, por lo que expulsa
mayor calor. Por lo que se define al calor de vaporización modificado como: (Çengel & Ghajar, 2011)
∗
ℎ𝑓𝑔 = ℎ𝑓𝑔 + 0.68 𝐶𝑝𝑙 (𝑇𝑠𝑎𝑡 − 𝑇𝑠 ) ( 41 )
propiedades relacionadas al flujo bifásico por lo que es importante aclarar que el calor de vaporización y
densidad de vapor, se deben evaluar a la temperatura de saturación de vapor. Por otro lado, las
𝑇𝑠𝑎𝑡 + 𝑇𝑠
𝑇𝑓 = ( 42 )
2
de sus sistemas se produce al interior de las superficies o tubos. Chato recomienda usar, para fluidos
bifásicos al interior de tubos horizontales la siguiente expresión: (Çengel & Ghajar, 2011)
1/4
𝑔 𝜌𝑙 (𝜌𝑙 − 𝜌𝑣 )𝑘𝑙3 3
ℎ𝑖𝑛𝑡𝑒𝑟𝑛𝑜 = 0.555 ( (ℎ𝑓𝑔 + 𝐶𝑝𝑙 (𝑇𝑠𝑎𝑡 − 𝑇𝑠 ))) ( 43 )
𝜇𝑙 (𝑇𝑠𝑎𝑡 − 𝑇𝑠 )𝐷 8
en los fluidos.
𝑉22 − 𝑉12
𝑄̇𝑒𝑛𝑡 𝑛𝑒𝑡𝑎 + 𝑊̇𝑓𝑙𝑒𝑐ℎ𝑎,𝑒𝑛𝑡 𝑛𝑒𝑡𝑎 = 𝑚̇ (ℎ2 − ℎ1 + + 𝑔(𝑧2 − 𝑧1 )) ( 44 )
2
𝑉22 − 𝑉12
𝑞𝑒𝑛𝑡 𝑛𝑒𝑡𝑎 + 𝑤𝑓𝑙𝑒𝑐ℎ𝑎,𝑒𝑛𝑡 𝑛𝑒𝑡𝑎 = ℎ2 − ℎ1 + + 𝑔(𝑧2 − 𝑧1 ) ( 45 )
2
Figura 14
Nota. Adaptado de Mecánica de Fluidos. Fundamentos y Aplicaciones (p. 206), por Y. Cengel, 2012, Mc
Graw Hill.
𝑃
Al aplicar la definición de entalpia ℎ = 𝑢 + 𝜌, para seguidamente reorganizar la ecuación de la
energía, se obtiene el primer miembro de la ecuación que representa la energía mecánica que entra al
sistema mientras que el segundo miembro representa la salida de dicha energía y las perdidas
𝑃1 𝑉12 𝑃2 𝑉22
𝑤𝑓𝑙𝑒𝑐ℎ𝑎,𝑒𝑛𝑡 𝑛𝑒𝑡𝑎 + + + 𝑔𝑧1 = + + 𝑔𝑧2 + (𝑢2 − 𝑢1 − 𝑞𝑒𝑛𝑡 𝑛𝑒𝑡𝑎 ) ( 47 )
𝜌 2 𝜌 2
La entrada de trabajo se denota como 𝑊𝑏𝑜𝑚𝑏𝑎 y la salida de esta como 𝑊𝑡𝑢𝑟𝑏𝑖𝑛𝑎 , las cuales
tendrán perdidas tratadas por separado una de otra, obteniendo la siguiente ecuación en términos de
carga.
𝑃1 𝑉12 𝑃2 𝑉22
+ + 𝑧1 + ℎ𝑏𝑜𝑚𝑏𝑎 = + + 𝑧2 + ℎ𝑡𝑢𝑟𝑏𝑖𝑛𝑎 + ℎ𝐿 ( 48 )
𝜌𝑔 2𝑔 𝜌𝑔 2𝑔
Perdidas Mayores. Las pérdidas de carga se expresan como perdidas mayores o menores dentro
de un sistema de bombeo. Las pérdidas de carga mayores se definen en base a las pérdidas de presión.
𝐿 𝜌𝑉 2
𝛥𝑃 = 𝑓 ( 49 )
𝐷 2
𝛥𝑃 𝐿 𝑉2
ℎ𝐿,𝑚𝑎𝑦𝑜𝑟 = =𝑓 ( 50 )
𝜌𝑔 𝐷 2𝑔
Donde f, representa el factor de fricción de Darcy y es aplicado con criterio dependiendo del
tipo de flujo, ya sea laminar o turbulento. En Figura 15, se puede observar el factor de fricción para flujo
Figura 15
Nota. Adaptado de Mecánica de Fluidos. Fundamentos y Aplicaciones (p. 332), por Y. Cengel, 2012, Mc
Graw Hill.
1 𝜀/𝐷 2.51
= −2.0 𝑙𝑜𝑔 ( + ) ( 51 )
√𝑓 3.7 𝑅𝑒√𝑓
rugosidad equivalente en la tubería, esta última siendo característica dependiente del tipo de material a
utilizar en la tubería, debido a esto en Tabla 3 se presentan algunas rugosidades equivalentes para
Tabla 3
Rugosidad 𝜺
Tuberías de Platico
(𝒎𝒎)
Polietileno
0.002
(PE)
Cloruro de Polivinilo
0.02
(PVC)
Nota. Adaptado de Hidráulica en tuberías a presión (p. 15), por Fagros, 2015.
Perdidas Menores. El coeficiente de perdida 𝐾𝐿 , depende de los accesorios dentro del sistema
manera.
𝑉2
ℎ𝐿,𝑚𝑒𝑛𝑜𝑟 = 𝐾𝐿 ( 52 )
2𝑔
Obteniendo así la perdida de carga total, definida por la suma entre las pérdidas de carga
mayores y menores.
𝐿 𝑉2
ℎ𝐿,𝑡𝑜𝑡𝑎𝑙 = (𝑓 + ∑ 𝐾𝐿 ) ( 54 )
𝐷 2𝑔
𝑃1 𝑉12 𝑃2 𝑉22 𝐿 𝑉2
+ + 𝑧1 + ℎ𝑏𝑜𝑚𝑏𝑎 = + + 𝑧2 + ℎ𝑡𝑢𝑟𝑏𝑖𝑛𝑎 + (𝑓 + ∑ 𝐾𝐿 ) ( 55 )
𝜌𝑔 2𝑔 𝜌𝑔 2𝑔 𝐷 2𝑔
69
Intercambiadores de Calor
Para lograr un intercambio de calor entre dos fluidos con temperaturas distintas, existen
dispositivos que permiten dicho intercambio sin mezclarse entre ellos mediante la transferencia de calor
Se pueden tener dos tipos de distribución para el movimiento de los fluidos en los
intercambiadores de calor, el flujo paralelo permite que se muevan en una misma dirección y que su
entrada sea por el mismo extremo, caso contrario del flujo en contraflujo, que se mueven en direcciones
contrarias y sus ingresos también son opuestos (Çengel & Ghajar, 2011).
En base a las necesidades que se presentan en diferentes áreas de procesos y operaciones que
están relacionados con: La ingeniería, medicina, calefacción, climatización, etc., se dispone de diferentes
tipos de intercambiadores que se los puede clasificar desde varios puntos de vista, según sea: su
Figura 16
Nota. Adaptado de Introducción a la Termodinámica con Algunas Aplicaciones de Ingeniería (p. 666), por
J. Rodríguez.
intercambiadores en donde uno de los fluidos de trabajo cambia de fase en el transcurso que circulan de
un punto a otro dentro del intercambiador de calor. En aplicaciones de refrigeración es necesario que
los dispositivos que permiten extraer o ceder calor, desde o hacia un medio, trabajen con un fluido en
donde el cambio de fase es necesario para que pueda trabajar adecuadamente el sistema.
Las aplicaciones más comunes y existentes en nuestro alrededor son los evaporadores y
condensadores que se utilizan en sistemas de refrigeración, sea que funcione como bomba de calor en
donde se busca generar calor hacia un medio o que el sistema trabaje de manera contraria permitiendo
energía térmica entre un fluido y un medio circundante, con el objetivo de absorber la energía de lo que
se busca refrigerar. El refrigerante es el que se encarga de recibir esa energía cuando pasa de un estado
líquido a vapor, por sufrir este cambio de estado e intercambio energético se denominan a estos
el medio circundante y el fluido que se encarga de ceder energía al medio exterior, pasa de vapor a
estado líquido, por producirse este cambio energético se conoce a estos sistemas como condensadores.
fenómeno de transferencia de calor están divididos por una pared sólida, que a menudo tienen
un fluido que se encuentra generalmente en depósitos cilíndricos verticales, también se los conoce
como serpentines mojados o intercambiadores de calor serpentín coraza. Los serpentines más comunes
Figura 17
Nota. Adaptado de Procesos de Transferencia de Calor (p. 819), por D. Kern, 1999. Mc Graw Hill.
una manera más económica y sencilla ya que consiste en un solo tubo en forma helicoidal.
Figura 18
Nota. Adaptado de Procesos de Transferencia de Calor (p. 822), por D. Kern, 1999. Mc Graw Hill.
73
Intercambiador de Calor de Carcasa y Tubos. Contiene una gran cantidad de tubos enlazados a
una carcasa, donde el intercambio de calor se produce cuando un fluido se mueve por dentro y el otro
Figura 19
Nota. Adaptado de Transferencia de Calor y Masa (p. 631), por Y. Cengel, 2012, Mc Graw Hill.
Son clasificados según el número de pasos que ejecutan por los tubos y la coraza, tomando en
cuenta que en todos los intercambiadores de calor de este tipo se forma una U en su coraza (Çengel &
Ghajar, 2011).
por ejemplo el calentamiento de glicerina y enfriamiento del agua. Son generalmente utilizados en la
parte química y petroquímica. En la parte farmacéutica, se encuentran de una manera más amplia,
debido a que de esta manera se pueda evitar una contaminación cruzada entre los dos fluidos (Çengel &
Ghajar, 2011).
74
Galindo L. (2013) define a la eficiencia en intercambiadores de calor como “la potencia que
transmite el intercambiador de dimensiones finitas dividido por la potencia máxima que transmitiría si el
energía (𝑈 = ∞).
obtener el criterio de potencia máxima a una longitud infinita las temperaturas 𝑇ℎ sería igual a 𝑇𝑐 , por lo
que la potencia máxima entregada por el intercambiador de calor en contraflujo vendría dada por la
ecuación 56.
Dónde la expresión descrita 𝑚̇ 𝐶𝑝 corresponde a la menor inercia térmica entre los fluidos de
trabajo o de estudio.
𝑄̇ ( 57 )
𝑒=
̇
𝑄𝑚𝑎𝑥
75
Dónde 𝑄̇ representa el calor cedido o perdido por el fluido caliente o frío respectivamente, de lo
cual podemos considerar al fluido con menor inercia térmica para dicha interpretación, debido a que
dicho fluido será aquel con mayor diferencia de temperaturas, obteniendo así la simplificación de la
ecuación 58:
(𝑇𝑒 − 𝑇𝑠 )𝑚𝑎𝑥 ( 58 )
𝑒=
(𝑇ℎ,𝑒 − 𝑇𝑐,𝑒 )
Sistemas de Refrigeración
cual Cengel & Boles (2012) lo define como la transferencia de calor producida desde una zona o región
de temperatura inferior a otra con una temperatura relativa superior. Los dispositivos o elementos que
producen dicho proceso de refrigeración se conocen como refrigeradores o bombas de calor, y los ciclos
de distintos materiales, trasladan energía térmica (calor) entre dos o más focos o ambientes conforme a
se puede almacenar distintos productos de interés como alimentos, compuestos químicos u otros.
Ciclos de Refrigeración
Existen dos tipos de ciclo de refrigeración, donde podemos notar al ciclo ideal y el ciclo real,
detallados a continuación.
Ciclo Ideal. Es el ciclo más utilizado en los sistemas de refrigeración, el cual contiene 4 etapas,
Figura 20
Nota. Adaptado de Termodinámica (p. 619, 620), por Y. Cengel, 2012, Mc Graw Hill.
ingresando al compresor que lo comprime hasta llegar isoentrópicamente a la presión del condensador,
para poder ingresar al condensador y ceder el calor hacia la periferia (Çengel & Boles, 2012).
presión de trabajo del evaporador, es decir, al estado 4 donde se transforma en vapor húmedo y con
77
baja calidad. El evaporador absorbe el calor del espacio a refrigerar, lo cual provoca que el refrigerante,
salga en forma de vapor saturado y vuelve al punto 1 para comenzar de nuevo el ciclo (Çengel & Boles,
2012).
Ciclo Real. Este ciclo difiere en las irreversibilidades en algunos de sus componentes, como la
fricción del fluido, que se debe a las caídas de presión y a la transferencia de calor en los alrededores ya
sea hacia el sistema o desde el sistema, en Figura 21 se presenta el diagrama T - s del ciclo real (Çengel &
Boles, 2012).
Figura 21
Nota. Adaptado de Termodinámica (p. 622), por Y. Cengel, 2012, Mc Graw Hill.
78
Para considerar en el campo real, se debe diseñar el sistema para que el refrigerante se
sobrecaliente ligeramente, antes de ingresar al compresor, lo cual garantiza que el refrigerante sea
totalmente vapor y no averíe el compresor. Adicionalmente, las caídas de presión producidas entre la
salida del evaporador y la entrada del compresor pueden llegar a ser muy representativas.
El compresor puede aumentar su potencial de trabajo, a causa del aumento del volumen
presión. En el ciclo real, no es fácil que después del condensador, el refrigerante sea liquido saturado
totalmente para poder ser estrangulado en la válvula que debe estar a una distancia muy corta del
evaporador e ingresar con baja entalpia, lo que produce que se extraiga más calor, tomando en cuenta
Un sistema de refrigeración por compresión mecánica está constituido básicamente por cuatro
• Compresor
• Condensador
• Válvula de expansión
• Evaporador
Compresor
comprensibles, pueden ser gases y vapores, sobre los cuales aporta un trabajo obteniendo como
Dentro del ciclo de refrigeración por compresión de vapor el compresor permite la compresión
del fluido de trabajo a condiciones de baja presión y temperatura, para descargarlo a propiedades de
alternante de pistón o émbolo para comprimir un gas y evacuarlo a alta presión (Figura 22).
Figura 22
Nota. Adaptado de Compresores Alternativos, por P. Fernández, 2012, Libros sobre Ingeniería Energética
([Link]
Compresor Rotativo. De igual forma son máquinas de desplazamiento rotativo que utilizan un
Figura 23
motor eléctrico y al compresor, los cuales van acoplados entre sí directamente y dentro de una cubierta
de acero, por lo cual no se puede realizar operaciones de mantenimiento. Un limitante de este tipo de
compresores es que deben ser alimentados exclusivamente con gas limpio, ya que no tienen la
capacidad de manejar gas con partículas líquidas o sólidas, por lo que es necesario realizar una limpieza
Figura 24
Compresor hermético
Nota. Adaptado de Compresor alternativo Modelo Serie AE3, por Tecumseh, 2010, Compresor
alternativo ([Link]
81
Generalmente en el exterior de la carcasa, llevan fijos tres tubos soldados, de los cuales dos son
del mismo diámetro y el otro de un diámetro menor, siendo el último el conectado a la descarga del
circuito mientras que de cualquiera de los otros se conecta a la aspiración, tal como se muestra en
Figura 24. Además, cuenta con una estructura de anclaje para sujeción a superficies planas y de un
Figura 25
([Link])
82
Condensador
transferencia de calor entre el refrigerante y un medio externo para así completar el proceso de
licuefacción del refrigerante una vez que sale del compresor y de eliminar el recalentamiento del gas
refrigerante.
• Tubo liso
Comúnmente utilizados en instalaciones pequeñas y son hechos con cobre. No poseen un flujo
forzado de aire ya que el aire que lo rodea al calentarse asciende y su volumen es sustituido por aire
Están conformados al igual que los condensadores de tubos lisos por serpentines de cobre, la
diferencia radica en la presencia de aletas separadas entre sí que aumentan el área de transferencia de
calor con el medio y siempre debe estar acompañado de una circulación de aire forzado. Las aletas
comúnmente son hechas de aluminio para no afectar propiedades del diseño en cuanto a peso.
Válvula de Expansión
Es un dispositivo cuyo objetivo es la restricción de flujo de un fluido para obtener una caída de
presión aprovechable para el proceso, la cual generalmente conlleva a una disminución de temperatura.
Dentro del campo industrial se utilizan 3 tipos distintos de elementos con el mismo objetivo, los cuales
Figura 26
Nota. Adaptado de Termodinámica (p. 617), por Y. Cengel, 2012, Mc Graw Hill.
propiedad conocida como coeficiente de Joule – Thomson, el cual según Cengel (2012) “es una medida
del cambio en la temperatura con la presión durante un proceso de entalpía constante”, resultando
igual si la temperatura permanece constante, mayor a cero si la temperatura disminuye y menor a cero
si la temperatura aumenta.
las más comunes son las de tipo termo-estáticas y las electrónicas, las cuales difieren en su
funcionamiento con respecto a la señal de entrada, la termo-estática recibe un estímulo físico en el gas
ese estímulo en forma de señal eléctrica mediante transductores para poder procesar dicha información
Campos de Aplicación
extremadamente bajas, por lo cual se han implementado sistemas de refrigeración en cascada, que
trabajen en serie para poder completar todo el ciclo de refrigeración. En el área de procesos científicos y
de ingeniería, siempre ha sido de vital importancia el estudio del ciclo de refrigeración, tal es el caso de
la separación del oxígeno y nitrógeno del aire, los cuales existen en forma gaseosa (Çengel & Boles,
2012).
Refrigerantes
Para que un sistema de refrigeración funcione de acuerdo al propósito por el cual fue diseñado y
construido, necesita de un refrigerante que permita transferir el calor desde un foco frio a un foco
Para lograr extraer calor de un espacio determinado y cederlo a otro, es indispensable el uso de
un refrigerante que estará sometido a cambios de estado, cuando pasa por una serie de procesos, los
que permiten al refrigerante volver al estado inicial dentro del ciclo de refrigeración. El refrigerante es
una sustancia capaz de absorber y ceder calor dentro de un sistema de refrigeración sin perder ninguna
de sus propiedades.
Figura 27
Refrigerantes Halogenados
hidrogeno, por lo que son químicamente muy estables al ser liberados a la atmosfera. Estos
años.
86
que les hace perjudícales a la capa de ozono, estos también contienen hidrogeno haciéndolos menos
HFC – Hidrofluorocarbono. Refrigerantes alternativos a los CFC y HCFC, estos contienen alto
potencial de calentamiento atmosférico, ya que al contener flúor se comportan como un gas de efecto
calentamiento atmosférico (PCA). Estos deben ser manipulados con precaución debido a su alta
propagación a la llama.
saturados. Su tiempo de vida es muy corta cuando se libera a la atmosfera y tiene valores muy bajos de
PCA.
Mezclas
Zeotrópica. A diferencia de la mezcla anterior, el cambio de fase del gas y el líquido varia en el
proceso de evaporación a una determinada presión, conociendo esa variación de temperatura como el
Refrigerantes no Orgánicos
Son refrigerantes alternativos que tienen valores muy bajos de PCA, no contienen carbono, su
Figura 28
Impacto Ambiental
Potencial de Agotamiento de Ozono – PAO. Es la capacidad que tiene una sustancia química
para causar daño a la capa de ozono. Se compara con el impacto de una masa de CFC-11. El valor de
referencia según la composición química se encuentra entre 0.01 y 1 para CFC y HCFC, respectivamente.
Calentamiento Global (PCG), que mide el efecto de gas invernadero. Se basa en la capacidad de atrapar
el calor de cada gas con respecto a la del CO2, así también como la rapidez que tienen de
Tiempo de Vida Media de la Sustancia. Se trata del tiempo que necesita la concentración para
que se disminuya en la atmosfera a la mitad. Si el tiempo de duración es larga, se lo considera como una
88
sustancia estable. Los gases fluorados tienen un periodo de vida muy extenso, por lo que el efecto que
• Clase 2L: Indica que son refrigerantes cuya inflamabilidad es inferior, con un máximo de
• Clase 3: Refrigerantes de alta propagación de la llama, cuyo LII (Límite Inferior de Inflamabilidad)
En Figura 29 se presenta como se clasifican los refrigerantes según los grupos de seguridad
Figura 29
Nota. Adaptado de Guía de Bolsillo de ASHRAE para Aire Acondicionado, Calefacción, Ventilación,
punto de burbujeo, rocío y nivel de seguridad. Los cuales servirán para seleccionar el refrigerante en un
Tabla 4
Nota. Esta tabla muestra los refrigerantes disponibles en el mercado adicionando el valor de media PAO
metodologías sistémicas y matemáticas para la representación de los procesos físicos. Varios conceptos
continuación.
Conceptos Básicos
Lazo de Control. Lazo de control es el conjunto formado por una unidad de medida, una unidad
de control, un elemento final de control y el propio proceso, generando así un bucle o lazo.
Un lazo de control puede ser abierto o cerrado. En el lazo de control abierto de la Figura 30, un
según su conveniencia para igualar el caudal de líquido de entrada con el de salida. Si existe una
variación regular de dicho parámetro será difícil para el operario de control mantener dicha condición,
por lo que deberá efectuar cambios con mayor frecuencia. Por otro lado, en el control de lazo cerrado
descrito en Figura 30, una vez hecha una primera calibración de parámetros de posición del vástago de
la válvula de control con la varilla de índice del flotador, propiamente dicho sistema regulará el nivel a
Figura 30
Nota. Adaptado de Instrumentación Industrial (p. 2), por A. Creus, 2010, Marcombo.
Campo de Medida. El campo de medida está definido por los valores de la variable medida o su
recepción o transmisión del instrumento y generalmente se lo expresa con dichos valores extremos.
(Creus, 2010)
Alcance. También conocido como span, el alcance consiste en la diferencia algébrica entre los
valor propio de la variable sujeta a medición. El error se lo clasificará dependiendo de las condiciones en
las que se desarrolle la experimentación, si sus condiciones son permanentes se conoce como error
93
estático, mientras que si las condiciones son variables y los instrumentos son susceptibles a las
variaciones del entorno se conoce como error dinámico, en el cual su valor dependerá de las
condiciones propias del proceso. Una relatividad en el error se la puede obtener mediante el error
medio del instrumento, que resulta de la media aritmética de los errores de la experimentación.
estadístico y resulta de una serie de mediciones, características de los equipos y condiciones varias.
Sensibilidad. Corresponde a la razón entre los valores del incremento de la lectura del
instrumento y el incremento de la variable que ocasiona dicha alteración, después de alcanzar un estado
de reposo.
al comparar una misma variable bajo las mismas condiciones, es decir exponer los mismos resultados
múltiples veces.
94
Histéresis. Corresponde a la máxima distancia que se evidencia en los valores de salida del
instrumento para un mismo valor medido en la variable, cuando se realiza una experimentación que
Carga. Es la cantidad del agente o fluido de control mediante el cual el proceso mantiene
Tipos de Entrada
En el sistema, se proporciona una relación causa y efecto, el cual será el proceso a controlar y se
pueda enlazar la salida con la entrada. Los tipos de entradas, en un sistema de control, se representa en
Figura 31.
Figura 31
Nota. Adaptado de introducción a los sistemas de control, aplicaciones y simulación con MATLAB (p. 2),
La entrada de tipo escalón, expresa una constante introducida en el sistema, para el caso de la
entrada tipo rampa, representa una variación continua en el tiempo; caso contrario en la entrada tipo
impulso, se lo realiza con una magnitud extremadamente grande, pero con una duración
Sensores y Transductores
Sensor. Dispositivo capaz de emitir una señal de salida transformable a partir de la energía
Figura 32
Nota. Adaptado de The measurement, instrumentation and sensors Handbook (p. 26), por J. Webster,
encuentran realizando la medición en: sensores pasivos y sensores activos. (Webster, 1999)
Sensores Pasivos o Moduladores. Son los que no agregan energía durante el proceso de
medición, pero pueden remover energía durante dicha operación. Como ejemplo de sensor pasivo se
describe la termocupla, la cual convierte una temperatura física a una señal de voltaje debido a que el
Sensores Activos o Generadores. Son aquellos que agregan energía durante el proceso de
medición. Un ejemplo son los radares o sistemas sonares, donde la distancia hasta objetos es medida
enviando señales de onda de radio o acústicas para ser reflejadas en un cuerpo obteniendo así la
Una segunda clasificación de los sensores se realiza considerando el tipo de señal que provee el
sensor, mediante este criterio se pueden clasificar en: sensores analógicos y sensores digitales.
Sensores analógicos brindan una señal que es continua tanto en su magnitud como en su
referencia espacial o temporal, esto quiere decir que la señal de salida que provee es continua y
directamente proporcional a la señal de entrada. Dentro de la naturaleza las variables físicas como
tienden a ser continuas y medibles con sensores analógicos, lo cual permite representarlas con una
Figura 33
Nota. Adaptado de The measurement, instrumentation and sensors Handbook (p. 37), por J. Webster,
Sensores digitales brindan una señal binaria (on/off ó 0/1), es decir que en dicha señal solo
existen valores discretos en el tiempo o espacio y dentro de este periodo de valores discretos la señal
puede ser representada solo con números discretos de valores de magnitudes. (Webster, 1999)
Un sistema binario es un sistema numérico en base dos, dentro del cual su representación más
identificar el número de vueltas de un sistema en rotación mediante el uso de un sensor magnético que
brinda una señal digital que representa un giro completo del sistema.
Figura 34
Nota. Adaptado de The measurement, instrumentation and sensors Handbook (p. 38), por J. Webster,
utilizadas generalmente.
98
Tabla 5
Magnitudes
Sensores
Caudal
Temperatura Presión Nivel
Flujo
Anemómetros de Potenciómetro +
RTD Potenciómetros + hilo caliente flotador
Resistivos
Termistores tubo Bourdon Galgas + voladizo Termistores
Termistores LDR
Condensador
Condensador
Capacitivos - variable + -
variable
diafragma
LVDT + diafragma
LVDT + flotador
Inductivos y Reluctancia LVDT + rotámetro
- Corrientes
electromagnéticos variable + Ley Faraday
Foucault
diafragma
Termopares
Generadores Piezoeléctricos - -
Piroeléctricos
Transductores. Dispositivo capaz de convertir una señal de particularidades físicas en una señal
correspondiente, pero de distintas particularidades físicas, por lo cual se los considera dispositivos
Dentro del campo habitual existen seis tipos de señales: eléctricas, magnéticas, mecánicas,
ópticas y químicas, por ende, cualquier dispositivo que transforme una señal de características
particulares en otra señal de otro tipo se lo debería considerar transductor, siendo su señal de salida
cualquier forma física aprovechable. En el ámbito práctico, los transductores son considerados a
aquellos que al convertir la energía su señal de salida es de naturaleza eléctrica. (Pallás Areny, 2005)
Dentro de un sistema, con todas sus partes e interconexiones pueden ser detallados
símbolo de la ejecución matemática a realizar en la entrada para obtener la salida (Figura 35 a). Cuando
existen varias entradas, el bloque, es cambiado por un círculo (Figura 35 b) y su salida resulta de la suma
algebraica de todas sus entradas, a esto se lo conoce como punto de suma. Adicionalmente, existe el
punto de reparto, que es generalmente utilizado cuando la señal es bifurcada para ser aplicada a varios
Figura 35
Nota. Adaptado de Introducción a los sistemas de control, aplicaciones y simulación con MATLAB (p. 4),
Transformada de Laplace
Transforma una función 𝑔(𝑡), del dominio del tiempo, en una función 𝐺(𝑠) en el dominio de s,
∞
𝐿{𝑔(𝑡)} = ∫ 𝑔(𝑡)𝑒 𝑠 𝑡 𝑑𝑡 = 𝐺(𝑠) ( 59 )
0
es un número complejo 𝑠 = 𝑠 + 𝑗𝜔, por lo tanto, la función 𝐺(𝑠) se representa en el plano cartesiano s
Figura 36
Nota. Adaptado de Introducción a los sistemas de control, aplicaciones y simulación con MATLAB (p. 26),
una ecuación algebraica de grado semejante. Las raíces resultantes del polinomio, se puede representar
de la siguiente manera:
(𝑠 + 𝑧0 )(𝑠 + 𝑧1 ) ⋯
𝐺(𝑠) = 𝐾 ( 60 )
(𝑠 + 𝑝0 )(𝑠 + 𝑝1 ) ⋯
Donde K representa la constante del sistema, las raíces del polinomio del numerador
sus raíces son llamados polos y son representados como una cruz para el plano s (Hernández, 2010).
El polo más próximo al eje 𝑗𝜔 se lo llama polo dominante, debido a que este influye más en el
Parte real. Un comportamiento, donde en el dominio del tiempo, se llega más rápidamente a un
valor final, para el dominio de s, la posición del polo dominante es un índice para la velocidad de
Parte imaginaria. Un comportamiento con respecto al dominio de s, la parte imaginaria del polo
SoC “System on a Chip”. Su traducción literal al español es “Sistema en un chip” y surge frente a
pequeños y autónomos posibles, con lo cual se crea un chip capaz de tener múltiples funciones
integradas.
102
que trabajan de manera colaborativa para tener funcionabilidad combinada, entonces un SoC se define
Un ejemplo ilustrativo de SoC se representa en Figura 37, el cual es el SoC SnapDragon 801 de
Qualcomm. Dentro de dicho elemento están presentes 8 componentes, los cuales son: Procesador (Krait
CPU), controlador de la memoria, la memoria, la GPU, el DSP, el módulo de conectividad (USB, LTE, WiFi
Figura 37
funcionamiento principal en los celulares y debido a su aplicación sus dimensiones son mínimas.
103
aplicaciones tecnológicas.
Las ventajas del dispositivo SoC es la versatilidad que posee debido a sus dimensiones, lo que lo
permite trabajar en distintos dispositivos, otra ventaja es la eficiencia que brinda ya que con su trabajo
colaborativo y dimensiones reducidas demanda un menor consumo energético. No obstante, tiene dos
desventajas puntuales, la primera radica en que en caso de algún fallo en cualquiera de sus
componentes todo el SoC queda inservible y la segunda es que debido al aumento en la densidad de
ESP32-WROOM-32U. “ESSP32” es un chip simple de 2.14 GHz con combo Wi-Fi y Bluetooth
diseñado con la tecnología TSMC ultra-low-power de 40 nm. Es diseñado para lograr el mejor
Es una placa de desarrollo o SoC con conectividad WiFi y Bluetooth basado en el módulo ESP32,
conectividad versátil.
arseniuro de galio en paralelo inverso acoplados a un fototransistor NPN de silicio en un paquete DIP de
5 pines. Está diseñado para aplicaciones que requieran detección o monitoreo de señales de AC.
104
Figura 38
Elemento H11AA1
Nota: Adaptado de Optocoupler, Phototransistor Output, Ac Input, with Base Connection, por Vischay,
2011.
Figura 39
Nota: Adaptado de Optocoupler, Phototransistor Output, Ac Input, with Base Connection, por Vischay,
2011.
IRM-20-5. Es un módulo miniatura de tipo fuente de energía de 20 W que está listo para
industrial. Permite un rango universal de voltaje de entrada de 85 ~ 305 VAC. Su estructura permite la
Tabla 6
Denominación Especificación
Figura 40
Elemento IRM-20-5
Figura 41
interruptor, el cual es activado mediante luz emitida por un LED que satura un componente
resistivas para operaciones de VCA o es utilizado para realizar un aislamiento entre dos etapas distintas
de un circuito.
Tabla 7
Denominación Especificación
Voltaje reverso 3V
altos y una alta capacidad de sobre corriente, se lo utiliza generalmente en el control de potencia en
Tabla 8
Denominación Especificación
Corriente en estado 40 𝐴
Figura 42
Elemento TPDV640
DAC. Es un dispositivo que permite convertir señales digitales a señales analógicas, de ahí que
sus siglas en ingles signifiquen “Digital to Analog Converter”. En Figura 43 se puede observar una señal
digital discontinua atravesando un filtro (DAC) para poder convertirse en una señal analógica continua
(Huircán, 2003).
Figura 43
Nota. Adaptado de Conversores Análogo-Digital y Digital-Análogo: Conceptos Básicos (p. 1), por J.
Huircán, 2003.
108
Figura 44
Nota. Ganancia de voltaje sin carga (A), Entrada no inversora ( V +) y Entrada inversora ( V −) Adaptado
termocuplas K, J, N, T o E. A su vez, tiene la capacidad de asimilar lecturas de hasta 1800 °𝐶 y −270 °𝐶.
En Figura 45 se representa un dispositivo MAX31855 con la descripción respectiva de cada uno de sus
Figura 45
Figura 46
Figura 47
Nota. Paquete libre de plomo - RoHS (+), empaquetamiento de componentes (T) y futuro producto (*)
Tiristores. Se constituye por 4 capas y dos terminales (ánodo y cátodo) que forman un
dispositivo, dichas capas son semiconductoras, generando la estructura pnpn como se muestra en la
interruptor ya que se encuentra apagado hasta que recibe voltaje para encenderse y poder conducir
(Floyd, 2008).
110
Figura 48
Tiristor (4 capas)
Nota. Adaptado de Dispositivos Electrónicos (p. 554), por T. Floyd, 2008, Pearson Educación.
Disipador térmico. Es una superficie aletada de aluminio cuyo propósito es la disipación del
calor generado por distintos elementos dentro de un ambiente cerrado, tal como la caja de control o la
aletas el dispositivo genera una mayor o menor transferencia de calor, siendo importante su correcta
Figura 49
Disipador de calor
Sistemas de Control
“Un sistema de control está formado por subsistemas y procesos (o plantas) unidos con el fin de
controlar las salidas de los procesos. En su forma más sencilla, un sistema de control produce una salida
o respuesta para una entrada o estímulo dado (Figura 50)” (Nise, 2006)
Figura 50
Nota. Adaptado de Sistemas de control para ingeniería (p. 2), por N. Nise, 2006, Compañía Editorial
Continental.
combustible. En dicho proceso los subsistemas se los identifica como válvulas de combustible y
controlado, además, existe el subsistema termostatos que funcionan como detectores midiendo la
temperatura de la habitación.
El uso de Sistemas de Control permite el movimiento de equipos con grandes dimensiones con
una precisión que no puede ser obtenida mediante otra forma, tal como lo hacen las enormes antenas
al dirigir sus captadores y emisores hacia regiones distintas del universo para realizar la comunicación de
Los requerimientos básicos satisfechos con la construcción de sistemas de control son los
• Amplificación de potencia
• Control remoto
constantes magnitudes, entre ellas la presión, el nivel, el caudal, la temperatura, la acidez, la velocidad,
la humedad, etc. El uso de instrumentos de control y medición permiten realizar esta función necesaria
de mantener las magnitudes, además de regular de manera idónea las constantes de un modo que
Un proceso cualquiera hace referencia a un sistema conformado o dispuesto de tal manera que
se cumpla el objetivo de manufacturar o cambiar las propiedades de un material mediante una serie de
calor, dentro del cual el controlador desarrolla las siguientes dos funciones para lograr el objetivo
Figura 51
Nota. Adaptado de Instrumentación Industrial (p. 482), por Creus, 2010, Alfaomega Marcombo.
El diagrama de bloques para el bucle de control descrito en Figura 51, se encuentra ilustrado en
Figura 52, del cual se puede identificar un lazo de control cerrado. En el cual la variable sujeta a control
es la temperatura del líquido (agua) que fluye en el intercambiador de calor y el elemento de control es
el vapor que fluye al interior de la tubería que pasa por el intercambiador de calor.
Figura 52
Nota. Adaptado de Instrumentación Industrial (p. 482), por Creus, 2010, Alfaomega Marcombo.
114
proceso, en el cual un operario detecta de manera sensorial la temperatura de salida del agua y según
su crítica procede a accionar la válvula de vapor para mantener condiciones deseadas de temperatura
en el agua de salida. Cabe destacar que debido a la inercia del proceso térmico no se obtienen
válvula, por lo que se necesita un cierto tiempo para realizar la calibración con la válvula del equipo.
Figura 53
Nota. Adaptado de Instrumentación Industrial (p. 483), por Creus, 2010, Alfaomega Marcombo.
corrección, se lo interpreta como una cadena cerrada de acciones, las cuales se las realiza múltiples
veces por el operador hasta obtener un resultado esperado. Este conjunto de operaciones dentro de un
Todos los procesos, según Creus (2010) presentan dos características fundamentales que deben
b) “el tiempo necesario para que la variable del proceso alcance un nuevo valor al ocurrir
un cambio de carga. Este retardo se debe a una o a varias propiedades del proceso:
El autor Ogata (2010), Dinámica de sistemas (p. 494) clasifica a los sistemas en lineales y no
lineales, dentro los cuales los de menor complejidad serán aquellos lineales representados por
Sistemas Lineales. Los sistemas lineales son analizados principalmente por sus propiedades de
linealidad a las cuales se puede aplicar el principio de superposición, el cual establece que “la respuesta
producida por la aplicación simultánea de dos funciones de excitación diferentes o entradas, es la suma
de dos respuestas individuales” (Ogata, 2010). Lo que produce que, dentro de un sistema lineal, su
respuesta a varias entradas puede computarse o calcularse analizando individualmente cada función de
entrada y posteriormente se suman los resultados, logrando que ecuaciones diferenciales con
Sistemas No Lineales. Por otro lado, los sistemas no lineales son aquellos interpretados
mediante ecuaciones no lineales, lo cual hace que el principio de superposición no sea aplicable.
116
Generalmente las ecuaciones físicas aplicadas a sistemas resultan ser lineales, no obstante, la
Debido a la dificultad para resolver matemáticamente los sistemas no lineales, se realiza una
linealización de los sistemas considerando una condición de operación. Ya conseguido un modelo lineal
a partir de un sistema no lineal, se aplican términos lineales para realizar el análisis y el diseño.
Dinámica de Sistemas
Según Ogata (2010), Dinámica de sistemas (p. 494), un sistema se define como “… una
combinación de componentes que actúan conjuntamente para alcanzar un objetivo específico”, lo cual
se puede aplicar a distintos ámbitos, expandiéndose a fenómenos dinámicos varios tales como el
En dinámica de sistemas se conoce como sistema dinámico a todo aquel en el cual su salida en
tiempo presente depende de una entrada o estímulo situado en un tiempo pasado, mientras que, un
sistema estático su salida depende únicamente de la entrada presente. Entonces un sistema estático
provee una salida constante si la entrada no posee variación alguna y caso contrario variaría su señal de
salida si la entrada varía. Dentro de un sistema dinámico la salida va a cambiar en función del tiempo
“La dinámica de sistemas trata del modelado matemático y el análisis de la respuesta de los
sistemas dinámicos” (Ogata, 2010) Toda interpretación de la realidad intenta explicarse bajo el
modelamiento matemático, es decir que casi cualquier fenómeno es susceptible a ser modelado bajo
conceptos matemáticos, lo cual lleva a la dinámica de sistemas a identificar y traducir los fenómenos o
condiciones de diseño en ecuaciones, ya sean sencillas como ecuaciones lineales o más complejas y
Elaboración de modelos. Un modelo matemático describe un sistema mediante sus leyes físicas,
que puede incluir parámetros desconocidos, que se pueden obtener realizando pruebas reales del
Es importante tener un modelo que sea lo más simplificado posible, cuando no se necesite una
ignorar estas variables, sus efectos, tienen respuestas pequeñas, los estudios del análisis del sistema
2. Escribir ecuaciones físicas en cada elemento, que deben ser combinados en base al esquema
3. Comprobar con resultados experimentales, una vez resuelto el modelo se verifica su validez.
Este proceso es repetitivo, hasta que exista una similitud entre las predicciones y los resultados
experimentales.
sistema cuyo modelo matemático es conocido. Una vez elaborado el modelo matemático del sistema
completo, que se constituye al unir los modelos matemáticos de cada componente se puede realizar el
análisis del mismo, cuando los parámetros del sistema varíen, puedan provocar varias soluciones, que
deben ser interpretados por el diseñador o también por la computadora (Ogata, 2010).
En el desarrollo para que un sistema cumpla con una acción especifica, se aplica el diseño que
matemático de manera que el procedimiento diseñado sea explicito para que un sistema funcione
Procedimientos de diseño. Inicialmente, es importante conocer los objetivos que tienen que
satisfacerse junto con sus especificaciones y también de la dinámica de sus componentes que se
introduce en el diseño. A medida que el diseño va tomando forma se puede notar que las
Ahora es importante resolver el diseño matemático del sistema que pone en evidencia las
soluciones del problema de diseño, es preferible, una vez halladas las soluciones del problema hacer una
simulación en computadora que nos demuestre el efecto producido por las entradas en el sistema
derivado. En la simulación, se deben tener soluciones satisfactorias para que el diseño sea considerado,
caso contrario, se deberá rediseñar hasta que se cumplan las soluciones. A partir de ese instante se
muestra en Figura 54 de manera que el sistema no llegue a la inestabilidad, se tiene también los valores
Figura 54
Nota. Adaptado de Sintonía experimental de controladores (p. 128), por A. Rosillo, 2008.
119
Al encontrar una relación del primer y segundo sobrepaso, que tiene un valor de 0.25, que se obtiene
cuando el valor del factor de amortiguamiento en los polos dominantes corresponde a 0.2176, mostrado
Figura 55
Nota. Adaptado de Sintonía experimental de controladores (p. 128), por A. Rosillo, 2008.
Tabla 9
Tipo de
Kp 𝝉i 𝝉d
Controlador
P 𝐾𝐶 - -
PI 𝐾𝐶 𝜏0 -
𝜏0 𝜏0
PID 𝐾𝐶
1.5 6
120
Modos de control
Considerando el ejemplo anterior, mantener un operador para cada lazo de control resultaría
algo absurdo, por lo cual las acciones o propósitos desarrolladas por dicho operador son realizadas por
En los distintos sistemas industriales existentes hoy en día se utilizan uno o varios sistemas de
a) De dos posiciones
c) Flotante
d) Proporcional (P)
control final se trabaja rápidamente entre una de dos posiciones fijas para un valor específico de la
variable a controlar, tal como se muestra en Figura 56. Además, se identifica una banda diferencial o
zona neutra en la cual el elemento de control permanece en una sola posición para un rango de valores,
Figura 56
Nota. Adaptado de Instrumentación Industrial (p. 488), por Creus, 2010, Alfaomega Marcombo.
Control Flotante
Posee similitud con el control todo-nada, pero mueve el elemento de control final a una
mediante una modificación o acoplamiento de un control todo-nada con una zona neutra con una
válvula motorizada de baja velocidad reversible, logrando así que la válvula permanezca inmóvil si la
variable objetivo se localiza dentro de la zona neutra, y al momento que dicho valor rebasa los límites la
válvula genera un movimiento en la dirección adecuada para modificar nuestra variable objetivo dentro
Figura 57
Nota. Adaptado de Instrumentación Industrial (p. 489), por Creus, 2010, Alfaomega Marcombo.
Al compararlo con un control todo-nada, el control de tipo flotante tiende a producir oscilaciones de la
variable controlada, las cuales pueden ser reducidas escogiendo adecuadamente velocidades del
Posee una ventaja frente al control todo-nada, la cual radica en que puede compensar cambios
importante o los cambios de carga son muy rápido no es recomendable utilizar este tipo de control.
Control Proporcional
Es el control en el cual la válvula varía su posición en la misma magnitud por cada unidad de
desviación, es decir existe una relación lineal continua con el valor de la variable sujeta a control y el
posicionamiento del elemento de control final, todo esto en el rango de la banda proporcional.
proporcional, en el cual su punto de consigna radica en 150 °C y su intervalo en el que actúa va desde los
100 hasta los 200 °C. Se observa además que mientras la temperatura se mantiene entre los 100 y 200
control, mientras que cuando está a 100 °C o menos la válvula se encuentra totalmente abierta y caso
Figura 58
Nota. Adaptado de Instrumentación Industrial (p. 490), por Creus, 2010, Alfaomega Marcombo.
De modo que mediante un control de tipo proporcional su salida de control está definida de la
siguiente manera:
Laplace:
Un control integral según (Hernández, 2010) “produce una acción de control que es
proporcional a la integral del error en el tiempo”, lo que significa que realiza la suma de los errores
cambio de los valores de salida del control 𝐼 es directamente proporcional a los valores de la señal de
error en la entrada
𝑑𝐼/𝑑𝑡 = 𝐾𝐼 ∙ 𝑒 ( 63 )
Entonces al ingresar una señal de error constante al controlador se obtiene una función de
salida de controlador del tipo recta pendiente, tal como se muestra en Figura 59. De donde se puede
interpretar que, al momento de tener una salida del controlador constante, la señal de error al ingreso
será resultante cero y cuando la salida del controlador posee una variación de manera constante la señal
Figura 59
Dentro del control industrial un modelo de control de tipo integral generalmente no se lo utiliza
de manera individual, sino se lo utiliza en combinación con el modo de control proporcional, resultando
Figura 60
Obteniendo como resultado un comportamiento ilustrado en Figura 60, donde se compara dos
𝐾𝐼
𝑓𝑢𝑛𝑐𝑖ó𝑛 𝑑𝑒 𝑡𝑟𝑎𝑛𝑠𝑓𝑒𝑟𝑒𝑛𝑐𝑖𝑎 = 𝐾𝑝 + ( 65 )
𝑠
El controlador PID al ser una combinación de los tres tipos de modos de control (proporcional,
integral y derivativo) ofrece la ventaja de no poseer una desviación en el error y a su vez decrece la
𝑑𝑒
𝑠𝑎𝑙𝑖𝑑𝑎 𝑑𝑒𝑙 𝑐𝑜𝑛𝑡𝑟𝑜𝑙𝑎𝑑𝑜𝑟 = 𝐾𝑝 ∙ 𝑒 + 𝐾𝐼 ∫ 𝑒 𝑑𝑡 + 𝐾𝐷 ( 66 )
𝑑𝑡
1
𝑓𝑢𝑛𝑐𝑖ó𝑛 𝑑𝑒 𝑡𝑟𝑎𝑛𝑠𝑓𝑒𝑟𝑒𝑛𝑐𝑖𝑎 = 𝐾𝑝 + 𝐾 + 𝑠 𝐾𝐷 ( 67 )
𝑠 𝐼
127
Capítulo III. Diseño y Construcción del Banco de Pruebas del Sistema De Refrigeración
por su buena conductividad térmica y además de su resistencia a la corrosión, que son características
principales para este tipo de sistemas (Muñoz & Zarantonelli, 2017). El tipo de tubo a utilizar es el tipo L,
que se detalla por su funcionamiento en Tabla 10 ya que es recocido y debido a su alta plasticidad que
posibilita su enrollamiento para que pueda ser doblado y facilite su manufactura, tomando en cuenta
que a nivel industrial esta tubería tiene una longitud de 15 metros (Scalofrios, 2013).
Tabla 10
Existen varios criterios generales de diseño para seleccionar una tubería adecuada de
refrigeración para tener una instalación optima, las líneas deben ser diseñadas lo mas cortas posibles y
Por otro lado, es necesario saber que el aceite que circula por las tuberías para lubricar y enfriar
al compresor debe ser trasladado por las líneas con una velocidad mínima de 5 m/s en tramos
Tabla 11
Parámetros Cantidad
Temperatura de evaporización 0 °𝐶
Temperatura de condensación 55 °𝐶
Tabla 12
Parámetros Cantidad
Nota. Adaptado de Nueva Enciclopedia de la Climatización (p. 301), por J. Ramírez, CEAC.
En donde de la Tabla 12, se selecciona que las pérdidas de carga serán admitidas a 1°C, cómo se
Figura 61
Nota. Adaptado de Nueva Enciclopedia de la Climatización (p. 300), por J. Ramírez, CEAC.
Para evaluar las pérdidas en una longitud de 10m se realiza la siguiente corrección en la caída de
presión. Tomando en cuenta que la tubería que se estableció en promedio donde se ubicara el sistema
0.1 ∗ 10
𝐽= = 0.5 𝑏𝑎𝑟
2
Tabla 13
Tabla 14
Velocidades en m/s
R717 15 a 20 15 𝑎 20
R134a
R22 4.5 a 20 10 a 18
R404a
Nota. Adaptado de Nueva Enciclopedia de la Climatización (p. 300), por J. Ramírez, CEAC.
En base a este criterio se considera evaluar en función del flujo necesario de refrigerante dentro
del sistema de 0.036 Kg/s, para obtener una velocidad adecuada en función de diferentes diámetros
Tabla 15
Velocidad del refrigerante dentro del sistema para diferentes diámetros de tubería
diferentes diámetros de tubería se tienen distintas velocidades, por lo que, el primer criterio de
selección de tubería se toma en consideración que para el refrigerante R404A, debe circular en las líneas
de aspiración a una velocidad de 4.5 a 20 m/s lo que se decide seleccionar tubería de 3/8” y para la línea
de descarga se selecciona la tubería de ½”, con el fin de que el refrigerante se mantenga en el rango
aspiración o succión con el evaporador, en primera instancia se seleccionan los diámetros de 3/8” para
Además, existe otro parámetro de selección en base al área de transferencia de calor como se
menciona a continuación.
Una vez seleccionado el tipo de tubería para el sistema, es importante que bajo características
de diseño se tomó la tubería de cobre de 3/8” ya que en esta tubería las perdidas por presión son
mínimas, debido a que la distancia entre la válvula de expansión y el evaporador debe ser pequeña,
Para el caso de la tubería del condensador, se tomó en cuenta que generalmente las áreas de
transferencia de calor son grandes y por ende se necesita una tubería de diámetro mayor, al considerar
seleccionar el tipo de refrigerante, aquí se toma en cuenta que en el evaporador se desea que circule
vapor sobrecalentado. Por tal motivo, se necesita que el punto de ebullición sea menor a la temperatura
evaporación sea aproximadamente 10°C menor. En este caso 10°C menor a la del agua en el reservorio
donde se encuentra el evaporador. Para definir esto se considera que la temperatura a la que va a llegar
el agua es de 10°C, por lo tanto, para este criterio de la Tabla 4, se seleccionan los refrigerantes cuyo
Tabla 16
22 -41 Si Cumple
141B 32 No cumple
Como se puede apreciar en Tabla 16, el único refrigerante que no cumple con estas condiciones
Seguridad. La seguridad es otro criterio que permite definir el tipo de refrigerante a usar, para
este diseño se toma en cuenta la salud de las personas, debido a que este equipo será manipulado por
estudiantes dentro del laboratorio de conversión de energía. Por lo tanto, se considera una categoría de
seguridad A1, con el fin de que el refrigerante a usar tenga baja propagación a la llama, que no sea
Para este análisis se parte desde el criterio anterior, descartando el refrigeranteR-141B, como se
Tabla 17
De la Tabla 17, se puede observar que existen refrigerantes como el R-290, R-600 A y R-406 A,
que se ubican en la categoría A3 y A2, por lo tanto, no son tóxicos, pero recaen la categoría de alta y
baja propagación a la llama, siendo negativos para este diseño y por lo tanto se descartan.
normativas donde mencionan que el uso de los refrigerantes o sustancias químicas deben ser
Como se mencionó en el Capítulo II, Existen indicadores para medir el impacto ambiental de los
El potencial de agotamiento de ozono (PAO), es el indicador que permite medir la capacidad que
tiene un refrigerante para perjudicar la capa de ozono, los valores están en el rango de 0 y 1.
Recordando que se tiene al refrigerante CFC-11 como referencia, teniendo un valor de 1, por lo tanto, es
generen gases de efecto invernadero, acelerando el calentamiento global. Recordando que los
Por lo tanto, el criterio de selección de los refrigerantes se toma un valor de PAO igual a cero y
para valores de PCA no muy elevados ya que mayor será el impacto atmosférico En Tabla 18 se
Tabla 18
R-22 y R-409 A, puesto a que su valor es mayor a 0 y para minimizar el efecto perjudicial a la capa de
Seguidamente, entre los refrigerantes que han cumplido los criterios anteriormente expuestos,
mediante la aplicación de matrices de decisión se selecciona el refrigerante R 404 a (Ver Apéndice 4).
136
Selección de Compresor
dimensionamiento energético.
Por lo tanto, como se puede observar en la Tabla 19, se procedió a calcular la energía calorífica
necesaria (Ver ecuación 21) en el agua para que un flujo másico, que será variado entre 0.6 a 0.05 kg/s,
aumente una determinada cantidad de temperatura (Δ𝑇). Además, la elección de los parámetros
deberá cumplir con una energía calorífica necesaria, dentro del rango de 1 Tonelada de refrigeración,
Tabla 19
En Figura 62, se muestra la variación de los parámetros con respecto a la energía necesaria para
Figura 62
55
50
45
40
35
Potencia
30
(kW)
25 ∆T: 20 C
20
∆T: 10 C
15
10 ∆T: 5 C
5
0
0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 0.65
Flujo Másico
(kg/s)
Luego de observar los resultados obtenidos, se presenta en Tabla 20, un resumen de los
resultados más cercanos a cumplir con el diseño previsto de 1 Tonelada de refrigeración, siendo este
Tabla 20
∆T = 10 C
20 10 4185 0.60 25110 25.11
20 10 4185 0.50 20925 20.93
20 10 4185 0.40 16740 16.74
20 10 4185 0.30 12555 12.56
20 10 4185 0.20 8370 8.37
20 10 4185 0.10 4185 4.19
20 10 4185 0.08 3348 3.35
20 10 4185 0.05 2093 2.09
balance de energía enfocado al intercambiador de calor externo al banco de pruebas. Para esto, se
𝐸̇𝑒𝑛𝑡𝑟𝑎𝑑𝑎 = 𝐸̇𝑠𝑎𝑙𝑖𝑑𝑎
La ecuación 68, permite calcular el flujo másico requerido por el aceite para determinada
Tabla 21
Temperatura
ENERGIA UTILIZADA PARA
(C) Calores Específicos Flujo Másico
DISMINUIR UN ∆T DE ACEITE
(J/kg K) (kg/s)
ENTRA SALE (kW)
35 45 1950.00 0.172
40 60 2006.00 0.08
En Tabla 22, se muestra un resumen de los parámetros de diseño finales para los fluidos agua y
aceite. Parámetros que fueron calculados y seleccionados a partir de Tabla 20 y Tabla 21.
Tabla 22
Temperatura
Flujo Másico (C) Se pretende enfriar y calentar
Fluido
(kg/s) respectivamente:
ENTRA SALE
dicha capacidad con respecto al fluido que se verá inmerso en el evaporador, se puede proceder a la
selección de la misma.
141
Aplicación de factor de seguridad del 10% según Norma ASHRAE (Pag 248)
del agua tiene que ser mayor a la del evaporador para poder perder calor y que la temperatura de
entrada del aceite menor a la del condensador para poder ganar calor (Barletta & Acevedo, 2021). De
acuerdo a los parámetros, tanto de agua como de aceite, se realiza la selección de temperaturas de
evaporación y condensación del compresor. Sabiendo que la temperatura de entrada del agua tiene que
ser mayor a la del evaporador para poder perder calor y que la temperatura de entrada del aceite
menor a la del condensador para poder ganar calor (Barletta & Acevedo, 2021).
Temperatura de condensación:
𝑇𝑐 = 55 °𝐶
Temperatura de evaporación:
𝑇𝑒 = 0 °𝐶
Temperatura de sobrecalentado:
𝑇𝑠𝑜𝑏𝑟𝑒𝑐𝑎𝑙𝑒𝑛𝑡𝑎𝑚𝑖𝑒𝑛𝑡𝑜 = 7 °𝐶
Temperatura de subenfriamiento:
𝑇𝑠𝑢𝑏𝑒𝑛𝑓𝑟𝑖𝑎𝑚𝑖𝑒𝑛𝑡𝑜 = 11 °𝐶
142
realizarán los cálculos necesarios tanto para, refrigerante 404a como refrigerante 134a, de manera que
Cálculos de unidad compresora para refrigerante R 404 A. Con los datos obtenidos se procede
al dimensionamiento del compresor para su posterior selección. Esto se realizará con la ayuda del
Figura 63
Nota. Adaptado de Handbook (p. 756), por ASHRAE, 2013, American Society of Heating, Refrigerating
Con el cálculo de propiedades termodinámicas en cada uno de los puntos del ciclo de
𝑄̇𝑒 = ℎ1𝑠 − ℎ4
𝐾𝐽
𝑄̇𝑒 = 103
𝐾𝑔
𝑄̇𝑐 = ℎ2 − ℎ3
𝐾𝐽
𝑄̇𝑐 = 134
𝐾𝑔
𝐾𝐽
𝑊̇𝑐𝑜𝑚𝑝 = 31
𝐾𝑔
COP de refrigeración
ℎ1 − ℎ4 372 − 269
𝐶𝑂𝑃𝑅𝑒𝑓 = =
ℎ2 − ℎ1 403 − 372
𝐶𝑂𝑃𝑅𝑒𝑓 = 3.3225
144
La producción especifica neta denota la cantidad de energía a absorber por unidad de masa,
por lo que su relación con la producción frigorífica requerida ayudará a encontrar el flujo másico de
𝐾𝐽 3600 𝑠
𝑄𝑇 = 3.685 ∙
𝑠 1ℎ
𝐾𝐽
𝑄𝑇 = 13266
ℎ
𝐾𝐽
13266 ⁄ℎ
𝑚𝑟 =
𝐾𝐽
103 ⁄𝐾𝑔
𝐾𝑔 𝑘𝑔
𝑚𝑟 = 128.80 = 0.036
ℎ 𝑠
El compresor deberá estar capacitado para recibir una cantidad de 128.80 𝑘𝑔/ℎ a la entrada del
su masa, por lo que se calcula el volumen de refrigerante aspirado por el compresor a su entrada.
𝑉𝑎 = 𝑚𝑟 ∙ 𝑉1𝑠
𝐾𝑔 𝑚3
𝑉𝑎 = 128.80 ∙ 0.03448
ℎ 𝐾𝑔
𝑚3
𝑉𝑎 = 4.44
ℎ
Para posteriormente calcular el volumen teórico que el compresor deberá ser capaz de
𝑉𝑎
𝜂𝑉 =
𝑉𝑏
𝑃𝐾 2.5
𝑟𝑝 = = = 4.16
𝑃𝑜 0.6
𝜂𝑉 = 1 − 0.05(𝑟𝑝 )
𝜂𝑉 = 1 − 0.05(4.16)
𝜂𝑉 = 0.79
𝑉𝑎
𝑉𝑏 =
𝜂𝑉
3
4.44 𝑚 ⁄ℎ
𝑉𝑏 =
0.79
3
𝑉𝑏 = 5.62 𝑚 ⁄ℎ
𝑃𝑡 = (ℎ2 − ℎ1𝑠 ) ∙ 𝑚𝑅
𝐾𝐽 𝐾𝑔 1ℎ
𝑃𝑡 = 31 ∗ 128.80 ∗
𝐾𝑔 ℎ 3600
𝑃𝑡 = 1.109 𝑘𝑊 = 1.49𝐻𝑃
146
Finamente, en Tabla 23, se detallan los cálculos obtenidos para el refrigerante R 404 a.
Tabla 23
Volumen de Potencia
Flujo Másico Relación de Eficiencia
Refrigerante barrido (𝒎𝟑 / requerida
(𝒌𝒈/𝒔) Compresión Volumétrica
𝒉) (HP)
R 404 a 0.036 4.16 0.79 5.62 1.49
El tipo de válvula de expansión a utilizar será la Válvula de Expansión Electrónica (VEE), la cual se
Una válvula de expansión se puede seleccionar bajo los criterios de la capacidad nominal de
Figura 64
Nota. Adaptado de Sanhua Series DPF-T/S Válvula de Expansión Electrónica (p. 2), por Sanhua, 2012.
Debido a que la capacidad nominal del compresor, es decir 1.5 HP, lo cual nos da 11.19 Kw y
Las dimensiones que nos ayudarán a saber las especificaciones de la válvula con los diámetros de
Figura 65
Nota. Adaptado de Sangua Series DPF-T/S Válvula de Expansión Electrónica (p. 3), por Sanhua, 2012.
Para enfriar agua y calentar aceite mediante el uso de un sistema de refrigeración, es necesario
que exista un intercambio de calor entre el fluido externo y el refrigerante que es el encargado de quitar
o ceder calor del medio en el que se encuentra. Por esa razón el evaporador y condensador a diseñar
combinación de convección y conducción por medio de la pared del tubo que divide a los dos fluidos.
En este caso se produce convección en una sola fase del fluido que circula por la coraza del
intercambiador y convección con cambio de fase del refrigerante que circula por el interior del
serpentín.
Por lo que, se requiere diseñar el evaporador para enfriar agua empleando refrigerante R404A
que circula por el interior del serpentín, con la finalidad de que el agua que ha perdido calor recircule
nuevamente por el sistema. En Tabla 24 se registra el rango de temperaturas y flujos másicos tanto para
Tabla 24
Temperatura de
20 0 °𝐶
entrada
Temperatura de
10 7 °𝐶
salida
Figura 66
transferencia de calor, para después determinar el área de transferencia de calor mediante la ecuación
de Fourier.
El esquema geométrico inicial del intercambiador de calor serpentín coraza se muestra en Figura
de calor de serpentín helicoidal desarrollado por R. Patil y B. Shende y R. Chosh. (Patil et al., 1982)
151
Figura 67
Parámetros Geométricos
En Figura 68, se muestran las dimensiones de la tubería de cobre tipo L, la medida nominal del
Figura 68
Nota. Adaptado de Manual Técnico de Cobre (p. 10), por Nacobre S.A.
152
Tabla 25
𝑑𝑜 = 0.0127 𝑚
𝑡𝑡𝑢𝑏𝑜 : 0.000889 𝑚
𝑑𝑖 = 𝑑𝑜 − 2 𝑡𝑡𝑢𝑏𝑜 = 0.01092 𝑚
𝐷ℎ = 0.3 𝑚
𝐷ℎ𝑒 = 𝐷ℎ + 𝑑𝑜 = 0.313 𝑚
𝐷ℎ𝑖 = 𝐷ℎ − 𝑑𝑜 = 0.287 𝑚
Paso de la hélice:
𝑝 = 1.5 𝑑𝑜 = 0.019 𝑚
153
𝐷ℎ 2
𝐿𝑠𝑒𝑟𝑝𝑒𝑛𝑡í𝑛 = 𝑁𝑣 √(2𝜋 ) + 𝑝2
2
𝜋 2
𝑉𝑠𝑒𝑟𝑝𝑒𝑛𝑡í𝑛 = 𝑑 𝐿
4 𝑜 𝑠𝑒𝑟𝑝𝑒𝑛𝑡í𝑛
𝜋
𝑉𝑠𝑒𝑟𝑝𝑒𝑛𝑡í𝑛 = (0.0127𝑚)2 0.943 𝑁𝑣 𝑚 = 0.0001194 𝑁𝑣 [𝑚3 ]
4
Volumen anular:
𝜋
𝑉𝑎𝑛𝑢𝑙𝑎𝑟 = (𝐷𝑖2 ) 𝑝 𝑁𝑣
4
𝜋
𝑉𝑎𝑛𝑢𝑙𝑎𝑟 = ((0.4𝑚)2 ) 0.019 𝑚 (𝑁𝑣 ) = 0.002 𝑁𝑣 [𝑚3 ]
4
El fluido que circula por el interior del serpentín experimenta un cambio de fase durante su
recorrido, el refrigerante ingresa al serpentín y en el transcurso que circula por el tubo cambia a estado
𝑇𝑟4 + 𝑇𝑟1
𝑇𝑝𝑟𝑒 =
2
154
0+7
𝑇𝑝𝑟𝑒 = = 3.5 °𝐶
2
del tubo, ya que no se tiene información inicial de los coeficientes de convección interno y externo, así
que una buena aproximación puede realizarse entre la temperatura que está dentro del tubo y la
Tomando en cuenta que la temperatura 𝑇𝑠 debe ser mayor que la 𝑇𝑠𝑎𝑡 para que se produzca la
evaporación, se tiene:
0 + 7 20 + 10
+
𝑇𝑠 = 2 2 = 9.25 °𝐶
2
𝑇𝑠 > 𝑇𝑠𝑎𝑡
𝑇𝑠𝑎𝑡 + 𝑇𝑠
𝑇𝑓 = = 4.625 °𝐶
2
En Tabla 26 y Tabla 27, se presentan las propiedades termo físicas del refrigerante evaluadas a
tabla de propiedades termo físicas para refrigerantes de la norma ASHRAE Capítulo 30 (Apéndice 1)
155
Tabla 26
Tabla 27
El calor especifico de vapor se toma a una temperatura promedio entre la temperatura de vapor
𝑇𝑟1 = 7 °𝐶
𝑇𝑟1 + 𝑇𝑠𝑎𝑡
𝑇𝑝𝑣𝑟 = = 3.5 °𝐶
2
𝑘𝐽
𝐶𝑝𝑣𝑅 = 1.02792
𝑘𝑔 𝐾
156
Se calcula la velocidad con al que circula el refrigerante por el interior del tubo:
𝑚𝑟 𝑚3
𝑄𝑟 = = 3.15 𝑥 10−5
𝜌𝑙𝑟 𝑠
𝑄𝑟 𝑚
𝑉𝑙𝑟 = 2 = 0.337 𝑠
𝑑
𝜋 𝑖
4
Número de Reynolds
𝜌𝑙𝑟 𝑉𝑙𝑟 𝑑𝑖
𝑅𝑒𝑟 = = 24 479.26
𝜇𝑙𝑟
Para encontrar el coeficiente de convección interno, se usa relación la descrita por Chaton
considerando que empiezan a generarse las primeras burbujas de vapor, ecuación 43.
1/4
9.81 𝜌𝑙𝑟 (𝜌𝑙𝑟 − 𝜌𝑣𝑟 ) (𝑘𝑙𝑟 )3 3
ℎ𝑖 = 0.555 ( (ℎ𝑓𝑔 ∙ 1000 + ∙ 𝐶𝑝𝑙𝑟 ∙ 1000 ∙ ((𝑇𝑠 − 𝑇𝑠𝑎𝑡 ))))
𝜇𝑙𝑟 (𝑇𝑠 − 𝑇𝑠𝑎𝑡 ) 𝑑𝑖 8
𝑊
ℎ𝑖 = 1 506.95
𝑚2 𝐾
Para calcular el coeficiente de convección externo se toma en cuenta que el fluido de trabajo
que circula por la coraza del evaporador es agua y no experimenta un cambio de fase.
Tabla 28
Diámetro equivalente
4 𝑉𝑓𝑙𝑢𝑗𝑜
𝐷𝑒𝑞 =
𝜋 𝑑𝑜 𝐿𝑠𝑒𝑟𝑝𝑒𝑛𝑡í𝑛
4 (0.002 𝑁𝑣 [𝑚3 ])
𝐷𝑒𝑞 = = 0.242 𝑚
𝜋 (0.0127 [𝑚]) (0.943 𝑁𝑣 [𝑚])
Velocidad de masa
𝑚𝑎
𝐺𝑎 = 𝜋
2 2 2
4 ((𝐷𝑖 )) − (𝐷ℎ𝑒 − 𝐷ℎ𝑖 ))
158
𝑘𝑔
0.08 𝑠 𝑘𝑔
𝐺𝑎 = 𝜋 = 0.704 2
2) ((0.313𝑚)2 − (0.287𝑚)2 ) 𝑚 𝑠
4 ((0.4𝑚) −
Número de Reynolds
𝐷𝑒𝑞 𝐺𝑎
𝑅𝑒𝑎 =
𝜇𝑎
𝑘𝑔
0.242 𝑚 (0.704 )
𝑚2 𝑠
𝑅𝑒𝑎 = = 149.568
𝑘𝑔
1.138 ∗ 10−3
𝑚𝑠
Debido al bajo número de Reynolds 149.568 < 2300, se considera que el agua se encuentra
en régimen laminar y se decide ocupar la expresión descrita por Chilton, Colburn, Gneraux y Vernons,
Para reemplazar los términos en ecuación 40, se debe tomar en cuenta las unidades en las que
se trabaja, por lo que en Tabla 29 se presenta las propiedades termofísicas en unidades inglesas y la
Tabla 29
∆𝑡 = 𝑇𝑎1 − 𝑇𝑝𝑠𝑒𝑟𝑝𝑒𝑛𝑡í𝑛
0.25
0.343 (62.3742 )(1) (7.667𝑥10−5 ) 51.35
ℎ𝑜 = 116 (( )( ))
1.138 0.5
𝐵𝑡𝑢 𝑊
ℎ𝑜 = 117.728 2
= 668.458 2
ℎ 𝑓𝑡 ℎ 𝑚 𝐶
Para el cálculo del coeficiente global de transferencia se calcula con la sumatoria de resistencias
1
𝑈𝑒 =
𝑑𝑜 d ln (𝑟𝑜 /𝑟𝑖 ) 1
+ o +
𝑑𝑖 ℎ𝑟 𝐾𝑐𝑢 ℎ𝑎
1
𝑈𝑒 =
0.0127 0.0127 ln (0.0127/0.01092) 1
+ +
(0.01092) (1996) 385 668.458
𝑊
𝑈𝑒 = 440.03
𝑚2 𝐾
se reemplazan los valores de las temperaturas al ingreso y a la salida del evaporador tanto como al
(20 − 7) − (10 − 0)
∆𝑇 = = 11.43 °𝐶
20 − 7
ln (10 − 0)
𝑄 = 𝑈 𝐴𝑜 ∆𝑇
𝑄
𝐴𝑜 =
𝑈𝑒 ∆𝑇
Por lo tanto, reemplazando la capacidad térmica para enfriar el agua, se tiene que el área de
3348
𝐴𝑜 = = 0.665 𝑚2
(479.937)(11.43)
El número de espiras queda expresado en función de la longitud del serpentín como se muestra
en la siguiente expresión:
𝐴𝑜
𝑁𝑣 =
𝐿𝑠𝑒𝑟𝑝𝑒𝑛𝑡í𝑛
𝜋 𝑑𝑜 ( )
𝑁𝑣
0.61
𝑁𝑣 = = 17.692 𝑣𝑢𝑒𝑙𝑡𝑎𝑠
0.943 𝑁𝑣
𝜋 (0.0127) ( 𝑁 )
𝑣
𝐻𝑒 = 𝑁𝑣 𝑝 + 𝑑𝑜
rango ya que se requiere que el equipo sea capaz de trabajar en las condiciones diseñadas. La longitud
𝐿𝑒 = 𝑁𝑣 𝐿𝑠𝑒𝑟𝑝𝑒𝑛𝑡í𝑛 1.15
𝐿𝑒 = 17.692(0.943)(1.15) = 19.179 𝑚
sumergido. En este caso se emplea para calentar aceite con la finalidad de que recircule nuevamente
por el sistema.
Tabla 30
Temperatura de entrada 35 68 °𝐶
Temperatura de salida 45 44 °𝐶
Figura 69
El esquema geométrico inicial del intercambiador de calor serpentín coraza se muestra en Figura
70. Los parámetros geométricos se calculan en base al procedimiento de diseño del evaporador.
Figura 70
Parámetros Geométricos
En Figura 71, se muestran las dimensiones de la tubería de cobre tipo L, la medida nominal del
tubo para el serpentín del condensador es de 1/2”. En Tabla 31 se presenta las dimensiones geométricas
Tabla 31
Dimensiones geométricas
Nota. Adaptado de Manual Técnico de Cobre (p. 10), por Nacobre S.A.
𝑑𝑜 = 0.015878 𝑚
𝑡𝑡𝑢𝑏𝑜 : 0.001016 𝑚
164
𝑑𝑖 = 𝑑𝑜 − 2 𝑡𝑡𝑢𝑏𝑜 = 0.01385 𝑚
𝐷ℎ = 0.4 𝑚
𝐷ℎ𝑒 = 𝐷ℎ + 𝑑𝑜 = 0.416 𝑚
𝐷ℎ𝑖 = 𝐷ℎ − 𝑑𝑜 = 0.384 𝑚
Paso de la hélice
𝑝 = 1.25 𝑑𝑜 = 0.02 𝑚
𝜋
𝑉𝑠𝑒𝑟𝑝𝑒𝑛𝑡í𝑛 = (0.015878𝑚)2 1.257 𝑁𝑣 𝑚 = 0.0002489 𝑁𝑣 [𝑚3 ]
4
Volumen anular
𝜋
𝑉𝑎𝑛𝑢𝑙𝑎𝑟 = ((0.5𝑚)2 ) 0.02 𝑚 (𝑁𝑣 ) = 0.004 𝑁𝑣 [𝑚3 ]
4
El fluido que circula por el interior del serpentín experimenta un cambio de fase durante su
recorrido, el refrigerante ingresa al serpentín como vapor sobrecalentado y en el transcurso que circula
165
por el tubo cambia a estado líquido, es decir, se produce un proceso de condensación y se analiza este
68 + 44
𝑇𝑝𝑟𝑒 = = 56 °𝐶
2
Tomando en cuenta que la temperatura 𝑇𝑠 debe ser menor que la 𝑇𝑐𝑜𝑛𝑑 para que se produzca
la condensación, se tiene:
68 + 44 35 + 45
+
𝑇𝑠 = 2 2 = 48 °𝐶
2
𝑇𝑠𝑎𝑡 + 𝑇𝑠
𝑇𝑓 = = 51.5 °𝐶
2
En las Tabla 32 y Tabla 33, se presentan las propiedades termofísicas del refrigerante evaluadas
Tabla 32
Tabla 33
El calor especifico de vapor se toma a una temperatura promedio entre la temperatura de vapor
𝑇𝑟2 = 68 °𝐶
𝑇𝑟2 + 𝑇𝑠𝑎𝑡
𝑇𝑝𝑣𝑟 = = 61.5 °𝐶
2
𝑘𝐽
𝐶𝑝𝑣𝑟 = 3.1312
𝑘𝑔 𝐾
167
Se calcula la velocidad con la que circula el refrigerante por el interior del tubo:
𝑚𝑟 𝑚3
𝑄𝑟 = = 4.055 𝑥 10−5
𝜌𝑙𝑟 𝑠
𝑄𝑟 𝑚
𝑉𝑣𝑟 = 2 = 0.269 𝑠
𝑑
𝜋 𝑖
4
Número de Reynolds
𝜌𝑣𝑟 𝑉𝑣𝑟 𝑑𝑖
𝑅𝑒𝑟 = = 30 824.595
𝜇𝑣𝑟
Para encontrar el coeficiente de convección interno, se usa relación la descrita por Chaton,
considerando que empiezan a generarse las primeras gotas de líquido refrigerante, ecuación 43.
1/4
9.81 𝜌𝑙𝑟 (𝜌𝑙𝑟 − 𝜌𝑣𝑟 ) (𝑘𝑙𝑟 )3 3
ℎ𝑖 = 0.555 ( (ℎ𝑓𝑔 ∙ 1000 + ∙ 𝐶𝑝𝑙𝑟 ∙ 1000 ∙ ((𝑇𝑠𝑎𝑡 − 𝑇𝑠 ))))
𝜇𝑙𝑟 (𝑇𝑠𝑎𝑡 − 𝑇𝑠 ) 𝑑𝑖 8
𝑊
ℎ𝑖 = 1 064.157
𝑚2 𝐾
Para calcular el coeficiente de convección externo se toma en cuenta que el fluido de trabajo
que circula por la coraza del evaporador es aceite y no experimenta un cambio de fase.
Tabla 34
Diámetro equivalente
4 (0.004 𝑁𝑣 [𝑚3 ])
𝐷𝑒𝑞 = = 0.233 𝑚
𝜋 (0.015878 [𝑚]) (1.257 𝑁𝑣 [𝑚])
Velocidad de masa
𝑘𝑔
0.17 𝑠 𝑘𝑔
𝐺𝑎𝑐 =𝜋 = 0.966 2
2) ((0.416𝑚)2 − (0.384𝑚)2 ) 𝑚 𝑠
4 ((0.5𝑚) −
Número de Reynolds
𝑘𝑔
0.233 𝑚 (0.966 )
𝑚2 𝑠
𝑅𝑒𝑎𝑐 = = 24.023
𝑘𝑔
9.364 ∗ 10−3 𝑚𝑠
169
Debido al bajo número de Reynolds 24.023 < 2300, se considera que el aceite se encuentra
en régimen laminar y se decide ocupar la expresión descrita por Chilton, Colburn, Gneraux y Vernons,
Tabla 35
∆𝑡 = 𝑇𝑠 − 𝑇𝑎𝑐1
0.25
0.0633 (54.1272 )(0.437) (3.889𝑥10−4 ) 55.4
ℎ𝑜 = 116 (( )( ))
9.364 0.625
𝐵𝑡𝑢 𝑊
ℎ𝑜 = 21.49 2
= 122.023 2
ℎ 𝑓𝑡 ℎ 𝑚 𝐶
170
Para el cálculo del coeficiente global de transferencia se calcula con la sumatoria de resistencias
1
𝑈𝐶 =
0.625 0.625 ln (0.625/0.01385) 1
+ + 122.023
(0.01385) (1 064.157) 385
𝑊
𝑈𝐶 = 107.78
𝑚2 𝐾
9, se tiene:
𝑄
𝐴𝑜 =
𝑈𝑐 ∆𝑇
Por lo tanto, reemplazando la capacidad térmica para calentar el aceite, se tiene que el área de
3348
𝐴𝑜 = = 2.082 𝑚2
(107.78)(14.92)
2.082
𝑁𝑣 = = 33.209 𝑣𝑢𝑒𝑙𝑡𝑎𝑠
1.257 𝑁𝑣
𝜋 (0.015878 ) ( 𝑁𝑣 )
171
rango ya que se requiere que el equipo sea capaz de trabajar en las condiciones diseñadas. La longitud
En Tabla 36 se enlista los parámetros finales del diseño de los serpentines tanto del evaporador
Tabla 36
Selección de Bombas
Tabla 37
Material PVC -
En Tabla 38, se presentan las características que poseerán los fluidos en el sistema de
alimentación.
Tabla 38
A 10 ℃ A 55 ℃
Densidad 𝑘𝑔/𝑚3
999.7 876
En Tabla 39, se presentan los datos para calcular las perdidas mayores en cada uno de los
sistemas de alimentación.
Tabla 39
𝑇𝑟𝑎𝑛𝑠𝑖𝑐𝑖ó𝑛 𝐿𝑎𝑚𝑖𝑛𝑎𝑟 -
Número de Reynolds
2039.94 26.06 -
En Tabla 40, se presentan los datos para calcular las perdidas menores en ambos sistemas de
alimentación.
Tabla 40
Accesorios Cantidad 𝑲𝑳 ∑ 𝑲𝑳
Codos (roscado) 5 0.9 4.5
Acople tipo T (roscado) 2 0.9 1.8
Reducción (1" - 3/4") 3 0.1 0.3
Unión (roscada) 1 0.08 0.08
Acoples Universales 2 0.5 1
Válvula de Bola, totalmente abierta 1 0.05 0.05
Válvula de Compuerta, 1/2" cerrada. 2 2.1 4.2
Total 11.93
Nota. Adaptado de Mecánica de Fluidos. Fundamentos y Aplicaciones (p. 350), por Y. Cengel, 2012, Mc
Graw Hill.
174
𝑃1 𝑉12 𝑃2 𝑉22 𝐿 𝑉2
+ + 𝑧1 + ℎ𝑏𝑜𝑚𝑏𝑎 = + + 𝑧2 + ℎ𝑡𝑢𝑟𝑏𝑖𝑛𝑎 + (𝑓 + ∑ 𝐾𝐿 )
𝜌𝑔 2𝑔 𝜌𝑔 2𝑔 𝐷 2𝑔
Se tiene como fin, calcular la carga que debe tener la bomba a seleccionar para el bombeo de
agua y aceite
Figura 72
Tabla 41
Para el punto 2 a la salida del intercambiador de calor, la presión (𝑃2 ) tendrá las siguientes
consideraciones: las cargas generadas por presiones serán insignificantes comparadas con las cargas
generadas por altura, y la presión manométrica en dicho punto será cercana a cero debido a los bajos
caudales a aplicarse en las prácticas de laboratorio. Entendiendo dichas consideraciones, se tienen las
Tabla 42
En Tabla 43 se muestran los resultados de la carga requerida por la bomba, tanto de aceite
como de agua, que son calculados aplicando la ecuación de la energía para sistemas de bombeo,
Ecuación 69.
𝑃2 − 𝑃1 𝑉22 − 𝑉12 𝐿 𝑉2
ℎ𝑏𝑜𝑚𝑏𝑎 = + + (𝑧2 − 𝑧1 ) + ℎ𝑡𝑢𝑟𝑏𝑖𝑛𝑎 + (𝑓 + ∑ 𝐾𝐿 ) ( 69 )
𝜌𝑔 2𝑔 𝐷 2𝑔
Tabla 43
Características Unidad
de alimentación con la disposición mostrada en Figura 72, teniendo una cabeza de succión positiva, se
garantiza una correcta alimentación en ambas bombas evitando la cavitación de las mismas, por tal
Luego de haber aprobado la selección de una bomba centrifuga para el flujo de agua (Véase
realiza verificando la necesidad que se requiere para flujo de agua según Tabla 43 y buscando la mejor
Finalmente, la bomba centrifuga SHIMGE modelo CPm130, se selecciona al cumplir con las
necesidades requeridas y ser adecuada para el transporte de agua a 10 ℃, teniendo como límite de
Figura 73
Luego de haber aprobado la selección de una bomba periférica para el flujo de aceite (Véase
realiza verificando la necesidad que se requiere para flujo de aceite según Tabla 43 y buscando la mejor
Finalmente, la bomba periférica EBARA modelo PRA 0.50 M permite el bombeo de aceite con
mercado que permite el suministro de altas presiones a cambio de bajos caudales, cumple con la
necesidad requerida para el bombeo de aceite, además, las presiones del fabricante con respecto a la
carga calculada poseen un amplio rango de trabajo seguro. Por lo tanto, se selecciona al cumplir con las
necesidades requeridas y ser adecuada para el transporte de aceite a 55 ℃, teniendo como límite de
trabajo una temperatura de hasta 80 ℃, además cumple con un comportamiento altamente resistente
Figura 74
Para realizar un análisis y elegir qué control es más adecuado para la temperatura en el lado frio
y lado caliente, se realiza el modelamiento del sistema para después acudir a las herramientas de
simulación y saber cómo se va a comportar la planta con las condiciones a las que fueron diseñadas.
vapor, donde la temperatura de sobrecalentamiento del refrigerante será la variable a controlar, la cual
variará con respecto a la temperatura ideal del refrigerante, siendo esta, la temperatura que permita la
Los componentes a analizarse dentro del proceso de refrigeración, en este caso evaporador y
condensador, serán los sistemas al cual se les aplicará un balance de energía dada por la ecuación 1.
𝑑𝑇𝑟
𝑚̇ℎ4 + 𝑄̇𝑟 − 𝑚̇ℎ1 = 𝑚𝐶𝑝
𝑑𝑡
𝑑𝑇𝑟
𝑚̇(ℎ4 − ℎ1 ) + 𝑄̇𝑟 = 𝜌𝐶𝑝 𝑉 ( 70 )
𝑑𝑡
180
Figura 75
1
𝑈𝐶 =
𝑑𝑜 d ln (𝑟𝑜 /𝑟𝑖 ) 1
+ o +
𝑑𝑖 ℎ𝑟 𝑘𝑐𝑢 ℎ𝑎𝑐
Para este análisis, el área y el coeficiente global de transferencia de calor serán necesarios. Por lo
refrigeración, en el evaporador:
𝑄̇𝑅 = 𝑈𝑒 𝐴𝑒 (𝑇𝐸𝑅 − 𝑇𝑟 )
Es necesario expresar las entalpias en función de la temperatura del refrigerante, por medio de
una regresión lineal como se muestra en la Figura 76 y Figura 77. Para realizar el ajuste de curvas se ha
tomado las propiedades termofísicas del refrigerante R404a de acuerdo con las temperaturas y
presiones de trabajo dentro del sistema diseñado, en este caso como la temperatura de condensación y
evaporación son de 55°C y 0°C, correspondientes a 2.5 Mpa y 0.6 Mpa, respectivamente. Se elige las
propiedades entre 0.6 Mpa y 2.5 Mpa de la tabla de Propiedades de líquido y vapor saturado del
Figura 76
300
0
-10 0 10 20 30 40 50 60
Temperatura [°C]
Figura 77
380
375
370 y = 0.2746x + 368.17
365 R² = 0.9155
360
-10 0 10 20 30 40 50 60
Temperatura [°C]
Tomando en cuenta que la entalpia en el punto 3 es igual a la entalpia del punto 4 a la entrada
ℎ3 = ℎ4 = ℎ4 𝑙𝑖𝑞.𝑠𝑎𝑡. + 𝑥(ℎ𝑓𝑔 ) ( 71 )
expresión:
( 73 )
ℎ4 𝑙𝑖𝑞.𝑠𝑎𝑡. (𝑇𝑟 ) = 197.4 + 1.5788 𝑇𝑟
del refrigerante.
𝐿 = 368.17 𝑘𝐽/𝑘𝑔
𝐺 = 0.2746 𝑘𝐽/𝑘𝑔𝐾
Donde las constantes L y G, describen una ecuación lineal, modelo que describe el valor de las
entalpias dentro de la curva de saturación (mezcla liquido – vapor). A su vez la entalpia ℎ4 𝑙𝑖𝑞 𝑠𝑎𝑡 , es una
ecuación que describe el comportamiento de entalpias del refrigerante en la curva de líquido saturado.
La ecuación 75, describirá la entalpia aproximada que tendrá el refrigerante en los puntos 3 y 4,
Sustituyendo:
𝑑𝑇𝑟
𝜌𝐶𝑝 𝑉 = 𝑈𝑒 𝐴𝑒 (𝑇𝐸𝑅 − 𝑇𝑟 ) + 𝑚̇ (ℎ4 𝑙𝑖𝑞 𝑠𝑎𝑡 (𝑇𝑟 ) + 𝑥(𝐿 − 𝐺𝑇𝑟 ) − ℎ1 (𝑇𝑟 ))
𝑑𝑡
Se tiene que:
𝑑𝑇𝑟
𝜌𝐶𝑝 𝑉 = 𝑈𝑒 𝐴𝑒 (𝑇𝐸𝑅 − 𝑇𝑟 ) + 𝑚̇(−(𝐿 − 𝐺𝑇𝑟 ) + 𝑥(𝐿 − 𝐺𝑇𝑟 ))
𝑑𝑡
Finalmente, obteniendo:
𝑑𝑇𝑟
𝜌𝐶𝑝 𝑉 = 𝑈𝑒 𝐴𝑒 (𝑇𝐸𝑅 − 𝑇𝑟 ) + 𝑚̇(𝑥 − 1)(𝐿 − 𝐺𝑇𝑟 ) ( 76 )
𝑑𝑡
𝑇𝑅 (𝜌, 𝐶𝑝 , 𝑉, 𝑈𝐸 , 𝐴𝐸 , 𝑚̇, 𝑥, 𝐿, 𝐺, 𝑡), por lo que es indispensable adimensionar la misma para evitar el uso de
todas ellas, y para lograrlo se dispone de las siguientes variables adimensionales (Flores, 2006).
𝑇𝐸𝑅 − 𝑇𝑟 ( 77 )
Θ𝑒 =
Δ𝑇𝑐
𝑡 ( 78 )
𝜏=
𝑡𝑐
Donde:
• 𝑇𝑐 temperatura característica
• 𝑡𝑐 tiempo característico
184
𝑡 = 𝑡𝑐 𝜏 ( 79 )
𝑇𝑟 = 𝑇𝐸𝑅 − Θ𝑒 Δ𝑇𝑐 ( 80 )
𝑑𝑇𝑟
𝜌𝐶𝑝 𝑉 = 𝑈𝑒 𝐴𝑒 (𝑇𝐸𝑅 − 𝑇𝑟 ) + 𝑚̇(𝑥 − 1)(𝐿 − 𝐺𝑇𝑟 )
𝑑𝑡
𝑑(𝑇𝐸𝑅 − Θ𝐸 Δ𝑇𝑐 )
𝜌𝐶𝑝 𝑉 = 𝑈𝑒 𝐴𝑒 (Θ𝑒 Δ𝑇𝑐 ) + 𝑚̇(𝑥 − 1)(𝐿 − 𝐺(𝑇𝐸𝑅 − Θ𝑒 Δ𝑇𝑐 ))
𝑑(𝑡𝑐 𝜏)
𝜌𝐶𝑝 𝑉
𝑡𝑐 = ( 82 )
𝑈𝑒 𝐴𝑒
𝑚̇𝐿
Δ𝑇𝑐 = ( 83 )
𝑈𝑒 𝐴𝑒
𝑑Θ𝑒 𝐺 Δ𝑇𝑐
− = Θ𝑒 + (𝑥 − 1) [1 − 𝑇𝐸𝑅 (1 − Θ )] ( 84 )
𝑑𝜏 𝐿 𝑇𝐸𝑅 𝑒
185
Δ𝑇𝑐 𝑚̇𝐿
𝛾= = ( 85 )
𝑇𝐸𝑅 𝑈𝑒 𝐴𝑒 𝑇𝐸𝑅
𝐺
𝛼= 𝑇 ( 86 )
𝐿 𝐸𝑅
𝑑Θ𝑒
− = Θ𝑒 + (𝑥 − 1)[1 − 𝛼(1 − 𝛾Θ𝑒 )] ( 87 )
𝑑𝜏
Figura 78
se tiene:
𝑑𝐸𝑠𝑖𝑠𝑡𝑒𝑚𝑎
𝐸̇𝑒𝑛𝑡𝑟𝑎 + 𝐸̇𝑔𝑒𝑛𝑒𝑟𝑎𝑑𝑎 − 𝐸̇𝑠𝑎𝑙𝑒 =
𝑑𝑡
𝑑𝑇𝐶
𝑄̇𝑅𝑒 + 𝑚̇𝑟 ∗ ℎ2 + (−𝑄̇𝑐 ) − 𝑚̇𝑟 ∙ ℎ3 = 𝑚𝑟 ∗ 𝐶𝑝 ∗ ( 88 )
𝑑𝑡
186
El calor aportado por la resistencia eléctrica se define por la ley de Ohm, en una relación entre el
voltaje y la resistencia, como:
𝑉𝑜2
𝑄̇𝑅𝑒 = 𝑃𝑒 = ( 89 )
𝑅
Donde:
𝑅𝑒 : 𝑅𝑒𝑠𝑖𝑠𝑡𝑒𝑛𝑐𝑖𝑎 𝑒𝑙é𝑐𝑡𝑟𝑖𝑐𝑎
El calor cedido por el serpentín del condensador hacia el entorno exterior sumergido (Aceite),
es:
𝑄𝑐̇ = 𝑈𝑐 ∗ 𝐴𝑐 ∗ (𝑇𝑐 − 𝑇∞ ) ( 90 )
1
𝑈𝐶 =
𝑑𝑜 d ln (𝑟𝑜 /𝑟𝑖 ) 1
+ o +
𝑑𝑖 ℎ𝑟 𝑘𝑐𝑢 ℎ𝑎𝑐
𝑑𝑇𝐶
̇ + 𝑚̇ 𝑟 (ℎ2 − ℎ3 ) − 𝑈𝑐 ∗ 𝐴𝑐 ∗ (𝑇𝑐 − 𝑇∞ ) = 𝑚𝑟 ∗ 𝐶𝑝 ∗
𝑄𝑅𝑒 ( 91 )
𝑑𝑡
Para analizar la energía del refrigerante en la entrada al condensador (punto 2), se supone que
ℎ2 = ℎ1 + 𝐶𝑝 (𝑇2 − 𝑇1 ) ( 92 )
como: ℎ1 = ℎ1 (𝑇𝑟 ).
187
𝑘−1
𝑇2 𝑃2 𝑘
( )=( )
𝑇1 𝑃1
𝑃2
𝑟𝑝 = ( 93 )
𝑃1
Por lo tanto, sustituyendo la ecuación 93, en ecuación 17, que define el proceso isentrópico, se
tiene:
𝑘−1
𝑇2 = 𝑇1 (𝑟𝑝 ) 𝑘 ( 94 )
Donde
𝑑𝑇𝐶
̇ + 𝑚̇ 𝑟 (ℎ1 (𝑇𝑟 ) + 𝐶𝑝 (𝑇2 − 𝑇1 ) − ℎ4
𝑄𝑅𝑒 (𝑇 )
𝑙𝑖𝑞.𝑠𝑎𝑡. 𝑟 − 𝑥(𝐿 − 𝐺 𝑇𝑟
)) − 𝑈𝐶 𝐴𝑐 (𝑇𝑐 − 𝑇∞ ) = 𝜌 𝐶𝑝 𝑉 ( 95 )
𝑑𝑡
La entalpia de vaporización entre los puntos 4 y 1 del sistema, es la diferencia entre la entalpia
ℎ𝑓𝑔 = ℎ1 − ℎ4𝑙𝑖𝑞.𝑠𝑎𝑡. ( 96 )
188
𝑘−1 𝑑𝑇𝐶
̇ + 𝑚̇ 𝑟 (ℎ𝑓𝑔 (𝑇𝑟 ) + 𝐶𝑝 (𝑇𝑟 (𝑟𝑝 )
𝑄𝑅𝑒 𝑘 − 𝑇𝑟 ) − 𝑥(𝐿 − 𝐺 𝑇𝑟 )) − 𝑈𝐶 𝐴𝑐 (𝑇𝑐 − 𝑇∞ ) = 𝜌 𝐶𝑝 𝑉
𝑑𝑡
𝑘−1 𝑑𝑇𝐶
̇ + 𝑚̇ 𝑟 ((𝐿 − 𝐺 𝑇𝑟 ) + 𝐶𝑝 (𝑇𝑟 (𝑟𝑝 )
𝑄𝑅𝑒 𝑘 − 𝑇𝑟 ) − 𝑥(𝐿 − 𝐺 𝑇𝑟 )) − 𝑈𝐶 𝐴𝑐 (𝑇𝑐 − 𝑇∞ ) = 𝜌 𝐶𝑝 𝑉
𝑑𝑡
𝑘−1 𝑑𝑇𝐶
̇ + 𝑚̇ 𝑟 (𝐶𝑝 (𝑇𝑟 (𝑟𝑝 )
𝑄𝑅𝑒 𝑘 − 𝑇𝑟 ) + (𝐿 − 𝐺 𝑇𝑟 )(1 − 𝑥)) − 𝑈𝐶 𝐴𝑐 (𝑇𝑐 − 𝑇∞ ) = 𝜌 𝐶𝑝 𝑉
𝑑𝑡
𝑘−1 𝐺 𝑑𝑇𝐶
̇ + 𝑚̇ 𝑟 (𝐶𝑝 (𝑇𝑟 (𝑟𝑝 )
𝑄𝑅𝑒 𝑘 − 𝑇𝑟 ) + (1 − 𝑇 ) (1 − 𝑥)𝐿) − 𝑈𝐶 𝐴𝑐 (𝑇𝑐 − 𝑇∞ ) = 𝜌 𝐶𝑝 𝑉 ( 97 )
𝐿 𝑟 𝑑𝑡
Igual que en el método del evaporador, para adimensionalizar la ecuación del condensador se
𝑇𝐶 − 𝑇∞ ( 98 )
𝜃𝐶 =
∆𝑇𝐶
𝑡 ( 99 )
𝜏=
𝑡𝑐
Donde:
𝑇𝐶 : 𝑇𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑐𝑎𝑟𝑎𝑐𝑡𝑒𝑟í𝑡𝑖𝑐𝑎
𝑡𝐶 : 𝑇𝑖𝑒𝑚𝑝𝑜 𝑐𝑎𝑟𝑎𝑐𝑡𝑒𝑟í𝑠𝑡𝑖𝑐𝑜
( 100 )
𝑇𝐶 = 𝑇∞ + 𝜃𝐶 ∆𝑇𝐶
𝑡 = 𝑡𝑐 𝜏 ( 101 )
189
Tomando en cuenta las ecuaciones 94, 100 y 101 al reemplazar en ecuación 97, se tiene:
𝑘−1 𝐺 𝑑𝑇𝐶
̇ + 𝑚̇ 𝑟 (𝐶𝑝 𝑇𝑟 ((𝑟𝑝 )
𝑄𝑅𝑒 𝑘 − 1) + (1 − 𝑇𝐸𝑅 − 𝜃𝑒 ∆𝑇𝐶 ) (1 − 𝑥)𝐿) − 𝑈𝐶 𝐴𝑐 (𝜃𝐶 ∆𝑇𝐶 ) = 𝜌 𝐶𝑝 𝑉
𝐿 𝑑𝑡
𝐶𝑝 𝑇𝑅 𝑘−1 𝐺 𝜃𝑒
̇ + ∆𝑇𝐶 (𝑈𝑒 𝐴𝑒 ) (
𝑄𝑅𝑒 ((𝑟𝑝 ) 𝑘 − 1) + (1 − 𝑥) (1 − 𝑇𝐸𝑅 (1 − ∆𝑇 ))) − 𝑈𝐶 𝐴𝑐 (𝜃𝐶 ∆𝑇𝐶 )
𝐿 𝐿 𝑇𝐸𝑅 𝐶
( 102 )
1 𝑑(𝑇∞ ) ∆𝑇𝐶 𝑑(𝜃𝐶 )
= 𝑡𝑐 (𝑈𝑒 𝐴𝑒 ) ( ( + ))
𝑡𝑐 𝑑(𝜏) 𝑑(𝜏)
𝑑(𝑇∞ )
=0
𝑑(𝜏)
𝐶𝑝 𝑇𝑟 𝑘−1 𝐺 𝜃𝑒
̇ + ∆𝑇𝐶 (𝑈𝑒 𝐴𝑒 ) (
𝑄𝑅𝑒 ((𝑟𝑝 ) 𝑘 − 1) + (1 − 𝑥) (1 − 𝑇𝐸𝑅 (1 − ∆𝑇 ))) − 𝑈𝐶 𝐴𝑐 (𝜃𝐶 ∆𝑇𝐶 )
𝐿 𝐿 𝑇𝐸𝑅 𝐶
( 103 )
∆𝑇𝐶 𝑑(𝜃𝐶 )
= (𝑈𝑒 𝐴𝑒 )
𝑑(𝜏)
1
Se multiplica a toda la ecuación 103 por ∆𝑇𝐶 (𝑈𝑒 𝐴𝑒 )
Ahora se define a la conductancia como la razón entre los coeficientes globales y áreas de
𝑈𝐶 𝐴𝑐
𝛽= ( 105 )
𝑈𝑒 𝐴𝑒
Sustituyendo las ecuaciones del evaporador 80 y 85, donde interviene la temperatura del
𝜃𝑒 ∆𝑇𝐶
𝑇𝑟 = 𝑇𝐸𝑅 (1 − )
𝑇𝐸𝑅
𝑇𝑟 = 𝑇𝐸𝑅 (1 − 𝛾 𝜃𝑒 ) ( 107 )
Se crea una nueva constante lambda, , que relaciona el calor especifico con la temperatura del
𝐶𝑝
= 𝑇 ( 108 )
𝐿 𝐸𝑅
̇
𝑄𝑅𝑒 𝑘−1 𝑑(𝜃𝐶 )
+ 𝑇𝑟 ((𝑟𝑝 ) 𝑘 − 1) + (1 − 𝑥) (1−∝ (1 − 𝛾 𝜃𝑒 )) − 𝛽𝜃𝐶 =
∆𝑇𝐶 (𝑈𝑒 𝐴𝑒 ) 𝑇𝐸𝑅 𝑑(𝜏)
𝑃𝑒 𝑘−1 𝑑(𝜃𝐶 )
+ (1 − 𝛾 𝜃𝑒 ) ((𝑟𝑝 ) 𝑘 − 1) + (1 − 𝑥) (1−∝ (1 − 𝛾 𝜃𝑒 )) − 𝛽𝜃𝐶 = ( 109 )
∆𝑇𝐶 (𝑈𝑒 𝐴𝑒 ) 𝑑(𝜏)
El gasto másico en el compresor viene dado por la siguiente relación (Flores, 2006, p. 25):
𝑘−1 ( 110 )
𝑚̇ 𝑟 = 𝑁𝑐𝑜𝑚𝑝 𝑉𝑃 𝜌 (1 + 𝐶𝑟 − 𝐶𝑟 𝑟𝑝 𝑘 )
191
2006, p. 25):
1
1 𝑃2 2
( 111 )
𝑚̇ 𝑟 = 𝐾 𝜌2 (𝑃1 ( − 1))
𝑃1
Donde:
igualan los gastos másicos del compresor (ecuación 110) y válvula de expansión (ecuación 117)
1
1
2
𝑁 𝑉𝑃 𝜌 = 𝐾 𝜌2 (𝑃1 (𝑟𝑝 − 1))
𝑁 𝑉𝑃 2 𝜌
( ) = (𝑟𝑝 − 1) ( 112 )
𝐾 𝑃1
𝜌 1
= ( 113 )
𝑃1 𝑅 𝑇𝑟
𝜌 1
= ( 114 )
𝑃1 𝑅 𝑇𝐸𝑅 (1 − 𝛾 𝜃𝑒 )
En donde para relacionar a la temperatura del refrigerante con las expresiones del compresor y
válvula de expansión, se reemplazan la ecuación 114 en 112:
𝑁 𝑉𝑃 2 1 ( 115 )
( ) = (𝑟𝑝 − 1)
𝐾 𝑅 𝑇𝐸𝑅 (1 − 𝛾 𝜃𝑒 )
192
Se reescribe la ecuación 115, creando una nueva constante épsilon (𝜀) que depende de las
𝑁 𝑉𝑃 2 1
𝜀=( ) ( 116 )
𝐾 𝑅 𝑇𝐸𝑅
𝜀
𝑟𝑝 = +1 ( 117 )
(1 − 𝛾 𝜃𝑒 )
Finalmente, se reemplaza la relación 117 en ecuación 109 y se tiene el modelo matemático para
𝑘−1
𝑃𝑒 𝜀 𝑘 𝑑(𝜃𝐶 )
+ (1 − 𝛾 𝜃𝑒 ) (( + 1) − 1) + (1 − 𝑥)(1−∝ (1 − 𝛾 𝜃𝑒 )) − 𝛽𝜃𝐶 = ( 118 )
∆𝑇𝐶 (𝑈𝑒 𝐴𝑒 ) (1 − 𝛾 𝜃𝑒 ) 𝑑(𝜏)
por lo que se trabaja con la ecuación 87 para aplicar el control. En Figura 79 se presentan los
Figura 79
En Figura 80 se presenta el esquema del evaporador con el modelo que gobierna al equipo.
Figura 80
𝑑Θ𝑒
− = Θ𝑒 + (𝑥 − 1)[1 − 𝛼(1 − 𝛾Θ𝑒 )]
𝑑𝜏
𝑇𝐸𝑅 − 𝑇𝑟
Θ𝑒 =
Δ𝑇𝑐
194
𝑚̇ 𝑟 𝐿
Δ𝑇𝑐 =
𝑈𝑒 𝐴𝑒
𝑚̇ 𝑟 𝐿
𝛾=
𝑈𝑒 𝐴𝑒 𝑇𝐸𝑅
𝐺
𝛼= 𝑇
𝐿 𝐸𝑅
Modelo Numérico
Aplicando la teoría clásica del control para analizar un ciclo ideal de refrigeración por
compresión de vapor ante una entrada escalón unitario y al realizar la transformada de Laplace en la
1 1
−(𝑠𝜃𝐸 (𝑠) − 𝜃𝐸 (0)) = (1 − 𝑥) [ − 𝛼 ( − 𝛾 𝜃𝐸 (𝑠))] − 𝜃𝐸 (𝑠)
𝑠 𝑠
Tomando en cuenta la condición inicial que 𝜃𝐸 (𝑠) = 0, entonces nos queda de la siguiente
manera:
(𝑥 − 1)(𝛼 − 1)
𝜃𝐸 (𝑠) =
𝑠2 + (𝛼𝛾(𝑥 − 1) + 1)𝑠
Se toma en cuenta, para estudio que es un ciclo ideal, con las condiciones en
Tabla 44
Estado 1 Estado 2
𝑻𝟏 = 𝟕 ℃ 𝑇2 = 68 ℃
Estado 3 Estado 4
𝑻𝟑 = 𝟒𝟒 ℃ 𝑇4 = 0 ℃
𝐶𝑝 = 2.095 𝑘𝐽/𝑘𝑔𝐾
𝑘 = 2.032
𝑊 𝑊
107.78 y de 440.03 del evaporador, con un área de 2.082 𝑚2 y de 0.665 𝑚2 para el
𝑚2 𝐾 𝑚2 𝐾
𝐿 = 368.17 𝐾𝐽/𝑘𝑔
𝐺 = 0.2746 𝐾𝐽/𝑘𝑔
196
∆𝑇𝐶 𝑚̇𝐿
𝛾= = = 2.7342 × 10−4
𝑇𝐸𝑅 𝑈𝑒 𝐴𝑒 𝑇𝐸𝑅
𝐺
𝛼= 𝑇 = 0.211076
𝐿 𝐸𝑅
El sistema de refrigeración posee dos componentes indispensables que permiten que el sistema
se comporte con las condiciones deseadas de trabajo; siendo estas, la válvula de expansión y compresor.
Figura 81
Para poder hallar la respuesta del sistema, con una entrada de escalón unitario, se toma en
comparan los tres tipos de control con respecto a la respuesta independiente del sistema cuando no se
tiene control.
Tabla 45
Controlador 𝑲𝒑 𝑲𝒊 𝒌𝒅
P 7.03056 − −
PI 6.72030 0.19080 −
Figura 82
De la Figura 82 se observa que los controladores P, PI y PID, tardan en estabilizarse con respecto
al comportamiento libre del sistema, por lo que, se opta por reajustar las constantes proporcionales,
integrales y derivativas.
198
En Figura 83 se realiza una segunda comparación entre los tres controles con la finalidad de
Tabla 46
Controlador 𝑲𝒑 𝑲𝒊 𝒌𝒅
P 19.25786 − −
PI 10.97473 0.41076 −
Figura 83
Comparando el control proporcional integral derivativo con los controles P y PI de la Figura 83,
se observa que el control PID se estabiliza en menor tiempo que los otros dos controladores, por esta
razón se selecciona el control PID para aplicar en el sistema de refrigeración, esto permite que el
199
en el evaporador disminuya de forma paulatina quitando el calor con la ayuda del refrigerante.
Por otro lado, como se mencionó en el modelado dinámico del condensador, este depende del
evaporador, por lo que el control independiente del lado caliente no es posible mediante el uso de este
sistema. De todas maneras, se quiere implementar un control PID con la ayuda de una niquelina que
será controlada por un control de fase, enviando pulsos de voltaje para calentar proporcionalmente al
condensador de la Tabla 72, cuyas constantes del control fueron 204, 7.2, 1.04 para las constantes
Figura 84
El Banco de Pruebas se encuentra conformado por cinco (5) sistemas o elementos que se
desarrollan y operan de manera integral, los cuales son: Estructura Metálica, Tanques de
200
Figura 85
fundamentales, tal como se muestra en Tabla 47. A su vez, las características de cada una de dichos
Tabla 47
1 Estructura Metálica
3 Sistema de Refrigeración
diagrama unifilar (Apéndice 6). En Tabla 48 se enlista el peso de los componentes que serán colocados
sobre la estructura.
Tabla 48
Peso de componentes
362.37 kg, correspondiente a la totalidad del banco de pruebas. Construida con perfil cuadrado de 1’’
con 2 mm de espesor, perfil cuadrado de 2” con 3 mm de espesor y una base de una unidad
condensadora. Los procesos seguidos para construir la estructura fueron: corte a 45° y 90°, soldadura
MIG, desbaste y afinado de uniones con amoladora, taladrado y finalmente, ensamble con sus
86.
Figura 86, las cuales representan el peso de los componentes con un factor que aumenta la carga en un
30%.
Figura 86
Figura 87
Tanques de Almacenamiento
condensadora y evaporadora. Tanto condensador como evaporador cuentan con dos tanques
Tanques Internos. El tanque interior servirá de almacenamiento para los fluidos de trabajo,
aceite y agua. Fabricados en acero inoxidable 430, debido a que permanecerán en contacto directo,
permanentemente, con el fluido de trabajo evitando la oxidación temprana del material y alargando su
vida útil. Para su proceso de fabricación se realizó una deconstrucción de la geometría cilíndrica para
obtener las dimensiones para cortar la plancha de acero inoxidable con oxicorte en maquinaria CNC,
obteniendo así dos geometrías después del corte, un rectángulo y un círculo. El rectángulo se utilizó
para el cuerpo del tanque, después del corte se procede a rolarlo próximo al diámetro requerido, luego
se procede a realizar la soldadura con proceso TIG y varilla de aporte 308L, para luego volver a rolarlo
204
con el fin de obtener el diámetro requerido. Luego ambos elementos se sueldan bajo el mismo proceso
del cuerpo y se procede a realizar los agujeros para los accesorios mediante un taladro de banco con
mesa adaptada para realizar agujeros en superficies curvas. Una vez realizados los agujeros se procede a
realizar la soldadura de los acoples y se colocaron tapones para realizar una prueba de fugas llenando el
tanque a tope con agua. Una vez culminado el proceso y con los tanques limpios se realizó un decapado
en los cordones de soldadura con el fin de evitar oxidación de estos con un Gel Decapante.
Figura 88
Tanques internos
Tanques Externos. El tanque exterior es el recipiente que abarcará al tanque interno. Fabricados
en acero A36. Para su proceso de fabricación se realizó una deconstrucción de la geometría cilíndrica
para obtener las dimensiones para cortar la plancha de acero A36 con oxicorte en maquinaria CNC,
obteniendo así dos geometrías después del corte, un rectángulo y un círculo. El rectángulo se utilizó
para el cuerpo del tanque, después del corte se procede a rolarlo próximo al diámetro requerido, luego
se procede a realizar la soldadura con proceso SMAW con electrodo revestido 6010, para luego volver a
rolarlo con el fin de obtener el diámetro requerido. Luego ambos elementos se sueldan bajo el mismo
205
proceso del cuerpo y se procede a realizar los agujeros para los accesorios mediante un taladro de
banco con mesa adaptada para realizar agujeros en superficies curvas. Para culminar el proceso se
realizaron los agujeros por donde se conectarían los tanques internos con el Sistema de Alimentación y
Figura 89
Tanques externos
Pintura y Aislamiento Térmico. Los tanques externos fueron limpiados y pintados con pintura
electrostática que brindan un gran recubrimiento superficial, durabilidad y secado rápido en el metal, a
su vez, el espacio anular generado entre ambos tanques (externo e interno) servirá como acopio de
Figura 90
Tanques de almacenamiento
Sistema de Refrigeración
El sistema de refrigeración consta de varios elementos que son necesarios tanto para el
detalla los elementos necesarios que se instalaron en el sistema y se los identifica en el Banco de
Figura 91
Sistema de refrigeración
Nota: 1. Llave de cierre manual Evaporador, 2. Visor de líquido, 3. Llave de cierre manual, 4. Válvula
antirretorno, 10. Llave de cierre manual Condensador entrada, 11. Llave de cierre manual Condensador
salida, 12. Válvula Rotolock del recibidor de líquido, 13. Recibidor de líquido, 14. Filtro deshidratador,
Compresor. El compresor se comporta como una bomba mecánica que permite mover el
refrigerante dentro del circuito de refrigeración, esta unidad aspira el vapor a la salida del evaporador y
En Figura 92, se muestra el compresor CAJ4519Z de la marca Tecumseh, junto con el kit de
arranque necesario para dar el primer impulso de trabajo, además de la válvula rotolock, necesaria para
realizar el vacío y carga de refrigerante, se ubican en la parte baja entren el evaporador y condensador
Figura 92
Manómetro de Presión de Alta y Baja. Con el objetivo de identificar las presiones a las que se
encuentra el refrigerante antes de ingresar al condensador y después de salir del evaporador, se coloca
dos manómetros, uno en la salida del condensador para medir la presión de alta y el otro antes de
ingresar al compresor para medir la presión de baja en el sistema. En Figura 93 se muestra los
manómetros de alta y baja que tienen un rango de 0 a 800 psi y de 0 a 500 psi, respectivamente.
Figura 93
Nota. 1. Manómetro de baja presión, 2. Entrada al compresor, 3. Manómetro de alta presión, 4. Salida del
condensador.
209
Presostatos de Alta y Baja. Este elemento mide una presión diferencial, se coloca en el sistema
para proteger al compresor, comparando las presiones cuando sobrepasa el límite de presión a la salida
del acumulador de líquido o cuando baja más de lo configurado en la línea de baja presión, abriendo o
cerrando un interruptor. En Figura 94 se muestra el presostato de alta y baja, de la marca Danfoss que
Figura 94
Nota. 1. Presostato de baja presión, 2. Línea de baja presión, 3. Presostato de alta presión, 4. Salida del
recibidor de líquido.
Recibidor de Líquido. El recibidor de líquido se coloca a la salida del condensador con el fin de
la válvula de expansión. Otra ventaja de colocar el acumulador es que permite almacenar el refrigerante
cuando sea necesario realizar mantenimiento del sistema. En Figura 95 se muestra el contenedor de
Figura 95
Recibidor de líquido
Nota. 1. Salida del acumulador, 2. Entrada del acumulador, 3. Acumulador de líquido, 4. Válvula Rotolock
Filtro Deshidratador. Debido a que por las tuberías del sistema de refrigeración no circula
solamente refrigerante sino también una pequeña cantidad de aceite que permite la lubricación del
compresor es necesario que la mezcla este totalmente seca cuando llega a la válvula de expansión y por
Visor de líquido. Se coloca un visor de líquido de diámetro 3/8” junto al filtro deshidratador
como se muestra en Figura 96, ya que permite observar en que condición circula el refrigerante antes de
Figura 96
Válvula Solenoide. Se coloca la válvula solenoide para apagar el compresor por el lado de baja
presión, la cual al llegar a una presión determinada se cierra completamente el paso de refrigerante
Figura 97
Válvula solenoide
Válvula de Cierre Manual. En el sistema se han colocado cuatro válvulas de cierre manual como
las que se muestran En Figura 98, con el fin de permitir el mantenimiento de una forma segura y
evitando que existan fugas de refrigerante, se colocan a la entrada y salida tanto del evaporador como
del condensador, esto permite cerrar el paso de refrigerante por el sistema cuando sea necesario
realizar una soldadura en caso de fugas y hacer vacío en esa línea sin desperdiciar refrigerante.
Figura 98
Nota. 1. Entrada del condensador, 2. Salida del condensador, 3. Entrada del evaporador, 4. Salida del
evaporador.
212
Válvula de Expansión Electrónica. Como se mencionó en el capítulo III, se selecciona una válvula
de expansión electrónica para refrigerante R404A, este se coloca antes del evaporador con el fin de que
se regule el paso de refrigerante, siendo el causante de pasar de una presión alta a una presión baja, la
Figura 99
para poder realizar la transferencia de calor entre el refrigerante y el medio circundante. Los serpentines
son sujetados a una estructura de canales ranurados y soportados en una base de perfiles en L de acero
A36 con pintura electrostática. En Figura 100 y Figura 100 se puede observar los serpentines para
Figura 100
Figura 101
Sistema de Alimentación
En el sistema de bombeo es importante mencionar todo el sistema está dividido en dos circuitos
separados, uno para agua y otro para aceite, con una bomba respectiva para cada fluido y anclados a la
estructura metálica.
Flujómetros para Sistemas de Alimentación. Para poder medir el flujo que pasa a través de
características.
Tabla 49
Brooks Instrument
Agua 0.1 − 0.8 1305D08A1A1A
Division
Brooks Instrument
Aceite 1 − 15 1305D10A3A1A
Division
Nota. Adaptado de Variable Area Flow Meters (Rotameters), por Brooks Instrument, 2013.
214
Figura 102
Nota. Flujómetro para sistema de alimentación de agua (izquierda), flujómetro para sistema de
Sistema de Alimentación para Agua. La bomba utilizada para impulsar el fluido agua es de la
Figura 103
Esta conecta su succión mediante un neplo de 1” a un acople universal, con el fin de poder
desmantelar la bomba, a continuación, está conectado a un acople tipo T, que se acopla por el extremo
de salida a una llave de paso para poder liberar el líquido del tanque, con una reducción de 1” a 1/2", y
por el otro extremo se conecta un neplo de 1” junto a un codo para lograr el ingreso al tanque mediante
otro neplo de 1”, que a su vez está conectado a una llave de paso tipo esférica para controlar el ingreso
Figura 104
Nota. 1. Neplos, 2. Llave de paso esférica, 3. Codo, 4. Reducción, 5. Llave de salida, 6. Conexión T, 7.
seguido de una conexión universal junto con un segmento de tubería que llega a la altura del neplo de
entrada del tanque para poder colocar una conexión Bypass. Se acopla mediante un acople tipo T que
de forma vertical continua con dos codos unidos mediante un neplo de 1” para poder reducir la presión
del fluido hasta llegar al flujómetro junto a una llave de paso para poder controlar el flujo de líquido que
216
ingresa a través un neplo metálico de 1”, a la salida del flujómetro se encuentra el acople rápido para
Para el Bypass, en el otro extremo del acople tipo T, se conecta también una llave de paso con
un neplo de 1” conectado hacia un codo, que se enlaza con tubería para ingresar al tanque, dentro del
tanque con la misma conexión se emplea un codo para evitar el desfogue del líquido dentro del tanque.
Figura 105
Conexión T, 8. Tubería 3/4", 9. Unión, 10. Neplos metálicos, 11. Acople Rápido
217
Sistema de Alimentación para Aceite. Se utilizar la bomba marca Ebara tipo PRA 046 M que se
Figura 106
Siendo esta, conectada de forma similar al sistema de alimentación del agua desde la salida del
Figura 107
A continuación, para la descarga de la bomba para aceite y el Bypass se procede con la misma
Figura 108
Nota. 1. Neplo metálico, 2. Unión, 3. Conexión T, 4. Codo, 5. Neplo 3/4", 6. Conexión Universal, 7.
Reducción 1” a 3/4”, 8. Neplo 1”, 9. Tubería 3/4”, 10. Llave de paso, 11. Flujómetro.
219
impresión de la placa (PCB) y el circuito de control físico; y el Software que incluye la programación en
su aplicación es sumamente fácil, posee algunos prototipos electrónicos con código abierto en los cuales
Proteus. Es un Software de CAD electrónico con un sistema de modelado virtual que combina la
simulación del circuito, sus componentes y los distintos microprocesadores para realizar la simulación
del microcontrolador que además resulta el esquemático del circuito y layout del PCB (Viveros, 2009).
Nextion - Editor HMI. Es un software que permite el desenvolvimiento del GUI de interfaz
hombre-máquina o como sus siglas en inglés HMI, logra funciones para animaciones mediante
diapositivas táctiles, botones, textos, gráficos, etc. Para su codificación de interacción de componentes,
Se utilizó el software Proteus, que ayuda a realizar las conexiones necesarias y el ruteo de estas,
para poder implementar todos los componentes del sistema de control, sus conexiones y
220
esquematización, con el fin de desarrollar un circuito de control capaz de satisfacer las necesidades del
banco de pruebas. Se presenta el diseño esquemático del circuito de control en la Figura 109.
Figura 109
Una vez esquematizado el circuito y realizado un modelo con el ruteo para la conexión del
circuito de control, mostrado en Figura 110, se realiza la impresión de placas PCB en fibra de vidrio con
recubrimiento de cobre.
Una vez terminado el diseño se realiza la impresión de placas en la cual el modelo del circuito
diseñado, mostrado en Figura 110, representa el circuito a ser impreso en la placa hecha de material
Figura 110
transfer y se sumerge en cloruro férrico con leves movimientos para remover las partículas de cobre
sobrantes en el circuito impreso. Una vez finalizado el proceso obtuvimos el resultado ilustrado en la
Figura 111.
222
Figura 111
El siguiente paso es realizar las perforaciones necesarias en el circuito para continuar con el
Circuito de Control
circuito de control de la Figura 109, tomando en cuenta las posiciones indicadas en el diseño de la placa
Figura 112
15k, 8. Resistencia 33k, 9. Bornera para disipador de calor, 10. Bornera para niquelina, 11. Fuente 5V,
12. 2N3904, 13. Relé, 14. Resistencia de Relé, 15. Led indicador, 16. Resistencia, 17. PC817, 18. Bornera
para actuadores, 19. Salida HMI, 20. Led Testigo, 21. Bornera para termocuplas, 22. MAX31855JASA, 23.
Microcontrolador ESP32
Para el CPU del circuito de control se seleccionó el SoC tipo microcontrolador ESP32 – WROOM
– 32U, el cual es un sistema integrado en un chip con múltiples funciones, basado en el chip ESP32, con
el cual podemos almacenar y desarrollar la programación requerida para el funcionamiento del sistema
224
realizar la interpretación de los datos obtenidos en la interfaz HMI, exportar los datos obtenidos
Figura 113
Tabla 50
Denominación Especificación
Core ESP32-D0WD
Temperatura de operación
-40 °C ~ 85 °C
recomendad
Nota: Transmisor receptor asíncrono universal (UART), Estándar de comunicaciones para transferir
información entre circuitos integrados (SPI: Serial Peripheral Interface), Seguridad digital de entrada y
salida (SDIO), Circuito inter-integrado, comunicación entre diferentes partes de circuito (I2C), Pulse
Width Modulation (PWM), Interconexión de circuitos de audio digital (I2S), Módulo receptor infrarrojo
226
(IR), General Purpose Imput Output, sistema de entrada y salida de propósito general (GPIO), Conversor
subsistemas, tal como se detalla en Figura 109, de manera integral y directa en la placa PCB, ilustrada
sus conexiones en Figura 110 con los cuales se realiza el sistema de control.
Los componentes que se utilizaron en el sistema de control son los descritos en Tabla 51
Tabla 51
Microprocesador ESP32 –
1 1 CPU
WROOM – 32U
A continuación, se detallan las conexiones del CPU con los distintos componentes.
Tabla 52
Figura 114
Figura 115
Tabla 53
IO32 - 8 R4 (U2)
IO33 - 9 R3 (U14)
IO26 - 10 R2 (U15)
IO27 - 11 R1 (U16)
229
Figura 116
Tabla 54
Figura 117
Tabla 55
IO13 - 15 ZC – 5 (U3)
Figura 118
Tabla 56
SCK - 5 (U6)
SCK - 5 (U7)
IO19 - 31
SCK - 5 (U8)
SCK - 5 (U9)
MISO - 7 (U6)
MISO - 7 (U7)
IO18 - 30
MISO - 7 (U8)
MISO - 7 (U9)
232
Figura 119
Tabla 57
IO16 - 27 3 (HMI)
IO17 - 28 2 (HMI)
233
Figura 120
Tabla 58
Componente Pin
8 (U11:A)
8 (U12:A)
MCP6002 8 (U12:B)
8 (U13:A)
8 (U13:B)
1 (U2)
1 (U14)
1 (U15)
1 (U16)
PC817
4 (U2)
4 (U14)
4 (U15)
4 (U16)
D4
D6
LED-RED
D8
D10
D3
D5
1N4007
D7
D9
TBLOCK - I2 1 (J1)
IRM - 20 - 5 4 (U10)
235
Tabla 59
Componente Pin
GND1 - 14 (U1)
GND3 - 38 (U1)
2 (U11:A)
4 (U11:A)
4 (U12:A)
MCP6002
4 (U12:B)
4 (U13:A)
4 (U13:B)
2 (U2)
2 (U14)
PC817
2 (U15)
2 (U16)
1 (Q1)
1 (Q2)
2N3904
1 (Q3)
1 (Q4)
TBLOCK - I2 2 (J3)
IRM - 20 - 5 3 (U10)
H11AA1 4 (U3)
236
Los sistemas eléctricos de potencia son aquellos que distribuyen la energía eléctrica para el
Figura 121
Nota: 1. HMI, 2. Cerradura, 3. Luz piloto led, 4. Paro de emergencia, 5. Selector metálico (2 posiciones),
Tabla 60
Elemento Especificaciones
Figura 122
Figura 123
Figura 124
Tabla 61
Característica Especificación
Apreciación 0.1 ℃
Dimensiones 71 x 28 x 71 mm
indeseadas como altibajos de tensión o corriente excesiva, protegiendo el sistema eléctrico en el que se
lo aplico (La bodega eléctrica, 2020). Para este proyecto se utilizaron 2 breakers en total, 1 breaker que
alimenta a todo el sistema eléctrico de fuerza y 1 breaker para el controlador de temperatura. Las
Figura 125
Nota: 1. Breaker para sistema eléctrico de fuerza y 2. Breaker para controlador de temperatura.
240
Tabla 62
Sistema eléctrico de
CHINT Nxb 3p-63a
fuerza
Controlador de
CHINT Nxb/Ebc 1p-6a
temperatura
Nxc-25a 3p
Bomba de Agua CHINT
120/220v
Nxc-25a 3p
Bomba de Aceite CHINT
120/220v
corriente eléctrica recorre la bobina del contactor produciendo el cierre de los contactos y permitiendo
la conexión eléctrica (Eminca, 2022). Para este proyecto se utilizaron 4 contactores en total, 1 contactor
calentamiento y 2 contactores para las bombas de agua y aceite. Las especificaciones de cada contactor
Figura 126
Nota: 1. Contactor para compresor del sistema de refrigeración, 2. Contactor para resistencia eléctrica
Tabla 63
Nxc-32a 3p
Compresor CHINT
120/220v
Nxc-32a 3p
Resistencia Eléctrica CHINT
120/220v
Nxc-25a 3p
Bomba de Agua CHINT
120/220v
Nxc-25a 3p
Bomba de Aceite CHINT
120/220v
activara cuando la conexión sea creada por un electroimán. A diferencia del contactor, el relé permite
que los contactos se cierran o abran según los requerimientos deseados (Arjona, 2013). Para este
proyecto se utilizó 1 Relé auxiliar. Las especificaciones del mismo se muestran en Tabla 64.
242
Figura 127
Tabla 64
Programación
Arduino, de esa manera se cargará la programación del control que gobernará el control de todo el
sistema.
Librerías de Arduino
Figura 128
• HardwareSerial.h: Librería utilizada para enviar datos hacia la PC, desde un puerto serial.
Se diseña mediante el software Nextion Editor HMI, mostrando su interfaz principal en Figura
129, en la cual se presenta la página a ser diseñada, las páginas a usar en el proyecto, el bloque de
atributos ayuda a mejorar los diferentes aspectos utilizados en el HMI como los botones, fondos, etc. A
continuación, en el bloque de programación se adhiere a los atributos para poder direccionar los
botones, elegir un tiempo de cambio de página y muchas funciones más; el bloque de compilaciones
indica los posibles errores producidos al diseñar el HMI, se presenta también las páginas cargadas en el
proyecto y finalmente el bloque de herramientas, donde tenemos varias opciones para implementar en
Figura 129
En el HMI, para poder cambiar mediante un toque se implementan botones que redirigen a una
Figura 130
Adicionalmente para poder solicitar datos en el HMI, se adiciona un texto que es redireccionado
a una interfaz que tiene diferentes tipos de entradas como texto, números, símbolos chinos, etc.
Figura 131
De igual manera, existen botones que generan alguna acción como encender y apagar
descrito anteriormente.
Figura 132
Finalizado este proceso, se procede a cargar al microcontrolador las funciones necesarias para
cada opción de control como son los datos ingresados, acciones de actuadores y toma de datos.
Una vez finalizado el diseño, se carga la información editada en el HMI y queda estructurado de
la siguiente manera:
Pantalla principal: Presenta los datos principales del proyecto, adicionalmente permite pulsar el
Figura 133
Pantalla previa al menú principal: Presenta el nombre e imagen del equipo y el botón de
siguiente.
Figura 134
Figura 135
Y dependiendo de la opción escogida, lleva a las diferentes páginas para ingresar los datos ya
sea en flujo paralelo o en contraflujo (Figura 136), los mismos introducidos mediante la entrada
Figura 136
Figura 137
Una vez introducidos los datos de temperatura de agua y de flujos másicos a ser manejados en
la práctica, sigue a la siguiente interfaz donde se debe seleccionar el intervalo de tiempo para recolectar
los datos y los botones de acción de bombas y sistema de refrigeración. Cuando el equipo esté en
funcionamiento, se selecciona el botón de toma de datos que a partir de ese instante enviará los datos
tomados por las termocuplas y presentados en pantalla hasta que se opima de nuevo el botón “TOMA
DE DATOS”.
Figura 138
correctamente.
Rotolock de la línea de alta, en este caso al acumulador de líquido, comprobando que el acople de la
válvula de aguja se conecte perfectamente a la entrada de los manómetros. A su vez la línea de servicio
del manifold debe estar conectada a los manómetros de salida del recipiente de nitrógeno. Se abre la
aproximadamente 250 psi, una vez regulada la presión de salida se abre la válvula del manómetro de
alta para que regule 100 psi a la conexión con el sistema de refrigeración para una primera verificación.
constante, caso contrario existen fugas y el sistema debe ser analizado para encontrar fugas o averías. A
continuación, se calibra en la válvula del manifold para una mayor presión y se realiza el mismo
procedimiento, para la última prueba a presión máxima de 250 psi se debe dejar presurizado el sistema
cerrar el suministro del tanque de nitrógeno a su válvula de regulación y posteriormente se abre la llave
Figura 139
Carga de refrigerante. Existen acciones previas que deben ser llevadas a cabo para poder iniciar
Se realiza el procedimiento de vacío en el sistema con el fin de remover la humedad del mismo
tanto por el lado de alta como por el lado de baja debido a la presencia de una válvula solenoide, para
ello se trabajará en conjunto con una máquina de vacío que será conectada a la manguera de servicio
del manifold y al mismo tiempo las mangueras de baja y alta presión se conectan a las válvulas rotolock
del acumulador de líquido y del compresor. Al momento de encender la máquina de vacío se deberá
tomar en cuenta que ambas llaves en el manifold tendrán que permanecer completamente abiertas, el
Figura 140
refrigerante que se encontrara en posición vertical donde solo liberara gas. Desenroscando de manera
leve la misma manguera en el manifold. Esto debido a que luego de hacer el vacío es posible que haya
Se procede a la carga de refrigerante en fase liquida, para lo cual el manifold tendrá la siguiente
disposición, manguera de servicio al tanque en posición vertical girado 180°, manguera de baja a la
entrada del compresor y manguera de alta conectada al acumulador de líquido, se abre la válvula del
manómetro de alta presión del manifold con el objetivo de que fluya liquido al acumulador y controlar la
abertura hasta 80 psi aproximadamente hasta que se estabilice el sistema, manteniendo la válvula del
manómetro de baja cerrada, por otro lado debido a que el equipo cuenta con una válvula solenoide que
al ser nueva se encuentra en posición abierta y por lo tanto la manguera de baja es suficiente para llenar
253
al sistema entero. Se procede a abrir la llave de baja permitiendo el ingreso paulatino de gas
refrigerante R 404 a al sistema. Finalmente, se prende el equipo para que trabaje durante breves
Figura 141
En este apartado se realizan las conexiones de los componentes físicos con los elementos
eléctricos para que permita la comunicación con el control de la máquina. En la Figura 142 se puede
Figura 142
Nota: 1. Conexión sensores de control PID, 2. Conexión de termocuplas, 3. Conexión a reles del sistema
Una vez diseñada la placa con los elementos eléctricos necesarios para el control, se
En la sintonización del PID con la herramienta PID Tuner de Matlab, se obtuvo que la constante
Figura 143
presión y temperatura enviarán señales eléctricas por medio de los pines 34 y 35 del microcontrolador
Figura 144
expansión con la herramienta de Arduino serial Plotter, se realizaron pruebas de funcionamiento a las
mantiene en un rango entre el 20% y 40%, estabilizando al sistema en 7.6°C, para garantizar que el
Figura 145
De la Figura 145 se identifica a la línea azul como el comportamiento PID, la línea roja
refrigerante a la salida del evaporador. Se puede observar que existe una diferencia de tiempo de
estabilización comparada con la simulación ya que realmente se estabiliza a los 11 s después de haberse
prendido el compresor.
refrigerante a 7 °C.
257
Tabla 65
5,9 35.7
6,1 29.1
6,2 28
6,3 27.7
6,4 27.5
6,7 26.8
6,9 26.6
7,2 26.5
7,5 24.8
proporcionales de voltaje y permitir que la niquelina caliente el aceite mientras usa un control PID para
control PID, se observa como varia la temperatura con la ayuda de la herramienta serial Plotter de
Arduino.
258
Figura 146
De la Figura 146 se puede observar que la línea de color azul representa la temperatura del
aceite a 45°C, la línea de color rojo es el comportamiento de la temperatura del aceite dentro del
El control PID permite encender la niquelina a partir de que el aceite alcanza los 33°C en el
condensador y empieza a elevar la temperatura a 45°C, al llegar a esta temperatura el control deja de
enviar pulsos de voltaje para apagar la niquelina hasta que la temperatura del aceite sea menor a 45°C y
almacena el aceite se encuentra asilado con poliuretano y posee un sistema bypass por donde retorna
fluido caliente, ocasionando que la temperatura del aceite no disminuya en gran proporción, además, el
de todo el volumen del fluido contenido en el tanque, por lo que para bajar la temperatura en el aceite
se debería esperar alrededor de 8 horas y en ese momento la niquelina volvería a encender para ser
259
controlada por el PID, por tal razón, no es posible aplicar el control cuando el sistema trabaja enfriando
entrada y salida del intercambiador de calor en cada uno de los fluidos de trabajo, agua y aceite.
Se procede con las pruebas de funcionamiento, las bombas son encendidas y calibradas a un
temperatura objetivo en evaporador, el equipo tiene que llegar a condiciones equilibradas, por lo tanto,
cuando se lo haya conseguido al verificar que las temperaturas no varían de forma brusca, el equipo
contraflujo, aplicando flujos masicos constantes para ambos fluidos de trabajo. Las mediciones serán
Tabla 66
Figura 147
20.2
20
19.8
19.6
19.4
19.2
19
0 10 20 30 40 50 60 70
Tiempo (min)
Figura 148
43
42
41
40
39
38
37
0 10 20 30 40 50 60 70
Tiempo (min)
Análisis: Tanto en Figura 147 y Figura 148 se evidencia el funcionamiento del intercambiador de
calor de tubos concéntricos en disposición contraflujo en sus extremos, es decir las temperaturas de
De los resultados obtenidos en Tabla 66, se realizó el análisis de eficiencia del intercambiador de
Tabla 67
̇
Potencia Máxima 𝑸𝒎𝒂𝒙 Calor cedido 𝑸̇𝒄 Eficiencia
ORD
(W) (W) (%)
Figura 149
tubo concéntrico.
pruebas para intercambiadores tipo tubo concéntrico y multitubos para flujo paralelo y contraflujo
utilizando como fluidos de trabajo aceites API y agua para el Laboratorio del DECEM” debido al diseño y
construcción sin normas certificadas en el intercambiador de calor (Mejía & Pérez, 2011).
agua, la cual es el objetivo establecido en el interfaz HMI, para distintas condiciones con su respectivo
Tabla 68
Toma de datos de la temperatura de ingreso del agua a temperatura setting de 19,5 °C
Figura 150
20.2
4.00%
ERROR (%)
20
19.8 3.00%
19.6 2.00%
19.4
19.2 1.00%
19 0.00%
0 10 20 30 40 50 60 70
TIEMPO (MIN)
Tabla 69
Figura 151
6.00%
ERROR (%)
17.50 5.00%
4.00%
17.00 3.00%
2.00%
16.50
1.00%
16.00 0.00%
0 10 20 30 40 50 60 70
TIEMPO (MIN)
Tabla 70
Temperatura setting: 14 °C
Temperatura
Ord Tiempo Error Error (%)
𝑻𝒄𝒊
1 0 14.50 0.00 0.00%
Figura 152
6.00%
15.00
ERROR (%)
5.00%
14.50 4.00%
14.00 3.00%
2.00%
13.50 1.00%
13.00 0.00%
0 10 20 30 40 50 60 70
TIEMPO (MIN)
Tabla 71
Figura 153
12.40 5.00%
ERROR (%)
12.20 4.00%
12.00
11.80 3.00%
11.60 2.00%
11.40
1.00%
11.20
11.00 0.00%
0 10 20 30 40 50 60 70
TIEMPO (MIN)
En Tabla 72 se resumen los valores obtenidos para las cuatro pruebas realizadas.
Tabla 72
Resumen de resultados y análisis obtenidos en adquisición de datos.
Con el objetivo de analizar el ámbito económico del presente proyecto se realiza un análisis
financiero en el que se desglosan los costos directos e indirectos que interfirieron dentro del proyecto.
Costos directos
Se establece como costos directos para el presente proyecto los honorarios a profesionales,
sistema de control.
Honorarios a profesionales
Dichos costos son asumidos por la Universidad de las Fuerzas Armadas “ESPE” debido a que son
profesionales que durante el presente período académico se encuentran nómina como docentes, se
Tabla 73
Cargo Nombre Tiempo (horas) Valor hora (USD) Costo total (USD)
A partir de los siguientes elementos los costos fueron asumidos voluntariamente por los
Remuneración a no profesionales
Tabla 74
Ord Función Tiempo (horas) Valor hora (USD) Costo total (USD)
1 Técnico de Refrigeración 12 4 48
2 Técnico Mecánico 12 4 48
3 Técnico de Control 12 4 48
Remuneración de estudiantes
El trabajo realizado por los estudiantes y la remuneración reconocida por el mismo se detalla en
Tabla 75.
Tabla 75
Remuneración de estudiantes
Ord Nombre Tiempo (horas) Valor hora (USD) Costo total (USD)
A continuación, se describen los gastos realizados en cada elemento o sistema del proyecto.
Estructura Metálica
Tabla 76
Tabla 77
Sistema de Refrigeración
Tabla 78
Sistema de Bombeo
Tabla 79
Sistema de Control
Tabla 80
Costos Indirectos
Dentro del proyecto se estiman como Costos Indirectos los gastos utilizados en: Combustible,
Tabla 81
Costos Indirectos
Gastos Totales
Tabla 82
1 Estructura 110.56
Tabla 83
Tabla 84
Análisis de costos
3 SUBTOTAL 16815.05
Conclusiones
• Los ciclos de refrigeración teóricos muestran a cuatro elementos fundamentales (compresor, condensador,
importantes que cumplen con los requerimientos que se necesitan en la práctica, como se mostró en
Capitulo III.
• A partir del estudio térmico de un sistema de refrigeración por compresión de vapor, se determina que es
posible controlar la temperatura cuando el sistema solamente trabaja para refrigerar o para calentar un
medio circundante por separado, aplicando el control de dos posiciones juntamente con el control PID que
determina que la temperatura del condensador será dependiente del comportamiento del evaporador, por
• La mejor estrategia de control para enfriar el agua es trabajar en conjunto con el control de dos posiciones,
que controla al compresor para encender cuando necesita bajar la temperatura o apagar cuando ha llegado
a la temperatura requerida en el evaporador, por otro lado, el control PID aplicado en el sistema de
refrigeración permite controlar el flujo de refrigerante en el evaporador para enfriar el agua garantizando
ingresar al compresor.
• Comparando los métodos de control para aplicarlos en el sistema de refrigeración se determina que el
control PID es el más adecuado, debido a que se estabiliza a los 7s, que comparado con los controles P y PI
• Al realizar el análisis de eficiencia del intercambiador de calor de tubo concéntrico, proporcionado por el
Laboratorio de Conversión de la Energía, se obtuvo una eficiencia promedio de 2.35%. Valoración que se
justifica debido al diseño y construcción sin normas certificadas del intercambiador de calor utilizado.
• Al analizar los datos obtenidos en la experimentación con el equipo Banco de Pruebas, se concluye que el
equipo permite obtener una temperatura de ingreso de agua con un error porcentual promedio de 2.95%
inversamente proporcional a la temperatura de setting en la interfaz HMI, es decir mientras menor sea la
Recomendaciones
• Realizar un estudio termo – físico computarizado de los cambios de fase en el refrigerante durante todo el
proceso dentro del circuito de refrigeración, con el fin de determinar con mayor precisión su
(condensador) para poder implementar un control PID adecuado, teniendo variables influyentes en las
• Aplicar un sistema de enfriamiento para disipar el calor generado en el compresor, debido a que al realizar
un trabajo constante tiende a sobrecalentarse, causando averías y reduciendo la vida útil del compresor.
• Al realizar las prácticas en el banco de pruebas se debe evitar trabajar a una temperatura mayor de 55 °C en
el aceite, debido a que afecta considerablemente el funcionamiento del sistema de refrigeración, debido a
que el refrigerante que circula por el serpentín del condensador eleva su temperatura y por ende su presión
ocasionando paros del sistema de refrigeración por la acción del presostato de alta presión como medida de
Referencias
Ambato, P., & Bravo, J. (2021). Recuperación energética, puesta a punto, automatización y adquisición
de datos del equipo “Túnel de Aire Multipropósito”, del Laboratorio de Conversión de Energía [Tesis
[Link]
ASHRAE. (2013). Handbook. American Society of Heating, Refrigerating and Air - Conditioning Engineers.
[Link]
compresor/
Barletta, G., & Acevedo, O. (2021). Buenas prácticas en los procesos de instalación y mantenimiento de
Çengel, Y. A., & Ghajar, A. J. (2011). Transferencia de Calor y Masa (4th ed.). McGraw-Hill.
Danfoss. (2016). Forma de carga recomendada para sistemas de refrigeración con Compresores
[Link]
286
content/uploads/2015/03/PRA/DataBook pra_prn_50.pdf
[Link]
[Link]
[Link] [Link]
Flores, B. A. (2006). Control térmico en un ciclo de refrigeración por compresión de vapor [Tesis de
[Link]
del-producto-1
Hernández, R. (2010). Introducción a los sistemas de control: Conceptos, aplicación y simulación con
[Link]
287
Incropera, F. P., & DeWitt, D. P. (1999). Fundamentos de Transferencia de Calor (4th ed.). Pearson
Educación de México.
Kern, D. Q. (1999). Procesos de Transferencia de Calor (1st ed.). Compañia Editorial Continental.
Kreith, F., Manglik, R. M., & Bohn, M. S. (2012). Principios de Transferencia de Calor (7th ed.). Cengage
Learning.
La bodega eléctrica. (2020). ¿Qué es y para qué sirve un breaker o interruptor de circuito?
[Link]
sobrecarga.
pdf/view/415789/MAXIM/[Link]
Mejía, J., & Pérez, J. (2011). Ingeniería conceptual, básica y de detalle de un banco de pruebas para
intercambiadores tipo tubo concéntrico y multitubos para flujo paralelo y contraflujo utilizando
como fluidos de trabajo aceites API y agua para el Laboratorio del DECEM [Tesis - Carrera de
[Link]
Moran, M. J., & Shapiro, H. N. (2011). Fundamentals of Engineering Thermodynamics (7th ed.). John
Muñoz, L. S., & Zarantonelli, F. (2017). Análisis, diseño y construcción de una bomba de calor residencial
[Link]
[Link]
Patil, R. K., Shende, B. W., & Ghosh, P. K. (1982). Designing a helical-coil heat exchanger.
Pitts, D., & Sissom, L. (1998). Heat Transfer (2nd ed.). McGraw-Hill.
[Link]
[Link]
289
Apéndices