CI41A Repaso Cont
CI41A Repaso Cont
CI41A Repaso Cont
ALDO TAMBURRINO
CI 41A - HIDRAULICA
ESCURRIMIENTO EN TUBERIAS
Consideremos una tubería de diámetro constante, rectilínea y horizontal. Como los filetes de
escurrimiento son paralelos, tenemos una distribución hidrostática de presiones. Si pB es la presión en un
punto ubicado en el eje de la tubería, la presión en los puntos A y B será:
pA = pB - Rγ
pC = pB + Rγ
2R
Consideremos ahora una tubería curvilínea. Debido a la curvatura de los filetes existe una fuerza
centrífuga que actúa sobre las partículas del fluido, haciendo que la distribución de presiones no sea la
hidrostática. Podemos usar esta variación de la presión para determinar el caudal que escurre por la tubería,
conocida la diferencia de presiones en las paredes externa e interna.
1 p1' = p1+(2R+l+h) γ
p2' = p2 + lγ
2R p1' = p2' + γmh
R2
l De donde resulta:
R1 2
(p1 - p2)/ γ = (γm/γ -1)h - 2R
2R h
1’
2’
1
UNIVERSIDAD DE CHILE CI 41A – HIDRAULICA Prof. ALDO TAMBURRINO
R1 v1 = R2 v2
p1 - p 2 2 2
v - v1
= 2 - 2R
γ 2g
El caudal está dado por Q = VA, donde V es la velocidad media de la sección. Para simplificar el
problema, consideremos V ≈ (v1 + v2)/2. Formando la suma por diferencia en la ecuación anterior y
reemplazando en ella la expresión para la diferencia de presiones encontrada anteriormente, obtenemos:
Utilizando R1 v1 = R2v2 , v1+v2 = 2V y V=Q/A podemos reescribir la ecuación anterior como:
γm v -v
- 1 h = 2 1 V
γ g
γm - 2 Q2
h = R 1 R 2
γ R 1 + R 2 g A2
g R1 + R2 γm
Q=A h
2 R1 - R2 γ
2
UNIVERSIDAD DE CHILE CI 41A – HIDRAULICA Prof. ALDO TAMBURRINO
∆x
D
FP1
τ0
τ0
(1) FP2
W
(2)
z1 α
x
z2
Datum
r r r r
∑F ext = FP + Froce + W
Fpx = Fp1 - F p2
Expandiendo Fp2 en series de Taylor y despreciando términos de orden superior, obtenemos la siguiente
expresión para la resultante de las fuerzas de presión según x:
∂ Fp ∂p
Fpx = - ∆x = - A ∆x
∂x ∂x
3
UNIVERSIDAD DE CHILE CI 41A – HIDRAULICA Prof. ALDO TAMBURRINO
donde A representa el área de la sección transversal de la tubería. La fuerza de roce está dada por:
Froce x = - τ 0 χ ∆x
W x = W sen α = γ A ∆x sen α
∂p
τ0 = (- + γ senα) R
∂x
B1 = B2 + Λf.
Λf dB
J= =-
∆x dx
∂p Λf
(- + γsen α) = γ
∂x ∆x
τ0 = γ R J
Notar que la expresión para t 0 es general. No se han hecho suposiciones de tipo de flujo o tipo de
4
UNIVERSIDAD DE CHILE CI 41A – HIDRAULICA Prof. ALDO TAMBURRINO
pared. Si limitamos el análisis para tuberías circulares, t 0 será una función de la velocidad media, V; del
diámetro de la tubería, D; del espesor de las asperezas, e; y las propiedades del pluido, µ y ρ. Notar que
hemos dejado de lado la presión. (Esto no es tan obvio. Euler gastó varios años tratando de ligar la fuerza
debido a la presión actuando sobre la pared con el esfuerzo de corte).
τ0 D
2
= φ(Re, )
ρV ε
o:
D
τ0 = ρ V 2 φ(Re, )
ε
de donde resulta
V2 D
J= φ(Re, )
gR ε
D 1 V2
J = 8 φ(Re, )
ε D 2g
El término 8 φ(Re, D/e) es nuestro conocido factor de fricción, el que denotamos como f, o sea:
1 V2
J=f
D 2g
La expresión anterior es la llamada fórmula de Darcy-Weisbach. f(Re, D/e) está dado, en forma general, en
el diagrama de Moody.
TAREA: Demostrar que la distribución del esfuerzo de corte en una tubería circular está dada por
τ = τ0 2r/D = τ0 (1 - 2y /D), donde y es la distancia desde la pared.
5
UNIVERSIDAD DE CHILE CI 41A – HIDRAULICA Prof. ALDO TAMBURRINO
Usando la ecuación anterior y la que liga t 0 con J es posible determinar f en forma analítica si se
conoce la distribución de velocidades. En efecto, reemplazando la expresión de Darcy-Weisbach en
τ 0 = γRJ y despejando f, obtenemos:
2
u
f =8 *2
V
Para el cálculo de V supongamos que la ley logarítmica de velocidades es válida en toda la región
del flujo y consideremos una tubería de pared hidrodinámicamente lisa.. Esto significa que estamos
despreciando la presencia de la subcapa viscosa. (TAREA: Estimar el error que se comete al despreciar la
subcapa viscosa). La ley logarítmica está dada por:
u 1 y u*
= ln( )+B
u* κ ν
D /2
1
V=
A ∫ u 2π rdr
0
donde A = πD2 /4 y se cumple la relación r+y = D/2, donde y es la distancia medida desde la pared.
Después de integrar resulta (revisar):
1 D u* 3
V = u * ln - + B
κ 2ν 2
6
UNIVERSIDAD DE CHILE CI 41A – HIDRAULICA Prof. ALDO TAMBURRINO
2
2
u 1
f =8 * = 8
V 1 ln Du* - 3
κ 2ν + B
2
Generalmente el dato que conocemos es V y no u*, por lo que conviene expresar la ecuación
anterior en términos de la velocidad media. Si también reemplazamos Re = VD/v, podemos obtener la
siguiente ecuación (TAREA):
1 1
= ( ln ( f Re ) + B κ - ln ( 32 ) - 1.5)
f 8 κ
1
= 0.884 ln( f Re ) - 0.913
f
Esta ecuación fue obtenida por Prandtl en 1935 y la expresó en términos del logaritmo decimal:
1
= 2.035 log ( f Re ) - 0.913
f
Debido a que Prandtl supuso que la ecuación de la capa límite interna de la región turbulenta era
válida en toda la región del flujo, la ecuación anterior tiene necesariamente que tener un error. Prandtl usó
los datos experimentales tomados por Nikuradse y ajustó las constantes numéricas, obteniendo:
1
= 2 log(Re f ) - 0.8
f
La ecuación anterior no tiene limitación en el número de Reynolds. Debemos recordar que en 1911, Blasius
propuso la siguiente relación empírica:
0,3164
f= 1/4
, Re ≤ 105
Re
7
UNIVERSIDAD DE CHILE CI 41A – HIDRAULICA Prof. ALDO TAMBURRINO
En forma similar puede obtenerse la ley de resistencia para régimen turbulento en pared
hidrodinámicamente rugosa. En este caso la ditribución de velocidades es:
u 1 y
= ln + B , ( κ , B ) = (0.4, 8.5)
u* κ k s
resultando:
1 D D
= 0,942 ln + 1.68 = 2,2 log
+ 1,68
f 2k s 2 k s
Al igual que en el caso de pared lisa, la ecuación anterior, no es completamente correcta. Al ajustar
las constantes numéricas con información experimental, resulta la siguiente expresión:
1 3,7D
= 2 log ( )
f ks
En este caso no tenemos una expresión analítica simple para la distribución de velocidades ya que la
constante aditiva depende de ks u* /v. Para paredes lisas teníamos:
1
2 log(Re f ) - = 0,8
f
3,7D 1
2 log( )- =0
ks f
Veamos cuál es la desviación de f para paredes lisas con respecto a f para paredes rugosas. Para
visualizar esto, manipulemos la ecuación de paredes lisas con el objeto de poder formar la de paredes
rugosas. Al hacerlo (TAREA), resulta:
8
UNIVERSIDAD DE CHILE CI 41A – HIDRAULICA Prof. ALDO TAMBURRINO
3,7D 1 k u
2 log( )- = 1,03 - 2 log s *
ks f ν
3,7 D 1
M = 2 log -
ks f
Insistamos que aunque la expresión anterior tiene la forma de la ecuación para tuberías rugosas, f
cubre todo el rango de rugosidades debido a la presencia de M. Usando los datos experimentales de
Nikuradse, podemos obtener la siguiente información:
Re D/ks f M ks u* /v
: : : : :
: : : : :
: : : : :
: : : : :
Nikuradse graficó los valores de M versus k s u* /v de la tabla anterior, obteniendo la curva indicada
en el gráfico siguiente:
3
2,5
2 Experiencias de Colebrook y White
1,5 (tuberías comerciales)
M 1
0,5
0
0,1 -0,5 1 10 100 1000
-1 Experiencias de Nikuradse
(tuberías con rugosidad uniforme artificial)
Ecuación de M para tuberías lisas u* ks/ ν
9
UNIVERSIDAD DE CHILE CI 41A – HIDRAULICA Prof. ALDO TAMBURRINO
Nikuradse realizó sus mediciones en tuberías con asperezas artificialmente añadidas, de tamaño
uniforme. Las tuberías comerciales no tienen una rugosidad uniforme, sino que fluctúa entre ciertos límites.
Para utilizar los resultados de Nikuradse en tuberías comerciales se procedió a la inversa, o sea, se midió la
caída de presión en tuberías comerciales (recordar ∂p/∂x →τ0 → f) y se determinó cuál sería la rugosidad
equivalente en una tubería de Nikuradse que diese la misma pérdida. La rugosidad equivalente se denota
con e. Así, para tuberías rugosas:
1 3,7D
= 2 log
f ε
3,29
M = 2 log 1 +
u* ε
ν
por lo que la ley friccional para tuberías en transición lisa-rugosa tiene la forma:
1 ε 2,51
= - 2 log +
f 3 ,7 D VD f
ν
La discrepancia entre los resultados de Nikuradse y Colebrook-White para M se debe a que las
tuberías comerciales al tener rugosidad no uniforme, presentan una distribución (en el sentido estadístico) de
rugosidades. Es así como algunas asperezas tienen un tamaño mayor que el promedio y ellas tienden a
perturbar la subcapa viscosa antes que si todas las asperezas tuviesen el mismo tamaño, haciendo que la
separación de la curva con respecto a la tubería lisa comience antes y gradualmente.
De este modo, hemos caracterizado con ecuaciones los factores de fricción para tuberías con
escurrimientos en régimen laminar y turbulento. No se ha hecho lo mismo para el régimen de transición
laminar-turbulento ya que es muy incierto y pequeñas variaciones de Re significan grandes cambios en f. En
general, se trata de evitar el diseñar en dicha zona.. Las expresiones aquí presentadas no son las únicas.
10
UNIVERSIDAD DE CHILE CI 41A – HIDRAULICA Prof. ALDO TAMBURRINO
Existe una enormidad de ellas, prácticamente todas determinadas experimentalmente, algunas de las cuales
introducen parámetros adicionales que dependen del material del que están hechas las tuberías.
MATERIAL ε
(mm)
Latón, cobre, plomo, vidrio 0,0015
Cemento-asbesto (rokalit) 0,0125
Acero sin remaches, fierro forjado 0,045
Fundición asfaltada 0,120
Fierro galvanizado 0,15
Dovelas de madera 0,18 – 0,90
Fundición 0,26
Concreto 0,30 – 3,0
Mortero seco 1,25
Acero remachado 0,9 – 9,0
Metal corrugado 20
11