Notas sobre las Leyes de Conservación
C.6 Mecánica de los Fluidos - Ingenierı́a Civil
Dr. Martı́n Romagnoli
4 de noviembre de 2020
Introducción
Los principios gobernantes en la Mecánica de Fluidos son las Leyes de
Conservación para la masa, la cantidad de movimiento (o momentum), y la
energı́a. Veremos cada una de estas leyes en este orden, pudiendo enunciarlas
en sus formas integrales, aplicadas a una región extendida del campo de flujo,
o en sus formas diferenciales, aplicadas en un punto o partı́cula de fluido.
Ambas formas son válidas, y puede una ser derivada de la otra (es decir, a
partir de las formas integrales es posible obtenerse sus formas diferenciales,
y viceversa).
Las formas integrales de las ecuaciones del movimiento de un fluido están
enunciadas en términos de la evolución de un volúmen de control y los flujos
de masa, momentum y energı́a a través de su superficie de control. En este
curso vamos a considerar, y lo veremos luego, volúmenes de control fijos
en el espacio e indeformables. Estas formas integrales son útiles en situa-
ciones donde no se necesite una descripción extremadamente detallada de
las condiciones del flujo (distribuciones de velocidad y presión) y en cambio
estemos interesados en obtener información sobre los valores resultantes, es
decir, de las integrales de dichas distribuciones.
Las formas diferenciales de las ecuaciones del movimiento de un fluido
son ecuaciones diferenciales en derivadas parciales no lineales acopladas para
las variables densidad, velocidad, presión y temperatura. De esta manera,
las formas diferenciales son más apropiadas para un análisis detallado del
flujo. En este curso básico, si bien vamos a comentar y presentar estas formas
diferenciales, en las aplicaciones prácticas utilizaremos las formas integrales.
Conservación de la masa o ecuación de continuidad
Sin considerar reacciones nucleares ni efectos relativistas, la masa no es
creada ni se destruye. El volumen ocupado por una cantidad especı́fica e
identificada de masa es denominado volumen material V (t). Dicho volumen
material se desplaza y se deforma en el interior de un fluido en movimiento
1
de manera que siempre contenga los mismos elementos de masa; ninguna
partı́cula sale del volumen material o ingresa a este. Esto implica que la
superficie material A(t) se debe mover con la velocidad del flujo u. De esta
manera, la conservación de la masa en un volumen material de un fluido
en movimiento puede ser enunciada de la siguiente manera, donde ρ es la
densidad del fluido: ˆ
d
ρ(x, t)dV = 0. (1)
dt
V (t)
Utilizando el Teorema de Transporte de Reynolds (TTR) (ver Notas
sobre Cinemática), con F = ρ y b = u, es posible escribir el término del
lado izquierdo de la igualdad anterior de la siguiente manera:
ˆ ˆ ˆ
d ∂ρ(x, t)
ρ(x, t)dV = dV + ρ(x, t)u(x, t) · ndA = 0. (2)
dt ∂t
V (t) V (t) A(t)
El problema que se nos presenta para la resolución de dicha ecuación es
que el movimiento y la evolución de V (t) y A(t) son determinados por el flu-
jo, el cual es habitualmente desconocido. Por lo tanto, para desarrollar una
ecuación integral que represente la conservación de masa en un volumen ar-
bitrario, V ∗ (t) y A∗ (t), la ecuación 2 debe ser modificada. Esta modificación
está normalmente motivada por la necesidad frecuente de conservar la masa
en un volumen que no es un volumen material, por ejemplo un volumen de
control fijo en el espacio.
Recordando el TTR para un volumen arbitrario V ∗ (t) y siendo F = ρ
obtenemos:
ˆ ˆ ˆ
d ∂ρ(x, t)
ρ(x, t)dV = dV + ρ(x, t)b · ndA. (3)
dt ∂t
V ∗ (t) V ∗ (t) A∗ (t)
Desplazando los términos de la derecha de la igualdad hacia la izquierda en
la ecuación 3:
ˆ ˆ ˆ
d ∂ρ(x, t)
ρ(x, t)dV − dV − ρ(x, t)b · ndA = 0. (4)
dt ∂t
V ∗ (t) V ∗ (t) A∗ (t)
El segundo paso es elegir un volumen de control arbitrario V ∗ (t) que coinci-
da instantáneamente con el volumen material V (t), de manera que V ∗ (t) =
V (t)´ y A∗ (t) = A(t) en el instante
´ t. Debemos tener muy presente que
d d
dt ρ(x, t)dV no es igual a dt ρ(x, t)dV , ya que la primera es la con-
V (t) V ∗ (t)
servación de la masa en un volumen material y la segunda es la derivada
total temporal de una integral sobre un volumen de control arbitrario. Sin
embargo, y por haber hecho coincidir V ∗ (t) = V (t) y A∗ (t) = A(t) en el
instante t, la integral de volumen ∂ρ
∂t en la ecuación 2 es igual a la de la
2
ecuación 4, y la integral de superficie de ρu · n sobre A(t) otorga el mismo
resultado sobre A∗ (t):
ˆ ˆ
∂ρ(x, t) ∂ρ(x, t)
dV = dV , (5)
∂t ∂t
V (t) V ∗ (t)
ˆ ˆ
ρ(x, t)u(x, t) · ndA = ρ(x, t)u(x, t) · ndA. (6)
A(t) A∗ (t)
Si de la ecuación 2 es posible escribir que:
ˆ ˆ
∂ρ(x, t)
dV = − ρ(x, t)u(x, t) · ndA, (7)
∂t
V (t) A(t)
entonces considerando las 2 igualdades anteriores (ecuaciones 5 y 6) es po-
sible observar que:
ˆ ˆ
∂ρ(x, t) ∂ρ(x, t)
dV = dV =
∂t ∂t
V ∗ (t) V (t)
ˆ ˆ
− ρ(x, t)u(x, t) · ndA = − ρ(x, t)u(x, t) · ndA. (8)
A(t) A∗ (t)
Reemplazando esta última expresión en el término central del lado iz-
quierdo de la igualdad en la ecuación 4 nos queda:
ˆ ˆ ˆ
d
ρ(x, t)dV + ρ(x, t)u(x, t) · ndA − ρ(x, t)b · ndA = 0, (9)
dt
V ∗ (t) A∗ (t) A∗ (t)
y reescribiendo:
ˆ ˆ
d
ρ(x, t)dV + ρ(x, t)(u(x, t) − b) · ndA = 0, (10)
dt
V ∗ (t) A∗ (t)
cuya única restricción para u y b es que deben ser observados en el mismo
marco de referencia. La expresión 10 es la ecuación de conservación de la
masa para un volumen de control arbitrario que está en movimiento. De
acuerdo a la elección de b este volumen puede ser estar siendo deformado,
acelerado, en movimiento uniforme o estacionario. En particular, cuando
b = u el volumen de control arbitrario se convierte en un volumen material
y la expresión 10 se reduce a la expresión 1.
Si el volumen de control es estacionario (b = 0), el primer término del
lado izquierdo de la igualdad puede escribirse de acuerdo al TTR según:
ˆ ˆ ˆ
d ∂ρ(x, t)
ρ(x, t)dV = dV + ρ(x, t)0 · ndA. (11)
dt ∂t
V ∗ (t) V ∗ (t) A∗ (t)
3
En la ecuación anterior el segundo término del lado derecho de la igualdad
es nulo, por lo que:
ˆ ˆ
d ∂ρ(x, t)
ρ(x, t)dV = dV , (12)
dt ∂t
V ∗ (t) V ∗ (t)
y reemplazando en la ecuación 10 el primer término del lado izquierdo de la
igualdad con b = 0 obtenemos la siguente:
ˆ ˆ
∂ρ(x, t)
dV + ρ(x, t)u(x, t) · ndA = 0. (13)
∂t
V ∗ (t) A∗ (t)
La ecuación 13 expresa la conservación de la masa en su forma integral
sobre un volumen de control arbitrario fijo en el espacio. Esta es una ecuación
de balance, y expresa que la variación de la masa en el interior del volumen
de control más el flujo neto de masa en su superficie de control debe ser
igual a 0.
Conservación de la cantidad de movimiento o ecua-
ción de momentum
En esta sección, la conservación de cantidad de movimiento equivalente
a la ecuación 13 es obtenida a partir de la 2da Ley de Newton, el principio
fundamental que gobierna el movimiento de los fluidos. Aplicada a un vo-
lumen material V (t) de superfice A(t), la 2da de Ley de Newton puede ser
enunciada de la siguiente manera:
ˆ ˆ ˆ
d
ρ(x, t)u(x, t)dV = ρ(x, t)gdV + f (n, x, t)dA, (14)
dt
V (t) V (t) A(t)
donde ρu es la cantidad de movimiento por unidad de volumen del fluido
en movimiento, g es la fuerza de cuerpo por unidad de masa que actúa sobre
fluido en V (t), f es la fuerza de superficie por unidad de area actuando en
A(t), y n es el versor normal saliente sobre A(t).
Considerando un proceso análogo al realizado para la ecuación de con-
servación de la masa, es posible obtener la ecuación de momentum para un
volumen de control arbitrario estacionario:
4
ˆ ˆ
∂
[ρ(x, t)u(x, t)]dV + [ρ(x, t)u(x, t)]u(x, t) · ndA =
∂t
V ∗ (t) ∗
ˆA (t) ˆ
ρ(x, t)gdV + f (n, x, t)dA. (15)
V ∗ (t) A∗ (t)
El lado derecho de la igualdad en la ecuación 15 no es más que la suma-
toria de fuerzas externas actuantes sobre el volumen de control arbitrario
estacionario. En resumen, la ecuación se lee de la siguiente manera: la va-
riación de momentum en el interior del volumen de control más el flujo neto
de momentum a través de la superficie de control es igual a la resultante de
fuerzas externas actuantes sobre el volumen de control.
Conservación de la energı́a
La primera ley de la termodinámica establece que la energı́a de un sis-
tema es conservada:
δq + δw = ∆e
donde δq es la cantidad de calor por unidad de masa entregada al sistema,
δw es el trabajo por unidad de masa realizado sobre el sistema y ∆e es el
incremento de energı́a interna por unidad de masa del sistema. La energı́a
interna es la manifestación del movimiento molecular aleatorio de los cons-
tituyentes del sistema. En en el caso del flujo de un fluido, debemos sumarle
la energı́a cinética de las partı́culas de fluido del movimiento macroscópico.
De esta manera, cuando aplicada a un volumen material V (t) de superficie
A(t), la conservación de la energı́a interna por unidad de masa e y la energı́a
cinética por unidad de masa |u|2 /2 puede ser escrita de la sgte manera:
ˆ ˆ ˆ ˆ
d 1 2
ρ(e + |u| )dV = ρg · udV + f · udA − q · ndA, (16)
dt 2
V (t) V (t) A(t) A(t)
donde el primer término del lado derecho de la igualdad es el trabajo realiza-
do sobre el fluido en V (t) por fuerzas de cuerpo, el segundo término del lado
derecho el trabajo realizado sobre el fluido en A(t) por fuerzas de superficie,
y el tercer término del lado derecho de la igualdad es el calor extraı́do del
V (t) a través de su A(t).
Expandiendo el término del lado izquierdo de la igualdad en la ecuación
16 a través de la utilización del TTR y considerando un proceso análogo al
realizado para las ecuaciones de conservación de la masa y de momentum, es
posible obtener la ecuación de conservación de la energı́a para un volumen
de control arbitrario estacionario:
5
ˆ ˆ
∂ 1 1
[ρ(e + |u|2 )]dV + [ρ(e + |u|2 )]u · ndA =
∂t 2 2
V ∗ (t) A∗ (t)
ˆ ˆ ˆ
ρg · udV + f · udA − q · ndA. (17)
V ∗ (t) A∗ (t) A∗ (t)
En el caso de un flujo sin fricción, es decir no existen esfuerzos tangencia-
les, el tensor de esfuerzos se reduce a sus fuerzas normales y f = −pdA/dA
(Ver ecuación 39, si µ = 0 entonces τxx = τyy = τzz = −p ). De esta manera
el segundo término del lado derecho de la igualdad puede escribirse:
ˆ ˆ ˆ ˆ
dA p
f · udA = −p · udA = −pu · ndA = − ρu · ndA
dA ρ
A∗ (t) A∗ (t) A∗ (t) A∗ (t)
Al mismo resultado se llega para el caso de un flujo viscoso paralelo en el
que la velocidad del fluido que atraviesa la superficie del volumen de control
es perpendicular en todos sus puntos a dicha superficie de control, si se
desprecian las variaciones hidrostáticas de presión en la sección transversal,
como es el caso del flujo en una cañerı́a, sea este laminar o turbulento (ya
lo veremos más adelante).
Reemplazando ésta última en la ecuación 17 y expresándola en el lado
izquierdo de la igualdad es posible obtener:
ˆ ˆ
∂ 1 1
[ρ(e + |u|2 )]dV + [ρ(e + |u|2 )]u · ndA+
∂t 2 2
V ∗ (t) A∗ (t)
ˆ ˆ ˆ
p
ρu · ndA = ρg · udV − q · ndA. (18)
ρ
A∗ (t) V ∗ (t) A∗ (t)
Reagrupando
ˆ ˆ
∂ 1 1 p
[ρ(e + |u|2 )]dV + [ρ(e + |u|2 + )]u · ndA+ =
∂t 2 2 ρ
V ∗ (t) A∗ (t)
ˆ ˆ
ρg · udV − q · ndA. (19)
V ∗ (t) A∗ (t)
Considerando que solamente existen fuerzas de cuerpo ocasionadas por el
campo gravitatorio, el primer término del lado derecho de la igualdad puede
reescribirse con ayuda del teorema de Gauss en una integral de superficie,
6
ˆ ˆ ˆ
ρg · udV = − ∇ · ρgzudV = − ρgzu · ndA
V ∗ (t) V ∗ (t) A∗ (t)
Con lo cual la ecuación 19 puede expresarse de la sgte forma:
ˆ ˆ
∂ 1 1 p
[ρ(e + |u|2 )]dV + [ρ(e + |u|2 + )]u · ndA
∂t 2 2 ρ
V ∗ (t) A∗ (t)
ˆ ˆ
+ ρ(gz)u · ndA = − q · ndA. (20)
A∗ (t) A∗ (t)
La ecuación 20 expresa la conservación de la energı́a en un volumen de
control arbitrario fijo y cuya única fuerza de cuerpo es la generada por el
campo gravitatorio, tanto para el caso de un flujo invı́scido (es decir, donde
los esfuerzos tangenciales son nulos) como para el caso de un flujo viscoso
paralelo (como lo es el escurrimiento en cañerı́as, sea laminar o turbulento).
Además, no existen adiciones ni sustracciones de energı́a al flujo mediante
bombas o turbinas, respectivamente.
Formas diferenciales de las leyes básicas de conser-
vación de masa y momentum para un flujo incom-
presible. Ecuaciones de Navier Stokes
Hasta aquı́ hemos desarrollado las leyes básicas en sus formas integrales
y, como hemos mencionado, serán muy útiles para la resolución de nuestros
problemas en la práctica. Sin embargo, la información que nos proveen éstas
ecuaciones no presenta un nivel de detalle adecuado si se pretende conocer
distribución del campo de velocidades y fuerzas (normales y tangenciales).
A tal fin es necesario hechar mano de las leyes básicas en sus formas dife-
renciales.
Partimos entonces de la ecuación de conservación de la masa en su forma
integral: ˆ ˆ
∂ρ(x, t)
dV + ρ(x, t)u(x, t) · ndA = 0, (21)
∂t
V ∗ (t) A∗ (t)
y transformando el segundo término del lado derecho en una integral de
volumen por intermedio del Teorema de Gauss:
ˆ ˆ
ρ(x, t)u(x, t) · ndA = ∇ · ρ(x, t)u(x, t)dV ,
A∗ (t) V ∗ (t)
7
es posible obtener:
ˆ ˆ
∂ρ(x, t)
dV + ∇ · ρ(x, t)u(x, t)dV = 0, (22)
∂t
V ∗ (t) V ∗ (t)
Reagrupando y empleando el Teorema del Valor Medio obtenemos la forma
general diferencial de la ecuación de continuidad:
∂ρ(x, t)
+ ∇ · ρ(x, t)u(x, t) = 0, (23)
∂t
la cual también podemos expresarla en notación ı́ndice1 :
∂ρ ∂
+ (ρui ) = 0. (24)
∂t ∂xi
Expandiendo el segundo término del lado izquierdo de la igualdad:
∂ρ ∂ ∂ ∂
+ (ρux ) + (ρuy ) + (ρuz ) = 0. (25)
∂t ∂x ∂y ∂z
Considerando ahora la derivada de un producto es posible expresar la ecua-
ción anterior de la siguiente manera:
∂ux ∂ρ ∂uy ∂ρ ∂uz ∂ρ ∂ρ
ρ + ux + ρ + uy + ρ + uz + = 0. (26)
∂x ∂x ∂y ∂y ∂z ∂z ∂t
La cual puede ser escrita de la siguiente forma:
∂ux ∂uy ∂uz ∂ρ ∂ρ ∂ρ ∂ρ
ρ( + + ) + ux + uy + uz + = 0. (27)
∂x ∂y ∂z ∂x ∂y ∂z ∂t
∂ux ∂uy ∂uz Dρ
+ρ( + )+ = 0. (28)
∂x ∂y ∂z Dt
Si consideramos que el fluido es incompresible, la derivada material es nula
(o sea, Dρ ∂ρ ∂ρ ∂ρ ∂ρ
Dt = ux ∂x + uy ∂y + uz ∂z + ∂t = 0) y la ecuación anterior se reduce
a la siguiente expresión:
∂ux ∂uy ∂uz
ρ( + + ) = 0. (29)
∂x ∂y ∂z
ux Axx Axy Axz
1
u = ui = uy tiene un único ı́ndice libre i. Aij = Ayx Ayy Ayz tiene dos
uz Azx Azy Azz
(Axx nx + Axy ny + Axz nz )
ı́ndices libres i y j. Aij nj = Aix nx + Aiy ny + Aiz nz = (Ayx nx + Ayy ny + Ayz nz ) tiene
(Azx nx + Azy ny + Azz nz )
un ı́ndice libre i y un ı́ndice auxiliar j. Si no queda claro ver capı́tulo 7.9, pág. 262 del
libro de Shames.
8
La expresión entre paréntesis no es más que la divergencia2 de la velocidad,
la cual puede ser expresada en notación vectorial como ∇ · u o en notación
ı́ndice según ∂u
∂xi .
i
De esta manera, como en la ecuación 29 la densidad no puede ser nula, la
ecuación de continuidad en su forma diferencial para el flujo de un fluido
incompresible es:
∂ui ∂ux ∂uy ∂uz
∇·u= = + + = 0. (30)
∂xi ∂x ∂y ∂z
Es decir, la divergencia de la velocidad debe ser nula si el flujo es incompre-
sible.
Del mismo modo partimos ahora de la ecuación de conservación de la
cantidad de movimiento en su forma integral:
ˆ ˆ
∂
[ρ(x, t)u(x, t)]dV + [ρ(x, t)u(x, t)]u(x, t) · ndA =
∂t
V ∗ (t) A∗ (t)
ˆ ˆ
ρ(x, t)gdV + f (n, x, t)dA. (31)
V ∗ (t) A∗ (t)
Transformando las integrales de superficie en integrales de volumen a través
del Teorema de Gauss,
ˆ ˆ
[ρ(x, t)u(x, t)]u(x, t) · ndA = ∇ · [ρ(x, t)u(x, t)u(x, t)]dV
A∗ (t) V ∗ (t)
ˆ
∂
= (ρui uj )dV
∂xi
V ∗ (t)
y ˆ ˆ ˆ
∂
f (n, x, t)dA = ni Tij dA = (Tij )dV ,
∂xi
A∗ (t) A∗ (t) V ∗ (t)
podemos escribirla de la siguiente manera:
ˆ
∂ ∂ ∂
{ (ρuj ) + (ρui uj ) − ρgj − (Tij )}dV = 0. (32)
∂t ∂xi ∂xi
V ∗ (t)
2 1
‚
La definición matemática del operador divergencia es ∇ · F = lı́m∆V →0 ∆V A
F · dA,
donde ∆V es cualquier volumen elemental en el espacio y A la superficie de este volumen
elemental.
9
Aplicando el Teorema del Valor Medio podemos deshacernos del opera-
dor integral y reescribirla de la siguiente manera:
∂ ∂ ∂
(ρuj ) + (ρui uj ) = ρgj + (Tij ). (33)
∂t ∂xi ∂xi
Trabajando sobre términos del lado izquierdo de la igualdad, los cuales
incluyen derivadas de productos, pueden ser expresados según:
∂ ∂ ∂uj ∂ρ ∂ ∂uj
(ρuj ) + (ρui uj ) = ρ + uj + uj (ρui ) + ρui =
∂t ∂xi ∂t ∂t ∂xi ∂ui
∂uj ∂ρ ∂ ∂uj
ρ + uj [ + (ρui )] + ρui .
∂t ∂t ∂xi ∂ui
La expresión entre corchetes [ ∂ρ ∂
∂t + ∂xi (ρui )] no es otra que la ecuación de
continuidad, por lo que reemplazando la ecuación 23 o 24 en la anterior y, a
su vez, ésta en la ecuación 33, es posible obtener la forma final:
∂uj ∂uj ∂
ρ + ρui = ρgj + (Tij ). (34)
∂t ∂ui ∂xi
D
Utilizando la definición de derivada total Dt , es posible escribirla:
Duj ∂
ρ = ρgj + (Tij ). (35)
Dt ∂xi
El ı́ndice j es libre y el ı́ndice i es auxiliar, en este caso j denota el carácter
vectorial de la ecuación. Las tensiones en un punto pueden ser representadas
completamente por las nueve componentes del tensor de tensiones Tij . El
primer ı́ndice indica la dirección de la normal a la superfice sobre la cual
es considerada la tensión, mientras que el segundo ı́ndice indica la dirección
sobre la cual actúa la tensión. Los elementos de la diagonal del tensor son
las tensiones normales, mientras que los elementos fuera de la diagonal son
las tensiones tangenciales. El tensor de esfuerzos es simétrico, esto es, es
válida la ley de complementaridad de esfuerzos, Tij = Tji . Expandiendo
el término del lado izquierdo y el segundo término del lado derecho de la
igualdad anterior obtenemos:
∂uj ∂uj ∂uj ∂uj ∂τxj ∂τyj ∂τzj
ρ( + ux + uy + uz ) = ρgj + + + . (36)
∂t ∂x ∂y ∂z ∂x ∂y ∂z
Se podrá observar aquı́ más claramente que el ı́ndice libre j en la ecuación
manifiesta el carácter vectorial de la ecuación de momentum. Por lo tanto,
cada una de las componentes de la ecuación anterior pueden ser escritas de
10
la siguiente forma:
∂ux ∂ux ∂ux ∂ux ∂τxx ∂τxy ∂τxz
ρ( + ux + uy + uz ) = ρgx + + +
∂t ∂x ∂y ∂z ∂x ∂y ∂z
∂uy ∂uy ∂uy ∂uy ∂τyx ∂τyy ∂τyz
ρ( + ux + uy + uz ) = ρgy + + + (37)
∂t ∂x ∂y ∂z ∂x ∂y ∂z
∂uz ∂uz ∂uz ∂uz ∂τzx ∂τzy ∂τzz
ρ( + ux + uy + uz ) = ρgz + + +
∂t ∂x ∂y ∂z ∂x ∂y ∂z
Se ha indicado anteriormente que el esfuerzo cortante τ en el caso de un
flujo paralelo esta dado por:
∂u
τ =µ (38)
∂n
donde n es la dirección coordenada perpendicular a la interfaz donde se
desarrolla el esfuerzo tangengial y µ es el coeficiente de viscosidad. Ésta es
la muy conocida ley de viscosidad de Newton. La relación entre tensiones y
deformaciones en un continuo es denominada ecuación constitutiva. La ecua-
ción constitutiva general para el caso de un flujo incompresible e isotrópico
puede escribirse de la siguiente forma:
∂ux ∂uy ∂ux
τxx = 2µ −p τxy = µ( + )
∂x ∂x ∂y
∂uy ∂uz ∂ux
τyy = 2µ −p τxz = µ( + ) (39)
∂y ∂x ∂z
∂uz ∂uy ∂uz
τzz = 2µ −p τzy = µ( + ),
∂z ∂z ∂y
donde p = − 31 (τxx + τyy + τzz ) es la presión termodinámica que hemos visto
en el capı́tulo sobre hidrostática. Los flujos que cumplen con esta última
ecuación se denominan Newtonianos (para el caso de un flujo paralelo, donde
uy = uz = 0 y la única componente de velocidad distinta de cero es ux , el
esfuerzo cortante en la ecuación 39 se reduce a τxy = µ ∂u ∂y , de esta manera
x
vemos que τxy se reduce a la ecuación 38).
Reemplazando los esfuerzos normales y cortantes de la ecuación 37 por la
ecuación constitutiva general (ecuación 39), después de agrupar y ordenar
términos (ayuda: cuando agrupen y ordenen términos debemos tener en
cuenta que estamos trabajando con un flujo incompresible, es decir ∇ · u =
∂ux ∂uy ∂uz
∂x + ∂y + ∂z = 0), obtenemos las ecuaciones de Navier-Stokes para un
flujo incompresible:
11
∂ux ∂ux ∂ux ∂ux ∂p
ρ( + ux + uy + uz ) = ρgx + (− + µ∇2 ux )
∂t ∂x ∂y ∂z ∂x
∂uy ∂uy ∂uy ∂uy ∂p
ρ( + ux + uy + uz ) = ρgy + (− + µ∇2 uy ) (40)
∂t ∂x ∂y ∂z ∂y
∂uz ∂uz ∂uz ∂uz ∂p
ρ( + ux + uy + uz ) = ρgz + (− + µ∇2 uz )
∂t ∂x ∂y ∂z ∂z
La misma ecuación (40) puede ser expresada en notación vectorial como
sigue:
Du
ρ = ρg + (−∇p + µ∇2 u) (41)
Dt
Las ecuaciones 30 y 40 forman un sistema de cuatro ecuaciones diferen-
ciales simultáneas con cuatro incógnitas, ux , uy , uz y p, las cuales podrı́a
pensarse en resolver en muchos problemas, si no fuera por su naturaleza no
lineal (la no linealidad está dada en la derivada sustancial o material don-
de se encuentran términos convectivos, ux ∂ux /∂x, etc. ) y la complejidad
general de las ecuaciones.
En el caso de un flujo invı́scido, donde se desprecian las fuerzas viscosas
(es decir, no existen tensiones tangenciales), el 3er término del lado derecho
de la igualdad es nulo y la ecuación 41 se reduce a la ecuación de Euler:
Du
ρ = ρg − ∇p (42)
Dt
Ecuación de Bernoulli
En las secciones precedentes se han presentado varias formas de las le-
yes de conservación de la materia, momentum y energı́a. La ecuación de
Bernoulli no es una nueva ley de conservación, sino que se deriva de la ecua-
ción de la conservación de la energı́a o de la ecuación de la conservación de
momentum.
La ecuación 20 describe la conservación de la energı́a tanto para el caso
de un flujo invı́scido (es decir, donde los esfuerzos tangenciales son nulos)
como para el caso de un flujo viscoso paralelo en cañerı́as (u = (ux , 0, 0)).
Considerando ahora que el flujo es permanente (∂/∂t = 0) y no hay trans-
ferencia de calor (q=0), la ecuación 20 puede ser expresada de la siguiente
manera:
ˆ
1 p
[ρ(e + |u|2 + + gz)]u · ndA = 0. (43)
2 ρ
A∗ (t)
12
Si además consideramos que el flujo es isentrópico (es decir, no hay va-
riación de la entropı́a, lo cual se cumple para el caso de un flujo invı́scido sin
transferencia de calor) e incompresible (Dρ/Dt = 0), la variación de energı́a
interna es nula (De/Dt = 0). Siendo que el volumen de control posee una
entrada (ent) y una salida (sal), luego de aplicar la ecuación de continui-
dad (el´ flujo neto de masa a través de la superficie de control es nulo, es
decir ρu · ndA = −ρQent + ρQsal = 0 donde Q es el caudal) y despreciar
la variación hidrostática de presión en las secciones de entrada y salida, la
ecuación anterior se reduce a:
1 p 1 p
( |u|2 + + gz)ent = ( |u|2 + + gz)sal . (44)
2 ρ 2 ρ
La ecuación 44 es normalmente conocida como la ecuación de Bernoulli
para el caso de un flujo permanente, no viscoso, incompresible y sin inter-
cambio de calor.
A la misma expresión es posible llegar partiendo de la ecuación de Euler
(42), la cual es válida para flujos incompresibles e invı́scidos (sin efectos
viscosos). Para eso primero vamos a reescribir la aceleración convectiva en
términos de el producto vectorial entre la velocidad y la vorticidad (ω), y el
gradiente de la energı́a cinética por unidad de masa:
(u · ∇)u = −u × ω + ∇(|u|2 /2) (45)
De esta manera, reemplazando la ecuación 45 en la ecuación 42 obtene-
mos:
Du ∂u ∂u
ρ =ρ + ρ(u · ∇)u = ρ − ρu × ω + ρ∇(|u|2 /2) = ρg − ∇p (46)
Dt ∂t ∂t
En el caso de un flujo incompresible, los efectos viscosos desaparecen
cuando el flujo es irrotacional (ω=0) o tiene una vorticidad uniforme (ω =
cte). Si consideramos que el flujo es permanente e irrotacional (∂/∂t = 0 y
ω=0), la ecuación 46 se reduce a:
ρ∇(|u|2 /2) = ρg − ∇p (47)
Si escribimos la aceleración de la gravedad como g = −∇(gz) y reem-
plazamos en la ecuación anterior es posible obtener:
p
∇(|u|2 /2 + gz + ) = 0 (48)
ρ
lo cual es equivalente a:
p
|u|2 /2 + gz + = constante. (49)
ρ
Notar que la ecuación 49 es idéntica a la ecuación 44.
13
Flujo entre placas paralelas (flujo de Poiseuillie)
Consideremos la situación que se describe en el esquema de la Figu-
ra 1, en la cual un flujo incompresible fluye entre placas paralelas al eje x
(asumimos a y como eje vertical). En esta situación, solamente tendremos
velocidad en la componente de la dirección del escurrimiento, por lo que
u = (ux (x, y), 0, 0) (uy = uz = 0). Además, el flujo debe ser independiente
∂
de la dirección z, por lo que ∂z =0. Como consecuencia, la ecuación 40 puede
ser expresada de la siguiente manera (considerando que la única fuerza de
cuerpo es la producida por el campo gravitatorio: gy = −g y gx = gz = 0):
∂ux ∂ux ∂p
ρ( + ux + 0 + 0) = (− + µ∇2 ux )
∂t ∂x ∂x
∂p (50)
ρ(0 + 0 + 0 + 0) = ρgy + (− + 0)
∂y
ρ(0 + 0 + 0 + 0) = 0 + (0 + 0).
Obserando las ecuaciones anteriores es posible señalar que para un flujo
paralelo e incompresible la distribución de presiones sobre el eje vertical y
es hidrostática.
∂
Si además de paralelo e incompresible, el flujo es permanente ( ∂t = 0) y
2
no viscoso (µ∇ ux = 0), la ecuación 50 se simplifica de la siguiente manera:
∂ux ∂p
ρux =−
∂x ∂x
(51)
∂p
−ρg = .
∂y
De esta manera, la ecuación 51 es válida entonces para un flujo paralelo,
incompresible, permamente y no viscoso.
Otra simplificación útil de la ecuación 50 es asumir que el flujo es per-
∂
mamente ( ∂t = 0), está desarrollado (es decir, el perfil de velocidades no
cambia en el sentido del escurrimiento, u = (ux (y), 0, 0), la componente x
de la velocidad ahora es solamente función de y; ∂u ∂x = 0), y se considera-
x
2
ran los efectos de la viscosidad (µ∇ ux 6= 0). A partir de lo mencionado,
la ecuación que describe el movimiento de un flujo paralelo, incompresible,
permanente, desarrollado y viscoso es la siguiente:
∂p
µ∇2 ux =
∂x
(52)
∂p
−ρg = .
∂y
14
Figura 1: Flujo entre placas paralelas.
Bibliografia
Kundu P., Cohen I.M. and Dowling D.R. (2016) Fluid Mechanics. 6th
Edition. 922 pag. Academic Press.
15