Flujo No Permanente - 07

FLUJO NO
PERMANENTE
Autor: Ing. William Chauca Nolasco

ECUACIONES DIFERENCIALES PARA
EL GOLPE DE ARIETE
Vamos a desarrollar las ecuaciones aplicando a
un segmento pequeo de fluido en una tubera.
Estas ecuaciones son las ecuaciones de cantidad

de movimiento y la ecuacin de continuidad .
Para ello las variables dependientes son las

alturas piezometricas H por encima de un eje de
referencia y la velocidad media v en una seccin
recta y las variables independientes son las
distancias x a lo largo de la tubera y el tiempo
t.

Haciendo un esquema para el anlisis de la cantidad de

movimiento, tenemos:
Fig. 01
D.C.L para la deduccin de la ecuacin del

movimiento
Para la deduccin no se tiene en cuenta la
relacin de Poisson. Para tuberas con juntas de
expansin estos no se consideran en la
deduccin.
Se considera que el rozamiento es proporcional

al cuadrado de la velocidad
Ecuacion de movimiento: Elemento de fluido

entre dos planos separados dx normales al eje
de la tubera. Se toma como cuerpo libre para la
aplicacin del segundo principio de Newton del
movimiento en la direccin axial (Figura 01)

Aplicando la ecuacin de movimiento para la grafica.
Desarrollando (1)

Desarrollando este termino
A esta ecuacin lo dividimos por una masa elemental de fluido

Sabiendo que:
Introduciendo la lnea de alturas piezometricas

Donde:
Luego desarrollando la ecuacin (2) teniendo en cuenta (3),

(4) y (*).
MULTIPLICANDO la ecuacin (5) POR -1
Por condiciones de flujo turbulento permanente.

Por hiptesis el coeficiente de rozamiento es el mismo que en flujo
permanente. Por lo tanto en (6).
Como el rozamiento se debe oponer al movimiento entonces
)
ECUACION DE CONTINUIDAD
La ecuacin de continuidad para flujo no permanente

La ecuacin de continuidad no permanente se aplica al
volumen de control de la figura anterior
no es una funcin de t desarrollando (9)
Dividiendo por la masa
Entonces:

Reordenando:
Agrupando
I II
I y II son las derivadas totales donde (I) se
refiere a la elasticidad de la tubera y a su
velocidad de deformacin con la presin. El
segundo termino (II) tiene en cuenta la
compresibilidad del liquido para la elasticidad
de la pared.
Fuerza de traccin en la
pared de una tubera
Podemos determinar lo siguiente:
a) La velocidad de variacin de las fuerzas de traccin por unidad de
longitud.
Si a la expresin (a) lo dividimos por t entonces.
b) Se determina la velocidad de variacin de la tensin.
Si a la (b) expresin lo dividimos por el modulo de elasticidad de Young

del material de la pared tenemos:
c) La velocidad de aumento de la deformacin unitaria.
Si a (c) lo multiplicamos por el radio obtenemos:
d) Velocidad de alargamiento radial:
Si a (d) lo multiplicamos por el permetro de la tubera nos da:
e) La velocidad de aumento de rea.
Por definicin del modulo volumtrico de elasticidad del flujo (II)
La velocidad de variacin de la densidad dividida por la densidad:

Utilizando la ecuacin (e1) y (12) entonces la ecuacin (11)
Se transforma en:
Factorizando:

Expresando la constante k con la relacin.

.(13)
Donde para las tuberas con juntas de expansin.
(13) en (12):
Como:

Desarrollando:

De esta expresin la variacin de la respecto a x o t es mucho menos

que la variacin de H respecto a x o t, por tanto se ha considerado a
contante.
Si la tubera estuviera en reposo y .
Entonces:
Donde la ecuacin final ser:
Donde encontramos la segunda ecuacin caracterstica denominada .
La ecuacin es la ecuacin de la continuidad para un lquido

compresible en una tubera elstica.
MTODO DE LAS CARACTERSTICAS
Lasecuaciones encontradas contienen 2 incgnitas . Se

