Physics · Libro 4 · Bachelor Year 2

Física universitaria — segundo año

Física universitaria — segundo año · Bachelor Year 2

2Cinemática de fluidos

Asómese uno al pretil de un puente y mire el río. Una hoja pasa a la deriva, se acelera entre las pilas, se remansa en la poza de aguas abajo, da una vuelta en un remolino y sigue. Podría describirse el río siguiendo esa hoja — su posición y su velocidad en cada instante — o quedándose quieto y registrando, en cada punto del río, con qué rapidez pasa el agua. La segunda descripción es la del físico de fluidos: un campo de velocidades. Este capítulo establece las dos descripciones, aprende a calcular la aceleración de una partícula a partir del campo, escribe la conservación de la masa como una ley local y clasifica los flujos según giren o no — el vocabulario que necesitarán las ecuaciones de Euler y de Navier.

Un río a su paso por un puente: agua rápida y lisa entre las pilas, remolinos a su abrigo — un campo de velocidades que se lee en puntos fijos, no una única trayectoria que se sigue.
Un río a su paso por un puente: agua rápida y lisa entre las pilas, remolinos a su abrigo — un campo de velocidades que se lee en puntos fijos, no una única trayectoria que se sigue.

2.1 La descripción como medio continuo

Definición 2.1 (Partícula de fluido; la hipótesis del medio continuo)

Una partícula de fluido es un volumen de fluido lo bastante pequeño para tratarse como un punto a la escala del flujo, y a la vez lo bastante grande para contener un número enorme de moléculas, de modo que su densidad ρ\rho, su presión PP, su temperatura TT y su velocidad media v\vect v sean promedios bien definidos (la escala mesoscópica, típicamente el micrómetro). La hipótesis del medio continuo supone que tal escala existe: el recorrido libre medio \ell de las moléculas debe ser mucho menor que la escala LL del flujo, /L1\ell/L \ll 1. En el aire a presión atmosférica 70nm\ell \approx 70\,\mathrm{nm}; en un líquido \ell es el tamaño molecular. El fluido queda entonces descrito por campos ρ(M,t)\rho(M, t), P(M,t)P(M, t), v(M,t)\vect v(M, t) definidos en todo punto.

Definición 2.2 (Descripciones lagrangiana y euleriana)

La descripción lagrangiana sigue cada partícula de fluido: su posición r(t)\vect r(t) y su velocidad  ⁣dr/ ⁣dt\dd\vect r/\dd t como funciones del tiempo y de su posición inicial. La descripción euleriana da, en cada punto fijo MM y en cada instante tt, la velocidad v(M,t)\vect v(M, t) de la partícula que pasa por MM en ese momento. Una trayectoria es el camino recorrido por una partícula; una línea de corriente en el instante tt es una curva tangente en cada uno de sus puntos a v(M,t)\vect v(M, t); un tubo de corriente es la superficie formada por las líneas de corriente que pasan por una curva cerrada. Un flujo es estacionario cuando los campos eulerianos no dependen del tiempo, v/t=0\partial\vect v/\partial t = \vect 0; entonces las líneas de corriente y las trayectorias coinciden y son fijas.

Ejemplo 2.3 (Dos descripciones de un mismo flujo)

El agua de una manguera de jardín de sección SS con caudal QQ se mueve a v=Q/Sv = Q/S en todas partes: el campo euleriano es uniforme y el movimiento lagrangiano de cada partícula es uniforme. En la boquilla que se estrecha el campo euleriano sigue siendo estacionario, pero vv crece a lo largo del eje: cada partícula se acelera al pasar, aunque nada cambie con el tiempo en ningún punto fijo. Eso es lo primero que la descripción euleriana debe aprender a expresar.

A la izquierda, un flujo convergente estacionario: las líneas de corriente son fijas, una partícula marcada sigue una de ellas y se acelera a medida que el tubo se estrecha. A la derecha, la descripción euleriana registra el vector velocidad en puntos fijos; la lagrangiana sigue una partícula a lo largo de su trayectoria.
A la izquierda, un flujo convergente estacionario: las líneas de corriente son fijas, una partícula marcada sigue una de ellas y se acelera a medida que el tubo se estrecha. A la derecha, la descripción euleriana registra el vector velocidad en puntos fijos; la lagrangiana sigue una partícula a lo largo de su trayectoria.

2.2 La derivada material

Teorema 2.4 (Derivada material (o particular))

La tasa de variación de una magnitud G(M,t)G(M, t) (escalar o vectorial) para la partícula de fluido que pasa por MM en el instante tt es

DGDt=Gt+(vgrad)G,en particulara=DvDt=vt+(vgrad)v,\frac{\mathrm DG}{\mathrm Dt} = \frac{\partial G}{\partial t} + (\vect v\cdot \operatorname{\vect{grad}})\,G , \qquad\text{en particular}\qquad \vect a = \frac{\mathrm D\vect v}{\mathrm Dt} = \frac{\partial\vect v}{\partial t} + (\vect v\cdot\operatorname{\vect{grad}})\,\vect v ,

donde (vgrad)=vxx+vyy+vzz(\vect v\cdot\operatorname{\vect{grad}}) = v_x\partial_x + v_y\partial_y + v_z\partial_z. El primer término es la tasa de variación local (en un punto fijo) y el segundo la convectiva (porque la partícula se desplaza a donde GG es distinto). Para la velocidad, (vgrad)v=grad(v2/2)+(curlv)v(\vect v\cdot \operatorname{\vect{grad}})\vect v = \operatorname{\vect{grad}}(v^2/2) + (\operatorname{\vect{curl}}\vect v)\wedge\vect v.

Demostración. En  ⁣dt\dd t la partícula pasa de MM a M+v ⁣dtM + \vect v\,\dd t, luego  ⁣dG=G(M+v ⁣dt,t+ ⁣dt)G(M,t)=tG ⁣dt+ ⁣dt(vxx+vyy+vzz)G\dd G = G(M + \vect v\dd t, t + \dd t) - G(M, t) = \partial_tG\,\dd t + \dd t\,(v_x\partial_x + v_y\partial_y + v_z\partial_z)G a primer orden (regla de la cadena para una función de cuatro variables). La identidad vectorial se comprueba sobre la componente xx: [grad(v2/2)]x=vxxvx+vyxvy+vzxvz[\operatorname{\vect{grad}}(v^2/2)]_x = v_x\partial_xv_x + v_y\partial_x v_y + v_z\partial_xv_z y [(curlv)v]x=vy(yvxxvy)+vz(zvxxvz)[(\operatorname{\vect{curl}}\vect v)\wedge\vect v]_x = v_y(\partial_yv_x - \partial_xv_y) + v_z(\partial_zv_x - \partial_xv_z) (con las componentes de curlv\operatorname{\vect{curl}}\vect v dadas más abajo en Definición 2.12); su suma vale vxxvx+vyyvx+vzzvxv_x\partial_xv_x + v_y\partial_yv_x + v_z\partial_zv_x.

Ejemplo 2.5 (Aceleración en una boquilla estacionaria)

Un flujo estacionario unidimensional v(x)=v0(1+x/L)v(x) = v_0(1 + x/L) a lo largo de una boquilla: el término local se anula y a=v ⁣dv/ ⁣dx=v02(1+x/L)/La = v\,\dd v/\dd x = v_0^2(1 + x/L)/L. Una partícula que entra a v0=2m/sv_0 = 2\,\mathrm{m}/\mathrm{s} y sale a 2v02v_0 tras un L=5cmL = 5\,\mathrm{cm} soporta a=80a = 80 a 160m/s2160\,\mathrm{m}/\mathrm{s}^{2}: de ocho a dieciséis gg, en un flujo donde nada depende del tiempo.

