PFC Victor Salguero Fernandez

UNIVERSIDAD CARLOS III DE MADRID
ESCUELA POLITÉCNICA SUPERIOR
DEPARTAMENTO DE INGENIERÍA TÉRMICA Y DE FLUIDOS
Valoración de colectores solares térmicos

planos de vacío para la producción de frío
PROYECTO FIN DE CARRERA
INGENIERÍA TÉCNICA INDUSTRIAL MECÁNICA
AUTOR: VÍCTOR SALGUERO FERNÁNDEZ
TUTORA: Dª MARÍA DEL CARMEN VENEGAS BERNAL
LEGANÉS, SEPTIEMBRE, 2008

Agradecimientos
En este pequeño texto, intentaré nombrar a todas las personas que me han
ayudado a la realización de este proyecto, ya que sin ellas no hubiera sido posible la
finalización del mismo.
En primer lugar a María Venegas Bernal, gracias por haberme introducido en la

ingeniería térmica y en las energías renovables ya que he encontrado la dirección de
mi profesión en la cual deseo trabajar, también por tu comprensión y por tu atención en
aquellos momentos en los no llegaba al nivel necesario en la realización del proyecto.
También se ha de agradecer la colaboración desinteresada de María del Carmen
Rodríguez, gracias a la cual he podido seguir adelante cuando María no podía
guiarme.
A mis compañeros David (Toxas), Nacho (Naxa Pop), Richi (Gusiluz), María
(Mari Paca), Pablo (recupérate pronto), Viki y Alfredo…muchas horas al lado vuestro,
codo con codo, tanto en los disgustos como en los buenos momentos, prácticas
exámenes….deciros que cada uno de vosotros me habéis aportado mucho en estos años
(risas, puntos de vista distintos, compañerismo, apoyo, amistad, decisión…) y que sin
vosotros la carrera hubiera sido de otra manera. Eh, que no me olvido de los
superiores: Víctor Sergio (Zidane), Esteban, Javi, y Manu, el año que viene con
vosotros de nuevo, gracias por abrir el camino jejeje.
A mis amigas Yolanda, Rebe, Laura (pimaaa), me habéis acompañado toda mi

vida, me habéis apoyado, animado, abrazado, noches sin dormir, fiestas, salidas,
recuerdos con vosotras…deciros que este es vuestro granito de arena en este proyecto,
gracias por haber estado ahí, muaks!! Al grupo de” ingeniero químicos (de palo)”:
Lau, Miriam y Piti, me habéis abierto vuestro corazoncito y me he sentido como uno
más a vuestro lado, gracias.
A mi tridente: María, Pablo (Melenudo) y Mari Carmen, personas en mi vida de

las cuales cada día sigo apreciando el día que os conocí, cada día me doy cuenta de los
afortunado que soy por teneros a mi lado, no me salen las palabras para describiros…

A mi Peque, casi cinco años a tu lado, llenos de amor, cariño, atención, dulzura,
de besos y caricias, de viajes para verte, de recuerdos y momentos a tu lado, por fin te
puedo decir he acabado y dedicarte el esfuerzo que he realizado, gracias porque has
estado a mi lado, TQ gorda.
Los últimos pero los más importantes, mi familia: Papá, Mamá y Memo. Papá y
Mamá gracias por la educación recibida de vosotros, por la oportunidad que me habéis
dado para estudiar, por vuestro apoyo, seguridad y esperanza porque terminara la
carrera, pues este es el resumen de ella, ya puedo decir que soy Ingeniero y todo por
vuestra dedicación a nosotros. Memo por aconsejarme en todo momento y por
compartir todo lo que eres conmigo. OS QUIERO MUCHO.
Al Jefe, por haberme puesto en mi camino a todas las personas que tengo,
gracias.

ÍNDICE

“Valoración de colectores solares térmicos planos de vacío para la producción de frío”
ÍNDICE
0 Introducción…………………………………………………………………………………1
I Importancia del proyecto…………………………………………………………1
II Objetivos……………………………………………………………………………5
III Descripción por capítulos………………………………………………………6
1 La máquina de absorción……………………………………………………….............8
1.1 Fundamentos termodinámicos de la máquina de absorción……………9
1.2 Componentes de la máquina ideal y su COP……………………………..13
1.3 Fluidos de trabajo………………………………………………………………18
1.4 Tipos de máquina de absorción……………………………………………..20
1.5 Consideraciones ambientales de la máquina de absorción................23
1.6 Descripción de una planta solar…………………………………………….24
1.7 Actualidad de la planta solar y la máquina de absorción………………25
2. El colector solar…………………………………………………………………………..35
2.1 Colectores solares planos……………………………………………………35
2.1.1 Descripción y funcionamiento de un colector plano…………36
2.1.2 Colector solar plano de vacío……………………………............40
2.1.3 Ecuaciones energéticas de un colector solar plano………….41
1)Rendimiento del colector solar plano……………………….41
2)Ecuación fundamental de colector…………………………..42
3)Pérdidas a través de la cubierta……………………………...43
4)Pérdidas a través de la parte posterior y laterales………..48
“Departamento de Ingeniería Térmica y Fluidos” pág:i

5)Pérdidas totales…………………………………………...........50
6)Calor útil transferido por el colector………………………...50
7)Rendimiento de colectores en régimen

permanente…………………………………………………………51
8) Pérdidas de calor por transferencia de calor de la

instalación…….………………………………….………………...53
9)Calor absorbido por el colector solar……………………….55
10)Calor absorbido por los componentes de la

instalación…………………………………………………………..55
11) Energía solar disponible…………………………………….56
2.2 Colectores concentradores…………………………………………………..58
2.2.1 Descripción y funcionamiento de un colector

concentrador………………………………………………………………..60
2.2.2 Tipos de colectores solares concentradores…………............60
2.2.2.1 Tubos evacuadores………………………………………60
2.2.2.2 Colectores parabólicos………………………………….66
3. Datos experimentales…………………………………………………………..............68
3.1 Descripción de la instalación………………………………………………68
3.1.1 Colectores solares planos de vacío……………………………...70
3.1.2 Intercambiador de calor……………………………………………72
3.1.3 Máquina de absorción LiBr-H2O…………………………............73
3.1.4 Tanque de almacenamiento……………………………………….76
3.1.5 Circuitos hidráulicos de la instalación…………………………..77
A) Circuito primario……………………………………………….77
“Departamento de Ingeniería Térmica y Fluidos” pág:ii

B) Circuito secundario……………………………………………78
C) Circuito terciario……………………………………………….79
D) Circuito cuaternario……………………………………………80
3.1.6 Fan-coils………………………………………………………………81
3.1.7 Fluidos de trabajo…………………………………………………...81
3.1.8 Sistema de regulación, bombeo y adquisición de datos…….82
3.2 Instalación solar térmica…………………………………………………….85
3.3 Máquina de absorción……………………………………………………….99
4. Simulación……………………………………………………………………………….100
4.1 Desarrollo del modelo………………………………………………………..100
4.1.1 Adquisición………………………………………………………….100
4.1.2 Caract Col1………………………………………………………….105
4.1.3 Instalación…………………………………………………………..110
4.2Resultados del modelo y comparación con los datos

experimentales……………………………………………………………………121
4.2.1 Temperaturas……………………………………………………….121
4.2.2 Rendimiento de la instalación…………………………………...125
4.2.3 Calor recibido por los colectores y pérdidas térmicas de la

instalación………………………………………………………………….129
4.3 Influencia de la superficie solar y del volumen del tanque de sobre la

instalación………………………………………………………………………….132
1) Variación de la superficie solar y volumen del tanque

constante…………………………………………………………………133
“Departamento de Ingeniería Térmica y Fluidos” pág:iii

2) Variación del volumen del tanque y superficie solar

constante…………………………………………………………………134
5. Estudio económico……………………………………………………………………..135
6. Estudio ambiental………………………………………………………………………140
6.1 Análisis de impacto sobre el aire………………………………………...141
6.2 Análisis de impacto sobre el agua………………………………...........142
6.3 Ruido…………………………………………………………………………..143
6.4 Cantidad no emitida de CO2………………………………………...........143
7. Conclusiones del proyecto…………………………………………………………...145
8.Bibliografía………………………………………………………………………………..147
“Departamento de Ingeniería Térmica y Fluidos” pág:iv

ÍNDICE DE FIGURAS
1. La máquina de absorción.
1.1 Modelo termodinámico de la máquina de absorción……………………..9
1.2 Máquina de absorción…………………………………………………………13
1.3 Generador………………………………………………………………………..16
1.4 Máquina de absorción de simple efecto…………………………………...20
1.5 Planta solar acoplada a una máquina de absorción……………………..24
1.6 Coste energético frente a capacidad frigorífica…………………………..27
1.7 Sistemas instalados en Europa……………………………………………...28
1.8 Tipos de máquinas empleadas en Europa………………………………...29
1.9 Esquema de la instalación……………………………………………………30
1.10 Colectores solares……………………………………………………………31
1.11 Croquis de la instalación……………………………………………………32
1.12 (a) Colectores solares, (b) calentador de agua, (c) depósitos y (d)
máquina de absorción……………………………………………………………..32
1.13 Colectores solares……………………………………………………………33
1.14 Máquinas de absorción……………………………………………………...33
1.15 Distribución de los colectores solares……………………………………34
1.16 Detalle de los colectores solares…………………………………………..34
2) El colector solar
2.1 Sistema abierto o por goteo………………………………………………….37
“Departamento de Ingeniería Térmica y Fluidos” pág:v

2.2 Configuraciones del sistema de soldadura………………………………38
2.3 Partes de un colector plano………………………………………………….39
2.4 Colector solar plano de vacío………………………………………………..40
2.5 Colector plano con dos cubiertas…………………………………………..43
2.6 Curvas de normalización……………………………………………………..53
2.7 Energía solar disponible………………………………………………………58
2.8 Esquema de tubos evacuadores…………………………………………….60
2.9 Esquema de captación de un sistema CPC……………………………….61
2.10 Tubo evacuador simple……………………………………………………...62
2.11 Tubo de vacío de flujo directo……………………………………………...63
2.12 Esquema de un tubo Heat Pipe…………………………………………….64
2.13 Funcionamiento de un tubo Heat Pipe……………………………………65
2.14. Colector Cilíndrico Parabólico…………………………………………….66
2.15 Reflector y refractor Fresnel………………………………………………..67
3) Datos experimentales.
3.1 Esquema de la instalación……………………………………………………69
3.2 Colector solar plano de vacío 400V plus, en el tejado…………………..70
3.3 Medidas del intercambiador de calor……………………………………….72
3.4 Intercambiador de calor Alfa Laval M6MF…………………………………72
3.5 Esquema de la máquina de absorción……………………………………..73
3.6 Máquina de absorción…………………………………………………………74
“Departamento de Ingeniería Térmica y Fluidos” pág:vi

3.7 Tanques de almacenamiento………………………………………………...76
3.8 Vista de la aplicación que muestra los valores registrados en la

instalación meteorológica, actualizados cada minuto………………………83
3.9 Radiación global frente al tiempo…………………………………………...95
3.10 Rendimiento de los colectores frente al tiempo………………………...96
3.11 Rendimiento de los colectores frente a la temperatura de salida del

colector……………………………………………………………………………….97
3.12 Calor útil frente a temperatura de salida de los colectores…………..98
3.13 Radiación global diaria recibida entre 30/06 y 11/07…………………...99
4) Simulación.
4.1 Esquema de la tubería……………………………………………………….113
4.2 Evolución de temperaturas simulación (superior) y experimental

(inferior) del día 30/06/2007……………………………………………………...121

(inferior) del día 05/07/2007……………………………………………………...122

(inferior) del día 11/07/2007……………………………………………………...123
4.5 Evolución temporal de ηs del día 30/06/2007…………………………….125
4.8 Evolución temporal de Qpl y Pt1 del día 30/06/2007……………………..130
4.9 Evolución temporal de Qpl y Pt1 del día 05/07/2007……………………..130
4.10 Evolución temporal de Qpl y Pt1 del día 11/07/2007……………………131
“Departamento de Ingeniería Térmica y Fluidos” pág:vii

5) Estudio económico.
5.1 Desglose del presupuesto de la instalación……………………………..139
“Departamento de Ingeniería Térmica y Fluidos” pág:viii

ÍNDICE DE TABLAS
1. La máquina de absorción.
1.1 COP y Rango de temperaturas de funcionamiento………………………22
1.2 Máquinas de absorción de baja potencia………………………………….26
3) Datos experimentales.
3.1 Especificaciones técnicas del panel 400V plus…………………………..71
3.2 Propiedades térmicas de los componentes del panel…………………..71
3.3 Características de las placas………………………………………………...72
3.4 Dimensiones del intercambiador……………………………………………72
3.5 Especificaciones de la máquina de absorción……………………………75
3.6 Dimensiones de los tanques de almacenamiento………………………..76
3.7 Propiedades de los materiales que conforman los tanques.………….77
3.8 Dimensiones de las tuberías del circuito primario………………………77
3.9 Propiedades de los materiales del circuito primario…………………….77
3.10 Dimensiones de las tuberías del circuito secundario………………….78
3.11 Propiedades de los materiales del circuito secundario……………….79
3.12 Dimensiones de las tuberías del circuito terciario……………………..79
3.13 Propiedades de los materiales del circuito terciario…………………...80
3.14 Dimensiones de las tuberías del circuito cuaternario…………………80
3.15 Propiedades de los materiales del circuito cuaternario………………81
3.16 Propiedades de los fluidos de la instalación…………………………….82
“Departamento de Ingeniería Térmica y Fluidos” pág:ix

3.17 Bombas de la instalación……………………………………………………82
3.18 Funcionamiento de la instalación………………………………………….84
3.19 Radiación global y temperatura de entrada y salida del agua que

circula por los colectores…………………………………………………………86
3.20 Calor útil y rendimiento de los colectores……………………………….91
3.21 Valores medios de Qu y ηc…………………………………………………..95
4) Simulación.
4.1 Valores geométricos de los colectores…………………………………..106
4.2 Valores de ηsT………….……………………………………………………...128
4.3 Resultados variación de la superficie solar……………………………..133
4.4 Resultados variación del volumen del tanque…………………………..134
5) Estudio económico.
5.1 Coste de la instalación………………………………………………………135
5.2 Precio de las bombas..………………………………………………………136
5.3 Coste de electricidad…………………………………………………………137
5.4 Gastos de explotación.………………………………………………………138
5.5 Desglose del presupuesto de la instalación……………………………..138
6) Estudio medioambiental.
6.1 Impactos sobre el aire………………….……………………………………141
6.2 Impactos sobre el agua………...……………………………………………142
6.3 Emisiones de CO2…………………………………………………………….143
“Departamento de Ingeniería Térmica y Fluidos” pág:x

0. INTRODUCCIÓN

0. Introducción.
I. Importancia del proyecto.
La contribución de la energía solar térmica al consumo energético mundial

sigue siendo muy escasa todavía, pese a que empiezan a percibirse ciertos
síntomas de cambio que permiten ser más optimistas de cara al futuro. Al
creciente interés de los ciudadanos por este tipo de soluciones hay que sumar
las ayudas e incentivos que se han puesto en marcha en muchos países del
mundo y la reducción de precios de los captadores solares en algunos
mercados, especialmente activos, como China o Japón. Una situación que
pone de manifiesto que estamos ante una tecnología madura que ha
experimentado un significativo avance durante los últimos años [36].
En la actualidad la capacidad de energía solar instalada en el mundo

supera a la de otras renovables con altos índices de desarrollo, como es el
caso de la energía eólica. Con una potencia instalada de 98,4 GW térmicos a
finales de 2004 (Datos del Solar Heat Worldwide 2004, considerando 41 países
que representan el 57% de la población mundial y el 85-90% del mercado
mundial de solar térmica), la solar térmica ha alcanzado unos niveles de
popularidad impensables hace tan solo unos años. Y no exclusivamente por lo
que a la producción de agua caliente se refiere, sino también en cuanto a la
calefacción de viviendas [36].
Europa representa tan solo el 9% del mercado mundial de energía solar

térmica con una potencia instalada de 10.000 MWth (MW térmico) a finales de
2004, o lo que es lo mismo, un total de 14 millones de m2 de captadores
solares en funcionamiento. El impulso que ha recibido esta industria durante los
últimos años, es lo que ha permitido dar un paso firme en el objetivo común de
alcanzar los 100 millones de m2 de superficie instalada que se pretenden
conseguir en el horizonte de 2010 [36].
“Departamento de Ingeniería Térmica y Fluidos” pág:1

España es el cuarto país europeo en el aprovechamiento de la energía

solar térmica, por delante de países como Italia, Francia o Gran Bretaña. Con
un 6% del total del mercado europeo, nuestro país ha alcanzado la madurez
tecnológica y comercial tras más de 20 años de experiencia [36].
Nuestro país aún se encuentra lejos de los objetivos nacionales fijados en

el Plan de Energías Renovables (PER), que plantea alcanzar una superficie
instalada de 4,9 millones de metros cuadrados para el año 2010. Para ello, la
entrada en vigor del Código Técnico de la Edificación, que obliga a instalar un
aporte de energía solar para agua caliente en todas las viviendas de nueva
construcción, junto a las medidas ya puestas en marcha con anterioridad,
darán un impulso definitivo a un mercado con excelentes perspectivas a medio
y largo plazo [36].
En la actualidad, el principal cliente de energía solar en España es el

usuario particular, que solicita la instalación de captadores solares de baja
temperatura para el consumo de agua caliente sanitaria. En segundo lugar se
encuentran los hoteles y restaurantes, en los que existe creciente interés por
este tipo de soluciones energéticas. Además de estos dos grupos de
consumidores, que son los que más aportan al total de la superficie instalada
en España, en general, se puede decir que existen buenos ejemplos en
múltiples sectores y para todo tipo de aplicaciones posibles, pudiendo
mencionar las instalaciones en centros educativos, centros deportivos, centros
sanitarios, albergues, campings… [36]
La energía solar térmica aprovecha la radiación del Sol para calentar un

fluido que, por lo general, suele ser agua o aire. La capacidad de transformar
los rayos solares en calor es, precisamente, el principio elemental en el que se
basa esta fuente de energía renovable [36].

Las aplicaciones de este tipo de energía se muestran a continuación [36]:
a) Producción de agua caliente sanitaria (ACS): es una alternativa

completamente madura y rentable. Esto es debido a que los
niveles de temperaturas que se precisan alcanzar (normalmente
entre 40 y 45 °C), coinciden con los más adecuados para el buen
funcionamiento de los sistemas solares estándar que se
comercializan en el mercado. Además, hacemos referencia a una
aplicación que debe satisfacer a lo largo de todo el año, por lo
que la inversión en el sistema solar se rentabilizará más
rápidamente que en el caso de otros usos solares, como la
calefacción, que sólo tienen utilidad durante los meses fríos.
Con los sistemas de energía solar térmica hoy en día podemos
cubrir el 100% de la demanda de agua caliente durante el verano
y del 50 al 80% del total a lo largo del año; un porcentaje que
puede ser superior en zonas con muchas horas de sol al año,
como por ejemplo el sur de España.
b) Sistemas de calefacción: se cubre parcialmente la demanda de

calefacción en una vivienda y se puede llegar a un ahorro
energético de un 25%.
c) Climatización de piscinas: La climatización del agua para piscinas
constituye otra aplicación interesante de la energía solar, tanto si
se trata de instalaciones cubiertas como a la intemperie. Estas
últimas merecen especial atención al existir en gran número y al
conseguir resultados más que satisfactorios con sistemas
sencillos y baratos.
De hecho, resulta bastante económico lograr una temperatura
estable y placentera en piscinas al aire libre. En primer lugar
porque, al circular el agua de la piscina directamente por los
captadores solares, no es necesario utilizar ningún tipo de
intercambiador de calor ni de sistema de acumulación. Y en

segundo lugar, porque la temperatura de trabajo suele ser tan

baja (en torno a los 30 °C) que permite prescindir de cubiertas,
carcasas o cualquier otro tipo de material aislante. De esta
manera, se consigue reducir el precio del captador sin excesivo
prejuicio en su rendimiento.
d) Refrigeración en edificios: el aprovechamiento de la energía solar

para producir frío es una de las aplicaciones térmicas con mayor
futuro, pues las épocas en las que más se necesita enfriar el
espacio coinciden con las que se disfruta de mayor radiación
solar. Además, esta alternativa a los sistemas de refrigeración
convencionales es doblemente atractiva porque permite
aprovechar las instalaciones solares durante todo el año,
empleándolas en invierno para la calefacción y en verano para la
producción de frío.
De las diversas fórmulas de aprovechar el calor solar para

acondicionar térmicamente un ambiente, la más viable en
términos de coste de la inversión y ahorro de energía es la
constituida por el sistema de refrigeración por absorción, utilizada
en el 60% de los casos. El funcionamiento de estos equipos se
basa en la capacidad de determinadas sustancias para absorber
un fluido refrigerante. Como absorbentes se utilizan
principalmente el amoniaco o el bromuro de litio, mientras que
como líquido refrigerante es el agua el más recomendado.

Los beneficios ambientales y económicos de la energía solar térmica, son

numerosos, por lo que se destacan los más importantes [36]:
- No emite gases contaminantes perjudiciales para la salud.

- No emite gases de efecto invernadero que provocan el cambio climático.
- No produce ningún tipo de desperdicio o residuo peligroso de difícil
eliminación.
- Su impacto sobre el medio ambiente es mínimo, y de producirse alguno
ocurre exclusivamente durante la fase de fabricación de los equipos.
- Este tipo de instalaciones no dejan huella ecológica cuando finaliza el
período de explotación.
- Es una energía que no corre peligro de agotarse a medio plazo, puesto
que su fuente productora es el Sol.
- No requiere costosos trabajos de extracción, transporte o
almacenamiento.
- Desde el mismo momento en que pongamos en marcha nuestra
instalación solar, la factura del gas o la electricidad destinada a la
producción de agua caliente sanitaria bajará. Esto se traduce en ahorros
medios de entre unos 75 a 150 euros al año en una economía familiar, en
función del combustible que se sustituya.
II Objetivos del proyecto.
Este proyecto se ha realizado parcialmente con una beca de colaboración

con cargo al proyecto de investigación “Técnicas innovadoras para la
producción eficiente de frío con energía solar y máquinas de absorción”
CCG06-UC3M/ENE-0689. Los objetivos a cumplir son:
- Revisión bibliográfica sobre el estado actual de la tecnología de energía

solar térmica de baja temperatura y sus aplicaciones para la producción
de frío.

- Análisis de los datos experimentales de los colectores planos de vacío y

de la máquina de absorción YAZAKI WFC-10.
- Adaptación de un modelo de una instalación de frío solar (planta solar y
máquina de absorción) empleando el calor procedente de los colectores
solares planos y analizando la influencia de la superficie solar y volumen
del tanque de almacenamiento.
- Realización de un estudio económico y medioambiental de la instalación.
III Descripción por capítulos.
Capítulo 1: La máquina de absorción.
Se describe la modelización termodinámica de la máquina de absorción,

a continuación se detallan los elementos constituyentes de este equipo y las
ecuaciones que rigen los balances másicos y energéticos, para la consecución
de la expresión del coeficiente de operación (COP).
Se informa sobre las propiedades y cualidades que debe poseer las

disoluciones empleadas en este tipo de tecnología. Se tipifican las máquinas de
absorción existentes en el mercado y se valora los efectos ambientales en el
uso de ellas.
Finalmente se realiza una perspectiva actual sobre el “frío solar” tanto a

nivel europeo como nacional y se destacan algunas instalaciones existentes.
Capítulo 2: El colector solar.
Se realiza un estudio sobre los captadores solares: planos, tubos

evacuadores y parabólicos. Se describe su funcionamiento y los componentes
de cada uno de ellos.
Se destacan las ecuaciones de diseño y de pérdidas térmicas de los

colectores solares planos, ya que éstas son empleadas en el modelo la
instalación.

