TP10 - Ecuación - Unidimensional de Calor

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 11

TRABAJO PRÁCTICO

N°10
ECUACIONES DIFERENCIALES EN DERIVADAS PARCIALES
AÑO 2020
Nociones generales
En el trabajo práctico 10, hallaremos soluciones particulares de algunas de las
EDP más importantes :
 (1) Ecuación de calor (flujo de calor unidimensional):
𝜕 2 𝑢 𝜕𝑢
𝑘 2= 𝑘>0
𝜕𝑥 𝜕𝑡
 (2) Ecuación de onda (por ejemplo: cuerda de guitarra):
𝜕 2𝑢 𝜕 2𝑢
𝑎2 2 = 2
𝜕𝑥 𝜕𝑡
 (3) Ecuación de Laplace (se presenta en problemas independientes del tiempo,
por ejemplo: potenciales electroestáticos):
2
𝜕2𝑢 𝜕2𝑢
∇ 𝑢=0→ 2+ 2=0
𝜕𝑥 𝜕𝑦
ECUACIÓN UNIDIMENSIONAL DE CALOR
Ecuación diferencial gobernante :
𝜕 2 𝑢 𝜕𝑢
𝑘 2= 𝑘>0
𝜕𝑥 𝜕𝑡

𝑢(𝑥, 𝑡) representa la temperatura en un punto 𝑥 a lo largo de una varilla delgada,


de longitud L, en algún instante 𝑡 suponemos que el flujo de calor se da en una
sola dirección.

Condiciones iniciales:
Si 𝑓(𝑥) denota la temperatura inicial en la varilla entonces la solución 𝑢(𝑥, 𝑡)
debe satisfacer 𝑢 𝑥, 0 = 𝑓(𝑥)
Condiciones de frontera (o de borde):
Para el caso de la varilla, la condición representativa en el extremo derecho puede
ser:
 La temperatura se mantiene constante en 𝑥 = 𝐿, esto es
𝑢 𝐿, 𝑡 = 𝑢0 .

 La frontera está aislada (no entra ni sale calor):


𝜕𝑢
ቤ =0
𝜕𝑥 𝑥=𝐿
Ejemplo
Considere una varilla delgada de longitud L, con temperatura inicial 𝑓(𝑥) y cuyos
extremos se mantienen a temperatura 0° para t > 0.

Planteo del problema:


Ecuación diferencial gobernante :
𝜕 2 𝑢 𝜕𝑢
𝑘 2= 𝑘>0
𝜕𝑥 𝜕𝑡
𝑡 > 0, 0<𝑥<𝐿
Condición inicial:
𝑢 𝑥, 0 = 𝑓(𝑥)
Condiciones de frontera (o de borde):
𝑢 0, 𝑡 = 0° , 𝑢 𝐿, 𝑡 = 0°
Resolución:
Queremos hallar 𝑢 𝑥, 𝑡 bajo las condiciones planteadas. Usamos el método de
separación de variables (Ver la presentación sobre este tema).
Proponemos la solución 𝑢 𝑥, 𝑡 = 𝑋 𝑥 𝑇(𝑡) y reemplazamos en la ecuación
diferencial gobernante:

𝑘𝑋 ′′ 𝑥 𝑇 𝑡 = 𝑋 𝑥 𝑇′(𝑡)

𝑋 ′′ 𝑥 𝑇′(𝑡)
= =𝛼
𝑋 𝑥 𝑘𝑇(𝑡)

𝑋´´ − 𝛼𝑋 = 0
ቊ ′
𝑌 − 𝛼𝑘𝑌 = 0
 Si 𝛼 = 0, entonces:
𝑋´´ = 0 → 𝑋 𝑥 = 𝐶1 𝑥 + 𝐶2

𝑇 ′ = 0 → 𝑇 𝑡 = 𝐶3
Si 𝑋 0 = 𝑋 𝐿 = 0 entonces 𝐶1 = 𝐶2 = 0, y 𝑢 𝑥, 𝑡 = 0

 Si 𝛼 = 𝜆2 > 0 (Ver en la presentación de separación de variables)

𝑋(𝑥) = 𝐶4 𝑐𝑜𝑠ℎ𝜆𝑥 + 𝐶5 𝑠𝑒𝑛ℎ𝜆𝑥


൝ 2
𝑇 = 𝐶6 𝑒 𝜆 𝑘𝑡

