2011

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 1

FLUJO POTENCIAL PLANO

FLUJO POTENCIAL
PLANO
POTENCIAL COMPLEJO DE VELOCIDADES
La teora de flujo potencial pretende describir el comportamiento
cinemtico de los fluidos basndose en el concepto matemtico de
funcin potencial, asegurando que el campo de velocidades (que es un
campo vectorial) del flujo de un fluido es igual al gradiente de una
funcin potencial que determina el movimiento de dicho fluido.

HUAYANEY MILLA MIGUEL

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA
16/12/2011

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 2

FLUJO POTENCIAL PLANO

UNIVERSIDAD NACIONAL DE ANCASH

SANTIAGO ANTUNEZ DE MAYOLO
FACULTAD DE INGENIERIA DE MINAS
GEOLOGIA Y METALURGIA
Escuela Profesional de Ingeniera De Minas

TRABAJO

TEMA: Flujo Potencial Plano

ALUMNOS:
HUAYANEY MILLA MIGUEL
DOCENTE:
Ing. BAYONA ANTUNEZ REMO CRISANTO

2011

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 3

FLUJO POTENCIAL PLANO

CAPITULO I

FLUJO POTENCIAL PLANO

1.1.

ANTECEDENTES
[]Los principios bsicos del movimiento de los fluidos se desarrollaron
lentamente a travs de los siglos XVI al XIX como resultado del trabajo de
muchos cientficos como Da Vinci, Galileo, Torricelli, Pascal, Bernoulli,
Euler, Navier, Stokes, Kelvin, Reynolds y otros que hicieron interesantes
aportes tericos a lo que se denomina hidrodinmica. Tambin en el
campo de hidrulica experimental hicieron importantes contribuciones
Chezy, Ventura, Hagen, Manning, Poiseville, Darcy, Froude y otros,
fundamentalmente durante el siglo XIX [] 1
Gracias a estos cientficos podemos comprender el resultado de sus
trabajos en la aplicacin a la mecnica de fluidos, que explican con
modelos matemticos, el comportamiento de los fluidos, para lo cual se
recurri a simplificar las propiedades de estos. As se hablaba de un fluido

http://www.monografias.com/trabajos12/mecflui/mecflui.shtml

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 4

FLUJO POTENCIAL PLANO

real. Esto hizo que los resultados no fueran siempre aplicables a casos
reales sino tambin a fluidos ideales, lo cual desarrollaron modelos
matemticos complejos que detallan el comportamiento de un flujo ideal
por medio del flujo potencial o potencial de velocidades que se desarrollo
estos modelos y se simplificaron algunos en la actualidad.
1.2.

DEFINICIN 2
La teora de flujo potencial pretende describir el comportamiento
cinemtico de los fluidos basndose en el concepto matemtico de
funcin potencial, asegurando que el campo de velocidades (que es un
campo vectorial) del flujo de un fluido es igual al gradiente de una funcin
potencial que determina el movimiento de dicho fluido:

Con

El signo menos en la ecuacin de arriba es slo una convencin de signos

sobre la definicin de . puede definirse sin el signo menos y la
formulacin que se obtendra sera la misma.
A un fluido que se comporta segn esta teora se le denomina fluido
potencial y a un flujo potencial.
Una de las primeras personas en aplicar esta formulacin para el flujo de
un fluido fue D'Alembert. l estudi la fuerza de resistencia producida por
un flujo de fluido sobre un cuerpo que se opona a ste en dos
dimensiones cuando este problema era completamente obscuro y

Definicin hecha por Wiki peda Enciclopedia Libre.

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 5

FLUJO POTENCIAL PLANO

Newton, a pesar de haberlo estudiado, no haba llegado a conclusiones

satisfactorias.
D'Alembert defini la funcin de corriente, , para describir la trayectoria
que tuviera cada partcula de un fluido a travs del tiempo. Esta funcin
corriente est determinada, en el plano, por dos variables espaciales y
para cada valor de la igualdad = (x, y) determina un lugar
geomtrico llamado lnea de corriente.
1.3.

EL POTENCIAL DE VELOCIDAD
El movimiento de los fluidos perfectos (no viscosos) se describe mediante
la ecuacin de Euler:

(1)

Este tipo de flujos es muy importante pues en muchas situaciones de

inters prctico, los efectos de la viscosidad de los fluidos reales quedan
limitados a las regiones del espacio (muchas veces pequeas) donde
tienen lugar fuertes gradientes de la velocidad (capas lmite o regiones
donde el flujo es turbulento), mientras que en el grueso del flujo los
efectos de la viscosidad son despreciables y el fluido se puede suponer
ideal. En las regiones materiales de flujo ideal no se crea ni se destruye
vorticosidad, de manera que si en un dado instante sta es nula, sigue
siendo nula en todo otro momento.
Una de las propiedades fundamentales de los fluidos invscidos es que
puede haber flujos que son permanentemente irrotacionales, es decir
que u = 0 en todos los puntos del fluido. Si adems el flujo es
incompresible, el campo de velocidad satisface las condiciones:
u = 0

u=0

(2)

La irrotacionalidad del campo de velocidad implica que u deriva de un

potencial de velocidad , esto es u = (3) y la incompresibilidad implica
entonces:

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 6

FLUJO POTENCIAL PLANO

2 2 2
= 2 + 2+ 2 =0
x
y
z
2

(4)

Por lo tanto satisface la ecuacin de Laplace, que junto con la (3)

determina u(r, t). Viceversa, la existencia de una funcin escalar que
satisface la ecuacin (4) implica la existencia de una funcin vectorial
u = que cumple las condiciones (2), como se puede verificar
fcilmente. Los flujos para los cuales se cumple la (3) se denominan flujos
potenciales. Claramente, la (3) implica que las lneas de corriente de un
flujo potencial son ortogonales a las superficies equipotenciales, definidas
por = cte.
En los flujos potenciales, la ecuacin de movimiento adopta esta forma:
2

