TEMA 6 enero 16

FLUJO
IDEAL 2017
D. Páez
Mecánica de Fluidos II
 TEMA 6 FLUJO IDEAL D. Páez 2017

CONTENIDO

FLUJO IDEAL O FLUJO POTENCIAL ................................................................. 3

FLUJO NO VISCOSO ............................................................................................. 3
Flujo Ideal o Potencial ............................................................................................ 3
Flujo Irrotacional..................................................................................................... 3
FUNCION DE CORRIENTE (ψ)........................................................................ 9
Ejemplo 1. Campo de velocidades irrotacional ................................................. 12
FUNCION POTENCIAL (Φ) ............................................................................. 17
Ejemplo 2. Determinación de la función Potencial .......................................... 18
Ejemplo 3. Determinación del vector velocidad ............................................... 20
Red de Flujo ........................................................................................................... 23
COORDENADAS POLARES .............................................................................. 25
Ejemplo 4. Uso de coordenadas polares ............................................................ 28
ECUACION DE EULER PARA FLUJO NO VISCOSO ............................... 33
Ejemplo 5. Aplicación de la ecuación de movimiento de Euler ..................... 37
FLUJOS POTENCIALES BIDIMENSIONALES ............................................ 42
Flujo Rectilíneo Uniforme ................................................................................... 42
Fuente y Sumidero ................................................................................................ 44
Vórtice Irrotacional ó Vórtice Libre ................................................................... 46
Ejemplo 6. Superficie libre de un vórtice irrotacional ...................................... 48
Flujo alrededor de una esquina ............................................................................ 51
SUPERPOSICION DE FLUJOS POTENCIALES BIDIMENSIONALES 52
Doblete ................................................................................................................... 52
Flujo alrededor de un cilindro ............................................................................. 58
Ejemplo 6. Flujo alrededor de un cilindro ......................................................... 59
Flujo alrededor de medio cuerpo ........................................................................ 62

1
 TEMA 6 FLUJO IDEAL D. Páez 2017

Ejemplo 7. Superposición de flujos ideales........................................................ 63

EJERCICIOS PROPUESTOS ............................................................................... 65
REFERENCIAS....................................................................................................... 70

2
 TEMA 6 FLUJO IDEAL D. Páez 2017

FLUJO IDEAL O FLUJO POTENCIAL

FLUJO NO VISCOSO
Se puede en algunos casos decir que el efecto de la viscosidad es

despreciable. Los efectos de viscosidad sólo se sienten dentro de la capa límite

muy cerca del contorno. Esta suposición permite encontrar soluciones analíticas

a problemas complejos.

El flujo ideal es un flujo no viscoso e irrotacional.

Flujo Ideal o Potencial

El flujo se puede denominar ideal si el campo de flujo cumple con las

siguientes condiciones:

1.- El flujo es conservativo

2.- El efecto de la viscosidad sobre el flujo es confinado a la zona de capa

límite

3.- No hay flujo a través del contorno

4.- El flujo es irrotacional

Flujo Irrotacional
Cada partícula de fluido puede experimentar cuatro tipos de movimientos

y deformaciones (traslación, deformación lineal, rotación y deformación angular).

3
 TEMA 6 FLUJO IDEAL D. Páez 2017

Se dice que un fluido es irrotacional si no existe una rotación neta de las

partículas de fluido, su rotación neta es cero.

𝜕𝜕
� 𝑑𝑑� 𝑑𝑑
𝑦 𝜕𝜕 B’

𝑑𝑑
𝑑𝑑 A’ 𝜕𝜕
� 𝑑𝑑� 𝑑𝑑
B 𝑑𝑑 𝜕𝜕
O
𝑑𝑑
𝑑𝑑
𝑣𝑣𝑣
𝑑𝑑
O
A
𝑢𝑢𝑢

𝑥
Figura 6.1 Rotación de una partícula de fluido

Estudiemos la rotación de una partícula (la partícula de fluido puede

deformarse), de acuerdo a la figura 6.1, se tiene que para una partícula en el plano

(x,y) de dimensiones: dx, dy. En un diferencial de tiempo OA y OB rotarán hasta

OA’ y OB’.

La velocidad angular en A’ está dada por la variación del ángulo de

rotación por unidad de tiempo:

𝒅𝒅
𝝎𝑶𝑶 = 𝐥𝐥𝐥 (𝟔. 𝟏)
𝒅𝒅→𝟎 𝒅𝒅

4
 TEMA 6 FLUJO IDEAL D. Páez 2017

donde:

𝜕𝜕
𝑑𝑑� 𝑑𝑑 𝜕𝜕
�
𝑑𝑑 ≅ 𝑡𝑡𝑡(𝛼) = 𝜕𝜕 = 𝑑𝑑
𝑑𝑑 𝜕𝜕
𝝏𝝏
𝒅𝒅 ≅ 𝒅𝒅 (𝟔. 𝟐)
𝝏𝝏

Sustituimos (6.2) en (6.1):

𝜕𝜕 𝜕𝜕
𝜔𝑂𝑂 = 𝑑𝑑 =
𝜕𝜕𝑑𝑑 𝜕𝜕

𝝏𝝏
𝝎𝑶𝑶 = (𝟔. 𝟑)
𝝏𝝏

De forma similar para la velocidad angular en B’:

𝒅𝒅
𝝎𝑶𝑶 = 𝐥𝐥𝐥 (𝟔. 𝟒)
𝒅𝒅→𝟎 𝒅𝒅

𝜕𝜕
𝑑𝑑� 𝑑𝑑 𝜕𝜕
�
𝜕𝜕
𝑑𝑑 ≅ 𝑡𝑡𝑡(𝛽) = = 𝑑𝑑
𝑑𝑑 𝜕𝜕

𝝏𝝏
𝒅𝒅 ≅ 𝒅𝒅 (𝟔. 𝟓)
𝝏𝝏

Sustituimos (6.5) en (6.4):

𝝏𝝏
𝝎𝑶𝑶 = (𝟔. 𝟔)
𝝏𝝏

5
 TEMA 6 FLUJO IDEAL D. Páez 2017

La rotación del elemento de fluido alrededor del eje Z, es el promedio de

las velocidades angulares, tomando sentido positivo anti horario:

𝜔𝑂𝑂 + 𝜔𝑂𝑂 1 𝜕𝜕 𝜕𝜕
𝜔𝑍 = = � − �
2 2 𝜕𝜕 𝜕𝜕

𝟏 𝝏𝝏 𝝏𝝏
𝝎𝒁 = � − � (𝟔. 𝟕)
𝟐 𝝏𝝏 𝝏𝝏

Se puede obtener en forma análoga las velocidades de rotación alrededor

de los ejes, 𝑥, 𝑦:

𝟏 𝝏𝝏 𝝏𝝏
𝝎𝒙 = � − � (𝟔. 𝟖)
𝟐 𝝏𝝏 𝝏𝝏

𝟏 𝝏𝝏 𝝏𝝏
𝝎𝒚 = � − � (𝟔. 𝟗)
𝟐 𝝏𝝏 𝝏𝝏

El vector rotación o el rotacional se puede escribir como:

𝝎 = 𝝎𝒊̂ + 𝝎𝒋̂ + 𝝎𝒌� (𝟔. 𝟏𝟏)

En forma vectorial:

𝒊 𝒋 𝒌
𝟏 𝟏
�⃗ = � 𝝏
𝝎 = 𝛁𝒙𝑽
𝝏 𝝏
� (𝟔. 𝟏𝟏)
𝟐 𝟐 𝝏𝝏 𝝏𝝏 𝝏𝝏
𝒖 𝒗 𝒘

6
 TEMA 6 FLUJO IDEAL D. Páez 2017

Donde las componentes del vector resultan:

𝟏 𝝏𝝏 𝝏𝝏
𝝎𝒙 = � − �
𝟐 𝝏𝝏 𝝏𝝏 𝒊

𝟏 𝝏𝝏 𝝏𝝏
𝝎𝒚 = � − �
𝟐 𝝏𝝏 𝝏𝝏 𝒋

𝟏 𝝏𝝏 𝝏𝝏
𝝎𝒌 = � − �
𝟐 𝝏𝝏 𝝏𝝏 𝒌

Recordemos que el vector velocidad está dado por:

�⃗ = 𝒖𝒊̂ + 𝒗𝒋̂ + 𝒘𝒌�

𝑽 (𝟔. 𝟏𝟏)

Por lo tanto:

𝟏
𝝎= 𝛁𝒙𝒙 = 𝟎 (𝟔. 𝟏𝟏)
𝟐

La ecuación (6.13) indica que el campo de flujo representa un flujo

irrotacional.

Para un flujo bidimensional el rotacional está dado por:

𝟏 𝟏 𝝏𝝏 𝝏𝝏
𝝎= �⃗ = � − �
𝛁𝒙𝑽 (𝟔. 𝟏𝟏)
𝟐 𝟐 𝝏𝝏 𝝏𝝏 𝒌

La vorticidad se expresa como dos veces el rotacional:

𝝃 = 𝟐𝟐 (𝟔. 𝟏𝟏)

7
 TEMA 6 FLUJO IDEAL D. Páez 2017

Recordemos que para un flujo incompresible la ecuación de continuidad

(Conservación de masa) se puede expresar vectorialmente como:

𝝏𝝏 𝝏𝝏 𝝏𝝏
𝛁. �𝑽⃗ = 𝟎 = + + =𝟎 (𝟔. 𝟏𝟏)
𝝏𝝏 𝝏𝝏 𝝏𝝏

8
 TEMA 6 FLUJO IDEAL D. Páez 2017

FUNCION DE CORRIENTE (ψ)

Es una función que es constante a lo largo de una línea de corriente. Para

flujo bidimensional, la función de corriente debe satisfacer las siguientes

condiciones, que satisfacen un flujo conservativo:

𝝏𝝏
=𝒖 (𝟔. 𝟏𝟏)
𝝏𝝏

𝝏𝝏
= −𝒗 (𝟔. 𝟏𝟏)
𝝏𝝏

En la ecuación (6.18) recordemos que u y v son las componentes en las

direcciones (x,y) respectivamente, del vector velocidad.

Se puede verificar la ecuación de continuidad en forma bidimensional a

partir de la ecuación (6.16)

𝝏𝝏 𝝏𝝏
�⃗ = 𝟎 =
𝛁. 𝑽 + =𝟎 (𝟔. 𝟏𝟏)
𝝏𝝏 𝝏𝝏

Sustituyendo las condiciones de la función de corriente:

𝜕𝜕 𝜕𝜕
𝜕� � 𝜕 �− �
�⃗ = 𝜕𝜕 𝜕𝜕
∇. 𝑉 +
𝜕𝜕 𝜕𝜕

𝝏𝟐 (𝝍) 𝝏𝟐 (𝝍)
𝛁. �𝑽⃗ = − =𝟎 (𝟔. 𝟐𝟐)
𝝏𝝏𝝏 𝝏𝝏𝝏

9
 TEMA 6 FLUJO IDEAL D. Páez 2017

La función de corriente satisface la ecuación de conservación para un

fluido incompresible.

La ecuación de una línea de corriente para flujo bidimensional está dada

por la siguiente condición (recordemos que el producto vectorial del vector

velocidad por el vector distancia a lo largo de una línea de corriente es cero):

