Iqa222 2021 S2 Ap2 Gi

UNIVERSIDAD TÉCNICA FEDERICO SANTA MARÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA Y AMBIENTAL

TRANSFERENCIA DE CALOR 2do SEMESTRE 2021
Profesor Paula Guerra

Ayudante Maximiliano Ortega
Actividad N°2 Ricardo Torres

Fernanda Troncoso
Integrantes Guillermo Córdova

Conocer, plantear y resolver los problemas asociados a la Simón Díaz
transferencia de calor en estado no estacionario por medio Nicolás Opazo
Conducción en
de herramientas de métodos numéricos, debido a que esto
importante para la modelación de procesosFecha industriales 01-11-2021sólidos
como la de los alimentos.
RESUMEN EJECUTIVO
En este informe se estudió la transferencia de calor en estado no estacionario, se identificó
cuáles son los números adimensionales y las variables asociadas, y en que procesos de la
industria se puede encontrar, de manera que se realizó la resolución de un problema
mediante la simulación con el software ANSYS Fluent y empleando métodos numéricos
para evaluar el efecto de distintas variables en el mecanismo de conducción.
Se estudió la transferencia de calor en estado no estacionario de un cilindro sólido que

posee una temperatura inicial de 278 K, el cual es sometido a un proceso de autoclave
donde circula vapor por el exterior de la lata hasta alcanzar una temperatura constante al
interior de este.
De esta forma se utilizaron materiales como el níquel y el aluminio para analizar la

diferencia debido a conductividades térmicas distintas, que se realizó mediante una
Universidad Técnica Federico Santa María
Departamento de Ingeniería Química y Ambiental
Transferencia de Calor 2021-2. Campus San Joaquín
modelación del comportamiento de la temperatura al interior del sólido a través del tiempo
a partir del software, y el método de diferencias finitas para un cilindro con dimensiones
definidas.
También se analizó como afecta la temperatura de la pared en la temperatura del interior

del sólido a través del tiempo para un mismo material, mediante la simulación en el
software.
Se concluye que tanto la conductividad térmica como la temperatura de la pared son

variables fundamentales que pueden modificar el tiempo de transferencia de calor. Además,
que, por su comportamiento gráfico, ambos métodos utilizados son válidos para realizar
estas aproximaciones, ya que se obtuvieron resultados similares para los tiempos
pertinentes. Las variaciones se pueden atribuir a la diferencia en las propiedades del
material, por el método utilizado, y la elección de los tramos de tiempo y dimensiones para
el cálculo numérico.
ÍNDICE
1. OBJETIVOS................................................................................................................3
2. PREGUNTAS..............................................................................................................3
2.1 Parte I: Teoría...........................................................................................................3
2.2 Parte II: Simulación ANSYS Fluent........................................................................7
4. REFERENCIAS.........................................................................................................14
5. ANEXOS...................................................................................................................15
Actividad 2: Conducción en sólidos Página | 1

ÍNDICE DE FIGURAS:
Figura 1: Gráfico de temperatura vs tiempo en el centro del cilindro (ANSYS Fluent)........7
Figura 2: Gráfico de temperatura vs tiempo en el centro del cilindro para el oro a diferentes
temperaturas de pared.............................................................................................................8
Figura 3: Esquema de elaboración propia del cilindro...........................................................9
Figura 4: Gráfico de temperatura vs tiempo del centro del cilindro mediante el método
diferencias finitas..................................................................................................................12
ÍNDICE DE TABLAS
Tabla 1: Datos de conductividad térmica (k) de distintos materiales.....................................4

Tabla 2: Datos de las Temperaturas de pared de la lata de alimento....................................15
Tabla 3: Datos de conductividad térmica, densidad y capacidad calorífica del aluminio....15
Tabla 4: Datos de conductividad térmica, densidad y capacidad calorífica del níquel........16
1. OBJETIVOS
 Analizar el mecanismo de transferencia de calor en estado no estacionario mediante

el uso de métodos gráficos y métodos de diferencias finitas.
 Describir conceptos asociados a la transferencia de calor en estado no estacionario.
 Simular la transferencia de calor en estado no estacionario mediante el software
ANSYS Fluent para evaluar el efecto de diferentes variables en el mecanismo.
2. PREGUNTAS
2.1 Parte I: Teoría

