Parte I.

Mecánica de fluidos:
Tema IV 01
ANÁLISIS DIMENSIONAL, SEMEJANZA, ESTUDIO SOBRE MODELOS

Análisis dimensional y leyes de semejanza. Parámetros adimensionales.

Teorema de Buckingham. Método de Ipsen. Estudio de modelos.
Números adimensionales de: Euler, Froude, Reynolds, Weber y Mach.
Modelado físico. El túnel aerodinámico.

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 1

Contenido
• Análisis dimensional y leyes de semejanza.
• Parámetros adimensionales.
• Teorema de Buckingham.
• Estudio de modelos.
• Números adimensionales de
™ Euler, Froude, Reynolds, Weber y Mach.
• Modelado físico.
™ Simulación parcial
• El túnel aerodinámico.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 2

1
 Objetivos docentes
Después de completar el Tema, el alumno deberá poder realizar las siguientes
actividades:
• Comprender más detalladamente las dimensiones, unidades y homogeneidad
dimensional de las ecuaciones.
• Conocer las aplicaciones y los beneficios del análisis dimensional.
• Conocer las distintas técnicas que se aplican en el análisis dimensional.
• Aplicar el teorema de Buckingham (variables repetitivas) para identificar
parámetros adimensionales.
• Entender los conceptos de similitud geométrica, cinemática y dinámica y
como aplicarlos al modelado experimental o al análisis teórico.
• Reconocer grupos adimensionales de uso corriente en la mecánica de fluidos.
• Conocer las principales características de un túnel aerodinámico.
• Resolver problemas relacionados con los contenidos del Tema.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 3

Bibliografía recomendada
• Çengel Cimbala, (1ra edic.) cap. VII, pp 269:305.
• Fox, (2da edic.) cap. VII, pp 305:322.
• Potter (3ra edic.), cap VI pp 209:332.
• Franzini (9na edic.) cap VII, pp 147:159.
• White, (5ta edic.), cap. V, pp 287:331.
• Crespo, (1ra edic.), cap. XV, pp 221:260.
• Barrero, (1ra edic.), cap. VI, pp 111:119.
• Streeter (8va edic.) cap IV, pp 162:179

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 4

2
 Resumen
Análisis dimensional, semejanza, estudio sobre modelos

Se estudiaron ya los métodos integrales y diferenciales para el análisis de flujos. Se presenta

aquí el tercer método de análisis: la experimentación apoyada por el análisis dimensional. El
objetivo fundamental aquí es ver si es viable el estudiar un problema mediante experimentación
sobre un modelo a escala, ya que muchas veces medir en un modelo real es difícil y caro.
Se discutirá como pueden escribirse de forma adimensional relaciones físicas familiares para
cumplir el principio de homogeneidad dimensional. Se presentará el Teorema π, una técnica
general para encontrar en forma sistemática parámetros adimensionales partiendo de las
variables dimensionales que rigen un cierto proceso físico. Se hará el análisis dimensional de las
ecuaciones de la mecánica de fluidos para obtener números adimensionales y finalmente se
presentan casos en donde la semejanza es sólo parcial, y que técnicas se aplican en ellos. .
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 5

Introducción a la hidrodinámica
Método de solución, técnicas básicas

Se pueden clasificar en dos: teóricas o analíticas, y experimentales.

