UNIVERSIDAD AUTÓNOMA DE NUEVO LEÓN

FACULTAD DE CIENCIAS QUÍMICAS

Producto Integrador
de Aprendizaje
Métodos Númericos

Equipo 4
Mayra Gabriela Castro Guzmán 1958446
David Arturo Flores de León 1961226
Ayrton Yahir De León Vázquez 2078637
Dana Paulina Herrera Villa 1949251
Rodolfo Valdez Morales 1971210
 ÍNDICE
Planteamiento del Códigos
01 Problema 04 Scilab

Descripción del Resultados y

02 fenómeno 05 discusión

Conclusiones
03 Métodologia 06

1
01 Planteamiento del problema
Se coloca una aleta cónica sobre una superficie [Link]és de un
tiempo lo suficientemente largo, la distribución de temperatura alcanza
el estado estacionario. El cono está rodeado de aire cuya temperatura
se puede considerar constante.

Fig 1. Aleta de enfriamiento Cónica 2

 Si en el balance de energía para este sistema se considera que la temperatura en
la dirección radial y angular es uniforme. Además, de considerar el estado
estacionario se llega la ecuación:

Ec. 1

Sujeta a las condiciones de frontera:

en x=0, T=Tw y
en x=L, se puede considerar que la punta está aislada, es decir, dT/dx=0

3
Definición de variables y datos:

Tabla 1. Condiciones Iniciales

Temperatura ambiente (T0) 25°C 298.15 K

Temperatura de la barra en la pared (Tw) 250°C 523.15 K

Coeficiente de transferencia de calor 35 W m^-2 K^-1

convectivo (h)

Conductividad térmica del cono (k) 180 W m^-1 K^-1

Longitud del cono (L) 7.5 cm 0.075 m

Diámetro del cono en la pared (D) 4 cm 0.04 m

B 0.2618 rad

4
a) Determine el perfil de temperatura en la dirección x. Es decir, resuelva la ecuación
diferencial para mostrar la función de la temperatura con respecto a x. Para la
solución utilice el método del disparo con Runge-Kutta de cuarto orden
(proponga un tamaño del paso h que considere adecuado para que la solución tenga
una buena aproximación al valor exacto).
b) Compare los resultados si se considera un cambio en el coeficiente de transferencia
de calor convectivo hc= 20, 50 y 70 W m^-1 K^-2
c) Para todos los casos calcule la pérdida total de calor por convección con la ecuación.

Ec. 2

5
01 Descripción del fenómeno

En la aplicación de efectos combinados de transferencia de calor por

conducción y convección, se utiliza una superficie de manera
extendida para aumentar la rapidez de transferencia de calor entre un
fluido continuo.

Tabla 2. Sistema reportado

T0= 100°C

Calor por convección 10 W*m^-1 *K^-1

Longitud 25 cm
 Sistema de estudio
Los sistemas de
convección y conducción
son combinados en el
área.

Al añadir una aleta a la

Perfil óptimo para la
superficie, el calor fluye
disipación de potencia
provocando la aceleración de
con el mínimo volumen
enfriamiento.

Aumenta la rapidez de Fig 3. Elemento diferencial

Contacto térmico transferencia de calor
entre el sólido y un fluido.

Fig 2. Diagrama del Fenómeno

 Obtención de las funciones
● Ecuación inicial Ec. 1

● L Ec. 3

Ec. 4
●

● Regla del producto Ec. 5

● Ec. 6
 Condiciones de frontera
En x = 0; T =y1= Tw = 523.15 K
En x = L; dT/dx=y2 = 0

Cuando: Cuando:

x= O m x= L= 0.075 m

y1= 523.15 K y2= O

y1= 305.393 K
(propuesto)
Con un tamaño de paso de h=-0.001
 Ecuaciones utilizadas para la resolución
por el método de Runge Kutta, 4to Orden
Ec. 7

Ec. 8

Ec. 9

Ec. 10