Capítulo 3: Datos experimentales.
Se describen los componentes de la instalación térmica solar y de la

máquina de absorción. Se analizan los datos experimentales obtenidos tanto
de los colectores solares planos como de la máquina.
Se obtiene un promedio de radiación global diaria recibida, rendimiento

diario de los colectores, el calor útil obtenido de la planta solar y se describe el
comportamiento de toda la instalación a lo largo de un período de
funcionamiento.
Capítulo 4: Simulación.
Se analiza el modelo matemático y se comparan los resultados obtenidos

con los experimentales.
Se caracteriza la influencia entre la superficie solar y el volumen del

tanque necesario para el funcionamiento de la instalación.
Capítulo 5: Estudio económico.
Se estima el coste de la instalación.
Capítulo 6: Estudio medioambiental.
Se realiza un análisis de los posibles impactos ambientales de la

instalación de frío solar y posibles medidas correctoras.
Capítulo 7: Conclusiones.
Se presenta de las conclusiones obtenidas al final del estudio.

1. LA MÁQUINA DE ABSORCIÓN

1. La máquina de absorción [1].
La máquina de absorción se define como una máquina frigorífica que

emplea como ciclo de trabajo el ciclo de absorción.
La absorción es un fenómeno químico que consiste en la separación de

uno o más componentes de una mezcla gaseosa con la ayuda de un solvente
líquido con el cual forma una disolución. Este proceso implica una
transferencia de masa del soluto A a través del gas B, que no se difunde y se
[8]
encuentra en reposo, hacia un líquido C, también en reposo . En el presente
proyecto el agua (H2O) es el soluto que se absorbe en bromuro de litio (LiBr).
La máquina de absorción fue creada por Nairne en 1774 y perfeccionada

en la segunda mitad del siglo XIX por Carré, sin sufrir ninguna modificación
hasta la década de 1950. Este hecho es debido a la producción de energía
eléctrica a precios muy reducidos.
Pero a partir de los años 70, debido al aumento del precio de la energía,
la industria ha incentivado y promovido el ahorro energético y más
concretamente en fuentes de calores residuales (p.ej: gases de escape de las
turbinas de gas empleados en calentar el agua en las plantas de
cogeneración). Estos focos de calor pueden ser empleados por las máquinas
de absorción, por lo que su comercialización aumentó en estos años sobre todo
en Estados Unidos (EEUU) y Japón en el sector del aire acondicionado.
Actualmente, la aplicación de esta tecnología se centra en cubrir la

demanda de refrigeración y aire acondicionado en los grandes edificios que es
un sector con una gran demanda energética [3] [5] [6].
“Departamento de Ingeniería Térmica y Fluidos” pág: 8

1.1 Fundamentos termodinámicos de la máquina de absorción.
El ciclo de absorción se constituye por una máquina térmica (máquina

directa de Carnot) y por una máquina frigorífica (máquina inversa de Carnot),
coincidiendo con el ciclo de compresión mecánica, representados en la Figura
1.1:
Foco caliente (Tg) Foco caliente (Tc)
Qg Qc
W
Qa
Qe
Foco frío (Te) Foco frío (Te)
Figura 1.1 Modelo termodinámico de la máquina de absorción[1].
El funcionamiento del modelo es el siguiente: la máquina térmica absorbe

energía calorífica, Qg, del foco caliente a temperatura Tg, transforma parte de
esta energía en trabajo mecánico W y transporta el resto, Qa, hasta el foco frío
a temperatura Te. La máquina frigorífica absorbe el trabajo mecánico W, y la
energía calorífica Qe de su foco frío a la temperatura Te, y la suma de ambos la
transporta hasta el foco caliente a la tempertatura Tc.

El parámetro termodinámico que define el comportamiento de la máquina

de absorción es el coeficiente de eficiencia energética o coeficiente de
operación (CEE o COP) [2] [3] [5] [6] [7], la expresión de este término es:
Qe
COP = ec [1.1]
Qg
[2] [3] [5] [6] [7]

Se procede al cálculo de COPmax o COPideal . Para ello se
asumen las siguientes hipótesis:
-. Las máquinas operan según un ciclo de Carnot (∆U=0).
-. Todos los procesos son reversibles (∆S=0).
-. Sistemas cerrados.
-. Signo de calor y trabajo:
Q>0: sistema recibe calor del exterior.
Q<0: sistema cede calor al exterior.
W>0: sistema realiza trabajo.
W<0: sistema recibe trabajo del exterior.
• Máquina térmica:
El rendimiento se define como:
W
η= ec [1.2]
Qg

Se aplica el primer principio de la termodinámica:
∆U = Q − W ⎫
⎪
Q = Qg − Qa ⎬ → W = Qg − Qa ec [1.3]
⎪
∆U = 0 (ciclo )⎭
Se aplica el segundo principio de la termodinámica:
δQ ⎫
∆S = ∫ ⎪ Q g Qa Q T
T ⎬→ − =0→ a = e ec [1.4]
∆S = 0 ( procesos reversibles )⎪⎭ Tg Te Q g Tg
Combinando la ec [1.2] con ec [1.3] y [1.4], se obtiene:
Te Tg − Te
η = 1− = ec [1.5]
Tg Tg
Tg − Te
W = Qg ⋅ ec [1.6]
Tg
• Máquina frigorífica:
El rendimiento se define como:
Qe
η= ec [1.7]
W
Se aplica el primer principio de la termodinámica:
∆U = Q − W ⎫
Q = Qe − Qc ⎪
⎪
⎬ → W = Qc − Qe ec [1.8]
W <0 ⎪
∆U = 0 (ciclo )⎪⎭

Se aplica el segundo principio de la termodinámica:
δQ ⎫
∆S = ∫ ⎪ Qe Qc Q T
T ⎬→ − =0→ c = c ec [1.9]
∆S = 0 ( procesos reversibles )⎪⎭ Te Tc Qe Te
Combinando la ec [1.7] con ec [1.8] y [1.9], se obtiene:
1 Tc
η= = ec [1.10]
Te Tc − Te
−1
Te
Tc − Te
W = Qe ⋅ ec [1.11]
Te
El trabajo en ambas máquinas debe ser el mismo, se igualan ec [1.6] y ec

[1.11], obteniéndose la expresión del COPideal:
Te Tg − Te
COPideal = ⋅ ec [1.12]
Tg Te − Te

1.2 Componentes de una máquina ideal y su COP.
La máquina de absorción ideal se compone de un generador, un

condensador, un evaporador y un absorbedor. Además de una válvula de
expansión entre el condensador y el evaporador, encargada de reducir la
presión entre estos dos equipos, otra de estrangulamiento con la misma
función que la anterior pero situada entre el generador y el absorbedor y una
bomba que impulsa la disolución desde el absorbedor hasta el generador, en la
Figura 1.2 se representa de forma gráfica los elementos de la misma:
Figura 1.2 Máquina de absorción [1].

Se asumen las siguientes suposiciones:
-. Equilibrio termodinámico en cada una de las trasformaciones.
-. Intercambiadores de calor de área infinita.
-. La temperatura de alimentación de la disolución diluida es la

misma que sale del generador.

Existen dos zonas bien diferenciadas, la de alta presión formada por el

generador y el condensador y la de baja presión constituida por el evaporador y
el absorbedor. Además de existir dos ciclos: el ciclo del refrigerante y el de la
disolución.
El funcionamiento de la máquina de absorción se describe a continuación.

Una fuente de calor proporciona el calor Qg la temperatura tg al generador,
separando el refrigerante de la disolución (punto1). El vapor del refrigerante a
la temperatura tg y presión Pc es conducido hacia al condensador donde
condensa a la temperatura tc y presión Pc (punto2). El refrigerante líquido sufre
una expansión isoentálpica en la válvula de expansión, donde se transforma
parcialmente en vapor, y disminuye su presión hasta Pe (punto 3). El
refrigerante hierve en el evaporador a la temperatura te y presión Pe (punto 4).
El refrigerante vapor, a la temperatura te y presión Pe, es transportado hasta el
absorbedor, donde es absorbido por la disolución a la temperatura ta y presión
Pe. Como el absorbedor y el condensador transfieren calor al medioambiente,
se suele cumplir que ta = tc.
La disolución caliente, a la temperatura tg y presión Pc, sale del generador

(punto 7), pasa a través de la válvula de estrangulamiento y retorna al
absorbedor a la temperatura tg y presión Pe (punto 8), donde se mezcla con la
disolución existente en el absorbedor. Ésta, a la temperatura tc y presión Pe
(punto 5), es bombeada al generador a la presión Pc, donde se calienta hasta la
temperatura tg (punto 6) continuando el ciclo.
En el condensador, el refrigerante se transforma a estado líquido cediendo

su calor latente Qc a un fluido exterior, en cambio en el evaporador el
refrigerante hierve a baja presión absorbiendo el calor Qe del fluido exterior que
refrigera. El conjunto absorbedor- bomba se denomina compresor
termoquímico porque realiza la misma función que el compresor mecánico ya
que el vapor procedente del evaporador entra en el absorbedor como vapor
saturado y sale en el generador como vapor recalentado al Igual que el fluido
de trabajo del ciclo de compresión mecánica.


El COP se define como:
Qe
COP = ec [1.1]
Qg
Aplicando el primer principio de la termodinámica en cada uno de los

componentes de la máquina de absorción, asumiendo las hipótesis:
-. Estado estacionario.
-. Se desprecian la diferencia de energía tanto cinética como

potencial.
Se obtienen las siguientes expresiones:
Q g = mr ⋅ (h1 − h7 ) + md ⋅ (h7 − h6 ) ec [1.13]
−Q c = mr ⋅ (h1 − h2 ) ec [1.14]
Q e = mr ⋅ (h4 − h3 ) ec [1.15]
−Q a = mr ⋅ (h4 − h8 ) + md ⋅ (h8 − h5 ) ec [1.16]
Wb = md ⋅ (h6 − h5 ) ec [1.17]

Suponiendo que la transformación en la válvula de expansión y de

estrangulamiento es isoentálpica y que el trabajo realizado por la bomba (Wb)
es despreciable en comparación con el calor suministrado al generador (h2=h3;
h5=h6; h7=h8), se obtiene:
(h4 − h2 )
COP = ec [1.18]
⎡ ⎛ md ⎞ ⎤
⎢(h1 − h7 ) + ⎜⎜ ⎟⎟ ⋅ (h7 − h5 )⎥
⎣ ⎝ mr ⎠ ⎦
Realizando un balance de masa en el generador, se obtiene el valor del

⎛m ⎞
cociente ⎜⎜ d ⎟⎟ , en función de la concentración de absorbente (X):
⎝ mr ⎠
Figura 1.3. Generador [1].

• Balance de masa en el generador
dm ⎫
= ∑ mei − ∑ m si ⎪⎪
dt
⎬ → mr1 + md 7 = md 6 → md 7 = md 6 − mr1 ec [1.19]
i i
= 0 (estado estacionario )⎪
dm
dt ⎪⎭
• Balance de masa de refrigerante
X ra = concentración de refrigerante en el absorbedor.

X rg = concentración de refrigerante en el generador.
mr 6 ⎫
X ra = → mr 6 = X ra ⋅ md 6 ⎪
md 6 ⎪ md 6 1 − X rg ec [1.20]
⎬ → mr1 + m r 7 = mr 6 → =
m m r1 X ra − X rg
X rg = r 7 → mr 7 = X rg ⋅ md 7 ⎪
md 7 ⎪⎭
Se cumple:
X ra + X a = 1⎫ md 6 1 − X rg m Xg
⎬→ = → d = ec [1.21]
X rg + X g = 1⎭ m r1 X ra − X rg mr Xg − Xa
Si mr1= 1 Kg, la ec [1.21] permite obtener el caudal de disolución

bombeado por cada kilogramo de refrigerante producido en el generador:
Xg
md = ec [1.22]
Xg − Xa

1.3 Fluidos de trabajo.
El sistema termodinámico de refrigeración por absorción emplea una

disolución, generalmente de dos o tres componentes. En el primer caso los
componentes de la disolución se denominan: refrigerante y absorbente.
Un refrigerante es una sustancia química que absorbe calor de un fluido

para que éste baje su temperatura. Las propiedades que debe poseer son:
a) Elevado calor latente de vaporización, a fin de disminuir en lo

máximo posible el caudal másico de refrigerante que ha de
circular por la instalación, con la consecuente reducción del
tamaño de la misma así como de las pérdidas de carga en
tuberías e intercambiadores de calor.
b) Presión de condensación a la temperatura de condensación no
muy elevada, debido a que si es elevada se dificulta la
estanqueidad de la instalación, facilitando las fugas de
refrigerante. Además aumenta el coste del compresor y del
equipo situado en la zona de alta presión (necesita una elevada
resistencia mecánica).
c) Presión de vapor a la temperatura de evaporación no debe ser
excesivamente baja debido a que si es menor que la
atmosférica se presentan problemas de estanqueidad y entrada
de aire disminuyendo la eficiencia del ciclo, compra de equipos
de purga de aire y aparición de la corrosión y formación de
cristales de hielo.
d) Volumen específico del vapor del refrigerante a la presión de
evaporación debe ser reducido, a fin de disminuir el tamaño del
compresor.
e) Temperatura de congelación baja para que no se solidifique en
las condiciones de trabajo normales.
f) Conductividad térmica elevada para disminuir las superficies de
intercambio en la medida de lo posible.


g) Viscosidad reducida a fin de disminuir las pérdidas de carga en

las tuberías e intercambiadores.
h) Comportamiento inerte frente al material que constituye la
instalación.
i) Debe ser barato, para disminuir costes y que el precio total de la
máquina sea lo más bajo posible.
El absorbente es una sustancia que tiene una gran afinidad con el

refrigerante. Esta afinidad permite que el refrigerante en estado vapor sea
absorbido por el absorbente aumentando la concentración de éste en la
disolución. Las propiedades que debe cumplir son:
a) Estabilidad en las temperaturas de trabajo del ciclo.

b) Fácilmente soluble en el refrigerante.
c) Punto de ebullición superior al del refrigerante para permitir su
fácil separación.
d) Temperatura de cristalización del eutéctico, formado con el
refrigerante, baja para la no formación de cristales.
e) Baja agresividad química para evitar la corrosión en la
instalación.
f) Barato y poco tóxico.
Por ser el fluido de trabajo una disolución, las características del ciclo y su
rendimiento, dependerán de sus propiedades físicas, de las cuales se
destacan:
a) Solubilidad total a las temperaturas y presiones del ciclo.

b) Absorción fácil del refrigerante por el absorbente.
c) Tendencia rápida al equilibrio.
d) Fácil separación del refrigerante puro.
e) Baja volatilidad del absorbente.

f) Baja viscosidad.
Las disoluciones más utilizadas comercialmente son el par amoníaco-
agua [NH3 (refrigerante) - H20 (absorbente)] y bromuro de litio-agua [LiBr
(absorbente) – H2O (refrigerante)]. El primero de ellos es empleado en
climatización en sistemas de gran potencia y de pequeña potencia mientras
que el segundo en sistemas de aire acondicionado de mediana y gran potencia.
1.4 Tipos de máquina de absorción.
Principalmente existen dos tipos de máquinas de absorción: simple efecto

y doble efecto.
La configuración de la máquina de simple efecto viene representada en la

Figura 1.4:
Figura 1.4 Máquina de absorción de simple efecto [1].

Gracias al recuperador de calor introducido entre el absorbedor y el

generador, la eficiencia del ciclo aumenta. Esto es debido a que se aprovecha
parte del calor que llega al absorbedor para calentar la disolución que se dirige
al generador, calor que no será necesario que aporte la fuente.

La fuente de calor puede ser de combustible fósil o de calor residual. Ésta

última es de gran interés, debido a que emplea un intercambiador de calor que
transfiere el calor residual a un fluido térmico intermedio y éste es bombeado
hacia el generador de la máquina de absorción para que se dé el proceso de
separación entre el refrigerante y el absorbente. En el presente proyecto se
estudia el calor que es capaz de transmitir el agua calentada a través de
colectores solares térmicos planos de vacío al generador para que la máquina
de simple efecto comience a funcionar.
El COP de la máquina de absorción de simple efecto es:
Qe (h4 − h2 )
COP = → COP = ec [1.23]
Qg + Wb ⎡ ⎛ md ⎞ ⎤
⎢(h1 − h7 ) + ⎜⎜ ⎟⎟ ⋅ (h7 − h6 ) + Wb ⎥
⎣ ⎝ mr ⎠ ⎦
A diferencia con ec [1.18], el efecto del recuperador de calor es apreciable

h6>h5 y se tiene en cuenta el trabajo de la bomba Wb. Esta expresión se
obtiene realizando los balances de energía en el generador, condensador,
evaporador y absorbedor que coinciden con las ecuaciones [1.13-1.17]
Las únicas variables no conocidas son la temperatura de entrada de la

disolución al generador (t6) y la de entrada del absorbedor (t8).Para ello el
recuperador de calor funciona con una eficiencia conocida (ε):
t 7 − t8
ε= → t 8 = t 7 − ε ⋅ (t 7 − t 5 ) ec [1.24]
t7 − t5

Además se cumple:
Q R = (m d − 1) ⋅ Cp g ⋅ (t 7 − t 8 )⎫ QR
⎬ → t6 = t5 + ec [1.25]
Q R = m d ⋅ Cp a ⋅ (t 6 − t 5 ) ⎭ m d ⋅ Cp a
Las máquinas de doble efecto se encuentran constituidas por dos

generadores de refrigerante, el generador de alta temperatura (GA) y el de baja
temperatura (GB), dos recuperadores de calor de la disolución (RA y RB), un
condensador y un subenfriador del líquido, una válvula de expansión , el
evaporador y el absorbedor. El equipo necesario resulta ser el mismo que en
las de simple efecto pero en este caso el generador se desdobla en dos.
La zona de alta presión está constituida por el evaporador y el

absorbedor, mientras que la de baja presión la constituyen el condensador y el
generador. Se pueden encontrar dos configuraciones en serie si se suministra
el mismo caudal de disolución sucesivamente al generador de baja y al de alta.
Si el caudal se reparte entre los dos generadores la distribución será en
paralelo.
Su principal uso está reservado a la recuperación de calor residual y el

efecto útil es el calor producido en el absorbedor, aunque existen casos en los
que no es posible obtener provecho económico simultáneo del condensador y
el absorbedor. Este tipo de máquinas se encuentran en desarrollo experimental
por lo que su presencia en el mercado es escasa.
Valores típicos de funcionamiento de las dos máquinas son:
Tabla 1.1 COP y Rango de temperaturas de funcionamiento [5].

Tipo de máquina COP Rango de temperaturas (ºC)
Simple efecto 0,6-0,7 75-90
Doble efecto 1,2-1,5 120-160

En el presente proyecto se estudiará la máquina de absorción de simple

efecto.
1.5 Consideraciones ambientales de la máquina de absorción.
La máquina de absorción es una alternativa para llegar al desarrollo

sostenible energético debido a:
- Emisión de CO2 baja, ya que pueden funcionar con fuentes de calor

residual o por medio de combustibles fósiles de baja emisión. De
esta manera se reduce el efecto invernadero que existe en la Tierra.
Actualmente se encuentra en fase de investigación el uso de la
energía solar como fuente de calor para la máquina de absorción
(Proyecto SACE) [3].
- Bajo potencial de destrucción del ozono, porque los refrigerantes

empleados son los denominados naturales (agua, amoníaco…), que
al emitirlos a la atmósfera no se combinan con los radicales libres de
oxígeno favoreciendo de esta manera la conservación de la capa de
ozono.
- La existencia de un creciente desarrollo del ciclo de absorción para

la consecución de máquinas capaces de trabajar con igual o mejor
eficiencia que las máquinas de compresión mecánica.
Al emplear un calor residual o renovable como fuente de energía de la

máquina de absorción se requiere evaluar el calor disponible y el flujo
energético de la máquina para que opere sin interrupción. Comparando ambos
flujos de energía se deducirá si la fuente de calor es capaz de accionar la
máquina o será necesario el empleo de otra.
Para evaluar el calor disponible y su fiabilidad es necesario conocer: las

temperaturas y caudales de la fuente, la continuidad del suministro o sus
variaciones temporales, que dependerán del régimen de carga del proceso que
lo genera.


1.6 Descripción de una planta solar [4].
Una planta solar es una instalación cuyo objetivo es emplear la energía

térmica procedente del Sol para fines energéticos (producción de energía
eléctrica, agua caliente para calefacción, ACS…). En este proyecto el uso de
este sistema, se encuentra justificado debido a que el generador de la máquina
de absorción es alimentado con agua caliente proveniente de la planta. A
continuación se representa un esquema de ésta:
Colectores
solares
Figura 1.5 Planta solar acoplada a una máquina de absorción [4].
Los componentes son:
1) Colectores solares: encargados de calentar el agua que circula a

través de ellos por medio de la radiación solar.

2) Sistema acumulador: almacenan el agua caliente proveniente de

los colectores. Generalmente son depósitos de gran capacidad.
Actúa cuando la radiación solar no permite alcanzar la temperatura
deseada para el funcionamiento de la máquina de absorción.
3) Sistema de energía auxiliar: equipo cuya función es suministrar

energía adicional cuando el agua procedente de los colectores no
posee la temperatura adecuada para el accionamiento de la máquina
de absorción. Generalmente es un calentador alimentado por un
combustible fósil o derivados.
Este tipo de instalaciones se han desarrollado en países cuya carga solar

es muy alta (España, Grecia, Italia…).
1.7 Actualidad de la planta solar y de la máquina de absorción.
El sector energético que emplea ambas tecnologías es el del aire

acondicionado. Éste representa un mercado en continuo ascenso debido a la
construcción de edificios comerciales y residenciales, originando una demanda
que exige condiciones de calidad en el ambiente en el que se desarrollan todas
las actividades de los usuarios. Gracias a los sistemas de acondicionamiento
de aire, se regulan tanto la temperatura como la humedad del mismo
cumpliendo los requerimientos de confort exigidos por las personas dentro del
edificio [6].
El empleo comercial de la energía solar para sistemas de aire

acondicionado es relativamente nuevo, concretamente se comienza en 1996
gracias a la investigación de Lamp y Ziegler y continúa en nuestros días
gracias al Proyecto SACE (Solar Air Conditioning in Europe) promovido por la
Unión Europea y al Solar Heating and Cooling Program financiado por la IEA
[3]
.

Destaca que casi todos los sistemas demandan una potencia de 50 kW y

debido a la existencia de diversas fuentes de calor, se pueden encontrar una
gran variedad de tecnologías frigoríficas que pueden ser combinadas con los
colectores solares. Sin embargo, recientes compañías han comenzado a
desarrollar máquinas por debajo de 50 kW promoviendo un desarrollo de
sistemas de baja y media potencia [6].
Este hecho tiene su fundamento en la aparición de colectores solares de

alta eficiencia realizados por las industrias medianas y pequeñas,
proporcionando una buena base para la realización de nuevas tentativas [6].
En la Tabla 1.2 se recoge una representación de máquinas de absorción

de la gama de baja potencia [9]:
Tabla 1.2. Máquinas de absorción de baja potencia. [9]

Fluidos Temperatura
Capacidad
Industria/ I+D País Refrigerante Tipo de de trabajo COP Estado
(kW)
trabajo (ºC)
Pink/Joanneum 1 En
Australia 10 Agua SE NH3/H2O 80 0,55
Research prueba
En
Yazaki Japón 35 Agua SE LiBr/H2O 87 0,71
venta
En
EAW Alemania 15 Agua SE LiBr/H2O 85 0,7
venta
CIimateWell En
Suecia 10 Agua SE+Almacenamiento LiCl/H2O 90 0,7
AB venta
En
Phoenix Alemania 10 Agua SE LiBr/H2O 85 0,7
prueba
En
Rotartica España 4,5 Agua SE LiBr/H2O 90 0,7
venta
Rinnai/Oska 2 Calentador En
Japón 6,7 Agua DE LiBr/H2O 1,2
gas gas venta
3 Calentador En
Robur Italia 15 Aire GAX NH3/H2O 0,9
de gas venta
1. Simple efecto; 2. Doble efecto; 3.Intercambiador Generador-Absorbedor.

