Modelo Enfriamiento de Un Sartén de Aluminio Mediante El Método de Diferencias Finitas
Modelo Enfriamiento de Un Sartén de Aluminio Mediante El Método de Diferencias Finitas
Modelo Enfriamiento de Un Sartén de Aluminio Mediante El Método de Diferencias Finitas
NRC:2186
Pablo Garcés
Tomás Jáuregui
Juan Vásquez
1. Abstract
2. Introducción
El concepto de diferencias finitas parte de una red nodal, la cual nos propone
sacar la temperatura de forma numérica en un punto especifico. Para esto se establecen
los puntos de interés mediante el trazado de un mallado en el objeto en cuestión, lo que
claramente se puede observar en la imagen 4.4 del texto “Introduction to Heat Transfer”
(Incropera & DeWitt, 2009).
1 δT δ 2 T δ 2 T
= +
α δt δ x 2 δ y 2
4. Resolución:
Figura 4: Plancha (de cocina) analizada por el método de diferencias finitas en estado transitorio.
El material del elemento usado es Aluminio. El mismo tiene dimensiones de
6cm×6cm con un grosor de 0.5cm. A esta se le realizo un corte de la siguiente manera:
Figura 6: Distribución nodal 1 considerando el eje de simetría con la línea entrecortada. Transferencia de calor en la parte
superior, inferior y lado derecho
Figura 7: Distribución nodal 2 considerando el eje de simetría con la línea entrecortada. Transferencia de calor en la parte
superior, inferior y lado derecho.
Estado transitorio
Bidimensional y unidimensional respectivamente
Eje de simetría en la mitad del elemento analizado
Por otro lado, se utilizaron como datos para la resolución de las ecuaciones a los
siguientes:
Δy=Δx=0.025[m]
Kg
ρ Al =2700
[ ]
m3
(Cengel & Boles, 2015)
W
k Al =237 [ ]
m℃
(Cengel & Boles, 2015)
J
c p=945[ ] (Cengel & Boles, 2015)
Kg ℃
W
h R 134 =2000
[ ]
m2 ℃
(Greco & Vanoli, 2005)
T ∞=15 ℃
T I , placa=100 ℃
Para los nodos internos se ocupó la ecuación 5.79. Los nodos internos son: 5, 8,
11, 14, 17, 20, 23, 26, 29, 32 y 35. A continuación se presenta la ecuación para los tres
primeros de dichos nodos. Para los resultados finales, se acomodó la ecuación base para
cada nodo restante, y se las inserto en el programa encargado de las iteraciones.
T P8 +1=Fo ( T 11
P
+ T 5P +T 7P +T P9 ) + (1−4 Fo ) T P8
T P9 +1=Fo[(2 ×T P8 )+T 10
P
+T P4 +(2× Bi ×T ∞)]+[1−( 4 × Fo)−(2 × Bi × Fo)]T P9
Para los nodos correspondientes a las esquinas del lado derecho se usó la
ecuación 5.92. Estos nodos son: 37, 39. Las ecuaciones son:
T P37+1=(2 × Fo)[T 36
P
+T P38+(2 × Bi ×T ∞)]+[1−(4 × Fo)−(4 × Bi × Fo)]T 37
P
T P39+1=(2 × Fo)[T 34
P
+T P38+(2 × Bi ×T ∞ )]+[1−(4 × Fo)−( 4 × Bi × Fo)]T 39
P
Debido a que se trabajó con un eje de simetría en el lado izquierdo del elemento
analizado, se desarrollaron las ecuaciones necesarias para los nodos pertenecientes a
estas zonas. Para los nodos 1 y 3 se obtuvo las ecuaciones de la siguiente manera,
considerando la convección solo en la parte superior del volumen de control:
(T 1P+1 −T 1P) P P P P
Δx (T 2 −T 1 ) Δy (T 6 −T 1 ) Δx P
ρ c p AΔx =k ( ) +k ( ) +h ( )(T ∞−T 1 )
Δt 2 Δy 2 Δx 2
(T 3P+1 −T 3P) P P P P P P
Δx (T 1 −T 2 ) Δy (T 5 −T 2 ) Δx (T 3 −T 2 )
ρ c p AΔx =k ( ) +k ( ) +k( )
Δt 2 Δy 2 Δx 2 Δy
Asumiendo que Δx=Δy, encontrando Fo, se obtuvo:
J
237[ ]
k s m℃ m2
∝= = =9.2887 ×10−5 [ ]
ρcp Kg J s
2700 [ ]× 945[ ]
m2
Kg ℃
W
2000 [ ]× 0.025[m]
h Δx m2 ℃
Bi= = =0.1831
k W
237 [ ]
m℃
Mediante Fo:
−5 m2
9.2887 × 10 [ ]× Δt
∝ Δt s
Fo= =
Δx 2 (0.025 [m])2
1−4 Fo ≥0 ⟶ 1.68[s ]≥ Δt
1−3 Fo ≥ 0⟶ 2.24 [s ]≥ Δt
El criterio que define nuestro Δt es el señalado con rojo, por lo que se escogió un
Δt =1.4[ s]. Por otro lado, se analizó la misma pieza por el método unidimensional y se
obtuvo las siguientes ecuaciones:
Para el nodo 1:
Para el nodo 7:
1−2 Fo ≥0 ⟶ 13.45[s ] ≥ Δt
El criterio que define nuestro Δt es el señalado con rojo, por lo que se escogió un
Δt =11.37 [s].
5. Resultados y Conclusiones
100
90
80
Temperatura del Nodo [C]
70
60
50
40
30
20
10
0
1 9 17 25 33 41 49 57 65 73 81 89 97 105113121129137145153161169177185193201
Tiempo [s]
NODO 1 NODO 39
75
70
65
60
55
50
45
40
35
30
25
20
15
10
5
0
1 9 17 25 33 41 49 57 65 73 81 89 97 105113121129137145153161169177185193201
Tiempo [s]
NODO 2 NODO 38
70
60
50
40
30
20
10
0
1 9 17 25 33 41 49 57 65 73 81 89 97 105113121129137145153161169177185193201
Tiempo [s]
NODO 19 NODO 20
En este caso, se comparó dos nodos aledaños, los cuales difieren en las
condiciones a las que están expuestas. El nodo 19 está en la superficie superior y por
debajo está el nodo 10 como un nodo interno. La tasa de perdida de calor en el 19 es
mayor, pero como son aledaños, la incidencia del uno sobre el otro es mucho más
fuerte, por lo que la decaída de temperatura del nodo 20 es ligeramente más lenta. La
tasa de perdida de calor en el 19 es mayor, pero como son aledaños, la incidencia del
uno sobre el otro es mucho más fuerte, por lo que la decaída de temperatura del nodo
20 es ligeramente más lenta.
Tabla 1: Caso unidimensional. Iteración de temperaturas en cada nodo en función del tiempo.
6. Adjuntos
Iteraciones Proyecto (Primera hoja Bidimensional y segunda hoja
unidimensional)
7. Referencias
Greco, A., & Vanoli, G. (2005). Flow-boiling of R22, R134a, R507, R404A and R410A
inside a smooth horizontal tube. International Journal of Refrigeration, 28(6),
872-880. doi:10.1016/j.ijrefrig.2005.01.008
Incropera, F. P., & DeWitt, D. P. (2009). Introduction to heat transfer. Hoboken, NJ:
Wiley.