CONDUCCIÓN

BIDIMENSIONAL
Método de Diferencias Finitas

Castilla Gleiver Edwin

Castro Fontalvo Rodolfo
García Arias Yeni
Marrugo Quintero Alejandra
Montero Guardia Luis

Transferencia de Calor

2018
 Conducción Bidimensional | Transferencia de Calor

INTRODUCCIÓN
Los fenómenos de transferencia de calor tienen un papel importante en muchos problemas
industriales y ambientales (Incropera & DeWitt, 1999). Por ejemplo, consideremos el área
vital de la producción y conversión de energía, así como en la industria química, metalúrgica,
entre otras. Estos problemas incluyen procesos de conducción, convección y radiación que se
relacionan con el diseño de sistemas.
Este tipo de problemas exigen un análisis permanente de los fenómenos que ocurren, los
cuales normalmente requieren tomar en consideración múltiples variables o parámetros según
la situación. En muchas ocasiones, el análisis está sujeto a la experiencia del personal y requiere
del mismo el manejo eficiente de la información, tal que permita hacer los ajustes
correspondientes a las variables del sistema con la pericia adecuada con fines de lograr su
máxima efectividad. Sin embargo, pueden surgir errores en las estimaciones lo cual conlleva a
posibles errores que afecten los resultados esperados, que, si corresponden a procesos
productivos, pueden llevar a pérdidas de recursos y tiempo (Gonzalez, Pardo, & Surga, 2014).
Existen varias maneras de obtener la formulación numérica de un problema de conducción de
calor, como los métodos de las diferencias finitas, de elementos finitos, de elementos frontera
y de balance de energía (o de volumen de control). Cada uno tiene sus propias ventajas y
desventajas y, en la práctica, se usa cada uno (Cengel, 2007).
El objetivo del trabajo, es emplear el método de diferencias finitas en dos placas de distintos
materiales, con el fin de apreciar el comportamiento de la transferencia de calor en conducción
bidimensional y estado estable, hallando la distribución de temperaturas de los nodos,
experimental y teóricamente, para posteriormente realizar una comparación y estudio entre
los resultados obtenidos del estudio.

1
 Conducción Bidimensional | Transferencia de Calor

CONTENIDO
INTRODUCCIÓN ..................................................................................................................................................... 1
PLANTEAMIENTO DEL PROBLEMA ................................................................................................................... 3
Consideraciones y Datos del Sistema................................................................................................................ 3
Condiciones de Frontera ...................................................................................................................................... 3
Propiedades Físicas y Termodinámicas.............................................................................................................. 3
ALGORITMO DE SOLUCIÓN ............................................................................................................................... 4
BALANCES DE CALOR............................................................................................................................................ 5
Evaluación de las Condiciones de Frontera...................................................................................................... 6
Nodo (I): Superficie plana con flujo de calor constante ........................................................................... 7
Nodo (II): Superficie plana con convección ................................................................................................. 8
Nodo (III): Superficie plana con temperatura constante .......................................................................... 9
Nodo (IV): Superficie plana con convección ............................................................................................... 9
Nodo (V): Esquina con flujo de calor constante y convección ............................................................. 10
Nodo (VI): Esquina con flujo de calor constante y convección ............................................................ 11
Nodo (VII): Esquina con convección y temperatura superficial constante ........................................ 12
Nodo (VIII): Esquina con convección y temperatura superficial constante ....................................... 13
DISEÑO EXPERIMENTAL ...................................................................................................................................... 14
Materiales ............................................................................................................................................................... 14
Montaje Experimental .......................................................................................................................................... 14
Metodología Experimental.................................................................................................................................. 15
REFERENCIAS ........................................................................................................................................................... 17

2
 Conducción Bidimensional | Transferencia de Calor

PLANTEAMIENTO DEL PROBLEMA

