Problema 2B.7 - BIRD

Problema 2B.
7
Flujo anular con cilindro interior moviéndose axialmente (ver Fig. 2B.7). Una varilla
cilíndrica de radio R se mueve axialmente con velocidad 𝑣𝑧 = 𝑣0 a lo largo del eje de
una cavidad cilíndrica de radio R como visto en la figura. La presión en ambos extremos
de la cavidad es la misma, por lo que el fluido se mueve a través de la región anular
únicamente debido al movimiento de la barra.
a) Encuentre la distribución de velocidades en la región anular estrecha.

b) Encuentre la tasa de masa a través de la región anular.
c) Obtenga la fuerza viscosa que actúa sobre la barra sobre la longitud L.
d) Demuestre que el resultado en c) se puede escribir como fórmula de “rendija
plana” multiplicada por una “corrección de curvatura”. Problemas de este tipo
surgen al estudiar el desempeño de troqueles de revestimiento de alambre.
Solución:
A diferencia del problema anterior, aquí el fluido fluye horizontalmente, por lo que la
gravedad no influirá en su velocidad. Si asumimos condiciones de contorno sin
deslizamiento, entonces la velocidad del fluido en la barra cilíndrica 𝑟 = 𝑘𝑅 es 𝑣0 y la
velocidad del fluido en la pared exterior 𝑟 = 𝑅 es 0.
Parte (a)
Elegimos un sistema de coordenadas cilíndricas con el eje 𝑧 positivo apuntando hacia la
derecha, la dirección del cilindro se mueve hacia adentro. Luego, el fluido fluye en la
dirección 𝑧 y varía en función de radio desde el eje del cilindro.
𝑣𝑧 = 𝑣𝑧 (𝑟)
Como resultado, las condiciones de contorno se describen de la siguiente manera:

B. C. 1: 𝑣𝑧 = 𝑣0 cuando 𝑟 = 𝑘𝑅
B. C. 2: 𝑣𝑧 = 0 cuando 𝑟 = 𝑅
La gravedad apunta hacia abajo, por lo que podemos decir que la presión no depende de
𝑧.
𝑝 ≠ 𝑝(𝑧)
Como 𝑣𝑧 = 𝑣𝑧 (𝑟), solo 𝜙𝑟𝑧 (la cantidad de movimiento 𝑧 en la dirección 𝑟 positiva) y
𝜙𝑧𝑧 (la cantidad de movimiento 𝑧 en dirección 𝑧 positiva) contribuye al balance de
momento.
Fuerza gravitatoria sobre

la capa
Figura 1: Este es el caparazón sobre el cual se hace el balance de la cantidad de

movimiento para el flujo en forma horizontal al anillo orientado.
Tasa de impulso 𝑧 en el capazón de en 𝑧: (2𝜋𝑟Δ𝑟)𝜙𝑧𝑧 |𝑧
Tasa de impulso 𝑧 fuera del caparazón en 𝑧 + 𝐿: (2𝜋𝑟Δ𝑟)𝜙𝑧𝑧 |𝑧+𝐿
Tasa de impulso 𝑧 en el capazón de en 𝑟: (2𝜋𝑟Δ𝐿)𝜙𝑟𝑧 |𝑟
Tasa de impulso 𝑧 fuera del caparazón en 𝑟 + Δ𝑟: [2𝜋(𝑟 + Δ𝑟)𝐿]𝜙𝑟𝑧 |𝑟+Δ𝑟
Componente de la fuerza gravitacional sobre el caparazón en dirección 𝑧: 0
Si suponemos un flujo constante, entonces el balance de cantidad de movimiento es
𝑇𝑎𝑠𝑎 𝑑𝑒 𝑚𝑜𝑣𝑖𝑚𝑖𝑒𝑛𝑡𝑜 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 − 𝑇𝑎𝑠𝑎 𝑑𝑒 𝑚𝑜𝑣𝑖𝑚𝑖𝑒𝑛𝑡𝑜 𝑑𝑒 𝑠𝑎𝑙𝑖𝑑𝑎
+ 𝐹𝑢𝑒𝑟𝑧𝑎 𝑑𝑒 𝑔𝑟𝑎𝑣𝑒𝑑𝑎𝑑 = 0
Considerando solo los componentes de 𝑧, tenemos
(2𝜋𝑟Δ𝑟)𝜙𝑧𝑧 |𝑧 − (2𝜋𝑟Δ𝑟)𝜙𝑧𝑧 |𝑧+𝐿 + 𝑟 − [2𝜋(𝑟 + Δ𝑟)𝐿]𝜙𝑟𝑧 |𝑟+Δ𝑟 = 0
Factorizamos el lado izquierdo
−2𝜋𝑟Δ𝑟(𝜙𝑧𝑧 |𝑧+𝐿 − 𝜙𝑧𝑧 |𝑧 ) − 2𝜋𝐿[(𝑟 + Δ𝑟)𝜙𝑟𝑧 |𝑟+Δ𝑟 − 𝑟𝜙𝑟𝑧 |𝑟 ] = 0
Dividimos ambos lados por 2𝜋Δ𝑟𝐿
Tomamos el límite como Δ𝑟 → 0

El segundo término es la definición de la primera derivada de 𝑟𝜙𝑟𝑧
Ahora sustituimos las expresiones por 𝜙𝑟𝑧 y 𝜙𝑧𝑧
Como 𝑝 y 𝑣𝑧 no dependen de 𝑧, los términos 𝑝 y 𝑝𝑣𝑧2 se cancelan
Entonces tenemos
De la Ley de la Viscosidad e Newton sabemos que 𝜏𝑟𝑧 = −𝜇(𝑑𝑣𝑧 ⁄𝑑𝑟), entonces
Integramos ambos lados de la ecuación diferencial con respecto a 𝑟
Dividimos ambos lados por −𝜇𝑟
Integramos ambos lados de la ecuación diferencial con respecto a 𝑟 una vez más
Aplicamos las condiciones de contorno ahora para determinar 𝐶1 y 𝐶2

