Movimiento de un objeto bajo la acción de una fuerza de rozamiento proporcional a la velocidad y al cuadrado de la velocidad.

En las páginas anteriores, hemos estudiando el movimiento de caída de un cuerpo bajo la acción de su peso y de una fuerza de rozamiento en los dos casos en los que es posible obtener expresiones analíticas para la evolución de la velocidad en función del tiempo y de la posición (altura) del móvil en función del tiempo.

La fórmula general de la fuerza de rozamiento de un cuerpo de densidad ρ que se mueve en el seno de un fluido de densidad ρf  y viscosidad η es

F r = 1 2 C d ρ f A v 2

Donde Cd se denomina coeficiente de arrastre, A es el área de la sección transversal al movimiento (en el caso de una esfera es πR2) y v es la velocidad.

El coeficiente de arrastre es una función del número de Reynolds, Re. Este número que es importante para definir la transición del flujo laminar al turbulento, se define como

Re= ρ f lv η

donde l representa la longitud del objeto medida a lo largo de su sección transversal (en el caso de una esfera es 2R).

Para un amplio intervalo de números Re, la dependencia la forma funcional del coeficiente de arrastre Cd se puede escribir.

C d 24 Re + 6 1+ Re +0.4

En esta página, vamos a considerar el caso de un cuerpo que se mueve a lo largo del eje X, bajo la acción de una fuerza de rozamiento que es una combinación de una fuerza proporcional a la velocidad y otra proporcional al cuadrado de la velocidad.

 

Fuerza de rozamiento proporcional a la velocidad

La ecuación del movimiento del cuerpo a lo largo del eje X es

m dv dt =6πηRv dv dt = 6πηR ρ 4 3 π R 3 v= 9η 2ρ R 2 v dv dt = 1 τ vτ= 2ρ R 2 9η

ρ es la densidad del cuerpo y η es la viscosidad del fluido

La solución de esta ecuación diferencial con la condición inicial t=0, v=v0 es

v=v0·exp(-t/τ)

Integrando con respecto del tiempo obtenemos la posición x del móvil en función del tiempo t

dx dt = v 0 exp( t/τ ) x= 0 t v 0 exp( t/τ )dt= v 0 τ( 1exp( t/τ ) )

Cuando t→∞, la posición final del objeto tiende a xv0τ

Fuerza de rozamiento proporcional al cuadrado de la velocidad

La ecuación del movimiento del cuerpo a lo largo del eje X es

m dv dt =0.2 ρ f π R 2 v 2 dv dt = 0.2 ρ f π R 2 ρ 4 3 π R 3 v 2 = 3 ρ f 20ρR v 2 dv dt = 1 τ v c v 2 v c = 30η ρ f R

ρf es la densidad y η es la viscosidad del fluido

La solución de esta ecuación diferencial con la condición inicial t=0, v=v0 es

1 v = 1 v 0 + t v c τ v= v 0 ( v c τ) ( v c τ)+ v 0 t

Cuando v0>>vc, en un tiempo corto t=τ la fuerza de rozamiento proporcional al cuadrado de la velocidad hace que ésta se reduzca a aproximadamente vc.

Integrando con respecto del tiempo obtenemos la posición x del móvil en función del tiempo t

x= 0 t v 0 v c τ v c τ+ v 0 t dt= v c τln( 1+ v 0 t v c τ )

La fuerza de rozamiento proporcional a la velocidad y al cuadrado de la velocidad

La ecuación del movimiento del cuerpo a lo largo del eje X es

dv dt = 1 τ v 1 τ v c v 2

Haciendo un cambio de variable, obtenemos una ecuación diferencial fácil de integrar

u= 1 v + 1 v c du dt = u τ

La solución de esta ecuación diferencial con la condición inicial t=0, v=v0 es

u= u 0 exp(t/τ) 1 v + 1 v c =( 1 v 0 + 1 v c )exp(t/τ) v= v 0 v c ( v 0 + v c )exp(t/τ) v 0

Cuando t→∞, la velocidad final del objeto tiende a cero

Integrando con respecto del tiempo, obtenemos la posición x del móvil en función del tiempo t

x= 0 t v 0 v c ( v 0 + v c )exp(t/τ) v 0 dt

Hacemos el cambio de variable

u=( v 0 + v c )exp(t/τ) v 0 du= u+ v 0 τ dt v 0 v c τ du u(u+ v 0 ) = v c τ ( 1 u 1 u+ v 0 )du= v c τln u u+ v 0

Deshaciendo el cambio de variable y evaluando la integral entre los limites 0 y t, obtenemos

x= v c τln( v 0 v c ( 1exp(t/τ) )+1 )

Cuando t→∞, la posición final del objeto tiende a

x v c τln( v 0 v c +1 )

Ejemplo

En el instante t=3 s

v= v 0 v c ( v 0 + v c )exp(t/τ) v 0 v= 10·8 (10+8)exp(3/9)10 =5.29m/s x= v c τln( v 0 v c ( 1exp(t/τ) )+1 ) x=8·9ln( 10 8 ( 1exp(3/9) )+1 )=21.83m

Cuando t→∞, la posición final del objeto tiende a

v0 x v c τln( v 0 v c +1 ) x8·9ln( 10 8 +1 )=58.39m

Actividades

Se selecciona el botón de radio:

Se introduce

Se pulsa el botón titulado Empieza

Se observa el movimiento del objeto a lo largo del eje X, hasta que alcanza la posición final.

Referencias

Saslow M. W., Hong Lu. Newton on objects moving in a fluid-the penetration length. Eur. J. Phys. 29 (2008) pp. 689-696.