A través de un montaje experimental y teórico determinar los perfiles térmicos de una placa
cuadrada que presenta transferencia de calor por conducción bidimensional, además presenta
transferencia de calor por convección gracias al contacto molecular ejercido por un fluido
generado por un secador. Utilizando el método de Diferencias Finitas, validaremos las
temperaturas obtenidas experimentalmente y las temperaturas calculadas para líneas de
corriente; así mismo se realizará para diferentes condiciones.

Ilustración 1. Sistema general de estudio

Consideraciones y Datos del Sistema

Para plantear las ecuaciones características de transferencia de calor por conducción

bidimensional, tendremos en cuenta las siguientes consideraciones globales:

 Sin generación de calor

 Estado estable
 Conductividad térmica constante
 Placa uniforme y cuadrada (∆𝑥 = ∆𝑦)
 Espesor constante y despreciable, debido a que ∆𝑧 ≪ ∆𝑥, ∆𝑦. Entonces, ∆𝑧 = 𝑊 = 1
 Áreas de transferencia de calor son ∆𝑦 ∙ 𝑊 y ∆𝑥 ∙ 𝑊 para el eje x y y, respectivamente.

Condiciones de Frontera

1. Temperatura constante: 2. Convección superficial: 3. Flujo de calor constante:

𝑇𝑠 𝜕𝑇 𝑞"𝑠
−𝐾 𝜕𝑥 | = ℎ[𝑇∞ − 𝑇(𝑡,0) ]
𝑥=0

Propiedades Físicas y Termodinámicas

Propiedad Material
k (W/m.K) Aluminio (Universidad de Cadiz, s.f.) Yeso (Especificar CDT, s.f.)
235 0,26
h (J/Kg.K) 29,04296778
Tabla 1. Propiedades físicas y termodinámicas

3
 Conducción Bidimensional | Transferencia de Calor

ALGORITMO DE SOLUCIÓN

4
 Conducción Bidimensional | Transferencia de Calor

BALANCES DE CALOR
Considerando un volumen de control alrededor de un nodo interior 𝑚, 𝑛

Ilustración 2. Transferencia de calor a un nodo interior desde sus cuatros nodos contiguos

Para condiciones bidimensionales, el intercambio de energía estará influido por la conducción

entre 𝑚, 𝑛 y sus cuatro nodos contiguos:

𝑞"(𝑚+1,𝑛→𝑚,𝑛) ∆𝑦𝑊 + 𝑞"(𝑚,𝑛+1→𝑚,𝑛) ∆𝑥𝑊 + 𝑞"(𝑚−1,𝑛→𝑚,𝑛) ∆𝑦𝑊 + 𝑞"(𝑚,𝑛−1→𝑚,𝑛) ∆𝑥𝑊

+ 𝑞̇ 𝑔𝑒𝑛 ∆𝑥∆𝑦𝑊 = 0
Div: ∆𝑉 = ∆𝑥∆𝑦𝑊
𝑞"(𝑚+1,𝑛) 𝑞"(𝑚,𝑛+1) 𝑞"(𝑚−1,𝑛) 𝑞"(𝑚,𝑛−1)
+ + + + 𝑞̇ 𝑔𝑒𝑛 = 0
∆𝑥 ∆𝑦 ∆𝑥 ∆𝑦
Si ∆𝑥 = ∆𝑦

𝑞"(𝑚+1,𝑛) + 𝑞"(𝑚,𝑛+1) + 𝑞"(𝑚−1,𝑛) + 𝑞"(𝑚,𝑛−1) + 𝑞̇ 𝑔𝑒𝑛 ∆𝑥 = 0

Por tanto:
4

∑ 𝑞(𝑖)→(𝑚,𝑛) + 𝑞̇ 𝑔𝑒𝑛 ∆𝑥∆𝑦 = 0

𝑖=1

Donde (𝑖) son los nodos vecinos.

∆𝑇𝑖,𝑗
Aplicando la ley de Fourier: 𝑞𝑖,𝑗 = −𝑘 ∆𝑥

Se obtiene la rapidez a la que se transfiere la energía por conducción del nodo (𝑖 → 𝑚, 𝑛):
Tabla 2.

