Clase EDP

DIFERENCIAS FINITAS: INTRODUCCIN
En la bsqueda de una descripcin cualitativa de un determinado

fenmeno fsico, normalmente el ingeniero plantea un sistema de
ecuaciones diferenciales ordinarias o parciales, vlidas para
determinada regin, e impone sobre dicho sistema condiciones
apropiadas de frontera e iniciales.
Modelo matemtico completo
Aparece la mayor dificultad,
pues solamente la forma ms
simple de ecuaciones, con
fronteras
geomtricamente
triviales, puede ser resuelta en
forma exacta con los mtodos
matemticos disponibles.
DIFERENCIAS FINITAS: INTRODUCCIN

Si se desea resolver este problema
sin
software,
es
necesario
presentarlo
de
una
manera
algebraica.
Proceso de discretizacin: el
conjunto infinito de nmeros que
representan la funcin o funciones
incgnitas
en
el
continuo,
es
reemplazado por un nmero finito de
parmetros incgnita, y este proceso
requiere alguna forma de aproximacin
Entre las diferentes formas de discretizacin posibles (elementos
finitos, volmenes finitos, etc.), una de las ms simples es mediante el
Mtodo de Diferencias Finitas.
DIFERENCIAS FINITAS: VALOR FRONTERA

En una solucin por este mtodo, las derivadas son reemplazadas por
aproximaciones en diferencias finitas, convirtiendo entonces un
problema de ecuaciones diferenciales en un problema algebraico
fcilmente solucionable por medios comunes (especialmente
matriciales).
VALOR EN LA FRONTERA
Consideremos el problema de encontrar la funcin (x) que satisface
la ecuacin diferencial:
Sujeta a las condiciones de frontera

Las ecuaciones anteriores
son
utilizadas
para
la
descripcin
analtica
de
muchos procesos fsicos, por
ejemplo:
Conduccin de calor a
travs de una pared plana
(TQ en 1-D)
Flujo
en
canales
y
tuberas
Deflexin transversal de
cables
Deformacin axial de
barras (ver Figura).
Entre otros

En la primera condicin de frontera, aplicada en x = 0, el valor de la
funcin (x) se especifica como0, tal como se muestra en la
siguiente ecuacin:
=0
en
x = x0
Una condicin de frontera de este tipo se denomina condicin de

frontera Dirichlet. En la segunda condicin, aplicable a la condicin
remanente de la frontera x = L, el valor de la funcin corresponde a la
ecuacin:
Este tipo de condicin de frontera se denomina condicin de frontera

Neumann .
DIFERENCIAS FINITAS EN 1-D

