30

Captulo 3

An
alisis Diferencial
En este captulo se desarrollaran las ecuaciones diferenciales que describen el movimiento de un
fluido real. En particular se desarrollan el principio de conservacion de masa y la segunda ley de
movimiento de Newton o principio de conservacion de cantidad de movimiento lineal (F~ = m~a)
para un elemento diferencial. Lo que se obtiene es un conjunto de ecuaciones diferenciales
parciales. Estas ecuaciones son difciles de resolver en forma analtica por lo que generalmente
se recurre a soluciones numericas. Este es el ambito de estudio de la Mecanica de Fluidos
Computacional o CFD, el cual no ser
a tratado en este curso. Sin embargo lo anterior, algunos
ejemplos sencillos pueden ser resueltos en forma analtica. Se analizaran ademas, algunos casos
particulares de flujo, para los cuales las ecuaciones se simplifican y de cuyo analisis se pueden
obtener conclusiones importantes. En una primera parte se analizara la cinematica de una
partcula elemental.

3.1

Cinem
atica

La cinematica estudia varios aspectos de un fluido en movimiento como velocidad, posicion y

aceleracion sin analizar las fuerzas necesarias para que se produzca dicho movimiento.
En una primera parte describiremos el movimiento en terminos del movimiento de una partcula
fluida y posteriormente se realizara un analisis macroscopico para la descripcion de un flujo.
La descripcion de cualquier propiedad del fluido puede ser descrita como una funcion de su
posicion. En particular se utilizan coordenadas espaciales (x, y, z por ejemplo) para identificar
las partculas de fluido y sus propiedades. Esta representacion se denomina representacion de
~ =V
~ (x, y, z). Como la
campo. As por ejemplo, el campo de velocidades vendra dado por V
representaci
on de una partcula puede ser diferente en tiempos diferentes la representacion debe
ser tambien una funcion del tiempo. Para el campo de temperaturas y velocidades por ejemplo
T = T (x, y, z, t)
~ =V
~ (x, y, z, t)
V
~ = u(x, y, z, t) + v(x, y, z, t)
V
+ w(x, y, z, t)k
El movimiento de una partcula puede ser descrito en terminos de la velocidad y la aceleracion.
Por definicion la velocidad de una partcula es la variacion temporal del vector posicion
~ = d~r .
V
dt

31

3.1 Cinematica

~ |. Si las propiedades de un flujo, en todos los puntos

La rapidez es el modulo de la velocidad |V
del espacio, permanecen invariantes en el tiempo, se dice que el flujo es permanente. En caso
contrario se llama nopermanente. Un campo de velocidades permanente estara dado por
~ =V
~ (x, y, z) .
V
Para el caso permanente se cumple (/t = 0).
Como se menciono anteriormente las propiedades de un fluido y las caractersticas del flujo se
pueden representar como una funcion de la posicion y del tiempo. Se desprende de lo anterior
dos formas de posibles de representacion:
1. La primera, utiliza el concepto de campo mencionado anteriormente. La descripcion del
flujo esta dada por la descripcion de las propiedades de este como una funcion de la posici
on
y del tiempo. De esta manera se obtiene informacion del flujo en terminos de que pasa
en un punto fijo del espacio en un tiempo t cuando el flujo pasa por el. Este metodo de
descripcion se denomina descripcion Euleriana. La velocidad queda representada por el
campo de velocidades dado por
~ =V
~ (x, y, z, t) .
V
2. El segundo metodo, denominado Lagrangiano, analiza una partcula generica del flujo para
analizar y caracterizar el flujo. En esta representacion la posicion x, y, z no son fijas sino
varan en el tiempo. Las coordenadas espaciales seran por lo tanto funciones del tiempo
y de una posicion preescrita xo , yo , zo en un instante to . Para la velocidad se tiene por lo
tanto
~ =V
~ (x(t), y(t), z(t), t) .
V
Para representar el flujo en forma grafica se utiliza el concepto
de lnea de corriente. Las lneas de corriente son las envolventes
de los vectores de velocidad de las partculas fluidas, es decir,
el vector de velocidad es siempre tangente a las lneas de corriente. Si el flujo es permanente (/t = 0) las lneas de corriente
estaran fijas en el tiempo y coincidiran con la trayectora de las
partculas. Si el flujo no es permanente (/t 6= 0) las lneas de
corriente seran solo una representacion instantanea del flujo.
Se llama tubo de corriente al conjunto de lneas de corriente que
pasan por el contorno de un area A, en un tiempo determinado.
Dado que la velocidad es tangente a las lneas de corriente, no
existira flujo a traves del manto de un tubo de corriente por lo
que se cumple que

A
Lneas y tubo de corriente

~ d~r = 0 ,
V
~ . De
donde d~r es el desplazamiento diferencial de una partcula fluida que tiene una velocidad V
la ecuacion anterior resulta
dy
dz
dx
=
=
,
u
v
w

(3.1)

que representan las ecuaciones para determinar las lneas de corriente.

C. Gherardelli

U. de Chile

32

3.1 Cinematica

3.1.1

Velocidad, Rotaci
on y Deformaci
on

Como se menciono en el captulo 1 un fluido es una substancia que se deforma al aplicar sobre
esta un esfuerzo de cizalle. Debemos esperar por lo tanto que una partcula fluida se encuentre
sometida a movimientos de traslacion, rotacion, deformacion lineal y angular como se muestra
en la figura 3.1. Este tipo de movimientos esta asociado a variaciones complejas de las diferentes
componentes de la velocidad (u, v, w) en todas las direcciones. Lo anterior nos indica que, en
general, (Vi /xj ) 6= 0 i, j. Analizaremos ahora cada uno de estos efectos por separado y su
relacion con la variacion de la velocidad seg
un los distintos ejes coordenados.

Figura 3.1: Superposicion de movimientos de una partcula fluida.

Traslaci
on
El movimiento mas sencillo al cual se puede encontrar sometida una partcula fluida es el
movimiento de traslacion. En la figura 3.2 se muestra una partcula que viaja con una velocidad constante desplazandose desde su posicion original una nueva posicion dada por los
puntos O0 A0 C 0 B 0 .

Figura 3.2: Movimiento de traslacion de una partcula fluida.

Deformaci
on lineal
Analizaremos la deformacion lineal seg
un el eje x, como se muestra en la figura 3.3. Nos interesa
por lo tanto analizar la variacion de la velocidad seg
un el mismo eje, es decir (u/x). Como
se muestra en esta figura, y debido a la diferencia de velocidad existente entre las lneas OB y
AC el elemento de fluido se deforma en un tiempo t. La variacion del volumen resulta


V =

u
x yzt .
x

C. Gherardelli

U. de Chile

33

3.1 Cinematica

Figura 3.3: Deformacion lineal de una partcula fluida.

El cambio de volumen por unidad de volumen es

1 d(V )
(u/x)t
u
= lim
=
.
t0
V t
t
x


El cambio de volumen, por unidad de volumen, seg

un todos los ejes es la superposicion de los
cambios seg
un cada eje y por lo tanto igual a
1 d(V )
u v w
~.
=
+
+
=V
V t
x y
z
~ , se encuentra asociada a la deformacion lineal
Vemos como la divergencia de la velocidad, V
de la partcula fluida. Como un cambio de volumen a masa constante significa una variacion de
la densidad, se debe cumplir que
~ = 0 para un flujo incompresible y
V
~ 6= 0 para un flujo compresible.
V
Rotaci
on
La velocidad angular de la lnea OA, OA , de la figura 3.4 queda definida por

.
t0 t

OA = lim

Para angulos peque

nos se tiene que la tangente del angulo se puede aproximar por el valor del
angulo. Analizando la figura y para peque
nos se tiene que
tan =

v
x xt

v
t ,
x

OA =

v
.
x

Analogamente la velocidad angular de la linea OB, OB , resulta

u

=
.
t0 t
y

OB = lim

C. Gherardelli

U. de Chile

34

3.1 Cinematica

Figura 3.4: Rotacion de una partcula fluida.

La velocidad angular en torno al eje z, z , se define como el promedio aritmetico de OA y

OB , es decir:

1
(OA OB )
2

1 v
u

2 x y

z =
=

(3.2)

Se puede observar de la ecuacion 3.2 que la partcula fluida rotara en torno al eje z como un
cuerpo rgido, es decir sin deformacion, solo si u/y = v/x. En otro caso la rotaci
on
estara asociada a una deformacion. Se ve ademas que cuando u/y = v/x la rotacion en
torno al eje z es cero.
Para los otros ejes se obtiene
1
2

