Introducción al flujo turbulento

1. Caracterı́sticas fı́sicas de los flujos turbulentos

Antes de definir en una forma precisa el significado de la turbulencia, resulta mas provechoso
para comprender su naturaleza analizar de qué formas se manifiesta y las caracterı́sticas esen-
ciales de esas formas. Las soluciones clásicas de las ecuaciones de Navier-Stokes, corresponden
a los flujos laminares, caracterizados por ser flujos ordenados y regulares. Pero los flujos que
se presentan en situaciones de la vida real no responden a esas soluciones porque manifiestan
un comportamiento completamente diferente. Las manifestaciones tı́picas de la turbulencia se
pueden resumir en tres efectos principales:

Inestacionaridad de pequeña escala

Un permanente y continuo proceso de mezcla

Un dominante efecto de la inercia del fluido

1.1. Inestacionaridad
Aún cuando se encuentre confinado y limitado por contornos estacionarios el flujo turbulento
presenta una progresiva sucesión de inestabilidades que causan e inducen fluctuaciones en las
componentes de la velocidad que hacen que el movimiento se torne inestacionario. Este tipo
de manifestaciones pueden ser observados en experimentos com la conocida experiencia de
Reynolds, en chorros fluidos y capas lı́mites.

1.2. Procesos de mezcla

La turbulencia introduce componentes de velocidad transversales el movimiento principal
que generan procesos de mezcla que aumentan notablemente el intercambio de masa, calor y
cantidad de movimiento en la masa fluida.

1.3. Efecto dominante de la inercia del fluido

La densidad del fluido tiene un papel mucho mas preponderante que el coeficiente de vis-
cosidad molecular en el establecimiento de las tensiones tangenciales efectivas. Una remarcable
caracterı́stica de la turbulencia es que la disipación friccional es mucho menos dependiente de
la viscosidad que de la densidad del fluido. Esto se debe a que la transferencia de energı́a
mecánica desde la componente estacionaria del flujo a la energı́a. Mientras que en los flujos
laminares las tensiones tangenciales son producidas exclusivamente por la viscosidad molecular
y esta domina sobre los efectos de inercia, en los flujos turbulentos la viscosidad molecular tiene
una muy débil influencia, y las tensiones tangenciales (de naturaleza diferente a la molecular)
son gobernadas casi exclusivamente por los efectos inerciales, asociadas a esa inestacionariedad
de pequeña escala presente en los flujos turbulentos.nética de la turbulencia es provocada por
efectos inerciales.
 Figura 1: Distribución de velocidades en flujos laminar y turbulento

Ciertos aspectos de los flujos turbulentos y su contraste con los laminares pueden observarse
mediante la realización de experimentos en un conducto con flujo totalmente desarrollado.
Algunas caracterı́sticas relevantes que pueden apreciarse de las observaciones son las siguientes:

Mientras un instrumento estándar para medir la velocidad como una sonda Pitot-estática
revela una distribución de velocidades en cada punto que parece no cambiar, la utilización
de un instrumento de respuesta rápida como un anemómetro de hilo caliente ( de muy
pequeño diámetro y consecuentemente baja inercia térmica), mostrará que la velocidad
en un punto variará en el tiempo, indicativo que la velocidad local en realidad es inesta-
cionaria.

Los flujos laminares tienen una distribución de velocidades tı́picamente parabólica, mien-
tras que los perfiles de velocidad medios en un flujo turbulento presentan un núcleo central
más aplanado, debido a un rápido intercambio de cantidad de movimiento de las partı́cu-
las, con un gradiente de velocidades sobre la pared mayor que provoca mayores tensiones
tangenciales sobre la misma.

Mientras que en los flujos laminares las tensiones tangenciales son producidas exclusi-
vamente por la viscosidad molecular y esta domina sobre los efectos de inercia, en los
flujos turbulentos la viscosidad molecular tiene una muy débil influencia, y las tensiones
tangenciales (de naturaleza diferente a la molecular) son gobernadas casi exclusivamente
por los efectos inerciales, asociadas a esa inestacionariedad de pequeña escala presente en
los flujos turbulentos.

Un alto grado de desorden en el flujo y movimientos inestacionarios e irregulares super-

