25.10.2014 Views

Golpe de Ariete

Golpe de Ariete

Golpe de Ariete

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

<strong>Golpe</strong> <strong>de</strong> <strong>Ariete</strong><br />

Resta<br />

Mejorar el mo<strong>de</strong>lo (complejizar(<br />

el cálculo)<br />

incorporando los efectos <strong>de</strong> fricción,<br />

topográficos y geometrías complejas<br />

Presentar como se analizan los dispositivos<br />

<strong>de</strong> alivio<br />

Die<strong>de</strong>rik Korteweg<br />

1848-1941<br />

Efecto <strong>de</strong> la fricción sobre el golpe<br />

<strong>de</strong> ariete<br />

Los coeficientes <strong>de</strong> fricción son<br />

consi<strong>de</strong>rados en general iguales tanto<br />

para flujos estacionarios como para no<br />

estacionarios<br />

Algunas <strong>de</strong> las discrepancias ocurren<br />

como consecuencia <strong>de</strong> usar velocida<strong>de</strong>s<br />

promedios en lugar <strong>de</strong> las verda<strong>de</strong>ras<br />

distribuciones <strong>de</strong> velocidad<br />

Pérdidas Localizadas<br />

Se supone que los coeficientes son los<br />

mismos que para el caso estacionario<br />

Pérdidas pequeñas localizadas conviene<br />

distribuir la pérdida aumentando el factor<br />

<strong>de</strong> fricción global <strong>de</strong> la cañería<br />

Perdidas localizadas importantes (válvulas<br />

parcialmente cerradas) conviene<br />

consi<strong>de</strong>rara a la válvula como la frontera<br />

entre dos cañerías<br />

Descripción simplificada <strong>de</strong>l fenómeno <strong>de</strong> golpe <strong>de</strong> ariete<br />

con fricción en la cañería<br />

Atenuación <strong>de</strong> la onda <strong>de</strong> choque<br />

∆p p c.<br />

inst<br />

a<br />

∆p<br />

fricción<br />

Variaciones <strong>de</strong><br />

presión<br />

producidas por el<br />

golpe <strong>de</strong> ariete<br />

en la válvula <strong>de</strong><br />

seccionamiento<br />

T=2L/a<br />

a) Mo<strong>de</strong>lo Sin<br />

disipación<br />

b) Mo<strong>de</strong>lo Con<br />

Disipación causada<br />

por fricción<br />

U<br />

U=0<br />

c) Mo<strong>de</strong>lo con<br />

disipación incluye<br />

fricción y transferencia<br />

<strong>de</strong> calor<br />

L<br />

T


Resonancia y <strong>Golpe</strong> <strong>de</strong> <strong>Ariete</strong><br />

Números adimensionales <strong>de</strong> interés<br />

El número <strong>de</strong> disipación Dn es la relación entre los tiempos <strong>de</strong> viaje <strong>de</strong> la onda y un término<br />

que representa el <strong>de</strong>caimiento viscoso<br />

n−1<br />

⎛ υ ⎞⎛<br />

L ⎞<br />

f n V0<br />

L<br />

Dn =<br />

⎜<br />

⎟<br />

2<br />

⎜ ⎟ Laminar<br />

Dn =<br />

⎝ r ⎠⎝<br />

a ⎠<br />

32 r a<br />

Turbulento<br />

f factor <strong>de</strong> fricción <strong>de</strong>l diagrama <strong>de</strong> Moody , n exponente <strong>de</strong> la velocidad en<br />

el término <strong>de</strong> pérdidas por fricción (n≈1.8<br />

1.8-2.0)<br />

Numero <strong>de</strong> propagación Pn: : relación entre el tiempo <strong>de</strong> tránsito en la línea con respecto a la<br />

frecuencia natural <strong>de</strong> la línea ω/2 (<strong>de</strong>pen<strong>de</strong><br />

<strong>de</strong> soportes, espesor,…)<br />

2 ω L<br />

Pn =<br />

π a<br />

Para análisis <strong>de</strong> la atenuación Dn


Método <strong>de</strong> las características<br />

Para que el sistema<br />