Ejemplo 2.6 (Local frente a convectivo)

La temperatura en una estación meteorológica baja durante la noche (T/t<0\partial T/ \partial t < 0: local), y baja también para una masa de aire arrastrada hacia el norte por el viento (vgradT<0\vect v\cdot\operatorname{\vect{grad}}T < 0: convectiva, o advectiva). Un termómetro montado en un globo que va a la deriva con el viento registra la suma, DT/Dt\mathrm DT/\mathrm Dt; la estación registra solo el primer término.

2.3 Conservación de la masa

Definición 2.7 (Caudal másico y caudal volumétrico)

El caudal másico a través de una superficie orientada SS es la masa que la atraviesa por unidad de tiempo, Dm=ρvn ⁣dSD_m = \sum\rho\,\vect v\cdot\vect n\,\dd S — el flujo de la densidad de corriente de masa j=ρv\vect j = \rho\vect v (unidad kgm2s1\mathrm{kg}\,\mathrm{m}^{-2}\,\mathrm{s}^{-1}); el caudal volumétrico es DV=vn ⁣dSD_V = \sum\vect v \cdot\vect n\,\dd S (m3/s\mathrm{m}^{3}/\mathrm{s}). Para una velocidad uniforme y normal a una sección plana de área SS, Dm=ρvSD_m = \rho vS y DV=vSD_V = vS.

Demostración. En  ⁣dt\dd t, el fluido que atraviesa el elemento  ⁣dS\dd S es el contenido en el cilindro oblicuo de base  ⁣dS\dd S y generatriz v ⁣dt\vect v\,\dd t, de volumen vn ⁣dS ⁣dt\vect v\cdot\vect n\,\dd S\,\dd t y masa ρ\rho veces ese volumen.

Teorema 2.8 (Conservación local de la masa)

En todo punto de un fluido,

ρt+div(ρv)=0,dondedivA=Axx+Ayy+Azz\frac{\partial\rho}{\partial t} + \operatorname{div}(\rho\vect v) = 0 , \qquad\text{donde}\qquad \operatorname{div}\vect A = \frac{\partial A_x}{\partial x} + \frac{\partial A_y}{\partial y} + \frac{\partial A_z}{\partial z}

es la divergencia de un campo vectorial: el flujo neto saliente de A\vect A a través de la superficie de un volumen pequeño, por unidad de volumen. De forma equivalente, Dρ/Dt+ρdivv=0\mathrm D\rho/\mathrm Dt + \rho\operatorname{div}\vect v = 0.

Demostración. Tómese la caja fija [x,x+ ⁣dx]×[y,y+ ⁣dy]×[z,z+ ⁣dz][x, x + \dd x]\times[y, y + \dd y]\times[z, z + \dd z]. Su masa ρ ⁣dx ⁣dy ⁣dz\rho\,\dd x\dd y\dd z varía a la tasa tρ ⁣dx ⁣dy ⁣dz\partial_t\rho\,\dd x\dd y \dd z. Entra masa por la cara situada en xx a la tasa jx(x) ⁣dy ⁣dzj_x(x)\dd y\dd z y sale por la cara situada en x+ ⁣dxx + \dd x a la tasa jx(x+ ⁣dx) ⁣dy ⁣dzj_x(x + \dd x)\dd y\dd z: salida neta xjx ⁣dx ⁣dy ⁣dz\partial_xj_x\,\dd x\dd y\dd z, y análogamente para yy y zz. La conservación de la masa — no se crea masa — da tρ=(xjx+yjy+zjz)\partial_t\rho = - (\partial_xj_x + \partial_yj_y + \partial_zj_z). La segunda forma se sigue de div(ρv)=vgradρ+ρdivv\operatorname{div}(\rho\vect v) = \vect v\cdot\operatorname{\vect{grad}}\rho + \rho\operatorname{div}\vect v.

Observación 2.9 (La divergencia, un primer encuentro)

La deducción muestra qué mide div\operatorname{div}: el flujo que sale de una caja pequeña por unidad de volumen. Sumando sobre las cajas que teselan un volumen finito VV, las caras interiores se cancelan dos a dos y solo queda la superficie exterior: VAn ⁣dS=VdivA ⁣dτ\iint_{\partial V}\vect A\cdot\vect n\,\dd S = \iiint_V \operatorname{div}\vect A\,\dd\tau — el teorema de la divergencia (de Ostrogradski), que se enunciará como es debido en Capítulo 11 y se usará en todo el electromagnetismo; se demuestra en el volumen de matemáticas del Año 3. Integrada sobre un volumen fijo, la ley local se lee  ⁣dmV/ ⁣dt=Dmsal\dd m_V/\dd t = -D_m^{\text{sal}}: la masa interior varía en lo que atraviesa la frontera.

A la izquierda, el balance de masa de una caja fija: el flujo neto saliente por sus seis caras, por unidad de volumen, es la divergencia de v. A la derecha, en un flujo estacionario el caudal másico es el mismo a través de cada sección de un tubo de corriente.
A la izquierda, el balance de masa de una caja fija: el flujo neto saliente por sus seis caras, por unidad de volumen, es la divergencia de ρv\rho\vect v. A la derecha, en un flujo estacionario el caudal másico es el mismo a través de cada sección de un tubo de corriente.

Corolario 2.10 (Flujo incompresible; tubos de corriente)

Un flujo es incompresible cuando la densidad de cada partícula se mantiene constante, Dρ/Dt=0\mathrm D\rho/\mathrm Dt = 0, lo que equivale a

divv=0.\operatorname{div}\vect v = 0 .

Un líquido homogéneo fluye de forma incompresible, y también un gas a velocidades muy inferiores a la del sonido (véase Capítulo 7). En un flujo estacionario el caudal másico ρvS\rho vS es el mismo a través de cada sección de un tubo de corriente; si además el flujo es incompresible y ρ\rho uniforme, vSvS es constante: el fluido se acelera donde el tubo se estrecha.

Demostración. Dρ/Dt=ρdivv\mathrm D\rho/\mathrm Dt = -\rho\operatorname{div}\vect v. Para el tubo, aplíquese el balance integrado al volumen comprendido entre dos secciones: no hay flujo por la superficie lateral (la velocidad es tangente) y el régimen es estacionario, luego lo que entra iguala a lo que sale.

Ejemplo 2.11 (La sangre)

La aorta (de sección 2.5cm22.5\,\mathrm{cm}^{2}) transporta 5L/min5\,\mathrm{L}/\mathrm{min} =83cm3/s= 83\,\mathrm{cm}^{3}/\mathrm{s}: v=33cm/sv = 33\,\mathrm{cm}/\mathrm{s}. El mismo caudal pasa por unos 101010^{10} capilares de sección total 2500cm22500\,\mathrm{cm}^{2}: v=0.3mm/sv = 0.3\,\mathrm{mm}/\mathrm{s} — bastante lento para los intercambios a través de las paredes, que son la razón de ser de la ramificación.

2.4 Vorticidad y flujos irrotacionales

Definición 2.12 (Vorticidad)

La vorticidad de un flujo es el campo vectorial