𝒅𝒅 𝒗
= (𝟔. 𝟐𝟐)
𝒅𝒅 𝒖

Representa la pendiente de una línea de corriente, y puede escribirse como:

𝒖𝒖𝒖 − 𝒗𝒗𝒗 = 𝟎 (𝟔. 𝟐𝟐)

A lo largo de una línea de corriente la función de corriente, ψ es

constante. Siendo ψ una función en (x,y) por lo tanto podemos escribir la

derivada de ψ como:

𝝏𝝏 𝝏𝝏
𝒅𝒅 = 𝒅𝒅 + 𝒅𝒅 (𝟔. 𝟐𝟐)
𝝏𝝏 𝝏𝝏

El valor de la derivada a lo largo de una línea de corriente es por lo

tanto igual a cero:

𝜕𝜕 𝜕𝜕
𝑑𝑑 = 𝑑𝑑 + 𝑑𝑑 = 0
𝜕𝜕 𝜕𝜕

10
 TEMA 6 FLUJO IDEAL D. Páez 2017

Sustituyendo las condiciones de la función de corriente, ecuaciones (6.17)

y (6.18) en la ecuación (6.23) se obtiene:

𝑑𝑑 = −𝑣𝑣𝑣 + 𝑢𝑢𝑢 = 0

La ecuación (6.23) satisface la condición de la pendiente de una línea de

corriente.

La función de corriente se puede relacionar con el cambio en flujo

volumétrico o caudal por unidad de ancho, en flujo bidimensional. En la figura

6.2 se muestra el flujo ideal entre dos líneas de corriente separadas un diferencial

de longitud, ds.

ψ +dψ
B
dq
ψ Líneas de corriente
udy
A
-vdx

Figura 6.2. Flujo entre líneas de corriente.

Sea 𝑑𝑑 el flujo volumétrico entre dos líneas de corriente, podemos escribir

en términos de las componentes de velocidad:

11
 TEMA 6 FLUJO IDEAL D. Páez 2017

𝒅𝒅 = |𝒖𝒖𝒖 − 𝒗𝒗𝒗| (𝟔. 𝟐𝟐)

Comparando (6.24) con la ecuación (6.23), se obtiene:

𝑑𝑑 = |𝑑𝑑|

El flujo volumétrico entre dos líneas de corriente se puede determinar

utilizando la siguiente expresión:

𝝍𝟐
𝒅𝒅 = � 𝒅𝒅 = 𝝍𝟐 − 𝝍𝟏 = ∆𝝍 (𝟔. 𝟐𝟐)
𝝍𝟏

Ejemplo 1. Campo de velocidades irrotacional

�⃗ = 4𝑥𝑥𝑖 + 2(𝑥 2 − 𝑦 2 )𝑗 que
Para el siguiente vector velocidad: 𝑉

representa un flujo bidimensional, incompresible y permanente, determine:

a) si el vector representa el flujo conservativo de un fluido incompresible.

b) si es un flujo irrotacional.
c) la función de corriente.
d) el flujo volumétrico entre dos puntos A(0,1) y B(1,1)

Solución:

�⃗ = 0
Para un flujo incompresible y conservativo se debe cumplir: ∇. 𝑉

Para un flujo bidimensional se tiene, según la ecuación (6.19):

12
 TEMA 6 FLUJO IDEAL D. Páez 2017

𝜕𝜕 𝜕𝜕
�⃗ =
∇. 𝑉 + =0
𝜕𝜕 𝜕𝜕

Las componentes del vector velocidad son:

u = 4xy

v = 2(𝑥 2 − 𝑦 2 )

y las derivadas correspondientes resultan:

𝜕𝜕
= 4𝑦
𝜕𝜕

𝜕𝜕
= −4𝑦
𝜕𝜕

Sustituimos en la ecuación de conservación, ecuación (6.19):

𝜕𝜕 𝜕𝜕
�⃗ =
∇. 𝑉 + = 4𝑦 − 4𝑦 = 0
𝜕𝜕 𝜕𝜕

Se verifica que:

�⃗ = 0
∇. 𝑉

por lo tanto el flujo es conservativo.

Para un flujo irrotacional bidimensional se debe cumplir que el rotacional

sea nulo, por lo que la ecuación (6.14) se debe verificar:

1 1 𝜕𝜕 𝜕𝜕
𝜔 = ∇𝑥𝑥 = � − � = 0𝑘
2 2 𝜕𝜕 𝜕𝜕 𝑘

13
 TEMA 6 FLUJO IDEAL D. Páez 2017

Se tiene las derivadas:

𝜕𝜕
= 4𝑥
𝜕𝜕

𝜕𝜕
= 4𝑥
𝜕𝜕

Al sustituir en la ecuación (6.14) se obtiene:

1
𝜔 = (4𝑥 − 4𝑥)𝑘 = 0𝑘
2

por lo tanto el flujo es irrotacional.

La función de corriente ψ es una función tal que satisfaga las siguientes

condiciones, ecuaciones (6.17) y (6.18):

𝜕𝜕
=𝑢
𝜕𝜕
𝜕𝜕
= −𝑣
𝜕𝜕
𝜕𝜕
= 𝑢 = 4𝑥𝑥
𝜕𝜕
𝜕𝜕
= −𝑣 − 2(𝑥 2 − 𝑦 2 )
𝜕𝜕
𝜕𝜕
Podemos tomar e integrar con respecto a y:
𝜕𝜕

𝜕𝜕
= 4𝑥𝑥
𝜕𝜕

14
 TEMA 6 FLUJO IDEAL D. Páez 2017

𝜓 = 2𝑥𝑦 2 + 𝑓(𝑥)

Derivamos con respecto a x e igualamos a la condición dada (6.18):

𝜕𝜕
= 2𝑦 2 + 𝑓´(𝑥) = −2(𝑥 2 − 𝑦 2 )
𝜕𝜕

Se obtiene:

𝑓´(𝑥) = −2𝑥 2

Integrando la expresión f’(x):

2
𝑓(𝑥) = − 𝑥 3 + 𝐶
3
Finalmente la función de corriente será:

𝟐
𝝍 = 𝟐𝟐𝒚𝟐 − 𝒙𝟑 + 𝑪
𝟑

donde C = constante de integración que se determina para la función de

corriente, ψ, conocida en un punto (x,y).

El flujo entre los puntos A y B, estará dado por:

𝑞 = 𝜓𝐵 − 𝜓𝐴

2
𝜓𝐴 = 2(0)(1)2 − 03 + 𝐶
3

2
𝜓𝐵 = 2(1)(1)2 − 13 + 𝐶
3

15
 TEMA 6 FLUJO IDEAL D. Páez 2017

2 2
𝑞 = 2(1)(1)2 − 13 + 𝐶 − �2(0)(1)2 − 03 + 𝐶�
3 3

2 4
𝑞 = �2(1)(1)2 − 13 � =
3 3

𝟒
𝒒 = 𝒎𝟑 /𝒔/𝒎
𝟑

16
 TEMA 6 FLUJO IDEAL D. Páez 2017

FUNCION POTENCIAL (Φ)

Si se tiene un flujo irrotacional entonces existe una función tal que:

𝝏𝝏
=𝒖 (𝟔. 𝟐𝟐)
𝝏𝝏

𝝏𝝏
=𝒗 (𝟔. 𝟐𝟐)
𝝏𝝏

La función potencial es constante a lo largo de las líneas equipotenciales,

que serán perpendiculares a las líneas de corriente.

𝝏𝝏 𝝏𝝏
𝒅𝒅 = 𝒅𝒅 + 𝒅𝒅 (𝟔. 𝟐𝟐)
𝝏𝝏 𝝏𝝏

𝒅𝒅 = 𝒖𝒖𝒖 + 𝒗𝒗𝒗 = 𝟎 (𝟔. 𝟐𝟐)

de donde se tiene la pendiente de las líneas equipotenciales es:

𝒅𝒅 𝒖
=− (𝟔. 𝟑𝟑)
𝒅𝒅 𝒗

Recordemos que la pendiente de la línea de corriente está dada por la

ecuación (6.21):

𝑑𝑑 𝑣
=
𝑑𝑑 𝜓 𝑢

17
 TEMA 6 FLUJO IDEAL D. Páez 2017

Comparando con la pendiente de las líneas equipotenciales, obtenemos

que las líneas de corriente y las líneas equipotenciales son perpendiculares:

𝒅𝒅 𝒖 𝟏
=− =− (𝟔. 𝟑𝟑)
𝒅𝒅𝝓 𝒗 𝒅𝒅
𝒅𝒅𝝍

La función potencial satisface la ecuación de conservación de masa (para

un flujo incompresible):

𝜕𝜕 𝜕𝜕
�⃗ = 0 =
∇. 𝑉 + =0
𝜕𝜕 𝜕𝜕

𝜕 𝜕𝜕 𝜕 𝜕𝜕
� �+ � �=0
𝜕𝜕 𝜕𝜕 𝜕𝜕 𝜕𝜕

Ecuación de Laplace:

𝟐 (𝝓)
𝝏𝟐 𝝓 𝝏𝟐 𝝓
𝛁 = + =𝟎 (𝟔. 𝟑𝟑)
𝝏𝒙𝟐 𝝏𝒚𝟐

Ejemplo 2. Determinación de la función Potencial

Para el ejemplo anterior determine la función potencial.

Solución:

Como el flujo es irrotacional entonces la función potencial existe:

𝑉 = 4𝑥𝑥𝑖 + 2(𝑥 2 − 𝑦 2 )𝑗

18
 TEMA 6 FLUJO IDEAL D. Páez 2017

u = 4x

𝜕𝜕
= 𝑢 = 4𝑥𝑥
𝜕𝜕

v = 2(𝑥 2 − 𝑦 2 )

𝜕𝜕
= 𝑣 = 2(𝑥 2 − 𝑦 2 )
𝜕𝜕

𝜕𝜕
podemos integrar ,
𝜕𝜕

𝜕𝜕
= 𝑢 = 4𝑥𝑥
𝜕𝜕

𝜙 = 2𝑥 2 𝑦 + 𝑓(𝑦)

derivamos con respecto a “y”

𝜕𝜕
= 2𝑥 2 + 𝑓´(𝑦) = 𝑣 = 2(𝑥 2 − 𝑦 2 )
𝜕𝜕

𝜕𝜕
= 2𝑥 2 + 𝑓´(𝑦) = 2(𝑥 2 − 𝑦 2 )
𝜕𝜕

𝑓´(𝑦) = −2𝑦 2 ahora integramos con respecto a y

2
𝑓(𝑦) = − 𝑦 3 + 𝐶
3

Finalmente la función potencial será:

19
 TEMA 6 FLUJO IDEAL D. Páez 2017

𝟐
𝝓 = 𝟐𝒙𝟐 𝒚 − 𝒚𝟑 + 𝑪)
𝟑

Ejemplo 3. Determinación del vector velocidad

La función de corriente de cierto flujo está dada por:

𝜓 = 2𝑥𝑥.

a) Hallar el vector velocidad.

b) Verificar la ecuación de conservación.
c) Determinar si el flujo es irrotacional
d) Determinar la función potencial, si existe.
e) Dibujar la línea de corriente que pasa por (3,4)

Solución:

𝜕𝜕 𝜕𝜕
=𝑢 , = −𝑣 para 𝜓 = 2𝑥𝑥, tenemos entonces:
𝜕𝜕 𝜕𝜕

𝜕𝜕 𝜕(2𝑥𝑥)
= = 2𝑥 = 𝑢
𝜕𝜕 𝜕𝜕

𝜕𝜕 𝜕(2𝑥𝑥)
= = 2𝑦 = −𝑣
𝜕𝜕 𝜕𝜕

El vector velocidad resultante: 𝑉 = 2𝑥𝑖 − 2𝑦𝑗

Se verifica la ecuación de conservación:

𝜕𝜕 𝜕𝜕 𝜕(2𝑥) 𝜕(−2𝑦)
�⃗ =
∇. 𝑉 + = + =2−2=0
𝜕𝜕 𝜕𝜕 𝜕𝜕 𝜕𝜕

20
 TEMA 6 FLUJO IDEAL D. Páez 2017

El rotacional está dado por:

1 𝜕𝜕 𝜕𝜕 1
𝜔 = � − � = (0 − 0) = 0
2 𝜕𝜕 𝜕𝜕 𝑘 2

El rotacional es cero por lo tanto el flujo es irrotacional.

Ya que el flujo es irrotacional podemos determinar la función potencial

𝜕𝜕
u = 2x , = 𝑢 = 2𝑥
𝜕𝜕

𝜕𝜕
v = −2𝑦 , = 𝑣 = −2𝑦
𝜕𝜕

𝜕𝜕
Integramos = −2𝑦
𝜕𝜕

𝜙 = −𝑦 2 + 𝑓(𝑥) ,

derivamos con respecto a x:

𝜕𝜕
= 𝑓´(𝑥) = 2𝑥, de donde se obtiene: 𝑓(𝑥) = 𝑥 2 + 𝐶
𝜕𝜕

𝜙 = −𝑦 2 + 𝑥 2 + 𝐶 Función Potencial.

Las líneas de corriente se definen como:

𝑑𝑑 𝑣 −2𝑦 𝑦
= = =−
𝑑𝑑 𝑢 2𝑥 𝑥

21
 TEMA 6 FLUJO IDEAL D. Páez 2017

𝑑𝑑 𝑦
=−
𝑑𝑑 𝑥

𝑑𝑑 𝑑𝑑
=−
𝑦 𝑥

C
ln(𝑦) = ln(𝑥 ) + ln(𝑐) = ln( )
x

Resolvemos y se obtiene la ecuación de las líneas de corriente:

𝐱𝐱 = 𝐂

Para el punto (3,4):

3∗4=C

12 = C

Sustituimos en la ecuación general de las líneas de corriente:

xy = 12

𝐲 = 𝟏𝟏/𝐱

Gráficamente para los reales positivos, obtenemos el siguiente gráfico que

representa la línea de corriente que pasa por el punto (3,4). Ver figura 6.3. La

ecuación no está definida en x =0, donde tiene una asíntota vertical.

22
 TEMA 6 FLUJO IDEAL D. Páez 2017

35
y
30

25

20

15

10

0
0 5 10 15
x

Figura 6.3. Líneas de corriente. Ejemplo 3.

Red de Flujo
Se denomina red de flujo al conjunto de líneas de corriente y de líneas

equipotenciales de un campo de flujo, tal que la separación entre líneas contiguas

sea la misma. En la figura 6.4 se muestra la red de flujo en el caso de una

compuerta de admisión inferior, en el caso de una presa de tierra y en el caso de

una presa impermeable.

23
 TEMA 6 FLUJO IDEAL D. Páez 2017

Φ1
Ψ2 Φ2

Ψ1

Figura 6.4. Redes de Flujo. Líneas punteadas: líneas equipotenciales,

líneas continuas: líneas de corriente. Tomado de:
http://www.ingenierocivilinfo.com/2010/09/construccion-de-la-red-de-flujo.html.

24
 TEMA 6 FLUJO IDEAL D. Páez 2017

COORDENADAS POLARES
Los flujos ideales pueden representarse en forma más sencilla en

coordenadas polares. Veamos cómo se representan las coordenadas polares en

términos de las coordenadas cartesianas, en la figura 6.5 se muestra la

representación de un punto tanto en el sistema cartesiano de coordenadas como

en el sistema polar de coordenadas.

y
P

x x
θ

Figura 6.5. Coordenadas Polares y Cartesianas.

Las coordenadas cartesianas se pueden hallar en términos de las

coordenadas polares (r, θ), según la figura 6.5:

𝒚 = 𝒓𝒓𝒓𝒓(𝜽) (𝟔. 𝟑𝟑)

𝒙 = 𝒓𝒓𝒓𝒓(𝜽) (𝟔. 𝟑𝟑)

En forma similar las coordenadas polares se pueden expresar en términos

de las coordenadas cartesianas (x,y)

𝑟2 = 𝑥2 + 𝑦2

25
 TEMA 6 FLUJO IDEAL D. Páez 2017

𝒓 = �𝒙𝟐 + 𝒚𝟐 (𝟔. 𝟑𝟑)

𝒚
𝜽 = 𝒂𝒂𝒂𝒂 � � (𝟔. 𝟑𝟑)
𝒙

El vector velocidad también puede expresarse en componentes polares,

según la figura 6.6:

�𝑽⃗ = 𝒗𝒓 + 𝒗𝜽 (𝟔. 𝟑𝟑)

𝒓 𝒏

Donde:

𝑣𝑟 = 𝑐𝑐𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚 𝑟𝑟𝑟𝑟𝑟𝑟 𝑑𝑑𝑑 𝑣𝑣𝑣𝑣𝑣𝑣 𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣

𝑣𝜃 = 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡 𝑑𝑑𝑑 𝑣𝑣𝑣𝑣𝑣𝑣 𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣
𝑽 = 𝒖𝒊 + 𝒗𝒋 (𝟔. 𝟑𝟑)

𝑢 = 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 ℎ𝑜𝑜𝑜𝑜𝑜𝑜𝑜𝑜𝑜 𝑑𝑑𝑑 𝑣𝑣𝑣𝑣𝑣𝑣 𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣

𝑣 = 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣 𝑑𝑑𝑑 𝑣𝑣𝑣𝑣𝑣𝑣 𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣

θ
y vθ

vr
V

v
θ
r u
θ x

Figura 6.6. Vector velocidad en componentes polares y cartesianas.

26
 TEMA 6 FLUJO IDEAL D. Páez 2017

La ecuación de conservación para un flujo incompresible:

𝟏 𝝏(𝒓𝒗𝒓 ) 𝟏 𝝏𝒗𝜽 𝝏𝝏
𝛁. �𝑽⃗ = + + (𝟔. 𝟑𝟑)
𝒓 𝝏𝝏 𝒓 𝝏𝝏 𝝏𝝏

Función de Corriente:

𝟏 𝝏𝝏
𝒗𝒓 = (𝟔. 𝟒𝟒)
𝒓 𝝏𝝏

𝝏𝝏
𝒗𝜽 = − (𝟔. 𝟒𝟒)
𝝏𝝏

Rotacional:

𝟏 𝝏(𝒓𝒗𝜽 ) 𝝏𝒗𝒓
𝝎𝒛 = � − � (𝟔. 𝟒𝟒)
𝟐𝟐 𝝏𝝏 𝝏𝝏

Función Potencial:

𝝏𝝏
𝒗𝒓 = (𝟔. 𝟒𝟒)
𝝏𝝏

𝟏 𝝏𝝏
𝒗𝜽 = (𝟔. 𝟒𝟒)
𝒓 𝝏𝝏

Ecuación de Laplace:

𝟐
𝟏 𝝏 𝝏𝝏 𝟏 𝝏𝟐 𝝓 𝝏𝟐 𝝓
𝛁 𝝓= �𝒓 �+ 𝟐 𝟐 + 𝟐 (𝟔. 𝟒𝟒)
𝒓 𝝏𝝏 𝝏𝝏 𝒓 𝝏𝜽 𝝏𝒛

27
 TEMA 6 FLUJO IDEAL D. Páez 2017

Ejemplo 4. Uso de coordenadas polares

Se tiene un flujo rectilíneo representado por la siguiente función potencial:

𝜙 = −𝐴𝐴

A = constante.

a) Determine las componentes del vector velocidad en coordenadas

polares y en coordenadas cartesianas.
b) Halle el rotacional en coordenadas polares.
c) Halle la función de corriente.

Solución:

a.1) A partir de la función potencial, 𝜙 ,podemos hallar el vector

velocidad en coordenadas cartesianas:

𝜙 = −𝐴𝐴

La función potencial está definida por las condiciones:

𝝏𝝏
=𝒖 (𝟔. 𝟐𝟐)
𝝏𝝏

𝝏𝝏
=𝒗 (𝟔. 𝟐𝟐)
𝝏𝝏

28
 TEMA 6 FLUJO IDEAL D. Páez 2017

La componente en la dirección x del vector velocidad, u, se obtiene

tomando la derivada de la función potencial con respecto a x, ecuación (6.26):

𝜕𝜕 𝜕(−𝐴𝐴)
𝑢= = = −𝐴
𝜕𝜕 𝜕𝜕

La componente del vector velocidad en la dirección 𝑦, 𝑣, de acuerdo a la

ecuación (6.27), será:

𝜕𝜕 𝜕(−𝐴𝐴)
𝑣= = =0
𝜕𝜕 𝜕𝜕

Por lo tanto el vector velocidad resultante es:

�⃗ = 𝐀𝐢 − 𝟎𝐲
𝐕

De donde obtenemos que A = velocidad horizontal y es constante.

a.2) El vector velocidad en coordenadas polares resultará de usar las

condiciones de la función potencial en coordenadas polares:

𝝏𝝏
𝒗𝒓 = (𝟔. 𝟒𝟒)
𝝏𝝏

𝟏 𝝏𝝏
𝒗𝜽 = (𝟔. 𝟒𝟒)
𝒓 𝝏𝝏

Escribimos la función potencial en coordenadas polares, sustituimos la

coordenada 𝑥 = 𝑟𝑟𝑟𝑟(𝜃)

ϕ = −Ax = −Arcos(θ)

29
 TEMA 6 FLUJO IDEAL D. Páez 2017

La componente radial del vector velocidad será, de acuerdo a la ecuación

(6.43):

𝜕𝜕 𝜕(−𝐴𝐴𝐴𝐴𝐴(𝜃))
𝑣𝑟 = = = −𝐴𝐴𝐴𝐴(𝜃)
𝜕𝜕 𝜕𝜕

La componente normal del vector velocidad será, de acuerdo a la ecuación

(6.44):

1 𝜕𝜕 1 𝜕(−𝐴𝐴𝐴𝐴𝐴(𝜃)) 1
𝑣𝜃 = = = [𝐴𝐴𝐴𝐴𝐴(𝜃)] = 𝐴𝐴𝐴𝐴(𝜃)
𝑟 𝜕𝜕 𝑟 𝜕𝜕 𝑟

Por lo tanto el vector velocidad en componentes polares estará dado por la

siguiente expresión:

�𝑽⃗ = −𝑨𝑨𝑨𝑨(𝜽)𝒓 + 𝑨𝑨𝑨𝑨(𝜽)𝜽

b) El rotacional en coordenadas polares, está dado por la ecuación (6.42):

