Función de flujo de Stokes
# Función de flujo de Stokes

En dinámica de fluidos, la función de corriente de Stokes se utiliza para describir las líneas de corriente y la velocidad de flujo en un flujo incompresible tridimensional con axisimetría. Una superficie con un valor constante de la función de corriente de Stokes encierra un tubo de corriente, en todas partes tangente a los vectores de velocidad del flujo. Además, el volumen dentro de este tubo de corriente es constante, y todas las líneas de corriente del flujo están situadas en esta superficie. El campo de velocidad asociado a la función de flujo de Stokes es un solenoidal, tiene divergencia cero. Esta función de flujo recibe su nombre en honor de George Gabriel Stokes.

## Coordenadas cilíndricas

Consideremos un sistema de coordenadas cilíndricas ( ρ , φ , z ), con el eje z la línea alrededor de la cual el flujo incompresible es axisimétrico, φ el acimut y ρ la distancia al eje z. Entonces las componentes de la velocidad del flujo $u_ρ$ y $u_z$ pueden expresarse en términos de la función de corriente de Stokes $\Psi$ mediante:

    ${\begin{aligned}u_{\rho }&=-{\frac {1}{\rho }}\,{\frac {\partial \Psi }{\partial z}},\\u_{z}&=+{\frac {1}{\rho }}\,{\frac {\partial \Psi }{\partial \rho }}.\end{aligned}}$

La componente azimutal de la velocidad $u_φ$ no depende de la función de flujo. Debido a la axisimetría, las tres componentes de la velocidad ( $u_ρ$ , $u_φ$ , $u_z$ ) sólo dependen de ρ y z y no del acimut φ.

El flujo de volumen, a través de la superficie limitada por un valor constante ψ de la función de corriente de Stokes, es igual a 2π ψ.

## Coordenadas esféricas

En coordenadas esféricas ( r , θ , φ ), r es la distancia radial desde el origen, θ es el ángulo cenital y φ es el ángulo azimutal. En flujo axisimétrico, con θ= 0 el eje de simetría rotacional, las cantidades que describen el flujo son de nuevo independientes del azimut φ. Las componentes de la velocidad del flujo $u_r$ y $u_θ$ están relacionadas con la función de corriente de Stokes $\Psi$ a través de:

    ${\begin{aligned}u_{r}&=+{\frac {1}{r^{2}\,\sin \theta }}\,{\frac {\partial \Psi }{\partial \theta }},\\u_{\theta }&=-{\frac {1}{r\,\sin \theta }}\,{\frac {\partial \Psi }{\partial r}}.\end{aligned}}$

De nuevo, la componente azimutal de la velocidad $u_φ$ no es una función de la función de Stokes ψ. El flujo de volumen a través de un tubo de corriente, limitado por una superficie de ψ constante, es igual a 2π ψ, como antes.

### Vorticidad

La vorticidad se define como:

    ${\boldsymbol {\omega }}=\nabla \times {\boldsymbol {u}}=\nabla \times \nabla \times {\boldsymbol {\psi }}$, donde ${\boldsymbol {\psi }}=-{\frac {\Psi }{r\sin \theta }}{\boldsymbol {\hat {\phi }}},$

con ${\boldsymbol {\hat {\phi }}}$ el vector unitario en la dirección.$\phi \,$

    Derivación de la vorticidad ${\boldsymbol {\omega }}$ utilizando una función de flujo de Stokes Consideremos la vorticidad definida por ${\boldsymbol {\omega }}=\nabla \times {\boldsymbol {u}}.$ A partir de la definición del gradiente en coordenadas esféricas: ${\begin{aligned}\omega _{r}&={1 \over r\sin \theta }\left({\partial \over \partial \theta }\left(u_{\phi }\sin \theta \right)-{\partial u_{\theta } \over \partial \phi }\right){\boldsymbol {\hat {r}}},\\\omega _{\theta }&={1 \over r}\left({1 \over \sin \theta }{\partial u_{r} \over \partial \phi }-{\partial \over \partial r}\left(ru_{\phi }\right)\right){\boldsymbol {\hat {\theta }}},\\\omega _{\phi }&={1 \over r}\left({\partial \over \partial r}\left(ru_{\theta }\right)-{\partial u_{r} \over \partial \theta }\right){\boldsymbol {\hat {\phi }}}.\end{aligned}}$ Primero se observa que las componentes $r$ y $\theta$ son iguales a 0. Segundo sustituye $u_{r}$ y $u_{\theta }$ en $\omega _{\phi }.$ El resultado es: ${\begin{aligned}\omega _{r}&=0,\\\omega _{\theta }&=0,\\\omega _{\phi }&={1 \over r}\left({\partial \over \partial r}\left(r\left(-{\frac {1}{r\sin \theta }}{\frac {\partial \Psi }{\partial r}}\right)\right)-{\partial \over \partial \theta }\left({\frac {1}{r^{2}\sin \theta }}{\frac {\partial \Psi }{\partial \theta }}\right)\right).\end{aligned}}$ A continuación se realiza la operación algebraica siguiente: ${\begin{aligned}\omega _{\phi }&={1 \over r}\left(-{\frac {1}{\sin \theta }}\left({\partial \over \partial r}\left({\frac {\partial \Psi }{\partial r}}\right)\right)-{\frac {1}{r^{2}}}{\partial \over \partial \theta }\left({\frac {1}{\sin \theta }}{\frac {\partial \Psi }{\partial \theta }}\right)\right)\\&={1 \over r}\left(-{\frac {1}{\sin \theta }}\left({\frac {\partial ^{2}\Psi }{\partial r^{2}}}\right)-{\frac {\sin \theta }{r^{2}\sin \theta }}{\partial \over \partial \theta }\left({\frac {1}{\sin \theta }}{\frac {\partial \Psi }{\partial \theta }}\right)\right)\\&=-{\frac {1}{r\sin \theta }}\left({\frac {\partial ^{2}\Psi }{\partial r^{2}}}+{\frac {\sin \theta }{r^{2}}}{\partial \over \partial \theta }\left({\frac {1}{\sin \theta }}{\frac {\partial \Psi }{\partial \theta }}\right)\right).\end{aligned}}$

