UNIVERSIDAD TÉCNICA DE AMBATO

FACULTAD DE INGENIERÍA CIVIL Y MECÁNICA

CARRERA DE INGENIERÍA MECÁNICA

TRABAJO EXPERIMENTAL PREVIO A LA OBTENCIÓN DEL TÍTULO

DE INGENIERO MECÁNICO

TEMA:

“ANÁLISIS TÉRMICO DE UN PROTOTIPO DE CÁMARA DE SECADO

ESTÁTICA POR CONVECCIÓN Y RADIACIÓN DE 1000 W PARA USO
SEMIINDUSTRIAL”

AUTOR: William Rene Criollo Camacho

TUTOR: Ing. Mg. Santiago Paúl Cabrera Anda

AMBATO – ECUADOR
Abril - 2021
 CERTIFICACIÓN

En mi calidad de Tutor del Trabajo Experimental, previo a la obtención del título de

Ingeniero Mecánico, con el tema: “ANÁLISIS TÉRMICO DE UN PROTOTIPO
DE CÁMARA DE SECADO ESTÁTICA POR CONVECCIÓN Y
RADIACIÓN DE 1000 W PARA USO SEMIINDUSTRIAL”, elaborado por el
Sr. William René Criollo Camacho, portador de la cedula de ciudadanía: C.I.
1804640587, estudiante de la Carrera de Ingeniería Mecánica, de la Facultad de
Ingeniería Civil y Mecánica.

Certifico:

• Que el presente trabajo experimental es original de su autor.

• Ha sido revisado cada uno de sus capítulos componentes.
• Esta concluido en su totalidad.

Ambato, Marzo 2021

SANTIAGO Firmado digitalmente por

SANTIAGO PAUL CABRERA ANDA
Nombre de reconocimiento (DN):

PAUL c=EC, o=BANCO CENTRAL DEL

ECUADOR, ou=ENTIDAD DE
CERTIFICACION DE INFORMACION-

CABRERA ECIBCE, l=QUITO,

serialNumber=0000497368,
cn=SANTIAGO PAUL CABRERA

ANDA
ANDA
Fecha: 2021.03.24 13:31:05 -05'00'

Ing. Mg. Santiago Paúl Cabrera Anda

TUTOR

ii
 AUTORÍA DE LA INVESTIGACIÓN

Yo, William René Criollo Camacho, con C.I. 1804640587 declaro que todas las
actividades y contenidos expuestos en el presente trabajo experimental con el tema:
“ANÁLISIS TÉRMICO DE UN PROTOTIPO DE CÁMARA DE SECADO
ESTÁTICA POR CONVECCIÓN Y RADIACIÓN DE 1000 W PARA USO
SEMIINDUSTRIAL”, así como también los análisis estadísticos, gráficos,
conclusiones y recomendaciones son de mi exclusiva responsabilidad como autor del
proyecto, a excepción de las referencias bibliográficas citadas en el mismo.

Ambato, Marzo 2021

William René Criollo Camacho

C.I. 1804640587
AUTOR

iii
 DERECHOS DE AUTOR

Autorizo a la Universidad Técnica de Ambato, para que haga de este trabajo

experimental o parte de él, un documento disponible para su lectura, consulta y
procesos de investigación, según las normas de la Institución.

Cedo los derechos en línea patrimoniales de mi Trabajo Experimental, con fines de

difusión pública, además apruebo la reproducción de este documento dentro de las
regulaciones de la Universidad, siempre y cuando esta reproducción no suponga una
ganancia económica y se realice respetando mis derechos de autor.

Ambato, Marzo 2021

……………………………………
William René Criollo Camacho
C.I.: 1804640587
AUTOR

iv
 APROBACIÓN DEL TRIBUNAL DE GRADO

Los miembros del Tribunal de Grado aprueban el informe del Trabajo Experimental,
realizado por el estudiante William René Criollo Camacho, de la Carrera de
Ingeniería Mecánica bajo el tema: “ANÁLISIS TÉRMICO DE UN PROTOTIPO
DE CÁMARA DE SECADO ESTÁTICA POR CONVECCIÓN Y
RADIACIÓN DE 1000 W PARA USO SEMIINDUSTRIAL”.

Ambato, Marzo 2021

Para constancia firman:

Digitally signed
by JAVIER
GARCIA DEL
ALEX SANTIAGO MAYORGA VALLE
PARDO Date:
2021-03-25 09:46-05:00 2021.03.27
08:10:59 ECT
----------------------------------------------- -------------------------------------------
Ing. Mg. Alex Santiago Mayorga Pardo Ing. Mg. Javier García del Valle
Miembro Calificador Miembro Calificador

v
 DEDICATORIA

Este proyecto va dedicado principalmente a Dios por darme salud y vida en estos
tiempos difíciles para el mundo entero.

A mis padres, Alfonso y Elena, por ser el pilar fundamental de ejemplo y apoyo para
seguir adelante en toda mi vida, padres los cuales nunca me han desamparado y me
han brindado todo lo necesario para lograr mis metas.

A mis hermanos que me han ayudado y apoyado en las circunstancias que se me han
presentado a lo largo de la vida.

A una persona especial quien ha estado conmigo en las buenas y malas desde que la
conocí, y me dio un hijo al cual amo con mi vida.

Willy

vi
 AGRADECIMIENTO

Agradezco a mis padres que con esfuerzo y sacrificio me supieron sacar adelante
para lograr esta meta ya que sin el apoyo de ellos no lo hubiera logrado. A mi familia
que me apoyaron moralmente para seguir siempre adelante.

Al Ingeniero Santiago Cabrera, tutor que supo compartir sus conocimientos los
cuales me ayudaron para el desarrollo del trabajo de titulación y culminación de
forma satisfactoria.

A la Universidad Técnica de Ambato, principalmente a la Carrera de Ingeniería

Mecánica por ayudarme a ser una persona preparada para ejercer la profesión de
manera adecuada.

A los docentes que me llenaron de conocimientos a lo largo de mis estudios.

Agradezco a mi hermana Ximena Criollo por ser una inspiración y ejemplo a seguir
como persona y como profesional que siempre estuvo presente brindándome apoyo
incondicional.

Agradezco a la madre de mi hermoso hijo por apoyarme y tenerme paciencia a lo

largo del tiempo que hemos compartido juntos.

Agradezco principalmente a mi Madre Elena que hizo todo lo posible para que este
objetivo se haga realidad, este propósito alcanzado es suyo por su apoyo
incondicional le agradezco de todo corazón.

A las personas que conocí al transcurso de mi vida estudiantil, a los verdaderos

amigos les agradezco por compartir tantas cosas conmigo.

vii
ÍNDICE GENERAL DE CONTENIDOS
CERTIFICACIÓN ..................................................................................................................ii
AUTORÍA DE LA INVESTIGACIÓN ................................................................................iii
DERECHOS DE AUTOR ..................................................................................................... iv
APROBACIÓN DEL TRIBUNAL DE GRADO ..................................................................v
DEDICATORIA ..................................................................................................................... vi
AGRADECIMIENTO .......................................................................................................... vii
ÍNDICE GENERAL DE CONTENIDOS .......................................................................... viii
ÍNDICE DE TABLAS........................................................................................................... xii
ÍNDICE DE FIGURAS........................................................................................................ xiii
RESUMEN............................................................................................................................. xv
ABSTRACT .......................................................................................................................... xvi
CAPÍTULO I .......................................................................................................................... 1
MARCO TEÓRICO .............................................................................................................. 1
1.1 Antecedentes Investigativos. ....................................................................................... 1
1.2 Objetivos ....................................................................................................................... 3
1.2.1 Objetivo general .................................................................................................... 3
1.2.2 Objetivos específicos ............................................................................................. 3
1.3 Fundamentación teórica .............................................................................................. 4
1.3.1 Transferencia de calor .......................................................................................... 4
1.3.2 Calor. .................................................................................................................... 10
1.3.3 Balance de energía para sistemas de flujo estacionario ................................... 11
1.3.4 Convección forzada interna ............................................................................... 11
1.3.5 Número de Reynols ............................................................................................. 11
1.3.6 Número de Nusselt .............................................................................................. 12
1.3.7 Longitud característica ....................................................................................... 12
1.3.8 Convección libre interna o convección interna ................................................ 13
1.3.9 Número de Grashof............................................................................................. 13
1.3.10 Número de Rayleigh.......................................................................................... 14
1.3.11 Maneras de transferencia de calor por radiación. ......................................... 15
1.3.12 Radiación infrarroja ......................................................................................... 18
1.3.13 Dinámica de fluidos computacional (CFD) ..................................................... 18
1.3.14 Estructura de Dinámica de Fluidos Computacional (CFD) .......................... 19
1.3.15 Definición de secadores ..................................................................................... 20
1.3.16 Tipos de secadores infrarrojos ......................................................................... 21

viii
CAPÍTULO II ...................................................................................................................... 25
METODOLOGÍA ................................................................................................................ 25
2.1 NIVELES O TIPOS DE INVESTIGACIÓN........................................................... 25
2.1.1 Descriptiva ........................................................................................................... 25
2.1.2 Correlacional ....................................................................................................... 25
2.1.3 Experimental ....................................................................................................... 26
2.1.4 Bibliográfica ........................................................................................................ 26
2.1.5 Exploratoria......................................................................................................... 26
2.2 PROCESO DEL TRABAJO DE INVESTIGACIÓN............................................. 27
2.3 DIMENSIONES GEOMÉTRICAS Y CARACTERISTICAS DEL MODELO .. 28
2.3.1 Dimensiones de la cámara de secado ................................................................. 28
2.3.2 Aislante (lana de vidrio) ..................................................................................... 28
2.3.3 Carcasa del secador ............................................................................................ 29
2.4 VERIFICACIÓN DEL MODELO PARA LA SIMULACIÓN ............................. 30
2.5 PASOS A SEGUIR PARA LA ENTRADA DE DATOS A LA SIMULACIÓN .. 31
2.5.1 Chequeo del modelo ............................................................................................ 31
2.5.2 Principio de la simulación .................................................................................. 31
2.5.4 Tipo de análisis .................................................................................................... 32
2.5.5 Parámetros para el tipo de análisis.................................................................... 33
2.5.6 Tipo de fluido....................................................................................................... 33
2.5.7 Condiciones iniciales ........................................................................................... 34
2.5.8 Dominio Computacional ..................................................................................... 34
2.5.9 Selección de materiales ....................................................................................... 35
2.5.10 Condiciones de frontera.................................................................................... 36
2.5.11 Cantidad de calor .............................................................................................. 37
2.5.12 Superficies radiantes ......................................................................................... 38
2.5.13 Calor por radiación ........................................................................................... 38
2.6 ANÁLISIS EXPERIMENTAL ................................................................................. 39
2.6.1 Instrumentos para la obtención de información .............................................. 39
2.7 OBTENCIÓN DE LA INFORMACIÓN ................................................................. 42
CAPÍTULO III ..................................................................................................................... 43
ANÁLISIS Y DISCUSIÓN DE RESULTADOS ............................................................... 43
3.1 Obtención de datos ..................................................................................................... 43
3.1.1 Parámetros de funcionamiento .......................................................................... 43
3.1.2 Cálculos térmicos ................................................................................................ 44

ix
 3.1.3 Cálculo del coeficiente de transferencia de calor para Convección Natural
Interior de la cámara ................................................................................................... 44
3.1.4 Calor generado por convección natural en el interior de la cámara .............. 47
3.1.5 Cálculo del coeficiente de transferencia de calor para Convección Forzada 48
3.1.6 Cálculo de calor generado por radiación al interior de la cámara ................. 50
3.2 CÁLCULO DE PÉRDIDAS POR CONVECCIÓN AL EXTERIOR DE LA
CÁMARA ......................................................................................................................... 51
3.2.1 Análisis del coeficiente de convección natural al exterior del secador ........... 51
3.2.2 Cálculo de la pérdida de calor por convección al exterior del secador
(radiación despreciable) .............................................................................................. 53
3.3 CÁLCULO DE LA RESISTENCIAS TÉRMICAS POR CONDUCCIÓN Y
CONVECCION EN LAS PAREDES DEL PROTOTIPO DE LA CÁMARA DE
SECADO ........................................................................................................................... 53
3.3.1 Resistencia térmica total por conducción ......................................................... 57
3.3.2 Resistencia térmica por convección natural al exterior del secador (radiación
despreciable) ................................................................................................................. 57
3.3.3 Cálculo del calor perdido total por las paredes del prototipo de secador ...... 58
3.4 CÁLCULO DE LA VARIACION DE LA ENERGIA INTERNA
ALMACENADA EN EL AIRE EN LA CÁMARA DEL SECADOR ........................ 58
3.4.1 Cálculo del calor que se necesita para elevar la temperatura del aire en el
interior de la cámara de secado (qair). ........................................................................ 58
3.4.2 Cálculo de la variación de la energía interna almacenada .............................. 59
3.5 BALANCE DE ENERGÍA ........................................................................................ 60
3.5.1 Cálculo del calor total de entrada requerido .................................................... 60
3.6 CÁLCULO DE LA EFICIENCIA DEL PROCESO DE SECADO (η) ................ 61
3.7 ETAPAS DE SIMULACIÓN TÉRMICA................................................................ 61
3.7.1 Validación del modelo......................................................................................... 62
3.8 ANÁLISIS DE LA DISTRIBUCIÓN (GLOBAL) DE TEMPERATURA REAL
AL INTERIOR DE LA CÁMARA DE SECADO ........................................................ 63
3.8.1 Temperaturas en el centro de la cámara de secado ......................................... 63
3.8.2 Temperaturas en la parte inferior de la cámara de secado ............................. 64
3.9 SIMULACIÓN DE DISTRIBUCIÓN DE TEMPERATURAS AL INTERIOR DE
LA CÁMARA EN ESTADO TRANSITORIO.............................................................. 66
3.9.1 Convección y radiación en el interior de la cámara en estado transitorio para
un tiempo de 180s ......................................................................................................... 67
3.9.2 Convección y radiación en el interior de la cámara en estado transitorio para
un tiempo de 420s ......................................................................................................... 70
3.9.3 Convección y radiación en el interior de la cámara en estado transitorio para
un tiempo de 688s ......................................................................................................... 74

x
 3.10 VALIDACIÓN DEL COMPORTAMIENTO DEL SISTEMA ........................... 77
3.11 SIMULACIÓN DE DISTRIBUCIÓN DE TEMPERATURAS AL INTERIOR
DE LA CÁMARA EN ESTADO ESTABLE ................................................................. 78
3.11.1 Distribución de temperaturas por convección y radiación en estado estable.
....................................................................................................................................... 78
CAPÍTULO IV ..................................................................................................................... 81
CONCLUSIONES Y RECOMENDACIONES ................................................................. 81
4.1 CONCLUSIONES...................................................................................................... 81
4.2 RECOMENDACIONES............................................................................................ 82
BIBLIOGRAFÍA.................................................................................................................. 83
ANEXOS ............................................................................................................................... 85