u w

z
x

1
x =
2

w v

y
z

y =

De las ecuaciones anteriores se puede ver que:

~ = 1 V
~

2
~ = 0 se llama irrotacional y representa un tipo especial de flujo como
Un flujo para el cual V
se vera mas adelante.
La vorticidad
~ de un flujo se define como
~ =V
~

~ = 2
Deformaci
on angular
Se ve de la figura 3.4 que las derivadas u/y y v/x pueden causar, ademas de la rotaci
on
de la partcula, una deformacion. La tasa de deformacion angular de una partcula se mide
por la rapidez de cambio del angulo que se forma entre las lneas OA y OB. Si OA gira a un
velocidad angular distinta a OB la partcula se esta deformando. Para el plano xy de la figura
la deformacion xy resulta:
C. Gherardelli

U. de Chile

35

3.1 Cinematica

1
(OA OB )
2

1 v
u
+
.
2 x y

xy =
=

La deformacion se representa mediante un tensor de deformacion, , cuya componente generica

ij esta dada por:
1
ij =
2

vj
vi
+
xi xj

Las componentes de la diagonal de este tensor representan la deformacion lineal por compresi
on
y/o traccion en los distintos ejes vista anteriormente y esta dada por:
ii =

vi
xi

Se vera mas adelante como este tensor de esfuerzos esta relacionado con los esfuerzos normales
y de corte.
Velocidad
Haciendo un desarrollo de Taylor del campo de velocidades y despreciando los termino de orden
2 y superiores se obtiene
vi (~x, t) = vi (~xo , t) +

vi
xj

xi i, j .
~
xo

El primer termino del lado derecho de la ecuacion anterior representa la traslacion por lo que el
segundo debe representar la rotacion y la deformacion. En forma matricial la ecuacion anterior
queda

u(~x)
u(~xo )
u/x u/y u/z
x

x) = v(~xo ) + v/x v/y v/z y .

