Ejemplo 2 Cilindro Sólido PDF

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 6

Ejemplo 2.

Perfil de temperaturas con generación de calor

Un cilindro sólido en el cual hay una generación de calor uniforme interna, q


(W/m3) está aislado en sus extremos. La temperatura en la superficie del
cilindro se mantiene contante, a Tw (K). El calor fluye solo en sentido radial.
a) Deduzca el modelo para el perfil de temperaturas en estado estacionario
cuando el sólido tiene una conductividad térmica constante e indique nivel
de descripción
 R R
r

 T T v T T   1   T  1  T  T
2 2 
C p  + vr + + vz  = k  r + 2 + 2  + SR
 t r r  z  
 r r  r  r  z
2

T
En estado estacionario =0
t
 2T  2T
Solo hay conducción en sentido radial = 0, =0
z 2
 2

Al ser sólido no hay flujo de calor convectivo

 T T v T T   1   T  1  T  T
2 2 
C p  + vr + + vz  = k  r + 2 + 2  + SR
 t r r  z   
 r r  r  r  z
2

Respuesta a)  1   T 
Nivel de descripción k r   = −q
de Gradiente Múltiple  r r  r 
Contiene un término    T   qr Condiciones
  r  = − En r = 0, dT
de dispersión  r  r   k
dr
=0
2
d T dT qr
Aplicando la regla de la r + =− En r = R, T = Tw
cadena en la dr 2 dr k
diferencial
d 2T dT qr Condiciones de frontera
r 2 + =− En r = 0, dT = 0
dr dr k dr
En r = R, T = Tw
b) Obtenga el perfil de temperatura considerando los parámetros
q = 15,000 W/m3 generación d calor
k = 52.3 W/m-K conductividad térmica níquel
Tw = 393 K Temperatura de la superficie
R = 0.5 m radio exterior del cilindro

Se va a resolver con la función dsolve de Matlab


%solución de una EDO de 2do orden
syms T(r) q k Tw R
DT = diff(T);
T = dsolve(r*diff(T,2) + DT + q*r/k == 0', T(R)== Tw, DT(0)== 0);
pretty(T)
q = 15000; % generación d calor (W/m3)
k = 52.3; % conductividad térmica niquel (W/m-K)
Tw = 393; % Temperatura de la superficie (K)
R = 0.5; % radio exterior del cilindro (m)
r = (0:0.01:R); %vector con valores de r

Te = eval(T); %vector con los valores evaluados de T en cada valor de r


% gráfica del perfil de temperatura

plot(r,Te)
Solución de b)
c) Comparar perfil de temperatura para cilindros de los siguientes materiales
níquel
plomo (k = 35 W/m-K) Al variar el parámetro k, no es necesario emplear
bronce (k = 120 W/m-K) en el script ciclos de ‘for’-’if’-’end’-’end’ porque las
variaciones son valores definidos
aluminio (k = 210W/m-K)
figure(2)
k = 52.3; % conductividad térmica níquel (W/m-K)
plot(r,Te)
hold on
k = 35; % conductividad térmica plomo (W/m-K)
Te2 = eval(T);
plot(r, Te2)
hold on
k = 120; % conductividad térmica bronce (W/m-K)
Te3 = eval(T);
plot(r, Te3)
hold on
k = 210; % conductividad térmica aluminio (W/m-K)
Te4 = eval(T);
plot(r, Te4)
xlabel('radio (m)')
ylabel('Temperatura (K)')
grid on
legend('niquel', 'plomo', 'bronce', 'aluminio')
Solución c)
Perfil de temperatura en el cilindro sólido, comparando diferentes materiales

También podría gustarte