ECUACIONES DIFERENCIALES PARA
INGENIERIA
Semana 13
CONTENIDOS
Ecuación de calor.
Ecuación de onda.
Ecuación de Laplace.
Logros esperados
• Resolver problemas relacionados a la ecuación del calor.
• Resolver problemas relacionados a la ecuación de la onda.
• Resolver problemas relacionados a la ecuación de Laplace.
Ecuación de calor
Ecuación de calor
Considere una varilla delgada de longitud 𝐿 con una temperatura inicial 𝑓 (𝑥) en
toda la varilla y cuyos extremos se mantienen a temperatura cero durante todo el
tiempo 𝑡 > 0, entonces la temperatura 𝑢(𝑥; 𝑡) en la varilla se determina del
problema con valores en la frontera:
𝜕 2 𝑢 𝜕𝑢
𝑘 2= , 0 < 𝑥 < 𝐿, 𝑡 > 0 … … . (1)
𝜕𝑥 𝜕𝑡
𝑢 0; 𝑡 = 0, 𝑢 𝐿; 𝑡 = 0, 𝑡>0
𝑢 𝑥; 0 = 𝑓 𝑥 , 0<𝑥<𝐿
𝐾
Donde 𝑘 = se le llama difusividad térmica y es una constante positiva.
𝜌𝛾
Ecuación de calor
Resuelva la ecuación de calor (1) sujeta a las condiciones dadas. Suponga una varilla
de longitud 𝐿.
𝑢 0; 𝑡 = 0, 𝑢 𝐿; 𝑡 = 0, 𝑡>0
𝐿
1, 0<𝑥<
𝑢 𝑥; 0 = 2
𝐿
0, <𝑥<𝐿
2
Solución
Ecuación de calor
Resuelva la ecuación de calor (1) sujeta a las condiciones dadas. Suponga una varilla
de longitud 𝐿.
𝑢 0; 𝑡 = 0, 𝑢 𝐿; 𝑡 = 0, 𝑡>0
𝑢 𝑥; 0 = 𝑥(𝐿 − 𝑥)
Solución
Ecuación de calor
Una barra de 100cm de longitud tiene sus extremos 𝑥 = 0, 𝑥 = 100 mantenidos a
0ºC. Inicialmente, la mitad de la barra esta a 60ºC, mientras que la otra mitad esta a
40ºC. Asumiendo un coeficiente de difusividad de 0,16 unidades c.g.s y un entorno
aislado, calcule la temperatura en cualquier posición de la barra en cualquier tiempo.
Solución
Ecuación de onda
Ecuación de onda
para una cuerda que vibra podemos especificar su desplazamiento inicial (o la
forma) 𝑓 (𝑥) así como su velocidad inicial 𝑔(𝑥). El desplazamiento vertical 𝑢(𝑥; 𝑡) de
la cuerda vibratoria de longitud 𝐿 que se muestra en la figura (a) se determina del
problema las dos condiciones iniciales:
𝜕𝑢
𝑢 𝑥; 0 = 𝑓(𝑥) y ∕ =𝑔 𝑥 , 0<𝑥<𝐿
𝜕𝑡 𝑡=0
El desplazamiento vertical 𝑢(𝑥; 𝑡) de la cuerda vibratoria de longitud 𝐿 que se
muestra en la figura (a) se determina del problema con valores en la frontera:
𝜕 2𝑢 𝜕 2𝑢
𝑎2 2 = 2 , 0 < 𝑥 < 𝐿, 𝑡>0
𝜕𝑥 𝜕𝑡
𝑢 0; 𝑡 = 0, 𝑢 𝐿; 𝑡 = 0, 𝑡>0
𝜕𝑢
𝑢 𝑥; 0 = 𝑓 𝑥 , ∕𝑡=0 = 𝑔 𝑥 , 0<𝑥<𝐿
𝜕𝑡
𝑇
Donde 𝑎2 = .
𝜌
Ecuación unidimensional de la onda
La ecuación de onda está dada por:
𝝏𝟐𝒖 𝝏𝟐𝒖
𝒂𝟐 𝟐 = 𝟐 … (𝟏)
𝝏𝒙 𝝏𝒕
𝒖 𝟎; 𝒕 = 𝟎, 𝒖 𝑳; 𝒕 = 𝟎, 𝒕 > 𝟎 … (𝟐)
𝝏𝒖 𝒙; 𝟎
𝒖 𝒙; 𝟎 = 𝒇 𝒙 , =𝒈 𝒙 𝟎 < 𝒙 < 𝑳 … (𝟑)
𝝏𝒕
donde 𝒖(𝒙, 𝒕) es la deflexión de la cuerda, la ecuación (2) son las
condiciones de frontera y la ecuación (3) son las condiciones iniciales.
𝟎 𝑳
Método de separación de variables
En el método de separación de variables o método del producto, se
determinan las soluciones de la ecuación de onda en (1) de la
forma:
𝑢 𝑥; 𝑡 = 𝐹 𝑥 𝐺(𝑡)
donde:
ሷ
𝑢𝑡𝑡 𝑥; 𝑡 = 𝐹(𝑥)𝐺(𝑡) y 𝑢𝑥𝑥 𝑥; 𝑡 = 𝐹 ′′ 𝑥 𝐺 (𝑡)
Reemplazando en la ecuación en (1), se tiene:
𝐹 𝑥 𝐺ሷ 𝑡 = 𝑎2 𝐹 ′′ 𝑥 𝐺 (𝑡)
Método de separación de variables
𝑭′′(𝒙) ሷ
𝑮(𝒕)
= 𝟐 = 𝒌 … (𝟒)
𝑭(𝒙) 𝒂 𝑮(𝒕)
Analizando las condiciones de frontera:
1) Si 𝑢 0; 𝑡 = 𝐹 0 𝐺 𝑡 = 0 → 𝑭 𝟎 = 0 … (∗)
2) Si 𝑢 𝐿; 𝑡 = 𝐹 𝐿 𝐺(𝑡) → 𝑭 𝑳 = 0 … (∗∗)
De la ecuación (4) analizamos, se tiene:
𝑭′′(𝒙)
= 𝒌 → 𝑭′′ 𝒙 − 𝒌𝑭 𝒙 = 𝟎 … (𝟓)
𝑭(𝒙)
Problema de Sturm-Lioville
Presentándose los siguientes casos en (5):
Caso 1: (𝒌 = 𝟎)
𝐹′′(𝑥) = 0 donde 𝐹 𝑥 = 𝐴𝑥 + 𝐵, usando las condiciones en (*) y (**), se
obtiene que 𝐴 = 𝐵 = 0, obteniéndose una solución trivial; 𝐹 𝑥 = 0, →
𝑢 𝑥; 𝑡 = 0
Caso 2: (𝒌 = 𝒑𝟐 , 𝒑 > 𝟎)
𝐹 ′′ 𝑥 − 𝑝2 𝐹 = 0 donde 𝐹 𝑥 = 𝐴𝑒 𝑝𝑥 + 𝐵𝑒 −𝑝𝑥 , usando las condiciones en
(*) y (**), se obtiene que 𝐴 = 𝐵 = 0, obteniéndose una solución trivial;
𝐹 𝑥 = 0, → 𝑢 𝑥; 𝑡 = 0
Problema de Sturm-Lioville
Caso 3: (𝒌 = −𝒑𝟐 , 𝒑 > 𝟎)
𝐹 ′′ 𝑥 + 𝑝2 𝐹 = 0, donde 𝐹 𝑥 = 𝐴𝑐𝑜𝑠(𝑝𝑥) + 𝐵𝑠𝑒𝑛(𝑝𝑥), usando las condiciones
en (*) y (**), se obtiene que 𝐴 = 0, 𝐵𝑠𝑒𝑛 𝑝𝐿 = 0 donde 𝐵 ≠ 0 y 𝑠𝑒𝑛 𝑝𝐿 = 0,
𝑛𝜋 𝑛𝜋
𝑝= obteniéndose una solución; 𝐹 𝑥 = 𝐵𝑛 𝑠𝑒𝑛( 𝑥)
𝐿 𝐿
De la ecuación (4) analizamos, se tiene:
ሷ
𝑮(𝒕)
= 𝒌 → 𝑮ሷ 𝒕 − 𝒌𝒂𝟐 𝑮 𝒕 = 𝟎 … (𝟔)
𝒂𝟐 𝑮(𝒕)
Problema de Sturm-Lioville
Presentándose el siguiente caso en (6):
1) 𝒌 = −𝒑𝟐 , 𝒑 > 𝟎)
𝐺ሷ 𝑡 + 𝑝2 𝑎2 𝐺(𝑡) = 0
G 𝑡 = 𝐷𝑐𝑜𝑠(𝑝𝑎𝑡) + 𝐸𝑠𝑒𝑛(𝑝𝑎𝑡),
𝑛𝜋𝑎 𝑛𝜋𝑎
𝐺 𝑡 = 𝐷𝑛 cos 𝑡 + 𝐸𝑛 𝑠𝑒𝑛 𝑡 , 𝑛 = 1; 2; 3; …
𝐿 𝐿
∞
𝑛𝜋𝑎 𝑛𝜋𝑎
𝐺 𝑡 = 𝐷𝑛 cos 𝑡 + 𝐸𝑛 𝑠𝑒𝑛 𝑡 … (7)
𝐿 𝐿
𝑛=1
𝑶𝒃𝒔𝒆𝒓𝒗𝒂𝒄𝒊ó𝒏
Para: 𝑘 = 𝑝2 𝑦 𝑘 = 0, se obtienen soluciones triviales.
Problema de Sturm-Lioville
Reemplazando en 𝑢 𝑥; 𝑡 = 𝐹 𝑥 𝐺(𝑡)
∞
𝑛𝜋𝑎 𝑛𝜋𝑎 𝑛𝜋
𝑢 𝑥; 𝑡 = 𝐷𝑛 cos 𝑡 + 𝐸𝑛 𝑠𝑒𝑛 𝑡 𝐵𝑛 𝑠𝑒𝑛( 𝑥)
𝐿 𝐿 𝐿
𝑛=1
∞
𝑛𝜋𝑎 𝑛𝜋𝑎 𝑛𝜋
𝑢 𝑥; 𝑡 = 𝐷∗ 𝑛 cos
∗
𝑡 + 𝐸 𝑛 𝑠𝑒𝑛 𝑡 𝑠𝑒𝑛( 𝑥)
𝐿 𝐿 𝐿
𝑛=1
𝑨𝒑𝒍𝒊𝒄𝒂𝒏𝒅𝒐 𝒍𝒂𝒔 𝒄𝒐𝒏𝒅𝒊𝒄𝒊𝒐𝒏𝒆𝒔 𝒊𝒏𝒊𝒄𝒊𝒂𝒍𝒆𝒔, 𝒕𝒆𝒏𝒆𝒎𝒐𝒔:
2 𝐿 𝑛𝜋 2 𝐿
𝑛𝜋
𝐷∗ 𝑛 = න 𝑓 𝑥 𝑠𝑒𝑛 𝑥 𝑑𝑥 𝐸∗𝑛 = න 𝑔 𝑥 𝑠𝑒𝑛 𝑥 𝑑𝑥
𝐿 0 𝐿 𝑛𝜋𝑎 0 𝐿
Ejemplo
Calcule la ley de vibraciones de una cuerda homogénea de longitud L, fijada en los
extremos, si en el momento inicial 𝑡 = 0 la cuerda tiene la forma de parábola
ℎ𝑥 𝐿 − 𝑥 , ℎ > 0 es constante, con velocidad inicial ausente.
Solución
Ejemplo
Resolver el problema siguiente:
2 2
2
𝜕 𝑢 𝜕 𝑢
𝑎 2
= 2, 0 < 𝑥 < 𝐿, 𝑡>0
𝜕𝑥 𝜕𝑡
𝑢 0; 𝑡 = 0, 𝑢 𝐿; 𝑡 = 0, 𝑡>0
3𝜋 𝜕𝑢
𝑢 𝑥; 0 = 𝑐𝑜𝑠 𝑥 , ∕𝑡=0 = 0, 0<𝑥<𝐿
𝐿 𝜕𝑡
Solución
Ejemplo
Resolver el problema siguiente:
2 2
2
𝜕 𝑢 𝜕 𝑢
𝑎 2
= 2, 0 < 𝑥 < 𝐿, 𝑡>0
𝜕𝑥 𝜕𝑡
𝑢 0; 𝑡 = 0, 𝑢 𝐿; 𝑡 = 0, 𝑡>0
𝜕𝑢 𝜋
𝑢 𝑥; 0 = 0, ∕𝑡=0 = 𝑐𝑜𝑠 𝑥 , 0<𝑥<𝐿
𝜕𝑡 𝐿
Solución
Ejemplo
Determine la solución de la ecuación de la onda (de longitud 𝐿 = 𝜋, extremos
fijos, y con 𝑎2 = 1) si su velocidad inicial es 0 y su deflexión inicial es 𝑓(𝑥)
donde:
1
a) 𝑓 𝑥 = 𝑠𝑒𝑛𝑥 − 𝑠𝑒𝑛(2𝑥)
2
𝜋
𝑥, 𝑠𝑖 0 < 𝑥 <
b) 𝑓 𝑥 = 𝜋 − 𝑥 3
𝜋
, 𝑠𝑖 <𝑥<𝜋
2 3
Solución
Ejemplo
Determine la solución de la ecuación de la cuerda vibrante (de longitud 𝐿 = 𝜋,
extremos fijos, y con 𝑐 2 = 1) si su velocidad inicial es 𝑔(𝑥) y su deflexión
inicial es 𝑓(𝑥) donde:
𝑎) 𝑓 𝑥 = 0; 𝑔 𝑥 = 0,1𝑠𝑒𝑛(2𝑥)
b) 𝑓 𝑥 = 0,1𝑠𝑒𝑛𝑥; 𝑔 𝑥 = −0,2𝑠𝑒𝑛𝑥
Solución
Ecuación de Laplace
Ecuación de Laplace
Suponga que deseamos encontrar la temperatura de estado estable 𝑢(𝑥, 𝑦) en
una placa rectangular cuyas aristas verticales 𝑥 = 0 y 𝑥 = 𝑎 están aislados, como se
muestra en la figura adjunta. Cuando no se escapa calor de las caras laterales de la placa,
resolvemos el siguiente problema con valores en la frontera:
Problema de Dirichlet
Un problema con valores en la frontera en el que se busca una solución de una ecuación
diferencial parcial de tipo elíptico tal como la ecuación de Laplace, ∇2 𝑢 = 0, dentro de una
región 𝑅 acotada (en el plano o en el espacio tridimensional) tal que 𝑢 tome los valores
prescritos en toda la frontera de la región se llama problema de Dirichlet. La solución del
problema de Dirichlet, para una región rectangular:
Principio de superposición
El problema de Dirichlet para un rectángulo se puede resolver con facilidad separando las
variables cuando se especifican condiciones homogéneas para dos fronteras paralelas. Sin
embargo, el método de separación de variables no se aplica a un problema de Dirichlet
cuando las condiciones en la frontera en los cuatro lados del rectángulo son no
homogéneas. Para salvar esta dificultad separamos el problema:
… … (∗)
en dos problemas, cada uno con condiciones homogéneas en la frontera, en lados
paralelos, como se muestra a continuación:
Principio de superposición
Suponga que 𝑢1 y 𝑢2 son las soluciones de los problemas 1 y 2, respectivamente. Si
definimos 𝒖 𝒙, 𝒚 = 𝒖𝟏 𝒙; 𝒚 + 𝒖𝟐 (𝒙; 𝒚), veremos que 𝑢 satisface todas las
condiciones en la frontera del problema original (*); por ejemplo,
Principio de superposición
y así sucesivamente. Además, 𝑢 es una solución de la ecuación de Laplace. En otras
palabras, al resolver los problemas 1 y 2 y sumar las soluciones, ya hemos resuelto el
problema original. Esta propiedad aditiva de las soluciones se llama principio de
superposición.
Solución 𝑢 = solución 𝑢1 del problema 1+ solución 𝑢2 del problema 2.
Donde:
Principio de superposición
La solución del problema 1 está dado por:
Principio de superposición
La solución del problema 2 está dado por:
Ejemplo
Resolver el problema siguiente problema de Laplace:
𝜕2𝑢 𝜕2𝑢
2
+ 2 = 0, 0<𝑥<𝜋 0<𝑦<𝜋
𝜕𝑥 𝜕𝑦
𝑢 0; 𝑦 = 𝑢 𝜋; 𝑦 = 𝑦,
𝑢 𝑥; 0 = 𝑢 𝑥; 𝜋 = 0
Solución
Ejemplo
Resolver el problema siguiente problema de Laplace:
𝜕2𝑢 𝜕2𝑢
2
+ 2 = 0, 0<𝑥<𝑎 0<𝑦<𝑏
𝜕𝑥 𝜕𝑦
𝑢 0; 𝑦 = 𝑢 𝑎; 𝑦 = 0,
𝑢 𝑥; 0 = 0; 𝑢 𝑥; 𝑏 = 𝑓(𝑥)
Solución
Ejemplo
Resolver el problema siguiente problema de Laplace:
𝜕2𝑢 𝜕2𝑢
2
+ 2 = 0, 0<𝑥<𝑎 0<𝑦<𝑏
𝜕𝑥 𝜕𝑦
𝑢 0; 𝑦 = 𝑢 𝑎; 𝑦 = 0,
𝑢𝑦 𝑥; 0 = 0; 𝑢 𝑥; 𝑏 = 𝑓(𝑥)
Solución
Ejemplo
Resolver el problema siguiente problema de Laplace:
𝜕2𝑢 𝜕2𝑢
2
+ 2 = 0, 0<𝑥<𝑎 0<𝑦<𝑏
𝜕𝑥 𝜕𝑦
𝑢 0; 𝑦 = 𝑢 𝑎; 𝑦 = 0,
𝑢 𝑥; 0 = 𝑓(𝑥); 𝑢 𝑥; 𝑏 = 0
Solución
Ejemplo
Resolver el problema siguiente problema de Laplace:
𝜕2𝑢 𝜕2𝑢
2
+ 2 = 0, 0<𝑥<𝑎 0<𝑦<𝑏
𝜕𝑥 𝜕𝑦
𝑢𝑥 0; 𝑦 = 𝑢𝑥 𝑎; 𝑦 = 0,
𝑢 𝑥; 0 = 𝑥; 𝑢 𝑥; 𝑏 = 0
Solución
Ejemplo
Resolver el problema siguiente problema de Laplace:
𝜕2𝑢 𝜕2𝑢
2
+ 2 = 0, 0<𝑥<1 0<𝑦<1
𝜕𝑥 𝜕𝑦
𝑢 0; 𝑦 = 0; 𝑢 1; 𝑦 = 1 − 𝑦
𝑢𝑦 𝑥; 0 = 𝑢𝑦 𝑥; 1 = 0,
Solución
Conclusiones
• Se resolvió problemas relacionados a la ecuación del calor.
• se resolvió problemas relacionados a la ecuación de la onda.
• Se resolvió problemas relacionados a la ecuación de Laplace.
REFERENCIAS
1] RE (2010). Ecuaciones diferenciales ordinarias : introducción / Ernesto Javier Espinosa
Herrera, comp (xiv, 478 p.). México, D.F. : Reverté.
[2] Cengel Ya, Palm III WJ (2014). Ecuaciones Diferenciales para Ingeniería y Ciencias (1st).
México DF: McGraw Hill Education.
[3] Zill DG, Wrights WS (2015). Ecuaciones Diferenciales con problemas con valores en la
frontera (8va. ed.). Cengage Learning: Cengage Learning.
[4] Xie WCh (2010). Differential equations for Engineers (3rd. ed.). Cambridge University
Press: Cambridge University Press.