0% encontró este documento útil (0 votos)
107 vistas18 páginas

Ecuaciones Diferenciales Parciales: Fundamentos

El documento describe las ecuaciones diferenciales parciales de segundo orden lineales y los métodos para resolverlas. Explica que estas ecuaciones se clasifican como hiperbólicas, parabólicas o elípticas dependiendo de los coeficientes. También presenta el método de separación de variables para resolver este tipo de ecuaciones, dividiéndolas en ecuaciones ordinarias separadas en x e y. Finalmente, aplica estos conceptos para modelar el flujo de calor a través de un alambre delgado usando la ecuación del calor
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
0% encontró este documento útil (0 votos)
107 vistas18 páginas

Ecuaciones Diferenciales Parciales: Fundamentos

El documento describe las ecuaciones diferenciales parciales de segundo orden lineales y los métodos para resolverlas. Explica que estas ecuaciones se clasifican como hiperbólicas, parabólicas o elípticas dependiendo de los coeficientes. También presenta el método de separación de variables para resolver este tipo de ecuaciones, dividiéndolas en ecuaciones ordinarias separadas en x e y. Finalmente, aplica estos conceptos para modelar el flujo de calor a través de un alambre delgado usando la ecuación del calor
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

Ecuaciones entre Derivadas Parciales (EDP)

Ecuación diferencial parcial lineal de segundo orden


Es una ecuación de la forma:

𝜕2𝑢 𝜕2𝑢 𝜕2𝑢 𝜕𝑢 𝜕𝑢


𝐴 2 +𝐵 +𝐶 2+𝐷 +𝐸 + 𝐹𝑢 = 𝐺 (1)
𝜕𝑥 𝜕𝑥𝜕𝑦 𝜕𝑦 𝜕𝑥 𝜕𝑦

Donde los coeficiente: 𝐴, 𝐵, 𝐶, 𝐷, 𝐸 𝑦 𝐹 son constantes o funciones de 𝑥 y de 𝑦.

Ecuación diferencial parcial lineal de segundo orden homogénea

Es cuando 𝐺 𝑥, 𝑦 = 0 en la ecuación (1).


𝜕2𝑢 𝜕2𝑢
+ =0
𝜕𝑥 2 𝜕𝑦 2

Ecuación diferencial parcial lineal de segundo orden no homogénea

Es cuando 𝐺 𝑥, 𝑦 ≠ 0 en la ecuación (1).


𝜕2𝑢 𝜕2𝑢
+ = 2𝑥𝑦
𝜕𝑥 2 𝜕𝑦 2
Clasificación de las ecuaciones
La ecuación diferencial parcial lineal de segundo orden
𝜕2𝑢 𝜕2𝑢 𝜕2𝑢 𝜕𝑢 𝜕𝑢
𝐴 2+𝐵 +𝐶 2+𝐷 +𝐸 + 𝐹𝑢 = 𝐺 (1)
𝜕𝑥 𝜕𝑥𝜕𝑦 𝜕𝑦 𝜕𝑥 𝜕𝑦

Donde los coeficiente: 𝐴, 𝐵, 𝐶, 𝐷, 𝐸 𝑦 𝐹 son constantes reales se clasifican en:

Hiperbólicas si: 𝐵2 − 4𝐴𝐶 > 0

Parabólica si: 𝐵2 − 4𝐴𝐶 = 0

Elípticas si : 𝐵2 − 4𝐴𝐶 < 0