Ec. 11
03 Metodología
ALGORITMO
1. Designar las variables iniciales de T0, Tw, h, k, L, D y β, así como sus
cantidades correspondientes
2. Proponer los valores iniciales para h (T paso), e indicar X1, y1 & y2, de
acuerdo a las condiciones de frontera del problema.
3. Proponer un valor para y1 (variable independiente).
4. Expresar las funciones de f2 y f1.
5. Evaluar las funciones en los valores iniciales de x1, y1 & y2.
 Metodología
ALGORITMO
6. Determinar el valor de las K1,j; K2,j; K3,j y K4,j; utilizando las ecuaciones de
Runge Kutta de 4to orden, de acuerdo a la función correspondiente, evaluados
con los valores iniciales.
7. Calcular las nuevas variables: xi+1; y1,i+1; y2,i+1
8. Por lo que ahora x1=xi+1; y1=y1,i+1; y2=y2,i+1
9. Realizar con el método de disparo, variando el valor de y1 propuesto al inicio,
hasta obtener la temperatura definida por las condiciones de frontera.
10. Agregar el comando csvWrite para transferir la matriz de resultados
obtenida a formato de excel.
04 Runge Kutta 4to Orden, con
método del disparo

Fig 4. Código Scilab Para RK 4to orden

 4. Expresar las funciones de f2 y f1.
5. Evaluar las funciones en los valores iniciales de x1, y1 & y2.

ALGORITMO Cálculo de la pérdida total de calor por convección

11. Expresar una nueva función para Q. (Es posible ingresar la función
después del paso 5)
12. Evaluar la función de Q con los valores iniciales de x1 y y1
13. Realizar la integral de la ecuación, evaluando los límites
correspondientes.
14. Evaluar función Q para cada cambio en el coeficiente de calor “h”.
Cálculo de la pérdida total de calor por
convección
Ec. 2

Fig 5. Código Scilab Para regla trapezoidal

05 Resultados y Discusiones
Tabla 3. Resultados de RK 4to orden

Fig 6. Perfil de temperatura en

dirección x
 Tabla 5.
Tabla 4. Comparación de Resultados Para
Resultados Para datos obtenidos RK 4to orden
RK 4to orden en Excel y en Scilab
Excel
Scilab
 Cambio de coeficiente de transferencia de
calor convectivo y considerando un cambio en
el coeficiente de transferencia de calor
convectivo hc= 20, 50 y 70 W m^-1 K^-2

Fig 7. Gráfica Perfil de temperatura

en dirección x (Scilab)

Fig 8. Gráfica Perfil de temperatura

en dirección x (Excel)
Para todos los casos calcule la pérdida total de
calor por convección con la ecuación.

Fig 9. Resultados de la pérdida de Calor (Scilab vs Excel)

Entre mayor sea nuestro coeficiente de transferencia de calor
convectivo, mayor será la cantidad de extracción de calor por
nuestra aleta cónica.
Comparación con datos obtenidos en Excel

Fig 10. Resultados

de la pérdida de
Calor (Excel)
06 Conclusiones

Cuanto mayor sea el coeficiente de transferencia de calor

convectivo, mayor será la eficiencia de la aleta de enfriamiento
cónica lo que lo hace depender del material de nuestra aleta de
enfriamiento

La forma cónica permite un mayor área de superficie en

comparación con una aleta de forma recta, lo que aumenta la
transferencia de calor. Además, la forma cónica puede mejorar
el flujo de aire alrededor de la aleta, lo que aumenta la
velocidad del flujo y promueve una mejor disipación del calor.
 Referencias

● Cengel, Y. A., & Ghajar, A. J. (2007). Transferencia de calor y

masa (Vol. 53). McGraw-Hill Interamericana.

● Bird, R. B., Stewart, W. E., & Lightfoot, E. N. (2020).

Fenómenos de transporte. Reverté.
ANEXOS
Anexo 1. Tabla con el procedimiento realizado para la solución, por Runge Kutta 4to orden
Anexo 2. Tabla con los resultados obtenidos de y1i por cada variación de h.
Anexo 3. Resultados de la integral para el cálculo de las pérdidas de calor variando el coeficiente de
calor
Anexo 3. Resultados de la integral para el cálculo de las pérdidas de calor variando el coeficiente de
calor