Flujo de fluidos estratificados - PreMAT
Flujo de fluidos estratificados - PreMAT
Flujo de fluidos estratificados - PreMAT
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
2.2. TEOREMA DE SQUIRE 34<br />
u(y)<br />
Figura 2.3. Estado inicial para un fluido con velocidad y <strong>de</strong>nsidad<br />
variables.<br />
Supongamos que en este estado inicial vale la aproximación <strong>de</strong> Boussinesq presentada<br />
en la sección 1.10.1. La componente según ê2 <strong>de</strong> esta ecuación es simplemente<br />
0 = − ∂p0<br />
− ρ0(y)g.<br />
∂y<br />
Apliquemos una perturbación tridimensional sobre este estado inicial: consi<strong>de</strong>remos<br />
los campos perturbados<br />
ũ = (u + ũ, ˜v, ˜w),<br />
˜ρ = ρ0 + ρ,<br />
˜p = p0 + p.<br />
La ecuación <strong>de</strong> continuidad ∇·ũ = 0 se traduce en una ecuación <strong>de</strong> continuidad<br />
para las perturbaciones:<br />
(2.10)<br />
∂ũ ∂˜v ∂ ˜w<br />
+ + = 0.<br />
∂x ∂y ∂z<br />
Tomando como <strong>de</strong>nsidad <strong>de</strong> referencia ρ0, tenemos que el estado perturbado también<br />
<strong>de</strong>be satisfacer la ecuación <strong>de</strong> Boussinesq<br />
∂ũ<br />
∂t<br />
ê2<br />
ρ<br />
+ (ũ · ∇)ũ = g − 1<br />
∇p.<br />
Luego <strong>de</strong> linealizada esta ecuación, tenemos sus tres componentes:<br />
(2.11)<br />
∂ũ<br />
∂t<br />
∂˜v<br />
∂t<br />
∂ ˜w<br />
∂t<br />
∂ũ ∂u<br />
+ u ∂x + ˜v ∂y<br />
+ u ∂˜v<br />
∂x<br />
+ u ∂ ˜w<br />
∂x<br />
ρ0<br />
ρ0<br />
∂p<br />
1 = − ρ0 ∂x ,<br />
ρ 1 ∂p<br />
= − g − ρ0 ρ0 ∂y ,<br />
1 ∂p<br />
= − ρ0 ∂z .<br />
Como los coeficientes en las ecuaciones (2.10) y (2.11) no <strong>de</strong>pen<strong>de</strong>n <strong>de</strong> x, z, t,<br />
po<strong>de</strong>mos consi<strong>de</strong>rar modos normales <strong>de</strong> la forma<br />
ũ(x, y, z, t) = û(y)e i(kx+mz−ωt) ,<br />
˜v(x, y, z, t) = ˆv(y)e i(kx+mz−ωt) ,<br />
˜w(x, y, z, t) = ˆw(y)e i(kx+mz−ωt) ,<br />
ρ(x, y, z, t) = ˆρ(y)e i(kx+mz−ωt) ,<br />
p(x, y, z, t) = ˆp(y)e i(kx+mz−ωt) .<br />
Po<strong>de</strong>mos suponer k, m > 0 sin per<strong>de</strong>r generalidad, ya que el signo <strong>de</strong> Re(ω)<br />
alcanzaría para <strong>de</strong>terminar el sentido <strong>de</strong> <strong>de</strong>splazamiento <strong>de</strong> la perturbación. Una<br />
ê1