INTRODUCCIÓN.

La viscosidad definida como la oposición del fluido a fluir, pues esta propiedad establece el
frotamiento interno entre las moléculas del fluido cuando se deslizan sobre otras; en la mayoría de
los fluidos tales como el agua, los esfuerzos a una velocidad de deformación es conoce como una
función lineal de la misma. Estos tipos de fluidos se conocen como fluidos Newtonianos, en donde
el coeficiente de viscosidad depende únicamente de la temperatura del fluido por lo que el esfuerzo
⃑ ⃑
( )
T
cortante (𝜏̿) se define de la siguiente manera: τ́ =−μ ( ∇ ⃑ ⃑
v + ∇ v ).
ij

En cambio, si hablamos de los fluidos no Newtonianos nos encontramos que estos no siguen la ley
de viscosidad de Newton, pues la relación entre el esfuerzo cortante y la velocidad de deformación
no es constante; por lo que el esfuerzo cortante tiene una variación de τ́ ij =−ŋ Δ́ .

Mediante un experimento casero, se moldea un viscosímetro de caída libre empleando canicas

donde se dejarán caer en dos diferentes fluidos; agua siendo un fluido Newtoniano y miel siendo
uno no Newtoniano. En el cual se hace un análisis al flujo alrededor de la canica mientras este
desciende por el tubo vertical por efecto de la gravedad, pues el objeto se acelera por la fuerza del
peso hasta que queda balanceada por la fuerza de flotación y arrastre viscoso.

Un parámetro muy importante que se debe tener en cuenta al momento de hacer el respectivo
análisis mencionado anteriormente es el número de Reynols (Re), el cual nos determina el régimen
de flujo.

ρ v prom D
ℜ=
μ
Para el estudio de la transferencia de cantidad de movimiento (análisis del flujo), es necesario
conocer y hacer uso de las siguientes ecuaciones.

ρD ⃗v
=−[ ⃗
∇ · ρ ⃗v ⃗v ]− ⃗
∇ p− [ ⃗
∇ · τ ] + ρ ⃗g ( 4 ) .
Dt
ρDv
Donde, es la velocidad de incremento de cantidad de movimiento por unidad de volumen en
Dt
otros términos acumulación, [ ⃗∇ · ρ ⃗v ⃗v ] es el transporte convectivo de cantidad de movimiento, ⃗
∇p
las pérdidas de presión, [ ∇
⃗ · τ ] el transporte conductivo por esfuerzos viscosos y ρ ⃗g es la fuerza
externa sobre el fluido, en este caso el peso.

Objetivos

 Aplicar los conceptos teóricos vistos en clase sobre transferencia de cantidad de

movimiento y reología para formular un modelo acorde con el sistema estudiado.
  Realizar el análisis del flujo alrededor de una esfera con las ecuaciones de balance de
cantidad de movimiento respectivas, determinando así, el perfil de velocidad para cada
fluido y su flux combinado.
 Determinar la expresión de la viscosidad de
los fluidos a partir del balance de fuerzas que actúan en el objeto (canica) mientras este
desciende por los fluidos.
 Analizar la influencia del diámetro del objeto (canica) sobre el flujo del sistema.

METODOLOGÍA

Se diseña una metodología apta para poder llevar a cabo con los objetivos planteados
anteriormente; dicha metodología está definida por tres fases:

1. Montaje y desarrollo del experimento y laboratorio de reología de viscosidad .(resultados)

