Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Analisis Dimencional

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 67

ANALISIS DIMENCIONAL

ERICK SANTIAGO ARENALES BAUTISTA

UNIVERSIDAD NACIONAL DE COLOMBIA – SEDE BOGOTÁ

FACULTAD DE INGENIERÍA

DEPARTAMENTO DE INGENIERÍA CIVIL Y AGRÍCOLA

23 DE JUNIO DE 2020
ANÁLISIS DIMENSIONAL

El análisis dimensional es un método para verificar ecuaciones y planificar experimentos sistemáticos. A


partir del análisis dimensional se obtienen una serie de grupos adimensionales, que van a permitir utilizar
los resultados experimentales obtenidos en condiciones limitadas, a situaciones en que se tengan
diferentes dimensiones geométricas, cinemáticas y dinámicas; y muchas veces en casos en que las
propiedades del fluido y del flujo son distintas de las que se tuvieron durante los experimentos.
La importancia del análisis dimensional viene dada por la dificultad del establecimiento de ecuaciones en
determinados flujos, además de la dificultad de su resolución, siendo imposible obtener relaciones
empíricas y teniendo que recurrir al método experimental.
Es importante considerar que si en un experimento en un modelo (a escala geométrica del prototipo), se
pueden obtener las escalas cinemáticas (relaciones de velocidades) y las escalas dinámicas (relaciones de
fuerzas), los resultados adimensionales que se obtienen para el modelo son también válidos para el
prototipo.

Homogeneidad dimensional

En toda ecuación física, cada termino deberá tener las mismas dimensiones: la ecuación debe ser
dimensionalmente homogénea; además la división de todos los términos por uno cualquiera de ellos haría
la ecuación adimensional, y cada cociente sería un grupo adimensional. Las dimensiones de las
magnitudes empleadas normalmente en Mecánica de Fluidos incluyen sólo una o más de las siguientes 4
dimensiones: M (masa), L (longitud), T(tiempo) y ϴ (temperatura):

Tabla 1: Dimensiones de magnitudes empleadas en mecánica de fluidos


Similaridad dinámica

Se define similaridad dinámica entre dos sistemas que sólo se diferencian en un factor de escala de
sus medidas, como la condición que lleva a determinar el comportamiento mecánico de uno de los
sistemas bajo ciertas condiciones, conociendo los valores de las propiedades del otro.
En mecánica de fluidos experimental, un sistema se conoce como prototipo y el otro es el real el cual
va entrar en operación. Es claro que el diseño de un nuevo vehı́culo terrestre, marino o aéreo o un
edificio, pasa primero por la construcción del prototipo con el que se va a predecir el comportamiento
del real en ciertas condiciones, que normalmente son simuladas en un túnel de viento o en una
cámara de agua.

Es una situación común que cientı́ficos e ingenieros que trabajan en mecánica de fluidos
experimental, se enfrenten a problemas donde el número de variables es alto y los modelos teóricos
no proporcionan la información suficiente para predecir cómo se comportarán los sistemas en
determinadas condiciones. La posibilidad de diseñar prototipos con el respaldo de la similaridad
dinámica, es una de las herramientas más poderosas y favoritas cuando toca caminar en la
oscuridad.
Como se mencionó anteriormente, es posible asegurar la similaridad dinámica entre dos sistemas
si éstos comparten la igualdad en valor de ciertos factores adimensionales, por ejemplo el
número de Reynolds, de Froude o el de Arquı́medes, etc. Entre las ganancias de conocer sus valores
está:

• Asegurar similaridad dinámica.

• Simplificar las expresiones matemáticas que relacionan las magnitudes en juego toman- do
casos límite.

• La presentación resumida de los datos experimentales.

• Existen dos casos prácticos donde estos factores adimensionales se pueden obtener:

• Dividiendo los términos de la ecuación o ecuaciones, que gobiernan el comportamiento


dinámico del sistema por un factor dimensional apropiado, de tal forma que los deje
adimensionales. En este caso, es necesario conocer los modelos matemáticos, las con-
diciones de frontera o incluso las ecuaciones constitutivas que describen la situación f´ısica
particular.

• Si no se conocen las ecuaciones, muy común en investigación, estos factores adimensionales


se pueden obtener a través de un análisis dimensional.

Parámetros adimensionales

Las magnitudes que intervienen en el movimiento de un fluido se pueden agrupar en tres tipos:
- magnitudes mecánicas del fluido
- magnitudes térmicas del fluido
- magnitudes del flujo
Tabla 2: Dimensiones y unidades de magnitudes empleadas en mecánica de fluidos

Es frecuente encontrar en los grupos adimensionales factores como el coeficiente de presión,


el coeficiente de arrastre CD, el número de Reynolds Re o el número de Froude Fr. Los parámetros
adimensionales asociados a las magnitudes anteriores vienen determinados por relaciones entre los
diversos efectos que se pueden considerar:

Tabla 3: Parámetros adimensionales más conocidos.


Existen muchos más de estos factores que tienen aplicación práctica en ingeniería. Una referencia más
específica se encuentra en el capítulo 7 del libro Cengel y Cimbala, Mecanica de fluidos.

Sistema Internacional y dimensión fı́sica

A finales de la década de los sesenta del siglo pasado, la comunidad cientı́fica se reunió en
Francia, con el fin de acordar un sistema de unidades de medidas universal que permitiera que
todos hablaran el mismo idioma en cuanto a comunicación de resultados experimentales,
reportes de investigación y avances tecnológicos. De esta forma nació el llamado Sistema
Internacional SI de medidas, en el cual se definen siete magnitudes físicas fundamentales y en
términos de éstas se construyen todas las demás que se llamarán magnitudes fı́sicas
derivadas. El mecanismo que permite construir magnitudes derivadas a partir de las
fundamentales, es la dimensión fı́sica. En la siguiente tabla se muestra las magnitudes
fundamentales seguida de su correspondiente dimensión fı́sica y de sus unidades en el sistema
internacional e inglés:
Magnitud fundamental Dimensión fı́sica2 Unidad SI Unidad Inglesa
Longitud L m ft
Masa M kg lbm
Tiempo T s s
Corriente eléctrica I A A
Cantidad de sustancia m mol mol
Temperatura termodinámica t K R
Intensidad luminosa i cd Cd

Tabla 4: Magnitudes fundamentales – Sistema Internacional.

La dimensión fı́sica de magnitudes fı́sicas muy famosas como la velocidad, la fuerza, la


aceleración, entre muchas otras, quedan en términos de las fundamentales. En la siguiente tabla,
se listan algunas de las magnitudes fı́sicas derivadas más usuales en los cursos de f´ısica de
primeros semestres, en las carreras de ingeniería y ciencias con su correspondiente dimensión fı́sica:
Tabla 5: Magnitudes derivadas del Sistema Internacional.
Magnitud Sı́mbolo Dimensión fı́sica Unidad SI
Velocidad [v] LT −1 m/s.
Aceleración [a] LT −2 m/s2 .
Fuerza [F ] MLT − 2 N (newton).
Energ´ıa [E] ML 2 T − 2 J (joule).
Potencia [p] ML 2 T − 3 W (watt).
Volumen [V ] L3 m3.
Posición [x] L m.
Desplazamiento [∆x] L m.
Distancia [d] L m.
Desplazamiento angular [∆θ] 1 rad (radián).
El hecho de reducir el mundo a siete magnitudes fundamentales tiene consecuencias como
que magnitudes fı́sicas que informan cosas distintas, puedan tener la misma dimensión fı́si-
ca. En la tabla (5) se observa que la posición, el desplazamiento y la distancia, tienen
dimensión de longitud. También ocurre que otras magnitudes fı́sicas, como la posición y el
desplazamiento angular, no tengan dimensión fı́sica. Para indicar la dimensión fı́sica de una
magnitud, se encierra el sı́mbolo que la representa entre paréntesis rectangulares y se iguala
a su correspondiente dimensión en términos de las fundamentales.
Una de las consecuencias más importantes que se dió con este tipo de organización, fue la
formalización de la rama de la fı́sica conocida como Análisis dimensional ; que hoy en dı́a
cuenta con reconocidos grupos de investigación en las universidades más prestigiosas del
mundo.
Uno de los principales teoremas del analisis dimensional es el de homogeneidad dimensional,
que declara que todos los términos de una ecuación de la fı́sica deben tener exactamente la
misma dimensión fı́sica. A continuación un par de ejemplos tradicionales:

Ejemplo: La ecuación que describe la posición de una partı́cula que se mueve en lı́nea
recta bajo la acción de una fuerza constante, es:

donde x es posición, v velocidad, a aceleración, t tiempo y 1/2 es una constante


sin dimensión. Se analiza la dimensión fı́sica de cada uno de los cuatro términos de
la ecuación, ası́:

L1 = L1 + L1T−1 × T 1 + 1 × L1T−2 × T 2
L1 = L1 + L1 + L1 .

Y se verifica que cada uno de los cuatro términos tiene dimensión de longitud.

Ejemplo: La ecuación de Bernoulli expresa la conservación de la energı́a por unidad de


volumen de un fluido que se mueve entre dos puntos dentro de una tuberı́a. Situación
que se describe con la ecuación:

donde P es la presión del fluido en un punto de la tuberı́a, ρ la densidad, v la rapidez,


g la aceleración gravitacional y h la altura. Igual que en el ejemplo anterior, se analiza
la ecuación término a término, ası́:
que, de igual forma, verifica que cada término tiene dimensión fı́sica de M 1 L−1 T −2 .

Matriz dimensional

En el ejemplo anterior de la ecuación de Bernoulli, se observa que la energı́a por unidad de volumen
en un punto es función de la presión, la rapidez, la altura y la densidad del fluido y también del valor
de la gravedad. De tal forma que es posible establecer una relación funcional de la forma:

f (P, v, h, ρ, g) = 0.

Con las dimensiones fundamentales de cada variable física se puede ensamblar un arreglo matricial
dimensional como el siguiente:
P v h ρ g
M 1 0 0 1 0
L -1 1 1 -3 1
T -2 -1 0 0 -2
En la literatura, el arreglo de dimensiones anterior se conoce como matriz dimensional. En la primera
fila se colocan las variables f´ısicas involucradas y en la primera columna las magnitudes f´ısicas
fundamentales y las otras entradas, se llenan con los exponentes de cada magnitud fundamental en
la correspondiente variable f´ısica relevante. El rango j de una matriz, es el tamaño de la mayor
submatriz cuadrada en la que su determinante no sea cero. Por ejemplo, calculando el determinante
formado por las tres primeras filas y columnas de la anterior matriz dimensional, se obtiene:

. .

De lo que se deduce que el rango de esta matriz es j = 3. Si por el contrario, el determinante fuera
cero, entonces se busca otra submatriz de 3 × 3. Si fuera el caso que los determinantes de todas
las submatrices de 3 × 3 fueran cero, entonces, es posible que el rango de esta matriz sea 2 y se
procede a calcular el determinante de todas las submatrices de 2 × 2.
Los problemas que no involucren fuerzas electromagnéticas ni reacciones quı́micas, se pueden
expresar en términos de las cuatro dimensiones fundamentales, masa M , longitud L, tiempo T
y temperatura t. Por otro lado, si los efectos térmicos se pueden despreciar, entonces el
análisis se reduce a las tres magnitudes fundamentales: longitud, tiempo y masa, si no se
pueden despreciar, pero es posible considerar la temperatura en combinación con la cons-
tante de Boltzmann ([k] [T ]) o con una constante de gas ([R] [T ]), es fácil verificar que la
dimensión fı́sica de la combinación en cualquiera de ambos casos es, L2 /T 2 y por lo tanto,
vuelve a sólo requerirse las tres magnitudes fundamentales L, T y M .

