Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 244 vistas

    Métodos de Diferencias Finitas (Ecuaciones Parabolicas) Crank-Nocolson ADI

    Cargado por

    Andy Carrera
    Este documento describe métodos numéricos para resolver ecuaciones diferenciales parciales (EDPs) como la ecuación de transmisión de calor. Explica el método de Crank-Nicolson para obtener una aproximación de segundo orden en el tiempo y el espacio. También describe el esquema ADI para resolver EDPs bidimensionales de manera eficiente usando matrices tridiagonales. El método ADI implica resolver la ecuación en dos pasos alternando entre direcciones implícitas y explícitas.

    Copyright:

    Attribution Non-Commercial (BY-NC)
    Descargue como DOC, PDF, TXT o lea en línea desde Scribd

    Métodos de Diferencias Finitas (Ecuaciones Parabolicas) Crank-Nocolson ADI

    Cargado por

    Andy Carrera
    244 vistas4 páginas
    Este documento describe métodos numéricos para resolver ecuaciones diferenciales parciales (EDPs) como la ecuación de transmisión de calor. Explica el método de Crank-Nicolson para obtener una aproximación de segundo orden en el tiempo y el espacio. También describe el esquema ADI para resolver EDPs bidimensionales de manera eficiente usando matrices tridiagonales. El método ADI implica resolver la ecuación en dos pasos alternando entre direcciones implícitas y explícitas.

    Título original

    Métodos de Diferencias Finitas (ecuaciones parabolicas) Crank-Nocolson ADI

    Derechos de autor

    © Attribution Non-Commercial (BY-NC)

    Formatos disponibles

    DOC, PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Este documento describe métodos numéricos para resolver ecuaciones diferenciales parciales (EDPs) como la ecuación de transmisión de calor. Explica el método de Crank-Nicolson para obtener una aproximación de segundo orden en el tiempo y el espacio. También describe el esquema ADI para resolver EDPs bidimensionales de manera eficiente usando matrices tridiagonales. El método ADI implica resolver la ecuación en dos pasos alternando entre direcciones implícitas y explícitas.

    Copyright:

    Attribution Non-Commercial (BY-NC)
    Descargue como DOC, PDF, TXT o lea en línea desde Scribd
    Descargar como doc, pdf o txt
    244 vistas4 páginas

    Métodos de Diferencias Finitas (Ecuaciones Parabolicas) Crank-Nocolson ADI

    Cargado por

    Andy Carrera
    Este documento describe métodos numéricos para resolver ecuaciones diferenciales parciales (EDPs) como la ecuación de transmisión de calor. Explica el método de Crank-Nicolson para obtener una aproximación de segundo orden en el tiempo y el espacio. También describe el esquema ADI para resolver EDPs bidimensionales de manera eficiente usando matrices tridiagonales. El método ADI implica resolver la ecuación en dos pasos alternando entre direcciones implícitas y explícitas.

    Copyright:

    Attribution Non-Commercial (BY-NC)
    Descargue como DOC, PDF, TXT o lea en línea desde Scribd
    Descargar como doc, pdf o txt
    Está en la página 1de 4

    Mtodos de Diferencias Finitas (ecuaciones parablicas)

    Mtodos Implcitos (mtodo de Crank-Nicolson)


    A pesar del mtodo implcito anterior ser estable y convergente, tiene el inconveniente de
    la aproximacin temporal ser apenas de 1ra orden mientras la aproximacin espacial es
    de 2da orden. El mtodo de Crank-Nicolson constituye una alternativa implcita que tiene
    2da orden en el tiempo y en el espacio.
    Para conseguir la aproximacin de 2da orden en el tiempo, el mtodo de Crank-Nicolson
    se basa en la aproximacin por diferencias finitas en un tiempo intermedio del rango t.
    As, la 1ra derivada temporal puede ser aproximada en el instante t
    l+1/2
    La 2da derivada espacial es tambin aproximada en el instante t
    l+1/2
    , haciendo el promedio
    de las aproximaciones en el inicio (t
    l
    ) y en el final (t
    l+1
    ) del paso de clculo t
    Substituyendo las aproximaciones en la EDP se obtiene:
    Molcula de clculo del mtodo Crank-Nicolson
    1 l l
    i i
    T T T
    t t
    +


    1 1 1 2
    1 1 1 1
    2 2 2
    2 2 1
    2
    + + +
    + +
    1
    + +
    + 1
    1
    ]
    l l l l l l
    i i i i i i
    T T T T T T T
    x x x
    1 1 1 1 2
    1 1 1 1
    2 2 2
    2 2
    2
    l l l l l l l l
    i i i i i i i i
    T T T T T T T T T T k
    k
    t t
    x x x
    + + + +
    + +
    1
    + +
    + 1

    1
    ]
    ( ) ( )
    1 1 1 1
    1 1 1 1
    2
    1 1 1
    1 1 1 1
    con 2 2
    2
    2 1 2 1
    l l l l l l l l
    i i i i i i i i
    l l l l l l
    i i i i i i
    k t
    T T T T T T T T
    x
    T T T T T T
    + + + +
    + +
    + + +
    + +

    1
    + + +
    ]

    + + + +
    Tal como en el caso del mtodo implcito ms simples, las condiciones de frontera pueden
    ser incorporadas en las ecuaciones del primero y del ltimo nodos interiores
    As, en el primero nodo (i = 1) se obtiene:
    En el ltimo nodo (i = m) se obtiene:
    A pesar de que las ecuaciones a las diferencias son ligeramente ms complejas do que
    las del mtodo simples, el sistema resultante es tambin tridiagonal por lo que puede ser
    solucionado de modo eficiente.
    EDPs en dos dimensiones espaciales
    La ecuacin de transmisin de calor puede ser aplicada a ms do que una dimensin.
    Para dos dimensiones se tiene:
    Mtodos explcitos e implcitos convencionales
    Puede obtenerse una solucin explcita substituyendo en la ecuacin anterior las
    derivadas por diferencias finitas tal como ha sido hecho para el caso unidimensional
    Este mtodo est condicionado por un criterio de estabilidad
    Los mtodos implcitos permiten obtener soluciones estables. Pero, la aplicacin directa
    de mtodos como el de Crank-Nicolson conducen a la resolucin de un sistema de m x n
    ecuaciones. Adems, la aplicacin directa de mtodos implcitos a dos o tres dimensiones
    implica la prdida de la propiedad tridiagonal, que permite la resolucin eficiente del
    sistema de ecuaciones
    ( )
    1 1
    0 0
    l l
    T f t
    + +

    ( )
    1 1
    1 1
    l l
    m m
    T f t
    + +
    + +

    ( )
    ( )
    ( )
    ( )
    1 1 1
    1 2 0 1 2 0
    2 1 2 1
    l l l l l l
    T T f t T T f t
    + + +
    + + + +
    ( )
    ( )
    ( )
    ( )
    1 1 1
    1 1 1 1
    2 1 2 1
    l l l l l l
    m m m m m m
    T T f t T T f t
    + + +
    + +
    + + + + +
    2 2
    2 2
    T T T
    k
    t
    x y
    _

    +


    ,
    1
    2 2
    , , 1, , 1, , 1 , , 1
    2 2 2 2
    2 2
    l l l l l l l l
    i j i j i j i j i j i j i j i j
    T T T T T T T T
    T T T
    k k
    t t
    x y x y
    +
    + +
    _
    _ + +


    + +



    ,
    ,
    2 2
    1
    8
    x y
    t
    k
    +

    El esquema A.D.I.
    El esquema A.D.I. (alternating-direction implicit) permite la resolucin de ecuaciones
    parablicas en dos dimensiones usando matrices tridiagonales. Para eso, cada intervalo
    de tiempo es ejecutado en dos pasos. En el 1er paso la EDP es aproximada por
    El 1er paso consiste entonces en pasar del instante t
    l
    para el instante t
    l+1/2
    a travs de una
    aproximacin centrada explcita en la direccin x y de una aproximacin centrada implcita
    en la direccin y
    Para una malla cuadrada (x = y), la ecuacin anterior puede ser escrita
    El sistema de ecuaciones resultante es tridiagonal, siendo por eso posible encontrar su
    solucin de modo eficiente
    El 2do paso consiste entonces en pasar del instante t
    l+1/2
    para el instante t
    l+1
    a travs de
    una aproximacin centrada implcita en la direccin x y de una aproximacin centrada
    explcita en la direccin y
    Para una malla cuadrada (x = y), la ecuacin anterior puede ser escrita
    El sistema de ecuaciones resultante es tambin tridiagonal
    Note que el esquema A.D.I. solo conduce a sistemas tridiagonales si es aplicado en la
    direccin que es implcita, i.e., en el 1ro paso segundo y y en el 2do paso segundo x
    1 2 1 2 1 2 1 2
    2 2
    , , , 1, , 1, , 1 , 1
    2 2 2 2
    2 2
    2
    l l l l l l l l
    i j i j i j i j i j i j i j i j
    T T T T T T T T
    T T T
    k k
    t t
    x y x y
    + + + +
    + +
    _
    + _ +


    + +



    ,
    ,
    1 2 1 2 1 2 1 2
    2 2
    , , , 1, , 1, , 1 , 1
    2 2 2 2
    2 2
    2
    l l l l l l l l
    i j i j i j i j i j i j i j i j
    T T T T T T T T
    T T T
    k k
    t t
    x y x y
    + + + +
    + +
    _
    + _ +


    + +



    ,
    ,
    ( ) ( )
    2
    1 2 1 2 1 2
    1, , 1, , , 1 , 1
    2 1 2 1
    k t
    x
    l l l l l l
    i j i j i j i j i j i j
    T T T T T T

    + + +
    + +
    + + + +
    1 2 1 2 1 2 1 2 1 1 1 1
    2 2
    , , , 1, , 1, , 1 , 1
    2 2 2 2
    2 2
    2
    l l l l l l l l
    i j i j i j i j i j i j i j i j
    T T T T T T T T
    T T T
    k k
    t t
    x y x y
    + + + + + + + +
    + +
    _
    + _ +


    + +



    ,
    ,
    ( ) ( )
    2
    1 2 1 2 1 2 1 1 1
    1, , 1, , , 1 , 1
    2 1 2 1
    k t
    x
    l l l l l l
    i j i j i j i j i j i j
    T T T T T T

    + + + + + +
    + +
    + + + +
    El mtodo A.D.I. es solo uno de los mtodos pertenecientes al grupo de los designados
    mtodos splitting (divisin de la aproximacin implcita en diversos pasos, dos para dos
    dimensiones, tres para tres dimensiones). Algunos permiten ultrapasar las dificultades de
    aplicacin del A.D.I. a problemas tridimensionales

    También podría gustarte