El método de Rabinowicz es un método que lo aplicamos para calcular
el flujo volumétrico realizando un cambio de variables para facilitar los
cálculos.
Para que se pueda aplicar se debe tener en cuenta que los ejes del
sistema en el que vamos a aplicarlo deben estar ubicados de manera
tal que el primer miembro de la integral (que Veremos más adelante)se
anule para que tenga validez la igualdad que se emplea.
Para el caso del ejercicio explicado en clase del fluido newtoniano que
circula por una cañería en sentido horizontal, por ejemplo,se ubican los
ejes de la siguiente manera;el eje z coincidente con el eje de la cañería
y el eje r coincidente con el radio
Recordando : para el sistema ya se ha determinado lo siguiente
que el perfil del esfuerzo de corte es :
p0 − pl
rz = .r y que
2L
El esfuerzo de corte en la pared es
𝑝 −𝑝
𝜏𝑤 = 02𝐿 𝑙 . 𝑅 = 𝜏𝑚𝑎𝑥
𝜏𝑟𝑧 𝑟
= y que la relación entre los mismos es ésta
𝜏𝑤 𝑅
(A)
Para calcular el caudal se puede usar esta expresión:
R R
dr 2
Q = vz .dA = vz .2 rdr = 2 vz .
A 0 0
2
Para este sistema , si el fluido es newtoniano esta resolución no
representa grandes dificultades, pero si el fluido no es newtoniano, se
complica la integral ,entonces se puede recurrir a un cambio de
variables , utilizaré como variable 𝜏𝑟𝑧 (despejando de A)
Recordemos antes lo que es una integral por partes
Si tengo el producto de dos funciones, y aplico las reglas de la derivada,
Tendré que d(f.g)=g.d(f) +f.d(g)….entonces
Si despejo de esa expresión g.d(f) será
g.d(f)=d(f.g)-f.d(g) y si aplico integral a ambos
miembros de la igualdad queda
∫ 𝑔. 𝑑 (𝑓) =ʃd(f.g)-ʃ f.d(g)
quiere decir que puedo reemplazar el primer miembro de la igualdad
por el segundo; pero como colocamos los ejes de manera tal que el
término ʃ d(f.g)=0
¿y porqué es cero? Porque esta integral sería la diferencia entre (f.g)en
el límite superior de la integral y (f.g) en el límite inferior , que en ambos
casos es cero,porque f en un caso y g en el otro son cero entonces
puedo reemplazar a la integral del producto de la función g por el
diferencial de f ; por la integral del diferencial de g por la función con el
signo cambiado
∫ 𝑔. 𝑑 (𝑓) =-ʃ f.d(g)
Aplicado a nuestro sistema
A = .r 2
porque dárea = 2 rdr y además
2
dr
rdr =
2
entonces si se aplica lo visto anteriormente y dado que se cumple que
Vz=0 para r=R
Vz=Vmax para r=0 entonces
𝑅 𝑅
𝑄 = 𝜋 ∫ 𝑣𝑧 𝑑𝑟 = 𝜋[𝑣𝑧 . 𝑟 2 ]𝑅0 − 𝜋 ∫ 𝑟 2 𝑑𝑣𝑧
2
0 0
R
Q = .0.R 2 − .vz max .0 − r 2 .dvz
0
𝑅
Se puede calcular 𝑄 = −𝜋 ∫0 𝑟 2 . 𝑑𝑣𝑧
Mi nueva variable es 𝜏𝑟𝑧
Necesito r=f(𝜏𝑟𝑧 ) r2 = f(𝜏𝑟𝑧 ) y 𝑑𝑣𝑧 = f(𝜏𝑟𝑧 )
Del recuadro del inicio se puede despejar
𝜏𝑟𝑧
𝑟= .𝑅 recordar que la variable es 𝜏𝑟𝑧
𝜏𝑤
𝜏𝑟𝑧 2 2
𝑟2 = ( ) .𝑅
𝜏𝑤
𝜕𝑣𝑧 𝜏𝑟𝑧
Para fluido newtoniano 𝜏𝑟𝑧 = −𝜇 ⇒ 𝑑𝑣𝑧 = . 𝑑𝑟 siendo
𝜕𝑟 −𝜇
𝜕𝜏𝑟𝑧
𝑑𝑟 = 𝑅
𝜏𝑤
Reemplazando en la expresión queda
𝑅 𝜏𝑤
2
𝜏𝑟𝑧 2 2 𝜏𝑟𝑧 𝜕𝜏𝑟𝑧
𝑄 = −𝜋 ∫ 𝑟 . 𝑑𝑣𝑧 = −𝜋 ∫ ( ) .𝑅 .( ). 𝑅
0 0 𝜏𝑤 −𝜇 𝜏𝑤
𝜏
𝜋𝑅3 𝜏𝑤 3 𝜋𝑅3 𝜏𝑟𝑧 4 𝑤
𝑄 = 3 ∫ 𝜏 𝑟𝑧 𝑑𝜏𝑟𝑧 = 3 [ ] ⇒
𝜇𝜏𝑤 0 𝜇𝜏𝑤 4 0
𝜋𝑅3
𝑄= 𝜏𝑤
4𝜇
Y se puede calcular el flujo volumétrico(comúnmente llamado
caudal)de esta manera.
Si quiero ahora expresarlo en función de la diferencia de presión
𝑝 −𝑝
𝜏𝑤 = 02𝐿 𝑙 . 𝑅 entonces quedará :
Q
=
( R )
3 4
max
=
R 4 ( p0 − pL )
rz
max
3
4 8 L
que es nuevamente la ecuación de Hagen-Pioseuille, válida para
régimen laminar,estado estacionario y fluido newtoniano,ydespreciando
los efectos de borde.
ACLARACIÓN:SE APLICA EN GENERAL ESTE MÉTODO CUANDO SIMPIFICA
LOS CÁLCULOS NO CUANDO LOS COMPLICA, ACÁ LO APLIQUÉ PARA
EJEMPLIFICAR SU USO