2. Análisis de la ecuación de movimiento de Navier-Stokes, incluir flux combinado ¿Es esta
ecuación apta para poder definir el perfil velocidad del flujo alrededor de la esfera?
3. Análisis y desarrollo de las fuerzas que actúan en el descenso de la canica por el fluido.
4. Resultados(resultados- determinar viscosidad con velocidad promedio de las tablas con los
resultados), Análisis de resultados (% error de viscosidad, diferencia del diámetro de
esferas, estudio de velocidad limite seguimiento de la esfera, y coeficiente de arrastre y
volver hallar la viscosidad

CONCLUSIONES ( De acuerdo a los objetivos)

VIDEO

Conceptos, ecuaciones, objetivos, explicación del experimento, explicar perfil de velocidad,

resultado conclusiones.

1.0. MONTAJE Y DESARROLLO DEL VISCOSIMETRO

Para dicho experimento se utilizaron los siguientes materiales:

Dos recipientes cilíndricos Una esfera (canica) grande

rectos transparentes de vidrio pesada y una pequeña
 Miel
Agua

1.1. Montaje del experimento

Inicialmente se debe llevar a cabo el montaje descrito a continuación:

Para los dos fluidos correspondientes los cuales fueron el agua y la miel se debe utilizar un
recipiente cilíndrico recto transparente seco y a temperatura ambiente, y en él se debe verter
dichos fluidos.

Recipiente lleno de agua Recipiente lleno de

1.2. Recolección de datos

Se deja caer la esfera desde una altura de 16,5 cm y se toma el tiempo que tarda cada esfera en
llegar al fondo del recipiente a través de los diferentes fluidos. Este procedimiento se realizó 5
veces con cada uno de los fluidos y con cada una de las esferas, se adjunta evidencia:

Para el fluido 1 (Agua): Podemos observar el desplazamiento que la canica realiza a través del
agua hasta tocar el fondo del recipiente.

Para el fluido 2 (Miel): Podemos observar el desplazamiento que la canica realiza a través de la
miel hasta tocar el fondo del recipiente.

1.3. Tabulación de datos

Se tabularon y analizaron los datos correspondientes a cada uno de los fluidos mencionados
anteriormente teniendo en cuenta las cifras significativas para una mejor evaluación.
2.0. ANÁLISIS DE LA ECUACIÓN DE MOVIMIENTO DEL FLUJO REPTANTE.

Se lleva a cabo el algoritmo estudiado para determinar el perfil de velocidad de nuestro sistema,
para ello se hace una representación gráfica con sus respectivas consideraciones y condiciones de
frontera para así tener la realización de esta fase de la metodología:

CONSIDERACIONES

- Estado estable. - El tensor de esfuerzos cortantes es

- Esfera lisa. simétrico.
- Fluidos incompresibles. - Se desprecian efectos de borde, no
- Fluido Newtoniano (agua) y deslizamiento
Fluido no Newtoniano (miel).
- Flujo laminar (reptante).
-

CONDICIONES DE FRONTERA.

 En r =R ; v r=0
 En r =R ; vθ =0
 En r =∞ ; v r =v ∞ cos ( θ ) ; vθ =−v ∞ sin ( θ )
 Fig.1. Líneas de flujo, fronteras y ejes referencias (r, θ, ϕ ).

Consideremos el flujo muy lento de un fluido incompresible alrededor de una esfera sólida (flujo
reptante), La esfera es de radio R y diámetro D. El fluido tiene una viscosidad µ y una densidad ρ, y
asciende verticalmente hacia la esfera con una velocidad uniforme v ∞ a lo largo del eje z.
 Fig.2. Perfil y líneas de velocidad del fluido reptante alrededor de la esfera.

Para la toma del volumen de control no se tiene en cuenta los efectos finales y se supone que el
flujo es infinito indicar lo del y se toma de una parte intermedia del recipiente de vidrio, la velocidad
del fluido alrededor de la esfera implica de dos de los componentes que no desaparecen v r y vθ ,
estos quedan en función del uno y del otro v r(𝑟, 𝜃) y vθ (𝑟, 𝜃) y v ϕ=0 indica que no hay rotación de
la esfera. Por otra parte, se ve claramente que, lejos de la superficie esférica, la distribución de
presión se transforma en la ecuación hidrostática p= p 0−ρgz , siendo la componente z en relación
con las coordenadas esféricas igual a z=r cos ( θ ). Por lo tanto, las ecuaciones deben de satisfacer
las condiciones límites para r = R y r =∞ .[1]

CARACTERIZACION DE LOS FLUIDOS

kg
AGUA: Fluido newtoniano, ρ cte=¿ 998 [ ] −3
3 a temperatura ambiente, μ cte=¿1.002x10 ¿Pa*s]
m
[2] y se le aplica el modelo de la Ley de la viscosidad de Newton. (Cengel, Y., & Cimbala)
 dv
Donde, τ = Esfuerzo cortante aplicado, μ=Viscocidad dinámica del fluido, = Gradiente de
dx
velocidad y velocidad de deformación.