+ ++ =

(5)

Donde indica el potencial de las fuerzas de volumen. En esta

aproximacin , y entonces u, quedan determinados por una ecuacin
lineal (la ecuacin de Laplace). Sin embargo, de acuerdo con la (5) p
depende de u en forma no lineal.
La solucin de un problema de flujo potencial 3 consiste pues en la
determinacin de dos magnitudes, y p, para lo cual disponemos de las
dos ecuaciones escalares (4) y (5). An as, el problema es sencillo slo
para flujos estacionarios. Ntese, sin embargo, que si las condiciones de
contorno no involucran la presin, la ecuacin de Laplace para se
puede resolver sin necesidad de conocer p; en tal caso, la (5) slo sirve
para calcular p una vez determinado .
Debido a que la ecuacin de Laplace es lineal, las combinaciones lineales
de soluciones son tambin soluciones de la misma. Se puede entonces
3

Flujo potencial tambin llamado o conocido como potencial de velocidades, y se utiliza mas en
fluidos ideales, por lo que presenta modelos matemticos complejos.

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 7

FLUJO POTENCIAL PLANO

construir el campo de velocidad de un problema de flujo potencial

superponiendo soluciones simples ya conocidas, de forma tal de
satisfacer las condiciones de contorno correspondientes a la solucin
buscada. Por otra parte, del mismo modo que en la Electrosttica cuando
se calcula el potencial elctrico de una distribucin de cargas, podremos
escribir el potencial de velocidad de un flujo potencial en forma de un
desarrollo multipolar. Tal desarrollo corresponde a la suma de los
potenciales elementales asociados a distribuciones de fuentes de fluido
de complejidad creciente (fuente nica, dipolo, etc.). De esta forma
veremos que el campo de velocidad de flujos simples, como el flujo
alrededor de un obstculo esfrico o cilndrico, se puede escribir como
una combinacin lineal del potencial de velocidad correspondiente a un
flujo uniforme y el potencial de velocidad de un dipolo.
Otra tcnica para resolver la ecuacin de Laplace consiste en buscar
soluciones en variables separadas, para lo cual es necesario utilizar
coordenadas que respeten las simetras del problema. Los problemas con
simetra cilndrica llevan a soluciones cuyas partes radiales son funciones
de Bessel; en aquellos que tienen simetra esfrica intervienen las
funciones de Legendre.

1.4.

CONDICIONES DE CONTORNO PARA FLUJOS POTENCIALES

La teora de la ecuacin de Laplace establece que la solucin est
determinada si se conoce sobre toda una superficie cerrada. Sin
embargo, es preferible en el presente contexto imponer condiciones
sobre la velocidad. Esto plantea el problema de averiguar qu condicin
sobre la velocidad es equivalente a asignar . Evidentemente, por lo
general no se podr asignar u (es decir, todas sus componentes), pues
seran tres condiciones, y no una sola como es asignar .

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 8

FLUJO POTENCIAL PLANO

Se puede demostrar que el campo u queda unvocamente determinado si

se asigna slo la componente normal de u sobre el contorno cerrado, es
decir, una condicin escalar.
En efecto, ntese que
(u) = u = u u

(6)

Pues u = 0. Entonces, si V es el volumen ocupado por el fluido,

. =

()

(7)

Cuando u es una funcin univaluada de la posicin, se puede utilizar el

teorema de la divergencia para transformar la integral de volumen. En
general, el volumen V se encuentra rodeado por fuera por la superficie
S2, y por dentro por la superficie S1, de modo que las normales n1 y
n2 son salientes del volumen (ver Fig. 1.4). Tenemos entonces:
=

. =

1 . 1 +

2 . 2 (8)

De aqu podemos deducir lo siguiente:

(a) si la componente normal de u es nula en todas las superficies del
contorno, entonces I = 0, lo cual puede suceder slo si u = 0 en todo el
volumen,
(b) si la componente normal de u no es nula en todo el contorno, es fcil
mostrar que su conocimiento es suficiente para determinar u en todo el
volumen. En efecto, supongamos que existan dos soluciones diferentes, u
= y u = , con la misma componente normal en el contorno. La
diferencia u u tiene entonces componente normal nula en el contorno
y entonces, por (a), tenemos que u u = 0 en todo el volumen, con lo
que queda demostrado que la solucin es nica.

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 9

FLUJO POTENCIAL PLANO

Fig. 1.4. Las superficies S1 y S2 limitan el flujo potencial que se desarrolla

en el volumen V. Fuente: www.google.com/potencialdevelocidades/ .

Si se asignaran, adems de la componente normal, las dos componentes

restantes, el problema quedara sobre determinado, y no tendra
solucin. En otras palabras, si se asigna un = u n sobre el contorno S, las
dos componentes tangenciales sobre S no son libres, sino que quedan
unvocamente determinadas al resolver la ecuacin de Laplace. Esto,
naturalmente, quiere decir que si parte del contorno (como suele ocurrir
frecuentemente) es una pared slida, sobre sta se tiene un = 0 en el
referencial de la pared. Pero las componentes tangenciales de u sobre la
pared, quedan determinadas por las condiciones sobre la totalidad del
contorno, de forma que por lo general no son nulas.
Es oportuno aqu un breve comentario acerca del hecho que la condicin
de contorno sobre una pared slida, esto es
un = 0

(9)

En el referencial de la pared, no determina las componentes tangenciales

de la velocidad. Del punto de vista matemtico esto no trae
inconsistencias, pues si los trminos viscosos de la ecuacin del

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 10

FLUJO POTENCIAL PLANO

movimiento son estrictamente nulos, no hay esfuerzos tangenciales. Pero

del punto de vista fsico es cuestionable, pues por pequeo que sea el
coeficiente de viscosidad, no es aceptable una discontinuidad tangencial
en la velocidad. Luego son inaceptables condiciones de contorno que
dejen libre la componente tangencial de u tan cerca como se quiera de
una pared slida. Volveremos sobre este punto cuando estudiemos los
flujos viscosos.

