Semana 10 - ECUACIONES HIPERBOLICAS

ECUACIONES HIPERBÓLICAS
Módulo 9
Métodos Numéricos para Ingeniería

2021 2
Videoconferencia 10
UPN, PASIÓN POR TRANSFORMAR VIDAS

Tema: ECUACIONES HIPERBÓLICAS
Logro
Al término de la sesión, el
estudiante resuelve
ejercicios y/o problemas de
EDP de tipo hiperbólico, en
base a participación activa,
aprendizaje autónomo y
con destreza en el manejo
del ordenador.
Tema: Límites, interpretación y métodos de cálculo
Temario
3. Clasificación 4. Ecuaciones
de una EDP Hiperbólicas
2. Derivadas
5. Método de
parciales de
Diferencias Finitas
una función
1. Introducción
Ecuaciones 6. Algoritmo
Hiperbólicas
1. Introducción
Muchos problemas en ciencia aplicada, física e ingeniería se modelan
matemáticamente mediante EDPs. Por ejemplo, el cambio de
temperatura en una placa y el movimiento vibratorio de una estructura.
A continuación estudiaremos los métodos de las diferencias finitas que

se basan en las fórmulas para aproximar las derivada parciales primera
y segunda de una función. Antes de ellos examinemos algunas
definiciones básicas de las EDPs.
2. Derivadas Parciales de una Función

Dada una función 𝑢 = 𝑢(𝑥, 𝑦) que depende de 𝑥 e 𝑦, la derivada parcial
de 𝑢:
a) Respecto de x en el punto (𝑥, 𝑦) se define por:

𝜕𝑢 𝑢 𝑥 + ∆𝑥, 𝑦 − 𝑢(𝑥, 𝑦)
= lim
𝜕𝑥 ∆𝑥→0 ∆𝑥
b) Respecto de y en el punto (𝑥, 𝑦) se define por:

𝜕𝑢 𝑢 𝑥, 𝑦 + ∆𝑦 − 𝑢(𝑥, 𝑦)
= lim
𝜕𝑦 ∆𝑦→0 ∆𝑦
Una ecuación que tiene derivadas parciales de una función desconocida

𝑢(𝑥, 𝑦) , de dos o más variables independientes, se llama ecuación
diferencial parcial o EDP.
 2u  2u
+ 2 xy 2 + u = 1 Orden 2 y lineal
x 2
y
 3u  2u
+ x 2 + 8u = 5y Orden 3 y Lineal
x 3
y
 2u u
+ xu = x Orden 2 y No lineal
x 2
y
Debido a la amplia aplicación a la ingeniería, centraremos nuestro

estudio en las EDPs lineales de 2° orden. Tales ecuaciones se puede
expresar de la forma general:
𝑨𝒖𝒙𝒙 + 𝑩𝒖𝒙𝒚 + 𝑪𝒖𝒚𝒚 + 𝑫 = 𝟎
Donde A, B y C son funciones de 𝑥 e 𝑦, y D es una función de 𝑥, 𝑦, 𝑢, 𝑢𝑥 y

𝑢𝑦 .
3. Clasificación de una EDP

Dependiendo de las relaciones entre A, B y C las EDP se pueden dividir en tres grupos:
Condición Categoría Ejemplo

Ecuación de Laplace:
𝐵2 − 4𝐴𝐶 < 0 Elíptica 𝑇𝑥𝑥 + 𝑇𝑡𝑡 = 0
∆𝑇 = 0
Ecuación del Calor:
𝐵2 − 4𝐴𝐶 = 0 Parabólica
𝑇𝑡 = 𝑘𝑇𝑥𝑥
Ecuación de la Onda:
𝐵2 − 4𝐴𝐶 > 0 Hiperbólica
𝑌𝑥𝑥 = 𝑐 2 𝑌𝑡𝑡
4. Ecuaciones Hiperbólicas
La Ecuación de Onda
Esta ecuación en DP Hiperbólica, esta dada por:
𝑢𝑡𝑡 𝑥, 𝑡 = 𝑐 2 𝑢𝑥𝑥 (𝑥, 𝑡) para 0 < 𝑥 < 𝑎 y 0 < 𝑡 < 𝑏
Con las condiciones de contorno:
u(0, t ) = 0 y u(a , t ) = 0 para 0  t  b,

