Capitulo VII - Ecuaciones Diferenciales Parciales 2

ECUACIONES DIFERENCIALES PARCIALES
7.1
INTRODUCCIN
Las Ecuaciones Diferenciales Parciales son aplicables a la transferencia de calor

y el mtodo ms empleado es el de CrancK Nicholson, quien emplea los nodos y
las diferencias hacia adelante, hacia atrs y las diferencias centrales para
calcular la temperatura en un tiempo y en un punto del objeto de estudio.
Muchos sistemas que se estudian en diversas reas de la ingeniera, se modelan
mediante las ecuaciones diferenciales parciales. Los modelos resultan de las
abstracciones que se realizan en campos como las transferencias de cantidad
de movimiento, calor y masa donde por lo general una variable dependiente
depende de 2 o ms variables independientes.
Objetivos:
Al finalizar esta unidad, a partir del anlisis del problema de la trasferencia de
calor en un sistema unidimensional, el estudiante estar en la capacidad de
identificar y determinar la ecuacin diferencial parcial EDP y aplicar los mtodos
de numricos para hallar su solucin utilizando el Excel y MatLab considerando
el PVI y PVF.
Contenidos Conceptuales
Modelacin de
fenmenos fsicos.
Aproximacin de las
EDP con ecuaciones
de diferencias.
7.2
Contenidos Procedimentales
Contenidos
Actitudinales
Analiza el problema de Trabaja en equipo

la transmisin de calor Comparte
en placas
informacin
Aplica los mtodos
Demuestra inters
numricos a la solucin en el tema.
de las EDP con PVI y
PVF.
ECUACIN GENERAL DE LA CONDUCCIN DE CALOR
Supngase un cuerpo slido de conductividad trmica k, peso especfico y

calor especfico
Cp
independientes de la temperatura T, en el cual fluye calor en
las tres dimensiones del espacio y puede generar o absorber calor debido a
algn fenmeno, por ejemplo de reaccin qumica.
Al efectuar un balance de calor en un elemento diferencial como el de la figura

de dimensiones
x , y , z se
tiene, de acuerdo con la ley de la continuidad:
Calor que entra

al elemento
diferencial
Acumulacin

de Calor
Donde
Calor que sale
del
elemento
q x y z
diferencial
(7.1)
puede ser positivo o negativo dependiendo de si el calor es
generado o absorbido por unidad de volumen y por unidad de tiempo en el

elemento diferencial.
Figura 60: Balance de Energa en un elemento diferencial.

Los flujos de calor representados por la ley de Fourier se sustituyen en la
ecuacin (7.1) y se tiene:
dT
k y z

dx
x
x x
dT

k x z
dy y y
y
dT
k x y
q x y z
dz
z
z
d
dT
x y z CpT k y z
dt
dx
k x z
dT
dx
k x y
dT
dx
(7.2)
Al dividir miembro a miembro entre x , y , z y hacerlos muy pequeos,

osea
x , y , z 0
, al aplicar la definicin de derivada se obtiene:

o
2T 2T 2T
T
q

2
2
2
t
Cp
Donde se ha sustituido a
q / Cp
(7.3)
con , la cual se llama coeficiente de difusividad
trmica.
7.3
ECUACIN DE CALOR EN DIFERENCIAS FINITAS
Una de las ecuaciones diferenciales parciales ms estudiadas es:

T
2T
2
t
x
(7.4)
Que describe la conduccin de calor en rgimen transitorio en una dimensin, la

difusin unidireccional de masa en rgimen transitorio. Por ejemplo, puede
describir la conduccin de calor en una barra aislada longitudinalmente durante
cierto periodo, tomado a partir de t = 0. La barra se considera suficientemente

delgada de longitud L muy grande en comparacin con su grosor.
Sean los extremos de la barra tomados como x = 0 y x = L. Sean adems
las condiciones siguientes:
T(x, 0) = f(x)
0<x<L
T 0,t g1 t
T L,t g 2 t
con t 0
(7.5)
(7.6)
Figura 61: Barra aislada longitudinalmente.

Estas expresiones conocidas como condiciones de frontera, dan los
valores de la temperatura T de la barra en sus extremos a cualquier tiempo t. La
ecuacin (7.4) y las condiciones (7.5) y (7.6) constituyen un problema de valor
en la frontera (PVF).
7.4
MTODO DE CRANK NICHOLSON

Si se tiene la EDP (7.4):
T
2T
2
t
x
Este mtodo consiste en combinar las aproximaciones de
T t
con
diferencias hacia delante apoyndose en la fila j y la aproximacin con diferencia

hacia
atrs apoyndose en la fila j+1, con los que se obtiene un algoritmo
implcito. Por ejemplo, al aproximarse T

delante y de
2T x 2
en el nodo i , j con diferencias hacia
con diferencias centrales se obtiene:

Ti , j 1 Ti , j
t
Al aproximar T
Ti 1, j 2Ti , j Ti 1, j
x 2
(7.7)
en el nodo (i,j+1) con diferencias hacia atrs y a
con diferencia centrales. Se llega a:

Ti , j 1 Ti , j
t
Ti 1, j 1 2Ti , j 1 Ti 1, j 1
x 2
(7.8)
Luego de sumar las ecuaciones (7.7) y (7.8) y rearreglar resulta:

Ti , j 1 Ti , j
Donde:
Ti 1, j 2Ti , j Ti 1, j Ti 1, j 1 2Ti , j 1 Ti 1, j 1
2
t x 2
(7.9)
2T x 2
Figura 62: Nodos usados en el mtodo de Cranck Nicholson.

Ejemplo de Aplicacin 7.1
Calcule la temperatura como una funcin de x y t en una barra aislada de
longitud unitaria ( pies ), sujeta a las siguientes condiciones iniciales y de frontera:
T x, 0 50sen x
CI :
T 1,t 50 F
CF 2 :
Y con
0 x 1
T 0,t 100 F
CF 1 :
t 0
1pie2 / h .
Solucin :
El problema de condiciones en la frontera (PVF) queda establecido como sigue:
T
2T
x2
t
T x, 0 50 sin x
: PVF
T 0,t 100 F
T L,t 50 F
Como ejemplo dividiremos a la barra en cuatro subintervalos (5 puntos o nodos

en la malla x), entonces
x 1 4 0, 25 .
El fenmeno se estudiar durante 1 hora y
se dividir este tiempo en 100 (101 nodos en la malla t), entonces
t 1 100 0, 01 .
Ntese en cuenta que en el nodo (0,0), la temperatura debera ser
50 sin 0 0
atendiendo a la condicin inicial, mientras que la condicin de frontera T (0,t)

establece que debera ser 100 F. Los puntos que representan stas
caractersticas se llaman puntos singulares; se acostumbra tomar en ellos un
valor de temperatura igual a la media aritmtica de las temperaturas sugeridas
por la condicin inicial y la condicin de frontera correspondientes.
La temeperatura tomada para el nodo (0,0) es 50 F. De igual manera el punto
(4,0),
cuya
temperatura
1 0, 01 / 0, 252 0,16
es
de
50 sin 50 /
. Se conocen tambin los nodos:

T1,0 50 sin 0, 25 35, 3553 F
T2,0 50 sin 0, 50 50 F
T3 ,0 50 sin 0, 75 35, 3553 F
2 25
F.
Calculando
Conociendo en resumen los nodos cuando t = 0:

T0 ,0 50 ; T1,0 35, 3553 ;T2,0 50 ; T3,0 35, 3553 ;T4 ,0 25
Resolviendo la ecuacin (7.9) para los siguientes nodos:

T1,1 T1,0
T0,0 2T1,0 T2 ,0 T0 ,1 2T1,1 T2 ,1

2
T2,1 T2,0
T1,0 2T2,0 T3 ,0 T1,1 2T2 ,1 T3 ,1

2
T3 ,1 T3,0
T2,0 2T3,0 T4 ,0 T2,1 2T3 ,1 T4 ,1

2
Los nodos (0,1); (4,1) son conocidos por las condiciones de frontera, sus valores
son:
T0 ,1 100 F
T4,1 50 F
Reemplazando estos valores y los calculados en las tres ecuaciones anteriores

tenemos:
1,16T1,1 0, 08 T2,1
45, 6985
0, 08 T1,1 1,16 T2 ,1 0, 08 T3,1 47, 6568
0, 08 T2,1 1,16T3,1 39, 6985
Este sistema ser resuelto con el mtodo de Thomas del capitulo II,
obtenindose:
T1,1 42.6093 ; T2,1 46, 6039 ; T3,1 37, 4369
Continuando de forma sucesiva se tienen como resultados finales:
Figura 63: Interfaz grfica del mtodo de Cranck Nicholson.

Ejemplo de Aplicacin 7.2
Resuelva el siguiente PVF:
CI :
CF1 :
CF 2 :
Y con
T x,0 20 F
0 x 1
T 1,t 100 F
t 0
T 0,t 100 F
1pie2 / h , L 1pie , tmax 1h
Solucin :
El problema sigue el procedimiento anterior y obtenemos los siguientes
resultados para un
x 0, 25 ; t 0, 01 :
Figura 63: Resultados para el ejemplo 7.2.

7.5 EJERCICIOS PROPUESTOS
7.5.1 Resuelva el siguiente PVF por los mtodos de explcitos e implcitos para
las condiciones dadas.
C 2C
EDP : 2 1pie / h
t x
CI T (x,0) 200 C L 1pie
0 x L
0
CF
T
(
0
,
t
)
100
C
0
12

