UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE INGENIERÍA QUÍMICA
Departamento de Operaciones Unitarias
Fenómenos de Transporte
Transporte de Momento
Guía de Aprendizaje GA04-EC02
Taller
Integrantes:
Dayana Marcela Ruiz Flores
Ximena Gabriela Mejía García
Julio Ricardo García Cuestas
Docente: MPhil. Ing. Rodolfo Espinoza
Grupo: 4
3-T-1 IQ
24 de abril del 2021
2C.4 Viscosímetro de cilindro descendente. Un viscosímetro de cilindro
descendente consta de un largo recipiente cilíndrico vertical (de radio R), cerrado
en ambos extremos, con un pedazo de metal cilíndrico sólido (de radio kR). El
pedazo de metal está equipado con aletas, de modo de modo que su eje coincide
con el del tubo.
La velocidad de descenso del pedazo metálico en el recipiente cilíndrico puede
observarse cuando este se encuentra lleno de fluido. Encontrar una ecuación que
proporcione la viscosidad del fluido en términos de la viscosidad terminal v o del
pedazo de metal y las diversas cantidades geométricas que se muestran en la
figura.
Desarrollo:
Sistema en estudio con las respectivas vistas y acotaciones:
Consideraciones y Suposiciones tomadas en cuenta
Sistema en estudio Flujo Fluido
Cilíndrico Estado Fluido
Simetría de eje estacionario Newtoniano
axial (z) Flujo Fluido de μ
Coordenadas (z, desarrollado constante
r, θ) Flujo Laminar Fluido
0 ≤ θ ≤ 2π Flujo reptante Incompresible (ρ
kR ≤ r ≤ R Flujo constante)
0≤z≤H unidireccional (z) No deslizante
No hay efectos Flujo de Couette
de borde
No hay
perturbaciones
de entrada ni de
salida
Área transversal
al flujo constante.
Áreas paralelas
al flujo constante
Sistema de Coordenadas Cilíndrico (r, θ, z)
Condiciones de fronteras:
Para r=R y para todo z, θ
Para r=kR y para todo z, θ
vmáx = vo
Con 0<k<1
0 ≤ vz ≤ vmáx
Ecuaciones diferenciales correspondientes:
Sistema en estado estacionario __
Flujo unidireccional __
Flujo desarrollado __
Vz = f(r) Vz = f(θ, z) __
Ecuación de continuidad
Como el fluido es incompresible
6𝑣𝜃
= 0 “Flujo desarrollado”
6𝜃
Ecuación de Movimiento (Navier Stokes)
∂ v r v θ ∂ v r v θ2 2 2
∂ 1 ∂ ( r v r ) 1 ∂ vr 2 ∂ vθ ∂ v r
ρ(
∂ vr
∂t
+ vr +
∂ r r ∂θ
− + vZ
r
∂ v r −∂ P
∂z
)=
∂r
+µ
[ (
∂r r ∂r )+ 2 − +
r ∂θ 2 r 2 ∂ θ ∂ z 2 ]
+ ρ gr
2 2
∂ 1 ∂ (r vθ ) 1 ∂ vθ 2 ∂ vr ∂ vθ
ρ(
∂ vθ
∂t
+v r
∂ vθ vθ ∂ vθ vr vθ
+
∂r r ∂θ
−
r
+ vZ
∂ v θ −1 ∂ P
∂z
)= 2
r ∂θ
+µ
[ (
∂r r ∂r )
+ 2 + +
r ∂θ2 r 2 ∂ θ ∂ z 2 ]
+ ρ gθ
∂ v v ∂ v z v θ2 2 2
∂vz ∂v 1 ∂ (r v z) 1 ∂ v z ∂ v z
ρ ( ∂t
+ vr z + θ
∂ r r ∂θ
− + vZ z =
r ∂z
−∂ P
)∂z
+µ
[
r ∂r
+ 2
r ∂θ
2
+
∂z
2 ]
+ ρ gz
Sistema en estado estacionario: no se observan movimientos en la dirección axial
ni en la dirección radial.
Flujo unidireccional: (el fluido en coordenada z)
Flujo desarrollado: la velocidad Vz no está en función de z
Componentes de velocidad y dependencia
vr = f(r, 𝜃, z, t) sin ondulación
v𝜃 = f(r, 𝜃, z, t) sin rotación
vz = f(r, 𝜃, z, t)
Simetría de eje axial
Flujo desarrollado
Estado estacionario
Condiciones de Contorno
CC1.
Fluido en contacto con la pared del cilindro fluido interior
Fluido no deslizante
Cilindro interior desciende
Para r=kR y para todo z, θ Vz = vmáx
CC2.
Fluido en contacto con la pared interior del cilindro exterior.
Cilindro exterior inmóvil
Para r=R y para todo z, θ
Ecuación 1, Ec cont.=0
∂p
Ecuación 2, = pgr
∂r
1 ∂p
Ecuación 3, = pgθ
r ∂θ
⃗g=9.81 m/ s2
g z|=|⃗g|
|⃗
gr|∴|⃗
|⃗g|⊥|⃗ gr|=0
gθ|=0
|⃗
La fuerza de gravedad aplicada se presenta en la coordenada z, r y θ están perpendicular a la
gravedad.
∂v
Ecuación 4,
−∂ p
∂z
+μ[ ( )]
1 ∂
r z + pgz
r ∂r ∂r
∂v
∂p
∂z
− pgz=μ
1 ∂
r z
r ∂r ∂r [ ( )]
Presión modificada, p=P−Pgz z
Signo menos porque la coordenada z esta en dirección a la gravedad.
∂v ( p 0− pg z O ) −( P H −pg z H )
−μ [ ( )] [
1 ∂
r z =
r ∂r ∂ r H ]
∂v p −p
−μ [ ( )] ( )
1 ∂
r z = 0 h
r ∂r ∂ r H
∂ ∂v p −p
( )( )
r z = 0 h r
∂r ∂r μH
∂ ∂v p −p
∂ r ( ∂ r ) ( μH )
z 0 h
r = rdr
∂v
∫ ∂∂r (r ∂ r )=X ∫ rdr
z
∂ vz x r 2
r = +C 1, dividendo por r
∂r 2
∂ v z xr C 1
= +
∂r 2 r
( xr2 + Cr ) dr
∂ v z= 1
C
∫ ∂ v =∫ ( xr2 + r ) dr
z
1
xr 2
v z= + C1 ln ( r ) +C2
4
p 0 − ph r 2
Usando la sustitución x=v z y sabiendo que v z= ( μH )4
k 2 R2
v 0=x +C 1 ln ( kR ) +C 2
4
k2 R 2
C 2=v 0 −x −C1 ln ( kR )
4
x R2 x k 2 R2
0= +C 1 ln R + v 0−
( ) −C 1 ln ( kR )
4 4
x R2 1
0=
4
(1−k 2 )+ C1 ln
k
+ v0 ()
x R2 2
−v 0 − (1−k )
4
=C 1
1
ln
k ()
x R2 ( x R2
1−k 2 )
( () )
v0 + −v −
x R2 2 2 0
4
− ( 1−k 2 ) ln ( r ) + v 0− x k R − 4 2
ln ( kR )−(1−k )
4 1 4 1
ln
k () ln
k
x R2 (
v 0+ 1−k 2 )
x ( 2 2 2) 4
r −k R + ¿
4 1
ln()
k
x R2 (
v0 + 1−k 2 )
x 4 kR
v 0+ ( r 2−k 2 R2 ) +
4
ln
1
()
( ( )) ln
r
k
Velocidad máxima.
Reescribiendo la ecuación:
xR2 (1−k 2)
v +
xr 2 0 4 kR xk 2 R2
4
+
1
ln
r (( )
+v 0 −
4 )
ln
k ()
xR2 (1−k 2 )
v+
xkR2 0 4 kr xk 2 R2
4
+
1
ln
kr (( )
+v 0 −
4 )
ln
k ()
v max=v 0
Velocidad mínima:
xR2 (1−k 2)
v +
xr 2 0 4 kR xk 2 R2
4
+
1
ln
r (( )
+v 0 −
4 )
ln
k ()
2 2
xR2 v 0 + xR (1−k ) kr xk 2 R2
4
+
1
ln
kr (( )
−v 0−
4 )
ln
k()
v min=0
Velocidad promedio, ⟨ v ⟩
caudal
⟨ v ⟩=
areatransversal
Caudal
P0−P H
v= ( Hμ )
R 2π
0=∫ ∫ v z dθdr
kR 0
R 2π
P0−P H
( Hμ )∫∫ v dθdr
kR 0
z
0 ≤ θ ≤2 π
0≤ z≤ H
xR2 kR
R
x
∫ v 0 + 4 ( r 2−k 2 R2 ) +
4
(1−k 2) ln
v0 +
r
dr
( )
1
kR
ln
k ()
v 0 + xR (1−k 2 )
3 2
x r kR
v0 r +
4 3 ( 2 2
+k R r + ) 1
r ln
r ( ( ) )
+r
ln
k ()
R
x R 3
v 0 + xR2 (1−k 2) kR
(
∫ v 0 R+ 4 3 + k R + 2 3
) 1
R ln
r ( ( ) )
+R
kR
ln
k ()
2 2
x k 3 R3 3 3 v 0+ xR (1−k ) x
[
− v 0 kR+
4 3 (
+k R +
ln
1
k
)
( kR ln ( 1 ) +kR )
() ]
Área transversal
2π R
AT =∫ ∫ rdrdθ
0 kR
2π R
r2
AT =∫ ∫ dθ
0 2 kR
R 2 k R2 2 π
AT =θ
2
−
2 0
dθ |
AT =π ( R2−k 2 R2 )
Velocidad promedio ⟨ v ⟩
2 2
x k 3 R3 3 3 v 0+ xR ( 1−k )
v 0 kR+
4 3 ( +k R +
1 ) (kR)
⟨ v ⟩=
ln
k ()
π ( R2−k 2 R2 )
Flujo masico:
m= pv
2 2
x k 3 R 3 3 3 v 0+ xR ( 1−k )
[
m= v 0 kR+
4 3 (+k R +
ln
1
k
) (kR ) p
() ]
Esfuerzo de corte:
∂vz
τ z=−μ
∂r
xR2
τ z=−μ
[
xr
2
−
v0 +
4
ln
(1−k 2)
1
k ()
kR
r
Esfuerzo de corte máximo para r =kR
1
( ) ]
xR 2
−μ
xR
2
[−
v0 +
4
ln
(1−k 2 )
1
k ()
kR
kR
xR2
( ) ]
τ rzmax =−μ
τ
τ rz = rz =ξ
[ xkR
2
−
v0 +
4
ln
(1−k 2)
1
k () ]
τ rzmax
Gráfica velocidad vs radio
r
16
14
12
10
0
0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1
Gráfica esfuerzo de corte vs radio