ω=curlv,curlA=(AzyAyz,  AxzAzx,  AyxAxy).\vect\omega = \operatorname{\vect{curl}}\vect v , \qquad \operatorname{\vect{curl}}\vect A = \Bigl(\frac{\partial A_z}{\partial y} - \frac{\partial A_y}{\partial z},\; \frac{\partial A_x}{\partial z} - \frac{\partial A_z} {\partial x},\; \frac{\partial A_y}{\partial x} - \frac{\partial A_x}{\partial y}\Bigr) .

Un flujo es irrotacional allí donde ω=0\vect\omega = \vect 0.

Proposición 2.13 (La vorticidad mide la rotación local)

Un elemento pequeño de fluido centrado en MM gira, como un todo, con la velocidad angular Ω=12ω(M)\vect\Omega = \tfrac12\vect\omega(M). En una rotación como sólido rígido v=ΩOM\vect v = \vect\Omega\wedge\vect{OM} la vorticidad es uniforme, ω=2Ω\vect\omega = 2\vect\Omega; en el flujo plano de cizalla v=kyex\vect v = ky\,\vect e_x vale ω=kez\vect\omega = -k\vect e_z: una ruedecilla de paletas soltada en un flujo de cizalla gira aunque las líneas de corriente sean rectas.

Demostración. Cerca de MM, v(M+ϵ)=v(M)+Gϵ\vect v(M + \vect\epsilon) = \vect v(M) + \mathbf G\vect\epsilon con G\mathbf G la matriz de jvi\partial_jv_i; sepárese en su parte simétrica (una deformación pura, que estira el elemento sin girarlo en conjunto) y su parte antisimétrica, cuya matriz 12(jviivj)\tfrac12(\partial_jv_i - \partial_iv_j) es la de ϵ12ωϵ\vect\epsilon \mapsto \tfrac12\vect\omega\wedge\vect\epsilon (Teorema 1.2). Rotación como sólido en torno a ez\vect e_z: v=Ω(y,x,0)\vect v = \Omega(-y, x, 0), ωz=xvyyvx=2Ω\omega_z = \partial_xv_y - \partial_yv_x = 2\Omega. Cizalla: ωz=y(ky)=k\omega_z = -\partial_y(ky) = -k.

Definición 2.14 (Potencial de velocidades; circulación)

En un flujo irrotacional la velocidad deriva de un potencial de velocidades: v=gradφ\vect v = \operatorname{\vect{grad}}\varphi (la circulación de v\vect v a lo largo de un camino depende entonces solo de sus extremos, ABv ⁣dl=φ(B)φ(A)\int_A^B\vect v \cdot\dd\vect l = \varphi(B) - \varphi(A)). Si además el flujo es incompresible, Δφ=divgradφ=0\Delta\varphi = \operatorname{div}\operatorname{\vect{grad}} \varphi = 0: el potencial obedece la ecuación de Laplace, exactamente igual que el potencial electrostático en una región sin cargas — un flujo potencial. La circulación de la velocidad a lo largo de una curva cerrada CC es Γ=Cv ⁣dl\Gamma = \oint_C\vect v\cdot\dd\vect l.

Demostración. Que un campo de rotacional nulo en una región simplemente conexa sea un gradiente, y que un gradiente tenga rotacional nulo (yzφ=zyφ\partial_y\partial_z\varphi = \partial_z\partial_y\varphi), corresponden al volumen de matemáticas del Año 2; la ecuación de Laplace se sigue de divv=0\operatorname{div}\vect v = 0.

Ejemplo 2.15 (Tres flujos planos)

(i) Flujo de estancamiento v=k(x,y)\vect v = k(x, -y): divv=0\operatorname{div}\vect v = 0, ω=0\vect\omega = \vect 0, φ=12k(x2y2)\varphi = \tfrac12k(x^2 - y^2); líneas de corriente xy=xy = const., hipérbolas — un chorro que choca contra una pared, cerca del eje. (ii) Vórtice puntual v=Γ2πreθ\vect v = \dfrac{\Gamma}{2\pi r}\vect e_\theta (en coordenadas polares planas): incompresible e irrotacional en todas partes salvo en r=0r = 0, φ=Γθ/2π\varphi = \Gamma\theta/2\pi (multivaluado), y la circulación sobre cualquier círculo en torno al centro vale Γ\Gamma — la vorticidad está concentrada en el eje. (iii) Vórtice de Rankine: un núcleo r<ar < a en rotación como sólido vθ=Ωrv_\theta = \Omega r (vorticidad 2Ω2\Omega) empalmado con el vórtice puntual vθ=Ωa2/rv_\theta = \Omega a^2/r en el exterior: el modelo de un tornado, del remolino de una bañera o del que queda tras la pila del puente — la velocidad es máxima en el borde del núcleo.

Tres flujos planos. A la izquierda, el flujo de estancamiento, irrotacional, con líneas de corriente hiperbólicas. En el centro, el vórtice puntual: líneas de corriente circulares y, sin embargo, vorticidad nula fuera del centro. A la derecha, el perfil del vórtice de Rankine, con rotación como sólido en el núcleo y un vórtice puntual en el exterior.
Tres flujos planos. A la izquierda, el flujo de estancamiento, irrotacional, con líneas de corriente hiperbólicas. En el centro, el vórtice puntual: líneas de corriente circulares y, sin embargo, vorticidad nula fuera del centro. A la derecha, el perfil del vórtice de Rankine, con rotación como sólido en el núcleo y un vórtice puntual en el exterior.

Observación 2.16 (Por qué importa la circulación)

Para un lazo plano pequeño de área  ⁣dS\dd S y normal n\vect n, la circulación vale v ⁣dl=ωn ⁣dS\oint\vect v\cdot\dd\vect l = \vect\omega\cdot\vect n\,\dd S — la vorticidad es la circulación por unidad de área, exactamente como la divergencia es el flujo por unidad de volumen. (Compruébese sobre la rotación como sólido: 2πrΩr=2Ωπr22\pi r\cdot\Omega r = 2\Omega\cdot\pi r^2.) Sumada sobre una superficie, esto da el teorema de Stokes, Cv ⁣dl=Scurlvn ⁣dS\oint_C\vect v\cdot\dd\vect l = \iint_S \operatorname{\vect{curl}}\vect v\cdot\vect n\,\dd S, enunciado en Capítulo 11. Un ala genera sustentación porque lleva una circulación en torno a sí; una línea de vórtice persiste en un fluido perfecto (teorema de Kelvin) — y por eso el remolino de detrás de la pila dura tanto.

Método 2.17 (Cómo leer un campo de velocidades)

Dado v(M,t)\vect v(M, t): (1) ¿es estacionario? (tv=0\partial_t\vect v = \vect 0); (2) ¿es incompresible? (divv=0\operatorname{div}\vect v = 0); (3) ¿es irrotacional? (curlv=0\operatorname{\vect{curl}}\vect v = \vect 0; si lo es, hállese φ\varphi); (4) hállense las líneas de corriente a partir de  ⁣dx/vx= ⁣dy/vy= ⁣dz/vz\dd x/v_x = \dd y/v_y = \dd z/ v_z a tt fijo, y las trayectorias a partir de  ⁣dr/ ⁣dt=v(r,t)\dd\vect r/\dd t = \vect v(\vect r, t); (5) calcúlese la aceleración tv+(vgrad)v\partial_t\vect v + (\vect v\cdot \operatorname{\vect{grad}})\vect v. En coordenadas polares planas (r,θ)(r, \theta), divv=1r(rvr)r+1rvθθ\operatorname{div}\vect v = \dfrac1r\dfrac{\partial(rv_r)}{\partial r} + \dfrac1r \dfrac{\partial v_\theta}{\partial\theta} y ωz=1r(rvθ)r1rvrθ\omega_z = \dfrac1r\dfrac{\partial(rv_\theta)} {\partial r} - \dfrac1r\dfrac{\partial v_r}{\partial\theta} (admitidas aquí; se deducen en Capítulo 11).