Si 𝑋 0 = 𝑋 𝐿 = 0 entonces 𝐶4 = 𝐶5 = 0, y 𝑢 𝑥, 𝑡 = 0
Analizando los valores 𝛼 = 0 𝑦 𝛼 > 0 de la constante junto con las condiciones
de borde, se observa que en ambos casos se obtiene una solución nula.

Si 𝛼 = −𝜆2 :
𝑋(𝑥) = 𝐴𝑐𝑜𝑠𝜆𝑥 + 𝐵𝑠𝑒𝑛𝜆𝑥
ቊ 2
𝑇 = 𝐶𝑒 −𝜆 𝑘𝑡
Analizando las condiciones de borde 𝑋 0 = 𝑋 𝐿 = 0, se obtiene:
𝑋 0 = 𝐴𝑐𝑜𝑠𝜆0 + 𝐵𝑠𝑒𝑛𝜆0 = 𝐴 = 0
𝑛𝜋
𝑋 𝐿 = 𝐵𝑠𝑒𝑛𝜆𝐿 = 0 → 𝐵 ≠ 0, 𝑒𝑛𝑡𝑜𝑛𝑐𝑒𝑠 𝜆 =
𝐿
En este caso se obtienen soluciones no triviales.
2 𝑘𝑡 2 𝑘𝑡
Luego 𝑢 𝑥, 𝑡 = 𝐵𝑠𝑒𝑛𝜆𝑥 𝐶𝑒 −𝜆 = 𝐷𝑠𝑒𝑛𝜆𝑥 𝑒 −𝜆 ,
𝑛𝜋
Pero 𝜆 = entonces se tiene una solución para cada 𝑛 ∈ Ν,
𝐿
2
𝑢𝑛 𝑥, 𝑡 = 𝐷𝑛 𝑠𝑒𝑛𝜆𝑛 𝑥 𝑒 −𝜆𝑛 𝑘𝑡

Por el principio de superposición, la solución es:

∞ ∞
2
𝑢 𝑥, 𝑡 = ෍ 𝑢𝑛 𝑥, 𝑡 = ෍ 𝐷𝑛 𝑠𝑒𝑛𝜆𝑛 𝑥 𝑒 −𝜆𝑛 𝑘𝑡 (∗)
𝑛=1 𝑛=1

Nos falta obtener la constante 𝐷𝑛 , usemos la condición inicial:



2
𝑢 𝑥, 0 = 𝑓 𝑥 → ෍ 𝐷𝑛 𝑠𝑒𝑛𝜆𝑛 𝑥 𝑒 −𝜆𝑛 𝑘0 =𝑓 𝑥
𝑛=1

𝑛𝜋
෍ 𝐷𝑛 𝑠𝑒𝑛 𝑥 = 𝑓(𝑥)
𝐿
𝑛=1

Observamos que el lado izquierdo es la representación en una serie de senos de


Fourier para la función 𝑓 𝑥 en el semintervalo 0 < 𝑥 < 𝐿.
Entonces debe suceder que:
2 𝐿 𝑛𝜋
𝐷𝑛 = 𝑏𝑛 = න 𝑓 𝑥 𝑠𝑒𝑛 𝑥 𝑑𝑥
𝐿 0 𝐿
Este es un valor numérico, como la función 𝑓 𝑥 está definida en forma genérica,
no hay ninguna cuenta más para resolver.
Por lo tanto, volviendo a (*)

2 𝐿 𝑛𝜋 𝑛𝜋 −
𝑛𝜋 2
𝑢 𝑥, 𝑡 = ෍ න 𝑓 𝑥 𝑠𝑒𝑛 𝑥 𝑑𝑥 𝑠𝑒𝑛 𝑥𝑒 𝐿 𝑘𝑡
𝐿 0 𝐿 𝐿
𝑛=1

2 𝐿 𝑛𝜋 𝑛𝜋 −
𝑛𝜋 2
𝑢 𝑥, 𝑡 = ෍ න 𝑓 𝑥 𝑠𝑒𝑛 𝑥 𝑑𝑥 𝑠𝑒𝑛 𝑥𝑒 𝐿 𝑘𝑡
𝐿 0 𝐿 𝐿
𝑛=1

Esta es la solución que estábamos buscando.

La solución debe depender de las variables 𝑥 𝑦 𝑡.

También podría gustarte