5
 Conducción Bidimensional | Transferencia de Calor

Nodo Calor transferido

𝑇(𝑚+1,𝑛) − 𝑇(𝑚,𝑛)
𝒎 + 𝟏, 𝒏 → 𝒎, 𝒏 (𝐾∆𝑦)
∆𝑥
𝑇(𝑚,𝑛+1) − 𝑇(𝑚,𝑛)
𝒎, 𝒏 + 𝟏 → 𝒎, 𝒏 (𝐾∆𝑥)
∆𝑦
𝑇(𝑚−1,𝑛) − 𝑇(𝑚,𝑛)
𝒎 − 𝟏, 𝒏 → 𝒎, 𝒏 (𝐾∆𝑦)
∆𝑥
𝑇(𝑚,𝑛−1) − 𝑇(𝑚,𝑛)
𝒎, 𝒏 − 𝟏 → 𝒎, 𝒏 (𝐾∆𝑥)
∆𝑦
Tabla 2. Calor transferido por los nodos vecinos acoplado a la ley de Fourier

𝑇(𝑖) −𝑇(𝑚,𝑛)
Donde el solo termino ∆𝑦 = ∆𝑥 representa el área de transferencia de calor y , la
∆𝑥=∆𝑦
aproximación en diferencias finitas del gradiente de temperatura en la frontera entre los dos
nodos.
Así la ecuación en diferencias finitas para un nodo interior que sea equidistante de sus cuatro
nodos vecinos con generación es:

𝑞̇ 𝑔𝑒𝑛 (∆𝑥)2
𝑇(𝑚+1,𝑛) + 𝑇(𝑚,𝑛+1) + 𝑇(𝑚−1,𝑛) + 𝑇(𝑚,𝑛−1) + − 4𝑇(𝑚,𝑛) = 0
𝐾
Sin generación de calor se obtiene:

𝑇(𝑚+1,𝑛) + 𝑇(𝑚,𝑛+1) + 𝑇(𝑚−1,𝑛) + 𝑇(𝑚,𝑛−1) − 4𝑇(𝑚,𝑛) = 0

Despejando se obtiene la ecuación del perfil térmico para un nodo interior con 4 nodos
contiguos:
𝑇(𝑚+1,𝑛) + 𝑇(𝑚,𝑛+1) + 𝑇(𝑚−1,𝑛) + 𝑇(𝑚,𝑛−1)
𝑇(𝑚,𝑛) =
4
Evaluación de las Condiciones de Frontera

Para el estudio de las condiciones de frontera tomamos cada uno de los puntos y se realiza el
balance correspondiente a este nodo. Para ello, tenemos en cuenta el siguiente diagrama:

Ilustración 3. Nodos establecidos en la placa

6
 Conducción Bidimensional | Transferencia de Calor

Nodo (I): Superficie plana con flujo de calor constante

Ilustración 4. Flujo de calor constante a un nodo en una superficie externa

Balance de calor:
∆𝑦 ∆𝑦
𝑞"(𝑚−1,𝑛) 𝑊 + 𝑞"(𝑚,𝑛−1) ∆𝑥𝑊 + 𝑞"(𝑚+1,𝑛) 𝑊 + 𝑞"𝑠∆𝑥𝑊 = 0
2 2
Donde: 𝑊 = 1
Aplicamos ley de Fourier:
Calor Trans. Ley de Fourier
𝒒"(𝒎−𝟏,𝒏) 𝑇(𝑚−1,𝑛) − 𝑇(𝑚,𝑛)
−𝐾
∆𝑥
𝒒"(𝒎,𝒏−𝟏) 𝑇(𝑚,𝑛−1) − 𝑇(𝑚,𝑛)
−𝐾
∆𝑦
𝒒"(𝒎+𝟏,𝒏) 𝑇(𝑚+1,𝑛) − 𝑇(𝑚,𝑛)
−𝐾
∆𝑥
Tabla 3. Aplicación de Ley de Fourier a la transferencia por nodo