2.1.1 Defina y explique con sus palabras los números adimensionales de Biot y
Fourier.
Número de Biot: Es un número adimensional que representa la relación que existe entre la
resistencia interior a la transferencia de calor convectiva de un sólido y la resistencia
exterior a la transferencia de calor conductiva alrededor del sólido.
h ∙ Lc
Bi=
k
( 1)
W
Donde h: coeficiente convectivo de transferencia de calor [ ] m2 K
, Lc : Longitud
W
característica [ m ] y k: conductividad térmica [ ]
m °C
.
El número de Biot es utilizado para elegir que método de resolución se utiliza para
problemas de transferencia de calor en estado no estacionario.
o Si Bi < 0,1 se utiliza el método de Lump – Capacitance

o Si Bi > 0,1 se utilizan cartas de temperatura transitoria
El método de Lump – Capacitance considera utilizar la relación (2) a partir del análisis de
un sistema donde puede despreciarse la conducción en el sólido al presentar una alta
resistencia a esta.
−∂ T h ∙ A s ∙ ( T 0−T ∞ )
=
∂t ρ∙ V ∙ C p
(2 )
Donde T 0: Temperatura inicial del sólido [ K ] , T ∞: Temperatura del medio [ K ] , A s: Área

W kg
superficial [ m 2 ], h: coeficiente convectivo de transferencia de calor
[ ]
2
m K
, ρ: densidad [ 3 ]
m
J
, c p: calor específico a presión constante [ ] y V : Volumen [m3 ].
Kg ∙ K
Número de Fourier: Es un número adimensional que representa la relación que existe
entre la velocidad de conducción de calor de un sólido y la velocidad de almacenamiento o
“absorción” de la energía de un sólido.