Exceptuando las dos últimas máquinas, las demás poseen un razonable

rango de temperaturas de trabajo bajo, entre 80-90ºC, el cual se puede obtener
con un tamaño de colector solar medio [9].
El ciclo de simple efecto domina sobre los demás, debido a que el

proceso de fabricación del aparato de absorción es más fácil y barato en
comparación con los de doble efecto. Los únicos inconvenientes al emplear
este tipo de máquina son: la utilización de una torre de refrigeración y la
relación de inversión inicial (torre de refrigeración, bomba, tuberías etc) y el
coste de operación (consumo de agua y electricidad, mantenimiento etc) [9].
Según el proyecto SACE (2003), el consumo de agua en este tipo de

sistemas en Europa es de 5,3 kg/(kWh)frigorífico. Aunque el precio del agua
puede ser insignificante, es un gasto elevado debido a que con 1,4 kg de agua
refrigerada se obtiene 1(kWh)frigorífico. Lo mismo ocurre con el gasto de eléctrico,
225 W/(kWh)frigorífico, cuya mayoría es empelado en la circulación del agua por
el circuito [9].
En la Figura 1.6, se representa el coste energético (€/kWh) frente al poder

frigorífico de los distintos sistemas de refrigeración existentes [10]:
Figura 1.6 Coste energético frente a capacidad frigorífica [10].

Por lo tanto, se deduce que las máquinas de absorción de simple efecto,

que emplean agua como refrigerante, se deben emplear en aplicaciones con
una carga frigorífica elevada (hospitales, hoteles…) y no en el uso doméstico o
en oficinas de pequeño o medio tamaño [9].
El consumo de energía de los edificios residenciales como de los

comerciales representa un 40% en la demanda energética europea. El
consumo de energía final en 2002 representó el 40,3 % debido a las
condiciones adversas del ambiente urbano donde se desarrolla la gran mayoría
de actividades laborables. Además el número de sistemas de aire
acondicionado se ha incrementado en 5 veces su valor a lo largo de estos
últimos 20 años, así como el espacio dedicado al aire-acondicionado 30
millones de m2 (1980) a 150 millones de m2 en el 2000. También ha aumentado
el consumo de energía de los equipos de aire-acondicionado desde 6 TJ en
1990, 40 TJ en 1996 y se estima 160 TJ en 2010 [3].
Esto ha significado un aumento del precio de la electricidad, así como

colapsos en las líneas de suministro eléctrico (sobre todo en verano) y un alto
[3]
índice de emisión de CO2 (20 veces más que en 1990) . Lo que conlleva a
encontrar nuevas alternativas para la producción de frío reduciendo el consumo
de energía y el coste de adquisición y mantenimiento de equipos.
Hoy en día se encuentran instalados en Europa 70 equipos de aire-

acondicionado que emplean colectores solares. La mayoría se encuentran en
Alemania y España [6].
Figura 1. 7 Sistemas instalados en Europa [6].


El tipo de máquinas empleadas en Europa se representa en la Figura 1.8

[6]
:
Figura 1.8 Tipos de máquinas empleadas en Europa [6].
El poder frigorífico de todos los sistemas es 6,3 MW con un área total de

colectores alrededor de los 17500 m2, destacando que el 59 % de los equipos
utilizan la máquina de absorción [6].
En cuanto a España, la demanda de energía producida en los inmuebles

por los sistemas de aire-acondicionado, ocupa un lugar importante en la
demanda total de energía. Gracias al nuevo Código Técnico la demanda
frigorífica de un edificio se encuentra regulada por lo que se propicia la
instalación de las nuevas innovaciones sobre los sistemas de refrigeración [5].
Debido a la abundante fuente de energía solar que recibe España (1200-

1500 kWh/m2 año) junto con las ordenanzas municipales y el nuevo CTE se
estima un aumento de la superficie de colectores solares de 0,014 m2/hab
(2003) a 0,15-0,30 m2/hab (2010). El uso fundamental de esta energía se
destinará a las máquinas de absorción, DHW (Domestic Hot Water), ACS
(Agua Caliente Sanitaria) y calefacción [5].

Se presentan a continuación algunos ejemplos destacables de

instalaciones tanto de carácter experimental como comercial:
Escuela de Tecnologías Energéticas Renovables de la Universidad

de Naresuan (Tailandia) [2].
Construida en 2005. El sistema consiste en cinco componentes:
1. Colectores solares de tubo: con una superficie total de 72 m2.

Consistente en la unión en serie de 5 unidades con 6 tubos cada una.
Además cada tubo posee 2,4 m2 de área de absorción.
2. Máquina de absorción de simple efecto LiBr/H2O de 10 Tn: modelo
Yazaki WFC SC-10. Sus condiciones de operación son 70-95ºC con un
consumo de agua caliente de 2,4 l/s produciendo 1,5 l/s de agua a 7ºC.
3. Torre de refirgeración: modelo BKC-30 RT con capacidad de 30Tn.
4. Depósitos de agua caliente: de 400 l de capacidad.
5. Recuperador de calor.
6. Unidades de ventilación: cuya cantidad son 4. Modelo Ecoaire WQW
32 V com un consumo total eléctrico de 9,4 KW.
Figura 1.9 Esquema de la instalación [2].

Figura 1.10 Colectores solares [2].
Departamento de Ingeniería Automática y Sistemas de la

[4]
Universidad de Sevilla (España) .
El equipo consiste:
1. Colectores solares: Colectores planos que trabajan entre 60-100ºC y

suministran una potencia nominal de 50 kW. La superficie total es 151
m 2.
2. Depósitos: con una capacidad de 2500 l, trabajando en paralelo.
3. Sistema auxiliar: consistente en un calentador de gas consumiendo
una potencia de 60 kW.
4. Máquina de absorción: cuyo par de trabajo es LiBr/H2O. Su rango de
temperatura de trabajo es 75-100ºC.

Figura 1.11. Croquis de la instalación [4].
Figura 1.12. (a) Colectores solares, (b) Calentador de gas, (c) Depósitos y (d) Máquina de
absorción [4].

Sistema de aire acondicionado solar en una fábrica de cosméticos

situada en Inofita Viotas (Grecia) [11].
Construido en 1999. Se compone de:
1. Dos máquinas de absorción que desarrollan una potencia frigorífica

de 350 kW cada una.
2. Superficie total de colectores planos de 2700 m2.
3. Torres de refrigeración.
4. Tres compresores que desarrollan una potencia de 350 kW.
Figura 1.13 Colectores solares [11].
Figura 1.14 Máquinas de absorción [11].


Sistema de refrigeración solar de un viñedo en Banyuls (Francia)

[11]
.
Operativa desde 1991.Se encuentra equipada con:
1. Dos máquinas de absorción de 52 kW, que acondicionan el agua

empleada en el sistema de ventilación cuyo caudal es de 250000 m3
de aire /h.
2. 130 m2 de tubos colectores.
Figura 1.15 Distribución de los colectores solares [11].
Figura 1.16 Detalle de los colectores solares [11]

2. El COLECTOR SOLAR

2. El Colector Solar.
Un colector solar es una especie de intercambiador de calor que

transforma la energía radiante en calor. La transferencia de energía se hace
desde una fuente radiante (sol), hacia un fluido (agua o aire generalmente) que
circula por los tubos o ductos del colector. El flujo de energía radiante que
finalmente intercepta el colector, proviene básicamente del rango visible del
espectro solar (longitudes de onda entre 300 y 3000 nm) y es por naturaleza
variable con el tiempo [12].
Los únicos problemas que presentan este tipo de dispositivos son: los
bajos y variables flujos de calor (en el mejor de los casos 1100 W/m2) y la
importancia relativa de la radiación. Esto último se justifica con la instalación de
depósitos que acumulan el agua proveniente de los colectores, la cual puede
ser empelada en momentos donde la radiación solar sea insuficiente [12].
2.1 Colectores solares planos.
Los colectores solares planos son diseñados para aplicaciones donde se

requiere que la energía sea liberada a bajas temperaturas, debido a que la
temperatura de este tipo de colectores, difícilmente pasa los 100ºC [13].
Las ventajas que se obtienen de este tipo de colectores con respecto a

los colectores de enfoque, es que éstos utilizan la energía solar directa y
difusa, no requieren movimiento continuo para dar seguimiento al sol,
prácticamente no necesitan mantenimiento y son mecánicamente de
construcción más simple que los colectores concentradores [13].
Las principales aplicaciones de estos dispositivos son en el campo del

calentamiento de agua a nivel doméstico e industrial, acondicionamiento
calorífico de edificios y secado de fruta y granos [13].

2.1.1 Descripción y funcionamiento de un colector plano.
Un colector plano se compone de:
1. Absorbedor o “placa negra”: unido a los tubos o ductos del

colector, se encarga de transformar la energía radiante en calor y
transmitirla al fluido [12]. Debe estar pintado de color oscuro con el
fin de aumentar su capacidad de absorción y reducir la reflexión
de los rayos solares, para que se alcance y mantenga una
temperatura de 100ºC. El pigmento empleado en la coloración de
este elemento debe ser selectivo, es decir, absorber de manera
eficaz la radiación solar e irradiar calor en pequeña proporción,
disminuyendo las pérdidas por radiación. Los absorbedores
poseen un coeficiente de absorción (α) de 0,95 y una emisividad
térmica de 0,1 [18].
Los absorbedores se encuentran formados por una base de metal

brillante y una capa de óxido metálico negro por encima. Los
materiales empleados más frecuentemente son: níquel y óxido de
cromo u óxido de cobre [18].
2. Cubierta transparente: situada en la parte superior de la placa de

absorción negra. Su misión principal es reducir las pérdidas de
[12]
calor por radiación y convección . La cubierta puede ser de
vidrio o plástico, aunque este último tiene la desventaja de
envejecer bajo la acción de los rayos ultravioleta y se hace
necesario reemplazarlo con relativa frecuencia [18].
Los valores típicos de los coeficientes de transmitancia ( τ ),

absorción, reflexión y transmisión son 8, 2, 10 y 88% en una
cubierta de vidrio [18].

Para reducir las pérdidas por radiación desde la superficie del

absorbedor, se realiza un aprovechamiento del “efecto
invernadero”. Este fenómeno consiste en que una capa de
material transparente permite el paso de radiaciones de longitud
de onda corta (por ejemplo: UV, espectro visible) pero no las de
longitud de onda larga (el infrarrojo), produciéndose una reflexión
de éstas últimas hacia el absorbedor incrementando su
temperatura [18].
Para atenuar las pérdidas por convección, se debe realizar una

unión estanca entre la cubierta y la caja, de esta manera se evita
la entrada de aire [18].
[18]
3. Circuito de fluido térmico : son los dispositivos que transportan
el fluido a calentar. Ha de estar en contacto directo con el
absorbedor, de forma que la transferencia de calor al fluido se
realice en las mejores condiciones. Existen tres tipos de
configuraciones:
a) Sistema abierto o por goteo: la placa que constituye el
absorbedor tiene forma ondulada de manera que el
fluido circula por los canales que origina la ondulación.
El fluido penetra en el colector por una tubería de
alimentación (parte superior) que está perforada y es
conducido hasta la tubería de recogida (parte inferior).
Es recomendado para líquidos, aunque presenta
pérdidas por evaporación y ataques de corrosión en el
material de la cubierta.
Figura 2.1 Sistema abierto o por goteo [18].


b) Sistema sándwich: el fluido circula entre dos placas de

forma que la de arriba es la que funciona como
absorbedor al recibir directamente la radiación solar
incidente.
c) Sistema de soldadura: el circuito y el absorbedor están
soldados tomando en general una forma de serpentín o
retícula. En el primer caso pueden originarse pérdidas
de carga del fluido y en el segundo la presencia de
burbujas de aire.
Figura 2.2 Configuraciones del sistema de soldadura [18].
4. Aislante: material localizado en la parte inferior del circuito de

[12]
fluido térmico, impiden la pérdida de calor por conducción .
Debe poseer un bajo coeficiente de conductividad. Ej: lana de
vidrio, poliestireno [18]…
[18]
5. Caja : aloja el conjunto de elementos que forman el colector y
asegura la estanqueidad. Puede construirse de madera, plástico
o chapa metálica.

En la Figura 2.3 se muestra un corte transversal de un colector plano, en

el se identifican los elementos descritos anteriormente:
Figura 2.3 Partes de un colector solar plano [12].
El funcionamiento es sencillo, la energía solar incidente atraviesa una o

varias capas de un material transparente adecuado, antes de alcanzar la placa
de absorción negra al cual está unido el tubo o ducto. En esta placa, es donde
la energía radiante es convertida en calor. Este calor, posteriormente es
transferido por conducción hacia el fluido de trabajo (agua, aire), que es el que
finalmente remueve la energía térmica del colector y la transfiere al tanque de
almacenamiento térmico o al espacio o producto que va a ser calentado, según
la aplicación que se le esté dando[13].

2.1.2 Colector solar plano de vacío.
El colector plano de vacío, posee un diseño muy similar al colector plano

convencional, pero fabricado con materiales que obtienen altos rendimientos, y
que poseen características especiales, tales como superficies selectivas, es
decir, materiales que absorben ciertas frecuencias como los infrarrojos, pero
que apenas los emiten, permitiendo absorciones del orden del 95%. Además,
las pérdidas de convección están minimizadas al máximo mediante el vacío de
la caja, que aumenta el rendimiento considerablemente. Con este sistema
pueden alcanzarse temperaturas de hasta 120 grados centígrados. La
desventaja de este sistema radica en el empleo de materiales más costosos y
el delicado proceso de vacío [17].
Figura 2.4 Colector solar plano de vacío [17].

2.1.3 Ecuaciones energéticas de un colector solar plano [18].
Se describen las principales ecuaciones que definen a un colector solar

plano:
1) Rendimiento del colector solar plano.
Se define el rendimiento de un colector solar plano (ηc) como: relación

existente entre la energía solar extraída por el fluido circulante y la energía
solar incidente sobre su superficie:
Qu
ηc = ec [2.1]
I gl ⋅ A
Donde: Qu= calor útil transmitido al fluido por la unidad de tiempo [W].
Igl = radiación global recibida por la unidad de tiempo y unidad de

área del captador [W/m2].
A = Área del captador [m2].
En un cierto intervalo de tiempo [t1, t2], será:
t2
∫ Q dt
t1
u
ηc = t2
ec [2.2]
A ⋅ ∫ I gl dt
t1
La ec [2.2] justifica porque ηc es un valor no constante, debido a que

depende de las condiciones de trabajo, que pueden influir en el valor de la
energía incidente y en la transferida al fluido. Así pues, en el rendimiento
influirá, por ejemplo, la velocidad del viento (pérdidas por convección), su

temperatura, la intensidad de la radiación solar, etc. El rendimiento de un

colector plano se puede cifrar entre un 30-50%.
Se ha de tener en cuenta las causas que influyen en el aprovechamiento

de la energía incidente, así como la transmisión de este calor al fluido. Por lo
tanto se debe evitar en lo posible: las corrientes de convección sobre el
absorbedor, las pérdidas de radiación y las de conducción a través de la
carcasa. Y se aumentarán los coeficientes de: transmitancia a través de la
cubierta transparente, absorbancia del absorbedor y el de transferencia al
fluido.
El principio básico del diseño de un colector es obtener el máximo calor

útil al menor precio (óptimo económico), aunque no siempre lleve consigo el
máximo rendimiento (óptimo energético).
2) Ecuación fundamental del colector.
La ecuación fundamental del colector se obtiene del balance energético:
A ⋅ ( I D ⋅ τ ⋅ α + I d ⋅ τ ⋅ α ) = Qu + Q p + Q c ec [2.3]
Donde: A = área del captador [m2].
ID = radiación solar directa incidente por unidad de área y tiempo

[W/m2].
Id = radiación solar difusa incidente por unidad de área y tiempo

[W/m2].
τ = transmitancia media de la cubierta transparente.
α = absorbancia media del captador.
Qu = calor útil transmitido al fluido, por unidad de tiempo [W].

Qp = pérdidas de calor por conducción, convección y radiación.

[W].
Qc = calor que se almacena en el colector a mayor temperatura

que el ambiente [W].
El valor de Qu puede deducirse si se determinan Qp y Qc. Los coeficientes

τ y α son conocidos, así como Igl. Evidentemente, Qu también puede
determinarse conociendo el caudal de fluido térmico, sus temperaturas de
entrada y salida en el colector y el calor específico Cp.
Para determinar las pérdidas en el captador, Qp, han de evaluarse las

pérdidas por radiación aplicando la ley de Stefan Boltzmann (E = σ·T4) y las
pérdidas por convección. En lo que se refiere al primer concepto, el captador
emitirá radiación en longitudes de onda prefijadas según la ley de Wien [λmax =
2897,9 (µm·K)/T(K)]. Así, por ejemplo, si se supone que la temperatura del
captador es 100ºC, la emisión corresponderá fundamentalmente a: λmax =
2897,9 (µm·K)/373 (K) → λmax ≈ 8 µm (correspondiente al infrarrojo).
3) Pérdidas a través de la cubierta.
La Figura 2.5 representa un colector plano donde se han considerado dos

cubiertas transparentes a temperaturas te y ti.
Figura 2.5 Colector plano con dos cubiertas [18].

Entre la placa (captador) y la cubierta inferior, a temperaturas tp y ti, se

transfiere una potencia térmica dada por la ecuación [2.4]:
σ ⋅ A ⋅ (t p4 − t i4 )
Ppi = A ⋅ h pi ⋅ (t p − t i ) + ec [2.4]
1 1
+ −1
εp εi
Donde: Ppi = potencia térmica transferida entre la placa y la cubierta

inferior [W].
A = área del captador [m2].
εi = poder emisivo de la cubierta inferior a temperatura ti ≈ 0,9.
εp = poder emisivo de la placa a temperatura tp.
σ = constante de Stefan- Boltzmann.
hpi = coeficiente de transferencia de calor por convección entre la

placa y la cubierta inferior, dado por [W/m2 ºC]:
1,14 ⋅ ∆T 0,310
h pi = ⋅ [l − 0,0018 ⋅ (t m − 10)] ec [2.5]
l 0,07
Siendo: ∆T = diferencia de temperatura entre placa y cubierta.
tm = temperatura media de placa y cubierta [ºC].
A continuación se define el coeficiente de radiación de la placa a la

cubierta inferior (hri) [W/m2 ºC]:
σ ⋅ (t p + t i ) ⋅ (t p2 + t i2 )
hri = ec [2.6]
1 1
+ −1
εp εi

Introduciendo esta expresión en la ec [2.4], la potencia térmica transferida

entre la placa y la cubierta inferior queda:
Ppi = A ⋅ (h pi + hri ) ⋅ (t p − t i ) ec [2.7]
Y definiendo R3 [m2 ºC/W]:, como resistencia térmica al paso del calor

desde la placa a la cubierta inferior:
= (h pi + hri )
1
ec [2.8]
R3
La potencia térmica transferida entre la placa y la cubierta inferior queda:
A ⋅ (t p − t i )
Ppi = ec [2.9]
R3
Siguiendo el mismo método, la potencia perdida entre la cubierta inferior,

que está a la temperatura ti y la cubierta exterior, a la temperatura te, puede
escribirse como:
A ⋅ (t i − t e )
Pie = ec [2.10]
R4
Donde:
σ ⋅ (t i + t e ) ⋅ (t i2 + t e2 )
hre = ec [2.11]
1 1
+ −1
εi εe

Si las cubiertas interior y exterior son del mismo material, los

correspondientes valores de ε serán iguales. Sin embargo puede resultar
interesante desde el punto de vista económico que la cubierta interior sea de
plástico y la exterior de vidrio.
Para evaluar la potencia perdida, desde la cubierta exterior, a temperatura

te, hasta el aire ta, se considera que el intercambio radiante tiene lugar entre las
temperaturas te y tcielo mientras que las pérdidas por convección tienen lugar
entre te y ta. Según esto, las citadas pérdidas serán:
σ ⋅ A ⋅ (t e4 − t cielo
4
)
Pea = A ⋅ hea ⋅ (t e − t a ) + ec [2.12]
1 1
+ −1
ε cielo εe
Si se considera que el poder emisivo del cielo es la unidad, al considerar

el cielo como cuerpo negro, la expresión anterior queda:
(
Pea = A ⋅ hea ⋅ (t e − t a ) + σ ⋅ ε e ⋅ A ⋅ t e4 − t cielo
4
) ec [2.13]
La cual se puede escribir de estas dos formas siguientes:
σ ⋅ ε e ⋅ A ⋅ (t e4 − t cielo
4
) ⋅ (t e − t a )
Pea = A ⋅ hea ⋅ (t e − t a ) + ec [2.14]
(t e − t a )
⎛
Pea = A ⋅ (t e − t a ) ⋅ ⎜⎜ hea +
(
σ ⋅ ε e ⋅ A ⋅ t e4 − t cielo
4
) ⎞⎟ ec [2.15]
(t e − t a ) ⎟
⎝ ⎠

Si se halla hrc:
σ ⋅ ε e ⋅ (t e4 − t cielo
4
)
hrc = ec [2.16]
(t e − t a )
Se obtiene:
A ⋅ (t e − t a )
Pea = A ⋅ (hea + hrc ) ⋅ (t e − t a ) = ec [2.17]
R5
Como en régimen permanente no hay acumulación de calor en las

cubiertas:
A ⋅ (t p − t i ) A ⋅ (t i − t e ) A ⋅ (t e − t a )
P= = = ec [2.18]
R3 R4 R5
Esta expresión se puede formular de la siguiente forma:
A ⋅ (t p − t a )
P= = A ⋅ U p ⋅ (t p − t a ) ec [2.19]
R3 + R 4 + R5
Siendo Up el coeficiente de pérdidas por las cubiertas. Este parámetro,

inverso a la suma de las resistencias térmicas, es del orden de 7 W/m2 ºC
cuando la placa lleva recubrimiento selectivo.
Para calcular R3 y R4, es necesario conocer la temperatura de las

cubiertas y como no son conocidas, la resolución de las ecuaciones anteriores
exige un método iterativo.

El coeficiente de pérdidas Up, depende de los coeficientes de convección

entre placas y cubierta (hpi y hie), y entre cubierta y aire (hea); del poder emisivo
de la placa a la temperatura considerada (εp), y del poder emisivo de las
cubiertas externa e interna (εe y εi). Cuando la temperatura de la placa varía
entre 40ºC y 140ºC, se puede utilizar la expresión propuesta por Hottel, Woertz
y Klein:
−1
⎡ ⎤
⎢ ⎥
⎢ N 1 ⎥ σ ⋅ (t p + t a ) ⋅ (t p2 + t a2 ) ec [2.20]
Up =⎢ + ⎥ +
⎢ 3,44 ⎡ (t p − t a )⎤
0 , 31
hea ⎥ 1 2 ⋅ N + f −1
+ −N
⎢ t ⋅⎢ N + f ⎥ ⎥ ε p + 0,05 ⋅ N ⋅ (1 − ε p ) εe
⎣ p ⎣ ⎦ ⎦
Donde: f = (1- 0,04·hea + 5·10-4· hea2)·(1+ 0,058·N)
N = número de cubiertas.
hea = coeficiente de transferencia de calor por efecto del viento

[W/m2 ºC]. Dado por:
hea = 5,7 + 3,8 ⋅ V (según McAdams); V = velocidad del viento [m / s ]) ec [2.21]
εe = emisividad de las cubiertas para el infrarrojo.
4) Pérdidas a través de la parte posterior y laterales.
Suponiendo que la placa tiene unas dimensiones l1 y l2 y área A, que la