u(x ,0) = f (x) para 0  x  a ,
ut (x ,0) = g(x) para 0  x  a ,
5. Diferencias Finitas
Discretización del Dominio
Primero se hace una partición del rectángulo
𝑅= 𝑥, 𝑡 : 0 ≤ 𝑥 ≤ 𝑎, 0 ≤ 𝑡 ≤ 𝑏
En una malla que consta de 𝒏 y 𝒎 rectángulos de tamaño de paso ℎ y 𝑘 para 𝑥 y 𝑡
respectivamente
𝑡
𝑡𝑖
𝑥
𝑥𝑖
Discretización de la variable dependiente
Segundo, se discretiza la condición inicial 𝑡 = 𝑡1 = 0

𝑢 𝑥𝑖 , 𝑡𝑗 = 𝑓 𝑥𝑖
Luego, se hace la discretización de la variable dependiente
𝑢𝑖,𝑗 ≈ 𝑢 𝑥𝑖 , 𝑡𝑗 : 𝑖 = 1, … , 𝑛 + 1; 𝑗 = 1, … , 𝑚 + 1
𝑡 𝑢𝑖𝑗
𝑡𝑖
𝑥
𝑥𝑖
Discretización de la ecuación:
Tercero, se utiliza diferencias centradas basadas en la malla anterior para aproximar
las derivadas 𝑢𝑡𝑡 (𝑥, 𝑡) y 𝑢𝑥𝑥 (𝑥, 𝑡):
𝑢 𝑥, 𝑡 + 𝑘 − 2𝑢 𝑥, 𝑡 + 𝑢(𝑥, 𝑡 − 𝑘)
𝑢𝑡𝑡 (𝑥, 𝑡) ≈
𝑘2
𝑢 𝑥 + ℎ, 𝑡 − 2𝑢 𝑥, 𝑡 + 𝑢(𝑥 − ℎ, 𝑡)
𝑢𝑥𝑥 (𝑥, 𝑡) ≈
ℎ2
Evaluando en los nodos de la malla, se tiene:

𝑢𝑖,𝑗+1 − 2𝑢𝑖,𝑗 + 𝑢𝑖,𝑗−1
𝑢𝑡𝑡 (𝑥𝑖 , 𝑡𝑗 ) ≈
𝑘2
𝑢𝑖+1,𝑗 − 2𝑢𝑖,𝑗 + 𝑢𝑖−1,𝑗
𝑢𝑥𝑥 (𝑥𝑖 , 𝑡𝑗 ) ≈
ℎ2
Reemplazando en la EDP de la onda:

𝑢𝑡𝑡 𝑥, 𝑡 = 𝑐 2 𝑢𝑥𝑥 (𝑥, 𝑡)
Se tiene:
𝑢𝑖,𝑗+1 − 2𝑢𝑖,𝑗 + 𝑢𝑖,𝑗−1 2
𝑢𝑖+1,𝑗 − 2𝑢𝑖,𝑗 + 𝑢𝑖−1,𝑗
2
=𝑐
𝑘 ℎ2
𝑐2𝑘2
𝑢𝑖,𝑗+1 − 2𝑢𝑖,𝑗 + 𝑢𝑖,𝑗−1 = 2 𝑢𝑖+1,𝑗 − 2𝑢𝑖,𝑗 + 𝑢𝑖−1,𝑗
ℎ
Es la ecuación en diferencias finitas que se usará como aproximación

