> restart:#"m11_p16"

Se trata de analizar el régimen estacionario de una instalación de agua caliente sanitaria que consta de un calentador, donde la corriente de agua recibe 20 kW, y de una tubería de 1 cm de diámetro y 10 m de longitud cuya transmitancia térmica global con el ambiente se supone que es K=200 W∙m 2∙K 1 referido al área del tubo. Se pide:
a) Temperatura de salida del calentador, T1 (en función del gasto másico de agua).
b) Perfil de temperatura a lo largo del conducto.
c) Temperatura de salida del conducto, T2, (en función del gasto másico de agua).
d) Gasto másico que daría máxima temperatura de salida del conducto y valor de ésta.
e) Comentar la viabilidad y conveniencia del punto de trabajo anterior.
Datos:

> read`../therm_eq.m`:read`../therm_const.m`:read`../therm_proc.m`:with(therm_proc):

> su:="H2O":dat:=[Q=20e3*W_,R=0.005*m_,L=10*m_,K=200*W_/(m_^2*K_)];

[Q = `+`(`*`(0.20e5, `*`(W_))), R = `+`(`*`(0.5e-2, `*`(m_))), L = `+`(`*`(10, `*`(m_))), K = `+`(`/`(`*`(200, `*`(W_)), `*`(`^`(m_, 2), `*`(K_))))]

Image

> ldat:=get_liq_data(su):dat:=op(dat),ldat,Const,SI1,SI2:

a) Temperatura de salida del calentador, T1 (en función del gasto másico de agua).

Sea T0 la de entrada y T1 la de salida.

> eqBE:=Q=m*c*DT;eqBE1:=T1=T0+Q/(m*c);

Q = `*`(m, `*`(c, `*`(DT)))
T1 = `+`(T0, `/`(`*`(Q), `*`(m, `*`(c))))

b) Perfil de temperatura a lo largo del conducto.

(Ver arriba.)

c) Temperatura de salida del conducto, T2, (en función del gasto másico de agua).

> eqBE2:=dH/dt=Qcond+Qconv+Qwall;eqBE2:=m*c*diff(T(x),x)=k*Pi*R^2*diff(T(x),x,x)-K*2*Pi*R*(T(x)-T0);eqBE2:=m*c*diff(T(x),x)=-K*2*Pi*R*(T(x)-T0);dsol1:=subs(eqBE1,dsolve({eqBE2,T(0)=T1},T(x)));eq2:=subs(T(x)=T2,x=L,dsol1);

`/`(`*`(dH), `*`(dt)) = `+`(Qcond, Qconv, Qwall)
`*`(m, `*`(c, `*`(diff(T(x), x)))) = `+`(`*`(k, `*`(Pi, `*`(`^`(R, 2), `*`(diff(diff(T(x), x), x))))), `-`(`*`(2, `*`(K, `*`(Pi, `*`(R, `*`(`+`(T(x), `-`(T0)))))))))
`*`(m, `*`(c, `*`(diff(T(x), x)))) = `+`(`-`(`*`(2, `*`(K, `*`(Pi, `*`(R, `*`(`+`(T(x), `-`(T0)))))))))
T(x) = `+`(T0, `/`(`*`(exp(`+`(`-`(`/`(`*`(2, `*`(K, `*`(Pi, `*`(R, `*`(x))))), `*`(m, `*`(c)))))), `*`(Q)), `*`(m, `*`(c))))
T2 = `+`(T0, `/`(`*`(exp(`+`(`-`(`/`(`*`(2, `*`(K, `*`(Pi, `*`(R, `*`(L))))), `*`(m, `*`(c)))))), `*`(Q)), `*`(m, `*`(c))))

d) Gasto másico que daría máxima temperatura de salida del conducto y valor de ésta.

Si m es pequeño se calienta mucho pero pierde mucho también; si m es grande se calienta menos pero pierde menos.

> m_Tm=T2(m)[max];m_Tm:=solve(diff(rhs(eq2),m)=0,m);m_Tm_:=evalf(subs(dat,m_Tm));T2m:=subs(m=m_Tm,rhs(eq2));T2m_:=evalf(subs(dat,T2m));

m_Tm = T2(m)[max]
`+`(`/`(`*`(2, `*`(K, `*`(Pi, `*`(R, `*`(L))))), `*`(c)))
`+`(`/`(`*`(0.1503e-1, `*`(kg_)), `*`(s_)))
`+`(T0, `/`(`*`(`/`(1, 2), `*`(exp(-1), `*`(Q))), `*`(K, `*`(Pi, `*`(R, `*`(L))))))
`+`(`*`(405.1, `*`(K_)))

Si representamos la temperatura tras el calentador, T1, y la del extremo de salida, T2, en función del gasto de agua:

> plot(subs(dat,SI0,{rhs(eqBE1),rhs(eq2),subs(dat,SI0,[[0,T0],[1,T0]])}),m=0..0.1,Temp=0..500,color=black);

Plot_2d

No es realista esa Tmax, que habrá de ser menor que la de ebullición a presión ambiente. En realidsad, para agua caliente sanitaria basta conque T2=45 ºC ( y la T1 tampoco puede ser mayor de unos 100 ºC para evitar la formación de vapor).

e) Comentar la viabilidad y conveniencia del punto de trabajo anterior.

Efectivamente, vemos que no es realista considerar gastos másicos pequeños; un gasto típico de un grifo es ya 0,1 kg/s (6 L/min).

>