Mariscos Final

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 76

UNIVERSIDAD NACIONAL DE TRUJILLO

Facultad de Ingeniería
Escuela Profesional de Ingeniería Mecatrónica

Diseño de una cámara de refrigeración para el almacenamiento de


mariscos

Proyecto de Investigación Formativa

Transferencia de Calor y Refrigeración

AUTORES : ALVAREZ CARBAJAL, Raúl


ARANDA ULLOA, Christian Antonio
SAAVEDRA PAREDES, Luis renato

DOCENTE :
Ing. Bengoa Seminario Juan Carlos

CICLO :

VII

Trujillo, Perú
2020
UNIVERSIDAD NACIONAL DE TRUJILLO

Tabla de Contenidos

Capítulo 1 Introducción ......................................................................................................... 2


Objetivos............................................................................................................................. 3
Objetivo General............................................................................................................. 3
Objetivos específicos ...................................................................................................... 3
Capítulo II Marco Teórico ...................................................................................................... 4
Capítulo III Materiales y Métodos.......................................................................................... 7
Características del producto y oficina: ........................................................................... 7
Características de los materiales de construcción para las paredes: ............................... 7
Diseño de la cámara frigorifica: ..................................................................................... 8
Dimensiones de cámara: ............................................................................................... 10
Cálculo de resistencias térmicas en las paredes: .......................................................... 10
Calor por paredes y techo: ............................................................................................ 26
Carga térmica................................................................................................................ 29
Carga de la oficina: ....................................................................................................... 30
Resistencias térmicas de las paredes de la oficina: ...................................................... 31
Refrigerante: Amoniaco ............................................................................................... 35
Cálculo de tuberías: ...................................................................................................... 40
Selección de equipo: ..................................................................................................... 42
Capítulo IV Resultados y Discusiones ................................................................................. 52
Capítulo V Conclusiones ...................................................................................................... 53
Capítulo VI Recomendaciones ............................................................................................. 54
Referencias Bibliográficas .................................................................................................... 55
Anexo ................................................................................................................................... 56

2
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Lista de tablas

Tabla 1Mínimo espesor de aislamiento recomendado ........................................................... 8


Tabla 2 Dimensiones del carro refrigerante ........................................................................... 9
Tabla 3 Corrección de temperatura por radiación solar ....................................................... 11
Tabla 4 Conductividad térmicas de algunos materiales aislantes ........................................ 22
Tabla 5.Renovación de aire diario por las aberturas de las puertas ...................................... 27
Tabla 6. Calor liberado por persona ..................................................................................... 29
Tabla 7.Calor disipado por motores eléctricos ..................................................................... 29
Tabla 8. Número de renovaciones de aire ............................................................................ 31
Tabla 9 Velocidades recomendadas en las líneas ................................................................. 41
Tabla 10.Selección de compresores según su potencia de trabajo ....................................... 42
Tabla 11.Modelos de diferentes acumuladores de succión .................................................. 47
Tabla 12.Modelos de presostato marca Danfoss según el tipo de presión ........................... 50
Tabla 13.Tres tipos de termostatos marca Danfoss y características ................................... 51

3
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Lista de figuras

Figura 1 Esquematización referencial de largo, ancho y alto ................................................. 9


Figura 2.Dimensiones de la cámara de refrigeración ........................................................... 10
Figura 3 Diagrama presión-entalpía del ciclo ideal de refrigeración ................................... 37
Figura 4 Esquema del sistema de refrigeración de 2 evaporadores y 2 compresores .......... 38
Figura 5.Compresor Josval MISTRAL 40-A-AS ................................................................. 42
Figura 6.Compresor Josval MC-AI-1000/A ......................................................................... 43
Figura 7.Evaporador DeltaFrio, modelo Casset 232 ............................................................ 43
Figura 8.Evaporador DeltaFrio, modelo DFC73.74 ............................................................. 44
Figura 9.Condensador García Camara, modelo CG30 ......................................................... 44
Figura 10 Válvula de expansión Danfoss TEA 20-12 .......................................................... 45
Figura 11.Filtro de secado DCL 167s................................................................................... 45
Figura 12.Características del tanque recibidor CW-644 ...................................................... 46
Figura 13.Tanque recibidor CW-644 ................................................................................... 46
Figura 14. Acumulador de succión 2842-AS-5117-7 ........................................................... 48
Figura 15.Visor de refrigerante VCYL 33S ......................................................................... 48
Figura 16.Tubo anti vibrador EVCYDEAC 21 S ................................................................. 49
Figura 17 Presostato de aceite KP15 .................................................................................... 49
Figura 18 Termostatos EKC 202ª ......................................................................................... 51

4
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Lista de ecuaciones
Ecuación 1 Resistencia térmica para la convección ........................................................... 4
Ecuación 2 Resistencia térmica para la conducción ........................................................... 4
Ecuación 3 Numero de Reynolds ........................................................................................ 4
Ecuación 4 Longitud critica ................................................................................................ 4
Ecuación 5 Numero de Nusselt (placa plana) ..................................................................... 4
Ecuación 6 Numero de Nusselt para conveccion ................................................................ 4
Ecuación 7 Calor de las paredes ......................................................................................... 5
Ecuación 8 Calor por infiltración de aire ............................................................................ 5
Ecuación 9 Calor sobre congelación ................................................................................... 5
Ecuación 10 Calor emitido por persona .............................................................................. 5
Ecuación 11 Efecto refrigerante.......................................................................................... 5
Ecuación 12 Flujo masico 1 ................................................................................................ 5
Ecuación 13 Trabajo de los compresores............................................................................ 5
Ecuación 14 Coeficiente de Funcionamiento ..................................................................... 5
Ecuación 15 Potencia del compresor .................................................................................. 6
Ecuación 16 Flujo de calor del evaporador......................................................................... 6
Ecuación 17 Flujo de calor del condensador ...................................................................... 6
Ecuación 18 Diámetro de las tuberías ................................................................................. 6
Ecuación 19 Flujo de calor por motores ............................................................................. 6
Ecuación 20 Numero de Rayleigh ...................................................................................... 6

1
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Capítulo 1
Introducción

Los mariscos son, en gastronomía, un animal marino invertebrado comestible en general

son los moluscos y los crustáceos, es decir, los mejillones, las almejas, los berberechos,

chipirones y semejantes. Estos merecen un análisis aparte del pescado, pues poseen

nutrientes diferentes y valiosas propiedades para nuestro cuerpo.

Por tratarse de un ingrediente de origen animal, los mariscos ofrecen proteínas de alta

calidad, pues contienen todos los aminoácidos esenciales. Son bajos en hidrato de carbono

y la mayor parte de ellos tiene un bajo contenido graso, no alcanzando al 2% de su

composición.

Sin embargo, muchos de ellos tienen una elevada proporción de colesterol, tal es el caso

de los chipirones, los mejillones o los langostinos.

Todos tienen una alta densidad nutritiva, es decir, concentran muchos nutrientes saludables

en su composición y aunque son ingredientes bajos en grasa, contienen valiosos minerales

y vitaminas.

Entre las vitaminas destaca la vitamina E con función antioxidante y vitaminas del

complejo B, dentro de las cuales merece una mención especial el ácido fólico,

indispensable en la dieta de la embarazada y para prevenir anemias nutricionales. También

encontramos en la mayoría de los mariscos vitamina A, ideal en esta temporada para cuidar

la piel y la visión.

2
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Objetivos

Objetivo General

• Realizar el cálculo, diseño y especificaciones de un sistema de refrigeración


destinado a la conservación y almacenamiento de mariscos y con el mismo
refrigerante hacer un sistema de aire acondicionado

Objetivos específicos

• Dimensionar el cuarto de refrigeración.


• Calcular la carga térmica.
• Determinar la temperatura del condensador y evaporador
• Determinar los estados del ciclo de refrigeración.
• Hallar el efecto refrigerante del sistema
• Calcular el flujo másico de refrigerante en el ciclo
• Seleccionar el compresor, evaporador, condensador y válvula de expansión
que se adapte a las condiciones del ciclo.
• Dimensionar las tuberías para el sistema
• Seleccionar componentes secundarios para el sistema de refrigeración

3
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Capítulo II
Marco Teórico

• Resistencia térmica para la convección.


1
𝑅1 =
ℎ ∗ 𝐴𝑟𝑒𝑎
Ecuación 1 Resistencia térmica para la convección

• Resistencia térmica para la conducción:


𝑒
𝑅4 =
𝑘 ∗ 𝐴𝑟𝑒𝑎
Ecuación 2 Resistencia térmica para la conducción

• Numero de Reynolds:
𝑉𝐿
𝑅𝑒 =
𝑣
Ecuación 3 Numero de Reynolds

• Longitud crítica para un Reynolds de 5 * 105:


𝑅𝑒 ∗ 𝑣
𝐿𝑐𝑟𝑖𝑡 =
𝑉
Ecuación 4 Longitud critica

• Número de Nusselt para el flujo sobre una placa plana (ecuacion combinada)
ℎ𝐿 1
𝑁𝑢 = = (0.037𝑅𝑒𝑙0.8 − 871)𝑃𝑟 3
𝑘
Ecuación 5 Numero de Nusselt (placa plana)

• Numero de Nusselt para convección natural:


2

1
ℎ𝐿 0.387 ∗ 𝑅a𝐿 6
𝑁𝑢 = = 0.825 + 8 … (21)
𝑘 9 27
0.492 16
[1 + ( 𝑃𝑟 ) ]
{ }
Ecuación 6 Numero de Nusselt para convección

4
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

• Calor que pasa a través de las paredes:


∆𝑇
𝑄𝑝𝑎𝑟𝑒𝑑 =
𝑅𝑇
Ecuación 7 Calor de las paredes

• Calor que ingresa por infiltración de aire


𝑄𝐼𝑛𝑓𝑖𝑙. 𝐴𝑖𝑟𝑒 = 𝑉𝑜𝑙𝐶𝑎𝑚 × 𝑁°𝑅𝑒𝑛𝑜𝑣𝑎𝑐𝑖𝑜𝑛𝑒𝑠 × 𝐶𝑎𝑙𝑜𝑟𝐸𝑙𝑖𝑚𝑖𝑛𝑎𝑑𝑜 … (7)
Ecuación 8 Calor por infiltración de aire

• Calor de sobre congelación viene dado por:


𝑄𝑓 = 𝑊𝑃𝑟𝑜𝑑 ∗ 𝐶𝑓 ∗ ∆𝑇
Ecuación 9 Calor sobre congelación

• Calor liberado por persona:

̇ = 𝑁°𝑃𝑒𝑟𝑠𝑜𝑛𝑎𝑠 ∗ 𝐶 𝑝𝑒𝑟
𝑄𝑝𝑒𝑟𝑠
Ecuación 10 Calor emitido por persona

• La ecuacion para determinar el efecto refrigerante es la siguiente (ERU):

𝐸𝑅𝑈 = ℎ1 − ℎ8
Ecuación 11 Efecto refrigerante

• Para el flujo másico 1 tenemos:


̇
𝑄𝑡𝑜𝑡𝑎𝑙 𝑐𝑎𝑚𝑎𝑟𝑎
̇ =
𝑚1
𝐸𝑅𝑈
Ecuación 12 Flujo masico 1

• El trabajo de compresión de viene dado por:

𝑊 = ℎ2 − ℎ1 … (𝑏𝑎𝑗𝑎)
𝑊 = ℎ4 − ℎ3 … (𝑎𝑙𝑡𝑎)
Ecuación 13 Trabajo de los compresores

• El coeficiente de funcionamiento se determina mediante:

𝐸𝑅𝑈
𝐶𝑂𝑃 =
𝑊
Ecuación 14 Coeficiente de Funcionamiento

5
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

• La potencia del compresor viene dada por:

𝑃𝑜𝑡𝑐𝑜𝑚𝑝𝑟𝑒𝑠𝑜𝑟 = 𝑚̇ ×̇ (ℎ2 − ℎ1 )
Ecuación 15 Potencia del compresor

• El flujo de calor en el evaporador viene dado por:

̇
𝑄𝑒𝑣𝑎𝑝𝑜𝑟𝑎𝑑𝑜𝑟 = 𝑚(̇ℎ1 − ℎ4 )
Ecuación 16 Flujo de calor del evaporador

• El flujo de calor en el condensador viene dado por:

̇
𝑄𝑐𝑜𝑛𝑑𝑒𝑛𝑠𝑎𝑑𝑜𝑟 = 𝑚(̇ℎ2 − ℎ3 )
Ecuación 17 Flujo de calor del condensador

• El diámetro en las tuberías viene dado por:

4 × 𝜈 × 𝑚̇
∅= √
𝜋 × 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑
Ecuación 18 Diámetro de las tuberías

• Para hallar el flujo de calor debido a los motores, nos regimos mediante la
siguiente ecuación:
𝑄̇𝑚𝑜𝑡𝑜𝑟𝑒𝑠 = # 𝑚𝑜𝑡𝑜𝑟𝑒𝑠 ∗ 𝐶𝑎𝑙𝑜𝑟𝑒𝑙𝑖𝑚.𝑚𝑜𝑡𝑜𝑟 ∗ 𝑃𝑜𝑡𝑚𝑜𝑡𝑜𝑟
Ecuación 19 Flujo de calor por motores

• Número de Rayleigh:
𝑔 ∗ 𝛽 ∗ ∆𝑇 ∗ 𝐿3
𝑅a𝐿 = ∗ 𝑃𝑟
𝑣2
Ecuación 20 Numero de Rayleigh

• Calor de respiración del marisco


𝐶𝑟𝑒𝑠 ∗ 𝑇𝑜𝑛 ∗ 4.184
̇ =
𝑄𝑟𝑒𝑠
24 ∗ 3600
Ecuación 21 Carlo de respiración

6
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Capítulo III
Materiales y Métodos

Características del producto y oficina:

• Temperatura inicial del producto: 20 ºC


• Temperatura exterior: 25 ºC
• Temperatura de almacenamiento: -25ºC
• Diferencia de Temperatura: 40ºC
• Temperatura de la oficina: 21 ºC
• Flujo total del producto: 1 toneladas = 1ton/h
𝐾𝐽
• Calor específico de la palta sobre el punto de congelación: 𝐶𝑒 = 3.3913 𝐾𝑔∗°𝐶 =
𝐾𝑐𝑎𝑙
0.81 (Véase anexo 1)
𝐾𝑔∗°𝐶

Características de los materiales de construcción para las paredes:

• Coeficiente de transmisión de calor para el material aislante para pared: poliestireno


𝐾𝑐𝑎𝑙∗𝑐𝑚 𝐽
expandido: 𝐾 = 3.10 ℎ∗𝑚2 ∗⁰𝐶 = 0.035 (Véase anexo 2).
𝑠∗𝑚∗⁰𝐶

• Coeficiente de transmisión de calor para el material aislante para techo:


𝐾𝑐𝑎𝑙∗𝑐𝑚 𝐽
corcho aglomerado: 𝐾 = 3.10 ℎ∗𝑚2 ∗⁰𝐶
= 0.035 𝑠∗𝑚∗⁰𝐶
(Véase anexo 2).

• Espesor de pared aislante: 180 mm = 0.18 m (Véase tabla 1)


𝐾𝑐𝑎𝑙∗𝑐𝑚
• Coeficiente de transmisión de ladrillo para construcción: 𝐾 = 62.00 =
ℎ∗𝑚2 ∗⁰𝐶
𝐽
0.72 (Véase anexo 2)
𝑠∗𝑚∗⁰𝐶
𝐽
• Coeficiente de transmisión de yeso para techo de construcción: 𝐾 = 0.149 𝑠∗𝑚∗⁰𝐶

(Véase anexo 2)
• Espesor de paredes, interior y exterior: 10 cm = 0.1 m

7
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Tabla 1Mínimo espesor de aislamiento recomendado

Fuente: Módulo de refrigeración industrial, Tecsup

Diseño de la cámara frigorifica:

Para esto se tomo en cuenta los siguientes puntos:


• Jaba de pesca
• Carro de refrigeracion
• Espacio para el personal.

Jaba de pesca
Medidas de la jaba contenedora:
• Largo: 60 cmLargo: 60 cm
• Ancho: 40 cmAncho: 40 cm
• Altura: 16 cmAltura: 16 cm

Por cada jaba, se contiene 120 Camarones, 30 de Cangrejos, 2 pulpos.

Carro de refrigeración
Dimensiones del carro de refrigeración
• Largo: 160 cm
• Ancho 120 cm
• Alto: 160 cm (10 pisos)

Capacidad:

Peso por unidad de los Mariscos:


• Peso de Camarones= 40 gr a kg = 0.04kg*120= 4.8kg
• Peso de Cangrejos= 150 gr a kg = 0.15kg*30= 4.5kg
• Peso de Pulpos= 1.50 kg *2= 3kg

Peso Total = 4.8kg + 4.5kg + 3kg = 12.3kg

8
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Peso total de mariscos por jaba = 12.3kg


Cantidad: 80 jabas/carro.
Peso total de carga: 12.3 x 8(jabas dentro del carro) x 10(niveles de la jaba) = 984 Kg

Figura 1 Esquematización referencial de largo, ancho y alto


Fuente: FAO documento técnico de pesca

Metros cuadrados que ocupa una jaba en el carro congelante:

60𝑥40 = 2400 𝑐𝑚2


En cada carro se colocará un total de 80 jabas, nuestro carro elegido soporta una carga
de 984 kg. Tenemos:
1𝑡𝑜𝑛 24ℎ 24ton
∗ =
ℎ 1𝑑𝑖𝑎 𝑑𝑖𝑎
La cámara frigorífica deberá tener una capacidad de 24 toneladas.

La cantidad de carros congelantes que necesitamos es:

24000𝑘𝑔
= 16 carros refrigerantes
1500𝑘𝑔

Tabla 2 Dimensiones del carro refrigerante

Medidas del carro refrigerante cm


Largo 160
Ancho 120
Altura 160
Fuente: Propia

Metros cuadrados por piso del pallet:


160𝑥120 = 19 200𝑐𝑚2 = 192 𝑚2

9
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Número de jabas por piso en el carro refrigerante:

Metros cuadrados por piso del pallet 19200


= = 𝟖 𝒋𝒂𝒃𝒂𝒔
Metros cuadrados que ocupa una jaba en el carro congelante 2400

Hacemos una distribución de 4 carros refrigerantes a lo largo y 4 a lo ancho, puesto


que el total es 16 carros refrigerantes donde cada uno contendrá 80 jabas 8 por piso con
10 pisos, tomando todo esto en consideración tenemos las siguientes medidas para la
cámara frigorífica.

Dimensiones de cámara:

Se consideraron como dimensiones dependiendo la cantidad de carros refrigerantes


por sus dimensiones: Largo = 6,4 m; Ancho = 4,8 m; Alto = 1,6 m

4,8

6,4
1,6

Figura 2.Dimensiones de la cámara de refrigeración


Fuente: Propia

Cálculo de resistencias térmicas en las paredes:

Debemos tener en cuenta una corrección en las temperaturas debido a la radiación


solar y a la ubicación de las paredes, descritas en la tabla 3.

10
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Tabla 3 Corrección de temperatura por radiación solar

Fuente: Módulo de refrigeración industrial, Tecsup

Para poder adaptarlo a nuestro diseño de la cámara de refrigeración se realizaron los


siguientes cambios Norte por Sur y Sur por Norte, ya que los datos en la tabla 3 son del
hemisferio Norte.

A lo largo de la pared Norte


Dimensiones: L=4.8m y h=1.6m
• Temperatura corregida:
T∞ = 25 + 2.2
T∞ = 27.2°C
• Temperatura promedio:
27.2 + 40
𝑇𝑝𝑟𝑜𝑚 =
2
𝑇𝑝𝑟𝑜𝑚 = 15.6 °C

• Viscosidad cinemática del aire a 15 ⁰C= 1.470*10-5 m2/s (anexo 11)


• Numero de Prandtl a 15 ⁰C: 0.7323 (anexo 11)
𝐽
• Conductividad térmica a 15 ⁰C: 0.02476 (anexo 11)
𝑠∗𝑚∗℃

• Coeficiente de expansión volumétrica:


1
𝛽=
293𝑘
11
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Hallando el h de la convección natural interna de la cámara, primero debemos hallar


el número de Rayleigh (Ecuación 20):
1
𝑔 ∗ 𝛽 ∗ ∆𝑇 ∗ 𝐿3 9.81 ∗ ∗ (27.2) ∗ 4.83
𝑅a𝐿 = ∗ 𝑃𝑟 = 293 ∗ 0.035 ∗ 0.7323
𝑣2 (1.470 ∗ 10−5 )2
𝑅a𝐿 = 9.75 ∗ 1012
Después hallamos el número de Nusselt natural (Ecuación 6):
2