Teorema Pi de Buckingham

Es muy conocido en f´ısica que sistemas en campos aparentemente muy diferentes, por ejemplo la
mecánica y el electromagnetismo, tienen comportamientos similares y están descritos por
ecuaciones que tienen exactamente la misma forma y por lo tanto la misma solución.
Por otro lado, los ingenieros obtienen provecho de esta situación de similaridad, que también se
repite en sistemas del mismo tipo o en el mismo campo, pero que tienen diferentes valores de escala
y que da paso a la construcción de prototipos o modelos a escala, que redundan en la posibilidad
real de diseño de experimentos que permitan comprobar hipótesis con una reducción drástica en
los costos de producción y protección de la salud de las personas que van operar tales dispositivos.
Con estos prototipos y la posibilidad de hacer experimentos en situaciones controladas y en
ambientes seguros, se obtiene conocimiento del comportamiento del correspondiente sistema a
escala real en las condiciones que finalmente va a trabajar.
La propiedad que verifica que dos sistemas iguales en forma, pero diferentes en escala, se comporten
de la misma manera, se conoce como similaridad dinámica y la clave para asegurarla es la
igualdad de valor de ciertos factores adimensionales que sólo dependen de la naturaleza de los dos
sistemas.

Números pi

Demostración:

El teorema Pi de Buckingham (1914) establece que si q1, q2, ..., qn son n variables f´ısicas
involucradas en un problema particular y si existe entre éstas una relación funcional de la
forma:

f (q1, q2, . . . , qn) = 0. (1)

Entonces, las n variables siempre se pueden combinar para formar exactamente (n−j) varia-
bles adimensionales independientes, donde j es el rango de la matriz dimensional. Cada
variable adimensional es llamada un número Π o factor adimensional.
El sı́mbolo Π es usado porque las variables adimensionales se pueden escribir como un
producto de las variables q1 , q2 , . . . , qn elevadas a alguna potencia. Ası́, la ecuación (1) se
puede escribir como la relación funcional:

φ (Π1, Π2, . . . , Πn−r) = 0. (2)

La ecuación que permite calcular estos Πi , es:

(3)

Donde VD son las variables fı́sicas relevantes en el problema y VF son las variables dinámicas
que hacen las veces de variables fundamentales (variables repetidas) y los ak se escogen de
manera tal que cada Πi sea adimensional. Finalmente, se puede establecer una relación entre
los Πi, de la forma:

Π1 = f (Π2, Π3, . . . , Πn−j).

La forma arbitraria de escoger las variables repetidas lleva a diferentes conjuntos de factores
adimensionales; sin embargo, en cada caso (n − j) son independientes y en palabras del
álgebra lineal, forman un conjunto completo.

Método de repetición de variables

La ecuación (3) provee un procedimiento formal de calcular los factores adimensionales


(generar los Π’s) de cualquier problema. Este procedimiento es conocido como el método de
repetición de variables y fue publicado incialmente por el cientı́fico ruso Dimitri Riabouchinski en
1911 y un año después se hizo popular por su inclusión en los trabajos del ingeniero Edgar
Buckingham. Buckingham presentó el método como un recetario de seis pasos, que permite al final
obtener los factores adimensionales:

Paso 1

Listar las n variables físicas relevantes del problema.

Paso 2

Escribir cada variable fı́sica en términos de su dimensión fı́sica.


Paso 3

Identificar las j magnitudes f´ısicas fundamentales en las variables físicas del paso 2. De acuer-
do al teorema Pi de Buckingham, el número k de factores adimensionales independientes, es
igual a:

k = n − j.

Paso 4

Escojer j variables fı́sicas que se repetirán para construir cada Π.

Nota 1: Debido a que las variables repetidas tienen el potencial de aparecer en cada Π,
conviene hacer una selección adecuada.

Nota 2: Es posible que la primera selección no lleve a relaciones que aporten a la


solución, entonces, volver al paso 3 y hacer una nueva selección.

Paso 5

Con la ecuación (3), calcular los Pis, uno a la vez, por la agrupación de las j
variables repetidas con una de las que quedaron, forzando a que el producto sea
adimensional. Es una costumbre que el primer Π, designado como Π1 , es el Π
dependiente, el único al lado izquierdo de la lista. Esta convención se tendrá en
cuenta a la conveniencia del ejercicio particular que se esté desarrollando.

Paso 6

Una vez obtenido todos los factores adimensionales, buscar relaciones matemáticas que
permitan conducir un experimento para determinar el valor de las constantes
adimensionales para posteriormente intentar predecir.

Ejemplo ilustrativo

Fuerza de ascenso (sustentación) sobre el ala de una aeronave

El siguiente ejemplo tiene como objetivo hacer énfasis en el poder de predicción que se logra
cuando se toman casos límites de los Pis.
Variables f´ısicas involucradas: densidad del aire ρ, viscosidad del aire µ, velocidad del aire v,
ángulo de ataque α del ala respecto al flujo del aire, velocidad del sonido c, longitud del perfil
del ala L y la fuerza de ascenso F , (ver figura 1). En total n = 7 variables físicas relevantes en
el problema. Las magnitudes fundamentales son longitud, masa y tiempo, por lo tanto, j = 3.
Se deben calcular 4 Pis (n − j = 4).
Al aplicar el teorema Π, se busca la forma funcional, Π1 = f (Π2, Π2, Π4). Procediendo, se
obtienen:

Figura 1.: Ala de aeronave dentro de un fluido de densidad ρ y viscosidad µ.

Los parámetro adimencionales que se dan como resultado son:

Conocidos como coeficiente de Lift, número de Reynolds, número de Mach y el ángulo de


ataque que en este caso es una variable fundamental en el problema. Por lo tanto,
CL= f (α, Re, Ma) .
Conocer esta forma funcional limita enormemente los experimentos. Basta con estudiar solo
algunos regı́menes para saber la dependencia de estos parámetros. De lo contrario el estudio
se debe empezar con las siete variables f´ısicas y probablemente en un tiempo razonable de
laboratorio nunca se hubiera encontrado las leyes correctas.
Para ejemplificar, si se trata de un fluido subsónico el número de Mach no es importante
y la ley de la sustentación sólo dependerı́a del ángulo de ataque y el número de Reynolds.
Pero en el limite de muy bajas velocidades donde el flujo de aire se puede considerar con
densidad constante y no turbulento, la dependencia de Re desaparece y la ley se reduce a la
forma simple:
CL = f (α),
es decir, la sustentación sólo depende del ángulo de ataque de la superficie respecto al flujo
del aire. Sólo los experimentos permiten conocer esta forma funcional y su lı́mite de validez
y cuándo entran en juego las otras variables en el problema como Re y Ma.
El ejemplo anterior es netamente ilustrativo, acontinuacion se realizarán varios ejercicios
donde se profundizará en cada detalle para poder entender el procedimiento general de la
aplicación del teorica Pi de Buckinham

Ejercicios Teorema Pi de Buckingham

1. Esfera que cae sin tener en cuenta la fricción del aire

En la mayorı́a de textos de mecánica de fluidos, se ejemplifica el método con el conocido


problema de un cuerpo que se arroja en un medio viscoso. En este momento, con el fin de
profundizar en la importancia que tiene el procedimiento en mecánica de fluidos y para
ilustrar la aparición de factores adimensionales famosos, se abordará este mismo
problema. Considere una esfera que cae en el vacı́o. La elevación y en un instante es función
del tiempo t, la rapidez inicial v0, la altura inicial y0 y el valor de la constante de gravedad g.
Se supone que son conocidas las dimensiones fundamentales de estas magnitudes f´ısicas.
Aplicando el método de repetición de variables:

Paso 1

De acuerdo al enunciado existen cinco variables físicas relevantes en este problema, por lo
tanto n = 5, así:

y = f (t, v0, y0, g) , n = 5.

Paso 2

Las variables del paso 1 se escriben en términos de sus dimensiones fundamentales. Se reco-
mienda escribir cada dimensión con su respectivo exponente, ya que, esto ayudará con los
procedimientos algebraicos posteriores.

y t v0 y0 g
{L } {T } {L T }
1 1 1 −1
{L } {LT −2}
1

Paso 3

Se identifica la presencia de sólo dos magnitudes fundamentales: el tiempo T y la longitud


L, por lo tanto j = 2.

Reducción j = 2.

El número de Πs que predice el teorema Pi de Buckingham, es:

Número esperado de ΠJ s, k = n − j = 5 − 2 = 3.
Paso 4

Como se identificó que j = 2, se deben escoger dos variables fı́sicas a repetir. Aunque, en la
literatura se establecen ayudas o guı́as para la selección de estas variables, la intensión en este punto
es acercarse lo más ingenuamente posible a la solución con el método. Ası́ que, sin mayor
refelxión se proponen como repetidas las siguientes variables:

Variables repetidas: v0 y y0 .

Con el uso continuado del método, se podrá observar que hay algo de arte en reconocer cuál es la
mejor selección.

Paso 5

Con la ecuación (3) se combinan las variables repetidas en productos con cada una de las otras
variables. Se calcula cada Pi, uno a la vez. El primero, Π1, se llama el Pi dependiente y en este caso
se forma con la variable dependiente y.

donde a y b son constantes que deben ser determinadas y que se ajustan para que la dimensión del Pi
sea cero. Cada factor que compone al Pi se expresa en términos de sus dimensiones
fundamentales.

Dimensiones de Π

Como las dimensiones fundamentales son por definición independientes, se iguala la suma de los exponentes
de cada una de ellas a cero, así:

Tiempo: T 0 = T −a , 0 = −a, a = 0.

Longitud: L0 = L 1 L aL b, 0 = 1 + a + b, b = −1 − a, b = −1.
Por lo tanto, el primer Pi, queda:
y
Π1 = . (4)
y0
De manera similar se crea el primer Pi independiente, es decir, Π2. Combinando los factores
que se repiten con la magnitud tiempo t, luego:
Tiempo: T 0 = T 1−a, 0 = 1 − a, a = 1.

Longitud: L0 = La+b, 0 = a + b, b = -1.


Por lo tanto:

(5)

Finalmente, el segundo Pi independiente (Π3), se obtiene combinando los factores repetidos


con la variable g. Procediendo de forma idéntica, se llega a:
gy0
Π3 = (6)
v02

Paso 6

Con los tres Pis encontrados, ecuaciones (4), (5) y (6), se pueden buscar relaciones que
permitan establecer una ley que describa la caída de la esfera.
Antes de terminar con este ejemplo, se puede observar que de forma natural se ensamblaron
parámetros adimensionales correspondientes a las magnitudes de posición y tiempo, y y t
respectivamente, así:
y v0t
z= , y τ=
y0 y0

z corresponderı́a a la versión adimensional de la posición y τ al correspondiente adimensional del


tiempo. En el caso de Π3, que aparentemente no es la forma adimensional de ninguna magnitud, se
puede reescribir como uno de los factores adimensionales más famosos usados en mecánica de fluidos,
el conocido número de Froude, ası́:

Pi 3 modificado:

En la mecánica de fluidos y termodinámica de sistemas abiertos este numero puede ser