1.5.

FLUJOS INCOMPRESIBLES BIDIMENSIONALES Y LA FUNCIN CORRIENTE

Muchos flujos de inters son esencialmente bidimensionales, es decir,

una de las tres coordenadas espaciales es ignorable y la correspondiente
componente de la velocidad es nula (o una constante). Los flujos planos y
los flujos con simetra axial son ejemplos de esta clase de flujos.
Consideremos los flujos planos. Es habitual elegir z como la coordenada
ignorable, de modo que se puede suponer que u z = 0 y que el flujo se
desarrolla en el plano (x, y). En estos casos es til introducir el concepto
de la funcin corriente.
Volvamos por un momento al caso general de un flujo incompresible en
tres dimensiones. En virtud de la incompresibilidad, podemos siempre
introducir una funcin vectorial A tal que
u=A

(10)

El potencial vectorial de velocidad A, asociado con la condicin de

incompresibilidad, es anlogo al potencial vectorial que se introduce en el
Electromagnetismo para describir el campo magntico en virtud de que
ste cumple con la condicin B = 0. Sin embargo, en general el inters
prctico de A es escaso, dado que estamos sustituyendo el campo de
velocidad a determinar por otro campo vectorial, con lo cual no se hace
un progreso significativo. Pero en el caso de los flujos planos, s existe una
ventaja, pues el campo de velocidad tiene slo las componentes ux y uy ,

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 11

FLUJO POTENCIAL PLANO

que adems no dependen de z. Por este motivo, la nica componente no

nula de A es Az , y de acuerdo a la (10) se cumple que

; =

(11)

En coordenadas polares (r, ) tendremos

; =

(12)

La funcin se denomina funcin corriente y representa, como hemos

visto, la componente perpendicular al plano del flujo del potencial
vectorial A. Cabe observar que es posible introducir la funcin corriente
para todos los flujos bidimensionales incompresibles (sean o no viscosos y
sean o no irrotacionales), pues basta que el flujo cumpla la condicin u
= 0. Por el contrario, el potencial de velocidad slo se puede introducir
para flujos irrotacionales (que por lo tanto son necesariamente invscidos,
barotrpicos y gobernados por fuerzas de volumen conservativas), pero
no est limitado a los flujos bidimensionales.
Se puede observar asimismo que si el flujo bidimensional incompresible
es tambin irrotacional, es decir, si cumple las condiciones (2), la funcin
corriente satisface la ecuacin de Laplace

2 2 2
=
+
+
2 2 2

(13)

Al igual que el potencial de velocidad . Para ver esto basta usar las
expresiones (6.11) de las componentes de u en la condicin u = 0. El
hecho que para estos flujos tanto como satisfacen la ecuacin de
Laplace permite (como veremos ms adelante) definir un nico potencial

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 12

FLUJO POTENCIAL PLANO

complejo w = + i y aprovechar las poderosas tcnicas matemticas

basadas en las propiedades de las funciones analticas de variable
compleja para encontrar las soluciones de esos problemas.
El significado de la funcin corriente se desprende de dos propiedades
que vamos a analizar ahora. En primer lugar, las lneas = cte. coinciden
con las lneas de corriente. Esto se puede ver tomando el producto
escalar u y usando las (11):

. = ux

+ uy
=0
x
y

(14)

Este resultado implica que las lneas de corriente son ortogonales en todo
punto a . Pero por otra parte es ortogonal en todo punto a las
lneas = cte., con lo cual queda demostrada la propiedad. Cuando,
adems, el flujo es tambin potencial, las lneas = cte. son ortogonales
a las lneas equipotenciales = cte.

Fig. 1.5. La variacin de la funcin corriente es igual al caudal volumtrico que

pasa por el tubo de corriente definido por las lneas de corriente = 1 y = 2
y de espesor unidad en la direccin z. Fuente: Idem.

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 13

FLUJO POTENCIAL PLANO

En segundo lugar, la cantidad = 2 1 representa el caudal

volumtrico que pasa por el tubo de corriente de seccin rectangular,
comprendido entre las lneas de corriente = 1 y = 2 y de espesor
unidad en la direccin z (ver Fig. 1.5). Para demostrar esta propiedad,
calculemos dicho caudal. Resulta:
2

. =
1

[ + ()] =
1

[
1

+ ( )()]

= 2 1

(15)

Donde hemos usado las (11) y el hecho que si dl (dx, dy, 0) entonces n
(+dy, dx, 0) . Se verifica as tambin que en un flujo incompresible, el
caudal volumtrico de un tubo de corriente es constante a lo largo del
tubo.

1.6.

FLUJOS POTENCIALES INCOMPRESIBLES ELEMENTALES

Describiremos a continuacin algunos flujos potenciales incompresibles

de particular inters para describir flujos ms complejos.
1.6.1. Flujo uniforme
Consideremos un flujo plano uniforme con u = u0 = cte., cuya
direccin forma un ngulo con el eje x. Es inmediato ver que el
potencial de velocidad es
= 0 . = 0 cos + 0 sin

(16)

Y la funcin corriente viene dada por

= 0 = 0 cos 0 sin

(17)

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 14

FLUJO POTENCIAL PLANO

Las equipotenciales y las lneas de corriente son familias de rectas

paralelas, perpendiculares entre s.
1.6.2. Fuente o sumidero puntiforme
Se describe por un campo de velocidad del tipo
0
=
, = 0
(18)
3
Donde r0 es el punto donde se encuentra la fuente (o sumidero),
o del cual salen (o entran) las lneas del campo de velocidad 4.
Mediante una traslacin de coordenadas que ubique el origen en
r0, y la eleccin de coordenadas esfricas (r, , ) centradas en el
nuevo origen, el campo en cuestin tiene como nica componente
no nula la componente radial, cuyo valor es
=

(19)

Que satisface las condiciones (6.1) excepto en r = 0.