Análogamente a los nodos, reemplazamos los datos de la Tabla 3, multiplicamos por 2 y

dividimos entre – 𝐾; y además tenemos en cuenta que ∆𝑥 = ∆𝑦.
De esta manera, obtenemos la transferencia que ocurre por un flujo de calor constante, la
cual está expresada de la siguiente manera:
2𝑞"𝑠
𝑇(𝑚−1,𝑛) + 2𝑇(𝑚,𝑛−1) + 𝑇(𝑚+1,𝑛) − 4𝑇(𝑚,𝑛) + ∆𝑥 = 0
𝐾
Así, la ecuación en diferencias finitas para un nodo externo con flujo de calor constante es:
𝑇(𝑚−1,𝑛) 𝑇(𝑚,𝑛−1) 𝑇(𝑚+1,𝑛) 𝑞"𝑠
𝑇(𝑚,𝑛) = + + + ∆𝑥
4 2 4 2𝐾

7
 Conducción Bidimensional | Transferencia de Calor

Nodo (II): Superficie plana con convección

Ilustración 5. Convección para un nodo en una superficie externa

Balance de calor:
∆𝑥 ∆𝑥
𝑞"(𝑚,𝑛+1) 𝑊 + 𝑞"(𝑚,𝑛−1) 𝑊 + 𝑞"(𝑚−1,𝑛) ∆𝑦𝑊 + 𝑞"𝑐𝑜𝑛𝑣 ∆𝑦𝑊 = 0
2 2
Donde: 𝑊 = 1; 𝑞"𝑐𝑜𝑛𝑣 = ℎ(−𝑇𝑚,𝑛 + 𝑇∞ )
Aplicamos ley de Fourier:
Calor Trans. Ley de Fourier
𝑇(𝑚,𝑛+1) − 𝑇(𝑚,𝑛)
𝒒"(𝒎,𝒏+𝟏) −𝐾
∆𝑦
𝑇(𝑚,𝑛−1) − 𝑇(𝑚,𝑛)
𝒒"(𝒎,𝒏−𝟏) −𝐾
∆𝑦
𝑇(𝑚−1,𝑛) − 𝑇(𝑚,𝑛)
𝒒"(𝒎−𝟏,𝒏) −𝐾
∆𝑥
Tabla 4. Aplicación de Ley de Fourier a la transferencia por nodo

Análogamente a los nodos, reemplazamos los datos de la Tabla 4, multiplicamos por 2 y

dividimos entre – 𝐾; y además tenemos en cuenta que ∆𝑥 = ∆𝑦.
De esta manera, obtenemos la transferencia que ocurre por un flujo de calor convectivo, la
cual está expresada de la siguiente manera:
2ℎ𝑇𝑚,𝑛 ∆𝑥 2ℎ𝑇∞ ∆𝑥
𝑇(𝑚,𝑛+1) + 𝑇(𝑚,𝑛−1) + 2𝑇(𝑚−1,𝑛) − 4𝑇(𝑚,𝑛) + − =0
𝐾 𝐾
Así, la ecuación en diferencias finitas para un nodo externo convección, es:
𝑇(𝑚,𝑛+1) 𝑇(𝑚,𝑛−1) 𝑇(𝑚−1,𝑛) ℎ𝑇∞ ∆𝑥/𝐾
𝑇(𝑚,𝑛) = + + −
2ℎ∆𝑥 2ℎ∆𝑥 ℎ∆𝑥 ℎ∆𝑥
(4 − 𝐾 ) (4 − 𝐾 ) (2 − 𝐾 ) (2 − 𝐾 )

8
 Conducción Bidimensional | Transferencia de Calor

Nodo (III): Superficie plana con temperatura constante

Ilustración 6. Temperatura superficial constante para un nodo en una superficie externa

Teniendo en cuenta que es temperatura superficial es constante y conocida, por lo tanto,

para cualquier nodo:

𝑇(𝑚,𝑛) = 𝑇𝑠
Nodo (IV): Superficie plana con convección