admita solución :<br />

⎡ dx<br />

⎢<br />

0<br />

<strong>de</strong>t ⎢<br />

⎢ 1<br />

⎢<br />

⎢⎣<br />

0<br />

2<br />

( ρ a )<br />

dt 0 0 ⎤<br />

0 dx dt<br />

⎥<br />

⎥<br />

0 0 1 ⎥<br />

⎥<br />

1 1/ ρ 0 ⎥⎦<br />

= 0<br />

Curvas características<br />

Las curvas características son líneas <strong>de</strong>l plano x-t x<br />

para las cuales el coeficiente angular dx/dt<br />

dt es igual a<br />

la velocidad <strong>de</strong> propagación <strong>de</strong> pequeñas<br />

perturbaciones.<br />

Lo que conduce a la<br />

obtención <strong>de</strong> la<br />

expresión <strong>de</strong> las<br />

curvas características<br />

dx<br />

= a<br />

dt<br />

dx<br />

= −a<br />

dt<br />

← Γ<br />

← Γ<br />

+<br />

−<br />

Sobre estas curvas se verifican lo que se conoce<br />

como relaciones <strong>de</strong> compatibilidad.<br />

Las relaciones <strong>de</strong> compatibilidad se obtiene a partir<br />

<strong>de</strong>l <strong>de</strong>terminante característico<br />

Determinantes característicos<br />

⎡ dx<br />

⎢<br />

⎢<br />

0<br />

⎢ 1<br />

⎢<br />

⎢⎣<br />

0<br />

2<br />

( ρ a )<br />

dt 0 0 ⎤<br />

0 dx dt<br />

⎥<br />

⎥<br />

0 0 1 ⎥<br />

⎥<br />

1 1/ ρ 0 ⎥⎦<br />

⎡∂U<br />

dU<br />

U x ⎤<br />

⎢ ⎥<br />

⎢∂<br />

∂<br />

dp<br />

p<br />

t ⎥<br />

⎢ ⎥ =<br />

⎢<br />

∂<br />

∂<br />

U<br />

p x ⎥<br />

⎢ ⎥ − f − g senθ<br />

⎢<br />

∂ ∂<br />

0<br />

2<br />

⎥<br />

⎣ ∂t<br />

⎦<br />

2D<br />

t<br />

Cálculo por el método <strong>de</strong> las<br />

características<br />

Γ − Γ+<br />

D<br />

⎡ dU dt 0 0 ⎤<br />

⎢<br />

0<br />

⎥<br />

⎢<br />

dp<br />

dx dt<br />

⎥<br />

<strong>de</strong>t⎢<br />

0 0 0 1 ⎥ = 0<br />

2<br />

⎢ U<br />

⎥<br />

⎢−<br />

f − g senθ 1 1/ ρ 0 ⎥<br />

⎣ 2D<br />

⎦<br />

⎡ dx<br />

⎢<br />

⎢<br />

0<br />

<strong>de</strong>t⎢<br />

1<br />

⎢<br />

⎢<br />

⎢ 0<br />

⎣<br />

dt dU 0 ⎤<br />

0 dp dt<br />

⎥<br />

⎥<br />

⎥<br />

1<br />

2<br />

U<br />

⎥<br />

− f − g sen θ 0 ⎥<br />

2D<br />

⎦<br />

2 0 0 1 ⎥ = 0<br />

( ρ a )<br />

2<br />

1 ⎛ U ⎞<br />

dU + dp +<br />

⎜ f + g senθ<br />

⎟ a dx = 0<br />

ρa<br />

⎝ 2g<br />

⎠<br />

2<br />

1 ⎛ U ⎞<br />

dU − dp −<br />

⎜ f + g senθ<br />

⎟ a dx = 0<br />

ρa<br />

⎝ 2g<br />

⎠<br />

+<br />

← Γ<br />

−<br />

← Γ<br />

x<br />

X − X<br />

D<br />

X − X<br />

D<br />

U −U<br />

+<br />

D<br />

U −U<br />

−<br />

D<br />

A<br />

B<br />

A<br />

B<br />