xi
ÍNDICE DE TABLAS
Tabla 1.1: Valores tradicionales del coeficiente de transferencia de calor por
convección.................................................................................................................... 7
Tabla 1.2: Emisividad de algunos materiales a 300 °K. ............................................... 9
Tabla 1.3: Tipos de radiación. .................................................................................... 16
Tabla 3.1: Parámetros del prototipo de secador industrial. ........................................ 43
Tabla 3.2: Propiedades del aire a T = 316°𝐾............................................................. 44
Tabla 3.4: Propiedades del aire a T = 298°𝐾............................................................. 51
Tabla 3.5: Datos de conducción de materiales. .......................................................... 56
Tabla 3.6 Temperaturas de media cámara.................................................................. 63
Tabla 3.7 Temperaturas en la parte inferior de la cámara. ......................................... 65
Tabla 3.8: Distribución de temperaturas en la mitad de la cámara. ........................... 69
Tabla 3.9: Distribución de temperaturas en la cara inferior de la cámara. ................. 70
Tabla 3.10: Temperatura media de la cámara a 420s. ................................................ 72
Tabla 3.11: Temperaturas en la cara inferior de la cámara a 420s. ............................ 73
Tabla 3.12: Temperaturas en la mitad de la cámara a 688s. ...................................... 75
Tabla 3.13: Temperaturas en la cara inferior de la cámara a 688s. ............................ 76
Tabla 3.14: Análisis experimental y simulado de las zonas estudiadas. .................... 77
Tabla 3.15: Error absoluto y relativo en la zona central. ........................................... 78
Tabla 3.16: Error absoluto y relativo en la cara inferior. ........................................... 78
Tabla 3.17: Distribución de temperaturas en la mitad de la cámara en estado estable.
.................................................................................................................................... 79
Tabla 3.16: Distribución de temperaturas en la cara inferior de cámara en estado
estable......................................................................................................................... 80

xii
ÍNDICE DE FIGURAS
Figura 1.1: Conducción. ............................................................................................... 4
Figura 1.2: Convección forzada. .................................................................................. 6
Figura 1.3: Convección natural. ................................................................................... 6
Figura 1.4: Radiación. ................................................................................................. 8
Figura 1.5: Niveles de ondas electromagnéticas. ....................................................... 15
Figura 1.6: Categorías que abarca CFD. .................................................................... 19
Figura 1.7: Secador por convección de bandejas. ...................................................... 22
Figura 1.8: Horno secador infrarrojo. ........................................................................ 23
Figura 1.9: Secador de túnel por infrarrojos. ............................................................. 24
Figura 2.1: Dimensiones internas de la cámara de secado. ........................................ 28
Figura 2.2: Aislante de la cámara de secado. ............................................................. 29
Figura 2.3: Carcasa de la cámara de secado. .............................................................. 30
Figura 2.4: Verificación del modelo. ......................................................................... 30
Figura 2.5: Tapa del modelo. ..................................................................................... 31
Figura 2.6: Creación de wizard. ................................................................................. 32
Figura 2.7: Tipo de análisis. ....................................................................................... 32
Figura 2.8: Parámetros de análisis. ............................................................................ 33
Figura 2.9: Selección de fluido. ................................................................................. 33
Figura 2.10: Condiciones iniciales. ............................................................................ 34
Figura 2.11: Dominio de cálculo. ............................................................................... 35
Figura 2.12: Selección de materiales. ........................................................................ 35
Figura 2.13: Condición de frontera (velocidad del aire). ........................................... 36
Figura 2.14: Condición de frontera (presión)............................................................. 37
Figura 2.15: Adecuación del calor por convección. ................................................... 37
Figura 2.16: Superficies radiantes. ............................................................................. 38
Figura 2.17: Adecuación del calor por radiación. ...................................................... 39
Figura 2.18: Controlador de temperatura. .................................................................. 40
Figura 2.19: Características del controlador. ............................................................. 40
Figura 2.20: Termocupla. ........................................................................................... 41
Figura 2.21: Prototipo del secador. ............................................................................ 41
Figura 3.1: Dimensiones de la cámara de secado. ..................................................... 54
Figura 3.2: Dimensiones de la carcasa del secador. ................................................... 55
Figura 3.3: Diagrama del balance de energía. ............................................................ 61
Figura 3.4: Volumen de la cámara. ............................................................................ 62
xiii
Figura 3.5: Ubicación central de la termocupla. ........................................................ 63
Figura 3.6: Ubicación inferior de la termocupla. ....................................................... 65
Figura 3.7: Convección y radiación para un tiempo de 180s. .................................... 67
Figura 3.8: Distribución de temperaturas para un tiempo de 180s. ........................... 67
Figura 3.9: Temperaturas del aire para un tiempo de 180s. ....................................... 68
Figura 3.10: Temperaturas en la mitad de la cámara para t=180s. ............................ 68
Figura 3.11: Temperaturas en la cara inferior de la cámara para t=180s. .................. 69
Figura 3.12: Convección y radiación para un tiempo de 240s. .................................. 70
Figura 3.13: Distribución de temperaturas para un tiempo de 240s. ......................... 71
Figura 3.14: Temperaturas del aire para un tiempo de 420s. ..................................... 71
Figura 3.15: Temperatura en la mitad de la cámara para t=240s. .............................. 72
Figura 3.16: Temperaturas en la cara inferior de la cámara para t=420s. .................. 73
Figura 3.17: Convección y radiación para un tiempo de 688s. .................................. 74
Figura 3.18: Distribución de temperaturas para un tiempo de 688s. ......................... 74
Figura 3.19: Temperaturas del aire para un tiempo de 688s. ..................................... 75
Figura 3.20: Temperaturas en la mitad de la cámara para t=688s. ............................ 75
Figura 3.21: Temperaturas en la cara inferior de la cámara para t=688s. .................. 76
Figura 3.22: Temperaturas en la mitad de la cámara para estado estable. ................. 79
Figura 3.23: Temperaturas en la cara inferior de la cámara para estado estable. ...... 79

xiv
 RESUMEN

El presente trabajo de investigación se centra en las distribuciones de temperaturas

en el interior de la cámara de secado con una fuente calorífica de radiación que en
este caso se utilizan cuatro lámparas infrarrojas de 250 W cada una, para la obtención
de dichas temperaturas se realizan análisis experimentales en el prototipo de la
cámara y también con un modelado del secador se obtuvo simulaciones realizadas en
el software, en estado transitorio con intervalos de tiempos de 180, 240 y 688
segundos.

La metodología planteada para el estudio fue la simulación en software verificando

primero el modelado de la cámara, luego se procede con los pasos para el ingreso de
datos como son: parámetros para el tipo de análisis, tipo de fluido (aire), condiciones
iniciales, selección de materiales, condiciones de frontera, cantidad de calor y
superficies radiantes, también mediante cálculos de los parámetros de transferencia
de calor en el interior de la cámara por convección natural, forzada y radiación, con
sus combinaciones alcanzan los 771 W, además se obtuvo la variación de energía
almacenada para el aire de 15,65 W. Aplicando balance de energía se determinan las
perdidas por convección al exterior que alcanzan los 25 W.

Los resultados combinados de los métodos de transferencia de calor ayudan al

cálculo de la eficiencia de la cámara en vacío (sin carga) y considerando solamente el
calentamiento del aire al interior nos da 6,4 por ciento. En los intervalos de tiempos
analizados experimentales y simulados se obtienen las distribuciones de temperaturas
en la zona media e inferior de la cámara alcanzando valores máximos de 47 y 58,75
grados centígrados respectivamente.

xv
 ABSTRACT

The present research work is focused on the temperature distributions inside the
drying chamber with a radiation heat source that in this case four infrared lamps of
250 W each are used. In order to obtain these temperatures, experimental analyses
are carried out in the prototype of the chamber and also with a modeling of the dryer
simulations were obtained in the software, in a transient state with time intervals of
180, 240 and 688 seconds.

The methodology proposed for the study was software simulation, first verifying the
modeling of the chamber, then proceeding with the steps for data entry such as:
parameters for the type of analysis, type of fluid (air), initial conditions, selection of
materials, boundary conditions, amount of heat and radiating surfaces, also through
calculations of the parameters of heat transfer inside the chamber by natural
convection, forced and radiation, with their combinations reaching 771 W, also
obtained the variation of stored energy for the air of 15.65 W. Applying energy
balance, the losses by convection to the outside are determined, which reach 25 W.

The combined results of the heat transfer methods help to calculate the efficiency of
the chamber in vacuum (without load) and considering only the heating of the air
inside gives us 6.4 percent. In the experimental and simulated time intervals
analyzed, the temperature distributions in the middle and lower part of the chamber
reach maximum values of 47 and 58.75 degrees Celsius, respectively.

xvi
 CAPÍTULO I

MARCO TEÓRICO

1.1 Antecedentes Investigativos.

En el presente trabajo necesitamos de diferentes antecedentes investigativos acerca
de los diferentes tipos de secadores industriales o semiindustriales construidos,
simulaciones térmicas y computacionales para lo cual detallamos las siguientes
investigaciones importantes para el desarrollo de esta tesis de simulación térmica de
un prototipo de cámara de secado referente a la tesis desarrollada por el señor Juan
Carlos Espinoza Freire, en el año 2013 bajo el tema: “ESTUDIO DEL PROCESO
DE SECADO CON LÁMPARAS INFRARROJAS Y SU INFLUENCIA SOBRE
LA DETERMINACIÓN DEL GRADO DE HUMEDAD EN CUBETAS DE
CARTÓN DE HUEVOS PRODUCIDAS EN LA FÁBRICA AVÍCOLA UBICADA
EN SAMANGA”, en la cual determinó y retiró humedad del cartón de las cubetas
para huevos, con lámparas infrarrojas como fuente de calor. Dichas cubetas de cartón
fueron sometidas un proceso de calentamiento por medio de lámparas infrarrojas
(radiación infrarroja) en un prototipo de secador con potencia calorífica de las
diferentes lámparas de 150W, 250W, hasta llegar a una potencia calorífica de 950W
logrando una temperatura interior de la cámara de 65 a 70°C, alcanzando con la
máxima carga tiempos de secado de entre 25 a 40 minutos. [1]

En la Universidad Tecnológica Equinoccial en la Facultad de Ciencias de la

Ingeniería Carrera de Ingeniería Automotriz, la tesis realizada por Guillermo Ortega
en 2016 bajo el tema “DISEÑO, COSTRUCCION E IMPLEMENTACION DE UN
PROTOTIPO DE UN HORNO DE SECADO DE PINTURA AUTOMOTRIZ Y
PRUEBAS DE PINTURA EN LAS PROBETAS AL FINAL DEL PROCESO”,
determinó y conllevó que un prototipo de horno de pintura automotriz que funciona
por energía infrarroja necesita de un análisis exhaustivo de los factores que
intervienen en el funcionamiento y constitución de la transmisión de energía, para lo
cual se debe determinar el número adecuado de lámparas infrarrojas para
introducirlas en el interior del prototipo, además se realizó un completo análisis

1
termodinámico. También se determinó el tiempo de encendido de las lámparas para
un secado uniforme de probetas, el aislante térmico con su espesor adecuado para
evitar fugas de calor hacia el exterior del horno. [2]

En el Centro de Innovación Aplicada en Tecnologías Competitivas (CIATEC), la

tesis realizada por Fernando Bravo en Diciembre del 2015 en “SIMULACION Y
ESCALAMIENTO INDUSTRIAL DE UN SECADOR INFRARROJO PARA EL
ACABADO DEL CUERO”, se determina como objetivo de la investigación,
desarrollar de manera virtual un modelo de maquina con la función de realizar el
proceso de sacado del cuero con un túnel de secado que tiene dentro paneles
infrarrojos con la finalidad de secar el agua del acabado del cuero, por lo cual se usó
un modelo en 3D de los paneles y se realizó un análisis en el software ANSYS CFX
conociendo las temperaturas en el transcurso del proceso y la cantidad de agua
evaporada, se empleó además el volumen finito es el método numérico utilizado, los
resultados se obtuvieron por medio de comparaciones de datos de evaporación y
perfiles térmicos experimentales. La elaboración de conjuntos de modelos virtuales
experimentales es la mejor opción para una elección adecuada del modelo escalado a
ser determinado y puesto en práctica. [3]

2
1.2 Objetivos
1.2.1 Objetivo general
➢ Analizar térmicamente un prototipo de cámara de secado estática por
convección y radiación de 1000 W para uso semiindustrial.
1.2.2 Objetivos específicos
➢ Realizar el balance de energía en la cámara de secado.
➢ Determinar la eficiencia térmica del proceso.
➢ Simular mediante software la distribución de temperatura al interior de la
cámara de secado en estado estable y transitorio.
➢ Comparación entre los resultados experimentales y los obtenidos de la
simulación.

3
1.3 Fundamentación teórica
1.3.1 Transferencia de calor
Es la ciencia que nos da a conocer el intercambio de energía entre sistemas, cuerpos,
elementos, materiales, dicho intercambio nos da como resultado de las diferencias
entre temperaturas. Teniendo en cuenta que al generar el calor puede haber perdidas
de calor entre superficies y fluidos en contacto respectivamente. [4]

Esta ciencia encierra tres mecanismos o formas diferentes de pasar la energía (calor)
de un cuerpo emisor al receptor. Estas son las siguientes: conducción, convección y
radiación. [5]

En el presente estudio los parámetros de transferencia de calor a tener en cuenta son:

convección y radiación, dados a conocer en el transcurso de la investigación.

1.3.1.1 Transferencia de calor por conducción

La transferencia de energía (calor) por conducción tiene dos procesos de realización.

El primero es interacción molecular, pues la partícula aumenta el movimiento

llegando a una temperatura de nivel alto el cual comparte o conduce calor a las
moléculas próximas que están a un nivel bajo de temperatura, y el segundo,
trasferencia de energía (calor) por conducción por medio de electrones libres, este
proceso es uno de los procesos importantes ya que se utiliza en solidos metálicos
puros que son los mejores conductores de calor. [6]

Figura 1.1: Conducción. [10]

4
La ecuación básica que nos permite calcular la transferencia de calor por conducción
es la siguiente:

dT
Q̇cond = −k. A. ( ) Ec.1.1
dx

Donde:

Q̇cond = Flujo de calor por conducción [W]

k = Conductividad térmica del material [W⁄m ∗ °C]

A = Area de la superficie [m2 ]

dT
(dx ) = Gradiente de temperatura [°C/m]

1.3.1.2 Transferencia de calor por convección