α ∙t k∙t
Fo= =
Lc ρ ∙C p ∙ L2c
2
(3 )
k m2
Donde α: difusividad térmica (α = [ ]
ρ ∙C p s ¿,
t: tiempo [ s ] , Lc : Longitud característica [ m ]
W kg
k: conductividad térmica [ ]
m °C
, ρ: densidad
m3 [ ] y c p: calor especifico a presión
J
constante [ Kg ∙ K ]
.
2.1.2 Defina que es la conductividad térmica, en que unidades se mide y señale 3

materiales con alta y 3 con baja conductividad.
La conductividad térmica (k) se puede definir como la capacidad que tiene cierto material
para transmitir el calor y corresponde a una propiedad física. La transferencia de calor será
mayor si el material tiene una alta conductividad térmica.
Algunos ejemplos de materiales con alta y baja conductividad térmica son los siguientes:
Tabla 1: Datos de conductividad térmica (k) de distintos materiales. (Young, 1992)
Alta conductividad térmica Baja conductividad térmica

Material W Material W
k
m∙ K [ ] k
m∙ K[ ]
Diamante 1000 Poliestireno 0,033
Plata 406 Poliuretano 0,02
Cobre 385 Corcho 0,04
2.1.3 Suponga que tiene una lámina de largo característico L y una esfera de diámetro
L, ambas del mismo material y a la misma temperatura inicial “cálida”. Estas se
sumergen en una gran cantidad de agua fría ¿Cómo se relaciona la superficie de
contacto con el tiempo en el que se alcanza el estado estacionario en el centro del
cuerpo? Explique con sus palabras para cada uno de los casos.
Tomando el primer caso como uno en el que se sumerge una lámina en agua, el cual posee
un largo característico L, y considerando una esfera con un diámetro L, también sumergida
en una gran cantidad de agua. La longitud característica de cada cuerpo va a depender de su
geometría ya que se define como la razón entre el volumen del sólido y el área del sólido,
por lo tanto, para distintos objetos la longitud característica se define a partir de la ecuación
(4) (Çengel, 2011).

V sólido
Lc =
A sólido
( 4)
∆x
Lc ( Placa Plana ) =
2
r
Lc ( Esféra )=
3
Donde Lc: Longitud Característica, ∆ x : Espesor de la placa plana, r: Radio de la esfera
En este caso se tiene una lámina (aplica para placa plana) con una longitud característica ya
dada, sin embargo, para la esfera solo tenemos su diámetro, por lo tanto, su longitud
característica viene dada por:
Lc ( Lámina ) =L
L
r
L ( Esféra )= =
( 2) L
=
c
3 3 6
Por lo tanto, en este caso para el ejercicio se tienen figuras con distinta longitud
característica.
Las constantes asociadas a la transferencia de calor en estado no estacionario como es el

número de Biot y el de Fourier, están directamente relacionadas con la longitud
característica de tal forma que:
h ∙ Lc Fo= α ∙ t
Bi= y
k L2c
En este caso, lo único que se mantiene para ambos cuerpos de estudio es la conductividad
térmica (k) y la difusividad térmica (α) al ser del mismo material, pero la longitud
característica es distinta, por lo tanto, las variables de Fourier y Biot también serán
distintas. Además, que el coeficiente convectivo (h) ya es distinto por esto mismo.
Finalmente se puede llegar a una relación de tal forma que:
hL ∙ L h ∙L
Bi ( Lámina )= y Bi( esféra)= e
k 6k
α ∙t 36 ∙ α ∙ t
Fo( Lámina)= 2 y Fo(esféra)=
L L2

Por lo tanto, el cálculo del tiempo que alcanza el estado estacionario en el centro del cuerpo
dependerá de los números adimensionales de Biot y Fourier, para todos los métodos de
resolución, ya sea Lump – Capacitance o las cartas de temperatura transitoria. Y para el
caso de este ejercicio cada número adimensional va a variar dependiendo del valor de su
longitud característica, lo cual está directamente relacionado en la distribución del calor de
este mismo.
2.1.4 Explique con sus palabras que es el estado no estacionario y mencione tres
procesos donde se observe este mecanismo de transferencia de calor.
El estado no estacionario dentro de la transferencia de calor se puede definir como un

estado en el cual las variables involucradas varían a lo largo del transcurso del tiempo y del
lugar en donde ocurre la transferencia. Es decir, varían a través del tiempo y la posición
entre un estado inicial y final.
Algunos ejemplos de procesos de transferencia de calor no estacionaria son:
o Enfriamiento o congelamiento de alimentos: Algunos procesos de las industrias

alimentarias utilizan técnicas de conservación para productos necesarios en los
procesos, siendo el enfriamiento a temperaturas muy bajas una de las más utilizadas,
debido a su efectividad al evitar el crecimiento de los microorganismos. En este
caso el factor del tiempo de refrigeración constituye un elemento de alta
importancia práctica, donde se utilizan métodos gráficos que relacionan las
contantes asociadas a la transferencia de calor en estado no estacionario (Nieto,
2015).
o Tostado de maní: Durante el proceso de tostado del maní, el factor del tiempo
resulta fundamental en la práctica para el almacenamiento, freído y enfriado del
maní.
o Esterilización de latas: Las latas son fundamentales para la industria de procesos,

ya que son capaces de almacenar alimentos para que luego, estos sean calificados
como alimentos no perecibles. Sin embargo, para que estos envases puedan cumplir
de buena manera su función es necesario que pasen por un proceso de esterilización,
en el cual la temperatura y el tiempo de duración de este son claves para poder
eliminar todo tipo de bacteria que pueda quedar presente en el material, puesto que
se debe conocer el tiempo exacto que se debe sumergir a las latas en agua a una alta
temperatura específica para realizar el procedimiento de manera eficaz. También
generalmente se trabaja con altas presiones para facilitar la obtención de las
condiciones aptas para el procedimiento (Guerra, 2021).

2.2 Parte II: Simulación ANSYS Fluent
2.2.1. Elabore un gráfico de temperatura en función del tiempo para distintos

materiales de la lata de alimento utilizando el software ANSYS Fluent, analizando lo
obtenido. Para la simulación utilice una temperatura de 400 [𝐾] para la pared del
sólido, que corresponde a la zona de la lata que se encuentre en contacto directo con el
vapor, y trabaje con níquel y aluminio como materiales del sólido.
A partir de los datos obtenidos por el software ANSYS Fluente, sobre la temperatura de la
pared del centro del cilindro respecto al tiempo para el níquel y el aluminio, teniendo una
temperatura de la pared del sólido de 400 [K] y considerando una temperatura inicial de
278[K]. De esta forma se obtuvo el gráfico de la figura 1 para un intervalo de tiempo de
t ∈[0 , 1400(s )].
418
398
Temperatura T [K]
378
358 Níquel
338 Aluminio
318
298
278
0 200 400 600 800 1000 1200 1400
Tiempo t [s]
Figura 1: Gráfico de temperatura vs tiempo en el centro del cilindro (ANSYS Fluent).
A partir del gráfico se puede observar que ambos materiales convergen o se estabilizan a
una temperatura de 400[K], sin embargo, la diferencia está en el tiempo que se demora en
converger. El níquel empieza a converger aproximadamente a los 1000[s], por otra parte, el
aluminio converge mucho antes, aproximadamente a los 300[s].

Una razón por la cual uno de los materiales converge antes que otro se puede encontrar en
las propiedades del material, si se analiza el coeficiente conductividad térmica (k) el níquel
W W
tiene un k = 90,7 [ ]
m∙K
, mientras que el aluminio tiene un k = 237
m∙K
. [ ]
El coeficiente de conductividad térmica es una propiedad física propia de cada material que
determina la capacidad de transmisión de calor. En vista de esta propiedad, se puede
afirmar que el aluminio es un mejor conductor de calor que el níquel, esto implica que el
tiempo requerido para que la temperatura en el centro de la lata iguale el valor de la
temperatura de pared es menor, ya que puede transmitir de mejor forma el calor que el
níquel. (Los datos de k para cada material se encuentran en la tabla 3 y 4 del Anexo B y C).
2.2.2. Elabore un gráfico de temperatura en función del tiempo haciendo uso de los
datos obtenidos del software ANSYS Fluent, para cada una de las temperaturas de
pared que se encuentran en el anexo utilizando como material del sólido oro. Comente
la diferencia entre lo graficado y por qué se genera.
En virtud de los datos obtenidos por el software ANSYS Fluent, se graficó la temperatura
del interior del cilindro respecto al tiempo para diferentes temperaturas de pared y
considerando una temperatura inicial de 278 [K]. Obteniendo el gráfico de la figura 2.
478
Temperatura T [K]
428
T=328[K
]
T=384[K
378 ]
T=425[K
]
T=477[K
328 ]
T=496K]
278
0 20 40 60 80 100 120 140 160 180 200
Tiempo t [s]
Figura 2: Gráfico de temperatura vs tiempo en el centro del cilindro para el oro a diferentes temperaturas de
pared.
A partir de la gráfica realizada, se puede observar que todas las curvas se comportan de
manera similar partiendo desde una temperatura inicial en el tiempo t = 0 [s] igual a