Si
queremos
resolver
el
problema de valor de frontera
unidimensional
presentado
anteriormente, por el mtodo de
diferencias finitas (ver Figura),
Figura
se
debe
comenzar
por
diferenciar
la
variable
independiente x, construyendo
un conjunto (o grilla o malla de
L+1 puntos discretos igualmente
espaciados (ver Figura).
Figura

APROXIMACIN DE DERIVADAS MEDIANTE DIFERENCIAS FINITAS
Forma alternativa para obtener aproximaciones de diferencia. Permite
deducir:
Frmulas de diferencia sistemticamente
Trminos de error de truncamiento.
Se pueden obtener las aproximaciones de diferencia hacia atrs,
centrada y hacia adelante.
La serie de Taylor para una funcin f evaluada en Xi+1 es:
f ' ' ( xi )
f ( xi 1 ) f ( xi ) f ' ( xi )( xi 1 xi )
( xi 1 xi ) 2 ... Rn
2!
Truncando en el trmino de la primera derivada y realizando los

cambios pertinentes se obtiene:
f ( xi 1 ) f ( xi )
'
f ( xi )
SOLUCIN DE EDP ELPTICAS

Una de las ecuaciones elpticas ms
representativas es la ecuacin de Laplace:
2T 2T
2 0
2
x
y
2T 2T
2 f x, y
2
x
y
Ecuacin empleada cuando NO

HAY fuentes o prdidas de calor
dentro del dominio. (ec. 1)
Ecuacin empleada cuando HAY
fuentes o prdidas de calor dentro
del dominio. (ec. 1)
Para la solucin por diferencias finitas se transforma la primera

ecuacin en una ecuacin algebraica de diferencias

Analicemos la transferencia de calor en una placa,
aplicando la Ley de Fourier con flujo en el eje x y en el
eje y :
Placa plana en donde se

da la transferencia de
calor

En Diferencias Finitas se trata el problema como una
malla de puntos discretos:
Malla de puntos
discretos
correspondiente a
la placa calentada
El anlisis de diferencias centradas presenta las

siguientes ecuaciones con respecto a la
variacin de la temperatura en las dimensiones x
y y:
2
x
2
Ti 1, j 2Ti , j Ti 1. j
x 2
2T Ti , j 1 2Ti , j Ti , j 1
2
y
y 2

Sustituyendo las ecuaciones 3 y 4 en la ecuacin 1,
asumiendo una malla en la que x = y y
reagrupando trminos se obtiene :
Ti 1, j Ti 1, j Ti , j 1 Ti , j 1 4Ti , j 0
Ecuacin Laplaciana en diferencias;
satisface todos los puntos interiores de la
placa.

De manera prctica se analiza el caso de la placa
calentada con condiciones de frontera fijas.

Aplicando la ecuacin Laplaciana a cada nodo:
Para el nodo (1,1)
T2,1 T0,1 T1, 2 T1, 0 4T1,1 0

Observando la figura se puede ver que To,1 =75 y T1,0 = 0,
por lo tanto la ecuacin se reduce a:
T2,1 T1, 2 4T1,1 75

De la misma manera se pueden obtener las ecuaciones en los otros 8 nodos
para obtener:
4T1,1 1T1, 2 0T1,3 1T2,1 0T2, 2 0T2,3 0T3,1 0T3, 2 0T3,3 T0,1 T1,0
1T1,1 4T1, 2 1T1,3 0T2,1 1T2, 2 0T2,3 0T3,1 0T3, 2 0T3,3 T0, 2
0T1,1 1T1, 2 4T1,3 0T2,1 0T2, 2 1T2,3 0T3,1 0T3, 2 0T3,3 T0,3 T1, 4
1T1,1 0T1, 2 0T1,3 4T2,1 1T2, 2 0T2,3 1T3,1 0T3, 2 0T3,3 T2,0
0T1,1 1T1, 2 0T1,3 1T2,1 4T2, 2 1T2,3 0T3,1 1T3, 2 0T3,3 0
0T1,1 0T1, 2 1T1,3 0T2,1 1T2, 2 4T2,3 0T3,1 0T3, 2 1T3,3 T2, 4
0T1,1 0T1, 2 0T1,3 1T2,1 0T2, 2 0T2,3 4T3,1 1T3, 2 0T3,3 T4,1 T3,0
0T1,1 0T1, 2 0T1,3 0T2,1 1T2, 2 0T2,3 1T3,1 4T3, 2 1T3,3 T4, 2
0T1,1 0T1, 2 0T1,3 0T2,1 0T2, 2 1T2,3 0T3,1 1T3, 2 4T3,3 T4,3 T3, 4

Lo cual se puede llevar a un sistema matricial de la
forma:
0
1
0
0
0
4 1
1 4 1
0
1
0
0
0
1 4 0
0
1
0
0
0 4 1
0
1
1
0
1
0
1 4 1
0
0
1
0
1 4 0
0
0
0
0
1
0
0 4
0
0
0
0
0
0
0
0
1
0
0
1
1
0
0
0
0
0
1
0
1
0
0
0
0
0
1
0
4 1
1 4
T1,1
T
1, 2
T1,3

T
2
,
1

T2, 2

T
2 ,3
T
3,1
T3, 2
T
3, 3
Ax b x A 1b
0 75
75
75 100
100
0 50
50
100 50

MTODO DE LIEBMANN
Es posible obtener sistemas de ecuaciones muy grandes, en los
cuales hay un nmero significativo de ceros (matices dispersas), en
estos casos no es eficiente aplicar mtodos tradicionales de
solucin matricial puesto que se gasta gran cantidad de memoria de
la computadora almacenando ceros, por esta razn los mtodos
aproximados representan un mejor procedimiento para obtener
soluciones. El mtodo ms comnmente empleado es el de GaussSeidel , el cual cuando se aplica a EDP tambin se conoce como
Mtodo de Liebmann.

Partiendo de esta tcnica se puede llegar a la siguiente
ecuacin:
Ti , j
Ti 1, j Ti 1, j Ti , j 1 Ti , j 1
4
La anterior ecuacin se resuelve de manera

iterativa para j=1 hasta n y i=1 hasta m

Algunas veces se utiliza la sobrerrelajacin para
acelerar la velocidad de convergencia, aplicando la
siguiente frmula despus de cada iteracin:
nuevo
i, j
nuevo
i, j
1 T
anterior
i, j
Donde
= Factor de ponderacin con valores entre 1 y 2

Igual que en el mtodo de Gauss-Seidel, las iteraciones
se repiten hasta que los valores absolutos de todos los
errores relativos porcentuales (a)i,j estn por debajo de
un criterio pre-especificado de terminacin s, dichos
errores porcentuales se estiman mediante la siguiente
ecuacin:
a i, j
anterior
Ti ,nuevo
T
j
i, j
nuevo
i, j
100%

CONDICIONES CON DERIVADA EN LA FRONTERA
Un caso especfico de condiciones de frontera es en donde se da la
derivada, conocida como condicin de frontera de Neumann.
Siguiendo con el ejemplo de la placa calentada se presenta el
siguiente esquema:
Esquema de un nodo frontera en

el extremo izquierdo de la placa
calentada

La ecuacin para el nodo frontera es la siguiente:
T1, j T1, j T0, j 1 T0, j 1 4T0, j 0

El punto imaginario nos permite incorporar la derivada de la
condicin de frontera en el problema, esto es representar la
primera derivada en la dimensin x en (0,j) por la diferencia
dividida finita:
Ti , j T1, j
2x
T1, j
T
T1, j 2x
x

Remplazando se obtiene:
T
2T1, j 2x
T0, j 1 T0, j 1 4T0, j 0
x
De la misma forma es posible establecer relaciones
similares para las condiciones de frontera con
derivadas en los otros extremos.

FRONTERAS IRREGULARES
En muchas ocasiones los problemas de ingeniera no presentan
geometras regulares como en el caso de la placa estudiada
anteriormente. La figura 4 muestra el caso de una placa calentada
con una frontera de forma irregular.
Malla de una placa

calentada con una frontera
de forma irregular.

Aproximacin de las primeras derivadas en la direccin x
i 1,i
Ti , j Ti 1, j
1x
i ,i 1
Ti 1, j Ti , j
2 x
De las primeras derivadas se obtienen las segundas

derivadas
T
T
2
T
T x i ,i 1 x

2
1x 2 x
x
x x
2

Sustituyendo las derivadas en la ecuacin anterior y
agrupando trminos se obtiene:
T
2
2
2
x
x
2
Ti 1, j Ti , j
Ti 1, j Ti , j
1 1 2 2 1 2
De la misma manera se obtiene una expresin para la

dimensin y
2T
2
2
2
y
y
Ti , j 1 Ti , j
Ti , j 1 Ti , j
1 1 2 2 1 2

Sustituyendo las ecuaciones en la ecuacin diferencial
inicial (Laplace) se obtiene:
2
0 2
y
Ti , j 1 Ti , j
Ti , j 1 Ti , j
2
2
1 1 2 2 1 2 x
Ti 1, j Ti , j
Ti1, j Ti , j
2
2
1
2
1 1
Esta ecuacin anterior se aplica a cualquier nodo que

sea adyacente a una frontera irregular tipo Dirichlet.
EDP Parablicas
Introduccin
Las
ecuaciones
parablicas a diferencia
de
las
elpticas,
determinan como una
incgnita vara tanto en el
espacio como en el
tiempo. Para la Ecuacin
de conduccin de calor
(caso de problemas de
propagacin), se toma
un diferencial de una
barra larga, delgada y
aislada.
q( X ) * Y * Z * t q( X X ) * Y * Z * t X * Y * Z * * C * T
Dividiendo por XYZ, aplicando Ley de Furier

y segn la definicin por lmite:
2T T
K 2
X
t
EDP Parablicas
Mtodos Explcitos o de diferencias progresivas
A travs de los
mtodos explcitos se
calculan los valores
en cada nodo para un
tiempo
posterior,
basndose en los
valores presentes del
nodo y sus vecinos.
EDP Parablicas
Las
ecuaciones
parablicas
estn
temporalmente
abiertas
en
los
extremos
mientras
que
las
elpticas
estn acotadas en
todas
las
dimensiones.
EDP Parablicas
Para la ecuacin de conduccin de calor:

2U U
2
X
t
Se aproxima la primera y segunda derivada por diferencias

finitas, hacia delante y centradas respectivamente.
2U UiL1 2UiL UiL1
2
X
X 2
U UiL1 UiL
t
t
UiL 1 UiL
Donde:
* t
L
L
L
*
(
U
2
*
U
U
i
1
i
i
1)
2
(X )
* t
2
( X )
EDP Parablicas
Mtodos Explcitos: Convergencia y estabilidad
Convergencia: Conforme a X y t tienden a cero, los

resultados de la tcnica por diferencias finitas se
aproximarn a la solucin verdadera.
Estabilidad: Los errores en cualquier etapa del clculo
no se amplifican, sino que se atenan conforme avanza
el clculo.
El mtodo es convergente y estable si 1/2 o
1 X 2
t
2
Se tendr un valor ptimo 1/6 al minimizar los errores
de truncamiento
EDP Parablicas
Mtodos Explcitos: Derivada en las condiciones de
frontera
En la frontera izquierda: i=0
U U1L U L1
dx
2x
U L1 U1L 2x
Para este punto

siguiente ecuacin
U
L 1
0
U
dx
en
U U mL 2 U mL
dx
2x
la
U (1 2)U U
L
1
L
0
En la frontera derecha: i=m+1
L
1
Se sustituye la aproximacin
en trminos de la 1ra
derivada
U
L
L
U 0L 1 U1L 2x
(1 2)U 0 U1
x
U mL U mL 2 2x
U
dx
Para este punto en la

siguiente ecuacin
UmL 11 UmL (1 2)UmL 1 UmL 2
Se sustituye la aproximacin
en trminos de la 1ra
derivada
U
UmL 11 UmL (1 2 )UmL 1 UmL 2x
EJEMPLO. EDP PARABLICA [2]

METODO EXPLICITO
Calcular la distribucin de temperatura en una barra larga y

delgada que tiene una longitud de 10 cm. K=0.835cm2/s.
T
2T
k 2
t
x
x 2cm
t 0.1s
t 0 : T 0 para 0 x 10
t 0 : T 100 en X 0 y T 50 en X 1
t
k 2
x
0.020875
EJEMPLO. EDP PARABLICA
Til1 Til Til1 Til1

nodo tiempo
0.1
X=2
0.2
0.1
X=4
X=6
X=8
Ecuacin
T1l 0 0.020875 0 2 0 100
T12 2.0875 0.020875 0 2 2.0875 100
T2l 0 0.020875 0 2 0 0
Resultado
2.0875
4.0878
0
0.2
T22 0 0.020875 0 2 0 2.0875
0.043577
0.1
T3l 0 0.020875 0 2 0 0
0.2
T32 0 0.0208751.0438 2 0 0
0.021788
0.1
T4l 0 0.020875 50 2 0 0
1.0438
0.2
T42 1.0438 0.020875 50 21.0438 0
2.0439
EDP Parablicas
Mtodo Implcito Simple
Aunque utilizan algoritmos ms complicados que los mtodos
explcitos, mejoran los problemas de estabilidad y no excluyen
informacin de importancia para la solucin.
La derivada espacial se aproxima en un nivel de tiempo posterior.
EDP Parablicas
Para el ejemplo de la barra visto anteriormente,
la segunda derivada se aproxima mediante:
2T Ti l 11 2Ti l 1 Ti l11
2
2
x
(x)
que tiene una exactitud de segundo orden.
Cuando esta ecuacin se reemplaza en la EDP
original, resulta una ecuacin con varias
incgnitas que no puede resolverse como en el
mtodo explcito.
EDP Parablicas
El sistema debe resolverse simultneamente pues con las condiciones de
frontera, las formulaciones implcitas dan como resultado un conjunto de
m ecuaciones lineales algebraicas con el mismo nmero de incgnitas. As,
el problema se reduce a la solucin de un sistema de ecuaciones
simultneas en cada punto en el tiempo.
Ti l 11 2Ti l 1 Ti l 11 Ti l 1 Ti l
k
2
(x)
t
que se puede expresa como:
Ti l11 (1 2 )Ti l 1 Ti l11 Ti l
Donde:
t
(x) 2
Esta ecuacin se aplica a todos los nodos interiores, excepto al primero y

al ltimo de los nodos, los cuales deben modificarse para considerar las
condiciones de frontera.
EDP Parablicas
Para el extremo izquierdo de la barra (i=0):
T0l 1 f o (t l 1 )
Donde f (t l 1 ) es una funcin que describe cmo cambia la temperatura

o
con el tiempo de la frontera.
Sustituyendo en la ecuacin de diferencias, se obtiene la ecuacin para el
primer nodo interior:
(1 2 )T1l 1 T2l 1 Ti l f o (t l 1 )
EDP Parablicas
De manera similar se obtiene la ecuacin para el ltimo nodo interior
(i=m):
Tml 11 (1 2 )Tml 1 Tml f m 1 (t l 1 )
Donde
f m 1 (t l 1 )
describe los cambios especficos de temperatura en el
extremo derecho de la barra. (i=m+1)

Cuando se escriben las ecuaciones de diferencias para todos los nodos,
se obtiene el sistema de ecuaciones a resolver. El mtodo tiene la ventaja
de que el sistema es tridiagonal.
EJEMPLO . EDP PARABLICA [2]

METODO IMPLICITO SIMPLE
Resolver el ejemplo anterior por este mtodo
1 2 T1l1 T2l1 T1l f0 t l1

Nodos interiores T l1 1 2 T l1 T l1 T l
i1
i
i1
i
Ultimo nodo
Tml11 1 2 Tml1 Tml fm1 t l1
Primer nodo
T1l 2.0875
l
0.020875 1.04175
0.020875
0
0
T
0
0.020875 1.04175
0.020875 T3l
0
0
0
0.020875 1.04175 T4l 1.04375
1.04175
0.020875
EJEMPLO . EDP PARABLICA

T1l
2.0047
l
0.0406
T
2
0
.
0209
T3
Tl
1.0023
4
1.04175
0.020875
0
0
0.020875 1.04175
0.020875
0
0
0.020875 1.04175
0.020875
0
0
0.020875 1.04175
T12
3.9305
2
0.1190
T
2
0
.
0618
T3
T2
1.9653
4
T12
4.09215
2
0.04059
T
2
0
.
02090
T3
T2
2.04069
4
EDP Parablicas
Aunque este mtodo es estable y convergente, presenta una
deficiencia:
la aproximacin en diferencias temporal tiene una
exactitud de primer orden; y la aproximacin en diferencias espacial
tiene una exactitud de segundo orden.
Adems, hay un lmite de
exactitud para el uso de pasos de tiempo grandes.
El mtodo de Richardson tiene una exactitud de segundo orden para el
espacio y para el tiempo, pero presenta serios problemas de
estabilidad.
El mtodo conocido como Crank- Nicholson ofrece un esquema
implcito que tiene una exactitud de segundo orden para el espacio y
para el tiempo y es incondicionamente estable.
EJERCICIOS
DE APLICACIN
EDP Parablicas
El Mtodo de Crank - Nicholson
Se desarrollan aproximaciones por diferencias en el punto medio del
incremento del tiempo.
As, la primera derivada temporal, para el caso de la barra, se aproxima en

tl+1/2 por:
l 1
l
T Ti Ti
t
t
EJERCICIOS
DE APLICACIN
EDP Parablicas
La segunda derivada en el espacio puede determinarse en el punto medio
promediando las aproximaciones por diferencias al principio (tl) y al final
(tl+1) del incremento del tiempo:
2T 1 Ti l1 2Ti l Ti l 1 Ti l11 2Ti l 1 Ti l 11

2
2
x
2
(x)
( x ) 2
EDP Parablicas
Sustituyendo y reagrupando:
Ti l 11 2(1 )Ti l 1 Ti l11 Ti l 1 2(1 )Ti l Ti l 1

Se determinan las condiciones de frontera
l 1
0
l 1
f o (t )
l 1
m 1
l 1
f m 1 (t )
para obtener versiones de la ecuacin de diferencias para los nodos

interiores primero y ltimo.
Para el primer nodo:
2(1 )T1l 1 Ti l 11 f o (t l ) 2(1 )T1l T2l f o (t l 1 )

Para el ltimo nodo:
Tml 11 2(1 )Tml 1 f m 1 (t l ) 2(1 )Tml Tml 1 f m 1 (t l 1 )
EJEMPLO . EDP PARABLICA [2]

METODO DE CRANK-NICHOLSON
Resolver ejemplo anterior por este mtodo
2.04175
0.020875
0
0
0.020875
2.04175
0.020875
0
0
0.020875
2.04175
0.020875
0
0
0.020875
2.04175
T1l
2.0450
l
0.0210
T2
T3
0.0107
Tl
1.0225
4
T1l
4.175
l
0
T2
0
T3
Tl
2.0875
4
EJEMPLO . EDP PARABLICA
2.04175
0.020875
0
0
0.020875
2.04175
0.020875
0
0
0.020875
2.04175
0.020875
0
0
0.020875
2.04175
T12
4.0073
2
0.0826
T
2
0
.
0422
T
3
T2
2.0036
4
T12
8.1801
2
0.0841
T2
0
.
0427
T3
T2
4.0901
4
Comparacin de los mtodos para la

solucin EDP parablicas
Explcito
Implcito
CrankNicolson
Solucin directa
Sistema de
ecuaciones
Sistema de
ecuaciones
Condicionalmente
estable
Incondicionalmente
estable
Incondicionalmente
estable
Segundo orden en
espacio O(x2) y
primer orden en
tiempo O(t)
Segundo orden en
espacio O(x2) y
primer orden en
tiempo O(t)
Segundo orden en
espacio y en tiempo
O(x2+t2)
EDP Parablicas
Ecuaciones parablicas en 2 dimensiones
La ecuacin parablicas que representan fenmenos de

propagacin y difusin, para dos dimensiones (X,Y) se
podr expresar: U
2U 2U
2
2
t
De igual forma se aproxima por diferencias finitas.

En el caso bidimensional para el mtodo explicito se
tendrn los siguientes criterios para la convergencia y
estabilidad.
1 X 2 Y 2
t
8
* t 1
2
(X )
4
EDP Parablicas
Ecuaciones parablicas en 2 dimensiones
Las tcnicas implcitas brindan mayor estabilidad pero

por para dos y tres dimensiones se obtienen sistemas
de mxn ecuaciones simultaneas, perdindose la
propiedad en la matriz de ser tridiagonales.
Se propone entonces otros mtodos (Mtodos de
separacin), como el esquema IDA, mtodo implcito de
direccin alternante
EDP Parablicas
Ecuaciones parablicas en 2 dimensiones: Esquema IDA
Resuelve ecuaciones diferenciales parablicas en 2

dimensiones espacial es usando matrices tri-diagonales
y para ello cada incremento se hace en dos pasos. En el
primero
Explicito
UiL,j1/ 2 UiL, j
t / 2
UiL1, j 2UiL, j UiL1, j
( X ) 2
UiL,j1/12 2UiL,j1/ 2 UiL,j 11/ 2

(Y )2
Para una malla cuadrada, la ecuacin anterior se expresa:

UiL,j 11/ 2 2(1 )UiL,j1/ 2 UiL,j1/12 UiL1, j 2(1 )UiL, j UiL1, j
EDP Parablicas
Ecuaciones parablicas en 2 dimensiones: Esquema IDA
En el segundo
Implcito
TiL, j1 Ti L, j1/ 2
t / 2
La
TiL1,1j 2TiL, j1 TiL1,1j
(X )2
ecuacin
entonces:
se
TiL, j11/ 2 2TiL, j1/ 2 Ti L, j11/ 2

( Y ) 2
expresa
TiL1,1j 2(1 )TiL, j1 TiL1,1j TiL, j11/ 2 2(1 )TiL, j1/ 2 TiL, j11/ 2
EJEMPLO. EDP PARABLICA [2]

METODO IDA
Utilice el mtodo IDA para encontrar la temperatura de la placa

suponga que la temperatura de la placa es cero en el tiempo 0 y
que se llevan instantaneamente a los valores de la figura. Emplee
un tamao de paso de 10s. K = 0.835cm2/s
100C
Placa de 40 x 40 cm
75C
50C
x = 10
= 0.0835
0C
2T 2T
T
k
2
2
t
y
x

(3,1)
(3,2)
(3,3)
(2,1)
(2,2)
(2,3)
(1,1)
(1,2)
(1,3)
Primer paso:
Til,j11/ 2 21 Til,j 1 / 2 Til,j11/ 2 Til, j1 21 Til, j Til, j1

Segundo paso:
Til11, j 21 Til,j 1 Til11, j Til,j11/ 2 21 Til,j 1 / 2 Til,j11/ 2

Primer paso:
Para los nodos (1,1) (1,2) y (1,3) se plantea el sistema
2.167
0.0835
0.0835
2.167
0.0835
1/ 2
T1,1
6.2625
1/ 2
0.0835 T1, 2 6.2625

14.6125
2.167 T11,3/ 2
Segundo paso:
Para los nodos (1,1) (2,1) y (3,1) se plantea el sistema
2.167
0.0835
0.0835
2.167
0.0835
1
T1,1
13.0639
1
0.0835 T1, 2 0.2577

8.0619
2.167 T11,3

Primer paso:
T31,/12
T31,/22
T31,/32
T21,/12
1 / 2
T1,1
T21,/22
T11, /2 2
T21,/32
T11,3/ 2
2.0181 2.2477 6.0256
0.1274 0.2900 4.1291

3.01060 3.2708 6.8692
Segundo paso:
T31,2
T31,3
T21,1 T21,2
1
1
T
T
1, 2
1,1
T21,3
1
T1,3
T31,1
3.7388 4.2359 11.3606
0.4782 0.8238 8.3207

5.5855 6.1683 13.1120
Ecuaciones hiperblicas
Aspectos bsicos
De forma general se puede decir que son las ecuaciones que tienen
variacin en la dimensin espacial y derivada de orden 2 con respecto al
tiempo.
,
Generalmente se emplean para problemas de vibraciones y

propagacin.
Para ilustrar los principios del mtodo de diferencias finitas aplicado a
la resolucin de ecuaciones diferenciales hiperblicas partiremos de una
ecuacin tpica, como la ecuacin de onda, que esta dada por:
2
2u
u
2
( x, t )
( x, t ) 0
2
2
t
x
Diferencias finitas aplicado a la ecuacin de onda
Sujeta a las condiciones de frontera:
u (0,t) = u(l,t) = 0,
para t>0
,
Lo interpretacin
fsica de estas condiciones es que para cualquier
tiempo t la posicin en los extremos es cero.
Las condiciones iniciales son:

u(x,0) = f(x)
u
( x, 0) g ( x)
x
para:
0 xl
Estas condiciones iniciales expresan la posicin y la velocidad inicial de

la onda, las cuales estn dadas por funciones que dependen de la
dimensin espacial x.
Para establecer el mtodo de diferencias finitas, se selecciona un entero
m>0 y el tamao de paso de tiempo k>0. Con h= l/m los puntos de red
(xi, tj) son:
xi = ih,
y
tj = jk
para cada i, = 0,1,2,m y j=0,1 En cualquier punto de red interior
(xi,tj) la ecuacin de onda se transforma en:
2
2u
u
2
( xi , t j )
( xi , t j ) 0
2
2
t
x
El mtodo de diferencias se obtiene usando el cociente de

diferencias centradas en las segundas derivadas parciales, dadas
por:
u ( xi , ti j ) 2u ( xi ,t j ) u ( xi , t j 1 ) k 2 4u
2u
(xj ,t j )
( xi , j )
2
2
4
t
k
12 t
Si se diferencia para la dimensin espacial:
u ( xi 1 , t j ) 2u ( xi ,t j ) u ( xi 1 , t j ) h 2 4u
2u
(x j ,t j )
( i , t j )
2
2
4
x
h
12 x
,
Si se sustituyen
las ecuaciones de diferencias en la original, se
obtiene la ecuacin de diferencias:
u ( xi , t j 1 ) 2u ( xi , t j ) u ( xi , t j 1 )
k
u ( xi 1 , t j ) 2u ( xi , t j ) u ( xi 1 , t j )
h
Si se reemplaza u (x,t) por el trmino de aproximaciones numricas

k
w con sus respectivos subndices (i, j) y se toma
h
La ecuacin queda:
wi , j 1 2 wi , j wi , j 1 2 wi 1, j 2 2 wi , j 2 wi 1, j 0
Esta ecuacin se resuelve para wi,j+1, o sea la aproximacin ms
avanzada para el paso del tiempo, obteniendo:
wi , j 1 2(1 2 ) wi , j 2 ( wi 1, j wi 1, j ) wi , j 1 0
,
De manera grfica se pueden observar los puntos

involucrados en el clculo:
t
tl+1
tl
xi-1
xi
xi+1
Esta ecuacin es aplicable para toda i=1,2,3.m-1 y para j=1,2,3.
Las condiciones de frontera, dan:
w(0,j)= w(m,j)=0, para cada j=1,2,3
,
Y la condicin inicial implica que:

w(i,0)=f(xi), para cada i=1,2,. m-1
Adems el punto w(i,1) se obtiene a partir de la condicin de
velocidad inicial y da:
wi ,1 wi ,0 kg ( xi ) (error 0(k))
Que se deduce de la aproximacin hacia delante de la derivada.
La ecuacin de forma matricial queda:
w1, j 1
w
2, j 1
m 1, j 1
2(1 )
2
0
M
0L L 0
2(1 2 )
O
O
0 L L 0
2
M
2(1 2 )
w1, j

w2, j
M

wm 1, j
w1, j 1
w2, j 1
wm 1, j 1
Si se desea mayor exactitud error 0(k3)

2
2
wi ,1 (1 2 ) f ( xi )
f ( xi 1 )
f ( xi 1 ) kg ( xi )
2
2
ALGORITMO PARA RESOLUCIN DE LA ECUACIN DE ONDA:
2
2u
u
2
(
x
,
t
)
( x, t ) 0
2
2
t
x
Para:
0 x l, 0 t T
Sujeta a las condiciones de frontera:

u (0,t) = u(l,t) = 0,
para
0t T
Y las condiciones iniciales:

u(x,0) = f(x), y
u
( x, 0) g ( x)
x
para:
0 xl
Entrada: extremo l, tiempo mximo T, constante , enteros m >=2, N>=2
Salida: aproximaciones wi,j a u(xi,tj) para toda i=0,.,m y j=0,.,N.
Paso 1.
Tome:
h l/m
k T / N
k / h
Paso 2
Para j=1,.N, tome
Paso 3
Tome:
w 0 ,0 f 0 ;
w m , 0 f l ;
w0, j 0
wm, j 0
Paso 4.
Para i=1,..m-1
(inicialize para t=0 y t=k)
Tome wi,0 = f(i,h);
2
2
wi ,1 (1 2 ) f ( xi )
f ( xi 1 )
f ( xi 1 ) kg ( xi )
2
2
Paso 5. Multiplicacin de matrices:

Para j=1,..N-1
Para i=1,m-1
Tome: wi , j 1 2(1 2 ) wi , j 2 ( wi 1, j wi 1, j ) wi , j 1 0
Paso 6.
Para j = 0, .N
Tome: t=jk;
Para i=0,..m
Tome x=ih;
SALIDA (x,t,wi,j)
Paso 7. FIN
APLICACION 1. EDP ELPTICA [1]

Ecuacin de Laplace
2u 2u
2 0
2
x
y
Condiciones de frontera
0, x 0
0, y 0
u x , y
1, x 1
1, y 1
Regin: [0, 1] [0, 1]

h= 0.25
APLICACION 1. EDP ELPTICA

u i1, j 2u i, j u i1, j u i, j1 u i, j u i, j1
0
2
2
x
y
u i1, j u i1, j u i, j1 u i, j1 4u i, j 0
Para el nodo (1,1)
u(2, 1) + u(0, 1) + u(1, 2) + u(1, 0) 4 u(1, 1) = 0
Para el nodo (3,2)

u(4, 2) + u(2, 2) + u(3, 3) + u(3, 1) 4 u(3, 2) = 0
u(0.25, 0) = u(0, 0.25) = 0

u(1, 0.5) = 1
u(2, 1) + u(1, 2) 4 u(1, 1) = 0
u(2, 2) + u(3, 3) + u(3, 1) 4 u(3, 2) = -1
SOLUCION: u = (1/7, 2/7, 1/2, 2/7, 1/2, 5/7, 1/2, 5/7, 6/7)
h= 0.05
h= 0.25
h= 0.01
APLICACIN 2. CRANK-NICHOLSON [3]

Se quiere solucionar la EDP de transferencia de calor (archivo heat3.m) usando
el mtodo de Crack-Nicholson, en el intervalo de x entre 0 y 10. Para resolver los
sistemas de ecuaciones implcitas en cada tiempo se usa la factorizacin LU y
se usa un x = 0.5 y un t = 0.05
condiciones de frontera en el
condiciones de frontera en el lado
lado derecho (archivo right.m):
izquierdo (archivo left.m):
t
2t
t
2t
T
te
T
te
1)
1)
1 t
1 t
2) T = 3
3) T = 0
2) T = -2;
3) T = 100.0;
4) T = 0;
Condiciones iniciales (archivo heatf.m):