profundidad del colector es l3 y el espesor del aislamiento posterior es la, las
pérdidas por conducción en la parte posterior serán:
⋅ (t p − t a )
k
Ppp = A ⋅ ec [2.22]
la
Siendo: Ppp = pérdidas por conducción [W].


A = área del colector [m2].
k = coeficiente de conductividad del aislamiento [W/m ºC].
Ra = resistencia térmica =la/k [m2 ºC/ W].
Si se supone que a lo largo de todo el lateral del colector existe un

aislante de igual naturaleza y espesor que en la parte posterior, las pérdidas
debidas a conducción en el lateral son:
⋅ (t p − t a )
k
Pla = 2 ⋅ (l1 + l 2 ) ⋅ l 3 ec [2.23]
la
Debiéndose incrementarse este valor en una cantidad constante por

efecto de borde. Este fenómeno queda recogido en la siguiente ecuación:
⋅ (t p − t a )
k
Pla = (l1 + l 2 ) ⋅ (2 ⋅ l 3 + l a ) ec [2.24]
la
Y las pérdidas totales serán:
⎡ (2 ⋅ l 3 + l a ) ⋅ (l1 + l 2 ) ⎤
⋅ (t p − t a ) ⋅ ⎢1 +
k
Ppp + la = A ⋅ ⎥ ec [2.25]
la ⎣ l1 ⋅ l 2 ⎦
En función de las resistencias térmicas:
A ⋅ ⋅(t p − t a )
Ppp + la = = A ⋅ U pp + la ⋅ (t p − t a ) ec [2.26]
R1 + R2

5) Pérdidas totales.
La potencia total perdida será la suma de las ecuaciones [2.19] y [2.26]:
PT = P + Ppp + la = A ⋅ U p ⋅ (t p − t a ) + A ⋅ U pp + la ⋅ (t p − t a ) = A ⋅ (U p + U pp + la ) ⋅ (t p − t a )
PT = A ⋅ U co ⋅ (t p − t a ) ec [2.27]
Donde Uco es designado como coeficiente global de pérdidas:
1 1
U co = + ec [2.28]
R1 + R2 R3 + R4 + R5
6) Calor útil transferido por el colector.
El calor útil transferido por el colector al fluido será:
[
Qu = A ⋅ I gl ⋅τ ⋅ α − U co ⋅ (t p − t a ) ] ec [2.29]
Al objeto de conseguir máximo rendimiento y, en consecuencia, máxima

economía de las instalaciones, se han de seleccionar el diseño y la forma de
trabajo del colector a fin de que se consiga maximizar la absorción de energía
en el captador y minimizar las pérdidas.
El uso de la ec [2.29] para determinar el calor útil, presenta el problema

de desconocer la temperatura de la placa tp, por lo que se realiza es sustituye
por la temperatura del fluido entrante en el colector (teco). Para ello, es preciso
aplicar un factor de corrección que conduce a la siguiente ecuación:

[
Qu = FR ⋅ A ⋅ I gl ⋅τ ⋅ α − U co ⋅ (t eco − t a ) ] ec [2.30]
El coeficiente FR se interpreta como la relación entre el calor transferido al

fluido a su paso por el colector y el calor que se podría transferir si la
temperatura de la placa fuera la misma que la del fluido a la entrada del
colector. Esta situación podría darse si el caudal másico fuera tan grande que
el incremento de temperatura a su paso por el colector se hiciera insignificante.
Este parámetro es afectado sólo por las características del colector y el flujo de
fluido a través del mismo.
7) Rendimiento de colectores en régimen permanente.
La ec [2.30] puede utilizarse para hallar el rendimiento de colector solar

dividiendo Qu por la radiación solar global interceptada:
Qu (t − t )
ηc = = FR ⋅ (τ ⋅ α ) − FR ⋅ U co ⋅ eco a ec [2.31]
I gl ⋅ A I gl
Para un colector dado, operando con flujo constante, los valores de A, FR

y Uco son aproximadamente invariables e independientes de las condiciones de
temperatura e irradiación. El coeficiente FR, puede ser expresado por:
⎛ ⎞
´
F ⋅U co
G ⋅ CP ⎜ G ⋅C ⎟
FR = ⋅ ⎜1 − e p ⎟⎟ ec [2.32]
U co ⎜
⎝ ⎠
Donde: G = caudal másico por unidad de área [kg/sm2]. Siendo el

coeficiente F´:
1
1 U co
F´ = ⋅ ec [2.33]
l´ 1 1
+
U co ⋅ (2 ⋅ a ⋅η al + D ) 2 ⋅ (D + d ) ⋅ hci

Donde: hci = coeficiente de transmisión entre el interior de la placa y el

fluido.
2·a = distancia entre los conductos que transportan el fluido por el

interior del colector.
D = ancho del conducto.
ηal = rendimiento de la aleta.
d = profundidad del conducto.
l´ = D + 2·a.
El rendimiento de un colector es representado a través de una línea recta

dada por la ec [2.31], en un gráfico en el cual las abcisas sean (teco-ta)/Igl y en
función de tm y las ordenadas ηc, dando lugar a las curvas de normalización de
un colector. Las condiciones para obtenerlas son: inclinación y orientación
hacia el Sur del colector y trabajando en mediodía, donde la radiación solar es
máxima. Las curvas de normalización son facilitadas por las firmas comerciales
y son empleadas para comparar colectores.
La intersección de esta recta con el eje de ordenadas es (FR ( τ ⋅ α ))n,

además de representar el rendimiento máximo debido a que la temperatura de
entrada del fluido en el colector (teco), coincide con la temperatura ambiente (ta).
El punto de corte con el eje horizontal representa una situación de baja
radiación solar, en este caso las pérdidas de calor son iguales a la radiación
solar interceptada, no hay calor útil transferido al fluido y la temperatura de éste
a la entrada del colector es igual a la temperatura de la placa. Esta temperatura
es la de equilibrio del colector.

A continuación, se muestra una curva de normalización de un colector

solar plano:
Figura 2.6 Curvas de normalización [18].
Según el gráfico superior, a medida que crece la temperatura de entrada

del fluido en el colector, disminuye el rendimiento del colector, lo que resulta
lógico, ya que para temperaturas altas de fluido y de placa absorbente
aumentan las pérdidas por radiación.
8) Pérdidas por transferencia de calor de la instalación solar.
En el instante de arranque de la instalación Qu = 0 y la radiación solar

incidente se le denomina umbral mínimo de Klein, cuya expresión es:
FR ⋅ U co ⋅ (t eco − t a )
[ ]
Qu = FR ⋅ A ⋅ I gl ⋅τ ⋅ α − U co ⋅ (t eco − t a ) I min,klein =
FR ⋅ (τ ⋅ α )
ec [2.34]

Dado que la curva de normalización se obtiene en condiciones

meteorológicas estándar, al mediodía, no se puede utilizar cuando el colector
trabaja a horas distintas, por lo que este parámetro es corregido dando lugar a
la ec [2.35]:
(t m − t a )h FR ⋅ U co ⋅ (t eco − t a )
I min,klein =
(t m − t a )12h FR ⋅ (τ ⋅ α ) ec [2.35]
Es un coeficiente adimensional igual al cociente entre el salto térmico

instantáneo y el salto térmico máximo al mediodía. La variación del producto
(τ ⋅ α ) suele ser muy pequeña durante el período de funcionamiento de los
sistemas solares, por lo que se puede realizar la siguiente hipótesis
(τ ⋅ α ) = (τ ⋅ α )n .
Las pérdidas de calor de la instalación solar, excluidos los colectores, se

desdoblan en pérdidas de los componentes: tanque de almacenamiento, y
otros equipos, por un lado y pérdidas en las tuberías por el otro. Para reducir
estas pérdidas estos dispositivos deben estar aislados térmicamente. Para
calcular estos parámetros se emplean las siguientes ecuaciones:
Q p , tan que = U tan que ⋅ Atan que ⋅ (t m − t a ) ec [2.36]
Q p ,tubo = U l ⋅ l tubo ⋅ (t m − t a ) ec [2.37]
Q p ,ins = Q p , tan que + Q p ,tubo ec [2.38]

9) Calor absorbido por el colector solar.
Es el calor necesario para elevar la temperatura de los componentes del

colector (captador, caja, aislamiento térmico, aire y cristal) desde tai
(temperatura ambiente inicial) hasta las temperaturas de trabajo:
Qa ,co = Qa , placa + Qa ,cubierta + Qa ,caja + Qa , aislante + Qa , aire + Qa , fluidotérmico ec [2.39]
Se detallan las ecuaciones necesarias para calcular cada término del

segundo miembro de la ec [2.39]:
Qa ,cubierta = mcubierta ⋅ C p ,cubierta ⋅ (t cubierta − t ai )

(t a + t m ) ec [2.40]
t cubierta =
2
Qa ,caja = m caja ⋅ C p ,caja ⋅ (t caja − t ai )

(t m + t cubierta ) ec [2.41]
t caja =
2
Qa , aire = m aire ⋅ C p , aire ⋅ (t aire − t ai ) ec [2.42]
Qa , fluidotérmico = m fluidotérmico ⋅ C p , fluidotérmico ⋅ (t m − t ai ) ec [2.43]
10) Calor absorbido por los componentes de la instalación solar.
Es el calor necesario para mantener a la temperatura tm de cada uno de

los componentes (tubería, fluido, bombas, válvulas, etc.) en todo instante:
Qa ,ins = Qa , material + Qa , aislante + Qa , fluidotérmico ec [2.44]

Se detallan las ecuaciones necesarias para calcular cada término del

segundo miembro de la ec [2.44]:
Qa , material = m material ⋅ C p , material ⋅ (t m − t ai ) ec [2.45]
Qa , aislante = m aislante ⋅ C p , aislante ⋅ (t aislante − t ai )

ec [2.46]
Qa , fluidotérmico = m fluidotérmico ⋅ C p , fluidotérmico ⋅ (t m − t ai ) ec [2.47]
11) Energía solar disponible.
Un balance energético horario sobre el colector, permite escribir:
I gl = Qu + I min + Pv ec [2.48]
La radiación solar interceptada se transforma en calor: una parte, Imin, se

invierte en incrementar la temperatura del colector y del sistema hasta alcanzar
la de proceso y compensar las pérdidas al ambiente. Otra serán las pérdidas
originadas por el viento, Pv, que no se pueden calcular ´a priori´ y su
determinación se debe hacer experimentalmente. Por último Qu es el calor útil
generado para el proceso.
Calculadas las pérdidas hora a hora y evaluada la radiación solar, la

energía solar disponible en cada hora será igual a la diferencia entre la
radiación global interceptada y las pérdidas. Si se denomina Imin a las pérdidas,
la energía solar disponible será:
Qu + Pv = I gl − I min
ec [2.49]
I min = I min, Klein + Q p ,ins + Qa ,co + Qa ,ins

La energía solar disponible máxima será:
Qdis = Qu + Pv = I gl − I min ec [2.50]
El rendimiento horario máximo del colector, hasta donde se pueda

calcular será:
Qdis
η h ,c = ec [2.50]
I gl
La integración de (Qu+Pv) a lo largo del día, cuantifica la energía

disponible diaria para el proceso, así como su rendimiento máximo diario:
Qu , max = ∑ Q hco , max ec [2.51]
Qu , max
η d ,c , max = ec [2.52]
∑I gl
Los rendimientos reales sólo se pueden obtener midiendo

experimentalmente el calor útil generado y la radiación global interceptada.

En la Figura 2.7, se representa la radiación solar global sobre un colector

plano y sobre ella la distribución horario Imin. La zona rayada es la energía solar
disponible para el proceso, incluidas las pérdidas originadas por el viento:
Figura 2.7 Energía solar disponible [18].
2.2 Colectores concentradores [13].
Existen muchas aplicaciones, sobre todo a nivel industrial, donde se

necesita que la energía sea liberada a altas temperaturas. Esto no se puede
lograr con los colectores solares planos debido a las características propias de
este tipo de colectores y a que la radiación solar es una energía de baja
intensidad.
En consecuencia, para obtener temperaturas altas (superiores a 100 °C),

es necesario incrementar la intensidad de la energía solar. Para ello, se
disminuye el área por donde ocurren las pérdidas de calor y se interpone un
dispositivo óptico entre la fuente de radiación (sol) y la superficie absorbedora,
que debe ser pequeña comparada con la del dispositivo óptico, dando lugar a
los colectores concentradores.

Con los colectores concentradores de energía solar, se pueden obtener

temperaturas:
1) 100 - 500°C si se usan colectores focales rudimentarios.

2) 500 - 1500 °C si el sistema óptico de los colectores tiene un buen
acabado.
3) 1500 - 3500 °C si el sistema óptico tiene un acabado perfecto.
Aunque con este tipo de colectores se pueden obtener altas temperaturas

de operación, estos presentan varios problemas técnicos desde el punto de
vista ingenieril comparados con los colectores solares planos:
a) Deben orientarse continuamente al sol de manera precisa mediante un

mecanismo apropiado debido a que este tipo de colectores utilizan
únicamente la energía solar directa.
b) El acabado de las superficies que constituyen el sistema óptico no sólo
debe ser de buena calidad, sino que debe mantener sus propiedades
por largos períodos de tiempo sin ser deterioradas por el polvo, lluvia y
medio ambiente, donde generalmente existen componentes oxidantes
y corrosivos.
c) La calidad de los materiales utilizados en el receptor (aislante térmico,
fluido de trabajo, tubos absorbedores y cubiertas) son mayores en este
tipo de colectores, debido a que es ahí donde se obtienen las altas
temperaturas.
El uso de este tipo de colectores no es muy generalizado debido a los

problemas de operación mencionados anteriormente y los costes elevados de
los materiales utilizados en este tipo de colectores focales. A causa de estos
hechos se trabaja con los colectores concentradores fijos, los cuales no poseen
las desventajas de los de enfoque, aunque sólo permiten incrementos
moderados de la intensidad de la radiación solar.

2.2.1 Descripción y funcionamiento de un colector concentrador.
Un colector concentrador se compone de:
1. Receptor: elemento encargado de absorber la radiación y

convertirla en otro tipo de energía (térmica, eléctrica…). En él se
encuentra tanto el absorbedor como el aislante térmico.
2. Concentrador o sistema óptico: dispositivo del colector que dirige
la radiación hacia el fluido a calentar.
Todas las partes mencionadas se encuentran dentro de una caja que
sirve como parte estructural del colector y que puede ser realizada de diversos
materiales como lámina metálica, madera o plástico.
El funcionamiento es igual al de un colector plano, la diferencia es que

con este tipo de colectores se consigue una mayor temperatura que con los
planos (por encima de 100ºC).
2.2.2 Tipos de colectores concentradores.
En la actualidad, se emplean dos tipos de colectores solares

concentradores: los tubos evacuadores y los parabólicos.
2.2.2.1 Tubos evacuadores [14].
Se encuentran formados por dos tubos concéntricos, entre los cuales se

ha producido el vacío a través de la extracción del aire localizado entre ellos.
En el interior se encuentra el absorbedor.
Figura 2.8 Esquema de tubos evacuadores [14].

Algunos colectores emplean un sistema denominado CPC (Colector

Parabólico Concéntrico) para aprovechar la radiación solar que incide entre dos
tubos. Este sistema consiste en una serie de reflectores que dirigen la luz que
cae entre tubo y tubo hacia la parte trasera de los mismos donde es también
aprovechada. Con ello los colectores reciben luz tanto de la parte delantera
como de la trasera. Con el sistema CPC se amplía la superficie efectiva de
captación por metro cuadrado para la tecnología de tubo de vacío factor que
sin embargo siempre estará por debajo de los colectores de placa plana.
Figura 2.9 Esquema de captación de un sistema CPC [14].
Existen dos clases de tubos evacuadores:
[14]
a) Tubos evacuadores simples : este sistema es únicamente utilizado en
calentadores solares termosifónicos. Son tubos evacuados ensamblados
directamente con el depósito acumulador y que por lo tanto contienen agua.
En la pared interior del tubo evacuado se sitúa una capa de color oscuro
de material absorbente. Cuando la radiación solar incide sobre la capa de
material absorbente, se transforma en calor y eleva la temperatura del agua
que está en contacto con él. El agua calentada se eleva por convección y
comienza a ascender siendo reemplazada por agua fría que a su vez se
calienta y reinicia el proceso.
Este tipo de tubo de vacío ofrece la ventaja de tener escasas pérdidas de

calor, y los inconvenientes de ser muy sensible a la presión y de no ofrecer
ninguna protección contra las bajas temperaturas (dilatación del agua al

congelarse y rotura de las tuberías de la instalación), Por lo tanto no es posible

su utilización en zonas con inviernos fríos, sin la inclusión de un calentador
eléctrico que caliente el agua del depósito cuando esta alcanza temperaturas
muy bajas.
Figura 2.10 Tubo evacuador simple [14].
[14]
b) Tubo de vacío de flujo directo : esta tecnología se emplea tanto para
colectores exentos como para calentadores solares compactos con depósito
integrado. Los colectores de tubo de vacío de flujo directo colocan en el interior
del tubo evacuado una plancha de material absorbente adecuado, que hace las
veces de absorbedor transformando la radiación solar en calor. El absorbedor
es recorrido en su superficie por un tubo con flujo directo en el que circula un
fluido que eleva su temperatura en contacto con él.
Los colectores de tubo de vacío de flujo directo, tienen la ventaja de poder

adoptar una posición tanto horizontal como vertical sin mermar su rendimiento,
ya que el tubo puede rotar sobre su eje, inclinándose el absorbedor de la
manera más adecuada. Además de ser utilizable en áreas frías ya que permite

usar las estrategias contra la congelación de uso general en la energía solar

térmica.
Figura 2.11 Tubo de vacío de flujo directo [14].

[14]
c) Tubo de vacío de Heat Pipe : esta tecnología de colectores solares
emplea un mecanismo denominado Heat Pipe (tubo de calor). Consiste en un
tubo cerrado en el cual se introduce un fluido de propiedades específicas.
Cuando el Sol incide sobre el absorbedor adosado al tubo, el fluido se evapora
y absorbe calor (calor latente). Como gas asciende sobre el líquido hasta lo alto
del tubo donde se sitúa el foco frio. Allí se licua (condensa) y cede su calor
latente al fluido que nos interesa calentar volviendo a caer al fondo del tubo por
gravedad. Este proceso se repite mientras dure la radiación del Sol o hasta que
el colector ha alcanzado una temperatura muy alta (de en torno los 130 grados
o más). El Heat Pipe o tubo de calor es considerado como un superconductor
térmico por lo eficaz de su funcionamiento.

Figura 2.12 Esquema de un tubo Heat Pipe [14].
Este tipo de colectores tienen la ventaja de no sufrir pérdidas por la noche

ya que el proceso de transferencia de calor no es reversible (es decir el fluido
caliente o el calor no puede pasar del acumulador al tubo y por lo tanto
perderse). Además cada tubo es independiente por lo que se puede cambiar en
pleno funcionamiento del sistema. Es altamente resistente a las heladas.
Dado que también pueden girar sobre su eje los tubos, existe la
posibilidad de que adopten posiciones verticales y horizontales al igual que
ocurre en los sistemas de flujo directo, aunque en este caso habrá que respetar
una inclinación mínima del largo del tubo para permitir que el fluido, una vez
licuado, descienda por gravedad. También se aplica el sistema CPC.

El funcionamiento de estos colectores, se detalla a continuación y se

muestra en la Figura 2.13:
1) La radiación solar incide en el absorbedor que se calienta y

transmite ese calor al tubo.
2) El calor recibido provoca que el fluido en el interior del tubo se

evapore y ascienda portando energía (calor latente)
3) El fluido evaporado cede su calor latente al fluido más frio que

circula por el exterior de la cabeza del tubo y al hacerlo se licua
4) El fluido de nuevo en estado líquido cae por gravedad al fondo del

tubo para reiniciar el proceso.
Figura 2.13 Funcionamiento de un tubo Heat Pipe [14].

2.2.2.2 Colectores parabólicos.
Este tipo de colectores poseen el mismo funcionamiento que el de los

colectores planos, la diferencia es que la sección de la placa solar es una
parábola, que dirige la radiación solar hacia un tubo o torre por donde circula el
fluido a calentar.
Existen dos configuraciones de colectores solares parabólicos:
a) Colectores cilíndricos parabólicos (CCP) [15]: se compone, básicamente, por

un espejo cilindro parabólico que refleja la radiación solar directa,
concentrándola sobre un tubo receptor colocado en la línea focal de la
parábola. Esta radiación concentrada provoca que el fluido que circula por el
interior del tubo se caliente, transformando así la radiación solar en energía
térmica en forma de calor sensible del fluido. La figura adjunta muestra una
fotografía de un CCP e ilustra su modo de funcionamiento:
Figura 2.14 Colector cilíndrico parabólico [15].
Los colectores cilindro parabólicos (CCP) resultan idóneos para trabajar

dentro del rango de temperaturas 125ºC-400ºC. Gracias a la concentración de
la radiación solar directa que incide sobre el plano de apertura del captador, se

consigue de forma eficiente elevar la temperatura del fluido de trabajo hasta

valores del orden de los 425ºC, pudiendo alimentar procesos industriales
dentro del rango de la media temperatura, convirtiéndolo en un captador ideal
para acoplarlo a una gran diversidad de procesos industriales.
[16]
b) Colectores de Fresnel : las superficies pueden ser planas o curvas,
estando diseñada cada una para reflejar individualmente la radiación sobre el
receptor; su gran ventaja radica en su facilidad de construcción, siendo muy
interesante su utilización en colectores solares de gran potencia y alta
temperatura. El refractor Fresnel utiliza un sistema de concentración parecido a
un montaje de lentes con un punto focal único.
Figura 2.15 Reflector y refractor Fresnel [16].

3. DATOS EXPERIMENTALES

3. Datos experimentales.
3.1 Descripción de la instalación.
La instalación térmica solar se encuentra situada en la azotea del edificio

Betancourt en la Escuela Politécnica Superior de la Universidad Carlos III de
Madrid, localizada en Leganés (España).
Se trata de un sistema de conversión térmica, el cual provee al edificio de:
- Agua Caliente Sanitaria (ACS).

- Calefacción.
- Aire acondicionado a través de una máquina de absorción LiBr-
H2O, empleando colectores solares planos y tubos evacuadores.
La instalación se compone de:
- 50 m2 de colectores solares planos.

- 50 m2 de colectores de tubos evacuadores.
- 20 m2 de colectores solares planos de vacío.
- Intercambiador de calor de placas.
- Dos tanques de almacenamiento o depósitos de agua.
- Bombas.
- Máquina de absorción LiBr- H2O [35 kW].
- Torre de refrigeración.
- Fan- coils.
- Sistema de control.
- Fluidos caloportadores.
- Estación meteorológica.

En la Figura 3.1 se presenta un esquema de la instalación:
Figura 3.1. Esquema de la instalación [38].
Esta instalación produce frío para la obtención de aire acondicionado

realizando el siguiente proceso, el calor que alimenta a la máquina de
absorción es captado por los colectores solares a una temperatura entre 65-
85ºC con una gran eficiencia. En verano, al mediodía, cuando la demanda de
aire acondicionado es elevada, el calor producido por los colectores es de 30
kW, que permite obtener una potencia frigorífica de 15 kW y un COP de 0,5. En
este proyecto se valoriza la producción de frío a través de colectores solares
planos de vacío.

3.1.1 Colectores solares planos de vacío.
Los colectores solares planos de vacío que se utilizan en la instalación

son de la marca THERMOSOLAR modelo 400V. La instalación posee 10
paneles de 2 m2 y una superficie total 20 m2.
Figura 3.2. Colector solar plano de vacío 400Vplus [22], en el tejado [40].
Las características de estos paneles se detallan a continuación [22]:
- Carcasa y marco en aluminio anticorrosivo con aleación de

magnesio.
- Vidrio de seguridad transparente con tratamiento antirreflejos.
- Creación de vacio in situ, consiguiendo una barrera aislante que

