Ecuaciones Diferenciales Parciales Parabolicas

Solucion numerica de ecuaciones en derivadas parciales
Ecuaciones diferenciales parciales parabolicas

Ing. Jes
us Javier Cortes Rosas
M. en A. Miguel Eduardo Gonzalez Cardenas
M. en A. Vctor D. Pinilla Moran *
2011
Resumen
Introducci
on. La ecuaci
on en derivadas parciales parabolicas. Ejemplo de aplicacion.
Conclusiones.
1.
Introducci
on
Uno de los ejemplos m

as representativos de las ecuaciones en derivadas parciales parabolicas es la
propagacion del flujo de calor en una direccion en funcion del tiempo.
De nuevo resulta muy importante establecer que este tipo de solucion esta sujeta de forma muy
importante a las caractersticas fsicas de los elementos del experimento, sobre todo en adelanto
a la introducci
on del nuevo concepto de estabilidad cuya aplicacion pudiera entenderse como la
realizacion de artilugios matem
aticos para manipular los resultados.
2.
La ecuaci
on en derivadas parciales parab
olicas
Las ecuaciones en derivadas parciales parabolicas [?] se presentan en problemas de propagacion. En

este tipo de problemas la soluci
on avanza indefinidamente a partir de valores iniciales conocidos,
satisfaciendo las condiciones de frontera conocidas, conforme la solucion progresa. Este tipo de
dominio de soluci
on, de extremo abierto, se muestra en la figura ??.
En estos modelos la variable dependiente de U y las variables independientes son X y t. La soluci
on
debe satisfacer la ecuaci
on diferencial parcial en todo el dominio abierto, as como las condiciones
iniciales y de frontera. Un ejemplo de ecuacion diferencial de tipo parabolico se encuentra al
considerar el problema unidimensional del flujo transitorio de calor que se define mediante la ecuacion diferencial numero ??
k
2
U (X, t) =
U (X; t)
(1)
t
C P X 2
*
Facultad de Ingeniera, UNAM. Profesores de tiempo completo del Departamento de Matem

aticas Aplicadas de
la Divisi
on de Ciencias B
asicas
Analisis numerico
Figura 1: Esquema general
Donde C es la capacidad calorfica del material, P es la densidad del material y k es la conductividad

termica del material.
Como se mencion
o, la soluci
on U (X, t) se propaga en un plano espacio-tiempo como el mostrado
en la figura ??. Si se considera una red espacio-tiempo como la que se muestra en la figura ??,
la solucion de la ecuaci
on (??) consiste en determinar la temperatura U en cada punto de la red
utilizada.
Para utilizar la red propuesta, en la ecuacion (??) deberan sustituirse las ecuaciones de derivaci
on
parcial numerica de acuerdo a lo establecido por el metodo de las diferencias finitas [?]. Se propone
el uso de las siguientes ecuaciones:
U
1
[Ui,j + Ui,j+1 ] + O(h)
(2)
=
t i,j
t
2 U
1
[ Ui1,j 2Ui,j + Ui+1,j ] + O(h2 )
(3)
=
2
X i,j
(X)2
Realizando la sustituci
on:
1
1
[Ui,j + Ui,j+1 ] =
[ Ui1,j 2Ui,j + Ui+1,j ]
t
(X)2
(4)
Despejando el termino inc

ognita Ui,j+1 resulta:
Ui,j+1 = Ui,j +
t
k
[Ui1,j 2Ui,j + Ui+1,j ]
(X)2 C P
(5)
Analisis numerico
Figura 2: Red planteada
Con la ecuaci
on (??) puede obtenerse la temperatura U en la coordenada Xi en el instante tj+1 en
funcion de las temperaturas en los puntos adyacentes de la red. Las temperaturas iniciales conocidas
y las temperaturas de la frontera suministran valores necesarios para iniciar los calculos que proceden
entonces rengl
on por rengl
on satisfaciendo los puntos extremos de cada renglon, las condiciones de
frontera dadas, hasta que satisfaga en forma aproximada alg
un estado final de temperatura (en
funcion del tiempo), conforme la soluci
on se aproxima a un estado estacionario. Esta progresi
on de
renglon por rengl
on en funci
on del tiempo, que contin
ua indefinidamente, ilustra la naturaleza de
extremo abierto del dominio de soluci
on de una ecuacion en derivadas parciales de tipo parab
olico.
Por otra parte, la soluci
on de la ecuaci
on (??) sera estable y no oscilatoria si:
k
t
0,25
C P (X 2 )
(6)
k
t
0,5
C P (X 2 )
(7)
O bien, la soluci
on ser
a estable si:
La referencia hacia la estabilidad del modelo matematico radica en que si los coeficientes antes
citados no tienen los valores necesarios para cumplir con las condiciones ?? y ?? los terminos de
la ecuacion (??) tendr
an valores tales que no se lograra que la ecuacion completa tienda a un
determinado valor.
Analisis numerico
3.
Ejemplo de aplicaci
on
El problema consiste en determinar la distribucion estacionaria de temperatura a lo largo de una