2.5 Ejercicios

Ejercicio 2.1

El recorrido libre medio del aire vale 70nm70\,\mathrm{nm} a 1bar1\,\mathrm{bar} y varía como 1/P1/P. ¿Es válida la descripción como medio continuo para (a) el aire en torno a un coche, (b) el aire en un microcanal de 1µm1\,\text{µ}\mathrm{m} a 1bar1\,\mathrm{bar}, (c) el aire a 100km100\,\mathrm{km} de altitud (P3×107barP \approx 3 \times 10^{-7}\,\mathrm{bar}) en torno a un satélite de 1m1\,\mathrm{m}? Dese /L\ell/L en cada caso.

Solución

Solución de Ejercicio 2.1.

(a) /L7×108\ell/L \sim 7 \times 10^{-8}: medio continuo. (b) 0.070.07: caso límite — las leyes del continuo empiezan a fallar (deslizamiento en las paredes). (c) =70×109/3×107=0.23m\ell = 70 \times 10^{-9} /3 \times 10^{-7} = 0.23\,\mathrm{m}, /L=0.23\ell/L = 0.23: no hay medio continuo; el gas es una lluvia de moléculas independientes (régimen molecular libre).

Ejercicio 2.2

Flujo unidimensional estacionario v(x)=v0(1+x/L)v(x) = v_0(1 + x/L) para 0xL0 \le x \le L. Aceleración de una partícula en xx; tiempo que tarda en ir de 00 a LL; compárese con L/v0L/v_0 y L/2v0L/2v_0.

Solución

Solución de Ejercicio 2.2.

a=v ⁣dv/ ⁣dx=v02(1+x/L)/La = v\,\dd v/\dd x = v_0^2(1 + x/L)/L (de v02/Lv_0^2/L a 2v02/L2v_0^2/L).  ⁣dx/ ⁣dt=v0(1+x/L)\dd x/\dd t = v_0(1 + x/L): t=(L/v0)ln2=0.69L/v0t = (L/v_0)\ln2 = 0.69\,L/v_0, entre L/2v0L/2v_0 y L/v0L/v_0.

Ejercicio 2.3

En una tubería de radio RR el perfil de velocidades es v(r)=vmax(1r2/R2)v(r) = v_{\max}(1 - r^2/R^2) a lo largo del eje. Caudal volumétrico; celeridad media v=DV/πR2\langle v\rangle = D_V/\pi R^2; valores para R=5mmR = 5\,\mathrm{mm}, vmax=0.4m/sv_{\max} = 0.4\,\mathrm{m}/\mathrm{s}.

Solución

Solución de Ejercicio 2.3.

DV=0Rvmax(1r2/R2)2πr ⁣dr=12πR2vmaxD_V = \int_0^Rv_{\max}(1 - r^2/R^2)\,2\pi r\,\dd r = \tfrac12\pi R^2v_{\max}; v=vmax/2=0.20m/s\langle v\rangle = v_{\max}/2 = 0.20\,\mathrm{m}/\mathrm{s}; DV=1.6×105m3/s=16mL/sD_V = 1.6 \times 10^{-5}\,\mathrm{m}^{3}/\mathrm{s} = 16\,\mathrm{mL}/\mathrm{s}.

Ejercicio 2.4

Para cada campo, dígase si el flujo es incompresible y si es irrotacional, y dese la vorticidad: (a) v=(ax,ay,0)\vect v = (ax, -ay, 0); (b) v=(Ωy,Ωx,0)\vect v = (-\Omega y, \Omega x, 0); (c) v=(ax,ay,0)\vect v = (ax, ay, 0); (d) v=(ky,0,0)\vect v = (ky, 0, 0); (e) v=(ax,ay,bz)\vect v = (ax, -ay, bz) — ¿para qué bb es incompresible?

Solución

Solución de Ejercicio 2.4.

(a) div=0\operatorname{div} = 0, ω=0\vect\omega = \vect 0: incompresible e irrotacional. (b) div=0\operatorname{div} = 0, ωz=2Ω\omega_z = 2\Omega: rotación como sólido. (c) div=2a0\operatorname{div} = 2a \ne 0, irrotacional (una fuente plana). (d) div=0\operatorname{div} = 0, ωz=k\omega_z = -k: cizalla. (e) div=b\operatorname{div} = b: incompresible si y solo si b=0b = 0; irrotacional para todo bb.

Ejercicio 2.5 ★★

Una manguera de jardín de diámetro interior 12mm12\,\mathrm{mm} llena un cubo de 10L10\,\mathrm{L} en 50s50\,\mathrm{s}. Celeridad en la manguera; celeridad en una boquilla de 3mm3\,\mathrm{mm}; si la boquilla se estrecha a lo largo de 4cm4\,\mathrm{cm}, estímese la aceleración media de una partícula que la atraviesa (en unidades de gg); tiempo que pasa en la boquilla.

Solución

Solución de Ejercicio 2.5.

DV=0.20L/sD_V = 0.20\,\mathrm{L}/\mathrm{s}. Manguera S=1.13cm2S = 1.13\,\mathrm{cm}^{2}: v=1.8m/sv = 1.8\,\mathrm{m}/\mathrm{s}; boquilla S=7.1mm2S = 7.1\,\mathrm{mm}^{2}: v=28m/sv = 28\,\mathrm{m}/\mathrm{s}. Aceleración media Δ(v2/2)/L=(2821.82)/0.081×104m/s21000g\Delta(v^2/2)/L = (28^2 - 1.8^2)/0.08 \approx 1 \times 10^{4}\,\mathrm{m}/\mathrm{s}^{2} \approx 1000\,g. Tiempo  ⁣dx/v\approx \int\dd x/v: unos milisegundos (4ms4\,\mathrm{ms} para un perfil de celeridad lineal).

Ejercicio 2.6 ★★

Flujo plano no estacionario v=(U,Vcosωt)\vect v = (U, V\cos\omega t) con UU, VV, ω\omega constantes. (a) Líneas de corriente en el instante tt. (b) Trayectoria de la partícula que está en el origen en t=0t = 0. (c) Dibújense ambas para V=UV = U en t=0t = 0 y en t=π/2ωt = \pi/2\omega; ¿por qué difieren? (d) Aceleración de una partícula.

Solución

Solución de Ejercicio 2.6.

(a)  ⁣dy/ ⁣dx=(V/U)cosωt\dd y/\dd x = (V/U)\cos\omega t: rectas de pendiente (V/U)cosωt(V/U)\cos\omega t, todas paralelas, que se inclinan con el tiempo. (b) x=Utx = Ut, y=(V/ω)sinωty = (V/\omega)\sin\omega t: la sinusoide y=(V/ω)sin(ωx/U)y = (V/\omega)\sin(\omega x/U). (c) En t=0t = 0 las líneas de corriente están a 4545{}^{\circ} y la trayectoria sale del origen a 4545{}^{\circ}; en t=π/2ωt = \pi/2\omega las líneas de corriente son horizontales mientras que la trayectoria es una sinusoide fija. Difieren porque el flujo no es estacionario: una línea de corriente es una instantánea, una trayectoria es una historia. (d) a=tv=(0,Vωsinωt)\vect a = \partial_t\vect v = (0, -V\omega\sin\omega t) — el campo es uniforme, luego el término convectivo se anula.