El flujo de masa que atraviesa una esfera centrada en el origen es
4 2 = 4 (20)

La cantidad Q es el flujo volumtrico o caudal erogado por la

fuente (o evacuado por el sumidero) y se denomina tambin
fuerza de la fuente (o sumidero). En trminos de Q, la (19) se
escribe como
=

4 2

(21)

El potencial de velocidad de la fuente (o sumidero) se expresa

como
=

(22)

1.6.3. Fuente o sumidero lineal

4

En el punto r0, naturalmente, est excluido la regin del flujo potencial.

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 15

FLUJO POTENCIAL PLANO

En este caso la fuente (o sumidero) est distribuida a lo largo de

una recta. Emplearemos un sistema de coordenadas cilndricas (r,
, z) cuyo eje es dicha recta. El flujo es plano y la nica
componente no nula del campo de velocidad es la componente
radial, que vale
=

(23)

Y que cumple con las condiciones (1) excepto en r = 0 (regin

excluida). La fuerza de la fuente se define ahora por unidad de
longitud y es igual al caudal volumtrico a travs de un tronco de
cilindro de altura unitaria (y radio arbitrario), de modo que
= 2. El potencial de velocidad es (recordemos que u r = /r
y u = r1/ )

= ln( ) =
ln( )
0
2 0

(24)

Donde r0 es una constante arbitraria5. La funcin corriente es

=

(25)

Por lo tanto las lneas de corriente son rectas radiales y las curvas
equipotenciales son circunferencias con centro en el origen (ver

El potencial de velocidad est definido a menos de una constante arbitraria.

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 16

FLUJO POTENCIAL PLANO

Fig. 1.6.3a).
Fig. 1.6.3. Campos de velocidad de: (a) una fuente lineal, (b) un vrtice lineal. Las
lneas de corriente se indican con lneas llenas y las equipotenciales con lneas de
trazos. Las diferencias y entre equipotenciales y lneas de corriente
consecutivas es la misma. Fuente: Idem.
1.6.4. Vrtice lineal
Consideremos un flujo plano cuyo campo de velocidad est dado
en un sistema de coordenadas cilndricas (r, , z) por
= 0 ; =

; = 0

(26)

Este campo se suele denominar campo de circulacin, puesto que

sus lneas de corriente son circunferencias con centro en el eje. Se
puede demostrar fcilmente (recurriendo a las expresiones de la
divergencia y del gradiente en coordenadas cilndricas) que
cumple con las condiciones (1) excepto sobre el eje r = 0. Sin
embargo, para demostrar que u = 0 es quizs ms intuitivo
aplicar el teorema de Stokes a un circuito cerrado formado por
dos arcos de circunferencia que abarcan un ngulo , de radios
son r1 y r2, y dos segmentos radiales que unen ambos arcos. La
circulacin es entonces
= 2 2 1 (1 ) = 0

(27)

Pues los tramos radiales no contribuyen por ser perpendiculares a

la velocidad. Resulta claro que toda curva cerrada que no rodea al
eje, se puede aproximar mediante tramos radiales y arcos de
crculo, por lo cual es fcil mostrar que para todas esas curvas se
tiene = 0 , lo cual en virtud del Teorema de Stokes implica que
u = 0 en todo el espacio, salvo en el eje mismo. Por otra parte, es
sencillo demostrar que la circulacin alrededor de toda lnea
cerrada que rodea (una sola vez) al eje vale
= 2

(28)

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 17

FLUJO POTENCIAL PLANO

La cantidad se conoce con el nombre de fuerza del vrtice. El

potencial de velocidad es
= =

(29)

La funcin corriente es
=

ln( )
2 0

(30)

En este caso las lneas de corriente son circunferencias con centro

en el origen y las equipotenciales son rectas radiales (Fig. 1.6.3b).
Notar la correspondencia entre este flujo y el de la fuente (o
sumidero) lineal: al pasar de uno a otro las funciones y
intercambian sus roles conservando su forma. Las lneas de
corriente radiales de la fuente lineal son idnticas a las
equipotenciales del vrtice lineal. Las lneas de corriente circulares
del vrtice correspondan a las equipotenciales de la fuente lineal.

1.6.5. Dipolo
Consideremos un sumidero S 1 en r1 y una fuente S 2 en r2 = r1 + d,
cuyos caudales volumtricos tienen el mismo valor absoluto Q (ver
Fig. 1.6.5. - 1). Si tomamos el lmite |d| 0, manteniendo
constante el vector
p = Qd

(31)

Obtenemos lo que se denomina flujo dipolar de momento dipolar

p.
Calculemos ahora el campo de velocidad del flujo dipolar en el
caso plano (dipolo lineal). Para simplificar tomemos origen de
coordenadas en S1. Entonces, los potenciales de velocidad
inducidos por la fuente y el sumidero en un punto P de
coordenadas r son, respectivamente
2 =

ln(
) ; 1 = ln( )
2
0
2
0

El potencial del flujo resultante es

(32)

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 18

FLUJO POTENCIAL PLANO

= 2 + 1 =

ln(
)
ln( ) =
ln
(33)
2
0
2
0
2

Fig. 1.6.5. 1. Un sumidero S1 en r1 y una fuente S2 en r2 = r1 + d, cuyos

caudales volumtricos tienen el mismo valor absoluto Q. Si tomamos el lmite
|d| 0, manteniendo constante el valor de p = Qd, se obtiene un dipolo.
Fuente: Idem.
Como d es muy pequeo podemos escribir

2 2 + 2
=

1
2

=1

+ (34)
2

Donde con indicamos trminos de orden superior. Por lo tanto

ln

= ln(1 2 ) + = 2 + (35)

Y en el lmite |d| 0 resulta

=
,
= lim (36)
0
2 2
Usando coordenadas polares en el plano del flujo y tomando p
como origen de , resulta

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 19

FLUJO POTENCIAL PLANO

cos
(37)
2

A partir de podemos calcular las componentes de la velocidad:

=

cos
1 sin
=
;

=
=

2 2

2 2

(38)

A partir de las componentes de la velocidad, usando las relaciones

(12) o la (19) podemos obtener la funcin corriente:
=

sin
=
2 2
2 2

(39)

Es fcil verificar que las lneas de corriente (Fig. 1.6.5. - 2) son una
familia de circunferencias tan- gentes a p y que pasan por el
origen, y

que las

equipotenciales

son

una

familia

de

circunferencias ortogonales a las primeras, que pasan por el

origen y son tangentes a la recta = 2 .

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 20

FLUJO POTENCIAL PLANO

Fig. 1.6.5. 2. Campo de velocidad de un dipolo lineal. Las lneas de corriente se

indican con lneas llenas y las equipotenciales con lneas de trazos. Fuente: Idem.
En el caso tridimensional se puede llevar adelante el clculo de
manera anloga. En coordenadas esfricas (r, , ) con origen en
el dipolo y eje polar en la direccin de p, resulta
=

.
cos
=
(40)
3
4
4 2

Las componentes de la velocidad son

=

cos
,
2 3

sin
,
2 3

= 0 (41)

Puesto que el flujo dipolar en tres dimensiones tiene simetra axial

alrededor del eje del dipolo, se trata en realidad de un flujo
bidimensional pues la coordenada es ignorable. Por lo tanto el
problema se puede tratar en un plano = cte. Introduciendo un
potencial vectorial de velocidad A por medio de la (10) y usando
las expresiones del rotor en coordenadas esfricas, es fcil ver que
la nica componente no nula de A es A . Por lo tanto, es posible

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 21

FLUJO POTENCIAL PLANO

tambin en este caso introducir una funcin corriente (r, ),

definida por
=

1
,
2 sin

1
(42)
sin

Usando las expresiones (42) de las componentes de la velocidad se

obtiene
=

sin2
4

43

La lnea de corriente y las curvas equipotenciales para el flujo

dipolar se pueden apreciar en la Fig. 1.6.5. 3.

Fig. 1.6.5. 3. Campo de velocidad de un dipolo. Las lneas de corriente se

indican con lneas llenas y las equipotenciales con lneas de trazos. Fuente: Idem.

1.7.

Fuerzas sobre un obstculo en un flujo potencial

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 22

FLUJO POTENCIAL PLANO

Trataremos ahora el problema de calcular la fuerza que se ejerce sobre

un obstculo de forma arbitraria en un flujo potencial. Supondremos que
a distancias grandes del obstculo el flujo es uniforme con una velocidad
U. La solucin se basa en el desarrollo multipolar del potencial de
velocidad, tomando en cuenta las condiciones de contorno sobre el
obstculo (nulidad de la componente normal de la velocidad, es decir un
= 0 en el referencial en que el obstculo est en reposo). A distancias
grandes del obstculo bastar con tomar en cuenta la contribucin del
primer trmino no nulo del desarrollo multipolar, pues es el que tiene la
potencia ms baja de 1/ r. Conviene tratar por separado los casos
bidimensional y tridimensional.

1.7.1. Caso bidimensional

El cilindro circular y el perfil de un ala de avin de longitud infinita
son ejemplos clsicos de este problema. Supondremos que el
objeto est en reposo, como sera el caso de estudiar el flujo en un
tnel de viento en el cual se impone la velocidad U, uniforme a
gran distancia del obstculo, y perpendicular a la coordenada
ignorable z. Vamos a suponer a priori que hay una circulacin no
nula alrededor del obstculo, sin preocuparnos por ahora de su
origen. Eso se ver ms adelante cuando estudiemos los efectos
de la viscosidad y la capa lmite.

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 23

FLUJO POTENCIAL PLANO

Fig. 1.7.1. Obstculo cilndrico cuya seccin tiene forma arbitraria, que est
inmerso en un flujo potencial uniforme y que tiene una circulacin atrapada.
Fuente: Idem.
Tomaremos un sistema de coordenadas cuyo origen est centrado
en el obstculo (Fig. 1.7.1).
Claramente, el potencial de velocidad tiene la forma
= + + (43)
Donde U es el potencial del flujo uniforme, es el potencial del
vrtice que da cuenta de la circulacin alrededor del obstculo, y
M es un desarrollo multipolar centrado en el obstculo, cuyos
trminos estn determinados por la condicin de contorno u n = 0
sobre la superficie del slido. En general, no habr ni fuentes ni
sumideros en el slido, de modo que
= + =

.
+ (44)
2

El primer trmino del miembro derecho de (44) es el potencial de

un dipolo de momento 2A (ver la ecuacin (36)) y los puntos

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 24

FLUJO POTENCIAL PLANO

suspensivos representan multipolos de orden superior que dan

contribuciones a M del orden de r2, r3, etc. El vector A es
constante y su valor depende de la forma del objeto y de su
orientacin con respecto de U. En el caso del cilindro circular los
trminos de orden superior son nulos, pues el trmino dipolar
alcanza para obtener la solucin exacta. Para un obstculo de una
forma cualquiera, sin embargo, no sern nulos y es necesario
calcularlos si se desea conocer con exactitud el campo de
velocidad. Veremos, sin embargo, que para calcular la fuerza
sobre el obstculo es suficiente emplear la aproximacin dipolar
de M, es decir A. De esta forma tenemos:
cos +

.
+ 2 (45)
2

Donde hemos tomado el eje polar en la direccin de U. La

velocidad es entonces
= = cos

+ sin +
+ 2 (46)
2

Vamos a plantear ahora el balance de la cantidad de movimiento

del fluido dentro de un volumen cilndrico 6 de radio r y altura
unidad, con r grande respecto de la dimensin lineal L del
obstculo (r >> L). Usando la ecuacin que se muestra tendremos
que la resultante R de las fuerzas que actan sobre el fluido
contenido en el cilindro (volumen fijo en el espacio) est dada por
=

[ . + ] (47)