barra delgada. La barra que se muestra en la figura ??. La barra esta hecha de aluminio y tiene
una longitud de un pie.
Figura 3: Barra del problema
Se supone que est

a perfectamente aislada en todas sus fronteras excepto en el extremo izquierdo.
Las propiedades correspondientes del material son:
k = 0,0370
Btu
s f t oF
Btu
Lb o F
Lb
P = 168 3
ft
C = 0,212
La barra se divide en doce incrementos iguales mediante las trece estaciones mostradas, de acuerdo
a la figura ??.
Inicialmente la barra se encuentra en estado de equilibrio termico a la temperatura de 100o F . La
temperatura en el extremo no aislado es reducida s
ubitamente a 0o F ; en ese instante la distribuci
on
de temperatura en la barra pasa a un estado transitorio. Este estado existe hasta el instante en que
la temperatura en todos los puntos de la barra se aproxime al estado final de equilibrio a 0o F . Las
Analisis numerico
condiciones de frontera establecidas matematicamente son:

U0,t = 0o F
U
X (l, t)
(8)
= 0 fFt
La condicion inicial es:

U (X, 0) = 100o F
(9)
2
k
En este sentido, los valores dados de k, C y P establecen una difusividad termica CP
de 0,00104 fst .
1
Con el n
umero de estaciones mostrado en la figura ??, el valor de X es 12
f t. Con todos estos
valores y haciendo referencia al criterio de estabilidad ya mencionado en la ecuacion ?? para una
solucion estable y no oscilatoria, se observa que el uso de un valor conveniente es t = 1,0 s para
satisfacer el requisito para dicha soluci
on.
En la tabla ?? se muestran las temperaturas en las estaciones 1, 3, 5, 7, 9, 11 y 13 en intervalos de

400s.
Cuadro 1: Temperaturas
Tiempo
100
500
900
1300
1700
2100
2500
2900
4.
1
0,00
0,00
0,00
0,00
0,00
0,00
0,00
0,00
3
28,51
9,13
3,27
1,17
0,42
0,15
0,05
0,01
5
53,49
17,64
6,32
2,26
0,81
0,29
0,10
0,03
7
72,60
24,95
8,94
3,20
1,14
0,41
0,14
0,05
9
85,26
30,56
10,95
3,92
1,40
0,50
0,18
0,06
11
92,16
34,09
12,21
4,37
1,56
0,56
0,20
0,07
13
94,31
35,29
12,64
4,53
1,62
0,58
0,20
0,07
Conclusiones
Si bien los resultados debieron ser resumidos para su presentacion, durante su obtencion es posible
constatar la practicidad del metodo de las diferencias finitas, sobre todo para procesos en los cuales
se desee una gran cantidad de puntos para lograr una buena aproximacion.
1
Referencias
[1] Gerald Curtis F. An
alisis numerico. 2a ed. Mexico 1991.
[2] Iriarte-Vivar Rafael. Metodos numericos. Mexico 1990.
1
Editado por Juan Carlos Marn Hel

u. Junio 2011

Ecuaciones Diferenciales Parciales Parabolicas

Cargado por

Copyright:

Formatos disponibles

Ecuaciones Diferenciales Parciales Parabolicas

Cargado por

Información del documento

Descripción 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

Ecuaciones Diferenciales Parciales Parabolicas

Cargado por

Copyright:

Formatos disponibles

Solucion numerica de ecuaciones en derivadas parciales

Ecuaciones diferenciales parciales parabolicas

Uno de los ejemplos m

Las ecuaciones en derivadas parciales parabolicas [?] se presentan en problemas de propagacion. En

Facultad de Ingeniera, UNAM. Profesores de tiempo completo del Departamento de Matem

Figura 1: Esquema general

Donde C es la capacidad calorfica del material, P es la densidad del material y k es la conductividad

Despejando el termino inc

Figura 2: Red planteada

El problema consiste en determinar la distribucion estacionaria de temperatura a lo largo de una

Figura 3: Barra del problema

Se supone que est

condiciones de frontera establecidas matematicamente son:

La condicion inicial es:

En la tabla ?? se muestran las temperaturas en las estaciones 1, 3, 5, 7, 9, 11 y 13 en intervalos de

Editado por Juan Carlos Marn Hel

También podría gustarte