puestos a un flujo regular y ordenado que produce intercambio y mezcla de masa, cantidad
de movimiento y energı́a entre las capas fluidas.

El conjunto de fenómenos antes descriptos es lo que caracteriza a lo que genéricamente

llamamos turbulencia, pudiendo resumirse los siguientes contrastes o diferencias entre un flujo
laminar y un flujo turbulento:

FLUJO LAMINAR FLUJO TURBULENTO

El movimiento es ordenado Es estacionario en conjunto en valores medios,
y en conjunto es realmente estacionario. pero contiene inestacionariedades de pequeña escala
Las fuerzas viscosas de origen Las fuerzas viscosas son
molecular son predominantes. relativamente importantes
Las fuerzas de inercia son irrelevantes Las fuerzas de inercia son dominantes.

Los primeros análisis teóricos detallados fueron encarados por O. Reynolds con relación a
su trabajo dos aspectos pueden ser analizados. Por un lado se puede demostrar, consideran-
do las ecuaciones diferenciales de la Mecánica de los Fluidos que las oscilaciones espaciales y
temporales presentes en las componentes de la velocidad afectan al movimiento promedio, y se
manifiestan como una resistencia a la deformación del fluido de un modo similar a la resistencia
que produce la fricción molecular interna a través de la viscosidad en los flujos laminares. La se-
gunda consideración teórica tiene que ver con las condiciones bajo las cuales una perturbación
del movimiento como los que ocurren en los flujos turbulentos puede incrementarse o amor-
tiguarse en el tiempo. Se abordará solamente el análisis del primer aspecto, considerando que el
movimiento es estacionario con fluctuaciones que tienen valores temporales medios. De acuerdo
a la evidencia experimental los movimientos sujetos a fluctuaciones en el espacio y el tiempo
pueden ser tratados considerándolos como integrado por un valor medio constante en el tiempo
y de oscilaciones positivas y negativas alrededor de ese valor medio. Los diferentes autores que
estudiaron el problema difieren en la elección de ese valor medio. Algunos consideran una media
espacial de valores promediados en un espacio definido de moderada extensión, en tanto que
otros toman los valores promediados en un simple punto del espacio durante un tiempo sufi-
cientemente largo, y algunos hacen uso de un concepto combinado. El presente análisis, debido
a Prandtl, considera a las componentes de la velocidad, presión, etc. en un punto del espacio
como compuestas de un valor medio temporal constante tomado como promedio durante un
perı́odo de tiempo suficientemente largo, y de fluctuaciones en el tiempo superpuestas a ese
valor medio. Un movimiento de estas caracterśticas se representa en la figura Nº2.
Si por ejemplo la componente de la velocidad según la dirección ƒx‚se designa con u(t), su
valor medio ū entre el instante to y el tiempo T se obtiene, aplicando el Teorema del Valor
Medio:

1 to +T
Z
ū = u(t)dt (1)
T to
 El tiempo T debe tomarse lo suficientemente grande para que ū sea independiente tanto de
T como de to . Movimientos con estas ceracterı́sticas existen en el mundo real, y se denominan
movimientos turbulentos con valor medio estacionario.

2. Reglas de cálculo para valores medios temporales de magni-

tudes fı́sicas
El tratamiento de ecuaciones con magnitudes fı́sicas utilizando sus medias temporales im-
plica el uso de ciertas reglas de cálculo elementales que determinan la validez de algunas opera-
ciones aritméticas básicas que se presentan en las mismas. Las mas frecuentes son:
0
El valor medio de las fluctuaciones en el tiempo u (t) es nulo:

Siendo
0
u (t) = u(t) − u (2)
0
El valor medio de u (t) es:
Z to +T
0 1 0
u = u (t) dt (3)
T to

llevando la (2) a la (3) y haciendo por comodidad:

Z T Z T Z T
1 1 1
u0 = [u(t) − u] dt = u(t) dt − u dt (4)
T 0 T 0 T 0

Las integrales de la (4) resultan:

Z T
1
u(t) dt = u (5)
T 0

Z T Z T
1 1 T
u dt = u dt = u =u (6)
T 0 T 0 T

de donde se concluye que el valor medio de las fluctuaciones de una magnitud es siempre
nulo:

u0 (t) = 0 (7)