1 1
ℎ𝐿 0.387 ∗ 𝑅𝐚𝑳 6 0.387 ∗ (9.75 ∗ 1012 )6
𝑁𝑢 = = 0.825 + 8 = 0.825 + 8
𝑘 9 27 9 27
0.492 16 0.492 16
[1 + ( 𝑃𝑟 ) ] [1 + (0.7323) ]
{ } { }
𝑘𝑁𝑢 0.02476 ∗ 48.35
ℎ= =
𝐿 4.8
𝐽
𝒉 = 𝟎. 𝟐𝟓
𝑠 ∗ 𝑚2 ∗ ℃

Finalmente hallando R1 utilizando la ecuación 1:


𝐽
h = 0.25 𝑠∗𝑚2 ∗℃

área = 4.8 x 1.6 = 7.68m2


1 1 ℃∗𝐬
𝑅1 = = = 0.52
ℎ ∗ 𝐴𝑟𝑒𝑎 0.25 ∗ 7.68 𝐉

Hallando R2 y R4 con ecuación 2 y:


e = 0.1 m
𝐽
k = 0.035 𝑠∗𝑚∗⁰𝐶

área = 7.68m2
𝑒 0.1 ℃∗s
𝑅4 = 𝑅2 = = = 37.2 ∗ 10−2
𝑘 ∗ 𝐴𝑟𝑒𝑎 0.035 ∗ 7.68 J

12
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Resistencia térmica para la capa de poliestireno expandido R3 con ecuación 2 y


e = 0.18 m
𝐽
k = 0.035 𝑠∗𝑚∗⁰𝐶

área = 7.68m2
℃∗s
𝑅3 = 0.669
J

Resistencia térmica por convección externa forzada al exterior de la cámara de


refrigeración R5.
Hallando número de Reynolds con ecuación 3 y:
V = 2 m/s
L =4.8 m
𝑣 = 1.516*10-5 m2/s
𝑉𝐿
𝑅𝑒 = = 633245.38
𝑣
En vista de que la longitud crítica es considerable frente a la longitud característica,
optaremos por usar la ecuación 5 con:
Re = 633245.38
Pr = 0.7323
ℎ𝐿 1
𝑁𝑢 = = (0.037𝑅𝑒𝑙0.8 − 871)𝑃𝑟 3 = 674.928
𝑘
Despejando en ecuación 5:
𝑘𝑁𝑢 0.02514 𝑥674.928
ℎ= =
𝐿 4.8
𝐽
ℎ = 3.53
𝑠 ∗ 𝑚2 ∗ ℃
Finalmente hallando R5 empleando ecuación 1 y:
𝐽
h = 3.52 𝑠∗𝑚2 ∗℃

área = 7.68 m2
1 ℃∗s
𝑅5 = = 36.99𝑥10−3
ℎ ∗ 𝐴𝑟𝑒𝑎 J

13
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Empleando la ecuación 7 para hallar el calor que pasa a través de la pared del norte:
℃∗s
∆𝑇 = 27.2; 𝑅𝑇 = 28𝑥10−3 J

∆𝑇 𝐽
𝑄̇𝑝𝑎𝑟𝑒𝑑,𝑁𝑜𝑟𝑡𝑒 = = 974.53
𝑅𝑇 𝑠

A lo largo de la pared Sur


Dimensiones: L=4.8m y h=1.6m
Con el coeficiente de convección natural previamente hallado:
• Temperatura corregida:
T∞ = 25℃
• Temperatura promedio:
25 + 40
𝑇𝑝𝑟𝑜𝑚 =
2
𝑇𝑝𝑟𝑜𝑚 = 32.5 ℃

• Viscosidad cinemática a 32.5 ⁰C= 1.632*10-5 m2/s (Interpolación de anexo 11)


• Numero de Prandtl a 34.15 ⁰C: 0.7275(Interpolación de anexo 11)
𝐽
• Conductividad térmica a 34.15⁰C: 0.2606𝑠∗𝑚∗℃ (Interpolación de anexo 11)

• Coeficiente de expansión volumétrica:


1
𝛽=
292𝑘

Hallando el h de la convección natural interna de la cámara, primero debemos hallar


el número de Rayleigh (Ecuación 20):
1
𝑔 ∗ 𝛽 ∗ ∆𝑇 ∗ 𝐿3 9.81 ∗ ∗ (25) ∗ 4.83
𝑅a𝐿 = ∗ 𝑃𝑟 = 292 ∗ 0.035 ∗ 0.7275
𝑣2 (1.632 ∗ 10−5 )2
𝑅a𝐿 = 724.9 ∗ 1020

14
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Después hallamos el número de Nusselt natural (Ecuación 6):


2

1 1
ℎ𝐿 0.387 ∗ 𝑅𝐚𝑳 6 0.387 ∗ (724.9 ∗ 1020 )6
𝑁𝑢 = = 0.825 + 8 = 0.825 + 8
𝑘 9 27 9 27
0.492 16 0.492 16
[1 + ( 𝑃𝑟 ) ] [1 + (0.7275) ]
{ } { }
𝑁𝑢 = 2100.78
𝑘𝑁𝑢 0.2606 ∗ 2100.78
ℎ= =
𝐿 4.8
𝐽
ℎ = 114.05
𝑠 ∗ 𝑚2 ∗ ℃

Finalmente hallando R1 utilizando la ecuación 1:


𝐽
h = 114.05 𝑠∗𝑚2 ∗℃

área = 7.68 m2
℃∗s
𝑅1 = 1.14x10−3
J
Hallando R2 y R4 con ecuación 2 y:
e = 0.1 m
𝐽
k = 0.035 𝑠∗𝑚∗⁰𝐶

área = 7.68m2
𝑒 ℃∗𝑠
𝑅4 = 𝑅2 = = 0.327
𝑘 ∗ 𝐴𝑟𝑒𝑎 𝐽
Resistencia térmica para la capa de poliuretano R3 con ecuación 2 y:
e = 0.18 m
𝐽
k = 0.035 𝑠∗𝑚∗⁰𝐶

área = 7.68m2
℃∗𝑠
𝑅3 = 0.67
𝐽

15
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Resistencia térmica por convección externa forzada al interior de la cámara de


refrigeración R5:
Hallando número de Reynolds con ecuación 3 y:
V = 2 m/s
L = 4.8 m
𝑣 = 1.655 * 10-5 m2/s
𝑉𝐿
𝑅𝑒 = = 580060.42
𝑣
Hallando número de Nusselt con ecuación 5 y:

Re = 580060.42 Pr = 0.7268
ℎ𝐿 1
𝑁𝑢 = = (0.037𝑅𝑒𝑙0.8 − 871)𝑃𝑟 3 = 574.53
𝑘
Despejando en la ecuación 5:
𝑘𝑁𝑢 0.02625 𝑥574.53
ℎ= =
𝐿 4.8
𝐽
ℎ = 3.142
𝑠 ∗ 𝑚2 ∗ ℃
Finalmente hallando R5 empleando ecuación 1 y:
𝐽
h =2.142 𝑠∗𝑚2 ∗℃

área = 7.68m2
℃∗𝑠
𝑅5 = 0.06
𝐽
Empleando la ecuación 7 para hallar el calor que pasa a través de la pared del sur:
∆𝑇 = 25
℃∗s
𝑅𝑇 = 36 𝑥10−3
J
∆𝑇 𝐽
𝑄̇𝑝𝑎𝑟𝑒𝑑,𝑆𝑢𝑟 = = 694.18
𝑅𝑇 𝑠

16
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

A lo largo de la pared Oeste


Dimesiones: a=6.4m y h=1.6m
• Temperatura corregida:
T∞ = 25 + 3.3
T∞ = 28.3°C
• Temperatura promedio:
28.3 + 40
𝑇𝑝𝑟𝑜𝑚 =
2
𝑇𝑝𝑟𝑜𝑚 = 34.15°C

• Viscosidad cinemática a 34.15 ⁰C= 1.655*10-5 m2/s (Interpolación de anexo 11)


• Numero de Prandtl a 34.15 ⁰C: 0.7268 (Interpolación de anexo 11)
𝐽
• Conductividad térmica a 34.15⁰C: 0.02625 𝑠∗𝑚∗℃ (Interpolación de anexo 11)

• Coeficiente de expansión volumétrica:


1
𝛽=
294𝑘

Hallando el h de la convección natural interna de la cámara, primero debemos hallar


el número de Rayleigh (Ecuación 20):
1
𝑔 ∗ 𝛽 ∗ ∆𝑇 ∗ 𝐿3 9.81 ∗ ∗ (28,3) ∗ 4.83
𝑅a𝐿 = ∗ 𝑃𝑟 = 294 ∗ 0.035 ∗ 0.7268
𝑣2 (1.655 ∗ 10−5 )2
𝑅a𝐿 = 797.16 ∗ 1020

Después hallamos el número de Nusselt natural (Ecuación 21):


𝑁𝑢 = 416.45
𝑘𝑁𝑢 0.02518 ∗ 416.45
ℎ= =
𝐿 3.4
𝐽
ℎ = 3.084
𝑠 ∗ 𝑚2 ∗ ℃

17
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Hallando R1 de la ecuación 1 con:


𝑊
ℎ = 3.084
𝑚2
∗ °𝐶
área = 10.24 m2:
℃∗𝑠
𝑅1 = 𝟑𝟏. 𝟔 ∗ 𝟏𝟎−𝟑
𝐽
Hallando R2 y R4 con ecuación 2 y:
e = 0.1 m
𝐽
k = 0.035 𝑠∗𝑚∗⁰𝐶

área = 10.24 m2
℃∗𝑠
𝑅4 = 𝑅2 = 𝟎. 𝟐𝟖
𝐽
Resistencia térmica para la capa de poliuretano R3 con ecuación 2 y:
e = 0.18 m
𝐽
k = 0.035 𝑠∗𝑚∗⁰𝐶

área = 10.24 m2
℃∗𝑠
𝑅3 = 𝟎. 𝟓𝟎𝟐
𝐽
Resistencia térmica por convección externa forzada al interior de la cámara de
refrigeración R5:
Hallando número de Reynolds con ecuación 3 y:
V = 2 m/s
L = a = 6.4m
𝑣 =1.655* 10-5 m2/s
𝑅𝑒 = 773413.9
Hallando número de Nusselt con ecuación 5 y:

Re = 773413.9 Pr = 0.7268
1
𝑁𝑢 = (0.037𝑅𝑒𝑙0.8 − 871)𝑃𝑟 3 = 925,86

18
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Despejando en la ecuación 5:
𝑘𝑁𝑢 0.035 𝑥925.86
ℎ= =
𝐿 6.4
𝐽
ℎ = 5.06
𝑠 ∗ 𝑚2 ∗ ℃

