Metodos Numericos para Ecuaciones Parciales Parabolicas

FACULTAD DE CIENCIAS MATEMÁTICAS
ESCUELA PROFESIONAL DE COMPUTACIÓN CIENTÍFICA
MÉTODOS NUMÉRICOS
Y PROGRAMACIÓN IV
Sesión 1
Ecuaciones Diferenciales Parabólicas
Alumno:
JEREMY KARSEN HOLGUÍN MORI
12140445
Ciudad Universitaria
2019
ECUACIONES DIFERENCIALES
PARCIALES NUMÉRICAS
Sesión 1 : Ecuaciones Diferenciales Parabólicas
1. Solución : ¿Es una E.D.P Hiperbólica?
4uxx + 6uxy + 9uyy = 0 (1)
Para que sea hiperbólica, debe cumplirse lo siguiente:
∆ = b2 − 4ac > 0
Donde :
b=6,a=4,c=9
∆ = b2 − 4ac = 62 − 4 · 4 · 9 = −108 < 0
Por lo tanto la E.D.P no es hiperbólica.
2. Solución : Resolver el siguiente problema.


∂u 2

 ∂t
= 0,9 ∂∂xu2

u(0, t) = 0 = u(1, t), t ∈]0, T ] (2)
u(x, 0) = sin(πx), x ∈ [0, 1]

Use: T = 0,3seg, con h = 0,1m, k = 0,1seg. Usar el método explícito. ¿Es buena la
aproximación?
2
La solucion a este problema es : u(x, t) = e−0,9π t sin(πx)
Daremos solución a este problema con la ayuda del código de Matlab, el cual nos muestra
los siguientes resultados:
1
Figura 1: Código del Método Explícito en MATLAB
Hemos obtenido la aproximación usando los datos que el problema nos pide, ahora anali-
zaremos la solucion exacta y la compararemos con los valores obtenidos cuando T = 0,3seg.
2
Figura 2: Comparación de los resultados obtenidos
Vemos que la aproximación es mala.
Figura 3: Gráfico de la solución aproximada
3
3. Solución :
Resolver el problema 2 usando T =0.3 seg, h =0.1 m, k =0.005 seg, usando el
Método Explícito. ¿La aproximación mejoró?
Figura 4: Aproximación con el Método Explícito
Podemos ver que el coeficiente lambda λ = r disminuyó de 9 a 0 · 45 , lo que nos lleva

a una mejora en la aproximación, pues para garantizar la estabilidad lamda debe estar
entre :
1
0≤λ=r≤
2
4
Figura 5: Gráfico de la Aproximacion por el Método Explícito
Nuestro gráfico también se ve mejor pues la variable temporal se particionó más

(k = 0 · 005 seg) lo que dió paso a que nuestro lambda λ = r mejore, pues:
c2 · k
λ=
h2
Por lo tanto, gracias a λ = r = 0 · 45 , nuestra aproximación mejoró.
4. Solución :
Muestre la distribución de la temperatura en la varilla cuando: t =0.052, t =0.1,
t =0.2. Elija en su criterio el número de nodos en ambas variables, o en su defecto, proponga
una solución satisfactoria para proveer el valor requerido de forma apropiada.
Convenientemente hacemos k =0.002 para encontrar los valores requeridos t = 0.052 ,

0.1 , 0.2.
5
Obtenemos los siguientes resultados para los tiempos requeridos:
Figura 6: Ejecutando el código con los nuevos datos
Figura 7: Valores obtenidos con Matlab
Figura 8: Gráfica para t = 0.052, 0.1, 0.2
6
5. Solución : Resuelva la ecuación del calor

∂u 2

 ∂t
= C 2 ∂∂xu2

u(0, t) = 0 = u(1, t), t > 0 (3)
u(x, 0) = sin(πx) + sin(3πx), x ∈ [0, 1]

Considere: C =1, T = 0,4seg, la cantidad de intervalos, m =10 (partición para x), n =100
(partición para el tiempo). Use el método explícito. ¿Es buena la aproximación?
Vemos que la aproximación es buena pues λ = 0.4, lo cual nos da la certeza

que el resultado converge a la solución exacta.
7
6. Solución : Repita el ejercicio anterior para:
(a) m = 20 (partición para x) y n = 100 (partición para t).
8
(b) m = 20 (partición para x) y n = 200 (partición para t).
9
(c) m = 50 (partición para x) y n = 500 (partición para t).
Las aproximaciones son malas.
10
7. Solución : ¿Qué observa en todos estos casos?
En el ejercicio anterior, note que no solo obtendrá una solución numérica
similar a la exacta cuando λ ≤ 12 .
En el ejercicio anterior todas la soluciones aproximadas divergen, por lo tanto obtenemos

una mala aproximación.
Esto se puede ver reflejado en el valor de lambda λ = r, pues en todos las casos λ > 1.
8. Solución :
Muestre el comportamiento de la temperatura de la varilla del problema 3 en
su punto medio a través del tiempo, elija según su criterio el número de nodos
en ambas variables.
Figura 9: Datos para hallar el comportamiento de la temperatura en el problema 5 en el

punto x =0.5
11

Metodos Numericos para Ecuaciones Parciales Parabolicas

Cargado por

Copyright:

Formatos disponibles

Metodos Numericos para Ecuaciones 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

Metodos Numericos para Ecuaciones Parciales Parabolicas

Cargado por

Copyright:

Formatos disponibles

FACULTAD DE CIENCIAS MATEMÁTICAS

ESCUELA PROFESIONAL DE COMPUTACIÓN CIENTÍFICA

1. Solución : ¿Es una E.D.P Hiperbólica?

4uxx + 6uxy + 9uyy = 0 (1)

Para que sea hiperbólica, debe cumplirse lo siguiente:

∆ = b2 − 4ac = 62 − 4 · 4 · 9 = −108 < 0

Por lo tanto la E.D.P no es hiperbólica.

2. Solución : Resolver el siguiente problema.

Vemos que la aproximación es mala.

Figura 3: Gráfico de la solución aproximada

Figura 4: Aproximación con el Método Explícito

Podemos ver que el coeficiente lambda λ = r disminuyó de 9 a 0 · 45 , lo que nos lleva

Nuestro gráfico también se ve mejor pues la variable temporal se particionó más

Por lo tanto, gracias a λ = r = 0 · 45 , nuestra aproximación mejoró.

Convenientemente hacemos k =0.002 para encontrar los valores requeridos t = 0.052 ,

Figura 6: Ejecutando el código con los nuevos datos

Figura 7: Valores obtenidos con Matlab

Figura 8: Gráfica para t = 0.052, 0.1, 0.2

Vemos que la aproximación es buena pues λ = 0.4, lo cual nos da la certeza

Las aproximaciones son malas.

En el ejercicio anterior todas la soluciones aproximadas divergen, por lo tanto obtenemos

Figura 9: Datos para hallar el comportamiento de la temperatura en el problema 5 en el

También podría gustarte