Si A y B son dos magnitudes de un flujo turbulento y “n” es una coordenada espacial

cualquiera se verifican las siguientes reglas para las operaciones elementales con valores
medios temporales :
Suma o diferencia de valores medios:

La suma de sus valores medios es igual al valor medio de la suma:

A+B =A+B (8)

Producto por una constante por un valor medio:

El producto de una constante por el valor medio de una magnitud es igual al valor medio
del producto de la constante por la magnitud:

const × A = const × A (9)

Derivada del valor medio según una dirección:

La derivada según una dirección del valor medio es igual al valor medio de la derivada:

∂A ∂A
= (10)
∂n ∂n

Valores medios de un producto

Si los valores medios de dos magnitudes A y B son nulos el valor medio del producto es
en general distinto de cero:

Z T
1
A= A(t) dt
T 0
Z T
1
B= B(t) dt
T 0
Z T
1
A·B = A(t) · B(t) dt 6= 0 (11)
T 0

Para ilustrar la regla anterior, si se supone que las magnitudes A y B varı́an en el tiempo
según una ley senoidal el gráfico siguiente muestra que el valor medio del producto es
distinto de cero:
3. Aplicación del Teorema de la Cantidad de Movimiento al
flujo turbulento.

Antes de trabajar con las ecuaciones de Navier-Stokes como lo hizo Reynolds para establecer
las ecuaciones que manejan los flujos turbulentos, se puede arribar a conclusiones de un modo
mas directo y obtener algunas conclusiones referidas a la mecánica de este tipo de flujo utilizando
el Teorema de la Cantidad de Movimiento.
Si las componentes de la velocidad según un sistema de referencia cartesiano x, y, z son
u, v, w respectivamente, y se está ante la presencia de un movimiento turbulento estacionario,
ellas se pueden expresar como la suma de un valor medio y sus correspondientes fluctuaciones
en el tiempo del siguiente modo:

u =u + u0
v = v + v0 (12)
w = w + w0

y similarmente para la presión:

p = p + p0 (13)

los valores medios se indican con los términos en barras y son los correspondientes a la parte
del flujo regular y perfectamente ordenado que pueden ser medidos por un observador con un
instrumento estándar como una sonda de presión Pitot-estática, y las fluctuaciones en el tiempo
se representan con las letras con tildes. El caudal en masa que pasa por una superficie elemental
perpendicular al eje “x” en un instante se expresa como:

dq = ρu dA dt (14)

la cantidad de movimiento elemental asociada a ese caudal en dirección “x” es:

dQx = ρu2 dA dt (15)

y las componentes en dirección “y” y “z” son:

dQy = ρuv dA dt (16)

dQz = ρuw dA dt (17)

El primer paso consiste en evaluar los valores medios de las (15), (16) y (17):

Z T
1
dQx = ρu2 dt
T 0
Z T
1
dQy = ρuv dt (18)
T 0
Z T
1
dQz = ρuw dt
T 0
 En las anteriores se hizo t0 = 0 por simplicidad, y se supone que el tiempo T es suficiente-
mente grande. Considerando la componente u y elevándola al cuadrado:

u2 = (u + u0 )2 = u2 + 2u · u0 + u02

y reemplazando en la primera de las (18):

Z T Z T Z T
1 1 1
dQx = ρu2 dt + ρu · u0 dt + ρu02 dt
T 0 T 0 T 0

Aplicando las reglas del punto 2 la segunda integral se anula por ser producto de una
constante u por el valor medio de una fluctuación que es cero, y luego de evaluar las restantes
integrales se tiene para el flujo de cantidad de movimiento medio elemental en dirección “x”:

dQx = (ρu2 + ρu0 2 ) dA (19)

0
Es importante señalar que u0 2 es el valor medio de los cuadrados de u , y no el cuadrado
del valor medio de u0 ,indicado como u0 2 que es nulo de acuerdo a las reglas y definiciones ante-
riormente vistas. La (19) es válida para un flujo incompresible, ya que cuando están presentes
cambios de densidad debe introducirse también a esta como suma de un valor medio y su fluc-
tuación. Aplicando un procedimiento similar a la segunda y tercera integrales de las (18) se
determinan las cantidades de movimientos medios elementales en las otras dos direcciones las
que resultan:

dQy = (ρ u · v + ρ u0 · v 0 ) dA (20)

dQz = (ρ u · w + ρ u0 · w0 ) dA (21)

Dividiendo las ecuaciones anteriores por el área elemental, se obtiene, desde un punto de
vista dimensional una relación fuerza/área, es decir una tensión:

dQx
= −(ρ u2 + ρu0 2 )
dA

dQy
= −(ρ u · v + ρu0 · v 0 ) (22)
dA

dQz
= −(ρ u · w + ρu0 · w0 )
dA

En las ecuaciones (22) se han cambiado los signos de los segundos miembros para tener
en cuenta que la ecuación de la cantidad de movimiento da la fuerza que actúa desde los
alrededores sobre el sistema, y ası́ obtener las tensiones que actúan sobre el elemento fluido.
Puede notarse que las componentes del cambio de la cantidad de movimiento elemental contiene
nuevos términos formados por productos y cuadrados de las fluctuaciones (u0 , v 0 , w0 ) sumados a
términos que provienen de productos de valores medios (u, v, w). Estos últimos son exactamente
los mismos que se obtienen cuando se considera el caso estacionario. El término ρu0 2 representa
tensiones normales, en tanto que ρ u0 v 0 y ρ u0 w0 denotan tensiones tangenciales en las direcciones
“y” y “z” respectivamente. El análisis permite ver claramente que los valores temporales medios
de las fluctuaciones de la velocidad introducen componentes de tensión adicionales que no
están presentes en los flujos laminares. Un estudio mas pormenorizado que conduce a idénticas
conclusiones se hará en el próximo punto.

4. Tensiones aparentes, tensor de tensiones de Reynolds

Las leyes de la viscosidad de Stokes y las ecuaciones de Navier-Stokes son aplicables a un
flujo turbulento si se usan las velocidades y magnitudes reales. Pero se mencionó anteriormente
que desde un punto de vista fı́sico resulta más práctico obtener expresiones a partir de valores
temporales medios, ya que estas son magnitudes que pueden ser medidas por un observador con
instrumentos convencionales como las sondas de presión. A continuación se derivarán las ecua-
ciones aplicables a movimientos turbulentos con valores temporales medios de magnitudes. Por
simplicidad y sin pérdida de generalidad se considerará un movimiento incompresible definido
por la condición ∇.V = 0 a partir de las ecuaciones de Navier–Stokes, las que se escriben:

 
∂u ∂u ∂u ∂u ∂p
ρ +u +v +w = ρKx − + µ∇2 u
∂t ∂x ∂y ∂z ∂x
 
∂v ∂v ∂v ∂v ∂p
ρ +u +v +w = ρKy − + µ∇2 v (23)
∂t ∂x ∂y ∂z ∂y
 
∂w ∂w ∂w ∂w ∂p
ρ +u +v +w = ρKz − + µ∇2 w
∂t ∂x ∂y ∂z ∂z

El primer paso consiste en sustituir todas las variables que aparecen en las ecuaciones (23)
por la suma de las medias temporales y de las componentes debidas a las fluctuaciones. Como
ejemplo, en la primera ecuación haciendo u = u + u0 el término u ∂u
∂x resulta:

(u + u0 ) u ∂u0 0 ∂u 0 ∂u
0
u + u0 ∂x = u + u + u + u (24)
∂(u + u0 ) ∂x ∂x ∂x ∂x

Tomando valores medios temporales y haciendo uso de las reglas antes vistas, el segundo
y tercer término del segundo miembro se anulan porque son productos de constantes (u y ∂u ∂x )
0 ∂u0
por valores medios de fluctuaciones (u y de ∂x ) que se hacen cero. El primer término u.∂u/∂x
permanece invariable al calcular la media, porque la expresión es constante con el tiempo.
Finalmente el último termino es el producto de dos magnitudes variables, y como se vio su
valor medio u0 .∂u0 /∂x no debe anularse necesariamente. Resumiendo al tomar valores medios
temporales del producto analizado se obtiene:

∂(u + u0 ) ∂u ∂u0
u + u0 · =u· + u0 ·
∂x ∂x ∂x

Realizando operaciones similares con todos los términos de las tres ecuaciones de NS se
arriba a las siguientes expresiones válidas para un movimiento turbulento con valores temporales
medios en sus magnitudes:
 !
