UNIVERSIDAD DE CHILE

DEPARTAMENTO DE INGENIERIA CIVIL

CI 3101, Mecánica de Fluidos
Prof. Aldo Tamburrino Tavantzis

ECUACIONES DE REYNOLDS
ALDO TAMBURRINO TAVANTZIS

Notación Indicial

Para ahorrar escritura, se usará la notación indicial:

r
x i = x = ( x1 , x 2 , x 3 ) = ( x , y , z )
r
u i = u = ( u1 , u 2 , u 3 ) = ( u ,v , w )

La repetición de índices repetidos indica suma sobre ese índice. De este modo, la
ecuación de continuidad para un fluido incompresible se escribe como:

∂u i
=0
∂x i

Al sumar sobre el índice i = 1, 2, 3 se tiene:

∂u1 ∂u 2 ∂u 3 ∂u ∂v ∂w
+ + = + + 0
∂x 1 ∂x 2 ∂x 3 ∂x ∂y ∂z

La ecuación del movimiento de un fluido newtoniano, incompresible en un campo

gravitacional está dada por la ecuación de Navier-Stokes:

∂u i ∂u i 1 ∂p̂ ∂ 2u i
+ uj =− +ν
∂t ∂x j ρ ∂x i ∂x 2j

En la deducción de la ecuación anterior no hubo restricción respecto al régimen de

flujo, por lo que puede ser aplicada tanto a flujos laminares como a turbulentos.

Descomposición de Reynolds

Si en un punto del flujo se mide cualquier variable e de un flujo turbulento en

función del tiempo, el registro que se obtiene es como el que se muestra a
continuación:

1
 UNIVERSIDAD DE CHILE
DEPARTAMENTO DE INGENIERIA CIVIL
CI 3101, Mecánica de Fluidos
Prof. Aldo Tamburrino Tavantzis

e
e’

Osborne Reynolds propuso descomponer la variable e en dos componentes: una

media temporal, e, y una fluctuación en torno a la media, e’. De este modo, se
puede escribir:

e( x i , t ) = e ( x i ) + e' ( x i , t )

donde el promedio temporal se define como:

1 T
∫
lim
e ( xi ) = e( x i , t )dt
T → ∞T 0

De la definición, es fácil ver que se cumple:

e' = 0
e =e
∂e ∂e
=
∂x i ∂x i

∫ edx = ∫ e dx
i i

Consideremos un flujo turbulento en el que sus propiedades estadísticas no varían

en el tiempo, en particular su valor medio temporal, su desviación estándar y otros
momentos de orden superior.

La descomposición de Reynolds aplicada a la ecuación de continuidad se traduce

en:

2
 UNIVERSIDAD DE CHILE
DEPARTAMENTO DE INGENIERIA CIVIL
CI 3101, Mecánica de Fluidos
Prof. Aldo Tamburrino Tavantzis

∂ ( u i + u i ′ ) ∂u i ∂u i ′
= + =0
∂x i ∂x i ∂x i

Promediando en el tiempo la ecuación anterior:

∂u i ∂u i ′
+ =0
∂x i ∂x i

De donde resulta:

∂u i ∂u i ′
=0 =0
∂x i ∂x i

O sea, la ecuación de continuidad se satisface tanto para las componentes medias

como para las fluctuaciones de la velocidad.

Las ecuación de continuidad instantánea podemos escribirla indistintamente con

el subíndice i o j, ya que ambos cubren (1, 2, 3), o sea:

∂u j
=0
∂x j

Multiplicando la ecuación anterior (que es un escalar) por el vector ui :

∂u j
ui =0
∂x j

Sumando esta ecuación a la de momentum, resulta:

∂u i ∂u i u j 1 ∂p̂ ∂ 2u i
+ =− +ν
∂t ∂x j ρ ∂x i ∂x 2j

Aplicando la descomposición de Reynolds:

′ ′
∂ui + ui ′ ∂ ( ui + u i )( u j + u j ) 1 ∂ ( p̂ + p̂´) ∂ 2 ( ui + ui ′ )
+ = − + ν
∂t ∂x j ρ ∂x i ∂x 2j

Tomando promedio temporal de esta ecuación:

3
 UNIVERSIDAD DE CHILE
DEPARTAMENTO DE INGENIERIA CIVIL
CI 3101, Mecánica de Fluidos
Prof. Aldo Tamburrino Tavantzis

′ ′
∂ui + ui ′ ∂ ( ui + ui )( u j + u j ) 1 ∂ ( p̂ + p̂' ) ∂ 2 ( ui + ui ′ )
+ =− +ν
∂t ∂x j ρ ∂x i ∂x 2j

O sea,

′ ∂ u ′u ∂ u i ′u j ′  ∂ 2u
∂u i ∂u i ′ ∂u i u j ∂u i u j 1  ∂p̂ ∂p̂ ′   ∂ 2ui ′ 
+ + + +
i j
+ =−  +  +ν i
+
∂t ∂t ∂x j ∂x j ∂x j ∂x j ρ  ∂x i ∂x i   ∂x 2j ∂x 2j 
 

Varios de los términos de la ecuación anterior son nulos. Por ejemplo:

∂ui u' j ∂ ui u' j ∂ui u' j

= = =0
∂x j ∂x j ∂x j

Considerando, además, que las propiedades estadísticas del flujo no varían en el

tiempo, la ecuación de momentum se reduce a:

∂u i u j ∂ u i ′u j ′ 1 ∂p̂ ∂ 2u i
+ =− +ν
∂x j ∂x j ρ ∂x i ∂x 2j

Expandiendo la derivada del producto, el primer término del lado izquierdo queda:

∂u i u j ∂u j ∂u i
= ui + uj
∂x j ∂x j ∂x j

Como se está considerando un fluido incompresible, la derivada del producto se

reduce a:

∂u i u j ∂u i
= uj
∂x j ∂x j

por lo que la ecuación de momentum se escribe como:

′ ′
∂u i ∂ u i u j 1 ∂p̂ ∂ 2u i
uj + =− +ν
∂x j ∂x j ρ ∂x i ∂x 2j

4
 UNIVERSIDAD DE CHILE
DEPARTAMENTO DE INGENIERIA CIVIL
CI 3101, Mecánica de Fluidos
Prof. Aldo Tamburrino Tavantzis

Con el objeto de dar a la expresión anterior la forma de la segunda ley de Newton,

por unidad de masa, pero aplicada a cantidades medias temporales, o sea:

a =∑f

el segundo término del lado izquierdo lo pasamos al lado derecho:

′ ′
∂u i 1 ∂p̂ ∂ 2u i ∂ u i u j
uj =− +ν −
∂x j ρ ∂x i ∂x 2j ∂x j

El efecto que tiene este simple paso algebraico no es menor, ya que podemos
∂ u i ′u j ′
interpretar de otra manera al término . Este término no es más que un
∂x j
flujo de momentum de las fluctuaciones turbulentas, pero al pasarlo al lado
derecho, podemos interpretarlo como una fuerza, resultante de la turbulencia del
flujo. Insisto, no es una fuerza turbulenta, lo interpretamos así ya que simplifica el
entendimiento de las ecuaciones y fenómenos turbulentos.

Siguiendo con el álgebra:

∂u i u j 1 ∂p ∂ 2 u j ∂ u i ′u j ′
=− +ν −
∂x j ρ ∂x i ∂x 2j ∂x j

∂u i u j 1 ∂p̂ ∂  ∂u i 
=− + ν − u i ′u j ′ 
∂x j ρ ∂x i ∂x j  ∂x j 
 

Multiplicando la ecuación de continuidad por la viscosidad:

∂u j
ν =0
∂x j

y derivando respecto a xi :

∂ ∂u j ∂ ∂u j
ν =ν =0
∂x i ∂x j ∂x j ∂x i

sumando a la ecuación de momentum:

5
 UNIVERSIDAD DE CHILE
DEPARTAMENTO DE INGENIERIA CIVIL
CI 3101, Mecánica de Fluidos
Prof. Aldo Tamburrino Tavantzis

∂u i u j 1 ∂p̂ 1 ∂  ∂u i  µ ∂ ∂u j
=− + µ − ρ u i ′u j ′  +
∂x j ρ ∂x i ρ ∂x j  ∂x j  ρ ∂x j ∂x i
 

agrupando términos:

∂u i u j 1 ∂p̂ 1 ∂   ∂u ∂u j  
=− +  µ i +  − ρ u i ′u j ′ 
∂x j ρ ∂x i ρ ∂x j   ∂x j ∂x i  
   

el término en paréntesis que multiplica la viscosidad dinámica es la tasa de

deformación, εij, y µεij corresponde a los esfuerzos de origen viscoso, los que
denotaremos como τvij. El subíndice v es para enfatizar el origen viscoso de estos
esfuerzos.

De este modo, la ecuación de momentum queda:

∂u i u j 1 ∂p̂ 1 ∂  ′ ′
=− + τ vij − ρ u i u j 
∂x j ρ ∂x i ρ ∂x j  

El término τ vij − ρ ui ′u j ′ lo interpretamos como el esfuerzo total actuando en un

flujo turbulento, cuando éste se representa en términos de sus cantidades medias
temporales. Tiene dos componentes, una de origen viscoso, τvij , y otra de origen
turbulento, − ρ u ′u ′ , por lo que podemos escribir un esfuerzo total τ :
i j ij

τ ij = τ vij + τ tij

donde τ tij = − ρ u i ′u j ′ corresponde al efecto turbulento y se denomina esfuerzos

aparentes de Reynolds o, simplemente, esfuerzos de Reynolds o esfuerzos
turbulentos.

Luego, las ecuaciones que rigen las propiedades medias temporales de un flujo
turbulento den un fluido newtoniano incompresible son:

Continuidad:

∂u j
=0
∂x j
Cantidad de movimiento:

6
 UNIVERSIDAD DE CHILE
DEPARTAMENTO DE INGENIERIA CIVIL
CI 3101, Mecánica de Fluidos
Prof. Aldo Tamburrino Tavantzis

∂u i u j 1 ∂p̂ 1 ∂  ′ ′
=− + τ vij − ρ u i u j 
∂x j ρ ∂x i ρ ∂x j  

El sistema de ecuaciones anteriores se conoce como ecuaciones de Reynolds.

Siendo más preciso en la definición del promedio temporal, o sea considerando un

T lo suficientemente grande para eliminar las fluctuaciones turbulentas, pero no
tan grande como para eliminar la variación de la dependencia que puedan tener
las características medias, es posible incluir la dependencia de las cantidades
medias con el tiempo, resultando:

∂u i ∂u i u j 1 ∂p̂ 1 ∂  ′ ′
+ =− + τ vij − ρ ui u j 
∂t ∂x j ρ ∂x i ρ ∂x j  

Las ecuaciones de Reynolds presentan un gran avance para el estudio de los

fenómenos turbulentos, pero no resuelve el problema, ya que τtij sigue siendo
desconocido. Considerando la simetría de los esfuerzos turbulentos, τtij = τtji,
resulta que tenemos 4 ecuaciones (una de continuidad y tres de momentum) y 10
incógnitas: u ,v , w , p̂ , u' 2 ,v' 2 , w' 2 , u' v' , u' w' ,v' w' . Esto es lo que se conoce como el
problema del cierre de la turbulencia y consiste en el hecho de tener más incógnitas
que ecuaciones. Luego, el próximo problema a resolver es determinar los esfuerzos
turbulentos τtij, pero esa es materia de la próxima clase.