La transferencia de energía (calor) por convección es la transmisión calorífica de un
fluido cuando este se mezcla con otra puede ser una superficie que no se produce por
contacto directo, a causa de movimientos de masa del mismo. La energía puede
transportarse de un punto a otro del espacio por el desplazamiento del mismo fluido,
el movimiento del fluido puede producirse a causas mecánicas externas puede ser por
un ventilador o una bomba la se le denomina convección forzada y la convección
natural o libre cuando un fluido circula a causa de densidades diferentes en el área
del fluido que resulta de la variación de temperaturas. [6]

La ecuación de newton de la rapidez es la que denomina la rapidez de la

transferencia de calor por convección, esta ecuación es:

Q̇conv = h. As . (Ts − Tf ) Ec.1.2

Donde:

Q̇conv = Flujo de calor por convección [W]

h = Coeficiente de transferencia de calor por convección [W⁄m2 ∗ °C]

5
As = Área de la superficie [m2 ]

Ts = Temperatura de la superficie [°C]

Tf = Temperatura del fluido alejado de la superficie [°C]

La transferencia de calor por convección se divide en dos aspectos que debemos

tomar en cuenta que es la convección forzada y la convección natural:

a) Convección forzada, esta ocurre cuando el fluido en forzado a circular por la

superficie en contacto esto sucede con la ayuda de diferentes mecanismos. [8]

Figura 1.2: Convección forzada.

b) Convección natural (libre), esta ocurre cuando el fluido circular libremente

sobre una superficie esto se debe a la diferencia de densidades del fluido o gas por la
diferencia de temperaturas. [8]

Figura 1.3: Convección natural.

6
 Tabla 1.1: Valores tradicionales del coeficiente de transferencia de calor por
convección. [8]

Tipo de convección h, [𝐖⁄ 𝟐 ]

𝐦 ∗ °𝐂

Convección libre de gases 2-5

Convección libre de líquidos 10-1000

Convección forzada de gases 25-250

Convección forzada de líquidos 50-20000

Ebullición y condensación 2500-100000

La resistencia de la transferencia de calor por convección es expresada por la

siguiente ecuación:

1
R conv = [°C⁄W] Ec.1.3
hc . A

Donde:

R conv = Resistencia Térmica por convección.

hc = Coeficiente de transferencia de calor por convección [W⁄m2 ∗ °C]

A = Área de la superficie.

1.3.1.3 Transferencia de calor por radiación

La transferencia de calor por radiación se diferencia bastante de la conducción y
convección ya que estas deben estar en contacto con superficies y con la intervención
de las diferencias de temperaturas en el entorno de contacto o ambiente de flujo, la
radiación se refiere en parte a un fenómeno relativo a cuerpos luminosos, no pretende
de la influencia de un medio y la energía calorífica puede ser transferida atravesando
el vacío absoluto. [5]

7
 Figura 1.4: Radiación.

La ley de Stefan-Boltzmann para cuerpos negros nos dice que la máxima razón de la
radiación que se puede exponer desde un medio o superficie a una temperatura
termodinámica Ts [K o R] se expresa en la siguiente formula. [8]

Q̇ emitida = σ. As . Ts4 Ec.1.4

La transferencia de calor por radiación para cuerpos reales (superficies) es más baja
que la expuesta por un cuerpo negro a las mismas condiciones (temperatura), donde
interviene la emisividad y área superficial de materiales expuestos en la (Tabla 1.2).
[8]

La ecuación es expresada de la siguiente manera:

Q̇rad = ε. σ. As . (Ts4 − Talred

4
) Ec.1.5

Donde:

Q̇rad = Flujo de calor por radiación [W]

ε = Emisividad de la superficie (material). Ver tabla 1.2

σ = Constante de Stefan-Boltzmann (5,67X10−8 [W⁄m2 ∗ K 4 ])

As = Área de la superficie [m2 ]

Ts4 = Temperatura de la superficie [°K]

4
Talred = Temperatura alrededor [°K]

8
 Tabla 1.2: Emisividad de algunos materiales a 300 °K. [8]

Material Emisividad

Hoja de aluminio 0,07

Aluminio anodizado 0,82

Cobre pulido 0,03

Oro pulido 0,03

Plata pulida 0,02

Acero inoxidable pulido 0,17

Pintura negra 0,98

Pintura blanca 0,90

Papel blanco 0,92-0,97

Pavimento de asfalto 0,85-0,93

Ladrillo rojo 0,93-0,96

Piel humana 0,95

Madera 0,82-0,92

Suelo 0,93-0,96

Agua 0,96

Vegetación 0,92-0,96

La resistencia de la transferencia de calor por radiación es expresada por la siguiente

ecuación:

1
R rad = [°C⁄W] Ec.1.6
hrad . As

9
Donde:

R rad = Resistencia térmica por radiación.

hrad = Coeficiente de transferencia de calor por radiación [W⁄m2 ∗ °C]

As = Área de la superficie [m2 ]

Para el cálculo del coeficiente de transferencia de calor por radiación se aplica la

siguiente ecuación:

2
hr = ε. σ. (Ts + Talred ). (Ts2 − Talred ) Ec.1.7

1.3.2 Calor.
Es el movimiento brusco de vibración molecular esto sucede cuando calentamos un
cuerpo o sustancia, las moléculas se mueven muy rápido produciendo energía (calor),
por el contrario, si enfriamos una sustancia o cuerpo el movimiento de las moléculas
se mueven lentamente hasta detenerse generando temperaturas bajas. [8]

Todos los cuerpos tienen calor en distintas partes con mayor o menor energía
calorífica, además podemos decir por los conocimientos adquiridos una frase que
dice lo siguiente “la energía no se crea ni se destruye solo se transforma”, el calor se
expresa con la siguiente ecuación. [8]

Q̇ = m. Cp. (T2 − T1 ) Ec.1.8

Donde:

m = Cantidad de masa de aire a ser calentado [kg]

kJ
Cp = Calor específico del fluido (aire) [ ⁄kg ∗ °C]

T2 = Temperatura caliente del fluido [°C]

T1 = Temperatura fría del fluido [°C]

10
1.3.3 Balance de energía para sistemas de flujo estacionario
El flujo de masa de mecanismos o sistemas de ingeniería es hacia adentro y hacia
afuera, este es considerado como volumen de control y se analiza en condición
estacionaria. Estacionario para un mejor entendimiento es sencillamente ningún
cambio con el tiempo en una ubicación específica, lo opuesto es considerado como
transitorio o no estacionario. Para un proceso de flujo estacionario la cantidad de
volumen que accede en un sistema en diferentes formas ya sea como: energía,
trabajo, transferencia de masa entre otros aspectos, tiene que ser la misma a la
cantidad de energía que sale del sistema. [8]

1.3.4 Convección forzada interna

Se considera para la convección de flujo interno parámetros hidrodinámicos como la
velocidad, también se considera aplicar un balance general de energía para así
establecer las variaciones de temperatura en dirección del flujo del fluido. Esto nos
ayuda para la distribución de flujo interno que nos presenta una geometría adecuada
para enfriar y calentar fluidos usadas en la tecnología de proceso químico, control de
ambiente y conservación de la energía. [9]

1.3.5 Número de Reynols

Es un número adimensional que depende de la velocidad de cierto fluido y sale por
cierta área de determinada configuración (diámetro de tubería), formándose el chorro
del flujo y se divide en dos tipos de flujo: laminar y turbulento. [10]

Re < 2000 → Flujo laminar

2000 ≤ Re ≤ 4000 → Flujo de transición

Re > 2000 → Flujo turbulento

V∗x
Re = [adimensional] Ec.1.9
v

Donde:

m
V = Velocidad del aire [ s ]

11
x = distancia entre las paredes (ingreso del aire) [m]

m2
v = viscosidad cinemática del aire [ ]
s

1.3.6 Número de Nusselt

Es el resultado experimental para la transferencia de calor sobre una placa plana con
una superficie que desprende calor constante, donde Nux en correlación de
convección para flujo turbulento es el valor local del número en una posición x
medida desde el inicio de la placa y se la expresa de la siguiente manera. [9]

Nux = 0.0296Re4/5 . Pr1/3 Ec.1.10

Donde:

Nux = Número de Nusselt local (flujo turbulento).

Re = Número de Reynolds.

Pr = Número de Prandtl.

De una forma más simplificada y para mejor uso del número de Nusselt en la
investigación nos queda de la siguiente manera:

Lc
Nu = hc . [adim] Ec.1.11
k

Donde:

hc = Coeficiente de transferencia de calor por convección forzada [W⁄m2 ∗ °C]

Lc = Longitud característica [m]

k = Constante de conductividad térmica [W⁄m ∗ °C]

1.3.7 Longitud característica

“Es la distancia que recorre una partícula dentro del volumen de un fluido con
respecto al área de la superficie en la capa limite utilizada”. [4]

Para calcular la longitud característica tenemos la siguiente ecuación:

12
 V
Lc = A [m] Ec.1.12
s

Donde:

V = Volumen [m3]

As = Área de la superficie [m2]

1.3.8 Convección libre interna o convección interna

1.3.8.1 Cavidades rectangulares
Disponemos de estudios experimentales y teóricos en forma extensa para el mejor
conocimiento de cavidad rectangular en la transferencia de calor por convección libre
o natural que nos dice, “dos de las paredes opuestas se mantienen a diferentes
temperaturas, mientras que las paredes restantes se aíslan de los alrededores”. [9]

También se da a conocer que existe entre las superficies caliente y fría un ángulo de
inclinación τ que puede variar que se detalla a continuación:

• 0° cavidad horizontal con calentamiento inferior.

• 90° cavidad vertical con calentamiento lateral.
• 180° cavidad horizontal con calentamiento superior. [9]

La siguiente ecuación expresa el flujo de calor a través de la cavidad.

Q = h(T1 − T2 ) Ec.1.13

1.3.8.2 Recintos cerrados

La transmisión de energía calorífica por convección a través de espacios cerrados
tiene un gran interés para el estudio ya que se complica porque el fluido encerrado no
está inmóvil. En un recinto vertical cerrado el fluido junto a la superficie o cara más
caliente sube y la superficie junto a la más fría baja. [8]

1.3.9 Número de Grashof

Este número físicamente puede interpretarse como adimensional que determina la
relación entre fuerzas de flotabilidad y las fuerzas viscosas en la corriente de flujo de
convección natural. El número de Grashof crítico tiene un valor aproximado de

13
4X108 para convección libre del aire que actúa sobre una placa plana vertical. Para
diferentes fluidos y también para diferentes niveles de turbulencia del ambiente se
pueden observar valor de entre 108 y 109. [4]

1.3.10 Número de Rayleigh

El número de Rayleigh es el producto del número Grashof y el número de Prandtl y
se expresa así:

Ra = Gr ∗ Pr [adim] Ec.1.14

Además, tenemos una ecuación para el cálculo del número de Rayleigh expresada a
continuación que facilitó el estudio en la cual nos dice que se aplica para recintos
cerrados. [8]

gβ(T1− T2 )L3c
Ra = . Pr [adim] Ec.1.15
v2

Donde:

Ra = Número de Rayleigh

g = Aceleración estándar de la gravedad [g=9.80665 m/s2]

β = coeficiente de expansión térmica [°K −1 ]

(T1− T2 ) = Diferencia de temperaturas

Lc = longitud característica [m]

m2
v = viscosidad cinemática [ ]
s

Pr = Número de Prandtl [adimensional]

1.3.10.1 Temperatura de película (promedio)

La temperatura de película o promedio se la calcula mediante las propiedades del
fluido dirigiendo a tablas de propiedades y se aplica la siguiente formula:

(Ts −T∞ )
Tf = Ec.1.16
2

14
Donde:

Tf = temperatuta de la película o promedio [°C]

Ts = temperatura de la superficie [°C]

T∞ = temperatura del ambiente [°C]

1.3.10.2 Coeficiente de expansión térmica

El coeficiente de expansión térmica se calcula con la temperatura de película o
promedio con la siguiente formula:

1
β=T Ec.1.17
f

1.3.11 Maneras de transferencia de calor por radiación.

1.3.11.1 Radiación térmica
Es la radiación más conocida como electromagnética formulada de un cuerpo como
consecuencia de una temperatura generada por el mismo. En la investigación
realizada a continuación se describe las diferentes generalidades de la radiación
térmica, sus características y materiales para dar a conocer su utilización e influencia
para la radiación. [4]

Figura 1.5: Niveles de ondas electromagnéticas. [11]

1.3.11.2 Tipos de radiación

La energía electromagnética por radiación se la clasifica en función a la longitud de
onda, cuando esta longitud se deduce como onda continua que se irradia a la
velocidad de la luz (3 X 108 m/s). Se considera a la radiación térmica es energía
15
electromagnética abarcando todo el espectro de longitudes de onda indicadas en
términos de micrómetro, μm (1 X 10-6 m). [12]

A continuación, la Tabla 1.3 muestra los tipos de radiación con una longitud de onda
aproximada:

Tabla 1.3: Tipos de radiación. [12]

Tipo De: Hasta:

Rayos Cósmicos 4 X 10-7 μm

Rayos Gamma 4 X 10-7 1.4 X 10-4 μm

Rayos X 1 X 10-5 2 X 10-2 μm

Rayos Ultravioletas 5 X 10-3 3.9 X 10-1 μm

Luz visible 3.9 X 10-1 7.8 X 10-1 μm

Radiación Solar 1 X 10-1 3.0 μm

Radiación Infrarroja 7.8 X 10-1 1 X 103 μm

Radiación Térmica 1 X 10-1 1 X 102 μm

Ondas Hertzianas 1 X 102 5 X 1010 μm

Ondas Radiofónicas 1 X 107 5 X 1010 μm

1.3.11.3 Radiación de un cuerpo negro

Todo cuerpo emite radiación en todas las direcciones a diferentes longitudes de onda
a una temperatura por encima de cero, esto depende de las condiciones del material y
la superficie, así como la temperatura del cuerpo. A un cuerpo negro se lo define
como un emisor y recibidor perfecto de radiación térmica incidente. [7]

1.3.11.4 Propiedades de la radiación térmica.

La radiación térmica es compartida o recibida a pequeñas distancias de la superficie
esto es para materiales grises. Por el contrario, otros materiales como el vidrio y agua

16
son materiales opacos a la radiación infrarroja y dejan que la radiación refleje sus
ondas a grandes profundidades antes que la radiación sea absorbida. Todo lo
mencionado se toma muy en cuenta para describir las propiedades relativas de la
radiación que son: la emisividad, la absortividad, la reflectividad y la transmisividad
de esos los diferentes materiales existentes puestos en estudio. [7]

1.3.11.4.1 La emisividad (ℇ)

“La emisividad de una superficie representa la razón entre la radiación emitida por la
superficie a una temperatura dada y la radiación emitida por un cuerpo negro a la
misma temperatura”. [7]

La emisividad es la distancia de la cual se encuentra cercana una superficie a un

