PFC - Marta Torralba Rubio
PFC - Marta Torralba Rubio
PFC - Marta Torralba Rubio
de aislamiento en tuberías
de sales fundidas
Para una Central de Captadores
Cilindro-Parabólicos (CCP)
Año 2014
Optimización del espesor de aislamiento
2
Optimización del espesor de aislamiento
UNIVERSIDAD DE SEVILLA
3
Optimización del espesor de aislamiento
RESUMEN
El tipo de fluido de trabajo empleado en este caso es una sal fundida compuesta por
un 40% de KNO3 y 60% de NaNO3, que permite alcanzar una temperatura de
operación de 550ºC. La sal fundida es un fluido de trabajo diferente de los comúnmente
empleados en centrales solares con tecnología cilindro-parabólica, debido al rango de
temperaturas de operación, porque presenta un punto de solidificación elevado, en torno
a los 200 ºC, que condiciona las pérdidas térmicas a lo largo de la extensa red de
tuberías necesaria en este tipo centrales.
La central bajo estudio en este proyecto es la primera que emplea la sal fundida
como fluido de trabajo y de almacenamiento con un sistema de captadores cilindro-
parabólicos. Fue inaugurada en julio del 2010, llamada Archimede, sitiada en Priolo
Gargallo, Sicilia (Italia). Está formada por 9 lazos en paralelo, de 6 captadores
constituidos por 8 módulos. La potencia eléctrica bruta de la central es de 5 MW,
resultando una potencia neta de 4,72 MW. El rendimiento (global) solar a eléctrico es
del 15,1 %.
5
Optimización del espesor de aislamiento
ABSTRACT
The thermal losses in a concentrating solar power plant (CSP) based on the
parabolic trough technology are analyzed and a method for the optimization of the
insulation thickness in the distribution and collection pipeline network is proposed and
implemented. The plants of this type are intended for converting the direct component
of solar radiation into electrical energy, by using parabolic trough collectors (PT) to
concentrate the solar energy onto a receiver tube placed in its focal line, and by
circulating a working fluid that increases its temperature as it flows through the
receiver.
The working fluid used in this case is a molten salt composed of 40% of KNO3 and
60% of NaNO3, which allows reaching an operating temperature of 550 °C. Compared
to the thermal oil currently used in PT plants, this molten salt mixture permits operating
at higher temperature, while its high freezing point (230 ºC approx.) makes it necessary
to make special provisions to avoid freezing at any point of the circuit (receivers,
pipelines, valves, pumps, etc.) and requires better insulation to reduce thermal losses.
The analysis of the thermal losses in the network of pipelines requires understanding
the mechanisms of heat transfer and its particularization to the specific operating
conditions of PT plants. The insulation does not eliminate thermal losses, simply
reduces it. After estimating the heat loss results as a function of the insulation thickness,
the optimum thickness is derived by establishing an economic balance between the
insulation costs and the costs associated to the thermal losses. The optimal thickness is
that which minimizes the total cost
The methodology has been applied to the Archimede plant in Priolo Gargallo, Sicily
(Italy). This is the first commercial PT plant that uses molten salt as working fluid and
storage medium. It was opened in July 2010, and consists of 9 parallel loops of 6
collectors each. The gross electrical output of the plant is 5 MW, resulting in a net
power of 4.72 MW. The global solar to electrical efficiency is 15.1%.
6
Optimización del espesor de aislamiento
ÍNDICE
1. OJBETIVO ......................................................................................................................................... 18
2. INTRODUCCIÓN.............................................................................................................................. 19
2.1. DESCRIPCIÓN GENERAL DE LA TECNOLOGÍA .............................................................. 19
2.2. HISTORIA ................................................................................................................................ 19
2.3. CENTRAL “ARCHIMEDE” .................................................................................................... 20
2.4. PROGRAMA DE CÁCULO “EES” ......................................................................................... 22
3. CARACTERÍSTICAS DE LA INSTALACIÓN ................................................................................ 23
3.1. FLUIDO DE TRABAJO ........................................................................................................... 23
3.1.1. ACEITE TÉRMICO ............................................................................................................. 23
3.1.2. AGUA LIQUIDO/VAPOR ................................................................................................... 25
3.1.3. SALES FUNDIDAS ............................................................................................................. 26
3.1.4. GASES A PRESIÓN ............................................................................................................ 29
3.2. CONFIGURACIÓN DEL CAMPO SOLAR ............................................................................ 30
3.3. TUBERÍAS ACERO INOXIDABLE ....................................................................................... 31
3.4. AISLAMIENTOS ..................................................................................................................... 34
3.4.1. CONDICIONES DE DISEÑO ............................................................................................. 34
3.4.2. LANA DE ROCA ................................................................................................................. 37
3.4.3. LANA CERAMICA ............................................................................................................. 40
3.4.4. INSTALACIÓN.................................................................................................................... 41
3.5. CAPTADOR ............................................................................................................................. 46
4. DIMENSIONAMENTO DE LA RED DE TUBERÍAS ..................................................................... 47
4.1. DATOS METEOROLÓGICOS ................................................................................................ 47
4.2. MULTIPLO SOLAR................................................................................................................. 49
4.3. DIMENSIONES DEL CAMPO SOLAR .................................................................................. 51
4.4. FLUIDO DE TRABAJO - SAL FUNDIDA ............................................................................. 54
4.5. TEMPERATURAS INICIALES ............................................................................................... 55
4.6. VELOCIDAD DEL FLUIDO Y CAUDAL MÁSICO .............................................................. 58
4.7. PÉRDIDAS DE CARGA .......................................................................................................... 60
4.8. DIÁMETRO ............................................................................................................................. 63
5. CÁLCULOS TÉRMICOS .................................................................................................................. 65
5.1. ESQUEMA ............................................................................................................................... 65
5.2. RESISTENCIAS TÉRMICAS .................................................................................................. 66
5.2.1. CONVECCIÓN FORZADA EN LA PARTE INTERIOR DEL CONDUCTO .................... 66
5.2.2. CONDUCCIÓN EN LA TUBERÍA DE ACERO ................................................................. 68
5.2.3. CONDUCCIÓN EN EL AISLAMIENTO ............................................................................ 70
5.2.4. CONVECCIÓN NATURAL (LIBRE) CON EL AIRE EXTERIOR .................................... 71
5.2.5. RADIACIÓN ENTRE LA SUPERFICIE EXTERIOR Y EL CIELO .................................. 73
5.2.6. RADIO CRÍTICO ................................................................................................................ 77
5.3. PROCESO ITERATIVO DE CÁLCULO ................................................................................. 78
7
Optimización del espesor de aislamiento
8
Optimización del espesor de aislamiento
INDICE DE FIGURAS
9
Optimización del espesor de aislamiento
Figura 37. Fotografía campo solar Archimede. Imagen del uso de liras.
Figura 39. Imagen entrada y salida de los lazos de campo solar de Arhimede
10
Optimización del espesor de aislamiento
Figura 47. Esquema de la retribución necesaria para dar viabilidad a cada tipo de
tecnología renovable
Figura 50. Fotografía tubería del campo solar sin aislamiento y tubería con
aislamiento- Archimede
11
Optimización del espesor de aislamiento
INDICE DE TABLAS
Tabla 10. Temperatura de salida del campo solar (Tramo 5)– Red caliente
Tabla 15. Pérdidas térmicas totales del año modelo – Red fría
Tabla 17. Temperatura de salida del tramo 1 (Tramo más alejado) – Red fría
Tabla 21. Influencia del espesor de aislamiento sobre la temperatura de salida con
12
Optimización del espesor de aislamiento
Tabla 22. Influencia del espesor de aislamiento sobre la temperatura de salida con
distinta irradiación – Red Fría
13
Optimización del espesor de aislamiento
NOMENCLATURA
Superficie exterior
Factor de fricción
Aceleración de la gravedad
Radiación incidente
14
Optimización del espesor de aislamiento
Tasa de interes
Factor de pérdidas
Presión de entrada
Presión de salida
Pr Prandtl fluido
Prandtl aire
Potencia perdida
15
Optimización del espesor de aislamiento
Resistencia de radiación
Reynolds exterior
Reynolds interior
Múltiplo solar
Temperatura media
Pérdida de carga
16
Optimización del espesor de aislamiento
Constante de Stefan-Boltzmann
17
Optimización del espesor de aislamiento
1. OJBETIVO
18
Optimización del espesor de aislamiento
2. INTRODUCCIÓN
Por otra parte las plantas con almacenamiento pueden generar mayor electricidad
para la misma potencia instalada, haciendo por tanto un uso más eficiente de la turbina.
Por el contrario, las plantas con almacenamiento requieren una mayor inversión.
2.2. HISTORIA
19
Optimización del espesor de aislamiento
Desde el año 2001 se puso en marcha un proyecto por parte de ENEA para
investigar el empleo de sales fundidas como fluido de trabajo y almacenamiento en
Central Cilindro-Parabólicas (CCP) [12].
20
Optimización del espesor de aislamiento
Probar la eficiencia óptica y térmica de los nuevos captadores solar, que emplean
CERMET y permiten alcanzar hasta 550ºC. Llamados HCEMS 11 [7]
21
Optimización del espesor de aislamiento
Comenzó a construirse en Julio del 2008 y se inauguró en Julio del 2010, con una
potencia eléctrica bruta de 5 MW y una potencia eléctrica neta de 4,72 MW. El campo
solar consta de 9 lazos con 6 captadores de 8 módulos cada uno. Con un sistema de
almacenamiento con capacidad de 7 horas equivalentes.
22
Optimización del espesor de aislamiento
3. CARACTERÍSTICAS DE LA INSTALACIÓN
En la actualidad los fluidos de trabajo más empleados son los orgánico sintéticos o
aceites sintéticos, entre los que destaca la mezcla compuesta por un 26,5% de oxido de
difenilo (C12H10O) y un 73,5% de bifenilo (C12H10). Esta mezcla presenta el mejor
compromiso entre economía y comportamiento aunque posee grandes dificultades en su
manejo, que condicionan los resultados de explotación de las plantas.
El aceite trabaja entre los 290ºC a la entrada del lazo y los 390ºC a la salida, lo que
supone un incremento de temperatura de 100ºC en el interior de los lazos.
Dichos aceites presentan las siguientes desventajas para su uso como fluido de
trabajo en sistemas de captadores cilindro-parabólicos (CCP):
23
Optimización del espesor de aislamiento
Actualmente se analizan diversos fluidos como candidatos para sustituir los aceites
térmicos en las plantas termosolares, como son: Sales Fundidas; Agua (con su
correspondiente cambio de estado líquido/vapor) y gases a presión. La tabla adjunta
reúne las ventajas y desventajas de dichas alternativas frente al aceite térmico. [5]
24
Optimización del espesor de aislamiento
La situación más lógica e ideal resultaría ser emplear agua (liq/vapor) como fluido,
por su carácter medioambiental y sin riesgo de incendio o contaminación, a diferencia
del aceite térmico que si lo plantea en caso de fuga.[5]
25
Optimización del espesor de aislamiento
La alternativa más atractiva hasta el momento para sustituir al aceite térmico, usado
prioritariamente en la actualidad, son las sales fundidas. Dichas sales fundidas son
empleadas actualmente en los sistemas de almacenamiento.
La tabla refleja las características técnicas que permiten la comparación entre los
dos sistemas para una planta de la misma potencia. Se llega a la conclusión de que las
sales tienen mejores propiedades cuando el sistema es directo por varios factores, se
requiere menor cantidad de sales, se almacenan a mayor temperatura y se eliminan las
pérdidas en el intercambiador de aceite y sal fundida, puesto que no sería necesario en el
sistema directo de almacenamiento.
26
Optimización del espesor de aislamiento
Las principales ventajas de las sales fundidas como fluido de trabajo son [5]:
Posibilidad de incrementar la temperatura del vapor por encima de los 400ºC del aceite.
Las sales son estables hasta temperaturas próximas a los 600ºC. Como consecuencia de
elevar la temperatura de salida del lazo se mejora la eficiencia del ciclo de Rankine de la
turbina con el vapor que se introduce. La eficiencia del ciclo aumenta de un 37,6% que
se obtiene con la temperatura máxima del aceite, del orden de 393ºC a un rendimiento
del 40% con el empleo de las sales fundidas a 500ºC. [12]
No son contaminantes. De hecho, las sales fundidas se usan como fertilizante de uso
común en la agricultura. Ofrecen mayor resistencia a las fugas, y en caso de producirse
una fuga, esta se solidifica antes de filtrarse en el terreno debido al alto punto de
solidificación que presentan las sales, en torno a 200ºC.
Las desventajas que presentan dichas sales fundidas son las siguientes:
Presentan un mayor índice de corrosión que el aceite sintético, que afectará del mismo
modo a la selección de los componentes del sistema de trasporte de las sales fundidas.
27
Optimización del espesor de aislamiento
Sal Solar (Solar salt): Fue seleccionada como la sal más práctica para producción
de energía en con tecnología de receptor central de torre, puesto que el límite superior
de temperatura de funcionamiento 600ºC permite utilizar las más avanzadas turbinas
del ciclo de Rankine. Presenta un coste más bajo que las sales de nitrato. Poseen la
mayor temperatura estable de operación. Presenta buenas propiedades a la hora de la
transmisión de energía por transporte, posee un alto coeficiente de transferencia de
calor, elevada densidad, bajas presiones de trabajo y elevado calor especifico. La
desventaja principal es su elevado punto de congelación que es de 220ºC que limita su
uso en CCP debido a los riesgos de congelación debido a la mayor extensión de la red
de tuberías que tienen este tipo de centrales.
28
Optimización del espesor de aislamiento
La sal fundida empleada en la planta Archimede, que es la planta CCP que se está
tomando como referencia, en el presente proyecto es Solar Salt, compuesta por un 60%
de NaNO3 nitrato sódico y un 40% de KNO3 nitrato potásico, y presenta la mayor
temperatura de funcionamiento. Cuyas propiedades en función de la temperatura, se
comentará en próximos apartados. [28]
29
Optimización del espesor de aislamiento
Las sales fundidas deben ser transportadas desde cada uno de los lazos (formados
por captadores) hasta el sistema de potencia mediante un conjunto de tuberías
interconectadas entre sí. Es decir, dos tuberías principales, una para la entrada del fluido
frio y otra para el retorno del fluido caliente. Cada lazo se conectara a ambas para la
toma y retorno del fluido.
Se dota al campo de una cierta inclinación que permita la limpieza y purga de las
tuberías, para evitar el ensuciamiento y consecuentes obstrucciones producidas por
posibles solidificaciones debidas al elevado punto de congelación que poseen las sales
fundidas, en torno a los 238ºC.
Los diámetros de las tuberías cambiaran a lo largo del recorrido por los distintos
tramos, porque a medida que la sal fundida va bifurcándose hacia los sucesivos lazos, el
diámetro de las tuberías se irá reduciendo para mantener la velocidad de la sal
constante.
30
Optimización del espesor de aislamiento
Las tuberías de sales fundidas deben estar preparadas para soportar elevadas
temperaturas y los efectos de la corrosión. Por ello se emplearán tuberías de acero
inoxidable austeníticos.
Los aceros inoxidables austeníticos (Serie 300) son los más utilizados por su amplia
variedad de propiedades, se obtienen agregando Níquel a la aleación, por lo que la
estructura cristalina del material se transforma en austenita y de ahí su nombre. El
contenido en Cromo varía de 16 al 28%, el de Níquel de 3.5 al 22% y el de Molibdeno
1.5 al 6%. Los tipos más comunes son el AISI 304, 304L, 316, 316L, 310 y 317. Sus
propiedades básicas son: excelente resistencia a la corrosión, excelente factor de
higiene-limpieza, fáciles de transformar, excelente soldabilidad, no se endurecen por
tratamiento térmico, se pueden utilizar tanto a temperaturas criogénicas como a elevadas
temperaturas. De las ventajas del acero inoxidable, aparte de las correspondientes a los
tubos de acero al carbono, cabe destacar, la buena resistencia a la acción de ciertos
agentes corrosivos como atmosferas industriales. Adecuados para ambientes húmedos y
elevadas temperaturas. El inconveniente es que su precio es más caro frente a los acero
al carbono.
31
Optimización del espesor de aislamiento
En la planta Archimede sobre la que se basa este proyecto se han empleado tuberías
de acero inoxidable austenítico AISI 316L que presenta las mejores propiedades contra
la corrosión y soportar la temperatura de 600ºC de operación,
Propiedades del AISI 316L (L= reducción del contenido de carbono): Altamente
resistente a la corrosión, notable en medios ácidos. Soldable. Dúctil a todas las
temperaturas. Buena resistencia a la oxidación hasta una temperatura de 900ºC (1652ºF)
pero no se garantizan las propiedades a la deformación por encima de 500ºC. Resiste a
la corrosión, no de una manera general, especialmente cuando se trata de una corrosión
por picaduras. Los elementos que producen este tipo de corrosión son el flúor, cloro,
bromo y yodo (halógenos). El cloro es el elemento más común de todos ellos presente
en ambientes marinos y agua potable. Sin embargo después de un periodo de tiempo
expuesto a 550ºC, el acero AISI 316L se oxida más (como se observa en la imagen), el
espesor de la pared del material habrá disminuido y por tanto resistirá menos ante un
mismo esfuerzo. Susceptible a la rotura por corrosión bajo tensiones. Su coeficiente de
dilatación térmica es menor. El coeficiente medio de expansión térmica del acero
inoxidable austenítico es de 17.46·10-6 /ºC para una temperatura máxima de 500ºC. [12]
32
Optimización del espesor de aislamiento
A la hora de dimensionar la red de tuberías hay que recordar del mismo modo que el
principal inconveniente que presenta el uso de sales fundidas como fluido de trasporte
de calor, es su alto punto de congelación, por ello se emplean sistemas de prevención
para que no se produzca la solidificación de las sales. Este sistema se compone de una
serie de resistencias distribuidas a lo largo de toda la red de tuberías, que mantienen la
temperatura por encima del punto de congelación. Mantener las tuberías de acero
calientes con una temperatura de 290ºC reduce las cargas térmicas y dilataciones en la
red de tuberías pero los consumos parasitarios de las resistencias térmicas disminuyen la
eficiencia de la central y ponen en entredicho su viabilidad. Para evitar el consumo de
estos sistemas de calentamiento auxiliares que disminuirían el rendimiento de la planta,
se emplean también recubrimientos aislantes por toda la red de tuberías que limitan la
transferencia de calor desde el fluido hacia el ambiente.
33
Optimización del espesor de aislamiento
3.4. AISLAMIENTOS
Los materiales aislantes se definen como tales por una propiedad física que expresa
la facilidad o dificultad con que el calor atraviesa el material, la conductividad térmica.
Se pueden clasificar según su composición en dos grandes grupos: minerales y
plásticos. [35]
Este proyecto que estudia del espesor óptimo del aislamiento necesario para una
planta termosolar cilindro-parabólico que emplea sales fundidas a 550ºC como fluido
caloportador, se podría emplear lana de roca y lana cerámica como materiales de
aislamiento. La lana cerámica resulta más costosa que la lana de roca, pero presenta
mejores propiedades con espesores menores.
Hay que tener en cuenta que la distancia mínima, incluyendo el aislamiento, entre
las superficies exteriores de las tuberías de entrada del fluido (fría) y de retorno
(caliente) ha de ser mayor de 100mm. Del mismo modo la distancia con el resto de
34
Optimización del espesor de aislamiento
objetos también tendrá que ser mayor de 100mm. Para que pueda considerarse
despreciable el flujo de calor por radiación entre las superficies, debido a que las
superficies se encontraran a distinta temperatura. En la instalación bajo estudio en este
proyecto, la distancia entre ambas tuberías es mayor de 100mm.
Presentan una muy buena durabilidad, con muy bajo coste de mantenimiento, dado
que la composición es una combinación de rocas y fibras inorgánicas.
35
Optimización del espesor de aislamiento
Pero aun sin riesgo de explosión como es el caso que se plantea en la planta de
Archimedes por emplear sales fundidas como fluido de trabajo, el elevado punto de
fusión de las lanas minerales permite evitar daños estructurales en la instalación en caso
de incendio.
Las lana de roca y cerámica son productos 100% reciclable por lo que permiten ser
empleados en proyectos que pretenden no ser perjudiciales para el medio ambiente
como es el caso de las plantas de energías renovables que se analiza en este proyecto.
36
Optimización del espesor de aislamiento
PROROX PS 970(680ºC)
PROROX WM 970(680ºC)
Las coquillas poseen mayores ventajas frente a las mantas armadas, puesto que
permiten mayor rapidez y mayor facilidad de instalación. Están formadas por una capa
37
Optimización del espesor de aislamiento
38
Optimización del espesor de aislamiento
PROROX WM 940 es una manta armada de lana de roca volcánica cosida por una
cara a una malla de acero galvanizado con hilos también de acero galvanizado. También
puede solicitarse con hilo y malla de acero inoxidable. Su temperatura máxima de
operación es 600ºC siendo su valor punta máximo de 750ºC.Ventajas que presenta son
el buen rendimiento térmico y gran resistencia a la humedad.
Otro modelo como la manta armada PROROX WM 960 permite soportar mayores
temperaturas, siendo la temperatura de trabajo de 660ºC y una temperatura máxima de
750ºC. Pero su valor se ve incrementado en 2 euros por metro cuadrado. Dicha
temperatura de operación no se alcanzara en ningún punto dado que la máxima
temperatura de las sales es de 600ºC antes de que se vuelvan inestables. No es necesario
el empleo de la manta armada PROROX WM 960 para plantas cilindro-parabólicas con
sales fundidas como fluido de trabajo.
Dado que las propiedades de los tres productos son semejantes, se tomará la más
económica, PROROX WM 940 debido a que cumple con los requisitos y permite un
ahorro en el coste del aislamiento.
39
Optimización del espesor de aislamiento
Se puede observar en la tabla que las coquillas presentan mejores propiedades que
las mantas armadas pero debido a la necesidad de realizar en este proyecto un estudio
detallado del espesor óptimo necesario para el aislamiento emplearemos mantas
armadas, las cuales permiten obtener varios espesores por superposición de capas. Otro
motivo es que los diámetros de tubería para el empleo de coquillas son menores de los
necesarios en la red de tuberías bajo estudio.
40
Optimización del espesor de aislamiento
Para la realización del análisis del aislamiento de este proyecto, se toma la lana
cerámica del fabricante UNIFRAX, denominada Fiberfrax Durablanket S. Que cumple
con todas características necesarias al menor coste. Siendo su temperatura máxima de
operación 1250ºC muy por encima de la temperatura de operación de las sales
fundidas.[18]
Realizando una comparativa entre la lana cerámica y la lana de roca cabe mencionar
que la lana cerámica presenta el beneficio adicional de soportar mayores temperaturas
de empleo que la lana de roca y una menor conductividad con el mismo espesor. De
manera que para el aislamiento de una tubería, se requeriría un espesor de lana cerámica
menor que de lana de roca, y con ello un menor peso de la estructura, que permitiría
disminuir los mecanismos de soporte necesarios. Ante el mismo espesor, la lana
cerámica representa un mayor coste de aislamiento que la lana de roca. Se puede
emplear una combinación de ambos materiales, aplicando inicialmente una capa de lana
cerámica y en sobre ella una capa de lana de roca. Todo ello se analiza en aparatados
posteriores de este proyecto.
3.4.4. INSTALACIÓN
Para realizar la instalación hay que cortar la manta y la malla de acero a la longitud
adecuada, para que encaje con el tubo con una presión ligera previa. Las uniones de
cierre deben alternarse a un ángulo de al menos 30 grados una de otra. Las juntas tope
de las mantas (longitudinales y circulares) deben unirse con alambre de acero (mín. de
0,5 mm) o asegurarse con ganchos.
41
Optimización del espesor de aislamiento
Los tubos de acero inoxidable y los tubos con una temperatura >350 deben por
preferencia aislarse con un producto en el que tanto la malla como el alambre de
costura sean de acero inoxidable, que puede ser solicitado al fabricante. Si las colchas se
ensamblan en multicapas, deben alternarse las uniones longitudinales y circulares
("unión de mampostería").
42
Optimización del espesor de aislamiento
acero. Las conexiones a los montajes, las tapas de cabezales y partes terminales, etc.,
deben de ser herméticas, mediante el sellador adecuado. Todos los componentes de
acero expuestos al ambiente corrosivo se deben limpiar, desgrasar y cubrir con un
terminado protector. [35]
43
Optimización del espesor de aislamiento
44
Optimización del espesor de aislamiento
45
Optimización del espesor de aislamiento
3.5. CAPTADOR
La energía solar se concentra gracias a los espejos que la reflejan al tubo receptor
situado en la línea focal del canal parabólico. El tubo receptor se compone de un tubo
metálico y una carcasa de vidrio concéntrica al tubo metálico. En el espacio entre la
carcasa y el tubo metálico se realiza el vacio.[34]
46
Optimización del espesor de aislamiento
El Sol puede considerarse desde la tierra, como un foco térmico a 5777K que nos
hace llegar hasta 1367 W/m2 en forma de radiación electromagnética. La energía solar
es abundante, de alta calidad, variable en el tiempo y en el espacio. La irradiación es la
radiación que incide durante un periodo de tiempo sobre una superficie medida en
Wh/m2 (Energía). La irradiación solar recibida por una superficie está influenciada por
el ángulo de incidencia que forman los rayos del Sol con la normal a la superficie.
Un año típico (modelo) está constituido por años de referencia. Un año de referencia
de datos meteorológicos es un conjunto de series temporales de valores de radiación
solar y otras variables meteorológicas de un año. Se compone de meses seleccionados
47
Optimización del espesor de aislamiento
entre los distintos años de medida disponibles (al menos 10 años) concatenados de
forma que componen un año completo.
Los datos meteorológicos horarios han sido obtenidos de SAM (System Advisor
Model) para Sevilla, España.
Las propiedades del aire que rodea a la red de tuberías se pueden obtener de
cualquier libro de termodinámica. Pero en este caso para la realización del proyecto, el
programa empleado EES, permite determinar las propiedades al aire a una presión y
temperatura determinada. Dichos parámetros serán empleados para el cálculo del
intercambio de calor por convección libre con el ambiente. La temperatura del aire en la
superficie será la temperatura media de película entre la temperatura ambiente y la
temperatura de la superficie de la tubería, del orden de 40ºC.
48
Optimización del espesor de aislamiento
El múltiplo solar (MS) se define como la potencia térmica total del campo solar
( dividida por la potencia térmica requerida por el ciclo de potencia para
producir su potencia eléctrica nominal ( é .
49
Optimización del espesor de aislamiento
La central de Archimede sobre la que se basa este proyecto, posee según sus
características técnicas, 7 horas de almacenamiento térmico, estando el campo solar
constituido por 9 lazos de 6 captadores.
En el análisis del múltiplo solar a lo largo del día, se observa que cuando el valor del
múltiplo solar en las primeras horas del día es menor que la unidad, para la operación de
la planta se recurre al almacenamiento térmico para producir la potencia térmica
requerida para el ciclo de potencia ( é ).
Al medio día se suceden las horas en la que se obtiene la mayor radiación del día
dado que el campo se encuentra sobredimensionado, el múltiplo solar es mayor que la
unidad. En el momento que el múltiplo solar es mayor que la unidad comienza el
almacenamiento térmico.
Próximo a las horas de la puesta de sol el múltiplo solar vuelve a ser menor que la
unidad y es necesario recurrir al almacenamiento térmico para producir la potencia
térmica necearía é La energía almacenada no se utiliza toda ella tras la puesta
50
Optimización del espesor de aislamiento
Este proyecto tiene como objetivo el análisis del óptimo aislamiento de la red de
tuberías que conduce el fluido caloportador por el campo solar hasta los diferentes lazos
y lo recoge a la salida de estos.
51
Optimización del espesor de aislamiento
Figura 37. Fotografía campo solar Archimede. Imagen del uso de liras.
De modo que la longitud de cada tramo de la red de tuberías bajo estudio será de
54m hasta la intersección, tanto para la red fría como para la red caliente. Dado que el
campo solar está formado por 9 lazos, y uno de ellos se encuentra situado en la zona
central, se divide el campo por simetría. Pudiendo distinguir 4 tramos de la red de
tuberías en los que se tienen distintos caudales debido a las intersecciones con cada uno
de los lazos. Un tramo final recoge todas las sales y las conduce hasta la isla de
potencia. [28]
54
52
Optimización del espesor de aislamiento
La red fría (azul) conduce el fluido de trabajo del tanque fría a los distintos lazos y
la red caliente (roja) recoge el fluido a la salida de los lazos y lo conduce a la isla de
potencia.
Figura 39. Imagen entrada y salida de los lazos de campo solar de Arhimede
53
Optimización del espesor de aislamiento
La sal fundida empleada en la central de Archimede como fluido de trabajo, tiene las
siguientes características técnicas recogidas en la ficha técnica de ENEA, que varían
con la temperatura.
Densidad: [21]
2090 0,636 º
1443 0,172 º
Viscosidad [21]:
0,443 1,9 10 º
º
54
Optimización del espesor de aislamiento
Es necesario conocer la temperatura en cada una de las partes del campo debido a
que las propiedades de la sal fundida dependen de la temperatura, como se ha
comentado en el apartado 4.4. Las sales fundidas tienen un rango de operación, fuera
del cual se producen problemas de descomposición o solidificación. Para mantener
dicho rango y obtener la máxima temperatura de operación es necesario disminuir las
pérdidas durante el recorrido de las sales por el campo solar. Para el cálculo de dichas
pérdidas es necesario conocer las temperaturas iniciales o temperaturas fijas de la
instalación. [28]
290º 5 Í
55
Optimización del espesor de aislamiento
La temperatura de salida del campo solar será algo menor que la de salida de los
captadores debido a las pérdidas ocasionadas en su recorrido por la red de tuberías
desde los captadores al tanque de almacenamiento, y que serán estudiadas en el
proyecto, en función del espesor y el tipo de aislamiento empleado.
La temperatura de las sales a la salida del lazo es de 550ºC, para ello es necesario
variar el caudal de sales en el interior de los lazos en función de la radiación solar, es
56
Optimización del espesor de aislamiento
550°
57
Optimización del espesor de aislamiento
Potencia útil = potencia procedente del sol (G_s) – pérdidas de calor hacia el
exterior (
.
1 3 / →
Potencia útil →
La velocidad en el interior de las tuberías oscila entre [1-3 m/s] dentro de dicho
intervalo, a mayor velocidad del fluido, se obtiene mayor potencia útil y menor
temperatura a la salida del captador.
↓ → ↓ ↑ → ↑
Por el contrario, ante una misma radiación solar, a mayor velocidad del fluido se
obtiene una mayor potencia útil y una menor temperatura a la salida de los captadores.
↑ → ↑ ↓ → ↓
Este hecho se explica debido a que una mayor velocidad ocasiona que el fluido no
esté el tiempo necesario circulando por el interior de los captadores y por tanto no logre
58
Optimización del espesor de aislamiento
1; 3 / →
/4
59
Optimización del espesor de aislamiento
Para la realización de los cálculos del presente proyecto se ha observado que dicha
variación se puede considerar lineal según la siguiente ecuación:
/ 0,356 0,0073 /
60
Optimización del espesor de aislamiento
2 2
La caída total será igual a la suma de las caídas de presión en cada tramo
∆ ∆
Para obtener la pérdida de presión primaria total se realizará el cálculo por tramos
horizontales de sección constante.
; →
. . → .
2
Existen multitud de tablas, curvas, ecuaciones, etc. para obtener el valor del factor
de fricción (f). Sin embargo, el más empleado es el “Diagrama de Moody”
61
Optimización del espesor de aislamiento
0,184
/
→ ó 2 10
.
2
Donde “K” es el factor de pérdidas de carga secundarias que depende del tipo de
accesorio. El valor de dicho parámetro se obtiene en función de la geometría del
elemento, dichos valores vienen proporcionados por el fabricante.
El valor del factor de pérdidas de carga correspondiente a los 4 codos de 90º que se
encuentran en cada tramo de la red de tuberías debido a la necesidad de instalar liras.
Codos roscados
° 1,5 → ° 4 ∗ 1,5 6
Pérdidas de carga asocias a las válvulas en la entrada de cada uno de los lazos, y por
tanto una en cada tramo.
0,15
El valor de las pérdidas secundarias será del 1Pa, es decir mucho menores que las
pérdidas primarias que son del orden de 2·104 Pa
El estudio de este proyecto se basa en los tramos de tuberías rectos despreciando por
tanto las pérdidas de carga secundarias ocasionadas debidas a otros elementos, por no
tener información detallada de la geometría y fabricante de los elementos empleado
como son las válvulas, estrechamientos, codos, bifurcaciones y expansiones.
62
Optimización del espesor de aislamiento
4.8. DIÁMETRO
Dimensiones y pesos de los tubos, según las normas ASME B 36.10; B 36.19.
Dadas las condiciones de presión y temperatura de la tubería se toma Sh10S.En la tabla
mostrada a continuación obtenida del catálogo de un fabricante, se observa que en negro
aparecen los espesores, expresados en milímetros y en color rojo el correspondiente
peso expresado en kg/m.
63
Optimización del espesor de aislamiento
64
Optimización del espesor de aislamiento
5. CÁLCULOS TÉRMICOS
5.1. ESQUEMA
El sistema bajo estudio en este proyecto es una tubería de acero inoxidable por la
que circula un fluido con una temperatura , que genera un coeficiente de película
interior . El fluido en su recorrido experimenta una pérdida de temperatura. [1]
La tubería se encuentra aislada con un espesor “x” compuesto por una capa de lana
cerámica y el doble de lana de roca. El sistema se encuentra en el exterior rodeado por
aire bajo condiciones ambientales a una temperatura exterior .
65
Optimización del espesor de aislamiento
Para analizar las pérdidas térmicas es necesario determinar las resistencias eléctricas
de los distintos mecanismos de transferencia de calor.
∆
∑
66
Optimización del espesor de aislamiento
ρvD 4 m
Re
μ πμD
Donde “v” es la velocidad media del fluido sobre la sección transversal del tubo, y
“D” es el diámetro del tubo. Para determinar el Reynolds del fluido se determinan las
propiedades del fluido de trabajo a una temperatura media de masas, que se aproxima
inicialmente a la temperatura de entrada del fluido en el conducto, puesto que el
objetivo es disminuir la diferencia de temperatura entre la entrada y la salida.
1
R
h A
67
Optimización del espesor de aislamiento
En primer lugar, es preciso determinar el valor del número de Reynolds, para poder
determinar el número de Nusselt con el que se determinará el coeficiente de película
interior. Se emplea la correlación de Gnielinski para determinar el valor del número de
Nusselt, que se aplica para convección forzada en flujo interno turbulento
completamente desarrollado con un valor del número de Reynolds comprendido entre
[3000; 5·106]. Con la condición de que el conducto sea suficientemente largo, es decir
(L / D)>10- Y la condición de número de Prandtl este comprendido en el siguiente
rango: 0.5 ≤ Pr < 2000. [2]
/8 1000
ó →
1 12,7 / 1
/
2 10 → 0,316
/
2 10 → 0,184
En los casos en los que la radiación es menor de 150W/m2 el caudal es del orden de
0,23kg/s lo que implica que el número de Reynolds alcance valores menores de 2300,
en los tramos más alejados, que supondrían un régimen laminar. En dicho casos el valor
del número de Nuseelt para régimen laminar completamente desarrollada con
temperatura superficial constate (temperatura ambiente) y Pr>0,6 será [2]:
3,66
El coeficiente de película interior depende del régimen de carga, puesto que depende
del caudal que circula por su interior en función de las variaciones meteorológicas. Por
lo que se pone de manifiesto la variación de las pérdidas térmicas con el régimen de
carga. La estimación de las pérdidas de generación anual se realizará con las
condiciones de operación en el año modelo a intervalos de una hora.
68
Optimización del espesor de aislamiento
Las temperaturas más altas se asocian con las energías moleculares más altas y
cuando las moléculas vecinas chocan, como lo hacen constantemente, ocurre una
transferencia de energía de las moléculas más energéticas a las menos energéticas. En
presencia de un gradiente de temperatura, la transferencia de energía por conducción
ocurre en la dirección de la temperatura decreciente.
D
ln D
R
2πL k
k 500º 21 /
∆
R
69
Optimización del espesor de aislamiento
Se trata de una resistencia térmica conductiva, que al igual que la tubería de acero
tendrá geometría cilíndrica porque el aislante se colocará recubriendo la tubería de
acero.
D
ln D
R
2πL k
D
ln D
R
2πL k
0,052 0,00048 T
0,14 /
Una menor conductividad implica una mayor resistencia térmica, y con ello una
70
Optimización del espesor de aislamiento
600º ; 96 0,14 /
600º 0,236 /
600º ; 96 ≅ 400º
La convección libre se origina cuando una fuerza actúa sobre un fluido en el que hay
gradientes de densidad. El efecto neto es una fuerza de empuje, que induce corrientes de
convección libre. El gradiente de densidad se debe al gradiente de temperatura, y la
fuerza de cuerpo se debe al campo gravitatorio.
71
Optimización del espesor de aislamiento
Como las velocidades de flujo de convección libre son por lo general mucho más
pequeñas que las que se asocian con la convección forzada, las transferencias de calor
por convección correspondientes también son más pequeñas.
1
R
h A
h T ; D
A D
g β D T T
Ra Pr
υ
72
Optimización del espesor de aislamiento
Nu C Ra
h D
Nu → h
k
T T
R
73
Optimización del espesor de aislamiento
Puede obtenerse una representación del valor de la temperatura efectiva del cielo en
la página web del Código Técnico de la Edificación en formato .MET, cuyos valores
para el caso de la Península Ibérica se encuentra en rango [-17,7; 25,7]ºC
correspondientes a un clima de referencia de un año tipo.
Correlación de Walton:
,
273,15 °
273,15 /
,
0,787 0,764 ln 1 0,0224 0,0035
0,00028
Se supone un valor fijo de la temperatura efectiva del cielo para la realización de los
cálculos de este proyecto por no tener datos más detallados de las condiciones
ambientales del año modelo. Para ello se considera una temperatura ambiente seca “ "
de 25ºC y una humedad relativa “HR” del 50%.
25º
14,31º
74
Optimización del espesor de aislamiento
El flujo de calor emitido por la superficie originado por la energía térmica, es menor
que el de un cuerpo negro a la misma temperatura. En el caso de un cuerpo negro la
emisividad es igual a la unidad
A T T
W
Siendo σ=5,67·10-8 la constante de Stefan Boltzmann [2].
m2 K4
75
Optimización del espesor de aislamiento
Hay que tener en cuenta a la hora de realizar los cálculos numéricos que para
determinar el calor por radiación el valor de la temperatura debe estar expresado en
Kelvin dado que la constante de Boltzmann está expresada en Kelvin y se trata de una
diferencia de temperaturas elevadas a la cuarta.
300K 0,08
300K 0,03
R 4 ε σ A ESTIMACIÓN
h ε σ A T T T T
76
Optimización del espesor de aislamiento
Q
0→ ≅ 1
En este proyecto el valor del radio crítico es del orden de 1cm de manera que si se
quiere aislar la red de tuberías, disminuyendo las pérdidas de calor, se debe tomar
espesores de aislamiento mayores a 1 cm.
↑ ; ↓ → ¿ ?
77
Optimización del espesor de aislamiento
5.3.1. DATOS
Se realizan los cálculos para varios espesores de aislamiento, y obtener así una curva
final que represente la variación de la temperatura del fluido en función del espesor de
aislamiento
Radiación solar →
Caudal →
Área exterior →
5.3.2. INCÓGNITAS
78
Optimización del espesor de aislamiento
de aislamiento.
5.3.3. HIPÓTESIS
→
CRITERIO 1
60 °
79
Optimización del espesor de aislamiento
CRITERIO 2
80
Optimización del espesor de aislamiento
5.3.5. DESARROLLO
j=1;
TRAMO j
Diámetros normalizados
≅ NO
j=j+1 _
j>5 NO
DIAMETROS TRAMOS
81
Optimización del espesor de aislamiento
82
Optimización del espesor de aislamiento
:
550º ; 54 ; ;
e=10 ; i=1;
e=e
5 /
; ;
NO
i=i+1
NO
i>5
e=e+10
NO
e>80
º ;
83
Optimización del espesor de aislamiento
SISTEMA DE ECUACIONES
R R
Incógnitas= ;
0
" A
" " 0
84
Optimización del espesor de aislamiento
" "
" T – T T T
" T – T T T
" " " /
Incógnitas ; ;
El valor del coeficiente de película exterior estará en el rango de [2-25] W/m2K por
tratarse de un proceso de convección libre de gases. [29]
2; 25
1º Ó → 5
85
Optimización del espesor de aislamiento
6. ANÁLISIS DE RESULTADOS
Para determinar el caudal nominal por tramo se realiza un estudio del rendimiento
de la planta con los datos proporcionados.
Número de captadores = 54
m ∆ → m 70,87 /
Sabiendo que el caudal en los nueve lazos del campo solar es el mismo, se puede
determinar el caudal nominal necesario por cada uno de los lazos
7,8 /
64
→ 1,34 /
4
86
Optimización del espesor de aislamiento
420º 1823 /
420º 0,00163
85533
El caudal en cada uno de los tramos de la red de tuberías irá variando conforme se
produzcan las intersecciones con los lazos, de manera que en las tuberías de entrada (red
fría) el caudal irá disminuyendo, y en las de retorno (red caliente) el caudal irá
aumentando tras cada intersección con un lazo.
Dada la configuración del campo solar, formado por nueve lazos divididos
simétricamente respecto a la isla de potencia, se simplifica el estudio térmico de la red
de tuberías considerando solo 5 tramos de 54m cada uno.
La longitud de los tramos es igual a la distancia entre los lazos considerando la lira.
L 54
TRAMO 1: es el más alejado del ciclo de potencia por lo que el caudal será el
menor de todos los tramos.
TRAMO 5: se trata del tramo que conduce todas las sales del campo solar.
TRAMO Caudal nominal [kg/s]
1 7,8
2 15,6
3 23,4
4 31,2
5 70,8
87
Optimización del espesor de aislamiento
Realizando los cálculos del apartado 4.6 y 4.7, manteniendo la velocidad constante
en el interior de los tramos, e igual a la velocidad nominal en el interior de los lazos. Se
obtienen los diámetros interiores que aparecen reflejados en la tabla y las pérdidas de
carga, primarias (P1) y secundarias (P2).
Caudal Diámetro
TRAMO v [m/s] Reynolds P1[Pa] P2[Pa] AP [Pa] AP [bar] f
[kg/s] interior [m]
Una vez determinado el diámetro en cada uno de los tramos con velocidad
constante, se procede a buscar en las tablas del fabricante (apartado 4.8), la tubería que
mejor se ajuste a los diámetros calculados, debido a que los diámetros se encuentran
normalizados. En esta búsqueda se tendrá en cuenta la presión máxima y la temperatura
de trabajo para determinar el espesor de la tubería. Estos valores se sobre dimensionarán
para prevenir cualquier anomalía en el funcionamiento que provoque un fallo en toda la
instalación. Considerando las tablas del fabricante de tuberías [4]
88
Optimización del espesor de aislamiento
Observaciones:
El número de Reynolds en todos los casos es mayor que 10000 lo que implica un
régimen turbulento completamente desarrollado que facilitará la transferencia de calor
en los captadores.
Las pérdidas de carga son del orden de 0,15 bares. Se comprueba que las pérdidas de
carga secundarias son menores que las primarias. El valor de las pérdidas secundarias es
aproximado dado que no se tiene detalle de las dimensiones de los elementos que las
provocan.
89
Optimización del espesor de aislamiento
TRAMO 1: 67 ; 3,05 ;
TRAMO 2: 83 ; 3,05 ; 2
TRAMO 3: 108 ; 3,05 ; 3
TRAMO 4: 108 ; 3,05 ; 4
TRAMO 5: 162 ; 3,4 ; 9
SIN AISLAMIENTO ‐ RED CALIENTE
TRAMO TRAMO 1 TRAMO 2 TRAMO 3 TRAMO 4 TRAMO 5 TOTAL
DIAMETRO (mm) 67 83 108 108 162
Pérdidas [MW] 491,63 603,72 723,27 710,06 1000,31 3528,99
Max Tª Superficie [ºC] 538,10 541,40 536,60 534,50 532,10 >60ºC
Max ΔT [ºC] 100,80 65,36 47,36 33,57 20,21
Min Tª salida[ºC] 449,2 484,6 459,1 448,2 439,3
Max Tª salida[ºC] 536,7 542 538,1 536,2 534,6 <550ºC
90
Optimización del espesor de aislamiento
RED FRÍA: Conduce la sal fundida desde la salida del tanque frío hasta la entrada
de los diferentes lazos
SIN AISLAMIENTO – RED FRÍA
TRAMO TRAMO 1 TRAMO 2 TRAMO 3 TRAMO 4 TRAMO 5 TOTAL
DIAMETRO (mm) 67 83 108 108 162
Pérdidas [MW] 96,84 115,75 142,09 144,23 209,27 708,17
Max Tª Superficie [ºC] 282,90 282,80 285,20 286,60 287,90 >60ºC
Max ΔT [ºC] 20,42 13,24 11,31 8,97 5,98
Min Tª salida[ºC] 230,10 250,50 263,70 275,10 284,00 <260ºC
En el caso sin aislamiento de la red fría, la temperatura del último lazo, que sería el
más desfavorecido (tramo 1) disminuye entre 7ºC y 60 ºC.
91
Optimización del espesor de aislamiento
Tª salida
captador 550ºC
92
Optimización del espesor de aislamiento
Max Pérdidas térmicas [W/m]
Espesor [cm] TRAMO 1 TRAMO 2 TRAMO 3 TRAMO 4 TRAMO 5
10 236,11 267,33 313,2 312,83 409,7
20 165,13 183,26 209,69 209,52 264,26
30 137,22 150,57 169,89 169,78 209,26
40 121,69 132,52 148,11 148,04 179,56
50 111,54 120,85 134,09 134,04 160,65
60 104,28 112,54 124,19 124,13 147,37
70 98,76 106,26 116,72 116,67 137,46
80 94,39 101,28 110,85 110,81 129,74
450,00
Tª salida 400,00
Perdidas térmicas(W/m)
93
Optimización del espesor de aislamiento
Pérdidas térmicas [MWh/año]
Espesor [cm] TRAMO 1 TRAMO 2 TRAMO 3 TRAMO 4 TRAMO 5 TOTAL
10 44,86 50,77 59,45 59,38 77,72 292,19
20 31,47 34,91 39,92 39,89 50,28 196,48
30 26,18 28,72 32,39 32,37 39,87 159,52
40 23,23 25,29 28,26 28,24 34,24 139,26
50 21,30 23,07 25,60 25,58 30,64 126,20
60 19,92 21,49 23,71 23,70 28,12 116,95
70 18,88 20,29 22,29 22,28 26,24 109,98
80 18,04 19,34 21,18 21,17 24,77 104,51
90,00
Pérdidas térmicas (MWh/año)
Tª salida 80,00
captador 550ºC 70,00
60,00
TRAMO 1
50,00
TRAMO2
40,00
TRAMO3
30,00
TRAMO4
20,00
TRAMO5
10,00
0,00
10 20 30 40 50 60 70 80
Espesor aislamiento (cm)
94
Optimización del espesor de aislamiento
Max Variación de temperatura [ºC]
Espesor [cm] TRAMO 1 TRAMO2 TRAMO3 TRAMO4 TRAMO5
10 10,87 6,19 4,80 3,58 2,08
20 7,65 4,26 3,24 2,42 1,35
30 6,37 3,51 2,63 1,97 1,08
40 5,66 3,09 2,30 1,72 0,92
50 5,19 2,82 2,08 1,56 0,83
60 4,86 2,63 1,93 1,44 0,76
70 4,61 2,48 1,81 1,36 0,71
80 4,40 2,37 1,72 1,29 0,67
12,00
Tª salida 10,00
captador 550ºC 8,00
TRAMO 1
AT (ºC)
6,00 TRAMO2
TRAMO3
4,00
TRAMO4
2,00 TRAMO5
0,00
10 20 30 40 50 60 70 80
Espesor aislamiento (cm)
95
Optimización del espesor de aislamiento
Tabla 10. TEMPERATURA DE SALIDA DEL CAMPO SOLAR (TRAMO 5) – RED CALIENTE
552,0
Tª salida 550,0
captador 550ºC 548,0
Tª sal TRAMO 1 (ºC)
546,0
544,0
542,0
MIN
540,0
MAX
538,0
536,0
534,0
532,0
10 20 30 40 50 60 70 80
Espesor de aislamiento (cm)
96
Optimización del espesor de aislamiento
Max Temperatura Superficie Exterior T sup [ºC]
Espesor [cm] TRAMO 1 TRAMO 2 TRAMO 3 TRAMO 4 TRAMO 5
10 96,25 98,89 102,10 102,10 107,20
20 74,23 72,49 74,23 74,22 77,18
30 63,01 63,86 64,97 64,97 66,92
40 59,12 59,71 60,51 60,50 61,91
50 56,87 57,31 57,92 57,91 58,99
60 55,42 55,77 56,25 56,25 57,11
70 54,42 54,71 55,09 55,09 55,80
80 53,69 53,93 54,26 54,26 54,85
120,00
Tª salida 110,00
captador 550ºC 100,00
Tª Superficie (ºC)
90,00 TRAMO 1
80,00 TRAMO2
70,00 TRAMO3
TRAMO4
60,00
TRAMO5
50,00
40,00
10 20 30 40 50 60 70 80
Espesor aislamiento (cm)
97
Optimización del espesor de aislamiento
Mínima Temperatura de salida [ºC]
Espesor [cm] TRAMO 1 TRAMO 2 TRAMO 3 TRAMO 4 TRAMO 5
10 539,1 543,8 541,1 539,7 538,8
20 542,3 545,7 543,9 543 542,4
30 543,6 546,5 545 544,3 543,9
40 544,3 546,9 545,6 545 544,6
50 544,8 547,2 546 545,5 545,1
60 545,1 547,4 546,3 545,8 545,5
80 545,4 547,5 546,5 546 545,8
70 545,6 547,6 546,7 546,2 546
550
Tª salida 548
captador 550ºC
Min Tª de salida (ºC)
546
544 TRAMO_1
542 TRAMO_2
540 TRAMO_3
TRAMO_4
538
TRAMO_5
536
534
10 20 30 40 50 60 70 80
Espesor de aislamiento (cm)
98
Optimización del espesor de aislamiento
Máxima Temperatura de salida [ºC]
Espesor [cm] TRAMO 1 TRAMO 2 TRAMO 3 TRAMO 4 TRAMO 5
10 548,8 549,3 549 548,9 548,8
20 549,2 549,5 549,4 549,3 549,2
30 549,3 549,6 549,5 549,4 549,3
40 549,4 549,7 549,5 549,5 549,4
50 549,4 549,7 549,6 549,5 549,5
60 549,5 549,7 549,6 549,6 549,5
80 549,5 549,7 549,6 549,6 549,6
70 549,5 549,7 549,6 549,6 549,6
549,8
Tª salida 549,6
captador 550ºC
549,4
Max Tª salida (ºC)
549,2 TRAMO_1
549 TRAMO_2
548,8 TRAMO_3
TRAMO_4
548,6
TRAMO_5
548,4
548,2
10 20 30 40 50 60 70 80
Espesor de aislamiento (cm)
99
Optimización del espesor de aislamiento
Tª entrada al
campo solar
290ºC
100
Optimización del espesor de aislamiento
Máxima Pérdidas térmicas [W/m]
Espesor [cm] TRAMO 1 TRAMO 2 TRAMO 3 TRAMO 4 TRAMO 5
10 79,7 89,9 105,0 105,4 137,1
20 57,2 63,2 71,9 72,1 89,9
30 48,3 52,7 59,1 59,3 72,2
40 43,3 46,9 52,1 52,3 62,7
50 40,0 43,1 47,6 47,8 56,5
60 37,6 40,4 44,4 44,5 52,2
70 35,8 38,3 41,9 42,1 49,0
80 34,4 36,7 40,0 40,1 46,4
160,0
140,0 Tª entrada al
campo solar
Pérdita témica (W/m)
120,0
290ºC
100,0 TRAMO 1
80,0 TRAMO 2
60,0 TRAMO 3
TRAMO 4
40,0
TRAMO 5
20,0
0,0
10 20 30 40 50 60 70 80
Espesor de aislamiento (cm)
101
Optimización del espesor de aislamiento
Tabla 15. PÉRDIDAS TERMICAS TOTALES DEL AÑO MODELO – RED FRÍA
Pérdidas térmicas Totales [MWh/año ]
Espesor [cm] TRAMO 1 TRAMO 2 TRAMO 3 TRAMO 4 TRAMO 5 SUMA
10 14,79 16,67 19,45 19,51 25,28 95,70
20 10,67 11,77 13,37 13,40 16,65 65,85
30 9,03 9,84 11,02 11,04 13,39 54,32
40 8,10 8,77 9,73 9,74 11,63 47,97
50 7,49 8,07 8,89 8,90 10,50 43,85
60 7,05 7,56 8,29 8,30 9,71 40,91
70 6,72 7,18 7,84 7,85 9,11 38,69
80 6,45 6,88 7,48 7,49 8,65 36,94
30,00
Tª entrada al
Pérdidas térmicas (MWh/año)
0,00
10 20 30 40 50 60 70 80
Espesor aislamiento (cm)
102
Optimización del espesor de aislamiento
Max Variación de temperatura [ºC]
Espesor [cm] TRAMO 1 TRAMO2 TRAMO3 TRAMO4 TRAMO5
10 3,64 2,09 1,64 1,24 0,72
20 2,66 1,48 1,13 0,86 0,47
30 2,26 1,24 0,93 0,71 0,38
40 2,03 1,11 0,83 0,62 0,33
50 1,88 1,02 0,75 0,57 0,30
60 1,77 0,96 0,70 0,53 0,28
70 1,69 0,91 0,67 0,50 0,26
80 1,62 0,87 0,64 0,48 0,25
4,00
Tª entrada al
3,50 campo solar
3,00 290ºC
2,50 TRAMO 1
AT (ºC)
2,00 TRAMO2
1,50 TRAMO3
TRAMO4
1,00
TRAMO5
0,50
0,00
10 20 30 40 50 60 70 80
Espesor aislamiento (cm)
103
Optimización del espesor de aislamiento
Tabla 17. TEMPERATURA SALIDA DEL TRAMO 1 (Tramo más alejado) – RED FRÍA
292 Tª entrada al
290 campo solar
290ºC
Tª sal TRAMO 1 (ºC)
288
286
284
MIN
282 MAX
280
278
276
10 20 30 40 50 60 70 80
Espesor de aislamiento (cm)
104
Optimización del espesor de aislamiento
Max Temperatura Superficial Exterior T sup [ºC]
Espesor [cm] TRAMO 1 TRAMO2 TRAMO3 TRAMO4 TRAMO5
10 67,16 68,10 69,25 69,28 71,03
20 57,59 58,05 58,66 58,67 59,68
30 54,51 54,81 55,19 55,19 55,85
40 53,06 53,26 53,53 53,53 54,00
50 52,22 52,37 52,58 52,58 52,93
60 51,69 51,81 51,97 51,97 52,25
70 51,32 51,42 51,55 51,55 51,78
80 51,06 51,14 51,25 51,25 51,44
75,00
Tª entrada al
70,00 campo solar
290ºC
Tª Superficie (ºC)
65,00
TRAMO 1
60,00
TRAMO2
55,00
TRAMO3
50,00 TRAMO4
TRAMO5
45,00
40,00
10 20 30 40 50 60 70 80
Espesor aislamiento (cm)
105
Optimización del espesor de aislamiento
Mínima Tª de salida [ºC]
Espesor [cm] TRAMO 1 TRAMO 2 TRAMO 3 TRAMO 4 TRAMO 5
10 280,7 284,3 286,4 288 289,3
20 283,4 286,1 287,5 288,7 289,5
30 284,5 286,7 288 288,9 289,6
40 285,1 287,1 288,2 289 289,7
50 285,5 287,4 288,4 289,1 289,7
60 285,8 287,5 288,5 289,2 289,7
80 286 287,7 288,6 289,2 289,7
70 286,1 287,8 288,6 289,3 289,8
292 Tª entrada al
290 campo solar
290ºC
Mínima Tª salida (ºC)
288
286 TRAMO_1
284 TRAMO_2
282 TRAMO_3
TRAMO_4
280
TRAMO_5
278
276
10 20 30 40 50 60 80 70
Espesor de aislamiento (cm)
106
Optimización del espesor de aislamiento
Máxima Tª de salida [ºC]
Espesor [cm] TRAMO 1 TRAMO 3 TRAMO 3 TRAMO 4 TRAMO 5
10 289 286,1 287,5 288,6 289,5
20 289,3 287,3 288,3 289,1 289,7
30 289,4 287,8 288,6 289,3 289,7
40 289,5 288 288,8 289,3 289,8
50 289,5 288,2 288,9 289,4 289,8
60 289,6 288,3 289 289,4 289,8
80 289,6 288,4 289 289,5 289,8
70 289,6 288,5 289,1 289,5 289,8
291
Tª entrada al
290 campo solar
Maxima Tª salida (ºC)
290ºC
289
TRAMO_1
288
TRAMO_2
287
TRAMO_3
286 TRAMO_4
TRAMO_5
285
284
10 20 30 40 50 60 80 70
Espesor de aislamiento (cm)
107
Optimización del espesor de aislamiento
El uso de aislamiento no elimina las pérdidas térmicas hacia el exterior tan sola las
reduce. Un mayor espesor de aislamiento supone una mayor resistencia térmica que
implica menores pérdidas térmicas.
↑
↑ ↑
→
↑
↓ ∆ → ↑
Por el contrario, el tramo1 tiene el menor caudal y el menor diámetro de toda la red,
de manera que se producirán menores pérdidas térmicas, mayor variación de
temperatura y menor temperatura superficial
q
→ T .
m m
∆ ∆ → T
q ≅q
→ ∆ ∆ → T
m m T ≅T
108
Optimización del espesor de aislamiento
CRITERIO TÉRMICO 1
é 1: 60°
CRITERIO TÉRMICO 2
∆ ∆ →
109
Optimización del espesor de aislamiento
Hay que recordar que las sales en su recorrido por la central atraviesan numerosas,
válvulas, bombas e intersecciones de tuberías, que resultan más difíciles de aislar debido
a su forma no uniforme, lo cual puede producir punto fríos donde por el alto punto de
solidificación de la sal se producen acumulaciones solidas que conducen a
obstrucciones de la red de tubería, y consecuentes roturas.
Para una muestra aleatoria de 19 horas (eje X) con distinta irradiación Gs (W/m2)
(representada mediante la curva negra). La variación de la irradiación afecta al caudal
de sales fundidas que circula por el interior de las tuberías, como se analizó en el
aparatado 4.6 de este proyecto. El diagrama de barras muestra el valor de la temperatura
de salida del tramo más desfavorable en función del espesor de aislamiento empleado.
La barra de color rojo representa el caso sin aislamiento, en azul el espesor de
aislamiento es de 10 cm, en verde 50 cm de aislamiento y por último en morado un
aislamiento de 80cm. El tramo más desfavorable depende del tipo de red de tuberías
analizada, en este proyecto se distinguen dos tipos de redes de distribución en función
de la temperatura de las sales que circulan por su interior, para la red caliente el tramo
más desfavorable es el tramo 5 y para la red fría es el tramo 1.
RED CALIENTE
La red caliente recoge las sales fundidas a la salida de los lazo a 550ºC y las
conduce hasta la isla de potencia.
570 1000
[ºC]
900
550 800
G_S (W/m^2)
700
Tª salida TRAMO 5
530
600
510 500
400
490 300
470 200
100
450 0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
Muestra horario del año modelo
Tabla 21. Influencia del espesor de aislamiento sobre la temperatura de salida con
distinta irradiación – Red Caliente
110
Optimización del espesor de aislamiento
RED FRÍA
Conduce las sales fundidas desde el tanque frío (290ºC) hasta entrada de los lazos.
300 1000
290 900
800
Tª salida TRAMO 1
280
700
G_s (W/m^2)
270 600
260 500
250 400
300
240
200
230 100
220 0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
Muestra horaria del año modelo
Tabla 22. Influencia del espesor de aislamiento sobre la temperatura de salida con
distinta irradiación – Red Fría
La temperatura de salida del campo solar influye sobre el rendimiento del ciclo de
potencia. Una mayor temperatura de salida implica un mayor rendimiento del ciclo de
potencia. El aumento de la radiación solar, durante la operación de la planta, aumenta el
caudal másico de sales fundidas, lo que implica una mayor temperatura de salida. La
variación del caudal de sales en función de la radiación solar se ha detallado en el
apartado 4.6 de este `proyecto.
Observando la diferencia entre la barra roja, caso sin aislamiento, y la barra morada,
caso con 10cm de aislamiento, se llega a la conclusión de que el uso de un mínimo
espesor de aislamiento del orden de 10 cm implica una mejora significativa, que
uniformiza la temperatura de salida del campo ante variaciones en el régimen de
operación (variación del caudal).
150 /
0,23 /
Las sales del tanque frio a una temperatura de 290ºC se introducen en el campo solar
y recorren toda la red de tuberías con un caudal por lazo de 0,23 kg/s. Se ha realizado un
análisis preliminar de la temperatura de la sal a la salida de los distintos tramos, tanto de
la red fría como de la red caliente, considerando una disminución de temperatura en el
interior del lazo de 10ºC
Temperatura de salida por Tramo
TRAMO TRAMO TRAMO TRAMO TRAMO
Esp TRAMO TRAMO TRAMO TRAMO TRAMO
CALIENTE CALIENTE CALIENTE CALIENTE CALIENTE
cm FRIO 5 FRIO 4 FRIO 3 FRIO 2 FRIO 1
1 2 3 4 5
20 288,4 285,6 282 277,4 270 257,6 253,8 252,2 251,4 250,9
25 288,6 286,1 282,8 278,6 271,7 258,8 255,3 253,8 253,1 252,7
30 288,7 286,4 283,4 279,5 273 260 256,7 255,3 254,7 254,3
35 288,8 286,7 283,9 280,1 273,9 261 257,8 256,6 256 255,6
40 288,9 286,9 284,2 280,7 274,7 261,8 258,8 257,6 257 256,7
45 289 287 284,5 281,1 275,3 262,5 259,6 258,5 257,9 257,6
50 289 287,1 284,7 281,4 275,8 263,1 260,3 259,2 258,7 258,4
55 289,1 287,3 284,9 281,7 276,3 263,7 260,9 259,9 259,3 259,1
60 289,1 287,3 285,1 282 276,7 264,1 261,5 260,4 259,9 259,6
65 289,1 287,4 285,2 282,2 277 264,5 262 260,9 260,4 260,2
70 289,1 287,5 285,3 282,4 277,3 264,9 262,4 261,4 260,9 260,6
75 289,2 287,6 285,4 282,6 277,6 265,2 262,8 261,8 261,3 261,1
80 289,2 287,6 285,5 282,8 277,8 265,5 263,1 262,2 261,7 261,4
112
Optimización del espesor de aislamiento
300
290
280
Tªde salida (ºC)
270
260
250
240
230
20 25 30 35 40 45 50 55 60 65 70 75 80
Espesor de aislamiento (cm)
En la operación real de la planta durante las horas de irradiación baja o nula, para
prevenir situaciones de disminución de temperatura cercas al punto de congelación de la
sal fundida, se utilizan resistencias térmicas instaladas a lo largo de toda la red de
tuberías que mantienen la sal fundida a 290ºC en todo su recorrido. Mantener las
tuberías de acero calientes con una temperatura de 290ºC recude las cargas térmicas y
dilataciones en la red de tuberías. Pero los consumos propios de las resistencias térmicas
disminuyen la eficiencia de la central y ponen en entredicho su viabilidad.
113
Optimización del espesor de aislamiento
7. ESTUDIO ECONÓMICO
e e
El coste de las pérdidas se sucede cada año, por lo que para calcular el valor neto
actual del gasto durante N años de vida de la instalación, con una inflación del coste de
la energía y un coste de oportunidad del dinero (interés que un banco da al invertir el
dinero) se debe multiplicar las pérdidas por el Coeficiente del VAN
114
Optimización del espesor de aislamiento
Proceso de cálculo para determinar el Coeficiente del VAN (Valor Actualizado Neto
de la inversión)
1
1
1 0,01
1 0,02
r=coste de oportunidad del dinero (interés que el banco nos hubiese dado por
invertir el dinero)
q n 3608 ; k 9
Siendo “n” las horas de operación de la planta en las que la irradiación solar del año
modelo es mayor de 150W/m2.
115
Optimización del espesor de aislamiento
Este coste se sucede cada año, por lo que para calcular el valor actual neto del gasto
durante N años se multiplica por el Coeficiente del VAN
Las pérdidas eléctricas, es decir, la electricidad que se deja de producir debido a las
pérdidas térmicas en la red de tuberías calientes es aproximadamente igual a 0,021
MWh para 80 cm de espesor y 0,057 MWh con 10 cm de espesor.
116
Optimización del espesor de aislamiento
Para la realización del estudio económico se han tomado las tarifas actualizadas de
la Orden IET/3586/2011, de 30 de diciembre 2011, y Orden IET/1491/2013 de 1 de
Agosto 2013 (que la prima de referencia)
6 €/ 60€/
8,7022 €/
117
Optimización del espesor de aislamiento
De manera que el valor de la energía para las dos opciones comentadas, será el
siguiente:
Tarifa Regulada, única para todos los periodos programación, expresada en céntimos de
euro por kilovatio hora.
29,89 €/ 8,7022 €/ ∙ 6%
30,4121 €/
6 €/ 28,1894 €/
8,7020 €/ ∙ 6%
34,71612 €/
Figura 47. Esquema de la retribución necesaria para dar viabilidad a cada tipo de
tecnologia renovable
118
Optimización del espesor de aislamiento
Este proyecto se basa en una central cilindro parabólica con almacenamiento de 7h,
por lo que según la tabla presenta 3950 horas equivalentes de referencia/año.
De las tablas del fabricante se determina la variación lineal del precio en función del
espesor empleado. [16]
119
Optimización del espesor de aislamiento
C €/ 3,77 1,233
C €/ 503,7
120
Optimización del espesor de aislamiento
La siguiente tabla recoge el coste en euros por metro cuadrado del los aislmientos.
espesor
CERÁMICA ROCA
(cm)
10 50,37 17,46
20 100,7 32,06
30 151,1 46,66
40 201,5 61,26
50 251,9 75,86
60 302,2 90,46
70 352,6 105,1
80 403 119,7
90 453,3 134,3
500
450
400
350
300
250
€
200
150
100
50
0
10 20 30 40 50 60 70 80 90
Espesor aislamiento(cm)
CERAMICA ROCA
El coste de la lana cerámica es mayor que la lana de roca. Para minimizar el coste
del aislamiento se emplea una combinación de ambos aislamiento.
121
Optimización del espesor de aislamiento
La primera tabla recoge el valor del coste de oportunidad, que hace referencia al
valor de la energía eléctrica que se deja de producir en la central debido a la pérdida de
temperatura en la red de tuberías. Considerando que la planta se acoge a la tarifa
regulada del régimen especial.
COSTE DE OPRTUNIDAD [€)
Esp 1AÑO 2AÑOS 3AÑOS 4AÑOS 5AÑOS 6AÑOS 7AÑOS 8AÑOS 9AÑOS 10AÑOS C. AISL
10 64820 131227 199263 268967 340378 413538 488493 565286 643961 724565 5738
20 43840 88753 134768 181910 230208 279690 330384 382320 435532 490045 12578
30 35715 72306 109794 148200 187548 227859 269159 311472 354822 399233 21090
40 31252 63269 96072 129678 164108 199382 235520 272545 310476 349338 31274
50 28369 57434 87211 117718 148973 180994 213799 247408 281842 317118 43132
60 26325 53296 80927 109236 138238 167951 198393 229581 261533 294269 56662
70 24784 50175 76189 102840 130145 158118 186777 216139 246220 277040 71865
80 23570 47719 72459 97805 123773 150377 177633 205557 234165 263476 88740
Coste total = Coste pérdidas + Coste aislamiento
espesor 1 AÑO 2 AÑOS 3 AÑOS 4 AÑOS 5 AÑOS 6 AÑOS 7 AÑOS 8 AÑOS 9 AÑOS 10 AÑOS
10 76296 142703 210739 280443 351854 425014 499969 576762 655437 736041
20 68996 113909 159924 207066 255364 304846 355540 407476 460688 515201
30 77895 114486 151974 190380 229728 270039 311339 353652 397002 441413
40 93800 125817 158620 192226 226656 261930 298068 335093 373024 411886
50 114633 143698 173475 203982 235237 267258 300063 333672 368106 403382
60 139649 166620 194251 222560 251562 281275 311717 342905 374857 407593
70 168514 193905 219919 246570 273875 301848 330507 359869 389950 420770
80 201050 225199 249939 275285 301253 327857 355113 383037 411645 440956
122
Optimización del espesor de aislamiento
Se observa que dado que los costes de las pérdidas térmicas en la red de tuberías se
suceden cada año, el coste de la energía eléctrica que se deja de producir se incrementa
a lo largo de los año, teniendo en consideración el Coeficiente del VAN, que actualiza
el valor del las pérdidas.
10 ñ
800000
700000
600000
500000
400000 C_perdidas
€
C_aislamiento
300000
C_TOTAL
200000
100000
0
10 20 30 40 50 60 70 80
Espesor de aislamiento (cm)
Tras a representación gráfica de los resultados se observa que el mínimo coste total
se produce con un espesor de aislamiento igual a 50 cm. Es decir, el óptimo espesor de
aislamiento para la red de tuberías es de 50 cm
Ó 50
123
Optimización del espesor de aislamiento
RED CALIENTE 50 cm
Tª salida lazo (ºC) 550
TRAMO TRAMO 1 TRAMO 2 TRAMO 3 TRAMO 4 TRAMO 5
Diámetro int. (mm) 67 83 108 108 162
Pérdidas (MW/año) 21,30 23,07 25,60 25,58 30,64
Max Tª sup. Ext. (ºC) 56,87 57,31 57,92 57,91 58,99
Max AT (ºC) 5,19 2,82 2,08 1,56 0,83
Max Tª salida (ºC) 549,4 549,7 549,6 549,5 549,5
Min Tª salida (ºC) 544,8 547,2 546 545,5 545,1
RED FRÍA (trasporta las sales fundidas desde el tanque frío de almacenamiento
hasta los lazos):
RED FRÍA 50 cm
Tª salida tanque frio (ºC) 290
124
Optimización del espesor de aislamiento
CRITERIO TÉRMICO 1:
60º
CRITERIO TÉRMICO 2:
∈ 260; 570 º
125
Optimización del espesor de aislamiento
La imagen adjunta muestra un esquema con los resultados térmicos obtenidos para
cada tramo de la red caliente, empleando el espesor óptimo de aislamiento de 50 cm.
126
Optimización del espesor de aislamiento
127
Optimización del espesor de aislamiento
ANALISIS ECONÓMICO
Espesor (cm) 50
Coste aislamiento (€) 86264
Coste lana cerámica (€) 27004
Coste lana de roca (€) 59260
Coste oportunidad (€) 317118
Coste total (€) 403382
128
Optimización del espesor de aislamiento
Figura 51. Fotografía tubería del campo solar sin aislamiento y tubería con aislamiento- Archimede
129
Optimización del espesor de aislamiento
130
Optimización del espesor de aislamiento
8. CONCLUSIONES
En la realización del estudio se han tomado hipótesis iniciales para las condiciones
ambientales y las condiciones de operación de la planta, tomado un año modelo en la
provincia de Sevilla y un régimen de carga obtenido de la documentación técnica de la
planta [9]. Las pérdidas térmicas del año modelo, se obtienen de un análisis horario
donde la radiación solar es superior a 150W/m2, debido a que en periodos prolongados
de radiación insuficiente sería necesario el uso de sistemas auxiliares de calentamiento
que no han sido objeto de análisis en el proyecto.
Para determinar el coste del aislamiento solo se han considerado los costes variables
en función de la cantidad de material y no los costes fijos asociados al transporte,
mantenimiento e instalación, que dependen del fabricante y no varían con el espesor de
aislamiento.
131
Optimización del espesor de aislamiento
En este proyecto es necesario emplear una capa inicial de la lana cerámica debido a
que la temperatura máxima de la lana de roca es de 500ºC y la red caliente presenta una
temperatura superficial de 550ºC sin aislamiento. La lana cerámica permite disminuir el
espesor de aislamiento necesario pero incrementa el coste del aislamiento como se
observa en la gráfica del apartado 7.3 del proyecto. Con la combinación de ambos
aislante obtenemos el menor coste con el menor espesor de aislamiento.
132
Optimización del espesor de aislamiento
9. BIBLIOGRAFÍA
[1] Universidad de Sevilla. Apuntes de Ingeniería Industrial asignatura de
Transmisión de Calor., Escuela Técnica Superior de Ingenieros (Septiembre 2013).
[10] ENEA. Archimede Solar Energy. Solar receiver tubes: The new frontiers
of molten salts. Disponible en:
http://www.enea.it/it/enea_informa/events/techitaly2012/Angelantoni_3ott2012.pdf
134
Optimización del espesor de aislamiento
[23] Kaynakli, O. Economic thermal insulation thickness for pipes and ducts:
A review study. Renewable and Sustainable Energy Reviews 2014; (30) 184–194
[26] Ucar, A; Balo,F. Determination of the energy savings and the optimum
insulation thickness in the four different insulated exterior walls. Renewable Energy
2010 (35) 88–94
[31] Ministerio de Agricultura, Alimentación y Medio Ambiente. Plan de
Energías Renovables 2011‐ 2020 (Vol. II): http://www.magrama.gob.es/es/cambio-
climatico/legislacion/documentacion/normativa-y-textos-legales/default.aspx [Última
visita: Junio 2014]
135
Optimización del espesor de aislamiento
[36] IDAE. OBJETIVOS DEL PLAN HASTA EL AÑO 2020. Plan de Energías
Renovables 2011- 2020.
136
Optimización del espesor de aislamiento
137