1 𝜕(𝑟𝑣𝜃 ) 𝜕𝑣𝑟
𝜔𝑧 = � − �
2𝑟 𝜕𝜕 𝜕𝜕

𝜕(𝑟𝑣𝜃 ) 𝜕(𝑟𝑟𝑟𝑟𝑟(𝜃))
= = 𝐴𝐴𝐴𝐴(𝜃)
𝜕𝜕 𝜕𝜕

𝜕(𝑣𝑟 ) 𝜕(−𝐴𝐴𝐴𝐴(𝜃))
= = 𝐴𝐴𝐴𝐴(𝜃)
𝜕𝜕 𝜕𝜕
1
𝜔𝑧 = (𝐴𝐴𝐴𝐴(𝜃) − 𝐴𝐴𝐴𝐴(𝜃) = 0
2𝑟

Se tiene que el rotacional es nulo, por lo que se verifica que el flujo es

irrotacional.

30
 TEMA 6 FLUJO IDEAL D. Páez 2017

c) La función de corriente estará dada por las siguientes condiciones en

coordenadas polares:

𝟏 𝝏𝝏
𝒗𝒓 = (𝟔. 𝟒𝟒)
𝒓 𝝏𝝏

𝝏𝝏
𝒗𝜽 = − (𝟔. 𝟒𝟒)
𝝏𝝏

Para determinar la función de corriente se puede resolver integrando con

respecto a la coordenada θ y luego derivando con respecto a r, o viceversa. En

este caso integramos con respecto a θ, la componente radial según la ecuación

(6.40) está dada por:

1 𝜕𝜕
𝑣𝑟 = = −𝐴𝐴𝐴𝐴(𝜃)
𝑟 𝜕𝜕

𝜕𝜕
= −𝐴𝐴𝐴𝐴𝐴(𝜃)
𝜕𝜕
𝜕𝜕
Integramos con respecto a θ y se obtiene una expresión de la función
𝜕𝜕

de corriente:

𝝍 = −𝑨𝑨𝑨𝑨𝑨(𝜽) + 𝒇(𝒓) 𝑰

Para hallar el valor de la función de r, f(r), derivamos con respecto a r:

𝜕𝜕
= −𝐴𝐴𝐴𝐴(𝜃) + 𝑓`(𝑟)
𝜕𝜕

E igualamos a la condición dada por la ecuación (6.41):

31
 TEMA 6 FLUJO IDEAL D. Páez 2017

𝜕𝜕
𝑣𝜃 = −
𝜕𝜕

𝜕𝜕
= −𝐴𝐴𝐴𝐴(𝜃) + 𝑓`(𝑟) = −𝐴𝐴𝐴𝐴(𝜃)
𝜕𝜕

