10.
2 PROBLEMA 2
Considere una esfera de radio R que rota a una velocidad angular en el eje z. La esfera est
sumergida e un n un volumen de fluido grande. Las escalas de velocidad y longitud en este
problema son y R respectivamente, sugiriendo que la ecuacin de Stokes ser buena
aproximacin si:
2
= 1
Figura 3- 1.-. Distribucin del flujo reptante a travs de la superficie de una esfera rotante.
La forma ms simple del sistema es que la esfera rote en el eje z. Para lo cual se tiene que la
ecuacin de continuidad es satisfecha si:
= = 0
Las condiciones lmite para la velocidad son:
= 0, = 0, = , r=R
= 0, = 0, = 0, r=
[Link]
Para verificar que V satisface la ecuacin de Stokes. Con los componentes y r de la ecuacin de
Stokes se reducen a:
= = 0
Las ecuaciones de movimiento en coordenadas esfricas, al eliminar trminos quedan de la
siguiente manera:
2 + 2
( + + + )
1 2 1 1 2
= + [ 2 2 ( 2 ) + 2 ( )+ 2 ] +
2 2
+ 2
( + + + + )
1
=
1 1 1 1 2 2
+ [ 2 ( 2 )+ 2 ( ( )) + 2 +
2 2 2
2
] +
2
+
( + + + + )
1
=
1 1 1 1 2
+ [ 2 ( 2 )+ 2 ( ( )) + 2
2 2
2 2
+ + 2 ] +
2
1 2 1 1
[ ( )+ 2 ( ( ))] = 0
2
1 2 1 1 2
2 = ( ) + ( ) + ( )
2 2 2 2 2
1 2 1 1
[ ( )+ ( ( ))] = 0
2 2
Por lo tanto queda:
2
=
2 1
=
1 2 1 1
[ ( )+ 2 ( ( ))] = 0
2
Al realizar el anlisis geomtrico del sistema, se puede proponer una solucin en base a puntos de
simetra, pues como se est trabajando con flujo reptante, la funcin de corriente en la
superficie de la esfrica es simtrica (existe una lnea de simetra en el centro de la esfera).
Esta solucin propuesta est dada por:
= ()
=
= (2 )
2 = 2
1 1
( 2 ) 2 = 0
2 2
Ahora
2 + 2 2 = 0
2 + 2 2 = 0
Se propone la solucin de la siguiente forma:
Ecuacin
3 2
3 3
+ 2
2
+ + = 0
Solucin
= 1 1 + 2 2 + 3 3
Por lo tanto para la ecuacin 2.6 es:
= 1 1 + 2 2
Resolviendo
( 1) + 2 2 = 0
2 + 2 = 0
Si
1 1 + 8
=
2
Obtenemos que n1=1 y n2=-2 y si las sustituimos en la ecuacin 2.7 tenemos que:
= 1 + 2 2
Sabemos entonces que la solucin propuesta es:
= [1 + 2 2 ]
De la condicin de frontera obtenemos que:
= 0, = 0, = , r=R
Si = r = R
= [1 + 2 2 ]
Ahora r=
0 = [1 + 2 ]
Tenemos que C1=0 y C2=R3
Sustituyendo nos queda el perfil:
3
=
2
Para el estudio del torque ejercido por la esfera. Siendo R d cuando es con el radio y R sin d
cuando es con la proyeccin.
El torque para un elemento diferencial de superficie dS est dado por:
= (, ) sin
Donde:
= 2 sin
Al evaluar encontramos que el anlisis de la Fuerza es:
| 2 sin sin
Diferencia de Torque:
= | 2 sin sin
2
= | 2 sin sin
0 0 0
3 sin
= ( )
3
= 3 sin (3 2 )
= 3 3 sin /
2
= 3 3 sin3
0 0
= 3 3 2 sin3
0
1
= 6 3 ( cos + cos3 )|
3 0
1 1 4
= 6 3 (1 + 1 ) = 6 3
3 3 3
El torque ejercido por el fluido sobre la esfera en un elemento diferencial, para un punto sobre la
superficie est a una distancia R sin del eje de rotacin es:
3
= 8
Se presenta una esfera a travs de la cual circula un fluido con un flujo de tipo reptante. Sus
caractersticas se muestran en la siguiente figura:
z
(, , )
Figura 3- 2.- Distribucin del flujo reptante (funcin de corriente) a travs de la superficie de una esfera..
Por otra parte, algunas consideraciones importantes son las siguientes:
El nmero de Reynolds es tan bajo que la velocidad y en la superficie de la esfera son cero.
Adems, no existe vorticidad en ningn sitio cercano a la esfera, por tanto, no existe conveccin.
Lejos de la esfera, al infinito, se supone, la velocidad de la corriente es solamente ( ), as, y
son representadas de la siguiente manera:
Figura 3- 3.- Velocidad de y en el infinito como funciones de .
Se entiende entonces que y son una funcin de .
= cos y = sen
As, las condiciones lmite para este caso estarn dadas por:
r = R; = 0
r = R; = 0
donde R es la superficie de la esfera.
Y en el infinito (un punto lo bastante lejos de la esfera):
r = ; = cos
r = ; = -sen
Al realizar el anlisis geomtrico del sistema, se puede proponer una solucin en base a puntos de
simetra, pues como se est trabajando con flujo reptante, la funcin de corriente en la
superficie de la esfrica es simtrica (existe una lnea de simetra en el centro de la esfera), y con
mayor razn en un punto muy lejano a la misma esfera (infinito).
Esta solucin propuesta est dada por:
= 2
donde es la funcin de corriente y 2 es la parte simtrica de la solucin.
Por otra parte tambin depender de r a travs del rango establecido r = R a r = , por tanto
ser necesario agregar una funcin () que habr de encontrarse.
La solucin entonces estar dada por:
(, ) = ()2
() se obtendr a partir de la sustitucin de = ()2 en la solucin encontrada sobre la
ecuacin de Navier-Stokes para la misma funcin de corriente, Ec. (1.2).