Ejercicio 2.7 ★★

Un río de 50m50\,\mathrm{m} de ancho y 2.0m2.0\,\mathrm{m} de profundidad fluye a 1.2m/s1.2\,\mathrm{m}/\mathrm{s}. (a) Caudal. (b) Entra en un desfiladero de 15m15\,\mathrm{m} de ancho y 4.0m4.0\,\mathrm{m} de profundidad: celeridad. (c) Un afluente añade 30m3/s30\,\mathrm{m}^{3}/\mathrm{s}; aguas abajo el río tiene 60m60\,\mathrm{m} de ancho y 2.5m2.5\,\mathrm{m} de profundidad: celeridad. (d) Un gas caliente fluye en régimen estacionario por una tubería de sección constante; entre dos puntos su temperatura se duplica a presión constante: ¿qué le ocurre a su celeridad?

Solución

Solución de Ejercicio 2.7.

(a) DV=50×2×1.2=120m3/sD_V = 50 \times 2 \times 1.2 = 120\,\mathrm{m}^{3}/\mathrm{s}. (b) 120/60=2.0m/s120/60 = 2.0\,\mathrm{m}/\mathrm{s}. (c) 150/150=1.0m/s150/150 = 1.0\,\mathrm{m}/\mathrm{s}. (d) ρ1/T\rho \propto 1/T a PP constante se reduce a la mitad; ρvS\rho vS se conserva: la celeridad se duplica.

Ejercicio 2.8 ★★

Vórtice de Rankine. vθ=Ωrv_\theta = \Omega r para r<ar < a, vθ=Ωa2/rv_\theta = \Omega a^2/r para r>ar > a. (a) Vorticidad en cada región (úsese la fórmula polar). (b) Circulación sobre un círculo de radio rr centrado en el eje, en ambas regiones; relaciónese con el flujo de la vorticidad. (c) Un tornado: a=50ma = 50\,\mathrm{m}, vmax=70m/sv_{\max} = 70\,\mathrm{m}/\mathrm{s}: Ω\Omega, Γ\Gamma y la celeridad a 500m500\,\mathrm{m}. (d) ¿Es irrotacional el flujo del núcleo? ¿Y el exterior? ¿Dónde está físicamente la “rotación”?

Solución

Solución de Ejercicio 2.8.

(a) ωz=1r ⁣d(rvθ) ⁣dr\omega_z = \frac1r\frac{\dd(rv_\theta)}{\dd r}: 2Ω2\Omega dentro, 00 fuera. (b) Dentro, Γ=2πrΩr=2Ωπr2\Gamma = 2\pi r\cdot\Omega r = 2\Omega\cdot\pi r^2, el flujo de la vorticidad; fuera, Γ=2πΩa2\Gamma = 2\pi\Omega a^2, constante, el flujo de todo el núcleo. (c) Ω=vmax/a=1.4rad/s\Omega = v_{\max}/a = 1.4\,\mathrm{rad}/\mathrm{s}; Γ=2πavmax=2.2×104m2/s\Gamma = 2\pi av_{\max} = 2.2 \times 10^{4}\,\mathrm{m}^{2}/\mathrm{s}; a 500m500\,\mathrm{m}: Ωa2/r=7m/s\Omega a^2/r = 7\,\mathrm{m}/\mathrm{s}. (d) El núcleo gira como un sólido; fuera, el flujo es irrotacional: las partículas dan vueltas en torno al eje sin girar sobre sí mismas. La vorticidad — la rotación — vive en el núcleo.

Ejercicio 2.9 ★★

Flujo de estancamiento v=k(x,y)\vect v = k(x, -y). (a) Potencial φ\varphi y comprobación de Δφ=0\Delta\varphi = 0. (b) Líneas de corriente; la función ψ=kxy\psi = kxy es constante sobre ellas — muéstrese que vx=ψ/yv_x = \partial\psi/\partial y, vy=ψ/xv_y = -\partial\psi/ \partial x (la función de corriente). (c) Campo de aceleraciones; muéstrese que vale grad(v2/2)\operatorname{\vect{grad}}(v^2/2) e interprétese. (d) Trayectoria de una partícula que parte de (x0,y0)(x_0, y_0): x=x0ektx = x_0\eu^{kt}, y=y0ekty = y_0\eu^{-kt}.

Solución

Solución de Ejercicio 2.9.

(a) φ=12k(x2y2)\varphi = \tfrac12k(x^2 - y^2), Δφ=kk=0\Delta\varphi = k - k = 0. (b)  ⁣dx/kx= ⁣dy/(ky)\dd x/ kx = \dd y/(-ky): xy=xy = const.; yψ=kx=vx\partial_y\psi = kx = v_x, xψ=ky=vy-\partial_x\psi = -ky = v_y. (c) a=(vxx+vyy)v=k2(x,y)=grad(12k2(x2+y2))=grad(v2/2)\vect a = (v_x\partial_x + v_y\partial_y)\vect v = k^2(x, y) = \operatorname{\vect{grad}}\bigl(\tfrac12k^2(x^2 + y^2)\bigr) = \operatorname{\vect{grad}} (v^2/2): en un flujo irrotacional estacionario la aceleración es el gradiente de la energía cinética por unidad de masa; apunta alejándose del punto de estancamiento — una partícula se frena al acercarse a la pared según yy y se acelera al alejarse según xx. (d) x˙=kx\dot x = kx, y˙=ky\dot y = -ky: x=x0ektx = x_0\eu^{kt}, y=y0ekty = y_0\eu^{-kt} (xyxy constante).

Ejercicio 2.10 ★★★

Un flujo en expansión. Flujo unidimensional v(x,t)=x/(t+τ)v(x, t) = x/(t + \tau) para t>0t > 0, con τ>0\tau > 0 constante. (a) Trayectorias: muéstrese que x(t)=x0(1+t/τ)x(t) = x_0(1 + t/\tau); aceleración de una partícula, tanto a partir de la trayectoria como de la derivada material. (b) La densidad es uniforme en cada instante, ρ(t)\rho(t): úsese la conservación de la masa para hallar ρ(t)\rho(t). (c) Compruébese siguiendo la masa comprendida entre dos partículas. (d) ¿Por qué es este campo un modelo unidimensional de un universo en expansión uniforme?

Solución

Solución de Ejercicio 2.10.

(a)  ⁣dx/ ⁣dt=x/(t+τ)\dd x/\dd t = x/(t + \tau) da x=x0(1+t/τ)x = x_0(1 + t/\tau), x¨=0\ddot x = 0; derivada material: tv+vxv=x/(t+τ)2+x/(t+τ)2=0\partial_tv + v\partial_xv = -x/(t + \tau)^2 + x/(t + \tau)^2 = 0. (b) tρ+x(ρv)=ρ˙+ρ/(t+τ)=0\partial_t\rho + \partial_x(\rho v) = \dot\rho + \rho/(t + \tau) = 0: ρ=ρ0τ/(t+τ)\rho = \rho_0\tau/(t + \tau). (c) La separación de dos partículas crece como (1+t/τ)(1 + t/\tau) y la masa entre ellas es fija: ρ1/(1+t/τ)\rho \propto 1/(1 + t/\tau). (d) v=Hxv = Hx con H=1/(t+τ)H = 1/(t + \tau): una ley de Hubble, expansión libre (sin aceleración), con la densidad decreciendo como el inverso del factor de escala.

