UNIVERSIDAD POLITÉCNICA DE PUEBLA

ORGANISMO PÚBLICO DESCENTRALIZADO

DEL ESTADO DE PUEBLA

INGENIERÍA MECATRÓNICA

TRANSFERENCIA DE CALOR

ACTIVIDAD

Conducción de calor en dos dimensiones

Solución analítica

ALUMNO:

ISAAC ALEJANDRO VALDEZ

ISRAEL DÍAZ CORTES

ASESORES:

DR. OBED CORTÉS ABURTO

JUAN C. BONILLA, PUEBLA NOVIEMBRE 2018

1
INTRODUCCION

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. Este proceso probablemente empírico, pierde
vigencia cuando se requiere un proceso analítico y sistemático fundamentado en la rigurosidad
científica. El empleo de una herramienta computacional como un simulador, permitiría que la
solución que se le dé al problema esté más ajustada a las necesidades, facilitando que tanto el
control del proceso como la toma de decisiones sean llevadas a cabo por métodos comprobables y
repetitivos. El objetivo del trabajo, es desarrollar una herramienta informática que permita obtener
la solución de las ecuaciones diferenciales en derivadas parciales de segundo orden de tipo
hiperbólico, que rigen el comportamiento de la transferencia de calor en superficies
bidimensionales, estrategia para el análisis de la superficie, el diagrama de flujo y seguidamente se
muestran los resultados y discusión de la simulación por el método de diferencias finitas, usando las
herramientas que nos ofrece el software MATLAB. Finalmente, se expresan las conclusiones del
estudio.

2
Contenido
DESARROLLO DE LAS OPERACIONES MATEMATICAS.......................................................................... 4
DONDE............................................................................................................................................. 4
A continuación de muestra el cambio para poder realizar las operaciones ............................... 4
A continuacion mostramos ya aplicadas las condiciones de barrera ......................................... 6
OBTENCION DE GRAFICAS POR MEDIO DEL SOFTWARE MATLAB ...................................................... 8
A continuacion los valores designados fueron los siguientes: ........................................................ 8
Se comenzo graficando los valores de n en graficas separadas. ........................................................ 9
Para n = 1 se obtiene la siguiente grafica........................................................................................ 9
Para n= 5 se obtiene la siguiente grafica. ....................................................................................... 9
Para n= 10 se obtiene la siguiente gráfica. ................................................................................... 10
Para n= 50 se obtiene la siguiente gráfica. ................................................................................... 10
Para n=100 se obtiene la siguiente grafica ................................................................................... 11
CONCLUSIONES ................................................................................................................................. 13
REFERENCIAS BIBLIOGRAFICAS ......................................................................................................... 13

3
DESARROLLO DE LAS OPERACIONES MATEMATICAS

𝑥=0 𝑇(𝑥, 𝑦) = 𝑇1 → 𝑇(0, 𝑦) = 𝑇1 → 𝜃(0, 𝑦) = 0

𝑥=𝐿 𝑇(𝑥, 𝑦) = 𝑇1 → 𝑇(𝐿, 𝑦) = 𝑇1 → 𝜃(𝐿, 𝑦) = 0
𝑦=0 𝑇(𝑥, 𝑦) = 𝑇2 → 𝑇(𝑥, 0) = 𝑇1 → 𝜃(𝑥, 0) = 0
𝑦=𝑊 𝑇(𝑥, 𝑦) = 𝑇2 → 𝑇(𝑥, 𝑤) = 𝑇2 → 𝜃(𝑦, 𝑤) = 1

DONDE

𝑇 − 𝑇1
𝜃=
𝑇2 − 𝑇1
A continuación de muestra el cambio para poder realizar las operaciones

𝜕𝜃 𝜕𝑥 𝜕𝜃 𝜕𝑦
= 𝑦 = 𝑥
𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑦

𝜕 2𝜃 𝜕 2𝜃 𝜕 2𝜃 𝜕 2𝑦
= 𝑦 =
𝜕𝑥 2 𝜕𝑥 2 𝜕𝑦 2 𝜕𝑦 2

2 2
1 𝑦𝜕 𝑥 𝜕 𝑦
( 2 = 𝑥 2=0 )
𝑥, 𝑦 𝜕𝑥 𝜕𝑦

4
 2 2
1 𝜕 𝑥 1 𝜕 𝑦
( + =0 )
𝑥 𝜕𝑥2 𝑦 𝜕𝑦2