2x 5 x
e
10
1)
T sin
2)
y1 = x; y2 = 10 -x; T = (x < 5.001).*y1 + (x > 5).*y2 ;
T 5x 3)2 sin 2x
4)
T=0
APLICACIN 2. CRANK-NICHOLSON
En la siguiente tabla se relacionan las diferentes soluciones segn las
condiciones de frontera e iniciales usadas
C. F.
derecha
C. F
izquierda
C. Inicial
Figura
Figura 1
Figura 2
Figura 3
Figura 4
Figura 1
Figura 3
Figura 2
Figura 4
APLICACIN 2. CRANK-NICHOLSON
Ahora se quiere solucionar la EDP de calor (archivo heat4.m) con las
condiciones de frontera de Neuman T(0,t) = T(10,t) = 0
q(x) puede ser (archivo q.m):

1) T e x 5
2)
T
2T
k 2 q x
t
x
T x 5 e 2 x 2 2 x 5 e 8 x 7
2
C. Inicial
q(x)
Figura
Figura 5
Figura 6
Figura 7
Figura 8
Figura 5
Figura 7
Figura 6
Figura 8
APLICACIN 3. ECUACIN DE LA ONDA [3]

Solucionar la ecuacin de la onda (archivo wave1.m) en un intervalo en x de o
a 10, con u igual a cero en las dos fronteras usando el mtodo de diferencias
finitas. Tanto x como t son 0.25. En t=0, u=0. d=0.1 y k = 4.
2u
u 2u
2d 2 ku 0
2
t x
t
APLICACIN 3. ECUACIN DE LA ONDA