de la EDP de la onda
2
𝑐𝑘
𝑢𝑖,𝑗+1 − 2𝑢𝑖,𝑗 + 𝑢𝑖,𝑗−1 = 𝑢𝑖+1,𝑗 − 2𝑢𝑖,𝑗 + 𝑢𝑖−1,𝑗
ℎ
𝑐𝑘
Se define 𝑟 = , entonces
ℎ
𝑢𝑖,𝑗+1 − 2𝑢𝑖,𝑗 + 𝑢𝑖,𝑗−1 = 𝑟 2 𝑢𝑖+1,𝑗 − 2𝑢𝑖,𝑗 + 𝑢𝑖−1,𝑗
De donde se tiene:
𝑢𝑖,𝑗+1 = 𝑟 2 𝑢𝑖+1,𝑗 − 2𝑟 2 𝑢𝑖,𝑗 + 𝑟 2 𝑢𝑖−1,𝑗 + 2𝑢𝑖,𝑗 − 𝑢𝑖,𝑗−1
𝑢𝑖,𝑗+1 = 2 − 2𝑟 2 𝑢𝑖,𝑗 + 𝑟 2 𝑢𝑖+1,𝑗 + 𝑢𝑖−1,𝑗 − 𝑢𝑖,𝑗−1

Idea gráfica:
𝒖𝒊,𝒋+𝟏
𝒖𝒊−𝟏,𝒋 𝒖𝒊,𝒋 𝒖𝒊+𝟏,𝒋

𝑡𝑗
𝒖𝒊,𝒋−𝟏
𝑥
𝑥𝑖
Discretización de las condiciones de contorno:

1. De 𝑢 0, 𝑡 = 0, se tiene:
𝑢1,𝑗 = 0
2. De 𝑢 𝑎, 𝑡 = 0, se tiene:
𝑢𝑛+1,𝑗 = 0
3. De 𝑢 𝑥, 0 = 𝑓(𝑥), se tiene:
𝑢𝑖,1 = 𝑓𝑖
4. De 𝑢𝑡 𝑥, 0 = 𝑔(𝑥), se tiene:
❑ 𝑢𝑖,2 = 𝑓𝑖 + 𝑘𝑔𝑖 . Al utilizar esta fórmula el error de truncamiento es

notorio. Si se desea utilizar se sugiere que el paso 𝑘 sea pequeño.
2 𝑟2
❑ 𝑢𝑖,2 = 1 − 𝑟 𝑓𝑖 + 𝑘𝑔𝑖 + 𝑓𝑖+1 − 𝑓𝑖−1 . Se sugiere utilizar esta fórmula
2
pues el error de truncamiento es pequeño.

En resumen, la forma cómo trabaja la fórmula de diferencia finita para la

ecuación de la onda se muestra a continuación:
𝑡
𝒖𝒊,𝒋+𝟏
𝒖𝒊−𝟏,𝒋 𝒖𝒊,𝒋 𝒖𝒊+𝟏,𝒋

𝑡𝑗
𝒖𝒊,𝒋−𝟏
𝑥
𝑥𝑖
6. Algoritmo
El algoritmo para solucionar numéricamente la ecuación de la onda por diferencias
finitas es el siguiente:
Paso 1: Ingrese a, b, c, h, k.
Paso 2: Calculo del 𝑛,m y r

𝑎 𝑏 𝑐𝑘
𝑛 = , 𝑚 = , r=
ℎ 𝑘 ℎ
Paso 3: Discretización del dominio

𝑥𝑖 = 𝑖 − 1 ∗ ℎ, para 𝑖 = 1: 𝑛 + 1
𝑡𝑗 = 𝑗 − 1 ∗ 𝑘, para 𝑗 = 1: 𝑚 + 1
𝑓𝑖 = 𝑓 𝑥𝑖 𝑦 𝑔𝑖 = 𝑔 𝑥𝑖 , para 𝑖 = 1: 𝑛 + 1
Paso 4: Discretización de las condiciones de contorno
𝑢1,𝑗 = 0 y 𝑢𝑛+1,𝑗 = 0, para 𝑗 = 1: 𝑚 + 1
𝑢𝑖,1 = 𝑓𝑖 y 𝑢𝑖,2 = 𝑓𝑖 +k𝑔𝑖 , para 𝑖 = 2: 𝑛
Paso 5: Calculo de las aproximaciones
𝒖𝒊,𝒋+𝟏 = 2 − 2𝑟 2 𝑢𝑖,𝑗 + 𝑟 2 𝑢𝑖+1,𝑗 + 𝑢𝑖−1,𝑗 − 𝑢𝑖,𝑗−1