Método de Separación de variables
𝜕2𝑢 𝜕2𝑢 𝑋′′ 𝑇′′
= =
𝜕𝑥 2 𝜕𝑡 2 𝑋 𝑇
Como el miembro izquierdo de la última ecuación es
𝑢 𝑥, 𝑡 = 𝑋 𝑥 𝑇(𝑡)
independiente de 𝑡 y es igual al miembro derecho, que
es independiente de 𝑥 , concluimos que ambos
𝜕𝑢 miembros de la ecuación son independientes de 𝑥 y de
= 𝑋′𝑇 t . En otras palabras, cada miembro de la ecuación
𝜕𝑥
debe ser constante. Para fines prácticos, resulta
𝜕2𝑢 conveniente escribir esta verdadera constante
= 𝑋′′𝑇 de separación como 𝜆:
𝜕𝑥 2
𝑋′′ 𝑇′′
𝜕𝑢 = = −𝜆
= 𝑋𝑇′′ 𝑋 𝑇
𝜕𝑡
𝑋 ′′ + 𝜆𝑋 = 0 𝑇 ′′ + 𝜆𝑇 = 0
𝜕2𝑢
= 𝑋𝑇′′
𝜕𝑡 2 Estas dos ecuaciones son dos ecuaciones ordinarias (EDO):
Se tendrán tres casos para 𝜆:

1. 𝜆 = 0
𝑋′′𝑇 = 𝑋𝑇′′
2. 𝜆 = −𝛼 2 < 0
3. 𝜆 = 𝛼 2 > 0
Caso I:Para 𝜆 = 0 𝑢 𝑥, 𝑡 = 𝑋 𝑥 𝑇(𝑡)
𝑇 ′′ + 𝜆𝑇 = 0
𝑋 ′′ + 𝜆𝑋 = 0
La primera solución para 𝜆 = 0 será:
𝑇 ′′ = 0
𝑋 ′′ = 0
Integrando:
Integrando: 𝑢 𝑥, 𝑡 = 𝑐1 𝑥 + 𝑐2 𝑐3 𝑡 + 𝑐4

𝑇′′𝑑𝑡 = 0𝑑𝑡
𝑋′′𝑑𝑥 = 0𝑑𝑥

𝑇 ′ = 𝑐1
𝑋′ = 𝑐1

𝑇 ′ 𝑑𝑡 = 𝑐1 𝑑𝑡
𝑋 ′ 𝑑𝑥 = 𝑐1 𝑑𝑥

𝑇 = 𝑐3 𝑡 + 𝑐4
𝑋 = 𝑐1 𝑥 + 𝑐2
Caso II:Para 𝜆 = −𝛼 2 Caso III:Para 𝜆 = 𝛼 2
𝑋 ′′ −𝛼 2 𝑋 = 0 𝑋 ′′ +𝛼 2 𝑋 = 0
𝑋 ′′ −𝛼 2 𝑋 = 0 𝑋 ′′ +𝛼 2 𝑋 = 0
𝐷 2 −𝛼 2 𝑋 = 0 𝐷2 +𝛼 2 𝑋 = 0
𝑟 2 −𝛼 2 = 0 𝑟 2 +𝛼 2 = 0

𝑟 = ±𝛼 𝑟 = ±𝛼𝑖

𝑋 = 𝑐5 𝑒 𝛼𝑥 + 𝑐6 𝑒 −𝛼𝑥 𝑋 = 𝑒 𝑎𝑖𝑥 𝑐9 𝑐𝑜𝑠𝑏𝑥 + 𝑐10 𝑠𝑒𝑛𝑏𝑥


Lo mismo para T:
𝑇 = 𝑐7 𝑒 𝛼𝑡 + 𝑐8 𝑒 −𝛼𝑡 𝑋 = 𝑒 0𝑥 𝑐9 𝑐𝑜𝑠𝛼𝑥 + 𝑐10 𝑠𝑒𝑛𝛼𝑥

𝑋 = 𝑐9 𝑐𝑜𝑠𝛼𝑥 + 𝑐10 𝑠𝑒𝑛𝛼𝑥


Lo mismo para T:
𝑢 𝑥, 𝑡 = 𝑐5 𝑒 𝛼𝑥 + 𝑐6 𝑒 −𝛼𝑥 𝑐7 𝑒 𝛼𝑡 + 𝑐8 𝑒 −𝛼𝑡
𝑇 = 𝑐11 𝑐𝑜𝑠𝛼𝑥 + 𝑐12 𝑠𝑒𝑛𝛼𝑥