pueden combinar con un multiplicador desconocido tal
como . Por lo que se debe tener un par de valores
cualquiera real y distinto. Entonces si realizamos un
agrupamiento preliminar de las ecuaciones L1 ^ L2.
Desarrollando la expresin y agrupando en forma conveniente.
(a) (b)
Este agrupamiento se hace con el fin de

obtener la derivada total
El primer trmino (a):

El segundo trmino (b):

Sabemos que:
Del clculo de derivadas totales podemos decir que (17) y (18) son
iguales:

Entonces
De (20) se ha hallado 2 valores reales de

que convierten las 2
ecuaciones diferenciales en derivadas parciales aun par de
ecuaciones diferenciales ordinarias.

Las 4 ecuaciones se pueden interpretar grficamente:
Se supone que se conoce V y H en R y S.

Es una curva de las ecuaciones (21) y (22) y solo es vlida a lo largo
de la caracterstica .La curva es un diagrama de las ecuaciones (23) y
(24) y solo es valida a lo largo de la curva .
Entonces cada ecuacin (21) y (23) contiene 2 incgnitas para un

punto conocido sobre su caracterstica.
En la interseccin (P) de y ambas ecuaciones se pueden resolver

para dar con Vp y Hp. En este punto las ecuaciones (22) y (24)
tambin se pueden resolver para dar x, t.
La solucin de obtiene a lo largo de la caracterstica partiendo de
condiciones inciales conocidas y hallando nuevas intersecciones de
manera que se hallan alturas y velocidades para tiempos posteriores.
Para los clculos de golpe de ariete en tuberas metlicas, V es muy

pequea comparado con a y se puede despreciar en las
ecuaciones (22) y (24)
Son rectas dichas caractersticas con pendiente
Se debe considerar que la tubera est construida en N tramos

iguales.
Red rectangular para los soluciones de las
ecuaciones caractersticas.
Se debe conocer H y V inicialmente en cada una de las secciones en

que se ha dividido la tubera.
Tambin deberamos tener en cuenta que la solucin se obtiene solo

sobre una regin limitada a menos que se de algn dato x=0, x=L de
alguna condicin extrema (c.c) como funcin del tiempo.
De la malla el intervalo de tiempo en
utilizando la retcula dado que x y t se conocen en cada interseccin y

no
.
consideramos la ecuacin
Escribiendo las ecuaciones (21) y (23) en notacin de diferencias

finita para la ecuacin (21).
Notacin en diferencias finitas.
En forma anloga para (23)

Sumando (25) y (26) se elimina Vp(I) y obtenemos HP(I)
De la ecuacin (25) y (26) al restar obtendremos Vp (I). (25) - (26):
Las Ecuaciones (27) y (28) se emplean en un circuito continuo para despejar Vp y Hp para
I=1 a N-1. Se deber utilizar las condiciones de contorno para despejar Hp(0) y Vp (0) as
como HP(N) y VP(N). Obtenidas estas , se pueden escribir los resultados, y despues se
sustituyen los valores de V(I) y H(I) por los de VP(I) y HP(I), se incrementa el tiempo y se
repite el proceso.
Condiciones de contorno
En el extremo aguas arriba o de izquierda la ecuacin (26) para la
caracterstica (para I=0) da una ecuacin en las 2 incgnitas Vp(0) y
Hp(0). Figura 1
Se necesita una condicin

exterior a la tubera para
resolver la ecuacin.
Fig. 1
Considerando que se sita un deposito en el extremo aguas arriba,
con altura HR por encima de la lnea central de la tubera.
Despreciando las perdidas menores en la entrada figura 2, se puede
aplicar la ecuacin de Euler no permanente. El ultimo termino es cero
excepto en la boca, y para esta distancia ds muy corta se puede
considerar pequea, de manera que se elimina, quedando.
Figura 2
S
Entrada a una
tubera desde un
deposito