cuerpo negro, por lo cual se da un valor de ℇ = 1, esta no es constante y depende del
material, la emisividad cambia por aspectos como la temperatura del cuerpo o
superficie, también influye la longitud de onda y la dirección en que se encuentra la
radiación expuesta. [7]

1.3.11.5 Espectrometría infrarroja

La espectrometría infrarroja o también llamada espectroscopia de infrarrojos es una
característica de la radiación electromagnética que se refiera a la absortividad. La
espectrometría infrarroja se refiere a que los enlaces químicos de las sustancias
tienen excitación energética potencial vibracional de las moléculas. [13]

Los usos y aplicaciones de la espectrometría infrarroja es muy amplia para

investigaciones científicas y para la industria, por lo que se destaca la manera rápida
y fiable de obtener datos de control y análisis, también con los instrumentos se puede
obtener datos de mediciones de campo ya que son pequeños y fácil de transportarlos,
con los avances tecnológicos computacionales y la favorable manipulación de
resultados se logra obtener medidas de gran precisión de muestras en una buscada
solución. A esto podemos mencionar la transformada de Fourier que es una técnica
de analizar el espectro infrarrojo y obtenerlo con una mayor rapidez. [13]

17
1.3.12 Radiación infrarroja
La radiación (energía) infrarroja es irradiada por medio de un fluido (aire) en forma
de ondas electromagnéticas como cualquier radiación térmica ya estudiada
anteriormente ya que es esta presente la temperatura absorbida por un cuerpo o
superficie cercana. La IR tiene mayor longitud de onda que la luz que podemos
apreciar (luz visible) y tiene menor frecuencia a esta, tiene también menor longitud
de onda a la radiación por microondas y mayor frecuencia que las microondas se
puede decir que son inversamente proporcionales, para lograr una mayor eficiencia
se debe decir que esto interviene directamente con la temperatura de la fuente que
debe lograr mayor temperatura para lograr mayor irradiación es este caso depende
una siempre de la otra. La radiación infrarroja se divide en tres regiones de espectro
infrarrojo. [14]

Las regiones cuyas longitudes de ondas son:

- Longitud de onda corta o cercana Entre: 0.78 a 2.5 μm

- Longitud de onda media Entre: 2.5 a 50 μm
- Longitud de onda larga o lejana Entre: 50 a 1000 μm

1.3.13 Dinámica de fluidos computacional (CFD)

La dinámica de fluidos computacional (CFD, siglas en ingles), es una forma
mecánica de fluidos a través de un software y se lo puede decir que es de manera
experimental, también podemos realizar estudios teóricos, estos dos estudios
dependen de las condiciones del problema a ser analizado. Es importante decir que
científicamente las dos técnicas son complementarias para la solución de un
problema a ser estudiado. Los datos experimentales que se analizan como los
resultados dan a conocer el entorno verdadero del fenómeno estudiado, el estudio
experimental implica un elevado costo. [15]

Existe también el estudio de modelación matemática (métodos numéricos), esta

técnica es menos costosa y eficiente obteniendo resultados en menor tiempo en
comparación al método experimental. Otro método es la solución analítica (métodos
analíticos), esta técnica es complicada ya que es difícil obtener la solución porque
contienen integrales, ecuaciones diferenciales, funciones especiales, entre otras. Por

18
el contrario, los métodos numéricos nos ayudan a obtener valores cercanos a la
solución del problema que queremos obtener. Los métodos numéricos abarcan tres
categorías de solución y dependen del fenómeno a ser estudiado y son: Método de
Diferencias Finitas (MDF), Método de Volumen Finito (MVF) y el Método de
Elementos Finitos (MEF) estos métodos nos ayudan a resolver ecuaciones varias
interferidas en la conservación de masa, conservación de energía, momento y
transporte de masa (especies químicas). [15]

Se establece que la dinámica de fluidos computacional (CFD) abarca a tres

categorías (Figura 1.6):

- Mecánica de fluidos, que enfatiza su estudio a fluidos en reposo y/o en

movimiento.
- Ciencias computacionales mediante un compilador.
- Matemáticas, ecuaciones diferenciales (ecuaciones gobernantes).

Figura 1.6: Categorías que abarca CFD. [15]

1.3.14 Estructura de Dinámica de Fluidos Computacional (CFD)

Por lo general, un problema necesita de una solución y en este caso (CFD) que
depende de una secuencia y estructura, esta secuencia fundamental consta de tres
partes principales que son: pre-proceso, procesamiento (solver) y postproceso. [15]

19
1.3.14.1 Pre-proceso de CFD
Esta parte es una de las más importantes ya que se caracteriza por detallar la entrada
de datos para su futura resolución y podemos describirlos para usuarios de softwares
comerciales de la siguiente manera. [15]

- Definir una geometría o dominio computacional conciso

- Generar una adecuada malla numérica
- Especificar los modelos que se pretenden representar
- Especificar las propiedades físicas y ambientales del medio a estudiar
- Establecer las condiciones iniciales y de frontera de la geometría

1.3.14.2 Procesamiento (Solver) de CFD

Esta segunda fase de dinámica de fluidos computacional (CFD) está constituida por
la parte media o central del software, el cual se encarga de resolver el modelo
matemático usando una computadora ya que está incluida en la solución del
problema, esta parte es cuando el controlador da la orden de ejecución y espera que,
de los resultados inicialmente propuestos, la obtención de los resultados puede variar
en el tiempo dependiendo de las restricciones y el modelo a resolver. [15]

1.3.14.3 Postproceso de CFD

En esta tercera fase de dinámica de fluidos computacional (CFD) podemos visualizar
la obtención de datos a gran escala descritos como en gráficos correspondientes para
el análisis de resultados y también ayuda a la mejor comprensión e interpretación del
fenómeno, lo cual nos conlleva a desarrollar conclusiones razonables de la solución
del fenómeno o problema de estudio, en esta etapa de postproceso también nos
permite representar y visualizar la malla, mapas, isolíneas, flujo de velocidad,
distribución de temperatura, etc. Para un mejor entendimiento además podemos
visualizar animaciones de diferente fenómeno del estudio a realzarse. [15]

1.3.15 Definición de secadores

Un secador es un equipo o máquina que ayuda al secado de materiales, sustancias,
etc. A través de la separación, eliminación, disipación de grados de humedad o

20
pequeñas fracciones de agua de los materiales mediante la evaporación. La corriente
de gas caliente es la fuente de calor de estos equipos o maquinas utilizadas para el
secado en general estas pueden funcionar de forma que el material puede permanecer
estático, en movimiento, en forma de fluido o de forma disuelta. [16]

1.3.16 Tipos de secadores infrarrojos

Para determinar los tipos de secadores iniciaremos dando a conocer que los
secadores tienen diferentes geometrías dependiendo de la forma del material a secar
o del uso industrial a utilizar, así como también la geometría depende del gas, fuente
de calor (tipo de calor) y los parámetros para el secado. [17]

1.3.16.1 Secador por convección o directo

Los secadores por convección puedes ser de dos métodos el de convección natural el
cual no se lo utiliza mucho en la industria del secado, en cambio el método de
convección forzada sobresale en este tipo de secadores ya que es más frecuente su
utilización esto es porque este utiliza aire caliente como fluido generado por la
quema de un combustible o por cualquier tipo de combustión con la ayuda de una
fuerza de aire por medio de ventiladores entre otros mecanismos que ayuden a
generas la convección forzada en general. [18]

Estos sistemas son generalmente sencillos y de fácil manipulación en el momento del

manejo, la industria agrícola es la que más utiliza este tipo de secadores por
convección ya que es de fácil acceso a adquirirlos por el bajo costo de
funcionamiento y porque no también diseñarlos y construirlos de acuerdo a la
necesidad de los usuarios, en la Figura 1.19, podemos observar un ejemplo de
secador por convección. [19]

21
 Figura 1.7: Secador por convección de bandejas. [19]

1.3.16.2 Secador por radiación

En estos secadores el principal agente que influye es la energía térmica infrarroja que
podemos observar en la Figura 1.10, como ya sabemos es la energía formada por
ondas electromagnéticas que se trasladan de un cuerpo con temperatura alta a uno de
menor temperatura. La transferencia de energía radiante a un material o cualquier
cuerpo húmedo nos ayuda con la evaporación del mismo ya que se produce
eléctricamente o por materiales refractarios calentados con fluidos provenientes de
algún gas, la energía utilizada es más costosa en comparación a los métodos de
secado por conducción y radiación. [19]

Los secadores infrarrojos utilizan lámparas infrarrojas, resistencias eléctricas o

materiales refractarios ya mencionados, suministrando eléctricamente la energía
radiante, su principal aplicación es en el cocido o horneado en la tecnología de
secado por capas de pintura y también estos secadores son utilizados para calentar
materiales con capas delgadas de masa. [16]

22
 Figura 1.8: Horno secador infrarrojo. [Autor]

1.3.16.3 Secador de infrarrojo a onda corta

Este equipo secador se desarrolló para servicio rápido de carrocerías o cualquier tipo
de taller de pintura con mayor capacidad productiva, se dispone de versiones
analógicas y digitales para un fácil manejo del funcionamiento para el secado de
cualquier proceso, pueden ser de distintos tipos que emiten ondas de: larga, media y
corta, con una diferencia principal de los tres que es en la penetración reflejada a la
película de pintura, la onda corta es la de mayor penetración que reduce el tiempo de
secado y reduce también el consumo de energía. [20]

1.3.16.4 Secador de túnel por infrarrojos

El principio de funcionamiento es la radiación por luz infrarroja cuyas ondas son
electromagnéticas y de espectro a la luz visible. La energía es llevada a través de la
onda y una vez impactada con la superficie esta se deposita en el objeto. Un
parámetro importante también es el ángulo de incidencia para su capacidad de
absorción, si el ángulo es más recto es decir 90 grados y el objeto es opaco la
cantidad de energía que absorbe es mayor. Estos hornos se utilizan para curar
recubrimientos sensibles al calor como plásticos, madera (MDF) entre otros
derivados, además reduce la utilización de hornos de curado continuos de
convección, reduce el tiempo de secado por lo que el costo de operación será menor.
Una de las desventajas principales de estos hornos es que cuando un objeto es más

23
brilloso la cantidad que refleja la luz es más alta por ende el calentamiento del objeto
se lo realizara lentamente. [21]

Figura 1.9: Secador de túnel por infrarrojos. [21]

24
 CAPÍTULO II

METODOLOGÍA
En este capítulo se estudia el comportamiento de la temperatura en la parte interior
de la cámara de secado, alcanzando los cálculos de los componentes de la
transferencia de calor por convección y radiación con diferentes flujos de calor
suministrado, además se crean los modelos en condiciones de estado estable y
transitorio para alcanzar los objetivos plantados en el proyecto, para lo cual se
procede a realizar una modelación computacional, concluyendo con un análisis del
proceder térmico computacional-experimental y encontrando la semejanza entre los
dos por medios gráficos y proximidades de temperaturas en las partes interiores de la
cámara a través de la comparación de varios parámetros térmicos expuestos en el
estudio.

2.1 NIVELES O TIPOS DE INVESTIGACIÓN

2.1.1 Descriptiva
Busca indicar las propiedades, características y mediciones de la distribución de
temperaturas dentro de la cámara de secado del proceso de los mecanismos de
transferencia de calor en estado estable y transitorio con datos tomados a través de
diferentes etapas por medio de la simulación y cálculos realizados. Estos son los
métodos más utilizados y confiables debido a que los valores se determinan en forma
cuantitativa.

2.1.2 Correlacional
Busca establecer relaciones entre el análisis térmico computacional con el
experimental (análisis global térmico dentro de la cámara de secado) con el fin de
llegar a lo planteado en los objetivos mediante cálculos matemáticos y
computacionales tienes o no relación.

25
2.1.3 Experimental
Experimentalmente se busca medir por medio del software la distribución de
temperatura en el interior de la cámara de secado a través de los diferentes ejes o
posiciones a tomar datos de la cámara simulada en el software.

2.1.4 Bibliográfica
Este método es uno de los más destacados o que tiene mayor incidencia en el
desarrollo de este trabajo experimental, ya que se basa en la recopilación de
información de fuentes bibliográficas provenientes de libros, tesis de grado, revistas,
artículos técnicos, manuales, etc. Estas fuentes bibliográficas fueron de gran
importancia en el desarrollo y entendimiento de esta investigación.

2.1.5 Exploratoria
Este tipo de investigación se la aplica también en la elaboración de trabajo
experimental ya que se analiza el uso y funcionamiento del software que se utiliza
comúnmente o de acuerdo al análisis térmico a realizar, con el fin de asegurar un
correcto análisis y una correcta medición de los parámetros encontrar, tratando de
disminuir en lo posible los errores presentes en los datos experimentales y
computacionales gracias a la utilización de softwares de simulación.

26
2.2 PROCESO DEL TRABAJO DE INVESTIGACIÓN

Búsqueda y elección de bibliografía

relacionada al tema de investigación

Selección de parámetros y configuración del

software para distribución de temperaturas.

Elaboración de partes y ensamblaje del modelo de la

cámara de secado estática para la simulación interna.

Validar y verificar el modelo para la

simulación.

Cálculos de los parámetros de transferencia

de calor en el interior de la cámara.

Simulación de la cámara de secado en el

software.

Análisis
experimental
vs simulación

Aprobación del sistema

simulado.

27
2.3 DIMENSIONES GEOMÉTRICAS Y CARACTERISTICAS DEL
MODELO
2.3.1 Dimensiones de la cámara de secado
El modelo de la cámara de secado mostrado en la Figura 2.1 es un prisma rectangular
que fue evaluada internamente por lo que a continuación se detallan las dimensiones
del modelo CAD:

- Largo = 900 mm
- Ancho = 600 mm
- Alto = 400 mm
- Espesor = 1,4 mm (material acero galvanizado)

Figura 2.1: Dimensiones internas de la cámara de secado. [Autor]

En la figura 2.1, además podemos observar cavidades, la cavidad del centro se

representa para la entrada de aire atmosférico y forzado simulando a un ventilador y
cuatro lamparás infrarrojas que alcanzarán los Watts necesarios para la simulación.

2.3.2 Aislante (lana de vidrio)

El aislante que se modeló va instalado en todas las caras de la cámara de secado
según el prototipo del secador que se simuló (no se realizó el cálculo del espesor),
tiene un espesor de 50mm separando la carcasa de la cámara, este aislante sirve para

28
reducir las pérdidas de calor por paredes, además vale recalcar que también tiene los
orificios ya detallados en el apartado anterior.

Figura 2.2: Aislante de la cámara de secado. [Autor]

Las propiedades del aislante para el ingreso de datos en el software se detallan en el

Anexo 5.

2.3.3 Carcasa del secador

la carcasa es la que cubre en su totalidad al aislante y la cámara de secado Figura 2.3,
el material es acero galvanizado de espesor 1.4 mm al igual que la cámara de secado
estos datos del material nos sirven para el ingreso de propiedades al software.