v(~
w(~x)
w(~xo )
w/x w/y w/z ~x
z
o

vi
La matriz x
se puede dividir en dos matrices, una antisimetrica, que representa la rotacion,
j ~
xo
y otra simetrica, que representa la deformacion, de la siguiente manera:

u/x

1
2

u
y

u/y u/z




v
v/y v/z = 12 u
0

x 

 y

w/x w/y w/z
w
12 v
12 u
z x
z

v
x

v/x

C. Gherardelli

u
 x
1 u
2  y +
1 u
2 z +

1
2

v
x 
w
x

u
y +
v
 y
1 v
2 z +

v
x

w
y

w
y

1
2
1
2

u
 z
v
z


1
2
1
2

w
x 
w
y

0


u
+
 z
v
z +
w
z

w
x 
w
y

U. de Chile

36

3.1 Cinematica

La velocidad queda por lo tanto para un fluido de la siguiente forma

u(~x)
u(~xo )
0
z y
x
x

x) = v(~xo ) + z
0
x y +  y .
v(~
w(~x)
w(~xo )
y x
0
z
z
Recordando los resultados de la mecanica del solido, donde solo se consideran movimientos de
traslacion y rotacion, la velocidad esta dada por
~~x = V
~~x +
~ ~r ,
V
o
~ es el vector de velocidad angular. En forma matricial esta ecuacion queda
donde

u(~x)
u(~xo )
0
z y
x

x) = v(~xo ) + z
0
x y .
v(~
w(~x)
w(~xo )
y x
0
z
Vemos que este resultado es un caso particular de la ecuacion para un fluido donde no existe
deformacion.

3.1.2

Aceleraci
on

La velocidad de una partcula cualquiera sera una funcion de la posicion as como del tiempo
~ =V
~ (~r, t) .
V
La aceleracion ~a es la variacion temporal de la velocidad
~a =

d~
d~
V (~r, t) = V
(x, y, z, t) .
dt
dt

Aplicando la regla de la cadena se obtiene

~a =

~ y V
~ z V
~
~ x V
V
+
+
+
x |{z}
t
y |{z}
t
z |{z}
t
t
u

~a =

~
~
~
V
V
V
u+
v+
w
x
y
z

{z

aceleracion
convectiva

~
V
t

(3.3)

| {z }

aceleracion
local

Se ve que existen dos efectos superpuestos en la aceleracion:

Aceleracion local: representa la variacion de la velocidad de una partcula en la posici
on
ocupada por esta, es decir, representa los efectos no permanentes existentes en un flujo.
Aceleracion convectiva: representa el hecho de que una propiedad asociada a una partcula
fluida puede cambiar debido al movimiento de esta de un punto en el espacio a otro.

C. Gherardelli

U. de Chile

37

3.1 Cinematica

Las componentes escalares de esta ecuacion en coordenadas cartesinas son:



ax =


ay =

u
u
u
u+
v+
w +
x
y
z

v
v
v
u+
v+
w +
x
y
z

aw =

w
w
w
u+
v+
w +
x
y
z


u
t

v
t


w
t

,


La ecuacion 3.3 puede reescribirse de la siguiente forma:

~a =



~
V
~.
~ V
+ V
t
| {z }
operador

~ es un operador matematico que, en el caso de la ecuacion anterior, se encuentra operando

V
~ . De lo anterior se puede decir que:
sobre la velocidad V

()  ~
+ V ()
t

es tambien un operador que, para el caso de la aceleracion, opera sobre la velocidad. Este
operador se denomina derivada material, sustancial o total y se representa por
D()
.
Dt

~a =

~)
D(V
.
Dt

El concepto de derivada total es aplicable a distintos parametros del flujo y no solo a la aceleracion. Para la temperatura T (x, y, z, t), que se diferencia de la velocidad por ser un campo
escalar, por ejemplo, la derivada total resulta
T
DT
~ T
=
+V
Dt
t
y para la presion p(x, y, z, t)
p ~
Dp
=
+ V p
Dt
t
donde vemos que el operador opera primero sobre el campo escalar.

C. Gherardelli

U. de Chile

38

3.2 Conservaci
on de la masa / Ecuaci
on de Continuidad

3.2

Conservaci
on de la masa / Ecuaci
on de Continuidad

El principio de conservacion de masa establece que la variacion

temporal de la cantidad de masa contenida en un volumen mas el
flujo de masa neto que pasa a traves de las paredes que encierran
el volumen es igual a cero. Para el desarrollo del principio de
conservacion de masa se considera una partcula diferencia como
la de la figura, la cual se supone inmersa en flujo con velocidad

~ = u + v
V
+ wk.
La cantidad de masa contenida en el volumen es igual a la masa
por unidad de volumen, es decir la densidad, por el volumen del
elemento diferencial.

Lneas y tubo de corriente

m = = x y z
La variacion temporal de m es, por lo tanto,
(m)

=
x y z
t
t

Figura 3.5: Volumen de Control diferencial.

El flujo de masa m
a traves de una superficie cualquiera A esta dado por
Z

m
=

~ dA
~ [kg/s]
V

Como convencion para definir el signo del producto punto se define la normal a la superficie
que encierra al volumen positiva hacia afuera del volumen. Lo anterior implica que para un
~ dA
~ sera negativo y para un flujo que sale
flujo que entra al volumen el signo del producto V
del volumen sera positivo. El flujo neto sera, por lo tanto, lo que sale menos lo que entra al
~ dA
~ representa el producto
volumen considerado. Otro punto a mencionar es que el producto V
escalar de la componente normal de la velocidad por el diferencial de area, es decir, Vn dA. La
componente tangencial al area dA no aporta flujo masico a traves de la superficie.
Para el elemento diferencial de la figura la ecuacion anterior se debe aplicar a cada una de las
caras del elemento y se debe realizar un desarrollo de Taylor de primer orden para la variable
V . Para la direccion se tiene que el flujo neto es


( u) x
( u) x
u +
y z u
y z =
x 2
x 2

C. Gherardelli

(u)
x y z
x


U. de Chile

39

3.3 Cantidad de Movimiento Lineal

El flujo neto en todas las direcciones es la suma de las componentes de cada direccion, es decir,


(u) ( v) ( w)
+
+
x y z
x
y
z


(m)

=
x y z +
t
t

(u) ( v) ( w)
+
+
x y z = 0
x
y
z


Agrupando los terminos y dividiendo por el volumen se obtiene finalmente:

( u) ( v) ( w)
+
+
+
=0
t
x
y
z

(3.4)

Esta ecuacion es la forma diferencial del principio de conservacion de masa para un sistema de
coordenadas cartesiano. En la Mecanica de Fluidos este principio se denomina Ecuaci
on de
Continuidad. Para un sistema de coordenadas cilndrico la ecuacion de continuidad esta dada
por:
1 (r ur ) 1 ( u ) ( uz )
+
+
+
=0
t
r
r
r
z

(3.5)

En forma vectorial la ecuacion de continuidad queda de la siguiente forma:

~ =0
+ V
t

(3.6)

Desarrollando el segundo termino se tiene

~
~ =0
+ V + V
t
Se puede apreciar que los dos primeros terminos representan la derivada total de la densidad
por lo que la ecuacion anterior se puede reescribir de la siguiente forma:
D
~ =0
+ V
Dt

(3.7)

~ = 0 se puede considerar, por lo tanto, que el flujo se

Para los flujos en que se cumple que V
comporta como un flujo incompresible.

3.3

Cantidad de Movimiento Lineal

3.3.1

Tensor de esfuerzos

Sobre un elemento diferencial de fluido act

ua una distribucion de esfuerzos seg
un todas las
direcciones como se muestra en la figura 3.6. Esta distribucion de esfuerzos se agrupa en un
:
tensor denominado tensor de esfuerzos

xx xy xz

yx yy yz .
zx zy zz
C. Gherardelli

U. de Chile

40

3.3 Cantidad de Movimiento Lineal

z s
zz

tzy

tzx

txz
x

txy

sxx

tyx

tyz

syy
y

Figura 3.6: Esfuerzos sobre un elemento de fluido.

El elemento generico de este tensor es i,j donde el subndice i representa la normal al plano
asociado con la tension y el subndice j representa la direccion de la tension. Por convencion se
adopta la siguiente convencion de signos:
La normal a la superficie es positiva hacia afuera del volumen que encierra la superficie.
Una componente de la tension es positiva cuando tanto el vector que representa la superficie
sobre la que act
ua la tension como la tension misma tienen sentidos coincidentes, es decir,
ambos positivos o ambos negativos.
Se cumple ademas que i,i = i,i que representa la componente de esfuerzos normal al plano i.
Se puede demostrar que este tensor es simetrico, es decir
xy = yx
xz = zx
zy = yz

xx xy xz

xy yy yz
xz yz zz
Se puede demostrar tambien que la suma de las tensiones normales es una invariante, es decir,
no depende de los ejes del sistema coordenado (x, y, z en este caso). Se define
1

= (xx + yy + zz )
3
como la tension volumetrica que es un escalar. Para el caso de fluidos ideales o no viscosos se
tiene que xx = yy = zz que se definio la presion como el negativo de la tension normal, es
decir

= xx = p .
C. Gherardelli

U. de Chile

41

3.3 Cantidad de Movimiento Lineal

Para el caso de fluidos viscosos se define la presion termodinamica como la tension volumetrica
con signo cambiado, es decir

= p

3.3.2

Ecuaciones de movimiento / Ecuaciones de Navier-Stokes

Para determinar las ecuaciones de movimiento para un fluido viscoso se analizara el elemento
diferencial de la figura 3.7 donde se han considerado solo las fuerzas seg
un la direccion y.
tzy dz
z 2

(tzy+
(syy-

syy dy
y 2

) dy dx

(txy-

) dx dz

dz

txy dx
x 2

(syy+

) dy dz

syy dy
y 2

) dx dz

dx

(txy+

txy dx
x 2

dy

) dy dz

(tzy-

tzy dz
z 2

) dy dx

Figura 3.7: Balance de fuerzas sobre un elemento de fluido.

Realizando un balance de fuerzas y desarrollando se obtiene el siguiente sistema de ecuaciones

para todas las direcciones:

xy
yy
zy
DVy
+
+
=
x
y
z
Dt
xx xy
xz
DVx
gx +
+
+
=
x
y
z
Dt
xz
yz
zz
DVz
gz +
+
+
=
z
y
z
Dt
gy +

(3.8)

La existencia de esfuerzos de corte esta asociada a las deformaciones a que esta sometido un
elemento diferencial de fluido. Debe por lo tanto existir una relacion entre el tensor de esfuerzos
y el tensor de deformaciones  analizado en el captulo 3.1

= f () .

Fluidos de Stokes
Se define un fluido de Strokes al fluido que cumple con las siguientes condiciones:
1. Tensor de esfuerzos
es una funcion continua del tensor de velocidad de deformacion  y
del estado termodinamico local.
es independiente de la traslacion y rotacion del elemento considerado.
2.
3. las propiedades del fluido son independientes del sistema de referencia utilizado.
4. El fluido carece de elasticidad.
5. El fluido es homogeneo, la funcion f no depende explicitamente de las coordenadas.
C. Gherardelli

U. de Chile

42

3.3 Cantidad de Movimiento Lineal

6. El fluido es isotropo, es decir, las propiedades son independientes de la direccion y las

y coinciden.
direcciones principales de
Fluido Newtoniano

Se define un fluido Newtoniano como un fluido de Stokes lineal, es decir, las componentes de
son funciones lineales de las componentes de .
Bajo las condiciones anteriores la relacion que se obtiene es la siguiente:
~ I
+ pI = 2 + V

(3.9)

donde
1
p = (xx + yy + zz )
3

(viscosidad dinamica) y (segundo coeficiente de viscosidad) son constantes de proporcionalidad.

I es el tensor identidad.
De acuerdo a la ec. 3.9 el elemento generico para los esfuerzos de corte esta dado por
vi
vj
ij =
+
xj
xi

Para un sistema cartesiano de coordenadas, las componentes escalares de la anterior resultan:

Vx Vy
+
y
x

xz

Vx Vz
+
=
z
x

yz

Vz
Vy
=
+
z
y

xy =


Para un flujo con un perfil de velocidades seg

un un solo eje (x) se obtiene:
xy =

Vx
y

que representa la Ley de Viscosidad de Newton vista en el captulo 1.

Para los esfuerzos normales la ecuacion 3.9 seg
un x queda:
xx + p = 2xx +

Vx
x

Sumando las tres componentes y recordando que xx + yy + zz = 3p se obtiene

2
= .
3
C. Gherardelli

U. de Chile

43

3.3 Cantidad de Movimiento Lineal

Reemplazando el resultado anterior en la ec. 3.9 se obtiene finalmente



= 2

2
~
p + V
3

{z

(3.10)

componente normal
de la deformacion
En notacion indicial la ecuacion anterior queda:
2
~
ij = 2ij p + V
3


ij

(3.11)

~ = 0) o es uniforme (V
~ = cte) se recupera lo visto en el captulo 2:
Si el flujo esta en reposo (V
ij = pij .
~ = 0) la componente normal del esfuerzo queda:
Si el flujo es incompresible ( V
ii = 2

Vi
p.
xi

Reemplazando el tensor de esfuerzos obtenido (ec. 3.10) en el sistema de ecuaciones 3.8 y considerando que

[%]ij =
[%] ,
xj
xi
se obtiene
Vi
Vi

+ Vj
t
xj

2
~
= gi +
(2)
p + V
xj
xi
3


Como
1
ij =
2

vj
vi
+
xi xj

se obtiene finalmente
Vi
Vi

+ Vj
t
xj

"

Vj
Vi
= gi +

+
xj
xi
xj

!#

2
~

p + V
xi
3


(3.12)

Las ecuaciones anteriores (3 ecuaciones escalares) representan las ecuaciones de movimiento

general para un fluido newtoniano y se denominan ecuaciones de Navier-Stokes. La ecuacion de
continuidad,

~ ) = 0,
+ (V
t
proporciona la ecuacion faltante para cerrar el sistema de ecuaciones. En el caso mas general
deben incluirse ademas la ecuacion de estado del fluido (f (p, , T ) = 0) y la dependencia de la
viscosidad con la temperatura y la presion ( = (T, p)). Estas ecuaciones no han sido resueltas
salvo en casos muy particulares y simples.
C. Gherardelli

U. de Chile

44

3.3 Cantidad de Movimiento Lineal

3.3.3

Flujo incompresible

La ecuacion de continuidad para un flujo incompresible esta dada por la siguiente relacion
~ = 0.
V
Veremos a continuacion como se modifican las ecuaciones de Navier-Stokes bajo esta condicion.
Para esto desarrollaremos la componente x (i = x) del segundo termino del lado derecho de la
ec. 3.12:
"

Vj
Vi

+
xj
xi
xj

!#

Vx Vx
+
x
x

Vx Vy
+
y
x

Vx Vz
+
z
x



2


Vx
2 Vx 2 Vx
Vx Vy
Vz

=
+
+
+
+
+
2
y 2
z 2
x x
y
z
x
|
{z
}
~ =0
V

"

2V

x
2
x

2V

x
2
y

2V

x
2
z

= 2 Vx .

"

Vj
Vi

+
xj
xi
xj

!#

= 2 Vi .

Realizando un desarrollo analogo seg

un los otros ejes coordenados y reemplazando en las
ecuacion 3.12 se obtienen las siguientes ecuaciones escalares en coordenadas cartesianas:
u
u
u
u

+u
+v
+w
t
x
y
z


v
v
v
v

+u
+v
+w
t
x
y
z


w
w
w
w

+u
+v
+w
t
x
y
z


p
2u 2u 2u
=
+ gx +
+ 2 + 2
x
x2
y
z

p
2v
2v 2v
=
+ gy +
+
+
y
x2 y 2 z 2

p
2w 2w 2w
=
+ gz +
+
+
z
x2
y 2
z 2

(3.13)

En coordenadas cilndricas el sistema de ecuaciones es el siguiente:

ur
u ur
u2
ur
ur

+ ur
+
+ uz
t
r
r
r
z

=
"

+
r

1
1 2 ur
2 u
2 ur
[rur ] + 2

+
r r
r 2
r2
z 2


1 p
+ g
r
"
#


1
1 2 u
2 ur
2 u
+
[ru ] + 2
+ 2
+
r r r
r 2
r
z 2


uz
uz
u uz
uz
p

+ ur
+
+ uz
=
+ gz
t
r
r
z
z
#
"


1
uz
1 2 uz
2 uz
+
r
+ 2
+
(3.14)
r r
r
r 2
z 2


u
u u
ur u
u
u
+ ur
+
+
+ uz
t
r
r
r
z

p
+ gr
r

C. Gherardelli

U. de Chile

3.3 Cantidad de Movimiento Lineal

45

En forma vectorial estas ecuaciones quedan de la siguiente forma:



~
V
~ V
~

+ V
t

C. Gherardelli

~
= p + ~g + 2 V

(3.15)

U. de Chile

46

3.3 Cantidad de Movimiento Lineal

Ejemplo Entre dos placas placas paralelas, fijas e infinitas

fluye un flujo laminar, incompresible y permanente en la direccion x. Se pide determinar el campo de velocidades y presion
que se estables bajo estas condiciones. Se pide ademas que determine el caudal volumetrico, la velocidad media del flujo y el
esfuerzo de corte en la pared w . Evaluar como cambian los resultados si la placa superior se mueve con una velocidad U en
la direccion x.

C. Gherardelli

Flujo entre dos placas

paralelas infinitas

U. de Chile

3.3 Cantidad de Movimiento Lineal

47

Ejemplo Por un tubo horizontal de radio R, diametro D =

2R, fluye un flujo laminar, incompresible y permanente en la direccion z. Se pide determinar el campo de velocidades y presion
que se establece bajo estas condiciones. Se pide ademas determinar la velocidad media del flujo,
el caudal volumetrico, el esfuerzo de corte y el factor de friccion f . Evaluar como cambian los
resultados si en el tubo original se introduce un tubo de radio R1 concentrico con el original,
por lo que el fluido pasara por un espacio anular.

C. Gherardelli

U. de Chile

48

3.3 Cantidad de Movimiento Lineal

3.3.4

Flujo turbulento

En esta seccion se veran las ecuaciones de Navier-Stokes para un flujo turbulento. Un flujo
turbulento se caracteriza por un movimiento aleatorio de las partculas fluidas con un comportamiento aleatorio de las variables del flujo como la velocidad, los esfuerzos de corte, etc.. Este
tipo de flujo se representa o modela por el valor medio (A) de la variable A mas una fluctuaci
on
0
(A ). Para la velocidad por ejemplo lo anterior queda expresado por
V =V +V0
donde
1
V =
T

tZ
0 +T

V (x, y, z, t)dt .
t0

T es el mayor periodo de la mayor fluctuacion. Lo anterior se puede ver en la figura XX.

Aplicando la definicion de promedio o media a la componente fluctuante (V 0 ) se obtiene:

V0 =

1
T

tZ
0 +T

(V V )dt

t0

t +T

tZ
0
0 +T
Z
1
V dt
V )dt

t0

t0

= V V =0,
es decir la media de las fluctuaciones es igual a cero.
Se define la intensidad de la turbulencia I como

(V 0 )2
I=
=
V

1
T

t+T
R

!1/2

(V 0 )2 dt

(3.16)

Se desarrollaran las ecuaciones de Navier-Stokes para las medias temporales de la velocidad (ya
que esta medida es facilmente cuantificable) y se vera el efecto de las fluctuaciones sobre estas.
Seg
un la coordenada x la ecuacion de Navier-Stokes es:
Vx
Vx
Vx
Vx

+ Vx
+ Vy
+ Vz
t
x
y
z


p
2 Vx 2 Vx 2 Vx
=
+ gx +
+
+
x
x2
y 2
z 2

En la ecuacion anterior se debe reemplazar Vi = Vi + Vi0 . Por ejemplo, el termino Vx (Vx /x)
queda

Vx (Vx /x) =

Vx + Vx0

= Vx
C. Gherardelli

"

(Vx + Vx0 )
x

Vx
V 0
V 0
Vx
+ Vx0 x
+ Vx x + Vx0
x
x
x
x
U. de Chile

49

3.3 Cantidad de Movimiento Lineal

Realizando todos los reemplazos y tomando la media temporal sobre toda la ecuacion (se propone
hacerlo como ejercicio) y considerando un flujo permanente se obtiene:
Vx
Vx
Vx
Vx
+ Vy
+ Vz
x
y
z

p
V 0
V 0
V 0
= +gx +2 Vx Vx0 x + Vy0 x + Vz0 x
x
x
y
z

.(3.17)

La ecuacion de continuidad para un flujo turbulento queda expresada mediante la siguiente

relacion:


Vx + Vx0
x

Vy + Vy0

Vz + Vz0
z

=0

Desarrollando y promediando en el tiempo se obtiene

Vx Vy
Vz
Vx0 Vy0 Vz0
+
+
+
+
+
=0
x
y
z
x
y
z
Como el promedio de las perturbaciones es, por definicion, igual a cero la ecuacion anterior se
puede escribir en forma separada para los promedios como para las fluctuaciones de la velocidad
de la siguiente manera:
Vx Vy
Vz
+
+
=0
z
x
y
Vx0 Vy0 Vz0
+
=0
+
y
z
x
De lo anterior se puede demostrar que
V 0
V 0
V 0
Vx0 x + Vy0 x + Vz0 x
x
y
z

(Vx0 )2 (Vx0 Vy0 ) (Vx0 Vz0 )

=
+
+
z
x
y

Sustituyendo en la ecuacion 3.17 se obtiene

Vx

Vx
Vx
Vx
+ Vz
+ Vy
z
y
x

p
+gx +2 Vx
x

(Vx0 )2

(Vx0 Vy0 ) (Vx0 Vz0 )

(3.18)
+
+
z
y

Comparando la ecuacion anterior con la ecuacion 3.15 se puede ver que la existencia de fluctuaciones en la velocidad genera esfuerzos en el fluido y estos afectan la velocidad media del
flujo. Estos esfuerzos se denominan esfuerzos aparentes o de Reynolds. Considerando todas las
direcciones se obtiene un tensor de esfuerzos denominado tensor de esfuerzos aparente:
0
0
0
(Vx0 )2 (Vx0 Vy0 ) (Vx0 Vz0 )
xx
xy
xz

0
0
0
0 = yx
yy
yz

= (Vy0 Vx0 ) (Vy0 )2 (Vy0 Vz0 )

0
0
0
zx
zy
zz
(Vz0 Vx0 ) (Vz0 Vy0 ) (Vz0 )2

C. Gherardelli

U. de Chile

50

3.3 Cantidad de Movimiento Lineal

Lo anterior se puede interpretar tambien como que el esfuerzo total en un flujo turbulento se
compone de un valor medio, asociado a la viscosidad del fluido, mas una fluctuacion, asociada
naturalmente a la turbulencia existente en el flujo, es decir:
= + 0

(3.19)

Escribiendo las ecuaciones de Navier-Stokes en estos terminos en notacion indicial resulta:


DVi
1 p
1 
T
=
+ gi +
ji + ji
Dt
xi
xj

(3.20)

Un punto importante a considerar es que la existencia de fluctuaciones introduce nuevas

incognitas y por lo tanto se requiere de nuevas ecuaciones para cerrar y solucionar
(numericamente) el sistema ecuaciones. Existen diversos modelos, llamados modelos de cierre,
que proporcionan estas ecuaciones adicionales. El estudio de estos modelos queda fuera del
alcance de este curso por lo que no ser
an tratados.

3.3.5

Fluido ideal

Para un fluido ideal se cumple que = 0 por lo que las ecuaciones de Navier-Stokes se reducen
a:


~
V
~ V
~

+ V
t

= p + ~g

(3.21)

Esta ecuacion se denomina Ecuaci

on de Euler y representa la forma diferencial de la segunda
ley de movimiento de Newton para un fluido ideal, donde no existen esfuerzos de corte. Para
un flujo permanente la ecuacion de Euler se reduce a:


~ V
~ = 1 p + ~g
V

El termino del lado izquierdo de la ecuacion anterior se puede desarrollar de la siguiente forma:


2
~ V
~ = V
V
2

~ ( V
~ ).
V

~ = 0, es decir que el flujo es irrotacional, la ecuacion de Euler queda

Si se cumple que V
V2

1
= p + ~g

Realizando un producto punto entre la ecuacion anterior y un elemento diferencial de movimiento

d~r seg
un una direccion arbitraria y suponiendo que ~g = g
z , se obtiene:
dp
V2
+d

C. Gherardelli

+ gdz = 0

U. de Chile

51

3.3 Cantidad de Movimiento Lineal

Para un flujo incompresible, la ecuacion anterior se puede integrar obteniendo finalmente:

1
p + gz + V 2 = cte.
2
La ecuacion anterior representa un balance de energa y se denomina ecuacion de Bernoulli. Se
ve que la energa de un fluido esta dada por la energa asociada a la presion p, la energa potencial
gz y la energa cinetica 1/2V 2 . Como la direccion d~r del desplazamiento fue elegida en forma
arbitraria, se puede concluir que la constante de integracion es la misma entre cualquier par de
puntos del fluido. En la siguiente seccion se vera mas en profundidad este punto.

3.3.6

Din
amica elemental

En esta seccion se analizara la ecuacion de cantidad de

movimiento lineal desde un sistema coordenado (
s, n
)
coincidente con la lnea de corriente, donde s es la
posicion a lo largo de la lnea de corriente. El analisis
se hara para un regimen permanente, por lo que las
lneas de corriente estaran fijas en el espacio. La velocidad de la partcula elemental en este sistema coordenado estara dada por:
Partcula fluida sobre una lnea de
corriente
~ =V
~ (s, t) s ,
V
dado que la velocidad es tangente a la lnea de corriente. El vector unitario s es, sin embargo,
una funcion tanto de s como n, es decir, s = s(s, n). La aceleracion de la partcula esta dada
por:
~a =

~
D(V s)
DV
D
s
DV
=
=
s + V
.
Dt
Dt
Dt
Dt

Como /t = 0 la ecuacion anterior queda

~
V
DV
= V
Dt
s


s + V

s
V
s

La derivada del vector unitario s resulta

s
n

s
= lim
= ,
s s0 s
R
donde n
es el vector normal a s y R el radio de curvatura de la lnea de corriente en el punto.

~a =

V2
V
.
n

s +
s }
|R{z }
| {z
componente componente
paralela a s normal a s
V

C. Gherardelli

U. de Chile

52

3.3 Cantidad de Movimiento Lineal

El termino de la aceleracion normal a s tiene su origen en el cambio de direccion de la velocidad

de una partcula al moverse sobre una trayectoria curva. Si la trayectoria es recta, es decir
R , este termino desaparace.
Analizaremos a continuacion la ecuacion de cantidad de movimiento lineal (F = ma) seg
un una
direccion de movimiento coincidente a la lnea de corriente y seg
un una direccion de movimiento
normal a la lnea de corriente. Para determinar las fuerzas externas en ambas situaciones, es
decir seg
un s y n
, se analizara el elemento diferencial de la figura 3.8.
Ecuaci
on de movimiento seg
un s; Ecuaci
on de Bernoulli
La ecuacion de cantidad de movimiento seg
un s es


p s
p s
V
ny p +
ny sny sin = snyV
s 2
s 2
s


Figura 3.8: Balance de fuerzas sobre una partcula fluida.

sin

p
V
= V
.
s
s

Se ve que para que exista movimiento debe existir un desbalance entre las fuerzas causadas por
la presion y el peso. Analizaremos a continuacion la ecuacion anterior a lo largo de la lnea de
corriente. El diferencial de la presion es


dp =

p
ds +
s


p
dn .
n


Sobre una lnea de corriente se cumple que n = cte o dn = 0, por lo que

p
dp
=
.
s
ds
Analogamente se tiene que
V

1 dV 2
V
=
s
2 ds

y
sin =

z
dz
=
.
s
ds

C. Gherardelli

U. de Chile

3.3 Cantidad de Movimiento Lineal

53

Reemplazando en la ecuacion de movimiento se obtiene

dz dp
1 d(V 2 )

=
.
ds ds
2
ds

Eliminando ds obtenemos
1
dz dp = d(V 2 )
2
o
1
dp + d(V 2 ) + dz = 0 .
2
Integrando sobre la lnea de corriente se obtiene
Z

dp 1 2
+ V + gz = C ,

(3.22)

donde C es una constante de integracion. Las ecuaciones anteriores son validas solo sobre una
lnea de corriente. Se ve que para poder integrar el primer termino de la ecuacion anterior es
necesario conocer la relacion existente entre la densidad y la presion.
Fluido incompresible.
Si la densidad es constante se obtiene
1
p + V 2 + gz = C .
2

(3.23)

La ecuacion anterior se denomina ecuacion de Bernoulli (1778) y tiene implcitas las siguientes
hipotesis
efectos viscosos despreciable,
flujo permanente,
flujo incompresible,
aplicable solo a una lnea de corriente.
La u
ltima de estas hipotesis significa que la constante de integracion sera diferente entre una
lnea de corriente a otra. La ecuacion de Bernoulli dice que para un flujo sin roce la energa total,
que es la suma de la energa cinetica, la energa potencial y la energa de presion, se mantiene
constante. Se ve que la ecuacion de Bernoulli, escrita en esta forma, tiene unidades de presi
on.
La constante C de la ecuacion de Bernoulli se denomina presion total pT , es decir
1
pT = p + V 2 + gz .
2
Por lo tanto, la presion total se mantiene constante sobre una lnea de corriente en un flujo ideal
( = 0). Vemos ademas que la presion total esta compuesta por
1 2
V = presi
on dinamica,
2
C. Gherardelli

U. de Chile

54

3.3 Cantidad de Movimiento Lineal

p = presion estatica y
gz = presion hidroestatica.
Dividiendo por g se obtiene
p
V2
+z+
= cte .
g
2g
Se puede apreciar que la ecuacion de Bernoulli se puede escribir en terminos de longitud. El
termino de elevacion z, que esta relacionado con la energa potencial se denomina altura topografica. El termino (P/g) se denomina altura de presion y representa la altura de la columna
de lquido necesaria para producir una presion p. (V 2 /2g) se llama altura de velocidad y representa la altura vertical necesaria para que si el fluido cae libremente, adquiera la velocidad
V.
Presi
on de estancamiento: se define como la presion que se obtiene al desacelerar un flujo
isoentropicamente (s=cte; proceso ideal) hasta el reposo. De la ecuacion de Bernoulli se ve que
la presion de estancamiento es igual a la presion total. Si el proceso de desaceleracion del flujo
no es ideal la presion que se obtiene es distinta a la de estancamiento y se denomina presi
on
de estancamiento local. La presion de estancamiento se obtiene de la conversion de la energa
cinetica y potencial, o de la presion dinamica e hidroestatica, en presion estatica y sera, por lo
tanto, mayor que esta.
Fluido compresible
Si suponemos ahora que el fluido es un gas ideal podemos utilizar la ecuacion de estado de los
gases ideales para expresar la dependencia de la densidad con la presion y la temperatura. De
la ecuacion de estado se obtiene
p
=
.
RT
Reemplazando en la ecuacion 3.22 se obtiene
Z

RT

1
dp
+ gz + V 2 = C ,
p
2

de donde se ve que debemos explicitar la forma en que vara la temperatura a lo largo de la

lnea de corriente. Para un flujo isotermico, es decir T = cte., se obtiene, integrando entre dos
puntos sobre una lnea de corriente
p1
1
1 2
V1 + gz1 + RT ln
= V22 + gz2 .
2
p2
2
Para un flujo isoentropico se cumple que
p
= cte.
k
Reemplazando, integrando entre dos puntos y reoordenando se obtiene


k
k1

p1 V12
+
+ gz1 =
1
2

k
k1

p2 V22
+
+ gz2 .
2
2

Esta ecuacion es equivalente a la ecuacion para un flujo incompresible salvo por el factor (k/k1)
que multiplica la presion y por el hecho de que las densidades son distintas.
C. Gherardelli

U. de Chile

3.3 Cantidad de Movimiento Lineal

55

Ecuaci
on de cantidad de movimiento seg
un n
.
Haciendo un desarrollo analogo al realizado en el punto anterior pero ahora seg
un n
se obtiene

dz
p
V 2

=
.
dn n
R

Esta ecuacion indica que la variacion en la direccion del flujo de la partcula esta acompa
nada de
una combinacion apropiada del gradiente de presion y el peso en la direccion normal a la lnea
de corriente. Si la partcula se mueve por una trayectoria recta (R ) la presion varia en
forma hidroestatica. Si por ejemplo despreciamos la gravedad o consideramos un flujo horizontal
obtenemos

V 2
p
=
,
n
R

que nos dice que la presion aumenta si uno se aleja del centro de curvatura, dado que n
apunta
hacia adentro del centro de curvatura y el termino del lado derecho de la ecuacion es positivo.
Para un s constante se tiene que ds = 0 y por lo tanto (p/n) = dp/dn. Por lo tanto, si
multiplicamos la ecuacion anterior por dn e integramos a traves de las lineas de corriente con
ds = 0 se obtiene
Z

dp
+
p

V2
dn + gz = cte. normal a la lnea de corriente.
R

Si el flujo es incompresible se tiene ademas que

Z

p+

V2
dn + gz = C .
R

Esta ecuacion nos dice que cuando una partcula viaja sobre una lnea de corriente curva (R < )
se requiere una fuerza neta adicional con direccion hacia el centro de curvatura para vencer los
efectos centrfugos asociados a la curvatura. Esta fuerza o diferencial de fuerza adicional es
proporcionada por la presion. La presion sera, por lo tanto, mayor en la parte externa que en
la parte interna de la curvatura.

C. Gherardelli

U. de Chile

3.3 Cantidad de Movimiento Lineal

56

Ejemplo Para el flujo ideal,

incompresible y permanente
de la figura describa la
variacion de la presion entre
los puntos (1) y (2) y (3) y
(4).

C. Gherardelli

U. de Chile

3.3 Cantidad de Movimiento Lineal

57

Aplicaciones Una de la aplicaciones importantes de las ecuaciones vistas anteriormente es la posibilidad de medir la velocidad
de un flujo a traves de la medicion de diferencias de presion. Una forma sencilla de lograr esto es la que se muestra en la figura.
Se pide evaluar la velocidad de una partcula
que pasa por punto 1, V1 si la linea de corriente de dicha partcula es horizontal y pasa
por el punto 2.

C. Gherardelli

U. de Chile

3.3 Cantidad de Movimiento Lineal

58

Ley de vaciado de un estanque Determinar la velocidad del lquido a la salida del estanque suponiendo
que el nivel del estanque se
mantiene constante y la viscosidad del lquido es despreciable. Rehaga el analisis para
el caso en que el estanque baja
de nivel.

C. Gherardelli

U. de Chile

3.3 Cantidad de Movimiento Lineal

59

Ejemplo A traves del Venturi de la figura circula

Kerosene (SG=0.85).
El
rango de flujos de Kerosene
se encuentra entre 0.005 y
0.050 m3/s. Determinar la
variacion de presiones asociada.

C. Gherardelli

U. de Chile

3.3 Cantidad de Movimiento Lineal

60

Ejemplo Para el sifon de

la figura determine la altura
maxima H tal que no se produzca cavitacion si el lquido
es agua y se encuentra a 60 F
(ps=0.256 psia). Patm=14.7
psia.

C. Gherardelli

U. de Chile

3.4 Flujo Potencial

3.4

61

Flujo Potencial

Se analizara en este captulo un tipo particular de flujo o escurrimiento denominado flujo potencial. Este tipo de flujo se denomina as ya que es posible definir una funcion potencial
mediante la cual se puede representar el campo de velocidades. La condicion necesaria para
~ = 0. Si
la existencia de la funcion potencial es que el flujo sea irrotacional, es decir, V
bien la condicion de irrotacionalidad en un flujo es difcil de encontrar existen, en algunos flujos,
zonas las cuales pueden ser tratadas como si el flujo fuese irrotacional. Para que una partcula
fluida, originalmente sin rotacion, comience a rotar se requiere de un esfuerzo de corte. Como se
vio anteriormente los esfuerzos de corte estan asociados a la viscosidad y los gradientes de
velocidad en la direccion normal al desplazamiento (V /n). Para fluidos de viscosidad baja,
como el aire por ejemplo, los esfuerzos de corte estaran asociados principalmente a la existencia
de gradientes de velocidad. En las regiones del flujo donde no existan gradientes de velocidad el
flujo podra ser considerado como irrotacional. De particular interes es el estudio de flujo alrededor de cuerpos solidos inmersos en un flujo, como un perfil alar por ejemplo. Sobre la pared
del cuerpo, y por el principio de adherencia, el fluido tendra una velocidad relativa al cuerpo
nula. A medida que uno se separa del cuerpo la velocidad del fluido aumenta aproximandose a la
velocidad de la corriente libre a partir de una cierta distancia, a partir de la cual practicamente
no existen gradientes de velocidad. La zona cercana al cuerpo es una zona de grandes gradientes
de velocidad y por lo tanto una zona donde los esfuerzos de corte son importantes. Esta zona
se denomina capa lmite y sera estudiada en el captulo 7. En la zona fuera de la capa lmite los
gradientes de velocidad desaparecen y con ellos los esfuerzos de corte, por lo que el flujo puede
ser considerado como irrotacional.

Figura 3.9: Flujo irrotacional y capa lmite sobre un cuerpo.

Ademas de la condicion de irrotacional se supondra que el fluido es incompresible ( = cte), el

flujo es permanente (/t = 0), y se analizaran solamente flujos bidimensional, es decir un flujo
donde las propiedades y caractersticas del flujo son independientes de una de las coordenadas
espaciales (2D).

3.4.1

Funci
on potencial

De la condicion de irrotacionalidad de un flujo se obtiene que

Vy
Vz
=
y
z
Vz
Vx
=
z
x

C. Gherardelli

U. de Chile

3.4 Flujo Potencial

62

Vy
Vx
=
x
y
Analizando estas relaciones se ve que las componentes de la velocidad se pueden expresar mediante una funcion escalar (x, y, z) tal que
Vx =

Vy =

Vz =

La funcion se denomina funcion potencial de velocidades y se cumple que

~ =
V
Reemplazando la relacion anterior en la ecuacion de continuidad para un flujo incompresible se
obtiene
~ = = 0
V
2 = 0
que se conoce como ecuacion de Laplace. En coordenas rectangulares la ecuacion anterior queda
2 2 2
+ 2 + 2 =0
x2
y
z
Una caracterstica importante de la ecuacion de Laplace es que es una ecuacion en derivadas
parciales lineal, lo que implica que si 1 y 2 son soluciones o satisfacen la ecuacion 2 = 0,
entonces 3 = 1 + 2 sera tambien solucion de la ecuacion de Laplace. Esta caracterstica
permite generar diferentes tipos de flujos a partir de otros conocidos superponiendo las funciones
potenciales respectivas. Esto se conoce como superposicion de flujos.

3.4.2

Funci
on de corriente

~ = 0, que en
La ecuacion de continuidad para un flujo incompresible y permanente es V
coordenadas cartesianas y para un flujo bidimensional resulta
Vx Vy
+
=0
x
y
Analizando la ecuacion anterior se ve que es posible definir una funcion = (x, y), llamada
funcion de corriente, tal que
Vx =

C. Gherardelli

U. de Chile

63

3.4 Flujo Potencial

Vy =

Reemplazando en la ecuacion de continuidad. se obtiene

y
x


=0

2
2

=0
xy yx
de donde vemos que satisface la ecuacion de continuidad. Obtenemos de esta manera nuevamente una reduccion del n
umero de funciones necesarias para representar el campo de velocidades. Se ve ademas, de la ecuacion anterior, que la funcion de corriente satisface tambien la
ecuacion de Laplace

2 = 0
Las lneas para las cuales la funcion de corriente es constante son las lneas de corriente. Diferenciando se obtiene

dx +
dy = 0
x
y

d =

Vy dx + Vx dy = 0
Esta ecuaci
on representa, como se vio anteriormente, la ecuacion para las lneas de corriente.
La variacion del valor de la funcion de corriente, entre
dos lneas de corriente, esta relacionado con el caudal que pasa entre ellas. La ecuacion de continuidad
aplicada a la figura queda
dq = Vx dy Vy dx
Introduciendo la funcion de corriente

q
Vydx

y
2

Vxdy
y
1

Caudal entre lneas de corriente.

dq =

dy +
dx = d
y
x

Integrando entre 1 y 2 se obtiene

Z2

q=

Z2

dq =
1

d
1

q = 2 1
C. Gherardelli

U. de Chile

64

3.4 Flujo Potencial

Se ve que la diferencia del valor de la funcion de corriente entre dos lneas de corriente es igual
al caudal volumetrico, por unidad de profundidad, que pasa entre las dos lneas.
Para una lnea de corriente se tiene que


dy
dx

=
=cte

Vy
Vx

que representa la pendiente de las lneas de corriente. La pendiente de las lneas equipotenciales,
es decir las lneas para las cuales = cte, resulta de igualar a cero el diferencial de , es decir
d =

dx +
dy = 0
x
y

Vx dx + Vy dy = 0

dy
dx

=
=cte

Vx
Vy

Multiplicando ambas pendientes se obtiene



dy
dx

dy
dx

=cte

dy
dx

= 1
=cte

o


dy
=
dx
=cte


1
=cte

lo cual indica que la interseccion de las lneas equipotenciales y las lneas de corriente ocurre
formando un angulo recto, es decir y son perpendiculares entre si. Esta condicion se utiliza
para representar un flujo graficamente mediante una malla formada por las lneas de corriente
y las equipotenciales.

3.4.3

Circulaci
on

La circulaci
on se define como la integral de lnea, sobre
una curva cerrada, de la componente tangencial de la
velocidad a lo largo de la curva, es decir,
I

ds

~ d~s
V

Aplicando el teorema de Stokes se obtiene ademas que

Circulacion.
Z

~ ) dA
~
( V

~ = 0, entonces no existira circulacion. Como se ver

Se ve que si el flujo es irrotacional, V
a
mas adelante la circulacion posee gran importancia en la teora de la sustentacion.
C. Gherardelli

U. de Chile

65

3.4 Flujo Potencial

En coordenadas cilndricas la ecuacion de Laplace para la funcion de corriente, la ecuacion de

continuidad y las componentes de la velocidad se expresan respectivamente por las siguientes
relaciones
2 =

r
r r
r


1 2
=0
r2 2

1
1 V
(r Vr ) +
=0
r r
r
Vr =

=
r
r

V =

3.4.4

1
=
r
r

Aplicaciones

Flujos simples
Se presentaran a continuacion algunos flujos bidimensionales sencillos y sus correspondientes
funciones de corriente y potenciales.
Flujo uniforme
El flujo mas sencillo es aquel que tiene lneas de corriente rectas y paralelas y donde la magnitud
de la velocidad es constante. Este tipo de flujo se llama flujo uniforme.

U
f

f
2

f
3

y
4

y
3

y
2

y
1

f
1

f
2

f
3

y
4

y
3

y
1

x
(a) Flujo uniforme paralelo a x

x
(b) Flujo uniforme inclinado

Figura 3.10: Flujo uniforme

Si la velocidad del flujo (U ) es paralela al eje x se tendra ademas que Vx = U y Vy = 0. De las

relaciones anteriormente vistas para la funcion potencial se obtiene que

=U
x

C. Gherardelli

U. de Chile

66

3.4 Flujo Potencial

= 0.
y
Integrando se obtiene
=Ux+C,
donde C es una constante de integracion que elegimos arbitrariamente igual a cero (C = 0)

=Ux.
Se ve que las lneas equipotenciales son paralelas al eje y. La funcion de corriente correspondiente
al flujo uniforme se obtiene a partir de

=U
y
y

=0
x

=Uy
que son lneas paralelas al eje x. y se pueden apreciar en la figura 3.10(a) para un flujo
uniforme paralelo al eje x. Si el flujo forma un angulo c/r al eje x se obtienen las siguientes
funciones de corriente y potencial respectivamente (figura 3.10(b))
= U (y cos x sin ) ,
= U (x cos + y sin ) .
Fuente y sumidero
Consideraremos ahora un fluido que fluye en forma
radial a partir de un punto y en todas las direcciones.
Si q es la razon volumetrica de fluido, por unidad de
profundidad, que sale de la fuente, por conservacion
de la masa se debe cumplir que
2 r Vr = q ,

y
q
q

de donde se puede despejar Vr

Fuente/sumidero.
Vr =

q
.
2r

C. Gherardelli

U. de Chile

67

3.4 Flujo Potencial

Como el flujo es radial se cumple ademas que V = 0.

q
=
r
2r
y
1
=0
r
de donde
=

q
ln r
2

que representa la ecuacion de una familia de crculos concentricos centrados en el origen.

Si q es positivo entonces el flujo es radial hacia afuera y se denomina fuente. Si q es negativo el
flujo es radial hacia adentro y se denomina sumidero. El caudal q se denomina intensidad de la
fuente o sumidero.
La funcion de corriente se obtiene de
1
q
= Vr =
r
2r
y

=y

cte

=0

f=

cte

que representa una familia de lneas radiales.

Lneas de corriente y
equipotenciales.

V
ortice libre o irrotacional
En este tipo de flujo las lneas de corriente son crculos concentricos1 como se muestra en la
figura.
Para este caso se tiene que Vr = 0 y V = V(r) . Las
funciones potencial y de corriente que se obtiene para
este caso son
=K

Vq
q

=y

f=

cte

cte

= K ln r ,
Vortice libre.
1

Dado que el flujo esta representado por un potencial de velocidades el flujo debe ser irrotacional. Esto puede
generar confusi
on con el tipo de flujo. Debe recordarse que la rotacionalidad esta relacionada con el cambio de
orientaci
on de una partcula fluida y no con la trayectoria seguida por la partcula.

C. Gherardelli

U. de Chile

68

3.4 Flujo Potencial

donde K es una constante. La velocidad V se obtiene de

V =

1
=
r
r

V =

K
.
r

Se ve que la velocidad vara inversamente proporcional con la distancia al centro, es decir, si r

entonces V y viceversa. Ejemplos que pueden ser aproximados mediante este tipo de flujo
son el tornado y el flujo de agua saliendo por un drenaje.
Se puede demostrar que
K=

donde es la circulacion sobre una curva que encierra el origen. La circulacion sera distinta de
cero ya que el origen representa una singuralidad dentro del flujo donde V . Sobre una
curva que no encierre al origen la circulacion sera cero ( = 0). Se obtiene por lo tanto
=

y
=

ln r .
2

Doblete
La combinacion de una fuente y un sumidero, de igual
intensidad, separados por una distancia infinitesimal
origina lo que se denomina doblete. Para la fuente y
el sumidero, separados por una distancia 2a la funcion
de corriente esta dada por

y
r1

q
= (1 2 ) .
2

r
q

-q

r2
q
2

Fuente y sumidero.

Expresando la funcion anterior en funcion del angulo se obtiene

q
2ar sin
= tan1
2
r2 a2


Para valores peque

nos de a la ecuacion anterior queda
=

qar sin
.
(r2 a2 )

El doblete se obtiene haciendo tender a 0 y q de tal forma que el producto (qa/) sea
constante. Para este caso se obtiene que
r2

1
r

2
a
r

C. Gherardelli

U. de Chile

69

3.4 Flujo Potencial

y=

cte

de donde
=

k sin
r

donde k = qa/ se denomina intensidad del doblete.

El potencial de velocidades asociado al doblete resulta
=

3.4.5

k cos
.
r

Lineas de corriente.

Superposici
on de Flujos

Como se menciono anteriormente los flujos potenciales estan gobernados por la ecuacion de
Laplace. Esto significa que se pueden combinar diferentes flujo potenciales para formar otros
de interes. Otro punto que se debe recordar es que a traves de una lnea de corriente no existe
flujo por lo que puede ser considerada como una pared solida. Lo anterior indica que si se logran
combinar distintos tipos de flujo de tal manera que una lnea de corriente tenga la forma de un
cuerpo particular, se puede analizar analticamente el flujo que se establece alrededor del cuerpo.
Este metodo se denomina superposicion. A continuacion se veran algunos ejemplos simples de
superposicion.
Fuente y flujo uniforme

La funcion de corriente y la funcion potencial para la

superposicion del flujo uniforme y la fuente esta dado
por
= uniforme + fuente

Punto de
estancamiento

Flujo uniforme y fuente.

= U r sin +

q
,
2

= U r cos +

q
ln r .
2

En alg
un punto del eje x (negativo) la velocidad de la fuente se anulara con la del flujo uniforme
y se formara, por lo tanto, un punto de estancamiento. Para la fuente se tiene que Vr = q/2r
por lo que el punto de estancamiento es tal que en x = b, U = q/2 r
b=

q
.
2U

(3.24)

Evaluando para r = b y = se obtiene

estancamiento =

C. Gherardelli

q
= bU .
2

U. de Chile

70

3.4 Flujo Potencial

Graficando estos resultados vemos como esta combinacion de flujos puede ser utilizada para analizar el
flujo sobre un cuerpo inmerso en un flujo uniforme.
Para esta combinacion el cuerpo es como el que muestra la figura, el cual se encuentra abierto aguas abajo.
Con la funcion de corriente conocida se puede obtener
el campo de velocidades en cualquier parte del flujo.

punto de
estancamiento

Lneas de corriente.
Vr =

1
q
= U cos +
r
2r

y
V =

= U sin
r

de donde el cuadrado del modulo de la velocidad resulta

2

Vr2

V =

V =U

V2

U q cos
=U +
+
r
2

b
b2
1 + 2 cos + 2
r
r

q
2r

2

Conocida la velocidad es posible determinar ademas el campo de presiones, utilizando la ecuaci

on
de Bernoulli entre dos puntos cualesquiera del flujo ya que el flujo es irrotacional. Por ejemplo,
entre un punto lejano del cuerpo, o de la fuente, donde V = U y p = p0 y despreciando las
variaciones de z se obtiene
1
1
p0 + U 2 = p + V 2
2
2
de donde se puede despejar la presion p
b
b2
1
p = p0 U 2 2 cos + 2
2
r
r

Doblete y flujo uniforme

La superposicion de un flujo uniforme con un doblete genera el flujo alrededor de un cilindro.
La funcion de corriente y la funcion potencial son, respectivamente, las siguientes
= U r sin

k sin
r

= U r cos +

k cos
,
r

donde k es la intensidad del doblete. Para que el cuerpo que se genera con esta superposici
on
sea un cilindro se debe cumplir que = cte para r = a, donde a es el radio del cilindro. Sobre la
C. Gherardelli

U. de Chile

71

3.4 Flujo Potencial

superficie del cilindro, o sobre la lnea de corriente que representa el cilindro, se cumple adem
as
que Vr = 0

1
k
Vr =
= U 2
r
r


r cos = 0 ,

de donde


k
r2

=0
r=a

k = U a2 .
y = cte

Reemplazando en y en se obtiene
!

=Ur

a2
1 2
r

=Ur

a2
1+ 2
r

sin ,

cos .
Lneas de corriente para un
cilindro.

Se ve que para r = a, = 0. En los puntos de estancamiento se cumple que V = 0

1
V =
= U r
r

a2
1+ 2
r

sin = 0

de donde sin = 0 o = .
Sobre la superficie del cilindro, es decir, para r = a, se tiene que V = 2U sin de donde las
velocidades maximas se obtiene para = /2
V,max = V ( = /2) = 2 U .
La distribucion de presiones en la superficie del cilindro (ps ) se obtiene utilizando la ecuaci
on
de Bernoulli y resulta
1
ps = po + U 2 (1 4 sin2 )
2
Integrando la presion ps sobre el manto del cilindro se puede obtener tanto la fuerza horizontal
o arrastre y la fuerza vertical o sustentacion a la cual esta sometido el cilindro. Para este caso,
y dada la simetra del flujo que se genera en torno al cilindro como se puede ver de la figura,
ambas fuerzas tendran un valor cero.
V
ortice, doblete y flujo uniforme
La funcion de corriente y la funcion potencial para esta superposicion de flujos son
a2
= Ur 1 2
r
C. Gherardelli

sin

ln r
2
U. de Chile

72

3.4 Flujo Potencial

a2
= Ur 1 + 2
r

cos

respectivamente, donde es la circulacion. Se puede ver que para r = a, = cte por lo que el
cuerpo generado es, al igual que el caso anterior, un cilindro de radio a. La diferencia es que el
cilindro generado por esta superposicion se encuentra girando en el sentido de giro del vortice
libre. La velocidad tangencial sobre la superficie (V,s ) toma ahora el siguiente valor
V,s =

= 2U sin +
.
r
2a

La forma que adquiere el flujo, y por lo tanto la forma que tienen las lneas de corriente, dependen
de la intensidad del vortice. La posicion de el/los puntos de estancamiento en la superficie de
cilindro se encuentran imponiendo la condicion V = 0

sin estanc. =

4U a

En la figura 3.11 se muestran las diferentes posibilidades que se pueden presentar, de acuerdo al

valor de 4U
a . Se ve que si 4U a > 1 entonces el punto de estancamiento no se encuentra sobre
el cilindro ya que sin estanc. > 1 no tiene solucion.

G
<1
4p Ua

G=0

punto de
estancamiento

G
=1
4p Ua

G
>1
4p Ua

Figura 3.11: Lneas de corriente y puntos de estancamiento para diferentes valores de (/4U a).

La presion sobre la superficie se encuentra utilizando al ecuacion de Bernoulli y resulta

1
2 sin
2
ps = p0 + U 2 1 4 sin 2 +
2 2 2
2
aU
4 a U
C. Gherardelli

U. de Chile

3.4 Flujo Potencial

73

La fuerza por unidad de longitud que se desarrolla sobre el cilindro se obtiene integrando ps
sobre el cilindro. Dada la simetra vertical del flujo sobre el cilindro el arrastre es cero. La
sustentacion, por unidad de longitud, resulta
FS = U .
Se puede apreciar que la fuerza de sustentacion apunta, para este caso, hacia abajo y que depende
de la densidad y velocidad del flujo libre y de la circulacion alrededor del cilindro. Si = 0
entonces se tendra que FS = 0. Para un cilindro girando en el sentido de giro del reloj, la fuerza
de sustentacion apuntara hacia arriba.

C. Gherardelli

U. de Chile

3.4 Flujo Potencial

74

Ejemplo Suponga que el flujo que se genera sobre un hangar de seccion semicircular
de diametro D = 6 m y largo L = 18 m se
puede aproximar por el flujo potencial que se
genera alrededor de un cilindro con [0, ].
Durante una tormenta el viento alcanza una
velocidad U = 100km/h y la temperatura
exterior es de 5 C. Si la presion dentro del
hangar es igual a p0 = 720 mmHg, se pide
que determine la fuerza neta sobre el hangar
que trata de levantarlo de sus fundaciones.

C. Gherardelli

U. de Chile

3.4 Flujo Potencial

75

Ejemplo Se pide determinar el campo de

velocidades que se obtiene mediante la superposicion de dos fuentes de igual intensidad q separadas por una distancia 2l. Para
un plano de simetra entre las dos fuentes se
pide determinar la velocidad del flujo a traves
del plano. Dado el resultado anterior, que
situacion real se podra representar mediante
esta superposicion? Para el plano de simetra
se pide determinar el campo de velocidades y
de presion sobre el plano.

C. Gherardelli

U. de Chile

3.4 Flujo Potencial

76

Ejemplo Un flujo potencial que fluye contra una placa plana se puede describir mediante la funcion de corriente = A xy donde
A es una constante. Este tipo de flujo permite describir aceptablemente el flujo en la
vecindad de un punto de estancamiento. Superponiendo una fuente de intensidad m en
el origen O se obtiene el flujo sobre una placa
plana con una protuberancia. Determine la
relacion entre la altura h, la constante A y la
intensidad de la fuente m.

C. Gherardelli

U. de Chile

3.4 Flujo Potencial

77

Ejemplo El flujo bidimensional de un fluido ideal e incompresible ( = 1000 kg/m3 )

en la vecindad de una esquina recta se puede
describir mediante la funcion de corriente
= 2r2 sin(2), donde tiene unidades de
m2 /s, r se expresa en metros. Determine
la funcion potencial correspondiente. Si la
presion en el punto (1) de la pared es 30 kPa,
cual es la presion en el punto (2). Asuma que
el plano xy es horizontal.

C. Gherardelli

U. de Chile