Ley de Stokes

La ley de Stokes se refiere a la fuerza de fricción experimentada por objetos esféricos moviéndose en el seno de un fluido viscoso en un régimen laminar de bajos números de Reynolds. En general la ley de Stokes es válida en el movimiento de partículas esféricas pequeñas moviéndose a velocidades bajas.

Fue derivada en 1851 por George Gabriel Stokes tras resolver un caso particular de las ecuaciones de Navier-Stokes.

La ley de Stokes es el principio usado en los viscosímetros de esfera en caída libre, en los cuales el fluido está estacionario en un tubo vertical de vidrio y una esfera, de tamaño y densidad conocidas, desciende a través del líquido. Si la bola ha sido seleccionada correctamente alcanzará la velocidad terminal, la cual puede ser medida por el tiempo que pasa entre dos marcas de un tubo.

A veces se usan sensores electrónicos para fluidos opacos. Conociendo las densidades de la esfera, el líquido y la velocidad de caída se puede calcular la viscosidad a partir de la fórmula de la ley de Stokes. Para mejorar la precisión del experimento se utilizan varias bolas. La técnica es usada en la industria para verificar la viscosidad de los productos, por ejemplo, la glicerina o el sirope.

La importancia de la ley de Stokes está ilustrada en el hecho de que ha jugado un papel crítico en la investigación de al menos 3 Premios Nobel.^[1]​

La ley de Stokes también es importante para la compresión del movimiento de microorganismos en un fluido, así como los procesos de sedimentación debido a la gravedad de pequeñas partículas y organismos en medios acuáticos.^[2]​ También es usado para determinar el porcentaje de granulometría muy fina de un suelo mediante el ensayo de sedimentación.

En la atmósfera, la misma teoría puede ser usada para explicar porque las gotas de agua (o los cristales de hielo) pueden permanecer suspendidos en el aire (como nubes) hasta que consiguen un tamaño crítico para empezar a caer como lluvia (o granizo o nieve). Usos similares de la ecuación pueden ser usados para estudiar el principio de asentamiento de partículas finas en agua u otros fluidos.

(velocidad límite)

La ley de Stokes puede escribirse como:

${displaystyle F_{d}=6pi R mu u}$

La condición de bajos números de Reynolds implica un flujo laminar lo cual puede traducirse por una velocidad relativa entre la esfera y el medio inferior a un cierto valor crítico. En estas condiciones la resistencia que ofrece el medio es debida casi exclusivamente a las fuerzas de rozamiento que se oponen al deslizamiento de unas capas de fluido sobre otras a partir de la capa límite adherida al cuerpo. La ley de Stokes se ha comprobado experimentalmente en multitud de fluidos y condiciones.

Si las partículas están cayendo verticalmente en un fluido viscoso debido a su propio peso puede calcularse su velocidad de caída o sedimentación igualando la fuerza de fricción con el peso aparente de la partícula en el fluido.

${displaystyle u_{s}={Bigl (}{frac {2}{9}}{Bigr )}{frac {( ho _{p}- ho _{f}) g r^{2}}{mu }}}$

En flujos de Stokes con un número de Reynolds muy bajo, la aceleración convectiva se puede considerar nula en los términos de la ecuación de Navier-Stokes. En ese caso las ecuaciones del flujo se igualan a las de un flujo incompresible y estacionario:^[3]​

${displaystyle {egin{aligned}& abla p=mu , abla ^{2}mathbf {u} =-mu , abla imes mathbf {oldsymbol {omega }} ,\& abla cdot mathbf {u} =0,end{aligned}}}$

${displaystyle {oldsymbol {omega }}= abla imes mathbf {u} }$

Usando algunas propiedades del cálculo de vectores, estas ecuaciones se pueden mostrar como resultado de una ecuación de Laplace para la presión y cada uno de los componentes del vector vorticidad:^[3]​

${displaystyle abla ^{2}{oldsymbol {omega }}=0}$

${displaystyle abla ^{2}p=0}$