29
 Figura 2.3: Carcasa de la cámara de secado. [Autor]

2.4 VERIFICACIÓN DEL MODELO PARA LA SIMULACIÓN

En este punto validaremos el modelo CAD ya ensamblado, esto es necesario ya que
el software nos pide la verificación del volumen interno para que la simulación este
correcta ver Figura 2.4.

Figura 2.4: Verificación del modelo. [Autor]

30
En la Figura 2.4 podemos observar que la geometría a analizar esta correcta para un
análisis térmico interno de la cámara de secado para poder alcanzar los objetivos
deseados.

2.5 PASOS A SEGUIR PARA LA ENTRADA DE DATOS A LA

SIMULACIÓN
2.5.1 Chequeo del modelo
Una vez verificado el modelo se procede con el primer paso que es chequear el
volumen interno de la cámara de secado. Crear lids (tapas) ver figura 2.5 en caso de
ser necesario si la geometría del sistema no está totalmente cerrada para el análisis.
Una vez que nos indique que esta el proceso exitoso continuar con los siguientes
pasos.

Figura 2.5: Tapa del modelo. [Autor]

2.5.2 Principio de la simulación

Crear un nuevo wizard para la entrada de los primeros datos para las condiciones del
estudio en este caso es para cualquier simulación.

31
 Figura 2.6: Creación de wizard. [Autor]

2.5.4 Tipo de análisis

En esta ventana seleccionamos el tipo de análisis, en este caso es interno.

Figura 2.7: Tipo de análisis. [Autor]

32
2.5.5 Parámetros para el tipo de análisis
Para un análisis por convección y radiación seleccionamos los siguientes parámetros
en la pestaña de tipo de análisis como se muestra en la figura 2.8.

Para análisis de convección es necesario activar gravedad (Gravity), y para radiación

activar (Radiation), la temperatura ambiente se la deja por default.

Figura 2.8: Parámetros de análisis. [Autor]

2.5.6 Tipo de fluido

Seleccionar los fluidos a utilizar en la simulación en este caso es aire (AIR).

Figura 2.9: Selección de fluido. [Autor]

33
2.5.7 Condiciones iniciales
Esta es la última ventana del wizard de las condiciones iniciales y las dejamos por
defecto y damos clic en finalizar para que el software procese el modelo y los datos
que se ingresaron.

Figura 2.10: Condiciones iniciales. [Autor]

2.5.8 Dominio Computacional

Una vez finalizado el wizard del software seleccionamos el dominio computacional
que es en este caso de 3D, esto nos sirve para realizar el análisis interno o externo de
todo el volumen del modelo a analizar.

34
 Figura 2.11: Dominio de cálculo. [Autor]

2.5.9 Selección de materiales

Los materiales del modelo de la cámara de secado son distintos y se los va
seleccionando uno a uno, en la figura 2.13 se observa los diferentes materiales.

Figura 2.12: Selección de materiales. [Autor]

35
2.5.10 Condiciones de frontera
En este paso se ingresa la velocidad de aire en el orificio de la chimenea de la cámara
de secado como se muestra en la figura 2.14.

Se ingresa el dato de 0,4 m/s, además se marca la opción de flujo completamente

desarrollado (fully developed flow).

Figura 2.13: Condición de frontera (velocidad del aire). [Autor]

En las condiciones de frontera insertar también la presión que sería la presión

ambiental (environment pressure) y seleccionar la respectiva cara de entrada de
presión y dejarlo por default como se muestra en la figura 2.14.

36
 Figura 2.14: Condición de frontera (presión). [Autor]

2.5.11 Cantidad de calor

En recursos de calor (heat resources) se debe aplicar la cantidad de calor que emiten
cada una de las lámparas infrarrojas que es de 250 W que tiene cada lampara para el
caso de análisis de convección.

Figura 2.15: Adecuación del calor por convección. [Autor]

37
2.5.12 Superficies radiantes
En las superficies radiantes (radiative surfaces) se ingresan los materiales de las caras
que irradian o están expuestas a radiación en este caso son las lámparas infrarrojas y
el acero galvanizado de la cámara de secado.

El material de las lámparas infrarrojas es GLASS PYREX, y el material del acero

galvanizado se lo definió como IRON SMOOTH OXIDED.

Figura 2.16: Superficies radiantes. [Autor]

2.5.13 Calor por radiación

Luego en RADIATION RESOURCES se introduce la cantidad de calor por
radiación emitido por las lámparas infrarrojas que en este caso es de 250 que tiene
cada lampara y marcar en la opción de DIFFUSIVE ya que las lámparas emiten
radiación a todas las direcciones.

38
 Figura 2.17: Adecuación del calor por radiación. [Autor]

2.6 ANÁLISIS EXPERIMENTAL

La recolección de datos experimentales se realiza en los laboratorios de la Facultad
de Ingeniería Civil y Mecánica, donde se encuentra el prototipo de la cámara de
secado infrarroja que se está utilizando para este estudio, donde se podrá obtener
datos que comprueben el análisis térmico realizado al interior de la cámara de
secado.

Los instrumentos utilizados para la obtención de datos son cuaderno de apuntes,

termocuplas tipo k y cronometro.

2.6.1 Instrumentos para la obtención de información

2.6.1.1 Controlador digital de temperatura
Mediante un sensor de temperatura se puede observar datos en la pantalla del
controlador ya sean componentes físicos o químicos extraíbles depende de la
utilización, en este caso se toma datos de temperatura.

39
 Figura 2.18: Controlador de temperatura. [Autor]

En controlador a utilizar en el presente estudio para la toma de datos de temperatura

tiene las siguientes características mostradas en la figura 2.21.

Figura 2.19: Características del controlador. [Autor]

El rango de temperatura se encuentra dentro de los parámetros a buscar.

2.6.1.2 Termocupla tipo k

La termocupla se encarga de transmitir mediante un cable especial los datos de
temperaturas al controlador para que no salga del rango de estudio.

40
 Figura 2.20: Termocupla. [Autor]

2.6.1.3 Prototipo de cámara de secado para pruebas

Para la obtención de datos de temperatura de utilizó un prototipo de cámara de
sacado mediante lámparas infrarrojas.

Figura 2.21: Prototipo del secador. [Autor]

El prototipo de cámara de secado tiene las siguientes características:

- Dimensiones exteriores = Alto: 0,5m. Ancho: 1m. Fondo: 0,7m

- Dimensiones interiores = Alto: 0,4m. Ancho: 0,9m. Fondo: 0,6m
- Temperatura máxima = 95°C
- Potencia calorífica de las lámparas = 1000 W
- Ventilador = 100 CFM
- Velocidad de circulación del aire = 0,4 m/s

41
2.7 OBTENCIÓN DE LA INFORMACIÓN
Se realizaron mediciones en dos diferentes posiciones de la cámara, con intervalos de
un minuto se midió las temperaturas, se midió en el centro de la cámara, y en la cara
inferior de la cámara, obteniendo diferentes datos en los puntos mencionados, se
realizó una tabulación para observar el comportamiento de los resultados realizados
en la simulación y los experimentales.

42
 CAPÍTULO III

ANÁLISIS Y DISCUSIÓN DE RESULTADOS

En el presente capítulo se detalla los cálculos experimentales utilizando algunas de
las ecuaciones propuestas en el capítulo I para el respectivo análisis matemático,
además se realiza el análisis computacional por etapas para una simplificada manera
de solución de las variables de distribución de temperatura al interior de la cámara de
secado, con estos cálculos se procede a la comparación experimental con los
obtenidos en la simulación realizada en el software mediante tablas y gráficas
necesarias.

3.1 Obtención de datos

3.1.1 Parámetros de funcionamiento
Se toma datos experimentales del prototipo secador industrial estático por radiación
infrarroja que se encuentra en los laboratorios de la Facultad de Ingeniería Civil y
Mecánica en la Carrera de Ingeniería Mecánica que son necesarios para la
comparación y simulación del presente estudio, se podría decir que esta simulación
es un complemento para dicho secador lo cual se realiza un análisis interno de la
cámara de secado, con el desarrollo de la investigación alcanzar el objetivo general,
los parámetros extraídos son los siguientes:

Tabla 3.1: Parámetros del prototipo de secador industrial. [Autor]

Temperatura alcanzada por la cámara de secado 65 °C

Velocidad de circulación del aire 0,4 m/s

Temperatura ambiente o inicial 20 °C

Humedad relativa 60 %

Además, se toma un dato adicional que es el tiempo de calentamiento que alcanza la

cámara de secado del prototipo de secador a 950 W, de forma experimental realizada
por medio de controladores de temperatura en el secador industrial mencionado

43
anteriormente y se obtuvo que la temperatura de 65° se alcanza en un tiempo de 22
min ver Anexo 1.

3.1.2 Cálculos térmicos

Para estos cálculos se considera factores en los cuales el calor es unidimensional,
esto es para que el resultado sea aplicable y razonable, además se admite que el
análisis de la cámara está expuesto por el lado, inferior y paredes laterales a una
temperatura constante.

3.1.3 Cálculo del coeficiente de transferencia de calor para Convección Natural

Interior de la cámara
Es necesario calcular el parámetro de la temperatura promedio, la cual es medida
relacionada con la temperatura ambiente y la temperatura que alcanza la superficie.

Con los datos de la tabla 3.1 calculamos la temperatura promedio (Tprom) necesaria
para obtener las propiedades del aire ver Anexo 2, se utiliza la ecuación 1.16.

(𝑇𝑠 + 𝑇∞ )
𝑇𝑓 =
2

65°𝐶 + 20°𝐶
𝑇𝑓 =
2

𝑇𝑓 = 42,5 °𝐶 + 273,15 = 315,16 °𝐾 ≈ 316°𝐾

Tabla 3.2: Propiedades del aire a T = 316°𝐾. [Autor]

Viscosidad cinemática v = 18,1x10-6 m2/s

Numero de Prandtl Pr = 0,704

Conductividad térmica k = 0,0274 W/m°K

Difusividad térmica α = 0,246 m2/s

3.1.3.1 Parámetros térmicos

A más de las propiedades del aire necesitamos calcular parámetros térmicos para
hallar el coeficiente de convección natural.

44
Con la temperatura de película o promedio y la ecuación 1.17 se calcula el
coeficiente de expansión térmica.

1
𝛽=
𝑇𝑓

1
𝛽=
316 𝐾

𝛽 = 0,003164 𝐾 −1

Es necesario calcular el número de Grashof para el cálculo del coeficiente de

convección natural con la siguiente ecuación:

𝑔𝛽(𝑇𝑠 −𝑇∞ )𝐿3𝑐

𝐺𝑟 = . [𝑎𝑑𝑖𝑚] Ec.3.1
𝑣2

Donde:

𝐺𝑟 = Número de Grashof

𝑔 = Aceleración estándar de la gravedad [g=9.80665 m/s2]

𝛽 = coeficiente de expansión térmica [°𝐾 −1]

(𝑇𝑠 − 𝑇∞ ) = Diferencia de temperaturas

𝐿𝑐 = longitud característica [m]

𝑚2
𝑣 = viscosidad cinemática [ ]
𝑠

Para calcular la longitud característica aplicamos la ecuación 1.12

𝑉
𝐿𝑐 =
𝐴𝑠

Donde:

V = Volumen de la cámara de secado [m3]

𝐴𝑠 = Área de la superficie [m2]

El volumen de la cámara de secado nos da como dato en la Figura 2.4 es V= 0,214

m3 y el área de la superficie se toma 𝐴𝑠 = 0,9𝑚 ∗ 0,6𝑚 según las dimensiones de la

45
cámara de secado modelada para la simulación entonces reemplazando en la
ecuación 1.12 tenemos:

0,214
𝐿𝑐 =
0,54

𝐿𝑐 = 0,40 𝑚

Reemplazamos en la ecuación de Grashof.

(9.80665 𝑚/𝑠2)(0,003164 𝐾 −1 )(338𝐾 − 293𝐾)(0,40𝑚)3

𝐺𝑟 = . [𝑎𝑑𝑖𝑚]
(18,1𝑥10−6 𝑚2 /𝑠)2

𝐺𝑟 = 272,86𝑥106

A continuación, calculamos el número de Rayleigh que es el producto del número

Grashof y el número de Prandtl y se expresa con la ecuación 1.14:

𝑅𝑎 = 𝐺𝑟 ∗ 𝑃𝑟 [𝑎𝑑𝑖𝑚]

𝑅𝑎 = (272,86𝑥106 ) ∗ (0,704)

𝑅𝑎 = 192,1 𝑥106

Para calcular el número de Nusselt considerando flujo laminar y se tiene la siguiente

ecuación:

1
0,387𝑅𝑎6
𝑁𝑢𝐿 = 0,825 + 8 Ec.3.2
9 27
0,492
16
(1 + ( )
( 𝑃𝑟 ) )

Donde:

𝑁𝑢𝐿 = número de Nusselt (Flujo laminar)

𝑅𝑎 = número de Rayleigh

𝑃𝑟 = número de Prandtl

46
 2

1
0,387(192,1 𝑥106 )6
𝑁𝑢𝐿 = 0,825 + 8
9 27
0,492 16
(1 + ( ) )
( (0,704) )

𝑁𝑢𝐿 = 74,181

De la ecuación 1.11 se despeja el coeficiente de convección natural así:

𝑁𝑢𝐿 ∗ 𝑘
ℎ𝑐𝑛 =
𝐿𝑐

Además, con el coeficiente de conducción térmica del material tomado del Anexo 3,
tenemos 𝑘 = 0,0274 𝑊⁄𝑚 ∗ 𝐾 y se reemplaza en la ecuación anterior.

(74,181) ∗ (0,0274 𝑊⁄𝑚 ∗ 𝐾 )

ℎ𝑐𝑛 =
(0,40𝑚)

ℎ𝑐𝑛 = 5,081 [𝑊⁄𝑚2 ∗ 𝐾 ]

3.1.4 Calor generado por convección natural en el interior de la cámara

Una vez que se determinó el coeficiente de convección natural al interior de la
cámara se calcula el calor perdido por el aire en el interior de la misma, con la
siguiente ecuación:

𝑄̇𝑐𝑛 = ℎ̅ ∗ 𝐴 ∗ 𝛥𝑇 Ec.3.3

Donde:

𝑄̇𝑐𝑛 = flujo de calor generado por el aire (𝑊)

ℎ̅ = coeficiente medio de convección (𝑊⁄𝑚2 ∗ 𝐾 )

𝐴𝑠 = área superficial en el fondo de la cámara (𝑚2 )

𝛥𝑇 = variación de temperatura (𝐾)

47
 𝑄̇𝑎𝑖𝑟𝑒𝑐𝑛 = 6,775 (𝑊⁄𝑚2 ∗ 𝐾 ) ∗ (0,54𝑚2 ) ∗ (338 − 293)𝐾