Solucionar la ecuacin (archivo wave2.m) en un intervalo en x de o a 10, con u
igual a cero en las dos fronteras usando el mtodo de diferencias finitas. Tanto
x como t son 0.25. =5
2u 2u
3
u
0
2
2
t
x
APLICACIN 4. ECUACINES DE TRANSPORTE [4]

Solucionar la ecuacin de transporte
En un intervalo [a,b], con
condiciones iniciales u(x,0) = u_0(x)
con los esquema de diferencias
finitas de (archivo transport.m) :
Euler
Lax-Friedrichs
Lax-Wendroff
Upwind
u
u
a
0
t
x
Euler
Lax-Friedrichs
Lax-Wendroff
Upwind

Clase EDP

Cargado por

Copyright:

Formatos disponibles

Clase EDP

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

Clase EDP

Cargado por

Copyright:

Formatos disponibles

DIFERENCIAS FINITAS: INTRODUCCIN

En la bsqueda de una descripcin cualitativa de un determinado

DIFERENCIAS FINITAS: INTRODUCCIN

DIFERENCIAS FINITAS: VALOR FRONTERA

Sujeta a las condiciones de frontera

DIFERENCIAS FINITAS: VALOR FRONTERA

DIFERENCIAS FINITAS: VALOR FRONTERA