𝑢 𝑥, 𝑡 = 𝑐9 𝑐𝑜𝑠𝛼𝑥 + 𝑐10 𝑠𝑒𝑛𝛼𝑥 𝑐11 𝑐𝑜𝑠𝛼𝑡 + 𝑐12 𝑠𝑒𝑛𝛼t


Principio de Superposición

𝑆𝑖 𝑢1 , 𝑢2 , … , 𝑢𝑘 son las soluciones de una ecuación diferencial lineal homogénea, entonces


la combinación lineal:

𝑢 = 𝑐1 𝑢1 + 𝑐2 𝑢2 + ⋯ + 𝑐𝑘 𝑢𝑘

𝐷𝑜𝑛𝑑𝑒 𝑙𝑎𝑠 𝑐𝑖 , 𝑖 = 1,2, … , 𝑘 𝑠𝑜𝑛 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒𝑠, 𝑒𝑠 𝑡𝑎𝑚𝑏𝑖é𝑛 𝑢𝑛𝑎 𝑠𝑜𝑙𝑢𝑐𝑖ó𝑛. :


La ecuación del Calor
Desarollaremos un modelo para flujo de calor a través de un alambre delgado, aislado, cuyos
extremos se mantienen a la temperatura constante 0℃ y cuya distribución inicial de
temperatura hay que especificar.
Sea L la longitud tiene área transversal A y coincide con el eje x en el intervalo [0, L].

Sección transversal de area A

∆𝑥

0 𝑥 𝑥 + ∆𝑥 𝐿
Suponiendo que:
 Dentro de la varilla, el flujo de calor tiene lugar sólo en la dirección x.

 La superficie lateral, o curva, de la varilla está aislada; esto es, no escapa calor de
su superficie

 No se está generando calor dentro de la varilla.

 La varilla es homogénea; esto es, su masa por unidad de volumen es constante.

 El calor específico y la conductividad térmica K del material de la varilla son


constantes.
Leyes y principios::
1.Calor: es energía en transferencia que solo es identificable cuando cruza la frontera
de un sistema.

2.Conducción de calor: La tasa del flujo de calor (la cantidad de calor por unidad
de tiempo que fluye a través de una unidad de área transversal en A) es proporcional a
𝜕𝑢
𝜕𝑥, el gradiente de temperatura en A. La constante de proporcionalidad k se
llama la conductividad térmica del material. En general, la conductividad térmica
puede variar de un punto a otro: k = k(x). En nuestro modelo es constante.
𝜕𝑄 𝜕𝑢
= −𝑘𝐴
𝜕𝑡 𝜕𝑥

3. Dirección del flujo de calor: La dirección del flujo de calor siempre va de puntos
de temperatura más alta a puntos de temperatura más baja.

4. Capacidad calórica específica: La cantidad de calor necesaria para elevar


la temperatura de un objeto de masa m en una cantidad ∆𝑢 es 𝑐𝑚∆𝑢, donde la
constante c es la capacidad calórica específica del material. La capacidad
calórica específica, al igual que la conductividad térmica, puede variar con la posición
𝑐 = 𝑐(𝑥). En nuestro modelo es constante.
∆𝑄𝑡 = 𝑐𝑚∆𝑢
Sección transversal de area A

∆𝑥

0 𝑥 𝑥 + ∆𝑥 𝐿
Sea 𝑄 el flujo de calor o ganancia de calor que fluye de 𝑥 a 𝑥 + ∆𝑥 a travez de
la seccion 𝐴 durante un intervalo de tiempo ∆𝑡, entonces la formula para la
conduccion de calor será:
𝜕𝑢
𝑄(𝑥) = −𝑘𝐴∆𝑡
𝜕𝑥
𝜕𝑢
𝑄 = −𝑘𝐴∆𝑡
𝜕𝑥
∆𝑄 = 𝐸𝑛𝑡𝑟𝑎𝑑𝑎 𝑑𝑒 𝑐𝑎𝑙𝑜𝑟 𝑒𝑛 𝑥 − 𝑆𝑎𝑙𝑖𝑑𝑎𝑑 𝑑𝑒 𝑐𝑎𝑙𝑜𝑟 𝑒𝑛 𝑥 + ∆𝑥

