Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 60 vistas

    Ejercicio 4 Edp Calor

    Cargado por

    Bryan Gustavo Paredes
    Este documento presenta la resolución de un problema de condiciones de contorno no homogéneas mediante el método de separación de variables. Primero se transforma el problema a uno con condiciones homogéneas mediante una función. Luego se determinan las funciones propias y valores propios del problema transformado. Finalmente, se obtiene la solución general como una serie de funciones propias y valores temporales.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd

    Ejercicio 4 Edp Calor

    Cargado por

    Bryan Gustavo Paredes
    60 vistas4 páginas
    Este documento presenta la resolución de un problema de condiciones de contorno no homogéneas mediante el método de separación de variables. Primero se transforma el problema a uno con condiciones homogéneas mediante una función. Luego se determinan las funciones propias y valores propios del problema transformado. Finalmente, se obtiene la solución general como una serie de funciones propias y valores temporales.

    Descripción original:

    Ecuaciones de calor

    Título original

    ejercicio 4 edp calor

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Este documento presenta la resolución de un problema de condiciones de contorno no homogéneas mediante el método de separación de variables. Primero se transforma el problema a uno con condiciones homogéneas mediante una función. Luego se determinan las funciones propias y valores propios del problema transformado. Finalmente, se obtiene la solución general como una serie de funciones propias y valores temporales.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    60 vistas4 páginas

    Ejercicio 4 Edp Calor

    Cargado por

    Bryan Gustavo Paredes
    Este documento presenta la resolución de un problema de condiciones de contorno no homogéneas mediante el método de separación de variables. Primero se transforma el problema a uno con condiciones homogéneas mediante una función. Luego se determinan las funciones propias y valores propios del problema transformado. Finalmente, se obtiene la solución general como una serie de funciones propias y valores temporales.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    Está en la página 1de 4

    Gustavo Paredes, Carlos Jaramillo

    Primero debemos transformar a un problema de condiciones de contorno


    homogéneas. Para ello, buscamos una función 𝑤(𝑥, 𝑡), la cual transformará las
    condiciones de contorno en homogéneas. Tomemos la siguiente función
    𝑤(𝑥, 𝑡) = 𝑎𝑡 2 , entonces, considerando el teorema de superposición se sigue
    que
    𝑢(𝑥, 𝑡) = 𝑣(𝑥, 𝑡) + 𝑤(𝑥, 𝑡) = 𝑣(𝑥, 𝑡) + 𝑎𝑡 2

    donde 𝑣(𝑥, 𝑡) es la solución del problema con condiciones de contorno


    homogéneas.

    Para determinar el coeficiente 𝑎 reemplazamos la expresión anterior en las


    condiciones de contorno
    𝑢𝑥 (0, 𝑡) = 𝑣𝑥 (0, 𝑡) + 0 = 0
    𝒗𝒙 (𝟎, 𝒕) = 𝟎.

    𝑢(𝜋, 𝑡) = 𝑣(𝜋, 𝑡) + 𝑎𝑡 2 = 𝑡 2
    De esta última expresión se sigue que 𝑎 = 1, para que

    𝒗(𝜋, 𝒕) = 𝟎.
    Entonces, nuestra solución tiene la forma

    𝒖(𝒙, 𝒕) = 𝒗(𝒙, 𝒕) + 𝒕𝟐 (*)

    Con estas expresiones, finalmente obtenemos el problema con condiciones de


    contorno homogéneas, reemplazando (*) en la ecuación diferencial original.
    Así

    (𝑣(𝑥, 𝑡) + 𝑡 2 )𝑡 = (𝑣(𝑥, 𝑡) + 𝑡 2 )𝑥𝑥


    𝑣𝑥 (0, 𝑡) = 0, 𝑣(𝜋, 𝑡) + 𝑡 2 = 𝑡 2 ,

    { (𝑣(𝑥, 𝑡) + 𝑡 2 )|𝑡=0 = 0.

    𝒗𝒕 + 𝟐𝒕 = 𝒗𝒙𝒙
    𝒗𝒙 (𝟎, 𝒕) = 𝟎, 𝒗(𝜋, 𝒕) = 𝟎,
    {
    𝒗(𝒙, 𝟎) = 𝟎.
    Para resolver esta ecuación diferencial, suponemos que la solución tiene la
    forma 𝑒 𝛼𝑥 y la reemplazamos en la ecuación anterior

    −𝛼 2 𝑒 𝛼𝑥 = 𝜆𝑒 𝛼𝑥 , 𝑒 𝛼𝑥 ≠ 0 −𝛼 2 = 𝜆 𝛼 = √−𝜆 = 𝑖√𝜆
    De donde se obtiene que la solución del problema de Sturm-Liouville tiene la
    forma

    𝑣(𝑥) = 𝐶1𝑐𝑜𝑠(√𝜆𝑥) + 𝐶2𝑠𝑒𝑛(√𝜆𝑥) (**)

    Para determinar los valores de los coeficientes es necesario utilizar las


    condiciones de contorno

    𝑣𝑥 (𝑥) = −𝐶1√𝜆𝑠𝑒𝑛(√𝜆𝑥) + 𝐶2√𝜆𝑐𝑜𝑠(√𝜆𝑥)

    Si 𝑥 = 0, se obtiene que 𝑣𝑥 (0) = 𝐶2√𝜆 = 0, de donde se concluye que 𝐶2 =


    0.

    Si 𝑥 = 𝜋, se obtiene que 𝑣(𝜋) = 𝐶1𝑐𝑜𝑠(√𝜆𝐿) = 0, 𝐶1 ≠ 0, de donde se


    concluye que 𝑐𝑜𝑠(√𝜆𝐿) = 0. Resolviendo esta ecuación trigonométrica, se
    obtiene que
    𝜋 1 + 2𝑛
    √𝜆𝐿 = + 𝑛𝜋 = ( )𝜋
    2 2
    De donde se obtiene que, los valores propios cumplen con la identidad

    𝟏+𝟐𝒏 𝟐 𝝅𝟐
    𝝀𝒏 = ( ) 𝟐, 𝒏 = 𝟎, 𝟏, 𝟐, 𝟑, … …
    𝟐 𝜋

    Regresando a la expresión (**) se obtiene que las funciones propias tienen la


    forma

    𝟏 + 𝟐𝒏
    𝒗𝒏 (𝒙) = 𝒄𝒐𝒔 (( ) 𝒙) , 𝒏 = 𝟎, 𝟏, 𝟐, 𝟑, … …
    𝟐

    Luego de haber determinado las bases de la serie, lo que debemos realizar es


    transformar todas las funciones que aparecen en el problema dado en términos
    de las bases

    2𝑡𝟏 = 2𝑡 ∑ 𝑩𝒏 𝒗𝒏 (𝒙)
    𝒏=𝟎
    ∞

    (𝟏, 𝒗𝒎 (𝒙)) = ∑ 𝑩𝒏 (𝒗𝒏 (𝒙), 𝒗𝒎 (𝒙))


    𝒎,𝒏=𝟎

    Por la propiedad de ortogonalidad de las bases se tiene que

    (1, 𝑣𝑛 (𝑥)) = 𝐵𝑛 (𝑣𝑛 (𝑥), 𝑣𝑛 (𝑥))


    De donde
    (1, 𝑣𝑛 (𝑥))
    𝐵𝑛 =
    (𝑣𝑛 (𝑥), 𝑣𝑛 (𝑥))

    𝜋 𝜋
    1 + 2𝑛 2 1 + 2𝑛
    (1, 𝑣𝑛 (𝑥)) = ∫ 𝑐𝑜𝑠 (( ) 𝑥) 𝑑𝑥 = ( ) 𝑠𝑒𝑛 (( ) 𝑥)|
    0 2 1 + 2𝑛 2 0
    2 1 + 2𝑛 𝟐
    =( ) 𝑠𝑒𝑛 (( )) = ( ) (−𝟏)𝒏
    1 + 2𝑛 2 𝟏 + 𝟐𝒏
    𝜋
    2
    1 + 2𝑛 1 𝜋
    (𝑣𝑛 (𝑥), 𝑣𝑛 (𝑥)) = ∫ 𝑐𝑜𝑠 (( ) 𝑥) 𝑑𝑥 = ∫ [1 + 𝑐𝑜𝑠((1 + 2𝑛)𝑥)]𝑑𝑥
    0 2 2 0
    𝜋
    1 𝐿 𝝅
    = [𝐿 + 𝑠𝑒𝑛((1 + 2𝑛)𝑥)| ] =
    2 (1 + 2𝑛)𝜋 0
    𝟐

    2
    (1 + 2𝑛) (−1)𝑛
    𝐵𝑛 = 𝜋
    2
    𝟒 𝟏
    𝑩𝒏 = ( ) (−𝟏)𝒏
    𝟏 + 𝟐𝒏 𝝅
    De donde se observa que, la solución de nuestro problema presenta la forma
    ∞ ∞
    𝟏 + 𝟐𝒏
    𝒗(𝒙, 𝒕) = ∑ 𝑻𝒏 (𝒕)𝒗𝒏 (𝒙) = ∑ 𝑻𝒏 (𝒕)𝒄𝒐𝒔 (( ) 𝒙) (∗∗∗)
    𝟐
    𝒏=𝟎 𝒏=𝟎

    La expresión anterior debe ser reemplazada en la ecuación diferencial parcial


    del problema transformado con condiciones de contorno homogéneas, además,
    de considerar la misma en las condiciones iniciales. En nuestro caso se obtiene
    que:
    ∞ ∞ ∞

    ∑ 𝑻𝒏 ′(𝒕)𝒗𝒏 (𝒙) + 𝟐𝒕 ∑ 𝑩𝒏 𝒗𝒏 (𝒙) = −𝝀𝒏 ∑ 𝑻𝒏 (𝒕)𝒗𝒏 (𝒙) ,


    𝒏=𝟎 𝒏=𝟎 𝒏=𝟎

    ∑ 𝑻𝒏 (𝟎)𝒗𝒏 (𝒙) = 0
    𝒏=𝟎

    De aquí, igualando las bases, se obtiene la ecuación diferencial ordinaria con


    valor inicial para la parte temporal, así tenemos que

    𝑻′𝒏 (𝒕) + 𝟐𝒕𝑩𝒏 = −𝝀𝒏 𝑻𝒏 (𝒕)


    𝑻𝒏 (𝟎) = 𝟎
    Para resolver el problema anterior, es necesario encontrar la solución de la
    parte homogénea, para luego sumarle con la solución de la parte no homogénea

    𝑻′𝒏 (𝒕) = −𝝀𝒏 𝑻𝒏 (𝒕)


    Integrando la expresión anterior resulta que
    𝑇𝑛′ (𝑡)
    = −𝜆𝑛
    𝑇𝑛 (𝑡)
    𝑙𝑛(𝑇𝑛 (𝑡)) = −𝜆𝑛 𝑡 + 𝐴1

    𝑻𝒏 (𝒕) = 𝜶𝒆−𝝀𝒏 𝒕
    Para resolver la parte no homogénea, consideremos a la siguiente función
    𝑻𝒏 (𝒕) = 𝒂𝒏 + 𝒃𝒏 𝒕, la cual debe ser reemplazada en la ecuación diferencial
    ordinaria no homogénea
    𝑏𝑛 + 2𝑡𝐵𝑛 = −𝜆𝑛 (𝑎𝑛 + 𝑏𝑛 𝑡)
    De esta expresión es necesario igual las bases, con lo que se obtiene el siguiente
    sistema de ecuaciones

    𝑏𝑛 = −𝜆𝑛 𝑎𝑛

    2𝐶𝑛 = −𝜆𝑛 𝑏𝑛 = 𝜆𝑛 2 𝑎𝑛
    De donde se obtiene que
    𝟐𝑩𝒏
    𝒂𝒏 =
    𝝀𝒏 𝟐
    𝟐𝑩𝒏
    𝒃𝒏 = −
    𝝀𝒏
    Con ello finalmente, por el teorema de la superposición, la solución general de
    la ecuación diferencial ordinaria es

    𝑻𝒏 (𝒕) = 𝜶𝒆−𝝀𝒏 𝒕 + 𝒂𝒏 + 𝒃𝒏 𝒕
    Ahora, necesitamos determinar el valor de 𝜶, para ello utilizamos la condición
    inicial

    𝑇𝑛 (0) = 0

    𝑇𝑛 (0) = 𝛼𝑒 −𝜆𝑛 0 + 𝑎𝑛 + 𝑏𝑛 0 = 0
    𝛼 = −𝑎𝑛
    Finalmente, la solución para la parte temporal es

    𝑻𝒏 (𝒕) = −𝒂𝒏 𝒆−𝝀𝒏𝒕 + 𝒂𝒏 + 𝒃𝒏 𝒕


    Este resultado debe ser reemplazado en la expresión (***), con lo que se sigue
    que
    ∞ ∞

    𝒗(𝒙, 𝒕) = ∑ 𝑻𝒏 (𝒕)𝒗𝒏 (𝒙) = ∑(−𝒂𝒏 𝒆−𝝀𝒏𝒕 + 𝒂𝒏 + 𝒃𝒏 𝒕)𝒗𝒏 (𝒙)


    𝒏=𝟎 𝒏=𝟎

    Ahora, para obtener la solución del problema original (del problema con
    condiciones de contorno no homogéneas) reemplazamos la serie anterior en la
    expresión (*)

    𝒖(𝒙, 𝒕) = 𝒗(𝒙, 𝒕) + 𝒕𝟐
    ∞

    𝒖(𝒙, 𝒕) = ∑(−𝒂𝒏 𝒆−𝝀𝒏 𝒕 + 𝒂𝒏 + 𝒃𝒏 𝒕)𝒗𝒏 (𝒙) + 𝒕𝟐 .


    𝒏=𝟎

    También podría gustarte