EDP Parciales - Ecuaciones Elípticas
EDP Parciales - Ecuaciones Elípticas
EDP Parciales - Ecuaciones Elípticas
Parciales.
Prof. Mara Angelica Vega U.
Modulo 14
Observaci
on 15.2 Este tema ha sido desarrollado, usando el texto Metodos
Numericos para Ingenieros. Steven C. Chapra - Raymond P.Canale
16.
Ecuaciones elpticas.
16.1.
La ecuaci
on de Laplace
4x
4x
y el segundo por
4y
4y
se obtiene
(248)
q
q
=0
x y
(249)
T
t
(250)
donde,
cal
qi es el flujo de calor en la direccion de la dimension i [ cm
2 s ]
k es el coeficiente de difusividad t
ermica
es la densidad del material
cm2
s
g
cm3
224
H
CV
donde,
H es el calor (cal)
V es el volumen en cm3
(251)
cal
donde, k 0 es el coeficiente de conductividad t
ermica [ scm
C ]. En ambos
casos k y k 0 son parametros que determinan, que tan bien el material conduce
calor.
La ley de Fourier se conoce tambien como ecuaci
on constitutiva, debido a que proporciona un mecanismo que define las interacciones internas
del sistema. Inspeccionando la ecuacion (250) indica que la ley de Fourier especifica que el flujo de calor perpendicular al eje i es proporcional a gradiente
o pendiente de la temperatura en la direccion i. El signo negativo asegura
que un flujo positivo en la direccion i resulta de una pendiente negativa de
alta a baja temperatura. En la figura 13 se ilustra la representacion grafica
de un gradiente de temperatura, recpordamos que el calor se transfiere hacia
abajo, desde una temperatura alta a una baja, el flujo en a) va de izquierda
225
Figura 13:
T
i
<0
Figura 14:
T
i
>0
se obtiene,
2T
2T
+
=0
x2
y 2
(252)
226
(253)
16.2.
Soluci
on Num
erica
(xm+1, yn+1 )
yj+1
(xi , yj+1 )
yj
(xi , yj )
yj1
(xi , yj1 )
x1
xi1
xi
(x0 , y0 )
xi+1 (xm+1 , y0 ) X
227
16.2.1.
La Ecuaci
on de Laplace.
(254)
(255)
(256)
(257)
Esta relacion, que se satisface por todos los puntos interiores de la placa, se
conoce como ecuaci
on laplaciana en diferencias. Debemos ademas especificar las condiciones de frontera en los extremos de la placa para obtener
una solucion u
nica. El caso mas simple es aquel donde la temperatura en la
frontera es un valor fijo. Este tipo de condicion se denomina condici
on de
frontera de Dirichlet. En la figura siguiente , se tiene una placa con este
tipo de condicion. En esta figura un balance en el nodo (1, 1) es de acuerdo
228
(258)
229
T21
T12
T11 +4T21
T31
T21
+4T31
T11
T21
T22
T32
+4T12
T22
T12
+4T22
T32
T22
+4T32
T31
T13
T12
T23
T33
4T13
T22
T23
T13 +4T23
T32
T23
75
50
75
50
= 175
T33
= 100
+4T33 = 150
(259)
16.2.2.
El m
etodo de Liebmann.
230
(260)
0 + 75 + 0 + 0
= 18,75,
4
(261)
aplicando sobrerrelajaci
on se obtiene,
T11 = 1,5(18,75) + (1 1,5)0 = 28,125
(262)
Para i = 2 , j = 1
T21 =
0 + 28,125 + 0 + 0
= 7,03125
4
(263)
Para i = 3 , j = 1
T31 =
50 + 10,54688 + 0 + 0
= 15,13672
4
(264)
(266)
32,51953 28,12500
|100 % = 13,5 %
32,51953
(267)
La novena iteraci
on da como resultado:
T11 = 4300061
(268)
16.3.
Condiciones en la frontera.
(269)
Observese que para aesta ecuacion se necesita un punto (1, j) que esta fuera
de la placa. Aunque parezca un pun to ficticio sirve para incorporar la derivada de la condicion de frontera en el problema, lo que se logra representando
la primera derivada en la dimension x en (0, j) por la diferencia dividida
finita
T1,j T1,j
T
=
x
24x2
(270)
despejando
T1,j = T1,j 24x
T
x
(271)
Ahora tenemos una relacion para T1,j que incluye la derivada, sustituyendo
en la relacion (269) se obtiene
2T1,j 24x
T
+ T0,j+1 + T1,j1 4T0,j = 0
x
234
(272)