𝑄̇𝑐𝑛 = 164,632 𝑊

3.1.5 Cálculo del coeficiente de transferencia de calor para Convección Forzada

Igualmente, con los datos requeridos en el apartado 3.1.3, calculamos la temperatura
promedio (Tprom).

𝑇𝑜 + 𝑇𝑓
𝑇𝑝𝑟𝑜𝑚 =
2

20°𝐶 + 65°𝐶
𝑇𝑝𝑟𝑜𝑚 =
2

𝑇𝑝𝑟𝑜𝑚 = 42,5 °𝐶 + 273,15 = 315,16 °𝐾 ≈ 316°𝐾

De la ecuación 1.9 se calcula el número de Reynolds.

𝑉∗𝑥
𝑅𝑒 = [𝑎𝑑𝑖𝑚]
𝑣

Para el cálculo y representación de x en la ecuación se tomará la distancia entre las

paredes laterales que es de 0,9m.

(0,4 𝑚⁄𝑠) ∗ (0,9 𝑚)

𝑅𝑒 = 2
(18,1𝑥10−6 𝑚 ⁄𝑠)

𝑅𝑒 = 19889,503 Flujo turbulento

El número de Reynolds es para flujo turbulento por lo que se utiliza la ecuación 1.10
que correlaciona el número de Nusselt.

𝑁𝑢 = 0.0296𝑅𝑒 4/5 . 𝑃𝑟 1/3

Aplicamos también la ecuación 1.11.

𝐿𝑐
𝑁𝑢 = ℎ𝑐
𝑘

48
Donde:

ℎ𝑐 = coeficiente de transferencia de calor por convección forzada [𝑊⁄𝑚2 ∗ °𝐶 ]

𝐿𝑐 = longitud característica [m]

𝑘 = constante de conductividad térmica [𝑊⁄𝑚 ∗ °𝐶 ]

Igualando la ecuación 1.10 y ecuación 1.11 obtenemos la siguiente igualdad.

𝐿𝑐
0.0296𝑅𝑒 4/5 . 𝑃𝑟 1/3 = ℎ𝑐
𝑘

Despejando hc y reemplazando datos en la igualdad se tiene:

4 1
[0.0296(98901,098)5 ∗ (0,704)3 ] [0,0274 𝑊⁄𝑚 ∗ °𝐶 ]
ℎ𝑐𝑓 =
0,40 𝑚

ℎ𝑐𝑓 = 4,955 𝑊⁄𝑚2 ∗ °𝐶

Con el dato del coeficiente de transferencia de calor por convección forzada se

calcula el flujo de calor generado por convección forzada con la siguiente ecuación:

𝑄̇𝑓𝑜𝑟𝑧 = 6 ∗ ℎ𝑐𝑓 ∗ 𝐴 ∗ 𝛥𝑇 Ec.3.4

Donde:

𝑄̇𝑓𝑜𝑟𝑧 = flujo de calor por convección forzada (W)

ℎ𝑐 = coeficiente de convección forzada (W/m2*°C)

𝐴= área de transferencia al fondo de la cámara (m2)

𝛥𝑇= variación de temperatura (°C)

49
 𝑄̇𝑓𝑜𝑟𝑧 = 6 (4,955 𝑊⁄𝑚2 ∗ °𝐶 ) ∗ 0,54 𝑚2 ∗ (65 − 20)°𝐶

𝑄̇𝑓𝑜𝑟𝑧 = 722,439 𝑊

3.1.6 Cálculo de calor generado por radiación al interior de la cámara

El calor por radiación es necesario calcular ya que la temperatura de operación en las
lámparas infrarrojas es alta en comparación a la temperatura ambiente.

Para este cálculo aplicamos la ecuación 1.5 con una emisividad del material para el
acero galvanizado 𝜀 = 0,28

𝑄̇𝑟𝑎𝑑 = 𝜀. 𝜎. 𝐴𝑠 . (𝑇𝑠4 − 𝑇𝑎𝑙𝑟𝑒𝑑

4 )

𝑄̇𝑟𝑎𝑑 = Flujo de calor disipado por radiación [𝑊]

𝜀 = Emisividad de la superficie (material)

𝜎 = Constante de Stefan-Boltzmann (5,67𝑋10−8 [𝑊⁄𝑚2 ∗ 𝐾 4 ])

𝐴𝑠 = Área de la superficie [𝑚2 ]

𝑇𝑠4 = Temperatura de la superficie [°𝐾]

4
𝑇𝑎𝑙𝑟𝑒𝑑 = Temperatura alrededor [°𝐾]

𝑄̇𝑟𝑎𝑑 = (0,28) (5,67𝑋10−8 [𝑊⁄𝑚2 ∗ 𝐾 4 ]) (0,9𝑚 ∗ 0,6𝑚). (3384 − 2934 )

𝑄̇𝑟𝑎𝑑 = 48,70 𝑊

3.1.6.1 Resumen de calores generados al interior de la cámara

Tabla 3.3: Resumen de calores. [Autor]

Convección natural + radiación 213,332 W

Convección forzada + radiación 771,139 W

50
Considerando que el calor de entrada es de 1000 W (producida por 4 focos), el 77 %
(771,13 W) se convierten en transferencia de calor por convección forzada más
radiación al interior de la cámara.

3.2 CÁLCULO DE PÉRDIDAS POR CONVECCIÓN AL EXTERIOR DE LA

CÁMARA
3.2.1 Análisis del coeficiente de convección natural al exterior del secador
Los parámetros a tomar en cuenta para el análisis son los siguientes:

- Relaciones de los efectos de radiación despreciables.

- Condiciones de entorno en estado estable.

Se toma los siguientes datos como la temperatura de la superficie al exterior del

secador y se aplica la ecuación 1.16 de la temperatura promedio con los siguientes
datos:

Ts = 30°C

𝑇∞ = 20°C

(𝑇𝑠 − 𝑇∞ )
𝑇𝑓 =
2

(30 + 20)
𝑇𝑓 =
2

𝑇𝑓 = 25°𝐶 = 298 °𝐾

Para obtener las propiedades del aire a la temperatura de 298 °𝐾 se toma los datos
del Anexo 2 e interpolando tenemos los siguientes valores:

Tabla 3.4: Propiedades del aire a T = 298°𝐾. [Autor]

Viscosidad cinemática v = 16,59x10-6 m2/s

Numero de Prandtl Pr = 0,709

Conductividad térmica k = 0,026 W/m°K

51
 Difusividad térmica α = 0,218x10-4 m2/s

Se calcula el coeficiente de expansión térmica con la ecuación 1.17 así

1
𝛽=
𝑇𝑓

1
𝛽=
298 𝐾

𝛽 = 0,003355 𝐾 −1

Con los datos de las propiedades del aire calculamos el número de Grashof con la
ecuación 3.5 y con el dato de la longitud característica ya calculada que es 𝐿𝑐 =
0,40 𝑚

𝑔𝛽(𝑇𝑠 −𝑇∞ )𝐿3𝑐

𝐺𝑟 = . [𝑎𝑑𝑖𝑚] Ec.3.5
𝑣2

(9.80665 𝑚/𝑠2)(0,003355 𝐾 −1 )(303𝐾 − 293𝐾)(0,40𝑚)3

𝐺𝑟 = . [𝑎𝑑𝑖𝑚]
(16,59𝑥10−6 𝑚2 /𝑠)2

𝐺𝑟 = 76,506𝑥106

A continuación, se calcula el número de Rayleigh con la ecuación 1.14

𝑅𝑎 = 𝐺𝑟 ∗ 𝑃𝑟 [𝑎𝑑𝑖𝑚]

𝑅𝑎 = (76,506𝑥106 ) ∗ (0,709)

𝑅𝑎 = 54,243 𝑥106

Cálculo del número de Nusselt considerando flujo laminar con la ecuación 3.6.

1
0,387𝑅𝑎6 Ec.3.6
𝑁𝑢𝐿 = 0,825 + 8
9 27
0,492 16
(1 + ( )
( 𝑃𝑟 ) )
52
 2

1
0,387(54,243 𝑥106 )6
𝑁𝑢𝐿𝐸 = 0,825 + 8
9 27
0,492 16
(1 + ( ) )
( (0,709) )

𝑁𝑢𝐿𝐸 = 50,926

Con la ecuación 1.11 se despeja el coeficiente de convección natural al exterior de la

cámara y resulta la ecuación:

𝑁𝑢𝐿𝐸 ∗ 𝑘
ℎ𝑒𝑥 =
𝐿𝑐

Con k que es dato tomado del Anexo 3.

(50,926 ) ∗ (0,0274 𝑊⁄𝑚 ∗ °𝐶 )

ℎ𝑒𝑥 =
(0,40𝑚)

ℎ𝑒𝑥 = 3,488 [𝑊⁄𝑚2 ∗ °𝐶 ]

3.2.2 Cálculo de la pérdida de calor por convección al exterior del secador

(radiación despreciable)
𝑄̇𝑝𝑒𝑟𝑑𝑖𝑑𝑎𝑠−𝑐𝑜𝑛𝑣𝑒𝑐𝑐𝑖ó𝑛 = ℎ̅ ∗ 𝐴 ∗ 𝛥𝑇

𝑄̇𝑝𝑒𝑟𝑑𝑖𝑑𝑎𝑠−𝑐𝑜𝑛𝑣𝑒𝑐𝑐𝑖ó𝑛 = 3,488 (𝑊⁄𝑚2 ∗ 𝐾 ) ∗ (0,54𝑚2 ) ∗ (303 − 293)𝐾

𝑄̇𝑝𝑒𝑟𝑑𝑖𝑑𝑎𝑠−𝑐𝑜𝑛𝑣𝑒𝑐𝑐𝑖ó𝑛 = 18,835 𝑊

3.3 CÁLCULO DE LA RESISTENCIAS TÉRMICAS POR CONDUCCIÓN Y

CONVECCION EN LAS PAREDES DEL PROTOTIPO DE LA CÁMARA DE
SECADO
Es necesario calcular las áreas de las caras de los componentes que conforman el
secador semiindustrial de la siguiente manera:

53
Para el cálculo de las respectivas áreas de la cámara de secado con el espesor del
acero galvanizado se toma las medidas descritas en la figura 3.1.

Figura 3.1: Dimensiones de la cámara de secado. [Autor]

A= 2[(0,9*0,6) +(0,9*0,4) +(0,6*0,4)]

A= 2,28 m2

A1= 2[(0,9028*0,6028) +(0,9028*0,4028) +(0,6028*0,4028)]

A1= 2,301 m2

Para el cálculo de las áreas restantes que son necesarias obtener para los cálculos
siguientes tomamos las dimensiones de la carcasa con el espesor del acaro
galvanizado ver figura 3.2.

54
 Figura 3.2: Dimensiones de la carcasa del secador. [Autor]

A2= 2[(1,0028*0,7028) +(1,0028*0,5028) +(0,7028*0,5028)]

A2= 3,124 m2

A3= 2[(1,0056*0,7056) +(1,0056*0,5056) +(0,7028*0,5056)]

A3= 3,149 m2

El espesor del aislante térmico a utilizar en este cálculo es de 50 milímetros de

espesor como se detalla en el capítulo II y el material es de lana de vidrio.

El secador semiindustrial está conformado por tres capas de materiales utilizados

para conformar las paredes del sistema en su totalidad, además de los cálculos
térmicos se debe calcular las resistencias térmicas que se describen a continuación:

- Dos capas de acero galvanizado de 1,4 mm de espesor (𝑘1 , 𝑘3 )

- Una capa de lana de vidrio de 50mm de espesor (𝑘2 )

La resistencia total se calcula con las siguientes ecuaciones:

𝑅𝑇𝑜𝑡𝑎𝑙−𝑐𝑜𝑛𝑑𝑢𝑐𝑐𝑖𝑜𝑛 = 𝑅1 + 𝑅2 + 𝑅3

𝐿1 𝐿2 𝐿3
𝑅𝑇𝑜𝑡𝑎𝑙−𝑐𝑜𝑛𝑑𝑢𝑐𝑐𝑖𝑜𝑛 = + +
𝑘1 𝐴1 𝑘2 𝐴2 𝑘3 𝐴3

55
Donde:

𝑅𝑇𝑜𝑡𝑎𝑙 = resistencia total

L = espesor de la pared

k = conductividad térmica del material

A = área de la superficie de la pared

Del Anexo 3 y Anexo 5 se obtiene los datos de conductividad térmica de los

materiales acero galvanizado y lana de vidrio respectivamente descritos en la tabla
3.5.

Tabla 3.5: Datos de conducción de materiales. [Autor]

Material Conductividad térmica (W/m.°C)

Acero galvanizado 46,5

Lana de vidrio 0,046

Datos:

𝐿1 = 𝐿3 = 1,4 mm = 1,4x10-3 m

𝐿2 = 50 mm = 0,05 m

𝑘1 = 𝑘3 = 46,5 W/m.°C

𝑘2 = 0,046 W/m.°C

𝐴1 = 2,301 m2

𝐴2 = 3,124 m2

𝐴3 = 3,149 m2

56
3.3.1 Resistencia térmica total por conducción
Reemplazamos los datos en la ecuación de la resistencia térmica total por
conducción.
(1,4𝑥10−3 𝑚) (0,05𝑚)
𝑅𝑇𝑜𝑡𝑎𝑙−𝑐𝑜𝑛𝑑 = +
𝑊 𝑊
(46,5 𝑚. °𝐶 ) (2,301𝑚2 ) (0,046 𝑚. °𝐶 ) (3,124𝑚2 )

(1,4𝑥10−3 𝑚)
+
𝑊
(46,5 𝑚. °𝐶 ) (3,149𝑚2 )

°𝐶
𝑅𝑇𝑜𝑡𝑎𝑙−𝑐𝑜𝑛𝑑 = 0,3479
𝑊

3.3.2 Resistencia térmica por convección natural al exterior del secador

(radiación despreciable)
Para el cálculo de esta resistencia no se considera la radiación y además se toma
como dato al coeficiente medio de convección natural ya calculado anteriormente y
se toma como referencia a la ecuación 1.3 para convección de la siguiente manera:

1
𝑅5 = [°𝐶⁄𝑊 ]
ℎ̅. 𝐴3

Donde:

𝑅5 = resistencia térmica por convección al exterior

ℎ̅ = coeficiente medio de convección natural al exterior

𝐴3 = área superficial exterior

Datos:

ℎ̅ = 6,775 (𝑊⁄𝑚2 ∗ °𝐶 )

𝐴3 = 3.149 m2

1
𝑅5 =
6,775 (𝑊⁄𝑚2 ∗ °𝐶 ) (3,149𝑚2 )

𝑅5 = 0,0468 °𝐶⁄𝑊

57
3.3.3 Cálculo del calor perdido total por las paredes del prototipo de secador
Una vez calculadas las resistencias térmicas es necesario calcular el calor disipado
por las paredes por lo que se emplea la siguiente formula:

𝑇𝑠 − 𝑇∞
𝑄̇𝑝𝑎𝑟𝑒𝑑 = Ec.3.7
∑𝑅

𝑇𝑠 − 𝑇∞
𝑄̇𝑝𝑎𝑟𝑒𝑑 =
𝑅𝑇𝑜𝑡𝑎𝑙−𝑐𝑜𝑛𝑑 + 𝑅5

Reemplazando las resistencias encontradas y reemplazando en la ecuación 3.8

tenemos:

(30 − 20)°𝐶
𝑄̇𝑝𝑎𝑟𝑒𝑑 =
(0,3479 + 0,0468) °𝐶⁄𝑊

𝑄̇𝑝𝑎𝑟𝑒𝑑 = 25 𝑊

3.4 CÁLCULO DE LA VARIACION DE LA ENERGIA INTERNA

ALMACENADA EN EL AIRE EN LA CÁMARA DEL SECADOR
3.4.1 Cálculo del calor que se necesita para elevar la temperatura del aire en el
interior de la cámara de secado (qair).
La temperatura máxima que alcanza el secador al interior de la cámara es un dato
experimental alcanzando 65°C y la temperatura del aire al exterior del secador es de
20°C. El calor necesario para elevar la temperatura del aire en el interior de la
cámara se expresa con la siguiente ecuación:

𝑞𝑎𝑖𝑟 = 𝑚. 𝐶𝑝. 𝛥𝑇

Donde:

𝑞𝑎𝑖𝑟 = cantidad de calor total de transferencia de calor

m = masa del sistema

58
Cp = calor específico del aire

𝛥T = diferencia de temperaturas

Con los datos obtenidos de la tabla del Anexo 2, tenemos a una temperatura de 45
grados centígrados (45+273= 318 K) los siguientes parámetros de la propiedad del
aire

Datos:

Volumen de la cámara de sacado (V) = 0,214 m3

Densidad del aire (ρ) = 1,112816 kg/m3

Calor específico del aire (Cp) = 1,00708 kJ/kg. °C

De esta forma reemplazamos los datos en la ecuación de calor total por transferencia,
sabiendo además que la masa es igual al producto de la densidad por el volumen y se
tiene que:

𝑞𝑎𝑖𝑟 = 𝑚. 𝐶𝑝. 𝛥𝑇

𝑘𝑔 𝑘𝐽
𝑞𝑎𝑖𝑟 = (0,214𝑚3 ) (1,112816 3
) (1,00708 ) (65 − 20)°𝐶
𝑚 𝑘𝑔. °𝐶

𝑞𝑎𝑖𝑟 = 10,772 𝑘𝐽

𝑞𝑎𝑖𝑟 = 10772,1 𝐽

De la prueba experimental se obtiene que, para alcanzar la temperatura de 65°C

(usada para el secado de cubetas de cartón para huevos), el tiempo de secado
requerido es de 11,47 minutos (688,2 s) ver Tabla 3.6.

3.4.2 Cálculo de la variación de la energía interna almacenada

𝑞𝑎𝑖𝑟
𝛥𝐸̇𝑎𝑙𝑚 =
𝑡

10772,1 𝐽
𝛥𝐸̇𝑎𝑙𝑚 =
688,2 𝑠

59
 𝛥𝐸̇𝑎𝑙𝑚 = 15,65 𝑊

3.5 BALANCE DE ENERGÍA

3.5.1 Cálculo del calor total de entrada requerido
El calor total requerido se calcula por la siguiente ecuación:

𝑄̇𝑇 = 𝑄̇𝑐𝑛 + 𝑄̇𝑓𝑜𝑟𝑧 + 𝑄̇𝑟𝑎𝑑 + 𝑄̇𝑝𝑒𝑟𝑑𝑖𝑑𝑎𝑠−𝑐𝑜𝑛𝑣𝑒𝑐𝑐𝑖𝑜𝑛 + 𝛥𝐸̇𝑎𝑙𝑚

Ec.3.8
+ 𝑜𝑡𝑟𝑎𝑠 𝑝𝑒𝑟𝑑𝑖𝑑𝑎𝑠 𝑛𝑜 𝑐𝑢𝑎𝑛𝑡𝑖𝑓𝑖𝑐𝑎𝑑𝑎𝑠

Donde:

𝑄̇𝑇 = flujo de calor total requerido (𝑊)

𝑄̇𝑐𝑛 = flujo de calor generado por el aire (𝑊)

𝑄̇𝑓𝑜𝑟𝑧 = flujo de calor por convección forzada (𝑊)

𝑄̇𝑟𝑎𝑑 = flujo de calor generado por radiación (𝑊)

𝑄̇𝑝𝑒𝑟𝑑𝑖𝑑𝑎𝑠−𝑐𝑜𝑛𝑣𝑒𝑐𝑐𝑖𝑜𝑛 = pérdidas por convección al exterior (W)

ΔĖalm = variación de energía almacenada

1000 W = (164,632 W) + ( 722,439 W) + (48,70 W) + (18,835 W)

+ (15,65 W) + otras pérdidas

otras pérdidas = 29,744 W

Las pérdidas no cuantificadas pueden producirse debido a pérdidas de la abertura en

el orificio del ventilador, defectos de aislamiento y apertura de puertas.

Pérdida por Pérdidas no

𝑄̇𝑐𝑛 y 𝑄̇𝑓𝑜𝑟𝑧 cuantificadas
convección
exterior
𝑄̇𝑟𝑎𝑑
Q̇perdidas−conveccion

𝛥𝐸̇𝑎𝑙𝑚

60
 Figura 3.3: Diagrama del balance de energía. [Autor]

3.6 CÁLCULO DE LA EFICIENCIA DEL PROCESO DE SECADO (η)

Q̇rad + ΔĖalm
𝛈= ∗ 100% Ec.3.9
Ėentra

48,7W + 15,65W
𝛈= ∗ 100%
1000 W

𝛈 = 𝟔, 𝟒%

Esta es la eficiencia de la cámara en vacío (sin carga), considerando solamente el

calentamiento del aire al interior. Estos valores se incrementarán en función de la
carga a ser secada, pudiendo ser para este tipo de cámara: cartón, probetas de aserrín
y materiales orgánicos.

3.7 ETAPAS DE SIMULACIÓN TÉRMICA

En esta sección de simulación térmica del sistema, permite realizar comparaciones de
los valores obtenidos por medio de las ecuaciones matemáticas aplicadas al proceso.
Además, permite obtener mediciones con el fin de conocer el comportamiento de la
distribución de temperatura al interior de la cámara de secado y su comportamiento
general.

Para esta simulación se tomó en cuenta las siguientes consideraciones de entrada

para que la validación del modelo sea lo más satisfactoria posible.

- Se simula todo el secador industrial con sus respectivos componentes y

especificaciones
- Se registran los datos de entrada como los parámetros y propiedades de los
materiales.
- El mallado se lo deja por defecto definido por el programa.
- Se requiere lo siguiente: al llegar a la temperatura deseada en la cámara de
secado con el tiempo para el equilibrio térmico, entonces los cálculos de los
parámetros serian correctos.

61
3.7.1 Validación del modelo
3.7.1.1 Volumen de la cámara de secado
En este diseño de cámara de secado es necesario obtener el volumen a ser analizado,
ya que el software a utilizar pide una geometría válida para el análisis mostrado en la
figura 3.4.

Figura 3.4: Volumen de la cámara. [Autor]

3.7.1.2 Cálculo del volumen de la cámara

V= A*B*C

Datos:

A; largo = 0,6 m

B; ancho = 0,9 m

C; altura = 0,4 m

Entonces:

V = A*B*C

V = (0,6)*(0,9)*(0,4)

V = 0,21 m3

62
3.8 ANÁLISIS DE LA DISTRIBUCIÓN (GLOBAL) DE TEMPERATURA
REAL AL INTERIOR DE LA CÁMARA DE SECADO
3.8.1 Temperaturas en el centro de la cámara de secado
Este análisis es experimental con datos tomados a un tiempo de 11,4 minutos, para el
cual se alcanza la temperatura de 65°C, a la cual se encuentra calibrado el
controlador de la cámara, como referencia se toma la mitad de la cámara de secado y
se extrae las siguientes temperaturas con sus tiempos respectivos por medio de una
termocupla tipo k, ver tabla 3.6.

Termocupla
Mitad de la
cámara

Figura 3.5: Ubicación central de la termocupla. [Autor]

Tabla 3.6 Temperaturas de media cámara. [Autor]

Tiempo (min) Tiempo (s) Temperatura (°C)

0 0 20

1 60 21,2

2 120 26,7

3 180 32,6

4 240 38,2

5 300 43,1

6 360 47,6

7 420 51,2

63
 8 480 54,4

9 540 57

10 600 59,4

11 660 61,7

11,47 688,2 65

T vs t
70
60
Temperatura (°C)

50
40
30
20
10
0
0 100 200 300 400 500 600 700 800
tiempo (s)

Gráfica 3.1: Distribución de temperaturas experimentales con respecto al tiempo en

la mitad de la cámara de secado. [Autor]

3.8.2 Temperaturas en la parte inferior de la cámara de secado

Al igual que el apartado anterior se obtuvo los datos de temperaturas, colocando la
termocupla tipo k en la parte inferior de la cámara de secado, ver tabla 3.7.

64
 Termocupla al
inferior de la
cámara

Figura 3.6: Ubicación inferior de la termocupla. [Autor]

Tabla 3.7 Temperaturas en la parte inferior de la cámara. [Autor]

Tiempo (min) Tiempo (s) Temperatura (°C)

0 0 20

1 60 22

2 120 23

3 180 27

4 240 33

5 300 38

6 360 42

7 420 46

8 480 50

9 540 53

10 600 59

65
 11 660 61

11,53 691,8 65

T vs t
70

60

50
Temperatura (°C)

40

30

20

10

0
0 100 200 300 400 500 600 700 800
tiempo (s)

Gráfica 3.2: Distribución de temperaturas experimentales con respecto al tiempo en

la parte inferior de la cámara de secado. [Autor]

3.9 SIMULACIÓN DE DISTRIBUCIÓN DE TEMPERATURAS AL

INTERIOR DE LA CÁMARA EN ESTADO TRANSITORIO

Para una mejor apreciación de las distribuciones de temperaturas en el interior de la

cámara de secado se analiza por tramos de tiempos de 180, 420 y 688 segundos,
además se muestran en diferentes vistas las simulaciones realizadas.

66
3.9.1 Convección y radiación en el interior de la cámara en estado transitorio
para un tiempo de 180s
Vista isométrica

Figura 3.7: Convección y radiación para un tiempo de 180s. [Autor]

Vista frontal

Figura 3.8: Distribución de temperaturas para un tiempo de 180s. [Autor]

67
Temperatura del flujo de aire a 0,4m/s de velocidad de entrada.

Figura 3.9: Temperaturas del aire para un tiempo de 180s. [Autor]

Temperaturas en diferentes puntos de la cámara (mitad de cámara)

Figura 3.10: Temperaturas en la mitad de la cámara para t=180s. [Autor]

68
 Tabla 3.8: Distribución de temperaturas en la mitad de la cámara. [Autor]

X [m] Y [m] Z [m] Temperatura (Fluido) [°C]

0,160145327 0,059978064 0,702800035 31,6606353
0,326953292 0,059978064 0,702800035 31,89011095
0,493761301 0,059978064 0,702800035 31,5302777
-0,006662704 0,059978064 0,702800035 31,59193796
-0,173470691 0,059978064 0,702800035 31,54323579
0,160145342 0,059978064 0,902800024 32,58372363
0,326953322 0,059978064 0,902800024 32,02031375
0,493761301 0,059978064 0,902800024 40,5293697
-0,006662656 0,059978064 0,902800024 31,95785872
-0,173470661 0,059978064 0,902800024 40,0284314
0,160145327 0,059978064 1,102800012 31,7929099
0,326953292 0,059978064 1,102800012 31,8474795
0,493761301 0,059978064 1,102800012 31,54221738
-0,006662704 0,059978064 1,102800012 31,83080049
-0,173470691 0,059978064 1,102800012 31,51421127

Temperatura en diferentes puntos de la cámara (cara inferior)

Figura 3.11: Temperaturas en la cara inferior de la cámara para t=180s. [Autor]

69
 Tabla 3.9: Distribución de temperaturas en la cara inferior de la cámara. [Autor]

X [m] Y [m] Z [m] Temperatura (Fluido) [°C]

0,520145357 -0,13922064 0,902800024 27,5326479
0,520145357 -0,13922064 1,102800012 26,83089783
0,520145357 -0,13922064 0,702800035 26,94443517
0,34014535 -0,13922064 0,902800024 28,51592199
0,34014535 -0,13922064 1,102800012 27,46874682
0,34014535 -0,13922064 0,702800035 27,6603037
0,160145357 -0,13922064 0,902800024 30,6527809
0,160145357 -0,13922064 1,102800012 28,04559001
0,160145357 -0,13922064 0,702800035 27,99125567
-0,019854635 -0,13922064 0,902800024 29,56017024
-0,019854635 -0,13922064 1,102800012 28,77599284
-0,019854635 -0,13922064 0,702800035 28,41704206
-0,199854627 -0,13922064 0,902800024 29,94164515
-0,199854627 -0,13922064 1,102800012 28,99468895
-0,199854627 -0,13922064 0,702800035 28,7778738

3.9.2 Convección y radiación en el interior de la cámara en estado transitorio

para un tiempo de 420s
Vista isométrica

Figura 3.12: Convección y radiación para un tiempo de 240s. [Autor]

70
Vista frontal

Figura 3.13: Distribución de temperaturas para un tiempo de 240s. [Autor]

Temperatura del flujo de aire a 0,4m/s de velocidad de entrada.

Figura 3.14: Temperaturas del aire para un tiempo de 420s. [Autor]

71
Temperaturas en diferentes puntos de la cámara (mitad de cámara)

Figura 3.15: Temperatura en la mitad de la cámara para t=240s. [Autor]

Tabla 3.10: Temperatura media de la cámara a 420s. [Autor]

X [m] Y [m] Z [m] Temperatura (Fluido) [°C]

0,160145311 0,052159418 0,9028 44,79678985

72
 Temperatura en diferentes puntos de la cámara (cara inferior)

Figura 3.16: Temperaturas en la cara inferior de la cámara para t=420s. [Autor]

Tabla 3.11: Temperaturas en la cara inferior de la cámara a 420s. [Autor]

X [m] Y [m] Z [m] Temperatura (Fluido) [°C]

