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.

2.5. INTERFASES DE ESPESOR NO NULO ENTRE CAPAS INFINITAS 40<br />

Figura 2.5. Perfiles <strong>de</strong> velocidad y <strong>de</strong>nsidad iniciales (2.23).<br />

De esta forma, se tiene que V es una medida típica <strong>de</strong> velocidad, y σ una <strong>de</strong> la<br />

variación relativa <strong>de</strong> <strong>de</strong>nsidad. La frecuencia <strong>de</strong> flotabilidad es<br />

N 2 (y) = σg<br />

d β′ (y),<br />

y el número <strong>de</strong> Richardson local es<br />

Ri(y) = Jβ′ (y)<br />

(u ′ σgd<br />

, don<strong>de</strong> J =<br />

(y)) 2 V<br />

La elección <strong>de</strong> estas escalas para U y ρ implican que J es el número <strong>de</strong> Richardson<br />

en y = 0. Para el estado inicial (2.23), si escribimos ρ1 = rρ2, tenemos<br />

V = u1 − u2<br />

2<br />

y σ =<br />

1 − r<br />

1 + r ,<br />

<strong>de</strong> modo que<br />

4(1 − r)gd<br />

J =<br />

.<br />

(1 + r)(u1 − u2) 2<br />

Con estas nuevas variables, la ecuación <strong>de</strong> Taylor-Goldstein en forma adimensionalizada<br />

es<br />

∂2 ˆ <br />

′′ ψ u Jβ′<br />

= − + α2 ˆψ,<br />

∂y2 u − c (u − c) 2<br />

don<strong>de</strong> c es la velocidad <strong>de</strong> fase adimensionalizada. Observemos que si ˆ ψ es solución<br />

<strong>de</strong> la ecuación anterior con un valor propio asociado c, entonces su conjugada compleja<br />

ˆ ψ ∗ es solución para el valor propio c ∗ . Luego, para cada modo que <strong>de</strong>cae en<br />

el tiempo existe uno que crece. Por lo tanto, una solución ˆ ψ es estable si y sólo si<br />

su valor propio asociado es real.<br />

Siguiendo [21] (sección 44) y [42], la curva <strong>de</strong> estabilidad marginal en el plano<br />

α − J está dada por J = α(1 − α), y obtenemos el diagrama presentado en la figura<br />

2.6.<br />

Queremos aproximar el estado inicial dado por (2.23) mediante dos capas homogéneas,<br />

<strong>de</strong> velocida<strong>de</strong>s u1, u2 y <strong>de</strong>nsida<strong>de</strong>s ρ1, ρ2, respectivamente. Para este<br />

estado inicial, que fue consi<strong>de</strong>rado en la sección anterior, obtuvimos que los modos<br />

neutralmente estables tienen número <strong>de</strong> onda<br />

k = g(1 − r2 )<br />

.<br />

r(u1 − u2) 2<br />

2 .

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

Saved successfully!

Ooh no, something went wrong!