Finalmente hallando R5 empleando ecuación 1 y:


𝐽
h =5.06 𝑠∗𝑚2 ∗℃; área = 10.24m2
℃∗𝑠
𝑅5 = 19.3𝑥10−3
𝐽
Empleando la ecuación 7 para hallar el calor que pasa a través de la pared del oeste:
℃∗s
∆𝑇 = 28.3; 𝑅𝑇 = 10.35 𝑥10−3 J

𝐽
𝑄̇𝑝𝑎𝑟𝑒𝑑,𝑂𝑒𝑠𝑡𝑒 = 2632.71
𝑠

A lo largo de la pared Este


Dimensiones: a=6.4m y h=1.6m
• Temperatura corregida:
T∞ = 25 + 3.3
T∞ = 28.3°C
• Temperatura promedio:
28.3 + 40
𝑇𝑝𝑟𝑜𝑚 =
2
𝑇𝑝𝑟𝑜𝑚 = 34.15°C

• Viscosidad cinemática a 34.15 ⁰C= 1.655*10-5 m2/s (Interpolación de anexo 11)


• Numero de Prandtl a 34.15 ⁰C: 0.7268 (Interpolación de anexo 11)
𝐽
• Conductividad térmica a 34.15⁰C: 0.02625 (Interpolación de anexo 11)
𝑠∗𝑚∗℃

• Coeficiente de expansión volumétrica:


1
𝛽=
294𝑘
19
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Hallando el h de la convección natural interna de la cámara, primero debemos hallar


el número de Rayleigh (Ecuación 20):
1
𝑔 ∗ 𝛽 ∗ ∆𝑇 ∗ 𝐿3 9.81 ∗ ∗ (28,3) ∗ 4.83
𝑅a𝐿 = ∗ 𝑃𝑟 = 294 ∗ 0.035 ∗ 0.7268
𝑣2 (1.655 ∗ 10−5 )2
𝑅a𝐿 = 797.16 ∗ 1020

Después hallamos el número de Nusselt natural (Ecuación 6):


𝑁𝑢 = 416.45
𝑘𝑁𝑢 0.02518 ∗ 416.45
ℎ= =
𝐿 3.4
𝐽
ℎ = 3.084
𝑠 ∗ 𝑚2 ∗ ℃
Hallando R1 de la ecuación 1 con:
𝑊
ℎ = 3.084
𝑚2
∗ °𝐶
Área = 10.24 m2:
℃∗𝑠
𝑅1 = 𝟑𝟏. 𝟔 ∗ 𝟏𝟎−𝟑
𝐽
Hallando R2 y R4 con ecuación 2 y:
e = 0.1 m
𝐽
k = 0.035 𝑠∗𝑚∗⁰𝐶

área = 10.24 m2
℃∗𝑠
𝑅4 = 𝑅2 = 𝟎. 𝟐𝟖
𝐽
Resistencia térmica para la capa de poliuretano R3 con ecuación 2 y:
e = 0.18 m
𝐽
k = 0.035 𝑠∗𝑚∗⁰𝐶

área = 10.24 m2
℃∗𝑠
𝑅3 = 𝟎. 𝟓𝟎𝟐
𝐽

20
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Resistencia térmica por convección externa forzada al interior de la cámara de


refrigeración R5:
Hallando número de Reynolds con ecuación 3 y:
V = 2 m/s
L = a = 6.4m
𝑣 =1.655* 10-5 m2/s
𝑅𝑒 = 773413.9
Hallando número de Nusselt con ecuación 5 y:
Re = 773413.9
Pr = 0.7268
1
𝑁𝑢 = (0.037𝑅𝑒𝑙0.8 − 871)𝑃𝑟 3 = 925,86
Despejando en la ecuación 5:

𝑘𝑁𝑢 0.035 𝑥925.86


ℎ= =
𝐿 6.4
𝐽
ℎ = 5.06
𝑠 ∗ 𝑚2 ∗ ℃
Finalmente hallando R5 empleando ecuación 1 y:
𝐽
h =5.06 𝑠∗𝑚2 ∗℃

área = 10.24m2
℃∗𝑠
𝑅5 = 19.3𝑥10−3
𝐽

Resistencia térmica de la puerta:

La puerta es de 5cm y consta de dos capas de madera (k=0.274 W/m. K, véase tabla
4) de 1cm y una capa central de poliestireno expandido (k=0.035) de 3cm de espesor.
Hallando la resistencia térmica de la madera de la ecuación 2 con:
e = 0.01 m
k = 0.274 W/m. K
área = 2 m2

21
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

𝑅𝑚𝑎𝑑𝑒𝑟𝑎 = 18.25 ∗ 10−3


Hallando la resistencia térmica del poliestireno expandido de la ecuación 2 con:
e = 0.03 m
k = 0.0361 W/m. K
área = 2 m2
℃∗s
𝑅𝑝𝑜𝑙.𝑒𝑥𝑝 = 415.51 ∗ 10−3
J
Finalmente, la resistencia térmica total será:
𝑅𝑝 = 2 ∗ 𝑅𝑚𝑎𝑑𝑒𝑟𝑎 + 𝑅𝑝𝑜𝑙.𝑒𝑥𝑝
℃∗s
𝑅𝑝 = 452 ∗ 10−3
J

Tabla 4 Conductividad térmicas de algunos materiales aislantes

Fuente: Código técnico de la edificación

Empleando la ecuación 7 para hallar el calor que pasa a través de la pared del este:
∆𝑇 = 28.3
℃∗s
𝑅𝑇 = 88.28 𝑥10−3
J
Esta última resultó de tratar las dos resistencias térmicas de la pared en serie, la
resistencia térmica de la puerta en paralelo más las resistencias de la pared.
𝐽
𝑄̇𝑝𝑎𝑟𝑒𝑑,𝐸𝑠𝑡𝑒 = 320.57
𝑠

22
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

A lo largo del techo


Dimensiones L=4.8m y h=6.4m
• Temperatura corregida:
T∞ = 25 + 8.3
T∞ = 33.3°C
• Temperatura promedio:
33.3 + 40
𝑇𝑝𝑟𝑜𝑚 =
2
𝑇𝑝𝑟𝑜𝑚 = 36.65°C

• Viscosidad cinemática a 36.65⁰C= 1.67*10-5 m2/s (Interpolación de anexo 11)


• Numero de Prandtl a 36.65⁰C: 0.2763(Interpolación de anexo 11)
𝐽
• Conductividad térmica a 36.65⁰C: 0.0276𝑠∗𝑚∗℃ (Interpolación de anexo 11)

• Coeficiente de expansión volumétrica:


1
𝛽=
296𝑘

Hallando el h de la convección natural interna de la cámara, primero debemos hallar


el número de Rayleigh (Ecuación 20), antes de eso como la superficie es horizontal se
tiene las siguientes consideraciones:
𝐴𝑠 1.6 ∗ 6.4
𝐿𝑐 = = = 0.64𝑚
𝑝 2 ∗ 1.6 + 2 ∗ 6.4
1
𝑔 ∗ 𝛽 ∗ ∆𝑇 ∗ 𝐿3 9.81 ∗ ∗ (33,3) ∗ 4.83
𝑅a𝐿 = ∗ 𝑃𝑟 = 296 ∗ 0.035 ∗ 0.2763
𝑣2 (1.67 ∗ 10−5 )2
𝑅a𝐿 = 𝟑𝟒𝟓. 𝟒𝟖 ∗ 𝟏𝟎𝟐𝟎
2

1
ℎ𝐿 0.387 ∗ 𝑅𝐚𝑳 6
𝑁𝑢 = = 0.825 + 8 = 1712.9
𝑘 9 27
0.492 16
[1 + ( 𝑃𝑟 ) ]
{ }

23
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

𝑘𝑁𝑢 0.0276 𝑥1712.9


ℎ= =
𝐿 0.64
𝐽
ℎ = 73.86
𝑠 ∗ 𝑚2 ∗ ℃

Hallando R1 de la ecuación 1 con:


𝑊
ℎ = 73.86
𝑚2
∗ °𝐶
área = 30.72m2
℃∗s
𝑅1 = 44 ∗ 10−3
J
Hallando R2 y R4 con ecuación 2 y:
e = 0.1 m
𝐽
k = 0.149 𝑠∗𝑚∗⁰𝐶

área = 30.72 m2
℃∗s
𝑅4 = 𝑅2 = 21.84 ∗ 10−3
J
Resistencia térmica para la capa de poliestireno expandido R3 con ecuación 2 y:
e = 0.07 m
𝐽
k = 0.035 𝑠∗𝑚∗⁰𝐶

área = 30.72 m2
℃∗s
𝑅3 = 65.1 ∗ 10−3
J
Resistencia térmica por convección externa forzada al interior de la cámara de
refrigeración R5:
Hallando número de Reynolds con ecuación 3 y:
V = 2 m/s
L = 4.8 m
𝑣 = 1.67 ∗ 10−5 m2/s
𝑅𝑒 = 574850.29

24
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Hallando número de Nusselt con ecuación 5 y:


Re = 574850.29
Pr = 0.2763
ℎ𝐿 1
𝑁𝑢 = = (0.037𝑅𝑒𝑙0.8 − 871)𝑃𝑟 3 = 409.12
𝑘
Despejando en la ecuación 5:
𝑘𝑁𝑢 0.035𝑥409.12
ℎ= =
𝐿 4.8
𝐽
ℎ = 2.98
𝑠 ∗ 𝑚2 ∗ ℃
Finalmente hallando R5 empleando ecuación 1 y:
𝐽
h = 2.98 𝑠∗𝑚2 ∗℃

área =30.72 m2
℃∗s
𝑅5 = 11𝑥10−3
J
Empleando la ecuación 7 para hallar el calor que pasa a través del techo:
℃∗s
∆𝑇 = 33.3; 𝑅𝑇 = 4.53𝑥10−3 J

𝐽
𝑄̇𝑡𝑒𝑐ℎ𝑜 = 7345.06
𝑠

25
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Calor por paredes y techo:

Para hallar el flujo de calor hacia la cámara de refrigeración sumamos los flujos de
calor de las cuatro paredes y el techo, resultando:
𝐽
𝑄̇𝑐𝑎𝑚𝑎𝑟𝑎 = 11 967.04
𝑠

Para determinar el calor que ingresa por infiltración de aire necesitaremos saber el
volumen de la cámara de refrigeración, el número de renovaciones y el calor eliminado.
𝑉𝑜𝑙𝐶𝑎𝑚 = 𝐴𝑛𝑐ℎ𝑜𝐶𝑎𝑚 × 𝐿𝑎𝑟𝑔𝑜𝐶𝑎𝑚 × 𝐴𝑙𝑡𝑜𝐶𝑎𝑚

