Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 1K vistas

    Esfera

    Cargado por

    Gabriela Gemio Lujan
    Este documento describe el flujo de un fluido alrededor de una esfera que gira lentamente. Se obtienen expresiones para la distribución de presión y velocidad en el fluido asumiendo flujo laminar. Se resuelve la ecuación de continuidad y las ecuaciones de cantidad de movimiento en coordenadas esféricas para derivar una ecuación diferencial que describe el perfil de velocidades en el fluido.

    Copyright:

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

    Esfera

    Cargado por

    Gabriela Gemio Lujan
    1K vistas15 páginas
    Este documento describe el flujo de un fluido alrededor de una esfera que gira lentamente. Se obtienen expresiones para la distribución de presión y velocidad en el fluido asumiendo flujo laminar. Se resuelve la ecuación de continuidad y las ecuaciones de cantidad de movimiento en coordenadas esféricas para derivar una ecuación diferencial que describe el perfil de velocidades en el fluido.

    Título original

    esfera

    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 describe el flujo de un fluido alrededor de una esfera que gira lentamente. Se obtienen expresiones para la distribución de presión y velocidad en el fluido asumiendo flujo laminar. Se resuelve la ecuación de continuidad y las ecuaciones de cantidad de movimiento en coordenadas esféricas para derivar una ecuación diferencial que describe el perfil de velocidades en el fluido.

    Copyright:

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

    Esfera

    Cargado por

    Gabriela Gemio Lujan
    Este documento describe el flujo de un fluido alrededor de una esfera que gira lentamente. Se obtienen expresiones para la distribución de presión y velocidad en el fluido asumiendo flujo laminar. Se resuelve la ecuación de continuidad y las ecuaciones de cantidad de movimiento en coordenadas esféricas para derivar una ecuación diferencial que describe el perfil de velocidades en el fluido.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    de 15

    FLUJO CERCA DE UNA ESFERA QUE GIRA

    LENTAMENTE

    Una esfera solida de radio R gira


    lentamente a velocidad angular
    constante Ω en un gran cuerpo de
    fluido en reposo. Obtener
    expresiones para la distribución
    de presión y de velocidad en el
    fluido. Se supone que la esfera
    gira lo suficientemente lento de
    modo que es adecuado utilizar la
    versión de flujo reptante.
    Solución

    DEFINICION DE CAMPOS
    𝒗𝒓 , 𝒗𝜽 , 𝒗𝝓 𝒗𝝓 = 𝒇(𝒓, 𝜽)
    𝑷𝒓 , 𝑷𝜽 , 𝑷𝝓
    𝑷𝝓 = 𝟎
    𝑷𝜽 = 𝒇(𝜽)
    𝑷𝒓 = 𝒇(𝒓)
    𝒈𝒓 , 𝒈𝜽 , 𝒈𝝓
    𝒈𝜽
    𝒈𝝓
    Ecuación de continuidad:
    COORDENADAS ESFÉRICAS (r, θ, ϕ):
    𝝏𝝆 𝟏 𝝏 𝟏 𝝏 𝟏 𝝏
    + 𝟐 (𝝆𝒓𝟐 𝒗𝒓 ) + (𝝆𝒗𝜽 𝑺𝒆𝒏𝜽) + (𝝆𝒗𝝓 ) = 𝟎
    𝝏𝒕 𝒓 𝝏𝒓 𝒓𝑺𝒆𝒏𝜽 𝝏𝜽 𝒓𝑺𝒆𝒏𝜽 𝝏𝝓
    𝝏𝝆
    =𝟎 𝑬𝑺𝑻𝑨𝑫𝑶 𝑬𝑺𝑻𝑨𝑪𝑰𝑶𝑵𝑨𝑹𝑰𝑶
    𝝏𝒕
    𝒗𝜽 = 𝟎
    𝒗𝒓 = 𝟎
    NOS QUEDA:
    𝟏 𝝏
    (𝝆𝒗𝝓 ) = 𝟎
    𝒓𝑺𝒆𝒏𝜽 𝝏𝝓
    SUPONIENDO DENSIDAD CONSTANTE:
    𝟏 𝝏𝒗𝝓
    =𝟎
    𝒓𝑺𝒆𝒏𝜽 𝝏𝝓
    𝝏𝒗𝝓
    =𝟎
    𝝏𝝓
    INTEGRANDO:
    𝒗𝝓 = 𝑪𝑶𝑵𝑺𝑻𝑨𝑵𝑻𝑬
    CONSIDERACIONES
    FLUJO REPTANTE: velocidades se igualan a cero debido a que nos
    indica que la esfera gira lo suficientemente lento.
    ECUACION DE CANTIDAD DE MOVIMIENTO:
    En función de los gradientes de velocidad para un fluido newtoniano de 𝝆 𝒚 𝝁
    constante:
    Componente r:
    𝝏𝒗𝒓 𝝏𝒗𝒓 𝒗𝜽 𝝏𝒗𝒓 𝒗𝝓 𝝏𝒗𝒓 𝒗𝜽 𝟐 + 𝒗𝝓 𝟐
    𝝆( + 𝒗𝒓 + + − )
    𝝏𝒕 𝝏𝒓 𝒓 𝝏𝜽 𝒓𝑺𝒆𝒏𝜽 𝝏𝝓 𝒓
    𝝏𝑷 𝟐 𝟐 𝝏𝒗𝜽 𝟐 𝟐 𝝏𝒗𝝓
    =− + 𝝁 (𝛁𝟐 𝒗𝒓 − 𝟐 𝒗𝒓 − 𝟐 − 𝟐 𝒗𝜽 𝒄𝒐𝒕𝜽 − 𝟐 )
    𝝏𝒓 𝒓 𝒓 𝝏𝜽 𝒓 𝒓 𝑺𝒆𝒏𝜽 𝝏𝝓
    + 𝝆𝓰𝒓 (𝑫)

    Componente 𝜽:
    𝝏𝒗𝜽 𝝏𝒗𝜽 𝒗𝜽 𝝏𝒗𝜽 𝒗𝝓 𝝏𝒗𝜽 𝒗𝒓 𝒗𝜽 𝒗𝝓 𝟐 𝒄𝒐𝒕𝜽
    𝝆( + 𝒗𝒓 + + + + )
    𝝏𝒕 𝝏𝒓 𝒓 𝝏𝜽 𝒓𝑺𝒆𝒏𝜽 𝝏𝝓 𝒓 𝒓
    𝟏 𝝏𝑷 𝟐 𝝏𝒗𝒓 𝒗𝜽 𝟐𝑪𝒐𝒔𝜽 𝝏𝒗𝝓
    =− + 𝝁 (𝛁 𝟐𝒗𝜽 − 𝟐 − 𝟐 − )
    𝒓 𝝏𝜽 𝒓 𝝏𝜽 𝒓 𝑺𝒆𝒏𝟐 𝜽 𝒓𝟐 𝑺𝒆𝒏𝟐 𝜽 𝝏𝝓
    + 𝝆𝓰𝜽 (𝑬)
    Componente ϕ:
    𝝏𝒗𝝓 𝝏𝒗𝝓 𝒗𝜽 𝝏𝒗𝝓 𝒗𝝓 𝝏𝒗𝝓 𝒗𝝓 𝒗𝒓 𝒗𝜽 𝒗𝝓
    𝝆( + 𝒗𝒓 + + + + 𝒄𝒐𝒕𝜽)
    𝝏𝒕 𝝏𝒓 𝒓 𝝏𝜽 𝒓𝑺𝒆𝒏𝜽 𝝏𝝓 𝒓 𝒓
    𝟏 𝝏𝑷 𝒗𝝓 𝟐 𝝏𝒗𝒓 𝟐𝑪𝒐𝒔𝜽 𝝏𝒗𝜽
    =− − 𝝁 (𝛁 𝟐 𝒗𝝓 − 𝟐 + 𝟐 + )
    𝒓𝑺𝒆𝒏𝜽 𝝏𝝓 𝒓 𝑺𝒆𝒏𝜽 𝒓 𝑺𝒆𝒏𝜽 𝝏𝝓 𝒓𝟐 𝑺𝒆𝒏𝜽 𝝏𝝓
    + 𝝆𝓰𝝓 (𝑭)

    En estas ecuaciones:

    𝟏 𝝏 𝟐 𝝏 𝟏 𝝏 𝝏 𝟏 𝝏𝟐
    𝛁𝟐 = (𝒓 ) + (𝑺𝒆𝒏𝜽 ) + ( )
    𝒓𝟐 𝝏𝒓 𝝏𝒓 𝒓𝟐 𝑺𝒆𝒏𝜽 𝝏𝜽 𝝏𝜽 𝒓𝟐 𝑺𝒆𝒏𝟐 𝜽 𝝏𝝓𝟐
    Por lo tanto simplificando:
    Componente r:
    𝝏𝑷
    𝟎=− + 𝝆𝓰𝒓
    𝝏𝒓
    Componente 𝜽:
    𝟏 𝝏𝑷
    𝟎=− + 𝝆𝓰𝜽
    𝒓 𝝏𝜽
    Componente ϕ:
    𝟏 𝝏 𝝏𝒗𝝓 𝟏 𝝏 𝟏 𝝏
    (𝒓𝟐 )+ 𝟐 ( (𝒗 ∗ 𝑺𝒆𝒏𝜽)) = 𝟎
    𝟐
    𝒓 𝝏𝒓 𝝏𝒓 𝒓 𝝏𝜽 𝑺𝒆𝒏𝜽 𝝏𝜽 𝝓

    En estas ecuaciones:

    𝟏 𝝏 𝝏𝒗𝝓 𝟏 𝝏 𝟏 𝝏 𝟏 𝝏𝟐 𝒗𝝓
    𝛁 𝟐 𝒗𝝓 = (𝒓 𝟐 ) + ( (𝒗 ∗ 𝑺𝒆𝒏𝜽)) + ( )
    𝒓𝟐 𝝏𝒓 𝝏𝒓 𝒓𝟐 𝝏𝜽 𝑺𝒆𝒏𝜽 𝝏𝜽 𝝓 𝒓𝟐 𝑺𝒆𝒏𝟐 𝜽 𝝏𝝓𝟐
    Por lo tanto:
    𝟏 𝝏 𝝏𝒗𝝓 𝟏 𝝏 𝟏 𝝏
    𝛁 𝟐 𝒗𝝓 = (𝒓𝟐 )+ 𝟐 ( (𝒗 ∗ 𝑺𝒆𝒏𝜽))
    𝟐
    𝒓 𝝏𝒓 𝝏𝒓 𝒓 𝝏𝜽 𝑺𝒆𝒏𝜽 𝝏𝜽 𝝓

    Reemplazando:
    𝟏 𝝏 𝟐
    𝝏𝒗𝝓 𝟏 𝝏 𝟏 𝝏
    (𝒓 )+ 𝟐 ( (𝒗 ∗ 𝑺𝒆𝒏𝜽)) = 𝟎
    𝒓𝟐 𝝏𝒓 𝝏𝒓 𝒓 𝝏𝜽 𝑺𝒆𝒏𝜽 𝝏𝜽 𝝓
    HALLANDO UNA ECUACIÓN PARA EL PERFIL DE VELOCIDADES:
    Cambio de variable:

    𝒗𝝓 = 𝒁(𝒓) ∗ 𝒔𝒆𝒏𝜽
    𝒗𝝓
    𝒁(𝒓) =
    𝒔𝒆𝒏𝜽
    Simplificando la ecuación:
    𝝏 𝝏𝒗𝝓 𝝏 𝟏 𝝏
    (𝒓𝟐 )+ ( (𝒗 ∗ 𝑺𝒆𝒏𝜽)) = 𝟎
    𝝏𝒓 𝝏𝒓 𝝏𝜽 𝑺𝒆𝒏𝜽 𝝏𝜽 𝝓

    Para el primer termino

    𝝏 𝝏𝒗𝝓 𝝏 𝒅(𝒁(𝒓) ∗ 𝒔𝒆𝒏𝜽)


    (𝒓𝟐 )= (𝒓𝟐 )
    𝝏𝒓 𝝏𝒓 𝝏𝒓 𝝏𝒓
    𝝏 𝝏𝒗𝝓 𝝏 𝒅𝒁(𝒓)
    (𝒓𝟐 ) = 𝒔𝒆𝒏𝜽 (𝒓𝟐 )
    𝝏𝒓 𝝏𝒓 𝝏𝒓 𝝏𝒓

    Para el segundo termino

    𝒗𝝓 ∗ 𝒔𝒆𝒏𝜽 = 𝒁(𝒓) ∗ 𝒔𝒆𝒏𝟐 𝜽


    𝝏
    (𝒗 ∗ 𝑺𝒆𝒏𝜽) = 𝟐 ∗ 𝒁 ∗ 𝒔𝒆𝒏𝜽 ∗ 𝑪𝒐𝒔𝜽
    𝝏𝜽 𝝓
    𝟏 𝝏
    ( (𝒗 ∗ 𝑺𝒆𝒏𝜽)) = 𝟐 ∗ 𝒁(𝒓) ∗ 𝑪𝒐𝒔𝜽
    𝑺𝒆𝒏𝜽 𝝏𝜽 𝝓
    𝝏 𝟏 𝝏
    ( (𝒗 ∗ 𝑺𝒆𝒏𝜽)) = −𝟐 ∗ 𝒁(𝒓) ∗ 𝑺𝒆𝒏𝜽
    𝝏𝜽 𝑺𝒆𝒏𝜽 𝝏𝜽 𝝓
    Reemplazando :
    𝝏 𝒅𝒁(𝒓)
    𝒔𝒆𝒏𝜽 (𝒓𝟐 ) − 𝟐 ∗ 𝒁(𝒓) ∗ 𝑺𝒆𝒏𝜽 = 𝟎
    𝝏𝒓 𝝏𝒓
    𝝏 𝒅𝒁(𝒓)
    (𝒓𝟐 ) − 𝟐 ∗ 𝒁(𝒓) = 𝟎
    𝝏𝒓 𝝏𝒓

    Resolviendo mediante el método de Euler


    𝝏 𝒅𝒁(𝒓)
    (𝒓𝟐 ) − 𝟐 ∗ 𝒁(𝒓) = 𝟎
    𝝏𝒓 𝝏𝒓

    𝒁(𝒓) = 𝑪𝟏 𝒓𝒎𝟏 + 𝑪𝟐 𝒓𝒎𝟐


    Parámetro de aproximación:
    𝒁(𝒓) = 𝒓𝒎
    𝒅𝒁(𝒓)
    = 𝒎 ∗ 𝒓𝒎−𝟏
    𝝏𝒓

    Reemplazando:
    𝝏 𝟐
    (𝒓 𝒎 ∗ 𝒓𝒎−𝟏 ) − 𝟐 ∗ 𝒓𝒎 = 𝟎
    𝝏𝒓
    Operando
    𝝏 𝒎−𝟏+𝟐
    𝒎 (𝒓 ) − 𝟐 ∗ 𝒓𝒎 = 𝟎
    𝝏𝒓
    𝝏 𝒎+𝟏
    𝒎 (𝒓 ) − 𝟐 ∗ 𝒓𝒎 = 𝟎
    𝝏𝒓
    Derivando:
    𝝏 𝒎+𝟏
    (𝒓 ) = (𝒎 + 𝟏)𝒓𝒎+𝟏−𝟏 = (𝒎 + 𝟏)𝒓𝒎
    𝝏𝒓
    𝒎 ∗ (𝒎 + 𝟏)𝒓𝒎 − 𝟐 ∗ 𝒓𝒎 = 𝟎
    Factorizamos 𝒓𝒎
    𝒎 ∗ (𝒎 + 𝟏) − 𝟐 = 𝟎
    Resolviendo
    𝒎𝟐 + 𝒎 − 𝟐 = 𝟎
    (𝒎 − 𝟏)(𝒎 + 𝟐) = 𝟎
    𝒎𝟏 = 𝟏
    𝒎𝟐 = −𝟐
    Reemplazamos en la formula:

    𝒁(𝒓) = 𝑪𝟏 𝒓𝟏 + 𝑪𝟐 𝒓−𝟐
    𝑪𝟐
    𝒁(𝒓) = 𝑪𝟏 ∗ 𝒓 +
    𝒓𝟐
    Reemplazando:
    𝒗𝝓
    𝒁(𝒓) =
    𝒔𝒆𝒏𝜽
    𝑪𝟐
    𝒗𝝓 = (𝑪𝟏 ∗ 𝒓 + 𝟐 ) ∗ 𝒔𝒆𝒏𝜽
    𝒓
    Condiciones frontera:
    𝒓=𝑹 𝒗𝝓 = 𝑹𝛀 ∗ 𝒔𝒆𝒏𝜽
    𝒓=∞ 𝒗𝝓 = 𝟎
    APLICANDO SEGUNDA CONDICION
    𝑪𝟐
    𝟎 = (𝑪𝟏 ∗ 𝒓 + 𝟐 ) ∗ 𝒔𝒆𝒏𝜽
    𝒓
    𝑪𝟐
    𝑪𝟏 = − =𝟎
    𝒓𝟑
    APLICANDO PRIMERA CONDICION
    𝑪𝟐
    𝑹𝛀 ∗ 𝒔𝒆𝒏𝜽 = (+ 𝟐 ) ∗ 𝒔𝒆𝒏𝜽
    𝑹

    𝑪𝟐 = 𝑹𝟑 𝛀

    REEMPLAZNDO

    𝑹𝟑 𝛀
    𝒗𝝓 = ( 𝟐 ) ∗ 𝒔𝒆𝒏𝜽
    𝒓

    También podría gustarte