1 IQ36A FENOMENOS DE TRANSPORTE, SEMESTRE 08-1 GUIA CAPITULO 6.2 Cap. 6.2: Ecuacin de Navier-Stokes.

PROBLEMA 6.2-1.- Flujo uniforme con superficie libre, a lo largo de un plano inclinado. Analizar el flujo gravitacional de un lquido en direccin x a lo largo de un plano inclinado, de ancho infinito, que forma un ngulo con la horizontal. Se supone que el lquido se mueve como una capa de altura uniforme a y est expuesto a la atmsfera en la superficie libre. a) Demostrar que las ecuaciones de Navier-Stokes se reducen a la nica expresin: gx + (d2vx/dy2) = 0 cuando se aplican las siguientes hiptesis: (1) Flujo laminar; (2) Velocidad slo en direccin x; (3) Fluido incompresible; (4) Flujo estacionario; (5) Superficie del lquido a presin atmosfrica. (Se recomienda usar una hoja con las ecuaciones completas y marcar cada trmino anulado con el nmero de la hiptesis que lo justifica, para estar as seguros de que se ha demostrado lo que se pide). b) Integrar la ecuacin para encontrar vx, usando las siguientes condiciones de borde: (i) Velocidad nula en la pared del plano inclinado. (ii) Esfuerzo tangencial nulo en la superficie de contacto con la atmsfera (esto implica despreciar una leve resistencia friccional presentada por el aire). Solucin: Solucin: a) Aplicacin de las hiptesis para anular trminos en las ecs. de Navier-Stokes: Hiptesis 1 (H1): Flujo laminar: La velocidad del fluido es paralela a la pared: vy = 0. H2: Velocidad slo en direccin x: Esto implica suponer que vz = 0. H3: Fluido incompresible: Se aplica la ecn. de continuidad v = 0. Tomando la forma desarrollada de este operador vectorial de la tabla, vemos que se reduce a: vx/x = 0 y, en consecuencia, tambin 2 2 vx/x = 0. H4: Flujo estacionario: /t = 0. H5: Superficie del lquido a presin atmosfrica. Esta hiptesis requiere algn anlisis: La ecn. (C), con las hiptesis anteriores y considerando que la aceleracin de gravedad no tiene componente en la direccin z, se reduce a: p/z = 0. Anlogamente, la ecn. (B) se reduce a: p/y = gy = constante. Como en la superficie, y = a, se tiene p = patm = constante, se deduce que p/x = 0 en todas partes.

2 H6: Ancho infinito: Esta hiptesis permite anular el trmino vx/z . Se obtiene as como resultado la ecuacin del enunciado: gx + (d2vx/dy2) = 0, en que la derivada parcial se transforma en derivada total porque vx slo depende de y. b) Condiciones de borde para la integracin: (i) (ii) Velocidad nula en la pared: y = 0: vx = 0. Esfuerzo tangencial nulo en la superficie libre: yx = 0 en y = a.
2 2