Ejercicio 2.11 ★★★

Una fuente en una corriente. Superpóngase un flujo uniforme v=Uex\vect v = U\vect e_x y una fuente plana v=q2πrer\vect v = \dfrac{q}{2\pi r}\vect e_r (qq es el caudal volumétrico por unidad de longitud según zz). (a) Potencial φ=Ux+(q/2π)lnr\varphi = Ux + (q/2\pi)\ln r; compruébese que ambos flujos son incompresibles e irrotacionales (salvo en el origen). (b) Punto de estancamiento. (c) Función de corriente ψ=Uy+qθ/2π\psi = Uy + q\theta/2\pi; la línea de corriente que pasa por el punto de estancamiento es rsinθ=(q/2πU)(πθ)r\sin\theta = (q/2\pi U)(\pi - \theta): dibújese. (d) Muéstrese que muy aguas abajo esa línea de corriente está en y=±q/2Uy = \pm q/2U: el flujo exterior a ella es el flujo alrededor de un semicuerpo romo de anchura q/Uq/U.

Solución

Solución de Ejercicio 2.11.

(a) v=gradφ=(U+qx/2πr2,qy/2πr2)\vect v = \operatorname{\vect{grad}}\varphi = (U + qx/2\pi r^2, qy/2\pi r^2); la fuente tiene vr=q/2πrv_r = q/2\pi r, div=1r ⁣d(rvr)/ ⁣dr=0\operatorname{div} = \frac1r\dd(rv_r)/\dd r = 0 y no depende de θ\theta, luego ωz=0\omega_z = 0; las sumas de campos así son de nuevo incompresibles e irrotacionales. (b) Sobre el semieje xx negativo, Uq/2πr=0U - q/2\pi r = 0: xs=q/2πUx_s = -q/2\pi U. (c) En el punto de estancamiento θ=π\theta = \pi, y=0y = 0: ψ=q/2\psi = q/2; la línea de corriente Ursinθ+qθ/2π=q/2Ur\sin\theta + q\theta/2\pi = q/2, es decir rsinθ=(q/2πU)(πθ)r\sin\theta = (q/2\pi U)(\pi - \theta) — una curva que sale de xsx_s y se abre aguas abajo. (d) θ0\theta \to 0: yq/2Uy \to q/2U; θ2π\theta \to 2\pi (rama inferior): yq/2Uy \to -q/2U. El fluido interior a esa línea de corriente procede de la fuente; fuera de ella, la corriente uniforme rodea un semicuerpo de anchura asintótica q/Uq/U.

Ejercicio 2.12 ★★★

Vórtice y sumidero. Flujo plano v=q2πrer+Γ2πreθ\vect v = -\dfrac{q}{2\pi r}\vect e_r + \dfrac{\Gamma}{2\pi r}\vect e_\theta para r>ar > a. (a) Muéstrese que es incompresible e irrotacional. (b) Líneas de corriente: espirales logarítmicas r=r0eqθ/Γr = r_0\eu^{-q\theta/\Gamma}. (c) Aceleración de una partícula: muéstrese que es puramente radial, v2/rer-v^2/r\,\vect e_r. (d) Tiempo que tarda una partícula en describir la espiral desde r0r_0 hasta aa; número de vueltas que da. (e) Valores: q=1.5×103m2/sq = 1.5 \times 10^{-3}\,\mathrm{m}^{2}/\mathrm{s}, Γ=3×102m2/s\Gamma = 3 \times 10^{-2}\,\mathrm{m}^{2}/\mathrm{s}, r0=0.5mr_0 = 0.5\,\mathrm{m}, a=2cma = 2\,\mathrm{cm}.

Solución

Solución de Ejercicio 2.12.

(a) divv=1rr(q/2π)=0\operatorname{div}\vect v = \frac1r\partial_r(-q/2\pi) = 0; ωz=1rr(Γ/2π)=0\omega_z = \frac1r\partial_r(\Gamma/2\pi) = 0. (b)  ⁣dr/vr=r ⁣dθ/vθ\dd r/v_r = r\,\dd\theta/v_\theta da  ⁣dr/r=(q/Γ) ⁣dθ\dd r/r = -(q/\Gamma)\dd\theta. (c) ar=vrrvrvθ2/r=(q2+Γ2)/4π2r3=v2/ra_r = v_r\partial_rv_r - v_\theta^2/r = -(q^2 + \Gamma^2)/4\pi^2r^3 = -v^2/r; aθ=vrrvθ+vrvθ/r=0a_\theta = v_r\partial_rv_\theta + v_rv_\theta/r = 0. (d)  ⁣dr/ ⁣dt=q/2πr\dd r/\dd t = -q/2\pi r: t=π(r02a2)/qt = \pi(r_0^2 - a^2)/q; vueltas =(Γ/2πq)ln(r0/a)= (\Gamma/2\pi q)\ln(r_0/a). (e) t=520s9mint = 520\,\mathrm{s} \approx 9\,\mathrm{min}; 1010 vueltas.

2.6 Problema: del desagüe al huracán

Problema 2.1

Problema de fin de semana — la misma cinemática a tres escalas: el remolino sobre un desagüe, un tornado y un ciclón visto desde un satélite

Parte I — La bañera. Una bañera se vacía por un orificio de radio a=2cma = 2\,\mathrm{cm}; el agua, de profundidad h=20cmh = 20\,\mathrm{cm}, sale con el caudal volumétrico DV=0.30L/sD_V = 0.30\,\mathrm{L}/\mathrm{s}. Lejos del orificio el flujo se modela como plano (independiente de la profundidad): v=vr(r)er+vθ(r)eθ\vect v = v_r(r)\vect e_r + v_\theta(r)\vect e_\theta.

  1. Usando la conservación de la masa a través de un cilindro de radio rr y altura hh, muéstrese que vr=q/2πrv_r = -q/2\pi r con q=DV/hq = D_V/h; valor de qq.
  2. Antes de quitar el tapón se había removido el agua, que giraba a 1.0cm/s1.0\,\mathrm{cm}/\mathrm{s} en r=50cmr = 50\,\mathrm{cm}. Admitiendo que cada anillo de agua conserva su momento angular al desplazarse hacia dentro (rvθrv_\theta constante), hállense vθ(r)v_\theta(r) y la circulación Γ\Gamma.
  3. Compruébese que el flujo es incompresible e irrotacional para r>ar > a.
  4. Líneas de corriente: muéstrese que son espirales logarítmicas y hállese cuántas vueltas da una partícula entre r=50cmr = 50\,\mathrm{cm} y r=ar = a.
  5. Tiempo que tarda una partícula en ir desde 50cm50\,\mathrm{cm} hasta el orificio, y su celeridad azimutal al llegar.
  6. Aceleración de una partícula en rr; compárese con gg en r=ar = a. (La superficie libre se hunde allí donde la presión debe equilibrarla — el embudo; capítulo siguiente.)
  7. La rotación de la Tierra da a todo fluido en reposo en la bañera una rotación “de fondo” a la tasa vertical local Ωsinλ5×105rad/s\Omega_\oplus\sin\lambda \approx 5 \times 10^{-5}\,\mathrm{rad}/\mathrm{s}. Compárese la celeridad azimutal que daría a 50cm50\,\mathrm{cm} con el 1cm/s1\,\mathrm{cm}/\mathrm{s} del removido, y dígase si el hemisferio decide el sentido de giro de una bañera que se vacía.