MIEL: Fluido no Newtoniano, ρ cte =? μ teorica=¿13.06¿Pa*s] enlazar con el estudio reológico y

se modela por medio de la Ley de potencias.

Donde, K es el índice de consistencia de flujo, y n es el índice de comportamiento de flujo. Cuando

n<1 el modelo predice un comportamiento pseudoplástico.

Se supone BORRAR TODO LO QUE HABIA

COMPONENTES DEL FLUX COMBINADO DE CANTIDAD DE MOVIMIENTO EXISTENTES EN

EL SISTEMA.

La ecuación del flux combinado de cantidad de movimiento:

ϕ ij =⃗
π ij + ρ ⃗v ⃗v
El flux molecular esta compuesto a su vez de:

π =⃗τ ij +¿ Pδ⃗ ij
⃗

Planteando la ecuacion general se representa de la siguiente manera:

π ij + ρ ⃗v ⃗v = τ⃗ij +¿ Pδ⃗ ij + ρ ⃗v ⃗v
ϕ ij =⃗

Analizamos la componentes del tensor de esfuerzo viscoso para mi condicion de velocidad, vθ (r , θ

) ; v r(r , θ) y v ϕ=0

[ ][ ]
τ rr τ rθ τ rϕ τ rr τ rθ 0
τ = τ θr τ θθ τ θϕ = τ θr τ θθ 0
τ ϕr τ ϕθ τ ϕϕ 0 0 0

Ley de viscosidad de newton (tomado de Apéndice B_Bird Stewart Lightfoot – 2ed. Educación B.1)
Resultado de los esfuerzos viscosos para mis condiciones de velocidad (tomado de Apéndice
B_Bird Stewart Lightfoot – 2ed. Educación B.1-18)

[
τ θθ=−μ 2(
1 ∂ vθ vr
r ∂θ r
+ ) ]
τ rr =−μ 2 [ ]
∂ vr
∂r

[ ( )
τ θr =τ rθ =−μ r
∂ ∂ vθ 1 ∂ vr
+
∂ r ∂ r r ∂θ ]
La igualación en esta expresión ( τ θr =τ rθ ) nos muestra la simetría que hay entre los dos tensores.

Procedemos a calcular el producto de Pδ⃗

[ ][ ]
1 0 0 P 0 0
Pδ=P 0 1 0 = 0 P 0
0 0 1 0 0 P

Dónde sabemos queδ⃗ es un vector identidad.

Ahora se procede a calcular el Flux molecular:

[ ][ ][ ]
τ rr τ rθ 0 P 0 0 P+ τ rr τ θr 0
π= τ θr τ θθ 0 + 0 P 0 = τ θr P+ τ θθ 0
0 0 0 0 0 P 0 0 P

Como último paso, se determina el flux convectivo:

[ ][ ]
vr vr vr vθ vr vϕ ρv r v r ρ vr vθ 0
ρvv=ρ v θ v r v θ v θ v θ v ϕ = ρvθ v r ρ vθ vθ 0
vϕ vr v ϕ vθ vϕ vϕ 0 0 0

Teniendo en cuenta el resultado del tensor de flujo molecular y el tensor flujo convectivo,
obtenemos entonces el flux combinado:

[ ][ ]
P+τ rr τ θr 0 ρv r v r ρ v r v θ 0
ϕ = τ θr P+τ θθ 0 + ρvθ v r ρ v θ v θ 0
0 0 P 0 0 0
Para obtener en la matriz el resultado final del flujo combinado
 [ ]
P+τ rr + ρvr v r τ θr + ρ v r v θ 0
ϕ = τ θr + ρvθ v r P+ τ θθ+ ρ v θ v θ 0
0 0 P

ϕ ij : tensor de flujo de cantidad de movimiento combinado

π ij: tensor de densidad de flujo de cantidad de movimiento molecular

⃗

ρ ⃗v ⃗v : tensor de densidad de flujo de cantidad de movimiento convectivo

τ⃗ij : Tensor de esfuerzos viscosos
P: Presion;

BALANCE DE CANTIDAD DE MOVIMIENTO DEL SISTEMA (NAVIER- STOKES)

Ecuación de movimiento para los fluidos agua y miel con μ y ρ constantes.[enlazar]

Componente en r:

ρ vr(∂ v r v θ ∂ v r v θ2 −∂ p
+
∂r r ∂θ
−
r
=
∂r )1 ∂2 2
r ∂r [
+ μ 2 2 ( r v r )+ 2
1 ∂
r sen θ ∂θ
∂v
(
sen θ r + ρ g r
∂θ )]
Componente en 𝜃:

(
ρ vr +
∂ r r ∂θ
+
r )
∂ v θ v θ ∂ v θ v r v θ −1 ∂ p
=
r ∂θ
+μ 2
r ∂r [ ( )
1 ∂ 2 ∂ vθ 1 ∂
r + 2
1 ∂
∂ r r ∂θ sen θ ∂ θ
( ( 2 ∂ vr
sen θ v θ ) + 2
r ∂θ )
+ ρ gθ
]
Ecuación de cantidad de movimiento en términos de τ .

Componente en r:

( ) [ ]
2
∂ vr vθ ∂ vr vθ −∂ p 1 ∂
2
1 ∂ τ
ρ vr + − = − 2 2 ( r 2 τ rr ) + 2 ( sen θ τ θr ) − θθ + ρ gr
∂r r ∂θ r ∂r r ∂r r sen θ ∂ θ r

Componente en 𝜃:

(
ρ vr +
∂ r r ∂θ
+
r
= )
∂ v θ v θ ∂ v θ v r v θ −1 ∂ p 1 ∂ 3
− 3
r ∂θ r ∂r
( [
r τ rθ ) +
1 ∂
rsenθ ∂ θ
( sen θ τ θθ) +
τ θr −τ rθ
r ]
+ ρ gθ

Ley de la viscosidad de Newton

 [ ( )
τ θr =τ rθ =−μ r
∂ ∂ vθ 1 ∂ vr
+
∂ r ∂ r r ∂θ ]
τ rr =−μ 2 [ ]
∂ vr
∂r

[
τ θθ=−μ 2(
1 ∂ vθ vr
r ∂θ r
+ ) ]
En los componentes r y 𝜃 de las ecuaciones de cantidad de movimiento no se descarta totalmente
el componente de convección, pero si el termino en función del tiempo ya que estos fluidos están
en estado estacionario; En los tensores de esfuerzos cortantes se omitieron los términos que
contienen (∇ ∙ v) por ser fluidos incompresibles.[1]

http://www.sc.ehu.es/sbweb/fisica/celeste/gravedad/gravedad.htm;

Para el fluido newtoniano (Agua), los perfiles de velocidad:

Como vemos, las ecuaciones de cantidad de movimiento resultantes para el componente r y 𝜃 se

encuentran en función de dos dimensiones (r, 𝜃) cada una, siendo más difícil obtener las
soluciones del perfil de velocidad. Para esto se introduce una función de corriente ψ en el flujo
asimétrico presente, que funciona como alternativa para hallar los perfiles de velocidad, teniendo
en cuenta de que esta se aplica para fluidos newtonianos, fluido en flujo laminar y numero de
Re<1.

La función de corriente establece que para los dos componentes de velocidad que no
desaparecen, se expresan como derivadas de ψ de tal forma que la ecuación de continuidad se
cumpla de manera automática. (ver tabla…) []. La ecuación diferencial que se establece para ψ , es
la ecuación D; Donde se utiliza el sistema de coordenadas esféricas en operaciones escalares.

La ecuación ψ para el flujo axisimétrico presente se asimila como:

( )
( )
2 2
∂ 2
( E ψ ) + 1 ∂ ψ , E ψ − 22 E ψ2 ∂ψ cos θ− 1 ∂ ψ Senθ
∂t rSen θ ∂ ( r , z ) r Sen θ ∂ r r ∂θ
4
E ψ =0
¿
Se establece las condiciones de frontera para la función de corriente correspondientes a las de
Navier Stokes dadas inicialmente.

Condiciones Limite

−1 ∂ ψ
C.L.1: en r=R, v r= =0
r sen θ ∂θ
2

+1 ∂ ψ
C.L.2: en r=R, vθ = =0
r sen θ ∂ r
Estas dos primeras condiciones limite describen la adherencia del fluido a la superficie de la esfera.