evita las pérdidas térmicas de la superficie selectiva absorbedora;
mayor durabilidad del panel ya que se evita la corrosión interna y la
condensación.
- Krypton insade; Inertización del aire restante mediante el llenado

con gas criptón mejorando aun más la conductividad térmica.
- Sistema Tichelmann de conductos y canalizaciones internas

permitiendo la conexión de hasta doce paneles consecutivos.

Las especificaciones técnicas y propiedades térmicas del panel se

recogen en la Tabla 3.1 y 3.2:
Tabla 3.1. Especificaciones técnicas del panel 400V plus.

Dimensiones (largo x alto x espesor) [mm] 2008 x 1008 x 75
Superficie total [m2] 2,03
Superficie útil [m2] 1,76
Peso total [kg] 48
Peso total con cubierta [kg] 60
Producción [kWh/a] 1021
Volumen total de fluido por panel [l] 1,3
Caudal recomendado [l/h] 100
Temperatura de estancamiento [ºC] 30
Radiación máxima a recibir [W/m2] 1000
Tabla 3.2. Propiedades térmicas de los componentes del panel.
Espesor ρ Cp Κ
Componente
(mm) (kg/m3) (J/kg·ºC) (W/m·K)
Placa absorbente AlOx [24] 0,3 3240 620 20
Cubierta de Cristal 4 2225 835 1,4
Caja de Aluminio 1 2770 875 177
Gas Criptón [25] ------ 3,708 248 0,01

3.1.2 Intercambiador de calor [19].
La transferencia de calor del circuito primario al secundario se realiza a

través del intercambiador de placas Alfa Laval M6MF. Las características
técnicas de este equipo son:
Tabla 3.3. Características de las placas.

Superficie máxima de
Material Número
transferencia [m2]
Placas Acero Inoxidable AISI 316 22 38
Tabla 3.4. Dimensiones del intercambiador de calor.

Modelo H W h
M6MF 920 320 140
Figura 3.3. Medidas del intercambiador de calor [19].
Figura 3.4. Intercambiador de calor Alfa Laval M6MF[19].

3.1.3 Máquina de absorción LiBr-H2O.
La máquina de absorción es una máquina de simple efecto de la marca

YAZAKI modelo WFC-10. Los fluidos de trabajo son LiBr (absorbente) y H2O
(refrigerante). Se encuentra situada en la nave de Ingeniería Térmica,
conectada con el agua caliente procedente de los tanques (generador), con la
torre de refrigeración y con los fan- coils.
Los principios termodinámicos que rigen este aparato se encuentran en el

capítulo 1 (pág 9-20).
Se compone de un generador, que se alimenta con agua caliente entre

65-85ºC, un evaporador que enfría el agua procedente de los fan-coils, un
condensador y un absorbedor conectados a una torre de refrigeración. En la
Figura 3.5 se representa un esquema de este equipo:
Figura 3.5. Esquema de la máquina de absorción.

Los tubos del generador son de acero inoxidable y los demás son de
cobre. El movimiento de la disolución se encuentra asegurado gracias a los
movimientos de convección natural de la misma, por lo que no es necesario el
uso de un motor para asegurar la circulación de la solución en el interior del
aparato.
En la Figura 3.6 se observa una fotografía real de la máquina de

absorción:
Figura 3.3. Máquina de absorción [39].

Las especificaciones de la máquina de absorción YAZAKI WFC-10, se

resumen en la siguiente tabla:
Tabla 3.5. Especificaciones de la máquina de absorción.

Modelo WFC-10
Capacidad Frigorífica [kW] 34,9

Producción de frío
Temperatura del agua enfriada
7 (salida), 12,5 (entrada)
[ºC]
Caudal [l/s] 1,67

Agua enfriada
Pérdida de carga [kPa] 37,27
Caudal [l/s] 4,05

Agua de refrigeración
Caudal [l/s] 2,38

Agua del generador
Electricidad Consumo [W] 210
Nivel de sonido Presión del sonido [dB(A)] 49
Control ON/OFF

3.1.4 Tanques de almacenamiento.
Los tanques de almacenamiento son dos y se encuentran situados en la

sala de calderas del edificio Betancourt:
Figura 3.7. Tanques de almacenamiento.
Las dimensiones y propiedades necesarias de los materiales que componen

los depósitos se recogen en las Tablas 3.6 y 3.7:
Tabla 3.6. Dimensiones de los tanques de almacenamiento.
Diámetro Interno Diámetro Altura Interna Altura Externa

Capacidad [l]
[m] Externo [m] [m] [m]
1,2 m 1,284 m 1,55 m 1,634 m 2000

Tabla 3.7. Propiedades de los materiales que conforman los tanques.
Espesor ρ Cp Κ
Material
[mm] [kg/m3] [J/kg·C] [W/m·K]
Acero
5 7832 434 63,9
galvanizado
Poliuretano 35 50 800 0,035
Estos depósitos se emplean para el funcionamiento del aire

acondicionado en el edificio. Poseen una configuración en serie.
3.1.5 Circuitos hidráulicos de la instalación.
A) Circuito primario.
El circuito primario cierra la conexión entre los colectores y el

intercambiador de calor. Está formado por tuberías de cobre sanitario
recubiertas por SH/Armaflex (aislante).
Las dimensiones de las tuberías que componen el circuito primario y las

propiedades sus elementos se detallan en la Tabla 3.8 y 3.9:
Tabla 3.8. Dimensiones de las tuberías del circuito primario.

Tramo Posición tubería Longitud [m] Características
Horizontal 5,44
Colectores-
Intercambiador de Vertical 21,05 De1 = 26 mm
calor
Azotea 2,56 Di1 = 24 mm
Horizontal 5,18 et1 = 1mm

Intercambiador de
Vertical 20,45 eais1 =20 mm
calor - Colectores
Azotea 0,85

Siendo: De1= diámetro exterior de la tubería del circuito primario.
Di1= diámetro interior de la tubería del circuito primario= De1-2et1.
et1= espesor de la tubería del circuito primario.
eais1 = espesor del aislante del circuito primario.
Tabla 3.9. Propiedades de los materiales del circuito primario.

ρ Cp Κ
Componente
[kg/m3] [J/kg·ºC] [W/m·K]
Cobre Sanitario 8960 380 401
SH/Armaflex 32 835 0,037
B) Circuito secundario.
El circuito secundario cierra la conexión entre el intercambiador de

calor y los tanques de almacenamiento. Está formado por tuberías de cobre
sanitario recubiertas por SH/Armaflex (aislante).
Las dimensiones de las tuberías que componen el circuito secundario y

las propiedades sus elementos se detallan en la Tabla 3.10 y 3.11:
Tabla 3.10. Dimensiones del circuito secundario.

Horizontal 2,51 De2 = 26 mm
Intercambiador de
calor- Tanques
Vertical 0,66 Di2 = 24 mm
Tanques - Horizontal 1,85
et2 = 1mm
Intercambiador de
calor Vertical 1,67 eais2 =20 mm

Tabla 3.11. Propiedades de los materiales del circuito secundario.

ρ Cp Κ
Componente
C) Circuito terciario.
El circuito terciario cierra la conexión entre los tanques y el generador

de la máquina de absorción. Está formado por tuberías de cobre sanitario
recubiertas por SH/Armaflex (aislante).
Las dimensiones de las tuberías que componen el circuito terciario y las

Tabla 3.12. Dimensiones del circuito terciario.


Tanques –
Generador

Generador-
Tanques eais3 =20 mm
Vertical 13,69

Tabla 3.13. Propiedades de los materiales del circuito terciario.

ρ Cp Κ
Componente
D) Circuito cuaternario.
El circuito cuaternario cierra la conexión entre el evaporador de la

máquina de absorción y los fan-coils. Está formado por tuberías de cobre
sanitario recubiertas por SH/Armaflex (aislante).
Las dimensiones de las tuberías que componen el circuito terciario y las

Tabla 3.14. Dimensiones del circuito cuaternario.


Evaporador – Fan
coil

Fan coils -
Evaporador eais4 =20 mm
Vertical 11,17

Tabla 3.15. Propiedades de los materiales del circuito cuaternario.

ρ Cp Κ
Componente
3.1.6 Fan-Coils.
En esta instalación se utiliza un Fan-Coil de la marca Jucker para
transferir el frío del circuito cuaternario al habitáculo a acondicionar, que es la
nave del departamento de Ingeniería Térmica.
3.1.7 Fluidos de trabajo.

Se detallan los fluidos que intervienen en la transferencia de calor en la
instalación:
-Circuito primario: disolución de etilenglicol y agua en una

composición de 20%/80%.
-Circuitos secundario, terciario y cuaternario: agua.

En la Tabla 3.16 se resumen las propiedades de los fluidos anteriores:
Tabla 3.16. Propiedades de los fluidos de la instalación.

ρ Cp Κ µ
Fluido
(kg/m3) (J/kg·ºC) (W/m·K) (Pa·s)
Etilenglicol[20] 1109 2400 0,26 20·10-3
Agua 1000 4191 668·10-3 389·10-6
3.1.8 Sistema de regulación, bombeo y adquisición de datos.
A) Sistema de bombeo.
En toda la instalación se pueden encontrar cuatro bombas, todas ellas de

la marca Grundfos, a continuación se detalla su ubicación y modelo. (Tabla
3.17).
Tabla 3.17. Bombas de la instalación.

Modelo UPS 25-60 UPS 25-60 UPS 32-80 UPS 32-120/F
Ubicación Terciario Terciario Secundario Primario
B) Sistema de regulación.
A lo largo del trazado de la instalación se pueden encontrar un total de

22 válvulas todo-nada, y un total de cuatro válvulas de regulación (de pantalla),
una en la salida de cada batería de colectores planos.

C) Equipos de adquisición de datos.
La instalación solar cuenta con una pequeña instalación meteorológica

que registra cada 2 segundos los siguientes parámetros: dirección del viento,
velocidad del viento, temperatura del aire, radiación global, radiación difusa,
humedad relativa, presión y lluvia. (Esta información se puede consultar vía
Internet en la página http://pc-134-222.uc3m.es/meteo.htm).
Figura 3.7. Vista de la aplicación que muestra los valores registrados en la instalación
meteorológica, actualizados cada minuto [21].
Existen numerosos puntos a lo largo de la instalación en los que se toman

medidas de temperatura cada dos segundos y hacen la media cada diez
minutos de forma automatizada. Los instrumentos empleados para la recogida
de este dato son termopares que se encuentran situados en la salida de los
colectores, las tuberías, los tanques de almacenamiento y en la máquina de

absorción. Se registran en un ordenador que coordina la operación de recogida

de todos los parámetros antes citados.
El funcionamiento de la instalación se resume en la Tabla 3.18:
Tabla 3.18. Funcionamiento de la instalación

Condiciones de Circuito de Circuito de Circuito de
funcionamiento colectores almacenamiento consumo
Radiación > 50 W/m2 SI
Radiación ≤ 50 W/m2 NO
Temperatura de salida de
los colectores >
SI
Temperatura del tanque +
5°C
Temperatura de salida de
los colectores ≤
NO
Temperatura del tanque +
5°C
Hora de funcionamiento
SI
del aire acondicionado
Otra hora NO
Temperatura externa ≤
NO
25°C
Temperatura externa >

SI
25°C
Temperatura del tanque >
Temperatura de salida del

SI
generador
+ 10°C
Temperatura del tanque≤
Temperatura de salida del NO

generador +10°

3.2 Instalación térmica solar.
La instalación solar se encuentra situada en la azotea del edificio

Betancourt de la Universidad Carlos III de Madrid en la población de Leganés.
En la Tabla 3.19 se recogen los datos experimentales de los parámetros

de intensidad global y la temperatura a la entrada y salida del agua que circula
por los colectores. Se representan tres días significativos del intervalo 30/06-
11/07 del año 2007, cuyo período de funcionamiento es 7:00-22:00. Esto último
tiene como objetivo comparar estos resultados experimentales con los de la
simulación en los mismos días, para proceder a una validación parcial del
modelo.

Tabla 3.19. Radiación global y temperatura de entrada y salida del agua que circula por los colectores.
Radiación global [W/m2] Temperatura entrada del colector [ºC] Temperatura salida del colector [ºC]
Tiempo [min]
30/06/2007 05/07/2007 11/07/2007 30/06/2007 05/07/2007 11/07/2007 30/06/2007 05/07/2007 11/07/2007
7:00 6,538 5,458 2,622 24,324 22,196 23,786 17,720 15,577 17,475
7:10 17,548 15,829 9,977 23,003 20,907 22,535 18,520 16,537 18,327
7:20 31,134 31,579 22,538 23,005 20,867 22,428 19,084 16,984 18,734
7:30 49,632 50,819 38,749 23,107 20,944 22,471 19,562 17,418 19,153
7:40 69,983 70,652 56,896 23,370 21,144 22,680 20,020 17,937 19,674
7:50 93,880 93,385 76,950 23,621 21,413 22,973 20,482 18,479 20,198
8:00 119,622 119,462 101,134 23,889 21,691 23,227 21,053 18,953 20,673
8:10 147,863 147,993 129,161 24,268 22,013 23,542 21,642 19,486 21,157
8:20 174,823 178,324 157,298 24,751 22,413 23,985 22,430 20,163 21,843
8:30 205,107 211,411 187,579 25,393 22,869 24,536 23,399 20,869 22,655
8:40 237,776 244,635 219,759 26,209 23,598 25,273 24,556 21,859 23,691
8:50 270,381 278,217 251,963 27,237 24,499 26,213 26,044 23,181 24,953
9:00 303,077 312,716 285,559 28,567 25,596 27,429 27,973 24,811 26,642
9:10 335,967 346,836 319,542 30,165 26,992 28,908 30,251 26,786 28,667
9:20 369,318 381,408 352,833 32,062 28,778 30,682 32,802 29,237 31,057
9:30 402,074 414,699 385,341 34,209 30,791 32,714 35,493 31,853 33,631
9:40 435,642 448,833 419,669 36,632 33,042 35,040 38,363 34,490 36,425
9:50 466,786 481,559 453,130 39,289 35,518 37,611 41,323 37,259 39,373
10:00 498,843 513,314 482,639 42,139 38,258 40,356 44,456 40,304 42,327
10:10 529,100 542,298 515,067 45,177 41,192 43,237 47,766 43,492 45,291
10:20 559,120 572,769 548,542 48,192 44,163 46,167 51,253 46,475 48,370

Continuación de Tabla 3.19.Radiación global y temperatura de entrada y salida del agua que circula por los colectores.
Tiempo [min] Radiación global [W/m2] Temperatura entrada del colector [ºC] Temperatura salida del colector [ºC]
30/06/2007 05/07/2007 11/07/2007 30/06/2007 05/07/2007 11/07/2007 30/06/2007 05/07/2007 11/07/2007
10:30 580,722 601,415 570,796 51,207 46,895 48,968 54,385 49,310 51,274
10:40 610,475 629,679 603,960 53,900 49,551 51,686 57,621 52,190 54,180
10:50 640,569 660,151 630,462 52,896 52,071 53,151 59,991 54,819 56,953
11:00 664,070 685,916 655,244 51,400 52,335 51,556 60,873 57,130 58,959
11:10 690,572 713,167 670,361 51,306 50,283 50,667 61,796 58,348 59,706
11:20 715,779 736,997 692,151 51,584 49,930 50,705 63,327 59,440 60,585
11:30 736,398 759,742 718,375 52,109 50,222 51,020 64,789 60,872 61,803
11:40 755,719 780,943 753,057 52,653 50,713 51,691 66,091 62,335 63,563
11:50 773,017 803,027 762,890 53,257 51,148 52,138 67,449 63,519 64,532
12:00 805,465 825,104 783,375 53,822 51,771 52,856 69,014 64,864 66,400
12:10 826,716 842,866 806,405 54,414 52,264 53,422 70,526 66,201 67,688
12:20 854,492 862,759 822,806 54,964 52,705 53,959 71,992 67,413 68,805
12:30 870,482 880,110 830,679 55,532 53,190 54,464 73,208 68,560 69,667
12:40 890,060 899,615 843,231 56,018 53,751 54,992 74,271 69,715 70,517
12:50 905,682 913,438 867,900 56,447 54,260 55,597 75,113 70,637 71,746
13:00 918,462 926,946 879,452 56,982 54,682 56,073 75,790 71,421 72,686
13:10 929,652 934,930 888,756 57,475 55,298 56,665 76,542 72,394 73,455
13:20 941,050 942,799 899,736 57,941 55,725 57,198 76,958 72,971 74,350
13:30 950,231 951,398 904,224 58,371 56,176 57,743 77,404 73,606 75,027
13:40 959,753 957,294 907,575 58,833 56,719 58,278 77,912 74,323 75,707
13:50 965,117 964,191 919,997 59,311 57,173 58,705 78,645 75,094 76,061
14:00 967,592 968,753 924,060 59,709 57,540 59,231 78,975 75,390 76,617
14:10 974,796 962,552 918,793 60,145 58,067 59,603 79,237 76,175 76,739
14:20 971,706 967,967 920,087 60,300 58,513 60,069 79,197 76,592 76,981

30/06/2007 05/07/2007 11/07/2007 30/06/2007 05/07/2007 11/07/2007 30/06/2007 05/07/2007 11/07/2007
14:30 970,054 966,338 917,177 60,281 59,004 60,475 79,309 76,892 76,895
14:40 965,993 967,686 911,288 60,604 59,255 60,932 79,683 76,862 77,212
14:50 967,355 958,813 901,679 60,724 59,650 61,278 79,912 76,859 77,409
15:00 959,033 957,064 890,515 60,831 59,911 61,644 79,321 77,156 77,487
15:10 955,448 949,391 891,813 61,166 60,001 61,917 79,278 77,134 77,495
15:20 948,047 941,505 882,047 61,336 59,903 62,334 78,961 76,926 77,754
15:30 936,498 930,555 872,739 61,506 60,002 62,535 78,620 76,691 77,297
15:40 921,060 915,870 872,284 61,742 60,250 62,653 78,558 76,782 76,592
15:50 909,117 897,856 862,311 61,856 60,477 62,920 78,138 76,748 75,870
16:00 897,164 887,241 845,826 61,962 60,457 63,158 77,881 75,979 75,724
16:10 879,054 871,284 834,054 61,989 60,573 63,310 77,242 75,391 75,538
16:20 865,378 853,943 822,987 62,038 60,655 63,492 76,665 74,906 75,246
16:30 847,030 843,686 818,753 62,103 60,639 63,722 76,112 74,347 75,304
16:40 825,890 821,020 799,826 62,113 60,626 63,857 75,488 73,667 74,887
16:50 803,274 799,164 774,649 62,045 60,589 64,024 74,768 73,086 74,547
17:00 782,151 774,699 749,886 62,064 60,713 64,136 74,245 72,476 74,308
17:10 757,806 752,311 722,365 62,088 60,515 64,194 73,629 71,572 73,531
17:20 733,217 730,632 704,084 61,978 60,288 66,870 72,787 70,286 72,960
17:30 705,793 707,411 587,612 61,807 60,680 68,705 71,631 69,372 70,982
17:40 682,749 678,268 655,883 63,880 63,615 68,971 70,795 68,896 71,815
17:50 653,274 648,542 637,883 65,685 65,246 69,462 70,534 69,008 71,330
18:00 625,803 617,144 613,100 66,506 66,038 69,640 70,205 69,110 70,795
18:10 599,358 586,231 597,645 66,537 66,237 69,272 69,418 68,477 70,457
18:20 567,599 561,599 556,234 66,151 65,730 68,717 68,199 67,311 69,413

30/06/2007 05/07/2007 11/07/2007 30/06/2007 05/07/2007 11/07/2007 30/06/2007 05/07/2007 11/07/2007
18:30 538,057 534,254 515,211 65,323 64,716 67,860 66,679 65,380 67,650
18:40 505,863 502,291 489,251 64,280 63,455 66,595 65,046 63,399 66,131
18:50 474,485 470,087 468,936 62,909 62,027 65,185 62,968 61,307 64,265
19:00 437,786 436,140 437,033 61,340 60,482 63,793 60,646 59,299 62,369
19:10 407,756 402,963 374,799 59,581 58,735 62,129 58,165 56,815 59,527
19:20 374,375 370,853 365,100 57,777 56,899 60,185 55,608 54,130 56,768
19:30 342,421 339,130 334,321 55,730 54,940 58,191 52,879 51,389 54,222
19:40 310,980 306,000 298,134 53,737 52,985 56,261 49,991 48,495 51,477
19:50 278,254 270,281 269,161 51,722 51,192 54,317 46,947 45,603 48,570
20:00 245,595 239,110 236,468 49,879 49,403 52,562 43,781 42,656 45,669
20:10 214,318 209,819 203,441 48,234 47,670 50,937 41,350 40,186 43,155
20:20 186,151 178,950 176,749 46,696 46,220 49,287 39,849 38,596 41,386
20:30 157,819 150,328 147,913 45,353 44,939 47,781 38,725 37,556 40,128
20:40 130,719 122,478 121,368 44,146 43,820 46,393 37,709 36,714 39,010
20:50 104,930 96,411 97,037 43,077 42,846 45,149 36,830 35,938 37,920
21:00 79,027 73,585 74,468 42,090 41,972 43,969 35,954 35,220 37,065
21:10 55,201 53,341 53,779 41,090 40,916 42,848 35,027 34,388 36,104
21:20 37,140 35,572 36,134 40,173 39,979 41,777 34,096 33,598 35,053
21:30 21,468 19,866 21,535 39,176 39,017 40,796 33,090 32,672 34,060
21:40 10,622 7,920 9,247 38,296 38,068 39,832 32,192 31,681 33,093
21:50 1,171 0,973 0,836 37,460 37,234 38,952 31,344 30,858 32,250
22:00 0 0 0 36,616 36,457 38,131 30,470 30,140 31,434

Para caracterizar dicha instalación se procede a calcular dos parámetros:
1. Qu: calor útil transmitido al agua que circula por los colectores [kJ]:
Qu = m& c ⋅ C p ⋅ ∆Tc ⋅ ∆t ec [3.1]
Donde: m& c = caudal másico de los colectores [kg/s]:
m& c = ρ ⋅ Qc ec [3.2]
Qc = caudal de agua que circula por los colectores = 0,15 [l/s].
ρ = densidad del agua = 1000 [kg/m3].
Cp = calor específico a presión constante del agua = 4,18 kJ/kg
∆Tc = diferencia de temperatura entre la salida y la entrada del

colector [K].
∆t = incremento de tiempo [s] = 10 [min]=600 [s].
2. Rendimiento de los colectores (ηc): cuya expresión es la dada por la ec

[2.2]:
t2
∫ Q dt
t1
u
ηc = t2
ec [2.2]
A ⋅ ∫ I gl dt
t1
Donde: Igl = radiación global recibida por la unidad de tiempo y unidad

de área del captador [W/m2].
A = Área de los captadores [m2] = 20 [m2].
Qu= calor útil [kW].
En la Tabla 3.20 se representan Qu y ηc, calculados a través de las

expresiones anteriores y con los datos recogidos en la Tabla 3.19.

Tabla 3.20. Calor útil y rendimiento de los colectores.

Qu [kJ] ηc [%]
Tiempo [min]
30/06/2007 05/07/2007 11/07/2007 30/06/2007 05/07/2007 11/07/2007
7:00 -2484,430 -2490,218 -2373,960 0 0 0
7:10 -1686,735 -1643,705 -1583,060 0 0 0
7:20 -1475,107 -1460,638 -1389,927 0 0 0
7:30 -1333,560 -1326,514 -1248,255 0 0 0
7:40 -1260,081 -1206,482 -1130,991 0 0 0
7:50 -1181,067 -1103,939 -1044,175 0 0 0
8:00 -1066,823 -1030,083 -961,009 0 0 0
8:10 -988,060 -950,566 -897,092 0 0 0
8:20 -872,809 -846,261 -805,622 0 0 0
8:30 -749,884 -752,274 -707,734 0 0 0
8:40 -621,925 -654,387 -595,251 0 0 0
8:50 -448,798 -496,106 -474,339 0 0 0
9:00 -223,455 -295,550 -295,802 0 0 0
9:10 32,336 -77,379 -90,590 0,802 0 0
9:20 278,312 172,750 140,918 6,280 3,774 3,328
9:30 483,021 399,351 344,997 10,011 8,025 7,461
9:40 651,241 544,672 520,767 12,457 10,113 10,341
9:50 765,108 655,016 662,943 13,659 11,335 12,192
10:00 871,677 769,511 741,705 14,562 12,493 12,806
10:10 973,968 865,260 772,405 15,340 13,296 12,497
10:20 1151,499 869,915 828,772 17,162 12,657 12,591

Continuación de Tabla 3.20. Calor útil y rendimiento de los colectores.

Qu [kJ] ηc [%]
Tiempo [min]
30/06/2007 05/07/2007 11/07/2007 30/06/2007 05/07/2007 11/07/2007
10:30 1195,410 908,542 867,650 17,154 12,589 12,667
10:40 1399,867 992,967 938,361 19,109 13,141 12,947
10:50 2669,133 1033,984 1430,063 34,723 13,052 18,902
11:00 3563,960 1803,873 2785,012 44,724 21,916 35,420
11:10 3946,577 3034,261 3400,395 47,624 35,455 42,271
11:20 4417,771 3577,675 3716,579 51,433 40,453 44,747
11:30 4770,316 4006,341 4056,795 53,983 43,944 47,060
11:40 5055,549 4372,350 4466,031 55,748 46,657 49,421
11:50 5339,146 4653,934 4662,489 57,557 48,296 50,930
12:00 5715,472 4925,452 5095,308 59,132 49,746 54,203
12:10 6061,224 5243,020 5366,826 61,097 51,837 55,460
12:20 6405,843 5533,034 5584,997 62,472 53,443 56,565
12:30 6649,681 5782,156 5719,372 63,659 54,748 57,376
12:40 6866,594 6005,486 5840,285 64,290 55,630 57,717
12:50 7021,855 6160,998 6075,315 64,609 56,207 58,334
13:00 7075,580 6297,009 6249,827 64,198 56,611 59,221
13:10 7173,090 6431,510 6316,511 64,299 57,326 59,226
13:20 7154,217 6488,003 6452,522 63,353 57,347 59,763
13:30 7160,256 6557,204 6502,473 62,794 57,435 59,927
13:40 7177,619 6622,504 6556,700 62,322 57,649 60,203
13:50 7273,619 6742,033 6529,272 62,804 58,270 59,142

Continuación de Tabla 3.20. Calor útil y rendimiento de los colectores.

Qu [kJ] ηc [%]
Tiempo [min]
30/06/2007 05/07/2007 11/07/2007 30/06/2007 05/07/2007 11/07/2007
14:00 7247,700 6714,982 6540,722 62,420 57,763 58,985
14:10 7182,400 6812,491 6446,734 61,401 58,979 58,471
14:20 7109,173 6801,293 6362,561 60,968 58,553 57,626
14:30 7158,621 6729,325 6177,355 61,497 58,031 56,127
14:40 7177,619 6623,888 6124,259 61,919 57,042 56,004
14:50 7218,762 6474,037 6068,647 62,186 56,268 56,087
15:00 6956,051 6487,751 5960,191 60,443 56,490 55,775
15:10 6813,749 6445,476 5860,542 59,429 56,576 54,762
15:20 6630,557 6404,207 5801,155 58,283 56,684 54,808
15:30 6438,556 6278,262 5553,417 57,293 56,223 53,027
15:40 6326,074 6219,127 5243,775 57,235 56,587 50,096
15:50 6125,140 6121,240 4871,601 56,146 56,813 47,079
16:00 5988,752 5839,152 4727,412 55,627 54,844 46,576
16:10 5738,245 5574,554 4600,083 54,398 53,317 45,961
16:20 5502,711 5361,417 4421,671 52,989 52,320 44,773
16:30 5270,197 5156,834 4357,251 51,850 50,936 44,348
16:40 5031,895 4906,076 4149,398 50,772 49,797 43,232
16:50 4786,547 4701,116 3958,656 49,657 49,021 42,585
17:00 4582,594 4425,194 3826,797 48,825 47,601 42,526
17:10 4341,776 4159,464 3512,374 47,745 46,074 40,519
17:20 4066,483 3760,994 2291,046 46,217 42,897 27,116

Continuación de Tabla 3.20.Calor útil y rendimiento de los colectores.
Qu [kJ] ηc [%]
Tiempo [min]
30/06/2007 05/07/2007 11/07/2007 30/06/2007 05/07/2007 11/07/2007
17:30 3695,693 3269,795 856,704 43,635 38,518 12,150
17:40 2601,316 1986,562 1070,094 31,751 24,407 13,596
17:50 1824,507 1415,217 702,827 23,274 18,185 9,182
18:00 1391,563 1155,399 434,203 18,530 15,601 5,902
18:10 1083,808 842,487 445,778 15,069 11,976 6,216
18:20 770,392 595,000 261,704 11,311 8,829 3,921
18:30 509,946 249,752 -79,027 7,898 3,896 0
18:40 288,000 -21,264 -174,386 4,744 0 0
18:50 22,144 -271,141 -345,878 0,389 0 0
19:00 -261,075 -445,149 -535,613 0 0 0
19:10 -532,468 -722,581 -978,749 0 0 0
19:20 -816,065 -1041,533 -1285,622 0 0 0
19:30 -1072,862 -1335,950 -1493,099 0 0 0
19:40 -1409,555 -1689,125 -1799,847 0 0 0
19:50 -1796,072 -2102,442 -2162,206 0 0 0
20:00 -2293,939 -2538,155 -2592,886 0 0 0
20:10 -2589,615 -2815,461 -2927,566 0 0 0
20:20 -2575,775 -2868,431 -2972,232 0 0 0
20:30 -2493,237 -2777,841 -2878,874 0 0 0
20:40 -2421,646 -2673,285 -2777,463 0 0 0
20:50 -2350,055 -2598,799 -2719,586 0 0 0
21:00 -2308,409 -2540,294 -2597,416 0 0 0
21:10 -2280,854 -2456,121 -2537,022 0 0 0
21:20 -2285,887 -2400,760 -2529,473 0 0 0
21:30 -2289,787 -2386,794 -2534,128 0 0 0
21:40 -2296,330 -2402,647 -2535,135 0 0 0
21:50 -2300,859 -2398,495 -2521,043 0 0 0

Se detallan, en la Tabla 3.21 los valores medios diarios tanto de Qu como

del rendimiento de los colectores:
Tabla 3.21 Valores medios diarios de Qu y ηc.

Radiación Global Qu diario Qu diario
ηc diario [%]
diaria [W/m2] [kW] [kJ]
30/06/2007 49895,070 431,968 259180,913 43,288
05/07/2007 50072,562 388,267 232959,900 38,770
11/07/2007 47644,877 355,420 213252,250 37,299
Promedio 49204,170 391,885 235131,021 39,786
Se representa gráficamente el comportamiento de la instalación solar en

estos tres días:
Figura 3.9. Radiación global vs Tiempo.

La radiación global posee la característica forma de campana, tal como se
detalla en el capítulo 2 (pág 58). En esta época del año amanece en torno a las
7:00, por tanto, la radiación captada es mínima. Al avanzar el día y debido al
movimiento de la Tierra respecto del Sol, llega mayor proporción de radiación y
se consigue un máximo en torno a los 960 W/m2, entre 12:45-16:15. Este
parámetro decae con la venida del atardecer y de la noche.

Figura 3.10. Rendimiento de los colectores frente Tiempo.
Entre las 7:00-9:20 el rendimiento de los colectores es nulo, debido a que

Qu es negativo ya que la radiación solar es tan baja que no es suficiente para
calentar el agua que circula por los colectores, por lo que se producen
pérdidas. A partir de este momento, se produce una subida del rendimiento
coincidiendo con el incremento paulatino de la radiación captada por los
colectores, en torno a las 10:50 se produce un rápido crecimiento de este
parámetro llegando a sus valores máximos, entre el 65-55%, en el intervalo
donde la radiación global es máxima (960 W/m2, 12:45-16:45). El rendimiento
de los colectores comienza a decaer bruscamente en torno a las 16:45 debido
al decrecimiento de la radiación solar, y a las 18:00 vuelve a ser nulo por el
mismo efecto descrito al principio.
Se concluye que el existe una relación directamente proporcional entre ηc

e Igl, a mayor radiación aumenta el rendimiento y viceversa.

Figura 3.11. Rendimiento de los colectores frente Temperatura de salida del colector.
Si el rendimiento de los colectores es nulo la temperatura de salida varía

entre los 15-25ºC, esto es debido a que a medida que se incrementa la
radiación solar, el agua se calienta pero no a la suficiente temperatura porque
se producen pérdidas en la instalación. A medida que este parámetro aumenta,
la temperatura también, llegando en un rendimiento del 13% a 55ºC. Se
observa que en el intervalo de máxima temperatura 75-80ºC, se alcanza el
mayor rendimiento 60-65%. Posterior a este punto las dos variables comienzan
a disminuir, destacando que la temperatura a rendimiento nulo decrece desde
los 65ºC hasta los 15ºC, este fenómeno es debido a que poco a poco el agua
de los colectores se va enfriando ya que la radiación es menor.
Se destaca que a medida que aumenta el rendimiento de los colectores,

la temperatura de salida del agua que circula por la instalación solar térmica
también crece, como era lógico. Por lo tanto existe una relación directamente
proporcional entre el rendimiento de los colectores y la temperatura de salida
del fluido que circula a través de ellos.

Figura 3.12. Calor útil frente a rendimiento de los colectores.
El calor útil negativo representa las pérdidas que existen en la instalación

y se corresponde con un rendimiento de colector nulo, esto significa que la
radiación captada no es lo suficientemente alta para elevar la temperatura del
agua que circula por los colectores. A medida que se eleva el rendimiento el
calor útil crece, llegando a un máximo entre 6500-6900 kJ en un intervalo de ηc
60-65%. Llegado a este punto, comienzan a decrecer ambas variables, debido
a una disminución de la radiación, y se regresa a la zona de pérdidas donde el
uso de la planta solar térmica es inadecuado.
Por lo tanto, existe una relación directamente proporcional entre el calor

útil transferido al fluido y el rendimiento de los colectores, como era de esperar.

Figura 3.13 Radiación global diaria recibida entre el 30/06 y 11/07.
El comportamiento de la radiación global es oscilante, posee un mínimo

de 40000 W/m2 y un máximo de 50000 W/m2. Este fenómeno se debe tanto al
movimiento de rotación y traslación de la Tierra respecto del Sol como a las
condiciones climatológicas de la zona en la cual se encuentra localizada esta
instalación.
3.3 Máquina de absorción.
La máquina de absorción posee una potencia de 35 kW y su par de trabajo

es una disolución LiBr- H2O. Es alimentada por el agua procedente del tanque,
la cual al poseer una temperatura inferior a 65ºC, no fue capaz de poner en
funcionamiento este equipo por lo que no se produjo frío.

4. SIMULACIÓN

4. Simulación.
4.1 Desarrollo del modelo.
Se realiza una simulación de la instalación de frío solar con la aplicación

informática MATHCAD 14. Los objetivos de esta simulación son:
- Predecir el comportamiento de la instalación de frío solar destinada a

suministrar la demanda de climatización de una vivienda de 90 m2.
- Analizar la influencia de la superficie solar y del volumen del tanque de
almacenamiento.
- Verificar con los datos experimentales del capítulo 3 si la simulación es
correcta.
El programa de cálculo es “Instalación”, una hoja de MATHCAD, este a su

vez consta de dos subprogramas “Adquisición” y “Caract Colector1”. A
continuación se realiza una descripción detallada de cada una de las partes
que constituyen esta simulación, para el día 30/06/2007.
4.1.1 Adquisición.
Este subprograma, es una hoja de MATHCAD cuya misión es la

adquisición de los datos recogidos por la estación meteorológica. Se divide en
cinco partes:
Antes se debe definir la longitud de los vectores, la variable C, será el

grado centígrado ya que la unidad que toma MATHCAD de temperaturas es el
grado Kelvin, y n es el día del año cuyo intervalo es [1,365]:
i := 0 .. 143 C := 1 ⋅K n := 181

1) Tiempos solares:
Se define la variable tiempo:
〈 〉
tiempo := DATOS 89 ⋅min ec [4.1]
Siendo DATOS<89>, la columna del archivo en el cual se encuentra el

