07.05.2013 Views

Flujo de fluidos estratificados - PreMAT

Flujo de fluidos estratificados - PreMAT

Flujo de fluidos estratificados - PreMAT

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

3.1. CONTROLES HIDRÁULICOS EN FLUJOS DISCRETAMENTE ESTRATIFICADOS 48<br />

alguna perturbación. Esta <strong>de</strong>finición generaliza a la presentada en el capítulo 1.<br />

Recor<strong>de</strong>mos que en el caso <strong>de</strong> un fluido homogéneo, el estado crítico representa un<br />

mínimo en la energía interna y en el momento en función <strong>de</strong>l tirante <strong>de</strong>l canal. En<br />

los <strong>fluidos</strong> <strong>estratificados</strong> ocurre algo similar.<br />

3.1.1. <strong>Flujo</strong> crítico para estratificaciones discretas con superficie libre.<br />

Supongamos que tenemos un sistema formado por n capas homogéneas, numeradas<br />

<strong>de</strong> arriba hacia abajo, <strong>de</strong> <strong>de</strong>nsida<strong>de</strong>s ρi, espesores yi, <strong>de</strong>splazándose con<br />

velocida<strong>de</strong>s ui, y con una superficie libre (figura 3.1). Asumiremos también que el<br />

ρ1 u1(x)<br />

ρ2<br />

u2(x)<br />

y1(x)<br />

y2(x)<br />

ρn un(x) yn(x)<br />

p = p0<br />

Figura 3.1. <strong>Flujo</strong> <strong>de</strong> n capas uniformes a superficie libre.<br />

flujo <strong>de</strong> cada capa es unidireccional, y que estamos en hipótesis <strong>de</strong> aguas someras.<br />

Consi<strong>de</strong>remos conocidos y fijos los caudales por unidad <strong>de</strong> ancho, qi.<br />

Sea Ei la energía interna <strong>de</strong> la i-ésima capa, entonces tendremos<br />

i−1<br />

n<br />

(3.1) Ei = q2 i<br />

2gy 2 i<br />

+<br />

j=1<br />

rjiyj +<br />

don<strong>de</strong> rji = ρj/ρi. Impongamos una perturbación infinitesimal, y tomemos un<br />

referencial solidario con dicha perturbación. Entonces, la altura <strong>de</strong> alguna <strong>de</strong> las<br />

capas <strong>de</strong>be variar, esto es, ∂yj/∂x = 0 para algún j. Asimismo, si no tenemos en<br />

cuenta los efectos <strong>de</strong> la fricción, la energía interna <strong>de</strong> cada capa <strong>de</strong>be permanecer<br />

constante:<br />

(3.2)<br />

<br />

∂E<br />

∂y<br />

ij<br />

= ∂Ei<br />

∂yj<br />

∂Ei<br />

∂x =<br />

n ∂Ei<br />

j=1<br />

y ∂y<br />

∂x =<br />

∂yj<br />

∂y1<br />

∂x<br />

∂yj<br />

∂x<br />

, . . . , ∂yn<br />

∂x<br />

j=i<br />

yj,<br />

= 0 ∀i = 1, . . . , n.<br />

Tenemos así un sistema <strong>de</strong> n ecuaciones, una para la energía interna <strong>de</strong> cada capa,<br />

con n incógnitas,<br />

∂E ∂y<br />

= 0,<br />

∂y ∂x<br />

⊤ don<strong>de</strong><br />

. Como estamos suponiendo que<br />

por lo menos una <strong>de</strong> las capas <strong>de</strong>be verificar ∂yj/∂x = 0, este sistema <strong>de</strong>be tener<br />

soluciones no triviales. Tenemos así una condición para estado crítico:<br />

<br />

<br />

(3.3)<br />

<br />

∂E <br />

<br />

∂y = 0.<br />

x

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

Saved successfully!

Ooh no, something went wrong!