Para las técnicas teóricas, y, dependiendo del punto de vista que se use
(Euleriano o Lagrangiano), los resultados surgen del correspondiente análisis.
Partícula fluida: Análisis diferencial
Estudio analítico: { Volumen de control: Análisis integral
Estudio experimental: Adquisición de datos

Técnica auxiliar útil en todos los casos: Análisis dimensional

El flujo en estudio debe satisfacer las leyes de

conservación de la mecánica en todos los casos.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 6

3
 Análisis dimensional y semejanza
Análisis dimensional
Es una de las técnicas básicas para el análisis de
problemas complejos en el área de
TERMOFLUIDODINÁMICA
Técnica matemática que se basa en el estudio de las dimensiones.
El análisis dimensional está relacionado con la semejanza, pero el
planteamiento es distinto.

La técnica de análisis dimensional consiste en reducir el número

de variables relacionadas con un fenómeno físico. Se intenta así,
encontrar una relación entre estas variables que permita
comprender mejor el fenómeno estudiado.

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 7

Análisis dimensional y semejanza

Análisis dimensional
Es muy útil en técnicas experimentales ya que permite definir y planificar
mejor el tipo de experimentación a realizar, ahorrando costos de ensayos al
minimizar la cantidad de experimentos a realizar. Permite analizar la viabilidad
de experimentar sobre un modelo a escala para hacer/optimizar el diseño de
un prototipo, permitiendo asegurar si hay condiciones de semejanza entre
ellos.
También se usa en el ámbito teórico para problemas en donde el fenómeno
físico o la geometría son complicados, las ecuaciones en forma integral no
puedan aportar los suficientes datos/detalles al resolverlas (cuando esto es
factible) y las ecuaciones diferenciales no puedan resolverse fácilmente.

En ambos casos, muchas veces se recurre a esta técnica, para

simplificar el problema, reteniendo sólo efectos que realmente son
importantes
17/05/2016 13:32 para el estudio
Mecánicaen cuestión.
de Fluidos – Miguel G. Coussirat Núñez 8

4
 Análisis dimensional y semejanza
Análisis dimensional
Magnitudes físicas:

* con dimensión:
• fundamentales: masa, distancia, temperatura...
• compuestas: densidad, velocidad, aceleración

* adimensionales: (compuestas, en virtud de su definición):

deformación específica, módulo de Poisson, densidad
relativa. Los ángulos son adimensionales (relación entre
arco y radio) y por eso se miden en radianes.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 9

Parámetros adimensionales
Fundamentos del análisis dimensional. Principio de Homogeneidad
dimensional:
Si una ecuación expresa correctamente una relación entre variables de un
proceso físico, ésta debe ser dimensionalmente homogénea.
Se puede deducir, usando este principio, si una relación entre variables es
correcta cuando uno no la recuerda.
Obsérvese que una integración o diferenciación de la ecuación que ya es
dimensionalmente homogénea puede cambiar sus dimensiones pero seguirá
siendo homogénea.
Ejemplo: Ecuación caída libre...
s = s0 + c0t + 0,5gt2 → dimensiones: [L]
∫sdt = s0t + 0,5c0t 2 + (1/6) gt3 → dimensiones: [Lt]
ds/dt = c0 + gt → dimensiones: [Lt-1]
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 10

5
 Parámetros adimensionales
Dimensiones fundamentales son propiedades fisicas que son
mensurables, por ejemplo longitud, tiempo, masa y temperatura.
Se usa un sistema de unidades para dar valores numéricos a
las dimensiones. Mills, “heat transfer”, Ed. Irwin, pg 45:
Los parámetros adimensionales (variables agrupadas) pueden
determinarse de dos formas:
• Pueden deducirse directamente de la ecuación integral/diferencial si
éstas son conocidas.
• Si las relaciones funcionales entre datos e incógnitas (ej.: ecuaciones
diferenciales) no son conocidas, los parámetros adimensionales
pueden determinarse haciendo un análisis dimensional de las variables
involucradas en el problema.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 11

Parámetros adimensionales
Objetivo: Reducir el número de variables independientes n (dimensionales)
de la que depende el problema a un número n-k, de variables (sin dimensión),
siendo j las dimensiones fundamentales y viendo que k ≤ j.
Técnicas de adimensionalización:
1. División por valores de referencia,
2. Producto de potencias,
3. Teorema π de Buckingham
A partir de la comprensión general del fenómeno se predicen los parámetros
físicos importantes y luego se agrupan en una serie de parámetros
adimensionales para estudiar con más detalle el fenómeno.
Magnitudes/Parámetros físicos importantes (n):
En mecánica de fluidos son: Velocidad, densidad, viscosidad, volumen, etc.
Dimensiones fundamentales (k): ...en sistemas
En mecánica de fluidos son: [M] ( o a veces [F] ), [L], [t], [T] eléctricos debe
incluirse la carga q
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 12

6
 Parámetros adimensionales
gS 0
División por valores de referencia: c02

Ejemplo:
Ecuación de caída libre en función del tiempo: s = s0 + c0t + 0,5 gt 2

Todos los términos de la ecuación son longitudes.

Variables dimensionales: s, t. Ctes dimensionales: g, c0 , s0. Ctes puras: 0,5 y 2. c0 t
S0

s ct gs s
Variables adimensionalizadas: s = ; t = 0 o bien s = 2 ; t = 0
* * * *

s0 s0 c0 c0t
6 4 4 4 447 4 4 4 4 48
s gs 2 2 gs 2
Variable t
s 0 s * = s 0 + c 0 0 t * + 0 , 5 20 t * ⇒ s * = 1 + t * + 0 , 5 20 t *
c0 c0 c0 dividiendo, no
{ puedo hacer t=0
α

⎛ gs ⎞
Se pasó de 5 cantidades: s = F(t, s0 , c0 , g) a sólo 3 cantidades s * = G ⎜⎜ t * , 20 ⎟⎟
⎝ c0 ⎠
En este caso la reducción k=j donde j número de dimensiones j [L][t] que hay en el problema, esto es, 5-3=2. Luego,
todas las variables tienen 2 dimensiones fundamentales [L][t] como máximo: [s,s0]= [L]; [t]=[t]; [g]=[Lt-2]; [c0]=[Lt-1].
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 13

Parámetros adimensionales
División por valores de referencia:
Ejemplo:
Ecuación de caída libre en función del tiempo: s = s0 + c0t + 0,5 gt 2 (Eq. 5.5)
Todos los términos de la ecuación son longitudes.
Variables dimensionales: s, t. Ctes dimensionales: g, c0 , s0. Ctes puras: 0,5 y 2.
gs tg
Variables adimensionalizadas: o bien : s** = 2
; t ** =
c0 c0
c0 gs0
( c0, g ) co2

c02 c0

c02

gt
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez t ** = 14
c0

7
 Parámetros adimensionales
División por valores de referencia:
Ejemplo:
Ecuación de caída libre en función del tiempo: s = s0 + c0t + 0,5 gt 2 (Eq. 5.5)
c0 .

c0
gs0

c0
gs0

c0

Fig 5.1c

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 15

Parámetros adimensionales
Producto de potencias:
Ejemplo: no se conoce la relación funcional entre variables. s = s0 + c0t + 0,5 gt
2

Si la intuición al realizar un experimento para medir como cambia s en el tiempo

nos indica que las variables importantes son:
s = F (t , s0 , c0 , g )
Si s tiene unidades de longitud [L] , luego s0, t, c0 y g deben combinarse de tal
forma que den unidades de longitud (homogeneidad dimensional).
Si la función F fuese de forma polinómica (ej.: ajuste puntos experimentales) y
a su vez cada monomio pudiese expresarse como producto de las cantidades t,
s0, c0, g, y siendo a, b, c y d coeficientes a determinar, F puede quedar
expresada en términos de [L] y/o [T]. En este caso debe quedar sólo la
dimensión longitud [L].

s = ∑ (cte)(t a s0 c0 g d ) ⇒ [ L]1 = [ t ]a [ L]b [ Lt−1 ]c [ Lt−2 ]d

b c

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 16

8
 Parámetros adimensionales
∴ s = ∑ (cte )(t a s0 c0 g d ) ⇒ [1 t ]0 = [1
]1[3 t ] a4
[ L4 [ Lt2−14
] b4 ]c [4 −2 d
]
b c
L2 Lt43
s ∑ (...)

Longitud 1 = b+c + d b =1 − a + d
Análisis
exponentes { Tiempo
{ 0 = a −c − 2d { c = a − 2d
Se tienen 2 ecuaciones con 4 incógnitas, se ponen dos incógnitas en función de las otras 2. Esto
es arbitrario y en función de que incógnitas se elija, cambiarán los grupos adimensionales.
∴ s = ∑ (cte)(t a s0 g d ) = ∑ (cte)(t a s0 s0 s0 c0 c0
1− a + d a −2d 1 −a d a −2 d
c0 gd )
⎡⎛ tc ⎞ a ⎛ s0 g ⎞ ⎤
d a
⎛ c0t ⎞ ⎛ gs0 ⎞
d

∴ s = ∑ (cte)⎢⎜ ⎟ (s0 )⎜ 2 ⎟ ⎥ = (cte )( s0 )∑ ⎜

⎜ 0
⎟ ⎜ ⎟ ⎜ ⎟ ⎜
⎟ ⎜ 2 ⎟ ⎟
⎢⎣⎝ s0 ⎠ ⎝ c0 ⎠ ⎥⎦ ⎝ s0 ⎠ ⎝ c0 ⎠
a d
⎡ c t ⎤ ⎡ gs ⎤ ⎛ c t gs ⎞
∴ = (cte)∑ ⎢ 0 ⎥ ⎢ 20 ⎥ → = G⎜⎜ 0 , 20 ⎟⎟ ⇒ s* = G (t*,α )
s s
s0 ⎣ s0 ⎦ ⎣ c0 ⎦ s0 ⎝ s0 c0 ⎠
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 17

Teorema π de Buckingham
El método de producto de potencias es de difícil aplicación cuando las
variables son muchas y las ecuaciones son pocas. El procedimiento alternativo
es el teorema π de Buckingham. La variable π significa un producto de
variables independientes que no tiene dimensiones.
Enunciado teorema π de Buckingham : “Si un proceso físico satisface el
principio de homogeneidad dimensional, se puede describir mediante
una relación entre sólo n-k variables adimensionales.
La reducción k es igual al máximo número de variables que no pueden
formar un grupo adimensional entre ellas y es siempre menor o igual al
número de dimensiones fundamentales j que describen estas variables
(k ≤ j)”, k se determina calculando el rango de una matriz.
No se debe omitir ninguna variable independiente importante, pues el análisis
evidentemente fallará debido a esa omisión.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 18

9
 Teorema π de Buckingham
Resumen
Se tienen F(x1,x2,...,xn)=0; xi variables dimensionales con i =1,...,n.
Se puede demostrar*Crespo, “Mecánica de Fluidos, pg 229” que pueden reorganizarse
como G(π1, π2,..., πn-k)=0; siendo los πn-k los n-k grupos adimensionales
independientes.
Siendo la reducción k como máximo, igual al número de dimensiones
fundamentales j para expresar todas las xi de forma dimensional.
La diferencia n-k se corresponde con el número de grupos
adimensionales independientes posibles π , esto es, los productos de
las k variables dimensionales elegidas entre las n variables xi , que se
repetirán en todos los grupos y elevadas a la potencia adecuada para que el
grupo πi no tenga dimensiones.
Las restante xi de las n variables x que aparece en cada grupo πi no se
eleva a ninguna potencia.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 19

Teorema π de Buckingham
Procedimiento
1. Haga una lista con los parámetros dimensionales del problema y
cuente su número, (n).
∴ Se tienen : F ( x1 , x2 ,..., xn ) = 0 ; xi con i = 1,.., n : dimensiona les

2. Haga una lista con las dimensiones primarias (j) de los n parámetros
dimensionales.

3. Calcule el rango de la matriz conformada por j filas y n columnas, el

que dará la reducción k. La variable k será el número de parámetros
repetititvos (variables repetitivas) en cada grupo π. Luego de acuerdo al
teorema π de Buckingham, el número esperado de grupos π es n-k.
G (π 1 , π 2 ,..., π n − k ) = 0 ;
Siendo π n − k los n − k grupos adimension ales independientes (la reducción).

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 20

10
 Teorema π de Buckingham
Procedimiento

4. Elija los k parámetros repetitivos que se usan para construir cada grupo
πi. Como estos parámetros aparecen en cada grupo, deben elegirse
adecuadamente (ver siguiente slide). Estos parámetros son los que están
elevados a las potencias necesarias para que el grupo sea adimensional,
luego de añadir cualquiera de los otros parámetros no seleccionados.

5. Se genera cada grupo π agrupando los k parámetros repetitivos

elegidos de los xi , con cada uno de los parámetros xi restantes de los n
parámetros dimensionales, que no se repetirán en cada grupo πi . Se
buscan los valores de las potencias a las que debe elevarse cada uno de
los k parámetros dimensionales, para que cada πi de los πn-k sea
adimensional.

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 21

Teorema π de Buckingham
Guía para elegir los parámetros que se repiten:
1. Nunca elija variables dependientes, ellas pueden luego aparecer en todos los grupos π.
2. Los parámetros a elegir no deben ser capaces de formar un grupo adimensional por si
mismo, pues puede que sea imposible generar el resto de los grupos π.
3. Los parámetros a elegir deben representar todas las dimensiones primarias observadas.
4. Nunca elija parámetros que ya son adimensionales.
5. Nunca elija dos parametros con las mismas dimensiones o con dimensiones que difieran en
sólo un exponente.
6. Elija constantes dimensionales sobre variables dimensional de modo tal que un grupo π
contenga la variable dimensional.
7. Elija parámetros comunes ya que pueden aparecer en cada uno de los grupos π.
8. Elija parámetros simples y no compuestos.

Guía para verificar los parámetros obtenidos:

1. Verifique que los grupos π.obtenidos sean efectivamente adimensionales
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 22

11
 Teorema π de Buckingham
Ejemplo 1:
Investigación dimensional de como varía la presión por unidad de
longitud en un flujo laminar, incompresible, viscoso a través de una
tubería recta de longitud l. Supongamos que las ecuaciones
diferenciales son muy complicadas para obtener esta relación
∴ Se busca una función que relacione las variables de interés.

¿Qué variables intervienen?:

caída de presión por unidad de longitud, caudal, viscosidad, diámetro
de la tubería y longitud de la tubería.
⎛ Δp ⎞
F⎜ , D , μ , V& ⎟ = 0
⎝ l ⎠
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 23

Teorema π de Buckingham
Aplicación del Teorema π de Buckingham:
Δp / l D V& μ
[M ] 1 0 0 1
⎛ Δp ⎞
F⎜ , D, μ , V& ⎟ = 0 [L ] −2 1 3 −1
⎝ l ⎠
[t ] −2 0 −1 −1
Matriz formada por el exponente a que está elevada cada dimensión fundamental en
cada una de las n variables dimensionales. Si se calcula el rango de esta matriz, (k)
éste vale 3, (en este caso k=j). Luego n-k =4-3=1 (número de grupos
adimensionales, πn-k ).
Ahora se eligen las k variables dimensionales que se repetirán en todos los grupos
adimensionales, esto es, k=3.
En mecánica de fluidos es típico usar c, ρ y D, pero podrían ser otras. En este caso
se usan Δp, D y V& .
Se forman así todos los grupos adimensionales πn-k, ciclando las tres variables
elegidas con las restantes variables dimensionales (1 grupo adimensional en total ).
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 24

12
 Teorema π de Buckingham
La expresión general de los coeficientes es:
⎛ Δp ⎞
a
⎛ Δp ⎞
F⎜ , D, μ ,V& ⎟ = 0 π i = ⎜ ⎟ DbV& c (variablesj )
⎝ l ⎠ ⎝ l ⎠
⎛ Δp ⎞
a
Luego, para el primer grupo, π1 → π1 = ⎜ ⎟ DbV& c μ
⎝ l ⎠
[ ]
M0 L0 t 0 = [ M L−2 t −2 ]a [L]b [L3 t -1 ]c ML-1t -1 = Ma+1 L−2a+b+3c−1 t −2a−c−1

μV&
Igualando exponentes queda: a = -1, b = -4, c = 1 ∴ π1 = Δp−1D−4V& 1μ =
Siendo la relación buscada: ΔpD4
μ V& ⎛ μ V&l ⎞ 1 μ V&l
G (π i ) = 0 → = π 1 → Δ p = G ⎜⎜ 4 ⎟⎟ =
Δp 4 ⎝ D ⎠ π1 D
4
D
l
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 25

Teorema π de Buckingham
Ejemplo 2:
Investigación dimensional de como varía la presión en un flujo
incompresible y viscoso a través de una tubería recta de longitud l.
Supongamos que las ecuaciones diferenciales son muy complicadas
para obtener esta relación
∴ Se escoge una opción EXPERIMENTAL

¿Qué variables intervienen?:

caída de presión, velocidad media, viscosidad, diámetro de la tubería,
longitud de la tubería, densidad del fluido y rugosidad de la tubería.

F (Δp, D , l , ε , c, ρ , μ , ) = 0
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 26

13
 Teorema π de Buckingham
Aplicación del Teorema π de Buckingham: F ( Δ p , c , μ , D , l , ρ , ε ) = 0
Δp D l ε c ρ μ
[M ] 1 0 0 0 0 1 1
[L ] −1 1 1 1 1 − 3 −1
[t ] −2 0 0 0 −1 0 −1
Matriz formado por el exponente a que está elevada cada dimensión fundamental. Si
se calcula el rango de esta matriz, (k) éste vale 3. Siendo n-k =7-3=4 ( número de
grupos adimensionales, πn-k ). Nótese que aquí k=j.
Luego se eligen 3 variables dimensionales que se repetirán en todos los grupos
adimensionales. En mecánica de fluidos es típico usar c, ρ y D, pero podrían ser otras.
Se forman así todos los grupos adimensionales πn-k, ciclando las tres variables
elegidas con las restantes variables dimensionales ( 4 grupos adimensionales en total,
agrupando las n-k=7-3=4 variables).
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 27

Teorema π de Buckingham

La expresión general de los coeficientes es: π1 = c a Db ρ c Δp

M 0 L0 t 0 = [ Lt −1 ]a [ L]b [ ML−3 ]c [ M L−1 t −2 ] = M c +1 La +b −3c −1 t − a −2

Δp
Igualando exponentes queda: a = -2, b = 0, c = -1 ∴ π1 = c−2 D0 ρ −1Δp =
ρc2
l ε μ
Luego: π2 = ca Db ρcl = ; π 3 = c a D b ρ cε = ; π 4 = c a D b ρ c μ =
D D cD ρ

Δp ⎛ l ε μ ⎞
Siendo la relación buscada: = G ⎜⎜ , , ⎟⎟
ρc ⎝ D D ρ cD
2
⎠
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 28

14
 Teorema π de Buckingham
Ejemplo 3:
La flexión δ de la punta de una viga en voladizo es función de la carga aplicada en la punta
P, la longitud de la viga l, el momento de inercia de la sección I y el módulo de elasticidad
del material E, esto es, δ =f(P, l, I, E). Reescriba esta función en forma adimensional .
δ P l I E
[L] [MLt-2] [L] [L4] [ML-1t-2]
Hay 5 variables dimensionales, n=5 y j=3 dimensiones primarias (MLT), puede
pensarse que los grupos π serían 2 ya que:
“...y k es como máximo el número de dimensiones fundamentales (j) para
expresar todas las xi”
δ P l I E
Luego k sería =3 y se tendrían 2 grupos π, [M ] 0 1 0 0 1
pero si se calcula el rango de la matriz es k=2.∴ [L ] 1 1 1 4 −1
n-k=3, deben obtenerse 3 grupos π. [t ] 0 −2 0 0 −2

Rango de la matriz, (reducción) k=2. Note que k es menor que el número de

dimensiones fundamentales j, (j=3).
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 29

Teorema π de Buckingham
Se ve que el rango de la matriz k es inferior al número de dimensiones del problema j. Esto
se debe a que [M] y [t] sólo aparecen en P y E y en ambos casos lo hacen de idéntica forma
[Mt-2] , indicando que es un caso especial en donde la reducción es k < j y no k=j.
Luego si n=5, xn=5 variables dimensionales, y j=3 dimensiones primarias (MLt), el rango de
la matriz es k=2.∴ n-k=3, debe obtenerse 3 grupos π. Para armar cada πi, si k=2., se
toman 2 xi (l y E) elevadas a las potencias a y b a determinar y una de las restantes xi que
quedan sin elevar a ninguna potencia.
La expresión para el 1er coeficiente es: π 1 = l a E b I ⎧ 0=b
b=0
π1 = l a E b I → [M 0 L0T 0 ] = [ L]a [ML−1T −2 ] [ L4 ] → ⎨0 = a − b + 4 ⇒
⎪ I
∴π1 = l −4 E 0 I =
b

⎪ 0 = −2b a = −4 l4
⎩
La expresión para el 2do coeficiente es: π 2 = l a E b P
⎧ 0 = b +1
b = −1
π 2 = l a E b P → [M 0 L0T 0 ] = [ L]a [ML−1T −2 ]
⎪ P
∴π 2 = l −2 E −1P = 2
b −2
[MLT ] → ⎨0 = a − b + 1 ⇒
⎪ 0 = −2b − 2 a = −2 El
⎩
La expresión para el 3er coeficiente es: π 3 = l E δ a b
⎧ 0=b
b=0 δ
π 3 = l a E bδ → [M 0 L0T 0 ] = [ L]a [ML−1T −2 ] [ L] → ⎨0 = a − b + 1 ⇒
⎪
∴π 3 = l −1E 0δ =
b

⎪ 0 = −2b a = −1 L
⎩
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 30

15
 Teorema π de Buckingham
Se obtiene finalmente δ ⎛ P I⎞
= f⎜ 2 , 4⎟
l ⎝ El l ⎠
Esta es una función complicada de 3 variables adimensionales y el análisis dimensional no
nos puede llevar más lejos.
Se puede “mejorar” esta ecuación conociendo más de la física del problema, ya que se
conoce que para deflexiones elásticas pequeñas δ es proporcional a la carga P, e
inversamente proporcional al momento de inercia I. Dado que P e I aparecen por separado
en la ecuación obtenida, significa que π3 es proporcional a π2, e inversamente proporcional
a π1, lo que permite escribir:
δ ⎛ P l4 ⎞ ⎛ Pl 3 ⎞
= (cte)⎜⎜ 2 ⎟⎟ → δ = (cte)⎜⎜ ⎟⎟
l ⎝ El I ⎠ ⎝ EI ⎠

El análisis dimensional, por si solo no es capaz de predecir esto. La teoría de la elasticidad

permite indicar que el valor de la cte es igual a 1/3.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 31

Teorema π de Buckingham
Ejemplo 4:

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 32

16
 Teorema π de Buckingham

n−k =
k

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 33

Teorema π de Buckingham

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 34

17
 Teorema π de Buckingham

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 35

Teorema π de Buckingham

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 36

18
 Método de Ipsen
El método basado en el teorema π, se le denomina a menudo el método de las
variables repetidas del análisis dimensional. Se seleccionan las variables
independientes, se añade una más y se obtiene un grupo π. El método tiene
algunas desventajas:
a) todos los grupos π contienen el mismo número de variables independientes
y podría perderse diversidad y con ello efectividad.
b) se debe comprobar que las variables repetidas no formen un grupo
adimensional entre ellas, lo cual en ocasiones es laborioso (relacionado con el
cálculo del menor complementario de la matriz variables dimensionales vs
dimensiones fundamentales).
Ipsen sugiere un método radicalmente diferente, un método paso a paso para
obtener todos los grupos π a la vez, sin necesidad de comprobaciones.
Simplemente se eliminan sucesivamente cada dimensión en la función deseada
mediante una multiplicación o una división.

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 37

Método de Ipsen
Ejemplo 5: Arrastre sobre cuerpo sumergido:
Se identifican los parámetros dimensionales relevantes y se escriben las
dimensiones de cada cantidad identificada.
F = f (L, c, ρ, μ)
[MLt ] =
-2
[ L] [ Lt −1 ] [ ML−3 ] [ M L−1 t −1 ]

Hay 3 dimensiones [MLt]. Se tratará de eliminar dimensiones de la ecuación,

para ello se toma una variable que tenga la dimensión [M], por ejemplo ρ y se
divide por ella a el resto de las variables que tengan dimensión de masa.
F ⎛ ρ μ⎞
= f ⎜⎜ L, c, , ⎟
ρ ⎝ ρ ρ ⎟⎠
⎡ MLt- 2 ⎤ −1 ⎡ ML ⎤ ⎡ M L t ⎤
−3 −1 −1

⎢ −3 ⎥
= [ L ] [ Lt ] ⎢ −3 ⎥ ⎢ −3 ⎥
⎣ ML ⎦ ⎣ ML ⎦ ⎣ ML ⎦
[ ] [
L4 t -2 = [ L] [ Lt −1 ] [1] L2 t −1 ]
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 38

19
 Método de Ipsen
Obsérvese que no se ha dividido ni L ni c, ya que no contienen [M]. La ecuación
puede parecer “extraña” pero no se ha perdido información.
F ⎛ μ⎞
∴ = f ⎜⎜ L, c, ⎟
ρ ⎝ ρ ⎟⎠
[L t ] =
4 -2
[ L] [ Lt −1 ] [ L2 t −1 ]

Hay ahora 2 dimensiones [Lt]. Se tratará de eliminar estas dimensiones de la

ecuación, para ello se toma una variable que tenga la dimensión [t], por ejemplo c y
se divide por ella el resto de las variables que tengan dimensión de tiempo. y a
las potencias adecuadas en cada término para eliminar la dimensión [t]
F ⎛ c μ⎞
= f ⎜⎜ L, , ⎟
ρc 2 ⎝ c ρc ⎟⎠
⎡ 4 -2 ⎤
L t ⎡ Lt −1 ⎤ ⎡ L2 t −1 ⎤
⎢ 2⎥
= [ L] ⎢ −1 ⎥ ⎢ −1 ⎥
⎢⎣ ( )
Lt −1 ⎥⎦ ⎣ Lt ⎦ ⎣ Lt ⎦
[ ]L2 = [ L] [1] [L]
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 39

Método de Ipsen
Se sigue sin perder información y sólo resta eliminar la dimensión L.
F ⎛ μ⎞
∴ 2 = f ⎜⎜ L, ⎟
ρc ⎝ ρ c ⎟⎠
[L ]2
= [ L] [ L]
Hay ahora 1 dimensión [L]. Se tratará de eliminar dimensiones de la ecuación, para
ello se toma la variable que tenga la dimensión [L], por ejemplo L y se divide por
ella el resto de las variables que tengan dimensión de longitud, y a las potencias
adecuadas en cada término para eliminar la dimensión [L]
F ⎛L μ ⎞ ⎫
= f ⎜⎜ , ⎟⎟ ⎪ ⎛ μ ⎞
ρc 2 L2 ⎝ L ρcL ⎠ ⎪
F
= f ⎜⎜ ⎟⎟
⎬⇒ ρc 2 L2 ⎝ ρcL ⎠
⎡ L2 ⎤ ⎡L⎤ ⎡ L⎤ ⎪
⎢ 2⎥ = ⎢⎣ L ⎥⎦ ⎢⎣ L ⎥⎦ ⎪
[1] = [1]
⎣L ⎦ ⎭
Se ha encontrado que el coeficiente de fuerza es sólo función del número de
Reynolds. No se han contado variables ni se ha identificado el rango de la matriz n-k.
Se operó con las variables adecuadas para eliminar las dimensiones primarias.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 40

20
 Método de Ipsen
Si se vuelve sobre el ejemplo 3, en donde el número de variables independientes
adimensionales es menor que el número de dimensiones primarias, pueden verse
las diferencias en la forma de operar.
δ = f ( P, L, I , E)
[L] = [ MLt-2 ] [ L] [ L4 ] [ M L−1 t −2 ]
Hay 3 dimensiones [MLt]. Se tratará de eliminar dimensiones de la ecuación, para
ello se toma una variable que tenga la dimensión [M], por ejemplo E y se divide
por ella el resto de las variables que tengan dimensión de masa.

⎛P E⎞
δ = f⎜ , L, I, ⎟
⎝E E⎠
⎡ MLt -2 ⎤ ⎡ M L−1 t −2 ⎤
[L] = ⎢ −1 − 2 ⎥
[ L] [ L4 ] ⎢ −1 − 2 ⎥
⎣ M L t ⎦ ⎣M L t ⎦
∴ [L] = [ ]
L [ L] [ L4 ] [1]
2

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 41

Método de Ipsen
Descartemos E y la dimensión [t] ha desaparecido junto con la dimensión [M],
Debe ahora eliminarse [L] dividiendo por ella el resto de las variables que
tengan dimensión de longitud y con las potencias adecuadas.

δ ⎛ P L I ⎞
⎛P ⎞⎫ = f⎜ 2 , , ⎟ δ ⎛ P I ⎞
δ = f ⎜ , L, I⎟⎪ L ⎝ EL L L4 ⎠ = f ⎜ 2 ,1, 4 ⎟
⎝E ⎠⎬ → ⇒L ⎝ EL L ⎠
⎡L⎤ ⎡ L2 ⎤ ⎡L⎤ ⎡L ⎤
4
2
[ ]
4 ⎪
[L] = L [ L] [ L ] ⎭ ⎢⎣ L ⎥⎦ ⎢ L2 ⎥
=
⎣ ⎦
⎢⎣ L ⎥⎦ ⎢ 4⎥
⎣L ⎦
[1] = [1] [1] [1]
δ ⎛ P I ⎞
∴ = f⎜ 2 , 4⎟
L ⎝ EL L ⎠
Se obtuvo idéntica respuesta que en el ejemplo 3. El hecho de que
desaparezcan [M] y [L] a la vez es la prueba de que en este caso hay sólo
dos variables independientes y no 3 como podría deducirse de la presencia de
[M], [L] y [t].
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 42

21
 Método de Ipsen
Ejemplo 4: El momento aerodinámico en el borde de ataque, MBA en un perfil
supersónico (Ma>1) es función de la longitud de la cuerda C, el ángulo de ataque α y
diversos parámetros del fluido (aire) a saber: la velocidad del flujo, c∞, la velocidad del
sonido a=cw,nc , y la relación de calores específicos γ. La experiencia indica que el
efecto viscoso (μ) es débil (capa límite delgada). Se pide reescribir esta relación
funcional en términos adimensionales utilizando el procedimiento de Ipsen.
M BA = f (C , α , c∞ , ρ, a, γ )
MBA [ML t ] = [L] [1] [ Lt
2 −2 -1 −3
] [ M L ][ Lt -1 ] [1]
c∞
α
Las variables α y γ ya son adimensionales y por
lo tanto forman dos de los grupos π buscados.

Hay 3 dimensiones [MLt]. Se tratará de eliminar dimensiones de la ecuación, para ello

se toma una variable que tenga la dimensión [M], por ejemplo ρ y se divide por ella el
resto de las variables que tengan dimensión de masa.

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 43

Método de Ipsen
M BA ⎛ ρ ⎞
∴ = f ⎜⎜ C , α , c∞ , , a, γ ⎟⎟
ρ ⎝ ρ ⎠
⎡ ML2 t −2 ⎤ ⎡ M L−3 ⎤ -1
⎢ −3 ⎥
= [ L] [1] [ Lt -1 ] ⎢ −3 ⎥
[ Lt ] [1]
⎣ M L ⎦ ⎣M L ⎦
[ ]
L5 t −2 = [ L] [1] [ Lt -1 ] [1] [ Lt -1 ][1]
Se descarta ρ quedando:
= f (C , α , c∞ , a, γ )
M BA
∴
ρ
[L t ] =
5 −2
[ L] [1] [ Lt -1 ] [ Lt -1 ] [1]
Se elimina ahora [t] dividiendo por a a las potencias adecuadas en cada término :
M BA ⎛ c a ⎞⎫
∴ = f ⎜ C, α , ∞ , , γ ⎟⎪
ρa 2 ⎝ ⎠⎪ M BA ⎛ c ⎞
a a = f ⎜ C, α , ∞ , γ ⎟
⎡ 5 −2 ⎤ ⎬→ ρa 2
⎝ a ⎠
⎡ Lt ⎤ ⎡ Lt ⎤
[ ]
-1 -1
= [ L] [1] ⎢ -1 ⎥ ⎢ -1 ⎥ [1] ⎪
Lt
⎢ 2⎥
L3 = [ L] [1] [1] [1]
⎢⎣ ( )
Lt −1 ⎥⎦ ⎣ Lt ⎦ ⎣ Lt ⎦ ⎪
⎭
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 44

22
 Método de Ipsen
Se elimina por último [L] dividiendo por C a las potencias adecuadas en cada término:

M BA ⎛C c ⎞⎫
= f ⎜ , α , ∞ , γ ⎟⎪ M BA ⎛ c ⎞
ρa 2C 3 ⎝C a ⎠⎪ = f ⎜ 1, α , ∞ , γ ⎟
∴ ⎬ → ρa 2C 3 ⎝ a ⎠
⎡L ⎤
3
⎡L⎤
⎢ 3⎥ = ⎢ ⎥ [1] [1] [1] ⎪ [1] = [1] [1] [1] [1]
⎣L ⎦ ⎣L⎦ ⎪⎭

⎛ ⎞
= f ⎜ α , ∞ , γ ⎟ ⇒ CM = f ( α , Ma, γ )
M BA c
∴
ρa 2C 3 ⎝ a ⎠

Se verá luego (flujo compresible) que por evidencia experimental c∞ tiene más
influencia en la variación de CM que a, por lo que en aerodinámica se define el
coeficiente de momento CM, como:
M BA
CM =
ρc∞2 C 3
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 45

Estudio de modelos
Se ha estudiado la homogeneidad dimensional y el Teorema π para,
llevar una relación físicamente homogénea a su forma adimensional usando
productos de potencias.
Aunque esta transformación sea matemáticamente correcta hay ciertas
dificultades ingenieriles previas que necesitan analizarse.
Hemos dado por cierto que es posible especificar todas las variables que
intervienen en el proceso. El primer análisis que debe hacerse es para
contestar las preguntas: ¿se han retenido todas?, ¿todas las variables
retenidas son importantes?,... cada variable retenida aumenta el costo del
análisis.
Una vez seleccionadas las variables y realizado el análisis dimensional,
debe analizarse la viabilidad de que los experimentos en un modelo a
escala sirvan para extrapolarlos a un prototipo. Debe garantizarse la
semejanza entre el modelo a ensayar y el prototipo a diseñar....¿Cómo?
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 46

23
 Estudio de modelos
Con suficientes ensayos, los datos obtenidos del modelo proporcionarán
la función adimensional buscada:
π1=f(π2, π3,....,πn-k)
Obtenida la ecuación y plasmada en tablas, gráficas o de forma analítica, se
está en posición de analizar la semejanza entre modelo y prototipo... para
responder el ¿cómo?.
Una definición formal del ¿cómo? podría ser:
“Las condiciones de flujo para un modelo de ensayo son
completamente semejantes a las de un prototipo si coinciden todos
los parámetros adimensionales correspondientes entre modelo y
prototipo”... y ver si se cumple o no...
Por ser difícil de obtener la semejanza completa, se habla de tipos
particulares de semejanza siendo las más comunes la geométrica, la
cinemática, la dinámica y la térmica.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 47

Estudio de modelos

Tipos de semejanza:
( x , y , z ) modelo homólogo ( x , y , z ) prototipo
• Geométrica
Lr = λ = L p / Lm ; λ : relación de escala lineal
r r r
cm wm um
• Cinemática Cr = r = r = r
cp wp up
t r = Lr / C r ; Ar = Lr / t r
2
14243 14243
escala de tiempos escala de aceleraciones
r
Fm
• Dinámica Fr = r
Fp
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 48

24
 Estudio de modelos
r
Semejanza dinámica: F
Fr = rm
Fp
Las fuerzas que pueden actuar en el modelo (Fm) y prototipo (Fp) son:
• fuerzas volumétricas FB : entre ellas la gravitatoria Fg =mg=ML/t2, las
fuerzas de inercia, FI=ma= ML/t2=ρc2L2 y las elásticas, FE=EvA=EvL2,
• fuerzas superficiales FS: debidas a la presión, Fp =ΔpA= ΔpL2, a la
viscosidad, Fv =μ(∂c/ ∂y )A=μ (c/L)L2 y a las tensiones superficiales si el
elemento fluido está en una interfase líquido-gas, FTS= σL.

En conservación de cantidad de movimiento puede tenerse:

ΣF =ma ⇒ ΣF -ma =0; ma=FI
∴ FB + Fs = Fg + FI + FE+ Fp + Fv + FTS = 0
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 49

Estudio de modelos
Semejanza dinámica:
Si se consideran las fuerzas de inercia, de presión, viscosas y gravitatorias semejantes
para modelo y prototipo puede escribirse: r r r r
FI ,m Fp,m Fv,m Fg ,m
r = r = r = r
FI , p Fp, p Fv, p Fg, p
r r r r r r
Pudiendo reordenarse como:
FI ,m FI , p FI ,m FI , p FI ,m FI , p
r = r ; r = r ; r = r
Fp , m Fp , p Fv ,m Fv , p Fg ,m Fg , p

Si las fuerzas ya señaladas fueran las únicas presentes se tiene:

( )
r r r r
FI = f Fp , Fv , Fg
Se verá el significado de esta relación cuando se definan los parámetros
adimensionales correspondientes
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 50

25
 Estudio de modelos
Semejanza analógica:
Puede darse el caso de que fenómenos de naturaleza intrínsecamente
diferente, respondan a una ecuación diferencial o integral de idéntica forma.
Ejemplo: Para la relación presión-caudal, se dedujo del análisis
dimensional una relación del tipo:
μV& 1 μV&l μl &
= π 1 → Δp = ⇒ Δp = V → Δp = RhidrV&
Δp 4 π D 4
π D 4
D 1 1
l
Si se recuerda la ley de Ohm: ΔV = Relect I
Puede verse que ambos fenómenos responden a una ley matemática similar.
Luego, se puede montar un circuito analógico en donde se definen
resistencias eléctricas equivalentes a las hidráulicas y se miden tensiones y
corrientes en vez de caudales y presiones.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 51

Números de Euler, Froude, Reynolds, Weber y Mach

Se ha observado a partir del estudio de numerosos

problemas de la termo-fluido-dinámica, que para poder
hacer comparaciones entre experimentos que resulten
similares, debe haber como mínimo semejanza dinámica
entre estos experimentos.
Se han definido así una serie de parámetros adimensionales
muy útiles y que comparan distintos aspectos de esta
semejanza.

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 52

26
 Números de Euler, Froude, Reynolds, Weber y Mach
Número de Reynolds
Re = Fuerzas inercia / Fuerzas viscosas.
L
r ρL3
T 2 = ρc L = ρcL
2 2
FI ma
Re = r = =
Fv ⎛ dc ⎞ ⎛c⎞ μcL μ
μ ⎜⎜ ⎟⎟ A μ ⎜ ⎟ L2
⎝ dy ⎠ ⎝L⎠
Semejanza entre modelo y prototipo Re m = Re p
Prototipo y modelo de turbomáquina con una Lr=λ=Lp/Lm=5, siendo np=1.000 y
considerando el mismo fluido en ambas máquinas Para que se cumpla que Rem=Rep ¿cuál
deberá ser la velocidad de giro del modelo? πD n πD n p p m m
Dp Dm
Re m = Re p ⇒ D p u p / ν p = D m u m / ν m ⇒ 60 = 60
νp νm
D p D pnp Dm Dm nm ν mD np 2
D 2

∴ = → nm = = np = 1.000 × 5 2 ⇒ n m = 25.000 rpm !

p p

νp νm ν p Dm
2
D 2
m
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 53

Números de Euler, Froude, Reynolds, Weber y Mach

El parámetro que cuantifica el comportamiento más o menos viscoso es un fluido
es el número adimensional de Reynolds
ρ : densidad, [kg/m3]
c : velocidad característica, [m/s] (en turbomáquinas
ρ c L c L L : longitud característica, [m]
c=Ω r)

Re = =
μ ν μν :: viscosidad
viscosidad dinámica, [kg/(m s)]
cinemática, [m /s] 2

L: para flujo en tuberías: diámetro (tubería).

L: para flujo sobre placas planas: longitud desde el borde de ataque.
L: para flujos sobre cuerpos sumergidos: longitud característica cuerpo, ej.:
diámetro (cilindros, esferas), cuerda (perfiles alares).
L: para capas límites: o bien la longitud a partir del borde de ataque del cuerpo o el
espesor de capa límite.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 54

27
 Números de Euler, Froude, Reynolds, Weber y Mach
Número de Froude
Fr = Fuerzas inercia / Fuerzas gravitacionales.
r L
FI ma 2 L L L c2 c
Fr = r = = T = 2
= 2
× = → Fr =
Fg mg g gT gT L gL gL

Semejanza entre modelo y prototipo Fr m = Fr p

Prototipo y modelo de barco con una Lr= λ = Lp/Lm=100, velocidad del prototipo,
cp de 20 nudos (1,8km/h=1nudo) y considerando el mismo fluido en ambas
máquinas, ¿cuál será la velocidad a la que debe mover el modelo para que se
cumpla que Frm=Frp?
cm cp Lm 1 1
Frm = Frp ⇒ = ⇒ cm = c p = cp = 20 = 2 nudos
g m Lm g p Lp Lp λ 100
∴ c m = 2 nudos
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 55

Números de Euler, Froude, Reynolds, Weber y Mach

Número de Euler
Eu = Fuerzas inercia / Fuerzas de presión.
r L
ρL3 2
FI ma T = ρL ∴ Eu = ρc → Eu = ρc =
2 2 2
c
Eu = r = =
Fp pA ΔpL2
ΔpT 2
Δp 2 Δp 2 Δp/ρ
Nota: Se agrega la raíz cuadrada y el 2 para que quede como relación de energía de presión y cinética. Además
la radicación sirve para dejar a la velocidad con exponente =1, como el Re
Semejanza entre modelo y prototipo Eu m = Eu p
Prototipo y modelo de barco con una Lr= λ = Lp/Lm=100, velocidad del prototipo 20
nudos (1,8km/h=1nudo). Si Fm=200gr y considerando el mismo fluido en ambas
máquinas, ¿qué fuerza se ejerce sobre el prototipo para que se cumpla que Eum=Eup?
r r 2
c 2p c m2 c 2p Fp Fm c p
Eu m = Eu p ⇒ = → Δp p = Δpm 2 → =
Δp p / ρ p Δpm / ρ m cm Ap Am c m2
2
r r L c 2p 2
⎛ 20 nudos ⎞
∴ F p = Fm = 0,2 kgf × 100 2
× ⎜ ⎟ = 2 × 10 4 × 10 2 kgf = 2 × 10 6 kgf
p

L c m2 ⎝ 2 nudos ⎠
2
m

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 56

28
 Números de Euler, Froude, Reynolds, Weber y Mach
Número de Weber
We = Fuerzas inercia / Fuerzas tensión superficial.
r L
FI ma ρL3 2 ρL3 ρL2 L ρc 2 L ρc 2 L c
We = r = = T = = = → We = =
F TS σL σ L σT 2
σT 2
σ σ σ / (ρL )

Semejanza entre modelo y prototipo We m = We p

Prototipo canal costero p/mareas de Lp=10m, con μ =1,07×10-3 kg/(ms), ρ =1.025kg/m3
σ=7,28 ×10-2 N/m, y se escala a Lr=Lp/Lm=1/100. El modelo se hace funcionar con agua
dulce, μ =1,0×10-3 kg/(ms), ρ=998kg/m3, σ=7,28 ×10-2 N/m. Para una velocidad de la
corriente marítima de 0,5m/s cual será la velocidad del modelo?.
ρ pu 2p Lp ρmum2 Lm ρ pσ m Lp 2
⎛ m ⎞ 998 m
Si se cumple Wem=Wep ⇒ = ⇒ um = u2p = ⎜ 0,5 ⎟ × ×100 = 4,93
σp σm ρmσ p Lm ⎝ s ⎠ 1.025 s
Semejanza total: Debe haber también igualdad de Re y Fr. Se comprueba que si
Rem=Rep→cm=48m/s,∴Wem≠Wep y no hay semejanza total. Aquí, la semejanza más importante a
cumplir es la de Fr. Un análisis para estas escalas muestra que es imposible cumplir similitud a la vez
de Re y Fr. Se empieza a hablar de simulación parcial (se tiene en cuenta efectos menores por
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 57
posteriores correcciones de los resultados). Esto se verá más adelante...

Números de Euler, Froude, Reynolds, Weber y Mach

Número de Mach
Ma = Fuerzas inercia / Fuerzas elásticas.

r ρL3 2
L
T = ρL = ρc → Ma = ρc =
2 2 2
FI ma c c
Ma = r = = =
FE EV A EV L2
E vT 2
E v E v E v / ρ c sonido

Semejanza entre modelo y prototipo Ma m = Ma p

Modelo con una Lr=λ = Lp/Lm=1, se desea conocer la velocidad de la corriente de ensayo
en un túnel de viento si la velocidad de diseño del avión para una altura de vuelo de 12km
con Tp=216,7K es de 650km/h, csonido=(γRT)0,5, γ=Cp/Cv
Para que se cumpla que Mam=Map

up um c km γ R ( 293K ) km
Ma m = Ma p ⇒ = ⇒ u m = u p s , m = 650 = 756
cs, p c s ,m cs, p h γ R ( 216,7K ) h
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 58

29
 Números de Euler, Froude, Reynolds, Weber y Mach
Número de Strohual
Sr =Tiempo característico del flujo /Período de la oscilación
Se supone una velocidad de entrada al vc que cambia en el tiempo
c ⎛ ωL * ⎞
c = c∞ cos ωt → Adimensionalizándola → c* = = cos⎜⎜ t ⎟⎟
c∞ ⎝ c∞ ⎠
ωL ωL fL
∴ Sr = ; ω = 2πf ⇒ Sr = =
c∞ 2πc∞ c∞
Semejanza entre modelo y prototipo Sr m = Sr p
¿Cómo afecta el viento a una chimenea de 2,75m de diámetro, al asumir que la frecuencia
de vibración natural del modelo y prototipo más importante es de 7 rad/s? ¿Qué velocidad
del viento induce vórtices con desprendimiento de 7 rad/s en el prototipo, si se define
Lr=Lp/Lm=3?, Otros datos: se considera el mismo fluido en ambos casos y se comprueba
el comienzo de las oscilaciones en el prototipo para una velocidad de 18,33m/s.
f p Lp f m Lm f L m m
Para que se cumpla que Srm=Srp ⇒ = ⇒ cm = c p m m = 18,33 × 3 = 54,99
cp cm f p Lp s s
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 59

Números de Euler, Froude, Reynolds, Weber y Mach

Teniendo en cuenta los distintos tipos de fuerzas que aparecen, se definen
ciertos parámetros adimensionales:
• Reynolds (Re = ρc L/µ ): relaciona fuerzas de inercia y viscosas;
• Weber (We = ρc2 L/σ ): relaciona fuerzas de inercia y tensiones superficiales
(hay superficie libre);

• Mach (Ma =c/csonido =(Δρ /ρ0)0.5 ): relaciona las fuerzas elásticas con las de
inercia; indica la compresibilidad de un fluido;

• Froude (Fr=c2/(gl) ): indica la relación entre las fuerzas gravitatorias y las de

inercia (ondas en superficies libres);
• Euler (Eu=ΔP/(ρc2) ) : relaciona fuerzas de presión y de inercia;

• Strohual ( Sr =fL/c = L/(c0t0) ): Flujo oscilatorio, relaciona la velocidad de la

oscilación con la velocidad media, f es la frecuencia de oscilación.
Ejemplo análisis N-S
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 60

30
 Números de Prandtl, Nusselt y Peclet
Asociados al estudio de transferencia de calor se han definido ciertos parámetros
adimensionales y se da su interpretación como el cociente de dos magnitudes:
• Prandtl (Pr = ν/α =μCp/kT ): relaciona difusión viscosa y difusión térmica;
α: Difusividad térmica, [m2/s]
[ M ][ L1 ][ t 1 ][ L2 ][ t -2 ][ T -1 ]
[Pr] = [ Adim] = Cp: Calor específico a presión cte. [J/(kgK)],
[ M ][ L][ t - 3 ][ T -1 ]
kT: Conductividad térmica, [w/(mK)]
• Nusselt (Nu = Lh/kT ): relaciona transf. calor por convección/transf. calor por
conducción molecular); 647 h
48
[ M ][ t ][ T -1 ][ L]
-3
h: Coeficiente de transferencia de calor, [w/(m2K)]
[Nu] = [ Adim] = -3 -1
[ M ][L][ t ][ T ] kT: Conductividad térmica, [w/(mK)]
144244 3
kT
• Peclet (Pe =Re×Pr =L c ρ Cp /kT ): relaciona transf. calor por convección con la transf.
calor por conducción;
[ L][ Lt -1 ][ ML-3 ][L2 t -2 T -1 ]
[Pe] = [ Adim] =
[1
M4 ][L]
42[ t -4
3
][ 4-1
T3 ]
kT
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez
Ejemplo análisis N-S 61

Números de Prandtl, Nusselt y Peclet

Asociados al estudio de transferencia de calor se han definido ciertos parámetros
adimensionales y se da su interpretación como el cociente de dos magnitudes:
• Grashof (Gr = βΔTgL3ρ2/μ2): relaciona flotación (buoyancy) y fuerzas inerciales;
β: Coeficiente de expansión volumétrica [(C,K)-1],
3 -3 2 -2
[L ][ML] [Lt ][T ][T][T -1 β =[ρ -1 (∂ρ/∂T)p=cte]
[Gr]= [Adim] =
[ML−1t-1]2 ΔT: Diferencia de temperaturas entre las “dos
fuentes” donde se establece el flujo. [C,K],

• Rayleigh (Ra = L3ρ2gβCpΔT /μkT = gΔρL3/νkTρ0= GrPr ): relaciona flotación

(buoyancy) y fuerzas inerciales;

[L3 ][ML-3]2[Lt-2 ][T-1][L2t -2T-1][T]

[Ra]= [Adim] = Cp: Calor específico a presión cte. [J/(kgK)],
[ML−1t -1][MLt−3T-1]
kT: Conductividad térmica, [w/(mK)]

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez

Ejemplo análisis N-S 62

31
 Números de Euler, Froude, Reynolds, Weber y Mach
Relaciones entre los números de Eu, Re y Fr
Si se consideran las fuerzas de inercia, de presión, viscosas y gravitatorias
semejantes para modelo y prototipo puede escribirse:
r r r r
FI ,m Fp,m Fv,m Fg ,m
r = r = r = r
FI , p Fp, p Fv, p Fg, p
Pudiendo reordenarse como:
r r r r r r
FI ,m FI , p FI ,m FI , p FI ,m FI , p
r = r ; r = r ; r = r ∴ Eu m = Eu p ; Re m = Re p ; Frm = Frp
F p , m Fp , p Fv ,m Fv , p Fg ,m Fg , p

Si las fuerzas ya señaladas fueran las únicas presentes se tiene:

( )
r r r r
FI = f Fp , Fv , Fg → Eu = f (Re, Fr)

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 63

Adimensionalización de las ecuaciones básicas

Se utilizó el Teorema π para analizar diferentes problemas, determinando los
parámetros adimensionales que gobiernan cada caso. Una técnica alternativa y
muy eficaz consiste en utilizar directamente las versiones diferenciales de las
ecuaciones básicas de conservación.

Aunque en los casos generales estas ecuaciones son de difícil resolución, revelan
los parámetros básicos dando pistas de cuando pueden despreciarse unos
términos respecto de otros.
Las condiciones de contorno e iniciales también deberán adimensionalizarse.
Las ecuaciones a adimensionalizar son las de conservación: de masa, de cantidad
de movimiento y de la energía.
∂ρ r
+ ∇ ⋅ ρc = 0
∂t
r r
⎛ ∂c r r⎞ r r
= ρ ⎜ + (c ⋅ ∇ )c ⎟ = −∇p + ρg + μ∇ 2 c
Dc
ρ
Dt ⎝ ∂t ⎠
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 64

32
 Adimensionalización de las ecuaciones básicas
Se redefinirán las variables dependientes e independientes para que sean
adimensionales, a través de dividirlas por constantes de referencia apropiadas para
el flujo.
Algunos ejemplos pueden ser: una velocidad de referencia del flujo, c0, o una
longitud de referencia L0. Es esencial usar constantes y no mezclar variables entre
ellas, esto es, hacer x/L0 o y/L0 y no x/y como nuevas variables.
Las lista de abajo es una lista completa de variables sugeridas:
1. Longitud de referencia, L0.
2. Velocidad de referencia del flujo, c0, (velocidad de corriente libre) o μ0/(ρ0L0)
si no hay una corriente libre (ej.: convección natural)
3. Tiempo de residencia t0= L0/c0, (en flujos estacionarios no hay tiempo
característicos, de allí esta definición).
4. otras propiedades de referencia de corriente libre: p0, μ0, ρ0, T0, k0 .
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 65

Adimensionalización de las ecuaciones básicas

Las variables adimensionales que se obtiene a partir de la división de las variables
dimensionales por las de la lista sugerida es:
r
xi p − p0 p − p∞ r* c T − T0 tc
xi = ,
*
p =
*
= , c = , T* = , t * = 0 = ft,
L0 ρ0c02
p0 − p∞ c0 Tw − T0 L0
r
μ ρ k C r g
μ * = , ρ * = , k * = , ∇* = L0∇, C*p = p , g * =
μ0 ρ0 k0 C p0 g
En lo que sigue, se utilizarán algunos de estos parámetros. Nótese que:
El orden de magnitud de las variables adimensionales es la unidad al haberles eliminado sus
dimensiones con longitudes, velocidades y otras variables de referencia, que sean
características del campo del flujo, por lo tanto t*~1, c*~1, x*~1, etc.
La definición de velocidad adimensional es vectorial y vale para las componentes también.
la presión adimensional se ha definido en términos de una diferencia de presiones ya que en el
caso de flujos incompresibles “el campo de velocidades no se afecta por la magnitud
absoluta de la presión, p, sino sólo por diferencias de presión”, como se vio en los
ejemplos de resolución de casos simples de las ecuaciones de Navier Stokes.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 66

33
 Adimensionalización de las ecuaciones básicas
Operando con estas variables en la ecuación de Continuidad se tiene:
∂ρ r ∂ρ 0 ρ * 1 * * r* c0 ρ 0 ∂ρ * ρ 0 c0 * * r *
+ ∇ ⋅ ρc = 0 → + ∇ ⋅ ρ ρ c c = 0 ⇒ + ∇ ⋅ρ c = 0
∂t ⎛ L t* ⎞ L L0 ∂t *
0 0
L0
∂⎜⎜ 0 ⎟⎟
⎝ c0 ⎠
c0 ρ 0 ⎛ ∂ρ * r ⎞ ⎛ ∂ρ * r ⎞
∴ ⎜⎜ * + ∇* ⋅ ρ *c * ⎟⎟ = 0 ⇒ ⎜⎜ * + ∇* ⋅ ρ *c * ⎟⎟ = 0
L0 ⎝ ∂t ⎠ ⎝ ∂t ⎠
Sólo hay un cambio de variables y no aparecen parámetros adimensionales. La
continuidad no es directamente influenciable por cualquiera de los parámetros del flujo.
Operando con estas variables en la ecuación de N-S se tiene (conviene reemplazar
aquí el t* por tf ). r r
⎛ ∂c r r⎞ r r
= ρ ⎜ + (c ⋅ ∇ )c ⎟ = −∇p + ρg + μ∇ 2 c
Dc
ρ
⎛ Dt ⎞⎝ ∂t ⎠
⎜ r ⎟
⎜ ∂c0 c * ⎛ r * 1 * ⎞ r * ⎟
ρ0 ρ * ⎜ + ⎜
⎛ t * ⎞ ⎜⎝
c0 c ⋅ ∇ ⎟
L0 ⎟⎠
c0 c
1 * 2 *
( * r* *
)
1 *2 r *
⎟ = − L ∇ ρ 0 c0 p + p0 + ρ 0 ρ gg + μ μ 0 L2 ∇ c0 c
⎜ ∂⎜ ⎟ ⎟ 0 0
⎜ ⎜f⎟ ⎟
⎝ ⎝ ⎠ ⎠
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 67

Adimensionalización de las ecuaciones básicas

⎛ ⎞
⎜ r* ⎟
* ⎜ ∂c0 c ⎛ r* 1 * ⎞ r* ⎟
ρ0 ρ ⎜
⎛ ⎞*
1
( r
)
1 2 r
+ ⎜⎜ c0c ⋅ ∇ ⎟⎟c0 c ⎟ = − ∇* ρ 0 c02 p * + p0 + ρ 0 ρ * gg * + μ * μ 0 2 ∇* c0 c *
⎜ ∂⎜ t ⎟ ⎝
L0 ⎠ L0 L0
⎟
⎜ ⎜f⎟ ⎟
⎝ ⎝ ⎠ ⎠
r ⎡ 78 ⎤
= 0 , p0 = cte
6
⎛ ∂c * c02 r * * r * ⎞ ⎢∇ ρ c p + ∇* ( p )⎥ + ρ gρ * gr * + μ 0 c0 μ *∇* 2 cr *
⎜
ρ 0 ρ ⎜ fc0 * +
*

∂
(
c ⋅ ∇ c ⎟⎟ = −)1
⎢
*
(
2
0 0
*
) 0
⎥ 0
L20
⎝ t L ⎠ L
⎣ ⎦
0 0

r
⎛ ∂c * ρ ρ *c 2 r r ⎞ ρ c2 r μc
( ) ( ) 2r
⎜⎜ ρ 0 fc0 ρ * * + 0 0 c * ⋅ ∇* c * ⎟⎟ = − 0 0 ∇* p* + ρ 0 gρ * g * + 02 0 μ *∇* c *
∂t
⎝ L0 ⎠ L0 L0
Dividiendo por las constantes del término convectivo se tiene:
⎛ ⎞ ρ 0c02 μ 0 c0
⎜ r* ⎟
ρ ρ ∂ ρ
( ) ( )
*
⎜ 0 fc c r r *⎟ L g r L20 2r
0
⎜ ρ ρ *c 2 ∂t * + c *
⋅ ∇ *
c ⎟ = − 0
∇ *
p *
+ 0
ρ * *
g + μ *∇* c *
ρ 0 ρ c0
* 2
ρ 0 ρ c0
* 2
ρ 0 ρ c0
* 2

⎜⎜ L ⎟⎟
0 0

⎝ 0 ⎠ L0 L0 L0
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 68

34
 Adimensionalización de las ecuaciones básicas
} p − p0 μ0
r
⎛ L0 f ∂c * r * * r * ⎞ ρ 0c02 * * ⎛ L0 g ⎞ r *
⎜⎜ (
+ c ⋅ ∇ c ⎟⎟ = − ) ∇ p + ⎜⎜ 2 ⎟⎟ g + ( )
L0
ρ 0 ρ c0
2r
μ *∇* c *
⎝ c0 ∂t ρ 0 ρ c0
* * 2 *
⎠ { ⎝ c0 ⎠
ρ
r
⎛ L0 f ∂c * r * * r * ⎞ p − p0 * * ⎛ L0 g ⎞ r * μ
⎜⎜
∂ *
(
+ c ⋅ ∇ c ⎟⎟ = − )ρ 2
∇ p + ⎜⎜ 2 ⎟⎟ g +
ρ
( )
2r
∇* c *
⎝ 0
c t ⎠ c0 ⎝ 0 ⎠
c L c
0 0

r
∂c * r * * r *
( 1 r
)
Sr * + c ⋅ ∇ c = − Eu ∇* p * + g * +
∂t Fr
1 *2 r *
Re
∇ c ( )
Nótese que:
Al las variables adimensionales t*~1, c*~1, x*~1, etc.(orden 1), en consecuencia los
r
( r
)
términos c * ⋅ ∇* c * , ∇* p * serán de orden 1 y a su vez son del mismo orden, por lo que la
magnitud relativa de unos términos de la ecuación sobre otros dependerá de la magnitud de
los parámetros adimensionales Re, Fr, Eu y Sr.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 69

Adimensionalización de las ecuaciones básicas

Un proceso similar se puede hacer con la ecuación de la energía, para flujo
incompresible, en donde Φ es la función disipación y se tiene:

⎛ ∂T r ⎞ Dp
+ (c ⋅ ∇ )T ⎟ = + ∇ ⋅ (k∇T ) + Φ
DT
ρC p = ρ⎜
Dt ⎝ ∂t ⎠ Dt
⎡ ⎛ ∂u ⎞2 ⎛ ∂v ⎞2 ⎛ ∂w ⎞2 ⎛ ∂v ∂u ⎞2 ⎛ ∂w ∂v ⎞2 ⎛ ∂u ∂w ⎞2 ⎤ ⎛ ∂u ∂v ∂w ⎞2 2
Φ = μ ⎢2⎜ ⎟ + 2⎜⎜ ⎟⎟ + 2⎜ ⎟ + ⎜⎜ + ⎟⎟ + ⎜⎜ + ⎟⎟ + ⎜ + ⎟ ⎥ + λ⎜⎜ + + ⎟⎟ ; λ = − μ
⎢⎣ ⎝ ∂x ⎠ ⎝ ∂y ⎠ ⎝ ∂z ⎠ ⎝ ∂x ∂y ⎠ ⎝ ∂y ∂z ⎠ ⎝ ∂z ∂x ⎠ ⎥⎦ ⎝ ∂x ∂y ∂z ⎠ 3

ρ *C *p
DT *
Dt *
= Ec
Dp
+
1
Dt Re Pr
∇* ⋅ k *∇*T * +
Ec
Re
Φ ( )
• Prandtl (Pr = ν/α =μCp/kT ): relaciona difusión viscosa y difusión térmica;
• Eckert (Ec = c02/ [Cp0(Tw-T0)] ): relaciona energía cinética y entalpía

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 70

35
 Modelado físico

El análisis dimensional aplicado a turbomáquinas tiene dos

importantes aplicaciones:

• La predicción del comportamiento de un prototipo a partir de

realizar estudios sobre modelos, esto es, una turbomáquina a
escala, según conveniencia (semejanza).

• La determinación del tipo más adecuado de turbomáquina,

sobre la base de máxima eficiencia, para un rango específico de
alturas, velocidades y caudales.

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 71

Modelado físico
Modelado:

• Si logro agrupar las variables independientes en grupos, sólo debo

cambiar el valor del grupo adimensional (por ejemplo, para cambiar el
número de Reynolds basta cambiar la velocidad y no simultáneamente
la longitud característica o la viscosidad del fluido).

• El análisis dimensional permite sugerir que variables pueden tener

poca importancia y sólo basta hacer unos pocos ensayos que
confirmen estas hipótesis.

• El análisis dimensional es el que en realidad proporciona las leyes

de escala (semejanza).

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 72

36
 Modelado físico
Modelado:
Se han escogido todos los grupos dimensionales satisfactoriamente pero
subsisten la siguientes dificultades:
¿Es posible especificar todas las variables que intervienen?,
¿qué efectos son los importantes? (temperatura, viscosidad, rugosidad,
etc...).

La experiencia y el buen juicio indican cuales son los parámetros a

retener. Un análisis estimativo de los valores de estos grupos
adimensionales, permiten cuantificar su importancia.

Cada grupo adimensional retenido aumenta el esfuerzo

necesario para obtener resultado y por ende, los costos.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 73

Modelado físico
Modelado:
El ensayo experimental puede realizarse sobre un prototipo o su
modelo a escala. El análisis estimativo implica buscar si hay semejanza
entre el modelo estudiado y el prototipo que se desea diseñar.
3-channel StreamLine with
Tri-axial wire probe

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 74

37
 Modelado físico
Modelado:

• Si los parámetros adimensionales coinciden para modelo y

prototipo podemos asegurar que hay semejanza completa entre
modelo y prototipo.

• Esto implicará hablar de semejanza geométrica, cinemática,

dinámica y térmica.

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 75

Modelado físico
Modelado: A tener en cuenta cuando se comparan modelos y prototipos:
• La velocidad en el modelo no debe ser tan baja que el flujo se
convierta en laminar o tan alta que los efectos de compresibilidad
sean notables.
• Que en el modelo la tensión superficial se vuelva importante.
• La rugosidad de un modelo respecto del prototipo debería reducirse
en la misma proporción de tal manera que se conservase la
semejanza geométrica.
Cuando no mantener alguna de estas relaciones es inevitable, se
habla de simulación parcial, esto es, tener en cuenta las
distorsiones de los resultados obtenidos en el modelo a la hora
de extrapolarlos al prototipo.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 76

38
 Modelado físico
Modelado: A tener en cuenta cuando se comparan modelos y prototipos:
Puede ocurrir que no se puedan igualar todos los parámetros adimensionales
simultáneamente, semejanza dinámica, entre modelo y prototipo ya que implicaría el
uso de: para el caso de igualar Rem=Rep y Frm=Frp y para una cierta escala
geométrica entre modelo y prototipo, tener un fluido con propiedades que no
existen (e.g.: viscosidad).
Lo mismo ocurriría en semejanza dinámica, en el caso de igualar Mam=Map y
Rem=Rep simultáneamente.
Ejemplo: Si para un prototipo de escala 25 veces el tamaño del modelo se usa una
presión de 70kPa y una T=268K, típica a 3km sobre el nivel del mar, cual es la
presión a la que se mantendrá la atmósfera en el dispositivo de ensayo para
condiciones a nivel del mar para el modelo?. El cálculo indica que la presión
requerida sería de 1.900kPa, lo que requiere una potencia muy grande para el
compresor.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 77

Estudio de modelos:
Modelado: A tener en cuenta cuando se comparan modelos y prototipos
Un problema menos evidente es cumplir con la semejanza geométrica en
términos de rugosidad idéntica entre modelo y prototipo. Cuando esto no se
cumple, debe hacerse una corrección por rugosidad que se basa en
extrapolaciones basadas en estudios del efecto de la rugosidad.

Puede ocurrir también que efectos que no sean tenidos en cuenta en el prototipo
aparezcan en el modelo debido a efectos de escala, como por ejemplo en el
caso de modelos a escala de puertos o ríos, en donde si se respeta la escala,
aparecen fuerzas de tensión superficial en el modelo, que no aparecían en el
prototipo.

Si el igualar dos o más parámetros adimensionales no es práctico, la información

se debe obtener de pruebas de modelos combinadas con cálculos teóricos. Se
hace por tanto un trabajo de SIMULACIÓN PARCIAL

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 78

39
 Estudio de modelos:
Simulación parcial
• EXAMPLE: For ship hydrodynamics, Fr similarity is
maintained while Re is allowed to be different.
• Why? Look at complete similarity, i.e.: Rem=Rep y
Frm=Frm

DDG-51 Destroyer

• To match both Re and Fr,

viscosity in the model test is a
1/25th scale model function of scale ratio! This is
not feasible:
If (Lm/Lp)=1/25→νm/νp= 8.0e-3
∴ νp= 1.0e-6×8.0e-3 ~ 8.0e-9 1.0e-7

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 79

Estudio de modelos:
Simulación parcial
Wanapum Dam on Columbia River
• Flows with free surfaces present unique challenges
in achieving complete dynamic similarity.

• For hydraulics applications, depth is very small in

comparison to horizontal dimensions. If geometric
similarity is used, the model depth would be so small
that other issues would arise
Physical Model at
Iowa Institute of Hydraulic Research
™ Surface tension effects (Weber number) would
become important.
™ Data collection becomes difficult.

• Distorted models are therefore employed, which

requires empirical corrections/correlations in order
to extrapolate model data to full scale.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 80

40
 Estudio de modelos:
Simulación parcial
Prototipo canal costero p/mareas de Lp=10m, con μ =1,07×10-3 kg/(ms), ρ
=1.025kg/m3 σ=7,28 ×10-2 N/m, y se escala a Lr=Lp/Lm=1/100. El modelo se hace
funcionar con agua dulce, μ =1,0×10-3 kg/(ms), ρ=998kg/m3, σ=7,28 ×10-2 N/m. Para
una velocidad de la corriente marítima de 0,5m/s cual será la velocidad del modelo?.
Si se cumple Wem=Wep , We = Fuerzas inercia / Fuerzas tensión superficial.

Semejanza entre modelo y prototipo We m = We p

ρ pu 2p Lp ρmum2 Lm ρ pσ m Lp ⎛ m⎞
2
998 m
Wem = Wep ⇒ = ⇒ um = u 2p = ⎜ 0,5 ⎟ × ×100 = 4,93
σp σm ρmσ p Lm ⎝ s ⎠ 1.025 s

• Velocidad calculada para el modelo conservando igualdad de We.

• El valor de We es en este caso: Wem=3,3 ×104 =Wep, ∴ efectos superficiales despreciables.
• La velocidad um 4,93 m/s para el modelo es poco práctica.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 81

Estudio de modelos:
Simulación parcial
Luego para revisar la velocidad del modelo debería analizarse que se cumplan también el
Re y el Fr:
Si se igualasen los Re, para estos datos se tendría una velocidad de 48m/s para el
modelo. Si se calcula el We para cada caso se vería que se tiene: Wep=3,5×104 y
Wem=3,2×106.
Estos números son altos y no habría efectos de tensión superficial, del mismo modo que si
se usa la velocidad definida para igualdad de We, como ya se analizó.

Verificando ahora la igualdad de Fr, Fr=c/(gL), en cada caso tendría Frm=48 (flujo
supercrítico Fr >1) y Frp =5 ×10-2 (flujo subcrítico Fr<1), luego Frm/Frp= 48/(5 ×10-2)=1
×103, ∴ no hay igualdad de Fr si hay igualdad de Re.

El que en un caso el flujo sea subcrítico y en el otro supercrítico introduce efectos de

escala muy significantes que no permitirían una extrapolación de los datos obtenidos en el
modelo hacia el prototipo.

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 82

41
 Estudio de modelos:
Simulación parcial
We = Fuerzas inercia / Fuerzas tensión superficial.

Se verá por tanto que es imposible cumplir a la vez el Re y Fr usando líquidos de

similares propiedades físicas. Los modelos de estudio de oleaje son definidos en
función de que se cumpla el Fr más que el Re.
Los cálculos de velocidad para el modelo basados en similaridad de Fr hace que las
velocidades para el modelo sean del orden de, um= 5 ×10-2 m/s con un Re=5 ×103.
A esta velocidad el Wem=3,4 lo que indica que las fuerzas de inercia aún son más
importantes que las superficiales pero no mucho más grandes, lo que conlleva a
definir una simulación parcial, que implica obtención de resultados con alguna
corrección a posteriori para tener en cuenta estos efectos de tensión superficial.

En casos en donde sólo sea posible la simulación parcial, se tienen en cuenta los
efectos “secundarios” que ya no son menores a partir de posteriores correcciones de
los resultados).
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 83

Estudio de modelos:
Simulación parcial
SIMULACIÓN PARCIAL: Si no pueden igualarse todos los parámetros para casos multi-
parámetricos, se debe seguir una jerarquización de éstos, es decir, deben considerarse los
más relevantes y “relajar” los de menor importancia mediante el demostrar que no es
importante recoger esos efectos a través de cálculos teóricos o la extrapolación de datos,
o bien corregir estos efectos a posteriori de ser esto posible.
Tipo de prueba Parámetros primarios (se Parámetros Parámetros irrelevantes
debe igualar) secundarios (normalmente)
Flujo a baja veloc., sin sup. libre (ej.: flujo en Re, α (ángulo de ataque) ε/l We, Fr, Ma, k o γ
tuberías o accesorios aerodinámicos) (k relac. Calores específicos)
Superficies de lámina, barcos. Fr Re, ε / l We, Ma, k o γ
Aliviaderos de presas, canales, puertos, ríos Fr Re, ε / l We, Ma, k o γ
Flujos a alta velocidad Ma, α , k (γ) Re, ε / l Fr, We
(k relac. Calores específicos)
Maquinaria de fluidos c/u Re, Ma, ε / l We, Fr, k o γ
(relación de 2 velocidades relevantes)

Flujos no estacionarios St Varios Varios

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 84

42
 Estudio de modelos:
Simulación parcial
CONSECUENCIAS:
No siempre pueden obtenerse grupos adimensionales π equivalentes a modelos y prototipo,
incluso cuando se obtiene semejanza geométrica, que también es difícil de cumplir.
Cuando esto ocurre se obtiene SIMILITUD INCOMPLETA o SIMILITUD PARCIAL.
Pese a esto, para el caso de túneles aerodinámicos, muchas veces se pueden extrapolar
resultados de modelos a prototipo si se tiene en cuenta que muchos coeficientes
experimentales dejan de ser función del Re cuando el Re es suficientemente alto como para
garantizar que la capa límite y la estela detrás del objeto sea completamente turbulenta.

Al preparar ensayos de modelos en túneles aerodinámicos debe tenerse en cuenta que el

bloqueo del túnel por la sección frontal del objeto sea como máximo del 7,5% de la
sección transversal del túnel.

Al preparar ensayos de vehículos terrestres si se tiene un piso “móvil”, éste debe ir a la misma
velocidad del viento generado por el ventilador del túnel.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 85

El túnel aerodinámico

Los túneles aerodinámicos son elementos esenciales para la

investigación, ya sea de estudio de modelos para diseño en ingeniería
como para pruebas de modelos en el ámbito de la investigación básica.

Circuito cerrado Circuito abierto

Diseño de túneles

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 86

43
 El túnel aerodinámico
Circuito cerrado

El centro de presiones (cdp) del vehículo es donde

se aplica la resultante de todas las fuerzas
aerodinámicas, es distinto del centro de gravedad
(cdg), que es donde se aplica la resultante de la
fuerzas de gravedad del vehículo. Si estos dos
centros divergen demasiado en su posición el
vehículo será inestable.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 87

El túnel aerodinámico
Circuito cerrado, VKI facilities

The VKI low speed wind tunnel L-1A has a free

jet test section of 3 m diameter and 4.5 m
length. The 4.2 m diameter contrarotating fans
of this facility are driven by a variable speed DC
motor of 580 kW, allowing a continuous
variation of velocity from 2 to 60 m/s. The
contraction ratio is 4 with a typical turbulence
level of 0.3 %.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 88

44
 El túnel aerodinámico
Circuito cerrado, VKI facilities

The VKI Wind Engineering facility L-1B

consists of a specially designed test
section interchangeable with that of the
low speed wind tunnel L-1A and hence
using the same return circuit, motor and
power supply. The rectangular test section
of L-1B is 2 m high, 3 m wide and 20 m long
with a roughened floor to allow the growth
of a turbulent boundary layer similar in
nature to the lower part of the atmospheric
neutral boundary layer.
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 89

El túnel aerodinámico
Circuito abierto, VKI facilities

The L-2A facility is a low-speed, open circuit wind tunnel of the suction type. It
incorporates an air inlet, a circular cross-section settling chamber fitted with
honeycomb and meshes, and an axisymmetric 8:1 contraction faired to an
octagonal test section of 0.28 m width and 1.3 m length. A two-directional
traversing mechanism is located in the test section.

17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 90

45
 El túnel aerodinámico
Circuito abierto, VKI facilities

The L-7 facility is a low speed, open circuit, miniature wind tunnel, with blowing fan, used for
laboratory training in boundary layer and three-dimensional velocity measurements. It is
specially designed to be of economical construction, to require minimal floor space, and to be
easily transportable. The facility is comprised of a centrifugal fan, driven by a 700 W, variable-
speed AC motor, a wide angle diffuser fitted with screens, a settling chamber fitted with a
honeycomb and a 9:1 contraction. Several test sections, with 0.16 m x 0.16 m cross section and
lengths up to 1 m are available. Maximum flow speed is 20 m/s, with a turbulence level of about
0.3 %.
The L-10, L-11, L-12 and L-13 are four identical facilities, of similar design to the L-7; each having
a smaller contraction ratio (5.76:1), a larger test section of 0.2 m x 0.2 m, and a maximum speed of
15 m/s. They are used mainly for laboratory training in several fields, including boundary layers,
17/05/2016 13:32 Mecánica de Fluidos – Miguel G. Coussirat Núñez 91
convective heat transfer, flow in cavities, and the calibration of small instruments.

El túnel aerodinámico
Circuito abierto, VKI facilities

The Mach 3.5 blow-down wind tunnel S-4 has an 8 cm x 10 cm test section. Air is
supplied at stagnation pressures from 3 to 18 bar and a supersonic ejector can be used
to decrease the exhaust pressure downstream of a variable geometry diffuser. A typical
unit Reynolds number for the S-4 is 5 x 107/m.

Depending on the stagnation pressure, the facility can achieve test times from 8 to 25
minutes. Model incidence can be adjusted in the range -10° to +10°, and models can be
injected rapidly into the air-stream when required.

The tunnel is equipped with shadow and Schlieren systems. Instrumentation includes
scani-valves
17/05/2016 13:32
for pressure measurements and three-component strain gauge balances.
Mecánica de Fluidos – Miguel G. Coussirat Núñez 92

46
 El túnel aerodinámico
Circuito abierto, VKI facilities

The VKI Longshot free piston tunnel is a short duration facility that can be operated with either
nitrogen or carbon dioxide, and it is designed to provide the attainment of very high Reynolds
number hypersonic flows.
It has a Mach 14 contoured nozzle of 0.43 m exit diameter and a 6° conical nozzle of 0.60 m exit
diameter which can be used throughout the Mach number range from 14 to 20. Typical Reynolds
numbers at Mach 14 range from 5 x 103 to 15 x 106/m. A high precision incidence mechanism for
pitch, roll, and yaw is mounted in the open-jet 4 m3 test section. Instrumentation includes a
force/moment balance, accelerometers, thin-film and coaxial thermocouples for heat flux
measurements, piezoresistive pressure transducers, and a Schlieren system. A total of 64
channels
17/05/2016of transient recorders with
13:32 a 50
Mecánica kHz sampling
de Fluidos rate are
– Miguel G. Coussirat controlled by a PC.
Núñez 93

47