𝑓`(𝑟) = 0

Integramos con respecto a la dirección radial, r, se obtiene:

𝑓(𝑟) = 𝐶 = 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐

Sustituimos en la expresión I de la función de corriente, y se obtiene la

expresión definitiva para la función de corriente. La constante C puede hallarse

conociendo una condición de borde.

𝝍 = −𝑨𝑨𝑨𝑨𝑨(𝜽) + 𝑪

32
 TEMA 6 FLUJO IDEAL D. Páez 2017

ECUACION DE EULER PARA FLUJO NO VISCOSO

Para un elemento diferencial la ecuación de movimiento se puede escribir

en forma vectorial como:

�⃗
𝑫𝑽 𝟏
�𝒂⃗ = =− 𝛁�⃗𝒑 + 𝒈𝛁
�⃗𝒉 (𝟔. 𝟒𝟒)
𝑫𝑫 𝝆

a = aceleración en un punto
V= velocidad en un punto
p = presión en un punto
h= altura de un punto

El vector aceleración es la derivada material del vector velocidad, en la

dirección x , y, z se puede escribir como:

𝝏𝝏 𝝏𝝏 𝝏𝝏 𝝏𝝏 𝟏 𝛛𝛛 𝛛𝛛
𝒂𝒙 = +𝒖 +𝒗 +𝒘 =− −𝒈 (𝟔. 𝟒𝟒)
𝝏𝝏 𝝏𝝏 𝝏𝝏 𝝏𝝏 𝝆 𝛛𝛛 𝛛𝛛

𝝏𝝏 𝝏𝝏 𝝏𝝏 𝝏𝝏 𝟏 𝛛𝛛 𝛛𝛛
𝒂𝒚 = +𝒖 +𝒗 +𝒘 =− −𝒈 (𝟔. 𝟒𝟒)
𝝏𝝏 𝝏𝝏 𝝏𝝏 𝝏𝝏 𝝆 𝛛𝛛 𝛛𝛛

𝝏𝝏 𝝏𝝏 𝝏𝝏 𝝏𝝏 𝟏 𝛛𝛛 𝛛𝛛
𝒂𝒛 = +𝒖 +𝒗 +𝒘 =− −𝒈 (𝟔. 𝟒𝟒)
𝝏𝝏 𝝏𝝏 𝝏𝝏 𝝏𝝏 𝝆 𝛛𝐳 𝛛𝛛

Recordemos que el vector velocidad en coordenadas cartesianas se puede

escribir como:

�𝑽⃗ = 𝒖𝒊 + 𝒗𝒋 + 𝒘𝒌

33
 TEMA 6 FLUJO IDEAL D. Páez 2017

Para flujo permanente se tiene que las variaciones con respecto al tiempo

son cero.

Para flujo irrotacional las condiciones para el rotacional nulo son:

1 𝜕𝜕 𝜕𝜕
𝜔𝑥 = � − � =0
2 𝜕𝜕 𝜕𝜕 𝑖
1 𝜕𝜕 𝜕𝜕
𝜔𝑦 = � − � =0
2 𝜕𝜕 𝜕𝜕 𝑗
1 𝜕𝜕 𝜕𝜕
𝜔𝑘 = � − � = 0
2 𝜕𝜕 𝜕𝜕 𝑘
De donde se tienen las siguientes igualdades:

𝝏𝝏 𝝏𝝏
= (𝟔. 𝟓𝟓)
𝝏𝝏 𝝏𝝏
𝝏𝝏 𝝏𝝏
= (𝟔. 𝟓𝟓)
𝝏𝝏 𝝏𝝏
𝝏𝝏 𝝏𝝏
= (𝟔. 𝟓𝟓)
𝝏𝝏 𝝏𝝏
Se sustituyen las ecuaciones (6.50) al (6.52) en las ecuaciones (6.47) al

(6.49), más la condición de flujo permanente, se tiene:

𝜕𝜕 𝜕𝜕 𝜕𝜕 1 ∂p ∂h
𝑎𝑥 = 0 + 𝑢 +𝑣 +𝑤 =− −𝑔 (6.47)
𝜕𝜕 𝜕𝜕 𝜕𝜕 𝜌 ∂x ∂x

𝜕𝜕 𝜕𝜕 𝜕𝜕 1 ∂p ∂h
𝑎𝑦 = 0 + 𝑢 +𝑣 +𝑤 =− −𝑔 (6.48)
𝜕𝜕 𝜕𝜕 𝜕𝜕 𝜌 ∂y ∂y

34
 TEMA 6 FLUJO IDEAL D. Páez 2017

𝜕𝜕 𝜕𝜕 𝜕𝜕 1 ∂p ∂h
𝑎𝑧 = 0 + 𝑢 +𝑣 +𝑤 =− −𝑔 (6.49)
𝜕𝜕 𝜕𝜕 𝜕𝜕 𝜌 ∂z ∂z

Se sustituye, por ejemplo en ax:

𝜕𝜕 1 𝜕𝑢2
𝑢 =
𝜕𝜕 2 𝜕𝜕
𝜕𝜕 1 𝜕𝑣 2
𝑣 =
𝜕𝜕 2 𝜕𝜕
𝜕𝜕 1 𝜕𝑤 2
𝑤 =
𝜕𝜕 2 𝜕𝜕
Se obtiene:

1 𝜕𝑢2 1 𝜕𝑣 2 1 𝜕𝑤 2 1 ∂p ∂h
𝑎𝑥 = + + =− −𝑔
2 𝜕𝜕 2 𝜕𝜕 2 𝜕𝜕 𝜌 ∂x ∂x
Agrupando se puede escribir:

1 𝜕(𝑢2 + 𝑣 2 + 𝑤 2 ) 1 ∂p ∂h
𝑎𝑥 = =− −𝑔
2 𝜕𝜕 𝜌 ∂x ∂x
Se tiene:

𝑉 = (𝑢2 + 𝑣 2 + 𝑤 2 )

representa la magnitud del vector velocidad en el punto (x,y,z), sustituyendo en la

aceleración en x:

1 𝜕(𝑉 2 ) 1 ∂p ∂h
𝑎𝑥 = =− −𝑔
2 𝜕𝜕 𝜌 ∂x ∂x

35
 TEMA 6 FLUJO IDEAL D. Páez 2017

Agrupando los términos y haciendo 𝑎𝑥 = 0:

1 𝜕(𝑉 2 ) 1 ∂p ∂h
+ +𝑔 =0
2 𝜕𝜕 𝜌 ∂x ∂x
1 p
𝜕 � 𝑉 2 + + 𝑔h�
2 𝜌
=0
𝜕𝜕
1 p
𝜕 � 𝑉 2 + + 𝑔h�
2 𝜌
=0
𝜕𝜕
Por lo tanto se tiene que en x, la cantidad:
1 p
� 𝑉 2 + + 𝑔h� = 𝑐𝑐𝑐 𝑒𝑒 𝑙𝑙 𝑑𝑑𝑑𝑑𝑑𝑑𝑑ó𝑛 𝑥
2 𝜌
El mismo análisis se puede hacer para las componentes del vector
aceleración ay y az obteniéndose la ecuación de Euler para flujo ideal:

𝑽𝟐 𝒑
+ + 𝒈𝒈 = 𝒄𝒄𝒄𝒄𝒄𝒄𝒄𝒄𝒄 (𝟔. 𝟓𝟓)
𝟐 𝝆

La ecuación de Bernoulli para cualquier punto del campo de flujo! Se

puede aplicar la ecuación de Euler entre cualquier punto del campo de flujo

siempre cuando el flujo sea no viscoso, incompresible e irrotacional.

Para flujos bidimensionales en el plano, podemos escribir la ecuación

como:

𝑽𝟐 𝒑 𝒅
+ = 𝒄𝒄𝒄𝒄𝒄𝒄𝒄𝒄𝒄 (𝟔. 𝟓𝟓)
𝟐 𝝆

36
 TEMA 6 FLUJO IDEAL D. Páez 2017

donde pd = presión dinámica = p+γh

Ejemplo 5. Aplicación de la ecuación de movimiento de Euler

El flujo en la vecindad de una esquina de 90° está descrito por la función

de corriente: 𝜓 = 2𝑟 2 𝑠𝑠𝑠(2𝜃) 𝜓(m2/s), r(m)

a) Determine la función potencial si existe.

b) Determine la presión en el punto B, si la presión en A es de 30KPa, use
una densidad de 1000kg/m3.
y

ψ2

B ψ1
r
0.5m ψ0
θ A x

1m
Figura 6.7. Ejemplo 5.
Solución:

a) Determinemos el vector velocidad:

1 𝜕𝜕 𝜕𝜕
𝑣𝑟 = , 𝑣𝜃 = −
𝑟 𝜕𝜕 𝜕𝜕

1 𝜕�2𝑟 2 𝑠𝑠𝑠(2𝜃)�
𝑣𝑟 = = 4𝑟𝑟𝑟𝑟(2𝜃)
𝑟 𝜕𝜕

37
 TEMA 6 FLUJO IDEAL D. Páez 2017

𝜕�2𝑟 2 𝑠𝑠𝑠(2𝜃)�
𝑣𝜃 = − = −4𝑟𝑟𝑟𝑟(2𝜃)
𝜕𝜕

𝑽 = 4𝑟𝑟𝑟𝑟(2𝜃)𝑟 − 4𝑟𝑟𝑟𝑟(2𝜃)𝜃

La función potencial existe si el flujo es irrotacional:

1 𝜕(𝑟𝑣𝜃 ) 𝜕𝑣𝑟
𝜔𝑧 = � − �
2𝑟 𝜕𝜕 𝜕𝜕

1 𝜕(𝑟(−4𝑟𝑟𝑟𝑟(2𝜃))) 𝜕(4𝑟𝑟𝑟𝑟(2𝜃))
𝜔𝑧 = � − �
2𝑟 𝜕𝜕 𝜕𝜕

1
𝜔𝑧 = − 4𝑠𝑠𝑠(2𝜃)(2𝑟) + 4(𝑠𝑠𝑠(2𝜃))
2𝑟

𝜔𝑧 = −4𝑠𝑠𝑠(2𝜃) + 4𝑠𝑠𝑠(2𝜃) = 0

El rotacional es cero, por lo tanto el flujo es irrotacional y la función

potencial existe. Se deben cumplir las condiciones:

𝜕𝜕
𝑣𝑟 =
𝜕𝜕
1 𝜕𝜕
𝑣𝜃 =
𝑟 𝜕𝜕
𝜕𝜕
= 𝑣𝑟 = −4𝑟𝑟𝑟𝑟(2𝜃)
𝜕𝜕

Integrando con respecto a r se obtiene:

𝜙 = −2𝑟 2 𝑐𝑐𝑐(2𝜃) + 𝑓(𝜃)

38
 TEMA 6 FLUJO IDEAL D. Páez 2017

1 𝜕𝜕
Derivando con respecto a θ, e igualando a la condición: 𝑣𝜃 = , se
𝑟 𝜕𝜕

tiene:
𝜕𝜕 𝜕(−2𝑟 2 𝑐𝑐𝑐(2𝜃) + 𝑓(𝜃))
= = −2𝑟 2 �−2𝑠𝑠𝑠(2𝜃)� + 𝑓´(𝜃)
𝜕𝜕 𝜕𝜕

= 4𝑟 2 𝑠𝑠𝑠(2𝜃) + 𝑓´(𝜃)

𝜕𝜕
= 𝑟𝑣𝜃 = 𝑟(4𝑟𝑟𝑟𝑟(2𝜃)) = 4𝑟 2 𝑠𝑠𝑠(2𝜃))
𝜕𝜕

Igualamos las dos expresiones:

4𝑟 2 𝑠𝑠𝑠(2𝜃) + 𝑓´(𝜃) = 4𝑟 2 𝑠𝑠𝑠(2𝜃))

𝑓´(𝜃) = 0

𝑓(𝜃) = 𝐶

Sustituimos 𝑓(𝜃) = 𝐶 en la expresión de la función potencial:

𝝓 = −𝟐𝒓𝟐 𝒄𝒄𝒄(𝟐𝟐) + 𝑪

b) Aplicamos la ecuación de Euler, ecuación (6.54), entre los puntos A y

B:

𝑉𝐴2 𝑝𝐴 𝑉𝐵2 𝑝𝐵
+ = +
2 𝜌 2 𝜌

Se obtiene una expresión para la presión en el punto B:

𝒑𝑩 𝑽𝟐𝑨 𝒑𝑨 𝑽𝟐𝑩
= + −
𝝆 𝟐 𝝆 𝟐

39
 TEMA 6 FLUJO IDEAL D. Páez 2017

Se calcula la magnitud de la velocidad en los puntos A y B, como:

𝑽 = �𝑽𝟐𝒓 + 𝑽𝟐𝜽

En el punto A, en coordenadas cartesianas (x=1m, y =0), corresponde a

polares:

𝑟2 = 𝑥2 + 𝑦2 = 1 + 0 = 1

𝒓 = 𝟏𝟏

𝑦 0
𝜃 = 𝑎𝑎𝑎𝑎𝑎 � � = 𝑎𝑎𝑎𝑎𝑎 � � = 0
𝑥 1

𝜽=𝟎

Sustituyendo las coordenadas (r,θ) en la velocidad:

𝑽𝑨 = −4(1)𝑐𝑐𝑐(0)𝑟 + 4(1)𝑠𝑠𝑠(0)𝜃 = −4𝑟

𝑉𝐴 = �(−4)2 + 0 = 4 𝑚/𝑠

El punto B, en coordenadas cartesianas (x=0, y =0.5m), corresponde a

polares:

𝑟 2 = 𝑥 2 + 𝑦 2 = 0 + (0.5)2

𝒓 = 𝟎. 𝟓𝟓

𝑦 0.5 𝜋
𝜃 = 𝑎𝑎𝑎𝑎𝑎 � � = 𝑎𝑎𝑎𝑎𝑎 � � = 𝑎𝑎𝑎𝑎𝑎(∞) =
𝑥 0 2

40
 TEMA 6 FLUJO IDEAL D. Páez 2017

𝝅
𝜽=
𝟐

Sustituyendo las coordenadas (r,θ) en la velocidad:

𝜋 𝜋
𝑉𝐵 = −4(0)𝑐𝑐𝑐( )𝑟 + 4(0.5)𝑠𝑠𝑠 � � = 2𝜃
2 2 𝜃

𝑉𝐵 = �(0)2 + 22 = 2 𝑚/𝑠

𝑽𝑩 = 𝟐 𝒎/𝒔

Sustituyendo en la expresión para la presión en B:

𝑝𝐵 16 30000 4
= + −
𝜌 2 𝜌 2

𝑝𝐵 = (6 ∗ 1000 + 30000)𝑃𝑃 = 36𝐾𝐾𝐾

𝒑𝑩 = 𝟑𝟑𝟑𝟑𝟑

Podemos hallar una expresión para la magnitud de la velocidad en forma

general para cualquier punto del campo de flujo:

𝑉 = �𝑉𝑟2 + 𝑉𝜃2

2 2
𝑉 = �(−4𝑟𝑟𝑟𝑟(2𝜃)) + (4𝑟𝑟𝑟𝑟(2𝜃))

2 2
𝑉 = 4𝑟�(𝑐𝑐𝑐(2𝜃)) + (𝑠𝑠𝑠(2𝜃)) = 4𝑟

𝑽 = 𝟒𝟒

41
 TEMA 6 FLUJO IDEAL D. Páez 2017

FLUJOS POTENCIALES BIDIMENSIONALES

A continuación estudiaremos algunos casos de flujos potenciales

bidimensionales, como son: el flujo rectilíneo uniforme, la fuente, el sumidero, el

vórtice irrotacional y el flujo alrededor de una esquina.

Flujo Rectilíneo Uniforme

El flujo rectilíneo uniforme, que forma un ángulo α con la horizontal, se

puede representar mediante el siguiente vector velocidad en coordenadas

cartesianas, a partir del cual se pueden hallar la función de corriente y la función

potencial.

�𝑽⃗ = 𝑼𝑼𝑼𝑼(𝜶)𝒊 + 𝑼𝑼𝑼𝑼(𝜶)𝒋 (𝟔. 𝟓𝟓)

U = magnitud del vector velocidad

La función potencial puede determinarse a partir de las ecuaciones:

𝝏𝝏
=𝒖 (𝟔. 𝟐𝟐)
𝝏𝝏

𝝏𝝏
=𝒗 (𝟔. 𝟐𝟐)
𝝏𝝏

Obteniéndose:

𝝓 = 𝑼𝑼𝑼𝑼(𝜶)𝒙 + 𝑼𝑼𝑼𝑼(𝜶)𝒚 (𝟔. 𝟓𝟓)

α constante.

42
 TEMA 6 FLUJO IDEAL D. Páez 2017

Similarmente la función de corriente puede hallarse haciendo uso de las

ecuaciones:

𝟏 𝝏𝝏
𝒗𝒓 = (𝟔. 𝟒𝟒)
𝒓 𝝏𝝏

𝝏𝝏
𝒗𝜽 = − (𝟔. 𝟒𝟏)
𝝏𝝏

Obteniéndose:

𝝍 = −𝑼𝑼𝑼𝑼(𝜶)𝒙 + 𝑼𝑼𝑼𝑼(𝜶)𝒚 (𝟔. 𝟓𝟓)

A partir de la ecuación (6.55) se puede determinar la pendiente de las líneas

de corriente:

𝑑𝑑 𝑣 𝑈𝑈𝑈𝑈(𝛼)
= = = 𝑡𝑡𝑡(𝛼)
𝑑𝑑 𝑢 𝑈𝑈𝑈𝑈(𝛼)

𝑑𝑑
= 𝑡𝑡𝑡(𝛼) = 𝑐𝑐𝑐𝑐𝑡𝑡𝑡𝑡𝑡
𝑑𝑑

La pendiente de las líneas de corriente es constante, por lo tanto se trata de

líneas rectas con pendiente igual a la tangente del ángulo de inclinación, α.

En la figura 6.8 se presenta el flujo rectilineo representado por líneas de

corriente rectilineas función del ángulo α

43
 TEMA 6 FLUJO IDEAL D. Páez 2017

U
α
x

Figura 6.8 Flujo Rectilíneo Uniforme

Fuente y Sumidero
La fuente se trata de un flujo radial hacia afuera desde una línea a través de

un origen perpendicular al plano x-y. Se tiene un flujo únicamente radial.

Se define q = flujo volumétrico por unidad de longitud.

𝒒 = 𝑽𝒓 𝒂𝒂𝒂𝒂 = 𝑽𝒓 (𝟐𝝅𝝅) (𝟔. 𝟓𝟓)

𝒒 = 𝒄𝒄𝒄𝒄𝒄𝒄𝒄𝒄𝒄 = 𝑽𝒓 (𝟐𝝅𝝅) = 𝟐𝟐𝟐 (𝟔. 𝟓𝟓)

m = constante y se denomina intensidad de la fuente.

A partir de la ecuación (6.59) se obtiene que la velocidad radial sea:

𝒎
𝑽𝒓 = (𝟔. 𝟔𝟔)
𝒓

Como se trata de un flujo radial únicamente, la componente del vector

velocidad en la dirección tangencial es nula:

44
 TEMA 6 FLUJO IDEAL D. Páez 2017

𝑽𝜽 = 𝟎 (𝟔. 𝟔𝟔)

El vector velocidad en coordenadas polares será para una fuente resulta:

𝒎
�⃗ =
𝑽 + 𝟎𝒋 (𝟔. 𝟔𝟔)
𝒓𝒓

A partir del vector velocidad, ecuación (6.62), se obtiene la función

potencial:

𝝓 = 𝒎 𝐥𝐥(𝒓) (𝟔. 𝟔𝟔)

De forma similar se puede determinar la función de corriente a partir del

vector velocidad:

𝝍 = 𝒎𝒎 (𝟔. 𝟔𝟔)

En la figura 6.9 se muestran las líneas de corriente y las líneas

equipotenciales para el caso de una fuente.

Para el flujo ideal denominado sumidero, el flujo es igualmente radial pero

en este caso la dirección del flujo es hacia adentro, hacia el origen. Por lo que

para el sumidero se tienen las mismas expresiones que para la fuente, donde la

intensidad es negativa, m =-m.

45
 TEMA 6 FLUJO IDEAL D. Páez 2017

Figura 6.9. Fuente representada por las líneas de corriente y líneas

equipotenciales.

Vórtice Irrotacional ó Vórtice Libre

El vórtice irrotacional como su nombre lo indica es un flujo. Se trata de un

campo de flujo para el cual las líneas de corriente son círculos concéntricos. El

vórtice es el flujo con velocidad radial nula y con una velocidad tangencial

función de la dirección radial:

𝑉𝑟 = 𝑓(𝑟)

En la figura 6.10 se representan las líneas de corriente para un vórtice

irrotacional.

46
 TEMA 6 FLUJO IDEAL D. Páez 2017

Vθ ψ

Figura 6.10. Vórtice representado por líneas de corriente.

Se tiene un flujo circular sin flujo radial, por lo que las componentes del

vector velocidad resultan:

𝑉𝑟 = 0
𝐶
𝑉𝜃 =
𝑟
𝐶 = 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐

El vector velocidad es:

𝑪
�𝑽⃗ = 𝟎𝒓 + (𝟔. 𝟔𝟔)
𝒓𝒋

La función potencial resulta:

𝝓=𝑪∗𝜽 (𝟔. 𝟔𝟔)

La función de corriente está dada por:

𝝍 = −𝑪 ∗ 𝒍𝒍(𝒓) (𝟔. 𝟔𝟔)

47
 TEMA 6 FLUJO IDEAL D. Páez 2017

𝐪
𝑪 = 𝒄𝒄𝒄𝒄𝒄𝒄𝒄𝒄𝒄 = (𝟔. 𝟔𝟔)
𝟐𝟐

Se define la circulación como, Γ, la integral lineal de la componente

tangencial de la velocidad alrededor de una curva cerrada. Esta trayectoria a lo

largo de un círculo con centro en el origen se puede escribir para el caso del

vórtice irrotacional como:

𝟐𝟐
𝚪 = � (𝑽𝜽 𝒓𝒓𝒓) = 𝟐𝟐𝟐𝑽𝜽
𝟎

𝑪
𝚪 = 𝟐𝟐𝟐𝑽𝜽 = 𝟐𝟐𝟐 =𝒒 𝒄𝒄𝒄𝒄𝒄𝒄 𝒖𝒖𝒖𝒖𝒖𝒖𝒖𝒖 (𝟔. 𝟔𝟔)
𝒓

Ejemplo 6. Superficie libre de un vórtice irrotacional

Se tiene un flujo a través de un orificio en el fondo de una tanque, como se

muestra en la figura 6.11. Si se forma un vórtice libre determine una expresión

para la superficie libre.

Solución:

La función potencial para el vórtice libre está dada por la ecuación (6.57),

con:
𝑞
𝐶=
2𝜋

48
 TEMA 6 FLUJO IDEAL D. Páez 2017

1 Z=0

Figura 6.11. Superficie de un vórtice libre.

𝒒
𝝓= 𝜽 (𝟔. 𝟔𝟔)
𝟐𝟐

La ecuación de energía entre los puntos 1 y 2:

𝑽𝟏 𝟐 𝒑 𝟏 𝑽𝟐 𝟐 𝒑 𝟐
+ + 𝒈𝒁𝟏 = + + 𝒈𝒁𝟐
𝟐 𝝆 𝟐 𝝆

Se tiene los puntos sobre la superficie libre por tanto la presión relativa es

cero:

𝑝1 = 𝑝2 = 0, 𝑍1 = 0

𝑉1 2 𝑉2 2
= + 𝑔𝑍2 𝑍2 = 𝑐𝑐𝑐𝑐 𝑑𝑑 𝑙𝑙 𝑠𝑠𝑝𝑒𝑒𝑒𝑒𝑒𝑒𝑒 = 𝑍𝑠
2 2

𝑉1 2 𝑉2 2
𝑍𝑠 = − 𝑍2 = 𝑍𝑠 = 𝑐𝑐𝑐𝑐 𝑑𝑑 𝑙𝑙 𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠
2𝑔 2𝑔

49
 TEMA 6 FLUJO IDEAL D. Páez 2017

La magnitud de la velocidad para el vórtice está dada por:

𝑉 = �𝑉𝑟 2 + 𝑉𝜃 2 = �0 + 𝑉𝜃 2 = 𝑉𝜃

𝑪 𝒒
𝑽𝜽 = =
𝒓 𝟐𝟐𝟐
𝒒
𝑽=
𝟐𝟐𝟐

En el punto 1, se toma 𝑟1 → ∞:

𝑞
𝑉1 = ≅0
2𝜋𝑟1

En el punto 2, se toma 𝑟2 = 𝑟:

𝑞 𝑞
𝑉2 = =
2𝜋𝑟2 2𝜋𝜋

En la ecuación de la superficie libre:

𝑉1 2 𝑉2 2 1 𝑞 2 𝑞2
𝑍𝑠 = − =0− � � =−
2𝑔 2𝑔 2𝑔 2𝜋𝜋 8𝑔𝜋 2 𝑟 2

𝒒𝟐
𝒁𝒔 = −
𝟖𝟖𝝅𝟐 𝒓𝟐

La ecuación de la superficie libre no es aplicable en el origen.

50
 TEMA 6 FLUJO IDEAL D. Páez 2017

Flujo alrededor de una esquina

El flujo alrededor de una esquina está representado por las siguientes

funciones de corriente y potencial:

𝝍 = 𝑼𝒓𝝅/𝜶 𝒔𝒔𝒔(𝝅𝝅/𝜶) (𝟔. 𝟕𝟕)

𝝓 = 𝑼𝒓𝝅/𝜶 𝒄𝒄𝒄(𝝅𝝅/𝜶) (𝟔. 𝟕𝟕)

Donde α es el ángulo que forma la esquina, es constante, r y θ variables.

A partir de la ecuación (6.70) o de la ecuación (6.71) se puede hallar el vector

velocidad para el flujo alrededor de una esquina.

Se muestran las líneas de corriente y equipotenciales en la figura 6.12.

φ3

φ2
φ1
φ0
ψ3

r ψ2
α ψ1
θ
ψ0

Figura 6.12. Flujo alrededor de una esquina.

51
 TEMA 6 FLUJO IDEAL D. Páez 2017

SUPERPOSICION DE FLUJOS POTENCIALES BIDIMENSIONALES

Se pueden describir flujos más complejos superponiendo las funciones

potenciales o las funciones de corriente de flujos simples.

Doblete
Un doblete es el flujo resultante de la superposición de una fuente y un

sumidero con igual intensidad cuando la distancia entre ellos tiende a cero.

Recordemos para la fuente:

𝜙𝑓 = 𝑚ln(𝑟1 ) y 𝜓𝑓 = 𝑚𝜃1

para el sumidero: 𝜙𝑠 = −𝑚ln(𝑟2 ) y 𝜓𝑠 = −𝑚𝜃2

𝑞
se tiene: 𝑚=
2𝜋

y P

r1ε r2
r

θ2
θ1 θ
x

-a a

Figura 6.13. Representación de una fuente y un sumidero equidistantes del

origen.

52
 TEMA 6 FLUJO IDEAL D. Páez 2017

La fuente y el sumidero se encuentran equidistantes del origen, a una

distancia “a”. Se tiene una fuente en (-a,0) y un sumidero en (a,0). Ver figura

6.13.

De acuerdo a la figura 6.13, por la ley del coseno, se escribe:

𝒓𝟏 𝟐 = 𝒓𝟐 + 𝒂𝟐 − 𝟐𝟐𝟐𝟐𝟐𝟐(𝜽) (𝟔. 𝟕𝟕)

𝑟2 2 = 𝑟 2 + 𝑎2 − 2𝑎𝑎𝑎𝑎𝑎(𝜋 − 𝜃)

𝒓𝟐 𝟐 = 𝒓𝟐 + 𝒂𝟐 + 𝟐𝟐𝟐𝟐𝟐𝟐(𝜽) (𝟔. 𝟕𝟕)

La función potencial resultante de superponer la fuente y el sumidero:

𝜙 = 𝜙𝑓 + 𝜙𝑠 = 𝑚 ln(𝑟1 ) − 𝑚ln(𝑟2 )

𝝓𝑫 = 𝒎 𝐥𝐥(𝒓𝟏 ) − 𝒎𝐥 𝐧(𝒓𝟐 ) (𝟔. 𝟕𝟕)

La ecuación (6.74) se puede escribir como:

𝑟1 𝑚 𝑟1 2
𝜙𝐷 = 𝑚 ln(𝑟1 ) − 𝑚l n(𝑟2 ) = 𝑚𝑚𝑚 � � = 𝑙𝑙 � �
𝑟2 2 𝑟2

Manipulando la ecuación:

𝑚 𝑟 2 + 𝑎2 − 2𝑎𝑎𝑎𝑎𝑎(𝜃)
𝜙𝐷 = 𝑙𝑙 � 2 �
2 𝑟 + 𝑎2 + 2𝑎𝑎𝑎𝑎𝑎(𝜃)

𝑚
= �𝑙𝑙�𝑟 2 + 𝑎2 − 2𝑎𝑎𝑎𝑎𝑎(𝜃)� − 𝑙𝑙�𝑟 2 + 𝑎2 + 2𝑎𝑎𝑎𝑎𝑎(𝜃)��
2

53
 TEMA 6 FLUJO IDEAL D. Páez 2017

𝑚 2𝑎𝑎𝑎𝑎𝑎(𝜃)
𝜙𝐷 = �𝑙𝑙 �(𝑟 2 + 𝑎2 ) �1 − 2 ��
2 (𝑟 + 𝑎2 )

2𝑎𝑎𝑎𝑎𝑎(𝜃)
− 𝑙𝑙 �(𝑟 2 + 𝑎2 ) �1 + ���
(𝑟 2 + 𝑎2 )

𝒎 𝟐𝟐𝟐𝟐𝟐𝟐(𝜽) 𝟐𝟐𝟐𝟐𝟐𝟐(𝜽)
𝝓𝑫 = �𝒍𝒍 �𝟏 − � − 𝒍𝒍 �𝟏 + �� (𝟔. 𝟕𝟕)
𝟐 (𝒓𝟐 + 𝒂𝟐 ) (𝒓𝟐 + 𝒂𝟐 )

La ecuación (6.75) se puede escribir en forma de series:

𝒙𝟐 𝒙𝟑 𝒙𝟒
𝐥𝐥(𝟏 + 𝒙) = 𝒙 − + − + ⋯ (𝟔. 𝟕𝟕)
𝟐 𝟑 𝟒

𝒙𝟐 𝒙𝟑 𝒙𝟒
𝐥𝐥(𝟏 − 𝒙) = −𝒙 + − + − ⋯ (𝟔. 𝟕𝟕)
𝟐 𝟑 𝟒

Se aplican las series (6.76) y (6.77) a la ecuación (6.75), se toma el primer

término de la serie:

𝑚 2𝑎𝑎𝑎𝑎𝑎(𝜃) 2𝑎𝑎𝑎𝑎𝑎(𝜃)
𝜙𝐷 = �𝑙𝑙 �1 − 2 � − 𝑙𝑙 �1 + ��
2 (𝑟 + 𝑎2 ) (𝑟 2 + 𝑎2 )

𝑚 2𝑎𝑎𝑎𝑎𝑎(𝜃) 2𝑎𝑎𝑎𝑎𝑎(𝜃) 𝑚 4𝑎𝑎𝑎𝑎𝑎(𝜃)

= �− 2 − � = �− �
2 (𝑟 + 𝑎2 ) (𝑟 2 + 𝑎2 ) 2 (𝑟 2 + 𝑎2 )

𝒓𝒓𝒓𝒓(𝜽)
𝝓𝑫 = −𝟐𝟐𝟐 � 𝟐 � (𝟔. 𝟕𝟕)
(𝒓 + 𝒂𝟐 )

Se llama la intensidad del doblete:

54
 TEMA 6 FLUJO IDEAL D. Páez 2017

𝝀 = 𝟐𝟐𝟐

En el límite cuando la fuente y el sumidero se acercan, es decir cuando a

tiende a cero, la ecuación (6.78) resulta:

𝑟𝑟𝑟𝑟(𝜃) 𝑐𝑐𝑐(𝜃)
𝜙𝐷 = lim −𝜆 � 2 � = −𝜆 � �
𝑎→0 (𝑟 + 𝑎2 ) 𝑟

La ecuación (6.79) representa la función potencial para el doblete:

𝒄𝒄𝒄(𝜽)
𝝓𝑫 = −𝝀 � � (𝟔. 𝟕𝟕)
𝒓

A partir de la función potencial podemos hallar el vector velocidad:

𝝏𝝏
𝒗𝒓 = (𝟔. 𝟒𝟒)
𝝏𝝏

𝟏 𝝏𝝏
𝒗𝜽 = (𝟔. 𝟒𝟒)
𝒓 𝝏𝝏

La velocidad radial:

𝑐𝑐𝑐(𝜃)
𝜕𝜕 𝜕 �−𝜆 � 𝑟 �� 1 𝜆𝜆𝜆𝜆(𝜃)
𝑣𝑟 = = = −𝜆𝜆𝜆𝜆(𝜃) �− 2 � =
𝜕𝜕 𝜕𝜕 𝑟 𝑟2

𝝀𝝀𝝀𝝀(𝜽)
𝒗𝒓 = (𝟔. 𝟖𝟖)
𝒓𝟐

La velocidad tangencial:

1 𝜕𝜕 1 𝜕 𝑐𝑐𝑐(𝜃) −𝑠𝑠𝑠(𝜃)
𝑣𝜃 = = �−𝜆 � �� = −𝜆 � �
𝑟 𝜕𝜕 𝑟 𝜕𝜕 𝑟 𝑟2

55
 TEMA 6 FLUJO IDEAL D. Páez 2017

𝝀𝝀𝝀𝝀(𝜽)
𝒗𝜽 = (𝟔. 𝟖𝟖)
𝒓𝟐

La función de corriente se puede obtener a partir de las componentes del

vector velocidad, ecuaciones (6.80) y (6.81):

𝟏 𝝏𝝏
𝒗𝒓 = (𝟔. 𝟒𝟒)
𝒓 𝝏𝝏

1 𝜕𝜕 𝜆𝜆𝜆𝜆(𝜃)
𝑣𝑟 = =
𝑟 𝜕𝜕 𝑟2

𝜕𝜕 𝜆𝜆𝜆𝜆(𝜃)
= , 𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖𝑖 𝑐𝑐𝑐 𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟 𝑎 𝜃:
𝜕𝜕 𝑟

𝜆𝜆𝜆𝜆(𝜃)
𝜓= + 𝑓(𝑟) 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑 𝑐𝑐𝑐 𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟 𝑎 𝑟:
𝑟

𝜕𝜕 𝜕 𝜆𝜆𝜆𝜆(𝜃) 𝜆𝜆𝜆𝜆(𝜃)
= � + 𝑓(𝑟)� = − + 𝑓`(𝑟)
𝜕𝜕 𝜕𝜕 𝑟 𝑟2

