Tesis 2
Tesis 2
TRABAJO DE INVESTIGACIÓN
PRESENTADO POR:
HUANCAYO - PERÚ
2019
AGRADECIMIENTOS
2
DEDICATORIA
Gracias a Dios
A mi abuela, por su amor y apoyo incondicional.
A mi abuelo que me cuida y guía desde el cielo.
3
ÍNDICE DE CONTENIDOS
AGRADECIMIENTOS.....................................................................................................II
DEDICATORIA.............................................................................................................. III
ÍNDICE DE CONTENIDOS...........................................................................................IV
INDICE DE TABLAS....................................................................................................VIII
INDICE DE FIGURAS...................................................................................................IX
RESUMEN.................................................................................................................... XI
ABSTRAC.................................................................................................................... XII
INTRODUCCIÓN.........................................................................................................XIII
CAPÍTULO I.................................................................................................................. 14
PLANTEAMIENTO DEL ESTUDIO...............................................................................14
1.1. Planteamiento y formulación del problema.........................................................14
1.1.1. Problema General.......................................................................................14
1.1.2. Problemas Específicos................................................................................14
1.2. Objetivos............................................................................................................ 15
1.2.1. Objetivo general..........................................................................................15
1.2.2. Objetivos específicos..................................................................................15
1.3. Justificación e importancia.................................................................................15
1.3.1. Justificación Técnica...................................................................................15
1.3.2. Justificación Económica..............................................................................15
1.3.3. Justificación Teórica....................................................................................15
1.3.4. Justificación Ambiental................................................................................15
1.3.5. Importancia de la investigación...................................................................16
1.4. Alcances de la investigación..............................................................................16
1.4.1. Alcance temporal........................................................................................16
1.4.2. Alcance geográfico......................................................................................16
1.5. Limitaciones.......................................................................................................16
1.5.1. Limitaciones bibliográficas..........................................................................16
1.5.2. Limitaciones de información........................................................................16
1.6. Descripción de variables y Operacionalización de variables..............................17
CAPÍTULO II................................................................................................................. 19
MARCO TEÓRICO.......................................................................................................19
2.1. Antecedentes del problema................................................................................19
2.1.1. Artículos Científicos:...................................................................................19
4
2.1.2. Antecedentes nacionales:...........................................................................20
2.1.3. Antecedentes internacionales.....................................................................21
2.1.4. Estado del Arte...........................................................................................24
2.2. Bases teóricas....................................................................................................26
2.2.1. Introducción a la refrigeración.....................................................................26
[Link] Evaporador:...................................................................................................32
[Link] Tubería de succión:.......................................................................................32
[Link] Compresor.....................................................................................................32
[Link] Línea de líquido:............................................................................................34
[Link] Condensador:................................................................................................34
[Link] Recibidor:......................................................................................................34
[Link] Válvulas de servicio:......................................................................................35
[Link] Elemento de expansión:................................................................................35
[Link] Sistemas complementarios:...........................................................................35
[Link] Sistema de deshielo....................................................................................35
[Link] Deshidratador:.............................................................................................35
[Link] Indicador de humedad y líquido:..................................................................36
2.2.2 Refrigerantes....................................................................................................37
2.2.3 Refrigerante R134a:.........................................................................................39
2.2.4 Operaciones en los sistemas de refrigeración..................................................42
2.2.5 Recuperación de refrigerantes.........................................................................47
2.3 Definición de términos básicos...............................................................................58
2.3.1 Estándares Internacionales..............................................................................59
2.3.2 Nomenclatura:..................................................................................................59
CAPÍTULO III................................................................................................................ 62
METODOLOGÍA........................................................................................................... 62
3.1 Identificación de la necesidad problema e identificación de requerimientos).....62
3.1.1 Lista de requerimientos u exigencias...............................................................63
3.1.2 Estructura de funciones....................................................................................67
3.1.3 Análisis de la solución......................................................................................74
3.1.4 Determinación de los conceptos de solución óptimos......................................78
[Link] Concepto de solución N° 1............................................................................78
[Link] Concepto de solución N° 2............................................................................78
3.1.5 Concluir cuál es el concepto solución adecuado para el diseño.......................79
3.2 Propuestas de solución al problema o necesidad...................................................81
3.3 Pruebas de solución...............................................................................................82
5
3.3.1 Diseño de la Solución.......................................................................................82
3.7 Diseño experimental...............................................................................................96
3.8 Simulación.............................................................................................................. 96
CAPÍTULO IV............................................................................................................... 99
ANÁLISIS Y RESULTADOS.........................................................................................99
4.1 Resultados.............................................................................................................. 99
4.2 Análisis de resultados (Explicación de los objetivos, experimentos, ensayos,
pruebas, evaluaciones, seguimientos, estadística).....................................................102
4.2.1 Diseño factorial completo...............................................................................103
4.2.2 Optimización de respuesta: Respuesta..........................................................108
4.2.3 CONCLUSIÓN GENERAL:............................................................................110
CAPÍTULO V.............................................................................................................. 111
ANÁLISIS ECONÓMICO............................................................................................111
5.1 Costo de Inversión................................................................................................111
5.1.1 Costos Directos..............................................................................................111
5.1.2 Costos Indirectos............................................................................................116
5.1.3 Costo Total.....................................................................................................117
CAPÍTULO VI............................................................................................................. 118
CONSTRUCCIÓN.......................................................................................................118
6.1 Proceso de Construcción......................................................................................118
CONCLUSIONES.......................................................................................................121
TRABAJOS FUTUROS...............................................................................................122
REFERENCIAS BIBLIOGRÁFICAS............................................................................123
ANEXOS:.................................................................................................................... 125
ANEXO 1.................................................................................................................... 125
ANEXO 2.................................................................................................................... 127
ANEXO 3.................................................................................................................... 130
ANEXO 5.................................................................................................................... 136
6
INDICE DE TABLAS
Tabla 1. Operacionalización de variables.............................................................17
Tabla 2. Descripción de variable independiente...................................................17
Tabla 3. Descripción de variable dependiente......................................................17
Tabla 4. Propiedades críticas de algunos refrigerantes........................................38
Tabla 5. Punto de ebullición de algunos refrigerantes..........................................38
Tabla 6. Propiedades químicas del R134a...........................................................39
Tabla 7. Diferencias entre R134A y R12...............................................................40
Tabla 8. Especificaciones técnicas del refrigerante R134A..................................41
Tabla 9. Colores de los cilindros contenedores de refrigerante............................48
Tabla 10. Gases refrigerantes usados para sustituir a los CFCs..........................56
Tabla [Link] tabla 10 y 11......................................................................56
Tabla 12. Aplicación de los refrigerantes más comunes.......................................57
Tabla 13. Estándares Internacionales...................................................................59
Tabla 14. Nomenclatura........................................................................................59
Tabla 15. Lista de exigencias...............................................................................63
Tabla 16. Matriz de Funciones de la Máquina Recuperadora de Refrigerante.....67
Tabla 17. Síntesis de Funciones de la Máquina Recuperadora de Refrigerante...68
Tabla 18. Matriz Morfológica Maquina recuperadora de Refrigerante...................74
Tabla 19. Evaluación del concepto de solución – Yi.............................................79
Tabla 20. Evaluación del concepto de solución- Xi..............................................79
Tabla 21. Datos del valor económico y valor técnico............................................80
Tabla 22. Refrigerante r-134a saturado, tabla de Temperatura............................83
Tabla 23. Refrigerante r-134a saturado, tabla de Temperatura............................83
Tabla 24. Diseño Factorial....................................................................................96
Tabla 25. Costos de equipos..............................................................................111
Tabla [Link] del sistema eléctrico................................................................112
Tabla [Link] de los accesorios de refrigeración............................................113
Tabla 28. Costos de herramientas......................................................................114
Tabla 29. Costos de los materiales consumibles................................................114
Tabla [Link] de mano de obra directa..........................................................115
Tabla 31. Costo de Diseño.................................................................................115
Tabla [Link] directos....................................................................................115
Tabla [Link]áneos.........................................................................................116
Tabla 34. Costos de asesoramiento...................................................................116
Tabla 35. Costos indirectos................................................................................116
Tabla 36. Costo total...........................................................................................117
Tabla 37. Descripción del proceso de construcción del recuperador..................119
7
INDICE DE FIGURAS
Figura 1. Prototipo de recuperación de refrigerante..............................................19
Figura 2. Llenado del tanque de refrigeración.......................................................20
Figura 3. Prototipo Recuperador...........................................................................20
Figura 4. Sujesor de Motor...................................................................................21
Figura 5. Máquina recuperadora de refrigerante...................................................22
Figura 6. Banco de pruebas..................................................................................22
Figura 7. Recuperadora CPS TR700....................................................................23
Figura 8. Recuperación de refrigerante Automotriz...............................................23
Figura 9. Equipo de recuperación de refrigerante.................................................24
Figura 10. Spotter.................................................................................................24
Figura 11. Equipo de recuperación eléctrico.........................................................25
Figura 12. Sistema de Refrigeración.....................................................................26
Figura 13. Refrigeración por compresión de vapor...............................................27
Figura 14. a) Ciclo de Carnot inverso. b) Diagrama T vs. S de Carnot.................27
Figura 15. Diagrama Presión Entalpía..................................................................28
Figura 16. Diagrama real P-h................................................................................29
Figura [Link] constitutivos de un sistema de refrigeración.......................31
Figura 18. Evaporación del refrigerante................................................................32
Figura 19. Tipos de compresores.........................................................................32
Figura 20. Compresor Abierto...............................................................................33
Figura 21. Compresor Semi-hermetico.................................................................33
Figura 22. Compresor Hermético..........................................................................33
Figura 23. Condensación del refrigerante.............................................................34
Figura 24. Intercambiadores de calor de flujo cruzado. a) Con aletas y b) Sin
aletas.................................................................................................................... 34
Figura 25. Filtros deshidratadores........................................................................35
Figura [Link] interna de un filtro deshidratador.................................................36
Figura 27. Visor de humedad y líquido.................................................................36
Figura 28. Gráfica presión-temperatura R134a.....................................................40
Figura 29. Diagrama Presión Entalpia HFC-134a (Unidades SI)..........................42
Figura 30. Bomba de vacío...................................................................................42
Figura 31. Proceso de vacion en el sistema.........................................................43
Figura 32. Fuga de refrigerante............................................................................44
Figura 33. Diagrama para cargar gases zeotròpicos (serie 400)..........................44
Figura 34. Diagrama para cargar gases puros ò azeotrópicos..............................45
Figura 35. Manómetro para carga, descarga y recuperación de refrigerante........45
Figura 36. Funcionamiento del Visor de Liquido...................................................46
Figura 37. Medición del consumo de corriente del compresor..............................47
Figura 38. Cilindro de Recuperación.....................................................................48
Figura 39. Diagrama para recuperar gas con compresor apagado.......................49
Figura 40. Diagrama para recuperar gas con compresor encendido....................50
Figura 41. Diagrama para recuperar gas en fase gaseosa...................................51
Figura 42. Diagrama para recuperar gas en conexión Push/Pull..........................52
Figura [Link] para recuperar refrigerante en refrigerador..........................54
Figura 44. Diagrama para recuperar refrigerante en un sistema de aire
acondicionado.......................................................................................................55
Figura 45. Black Box.............................................................................................67
8
Figura 46. Caja Blanca de la Máquina Recuperadora de Refrigerante.................69
Figura 47. Fijación de procesos Técnicos.............................................................72
Figura 48. Representación final las funciones para la máquina recuperadora de
Refrigerante.......................................................................................................... 73
Figura 49. Alternativa de solución N°1..................................................................78
Figura [Link] de solución N°2...................................................................78
Figura 51. Solución lineal Valor Técnico Vs Valor Económico..............................80
Figura 52. Propuesta de Solución.........................................................................81
Figura 53. Grafica P-h del R-134a. Identificación de los puntos de operación......84
Figura 54. Grafica de eficiencia promedio para compresores reciprocantes.........88
Figura [Link] Guía rápida de Selección...............................................89
Figura 56. Compresor hermético de 1/5 hp...........................................................89
Figura 57. Diagrama Ph para el sistema de recuperación....................................90
Figura 58. Catalogo para selección del condensador...........................................91
Figura [Link] CDE 2782.......................................................................91
Figura 60. Resultado del Yate [1]..........................................................................99
Figura 61. Resultado del Yate a..........................................................................100
Figura 62. Resultado del Yate b..........................................................................100
Figura 63. Resultado del Yate b..........................................................................101
Figura 64. Diagrama de Pareto de efectos estandarizados................................102
Figura 65. Grafica normal de efectos estandarizados.........................................104
Figura 66. Graficas de residuos para respuesta.................................................105
Figura 67. Gráfica de cubos (medios ajustadas).................................................106
Figura 68. Grafica de efectos principales para Respuesta..................................107
Figura 69. Gráfica de interacción para respuesta...............................................107
Figura 70. Grafica de predicción.........................................................................108
Figura 71. Proceso de Construcción...................................................................118
9
RESUMEN
10
ABSTRAC
This technological research work includes the design and construction of the
TARIKAPUY refrigerant recovery machine R-134a Hydrofluorocarbons (HFC)
according to the suction / compression technique, whose main function is to
recycle the R134a refrigerant used in commercial refrigeration equipment from
refrigerators to air conditioning systems, which are excessively evacuated to the
environment after each corrective maintenance.
The machine aims to become a learning tool for current and new generations of
students of the Professional Academic School of Mechanical Engineering of the
Huancayo Continental University as a learning module; and at the same time
encourage commercial refrigeration companies to incorporate the use of
refrigerant recovery equipment in good refrigeration practices, generating a culture
of environmental protection avoiding the excessive release of depleting refrigerant
gases from the ozone layer to the atmosphere.
The machine has been designed based on thermodynamic calculations mainly for
the selection of the heart of the refrigerant recovery equipment, the reciprocating
compressor of a Tecumseh brand piston of 1 / 5HP of power. The specific
objectives are to determine how the recovery of refrigerant R134a in liquid and
vapor phase would be carried out.
Finally, a software analysis was performed to verify the machine's operating
parameters in the SolidWorks environment to specify the optimal operation of the
materials used by corroborating them by calculating the compressor selection,
using the Mollier diagram.
Finally, a software analysis was performed to verify the machine's operating
parameters in the CoolPack environment with the Engineering Equation Solver
extension, M-code programming (MATLAB) and finally the LabView software to
verify the optimal operation of the materials used Corroborating them by
calculating the compressor selection, using the Mollier diagram.
Keywords: R134a refrigerant, mechanical design, refrigeration, Mollier diagram.
11
INTRODUCCIÓN
Estudios realizados en el año 1974, advirtieron que el daño que estaba sufriendo
la capa ozono la cual protege el planeta de los rayos ultravioleta, debido a los
gases emanados desde la tierra. Permitieron detectar el peligro inesperado
causado por la emisión de clorofluorocarbonos conocidos como CFCs, los cuales
son usados en sistemas de refrigeración por compresión de vapor y ciclos
termodinámicos de potencia. Los resultados fueron contrastados hacia los años
80, por tal motivo en el año 1987 se firma el acuerdo internacional ambiental
Protocolo de Montreal con el objetivo de regular, reducir la producción y consumo
de sustancias agotadoras de la capa de Ozono (SAOs).[CITATION Jac09 \l 10250
]
La máquina recuperadora de refrigerante es un equipo frigorífico que condensa el
refrigerante extraído en fase gaseosa de cualquier maquina u equipo de
refrigeración y lo separa del aceite que pudiera llevar, realizando este trabajo en
varias fases (liquido subenfriado, liquido saturado, vapor saturado y vapor
sobrecalentado) [CITATION Fra12 \l 10250 ]
Entonces la maquina recuperadora consiste en en un sistema abierto el cual se
rige bajo el método de recuperación en fase de vapor, extrayendo el refrigerante
en su fase gaseosa, comprimiéndolo con la ayuda de un compresor reciprocante,
a través del área de succión o baja presión para luego ser condensado y
almacenado en un tanque de recuperación. Todo esto diseñado y construido de
tal manera que sus elementos sean accesibles para un mantenimiento y visibles
para que el operador pueda identificar cada uno de sus componentes. [CITATION
ISO \l 10250 ]
Al respecto se considera que el diseño de una maquina recuperadora de
refrigerante es una herramienta de aprendizaje muy importante para las actuales y
nuevas generaciones de estudiantes de la Escuela Académica Profesional de
Ingeniería Mecánica de la Universidad Continental; y al mismo tiempo incentivar al
colectivo de técnicos en refrigeración industrial comercial a incorporar el uso de
un equipo de recuperación de refrigerante a las buenas prácticas de refrigeración
generando así una cultura de protección al medio ambiente y evitando la
liberación desmedida de gases refrigerantes HCFs a la atmosfera.
El primero de enero del 2019, entro en vigor la encomienda de Kigali al protocolo
de Montreal, iniciativa para reducir “drásticamente” sustancias agotadoras de la
capa de ozono. El acuerdo contempla que, durante los próximos 30 años, se
reducirá en más de un 80 por ciento la producción y consumo de
hidrofluorocarbonos.[CITATION ONU03 \l 10250 ]
Justamente el presente trabajo de investigación tecnológica tiene como objetivo
central diseñar una maquina recuperadora de refrigerante R134a mediante la
técnica Succión/Compresión (considerando como indicadores el nivel de
recuperación de refrigerante en fase vapor, el nivel de recuperación de
refrigerante en fase líquido y el nivel de recuperación compresión/succión nivel de
en la maquina recuperadora de refrigerante R134a respectivamente.
12
CAPÍTULO I
PLANTEAMIENTO DEL ESTUDIO
13
1.2. Objetivos
1.2.1. Objetivo general
Diseñar la máquina TARIKAPUY recuperadora de refrigerante
R-134a según la técnica Succión/Compresión.
14
calentamiento global del planeta. Haciendo de la máquina eco
amigable con el medio ambiente.
1.5. Limitaciones
1.5.1. Limitaciones bibliográficas
Al pertenecer el tema de estudio a una investigación tecnológica
los recursos bibliográficos serán escasos.
15
1.6. Descripción de variables y Operacionalización de variables
16
recuperado con recuperación de Engineering Equation Solver
niveles de pureza refrigerante en la Flujo másico=Kg/h
aceptables. ARI700 fase Liquida.
Nivel de
recuperación Engineering Equation Solver
compresión/succión. Flujo másico=Kg/h
Fuente: Elaboración propia, 2019
17
CAPÍTULO II
MARCO TEÓRICO
18
Figura 1. Prototipo de recuperación de refrigerante
Fuente: Elaboración Carlos Orosco, 2004.
19
pero no
necesariamente con un alto grado de pureza”[CITATION
Jac09 \l 10250 ].
20
CESAR ANDRES ABARCA ERAZO, JOSE CARLOS ARTIGA
PORTILLOS Y YAMILETH MARISELA RIVERA CUBIAS
(2017). DISEÑO Y CONSTRUCCION DE UN EQUIPO DE
RECUPERACIÓN DE REFRIGERANTE R-134a. “Debido al
gran impacto que genera la liberación indiscriminada de gases
refrigerantes al ambiente, el uso de un recuperador de
refrigerante debe ser de carácter obligatorio dentro de las
buenas prácticas en refrigeración para un manejo adecuado y
responsable de estos gases refrigerantes presentes en diversos
sistemas; ya sea por mantenimiento o por la depuración de un
sistema, debido a un cambio de tecnología más reciente”.
[ CITATION Ces17 \l 10250 ]
21
Figura 6. Banco de pruebas.
Fuente: Elaboración Mero Jefferson, 2012.
22
Fuente: Elaboración Muñoz José, 2013.
23
Fuente: ICOR International
24
2.2. Bases teóricas
2.2.1. Introducción a la refrigeración
Una de las principales áreas de aplicación de la termodinámica es
la refrigeración, que es la transferencia de calor de una región de
alta temperatura hacia una región de baja temperatura. Los ciclos
en los que operan con mayor frecuencia, se denominan ciclos de
refrigeración por compresión de vapor, donde el refrigerante se
evapora y condensa alternadamente, para luego comprimirse en la
fase de vapor[ CITATION Ala99 \l 10250 ].
25
Figura 13. Refrigeración por compresión de vapor
Fuente: Revista Fic-Frio, 2019.
26
Considere un ciclo de Carnot invertido ejecutado dentro de la campana de
saturación de un refrigerante, como el que se muestra en la figura 4 a) y b).
1-2: Se transfiere (absorción) calor reversiblemente desde la región fría TL, de
forma isoterma donde el refrigerante experimenta cambios de fase[ CITATION
Gru13 \l 10250 ].
2-3: Se comprime el refrigerante isoentrópicamente, hasta que alcanza la
temperatura máxima TH[ CITATION Gru13 \l 10250 ].
3-4: Se transfiere calor reversiblemente a la región caliente a TH, de forma
isoterma, donde el refrigerante experimenta cambios de fase (vapor a líquido)
[ CITATION Gru13 \l 10250 ].
4-1: Se expande el refrigerante isoentrópicamente hasta, alcanzar la temperatura
mínima TL[ CITATION Gru13 \l 10250 ].
27
El proceso de expansión (proceso de 3-4):
Esto se produce en la válvula de expansión, cuya misión es disminuir la
presión y dosificar la cantidad de fluido refrigerante que debe ingresar al
evaporador a entalpia constante hasta el punto 4.
El proceso de evaporación (proceso de 4 -1):
La función que realiza el evaporador es extraer el calor del medio
refrigerado por medio de la evaporación del refrigerante. El evaporador
recibe refrigerante a baja presión y baja temperatura y lo entrega a baja
presión y alta temperatura”.[ CITATION Yun96 \l 10250 ]
28
(b-c’): El vapor es comprimido hasta una presión considerablemente
mayor que la presión de condensación, para forzar la salida de vapor a
través de la válvula de descarga.
(c’-2’): Caída de presión en la línea de descarga.
(3-3’): Zona de sub. Enfriamiento.
(4-4’): Refrigerante a menor temperatura y mayor concentración de
refrigerante líquido, aumentando el valor del efecto
refrigerante”[ CITATION Yun96 \l 10250 ].
[Link].3.2 Relación entre la presión y la temperatura:
Cuando se aumenta la temperatura debido a un suministro
de energía, la presión absoluta se aumentará en proporción
directa al aumento de la temperatura absoluta. Si el gas es
enfriado, la presión absoluta del gas disminuirá en
proporción directa a la disminución de la temperatura
absoluta[ CITATION Yun96 \l 10250 ].
Durante todo el ciclo de refrigeración la presión y
temperatura guardan una estrecha relación, según se
expone a continuación:
29
ciclo”.[ CITATION Fra06 \l 10250 ]
30
[Link] Evaporador:
“Superficie de transferencia de calor a través de la cual se absorbe el calor
del espacio refrigerado (el refrigerante se vaporiza)”[ CITATION Tor11 \l
10250 ].
[Link] Compresor
“Saca el vapor del evaporador para mantener la presión deseada de
vaporización a la temperatura deseada”[ CITATION Tor11 \l 10250 ].
a. Reciprocantes
b. Paleta
c. Tipo scroll
d. Tipo Tornillo
e. Turbocompresor
31
Figura 19. Tipos de compresores
Fuente: TECUMSEH, 2001
[Link].1 Abiertos
Pueden ser accionados por todo tipo de motor desde eléctricos
hasta los motores de combustión interna.
[Link].2 Semi-herméticos
Tiene al motor y compresor alojados en el mismo cuerpo dividido
en dos partes.
32
[Link].3 Herméticos
Son los compresores libres de mantenimientos, el mercado actual
de la refrigeración industrial tiene como favorito a este tipo de
compresor.
[Link] Condensador:
Es el elemento de la instalación que se encarga de pasar el estado de
vapor del fluido refrigerante a estado líquido. El fluido refrigerante ingresa
al condensador en estado gas (vapor recalentado) y sale en estado líquido
a la temperatura que se condenso o incluso a una temperatura menor si se
produce subenfriamiento[ CITATION Fra06 \l 10250 ].
Fuente: SEMARNAT, 2019.
33
[Link] Recibidor:
“Almacena el líquido condensado, de manera que este se pueda
suministrar continuamente, conforme a la necesidad del
evaporador”[ CITATION Tor11 \l 10250 ].
[Link] Deshidratador:
“La humedad es uno de los enemigos básicos de un sistema de
refrigeración y el nivel de humedad debe ser mantenido al mínimo con el
fin de evitar alteraciones en el funcionamiento del sistema o el deterioro
del compresor”[ CITATION Min17 \l 10250 ]. Los deshidratadores (Figura
18) o secadores, tal como se denominan comúnmente por cumplir su
función, está provisto de una malla a la entrada en forma cilíndrica y otra
malla a la salida en forma circular como se puede observar en la figura 9.
Están constituidos por una envoltura rellena con un secante o agente de
secado provista de un filtro en cada extremo[ CITATION Min17 \l 10250 ].
34
35
Figura 27. Visor de humedad y líquido.
Fuente: [EMERSON], 2019.
36
2.2.2 Refrigerantes
Los refrigerantes son los fluidos con elevado punto de ebullición y bajo punto de
fusión.
37
[Link]. Propiedades críticas de los refrigerantes:
Los refrigerantes normalmente tienen temperaturas críticas por encima de
90 C, pero hemos seleccionado algunos en la siguiente tabla (R404A,
R410A y el R744) cuyas temperaturas críticas están por debajo de ese
valor.[CITATION Uni05 \l 10250 ]
Tabla 4. Propiedades críticas de algunos refrigerantes
38
Fuente: ARI, 2019.
[Link] Componentes:
Tabla 6. Propiedades químicas del R134a
39
Tabla 7. Diferencias entre R134A y R12
40
Fuente: DuPont™ Suva® Refrigerantes, 2019
DESVENTAJAS
No se mezcla bien con los aceites minerales o alquilbencenos por lo
cual se tienen que sustituir (95% mínimo) por aceite polyester que es
más caro.
Tiene un deslizamiento de temperatura mínimo de 0.5 C[ CITATION
Doc19 \l 10250 ].
41
[Link] Diagrama presión entalpia del refrigerante r134a
Cada refrigerante tiene su propio gráfico de presión - entalpía. En el anexo 1, se
presenta en forma esquemática un diagrama de ese tipo. Allí la presión se
encuentra en el eje y, en escala logarítmica; en el eje horizontal se dibuja la
entalpía por kilogramo del refrigerante. Las líneas horizontales de la gráfica son
de presión constante (isóbaras) y las verticales son entonces de entalpía
constante (isoentálpicas).[ CITATION Esc13 \l 10250 ]
42
43
requerido, pero éste no se mantiene indica que aún existe una
fuga”[ CITATION Leo14 \l 10250 ].
Para localizar una fuga, el método más sencillo es cargar el sistema con
nitrógeno tanto en el lado de alta como en el de baja presión. Con una
mezcla de agua con jabón y ayuda de un pincel, se recorren las partes
44
en fase líquida se hace con el sistema apagado y la carga en fase de
vapor, se hace con el equipo funcionando[ CITATION Cet16 \l 10250 ].
Fuente: SEMARNAT, 2019
45
[Link] Determinación de la carga correcta de refrigerante:
Es primordial que el sistema cuente con la carga correcta de refrigerante,
lo que se determina considerando los siguientes aspectos:
46
Figura 37. Medición del consumo de corriente del compresor
Fuente: Frio Polar, 2010.
47
Que los cilindros tengan válvulas separadas para líquido y gas y estén
dotados de un dispositivo de alivio de la presión. En la figura 35 se puede
ver un cilindro típico de recuperación[CITATION Uni05 \l 10250 ].
48
[Link] Tecnologías de recuperación
Los cilindros recuperadores para refrigerante deben estar completamente
vacíos antes de proceder a cargarlos con gas[ CITATION Sec06 \l 10250 ].
Esto evita que el gas recuperado se contamine con aire, humedad o
remanentes del gas refrigerante que estuvo contenido anteriormente
[ CITATION Sec06 \l 10250 ]. Como antes se mencionó, se debe hacer un
vacío al tanque recuperador de al menos 1000 micrones [ CITATION
Sec06 \l 10250 ].
Para poder acelerar la recuperación de gas, se debe mantener frío el
tanque recuperador durante todo el proceso [ CITATION Sec06 \l 10250 ].
Esto se puede lograr colocándolo en una cubeta con hielo. Mientras más
frío esté el tanque, la presión del gas disminuye, pero si el equipo de
donde se está recuperando el gas está a una temperatura ambiente,
entonces el proceso de recuperado es más lento [ CITATION Sec06 \l
10250 ].
Antes de comenzar la recuperación de gas debe revisarse la posición de
todas las válvulas y, si aplica, constatar el nivel del aceite del compresor
de la recuperadora [ CITATION Sec06 \l 10250 ]. Es aconsejable recuperar
el gas refrigerante líquido en un tanque recibidor [ CITATION Sec06 \l
10250 ]. Debe recuperarse el líquido primero y después el vapor.
Recuperar el refrigerante en fase gaseosa deja aceite en el sistema,
minimizando la pérdida del mismo [ CITATION Sec06 \l 10250 ].
Si el compresor del sistema no funciona, hay que entibiar el cárter del
compresor [ CITATION Sec06 \l 10250 ]. Esto ayuda a liberar el
refrigerante atrapado en el aceite [ CITATION Sec06 \l 10250 ]. Figura 39.
49
Se tienen que instalar dos válvulas removibles, una para alta presión y otra para
baja [ CITATION Sec06 \l 10250 ]. El refrigerante migra y se condensa en el
tanque recuperador [ CITATION Sec06 \l 10250 ]. Con este método se recupera el
80% del gas y es aprobado por la EPA [ CITATION Sec06 \l 10250 ].
Recuperar el refrigerante de ambos lados, alta y baja, para poder
lograr un vacío completo [ CITATION Sec06 \l 10250 ]. Esta acción
también ayuda a acelerar el proceso de recuperación de gas
[ CITATION Sec06 \l 10250 ].
Si el compresor del sistema funciona, se debe encender y recuperar el
gas del lado de alta presión [ CITATION Sec06 \l 10250 ].
50
El refrigerante, en fase de vapor, es normalmente aspirado por la
succión de la máquina recuperadora y, una vez condensado, es
enviado al tanque recuperador [ CITATION Sec06 \l 10250 ].
Hay dos formas de conectar la máquina para recuperar vapor, según sea
el caso:
En el juego de manómetros de los dos lados del compresor. (Sistemas
comerciales medianos) [ CITATION Sec06 \l 10250 ].
Sólo del lado de baja, donde hay que instalar una válvula pinchadora
para extraer el refrigerante, y la cantidad a recuperar es pequeña
[ CITATION Sec06 \l 10250 ]. (Refrigeradores domésticos, aires
acondicionados de baja capacidad, congeladores pequeños)
[ CITATION Sec06 \l 10250 ].
51
En el caso de los enfriadores de agua industriales,
mantener agua circulando (lo cual, adicionalmente, previene
que se congele) [ CITATION Sec06 \l 10250 ].
En caso de que la máquina de recuperación no tenga un sistema
de evaporación, se debe proteger contra la llegada de refrigerante
líquido utilizando el juego de manómetros para ir dosificando,
mediante las válvulas de operación, su ingreso desde el sistema a
la máquina (utilizándolas como si fueran un dispositivo de
expansión) durante las etapas iniciales de la recuperación
[ CITATION Sec06 \l 10250 ].
El refrigerante líquido puede ser recuperado por técnicas de
decantación, separación o “push/pull” (succión y retroalimentación),
con el consiguiente arrastre de aceite [ CITATION Sec06 \l 10250 ].
52
El tanque recuperador succionará el refrigerante líquido (movimiento pull) de la
Unidad de Calefacción, Ventilación, Aire Acondicionado y Refrigeración (HVAC&R
– Heating, Ventilation, Air Conditioning and Refrigeration) desactivada, cuando la
máquina recuperadora haga disminuir la presión del cilindro [ CITATION Sec06 \l
10250 ]. El vapor succionado del tanque recuperador por la máquina recuperadora
será entonces empujado de vuelta (movimiento push), es decir, comprimido hacia
el lado que corresponde al vapor en la unidad de HVAC&R desactivada
[ CITATION Sec06 \l 10250 ].
Una vez que la mayoría del refrigerante haya sido cargado del sistema al tanque
recuperador, la maquina recuperadora comenzará a ciclar, controlada por su
presostato de baja presión de succión, removiendo el resto del refrigerante en
forma de vapor [ CITATION Sec06 \l 10250 ]. Cuando la máquina de recuperación
ya no continúe ciclando y se detenga por completo, eso indica que se ha
recuperado todo el refrigerante posible del sistema [ CITATION Sec06 \l 10250 ].
No se debe utilizar el método “push/pull”:
Si el sistema o equipo tienen una carga menor de 9 kilos ó 20 libras, de
gas refrigerante [ CITATION Sec06 \l 10250 ].
Si el equipo es una bomba de calor u otro sistema en donde el refrigerante
líquido pudiera quedar aislado [ CITATION Sec06 \l 10250 ].
Si el equipo tiene un acumulador entre los puertos de servicio, utilizados
para recuperar líquido [ CITATION Sec06 \l 10250 ].
Si ha ocurrido una migración de refrigerante líquido, y se desconoce su
ubicación [ CITATION Sec06 \l 10250 ].
Si el diseño de la tubería en el equipo no permite crear una columna sólida
de líquido [ CITATION Sec06 \l 10250 ].
Si se utiliza el método “push/pull:
Se necesita una mirilla, para poder saber que se terminó de recuperar todo
el líquido [ CITATION Sec06 \l 10250 ].
Tener una tercera manguera lista, ya que será necesaria [ CITATION
Sec06 \l 10250 ].
Después de haber retirado todo el líquido, se deben reconfigurar las
mangueras para recuperar vapor, ya que este método no hace un vacío en
el sistema [ CITATION Sec06 \l 10250 ].
53
10250 ].Para esto, hay que instalar una válvula punzonadora en el
sistema, siguiendo las instrucciones del fabricante, y utilizar una unidad de
recuperación para extraer el refrigerante de la unidad mediante un injerto
de toma de línea al igual que con un sistema mayor [ CITATION Uni05 \l
10250 ].Las válvulas punzonadoras nunca deben dejarse instaladas de
modo permanente, sino que hay que retirarlas después de su uso si están
instaladas en el tubo de proceso [ CITATION Uni05 \l 10250 ].En la gráfica
4-7, la unidad de recuperación está conectada al refrigerador mediante
una válvula punzonadora típica [ CITATION Uni05 \l 10250 ].Debido a que
la carga de refrigerante es pequeña, sólo hace falta recuperar gas
[ CITATION Uni05 \l 10250 ].Se recomienda instalar válvulas
punzonadoras en ambos lados de presión [ CITATION Uni05 \l 10250 ].
Fuente: SEMARNAT, 2019
54
Figura 44. Diagrama para recuperar refrigerante en un sistema de aire
55
Tabla 10. Gases refrigerantes usados para sustituir a los CFCs.
56
[Link]. Tabla de aplicación de los gases más comunes
Tabla 12. Aplicación de los refrigerantes más comunes.
57
2.3 Definición de términos básicos
58
que el equipo de vacío para la recuperación final y no superior a 100 kPa
(absoluta)[ CITATION ISO \l 10250 ].
2.3.2 Nomenclatura:
SIMBOLO NOMBRE
T Temperatura
S Entropía
TL| Baja Temperatura
TH Alta Temperatura
H Entalpia
P Presión
Q Carga térmica
Qe Carga térmica por estructura
Qc Carga por cambio de aire
Qv Cargas térmicas varias
Qt Carga térmica total
U Coeficiente global de transferencia de calor
A Superficie de transferencia de calor
∆T Diferencial de temperatura
59
1/Mi Coeficiente de convección de pared interior
ex Espesor de la pared
Kx Coeficiente de conductividad térmica
1/Mo Coeficiente de convección de pared exterior
C Conductancia térmica
K Conductividad térmica
X Espesor de material en pulgadas
R Resistencia que posee cada material
Hr Humedad relativa
Vi Volumen interior de la cámara
Fv Factor de cambio de aire en un promedio de 24 hr del espacio
refrigerado
Fh Factor por cambios de aire
q Carga térmica por producto
q1 Carga por producto antes del punto de congelación
q2 Carga del producto al fundirse el hielo
q3 Calor sensible de fusión
m masa del producto
Mtg Diferencia de entalpías de saturación
tr Temperatura de rocío
Ti Temperatura interna
Tevap Temperatura de evaporación
Tr Temperatura del espacio refrigerado
Dt Diferencia de temperatura entre el evaporador y espacio refrigerado
Tcond Temperatura de condensación
Tm Temperatura del medio ambiente
Pevap Presión de evaporación
Pcond Presión de condensación
RC Relación de compresión
Qf Potencia frigorífica
m Caudal Masico
Vc Caudal volumétrico
y Razón de compresión
ŋm Rendimiento mecánico del compresor
ŋv Rendimiento volumétrico del compresor
ŋT Rendimiento total del compresor
Wc Potencia del compresor
Qc Potencia calorífica evacuada en el condensador
V Caudal volumétrico de refrigerante entrando en la válvula de expansión
COP Coeficiente de rendimiento del sistema de refrigeración
Fuente: Elaboración propia
60
61
CAPÍTULO III
METODOLOGÍA
Fecha:2/05/19
Características D Exigencias Responsable
E
GEOMETRÍA
Donde:
a. Compresor hermético
b. Moto ventilador
c. Separador de aceite
d. Cañería de cobre
e. Válvulas de alta y baja presión
f. Condensador
g. Presostato
h. Válvula de Gusanillo
i. Manómetros
63
D 2. La unidad de recuperación de Y. Martinez
GEOMETRÍA refrigerante R134a debe tener
como componentes esenciales
al compresor y condensador, su
disposición será en forma de
espejo uno al frente del otro.
FABRICACIÓN
E 8. El consumo de corriente Y. Martinez
eléctrica debe ser mínimo.
FABRICACIÓN
E Y. Martinez
9. La potencia del compresor
hermético debe ser de ½ hp
FABRICACIÓN
E 10. La máquina recuperadora de Y. Martinez
refrigerante R-134a debe
funcionar con suministro
monofásico.
64
FABRICACIÓN
D 12. El refrigerante recuperado debe Y. Martinez
cumplir con los estándares de
calidad adecuados para su
reutilización.
ERGONOMÍA
D 13. La aplicación máquina- Y. Martinez
ambiente- humano, no debe
infligir daños en el técnico al
momento de realizar el servicio
de mantenimiento.
USO
E 15. El servicio que prestará la Y. Martinez
maquina recuperadora de
refrigerante será para trabajos
de recuperación de refrigerante
en maquinaria pesada.
MONTAJE
E 16. La máquina recuperadora de Y. Martinez
refrigerante R134 será de fácil
acoplamiento, de manera que
las actividades del
mantenimiento no demanden
mucho tiempo.
D 17. El estado de tecnología de estas Y. Martinez
máquinas recuperadoras de
FABRICACIÓN refrigerante, se obtuvo a través
de catálogos de fabricantes de
recuperadoras de refrigerante
de con capacidades de
recuperación desde CFC hasta
HCFC. No ha sido posible
encontrar planos, pero si un
poco de bibliografía referido al
tema.
65
E [Link] fecha de entrega del
diseño de la maquina
USO recuperadora de refrigerante
R134a, que comprende los Y. Martinez
cálculos efectuados y los
planos, será para fines del
mes de noviembre de 2019.
D: Deseo
E: Exigencia
66
3.1.2 Estructura de funciones
En el siguiente apartado se implementará la Caja Negra: “Black Box” la cual
consiste en un conjunto de entradas y salidas fundamentada en la conservación
de la energía[ CITATION Cri17 \l 10250 ].
67
[Link] Abstracción: Caja negra (Black-Box) de la maquina
recupera de refrigerante.
• Temperatura de
Descarga.
• Temperatura • Flujo masico de
de Succión. Máquina refrigerante.
Recuperadora • Potencia del
• Presiones de compresor.
Evaporación y de Refrigerante • Estados
condensación. Termodinámicos
• COP.
Varia
Varia ble
ble Depe
Inde ndien
pendi Figura 45. Black Box te
ente
Fuente: Elaboración Propia, 2019.
ELEMENTO FUNCIÓN
Bombea y comprime el fluido
Compresor refrigerante, creando una zona de
alta y baja presión.
Intercambiar calor; en el cual el
Condensador refrigerante succionado ingresa en
forma de vapor y sale en forma
líquida.
Separador de Aceite Separa el aceite del refrigerante y lo
devuelve al cárter del compresor.
Dispositivo encargado de realizar la
Moto-Ventilador convección forzada y acelerar el
proceso de condensación.
Filtro deshidratador Material filtrante que remueve la
humedad del sistema.
Manómetro de Baja Presión Dispositivo diseñado para medir la
68
presión de succión del compresor.
Manómetro de Alta Presión Dispositivo diseñado para medir la
presión de descarga dl compresor.
Visor de Liquido-Humedad Permite observar la condición del
refrigerante.
Abre o cierra el paso del flujo de
Válvula solenoide refrigerante es accionada
eléctricamente.
Dispositivo de seguridad, cierra o
Presostatos de Alta y Baja abre circuito eléctrico según la
lectura de presión del fluido en el
sistema.
Fuente: Elaboración Propia,2019.
69
[Link] Caja Blanca de la Máquina Recuperadora de Refrigerante
ii. Verificar que las dos válvulas de servicio estén cerradas y proceder a
conectar las mangueras de la siguiente manera:
a) Uno de los extremos de la manguera de color azul a la válvula de servicio
del mismo color, es decir, en el lado de baja presión, roscando la tuerca
tipo Flare.
b) Uno de los extremos de la manguera de color rojo a la válvula de servicio
del mismo color, es decir, en el lado de alta presión, roscando la tuerca
tipo Flare.
vii. Manteniendo aún cerrada la válvula de servicio de alta presión (color rojo)
y abrir lentamente la válvula de servicio de baja presión (color azul). En
este instante comenzará a ingresar el refrigerante R-134a en el sistema.
Se podrá visualizar en el manómetro de baja, como la presión comienza a
ascender conforme ingresa el gas refrigerante al sistema.
xi. Una vez recuperado en su totalidad el refrigerante del sistema, y que este
se haya almacenado en el tanque, apague el equipo mediante el
interruptor selector.
xii. Como último paso, se debe realizar la purga del aceite acumulado en la
trampa de aceite para garantizar que éste retorne hacia el cárter del
compresor.
72
[Link] Fijar los procesos técnicos
73
[Link] Representación de la estructura de funciones
SUBFUNCIONES
01 02 03 04 05 06
FUNCION ÓPTIMA
01+02+03 04 05+06
REPRESENTACIÓN
SISTEMA
HVAC&R
SIN
REPARAR
UNIDAD
DE
PROCESO
REFRIGERANTE RECUPERADO
Figura 48. Representación final las funciones para la máquina recuperadora de Refrigerante
74
75
3.1.3 Análisis de la solución
[Link] Determinación del concepto de soluciones (matriz morfológica).
SELECCIÓN
DE
MANGUERAS
EXTRACCION
DEL
REFRIGERANTE
SISTEMA DE
VENTILACION
ALIMENTACIO
N DE
ENERGIA
ELECTRICA
SISTEMA DE
CONDENSACI
ON
77
SISTEMA DE
SEPARACIÓN
DE ACEITE
SOPORTE O
BASTIDOR
SISTEMAS
DE
SEGURIDAD
78
ENVASADO
CONCEPTOS
DE SOLUCIÓN
79
3.1.4 Determinación de los conceptos de solución óptimos
Luego de realizar la matriz morfológica se determinó dos posibles soluciones
(Conceptos de soluciones)
80
3.1.5 Concluir cuál es el concepto solución adecuado para el diseño
g1 . p1 + g2 . p2 +…+ g n . pn
y i= ≤1
( g1 + g2 +…+ g n ) . pmax
81
8 Calidad 3 2 6 3 9 4 12
Puntaje máximo Ʃ 29 65 75 116
Valor económico xi 0,75 0,86 1
g 1 . p1+ g 2 . p2 +…+ gn . p n
x i= ≤1
( g1 + g2 +…+ g n) . pmax
Solución Ideal
1.2
1
Valor Económico Y
0.8
0.6
0.4
0.2
0
0.7 0.75 0.8 0.85 0.9 0.95 1 1.05
Valor técnico X
82
Fuente: elaboración propia, 2019.
DEFINICION DE LA SOLUCION:
La solución elegida presenta las siguientes características: es ergonómica,
facilidad para su transporte y manipulación. La selección de sus
componentes mediante cálculos termodinámicos y el uso de tablas se
lograrán contrastar con la ayuda del software Cool Pack y la extensión
EEQ (Engineering Equation Solver); los cuales juegan un papel
fundamental al momento de optimizar el funcionamiento de la máquina.
83
Como se mencionó previamente se considera como
punto inicial la presión de succión del sistema como la
presión a la cual se encuentra el refrigerante en el equipo
al cual se le realizara la recuperación.
84
Tabla
Tabla 23. Refrigerante 22. Refrigerante
r-134a r-134a
saturado, tabla saturado, tabla de Temperatura
de Temperatura
85
Figura 53. Grafica P-h del R-134a. Identificación de los puntos de operación
Estado 3:
Para este estado se necesita conocer la temperatura ambiente a la cual se
realizará la recuperación ya que la presión de condensación está directamente
ligada a la temperatura del refrigerante.
Figura 54: Grafica P-h del R-134a. Identificación de los puntos de operación
del sistema recuperador de refrigerante.
T 3=Tamb+∆ T
Siendo:
∆ T =20 ° F
Tamb=86 ° F
86
Entonces:
r 3=106.0 ° F
BTU
P 3=116.83
lbm
87
Tabla 25. Refrigerante r-134a sobrecalentado25
88
Una vez determinadas las características termodinámicas de cada uno de los
estados que intervienen en el proceso de recuperación, se obtienen otros
parámetros de importancia para el dimensionamiento del equipo.
𝑊𝑐𝑜𝑚𝑝 = ℎ2 − ℎ1
𝑊𝑐𝑜𝑚𝑝 = (119.18 𝐵𝑇𝑈⁄𝑙𝑏𝑚) − (107.24 𝐵𝑇𝑈⁄𝑙𝑏𝑚)
𝑊𝑐𝑜𝑚𝑝 = 11.94 𝐵𝑇𝑈⁄𝑙𝑏𝑚
CALOR RECHAZADO:
𝑄 = ℎ2 − ℎ3
𝑄 = (119.18 𝐵𝑇𝑈⁄𝑙𝑏𝑚) − (116.83 𝐵𝑇𝑈⁄𝑙𝑏𝑚)
𝑄 = 2.35 𝐵𝑇𝑈⁄𝑙𝑏𝑚
𝑃𝑇 = (𝑚̇)(𝑊𝑐𝑜𝑚𝑝)
Y conociendo que:
1 𝐻𝑃 = 2544.4 𝐵𝑇𝑈⁄ℎ𝑜𝑟𝑎
Tenemos:
Para conocer la potencia nominal que se requiere como unidad motora del
sistema se utiliza la eficiencia total de un compresor reciprocante, y se toman
los datos de la figura 55.
89
Figura 54. Grafica de eficiencia promedio para compresores reciprocantes
Fuente:
RELACIÓN DE COMPRESIÓN:
𝜂𝑇 = 0.75
Es decir:
𝜂𝑇 = 75 %
𝑛 = 1750 𝑟𝑝𝑚
POTENCIA UTIL:
𝐵ℎ𝑝 =𝑃𝑇/𝜂𝑇
𝐵ℎ𝑝 =0.141 𝐻𝑃/0.75
𝐵ℎ𝑝 = 0.188 𝐻𝑃
90
𝐵ℎ𝑝 = 0.20 𝐻𝑃 = 1/5 HP
91
[Link].1.3 Selección del Condensador
Para seleccionar el condensador se debe conocer
el calor que se debe despejar, la temperatura de
condensación y la diferencia de temperatura con
el ambiente. Para hallar el calor a despejar se
deben hallar las propiedades en el estado 2 de
acuerdo a la Figura 57. Como no se tiene ningún
dato más que la presión de 2, se asumirá que el
compresor realiza un proceso isoentrópico, es
decir s1=s2, luego con esa propiedad, y haciendo
uso de tablas termodinámicas se puede hallar la
entalpía isoentrópica en el estado 2. De
considerarse el caso en que el compresor no
fuese isoentrópico, la entalpía alcanzada en 2
sería menor puesto que el trabajo entregado al
compresor no sería aprovechado solamente por
el fluido refrigerante sino también se perdería
como calor cedido al medio ambiente y calor por
fricción. Es por ello que al asumir un proceso de
compresión isoentrópico la variación de entalpía
de 2 a 1 será mayor, y por ello el calor que
deberá evacuar el condensador también será
mayor, es decir se establecen condiciones críticas
de operación.
kg kj
vcond=mt∗( h 2−h 3 ) =0.0067 ∗( 303.47−115.8 ) =1.26 kw
s kg
Se opta por trabajar con condensadores de la
marca brasilera Elgin. De acuerdo a su catálogo,
las capacidades están dadas para una
92
temperatura de condensación de 45ºC y una
diferencia de temperaturas de 35ºC, en el lugar
de operación del equipo, Lima, la temperatura
ambiental es menor lo cual favorece a una mejor
transferencia de calor. Las capacidades están
dadas en kcal/h por lo cual se realiza la
conversión de unidades:
93
Fuente:[Link]
atalogo%202013%20(Separado)/6.0%20-%[Link]
Con los Datos obtenidos de las presiones de alta y baja, así como
de las temperaturas nos ubicamos en el Diagrama de Moliere:
Las presiones de alta y baja las localizamos en la parte
derecha elegimos adecuadamente los puntos y trazamos
dos líneas horizontales las cuales cortaran la gráfica de
(Presión vs Entalpia).
94
De los puntos localizados serán enumerados desde la parte
izquierda respectivamente en sentido anti horario 1, 2, 3, 4.
Nos servirán para encontrar las entalpias las cuales serán
usadas en nuestros cálculos posterior mente. Las
localizaremos trazando líneas verticales a los puntos.
Qevaporador=(h 3−h2)
402 kJ 201 kJ
Qevaporador=( − )
kg kg
95
Qevaporador=201 kJ /kg
ṁ=PF / Qevaporador
PF: Potencial Frigorífico. Se obtiene del Abaco psicométrico
Kg/s
También como:
ṁ:flujo de masa=26 gramos /seg
Qevaporador=201 kJ /kg
PF: Potencial frigorífico
ṁ=PF / Qevaporador
PF=ṁ x Qevaporador
PF=0.026 kg /segx 201 kJ /kg
PF=5.226 kJ /seg
c. Trabajo de compresión
Qcompresor=(h 4−h3)
434 kJ 402 kJ
Qcompresor=( − )
kg kg
Qcompresor=402 kJ /kg
d. Potencia eléctrica del compresor
Pec=Qcompresorx ṁ
402 kJ
Pec= x 0.026 kg /seg
kg
Pec=10.452 kJ / seg
e. Caudal volumétrico
V =ṁ xve
Ve: Volumen especifico R134a
ve :volumen especifico R 134 a=0.235 m3 /k
0.026 kg
V= 0.235 m3 /kg
segx
96
97
f. V =0 .00611 m 3/seg∨¿ Calor a disipar del condensador
Qcondensador=h 4−h 1
434 kJ 201 kJ
Qcondensador=( − )
kg kg
Qcondensador=233 kJ /kg
94
Qev=5.226 kJ /s
m. Temperatura de entrada y salida de la válvula de expansión
T 1=0.1° C
T 2=8.5 ° C
n. COP del sistema
Qevaporador
COP=
Wcompresor
5.226
COP=
2.0262
COP=2.46766
Finalmente:
o. Capacidad de recuperación:
Qo Qo kg
ṁ= = =0.0960
ε (h 4−h 1) s
ṁ=345.6 kg/h
p. Tiempo de recuperación:
m
t=
ṁr
95
Diseño factorial 22 o 2x2
Tabla 24. Diseño Factorial
3.8 Simulación
3.8.1 Simulación con aplicación grafica visual(LabView-Student)
Ver código de programación en el anexo 5 del presente trabajo de investigación.
96
97
98
CAPÍTULO IV
ANÁLISIS Y RESULTADOS
4.1 Resultados
Primer resultado obtenido para el coeficiente de performance en base al programa
CoolPack-Engineering Equation Solver
Fuente: CooLPack
99
Figura 61. Resultado del Yate a
. Fuente: CooLPack
Fuente: CooLPack
100
Figura 63. Resultado del Yate b
Fuente: CooLPack
101
4.2 Análisis de resultados (Explicación de los objetivos,
experimentos, ensayos, pruebas, evaluaciones, seguimientos,
estadística)
*Respuesta COP
AB
0 5 10 15 20
Efecto estandarizado
102
4.2.1 Diseño factorial completo
Resumen del diseño
Factores 2 Diseño de la base: 2; 4
:
Corridas: 8 Réplicas: 2
Total 7 0.139092
103
El resumen de modelo muestra niveles altos de predicción, los
cuales son contrastados con la tabla de Anova y una desviación
estándar igual a 0.0173421.
Coeficientes codificados
Términ Efecto Coef EE del Valor T Valor p FIV
o coef.
A - - 0.00 - 0. 1.
0.119 0.059 613 9.7 00 0
50 75 4 1 0
B - - 0.00 - 0. 1.
0.218 0.109 613 17. 00 0
00 00 78 0 0
Estructura de alias
Facto Nombre
r
A A
B B
104
Gráfica normal de efectos estandarizados
(la respuesta es Respuesta; α = 0.05)
99
Tipo de efecto
Significativo
95 Factor Nombre
90 A A
B B
80 AB
70
Porcentaje
60
50 A
40
30
20 B
10
1
-20 -15 -10 -5 0 5 10
Efecto estandarizado
Fuente: Minitab 18
105
Gráficas de residuos para Respuesta
Gráfica de probabilidad normal vs. ajustes
99
0.02
90
0.01
Porcentaje
Residuo
50 0.00
-0.01
10
-0.02
1
-0.030 -0.015 0.000 0.015 0.030 2.3 2.4 2.5 2.6
Residuo Valor ajustado
1.5 0.01
Frecuencia
Residuo
1.0 0.00
-0.01
0.5
-0.02
0.0
-0.02 -0.01 0.00 0.01 0.02 1 2 3 4 5 6 7 8
Residuo Orden de observación
106
Gráfica de cubos (medias ajustadas) de Respuesta
2.2920 2.2570
1
2.5945 2.3905
-1
-1 1
A
107
Gráfica de interacción para Respuesta
Medias ajustadas
-1 1
B*A A
2.6
-1.0
1.0
2.5
Media de Respuesta
2.4
2.3
A* B B
2.6
-1.0
1.0
2.5
2.4
2.3
-1 1
A B
2.45
Media de Respuesta
2.40
2.35
2.30
-1 1 -1 1
108
Óptima A B
Alto 1.0 1.0
D: 0.9387
Act [-1.0] [-1.0]
Predecir Bajo -1.0 -1.0
Respuest
Máximo
y = 2.5945
d = 0.93869
Solución
Solución A B Respuesta Deseabilidad
Ajuste compuesta
1 -1 -1 2.5945 0.938692
109
Variabl Valor de
e configuración
A -1
B -1
Configuración
Variabl Valor de
e configuración
A -1
B 1
Predicción
Ajust EE de ajuste IC de 95% IP de 95%
e
110
2.292 0.0122627 (2.25795; 2.32605) (2.23303; 2.35097)
Configuración
Variabl Valor de
e configuración
A 1
B -1
Predicción
Ajuste EE de ajuste IC de 95% IP de 95%
111
CAPÍTULO V
ANÁLISIS ECONÓMICO
5.1 Costo de Inversión
5.1.1 Costos Directos
Los costos directos se derivan de los materiales, mano de obra,
herramienta, maquinaria o labores directamente involucradas en la
construcción física del proyecto. A continuación, se detalla los
costos directos:
SOLDADURA TIG
Modelo: TIG 170 (60%) A 4mm 5kg KVA5
1 S/0.00 S/0.00
POTENZA TIG 170 HF
Marca: WELD WELL
SOLDADURA ELECTRICA
Capacidad: 200A 1 S/0.00 S/0.00
Marca: WELD WELL
TOTAL S/0.00
Fuente: elaboración propia, 2019
112
Tabla [Link] del sistema eléctrico
VALOR
DESCRIPCIÓN CANTIDAD VALOR TOTAL
UNITARIO
Cable Surperflex 4x12 AWG 1Kv C/M CA 3 S/8.00 S/24.00
Cable Flexible # 16 5 S/10.00 S/50.00
Cable Flexible # 18 3 S/12.00 S/36.00
Cable flexible # 12 2 S/6.00 S/12.00
Luz piloto 11 S/2.00 S/22.00
Selectores 10 S/2.50 S/25.00
Adhesivos plásticos pequeños 2 S/1.00 S/2.00
Bloque de pulsador 1 S/5.00 S/5.00
Bornera plástica 6A 18 AWG 5 S/2.50 S/12.50
Bornera plástica 15A 14 AWG 1 S/1.50 S/1.50
Cinta ordenadora 1 S/5.50 S/5.50
Bornera para Relay 8 Pin Plano PTF-08AE 7 S/1.50 S/10.50
Cable multipar 5 S/0.40 S/2.00
Toma de 3 polos 50A 220V 1 S/6.00 S/6.00
Enchufe trifásico 50 Amp. 250V 1 S/5.50 S/5.50
Botonera on/off 15 Amp. 1 S/4.50 S/4.50
Interruptor iluminado 3 S/1.50 S/4.50
Riel de acero 1 S/2.50 S/2.50
Terminal hembra de alta temperatura 1 S/4.80 S/4.80
TOTAL S/235.80
Fuente: elaboración propia, 2019
113
Tabla [Link] de los accesorios de refrigeración
VALOR
DESCRIPCIÓN CANTIDAD VALOR TOTAL
UNITARIO
Tubo flexible de 3/8" 13 S/3.70 S/48.10
Tubo flexible de 1/4" 3 S/2.90 S/8.70
Tubo flexible 1/2" 3 S/5.80 S/17.40
Tubo flexible de 5/8 8 S/5.20 S/41.60
Válvula de paso de 3/8" 5 S/9.00 S/45.00
Filtro secador de 3/8" 1 S/5.50 S/5.50
Filtro secador de 1/4" 1 S/6.80 S/6.80
Tuercas de bronce de 3/8" 35 S/0.60 S/21.00
Tuercas de bronce de 1/4" 6 S/0.70 S/4.20
Tee 3/8 4 S/4.90 S/19.60
Conector macho 3/8 x 3/8 8 S/2.60 S/20.80
Codos 5/8 4 S/1.40 S/5.60
Codos ½ 4 S/1.30 S/5.20
Codos 3/8 1 S/0.70 S/0.70
Válvulas de carga y 5
S/2.50 S/12.50
accesorios
Manómetro de alta QG 2
S/4.70 S/9.40
R134A
Manómetro de baja QG 1
S/4.70 S/4.70
R134A
Acople de 1/4 a 1/8 3 S/2.60 S/7.80
Tapón de ¼ 1 S/0.60 S/0.60
Unión de 5/8 soldable 2 S/0.50 S/1.00
Aislante para cañería 5/8 2 S/3.00 S/6.00
Mangueras para R134A 1 S/7.20 S/7.20
Tubo flexible de 3/8" 13 S/3.70 S/48.10
TOTAL S/299.40
Fuente: elaboración propia, 2019
[Link] Bienes:
Tabla 28. Costos de herramientas
VALOR
DESCRIPCIÓN CANTIDAD VALOR TOTAL
UNITARIO
Lima plana 1 S/6.00 S/6.00
Broca de 1/8 1 S/1.80 S/1.80
Broca de 5/32 1 S/3.00 S/3.00
Broca de Molibdeno HSS SAE 1
S/6.00 S/6.00
¼
Pistola para silicona 1 S/7.00 S/7.00
Prensa 1 S/0.00 S/0.00
Arco de Sierra 1 S/0.00 S/0.00
TOTAL S/23.8
Fuente: elaboración propia, 2019
VALOR VALOR
DESCRIPCIÓN CANTIDAD
UNITARIO TOTAL
TOTAL S/82.50
Fuente: elaboración propia, 2019
115
[Link] Personal:
[Link].1 Mano de obra directa
Es la fuerza laboral que se encuentra en contacto directo con la
fabricación del proyecto. Esta ha sido responsable del grueso
manejo del trabajo por horas y es el costo más importante de
controlar y medir. A continuación, se detalla la mano de obra
directa:
Tabla [Link] de mano de obra directa
VALOR/HOR VALOR
DESCRIPCIÓN HORAS
A TOTAL
Construcción de estructura 30 S/6.00 S/180.00
Servicio de doblez 1 S/2.60 S/2.60
Servicio de pintura 10 S/6.00 S/60.00
TOTAL S/242.60
Fuente: elaboración propia, 2019
DESCRIPCIÓN VALOR
TOTAL
Martinez Pucuhuayla Yojan Andre S/ 600.00
TOTAL S/ 600.00
Fuente: elaboración propia, 2019
116
otros: los gastos de administración, organización, supervisión, construcción de
instalaciones generales necesarias para realizar conceptos de trabajo, el
transporte, imprevistos y, en su caso, prestaciones laborales. A continuación, se
detalla los costos indirectos:
Tabla [Link]áneos
DESCRIPCIÓN VALOR
TOTAL
Bibliografía S/10.00
Internet S/40.00
Impresiones S/80.00
TOTAL S/ 235.00
Fuente: elaboración propia, 2019
DESCRIPCIÓN VALOR
TOTAL
Asesoramiento S/ 850,00
Misceláneos S/ 235,00
SUBTOTAL S/ 1085,00
Fuente: elaboración propia, 2019
117
Tabla 36. Costo total
TOTAL S/ 2569,10
Fuente: elaboración propia, 2019
118
CAPÍTULO VI
CONSTRUCCIÓN
6.1 Proceso de Construcción
P-10 Corte de tubería de cobre para las conexiones mecánicas entre los
elementos.
121
Marca Descripción
A-1 Compra y almacenamiento de materiales y materia prima en el
taller de trabajo.
P-1 Operación de trazar las medidas según diseño, en la materia prima.
P-2 Cortar el material para la estructura principal del recuperador.
I-1 Verificación de medidas finales en pieza para la estructura principal.
P-3 Soldar y ensamblar las piezas antes cortadas para armar la
estructura.
P-10 Corte de tubería de cobre para las conexiones mecánicas entre los
elementos.
123
CONCLUSIONES
124
TRABAJOS FUTUROS
125
126
REFERENCIAS BIBLIOGRÁFICAS
127
capacidad, usando refrigerante R404A para proceso de carga, descarga,
recuperación de refrigerante del sistema y simulacion de fallas. . Sangolqui :
Escuela Politecnica del Ejercito, 2012.
13. Mery, Castro Clavijo Luz. Recuperación del gas refrigerante R-22
clorofluorocarbonado, a traves de la unidad recuperadora CPS TR700 y el modulo
CRZRM PRO-SET. Acacias, Colombia : Universidad Nacional Abierta y a
Distancia, 2017.
14. Vizhñay, Jose Fernando Muñoz. Estudio del uso, manejo y reciclaje de
refrigerante automotriz . Cuenca,Ecuador : Universidad de Azuay, 2013.
15. José, Alarcón Creus. Teoria de la refrigeración . España : Marcombo S.A,
1999. 9788426716224.
16. Cengel, Yunus A. Transferencia de calor y masa. Mexico : McGraw-Hill
Interamericana, 2007. 978-970-10-6173-2.
17. —. Termodinamica. Mexico : Mc Graw-Hill, 1996. 978-0-07-352932-5.
18. 1MI131, Grupo. [Link]. [En línea] 21 de Noviembre de 2013. [Citado el:
20 de Julio de 2019.] [Link]
19. Manuel, Franco Lijó Juan. Manual de refrigeración . España : Reverté, 2006.
978-84-291-8011-4.
20. Benito, Torres. SCRIB. [En línea] Danfoss, 13 de Agosto de 2011. [Citado el:
28 de Julio de 2019.] [Link]
21. ConFrio. [Link]. [En línea] 11 de Junio de 2013. [Citado el: 19 de Mayo
de 2019.] [Link]
22. Cristina, Escudero Salas. Maquinas y equipos termicos, instalacion y
mantenimiento, instalacions frigorificas de climatizacion, instalaciones de
produccion de calor. Madrid : Paraninfo, 2013. 9788497322614.
23. DocPlayer. [Link]. [En línea] 23 de Setiembre de 2019. [Citado el: 1 de
Julio de 2019.] [Link]
[Link].
24. Ozono, Unidad Tecnica. Acoplasticos. [En línea] 2 de Febrero de 2019.
[Citado el: 29 de Junio de 2019.]
[Link]
%20de%20buenas%20pr%C3%A1cticas%20en%20refrigeraci%C3%B3n.
%20Unidad%20t%C3%A9cnica%20Ozono.%20Minambiente.%[Link].
25. Naturales(SEMARNAT), Secretaria del Medio Ambiente y Recursos.
Buenas prácticas en sistemas de refrigeracion y aire acondicionado. SEMARNAT.
[En línea] 2 de Enero de 2006. [Citado el: 20 de Mayo de 2019.]
[Link]
.pdf.
128
26. Gas-Servei. Gasservei. [En línea] 21 de Setiembre de 2019. [Citado el: 2 de
Junio de 2019.] [Link]
hfc/gasficha/r-134a/.
27. Celcius, Cero Grados. Cero Grados Celcius. [En línea] 17 de Enero de 2018.
[Citado el: 22 de Mayo de 2019.] [Link]
28. Rosario, Leopoldo. [Link]. [En línea] 28 de Noviembre de 2014.
[Citado el: 20 de Mayo de 2019.] [Link]
q=fugas+en+refrigerantes.
29. MOTOREX. MOTOREX. [En línea] 10 de Diciembre de 2018. [Citado el: 10 de
Junio de 2019.] [Link]
30. Cetis1466a3. Blogger. [En línea] 1 de Junio de 2016. [Citado el: 14 de Mayo
de 2019.] [Link]
31. TECUMSEH. Guia de refrigeracion . Brazil : FicFrio, 2001.
32. Samir, Saydaoui. Practicas de las maquinas frigorificas. México : Alfaomega
GrupoEditor, SA. de C.V., 2013. 978-607-707-156-3.
33. Antonio, Lopez Gomez. Las instalaciones frigorificas en las industrias
agroalimentarias. Madrid : A. Madrid Vicente, 1994. 84-87440-51-7.
34. Benjamín, Barriga Gamarra. Métodos de Diseño en Ingeniería Mecánica.
[Pdf] Lima : SCRIB, 2016.
35. [Link]. Quimobasicos S.A.C de C.V. [En línea] 18 de
Noviembre de 2018. [Citado el: 29 de Abril de 2019.]
[Link]
36. Santos, Noé Ramírez. MUNDO HVAC&R. [En línea] 02 de Febrero de 2013.
[Citado el: 24 de Julio de 2019.] [Link]
37. Montes, Ciro Espinosa. Metodología de Investigación Tecnológica.
Huancayo : Soluciones Gráficas S.A.C., 2014. 978-612-00-1667-1.
38. Diego, González Ruiz Jaime. Repository UPB. [En línea] 2 de Abril de 2012.
[Citado el: 20 de Mayo de 2019.]
[Link]
%[Link]?sequence=1.
39. Luis, Paredes Yataco Ismael. Diseño de un sistema de refrigeración con dos
temperatura de evaporación empleando un solo compresor. Lima : Pontificia
Universidad Católica del Perú, 2015. 20090288.
129
ANEXOS:
ANEXO 1
MATRIZ DE OPERACIONALIZACIÓN DE VARIABLES
DEFINICIÓN Variable que expresa la compresión de Variable que expresa la cantidad de refrigerante
OPERACIONAL vapor para un fluido frigorífero. recuperado.
125
VARIABLE RECUPERADORA DE REFRIGERANTE
INDEPENDIENTE
DIMENSIONES INDICADORES ACTIVIDADES DE INVESTIGACIÓN O TRATAMIENTO
Proceso de aplicación de la Temperatura °C/F
técnica Push Pull con la
finalidad de recuperar
refrigerante R134a. Presión Psi/bar
126
ANEXO 2
MATRIZ DE COHERENCIA
127
según la técnica coeficiente de PARA UN SISTEMA DE 6. Diseño conceptual
Succión/ Compresión? performance REFRIGERACIÓN POR 7. Recolección de
mayor a 3. COMPRESIÓN DE VAPOR antecedentes
DE 1HP DE CAPACIDAD, 8. Aplicación del
¿Cómo se realizaría la USANDO REFRIGERANTE
recuperación en fase R-404A PARA PROCESO Instrumento
liquida en la maquina DE CARGA, DESCARGA, (Operacionalización
recuperadora de RECUPERACIÓN DE de variables.)
refrigerante R134a? REFRIGERANTE DEL 9. Control de la calidad
SISTEMA Y SIMULACIÓN de los datos
OBJETIVO DE FALLAS PARA EL 10. Verificación de los
GENERAL LABORATORIO DE
requerimientos
Diseñar la máquina CONVERSIÓN DE
ENERGÍA DEL D.E.C.E.M. 11. Estudio del estado del
recuperadora de
refrigerante R-134a Tesis Descriptiva, Escuela arte
según la técnica Succión/ Politécnica del Ejercito, 12. Análisis Técnico
Compresión. Sangolquí, Ecuador, 2012. Económico
OBJETIVOS CASTRO, Luz. 13. Redacción del
ESPECÍFICOS RECUPERACION DEL GAS Proyecto de
Determinar cómo se REFRIGERANTE R-22
investigación
realizará el diseño para la CLOROFLUOROCARBONO,
A TRAVES DE LA UNIDAD 14. Elaboración del
recuperación en las
RECUPERADORA CPS Proyecto de
etapas
succión/compresión en la TR700 Y EL MODULO investigación
máquina recuperadora de CRXRM PRO-SET. Tesis 15. Revisión y
refrigerante R134a según Descriptiva, Universidad modificación del
la técnica Nacional Abierta y a Distancia Proyecto de
Succión/Compresión. UNAD, ACACIAS, Colombia,
investigación
2017.
Determinar el coeficiente TEORÍA BÁSICA
de Performance del Segunda ley de la Termodinámica.
Ciclo invertido de Carnot.
128
sistema.
DEFINICIÓN CONCEPTUAL
Conjunto de actividades y procesos
realizados en la máquina frigorífica,
aplicando la técnica
Succión/Compresión con la finalidad
de recuperar refrigerante r134a, que
condensa el refrigerante extraído en
fase liquida de cualquier máquina y lo
separa del aceite que pudiera llevar,
realizando así este trabajo en las fases
liquida y vapor.
DEFINICIÓN OPERACIONAL
Variable que expresa la compresión de
vapor para un fluido frigorífero.
129
ANEXO 3
PLANOS CONSTRUCTIVOS Y ESQUEMA DE FUNCIONAMIENTO DEL EQUIPO RECUPERADOR DE REFRIGERANTE
130
131
132
133
ANEXO 4
DIAGRAMA DE MOLLIERE PARA EL REFRIGERANTE R134
134
ANEXO 5
ESTRUCTURA DE PROGRAMACIÓN REALIZADA PARA EL SISTEMA DE
RECUPERACIÓN DE REFRIGERANTE
$WARNING OFF
"Función para determinar la capacitancia térmica mínima entre dos fluidos"
FUNCTION Cmin(C_refrigerante,C_secundario)
IF (C_refrigerante>C_secundario) THEN
Cmin:=C_secundario
ELSE
Cmin:=C_refrigerante
ENDIF
END
"Función para determinar la relación de capacitancias térmicas entre dos fluidos"
FUNCTION Cr(C_refrigerante,C_secundario)
IF (C_refrigerante>C_secundario) THEN
Cr:=C_secundario/C_refrigerante
ELSE
Cr:=C_refrigerante/C_secundario
ENDIF
END
"Procedimiento para direccionar las constantes adecuadas para la correlación de Zukauskas"
PROCEDURE Cm(Re:C,m)
IF (10<Re) and (Re<=100) THEN
C=0.9
m=0.4
ENDIF
IF (100<Re) and (Re<=1000) THEN
C=0.683
m=0.466
ENDIF
IF (1000<Re) and (Re<=200000) THEN
C=0.35
m=0.6
ENDIF
END
"Coeficientes para la eficiencia volumétrica"
a0=9.02017630E-01
a1=5.25041734E-05
a2=1.01302088E-03
a3=4.45451471E-04
a4=8.07560303E-07
a5=7.30516887E-02
a6=1.09549289E-03
a7=1.81500280E-05
a8=1.51540591E-03
a9=8.80285574E-05
"Coeficientes para la eficiencia combinada"
b0=1.07581877E+00
b1=7.34813961E-02
b2=7.93571514E-03
b3=7.68970614E-03
b4=1.11433756E-05
135
b5=1.05285567E+00
b6=1.65708127E-01
b7=1.96106652E-04
b8=9.64519165E-03
b9=3.74398417E-04
"DATOS GEOMÉTRICOS DEL CONDENSADOR"
A_ref_c=2.87 "área lado refrigerante parte externa tubos"
Din_c=13*convert(mm,m) "diámetro interno de los tubos en el condensador"
Dex_c=16*convert(mm,m) "diámetro externo de los tubos en el condensador"
Nt_c=10 "número de tubos en el condensador"
Dcor_c=183*convert(mm,m) "diámetro interno de la coraza"
Lc_c=19.5*convert(mm,m) "distancia entre centros"
Lt_c=3.5*convert(mm,m) "distancia entre tubos"
np_c=2 "número de pasos por coraza"
Ltubo_c=800*convert(mm,m) "longitud del tubo"
Nb_c=4 "número de baffles"
Nt_fila=4 "número de filas verticales de tubos"
K_tubo=400 ”conductividad del material del tubo"
Lcor_c=0.87 "longitud de la coraza"
Rt_i=0.00062 "resistencia interna asociada al ensuciamiento en tubos"
Dcor=0.195 “diámetro de la coraza”
"DATOS GEOMÉTRICOS DEL EVAPORADOR"
A_ref_e=1.81 "área lado refrigerante parte externa tubos"
Din_e=0.00822 "diámetro interno de los tubos en el evaporador"
Dex_e=0.00952 "diámetro externo de los tubos en el evaporador"
136
nu_suc=volume(refrigerant$,P=P_evap,T=T_7)
T_1=(T_7+273.15)*(rp)^((k-1)/k)-273.15
W_comp=(P_evap*nu_suc/(1-k))*((rp)^((1-k)/k)-1)*10^2 "trabajo teórico del compresor"
m_ref=(V_G*eta_v*N)/(60*nu_suc)
eta_comb=b0-b1*GR+b2*GR^2-b3*N+b4*N^2+b5*P_evap-b6*P_evap^2-
b7*GR*N+b8*GR*P_evapb9*
N*P_evap
Pot_comp*eta_comb=m_ref*W_comp
“LÍNEA DE DESCARGA”
"Balance energético para estimar la temperatura de superfice, Tw_l_des"
Tm_l_des=(T_1+T_2)/2 "temp. media del refrigerante"
Tw_l_des=(hi_l_des*Tm_l_des+ho_l_des*T_amb)/(hi_l_des+ho_l_des)
"Correlación Gnielinski, flujo interno lado refrigerante"
hi_l_des*Din_l_des/k_ref_des=(((ff_des/8)*(Re_ref_des-
1000)*Pr_ref_des)/(1+12.7*(ff_des/8)^0.5*(Pr_ref_des^(2/3)-1)))
ff_des=(0.790*ln(Re_ref_des)-1.64)^(-2)
k_ref_des=conductivity(refrigerant$,P=P_cond,T=Tm_l_des)
Cp_ref_des=CP(refrigerant$,P=P_cond,T=Tm_l_des)
rho_ref_des=density(refrigerant$,P=P_cond,T=Tm_l_des)
mu_ref_des=viscosity(refrigerant$,P=P_cond,T=Tm_l_des)
mu_w_des=viscosity(refrigerant$,P=P_cond,T=Tw_l_des)
Pr_ref_des=Cp_ref_des*mu_ref_des*1000/k_ref_des
Re_ref_des=4*m_ref/(pi*Din_l_des*mu_ref_des)
"Cilindro horizontal caso para el aire, multiplicada por un factor de conversión de unidades"
ho_l_des=(1.13*((Tw_l_des-T_amb)/Dex_l_des)^0.25)*(41.816/36) " análisis con
aproximaciones"
"Cilindro vertical caso para el aire"
ho_sep=(1.22*((Tw_l_des-T_amb)/D_sep)^0.25)*(41.816/36) " análisis con aproximaciones"
"Transferencia de calor al ambiente"
h_1=enthalpy(refrigerant$,P=P_cond,T=T_1)
h_2=enthalpy(refrigerant$,P=P_cond,T=T_2)
A_l_des=pi*Dex_l_des*L_l_des "superficie de la línea de descarga"
A_sep=pi*D_sep*L_sep+2*pi*D_sep^2/4 "superficie del separador de aceite"
Q_l_des=ho_l_des*A_l_des*(Tw_l_des-T_amb) "pérdida en línea de descarga"
Q_sep=ho_sep*A_sep*(Tw_l_des-T_amb) "pérdida en el separador de aceite"
Q_amb_l_des=Q_l_des+Q_sep "desde descarga hasta el condensador"
Q_amb_l_des=m_ref*(h_1-h_2)*1000
“CONDENSADOR”
P_cond=pressure(refrigerant$,T=T_cond,x=0)
rho_sc=density(water,T=(T_aec+T_asc)/2,P=1.9)
m_aec=C_aec*rho_sc*(1/3600)
"ZONA DE DESRECALENTAMIENTO"
"Lado refrigerante"
mu_v_drc=viscosity(refrigerant$,T=(T_cond+T_2)/2,P=P_cond)
mu_w_drc=viscosity(refrigerant$,T=Tw_drc,P=P_cond)
k_v_drc=conductivity(refrigerant$,T=(T_cond+T_2)/2,P=P_cond)
rho_v_drc=density(refrigerant$,T=(T_cond+T_2)/2,P=P_cond)
Cp_ref_drc=cp(refrigerant$,T=(T_2+T_cond)/2,P=P_cond)
h_2v=enthalpy(refrigerant$,P=P_cond,x=1)
h_2l=enthalpy(refrigerant$,P=P_cond,x=0)
Pr_v_drc=Cp_ref_drc*mu_v_drc*1000/k_v_drc
C_ref_drc=m_ref*Cp_ref_drc
"Lado agua"
mu_sc_drc=viscosity(water,T=(T_i+T_asc)/2,P=1.9)
rho_sc_drc=density(water,T=(T_i+T_asc)/2,P=1.9)
k_sc_drc=conductivity(water,T=(T_i+T_asc)/2,P=1.9)
Cp_sc_drc=CP(water,T=(T_i+T_asc)/2,P=1.9)
137
Pr_sc_drc=Cp_sc_drc*mu_sc_drc*1000/k_sc_drc
C_sc_drc=m_aec*Cp_sc_drc
"Balance de energía entre las corrientes y análisis térmico e-NTU"
m_ref*(h_2-h_2v)=m_aec*Cp_sc_drc*(T_asc-T_i)
Q_drc=m_ref*(h_2-h_2v)
m_ref*(h_2-h_2v)=epsilon_drc*Cmin_drc*(T_2-T_i)
Cmin_drc=Cmin(C_ref_drc,C_sc_drc)
Cr_drc=Cr(C_ref_drc,C_sc_drc)
epsilon_drc=2*(1+Cr_drc+(1+Cr_drc^2)^(1/2)*((1+exp(-NTU_drc*(1+Cr_drc^2)^(1/2)))/(1-
exp(-
NTU_drc*(1+Cr_drc^2)^(1/2)))))^(-1)
NTU_drc=U_drc*A_drc/(Cmin_drc*1000)
U_drc=1/((1/hi_drc)*(Dex_c/Din_c)+(Dex_c*ln(Dex_c/Din_c))/(2*K_tubo)+(1/ho_drc)
+Rt_i*(Dex_c/D
in_c))
"Correlación de COLBURN para flujo interno"
{hi_drc=0.023*(Re_sc_drc)^(4/5)*Pr_sc_drc^(1/3)*(K_sc_drc/Din_c)}
"Correlación DITTUS-BOELTER, flujo interno"
{hi_drc=0.023*Re_sc_drc^0.8*Pr_sc_drc^0.4*(K_sc_drc/Din_c)}
"Correlación de GNIELINSKI flujo interno"
hi_drc=((f_drc/8)*(Re_sc_drc-1000)*Pr_sc_drc)/(1+12.7*(f_drc/8)^(1/2)*(Pr_sc_drc^(2/3)-
1))*(K_sc_drc/Din_c)
A_paso_drc=(pi*Din_c^2/4)*(2*Nt_c/Np_c)
V_drc=m_aec/(rho_sc_drc*A_paso_drc)
Re_sc_drc=V_drc*rho_sc_drc*Din_c/mu_sc_drc
f_drc=1/(0.79*ln(Re_sc_drc)-1.64)^2
"Correlación de KERN flujo externo considerando un arreglo escalonado"
{ho_drc=(k_v_drc/De_drc)*0.36*(De_drc*Gs_drc/mu_v_drc)^(0.55)*(Cp_ref_drc*1000*mu_v
_drc/k_v
_drc)^(1/3)*(mu_v_drc/mu_w_drc)^0.14
Gs_drc=m_ref/As_drc
As_drc=(1-Dex_c/Lc_c)*Dcor_c*(Lcor_c/(Nb_c+1))
De_drc=4*(((Lc_c^2)*sqrt(3))/4-((pi*Dex_c^2))/8)/(pi*Dex_c/2)
Re_drc=De_drc*Gs_drc/mu_v_drc}
"Correlación de ZUKAUSKAS flujo externo considerando arreglo escalonado"
ho_drc=C_drc*(Re_drc)^(m_drc)*Pr_v_drc^(0.36)*(K_v_drc/Dex_c)*(mu_v_drc/mu_w_drc)^0
.25
Apaso=(Lcor_c/(Nb_c+1))*(Dcor_c-(Dcor_c/Lc_c-1)*Dex_c)
V_ref_drc=m_ref/(rho_v_drc*Apaso)
Re_drc=V_ref_drc*rho_v_drc*Dex_c/mu_v_drc
CALL Cm(Re_drc:C_drc,m_drc)
"Balance térmico auxiliar en la determinación de la temperatura de pared"
(Tw_drc-(T_asc+T_i)/2)*hi_drc=((T_2+T_cond)/2-Tw_drc)*ho_drc
"ZONA DE CONDENSACIÓN"
"Lado agua"
T_media=(T_i+T_ii)/2
rho_sc_cond=density(water,T=T_media,P=1.9)
Cp_sc_cond=CP(water,T=T_media,P=1.9)
mu_sc_cond=viscosity(water,T=T_media,P=1.9)
k_sc_cond=conductivity(water,T=T_media,P=1.9)
Pr_sc_cond=Cp_sc_cond*mu_sc_cond*1000/k_sc_cond
Cmin_cond=m_aec*Cp_sc_cond
"Lado refrigerante"
mu_l_cond=viscosity(refrigerant$,P=P_cond,x=0)
rho_l_cond=density(refrigerant$,P=P_cond,x=0)
k_l_cond=conductivity(refrigerant$,P=P_cond,x=0)
Cp_l_cond=CP(refrigerant$,P=P_cond,x=0)
138
rho_v_cond=density(refrigerant$,P=P_cond*0.99,x=1)
hfg_cond=(h_2v-h_2l)
"Balance de energía"
m_ref*(h_2v-h_2l)=m_aec*Cp_sc_cond*(T_i-T_ii)
Q_cond=m_aec*Cp_sc_cond*(T_i-T_ii)
m_ref*(h_2v-h_2l)=epsilon_cond*Cmin_cond*(T_cond-T_ii)
epsilon_cond=1-exp(-NTU_cond)
NTU_cond=U_c_cond*A_cond/(Cmin_cond*1000)
U_c_cond=1/((1/hi_cond)*(Dex_c/Din_c)+(Dex_c*ln(Dex_c/Din_c))/(2*K_tubo)+(1/ho_cond)
+Rt_i*(D
ex_c/Din_c))
"Correlación de NUSSELT flujo externo condensación sobre paquete de tubos vertical"
ho_cond*(variable1^0.25)=0.729*(variable2^0.25)
variable1=(Nt_fila)*mu_l_cond*Dex_c*(T_cond-Tw_c)
variable2=(rho_l_cond*(rho_l_cond-rho_v_cond)*g#*k_l_cond^3*hfg_mod*1000)
hfg_mod=hfg_cond*(1+0.68*Ja_cond)
Ja_cond=Cp_l_cond*(T_cond-Tw_c)/hfg_cond
"Correlación de COLBURN flujo interno"
{hi_cond=0.023*(Re_sc_cond)^(4/5)*Pr_sc_cond^(1/3)*(K_sc_cond/Din_c)}
"Correlación DITTUS-BOELTER flujo interno"
{hi_cond=0.023*Re_sc_cond^0.8*Pr_sc_cond^0.4*(K_sc_cond/Din_c)}
"Correlación de GNIELINSKI flujo interno"
hi_cond=((f_cond/8)*(Re_sc_cond-1000)*Pr_sc_cond)/
(1+12.7*(f_cond/8)^(1/2)*(Pr_sc_cond^(2/3)-
1))*(K_sc_cond/Din_c)
V_cond=m_aec/(rho_sc_cond*A_paso_drc)
Re_sc_cond=V_cond*rho_sc_cond*Din_c/mu_sc_cond
f_cond=1/(0.79*ln(Re_sc_cond)-1.64)^2
"Balance térmico auxiliar en la determinación de la temperatura de pared"
ho_cond*(T_cond-Tw_c)=hi_cond*(Tw_c-T_media) "aproximación"
"ZONA DE SUBENFRIAMIENTO"
"Lado refrigerante"
mu_v_sub=viscosity(refrigerant$,T=(T_cond+T_3)/2,P=P_cond)
mu_w_sub=viscosity(refrigerant$,T=Tw_sub,P=P_cond)
k_v_sub=conductivity(refrigerant$,T=(T_cond+T_3)/2,P=P_cond)
rho_v_sub=density(refrigerant$,T=(T_cond+T_3)/2,P=P_cond)
Cp_ref_sub=cp(refrigerant$,T=(T_3+T_cond)/2,P=P_cond)
Pr_v_sub=Cp_ref_sub*mu_v_sub*1000/k_v_sub
C_ref_sub=m_ref*Cp_ref_sub
"Lado agua"
mu_sc_sub=viscosity(water,T=(T_ii+T_aec)/2,P=1.9)
rho_sc_sub=density(water,T=(T_ii+T_aec)/2,P=1.9)
k_sc_sub=conductivity(water,T=(T_ii+T_aec)/2,P=1.9)
Cp_sc_sub=CP(water,T=(T_ii+T_aec)/2,P=1.9)
Pr_sc_sub=Cp_sc_sub*mu_sc_sub*1000/k_sc_sub
C_sc_sub=m_aec*Cp_sc_sub
"Balance de energía entre las corrientes y análisis térmico e-NTU"
Q_sub=m_aec*Cp_sc_sub*(T_ii-T_aec)
h_3=((m_ref*h_2l-epsilon_sub*Cmin_sub*(T_cond-T_aec))/m_ref)
T_3=temperature(refrigerant$,P=P_cond,h=h_3)
T_ii=m_ref*(h_2l-h_3)/(m_aec*Cp_sc_sub)+T_aec
Cmin_sub=Cmin(C_ref_sub,C_sc_sub)
Cr_sub=Cr(C_ref_sub,C_sc_sub)
epsilon_sub=2*(1+Cr_sub+(1+Cr_sub^2)^(1/2)*((1+exp(-NTU_sub*(1+Cr_sub^2)^(1/2)))/(1-
exp(-
NTU_sub*(1+Cr_sub^2)^(1/2)))))^(-1)
NTU_sub=U_sub*A_sub/(Cmin_sub*1000)
139
U_sub=1/((1/hi_sub)*(Dex_c/Din_c)+(Dex_c*ln(Dex_c/Din_c))/(2*K_tubo)+(1/ho_sub)
+Rt_i*(Dex_c/
Din_c))
"Correlación de GNIELINSKI flujo interno"
hi_sub=((f_sub/8)*(Re_sc_sub-1000)*Pr_sc_sub)/(1+12.7*(f_sub/8)^(1/2)*(Pr_sc_sub^(2/3)-
V_sub=m_aec/(rho_sc_sub*A_paso_drc)
Re_sc_sub=V_sub*rho_sc_sub*Din_c/mu_sc_sub
f_sub=1/(0.79*ln(Re_sc_sub)-1.64)^2
"Correlación de COLBURN flujo interno"
{hi_sub=0.023*(Re_sc_sub)^(4/5)*Pr_sc_sub^(1/3)*(K_sc_sub/Din_c)}
"Correlación DITTUS-BOELTER flujo interno"
{hi_sub=0.023*Re_sc_sub^0.8*Pr_sc_sub^0.4*(K_sc_sub/Din_c)}
"Correlación de KERN flujo externo considerando arreglo escalonado"
{ho_sub/
(k_v_sub/De_drc)=0.36*(De_drc*Gs_drc/mu_v_sub)^(0.55)*(Cp_ref_sub*1000*mu_v_sub/k_v
_sub)^(1/3)*(mu_v_sub/mu_w_sub)^0.14}
"Correlación de ZUKAUSKAS flujo externo considerando arreglo escalonado"
ho_sub*Dex_c/K_v_sub=C_sub*(Re_sub)^(m_sub)*Pr_v_sub^(0.36)*(mu_v_sub/mu_w_sub)^
0.25
V_ref_sub=m_ref/(rho_v_sub*Apaso)
Re_sub=V_ref_sub*rho_v_sub*Dex_c/mu_v_sub
CALL Cm(Re_sub:C_sub,m_sub)
"Balance térmico auxiliar en la determinación de la temperatura de pared"
(Tw_sub-(T_aec+T_ii)/2)*hi_sub=((T_3+T_cond)/2-Tw_sub)*ho_sub
A_ref_c=A_cond+A_drc+A_sub "ecuación de cierre"
Q_dot_cond=Q_drc+Q_cond+Q_sub "transferencia total en el condensador"
GS=7.51825445E-02+2.89126527E-01*P_cond-4.81205835E-01*P_evap
T_3=T_cond-GS
“VÁLVULA DE EXPANSIÓN”
T_4=-2.08016600E+00+8.99186897E-01*T_3+1.19696088E-01*T_amb
T_7=-1.81611023E-01+1.65695519E-01*T_amb+8.33658691E-01*T_6
rho_4=density(refrigerant$,T=T_4,P=(P_cond*0.96))
h_4=enthalpy(refrigerant$,T=T_4,P=(P_cond*0.96))
ka=0.0000024431 + 4.85700000E-08*T_evap
DPP=((P_cond*0.96)-P_evap)
DP=DPP*convert(bar,Pa)
m_cat=ka*(sqrt(rho_4*DP))
m_ref*(GR_max-SS)=m_cat*(GR-SS)
GR_max=-0.75 + 1.75*SS
T_6=GR+T_evap
“EVAPORADOR”
P_evap=pressure(refrigerant$,T=T_evap,x=1)
h_evap=enthalpy(refrigerant$,P=P_evap,x=1)
CALL BRINEPROP('Density','PG',conc,T_aee:rho_se)
m_aee=C_aee*rho_se*(1/3600)
"ZONA DE EVAPORACIÓN"
"Lado refrigerante"
h_5=h_4
x_5=quality(refrigerant$,h=h_5,P=P_evap)
mu_v_evap=viscosity(refrigerant$,x=1,P=P_evap)
k_v_evap=conductivity(refrigerant$,x=1,P=P_evap)
rho_v_evap=density(refrigerant$,x=1,P=P_evap)
Cp_ref_evap=cp(refrigerant$,x=1,P=P_evap)
h_l_evap=enthalpy(refrigerant$,x=0,P=P_evap)
mu_l_evap=viscosity(refrigerant$,x=0,P=P_evap)
k_l_evap=conductivity(refrigerant$,x=0,P=P_evap)
rho_l_evap=density(refrigerant$,x=0,P=P_evap)
140
Cp_l_evap=cp(refrigerant$,x=0,P=P_evap)
Pr_l_evap=Cp_l_evap*mu_l_evap*1000/k_l_evap
"Lado secundario, propylene glycol, [15.2<Conc<57.0]"
CALL BRINEPROP('ThermalC','PG',conc,Tw_evap:k_wbrine_evap)
CALL BRINEPROP('SpecHeat','PG',conc,Tw_evap:Cp_wbrine_evap)
CALL BRINEPROP('DynVisc','PG',conc,Tw_evap:vis_wbrine_evap)
mu_wbrine_evap=vis_wbrine_evap/1000
Pr_wbrine_evap=Cp_wbrine_evap*mu_wbrine_evap*1000/k_wbrine_evap
CALL BRINEPROP('Density','PG',conc,(T_o+T_ase)/2:rho_brine_evap)
CALL BRINEPROP('DynVisc','PG',conc,(T_o+T_ase)/2:vis_brine_evap)
CALL BRINEPROP('ThermalC','PG',conc,(T_o+T_ase)/2:k_brine_evap)
CALL BRINEPROP('SpecHeat','PG',conc,(T_o+T_ase)/2:Cp_brine_evap)
mu_brine_evap=vis_brine_evap/1000
Cmin_brine_evap=m_aee*Cp_brine_evap
Pr_brine_evap=Cp_brine_evap*mu_brine_evap*1000/k_brine_evap
"Balance de energía entre las corrientes y análisis térmico"
m_ref*(h_evap-h_5)=m_aee*Cp_brine_evap*(T_o-T_ase)
Q_evap=m_aee*Cp_brine_evap*(T_o-T_ase)
m_ref*(h_evap-h_5)=epsilon_evap*Cmin_brine_evap*(T_o-T_evap)
epsilon_evap=1-exp(-NTU_evap)
NTU_evap=U_evap*A_evap/(Cmin_brine_evap*1000)
U_evap=1/((1/h_e*(Dex_e/Din_e))+(Dex_e*ln(Dex_e/Din_e))/(2*K_tubo)+(1/ho_evap)+Rt_o)
"Correlación de CHEN, componente de ebullición nucleada + componente convectivo"
h_e=h_mic*S+h_mac*F
"Correlación DITTUS-BOELTER término convectivo"
h_mac=0.023*Re_l_evap^0.8*Pr_l_evap^0.4*(k_l_evap/Din_e)
"Correlación FORSTER-ZUBER término ebullición nucleada"
h_mic=0.00122*((k_l_evap^(0.79)*Cp_l_evap^(0.45)*rho_l_evap^(0.49))/
(sigma_e^(0.5)*mu_l_evap^(
0.29)*(landa_evap)^(0.24)*rho_v_evap^(0.24)))*factor
factor=deltaT_sat^0.24*deltaP_sat^0.75*(1000^0.45/1000^0.24)
sigma_e=SURFACETENSION(refrigerant$,T=T_evap)
landa_evap=h_evap-h_l_evap
Gs_ref=m_ref/As_ev
As_ev=(pi*Din_e^2/4)*(2*Nt_e/Np_e) "superficie total"
Re_l_evap=Gs_ref*(1-x_5)*Din_e/mu_l_evap
Re_bl=Re_l_evap*F^1.25 "número de Reynolds bifásico"
F=2.35*(1/(X_tt)+0.213)^0.736 "factor de mejoramiento"
X_tt=((1-x_5)/x_5)^0.9*(rho_v_evap/rho_l_evap)^0.5*(mu_l_evap/mu_v_evap)^0.1
S=1/(1+0.00000253*Re_bl^1.17) "factor de supresión"
T_brine=(T_ase+T_o)/2
T_sat=T_evap
P_sat=P_evap
P_w=pressure(refrigerant$,T=Tw_evap,x=1)
deltaT_sat=(Tw_evap-T_sat)
deltaP_sat=(P_w-P_sat)*convert(bar,Pa)
"Correlación de KERN, coeficiente de transferencia externo considerando arreglo escalonado"
{ho_evap=(k_brine_evap/De_evap)*0.36*(De_evap*Gs_evap/mu_brine_evap)^(0.55)*(Cp_brin
e_evap*
1000*mu_brine_evap/k_brine_evap)^(1/3)*factor1
factor1=(mu_brine_evap/mu_wbrine_evap)^0.14
Gs_evap=m_aee/As_evap
As_evap=(1-Dex_e/Lc_e)*Dcor_e*(Lcor/(Nbaffles+1))
De_evap=4*(((Lc_e^2)*sqrt(3))/4-((pi*Dex_e^2))/8)/(pi*Dex_e/2)}
"Correlación de ZUKAUSKAS, coeficiente de transferencia externo considerando arreglo
escalonado"
141
ho_evap=C_evap*(Re_brine_evap)^(m_evap)*Pr_brine_evap^(0.36)*(K_brine_evap/Dex_e)*(P
r_brine_
evap/Pr_wbrine_evap)^0.25
V=m_aee/(rho_brine_evap*A_paso)
A_paso=(Lcor/(Nbaffles+1))*(Dcor_e-(Dcor_e/Lc_e-1)*Dex_e)
Re_brine_evap=V*rho_brine_evap*Dex_e/mu_brine_evap
CALL Cm(Re_brine_evap:C_evap,m_evap)
"Balance térmico auxiliar en la determinación de la temperatura de pared"
ho_evap*(T_brine-Tw_evap)=h_e*(Tw_evap-T_evap) "aproximación"
"ZONA DE RECALENTAMIENTO"
"Lado refrigerante"
h_6=enthalpy(refrigerant$,T=T_6,P=P_evap)
mu_v_rec=viscosity(refrigerant$,T=(T_evap+T_6)/2,P=P_evap)
k_v_rec=conductivity(refrigerant$,T=(T_evap+T_6)/2,P=P_evap)
rho_v_rec=density(refrigerant$,T=(T_evap+T_6)/2,P=P_evap)
Cp_ref_rec=CP(refrigerant$,T=(T_evap+T_6)/2,P=P_evap)
Pr_ref_rec=Cp_ref_rec*mu_v_rec*1000/k_v_rec
C_ref_rec=m_ref*Cp_ref_rec
"Lado secundario, propylene glycol, [15.2<Conc<57.0]"
CALL BRINEPROP('Density','PG',conc,(T_aee+T_o)/2:rho_brine_rec)
CALL BRINEPROP('DynVisc','PG',conc,(T_aee+T_o)/2:vis_brine_rec)
CALL BRINEPROP('DynVisc','PG',conc,Tw_rec:vis_wbrine_rec)
mu_wbrine_rec=vis_wbrine_rec/1000
mu_brine_rec=vis_brine_rec/1000
CALL BRINEPROP('ThermalC','PG',conc,(T_aee+T_o)/2:k_brine_rec)
CALL BRINEPROP('SpecHeat','PG',conc,(T_aee+T_o)/2:Cp_brine_rec)
C_brine_rec=m_aee*Cp_brine_rec
Pr_brine_rec=Cp_brine_rec*mu_brine_rec*1000/k_brine_rec
T_brine2=(T_aee+T_o)/2
"Balance de energía entre las corrientes y análisis térmico"
m_ref*(h_6-h_evap)=m_aee*Cp_brine_rec*(T_aee-T_o)
Q_rec=m_ref*(h_6-h_evap)
m_ref*(h_6-h_evap)=epsilon_rec*Cmin_rec*(T_aee-T_evap)
Cmin_rec=Cmin(C_ref_rec,C_brine_rec)
Cr_rec=Cr(C_ref_rec,C_brine_rec)
epsilon_rec=2*(1+Cr_rec+(1+Cr_rec^2)^(1/2)*((1+exp(-NTU_rec*(1+Cr_rec^2)^(1/2)))/(1-
exp(-
NTU_rec*(1+Cr_rec^2)^(1/2)))))^(-1)
NTU_rec=U_rec*A_rec/(Cmin_rec*1000)
U_rec=1/((1/hi_rec*(Dex_e/Din_e))+(Dex_e*ln(Dex_e/Din_e))/(2*K_tubo)+(1/ho_rec)+Rt_o)
"Correlación de GNIELINSKI, flujo interno convección forzada"
hi_rec=((f_f/8)*(Re_ref_rec-1000)*Pr_ref_rec)/(1+12.7*(f_f/8)^(1/2)*(Pr_ref_rec^(2/3)-
1))*(K_v_rec/Din_e)
f_f=1/(0.79*ln(Re_ref_rec)-1.64)^2
V_ref_rec=m_ref/(rho_v_rec*As_ev)
Re_ref_rec=V_ref_rec*rho_v_rec*Din_e/mu_v_rec
"Correlación de COLBURN flujo interno"
{hi_rec=0.023*(Re_ref_rec)^(4/5)*Pr_ref_rec^(1/3)*(K_v_rec/Din_e)}
"Correlación DITTUS-BOELTER flujo interno"
{hi_rec=0.023*Re_ref_rec^0.8*Pr_ref_rec^0.4*(k_v_rec/Din_e)}
"Correlación de ZUKAUSKAS, coeficiente de transferencia externo considerando arreglo
escalonado"
ho_rec=C_rec*(Re_brine_rec)^(m_rec)*Pr_brine_rec^(0.36)*(K_brine_rec/Dex_e)*(mu_brine_
rec/mu_
wbrine_rec)^0.25
V_rec=m_aee/(rho_brine_rec*A_paso)
Re_brine_rec=V_rec*rho_brine_rec*Dex_e/mu_brine_rec
142
CALL Cm(Re_brine_rec:C_rec,m_rec)
"Balance térmico auxiliar en la determinación de la temperatura de pared"
ho_rec*((T_aee+T_o)/2-Tw_rec)=hi_rec*(Tw_rec-(T_evap+T_6)/2) "aproximación"
A_evap+A_rec=A_ref_e "ecuación de cierre"
Q_dot_evap=Q_evap+Q_rec "transferencia de calor total"
COP_G=Q_dot_evap/(Pot_comp+Pot_ba+Pot_bm)
"Caídas de presión a través de los intercambiadores de calor"
Dp_e=2261.95*C_aee^2-23.02
Dp_c=3143.49*C_aec^2-34.35
"Potencia consumida por la bomba del agua-glicol"
CALL BRINEPROP('DynVisc','PG',conc,T_aee:vis_brine)
mu_brine=vis_brine/1000
Re_brine=(4*rho_se*C_aee)/(pi*D_tubo*mu_brine)
f_tub_brine=0.25/(((log10((epsilon/(3.7*D_tubo))+(5.74/Re_brine^0.9)))))^2
Vel_brine=4*C_aee/(pi*D_tubo^2)
Dp_evap=f_tub_brine*(1/D_tubo)*(Vel_brine^2/(2*g#))*L_tubo+Dp_e
Pot_ba=((Dp_evap*C_aec/eta_ba)/3600)/1000
"Potencia consumida por la bomba del agua"
mu_agua=viscosity(water,T=T_aec,P=1.9)
rho_agua=density(water,T=T_aec,P=1.9)
Re_agua=(4*rho_agua*C_aec)/(pi*D_tubo*mu_agua)
f_tub_agua=0.25/(((log10((epsilon/(3.7*D_tubo))+(5.74/Re_agua^0.9)))))^2
Vel_agua=4*C_aec/(pi*D_tubo^2)
Dp_cond=f_tub_agua*(1/D_tubo)*(Vel_agua^2/(2*g#))*L_tubo+Dp_c
Pot_bm=((Dp_cond*C_aec/eta_bm)/3600)/1000
143