] cos ML Sh ML + i sen ML Ch ML
Siendo la carga termica en funci6n del tiempo para ese sistema la aplicaci6n de la ec. K.31 enx=O: qt =
=
q(O,t) =
(1<.32)
~ (-T-T) -1m [_ _ _ _A._M---..;..---..;.._F_ _ _ _ e i(wt+
q.L(L, -r) T(L,t--r) d-r
q(L,t)
+
0
0
aplicando las transformadas de Laplace y sus propiedades (ec. K.73) podemos expresar el anterior conjlmto de ecuaciones como: h(O,s) = -r q,o(O,s) u(O,s)
+
-r q,L(O,s) u(L,s) (K.90)
h(L,s) = -r q,o(L,s) u(O,s)
+
-r q,L(L,s) u(L,s)
siendo la nomenclatura utilizada : h(O,s)
Transformada de Laplace de Ia funcion flujo de calor en x=O
h(L,s)
Transformada de Laplace de la funcion flujo de calor en x=L
u(O,s)
Transformada de Laplace de la funcion temperatura en x=O
u(L,s)
Transformada de Laplace de Ia funcion temperatura en x=L
1: q.o(O,t)
Transfonnada de Laplace de la respuesta en flujo de calor en x=O ante una excitacion delta de Dirac en x=O.
'tqJ,(O,t)
Transformada de Laplace de la respuesta en flujo de calor en x=o ante una excitacion delta de Dirac en x=L.
'tq,o(L,t)
Transfonnada de Laplace de la respuesta en flujo de calor en x=L ante una excitacion delta de Dirac en x=O. Transfonnada de Laplace de Ia respuesta en flujo de calor en x=L ante una excitacion delta de Dirac en x=L.
500
Apendice K
T
T
o
t
t- - - - - -
O'
o
-1-. L' x
t
Fig. K.12. Pared unicapa sometida a distribuc;ones aleatorias en temperatura en ambas caras.
Expresi6n que en forma matricial resulta :
h(O,S)] [h(L,s) =
['t q,o(O,s) 't qJ,(O,S)] ("(O,S)] 't q,o(L,s) 't qJ,(L,s) u(L,s)
(K.91)
Recordando que la transformada de Laplace de Ia distribuci6n de temperaturas resultante para una excitaci6n delta en x = 0 venia dada por K.84 :
Sh (L - x)
~
s (1,"
(K.92)
U II-Dirac T,O(X,S)
:y aplicando Fourier:
q
= _
1 aT
ax
't
= -
1
au ax
La transformada de Laplace del flujo de calor en x=L ante una excitaci6n en temperaturas 6(t) en x=o sera: 501
Manual de climatizaci6n. Torno II
't
_ A
o(L,s)
q.
[a u
a-Dirac T,O(X,S)
ax
J.~
1
~~
-L
Sh
ffl
(K.93)
Y la transfonnada de Laplace del flujo de calor en x=O ante una excitaci6n en temperaturas ()(t) en x=O sera:
't
q.
o(O,s)
_ A[0 tI a-Dirac T,O(X,S) ox
1 0
~~ . J. ~ ~ f ~ :"l Sh
(K.94)
Amllogamente a la expresi6n K.84 se puede calcular la transfonnada de Laplace de Ia distribuci6n de temperaturas resultante ante una excitaci6n delta en x = L :
Sh
L [;1
~
ffl 1
(K.95)
Sh
La transfonnada de Laplace del flujo de calor en x=L ante una excitaci6n en ternperaturas Il(L) en x=L sera:
'l
"t q,rJL,s)
-
II.
[
0 U a-Dirac
.,lJx,8)
T
ox
'l
L = -
II.
~
~~
Ch
.
Sh
~~ s
1
(X "
r. ~ s 1 (K.96)
l «"J 502
Apendice K
Y la transfonnada de Laplace del flujo de calor en x=O ante una excitacian en temperaturas o(t) en x::::L sera:
(O,s)
1:
~
=
_
A
au. (Xs)i &-Dlrac T,L , [ ox L Sh
Por tanto el procedimiento basico seria :
~ H.l
(K.97)
1) Calcular las transfonnadas de Laplace de las funciones excitacian en temperatura a ambos lados del muro unicapa u(O,s)
u(L,s)
2) Calcular mediante Ia matriz representativa de Ia capa, las transfonnadas del flujo de calor en x = 0 y x = L (aplicacian de la ec. K.91 con la ayuda de las ec. K.93, K.94~ K.96 y K.97). h(D,s) h(L,s) 3) Obtener las antitransformadas de las anteriores expresiones, resultando e1 flujo de calor en funcian del ticmpo en la cara x=D y x=L. q(O,t)
q(L,t)
Como es evidente el problema radica fundamentalmente en calcular las antitransfonnadas de las expresiones resultantes.
K. 3. 3. Resoillcioll del problema mediante transformada de Laplace para lin mllro multicapa COli excitacioll ell temperatllras en las dos caras Hasta ahora hemos planteado la resoluci6n del problema para un muro con una cap a de material y con excitaciones de temperatura cualesquiera en las superficies de dicha capa. no obstante, en el problema real existiran varias capas de matcrial y el cambio de temperaturas se renlizara en los medios que bafian las paredes extema e intcnla del muro, y no en las superficies exterior e interior de la pared. 503
Manual de climatizaci6n. Tomo II
Operando eon Ia eeuaei6n matricial K.91 se puede cxpresar las trnsrormadas de la temperatura y el flujo de calor en una cara de la plaea, en funei6n de la otra, nsi :
U(O,S)] = [A(S) B(S)] [U(L,S)] [h(O,s) C(s) D(s) h(L,s)
(K.98)
y teniendo en cuenta K.93, K.94, K.96 y K.97, resulta : 't q,L(L,s)
A(s)
't q,o(L,s)
B(s)
1
+
't
q.
o(L,s)
(K.99)
C(s)
't q,L(L,s) 't q,L(O,s)
-
't q,o(O,s)
't q,o(L,s)
D(s)
+
't q. o(O,s)
't
q.
o(L,s)
Pudiendo conc111ir que Ia matriz de coefieientes A(s), B(s), C(s) YD(s) representa las caraeteristicas de Ia capa de material que posee un espesor L, una conductividad A. y una difusividad tcrmiea a". La capa rcprescntativa de la eonvecei6n interior y exterior se obtiene realizando Ia transformada de Laplace de las eeuaciones que determinan el fllljO de calor intcremnbiado, tomando los coefieientes de convecei6n eomo constantes, asi : . 504
Apendice K
q e = qpext = he e (Te - Tpm)
(K.lOO)
y con las propiedades de las transformadas de Laplace (ec.K.73) :
h ;Cs)
=
h(L,s)
=
he i [u(L,s) - u ,(s)]
h e(s) = h(O,s) = he e
[u e(s)
(K.lOl)
- u(O,s)]
siendo: hj(s) - Transformada de Laplace de la funcion flujo de calor que absorbe el aire interior (carga termica) he(s) - Transformada de Laplace de la funcion flujo de calor que entTa en In pared desde el aire exterior. Uj(s) - Transfonnada de Laplace de Ia funcion temperatura del aire que bafia la capa interior del muro. ue(s) - Transfonnada de Laplace de la funcion temperatura del aire que bana la capa exterior del muro. Las expresioncs K. 101 expresadas en forma matricial : 1 he e
[U(O,S)]
(K.I02)
h(O,s)
1 Analogamcnte, en el caso de tener una camara de aire caracterizada por su resistcncia al paso de calor (R j ), y por no poseer inercia termica, In podemos reprcsentar mediante la malriz earacteristica :
A(S) B(S)] [C(s) D(s)
=
[1
R/]
0
1
(K.103)
505
l\fanual de climatizaci6n. Torno II Observamos que en esta demostracion, en comparacion con el apartado K.2, existe que ha desaparecido, ello es debido a que estamos despej undo las condiciones exteriores en funcion de las interiores (al contrario que hicimos en el caso anterior). W1 signo (-)
Para estudiar el problema de Ia existencia de varias capas de material se abordanl el problema mediante superposicion de capas, igual como se hizo ante variaciones senoidales, asi:
I
(S)]
tl
I
1
.
(K.I04)
hl(s)
que se puede resumir en : u e(S)]
l
h e(S)
=
[AA (s) BB (S)] Itl t(S)] cc (s)
DD (s)
(K.lOS)
h I(S)
Interior
Exterior
u (s)
u.(s)
h (s)
h.(s)
e
I
e
I
1
2
n
Fig. K.13. Superposicion de varias capas.
506
Apendice K b02
EI detenl1inanle de Ia matriz de cada capa es Ia unidad (A(s) D(s) - B(s) C(s) = 1, ec. K.99), y por tanto el determinante de Ia matriz producto tambicn sent la unidad (AA(s) DD(s) - BB(s) CC(s) = I)} por 10 que podemos expresar la anterior matriz C0l110 :
h e(S)] = [ h ,(s)
[DD (s)/BB (s)
] [u e(O)S)]
- IIBB (s)
1IBB (s)
- AA (s)/BB (s)
(K.I06)
u ;(L,s)
EI procedimiento para calcular el flujo de calor en una de las caras del muro, en funcion de Ia evolucion de las temperaturas en los medios que banan Ia pared, seria semejante al visto para una sola cap a (ver apartado K.3.2.2), y que a continuacion se especifica : 1) Calcular las transformadas de Laplace de las funciones temperatura exterior (Te(t» y temperatura interior (Ti(t» en los medios que hanan las paredes del muro:
2) Ca1cular Ia matriz representativa de cada capa (capa gencrica "j") a partir de Ia ec. K.99.
3) Obtener el prodllcto matricial de todas las capas, nlanteniendo el orden indicado en K.l 04, estimandose Ia matriz de coeficientes representativa de dicho muro:
AA(s)
BB(s)
CC(s)
DD(s)
4) Realizar e1 producto matricial establecido en K. 106, con 10 que obtcndriamos las funcioncs:
5) Caicular Ia antitransfonl1ada de hies) y he(s), resllltando los flujos de calor en el interior y exterior respectivamente.
507
Manual de clirnatizaci6n. Torno II
Como es evidente el problema radica fundamentalmente en obtener la antitransformada del flujo de calor.
K.3.4.
Factores de resplIesta
En la aplicaci6n del metodo propuesto a nuestro caso particular, observamos una dificultad aiiadida, y es que la funci6n temperatura sol-aire no es una funci6n continua, sino que conocemos los valores de dicha variable a intervalos horarios, ya que tanto la temperatura ambiente, como el valor de la radiaci6n solar incidente, son conocidas con dicha cadencia de tiempo. (EI procedimiento seria amilogo con cualquier otro intervalo de tiempo). Podemos generar una funcian continua de la temperatura sol-aire, uniendo mediante rectas los valores conocidos en cada instante de tiempo (horas), producicndose una funcian continua tal y como indica la figura K.14.
Temperatura °C 60~------------------------------------------------~
,
50
-
-
-,
..
..
..
"."
..
..
..
t"
..
..
... ;
•
..
-
-
-;
•
..
..
-
...
,
-
...
..
I
","
..
sb-i.:air~
-- :- ---:- --- ~ ----: ----:-
,-
-
""
......... -. ...... •
..
..
I
-
Temperatura : 40
I
..
~
..
I
..
..
..
....
I
..
-
-
... -
..
-
..
.
-, - ............ . ,
•
30
20
.
......
10
'"
.
......
.
. . . ... .
.. .. .. .. .. . .. . I
,
,
.. .
,
..
.. ..
0~--~---+---+---4--~----~--+---+---~--~--~---4
o
2
4
6
8
10
12
14
16
18
20
Horas
Fig. K.14. Funcion continua de La temperatura sol-aire
508
22
24
Apendice K
La funci6n generada puede ser asimismo reconstruida mediante suma de fUl1ciones triangulo de altura unidad y anchura el doble del intervalo propuesto (es decir, 2h.), multiplicada por el valor de Ia temperatura sol-aire en cada instante. Ver figura K.I5 Temperatura °c 60~----------------------------------------------~
50
Ternp~rat~ra : SOl.:a ir-~
40
- - - - - -: - - - -
30
20 10
o
2
4
6
6
10
12
14
16
16
20
22
24
Horas
Fig. K.JS. Reconstruccion de IofunciOn continua TSQmediante slim a (Iefunciones tridngulo.
Dicha descomposici6n en suma de triangulos se puede realizar tanto con la temperatura sol-aire como con la temperatura interior, por 10 tanto: CD
(K.I07)
CD
(K.l08) siendo (n-I) at
~
t
~
n at
De acuerdo con los pasos establecidos en el apartado anterior para resolver este tipo de problemas, deberemos en primer lugar calcular las transformadas de Laplace de las excitaciones a ambos lados del muro, asi : 509
Manual de climatizaci6n. Torno II
(K.I09)
y tcniendo en cuenta 1a propiedad de que Ia transformada dc una suma es Ia suma de
trans formadas .:
alser los valores de Tsa(n-j) constantes, y teniendo en cuenta las propicdades de las transformadas: 00
U e(s)
=
L
Tsa(n - j) S£ [ftrilJnguIoJ.(t)]
(K.IIO)
j-O
Analogamente para la transfonnada de la temperatura interior tendriamos :
A continuaci6n debemos ca1cular para el muro en concreto estudiado la matriz caracteristica del muro AA(s)
CC(s)
BB(s)
DD(s)
de acuerdo con las ecuaciones K.99, K.I02, K.I03, y el producto matricial establecido en K.I04. Posterionnente, y de acuerdo con K.l 06 obtenemos las transformadas dc los flujos de calor en el interior y exterior mediante : h (s)
e
=
DD(s)
BB(s)
U e(s)
-
1 BB(s)
uls) (K.112)
h.(s) I
510
1 BB(s)
U e(s)
- AA (s) u .(s) BB(s) I
Ap~ndice
K
en la cual sustituycndo los valores de las transfonnadas de las funciones excitacion (ee. K. 110 YK.lll ) oblcnemos las transfonnadas de los flujos de calor buscados : CD
DD(S)E T (n-J) h (s) = •
BB (s)
e.
1
-
sa
~
rf ...,,, l .(/)]t' t"ungu OJ
CD
~ T.(n- j)~ rftrt " l .(/)] BB (s) ~ I t't ungu OJ
~
1 T (n- j)$i.[itr:'" 10 .(/)]BB (8) ~ sa 'ungu J
h .(8) = ,
j-O
AA(s)~ , - BB(s) ~ T,(n- })Si.
[flrldngulojt)]
Y aplicando antitransfonnadas tendremos el problema resuelto.
~
I) = Sf-1[DD (s) T (n _ ') q e( BB (s) ~ sa ]
B;(S)
q,(/) =
~
T/(n - j) Si.
~-l[ BB(s) 1 ~ ~
~[
(I)]
[flrldngulojt)]]
T (n - j) Sf. sa
.f
J trtangu1oJ'
[J:trl·..(.ungu10J.(/)] -
AA(s)~, ] ~ T,(n - J) Si. [flrldngulojt)]
- BB(s)
511
Manual de clirnatizacion. Torno II
Pnra cn1cular Ins nntitransformndns en primer lugnr utiliznremos In propiedad de que la antitransfomlada de una suma de funciones es la sum a de Ins antitransformadas :
. -I[ DD(s) BB (s)
. qit ) = L
Tsa(n - J) ~
j-O
co
+L
Ti(n - j)
~-
I[ - 1 BB(s)
j-O
~
~ [
ftrtanguloJ( t)
J] +
[ ftrltlngulojt) J] (K.113)
Lj..o
qlt)
Tsa(n - j)
00
T(n _ ') t J
+L j-O
siendo (n-I)
~t ~
~-I[ B;(S) ~ [ftrldngulojt)J] •
~_I[ - AA (S)
~ [ftrldngulojt)J)
BB(s)
t ~ n ~t
Y si somos capnces de obtener las antitransformadns de tres funciones respuestn del muro nnte entradas triangulo, tales como :
XG) - Respuesta en la superficie exterior en flujo de cnlor nnte unn excitnci6n triangulo en temperatura en In superficie exterior (instante j desde In perturbaci6n).
XU) =
g_l[DD (s) BB (s)
g [/,
.(/)]]
trlanguloJ
(K.114)
YG)" Respuesta en Ia superficie interior en flujo de cnlor ante una excitaci6n
triangulo en temperatura en la superficie exterior (instante j desde Ia perturbaci6n).
(K.llS)
512
Apendice K
ZG) - Respuesta en la superficie interior en flujo de calor ante una excitaci6n triangulo en temperatura en la superficie interior (instante j desde la perturbaci6n). ZU) =
~-1[ -BB(s) AA (s) ~ [1: .(t)]] ~angu~J
(K.116)
Podemos obtener los flujos de calor como : 01)
01)
E }-o
T sa(n - j) XU) -
E
Tt(n - j) YU)
}-o
(K.117) 01)
01)
qlt)
E }-o
Tsa(n - j) YU)
+
E
Tln - j) ZU)
}-o
siendo (n·l) .At s: t s: n .At Que para flujos de calor a intervalos horarios resulta : 01)
01)
qe(n)
E }-o
Tsa(n - j) XU) -
E
Tt(n - j) YU)
1-0
(K.llS) 01)
qln) =
E }-o
Tso(n - j) YU)
+
E
Tln - j) ZU)
}-o
Los valores de las funciones XG), YG) y ZG) son conocidos como factores de respuesta del muro considerado, y si observamos son completamente independientes de las excitaciones en temperatura a ambos lados del cerramiento, por 10 que dichos coeficientes definen el comportamiento del muro ante cualquier excitaci6n en temperaturas a ambos lados de la pared. El calculo de dichos coeficientes para un muro concreto se realiza a partir de su defmici6n, y teniendo en cuenta que Ia funci6n triangulo se puede descomponer a su vez en Ia suma de tres funciones ramp a, tal y como se indica en la figura K.16 513
Manual de climatizaci6n. Tomo II
A- --
I.
I
I
I
I
-1
0
-1
+1
0
+1
-1
0
+1
Fig. K.J6. Descomposicion de Ia funciOn tridngulo como suma de !res funciones I'ampa.
(K.119)
siendo:
t~(j-l)At
~
f1j(t)=O
t>(j-l)At
~
f1j(t)=(t+At)/L\t
funci6n rampa con existencia a partir del intervalo (j-l) At
j At
~
f 2j (t) = 0
t > j At
~
f2j (t) = t I At
t
~
funci6n rampa con existencia a partir del intervalo j At
514
Apendice K
0+1) ~t
~
f3J (t) = 0
t> 0+ 1) ~t
~
f3,j (t) = (t - ~t) /
t
~
~t
funci6n rampa con existencia a partir del intervalo G+ 1) ~ t Desarrollaremos el calculo para la funci6n X(t), siendo analogo el planteamiento para las otras dos. La expresi6n K.l 14 definia la funci6n X(t) en la que sustituyendo la funci6n triangulo por suma de las tres rampas, tenemos :
XU)
=
g_I[DD (s) BB(s)
g
[/'1 (t)
- 2
J
12 (t)
+
13 (t)]] J
J
(K.120)
La funci6n rampa de pendiente unidad se defme de forma generica por : f(t) = t siendo su transformada de Laplace 1/s2, por 10 que la expresi6n general, a la que dcbemos calcular su antitransfonnada sera de la forma indicada en la ec. K.121, la cual se puede descomponer en :
H(s) S2
(K.121)
G(s)
s+X
m
siendo Co' Cl' ~ constante y - xm las raices de la ecuaci6n B(8) 1) EI valor de Co se obtenia sin mas que multipIicar por expresi6n cuando s tiende a 0
d
s2
m
S+X
m
0
s2 y IIevar al limite la
~ C -:H(O) o
g(O)
2) EI valor de C 1 se obtenia sin mas derivar con respecto asia expresi6n anterior y llevar al limite cuando s tiende a 0 515
Manual de c1irnatizacion. Torno II
d [H (s) ]
OCI
---=O+C+Ld d8 G(s) ,1 m.l m
28 (s + x m) - 8 2 (s + X
)2 m
d[H
(8) ] =}C=--1 ds G(s) s~o
3) Finalmente el valor de dj se obtenia multiplicando Ia exprcsion general por (s+x) y IlcvtmdoJa al limite cuando s ticnde a Xj' (en dicho limite se debera aplicar L'Hopital pues da una indetcrminacion).
d L __ m_(8+x.)+d.
H(s) Co C1 ---(s+x) =-(8+X)=+-(s+x.)+ 8 2 G(s) J S2 J S J m.l
mOj S + Xm
J
J
luego:
.d .= l'1m )
s"'x· J
[H( s)(s +X.})] = l' 2
s G(s)
1m
H('S) d
-(s+x)+H(s) J
.
'S 1--------------s....x· dG( ) J
2s G(s)
+
s
S2
ds
Finalmcnte hacicndo uso de las relaciones de transformadas simples dc Ia pag. 495
(K.122)
siendo: - Grado dc1 polinomio H(s) < G(s) - O'(X ) = &O(s)/&s Is=xm
-xm
raices de G(s) = 0, siendo - m = 1,2,.... contador
516
IXm+d>\xml
Apendice K
Si denominamos
~ (t) = ~-1
DD (s) BB(s)
[
1
(K.123)
S2
la expresion K.120 se transforma en : ~ (t +
X(t)
at) - 2
~ (I) + ~ (t
- at)
(K.124)
at) - 2 ~(at)
(K.125)
que particularizada para intervalos de tiempo resulta : j
j =I j
~
2
XU)
X(O)
0
XCI) = ~(U+l)
= ~(2
= ~ (at)
at) - 2 ~U at)
+ ~(U-l)
at)
Ya que 141 funcion f2•j tiene existencia a partir de At Yf 3,j tiene existencia a partir de 21lt
Amllogamentc definimos 141 funcion :
1l (t) = $i-I [
1 s 2 BB(s)
1
(K.126)
y obtenemos :
YeO) = 1l(at)
j =0
j =1 j
~
2
YU)
Y(l) =
=
,,(2 at) - 2 1l (a t)
1l(U+l)·at) - 2 1lU at)
+
(K.127)
1l(U-l) at)
Por ultimo
J.l(t) =
~-1
[ 8
2
AA (s) ] BB(s)
(K.128)
517
Manual de climatizaci6n. Tomo II y obtcnemos :
j
=
0
j =1 ~
j
2
Z(l)
Z(O) =
=
Jl(a t)
Jl(2 at) - 2 Jl(at)
(K.129)
ZU) = Jl(U+l) at) - 2 JlU at) + J.l(U-l) at)
Las antcriorcs funciones de respuesta ante excitaciones triangulo, 0 los factores de respuesta quedan representados en la figura K.17 Por 10 tanto la fonna de calcular los factores de respuesta de un muro concreto sera: I) Obtener las matrices caracteristicas de cada capa de material, ec. K. 99, K.l 02 Y K.I03. 2) Realizar el producto matricial de acuerdo con K.I 04y cnleuInr In mntriz
representntiva del muro, ce. K.IOS 3) Caleular las nntitransformndas de las expresiones K.l23, K.l26 Y K.12S, ulilizando In eeuaci6n K.122 4) Obtener los diferentes factores de respuesta mediante las eeuneiones K.125, K.127 y K.129
Por ultimo debemos senalar que en general se utiliza como intervalo de tiempo unn hora (At = 1 ~.)
Temperatura
Muro
Temperatura
1°C /\ I
1
\
\
--~---
1°C
1\
1
-
Exterior
/
\ \
-~---
Interior
Fig. K.l7. Representacion grdfica de los factores de respuesta.
SIS
Apendice K
K.3.4.]. Propiedades de los factores de respuesta 1) Propiedad de la sumatoria Si las tcmperaturas interior y exterior son constantes, el fllUO de calor transmitido 10 sera en estado eSlacionario, y por tanto dicho flujo de calor sent igual a :
De otra parte al mantenerse constantes las temperaturas, estas pueden extraerse de la sumatoria de la ec. K.ll 8, con 10 que: co
q,
=
K (Tsa - T) = Tsa(n - j)
E }.o
YU)
+
Tln - j)
j ..
co
qe
=
K (Tla - T)I
Tsa(n - j)
E
E
ZU)
o co
XU) - Tln - j)
E
YU)
}..o
j-O
de donde se deduce : co
K =
GO
E }..o
E }..o
Y(j)
GO
XU)
-E }..o
ZU)
(K.130)
" La smna de los factores de respuesta es el coeficiente global de transmisi6n de calor"
2) Propiedad del ratio comun La soluci6n para los factores de respuesta era del tipo
s.e-1[
H(s) S2
G(s)
] _ ReO)
t
+
[~
G(O)&S
[H(S)] G(s)
+ 8.0
E m.l
H(x m) X 2 G '(X ) m
Los dos primeros tenninos quedan compensados cuando estimamos los factores de respuesta para j22 (ec. K.125, K.l26 y K.l27), ya que son smna de dos antitransformadas en los instantes 'J+ 1" y. 'j-l ", menos dos veces la antilransformada en "j"; y por tanto su suma total es cero paraj;;::2.
519
Manual de climatizacion. Torno II
Respecto a la sumatoria recordemos que en nuestro caso los xm cnm las soluciones del polinomio 88(5), siendo dichas raices sicmpre negativas, y cumpliendose IXm+ll>lxml, lucgo a partir de un cierto intervalo de tiempo "p", la soluci6n sera x1 unicamente funci6n de la primera raiz (e t ), (los detnas tcrminos de Ia sumatoria seran dcspreciables), y por 10 tanto los coeficientes senin de la forma:
Vj > p
YU) ~ Cte e Xl
j
pudiendo definirse un ratio comun para cualquier j > p
YU)
c
(K.131)
YU - 1)
Dicho ratio comun es el mismo para los factores de respuesta X, Y, y Z, ya que los valores de ~n son todos iguales, por ser el denominador de la ecuaci6n transformada de Laplace igual en todos los casos (8S(s)).
3) Expresion finita para el calculo del flujo de calor
Haciendo uso del ratio comun, la expresi6n del flujo de calor en un instante dado, se puedc escribir como:
qj{n) = Y{O) Tsa
Y(p) T 10(n-p) +
+
520
+
Y{I) Tso{n-l)
c Y(p) Tla(n-p-l)
Z{O) Tt{n)
Z(P) T;
+
+
+
Z(l) T t{n-l)
+
c Z(P) T;(n-p-l)
+
+
Y(2) Tso(n-2)
+ •..• +
c 2 Y(p) T sa(n-p-2)
+ ••••
Z(2) Tt {n-2) - .... + C2
Z(P) T;
-
Apendice K
y la del intervalo de tiempo anterior:
+
yep) Tsa
+
c Y(p) Tsa(n-p-2)
+
Z(O) T,(n-l)
Z(P) Tln-p-l)
+
+
Z(l) Tln-2)
+
c Z(P) Tln-p-2)
+
c
2
Y(p) Tsa(n-p-3)
+ .... -
Z(2) Tln-3) +
c 2 Z(P) T t(n-p-3) -
Por 10 que podemos expresar el flujo de calor que atraviesa la pared interior (en el exterior seria semejante), en un numero de terminos finitos sin mas que multiplicar la segunda expresion por el ratio comun y res tar ambas ecuaciones : p
qi(n)
=
c qln-l)
+
YeO) Tsa(n)
+
L
Tsa(n-j) [YU) - c YU-I)] -
/-1 p
+
Z (0) Tln)
+
L
Ti(n -j)
[ZU) - c Z U-1)]
/-1
siendo p el intervalo a partir del cual el cociente de los [actores de respuesta se puede considerar constante, e igual al ratio eomun. EI procedimiento puede ampliarse para recoger mas historia en flujo de calor qj(n-j) y menos en cllanto a temperaturas exterior e interior Tsa(n-j) y Tj(n-j).
K.3.4.2. Factores de respuesta de algunas paredes consideradas tipo A continuacion se facilitan los resultados del caleulo de los faetores de respuesta de tres tipos de muros representativos de la construccion espanola. Tambicn se han anadido otros datos como: - Peso .. Analisis en flujo estacionario · Coef. global de transmision de calor "K" .. Analisis con variacion senoidal de la temperatura exterior (24h) · Amortiguamiento 0 impedancia IBB I · Dcsfase en horas de la onda en flujo de calor 521
Manual de cHmatizaci6n. Torno II
Scm
Muro peso medio Fig. K.18. Muro peso medio
Peso = 307 kg/m2 Amilisis en flujo estacionario : · Coef. global de transmisi6n de calor K = 0,60154 W/m 2 K Amilisis en flujo de calor senoidal exterior (periodo de 24h): · Amortiguamiento IBBI = 5,400645 °C m2 /W · Desfase en horas = 8,696657 horas Coeficicnte de convecci6n exterior 16,67 W/m2 K Coeficiente de convecci6n interior 9,09 W/m2 K ESJ>C?sor m
C°-WJrui~dad
Ladrillo macizo (1 12 pic)
0,12
0,87
1800
1380
Enlbscado de cementa
0,015
1,4
2000
050
Policstireno expandido
0,04
0,033
25
837
Ladrillo hueco sencillo
0,04
0,49
1200
920
Enlucido yeso
0,015
0,3
800
920
Dtn&id~d "81 m
C1'~g~P
Exterior
Interior
522
Apendice K
Orden
Coeficienies Y
Coeficientes Z
Coeficientes X
1
0,000
-5.0386
10.855
2
0.88018E-3
1,7154
-2.846
3
0.76351E'-2
0.94904
-1.2021
4
0.19225E-l
0,61076
-0.78141
5
0,29196E-I
0,39396
-059939
6
0,35409E-l
0,25481
-0,50186
7
0,38279E-l
0,16540
-0.43797
8
0,38769E-1
0,10789
-038920
9
0,37741E-l
0,70844E-l
-034852
10
0,35829E-l
0,46926E-l
-0,31309
11
0,33458E-1
0,31443E-l
-0.28166
12
0,30901E-l
0,21382E-l
-0,25354
13
0,28324E-l
0,14810E-l
-0,22829
14
0,25827E-l
0,10488E-l
-0,20559
15
0,23466E-1
0,76183E-2
-0,18516
16
0,21266E-l
0,56898E-2
-0.16676
17
0,19239E-l
0,43734E-2
-0,15020
18
0,17383E-l
0,34575E-2
-0,13529
19
0,15692E-l
0,28054E-2
-0,12186
20
0,14156E-l
0,23289E-2
-0,10976
21
0,12765E-l
0,19708E-2
-0.98864E-l
22
0,11507E-l
0,16937E-2
-0.89051 E-l
23
0,10371E-l
0,14734E-2
-0,80212E-l
24
0,9341 1E-2
0, 12937E-2
-0,72250E-l
RELACION SOLUCIONES (ratio cornun)
c= 0,90074
523
Manual de climatizaci6n. Tomo II Coeflclente X
12r-------------------------------------------_________ •
10
. . . . . .: .. .... :. . . . . .
•
I
•
•
. . . . . . ' . . . . . . '. . . . . . .' . . . . . . ' . . . . . . , . . . . . .' ..
~
8
............: .......... :........... : ........... :......... -:.............: ........ . :........... : ...... .
6
. . . . . ., .... .. :. . . . . . ~ . . . . . . : . . . . . . :. . . . . . ~ . . . . . . :.. . ... ; .. . ...: ..
4
. . . . . .: . . . . . . :. .. ... : . . . . . . : . . . . . , ~ . . . . . ~ . . . . . . :. . . . . . : . . . . . .: . .
2
.........' .......................
;
..........
w' . . . . . . . . .
_:_
......
..
..
_: . . . . . . . . . .
_:_
....
......
:
..............
..
O~lr~:=='==~·~·--~:~:--~:~:--~:
V'" .: ...... ;...... : ......;...... ~ ..... : ..... .:..... ": ......:..
-2
4r-~~~~~~--~~~~----~--~~~_4~~~~~~,~
o
5
10
15
20
25
30
40
35
45
Orden Coeficiente Y
O,05r-----------------------------------------------------~ 0,04
........
... ' ............ '............
..
..
.
,
..
..
..
. . . . . - ....
....
.... ' . . . . . . . . . . .
~
..................
..
......
"
........
..
e' ..
..
..
0,03
.
.......................
0,02
.. ,'"'
..........
. ................... ..
",
. ..........
...... ;
.. -,
. ..
0,01
0~~4-~~--~~~~~~--~--~====c=~~ 5 o 10 15 20 30 40 45 25 35 Orden Coeficlente Z
2~----------------------------------------------------~
k···.·····.'······.·····:······,·················:··
1
O~----~~--~----~----~----~----~----~----~----~
.
.
.
,
-1
. . . . . ., . . . . . . " . . . . . , . . . . . ',' . .....' . . . . . " . . . . . . " ... .. ; . . . . . " ..
-2
........ . . . . . . . . . . . . . . . . . . . , . . . . . . - . . . . . . . . . . . . ,. . . . . . . . . . . - ... .
-3
. . . . . .: . . . . . . :. . . ... : . . . . . . : . . . . . . :. . . . . . .:. '"
-4
.... .......: ............ :...........
-5
. . . . . ., . .... ',' , , , .. , ' , .. . '.' . , . , "
~
.. :- . .... : . . . . . .' ..
........ .... :- .......... :......... . .: ........... :..... . . . : ........ ...: . . ' . , .. "
... , ',' . . . . . . ' . , .. ',' .
-6~~~~~~~~--~~~~~~~~--~~~~~~~~--~
o
5
·1 0
15
20
25
30
35
40
45
Orden
Fig. K.19. Representaeion graftea de los eoeJieientes 524
y,zX para el muro de peso medio
Apendice K Muro peso ligero
1,5 em
5 em
Fig. [(,20 Muro peso ligero
Peso = 200 kg/m2 Annlisis en flujo estacionario : · Coer. global de transmision de calor K = 0,37454 W/m2 K Annlisis en flujo de calor senoidal exterior (periodo de 24 h.) : · Amortiguamiento IBB I = 4,533943 °C m2 IW · Desfase en horas = 6,552024 horas Coeficiente de conveccion exterior 16,67 W/m2 K Coeficiente de conveccion interior 9,09 W/m 2 K
Espesor m
Exterior
Conductividad W/mK
Dcnsidad kglm3
Caloresp. J/kgK
0
Enfoscado de cemento
0,015
1,4
2000
1050
Ladrillo hueco doble
0,09
0,49
1200
920
Poliurctano expandido
0,05
0,023
40
837
Ladrillo hueco sencillo
0,04
0,49
1200
920
Enlucido yeso
0,015
0,3
800
920
Interior 525
Manual de clirnatizacion. Torno II III
Orden
Coeficicntcs Y
Cocficicntcs X
1
0,29565E-4
-5,0382
10,001
2
0,36765E-2
1,7282
-3,99R6
3
0,19647E-l
0,98292
-1,5334
..
O,35209E-l
0,65145
-O,97R78
5
0,42194E-l
0,43239
-0,72226
6
0,42732E-I
0,28739
-0,55082
7
O,39716E-I
0,19133
-0,42306
8
0,35093E-l
0,12761
-0,32551
9
O,30010E-l
0,85288E-I
-0,25063
10
0,25087E-I
0,57139E-I
-0,19306
11
0,2063 lE-1
O,38384E-l
-0,14876
12
0,16758E-l
0,25863E-l
-0,11465
13
O,13485E-l
O,17486E-I
-0,88390E-l
14
0,10770E-l
O,l1867E-l
-O,68156E-l
15
O,85513E-2
O,80872E-2
-O,52563E-I
16
0,67575E-2
O,55365E-2
-0,40543E-I
17
0,53 193E-2
O,38090E-2
-0,3 1276E-I
18
0,41739E-2
0,i-6344E-2
-0,24130E-l
19
O,32667E-2
0,18323E-l
-0, 18618E-I
20
O,25511E-2
O,12820E-2
-0,14366E-l
21
0,1 9886E-2
O,90246E-3
-O,11086E-l
22
O,15478E-2
O,63929E-3
-0,85555E-2
23
O,12032E-2
0,45575E-3
-O,66028E-2
24
0,93437E-3
O,32698E-3
-0,50961 E-2
RELACION SOLUCIONES (ratio comun)
526
Cocficicntcs Z
c = 0,77203
Apendice K
Muro pesado
em
Fig. K.21.
MUTO pesado
Peso = 494 kg/m2 Amilisis en flujo estacionario : · Coef. global de transmisi6n de calor K = 0.,43166 W/m 2 K Analisis en flujo de calor senoidal exterior (periodo de 24 h.) : · Amortiguamiento
IBB\ = 22,58752 °C m2 /W
· Dcsfase en horas
=
13,0.0.8 horas
Coeficicntc de convecci6n exterior 16,67 W/m2 K Cocficicnte de convecci6n interior 9,0.9 W/m2 K Espesor ConduCtividad W/mK m
Densidad kg/ml
Calor esp.
180.0.
1380.
40.
837
J/kgK
ExteriQr Ladrillo macizo (I pie)
0.,24
0.,87
Poliuretano expandido
0.,0.4
0.,0.23
Ladrillo hueco sencillo
0.,0.4
0.,49
120.0.
920.
Enlucido yeso
0.,0.15
0.,3
80.0.
920.
Interior
527
Manual de climatizaci6n. Torno II
Orden
Coeficicntes Y
Coeticientes X
1
0,000
-5,0383
10,855
2
0,59685E-5
1,7244
-2,8459
3
0,45990E-4
0,97185
-1,1960
4
0,44108E-3
0,63737
-0,75208
5
0,15585E-2
0,41852
-0,53669
6
0,33162E-2
0,27518
-0,41057
7
0,53458E-2
0,18124
-0,32890
8
0,73115E-2
0,11965
-0,27277
9
0,90112E-2
0,79234E-l
-0,23266
10
0,10362E-l
0,52700E-l
-0,20318
11
0,11355E-l
0,35263E-l
-0,18099
12
0,12026E-l
O,23790E-l
-0,16389
13
0,12423E-l
O,16230E-l
-0,15041
14
0,12600E-l
0,11238E-l
-0,13950
15
0,12606E-l
0,79341E-2
-0,13047
16
0,12483E-l
0,57387E-2
-0,12280
17
0,12265E-l
0,42728E-2
-0,11615
18
0,11981E-l
0,32874E-2
-0,11026
19
0,11652E-l
0,26188E-2
-0,10495
20
0,11293E-l
0,21595E-2
-0,10010
21
0,10918E-l
O,18389E-2
-0,95625E-l
22
0,10535E-l
0,16104E-2
-0,91449E-I
23
0,10151E-l
0,14432E-2
-0,87526E-l
24
0,~7696E-2
0,13172E-2
-0,83823E-l
RELACION SOLUCIONES (ratio comun)
528
Coeficientcs Z
c=0,95928
ApendiceL Comportamiento y ciasijicacion de los cerramientos verticales (paredes) y horizontales (techos) L.l. Tipos de cerramientos estudiados para paredes Como hemos visto (ec.7.20) Ia temperatura equivalente para un muro viene dada por la expresion : yu) (L.l) T'eI/(n) = T,a(n - j) -K
E GO
j-O
De ella se' deduce que los tipos de cerramientos en cuanto a la detenninacion de la temperatura seca equivalente exterior, se podnin clasificar por su caracteristica YG)IK. Se han estudiado 14 tipos demuros tipicos de construccion espanola, en los cuales se han mantenido constante los coeficientes de conveccion interior y exterior en los valores de 9,09 W/m2 °C y 16,67 W/m2 °C respectivamente, segUn norma ~BE-CT -79 (ref 6). Se facilitan los datos constructivos de los 14 tipos de muros estudiados, asf como : - Su peso (kg/m2) - Su comportamiento en estado estacionario
K Coef global de transmision de calor (W/m2 °C) - Su comportamiento en regimen variable senoidal de la temperatura exterior con un periodo de 24 h.
. IBB\ Amortiguamiento (m2 °CIW) . Desfase en horas. MURON.l
492 kg/m2 1,320 W/m 2 °C IBBI = 4,989 m2 °CIW Desfa'se -= 11,787 horas
Peso K
= =
24 em r
...
Fig. L.l
529
Manual de climatizaci6n. Tomo II ________ __________ ~
..
----'~--~--~----A~--------------
Espesor . Conductividad m W/m °C
Cnmposici6n Exterior Ladrillo macizo (1 pie) Camara de aire Ladrillo hueco sencillo Enlucido yeso . Interior
Densidad kg/m3
Calor. esp. J/kgOC
0,24 0,87 1800 0,05 Resist. termica = 0,18 m2 °CIW 0,04 0,49 1200 0,015 0,3 800
1380 920 920
MURON.2 =494 kglm2 = 0,432 W/m2 °C K = 22,587 m2 °CIW Desfase = 13,008 horas
Peso
IBsl
4 em 4 em
1.5 en
Fig. L2
Espesor
Composici6n
m
Exterior Ladrillo macizo (1 pie) Poliuretano expandido Ladrillo hueco sencillo Enlucido yeso Interior
0,24 0,04 0,04 0,015
Conductividari
Densidad
W/moC
kg/m3
0,87 0,023 0,49 0,3
1800 40 1200 800
Calor. esp. J/kgOC 1380 837 920 920
MURON.3
Peso K
= 276 kglm2
IBBI
= 1,429 m2
=
1,614 W/m2 °C
°CIW Desfase = 7,140 horas
.. ...
12 em
.. 5 em 4 em r
Fig.L.3
530
... 1.5 em
Apendice L
Espesor
Composici6n
m Exterior Ladrillo macizo ( 112 pie) Camara de nire Ladrillo hueco sencillo Enlucido yeso Interior
0,12 0,05 0,04 0,015
Conductividad W/m °C
Densidad kglm3
Calor. esp.
J/kgOC
0,87 1800 1380 2 Resist. termica = 0,18 m °CIW 0,49 1200 920 0,3 800 920
MURON.4 Peso
=307kg/m2
K
= ,602 W/m 2 °C = 5401 m2 °elW ,
\BB\
°
1 5 ern
Desfase"'= 8,697 horas Fig. L.4
Espesor m
Composici6n Exterior Ladrillo macizo (112 pie) Enfoscado cementa Poliestireno expandido Ladrillo hueco sencillo Enlucido yeso Interior
0,12 0,015 0,04 0,04 0,015
Conductividad W/moC
Densidad kg/m3
0,87 1,4 0,033 0,49 0,3
1800 2000 25 1200 800
o MURON.5 Peso
K
\BBI
= 282 kg/m2 = 1,538 W/m2 °C
:-?,X
IX
x ~
0 ·· ~ ·X . )1:
y
)t
12em
1380 1050 837 920 920
1;0 :
~
,Bj
1>< .
~ k~ ~,? -'"
J/kgOC
B
~.
0x~ ·
Calor. esp.
J::k 5 em 4 em
... 1 5 em
2
= 1,417 m °elW Desfase = 7,034 horas
Fig. L.5
531
Manual de clirnatization. Torno II
Espesor rn
Composition
Exterior Ladrillo perforado (1/2 pie) Enfoscado cementa Camara de aire Ladrillo hueco sencillo Enlucido yeso Interior
0,12 0,015 0,05 0,04 0,015
Conductividad W/m °C
Densidad kg/mJ
Calor. esp. J/kgOC
0,76 1,4
1600 1135 2000 1050 2 Resist. tennica = .0,18 m °CIW 0,49 1200 920 0,3 800 920
MURON.6
Peso K
IBBI
= 283,6 kg/m2 =
°,453 W/m
2
°C
2
°CIW Desfase = 8,352 horas = 6,241 m
Fig. L.6
Espesor
Composition
m Exterior Ladrillo perforado (112 pie) Enfoscado cemento Poliuretano expandido Ladrillo hueco sencillo Enlucido yeso Interior
0,12 0,015 0,04 0,04 0,015
Conductividad W/m °C . 0,76 1,4 0,023 0,49 0,3
Densidad kg/mJ
Calor. esp. J/kgOC
160()
1135 1050 837 920 920
~2000
40 1200 800
MURON.7
Peso K
= 200 kg/m2
IBB\
= 4534 m2 ,
°
= ,375 W/m2 °C
°CIW Desfase = 6,552 horas 532
em
Fig.L7
Apendice L Espesor
Composicion
m
Exterior Enfoscado cemento Ladrillo hueco doble Poliuretano expandido Ladrillo hueco sencillo Enlucido yeso Interior
0,015 0,09 0,05 0,04 0,015
Conductividad W/moC
Densidad
1,4 0,49 0,023 0,49 0,3
2000 1200 40 1200 800
Calor. esp. J/kgOC
kg/m3
1050 920 837 920 920
MURON.8
IBBI
= 234 kg/m2 = 1,356 W/m2 °C 2 = 1,258 . m °elW
Desfase
=
Peso K
5cm
6,298 horas
1,5cm
Fig. L.8
Composicion
Espesor m
Exterior Enfoscado cemento Ladrillo hueco Camara de aire Ladrillo hueco sencillo Enlucido yeso Interior
0,015 0,012 0,05 0,04 0,015
Conductividad W/moC
Densidad
Calor. esp. J/kgOC
kg/m3
1,4 2000 1050 0,49 1200 920 Resist. termica = 0,18 m2 CIW 0,49 1200 920 0,3 800 920 0
MURON.9
Peso K
= 236,5 kg/m2
0563 , W/m2 °C IBBI = 3,680 m2 °CfW Desfase = 7,550 horas =
Fig. L.9
533
Manual de climatizaci6n. Tomo II
Espesor m
Composici6n
Exterior Enfoscado cementa Ladrillo hueco Fibra de vidrio Ladrillo hueco sencillo Enlucido yeso Interior
0,015 0,012 0,05 0,04 0,015
Conductividad W/m °C
Densidad kg/m3
1,4 0,49 0,041 0,49 0,3
2000 1200 50 1200 800
Calor. esp. J/kgOC 1050 920 837 920 920
MURON.I0
Peso
= 322 kglm2
1,701W/m2 °C IBBI = 1,424 m2 °elW Desfase = 7,771 horas
K
=
15 em
1,5em
Fig. L.IO
Composici6n
Espesor m
Exterior Enfoscado cemento Bloque hueco hormigon Enlucido yeso Interior
0~015
0,20 0,015
Conductividad W/m °C 1,4 0,56 0,3
Densidad kg/ml 2000 1400 800
Calor. esp. J/kgOC 1050 1050 920
MURON.ll
Peso K
= 330 kg/m2
°,
820 W/m2 °C = 5,783 m2 °elW Desfase = 11,886 horas
IBBI
=
'
1,5 c
15em
Fig. L.II
534
Apendice L
Espesor m
Composici6n
Exterior Enfoscado cementa Bloque hormigon ligero Cfunara de aire Ladrillo hueco sencillo Enlucido yeso Interior
0,015 0,24 0,05 0,04 0,015
Conductividad W/m °C
Densidad kg/m3
Calor.-esp.
J/kgOC
1,4 0,33
2000 1050 1000 1050 2 Resist. termica = 0,18 m °C/W 920 0,49 1200 920 0,3 800
MURON.12
Peso K
=900kg/m2 = 1,662 W/m2 °C
IBBI
= 3,325 m2 °C/W
Desfase
= 10,038 horas Fig. L.12
Espesor m
Composici6n
Exterior Mamposteria de granito Cfunara de aire Ladrillo hueco sencillo Enlucido yeso Interior
0,3 0,05 0,04 0,015
Conductividad W/m °C
Densidad kg/m3
Calor. esp.
J/kgOC
2800 920 2 Resist. termica = 0,18 m °C/W 920 0,49 1200 920 0,3 800
2,5
MURON.13
Peso K
= 854 kg/m2
IBBI
=
2
Desfase
=
= 0339 W/m2 ,
°C 18991 m °C/W , 9,370 horas
Fig. L.13
535
Manual de climatizaci6n • Tomo II
Espesor
Composici6n
m
Exterior Mamposteria de granito Poliuretano Enlucido yeso Interior
0,3 0,06 0,015
Conductividad W/moC
Densidad kg/m3
2,5 0,023 0,3
2800 40 800
Calor. esp. .JJkgoC 920 837 920
MURON.14
kglm2
Peso
= 541
K
IBBI
= ,637 W/m2 °C = 10266 , m2 °CIW
Desfase
= 9,963
°
horas
1,5 em
Fig.L.14
Composici6n
Exterior Hormigon con mallazo Poliestireno expandido Honnigon en mas a Enlucido yeso Interior
Espesor m 0,07 0,04 0,15 0,015
Conductividad W/moC
Densidad kg/m3
1,63 0,033 1,6 0,3
240O· 25 2400 800
Calor. esp. JJkgoC 1050 837 1050 920
En In figura L.15 se han representado los valores Y(j)/K de los diferentes rnuros, especificandose su peso en kglrn2 . En ella se puede observar una cierta clasificaci6n en cuanto a su comportarniento termico en funcion del peso por metro cuadrado, aunque se pueden sefialar diferencias significativas en ciertos muros, como es el caso del muro nOll. Del analisis de la figura L~ 15 se puede clasificar los muros, respecto a su comportamiento tennico en : ligeros (peso sobre 200 kglm2), medios (peso sobre 300 kglm2) y pesados (500 kg/rn2 0 superiores). Habiendose tornado como representativos:
536
Apendice L
Muro Iigero
MuroN° 7
200 kg/m2
Muromedio
MuroN° 6
283,6 kg/m2
Muropesado
MuroN° 2
494 kg/m2
COEFICIENTE Y(i)/K
0,14
~eso Kg/m 2: (Numero m:uro)
0,12 282,(5) - " . - - 236".& (9} - - , -' - ' .- , , . -' , , - - - . : , . - . . , ! , 276 {3) . : 20 K 2,.: ----"'-"___ 322 (10) , 0 g/m, lIgera .
0,10
.
-
..
,
,
'.
,
,
.
,
..
.
. . -', , . . - , -'. , - . . , -', , , . - - :2 . , , , , : , - , , , . '- - . , - . '283.6 (6) , '283.6 Kg/ni Media ' ,
0,08
:
307 (4) :
:
:. - 900'(12):' . - . , .
0,06
:
:
-:494.Kg/m2~
330 (11):
:
, Pesada
'
:
~
- - . . . . :- . - . . . -
,
,
0,04
0.0: lJ~~E~!§i;;;-
":'-::'::-,': - : -:' :-'~::-'~'
°
6
12
18
24
30
-J'.
36
42
48
HORAS
Fig. L.IS. Caracterlstica YO)/K de diferentes muros
L.2. Tipos de cerramientos e~tudiados para tecltos Anfllogamente at caso de muros, la temperatura equivalente para un. techo viene dada por la ec. 7.20 : Q)
T (n) seq
=
L 1-0
T,a(n - j) yU)
K
(L.2)
De ella se deduce que los tipos de cerramientos en cuanto a la determinacion de la temperatura seca equivalente exterior, tambien se podnin clasificar por su caracteristica YG)IK. 537
Manual de clirnatizacion . Torno II
---------------------------------------------------~.~ Se han estudiado 14 tipos de techos tipicos de construcci6n espanola, en los cuales se han mantenido constante los coeficientes de convecci6n interior y exterior en los valores de 5,88 W/m2 °C y 20,00 W/m 2 °C respectivamente, segtin norma NBE-CT-79 (ref. 6). Se facilitan los datos constructivos de los 14 tipos de techos estudiados, asf como:
- Su comportamiento en estado estacionario . K Coef global de transmisi6n de calor (W1m2 ° C) - Su comportamiento en regimen variable senoidal de la temperatura exterior con un periodo de 24 h. . IBBI Amortiguamiento (m2 °CIW) . Desfase en horas.
TECHON.l
Peso K
= 538 kglm2 = I ,169 W/m2 °C
= 6,304 m2 °CIW IBBI Desfase = 12,086 horas
Fig. L.16
Composicion
Espesor m
Exterior Grav ill a Impermeabilizante Hormig. ligero de pendiente Bovedilla hormig6n normal Guamecido yeso Interior 538
Conductividad W/m °C
0,05 0,01 0,10 0,22+0,04 0,02
0,81 0,19 0,35 1,54 0,3
Densidad kg/ml
1700 1100 1000 1254 800
Calor. esp. J/kgOC
920 1680 1050 1050 920
Apendice L
TECHON.2 Peso K
IBBI
= 495 kgfm2
1 em
1,224 W/m2 °C = 4975 m2 °CIW ,
=
2:.:c~m~===~~~~meJ!=!:=~!I1~
Desfase = 11,155 horas
Fig. L.17
Composici6n
Espesor m
Exterior Gravilla Impermeabilizante Hormig. ligero de pendiente Bovedilla hormigan normal Guarnecido yeso Interior
Conductividad W/m °C
0,05
0,81
0,01 0,10
0,19 0,35
0,16+0,04 0,02
l,54 0,3
Densidad kg/m3 1700 1100 1000 1415 800
Calor. esp.
J/kgOC 920 1680 1050 1050 920
TECHON.3 Peso K
= 574 kglm2 = 1,126 W/m2 °C
IBBI
=
7696 m2 °CIW ,
Desfase
=
12,881 horas Fig. L.18
Composici6n
Espesor Conductividad m W/m °C
Exterior Gravilla Impermeabilizante Honnig. ligero de pendiente Bovedilla hormigon normal Guarnecido yeso Interior
0,81 0,05 ·0,19 0,01 0,35 0,10 0,27+0,04 1,54 0,3 0,02
Densidad kg/m3 1700 1100 1000 1168 800
Calor. esp.
J/kgOC 920 1680 1050 1050 920
539
Manual de climatizacion. Torno II
TECHON.4
Peso
= 453 kglm2
K
= 1,260 W/m2 °C
IBBI
= 4 ,708 m2
Desfase =
°CIW 10,417 horas Fig. L.19
Composition
Espesor Conductividad W/moC m
Exterior Impermeabilizante Honnig. ligero de pendiente Bovedilla hormigon normal Guarnecido yeso Interior
0,01 0,10 0,22+0,04 0,02
1100 1000 1254 800
Calor. esp. J/kgOC 1680 1050 1050 920
1cm
TECHON.5
K
= 553 kglm2 = 1,191 W/m2 °C
IBBI
=. 6,548 m2
Peso
0,19 0,35 1,54 0,3
Densidad kg/m3
°CIW Desfase = 12,639 horas
2cm
Fig. L.20
Composicion
Espesor Conductividad W/m °C m
Exterior Pavimento Impermeabilizante Hormig. Iigero de pendiente Bovedilla hormigon normal Guanlecido yeso Interior 540
0,05 0,01 0,10 0,22+0,04 0,02
1,1 0,19 0,35 1,54 0,3
Densidad kg/m3
Calor. esp. J/kgOC
2000 1100 1000 1254 800
1380 1680 1050 1050 920
Apendice L
TECHON.6 1cm
Peso K IBBI Desfase
=
510 kg/~2
2 = 1249W/m °C , 2 = 5, 166 m °CIW
=
2cm
11,709 horas Fig. L.2l
Composici6n
Espesor Conductividad m W/m °C
Exterior Pavimento Impermeabilizante Honnig. ligero de pendiente Bovedilla honnigon normal Guamecido yeso Interior
0,05 0,01 0,10 0,16+0,04 0,02
1,1 0,19 0,35 1,54 0,3
Densidad kg/mJ 2000 1100 1000 1415 800
Calor. esp.
J/kgOC 1380 1680 1050 1050 920
TECHON.7
Peso
= 589 kg/m2
K
= 1, 147 W/m2 °C
= 7,994 m2 °CIW IBBI Desfase = 13,433 horas Fig. L.22
Composici6n
Espesor Conductividad m W/m °C
Exterior Pavimento Impermeabilizante Hormig. ligero de pendiente Bovedilla honnigon normal Guamecido yeso Interior
0,05 0,01 0,10 0,27+0,04 0,02
1,1 0,19 0,35 1,54 0,3
Densidad kg/mJ 2000 1100 1000 1168 800
Calor. esp.
J/kgOC 1380 1680 1050 1050 920
541
Manual de climatizacion . Torno II
TECHON.8
Peso
= 443,5 kglm2 = 0,431 W/m 2
°C = 19220 °CIW , IBBI Desfase = 10,553 horas
K
m2
Fig.L.23
Composicion
Espesor Conductividad m W/m °C
Exterior Pavimento Poliestireno expand. Tipo IV Bovedilla honnig6n nonnal Guamecido yeso Interior
0,05 0,06 0,22+0,04 0,02
1,1 0,033 1,54 0,3
Densidad kglm3
2000 25 1254 800
Calor. esp.
J/kgOC 1380 837 1050 920
TECHON.9
Peso K IBBI
= 443,5 kglm2 = =
0,431 W/m 2 °C 10741 m2 °CIW ,
Fig.L.24
Composicion
Espesor
Conductividad
m Exterior Pavimento Bovedilla honnigon nonnal Poliestireno expand. Tipo IV Guamecido yeso Interior
542
0,05 0,22+0,04 0,06 0,02
Densidad kglm3
Calor. esp.
W/m °C 1,1 1,54 0,033 0,3
2000 1254 25 800
1380 ·
J/kgOC
1050 837 920
Apendice L
1 ECHO N.IO 337 kg/m2 = 0,751 W/m2 °C = 6,287 m2 °CIW Desfase = 12,512 horas
Peso K
=
IBBI
2cm
Fig. L.25
Composlci6n
Espesor m
Exterior Pavimento Bovedilla poliestireno Guarnecido yeso Interior
Conductividad W/moC
0,05 0,22+0,04 0,02.
1,1 0,26 0,3
Densidad kglm3
Calor. esp.
J/kgOC
2000
1380
850 800
1000 920
TECHO N.II
Peso K
= 317,5 kg/m2
= 0,367 W/ril2 °C
IBBI
=
Desfase
=
21 ,330 m2 °CIW 12,351 horas
2cm
Fig. L.26
Composici6n
Espesor
Exterior Pavimento Poliestireno expand. Tipo IV Honnigon ligero macizo Guarnecido yeso Interior
m
Conductividad W/m °C
Densidad kglm3
0,05 0,06 0,20 0,02
1,1 0,033 0,35 0,3
2000 25 1000 800
Calor. esp.
J/kgOC 1380 837 1050 920
543
Manual de climatizaci6n • Tomo II
TECHON.12 Peso K
IBBI Desfase
317,5 kg/m2 =0,367W/m2 °C
=
2 = =
15242 °CIW 11:467 m heras
~1~1~~il~~~~
.1 em
Fig. L.27
Espesor m
Composici6n
Exterior Pavimento Honnigon ligero macizo Poliestireno expand.Tipo IV Guarnecido yeso Interior
Conductividad W/moC
Densidad kglm3
1,' 0,35 0,033 0,3
2000 1000 25 800
0,05 0,20 0,06 0,02
Calor. esp.
J/kgOC 1380 1050 837 920
TECHON.13 263,5 kg/m2 = 0378 W/m2 °C
=
1~,899
2
m °CIW Desfase = = 10,285 horas
1 em
:1:~1111111111 ~em
Fig.L.28
Composici6n Exterior Gravilla Poliestireno Impermeabilizante Hormigon ligero macizo Guarnecido yeso Interior 544
Espesor m 0,05 0,06 0,01 0,15 0,02
Conductividad W/m °C
Densidad kglm3
Calor. esp.
0,81 0,033 0,19 0,35 0,3
1700 25 1100 1000 800
920 837 1680 1050 920
J/kgOC
Apendice L
TECHON.14
489,5kglm2 = 0395 W/m2 °C , K m2 °CIW , IBBI = 32616 Desfase = 11,684 horas Peso
=
Fig. L.29
Composici6n
Espesor Conductividad m W/m °C
Exterior Gravilla Poliestireno Impenneabilizante Bovedilla honnigon nonnal Honnigon ligero macizo Guamecido yeso Interior
0,05 0,06 0,01 0,22+0,04 0,05 0,02
0,81 0,033 0,19 1,54 0,35 0,3
Densidad kg/m3
Calor. esp. J/kgOC
1700 25 1100 1254 1000 800
920 837 1680 1050 1050 920
En la figura L;30 se representan los valores Y(j)/K de los diferentes cerramientos, especificnndose su peso en kg/m2 . Se puede observar una cierta clasificacion de su comportamiento tennico en funcion del peso por metro cuadrado, aunque existen ciertas discrepancias como las que se pueden detectar con los techos nOlI y 12 (misma composicion con diferente orden). Analizandose la fig. L.30 se clasifican en funcion de: POSICI6N DEL AISLANTE
. Techo invertido (aislante hacia el exterior) . Techo nonnal (aislante hacia el interior)
PESO
. Muy ligero ~ 200 kglm2 . Ligero ~ 350 kglm2 . Pesado ~ 550 kglm2
Asi se producen los siguientes gropos : - Techo ligero nonnalo invert~do muy ligero - Techo pes ado nonnal 6 invertido ligero - Techo invertido pesado
=* =* =*
Techo N° 10 Y 13 Techo N° 5 Y 11 Techo N° 14 545
Manual de clirnatizacilm. Torno II
COEFICIENTE Y(I)/U 0,06~--------------------------------------------~ 2 (9) 44~.5 . . (N;Techo) Pe,so Kg/m (10). 337 Bovedilla poliestireno (.13) .2.63,5. : . . 0,05 (12) 317,5: . . . (1) 538 . (5) 553 Bovedilla hormig6n normal . '(11·) 817',5 . . • . . . - .• - . - . - • - . - - . 0,04 (3) 5 7 4 ' . , (7) 589 (~L4:4~,?
0,03
___ . : .. _ . _ : _ . _ . _
(14) 489,5 Techo invertido
0,02 0,01
o
o
6
12
18
24
30
36
42.
'Fig. L.30. Caracteristica YO)/K de diferentes tech os
546
48
HORAS
ApendiceM Temperatllra equivalente en muros eortilla, tee/IDS ventilados y tee/lOS con eieloraso M.I.
MUTOS
cortina
M.I.I. Planteamiento del problema En Ia figura M.1 y M.2 se esquematizan dos secciones de este tipo de muros, los cuales estan constituidos por una chapa metalica exterior, una camara de aire ventilada, y una sucesi6n de capas de material que conforman el muro. En este tipo de cerramientos -Ia chapa metalica exterior posee una inercia termica despreciable, por 10 que actua en la pnlctica como una cortina respecto al muro real, entendido este como la sucesi6n de capas de material desde la ~amara de aire ventilada hasta el interior.
TECHON.l 156,8 kg/m2 = 0722 W/m2 °C , = 2,985 m2 °CIW Desfase = 5,579 horas
Peso K
=
IBBI
~ O~.,. -
c..
12 em
' ..... 2cm
~
110
lA 1 5 em
..
['0
Fig. M.l
Composici6n
Espesor m
Exterior Chapa metalica Camara de aire ventilada Poliuretano proyectado Ladrillo hueco Guarnecido yeso Interior
Conductividad W/m °C
Densidad kg/m3
0,023 0,49 0,3
40 1200 800
Calor. esp. J/kg °C
0,0015 0,02 0,12 0,015
837 920 920
547
Manual de climatizacion. Tomo II
TECHON.2 kg/m2 = 0,372 W/m2 °C K , m2 °CIW IBBI = 19433 Desfase = 9,260 horas
Peso
= 207,3
Fig. M.l
Composicion
Exterior Chapa metalica Camara de aire ventilada Poliuretano proyectado Ladrillo hueco Fibra de vidrio Ladrillo hueco Guarnecido yeso Interior
Espesor m
Conductividad W/moC
Densidad kg/mJ
0,023 0,49 0,041 0,49 0,3
40 1200 50 1200 800
Calor. esp. J/kgOC
0,0015 0,02 0,12 0,05 0,04 0,015
837 920 837 920 920
Pretendcmos en este apendice obtener una expresi6n de la temperatura sol-aire a considcrar, la cual sea equivalente 'a la ecuaci6n 7.16 para muros convencionalcs.
(M.l)
A partir de esta expresi6n seria de aplicaci6n la ce. 7.20 para el calculo de In temperatura equivalente del muro.
M.I.2. Flujos de calortransferidos En la figura M.3 se esquematizan los intercambios de calor que se suceden en un muro tipo cortina. 548
Apendice M
Interior
Fig. M.3. lntercambios de calor en un muro tipo cortina.
La temperatura sol-a ire se definia como aquella temperatura que deberiamos en el ambiente exterior, para que manteniendose constante, e igual al real el coeficiente de convecci6n y la temperatura de la pared exterior, nos contabilizara el calor absorbido por radiaci6n solar en dicha superficie, y por tanto pudieramos evaluar el flujo de calor entrante al muro mediante la expresi6n : co~siderar
q elft = he c (Tsa - Tpe)
(M.2)
siendo: qent - Calor entrante en el muro por la sup. exterior (W1m2) hCe - Coeficiente de convecci6n exterior, en nuestro caso el que ofrece una camara de air~ ventilada, (segitn norma NBE-CT-79, ref 6, sera 9,09 W/m2 °C).
Tsa - Temperatura sol-aire (OC) Tpe
-
Temperatura pared exterior del muro (OC)
EI flujo de calor entrante en el instante "nu (de acuerdo con el apendice K), se puede obtener mediante la ec. K.llS. Aplicando dicha expresi6n en funci6n de la temperatura superficial de la pared exterior y la temperatura seca interior, tenemos : CD
00
qelft
=
E
rO
Tpe(n-j) X(n) -
E
.
T lL(n-j) Yen)
(M.3)
foO
549
Manual de climatizaci6n. Torno II
donde: TsLCn-j)
- Temperatura seca interior en el instante (n-j).
XCn),Y(n) - Cocficicntes "X" y "Y" de Ia funcion de transferencia del muro Cdesde Ia capa de poliuretano al interior, ya que se considera la variaci6n de temperatura sobre Ia pared exterior, y no sobre el ambiente exterior). Asumicndose la temperatura interior eonstante con el tiempo, sera de aplicaci6n Ia simplificaci6n propuesta en K.130 :
qent
=
(i:
Tpe(n-j) X(n)) - K' TsL
(M.4)
}-o
siendo: K' - Coeficientc global de transmision de calor desde Ia capa de poliuretano hasta el intcrior (W/m2 °C).
K' =
Lj
-
1
Espesores de las capas (m)
Ai - Conductividades de las capas (W/m °C) hCj - Cocficiente de conveccion interior (9,09 W/m 2 °C) Norma NBE-CT-74 (ref 6). Los coeficientes X(n) para los muros cortina "A" y "B" representados en las figuras M.l y M.2, tomando como intervalo de tiempo una hora, son:
Muro cortina A X(O) = 1,1407
X(O) = 1,1407
X(l) =.0,12354
X(l) =·0,12357
X(2) =·0,035815
X(2) =·0,036348
X(3) =-0,027794
X(3) =-0,026666
............. 550
Muro cortina B
............
Apendice M
Pudicndose concIuir : 1) Los valores de los coef. X(n) en ambos muros son muy parecidos. 2) EI peso por el que se pondera la temperatura de la pare~ en el instante de calculo es 9,23 veces el peso por el que se pondera Ia temperatura de Ia pared en Ia hora prccedente, y 31,85 veces el peso por el que se pondera Ia temperatura de la pared
dos horas antes, etc .. 3) Con los anteriores pesos de ponderacion, y debido a la fluctuacion limitada de Ia
temperatura de la pared exterior, podemos afirmar que el suponer constante la temperatura de Ia pared exterior, e igual al momento actual, no supone un error apreciable en el calculo del calor entrante por la pared exterior mediante la ec. MA. 4) Admitiendo 10 anterior, y aplicando la ec. K.130, podemos transformar la ec. M.4
en:
q ent = K' (Tpe - T sL )
(M.S)
Igualando las expresiones M.2 y M.5, podemos despejar la temperatura sol-aire como:
Tsa
T se
he c
(M.6)
La temperatura de la pared exterior se podra estimar realizando un balance de energia sobre dicha superficie, asi : Calor entrante en 1a pared = Calor transfcrido por radiaci6n + Calor tnmsferido por convecci6n
(M.7) siendo: .. Calor entrante en la pared
q ent -- K I (Tpe - T sL )
551
Manual de c1imatizacion. Torno II
- Calor transferido por conveccion
(M.8) - Calor transfe.-ido por radiacion (pianos paralelos e infinitos)
a 1 u em
+
1
(M.9)
- 1
Up
6 Iinealizando :
1
1
1
-- +
U
Cl (Tem - Tpe)
(M.IO)
- 1
em
sicndo: CI
=
~p
a [(Tcm +273, 15)+(T pe +273,15)] [(Tcm+273,15)2 +(Tpc+273~15)2] - Coeficicnte de absorci6n chapa metalica - Coeficiente de absorci6n pared exterior .
Tcm - Temperatura chapa metalica (OC) Tai a
- Temperatura aire camara ventilada (OC) - Cte. Stefan-Boltzman (5,6697 10-8 W/m2 K4)
Aplicando el balance energctico (ec. K.7) :
K I (T - T_,.) he (T pe
iSLI
e
i -T
a
\
peJ
+
1
1
-- +
u em
1
--1
Cl (T
em
-T \
peJ
(M.II)
up
Como observamos Ia temperatura de Ia pared exterior (T pe) depende a su vez de dos nuevas temperaturas, Ia temperatura del aire en Ia eanlara ventilada (Tai) Y Ia temperatura que alcanza la chapa metalica (TcnJ Estas temperaturas se podran obtener mediante balances de energia sobre el aire que se encuentra en Ia camara, y sobre Ia chapa metalica, considerando en ambos casos inercia termica despreciable. 552
Apendice M
BALANCE DE ENERGIA TOTAL EN EL AlRE DE LA CAMARA
Calor total cedido por convecci6n chapa metalica = Calor total absorbido por convecci6n pared exterior + + Calor removido por movimientode aire en la camara - Calor total transferido por conveccion desde la chapa metiilica
he e Aneh. All. (Tem - Tat)
(M.12)
- Calor total transferido por convecci6n hacia la pared exterior
(M.13) - Calor total removido por el aire Va
Aneh . Esp. P a Cp a (Tal - T,e)
(M.14)
siendo: Velocidad aire en la camara (mls) Esp - Espesor camara de aire ventilada (m) Anch. - Anchura camara aire ventilada (m) All. - Altura camara aire ventilada (m) 3 Pa - Densidad del aire (~ 1,2 kglm ) CPa - Calor especifico aire (= 1000 J/kgOC)
Va
-
Aplicando el balance de energia al aire de la camara:
BALANCE DE ENERGIA EN LA CHAPAMETALICA EXTERIOR (por m2 sup.)
Calor absorbido por radiaci6n solar = Calor radiaci6n al exterior + Calor convecci6n al exterior + + Calor radiaci6n a la pared exterior. + Calor convecci6n a la camara de aire 553
Manual de clirnatizacion. Torno II
(M.16) - Calor absorbido por radiacion 'solar
sicndo:
IT (11, y) - Radiaci6n total incidcnte sobre la chapa metalica - Calor transferido por conveccion al exterior
q c1
=
he e (Tem - T se)
(M.17)
- Calor transferido por radiacion al exterior. Temp. exterior (T so) y con una
superficic mucho mayor que la pared del edificio.
6 lincalizando
qr1
(X
em
C2
(Tem - Ts)
(M.I8)
sicndo: C2 = ~ C~,lor ''',1
0
[(rcm +273,15)+(T se +273,15)] [(Tcm +273,15)2 +(T se+273,15)2]
transferido por conveccion al aire de la camara
q e2
=
he e (Tem - Tai )
(M.19)
....
- Calor transferido por radiacion a la pared exterior
1 1 (X
554
em
1
Cl (Tem - Tpe) - 1
(M.20)
Apendice M Aplicando eI balance propuesto en M.16 : ex em I 1"(-n. I , y) = he e (Tem - T se)
+
tx em
C2 (Tem - Tse) + he C. (T - Tai) + "em
l _ _ Cl (T -T) 1 em pe --+--1
+ __
1
ex em
(M.21)
ex p
M.l. 3. Expresioll de la temperatura sol-aire, y la temperatura
equivalente ell muros cortina Se han establccido tres ecuaciones, M.ll, M.15 YM.21" con tres incognitas, Tai' Telll' y TpC' de donde se puedc obtener Ia temperatura de Ia pared exterior, la cual sustituida en la ec. M.6, nos proporciona la temperatura sol-aire :
Tsa = Ai T
se
+
A2
txcm
I/Tl,Y) he
+
A3 TsL
(M.22)
c
en la que los coeficientes AI, A2 YA3 dependen de : - EI coeficientc de convecci6n exterior (hce = 16,67 W/m2 °C) - El coeficiente de convecci6n camara ventilada (hce = 9,09 W/m 2 °C) • -
EI coef. global transmisi6n calor pared exterior - aire interior (K') EI coeficiente de absorci6n chapa metalica (ex em) EI coeficiente de absorcion pared exterior (ex p) Temperatura interior (T sL ) Temperatura exterior (Tse = 29,2°C) La"radiaci6n solar incidente (IT(ll,Y» La velocidad del aire en la camara ventilada (va) EI espesor de la camara ventilada (Esp. = 0,03 rn.) La altura la pared exterior, 6 distancia entre aberturas (AlL)
555
Manual de clirnatizacion. Torno II
Se ha estudiado lID amplio ran go de variacion de dichas variables para los dos muros objeto de estudio, resultando :
Clem Cl
p
TsL(OC) . 2 IT(W/m) va(m/s) Alt.(m)
0,9 0,9 25 500 1 3
0,75 0,75 25 500 1 3
0,6 0,6 25 500 1 3
0,9 0,9 20 500 1 3
0,9 0,9 25 100 1 3
0,9 0,9 25 500 3 3
0,9 0~9
25 500 1 20
Muro A
Al A2 A3
1,008 1,000 0,223 0,216 -0,008 0,0
0,992 0,208 0,008
1,008 1,006 1,019 0,989 0,223 0,221 0,158 0,340 -0,008 -0,005 -0,019 0,011
Muro B
Al A2 A3
1,004 1,000 0,222 0,216 -0,004 0,0
0,996 0,209 0,004
1,004 1,003 1,009 0,995 0,222 0,220 0,156 0,341 -0,004 -0,003 -0,009 0,005
TABLA M.I. Valores de AI, A2, Y A3 para muros corfil,a
De 10 cual se conc1uye : 1) EI cocficicnte A3 es pnlcticamente nulo. 2) EI coeficicnte Ales pnlcticamente la unidad. 3) EI coeficiente A2 varia entre 0,34 y 0,16, fluctuando en gran medida con Ia vclocidad que alcanza el aire dentro de la camara ventilada, y la distancia entre aperturas en la pared exterior. Dicho coeficiente se puede aproximar a 0,25. 4) Por todo 10 anterior, la temperatura sol-aire en este tipo de muros se puede obtener mediante la expresion :
he c
556
(M.23)
Apendice M
considerando : a) Muro constituido por la sucesion de capas de material desde la camara de aire ventilada hasta el interior. b) EI coeficiente de conveccion exterior sera el de la camara de aire ventilada (9,09 W/m2 °C)
c) Se debera con~iderar como cocficiente de absorcion exterior de la pared, el de la chapa metalica exterior. 5) La temperatura seea equivalente se obtendnl mediante la sustitucion de)a cc. M.25 en 7.20, es decir :.
Tseq(n)
=
~ [T (n-") 0 25 J,
L." ~o
se
+
IX em
IrC:'Il'V)] YU)
h
Cc
K
(M.24)
La temperatura seea equivalente de los muros cortina (tabla 7.21) se han obtenido con las condiciones expuestas en este apendice.
M.2. Tee/lOS ventilados M.l.I. Planteamiento del problema En las figuras M.4, M.5 YM.6 se esquematizan diferentes configuraciones de techos ventilados, los cuales estall constituidos basicamente por lma serie de capas de material entre el exterior y una camara de aire ventilada, para a continuacion pasar a otra serie de capas hasta encontrar la superficie interior. Ell este tipo de cerramientos sf se debe considerar que ambas capas de materiales, tanto desde la camara de aire ventilada hasta eI exterior, como desde esa misma camara de aire al interior, poseen una inercia termica no despreciable, no siendo por tanto de aplicacion Ia base teorica desarrollada para muros tipo cortina. La temperatura que aIcanzan las superficies, es inferior a la que se alcanzaba en el caso de poseer una chapa metalica al exterior (muro cortina), por 10 cual asumiremos que los intercambios de energia por radiacion entre paredes son despreciables (0 bien ya incluidos dentro del coeficiente de conveccion). 557
Manual de climatizaci6n. Torno II
En este apendice se va a intentar desarrollar una expresion directa para estimar la temperatura equivalente en este tipo de eerramientos, dicha temperatura se define como aquella que deberia considerarse en el ambiente exterior, para que suponiendo el techo unicamente constituido por las capas de material que van desde la camara de aire ventilada hasta el interior, y manteniendo los coeficientes de conveccion existentes, obtuvicrnmos matematicamente el flujo de calor entrante como carga en el interior bajo Ia suposicion de regimen permanente de transmision de calor. . (M.2S) siendo:
K2 - Coeficiente global de transmision de calor des de Ia camara de aire ventilada hasta el interior (W/m2 °C).
1
K2 = 1
hc c Lj Aj hCj
hCe
capas +
E 1-1
L _I +
At
1
hC t
Espesores de las capas (m) - Conductividades de las capas (W/m °C) - Coeficiente de conveccion interior (5,88 W/m2 °C) Norma NBE-CT-74 (ref 6).
'-
-
Coef. de conveccion camara de aire (5,88 W/m 2 °C) Norma NBE-CT-74 (ref 6).
Tseq(n) - Temperatura seca equivalente exterior en el instante "n".COC) TsL
- Temperatura seca interior local (OC)
20m
Fig.M.4
558
Apendice M Espesor
Composicion
m Exterior Gravilla Doble tablcro ladrillo hueco Camara de aire ventilada Aislante Bovedilla honnigon nonnal Guamecido yeso Interior
Conductividad W/m °C
Densidad kg/m3
Calor. esp. J/kg °C
0,05 0,08
0,81 0,49
1700 1200
920 920
0,05 0,22+0,04 0,02
0,044 0,87 0,3
50 1254 , 800
837 1050 920
2 em
Fig.M.5
Espesor
Composicion
m Exterior Pavimento Doble tablero ladrillo hueco Camara de aire ventilada Aislante Bovedi1Ia honnigon nonnal Guamecido yeso Interior
Conductividad W/m °C
0,05 0,08 0,05 0,22+0,04 0,02
Densidad kglm 3
Calor. esp. J/kgOC
0,49
2000 1200
1380 920
0,044 0,87 0,3
50 1254 800
837 1050 920
1,1
2 em
Fig.M.6
559
Manual de climatizacion. Tomo II Composition
Espesor m
Exterior Impenneabilizante Doble tablero ladrillo hueco Camara de aire ventilada Aislante Bovedilla honnigon normal Guamecido yeso Interior
Densidad kg/m3
Conductividad W/m °C
Calor. esp. J/kg °C
0,01 0,08
0,19 0,49
1100 1200
1680 920
0,05 0,22+0,04 0,02
0,044 0,87 0,3
50 1254 800
837 1050 920
Se pueden plantear muchas variaciones sobre el techo de la figura M.6 en base a la existencia de una capa exterior con poca inercia.
M.2.2. Filijos de calor transferidos En la figura M.1 se tiene una representacion esquematica de este conjunto de techos.
IT (¥,tz)
•
T
se
Exterior
Camara aire ventilada
q.
I
Interior Fig. M. 7. Esquema techos ventilados
560
Apendice M
Para la estimaci6n del flujo de calor entrante (qent), es necesario conocer la relaci6n existente entre los flujos de calor que alcanzan la camara de aire ventilada, y que en definitiva producen una temperatura de, equilibrio del aire de la misma, la cual determinant el flujo calorifico finalmente transferido por el techo, (carga en nuestra instalacion).
I
FLU,IO DE CALOR DESDE LAS CAPAS SUPERlORES
EI flujo de calor que alcanza la camara de aire ventilada desde eI extcrior se puede obtener como si se tratara de un cerramiento normal, por tanto sera de aplicaci6n la ec. 7.10 y 7.21, es dccir :
(M.26) donde:
(M.27)
K1
-
Coeficiente global de transmisi6n de calor desde el exterior a Ia camara de aire ventilada (W/m2 °C).
1
Kl 1 -- +
he e Li
capas
L
L
_t
+ --
1-1
At
he c
1
- Espesores de las capas (m)
Ai - Conductividades de las capas (W/m °C) hCe - Coeficiente de convecci6n exterior (20,0 W/m2°C) Norma NBE-CT-74 (ref 6). hCe
-
Coer. de convecci6n camara de aire (5,88 W/m 2o C) Norma NBE-CT-74 (ref 6).
561
Manual de climatizacion. Tomo II
Tscq1(n) -
Temperatura seca equivalente exterior en el instante "n" del trozo de cerramiento comprendido entre la camara de aire y el exterior.(OC)
Tai
-
Temperatura seca del aire en la camara (OC)
Y I (j)
~
Coef. Y de la funci6n de transferencia del lramo de cerramiento comprendido entre la camara de aire y el exterior (W/m2 °C)
1'}"( 11, y) - Radiaci6n solar incidente (W1m2) a. te
-
Coeficiente de absorci6n techo exterior
FLUJO DE CALOR HACIA LAS CAPAS INFERlORES EI fllUO de calor que atraviesa la superficie de la camara de aire ventilada y se dirige hacia el interior del local q' se puede obtener mediante Ia ec. K.118, por tanto:
q 'en) =
L 1'-0
Taln-j) X 2(n) -
L
T sL(n-j) Y2(n)
(M.28)
1'-0
Xin),Yin) - Cocficientes "X" y "Y" de la funci6n de transferencia del muro (des de la camara de aire ventilada hasta el interior, W/m2 °C) Asumicndose la temperatura interior de nuestra instalaci6n constante con el tiempo, sera de aplicaci6n la simplificacion propuesta en K.130 :
(M.29)
sicndo:
K2 - Coeficiente global de traqsmisi6n de calor desde la camara de aire ventilada hasta el interior (W/m2 °C).
1
K2 =
1
-- +
hc c
562
capas
L
i-I
L
1
_I + -
Ai
hc,
Apendice M
hCi - Coeficiente de conveccion interior (5,88 W/m 2 °C) Nonna NBE-CT -74 (ref. 6). Los coeficicntes Xz{n) para los cerramientos estudiados toman los valores : X(O) = 0,89729 X( 1) =-0,17043 X(2) =-0,013187
X(3) =-0,010047
Pudi6ndose conduir : 1) EI factor por el que se pondera la temperatura del aire de la camara ventilada en cl instantc de calculo es 5,265 veces el factor de ponderacion de Ia temperatura de dicho aire en la hora precedente, y 68,043 veces el factor por el que se pondera esa misma temperatura dos horas antes, etc .. 2) Con los antcriores factores de ponderacion, y dcbido ala fluctuacion limitada de la temperatura del aire en la camara, podemos afinnar que el suponer constante la temperatura de dicho aire, e igual al momento actual, no supone un error apreciable en el calculo del flujo de calor entrante desde la camara de aire hacia el interior. 3) AdmiLiendo el mismo razonamiento, y aplicando la ee. K.130, podemos transfonnar la ec. M.29 en :
(M.30)
FLUJO DE CALOR REMOVIDO POR EL AIRE EI flujo de calor que es removido por el movimiento del aire de la camara ventilada por metro cuadrado de superficie exterior, es funcion de In superficie de aberturas cxistentes (S), y de Ia velocidad inducida del aire (por circulacion natural 0 forzada). Asi:
(M.31)
563
Manual de climatizaci6n. Tomo II donde: - Supcrficie aberturas camara ventilada (Se divide por 2, por suponer la mitad de superficie extrayendo aire, m2 ). A - Superficie total de techo (m2) va - Vclocidad inducida al aire (mls) Pa - Densidad del aire (~ 1,2 kglm3) CPa - Calor especifico aire (= 1000 J/kgOC) RALANCE DE ENERGIA Rcalizando el balance de energias propuesto sobre el aire de la camara, y supuesto este con inercia tcrmica despreciable : I
q;=qaj+
qI
(M.32)
Sustituyendo :
y despejando 141 temperatura del aire en Ia camara:
T at.(n)
(M.33)
Una vez estimada la temperatura del aire, y despreciando posibles Intercambios de calor por radiaci6n dentro de Ia camara ventilada (0 asumiendo que ya han sido contabilizados con los coeficientes de convecci6n), la temperatura seca equivalente buscada se obtendra aplicando Ia ec. 7.21 para el tramo de cerramiento existente entre la camara de aire ventilada y el exterior, es decir 564
Apendice M
1
v P Cp S/2 y 2(i) j(n-i)+K T_ r + a a a T (n-i) _'_ 2 seq A se K ~ 2 Tseq(n)=-------------------VaP aCp a S/2 co
[
~ KIT
L..,;
0lI.J,
Kl+K2
+
(M.34)
A
Como obscrvamos en Ia cc. M.34, en este tipo de techos dcberemos estimar dos veces una temperatura equivaicnte, 10 cllal complica ciertamente el proceso matematico. Finalmente scfialemos que deberemos afiadir dos nuevas variables: 1)
La temperatura seca deliocal (g~eralmente TsL=25°C)
2) Los metros cubicos por segundo y metro cuadrado de techo que se indllcen. va S/(2A). Dcberemos utilizar un valor consen'ativo (menos aire inducido con 10 que
se alcanzara mayor temperatura en el aire de la camara ventilada). Los valores standard que se proponen son : S/A = 0,01 m2 aberturas/m2 superficie techo va = 0,5 mls luego: 3 Wlln
2
°c
La temperatura equivalente de los techos ventilados(tabla 7.24) sc ha obtenido con las consideraciones expuestas en este apendice.
M.3. Teclzos con cieloraso Una fonna constructiva muy tipica de techos es sustituir la capa interior de guarnecido de yeso por una capa de escayola no aplicada directamente sobre la bovedilla, sino a una cierta distancia de la misma, y que en general oscila entre lOy 30 cm. Esto crea una camara de aire por Ia que se suelen conducir las Iineas electricas y en algunos casos los conductos de aire acondicionado si hay suficiente espacio. Desde un punto de vista temlico dicha camara de aire aport a en primer lugar una cicrta resistencia a Ia entrada de calor, y que como sabemos por la tabla 7.2 se puede 565
Manual de climatizacion. Torno II
establceer en 0,21 m2 °eIW. Esto haee disminuir en eonseeueneia el eoefieiente global de intereambio tcrmico "K" (vcr tabla M.2). Dc otra parte dicha resistencia tonnica modi fica la distribucion de temperaturas existente, por 10 que tambicn afecta a In energia almacenada en las diferenles capas del tceho, y en conseeuencia a la inercia que presenta el cerramiento. Es decir, se modifican asimismo las tcmpcraturas equivalcntes a considcrar. En Ia figura M.8 se han representado las temperaturas equivalentes para cI mes de Agosto, de los techos que hemos establecido como standard, el nonnal ligero, normal pes ado y el invertido pesado. Todos se han calculado con y sin la camara de aire comentada y en condiciones standard de calculo (tabla 7.22), excepto Ia consideracion de color oscuro.
Techo normalligero
Pavimento Bovedilla poliestireno Guamecido de yeso
0,05 m 0,22+0,04 m 0,02m
Techo normal pes ado
Pavimento Impenneabilizante Boved. hormigon nonnal Guarnecido de yeso
566
Pavimento Bovedilla poliestireno Camara de aire Escayola
0,05m 0,22+0,04 m 0,2 m 0,02m
Techo normal pesado can cieloraso 0,05m 0,01 m 0,22+0,04m 0,02m
Techo invertido pesado
Gravilla Poliestireno Impermeabilizante Boved. honnigon nonnal Honnigon ligero macizo Guamecido de yeso
Techo normalligero con cieloraso
0,05 In 0,06m 0,01 m 0,22+0,04 m 0,05 m 0,02m
Pavimento Impenneabilizante Boved.honnigon nonna Camara de aire Escayola
0,05111 0,01 m 10,22+0,04m 0,2m 0,02m
Techo invertido pesado con cieloraso
Gravilla Poliestireno Impenneabilizante Boved. honnigon nonna Honnigon ligero nOnna Camara de aire Escayola
0,05 n1 0,06m 0,01 m 10,22+0,04 m 10,05 m 0,2m 0,02111
Apendice M
Dc la observnci6n de la figura podemos extraer las siguientes conclusiones : 1) La cxistencia dcl cicIo raso hace disminuir tanto el coef. global dc transmisi6n de calor "K", como In temperatura equivalente maxima a considcrar~ oscilando entre 1°C yO, I °C en funci6n de la inercia del cerramicnto. 2) Con cieloraso Ia osciIaci6n a considerar de temperatura equivalcnte es mcnor. 3) EI flujo de calor ganado por metro cuadrado (T sL situaciones se observa en Ia tabla M.2
K W/m 2 °C
= 25°C) en diferentcs
qrnax W
Teq 18 h.
q18 h.
°C
°C
W
Tcqmax
Ligero sin cicIo raso
0,751
41,1
12,091
34,8
7,360
Ligero con cicIo raso Pcsado sin cicIo raso Pcsado con cielo raso lnvertido sin cicio raso Invertido con cicIo raso
0,649 1,191 0,953 0,395 0,365
40,2 39,2 38,5 38;1 38,0
9,865 16,912 12,866 5,174 4,745
34,5 35,3 35,3 36,2 36,6
6,166 12,267 9,816 4,424 4,234
TABlA M.2. Comparacion tee/lOs con)' sill cieloraso
Del anal isis dc los resultados podemos concluir : - Tanto cl flujo de calor como el coeficiente global de intercambio sc modi fica por la existencia del cieloraso. - La temperatura equivalente maxima (::::23 h. solares) es distillta. - La temperatura equivalente a las 18h. solares (:::: 16 h. civiles) practicamente no se modifica por la existeneia del cielo raso. Por 10 tanto, y de forma aproximada se aconseja utilizar las mismas temperaturas equivalentcs del techo sin la existeneia del cieloraso, pero modificado cl coeficiente global de calor. Para un resultado mas exacto se deberian desarrollar las temperaturas equivalentes para esta nueva situaci6n.
567
Manual de climatizaci6n. Torno II
Horizontal · 24 Agosto
Temp. equivalente (Oe)
42 41 40 39 38 37 36 35 34 33 32 31 30
-
UG-S.C.
-t- UG-C.C. "'PES-S.C. ---- PES-C.C.
*"lNV-S.C.
-
..
-
-."
' I
.,
..
..
"."
..
-+-INV-S.C. I
.. ... .. ..
f ~
.. ..
- - -. - - -. - - -. - - - .- - - ~ - - -. - - -. - - - .- - - ~ - - -, - - '. - . 0
2
4
6
8
10
12
14
16
18
20
22
s.c. - (SIN
CAMARA AlAE) C.C. - "(CON CAMARA AlAE)
24
Hora solar (h)
Fig. M.B. Comparacion temp. elJuivaientes techos con y sin cieloraso.
568
ApendiceN Coejicientes "nvi" para elementos adicionales en acristalamientos
N.J. Elementos jijos
N.J.l. Carpinteria metaiica La construccion tipica de ventanales metalicos, comporta una superficie acristalada de alrededor del 960/0, restando el 4% para los propios marcos del ventanal, luego la energia que realmente ha atravesado la superficie acristalada sera : (N.I) donde: ITRcm(ll,Y) - Energia que ha atravesado una ventana con carpinteria metalica (W/m2), tomando como area total de la ventana eI hueco de la misma. - Energia que ha atravesado una superficie completamente acristalada con vidrio comun sin marco. Tabla 7.25 (W/m2 ) No se ha modificado la proporcion de energia radiante respecto al vidrio comun.
nv
=
0,96
mr
=
97,3 %
(N.2)
N.J.2. Carpinteria de madera Los ventanales de madera conllevan lID marco muy superior al existente en ventanas con marco metalico, asi en general podemos suponer que el marco de madera ocupa sobre el 15% del hueco total de ventana, 0 10 que es 10 mismo la superficie acristalada alcanza el 85% del hueco, asi en este caso : (N.3) 569
Manual de climatizacion. Torno II
donde: ITRcM (l1,Y) - Energia que atraviesa una ventana con carpinteria de madera. (W/m2), tomando como area total de In ventana el hueco de Ia misma. Igualmente al caso anterior no se ha modificado Ia proporci6n de cnergia radiante respecto al vidrio comun.
mr
nv = 0,85
N.1. 3.
97,3 %
(N.4)
Vidr;os de distinto espesor
EI procedimiento de calculo es el desarrollado en el capitulo 7.2.2.2.1. aunque evidentemente con distinto espesor a1 considerado standard (4 mm). Como observamos en Ia ec. 7.59 ia cantidad de energia fmalmente transferida dependc de como se distribuye el porcentaje que sea radiaci6n directa y difusa. Con vistas a poder estimnr un coeficiente "nv", supondremos dos casos extremos, de una parte que la energia dirccta sea eI 80% del total (caso de sup. horizontales) y de otra que Ia energia difusa sea el 100% (orientaci6n norte en el hemisferio norte). Con estos dos casos extremos se puede evaluar el cociente :
nv
I TR
esp( 1'\ , Y)
(N.S)
ITR(ll, y)
ITR esp( 11, y) -
Energia que ha atravesado una ventana con un cristal cornun de espesor "esp" (m).
En la figura N.l se observa dicho coeficiente en funcion del angulo, para 4 espesores diferentes y para los dos tipos de radiaci6n solar comentados (todo radiaci6n difusa, 0 s610 cI 20%). Asimismo en la figura N.2 se representa "mr" (tanto por cien de energia radiante, calculado mediante 7.60).
570
Apendice N Valor de nv
0,98 - ' - ~ .. -
.~
- - -
.
- .
.. ..
'"
...
.. .. .. .. .. .. .. .. .. ..
..
Espesor 6 mm SOO.4 Rad. directa Espesor 6 mm 100 %Aad. difusa
0,96 ----:--............-'.:....':--...............:.-:..~ . .:.'",:,-";"_':":"-.:..':,,'-..,:'__':...:..':.."~~:""';"1::::::~~~1-=--=--=--=--=--=-- Espesor Espesor ...4.-~01----- Espesor A..---'''----- Espesor
0,92
.
.
I
I
..................
.
, , - .......................... ,
•
•
. '"
•
8 mm 80% Rad. directa 8 mm 100 %Rad, difusa 10 mm SOOk Rad, directa 10 mm 100 %Rad, difusa
......... . •
0.9~----------------~---------~-~ 60 70 90 o 10 20 30 40 50 80 Angulo de incidencia
Fig. N.I. Coefidente "nv" para cristal comun y diferentes espesores.
De la observaci6n de las figuras N.I y N.2, podemos establecer las siguientes relaciones entre el factor "nv", 1a proporci6n de energia radiante y el espesor del cristal, (siendo para todo tipo de radiaci6n practicamente el mismo).
nv = 1 - 10 (esp - 0,004) mr = 97,3 - 700. (esp - 0,004)
(N.6) (N.?)
donde: esp - Espesor del cristal (m). 100
r-----
Espeso,
~:~::~~
4 mm 100% Aed. di1utle
: ~~ ~gg::: ~:~: ~:~~::
~......---....,
Esp_o, 10 mm 100% Red . difuse
, eap."O' 4 min 80% Red . dirltCta • , esp.so. 6 mm 80% Red . direct. ::..,..:...---~ , , , " , , ,'. ~s~.!-O~ f8!"!,,,.8~~ ~e~.d~r~~ ,~ , ! , , ~sp . .,o, 0 mrn 8".,.. ,ned. dor~.. .
80~----------------------~ 70 80 40 90 30 50 60 o 10 20 Angulo de incidencia (0)
Fig. N.l. Tanto por den de energia radiante funcion espesor cristal
571
Manual de clirnatizaci6n. Torno II
N.l.4. Vidrio doble EI proceso de reflexiones y absorciones de energia radiante se esquematiza en la figura N.3. En la misma debemos de tener en cuenta que los coeficientes de reflexi6n y transmisi6n utilizados para ambos cristales vienen dados mediante las ecuaciones 7.41 Y7.42 respectivamcnte, y los subindices indican el vidrio referenciado, "I" exterior y "2" interior.
Radiaci6n reflejada
Radiaci6n incidente
Interior Radiaci6n transmitida
Fig. N.3. Proceso de rej1exion y transmisiOn en vidrio doble.
Como sicmprc los coeficientes globales de transmisi6n y reflexi6n para la radiaci6n dirccta se obtendnln por sumas de series, y el coeficiente de absorci6n por difcrencia de los anteriores respecto a la unidad. - Coeficiente global de transrnisi6n :
y sumando la serle : t
572
=
(N.8)
Apendice N
donde: PI - Coer. reflexion crist'll" 1", exterior (ec. 7.42) P2 - Coer. rcflexion crist'll "2", interior (ec. 7.42) 't 1 - Coer. transmisi6n crist'll" 1", exterior (ec. 7.41) 't 2 - Coef. transmisi6n crist'll "2", interior (ec. ? 41 )
- Coeficiente global de reflexi6n :
y sumando Ia serie :
p
(N.9)
- Coeficiente de absorci6n :
(N.lO)
Respecto a la radiacion difusa recordemos que su coeficiente de transmisi6n era el 89% del que se tiene ante incidencia nonnal (ec. 7.57), y que el coeficiente de absorci6n se mantenia pnicticamente constante (fig. 7.23), por tanto la energia transmitida ante Ia existencia de un doble crist'll se evaluara como: 1 TR doble( 11 , y) = 1TR doble radiante( 11 , y) + 1TR doole conv-cond. ( 11 , Y)
1TRdobleradiante(l1 ,V)
=1 D(Tl ,V)
't 1't 2
(N.ll)
, 't 01 't 02
+1/11, y) 0,89 - - - I - PIP 2 1 - P 01 P 02
573
.Manual de climatizaci6n. Torno II don de se ha supuesto el mismo porcentaje de energia absorbida que se transmite al exterior (65%) y al interior (35%) que en el caso de vidrio simple, ya que esto es funcion de los coeficientes de conveccion y de las ternperaturas existentes en ambos medios. Como en los anteriores casos definimos el coeficiente "nv" mediante el cociente :
I TR
nv
doble('"
Y)
(N.12)
I TR(", Y)
dicha relacion es funcion de la proporcion de energia recibida en forma directa y dirusa, y del angulo de incidencia, en la figura N.4 se contempla 'su variacion.
nv 0,9
~
Doble crista!. Espesor de cada uno 4 mm.
-
~
. -
•
---
-.-·
...
...
... •
........... :
-
...
-
...
.. - .. - .. I
- ........ - .......................... - .... -
...
..................... - ................. I
..
-
0,4
,
,
. , . , . . ; ..... - .. -.- .......... -; ..........
............
•
..
- ...... -,- ........ : ........ -." ........
0,5
·
..
. I
-
.
- ........ - ........ I
I
,
. I
. I
. I
0,3 0,2 0,1
.......... -:- ........ : . - .... -:- ..... - . : - - .....
<- ......... ; ......... -: .... - . . . : ............
o~--~----~----~----~--~----~----~----~--~
o
10
20
30
40
50
60
70
80
90
Angulo de incidencia (0)
Fig. N.4. Coeficiente "nv" para cristal doble.
Como observamos el cornportamiento es diferente para radiacion difusa 0 radiacion directa, aWlque para pequefios angulos de incidencia de la radiacion directa (menores de 50°) el factor "nv" es sensiblemente igual. 574
Apendice N
En el ealeulo de eargas termieas siempre estamos intentando ealeular las demandas tenniens superiores, y estas se daran cuando la energia transmitida sea elcvada, es decir euando eI angulo de ineideneia no sea exeesivamente elevado, es por 10 que se torna un valor constante de nv = 0,92. La proporci6n de energia radia.nte se puede estimar mediante Ia expresion N.13, la eual se representa en In fig. N.5:
mr
I TR
100
I TR mr (%)
doble radianti 11 , Y)
doble radiante( 11, y) +
I TR
(N.l3)
doble conv-condJ 11 , Y)
Vidrio dobre de Igual espesor ambos cristales (4 mm)
100~--------------------------------------~
95
. - - .. - .. ,
, - ............. -
..
,
90 . -"-,- ---. ---
, •
, .......... - .... . . I
•
85 80 75
,
,
,Rad.,eO% directa,
- - - -:" " - - : - - - -:- . - -
,
,
:Rad."10~% ciif~;~'~' "-:"-
70 65
•
*
•
- - " -, - . . . ; - - . - " " . -
• ~
•
,
•
•
. - - . ,- . . . -, . . " - ,- - - . -, - .
Angulo de incidencia
e)
Fig. N.S. Proporcion de energla radiante. Cristal dohk.
Dc la simple observacion de la figura N.5 se concluye que In proporeion de energia radinnte es del 91%, y unieamente con angulos de incidencia grandes (>60°) es apreciable Ia proporci6n de energia intercambiada por conducci6n - convecci6n.
N. 1. s. Vidrios especia/es En la actualidad existen muchas variantes de vidrios utilizados en ventanas, tanto debido a su diferente composici6n, como a tratamientos especiales de su superficie. La modifieaci6n de la composici6n del cristal (vidrios atermieos) aumenta en general su coeficiente de extinei6n, 10 que provoca mayor energia absorb ida y por tanto menor energia transfcrida al interior. 575
Manual de climatizacion. Tomo II
EI tratamiento superficial del vidrio aumenta 141 energia reflejada, 10 que se traduce en una menor transmisi6n de radiaci6n 411 interior. Ante la utilizaci6n de vidrios especiales deberemos conocer (por catalogo), su influencia sobre 141 radiaci6n total tmnsmitida en comparaci6n con cI vidrio comllll. En las tab las N.t y N.2 sc racilitan los ractores solares (referidos a radi41ci6n directa) de los principales tipos de vidrios utilizados en Espana, extraidos de 141 referencia 14. A partir de este valor y mediante 141 ee. 7.60 se estima el valor de "nvt Produrto
Crista nola Planilux Cristanola Planilux Cristanola Planilux Cristanola Planilux Cristanola Planilux Cristanola Planilux Cristanola Planilux Cristanola Planilux Cristanola Planilux Parsol Parsol Parsol Parsol Parsol Parsol Parsol Parsol Parsol Parsol Parsol Parsol Pink-Rosa Pink-Rosa Pink-Rosa Pink-Rosa Rellcctasol K Rctlcctasol G Antclio Antclio Antclio Antclio Antclio Antclio Antelio AnteJio AnteJio Cool-Lite SS 108
576
Espesor (nun)
Color
Fador solardlncta
nVI
2 3 4 5 6 8
Transparente Transparente Transparente Transparente Transparente Transparente Transparente Transparente ' Transparente Oris Oris Oris Oris Bronee Bronce Bronee Bronee Verde Verde Verde Verde Rosa Rosa Rosa Rosa
0,88 0,88 0.88 0,87 0,85 0,82 0,80 0,79 0,79 0,69 0,64 0,60 0,48 0,68 0,63 0,59 0,47 0,65 0,61 0,57 0.48 0,81 0,80 0.78 0.70 0,52 0,42 0,60 0,60 0,59 0.57 0,65 0,65 0,46 0,44 0,41 0,17
1 1 1 0,99
10 12 15 4 5 6 10 4 5 6
10 4 5 6
10 4 5 6 10 6 6 5 6 8
10 6 8 5 6 8 6
Clair Clair Clair Clair Argent Argent Havane Havane Havane Plata
0~7
0,93 0~1
0,90 0,90 0,78 0,73 0.68 0,55 0,77 0,72 0.67 0,53 0,74 0,69 0,65 055 0,92 0,91 0,89 0.80 059 0,48 0,68 0,68 0.67 0,65 0,74 0,74 0,52 0,50 0.47 0,19
Apendice N
Cool-Lite SS 114 Cool-Lile SS]:W Co()I-Lite SS 132 Cool-Lite TS120 Cool-Lite 5CI08 Cool-Lite SCI 14 Cool-Lite 5C 120 Coot-Lite TO 108 Cool-Lile TE 11 0 Cool-I.ile TE 115 Cool-Lite TB 1)0 Cool-Lite TB140 Cool-Lite SS208 Cool-Lite 55214 Cool-Lite SS220 Cool-Lite SS232 Cool-Lite SS308 Cool-Lite 5S314 Cool-Lite SS320 Cool-Lite 5S332 Cool-Lite 5R 132 Cool-Lite SS408 Cool-Lite SS414 Cool-I -ite SS420 Cool-Lile S5432 Cool-Lite 5S508 Cool-Lite SS514 Cool-I.ite 5S520 Cool-Lite 55532 Cool-Lite 55608 Cool-I.ite S5614 Cool-l.ite 5S620 Cool-Lite SS632 Cool-I -ite PB 108 Cool-Lite PHl]4 Cool-Lite PI1120 Cool-Lite PH 135 Cool-Lite PR408 Cool-Lite PB414 Cool-Lite PB420 Cool-Lite PH435 Cool-Lite L8120 Cool-Lite L8] 35 Cool-Lite DB 120 Cool-Lite DB 135
6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6
Plata Plata Plata Plata Cromo Cromo Cromo Oro Ocre Oere Azul Azul Bronec Bronec Bronec Bronec Oris Oris Oris Oris Oris Verde Verde Verde Verde Pink-Rosa Pink-Rosa Pink-Rosa Pink-Rosa Ambar Ambar Ambar Ambar Pastel Azul Pastel Azul Pastel Azul Pastel Azul Pastel Aguamarina Pastel Aguamarina Pastel Aguamarina Pastel Aguamarina Pastel Bronec Claro Pastel Bronec Claro Pastel Bronce Oseuro Pastel Bronce Oscuro
023 0.29 0.41 0.29 0.17 0,25 0,31 0.16 · 0,22 0.24 0,37 0,46 0,20 0.23 0,26 0,36 0,20 0,23 0,26 0,37 0,43 0.19 0,22 0,26 0,35 0.18 0,23 0.28 0,43 0.19 0,23 0.26 0,37 0,21 027 0,33 0,46 0,22 0,25 0.27 0,37 0,31 0,46 0,30 0,45
026 03) 0.47 033 0.19 029 035 0.18 0,25 027 0,42 0,52 0,23 0,26 0,30 0,41 0,23 0,26 0,30 0,42 0,49 0,22 0,25 0,30 0.40 020 0,26 0,32 0,49 0,22 0,26 0,30 0,42 024 0,..11 038 052 0,25 0,28 0,31 0.42 0,35 0,52 0,34 0,51
TAB/AN.I. Fuctores solares (nvJ de vidrios simples comerciales
577
Manual de climatizacion. Torno II Producto
Espcsor (nma)
Climalit I Planilux Ctimalit I Planilux Climalit I Planilux Climalit / Planilux Climalit I Planilux Climalit I Planilux ClillUllit / Eko Climalit / Eko Climalit / Phmithcrm Climalit I Parsol Climalit / Parsol Climalit / Parsol Climalit / Parsol Climalit I Parsol Climalit / Parsot Climalit / Parso! Climalit ! Parsol Climalit / Parsol Climalit / Parsol Climalit / Pink-Rosa Climalit / Pink-Rosa Climulit / Pink-Rosa Climalit I Pink-Rosa Climalit I Refleetasol K Climalit / Rcfleetasol K Climalit I Antclio Climalit / Antclio CHmalit / Antelio Climalit / Antelio Climalit / Antelio Climalit / Antelio Climalit / Antelio Climalit I Antelio Climalit I Antelio Climalit / Cool-I,ite SS 108 Climalit I Cool-Lite SS114 Climalit I Cool-Lite SS 120 Climalit I Cool-Lite SS132 Climalit / Cool-Lite TS120 Climalit I Cool-Lite SC 108 Climalit I Cool-Lite SC114 Climalit / Cool-Lite SC120 CUmalit / Cool-Lite TO 108 CUmalit / Cool-Lite TE 11 0 Climalit / Cool-Lite TE 115 Climalit / Cool-Lite TB 130
578
4/4 SIS
6/5 6/6 8/8 10/8 4/4
6/6 6/6 4/4
515 6/6 8/10 4/4 5/5 6/6 8/10 6/6 8/10 4/4
515 6/6 8/10 6/6
Color
Transparentc Transparente Transparentc Transparcnte Transparentc Transparcnte Transparclltc Transparentc Transparentc Oris Oris Oris Oris Bronee Bronee Bronee Bronee Verde Verde Rosa Rosa Rosa Rosa
6/6 5/5 6/6 8/8 8/10 5/5 6/6 8/8 6/6 8/8 6/6 6/6 6/6
6/6 6/6 6/6
6/6 6/6
6/6 6/6 6/6 6/6
Clair Clair Clair Clair Havane Havane Havane Argent Argent Plata Plata Plata Plata Plata Cromo Cromo Cromo Oro Oere Ocre Azul
Fm:tor SOilll''''rt'C'la
0.76 0.75 0.73 0.72 0.68
0,66 0.74 0.71 0,65 0.58 0.54 0.49 0,38 0.58 0,54 0,48 0,37 0,47 0,36 0,72 0,69 0.67 0,59 0,45 0.31 0,52 0.51 0.50 0,48 0.36 0.34 0,30 0,58 0,57 0,13 0,17 0.22 0.37 0,22 0.12 0.19 0,24 0,11 0,15 0,17 031
ltV,
0,86 O.R5 0,83 0,82 0.77 0.75 0,84 0.81 0,74 0,66 0,61 0,56 0,43 0,66 0,61 0,55 0,42 0,53 0041
0,82 0,78 0,76 0,67 0.51 0,35 0,59 0,58 0,57 0,55 0,41 0,39 0,34 0,66 0,65 0,15 0,19 0.25 0.42 0,25 0,14 0,22 0,27 0,13 0,17 0,19 0,35
Apendice N
Climulit / Cool-Lite TR 140 Climalit I Cool-I.ile SS208 Climalit I Cool-I.ile SS214 Climalit .I Cool-J ,ite SS220 Climalit / Cool-Lite SS232 Climnlit ! Cool-I.ile SS308 Climalit I Cool-I.ite SS314 Clilllnlit ! Cool-Litc SS320 Climalit I Cool-l.ite SS332 Climalit / Cool-Lite SR 132 Climalit / Cool-Lite SS408 Climalit.l Cool-Lite SS414 Climalit / Cool-Lite SS420 Climalit I Cool-Lite SS432 Climalit / Cool-Lite SS508 Climalit / Cool-Litc SS514 Climalit I Cool-Lite SS520 Climalit I Cool-Lite SS532 ClimaIit / Cool-Lite SS608 Climalit / Cool-Lite SS614 Climalit / Cool-Lite SS620 Climalit / Cool-Lite SS632 Climalit / Cool-l.ite PJ3108 Climalit I Cool-Lite PB 114 Climalit i Cool-Lite PA 120 Climalit I Cool-Lite PB 135 Climalit / Cool-Lite PB408 Climalit / Cool-J lite PB414 Climalit / Cool-I,ite PB420 Climalit I Cool-I.ite PB435 Climalit / Cool-I.ite LB 120 Climalit / Cool-I.ite LB 135 Climalit I Cool-Lite DB 120 Climalit / Cool-I.ile DB 135
6/6
6/6 6/6 6/6
6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6 6/6
Azul Bronec Bronee Bronee Bronce Oris Oris Oris Oris Oris Verde Verde Verde Verde Pink-Rosa Pink-Rosa Pink-Rosa Pink-Rosa Ambar Ambar Ambar Ambar Pastel Azul Pastel Azul Pastel Azul Pastel Azul Pastel Aguamarina Pastel Aguamarina Pastel Aguamarina Pastel Aguamarina Pastel Bronee Claro Pastel Bronee-Claro Pastel Bronee Oscuro Pastel Bronec Oscuro
0.39 0.12 0.15 0,18 0.26 0,12 0.15 0.18 0.27 0.34 0,12 0.15 0.18 026 0,13 0,16 0,21 0,34 0.12 0,15 0,18 027 0,15 021 0,25 0.37 0,13 0.17 0,18 027 0,24 0,37 0,24 0.37
0.44 0.14 0,17 0,20 0,30 0.14 0.17 0.20 031 0,39 0,14 0,17 0.20 0.30 0,15 0,18 0,24 039 0,14 0,17 0,20 0.31 0,17 0,24 0,28 0,42 0,15 0,19 0,20 0,31 0,27 0,42 0,27 0,42
TABLA N.l. Factores solares (nv) de vidrios dobles comerciules.
Para vidrios de color si desconocemos Ia transmisi6n de energia solar que producen podemos de forma aproximada utilizar Ia tabla N.3.
579
Manual de climatizaci6n. Tomo II
Material
nv
Vidrio color Ambar Bronce Rosa Azul Oris Verde Opalescente claro Opalescente oscuro
0,85 0,77 0,92 0,80 0,78 0,74 0,43 0,37
TABlA N.3. Factores so/Qres standard (nv) para vidrios de color
En principio y a falta de datos mas concretos de ensayos supondremos que Ia energi a transmitida en forma radiante se mantiene en un 90%.
N. 2. Elementos moviles
N.2.1. Cortinas interiores Los intercambios energeticos que se producen ante la existencia de cortinas interiores se representan en la figura N.6.
TRANSMISION
ABSORCION
•..... ,...... Tallil TllllT j
C
j
C
I
j
Fig. N.6. Intercamhios energencos con cortinas interiores
580
Apendice N
La mayor incertidumbre que se tiene para detenninar el intercambio energctico producido es la cantidad de energia reflejada, absorb ida y transmitida por Ia cortina. a i-Coer. absorci6n clcmento interior 't i-Coer. transmisi6n elemento interior Pi - Coer. rcflcxi6n elemento interior
recordando :
(N.14) Una vez estimadas dichas propiedades tt~nnicas podemos obtener la energia radiante total introducida al ambiente mediante la suma de la serie :
(N.IS)
de otra parte Ia energia absorb ida tanto por el cristal como por Ia cortina y que se transfiere al interior por convecci6n - conducci6n se estima por :
I TRconv-conJ...ll , y )=
se ha utilizado e1 coeficiente 0,25 en vez de 0,35, ya que la cortina impone la existencia de una camara de aire semi-ventilada, la eual aporta una eierta resisteneia tenniea al paso de calor por condueci6n - convecci6n desde el exterior.
581
Manual de climatizaci6n. Torno II
La cncrgia total transmitida sera suma de ambos cfcetos
EI factor "nv" se obtcndnl por :
nv
1TR
cort.-inte. ( 11 , Y)
(N.l8)
I TR(11 , y) sicndo: - Energia total transmitida con vidrio simple (ce. 7.59) ITR cort.-intc. (11, y) -
Energia total transmitida con vidrio simple y cortinas interiores (ec. N.l7)
Sustituyendo y aproximando (se considera que el cristal presenta tanto para la radiacion directa como para la difusa el mismo coefieiente de transmisi6n 't c) :
+
0,25
(X.
(N.19)
nv 't
c
+
0,35
(X.
Y el porcentaje de energia que es transmitida en forma radiante sobre el total :
mr
100
1TR 1TR
raj. 11 , Y)
100
(N.20)
cort.-inte. ( 11 , Y)
1 - Pi P c 582
Apendice N
La caracterizaci6n de las propiedades termicas de las cortinas utilizadas es funci6n basicamente del color (coef. reflexi6n del hilo), y del espaciado del mismo '11 confomlar 1'1 tela (coer. transmisi6n directo). Una vez estimados con exactitud dichos panlmetros tendremos que establecer la distribuci6n fisica de los hilos al conformar la tela, 1'1 direcci6n de 1'1 ra,diaci6n recibida, y suponiendo reflexi6n difusa poder, establecer definilivamente cual es para la lela (6 cortina) los coeficientes globales de transmisi6n, reflexi6n y absorci6n. Ante el desconocimiento ponnenorizado de detalles de tejidos y telas en una cortina determinada, se asume 1'1 siguiente clasificaci6n standard : - En funci6n del color del hilo · Clams (coer. reflexi6n hilo = 0,6) · Medias (coer. reflexi6n hilo = 0,4) · Oscuras (coer. reflexi6n hilo = 0,2) - En funci6n del espaciado del tejido · Espaciada (se observan con nitidez los detalles a su traves) · Media (se ven objetos pero no detalles) · Tupida (los objetos no son visibles) En la tabla N.4 se facilitan tanto las propiedades radiantes asignadas a cada cortina concreta, como el factor "nv" y "mr" aplicando las ec. N.19 y N.20, con cristal de 4 mm.
('t cristal= 0,855
ex * = 0,067
Pc.cristal= 0,04)
583
Manual de clirnatizaci6n. Torno II
CORTINAS
COLOR Claro
I
Medio
I
Oscuro
Espaciada
'tj = 0,57 Pj = 0,35 (Xj = 0,08 nv = 0,67 mr=84%
'tj=0,5 Pi = 0,24 (Xj = 0,26 nv = 0,77 mr=64%
'tj = 0,44 Pi = 0,13 "j = 0,43 nv = 0,87 mr=49%
Tela Media
'ti = 0,42 Pi = 0,46 tti = 0,12 nv = 0,56 mr=75%
' 't i 0,3 Pi = 0,32 (Xi = 0,38 nv = 0,69 mr=42%
'tj = 0,22 Pi=0,17 " j = 0,61 nv = 0,84 mr 260/0
Tupida
'tj = 0,22 Pj = 0,57 " j = 0,21 nv = 0,45 mr=48%
'tj=0,15 Pi = 0,38 ttj = 0,47 nv = 0,64 mr=23%
'tj = 0,08 pj=0,19 (Xj = 0,73 nv = 0,82 mr= 10%
TABLA N.4. Valores de "nv" y "mr" para cortinas
N. 2. 2. Persiallas Hay basicamente dos tipos de ventanas utilizadas en: Ia constnlcci6n, las enrrollables y las del tipo venecianas (laminillas con inc1inaci6n graduable). Las primeras se usan siempre al exterior, y en caso de ser utilizadas no dejan pasar radiaci6n solar en el tramo de puerta 0 ventana que abarcan ..Por tanto en ese supuesto Ia radiaci6n solar transferida debera eonsiderarse nula, aunque no la cantidad de calor transmit ida, pues se comportara como un muro sin inercia tetmica con un determinado coeficiente global de intercambio (recordar que existe una camara de aire entre In propia persiana y el cristal). En general para estimaei6n de la carga termica se supone que dichas persianas estan enrolladas y por 10 tanto es como sino existieran. 584
Apendice N
Las de tipo venecianas se pueden colocar tanto al exterior como al interior, siendo regulable el angulo que fonnan las laminillas conla horizontal. Una regulacion adecuada sera establecida por el propio usuario, de fonna que la radiacion directa sea completamente interceptada, transmitiendose por tanto unicamente radiacion difusa, y siendo esta, funcion del coeficiente de reflexion del material, del angulo de inclinacion de las laminas ("'), y del angulo de incideneia de la radiacion solar con respecto a la ventana (0), ver fig. N. 7. El angulo de incidencia de la radiaeion solar con respecto a la "entana se obtiene mediante (ver ee. 7.76) :
tg
a
(N.21)
donde: hs - Altura solar (0) IBvl - Diferencia entre el aeimut solar y el acimut de la superficie acristalada (0). En valor absoluto
REALIDAD
SIMPLIFICACION
\
\
a,
Fig. N. 7. Intercambios energeticos con persianas venecianas.
585
Manual de climatizaci6n. Tomo II
Como deducimos se pueden presentar una gran variedad de coeficientcs de transmisi6n, reflexi6n y absorci6n. Ante ello se toman unos valores standard en funci6n del color de las laminillas : clara reflectante, color medio no refiectante 6 oscura no reficctante. Clara-reOectante
Media-no reflectante
,;·=015 I ' PI· = 050 ' ex·1 = 0 ' 35
Oscura-no reflectante
,;·=012 I ' PI· = 0' 35 (X·=047 I '
,;·=010 I ' Pi = 0,20 (Xi =0,70
N.2.2.1. Persianas yenecianas al exterior EI comportamiento respecto a la energia total transferida al interior difiere en funci6n de la posici6n de la persiana, asi si se encuentran al exterior In energin absorb ida por las mismns no nlcanzani ellocru, y imicamente incidini en el cristal radiaci6n difuso, (ver fig. N.8).
persiana
cristal 55% • .
..... 45%
Fig. N.B. Intercambio calor con persianas venecianas exteriores 586
Apendice N De acuerdo con cl esquema de la figura N.S podemos obtener la energia radiante incidente sobre el cristal como :
radiacion que sera toda ella en forma difusa
(N.22)
A la que dcbcrcmos aplicar todo el razonamiento establecido en ei apartado 7.2.2.2.1, resultando una cnergia total transmitida :
(N.23)
Como siempre cl factor "nv" se obtiene por :
nv
I TR
pers.-ext. ( 11 , Y)
(N.24)
ITR(ll, y) y la proporcion de energia radiante sera:
mr
(N.23)
Sustituyendo las propiedades radiantes de las persianas consideradas, y utilizando los valores standard : 'to
= 0,855
a * = 0,067
Y
Pc = 0,04 587
Manual de climatizaci6n. Torno II
sc obticnc: Clara-reOectante
Media-no reOectante
nv = 0,14
nv = 0,11
nv = 0,09
mr=96,2%
mr=96,2%
mr= 96,2%
Oscura-no reflectante
N.2.2.2. Persianas yenecianas a1 interior Si la persiana se coloca en el interior la cantidad de energia absorbida por la misma se convertira rapidamente (inercia nula) en energia transferida por conveccion, y sobre el cristal incidini tanto radiacion directa como difusa, ver fig. N. 9.
cristal
persiana
Fig. N.9.lntercambio calor con pers;anas venecianas interiores
De acuerdo con el esquema de la figura N.9 podemos obtener la energia radiante que atravicsa cl cristal y persiana por :
588
Apendice N en la que se supone id6ntico coeficiente de reflexi6n para la radiaci6n direct a y difusa, realmente a partir de la primera reflexi6n en la persiana se considera toda la radiaci6n difusa). Sumando Ia serie :
La energia transferida al interior por conducci6n - convecci6n se compone de dos tenninos, el primero la energia absorbida por la persiana (se supone inercia t<~nnica nula), y el segtUldo sera lUla parte de la energia absorbida por cl cristal. Al fonnarse una camara de ai,-e semi-ventilada entre la persiana y el propio cristal se modifican cI coeficiente de convccci6n interior del cristal, 10 que hani variar el porcentaje de energia que se transfiere al exterior e interior respecto a la existencia (mica del cristal. En este caso concreto supondremos que se transfiere al interior el 25%. Con todo ello la energia transferida por conducci6n - convecci6n sera :
I TR conv-coni Tl ,Y) =['t c I D( Tl , Y) + 0,89
(1 - 't i - P, -
t' j
P , P c - P12 P c -
't t
't 0
I i Tl , Y) ] 2
P,2 P c
-
2
P ,3 P c
- ....... ) +
es dccir:
589
Manu8'1 de climatizaci6n. Tomo II Siendo la energia total transferida :
I TR
perl.-inteJ 11 , y) =
I TR
l'aj..l1 , y) +
I TR
conv-conJ.l1, y )
(N.28)
Como siempre el factor "nv" y "mr" se obtiene por las correspondientes ecuaciones N.18yN.20.
Sustituyendo las propiedades radiantes de las persianas consideradas, y utilizando los valores standard: a* = 0,067
Pc = 0,04
se tiene : Media-no reflect ante
Oscura-no reOectante
nv = 0,50
nv = 0,65
mr=28,8%
mr= 17,9%
nv = 0,80 mr= 12,2%
Clara-ret1ectante
N.3. Existencia de varios accesorios EI procedimiento correcto ante Ia existencia de varios accesorios seria plantear el intercambio existente y posteriormente evaluar los factores "nv" y "mr" como hemos hecho hasta este instante. No obstante si despreciamos Ia multiple reflexi6n entre los accesorios utilizados el factor global "nv" se puede establecer como producto de todos los correspondientes a los accesorios utilizados. Respecto al factor "111r" coincidinl aproximadamente con elmenor de cada uno de los accesorios.
590
Apendice 0 Illercia (Ie los cerramientos. Carga installtanea (lebi(la a aportes energeticos por radiacion
0.1. Allalisis del problema La encrgia aportada por los difcrentes tipos de cargas existentes en aire acondicionado Ioes en divers as fonnas de transmision de cal Of, ver tabla 0.1 :
TIPO DE CARGA (s610 calOl~ sensible)
FORMA DE TRANSMISION
Radiacion Radiacion sup. acrist. sin accesorios interiores * Radiacion sup. acrist. con accesorios interiores * Luces fluorescentes Lamparas de incandescencia Personas ** Infiltracion Ventilacion Maquin'aria y utensiJios
Conveccion
95% 50% 500/0 80% 40%
20-800/0
5% 50% 50% 20% 60% 100% 100% 80-20%
* Se entiende por accesorios interiores cortinas 0 persianas ** Los porcentajes pueden cambiar sustancialmente en funci6n de Ia actividad que desarrolle Ia persona, ver apartado 5.2. TABlA 0.1. Porcentaje (Ie energia radiallte segl'in tipo (Ie cargu.
La energia transmitida en fonna radiante es absorb ida por los cerramientos y objetos del local sobre los que incide; estos al absorber dicha energia elevan su temperatura y posteriormente. por conveccion-conduccion Ia devuelven al local, 10 que se traduce en carga de refrigeracion, (hay que tener en cuenta que parte de dicha energia habra sido transmitida por conduccion a traves del muro). EI proceso no es instantaneo, pudiendo devolver la energia absorbida en el transcurso de varias horas. 591
Manual de c1imatizacion. Torno II
Del proceso expuesto se deduce que no toda la energia radiantc que alcanza cllocal es carga tcrmica en dicho instantc, sino que existira una atcnuaci6n de la misma y un decalaje en cl ticmpo, funci6n del tipo de cerramientos utilizados.
0.2. So!ucioll matematica Para abordar el problema recurriremos al concepto de temperatura radiaci6n-airc al que estan sometidos los cerramientos, que es amllogo al concepto temperntura sol-a ire definido en el apartado 7.2.1.2, y que record amos :
T sa
T se
+
" Ihcr{11,Y)
(0.1)
e
el cual transformamos en :
T ra
TsL
" I {11,Y) T
+
hc,
(0.2)
cuando incide una radiaci6n en un cerramiento desde el interior. Como recordamos del apendice K, el flujo de calor que se transmite a traves de la superficie interior de un cerramiento se puede evaluar mediante la expresi6n K.llS, el cual recordamos : 00
qi(n)
E Tsa(n-j)
YU)
+
i-O
E
T sL(n-j) ZU)
(0.3)
i-O
y aplicando el concepto de temperatura radiante : 00
00
q,(n)
E i-O
Tsa(n-}) YU)
+
E
Tra(n-j) 'ZU)
(0.4)
i-O
Si queremos analizar la transmisi6n de calor existente unicamente dcbida a Ia radiaci6n incidente sobre el cerramiento por el interior del local, deberemos considerar la temperatura sol - aire exterior constante e igual a la interior, por 10 tanto: 592
Apendice 0
q,(n) =
E TsL(n-j) co
~
E IT co
Y(j)
+
~
sL +
(X
1 (n -j) T
~i
1Z(j)
(0.5)
Tenicndo en cuenta la constaneia dc TsL' y recordando las propiedadcs de la sumatoria de los faetores de rcspuesta (cc. K.130) : co
K
co
E
E
Y(j)
j-O
X(j)
-
j-O
E
Z(j)
(0.6)
j-O
Sustituyendo y cancelando tcnninos iguales :
E j-O
(X IT(n-j) - - - Z(j)
hc,
(0.7)
La energia reflejada por un eerramiento concreto sera absorb ida por los demas, aunque es de reeonoeer que los eerramientos en general se pueden eomportar dc forma diferentc (distintos factores Z(j».
0.2.1. Carga por radiacion a traves de superficies acristaladas La radiaci6n solar que atraviesa la sup. aeristalada en forma de radiaci6n ineide fundamcntalmente en el suelo, y este presenta un eoefieiente de absorei6n generahnente alto (entre 0,8 y 0,9), 10 cual impliea que este sera cl eerramiento predominantc, siendo praeticamente insignifieante el difcrente eomportamiento de los dcmas eerramientos dcbido a la esc as a energia radiante que reciben. Por ello se puede aproximar la ee. 0.7 a : ~
qln)
L..J. j-O
I (n-j) 1"-
Z(j)
hc,
(0.8)
donde: Z(j)
.. Coef. "j" de la funei6n de transfereneia "Z" para suclos (W/m2 °C)
IT(n-j) - Radiaei6n ineidente sobre el suelo (W/m2) 593
Manual de climatizaci6n. Torno II
Coef. de conveccion interior, que segilll la norma NBE-CT-74 (rcf.6), para aire acondicionado (flujo de calor ascendente) presenta un \'alor de 11 W/m 2 °C. Flujo de calor absorbido por el cerramicnto (W1m2) EI flujo de calor total absorbido por el suclo se obtendnl multiplicando la anterior expresion por el ntunero de metros cuadrados sobre los que inc ide la radiacion :
lr(n-}) --Zlj)
(0.9)
hc,
donde: Qj(n) - Flujo de calor total absorbido por el cerramiento (W) Ar
- Area de suelo expuesto al sol (m2)
Y recordando que toda Ia radiacion que atraviesa Ia superficie acristalada se evailla por: (0.10) resulta:
Lo
ITr(n-j) ZU)
hc,
j ..
(0.11)
Finalmente la energia que es carga en cada instante sera 141 energia total incidente (suma de Ia energia radiante mas 141 conveetiva conductiva), menos Ia energia absorbida, y teniendo en cuenta eI criterio de signos adoptado para Qi (entrante en el cerramiento <0) :
Qcarga(n) = PT(n) - [ -
~o
mr PT(n-j)
L
j ..
(0.12)
100 he t
donde: PT(n-j) - Energia total que atraviesa 141 sup. mr - Energia en forma radiante (0/0) 594
acristalada~ instante
"n_j" (W)
Apendice 0
0.2.2. Caracterizacion ante illtercambio de calor por radiacion solar a traves de sliperficies acristaladas Dc Ia exprcsion 0.12 podemos conduir que eI panimetro que define In distinta carga aportada en cuanto al comportamiento del cerramientq es la sucesion "Z(j)/hct, cl eual unicamente depende del tipo de suelo utilizado para esta carga.
Se han estudiado 5 tipos de suelos tipicos de eonstnlccion espanola, manteniendo constantes los coeficientes de conveccion interior y exterior en los valores de 11 W/m2 °C y 11 W/m 2 °C respectivamente, segiIn nonna NBE-CT-79 (ref. 6), yo que el flujo de calor en este caso es ascendente (se catienta ellocaI desde el suelo en el que incide Ia radiacion solar) .. rooos los suelos analizados han sido suponiendo la existencia 0 no tanto de parquet como de moqueta, para los cunles se han supuesto los siguientes espesores y propiedades tcnnicas:
m
Conductividad W/m °C
Densidad kg/m3
Calor. esp. J/kg °C
Parquet
0,01
0,14
600
2810
Moqueta
0,01
0,05
1000
1500
Material
Espesor
Se facilitan los datos eonstructivos de los 5 tipos de suelos estudiados, asi como: - Su peso (kg/m2) - K Coef. global de transmision de calor (W1m2 C) 0
SUELO N.t
Peso = 596,0 kg/m2 K
= 2,400 W/m 2 °C
Fig. 0.1
595
Manual de climatizacion. Tomo II Composition
Espesor m
Conductividad W/m °C
Densidad kg/m3
Calor. esp. J/kg °C
0,05 0,20 0,02
1,1 1,63 0,3
2000 2400 800
1380 1050 920
Interior Pavimento Hormigon armado Guarnecido yeso Exterior (piso inferior) SUELON.2
1 em
Peso K
= =
kg/m2
553,0 I,248 W/m2 °C 2em
Fig. 0.2
Composition
Interior Pavimento Impermeabilizante Hormig. ligero de pendiente Bovedilla hormigon normal Guarnecido yeso Exterior (piso inferior)
Espesor Conductividad Densidad Calor. esp. W/m °C kg/m3 J/kg °C m
0,05 0,01 0,10 0,22+0,04 0,02
1,1
0,19 0,35 1,54 0,3
2000 1100 1000 1254 800
SUELON.3
Peso K
442,0 kg/m~ = 2,161 W/m 2 °C
=
2cm
Fig. 0.3
596
1380 1680 1050 1050 920
Apendice 0 Espesor Conductividad Densidad Calor. esp. W/moC m kglm3 J/kg °C
Composicion Interior Pavimento Bovedilla homligon nonnal Guamecido yeso Exterior (piso inferior)
0,05 0,22+0,04 0,02
1,1 1,54 0,3
2000 1254 800
1380 1050 920
SUELON.4 Peso
= 337,0 kglm2
K
=
°,
2
773 W/m °C
2cm
Fig. 0.4
Composicion Interior Pavimento Bovedilla poliestireno Guamecido yeso Exterior (piso inferior)
Espesor Conductividad Densidad Calor. esp. W/moC kglm3 J/kg °C m 0,05 0,22+0,04 0,02
1,1 0,26 0,3
2000 850 800
1380 1000 920
SUELON.S Peso K
= 316,0 kglm2 = 1,156W/m2 °C Fig. 0.5
Composicion Interior Pavimento Honnigon ligero macizo Guamecido yeso Exterior (piso inferior)
Espes 9r m 0,05 0,20 0,02
Conductividad Densidad Calor. esp. W/moC J/kgOC kglm3 1,1 0,35 0,3
2000 1000 800
1380 1050 920
597
Manual de climatizaci6n. Torno II
En la figura P.6 se representan los valores ZG)lhci de los diferentes suelos, (hasta 15 combinaciones con la consideraci6n de moquetas y parquets). Como observamos aunque cada suelo en concreto se comporta de una manera distinta, es evidente que se pueden c1asificar en funci6n de su terminaci6n, ya sea pavimento normal, parquet 6 moqueta.
Coeficiente ZO)/hc j 0,4~------------------------------------------------~
0,2 . - - - - -.- - - - - -
-I
-
-
-
-
-
-
...
-
-
-
-
-
-.-
-
-
-
-
-
...
-
-
-
-
-
-
po
-
-
-
-
-
-0"
-
-
-
--
-0,2 - -Suelo~ -con- ~ta - ~ - - - - - - ~ - - - - -
-0, ,4
<- - - - - -
Suelos con parquet
-0,6
Suelo~ de pa~imento :
-0,8
-1
~------------------------------------------------~
o
1
2
3
4
5
6
7
8
Instante "j"
Fig. 0.6 Caracteristica ZO)lhci de diferentes sue/os
EI conjunto de coeficientes ZG)/hci representativo de los tres tipos de terminaciones de suelos utilizados se facilitan en la tabla 0.2. El suelo concreto tornado como standard es el n03.
598
Apendice 0
N° DE ORDEN
PAVIMENTO
PARQUET
MOQUETA
0
-0)72872
-0,519755
-OJ() I () I G
I
0,140868
0,096267
0,085189
2
0,073345
O,()16010
0,014780
.,'"
0,051081
0,O2()O51
0,010709
4
0,038840
0,020()23
O,O()8690
5
0,031095
0,017104
0,007390
6
0,025811
0,014632
0,006471
7
0,022001
0,012802
0,005784
8
0,019121
0,011390
0,005249
9
0,016850
0,010259
0,004817
10
0,014991
0,009322
0,004456
11
0,013426
0,008525
0,004146
12
0,012078
0,007831
0,003873
13
0,010898
0,007215
0,003629
14
0,009853
0,006663
0,003407
15
0,008919
0,006161
0,003202
16
0,008082
0,005704
0,003013
17
0,007327
0,005284
0,002837
18
0,006645
0,004897
0,002673
19
0,006028
0,004541
0,002519
20
0,005469
0,004211
0,002374
21
0,004963
0,003905
0,002238
22
0,004504
0,003623
0,002110
23
0,004087
0,003361
0,001989
24
0,003709
0,003117
0,001876
25
0,003367
0,002893
0,001769
26
0,003055
0,002683
0,001668
599
Manual de climatizacion. Torno II
27
0,002773
0,002490
0,001573
28
0,002517
0,002J 10
0,001483
29
0,002284
0,002143
0,001398
30
0,002073
0,001988
0,001319
31
0,OQ1882
0,001845
0.001244
32
0,001708
0,001712
0,001173
33
0,001550
0,001588
0,001106
34
0,0014.07
0,001473
0,001043
35
0,001277
0,001367
0,000983
36
0,001159
0,001268
0,000927
37
0,001052
0,001177
0,000875
38
0,000955
0,001092
0,000825
39
0,000866
0,001013
0,000778
0,000~40
0,000733
,
40
0,000786
41
0,000714
0,000872
0,000691
42
0,000648
0,000809
0,000652
43
0,000588
0,000751
0,000615
44
0,000534
0,000696
0,000580
45
0,000484
0,000646
0,000547
46
0,000440
0,000600
0,000516
47
0,000399
0,000556
0,000486
ReI Z(i+ 1)/Z(j)
0,9074623
0,9276434
0,9430648
TABlA 0.2. Coef. ZOJ//tcipara los Ires t!pos de termillaciones de slle/os.
Algunos autores suclcn indicar pequeuas diferencias de comportamiento en funcion del peso del suc1o, aunque de 10 observado podemos afirmar que Ia incertidumbre en la determinacion del coeficiente de conveccion, 0 incluso las propiedades tcmlicas de parquets 0 moquetas, conlleva mucho mayor error que la precision que se quiere ganar con la anterior distincion de tipos de suelos. 600
Apendice 0
0.2. 3. Carga por radiacion desdefocos interiores Dentro de este apart ado se encontranin todos los demas tipos de cargas comentadas (luces fluorescentes e incandesccnles, personas y maqllinas). En este caso se puede realizar el mismo razonamiento anterior, pero hay que tener en client a que la energia radiante se considera difllsa, y por 10 tanto incidira en todas las superficies interiores del local. Para un calclilo riguroso deberiamos conocer las dimensiones del mismo y la posicion relativa del foco de calor, con el fin de detenllinar los factores de fonna de las distintas superficies. Una primera aproximacion podria establccerse suponiendo que un 30% sera absorbida por el suelo, un 30% por el techo y un 40% por las paredes. Con estos porcentajes la carga tcnnica definitivamente se obtendria mediante:
+" 0,3 mr100PT{n-j) Z he 00
Q
(n) =PT{n)
carga
L.J
j ..o
sIlelo
u)
+
i sIlelo
(0.13)
+;... L.J j-O
u) +;...
0 3 mr PT(n-j) Z 0,4 mr PT(n-j) Z (i) ' 100 h techo.l..J 100 he parerJV C t techo
i-O
t pared
0.2.4. Caracterizacion ante intercambio de calor por radiacion desde focos internos Dc Ia expresion 0.13, concluimos que el parametro que define el comportatlliento del local es la sucesion Z(j)/hCj para los distintos cerramientos del local, y ponderados por la proporcion de cerramiento que es "vista" desde la fuente de calor, y que se ha considerado un 300/0 el suelo, un 30% el techo y finalmente un 40% las paredes.
EI comportamiento del suelo ya ha sido estudiado en el apartado anterior, y se ha comprobado que no depende pnicticamente del peso 0 constitucion del mismo, sino de su acabado interior. Respecto al suelo seguiremos adoptando la clasificacion anterior: Tenninacion de Pavimento, Parquet, 0 Moqueta-alfombra.
601
Manual de dirnatizaci6n. Torno II
Las paredes han sido estlldiadas en el apcndice K y L, y en ellos podemos comprobar que para el cocficicnte YO) del muro era mll)' importante la inercia del mismo (cs csla la clasificacion adopt ada para excitacion extcrior), mientras quc los 20) sc manticnen pniclicamentc constantcs cuando el acabado interior es el mismo (vcr apartado K.3.4.2), en cl cllal para las tres parcdes consideradas tipo (en cuanto a coer. yO», se observa pnlclicamentc el mismo valor de 20). Los acabados interiores de paredcs son nomlalmente en enlucido de yeso, 0 en madera, (olros tipos de acabados, tales como telas, homligon~ etc ... aunque cxisten son atipicos). Con estos dos tipos de acabados y Ia pared considerada ligera (muro n07), se obtienen los siguientes coeficientes 2U) :
N° DE ORDEN
Tipo aeabado
0
1
2
3
Enlucido yeso
-5,038222
1,728167
0,9829197
0,6514493
Madera (1 em)
-4,822867
1,542535
0,6903011
0,5173692
TABLA 0.3. Clltltro primeros coef. ZO) para diversos acahados de paredes
Como observamos en Ia tabla 0.3 Ia variacion es tambicn escasa para los acabados interiores considerados "standard". Los tcchos fueron igualmente estudiados en los apendices K y L y podemos extraer la misma conclusion que con las paredes respecto a Ia nula importancia de la inercia de Ia misma respecto al coer 20), y por el contrario una incidencia real respecto at acabado intcrior. Los acabados interiores de techos son normalmente en enlucido de yeso, en cielo raso, o en madera, (olros tipos de acabados, tales como leIas, honnigon, etc ... aunque existen son atipicos). Con estos tres tipos de acabados y el techo mas representativo de Ia constntccion espanola (n° 9 sin aislamiento) se obtienen los siguientes coeficientes ZU):
602
Apendice 0
N° DE ORDEN
Tipo aeahado
0
1
2
3
Enlueido ycso
-4,02608
0,5448628
0,2283137
0,16318
Mader~l
-3,775254
0,789075
0,2189659
0,1333951
-3,324149
1,14616
0,1897238
0,077591
(l em)
Cicio raso
TABlA 0.4. Cuatro primeros coef. ZO) para diversos «cab at/os de tee/lOs
Como observamos en la tabla 0.4 la variaci6n es tambien esc as a para los acabados interiores considerados "standard". Con Jas considcraciones anteriores podemos resumir : - EI rctardo a considerar para cargas intemas (luces, ocupantes, etc .. )~ vendni dado
por Ia porcion de energia disipada en forma radiante "mr", y por el tipo de acabado del suelo. Siendo de segundo orden el acabado de techo y paredes, asi como la inercia termica del tipo de construccion. • Respecto al tipo de acabado de suelos se toman el pavimento (terrazo nOffilal), el parquet (madera), y Ia moqueta (0 alfombras). • Tanto para paredes como para techos asumiremos el acabado en enlucido de yeso por ser este el mas frecuente. EI conjunto de coeficientes
L1'-0
0,3
Z sueloU ')
he i
suelo
+L / .. 0
0,3
V)
Ztecho
he t
techo
+L / .. 0
0,4
ZpareP)
he i
(0.14)
pared
representativo de los tres tipos de terminaciones de suelos utilizados se facilitan en Ia tabla 0.5. 603
Manual de climatizacion. Torno II
N° DE ORDEN
PAVIMENTO
0
604
PARQUET
MOQUETA
-0,663631
-0,575207
-0,507637
1
0,151677
0,134759
0,130558
2
0,080885
0,066731
0,058671
3
0,054954
0,045460
0,039641
4
0,039050
0,032140
0,027614
5
0,028683
0,023376
0,019692
6
0,021722
0,017482
0,014386
7
0,016939
0,013449
0,010787
8
0,013576
0,010643
0,008314
9
0,011151
0,008651
0,006587
10
0,009356
0,007206
0,005360
11
0,007990
0,006131
0,004470
12
0,006920
0,005309
0,003808
13
0,006060
0,004663
0,003303
14
0,005353
0,004143
0,002908
15
0,004758
0,003712
0,002590
16
0,004251
0,003349
0,002329 .
17
0,003811
0,003036
0,002108
18
0,003426
0,002764
0,001920
19
0,003087
0,002522
0,001756
20
0,002785
0,002307
0,001611
21
0,002515
0,002114
0,001481
22
0,002273
0,001939
0,001365
23
0,002056
0,001781
0,001260
24
0,001861
0,001636
0,001165
Apendice 0
25
0,001684
0,001504
0,001078
26
0,001525
0,001384
0,000999
27
0,001381
0,001274
0,000926
28
0,001251
0,001173
0,000859
29
0,001133
0,001080
0.000797
30
0,001027
0,000995
0,000741
31
0,000930
0,000916
0,000688
32
0,000843
0,000844
0,000640
33
0,000764
0,000778
0,000595
34
0,000692
0,000718
0,000554
35
0,000627
0,000662
0,000516
36
0,000569
0,000610
0,000481
37
0,000515
0,000563
0,000448
3'8
0,000467
0,000519
0,000418
0,000423
0,000479
0,000390
40
0,000384
0,000442
0,000363
41
0,000348
0,000408
0,000339
42
0,000315
0,000376
0,000317
43
0,000286
0,000347
0,000296
44
0,000259
0,000321
0,000276
45
0,000235
0,000296
0,000258
46
0,000213
0,000273
0,000242
47
0,000193
0,000252
0,000226
39
ReI Z(j+ 1)/Z(j)
.
0,850767
0,856822
0,861448
TABlA 0.5. Coef. ZOJ1hcipara!ocos internos de calor por radiacion.
605 ,
ApendiceP Cliicllio de cargas termicas en cllle/accioll con varillcioll Iloraria de temperattlras
P.l.lntroduccion EI estudio de cargas tennicas para refrigeracion se ha establecido en base a una cierta oscilacion de las variables ambientaIes, tal como temperatura seca, radiaci6n solar, etc ... , en cambio en la determinaci6n de la potencia necesaria para calefacci6n se ha aconsejado tomar unas condiciones constantes de dichas variables (temperatura seca constante e igual a Ia minima y radiacion solar nula). Se plantea pues Ia necesidad de justificar dicha forma de actuar y cuantificar en la medida de 10 posible Ia desviacion que clio diera Iugar.
P.2. Variacion aconsejada de las variables ambientales En el apendice G se discuti6 como obtener Ia evoluci6n de estas variables para las condiciones criticas de inViemo, (es de sefialar que faHan estudios y nonnativas al respecto para una generalizacion de dichas tluctuaciones). Asi con rcspecto a la temperatura seca se proponia una oscilacion de la forma:
T,
min =
O~
[1
+
cos (FII (HS
+
F3)
+
F2 1)
1
+
T,
miN IN.
(p.I)
siendo: OmO
- Oscitacion temperatura diaria minima (OC), se aconseja para Valencia OmO = OMO/2, donde OMD es Ia oscilaci6n media diaria (tabla 6.1)
Ts min inv
-
HS
- Hora solar (h)
Fl j
-
Temperatura seca minima de proyecto en inviemo con un determinado nivel percentil (tabla 6.7)
Funciones defmidas en el apendice G (funcion de Ia hora solar y Ia hora de salida del sol) 607
Manual de climatizac;ion. Torno II La humedad rclativa maxima en inviemo se considera en Valencia del 95% (Ia cual se produce a la hora de salida del sol); como siemprc faltan estudios para poder extrapolar a otras ciudades de Espana. Se considcra constante la humedad especifica a 10 largo del dia, tal y como se senalaba en Ia fig. 6.2. Con estos condicionantes se obtuvo para Valencia y el mes de Enero (minima temperatura), tanto la evoluci6n diaria de la temperatura seca con un nivel percentil del 97,5%, como la humedad relativa, tal y como se representa en Ia fig. G.6 y que a continuaci6n se reproduce.
.
Temperatura seca (OC) Humedad relativa (%) 10.-----------------------------------------,100
: :~:~ ~: ~ ~: ~: --: ~ ,:~ ,:>...;..: : ::..::::::.::::. .::_.:-:. :-.. ~:~: ...:: . . . : . . . ~ . . . : . . . :. . . .: . .
7
~~
~
;..~
70
• . . . ' . . . • . . . 60
5
''- .:. . . . ' . . . ' . . . '. . . .' . . . '. . . .' . . . ' . . .
~
4
. . . : . . . ~ . . . : . .4.
,
:....: .
'
3
,
,
' • • •' .
-
'.
~_ _ _...
~
~
... '
.....
•
o
-
•
;. .
2
o
2
4
6
6
t ....
w
...'
•
10
12
.......
I
14
'
..
•
16
~
..
.!
'. . .
J
,
,
•• ' ••
•••
!.
50
30
: . -. : . . . ; . . - 20 ..
:.
. . :. . . : - . . 10
•
16
20
22
0
24
Hora solar
Fig. P.I. Fluctuacion de fa temperatura secay humedad relativa en las condiciones de proyecto para ;nv;erno en Valencia, mes de Enero. (N;vel percentil del 97,5%).
La radiaci6n solar el dia critico de inviemo se considera formada unicamente por radiaci6n difusa, siendo esta de fonna aproximada el 50% de la que Ie corresponderia ese dia de inviemo con cielo despejado. Con estas consideraciones se obtuvieron los valores senalados en la tabla P.l para Valencia, mes de Enero y que se represent an gnificamente en la fig. P.2. Hora solar (h)
8
9 10 11 12 13 14 15 16
Horizontal
22 36 42 45 45 45 42 36 22
Cualquier orientaci6n
13 22 25 27 27 27 25 22 13
TABlA P.I. Radiacion mInima, toda di/usa en Wlm 2 , mes de Enero. Latitud 40
608
0
Apendice P
Radiaci6n solar (w/m2) 50r-------------------------------------------------~ . . . -: . - Radiaei6n :sobre- ~ . . . .: . ~-"""t.
40
_ •..' ..
~u~~r~j~i~ :~~ri~~n;t~1
. . .' . _ . _' ..
.-.~-.:-~-~--
'
30
20
10
- - . -: . - .. : - - .. :- - -
~
- - - -: . - - -: - - - - :- - - -:
. - -: - - - -: - - - - :- - - -
O~------------~--------------------~------------~ 6 10 12 14 16 18 20 22 24 2 4 8 o
Hora solar
Fig. P.l. Fb,ctuacion de la radiacion solar incidente en Valellcia sobre sup. horizontal y vertical el dia de proyecto en invierno
P.3. Temperatllra equivalente en invierno La estimacion de la temperatura equivalente, ya sea para una pared 0 para un techo se deberia estimar en base a la ec. 7.20 0 7.21, es dccir deberemos conocer en primer lugar la temperatura sol-aire.
Tsa
(P.2)
En la fig. P.3 se representa la temperatura sol-aire para techos oscuros (a = 0,9 y hCe = 20 W/m2°C) y para paredes oscuras (a = 0,9 y hCe = 16,67 Wltn2 °C). Para colores medios 0 elaros estaria entre la evolucion representada y la temperatura seca exterior.
609
Manual de climatizaci6n. Tomo II
Temperatura (OC) 10r-----------------~-----------------------------------,
9 - - - -:- - - - ~ - - - -Jemp.-solaifetecho ':-; •. "::_0_0 .• - -:-.
8
F -:::---: :::'--:-
: :_::': :-, ~ :::::- -::-::: >;~.;.; <" ~.\ ,
,
::,'
,
,
- -; - ... :- - - -.: - ...
"
"
5 4
3 2 I
1
·1
0
0
2
4
6
8
10
• ..
12
....
;'''
~
14
..
-,.
16
I ..
•
..
-:
1 ..
18
..
..
". ~
20
22
24
Hora solar
Fig. P.3. Fluctuacion de fa temperatura sol-aire en Valencia sobre superficie horizontal y vertical el dia de proyecto en invierno
Para dctcnninar la temperatura equivalente (ee. 7.21), dcberemos dcterminar el comportamiento de los cerramientos, los cuaIes como sabemos (apcndice N) se pueden clasifiear en tres grupos, ligeros, medios y pesados.
T
seq
;... T (n-j) yU) L..J J-O
SQ
K
(p.3)
Aplicando la cc. P.3 sobre los cerramientos standard podemos obtener Ia temperatura equivalcntc, la cual se representa en las fig. P.4 para paredes y P.5 para techos. En ellas se reflcjan asimismo las temperaturas equivalentes a considerar con radiacion solar nuia (lineas a trazos).
610
Apendice P Paredes oscuras
Temperatura (OC)
10~----------------------------------------------~ :
8
temp. sol"aire par~des
" , / " ' " "Titlriperlltura sei::a" ainDlente
6
4 Pare~ pesadll:
2 Pared: ligera
O~-----------------------------------------------------J
o
4
2
6
8
10
12
14
16
18
20
22
24
Hora solar (h)
Fig. P.4. Temperatura equivalente para paredes standard oscuras en Valencia, con y sin contar el aporte debido a La radiacion solar Techos oscuros
Temperatura (oG)
'O~--------------------------------------------------~ SOI:"aire tec~os
8
-""".
-
. - .
P~red pes~da .
,
~
,
,
temperafura seca "amtiientti "
Pared media ,
: ; Ji i ;
4
~ ; : : :: : :
::1 ! .
~ ~ -
~. .
! :':'. -: .. , : t,.'
! .' .::.: .: .: ::: .: .:
::!.:.:.~
.'
2 .
o~--------------------------------------------------~
o
2
4
6
8
10
12
14
16
18
20
22
24
Hora solar (h)
Fig. P.S. Temperatura equivalente para techos standard oscuros en Valencia, con y sin contar el aporte debido a La radiacion solar
611
Manual de climatizacion. Torno II
Oc 141 obscn'acion de las figuras PA y P.5 podemos conc1uir : 1) La temperatura equivalente exterior varia a 10 largo del dia, sobre todo en paredes ligcras. 2) La minima temperatura exterior equivalente a considerar no se produce a la hora de salida del sol, sino unas horas mas tarde dependiendo de 141 inercia del muro. 3) La importancia relativa de la radiacion solar sobre 141 temperatura equivalente a considerar es muy escasa, pudiendo pnicticamente considerarse nula. 4) La corrcccion a establecer para obtener Ia temperatura equivalente minima deberia eontemplar la temperatura minima exterior mas un porcentaje de la oscilaeion minima diaria COmO), siendo dieho poreentaje funeion de la inercia del cerramiento y. valido para eualquier localidad Cobservar las ee. P.I, P.2 y P.3 con radiacion solar nula). De 141 observacion de las figuras PAy P.5 se conc1uye Pared ligera Pared media Pared pesada Techo
Tseq,mfn = Ts,min,inv + 0,20 OmD Tseq,min = Ts,min,inv + 0,34 OmD Tseq,min = Ts,min,inv + 0,45 OmD Tseq,min = Ts.min,inv + 0,43 0111D
5) La temperatura equivalente minima a 141 hora de salida del sol se puede estimar mediante: Pared ligera Pared media Pared pesada Teeho
Tseq,min = Ts,min,inv + 0,40 OmD Tseq,min = Ts,min,illv + 0,52 OmD Tseq,mfn Ts,mfn,inv + 0,54 OmD Tseq,min = Ts,mfn,inv + 0,55 OmD
6) En comparacion con el procedimiento "normal" de ealculo de eargas (se toma como temperatura la minima de proyecto Ts,min,inv), afirmamos que 141 earga por transmision a traves de paredes y techos esta mayorada, siendo el error relativo cometido 612
Apendice P
(~
erc1ativo
)
100
A K (T .
i + Por OmD - T ,T) - A K (T . t - T ,T) s,nun. ny SL ',nun. ny 'L A K (T . t + Por OmD - TNT) ',nun. nv
=
0
:i.L
e relativo (%)
Por O,nD
100
T I, nun,"my
+
Por 01nD - T sL
(P.4)
donde: ercJativo Por OmD
- Error relativo en ~{,. - Por~entaje de Ia OmD a contabilizar (tanto por lIDO) - Oscilaci6n minima diaria (0 C)
Ts.min.inv
-
TsL
-
Temperatura minima de proyecto con un determinado nivel percentil (tabla 6.7 en °C). Temperatura seca del local en inviemo
Para dar un orden de magnitud aplicaremos para Valencia la anterior expresi6n : Por= 0,25 OmD= 5°C Ts,min.inv = 1°C (NP =
97,5%)
TsL = 20°C e relatiYo(%) = 100
Por OmD
= 100
TI,nun,'my • + PorOmD -T' sL
0,25 5 1 +0,25 5-20
7%
Luego Ia carga por este concepto seria un 7% menor, por 10 que su incidcncia en el cnlculo global de cargas es despreciable. .
P.4. Carga por radiacion solar a travis de superficies acristaladas En general como hemos seiialado dicha carga se considera nula, no obstante una cierta radiaci6n solar a 10 largo del dia existe como hemos visto en la tabla P.1. 613
Manual de clirnatizacion. Torno II Dicha energia debcria atravesar un vidrio comun, siendo toda ella di fusa, por 10 que Ia carga seria :
Carga = 0,89
't 0
IT = 0,89 0,855 IT = 0,76 IT
(P.S)
A este valor dcbcriamos posteriormente afectarlc del correspondicnte cocficiente nVi (tabla 7.26) en funcion de los materiales adicionales utilizados.
Carga
=. nV
t
0,76 IT
(P.6)
Por todo 10 cual para ventanas en paredes. estariamos hablando en el caso mas favorable de :
Carga
=
nv, 0,76 IT = 1 0,76 27 = 20 Wlm
2
El error cometido al despreciar este aporte calorifico es practicamente nulo, ya que: 1)' EJ momento de maxima carga se suele producir en las horas de la salida 'del sol, y en esos momentos dichacarga es efectivamente nula. 2) Si tenemos en cuenta cortinas, doble cristal 0 alglm elemento accesorio se nos convicrte en un aporte de unos lOW1m2 de ventana en las horas centrales del dia, valor que es c1aramente despreciable.
P.5. Metodo de calculo El metodo de calculo que podriamos establecer es por 10 demas coincidente con el de cargas de refrigeracion, atmque contabilizando como qued6 dicll0 el minimo de personas, luces, etc, es decir cl minimo de aportes intemos de energia. La discrepa.ncia que podemos finalmente observar no seria superior a un 5% con el metoda dcscrito en el capitulo 7.3, por 10 que no se suele contemp]ar, aunque insistimos una corrccci6n de la temperatura exterior para el calculo de cargas a traves de paredes y techos seria conveniente.
614
Apendice P
ANEXO. PROGRAMA DE ORDENADOR
El paquetc de programas que se facilita con ellibro desarrolla todo el complicado proceso matemalico explicado en los diversos apartados del mismo, con el fin de poder adaptar los resultados a un determinado equipo. A continuaci6n se describen las principales caracteristicas, asi como los datos de partida y resultados obtenidos. - Arrancar con la palabra mpAlRE. - Pennite un almacenamiento en disco duro con copy *.* - En cualquier punto del programa pulsando la tecla ESC el programa vuelve a la pantalla inicial. CARACTERISTICAS DEL PAQUETE DE PROGRAMAS
I
mpAIRE
I
EI paquete se com pone de los siguientes programas :
B> 1 Lectura en el diagrama sicrometrico ~ 2 Transformaciones sicrometricas
~ 3 Estudio detallado de serpentines [> 4 Humectadores 0 torres refrigeracion
D>
~
5 (> 6 [Y 7 B> 8 l> 9
Confort humano EI ambiente exterior Hora civil < .. > Hora solar Tablas cargas Aire Acondicionado AnaHsis trans. calor cerramientos
615
Manual de climatizaci6n. Tomo II
Caracterfsticas del programa ·
C> [>
Se estudian : - Recintos cerrados (temp. uniforme paredes) - Recintos cerrados (una sup. a temp. muy distinta) - Recintos abiertos Datos de entrada: - Constituci6n ffsica de la persona - ripo, color y cantidad de vestidos - Actividad desarrollada por la persona - Temperatura seca aire (OC) - Humedad relativa aire (%) - Temperatura de las paredes (OC) - Velocidad del aire (m/s) Otros datos para recintos abiertos : - Radiaci6n solar (w/m2) - Coef. absorci6n vestidos a long. onda corta - Temperatura suelo (OC) - Factor de forma persona - cielo Otros datos recinto cerrado (una sup. a temp. diferente) - Temperatura de esa superficie (OC) .. Area de esa superficie (m 2) - Factor de forma persona - superficie caliente - Coef. absorci6n esa superficie
[> 616
Datos de salida: - Cantidad y concepto de calor intercambiado - Indices termicos del ambiente - Temperatura efectiva e fndice de esfuerzo termico - Grado de confort
Apendice P
Caracterfsticas del programa :
Se puede calcular : - Posicion del sol (altura y azimut) - Radiacion solar maxima (segun modelo C de Iqbal) - Evolucion temp. seca y humedad relativa (Norma Une 100-014-84)
C>
POSICION DEL SOL ..
Datos de entrada: - Latitud del lugar Datos de salida : - Altura y azimut solar para todos los 21 de Mes y para todas las horas solares
C>
EVOLUCION TEMP. SECA V HUMEDAD RELATIVA (NORMA UNE 100-014-84) Datos de entrada: - Variacion anual de temperatura (OC) - Temp. seca maxima proyecto (OC) - Variacion diaria temperatura (OC) - Temp. humeda de proyecto (OC) - Cota sobre fa estaci6n (m) metereol6gica Datos de salida: - Temperatura seca y humedad relativa para tad as los meses y horas sol res 617
Manual de clirnatizaci6n. Torno II
Caracterfsticas del programa ·
SOLAR MAXIMA [> RADIACION (MODELO "CII DE IQBAL) Datos de entrada: - Altura sobre el nivel del mar (m) - Tipo de atm6sfera · Elegir una 0 imponer coef. anstrong - Coef. reflexi6n alrededores · Elegir alrededores · Imponer uno - Cantidad de vapor de agua en el aire · Imponer uno · Calculado por programa (en base a la temp. seca y humedad relativa, estimada en base al anterior programa) _ Latitud (0) Datos de salida: - Radiaci6n solar maxima para todos los meses, todas las horas solares y orientaciones posibles
[> Trabaja en Kcal/h m
2
0 W/m
2
[> Se puede especificar caulquier mes y orientacion determinada de calculo
618
Apendice P
Caracterfsticas del programa :
/7
Hora civil
<---->
Hora solar
[> Convierte la hora civil en hora solar
0
I
viceversa
para cualquier localidad
[> Datos entrada para pasar de hora civil a hora solar: - Hora civil (Horas, minutos y segundos) - Longitud del lugar (Grados y minutos) - Longitud de referencia del pais (Grados) (Para Espana -30 0 Verano y -15 0 Invierno) - Mes y dfa del ano
[> Datos entrada para pasar de hora solar a hora civil : - Hora solar (Haras, minutos y segundos) - Langitud del lugar (Grados y minutos) - Langitud de referencia del pais (Gradas) (Para Espana -30 0 Verano y -15 0 Invierno) - Mes y dfa del ana
619
Manual de clirnatizaci6n. Torno II
Caracterfsticas del programa : 8 Tablas cargas Aire Acondicionado
t> t>
Se desarrolian todas las tablas necesarias para calcular las cargas en aire acondicionado. Siendo los programas implementados :
[> Estimacion
cargas solares: . · Energfa maxima transmitida a traves de un cristal · Carga instantanea para el equipo de aire acondicionado por radiacion solar a taves de un cristal
l? Cargas internas teniendo en cuenta la inercia del local · Personas · Luces .Otras
[> Cargas por transmision
de calor en cerramientos. Obtencion de la temperatura equivalente, tanto para techos como para paredes.
620
Apendice P
Caracterfsticas del programa : 8 Tablas cargas Aire Acondicionado ENERGIA SOLAR MAXIMA 0 CARGA POR RADIACION SOLAR
[> [> .. [> [>
Trabaja en Kcal/h 6 W Se estima para las orientaciones basicas 6 sobre cualquier orientaci6n especificada Se calcula para todo el ano
0
para cualquier mes
Datos de entrada: - Altura sobre el nivel del mar (m) - Tipo de atm6sfera · Elegir una 0 imponer coet. anstrong - Coet. reflexi6n alrededores · Elegir alrededores · Imponer uno - Cantidad de vapor de agua en el aire · Imponer uno · Calculado por programa (en base a la temp. seca y humedad relativa, estimada en base al programa sobre el ambiente exterior _ Latitud (0) . - Tanto por cien de hueco ocupado por el cristal (*) - Espesor del cristal (mm) (*) - Coet. transm. a radiaci6n perpendicular (93.3% cristal normal) (*) - Indice de refracci6n (*) (*) Se puede indicar cristal normal y tipo de carpinteria
En caso de estimar la carga se debe ariadir : - Existencia 0 no de persianas 0 cortinas interiores - Tipo de construcci6n (elegir standard, 3 tipos, o especificar uno determinado) - Horas de funcionamiento del sistema de aire acondicionado
621
Manual de climatizacion. Torno II
Caracterfsticas del programa : 8 Tablas cargas Aire Acondicionado INERCIA CARGAS INTERNAS
Se estima el % de energfa que es carga para el equipo de aire acondicionado
[>
Datos de entrada: - Tipo de carga a considerar .Ocupantes · Luz incandescente · Luz fluorescente · Otras (especificar % energfa radiante) - Tipo de cerramientos : · Standard (3 modelos) · Especificar uno cualquiera determinadado (observar datos requeridos en el programa de analisis trans. calor cerramientos regimen cualquiera) - Horas de funcionamiento del equipo de aire acondicionado
622
Apendice P
Caracterfsticas del programa : 8 Tablas cargas Aire Acondicionado TEMPERATURAS EQUIVALENTES
[> [> [> [> [> [>
Se calcula tanto para techos como paredes Para todas las orientaciones 0 una determinada Para todos los meses 0 uno concreto Para cerramientos standard 6 un cerr~miento especffico (ver en este caso datos de entrada en el programa de analisis trans. calor cerramientos regimen cualquiera) Para 3 colores superficie exterior, claro, medio y oscuro Datos de entrada: - Datos referentes a la temperatura exterior (ver programa sobre el ambiente exteri,ar) - Datos de radiaci6n solar (ver programa sobre el ambiente exterior) Resultados: - Temperatura equivalente a considerar eC)
623
Manual de climatizacion. Torno II
Caracterfsticas del programa ·
Se estudian los siguientes procesos : - Regimen estacionario - Regimen senoidal puro exterior (24 h) - Regimen cualquiera (Func. de transferencia) Los resultados se puden obtener en : W 6 Kcal/h
[>
Los datos requeridos son: 2 -Coet. convecci6n exterior e interior (W/m K) - Espesores distintas capas (m) - Conductividades capas hW/mK) - Densidades capas (Kg/m) (*) - Calor especffico capas (J/Kg K) (*) 2 - En caso capa de aire la resistencia termica (m K/W)
(*) No requerido en regimen estacionario
[>
624
Los resultados obtenidos son : - Regimen estacionario · Coef. global - Regimen senoidal puro · Coef. global · Amortiguamiento · Desfase · Peso - Regimen cualquiera · Coef. global · Amortiguamiento · Desfase · Peso · 48 coef. funci6n transf. X,Y,Z · Relac6n entre los ultimos coef.
Terminos principales en la nomenclatura utilizada Signos latinos A AH As
- Area superficial (m2) - Angulo horario (0) - Acimut solar (0)
B CA
- Coeficiente traspaso de masa (kg/s Pa m2) - Calidad de aire (decipoIs)
Cp Cf
- Calor especifico (kJ/kg °C 6 kcal/kg °C) - Calor de cambio de fase (liquido-vapor,kJ/kg 6 kcal/kg)
c
- Ratio comun de Ia funci6n de transferencia de un cerramiento
D
- Diametro tuberia (m)
DEC - Declinaci6n C) E - Energia (J 6 kcat) Eo
- Correcci6n para Ia excentricidad de la 6rbita solar
esp FS GD HCR HS
- Espesor (m) - Factor 6 Coeficiente de simultaneidad - Grados-dia - Hora civil (It) - Hora solar (h)
HSS - Hora de salida del sol (b) h
- Entalpia (kJ/kg a.s. 6 kcal/kg. a.s.)
hc
- Coeficiente de convecci6n (W/m2 °C 6 kcal/h m2 °C)
ho hs
- Coeficiente de traspaso de masa (kg/s m2 . kg/kg as) - Altura solar (0) - Radiaci6n solar (W1m2 6 kcallh m2) - Angulo de incidencia (0)
I
K Kc
- Coef.global de transmisi6n de calor (W/m2 °C 6 kcaIIh m2 °C) - Coeficiente de extinci6n (cm- I )
KI j
-
Relaci6n entre carga instantimea y potencia temlica, i= 1 (Ocupantes) i=2 (Hum. incandescente) i=3 (Hum. fluorescente) 625
L
- Longitud (m)
LA T - Latitud C) Le - N° de Lewis MN - Energia neta/unidad tiempo producida por las personas (W 6 Met) m - Masa (kg) mr - Porcentaje de energia en forma radiante (tanto por cien) ND - N° del dia segiln el ano Juliano np - N°de personas 0 ocupantes nr - Indice de refraccion nv - Coeficiente de materiales adicionales para acristalamientos P - Presion (Pa) PT - Potencia termica (W 0 kcallh) Q - Potencia calorifica (W 6 kcalJh)· Qp q r S ST T TSV t V v ve W XU) VU)
- Potencia calorifica disipada por persona (WIper. 0 kcallh per.) - Potencia tcrmica por unidad de superficie (W1m2 6 kcallh m2) - Radio (m) - Superficie (m2) - Superficie exterior total (m2) - Temperatura (OC) - Tiempo solar verdadero (h) - Ticmpo (h 6 s) - Volumen (m3) - Vclocidad (mls) - Volumen especifieo del aire hu.medo (m3/kg a.s.) - Humedad especifiea (kg/kg a.s.) - Coer. "X" de la func. de transferencia (W/m 2 °C 6 kcallh m2 °C) - Cocf. "V" de la fune. de transfereneia (W/m2 °C 0 kcallh m2 °C) ZU) - Coef. "Z" de la func. de transferencia (W/m2 °c 6 kcallh m2 °C) z - Altura sobre el nivel del mar (m)
Signos griegos A t
626
- Incremento 0 diferencia - Transmisividad 6 coef. de transmisi6n
P
- Densidad (kglm3)
Pg Pc
Coef de reflexi6n del suelo (tanto por uno) - Coef. de reflexi6n del cristaI (tanto por uno)
Humedad rclativa del aire (%) - Coef. de absorci6n - Difusividad ternlica del material (a"=A/p Cp) (m2/s) - Conductividad del material (W/mOC 6 kcallh m °C)
a a" A
-
-
o
Coef. de emisi6n 6 emisividad - Cte. de Stefan-Bolzman 5,6697 10-8 W/n12 K4 6 4,965 10-8 kcallh m2 K4
U)
-
€