𝜕𝑢(𝑥, 𝑡) 𝜕𝑢 𝑥 + ∆𝑥, 𝑡
∆𝑄 = −𝑘𝐴∆𝑡 − −𝑘𝐴∆𝑡
𝜕𝑥 𝜕𝑥

𝜕𝑢 𝑥 + ∆𝑥, 𝑡 𝜕𝑢(𝑥, 𝑡)
∆𝑄 = 𝑘𝐴∆𝑡 − (1)
𝜕𝑥 𝜕𝑥
Por el principio de capacidad calorífica, la ganancia de calor debe ser igual:
∆𝑄 = 𝑐𝑚∆𝑢(𝑥, 𝑡)
∆𝑄 = 𝑐𝑚 𝑢 𝑥, 𝑡 + ∆𝑡 − 𝑢(𝑥, 𝑡) (2)

Igualando las ecuaciones (1) y (2):

𝜕𝑢 𝑥 + ∆𝑥, 𝑡 𝜕𝑢(𝑥, 𝑡)
𝑘𝐴∆𝑡 − = 𝑐𝑚 𝑢 𝑥, 𝑡 + ∆𝑡 − 𝑢(𝑥, 𝑡) (3)
𝜕𝑥 𝜕𝑥

Para deteminar la masa de la seccion ∆𝑥, tenemos que el volumen de un prisma es igual
al area de la base por la altura.
𝑉 = 𝐴∆𝑥
𝑉 = 𝐴𝑏 ℎ
𝑃𝑒𝑟𝑜:
𝐴𝑏 = 𝐴
𝑚
𝜌=
𝑉
ℎ = ∆𝑥 𝑚 = 𝜌𝑉

𝑚 = 𝜌𝐴∆𝑥
𝑆𝑢𝑠𝑡. 𝑒𝑛 𝑙𝑎 𝑒𝑐. 3 :
𝜕𝑢 𝑥 + ∆𝑥, 𝑡 𝜕𝑢(𝑥, 𝑡)
𝑘𝐴∆𝑡 − = 𝑐𝜌𝐴∆𝑥 𝑢 𝑥, 𝑡 + ∆𝑡 − 𝑢(𝑥, 𝑡)
𝜕𝑥 𝜕𝑥
𝜕𝑢 𝑥 + ∆𝑥, 𝑡 𝜕𝑢(𝑥, 𝑡)
𝑘𝐴∆𝑡 − = 𝑐𝜌𝐴∆𝑥 𝑢 𝑥, 𝑡 + ∆𝑡 − 𝑢(𝑥, 𝑡) ÷𝐴
𝜕𝑥 𝜕𝑥

𝜕𝑢 𝑥 + ∆𝑥, 𝑡 𝜕𝑢(𝑥, 𝑡)
𝑘 − 𝑢 𝑥, 𝑡 + ∆𝑡 − 𝑢(𝑥, 𝑡)
𝜕𝑥 𝜕𝑥
=
𝑐𝜌 ∆𝑥 ∆𝑡

Tomando los limites cuando ∆𝑥 → 0 𝑦 ∆𝑡 → 0

𝜕𝑢 𝑥 + ∆𝑥, 𝑡 𝜕𝑢(𝑥, 𝑡)
𝑘 − 𝑢 𝑥, 𝑡 + ∆𝑡 − 𝑢(𝑥, 𝑡)
𝜕𝑥 𝜕𝑥
lím = lím
𝑐𝜌 ∆𝑥→0 ∆𝑥 ∆𝑡→0 ∆𝑡