determinante. Para explorar un poco lo mencionado, se observa que la posición adimensional es
función de y y del tiempo t,

z = f (y, t),

entonces, la primera derivada respecto al tiempo normalizado es:


por lo tanto,
∗ 1
z= ẏ (7)
v0
La segunda derivada, relacionada con la aceleración, también se obtiene con la regla de la
∗∗
cadena. Antes de proceder, se hace énfasis en que según la ecuación (7), z es función de ẏ
y de el tiempo, as´ı:

entonces,

(8)
La segunda ley de Newton aplicada al problema de la caı́da libre lleva a la ecuación,

ÿ = −g,

que con 2-7 y 2-8 se reescribe como,

v20 ∗∗
z = −g.
y0

Relacionando con el número de Froude se llega a la versión adimensional de la caı́da libre:

(9)

Al integrar dos veces la ecuación (9), se obtiene:


1 2
z(τ ) = 1 + τ − τ . (10)
2Fr2
El número de Froude combina los tres parámetros del problema original, y0 , v0 y g. La ecua-
ción (10) es una familia de curvas que resumen muchas combinaciones de dichos parámetros.
Figura 2.: Gráfica de z vs τ para distintos valores del número de Froude al cuadrado.

La reducción paramétrica de tres variables a una, que se logró con el método de repetición y
los Pis encontrados, hace mucho más fácil estudiar el sistema; pues el esfuerzo experimental
sólo se debe concentrar en las curvas z contra τ como las mostradas en la figura 2. En
caso contrario, se debe buscar la relación de la altura contra el tiempo dejando el valor de
g constante, luego altura contra gravedad dejando el tiempo constante y finalmente tiem-
po contra gravedad dejando la altura constante. Ahora, dejar la gravedad constante no es
problema pero cuando toque variarla la situación se complica mucho.
El poder predictivo que se logra con el procedimiento anterior se puede ejemplificar en la
siguiente situación; considere que en la superficie de marte se lanza un objeto con rapidez
de 2 m/s desde una altura de 2 m. La gravedad marciana es de 3,711 m/s2 . ¿Cuánto tiempo
tarde en llegar al piso?
El cuadrado de Froude es,

En la gráfica 2 se interpola que el valor de τ correspondiente a un Fr2 de 0,54 es de


aproximadamente 1,7, así:
Luego el tiempo que demora la caída es de aproximadamente 1,7 s, que al contrastar con el
valor que se obtiene de las ecuaciones del movimiento uniformemente acelerado, es práctica-
mente el mismo.

En general, el análisis dimensional arroja de manera natural y continuada una serie de estos
factores que por su importancia en el diseño de dispositivos, han tomado el nombre de per-
sonajes destacados de las ciencias y la ingenier´ıa, tales como Froude, Reynolds, Arquímedes,
entre otros.
Continuando con la discusión anterior, el Pi tres original no está de ninguna manera errado;
de hecho tiene cierto parecido con la famosa relación v 2f = v20 + 2a∆y de la caı́da libre,
pero si no representa ventaja matemática alguna, es mejor dejarlo en términos de un factor
adimensional famoso como el número de Froude. En el texto de Cengel se mencionan
algunas reglas para manipular los Pis de manera tal, que lleven a este tipo de factores. Pero
lo anterior no entra en los propósito de este trabajo y por lo tanto no se hará mayor énfasis
en esto.
Como se ha insistido, todos los términos de una ecuación deben tener la misma dimensión
f´ısica, teorema de homogeneidad dimensional. Pero dimensiones o unidades de medida son
construcciones humanas muy convenientes si se desea crecer en tecnología y lograr relaciones
de intercambio de mercancı́as, ciencia y otros. Debido a que ningún sistema de unidades es
mejor que otro, se puede plantear como una necesidad, que cualquier problema f´ısico debe
ser expresable en forma completamente adimensional, es decir, tener una formulación sin
dimensiones. Y los parámetros usados para construir variables adimensionales dependientes
o independientes, tendrán que aparecer en las ecuaciones y en las condiciones de frontera.
Si el procedimiento se hace de manera apropiada, habrá una reducción en la dependencia
paramétrica de la formulación, generalmente, por el número de unidades independientes.
Esta reducción paramétrica se llama similitud y las similitudes facilitan en gran medida la
correlación de los datos experimentales, permitiendo obtener una solución al problema y de
paso, convirtiéndose en la base de la similaridad dinámica anteriormente mencionada.

2. Burbuja de jabón
Intuitivamente se comprende que la estabilidad de una burbuja de jabón está
directamente relacionada con la diferencia de presiones entre el interior y exterior de
ésta. Ir un poco más allá de esta conjetura intuitiva no es fácil, si no se conoce algo de los
temas de tensión superficial y presión. La intensión de este ejercicio, es encontrar una
expresión matemática simple que permita una explicación que integre otros elementos y que
pueda en principio dar origen a experimentos que permitan validarla.
La figura 3, muestra las variables que a primera vista tienen relevancia en el problema:
∆P : La diferencia de presiones entre el interior y el exterior de la burbuja es la causa
que la estructura se infle, dando paso a las siempre curiosas esferas hechas de jabón.
Figura 3: Burbuja de jabón. A primera vista las magnitudes fı́sicas que intervienen tienen que
ver con la geometrı́a, las fuerzas de tensión superficial y la diferencia de
presiones.

γ: Se entiende que la superficie de la pelı́cula de jabón debe tener cierta resistencia


a rasgarse y que excedida ésta la burbuja explota. Tal resistencia es conocida
como tensión superficial y sin mayor rigor y de manera un poco ingenua se tomará
como una variable fı́sica relevante de la cual no se dará más explicación que sus
unidades SI y por lo tanto su dimensión fı́sica.

R: El radio de la esfera está directamente relacionado con el espesor de la pelı́cula que a


su vez tiene relación con la tensión superficial. Se aprecia fácilmente, que si la burbuja
se infla continuamente el radio aumenta hasta explotar.

De acuerdo a lo anterior, el número de variables fı́sicas relevantes en el problema es n = 3.


La dimensión fı́sica de cada una de las tres variables relevantes en el problema, son:

Magnitud R ∆P γ
Unidad SI m Pa N/m
Dimensión fı́sica {L } {M L T }
1 1 −1 −2
{M 1 T −2 }

Por simple inspección se verifica que son 3 las mágnitudes fı́sicas fundamentales presentes
en las variables relevantes: el tiempo, la longitud y la masa, j = 3. El número de Pis, n − j,
ser´ıa cero. Sin embargo, el tiempo se puede descartar ya que la presencia de una burbuja de
jabón se puede considerar una situación estática, y por lo tanto j = 2. Ası́, el número de Pis
a calcular es:

n − j = 3 − 2 = 1.

Se propone tomar como variables repetidas el radio R y la tensión superficial γ. Sólo hay
que calcular un Pi, entonces:

Π = ∆P Ra γ b = M 1 L−1 T −2 × La × M 1 T −2 = M 1+b L−1+a T −2−2b . (1)


Asegurar adimensionalidad lleva al siguiente sistema lineal:

1 + b =0
−1 + a = 0
−2 − 2b = 0

cuya solución es:

b = −1, a = 1.

Por lo tanto la forma del Π es:


∆P · R
Π= .
γ

En la eventualidad de un sólo Π, inmediatamente se busca relacionar las variables presentes


intentando una ecuación, que salvo una constante adimensional por determinar, permita
obtener la expresión matemática deseada. En este caso, al despejar la diferencia de presiones:

ecuación que tiene la forma de la conocida Ley de Laplace.

3. Fluido fluyendo por una tuberı́a

Para encontrar una expresión matemática que dé una primera aproximación a la descripción
del comportamiento de un fluido (como el agua) que se mueve dentro de una tuberı́a de
diámetro d, longitud l, conociendo que en la superficie interna se formó una rugosidad de
espesor e.

Figura 4: Fluido que se mueve dentro de una tuberı́a de diámetro d y longitud l.


Con el teorema Pi se pretende encontrar una expresión matemática que permita describir
algo de un fluido que se mueve a través de una tuberı́a, pero sin entrar en demasiados detalles
técnicos ni hacer uso de los principios de la fı́sica. Para resolverlo se implementa el método
de repetición de variables.
A partir de la información dada en el enunciado, las variables fı́sicas relevantes son el diámetro
d, la longitud l y el espesor de la rugosidad e de la superficie interior de la tuberı́a, la densidad
ρ , la viscosidad µ y la rapidez v del fluido. También se debe considerar el cambio de presión
∆P que experimenta el fluido a medida que avanza. Lo anterior queda bien expresado por la
relación funcional:

f (∆P, d, l, e, v, ρ, µ) = 0.

El conteo arroja siete variables físicas que son relevantes en el problema, luego:

Número de variables relevantes: n = 7.

Se reescriben en términos de las magnitudes fundamentales:

∆P d l e v ρ µ
Pa m m m m/s kg/m3 Pa · s
{M L T }
1 −1 −2
{L1} {L1} {L1} {L 1T −1 } {M 1L−3} {M 1 L −1 T −1 }

Por simple inspección se ve que sólo tres magnitudes fundamentales están presentes en el
problema, luego:

Magnitudes fundamentales: j = 3.

La matriz dimensional tiene rango r = 3. Como son n = 7 variables en el problema, el


número de factores adimensionales deben ser:

Número de Pis a calcular: (k = n − j) k = 4.

De las 7 variables se seleccionan 3 cualquiera, pero se sugiere que estas variables repetidas
tengan diferente dimensión fı́sica, es decir, no se seleccionará d, l y e. Además es importante
que entre las elegidas estén presentes las dimensiones fundamentales M , L y T .
La selección que se propone tiene en cuenta dos parámetros de la tuberı́a y uno del fluido
y se escogieron de forma completamente arbitraria pero siguiendo los consejos del párrafo
anterior. Entonces:

Variables repetidas: v, d y ρ.
Se recalca que la selección fue completamente arbitraria y que cualquier otra es igualmente
viable. Cada factor adimensional se construye combinando las tres variables que se repiten
con una de las que no se seleccionaron. Por ejemplo, tomando el primer factor adimensional,
como:

Π1 = ∆P vadbρc.

Los exponentes a, b y c se obtienen del requerimiento de dejar a Π1 sin dimensiones. Por lo


tanto,

Igualando ´ındices, se deduce que a = −2, b = 0, c = −1, de tal forma que:

Con procedimientos similares, se calculan los otros pis,

Por lo tanto, la relación buscada, es:


A primera vista se puede reconocer que de nuevo aparece el número de Reynolds en uno de
los Pis y también que, desde la ecuación de Bernoulli, ρv 2 tiene las mismas unidades de ∆P .

4. Fı́sica Nuclear - Parámetro de impacto

A finales del siglo XIX, Thompson, debido al descubrimiento que él mismo hizo del electrón,
desarrolló el modelo del pudı́n de uvas pasas como una primera aproximación para describir
cómo estaba conformada la materia. Continuando con este trabajo, ya en el siglo XX, Ruther-
ford puso a prueba la propuesta de Thompson en un épico experimento.
En 1898 Marie y Pierre Curie descubrieron un elemento, que llamaron Radio, que emit´ıa ra-
diación en tres formas: rayos beta (electrones), rayos alfa (núcleos de Helio) y rayos gamma,
y era notablemente más activo que los ya conocidos uranio y torio descubiertos por Becque-
rel y Marie Curie respectivamente. En 1907, Rutherford utilizaba el radio como fuente de
part´ıculas alfa para sus famosos experimentos. Hac´ıa pasar a las part´ıculas alfa que emit´ıa
la fuente por una rejilla (colimador) para obtener un haz muy estrecho, que impactaba di-
rectamente sobre una delgada lámina de oro que lo dispersaba haciéndolo incidir, una vez
emerg´ıa, sobre una superficie de sulfuro de zinc, ver figura 5.