para 𝑖 = 2: 𝑛 y 𝑗 = 2: 𝑚
Paso 6: Imprimir la solución aproximada 𝒖
Paso 7: Graficar la solución aproximada 𝒖

Ejemplo:
Dada la ecuación diferencial de la onda:
𝑢𝑡𝑡 𝑥, 𝑡 = 4𝑢𝑥𝑥 𝑥, 𝑡 para 0 < 𝑥 < 1 y 0 < 𝑡 < 0.5,
con las condiciones de contorno:

u (0, t ) = 0 y u (1, t ) = 0, para 0  t  0.5,
u ( x, 0) = sin(  x) , para 0  x  1,
ut ( x, 0) = 0 , para 0  x  1,
Elabore el algoritmo para solucionar numéricamente esta ecuación con su
respectiva implementación en OCTAVE:, h=0.1 y k=0.05

Solución:
Paso 1: Ingrese a, b, c, h, k. Paso 4: Discretización de las condiciones de contorno
𝑢1,𝑗 = 0 y 𝑢𝑛+1,𝑗 = 0, para 𝑗 = 1: 𝑚 + 1

Paso 2: Calculo del 𝑛, 𝑚 y r
𝑎 𝑏 𝑐𝑘
𝑛 = ℎ , 𝑚 = 𝑘 , r= ℎ 𝑢𝑖,1 = 𝑓𝑖 y 𝑢𝑖,2 = 𝑓𝑖 +k𝑔𝑖 , para 𝑖 = 2: 𝑛
Paso 3: Discretización del dominio Paso 5: Calculo de las aproximaciones

𝒖𝒊,𝒋+𝟏 = 2 − 2𝑟 2 𝑢𝑖,𝑗 + 𝑟 2 𝑢𝑖+1,𝑗 + 𝑢𝑖−1,𝑗 − 𝑢𝑖,𝑗−1
𝑥𝑖 = 𝑖 − 1 ∗ ℎ, para 𝑖 = 1: 𝑛 + 1 Para 𝑖 = 2: 𝑛 y 𝑗 = 2: 𝑚
𝑡𝑗 = 𝑗 − 1 ∗ 𝑘, para 𝑗 = 1: 𝑚 + 1
𝑓𝑖 = sin(𝜋𝑥𝑖 ) 𝑦 𝑔𝑖 = 0, para 𝑖 = 1: 𝑛 + 1 Paso 6: Imprimir la solución aproximada 𝒖
Paso 7: Graficar la solución aproximada 𝒖

Consultas adicionales
Realice consultas adicionales a

través del chat o solicita al
docente activar el micrófono
para participar.
GRACIAS

Semana 10 - ECUACIONES HIPERBOLICAS

Cargado por

Copyright:

Formatos disponibles

Semana 10 - ECUACIONES HIPERBOLICAS

Cargado por

Información del documento

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

Semana 10 - ECUACIONES HIPERBOLICAS

Cargado por

Copyright:

Formatos disponibles

ECUACIONES HIPERBÓLICAS

Métodos Numéricos para Ingeniería

UPN, PASIÓN POR TRANSFORMAR VIDAS

A continuación estudiaremos los métodos de las diferencias finitas que

2. Derivadas Parciales de una Función

a) Respecto de x en el punto (𝑥, 𝑦) se define por:

b) Respecto de y en el punto (𝑥, 𝑦) se define por:

Una ecuación que tiene derivadas parciales de una función desconocida

Debido a la amplia aplicación a la ingeniería, centraremos nuestro