∂u0 ∂u0 ∂u0
 
∂u ∂u ∂u ∂p
ρ u· +v· +w· = ρKx − + µ∇2 u − ρ u0 · + v0 · + w0 ·
∂x ∂y ∂z ∂x ∂x ∂y ∂z

!
∂v 0 ∂v 0 ∂v 0
 
∂v ∂v ∂v ∂p
ρ u· +v· +w· = ρKy − + µ∇2 v − ρ u0 · + v0 · + w0 · (25)
∂x ∂y ∂z ∂y ∂x ∂y ∂z

!
∂w0 ∂w0 ∂w0
 
∂w ∂w ∂w ∂p
ρ u· +v· +w· = ρKz − + µ∇2 w − ρ u0 · + v0 · + w0 ·
∂x ∂y ∂z ∂z ∂x ∂y ∂z

Estas ecuaciones luego de algunas sustituciones y operaciones algebraicas se pueden escribir

en forma más compacta como:

!
Du ∂p ∂(u0 )2 ∂(u0 · v 0 ) ∂(u0 · w0 )
ρ = ρKx − + µ∇2 u − ρ + +
Dt ∂x ∂x ∂y ∂z
!
Dv ∂p ∂(u0 · v 0 ) ∂(v 0 2 ) ∂(v 0 · w0 )
ρ = ρKy − + µ∇2 v − ρ + + (26)
Dt ∂y ∂x ∂y ∂z
!
Dw ∂p ∂(u0 · w0 ) ∂(v 0 · w0 ) ∂(w0 2 )
ρ = ρKz − + µ∇2 w − ρ + +
Dt ∂z ∂x ∂y ∂z

Una comparación de las ecuaciones anteriores con las correspondientes al flujo laminar,
muestra que ambas son formalmente idénticas, a excepción de aquellos términos que surgen de
productos de las fluctuaciones de las componentes de la velocidad. El análisis dimensional de
esos términos da como unidad fuerza/superficie, y en consecuencia esto sugiere que una forma
apropiada es considerar al flujo turbulento con valores medios temporales como un flujo laminar
con la presencia de fuerzas adicionales que tienen su origen en la turbulencia. Estas fuerzas son
llamadas fuerzas aparentes y tienen las siguientes expresiones:

!
∂(u0 )2 ∂(u0 .v 0 ) ∂(u0 .w0 )
(d fx )ap = −ρ + + .dVol
∂x ∂y ∂z
!
∂(v 0 u0 ) ∂(.v 0 )2 ∂(v 0 .w0
(d fy )ap = −ρ + + .dVol (27)
∂x ∂y ∂z
!
∂(w0 .u0 ) ∂(w0 .v 0 ) ∂(.w0 )2
(d fz )ap = −ρ + + .dVol
∂x ∂y ∂z

Las 9 tensiones asociadas a estas fuerzas se denominan tensiones aparentes o tensiones

de Reynolds, y conforman las componentes de un tensor denominado tensor de tensiones
aparentes de Reynolds,

 0 0
ρ(u u ) ρ(u0 v 0 ) ρ(u0 w0 )
  
τxx τxy τxz
τyx τyy τyz  = ρ(u0 u0 ) ρ(u0 v 0 ) ρ(u0 w0 ) (28)
τzx τzy τzz ap ρ(u0 u0 ) ρ(u0 v 0 ) ρ(u0 w0 )
5. Manifestación de la tensión aparente, viscosidad turbulenta
o de remolino.
Aunque ası́ planteado el concepto de tensión aparente del flujo turbulento pareciera en
cierto modo arbitrario y carente de sustento fı́sico, un análisis comparativo de los mecanismos
que generan las tensiones debidas a la fricción en el flujo laminar permite descubrir algunas
analogı́as entre ambos tipos de tensiones. En los flujos laminares existe un movimiento molecular
aleatorio que macroscópicamente se observa como un movimiento de moléculas perfectamente
ordenado sobre el que se superpone un movimiento o fluctuaciones moleculares al azar. Esto
sumado a las fuerzas de cohesión molecular que son acciones microscópicas se manifiesta o se
hacen sentir como una tensión tangencial entre las capas fluidas. La relación entre la tensión
tangencial y los gradientes del campo de velocidades lleva a introducir el concepto de viscosidad
absoluta o dinámica que es de origen esencialmente molecular.
En un flujo turbulento existe también un flujo medio temporal bien ordenado al que se
superpone una fluctuación al azar de partı́culas o elementos fluidos de carácter macroscópico,
cuyo efecto es similar a la acción molecular de los flujos laminares. Esto provoca también un
efecto macroscópico pero de naturaleza diferente y de carácter mas grosero que se manifiesta
como un valor medio temporal que llamamos tensión aparente. A igual que en el flujo laminar,
la tensión aparente se puede relacionar con el campo de velocidades y obtener una propiedad
media temporal análoga a la viscosidad molecular. Dicha analogı́a conduce a la definición de
un coeficiente de intercambio turbulento de cantidad de movimiento llamado también
viscosidad turbulenta o de remolino
Para determinar la viscosidad turbulenta se puede partir de las ecuaciones de NS, con-
siderando un flujo turbulento estacionario y bidimensional cuya componente “x”se escribe:

!
∂u ∂u ∂p ∂u ∂(u0 .u0 ) ∂(v 0 u0 )
ρu +ρv =− +µ −ρ + (29)
∂x ∂y ∂x ∂y ∂x ∂y

Usualmente se ignora la variación longitudinal de la intensidad turbulenta, ρu0 u0 , y la (29)

se escribe:

 
∂u ∂u ∂p ∂ ∂u ∂p ∂
ρu +ρv =− + µ − ρv 0 u0 =− + (τl + τt ) (30)
∂x ∂y ∂x ∂y ∂y ∂x ∂y

Siendo τl =µ ∂u 0 0
∂y la fricción laminar o molecular τt = −ρv u la fricción turbulenta. Por ana-
logı́a con la ecuación que determina la fricción laminar, se define un coeficiente de intercambio
turbulento de cantidad de movimiento viscosidad turbulenta como:
0 0 0 0
ρv u
τt = −ρv u =εt ∂u
∂y de donde resulta εt = − ∂u/∂y

Debe señalarse que mientras la viscosidad molecular es una propiedad constitutiva que
depende de las propiedades termodinámicas locales del fluido, la viscosidad turbulenta viene
0 0
determinada por las condiciones locales del flujo mismo ( u , v y ∂u/∂y) y la densidad ρ del
fluido. En un flujo turbulento próximo a una pared sólida estacionaria, se pueden identificar
varias regiones gobernadas por diferentes fenómenos. En la zona muy cercana a la pared denom-
inada subcapa laminar, la pared amortigua las fluctuaciones de las partı́culas y la fricción
molecular es predominante. En esta región el perfil de velocidades es lineal lo que implica que:

∂u τ
= =cte (31)
∂y µ
 Entre esta y el flujo irrotacional se desarrolla el flujo turbulento propiamente dicho donde
resulta posible diferenciar tres regiones con distinto comportamiento, denominadas región
interna, región de transición y región externa.
Figura Nº4- Regiones de una capa lı́mite turbulenta.

6. Concepto fı́sico de las tensiones de Reynolds

El mecanismo fı́sico asociado a las tensiones de Reynolds puede ser explicado del siguiente
modo Un flujo de corte turbulento tiene une distribución de velocidades media temporal que
no es uniforme. Debido a la turbulencia de la corriente, pequeñas masas fluidas M1 y N1
intercambian su posición en dirección transversal a la corriente principal, y cada una de ellas
tienda a conservar su cantidad de movimiento. Luego del intercambio M2 tiende a acelerar
al fluido en el que está inmerso y que se mueve mas lentamente, e inversamente N2 trata de
decelerar al fluido en que está inmerso y que se mueve mas rápidamente. El resultado neto de
este intercambio de cantidad de movimiento es una tensión tangencial τxy en la coordenada “y”
actuando en dirección “x”.
Bibliografı́a
1. Fluid Dynamics- MIT Courses- Ascher H. Shapiro.

2. Aerodynanics Theory – W. F. Durand, Vol III, Division G-The mechanicas of Viscous Flow-
L. Prandtl.

3. La mecánica de los Fluidos -Irving. H. Shames