Figura 5.: Experimento de Geiger y Marsden.

Despúes de los sorprendentes resultados, Rutherford tuvo que cambiar el planteamiento ori-
ginal y para explicar lo que estaba pasando, propuso un modelo planetario como la estruc-
tura del átomo. Explicaba que la fuerte desviación que experimentaban algunos proyectiles
(partı́culas alfa) se debı́a a lo cerca que pasaban de los pesados núcleos del material de la hoja
de oro que utilizó como blanco. Luego, para Rutherfor era claro que el ángulo de dispersión
estaba directamente relacionado con la distancia más cercana del proyectil, partı́cula alfa,
al núcleo dispersor, distancia conocida como parámetro de impacto b, que se aprecia en la
figura 6. Entonces, se propuso encontrar la distribución del ángulo de dispersión .
Figura 6: Esquema de la dispersión de Rutherford. Partı́culas alfa se disparan contra un
núcleo pesado y debido a la interacción entre ellos se dispersan un ángulo Φ.

Después de muchos experimentos y aplicación rigurosa de las leyes de la fı́sica y los principios
de conservación del momento y la energı́a, Rutherford obtuvo la siguiente expresión para cada
valor del parámetro de impacto:

. (2)
El objetivo es encontrar una ecuación para el parámetro de impacto usando el teorema Pi de
Buckingham y con muy poco respaldo de teorías f´ısicas.
Tomando como magnitudes relevantes en el problema la velocidad del proyectil v, el paráme-
tro de impacto b, el ángulo de dispersión φ, la fuerza eléctrica entre las dos partı́culas Fe y la
masa del proyectil mα, n = 5.

f (v, b, φ, Fe, mα) = 0.

El conteo arroja cinco variables físicas que son relevantes en el problema, luego:

Número de variables relevantes: n = 5.

Las cinco variables se reescriben en términos de las magnitudes fundamentales, ası́:

v b φ Fe mα
m/s m rad N kg
{M L } {L } {1} {M L T }
1 −1 1 1 1 −2
{M 1}

Las dimensiones fı́sicas fundamentales que están presentes son el tiempo (T), la masa (M) y
la longitud (L), j = 3. En este caso tenemos n − j = 2. Luego, de acuerdo al teorema Pi
de Buckingham, existen dos Π’s:

Número de Pis a calcular: (k = n − j) k = 2.

Para construirlos se deben elegir tres variables dimensionales como variables repetidas que
van a sustituir a las dimensiones f´ısicas fundamentales presentes en el problema. Se propone:

Variables repetidas: v, b, Fe.

Se puede considerar el ángulo de dispersión como un Π, ya que es una magnitud sin dimensión
f´ısica, as´ı:

Π1 = φ. (3)

Cálculo del segundo Π:

Π2 = mα v a bb Fe c (4)
en dimensiones:

Asegurando adimensionalidad del Π2,

c = −1
a=2 (5)
b = −1

reemplazando la ecuación 5 en la 4,

(6)

La fuerza eléctrica a un valor tı́pico está dada por:

ke(2e)(Ze) 2 Z ke e2
Fe = = , (7)
b2 b2
reemplazando en 6, tenemos:
mαv2b
Π2 = . (8)
2Zkee2
El teorema de Buckingham establece que existe una relación entre los dos Π’s, es decir:

Π2 = f (Π1).

Teniendo en cuenta lo anterior y despejando b de 3-8, finalmente se obtiene:


2Zkee2
b= f (φ). (9)
mα v2
Expresión que aparte de estar de acuerdo con la de Rutherford, cumple el objetivo principal.
Nota: A pesar que no se incluyó la carga eléctrica como una magnitud fundamental, final-
mente se tuvo en cuenta al tomar la fuerza en el caso particular de interacciones eléctricas, es
decir, la ley de Coulomb.
5. Caida de una pluma y una esfera

Figura 7: Una pluma y una esfera se sueltan al mismo tiempo y desde la misma altura,
¿cuál llega al piso primero?.

Situación de clase: Se dejan caer una pluma y una esfera (metálica o de vidrio) desde la
misma altura y al mismo tiempo, figura 7, ¿cuál llega primero al piso?
La experiencia demuestra que la esfera cae primero. Sin mayor reflexión, se puede argumentar
que la forma y el peso de los objetos y la interacción con el aire tienen que ver con el resultado
final. Pero pasar de la argumentación a una explicación más técnica parece en principio lejana
y seguramente llena de principios fı́sicos que sólo una persona con formación en ciencias o
ingenierı́a lo podrı́a lograr.

Solución con el teorema Π


¿Qué parámetros fı́sicos influyen en la caı́da?

Geométricos: En ambos casos, la pluma y la esfera, se debe tener en cuenta


una longitud tı́pica el radio para la esfera y la longitud para la pluma. Para
lograr una descripción general se dejará como variable geométrica R. Al final,
para cada caso par- ticular se interpretará como un radio o como una longitud,
según sea el caso. Variable f´ısica: R.

Aspectos del fluido: Tanto la densidad (ρ) como la viscosidad (µ) del aire, juegan un
papel importante. Variables f´ısicas: ρ y µ.
Factores dinámicos: En este problema la velocidad y la fuerza juegan un papel
importante en la cinemática de cada objeto. Variables fı́sicas: v, F .

De acuerdo al análisis anterior, el número de variables fı́sicas que intervienen en el movimiento


de la pluma y la esfera son cinco: R, ρ, µ, v y F .

Número de variables relevantes: n = 5.


La dimensión fı́sica de cada una de las cinco variables relevantes en el problema, son:
R ρ µ v F
m kg/m3 kg/ms m/s kg m/s2
{L1} {M 1L−3} {M 1 L −1 T −1 } {L 1T −1 } {M 1 L 1 T −2 }
Observando las dimensiones f´ısicas se concluye que las magnitudes f´ısicas fundamentales,
son la masa M , el tiempo T y la longitud L. Luego j = 3, que es el orden de la matriz
dimensional. De acuerdo al teorema Pi de Buckingham, con este análisis se deben obtener dos
Πs:

Número de Pis esperados: k = n − j = 2.

Como primer intento se escogen como variables repetidas: el radio, la velocidad y la densidad:

Variables repetidas: R, v, ρ.

Construcción de los Pis.

Para Π1 se selecciona la viscosidad y por lo tanto:

Π1 = µRavbρc.

En términos de las dimensiones fundamentales:


. a . b . c
Π
1 = M L T
1 −1 −1
× L1 × L1 T −1 × M 1 L−3 = M 1+c L−1+a+b−3c T −1−b .

Forzando el Pi para que sea adimensional, se llega al siguiente sistema:

−1 + a + b − 3c = 0,
−1 − b = 0,
1 + c = 0.

Con solución:

a = −1, b = −1, c = −1.

Al sustituir,
µ 1
Π1 = = . (10)
Rvρ Re
Factor adimensional relacionado directamente con el número de Reynolds.
Se procede de igual manera para calcular el segundo Π. Para Π2 se selecciona el paráme-
tro faltante F :

Π2 = FRavbρc.

En términos de las dimensiones fundamentales:


a . b . c
Π2 = M 1 L1 T −2 × L1 × L1 T −1 × M 1 L−3 = M 1+c L1+a+b−3c T −2−2b .

Forzar la adimensionalización del Pi lleva al sistema lineal:


1 + a + b − 3c = 0,
−2 − 2b = 0,
1 + c = 0.

Con solución:

a = −2, b = −2, c = −1.

Al sustituir,
F
Π2 = . (3-11)
R2 v2 ρ

De acuerdo al teorema Π existe una relación entre los Pis del ejercicio, pero buscar la ade-
cuada requiere de paciencia y habilidad que se va ganando con el uso continuado del método.
En este caso, el primero está relacionado con el número de Reynolds. El segundo, permite
visualizar una ley de fuerzas. Al relacionar los dos Pis, se obtiene:

(12)
Si el número de Reynolds no es importante en el problema, la fuerza se ve afectada por el
cuadrado de la velocidad, haciendo que ésta crezca muy rápido de tal forma que el sistema (la
pluma) se sustenta en el aire, así:

F ∼ ρR2v2,
que explicarı́a por qué la pluma queda como suspendida en el aire cayendo lentamente.
Si por el contrario, el Reynolds toma importancia, la ecuación 12 se multiplica y divide por µ y
se reorganizan los factores, as´ı:

Reemplazando Re = ρRv/µ, se obtiene:

F ∼ µReRv.

que tiene la forma de la conocida ley de Stokes, en la que la fuerza de fricción sobre un cuerpo
que se mueve dentro de un fluido viscoso es directamente proporcional a la velocidad, representando
bien lo que le sucede a la esfera.

6. Caida libre

Paso 1
Variables f´ısicas relevantes: altura y, masa m, gravedad g, tiempo t. De acuerdo a lo anterior,
el número de variables fı́sicas relevantes en el problema es n = 4.

Paso 2
La dimensión fı́sica de cada una de las cuatro variables relevantes en el problema, son:

Magnitud y m v t
Unidad SI m kg m/s s
Dimensión fı́sica {L1} {M 1} {L1T−1} {T1}

Paso 3
Por simple inspección se verifica que son 3 las mágnitudes fı́sicas fundamentales presentes
en las variables relevantes del paso 2: tiempo, longitud y masa, j = 3. Ası́, el número de Pis
a calcular son:

n − j = 4 − 3 = 1.

Paso 4
Se propone como variables repetidas: la masa m, la gravedad g y el tiempo t. Se deja como
no repetida la altura y.

Cálculo del Pi
Sólo hay que calcular un P i, entonces:
(1)

Que para asegurar adimensionalidad lleva al sistema:

a=0
1 + b =0
−2 + c = 0

cuya solución es:

a = 0, b = −1, c = −2.

Por lo tanto la forma del Π es:

Paso 6
En la eventualidad de un sólo Π, se busca relacionar las variables presentes intentando una
ecuación, que salvo una constante adimensional por determinar, permita obtener la expresión
matemática deseada. En este caso, al despejar la altura:

y = Πgt2.
7. Pendulo simple

Paso 1
Variables f´ısicas relevantes: longitud de la cuerda L, masa m, gravedad g, el periodo T . De
acuerdo a lo anterior, el número de variables fı́sicas relevantes en el problema es n = 4.

Paso 2
La dimensión fı́sica de cada una de las cuatro variables relevantes en el problema, son:

Magnitud L m g T
Unidad SI m kg m/s2 s
Dimensión fı́sica {L1} {M 1} {L1T−2} {T1}

Paso 3
Por simple inspección se verifica que son 3 las mágnitudes fı́sicas fundamentales presentes
en las variables relevantes del paso 2: tiempo, longitud y masa, j = 3. Ası́, el número de Pis
a calcular son:

n − j = 4 − 3 = 1.

Paso 4
Se propone como variables repetidas: la masa m, la gravedad g y la longitud L. Se deja como
no repetida el periodo T .

Cálculo del Pi
Sólo hay que calcular un P i, entonces:

. (1)
Que para asegurar adimensionalidad lleva al sistema:

1 − 2c = 0
b=0
a + c =0

cuya solución es:

a = −1/2, b = 0, c = 1/2.

Por lo tanto la forma del Π es:

Paso 5
En la eventualidad de un sólo Π, se busca relacionar las variables presentes intentando una
ecuación, que salvo una constante adimensional por determinar, permita obtener la expresión
matemática deseada. En este caso, al despejar el periodo:

8. Aceleración Centrípeta

Paso 1
Variables fı́sicas relevantes: masa m, rapidez v, radio R, aceleración a. De acuerdo a lo
anterior, el número de variables fı́sicas relevantes en el problema es n = 4.

Paso 2
La dimensión fı́sica de cada una de las cuatro variables relevantes en el problema, son:

Magnitud m v R a
Unidad SI kg m/s m m/s2
Dimensión fı́sica {M 1} {L1T−1} {L1} {L1T−2}

Paso 3
Por simple inspección se verifica que son 3 las mágnitudes fı́sicas fundamentales presentes
en las variables relevantes del paso 2: tiempo, longitud y masa, j = 3. Ası́, el número de Pis
a calcular son:
n − j = 4 − 3 = 1.

Paso 4
Se propone como variables repetidas: la masa m, la rapidez v y el radio R. Se deja como no
repetida la aceleración a.

Cálculo del Pi
Sólo hay que calcular un P i, entonces:

(1)

Que para asegurar adimensionalidad lleva al sistema:

a=0
1 + b + c =0
−2 − 2b = 0

cuya solución es:

a = 0, b = −2, c = 1.

Por lo tanto la forma del Π es:

Π = av−2R1.

Paso 5
En la eventualidad de un sólo Π, se busca relacionar las variables presentes intentando una
ecuación, que salvo una constante adimensional por determinar, permita obtener la expresión
matemática deseada. En este caso, al despejar la aceleración:
9. Trabajo y energía cinetica
10.
Paso 1
Variables f´ısicas relevantes: Fuerza F , masa m, rapidez v, distancia d. De acuerdo a lo
anterior, el número de variables fı́sicas relevantes en el problema es n = 4.

Paso 2
La dimensión fı́sica de cada una de las cuatro variables relevantes en el problema, son:

Magnitud F m v d
Unidad SI N kg m/s m
Dimensión fı́sica {M 1 L1 T −2 } {M 1 } {L1 T −1 } {L1 }

Paso 3
Por simple inspección se verifica que son 3 las mágnitudes fı́sicas fundamentales presentes
en las variables relevantes del paso 2: tiempo, longitud y masa, j = 3. Ası́, el número de Pis
a calcular son:

n − j = 4 − 3 = 1.

Paso 4
Se propone como variables repetidas: la masa m, la rapidez v y la distancia d. Se deja como
no repetida la fuerza F .

Cálculo del Pi
Sólo hay que calcular un P i, entonces:
Que para asegurar adimensionalidad lleva al sistema:

1 + a =0
1 + b + c =0
−2 − 2b = 0

cuya solución es:

a = −1, b = −2, c = 1.

Por lo tanto la forma del Π es:

Π = Fm−1v−2d1.

Paso 5
En la eventualidad de un sólo Π, se busca relacionar las variables presentes intentando una
ecuación, que salvo una constante adimensional por determinar, permita obtener la expresión
matemática deseada. En este caso, al despejar el producto Fuerza por distancia (Trabajo):

Fd = Πmv2.
TEORÍA DE MODELOS

En los ensayos experimentales del flujo en un determinado prototipo, a veces no es


posible realizar los ensayos con el propio prototipo, por su tamaño o por la dificultad de
reproducir las condiciones reales de flujo, con lo que se realizan los ensayos con modelos a
escala (geométricamente semejantes),.

Por ejemplo, en el estudio experimental de hélices marinas, se realizan dos tipos


fundamentales de ensayos con hélices modelo a escala reducida: los de autopropulsión en un
canal de agua dulce y los de cavitación en un túnel de cavitación cerrado con agua caliente y a
depresión; en la evaluación del comportamiento de una carena, se realizan ensayos de arrastre
con modelos a escala reducida en agua dulce, tanto en canales de aguas tranquilas como en
piscinas con generadores de olas.

La teoría de modelos permite obtener las condiciones de ensayo del modelo a partir de
las condiciones de flujo del prototipo y las magnitudes del prototipo a partir de las medidas
experimentales del modelo.

SEMEJANZA

Prototipo, modelo y sus respectivos flujos considerados, están relacionados entre si


por tres tipos de semejanza: geométrica, cinemática y dinámica.

Semejanza geométrica, con un factor de escala de longitudes constante entre modelo y prototipo
(NL):

Semejanza cinemática del campo de velocidades, con un factor de escala de velocidades


entre modelo y prototipo:

la relación entre los dos factores de escala: de longitudes y de velocidades, viene determinada
por el factor de escala de tiempos:
Semejanza dinámica de los campos de las distintas fuerzas que pueden intervenir en el flujo,
con un factor de escala de fuerzas, que debe ser constante, entre modelo y prototipo:

El factor de escala de fuerzas, es el que va a permitir establecer las condiciones del flujo en el
ensayo del modelo a partir de las condiciones del flujo en el prototipo, y obtener ““fuerzas,
potencias y rendimientos”” del prototipo a partir de sus correspondientes valores
experimentales en el modelo.

Los campos de fuerzas que pueden aparecer en la interacción de un fluido y un objeto, pueden ser:

Fuerzas de inercia, determinadas por la variación temporal de la cantidad de movimiento, y


cualitativamente podemos expresarlas por:

Fuerzas de rozamiento por viscosidad, determinadas por el campo de tensiones, que a su vez
viene determinado por la viscosidad y el campo de velocidades, y podemos expresarlas cualitativamente
por:

Fuerzas gravitatorias, determinadas por la posición en el campo gravitatorio, expresadas por:

Fuerzas de presión, determinadas por el campo de presiones:

Fuerzas de elasticidad, determinadas por la compresibilidad del fluido, o bien por la velocidad de
pequeñas perturbaciones en el seno del fluido:

Fuerzas de tensión superficial, determinadas por:

Para la semejanza dinámica total, el factor de escala de fuerzas debe ser constante,
independientemente del campo de fuerzas considerado:
De la relación anterior se obtienen los siguientes resultados:
(1) relación entre fuerzas de inercia y fuerzas viscosas

ws decir, para la semejanza de los campos de fuerzas de inercia y de esfuerzos viscosos, entre
modelo y prototipo, el numero de REYNOLDS del modelo debe ser el mismo que el del
prototipo.
(2) relación entre fuerzas de inercia y fuerzas gravitatorias

es decir, para la semejanza de los campos de fuerzas de inercia y fuerzas gravitatorias, entre
modelo y prototipo, el numero de FROUDE del modelo debe ser el mismo que el del prototipo.

(3) relación entre fuerzas de inercia y fuerzas elásticas.


es decir, para la semejanza de los campos de fuerzas de inercia y de compresibilidad, entre
modelo y prototipo, el numero de MACH del modelo debe ser el mismo que el del prototipo.

(4) relación entre fuerzas de inercia y fuerzas de tensión superficial:

es decir, para la semejanza de los campos de fuerzas de inercia y de fuerzas de tensión superficial, entre modelo y
prototipo, el numero de WEBER del modelo debe ser el mismo que el del prototipo.

(5) relación entre fuerzas de inercia y fuerzas de presión

es decir, para la semejanza de los campos de fuerzas de inercia y de fuerzas de presión, entre
modelo y prototipo, el numero de EULER del modelo debe ser el mismo que el del prototipo.
Para la semejanza dinámica total entre modelo y prototipo, el factor de escala de fuerzas debe
ser constante, y con ello se tiene las siguientes igualdades entre los parámetros
adimensionales de modelo y prototipo:
Rem=Rep Frm=Frp Mam=Map Wem=Wep Eum=Eup ...

El cumplimiento simultaneo de todas las igualdades anteriores, lleva al absurdo de que el


factor de escala de longitudes sea 1, es decir que el modelo es el propio prototipo; lo que
quiere decir que es imposible ensayar un modelo a escala, y que se conserve la semejanza en
todos los campos de fuerza entre modelo y prototipo; no obstante, en los ensayos con modelos
lo que se hace es considerar cual es el parámetro adimensional controlante, es decir el campo
de fuerzas más importante (a parte del de inercia).

Por ejemplo, el ensayo de arrastre de un modelo de carena en un canal de experiencias


hidrodinámicas, se realiza de tal forma que el numero de FROUDE del modelo sea el del
prototipo, pues en el movimiento de un buque, las fuerzas gravitacionales del campo de olas
generado predomina sobre los esfuerzos viscosos sobre la propia carena. Así si la velocidad
del buque es de 15 nudos, y el modelo se construye a escala 1:100, la velocidad con la
que se debe mover el modelo en el canal será :

NUMEROS ADIMENSIONALES CONTROLANTES EN SEMEJANZA.

A continuación se establecen una pautas cualitativas en el análisis dimensional de


diversos flujos, para establecer que números adimensionales tienen mayor importancia y
cuales se pueden no considerar (sus efectos).
(1) En flujo incompresible, estacionario y sin superficie libre, el número adimensional controlante es el Re. Es
decir se podrá considerar semejanza si el número de Reynolds del modelo es el mismo que el del prototipo.

(2) En flujo incompresible, estacionario y con superficie libre, los números controlantes son el Re y el Fr.
Pero normalmente no es posible mantener los dos números simultáneamente 8 iguales; con lo que en este caso se
considera controlante el número de Froude.

(3) En flujo compresible, estacionario y sin superficie libre, los números controlantes son el Re y el Ma; por
no poder mantener simultáneamente las igualdades, se considera el mas controlante el número de Mach.

(4) En general si el Re es muy grande, se puede no considerar sus efectos, siendo controlantes el Ma o el Fr. No
obstante el Re marca en que tipo de flujo se esta, pero si se esta en flujo turbulento completamente desarrollado,
la influencia del Re es practicamente constante.

(5) Si el flujo es no estacionario se debe considerar además el número de Strouhal. Es importante destacar que
se pueden tener dos tipos de procesos no estacionarios:

- flujo no estacionario con condiciones de contorno estacionarias.


- flujo no estacionario con condiciones de contorno no estacionarias.

El primer caso de condición de contorno estacionaria, se tiene en la formación de la


estela de torbellinos de Von Karman, en el flujo externo de la interacción de un flujo con un
objeto. Aguas arriba del objeto, el flujo puede ser bastante uniforme, pero aguas abajo se
origina una estela que integra una serie de vortices contrarrotantes, que van creciendo en el
sentido del flujo, hasta alcanzar un tamaño crítico, y son arrastrados por la propia corriente: el
crecimiento y desprendimiento de vórtices es alternativo y no estacionario. La frecuencia de
desprendimiento de vórtices es función de Re, y viene determinada por el número de Strouhal:
El segundo caso de condición de contorno no estacionaria, se tiene cuando algunos de
los contornos geométricos con los que interacciona el flujo, tiene posiciones no estacionarias
(dependen del tiempo).

Un ejemplo clásico es la destrucción del puente de TACOMA por la interacción del


viento. En condiciones normales, la estructura esta estacionario (no se mueve), pero para
determinadas condiciones del viento, la frecuencia de los torbellinos coincidió con una
frecuencia natural de la estructura, que empezó a vibrar (contorno no estacionario), generando
nuevas frecuencias y amplitudes de los torbellinos de Von Karman, que resonaron con las
frecuencias naturales de la estructura, dando lugar a un proceso de auto excitación, que llevo a
la destrucción del puente.