𝑘 𝜕 2 𝑢(𝑥, 𝑡) 𝜕𝑢(𝑥, 𝑡)
=
𝑐𝜌 𝜕𝑥 2 𝜕𝑡

La difusividad térmica se define como 𝛽 = 𝑘 𝑐𝜌 , por tanto la ecuación se puede


expresar como:
𝜕 2 𝑢 𝜕𝑢
𝛽 2= Ecuación del flujo de calor en una dimensión.
𝜕𝑥 𝜕𝑡
Solución del problema de valores en la frontera de la ecuación del calor

𝜕 2 𝑢 𝜕𝑢
𝛽 2= 0 < 𝑥 < 𝐿, 𝑡>0
𝜕𝑥 𝜕𝑡
𝑉𝑎𝑙𝑜𝑟𝑒𝑠 𝑒𝑛 𝑙𝑎 𝑓𝑟𝑜𝑛𝑡𝑒𝑟𝑎:
𝑢 0, 𝑡 = 0 𝑢 𝐿, 𝑡 = 0
𝑢 0, 𝑡 = 0, 𝑢 𝐿, 𝑡 = 0, 𝑡>0

𝑥
𝑢 𝑥, 0 = 𝑓 𝑥 , 0<𝑥<𝐿 𝐿
0
Utilizando el método de separación de variables
𝑢 𝑥, 𝑡 = 𝑋 𝑥 𝑇(𝑡)
𝛽𝑋′′𝑇 = 𝑋𝑇′
𝜕𝑢
= 𝑋′𝑇 𝑋′′ 𝑇′
𝜕𝑥 = = −𝜆
𝑋 𝛽𝑇
𝜕2𝑢
= 𝑋′′𝑇
𝜕𝑥 2
𝑋 ′′ + 𝜆𝑋 = 0
𝜕𝑢
= 𝑋𝑇′ 𝑇′ + 𝜆𝛽𝑇 = 0
𝜕𝑡
Evaluando las condiciones en las fronteras: Evaluando los condiciones en la frontera:
𝑢 𝑥, 𝑡 = 𝑋 𝑥 𝑇(𝑡)
𝑋(𝑥) = 𝑐1 𝑥 + 𝑐2
𝑢 0, 𝑡 = 0, 𝑢 𝐿, 𝑡 = 0, 𝑡>0
𝑋(0) = 𝑐2
𝑢 0, 𝑡 = 𝑋 0 𝑇 𝑡 = 0

𝑢 𝐿, 𝑡 = 𝑋 𝐿 𝑇 𝑡 = 0 0 = 𝑐2 𝑐2 = 0
𝐷𝑒 𝑚𝑜𝑑𝑜 𝑞𝑢𝑒:
Puesto que estas igualdades deben ser 𝑋(𝑥) = 𝑐1 𝑥
válidas para todo tiempo 𝑡, entonces:
𝑋 0 =0 𝑦 𝑋 𝐿 =0 𝑋(𝐿) = 𝑐1 𝐿

𝑋 ′′ + 𝜆𝑋 = 0
0 = 𝑐1 𝐿 𝑐1 = 0
Caso I:Para 𝜆 = 0
𝑋 ′′ + 𝜆𝑋 = 0 𝑐1 = 𝑐2 = 0, 𝑝𝑜𝑟 𝑡𝑎𝑛𝑡𝑜:

𝑋 ′′ = 0 𝑋(𝑥) = 0
𝐿𝑎 𝑠𝑜𝑙𝑢𝑐𝑖ó𝑛 𝑠𝑒𝑟á:
Integrando:
𝑢 𝑥, 𝑡 = 𝑋 𝑥 𝑇(𝑡)
𝑋 ′ = 𝑐1
𝑢 𝑥, 𝑡 = 0 𝑆𝑜𝑙𝑢𝑐𝑖ó𝑛 𝑡𝑟𝑖𝑣𝑖𝑎𝑙 𝑝𝑎𝑟𝑎 𝜆 = 0
𝑋 = 𝑐1 𝑥 + 𝑐2
Caso II:Para 𝜆 = −𝛼 2 < 0 𝑐5 + 𝑐6 = 0
𝑋 ′′ −𝛼 2 𝑋 = 0 𝑐5 𝑒 𝛼𝐿 + 𝑐6 𝑒 −𝛼𝐿 = 0

