Propiedades Del Campo de Velocidad MF
Propiedades Del Campo de Velocidad MF
Propiedades Del Campo de Velocidad MF
Las propiedades cinemticas del campo de velocidad son determinadas por su divergencia, v, y por el rotor, v. Se adoptar un sistema de coordenadas con en el eje x positivo hacia el este, eje y positivo hacia el norte y el eje z hacia arriba en direccin de la lnea de accin de la fuerza gravitacional. Como el plano xy es tangente a la superficie de la tierra, las coordenadas xyz se conocen como coordenadas del plano tangente. Pero se aplica slo localmente, en un punto. COORDENADAS NATURALES. Es un sistema de coordenadas til para describir la dinmica del movimiento de un fluido. Aunque se pueden formular en tres dimensiones, veremos slo los aspectos en dos dimensiones. Usando las coordenadas naturales, se puede distinguir entre trayectoria y lnea de corriente de una parcela de fluido. Trayectoria: Es la curva descrita por las posiciones sucesivas de una parcela de fluido en movimiento. En un instante dado, el vector velocidad es tangente a la trayectoria. Lnea de corriente: Es una lnea cuya tangente en cualquier punto del fluido en movimiento es paralela a la velocidad instantnea. Las trayectorias y las lneas de corriente generalmente no coinciden, excepto en el caso de flujo estacionario. En la figura se ilustran ambas lneas
ww w.
Fi s
ic
aA
.c
om
VA FIGURA En la horizontal se usa una coordenada s a lo largo de la trayectoria o de la r lnea de corriente. La velocidad v que es tangente a la lnea en P define un ) vector unitario t tangente en P. Para mantener el sistema ortogonal derecho, ) se elige una coordenada normal n en P, con vector unitario n . La tercera ) coordenada puede ser la vertical local z con vector unitario k . VA FIGURA
Considerando una parcela de fluido que en un pequeo dt se mueve ds desde r P1 a P2, y donde v gira un ngulo d , que por convencin se considera positivo si la rotacin es antihorario. VA FIGURA
ww w.
Fi s
ic
aA
) ) ) ) ) ) ) ) ) t n = k, n k = t ; k t = n
.c
om
d ds es una derivada material, es decir que describe los cambios de la r direccin de v siguiendo la trayectoria de la parcela de fluido. Si la curva fuera una lnea de corriente, los puntos P1 a P2 seran para un mismo instante. En este caso se tiene: 1 = K s curvatura de la lnea de corriente, y Rs el radio de curvatura s Rs de la lnea de corriente. Segn la convencin > 0 si d > 0 K, Ks, R, Rs < 0 si d < 0
)
) ) ) ) )) ) ) ) dt ) dt dt ) dt ) d (t t ) t = 0 2t =0 t + t t =1 =0 t d d d d d
ww
w.
Fi s
El vector unitario t no es constante. Se puede demostrar que su derivada es ) ) ) normal a t . Por ejemplo si t = t ( ) , como es unitario
ic
)
aA
)
.c
om
Se deja como ejercicio demostrar que si t cambia a t +d t cuando gira un d , entonces ( pag, 103- 104 de Riegel)
) dt ) =n d ) t ) =n
Interesan tambin las derivadas de t respecto de n y s, es decir normal y a lo largo de la lnea de corriente o trayectoria. Por la regla de la cadena
) ) ) dt dt d ) dt = Kn trayectorias ) = = Kn ds d ds ds ) t ) = K s n lneas de corriente s De la misma forma: ) dt d ) = n trayectorias dn dn ) ) t n lneas de corriente = n n Tambin interesa la derivada de t respecto al tiempo:
) ) ) dt dt ds dt ) = =V = VK n dt ds dt ds )
) t ) = n t t
ww
w.
Fi s
) dt ) = Kvn dt
aA
ic
lnea de corriente
.c
trayectorias
om
DIVERGENCIA. Es una importante propiedad del campo de velocidad. La divergencia del campo de velocidad en el espacio es: u v w r r + + divv = v = x y z La velocidad instantnea, si vH es el viento horizontal, se puede escribir en coordenadas naturales como: r v H = v H t El operador nabla en coordenadas naturales es: H = t + n s n
r H vH = t s v r H vH = + s
(v H t ) (v H t ) n v H t = t + n n s n t t v t + n n s +
Reemplazando las definiciones de s y n se llega a: t t v H r H vH = + vH s n Indica que la divergencia horizontal del viento se debe a la contribucin de dos efectos: 1. Al trmino v/s llamado divergencia longitudinal, con v/s > 0 si el viento aumenta corriente abajo a lo largo de las lneas de corriente.
ww
w.
Fi s
ic
aA
.c
om
2. Al trmino v/n llamado divergencia transversal, con v/n > 0 si la lnea de corriente diverge en direccin normal al flujo. Estos efectos se ilustran en la figura:
Es posible que el flujo sea no divergente, aunque las lneas de corriente parezcan indicar convergencia o divergencia. El flujo es no divergente cuando el campo de viento se ajusta de manera que v/s se balancea con v/n. Geomtricamente, la divergencia se asocia con la tasa de cambio del rea por unidad de rea entre lneas de corriente, esto es:
ww w.
Si el movimiento es tal que el rea no cambia su valor numrico, aunque cambie su forma, el flujo es no divergente, como se ilustra en la figura:
Fi s
1dA r v = Adt
ic
aA
.c
om
VORTICIDAD. Matemticamente es el rotor del vector velocidad. En dinmica de fluidos, para el campo de velocidad, se llama vorticidad. Es una extensin del concepto de velocidad angular de una parcela de fluido que rota en torno a algn eje. r j Si v = ui + v + wk , la expresin matemtica de la vorticidad es definida por el vector q, v w w u v u r r q = v = j y z i + z x + x y k La vorticidad es un campo vectorial y se pueden dibujar curvas que son tangentes al vector q en cada punto. Se define lnea de vrtice a una curva tangente en cada punto al vector vorticidad en ese punto. De particular inters es la tendencia de la parcela de fluido geofsico a rotar en torno a la vertical local. Esta se representa por la componente vertical de la vorticidad, se denota con el smbolo v u r = k v = x y En meteorologa cuando se habla de vorticidad se refiere frecuentemente slo a la componente vertical . Como v o vH es el viento relativo, es decir el viento relativo a un punto sobre la superficie de la Tierra que gira, corresponde a la vorticidad relativa. A menos que se diga lo contrario se hablar de vorticidad refirindose a . La vorticidad tiene signo; en el hemisferio sur se tiene (NOTA: sentido ciclnico se refiere al mismo sentido de la rotacin terrestre):
ww w.
Fi s
ic
aA
.c
om
r v H = vt
r H v H = t + n (vt) n s Desarrollando esta expresin y usando las definiciones apropiadas, se llega a la relacin: v r H v H = vK S k n r y para la componente vertical de la vorticidad, calculando = k H v H ,
ww w.
= vK S
El primer trmino se llama de curvatura y el segundo se llama trmino de la cortante (del viento). Por lo tanto la vorticidad (o la rotacin de la parcela de fluido en torno a un eje a travs de su centro de masa) se debe a la
Fi s
ic
aA
v v v = n n RS
.c
superposicin de dos efectos: el efecto de la curvatura de la lnea de corriente y el de la cortante de la velocidad normal al flujo. Esto explica porque, como se ilustra en el esquema de ms abajo: 1. Un flujo en lnea recta puede tener vorticidad; este flujo no tiene curvatura, pero si existe una variacin de la velocidad normal a la direccin del flujo, la vorticidad no es cero. 2. Un flujo con curvatura puede ser irrotacional ( = 0) cuando el efecto de curvatura se balancea exactamente con el efecto de la cortante.
En un esquema de un flujo circular como un ejemplo de flujo curvado, si la parcela en direccin A tiene > 0 en el hemisferio sur, pero si la velocidad se ajusta adecuadamente, la parcela segn la direccin B tendr < 0, resultando = 0 en total. HS
ww
w.
Fi s
ic
aA
.c
om
En el hemisferio sur, la distribucin de velocidad y la curvatura alrededor de los sistemas de baja (alta) presin es tal que < 0 ( > 0), que llamamos vorticidad ciclnica (anticiclnica). En el HN cambian de signo. Esquema:
FUNCIN CORRIENTE. Si en un flujo horizontal se tiene div v = 0, se dice que es no divergente. Para estos flujos las componentes de la velocidad son tal que: u v u v r v = 0 = + = x y x y
Matemticamente, esta relacin requiere que la diferencial vdx udy sea exacta, es decir existe una funcin = (x,y) tal que:
ww
w.
Fi s
ic
aA
.c
om
d =
dx + dy = vdx udy x y
La funcin se llama funcin corriente, sus dimensiones son m2/s. La velocidad se puede escribir ahora como: r v = ui + v = j i + j y x r v = k La divergencia y la vorticidad en trminos de la funcin corriente, se expresan de la siguiente forma: 2 2 r + + 0 v = = xy yx x y y x r = k v = 2 Las isolneas de = cte son siempre tangente al vector velocidad, as que por definicin las isolneas = cte son lneas de corriente. Pero no se debe olvidar que las lneas de corriente existen siempre, an para flujos divergentes, en cambio las lneas = cte slo existen cuando el flujo es no divergente. VELOCIDAD POTENCIAL. Considerando ahora el caso de un flujo irrotacional, es decir con = 0, se tiene: v u r = 0 = k v = x y
ww w.
Fi s
ic
aA
.c
om
de donde se deduce que las componentes de la velocidad satisfacen la relacin v u = x y Este caso se da cuando el campo de viento es tal que existe un balance entre los efectos de curvatura y de la cortante del viento. Matemticamente se obtiene un flujo irrotacional cuando el campo de velocidad horizontal se puede obtener en trminos de una funcin escalar tal que v = , as se tiene: r v = = i + j x y ,v = entonces : u = x y
De esta forma se cumple que = 0, ya que: 2 2 r 0 = k v = k = xy yx Para la divergencia se obtiene: r v = = 2 Las isolneas de = cte se llaman equipotenciales y son perpendiculares al flujo (que es irrotacional). Se dibujan en las figuras.
ww w.
Fi s
ic
aA
.c
om
En el caso especial en que flujo sea no divergente e irrotacional, se pueden representar en trminos slo de de . En este caso: r v = 2 = 0
= 2 = 0
es decir las funciones corriente y potencial satisfacen la ecuacin de Laplace, por lo que se llama un flujo laplaciano. En este caso y son funciones armnicas, ortogonales que satisfacen las ecuaciones de Cauchy-Riemann: u v u v = ; = x y y x La mayora de los flujos reales se pueden aproximar con distintos grados de realismo a los tres casos anteriores. Pero en general los flujos (especialmente el atmosfrico) no son ni puramente divergente ni puramente irrotacional, y adems interesa tener en cuentas las divergencias y rotaciones del flujo. Para hacer esto podemos elegir una funcin F vectorial que se escriba en la forma:
r r F = h +
ww w.
donde h se llama potencial vectorial y potencial escalar de F. Escrita F de esta forma, el primer trmino contiene todo el rotor de F y el segundo toda su divergencia. Si se elige h = z , entonces:
r h = k
Fi s
ic
aA
.c
om
que se conoce como el teorema de Helmholtz, que dice que el campo de viento se puede descomponer en sus partes rotacional y divergente. De esta forma se puede considerar el campo de velocidad puramente irrotacional, puramente no divergente, o ambos o ninguno. Tambin se llaman a y funcin corriente y velocidad potencial, pero su forma funcional difieren de y .
El trmino k se llama parte rotacional del viento y se denota por v, y el trmino se llama parte divergente del viento y se denota por v, as se escribe:
r r r v = v + v
LINEAS DE CORRIENTE. En todo punto v es tangente a la lnea de corriente, y en dos dimensiones dR r dR tambin es tangente. Adems t = , y como v es paralelo a dR, entonces ds r r v dR = 0 , entonces: j j (ui + v) (dxi + dy) = 0 (udy vdx)k = 0 dy v( x, y, t ) = dx v( x, y, t )
ww w.
es la ecuacin diferencial de las lneas de corriente. Si en dos dimensiones se escribe en la forma vdx udy = 0 y si el flujo es no divergente, entonces existe tal que u = -/y, v = /x, reemplazando se tiene
Fi s
ic
aA
.c
om
dx + dy = d = 0 x y = cte. lo que indica que la funcin coincide con las lneas de corriente. TRAYECTORIAS. Es la curva que se forma al unir las posiciones sucesivas de una parcela de fluido en movimiento, que se puede definir matemticamente por la ecuacin r dr r r = v (r , t ) dt con r y v en tres dimensiones. Cuando se determina r(t) para alguna condicin inicial t = to, entonces se conoce la trayectoria de la parcela de fluido. En la mayora de los casos, la solucin de la ecuacin diferencial es dificil, y debe simplificarse con suposiciones o usar mtodos de solucin numrica. Ejemplo. Para un flujo con componente zonal constante y meridional con perfil sinusoidal de la forma: u = U = cte > 0 , v = vo cos k ( x ct ) La ecuacin de las lneas de corriente es: dy v vo = = cos k ( x ct ) dx u U Integrando se obtiene: ys = vo sen k ( x ct ) + yo Uk vo, k y c ctes.
ww w.
Fi s
ic
aA
.c
om
con yo cte de integracin. Las lneas de corriente son sinusoidales, de amplitud vo/Uk.
que corresponde a las lneas de corriente, y por comparacin se tiene que = Uyo tambin es lnea de corriente. El clculo de la trayectoria da: dx = u =U dt dy = v = vo cos k ( x ct ) dt x(t ) = Ut + xo
ww
vo sen k ( x ct ) Uk U
w.
Fi s
ic
aA
vo sen k ( x ct )dx Uy k
.c om
Si x = 0 para t = 0, entonces xo = 0, y x(t ) = Ut Se reemplaza en dy/dt: dy = vo cos k (Ut ct ) = vo cos k (U c)t dt vo y (t ) = sen k (U c)t + yo k (U c) Para y = 0 en t = 0, entonces yo = 0, y las ecuaciones paramtricas de la trayectoria si la parcela de fluido est en el origen cuando t = 0 son: x(t ) = Ut y (t ) = vo sen k (U c)t k (U c)
om
y ( x) =
vo c sen k 1 x k (U c) U
de donde se observa que la trayectoria no coincide con las lneas de corriente. En la figura se dibujan para el caso en que la lnea de corriente se mueve con velocidad de fase c = U/2.
ww
w.
Fi s
ic
aA
.c