El primer subndice del esfuerzo tangencial designa la normal a la superficie en la cual acta (la superficie libre es un plano perpendicular al eje y). El segundo subndice designa la direccin del esfuerzo (que acta en la direccin negativa del eje x). La frmula se encuentra en la tabla del esfuerzo tangencial:
= = - ( vx v y + y x

xy

yx

Como vy = 0, la condicin (ii) se reduce a : vx/y = 0 en y = a. Integrando dos veces la ecuacin diferencial, se obtiene: 2 vx = - (gx y / 2 ) + c1 y + c2 en que c1 y c2 son las constantes de integracin. Mediante las dos condiciones de borde, se obtiene finalmente: vx = - (gx y2 / 2 ) + (gx a y / ) PROBLEMA 6.2-2.- Medicin de la viscosidad Un aparato para medir viscosidades de lquidos consta de dos vasos cilndricos concntricos. Se coloca el lquido en el espacio anular entre ambos cilindros. El vaso interior gira a velocidad angular conocida; se mide el torque Tw que resulta de la friccin en la pared de uno cualquiera de los vasos; del torque se calcula el esfuerzo tangencial; se aplican finalmente las ecs. de Navier-Stokes para obtener la viscosidad del lquido. Para este fin, aplicamos las siguientes hiptesis para simplificar las ecuaciones de Navier-Stokes en flujo laminar: flujo estacionario, eje z vertical hacia arriba, trayectorias concntricas en planos z = const., simetra axial.

3 Notar adems que una de las condiciones de borde es: en r = Ri, vi = Ri para todo z. Esto requiere postular que v es independiente de z (de lo contrario, no podramos encontrar una solucin particular simple). En consecuencia: v/z = 0 a) Con las hiptesis indicadas, obtener una relacin entre el torque en el eje de rotacin y la viscosidad del lquido. b) Aplicar la solucin analtica para calcular la viscosidad de un lquido cuando se mide un torque igual a 0,28 (N cm). Los datos del viscosmetro son: radio del rotor Ri = 2 cm; radio del vaso exterior fijo Re = 2,2 cm; altura del lquido H = 5,4 cm; velocidad de rotacin = 200 RPM. Solucin: a) En primer lugar, se establece una relacin entre el torque en el eje en la pared del cilindro interior (Twi) y el esfuerzo tangencial en la pared. En la pared interior (superficie normal al eje r) acta un esfuerzo tangencial en la direccin tangencial ; por lo tanto, la designacin general es r, y por actuar en la pared interior, es ri. Sobre un elemento de rea dAi, el elemento de fuerza es ri dAi y el elemento de torque es un vector en la direccin axial z, de magnitud igual a la fuerza elemental por el brazo Ri, esto es: dTwi = ri dAi Ri. El mismo torque acta sobre cada uno de los elementos dAi que componen la superficie total del cilindro interior Ai = 2 Ri H. Como todos estos vectores tienen la misma direccin axial (z), se suman para dar el torque total: Twi = ri Ai Ri (1) Si se repite el mismo anlisis para la pared del cilindro exterior de radio Re, se obtiene en forma anloga: Twe = re Ae Re. Como el sistema est en equilibrio, girando con velocidad angular constante, ambos torques deben ser iguales. Para resolver las ecuaciones de Navier-Stokes, se aplican las hiptesis del enunciado, que conducen a las siguientes condiciones matemticas: - flujo estacionario: /t = 0; - eje z vertical hacia arriba, esto es, las componentes de la aceleracin de gravedad son: gz = -g; gr = g = 0; - trayectorias concntricas en planos z = const., esto es: vr = vz = 0. - simetra axial, esto es: / = 0. - v independiente de z: v/z = 0 La ecuacin de Navier-Stokes (ecuacin E) se reduce a: d 1 d (rv ) = 0 dr r dr

(2)

4 Las condiciones de borde son: Para r = Ri: v = Ri, para todo z. Para r = Re: v = 0, para todo z. La integracin conduce a:
2 2 Ri Re r2 v = (3) 2 2 r Re Ri El esfuerzo tangencial est definido por: ( v / r ) 1 v r ) (en este caso: vr = 0) (4) r = r = - (r + r r Reemplazando v segn ec. (3), se encuentra para el esfuerzo tangencial en la pared interior y en la pared exterior: 2 2 2 ri = 2 Re / (Re - Ri ) 2 2 2 re = 2 Ri / (Re - Ri ) Reemplazando en ec.(1), se obtiene: 2 2 2 2 (5) Twi = Twe = 4 H Ri Re / (Re - Ri )

b) Reemplazando los datos, se calcula: = 0,342 Pa s.

PROBLEMA 6.2-3.- Flujo vertical descendente a lo largo de la pared exterior de un tubo vertical.
Un tubo cilndrico vertical, de radio exterior R, est lleno de lquido, que rebalsa por el extremo superior y cae como una pelcula de espesor uniforme e adherida a la pared exterior del tubo. Modelar slo este flujo exterior, como si el cilindro tuviera longitud infinita. Suponer flujo laminar, con velocidad slo en la direccin axial hacia abajo, estacionario, incompresible, con simetra axial. En la superficie exterior de la pelcula (que est en contacto con el aire), suponer que la presin es atmosfrica y que el esfuerzo tangencial es despreciable. A partir de la ecuacin de Navier-Stokes, obtener la velocidad axial vz en funcin de la posicin radial.