Igualando a la velocidad tangencial:

𝝏𝝏
𝒗𝜽 = − (𝟔. 𝟒𝟒)
𝝏𝝏

𝜕𝜕 𝜕 𝜆𝜆𝜆𝜆(𝜃) 𝜆𝜆𝜆𝜆(𝜃) 𝜆𝜆𝜆𝜆(𝜃)

= � + 𝑓(𝑟)� = − + 𝑓`(𝑟) = −𝑣𝜃 = −
𝜕𝜕 𝜕𝜕 𝑟 𝑟2 𝑟2

𝜆𝜆𝜆𝜆(𝜃) 𝜆𝜆𝜆𝜆(𝜃)
− + 𝑓`(𝑟) = −
𝑟2 𝑟2

𝑓`(𝑟) = 0

56
 TEMA 6 FLUJO IDEAL D. Páez 2017

𝑓(𝑟) = 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 = 𝐶

La ecuación (6.82) representa la función de corriente para el doblete:

𝝀𝝀𝝀𝝀(𝜽)
𝝍= +𝑪 (𝟔. 𝟖𝟖)
𝒓
𝒒
𝝀 = 𝟐𝟐𝟐 = 𝒂 (𝟔. 𝟖𝟖)
𝝅

En la figura 6.14 se muestran las líneas de corriente y equipotenciales para

un doblete.

Figura 6.14. Líneas de corriente y equipotenciales para un Doblete

57
 TEMA 6 FLUJO IDEAL D. Páez 2017

Flujo alrededor de un cilindro

Se obtiene superponiendo: Flujo Uniforme horizontal +Doblete

En la figura 6.15 se muestra el flujo resultante. La función potencial

resultante será:

𝜙 = 𝜙𝑈𝑈𝑈𝑈𝑈𝑈𝑈𝑈 + 𝜙𝐷𝐷𝐷𝐷𝐷𝐷𝐷

𝐾
𝜙 = 𝑈. 𝑥 + cos(𝜃)
𝑟

Hacemos: 𝐾 = 𝑈𝑅 2

donde U= velocidad uniforme, K= intensidad del doblete

𝑼𝑹𝟐
𝝓 = 𝑼𝑼 + 𝐜𝐜𝐜(𝜽) (𝟔. 𝟖𝟖)
𝒓

De forma similar la función de corriente resulta en:

𝐾
𝜓 = 𝑈𝑈 − 𝑠𝑠𝑠(𝜃)
𝑟

𝑼𝑹𝟐
𝝍 = 𝑼𝑼 − 𝒔𝒔𝒔(𝜽) (𝟔. 𝟖𝟖)
𝒓

Estas expresiones (6.84) y (6.85) son válidas para un radio r>a, en otras

palabras se pueden aplicar fuera del contorno sólido. Cuando el flujo se acerca al

cuerpo se tiene en R = a un punto de velocidad cero.

58
 TEMA 6 FLUJO IDEAL D. Páez 2017

ψ
U

R=a
x

Figura 6.15. Flujo alrededor de un cilindro.

Ejemplo 6. Flujo alrededor de un cilindro

Para el flujo alrededor de un cilindro:

a) Determine el vector velocidad en coordenadas polares.

b) Determine la magnitud del vector velocidad en la superficie del
cilindro.
c) Compruebe la ecuación de conservación.

Solución:

A partir de las ecuaciones de función potencial o función de corriente

podemos hallar las componentes del vector velocidad.

59
 TEMA 6 FLUJO IDEAL D. Páez 2017

𝑈𝑈2
𝜙 = 𝑈. 𝑥 + cos(𝜃)
𝑟

En coordenadas polares la función potencial resulta en:

𝑈𝑈2
𝜙 = 𝑈𝑈𝑈𝑈𝑈(𝜃) + cos(𝜃)
𝑟

𝜕𝜕 1 𝜕𝜕
𝑣𝑟 = 𝑣𝜃 =
𝜕𝜕 𝑟 𝜕𝜕

a) Hallamos las correspondientes derivadas:

𝑈𝑈2
𝜕(𝑈𝑈𝑈𝑈𝑈(𝜃) + cos(𝜃)) 𝑅2
𝑣𝑟 = 𝑟 = 𝑈𝑈𝑈𝑈(𝜃) �1 − 2 �
𝜕𝜕 𝑟

𝑈𝑈2 𝑅2
1 𝜕 �𝑈𝑈𝑈𝑈𝑈(𝜃) + 𝑟
cos(𝜃)� 1 𝜕 �𝑈𝑈𝑈𝑈(𝜃) �𝑟 + 𝑟 ��
𝑣𝜃 = =
𝑟 𝜕𝜕 𝑟 𝜕𝜕

1 𝑅2 𝑅2
= �−𝑈𝑈𝑈𝑈 (𝜃) �𝑟 + �� = −𝑈𝑈𝑈𝑈(𝜃) �1 + 2 �
𝑟 𝑟 𝑟

Obtenemos las coordenadas del vector velocidad:

𝑅2
𝑣𝑟 = 𝑈𝑈𝑈𝑈(𝜃) �1 − 2 �
𝑟

𝑅2
𝑣𝜃 = −𝑈𝑈𝑈𝑈(𝜃) �1 + 2 �
𝑟

b) La magnitud del vector velocidad:

60
 TEMA 6 FLUJO IDEAL D. Páez 2017

𝑉 = �𝑉𝑟2 + 𝑉𝜃2

Sustituimos las expresiones para las componentes del vector velocidad:

2 2
𝑅2 𝑅2
𝑉 = ��𝑈𝑈𝑈𝑈(𝜃) �1 − � � + �−𝑈𝑈𝑈𝑈(𝜃) �1 + ��
𝑟2 𝑟2

En la superficie del cilindro r = R

2 2
𝑅2 𝑅2
𝑉 = ��𝑈𝑈𝑈𝑈(𝜃) �1 − � � + �−𝑈𝑈𝑈𝑈(𝜃) �1 + ��
𝑅2 𝑅2

2
𝑅2
𝑉 = ��−𝑈𝑈𝑈𝑈(𝜃) �1 + ��
𝑅2

𝑉 = 2𝑈𝑈𝑈𝑈(𝜃)

𝑉𝜃 = −2𝑈𝑈𝑈𝑈(𝜃)

𝑉𝑟 = 0

En la superficie del cuerpo se tiene solamente velocidad tangencial, ya que

no hay flujo a través del cilindro! La velocidad radial es cero. Ver figura 6.15.

61
 TEMA 6 FLUJO IDEAL D. Páez 2017

Flujo alrededor de medio cuerpo

El flujo alrededor de medio cuerpo resulta de superponer dos flujos: Flujo

Uniforme +Fuente.

La función potencial resultante de la superposición de los flujos, resulta:

𝜙 = 𝜙𝑈𝑈𝑈𝑈𝑈𝑈𝑈𝑈 + 𝜙𝐹𝐹𝐹𝐹𝐹𝐹

𝜙 = 𝑈𝑈 + 𝑚ln(𝑟)

La escribimos en coordenadas polares:

𝝓 = 𝑼𝑼𝑼𝑼𝑼(𝜽) + 𝒎𝐥 𝐧(𝒓) (𝟔. 𝟖𝟖)

La función de corriente resultante será:

𝜓 = 𝜓𝑈𝑈𝑈𝑈𝑈𝑈𝑈𝑈 + 𝜓𝐹𝐹𝐹𝐹𝐹𝐹

𝜓 = 𝑈𝑈 + 𝑚𝑚

La escribimos en coordenadas polares:

𝝍 = 𝑼𝑼𝑼𝑼𝑼(𝜽) + 𝒎𝛉 (𝟔. 𝟖𝟖)

Para el flujo alrededor de medio cuerpo la distancia a la cual se tiene el

contorno del cuerpo y la velocidad es nula, se denomina, b = punto de

estancamiento con velocidad cero. Ver figura 6.16.

62
 TEMA 6 FLUJO IDEAL D. Páez 2017

U
r

Figura 6.16. Flujo alrededor de medio cuerpo.

Ejemplo 7. Superposición de flujos ideales.

Determine el flujo resultante de superponer los flujos ideales

representados por las funciones de corriente:

𝐶 𝐶
𝜓1 = − ln(𝑟) vórtice, 𝑘 =
2𝜋 2𝜋

𝑞 −𝑞
𝜓2 = θ sumidero, 𝑚 =
2𝜋 2𝜋

Para los flujos ideales bidimensionales dados podemos superponer las

funciones de corriente para determinar el flujo resultante:

𝜓 = 𝜓1 +𝜓2

𝐶 𝑞
𝜓=− ln(𝑟) + θ
2𝜋 2𝜋

El vector velocidad resultante estará definido por:

63
 TEMA 6 FLUJO IDEAL D. Páez 2017

1 𝜕𝜕 𝜕𝜕
𝑣𝑟 = , 𝑣𝜃 = −
𝑟 𝜕𝜕 𝜕𝜕

𝐶 𝑞
𝜕𝜕 𝜕 �2𝜋 ln(𝑟) + 2𝜋 𝜃� 𝑞
= =
𝜕𝜕 𝜕𝜕 2𝜋

La componente tangencial será:

1 𝑞
𝑣𝑟 = −
𝑟 2𝜋

𝐶 𝑞
𝜕𝜕 𝜕 �2𝜋 ln(𝑟) + 2𝜋 𝜃� 𝐶 1
= =
𝜕𝜕 𝜕𝜕 2𝜋 𝑟

La componente normal será:

𝐶 1
𝑣𝜃 = −
2𝜋 𝑟

1 𝑞 𝐶 1
�⃗ = �−
𝑉 � + �− �
𝑟 2𝜋 𝑟 2𝜋 𝑟 𝜃

Este vector velocidad debe satisfacer continuidad y la condición de

irrotacionalidad.

64
 TEMA 6 FLUJO IDEAL D. Páez 2017

EJERCICIOS PROPUESTOS
6.1 Para el campo bidimensional dado por:
−𝐾𝐾
𝑢=
𝑥2+ 𝑦2
𝐾𝐾
𝑣=
𝑥2 + 𝑦2
Donde K = constante.
a) Indique si el vector velocidad representa un flujo conservativo.
b) Hallar el rotacional en coordenadas cartesianas y en coordenadas
polares.
c) Exprese el vector velocidad en sus componentes tangenciales y
normales (coordenadas polares).

𝑅2
6.2 Demuestre que el potencial: 𝜙 = 𝑈 �𝑟 + � cos(𝜃) y la función de
𝑟
𝑅2
corriente 𝜓 = 𝑈 �𝑟 + � sen(𝜃)
𝑟

a) Satisface la ecuación de continuidad

b) Representa un flujo irrotacional
El campo de velocidades satisface la condición de frontera:
𝑉𝑟 = 0 𝑒𝑒 𝑟 = 𝑅
𝑉𝑟 = 𝑈𝑈𝑈𝑈(𝜃) 𝑐𝑢𝑢𝑢𝑢𝑢 𝑟 → ∞
𝑉𝜃 = −𝑈𝑈𝑈𝑈(𝜃) 𝑐𝑐𝑐𝑐𝑐𝑐 𝑟 → ∞
6.3 Si se tiene un flujo representado por la función de corriente:
𝑞
𝜓 = 2 sen(2𝜃)
𝑟
a) Determine el vector velocidad
b) Indique si se trata de un flujo conservativo
c) Indique si representa un flujo irrotacional

65
 TEMA 6 FLUJO IDEAL D. Páez 2017

d) Determine la función potencial si existe

e) Si la presión en el punto (2,3) es de 180N/m2, determine la presión en el
punto (4,8). La densidad del fluido es de ρ = 1.2Kg/m3. Use q = 1m3/s/m

6.4 Se tiene la siguiente función: 𝜙 = −𝐾(𝑥 2 − 𝑦 2 )

a) Indique que tipo de flujo representa
b) Determine la función de corriente

6.5 Halle el flujo resultante de la superposición de un vórtice y de una fuente,

ambos ubicados en el origen, junto a un flujo uniforme horizontal hacia la
derecha.
a) Determine la función de corriente
b) Determine la función potencial
c) Halle una expresión para la magnitud del vector velocidad

6.6 Verificar si el vector:

�⃗ = 𝑎(𝑥 2 − 𝑦 2 )𝑖 − 2𝑎𝑎𝑎𝑗
𝑉
define el flujo bidimensional de un fluido ideal. En caso afirmativo:
a) Establecer las ecuaciones de la función de corriente y de la función
potencial si φ=0, ψ =0 en (1,1).
b) ¿Cuál es el caudal por unidad de ancho que fluye entre las líneas de
corriente que pasan por los puntos (3,2) y (6,4)

q
6.7 Dadas las funciones potenciales: ϕ1= -Ay, y ϕ2 = ( )ln(r), donde r2 =
2π

x2 + y2.

66
 TEMA 6 FLUJO IDEAL D. Páez 2017

a) Determinar la función potencial resultante de superponer los flujos

representados por las funciones potenciales dadas.

b) Hallar las componentes de la velocidad en coordenadas cartesianas y en

coordenadas cilíndrico polares en (4,2). Verifique que la magnitud de la
velocidad en ese punto es la misma en ambas coordenadas.
c) Determinar la función de corriente resultante.

5.8 Un flujo bidimensional está caracterizado por la función de corriente:

𝜓 = 2𝑦 − 3𝑥 −3xy
Determinar:
a) El rotacional.
b) El potencial, si existe.
c) Las componentes del vector velocidad.
d) El vector velocidad en el punto (2,5).

6.9 El flujo alrededor de un esquina, con ángulo = 135º, tiene una función de
q
corriente: ψ = sen(2θ)
r2

Determine:
a) Las componentes del vector velocidad.
b) La función potencial, si existe.
c) La presión en el punto (4,8), si la presión en (2,3) es 180N/m2. La
densidad del fluido es de ρ = 1.2kg/m3.

6.8 La función potencial de un flujo es ϕ = 9 + 6𝑥 − 4𝑦 + 7𝑥𝑥

a) Determine la función de corriente.

67
 TEMA 6 FLUJO IDEAL D. Páez 2017

b) Pruebe que se trata de un flujo ideal.

d) Determine el caudal entre (2,3) y (5,7).
e) Dibujar esas líneas de corriente y un conjunto de líneas equipotenciales.
f) Si la presión en el primer punto es: 100N/m2, determine la presión en
el segundo punto. La densidad del fluido es 1,3.

6.9 La función de corriente de un flujo es ψ= 6x2y – 2y3.

a) Hallar las expresiones para las componentes de la velocidad.
b) Determinar el modulo y ángulo del vector velocidad en el punto (2,3).
c) Hallar el rotacional.
e) Hallar la función potencial si existe.

6.10 Las líneas de corriente de cierto flujo permanente e incompresible están

expresadas por la ecuación: ψ= x2 – 2xy –y2. Hallar:
a) El vector velocidad.
b) Verificar si se cumple la ecuación de continuidad.
c) Determinar la vorticidad.
d) Determinar el potencial si existe.

6.11 Un flujo bidimensional está caracterizado por la función de corriente:

ψ= 2y – 3x – 3xy.
a) ¿Cuáles son las componentes de la velocidad en el punto (3,-1).
b) ¿Cuál es la vorticidad y cuál el potencial si existe?

68
 TEMA 6 FLUJO IDEAL D. Páez 2017

q
6.12 Superponer el flujo potencial de una fuente: ϕ1 = − ln R , con el de un
2π
Г
vórtice libre: ϕ2 = − θ, es decir ϕ = ϕ1 + ϕ2 , para hallar las
2π

expresiones de:
a) La vorticidad del flujo resultante.
b) Las velocidades uR Y vθ del flujo resultante.
c) La función de corriente ψ del flujo resultante.
d) La velocidad en R= 4 m, θ= π/4, si q= 8 m2/s y Г= 4 m2/s.

6.13 Se tiene un flujo ideal que resulta de la superposición de las funciones

potenciales: Φ1 = -Ax, y Φ2 = (q/2π)ln(r), A = 1, q = 3 m3/s.
Determine:
a) El vector velocidad resultante, tanto en coordenadas cartesianas como
polares.
b) La magnitud del vector velocidad en el punto (x = 1, y = 0.5).
c) Hallar la función de corriente, ψ, del flujo superpuesto.
d) Hallar el caudal que circula entre las líneas de corriente que pasan por el
punto (1, 0.5) y el punto (2, 1).

6.14 Para el flujo correspondiente al campo de velocidades:

u= 3y + 2, v= 2x - 1
a) Verificar la ecuación de continuidad.
b) Hallar el rotacional. Determine si el flujo es irrotacional o no.
c) Hallar el potencial, si existe.
d) Determinar la función de corriente.
e) Si la presión en el punto (4,5) es de 25 KPa determinar la presión en el
punto (8,7).

69
 TEMA 6 FLUJO IDEAL D. Páez 2017

REFERENCIAS
- Aguirre, J., et al. (1987). Mecánica de Fluidos Fundamental. Mérida. Venezuela.
Consejo de Publicaciones de la Universidad de Los Andes.
- Boxer, G. (1994). Mecánica de Fluidos. Wilmingon. DE: Addison-Wesley
Iberoamericana.
- Departamento de Ingeniería Civil, UNAH. (2016). Mecánica de Fluidos.
Recuperado desde:
https://mecanicadefluidosunah.files.wordpress.com/2015/11/moody.jpg.
7/10/2016

- Munson, B., et al. (1990). Fundamentals of Fluid Mechanics. New York. NY:
John Wiley and Sons.
- Olson, R. y Wright, S. (1990). Essentials of Engineering Fluid Mechanics. New
York. NY: Harper and Row.
- Redes de Flujo. Tomado de:
http://www.ingenierocivilinfo.com/2010/09/construccion-de-la-red-de-flujo.html

- Rouse, H. (1978). Elementary Mechanics of Fluids. NewYork, NY: Dover

Publications.
- White, F. (2008). Mecánica de Fluidos. Sexta Edición. México. McGraw
Hill.

70