Ejemplo 2 Cilindro Sólido PDF
Ejemplo 2 Cilindro Sólido PDF
Ejemplo 2 Cilindro Sólido PDF
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
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
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