Solucin: Es ms fcil visualizar el flujo como un cilindro macizo vertical, de longitud infinita y radio R. Adherida a la pared de este cilindro, cae una pelcula de lquido de espesor uniforme e. En coordenadas cilndricas de eje vertical hacia arriba, la pelcula de lquido ocupa el espacio comprendido entre r = R y r = R+e.
Se aplica la ecuacin de Navier-Stokes para encontrar la funcin vz = vz(r), bajo las siguientes hiptesis: (1) Flujo laminar axial: vr = v = 0; (2) Flujo estacionario: /t = 0; (3) Flujo incompresible, lo que significa que se aplica la ecuacin de continuidad: v = 0, que, junto con

5 (1), implica: vz/z = 0; (4) Longitud infinita (sin efectos de extremos); (5) Simetra axial: / = 0; (6) Eje vertical: gr = g = 0; gz = -g. Como condiciones de borde, se tiene: (7) Presin atmosfrica en la superficie de la pelcula lquida (esto es, para r = R+e); (8) Esfuerzo tangencial nulo en la misma superficie; (9) Velocidad nula en la pared del cilindro (r = R). La ecn. (D) de Navier-Stokes se reduce a: p/r = 0. Con condicin (5), resulta: p = p(z). Como p es constante e igual a patm en r = R+e, se concluye que p = patm en todo el interior del flujo. Por tanto: p/z = 0. La ecn. (E) es: 0 = 0. En ecn. (F), de (3) y (5) resulta que: vz = vz(r) y la derivada parcial es derivada total dvz/dr. As, la ecn. (F) es: dv z d g + ( r )= 0 r dr dr con las condiciones de borde: (i) en r = R: vz = 0; (ii) en r = R + e: rz = 0.
De la tabla de componentes del tensor esfuerzo tangencial, con las condiciones de este problema, resulta que la condicin (ii) implica: dvz/dr = 0. Por integracin: vz =
1 [g (r2 R2) 2g (R+e)2 ln(r/R)] 4

PROBLEMA 6.2-4.- Flujo vertical descendente en el espacio anular entre dos cilindros verticales concntricos, cuando el cilindro interior gira en torno a su eje.
Considerar un eje cilndrico vertical de radio 2 cm que gira con velocidad angular 450 RPM dentro de un cilindro fijo de radio 2,2 cm. En el espacio anular entre ambos hay un lubricante de densidad 950 Kg/m3 y viscosidad 6 cp. El lubricante gira arrastrado por el eje; tambin tiene un desplazamiento axial hacia abajo, con gradiente de presin constante dp/dz = 3500 Pa/m. a) Aplicar la ecuacin de Navier-Stokes al movimiento del fluido en el espacio anular, suponiendo flujo laminar, estacionario, incompresible, con componente radial de la velocidad nula en todas partes, longitud infinita, simetra axial.

NOTAR que las ecuaciones para las componentes axial y tangencial resultan independientes entre s. Por lo tanto, se pueden integrar separadamente, obtenindose las componentes axial y tangencial de trayectorias en espiral concntricas con el eje.

b) Calcular la magnitud del vector velocidad para la posicin r = 2,08 cm. Solucin: a) Modelacin matemtica: Al introducir las hiptesis del enunciado, la ecn. (D) de Navier-Stokes en coordenadas cilndricas se reduce a :
2 v 1 p = r r r

que permite calcular la distribucin radial de la presin despus de haberse resuelto para v (pero esto no forma parte de la pregunta). La ec. (E) de Navier-Stokes con las condiciones de borde indicadas para v constituye un problema idntico al Problema 6.2-2 (el viscosmetro rotatorio) y puede utilizarse dicho resultado en forma directa (en todo caso, la solucin se encuentra por integracin de la ec. (E) con las condiciones de borde: v = 0 en Re, v = Ri en Ri). Notar que la hiptesis de longitud infinita lleva a postular que v/z es nulo en todas partes. El resultado del Problema 6.2-2 es: v =
2 r2 R i2 R e 2 r Re R i2