Ilustración 7. Convección para un nodo en una superficie externa

Balance de calor:
∆𝑥 ∆𝑥
𝑞"(𝑚,𝑛+1) 𝑊 + 𝑞"(𝑚,𝑛−1) 𝑊 + 𝑞"(𝑚+1,𝑛) ∆𝑦𝑊 + 𝑞"𝑐𝑜𝑛𝑣 ∆𝑦𝑊 = 0
2 2
Donde: 𝑊 = 1; 𝑞"𝑐𝑜𝑛𝑣 = ℎ(−𝑇𝑚,𝑛 + 𝑇∞ )
Aplicamos ley de Fourier:
Calor Trans. Ley de Fourier
𝑇(𝑚,𝑛+1) − 𝑇(𝑚,𝑛)
𝒒"(𝒎,𝒏+𝟏) −𝐾
∆𝑦
𝑇(𝑚,𝑛−1) − 𝑇(𝑚,𝑛)
𝒒"(𝒎,𝒏−𝟏) −𝐾
∆𝑦

9
 Conducción Bidimensional | Transferencia de Calor

𝑇(𝑚+1,𝑛) − 𝑇(𝑚,𝑛)
𝒒"(𝒎+𝟏,𝒏) −𝐾
∆𝑥
Tabla 5. Aplicación de Ley de Fourier a la transferencia por nodo

Análogamente a los nodos, reemplazamos los datos de la Tabla 5, multiplicamos por 2 y

dividimos entre – 𝐾; y además tenemos en cuenta que ∆𝑥 = ∆𝑦.
De esta manera, obtenemos la transferencia que ocurre por un flujo de calor convectivo, la
cual está expresada de la siguiente manera:
2ℎ𝑇𝑚,𝑛 ∆𝑥 2ℎ𝑇∞ ∆𝑥
𝑇(𝑚,𝑛+1) + 𝑇(𝑚,𝑛−1) + 2𝑇(𝑚+1,𝑛) − 4𝑇(𝑚,𝑛) + − =0
𝐾 𝐾
Así, la ecuación en diferencias finitas para un nodo externo convección, es:
𝑇(𝑚,𝑛+1) 𝑇(𝑚,𝑛−1) 𝑇(𝑚+1,𝑛) ℎ𝑇∞ ∆𝑥/𝐾
𝑇(𝑚,𝑛) = + + −
2ℎ∆𝑥 2ℎ∆𝑥 ℎ∆𝑥 ℎ∆𝑥
(4 − 𝐾 ) (4 − 𝐾 ) (2 − 𝐾 ) (2 − 𝐾 )

Nodo (V): Esquina con flujo de calor constante y convección

Ilustración 8. Flujo de calor constante y convección para un nodo en una esquina

Balance de calor:
∆𝑦 ∆𝑥 ∆𝑦 ∆𝑥
𝑞"(𝑚+1,𝑛) 𝑊 + 𝑞"(𝑚,𝑛−1) 𝑊 + 𝑞"𝑐𝑜𝑛𝑣 𝑊 + 𝑞"𝑠 𝑊=0
2 2 2 2
Donde: 𝑊 = 1; 𝑞"𝑐𝑜𝑛𝑣 = ℎ(−𝑇𝑚,𝑛 + 𝑇∞ )
Aplicamos ley de Fourier:
Calor Trans. Ley de Fourier
𝑇(𝑚+1,𝑛) − 𝑇(𝑚,𝑛)
𝒒"(𝒎+𝟏,𝒏) −𝐾
∆𝑥
𝑇(𝑚,𝑛−1) − 𝑇(𝑚,𝑛)
𝒒"(𝒎,𝒏−𝟏) −𝐾
∆𝑦
Tabla 6. Aplicación de Ley de Fourier a la transferencia por nodo

10
 Conducción Bidimensional | Transferencia de Calor

Análogamente a los nodos, reemplazamos los datos de la Tabla 6, multiplicamos por 2 y