𝑉𝑜𝑙𝐶𝑎𝑚 = 4.8 × 6.4 × 1.6

𝑉𝑜𝑙𝐶𝑎𝑚 = 49.15 𝑚3

El número de renovaciones de aire diario es un valor tabulado (véase tabla 6):


𝑁°𝑟𝑒𝑛𝑜𝑣 = 9.85

De igual manera, el calor eliminado es un valor tabulado en función de la humedad


relativa y de la temperatura del aire de entrada (anexo 3); para nuestro caso una humedad
relativa de 90% y la temperatura del aire de entrada a 25°C, tenemos:
𝐾𝑐𝑎𝑙
𝐶𝑎𝑙𝑜𝑟𝐸𝑙𝑖𝑚𝑖𝑛𝑎𝑑𝑜 = 15,42 = 72.82 𝐾𝐽/𝑚3
𝑚3

Finalmente tenemos, según la ecuación 7, que el calor que ingresa por infiltración de
aire es:
𝑄𝐼𝑛𝑓𝑖𝑙. 𝐴𝑖𝑟𝑒 = 𝑉𝑜𝑙𝐶𝑎𝑚 × 𝑁°𝑅𝑒𝑛𝑜𝑣𝑎𝑐𝑖𝑜𝑛𝑒𝑠 × 𝐶𝑎𝑙𝑜𝑟𝐸𝑙𝑖𝑚𝑖𝑛𝑎𝑑𝑜

𝑄𝐼𝑛𝑓𝑖𝑙. 𝐴𝑖𝑟𝑒 = 35 254.16𝐾𝐽

Como se desea llegar a la temperatura de almacenamiento en una hora (3600 s):


𝑄𝐼𝑛𝑓𝑖𝑙.̇ 𝐴𝑖𝑟𝑒 = 35 254.16𝐾𝐽/(3600 𝑠𝑒𝑔)

𝑄𝐼𝑛𝑓𝑖𝑙.̇ 𝐴𝑖𝑟𝑒 = 9.81 𝐾𝐽/𝑠𝑒𝑔

26
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Tabla 5.Renovación de aire diario por las aberturas de las puertas

Fuente: Módulo de refrigeración industrial, Tecsup.

Para hallar el calor de sobre congelación, obtenemos el valor de Cf de anexo 1.


𝐾𝐽
𝐶𝑓 = 0.43 𝐾𝑔∗°𝐶.

Finalmente, tenemos que calor de sobre congelación, aplicando la ecuación 9:


Wprod = 1000 kg
𝐾𝐽
𝐶𝑓 = 0.43 𝐾𝑔∗°𝐶;

∆𝑇 = 40°𝐶
𝑄𝑓 = 𝑊𝑃𝑟𝑜𝑑 ∗ 𝐶𝑓 ∗ ∆𝑇
𝑄𝑓 = 17200 KJ
Pero, en el lapso de una hora:
17200KJ
𝑄𝑓̇ =
3600𝑠
𝐾𝐽
𝑄𝑓̇ = 4.77
𝑠

Para el calor de respiración del marisco se consideró valor de tabla. (Véase anexo 4).
Determinando el flujo de calor debido a la respiración del marisco, según la ecuación 21

27
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

𝑊
Toneladas = 1; Cres = 12 𝑡𝑜𝑛𝑒𝑙𝑎𝑑𝑎𝑠
𝐶𝑟𝑒𝑠 ∗ 𝑇𝑜𝑛 ∗ 4.184
̇ =
𝑄𝑟𝑒𝑠
24 ∗ 3600
𝐽
̇ = 0.00058
𝑄𝑟𝑒𝑠
𝑠

Calor liberado por personas:


El calor liberado por persona está en función de la temperatura de nuestra cámara de
refrigeración (-25°C) (véase tabla 7):
𝐾𝑐𝑎𝑙
𝐶𝑝𝑒𝑟 = 363

Hallando el flujo de calor debido al calor liberado por personas, empleamos la


ecuación 10 con:
N° personas = 2
Cper = 363kcal/h
̇ = 𝑁°𝑃𝑒𝑟𝑠𝑜𝑛𝑎𝑠 ∗ 𝐶 𝑝𝑒𝑟
𝑄𝑝𝑒𝑟𝑠
𝐽
̇ = 726
𝑄𝑝𝑒𝑟𝑠
𝑠

Calor suministrado por motores:


El calor liberado por motor está en función del tipo de cámara que se tenga (véase
tabla 8), para nuestro caso:
𝐾𝑐𝑎𝑙
𝐶𝑎𝑙𝑜𝑟𝐸𝑙𝑖𝑚.𝑀𝑜𝑡𝑜𝑟 = 640
𝐻𝑃 ∗ ℎ

Finalmente, hallamos el flujo de calor suministrado por los 3 motores, mediante la


ecuación 19, resultando:
𝑄̇𝑚𝑜𝑡𝑜𝑟𝑒𝑠 = # 𝑚𝑜𝑡𝑜𝑟𝑒𝑠 ∗ 𝐶𝑎𝑙𝑜𝑟𝑒𝑙𝑖𝑚.𝑚𝑜𝑡𝑜𝑟 ∗ 𝑃𝑜𝑡𝑚𝑜𝑡𝑜𝑟
𝐾𝐽
̇
𝑄𝑚𝑜𝑡𝑜𝑟𝑒𝑠 = 1.08
𝑠

28
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Tabla 6. Calor liberado por persona

Fuente: Módulo de refrigeración industrial, Tecsup.

Tabla 7.Calor disipado por motores eléctricos

Fuente: Módulo de refrigeración industrial, Tecsup.

Carga térmica

Carga térmica total:


̇ = 𝑄𝑐𝑎𝑚𝑎𝑟𝑎
𝑄𝑡𝑜𝑡𝑎𝑙 ̇ + 𝑄𝐼𝑛𝑓𝑖𝑙. ̇ + 𝑄̇ + 𝑄 ̇ + 𝑄 ̇ + 𝑄 ̇ ̇
𝐴𝑖𝑟𝑒 𝑓 𝑟𝑒𝑠 𝑝𝑒𝑟𝑠 𝑚𝑜𝑡𝑜𝑟𝑒𝑠 + 𝑄𝑙𝑢𝑐𝑒𝑠

̇ = 2.94 ∗ A ∗ w ∗ n
𝑄𝑙𝑢𝑐𝑒𝑠
̇ = 11 967.04 + 93.34 𝑥 103 +4.77 𝑥 103 + 0.00058 + 726 + 1.08 𝑥 103 + 30
𝑄𝑡𝑜𝑡𝑎𝑙
𝐽
̇ = 111 913.04
𝑄𝑡𝑜𝑡𝑎𝑙
𝑠

29
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Cárga térmica de diseño


El compresor trabaja 16 horas de las 24 del día.
24 / 16 = 1.5
𝐽
𝑄̇𝑐𝑎𝑟𝑔𝑎 𝑡é𝑟𝑚𝑖𝑐𝑎 = 1.5 𝑥 111 913.04 = 167 870.01
𝑠
𝐽
𝑄̇𝑐𝑎𝑟𝑔𝑎 𝑡é𝑟𝑚𝑖𝑐𝑎 𝑓𝑎𝑐𝑡𝑜𝑟 𝑑𝑒 𝑠𝑒𝑔𝑢𝑟𝑖𝑑𝑎𝑑 = 1.1 𝑥 167 870.01 = 184 656.51
𝑠

Carga de la oficina:

En este caso como no se nos dio como tarea diseñar la oficina se optó por despreciar el
calor que pasan a través de las paredes y solo considerar el calor que emiten las personas,
junto con el calor de la renovación del aire, así tenemos:

Calor por infiltración de aire:


Supongamos que la oficina tiene un volumen de 60𝑚3
𝑟𝑒𝑛𝑣 𝑟𝑒𝑛𝑣
𝑁°𝑟𝑒𝑛𝑜𝑣 = 5.6 = 134.4
ℎ 𝑑𝑖𝑎
𝐾𝑐𝑎𝑙
𝐶𝑎𝑙𝑜𝑟𝐸𝑙𝑖𝑚𝑖𝑛𝑎𝑑𝑜 = 8.37 3 = 41.67 𝐾𝐽/𝑚3
𝑚
𝑄𝐼𝑛𝑓𝑖𝑙. 𝐴𝑖𝑟𝑒 = 𝑉𝑜𝑙𝐶𝑎𝑚 × 𝑁°𝑅𝑒𝑛𝑜𝑣𝑎𝑐𝑖𝑜𝑛𝑒𝑠 × 𝐶𝑎𝑙𝑜𝑟𝐸𝑙𝑖𝑚𝑖𝑛𝑎𝑑𝑜
𝑄𝐼𝑛𝑓𝑖𝑙. 𝐴𝑖𝑟𝑒 = 336 026,82 𝑘𝐽

Como se desea llegar a la temperatura de confort en una hora (3600 s):


336 026,82
𝑄𝐼𝑛𝑓𝑖𝑙. 𝐴𝑖𝑟𝑒 = 𝐾𝐽
3600𝑠
𝑄𝐼𝑛𝑓𝑖𝑙. 𝐴𝑖𝑟𝑒 = 93.32𝐾𝐽/𝑆

Calor liberado por personas:


El calor liberado por persona está en función de la temperatura de la oficina (véase
tabla 6), para nuestro caso:
𝐾𝑐𝑎𝑙
𝐶𝑝𝑒𝑟 = 121

30
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Hallando el flujo de calor debido al calor liberado por personas, empleamos la ecuación
11 con:
N° personas = 10; Cper = 121 kcal/h
̇ = 𝑁°𝑃𝑒𝑟𝑠𝑜𝑛𝑎𝑠 ∗ 𝐶 𝑝𝑒𝑟
𝑄𝑝𝑒𝑟𝑠
𝐽
̇ = 1210
𝑄𝑝𝑒𝑟𝑠
𝑠

Tabla 8. Número de renovaciones de aire

Fuente: scalofrios.es

Resistencias térmicas de las paredes de la oficina:

Para estos cálculos solo consideramos transferencia de calor por conducción térmica y
que la oficina está hecha de ladrillo (k=0.72) y el techo de yeso (0.149).

Pared norte
Dimensiones L=4.8m y h=1.6m
• Temperatura corregida:
T∞ = 25 + 2.2
T∞ = 27.2°C
• Temperatura promedio:
27.2 + 21
𝑇𝑝𝑟𝑜𝑚 =
2
𝑇𝑝𝑟𝑜𝑚 = 24°C