La ec. (F) de Navier-Stokes, mirada separadamente, equivale al problema de un flujo axial a lo largo del espacio anular entre dos cilindros concntricos fijos. Este problema se formula en forma anloga al del flujo laminar en una tubera cilndrica (resuelto en clase), pero ahora con las condiciones de borde: En r = Ri, vz = 0; En r = Re, vz = 0; El resultado por integracin es:
2 B (R e R i2 ) ln (r / R i ) B (r 2 R i2 ) + 4 4 ln (R i / R e ) 1 p en que B = + g (una constante) z

vz =

b) Con las expresiones anteriores para v y vz, que son perpendiculares entre s, se calcula v =

[v

+ v2 z

1/ 2

. Con los datos numricos, resulta: v = 1,395 m en el punto (r = 2,08 cm).

7 PROBLEMA 6.2-5.- Slido que cae dentro de un tubo vertical. Un tubo cilndrico vertical de radio Re contiene un lquido de densidad y viscosidad . En el interior de este tubo se deja caer un objeto cilndrico macizo, de densidad s y radio Ri. El cilindro interior cae por efecto de la fuerza de gravedad y experimenta friccin con el lquido que lo rodea. Cuando ambas fuerzas se equilibran, el cilindro contina cayendo con velocidad uniforme V. Modelar este fenmeno para la condicin estacionaria indicada, con las siguientes hiptesis: el tubo y el cilindro tienen largo infinito; el cilindro mantiene una posicin concntrica con el tubo; el movimiento del fluido es laminar y slo tiene componente axial. Determinar la velocidad uniforme V. Solucin: Seguir el siguiente esquema: a) Resolver la ecn. de Navier-Stokes para obtener la velocidad del fluido: v = v(r,,z) en coordenadas cilndricas, suponiendo longitud infinita para el objeto. b) Calcular el esfuerzo tangencial en la pared del objeto y, con esto, calcular la fuerza de resistencia en la superficie total del objeto. c) Calcular la velocidad de cada del objeto con la condicin fuerza de gravedad + fuerza de resistencia = 0. Etapa a): Con las hiptesis indicadas, el problema es similar a la segunda parte del Problema 1b-3, flujo axial a lo largo del espacio anular entre dos cilindros concntricos, con la nica diferencia de que ahora las condiciones de borde son: En r = Ri, vz = -V (en que V es la velocidad de cada del objeto); En r = Re, vz = 0;
dv dp La ecn. (F) es: B + d (r z ) = 0 en que B = 1 +g r dr dr dz
Con las condiciones de borde indicadas, se obtiene:

Br2 vz = + c1 ln r + c 2 en que las constantes de integracin son: 4

rz = -

dv z B r c1 =- + dr 2 r

La fuerza fricional total es: Fz = rz (en Ri) 2 Ri L, para una longitud L del objeto.
Etapa c) Igualando Fz con el peso del objeto, W = ( Ri2 L s g), se obtiene:

B R 2 R 2 s g R i i ln i + B (R 2 R 2 ) V= e i 2 Re 4 2
PROBLEMA 6.2-6.- Forma de la superficie libre en un lquido que gira con velocidad angular constante.

Un recipiente cilndrico vertical, de radio R, contiene lquido hasta el nivel H sobre el fondo. El recipiente gira con velocidad angular constante. Debido al efecto centrfugo del movimiento, el lquido tiende a tomar un nivel ms alto en el borde del recipiente y ms bajo en el centro (eventualmente, si es suficientemente grande, el lquido se va a derramar sobre el borde del recipiente).