278[K], a medida que transcurre el tiempo existirá una transferencia de calor desde las
paredes hacia el centro del cilindro. La temperatura del centro del cilindro a medida que
transcurre el tiempo tenderá a acercarse a la temperatura de las paredes de este. Cuando la
temperatura de las paredes es mayor, el tiempo que debe transcurrir para que la temperatura
del centro llegue a un valor fijo es mayor, esto se puede observar en el gráfico, ya que la
curva de T = 496[K] empieza a converger aproximadamente a los 180 [s]. Por otro lado,
cuando la temperatura de las paredes es menor, el tiempo necesario para que el centro
converja a dicha temperatura es menor, puesto que existe una diferencia de temperatura
inferior. En el caso de la curva de T = 328[K], se puede observar mediante el gráfico que
esta empieza a dirigirse hacia una cifra estable aproximadamente a los 120 [s].
Estas diferencias entre el tiempo que le toma a cada curva estabilizarse en un número,
dependen netamente de la diferencia de temperaturas entre la pared y el centro. Sin
embargo, se puede notar que los tiempos de convergencia para cada temperatura son
bastante bajos, esta situación se produce ya que el material utilizado para realizar el
experimento es oro, un elemento que particularmente posee una conductividad térmica alta
W
(314 [ ]). Esto explica que, a pesar de que existen diferencias entre las curvas
m∙ K
explicitadas en el gráfico, estas no son tan amplias en relación con el tiempo que tardan en
coincidir en un valor final, pero lo podrían llegar a ser en el caso de que el material del
sólido fuera otro, específicamente uno que posea una menor conductividad térmica.
2.2.3 Realice el gráfico de temperatura en el centro del cilindro en función del tiempo
para el caso en estudio, utilizando el método de diferencias finitas, para ello haga uso
de una discretización adecuada para un sólido cilíndrico. El radio y largo del sólido es
de 0,1 [𝑚] y 1,0 [𝑚] respectivamente, mientras que la temperatura de pared es de
400[𝐾]. Por último, utilice un rango de tiempo de 200 [𝑠] y como material para el
sólido níquel y aluminio. Luego compare y analice con los gráficos obtenidos en la
pregunta 1 con el software ANSYS Fluent.
Se tiene un cilindro con las dimensiones de la figura 3 (Diagrama de elaboración propia

creado mediante la herramienta AutoCAD).

Figura 3: Esquema de elaboración propia del cilindro a analizar