tiempo en minutos. Se debe corregir el tiempo local teniendo en cuenta la
excentricidad de la órbita terrestre, ya que se trata de una elipse cuyo centro no
es el sol, de acuerdo con U. Eicker, (2003)/Pág: 19/(2.8) :
ec [4.2]
Por lo que queda la corrección global: zona horaria local + órbita + zona
horaria legal + horario estacional como:
φ := 40.41 ⋅deg Latitud
L := 3.5 ⋅deg Longitud
Lstd := 0 ⋅deg Longitud del meridiano de referencia para la zona horaria
⎡ min ⋅ L − L + E⎤
∆h := −⎢4 ⋅ ( std ) ⎥ ec [4.3]
⎣ deg ⎦
∆h = 17.2586 ⋅min
tiempo := tiempo + ∆h ec [4.4]

2) Velocidad y dirección del viento: Se definen estas dos variables cómo:
〈 〉 m
Vmed_viento:= DATOS 61 ⋅ ec [4.5]
sec
〈 〉 ec [4.6]
Dir_viento:= DATOS 62
3) Temperaturas: la temperatura media del aire posee la siguiente expresión:
〈 〉
tmed := DATOS 63 ⋅C ec [4.7]
4) Insolaciones: La insolación global incidente es:
〈 〉 watt ec [4.8]
Imed := DATOS 67 ⋅
2
m
5) Cálculo de la radiación sobre el plano inclinado de los colectores:
En primer lugar se necesita saber la latitud e inclinación de los colectores:
Latitud: φ := 40.41 ⋅ deg
Inclinación del panel:

β i := 45 ⋅ deg

A continuación se calcula la Constante Solar (G0ni), es decir la irradiación

solar exterior de la atmósfera terrestre, y la declinación (δi), es la inclinación del
eje de rotación de la Tierra respecto al plano de la elíptica, de acuerdo con U.
Eicker (2003)/Pág: 17/(2.4) :
⎛
G0n := 1353 ⋅⎜ 1 + 0.033 ⋅cos ⎜
⎛ 360 ⋅n ⋅deg⎞ ⎞ ⋅ watt ec [4.9]
i ⎝ ⎝ 365 ⎠ ⎠ m2
δi := 23.45 ⋅sin⎜ 360 ⋅

⎛ 284 + n ⎞
⋅deg ⋅deg ec [4.10]
⎝ 365 ⎠
Cuyos valores son: G0ni= 1308 W/m2 y δi= 23,18º.
Se muestran los ángulos horarios de salida y puesta de sol sobre el plano

horizontal (ωszi) e inclinado (ωsTi), de acuerdo con U.Eicker (2003)/Pág:
21/(2.14) :
ω sz := acos ( −tan( φ ) ⋅tan( δ i) ) ec [4.11]

i
ω sT := acos ( −tan( φ − β i) ⋅tan( δ i) ) ec [4.12]

i
Finalmente se le asigna un ángulo horario a cada serie de datos tomados.

Hallados estos parámetros se procede al cálculo de la radiación total y difusa
en el plano inclinado para cada serie de datos, para ello se obtiene en primer
lugar el índice de transparencia atmosférica (kt), la cual determina la oposición
de la atmósfera al paso de la radiación solar debido a los obstáculos que se
pueden encontrar en ella (partículas principalmente), de acuerdo con U. Eicker
(2003)/Pág: 27/(2.24) :

ec. [4.13]
La componente difusa (Id), aquella componente de la radiación solar que

llega al colector rebotando en diferentes obstáculos, se calcula según la
correlación de Orgill y Hollands (Canadá 1977):
Id ( kT , I , i ) := ⎛ 1.0 − 0.249 ⋅kT if 0 ≤ kT < 0.35 ⎞ ⋅I

⎜
⎜ 1.557 − 1.84 ⋅kT if 0.35 ≤ kT < 0.75 ⎟ ec [4.14]
⎜ 0.6584 ⋅kT − 0.3168 if 0.75 < kT ≤ 2 ⎟
⎜
⎝ 1 otherwise ⎠
Para obtener la componente directa se calcula la relación Rb= In/IT

(In=radiación sobre el plano horizontal e IT= radiación sobre el inclinado):
cos( φ − β i) ⋅cos( δ i) ⋅cos(ω ) + sin( φ − β i) ⋅sin( δ i)

Rb (ω , i) := if ω < ω s ec [4.15]
cos(φ ) ⋅cos( δ i) ⋅cos(ω ) + sin(φ ) ⋅sin( δ i) i
0 otherwise
Finalmente se obtiene la radiación total solar sobre el plano inclinado

(Itilted) para cada serie de datos y la energía total por metro cuadrado tanto del
plano inclinado como la total:
⎛ 1 + cos ( β i) ⎞ ⎛ 1 − cos ( β i) ⎞
ITilted ( I , Id , Rb , i ) := ( I − Id ) ⋅Rb + Id ⋅⎜ + I ⋅0.2 ⋅⎜
⎝ 2 ⎠ ⎝ 2 ⎠
ec [4.16]

Itilted := for i ∈ 1 .. 141

kT1 ← kT ( ω i −1 , ω i , Imedi , i )
Id1 ← Id ( kT1 , Imedi , i )
ec [4.17]
Rb1 ← Rb ( ω i , i )
Itiltedi ← ITilted ( Imedi , Id1 , Rb1 , i )
Itilted143 ← 0
Itilted
kWh
Etilted := ∑ ( Itiltedi ⋅600 ⋅sec) Etilted = 7.0517 ⋅
m
2
ec [4.18]
Emed := ∑ ( Imedi ⋅600 ⋅sec) Emed = 8.3161 ⋅

kWh
2
ec [4.19]
i m
4.1.2 Caract Colector1.
Este subprograma, es una hoja de MATHCAD en la cual se han recogido

las dimensiones y propiedades de los materiales que componen la instalación.
Se compone de tres partes:
1) Calor específico (Cp) y densidad de los materiales y del agua (ρ): se detalla
el número de colectores (Ncol):

2) Características de los colectores: se definen las variables geométricas

detalladas en la Tabla 4.1:
Tabla 4.1. Variables geométricas de los colectores.

Nombre Símbolo
Área de un colector aco
Área total Aco
Área efectiva de un colector aef
Área efectiva total Aef
Largo del colector lco
Ancho del colector anco
Altura del colector hco
Espesor de la caja ealum
Se calcula la masa total de cada componente de los colectores:
-Cubierta de cristal de 4 mm de espesor:
ec [4.20]
mcov := A co⋅ 0.004⋅ m⋅ ρcristal
-Placa absorbente de AlOx de 0.3 mm de espesor:
mpl := A ef ⋅ 0.0003⋅ m⋅ ρAlOx

ec [4.21]
-Aislamiento térmico de gas Criptón, caras laterales y cara posterior:
maiscol := ( h co − 0.0003m) ⋅ lco⋅ an co⋅ ρaiscol⋅ Ncol

ec [4.22]

Caja de aluminio de 1mm de espesor
mbox := ⎡⎣2⋅ h co⋅ ( lco + an co) + lco⋅ an co⎤⎦ ⋅ ealum⋅ ρalum⋅ Ncol
ec [4.23]
-Fluido caloportador contenido en el colector:
kg
mfco := A co⋅ 0.64⋅
2 ec [4.24]
m
Se definen los parámetros que caracterizan el funcionamiento de

operación de los colectores:
-Gasto másico nominal:
ec [4.25]
gmNom := 0.012⋅ ⎛ aco⋅ ⎞
kg
⎜ 2
⎝ sec ⋅ m ⎠
GmNom := gmNom⋅ Ncol
ec [4.26]
-Parámetros de la curva de normalización referida a la temperatura media del fluido en el colector:
Fr τα := 0.82 ec [4.27]
watt
FrUo := 2.2⋅ ec [4.28]
2
m ⋅C
3) Características de las tuberías: sólo se tiene en cuenta los circuitos de

hidráulicos de ida, recogidas en la Tabla 3.8-3.13:
Circuito primario:
-Tuberías de cobre sanitario de 24*26 mm de diámetro:
0.026
2
− 0.024
2 ec [4.29]
2
m tub1 := ρ cobre ⋅ π ⋅ ⋅ m ⋅ 29.05 m
4
-Aislamiento de las tuberías, SH/Armaflex, espesor 20 mm:
⎛ 0.066 2 − 0.026 2 ⎞ 2 ec [4.30]

maistub1 := π ⋅ ⎜ ⋅ m ⋅ 29.05 ⋅ m⋅ ρ aistub
⎝ 4 ⎠

-Fluido caloportador:
2
0.024 2 ec [4.31]
mf1 := π ⋅ ⋅ m ⋅ 29.05 m⋅ ρ agua
4
Circuito secundario:
0.026
2
− 0.024
2
2
ec [4.32]
m tub2 := ρ cobre ⋅ π ⋅ ⋅ m ⋅ 3.17 ⋅ m
4
-Aislamiento de las tuberías, SH/Armaflex, espesor 20mm:
⎛ 0.066 2 − 0.026 2 ⎞ 2 ec [4.33]

maistub2 := π ⋅ ⎜ ⋅ m ⋅ 3.17 ⋅ m⋅ ρ aistub
⎝ 4 ⎠
2
0.026 2 ec [4.34]
m f2 := π ⋅ ⋅ m ⋅ 3.17 m ⋅ ρ agua
4
Circuito terciario:
2 2
0.035 − 0.033 2
m tub3 := ρ cobre ⋅ π ⋅ ⋅ m ⋅ 17.78 ⋅ m ec [4.35]
4
-Aislamiento de las tuberías, SH/Armaflex, espesor 20mm:
⎛ 0.0752 − 0.0352 ⎞ 2 ec [4.36]

maistub3 := π ⋅ ⎜ ⋅ m ⋅ 17.78⋅ m⋅ ρaistub
⎝ 4 ⎠
2 ec [4.37]
0.033 2
mf3 := π ⋅ ⋅ m ⋅ 17.78 m⋅ ρ agua
4

4) Características del tanque de almacenamiento: recogidas en la Tabla 3.6:
V tan := 2 m
3 ec [4.38]
-Volumen total almacenado:
-Carcasa de acero cilíndrica de 42 mm de espesor

h tane := 1.634 m
D tane := 1.284 m h tani := 1.55 m
D tani := 1.2 m
⎡ π ⋅ ⎛ D 2 − D 2 ⎞ + π ⋅ (D ⋅ h ⎤
mtan := ⎢ 4 ⎝ tane tani ⎠ tane tane − D tani ⋅ h tani )⎥ ⋅ 0.042 m ρ acerog ec [4.39]
⎣ ⎦
-Aislamiento de poliuretiano con un espesor de 35 mm:
h aistan := 1.704 m
D aistan := 1.354 m
ec [4.40]
maguatan := Vtan⋅ ⎛ 1000⋅

kg ⎞ ec [4.41]
⎜ 3
⎝ m ⎠

4.1.3 Instalación.
Este programa se trata de una hoja de MATHCAD, que emplea los datos
recogidos en “Adquisición” y “Caract Colector1”. El objetivo es predecir el
comportamiento de la instalación de frío solar ante una carga térmica
establecida para una vivienda. Comprende tres partes:
1) Datos de partida: se definen los datos iniciales:
Siendo: DATOS= archivo de donde proceden los datos experimentales,

empleados en “Adquisición”.

