title: Función de flujo de Stokes
text: # 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.$
  