Guía para Práctica de Flujo Transitorio Pozo de Oscilación

Escuela de Ingeniería de Antioquia. Ingeniería civil.
Hidráulica. 2016-1. Profesor Juan Fernando Barros Martínez

Jueves 11 de febrero
Flujo Transitorio. Golpe de Ariete y Pozo de Oscilación

Fundamento teórico
Golpe de ariete en conductos por gravedad
El golpe de ariete es un fenómeno transitorio en el que la tubería no se considera rígida y el
líquido se trata como compresible. Este fenómeno se produce en los conductos al cerrar o abrir
una válvula, al poner en marcha o parar una máquina hidráulica, o al disminuir bruscamente el
caudal. Un flujo es transitorio cuando sus características hidráulicas, presión (p), velocidad (v) y
caudal (q), dependen del espacio y del tiempo. Si el flujo se trata como unidimensional, entonces
p=p(x, t), v=v(x, t), q=q(x, t).
La descripción del fenómeno puede consultarse en el anexo al final de este documento (que
corresponde a la sección 15 del libro Mecánica de Fluidos y Máquinas Hidráulicas de Claudio
Mataix).
Fórmulas de presión máxima o sobrepresión
La sobrepresión depende del tiempo de cierre de la válvula (tc). El cierre puede ser:
T T L (1)
Rápido, si t c < , =2
2 2 c
T (2)
Lento, si t c >
2
Donde T es el periodo de la onda de presión que recorre el tramo entre la válvula y el depósito, L
es la longitud del tramo de tubería y c es la velocidad de la onda en el agua.
Para cierres instantáneos (tc = 0), las fórmulas de Jowkosky expresan:
∆ P=ρcv (cierre total) (3)
'
Δ P=ρc (v−v ) (cierre parcial) (4)
donde ΔP es el aumento de presión, ρ la densidad, v la velocidad inicial, v’ la velocidad final y c es

la velocidad de la onda expresada como
√ Eo (5)
ρ
c=
√ 1+
Eo D
Eδ
donde:
c = velocidad de la onda en m/s
1
Eo = módulo volumétrico del fluido en N/m2

D = diámetro de la tubería en m
E = módulo elástico del material de la tubería en N/m2
δ= espesor de la tubería en m
El golpe de ariete aumenta directamente con:
La longitud de la tubería
La velocidad del flujo
La rapidez del cierre
Ecuación dinámica del pozo de oscilación

El tanque de oscilación es un dispositivo de alivio utilizado para reducir el efecto producido por el
golpe de ariete. Tomando un tramo dx de la tubería de conducción, puede verse que está
sometido a tres fuerzas en dirección x (figura 1):
Figura 1. Sistema de fuerzas en un elemento diferencial de la tubería.

El Peso del agua = dw sen α= γ Ac dx sen α
La fuerza debida a la diferencia de presiones = -dp Ac
La fuerza debida a la fricción = -γ dhfc Ac
donde dw es el peso de un elemento diferencial, α es el ángulo que el eje de la tubería hace con la
horizontal, Ac es el área de la tubería y dx, dp, dhfc son los diferenciales de longitud en dirección
de la tubería, presión y pérdidas por fricción a lo largo de la tubería.
Aplicando la Segunda ley de Newton:
2
γdxAc ∂ v (6)
γAc dx sen α−dp Ac−γd hfc Ac=
g ∂t
Reemplazando dx sen α = dh y simplificando,

dp 1 ∂v
dh− −d hfc − dx=0
γ g ∂t
En el sistema tanque-tubería-almenara (figura 2) no hay ondas de presión pues el golpe de ariete
se anula con la presencia del pozo de oscilación (almenara). Puede considerarse que líquido y
dv ∂ v dv
tubería son indeformables y por tanto =0 , entonces =
dx ∂ t dt
Figura 2. Sistema tanque-tubería-almenara.

Integrando la ecuación anterior según los límites que se señalan en la figura
H a−H 1 H a +z hfc lc
dp 1 dv
∫ dh− ∫ 2
−∫ d h fc−
γ 0 g dt ∫
dx=0
0 v 0
H 1−
2g
Haciendo hfc = cv2 y simplificando
v
2
2 lc dv (7)
z+ +c v + =0
2g g dt
1
Llamando K= +c ,
2g
3
2 lc dv (8)
z +k v + =0
g dt
y para tener en cuenta la dirección de la velocidad:
lc dv (9)
z +kv |v|+ =0
g dt
Esta expresión se conoce como la ecuación dinámica del pozo de oscilación. Para completar el
sistema de ecuaciones se acompaña la ecuación anterior con la de continuidad
dz (10)
vAc=Ap +Q
dt
donde Ap es el área transversal del pozo de oscilación y Q es el caudal en la tubería aguas abajo de
la almenara, siempre positivo o cero (en caso de cierre total).
Solución teórica de las ecuaciones del pozo de oscilación