Parte II — El tornado como vórtice de Rankine. Radio del núcleo a=60ma = 60\,\mathrm{m}, viento máximo vmax=80m/sv_{\max} = 80\,\mathrm{m}/\mathrm{s}, densidad del aire ρ=1.2kg/m3\rho = 1.2\,\mathrm{kg}/\mathrm{m}^{3}.

  1. Escríbase vθ(r)v_\theta(r) en ambas regiones y dense Ω\Omega y Γ\Gamma.
  2. Vorticidad en el núcleo y fuera de él; represéntense vθv_\theta y ωz\omega_z frente a rr.
  3. Aceleración de una partícula de aire en el borde del núcleo; compárese con gg.
  4. Muéstrese que una ruedecilla de paletas colocada fuera del núcleo deriva alrededor del tornado sin girar sobre su propio eje, mientras que una colocada en el núcleo da una vuelta sobre sí misma por cada vuelta en torno al eje.
  5. Energía cinética del núcleo, por metro de altura.
  6. Energía cinética del flujo exterior por metro de altura, entre aa y R=5kmR = 5\,\mathrm{km}; ¿por qué depende de RR solo logarítmicamente?
  7. Con una altura de 1km1\,\mathrm{km}, energía cinética total; compárese con la energía liberada por 1t1\,\mathrm{t} de TNT (4.2GJ4.2\,\mathrm{GJ}).
  8. En un tornado el aire también fluye hacia dentro y hacia arriba. Si el núcleo se alimenta de una entrada radial a 5m/s5\,\mathrm{m}/\mathrm{s} por su superficie lateral a lo largo de los primeros 300m300\,\mathrm{m} de altura, ¿qué celeridad vertical debe alcanzar el aire en lo alto de esa capa (uniforme sobre la sección del núcleo)?

Parte III — El ciclón. Un ciclón tropical se modela como un vórtice de Rankine de radio de núcleo a=40kma = 40\,\mathrm{km} y vmax=50m/sv_{\max} = 50\,\mathrm{m}/\mathrm{s}. La rotación de la Tierra se deja sentir a través de la vorticidad planetaria f=2Ωsinλ5×105s1f = 2\Omega_\oplus\sin\lambda \approx 5 \times 10^{-5}\,\mathrm{s}^{-1} a la latitud 2020{}^{\circ}: un anillo de aire en reposo respecto del suelo ya lleva, en el sistema inercial, la celeridad azimutal fr/2fr/2 en torno a cualquier eje vertical.

  1. Vorticidad del núcleo y circulación Γ\Gamma del ciclón; compárese la vorticidad del núcleo con ff.
  2. Un anillo de aire de radio r0=500kmr_0 = 500\,\mathrm{km}, inicialmente en reposo respecto del suelo, es atraído hacia el centro por la baja presión. Admitiendo que se conserva su momento angular por unidad de masa respecto del eje, r(vθ+fr/2)r(v_\theta + fr/2), hállese su velocidad vθv_\theta respecto del suelo en r=100kmr = 100\,\mathrm{km} y en r=40kmr = 40\,\mathrm{km}.
  3. ¿En qué sentido gira el ciclón en el hemisferio norte, y por qué? ¿Por qué el sentido está fijado en un ciclón y no en una bañera?
  4. ¿Por qué no puede arrastrarse así el aire hasta el eje mismo — qué limita el viento y deja un ojo?
  5. El aire converge hacia el ciclón a 5m/s5\,\mathrm{m}/\mathrm{s} a través de un cilindro de radio 500km500\,\mathrm{km} y altura 1km1\,\mathrm{km}: caudal másico arrastrado (ρ=1.2kg/m3\rho = 1.2\,\mathrm{kg}/\mathrm{m}^{3}).
  6. El ciclón se desplaza en bloque hacia el oeste a 6m/s6\,\mathrm{m}/\mathrm{s}. Escríbase el campo de velocidades euleriano visto desde el suelo en función del campo v0\vect v_0 del vórtice en su propio sistema; ¿es estacionario el flujo en el sistema del suelo? ¿En qué lado de la tormenta son más fuertes los vientos?

Parte IV — Caudales y un resumen.

  1. Llueve a 50mm/h50\,\mathrm{mm}/\mathrm{h} sobre una cuenca de 10km210\,\mathrm{km}^{2} y toda el agua llega a un río de sección 40m240\,\mathrm{m}^{2}: celeridad media del río en crecida.
  2. La sangre: la aorta (2.5cm22.5\,\mathrm{cm}^{2}) transporta 5.0L/min5.0\,\mathrm{L}/\mathrm{min}; los capilares tienen una sección total de 0.25m20.25\,\mathrm{m}^{2}: celeridades en cada caso; tiempo que tarda la sangre en recorrer un capilar de 1mm1\,\mathrm{mm}.
  3. Un depósito de sección SS se vacía por un orificio de sección ss en el fondo; la celeridad de salida es 2gh\sqrt{2gh} (hh es el nivel del agua — admitido, se deduce en el capítulo siguiente). Escríbase la conservación de la masa y hállense h(t)h(t) y el tiempo de vaciado para S=1m2S = 1\,\mathrm{m}^{2}, s=1cm2s = 1\,\mathrm{cm}^{2}, h0=1mh_0 = 1\,\mathrm{m}.
  4. Reúnanse los tres flujos de las partes I–III en una tabla: radio, celeridad, vorticidad del núcleo, circulación — y la magnitud que, conservada durante la convergencia, los pone a girar.
Solución

Solución de Problema 2.1.

1. A través de un cilindro de radio rr: 2πrhvr=DV2\pi rh|v_r| = D_V, luego vr=q/2πrv_r = -q/2\pi r, q=DV/h=1.5×103m2/sq = D_V/h = 1.5 \times 10^{-3}\,\mathrm{m}^{2}/\mathrm{s}.

2. rvθ=0.5×0.01=5×103m2/srv_\theta = 0.5 \times 0.01 = 5 \times 10^{-3}\,\mathrm{m}^{2}/\mathrm{s}: vθ=Γ/2πrv_\theta = \Gamma/2\pi r con Γ=3.1×102m2/s\Gamma = 3.1 \times 10^{-2}\,\mathrm{m}^{2}/\mathrm{s}.

3. divv=1rr(rvr)=0\operatorname{div}\vect v = \frac1r\partial_r(rv_r) = 0, ωz=1rr(rvθ)=0\omega_z = \frac1r\partial_r(rv_\theta) = 0.

4.  ⁣dr/r=(q/Γ) ⁣dθ\dd r/r = -(q/\Gamma)\dd\theta: r=r0eqθ/Γr = r_0\eu^{-q\theta/\Gamma}; vueltas =(Γ/2πq)ln(r0/a)=3.3×3.211= (\Gamma/2\pi q)\ln(r_0/a) = 3.3 \times 3.2 \approx 11.

5. t=π(r02a2)/q=520s9mint = \pi(r_0^2 - a^2)/q = 520\,\mathrm{s} \approx 9\,\mathrm{min}; vθ(a)=0.25m/sv_\theta(a) = 0.25\,\mathrm{m}/\mathrm{s}.

6. a=(v2/r)er\vect a = -(v^2/r)\vect e_r; en r=ar = a: v2=0.252+0.0122v^2 = 0.25^2 + 0.012^2, a=3.1m/s20.3ga = 3.1\,\mathrm{m}/\mathrm{s}^{2} \approx 0.3g.

