C. Mira, "Diseño, Construcción y Prueba de Un Prototipo Calentador Solar de Agua Con Tubos de Calor.,"
C. Mira, "Diseño, Construcción y Prueba de Un Prototipo Calentador Solar de Agua Con Tubos de Calor.,"
C. Mira, "Diseño, Construcción y Prueba de Un Prototipo Calentador Solar de Agua Con Tubos de Calor.,"
UNIVERSIDAD EAFIT
DEPARTAMENTO DE INGENIERÍA MECÁNICA
ÁREA DE DISEÑO
MEDELLÍN
2007
DISEÑO, CONSTRUCCIÓN Y PRUEBA DE UN PROTOTIPO CALENTADOR
SOLAR DE AGUA CON TUBOS DE CALOR
Asesor:
Luis Santiago París
Ingeniero Mecánico
UNIVERSIDAD EAFIT
DEPARTAMENTO DE INGENIERÍA MECÁNICA
ÁREA DE DISEÑO
MEDELLÍN
2007
AGRADECIMIENTOS
Durante el desarrollo del proyecto recibí ayuda de muchas personas a los cuales
agradezco profundamente por el conocimiento que me compartieron, el tiempo
que me brindaron para aportar con su trabajo y el constante apoyo que sentí de
ellos.
Gracias, a mi familia por su apoyo en esta y otras ocasiones. Ellos soportaron gran
parte de la carga emocional del proyecto.
pág
TABLA DE CONTENIDO 6
LISTA DE ILUSTRACIONES 10
LISTA DE TABLAS 13
LISTA DE ANEXOS 15
1 INTRODUCCIÓN 13
2 OBJETIVOS 16
2.1 GENERAL 16
2.2 ESPECÍFICOS 16
3 MARCO TEÓRICO 17
3.1 LA ENERGÍA SOLAR 17
3.1.1 Radiación Extraterrestre 17
3.1.2 Radiación solar en la superficie terrestre 19
3.1.3 Geometría de la radiación solar 24
3.1.4 Medición de la radiación solar 27
3.2 CALENTADORES SOLARES 28
3.2.1 Componentes un sistema de calentamiento solar 28
3.2.2 Colector de Placa Plana 30
3.2.3 Colector de tubos de vacío 31
3.2.4 Colectores Concentradores 32
3.2.5 Superficies para absorción 33
3.2.6 Características de medios transparentes 35
3.2.7 Eficiencia óptica de un colector solar de placa plana 37
3.3 TUBOS DE CALOR 39
3.3.1 Ventajas y aplicaciones 40
3.3.2 Componentes de los tubos de calor 41
3.3.3 Límites de transferencia de calor 43
3.3.4 Variaciones de los tubos de calor 47
3.3.5 Modelo de cálculo de la resistencia térmica equivalente 50
3.4 CALENTADORES SOLARES CON TUBOS DE CALOR 52
4 DISEÑO DE UN PROTOTIPO DE CALENTADOR SOLAR CON TUBOS DE
CALOR 56
4.1 PLANTEAMIENTO DEL PROBLEMA DE DISEÑO 56
4.1.1 Condiciones de operación 56
4.1.2 Especificación de diseño el prototipo de calentador solar 59
4.2 ANÁLISIS FUNCIONAL 60
4.2.1 Caja Negra del prototipo de calentador solar 61
4.2.2 Estructura funcional del prototipo de calentador solar. 62
4.2.3 Síntesis Funcional del prototipo de calentador solar. 62
4.3 MATRIZ MORFOLÓGICA 65
4.4 EVALUACIÓN PREVIA 66
4.5 CONCEPTOS SOLUCIÓN 68
4.6 EVALUACIÓN DE ALTERNATIVAS 68
4.6.1 Árbol de objetivos 68
4.6.2 Valores de cumplimiento 69
4.6.3 Matriz de evaluación 70
4.7 CONFIGURACIONES DE FORMA 71
4.7.1 Intercambiador 71
4.7.2 Placa 72
4.7.3 Cubierta 74
4.8 SELECCIÓN DE MATERIALES 75
4.8.1 Placa Colectora 76
4.8.2 Cubierta transparente 78
4.8.3 Soporte 80
4.8.4 Intercambiador 82
4.8.5 Aislamientos 83
4.9 DIMENSIONAMIENTO DE LA SOLUCIÓN 85
4.9.1 Cálculo aislamiento tanque 85
4.9.2 Estimación de las pérdidas de calor en el panel 87
4.9.3 Capacidad de transferencia de calor de la placa 97
4.9.4 Cálculo convección natural en el intercambiador 102
5 DISEÑO DE TUBOS DE CALOR PARA EL PROTOTIPO DE
CALENTADOR SOLAR 105
5.1 PLANTEAMIENTO DEL PROBLEMA DE DISEÑO 105
5.2 ANÁLISIS FUNCIONAL 107
5.2.1 Caja Negra del tubo de calor 107
5.2.2 Síntesis Funcional del tubo de calor 107
5.3 EVALUACIÓN DE ALTERNATIVAS DE SOLUCIÓN 110
5.3.1 Fluido 110
5.3.2 Contenedor 111
5.3.3 Estructura Capilar 113
5.4 DIMENSIONAMIENTO DE LOS TUBOS DE CALOR 114
5.4.1 Geometría 114
5.4.2 Resistencia Mecánica 115
5.4.3 Cantidad de Fluido 116
5.4.4 Límites de Operación 118
6 ASPECTOS IMPORTANTES EN LA CONSTRUCCIÓN DE LOS TUBOS
DE CALOR PARA EL PROTOTIPO DE CALENTADOR SOLAR 120
6.1 LIMPIEZA 120
6.2 ENSAMBLE 126
6.3 CIERRE 127
6.4 SOLDADURA 128
6.5 CARGA 129
6.6 SELLADO 131
7 EVALUACIÓN DEL PROTOTIPO DE CALENTADOR SOLAR CON TUBOS
DE CALOR 132
7.1 INSTRUMENTOS DE MEDICIÓN Y REGISTRO 132
7.1.1 Medición de la radiación solar 132
7.1.2 Medición de temperaturas y calibración de termocuplas 133
7.2 PRUEBA EQUILIBRIO DE TEMPERATURA 134
7.3 PRUEBA CIRCULACIÓN NATURAL 136
7.4 PRUEBA CIRCULACIÓN FORZADA 140
7.5 PRUEBA VARIACIÓN DE INCLINACIÓN EN UN TUBO 143
7.6 COSTOS 147
8 CONCLUSIONES 149
9 RECOMEDACIONES Y SUGERENCIAS 153
10 BIBLIOGRAFÍA 155
10.1 CLÁSICA 155
10.2 INTERNET 157
LISTA DE ILUSTRACIONES
pág
Ilustración 1 Variación de la radiación solar extraterrestre durante el año 18
Ilustración 2 Espectro de emisión del sol 19
Ilustración 3 Balance Radiación Solar 21
Ilustración 4 Masa de aire 22
Ilustración 5 Esquema Latitud 25
Ilustración 6 Declinación Solar 25
Ilustración 7 Ángulos posición superficie 26
Ilustración 8 Ecuación del tiempo 27
Ilustración 9 Esquema de un piranómetro 28
Ilustración 10 Componentes de un calentador solar de agua 30
Ilustración 11 Colector de Placa Plana 31
Ilustración 12 Colector de Tubos de Vacío 32
Ilustración 13 Principio de concentración 33
Ilustración 14 Distribución energía incidente 33
Ilustración 15 Ley de Kirchhof 35
Ilustración 16 Coeficiente de transmisión de un medio transparente 36
Ilustración 17 Radiación absorbida por la placa 38
Ilustración 18 Esquema del funcionamiento de un tubo de calor 39
Ilustración 19 Estructuras Capilares para tubos de calor 42
Ilustración 20 Esquema de operación de un termosifón 48
Ilustración 21 Esquema de operación de un tubo de calor de circuito cerrado 49
Ilustración 22 Resistencia térmica equivalente de un tubo de calor 50
Ilustración 23 Resistencia térmica equivalente simplificada 51
Ilustración 24 Calentador solar de agua compacto con tubos de calor 54
Ilustración 25 Colector Solar con Tubos de Calor Evacuados 55
Ilustración 26 Mapa de radiación solar global en Colombia 58
Ilustración 27 Caja negra de un calentador solar 62
Ilustración 28 Estructura funcional de un prototipo de calentador solar 62
Ilustración 29 Síntesis funcional de un prototipo de calentador solar 64
Ilustración 30 Árbol de objetivos 69
Ilustración 31 Configuraciones de forma para el intercambiador 72
Ilustración 32 Configuraciones de forma para la placa absorbedora 73
Ilustración 33 Configuraciones de forma para la cubierta 74
Ilustración 34 Pérdida de estabilidad dimensional de las cubiertas individuales 75
Ilustración 35 Esquema para el cálculo del aislamiento del tanque 86
Ilustración 36 Esquema general de los flujos de calor 89
Ilustración 37 Pérdidas por conducción en la parte inferior 89
Ilustración 38 Pérdidas por convección en la parte superior 91
Ilustración 39 Pérdidas por radiación en la parte superior 93
Ilustración 40 Gráfico eficiencia del colector 97
Ilustración 41 Esquema placa para análisis térmico 98
Ilustración 42 Flujos de calor para tramo de la placa analizada 98
Ilustración 43 Distribución de temperatura 102
Ilustración 44 Esquema transferencia de calor por convección 103
Ilustración 45 Caja negra del tubo de calor 107
Ilustración 46 Síntesis funcional del tubo de calor 109
Ilustración 47 Esquema cálculo volumen malla 116
Ilustración 48 Estructura Malla 118
Ilustración 49 Enjuague de los tubos 121
Ilustración 50 Secado de los tubos 121
Ilustración 51 Dispositivo para el lavado de los tubos 122
Ilustración 52 Esquema dispositivo para el lavado de los tubos 122
Ilustración 53 Inmersión de la malla en solución 124
Ilustración 54 Enjuague de la malla con abundante agua 124
Ilustración 55 Proceso secado malla 125
Ilustración 56 Verificación de limpieza malla 125
Ilustración 57 Inmersión de capilares y termopozos en solución 126
Ilustración 58 Ensamble malla y tubo 126
Ilustración 59 Esquema copa para cierre de tubos de calor 127
Ilustración 60 Proceso cierre tubo 128
Ilustración 61 Soldadura capilares 128
Ilustración 62 Soldadura termopozos 129
Ilustración 63 Pesaje del tubo 129
Ilustración 64 Vacío en el tubo 130
Ilustración 65 Llenado del tubo 130
Ilustración 66 Circuito para la medición de la radiación solar. 133
Ilustración 67 Esquema prueba equilibrio de temperatura 135
Ilustración 68 Resultado prueba equilibrio de temperatura 135
Ilustración 69 Esquema prueba de circulación natural 136
Ilustración 70 Evolución temperaturas en prueba con circulación natural 137
Ilustración 71 Eficiencia en la prueba con circulación natural 138
Ilustración 72 Distribución temperaturas del tubo en prueba circulación natural 139
Ilustración 73 Esquema prueba circulación forzada 140
Ilustración 74 Evolución temperaturas en prueba con circulación forzada 141
Ilustración 75 Eficiencia en prueba con circulación forzada 142
Ilustración 76 Distribución temperatura tubo en prueba circulación forzada 143
Ilustración 77 Esquema prueba tubo con diferentes inclinaciones 144
Ilustración 78 Montaje prueba para un tubo 144
Ilustración 79 Distribución de temperatura con una inclinación de 12° 145
Ilustración 80 Distribución de temperatura con una inclinación de 18° 145
Ilustración 81 Distribución de temperaturas con una inclinación de 25° 146
Ilustración 82 Degradación de la malla 147
LISTA DE TABLAS
pág
Tabla 1 Coeficientes para el modelo de Hottel 23
Tabla 2 Radiación solar en Medellín 57
Tabla 3 Especificación de diseño del prototipo de calentador solar 59
Tabla 4 Subfunciones de la síntesis funcional 63
Tabla 5 Matriz Morfológica 65
Tabla 6 Criterios para la evaluación previa de principios de solución 66
Tabla 7 Evaluación previa de los principios de solución 67
Tabla 8 Concepto de Solución 1 68
Tabla 9 Concepto de Solución 2 68
Tabla 10 Criterios de cumplimiento con los objetivos 69
Tabla 11 Matriz de evaluación de los conceptos de solución 70
Tabla 12 Valores para la calificación de aspectos cualitativos 76
Tabla 13 Matriz de evaluación para la selección del material en la placa 78
Tabla 14 Matriz de evaluación para la selección del material de la cubierta. 80
Tabla 15 Matriz de evaluación para la selección del material del soporte 81
Tabla 16 Matriz de evaluación para la selección del material del intercambiador 83
Tabla 17 Matriz de evaluación para la selección del material del aislamiento 84
Tabla 18 Dimensiones tanque de almacenamiento 85
Tabla 19 Parámetros para el cálculo del aislamiento del tanque 86
Tabla 20 Valores de entrada para cálculo de las pérdidas térmicas 88
Tabla 21 Parámetros aislamiento parte inferior panel 90
Tabla 22 Pérdidas de calor en la parte inferior 90
Tabla 23 Parámetros para cálculo de coeficiente por convección natural 92
Tabla 24 Resistencia a la convección entre la placa y la cubierta 92
Tabla 25 Resistencia a la convección entre la cubierta y el ambiente 93
Tabla 26 Resistencia térmica a la radiación entre la placa y la cubierta 94
Tabla 27 Resistencia térmica a la radiación entre la placa y el ambiente 96
Tabla 28 Resultados eficiencia térmica 96
Tabla 29 Flujos de calor en la placa 99
Tabla 30 Parámetros de entrada para la distribución de temperatura de a placa 101
Tabla 31 Parámetros y resultados cálculo convección natural 104
Tabla 32 Especificación de diseño del tubo de calor 106
Tabla 33 Relación entre funciones parciales y componentes del tubo de calor 108
Tabla 34 Matriz Selección fluido de trabajo 111
Tabla 35 Matriz selección material contenedor tubo de calor. 112
Tabla 36 Matriz selección mecha 114
Tabla 37 Dimensiones tubo de calor 115
Tabla 38 Parámetros de entrada y resultados cálculo resistencia del tubo 115
Tabla 39 Parámetros y resultados cálculo volumen malla 117
Tabla 40 Propiedades estructura capilar 119
Tabla 41 Límites de operación de los tubos de calor 119
Tabla 42 Partes del dispositivo para el lavado de los tubos 123
Tabla 43 Eficiencia prueba convección natural 138
Tabla 44 Eficiencia en prueba con circulación forzada 142
Tabla 45 Costos materiales del prototipo 148
LISTA DE ANEXOS
pág
ANEXO PLANOS 161
1 INTRODUCCIÓN
13
El calentamiento solar de agua puede llevarse a cabo con distintos tipos de
equipos, los cuales tienen un funcionamiento basado en el mismo principio básico.
Un cuerpo oscuro absorbe energía de la luz solar y transfiere este calor absorbido
a un líquido.
Los tubos de calor han sido un tema de investigación en países como China,
Estados Unidos, Rusia desde hace muchos años. Allí se han desarrollado
importantes aplicaciones industriales y de consumo particular; además, se
producen y comercializan habitualmente calentadores solares con tubos de calor.
Las aplicaciones que sacan provecho de la energía solar son muy comunes en
países desarrollados con condiciones climáticas adversas. Resulta curioso que en
las zonas tropicales, donde el nivel de radiación es relativamente alto y constante
a lo largo del año, este tipo de tecnologías tengan aplicaciones tan reducidas.
Existen en Colombia otras fuentes de energía tradicional cuya utilización es más
simple. Sin embargo, puede resultar más costosa su explotación, y aún hoy hay
áreas remotas donde no hay cubrimiento de los servicios públicos.
14
Colombia, específicamente en la ciudad de Medellín, y se pretende evaluar el
desempeño del mismo.
15
2 OBJETIVOS
2.1 GENERAL
2.2 ESPECÍFICOS
16
3 MARCO TEÓRICO
El flujo de energía proveniente del sol sobre el borde exterior de la atmósfera varía
de manera inversa proporcional al cuadrado de la distancia entre el sol y la tierra,
la cual tiene variaciones a lo largo del año. Al valor promedio de dicho flujo por
unidad de área se le llama constante solar, su valor medio es de 1353W/m2 el cuál
fluctúa levemente por la variación periódica de la distancia entre la tierra y el sol.
La radiación extraterrestre se define como el flujo de energía que incide sobre una
superficie perpendicular a la dirección de propagación de la radiación fuera de la
atmósfera. Para un día particular el flujo puede encontrarse como (ITESO@,
2006):
17
360n
I on = I SC 1 + 0.033 cos
365 (1)
1400
Radiación Extraterrestre W/m2
1380
1360
1340
1320
1300
1280
1260
ro
zo
e
ril
lio
re
e
er
ay
ni
br
br
br
st
Ab
ne
ub
Ju
ar
Ju
go
br
m
M
M
E
ct
ie
Fe
ie
ie
A
ov
pt
ic
D
Se
Mes
1
Cuerpo Negro: Cuerpo cuya superficie absorbe toda la radiación electromagnética incidente
18
Ilustración 2 Espectro de emisión del sol
WIEDER, 2007
19
menor que la constante solar debido a la absorción y a la dispersión de la
radiación que origina la interacción de los fotones con la atmósfera. Entonces,
aparecen clasificaciones para la radiación en la superficie terrestre.
• Radiación directa (Ib,Hb): Se recibe directamente del sol sin sufrir ninguna
dispersión atmosférica.
• Radiación difusa (Id,Hd): se recibe del sol después de ser afectada por la
dispersión atmosférica.
• Radiación Terrestre: Proviene de objetos terrestres por efectos reflectivos.
• Radiación Total: Es la suma de las radiaciones directa, difusa y terrestre.
• Radiación Global: Es la suma de las radiaciones directa y difusa. No se
considera la radiación terrestre (ITESO@,2006).
20
Ilustración 3 Balance Radiación Solar
Existen modelos para el cálculo de la radiación directa para días soleados con
cielo claro, los cuales tienen en cuenta la absorción de radiación de la atmósfera.
Uno de los más comunes es el de la masa de aire o ley de Bourger, el cual
establece que (KREITH Y OTROS, 1978):
I b = I o e − km (2)
21
masa que atravesaría si el sol estuviera directamente perpendicular al plano del
piso. Cuando el sol esta directamente vertical m=1, cuando la altitud solar es igual
a 30º, m=2. La masa de aire se define como:
m = sec θ z (3)
θz
α
τ atm =
Ib
(
= 0.5 e − 0.65 m ( z ,α ) + e −0.095 m ( z ,α ) ) (4)
Io
P( z )
m( z , α ) = m(0, α ) (5)
P(0)
22
[
m(0, α ) = 1229 + (614 senα ) 2 ] 1
2
− 614 senα (6)
Sin embargo estas ecuaciones solo aplican con cielo claro en una atmósfera
estándar ausente de contaminación. Hottel hizo modificaciones a estas
ecuaciones para tener en cuenta las partículas de vapor (KREITH Y KREIDER,
1978):
23
Existen otros modelos para la determinación de la radiación global diaria, como la
ecuación de regresión de tipo Angstrom.
H S
= a+b (9)
Ho So
Donde H es la radiación total promedio para el día sobre una superficie horizontal,
Ho es la radiación solar extraterrestre, S son las horas reales de brillo solar y So las
horas posibles de iluminación solar para el día. Las constantes a y b dependen del
lugar, el clima y la vegetación y se encuentran en los siguientes rangos:
• 0.1 ≤ a ≤ 0.54
• 0.1 ≤ b ≤ 0.73
I T = I n cos θ (10)
24
Ilustración 5 Esquema Latitud
Φ
Línea del Ecuador
• Declinación Solar (δ): Posición angular del sol al medio día, respecto al
plano del Ecuador, es positivo hacia el norte. La causa de este ángulo es la
inclinación que tiene la tierra respecto al eje terrestre de rotación de 23.45º.
Esta inclinación determina la existencia de las estaciones y las distintas
duraciones de los días. (-23.5º ≤ δ ≤ 23.5º).
Mar 21 23.5º
Jun 21 Dic 21
SOL
Sep 21
25
Ilustración 7 Ángulos posición superficie
26
Ilustración 8 Ecuación del tiempo
20
15
10
0
0 50 100 150 200 250 300 350 400
-5
-10
-15
-20
Día n
27
de estar expuestos a la radiación, hay una diferencia de temperatura entre ellos.
La hipótesis de trabajo de un piranómetro, es que la irradiancia es directamente
proporcional a la diferencia de temperatura entre ambos sensores. Para evitar
ruido en las lecturas, causado por el viento y otros factores meteorológicos, el
sensor expuesto a la radiación suele estar protegidos por un hemisferio de vidrio.
Este hemisferio, de características ópticas adecuadas, permite el paso de la
radiación, pero evita el enfriamiento por viento, lo cual puede alterar la lectura.
ITESO@,2006
28
disposición espacial. Sin embargo se pueden identificar ciertos elementos
comunes a la mayoría de los sistemas de calentamiento.
• Colector solar térmico: Capta la energía calorífica del sol. En otras palabras
transforma la energía solar incidente en otra forma de energía útil. Existen
diferentes tipos de colectores solares pero los más comunes son los de
placa plana y los de tubos de vacío.
• Acumulador: Almacena el agua que se pretende calentar. Consiste de un
tanque aislado térmicamente. En cuyo interior puede encontrarse el
intercambiador de calor.
• Intercambiador de calor: Permite transmitir el calor que absorbe el líquido
caloportador2 en el colector al agua contenida en el acumulador.
• Sistema de bombeo: Garantiza la circulación del fluido caloportador entre el
colector y el intercambiador de calor.
• Sistema de Tuberías: conduce los flujos de los líquidos y permite las
uniones entre las partes del sistema.
• Sistema de Respaldo: Consiste en un elemento adicional, como resistencia
eléctrica, que entra en funcionamiento cuando la captación de radiación
solar es insuficiente para suplir la demanda.
• Elementos de control y seguridad: Hacen posible mantener ciertas variables
del sistema en niveles establecidos.
2
Caloportador: Transporta calor
29
Ilustración 10 Componentes de un calentador solar de agua
GARBITEK@,2006
El tipo de colector más popular es el de placa plana. Consiste en una caja aislada
y hermética para soportar las adversidades del clima. En la parte superior tiene
una cubierta transparente de plástico o vidrio que favorece el efecto invernadero al
permitir que pasen los rayos del sol y evitar pérdidas de calor por convección
(GARBITEK@,2006). Una placa absorbente se ubica debajo de la cubierta. Esta
placa tiene un recubrimiento oscuro que incrementa el poder de absorción.
30
Ilustración 11 Colector de Placa Plana
EVE@, 2006
31
Ilustración 12 Colector de Tubos de Vacío
GARBITEK@,2006
32
Ilustración 13 Principio de concentración
Energía
Absorbida
α
Energía
Transmitida
τ
33
• Coeficiente de reflexión o reflectancia (ρ): Fracción de energía que se
reflecta desde el material.
α +τ + ρ = 1 (14)
α + ρ =1 (15)
E (T )
ε = (16)
Eb (T )
Eb (T ) = σT 4
(17)
σ = 5.6697 * 10 W / m K
8 2 4
34
Ilustración 15 Ley de Kirchhof
εEA
αIA
EA = GA (19)
35
La reflectancia de un medio transparente se expresa como una función del índice
de refracción, mediante la siguiente expresión (LUNDE, 1980):
n −1
2
ρ = (21)
n + 1
1 ρ (1-ρ)2ρ (1-ρ)2ρ3
(1-ρ) (1-ρ)ρ2
(1-ρ)ρ
Τr
∞
τ r = (1 − ρ ) + (1 − ρ ) ρ 2 + (1 − ρ ) ρ 4 + K = (1 − ρ ) ∑ρ
2 2 2 2 2i
(22)
i =0
1− ρ
τr = (23)
1+ ρ
Para un número n de cubiertas se tiene que:
36
1− ρ
τr = (24)
1 + (2n − 1) ρ
1 1 − ρ1 1 − ρ2
τr = +
2 1 + (2n − 1) ρ1 1 + (2n − 1) ρ 2
(25)
37
Ilustración 17 Radiación absorbida por la placa
Radiación
incidente
τ
τ(1-α)
τα
τα(1-α)ρd τα(1-α)2ρd2
τα
∞ i
38
3.3 TUBOS DE CALOR
3
Termopozo: Funda hecha de material conductivo que sirve para separar el sensor de temperatura
del medio de medición.
39
3.3.1 Ventajas y aplicaciones
Los tubos de calor suelen ser utilizados en situaciones en las que: la fuente de
calor y el sumidero deben estar separadas, es necesario incrementar la
conductividad de un sólido, o se quiere esparcir el calor en una superficie. Su
capacidad de transferencia de calor puede ser muchas veces mayor que la del
conductor sólido. Esto se traduce en una resistencia térmica pequeña y permite la
separación del evaporador y el condensador sin un aumento significativo en la
caída de temperatura a lo largo del tubo (PETERSON, 1994).
El ciclo bifásico en el que operan los tubos de calor favorece cortos tiempos de
respuesta. Estos son menores que otros aparatos para la transferencia de calor y
no son función de la longitud.
Los tubos de calor son partes inmóviles y son casi libres de mantenimiento. Sin
embargo, debe prestarse atención en su diseño y manufactura para lograr un buen
desempeño. Hay dos aspectos de la manufactura de especial cuidado: la limpieza
y el sellado. Un defecto en la primera puede generar gases no condensables que
significan la pérdida de funcionalidad de un tramo del tubo; y un defecto de la
segunda se traduce en un fallo definitivo del tubo. También existen condiciones de
uso que pueden reducir la vida útil de los tubos de calor como: vibración, impactos
fuertes, choque térmico y ambientes corrosivos.
40
3.3.2 Componentes de los tubos de calor
41
Ilustración 19 Estructuras Capilares para tubos de calor
42
y buena estabilidad térmica química. Esta última es muy importante, ya que
se refiere a la posible degradación del fluido al someterlo a ciertas
temperaturas. Los calores de vaporización altos permiten el transporte de
grandes cantidades de calor con bajos flujos másicos y así mantienen una
caída de presión pequeña dentro del tubo. La buena conductividad térmica
minimiza el gradiente de temperatura radial y reduce la posibilidad de
núcleos de ebullición. La baja viscosidad minimiza la resistencia al flujo y la
alta tensión superficial evita que partículas de líquido sean arrastradas por
el vapor.
2 ⋅ σ ⋅ cos θ
∆Pc = (27)
rc
43
presión capilar. Si esto no es así, se destruye la interfase líquido – vapor,
produciéndose secado y sobrecalentamiento del evaporador (APPLEBY Y
OTROS; 2003). El flujo de calor límite puede aumentarse con una reducción
de los poros de la malla; sin embargo debe tenerse cuidado ya que esto
reduce la permeabilidad.
El límite capilar se encuentra como:
2 ⋅ σ ⋅ cosθ
± ρ l ⋅ g (d v ⋅ cosψ + L ⋅ senψ )
qc =
rc
(28)
C ⋅ ( f v Re v ) ⋅ µ v µL
+ Leff
2 ⋅ (rb ,v ) ⋅ Av ⋅ ρ v ⋅ h fg K ⋅ Aw ⋅ h fg ⋅ ρ l
2
44
2 ⋅ rb ,v ⋅ q c
( f v Re v ) = 0.038
C = 1.00
A ⋅ µ
v v fg ⋅ h
2 ⋅ rb,v ⋅ q c
Re v = (31)
Av ⋅ µ v ⋅ h fg
qc
Ma v = (32)
Av ρ v h fg Rv ⋅ Tv ⋅ γ v
• Límite de arrastre: Los flujos de vapor y líquido dentro del tubo son en
direcciones opuestas, la interacción de estos genera una fuerza de arrastre
sobre el líquido debido a una mayor velocidad en el vapor. La tensión
superficial se encarga de balancear esta fuerza, pero existe un flujo de calor
límite con el que se alcanzan velocidades de vapor capaces de arrastrar
gotas de líquido. La cantidad de líquido que llega al evaporador disminuye y
este se seca y sobrecalienta (SILVERSTEIN, 1992). El modelo de cálculo
para el límite de arrastre es (PETERSON, 1994):
45
1
σ ⋅ ρ v2 2
qa = Av ⋅ h fg ⋅
(33)
2 ⋅ r
b,w
2 ⋅ π ⋅ Le ⋅ keff ⋅ Tv 2 ⋅ σ
qe = r − ∆Pc ,m (35)
h fg ⋅ ρ v ⋅ ln ri n
r
v
46
vapor y la velocidad del mismo. La velocidad máxima del vapor es la del
sonido, al alcanzarla el flujo de calor no puede aumentar más y llega a un
límite. El flujo de calor máximo que puede transferirse a cierta temperatura,
se encuentra utilizando la ecuación (PETERSON, 1994).
qs = 0.474 ⋅ h fg ⋅ AV ⋅ ρV P (36)
d i2 ⋅ h fg ⋅ AV ⋅ PV ⋅ ρV
qv = (37)
64 ⋅ µV ⋅ Le
47
del condensador al evaporador deslizándose por la pared del contenedor,
en la pared se forma una película de líquido. Por esta razón, los
termosifones bifásicos están restringidos a tener el evaporador en la parte
inferior y el condensador en la parte superior.
Condensador
Sumidero de
Calor
Flujo de Vapor
Zona Adiabática
Película de
condensado
descendiendo
Evaporador
Fuente de
Calor
Pozo de líquido de
evaporación
48
Ilustración 21 Esquema de operación de un tubo de calor de circuito cerrado
Flujo de Vapor
Condensador
Evaporador
Fuente de Sumidero
Calor de Calor
Flujo de Líquido
En los LHP la reserva esta conectada al evaporador, por una unión capilar,
permitiendo una regulación automática de la conductancia del tubo acorde
con el flujo de calor. Los LHP no necesitan precalentamiento antes la
operación (APPLEBY y OTROS).
49
pequeña para ser obligada a retornar al evaporador (GIRALDO y OTROS,
1988).
50
Ilustración 23 Resistencia térmica equivalente simplificada
d o, p
ln
di, p
Rp =
2πLK p
(39)
Donde do,p y di,p son el diámetro externo e interno del tubo respectivamente, L es
la longitud del evaporador o el condensador y Kp es la conductividad de la pared.
51
La resistencia térmica equivalente de la combinación líquido mecha se calcula de
manera similar utilizando una conductividad térmica efectiva del líquido y la
mecha.
d
ln o ,m
d
Rp =
i,m
(40)
2πLK eff
Así la resistencia térmica equivalente del tubo de calor queda definida como:
d o, p d
ln ln o, m
1 1 1 di, p
+ d i ,m
Req = + (41)
2π Le Lc K p K eff
Los tubos de calor tienen características muy importantes que hacen atractiva y
práctica su utilización en sistemas de calentamiento solar de agua. Algunas de
estas características ventajosas son: su capacidad de transferir calor latente a
altas tasas y grandes distancias y la corta diferencia de temperatura necesaria
entre la zona fría y la zona caliente lo que implica una baja degradación de la
calidad de la energía. (MATHIOULAKIS Y OTROS, 2000).
52
medio de un fluido caloportador y un sistema de bombeo, al tanque de
almacenamiento.
En la segunda categoría están los sistemas en los cuales el condensador del tubo
de calor esta dentro del acumulador. En este caso el calor se transfiere, por medio
del tubo de calor, desde el colector al fluido sin necesidad de otros fluidos o
sistemas de bombeo. Estos sistemas generalmente emplean termosifones
bifásicos y el tanque se ubica en la parte superior del colector. (MATHIOULAKIS Y
OTROS,2000).
53
Ilustración 24 Calentador solar de agua compacto con tubos de calor
4
Loop: lazo cerrado
54
Ilustración 25 Colector Solar con Tubos de Calor Evacuados
THERMOMAX@,2006
55
4 DISEÑO DE UN PROTOTIPO DE CALENTADOR SOLAR CON TUBOS DE
CALOR
El uso de la radiación solar como energía térmica parece ser la aplicación más
directa de esta energía amigable con el ambiente, disponible y gratuita. La
radiación es uno de los tres mecanismos básicos de transferencia de calor. Utilizar
este calor radiante desde el sol para aumentar la energía interna de una masa de
agua no requiere ningún tipo de transformación en la forma de la energía. Se
necesita diseñar un mecanismo que permita transferir el calor y prevenir de las
pérdidas del mismo.
56
La capacidad del colector esta determinada por la intensidad de radiación que este
reciba, y esta a su vez depende, entre otros factores, del lugar geográfico en el
que se ubique. El prototipo que se diseña, se ubica en Medellín, ciudad que tiene
las siguientes coordenadas geográficas:
• Latitud: 6º 13’
• Longitud: -75º 36’
• Altitud: 1490m
H 4,78 4,94 4,76 4,76 4,51 4,96 5,47 5,33 4,98 4,39 4,34 4,37 4,80
H/Ho 0,51 0,50 0,46 0,46 0,44 0,50 0,55 0,52 0,48 0,44 0,46 0,48 0,47
Hb 2,84 2,88 2,62 2,58 2,38 2,83 3,30 3,13 2,82 2,35 2,41 2,50 2,72
Hb/H 0,59 0,58 0,55 0,54 0,53 0,57 0,60 0,59 0,57 0,54 0,56 0,57 0,57
Hd 1,94 2,06 2,15 2,18 2,13 2,14 2,17 2,20 2,17 2,03 1,93 1,87 2,08
Hd/H 0,41 0,42 0,45 0,46 0,47 0,43 0,40 0,41 0,44 0,46 0,44 0,43 0,43
So 11,70 11,80 12,00 12,10 12,30 12,40 12,30 12,20 12,00 11,90 11,70 11,60 12,00
S 5,74 5,54 4,77 4,77 4,55 5,83 6,87 6,19 5,30 4,29 4,67 5,03 5,30
S/So 0,49 0,47 0,40 0,39 0,37 0,47 0,56 0,51 0,44 0,36 0,40 0,43 0,44
57
La determinación teórica de la radiación tiene cierta complejidad y un grado
considerable de incertidumbre. Por lo cual, se utilizan este tipo de modelos que
combinan mediciones experimentales y cálculos teóricos. En promedio la
irradiación solar en Medellín, para una superficie horizontal, es de 4,80kWh/m2/día
o 1800kWh/m2/año y las horas de brillo son 5.30 horas diarias o 2000 horas
anuales. La radiación instantánea promedio para una superficie horizontal es
cercana a 400W/m2.
58
En Colombia, en general los valores de radiación promedio son similares a los de
Medellín. Los valores extremos son 2.5kWh/m2 por día y 7.0kWh/m2 por día. La
radiación promedio para el territorio nacional se muestra en la ilustración.
Rm Número de usuarios
cuadrado de captación comerciales
Calentadores solares
Rm Tprom para uso doméstico 40ºC<= T <= 60ºC
comerciales
Valor promedio en
D Alta eficiencia térmica 50%
calentadores solares
Pérdida de calor durante la
Rf Aislamiento Efectivo ΔT <= 10ºC
noche
Cantidad de Agua disponible
Materia
59
Rm: Requisito Mínimo
Rf: Requisito Fijo
D: Deseo
Tipo Requerimiento Intervalos Criterio/Comentario
D Indicador de Temperatura --
Señal Facilidad para la instalación
Rf de instrumentos de medición -- Prototipo para pruebas
durante las pruebas
Rf Control de nivel de agua --
Manufactura Apariencia Seguridad Ambientales
Recursos y Tecnología
D Construcción Simple --
disponibles.
Costos Mantenimiento
Recursos y Tecnología
Rm Bajo Costo C <= $ 700.000
disponibles.
60
de partida para el diseño, ya que permite obtener un conjunto de funciones con
sus interrelaciones, a las cuales se les debe encontrar portadores en la etapa de
generación de alternativas.
61
Ilustración 27 Caja negra de un calentador solar
Flujo de Materia
Flujo de Energía
Flujo de Información
SUMINISTRAR
AGUA CALIENTE
Transferir Entregar
Calor al Agua
FUNCIONES CLAVES
Agua Caliente EN EL DISEÑO
62
Tabla 4 Subfunciones de la síntesis funcional
Esquema Subfunción Descripción Reto de diseño
Pérdidas Involucra un conjunto
de elementos que SI: Es una
Captar
Radiación permiten el flujo de función
Captar Radiación Calor radiación
solar calor eficiente a determinante en
Solar solar
través de una el desempeño.
superficie.
Pérdidas SI: Se debe
seleccionar el
Es la función
Transportar tipo,
desempeñada por
Calor Transportar Calor Calor Calor componentes y
los tubos de calor
forma del tubo
de calor.
Calor
El calor captado y
SI: Es una
Transferir transportado debe
Agua Agua función
Transferir Caliente calor al ser entregado al
Fría determinante en
Calor al agua agua agua, a través de
el desempeño.
una superficie
Agua Caliente
El agua que ingresa
al sistema y la que
Agua Agua NO: El diseño se
se calienta debe
Fría Almacenar
Fría Almacenar centra en el
almacenar para
Agua Pérdidas agua panel de
lograr su
colección.
disponibilidad para el
Agua Caliente consumo.
63
La síntesis funcional se muestra en la siguiente figura:
Pérdidas
Cargas
Transferir
Calor al agua
Nivel
Agua
Agua Caliente
Fría disponible
Agua Agua Agua
Agua Fría Ingresar Fría Sostener Fría Almacenar Caliente Entregar
Agua Agua Agua
Caliente
Agua Caliente
Pérdidas
Cargas Transmitidas
4.3 MATRIZ MORFOLÓGICA
Captar
Radiación
Solar
Placa plana para absorción Espejo con forma parabólica
con cubierta protectora. para aplicación del principio
de concentración.
Transportar
Calor
Tubo de calor recto con
Termosifón bifásico con estructura capilar para Tubo de calor de
condensador vertical facilitar la distribución del circuito cerrado
líquido en el evaporador
Transferir Calor Área del condensador
Aletas
al Agua recto del tubo de calor.
Almacenar
Tanque integrado al
mecanismo de captación
65
4.4 EVALUACIÓN PREVIA
66
Tabla 7 Evaluación previa de los principios de solución
Solución
Solución 1 Solución 2 Solución 3
Función
CAPTAR
Placa Plana Principio de concentración
RADIACIÓN SOLAR
Complejidad + -
-
Compatibilidad con +
Tiene un costo muy
especificaciones Bajo costo
elevado
-
Capacidad Técnica + No se posee la tecnología
para su producción
=
+
Aporte Conocimiento No es una solución
innovadora
TRANSPORTAR Termosifón Tubo de calor de
Tubo de calor
CALOR Bifásico circuito cerrado
Complejidad + = -
-
Compatibilidad con
+ + Incrementa los costos
especificaciones
del colector
+ = =
Se tiene Puede desarrollarse la Puede desarrollarse la
Capacidad Técnica
experiencia en su capacidad técnica para capacidad técnica para
fabricación producirlo producirlo
Aporte Conocimiento + + +
TRANSFERIR Área del
Aletas
CALOR AL AGUA condensador
Complejidad + =
Compatibilidad con -
+
especificaciones Incrementa costos
Capacidad Técnica + =
-
-
Solución
Aporte Conocimiento Solución ampliamente
ampliamente
conocida
conocida
Tanque separado
ALMACENAR Tanque integrado horizontal o vertical e
intercambiador de calor
Complejidad = =
=
Compatibilidad con
Puede ser más +
especificaciones
costoso.
+
Capacidad Técnica =
Aporte Conocimiento + =
Principio de solución pertinente
Principio de solución descartado
67
4.5 CONCEPTOS SOLUCIÓN
68
aceptable, debe satisfacer de una forma superior los deseos de la especificación.
Por esta razón los deseos corresponden en gran medida con los objetivos de
segundo nivel.
69
Los criterios cualitativos son traducidos numéricamente así: una clasificación alta
es igual a 1, una clasificación media es igual a 0.6, y una clasificación baja es igual
a 0.1.
70
4.7 CONFIGURACIONES DE FORMA
4.7.1 Intercambiador
El fluido a calentar pasa por el intercambiador donde entra en contacto con los
condensadores de los tubos de calor. Las propuestas para la geometría del
intercambiador son:
71
Ilustración 31 Configuraciones de forma para el intercambiador
Configuración de forma 1
Configuración de forma 2
Configuración de forma 3
4.7.2 Placa
72
• Doblado de la placa para que adquiera la forma del tubo: Los tubos pueden
ubicarse encima o debajo de la placa. El doblez de la semicircunferencia
incrementa el área de contacto con los tubos, lo que favorece la
transferencia de calor. Esta configuración implica un proceso adicional de
manufactura que puede incrementar los costos.
• Placa plana: No se realiza ningún doblez en la placa y el tubo se une
mediante abrazaderas de un material conductor. El área para la
transferencia de la placa es muy pobre y puede afectar el desempeño del
sistema.
Configuración de forma 1
Configuración de forma 2
73
4.7.3 Cubierta
Configuración de forma 1
Configuración de forma 2
74
Se estima que 9 tubos de calor tienen la capacidad de retirar la energía incidente
de 400W/m2. Además, así se facilita la posibilidad de separar un tubo para
pruebas y utilizar botellas plásticas como cubierta.
75
Los criterios de evaluación tienen valores de 0 a 1. Cuando son criterios
cualitativos, pueden tomar valores de alto, medio y bajo. Si la característica es
positiva una calificación de alto equivale a 1, una calificación de medio equivale a
0.6 y una calificación de bajo equivale a 0.1. Si la característica es negativa, la
calificación es inversa, alto equivale a 0.1, medio equivale a 0.6 y bajo equivale a
0.1.
76
La conductividad, la densidad y la capacidad calorífica se relacionan mediante la
difusividad térmica. Esta es una medida de que tan rápido un material puede
transferir el calor desde una fuente. Esta tiene en cuenta la conductividad, pero
también la capacidad calorífica y la densidad, ya que el material no solo transfiere
calor sino que también se calienta (LIENHARD y OTROS, 2004).
K
α= (1)
ρC p
Los materiales que normalmente se utilizan en los colectores solares son (LUNDE,
1980):
• Aluminio.
• Acero.
• Cobre.
• Latón.
Se construye una matriz de evaluación para la selección del material para la placa
colectora.
77
Tabla 13 Matriz de evaluación para la selección del material en la placa
Conductividad Térmica
Capacidad Calorífica
Difusividad Térmica Costo
Densidad (kg/m3)
(W/ºCm)
(J/kgºC)
Valor (m2/s)
Calificación
Calificación
Calificación
Calificación
Ponderada
Ponderada
Material Total
Costo
Valor
Peso
Peso
Aluminio 2722 921 206 8,22E-05 0,76 0,6 0,46 Medio 0,6 0,4 0,24 0,70
Acero 7860 448 45 1,28E-05 0,12 0,6 0,07 Medio 0,6 0,4 0,24 0,31
Cobre 8902 393 377 1,08E-04 1,00 0,6 0,60 Alto 0,1 0,4 0,04 0,64
Latón 8522 385 109 3,32E-05 0,31 0,6 0,18 Bajo 1 0,4 0,4 0,58
• Transmitancia (0.2).
• Resistencia a la intemperie (0.15).
• Resistencia a la temperatura (40ºC) (0.15).
• Fácil consecución (0.2)
• Costo (0.3).
78
El material tradicional en el que se fabrican estas cubiertas es el vidrio, sin
embargo pueden explorarse otras posibilidades en los polímeros. Los posibles
materiales para la cubierta son:
79
UV. Se estima su duración entre 2 ó 3 años para láminas flexibles, siendo
superior a 6 años para láminas rígidas (INFOAGRO@, 2007).
4.8.3 Soporte
80
Además el soporte esta expuesto a la intemperie. Las principales características
que se buscan en el material para el soporte son:
81
El material seleccionado es la madera plástica teniendo en cuenta el uso creciente
que tiene este material, su bajo peso y resistencia a la intemperie.
4.8.4 Intercambiador
Se selecciona el PVC, a pesar que tiene una baja resistencia a las temperaturas.
Su temperatura límite de operación es igual a la temperatura límite de operación
del sistema. Sin embargo, el sistema solo alcanza esta temperatura en casos
82
puntuales, y en situaciones normales el agua se calienta hasta 50ºC. Un material
más apropiado es el PE, pero este es de difícil consecución en el medio.
4.8.5 Aislamientos
Los aislamientos influyen notoriamente en la eficiencia del dispositivo. Su función
es disminuir las pérdidas térmicas, de esta manera se aumenta el calor útil. El
dispositivo esta expuesto a la lluvia y por eso la impermeabilidad del aislante se
83
debe considerar. Las propiedades que se buscan en el material de los
aislamientos son:
• Baja conductividad térmica (0.3).
• Impermeabilidad (0.3).
• Costos (0.4).
• Fibra de vidrio
• Espuma de poliuretano.
84
4.9 DIMENSIONAMIENTO DE LA SOLUCIÓN
85
Ilustración 35 Esquema para el cálculo del aislamiento del tanque
Tamb,h
a Ta T1 T2
h
L
L
Ta T1 T2 Tamb
t L 1
K t A K a. A hA
• Ta − Tamb Ta − Tamb
Q= =
Rt + R a + Rc t L 1
+ +
K t At K a Ac hAc
(4)
86
Si se supone este flujo de calor constante durante toda la noche se puede
encontrar la energía total cedida por el agua.
•
Q p = Q* Z
(5)
Donde Z es el tiempo sin radiación, que en este caso se supone 16h es decir
57600 segundos. Dando lugar a un calor de pérdida de 2468.12kJ. La disminución
de temperatura en el agua se calcula como:
Qp
∆T =
mC p
(6)
87
Qtrans U p ⋅ (T p − Ta )
η= = (τα ) − (9)
(τα ) Qi
Para efectos de calcular las resistencias térmicas se hace una suposición de los
valores de las temperaturas de la placa (Tp), la cubierta (Tc) y el ambiente (Ta),
teniendo en cuenta valores comunes en colectores solares. Además se supone en
flujo de radiación incidente de 700W/m2. Estos y otros valores de entrada para los
cálculos se especifican en la tabla.
En la parte inferior hay un aislamiento térmico, por lo tanto hay una resistencia a la
transferencia de calor por conducción. Las pérdidas por convección y radiación en
la parte inferior son pequeñas, por lo tanto, no se consideran.
88
Ilustración 36 Esquema general de los flujos de calor
Qp,1
Ta
R3,c R3,r
Tc
R1
Ta
Qp,2
L R1
Ta
Qp,2
1 ka
U p, 2 = = (10)
R1 L
89
Donde ka es la conductividad térmica del aislante y L es el espesor del mismo. Y
las el calor de pérdida en la parte inferior se expresa como:
Q p , 2 = U p , 2 (T p − Ta )
(11)
Así los valores del coeficiente global de pérdidas y el calor de perdida en la parte
inferior son:
90
Ilustración 38 Pérdidas por convección en la parte superior
Qc
Ta
R3,c
Tc
R2,c d
Tp
1
Rc =
h (12)
Nu = 0.152Gr 0.281
10 4 ≤ Gr ≤ 10 7 (13)
hd
Nu = (14)
k
91
gβ ∆Td 3
Gr =
ν2 (15)
1
β=
T (16)
92
Entre la cubierta y el medio ambiente se presenta un fenómeno de convección
forzada cuando hay viento pasando sobre la cubierta. En este caso el coeficiente
de convección se calcula como (PARIS Y OTROS, 1980):
Dónde V es la velocidad del viento. En este caso se supone una velocidad del
viento de 3m/s y se obtiene un coeficiente de convección entre la placa y el
ambiente de 17.10W/m2K y la resistencia térmica es de 0.058m2K/W. Estos
resultados se resumen en la tabla.
Ahora se propone establecer el modelo para el cálculo de los flujos de calor por
radiación.
Qr
Ta
R3,r
Tc
R2,r
Tp
El flujo de calor por radiación entre la placa y la cubierta se encuentra como el flujo
por radiación entre superficies paralelas (WIEDER, 1982).
93
(
J = ε p − cσ T p2 − Tc2 ) (18)
Donde,
1 1 1
= + −1
ε p−c ε p εc
J = h2.r (T p − Tc )
(19)
Donde,
3
h2, r = 4ε p −cσ T
T p + Tc
T=
2
94
La cubierta establece un flujo de calor por radiación con el cielo que se expresa
como (WIEDER,1982):
(
J = ε cσ Tc − Tcielo
4
) (20)
El flujo radiante entre la placa y el cielo también puede escribirse similar a un flujo
convectivo.
Donde,
3
hr = 4ε c T
Tc + Tcielo
T=
2
También puede encontrarse el flujo referido a la temperatura ambiente.
J = hr
(Tc − Tcielo ) (T − T )
(Tc − Ta ) c a (23)
Entonces es coeficiente de radiación entre la placa y el ambiente es el siguiente.
h3, r = 4ε c T
3 (Tc − Tcielo )
(Tc − Ta ) (24)
95
Tabla 27 Resistencia térmica a la radiación entre la placa y el ambiente
VARIABLE NOMBRE ECUACIÓN VALOR UNIDADES
Tcielo Temperatura Cielo 0,0552(Ta)^(1,5) 288,48 K
hr3 Coeficiente Radiación εc*σ(Tc^2+Tcielo^2)(Tc+Tcielo)(Tc- 11,95 W/m2K
Equivalente Tcielo)/(Tc-Ta)
R3,1 Resistencia Térmica 1/hr2 0,084 m2K/W
Radiación
1 1
Re = +
1 1 1 1
+ +
R 2, c R2 , r R3,c R3,r
(25)
Puede hallarse la eficiencia teórica del colector como una función del calor
incidente, la temperatura de la placa y la temperatura ambiente.
96
Qtrans U p (T p − Ta )
η= = (τα ) − (26)
Qi Qi
8.13(T p − Ta )
η = 0,83 −
Qi
90%
80%
70%
60%
Eficiencia
50%
40%
30%
20%
10%
0%
0 0,02 0,04 0,06 0,08 0,1 0,12
(Tp-Ta)/Qi (ºC m2 /W)
97
coeficiente de pérdida por unidad de área. Estas transferencias de calor ocurren
por la parte superior de la placa y la parte inferior se supone aislada. Por simetría,
el análisis puede hacerse únicamente en una porción de la placa, tomando un
tramo a la izquierda de un tubo, con longitud igual a la mitad de la separación
entre tubos.
qa Up, Ta
qp
qb
qs qe t x
dx
x=0 x=L
• •
E = S …….
98
qa + qe = q s + q p
(27)
dT
qs Calor que entra en x=x+Δx Kt
dx x + ∆x
dT dT
(τα )Qi ∆x + Kt = Kt + U p (T − Ta )∆x
dx x + ∆x
dx x (28)
d 2T U p
2
− (T − Ta ) + (Tα ) Qi = 0
dx Kt Kt (29)
Se define,
Up
m=
Kt (30)
99
Y se sustituye T por la temperatura en exceso θ=T-Ta. La ecuación queda escrita
como,
d 2θ (Tα )Qi 0
2
− m 2θ + m 2 =
dx Up
(Tα )Qi
θ = C1 e − mx + C 2 e mx + Up
(31)
dθ
=0
dx x =0
θ x=L
= Tb − Ta
(32)
(τα )Qi
T − Ta −
Up cosh( mx)
=
Tb − Ta −
(τα )Qi cosh(mL)
Up
(33)
100
Q=
1L
L ∫0
[ ]
(τα )Qi − U p (T − Ta ) dx
Tb − Ta − (τα )Qi
Up
Q=− tanh( mL)
mL Up
(34)
101
Ilustración 43 Distribución de temperatura
70,0000008
70,0000007
70,0000006
70,0000005
70,0000004
T(ºC)
70,0000003
70,0000002
70,0000001
70
69,9999999
0 0,01 0,02 0,03 0,04 0,05 0,06
x(m)
La propuesta es que el calentador solar pueda ser utilizado con o sin una bomba
que haga circular el fluido entre el tanque y el intercambiador. Cuando no se utiliza
una bomba, la transferencia de calor dentro del intercambiador ocurre por
convección natural y la circulación ocurre por el efecto termosifón.
102
calculan los coeficientes utilizando las relaciones empíricas encontradas en la
literatura.
Tagua
Tc
dext
Lc
h ⋅ d ext
Nu D = (34)
k
103
Los parámetros de entrada y los resultados del cálculo de la transferencia de calor
por convección natural en el condensador se presentan en la tabla.
El flujo de calor por tubo que se consigue con el mecanismo de convección natural
es de 51.93W. Este flujo de calor es apenas suficiente para evacuar el calor
absorbido, para lograr un mejor desempeño debe considerarse el uso de una
bomba o tener una mayor longitud de condensador.
104
5 DISEÑO DE TUBOS DE CALOR PARA EL PROTOTIPO DE CALENTADOR
SOLAR
105
Tabla 32 Especificación de diseño del tubo de calor
Rm: Requisito Mínimo
Rf: Requisito Fijo
D: Deseo
Tipo Requerimiento Intervalos Criterio/Comentario
Rf Geometría Longitud del evaporador Le = 1m Diseño del calentador solar
Rf Longitud del condensador Lc=0.1m Diseño del calentador solar
Latitudes extremas de
Rm Angulo de operación 0<= φ <= 20º
Colombia
Rf Forma del contenedor Recto Diseño del calentador solar
Rf Sección transversal Circular Diseño del calentador solar
Fuerza
Resistente a presiones
Rf Pmáx = 198 kPa Presión de saturación a 120ºC
internas
térmica
Arranque en frío cada
Rf -- Condición de operación
mañana.
Temperatura del panel sin
Rf Temperatura de arranque T = 20ºC
radiación.
Seguridad Manufactura Ambientales
Condición de operación. El
condensador está
Rf Resistente a la corrosión --
permanentemente inmerso en
el agua.
D Fácil instalación --
Recursos y Tecnología
Rf Bajo Costo --
disponibles.
106
5.2 ANÁLISIS FUNCIONAL
Energía de Energía de
entrada TRANSPORTAR salida
ENERGÍA TÉRMICA
Pérdidas
107
Tabla 33 Relación entre funciones parciales y componentes del tubo de calor
Componente Funciones Relacionadas Propiedades Componente
Contenedor • Recibir calor. • Resistencia química al medio de
• Conducir calor. la fuente y del sumidero de calor.
• Contener Vapor. • Conductividad térmica alta.
• Contener Líquido. • Resistencia a presiones de
• Entregar calor operación.
Malla • Transportar líquido. • Capilaridad.
• Conducir calor. • Permeabilidad.
• Contener Líquido. • Mojabilidad.
• Conductividad térmica alta.
Fluido • Evaporar líquido. • Estado líquido a temperaturas de
• Transportar vapor. operación.
• Contener vapor. • Alta entalpía de vaporización.
• Baja viscosidad del líquido.
• Alta tensión superficial.
108
Ilustración 46 Síntesis funcional del tubo de calor
Energía
Térmica Energía
Energía Energía Térmica
Entrada Térmica Térmica
Recibir Conducir Evaporar Transportar Vapor Contener Vapor Condensar Conducir Entregar Salida
calor Calor líquido vapor vapor Vapor calor calor
Líquido
Líquido
Flujo de energía
Flujo de masa
109
5.3 EVALUACIÓN DE ALTERNATIVAS DE SOLUCIÓN
Ahora se procede a evaluar entre cual es el material más apropiado para cada uno
de los componentes teniendo en cuenta la aplicación y las restricciones del
diseño.
5.3.1 Fluido
Las propiedades que el fluido de trabajo debe de tener asociadas a las funciones
parciales ya se han mencionado. Estas propiedades son función de la temperatura
y por ende la selección del fluido de calor depende de las temperaturas de
operación. Para temperaturas medias, como es el caso de la presente aplicación,
se sugieren materiales como el agua, el etanol y el metanol.
ρ L ⋅ σ ⋅ h fg
Nl =
µ (1)
Los fluidos con mayor número de mérito resultan más adecuados en las
aplicaciones de tubos de calor. La evaluación del fluido, se realiza teniendo en
cuenta el número de mérito y otros criterios de evaluación.
110
Tabla 34 Matriz Selección fluido de trabajo
CRITERIO DE CALIFICACIÓN Agua Etanol Metanol
Densidad kg/m3 998 810 792
Tensión Superficial N/m 0,07198867 0,022 0,022
Calor de Vaporización J/kg 2,44E+06 8,39E+05 1,13E+06
Número de Viscosidad Pa.s 8,90E-04 1,17E-03 6,10E-04
Mérito Valor (J/m2 s) 1,97E+11 1,28E+10 3,22E+10
Calificación 1 0,06 0,16
Factor de Peso 0,3 0,3 0,3
Calificación Ponderada 0,3 0,02 0,05
Valor Alta Media Media
Calificación 1 0,6 0,6
Manejabilidad
Factor de Peso 0,2 0,2 0,2
Calificación Ponderada 0,2 0,12 0,12
Valor Alta Alta Alta
Facilidad de Calificación 1 1 1
Consecución Factor de Peso 0,2 0,2 0,2
Calificación Ponderada 0,2 0,2 0,2
Valor Bajo Media Media
Calificación 1 0,6 0,6
Bajo Costo
Factor de Peso 0,3 0,3 0,3
Calificación Ponderada 0,3 0,18 0,18
Total 1 0,52 0,55
5.3.2 Contenedor
El contenedor esta en contacto directo con el medio del que se extrae calor o al
que se entrega. Por ello, su resistencia al ambiente de operación es importante.
En este caso el contenedor debe ser resistente a la corrosión, ya que el
condensador esta sumergido en el agua que se pretende calentar, la cual puede
tener compuestos oxidantes.
111
La evaluación para la selección del material del contenedor incluyendo los criterios
mencionados y otros se presenta a continuación.
Conductividad (W/mK)
206 45 377
3
Densidad (kg/m )
2722 7860 8902
Valor
Calificación
Esfuerzo de cedencia Factor de Peso 0,2 0,2 0,2
Calificación Ponderada 0 0 0
112
El contenedor del tubo de calor se fabrica en tubería flexible de cobre con ½
pulgada de diámetro exterior.
113
Tabla 36 Matriz selección mecha
Malla Espiral Fibra de
CRITERIO DE CALIFICACIÓN Resorte
metálica labrada vidrio
Valor Media Alta Alta Media
Calificación 0,6 0,60 1,00 0,60
Permeabilidad Factor de Peso 0,125 0,125 0,125 0,125
Calificación
0,075 0,075 0,125 0,075
Ponderada
Valor Media Baja Baja Media
Calificación 0,6 0,1 0,1 0,6
Capilaridad Factor de Peso 0,5 0,5 0,5 0,5
Calificación
0,3 0,05 0 0,3
Ponderada
Valor Alta Alta Alta Baja
Calificación 1 1 1 0,1
Conductividad térmica Factor de Peso 0,1 0,1 0,1 0,1
Calificación
0,1 0,1 0,1 0,01
Ponderada
Valor Media Baja Baja Media
Calificación 0,6 0,6 0,6 0,6
Facilidad de Consecución Factor de Peso 0,2 0,2 0,2 0,2
Calificación
0,12 0,12 0,12 0,12
Ponderada
Valor Media Baja Baja Baja
Calificación 0,6 0,1 0,1 0,1
Facilidad Manufactura /
Instalación Factor de Peso 0,2 0,2 0,2 0,2
Calificación
0,12 0,02 0,02 0,02
Ponderada
Valor Medio Alto Medio Alto
Calificación 0,6 0,1 0,6 0,1
Bajo Costo Factor de Peso 0,2 0,2 0,2 0,2
Calificación
0,12 0,02 0,12 0,02
Ponderada
Total 0,835 0,385 0,485 0,545
5.4.1 Geometría
La geometría del tubo esta definida por el diseño global del prototipo de calentador
solar y por la facilidad de consecución.
114
Tabla 37 Dimensiones tubo de calor
Dimensión Nombre Valor Unidades Criterio
de Diámetro externo 0.5 pulg Dimensiones tubería
tubo de calor flexible de cobre comercial
di Diámetro interno 11 mm Dimensiones tubería
contenedor flexible de cobre comercial
Le Longitud 1 m Área de colección de
Evaporador 1m*1m
Lc Longitud 0.1 m Diámetro interno tubería
condensador para el intercambiador
La Longitud zona 0.2 m Posibilitar reparaciones
adiabática futuras.
∆P
t= ri FS
S yp
(2)
115
Del análisis anterior se concluye que la tubería flexible de cobre de ½" de diámetro
exterior, tiene el espesor suficiente para soportar la diferencia de presiones
máximas.
La cantidad de fluido a cargar dentro de los tubos debe ser tal que logre mojar la
malla y garantice la existencia de una pequeña piscina de fluido en la parte
inferior. Para calcular el volumen de agua absorbido por la malla se encuentra el
espesor de esta y se halla el volumen de un tubo concéntrico al contendor con ese
espesor.
ri
2dw
π 2
V f = 2π ⋅ d w ⋅ d i ⋅ L + ⋅ di ⋅ Lp (3)
4
116
Donde dw es el diámetro del alambre con el cual está tejida la malla, di es el
diámetro interno del contendor, L la longitud total del tubo de calor y Lp la longitud
de la piscina.
Los parámetros y resultados para el cálculo del volumen del fluido son los
siguientes.
Finalmente se decide cargar los tubos con un volumen de fluido un poco mayor
igual a 38ml, lo que equivale aproximadamente a una masa de 38g en las
condiciones de carga.
117
5.4.4 Límites de Operación
2 ⋅ π ⋅ d i − 4nd w
w=
n (4)
Ilustración 48 Estructura Malla
4dw
118
Tabla 40 Propiedades estructura capilar
PROPIEDAD ECUACIÓN
1
Número de malla N=
dw + w
d +w
Radio Capilar rc = w
2
1.05 ⋅ π ⋅ N ⋅ d w
Porosidad ε = 1−
4
dw ⋅ε 2
2
Permeabilidad K=
122(1 − ε ) 2
Conductividad Efectiva k [(k + k w ) − (1 − ε )(kl − k w )]
keff = l l
kl: Conductividad del líquido
kw: Conductividad material malla (kl + k w ) + (1 − ε )(kl − k w )
Los límites se evalúan para el tubo horizontal con una temperatura de operación
de 90ºC y se obtienen los siguientes resultados.
119
6 ASPECTOS IMPORTANTES EN LA CONSTRUCCIÓN DE LOS TUBOS DE
CALOR PARA EL PROTOTIPO DE CALENTADOR SOLAR
6.1 LIMPIEZA
120
• Llenado de los tubos con una solución de soda cáustica (NaOH) al 2%
masa. La solución permanece dentro de los tubos durante 15 minutos.
• Enjuague del interior de los tubos con abundante agua.
• Llenado de los tubos con una solución de ácido clorhídrico (HCl) al 12%
masa. La solución permanece dentro de los tubos durante 15 minutos.
• Enjuague del interior de los tubos con abundante agua.
121
de accesorios se diseña un dispositivo para lavado. El dispositivo permite llenar
completamente de solución tres tubos cada vez. Además facilita el enjuague,
porque dispone de entradas y salidas para hacer correr el agua. El uso del
dispositivo es seguro, ya que se controla la emisión de vapores del ácido y el
contacto con los reactivos es mínimo.
122
Tabla 42 Partes del dispositivo para el lavado de los tubos
ESQUEMA NOMBRE ACCESORIO CANTIDAD
Unión hembra – hembra
4
½” IPS
T con tres extremos
6
macho ½”
Adaptador ½” PVC 6
Válvula en T ½” 2
Tapón ½” IPS 1
Adaptador ½” para
2
conexión manguera
La malla también debe limpiarse por las mismas razones que los tubos, y para
garantizar el flujo a través de ella. La limpieza de la malla difiere con la de los
tubos de calor en la concentración de las soluciones utilizadas. La malla es más
susceptible al ataque de las sustancias utilizadas, por eso las soluciones tienen
bajas concentraciones. Los pasos para este procedimiento son:
123
Ilustración 53 Inmersión de la malla en solución
124
Ilustración 55 Proceso secado malla
Los capilares son tramos cortos de tubería de cobre con diámetro de 1/8 pulgada,
y los termopozos son tramos cortos con un extremo cerrado de tubería de cobre
de 3/8 pulgada. La limpieza de estos elementos, por su tamaño, es sencilla. Estos
se pueden sumergir fácilmente en un recipiente.
Las concentraciones y lo tiempos de permanencia son los mismo que para los
tubos de calor.
125
• Inmersión en una solución de soda cáustica (NaOH) al 2% masa durante
15 minutos.
6.2 ENSAMBLE
126
6.3 CIERRE
127
Ilustración 60 Proceso cierre tubo
6.4 SOLDADURA
128
Ilustración 62 Soldadura termopozos
6.5 CARGA
129
• Realizar vacío dentro del tubo.
130
6.6 SELLADO
Las válvulas de carga de los tubos no garantizan una conservación del vacío por
largos períodos, entonces se hace necesario retirarlas y sellar los tubos de
manera definitiva. Los pasos para el sellado son:
• Estrangular el capilar del tubo.
• Cortar el capilar con soldadura de arco eléctrico.
• Probar el tubo de calor.
• En caso que el proceso de carga no resultase exitoso, se retira el capilar
calentando la soldadura de nuevo y se instala uno nuevo.
131
7 EVALUACIÓN DEL PROTOTIPO DE CALENTADOR SOLAR CON TUBOS
DE CALOR
Las pruebas de desempeño técnico del prototipo requieren del uso de diferentes
instrumentos de medición. Las variables más importantes a registrar son: la
intensidad de radiación, las temperaturas y flujos másicos. En esta sección se
discute alrededor de los instrumentos utilizados para la medición de algunas de
estas variables.
132
que ingresa al sistema. Este instrumento utiliza una celda solar, aprovechando el
hecho que para bajas resistencias la corriente de la celda es directamente
proporcional a la energía que incide (SOLAREX CORP, 1979). Esta celda se
conecta de acuerdo con el circuito mostrado en la figura.
mA
+
1 kΩ
1Ω
-
133
Los puntos de temperatura que se utilizan para la calibración son tres:
134
Ilustración 67 Esquema prueba equilibrio de temperatura
Tc Te
4
Te - Tc (ºC)
0
0 100 200 300 400 500 600 700
-1
-2
Radiación Aproximada en W/m2
135
Los resultados de las pruebas permiten intuir que, para el tubo evaluado, la
diferencia de temperatura entre condensador y evaporador aumenta cuando la
intensidad de radiación solar aumenta. Al tubo le resulta más difícil mantener el
equilibrio cuando esta recibiendo más calor de la fuente.
TANQUE
136
la intensidad de radiación solar. Las temperaturas de los condensadores y el
tanque, en cambio, mantuvieron una pequeña tendencia de crecimiento durante
toda la prueba. Esto se debe a la alta inercia térmica que agrega el agua que se
pretende calentar.
120,0
100,0
80,0
T ºC
60,0
40,0
20,0
0,0
9:21:36 10:33:36 11:45:36 12:57:36 14:09:36
Hora
137
calor útil para la eficiencia se calcula encontrando el incremento en la energía
interna para la masa en una unidad de tiempo de una hora.
Eficiencia Real
tanque inicial
Temperatura
Temperatura
Temperatura
Hora Inicial
Hora Final
Factor del
Radiación
placa (ºC)
Eficiencia
(ºCm2/W)
Calor Útil
final (ºC)
colector
Teórica
(W/m2)
(W/m2)
(ºC)
10:21:02 11:20:17 418 23,8 24,1 67,9 24,1 0,096 5,8% 2,5%
11:20:17 12:26:52 610 24,1 25,4 77,4 93,0 0,081 15,2% 14,8%
12:26:52 13:21:03 738 25,4 25,4 96,6 0,0 0,093 0,0% 4,7%
13:21:03 14:14:04 648 25,4 25,8 87,3 35,9 0,091 5,5% 5,9%
14,0%
12,0%
10,0%
Eficiencia
8,0%
6,0%
4,0%
2,0%
0,0%
0,073 0,078 0,083 0,088 0,093 0,098
(Tp-Ta)/Qi (m2K/W)
138
La eficiencia del prototipo de calentador solar es algo baja, ya que se está
elevando la temperatura de la placa. Existe una deficiencia en la manera como se
retira el calor de la placa. En promedio, durante la prueba, la eficiencia fue del
7.6% y se alcanzaron eficiencias máximas del 17.5%.
90,00
80,00
70,00
60,00
T (ºC)
50,00
40,00
30,00
20,00
10,00
0,00
0 200 400 600 800 1000 1200 1400
Evaporador Zona adiabática Cond
Longitud (mm)
120 W/m2 207 W/m2 312 W/m2 404 W/m2 514 W/m2 608 W/m2
139
7.4 PRUEBA CIRCULACIÓN FORZADA
Se realizan pruebas con circulación forzada con el fin de aumentar los coeficientes
de transferencia en la parte superior y así facilitar la evacuación del calor de los
tubos.
TANQUE
La captura de datos tiene una duración de tres horas, en este tiempo las
temperaturas tienen el siguiente comportamiento. La temperatura media en el
140
tanque pasa de 24,4°C a 34,8°C, esto es equivalente para los 67kg a 2917,50kJ
de energía y en promedio durante la prueba a un flujo de calor de 301,2W. Este
valor es alto si se tiene en cuenta que la radiación durante la prueba fue de
aproximadamente 680W/m2, lo que se traduce en una eficiencia del 44%. Sin
embargo existe cierta distorsión, ya que la potencia de la bomba es más grande
de la requerida y una parte considerable de esta se traduce en calor que va al
agua.
100,0
80,0
Temperatura (°C)
60,0
40,0
20,0
0,0
10:40:48 11:09:36 11:38:24 12:07:12 12:36:00 13:04:48 13:33:36 14:02:24
Hora
141
Tabla 44 Eficiencia en prueba con circulación forzada
Caudal 3 LPM
Flujo Másico 0,05 kg/s
Temperatura inferior
superior tanque (°C)
Temperatura Placa
Eficiencia Teórica
ambiente (°C)
Temperatura
Temperatura
tanque (°C)
Hora Inicial
Hora Final
Eficiencia
Calor útil
(°C)
10:45:05 11:01:30 627,1 27,6 25,7 24,2 75,0 0,08 324,7 51,8% 16,3%
11:01:30 11:30:04 780,5 26,8 26,3 27,2 94,2 0,09 -188,3 -24,1% 6,7%
11:30:04 12:00:10 769,7 27,3 29,6 29,7 102,4 0,10 -12,5 -1,6% -3,2%
12:00:10 12:30:19 773,8 27,5 32,2 31,2 105,0 0,10 203,7 26,3% -5,5%
12:30:19 13:00:38 478,7 30,1 34,9 34,0 90,1 0,13 197,4 41,2% -27,8%
13:00:38 13:32:39 647,4 30,9 36,0 35,9 86,7 0,09 21,0 3,2% 6,8%
80,0%
70,0%
60,0%
Eficiencia
50,0%
40,0%
30,0%
20,0%
10,0%
0,0%
0,00 0,02 0,04 0,06 0,08 0,10 0,12 0,14
(Tp-Ta)/Qi
142
parte, a la distorsión impuesta por la bomba. En cuanto a eficiencia los datos no
son concluyentes.
100,0
80,0
Temperatura (°C)
60,0
40,0
20,0
0,0
0 200 400 600 800 1000 1200 1400
Cond
Evaporador Zona adiabática
332 W/m2 500 W/m2 595 W/m2 783 W/m2 Longitud (mm)
143
Ilustración 77 Esquema prueba tubo con diferentes inclinaciones
Intercambiador
144
Ilustración 79 Distribución de temperatura con una inclinación de 12°
110
100
90
80
Temperatura (ºC)
70
60
50
40
30
20
10
0
0 200 400 600 800 1000 1200 1400
Evaporador Zona adiabática Cond
110
100
90
80
70
60
T (ºC)
50
40
30
20
10
0
0 200 400 600 800 1000 1200 1400
Evaporador Zona adiabática Cond
145
Ilustración 81 Distribución de temperaturas con una inclinación de 25°
100,0
90,0
80,0
70,0
60,0
T (ºC)
50,0
40,0
30,0
20,0
10,0
0,0
0 200 400 600 800 1000 1200 1400
Evaporador Zona adiabática Cond
Después de la carga se realiza una purga para retirar los gases no condensables y
se hace un chequeo para evaluar su efectividad. Es decir, que en el momento de
sellar los tubos no hay presencia de gases no condensables. Estos se generan,
146
posiblemente, en reacciones dentro del tubo de la malla con algunos iones diluidos
en el agua.
No se conoce ninguna experiencia previa del uso de este tipo de malla, por lo cual
en el momento de construcción se deja un testigo de esta sumergida en agua. Con
el paso del tiempo el agua se enturbia con partículas blancas, manifestación de la
degradación de la malla.
7.6 COSTOS
Los costos finales de los materiales necesarios para la construcción del prototipo
se especifican en la tabla. Como puede observarse el ítem que tiene más peso
sobre el costo total es el cobre, seguido por la madera plástica y el material para
los aislamientos. Los elementos más costosos del colector son los tubos de calor,
el soporte y los aislamientos. Un incremento en el precio de los materiales para su
fabricación tiene una notoria influencia en el costo total.
El costo unitario de los materiales de los tubos de calor es cercano a $13.000 por
tubo. Adicionalmente, está el costo de los procesos de manufactura, en la
fabricación de los tubos se invierten aproximadamente 30h lo que es equivalente a
3,3h por tubo. Esta es una clara manifestación de las complicaciones que tiene la
147
tecnología de fabricación de tubos de calor y que esta se desarrolla aún a nivel
experimental.
148
8 CONCLUSIONES
Se desarrolló un equipo que combina la tecnología de los tubos de calor con la del
uso de la energía solar. El proceso consistió: de la adquisición de conocimiento
general sobre ambos temas, la identificación de los requerimientos, el diseño
metódico, la construcción y finalmente una breve evaluación del prototipo
resultante.
Se diseñó un prototipo de calentador solar con tubos de calor que tiene un área de
colección de un metro cuadrado y calienta una masa de agua de 67kg. En el panel
se ubican 9 tubos de calor con una longitud de 1.3m, los límites de transferencia
de calor para los tubos diseñados son todos superiores a 190W y cada tubo debe
transferir 50W. Los modelos matemáticos utilizados indican que el panel tendría
una eficiencia teórica del 40%.
149
En el prototipo se utilizaron varios materiales plásticos, por la facilidad para su
procesamiento y costo moderado. Se hicieron las respectivas consideraciones
sobre la coherencia de sus propiedades con las condiciones de operación y se
concluyó que esta combinación es innovadora, técnicamente viable y
económicamente favorable.
150
Se realizaron pruebas en el prototipo para evaluar su desempeño. En la prueba
con convección natural el equipo tuvo una eficiencia del 7.6% en promedio. Esta
eficiencia es baja considerando que para las mismas condiciones la eficiencia
teórica es del 10%, lo que implica que las pérdidas de calor son mayores a las
consideradas inicialmente y / o una capacidad de transporte térmico baja.
151
El prototipo construido actualmente no es competitivo comparado con valores de
rendimiento de calentadores solares de agua comerciales.
152
9 RECOMEDACIONES Y SUGERENCIAS
Realizar pruebas de resistencia del extremo del tubo que se cierra laminando el
material. Debe determinarse con exactitud el espesor final y la existencia o no de
concentradores de esfuerzo. Para tubos que tengan temperaturas de operación
más altas, este podría ser un punto de falla.
153
pero fabricadas de acero inoxidable. En este aspecto es necesario realizar una
exploración para determinar la facilidad de consecución.
154
10 BIBLIOGRAFÍA
10.1 CLÁSICA
BEJAN, Adrian y KRAUS, Allan. “Heat transfer handbook”. Hoboken: John Wiley &
Sons, 2003. ISBN: 0-471-39015-1
155
FAGHRI, Amir. “Heat Pipe Science and Technology”. Londres: Taylor & Francis,
1995. ISBN 1-56032-383-3
LUNDE, Peter. “Solar Thermal Engineering, Space Heating and Hot Water
Systems”. Nueva York: John Wiley & Sons, 1980. ISBN: 0-471-08233-3.
156
RODRIGUEZ, Humberto y GONZALEZ, Fabio. “Manual de radiación solar en
Colombia”. Bogotá:, 1992. ISBN 958-95506-0-6.
SILVERSTEIN, Calvin C. “Design and Technology of Heat Pipes for Cooling and
Heat Exchange”. Nueva York: Taylor & Francis, 1992. ISBN:0-89116-859-1.
WIEDER, Sol. “An introduction to Solar Energy for Scientist and Engineers”. Nueva
York: John Wiley & Sons, 1982. ISBN: 0-471-06048-8.
10.2 INTERNET
ENCARTA@
Enciclopedia Virtual Encarta. MSN, 2005 http://es.encarta.msn.com/
GARBITEK@
Fabricante de calentadores solares. Visitado en marzo de 2007.
http://www.garbitek.com
GREEN PAGES@
157
Uso racional de los recursos naturales – Visitado en agosto de 2006.
http://www.eco-web.com
INFOAGRO@
Información sobre montaje de invernaderos – Visitado en abril 2007.
http://canales.hoy.es/canalagro/datos/industria_auxiliar/plasticos5.htm
ITESO@
Instituto tecnológico y de estudio superiores de occidente. Notas sobre el curso de
energía solar – Visitado en marzo 2007. www.pti.iteso.mx
MYSOLAR@
Fabricante de productos que utilizan energía solar - Visitado en septiembre de
2006. http://www.mysolar.com/
THERMOMAX@
Fabricante de aplicaciones solares – Visitado en marzo de 2007.
http://www.thermomax.com
158
UNIVERSIDAD DE CANTABRIA@
Procesos termosolares en baja, media y alta temperatura – Visitado en marzo
2007. http://www.termica.webhop.info/
WIKIPEDIA@
Enciclopedia Virtual – Agosto de 2006. http://es.wikipedia.org
159
ANEXO PLANOS
160