Donde n es la normal a la superficie lateral S del cilindro (no hay

flujo a travs de las caras inferior y superior), hemos ignorado los
efectos de la viscosidad, y hemos supuesto que el flujo es
estacionario y que no hay fuerzas externas de volumen. La (47) es

De esta forma evitamos el problema de integrar el campo de presin sobre la superficie del
obstculo.

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 25

FLUJO POTENCIAL PLANO

la fuerza que el obstculo ejerce sobre el fluido. Por lo tanto, la

fuerza que el fluido ejerce sobre el obstculo 7es
= =

. + =

( + )

(48)

La presin est dada por la ecuacin de Bernoulli

1
= + 2 2 (49)
2
Usando la expresin de la velocidad (46) en (49) resulta
= +

sin + (50)
2

Donde los puntos suspensivos indican trminos del orden de r2 o

superior, que no contribuyen a F en el lmite r . En cuanto al
trmino u ur de la integral (49), tenemos que
= 2 cos2 + 2 cos sin +

cos +
2

(51)
Donde los puntos suspensivos indican trminos del orden de r2 o
superior, que no contribuyen a F en el lmite r . Por lo tanto,
usando la (6.92) donde omitimos el trmino p er (dado que su
contribucin a F es nula), y la (51), el integrando de la (48) queda
como
+

sin
2

+ 2 cos sin +
cos
2

= 2 cos +
(52)
2
= 2 cos2 +

Para que el obstculo permanezca inmvil debemos pensar que mediante algn agente externo

(por ejemplo, un soporte) ejercemos sobre el objeto una fuerza que equilibra a F.

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 26

FLUJO POTENCIAL PLANO

A menos de trminos del orden de 1/r. La integral se puede