2 2
1𝑑 𝑥 1 𝑑 𝑥
2
= − 2
= 𝜆2
𝑦 𝑑𝑦 𝑦 𝑑𝑥

2
1 𝑑 𝑦
= 𝜆2
𝑦 𝑑𝑦2

𝑑2 𝑦
− 𝜆2 𝑦 = 0
𝑑𝑦 2
𝜃(𝑥, 𝑦) = (𝐶1 cos 𝜆𝑥 + 𝐶2 Sen 𝜆𝑥) (𝑦 = 𝐶3 𝑒 −𝜆𝑦 + 𝐶4 𝑒 𝜆𝑦 )

𝐶1 = 0 𝐶2 Sen 𝜆𝑥 (𝐶3 𝐶4 ) = 0
𝐶3 = −𝐶4

(𝐶2 Sen 𝜆𝑥)(− 𝐶3 𝑒 −𝜆𝑦 + 𝐶4 𝑒 𝜆𝑦 )

𝐶2 𝐶4 Sen 𝜆𝐿(𝑒 𝜆𝑦 − 𝑒 −𝜆𝑦 ) = 0

∞
𝑛𝜋𝑥 𝑛𝜋𝑦
𝜃 = ∑ 𝐶𝑛 (sinh )(sinh )
𝐿 𝐿
𝑛=1

∫ 𝑔𝑚 𝑥 𝑔𝑛 (𝑥)𝑑𝑥 = 0 𝑚≠𝑛
𝑎

5
 ∞

𝑓(𝑥) = ∑ 𝐴𝑛 𝑔𝑛 (𝑥)
𝑛=1
∞

𝑓(𝑥)𝑔𝑚 (𝑥) = 𝑔𝑛 (𝑥) ∑ 𝐴𝑛 𝑔𝑛 (𝑥)

𝑛=1

𝑏 𝑏 ∞

∫ 𝑓(𝑥) 𝑔𝑚 (𝑥)𝑑𝑥 = ∫ 𝑏𝑛 (𝑥) ∑ 𝐴𝑛 𝑔𝑛 (𝑥) 𝑑𝑥

𝑎 𝑎 𝑛=1

𝑏 ∞ 𝑏

∫ 𝑓(𝑥) 𝑔𝑚 (𝑥)𝑑𝑥 = ∑ 𝐴𝑛 ∫ 𝑔𝑚 (𝑥) 𝑑𝑥

𝑎 𝑛=1 𝑎

𝑏 𝑏

∫ 𝑓(𝑥) 𝑔𝑚 (𝑥)𝑑𝑥 = ∫ 𝑔𝑚 (𝑥)𝑑𝑥

𝑎 𝑎
𝑏
∫𝑎 𝑓(𝑥)𝑔𝑚 (𝑥)𝑑𝑥
𝐴𝑚 𝑏
∫𝑎 𝑔𝑚 (𝑥)𝑑𝑥

A continuacion mostramos ya aplicadas las condiciones de barrera

∞
𝑛𝜋𝑤 𝑛𝜋𝑤
∑ 𝐶𝑛 sin sinh
𝐿 𝐿
𝑛=1

𝐿 𝑛𝜋𝑤
∫0 sin 𝐿 1 − 𝑐𝑜𝑠𝑛𝜋
𝐴𝑚 = 𝐿 𝑛𝜋𝑤 = 𝑛𝜋 1
∫0 sinh 𝐿
2 − 4 𝑠𝑒𝑛2𝜋𝑛