0,520145357 -0,147842586 0,902800024 41,34014267
0,520145357 -0,147842586 1,102800012 37,46693871
0,520145357 -0,147842586 0,702800035 37,43048545
0,34014535 -0,147842586 0,902800024 44,03624619
0,34014535 -0,147842586 1,102800012 39,87089089
0,34014535 -0,147842586 0,702800035 40,47959917
0,160145357 -0,147842586 0,902800024 46,03707342
0,160145357 -0,147842586 1,102800012 41,69317006
0,160145357 -0,147842586 0,702800035 41,71422655
-0,019854635 -0,147842586 0,902800024 50,75186494
-0,019854635 -0,147842586 1,102800012 44,97149319
-0,019854635 -0,147842586 0,702800035 44,08501097
-0,199854627 -0,147842586 0,902800024 55,7057369
-0,199854627 -0,147842586 1,102800012 45,37138441
-0,199854627 -0,147842586 0,702800035 45,35066386

73
3.9.3 Convección y radiación en el interior de la cámara en estado transitorio
para un tiempo de 688s
En esta prueba se utilizó un tiempo estimado de 11,46 minutos (688s), ya que en la
prueba experimental en el laboratorio se llegó a 65 grados centígrados y en la
simulación de la cámara en el software obtenemos las diferentes imágenes de
distribuciones de temperaturas descritas a continuación en las diferentes vistas para
una mejor apreciación.

Vista isométrica

Figura 3.17: Convección y radiación para un tiempo de 688s. [Autor]

Vista frontal

Figura 3.18: Distribución de temperaturas para un tiempo de 688s. [Autor]

74
Temperatura del flujo de aire a 0,4m/s de velocidad de entrada.

Figura 3.19: Temperaturas del aire para un tiempo de 688s. [Autor]

Temperaturas en diferentes puntos de la cámara (mitad de cámara)

Figura 3.20: Temperaturas en la mitad de la cámara para t=688s. [Autor]

Tabla 3.12: Temperaturas en la mitad de la cámara a 688s. [Autor]

X [m] Y [m] Z [m] Temperatura (Fluido) [°C]

0,160145327 2,19369E-05 0,752799988 47,27098469
0,460145295 2,19369E-05 0,752799988 46,78775746

75
 -0,139854714 2,19369E-05 0,752799988 46,95083269
0,160145327 2,19369E-05 1,05279994 47,27627047
0,460145295 2,19369E-05 1,05279994 46,7378967
-0,139854714 2,19369E-05 1,05279994 47,01415058

Temperatura en diferentes puntos de la cámara (cara inferior)

Figura 3.21: Temperaturas en la cara inferior de la cámara para t=688s. [Autor]

Tabla 3.13: Temperaturas en la cara inferior de la cámara a 688s. [Autor]

X [m] Y [m] Z [m] Temperatura (Fluido) [°C]

0,520145357 -0,13922064 0,902800024 48,95335398
0,520145357 -0,13922064 1,102800012 51,44040163
0,520145357 -0,13922064 0,702800035 51,74113343
0,34014535 -0,13922064 0,902800024 58,41055896
0,34014535 -0,13922064 1,102800012 53,83750471
0,34014535 -0,13922064 0,702800035 53,96134729
0,160145357 -0,13922064 0,902800024 58,75555008
0,160145357 -0,13922064 1,102800012 55,26061453
0,160145357 -0,13922064 0,702800035 55,31706205
-0,019854635 -0,13922064 0,902800024 64,59551183
-0,019854635 -0,13922064 1,102800012 58,57805581
-0,019854635 -0,13922064 0,702800035 58,14373573
-0,199854627 -0,13922064 0,902800024 52,72717071

76
 -0,199854627 -0,13922064 1,102800012 58,9960443
-0,199854627 -0,13922064 0,702800035 58,80808125

3.10 VALIDACIÓN DEL COMPORTAMIENTO DEL SISTEMA

Se realiza el análisis para comparar entre la solución obtenida computacionalmente y
la obtenida experimentalmente o sistema real, para la obtención de la información se
utiliza el error absoluto para cuantificar variaciones con respecto al valor real o
verdadero que es el valor de las temperaturas obtenidas de la simulación realizada,
para lo cual se utiliza las siguientes ecuaciones:

E𝐚𝐛𝐬 = Medidaexperimental − Medidasimulada Ec.3.10

E𝐚𝐛𝐬
E𝐫 = ∗ 100% Ec.3.11
Medidasimulada

Tabla 3.14: Análisis experimental y simulado de las zonas estudiadas. [Autor]

Análisis Tiempo (s) Temperatura Temperatura

Central (°C) Cara inferior (°C)

Experimental 180 32,6 27

420 51,2 46

688 65 65

Simulado 180 32,58 30,65

420 44,79 46,04

688 47 58,75

Para los intervalos de tiempos mostrados en las tablas 3.15 y 3.16 se procede a
calcular el error absoluto y relativo de las zonas estudiadas.

77
 Tabla 3.15: Error absoluto y relativo en la zona central. [Autor]

Tiempo (s) Error absoluto (°C) Error relativo (%)

180 0,02 0,06

420 6,41 14

688 18 38

Tabla 3.16: Error absoluto y relativo en la cara inferior. [Autor]

Tiempo (s) Error absoluto (°C) Error relativo (%)

180 3,65 12

420 0,04 0,08

688 6,25 10,6

3.11 SIMULACIÓN DE DISTRIBUCIÓN DE TEMPERATURAS AL

INTERIOR DE LA CÁMARA EN ESTADO ESTABLE

3.11.1 Distribución de temperaturas por convección y radiación en estado

estable.
En esta prueba de simulación como es en estado estable no depende del tiempo y el
software calcula por medio de iteraciones solamente y se aprecian en las figuras las
distribuciones de temperaturas en diferentes vistas mostradas a continuación:

78
Temperatura en diferentes puntos de la cámara (mitad de cámara)

Figura 3.22: Temperaturas en la mitad de la cámara para estado estable. [Autor]

Tabla 3.17: Distribución de temperaturas en la mitad de la cámara en estado estable.

[Autor]

X [m] Y [m] Z [m] Temperatura (Fluido) [°C]

0,160145327 0,050779354 0,752799988 138,1675554
0,460145295 0,050779354 0,752799988 135,2464754
-0,139854714 0,050779354 0,752799988 135,6131442
0,160145327 0,050779354 1,05279994 133,134392
0,460145295 0,050779354 1,05279994 137,4264795
-0,139854714 0,050779354 1,05279994 134,3410406
Temperatura en diferentes puntos de la cámara (cara inferior)

Figura 3.23: Temperaturas en la cara inferior de la cámara para estado estable.

[Autor]

79
 Tabla 3.16: Distribución de temperaturas en la cara inferior de cámara en estado
estable. [Autor]

X [m] Y [m] Z [m] Temperatura (Fluid) [°C]

0,510145307 -0,13922064 0,902800024 125,7323429
0,510145307 -0,13922064 1,102800012 119,6798925
0,510145307 -0,13922064 0,702800035 118,5104388
0,310145319 -0,13922064 0,902800024 120,0697236
0,310145319 -0,13922064 1,102800012 120,2257472
0,310145319 -0,13922064 0,702800035 119,8093775
0,110145316 -0,13922064 0,902800024 119,442535
0,110145316 -0,13922064 1,102800012 118,0804695
0,110145316 -0,13922064 0,702800035 121,2637549
-0,089854687 -0,13922064 0,902800024 119,8726118
-0,089854687 -0,13922064 1,102800012 119,0148017
-0,089854687 -0,13922064 0,702800035 119,7426074

En estado estable las temperaturas fluctúan entre 120°C y 135°C al interior de la

cámara esto podría permitir su uso para otras aplicaciones como para el secado de
pintura con tiempos más elevados de proceso.

80
 CAPÍTULO IV

CONCLUSIONES Y RECOMENDACIONES
4.1 CONCLUSIONES
• Aplicando balance de energía se determinaron perdidas por convección hacia
el exterior que alcanzan aproximadamente 25 W, se analizó la transferencia
de calor interna por convección natural, forzada y radiación y sus
combinaciones que alcanzan los 771 W, la variación de energía almacenada
para el aire es de 15,65 W

• Al combinar los resultados de los métodos de transferencia de calor por

convección forzada con radiación y convección natural con radiación, la
eficiencia de la cámara en vacío (sin carga) y considerando solamente el
calentamiento del aire al interior nos da 6,4%.

• La simulación realizada en el software ayudó a comprobar los valores

calculados mediante ecuaciones y los valores transmitidos por el programa, se
analizó en la zona media e inferior de la cámara la distribución de
temperaturas en función de diferentes intervalos de tiempo de 180, 420 y 688
segundos, alcanzando valores máximos de 47°C en el centro y 58,75°C en la
zona inferior.

• Con la simulación en el software se proporciona una mayor proximidad de la

variación de distribuciones de temperaturas en la cámara, mostrando que los
puntos analizados computacionalmente y los realizados experimentalmente
varían de acuerdo al tiempo y a las condiciones iniciales del programa.

• De acuerdo a los datos obtenidos experimentalmente el tiempo mínimo para

alcanzar la temperatura de 65 grados centígrados al interior de la cámara de
secado es de 11,4 minutos (688s), mientras que en la obtención de datos de la
simulación realizada en el estudio las distribuciones de temperaturas en los
puntos analizados (medio e inferior) son menores. La diferencia entre los

81
 valores experimentales y simulados se debe principalmente a la dificultad de
medición usando termocupla convencional.

4.2 RECOMENDACIONES

• Es importante que la simulación incluya diversos parámetros para la

convección natural y forzada que permitan que la solución se aproxime a la
real.

• Se recomienda utilizar las termocuplas necesarias instaladas adecuadamente

al sistema, para la obtención de datos sin caídas de temperaturas al interior
del modelo a ser analizado.

• Para un mejor análisis del prototipo se recomendaría utilizar lámparas

infrarrojas tipo tubo para una mejor distribución de temperaturas y así
alcanzar una mejor eficiencia del secador.

• Sería importante complementar el estudio con simulaciones que incluyan

carga para el secado que permitirían verificar la calidad del proceso.

82
BIBLIOGRAFÍA
[1] J. C. Espinoza Freire, “‘ESTUDIO DEL PROCESO DE SECADO CON
LÁMPARAS INFRARROJAS Y SU INFLUENCIA SOBRE LA
DETERMINACIÓN DEL GRADO DE HUMEDAD EN CUBETAS DE
CARTÓN DE HUEVOS PRODUCIDAS EN LA FÁBRICA AVÍCOLA
UBICADA EN SAMANGA,’” Universidad Técnica de Ambato, 2013.

[2] G. F. Ortega Sánchez, “‘DISEÑO, CONSTRUCCIÓN E

IMPLEMENTACIÓN DE UN PROTOTIPO DE UN HORNO DE SECADO
(CURADO) DE PINTURA AUTOMOTRIZ Y PRUEBAS DE PINTURA EN
LAS PROBETAS AL FINAL DEL PROCESO,’” UNIVERSIDAD
TECNOLÓGICA EQUINOCCIAL, 2016.

[3] F. Bravo Bravo, “‘SIMULACION Y ESCALAMIENTO INDUSTRIAL DE

UN SECADOR INFRARROJO PARA EL ACABADO DEL CUERO,’”
CIATEC, 2015.

[4] J. P. Holman, Tranferencia de calor, Octava edi. Madrid: Mc Graw Hill, 1999.

[5] K. Donald, Procesos de tranferencia de calor, Trigésima. México:

CONTINENTAL, 1999.

[6] J. R. Welty, Fundamentos de tranferencia de momento, calor y masa, Primera

ed. Corvallis, Oregón: LIMUSA, 1960.

[7] Dassault Systemes, “Conducción,” 2018. https://n9.cl/fifu%0A (accessed Sep.

18, 2020).

[8] Y. A. Cengel, Transferencia de calor y masa, Cuarta Edi. México: Mc Graw

Hill, 2011.

[9] F. P. Incropera, FUNDAMENTOS DE TRANFERENCIA DE CALOR,

Cuarta Edi. México: Pearson, 1999.

[10] R. L. MOTT, Mecánica de Fluidos Aplicada, Cuarta Edi. 2004.

[11] Dassault Systemes, “Radiation,” 2018.

http://help.solidworks.com/2018/spanish/SolidWorks/cworks/c_Radiation.htm
(accessed Sep. 18, 2020).

83
[12] A. J. Chapman, TRANSMISIÓN DE CALOR, Tercera Ed. Madrid: GREFOL,
1984.

[13] G. Pérez, “ESPECTROMETRÍA.COM,” 2007.

https://www.espectrometria.com/espectrometra_infrarroja (accessed Sep. 18,
2020).

[14] M. del C. Cañote V., “ESPECTROMETRÍA DE ABSORCIÓN

INFRARROJA,” vol. 0, 2011, [Online]. Available:
https://www.scribd.com/presentation/69212797/INFRARROJO.

[15] Xamán J., Dinámica De Fluidos Computacional Para Ingenieros, Primera ed.
EE.UU: Palibrio, 2015.

[16] A. Morillo Velázques, “EVAPORADORES Y SECADORES,” Equipos,

máquinas e Instal. Ind., [Online]. Available: www.docsity.com.

[17] Camen, “SISTEMAS DE SECADO Y SECADORES INDUTRIALES,” 2008.

camen-q.com/blog/sistemas-de-secado-y-secadores-industriales.html
(accessed Oct. 08, 2020).

[18] M. Peralta, “Secadores Industriales,” Academia.edu, 1998, [Online].

Available: https://www.academia.edu/9626426/Secadores_Industriales.

[19] P. Fito Maupoey, A. M. Andrés Grau, J. M. Barat Baviera, and A. M. Albors

Sorolla, Introducción al secado de alimentos por aire caliente. España: Byprint
Percom, 2001.

[20] LAGOS S.A, “Equipos Infrarrojos Lagos,” 1999.

https://cabinaslagos.com/secadores-infrarrojos-ir-taller-industria-pintura/
(accessed Oct. 11, 2020).

[21] POWDERTRONIC, “HORNOS INFRARROJOS.”

https://powdertronic.com/hornos-infrarrojos/ (accessed Oct. 12, 2020).

84
ANEXOS

Anexo 1. Temperatura de calentamiento del prototipo por convección forzada y

radiación respecto al tiempo

Fuente: [1]

85
Anexo 2: Tabla de propiedades térmicas del aire

86
Anexo 3: Valores de propiedades para metales (acero galvanizado)

87
Anexo 4: Propiedades del gas ideal del aire

88
Anexo 5: Propiedades térmicas de la lana de vidrio

89
Anexo 6: Controladores de funcionamiento del prototipo de secador

90
Anexo 7: Entrada de aire al prototipo por ventilación a 0,4m/s (convección
forzada)

91
Anexo 8: Instalación de las lámparas infrarrojas de 250 W

92
Anexo 9: Ubicación de las lámparas infrarrojas en el interior de la cámara.

93
Anexo 10: Termocuplas instaladas en la zona media e inferior de la cámara de
secado

94