calcular de inmediato y resulta
= 0 ; = (53)
Por lo tanto, el arrastre es nulo, un resultado que refleja la
ausencia de disipacin de energa en los flujos potenciales
estacionarios. En cuanto a la sustentacin (Fy), depende
nicamente de la presencia de una circulacin de la velocidad
alrededor del obstculo. La sustentacin, tambin llamada fuerza
de Magnus, se puede expresar en forma vectorial utilizando un
vector circulacin paralelo a la direccin z, de modo que
= ; = (54)
Para aplicar correctamente la frmula (6.96) hay que respetar la
convencin de calcular recorriendo una curva cerrada alrededor
del obstculo en sentido anti horario ( creciente, por lo tanto es
negativo.
La fuerza de Magnus es responsable de la sustentacin del ala de
un avin, y la geometra del perfil alar es lo que origina la
aparicin de la circulacin necesaria. Esa fuerza tambin es la base
de los mecanismos de propulsin de embarcaciones y aviones y de
la sustentacin de los helicpteros mediante hlices. Tambin se
la invoca para explicar la trayectoria de una pelota de tenis o de
ping-pong cuando se le ha impartido efecto, e igualmente en el
caso de la pelota de ftbol cuando se la patea con chanfle. Esto
es cierto, pero cabe observar que no puede haber circulacin
atrapada en el flujo potencial alrededor de un obstculo todas
cuyas dimensiones son finitas. En efecto, sea cual fuere el circuito
C de integracin que se elija, es siempre posible encontrar una
superficie S que se apoye sobre C y que est ntegramente
contenida en la regin del flujo potencial. Pero la circulacin de u
alrededor de C es igual al flujo de la vorticosidad = a

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 27

FLUJO POTENCIAL PLANO

travs de S, y en un flujo potencial, = 0. Por consiguiente la

circulacin de u al- rededor de C resulta idnticamente nula. Lo
que produce la circulacin en los casos menciona- dos (pelotas de
tenis, ping-pong, ftbol, etc.) es el arrastre del fluido cercano a la
superficie del objeto que est en rotacin, debido a la viscosidad
del fluido8. En esos casos se puede tambin aplicar la ecuacin de
Bernoulli, suponiendo la continuidad de la presin a travs de la
capa lmite prximo a la superficie del cuerpo, donde el flujo no es
potencial. De esta forma se encuentra que la fuerza sobre el
cuerpo es perpendicular a U, al igual que para la geometra
cilndrica.
1.8.

Comentarios sobre los flujos potenciales

Antes de concluir este Captulo es oportuno comentar algunas cuestiones
que seguramente intrigan al lector. Parte de ellas se estudiarn en
Captulos posteriores, pero no est dems una discusin preliminar.
Claramente la teora ideal que presentamos aqu, si bien predice
correctamente la sustentacin de un ala, no est de acuerdo con la
experiencia en ciertos aspectos y no aclara satisfactoriamente otros. Por
ejemplo predice que el arrastre sobre un cuerpo que se mueve en el seno
de un fluido es nulo, lo cual est en violenta contradiccin con la
experiencia cotidiana.
Cmo conciliar este fracaso con la correcta prediccin de la
sustentacin? Adems, cmo se origina la circulacin que produce la
sustentacin del ala? De acuerdo con los teoremas de la circulacin, si el
ala parte del reposo y el aire est inmvil, y entonces claramente la
circulacin de la velocidad alrededor del perfil es nula y no hay
sustentacin, parecera a primera vista que tendra que seguir siendo nula
cuando el ala se pone en movimiento.

Por este motivo es esencial la rotacin de la pelota para producir el efecto.

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 28

FLUJO POTENCIAL PLANO

La pista que permite aclarar estas cuestiones parte de considerar los

efectos de la viscosidad, que hasta aqu hemos ignorado. Pero la forma
como interviene la viscosidad en muchos flujos es muy sutil y para nada
obvia como se ver en seguida.
1.8.1. Flujos viscosos
Para comenzar a entender el rol de la viscosidad conviene regresar
por un momento a la ecuacin de Navier-Stokes

+ . =
+ 2 (55)

Esta ecuacin difiere de la ecuacin de Euler, sobre la cual se basa

la teora de los fluidos idea- les, por la presencia del trmino
viscoso v2u (aqu = es la viscosidad cinemtica).
Por de pronto, la observacin de los flujos reales (o sea, viscosos)
muestra que sobre un con- torno rgido, no slo la componente
normal, sino tambin la componente tangencial de u deben ser
iguales a la velocidad del contorno 9. Si el contorno est en reposo,
debemos entonces tener
u=0

sobre el contorno

(56)

Esta condicin sobre la componente tangencial de la velocidad se

suele llamar condicin de no deslizamiento, y vale para todo fluido
con v 0, por pequea que sea la viscosidad.
En este contexto, podemos observar que usando la identidad
vectorial
2A = ( A) ( A)

(57)

Y recordando que u = y que en un flujo incompresible u =

0, tenemos que
2 u = (58)

El hecho que las parcelas del fluido que estn en contacto con la superficie de un slido se
tienden a adherir a l se puede entender cualitativamente en trminos de las interacciones
intermoleculares.

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 29

FLUJO POTENCIAL PLANO

Por lo tanto 2u = 0 en un flujo cuya vorticosidad es nula (en

realidad basta que w = 0) y la ecuacin de Navier-Stokes
coincide con la ecuacin de Euler, pese a la existencia de esfuerzos
tangenciales. Este hecho, sobre el cual volveremos ms adelante,
explica porqu en un fluido real es posible tener un flujo potencial
en la mayor parte del volumen, salvo en determinadas regiones
donde por algn motivo hay vorticosidad. Claramente, cerca de
superficies slidas, cuando la teora ideal predice una componente
tangencial no nula de la velocidad relativa, los efectos de la
viscosidad son seguramente importantes, pues es precisamente
debido a la viscosidad que se consigue que dicha componente se
anule. En tales regiones hay vorticosidad y por lo tanto el flujo no
es potencial y tampoco, obviamente, se puede aplicar la ecuacin
de Euler.
1.8.2. El nmero de Reynolds
Consideremos un fluido viscoso en movimiento, y sea U el valor
tpico de la velocidad. Adems, sea L la escala espacial
caracterstica del flujo 10. Una vez elegidos L y U, podemos formar
la cantidad
=

(59)

Que es un nmero puro que se conoce como nmero de Reynolds.

El motivo de la importancia del nmero de Reynolds se puede
apreciar si se observa que las derivadas de las componentes de u
son tpicamente del orden de magnitud de U / L (suponiendo, se
entiende, que las componentes de u tengan variaciones del orden
de U sobre distancias del orden de L). Tpicamente, las derivadas
de las componentes de u tendrn ellas mismas variaciones del
orden de U / L sobre distancias del orden de L, luego derivadas
10

Un ejemplo puede permitir aclarar esto. En el caso del aquietamiento de la taza de t, un

valor razonable para L sera 4 cm, y para U, 5 cm/s.

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 30

FLUJO POTENCIAL PLANO

segundas como 2ui/x2,, etc., sern del orden de magnitud de U

/ L2. De esta forma podemos estimar el orden de magnitud de los
trminos de inercia y viscoso de la ecuacin de Navier-Stokes
como
Trmino de inercia: ( . ) = ( 2/)
Trmino viscoso: 2 = ( 2 /) (60)
Siempre y cuando estas estimaciones sean correctas, se deduce
que

2
=
= (61)

2
Por consiguiente, el nmero de Reynolds es importante porque
permite estimar groseramente la importancia relativa de dos
trminos importantes de la ecuacin de Navier-Stokes. No debe
sor- prender, entonces, que flujos caracterizados por un nmero
de Reynolds grande y flujos caracterizados por un nmero de
Reynolds pequeo tengan caractersticas muy diferentes.

1.8.3. Flujos de alto nmero de Reynolds

El caso Re >> 1 corresponde a lo que podramos llamar flujos de
baja viscosidad. El resultado (61) sugiere que cuando Re >> 1 los
efectos de la viscosidad deberan ser, en general, despreciables. El
flujo alrededor de un perfil alar delgado, cuando el ngulo de
ataque es pequeo, es un caso de esta clase. An as, sin embargo,
los efectos de la viscosidad se vuelven importantes en capas lmite
delgadas, donde existen gradientes de velocidad muy grandes.
Debido a ello, en esas capas delgadas, el trmino viscoso resulta
mucho mayor de lo que indica la estimacin (61), que se basa
sobre las escalas caractersticas del grueso del flujo. Veremos ms
adelante que el espesor tpico de tales capas lmite est dado
por = 1

(62)

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 31

FLUJO POTENCIAL PLANO

Por lo tanto, cuanto mayor sea el nmero de Reynolds del grueso

del flujo, tanto ms delgada resulta la capa lmite.
Est claro que para poder aplicar al grueso del flujo la teora de
flujos ideales que desarrollamos en este Captulo, es necesario que
el nmero de Reynolds sea grande (de lo contrario la capa lmite
abarcara el grueso del flujo). Sin embargo esta condicin no es
suficiente. Veremos, en efecto, que el fenmeno de la separacin
de la capa lmite, que consiste en que la capa lmite se
desprende de la superficie del slido, lleva a un estado de cosas
muy diferente del que predice la teora ideal. Existe, adems, otra
complicacin en los flujos de alto nmero de Reynolds. Ocurre
muchas veces, en efecto, que las soluciones estacionarias de la
ecuacin

de

Euler

son

inestables

frente

pequeas

perturbaciones, y de resultas de eso el flujo deja de ser

estacionario y se vuelve turbulento. En realidad, fue precisamente
en este contexto que Reynolds introdujo el parmetro a
dimensional que hoy lleva su nombre.

Fig. 1.8.3. - 1. Arrastre sobre un cilindro circular en un flujo uniforme.

Fuente: Idem.

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 32

FLUJO POTENCIAL PLANO

La Fig. 1.8.3. 1. Muestra el coeficiente de arrastre sobre un cilindro

circular, en funcin del nmero de Reynolds. Se puede observar que el
arrastre sigue siendo importante en todo el intervalo entre Re = 102 y
107, que equivale a disminuir la viscosidad por un factor 105. Adems,
como lo indican los esquemas, el patrn general del flujo no muestra
signos de tender v 0 (o sea Re ). Esto se debe a que para lograr esa
forma, cerca de la superficie del cilindro la velocidad del flujo principal
debera disminuir muy rpidamente del lado de atrs del mismo, y por lo
tanto habra un fuerte gradiente de presin en sentido contrario al flujo
(lo que se llama un gradiente adverso). Pero ocurre que la capa lmite
adherida a la superficie del cilindro no puede soportar ese gradiente
adverso. El motivo es que debido al gradiente adverso, muy cerca de la
superficie se produce una inversin del flujo, que hace que la capa lmite
(en la cual el flujo es rotacional) se separe de la superficie del cilindro
provocando una importante estela detrs del mismo, en la cual el flujo es
vorticoso. Si bien el carcter de la estela cambia al aumentar el nmero
de Reynolds, sta nunca desaparece. La presencia de la estela es
responsable del arrastre que se observa en los experimentos.

Fig. 1.8.3. 2. Generacin de la sustentacin de un ala. El patrn (1) del flujo en

el instante que el ala se pone en movimiento a partir del reposo no tiene
circulacin. Para un perfil alar diseado para tener sustentacin, cuando el ala se

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 33

FLUJO POTENCIAL PLANO

comienza a mover el flujo tiene un punto de estancamiento sobre la superficie

superior, antes del borde de fuga, y la velocidad es singular en el filo del borde
de fuga. Tal situacin se mantendra siempre si no fuera por los efectos de la
viscosidad, que hacen que la capa lmite de la parte de abajo del ala se separe
justo en el borde de fuga. El vrtice que as se produce interacta con la regin
no viscosa del flujo haciendo que el punto de estancamiento superior se
desplace hacia el punto de separacin. El nuevo flujo (2) no viscoso (fuera de la
capa lmite) tiene ahora una circulacin, cuya magnitud est determinada por la
geometra del perfil alar, junto con la condicin que el punto de estancamiento
posterior coincida con el borde de fuga. El pro- ceso por el cual se genera la
sustentacin tiene lugar gracias a que se descarga vorticosidad en el fluido (3).
De acuerdo con el Teorema de Helmholtz, la circulacin alrededor del circuito (a
e c d a) que rodea el ala est compensada exactamente por la circulacin
alrededor del circuito (a b c e a) que rodea la estela y el vrtice pero no el ala.
A diferencia del cilindro, un perfil alar est especialmente diseado para
evitar la separacin de la capa lmite. La caracterstica esencial que
permite lograr esto es el perfil alargado y afilado de la parte posterior,
que sirve para reducir el gradiente adverso de presin, haciendo que la
presin aumente en forma muy gradual en la parte del ala donde el
espesor disminuye. Por este motivo, la capa lmite no se separa sino muy
cerca del borde de fuga, y entonces se produce una estela sumamente
delgada. En este caso, s es posible obtener valores muy pequeos del
arrastre para cuerpos bien perfilados (aerodinmicos) el coeficiente de
arrastre es tpicamente del orden de Re1/2 para Re . Este estado
de cosas se mantiene siempre y cuando el ngulo de ataque sea pequeo.
Si el ngulo de ataque supera cierto lmite, el aumento de presin sobre
la parte superior del ala deja de ser suficientemente gradual, y se produce
la separacin de la capa lmite muy adelante del borde de fuga. Cuando
eso ocurre se dice que el ala est en prdida (de sustentacin), en efecto,

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 34

FLUJO POTENCIAL PLANO

la sustentacin cae bruscamente a la par que aumenta fuertemente el

arrastre.
Las propiedades de la capa lmite tambin permiten entender el origen de
la circulacin necesaria para la sustentacin del ala. Cuando el ala y el aire
que la rodea estn en reposo, la velocidad y la vorticosidad son nulas en
todas partes. Los experimentos muestran que al ponerse en movimiento
el ala, se genera vorticosidad en la capa lmite (Fig. 1.8.3. 2.); esta
vorticosidad queda confinada en esa capa y en una delgada estela que
queda atrs del ala, y que termina en un vrtice que se suele llamar
vrtice de arranque. Consideremos ahora un circuito material (a b c d a)
suficientemente grande como para estar ntegramente afuera de esas
regiones desde el inicio del movimiento. Puesto que el estado inicial es de
reposo, la circulacin alrededor de ese circuito es nula inicialmente y por
el Teorema de Kelvin sigue siendo nula en todo tiempo posterior. Si ahora
dividimos el circuito en dos partes, una de los cuales (a b c e a) rodea el
ala pero no la es- tela y el vrtice de arranque, y la otra (a e c d a) rodea la
estela y el vrtice pero no el ala, clara- mente la circulacin alrededor de
estos dos circuitos debe ser de igual valor absoluto y signo opuesto.
Ahora bien, el diseo del perfil alar est concebido de modo que la
vorticosidad que se desprende del ala es positiva, para que quede
alrededor del ala una circulacin negativa, que proporcione la
sustentacin. El desprendimiento de vorticosidad contina hasta que la
circulacin alrededor del perfil alar alcanza el valor necesario para que el
flujo principal (irrotacional) no presente singularidad en el borde de fuga.
Este valor es precisamente el K de Kutta- Joukowski. A partir de ese
momento no hay ms variacin neta de la vorticosidad de las capas lmite
de la superficie superior e inferior del ala, y se mantiene el valor final K
de la circulacin.

FACULTAD DE INGENIERIA DE MINAS, GEOLOGIA Y METALURGIA 35

FLUJO POTENCIAL PLANO

BIBLIOGRAFIA

Pginas web usadas:

http//www.google.com.pe (se utilizo el buscador como potencial de

velocidad)

http//www.wikipedia.com (potencial de flujo)