dividimos entre – 𝐾; y además tenemos en cuenta que ∆𝑥 = ∆𝑦.
De esta manera, obtenemos la transferencia que ocurre por un flujo de calor constantes y uno
convectivo, la cual está expresada de la siguiente manera:
ℎ𝑇𝑚,𝑛 ℎ𝑇∞ 𝑞"𝑠
𝑇(𝑚+1,𝑛) + 𝑇(𝑚,𝑛−1) − 2𝑇(𝑚,𝑛) + ∆𝑥 − ∆𝑥 − ∆𝑥 = 0
𝐾 𝐾 𝐾
Así, la ecuación en diferencias finitas para un nodo en una esquina externa con convección y
flujo de calor constante, es:
ℎ𝑇∞ ∆𝑥
𝑇(𝑚+1,𝑛) 𝑇(𝑚,𝑛−1) 𝐾 𝑞"𝑠∆𝑥/𝐾
𝑇(𝑚,𝑛) = + − −
ℎ ℎ ℎ ℎ
(2 − 𝐾 ∆𝑥) (2 − 𝐾 ∆𝑥) (2 − 𝐾 ∆𝑥) (2 − 𝐾 ∆𝑥)

Nodo (VI): Esquina con flujo de calor constante y convección

Ilustración 9. Flujo de calor constante y convección para un nodo en una esquina

Balance de calor:
∆𝑦 ∆𝑥 ∆𝑦 ∆𝑥
𝑞"(𝑚−1,𝑛) 𝑊 + 𝑞"(𝑚,𝑛−1) 𝑊 + 𝑞"𝑐𝑜𝑛𝑣 𝑊 + 𝑞"𝑠 𝑊=0
2 2 2 2
Donde: 𝑊 = 1; 𝑞"𝑐𝑜𝑛𝑣 = ℎ(−𝑇𝑚,𝑛 + 𝑇∞ )
Aplicamos ley de Fourier:
Calor Trans. Ley de Fourier
𝑇(𝑚−1,𝑛) − 𝑇(𝑚,𝑛)
𝒒"(𝒎−𝟏,𝒏) −𝐾
∆𝑥
𝑇(𝑚,𝑛−1) − 𝑇(𝑚,𝑛)
𝒒"(𝒎,𝒏−𝟏) −𝐾
∆𝑦
Tabla 7. Aplicación de Ley de Fourier a la transferencia por nodo

Análogamente a los nodos, reemplazamos los datos de la Tabla 7, multiplicamos por 2 y

dividimos entre – 𝐾; y además tenemos en cuenta que ∆𝑥 = ∆𝑦.

11
 Conducción Bidimensional | Transferencia de Calor

De esta manera, obtenemos la transferencia que ocurre por un flujo de calor constantes y uno
convectivo, la cual está expresada de la siguiente manera:
ℎ𝑇𝑚,𝑛 ℎ𝑇∞ 𝑞"𝑠
𝑇(𝑚−1,𝑛) + 𝑇(𝑚,𝑛−1) − 2𝑇(𝑚,𝑛) + ∆𝑥 − ∆𝑥 − ∆𝑥 = 0
𝐾 𝐾 𝐾
Así, la ecuación en diferencias finitas para un nodo en una esquina externa con convección y
flujo de calor constante, es:
ℎ𝑇∞ ∆𝑥
𝑇(𝑚−1,𝑛) 𝑇(𝑚,𝑛−1) 𝐾 𝑞"𝑠∆𝑥/𝐾
𝑇(𝑚,𝑛) = + − −
ℎ ℎ ℎ ℎ
(2 − 𝐾 ∆𝑥) (2 − 𝐾 ∆𝑥) (2 − 𝐾 ∆𝑥) (2 − 𝐾 ∆𝑥)

Nodo (VII): Esquina con convección y temperatura superficial constante

Ilustración 10. Convección y temperatura superficial constante para un nodo en una esquina