2) Cálculo de temperaturas:
Colector.
La ecuación 4.42, muestra las pérdidas del colector en caso de que se
formara una hipotética capa de hielo sobre la superficie de los colectores. En
este caso será cero, debido a que la estación es verano y la temperatura media
del ambiente no va a estar por debajo de 0ºC. Se muestra la ecuación porque
es necesaria para el cálculo de temperaturas[29]:
Qhielo :=
⎛ ⎛
Aco ⋅ 0.003 ⋅ m⋅ ⎜650 ⋅ exp⎜0.227 ⋅ min⎜
⎛ tmed ⎞⎞⎞ ⋅ kg ⋅329 ⋅ kjoule if min( tmed) < 0 ⋅C
⎝ ⎝ ⎝ C ⎠⎠⎠ m3 kg
0 ⋅ kjoule otherwise
ec [4.42]
Donde: 0,003m = estimación del espesor de la capa de hielo que se

formaría.
329 kJ/kg = calor latente de cambio de fase del hielo
650 exp (0.227·min (tmed/C) kg/m3= densidad del hielo.
Por lo tanto el calor que reciben los colectores solares es:
Qpl ( i , hielo ) := Aco ⋅ Itilted ⋅ 600 ⋅ sec ⋅ hielo ec [4.43]

i
Donde: 600sec es debido a que los datos se toman cada diez minutos.
hielo = es una variable que indica si hay o no presencia de hielo.

Si toma el valor 1 es que no existe hielo y si es 0 es que existe
hielo y toda la radiación solar se emplea en derretirlo.

Las pérdidas ópticas, por radiación y convección en los colectores

presenta la siguiente expresión:
( )
FrUo ec [4.44]
Pco ( tave , i) := Aco ⋅ ⋅ 600 ⋅ sec ⋅ tave − tmed
Fr τα i
Siendo: tave= temperatura media entre la temperatura de entrada (tinco) y

salida del colector (toutco).
tmedi = temperatura media del aire exterior.
Las pérdidas por efecto del viento presentan la siguiente ecuación, tanto
para la superficie del colector como para la caja de aluminio[29]:
sec watt ⎛ ave ⎞

t + tmed
i
Pcov ( tave , i) := Aco ⋅ 3 ⋅ Vmed ⋅ ⋅ ⋅ ⎜ − tmed ⋅ 600 ⋅ sec ec [4.45]
m ⋅C ⎝
i m 2 2 i⎠
⎡ 3 ⋅ tave + tmed
sec watt ⎛ ⎞ ⎤ ec [4.46]
Pbox ( tave , i) := ⎢Aco ⋅ 3 ⋅ 1.2 ⋅ Vmed ⋅ − tmed ⋅ 600 ⋅ sec⎥
i
⋅ ⋅⎜
⎢ i m
m ⋅C ⎝
2 4 i⎠ ⎥
⎣ ⎦
Siendo: Vmedi = velocidad del viento.
1,2 = el área de la caja de aluminio es 1,2 veces que la cubierta.
tAVE = temperatura tave en el instante i-1.
Los índices fraccionarios son una aproximación.
Finalmente se tienen en cuenta las pérdidas debidas a los calores

almacenados por las diferentes partes que forman el colector:
HCE co ( t ave , tAVE , i) := mpl ⋅ Cp AlOx ⋅ ( tave − tAVE) ...

⎛ tave + tmed i tAVE + tmed i−1 ⎞
+ mcov ⋅ Cp cristal ⋅ ⎜ − ...
⎝ 2 2 ⎠
⎛ 3 ⋅ tave + tmed i 3 ⋅ tAVE + tmed i−1 ⎞
+ mbox ⋅ Cp alum ⋅ ⎜ − ...
⎝ 4 4 ⎠
⎛ 7 ⋅ tave + tmed i 7 ⋅ tAVE + tmed i−1 ⎞
+ maiscol ⋅ Cp aiscol ⋅ ⎜ − ...
⎝ 8 8 ⎠
+ mfco ⋅ Cp agua ⋅ ( tave − tAVE)
ec [4.47]

Circuito primario (colectores solares- intercambiador de calor):
En este tramo, el diámetro externo es de 26 mm posee un espesor 1 mm

y el aislante posee 20 mm de espesor.
Las pérdidas por efecto del viento se calculan de la siguiente forma:

aire
tmed
aislante
tubería
toutco
agua salida colector
Figura 4.4. Esquema de la tubería.
Las resistencias térmicas son:
R tt = R t , hi + R t , ci + R t , ce + R t , he
1
R t , hi =
hint ⋅ Aint
⎛r ⎞
ln ⎜⎜ ext ⎟⎟
R t , ci = ⎝ rint ⎠
2 ⋅ π ⋅ L ⋅ k cobre
⎛r ⎞
ln ⎜⎜ aistub ⎟⎟
= ⎝ rext ⎠
R t , ce
2 ⋅ π ⋅ L ⋅ k aistub
1
R t , he =
h ext ⋅ Aext

Siendo: Rt,hi = resistencia térmica de convección interna.
Rt,ci = resistencia térmica de conducción tubería de cobre.
Rt,ce = resistencia térmica de conducción del aislante.
Rt,he = resistencia térmica de convección externa.
Se procede al cálculo de los coeficientes de convección interna y externa

(hint y hext):
A) Flujo interno:
hint se obtiene del número de Nusselt:
watt
Nu⋅ 0.668 ⋅
m⋅ C
hint := ec [4.48]
0.024 ⋅ m
Donde: 0.668 W/m C = coeficiente de conducción del agua.
0.024m = diámetro interior del tubo.
Se emplea la siguiente ecuación, de acuerdo con F.Incropera (Capítulo

8/pág 445):
0.8 0.3
Nu := 0.023 ⋅ Rey ⋅ Pr
(vint ⋅ρ agua ⋅0.024 ⋅m) ec [4.50]

Rey :=
−6
389 ⋅ 10 ⋅ Pa ⋅ sec
−6 ec[4.51]
389 ⋅ 10 ⋅ Pa ⋅ sec ⋅ Cpagua
Pr :=
watt
0.668 ⋅
m⋅ C

B) Flujo externo:
hext se obtiene del número de Nusselt:
− 3 watt
Nuex ⋅ 26.3 ⋅ 10 ⋅
i m⋅ C ec [4.52]
hex :=
i 0.066 ⋅ m
Donde: 26.3*10-3 W/m C = coeficiente de conducción del aire.
0.066m = diámetro exterior del tubo.
Se emplea la siguiente ecuación, de acuerdo con F.Incropera (Capítulo

7/pág 369):
Nu ex := 0.193 ⋅ Re ex
i ( i )0.618 ⋅ Pr ex 0.33 ec [4.53]
kg
Vmed ⋅ 1.16 ⋅ ⋅ 0.066 ⋅ m
i 3
m
Reex :=
i −7
184.6 ⋅ 10 ⋅ Pa ⋅ sec ec [4.54]
−7 joule
184.6 ⋅ 10 ⋅ Pa ⋅ sec ⋅ 1007 ⋅
kg ⋅ C
Prex := ec [4.55]
− 3 watt
26.3 ⋅ 10 ⋅
m⋅ C

Se define un coeficiente global de transferencia de calor del tubo y se

calcula el área exterior del tubo:
1
Utub :=
i ⎛ 0.013 ⎞
0.012 ⋅ m⋅ ln ⎜ 0.033 ⋅ m⋅ ln ⎛⎜
0.033 ⎞
1
+
⎝ 0.012 ⎠ + ⎝ 0.013 ⎠ + 0.033
− 5 watt watt watt hint ⋅ 0.013
hex + 10 ⋅ 63.9 ⋅ 0.035 ⋅
i
C ⋅m
2 m⋅ C m⋅ C
ec [4.56]
Aext := π ⋅ 0.066 ⋅ m ⋅ 29.05 m ec[4.57]
Por lo tanto la energía disipada por el tubo (Pins(tave,i)), la temperatura en

la pared del aislante del tubo (taistub(tave,i)) y los calores almacenados por cada
uno de los elementos del circuito son:
(
Pins ( tave , i) := Aext ⋅ Utub ⋅ 600 ⋅ sec ⋅ tave − tmed
i i ) ec [4.58]
Pins ( tave , i)
taistub ( tave , i) := tmed + ec [4.59]
i ⎛h + 10− 5 ⋅ watt ⎞ ⋅A ⋅600 ⋅sec
⎜ exi 2
ext
⎝ C ⋅m ⎠
HCEins ( tave , tAVE, i) := mtub1 ⋅ Cpcobre ⋅ ( tave − tAVE) ...

⎛ tave + taistub ( tave , i) tAVE + taistub ( tAVE, i) ⎞
+ maistub1 ⋅ Cpaistub ⋅ ⎜ − ..
⎝ 2 2 ⎠
+ mf1 ⋅ Cpagua ⋅ ( tave − tAVE)
ec[4.60]

Tanque.
Se define un coeficiente global de transferencia de calor y el área exterior

del tanque:
watt ec [4.60]
Utan := 0.68 ⋅
2
m ⋅C
(
π ⋅ Daistan )2 ec [4.61]
Aexttan := ⎡⎣π ⋅ ( Daistan ) ⋅ haistan⎤⎦ +
2
Por lo tanto la energía disipada por el tubo (Ptan(tave,i)), la temperatura en

la pared del aislante del tanque (taistan(tave,i)) y los calores almacenados por
cada uno de los elementos del circuito son:
(
Ptan ( ttan , i) := Aexttan ⋅ Utan ⋅ 600 ⋅ sec ⋅ ttan − tmed
i ) ec [4.58]
Ptan ( ttan , i)
taistan ( ttan , i) := tmed +
i Utan ⋅ Aexttan ⋅ 600 ⋅ sec
ec [4.59]
HCEtan ( ttan , tTAN , i) := ( maguatan ⋅ Cpagua + mtan ⋅Cpacerog ) ⋅( ttan − tTAN) ...
⎛ ttan + taistan ( ttan , i) tTAN + taistan ( tTAN , i) ⎞
+ maistan ⋅Cpaistub ⋅⎜ −
⎝ 2 2 ⎠
ec[4.60]
Siendo: ttan= temperatura del tanque.
tTAN = la temperatura ttan en el instante i-1.

Finalmente se define Qs, el calor cedido al tanque como:
ec[4.61]
Siendo: tintan= temperatura en la entrada del tanque.
Circuito secundario y terciario.
Los calores acumulados por cada componente de estos circuitos son:
ec[4.62]
Calor útil (Qu).
Qu es el calor necesario para que la máquina funcione. Se ha simulado de

la siguiente forma:
Qe
Qu = ec[4.63]
COP

Siendo Qe= demanda frigorífica, hallada según [38]:
Qe = U ⋅ A ⋅ (t med − t confot ) ec[4.64]
Donde: U*A=360 W/ºC
tconfort= 24ºC.
tmed= temperatura del aire de los datos experimentales.
COP = 0,36, valor medio experimental de la máquina de

absorción.
Realizado esto se inserta la hoja de datos en el MATHCAD:
Sistema de ecuaciones para el cálculo de temperaturas.
El sistema de ecuaciones se encuentra formado por aquellas ecuaciones

de balances energéticos en cada una de las partes de la instalación:
colector(1), circuito primario(2) y el tanque(3):
(1)
(2)
(3)
Se especifican las condiciones iniciales para resolver el sistema de

ecuaciones:
ttan := 40 ⋅ C tinco := tmed1

tave := tmed1 tintan := tmed1

El sistema de ecuaciones se resuelve con el siguiente algoritmo:
temps := temp0 , 0 ← tmed0

temp0 , 1 ← tmed0
temp0 , 3 ← tmed0
temp0 , 2 ← 60 ⋅ C
Qs ← 0 ⋅ joule
hielo ← 0
for j ∈ 1 .. 141
(
t ← temps j , temp
j −1 , 0
, temp
j −1 , 2
, hielo )
temp ← t0
j,0
temp ← t1
j,1
temp ← t2
j,2
temp ← t3
j,3
Qs ← Qs + Aco ⋅ Fr τα ⋅ Itilted ⋅ 600 ⋅ sec
j
hielo ← 1 if Qs > Qhielo
temp i := 0 .. 141
( 〈0〉
tave := temps i
i
) ( 〈2〉
ttan := temps i ) (
〈1〉
tintan := temps i )
i i
( 〈3〉
tinco := temps i
i
) toutco := 2 ⋅ tave − tinco
i i i
3) Resultados: se obtienen la simulación de las temperaturas a lo largo del día,

así como el rendimiento total e instantáneo de la instalación y los valores
totales de los distintos parámetros definidos en el programa, los cuales se
comentarán el apartado siguiente.

4.2 Resultados de la simulación y comparación con los datos experimentales.
Antes de proceder al análisis de los resultados obtenidos, se ha de

destacar que se procede a una validación parcial del modelo, debido a que el
período de tiempo analizado es corto. Los días analizados son los mismos que
los del capítulo 3.
4.2.1Temperaturas.
Se muestran los resultados de la simulación en el cálculo de las

temperaturas de entrada y salida del colector junto con los experimentales
recogidos en los días 30/06, 05/07 y 11/07.
a) 30/06/2007
Figura 4.5. Evolución de temperaturas simulación (superior) y experimental (inferior).

b) 05/07/2007

c) 11/07/2007

Donde: tinco = temperatura de entrada de los colectores.
toutco = temperatura de salida de los colectores.
tave = temperatura media de los colectores.
tmedi = temperatura del aire exterior.
ttan = temperatura del tanque.
En las Fig 4.2, 4.3 y 4.4,basándose en las gráficas experimentales, se

observa un descenso de temperaturas inicial, hasta situarse en un valor
aproximado de 18ºC, debido a que desde 0:00 hasta aproximadamente las
7:00, es el intervalo más frío de las noches de verano. Posterior a este
momento, y gracias a la radiación solar, tinco aumenta hasta llegar a un máximo
de 70ºC en torno a las 18:00, a partir del cual vuelve a bajar por la disminución
de la radiación solar hasta un valor de 35ºC.
Por otra parte, toutco es baja al inicio del día, causado porque el colector
plano no calienta el agua que lo atraviesa. Crece paralelamente a tinco hasta
que a partir de las 8:00, el colector recibe la suficiente radiación solar para que
el agua salga a mayor temperatura, alcanzando su máximo en 80ºC a las 16:00
aproximadamente. A partir de este instante comienza a decaer debido a la
disminución de la radiación solar.
Comparando con los resultados de la simulación, se demuestra que el

modelo matemático comete un error del 10% en la predicción de los valores de
tinco, toutco y tave, debido a que la carga térmica en la simulación es para una
vivienda mientras que experimentalmente se alimenta al generador de la
máquina , sin producir frío. Aún así, los tres parámetros poseen la misma
tendencia que los experimentales, por lo que se puede concluir que el modelo
posee un ajuste bueno con lo que ocurre en la realidad. Destacar que tinco,
toutco y tave poseen un mismo comportamiento hasta el comienzo de
funcionamiento debido a que todo el sistema se encuentra expuesto a la misma
temperatura ambiente.

En la Fig 4.2 , hay un descenso ttan porque previamente se calienta la

disolución hasta llegar a 65ºC . Se observa que a partir de las 14:00, ttan
alcanza una temperatura mayor de 65ºC, por lo que a partir de este instante se
produciría frío.
4.2.2 Rendimiento de la instalación solar (ηs).
El rendimiento instantáneo y total de la instalación se define como:
Qs
ηs = ec [4.65]
Q pl
Q sT
η sT = ec [4.66]
Q plT
En la Fig 4.5, 4.6, 4.7, se representa de manera gráfica el rendimiento de

la instalación a lo largo del día:
Figura 4.5. Evolución temporal del ηs el día 30/06/2007.




El rendimiento experimental es el rendimiento de los colectores (ηc),

calculado en el capítulo 3 pág 90-95, esta aproximación se ha realizado por la
imposibilidad de cuantificar las pérdidas térmicas de la instalación y obtener el
calor cedido de los colectores de manera experimental, al no poseer los datos
de caudal y temperaturas de entrada y salida del tanque. Las gráficas se
encuentran comentadas en el capítulo 3, pág 96.
En la Fig 4.5, 4.6, 4.7, basándose en las gráficas de la simulación, se

observa que a medida que transcurre el día el rendimiento de la instalación
aumenta hasta llegar a un máximo de 0,6 en las horas donde la radiación solar
es máxima, entre 12:00 y 14:00. A partir de este momento comienza a decaer
por la reducción de la intensidad solar y a las 18:00, la instalación solar ya no
produce calor suficiente por lo que deja de funcionar dando lugar a un
rendimiento nulo.
Comparando con las experimentales, el rendimiento de la instalación

simulado posee la misma tendencia que éstos, es decir aumenta a medida que
la radiación solar crece y disminuye cuando decrece la intensidad solar.
Se recogen en la Tabla 4.2 el rendimiento solar total diario de la

instalación:
Tabla 4.2. Valores de ηsT.

DÍA ηsT ηc
30/06/2007 0,275 0,432
05/07/2007 0,248 0,387
11/07/2007 0,332 0,373
PROMEDIO 0,285 0,397
Se demuestra que el rendimiento de la instalación se encuentra alrededor

del 29%, es menor que el experimental debido a que en éste no se han
contabilizado las pérdidas térmicas.

4.2.3 Calor recibido por los colectores (Qpl) y las pérdidas térmicas del la
instalación (Pt1).
Se muestran los resultados de la simulación en el cálculo de Qpl y Pt1, con

el objetivo de caracterizar cualitativamente en qué momentos se producen más
pérdidas y determinar cuáles son los intervalos críticos de funcionamiento, es
decir aquellos en los que la instalación no produciría calor para que funcionara
la máquina de absorción.
Pt1 posee la siguiente ecuación:
ec [4.67]
i i i (
Pt1 := Pco + Pcov + Pbox + Pins + HCEins + HCEco + HCEinss
i i i i i )
Donde: Pco= pérdidas dadas por la curva de normalización del colector,

incluyendo los efectos ópticos y de transferencia de calor.
Pcov = pérdidas de la cubierta.
Pbox = pérdidas en la caja de aluminio.
Pins = pérdidas disipada por los tubos del primario.
HCEco = calor acumulado por los elementos del colector.
HCEins = calor acumulado por los elementos del primario.
HCEinss = calor acumulado por los elementos del secundario.

En la Fig 4.8, 4.9, 4.10, se representa la evolución temporal de estos dos

parámetros:
Figura 4.8. Evolución temporal Qpl y Pt1 del día 30/06/2007.

Ambos parámetros crecen al aumentar la radiación solar hasta

aproximadamente las 10:00, donde Pt1 alcanzan valores próximos a 600 W/m2.
Posterior a este punto, Pt1 se mantiene por debajo de Qpl, esta separación
indica la producción de calor (Qs), el cual es transferido al tanque. A partir de
las 18:00 disminuyen estos valores debido a la disminución de la radiación
solar y de la temperatura ambiente.
Cabe destacar en el intervalo 12:00-16:00, donde Qpl alcanza sus valores

máximos y las pérdidas térmicas son bajas, por lo que la producción de Qs es
máxima en este período.

4.3 Influencia de la superficie solar y del volumen del tanque sobre la

instalación.
Se procede a analizar cómo afecta el área de colectores (Aco) y la

capacidad del tanque a la instalación. Para ello, se estudiarán dos casos:
1. Variación de la superficie solar, manteniendo constante el volumen del

tanque.
2. Variación del volumen del tanque manteniendo constante la superficie
solar.
Las variables a tener en cuenta para detectar esta influencia son: QsT, ηsT,
Qpl y ttan. Se realizan en un día típico, en este caso con las condiciones del día
05/07/2007. No se modifican las características geométricas de los elementos
de la instalación ni los materiales que la forman. Se tiene como condición de
diseño que : 65º C ≤ t tan ≤ 85º C , ya que este rango de temperaturas permite el
funcionamiento de la máquina de absorción.

1) Variación de la superficie solar y volumen del tanque constante.
Aco se modifica según el número de colectores (Ncol) y el volumen del

tanque se mantiene constante Vtan=2m3. En la Tabla 4.3, se recogen los
resultados de este caso:
Tabla 4.3. Resultados de variación de la superficie solar.
Ncol 5 10 15 20 25
ηsT 0,270 0,249 0,228 0,210 0,193
QsT [MJ] 69,770 129,008 177,435 216,939 249,614
QsT/Aco
[MJ/m2]
6,874 6,355 5,827 5,343 4,918
Qpl [MJ] 258,700 517,399 776,099 1034,798 1293,495

ttan [ºC] 61,59 67,398 72,91 77,50 81,44
A medida que aumenta el número de colectores, el rendimiento total de la

instalación baja. Esto es debido a que la relación QsT/Aco disminuye, la causa
de este fenómeno es el aumento de pérdidas en los colectores solares. Todas
las propuestas cumplen con la condición de diseño, aunque el óptimo para las
condiciones analizadas se encuentra con cinco colectores solares ya que
posee un rendimiento de 0,270, pero no se puede asegurar que satisfaga la
demanda frigorífica real, ya que Qu se ha simulado.
Cabe destacar que ttan aumenta a medida que la superficie solar se

incrementa. Atendiendo a esto último el máximo número de colectores es 25,
para un volumen de tanque de 2000l, ya que si se supera la condición de
diseño se supera.

2) Variación del volumen del tanque y superficie solar constante.
Se modifica Vtan y se mantiene constante el número de colectores Ncol=10.

En la Tabla 4.4, se recogen los resultados de este caso:
Tabla 4.4. Resultados de la variación del volumen del tanque.
Vtan [m3] 0,5 1 2 4 6
ηsT 0,198 0,227 0,249 0,264 0,260
QsT [MJ] 102,720 117,372 129,008 136,689 139,685
QsT/Aco
5,060 5,782 6,355 6,733 6,877
[MJ/m2]
Qpl [MJ] 517,399 517,399 517,399 517,399 517,399
ttan [ºC] 79,42 72,60 67,52 64,13 62,84
Se observa un incremento del rendimiento a medida que Vtan es mayor. La

causa de este efecto es el crecimiento de Qs/Aco, ya que se almacena más
energía térmica en el tanque. El valor de ttan disminuye debido a un aumento de
la cantidad de agua en el tanque.
Se destaca que el volumen máximo para 10 colectores es de 3 m3, ya que

la temperatura del tanque es el mínimo exigido 65ºC, este cálculo se ha
realizado en la simulación.

5. ESTUDIO ECONÓMICO

5. Estudio económico.
La instalación completa se compone de: una instalación solar térmica y

una máquina de absorción, además de los fan coils, el sistema de tuberías y
los equipos de regulación, adquisición de datos y bombeo.
En la realización de este presupuesto se ha considerado:
- Al ser una instalación experimental se encuentra

sobredimensionada.
- La longitud de las tuberías es mayor de lo convencional debido a
los distintos dispositivos de medida a lo largo de su recorrido y
por la gran altura del edificio Betancourt.
- Se incluye una instalación de gas natural, para que la máquina de
absorción siga funcionando en el caso de no existir radiación
solar suficiente para satisfacer la demanda energética [4].
Debido a esto, este estudio económico no representa los gastos reales de
una instalación estándar para un edifico comercial o familiar, o para una casa
unifamiliar.
En la Tabla 5.1 se encuentran representado el precio de los elementos de

la instalación completa:
Tabla 5.1. Costes de la instalación.

Concepto Precio [p.u] Nº unidades Coste total [€]
[26]
Colectores 915 €/panel 10 9150,00
[27]
Tubo 26 mm 6,05 €/m 62,22 376,43
[27]
Tubo 35 mm 9,19 €/m 86,78 797,51
[28]
Armaflex 26 mm 0,81 62,22 50,40
[28]
Armaflex 35 mm 0,84 86,78 72,90
[29]
Intercambiador Calor 250 1 250,00
[29]
Tanque acero (2000 l) 2000 1 2000,00
[29]
Vaso de expansión 100 1 100,00
Fan coil [29] 200 1 200,00
[27]
Codos 0,44 €/codo 111 48,84
[30]
Válvulas (todo-nada) 2,9 €/ud 22 63,80
[30]
Válvulas (regulación) 271,58 €/ud 4 1086,32
[30] [31]
Aguaglicol 18 €/ud 2,73 49,14

