Diseño y Calculos Aislamiento AISLAM GT3 - 07 (IDAE)
Diseño y Calculos Aislamiento AISLAM GT3 - 07 (IDAE)
Diseño y Calculos Aislamiento AISLAM GT3 - 07 (IDAE)
Guía técnica para el diseño y cálculo del aislamiento térmico de conducciones, aparatos y
equipos
CONTENIDO
Esta publicación ha sido redactada por la Asociación Técnica Española de Climatización y
Refrigeración (ATECYR) para el Instituto para la Diversificación y Ahorro de la Energía (IDAE),
con el objetivo de promocionar la eficiencia en el uso final de la energía en los edificios
(contiene un CD con programa informático de cálculo).
........................................................
IDAE
Instituto para la Diversificación y Ahorro de la Energía
C/ Madera, 8
E-28004-Madrid
[email protected]
www.idae.es
Madrid, febrero de 2007
ÍNDICE
Presentación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
2 Transmisión de calor . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.1 Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
de transferencia de calor . . . . . . . . . . . . . . . . . . . . . . . . 16
al elemento no aislado . . . . . . . . . . . . . . . . . . . . . . . . 23
de intercambio de calor . . . . . . . . . . . . . . . . . . . . . . . . 23
5 Ejemplos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
en estado estacionario . . . . . . . . . . . . . . . . . . . . . . . . 34
a tubería desnuda . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
en un tramo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
P R E S E N TA C I Ó N
Así, para justificar que una instalación cumple las exigencias que se establecen en el RITE
podrá optarse por una de las siguientes opciones:
Por esta razón, el IDAE con el fin de facilitar a los agentes que participan en el diseño y di
mensionado, ejecución, mantenimiento e inspección de estas instalaciones, ha promovido
la elaboración de una serie de guías técnicas de ahorro y eficiencia energética en climatiza
ción, que desarrollen soluciones alternativas.
En concreto, la que nos ocupa, titulada “Guía técnica. Diseño y cálculo del aislamiento tér
mico de conducciones, aparatos y equipos”, junto con el programa informático AISLAM,
pretende ser un procedimiento alternativo, de acuerdo con lo establecido en la IT 1.2.4.2,
para el cálculo de los espesores de aislamiento.
NOTA: En este documento, todas las menciones al Reglamento de Instalaciones Térmicas en los Edificios se refieren al
último borrador disponible.
5
1
Objeto y
campo de aplicación
7
2
Transmisión
de calor
Conducción
9
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
dT
q / A = −k
Particularizada la anterior ecuación al caso de una capa
dn
cilíndrica en que sus superficies tengan una diferencia
[1] de temperaturas ΔT, tenemos:
⎛ 1 1 ⎞
2.2.1 Resistencias térmicas conductivas ⎜ − ⎟
= ⎝ int ext ⎠
en caso de placas planas r r
Rcond_esférica
4π k
Particularizada la anterior ecuación al caso de una placa [7]
plana en que sus superficies tengan una diferencia de
temperaturas ΔT, tenemos:
2.3 RESISTENCIAS TÉRMICAS POR CONVECCIÓN
q/ A =
ΔT
esp / k La ecuación que rige el intercambio de calor por convec
[2] ción es la conocida ecuación de Newton, la cual considera
que la densidad de flujo de calor por unidad de área es
Donde se define la resistencia térmica por conducción proporcional a la diferencia de temperaturas entre la su
de una placa plana como: perficie y la temperatura del fluido (líquido o sólido).
esp
Rcond _ plana = q / A = hconv (ΔT )
k
[3] [8]
10
Transmisión de calor
q/ A =
ΔT
1/ hconv
sente el fluido, así como de su estado, e incluso del
mismo gradiente de temperaturas (pared-fluido).
[9]
Respecto al movimiento se debe diferenciar entre movi
miento provocado (forzado) por un elemento (bomba, Donde se define la resistencia térmica por convección de
ventilador) o por el ambiente (velocidad viento), y movi una placa plana como:
miento natural (debido a la diferencia de temperaturas
dentro del fluido que a su vez provoca diferencia de den 1
sidades y por tanto desplazamiento). Rconv_plana =
hconv [10]
Respecto a su estado, cabe diferenciar el caso de gases,
líquidos o fluidos que en las condiciones de trabajo pre
senten cambios de fases (tuberías bifásicas).
2.3.2 Resistencias térmicas convectivas
Como órdenes de magnitud se pueden señalar: en caso de capa cilíndrica
q/ H =
ΔT
• Con movimiento forzado del orden de varias
1/(2 π r hconv )
decenas (10-100 W/m2 K). [11]
Esta gran variedad de coeficientes de convección hace Particularizada la anterior ecuación al caso de una capa
que el comportamiento al paso de calor en el caso de ga esférica en que tengamos una diferencia de temperatu
ses sea muy diferente con respecto a los demás. En ras ΔT entre la superficie y el fluido:
otras palabras, la resistencia térmica que ofrece un líqui
q=
ΔT
1/(π r 2 hconv )
do o un fluido en cambio de fase es despreciable frente
a la que ofrece un gas. En el apéndice I se dan las corre
laciones más usuales. [13]
11
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
Donde se define la resistencia térmica por convección de lo que una buena aproximación será suponerla igual a la
una capa esférica como: temperatura del aire. Por tanto, la expresión del flujo de
calor se puede expresar (linealizando la ecuación) como
1 un coeficiente de convección equivalente de radiación
Rconv_esférica = mediante:
π r 2 hconv
q / A = hrad (ΔT )
[14]
[17]
2.4 RESISTENCIAS TÉRMICAS POR RADIACIÓN Donde ΔT representa la diferencia de temperaturas en
tre la pared y el medio (aire) dicha diferencia se
La ecuación que rige el intercambio de calor por radia contabiliza en °C ya que es lo mismo que en Kelvin.
ción es la conocida ecuación de Stefan-Boltzman, la cual
considera que la densidad de flujo de calor por unidad Evidentemente el valor del coeficiente de convección
de área es proporcional a la diferencia a la cuarta poten equivalente en radiación será (suponiendo la tempera
cia de temperaturas (en Kelvin) entre superficies. tura del resto de superficies igual a la del aire):
q / A = C(TK sup1
4
− TK sup2
4
) Finalmente al haber expresado el intercambio radiante
[15] con la misma ecuación formal que la convección, la re
sistencia térmica para placa plana, capa cilíndrica o
En este caso la constante de proporcionalidad C presen esférica serán del mismo tipo que las convectivas, sin
ta diferentes valores en función de las propiedades más que sustituir el coeficiente de película en convec
radiantes de las superficies (coeficiente de emisión), y de ción por el correspondiente en radiación.
la forma del recinto (distancias y ángulos entre superfi
cies), siendo en general difícil su determinación exacta.
12
Transmisión de calor
1
Rconv −rad _ plana =
hconv _ rad [21]
1
Rconv_rad_cilíndrica =
2 π r hconv_rad [22]
1
Rconv_rad_esférica =
π r 2 hconv_rad [23]
Recordemos que en el caso de líquidos (o fluidos en cambio de fase) el calor intercambiado por radiación es despre
ciable y, por tanto, siguen siendo válidas las anteriores expresiones sin más que asignar un valor nulo al coeficiente
de película equivalente de radiación (hrad = 0).
Resaltemos que asumimos el estado estacionario, y no presencia de radiación de longitud de onda corta (exposición so
lar), por lo tanto, mediante el uso de estas expresiones no obtendremos el calor real transferido por los muros exteriores
de un edificio, ya que no consideramos ni inercia térmica, ni radiación solar.
De forma general se deberá contar la posibilidad de existencia de intercambio convectivo y radiante a ambas partes del elemento.
Las anteriores afirmaciones se resumen para placas planas en la constancia del flujo de calor por unidad de área, es decir:
q/ A=
ΔT int ΔT i ΔT ext
Rconv − rad _ plana,int R cond _ plana,i R conv −rad _ plana ,ext
= =
[24]
De donde se obtiene, simplemente sumando numeradores y denominadores (propiedad de las fracciones), y contabili
zando todas las capas):
T int − Text
q/ A =
Rconv −rad _ plana,int + ∑ Rcond _ plana ,i + R conv− rad _ plana ,ext
capas
material
[25]
13
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
Capa 1
Interior Exterior
Capa 3
Capa 2
Capa 1
Capa 2
Ti Te
Capa 3
Exterior Te
Para capa cilíndrica tenemos la constancia del flujo de calor por unidad de longitud, es decir:
De donde se obtiene, simplemente sumando numeradores y denominadores (propiedad de las fracciones), y contabili
zando todas las capas:
Tint − Text
q/H =
R conv−rad_cilíndrica,int + ∑ Rcond_cilíndrica,i + Rconv−rad_cilíndrica,ext
capas
material [28]
Las superficies interiores de la tubería están a la misma temperatura y, por tanto, aunque el fluido sea un gas el intercam
bio de calor por radiación es despreciable (o nulo). Consecuentemente, en el interior sólo se considerará el intercambio
Tint − T ext T −T
q/H = = int ext
⎛ ⎞ Rtotal,cilíndricas
ln⎜ r i+1 ⎟
1 ⎝ ri ⎠ 1
+ ∑ +
2π rint hconv,int
capas
2π k i 2π rext hconv_rad,ext
material [29]
14
Transmisión de calor
Tuberías
Exterior
Capa 2
Te
Capa 1
Interior
Di
Ti
Finalmente, para capa esférica tenemos la constancia del flujo de calor, es decir:
Tint − Text
q=
R conv−rad_esférica,int + ∑ Rcond_esférica,i + Rconv−rad_esférica,ext
capas
material
[31]
Igualmente al caso de tuberías, las superficies interiores de las esferas están a la misma temperatura y, por tanto, aunque el flui
do sea un gas el intercambio de calor por radiación es despreciable (o nulo). Consecuentemente, en el interior sólo se considerará
el intercambio de calor por convección.
Tint − T ext T −T
q = = int ext
⎛ 1 1 ⎞ Rtotal,esféricas
⎜ − ⎟
1 ⎝ ri ri+1 ⎠ 1
+ ∑ +
π rint2 hconv,int
capas
4π k i π rext
2
hconv_rad,ext
material [32]
15
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
Esferas
Exterior
Capa 2
Te
Capa 1
Interior
Di
Ti
• Paredes
Tint − T ext
q/ A = = U plana (T int − Text )
Rtotal , placas
[33]
1 1
U plana =
R total , planas 1 espi 1
=
• Tuberías
Tint − Text
q/H =
Rtotal,cilíndricas
(
= U cilíndrica 2π r Tint − Text )
[35]
1
1 2π r
U cilíndrica = =
2π r Rtotal,cilíndricas
⎛r ⎞
ln⎜ i+1 ⎟
1 ⎝ ri ⎠ 1
+ ∑ +
2π rint h conv,int capas 2π k i 2π rext hconv_rad,ext
material
[36]
16
Transmisión de calor
• Esferas
Tint − Text
q=
Rtotal,esféricas
(
= U esférica π r 2 Tint − Text )
[37
1
1 π r2
U esférica =
π r 2 Rtotal,esféricas
⎛1 1 ⎞
⎜ − ⎟
1 ⎝ ri ri+1 ⎠ 1
+ ∑ +
π rint2 hconv,int
capas
4π k i π rext
2
hconv_rad,ext
material
[38]
Observemos que tanto para tuberías como para esferas el coeficiente global de transferencia de calor no es constante, y
depende de la superficie de referencia (radio r).
Realmente, las cámaras de aire sólo son utilizables en configuraciones planas (raramente se utilizan en algún caso en tu
berías).
En dicha cámara de aire existe convección y radiación, tal como hemos visto. No obstante, en la práctica se le asigna a di
cha capa una cierta resistencia térmica global en función de su espesor y disposición. Es de destacar que la resistencia
térmica que se asigna no es proporcional al espesor, e incluso a partir de un cierto espesor el aumento del mismo no con
lleva una mayor resistencia, sino que disminuye (pensemos que en esos casos el movimiento del aire dentro de la cámara
es más libre).
En la siguiente tabla se facilitan los valores típicos (extraídos del código técnico de la edificación):
• Horizontales • Verticales
500 mm2 < Aapertura < 1500 mm2 por m2 500 mm2 < Aapertura < 1500 mm2
de superficie por m de longitud horizontal
17
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
Se considera la cámara de aire como un ambiente interior. Se toma el coeficiente de convección correspondiente y se des
precia el resto de cerramiento al exterior.
Se trata de la existencia de elementos no homogéneos dentro de las capas analizadas, y que en general favorecen el intercam
bio de calor por poseer mayor conductividad. Como caso típico podemos observar existencia de pilares en paredes, o
existencia de bridas o válvulas en tuberías.
Contabilizar estos elementos para el cálculo del flujo de calor es realmente complejo, ya que se trata de configuraciones en
donde existe flujo de calor bidimensional o tridimensional.
• Aumentar una cierta cantidad de calor (disminuyendo la resistencia térmica global o directamente con un por
centaje).
Todos estos procedimientos son muy difíciles de cuantificar y asignar correctamente a un caso, por lo que se debe recu
rrir a la experiencia del calculista.
En la práctica pueden existir dos tipos de resistencias térmicas aún no consideradas, las resistencias térmicas de con
tacto entre capas de diferente material y la resistencia térmica por depósitos de materiales en las superficies exteriores
o interiores.
Las primeras son debidas a que algunas combinaciones (chapa metálica sobre superficie constructiva rugosa, por ejem
plo) dejan pequeñas oquedades (en general de aire) que presentan resistencia térmica al paso de calor.
Las segundas son debidas a depósitos de material (caso típico de carbonatos cálcicos dentro de tuberías con circulación
de agua caliente).
Di
2.8.6 Caso de tuberías enterradas Ti
18
Transmisión de calor
terreno como una resistencia térmica del mismo, suponiendo el intercambio de calor entre la tubería y la temperatura del
suelo, es decir:
Tint − Tsuelo
q/H =
⎛ ⎞
ln⎜ ri+1 ⎟
1 ⎝ ri ⎠
+ ∑ + Rterreno
2π rint hconv,int
capas
2π k i
material [39]
La resistencia del terreno depende evidentemente de la profundidad a la que está enterrada la tubería y a la conductivi
dad del terreno, pudiéndose aproximar mediante:
1 ⎛2 p⎞ ⎛2 p
⎡ 2 ⎤
Rterreno ln ⎢⎜⎜
⎞
2π k terreno ⎢⎝ Dext ⎟⎠ ⎜D
= ⎟+ ⎜ ⎟ −1⎥
⎟ ⎥
⎝ ext ⎠
⎣ ⎦ [40]
⎛ Hr −10 ⎞
(
k terreno = rterreno
0,266
⎝
)
− 6,05 ⎜ 1− terreno
80
⎟
⎠
[41]
Donde la densidad del terreno está en el rango (1.200 y 1.800 kg/m3) y la humedad relativa del terreno (entre 5% y 30%).
Finalmente, la temperatura del terreno (a nivel del suelo) varía aproximadamente de forma cosenoidal, mediante la expresión:
⎛ 2π
Tsuelo = Tmedia ,suelo − Va cos⎜ (dia − dia 0 )⎞⎟
⎝ 365 ⎠
[42]
Donde el rango de valores posibles es:
Tmedia,suelo = Temperatura media, dependiendo de la zona climática de la localidad (Zona E 15 °C, Zona D 16 °C, Zona C 17 °C,
Zona B 18 °C y Zona A 19 °C).
dia0 = día valor mínimo (Zona E 34, Zona D 33, Zona C 32, Zona B
32 y Zona A 31).
c
Di
Ti
2.8.7 Caso de tuberías enterradas con caja o canalización
Como en el caso anterior se trata de flujo de calor bidimensional,
aunque se aproximará a un flujo unidimensional con una resisten
cia térmica de la canalización, así:
19
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
• Canalización con un material (en general arena El flujo de calor se calcula a través de cada pared, tomado
o similar) como placas planas, y con los coeficientes de convección
que se especifican para este caso en el apéndice I.
ln(requ /rext )
Rcanal = Como en tuberías, en el interior sólo se contabilizará el in
2π k [43] tercambio de calor por convección, ya que por radiación
es despreciable (las paredes interiores se encuentran a la
Donde el radio equivalente asignado será: misma temperatura).
1,07 c
requ =
2
Simplemente resaltar que el resultado final que nos intere
sa es el calor perdido por metro de conducto, expresando
[44] por tanto el resultado como:
En todos los casos se ha supuesto la existencia de una Para estimar la temperatura final a lo largo de un tramo
única tubería. Si existen varias y están a diferentes tem (longitud dada), o el flujo total de calor intercambiado
peraturas existen interacciones entre ellas, y las ver apartado 5.8.
anteriores expresiones no son válidas, consultar ASHRAE
2004 (Systems & Equipment cap. 11).
Depósito horizontal
Interior
Ti
Altura
Anchura
Longitud
20
Transmisión de calor
q q q
q= (π rint ) + (π r int ) + ( L)
2 2
[47]
Una muy buena aproximación se obtiene (debido al gran diámetro del depósito y
a la incertidumbre de cómo se produce la convección en la parte inferior del depó
sito) considerando que la resistencia térmica por metro cuadrado es la misma
independientemente de la superficie que consideremos (lateral, superior, infe
rior), e igual a la que se produce en el lateral del depósito. Así, la ec. [47] se
reconvierte en:
q
H cilindro
q= π rint2 + 2π rint L + π rint2
( )
2 π rint [48]
21
3
Estimación
del espesor de
aislante
Definir el espesor de aislante en una determinada insta 3.2 PARA PERDER UN PORCENTAJE DE CALOR
lación puede ser función de varios criterios técnicos. A CON RESPECTO AL ELEMENTO NO AISLADO
continuación se irán definiendo los mismos e indicando
la forma de estimar el espesor de aislamiento. Es quizás uno de los criterios más acertados, pues es
una forma relativa de establecer la bondad del sistema.
En todos los casos el procedimiento suele ser iterativo, ya No obstante, hay que destacar que en tuberías cambian
que los coeficientes de convección, y el de radiación de los coeficientes de convección exteriores, y esto para tu
penden en general de la temperatura de las superficies berías pequeñas es crítico. Por lo tanto, este criterio
(exterior e interior) y éstas a su vez dependen del flujo de puede parecer acertado para paredes y para tuberías
calor transferido que es función del espesor utilizado. Fi con un diámetro superior a 10 cm (en tuberías de menor
nalmente, en el caso de tuberías el coeficiente de diámetro este criterio no es adecuado).
convección puede depender del diámetro exterior de la tu
bería, y éste a su vez depende del espesor de aislamiento. El proceso de cálculo es idéntico al anterior, no obstan
te, el cálculo se debe realizar dos veces, uno sin la
existencia de aislamiento (probablemente con su nece
saria iteración), y otro partiendo del flujo de calor que
3.1 PARA INTERCAMBIAR UN FLUJO DE CALOR DADO finalmente se desea intercambiar, el cual se obtiene del
flujo de calor anterior, y el porcentaje asignado. Partien
Es el caso más sencillo, y el valor asignado a la densidad do de este valor se obtiene el espesor de aislamiento
de flujo de calor suele ser fijado por la experiencia. Es (como en el apartado 3.1).
una práctica habitual, aunque sus resultados pueden
ser muy alejados de valores óptimos de diseño.
Un caso práctico es asignar un valor de pérdidas en pa 3.3 PARA LIMITAR UNA RESISTENCIA TÉRMICA
redes de cámaras frigoríficas en función de la O UN COEFICIENTE GLOBAL DE INTERCAMBIO
temperatura interior de la misma (entre 6 y 7 W/m2 para DE CALOR
cámaras de congelación y entre 8 y 9 W/m2 para cáma
ras de refrigeración). La diferencia viene establecida Las ecuaciones a utilizar según la geometría analizada
fundamentalmente por el distinto coste de producción son (33,34,35), junto con (26,29,32), de las cuales se
de frío en función del nivel térmico requerido. deben despejar el espesor de aislamiento.
Para el caso de tuberías (limitar W/m) y esferas (limitar Señalar que en el caso de tuberías y esferas el procedi
W), señalar que no se suele utilizar este criterio. miento es necesariamente iterativo, ya que el radio
exterior aparece en dos términos de la correspondiente
Finalmente, especificar que en el caso de tuberías y esfe ecuación.
ras el procedimiento es necesariamente iterativo, ya que
el radio exterior aparece en dos términos de la corres Para capas planas, y como ejemplo en el Código Técnico
pondiente ecuación (resistencia de su capa y resistencia de la Edificación, vienen fijados unos valores máximos
convectiva-radiativa exterior). del coeficiente global de transferencia de calor (el inverso
23
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
en este caso corresponde con la resistencia térmica) para cerramientos en función de la zona climática y su tipo, el cual
se reproduce a continuación:
Zona Climática A B C D E
Muro 0,94 0,82 0,73 0,66 0,57
Suelo 0,53 0,52 0,50 0,49 0,48
Cubierta 0,50 0,45 0,41 0,38 0,35
Medianera 1,00 1,00 1,00 1,00 1,00
Particiones interiores 1,20 1,20 1,20 1,20 1,20
En tuberías y esferas no se suele utilizar este criterio. Señalar que el procedimiento es necesariamente iterativo, ya que
el radio exterior aparece en dos términos de la correspondiente ecuación.
En la práctica se trata de tuberías que transportan fluidos calientes (geometría cilíndrica), o depósitos que los contienen
(geometría cilíndrica en las paredes y plana en las superficies superior e inferior).
La estimación del necesario aislamiento se realiza igualando el flujo de calor total transferido al correspondiente entre la su
perficie que se quiere proteger (la exterior) y el ambiente exterior.
Tsup.ext − Text
• Caso de placas planas
Tint − Text
1 1 espi 1
=
1 r
=
ln ⎛⎜ i+1 ⎟⎞
2π rext hconv _ rad,ext 1 ⎝ ri ⎠ 1
2π rint h conv ,int capas 2π k i 2π rext hconv _ rad ,ext
+ ∑ +
material [50]
Observemos que en este caso es necesaria la iteración por estar el radio exterior en varios términos de la ecuación.
1 ⎛1 1 ⎞
=
24
Estimación del espesor de aislante
En el supuesto de querer mantener una temperatura su Para obtener el espesor económico expresamos todos
perficial interior se procedería de igual forma pero los costes en función del m2 de superficie de aisla
igualando al flujo de calor interior. (Realmente sólo tie miento. El coste de inversión en aislamiento CI (€/m2)
ne sentido en el caso de placas planas, ya que en se puede poner en función de dos términos, uno de
tuberías y esferas es imposible el contacto accidental en pendiente de la cantidad de material utilizado y otro
el interior de las mismas). del coste de instalación (independiente del espesor),
por tanto:
CI = aa + bb esp
3.5 EVITAR CONDENSACIONES SUPERFICIALES [52]
Este caso es semejante al anterior (se utilizan por tanto Si se conoce el precio del aislamiento para dos espeso
las mismas ecuaciones). Simplemente se trata de impo res dados es inmediato obtener dichas constantes:
CI1 − CI 2
ner una temperatura superficial que sea igual (o
bb =
superior) a la temperatura de rocío del ambiente, y con
esp1 − esp 2
ello que no se produzcan condensaciones superficiales.
CI1 − CI 2
(Ver apéndice II para estimar la temperatura de rocío).
aa = CI 1 − esp1
esp1 − esp2
La posibilidad de condensación superficial siempre se
da en el “lado caliente”, es decir, en paredes en la parte
más caliente, y en tuberías, únicamente si por ellas cir [53]
cula un fluido a temperatura inferior a la de rocío del aire
que la circunda exteriormente. El coste de explotación durante el primer año será en
función del flujo de calor transferido por m2, del tiempo
Es de señalar que en paredes la no existencia de con anual de funcionamiento t (s) y del coste de la energía
densaciones exteriores no garantiza que no puedan térmica producida s (€/J)
q
existir condensaciones dentro de la misma. En el su
t s
A
puesto de que las superficies exteriores sean
impermeables al paso de vapor sí que se garantiza la
no condensación interior. Esta situación se da en pa [54]
neles aislantes en cámaras frigoríficas.
Este coste se sucede cada año, por lo que para calcu
En tuberías por las que circula un fluido a temperatura lar el valor actual neto del gasto durante N años de
inferior a la de rocío del ambiente, la capa exterior siem vida de la instalación, con una inflación del combusti
pre debe ser impermeable al paso de vapor de agua, por ble utilizado i (%) y un coste de oportunidad del
lo que la presencia de dicha capa imposibilita el paso de dinero d (%), (interés que un banco nos hubiera dado
vapor y por tanto las posibles condensaciones dentro de por invertir ese dinero), se debe multiplicar por el VAN
la misma. (Valor Actual Neto).
⎛ 1 +i ⎞
N
⎟ −1
1+d ⎠
⎜
VAN = ⎝
3.6 EN FUNCIÓN DEL ESPESOR ECONÓMICO
⎛ 1 +i ⎞
⎟ − 1
Si i = d VAN = N
⎝1 + d ⎠
En este supuesto se trata de determinar aquel espesor ⎜
que minimice el coste total de la instalación teniendo en [55]
25
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
⎛ ⎞
VAN s k aisl 1 esp i 1
⎜ ⎟
especon = t (Text − Tint )− k aisl ⎜
bb ⎜ hconv _ rad ,int capas k i hconv _ rad ,ext
+ ∑ + ⎟
⎟
⎜ excepto ⎟
⎝ aislante ⎠ [57]
En la anterior expresión (si se calcula con precisión) es igualmente necesaria la iteración, ya que los coeficientes de
convección-radiación interior y exterior dependen del salto de temperaturas entre las paredes y los medios, y este
salto es función del espesor de aislamiento considerado.
En correspondencia con la situación anterior expresamos todos los términos en función del metro lineal de tubería.
El coste de inversión en aislamiento CI (€/m) se puede expresar en dos términos, uno dependiente de la cantidad de mate
rial utilizado (función del perímetro total) y otro en función del coste de instalación (independiente del espesor), por tanto:
CI = 2π (r + esp)(aa + bb esp)
[58]
Si se conoce el precio del aislamiento para dos espesores dados es inmediato obtener dichas constantes:
⎡ CI1 CI
(2π (r + esp1 ))− 2
(2π (r + esp2 ))⎥⎦
⎤
bb =
⎢⎣
esp1 − esp2
⎡ CI1 CI
2 r esp 2 r esp
− 2 ⎤
CI
aa = 1 esp1
⎢ ( π ( + 1 )) ( π ( + ))
2 ⎥
(2π (r + esp1 ))
⎣ ⎦
esp1 − esp 2
−
[59]
El coste de explotación del primer año (o funcionamiento) será en función del flujo de calor perdido por m, el tiempo anual
de funcionamiento t (s) y el coste de la energía térmica producida s (€/J).
q
t s
H [60]
Este coste se sucede cada año, por lo que utilizaremos el valor actual neto del gasto durante N años de vida de la insta
lación (VAN), como en el caso anterior, para obtener el coste de funcionamiento o explotación:
esp
Rcond _ plana =
k [61]
Y, finalmente, el coste total será CI+CF. Obteniendo el mínimo (derivando e igualando a cero) de dicha función de coste
se obtiene el espesor económico, aunque en este caso no es posible expresar el espesor económico de forma explícita en
función de las demás variables, y se hace necesaria la iteración para obtener la solución.
26
Estimación del espesor de aislante
Análogamente a los anteriores supuestos expresamos todos los términos en función del flujo de calor de la esfera.
El coste de inversión en aislamiento CI (€) se puede expresar en dos términos, uno dependiente de la cantidad de material
utilizado (función del área superficial total) y otro en función del coste de instalación (independiente del espesor), por tanto:
En caso de conocerse el precio del aislamiento para dos espesores dados es inmediato obtener dichas constantes:
⎡ CI1 CI 2 ⎤
π (r + esp1 )
( 2
π (r + esp 2 ) ⎥⎦
) ( 2
)
−
bb =
⎢
⎣
esp1 − esp2
⎡CI 1 CI 2 ⎤
CI 1 π (r + esp1 ) π (r + esp 2 ) ⎦
2 − 2 ⎥
aa =
⎢
⎣ ( ) ( esp1
)
(π (r + esp ) )
1
2 −
esp1 − esp2
[63]
El coste de explotación del primer año (o funcionamiento) será en función del flujo de calor perdido, el tiempo anual de
funcionamiento t (s) y el coste de la energía térmica producida s (€/J)
tqs
[64]
Este coste se sucede cada año, por lo que utilizaremos el valor actual neto del gasto durante N años de vida de la
instalación (VAN), como en los supuestos anteriores, para obtener el coste de funcionamiento o explotación:
CF = t q s VAN
[65]
Y, finalmente, el coste total será CI+CF. Obteniendo el mínimo (derivando e igualando a cero) de dicha función de coste se
obtiene el espesor económico, aunque en este caso es necesaria la iteración para su determinación.
En tuberías, es interesante conocer el tiempo que tardará (sin movimiento de fluido) en congelarse el agua de su in
terior partiendo de una determinada temperatura inicial, o planteado de forma alternativa, qué espesor de
aislamiento debemos utilizar para que se congele un determinado porcentaje de agua en un determinado tiempo (por
ejemplo, 8 horas nocturnas) sin movimiento del fluido y en unas condiciones dadas de temperatura exterior.
Supuesta toda la tubería (materiales y fluido) a una misma temperatura, podemos afirmar que el calor intercambiado se
utilizará en disminuir su temperatura (si es superior a la temperatura de congelación, normalmente asumida en 0 °C) o
en congelar el agua interior (si estamos a 0 °C). Así, la ecuación diferencial que marca el proceso será:
27
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
• Si Tagua > 0 °C
T
agua −
Text
−(magua Cpagua +
mtubería Cptubería )dTagua =
dt
Rtotal,cilíndricas [66]
• Si Tagua = 0 °C
0 − Text
dmhieloCf agua =
Rtotal,cilíndricas
[67]
En la práctica se suele despreciar el calor necesario para enfriar el material de la tubería (en todo caso esta afirmación su
pone una posición conservadora, en realidad, el tiempo para alcanzar esa temperatura será algo mayor).
El proceso no se realiza en condiciones constantes de resistencia térmica del sistema, ya que de una parte la resistencia
térmica (exterior) varía en función de la temperatura del agua, así en el primer instante la temperatura superficial será su
perior a la que se alcance en los instantes finales (la temperatura interior va disminuyendo), y de otra parte cuando se
está formando hielo, éste presenta una cierta resistencia térmica al paso de calor, por lo que la misma variando en fun
ción del porcentaje de hielo formado.
Como la mayor parte del intercambio de calor se realizará en el proceso de cambio de estado (debido al alto calor de
cambio de estado), y despreciando la resistencia térmica que ofrece el hielo (posición conservadora), operaremos con
la resistencia térmica que ofrece el sistema (Rtotal,cilíndricas) cuando la temperatura del agua sea de 0 °C.
Integrando las anteriores ecuaciones (con la consideración de resistencia térmica constante) obtenemos:
πDint
2 ⎡T − T ⎤
t enfriamiento =
R
total,cilíndricas ρ agua Cpagua ln⎢ ext agua,0 ⎥
4
⎢⎣ Text ⎥⎦
[68]
Por πDint
2
ρ
hieloCf agua Rtotal,cilíndrica
t congelación =
100 4
−Text
[69]
Siendo:
tenfriamiento = tiempo que se tarda en alcanzar 0 °C desde una temperatura inicial de Tagua.
tcongelación = tiempo que tarda en congelar un porcentaje (Por) de agua contenida en una tubería cuando en el instante ini
cial se encuentra toda la tubería en estado líquido a 0 °C.
Evidentemente el tiempo total será la suma de ambos, y podemos obtener la resistencia térmica (y por tanto el espesor de
aislamiento) que es necesaria imponer para afirmar que tras un tiempo “ttotal” el porcentaje de agua congelada será “Por”.
t total
Rtotal,cilíndrica =
⎡ Por ⎤
⎢
πD 100
2 ρ
hieloCf agua ⎡ T −T ⎤ ⎥
⎢
int
+
ρ
agua Cpagua ln⎢ ext agua,0 ⎥⎥
4 ⎢ −Text ⎢⎣ Text ⎥⎦⎥
⎢⎣
⎥⎦
[70]
28
Estimación del espesor de aislante
Siendo:
3.8 EN FUNCIÓN DE PRESENTAR UNA DIFERENCIA DE TEMPERATURAS A LO LARGO DE UNA TUBERÍA O CONDUCTO
El calor intercambiado a lo largo de una tubería (o conducto) será utilizado por el fluido interior en modificar su tempera
tura. Si limitamos la máxima diferencia de temperatura del fluido estaremos limitando el máximo flujo de calor
intercambiado, y con ello el espesor de aislamiento a imponer.
.
Tfluido − Text
−m fluido Cp fluido dTfluido = dH
Rtotal,cilíndricas
[71]
Como en el caso anterior, una vez obtenida la resistencia térmica total es posible evaluar el espesor de aislamiento, como
en el caso 3.3 (recordemos que se trata de un proceso iterativo).
Evidentemente despejando de la anterior expresión la temperatura del fluido al final de una distancia H (Tfluido,sal)
obtenemos:
.
−H
πDint
2
ρ fluido v fluidoCp fluidoR R total,cilíndricas
(
Tfluido,sal = Text + Tfluido,ent − Text e ) 4
[73]
.
πD 2
q= ρ fluidov fluidoCp fluido (Tfluido,ent − Tfluido,sal )
int
4 [74]
Lógicamente si los conductos son rectangulares el área a considerar no se corresponde con un círculo (πD2/4 ), sino con
un rectángulo (altura x anchura).
29
4
Estudio de
condensaciones
interiores
El estudio de condensaciones interiores se realiza para El flujo de vapor que atraviesa una placa plana por me
placas planas y se aplica a los cerramientos en construc tro cuadrado depende de la diferencia de presiones
ción, tal y como se especifica en el Código Técnico. El parciales de vapor a ambos lados de la placa y de la re
análisis se puede realizar bajo dos variables: la presión sistencia al paso de vapor de la misma, la cual es función
parcial de vapor y la temperatura de rocío. de su resistividad al paso de vapor y de su espesor, así:
m vapor =
rvi espi RVi
En el primer caso se trata de comparar la presión parcial
de vapor existente en cada punto del interior de la pared =
31
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
Si no existe condensación el flujo de vapor será constan ec. [24]), y con dicha temperatura podemos obtener la
te a lo largo de todas las capas, y ésta será la condición presión parcial de vapor en saturación (con Hr=100%,
que utilizamos para obtener la presión de vapor en cada ver apéndice III).
punto de unión de materiales. Por último, dentro de un
material homogéneo la distribución es rectilínea. (La ex Como hemos señalado, si dicha presión parcial en satu
presión final adoptada es semejante a la ec. [24] obtenida ración es inferior a la presión parcial de vapor estimada
para transmisión de calor en placas planas). anteriormente condensará.
32
5
Ejemplos
1
Rconv _ plana =
5.1 PAREDES. ESPESOR DEL AISLAMIENTO
PARA TENER UN COEFICIENTE GLOBAL hconv
DE TRANSFERENCIA DE CALOR
Al tratarse de placas planas las resistencias térmicas de Finalmente podemos despejar el espesor de aislamiento
las diferentes capas se obtendrán según la ec. [3] y la [10]: requerido, ec. [3]:
33
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
Tint − Text 22 − 0
convección son idénticos al caso anterior, añadiendo la
q/ A = = 31,68 W/m2
Rtotal, plana 0,695
conductividad del enfoscado de cemento (1,4 W/m2 K).
Con dichos datos se obtienen las resistencias térmicas =
de las diferentes capas (utilizando las ec. [3] y [10].
q/ A=
ΔT int
Rconv −rad _ plana,int
Datos diferentes capas material =
Nombre Espesor Conductividad Resist. térmica
(mm) (W/m K) (m2 K/W) ΔT i ΔT ext
Rcond _ plana ,i R conv− rad _ plana, ext
Interior 0,130 = =
Enlucido yeso 15,00 0,300 0,050
Ladrillo hueco 65,00 0,490 0,133
Cámara de aire 50,00 0,18 Sus resultados se representan gráficamente en la figura
Enfoscado inferior.
15,00 1,4 0,011
de cemento
Ladrillo perforado 115,00 0,760 0,151
Exterior 0,040
16
Temperatura Tx (°C)
14
12,101
12
10
8 6,399
6,06
6
4
1,267
2 0
0
0 20 40 60 80 100 120 140 160 180 200 220 240 260
Distancia x (mm)
34
Ejemplos
En este caso debemos añadir la resistividad al paso de vapor de los diferentes materiales, y calcular la resistencia al paso
de vapor de las distintas capas de material de acuerdo con la ec. [75]
Rvapor , j = esp j rv j
Y su coeficiente global de transmisión de calor (inverso A partir de la presión parcial de vapor podemos calcular
resistencia térmica) es 0,58 W/m2 K. igualmente la temperatura de rocío (apéndice III, ec
[A.III.7]) para dichos puntos, obteniéndose:
Asimismo debemos obtener la presión parcial de vapor
(y su temperatura de rocío) a ambas partes del muro: Pv (Pa) Tr (°C) Ts (°C)
Interior 22
• Exterior
1587 13,89 20,34
Text = 0 °C Hrext = 80% (apéndice III)
> Pvs = 611 Pa
Text = 0 °C y ec. [A.III.5] —— Enlucido yeso
Hrext = 80% y ec. [A.III.6] —— > Pvext = 489 Pa 1444 12,44 19,70
Pvext = 489 Pa y ec. [A.III.7] ——> Trext = -2,69 °C Ladrillo Hueco
1133 8,81 18,00
• Interior
Lana de vidrio
Tint = 22 °C Hrext = 60% (apéndice III
análogamente) —— > Trint = 13,89 °C Pvint = 1587 Pa 1062 7,85 2,03
Ladrillo Perforado
Seguidamente deberemos ir obteniendo las diferentes 489 -2,69 0,51
presiones parciales en los puntos de unión de materia Exterior 0
les mediante la ec. [78]
RV i RV j ∑ RVi
= = ratura seca (obtenida por el mismo procedimiento que
en el ejemplo 5.2) y la temperatura de rocío tenemos:
capas
material
35
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
12
10 8,81
7,847
8
6
2,025
4
2
0 0,000
0
-2 -2,689
0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200 210 220
Distancia x (mm)
Y como observamos se cruzan, es decir, SÍ que existirán en dichas condiciones condensaciones dentro del muro, aun
que no condensa en las superficies, ni interior ni exterior.
5.4 PAREDES. CÁLCULO DEL ESPESOR DE AISLAMIENTO PARA TRANSFERIR UN FLUJO DE CALOR DADO
Se desea estimar el espesor de panel de poliuretano a utilizar en una cámara frigorífica para transferir en
estado estacionario 7 W/m2. La cámara se encuentra a -20 °C estando en el interior de una nave que se su
pone a 30 °C (condiciones de verano).
Considerar una altura de las paredes de 3 m y superficie no metálica (coef. Emisión 0,9).
Los coeficientes de convección_radiación a utilizar tanto en el interior como en el exterior, dependen de las temperaturas
superficiales (ver apéndice I, ec. [A.I.1] o [A.I.2], ec. [18] y [20]), y estas a su vez dependen del flujo de calor y del espesor
de aislamiento utilizado, por lo que es necesaria la iteración (se puede iniciar con un espesor de aislamiento cualquiera).
Condiciones Interiores
Tint = -20 °C
Altura: = 3 m
Coef.emisión: = 0,9
36
Ejemplos
= 3,34W / m 2 K
Condiciones Exteriores
Text = 30 °C
Altura: = 3 m
Coef. emisión: = 0,9
Coef. convección turbulenta 1,71 W/m2 K ec. [A.I.2]
= 5,66 W / m 2 K
Coeficiente convección_radiación exterior calculado (recinto cerrado): hext = 7,37 W/m2 K ec. [20]. La distribución de tem
peraturas en el interior (en este caso en las superficies) se calcula con la ec. [24].
Por lo tanto el resultado final obtenido es (recordamos que después de realizar varias iteraciones):
37
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
Recordemos que el interior de la cámara se encuentra a -20 °C, y está dentro de una nave por lo que se supone una
temperatura media ambiente a lo largo del año de 22 °C.
10 cm 41 €/m2
12,5 cm 47 €/m2
La función a minimizar es el coste total CT, suma del coste de instalación CI (ec. [52]) más el coste de explotación en eu
ros constantes CF (ec. [56]).
Una aproximación se puede obtener mediante la ec. [57]. No obstante, se deben realizar igualmente iteraciones, ya
que ésta depende de los coeficientes de convección que a su vez dependen de las temperaturas superficiales y és
tas dependen del espesor de aislamiento utilizado. Si se resuelve a mano es mejor imponer unos coeficientes de
convección_radiación, y luego comprobar que son del mismo orden que los finalmente existentes.
Desde un punto de vista exacto el procedimiento es análogo al ejemplo anterior; únicamente algo más tedioso por
que hay que calcularlo para diferentes espesores y finalmente elegir aquel que minimice la función coste total.
Condiciones Interiores
Tint = -20 °C
Altura: = 3 m
= 3,33W / m 2 K
38
Ejemplos
Condiciones Exteriores
Text = 22 °C
Altura: = 3 m
Coef. emisión: = 0,9
Coef. convección turbulenta ec. [A.I.2]
= 5,22 W / m 2 K
esp
Rcond _ plana =
k
1
Rconv _ plana =
hconv
Resistencia total = 0,193+6,458+0,145 = 6,796 m2 K/W
39
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
CI1 − CI 2 41− 47
bb = = 240
esp1 − esp2 0,1− 0,125
=
CI 1 − CI 2
aa = CI 1 − esp1 = 41− 240 x 0,1 = 17
esp1 − esp2
⎛ 1 +i ⎞ ⎛ 1+ 0,06 ⎞
N 20
⎟ −1 ⎜
⎝1 + d ⎠ ⎝ 1+ 0,04 ⎠
⎜ ⎟ −1
VAN = = 24,11
⎛1+i ⎞ ⎛ 1+ 0,06 ⎞
=
⎟ −1 ⎟ −1
1 d ⎝ 1+ 0,04 ⎠
⎜ ⎜
⎝ + ⎠
q 0,033
CF = t s VAN = 7920 x 3600 x 6,19 x x 24,11 = 39,04 ∈ / m
2
A 1000 x 3600
El espesor económico se comprueba que coincide con el obtenido mediante la ec. [57]
⎛ ⎞
VAN s k aisl 1 espi 1
⎜ ⎟
especon = t (T − Tint )− k aisl ⎜
bb ⎜ hconv _ rad,int capas k i hconv _ rad ,ext
+ ∑ ⎟=
ext +
⎟
⎜ excepto ⎟
⎝ aislante ⎠
0,033
24,11 x 0,024
= 7920 x3600 1000x3600 (22 − (−20)) − 0,024⎛⎜ 1 + 1 ⎟⎞ = 0,159 m
240 ⎝ 5,18 6,9 ⎠
40
Ejemplos
5.6 TUBERÍAS. ESPESOR PARA PERDER UN TANTO POR CIEN RESPECTO A LA TUBERÍA DESNUDA
Apéndices
Calcular el espesor de aislamiento de lana de vidrio que es necesario utilizar en una tubería horizontal de acero
DN40, por la que circula agua a 90 °C en un recinto cerrado que se encuentra a 25 °C, para que la pérdida de calor
sea el 10% de la que se produciría con la tubería desnuda.
En cualquier caso, el coeficiente de convección interior es muy grande, lo que produce una resistencia térmica despreciable.
En este ejemplo, en primer lugar, hay que calcular el flujo de calor que cedería esa tubería desnuda, para la cual se pue
de suponer que la temperatura exterior de la tubería coincide con la interior (ya que al ser de acero la resistencia térmica
es despreciable).
El diámetro interior de la tubería es de 41,9 mm, teniendo un espesor de 3,2 mm, lo que da un diámetro exterior de
(41,9+2x3,2) = 48,3 mm.
El coeficiente de convección exterior a utilizar viene dado por las ec. [A.I.10], ya que:
ΔT 90 − 25
hconv = 1,254 = 1,254 = 7,57W / m 2 K
D 0,0483
= 7,45W / m 2 K
Obteniéndose un coeficiente total de 15,02 W/m2 K, lo que produce una resistencia térmica de ec. [12]
1 1
Rconv_cilíndricas = = = 0,439 m K/W
2 π r hconv 2 π (48,310−3 /2)15,02
90 − 25
q/ H = = 148 W / m
ΔT
1/(2 π r hconv ) 1/(2 π 24,1510 −3 15,02)
=
Seguidamente debemos de calcular el espesor de aislamiento que impuesto en esa tubería hiciera perder ese flujo de ca
lor, teniendo en cuenta que el aislamiento añade una resistencia tal como la ec. [5], y modifica la resistencia convectiva
(de acuerdo con la ec. [A.I.10]), al modificarse la temperatura superficial.
Una vez resuelto el proceso iterativo la solución final es espesor igual a 43,8 mm.
41
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
29,02 − 25
hconv = 1,254 = 1,254 = 2,91W / m 2 K
ΔT
D 0,1361
= 5,52W / m 2 K
Obteniéndose un coeficiente convección_radiación exterior calculado (recinto cerrado): hext 8,44 W/m2 K, y finalmente
produce una resistencia térmica de ec. [12]
1 1
Rconv_cilíndricas = = = 0,277 m K/W
2 π r hconv 2 π (136,110−3 /2) 8,44
Tint − T ext T −T
q/H = = int ext
⎛r ⎞ Rtotal,cilíndricas
ln⎜ i+1 ⎟
1 ⎝ ri ⎠ 1
+ ∑ +
2π rint hconv,int
capas
2π k i 2π rext hconv_rad,ext
material
90 − 25 90 − 25
q / H = = 14,80 W / m
0 + 0,001+ 4,122 + 0,277 4,395
=
42
Ejemplos
5.7 TUBERÍAS. CÁLCULO ESPESOR PARA QUE NO EXISTAN CONDENSACIONES SOBRE UNA TUBERÍA
Calcular el espesor de aislamiento (de conductividad 0,03 W/m K) para una tubería vertical de acero DN40,
que es necesario utilizar para que no condense en la superficie externa, cuando por ella circula agua fría a 8
°C, estando el exterior a 25 °C, con una humedad relativa del 90%, una velocidad del aire de 3 m/s y un coef.
de emisión de 0,9 (pintada).
La temperatura exterior es de 25 °C, y con una humedad relativa del 90%, obtenemos (con el apéndice III) una tempera
tura de rocío de 23,3 °C.
El resultado obtenido (después de iterar) es de 6,3 mm, con el que pasamos a realizar los cálculos.
El coeficiente de convección interior se considera muy grande, por lo que su resistencia térmica es despreciable.
La resistencia térmica de una capa de aislante de 6,3 mm ofrece una resistencia térmica de ec. [4]
Condiciones Exteriores
Text = 25 °C
Velocidad: = 3 m/s
Coef.emisión: = 0,3
Coef. convección turbulenta (v D=3 0,0609 = 0,1827 > 0,00855 m2/s) ec. [A.1.14]
(D = 41,9+2*3,2+2*6,3=60,9 mm)
43
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
v0,9 30,9
hconv = 8,9 8,9 = 31,65 W / m 2 K
D 0,0609
0,1
= 0,1
= 5,36W / m 2 K
Coeficiente convección exterior calculado (al ambiente exterior): hext = 37,01 W/m2 K
1 1
Rconv_cilíndricas = = = 0,141 m K/W
2 π r hconv 2 π (60,910−3 /2) 37,01
25 − 8 17
q/ H = = 12,43W / m
0 + 0,001+ 1,230 + 0,141 1,371
=
Como comentario destacar que si la misma tubería estuviera en el interior de un recinto (sin una velocidad de aire impues
ta, pero con la misma temperatura), el espesor necesario sería 23,7 mm (mucho mayor, pues el coeficiente de
convección_radiación en este caso sería de 7,94 W/m2 K, lo cual haría una mayor resistencia térmica superficial y para
mantener la misma temperatura superficial se debería aumentar la resistencia interior, lo que originaría finalmente un
mayor espesor de aislamiento).
44
Ejemplos
5.8 TUBERÍAS. CÁLCULO ESPESOR AISLAMIENTO EN TUBERÍA ENTERRADA, PARA PERDER UN FLUJO DE CALOR DADO
Apéndices
Calcular el espesor de aislamiento de lana de vidrio que es necesario utilizar para perder 18 W/m en una tubería de ace
ro DN40 enterrada a 1 m de profundidad, por la que circula agua a 90 °C. Suponer una temperatura de suelo de 15 °C, y
que la tubería se asienta sobre una canalización de arena (k = 1W/m K) de 350 mm de lado.
El proceso, como siempre, es iterativo, por lo que únicamente comprobamos la solución obtenida, y que en esta ocasión
es 34,4 mm.
En primer lugar hay que estimar la resistencia térmica del aislante, ec. [4].
2π k 2π 1
=
Finalmente la resistencia del terreno, ec. [40], para lo cual debemos primeramente calcular la conductividad del mismo, ec. [41].
⎛ Hr − 10 ⎞ ⎛ 10 −10 ⎞
k terreno = ρ terreno
( 0,266
− 6,05 ⎜1 − terreno
) ⎟ = 1400
0, 266
− 6,05 ⎜1 −
( ⎟ = 0,82 W / m K
)
⎝ 80 ⎠ ⎝ 80 ⎠
Y
1 ⎛2 p⎞ ⎛2p
⎡ 2 ⎤
Rterreno ln ⎢⎜⎜
⎞
2π k terreno ⎢⎝ Dext ⎟⎠ ⎜D
= ⎟+ ⎜ ⎟ ⎥
⎟ −1⎥ =
⎝ ext ⎠
⎣ ⎦
1 2x1 2x1
2
ln ⎢⎜ ⎟ −1 ⎥ = 0,459 m K /W
⎡ ⎤
⎛ ⎞ ⎛ ⎞
2π 0,82 ⎢⎝ 2x0,18725 ⎠ ⎝ 2x0,18725 ⎠
= ⎟+ ⎜
⎥
⎣ ⎦
La resistencia total será:
90 −15
q/ H = = 18W / m K
4,169
45
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
89,99
Lana de vidrio 34,4 0,040 3,524
26,56
5.9 TUBERÍAS. CÁLCULO ESPESOR AISLAMIENTO PARA QUE CONGELE UN TANTO POR CIEN DEL AGUA
CONTENIDA EN UN DETERMINADO TIEMPO
Calcular el espesor de aislamiento de lana de vidrio que es necesario utilizar en una tubería horizontal de acero DN40,
en la que en el instante inicial se tiene agua a 20 °C y el ambiente está a -15 °C con una velocidad de aire de 3 m/s
(pintada, coef. emisión=0,9). Se desea que como máximo se congele el 30% en 8 horas.
Como siempre el proceso es iterativo, y analizamos la solución final, que es 16,2 mm.
La resistencia que debe ofrecer la tubería viene dada por la ec. [70], así:
t total
Rtotal,cilíndrica = =
⎡ Por ⎤
2 ⎢
πDint ρ hieloCf agua ⎡T − T ⎤⎥
⎢ 100 + ρ agua Cpagua ln⎢ ext agua,0 ⎥⎥
4 ⎢ −T ext ⎢⎣ Text ⎥⎦⎥
⎢⎣ ⎥⎦
8 x 3600
= = 2,16 m K /W
⎡ 30 ⎤
π (41,910−3 ) 2 ⎢100 920 x 333800 ⎡−15 − 20 ⎤⎥
⎢ +1000 x 4190ln⎢ ⎥⎥
4 ⎢ 15 ⎣ −15 ⎦⎥
⎢
⎣ ⎦⎥
Condiciones Exteriores
Text = -15 °C
Velocidad: = 3 m/s
Coef.emisión: = 0,9
Coef.convección turbulenta (ya que v D = 3 x 0,0807=0,2421 > 0,00855 m/s) ec. [A.1.14]
46
Ejemplos
v0,9 30,9
hconv = 8,9 8,9 = 30,77 W / m 2 K
D 0,0807
0,1
= 0,1
Apéndices
= 3,53W / m 2 K
Coeficiente convección exterior calculado (al ambiente exterior): hext = 34,3 W/m2 K
1 1
Rconv_cilíndricas = = = 0,115 m K/W
2 π r hconv 2 π (80,710−3 /2) 34,3
ln(rext /rint ) ln((80,7/2)/24,15)
Rcapa_cilíndrica = = = 2,042 m K/W
2π k 2 π 0,04
47
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
⎛ 1 1 ⎞ ⎛ 1 1 ⎞
⎜ − ⎟ ⎜ − ⎟
⎝ rint rext ⎠ ⎝ 4/2 4/2 + 0,04 ⎠
Rcond_esférica = = = 0,033 m 2K/W
4π k 4 π 0,024
1 1
Rconv_esférica = = = 0,005 m 2 K/W
π r hconv π (4/2 + 0,04) 16
2 2
de calor, ec. [32]
Y un flujo
Coef.convección ec [A.1.16]
v0,8 50,8
hconv = (3,76 − 0,00497 T ) 3,76 0,00497 x 16 = 15,06 W / m 2 K
D 0,54
0,2
= ( − ) 0, 2
48
Ejemplos
En el interior el intercambio de calor por radiación es despreciable, luego la resistencia térmica interior será 0,066 m2 K/W.
Apéndices
El aislamiento ofrecerá una resistencia térmica de:
esp 0,025
Rcond _ plana = = 0,625 m K /W
2
k 0,04
=
25 − 22,35
hconv = 1,174 = 1,174 = 1,70 W / m 2 K
ΔT
L 0,6
= 1,78 W / m 2 K
2a+2b+2esp=2*0,5+2*0,6+4*0,025=2,30 m2
0,979
Rlineal = = 0,426 W / m K
2,30
qlineal 21,14
qsup i.int = = 9,609 W / m
2
Y la temperatura del aire a la salida del conducto se obtendrá mediante la ec. [73], (con el área del conducto).
49
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
.
−H
(
Tfluido,sal = Text + Tfluido,ent − Text e ) anc alt ρ fluido v fluidoCp fluido R total,cilíndricas
=
−30
0,5*0,6*1,2x5x1000x0,426
= 25 + (16 − 25) e =16,346º C
.
q = anc alt ρ fluido v fluidoCp fluido (T fluido ,ent − T fluido ,sal )= 0,5x0,6x1,2x5x1000 (16,346 −16) = 622 W
El proceso, como siempre, es iterativo, ya que los coeficientes de convección a utilizar dependen de la distribución de
Resistencia térmica del material del depósito nulo (conductividad muy alta)
rext =1000/2+40=540 mm
rint = 1000/2=500 mm
Coef. convección laminar (D3�T = 1,083(24,73-20) = 5,96 < 10) ec. [A.1.8]
50
Ejemplos
24,73 − 20
hconv = 1,324 = 1,32 4 = 1,91W / m 2 K
ΔT
D 1,08
Apéndices
= 5,27 W / m 2 K
Presentando una coeficiente convección_radiación exterior de 7,18 W/m2 K y una resistencia térmica, ec. [12]
1 1
Rconv _ cilindricas = = 0,041 m K / W
2 π r hconv 2 π (1080 10 / 2) 7,18
= −3
q T int − Text 60 − 20
=115,18 W / m
H R total 0,347
= =
cilindro
q
H cilindro 115,18
q= π rint2 + 2π rint L + π rint2 = 2 π 0,5 + 2π 0,5 x 2 = 288 W
2
2 π rint
( 2 π 0,5
) ( )
Produciéndose una distribución de temperaturas en el primer instante:
51
Apéndices
hconv = 1,324
I.2.1 Interior de edificios ΔT
I.2.1.1 Verticales H [A.I.3]
El régimen de circulación viene dado en función del pa • H3�T > 10 m3 °C régimen turbulento
hconv = 1,523 ΔT
rámetro H3�T
Donde: [A.I.4]
H es la altura de la pared (m)
�T es el valor absoluto de la diferencia de tempera • Coeficiente de convección superior cuando la pared
turas entre la pared y el aire (°C) fría está hacia arriba, o coeficiente de convección
inferior cuando la cara caliente esta hacia abajo (la
– H3�T ≤ 10 m3 °C régimen laminar gravedad no favorece el movimiento)
hconv = 1,324
ΔT • H3�T ≤ 10 m3 °C régimen laminar
H
hconv = 0,594
[A.I.1] ΔT
– H3�T > 10 m3 °C régimen turbulento H [A.I.5]
hconv = 1,743 ΔT
[A.I.2]
53
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
I.2.2.1 Verticales y Horizontales El régimen de circulación viene dado en función del pa
rámetro D3�T.
El régimen de circulación viene dado en función del pa
rámetro v H Donde:
D es el diámetro exterior de la tubería (m)
Donde: �T es el valor absoluto de la diferencia de tempera
H es la altura de la pared (o la anchura para superfi turas entre la pared y el aire (°C)
cies horizontales) (m)
v es la velocidad del aire (m/s) • D3�T ≤ 10 m3 °C régimen laminar
v
[A.I.11]
4
hconv = 5,765
H
I.3.1.3 Caso de conductos de aire
[A.I.7]
En este caso existen dos superficies planas verticales y
I.3 Tuberías (coeficiente de convección exterior dos horizontales, siendo de dimensiones relativas varia
con aire) bles en función de la relación altura/anchura y su
régimen de circulación. Una vez ponderados dichos coe
I.3.1 Interior de edificios ficientes, en la práctica se toma el valor medio que se
expresa a continuación.
I.3.1.1 Verticales
Donde:
El régimen de circulación viene dado en función del pa H es la anchura del conducto (m)
rámetro D3�T �T es el valor absoluto de la diferencia de tempera
turas entre la pared y el aire (°C)
Donde:
D es el diámetro exterior de la tubería (m) • H3�T ≤ 10 m3 °C régimen laminar
�T es el valor absoluto de la diferencia de tempera
hconv = 1,17 4
turas entre la pared y el aire (°C) ΔT
• D3�T ≤ 10 m3 °C régimen laminar H [A.I.12]
hconv = 1,324
ΔT I.3.2 Exterior de edificios
D [A.I.8] I.3.2.1 Verticales y Horizontales
• D3�T > 10 m3 °C régimen turbulento El régimen de circulación viene dado en función del pa
hconv = 1,743 ΔT
rámetro v D.
[A.I.9] Donde:
D es el diámetro exterior de la tubería (m)
v es la velocidad del aire (m/s)
54
Apéndices
hconv = + 3,14
D D [A.I.13] II.1 Introducción
• v D > 0,00855 m2/s régimen turbulento En los cálculos usados en el código técnico de la edifica
v 0,9
ción se han fijado unos coeficientes equivalentes de
hconv = 8,9
convección_radiación en función de su posición relativa
D 0,1
y tipo de cerramiento.
[A.I.14]
1057(1,352 + 0,019T )v
• Cara caliente hacia abajo o fría hacia arriba
0,8
hconv =
hint = 5,88 W/m2 K (resistencia térmica 0,17 m2 K/W)
v 0,8
Parte superior para suelos en refrigeración
D
0,2 Parte inferior para suelos en refrigeración
[A.I.16] Cara fría hacia abajo
Donde:
D es el diámetro interior de la tubería (m) (o diáme II.3 Exterior de edificios
tro hidráulico en caso de conducto rectangular
Dh = 2 anchura x altura / (anchura + altura)) • Verticales: (Paredes)
v es la velocidad del aire (m/s) hint = 25 W/m2 K (resistencia térmica 0,04 m2 K/W)
T es la temperatura en °C
• Horizontales:
Expresión de ASHRAE Fundamentals 2005 • Cara caliente hacia arriba o fría hacia abajo
hint = 25,0 W/m2 K (resistencia térmica 0,04 m2 K/W)
55
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
La temperatura de rocío de un ambiente se define como aquella temperatura en que el vapor de agua existente se encon
traría en saturación. En otras palabras, si existe en dicho ambiente una superficie a dicha temperatura (o inferior)
aparecerá agua líquida condensada sobre la misma. Para obtenerla nos hará falta conocer la relación existente para el
agua entre la temperatura y la presión de saturación.
ASHRAE propone dos relaciones en función del valor de la temperatura (realmente define la presión parcial de vapor en
saturación en función de su temperatura, luego es necesaria la iteración si la variable de entrada es la presión parcial de
vapor), así:
• T < -0,0606 °C
⎛ − 5674,359
+ 6,3925247 − 9,677843 10−3 TK + 6,2215701 10 −7 TK 2
⎞
Pvs = exp⎜ TK
⎜ + ⎟
• T > -0,0606 °C
⎛ − 5800,2206
+ 1,3914993 − 4,8640239 10 −2 TK + 4,1764768 10 −5 TK 2 − ⎟
⎞
Pvs = exp⎜ TK
⎜
⎟
Donde
TK = Temperatura (K)
T = Temperatura (°C)
Pvs = Presión vapor en saturación (Pa)
Una aproximación más fácil de utilizar par el rango (0 °C a 200 °C) es:
2148,496−10,2858TK
Pvs =10 35,85−TK
[A.III.4]
56
Apéndices
En primer lugar con las anteriores expresiones obtenemos Las principales normas y manuales de consulta son:
la presión de saturación del agua a la temperatura seca.
Pvs = f (Ts)
ISO 7345
Aislamiento térmico. Magnitudes físicas y defini
[A.III.5] ciones
Hr
Pv = Pvs
100
ISO 12241
Aislamiento térmico para equipos de edificación
[A.III.6] e instalaciones industriales. Método de cálculo
Finalmente con dicha presión parcial obtenida y las an DTIE 12.01
teriores ecuaciones de ASHARE, calculamos la A. Viti. Cálculo del aislamiento térmico de con
temperatura de saturación a esa presión parcial. Tempe ductos y equipos. Atecyr
ratura, que como hemos definido, será la temperatura
de rocío (este cálculo deberá ser iterativo por la forma ASHRAE 2005
de las expresiones, aunque si utilizamos la expresión Fundamentals
aproximada dicha iteración no existe).
Tr = f (Pv)
ASHRAE 2004
Systems and Equipment Handbook
[A.III.7]
CÓDIGO TÉCNICO DE LA EDIFICACIÓN
R.D. 314/2006, de 17 de marzo
57
Guía técnica
Diseño y cálculo del aislamiento térmico de conducciones, aparatos y equipos
58
Apéndices
Subíndices
aire Aire
agua Agua
canal Canal
cilíndrica Capa cilíndrica
cond Conducción
conv Convección
conv_rad Convección + Radiación
econ Económico
equ Equivalente
ext Exterior
esférica Capa esférica
final Final
hielo Hielo
i Capa genérica de material
inicial Inicial
int Interior
media Media
plana Capa plana
rad Radiación
r Radio
suelo Suelo
sup Superficie
terreno Terreno
59
Títulos publicados de la serie ”Ahorro y Eficiencia Energética en Climatización”
Guía nº 1:
Guía técnica.
Mantenimiento de instalaciones térmicas
Guía nº 2:
Guía técnica.
Procedimientos para la determinación
del rendimiento energético de plantas
enfriadoras de agua y equipos autónomos
de tratamiento de aire
Guía nº 3:
Guía técnica.
Diseño y cálculo del aislamiento térmico
de conducciones, aparatos y equipos.
Incluye CD-ROM con programa AISLAM
9 788496 680081