La tercera condición implica que v ∞ (velocidad limite o uniforme) lejos de la esfera (esto puede
verse al reconocer que v r=v ∞ cos ( θ ) ; vθ =−v ∞ sin ( θ ) lejos de la esfera)[1]

C.L.3: en r→ ∞,
−1 ∂ψ
v r= =v ∞ cos θ
r sinθ ∂ θ
2

δψ 2 2 2
=r v r sin θ=r ( v ∞ cos θ ) sin θ=r v ∞ ∈θ cos θ
∂θ
−1 2 2
ψ= v ∞ r Sen ( θ )
2
Sugiere esta condición de frontera que ψ tiene la forma
2
ψ (r , θ)=f (r ) Sen ( θ )
Se igualan los términos de la ecuación D (ver tabla…) con la solución propuesta descrita
anteriormente:
4
E ψ =0
¿
Para desarrollar la ecuación diferencial anterior, se aplica el método de separación de variables, ya
que el problema cede fácilmente a una solución de producto.

Solución del modelo en el primer operador E 2 de la ecuación D

 ( )
2
2 ∂ sen θ ∂ 1 ∂
E= + 2
∂r 2
r ∂ θ sen θ ∂θ

E2 (E 2 ψ )=0

[
∂2 sen θ ∂
∂r
2
+
1 ∂
r ∂ θ sen θ ∂θ
2 (
]([
∂ 2 sen θ ∂
∂r
2
+ 2 )
1 ∂
r ∂θ sen θ ∂ θ
2
] f (r )Sen ( θ ) )=0 ( )
Primero se determina

( E2 ψ )=0

( )
2
∂ sen θ ∂ 1 ∂ 2
[ + 2 ]f (r ) Sen ( θ ) =0
∂r 2
r ∂ θ sen θ ∂θ

[sen θ
2 ∂2
∂ r2
f ( r ) + f ( r )
sen θ ∂
r 2 ∂θ sen θ ∂ θ (
1 ∂ sen 2 θ
]=0 )
2
∂ f (r) 2
sen θ
[sen¿ ¿ 2θ 2
−2 2
f ( r )]=0 ¿
∂r r

[ ]
2 2
∂ 2
− 2 f ( r ) sen θ=0→ Segundo operador ( E ψ )=0
2
2
∂r r
Sustituyendo el cuadrado
2 2
E (E ψ )=0

[ ][ ]
2
∂ 2 2 ∂2 2
2
− 2 2
− 2 f ( r ) sen θ=0
∂r r ∂r r

∂ ∂2 f ( r )
Se utiliza una función de ensayo de la forma 𝑓(𝑟) = r n ; 𝑓(𝑟) = nr n−1 ; =n(n−1)r n−2
∂r ∂r
2

para E 4 ψ =0

[ n ( n−1 ) r n−2−
2 n
r 2 ][ 2
r n ( n−1 ) r n−2− 2 r n sen2 θ=0
r ]
¿
¿
Se encuentra que n puede tener los valores -1,1,2 y 4. Por consiguiente, f ( r ) tiene la forma

A
f ( r )= + Br +C r 2 + D r 4 =ψ (r ,θ)
r
 −1
Para cumplir la tercera condición límite, D debe ser cero, y C debe ser v . Las dos primeras
2 ∞
1 3 −3
condiciones límite dan A= V ∞ R y B= V R , entonces la distribución de la velocidad es
4 4 ∞
A 1 2
ψ (r , θ)= + B r − v∞ . r
r 2

[ ( ) ] cos θ
3
3R 1 R
v r=v ∞ 1− +
2 r 2 r

[ ( ) ] sinθ
3
3R 1 R
vθ =−v ∞ 1− −
4 r 4 r

√[ ( )] [ ( ) ] sin θ
3 3
3R 1 R 3R 1 R
v=√ v r + v θ =v ∞
2 2
1− + cos θ− 1− −
2 r 2 r 4 r 4 r
https://slideplayer.com/slide/17916995/

https://slidetodoc.com/analysis-of-creeping-flows-p-m-v-subbarao/

Para el fluido no Newtoniano (miel), el perfil de velocidad:

La miel al ser un fluido no newtoniano de tipo pseudoplastico, también llamado reofluidizante