Determinar la forma de la superficie libre del lquido y, en particular, la altura Hw que toma el lquido en contacto con la pared y la altura H0 del lquido en el eje del cilindro. Para este fin, aplicar la ecuacin de Navier-Stokes en la forma siguiente: a) Resolver la ec. de Navier-Stokes suponiendo que el movimiento ocurre en rgimen laminar, estacionario, con simetra axial. Se supone que el lquido se mueve en trayectorias circulares concntricas con el eje z. Para integrar la ecuacin y obtener una expresin para la velocidad tangencial v, postular como condiciones de borde: Para r = R: v = R, para todo z. Para r = 0: v = 0, para todo z. Estas condiciones de borde requieren postular que v es independiente de z (de lo contrario, no podramos encontrar una solucin particular simple). En consecuencia, v/z = 0, b) Conocida la velocidad tangencial, integrar para encontrar la presin p(r,z). Notando que la superficie libre del lquido es una superficie de presin constante (igual a la presin atmosfrica), se utiliza la condicin de borde: para z = H0, p = patm. La ecuacin de la superficie libre se deduce de la ecuacin de la presin , con la condicin p = const = patm. c) Finalmente, para encontrar Hw, se aplica la condicin: el volumen de lquido en reposo es igual al volumen de lquido en rotacin (integrar la ecuacin de la superficie libre para encontrar el volumen en rotacin).

9 Solucin: La superficie libre es una superficie que se encuentra a presin constante e igual a la presin atmosfrica. Por lo tanto, es necesario resolver la ecuacin de Navier-Stokes para la velocidad y la presin y, con esa solucin, buscar la forma de las superficies de presin constante. a) Solucin de la ec. de Navier-Stokes para la velocidad v: se aplican las siguientes hiptesis: H1: flujo estacionario: /t = 0; H2: trayectorias concntricas en planos z = const., esto es: vr = vz = 0. H3: simetra axial, esto es: / = 0. H4: fluido incompresible: Se aplica la ecn. de continuidad v = 0. H5: eje z vertical hacia arriba, esto es, las componentes de la aceleracin de gravedad son: gz = -g; gr = g = 0 Con esto, la ecuacin de Navier-Stokes (ecuacin E) se reduce a: 2 1 (rv ) + v = 0 r r r z 2

(1)

Adems, se indic en el enunciado la necesidad de postular v/z = 0, con lo cual la ec. (1) se reduce a: d 1 d (rv ) = 0 (2) dr r dr Con las condiciones de borde indicadas en el enunciado, la integracin conduce a: v = r (3) (este resultado implica que el lquido gira como si fuera un cuerpo slido).

b) Solucin de la ecuacin de Navier-Stokes para la presin: con la solucin encontrada en ec. (3) para v, se integran ahora las ecs. (D) y (F) para encontrar la presin p. Con las hiptesis indicadas, estas ecuaciones se reducen a:

2 v 1 p = r r

( D)

1 p 0= (F) g z La integracin de la ec. (D) lleva a: p = 2 r2 + f1(z) (donde f1(z) es una funcin arbitraria de z) Anlogamente, la ecn. (F) lleva a: p = - g z + f2(r) (donde f2(r) es una funcin arbitraria de r) De aqu resulta que la presin p est dada por: p = 2 r2 - g z + a (donde a es una constante arbitraria) La constante a se determina con la condicin de borde: Para r = 0: p = patm: z = H0. Se obtiene: a = patm + g H0. Con esto, la solucin es:

10 p = 2 r2 + g (H0 z) + patm (4)

c) Ecuacin de la superficie libre: Definiendo como (rs, zs) las coordenadas de la superficie libre, la ecuacin se obtiene haciendo p = patm en ec. (4): 2 rs2 + g (H0 zs) = 0 (5) Para determinar el valor de H0, se iguala el volumen ocupado por el lquido en estado de reposo (con altura H, dato) con el volumen bajo la superficie libre, obtenido por integracin: R2 H =

z s 2 r dr

(donde zs se reemplaza desde ec. (5) como funcin de r) (6)

El resultado es: H0 = H (2 R2 / 4 g)

Finalmente, se obtiene Hw mediante ec. (5), calculando el valor de zs para rs = R. Introduciendo en ec. (5) el valor de H0 obtenido en ec.(6), se llega a: Hw = H + (2 R2 / 4 g) (7)