Balance de calor:
∆𝑥 ∆𝑦 ∆𝑦
𝑞"(𝑚−1,𝑛) 𝑊 + 𝑞"(𝑚,𝑛+1) 𝑊 + 𝑞"𝑐𝑜𝑛𝑣 𝑊=0
2 2 2
Donde: 𝑊 = 1; 𝑞"𝑐𝑜𝑛𝑣 = ℎ(−𝑇𝑚,𝑛 + 𝑇∞ )
Aplicamos ley de Fourier:
Calor Trans. Ley de Fourier
𝑇(𝑚−1,𝑛) − 𝑇(𝑚,𝑛)
𝒒"(𝒎−𝟏,𝒏) −𝐾
∆𝑥
𝑇(𝑚,𝑛+1) − 𝑇(𝑚,𝑛)
𝒒"(𝒎,𝒏+𝟏) −𝐾
∆𝑦
Tabla 8. Aplicación de Ley de Fourier a la transferencia por nodo

Análogamente a los nodos, reemplazamos los datos de la Tabla 8, multiplicamos por 2 y

dividimos entre – 𝐾; y además tenemos en cuenta que ∆𝑥 = ∆𝑦 y 𝑇(𝑚−1,𝑛) = 𝑇𝑠 .
De esta manera, obtenemos la transferencia que ocurre por un flujo convectivo y temperatura
superficial constante, la cual está expresada de la siguiente manera:

12
 Conducción Bidimensional | Transferencia de Calor

ℎ𝑇𝑚,𝑛 ℎ𝑇∞
𝑇(𝑚−1,𝑛) + 𝑇(𝑚,𝑛+1) − 2𝑇(𝑚,𝑛) + ∆𝑥 − ∆𝑥 = 0
𝐾 𝐾
Así, la ecuación en diferencias finitas para un nodo en una esquina externa con estas
condiciones es:
𝑇𝑠 𝑇(𝑚,𝑛+1) ℎ𝑇∞ ∆𝑥/𝐾
𝑇(𝑚,𝑛) = + −
ℎ ℎ ℎ
(2 − 𝐾 ∆𝑥) (2 − 𝐾 ∆𝑥) (2 − 𝐾 ∆𝑥)

Nodo (VIII): Esquina con convección y temperatura superficial constante

Ilustración 11. Convección y temperatura superficial constante para un nodo en una esquina

Balance de calor:
∆𝑥 ∆𝑦 ∆𝑦
𝑞"(𝑚+1,𝑛) 𝑊 + 𝑞"(𝑚,𝑛+1) 𝑊 + 𝑞"𝑐𝑜𝑛𝑣 𝑊=0
2 2 2
Donde: 𝑊 = 1; 𝑞"𝑐𝑜𝑛𝑣 = ℎ(−𝑇𝑚,𝑛 + 𝑇∞ )
Aplicamos ley de Fourier:
Calor Trans. Ley de Fourier
𝑇(𝑚+1,𝑛) − 𝑇(𝑚,𝑛)
𝒒"(𝒎+𝟏,𝒏) −𝐾
∆𝑥
𝑇(𝑚,𝑛+1) − 𝑇(𝑚,𝑛)
𝒒"(𝒎,𝒏+𝟏) −𝐾
∆𝑦
Tabla 9. Aplicación de Ley de Fourier a la transferencia por nodo

Análogamente a los nodos, reemplazamos los datos de la Tabla 9, multiplicamos por 2 y

dividimos entre – 𝐾; y además tenemos en cuenta que ∆𝑥 = ∆𝑦 y 𝑇(𝑚+1,𝑛) = 𝑇𝑠 .
De esta manera, obtenemos la transferencia que ocurre por un flujo convectivo y temperatura
superficial constante, la cual está expresada de la siguiente manera:
ℎ𝑇𝑚,𝑛 ℎ𝑇∞
𝑇(𝑚+1,𝑛) + 𝑇(𝑚,𝑛+1) − 2𝑇(𝑚,𝑛) + ∆𝑥 − ∆𝑥 = 0
𝐾 𝐾

13
 Conducción Bidimensional | Transferencia de Calor