31
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

• Viscosidad cinemática del aire a 24 ⁰C = 1.5528*10-5 m2/s (Interpolación de


anexo 11)
• Número de Prandtl a 24 ⁰C: 0.73 (Interpolación de anexo 11)
𝐽
• Conductividad térmica a 24 ⁰C: 0.0254 𝑠∗𝑚∗℃ (Interpolación de anexo 11)

Tenemos en cuenta que el espesor de las paredes es 0.1m y aplicando la ecuación 2:


e = 0.1 m
𝐽
k = 0.72 𝑠∗𝑚∗⁰𝐶

área = 7.68 m2
℃∗s
𝑅𝑐𝑜𝑛𝑑 = 4.76 ∗ 10−3
J
℃∗s
∆𝑇 = 6.2; 𝑅𝑇 = 9.26𝑥10−3 J

𝐽
𝑄̇𝑝𝑎𝑟𝑒𝑑,𝑁𝑜𝑟𝑡𝑒 = 669.54
𝑠

Pared sur
Dimensiones L=4.8m y h=1.6m
• Temperatura corregida:
T∞ = 25 + 0
T∞ = 25°C
• Temperatura promedio:
25 + 40
𝑇𝑝𝑟𝑜𝑚 =
2
𝑇𝑝𝑟𝑜𝑚 = 32.5°C

• Viscosidad cinemática del aire a 32.5 ⁰C= 1.628*10-5 m2/s (Interpolación de


anexo 11)
• Número de Prandtl a 32.5 ⁰C: 0.7281 (Interpolación de anexo 11)
𝐽
• Conductividad térmica a 32.5 ⁰C: 0.0261 (Interpolación de anexo 11)
𝑠∗𝑚∗℃

32
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Tenemos en cuenta que el espesor de las paredes es 0.1m y aplicando la ecuación 2:


e = 0.1 m
𝐽
k = 0.72 𝑠∗𝑚∗⁰𝐶

área = 7.68m2
℃∗s
𝑅𝑐𝑜𝑛𝑑 = 4.76 ∗ 10−3
J
℃∗s
∆𝑇 = 4; 𝑅𝑇 = 9.26𝑥10−3 J

𝐽
𝑄̇𝑝𝑎𝑟𝑒𝑑,𝑠𝑢𝑟 = 131.97
𝑠

Pared oeste
Dimensiones a=6.4, h=1.6
• Temperatura corregida:
T∞ = 25 + 3.3
T∞ = 28.3°C
• Temperatura promedio:
28.3 + 40
𝑇𝑝𝑟𝑜𝑚 =
2
𝑇𝑝𝑟𝑜𝑚 = 34.1°C
• Viscosidad cinemática del aire a 34.1 ⁰C= 1.562*10-5 m2/s (Interpolación de
anexo 11)
• Número de Prandtl a 34.1 ⁰C: 0.7296 (Interpolación de anexo 11)
𝐽
• Conductividad térmica a 34.1⁰C: 0.02551 𝑠∗𝑚∗℃ (Interpolación de anexo 11)

Tenemos en cuenta que el espesor de las paredes es 0.1m y aplicando la ecuación 2:


e = 0.1 m
𝐽
k = 0.72 𝑠∗𝑚∗⁰𝐶

área = 10.24m2
℃∗s
𝑅𝑐𝑜𝑛𝑑 = 10.26 ∗ 10−3
J

33
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

℃∗s
∆𝑇 = 7.3; 𝑅𝑇 = 8.6𝑥10−3 J

𝐽
𝑄̇𝑝𝑎𝑟𝑒𝑑,𝑜𝑒𝑠𝑡𝑒 = 429.31
𝑠

Pared este
Dimensiones a=6.4, h=1.6
• Temperatura corregida:
T∞ = 25 + 3.3
T∞ = 28.3°C
• Temperatura promedio:
28.3 + 40
𝑇𝑝𝑟𝑜𝑚 =
2
𝑇𝑝𝑟𝑜𝑚 = 34.1°C
• Viscosidad cinemática del aire a 34.1⁰C= 1.562*10-5 m2/s (Interpolación de
anexo 11)
• Número de Prandtl a 34.1 ⁰C: 0.7296 (Interpolación de anexo 11)
𝐽
• Conductividad térmica a 34.1 ⁰C: 0.02551 𝑠∗𝑚∗℃ (Interpolación de anexo 11)

Tenemos en cuenta que el espesor de las paredes es 0.1m y aplicando la ecuación 2:


e = 0.1 m
𝐽
k = 0.72 𝑠∗𝑚∗⁰𝐶

área = 10.24m2
℃∗s
𝑅𝑐𝑜𝑛𝑑 = 10.26 ∗ 10−3
J
℃∗s
∆𝑇 = 7.3; 𝑅𝑇 = 8.6𝑥10−3 J

𝐽
𝑄̇𝑝𝑎𝑟𝑒𝑑,𝑒𝑠𝑡𝑒 = 429.31
𝑠

34
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Refrigerante: Amoniaco

Para analizar el ciclo de refrigeración del AMONÍACO, se debe primero, tener las
siguientes consideraciones respecto a tanto la temperatura de condensación como la de
evaporación.
Tevaporación = Tfinal – (5 - 8) °C
Tcondensación = Tambiente + (10 - 20) °C
𝑇𝑒𝑣𝑎𝑝𝑜𝑟𝑎𝑐𝑖ó𝑛 1 = −25 − 5 = −30 °𝐶
𝑇𝑐𝑜𝑛𝑑𝑒𝑛𝑠𝑎𝑐𝑖ó𝑛 = 25 + 15 = 40 °𝐶
𝑇𝑒𝑣𝑎𝑝𝑜𝑟𝑎𝑐𝑖ó𝑛 2 = 21 − 5 = 16 °𝐶

