Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 53 vistas

    Examen Extraordinario

    Cargado por

    LC Sosa
    Este documento presenta un análisis del problema de cavidad impulsado por una placa a través de 6 secciones. Se describe la teoría del flujo potencial, la vorticidad, la discretización de las ecuaciones de Navier-Stokes mediante diferencias finitas, y se plantea el problema específico de una cavidad con una placa en movimiento. Los resultados muestran la función de corriente y el campo de velocidades obtenidos mediante un programa de código.

    Copyright:

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

    Examen Extraordinario

    Cargado por

    LC Sosa
    53 vistas11 páginas
    Este documento presenta un análisis del problema de cavidad impulsado por una placa a través de 6 secciones. Se describe la teoría del flujo potencial, la vorticidad, la discretización de las ecuaciones de Navier-Stokes mediante diferencias finitas, y se plantea el problema específico de una cavidad con una placa en movimiento. Los resultados muestran la función de corriente y el campo de velocidades obtenidos mediante un programa de código.

    Descripción original:

    dinamica de fluidos computacionales

    Título original

    Examen extraordinario

    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 un análisis del problema de cavidad impulsado por una placa a través de 6 secciones. Se describe la teoría del flujo potencial, la vorticidad, la discretización de las ecuaciones de Navier-Stokes mediante diferencias finitas, y se plantea el problema específico de una cavidad con una placa en movimiento. Los resultados muestran la función de corriente y el campo de velocidades obtenidos mediante un programa de código.

    Copyright:

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

    Examen Extraordinario

    Cargado por

    LC Sosa
    Este documento presenta un análisis del problema de cavidad impulsado por una placa a través de 6 secciones. Se describe la teoría del flujo potencial, la vorticidad, la discretización de las ecuaciones de Navier-Stokes mediante diferencias finitas, y se plantea el problema específico de una cavidad con una placa en movimiento. Los resultados muestran la función de corriente y el campo de velocidades obtenidos mediante un programa de código.

    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 11

    Problema de Cavidad

    Luis Cuauhtémoc Sosa Manzo

    June 25, 2018

    1 Flujo de la Cavidad impulsado por una placa


    (Lid Driven Cavity Flow)
    El ujo de la cavidad impulsado por una placa es probablemente uno de los
    problemas de uidos más estudiamos en la dinámica de uidos computacionales.
    La simplicidad de la geometría hace al problema fácil de computar y aplicar
    condiciones de frontera. Aunque el problema parece simple en muchas formas,
    el ujo en la cavidad retiene toda la física del ujo con vórtices contrarrotanques
    que aparecen en la esquina de la cavidad.
    Los problemas de cavidad sirven como un punto de referencia para problemas
    de métodos numéricos en términos de exactitud, eciencia numérica, etc. En
    la literatura es posible encontrar numerosos estudios de ujo de la cavidad. El
    estudio numérico en materia de ujo de cavidad puede ser agrupado básicamente
    en tres categorías;

    1. En la primer categoría de estudio, la solución estable del ujo en la cavidad


    es buscada. En este tipo de estudios la solución numérica de la ecuaciones
    estables de Navier-Stokes para ujo incompresible son presentadas a distinto
    número de Reynolds.
    2. En la segunda categoría de estudio, la bifurcación del ujo en una cavidy
    desde un régimen estable hasta uno inestable es estudiada. En estos estudios
    el análisis de estabilidad hidrodinámica es realizado a un número de Reynolds
    en el cual la bifucación de Hopf se presente en el ujo.
    3. En la tercera categoría, es estudiada la transición de estable a inestable de
    un ujo a través de Simulación Numérica Directa (DNS) y el número de
    Reynolds en la transición es presentado

    1
    2 Teoría de Flujo Potencial
    La teoría de ujo potencial pretende describir el comportamiento cinemático de
    los uidos basándose en el concepto matemático de función potencial, asegu-
    rando que el campo de velocidades (que es un campo vectorial) del ujo de un
    uido es igual al gradiente de una función potencial que determina el movimiento
    de dicho uido:
    ~ = −∇φ
    V
    donde el campo de velocidades queda denido como:
    ~ = uî + v ĵ + wk̂
    V

    El signo menos en la ecuación de arriba es sólo una convención de signos sobre


    la denición de φ. φ Puede denirse sin el signo menos, y la formulación que se
    obtendría sería la misma. A un uido que se comporta según esta teoría se le
    denomina uido potencial, que da lugar a un ujo potencial.
    Una de las primeras personas en aplicar esta formulación para el ujo de un
    uido fue D'Alembert. Estudió la fuerza de resistencia producida por un ujo
    de uido sobre un cuerpo que se oponía a éste en dos dimensiones cuando la
    solución a este problema era completamente desconocida y Newton, a pesar de
    haberlo estudiado, no había llegado a conclusiones satisfactorias.
    D'Alembert denió la función de corriente, ψ , para describir la trayectoria
    que tuviera cada partícula de un uido a través del tiempo. Esta función cor-
    riente está determinada, en el plano, por dos variables espaciales y para cada
    valor de ψ la igualdad ψ = ψ(x, y) determina un lugar geométrico llamado línea
    de corriente.
    Primeramente deniremos la función corriente en el plano, para luego ex-
    plicar sus características. La función ψ se dene como aquella que cumple con
    las siguientes condiciones:
    ∂ψ ∂ψ
    =u , = −v
    ∂y ∂x

    2
    3 Vorticidad
    La vorticidad es una magnitud física empleada en mecánica de uidos para
    cuanticar la rotación de un uido.
    Matemáticamente la vorticidad es el campo vectorial denido por el rota-
    cional del campo de velocidades:
    ~
    ω =∇×V

    La presencia de vorticidad en un uido siempre implica la rotación de las partícu-


    las uidas, acompañada o no de alguna deformación transversal. En un u-
    ido real su existencia está íntimamente ligada a las tensiones tangenciales. La
    ecuación que permite estudiar la cinética de este campo (llamada ecuación de
    transporte de vorticidad) se obtiene tomando el rotacional a ambos lados de la
    ecuación de momentum de las Ecuaciones de Navier-Stokes.
    La vorticidad se origina fundamentalmente en los contornos sólidos debido
    a que los uidos no son capaces de deslizar sobre ellos, y luego se propaga al
    interior del uido siguiendo la ley de variación descrita por la Ecuación 2. El
    primer término corresponde a la variación de vorticidad por deformación de
    las líneas vorticosas. Este fenómeno ocurre tanto en uidos viscosos como no
    viscosos, sin embargo es un hecho notable que cuando el uido es no viscoso
    (ideal) esta es la única forma en que la vorticidad puede variar. Tal como
    lo demostró Kelvin en uno de sus teoremas, esta variación ocurre siempre de
    manera que el ujo de vorticidad asociado a una supercie abierta que se mueve
    con el uido permanece constante, lo cual también implica que la variación de
    la circulación ? de la velocidad a lo largo del contorno de esa misma supercie
    sea nula:
    dΓ
    =0
    dt

    3
    4 Discretización de las ecuaciones
    El método de las diferencias nitas fue el primer método aplicada para obtener
    la solución numérica de unas ecuaciones diferenciales. La primera persona en
    emplearlo fue Euler (probablemente en 1768). Este método emplea el desarrollo
    en seres de Taylor para la discretización de las derivadas de las variables de
    ujo; por ejemplo, si se quiere calcular la derivada de una función escalar U (x)
    evaluada en un punto x0 , tenemos:

    ∂U ∆x2 ∂2U
    U (x0 + ∆x) = U (x0 ) + ∆x + + + ... (1)
    ∂x 2 ∂x2

    Con esto, se tiene que la primera derivada de U (x) puede aproximarse por:

    ∂U U (x0 + ∆x) − U (x0 )


    = + O(∆x) (2)
    ∂x ∆

    Esta aproximación es de primer orden, porque el truncamiento del error es


    proporcional al término de mayor orden del resto de la serie. Este mismo proced-
    imiento se puede emplear para obtener aproximaciones más exactas del método
    de las diferencias nitas, sin más que aproximar con más términos de la serie de
    Taylor que se obtiene. Una ventaja importante del método de las diferencias ni-
    tas es su simplicidad teórica y la posibilidad de aumentar la precisión del método
    sin más que aumentar el orden de aproximación de las derivadas. Por otro lado,
    este método requiere una malla estructurada, con lo cual el campo de aplicación
    se ve reducido a algunos problemas. Normalmente con geometrías sencillas.
    Además, este método puede aplicarse directamente en coordenadas curvilíneas,
    teniéndose que transformar las ecuaciones de Navier-Stokes de curvilíneas a
    cartesianas para la resolución del problema.

    5 Formulación de las diferencias nitas de ecua-


    ciones diferenciales
    Las diferencias nitas de ecuaciones diferenciales plantean que se pueden reem-
    plazar los términos de derivadas por diferencias, en donde las derivadas se en-
    tienden como:
    ∆f f (x + ∆x) − f (x)
    f 0 (x) = lim = lim (3)
    ∆x→∞ ∆x ∆x→∞ ∆x

    4
    Recordando que la derivada representa a la pendiente de una recta tangente
    a la curva en un punto, entonces la derivada puede aproximarse por:

    f (x + ∆x) − f (x)
    f 0 (x) ≈ (4)
    ∆x

    6 Planteamiento del problema


    Se plantea conocer la Función de corriente y el campo de velocidades formado
    por un uido a través de un problema de cavidad, donde una placa en la parte
    superior se mueve a una velocidad u0 = 5 ms
    Se requiere hacer uso de las condiciones de frontera de Newman y de Dirichlet
    para delimitar las paredes tanto internas como externas para la resolución del
    problema.

    5
    7 Diagrama de ujo y Código

    Figure 1: Diagrama de ujo

    6
    Figure 2: Primera parte del código

    7
    Figure 3: Segunda parte del código

    8
    8 Resultados
    A continuación, se presentan los resultado de Función de corriente y Velocidades
    obtenidas del programa.

    Figure 4: Función de corriente

    9
    Figure 5: Campo de velocidades

    10
    9 Conclusiones

    En los resultados obtenidos en Octave, se aprecia el comportamiento de la Fun-


    ción de corriente y las velocidades en la cavidad. Cabe mencionar que el prob-
    lema es una cavidad muy pequeña y se maneja una velocidad de desplazamiento
    de la placa relativmente grande, además de que la distacia entre nodos es pe-
    queña por lo cual las condiciones de velocidad cero, son poco notorias en la
    cavidad.

    10 Bibliografía

    Anderson, J. (1995). Computational Fluid Dynamics. Singapore: Mc Graw


    Hill.

    Frank M. White (2008). Mecánica de uidos. McGraw Hill.

    Hong G. Im (2001). A Finite Dierence Code for the Navier-Stokes Equa-


    tions in Vorticity/Stream Function Formulation. University of Michigan.

    11

    También podría gustarte