Ejercicios

1. Ensayo de modelos de carenas: números de Froude y de Reynolds. La resistencia al avance de una carena,
viene determinada por las fuerzas de superficie debidas a la presión y a la viscosidad, sobre la propia
superficie mojada de la carena. En cuanto a la presión, se tienen presiones hidrostáticas, generadas por el
campo de olas. Para mantener semejanza dinámica total, entre modelo y prototipo, es necesario que los
campos de fuerzas sean semejantes, lo que lleva: en el caso de las fuerzas de rozamiento viscoso, a que el Re
de modelo y prototipo sean iguales, y en el caso de las fuerzas de presión hidrostática, a que el Fr de modelo
y prototipo sean iguales. Por la imposibilidad experimental de que simultáneamente se cumpla que FrM=FrP
y ReM=ReP, los ensayos, se deben realizar con una única condición, que en el caso de carena de buques, es el
Fr, por ser el campo gravitatorio el predominante, en la resistencia al avance.

En el ensayo de un modelo de carena, se realizan ensayos en el canal de experiencias,


arrastrando el modelo a determinada velocidad (marcada por el Fr), y se mide la resistencia al
avance.

DETERMINE para cada punto de funcionamiento de la carena prototipo:


1. Coeficientes de resistencia para el modelo: total, friccional y residual.
2. Coeficientes de resistencia para el prototipo: residual, friccional y total.
3. Resistencia al avance del prototipo para las distintas velocidades del buque.

DATOS: Prototipo: Eslora: LP = 103,77 m; Superficie mojada: AmP = 3103 m2;


Modelo: factor de escala de longitudes: Z = 13,02

1 2 3 4 5 6 7 8 9
vM(m/s) 1,568 1,640 1,711 1,782 1,853 1,925 1,996 2,067 2,139
RTM(N) 97,9 107,3 118,1 130,5 144,8 161,9 182,6 207,9 238,4

propiedades físicas: aceleración gravedad: g = 9,80665 m/s2


densidad agua de mar (15ºC): qP=1025,9 kg/m3
viscosidad cinemática agua de mar: vP=1,1882
10-6 m2/s densidad agua del canal (17,6ºC): qM
= 998,6 kg/m3
viscosidad cinemática agua del canal : vM = 1,0653 10-6 m2/s

RESOLUCIÓN: Ecuaciones de cálculo de coeficientes de resistencia:


La eslora y la superficie mojada del modelo son:

(1) COEFICIENTES PARA MODELO: para cada uno de los puntos de funcionamiento, se obtiene el coeficiente
de resistencia total, a partir de la medida de la resistencia al avance; el coeficiente de resistencia friccional, por la
ecuación ITTC-57; y finalmente el coeficiente residual por diferencia entre ambos:

punto de funcionamiento 1:

Análogamente, se obtienen los valores para el resto de puntos de funcionamiento:


1 2 3 4 5 6 7 8 9
vM(m/s) 1,568 1,640 1,711 1,782 1,853 1,925 1,996 2,067 2,139
RTM(N) 97,9 107,3 118,1 130,5 144,8 161,9 182,6 207,9 238,4
CRTM 4,360E-03 4,368E-03 4,417E-03 4,499E-03 4,617E-03 4,783E-03 5,018E-03 5,328E-03 5,705E-03
ReM 1,173E+07 1,227E+07 1,280E+07 1,333E+07 1,386E+07 1,440E+07 1,493E+07 1,546E+07 1,600E+07
CRFM 2,919E-03 2,896E-03 2,875E-03 2,856E-03 2,837E-03 2,819E-03 2,801E-03 2,785E-03 2,769E-03
CRRM 1,441E-03 1,472E-03 1,541E-03 1,644E-03 1,780E-03 1,965E-03 2,217E-03 2,542E-03 2,936E-03

(2) COEFICIENTES PARA PROTOTIPO: para cada uno de los puntos de funcionamiento, la determinación de
los coeficientes de resistencia es:

Coeficiente de resistencia RESIDUAL: por dependen del Fr, y ensayar con FrM=FrP, los
coeficientes residuales serán iguales. La relación entre las velocidades de prototipo y modelo,
vienen determinadas por esta igualdad del numero de Froude:

Coeficiente de resistencia FRICCIONAL: por depender del Re, y no poder ensayar a Re


constante:

se obtiene a partir de la determinación del Re.

Coeficiente de resistencia TOTAL: es la suma de los dos anteriores. Para el punto de


funcionamiento 1:

Análogamente, para el resto de puntos de funcionamiento del prototipo:

1 2 3 4 5 6 7 8 9
VP (m/s) 5,658 5,918 6,174 6,430 6,686 6,946 7,202 7,458 7,718
CRRP 1,441E-03 1,472E-03 1,541E-03 1,644E-03 1,780E-03 1,965E-03 2,217E-03 2,542E-03 2,936E-03
ReP 4,941E+08 5,168E+08 5,392E+08 5,616E+08 5,839E+08 6,066E+08 6,290E+08 6,514E+08 6,741E+08
CRFP 1,674E-03 1,664E-03 1,655E-03 1,646E-03 1,638E-03 1,630E-03 1,623E-03 1,615E-03 1,608E-03
CRTP 3,115E-03 3,136E-03 3,196E-03 3,290E-03 3,419E-03 3,595E-03 3,839E-03 4,158E-03 4,544E-03
(3) RESISTENCIA TOTAL DEL PROTOTIPO: Para cada punto de funcionamiento, conociendo su coeficiente
de resistencia total, se tiene:

para el punto de funcionamiento 1:

Análogamente para el resto de puntos de funcionamiento:


1 2 3 4 5 6 7 8 9
VP (m/s) 5,658 5,918 6,174 6,430 6,686 6,946 7,202 7,458 7,718
RTP (kN) 158,710 174,783 193,926 216,521 243,252 276,076 316,978 368,14 430,85

EN LA TABLA SIGUIENTE, SE RESUMEN TODOS LOS RESULTADOS:

1 2 3 4 5 6 7 8 9

M O D E L O
vM(m/s) 1,568 1,640 1,711 1,782 1,853 1,925 1,996 2,067 2,139
RTM(N) 97,9 107,3 118,1 130,5 144,8 161,9 182,6 207,9 238,4
FrM 0,176 0,184 0,192 0,202 0,210 0,217 0,226 0,235 0,243
ReM 1,173E+07 1,227E+07 1,280E+07 1,333E+07 1,386E+07 1,440E+07 1,493E+07 1,546E+07 1,600E+07
CRTM 4,360E-03 4,368E-03 4,417E-03 4,499E-03 4,617E-03 4,783E-03 5,018E-03 5,328E-03 5,705E-03
CRFM 2,919E-03 2,896E-03 2,875E-03 2,856E-03 2,837E-03 2,819E-03 2,801E-03 2,785E-03 2,769E-03
CRRM 1,441E-03 1,472E-03 1,541E-03 1,644E-03 1,780E-03 1,965E-03 2,217E-03 2,542E-03 2,936E-03

P R O T O T I P O
VP(NUD0S) 11,0 11,5 12,0 12,5 13,0 13,5 14,0 14,5 15,0
VP (m/s) 5,658 5,918 6,174 6,430 6,686 6,946 7,202 7,458 7,718
RTP(kN) 158,71 174,78 193,93 216,52 243,25 276,08 316,98 368,14 430,85
PARRASTRE(kW) 897,96 1034,31 1197,27 1392,24 1626,44 1917,63 2282,94 2745,77 3325,36
FrP 0,176 0,184 0,192 0,202 0,210 0,217 0,226 0,235 0,243
ReP 1,173E+07 1,227E+07 1,280E+07 1,333E+07 1,386E+07 1,440E+07 1,493E+07 6,514E+08 6,741E+08
CRTP 3,115E-03 3,136E-03 3,196E-03 3,290E-03 3,419E-03 3,595E-03 3,839E-03 4,158E-03 4,544E-03
CRFP 1,674E-03 1,664E-03 1,655E-03 1,646E-03 1,638E-03 1,630E-03 1,623E-03 1,615E-03 1,674E-03
CRRP 1,441E-03 1,472E-03 1,541E-03 1,644E-03 1,780E-03 1,965E-03 2,217E-03 2,542E-03 2,936E-03
2. Fuerza de arrastre: números de Reynolds y de Euler. La fuerza de arrastre sobre un dirigible debe
determinarse mediante ensayos con un modelo a escala reducida. Los ensayos se realizan en un túnel de
agua, en donde el modelo de dirigible se fija a un soporte instrumentado con un sensor de fuerza; la
velocidad del flujo a agua que incide sobre el modelo, se controla con una bomba de velocidad variable que
es la que origina el flujo hacia el modelo.

DETERMINE: 1. Cualitativamente los números adimensionales controlantes.


2. Justifique porque se tiene que ensayar el modelo en un túnel de agua.
3. Velocidad a la que se debe ensayar el dirigible modelo en el túnel de agua.
4. Potencia de arrastre del dirigible prototipo.

DATOS: Dirigible: Lp=20m; Up=36 km/h;


aire a 25ºC y 1 bar: viscosidad = 15,58 10-6 m2/s; densidad
= 1,2 kg/m3 Modelo: Lm=1m; agua, viscosidad = 10-6 m2/s;
3
densidad = 1000 kg/m
Fuerza de arrastre a la velocidad de ensayo = 1649 N

RESOLUCIÓN:

(1) PARÁMETROS ADIMENSIONALES CONTROLANTES.

Los números adimensionales más importantes en un proceso de flujo son: Re, Eu, Ma y Fr.
Debido a la baja velocidad del submarino con respecto a la velocidad del sonido, los efectos
de compresibilidad no son importantes, y en número de MACH no es controlante; por no haber
superficies libre, el número de Froude tampoco es controlante; con lo que los números
adimensionales que se tienen en cuenta para la similitud dinámica son exclusivamente el
número de REYNOLDS y el número de EULER.

Se puede llegar analíticamente a esta conclusión, partiendo de que la fuerza de arrastre


depende de las siguientes variables: FD=f(U,L,q,µ), y aplicando el teorema de BUCKINGHAM, se
obtienen dos parámetros adimensionales, que están relacionados, respectivamente, con los
números de REYNOLDS y de EULER.

(2) ENSAYO EN UN TUNEL DE AGUA.

Para la semejanza dinámica entre modelo y prototipo, sus números de Reynolds deben ser iguales:

en un dirigible la longitud característica en su ESLORA.

Con lo que la velocidad a la que se debe mover el modelo es:


Numéricamente, con los datos: Up = 36(1000/3600)= 10 m/s
Lp = 20 m;
Lm = 1 m
vp = 10-6 m2/s;
vmagua = 15,58 10-6 m2/s

Si el ensayo se realiza en un túnel de viento, la velocidad del aire en el túnel debería ser de:

Esta velocidad tiene dos problemas: uno el experimental, porque se tiene que trabajar con
grandes velocidades de flujo de aire y sobre todo porque se tiene números de Mach
relativamente altos (Ma>0,3), y se deberían considerar los efectos de la compresibilidad. Para
evitar estos dos problemas, se ensaya en un túnel de agua, debido a la que la viscosidad
cinemática del agua es unas 15 veces menor que la del aire, y con ello la velocidad de
ensayo en agua baja en la misma proporción.

