Ver el código de «Teoría de flujo potencial»


La '''teoría de flujo potencial''' pretende describir el comportamiento [[cinemática|cinemático]] de los [[fluido]]s basándose en el concepto matemático de [[Potencial (física)|función potencial]], asegurando que el campo de velocidades (que es un [[campo vectorial]]) del flujo de un fluido es igual al [[gradiente]] de una función potencial que determina el movimiento de dicho fluido:

<math> \vec{V} = - \nabla\phi</math> 

donde el campo de velocidades queda definido como

<math>\vec{V}= u \hat{i} + v \hat{j} + w \hat{k}</math> 

El signo menos en la ecuación de arriba es sólo una convención de signos sobre la definición de <math>\phi</math>. <math>\phi</math> Puede definirse sin el signo menos, y la formulación que se obtendría sería la misma. A un fluido que se comporta según esta teoría se le denomina ''fluido potencial'', que da lugar a un ''flujo potencial''. 

Una de las primeras personas en aplicar esta formulación para el flujo de un fluido fue [[D'Alembert]]. Estudió la fuerza de resistencia producida por un flujo de fluido sobre un cuerpo que se oponía a éste en dos dimensiones cuando la solución a este problema era completamente desconocida y [[Isaac Newton|Newton]], a pesar de haberlo estudiado, no había llegado a conclusiones satisfactorias.

D'Alembert definió la función de corriente, <math>\psi</math>, para describir la trayectoria que tuviera cada [[partícula de fluido|partícula de un fluido]] a través del tiempo. Esta función corriente está determinada, en el plano, por dos variables espaciales y para cada valor de <math>\psi</math> la igualdad <math>\psi=\psi(x,y)</math> determina un [[lugar geométrico]] llamado línea de corriente.

== La Naturaleza de <math>\psi</math> y su relación con <math>\phi</math> ==

Primeramente definiremos la función corriente en el plano, para luego explicar sus características. La función <math>\psi</math> se define como aquella que cumple con las siguientes condiciones:


<math>
\frac{\partial \psi}{\partial y}=u 
</math>
y 
<math>
\frac{\partial \psi}{\partial x}=-v
</math>

Las líneas de corriente determinan la trayectoria de una partícula de fluido que se encuentra sobre estas. Así, por ejemplo, si una partícula de fluido se halla sobre la línea ''equipotencial'' de <math>\psi=3</math>,  tendrá una trayectoria que se situará exactamente sobre el lugar geométrico que determinará la igualdad <math>\psi(x,y)=3</math> (línea de corriente=trayectoria es debido a que contemplamos un movimiento plano independiente de t "ψ(x,y)"). Esta propiedad de las líneas de corriente exige que las funciones <math>\psi</math> y <math>\phi</math> estén "sincronizadas", ya que la velocidad en cualquier punto del flujo de fluido será siempre tangente a la trayectoria de la línea de corriente. Fácilmente se puede demostrar que la familia de curvas determinadas por la función corriente y la función potencial de velocidades forman una red ortogonal, como se verá a continuación:

Partimos del diferencial total de la función <math>\phi</math>:

<math>d\phi=\frac{\partial \phi}{\partial x} dx+\frac{\partial \phi}{\partial y} dy</math>

Así en cualquier curva equipotencial <math>\phi=\mbox{constante}</math> se cumplirá que 

<math>\frac{\partial \phi}{\partial x} dx+\frac{\partial \phi}{\partial y} dy = 0</math>

Esto implica que: 

<math> (\frac{dy}{dx})_{\phi=cte}=\frac{-u}{v}</math>

La misma propiedad se aplica a cada línea de corriente:

<math>\frac{\partial \psi}{\partial x} dx+\frac{\partial \psi}{\partial y} dy = 0</math>

y 


<math> (\frac{dy}{dx})_{\psi=cte}=\frac{v}{u}</math>

por lo cual de determina que:


<math> (\frac{dx}{dy})_{\psi=cte}=-(\frac{dy}{dx})_{\phi=cte} </math>

Esta propiedad de ambas funciones permite intercambiarlas para generar otros patrones de flujo y, como las líneas de corriente no pueden cortarse entre sí, no existe ningún caudal que las atraviesa perpendicular a estas. Esto permite suponer a las líneas de corriente límites materiales, es decir, paredes u obstáculos que restringen o determinan el flujo que se desea estudiar. Esto ya lo supuso D'Alembert al estudiar el efecto de empuje de un flujo corriendo sobre un objeto que lo obstaculiza. Estudiando las propiedades de la función corriente y la función potencial determinó que podían superponerse para generar así un patrón de fluido que combinara diversos movimientos. Superponiendo una fuente y un sumidero de igual caudal obtuvo una circunferencia, que combinó con un flujo uniforme para modelar el flujo de fluido sobre un cilindro de longitud infinita. Una vez obtenido esto, demostró que la suma de las presiones sobre el cilindro se anulaban, lo que hacía que la fuerza resultante sobre el cilindro fuera cero. Esta es la llamada [[paradoja de D'Alembert]].

== Algunos patrones de flujo simple ==

=== Flujo uniforme ===

Un flujo el cual sea uniforme en una misma dirección cumple con que <math>\vec{V}=U_{0}\hat{n}</math>
donde <math>\hat{n}</math> es la dirección del flujo. Si tomamos esta dirección como la del eje <math>x</math>, obtendremos que el campo de velocidades estará dado por:

<math>\vec{V}=U_{0}\hat{i}</math>

Con lo cual podremos encontrar la función potencial integrando:

<math>
\phi=-\int U_{0} dx +f(y)
</math>

Sabiendo que <math>\frac{\partial \phi}{\partial y}=0</math> y que no existe componente vertical de la velocidad en ningún punto del flujo tenemos:

<math>
\phi=-\int U_{0} dx + C_1
</math>
 
En esta ecuación podemos tomar la constante de integración igual a cero. Quedando las funciones potencial y corriente como:

<math>
\mathrm{\phi} = -U_{0} x
</math>

<math>
\psi = -U_{0} y
</math>

=== Fuentes y sumideros ===
Una fuente o un sumidero de algún fluido tiene la particularidad de que el flujo sólo sale o entra, lo que implica que el vector velocidad para cada punto del flujo será colineal al origen para ambos casos. Es mucho más sencillo hallar esta función potencial usando coordenadas polares. Así:



<math>
v_r = \frac{Q}{2 \pi r}
</math>



<math>
v_\theta=0
</math>


Donde <math>Q</math> es el caudal que sale si es positivo o entra si es negativo. Para hallar la función potencial integramos:

<math>
-\phi=\int \frac{Q}{2 \pi r} dr +f(\theta)
</math>

Como la velocidad en <math>\theta</math> es igual a cero sólo queda una constante de integración la cual podemos hacer cero; entonces:

<math>\phi = -\frac{Q ln(r)}{2 \pi}</math>

Para obtener la función corriente podemos realizar un procedimiento análogo considerando la forma del operador [[gradiente]] en [[coordenadas polares]]:

<math>
-\frac{1}{r}\frac{\partial \psi}{\partial \theta} = \frac{Q}{2 \pi r}
</math>

<math>
\int \frac{Q}{2 \pi} d\theta = - \psi + f(r)
</math>

entonces:

<math>
\psi=\frac{Q}{2 \pi} \theta
</math>

== Véase también ==
* [[Divergencia]]
* [[Rotacional]]

{{Control de autoridades}}
[[Categoría:Mecánica de fluidos]]
[[Categoría:Ingeniería mecánica]]
[[Categoría:Ingeniería civil]]