Diferencias Finitas

1
Tareas Taller de integracin I Tarea trmica

Jos Quintanilla Acevedo1, Nicols Vicencio Mora2
Facultad de Ingeniera, Universidad de Talca, Curic
1
[email protected]
2
[email protected]
Ingeniera en Mecatrnica
Chile
AbstractEn el presente documento se detallaran los pasos
para la modelacin y simulacin de una barra sometida a un
flujo de calor, tarea correspondiente a una de las tareas del
mdulo integracin I de ingeniera en Mecatrnica, en donde se
ver un clculo asociado estado estacionario y luego el mtodo
numrico de diferencias finitas para describir temperaturas en el
tiempo.
1 metro de largo.
1 cm de espesor.
10 cm de ancho.
El problema se expresa en la siguiente figura:
I. INTRODUCCIN
n proceso o sistema en ingeniera puede estudiarse en

forma experimental o en forma analtica, la primera
implica realizacin de pruebas y toma de decisiones y la
segunda implica la realizacin de clculos, el proceso
experimental suele ser caro para sistemas complejos, en
cambio el procedimiento analtico, el cual por lo dems pude
ser numrico, tiene la ventaja que es ms rpido y barato en
algunos casos.
Lo que hay que tener en cuenta es que los resultados
obtenidos analticamente y numricamente estn sujetos a la
exactitud de las suposiciones, de las aproximaciones y de las
idealizaciones tomadas en cuenta para realizar los clculos.
En este caso se modelar un sistema en el que se produce
transferencia de calor, calor que se define como la forma de
la energa que se puede transferir de un sistema a otro como
resultado de la diferencia en la temperatura.
El problema que se tratar en este documento se relaciona con
la transferencia de calor, ciencia que determina el tiempo en
donde hay un intercambio de temperaturas o razones de
transferencia de energa, para el cual se resolver de manera
discreta por medio del mtodo de diferencias finitas.
II. PROBLEMA A RESOLVER
El sistema a resolver corresponde a una placa de plomo que se
encuentra en un ambiente con temperatura de 25 C, a la cual
se le aplica un flujo de calor constante de 15 W en una cara de
la placa.
Se pretende describir la temperatura de una cantidad de puntos
de la placa, y el comportamiento de sta en el tiempo.
Las dimensiones de la placa son:
Fig. 1 Esquema del problema a resolver.
Para abordar este problema se deben reconocer los fenmenos

que se estn produciendo, la transferencia de calor se puede
presentar en 3 formas:
Conduccin: Es un proceso de transmisin de calor basado en
el contacto directo entre los cuerpos, sin intercambio de
materia, por el que el calor fluye desde un cuerpo de mayor
temperatura a otro de menor temperatura que est en contacto
con el primero
Radiacin: La radiacin a la transmisin de calor entre dos
cuerpos los cuales, en un instante dado, tienen temperaturas
distintas, sin que entre ellos exista contacto ni conexin por
otro slido conductor. Es una forma de emisin de ondas
electromagnticas que emana todo cuerpo que est a mayor
temperatura que el cero absoluto.
Conveccin: se caracteriza porque se produce por medio de un
fluido (lquido o gas) que transporta el calor entre zonas con
diferentes temperaturas
Con los datos proporcionados para el problema, se pueden
identificar transferencia de calor por medio de dos formas,
conduccin y conveccin, la primera se produce en el contacto
directo de una fuente de flujo de calor de 15 W en un extremo,
y la segunda es la conveccin con el aire del ambiente a 25C.
III. ENFOQUE ANALTICO

Para la solucin de los distintos anlisis, se deben considerar
las propiedades o parmetros del sistema, las cuales se
2
presentan a continuacin:
Flujo Unidimensional:
se toma en cuenta que las magnitudes de las medidas de la
placa son comparablemente grandes porque el largo de 1
metro es una medida demasiado mayor a su ancho y espesor
de 0,1 y 0,01 metros, por lo cual es entendible que la direccin
sea a lo largo de la placa que ser representada como la
direccin x. Lo que corresponde a dejar la extensa ecuacin
de calor por conduccin as:
Constante Trmica Plomo:

*
Calor especfico Plomo:

[
Densidad Plomo:
[
No hay generacin interna de calor.
Para aletas de seccin transversal uniforme se tiene la

distribucin de temperaturas es:
Difusidad Trmica de plomo:

*
(Cengel, 2007)
Coeficiente de conveccin aire :
*
Donde los elementos de la ecuacin son iguales a :
+
( )
Donde
Flujo entrada al material:

[ ]
es la temperatura de la base de la aleta.
+
La transferencia de calor en la aleta est dado por la siguiente
expresin:
Temperatura aire o ambiente:
Para la estimacin analtica del problema se tomar como

referencia el anlisis de superficies extendidas o bien llamadas
aletas, estas son de uso comn en la prctica para mejorar la
transferencia de calor para incrementar la razn de
transferencia.
Se elige este anlisis ya que el caso es muy cercano a la
aplicacin de estas generalmente solo una superficie tiene
consigo un flujo de calor el cal de debe disipar el cual es
transferido a la aleta, cosa que sucede en el problema, ya que
se tiene un flujo de calor constante hacia la placa.
Como se trata de este tipo de anlisis se deben tener
suposiciones claras, cosa comn en el anlisis de transferencia
de calor, en este tipo de anlisis se presentan las suposiciones
siguientes:
En el anlisis de aletas se considera una operacin
estacionaria sin generacin de calor en la aleta y se supone
que la conductividad trmica del material permanece constante
(Cengel, 2007).
Donde
Tambin
transversal y
analizar.
es una constante y corresponde al rea

es el permetro de la pieza y es el largo a
Con estas relaciones entregadas por la literatura[1], es posible

encontrar una solucin analtica para encontrar temperaturas
en puntos de la placa. Se puede encontrar temperatura en la
base de la aleta:
3
Reemplazando los paramentos del problema antes
mencionados, con un permetro de aleta de 22 cm, y L=1
debido a las dimensiones, se tiene que:
Despejando desde
La temperatura en la base es:
Ya con esta temperatura adquirida se est en condiciones de

encontrar la temperatura en distintos puntos de la barra ya que
se presenta la relacin antes mencionada:
Por ejemplo para el final de la barra con x=1:
(
Temperatura en C
78.75
0.2
36
0.5
26.02
0.8
25.09
Para discretizar esta ecuacin, se debe recurrir a algn mtodo

numrico,en donde se acostumbra el reemplazo de ecuaciones
diferenciales por ecuaciones algebraicas.
Existe un popular mtodo que reemplaza las derivadas por
diferencias, mtodo llamado diferencias finitas, este mtodo
puede ser aplicado para problemas de transferencia de calor en
estado estacionario y de forma diferente para problemas de
rgimen transitorio. Para problemas de estado estacionario se
aplica una doferenciacion del problema en el sentido de las
variables espaciales, que es valida para cualquier instante del
estado estacionario, sin embargo en rgimen transitorio, se
produce una diferenciacin con respecto al tiempo y al
espacio.
Esto quiere decir que en la aleta ms exactamente en su punto

final la temperatura es la temperatura del ambiente.
Si se hace un pequeo barrido de puntos se podr comprobar
que a medida el largo avanza hay menos temperatura.
Distancia "x"
estacionario implica que no hay cambio con el tiempo en

todos los puntos del sistema, al contrario de un anlisis
transitorio que implica la variacin con el tiempo de las
variables, durante la transferencia de calor transitoria, la
temperatura suele variar con el tiempo como con la posicin.
Para los efectos de simulacin se intentara aproximar la
ecuacin de calor unidimensional:
Fig. 2 formulacin en diferencias finitas en r. transitorio.
Ahora para proceder con el mtodo hay que tener en cuenta

que se tiene una segunda derivada, por lo cual se debe realizar
dos veces la diferenciacin.
Se debe aplicar el mtodo en la primera derivada, definiendo
como el nmero del nodo a analizar.
25
(
IV. SOLUCIN NUMRICA

Los problemas de transferencia de calor casi siempre se
clasifican como estacionarios o transitorios, el termino
Ahora cada gradiente se expresa, como funcin de las

temperaturas nodales.
( )
Quedando finalmente para un sistema unidimensional en la
coordenada x, o sea
:
)
(
) (
Al sustituir los dos gradientes en la primera aproximacin

queda lo siguiente:
Ya discretizado el espacio, la ecuacin debe discretearse en el

tiempo, esto se har con el mtodo explicito, el cual representa
a la derivada con respecto al tiempo en forma de diferencia
hacia adelante. El entero se introduce en la diferenciacin
para lo que se debe tomar en cuenta que
. Tomando
en cuenta esto la aproximacin por diferencias finitas para la
derivada respecto del tiempo se expresa como:
Esta ltima formula es explicita ya que las temperaturas

nodales desconocidas se determinan desde las temperaturas
nodales conocidas del tiempo anterior, adems es factorizada
as, por que ms adelante se ver una interesante propiedad de
estabilidad en la cual se deben tener un coeficiente de .
Se definir un espaciamiento de
, un nmero par
que permitir conocer la temperatura en la ltima parte de la
barra.
Como el valor
no est disponible en un principio se debe
aplicar un balance de energa alrededor el nodo 1. Este balance
viene dado por:
Por lo que el balance se expresa de la siguiente forma:

El superndice
se sirve para denotar la dependencia con
respecto al tiempo y la derivada respecto al tiempo se expresa
como la diferencia de las temperaturas asociadas con los
( ), donde
tiempos
(
) y
. Lo que conlleva a que temperaturas deben ser
medidas o muestreadas con un intervalo
de desfase entre
cada tiempo de muestra.
Tomando en cuenta el tiempo la relacin se debe incluir a la
ecuacin final:
Despejando el valor de inicio:

(
Para realizar el cdigo se debe seguir el siguiente

procedimiento:
Quedando finalmente:
Equivalente a:
Ahora despejando la variable de temperatura:
Desde esta ecuacin se define el nmero discreto de Fourier,

que es adimensional:
Evaluando dentro de un ciclo o bucle de iteraciones.
5
Ahora para eleccin de
normalmente se tienen que tener
en cuenta la presicion requerida, sin embargo se debe ser
cuidadoso ya que al hacer la eleccin, se debe tener en cuenta
el cual debe ser elegido mediante un criterio de
estabilidad, para evitar resultados errneos el valor de
debe mantenerse por debajo de cierto limite, el cual depende
de
y algunos de los parmetros.
El criterio se determina requiriendo que el coeficiente
asociado con el nodo de inters en el tiempo anterior mayor o
igual a 0, esto se hace reuniendo todos los trminos que
incluyen .
Entonces para la forma de diferencias finitas de un nodo
unidimensional, que es el caso tratado se debe cumplir que:
(
(
(
Para forzar al lmite el sistema, para describir temperaturas, se

elegir un numero de Fourier de 0.5, para el cual se tendr el
tiempo superior de :
(
)
(
[ ]
TABLA I
RESULTADOS SOLUCIN EXPLICITA FO=1/4.
Lo que significa:
forma continua con el tiempo.

Para mejorar estos resultados y eliminar la deformacin
expuesta se debiera reducir e valor de Fo, se probar un
numero de Fourier reducido a la mitad Fo=1/4, para lo cual se
se usara el mismo espaciamiento en x:
[ ]
Definiendo entonces el nmero de Fourier y el flujo de calor

como sigue:
(s)
0
414
828
1242
1656
2070
2484
25
67.49
88.74
104.67
117.95
129.57
140
25
25
35.62
46.25
56.21
65.50
74.22
25
25
25
27.66
31.64
36.29
41.27
25
25
25
25
25.66
26.99
28.86
25
25
25
25
25
25.17
25.58
25
25
25
25
25
25
25.04
Como se observa el paso de tiempo se reduce a la mitad, no

as las temperaturas en el tiempo, el cambio solo fue que el
algoritmo est siendo ms fino en su respuesta, lo que prueba
que para un nmero pequeo de Fo la precisin de los clculos
mejora, pero se debe aumentar la cantidad de puntos a analizar
si se quiere analizar un lapso de tiempo especfico.
Y la frmula para el primer nodo:

(
V. CONCLUSIN
Arroja los siguientes resultados en cdigo C.

TABLA I
RESULTADOS SOLUCIN EXPLICITA FO=1/2.
(s)
0
828
1656
2484
3312
4140
4968
25
109.99
109.99
152.48
152.48
184.35
184.35
25
25
67.49
67.49
99.36
99.36
125.92
25
25
25
46.25
46.25
67.49
67.49
25
25
25
25
35.62
35.62
48.90
25
25
25
25
25
30.31
30.31
25
25
25
25
25
25
27.66
Se puede observar que el clculo de temperaturas es igual en

tiempos sucesivos para el mismo nodo, esto corresponde a un
vicio del algoritmo al escoger el nmero de Fourier mximo
permitido por el criterio de estabilidad antes mencionado,
comportamiento que no se asemeja al comportamiento real de
transferencia de calor, ya que ante un flujo constante al menos
los nodos ms cercanos a la conduccin deben cambiar de
La tcnica numrica aproximada de diferencias finitas sirve

frente a casos con complejidades geomtricas en problemas
de transferencia de calor, aun as en este documento se llev a
cabo un sistema muy bsico de transferencia de calor, en
primer lugar se conoci el concepto de cmo se tratan los
ejercicios en forma estacionaria, los cuales entregan datos solo
cuando el sistema no cambia sus variables en el tiempo, y
aunque en el caso analtico se estudi la conveccin, fue una
manera de expresar con algn mtodo la transferencia de calor
en estado estable de la barra propuesta y dejar a entrever que
los datos entregados en modo estable muchas veces
corresponden a datos mximos de transferencia o temperatura
que suelen ser tomados en cuenta en clculos de productos o
prestaciones con el simple trabajo de reemplazar valores en
correlaciones que estn plasmadas en los textos.
En segunda instancia se hicieron suposiciones ms estrictas
respecto a la barra, asignndola solo con transferencia de calor
por conduccin, lo cual fue hecho para expresar de manera
ms sencilla el mtodo de diferencias finitas que de todos
modos es aplicable para sistemas donde, no solo exista
conduccin, si no tambin conveccin y radiacin o
combinaciones, solo que habra que reescribir la
diferenciacin en el espacio en caso de un anlisis estacionario
6
o una diferenciacin en el sentido del espacio y el tiempo, que
es lo que se hizo en las ltimas pginas.
Se pudo aprender mucho de la literatura, nunca se haba
mencionado el mtodo de diferencias finitas en los cursos de
termodinmica y transferencia de calor, por lo cual fue una
buena ocasin para conocer el mtodo y saber todo lo que est
detrs de los simuladores mecnicos, como por ejemplo el
software Ansys que a veces trabaja con estos mtodos pero
aplicados en forma no tan solo unidimensional si no que hasta
en tres dimensiones de nodos, por ultimo destacar que todos
los mtodos discretos que se han visto en las distintas tareas
son aplicables, precisamente por su condicin discreta a
sistemas de control y de soluciones tecnolgicas.
.
VI. BIBLIOGRAFA
Cengel, Y. A. (2007). transferencia de calor y masa.
McGRAW-HIL.
incropera, F. (1999). Fundamentos de transferencia de calor.
Pearson.

Diferencias Finitas

Cargado por

Copyright:

Formatos disponibles

Diferencias Finitas

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

Diferencias Finitas

Cargado por

Copyright:

Formatos disponibles

1

Tareas Taller de integracin I Tarea trmica

n proceso o sistema en ingeniera puede estudiarse en

Fig. 1 Esquema del problema a resolver.

Para abordar este problema se deben reconocer los fenmenos

III. ENFOQUE ANALTICO

Constante Trmica Plomo:

Calor especfico Plomo:

No hay generacin interna de calor.

Para aletas de seccin transversal uniforme se tiene la

Difusidad Trmica de plomo:

Donde los elementos de la ecuacin son iguales a :

Flujo entrada al material:

es la temperatura de la base de la aleta.

Temperatura aire o ambiente:

Para la estimacin analtica del problema se tomar como

es una constante y corresponde al rea

Con estas relaciones entregadas por la literatura[1], es posible

La temperatura en la base es:

Ya con esta temperatura adquirida se est en condiciones de

Para discretizar esta ecuacin, se debe recurrir a algn mtodo

Esto quiere decir que en la aleta ms exactamente en su punto

estacionario implica que no hay cambio con el tiempo en

Fig. 2 formulacin en diferencias finitas en r. transitorio.

Ahora para proceder con el mtodo hay que tener en cuenta

IV. SOLUCIN NUMRICA

Ahora cada gradiente se expresa, como funcin de las

Al sustituir los dos gradientes en la primera aproximacin

Ya discretizado el espacio, la ecuacin debe discretearse en el

Esta ltima formula es explicita ya que las temperaturas

Por lo que el balance se expresa de la siguiente forma:

Despejando el valor de inicio:

Para realizar el cdigo se debe seguir el siguiente

Ahora despejando la variable de temperatura:

Desde esta ecuacin se define el nmero discreto de Fourier,

Para forzar al lmite el sistema, para describir temperaturas, se

forma continua con el tiempo.

Definiendo entonces el nmero de Fourier y el flujo de calor

Como se observa el paso de tiempo se reduce a la mitad, no

Y la frmula para el primer nodo:

Arroja los siguientes resultados en cdigo C.

Se puede observar que el clculo de temperaturas es igual en

La tcnica numrica aproximada de diferencias finitas sirve

También podría gustarte