Una condicin de frontera de este tipo se denomina condicin de

Este tipo de condicin de frontera se denomina condicin de frontera

DIFERENCIAS FINITAS EN 1-D

DIFERENCIAS FINITAS EN 1-D

Truncando en el trmino de la primera derivada y realizando los

DIFERENCIAS FINITAS EN 1-D

DIFERENCIAS FINITAS EN 1-D

DIFERENCIAS FINITAS EN 1-D

SOLUCIN DE EDP ELPTICAS

Ecuacin empleada cuando NO

Para la solucin por diferencias finitas se transforma la primera

SOLUCIN DE EDP ELPTICAS

Placa plana en donde se

SOLUCIN DE EDP ELPTICAS

SOLUCIN DE EDP ELPTICAS

El anlisis de diferencias centradas presenta las

SOLUCIN DE EDP ELPTICAS

SOLUCIN DE EDP ELPTICAS

SOLUCIN DE EDP ELPTICAS

Para el nodo (1,1)

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

T2,1 T1, 2 4T1,1 75

SOLUCIN DE EDP ELPTICAS

SOLUCIN DE EDP ELPTICAS

SOLUCIN DE EDP ELPTICAS

SOLUCIN DE EDP ELPTICAS

La anterior ecuacin se resuelve de manera

SOLUCIN DE EDP ELPTICAS

SOLUCIN DE EDP ELPTICAS

SOLUCIN DE EDP ELPTICAS

Esquema de un nodo frontera en

SOLUCIN DE EDP ELPTICAS

T1, j T1, j T0, j 1 T0, j 1 4T0, j 0

SOLUCIN DE EDP ELPTICAS

SOLUCIN DE EDP ELPTICAS

Malla de una placa

SOLUCIN DE EDP ELPTICAS

De las primeras derivadas se obtienen las segundas

SOLUCIN DE EDP ELPTICAS

De la misma manera se obtiene una expresin para la

SOLUCIN DE EDP ELPTICAS

Esta ecuacin anterior se aplica a cualquier nodo que

Dividiendo por XYZ, aplicando Ley de Furier

Para la ecuacin de conduccin de calor:

Se aproxima la primera y segunda derivada por diferencias

Convergencia: Conforme a X y t tienden a cero, los

Para este punto

En la frontera derecha: i=m+1

Para este punto en la

UmL 11 UmL (1 2 )UmL 1 UmL 2x

EJEMPLO. EDP PARABLICA [2]