7. Ωsinλr=5×105×0.5=2.5×105m/s\Omega_\oplus\sin\lambda\,r = 5 \times 10^{-5} \times 0.5 = 2.5 \times 10^{-5}\,\mathrm{m}/\mathrm{s}, cuatrocientas veces menos que el removido: el hemisferio no decide nada en una bañera corriente (sí lo hace en un depósito dejado quieto un día).

8. vθ=Ωrv_\theta = \Omega r (r<ar < a), Ωa2/r\Omega a^2/r (r>ar > a); Ω=80/60=1.33rad/s\Omega = 80/60 = 1.33\,\mathrm{rad}/\mathrm{s}, Γ=2πavmax=3.0×104m2/s\Gamma = 2\pi av_{\max} = 3.0 \times 10^{4}\,\mathrm{m}^{2}/\mathrm{s}.

9. ωz=2Ω=2.7s1\omega_z = 2\Omega = 2.7\,\mathrm{s}^{-1} en el núcleo, 00 fuera; vθv_\theta crece linealmente hasta 80m/s80\,\mathrm{m}/\mathrm{s} y luego cae como 1/r1/r; ωz\omega_z es un escalón.

10. a=vmax2/a=6400/60=107m/s211ga = v_{\max}^2/a = 6400/60 = 107\,\mathrm{m}/\mathrm{s}^{2} \approx 11g, hacia dentro.

11. Fuera, la vorticidad es nula: un elemento pequeño se traslada a lo largo de su círculo sin girar sobre sí mismo — la ruedecilla de paletas conserva su orientación. En el núcleo (rotación como sólido) gira con el fluido, una vuelta por revolución.

12. Enuˊcleo=0a12ρΩ2r22πr ⁣dr=14πρΩ2a4=2.2×107J/mE_{\text{núcleo}} = \int_0^a\tfrac12\rho\Omega^2r^2\,2\pi r\,\dd r = \tfrac14\pi\rho\Omega^2a^4 = 2.2 \times 10^{7}\,\mathrm{J}/\mathrm{m}.

13. Eext=aR12ρ(Ωa2/r)22πr ⁣dr=πρΩ2a4ln(R/a)=8.7×107×4.4=3.8×108J/mE_{\text{ext}} = \int_a^R\tfrac12\rho(\Omega a^2/r)^2\,2\pi r\,\dd r = \pi\rho\Omega^2a^4\ln(R/a) = 8.7 \times 10^7 \times 4.4 = 3.8 \times 10^{8}\,\mathrm{J}/\mathrm{m}; el integrando ρv2r1/r\rho v^2r \propto 1/r, de ahí el logaritmo.

14. 4.0×108J/m×1000m=4×1011J4.0 \times 10^{8}\,\mathrm{J}/\mathrm{m} \times 1000\,\mathrm{m} = 4 \times 10^{11}\,\mathrm{J}: unas 100t100\,\mathrm{t} de TNT.

15. Entrada 2πa×300×5=5.7×105m3/s2\pi a \times 300 \times 5 = 5.7 \times 10^{5}\,\mathrm{m}^{3}/\mathrm{s} a través de πa2=1.1×104m2\pi a^2 = 1.1 \times 10^{4}\,\mathrm{m}^{2}: vz=50m/sv_z = 50\,\mathrm{m}/\mathrm{s}.

16. ω=2vmax/a=2.5×103s1=50f\omega = 2v_{\max}/a = 2.5 \times 10^{-3}\,\mathrm{s}^{-1} = 50f; Γ=2πavmax=1.3×107m2/s\Gamma = 2\pi av_{\max} = 1.3 \times 10^{7}\,\mathrm{m}^{2}/\mathrm{s}.

17. r0fr0/2=r(vθ+fr/2)r_0\cdot fr_0/2 = r(v_\theta + fr/2): vθ=f(r02r2)/2rv_\theta = f(r_0^2 - r^2)/2r; 60m/s60\,\mathrm{m}/\mathrm{s} a 100km100\,\mathrm{km} y 155m/s155\,\mathrm{m}/\mathrm{s} a 40km40\,\mathrm{km}.

18. f>0f > 0 en el hemisferio norte da vθ>0v_\theta > 0: sentido antihorario visto desde arriba (ciclónico). Sobre 500km500\,\mathrm{km} la rotación planetaria es la rotación inicial dominante; en una bañera lo es el removido.

19. El rozamiento con el mar retira momento angular, y el gradiente de presión no puede equilibrar la aceleración centrípeta de un anillo cada vez más rápido: el aire asciende antes de llegar al eje, en la pared del ojo, y deja un ojo en calma.

20. ρ2πrHv=1.2×2π×5×105×103×5=1.9×1010kg/s\rho\,2\pi rHv = 1.2 \times 2\pi \times 5 \times 10^5 \times 10^3 \times 5 = 1.9 \times 10^{10}\,\mathrm{kg}/\mathrm{s}.

21. v(M,t)=Uex+v0(M+Utex)\vect v(M, t) = -U\vect e_x + \vect v_0(M + Ut\,\vect e_x): el patrón es advectado y el campo no es estacionario en el sistema del suelo. La traslación se suma al remolino en el lado donde vθv_\theta apunta hacia el oeste — el lado norte (el de la derecha) de una trayectoria hacia el oeste.

22. 0.05×107/3600=139m3/s0.05 \times 10^7/3600 = 139\,\mathrm{m}^{3}/\mathrm{s}; v=139/40=3.5m/sv = 139/40 = 3.5\,\mathrm{m}/\mathrm{s}.

23. DV=83cm3/sD_V = 83\,\mathrm{cm}^{3}/\mathrm{s}: aorta 0.33m/s0.33\,\mathrm{m}/\mathrm{s}; capilares 0.33mm/s0.33\,\mathrm{mm}/\mathrm{s}; 3s3\,\mathrm{s} por milímetro.

24. S ⁣dh/ ⁣dt=s2ghS\,\dd h/\dd t = -s\sqrt{2gh}: h=h0sSg/2t\sqrt h = \sqrt{h_0} - \frac sS\sqrt{g/2}\,t, vacío en T=(S/s)2h0/g=104×0.45=4500s75minT = (S/s)\sqrt{2h_0/g} = 10^4 \times 0.45 = 4500\,\mathrm{s} \approx 75\,\mathrm{min}.

25. Bañera: r0.5mr \sim 0.5\,\mathrm{m}, v1cm/sv \sim 1\,\mathrm{cm}/\mathrm{s}, Γ=3×102m2/s\Gamma = 3 \times 10^{-2}\,\mathrm{m}^{2}/\mathrm{s}; tornado: 60m60\,\mathrm{m}, 80m/s80\,\mathrm{m}/\mathrm{s}, ω=2.7s1\omega = 2.7\,\mathrm{s}^{-1}, Γ=3×104m2/s\Gamma = 3 \times 10^{4}\,\mathrm{m}^{2}/\mathrm{s}; ciclón: 40km40\,\mathrm{km}, 50m/s50\,\mathrm{m}/\mathrm{s}, ω=2.5×103s1\omega = 2.5 \times 10^{-3}\,\mathrm{s}^{-1}, Γ=1.3×107m2/s\Gamma = 1.3 \times 10^{7}\,\mathrm{m}^{2}/\mathrm{s}. En los tres, el momento angular por unidad de masa rvθrv_\theta del fluido que converge se conserva: la convergencia los pone a girar.

Términos definidos en este capítulo

Ver los 393 términos del glosario