𝑋 ′′ −𝛼 2 𝑋 = 0
1 1
𝐷 2 −𝛼 2 𝑋 = 0 ∆= = 𝑒 −𝛼𝐿 −𝑒 𝛼𝐿
𝑒 𝛼𝐿 𝑒 −𝛼𝐿
2 2
𝑟 −𝛼 = 0
Como el sistema de ecuaciones es homogéneo y
𝑟 = ±𝛼 ∆≠ 0, entonces:
𝑐5 = 𝑐6 = 0
𝑋(𝑥) = 𝑐5 𝑒 𝛼𝑥 + 𝑐6 𝑒 −𝛼𝑥
𝑋(𝑥) = 0
Evaluando los condiciones en la frontera:
𝐿𝑎 𝑠𝑜𝑙𝑢𝑐𝑖ó𝑛 𝑠𝑒𝑟á:
𝑋(0) = 𝑐5 𝑒 0 + 𝑐6 𝑒 0 𝑢 𝑥, 𝑡 = 𝑋 𝑥 𝑇(𝑡)
𝑢 𝑥, 𝑡 = 0
𝑐5 + 𝑐6 = 0
𝑆𝑜𝑙𝑢𝑐𝑖ó𝑛 𝑡𝑟𝑖𝑣𝑖𝑎𝑙 𝑝𝑎𝑟𝑎 𝜆 = −𝛼 2 < 0
𝑋(𝐿) = 𝑐5 𝑒 𝛼𝐿 + 𝑐6 𝑒 −𝛼𝐿

𝑐5 𝑒 𝛼𝐿 + 𝑐6 𝑒 −𝛼𝐿 = 0
Caso III:Para 𝜆 = 𝛼 2 > 0 𝑋 𝑥 = 𝑐2 𝑠𝑒𝑛𝛼𝑥
𝑋 ′′ +𝛼 2 𝑋 = 0 𝑋 𝐿 = 𝑐2 𝑠𝑒𝑛𝛼𝐿
𝑋 ′′ +𝛼 2 𝑋 = 0 0 = 𝑐2 𝑠𝑒𝑛𝛼𝐿
𝐷 2 +𝛼 2 𝑋 = 0 𝑐2 𝑠𝑒𝑛𝛼𝐿 = 0 (2)
𝑟 2 +𝛼 2 = 0 𝑐2 = 0 ó 𝑠𝑒𝑛𝛼𝐿 = 0
𝑟 = ±𝛼𝑖 𝑐2 = 0 nos conduce a otra solución trivial
𝑢 𝑥, 𝑡 = 0
𝑋 = 𝑒 𝑎𝑖𝑥 𝑐1 𝑐𝑜𝑠𝑏𝑥 + 𝑐2 𝑠𝑒𝑛𝑏𝑥
Pero la ecuación (2) se verifica para:

𝑋 = 𝑒 0𝑥 𝑐1 𝑐𝑜𝑠𝛼𝑥 + 𝑐2 𝑠𝑒𝑛𝛼𝑥 𝑐2 ≠ 0 𝑐𝑢𝑎𝑛𝑑𝑜 𝑠𝑒𝑛𝛼𝐿 = 0

𝑋 𝑥 = 𝑐1 𝑐𝑜𝑠𝛼𝑥 + 𝑐2 𝑠𝑒𝑛𝛼𝑥 (1) 𝛼𝐿 = 𝑎𝑟𝑐𝑠𝑒𝑛(0)