Resolviendo el sistema de ecuaciones obtenemos:
Armando con las constantes de integración, se conoce la distribución de velocidades
Por lo tanto,
Parte b
El caudal másico 𝑤 es, suponiendo una densidad constante 𝜌,
El caudal volumétrico 𝑑𝑉/𝑑𝑡 es la velocidad media multiplicada por el área de la

sección transversal
La velocidad promedio se obtiene integrando la velocidad sobre el área por la que fluye
el fluido a través y luego dividiendo por esa área
Hacemos una sustitución para resolver la integral
El caudal másico se convierte en

Por lo tanto,
Parte c
La tensión viscosa 𝜏𝑟𝑧 representa físicamente la fuerza en la dirección 𝑧 en una unidad
de área perpendicular a la dirección 𝑟. Evaluando en 𝑟𝑧 en 𝑟 = 𝑘𝑅 y multiplicándolo
por la superficie del área del cilindro interior, obtenemos la fuerza viscosa que actúa a lo
largo de su longitud. El punto final a tener en cuenta es que debido a que el fluido está
actuando desde un radio mayor 𝑟 en el cilindro interior, que tiene un radio menor 𝑘𝑅,
colocamos un signo menos delante de 𝜏𝑟𝑧
De la Ley de la viscosidad de Newton, tenemos
Ahora podemos encontrar la fuerza

Por lo tanto,
El signo menos aquí tiene sentido porque la fuerza viscosa (de fricción) del fluido se
opone a la del movimiento del cilindro, que se mueve en la dirección 𝑧 positiva.
Parte d
Para hacer que el anillo sea una hendidura plana, dejamos que el radio del cilindro
interior tienda hacia el radio del cilindro exterior; es decir, 𝑘 es muy ligeramente menor
que 1.
Sustituimos esto en la fórmula de 𝐹𝑧
El desarrollo de la serie de Taylor, para ln(1 − 𝜀) es el siguiente
Sustituimos esta fórmula en 𝐹𝑧
Usa la división larga para obtener una serie para la fracción

Por lo tanto,
La tensión viscosa 𝜏𝑥𝑧 para flujo laminar en hendidura con una pared que se mueve a
velocidad 𝑣0 se obtuvo de
Si orientamos la rendija horizontalmente de modo que la presión y la gravedad no

influyan en el flujo del fluido como es el caso en este problema, entonces 𝜏𝑥𝑧 se
simplifica
Dirección de fuerza
gravitacional
Figura 2: Este es el flujo de Couette en una rendija horizontal con una pared móvil en
𝑥 = 𝐵.
La fuerza viscosa que actúa sobre la pared móvil es 𝐹𝑧 = 𝜏𝑥𝑧 |𝑥=𝐵 ∙ 𝑊𝐿. No hay signo
menos aquí porque el fluido que actúa sobre la pared tiene una coordenada 𝑥 más baja
que la pared en movimiento en 𝑥 = 𝐵.
Comparando esta fórmula con el resultado encuadrado, vemos que son equivalentes si
𝜀 = 2𝐵, la rendija ancha y 𝑊 = 2𝜋. 𝑊debería der una distancia, pero es solo un ángulo;
por lo tanto, la fracción en el resultado no es realmente una fórmula de rendija plana. La
serie restante entre paréntesis entonces debe ser una corrección de curvatura
Fórmula de rendija
Corrección de curvatura
plana

Problema 2B.7 - BIRD

Cargado por

Copyright:

Formatos disponibles

Problema 2B.7 - BIRD

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Problema 2B.7 - BIRD

Cargado por

Copyright:

Formatos disponibles

Problema 2B.

a) Encuentre la distribución de velocidades en la región anular estrecha.

Como resultado, las condiciones de contorno se describen de la siguiente manera:

Fuerza gravitatoria sobre

Figura 1: Este es el caparazón sobre el cual se hace el balance de la cantidad de

Tomamos el límite como Δ𝑟 → 0

Ahora sustituimos las expresiones por 𝜙𝑟𝑧 y 𝜙𝑧𝑧

Como 𝑝 y 𝑣𝑧 no dependen de 𝑧, los términos 𝑝 y 𝑝𝑣𝑧2 se cancelan

De la Ley de la Viscosidad e Newton sabemos que 𝜏𝑟𝑧 = −𝜇(𝑑𝑣𝑧 ⁄𝑑𝑟), entonces

Integramos ambos lados de la ecuación diferencial con respecto a 𝑟

Dividimos ambos lados por −𝜇𝑟

Aplicamos las condiciones de contorno ahora para determinar 𝐶1 y 𝐶2

Armando con las constantes de integración, se conoce la distribución de velocidades

El caudal volumétrico 𝑑𝑉/𝑑𝑡 es la velocidad media multiplicada por el área de la

Hacemos una sustitución para resolver la integral

El caudal másico se convierte en

De la Ley de la viscosidad de Newton, tenemos

Ahora podemos encontrar la fuerza

Sustituimos esto en la fórmula de 𝐹𝑧

El desarrollo de la serie de Taylor, para ln(1 − 𝜀) es el siguiente

Sustituimos esta fórmula en 𝐹𝑧

Usa la división larga para obtener una serie para la fracción

Si orientamos la rendija horizontalmente de modo que la presión y la gravedad no

También podría gustarte