𝑨𝒖𝒙𝒙 + 𝑩𝒖𝒙𝒚 + 𝑪𝒖𝒚𝒚 + 𝑫 = 𝟎

Donde A, B y C son funciones de 𝑥 e 𝑦, y D es una función de 𝑥, 𝑦, 𝑢, 𝑢𝑥 y

3. Clasificación de una EDP

Condición Categoría Ejemplo

u(0, t ) = 0 y u(a , t ) = 0 para 0  t  b,

Discretización de la variable dependiente

Segundo, se discretiza la condición inicial 𝑡 = 𝑡1 = 0

Evaluando en los nodos de la malla, se tiene:

Reemplazando en la EDP de la onda:

Es la ecuación en diferencias finitas que se usará como aproximación

𝑢𝑖,𝑗+1 − 2𝑢𝑖,𝑗 + 𝑢𝑖,𝑗−1 = 𝑟 2 𝑢𝑖+1,𝑗 − 2𝑢𝑖,𝑗 + 𝑢𝑖−1,𝑗

𝑢𝑖,𝑗+1 = 𝑟 2 𝑢𝑖+1,𝑗 − 2𝑟 2 𝑢𝑖,𝑗 + 𝑟 2 𝑢𝑖−1,𝑗 + 2𝑢𝑖,𝑗 − 𝑢𝑖,𝑗−1

𝑢𝑖,𝑗+1 = 2 − 2𝑟 2 𝑢𝑖,𝑗 + 𝑟 2 𝑢𝑖+1,𝑗 + 𝑢𝑖−1,𝑗 − 𝑢𝑖,𝑗−1

𝒖𝒊−𝟏,𝒋 𝒖𝒊,𝒋 𝒖𝒊+𝟏,𝒋

Discretización de las condiciones de contorno:

❑ 𝑢𝑖,2 = 𝑓𝑖 + 𝑘𝑔𝑖 . Al utilizar esta fórmula el error de truncamiento es

pues el error de truncamiento es pequeño.

En resumen, la forma cómo trabaja la fórmula de diferencia finita para la

𝒖𝒊−𝟏,𝒋 𝒖𝒊,𝒋 𝒖𝒊+𝟏,𝒋

Paso 2: Calculo del 𝑛,m y r

Paso 3: Discretización del dominio

Paso 4: Discretización de las condiciones de contorno

𝑢1,𝑗 = 0 y 𝑢𝑛+1,𝑗 = 0, para 𝑗 = 1: 𝑚 + 1

𝑢𝑖,1 = 𝑓𝑖 y 𝑢𝑖,2 = 𝑓𝑖 +k𝑔𝑖 , para 𝑖 = 2: 𝑛

Paso 5: Calculo de las aproximaciones

𝒖𝒊,𝒋+𝟏 = 2 − 2𝑟 2 𝑢𝑖,𝑗 + 𝑟 2 𝑢𝑖+1,𝑗 + 𝑢𝑖−1,𝑗 − 𝑢𝑖,𝑗−1

Paso 6: Imprimir la solución aproximada 𝒖

Paso 7: Graficar la solución aproximada 𝒖

𝑢𝑡𝑡 𝑥, 𝑡 = 4𝑢𝑥𝑥 𝑥, 𝑡 para 0 < 𝑥 < 1 y 0 < 𝑡 < 0.5,

con las condiciones de contorno:

Elabore el algoritmo para solucionar numéricamente esta ecuación con su

respectiva implementación en OCTAVE:, h=0.1 y k=0.05

Paso 1: Ingrese a, b, c, h, k. Paso 4: Discretización de las condiciones de contorno

𝑢1,𝑗 = 0 y 𝑢𝑛+1,𝑗 = 0, para 𝑗 = 1: 𝑚 + 1

Paso 3: Discretización del dominio Paso 5: Calculo de las aproximaciones

Paso 7: Graficar la solución aproximada 𝒖

Realice consultas adicionales a

También podría gustarte

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.