Las ecuaciones del pozo de oscilación (9 y 10) pueden resolverse analíticamente sólo para el caso
de un cierre total instantáneo e ignorando la fricción, resultando para un período de oscilación
teórico
√ lc Ap (11)
T =2 π
g Ac
Método numérico de Scimemi

Expresando las ecuaciones 9 y 10 en forma de diferencias finitas
lc Δ v (12)
z +kv |v|+ =0
g Δt
Δz (13)
v Ac=Ap +Q
Δt
Se procede a calcular elementos finitos Δv y Δz para cada incremento Δt. Siendo i el número de
orden con el cual se mide el tiempo de funcionamiento del pozo en incrementos Δt, se tiene:
−g Δ t
Δ v i= (z i +k |v i−1| v i−1)
lc
Δt
Δ z i= (v Ac−Q)
Ap i
4
Haciendo
gΔt g Δt k Ac Δ t Q Δt
C 1= ,C 2= , C 3= , C 4=
lc lc Ap Ap
Δ v i=−ci z 1−c 2|v i−1| v i−1
Δ z i=c3 vi −c 4
Procedimiento de cálculo modificado

Las ecuaciones del pozo de oscilación presentadas hasta el momento, consideran un factor de
fricción (c) constante, independiente de la velocidad. Para tener en cuenta la variación de este
factor, se realizará un trabajo experimental en el sistema donde se llevará a cabo el golpe de
ariete. El trabajo consiste en registrar las pérdidas más la cabeza de velocidad desde el tanque de
alimentación hasta el pozo de oscilación para distintos caudales. Como las pérdidas se expresan en
la forma hf = cv2, puede obtenerse mediante regresión una ecuación relacionando z (es decir las
pérdidas más la cabeza de velocidad) y v, ya que z= (1/(2g) +c) v 2. Al hacer la regresión resultará
una ecuación z= c*vn de manera que k = c*v n-2. Con esta consideración el valor de C 2 se calcula
según la ecuación
g Δt K g Δt ¿ n −2
C 2= = ( c |v| )
lc lc
El procedimiento a seguir es el siguiente:
1. Se selecciona Δt (Δt < (1/40) T)
2. Se calcula C1, C3, C4
Con i = 0
g Δt ¿ n−2
3. Se calcula C 2=
lc
( c |v i| )
4. Se calcula Δ z i=C 3 v i−C 4
con i 1
5. Se calcula z i=z i−1+ Δ z i−1
6. Se calcula Δ v i=−C 1 z i−C 2|v i−1|v i−1
7. Se calcula v i=v i−1 + Δ v i
g Δt ¿ n−2
8. se calcula C 2=
lc
( c |v i| )
9. Se calcula Δ z i=C 3 v i−C 4
con i=i+1 se vuelve al paso 5
Práctica de laboratorio
5
1. Para diferentes caudales tomar lectura de h en el vertedero y Zo (entre el nivel del tanque y el
nivel en el pozo de oscilación).
2. Para un caudal determinado hacer el cierre rápido de la válvula y tomar los tiempos entre los
niveles máximos y mínimos en el pozo de oscilación marcando los niveles.
Resultados
1. Obtener por regresión la ecuación z= c*vn
2. En un diagrama de tiempo y altura (tiempo en la abscisa y Z en la ordenada), grafique la
oscilación del nivel en el pozo de oscilación y ubique los datos registrados en el experimento.
Precaución: tenga en cuenta que el valor inicial de z es negativo y corresponde a la diferencia
de nivel entre la superficie en el depósito y en el pozo de oscilación antes del cierre.
6
Anexo:
7
8
9
10

Guía para Práctica de Flujo Transitorio Pozo de Oscilación

Cargado por

Copyright:

Formatos disponibles

Guía para Práctica de Flujo Transitorio Pozo de Oscilación

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Guía para Práctica de Flujo Transitorio Pozo de Oscilación

Cargado por

Copyright:

Formatos disponibles

Escuela de Ingeniería de Antioquia. Ingeniería civil.

Hidráulica. 2016-1. Profesor Juan Fernando Barros Martínez

Flujo Transitorio. Golpe de Ariete y Pozo de Oscilación

∆ P=ρcv (cierre total) (3)

donde ΔP es el aumento de presión, ρ la densidad, v la velocidad inicial, v’ la velocidad final y c es

Eo = módulo volumétrico del fluido en N/m2

Ecuación dinámica del pozo de oscilación

Figura 1. Sistema de fuerzas en un elemento diferencial de la tubería.

Reemplazando dx sen α = dh y simplificando,

Figura 2. Sistema tanque-tubería-almenara.

Haciendo hfc = cv2 y simplificando

y para tener en cuenta la dirección de la velocidad:

Solución teórica de las ecuaciones del pozo de oscilación

Método numérico de Scimemi

Procedimiento de cálculo modificado

También podría gustarte

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.