Punto 1:
𝑇1 = −30°𝐶 (𝑇𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑑𝑒𝑙 𝑒𝑣𝑎𝑝𝑜𝑟𝑎𝑑𝑜𝑟 1)
𝑃1 = 465.99𝑘𝑃𝑎 (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝐽
ℎ1 = 1445.8 (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝑔
𝑚3
𝑣1 = 0.2680 (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝑔
𝐾𝐽
𝑠1 = 5.31 (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝑔𝐾
Punto 2:
𝑇2 = 34.13°𝐶 (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝑃2 = 𝑃3 = 752.79𝑘𝑃𝑎 (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝐽
ℎ2 = 1509.30 (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝑔
𝐾𝐽
𝑠2 = 𝑠1 = 5.31 (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝑔𝐾
0.1772𝑚3
𝑣2 = (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝑔

35
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Punto 3:
𝑇3 = 16°𝐶 (𝑇𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑑𝑒𝑙 𝑒𝑣𝑎𝑝𝑜𝑟𝑎𝑑𝑜𝑟 2 )
𝑃3 = 752.79𝑘𝑃𝑎(𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝐽
ℎ3 = 1458.5 (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝑔
0.1693𝑚3
𝑣3 = (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝑔
𝐾𝐽
𝑠3 = 5.138 (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝑔𝐾

Punto 4:
𝑇4 = 67.33°𝐶 (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝑃4 = 𝑃5 = 1554.33𝑘𝑃𝑎(𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝐽
ℎ4 = 1558.26 (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝑔
𝐾𝐽
𝑠4 = 𝑠3 = 5.138 (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝑔𝐾
0.0951𝑚3
𝑣4 = (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝑔

Punto 5:
𝑇5 = 40°𝐶 (𝑇𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑑𝑒𝑙 𝑐𝑜𝑛𝑑𝑒𝑛𝑠𝑎𝑑𝑜𝑟 )
𝑃5 = 1554.33𝑘𝑃𝑎(𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝐽
ℎ5 = 371.7 (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝑔
0.001726𝑚3
𝑣5 = (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝑔
𝐾𝐽
𝑠5 = 1.3591 (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝑔𝐾

36
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Punto 6:
𝐾𝐽
ℎ6 = ℎ5 = 371.7 (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝑔

Punto 7:
𝑇7 = 𝑇3 = 16°𝐶 (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)

𝑃7 = 𝑃3 = 752.79𝑘𝑃𝑎(𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)

𝐾𝐽
ℎ7 = 255.9 (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝑔

𝐾𝐽
𝑠7 = 0.9788 (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝑔𝐾

0.001623𝑚3
𝑣7 = (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝑔

Punto 8:
𝐾𝐽
ℎ8 = ℎ7 = 255.9 (𝑇𝑎𝑏𝑙𝑎𝑠 𝑡𝑒𝑟𝑚𝑜𝑑𝑖𝑛á𝑚𝑖𝑐𝑎𝑠)
𝐾𝑔

40 °𝐶

16°𝐶

2.2 °𝐶

Figura 3 Diagrama presión-entalpía del ciclo ideal de refrigeración


Fuente: Stoecker-Refrigeracion
37
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

40 °𝐶

16 °𝐶

2.2 °𝐶

Figura 4 Esquema del sistema de refrigeración de 2 evaporadores y 2 compresores


Fuente: Stoecker-Refrigeracion

Calculando el ERU según ecuación 11:


𝐸𝑅𝑈 = ℎ1 − ℎ8
𝐸𝑅𝑈 = 1189.9 𝐾𝐽/𝑘𝑔
Hallando flujo másico 1, según ecuación 12:
̇
𝑄𝑡𝑜𝑡𝑎𝑙 𝑐𝑎𝑚𝑎𝑟𝑎
̇ =
𝑚1
𝐸𝑅𝑈
𝐽
𝑄𝑡𝑜𝑡𝑎𝑙 𝑐𝑎𝑚𝑎𝑟𝑎 = 111 913.04
𝑠
𝑚̇ = 0.941𝐾𝑔/𝑠
𝑚̇1= 𝑚̇2= 𝑚̇7= 𝑚̇8
Hallamos el segundo flujo másico con balance de energía:
𝐽 𝐾𝐽
𝑚̇5 ℎ5 + 𝑚̇2 ℎ2 + 111 913.04 = 𝑚̇3 ℎ3 + 𝑚̇7 ℎ7
𝑠 𝑠

38
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Sabemos que por balance de masas tenemos:


𝑚̇5 + 𝑚̇2 = ̇ 𝑚7 + 𝑚̇3
De esto sacamos que:
𝑚̇5 =̇ 𝑚̇3
Por el balance de energía:
𝑚̇5 =̇ 𝑚̇3 = 0.449 𝐾𝑔/𝑠
Hallando trabajo de compresión de baja, según ecuación 13:
𝑊1 = ℎ2 − ℎ1 = 1509.3 − 1445.8 = 63.5𝐾𝐽/𝐾𝑔
Hallando trabajo de compresión de alta, según ecuación 13:
𝑊 = ℎ4 − ℎ3 = 99.76𝐾𝐽/𝐾𝑔
Calculando la potencia:
𝐾𝑔 𝐾𝐽 𝐾𝐽
𝑝𝑜𝑡𝑒𝑛𝑐𝑖𝑎 𝑐𝑜𝑚𝑝𝑟𝑒𝑠𝑜𝑟 𝑑𝑒 𝑏𝑎𝑗𝑎 = 𝑃𝑏 = 0.119 ∗ 63.5 = 7.56 = 10.25𝐻𝑃
𝑠 𝐾𝑔 𝑠
𝐾𝑔 𝐾𝐽 𝐾𝐽
𝑝𝑜𝑡𝑒𝑛𝑐𝑖𝑎 𝑐𝑜𝑚𝑝𝑟𝑒𝑠𝑜𝑟 𝑑𝑒 𝑎𝑙𝑡𝑎 = 𝑃𝑎 = 0.267 ∗ 99.76 = 26.64 = 35.71𝐻𝑃
𝑠 𝐾𝑔 𝑠

Aplicando un factor de seguridad de 1.2 a los 2 compresores tenemos un compresor de


baja de 15 hp y 40 hp.

Calculando el flujo de calor en el condensador, según ecuación 17 si:


h4= 1558.26 kJ/kg; h5 = 371.71 kJ/kg
𝑚̇5 =̇ 𝑚̇3 = 0.249 𝐾𝑔/𝑠
𝐾𝐽
̇
𝑄𝑐𝑜𝑛𝑑𝑒𝑛𝑠𝑎𝑑𝑜𝑟 = 295.45
𝑠

39
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Cálculo de tuberías:

Considerando una instalación comercial mediana (véase tabla 8) y para el diámetro


(véase ecuación 18):

Líneas de succión:
𝑚3
Flujo másico 𝐾𝑔/𝑠 Volumen especifico 𝐾𝑔 Velocidad m/s

𝑚̇ = 0.941 𝑣1 = 0.2680 7.5

∅ = 6.2cm

𝑚3
Flujo másico 𝐾𝑔/𝑠 Volumen especifico 𝐾𝑔 Velocidad m/s

0.249 𝑣3 = 0.1693 7.5

∅ = 7.1cm

Línea de descarga:
𝑚3
Flujo másico 𝐾𝑔/𝑠 Volumen especifico 𝐾𝑔 Velocidad m/s

𝑚̇ = 0.119 𝑣2 = 0.1772 10

∅ = 5cm

𝑚3
Flujo másico 𝐾𝑔/𝑠 Volumen especifico 𝐾𝑔 Velocidad m/s

0.249 𝑣4 = 0.0951 10

∅ = 5.2cm

40
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Línea de líquido:
𝑚3
Flujo másico 𝐾𝑔/𝑠 Volumen especifico 𝐾𝑔 Velocidad m/s

0.249 𝑣5 = 0.001726 1

∅ = 2cm

Tabla 9 Velocidades recomendadas en las líneas

Fuente: Diapositivas proporcionadas en clase

41
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Selección de equipo:

Compresor:

Se selecciona un compresor con una potencia de 15 y 40 HP, para discernir entre el tipo
de compresor a emplear, según su potencia escogemos el tipo basándonos en la tabla 9.

Tabla 10.Selección de compresores según su potencia de trabajo

Fuente: Guía práctica: Selección de los elementos de una instalación frigorífica

CV quiere decir caballos de vapor y tienen un valor similar a los HP, por lo que será
para el compresor de 15 hp uno de pistón y para el de 40 hp un compresor de tornillo.
Elegimos el compresor según el catálogo de Compresores Josval, véase anexo 5 y 6.
Finalmente, nuestro compresor será: Compresor Josval, modelo MISTRAL 40-A AS y
modelo MC-AI-1000/A

Figura 5.Compresor Josval MISTRAL 40-A-AS


Fuente: Imagen de catálogo de compresores Josval

42
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Figura 6.Compresor Josval MC-AI-1000/A


Fuente: Imagen de catálogo de compresores Josval

Evaporador:
Se selecciona un evaporador de acuerdo con la potencia del compresor. Nos basamos
en el catálogo de Delta Frio (véase anexo 7) para la elección de nuestro modelo de
evaporador. Según la tabla adjunta en el anexo 7, para nuestra elección nos basamos en la
temperatura de evaporación y la potencia en kcal/h.
Finalmente, nuestro evaporador será:

Figura 7.Evaporador DeltaFrio, modelo Casset 232


Fuente: DeltaFrio

43
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Figura 8.Evaporador DeltaFrio, modelo DFC73.74


Fuente: DeltaFrio

Condensador:
Basamos nuestra elección en el catálogo de condensadores García Cámara (véase anexo
8), según la potencia del compresor de alta (40 HP).
Por lo que, nuestro condensador será: Condensador García Cámara, modelo CG30.

Figura 9.Condensador García Camara, modelo CG30


Fuente: Catálogo general García Camara

Válvula de expansión:
Para la selección de las válvulas de expansión, buscamos uno que sea adecuado para el
refrigerante elegido, en nuestro caso, amoniaco.

44
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Según el anexo 10, nos encontramos en el rango -20 °C a 30 °C y finalmente la válvula


de expansión elegida será: Válvula de expansión Dnfoss, modelo TEA 20-12.

Figura 10 Válvula de expansión Danfoss TEA 20-12


Fuente: Catálogo Danfoss, 2016

Filtro de secado:
Para su selección hay que tener en cuenta el diámetro de línea de líquido, también su
presión y temperatura de trabajo (véase anexo 12).
Finalmente, basándonos en el diámetro adecuado, se elige filtro de secado a: Filtro de
secado Danfoss, modelo DCL 167s.

Figura 11.Filtro de secado DCL 167s


Fuente: Total Parts

45
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Tanque recibidor:
Los sistemas refrigeración están expuestos a cambios de temperatura externa o a
variaciones en la carga térmica, y para ese tipo de escenarios, se utiliza el tubo recibidor,
para almacenar el exceso de refrigerante producido por estas condiciones mencionadas.

También provee en forma constante e ininterrumpida refrigerante en estado líquido a


la entrada de la válvula de expansión. Es la de almacenar toda la carga de refrigerante de
la instalación, debido a tareas de mantenimiento o reparación. Para la selección de este
componente, se debe de tener de referencia la presión y temperatura de trabajo,
dimensiones de conexión y capacidad.
El tanque que elegimos para nuestra cámara refrigerante será: Tanque recibidor NR
RefrigeraCION, modelo CW-644 (figura 15), cuyas características son:

Figura 12.Características del tanque recibidor CW-644

POTENCIA 3-5HP

Fuente: propio

Figura 13.Tanque recibidor CW-644


Fuente: NR REFRIGERACION

46
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Acumulador de succión:
Acumuladores de succión para proteger al compresor, evitando el regreso de
refrigerante líquido al mismo. El acumulador de succión se colocará en la línea de líquido.

Por lo que hay que tener en cuenta el tamaño de la tubería, también hay que considerar
la presión, temperatura.

Tabla 11.Modelos de diferentes acumuladores de succión

Fuente: IMCOSA

Basándonos en la tabla 10 y en las medidas de diámetro estándares, se seleccionó como


acumulador de succión: Acumulador de succión IMCOSA, modelo 2842-AS-5117-7

47
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Figura 14. Acumulador de succión 2842-AS-5117-7


Fuente: IMCOSA

Visor de refrigerante:
Control inmediato y directo de la circulación, del estado y cantidad de humedad del
refrigerante en fase líquida o bifásica, en las instalaciones de refrigeración y
acondicionamiento de aire; para seleccionar un visor adecuado, hay que tener en cuenta el
diámetro de línea de líquido (véase anexo 13). Finalmente, para nuestro caso, se seleccionó:
Visor de refrigerante Carly, modelo VCYL 33S.

Figura 15.Visor de refrigerante VCYL 33S


Fuente: Carly

Separador de aceite.
Cumple la función de regresar el aceite al compresor, manteniendo así el refrigerante
libre de contaminación, también ahorramos energía y costo de operación. Para su selección
nos basaremos en el diámetro de descarga, y el separador de aceite seleccionado es:
Separador de aceite Emerson, modelo 300-SAS-4

Tubo anti vibrador


Para evitar las vibraciones mecánicas causadas por el compresor en el momento en que
realiza trabajo, se acoplo tubos anti-vibradores. Según la dimensión del acoplamiento al

48
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

compresor y la dimensión del tubo de descarga. Véase la tabla del anexo 15. Finalmente,
el tubo anti-vibrador seleccionado será: Tubo anti vibrador Carly, modelo EVCYDEAC21
S

Figura 16.Tubo anti vibrador EVCYDEAC 21 S


Fuente: Mercado Libre

Presostato de aceite:
• Presostato de alta presión: Este elemento actúa cuando la presión de condensación
es muy alta, protegiendo el sistema contra la sobrecarga o una posible rotura de las tuberías,
parando el compresor. Se conecta en la línea de descarga, a la salida del compresor.
• Presostato de baja presión: Actuará cuando la presión en la línea de aspiración sea
inferior al nivel establecido. Esta caída de presión puede causar una temperatura de
evaporación muy baja, dañando los productos almacenados o la propia cámara.
Para seleccionar un presostato hay que tener en cuenta la presión (véase tabla 13),
finalmente, el presostato elegido es: Presostato de aceite Danfoss, modelo KP15. (figura
19).

Figura 17 Presostato de aceite KP15

49
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Tabla 12.Modelos de presostato marca Danfoss según el tipo de presión

Fuente: Catálogo de presostatos Danfoss, 2013.

Termostatos.
Permite controlar la temperatura del recinto y por ende de los productos. Este
dispositivo se instala al interior de la cámara. El instante en que la planta frigorífica debe
empezar a funcionar o pararse, viene determinado por este sensor. Parando o poniendo en
marcha el motor del compresor, en función del nivel térmico. Para un correcto
funcionamiento de la cámara, es imprescindible que dichos valores estén bien definidos.
Finalmente, el termostato empleado será: Termostatos Danfoss, modelo EKC 202A.
(figura 15)

50
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Figura 18 Termostatos EKC 202ª

Tabla 13.Tres tipos de termostatos marca Danfoss y características

Fuente: Catálogo termostatos Danfoss, 2016.

51
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Capítulo IV
Resultados y Discusiones
o Carga térmica total: 184 656.51 KJ/s (cámara)-121.629 KJ/s(oficina)

o Flujo másico de refrigerante: 0.119Kg/s y 0.249 Kg/s

o Potencia del compresor de baja: 15HP

o Potencia del compresor de alta: 40HP

Elementos seleccionados:

o Compresor: Marca Josval, Modelo Mistral 40-A AS y MC-AI-1000/A

o Condensador: Marca García Cámara, modelo CG30

o Evaporador: Marca Delta Frío, modelo DF Casset 232 y DFC73.74

o Válvula de expansión: Marca Danfoss, modelo TEA-20-12.

o Filtro de secado: Marca Danfoss, modelo DCL 167s

o Tanque recibidor: Marca NR, modelo CW-644

o Acumulador de succión: Marca Imcosa, modelo 2844-AS-5134-7

o Visor de refrigerante: Marca Carly, modelo VCYL 33S

o Separador de aceite: Marca Emerson, Modelo 300-SAS-3

o Válvula de seguridad solenoide: Marca Danfoss, Modelo E14S250

o Tubo anti vibrador: Marca Carly, Modelo EVCYDEAC 2 S

o Presostato de aceite: Marca Danfoss, Modelo KP15

o Termostatos: Marca Danfoss, Modelo EKC 202A

52
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Capítulo V
Conclusiones
- Se consiguió proveer una solución de mantenimiento y almacenamiento de mariscos
realizando los cálculos respectivos, dimensionamiento adecuado y correcta elección
de los elementos de la cámara de refrigeración.
- Un cuarto de 60 metros cúbicos fue diseñado para mantener 1tonelada de maricos.
- En base a las dimensiones del cuarto y a las distintas fuentes de calor que pueden
afectar el proceso de enfriamiento se llegó a una carga térmica de 184 656.51KJ/s y
para la oficina la carga térmica es de 121.629 KJ/s.
- Se determinó la temperatura de evaporación 1, de 2.2°C, de evaporación 2 , de 16°C y
condensación, de 40°C, para hallar los estados del ciclo de refrigeración.
- Con ayuda de las temperaturas halladas y tablas termodinámicas del amoníaco se
obtuvieron los valores de presión, temperatura, entalpía, entropía y volumen específico
de cada estado termodinámico.
- El flujo másico del sistema resultó de 0.119 Kg/s y 0.249 Kg/s
- Se dimensionaron las tuberías resultando con diámetros de 7 y 8.5cm las líneas de
succión, 5 y 6 cm las líneas de descarga y 2cm la línea de líquido.
- Se seleccionaron los siguientes elementos secundarios para proveer seguridad en el
ciclo de refrigeración:
o Filtro de secado: Marca Danfoss, modelo DCL 167s
o Tanque recibidor: Marca NR, modelo CW-644
o Acumulador de succión: Marca Imcosa, modelo 2844-AS-5134-7
o Visor de refrigerante: Marca Carly, modelo VCYL 33S
o Separador de aceite: Marca Emerson, Modelo 300-SAS-3
o Válvula de seguridad solenoide: Marca Danfoss, Modelo E14S250
o Tubo anti vibrador: Marca Carly, Modelo EVCYDEAC 2 S
o Presostato de aceite: Marca Danfoss, Modelo KP15
o Termostatos: Marca Danfoss, Modelo EKC 202A

53
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Capítulo VI
Recomendaciones

o Hacer los cálculos para un ciclo real, para obtener una mayor precisión.
o Considerar subenfriamiento y sobre calentamiento, con el fin de proteger el
compresor.
o Emplear refrigerantes más comerciales, debido a que suelen ser más frecuentes los
equipos para los refrigerantes comerciales, y es muy dificultoso encontrar equipos
para refrigerantes ya no empleados.
o Procurar hacer la selección de los elementos de la cámara de refrigeración de un
mismo proveedor (marca), esto ya que en muchas ocasiones algunos elementos de
marcas distintas son incompatibles entre ellos, es mejor adquirirlos todos de una
misma entidad.
o Instalar una o dos láminas de polietileno entre la superficie superior de la capa de
aislamiento y la superficie inferior de las losas de concreto, para evitar la
penetración de vapor de agua.
o Para el buen funcionamiento del ciclo del aire y de la temperatura, la carga debe ir
paletizada para que circule el aire por debajo, ser uniforme y cubrir toda la superficie
del suelo. En caso de usar cajas estas deben tener ventilación en forma de agujeros.

54
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Referencias Bibliográficas

• Catálogos de siñadores
• Teoria de congelacion y refrigeracion de pescados y marisos..
• Castillo, J. F. (2015). Diseño de camara frigorifica.
• Ceballos Maya, J. (2016). Diseño óptimo de una cámara de conservación de
productos congelados. Sevilla.
• Fundación Española de la Nutrición. (18 de 07 de 2018). Fen.org.es. Obtenido de
http://www.fen.org.es/mercadofen/pdfs/frambuesa.pdf
• Martínez Salinas, J. J. (2016). DISEÑO DE UNA CÁMARA FRIGORÍFICA PARA
LA CONSERVACIÓN DE POLLOS EN LA AGENCIA DE VENTA PIO LINDO
DE LA CIUDAD DE EL ALTO. La Paz.
• Universidad de Lima. (06 de 02 de 2012). Fresno.ulima.edu.pe. Obtenido de
http://fresno.ulima.edu.pe/sf%5Csf_bdfde.nsf/imagenes/9A5B9CD541FA172005
2573540070AE16/$file/10-25-ayta.pdf

55
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Anexo
ANEXO 1: Calor específico de alimentos

Fuente: Módulo de refrigeración industrial, bdigital.

56
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

ANEXO 2: Coeficiente de transmisión de calor para diferentes tipos de materiales

Fuente: Módulo de refrigeración industrial, Tecsup

57
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

ANEXO 3: Ganancia de calor por infiltración de aire.

Fuente: Módulo de refrigeración industrial, Tecsup

58
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

ANEXO 4: Calor de respiración de vegetales y frutas

Fuente: Calor de respiración de frutas y vegetales, Carlos E. Orrego.

59
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

ANEXO 5: Tipos de compresores de tornillo de marca Josval

Fuente: Catálogo de Compresores Josval, 2014

60
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

ANEXO 6: Tipos de compresores de pistón de marca Josval

Fuente: Catálogo de Compresores Josval, 2014

61
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

ANEXO 7: Tipos de evaporadores para bodegas climatizadas marca delta frio

Fuente: Catálogo de evaporadores Delta Frio

62
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

ANEXO 8: Tipos de evaporadores para bodegas climatizadas marca delta frio

Fuente: Catálogo de evaporadores Delta Frio

63
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

ANEXO 9: Modelos de condensadores de marca García Camara

Fuente: Catálogo de condensadores García Camara, 2017.

64
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

ANEXO 10: Válvulas de expansión termostáticas para amoníaco, tipo TEA

Fuente: Folleto técnico de válvulas de expansión marca Danfoss.

65
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

ANEXO 11: Propiedades del aire a la presión de 1 atm

Fuente: Transferencia de calor y masa, Yunus – Cengel, Tercera edición

66
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

ANEXO 12: Filtros de secado

Fuente: Catálogo de revista de filtros de secado Danfoss

67
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

ANEXO 13: Diferentes modelos de visores de refrigerantes según diferentes

diámetros

Fuente: Carly

68
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

ANEXO 14: Tabla de selección de capacidades para separadores de aceite

Fuente: IMCOSA

69
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

ANEXO 15: Modelos de tubos antivibradores marca Carly

Fuente: Catálogo de tubos antivibradores marca Carly, 2015.

70
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

Test de comprobación
• ¿Cuáles son las razones desde el punto de vista técnico, económico, ambiental y
social de la selección del tipo de refrigerante a usar en su sistema de
refrigeración? (comparar mínimo con otros dos tipos de refrigerantes).

Las soluciones que usamos de refrigerante basadas en amoniaco van a costar entre un
10% y un 20% respecto a los CFC (clorofluorocarbonos), el cual se suele usar en los aires
acondicionados y es uno de los refrigerantes más usados . El amoniaco es un refrigerante
mucho más eficiente que los CFC (clorofluorocarbonos), por lo que un sistema con
amoniaco consumirá menos electricidad y por tanto menos costos de operación. El
amoniaco disfruta rendimientos termodinámicos en torno a un 3-10% superiores a otros
refrigerantes, por lo bajo las mismas premisas de capacidad frigorífica tiene un consumo
energético menor. Es un refrigerante seguro para el medio ambiente y necesario para
impulsar la sostenibilidad, con un potencial de agotamiento del ozono de 0 y un potencial
calentamiento global de 0.
Lo mismo pasa con los derivados de los clorofluorocarbonos, la principal ventaja del
amoniaco es su bajo costo y su amabilidad con el ambiente, ya que es de los que menos
contaminan

• ¿Cuáles son las razones desde el punto de vista técnico, económico, ambiental y
social de la selección del tipo de compresor (tornillo, alternativo, scroll, etc.) de
su sistema de refrigeración? (comparar mínimo con otros dos tipos de
compresores).

o Para nuestro sistema de refrigeración trabajamos con dos compresores para obtener un
mejor rendimiento ya que los dos compresores trabajarían menos utilizando menor
potencia eléctrica reduciendo costo a largo plazo con una mejor condición de
operación desde el punto de vista mecánico y de mantenimiento
o Nuestro tipo de compresores elegidos son de tipo tornillo y otro de pistón ya que estos
dos compresores se complementan entre si ya que el del tipo tornillo es más usado de
71
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO

manera industrial y el de pistón en sistema de refrigeración con menor impacto por eso
el compresor de tipo pistón es usado en la primera fase de compresión y el de tornillo
en la segunda ya que tiene una mayor durabilidad y eficiencia mecánica
o Al comparar los compresores elegidos con uno del tipo scroll podemos deducir que el
de tipo scroll es un compresor para mayores rangos de trabajo y con un costo mayor
de uno de tipo tornillo, pero para nuestras especificaciones de sistema de refrigeración
calculadas no sería muy factible utilizar uno de ellos
o Por otro lado al comparar los compresores elegidos con uno de tipo rotativo,
encontramos que el de tipo rotativo tiene muy poca versatilidad lo que hace que su
uso sea limitado, una vez instalado este tipo de compresor es difícil su traslado lo que
hace que sea dependiente del equipo al que se encuentre conectado, otra desventaja de
usar este tipo de compresor seria que sus reparaciones son muy complicadas ya que si
sufren de algún daño o desperfecto la reparación puede salirte muy costosa o incluso
llegando a ser necesario un cambio completo.

72
INGENIERÍA MECATRÓNICA

También podría gustarte