= a<br />

= a<br />

1<br />

( TD<br />

−TA<br />

)<br />

( T −T<br />

)<br />

( PD<br />

− PA<br />

) + U<br />

A<br />

U<br />

D( X<br />

D<br />

− X<br />

A<br />

)<br />

ρ a 2 d<br />

1<br />

B<br />

A<br />

D<br />

B<br />

f<br />

f<br />

+ a g senθ<br />

= 0<br />

( PD<br />

− PB<br />

) + U<br />

B<br />

U<br />

D<br />

( X<br />

D<br />

− X<br />

B<br />

) + a g senθ<br />

= 0<br />

ρ a 2 d<br />

a1<br />

U=0<br />

TIEMPO<br />

L1/a1+<br />

L2/a2<br />

L1/a1<br />

Condición Frontera<br />

P=cte<br />

Condición p1=p2<br />

+ cons. masa<br />

∆U0<br />

Condición Frontera U=0<br />

DISTANCIA


P=cte<br />

U=0<br />

a1<br />

Dispositivos <strong>de</strong> Protección frente<br />

al golpe <strong>de</strong> ariete<br />

TIEMPO<br />

Igualdad <strong>de</strong><br />

Presiones a<br />

ambos lados<br />

Conservación<br />

<strong>de</strong> la masa<br />

U=0<br />

D<br />

∆t<br />

=<br />

∆x<br />

a<br />

A<br />

X D − X A = a ( TD<br />

−TA<br />

)<br />

X D − X B = a ( TD<br />

−TB<br />

)<br />

1<br />

f<br />

U D −U<br />

A + ( PD<br />

− PA<br />

) + U A U D<br />

( X D − X A<br />

) + a g senθ<br />

= 0<br />

ρ a 2 d<br />

1<br />

f<br />

U D −U<br />

B − ( PD<br />

− PB<br />

) + U B U D<br />

( X D − X B<br />

) + a g senθ<br />

= 0<br />

ρ a 2 d<br />

B<br />

C<br />

DISTANCIA<br />

Estudio Técnico <strong>de</strong> prevención <strong>de</strong><br />

golpe <strong>de</strong> ariete<br />

Especificación <strong>de</strong>l régimen transitorio estudiado:<br />

• régimen inicial y origen <strong>de</strong>l transitorios<br />

Datos necesarios para llevar a<strong>de</strong>lante el estudio:<br />

• Reservorios, bombas, Purgas, Válvulas unidireccionales, Válvulas<br />

<strong>de</strong> cierre, Puntos muertos, Dispositivos <strong>de</strong> regulación<br />

Cálculo<br />

Evaluación <strong>de</strong>l riesgo <strong>de</strong> golpe <strong>de</strong> ariete<br />

Presentación <strong>de</strong> resultados<br />

Resumiendo<br />

Hemos visto una <strong>de</strong>scripción cualitativa <strong>de</strong>l golpe <strong>de</strong> ariete y dimos la<br />

explicación <strong>de</strong>l orígen <strong>de</strong> las sobrepresiones y subpresiones.<br />

Analizamos las mismas en el marco <strong>de</strong> las leyes <strong>de</strong> conservación <strong>de</strong> la<br />

mecánica <strong>de</strong> fluidos que conocemos.<br />

Logramos expresiones que permiten una estimación rápida <strong>de</strong> la<br />

importancia <strong>de</strong> las sobrepresiones y subpresiones que se pue<strong>de</strong>n<br />

<strong>de</strong>sarrollar en este fenómeno que consi<strong>de</strong>ran la maniobra <strong>de</strong> cierre.<br />

e.<br />

Efectuamos un análisis <strong>de</strong> casos simples en el plano x-t x t que sirve <strong>de</strong><br />

base para consi<strong>de</strong>rar los efectos <strong>de</strong> fricción y las singularida<strong>de</strong>s s <strong>de</strong>l<br />

conducto por el método <strong>de</strong> las características<br />

Presentamos algunos <strong>de</strong> los elementos <strong>de</strong> alivio <strong>de</strong>l golpe <strong>de</strong> ariete

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!