Así, la ecuación en diferencias finitas para un nodo en una esquina externa con estas
condiciones es:
𝑇𝑠 𝑇(𝑚,𝑛+1) ℎ𝑇∞ ∆𝑥/𝐾
𝑇(𝑚,𝑛) = + −
ℎ ℎ ℎ
(2 − 𝐾 ∆𝑥) (2 − 𝐾 ∆𝑥) (2 − 𝐾 ∆𝑥)

DISEÑO EXPERIMENTAL
Materiales

Para la realización de esta práctica, se requirieron los siguientes materiales:

 Placa plana de aluminio galvanizado de 18 cm x 18 cm
 Placa plana de drywall de 18 cm x 18 cm
 Lamina de poliestireno
 Termopares
 Plancha digital
 Soportes
 Plancha de ropa.
 Secador de cabello

Montaje Experimental

14
 Conducción Bidimensional | Transferencia de Calor

Metodología Experimental

Para llevar a cabo el montaje de esta práctica y posteriormente el desarrollo de esta, se hizo
necesario seguir una serie de pasos y condiciones para terminarla de manera exitosa. Los
pasos fueron los siguientes:
PASOS:
1. Para comenzar, tomamos dos placas planas de 18 cm x 18 cm cada una de un material
diferente, drywall y aluminio.
2. Para cada placa se tomaron láminas de poliestireno (aislante), las cuales fueron pegadas
a las placas de ambos lados; y de esta manera no se presentarán alteraciones en las
temperaturas por la transferencia de calor.
3. En cada lámina de aislante fue trazada una cuadrícula para simular las particiones de la
superficie.
4. Cada lámina fue calentada con la plancha de ropa.
5. Posteriormente, se introdujo unos termopares en cada punto de las placas para así medir
las temperaturas de dichos puntos.
6. Las condiciones de frontera que se simularon experimentalmente se describieron al
inicio, en el planteamiento del problema.

Ilustración 12. Cuadrícula trazada para el aislante de cada placa, con su respectiva enumeración en los nodos

ANALISIS RESULTADOS
Para obtener el perfil térmico de la placa, se debe resolver el sistema de ecuaciones lineales
que se forman al plantear los balances de calor para cada nodo, para resolver este sistema es
necesario utilizar un método numérico, en este caso, se utilizó el método de referencia circular
con ayuda de la herramienta informática de Microsoft Excel.
Antes de explicar en qué consiste el método de cálculo iterativo, es necesario aclarar en qué
consiste las referencias circulares; Estas, surgen cuando existen ecuaciones que se refieren a

15
 Conducción Bidimensional | Transferencia de Calor

otras y estas dependen unas de otras. Las referencias circulares devuelven un error muy típico
que condiciona el resto de resultados de la hoja de cálculo sobre la que se trabaja. Ahora, esa
referencia circular se debe romper de una u otra forma, y es ahí donde entra la opción de
habilitar cálculo iterativo, que permite calcular tantas veces como se le haya indicado hasta
alcanzar un resultado con la precisión definida.
Por otro lado, se sabe que nuestro proceso es conductivo y convectivo (para ello solo
empleamos la velocidad para calcular el coeficiente de convección); por lo tanto, para
comparar los valores de temperatura experimentales con los valores de temperatura teóricos,
se creó una matriz de 50 x 50, de tal modo que se obtuvieran 2500 particiones en la placa.

16
 Conducción Bidimensional | Transferencia de Calor

REFERENCIAS
Cengel, Y. A. (2007). Transferencia de Calor y Masa (Tercera ed.). México: McGraw - Hill
Interamericana.
Gonzalez, T. C., Pardo, M., & Surga, R. (2014). Simulación Bidimensional de la Ecuación de
Calor. (U. d. Guajira, Ed.) Ciencia e Ingeniería, 1. Obtenido de
[Link]
Incropera, F., & DeWitt, D. (1999). Fundamentos de Transferencia de Calor (Cuarta ed.).
Prentice Hall.

17