UNIVERSIDAD DE CARTAGENA

FACULTAD DE INGENIERIA
PROGRAMA DE INGENIERIA QUIMICA
TRANSFERENCIA DE CALOR
TAREA N°4

Castro R. José a; Espinosa F., María Patricia a; Escorcia A. Henrya.

Mueses Miguel b.
a
Estudiantes de Ingeniería Química, Universidad de Cartagena
b
Ingeniero químico Ph.D, Docente de la Universidad de Cartagena

Un ladrillo de construcció n mostrado en la figura 1a se encuentra en proceso de

enfriamiento después de que el horno donde se calcinó es apagado. El ladrillo está en
contacto con la superficie del horno a 500ºF. Uno de los operarios del horno derramó
un material polimérico sobre dos de los extremos del ladrillo (véase figura 1b). El
material contaminó estas á reas por tanto se puede considerar que se encuentran
aisladas.

Una corriente de aire a 200ºF proveniente de una zona má s fría del horno, pasa a
través del ladrillo. El coeficiente de convecció n es h = 7 BTU/h.ft 2.ºF, la conductividad
térmica del ladrillo es k = 0.58 BTU/h.ft.ºF, la altura y el ancho del ladrillo son W = H =
12 in y el espesor es δ = ½ in.

Formule las ecuaciones que permiten describir el perfil térmico del ladrillo en cada
una de sus secciones y calcule la transferencia de calor total del sistema.
Fig. 1a Esquema del ladrillo
 Fig. 1b Vista superior del ladrillo

DESARROLLO

Para desarrollar este problema es importante tener en cuenta las siguientes

suposiciones:

 La corriente de aire que llega al ladrillo tiene una temperatura representada

por T∞ y un coeficiente de convecció n hc; estos valores se consideran
invariables y por tanto es posible afirmar que la transferencia de calor se da en
condiciones de estado estable. La condició n de estado estable implica que la
entrada de energía al sistema, por unidad de tiempo, debe ser igual a la energía
neta de salida.
 No hay generació n interna de calor
 La conductividad térmica del ladrillo es constante.
 Los efectos de radiació n desde la superficie son insignificantes.
  La transferencia de calor se da de forma bidimensional, razó n por la cual se
debe acudir al método de diferencias finitas para poder obtener una
distribució n de temperatura a lo largo del ladrillo segú n las condiciones dadas.
 La superficie del ladrillo está sometida a un flujo de calor entrante de valor
constante (q’’).

Retomando las figuras 1a y 1b, se puede observar que dicho ladrillo puede ser
analizado en 3 partes diferentes: la base (en contacto directo con la superficie del
horno), la parte del ladrillo que se encuentra en posició n vertical y la parte superior
del ladrillo (que contiene los extremos aislados). Como ya se dijo anteriormente, el
método má s adecuado para llegar a obtener los perfiles de temperatura es mediante el
método de diferencias finitas.

Existe una particularidad en el sistema propuesto por el ejemplo, ésta corresponde a

una transferencia de calor tanto a lo largo del eje x coordenado como a lo largo del eje
y coordenado y a lo largo del eje z. Por esta razó n se propone entonces utilizar el
método de diferencias finitas tridimensional.

Genéricamente, se sabe que el método de diferencias finitas consiste en aproximar las

ecuaciones diferenciales por ecuaciones en diferencias obtenidas usualmente a partir
del truncamiento de series de Taylor. El conjunto de ecuaciones en diferencias
resultante constituye un sistema de ecuaciones lineales, el cual se puede resolver
numéricamente mediante cualquier tipo de método apropiado para ello y de esta
forma obtener una solució n muy aproximada. El método de diferencias finitas es una
aproximació n para encontrar solució n numérica a las ecuaciones que gobiernan el
modelo matemá tico de un sistema continuo.

Para el caso del sistema tridimensional se debe establecer que Δx = Δy = Δz segú n lo

cual el conjunto de puntos discretos (nodos) quedaría representado de la siguiente
forma:
 Fig. 2 Malla de puntos – arreglo espacial de orden tridimensional

Para la malla escogida se debe cumplir que: u (xi , yj , zk) = u (iΔx , jΔy, kΔz). A
continuació n se indica la deducció n para la malla de puntos nodales escogida:

Se requiere de tres coordenadas (x, y, z) para expresar la distribució n de

temperaturas en el ladrillo, para tal fin se ha considerado el espacio como un conjunto
de celdas cú bicas de dimensiones Δx = Δy = Δz de tal forma que sus nodos extremos
coinciden con sus límites.

Un nodo interno de coordenadas (i, j, k) se representaría como:

Fig. 3 Arreglo de celdas cúbicas para nodos internos

Se considera un nodo interno (i, j, k) de una cá mara cú bica sin generació n de calor, en
estado estacionario. Se identifican las celdas como puntos nodales, de la siguiente
forma:

 Nodo 0 = (i, j, k)  Nodo 4 = (i, j – 1, k)

 Nodo 1 = (i + 1, j, k)  Nodo 5 = (i, j, k + 1)
 Nodo 2 = (i – 1, j, k)  Nodo 6 = (i, j, k – 1)
 Nodo 3 = (i, j + 1, k)