(disminución de la viscosidad a medida que aumenta la deformación o esfuerzo cortante, según los
resultados del estudio reológico de la figura 232356; este es característico por tener un numero de
Re<1 y una viscosidad aparente, donde se le atribuye por su variación. Al determinar el perfil de
velocidad se aplica el modelo de Ley de potencias sobre las mismas consideraciones que al agua.

Modelo de la ley de potencias:

n−1
τ =−K |γ̇ | ( γ̇ )
Donde γ̇ es el tensor de velocidad de la deformación de la ley de la viscosidad de Newton
analizados inicialmente (ver ecuación 265,6). [1]

| ( ) |( ( ) )
n−1
∂ ∂ vθ 1 ∂ vr ∂ ∂ vθ 1 ∂ vr
τ θr =τ rθ =−K r + r + =−K ¿
∂r ∂r r ∂θ ∂r ∂ r r ∂θ

| | ( )
n−1
∂v ∂ vr
τ rr =−K 2 r 2 =−K ¿
∂r ∂r

| | ( )
n−1
1 ∂ vθ v r 1 ∂ vθ vr
τ θθ =−K 2( + ) 2( + ) =−K ¿
r ∂θ r r ∂θ r
Luego se reemplaza los tensores de esfuerzos cortantes del modelo de la ley de potencias
correspondientes, en cada componente de las ecuaciones de cantidad de movimiento de Navier-
Stokes en términos de τ .(ver ecuaciones 1,2)

Componente en r:

( )
2
∂ vr vθ ∂ vr vθ
ρ vr + − =¿
∂r r ∂θ r
−∂ p
−¿
∂r
−∂ p
−¿
∂r

Componente en 𝜃:

( vr
∂ vθ vθ ∂ vθ vr vθ
+
∂r r ∂θ
+
r
=¿ )
−1 ∂ p
−¿
r ∂θ
Se deja expresada las ecuaciones diferenciales para hallar los perfiles de velocidad en el fluido no
Newtoniano, ya que igualmente en la solución del fluido Newtoniano La solución matemática es
muy compleja.

4.0. FUERZAS QUE ACTUAN EN EL DESCENSO DE LA CANICA EN EL FLUIDO

Distribución de la fuerza de la presión:

Luego de obtener las distribuciones de la velocidad v r, vθ , se procede a calcular la presión que

actúa sobre la superficie de la esfera, los perfiles de velocidad se sustituyen en las componentes
resultantes de la ecuación de cantidad de movimiento Navier- Stokes (ver ecucion …) siendo
Ƥ=( p+ ρgz )
Componente en r:

( [ ( )] [ ( )] [ ( )] ( [
3R 1 R 3
3 R 1 R 3
3R 1 R 3
3R 1
∂ v ∞ 1− + cos θ −v ∞ 1− − sin θ ∂ v ∞ 1− + cos θ −v ∞ 1− −
2 r 2 r 4 r 4 r 2 r 2 r 4 r 4
ρ vr + −
∂r r ∂θ r

[ ( [ ( )]
)]
3
3R 1 R
∂ v ∞ 1− +

( [ ( )])
2 r 2 r
2 3
−∂ p 1 ∂ 3R 1 R 1 ∂
+ μ 2 2 r 2 (v ∞ 1− + )+ 2 sen θ + ρ gr
∂r r ∂r 2 r 2 r r sen θ ∂θ ∂θ
 ( )( )
μ v∞ 3
R
∫ ∂ Ƥ=∫ 3 R
2
r
3
cos θ ∂r

( )( )
−3 μ v ∞
2
R
P= 2
cos θ+ C r
2 R r
Se realiza el mismo procedimiento para el componente θ , y se aplica la C.L.3 que r→ ∞, donde la P
tiende a p0 siendo C θ=C r =p 0

( )( )
−3 μ v ∞
2
R
P= 2
cos θ+ C θ
2 R r

( R (r )
−3 μ v
)
2
R ∞
p+ ρgz= cos θ+¿ p ¿
2 0
2

( )( )
3 μv∞
2
R
p= p 0−ρgz− 2
cos θ
2 R r
La distribución de la presión resultante indica que hay una presión inicial que es la del plano z=0 y
(−ρgz ) el efecto del peso del fluido, donde se considera que z es la altura en el eje z alejado de la
esfera, siento esta valida para los flujos reptantes con numero de Reynolds inferiores a 1.

Balance de fuerzas que actúan en el sistema.

Cuando la esfera desciende verticalmente en el fluido en reposo, adquiere un movimiento

acelerado hasta alcanzar una velocidad de caída constante en condición de estado estacionario: la
velocidad limite o uniforme ( v ∞ ¿ , donde la suma de todas las fuerzas sobre la esfera debe dar cero
[1]. En la figura. Se representa el diagrama de fuerzas.

Balance de fuerzas sobre la esfera → ∑ F=Fe+ Fr −W =0

Donde, la fuerza de empuje o de flotación (principio de Arquímedes) es Fe y la fuerza de
resistencia Fr (Ley de Stokes), son el resultado de la fuerza total Ff que actúa sobre la esfera.
Siendo la fuerza normal Fn, la que actúa perpendicular a la superficie de ella y la fuerza tangencial
Ft o también llamada resistencia de fricción es donde el esfuerzo cortante actúa tangencialmente a
la esfera [1].

A continuación, se determina la fuerza total ejercida por el fluido

Ff =Fn+ Ft
Debido a la simetría alrededor del eje z, la fuerza resultante estará en la dirección z. Por
consiguiente, la fuerza se obtiene integrando las componentes z de las fuerzas normal y tangencial
sobre la superficie de la esfera donde los limites son 0 ≤ θ ≤ π , 0 ≤ ϕ ≤ 2 π ; en función del perfil de
presión p y el tensor de esfuerzos cortantes τ rθ según corresponda.

Expresión de la fuerza normal:

2π π
Fn=∫ ∫ ¿ ¿ ¿
0 0

2π π
Fn=∫ ∫ ¿ ¿ ¿
0 0

4 3
Fn=Fe+resistencia de forma → Fn= π R ρg+ 2 πμ R v ∞
3
Expresión de la fuerza tangencial:
2π π
Ft =∫ ∫ ¿¿ ¿
0 0

siendo τ rθ =−μ r
[ ( )
∂ ∂ vθ 1 ∂ vr 3 μ v∞
+
∂r ∂ r r ∂ θ
=
2 R ] ( )
senθ

Ft =Resistencia de fricción→ Ft=4 πμ R v ∞

Entonces la expresión resultante para la fuerza total del fluido sobre la esfera (Ver ecuación), es:

4 3
Ff = π R ρg+2 πμ R v ∞ +4 πμ R v ∞
3
Ff =Fe+resistencia de forma+ Resistencia de fricción
4 3
Ff = π R ρg+6 πμ R v ∞ → Ff =Fe + Fr
3
A partir del balance de fuerzas del sistema (ver ecuación ), se determina la relación para obtener la
viscosidad de los fluidos al medir la velocidad limite.

4 4
∑ F= 3 π R3 ρg+ 6 πμ R v ∞− 3 π R 3 ρ s g=0
Donde ρ es la densidad del fluido y ρ s la densidad de la esfera sólida, la ecuación de la viscosidad
según la ley de Stokes, y velocidad limite corresponden:
2 2
2 R ( ρs− ρ)g 2 R ( ρs− ρ) g
μ= → v∞=
9 v∞ 9 μ
Estas ecuaciones dependen de que el número de Reynolds sea menor a 1, también de que
velocidad limite sea la velocidad constante desde una cierta distancia descendiendo, donde las
fuerzas que actúan sobre esta se equilibren.

es que el fluido sea infinito, pero este no es el caso ya que está limitado por las paredes

del recipiente cilíndrico.

Por ende, para limitar la expresión 3.2-1 a la condición de fluido finito y poder medir la

velocidad media transformándola a velocidad limite, se utiliza la corrección de Ladenburg

que establece:

𝑣𝑙𝑖𝑚 = (1 + 2,4

𝑅𝑒𝑠𝑓𝑒𝑟𝑎

𝑅𝑟𝑒𝑐𝑖𝑝𝑖𝑒𝑛𝑡𝑒

) ∗ 𝑣𝑚𝑒𝑑𝑖𝑎 (3.2-2)

Otra manera de determinar la 𝑣𝑙𝑖𝑚, es realizando un seguimiento de la velocidad a través del

tiempo, hasta que esta se estabilice. (Tabla 5)