PVF 1
t 0
0
T (0,t) 20 C 12 t 60
CF2 T (L,t) 1000 C 0 t 12
T (L,t) 200 C 12 t 60
7.6 REFERENCIAS BIBLIOGRFICAS
NIEVES,
Antonio-
METODOS
NUMERICOS
Aplicados
la
Ingeniera Primera Edicin, Editorial CECSA, pag. 533 - 588, Mxico

1996.
CARRASCO, Luis METODOS NUMERICOS APLICADOS A LA

INGENIERA. Segunda Edicin, Ediciones RFG, pg. 693 - 773, Per
2007.
NAGLE, SAFF, SNIDER- ECUACIONES DIFERENCIALES.
Cuarta
Edicin, Editorial Pearson Education, pag. 576 - 660, Lima Mxico

2005.

Capitulo VII - Ecuaciones Diferenciales Parciales 2

Cargado por

Copyright:

Formatos disponibles

Capitulo VII - Ecuaciones Diferenciales Parciales 2

Cargado por

Información del documento

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

Capitulo VII - Ecuaciones Diferenciales Parciales 2

Cargado por

Copyright:

Formatos disponibles

ECUACIONES DIFERENCIALES PARCIALES

Las Ecuaciones Diferenciales Parciales son aplicables a la transferencia de calor

Analiza el problema de Trabaja en equipo

ECUACIN GENERAL DE LA CONDUCCIN DE CALOR

Supngase un cuerpo slido de conductividad trmica k, peso especfico y

independientes de la temperatura T, en el cual fluye calor en

Al efectuar un balance de calor en un elemento diferencial como el de la figura

tiene, de acuerdo con la ley de la continuidad:

Calor que entra

Calor que sale

puede ser positivo o negativo dependiendo de si el calor es

generado o absorbido por unidad de volumen y por unidad de tiempo en el

Figura 60: Balance de Energa en un elemento diferencial.

Al dividir miembro a miembro entre x , y , z y hacerlos muy pequeos,

, al aplicar la definicin de derivada se obtiene:

con , la cual se llama coeficiente de difusividad

ECUACIN DE CALOR EN DIFERENCIAS FINITAS

Una de las ecuaciones diferenciales parciales ms estudiadas es:

Que describe la conduccin de calor en rgimen transitorio en una dimensin, la

cierto periodo, tomado a partir de t = 0. La barra se considera suficientemente

Figura 61: Barra aislada longitudinalmente.

MTODO DE CRANK NICHOLSON

Este mtodo consiste en combinar las aproximaciones de

diferencias hacia delante apoyndose en la fila j y la aproximacin con diferencia

atrs apoyndose en la fila j+1, con los que se obtiene un algoritmo

implcito. Por ejemplo, al aproximarse T

en el nodo i , j con diferencias hacia

con diferencias centrales se obtiene:

en el nodo (i,j+1) con diferencias hacia atrs y a

con diferencia centrales. Se llega a:

Luego de sumar las ecuaciones (7.7) y (7.8) y rearreglar resulta:

Figura 62: Nodos usados en el mtodo de Cranck Nicholson.

Como ejemplo dividiremos a la barra en cuatro subintervalos (5 puntos o nodos

El fenmeno se estudiar durante 1 hora y

se dividir este tiempo en 100 (101 nodos en la malla t), entonces

Ntese en cuenta que en el nodo (0,0), la temperatura debera ser

atendiendo a la condicin inicial, mientras que la condicin de frontera T (0,t)

. Se conocen tambin los nodos:

T3 ,0 50 sin 0, 75 35, 3553 F

Conociendo en resumen los nodos cuando t = 0:

Resolviendo la ecuacin (7.9) para los siguientes nodos:

T0,0 2T1,0 T2 ,0 T0 ,1 2T1,1 T2 ,1

T1,0 2T2,0 T3 ,0 T1,1 2T2 ,1 T3 ,1

T2,0 2T3,0 T4 ,0 T2,1 2T3 ,1 T4 ,1

Reemplazando estos valores y los calculados en las tres ecuaciones anteriores

Continuando de forma sucesiva se tienen como resultados finales:

Figura 63: Interfaz grfica del mtodo de Cranck Nicholson.

1pie2 / h , L 1pie , tmax 1h

Figura 63: Resultados para el ejemplo 7.2.

Ingeniera Primera Edicin, Editorial CECSA, pag. 533 - 588, Mxico

CARRASCO, Luis METODOS NUMERICOS APLICADOS A LA

NAGLE, SAFF, SNIDER- ECUACIONES DIFERENCIALES.

Edicin, Editorial Pearson Education, pag. 576 - 660, Lima Mxico

También podría gustarte