Se le aplica un balance de energía al nodo 0:

ΣQ = 0

Q1 →0 +Q2 → 0+ Q3 →0 +Q4 → 0 +Q5 →0 +Q6 → 0=0

k ∆y∆z ( T ∆−Tx )+ k ∆ y ∆ z ( T ∆−Tx )+ k ∆ x ∆ z ( T ∆−Ty )+ k ∆ x ∆ z ( T ∆−Ty )+ k ∆ x ∆ y ( T ∆−Tz )+ k ∆ x ∆ y ( T

1 0 2 0 3 0 4 0 5 0 6

Desarrollando:

k ∆ y ( T 1−T 0 ) +k ∆ y ( T 2−T 0 ) +k ∆ z ( T 3 −T 0 ) + k ∆ z ( T 4 −T 0 ) +k ∆ x ( T 5 −T 0 ) + k ∆ x ( T 6−T 0 ) =0

Dividiendo por k ∆ y ; y teniendo en cuenta que ∆ x=∆ y=∆ z :

( T 1−T 0 ) + ( T 2−T 0 ) + ( T 3−T 0 )+ ( T 4−T 0 ) + ( T 5−T 0 ) + ( T 6−T 0 )=0

T 1+T 2 +T 3+T 4 +T 5 +T 6−6 T 0=0

T 1 +T 2+ T 3 +T 4 +T 5+ T 6
T 0=
6

Ahora se procede a realizar el mismo aná lisis anterior pero para un punto en una
esquina con convección:
 Fig. 4 Arreglo de celdas cúbicas para nodos en esquina externa con convección

Se identifican las celdas como puntos nodales, de la siguiente forma:

 Nodo 0 = (i, j, k)  Nodo 2 = (i, j, k – 1)

 Nodo 1 = (i + 1, j, k)  Nodo 3 = (i, j + 1, k)

Se le aplica un balance de energía al nodo 0:

ΣQ = 0

Q1 →0 +Q2 → 0+ Q3 →0=0

k
∆y
2
T −T 0
∆z 1
∆x (+k
2
∙
2)
∆ y ∆ x T 2−T 0
∆z
+k
∆x
2 (
T −T 0
∆z 3
∆y
+h
∆x
2 ) (
∆ z ( T ∞ −T 0 )+ h
∆y
2 )
∆ z ( T ∞ −T 0 )=0

2 2
∆y ∆x ∆x ∆x ∆x
k ( T 1−T 0 ) +k ( T 2−T 0 ) +k ( T 3−T 0 )+ h ( T ∞ −T 0 ) + h ( T ∞ −T 0 ) =0
2 4 2 2 2

Multiplicando por 4 /k ∆ y ; teniendo en cuenta que ∆ x=∆ y=∆ z :

4h ∆x
2 ( T 1 −T 0 ) + ( T 2−T 0 ) +2 ( T 3−T 0 ) + ( T ∞ −T 0 )=0
k

4 h∆ x 4h ∆x
2 T 1−2T 0 +T 2 −T 0 +2 T 3−2 T 0 + T∞− T 0=0
k k
 4h∆ x 4 h∆ x
2 T 1 +T 2+ 2T 3−5 T 0+ T ∞− T 0 =0
k k

2 T 1 +T 2+ 2T 3− 5+ ( 4h ∆x
k )
T 0+
4h∆ x
k
T ∞=0

Tomando todos los flujos convectivos, la ecuació n queda expresada como:

2 T 1 +T 2+ 2T 3− 5+ ( 6h∆ x
k )
T 0+
6h∆ x
k
T ∞=0

Se procede ahora a realizar la deducció n para los puntos laterales donde se

presentan los fenó menos de conducció n y convecció n al tiempo.

Fig. 5 Arreglo de celdas cúbicas para nodos laterales

Se le aplica un balance de energía al nodo 0:

ΣQ = 0

Q1 →0 +Q2 → 0+ Q3 →0 +Q4 → 0=0

k
∆y
2
∆z 1
∆x(
T −T 0
+k
2
∙ )
∆ y ∆ x T 2−T 0
2 ∆z
+k
2
∙(
∆ x ∆ y T 3−T 0
2 ∆z )+k
∆x
2
T −T 0
∙∆ z 4
∆y
+h(∆x
2 )
∆ z ( T ∞−T 0 ) + h(∆x
2
∆ )
Dividiendo por 4 /k ∆ y ; teniendo en cuenta que ∆ x=∆ y=∆ z

2h∆ z 2h∆ z 2h ∆ x 2h∆ x

2 ( T 1 −T 0 ) + ( T 2−T 0 ) + ( T 3−T 0 ) +2 ( T 4−T 0 ) + ( T ∞−T 0 ) + ( T ∞−T 0 ) + ( T ∞ −T 0 ) + ( T ∞−T 0
k k k k
Desarrollando:

2h∆ z 2h ∆ z 2h∆ z 2h ∆ z 2h ∆ x 2h∆ x

2 T 1−2T 0 +T 2 −T 0 +T 3−T 0 +2 T 4 −2 T 0 + T ∞− T 0+ T∞− T + T ∞− T 0+
k k k k 0 k k

8h∆ z 8h∆ z
2 T 1 +T 2+ T 3 +2 T 4−6 T 0− T ∞− T 0=0
k k

2 T 1 +T 2+ T 3 +2 T 4− 6+ ( 8h∆ z
k )
T0+
8h∆ z
k
T ∞−¿ 0

Se procede ahora a realizar la deducció n para llegar a obtener la ecuació n

correspondiente a los arreglos de celdas cúbicas ubicados por encima de los
arreglos laterales, para de esta forma ir realizando el barrido de distribució n de
temperaturas a lo largo del ladrillo:

Se realiza mediante una figura aná loga a la figura 5.

Se aplica un balance de energía para la celda cú bica de coordenadas (i, j, k); nodo 0:

ΣQ = 0

Q1 →0 +Q2 → 0+ Q3 →0 +Q4 → 0 +Q 5 →0 =0

k
∆y
2
∆z (
T 1−T 0
∆x
+k )
∆ y ∆ x T 2−T 0
2
∙
2 ∆z
+k( 2
∙
2 )
∆ x ∆ y T 3−T 0
∆z
+k
∆x
2
∙∆ z (
T 4 −T 0
∆y
+k )
∆x
2
∙ ∆z
T 5−T 0
∆y (
+4 h
∆
2 ) ( )
Dividiendo por 4 /k ∆ y ; teniendo en cuenta que ∆ x=∆ y=∆ z :

∆x
2 ( T 1 −T 0 ) + ( T 2−T 0 ) + ( T 3−T 0 ) +2 ( T 4−T 0 ) +2 ( T 5−T 0 ) +8 h ( T ∞−T 0 )=0
k

Desarrollando:

∆x ∆x
2 T 1−2T 0 +T 2 −T 0 +T 3−T 0 +2 T 4 −2 T 0 +2 T 5−2T 0 +8 h T ∞−8 h T =0
k k 0

∆x ∆x
2 T 1 +T 2+ T 3 +2 T 4 +2T 5−8 T 0 +8 h T ∞−8 h T =0
k k 0
 (
2 T 1 +T 2+ T 3 +2 T 4 +2T 5−8 1+h
∆x
k )
T +8 h
0
∆x
k ∞
T =0

Se procede a realizar la deducció n correspondiente a la primera fila interna

considerando los efectos de convecció n:

Fig. 6 Arreglo de celdas cúbicas para la segunda fila interna (nodos cara interna)

Se aplica un balance de energía para la celda cú bica de coordenadas (i, j, k); nodo 0:

ΣQ = 0

Q1 →0 +Q2 → 0+ Q3 →0 +Q4 → 0 +Q5 →0 +Q6 → 0=0

k ∆y∆z ( T 1−T 0
∆x )+k ∆ y ∆ z (
T 2−T 0
∆x )
+k ∆ x∆ z
T 3−T 0
∆y (
+k ∆ x ∆ z )
T 4 −T 0
∆y (
+k ∆ x ∆ y )
T 5−T 0
∆z
+k ∆ x∆ y
T6
( ) (
Dividiendo entre k ∆ y ; y conociendo que ∆ x=∆ y=∆ z :

h∆ z
( T 1−T 0 ) + ( T 2−T 0 ) + ( T 3−T 0 )+ ( T 4−T 0 ) + ( T 5−T 0 ) + ( T 6−T 0 ) + k ( T ∞−T 0 ) =0

Desarrollando:

h∆ z h∆ z
T 1−T 0+T 2−T 0 +T 3−T 0+ T 4−T 0 +T 5 −T 0 +T 6−T 0 + T ∞− T 0=0
k k

h∆ z h∆ z
T 1+T 2 +T 3+T 4 +T 5 +T 6−6 T 0 + T∞− T 0 =0
k k
 (
T 1+T 2 +T 3+T 4 +T 5 +T 6− 6+
h∆ z
k )T0+
h∆ z
k
T ∞=0

Ahora, para la segunda fila interna (nodos situados en la cara interna del sistema
tridimensional), bajo el mismo balance de energía al nodo 0 y un aná lisis similar al
anterior:

ΣQ = 0

Q1 →0 +Q2 → 0+ Q3 →0 +Q4 → 0 +Q5 →0 +Q conv→ 0=0

Se llegó a establecer que su ecuació n correspondiente era:

(
T 1+T 2 +T 3+T 4 +T 5 +T 6−2 3+
h ∆z
k )T 0+
2h∆ z
k
T ∞=0