Fuerzas adicionales como la gravedad o la flotabilidad no han sido tomados en cuenta, pero pueden ser fácilmente añadidos a la ecuación ya que son lineales, así que se puede aplicar la superposición lineal a las soluciones.

Para el caso de una esfera en un campo de velocidades, es ventajoso usar el sistema de coordenadas cilíndrico (${displaystyle r,varphi ,z}$). El eje (${displaystyle z}$) pasa por el centro de la esfera y está alineado con la dirección del flujo, mientras que (${displaystyle r}$) es el radio medido perpendicular al eje (${displaystyle z}$). El origen es el centro es de la esfera. Debido a que el flujo es asimétrico respecto al eje z, éste es independiente del azimut (${displaystyle varphi }$).

En el sistema de coordenadas cilíndrico, el flujo incompresible puede ser descrito por la función del flujo de Stokes (${displaystyle psi }$), la cual está en función de (${displaystyle r}$) y (${displaystyle z}$):^[4]​^[5]​

${displaystyle v=-{frac {1}{r}}{frac {partial psi }{partial z}},qquad w={frac {1}{r}}{frac {partial psi }{partial r}}}$

con (${displaystyle v}$) y (${displaystyle w}$) como componentes del flujo de velocidad en la dirección (${displaystyle r}$) y (${displaystyle z}$), respectivamente. La componente de la velocidad acimutal en la dirección φ es cero, en el caso simétrico. El flujo de volumen, a través de un tubo limitado por una superficie de valor constante (${displaystyle psi }$), es igual a (${displaystyle 2pi psi }$) y es constante.^[4]​

Para el caso de un flujo simétrico por los ejes, el único componente no nulo del vector vorticidad (${displaystyle omega }$) es el azimutal (${displaystyle varphi }$), el componente (${displaystyle omega _{varphi }}$).^[6]​^[7]​

${displaystyle omega _{varphi }={frac {partial v}{partial z}}-{frac {partial w}{partial r}}=-{frac {partial }{partial r}}left({frac {1}{r}}{frac {partial psi }{partial r}} ight)-{frac {1}{r}} {frac {partial ^{2}psi }{partial z^{2}}}}$

El operador de Laplace, aplicado a la vorticidad (${displaystyle omega _{varphi }}$), aplicado en el sistema cilíndrico con simetría en los ejes:^[7]​

${displaystyle abla ^{2}omega _{varphi }={frac {1}{r}}{frac {partial }{partial r}}left(r {frac {partial omega _{varphi }}{partial r}} ight)+{frac {partial ^{2}omega _{varphi }}{partial z^{2}}}-{frac {omega _{varphi }}{r^{2}}}=0.}$

De las dos ecuaciones anteriores, y con las apropiadas condiciones de contornos, para un campo de velocidad uniforme y paralela (${displaystyle V}$) en la dirección (${displaystyle z}$) y en una esfera de radio (${displaystyle R}$), la solución resulta ser^[8]​

${displaystyle psi =-{frac {1}{2}} V r^{2} left[1-{frac {3}{2}}{frac {R}{sqrt {r^{2}+z^{2}}}}+{frac {1}{2}}left({frac {R}{sqrt {r^{2}+z^{2}}}} ight)^{3} ight].}$

La fuerza viscosa por unidad de área σ, ejercida por el flujo en la superficie de la esfera, está en la dirección z sobre toda la esfera. Más exactamente, tiene el mismo valor en cualquier punto de la esfera:

${displaystyle {oldsymbol {sigma }}={frac {3 mu V}{2 R}} mathbf {e} _{z}}$

con (${displaystyle mathbf {e} _{z}}$) el vector unitario en la dirección (${displaystyle z}$). Para otras formas que no sean la esférica, (${displaystyle sigma }$) no es constante a lo largo de la superficie del cuerpo. Integrando la fuerza viscosa por unidad de área (${displaystyle sigma }$) sobre la esfera resulta la fuerza de fricción (${displaystyle F_{d}}$) de acuerdo con la ley de Stokes.

---