Como resultado, del cálculo se encuentra que el vector vorticidad es igual a:

    ${\boldsymbol {\omega }}={\begin{pmatrix}0\\[1ex]0\\[1ex]\displaystyle -{\frac {1}{r\sin \theta }}\left({\frac {\partial ^{2}\Psi }{\partial r^{2}}}+{\frac {\sin \theta }{r^{2}}}{\partial \over \partial \theta }\left({\frac {1}{\sin \theta }}{\frac {\partial \Psi }{\partial \theta }}\right)\right)\end{pmatrix}}.$

### Comparación con las cilíndricas

Los sistemas de coordenadas cilíndricas y esféricas se relacionan a través de

    $z=r\,\cos \theta \,$ y $\rho =r\,\operatorname {sen} \theta .\,$

## Definición alternativa con signo opuesto

Como se explica en el artículo general función de corriente, las definiciones que utilizan una convención de signo opuesto -para la relación entre la función de corriente de Stokes y la velocidad del flujo- también están en uso.

## Divergencia cero

En coordenadas cilíndricas, la divergencia del campo de velocidades u se convierte en:

    ${\begin{aligned}\nabla \cdot {\boldsymbol {u}}&={\frac {1}{\rho }}{\frac {\partial }{\partial \rho }}{\Bigl\rho \,u_{\rho }{\Bigr )}+{\frac {\partial u_{z}}{\partial z}}\\&={\frac {1}{\rho }}{\frac {\partial }{\partial \rho }}\left(-{\frac {\partial \Psi }{\partial z}}\right)+{\frac {\partial }{\partial z}}\left({\frac {1}{\rho }}{\frac {\partial \Psi }{\partial \rho }}\right)=0,\end{aligned}}$

como era de esperar para un flujo incompresible.

Y en coordenadas esféricas:

    ${\begin{aligned}\nabla \cdot {\boldsymbol {u}}&={\frac {1}{r\,\sin \theta }}{\frac {\partial }{\partial \theta }}(u_{\theta }\,\sin \theta )+{\frac {1}{r^{2}}}{\frac {\partial }{\partial r}}{\Biglr^{2}\,u_{r}{\Bigr )}\\&={\frac {1}{r\,\sin \theta }}{\frac {\partial }{\partial \theta }}\left(-{\frac {1}{r}}{\frac {\partial \Psi }{\partial r}}\right)+{\frac {1}{r^{2}}}{\frac {\partial }{\partial r}}\left({\frac {1}{\sin \theta }}{\frac {\partial \Psi }{\partial \theta }}\right)=0.\end{aligned}}$

## Líneas de corriente como curvas de función de corriente constante

Del cálculo se sabe que el vector gradiente $\nabla \Psi$ es normal a la curva $\Psi =C$. Si se demuestra que en todas partes ${\boldsymbol {u}}\cdot \nabla \Psi =0,$ utilizando la fórmula para ${\boldsymbol {u}}$ en términos de $\Psi ,$ entonces se demuestra que las curvas de nivel de $\Psi$ son líneas de corriente.

### Coordenadas cilíndricas

En coordenadas cilíndricas,

    $\nabla \Psi ={\partial \Psi \over \partial \rho }{\boldsymbol {e}}_{\rho }+{\partial \Psi \over \partial z}{\boldsymbol {e}}_{z}$.

y

    ${\boldsymbol {u}}=u_{\rho }{\boldsymbol {e}}_{\rho }+u_{z}{\boldsymbol {e}}_{z}=-{1 \over \rho }{\partial \Psi \over \partial z}{\boldsymbol {e}}_{\rho }+{1 \over \rho }{\partial \Psi \over \partial \rho }{\boldsymbol {e}}_{z}.$

De modo que

    $\nabla \Psi \cdot {\boldsymbol {u}}={\partial \Psi \over \partial \rho }(-{1 \over \rho }{\partial \Psi \over \partial z})+{\partial \Psi \over \partial z}{1 \over \rho }{\partial \Psi \over \partial \rho }=0.$

### Coordenadas esféricas

Y en coordenadas esféricas

    $\nabla \Psi ={\partial \Psi \over \partial r}{\boldsymbol {e}}_{r}+{1 \over r}{\partial \Psi \over \partial \theta }{\boldsymbol {e}}_{\theta }$

y

    ${\boldsymbol {u}}=u_{r}{\boldsymbol {e}}_{r}+u_{\theta }{\boldsymbol {e}}_{\theta }={1 \over r^{2}\sin \theta }{\partial \Psi \over \partial \theta }{\boldsymbol {e}}_{r}-{1 \over r\sin \theta }{\partial \Psi \over \partial r}{\boldsymbol {e}}_{\theta }.$

De modo que

    $\nabla \Psi \cdot {\boldsymbol {u}}={\partial \Psi \over \partial r}\cdot {1 \over r^{2}\sin \theta }{\partial \Psi \over \partial \theta }+{1 \over r}{\partial \Psi \over \partial \theta }\cdot {\Big-{1 \over r\sin \theta }{\partial \Psi \over \partial r}{\Big )}=0.$
