Anterior

La gota de lluvia que cae a través de una nube

Una gota de agua de lluvia cae a través de una nube de pequeñas gotitas. A medida que cae, incrementa su masa al chocar inelásticamente con las pequeñas gotitas. El problema consiste en determinar la posición x y velocidad v de la gota en función del tiempo t, conocida la masa inicial m0, la velocidad inicial v0 y la altura inicial x0 en el instante t=0.

 

 

 

 

La masa de la gota

Hemos de hacer una suposición acerca de la forma en que la masa de la gota se incrementa con el tiempo. Si la gota va absorbiendo las pequeñas gotitas que encuentra en su camino, entonces

dm dt (área)×(velocidad)π r 2 v= ρ n π r 2 v=k m 2/3 v

El valor de la constante de proporcionalidad k es

k= ρ n π ( ρ a 4 3 π ) 2/3

En general, supondremos que la razón del incremento de la masa de la gota con el tiempo es de la forma

dm dt =k m α v

Como la velocidad v=dx/dt

dm dt =k m α dx dt

Integramos esta ecuación con las condiciones iniciales para x=0, m=m0

m 0 m m α dm= 0 x k·dx m 1α m 0 1α =(1α)kx m= ( (1α)kx+ m 0 1α ) 1/(1α)

Esta ecuación nos proporciona la masa m de la gota en función de la posición x.

Ecuaciones del movimiento

Sobre la gota de masa m actúa una única fuerza que es su peso mg. La segunda ley de Newton aplicada a este sistema de masa variable se escribe

d(mv) dt =mgm dv dt +v dm dt =mg

Cuando g=0

Empezaremos por el caso más simple, aquél en el que la aceleración de la gravedad es cero. Podría ser el caso de un objeto que pasase a través de la materia interestelar.

Como la fuerza exterior es nula, el momento lineal se conserva, al aumentar la masa disminuye la velocidad de la gota

m0v0=mv

dx dt = m 0 v 0 m = m 0 v 0 ( (1α)kx+ m 0 1α ) 1/(1α)

Integramos

0 x ( (1α)kx+ m 0 1α ) 1/(1α) dx= 0 t m 0 v 0 dt 1 k(2α) { ( (1α)kx+ m 0 1α ) (2α)/(1α) m 0 (2α) }= m 0 v 0 t

Expresamos x en función del tiempo t

x= 1 k(1α) m 0 (1+α) { ( k(2α) m 0 (1+α) v 0 t+1 ) (1α)/(2α) 1 }

Calculamos ahora la velocidad v en función del tiempo t

v= m 0 v 0 m = m 0 v 0 ( (1α)kx+ m 0 1α ) 1/(1α) = v 0 ( k(2α) m 0 (1+α) v 0 t+1 ) 1/(2α)

Integrando de nuevo, obtenemos la posición x de la gota en función del tiempo t.

x= 0 t v 0 ( k(2α) m 0 (1+α) v 0 t+1 ) 1/(2α) dt = 1 k(1α) m 0 (1+α) { ( k(2α) m 0 (1+α) v 0 t+1 ) (1α)/(2α) 1 }

Cuando α=2/3 las expresiones de la masa m de la gota, la velocidad v y la posición x en función del tiempo t son:

m= ( 1 3 kx+ m 0 1/3 ) 3 v= v 0 ( 4 3 k m 0 1/3 v 0 t+1 ) 3/4 x= 3 m 0 1/3 k { ( 4 3 k m 0 1/3 v 0 t+1 ) 1/4 1 }

Cuando g≠0

Las ecuaciones que tenemos que resolver son

dm dt =k m α v m dv dt +v dm dt =mg

Con la notación

m ˙ = dm dt v ˙ = dv dt v ˙ ˙ = d v ˙ dt  = d 2 v d t 2

Las ecuaciones anteriores se escriben

m ˙ =k m α v m v ˙ +v m ˙ =mg

En general, la aceleración de la gota dv/dt no es constante, para que fuese constante se debería cumplir que

m ˙ m v=gc

donde c es una constante

Eliminado la derivada de m en las dos ecuaciones que describen el movimiento de la gota, obtenemos una ecuación diferencial de primer orden en v.

m ˙ =k m α v v ˙ +k m α1 v 2 =g m α+1 =k v 2 g v ˙

Derivamos respecto del tiempo

(α+1) m α m ˙ =k 2v v ˙ (g v ˙ ) v 2 ( v ˙ ˙ ) g v ˙ (α+1)kv=k 2v v ˙ (g v ˙ ) v 2 ( v ˙ ˙ ) g v ˙ v ˙ ˙ = (g v ˙ ) v ( (1α)g(3α) v ˙ )

Esta ecuación diferencial no tiene solución analítica conocida

La aceleración es constante cuando el término entre paréntesis es cero

v ˙ ˙ =0 v ˙ = 1α 3α g

Cuando α=2/3, la aceleración es constante e igual a 1/7 de la aceleración de la gravedad g

dv dt = 1 7 g

Procedimiento numérico

Se resuelve la ecuación diferencial de segundo orden por el procedimiento de Runge-Kutta

d 2 x dt =k m α1 v 2 +g   m 1α =(1α)kx+ m 0 1α

Con las condiciones iniciales t=0, v=0, m=m0

La masa inicial m0 en gramos es el producto de la densidad del agua 1.0 g/cm3 por el volumen de una esfera de radio r0 en cm

m 0 = 4 3 π r 0 3

El valor de la constante de proporcionalidad k es

k= ρ n π ( ρ a 4 3 π ) 2/3 = 10 6 π ( 4 3 π ) 2/3 =1.21· 10 6

Donde ρn≈10-6 g/cm3 es la densidad de la niebla, y ρa=1.0 g/cm3 es la densidad del agua. La constante de proporcionalidad k es por tanto, del orden de 10-6.

La aceleración de la gravedad es g=980 cm/s2.

Actividades

Se introduce

Se pulsa el botón titulado Empieza

Se observa la caída de la gota de agua en la parte izquierda del applet

Se representa la aceleración dv/dt de la gota en función del tiempo t, observando que tiende hacia el valor límite g/7=1.4 m/s2.

Vemos como la gota cambia su tamaño a medida que absorbe las pequeñas gotas suspendidas en el aire y que forman la niebla.

Referencias

Adawi I. Comments on the raindrop problem. Am. J. Phys. 54 (8) August 1986, pp. 739-740

Anterior