Y para flujo de retorno.
Entonces entre la ecuacin (26 y nos dar HP (0) y VP (0).
Ahora para el extremo aguas abajo o de la derecha de la tubera

figura 3.
Figura 3 Vlvula
Se da la caracterstica mediante la ecuacin (25) para I=N por lo que
tambin se necesita una condicin externa (vlvula) Figura 4
Figura 4. vlvula en el extremo aguas debajo de una

tubera
Para un flujo permanente a travs de una vlvula puede ser

considerada como un orificio

A= rea de la tubera
= Velocidad en rgimen permanente en la tubera.
H0= Altura a travs de la vlvula
= El rea de orificio por el coeficiente de descarga.
En forma general (otro orificio)
(31) entre (30)

En (32) es el coeficiente sin dimensiones de abertura de la vlvula,
=1 para flujo permanente y se reduce a cero cuando la vlvula se
cierra.
Entre la ecuacin (32) y (25) se obtendr y
: Generalmente est en funcin de t y se puede expresar:
= Es el tiempo de cierre.
ALTERNATIVA DE LAS ECUACIONES CARACTERISTICAS

Flujo No Permanente - 07

Cargado por

Copyright:

Formatos disponibles

Flujo No Permanente - 07

Cargado por

Información del documento

Descripción original:

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

Flujo No Permanente - 07

Cargado por

Copyright:

Formatos disponibles

FLUJO NO

Autor: Ing. William Chauca Nolasco

Estas ecuaciones son las ecuaciones de cantidad

Para ello las variables dependientes son las

Haciendo un esquema para el anlisis de la cantidad de

D.C.L para la deduccin de la ecuacin del

Se considera que el rozamiento es proporcional

Ecuacion de movimiento: Elemento de fluido

Aplicando la ecuacin de movimiento para la grafica.

Desarrollando este termino

A esta ecuacin lo dividimos por una masa elemental de fluido

Introduciendo la lnea de alturas piezometricas

Luego desarrollando la ecuacin (2) teniendo en cuenta (3),

MULTIPLICANDO la ecuacin (5) POR -1

Por condiciones de flujo turbulento permanente.

Como el rozamiento se debe oponer al movimiento entonces

La ecuacin de continuidad para flujo no permanente

no es una funcin de t desarrollando (9)

Dividiendo por la masa

Si a la expresin (a) lo dividimos por t entonces.

b) Se determina la velocidad de variacin de la tensin.

Si a la (b) expresin lo dividimos por el modulo de elasticidad de Young

Si a (c) lo multiplicamos por el radio obtenemos:

d) Velocidad de alargamiento radial:

Si a (d) lo multiplicamos por el permetro de la tubera nos da:

e) La velocidad de aumento de rea.

Por definicin del modulo volumtrico de elasticidad del flujo (II)

La velocidad de variacin de la densidad dividida por la densidad:

Expresando la constante k con la relacin.

De esta expresin la variacin de la respecto a x o t es mucho menos

Si la tubera estuviera en reposo y .

Donde la ecuacin final ser:

Donde encontramos la segunda ecuacin caracterstica denominada .

La ecuacin es la ecuacin de la continuidad para un lquido

Lasecuaciones encontradas contienen 2 incgnitas . Se

Desarrollando la expresin y agrupando en forma conveniente.

Este agrupamiento se hace con el fin de

El primer trmino (a):

De (20) se ha hallado 2 valores reales de

Las 4 ecuaciones se pueden interpretar grficamente:

Se supone que se conoce V y H en R y S.

Entonces cada ecuacin (21) y (23) contiene 2 incgnitas para un

En la interseccin (P) de y ambas ecuaciones se pueden resolver

Para los clculos de golpe de ariete en tuberas metlicas, V es muy

Son rectas dichas caractersticas con pendiente

Se debe considerar que la tubera est construida en N tramos

Se debe conocer H y V inicialmente en cada una de las secciones en

Tambin deberamos tener en cuenta que la solucin se obtiene solo

utilizando la retcula dado que x y t se conocen en cada interseccin y

Escribiendo las ecuaciones (21) y (23) en notacin de diferencias

En forma anloga para (23)

Sumando (25) y (26) se elimina Vp(I) y obtenemos HP(I)

De la ecuacin (25) y (26) al restar obtendremos Vp (I). (25) - (26):

Se necesita una condicin

Y para flujo de retorno.

Entonces entre la ecuacin (26 y nos dar HP (0) y VP (0).

Ahora para el extremo aguas abajo o de la derecha de la tubera

Figura 4. vlvula en el extremo aguas debajo de una