En este caso, se deben definir los siguientes supuestos para la resolución del problema, los
cuales son:
 Conducción en dirección r
 No hay generación de calor
 Capacidad Calorífica constante
 Despreciar efectos de la viscosidad
 Fluido incompresible (Densidad constante)
Para ello se definen las condiciones de contorno y de borde del sólido, que son:
C.I.: t ≤ 0 , T =T 0 ∀ r
C.B. 1: x=0 , T =T ∞ ∀ t >0
C.B. 2: x=L , T=T ∞ ∀ t> 0
Donde x se define en coordenadas cartesianas desde una pared del cilindro hacia la otra
pared y r se define en coordenadas cilíndricas desde el centro del cilindro hasta su radio
máximo (ver figura 3). Además L=0,2[m], T 0=278 [ K ] y T ∞=400[ K ] como se define en
nuestro enunciado.
Entonces, para la resolución del problema se utiliza la ecuación de conservación de energía,

en coordenadas cilíndricas de Navier Stokes (5) (Bird, 2006).

2 2
1 ∂2 T ∂2 T ∂ vz 2
k
1 ∂
( ( )r
∂T
+ 2
r ∂r ∂r r ∂θ ∂z 2
+ 2 + 2 μ[
∂ vr
∂r
¿ ¿ 2+
1 ∂ vθ vr
) ( ) (
r ∂θ r
+ +
∂z ) ( )
]+ ¿ μ r
∂r r
+
r ∂θ [(
∂ ∂ vθ 1 ∂ v r
( )
+
1 ∂ vz ∂ v
+
r ∂θ ∂z )(
∂T ∂ T v ∂T ∂T
¿ ρC (
∂z )
θ
p +v r + +v z
∂t ∂r r ∂ θ
(5)
Como nuestro cilindro sólo presenta conducción solo en dirección r y se desprecian los
efectos de la viscosidad para efectos de cálculo, la ecuación de energía queda de la
siguiente forma:
∂ ∂T ∂T
k
∂r ∂r( ) =ρ C p
∂t
( 6)
En el caso del cálculo de forma numérica con diferencias finitas, es necesario separar la
derivada en diferentes nodos de tal forma que se calcule el efecto la temperatura en cada
uno de esos nodos, por lo tanto, las variables que son necesarias para poder dividirlas en
nodos son la variación de radio ∆ r y la variación de tiempo ∆ t , por lo que cada derivada
queda de la forma:
t +1 t
∂T T m −T m
=
∂t ∆t
(7 )
t t
∂ ∂T 1 ∂T ∂T
( )
∂r ∂r
=
∆r ∂r ( | m+
1
2
−
∂r | )
m−
1
2
t t t t
∂ ∂T 1 T m +1−T m T m −T m −1
( )
∂r ∂r
=
∆r (
∆r
−
∆r )
(8)
Por lo tanto, reemplazando (7) y (8) en (6) generamos la siguiente ecuación diferencial:
T tm+1−T tm t t t t
1 T m +1−T m T m−T m−1
ρ Cp
∆t
=k
∆r ( ∆r
−
∆r )
k
Como la difusividad térmica se define como α = , entonces se despeja la temperatura
ρ ∙C p
en un tiempo t + 1 para un nodo m

t t t t
α T m+1−T m T m −T m−1
T t+1 t
m =T m +
∆r ( ∆r
−
∆r
∆t )
(9)
Finalmente, con la expresión (9) se puede analizar la temperatura que varía para cada nodo
m y en cada tiempo t de manera numérica para cada material de aluminio y níquel, lo que
se define como método de diferencias finitas.
Como para cada material se tienen distintas propiedades caloríficas, la difusividad de cada
uno va a tender a variar, y se tiene el siguiente criterio para obtener el número adimensional
de Fourier, de forma que:
∂t ∆t
Fo=α 2 ≅ α 2
∂r ∆r
( 10 )
Para que el método de diferencias finitas sea válido y pueda converger se debe cumplir que
el número de Fourier debe tener un valor menor a 1/4. (Incropera, 1999)
Para el Níquel se trabajó con una discretización de ∆ t=4 [ s ] y ∆ r =0,02[m] así obteniendo
un Fo = 0,23 cumpliendo con la restricción antes mencionada, por otro lado, para el
aluminio se trabajó con un ∆ t=1 [ s ] y ∆ r =0,02[m], obteniendo un Fo = 0,21 de igual
forma cumpliendo con la restricción.
Fo ( Níquel )=0,23
Fo ( Aluminio )=0,21
Generando el gráfico de la figura 4 a partir del cálculo de cada temperatura de cada nodo en
el intervalo de x ∈(0 ; 0,2 [ m ] ) y teniendo en cuenta que ∆ r =∆ x.

418
398
378
Níquel
Temperatura [K]
358
Aluminio
338
318
298
278
0 20 40 60 80 100 120 140 160 180 200
Tiempo [s]
Figura 4: Gráfico de temperatura vs tiempo del centro del cilindro mediante el método diferencias finitas.
Finalmente, se obtuvo el gráfico para el níquel y aluminio en un intervalo de tiempo de

200[s] del enunciado mediante un método numérico. Si se compara con el gráfico de la
figura 1 obtenido mediante ANSYS Fluent, se observa que ambas curvas poseen el mismo
tipo de comportamiento en el tramo de tiempo correspondiente.
Analizando cuantitativamente, se puede decir, que el aluminio en el gráfico de la figura 1 a

los 200[s] converge a una temperatura de aproximadamente 398 [K] al igual que el gráfico
de la figura 4 mediante el método de diferencias finitas.
Por otro lado, para el caso del níquel la temperatura de convergencia en un tiempo igual a
200[s] es de 348[K] utilizando el primer método, por lo que se observa una pequeña
diferencia con el número obtenido utilizando la metodología numérica el cual es 351 [K] en
el centro del cilindro aproximadamente.
Si se analizan estas cifras obtenidas, se puede notar que ambos métodos llegan a valores
similares, sin embargo, las pequeñas diferencias que se producen son atribuibles al método
en sí que se utilizó para obtenerlas, esto debido a que es necesario tener en cuenta algunas
reglas que pueden afectar el resultado final, tales como que el número de Fourier utilizado
para el cálculo debe ser menor a 0,25, y que la fórmula utilizada para calcular las
temperaturas (ecuación 9) requiere definir intervalos de tiempo y largo arbitrarios que se
adecúen de buena manera a las especificaciones del sistema, las cuales son un radio de 0,1
[m], un tiempo de análisis de 200 [s] y una temperatura de la pared de 400 [K].
Otro factor que se debe tomar en cuenta para el método de diferencias finitas es que
también toma en cuenta las propiedades del material utilizado, la capacidad calorífica, el

coeficiente de conductividad térmica y la densidad, para obtener estas características y

llevar a cabo la metodología se realizó una búsqueda bibliográfica y de esta manera se
obtuvo su valor, mientras que en el programa ANSYS Fluent simplemente se elegía el
material con el que se trabajaría, por lo que tiene este sus propiedades incluidas
implícitamente, aquí es donde se pudo haber producido alguna disyuntiva en estos valores
numéricos que si bien, deben ser muy similares entre sí, pueden causar diferencias en los
resultados finales.
Otra diferencia observada a partir de las gráfica, es que los intervalos de tiempo que se
estudiaron son distintos, ya que para elaborar el primer gráfico se utilizó un tiempo mucho
mayor que llega a los 1400 [s], en cambio para el segundo solo se llegó a los 200 [s], aun
así, si se observa de manera detallada el comportamiento del segundo gráfico y se tienen en
cuenta los valores entregados para la temperatura en un tiempo igual a 200 [s], es posible
manifestar que sería correcto extrapolar las curvas y predecir que estas tendrán la misma
conducta que el otro gráfico al llegar a un tiempo de 1400 [s]. Por lo tanto, ambos métodos
serían correctos en predecir valores de temperatura en el centro del cilindro para distintos
tiempos y ambos se pueden utilizar de buena forma para lograr el objetivo propuesto.
3. CONCLUSIONES
Dentro de la transferencia de calor en estado no estacionario de un sólido, se pueden

encontrar variables asociadas al material de este como la conductividad térmica, la
difusividad térmica, etc., además de variables asociadas a la geometría del objeto como la
longitud característica, esto en conjunto relaciona los números adimensionales de Biot y de
Fourier, los cuales son los principales factores que determinan la distribución de calor a lo
largo del tiempo.
A partir de la simulación realizada se obtuvo del gráfico de la figura 1, en vista de este se

puede concluir que la conductividad térmica (k) de un sólido afecta directamente en el
proceso de transferencia de calor, mientras mayor sea este k, la temperatura en el centro del
objeto tardará menos en converger en un valor definitivo y el material en cuestión será un
mejor conductor de calor. La temperatura de la pared del cilindro también resulta ser un
factor determinante en la transferencia de calor en el tiempo, en vista del gráfico 2 se
observaron diferencias en la convergencia de la temperatura del interior cuando se cambia
esta variable.
Por otra parte, fue posible notar que utilizando tanto el método gráfico obtenido por el
software ANSYS Fluent como el método de diferencias finitas se obtuvieron valores
similares para la temperatura en el centro del cilindro sólido a cierto tiempo. Ambos
métodos se ajustan de manera correcta al sistema presentado y se pueden utilizar de buena

manera para realizar cálculos de temperatura. Las leves variaciones que existen entre ellos
se pueden atribuir a diferencias decimales entre los valores investigados y los implícitos en
el programa para las propiedades del material en cuestión como también a los intervalos de
tiempo y largo utilizados en el método numérico, los cuales pueden causar una alteración.
4. REFERENCIAS
 Bird, R. B., Stewart, W. E., Lightfoot, E. N. (2006) “Transport Phenomena” 2ed.

John Wiley & Sons, Inc.
 Çengel, Y. A., & Ghajar, A. J. (2011). Transferencia de calor y masa: Fundamentos y

aplicaciones. 4a edición. México D.F. McGraw Hill.
 ChemGlobe. (2014). Aluminio (Al) [13] — Elemento químico — Tabla periódica.
ChemGlobe - Tabla periódica de los elementos. https://chemglobe.org/tabla-
periodica/elemento/aluminio/
 ChemGlobe. (2014). Níquel (Ni) [28] — Elemento químico — Tabla periódica.
ChemGlobe - Tabla periódica de los elementos. https://chemglobe.org/tabla-
periodica/elemento/niquel/
 Franco García, A. (2010). Determinación del calor específico de un sólido. Física con
ordenador. Recuperado 29 de octubre de 2021, de
http://www.sc.ehu.es/sbweb/fisica/estadistica/otros/calorimetro/calorimetro.htm
 Guerra, P. (2021). Conducción en sólidos Parte 1. Material del Aula, Universidad
Técnica Federico Santa María, Santiago de Chile.
 Incropera, Frank P. (1999). Fundamentos de transferencia de calor, 4a. ed.
PRENTICE HALL, México.
 Nieto, A. (2015, 20 julio). El Frío en la Conservación de Alimentos -. Mundo Hvacr.

https://www.mundohvacr.com.mx/2008/05/el-frio-en-la-conservacion-de-alimentos/
 Young, H. D. (1992). University physics (7.a ed.). Reading, Mass: Addison-Wesley

Pub. Co.
5. ANEXOS
ANEXO A
“Datos de las Temperaturas de pared de la lata de alimento para el cálculo con el
software”
A continuación, se presenta una tabla con los datos correspondientes de las temperaturas de
pared de lata de alimentos propuestas para el cálculo en el inciso 2.2.1 con el software
ANSYS Fluent.

Tabla 2: Datos de las Temperaturas de pared de la lata de alimento.
T1 (K) T2 (K) T3 (K) T4 (K) T5 (K)

328 384 425 477 496
ANEXO B
“Datos de constantes de conductividad, densidad y capacidad calorífica para el
aluminio”
A continuación, se presenta una tabla con los datos correspondientes de las contantes de
conductividad térmica (k), densidad (ρ) y capacidad calorífica a presión constante (C p) del
aluminio con su respectiva fuente.
Tabla 3: Datos de conductividad térmica, densidad y capacidad calorífica del aluminio.
Variable Valor Unidad Fuente
W
Conductividad térmica (k) 205 (Young, 1992)
m° C
kg
Densidad (ρ) (a 293 [K]) 2702 (ChemGlobe, 2014)
m3
J
Capacidad Calorífica (C p) 900 (ChemGlobe, 2014)
kg ∙ K
ANEXO C
“Datos de constantes de conductividad, densidad y capacidad calorífica para el
níquel”
A continuación, se presenta una tabla con los datos correspondientes de las contantes de
conductividad térmica (k), densidad (ρ) y capacidad calorífica a presión constante (C p) del
níquel con su respectiva fuente.
Tabla 4: Datos de conductividad térmica, densidad y capacidad calorífica del níquel.

Variable Valor Unidad Fuente
W
Conductividad térmica (k) 90,7 (ChemGlobe, 2014)
m° C
kg
Densidad (ρ) (a 293 [K]) 8900 (ChemGlobe, 2014)
m3
J
Capacidad Calorífica (C p) 440 (ChemGlobe, 2014)
kg ∙ K

Iqa222 2021 S2 Ap2 Gi

Cargado por

Copyright:

Formatos disponibles

Iqa222 2021 S2 Ap2 Gi

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Iqa222 2021 S2 Ap2 Gi

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD TÉCNICA FEDERICO SANTA MARÍA

DEPARTAMENTO DE INGENIERÍA QUÍMICA Y AMBIENTAL

Profesor Paula Guerra

Actividad N°2 Ricardo Torres

Integrantes Guillermo Córdova

Se estudió la transferencia de calor en estado no estacionario de un cilindro sólido que

De esta forma se utilizaron materiales como el níquel y el aluminio para analizar la

También se analizó como afecta la temperatura de la pared en la temperatura del interior

Se concluye que tanto la conductividad térmica como la temperatura de la pared son

Actividad 2: Conducción en sólidos Página | 1

Tabla 1: Datos de conductividad térmica (k) de distintos materiales.....................................4

 Analizar el mecanismo de transferencia de calor en estado no estacionario mediante

2.1 Parte I: Teoría

Actividad 2: Conducción en sólidos Página | 2

o Si Bi < 0,1 se utiliza el método de Lump – Capacitance

Donde T 0: Temperatura inicial del sólido [ K ] , T ∞: Temperatura del medio [ K ] , A s: Área

Actividad 2: Conducción en sólidos Página | 3

2.1.2 Defina que es la conductividad térmica, en que unidades se mide y señale 3

Tabla 1: Datos de conductividad térmica (k) de distintos materiales. (Young, 1992)

Alta conductividad térmica Baja conductividad térmica

Actividad 2: Conducción en sólidos Página | 4

Donde Lc: Longitud Característica, ∆ x : Espesor de la placa plana, r: Radio de la esfera

Las constantes asociadas a la transferencia de calor en estado no estacionario como es el

Finalmente se puede llegar a una relación de tal forma que:

Actividad 2: Conducción en sólidos Página | 5

El estado no estacionario dentro de la transferencia de calor se puede definir como un

Algunos ejemplos de procesos de transferencia de calor no estacionaria son:

o Enfriamiento o congelamiento de alimentos: Algunos procesos de las industrias

o Esterilización de latas: Las latas son fundamentales para la industria de procesos,

Actividad 2: Conducción en sólidos Página | 6

2.2 Parte II: Simulación ANSYS Fluent

2.2.1. Elabore un gráfico de temperatura en función del tiempo para distintos

Actividad 2: Conducción en sólidos Página | 7

Actividad 2: Conducción en sólidos Página | 8

Se tiene un cilindro con las dimensiones de la figura 3 (Diagrama de elaboración propia

Actividad 2: Conducción en sólidos Página | 9

Figura 3: Esquema de elaboración propia del cilindro a analizar

Entonces, para la resolución del problema se utiliza la ecuación de conservación de energía,

Actividad 2: Conducción en sólidos Página | 10

Actividad 2: Conducción en sólidos Página | 11

Actividad 2: Conducción en sólidos Página | 12

Finalmente, se obtuvo el gráfico para el níquel y aluminio en un intervalo de tiempo de

Analizando cuantitativamente, se puede decir, que el aluminio en el gráfico de la figura 1 a

Actividad 2: Conducción en sólidos Página | 13

coeficiente de conductividad térmica y la densidad, para obtener estas características y

Dentro de la transferencia de calor en estado no estacionario de un sólido, se pueden

A partir de la simulación realizada se obtuvo del gráfico de la figura 1, en vista de este se

Actividad 2: Conducción en sólidos Página | 14

 Bird, R. B., Stewart, W. E., Lightfoot, E. N. (2006) “Transport Phenomena” 2ed.

 Çengel, Y. A., & Ghajar, A. J. (2011). Transferencia de calor y masa: Fundamentos y

 Nieto, A. (2015, 20 julio). El Frío en la Conservación de Alimentos -. Mundo Hvacr.

 Young, H. D. (1992). University physics (7.a ed.). Reading, Mass: Addison-Wesley

Actividad 2: Conducción en sólidos Página | 15

Tabla 2: Datos de las Temperaturas de pared de la lata de alimento.

T1 (K) T2 (K) T3 (K) T4 (K) T5 (K)

Tabla 3: Datos de conductividad térmica, densidad y capacidad calorífica del aluminio.

Variable Valor Unidad Fuente

Tabla 4: Datos de conductividad térmica, densidad y capacidad calorífica del níquel.