𝛼𝐿 = 𝑛𝜋
Evaluando los condiciones en la frontera: 𝑛𝜋
𝛼= 𝑛 = 1,2,3, …
0 𝐿
𝑋(0) = 𝑐1 𝑐𝑜𝑠0 + 𝑐2 𝑠𝑒𝑛0 𝑃𝑜𝑟 𝑡𝑎𝑛𝑡𝑜 𝑙𝑎 𝑒𝑐. 𝑋 ′′ + 𝜆𝑋 = 0 tiene soluciones
No triviales para:
0 = 𝑐1 𝑐1 = 0
2 2
𝜆𝑛 = 𝛼 𝑛 = 𝑛 𝜋
2
𝑛 = 1,2,3, …
Sust.: en (1): 𝐿 2
Los valores 𝜆𝑛 y las soluciones correspondientes:
𝑛𝜋𝑥 𝑛 = 1,2,3, …
𝑋 𝑥 = 𝑐2 𝑠𝑒𝑛 ,
𝐿

son los valores propios y las funciones propias, respectivamente, del problema planteado
en 𝑋 ′′ + 𝜆𝑋 = 0 .
Pero:
Resolviendo:
𝑢 𝑥, 𝑡 = 𝑋 𝑥 𝑇(𝑡)
𝑇′ + 𝜆𝛽𝑇 = 0
𝐷 + 𝜆𝛽 𝑇 = 0 𝑛 2 𝜋2 𝑛𝜋𝑥
−𝛽 𝑡
𝑢𝑛 𝑥, 𝑡 = 𝑐2 𝑐3 𝑒 𝐿2 𝑠𝑒𝑛
𝑟 + 𝜆𝛽 = 0 𝐿

𝑛2 𝜋2 𝑛𝜋𝑥
−𝛽 𝑡
𝑟 = −𝜆𝛽 𝑢𝑛 𝑥, 𝑡 = 𝐴𝑛 𝑒 𝐿2 𝑠𝑒𝑛
𝐿
−𝜆𝛽𝑡
𝑇 𝑡 = 𝑐3 𝑒
Por el principio de superposición, la solución
Pero: sera :
𝑛2 𝜋 2 ∞
𝜆= 2 𝑛 2 𝜋2 𝑛𝜋𝑥
−𝛽 𝑡
𝐿 𝑢 𝑥, 𝑡 = 𝐴𝑛 𝑒 𝐿2 𝑠𝑒𝑛 3
𝑛2 𝜋2
𝐿
−𝛽 𝑡 𝑛=1
𝑇 𝑡 = 𝑐3 𝑒 𝐿2
Para la distribución inicial de temperatura 𝑡 = 0 se tendrá:

𝑛 2 𝜋2 𝑛𝜋𝑥
−𝛽 0
𝑢 𝑥, 0 = 𝐴𝑛 𝑒 𝐿2 𝑠𝑒𝑛
𝐿
𝑛=1


𝑛𝜋𝑥
𝑢 𝑥, 0 = 𝑓(𝑥) = 𝐴𝑛 𝑠𝑒𝑛
𝐿
𝑛=1

Esta solución es la expansión de medio intervalo de una serio de Fourier en seno.

Comparando, 𝐴𝑛 = 𝑏𝑛 , 𝑛 = 1,2,3, … , 𝑑𝑒 𝑚𝑜𝑑𝑜 𝑞𝑢𝑒:

𝐿
2 𝑛𝜋𝑥
𝐴𝑛 = 𝑓(𝑥)𝑠𝑒𝑛 𝑑𝑥
𝐿 0 𝐿

Sust. 𝐴𝑛 en la ec. (3) ,tendremos la solución al problema de valores en la frontera,


como la serie infinita:


2 𝐿
𝑛𝜋𝑥 𝑛 2 𝜋2 𝑛𝜋𝑥
−𝛽 𝑡
𝑢 𝑥, 𝑡 = 𝑓(𝑥)𝑠𝑒𝑛 𝑑𝑥 𝑒 𝐿2 𝑠𝑒𝑛
𝐿 0 𝐿 𝐿
𝑛=1

También podría gustarte