(3) VELOCIDAD DEL FLUJO QUE INCIDE SOBRE EL MODELO:

La velocidad del flujo uniforme que se hace incidir sobre el modelo, viene determinada por la
igualdad de números de Reynolds entre modelo y prototipo:

(4) POTENCIA DE ARRASTRE:

La potencia de arrastre, viene determinada por: PD = FD U.

El ensayo suministra el valor de la fuerza de arrastre sobre el modelo. La fuerza de arrastre


sobre el prototipos se determina por análisis dimensional, que establece que para la semejanza
dinámica los números de EULER, deber ser iguales:
Numéricamente con los datos: FDm=1649 N; qp = 1,2 kg/m3 qm = 1000 kg/m3
Up = 10 m/s Um = 12,84 m/s
Lp = 20 m Lm = 1 m

Con lo que la potencia de arrastre, estimada para el prototipo del dirigible es:

PDp= FDp Up = (480,1 N)(10m/s) = 4,8 kW

La potencia de arrastre involucrada en el modelo a escala reducida es:

PDm=FDmUm=(1649N)(12,84m/s) =21,17kW

No debe extrañar, que tanto la fuerza de arrastre, como la potencia de arrastre, sean mayores
en el modelo a escala reducida, que en el prototipo, debido a que el modelo se ensaya en agua
y el prototipo se mueve en aire.

3. Cavitación: número de Euler (número de cavitación) En el flujo de líquidos, si por aumento de velocidad,
la presión baja hasta valores próximos a la presión de vapor (a la temperatura del flujo), el líquido cambia
de fase, se generan burbujas de vapor, que son arrastradas por la corriente hacia zonas de alta presión, en
donde las implosiones generan pulsos de alta frecuencia y alta presión. A este fenómeno se denomina
cavitación. Consideremos el flujo entorno a un arpón de pesca submarina : conforme el disparo se realiza
a menos profundidad, la velocidad para condiciones de cavitación es menor. Se realizan un ensayo en
donde se determina que a una profundidad de 2 m, si la velocidad es menor de 44 m/s no hay cavitación.

DETERMINE: 1. Número de cavitación mínimo para no cavitación.


2. Velocidad máxima para no cavitación a una profundidad de10 m.
3. A que profundidad deja de cavitar, si el arpón se dispara verticalmente hacía abajo, desde un
profundidad de 2 m a la velocidad del apartado anterior

DATOS: agua de mar: densidad=1025 kg/m3; presión de vapor = 23 mbar


Arpón: a 2 m de profundidad la velocidad de no
cavitación es 44 m/s CD=0,4 (suponerlo constante)
Peso = 0,250 kg; densidad = 2300 kg/m3;
Área recta máxima =2 cm2
Presión atmosférica: 990 mbar

RESOLUCIÓN:

Se define como número de cavitación, al número de Euler, derivado de considerar la


diferencia, entre la presión de vapor y la presión no perturbada por la interacción del objeto
que se mueve en el seno del líquido. En nuestro caso:

en donde ‘’h’’ es la profundidad


(1) NÚMERO DE CAVITACIÓN MÍNIMO:

Con el resultado experimental, de que a 2 m de profundidad la cavitación se inicia a


velocidades mayores de 12 m/s; el número de cavitación mínimo para no cavitación es:

y se denomina número de cavitación crítico

(2) VELOCIDAD MÁXIMA PARA NO CAVITACIÓN A UNA PROFUNDIDAD DE 10 m:

la condición de no cavitación es que el número de cavitación sea mayor que el crítico:

a 10 metros de profundidad, la velocidad debe ser:

(3) PROFUNDIDAD A LA QUE DEJA DE CAVITAR, CON UN DISPARO VERTICAL DESCENDENTE, EN LAS
CONDICIONES: h 0= 2m; U0 = 57,1 m/s.

La ecuación de equilibrio de fuerzas sobre el arpón es:


con todo se tiene la ecuación diferencial:

D2

como lo que interesa obtener en la ecuación h=h(v); en la ecuación se puede expresar el elemento de
tiempo, en función de la velocidad y de la variación de profundidad: dt=dh/v; con todo se tiene la
ecuación diferencial:

en donde las constantes son:

cuya integración da la ecuación h=h(v):

obteniedose también la función v=v(h)

Numericamente:

Con lo que la velocidad en función de la profundidad es:

el número de cavitación en función de la profundidad es:

Con las dos últimas ecuaciones se construye la tabla : h -v -Ca:

profundidad velocidad Nº cavitación


m m/s -
2,00 57,10 0,070
2,50 52,65 0,086
3,00 48,56 0,105
3,10 47,78 0,109
3,20 47,01 0,114
3,30 46,26 0,118
3,50 44,79 0,128

La zona de cavitación esta en los puntos de número de cavitación menor de 0,118 (Ca<Ca*);
es decir, el arpón sale en cavitación (Ca=0.070) y al alcanzar la profundidad de 3,3 m deja de
cavitar.

4. Campo de olas en dinámica de buques: número de Froude. En la dinámica de buques, se


estudia el efecto de un campo de olas, sobre el comportamiento de la nave. Para ello se
realizan ensayos en canales que incorporan un generador de olas; una vez establecido el
campo de olas, se arrastran las carenas modelo y se miden los esfuerzos. Para un
determinado buque, se quiere obtener el par escorante al que le sometería un campo de olas
de 8 metros de altura y un periodo de 20 segundos. Para ello se construye un modelo que se
ensaya en el canal de dinámica del buque. A partir de los datos.

DETERMINE: 1. Números adimensionales que controlan la hidrodinámica de vehículos marinos.


2. La velocidad a la que hay que mover el modelo.
3. El campo de olas que hay que generar: amplitud y periodo.
4. El factor de escala de momentos escorantes entre modelo y prototipo.

DATOS: Buque prototipo: Eslora = 200 m; velocidad 14 nudos;


agua de mar: 1025kg/m3
Carena modelo: Eslora = 4m; agua dulce: 1000 kg/m3

RESOLUCIÓN:
(1) NUMEROS ADIMENSIONALES CONTROLANTES:

En hidrodinámica de vehículos marinos, los números adimensionales controlantes son el


número de Froude (Fr) y el número de Strouhal (St).

En el movimiento de buque en un campo de olas externo el número adimensional controlante, es el


que controla
las amplitudes, periodos y frecuencias del oleaje, es decir el número de Strouhal:
En el movimiento de un buque, el campo de olas generado en su movimiento, por
modificación de la superficie libre, es de tal importancia, que la mayor parte de la resistencia
al avance, es debida a este campo de olas. Por ello, el parámetro adimensional que controla
los fenómenos de superficies libres en campo gravitacional, es el número de Froude, que sería
controlante.

En principio los efectos de la viscosidad también deberían considerarse, es decir el número de


Reynolds sería también controlante; con lo que los ensayos de modelo, se deberían realizar de
tal forma que se cumpliese: Frm=Frp y Rem=Rep; para asegurar la semejanza dinámica. Pero
esta doble igualdad, lleva a un absurdo:

La igualdad de Froude lleva a que la velocidad de arrastre del modelo sea:

La igualdad de Reynolds lleva a que la velocidad de arrastre del modelo sea

Para que las dos expresiones, coincidan, se tiene que:

En nuestro caso, la escala geométrica es Lm/Lp = 1/50, con lo que la viscosidad del fluido
utilizado en el canal de ensayos, debería ser 1000 veces menor que la del agua de mar(
(1/50)3/2=1/354). Evidentemente, la gran cantidad de fluido que se tiene que utilizar en el
canal de ensayos, hace inviable la búsqueda de un fluido de ensayo de tal viscosidad, por lo
que evidentemente se usa agua dulce como fluido de ensayo.

Así pues, usando agua dulce como fluido de ensayo, la igualdad de Renolds, lleva a
velocidades de ensayos mucho mayores que las correspondientes a considerar la igualdad de
Froude.

No obstante, este problema se resuelve, debido a que en el movimiento de un buque, el Re es


poco controlante, frente al Fr, con lo que el ensayo en el canal de dinámica del buque, se
realiza a Froude constantes.

(2) VELOCIDAD A LA QUE HAY QUE ARRASTRAR EL MODELO.

Por la igualdad de Fr en el modelo y el prototipo, la velocidad de ensayo es:


(3) CAMPO DE OLAS:

El modelo y el prototipo deben moverse en campos de olas semejantes cinemáticamente, con


lo que sus números de Strouhal deben ser iguales:

Con lo que el periodo de las olas a generar es:

Y la amplitud de las olas generadas debe ser:

(4) FACTOR DE ESCALA DE MOMENTOS:

Por el Buckingham se puede obtener que hay el parámetro adimensional, con el que esta
relacionado el momento de escora es:

Con lo que por semejanza dinámica entre modelo y prototipo, ambos parámetros
adimensionales de momento, deben ser iguales:

Lo que lleva a obtener el factor de escala de momentos:


3.5 Frecuencias características: número de Strouhal. En el diseño de aletas de submarinismo, se intenta realizar un
prototipo, a partir de los datos obtenidos de la propulsión de un delfín. Al tipo de propulsión que utilizan los peces, se
le denomina propulsión por estela de remolinos, en donde el empuje conseguido depende de la amplitud del coleteo, de
la velocidad de avance y de la densidad del agua. A partir de los datos:

DETERMINE: 1. Parámetros controlantes de la propulsión por estela de remolinos.


2. Amplitud con la que deben moverse las aletas del submarinista, para tener un flujo
dinámicamente semejante al coleteo del delfín.
3. Frecuencia de aleteo del submarinista.

DATOS: Delfín: velocidad = 5 m/s; coleteo: amplitud = 1m;


frecuencia = 1,5 Hz
Submarinista: velocidad = 1 m/s
Hidrodinámica: CDsub = 3 CDdelfín; Asub = 0,5 Adelfín

RESOLUCIÓN:

(1) PARAMETROS CONTROLANTES:

En procesos oscilatorios y/o no estacionarios, el parámetro controlante, es el factor de escala de


tiempos, o número de STROUHAL:

en donde “f” es una frecuencia característica,“T” su periodo, L una longitud característica, U una
velocidad característica y L/U el tiempo de residencia.

El otro parámetro controlante, viene determinado por la función T=T(L,U,q); se obtiene por el
teorema de Buckingham, se denomina coeficiente de empuje KT, y deriva del número de EULER:

en donde: “T” es la fuerza de empuje suministrada por el sistema de propulsión (del inglés
“thrust”)

“L” es una longitud característica. En el caso de hélices marinas, es el diámetro


de la hélice, y en este caso es la amplitud de los coleteos-aleteos.

(2) AMPLITUD DE LAS OSCILACIONES:

De la igualdad de coeficientes de empuje de modelo y prototipo, se obtiene el factor de escala de


longitudes:

evidentemente las densidades son iguales, porque tanto el delfín como el submarinista se mueven
en el mar.
El factor de escala de fuerzas Td/Ts, se obtiene a partir de las fuerzas de arrastre que se oponen al
movimiento:

Con las ecuaciones [1] y [2], se obtiene el factor de escala de longitudes:

Como la amplitud de los coleteos del delfín es de 1m, la amplitud de los aleteos del submarinista,
será de:

(3) FRECUENCIA DE LOS ALETEOS DEL SUBMARINISTA:

Por la igualdad de números de Strouhal:

La frecuencia del movimiento de las aletas del submarinista debe ser de 0,245 Hz, es decir un
aleteo cada 0,245 segundos o bien 1/0,245 = 4 segundos /aleteo, que es el periodo del movimiento
de las aletas.

a. Aplicación de la teoría de modelos: diseño de un experimento para el


estudio de impresión por chorro de tinta.

En la impresión por chorro de tinta, se tienen impactos de gotas muy pequeñas sobre
superficies planas. La sección de contacto entre la gota de tinta y la superficie del papel, es
una función del tiempo. Debido al tamaño de las gotas (volumen del orden de nanolitros y
diámetros del orden de décimas de milímetros ), es imposible la visualización experimental de
la dinámica del proceso; por ello se planifica un experimento, con gotas mas grandes (del
orden de milímetros), y se mide el tamaño de contacto en diversos instantes correlativos de
tiempo. A partir de los resultados experimentales para las gotas modelo:

DETERMINE: 1. Números adimensionales controlantes.


2. Escala de longitudes, velocidades y tiempos.
3.Diámetro y velocidad de impacto de las gotas modelo.
4. Evolución temporal del tamaño de contacto de las gotas de tinta.
DATOS: Gotas de tinta: tamaño: Dp=0,2 mm ; Velocidad de impacto: Vp=1m/s
densidad: q=1040 kg/m3; viscosidad: µ=5cP; tensión superficial en
aire: o = 0,04 N/m
Líquido de ensayo: densidad = 1000 kg/m3; viscosidad: 18cP,
tensión superficial en aire = 0,05N/m

Experimentales: Φ vs t ; Φ= diámetro del área mojada;

t=tiempo

Φ(mm) 0 2,1 4,2 5,9 7,6 7,1 6,2 5,3 4,5 4,1 3,2 2,7
t(ms) 0 0,5 1 2 3 4 5 6 7 8 9 10

tiempo

RESOLUCIÓN:
(1) NÚMEROS ADIMENSIONALES CONTROLANTES:

Las variables que intervienen son: Diámetro de las gotas: D [D] = L


Velocidad de impacto de las gotas: U [U] = LT-1
Densidad gotas: q [q] = ML-3
Viscosidad gotas: µ [µ] = ML-1T-1
Tensión superficial gota-aire: o [o] = MT-2

Por el teorema de Buckingham, con 5 variables y 3 dimensiones distintas, se obtienen 5-3=2


parámetros adimensionales: como variables independientes se consideran la viscosidad y la
tensión superficial.

Se procede a hallar el primer numero adimensional:

número adimensional = número de REYNOLDS:

Ahora se procede a hallar el segundo número adimensional:


número adimensional = número de WEBER:

(2) FACTORES DE ESCALA: LONGITUDES (NL), VELOCIDADES (NV ) y TIEMPOS


(NT) :

Para semejanza entre las gotas modelo y las gotas prototipo, se debe cumplir:

Obteniéndose:

El factor de escala de tiempos viene determinado, a partir de los factores de escala de velocidades y
longitudes:
Con los datos numéricos: ρp = 1044 kg/m3; ρm = 1000 kg/m3

µp = 5 cP µm = 17,3 cP

σp = 0,04 N/m σm = 0,05

se obtienen los siguientes factores de escala:

Por lo tanto:

(3) TAMAÑO Y VELOCIDAD DE IMPACTO DE LAS GOTAS MODELO:

A partir del tamaño y velocidad de impacto de las gotas de tinta (prototipo), y de los factores
de escala de longitudes y de velocidades, se obtienen el tamaño y la velocidad de las gotas de
líquido, con las que se realizan los experimentos:

Dm = NL Dp = (10) (0,2 mm) = 2 mm

Vm = NV Vp = (0,361)(1 m/s) = 0,361 m/s

(4) EVOLUCIÓN TEMPORAL DEL TAMAÑO DE LA SUPERFICIE MOJADA POR LA TINTA:

A partir de la tabla de valores experimentales para las gotas modelo, se obtiene una tabla
para valores de las gotas de tinta, a partir de los factores de escala de longitudes y de
tiempos:

Φp = Φm/NL = $m/10

tp = tm / NT = tm / 27,7
Φmodelo(mm) 0 2,1 4,2 5,9 7,6 7,1 6,2 5,3 4,5 4,1 3,2 2,7
tmocledo(ms) 0 0,5 1 2 3 4 5 6 7 8 9 10
Las unidades de tiempo para el modelo son milisegundos

$prototipo(mm) 0 0,21 0,42 0,59 0,76 0,71 0,62 0,53 0,45 0,41 0,32 0,27
tprototipo(µs) 0 18,05 36,1 72,2 108,3 144,4 180,5 216,6 252,7 288,8 324,9 361
Las unidades de tiempo para el prototipo son microsegundos

EVOLUCIÓN DEL TAMAÑO


0,8 DE LA ZONA MOJADA POR
LA TINTA, EN FUNCIÓN
0,7 DEL TIEMPO
Tamaño
zona mojada 0,6
(milímetros)
0,5

0,4

0,3

0,2 Tamaño gota de tinta antes del impacto = 0,2 mm

0,1

0
0 50 100 150 200 250 300 350 400

Tiempo (microsegundos)

b. Reglas de semejanza en bombas y ventiladores: El flujo en una bomba depende


fundamentalmente del caudal, de la velocidad de giro y del tamaño del rodete; además de las
propiedades del fluido (viscosidad y densidad). Los parámetros de funcionamiento de interés son:
el aumento de presión “ “Op”” (altura o carga en unidades de longitud: H = Op/qg), la potencia
consumida ”“P”” y el rendimiento “ “y”” (relación entre la potencia entregada al líquido y la consumida:
y=(Q Op)/P = qQgH/P). La energía especifica (por unidad de masa) que se transfiere al líquido es el
termino ”“gH””, y es el parámetro que normalmente se utiliza en análisis dimensional en bombas: A
partir de la información anterior y de los datos:

DETERMINE: 1. Parámetros adimensionales, a partir del teorema de Buckinham.


2. Parámetro adimensional independiente del tamaño de la bomba: velocidad específica.
3. Reglas de semejanza entre máquinas geométrica y dinámicamente semejantes:
4. Punto de funcionamiento cuando la velocidad de giro aumenta un 20%.

DATOS: Bomba: D = 127 mm; Q = 18 L/s; q =1000 kg/m3; N = 1450 rpm; H = 6,8 m; y = 63%

centrifuga o radial mixta axial

Ns = 0,2 0,4 0,6 0,8 1,0 2,0 3,0 4,0

Clasificación de bombas en función de la velocidad específica


RESOLUCIÓN:

(1) PARÁMETROS ADIMENSIONALES:

En el flujo de bombas, las variables son:


Energía específica gH [gH] = L2 T-2
Caudal Q [Q] = L3 T-1
Potencia consumida P [P] = M L2
T-3
Diámetro del rodete D [D] = L
Velocidad de giro x [x] = T-1
Densidad del líquido q [q] = M L-3
Viscosidad del líquido µ [µ] = M L-1 T-1

Consideraremos como variables independientes: la densidad, la velocidad de giro y el diámetro;


con lo que los 4 parámetros adimensionales que se tienen (7 variables – 3 dimensiones = 4
parámetros) son:

El parámetro adimensional del aumento de energía específica del líquido, se denomina


coeficiente de altura, y determina que para una misma bomba, si la velocidad de giro se
duplica, la altura de bombeo de cuadruplica; y que la altura depende del cuadrado del tamaño
del rodete.

El parámetro adimensional de la potencia consumida, se denomina coeficiente de potencia, y


determina que la potencia depende del cubo de la velocidad de giro y de la quinta potencia
del tamaño del rodete.

El parámetro adimensional del caudal, se denomina coeficiente de caudal, y determina que el


caudal aumenta linealmente con la velocidad, y cúbicamente con el tamaño:
El cuarto parámetro adimensional es la inversa del número de Reynolds:

(2) PARÁMETRO ADIMENSIONAL INDEPENDIENTE DEL TAMAÑO: VELOCIDAD ESPECIFICA.

En bombas, hay un parámetro adimensional, que se obtiene de los anteriores, que no depende
del tamaño de la máquina, sino de su morfología: es la velocidad especifica Ns, que se obtiene
por la relación entre los parámetros adimensionales (n3)1/2 / (n1)3/4:

El valor de la velocidad característica, permite clasificar las bombas:

Bombas centrífugas o radiales Ns bajos (0,2 a 0,5) bajos caudales y altas cargas
Bombas de flujo mixto Ns intermedio (0,5 a 3) caudales y alturas intermedias
Bombas axiales Ns alto (3 a 4) altos caudales y bajas cargas

En la ecuación de la velocidad característica, las variables se expresan en el S.I., es decir la


velocidad de giro en radianes / segundo, el caudal en m3/s, la altura o carga en metros y la
aceleración gravitatoria en m/s2. No obstante es habitual, utilizar dos tipos de sistemas
técnicos, el europeo ( en donde la velocidad de giro se expresa en Herzios
(Hz=revoluciones/segundo); y el americano (en donde la velocidad de giro se expresa en rpm
(revoluciones/minuto), el caudal en gpm (galones / minuto), y en lugar de la energía especifica,
se pone la carga en pies; las relaciones entre las distintas velocidades características es:

(3) REGLAS DE SEMEJANZA: si se disponen de dos bombas geométrica y dinámicamente semejantes, y los
viscosos son despreciables (por trabajas a Re altos), se tiene que el coeficiente de caudal es el mismo en las dos
bombas; y además como los coeficientes de altura y de potencia solo dependen (despreciando efectos viscosos)
del coeficiente de caudal, se tiene también las igualdades correspondientes de los coeficientes de altura y de
potencia:

por similitud geométrica y cinemática: CQ1=CQ2


(1) q
por trabajar a Re muy altos:

(2)

(3)

Las expresiones [1], [2] y [3], se denomina leyes de semejanza de las bombas.

(4). PUNTO DE FUNCIONAMIENTO CON UN AUMENTO DEL 20% DE VELOCIDAD DE


GIRO.

En este caso, se trata de la misma bomba y la semejanza geométrica es total: en este caso D1=D2.

Si lo que varia es la velocidad de giro, el nuevo punto de funcionamiento viene determinado


por las leyes de semejanza anteriores:

Punto de funcionamiento “1”: N1 = 1450 rpm

Punto de funcionamiento “2”: N2 = 1,2 (1450)=1740 rpm


BIBLIOGRAFIA

• CENGEL, C. “Mecánica de fluidos fundamentos y aplicaciones" Ed. McGraw Hill ( 2006


)
• SHAMES , I. “La mecánica de los fluidos" Ed. McGraw Hill ( 1995)
• WHITE , F. “Mecánica de fluidos" Ed. McGraw Hill ( 2008 )
• CRESPO , A. “Mecánica de fluidos” Ed. Thomson ( 2006 )
• BARRERO RIPOLL, A., PÉREZ-SABORID SÁNCHEZ-PASTOR, M. ‘’Fundamentos y
aplicaciones de la Mecánica de Fluidos’’ Ed. McGraw Hill ( 2005 )
• GILES. R.V. “Mecánica de fluidos e hidráulica” Ed. McGraw Hill ( 1994 )
• STREETER , V. “Mecánica de los fluidos" Ed. McGraw Hill ( 1989 ) 532 STR mec
• WEBBER , N.B. “Mecánica de fluidos para ingenieros" Ediciones Urmo ( 1969

También podría gustarte