Continuación Tabla 5.1. Costes de la instalación

Concepto Precio [p.u] Nº unidades Coste total [€]
2 [29]
Bombas 1446,3 ----- 1446,30
Caldera Gas Natural3 [32] [37] 1162,25 1 1162,25
Regulador Gas Natural [29] 13,71 1 13,71
Llave Gas Natural [29] 5,02 1 5,02
Máquina absorción4 18000 1 18000,00
Subtotal 34872,61 €
Mano de obra5 10461,78 ---- 10461,78
Ayudas estatales6 - 10461,78 ---- -10461,78
Ayudas autonómicas7 -225 €/m2 20 -4500
Costes fijos electricidad8 67,41 ---- 67,41
Costes fijos Gas Natural9 51,24 ---- 51,24
Total 30491,26 €
1 Para el cálculo de los litros de etilenglicoglicol se suponen las tuberías del

circuito primario y los colectores solares llenos (1,3 l/panel) además de un 5%
más en el vaso de expansión, de lo cual el 20% es Etilenglicol y lo restante
agua que se desprecia su coste:
⎡⎡ l ⎤ ⎡ l ⎤ ⎤
⎢ ⎢10 paneles ⋅ 1,3 + 10 paneles ⋅ 1,3 ⋅ 0, 05 ⎥ ⋅ 0,2 = 2,73 l
⎢⎣ ⎣ panel ⎥⎦ ⎢
⎣ panel ⎥⎦ ⎥⎦
2 El precio de las bombas se encuentran en la tabla 5.2:
Tabla 5.2. Precio de las bombas

Modelo Nº unidades Coste unitario [€]
Grundfos UPS 25-60 2 216,06
3 Caldera mural VAILLANT atmoTEC pro 20640, de 20640 kcal/h, con un 90%
[37]
de rendimiento, potencia 24 kW y temperatura máxima de 130 ºC . Se ha
elegido esta caldera porque cumple con las especificaciones de la instalación y
posee un precio medio.

4 El precio de la máquina de absorción ha sido facilitado por Doña María del

Carmen Venegas Bernal, profesora del departamento de Ingeniería Térmica y
de Fluidos de la Universidad Carlos III.
5 La mano de obra de una instalación se presupuesta como el 30-35% del

coste total, en este caso se supone en un 30%.
6 Según IDAE, al ser una instalación solar térmica de 20 m2 le corresponde una

subvención del 30% del coste total [33].
7 Según la Comunidad de Madrid, la subvención para colectores solares

térmicas equivale a 225 €/m2 [33].
8 En la Tabla 5.3 se recoge el consumo eléctrico de aquellos elementos de la

instalación que lo necesitan:
Tabla 5.3. Costes de electricidad.

Consumo
Elemento Consumo [kWh] Consumo [kWh/mes]
[kWh/día]
Bomba 1º 0,32 7,68 230,4
Bomba 2º 0,145 3,48 104,4
(*)
Bomba 3º 0,1 2,4 72
Máquina de abs 0,21 5,04 151,2
Fan coil 0,1 2,4 72
Total 24,6 738
(*) Es el consumo de ambas bombas.
Estos elementos funcionan todo el día. Aplicando la tarifa 3.0.1 del BOE 156
28/06/2008 y despreciando la dada de alta a la red, se obtiene un coste
eléctrico de:
kWh € mes
630 ⋅ 0,107 = 67,41 €
mes kWh

9 Sin contar los gastos de dada de alta en la red de Gas Natural, el coste de
gas natural asciende, aplicando la tarifa T3 [34]:
⎡ € ⎤
⎢0,15 kWh ⋅ ⋅0,037 kWh ⎥ + 44,17 + 16% IVA = 51,24 €
⎣ ⎦
En la Tabla 5.4 se representa el gasto de explotación de la instalación

anualmente:
Tabla 5.4. Gastos de explotación

Concepto Precio [€/año]
Mantenimiento de la instalación solar (*) 100
Electricidad (**) 947,52
Gas natural (**) 731,01
Total 1778,53
(*) Consultado con un técnico del mantenimiento de la instalación.
(**) Resulta de multiplicar por doce meses los 78,96 €/mes de electricidad.
[35]
(***) Hogares con calefacción de gas natural: 116,13€/año la revisión y los
doce meses de 51,24 €/mes de gasto de Gas Natural.
En resumen, se muestra en la Tabla 5.5 un desglose del presupuesto de

toda la instalación, en función de los subsistemas que la componen y en la
Figura 5.1 se ha realizado una representación gráfica de la misma:
Tabla 5.5 Desglose del presupuesto de la instalación.

Precio [€] Porcentaje (%)
Instalación solar 15770,59 44,80
Caldera Gas 1232,22 3,50
Máquina de absorción 18000,00 51,13
Fan coil 200,00 0,57

Figura 5.1. Desglose del presupuesto de la instalación.
Tras haber realizado el estudio económico se puede concluir:
- El coste total de la instalación es de 30502,81 € debido a que esta

instalación es experimental. Es un precio elevado para que se
pueda adquirir por una familia no así por una industria,
- Las ayudas recibidas del Estado y de la Comunidad de Madrid
disminuyen significativamente el coste total ahorrándose casi
5000 €, aún así las subvenciones son bajas.
- Los gastos de explotación es elevado, sobre todo en los aspectos
de electricidad y gas por lo que se recomienda un uso más
racional de la instalación.
- La mitad del presupuesto es destinado a la máquina de absorción,
esto es debido a que es una tecnología empleada en grandes y
medianas potencias, resultando competitiva con los sistemas de
compresión mecánica al utilizar calores residuales además se
encuentra en investigación. por lo que su adquisición requiere
una inversión alta.

6. ESTUDIO MEDIOAMBIENTAL

6. Estudio ambiental.
Se denomina Evaluación de Impacto Ambiental (EIA) a todo el

procedimiento administrativo cuyo fin es establecer si un determinado proyecto
es aceptable o no desde un punto de vista ambiental. Con esta herramienta de
trabajo se pretende, de una forma metodológica, la descripción de distintas
alternativas a un proyecto, la identificación de las distintas acciones que posee
éste en el entorno natural así como la valoración de los diversos elementos
ambientales. Realizada dicha valoración, se procede a la identificación y
cuantificación de los diversos impactos ambientales recomendando medidas
protectoras y correctoras de los mismos [23].
Dentro de los diversos documentos de la EIA se encuentra el denominado

Estudio de Impacto Ambiental (EstIA) que es el documento técnico exigido por
la Administración en el Reglamento de Evaluación de Impacto Ambiental (Real
Decreto 1131/1988), con el fin de conocer de forma detallada cómo afectaría
un proyecto dado al entorno [23].
De los múltiples apartados que contiene el EstIA, para la realización del

presente proyecto, se va a centrar en los aspectos más relevantes del
inventario ambiental procediendo a una valoración cualitativa de los impactos
ambientales que se puedan encontrar, tanto en la instalación solar como en la
máquina de absorción, y recomendando algunas medidas para la minimización
de los mismos.
Los aspectos del inventario ambiental que se van a considerar para el

análisis de impactos van a ser:
- Aire.
- Agua.
- Ruido.

6.1. Análisis de impactos sobre el aire.
En la Tabla 6.1 se detallan los impactos y las posibles medidas

correctoras de la instalación:
Tabla 6.1. Impactos sobre el aire.

Fase Impacto Medidas
Generación de Realización
partículas y sólo de los
polvo en el orificios
anclaje de los indispensable
Instalación Montaje paneles. para la
solar instalación de
los soportes
AIRE
de los
paneles.
Explotación Ninguno Ninguno
Montaje Ninguno Ninguno
Transporte de Colocación
legionella de filtros
Máquina de debido a la esterilizantes
absorción Explotación emisión de en la salida
gotas de agua de la torre.
por la torre de
refrigeración.
La instalación solar térmica en la fase de explotación no contamina la

atmósfera debido a que la fuente de calor empleada para la elevación de la
temperatura del agua es la energía térmica proveniente de la radiación solar,
por lo que no se genera ninguna emisión de gases de efecto invernadero (CO,
CO2), ni de agentes que producen la lluvia ácida (SOx y NOx).

La máquina de absorción no ejerce ningún impacto en la fase de

construcción ni en la de explotación debido a que el refrigerante empleado es el
[1]
agua cual posee un potencial destructor del ozono nulo .El absorbedor
empleado es LiBr, el cual se encuentra en ambas fases en estado líquido por lo
que no existirá fuga en estado gaseoso durante el funcionamiento normal del
aparato.
6.2 Análisis de impacto sobre el agua.
En la Tabla 6.2 se detallan los impactos y las posibles medidas

correctoras de la instalación:
Tabla 6.2. Impactos sobre el agua.

Fase Impacto Medidas
Instalación Montaje Ninguno Ninguno
solar Explotación Ninguno Ninguno
Montaje Ninguno Ninguno
AGUA
Control de la
corrosión y
Contaminación fisuras de los
del agua de conductos del
Máquina de
refrigeración absorbedor,
absorción Explotación
por LiBr. para que no
existan
roturas de los
mismos.
La instalación solar térmica posee una disolución de 80% agua y un 20%

etilenglicol (anticongelante) la cual se emplea como fluido caloportador, al
encontrarse en un circuito cerrado sin comunicación con la red de tuberías no
existe riesgo de contaminación de aguas del edificio.

6.3 Ruido.
Los elementos de la instalación que generan ruido durante su

funcionamiento son: la máquina de absorción, el calentador de gas y los fan
coils.
El calentador de gas se encuentran en la sala de calderas del edificio

Betancourt, es decir en la parte inferior del edificio, la cual se encuentra alejada
de las zonas de trabajo del mismo por lo que no requiere un aislamiento
acústico que insonorice esta zona pero si uno que reduzca los sonidos emitidos
por estos equipos.
Los fan coils y la máquina de absorción, se encuentran situados en la

parte superior de la nave de Ingeniería Térmica, se encuentran en una zona
habilitada para el trabajo por lo que el nivel de ruido exterior no podrá ser
superior a 60 dB(A) e interior a 40 dB (A) según el Decreto 78/1999, por el que
se regula el régimen de protección contra la contaminación acústica de la
Comunidad de Madrid. Aplicando lo anteriormente descrito, los fan coils se
deben aislar para reducir su emisión de ruido si supera estos límites.
6.4 Cantidad no emitida de CO2.
Se muestra a continuación la cantidad promedio no emitida de CO2, al

emplear una instalación de frío solar frente al sistema de compresión mecánica,
para una casa de 90 m2:
Tabla 6.3. Emisiones de CO2.

Frío Solar Compresión mecánica
30/06 05/07 11/07 30/06 05/07 11/07
QuT (kWh)(*) 12,34 36,15 21,17 12,34 36,15 21,17
Kg CO2 (**) 2,71 7,95 4,65 6,41 18,80 11,00
Promedio
5,10 12,07
(kg CO2)
(*) Valores obtenidos en la simulación.
(**) Se producen 0,22 kg CO2 en la instalación de frío solar por el gas natural[38]
Se produce 0,52 kg CO2/kWh en la producción de electricidad[38].

Se observa una reducción del 50% de emisión, por lo que se refleja de

forma cuantitativa, que las instalaciones de frío solar generan menos CO2, que
los ciclos de compresión mecánica. Desde un punto medioambiental, este tipo
de sistemas son más limpios que los eléctricos debido al empleo de la
radiación solar como fuente de energía.

7. CONCLUSIONES

7. Conclusiones.
Como resultado del estudio de este proyecto, se ha llegado a las

siguientes conclusiones.
1. Según la bibliografía consultada, las instalaciones de frío solar

contribuyen a satisfacer la demanda frigorífica durante los meses de
verano. Es una tecnología en auge y de uso entre las grandes y
medianas empresas, aunque se está trazando una línea para las
viviendas unifamiliares.
2. La instalación analizada recibe un promedio diario de 49204,170

W/m2, obteniéndose un calor medio procedente de los colectores de
235131,021 KJ con un rendimiento promedio de colectores del
39,786%.El rendimiento de los colectores alcanza su máximo
aproximadamente entre las 10:50-16:45, con valores situados entre el
65-55% coincidiendo con la radiación solar máxima (960W/m2). El
intervalo de temperaturas de salida del agua se encuentra entre 65-
80ºC, apto para el funcionamiento de la máquina de absorción.
3. Las temperaturas obtenidas en la simulación poseen la misma

tendencia que los experimentales y corroboran el intervalo
experimental de trabajo de los colectores 65ºC-80ºC, aunque con la
simulación se llegan a alcanzar los 85ºC.
4. La instalación solar posee un rendimiento medio del 28,5%, llegando a

alcanzar un máximo del 60% en su período de funcionamiento. Este
parámetro disminuye al aumentar la superficie solar y se incrementa al
aumentar el volumen del tanque.
5. La instalación posee un coste total de 30491,26€ y un gasto de

explotación de 1778,53€/año, lo que requiere una inversión inicial
elevada que puede ser sufragada por grandes y medianas empresas.

6. La instalación de frío solar posee prácticamente un impacto

medioambiental nulo, no emite gases de efecto invernadero y protege
la capa de ozono.
7. En una vivienda de 90 m2, se reduciría las emisiones de CO2 en un

50% respecto a los sistemas de compresión mecánica y supondría
una ahorro energético de 75 €/mes. Disminuyendo de esta formal
dependencia de la energía eléctrica en la casa unifamiliar.
Trabajos futuros
Se afirma que los resultados de la simulación realizan una primera

aproximación válida, aunque es necesario que el modelo sea más realista
incluyendo los circuitos de vuelta de cada uno de los circuitos de la instalación.
Se recomienda la realización de simulaciones, partiendo de la base de

este proyecto, con una variable que relacione Vtan/Ncol de forma que se puedan
extrapolar los resultados a otro tipo de instalaciones. Además de proyectos que
caractericen de forma experimental la máquina de absorción para determinar
su comportamiento y así alcanzar una comparación entre ésta y las
simulaciones que se realicen.

8. BIBLIOGRAFÍA

8. Bibliografía
[1] Izquierdo Millán, Marcelo. “Sistemas de ciclos de absorción para la

climatización y refrigeración”.Curso de Postgrado (1998), Escuela Politécnica
Superior de la Universidad Carlos III de Madrid.
[2] A. Pongtornkulpanicha, S. Thepa, M. Amornkitbamrung, C. Butcher.

“Experience with fully operational solar-driven 10-ton LiBr/H2O single-effect
absorption cooling system in Thailand”. Renewable Energy 33 (2008) 943-949.
[3] Constantinos A. Balarasa, Gershon Grossmanb, Hans-Martin Henningc,

Carlos A. Infante Ferreirad, Erich Podessere, Lei Wangd y Edo Wiemkenc.
“Solar air conditioning in Europe—an overview”. Renewable and Sustainable
Energy Reviews 11 (2007) 299–314.
[4] Darine Zambrano, Carlos Bordons, Winston Garcia-Gabin, Eduardo F.

Camacho. “Model development and validation of a solar cooling plant”. I n t e
rna t i o n a l Journa l O f R e f r i g e r a t i o n xxx (2007) 1 – 1 3.
[5] García Casals, Xavier. “Solar absorption cooling in Spain: Perspectives and
outcomes from the simulation of recent installations”. Renewable Energy 31
(2006) 1371–1389.
[6] Hans-Martin Henninga. “Solar assisted air conditioning of buildings – an

overview”. Applied Thermal Engineering 27 (2007) 1734–1749.
[7] N.E.Wijeysundera. “Simplified models for solar-powered absorption cooing

systems”. Renewable Energy 16 (1999) 679-684.
[8]http://es.wikipedia.org/wiki/Absorci%C3%B3n_(qu%C3%ADmica),26/02/2008
[9] D.S Kim y C.A. Infante Ferreira “Air-cooled Solar Absorption Air
Condiotioning, Final Report”.Delft University of Technology, 2005.
[10] Altener Project Number 4.1030/Z/02-121/2002. “Promoting Solar Air

Conditioning: Technical overview of active techniques”.Año 2002.

[11] Edo Wiemken “Experiences from design studies, simulation studies and
installations” International Energy Agency Solar Heating and Cooling
Programme.Task 38: Solar Air-Conditioning and Refrigeration. Año 2007.
[12] John A. Duffie y William A. Beckman “Solar Engineering of thermal

processes”. Editorial John Wiley & Sons, Inc. Año 1991.
[13]http://www.quimica.urv.es/~w3siiq/DALUMNES/99/siiq51/Colc.html,2/03/20
08.
[14]http://www.sitiosolar.com/colectores%20de%20tubo%20de%20vacio.htm,
2/03/2008.
[15]http://www.psa.es/webesp/areas/cilindroparabolica/general.html, 2/03/2008.
[16]http://personales.ya.com/universal/TermoWeb/EnergiasAlternativas/solar/P
DFs/Solar4.pdf, 2/03/2008.
[17] http://www.natureduca.com/energ_alternat_colectsolar2.php, 2/03/2008.
[18] Izquierdo Millán, Marcelo. “Energías Renovables, Módulo: Energía Solar”.

Departamento de Ingeniería Térmica de la Escuela Politécnica Superior de la
Universidad Carlos III de Madrid.
[19].http://www.alfalaval.com/solution-finder/products/gasketed-industrial-range-
phe/Documents/M6.pdf, 15/07/2008.
[20]http://webserver.dmt.upm.es/~isidoro/bk3/Appendices/Propiedades%20term
icas%20de%20la%20materia.pdf, 15/07/2008.
[21] http://pc-134-222.uc3m.es/meteo.htm, 15/07/2008.
[22] Catálogo de colectores THERMOSOLAR.
[23] Garmendia Salvador A., Salvador Alcaide, A., Crespo Sánchez, C.,
GarmediaSalvadorL., “Evaluación de Impacto Ambiental” Editorial Pearson-
Prentice Hall, Año 2005.

[24]http://www.goodfellow.com/csp/active/STATIC/S/Nitruro_de_Silicio-
Nitruro_de_Aluminio-Oxido_de_Aluminio.HTML, 22/08/2008.
[25] http://enciclopedia.us.es/index.php/Kript%C3%B3n, 22/08/2208.
[26] Tarifa Thermosolar 2006.
[27] http://www.preoc.es/precios/U/U28/U28AF.asp, 27/08/2008.
[28] http://www.preoc.es/precios/U/U15/U15AM.asp, 27/08/2008.
[29] Alcañiz Martín, Sergio, “Modelado y estudio económico y medioambiental

de la instalación solar térmica de la UC3M con aplicación a calefacción”,
Proyecto Fin de Carrera, Ingeniería Industrial, Año 2005.
[30] http://www.costonet.com.mx, 27/08/2007.
[31] http://www.quimicosballester.com.ar/, 27/08/2008.
[32] http://www.preoc.es/precios/D/D30/D30OD.asp, 27/08/2008.
[33] http://www.soliclima.com/pt/subvencoes.html, 27/08/2008.
[34]http://www.hc.es/index.asp?MP=4&MS=50&MT=65#tarifas&resolucion=128
0*768, 27/08/2008.
[35]http://www.gasnatural.com/servlet/ContentServer?gnpage=1-1-
1&centralassetname=1-1-3-1-1-0-0, 27/08/2008.
[36]http://www.idae.es/uploads/documentos/documentos_2.1.7_Energia_Solar_
Termica_7b7f5630.pdf, 4/09/2008.
[37] http://www.calderasmadrid.com/calderas/turbotecplus.html, 6/09/2008.
[38] M.C. Rodríguez, P. Aumente, M. Izquierdo Millán, A. Lecuona Neumann,

R. Salgado Mangual “Energy and carbon emission savings in Spanish housing
air-conditioning using solar driven absorption system” Applied Thermal
Engineering 28 (2008) 1734-1744, Año 2008.

[39] http://termica.uc3m.es/Solar/yazaki.jpg, 9/09/2008.
[40] http://termica.uc3m.es/Solar/vacio.jpg, 9/09/2008.
[41] Incropera, Frank P. “Fundamentos de transferencia de calor”, Editoria

Pearson, Año 1999.

PFC Victor Salguero Fernandez

Cargado por

Copyright:

Formatos disponibles

PFC Victor Salguero Fernandez

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

PFC Victor Salguero Fernandez

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD CARLOS III DE MADRID

ESCUELA POLITÉCNICA SUPERIOR

DEPARTAMENTO DE INGENIERÍA TÉRMICA Y DE FLUIDOS

Valoración de colectores solares térmicos

PROYECTO FIN DE CARRERA

INGENIERÍA TÉCNICA INDUSTRIAL MECÁNICA

AUTOR: VÍCTOR SALGUERO FERNÁNDEZ

TUTORA: Dª MARÍA DEL CARMEN VENEGAS BERNAL

LEGANÉS, SEPTIEMBRE, 2008

En primer lugar a María Venegas Bernal, gracias por haberme introducido en la

A mis amigas Yolanda, Rebe, Laura (pimaaa), me habéis acompañado toda mi

A mi tridente: María, Pablo (Melenudo) y Mari Carmen, personas en mi vida de

I Importancia del proyecto…………………………………………………………1

III Descripción por capítulos………………………………………………………6

1.1 Fundamentos termodinámicos de la máquina de absorción……………9

1.2 Componentes de la máquina ideal y su COP……………………………..13

1.3 Fluidos de trabajo………………………………………………………………18

1.4 Tipos de máquina de absorción……………………………………………..20

1.5 Consideraciones ambientales de la máquina de absorción................23

1.6 Descripción de una planta solar…………………………………………….24

1.7 Actualidad de la planta solar y la máquina de absorción………………25

2.1 Colectores solares planos……………………………………………………35

2.1.1 Descripción y funcionamiento de un colector plano…………36

2.1.2 Colector solar plano de vacío……………………………............40

2.1.3 Ecuaciones energéticas de un colector solar plano………….41

1)Rendimiento del colector solar plano……………………….41

2)Ecuación fundamental de colector…………………………..42

3)Pérdidas a través de la cubierta……………………………...43

4)Pérdidas a través de la parte posterior y laterales………..48

“Departamento de Ingeniería Térmica y Fluidos” pág:i

6)Calor útil transferido por el colector………………………...50

7)Rendimiento de colectores en régimen

8) Pérdidas de calor por transferencia de calor de la

9)Calor absorbido por el colector solar……………………….55

10)Calor absorbido por los componentes de la

11) Energía solar disponible…………………………………….56

2.2 Colectores concentradores…………………………………………………..58

2.2.1 Descripción y funcionamiento de un colector

2.2.2 Tipos de colectores solares concentradores…………............60

2.2.2.1 Tubos evacuadores………………………………………60

2.2.2.2 Colectores parabólicos………………………………….66

3.1 Descripción de la instalación………………………………………………68

3.1.1 Colectores solares planos de vacío……………………………...70

3.1.2 Intercambiador de calor……………………………………………72

3.1.3 Máquina de absorción LiBr-H2O…………………………............73

3.1.4 Tanque de almacenamiento……………………………………….76

3.1.5 Circuitos hidráulicos de la instalación…………………………..77

“Departamento de Ingeniería Térmica y Fluidos” pág:ii

3.1.7 Fluidos de trabajo…………………………………………………...81

3.1.8 Sistema de regulación, bombeo y adquisición de datos…….82

3.2 Instalación solar térmica…………………………………………………….85

3.3 Máquina de absorción……………………………………………………….99

4.1 Desarrollo del modelo………………………………………………………..100

4.1.2 Caract Col1………………………………………………………….105

4.2Resultados del modelo y comparación con los datos

4.2.2 Rendimiento de la instalación…………………………………...125

4.2.3 Calor recibido por los colectores y pérdidas térmicas de la

4.3 Influencia de la superficie solar y del volumen del tanque de sobre la

1) Variación de la superficie solar y volumen del tanque