2
𝐴𝑚 = (1 − (−1)𝑛
𝑛𝜋

6
 ∞
2 𝑛𝜋𝑥
1= ∑ (1 − (−1)𝑛 ) sin
𝑛𝜋 𝐿
𝑛=1

2
(1 − (−1)𝑛 )
𝐶𝑛 = 𝑛𝜋
𝑛𝜋𝑤
sinh 𝐿

∞ 𝑛𝜋𝑦
2 1 + (1)𝑛+1 𝑛𝜋𝑥 sinh 𝐿
𝜃(𝑥, 𝑦) = ∑ sin
𝜋 𝑛 𝐿 sinh 𝑛𝜋𝑤
𝑛=1 𝐿

7
OBTENCION DE GRAFICAS POR MEDIO DEL SOFTWARE MATLAB

Para la elaboracion y obtencion de las graficas solicitadas por medio del software matlab se
desarrollo el respectivo codigo en MATLAB el cual nos permite realizar el calculo
respectivamente para cada valor asiganado en las siguientes actividades.

A continuacion los valores designados fueron los siguientes:

1. n=1
2. n=5
3. n=10
4. n=50
5. n=100

Los cuales se graficaron en el software MATLAB a apartir de la ecuacion final obtenida

mediante el analisis matematico que se ha descrito previamente; la ecuacion
correspondiente que sera ingresada en el editor de MATLAB es la siguiente:

∞ 𝑛𝜋𝑦
2 1 + (1)𝑛+1 𝑛𝜋𝑥 sinh 𝐿
𝜃(𝑥, 𝑦) = ∑ sin
𝜋 𝑛 𝐿 sinh 𝑛𝜋𝑤
𝑛=1 𝐿

Una vez que ya esta definida la ecuacion a graficar se procede a ingresar la programacion
adecuada para poder obtener la graficas de los diferentes valores que se le puede dar a “n”,
que en este caso son los mencionados anteriomente.

Para los valores que se han establecido para “n”, se han definido las siguientes
caracteristicas:
Donde “x” que es nuestra “L” tendra un valor de 5 y donde nuestre eje “Y” que sera “W”
obtendra un valos especifico de 10, con un espaciamiento de 0.1 en ambos lados para la
obtencion de mejores resultados.

8
Se comenzo graficando los valores de n en graficas separadas.

Para n = 1 se obtiene la siguiente grafica.

Figura numero 1.- Donde N= 1

Para n= 5 se obtiene la siguiente grafica.

Figura numero 2.- Donde N= 5

9
Para n= 10 se obtiene la siguiente gráfica.

Figura numero 3.- Donde N= 10

Para n= 50 se obtiene la siguiente gráfica.

Figura numero 4.- Donde N= 50

10
Para n=100 se obtiene la siguiente grafica

Figura numero 5.- Donde N= 100

11
De esta manera es como quedarían las gráficas independientes para cada valor de “n”
asignado.
A continuación se procede a mostrar las líneas de flujo de calor en dos dimensiones de
una placa rectangular que es el objetivo de estudio en este tema.

Figura numero 6.- Líneas de flujo de calor en la placa rectangular previamente

..definida

12
CONCLUSIONES

Se desarrolló un software en Matlab, que permitió simular la solución de la ecuación

diferencial en derivadas parciales de segundo orden de tipo hiperbólico, que permite
analizar el comportamiento de la transferencia de calor en las dimensiones dadas,
considerando el comportamiento de los nodos de los bordes, las esquinas y los internos.

Se desarrolló una interfaz gráfica para hacer más interactivo el proceso de análisis bajo
corridas sucesivas, permitiendo modificar los parámetros de simulación y las características
del material bajo estudio, de una manera más práctica.

El software permite variar dimensiones de la lámina, las temperaturas de frontera, externa

e interna (condición inicial), número de iteraciones, el coeficiente de transferencia de calor
por convección, propiedades térmicas y físicas del material, proporcionando los resultados
por pantalla, pudiendo ser almacenados en medios digitales.

Esta herramienta facilita el análisis del proceso de transferencia de calor en láminas con
diversas propiedades térmicas y físicas, permitiendo el estudio con distintos número de
iteraciones lo cual facilita la visualización del fenómeno físico de la propagación de calor en
un material, pudiendo establecer comparaciones para apoyar la toma de decisiones en
proceso industriales.

REFERENCIAS BIBLIOGRAFICAS
 Incropera, F., DeWitt, D. 1999. “Fundamentos de transferencia de calor”. Pp. 173-188, Pp
248-257. Prentice Hall. Cuarta Edicion.
 Mathworks, Inc.1997.“Partial Differential Equation Toolbox User’s Guide”,[On line].
Disponible
en:htpp/www.mathenath.chalmers.se/~mohammad/teaching/PDE1_TM/pde1_lab/la
b2/toolboxguide.pdf
 Zozaya, A. “Método de las diferencias finitas y su aplicación a problemas de electrostática”.
Universidad de Carabobo. Facultad de Ingeniería. 2004. pp.
 http://avalon.utadeo.edu.co/dependencias/publicaciones/alimentica2/libro4a.pdf

13