Análisis Diferencial en Mecánica de Fluidos

2.
Análisis diferencial en Mecánica de Fluidos 1
UNIVERSIDAD DE OVIEDO
Escuela Politécnica Superior de Ingeniería de Gijón
Ingenieros Industriales
Curso 2004-2005
Apuntes de Mecánica de Fluidos
ANÁLISIS DIFERENCIAL EN MECÁNICA DE FLUIDOS.
1. Cinemática y dinámica de fluidos.
2. Ecuaciones de constitución.
3. Ecuaciones de conservación
4. Problemas resueltos.
Julián Martínez de la Calle

Área de Mecánica de Fluidos
Gijón noviembre 2004
_________________________________________________________________________________________________________________
Apuntes de Mecánica de Fluidos JMC 04
2. Análisis diferencial en Mecánica de Fluidos 2
ANÁLISIS DIFERENCIAL EN MECÁNICA DE FLUIDOS
1. CINEMÁTICA Y DINÁMICA DE FLUIDOS
1.1. Métodos de análisis en Mecánica de Fluidos: Lagrangiano y Euleriano. 3

1.2. Tipos de análisis en Mecánica de Fluidos: Diferencial, Dimensional e Integral. 4
1.3. Cinemática de Fluidos: propiedades del vector velocidad. 5
1.3.1. Derivada sustancial o material. 5
1.3.2. El operador NABLA: gradiente, divergencia y laplaciana. 7
1.4. Dinámica de Fluidos: fuerzas macroscópicas sobre los fluidos. 10
1.4.1. Fuerzas de volumen. 10
1.4.2. Fuerzas de superficie. 11
1.4.3. Tensor de tensiones. 11
1.5. Tipos de flujo. 14
2. ECUACIONES DE CONSTITUCIÓN
2.1. Comportamiento mecánico: tensor de velocidad de deformación. 15

2.2. Fluidos Stokesianos: tensor de tensiones viscosas. 18
2.3. Fluidos Newtonianos. 18
2.3.1. Ec. de NAVIER-POISSON. 18
2.3.3. Tensor de tensiones de un fluido newtoniano. 19
2.4. Comportamiento térmico. 20
2.4.1. Ecuaciones de Estado.
2.4.2. Ecuaciones de transmisión de calor.
3. ECUACIONES DE CONSERVACIÓN.
3.1. Ecuación diferencial de conservación de masa: ecuación de continuidad. 21

3.2. Ecuación diferencial de conservación de cantidad de movimiento: ecuación de movimiento de CAUCHY. 24
3.2.1. Fluido no viscoso: ecuación de EULER. 26
3.2.2. Flujo no viscoso en línea de corriente: ecuación de BERNOULLI. 28
3.2.3. Fluido newtoniano: ecuaciones de NAVIER-STOKES. 29
3.2.4. La función de corriente y la función potencial de velocidad. 30
3.3. Ecuación diferencial de conservación de la energía: ecuación de la energía. 32
3.3.1. Ecuación de energía interna. 33
3.3.2. Ecuación de entalpía. 34
3.3.3. Ecuación de entropía. 34
3.4. Condiciones de contorno. 35
4. PROBLEMAS RESUELTOS.
4.1. Métodos de análisis: Euleriano y Lagrangiano. 38

4.2. Aplicación de la ecuación de continuidad: criterios de incompresibilidad. 40
4.3. Aplicación de las ecuaciones de continuidad y de BERNOULLI: descarga de depósitos. 42
4.4. Aplicación de la ecuación de BERNOULLI i: flujo no viscoso entre discos horizontales. 45
4.5. Aplicación de la ecuación de BERNOULLI i no estacionario: oscilaciones en un tubo en U. 48
4.6. Aplicación de la ecuación de BERNOULLI con aceleración de arrastre: bomba rotativa. 50
4.7. Aplicación de las ecuaciones de NAVIER-STOKES: flujo de COUETTE-POISEUILLE. 53
4.8. Aplicación de la ecuación de NAVIER-STOKES: flujo viscoso entre discos horizontales. 55
4.9. Aplicación de la ecuación de Energía: distribución de temperaturas en flujo de Poiseuille. 59
ANEXOS.
A1: Ecuaciones de conservación en coordenadas cartesianas, cilíndricas y esféricas
A2: Newton, Berrnoulli, Euler, Lagrange, Laplace, Poisson, Navier, Cauchy, Stokes, Fraude, Mach, Reynolds,
Weber, Prandtl, von Karman.
_________________________________________________________________________________________________________________
1. CINEMÁTICA Y DINÁMICA DE FLUIDOS

1.1. Métodos de análisis en Mecánica de Fluidos: Lagrangiano y Euleriano.
1.2. Tipos de análisis en Mecánica de Fluidos: Diferencial, Dimensional e Integral.
1.3. Cinemática de Fluidos: propiedades del vector velocidad.
1.3.1. Derivada sustancial o material.
1.3.2. El operador NABLA: gradiente, divergencia y laplaciana.
1.4. Dinámica de Fluidos: fuerzas macroscópicas sobre los fluidos.
1.4.1. Fuerzas de volumen.
1.4.2. Fuerzas de superficie.
1.4.3. Tensor de tensiones.
1.5. Tipos de flujo.
1.1. MÉTODOS EN MECÁNICA DE FLUIDOS: LAGRANGIANO Y EULERIANO.
El ámbito general de la Mecánica de Fluidos, es la interacción entre fluidos y su entorno. Además el

fluido esta constituido por una sucesión continua de partículas que interaccionan entre si y entre los contornos.
La partícula fluida esta formada por una sucesión continua de puntos materiales, que integran un volumen
infinitesimal; que en el proceso de fluir, se deforma por la interacción con el resto de fluido, permaneciendo su
masa y su volumen elementales constantes, es decir, la densidad en toda la extensión de la partícula fluida es
constante.
Como metodología de estudio se dispone de dos alternativas:
- La identificación de cada partícula fluida y su seguimiento en el tiempo; es decir, hay que

determinan la posición de la partícula en función del tiempo, además de conocer las magnitudes
asociadas a cada partícula. Este el método de LAGRANGE, que es el usado normalmente en
Mecánica de Sólidos.
- En Mecánica de Fluidos, es suficiente conocer el valor de las propiedades en los diversos puntos
del campo fluido a lo largo del tiempo, con independencia de la partícula que lo ocupa en un
instante determinado; ésta es la base del método de EULER, en el que las magnitudes de las
propiedades de una partícula fluida en un determinado instante, vienen dadas por los valores de las
propiedades en el punto que es ocupado por la partícula en el citado instante. En el método
Euleriano, se deben determinar los “campos de propiedades”; así, el campo de presiones, es la
expresión espacial y temporal de la presión: p=p(x,y,z,t), con lo que una partícula que en un
instante “t0”, ocupe una posición (x0, y0, z0), tiene una presión dada por el campo de presiones:
p=p(x0,y0,z0,t0)
Al movimiento de un fluido se le denomina flujo, y en su análisis es interesante tener algún tipo de

representación. Cada método de análisis utiliza diferentes procedimientos de representación:
- En el método lagrangiano, se definen las trayectorias de las partículas como lugar geométrico de las
diferentes posiciones temporales de las partículas. La trayectoria de una partícula es el lugar geométrico
de las posiciones sucesivas, a lo largo del tiempo, de la partícula, que en el instante inicial (t0) estaba en la
r
posición inicial ( r0 ).
A(t2)
A(t1)
A(t3)
A(t0)
Fig.1.1.: Trayectoria de una partícula A a lo largo del tiempo: Método Lagrangiano
_________________________________________________________________________________________________________________
- En el método euleriano, se definen las líneas de corriente, que son las envolventes tangenciales de los
vectores velocidad de todas las partículas en un instante determinado.
v (x1,y1,z1,t0)
B(t0)
v (x0,y0,z0,to)
(x1,y1,z1)
A(t0) (x0,y0,z0)
Fig. 1.2. Línea de Corriente en un instante dado (t0) en distintos puntos ocupados por distintas partículas A,B,...
El vector unitario de la velocidad, es el vector unitario de la dirección tangencial, en determinado en un

punto de la línea de corriente; con lo que las ecuaciones diferenciales que dan lugar a las líneas de
corriente son:
dx dy dz
= =
u ( x, y, z, t ) v(x, y, z, t) w(x, y, z, t) (1.)
1.2. TIPOS DE ANÁLISIS EN MECÁNICA DE FLUIDOS.
La dinámica de fluidos trata del movimiento de los fluidos, a lo que se denomina flujo, y de sus
interacciones con el entorno. En el estudio de flujos hay que analizar el estado de movimiento del fluido,
definido por las ecuaciones de conservación (leyes fundamentales en el movimiento de fluidos), por las
ecuaciones de constitución (leyes del comportamiento del fluido) y por las condiciones de contorno (impuestas
por la geometría y el entorno).
Las ecuaciones de conservación y de constitución, junto con las condiciones de contorno, aplicadas a cada
una de las partículas del fluido, dan sistemas de ecuaciones diferenciales, cuya resolución lleva a definir el flujo,
en cuanto al campo de velocidades (cinemática) y al campo de fuerzas (dinámica). Este tipo de análisis
diferencial, da sistemas de ecuaciones simultaneas en derivadas parciales, que son de difícil o imposible
resolución; aunque pueden encontrarse soluciones analíticas con hipótesis restrictivas y en determinados casos,
en donde se pueden obtener soluciones parciales por cálculo numérico, utilizando las técnicas actuales de
simulación que constituye la mecánica de fluidos computacional (CFD: computational fluid mechanics), en
donde las derivadas se sustituyen por relaciones algebraicas en un número finito de puntos del flujo (mallado).
Si lo que se pretende obtener, no es el estado de movimiento del fluido, sino sus efectos sobre una
determinada región del flujo, se puede establecer otro tipo de análisis que evalúe las características globales del
flujo: caudales, fuerzas, momentos, potencias,... A la región de estudio, en donde se consideran las
interrelaciones entre entorno y flujo, se le denomina volumen de control; las modificaciones sobre el entorno
que introduce el flujo en su entrada-residencia-salida del volumen de control, o que el entorno introduce en las
propiedades del flujo, vienen determinadas por las ecuaciones integrales de conservación aplicadas al sistema
aislado entorno-volumen de control. Este método de análisis integral, se fundamenta en las ecuaciones
integrales que dan las velocidades de variación de las propiedades del fluido a su paso por el volumen de control.
Cuando el flujo es complejo y el análisis diferencial no aporta soluciones (por ser insuficientes las
ecuaciones o porque la resolución de los sistemas en derivadas parciales no es posible); y debido a que el análisis
integral da resultados globales; es necesario recurrir a un análisis experimental, en donde los resultados se
obtienen a partir de las magnitudes medidas en los experimentos. En este método de análisis aparecen dos
problemas propios: el gran número de variables que intervienen en la descripción del flujo y la imposibilidad, en
ciertos casos, de ensayar en condiciones reales. Para abordar estos problemas, se dispone del análisis
_________________________________________________________________________________________________________________
dimensional que permite reducir el número de variables y la teoría de modelos, con la que se correlacionan los
resultados experimentales de un modelo con los que tendría su prototipo.
El análisis diferencial puede ser utilizado para cualquier tipo de flujo, pero la dificultad de establecer y
resolver sistemas de ecuaciones diferenciales, limita el método; también el análisis experimental puede aplicarse
a cualquier flujo, pero las dificultades inherentes a las técnicas experimentales, presupuesto y universalidad, son
las que limitan el método; en cuanto al análisis integral, si que aporta resultados en el estudio técnico de flujos,
pero siempre de magnitudes globales.
El análisis diferencial comenzó con EULER y LAGRANGE en el siglo XVIII, el análisis dimensional tuvo
sus primeros pasos con RAYLEIG a finales del siglo XIX, y el análisis integral, aunque propuesto por EULER, se
desarrollo a mediados del siglo XX. En la actualidad las potentes técnicas de cálculo numérico, implementadas
en ordenadores cada vez más rápidos, han hecho posible, el resurgimiento del análisis diferencial, en cuanto a la
posibilidad de resolución de flujos cada vez más complejos, en donde se consideran los efectos viscosos. En
cuanto al análisis experimental, el desarrollo de sensores específicos (piezoeléctricos de presión,
extensiométricos de fuerza,...) y de técnicas cada vez menos intrusivas (velocimetría laser-doppler, hilo
radiante,...), esta aportando medidas cada vez más fiables.
1.3. CINEMÁTICA DE FLUIDOS: PROPIEDADES DEL VECTOR VELOCIDAD.

El término cinemática esta asociado a las propiedades del campo de velocidades. En la descripción
r
Euleriana del flujo, la velocidad local del fluido v(x, y, z, t) es la variable fundamental. Se pueden tener dos
casos extremos:
- flujo estacionario, cuando la velocidad es independiente del tiempo, con lo que en un determinado
punto, la velocidad (y en general cualquier propiedad) no varía con el tiempo, es decir la velocidad
solo depende de la posición:
r
r r ∂v
v = v(x, y, z) ⇒ =0
∂t (2.)
- flujo uniforme, cuando la velocidad no depende de la posición, con lo que en un determinado

instante, la velocidad (y en general cualquier propiedad) es la misma en todos los puntos del campo
fluido, es decir la velocidad solo depende del tiempo, y su gradiente es nulo.
r r r
v = v(t) ⇒ ∇v = 0
(3.)
1.3.1. DERIVADA SUSTANCIAL O MATERIAL.
Consideremos una propiedad intensiva del campo fluido, por ejemplo, la presión termodinámica. El
campo de presiones, es tetradimensional: p=p(x,y,z,t), y para su obtención, hay que resolver la ecuación
diferencial:
∂p ⎛ ∂p ∂p ∂p ⎞
dp = p(x + dx, y + dy, z + dz, t + dt) − p(x, y, z, t) = ⋅ dt + ⎜ ·dx + ·dy + ·dz ⎟
∂t ⎝ ∂x ∂y ∂z ⎠
para lo cual, hay que determinar las 4 variaciones parciales de la presión: la temporal y las tres espaciales.
Una vez conocido el campo de presiones, si se sigue a una partícula, la magnitud de su presión, viene
dada por el valor del campo de presiones, que es función de la posición en que se encuentre en un determinado
instante.
La variación temporal del campo de presiones, viene dada por:
dp ∂p dx ∂p dy ∂p dz ∂p
= + ⋅ + ⋅ + ⋅
dt ∂t dt ∂x dt ∂y dt ∂z
_________________________________________________________________________________________________________________
en las expresiones anteriores, los espacios recorridos elementales: dx. dy, dz, representan los incrementos
espaciales, en los que se conoce el campo fluido. Si los citados incrementos, los marcan las partículas en su
movimiento a una determinada velocidad1; con lo que los respectivos incrementos espaciales son:
dx = u·dt; dy = v·dt; dz = w·dt

r r r r
en donde u,v,w son las componentes cartesianas del vector velocidad: v = u⋅ i + v⋅ j + w⋅k
con todo, se tiene, que para un observador que se mueve con la velocidad del fluido, la variación temporal de la
presión es:
dp ∂p ∂p ∂p ∂p
= + u⋅ + v⋅ + w⋅
dt drr = vr ⋅dt ∂t ∂x ∂y ∂z (4.)
en notación compacta, se tiene:

dp ∂p r
= + v ⋅ ∇p
dt r
dr r
=v ∂t
dt (5.)
∂ r
en donde el operador + v ⋅ ∇ , se denomina derivada sustancial o material, y representa la velocidad de
∂t
cambio de una determinada propiedad de una partícula de fluido en su movimiento. Se suele denotar, o bien por
dp D
r , o bien por . Así, en el caso consideradode la presión termodinámica, su derivada material es:
dt dr = vr Dt
dt
Dp ∂p r
= + v ⋅ ∇p
Dt ∂t (6.)
∂p
en donde la variación local es: , que evalúa la variación temporal de la presión en un determinado punto
∂t
r
y la variación convectiva es: v ⋅ ∇p , que evalúa la variación de la presión (por unidad de tiempo) de un punto a
otro, al moverse el fluido por un gradiente de presiones.
el termino ∇p se denomina gradiente de presión, y su expresión depende del sistema de coordenadas:
∂p r ∂p r ∂p r
coordenadas cartesianas: ∇p = ⋅ i + ⋅ j + ⋅k
∂x ∂y ∂z
∂p r 1 ∂p r ∂p r
coordenadas cilíndricas: ∇p = ⋅ u r + ⋅ + ⋅ uz
∂r r ∂θ ∂z
Si consideramos una propiedad vectorial, como por ejemplo la misma velocidad, su derivada
material, viene determinada por:
r r r r r r
Dv ∂v ⎛ ∂v ∂v ∂v ⎞ ∂v r r
Coordenadas cartesianas: = +⎜u + v + w ⎟ = + v ⋅ ∇v
Dt ∂t ⎝ ∂x ∂y ∂z ⎠ ∂t
La aceleración del fluido, es la derivada material de la velocidad, y es suma de una aceleración local y
una aceleración convectiva:
r r
r Dv ∂v r r
a= = + v ⋅ ∇v
Dt ∂t (7.)
r
∂v
aceleración local
∂t
r r ⎛ ∂u ∂u ∂u ⎞ r ⎛ ∂v ∂v ∂v ⎞ r ⎛ ∂w ∂w ∂w ⎞ r
aceleración convectiva v ⋅ ∇v = ⎜ u + v + w ⎟ ⋅ i + ⎜ u + v + w ⎟ ⋅ j + ⎜ u +v +w ⎟⋅k
⎝ ∂x ∂y ∂z ⎠ ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂x ∂y ∂z ⎠
1
es como si el campo fluido se estuviera marcando por un observador que se mueve con la velocidad del fluido
_________________________________________________________________________________________________________________
El gradiente de velocidad, es un tensor, cuyas componentes cartesianas son:
⎛ ∂ ⎞ ⎛ ∂u ∂v ∂w ⎞
⎜ ∂x ⎟ ⎜ ∂x ∂x ∂x ⎟⎟
⎜ ⎟ ⎜
∂ ∂u ∂v ∂w ⎟
v = ⎜⎜ ⎟⎟ ( u v w ) = ⎜⎜
r
∂y ∂y ∂y ∂y ⎟
⎜ ⎟ ⎜ ⎟
⎜ ∂ ⎟ ⎜ ∂u ∂v ∂w ⎟
⎜ ⎟ ⎜ ⎟
⎝ ∂z ⎠ ⎝ ∂z ∂z ∂z ⎠ (8.)
Si el sistema de referencia es no inercial, para obtener la aceleración absoluta, hay que añadir a la
aceleración relativa, calculada anteriormente, la aceleración de arrastre del sistema de referencia no inercial2,
respecto a un sistema de referencia inercial:
r r
r Dv ∂v r r
a relativa = = + v ⋅ ∇v
r r r Dt ∂t
aabsoluta = a relativa + aarrastre r
r r dω r r r r r r
aarrastre = a0 + SRNI × r + ωSRNI × ( ωSRNI × r ) + 2ωSRNI × v
dt
en donde la aceleración de arrastre, tiene 4 sumandos, respectivamente:

r
a0 = aceleración lineal del sistema de referencia no inercial (SRNI) respecto al inercial
r
dωSRNI r
×r = aceleración tangencial de la partícula
dt
r r r
ωSRNI × ( ωSRNI × r ) = aceleración centrípeta de la partícula
r r
2ωSRNI × v = aceleración de CORIOLIS.
r r
en ωSRNI es la velocidad de giro del SRNI, respecto al sistema de referencial inercial, y r es su vector de posición
en el SRNI
1.3.2. EL OPERADOR NABLA.
El operador NABLA: ∇ , se utiliza para compactar las expresiones cinemáticas. Es un operador

vectorial, que se puede aplicar a un escalar o a un vector, representando su gradiente.
NABLA es un operador vectorial, que se puede operar con otro vector o con un tensor. Así el producto
escalar del operador NABLA representa la divergencia: si el producto escalar es con otro vector, se tiene la
divergencia del vector, y si el producto escalar es con un tensor se tiene la divergencia del tensor. El producto
vectorial del operador NABLA con un vector, representa el rotacional del vector:
∂ r ∂ r ∂ r
OPERADOR NABLA en coordenadas cartesianas: ∇ = ⋅ i + ⋅ j + ⋅k
∂x ∂y ∂z
∂ r 1 ∂ r ∂ r
OPERADOR NABLA en coordenadas cilíndricas: ∇= ⋅ ur + ⋅ + ⋅ uz
∂r r ∂θ ∂z
∂ r 1 ∂ r 1 ∂ r
OPERADOR NABLA en coordenadas esféricas: ∇= ⋅ ur + ⋅ uθ + ⋅ uϕ
∂r r ∂θ rsenθ ∂ϕ
2
Un sistema de referencia NO INERCIAL, se mueve con velocidad no constante, respecto a un sistema de
referencia INERCIAL (a velocidad constante).
_________________________________________________________________________________________________________________
GRADIENTE: la aplicación del operador NABLA a una magnitud, representa el gradiente de la magnitud a lo
largo del espacio: es decir, su variación de un punto a otro.
Cuando el operador ∇ se aplica a un escalar se obtiene un vector que representa la variación del
escalar desde un punto a otro del flujo y que se denomina gradiente de la magnitud escalar.
Cuando el operador ∇ se aplica a un vector se obtiene un tensor que representa la variación del vector
desde un punto a otro del flujo, que se denomina gradiente de la magnitud vectorial.
cartesianas = i+ j+ k
∂x ∂x ∂z
GRADIENTE DE PRESIÓN: ∇p =
∂p r 1 ∂p r ∂p r
cilíndricas = ur + uθ + u z
∂r r ∂θ ∂z
∂T r ∂T r ∂T r
∂x ∂x ∂z
GRADIENTE DE TEMPERATURA: ∇T =
∂T r 1 ∂T r ∂T r
cilíndricas = ur + uθ + uz
∂r r ∂θ ∂z
∂ρ r ∂ρ r ∂ρ r
∂x ∂x ∂z
GRADIENTE DE DENSIDAD: ∇ρ =
∂ρ r 1 ∂ρ r ∂ρ r
cilíndricas = ur + uθ + u z
∂r r ∂θ ∂z
⎛ ∂ ⎞ ⎛ ∂u ∂v ∂w ⎞
⎜ ⎟ ⎜ ⎟
⎜ ∂x ⎟ ⎜ ∂x ∂x ∂x ⎟
r ⎜ ∂ ⎟ ⎜ ∂u ∂v ∂w ⎟
GRADIENTE DE VELOCIDAD: ∇v = ⎜ ⎟ ( u v w ) = ⎜
∂y ∂y ∂y ∂y ⎟
⎜ ⎟ ⎜ ⎟
⎜ ⎟∂ ∂
⎜ u ∂v ∂w ⎟
⎜ ⎟ ⎜ ⎟
⎝ ∂z ⎠ ⎝ ∂z ∂z ∂z ⎠
DIVERGENCIA: el producto escalar del operador ∇ con un vector, es un escalar que se denomina
divergencia del vector. El producto escalar del operador NABLA con el vector velocidad, es la divergencia de
la velocidad; que representa la velocidad de variación de la densidad por unidad de densidad3 (que también es
igual a la velocidad de variación del volumen por unidad de volumen)
∂u ∂v ∂w
cartesianas = + +
r ∂x ∂x ∂z
DIVERGENCIA DE VELOCIDAD: ∇⋅v =
1 ∂ ( r ⋅ v r ) 1 ∂v θ ∂v z
cilíndricas = + +
r ∂r r ∂θ ∂z
El producto escalar del operador ∇ con un tensor, es un vector que se denomina divergencia del
tensor; así la divergencia del tensor de tensiones es el vector de tensiones de contacto (fuerzas de superficie
por unidad de área) sobre una determinada partícula fluida:
⎛ ∂τ ∂τ ∂τ ⎞ r ⎛ ∂τ ∂τ ∂τ ⎞ r ⎛ ∂τ ∂τ ∂τ ⎞ r
∇ ⋅ τ = ⎜ xx + yx + zx ⎟ ⋅ i + ⎜ xy + yy + zy ⎟ ⋅ j + ⎜ xz + yz + zz ⎟ ⋅ k
⎝ ∂x ∂y ∂z ⎠ ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂x ∂y ∂z ⎠
r −1 dρ
Ver ecuación de continuidad: ∇ ⋅ v =
3
ρ dt
_________________________________________________________________________________________________________________
ROTACIONAL: el producto vectorial del operador ∇ con un vector, es un vector que se denomina rotacional.
El producto vectorial del operador NABLA con el vector velocidad es el rotacional de la velocidad que
representa la velocidad angular local de una partícula fluida (el doble), y que se denomina vorticidad:
r r r
ROTACIONAL DE VELOCIDAD: ∇ × v = Ω = VORTICIDAD =2· ω
r ⎛ ∂w ∂v ⎞ r ⎛ ∂u ∂w ⎞ r ⎛ ∂v ∂u ⎞ r
Cartesianas ∇ × v = ⎜ − ⎟⋅ i + ⎜ − ⎟⋅ j + ⎜ − ⎟⋅k
⎝ ∂y ∂z ⎠ ⎝ ∂z ∂x ⎠ ⎝ ∂x ∂y ⎠
r ⎛ 1 ∂v z ∂v θ ⎞ r ⎛ ∂v r ∂v z ⎞ r ⎛ 1 ∂ ( r ⋅ v θ ) 1 ∂v r ⎞ r
Cilíndricas ∇ × v = ⎜ − ⎟ ⋅ ur + ⎜ − ⎟ ⋅ uθ + ⎜ − ⎟ ⋅ uz
⎝ r ∂θ ∂z ⎠ ⎝ ∂z ∂r ⎠ ⎝r ∂r r ∂θ ⎠
LAPLACIANA: la divergencia de un gradiente, se denomina laplaciana; que cuando se aplica a un escalar da

lugar a otro escalar, y cuando se aplica a un vector da otro vector. Como ejemplos significativos consideraremos
la laplaciana de temperatura y la laplaciana de velocidad
LAPLACIANA DE TEMPERATURA: ∇ ⋅ ( ∇T ) = ∇ 2 T
∂ 2T ∂ 2T ∂ 2T
Cartesianas ∇ 2T = + +
∂x 2 ∂y 2 ∂z 2
1 ∂T ∂ 2 T 1 ∂ 2 T ∂ 2 T
Cilíndricas ∇ 2T = + + +
r ∂r ∂r 2 r 2 ∂θ 2 ∂z 2
1 ∂ ⎛ 2 ∂T ⎞ 1 ∂⎛ ∂T ⎞ 1 ∂ 2T
Esféricas ∇2T = ⎜r ⎟+ 2 ⎜ senθ ⎟+ 2 2
r ∂r ⎝ ∂r ⎠ r senθ ∂r ⎝
2
∂θ ⎠ r sen θ ∂ϕ2
r r
LAPLACIANA DE VELOCIDAD: ∇ ⋅ ( ∇v ) = ∇ 2 v
r ⎛ ∂2u ∂2u ∂2u ⎞ r ⎛ ∂2 v ∂2 v ∂2 v ⎞ r ⎛ ∂2w ∂2w ∂2 w ⎞ r

Cartesianas: ∇2 v = ⎜ 2 + 2 + 2 ⎟ ⋅ i + ⎜ 2 + 2 + 2 ⎟ ⋅ j + ⎜ 2 + 2 + 2 ⎟ ⋅ k
⎝ ∂x ∂y ∂z ⎠ ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂x ∂y ∂z ⎠
r ⎛ 1 ∂v r ∂ 2 v r 1 ∂ 2 v r ∂ 2 v r v r 2 ∂v θ ⎞ r
∇2 v = ⎜ + 2 + 2 + 2 − 2 − 2 ⎟ ⋅ ur +
⎝ r ∂r ∂r r ∂θ2 ∂z r r ∂θ ⎠
⎛ 1 ∂v θ ∂ 2 v θ 1 ∂ 2 v θ ∂ 2 v θ v θ 2 ∂v r ⎞ r
Cilíndricas: +⎜ + 2 + 2 + 2 − 2 + 2 ⎟ ⋅ uθ +
⎝ r ∂r ∂r r ∂θ2 ∂z r r ∂θ ⎠
⎛ 1 ∂v z ∂ 2 v z 1 ∂ 2 v z ∂ 2 v z ⎞ r
+⎜ + 2 + 2 + 2 ⎟ ⋅ uz
⎝ r ∂r ∂r r ∂θ2 ∂z ⎠
r ⎡ 2⎛ ∂v v 1 ∂v ϕ ⎞ ⎤ r
∇ 2 ⋅ v = ⎢∇ 2 v r − 2 ⎜ v r + θ + θ + ⎟ ⎥ ⋅ ur +
⎣ r ⎝ ∂θ tgθ senθ ∂ϕ ⎠ ⎦
⎡ 2 ∂v 1 ⎛ ∂v ⎞ ⎤ r
Esféricas: + ⎢∇ 2 v θ + 2 r − 2 2 ⎜ v θ + 2 cos θ ϕ ⎟ ⎥ ⋅ u θ +
⎣ r ∂θ r sen θ ⎝ ∂ϕ ⎠ ⎦
⎡ 1 ⎛ ∂v ∂v ⎞ ⎤ r
+ ⎢∇ 2 v ϕ − 2 2 ⎜ v ϕ − 2senθ r − 2 cos θ θ ⎟ ⎥ ⋅ uϕ
⎣ r sen θ ⎝ ∂ϕ ∂ϕ ⎠ ⎦
_________________________________________________________________________________________________________________
1.4. DINÁMICA DE FLUIDOS: FUERZAS MACROCÓPICAS.
En el Análisis Diferencial de Fluidos, hemos considerado como volumen de control a la partícula

fluida, que es una porción de fluido de dimensiones infinitesimales y arbitrarias. El tamaño esta en relación a las
dimensiones del equipo de medida y del tiempo de medición. En todo caso los valores de las magnitudes son
medias temporales y espaciales, referidas a un intervalo de tiempo elemental y al conjunto de moléculas que
integran la partícula fluida.
Para su análisis, la partícula fluida se aísla del fluido que la envuelve mediante las superficies de
contacto partícula-fluido, y en este estado de equilibrio de la partícula aislada, se analizan las fuerzas que la
mantienen en equilibrio; estas fuerzas se dividen en tres tipos: fuerzas de volumen, fuerzas de superficie y
fuerzas de inercia.
1.4.1. FUERZAS DE VOLUMEN.
Las fuerzas de volumen, son fuerzas débiles de largo alcance, actúan sobre cada elemento de volumen
del fluido, y son debidas a campos de fuerzas externos. La fuerza elemental, que ejerce el campo externo, sobre
el elemento de volumen, por unidad de volumen es:
r
dFV r
= fV
dV
El campo externo de fuerzas, más habitual, es el gravitatorio terrestre, en donde la fuerza por unidad
r
r
de volumen, viene determinada por la densidad de la partícula fluida y la aceleración gravitacional ( g = − gk )
r r
f V = ρg
Si la partícula fluida, tiene una determinada densidad de carga (ρq=dq/dV), y esta inmersa en un campo
r
eléctrico de intensidad E , la fuerza elemental por unidad de volumen es:
r r
fV = ρ q E
Si la partícula fluida, tiene una determinada densidad de carga (ρq=dq/dV), y esta inmersa en un campo
r
magnético de intensidad B , la fuerza elemental por unidad de volumen es:
r r r
fV = ρ q ⋅ v × Β
En general, si la partícula cargada, esta en un campo electromagnético, se tiene la fuerza de LORENTZ:
r r r r
(
fV = ρ q E + v × Β )
Genéricamente, consideraremos como expresión de la fuerza elemental por unidad de volumen, debida
a la acción de un campo de fuerzas externos:
r
dFV r r
= fV = ρ ⋅ g
dV (9.)
r
Si el sistema de referencia es no inercial, a la aceleración absoluta g , hay que restarle la aceleración de
r
arrastre a SRNI (del sistema de referencia no inercial respecto a un sistema de referencia inercial), para poder
establecer el equilibrio de todas las fuerzas, en el sistema de referencia no inercial.
r
dFV r r r
= fV = ρ ⋅ ( g - aSRNI )
dV
r
r r dωSRNI r r r r r r
a SRNI = a0 + × r + ωSRNI × ( ωSRNI × r ) + 2ωSRNI × v
dt
_________________________________________________________________________________________________________________
1.4.2. FUERZAS DE SUPERFICIE.
Las fuerzas de contacto, que sobre las superficies de una partícula fluido en movimiento, ejerce el fluido
que la rodea, se denominan fuerzas de superficie. Son debidas a la agitación molecular y a la interacción
molecular, por lo que son apreciables sólo a distancias del orden de la interacción molecular.
Se puede demostrar que la resultante de las fuerzas elementales de superficie (por unidad de volumen),
sobre una partícula fluida, viene determinado por la divergencia de un tensor asociado al punto por el que pasa la
partícula, y que se denomina tensor de tensiones:
r
dFS r
= fS = ∇ ⋅ σ
dV (10.)
1.4.3. TENSOR DE TENSIONES.
En un determinado elemento de área de la partícula fluida, las fuerzas de superficie son proporcionales
al área de contacto con el resto del fluido; y además por el principio de acción-reacción, la fuerza que el resto de
fluido hace sobre la partícula a través del citado elemento de área, está equilibrada por la correspondiente fuerza
que la partícula hace sobre el resto de fluido.
A la fuerza de contacto elemental por unidad de área de contacto, se le denomina tensión. Sobre un
determinado elemento de área, la tensión (resultante elemental de las fuerzas de superficie, por unidad de área de
contacto), se pueden descomponer en tres componentes elementales ortogonales: una componente en la dirección
normal (del elemento de área) y dos en direcciones tangenciales, es decir una tensión normal y dos tensiones
tangenciales.
Consideremos un punto material como intersección de 3 planos ortogonales, cada uno de los cuales esta
marcando una dirección normal y dos tangenciales, con lo que se tendrá en conjunto, en el citado punto material,
3 tensiones normales y 6 tensiones tangenciales, que se agrupan en el denominado tensor de tensiones:
⎛ σ xx τ yx τ yz ⎞
⎜ ⎟
σ = ⎜ τ xy σ yy τzy ⎟
⎜τ τ yz σzz ⎟⎠
⎝ xz
σzz
z τyz
τxz
yπz=cte. τzx
πx=cte. τyx
x πy=cte. σxx
τzy
σyy
τxy
Fig. 1.3. Tensor de tensiones en un punto.
_________________________________________________________________________________________________________________
Consideremos una partícula fluida elemental, en coordenadas cartesianas, de tamaño dx·dy·dz, en cuyo
centro de gravedad, se tiene el tensor de tensiones σ . El balance de fuerzas elementales de superficie, en la
dirección “x”, es:
(1) sobre las caras x=cte, si la tensión normal en el cdg es σxx , en la cara anterior, la fuerza de superficie en la
⎛ ∂σ dx ⎞
dirección “x” es: ⎜ σxx + xx ⎟ ( dy·dz ) y de sentido positivo (consideramos que el fluido que rodea la citada
⎝ ∂x 2 ⎠
cara anterior es más rápido que la propia cara), y en la cara posterior, la fuerza de superficie en la dirección “x”
⎛ ∂σ −dx ⎞
es ⎜ σxx + xx ⎟ ( dy·dz ) y de sentido negativo ( suponemos que el fluido que rodea la citada cara posterior
⎝ ∂x 2 ⎠
es más lento que la propia cara).
dy ⎛ ∂τ xx −dx ⎞
⎜ τ xx + ⎟ ( dy·dz )
z dx ⎝ ∂x 2 ⎠
dz
cdg τ xx
dx
y
⎛ ∂τ xx dx ⎞
⎟ ( dy·dz )
x
⎜ τ xx +
⎝ ∂x 2 ⎠
Fig.1.4. Partícula fluida. Fig. 1.5. Fuerzas “x” en caras x=cte,
(2) Análogamente sobre las caras y=cte, z=cte.
⎛ ∂τ xy dy ⎞
⎜ τ xy + ⎟ ( dx·dz )
⎝ ∂y 2 ⎠
⎛ ∂τ xz −dz ⎞
⎜ τ xz + ⎟ ( dx·dy )
⎝ ∂z 2 ⎠
τ xz
τ xy
⎛ ∂τ xz dz ⎞
⎜ τ xz + ⎟ ( dx·dy )
⎝ ∂z 2 ⎠
⎛ ∂τ xy −dy ⎞
⎜ τ xy + ⎟ ( dx·dz )
⎝ ∂y 2 ⎠
Fig. 1.6. Fuerzas “x” en caras y=cte. Fig. 1.7. Fuerzas “x” en caras z=cte.
⎛⎡ ∂σ xx dx ⎞ ⎛ ∂σ xx −dx ⎞ ⎤ ⎞
( dFS )x = ⎜ ⎢⎛⎜ σxx + ⎟ − ⎜ σ xx + ⎟ ⎟ ( dy ⋅ dz ) +
⎝ ⎣⎝ ∂x 2 ⎠ ⎝ ∂x 2 ⎠ ⎦⎥ ⎠
⎛ ⎡⎛ ∂τ dy ⎞ ⎛ ∂τ −dy ⎞ ⎤ ⎞
+ ⎜ ⎢ ⎜ τ xy + xy
⎜ ⎝ ⎟ − ⎜ τ xy + xy ⎟ ⎥ ⎟ ( dx ⋅ dz ) +
⎝⎣ ∂y 2 ⎠ ⎝ ∂y 2 ⎠ ⎦ ⎟⎠
⎛ ⎡⎛ ∂τ dz ⎞ ⎛ ∂τ xz −dz ⎞ ⎤ ⎞ ⎛ ∂σxx ∂τ xy ∂τ xz ⎞
+ ⎜ ⎢ ⎜ τ xz + xz ⎟ − ⎜ τ xz + ⎟ ⎥ ⎟ ( dx ⋅ dy ) = ... = ⎜ + + ⎟ ( dx·dy·dz )
⎝ ⎣⎝ ∂z 2 ⎠ ⎝ ∂z 2 ⎠ ⎦ ⎠ ⎝ ∂x ∂y ∂z ⎠
Análogamente, las componentes “y” y “z”, de las fuerzas de superficie sobre la partícula fluida son:
_________________________________________________________________________________________________________________
⎛ ∂τ yx ∂σ yy ∂τ yz ⎞
( dFS )y = ⎜ + + ⎟ ( dx·dy·dz )
⎝ ∂x ∂y ∂z ⎠
⎛ ∂τzx ∂τzy ∂σzz ⎞

( dFS )z = ⎜ + + ⎟ ( dx·dy·dz )
⎝ ∂x ∂y ∂z ⎠
Con todo lo anterior, la resultante de las fuerzas de superficie sobre la partícula fluida es:
r ⎡ ⎛ ∂σ ∂τ ∂τ ⎞ r ⎛ ∂τ ∂σ ∂τ ⎞ r ⎛ ∂τ ∂τ ∂σ ⎞ r ⎤
dFS = ⎢ ⎜ xx + xy + xz ⎟ ⋅ i + ⎜ yx + yy + yz ⎟ ⋅ j + ⎜ zx + zy + zz ⎟ ⋅ k ⎥·dV
⎣ ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂x ∂y ∂z ⎠ ⎦
con lo que la fuerza elemental de superficie, por unidad de volumen se puede expresar por:
r ⎛ σ xx τ yx τzx ⎞
dFS ⎛ ∂ ∂ ∂ ⎞ ⎜ ⎟
= ⋅ τ xy σ yy τ zy ⎟ = ∇ ⋅ σ
dV ⎜⎝ ∂x ∂y ∂z ⎟⎠ ⎜⎜
⎝ τxz τ yz σ zz ⎟⎠
Por el equilibrio dinámico de la partícula fluida, los momentos de las fuerzas elementales tangenciales
deben ser nulos, lo que lleva a la igual de las tensiones tangenciales cruzadas; es decir:
τxy = τyx
τxz = τzx
τyz = τzy
TENSOR DE LA PRESIÓN TERMODINÁMICA: si el fluido esta en reposo, o se considera como fluido ideal
(los coeficientes de transportes son idénticamente nulos), las únicas interacciones moleculares, son debidas a la
agitación térmica molecular, que tiene dos características fundamentales: no hay dirección privilegiada de los
esfuerzos, y estos son exclusivamente normales y de compresión; lo que lleva a que el tensor de tensiones en un
determinado punto, sea un tensor esférico, con tensiones exclusivamente normales, iguales a la presión
termodinámica:
⎛ −p 0 0 ⎞ ⎛ 1 0 0⎞
σ ideal = ⎜ 0 − p 0 ⎟ = − p ⋅ ⎜⎜ 0 1 0 ⎟⎟ = − p ⋅ 1
⎜ ⎟
⎜ 0 0 p ⎟⎠ ⎜ 0 0 1⎟
⎝ ⎝ ⎠
TENSOR DE TENSIONES VISCOSAS: cuando hay flujo, el fluido esta en movimiento, el tensor de tensiones
se puede descomponer en dos: uno que esférico, correspondiente a la presión termodinámica, y otro, que se
denomina tensor de tensiones viscosas, y que determina la diferencia entre las tensiones en un determinado
punto y las correspondientes a la presión termodinámica:
⎛ − p + τ xx τ yx τ zx ⎞ ⎛ − p 0 0 ⎞ ⎛ τ xx τ yx τ zx ⎞
⎜ ⎟ ⎜ ⎟
σ = ⎜ τ xy − p + τ yy τ zy ⎟ = ⎜⎜ 0 − p 0 ⎟⎟ + ⎜ τ xy τ yy τ zy ⎟ = − p ⋅ 1 + τ
⎜ τ
⎝ xz τ yz − p + τzz ⎟⎠ ⎜⎝ 0 0 0 ⎟⎠ ⎜⎝ τ xz τ yz τ zz ⎟⎠
en donde: σ es el tensor de tensiones; “p” es la presión termodinámica y τ el tensor de tensiones viscosas.
La resultante de las fuerzas de contacto sobre toda la partícula fluida, (Ec. 10.), se suele expresar, como
suma de la contribución de la presión termodinámica y las tensiones viscosas:
r
dFs r
= fS − ∇p + ∇ ⋅ τ
dV (11.)
_________________________________________________________________________________________________________________
1.5. TIPOS DE FLUJOS

Para poder acotar el estudio del movimiento de un fluido, se establecen las pertinentes restricciones, que
determinan los siguientes tipos de flujos:
- Flujo estacionario: ( ∂ #/ ∂t = 0 ∧ ∇ # ≠ 0) , en un determinado punto las propiedades del fluido no

varían con el tiempo, aunque puedan variar de un punto a otro (gradiente no nulo).
- Flujo uniforme: ( ∂ #/ ∂t ≠ 0 ∧ ∇ # = 0) en un determinado instante todas las partículas tienen la misma
velocidad en cualquier posición (gradiente nulo)
r
- Flujo transitorio o no estacionario: ( ∂ # / ∂t ≠ 0 ∧ ∇⋅# ≠ 0) , las propiedades del fluido varían con el
tiempo en cada punto y de un punto a otro.
r r
- Flujo irrotacional: ( ∇xv = 0 ≡ Ω = 0 ), el vector rotacional de velocidad es nulo y con ello la
vorticidad es nula.
r r
- Flujo rotacional: ( ∇xv ≠ 0 ≡ Ω ≠ 0 ), de vorticidad no nula.
r
- Flujo incompresible: ( ρ = cte. ≡ ∇ ⋅ v = 0 ); la densidad es constante en todos lo puntos y a lo largo
del tiempo, lo que lleva a que la divergencia de la velocidad sea nula, lo que suele expresar como flujo
adivergente.
- Flujo compresible: ( ρ ≠ cte. ), la densidad varia a lo largo del tiempo y del espacio.
- Flujo no viscoso: ( µ = 0 ), no hay transporte de cantidad de movimiento entre las partículas del flujo.
- Flujo viscoso: ( µ ≠ 0 ), hay interacción entre las partículas que constituyen el fluido, manifestándose
como intercambios de cantidad de movimiento, que dan lugar a fenómenos de disipación de energía, que
se denomina disipación viscosa.
- Flujo ideal: (µ=0 ∧ κ =0 ∧ δ=0 ), no hay interacción entre las partículas que constituyen el fluido, ni
de transporte de cantidad de movimiento (viscosidad), ni de transporte de calor (conductividad), ni de
transporte de masa (difusividad)
- Flujo laminar: en donde las fuerzas viscosas predominan sobre las de inercia; y en la interacción viscosa
con otras partículas, una determinada partícula de fluido no cambia su trayectoria, siendo arrastrada por la
acción del resto de partículas: frenada por partículas más lentas y acelerada por partículas más rápidas.
- Flujo turbulento: en donde las fuerzas de inercia predominan sobre las fuerzas viscosas; y en la
interacción viscosa con otras partículas, una determinada partícula es desplazada de su trayectoria por los
intercambios de cantidad de movimiento de otras partículas, además de ser arrastrada.
- Flujo interno: en donde el fluido esta confinado por límites; a partir de los campos de velocidades y de
presiones, se obtienen las fuerzas sobre los límites y la perdida de energía del fluido a su paso entre los
límites; un caso típico es el estudio de flujo en tuberías.
- Flujo externo: en donde el fluido rodea a un objeto; a partir de los campo de velocidades y presiones, se
obtienen se obtienen las fuerzas del fluido sobre el objeto; un caso típico es el estudio de perfiles
aerodinámicos.
_________________________________________________________________________________________________________________
2. ECUACIONES DE CONSTITUCIÓN
2.1. Comportamiento mecánico: tensor de velocidad de deformación.
2.2. Fluidos Stokesianos: tensor de tensiones viscosas.
2.3. Fluidos Newtonianos.
2.3.1. Ec. de NAVIER-POISSON.
2.3.3. Tensor de tensiones de un fluido newtoniano.
2.4. Comportamiento térmico.
2.4.1. Ecuaciones de Estado.
2.4.2. Ecuaciones de transmisión de calor.
2.1. COMPORTAMIENTO MECÁNICO: Tensor de velocidad de deformación.
En función de las hipótesis restrictivas, con las que se analiza el comportamiento de los fluidos reales,
se tienen las Ecuaciones de Constitución, que son inherentes a cada fluido analizado. El comportamiento
especifico de un determinado fluido, viene determinado por su comportamiento mecánico y su comportamiento
térmico.
El comportamiento mecánico del fluido, viene determinado por la relación entre las tensiones a las que
esta sometido y las velocidades de deformación que se producen por la acción de las tensiones mecánicas. Este
es el comportamiento inherente de los fluidos, es decir, las débiles fuerzas intermoleculares, hacen que cualquier
esfuerzo tangencial, deforme continuamente el fluido, originando el movimiento de las partículas o flujo. La
velocidad de deformación viene determinada por la magnitud del esfuerzo tangencial y de la capacidad de
transporte de cantidad de movimiento entre partículas, que es la propiedad mas importante, inherente al fluido, y
que se denomina viscosidad.
En el método Euleriano, en cada punto del flujo, la velocidad de deformación viene determinada por el
campo de velocidades. Cada punto del flujo, tiene asociado un valor del tensor de velocidades de deformación,
que marca la deformación unitaria de una partícula fluida a su paso por el citado punto. La deformación de una
determinada partícula en su movimiento por el campo fluido, viene determinada por las posibles variaciones
espaciales de la velocidad de cada uno de los puntos que la integran, es decir del gradiente de velocidad, que al
ser una magnitud tensorial (9 variaciones posibles), marca la misma condición tensorial a la velocidad de
deformación.
Se tienen dos tipos de deformación: las debidas a alargamientos o contracciones, provocadas por los
gradientes de las componentes de la velocidad en sus respectivas direcciones, y que se determinan por la
velocidad de la variación unitaria (por unidad de longitud); y las debidas a giros, provocados por los gradientes
de las componentes de la velocidad en direcciones perpendiculares a la propia componentes, y que se determinan
por la velocidad de variación angular.
r r
Consideremos, un caso muy simple, en donde v = v(y) ⋅ j , es decir, la única componente de la
velocidad, es en la dirección “y”, y además esa componente sólo varia en la propia dirección “y”. Si
consideramos una partícula elemental dx·dy·dz, al cabo de un tiempo elemental, se ha deformado (en este caso
sólo en la dirección “y”), teniendo que su velocidad de deformación unitaria (dilatación o contracción por
unidad de longitud y de tiempo) viene dada por:
⎛⎛ ∂v ⎞ ⎞
⎜ ⎜ v + dy ⎟ ⋅ dt − v ⋅ dt ⎟ / dy
dl / dy ⎝ ⎝ ∂y ⎠ ⎠ ∂v
ε& yy = yy = =
dt dt ∂y
_________________________________________________________________________________________________________________
que además representa la velocidad del aumento (o disminución) unitario de volumen :
dV(t + dt) dl yy ⋅ dx ⋅ dz
dV(t) dx ⋅ dy ⋅ dz dl / dy
= = yy = ε& yy
dt dt dt
dy dlyy
z
dz
dx
y
x
Fig.1.8. Dilatación cúbica: debida a la deformación lineal.
Si se tiene un campo de velocidades genérico: u=u(x,y,z), v=v(x,y,z), w=w(x,y,z), se obtienen las

correspondiente velocidades de dilataciones lineales unitarias, dadas por:
∂u ∂v ∂w
ε& xx = ε& yy = ε& zz =
∂x ∂y ∂z (12.)
Con lo que, la velocidad de dilatación cúbica, viene determinada por, la suma de las dilataciones
posibles en cada una de las tres direcciones; que es la divergencia de la velocidad:
1 dV ∂u ∂v ∂w r
= + + = ∇⋅v
V dt ∂x ∂y ∂z (13.)
Así un fluido de densidad constante, por conservación de masa, no hay variación del volumen, y por lo
tanto su flujo es adivergente.
Cada una de las tres dilataciones cúbicas, son la diagonal principal del tensor gradiente de velocidad; es
decir, el citado tensor esta marcando la dilatación cúbica que experimenta una partícula, cuando pasa por el
citado punto.
Hasta ahora hemos considerado deformaciones puramente lineales de dilatación o de contracción,

debidas a las variaciones de cada una de componentes del vector velocidad, en sus respectivas direcciones:
∂u ∂v ∂w
, , .
∂x ∂y ∂z
Consideremos el efecto de deformación, que tienen las variaciones cruzadas de las componentes de la
∂u ∂u ∂v ∂v ∂w ∂w
velocidad, es decir : , , , , , . Para lo cual, analicemos el caso más simple, en donde el vector
∂y ∂z ∂x ∂z ∂x ∂y
_________________________________________________________________________________________________________________
r r r
velocidad sea: v = u(y) ⋅ i + v(x) ⋅ j ; obteniéndose, que la deformación angular por unidad de tiempo, viene dada
dα 1 ⎛ ∂u ∂v ⎞
por: ε& xy = = ⎜ + ⎟ , que representa la velocidad de deformación angular, en un plano z = cte.
dt 2 ⎝ ∂y ∂x ⎠
∂u
y dydt
∂y
dα
∂v
∂u v+ dx
u + dy ∂x t+dt
v ∂y ∂v
dxdt
∂x
∂u
u+ dy
∂y
t
∂v
v+ dx
∂x v·dt
v
u u
u·dt
x
∂u ∂v
Fig.1.9. Velocidad de deformación angular, debida a gradientes velocidad cruzados: ,
∂y ∂x
Análogamente, para los gradientes cruzados, sin variación de “x”, se tiene que la velocidad de
1 ⎛ ∂v ∂w ⎞
deformación angular en un plano x = cte, es igual a: ε& yz = ⎜ + ⎟ ; y para los gradientes cruzados, sin
2 ⎝ ∂z ∂y ⎠
variación de “y”, se tiene que la velocidad de deformación angular en un plano y = cte, es igual a
1 ⎛ ∂u ∂w ⎞
ε& xz = ⎜ + ⎟
2 ⎝ ∂z ∂x ⎠
Con todo, se tiene que el tensor gradiente de velocidad, en un determinado punto, provoca que las
partículas que pasan por el citado punto, se deformen con una determinada velocidad, tanto longitudinal como
angularmente. El tensor, que marca las velocidades de deformaciones, es el tensor de velocidades de
deformación, y viene determinado por el tensor gradiente de velocidad.
En coordenadas cartesianas, el tensor de velocidad de deformación es:
⎛ ∂u 1 ⎛ ∂u ∂v ⎞ 1 ⎛ ∂u ∂w ⎞ ⎞
⎜ ⎜ + ⎟ ⎜ + ⎟⎟
⎜ ∂x 2 ⎝ ∂y ∂x ⎠ 2 ⎝ ∂z ∂x ⎠ ⎟
⎜ 1 ⎛ ∂u ∂v ⎞ ∂v 1 ⎛ ∂v ∂w ⎞ ⎟
ε& = ⎜ ⎜ + ⎟ + ⎟
⎜ 2 ⎝ ∂y ∂x ⎠ ∂y 2 ⎜⎝ ∂z ∂y ⎟⎠ ⎟
⎜ ⎟
⎜ 1 ⎛ ∂u + ∂w ⎞ 1 ⎛ ∂v ∂w ⎞ ∂w ⎟
⎜ ⎜ ⎟ ⎜ + ⎟ ⎟
⎝ 2 ⎝ ∂z ∂x ⎠ 2 ⎝ ∂z ∂y ⎠ ∂z ⎠ (14.)
_________________________________________________________________________________________________________________
2.2. FLUIDOS STOKESIANOS: TENSOR DE TENSIONES VISCOSAS.
STOKES, estableció que la diferencia de tensiones, entre un fluido viscoso y un fluido ideal, venia
determinada por una función tensorial del tensor de velocidad de deformación, que se denomina función de
tensiones viscosas (f); con ello el tensor de tensiones para un fluido Stokesiano esta integrado por dos términos:
el debido a la presión termodinámica y el debido a la viscosidad:
σ = − p ⋅ 1 + f ( ε& )
(15.)
Como se había visto anteriormente, el tensor de tensiones en un determinado punto, viene dado por los
esfuerzos normales y tangenciales, provocados por la interacciones entre partículas:
⎛ σ xx τ yx τ yz ⎞
⎜ ⎟
σ = ⎜ τ xy σ yy τzy ⎟
⎜ σzz ⎠⎟
⎝ τ xz τ yz (16.)
Las 3 componentes normales, se denotan por σxx, σyy, σzz. Las 6 componentes tangenciales, se denotan
por: τxy,τyx,τxz, τzx,τyz, τzy; siendo respectivamente iguales: τxy=τyx, τxz=τzx, τyz=τzy. Con lo que se tiene 6
tensiones distintas: 3 normales y 3 tangenciales.
El tensor de tensiones viscosas, es la diferencia entre el tensor de tensiones y el tensor esférico,

correspondiente a la presión termodinámica; se denota por τ , y tiene 3 componentes normales: τxx, τyy, τzz; y 6
componentes tangenciales, que coinciden con la del tensor de tensiones.
⎛ τ xx τ yx τ yz ⎞
⎜ ⎟
τ = ⎜ τ xy τ yy τzy ⎟
⎜ τ zz ⎠⎟
⎝ τ xz τ yz (17.)
2.3. FLUIDOS NEWTONIANOS.
El conocimiento de la función tensorial “f” de la Ec. 15., permitiría la determinación del campo de
tensiones viscosas a partir del campo de deformaciones, que a su vez depende del campo de velocidades. El
comportamiento más simple, es que la función sea lineal, en donde las tensiones viscosas sean proporcionales a
las velocidades de deformación; este es el comportamiento experimental dado por NAVIER y POISSON, para el
comportamiento reológico de un gran número de líquidos y de gases, que se denominan fluidos newtonianos.
2.3.1. ECUACIÓN DE NAVIER-POISSON.
En un fluido newtoniano, cada componente del tensor de tensiones viscosas, es función lineal de cada
componente del tensor de velocidades de deformación, pudiendo tener genéricamente 9x9=81 coeficientes de
proporcionalidad; pero si el medio es isótropo, los coeficientes se reducen exclusivamente a 3: uno para la
dirección normal, o coeficiente de viscosidad normal λ, y dos para cada dirección tangencial, o coeficientes de
viscosidad tangencial µ y µ’. Con lo que en un medio isótropo el tensor de tensiones viscosas se puede expresar
como:
r
τ = µ ⋅ ε& + µ '⋅ ε& + λ ( ∇ ⋅ v ) ⋅ 1
por simetría del propio tensor, se tiene que µ=µ’, con lo que µ+µ’=2µ; quedando como expresión del tensor de
tensiones para un fluido newtoniano, la denominada ecuación de NAVIER-POISSON:
r
τ = 2µ ⋅ ε& + λ ( ∇ ⋅ v ) ⋅ 1
(18.)
_________________________________________________________________________________________________________________
2.3.2. TENSOR DE TENSIONES EN UN FLUIDO NEWTONIANO.
El tensor de tensiones, lo hemos considerado suma de dos tensores: uno esférico debida a la acción de la
presión termodinámica, y otro debido a los esfuerzos viscosos, que en el caso de un fluido newtoniano, viene
determinado por la Ec. 18 de NAVIER-POISSON; con lo que el tensor de tensiones en este caso es:
r r
σ = − p ⋅ 1 + τ = ⎡⎢ − p ⋅ 1⎤⎥ + ⎡⎢ 2µ ⋅ ε& + λ ( ∇ ⋅ v ) ⋅ 1⎤⎥ = ⎡⎣ − p + λ ( ∇ ⋅ v )⎤⎦ ⋅ 1 + 2µ ⋅ ε&
⎣ ⎦ ⎣ ⎦ (19.)
Su expresión en coordenadas cartesianas es:
⎡ ⎛ ∂u ∂v ∂w ⎞ ∂u ⎛ ∂u ∂v ⎞ ⎛ ∂u ∂w ⎞ ⎤
⎢− p + λ⎜⎜ + + ⎟⎟ + 2µ µ⎜⎜ + ⎟⎟ µ⎜ + ⎟ ⎥
⎢ ⎝ ∂x ∂y ∂z ⎠ ∂x ⎝ ∂y ∂x ⎠ ⎝ ∂z ∂x ⎠ ⎥
= ⎢ ⎛ ∂u ∂v ⎞ ⎛ ∂u ∂v ∂w ⎞ ∂v ⎛ ∂v ∂w ⎞ ⎥
σ = ⎢ µ⎜⎜ + ⎟⎟ − p + λ⎜⎜ + + ⎟⎟ + 2µ µ⎜⎜ + ⎟⎟ ⎥
⎢ ⎝ ∂y ∂x ⎠ ⎝ ∂x ∂y ∂z ⎠ ∂y ⎝ ∂z ∂y ⎠ ⎥
⎢ ⎛ ∂u ∂w ⎞ ⎛ ∂v ∂w ⎞ ⎛ ∂u ∂v ∂w ⎞ ∂w ⎥
⎢ µ⎜ + ⎟ µ⎜⎜ + ⎟⎟ − p + λ⎜⎜ + + ⎟⎟ + 2µ ⎥
⎣⎢ ⎝ ∂z ∂x ⎠ ⎝ ∂z ∂y ⎠ ⎝ ∂x ∂y ∂z ⎠ ∂z ⎦⎥
(20.)
A partir de las componentes de la diagonal principal del tensor de tensiones , se puede establecer la
distinción entre presión termodinámica (p) y presión normal media ( p ):
- Presión termodinámica: es la propiedad termodinámica que genera un estado de tensión definido por un
tensor esférico de esfuerzos; es función del recorrido libre medio de las partículas, de la densidad de partículas y
de la agitación térmica de las partículas (fluctuaciones de posición debido a la temperatura).
- Presión normal media: es el valor medio de las tensiones normales sobre una partícula:
σx + σ y + σz
p=−
3
⎛ ∂u ∂v ∂w ⎞ ∂u
σ x = − p + λ⎜⎜ + + ⎟⎟ + 2µ
⎝ ∂x ∂y ∂z ⎠ ∂x ⎛ 2 ⎞⎛ ∂u ∂v ∂w ⎞
⇒ p = p − ⎜ λ + µ ⎟⎜⎜ + + ⎟
⎛ ∂u ∂v ∂w ⎞ ∂v ⎝ 3 ⎠⎝ ∂x ∂y ∂z ⎟⎠
σ x = − p + λ⎜⎜ + + ⎟⎟ + 2µ
⎝ ∂x ∂y ∂z ⎠ ∂y
⎛ ∂u ∂v ∂w ⎞ ∂w
σ z = − p + λ⎜⎜ + + ⎟⎟ + 2µ
⎝ ∂x ∂y ∂z ⎠ ∂z (21.)
al término (λ+2µ/3) se le denomina viscosidad bruta : κ=λ+2µ/3; con lo que la ecuación de la presión normal
r
media puede ponerse como: p = p − κ ⋅ ( ∇ ⋅ v )
.La igualdad entre ambas presiones (termodinámica y normal media), se verifica en los siguientes casos:
- viscosidad bruta nula: κ=0 ⇒ p = p ; obteniéndose por la condición de STOKES:
2 2
κ=0 ∧ κ = λ+ µ ⇒ λ=− µ ≡ CONDICIÓN DE STOKES
3 3 (22.)
r
- fluido incompresible: ρ=cte. ≡ ∇ ⋅ v =0 ⇒ p=p por la condición de adivergencia.
_________________________________________________________________________________________________________________
2.4. COMPORTAMIENTO TÉRMICO.
El comportamiento térmico del fluido viene determinado por las ecuaciones de estado y por la relación
entre flujo de calor y gradiente térmico.
2.4.1. ECUACIONES DE ESTADO.
Las ecuaciones de estado, son las relaciones entre las diversas variables de estado del fluido: presión,
temperatura, densidad, energía interna, entalpía, entropía,...:
- ecuación térmica de estado: f(p,T,ρ) = 0

- ecuación calórica de estado: û = û (T, p) = û (T, ρ) = û (p, ρ)
- ecuación entálpica de estado: h = h (T, p) = h (T, ρ) = h (p, ρ)
- ecuación entrópica de estado: s = s(T, p) = s(T, ρ) = s(p, ρ)
Cada una de las funciones que aparecen en las ecuaciones de estado anteriores, vienen determinadas por
las características intrínsecas de cada tipo de fluido. Así para el caso hipotético de un gas ideal, las ecuaciones de
estado se simplifican notablemente:
ecuación térmica: p=ρRT R = constante del gas4 = Ru/M

ecuación calórica: dû=cv dT cv = calor especifico a volumen constante
ecuación entálpica: dh = cp dT cp = calor especifico a presión constante
ecuación entrópica: ds = cp dT/T-Rdp/p
Si los calores específicos son constantes, se tienen las siguientes expresiones:
û - û 0 = cv (T-T0 )
h – h0 = cp (T-T0)
T p T ρ
s – s0 = c p ln − R ln = c v ln + R ln
T0 p0 T0 ρ0
Si el proceso es isentrópico (adiabático sin efectos disipativos), se tiene la relación entre las
variables, p,ρ,T:
γ
γ
p ⎛ ρ ⎞ ⎛ T ⎞ γ−1
=⎜ ⎟ =⎜ ⎟
p0 ⎝ ρ0 ⎠ ⎝ T0 ⎠
2.4.2. ECUACIONES DE TRANSMISIÓN DE CALOR.
La relación entre el flujo de calor y el gradiente térmico, viene determinada por el tipo de transmisión
de calor que tenga lugar: conducción, convección o radiación. Por simplicidad, considerando exclusivamente
transmisión de calor por conducción, se puede tomar como primera aproximación la ley de FOURIER de
conducción:
r
q = − κ ⋅ ∇T
(23.)
r
en donde “κ” es el coeficiente de conductividad térmica del fluido (W/mK), “ q ” el vector flujo de calor por
∂T r ∂T r ∂T r
unidad de área (W/m2), y “ ∇T ” el gradiente de temperatura (en coordenadas cartesianas: i+ j+ k)
∂x ∂y ∂z
4
Ru es la constante universal de los gases, cuyo valor es en el S.I. de 8310 J/kgK
_________________________________________________________________________________________________________________
3. ECUACIONES DE CONSERVACIÓN.
3.1. Ecuación diferencial de conservación de masa: ecuación de continuidad.
3.2. Ecuación diferencial de conservación de cantidad de movimiento: ecuación de movimiento de CAUCHY.
3.2.1. Fluido no viscoso: ecuación de EULER.
3.2.2. Flujo no viscoso en línea de corriente: ecuación de BERNOULLI.
3.2.3. Fluido newtoniano: ecuaciones de NAVIER-STOKES.
3.2.4. La función de corriente y la función potencial de velocidad.
3.3. Ecuación diferencial de conservación de la energía: ecuación de la energía.
3.3.1. Ecuación de energía interna.
3.3.2. Ecuación de entalpía.
3.3.3. Ecuación de entropía.
3.4. Condiciones de contorno.
3.1. ECUACIÓN DIFERENCIAL DE CONSERVACIÓN DE MASA: ECUACIÓN DE

CONTINUIDAD.
Los principios generales, que son válidos para cualquier tipo de entidad material, son una expresión
matemática de las leyes de conservación. En el caso del análisis diferencial en Mecánica de Fluidos se
consideran las siguientes leyes de conservación: conservación de masa, conservación de cantidad de movimiento
y conservación de energía. Consideraremos como entidad, la de una partícula fluida, que se aísla del resto del
fluido, y se le aplican las leyes de conservación.
Analizaremos en primer lugar la conservación de masa: utilizando el método “euleriano”,

consideraremos que la partícula es indeformable y que su volumen elemental ( dV = dx ⋅ dy ⋅ dz , en coordenadas
cartesianas) es siempre el mismo y esta siempre en la misma posición; se establece el siguiente balance de masa
entre dos instantes de tiempo “t” y “t+dt”:
masa ( t ) = ρ ⋅ dV
∂ρ
⎛ ∂ρ ⎞ dm = ⋅ dV ⋅ dt
masa ( t + dt ) = ⎜ ρ + dt ⎟ ⋅ dV ∂t
⎝ ∂t ⎠
La variación de masa en el volumen considerado durante el intervalo de tiempo “dt”, es debida al flujo
másico por las caras del elemento de volumen en el tiempo “dt”
⎡ ⎤
⎢ ⎛ & ∂m
dm = ⎢⎜ m −
& ⎞
dx ⎟
⎛
+ ⎜⎜ m
& −
∂m
& ⎞
dy ⎟⎟
⎛
+ ⎜m& −
∂m
& ⎞
dz ⎟ ⎥ dt = ... = − ⎡ ∂ (ρu ) + ∂ (ρv ) + ∂ (ρw ) ⎤ ⋅ dV ⋅ dt
∂x ∂y ⎠ caras ⎝ ∂z ⎠ caras ⎥ ⎢ ∂t ∂t ∂t ⎥⎦
⎝ ⎠ caras ⎝ ⎣
⎣⎢ dy ⋅ dz dx ⋅ dz dx ⋅ dy ⎦⎥
con las dos expresiones de la variación de masa de la partícula fluida considerada, se tiene:
∂ρ r
+ ∇ ⋅ (ρv ) = 0
∂t (24.)
Ecuación que se denomina de continuidad, porque en la ecuación de conservación de masa sólo se

requiere la derivabilidad de las funciones que dan la densidad y las componentes de la velocidad, es decir se
requiere su continuidad. Las funciones son continuas, porque estamos considerando como modelo del fluido, el
formado por una sucesión continua de partículas, es decir es un medio continuo.
_________________________________________________________________________________________________________________
w dy
v dx
∂ (ρ v dx dz ) ⎤
dz
⎡
u (ρ v dx dz )dt ⎢(ρ v dx dz ) + dy ⎥ ⋅ dt
⎣ ∂y ⎦
(masa entra − masa sale )caras dx⋅dz = (ρv dx dz )dt − ⎡⎢(ρv dx dz ) + ∂(ρv dx dz ) dy⎤⎥ dt = − ∂(ρ v ) ⋅ dx dy dz ⋅ dt = − ∂(ρv ) ⋅ dV ⋅ dt
⎣ ∂y ⎦ ∂y ∂t
Fig.1.10. Balance de masa en las caras dx dz, por variación del flujo másico de una cara a la otra
La ecuación de continuidad también se puede expresar en función de la derivada total de la densidad, al

r
descomponer la divergencia de ρ v en dos términos, y reagrupar la variación local de la densidad con su
variación convectiva, obteniendo:
∂ρ r ∂ρ r r ⎛ ∂ρ r ⎞ r dρ r
+ ∇ ⋅ ( ρv ) = + ( ρ ⋅ ( ∇ ⋅ v ) + v ⋅ ∇ρ ) = ⎜ + v ⋅ ∇ρ ⎟ + ρ ⋅ ( ∇ ⋅ v ) = + ρ ⋅ (∇ ⋅ v )
∂t ∂t ⎝ ∂t ⎠ dt
dρ r
+ ρ ⋅ (∇ ⋅ v ) = 0
dt (25.)
r r r r
En coordenadas cartesianas, y siendo el vector velocidad: v = u i + v j + wk , la ecuación de
continuidad es:
∂ρ ∂ (ρu ) ∂ (ρv) ∂ (ρw )
+ + + =0
∂t ∂x ∂y ∂z (26.)
r r r r
En coordenadas cilíndricas, y siendo el vector velocidad: v = v r ur + v θu θ + v z k , la ecuación de
continuidad es:
∂ρ 1 ∂ (r ρ v r ) 1 ∂ (r ρ v θ ) ∂ (ρv z )
+ + + =0
∂t r ∂r r ∂θ ∂z (27.)
Es interesante ver el efecto que tiene la divergencia de la velocidad sobre el flujo; de la Ec. 25.:
r 1 dρ
∇⋅v = −
ρ dt (28.)
Si la divergencia es positiva, la densidad disminuye y se tiene una expansión; si la divergencia es

negativa la densidad aumenta y se tiene una compresión; y si la divergencia es nula la densidad es constante.
_________________________________________________________________________________________________________________
En flujo estacionario compresible queda la ecuación diferencial:
∂ (ρ u) ∂ (ρ v) ∂ (ρ w)
cartesianas : + + =0
r ∂x ∂y ∂z
∇ ⋅ (ρ v ) = 0 ⇒
1 ∂ (r ρ v r ) 1 ∂ (r ρ v θ ) ∂ (ρ v z )
cilíndricas : + + =0
r ∂r r ∂θ ∂z (29.)
En flujo incompresible queda la ecuación diferencial:
∂u ∂v ∂w
cartesianas : + + =0
r ∂x ∂y ∂z
∇⋅v = 0 ⇒
1 ∂ (r v r ) 1 ∂ (r v θ ) ∂v z
cilíndricas : + + =0
r ∂r r ∂θ ∂z (30.)
La ecuación de continuidad para flujo incompresible es lineal y por tanto con solución. Por ello es
interesante conocer en que condiciones es aplicable; estrictamente es válida cuando la densidad del fluido es
constante; no obstante en la mayor parte de las aplicaciones de líquidos, la densidad prácticamente no varía, así
como en gases a baja velocidad.
Un criterio, para poder establecer cuando es aplicable la incompresibilidad del fluido, viene marcada
por la relación entre la velocidad del fluido y la velocidad de pequeñas perturbaciones en el seno del fluido (que
se denomina velocidad sónica); la relación es el número de Mach, y el criterio es considerar flujo incompresible
a bajos números de Mach, tomando normalmente el límite en 0,3; es decir si Ma ≤ 0,3 se puede considerar la
hipótesis de incompresibilidad.
∂ (ρv) ∂v ∂ρ ∂v δρ δv
INCOMPRESIBILIDAD ⇒ ≈ρ ≡ v << ρ ≡ <<
∂y ∂y ∂y ∂y ρ v
δp
VELOCIDAD SÓNICA: a2 ≈ ; ECUACIÓN ENERGÍA: δp ≈ −ρv δv
δρ
δρ 1 δp
2
ρ a2 ρ v2 ⎛v⎞
CRITERIO DE INCOMPRESIBLIDAD: << 1 ≡ = = ⎜ ⎟ = Ma 2 << 1 CRITERIO PRÁCTICO: Ma < 0,3
δv δp 1 a2 ⎝a⎠
v ρ v2
_________________________________________________________________________________________________________________
3.2. ECUACIÓN DIFERENCIAL DE CONSERVACIÓN DE CANTIDAD

DE MOVIMIENTO: ECUACIÓN DE MOVIMIENTO DE CAUCHY.
Considerando la primera ley del movimiento de NEWTON aplicadas a una partícula fluida en el seno de
un campo fluido o flujo, se pueden establecer el principio de conservación de cantidad de movimiento para
una partícula fluida: en una partícula en equilibrio su cantidad de movimiento se conserva; ello permite
establecer como nula la resultante de las fuerzas que actúan sobre la partícula.
Una partícula fluida es una porción de fluido de dimensiones infinitesimales y arbitrarias; el tamaño
esta en relación a las dimensiones del equipo de medida y del tiempo de medición. En todo caso los valores de
las magnitudes son medias temporales y espaciales, referidas a un intervalo de tiempo elemental y al conjunto de
moléculas que integran la partícula fluida.
Para su análisis, la partícula fluida se aísla del fluido que la envuelve mediante las superficies de
contacto partícula-fluido, y en este estado de equilibrio de la partícula aislada, se analizan las fuerzas que la
mantienen en equilibrio; estas fuerzas se dividen en tres tipos: fuerzas de volumen, fuerzas de superficie y
fuerzas de inercia.
FUERZA DE VOLUMEN: en función de que la masa de fluido (contenida en el volumen de la partícula) esta en
una determinada posición de un campo de fuerzas; lo más usual es que el campo de fuerzas sea central,
y que sea el campo gravitatorio. La evaluación de estas fuerzas es simple si derivan de un campo
central, del que se conoce su vector aceleración, y que genéricamente se denomina; en el caso de campo
r r
gravitatorio, éste vector tiene únicamente componente vertical: g = − g k . A estas fuerzas se les
denomina fuerzas másicas o fuerzas de volumen. La expresión diferencial de las fuerzas de volumen de
un campo central sobre una partícula fluida de volumen elemental dV y de masa dm es:
r r
dFv = g dm r r
⇒ dFv = ρ g dV
dm = ρ dV (31.)
FUERZAS DE SUPERFICIE: las fuerzas de contacto, que sobre las superficies de la partícula, ejerce el fluido
que la rodea, se denominan fuerzas de superficie y son debidas a los esfuerzos en las superficies de contacto
partícula fluido; los esfuerzos son debidos a la presión termodinámica y a los esfuerzos viscosos que aparecen en
el movimiento del fluido con gradiente de velocidad.
− p + τ xx τ yx τ zx −p 0 0 τ xx τ yx τ zx
tensor de tensiones: σ = τxy − p + τ yy τzy = 0 − p 0 + τ xy τ yy τ zy = − p ⋅ 1 + τ
τ xz τ yz − p + τ zz 0 0 0 τ xz τ yz τ zz
en donde “p” es la presión termodinámica y τij las tensiones viscosas
La resultante de las fuerzas de contacto sobre toda la partícula fluida viene determinada por el gradiente de
presión y por el gradiente del tensor de tensiones viscosas:
r
dFs = dV ⋅ ⎛⎜ − ∇p + ∇ ⋅ τ ⎞⎟
⎝ ⎠ (32.)
_________________________________________________________________________________________________________________
FUERZAS DE INERCIA: las fuerzas de inercia que el fluido ejerce sobre su entorno, vienen dada por su masa y
por su aceleración; y la fuerza de inercia de reacción correspondiente (3ª ley de Newton del movimiento) del
entorno del flujo sobre la partícula fluida será:
r r
dv
dFi = −ρdV ⋅
dt (33.)
Al estar la partícula en equilibrio, la resultante de las fuerzas que actúan sobre ella es nula, con lo que
combinando las ecuaciones anteriores se tiene:
r
dFvolumen = ρ dV g
⎛ r r =⎞
dFvolumen + dFsup erficie + dFinercia = 0 dFsup erficie = ⎜⎜ − ∇p + ∇ ⋅ τ ⎟⎟ dV
≡
⎝ ⎠
r
dv
dFinercia = − ρ dV
dt
r r
r ⎛ ⎞ ⎡ r ⎛ ⎞ dv ⎤
∑ dF = ρ dV g + ⎜⎝ − ∇p + ∇ ⋅ τ ⎟⎠dV − ρdV dt = dV ⎢⎣ρ g + ⎜⎝ − ∇p + ∇ ⋅ τ ⎟⎠ − ρ dt ⎥⎦ = 0
dv
⇒
∑ dF = ρ gr + ⎛⎜ − ∇p + ∇ ⋅ τ ⎞⎟ − ρ dvr = 0
dV ⎝ ⎠ dt (34.)
Con lo que para una partícula fluida, la expresión de la 1ª ley de Newton del movimiento o ecuación de
conservación de cantidad de movimiento, y que en Mecánica de Fluidos se denomina ecuación de movimiento
de CAUCHY, es:
r
r = dv
ρ g −∇p + ∇ ⋅ τ = ρ
dt
(35.)
En coordenadas cartesianas los términos de la ecuación de movimiento de CAUCHY son:
gradiente de presión: ∇p = i+ j+ k
∂x ∂y ∂z
divergencia del tensor de tensiones viscosas:
⎛ ∂τ ∂τ ∂τ ⎞ r ⎛ ∂τ ∂τ ∂τ ⎞ r ⎛ ∂τ ∂τ ∂τ ⎞ r
∇ ⋅ τ = ⎜ xx + yx + zx ⎟ ⋅ i + ⎜ xy + yy + zy ⎟ ⋅ j + ⎜ xz + yz + zz ⎟ ⋅ k
⎝ ∂x ∂y ∂z ⎠ ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂x ∂y ∂z ⎠
r r
dv ∂v r r
vector aceleración ( ) = aceleración local ( ) + aceleración convectiva ( v ⋅ ( v ) )
dt ∂t
r r r r r
dv ∂v ∂v ∂v ∂v ⎛ ∂u ∂u ∂u ∂u ⎞ r ⎛ ∂v ∂v ∂v ∂v ⎞ r ⎛ ∂w ∂w ∂w ∂w ⎞ r
= +u +v +w = ⎜⎜ +u +v + w ⎟⎟ ⋅ i + ⎜⎜ +u +v + w ⎟⎟ ⋅ j + ⎜⎜ +u +v +w ⎟⋅k
dt ∂t ∂x ∂y ∂z ⎝ ∂t ∂x ∂y ∂z ⎠ ⎝ ∂t ∂x ∂y ∂z ⎠ ⎝ ∂t ∂x ∂y ∂z ⎟⎠
_________________________________________________________________________________________________________________
La ecuación vectorial de movimiento de CAUCHY, da lugar a tres ecuaciones diferenciales, una por cada
componente; teniendo en coordenadas cartesianas las siguientes ecuaciones de movimiento de CAUCHY:
∂p ⎛⎜ ∂τ xx ∂τ yx ∂τ zx ⎞ ⎛ ⎞
⎟ + ρg x = ρ⎜ ∂u + u ∂u + v ∂u + w ∂u ⎟
− + + + ⎜ ∂t ⎟
∂x ⎜⎝ ∂x ∂y ∂z ⎟
⎠ ⎝ ∂x ∂y ∂z ⎠
∂p ⎛⎜ ∂τ xy ∂τ yy ∂τ zy ⎞ ⎛
⎟ + ρg y = ρ⎜ ∂v + u ∂v + v ∂v + w ∂v ⎟
⎞
− + + + ⎜ ∂t
∂y ⎜⎝ ∂x ∂y ∂z ⎟
⎠ ⎝ ∂x ∂y ∂z ⎟⎠
∂p ⎛⎜ ∂τ xz ∂τ yz ∂τ zz ⎞ ⎛
⎟ + ρg z = ρ⎜ ∂w + u ∂w + v ∂w + w ∂w ⎟
⎞
− + + + ⎜ ∂t
∂z ⎜⎝ ∂x ∂y ∂z ⎟
⎠ ⎝ ∂x ∂y ∂z ⎟⎠ (36.)
3.2.1. FLUIDO NO VISCOSO: ECUACIÓN DE EULER.
Para poder utilizar la ecuación de movimiento de CAUCHY, es necesario conocer los términos de las
tensiones viscosas; que están relacionadas con la velocidad de deformación. La relación entre tensiones y
velocidades de deformación depende de la propia constitución del fluido correspondiente.
El caso más simple es cuando el fluido es no viscoso, y son idénticamente nulas todas las tensiones
viscosas, con lo que la ecuación de movimiento de CAUCHY se simplifica, quedando la ecuación de EULER del
movimiento de un fluido no viscoso:
r
r dv
− ∇ p + ρg = ρ
dt (37.)
En la ecuación de EULER, es interesante introducir el término de vorticidad local del fluido, que forma
parte del vector aceleración:
r r
dv ∂v r r
vector aceleración: = + v ⋅ ( ∇v )
dt ∂t
r r ⎛ vr 2 ⎞ r r
aceleración convectiva: v ⋅ ( ∇v ) = ∇ ⎜ ⎟+ Ω x v
⎜ 2 ⎟
⎝ ⎠
r r r
vorticidad: Ω = 2ω = ∇xv
r 1 ⎛ ∂w ∂v ⎞ r 1 ⎛ ∂u ∂w ⎞ r 1 ⎛ ∂v ∂u ⎞ r
vector velocidad angular. ω = ⎜ − ⎟⋅ i + ⎜ − ⎟⋅ j + ⎜ − ⎟⋅k
2 ⎝ ∂y ∂z ⎠ 2 ⎝ ∂z ∂x ⎠ 2 ⎝ ∂x ∂y ⎠
Con todo lo anterior, la ecuación de EULER queda como:

r
( )
∂v 1 r r 2 ∇p r
+ ∇⋅ v +
∂t 2 ρ
r r
− g + Ω x v =0 ( ) (38.)
La dificultad de trabajar con flujos con vorticidad, es inherente a la dificultad de la ecuación anterior; no
obstante con determinadas condiciones el término que incluye la vorticidad es nulo, aún no siendo nula la
verticidad.
_________________________________________________________________________________________________________________
Para poder establecer las condiciones en donde el término de vorticidad es nulo, multipliquemos
r
escalarmente la ecuación anterior por un vector desplazamiento “ dr ”, arbitrario; con lo queda la ecuación:
r
( )
⎡ ∂v 1 r 2 ∇p r r ⎤ r
r
(
⎢ ∂t + 2 ∇ ⋅ v + ρ − g + Ω x v ⎥ ⋅ dr = 0 )
⎣ ⎦ (39.)
En donde el término que incluye la vorticidad, es nulo cuando:
r r r
( Ω x v ) ⋅ dr = 0
(40.)
Lo que tiene lugar bajo una de las siguientes condiciones:

r r
- v = 0 : no hay flujo, estamos en fluidoestática y la ecuación de EULER queda: ∇p = ρg
r
( )
r ∂v 1 r2 r
- Ω = 0 : flujo irrotacional, la ecuación de Euler queda: ρ + ρ∇ ⋅ v + ∇p − ρg = 0
∂t 2
r r r
- dr es perpendicular al vector Ω x v : se denomina flujo de BELTRAMI
r r
- d r es paralelo al vector velocidad v : la ecuación de Euler queda como:
r
∂v r 1
∂t 2
( )
r 2 r ∇p r
⋅ dr + ∇ v ⋅ dr +
ρ
r r
⋅ dr − g ⋅ dr = 0
(41.)
Esta última condición es la que conduce a resultados más útiles, y se tiene en los puntos del flujo que
cumplen la condición de que su vector desplazamiento es paralelo al vector velocidad, es decir puntos de las
r r r r
líneas de corriente, en donde v = v ⋅ u LC y dr = ds ⋅ u LC ; en donde “ds” es el módulo del vector desplazamiento a
lo largo de una línea de corriente
r y “v” el módulo del vector velocidad. Considerando además, exclusivamente
r
campo gravitatorio ( g = −gk ), los términos de la ecuación anterior son:
r r
∂v r ∂ ( v ⋅ u LC ) r ∂v
⋅ dr = ⋅ ( ds ⋅ u LC ) = ⋅ ds
∂t ∂t ∂t
1
2
( )
r2 r 1
2
r 1 ⎛ dv 2 r ⎞
∇ v ⋅ dr = ∇ ( v 2 ) ⋅ ( ds ⋅ u LC ) = ⎜
2 ⎝ ds ⎠
r 1
⋅ u LC ⎟ ⋅ ( ds ⋅ u LC ) = dv 2
2
dp r
⋅u
∇p r ds LC r dp
⋅ dr = ⋅ ( ds ⋅ u LC ) =
ρ ρ ρ
r r r r
r r
( )( )
− g ⋅ dr = − − gk ⋅ dx ⋅ i + dy ⋅ j + dz ⋅ k = g ⋅ dz
Obteniendo la Ec. de EULER para puntos de una línea de corriente:
∂V 1 dp
⋅ ds + dv 2 + + g ⋅ dz = 0
∂t 2 ρ (42.)
_________________________________________________________________________________________________________________
3.2.2. FLUJO NO VISCOSO EN PUNTOS DE LÍNEA DE CORRIENTE: EC. DE BERNOULLI.
En la ecuación de Euler en puntos de una línea de corriente, la única condición restrictiva es considerar
flujo no viscoso. La ecuación es de gran utilidad pues es posible integrarla entre dos puntos de una misma línea
de corriente, pues dos términos son diferenciales exactas; obteniéndose la ecuación de BERNOULLI para flujo no
estacionario y no viscoso:
∂v v 2 − v12
∫ ∫
2 2
dp
⋅ ds + 2 + + g(z 2 − z1 ) = 0
1 ∂t 2 1 ρ (43.)
La ecuación de BERNOULLI, para flujo estacionario e incompresible (además de no viscoso), toma la

forma:
p2 − p1 v 22 − v12
+ + g(z 2 − z1 ) = 0
ρ 2 (44.)
Es decir, a lo largo de una línea de corriente, permanece constante la suma de los tres términos:
p v2
+ + gz = cte.
ρ 2 (45.)
La constante de la ecuación de Bernoulli puede variar de una línea de corriente a otra, a no ser que
además se tenga la condición de irrotacionalidad, con la que la ecuación se cumple independientemente de la
r
dirección del vector desplazamiento (d r ), y la constante de la ecuación de Bernoulli es la misma en todo el
flujo.
En la ecuación anterior, todos los sumandos, son dimensionalmente términos de energía específica
(energía por unidad de masa), representando cada término: el trabajo de flujo, la energía cinética y la energía
potencial, en un determinado punto de la línea de corriente.
La Ec. de BERNOULLI, también se suele expresar en términos de presión:
1
p + ρv 2 + ρgz = cte.
2 (46.)
representando cada término:
p = presión absoluta o termodinámica

1 2
ρv = presión dinámica
2
ρgz = presión hidrostática
1
p + ρv 2 = presión de estancamiento
2
p + ρgz = presión piezométrica
_________________________________________________________________________________________________________________
3.2.3. FLUIDO NEWTONIANO: ECUACIONES DE NAVIER-STOKES.
Para un fluido newtoniano las tensiones viscosas son proporcionales a las velocidades de deformación,
r
y viene determinadas por la ecuación de NAVIER-POISSON (Ec.18.): τ = 2µ ⋅ ε& + λ ( ∇ ⋅ v ) ⋅ 1 . Con lo que en la
ecuación de movimiento de CAUCHY, la divergencia del tensor de tensiones viscosas, para el caso de un fluido
newtoniano, se puede obtener a partir de la expresión de su tensor de tensiones, quedando:
[ ( )]
∇ ⋅ τ = (λ + µ ) ∇ ∇ ⋅ v + µ∇ 2 v
(47.)
Con lo que la ecuación de movimiento para un fluido newtoniano, que se denomina ecuación de
NAVIER-STOKES queda como:
[ ( )]
ρg − ∇p + (λ + µ ) ∇ ∇ ⋅ v + µ∇ 2 v = ρ
dv
dt (48.)
Con la condición de STOKES: λ = -2µ/3, es decir (λ+µ)=µ/3, queda como:
ρg − ∇p +
µ
3
[ ( )]
∇ ∇ ⋅ v + µ∇ 2 v = ρ
dv
dt (49.)
r
En donde ∇ 2 v es el vector “laplaciana de velocidad”, y que puede descomponerse en función de la
r r r
divergencia y de la vorticidad: ∇ 2 v = ∇ ( ∇ ⋅ v ) − ∇xΩ ; con lo que la ecuación vectorial de NAVIER-STOKES
para el movimiento para un fluido newtoniano se puede rescribir como:
r r
r 4 r dv
ρ ⋅ g − ∇p + µ ⋅ ∇ (∇ ⋅ v ) − µ ⋅ ∇xΩ + = ρ ⋅
3 dt (50.)
La ecuación de NAVIER-STOKES, con la restricción de fluido no viscoso: µ=0, lleva evidentemente a la

r
r dv
ecuación de EULER: −∇p + ρ ⋅ g = ρ ⋅
dt
r
Las restricciones simultaneas de flujo incompresible (adivergente: ∇ ⋅ v = 0) y de flujo irrotacional
r
(vorticidad nula: ∇xv = 0), también llevan a la ecuación de EULER, aunque el fluido sea viscoso.
r r
En la ecuación de NAVIER-STOKES, con la restricción de flujo incompresible (adivergente: ∇ ⋅ v = 0) ,
el segundo término es nulo por obligar la incompresibilidad, a que el flujo sea adivergente (divergencia del
vector velocidad nula); con lo que se tiene la ecuación de NAVIER-STOKES para el flujo de un fluido
newtoniano incompresible:
r
r r dv
ρ ⋅ g − ∇p + µ ⋅ ∇ 2 v + = ρ ⋅
dt (51.)
_________________________________________________________________________________________________________________
La ecuación vectorial anterior, puede descomponerse en sus tres componentes cartesianas, teniendo las
siguientes ecuaciones escalares:
⎛ ∂ 2u ∂ 2u ∂ 2u ⎞
−
∂p
+ µ⎜ + + ⎟ + ρg = ρ⎛⎜ ∂u + u ∂u + v ∂u + w ∂u ⎞⎟
⎜ ⎜ ∂t ∂z ⎟⎠
∂z ⎟⎠
x
∂x 2 2 2 ⎝ ∂x ∂y
⎝ ∂x ∂y
⎛ ∂2v ∂2v ∂2v ⎞

−
∂p
+ µ⎜ + + ⎟ + ρg = ρ⎛⎜ ∂v + u ∂v + v ∂v + w ∂v ⎞⎟
⎜ y ⎜ ∂t ∂z ⎟⎠
∂y
⎝ ∂x
2
∂y 2
∂z ⎟⎠
2 ⎝ ∂x ∂y
⎛ ∂ 2w ∂2w ∂2w ⎞
−
∂p
+ µ⎜ + + ⎟ + ρg = ρ⎛⎜ ∂w + u ∂w + v ∂w + w ∂w ⎞⎟
⎜ ⎜ ∂t ∂z ⎟⎠
∂z ⎟⎠
z
∂z 2 2 2 ⎝ ∂x ∂y
⎝ ∂x ∂y (52.)
Estas son las que normalmente se denominan ecuaciones de Navier-Stokes, y fueron obtenidas de forma
independiente por Louis M. NAVIER (1785-1836) y por George G. STOKES (1819-1903). Son ecuaciones
diferenciales en derivadas parciales de 2º orden no lineales. Las tres ecuaciones de Navier-Stokes, incluyen 4
incógnitas: la presión (p) y las tres componentes de la velocidad (u,v,w). La ecuación que completa el sistema
homogéneo de ecuaciones diferenciales, es la ecuación de continuidad (du/dx+dv/dy+dw/dz=0). Esto en el caso
de flujo incompresible, en donde la densidad es constante; en el caso de flujo compresible, aparecen tres nuevas
incógnitas, la densidad, la temperatura y la energía interna, necesitando para completar el sistema homogéneo, 3
nuevas ecuaciones, que como veremos posteriormente son la ecuación de energía, la ecuación térmica de estado
y la ecuación calórica de estado.
3.2.4. LA FUNCIÓN DE CORRIENTE Y LA FUNCIÓN POTENCIAL.
LA FUNCIÓN DE CORRIENTE
Si consideremos un flujo estacionario incompresible bidimensional, la ecuación de continuidad se

reduce a:
∂u ∂v
+ =0
∂x ∂y
Con lo que se puede introducir una función ψ, definida por:
∂ψ ∂ψ r ∂ψ r ∂ψ r
u= ∧ v=− ≡ v= i− j
∂y ∂x ∂y ∂x
r r ⎛ ∂ 2ψ ∂ 2ψ ⎞ r r
Con lo que la expresión de la vorticidad es: Ω = ∇xv = ... = − ⎜ 2 + 2 ⎟ ⋅ k = −∇ 2ψ ⋅ k
⎝ ∂x ∂y ⎠
De la ecuación de movimiento se puede establecer:
∂ψ ∂ (∇ 2 ψ ) ∂ψ ∂ (∇ 2 ψ ) µ 2 2
− = ⋅ ∇ (∇ ψ )
∂y ∂x ∂x ∂y ρ (53.)
que es una ecuación diferencial de 4º orden, con necesidad de solución por análisis numérico.
_________________________________________________________________________________________________________________
Si se añade además la restricción de flujo irrotacional, la ecuación anterior se reduce a la ecuación de

LAPLACE en dos dimensiones:
∂ 2ψ ∂ 2ψ
+ = ∇ 2ψ = 0
2 2
∂x ∂y
Se puede obtener que en una línea de corriente no hay cambio en la función ψ, por lo que a la citada
función se le denomina función de corriente:
dx dy
ecuación línea de corriente: = ≡ − v ⋅ dx + u ⋅ dy = 0
u v
∂ψ ∂ψ
introduciendo la función de corriente: dx + dy = 0 = dψ ⇔ ψ = cte.
∂x ∂y
LA FUNCIÓN POTENCIAL
Consideremos como única restricción que el flujo es irrotacional, con ello se tiene que la vorticidad es
nula y se obtiene que el vector velocidad es el gradiente de una función escalar5, a la que se denomina función
potencial de velocidad:
r r r
Ω = ∇xv = 0 ⇔ v = ∇φ
Conocida la función potencial de velocidad φ=φ(x,y,z,t), se obtienen fácilmente las componentes del
vector velocidad:
∂φ ∂φ ∂φ
u= v= w=
∂x ∂y ∂z
El lugar geométrico de puntos del flujo con igual función potencial, se denomina superficie
equipotencial. En el caso particular de flujo bidimensional, el lugar geométrico esta contenido en el plano en el
que se fluye el fluido, son las líneas equipotenciales, y al ser bidimensional existe también la función de
corriente, verificándose que las líneas de corriente son ortogonales a las líneas equipotenciales.
Si además de irrotacional, se considera flujo incompresible, se obtiene que la función potencial

cumple la ecuación de LAPLACE en tres dimensiones:
r r ∇2 φ = 0
∇xv = 0 ⇒ v = ∇φ
r ⇒ ∇ ⋅ (∇φ) = 0 ≡ ∂ 2 φ ∂ 2 φ ∂ 2 φ
∇⋅v = 0 + + =0
∂x 2 ∂y 2 ∂z 2
5
Por cálculo vectorial, un vector con rotacional nulo es el gradiente de una función escalar.
_________________________________________________________________________________________________________________
3.3 ECUACIÓN DIFERENCIAL DE CONSERVACIÓN DE ENERGÍA: ECUACIÓN DE ENERGÍA.
El principio de conservación de energía (PRIMER PRINCIPIO DE TERMODINÁMICA) aplicado a

una partícula fluida, establece que la energía total de la partícula fluida es constante, siempre que no existan
aportes energéticos por transferencia de calor o de trabajo.
Siguiendo el criterio termodinámico de signos, se consideran como positivos el trabajo desarrollado por
la partícula y el calor aportado a la partícula, y como negativos el trabajo consumido por la partícula y el calor
cedido por la partícula; con todo la ecuación de conservación de energía es:
dE d(uˆ + v 2 / 2 + gz)
= ρdV
dt dt
δQ &
dE − δQ + δW = 0 ≡ = Q = dV [ ∇ ⋅ ( κ ∇T)]
dt
δW & ⎡ ⎤
= W = dV ⎢ −∇ ⋅ ⎛⎜ v ⋅ τ ⎞⎟ + v ⋅ ∇p ⎥
r = r
dt ⎣ ⎝ ⎠ ⎦
d(uˆ + v 2 / 2 + gz) r r =
ρ + v ⋅ ∇p = ∇ ⋅ ( κ∇T) + ∇ ⋅ ( v ⋅ τ )
dt (54.)
En donde:
- la energía total de la partícula viene dada por la suma de la energía interna, la energía cinética y la
energía potencial: E = U + Ec + Ep = mû + mv 2 / 2 + mgz
- la transferencia de calor (por unidad de tiempo) entre partícula y su entorno por conducción viene
determinada por el gradiente de temperatura ( ∇T ) y por la conductividad térmica (κ)
- el trabajo (por unidad de tiempo) intercambiado entre partícula y su entorno tiene dos términos, el
debido a las fuerzas de presión (trabajo de flujo) y el debido a los esfuerzos viscosos.
El trabajo debido a los esfuerzos viscosos, se puede expresar como suma de dos términos, introduciendo
el concepto de función de disipación viscosa de RAYLEIG Φ:
r r =
&
W viscosidad = ∇ ⋅ ( v ⋅ τ ) = v ⋅ ( ∇ ⋅ τ ) + Φ
(55.)
r = r =
Φ = ∇ ⋅ ( v ⋅ τ) − v ⋅ (∇ ⋅ τ)
(56.)
En coordenadas cartesianas para un fluido newtoniano, la función de disipación viscosa es:
⎡ 2
⎛ ∂v ⎞
2 2
⎛ ∂v ∂u ⎞
2
⎛ ∂w ∂v ⎞
2 2⎤
⎛ ∂u ⎞ ⎛ ∂w ⎞ ⎛ ∂u ∂w ⎞ ⎥
Φ = µ ⎢2⎜ ⎟ + 2⎜⎜ ⎟⎟ + 2⎜ ⎟ + ⎜
⎜ + ⎟
⎟ + ⎜
⎜ + ⎟
⎟ +⎜ + ⎟
⎢ ⎝ ∂x ⎠ ⎝ ∂y ⎠ ⎝ ∂z ⎠ ⎝ ∂x ∂y ⎠ ⎝ ∂y ∂z ⎠ ⎝ ∂z ∂x ⎠ ⎥
⎣ ⎦
En la ecuación de disipación viscosa todos los términos son cuadráticos, por lo que su valor siempre es
positivo, es decir en flujo viscoso parte de su energía disponible se disipa por las irreversibilidades de los
fenómenos de transporte de cantidad de movimiento entre partículas; lo que esta de acuerdo con el segundo
principio de Termodinámica de que los procesos reales son irreversibles con degradación de energía y su
consiguiente aumento de entropía del universo.
_________________________________________________________________________________________________________________
3.3.1. ECUACIÓN DE ENERGÍA INTERNA.
Introduciendo el término de función de disipación viscosa en la ecuación de energía (29), y utilizando la

ecuación de Navier-Stokes para un fluido newtoniano (23) multiplicada escalarmente por el vector velocidad,
r
para que desaparezca el término v ⋅ (∇ ⋅ τ) , se obtiene una expresión de la ecuación de energía en donde no
aparecen las energías cinética ni potencia, solo la energía interna:
dû r
ρ = − p(∇ ⋅ v ) + κ ⋅ ∇ 2 T + Φ
dt
(57.)
La ecuación anterior permite obtener la energía interna6 en función del flujo de calor por conducción, del
trabajo de expansión o compresión y de la disipación por viscosidad. Esta ecuación de energía, es valida para un
fluido newtoniano en condiciones muy generales de flujo transitorio, compresible, viscoso y conductor de calor;
solo se desprecian la transferencia de calor por radiación y por fuentes internas.
En la ecuación de energía anterior, la derivada total de la energía interna se puede expresar como suma de
duˆ ∂uˆ r
dos términos: el de variación local y el de variación convectiva: = + v ⋅ ∇û .
dt ∂t
r r ∂û ∂û ∂û

El termino convectivo v ⋅ ∇û = u +v +w , representa el transporte de calor por convección.
∂x ∂y ∂z
Como primera simplificación restrictiva para el manejo de la ecuación de energía, se suele considerar que
la energía interna es proporcional a la temperatura, a través del calor específico a volumen constante:
⎛ ∂û ⎞
cv = ⎜ ⎟
⎝ ∂T ⎠ v
considerando además la constancia de los siguientes parámetros del fluido: calor especifico a volumen constante,
viscosidad dinámica, coeficiente de conductividad térmica y densidad: cv, µ,κ,ρ , se tiene la forma más simple de
la ecuación de la energía:
dT r
ρ ⋅ cv + p(∇ ⋅ v ) = κ ⋅ ∇ 2 T + Φ
dt (58.)
En el caso particular de flujos muy lentos o en reposo, se pueden despreciar los términos disipativos y
convectivos, con lo que se tiene la siguiente expresión que permite obtener el campo de temperaturas; que es la
ecuación de conducción de calor para sólidos y fluidos en reposo:
∂T
ρ ⋅ cv = κ ⋅ ∇ 2T
∂t (59.)
6
La energía interna esta asociada a la agitación térmica y a la composición química; en el caso de un gas ideal no hay interrelaciones entre
los átomos, y la energía interna solo depende de la temperatura termodinámica.
_________________________________________________________________________________________________________________
3.3.2. ECUACIÓN DE ENTALPÍA.
La entalpía es la suma de la energía interna y el trabajo de flujo: h = u + p / ρ ; con lo que la variación

temporal de la entalpía es:
dh dû 1 dp p dρ
= + −
dt dt ρ dt ρ 2 dt
el término de variación temporal de la densidad, se puede expresar en función de la divergencia de la velocidad ,

a partir de la ecuación de continuidad:
r 1 dρ
∇⋅v = −
ρ dt
con lo que se tiene:
dh duˆ 1 dp p r
= + + (∇ ⋅ v )
dt dt ρ dt ρ
y finalmente con la ecuación de energía interna (48), se llega a la ecuación de entalpía:
dh dp
ρ = + κ ⋅ ∇ 2T + Φ
dt dt
(60.)
3.3.3. ECUACIÓN DE ENTROPÍA.
El Segundo Principio de Termodinámica se puede establecer como:
ds 1 duˆ p r
= + ⋅ (∇ ⋅ v )
dt T dt ρT
(61.)
Con lo que utilizando la ecuación de energía interna (48), se llega a la ecuación de la entropía:
ds
ρT = κ ⋅ ∇ 2T + Φ
dt
(62.)
El termino de disipación de energía es siempre positivo, con lo que genera siempre aumento de entropía:
es lo inherente al Segundo Principio de Termodinámica: las irreversibilidades hacen aumentar la entropía; el
término de transmisión de calor por conducción, aumenta la entropía si el flujo de calor es positivo (es decir se
calienta el flujo) y disminuye la entropía si en flujo de calor es negativo ( es decir se enfría el flujo)
_________________________________________________________________________________________________________________
3.4. CONDICIONES DE CONTORNO.
A partir de las ecuaciones de conservación para una partícula fluida se han obtenido las ecuaciones:
∂ρ r
(1) ecuación de continuidad: + ∇ ⋅ (ρv ) = 0
∂t
r
r r r dv
(2)(3)(4) ecuación de movimiento: ρ ⋅ g − ∇p + (λ + µ) ⋅ ( ∇ (∇ ⋅ v ) ) + µ ⋅ ∇ 2 v = ρ ⋅
dt
dû r
(5) ecuación de energía: ρ + p( ∇ ⋅ v ) = ∇ ⋅ ( κ ⋅ ∇ T) + Φ
dt
Las ecuaciones de continuidad y de energía son ecuaciones diferenciales escalares y la ecuación de

movimiento es vectorial, por lo que entre todas aportan 5 ecuaciones diferenciales escalares. En cuanto a las
incógnitas se tienen: la densidad (ρ), las componentes del vector velocidad (u,v,w), la presión (p), la temperatura
(T) y la energía interna (û), es decir se tienen 7 incógnitas, por lo que para poder tener un sistema homogéneo de
ecuaciones es necesario disponer de 2 ecuaciones adicionales; estas ecuaciones son las ecuaciones de estado de
constitución del propio fluido considerado:
(6) ecuación térmica de estado: ρ = ρ(p,T)

(7) ecuación calórica de estado: û = û(p,T)
Con todo lo expuesto anteriormente, se dispone de un sistema homogéneo de 7 ecuaciones diferenciales

con 7 incógnitas, cuya resolución es posible, con las condiciones de contorno apropiadas para cada caso, y
normalmente con técnicas numéricas, siendo posible solo para casos muy concretos la solución analítica.
Con la restricción de flujo incompresible y propiedades constantes, se tiene solo 5 incógnitas: la

presión, las tres componentes de la velocidad y la temperatura; siendo suficientes las ecuaciones de continuidad,
movimiento (3) y energía:
r
(1) continuidad: ∇⋅v = 0
r
r r dv
(2)(3)(4) movimiento: ρ ⋅ g − ∇p +µ ⋅ ∇ 2 v = ρ ⋅
dt
dT
(5) energía: ρ ⋅ cv = κ ⋅ ∇ 2T + Φ
dt
además la ecuación de energía esta desacoplada, es decir en las cuatro ecuaciones aportadas por la continuidad y
por la cantidad de movimiento, sólo aparecen 4 incógnitas: presión y componentes de la velocidad, por lo que es
posible su resolución; si se requiere el campo de temperaturas, se obtiene a partir de la ecuación de energía,
previo conocimiento del campo de velocidades.
_________________________________________________________________________________________________________________
La solución de los sistemas de ecuaciones diferenciales anteriores, están condicionadas por las
condiciones de contorno apropiadas, que dependen de cada caso, y vienen determinadas por los valores de las
propiedades en el instante inicial, por la geometría de las paredes y por las condiciones en las entradas y en las
salidas.
En las paredes impuestas por la geometría en la que esta confinado el flujo, se tiene la condición de no
deslizamiento ni de cambio de temperatura, es decir: en las partículas que “tocan” una pared se ponen a la
velocidad de la pared y a su temperatura: velocidad del fluido en la pared = velocidad de la pared; y temperatura
del fluido en la pared = temperatura de la pared.
Un caso muy particular de condición de contorno impuesta por una pared, es el caso de los flujos que se
consideran no viscosos, en donde no se cumple la condición de no deslizamiento; siendo la única condición de
contorno establecida por la pared, que el flujo no la atraviese, es decir que sean iguales las velocidades normales
de la pared y del fluido, no pudiendo decir nada sobre la velocidad tangencial del flujo cerca de la pared.
En las entradas y salidas se deben conocer las distribuciones de velocidad, temperatura y presión.
Las condiciones de contorno más complejas se tienen cuando existe superficie libre, en la interfase
líquido-líquido o líquido-gas; en donde se cumple la condición cinemática de contorno, de igualdad de
velocidades perpendiculares a la superficie de separación (no debe haber huecos entre el líquido y el gas); así
como el equilibrio de tensiones en la superficie libre (excepto por los efectos de tensión superficial), es decir
igualdad de tensión normal o presión e igualdad de tensión tangencial. Además debe cumplirse la condición de
igualdad de temperaturas en todos los puntos de la superficie libre.
Por la distinta viscosidad de cada fluido, son distintos los gradientes de velocidad de cada fluido en la
superficie libre, aunque los esfuerzos tangenciales deben ser iguales, con lo que el perfil de velocidades (que
incluye la propia superficie libre) si que debe ser una función continua, pero no es derivable en los puntos de la
superficie libre:
⎛ ∂u ⎞ τL
⎜⎜ ⎟⎟
⎝ ∂y ⎠ L µL µ
= = G
⎛ ∂u ⎞ τG µL
⎜⎜ ⎟⎟ µG
⎝ ∂y ⎠ G
Normalmente el aumento de presión debido al efecto de la tensión superficial es despreciable, excepto

⎛ 1 1 ⎞
cuando los radios de curvatura son muy pequeños: ∆p = σ ⎜⎜ + ⎟ ; así en el caso de una gota de
⎟
⎝ RL RG ⎠
líquido en el seno de un gas o de otro líquido, como los radios son pequeños y además iguales, se tiene que la
2σ
sobrepresión que tiene lugar entre puntos separados por la superficie libre es: ∆p =
R gota
_________________________________________________________________________________________________________________
4. PROBLEMAS RESUELTOS.
4.1. Métodos de análisis: Euleriano y Lagrangiano
4.2. Aplicación de la ecuación de continuidad: criterios de incompresibilidad.
4.3. Aplicación de las ecuaciones de continuidad y de Bernoulli: descarga de depósitos.
4.4. Aplicación de la ecuación de Bernouilli: flujo no viscoso entre discos horizontales.
4.5. Aplicación de la ecuación de Bernouilli no estacionario: oscilaciones en un tubo en U.
4.6. Aplicación de la ecuación de Bernouilli con aceleración de arrastre: bomba rotativa.
4.7. Aplicación de las ecuaciones de Navier-Stokes: flujo de Couette-Poiseuille.
4.8. Aplicación de la ecuación de Navier-Stokes: flujo viscoso entre discos horizontales.
4.9. Aplicación de la ecuación de Energía: distribución de temperaturas en flujo de Poiseuille.
_________________________________________________________________________________________________________________
4.1. Métodos de análisis: Euleriano y Lagrangiano.

Para determinar la aceleración de una partícula, se puede utilizar el método Lagrangiano, en donde la
aceleración de la partícula se obtiene por la derivada segunda de su vector de posición, respecto al tiempo.
Utilizando el método Euleriano la aceleración de una partícula que se mueve en un campo de velocidad, es una
función del tiempo y de la posición, y es suma de la aceleración local y de la convectiva. Se considera un flujo
unidimensional, estacionario e incompresible a través de una tobera convergente. A partir de los datos:
DETERMINE: (1) Aceleración por el método Euleriano.

(2) Aceleración por el método Lagrangiano.
r ⎛ x⎞ r
DATOS: campo de velocidades: v = u0 ⎜1 + ⎟ ⋅ i
⎝ A⎠
Q “x”
x
A
RESOLUCIÓN:
r r r
En el método Euleriano, las partículas se mueven por un campo de velocidad: v = v(r , t) , en una
determinada posición y en un instante de tiempo, la aceleración de la partícula, que en ese instante, esta en la
posición determinada, es:
r r
r dv ∂v r r
a= = + v ⋅ ∇v
dt ∂t
r
∂v
es la aceleración local, y viene determinada, en un determinado punto (local), por la variación de
∂t
la velocidad con el tiempo; si el flujo es estacionario, en un determinado punto las propiedades no varían con el
tiempo (no hay variaciones locales), y en particular la velocidad en ese punto es la misma a lo largo del tiempo,
con lo que la aceleración local será nula.
r r
v ⋅ ∇v es la aceleración convectiva, y viene determinada, en un determinado instante, por el gradiente
r
de velocidad ( )
En el método Lagrangiano, se parte del conocimiento del vector de posición de una determinada
r r
partícula a lo largo del tiempo: rp = rp (t) ; y se obtiene su aceleración por la derivada segunda del vector de
r
r d 2r
posición con respecto al tiempo: a=
dt 2
(1) ACELERACIÓN DE LAS PARTICULAS EN EL MÉTODO EULERINO.

En el problema, el flujo es estacionario y unidimensional:
r
∂v
Estacionario ⇒ =0
∂t
r ⎛ x⎞ r
Unidimensional ⇒ v = u0 ⎜1 + ⎟ ⋅ i [1]
⎝ A⎠
_________________________________________________________________________________________________________________
r r r r
La aceleración a es puramente convectiva : a = v ⋅ v
⎛ ∂ ⎞ ⎛ ∂u ∂v ∂w ⎞
⎜ ⎟ ⎜ ∂x ⎟⎟
⎜ ∂x ⎟ ⎜ ∂x ∂x
∂ ∂u ∂v ∂w ⎟
v = ⎜⎜ ⎟⎟ ( u v w ) = ⎜⎜
r
El gradiente de velocidad es:
∂y ∂y ∂y ∂y ⎟
⎜ ⎟ ⎜ ⎟
⎜ ∂ ⎟ ⎜ ∂u ∂v ∂w ⎟
⎜ ⎟ ⎜ ⎟
⎝ ∂z ⎠ ⎝ ∂z ∂z ∂z ⎠
en el problema: u=u0(1+x/A), v = w = 0;
∂u u 0
= ;
∂x A
⎛ u0 ⎞
⎜A 0 0⎟
r ⎜ ⎟
con lo que el gradiente de velocidad es: v=⎜ 0 0 0⎟
⎜ 0 0 0⎟
⎜ ⎟
⎝ ⎠
⎛ u0 ⎞
⎜A 0 0⎟
r r r ⎛ ⎛ x⎞ ⎞⎜ ⎟ u2 ⎛ x⎞ r
y la aceleración convectiva es: a = v ⋅ v = ⎜ u 0 ⎜ 1 + ⎟ 0 0 ⎟ ⎜ 0 0 0 ⎟ = 0 ⎜1 + ⎟ ⋅ i
⎝ ⎝ A⎠ ⎠⎜ A ⎝ A⎠
0 0 0⎟
⎜ ⎟
⎝ ⎠
r u2 ⎛ x⎞ r
a = 0 ⎜1 + ⎟ ⋅ i [2 ]
A ⎝ A⎠
(2) ACELERACIÓN DE LAS PARTÍCULAS EN EL MÉTODO LAGRANGIANO.

Consideremos una partícula que en el instante inicial (t=0), esta situada en el inicio de la tobera (x=0); la
posición de esa partícula a lo largo del tiempo es: x p = x p ( t ), y se determina a partir del campo de
velocidades:
⎛ x ⎞
∫ ∫
dx p xp dx p t
dx p = u ⋅ dt = u 0 ⎜ 1 + p ⎟ ⋅ dt ⇒ = u 0 ⋅ dt ⇒ = u 0 ⋅ dt
⎝ A⎠ xp xp =o 1 + xp / A t =0
1+
A
⎛ ⎛
⎜ u 0 t ⎞⎟ ⎞
⎛ xp ⎞ ⎜ ⎜
A ⎟⎠ − 1⎟
⎟
A ⋅ ln ⎜ 1 +
⎜ ⎟ = u0 ⋅ t
A ⎟⎠
⇒ xp = A ⎜ e ⎜
⎝
⎟
⎝ ⎜ ⎟
⎜ ⎟
⎝ ⎠
r d2x r u 2 u0 t r
con lo que la aceleración de la partícula será: a p = 2p ⋅ i = ... = 0 ⋅ e A ⋅ i [3]
dt A
Evidentemente las dos expresiones deben dar el mismo valor de la ACELERACIÓN (¡compruébelo!).
_________________________________________________________________________________________________________________
4.2. Aplicación de la ecuación de continuidad: criterios de incompresibilidad.

La condición estricta de incompresibilidad, es que la densidad sea constante; no obstante, bajo
determinadas condiciones del flujo, es posible asumir la hipótesis de incompresibilidad. Uno de los criterios, es
que el número de Mach, sea relativamente pequeño, tomando normalmente como límite Ma<0,3. Así en el flujo
en un ventilador, aunque el aire tiene un módulo de compresibilidad bajo, es posible suponer su flujo
incompresible, si el nº de Mach en la punta del alabe , no supera el límite de 0,3.
DETERMINE: La máxima velocidad de giro del ventilador para asumir la hipótesis de incompresibilidad.
DATOS: Diámetro del ventilador = 1 m

Suponer aire ideal a una temperatura de 32 ºC; Raire ideal = 287 J/kgK; cp/cv = 1,4
RESOLUCIÓN:
∂ρ r
Ec. Continuidad + ∇ ⋅ ( ρv ) = 0
∂t
∂ρ
Hipótesis: flujo estacionario ⇒ =0
∂t
r d (ρ u )
Flujo unidimensional ⇒ ∇ ⋅ (ρ v ) =
dx
Con lo que la ecuación de continuidad queda en flujo compresible (con las hipótesis anteriores), como:
d ( ρu ) du dρ
=0 ≡ ρ +u =0 [1 ]
dx dx dx
du
Si el flujo, se hubiese considerado además incompresible, la ecuación de continuidad sería: ρ = 0 [2 ]
dx
du dρ
Con la condición ρ >> u , la ecuación [1] se convierte “casi” en la [2]; es decir si se cumple la
dx dx
desigualdad, aunque la densidad varíe, puede considerarse el flujo como “cuasi-incompresible”.
u dρ
Para flujo unidimensional, la desigualdad anterior, queda como: << 1 [3]
ρ du
La variación de densidad se pueden expresar en función de la velocidad que tienen pequeñas perturbaciones en el
flujo, que se denomina velocidad sónica7 “a”:
dp dp
a= ⇒ dρ =
dρ 2
a
7
Estrictamente las variaciones de presión debidas a las variaciones de densidad, dependen del tipo de proceso
que tenga lugar; en este caso el movimiento de las pequeñas perturbaciones en el seno de un fluido, es rápido y
prácticamente no hay tiempo para la transmisión de calor (además las irreversibilidades son también
despreciables), es decir se puede considerar que el proceso es isentrópico, definiéndose por tanto la velocidad
⎛ ∂p ⎞
sónica como: a= ⎜ ⎟
⎝ ∂ρ ⎠ s
_________________________________________________________________________________________________________________
La variación de velocidad se puede expresar en función de la variación de presión (Ec. de BERNOULLI sin
variaciones de energía potencial):
⎛ u2 ⎞ dp
dp + d⎜ ⎟=0 ⇒ du = −
⎜ 2 ⎟ ρu
⎝ ⎠
Con las dos últimas expresiones se puede rescribir la desigualdad [3], quedando:
dp
2
u a2 ⎛u⎞
= ⎜ ⎟ = Ma 2 << 1 ≡ Ma 2 << 1
ρ −dp ⎝ a ⎠
ρu
La relación entre la velocidad de una partícula del flujo, y la velocidad de pequeñas perturbaciones en el seno del
u
fluido, es el número de MACH, “Ma”: Ma =
a
Con todo, se tiene que el criterio de flujo “cuasi-incompresible” es Ma2<<1. Normalmente se toma Ma<0,3
(Ma2<0,09), como criterio de incompresibilidad.
En el caso de un ventilador, el mayor Ma del flujo, se obtiene en las partículas de aire que tocan las puntas de los
alabes; en donde prácticamente adquieren la velocidad tangencial de la punta (ωR).
La velocidad de giro máxima del ventilador, para poder considerar el flujo como incompresible, viene
determinada por la condición de que el número de Ma en cualquier punto del flujo sea menor de 0,3; y como las
partículas más rápidas son las de las puntas, se tiene que el flujo se podrá considerar como incompresible, si el
Ma en la punta del alabe no supera el valor de 0,3:
ωR 0,3 a
Ma punta = < 0,3 ⇒ ω< [4]
a R
La velocidad del sonido en el aire, se puede determinar considerándolo como gas ideal, con lo que se tiene:
γ aire = c p / c v = 1, 4
a aire ideal = γRT = ... R aire = 287 J / kgK ... a = 1,4 ⋅ 287 ⋅ 305,15 = 350,156 m / s
T = 32 + 273,15 = 305,15 K
con lo que da la ecuación [4] se obtiene el máximo valor de la velocidad de giro, para que el flujo se pueda
considerar como incompresible:
0,3 a a = 350,156 m / s 0,3 ⋅ 350,156 −1

ω< = ... ... ω< = 210,094 s
R R = 0,5m 0,5
60
N < 210,094 = 2006, 25 rpm
2π
El ventilador de 1 m de diámetro, debe girar a una velocidad menor de 2006,25 rpm, para poder considerar la
hipótesis de “cuasi-incompresibilidad” en todo el flujo.
_________________________________________________________________________________________________________________
4.3. Aplicación de las ecuaciones de continuidad y de BERNOULLI: descarga de depósitos.

Uno de los problemas clásicos de hidrodinámica, es la descarga de un deposito. Se suelen tener tres
casos: a) determinación de la velocidad de salida, por un orificio inferior, considerando que el nivel del deposito
es constante; b) tiempo de descarga de un deposito y c) forma del deposito para que la velocidad de descenso del
nivel sea constante (reloj de agua o clepsidra).
DETERMINE para cada caso: 1. velocidad del chorro para flujo no viscoso.
2. tiempo de descarga.
3. forma del deposito de revolución.
4. valores de 1.), 2.) y 3.) a partir de los datos numéricos.
DATOS: nivel del deposito = H.

nivel inicial del deposito =h0 ; diámetro del deposito = D1; diámetro del orificio = D2.
velocidad de descenso del nivel de agua = v; diámetro de salida = D2
numéricos: (a) H=3m
(b) h0= 3m; D1=2m; DO=20mm
(c) v=2 mm / min; D2=4 mm
RESOLUCIÓN:
a) Deposito de nivel constante

1
vd
BERNOULLI (1) (2) p1 + 12 ρ v12 + ρg z1 = p 2 + 12 ρ v 22 + ρg z 2
p1 = p atm
(1)
v1 = 0 (H = cte)
⇒ p atm + 12 ρ 0 2 + ρg z1 = p atm + 12 ρ v d2 + ρg z 2
p 2 = p atm
(2)
v 2 = v desc arg a = v d
Con lo que la velocidad de descarga es: v d = 2g (z 2 − z 1 ) = 2gH v d = 2gH [1]
Que es la ecuación de TORRICELLI. Numéricamente, si el nivel del deposito es de H=3m, la velocidad de

descarga, es:
v d = 2gH = 2 ⋅ 9,8 ⋅ 3 = 7,668 m / s
COEFICIENTE DE CONTRACCIÓN: existe una contracción del chorro de salida, con respecto al área del
orificio de salida; la relación entre las dos áreas, se denomina coeficiente de contracción de la vena líquida:
Cc=AC/A0, que suele ser del orden de 0,65.
_________________________________________________________________________________________________________________
COEFICIENTE DE VELOCIDAD: experimentalmente la velocidad de salida en la contracción, es menor que la

dada por la ecuación de Torricelli, según el coeficiente de velocidad: Cv=(vd)experimental / (vd)Torricelli; que suele ser
del orden de 0,95.
COEFICIENTE DE DESCARGA: el caudal que atraviesa la sección contraída es:
( )
Q = A c ⋅ v d = (A 0 C c ) C v 2gH = (C c C v )A 0 2gH
Al producto CcCv se le denomina coeficiente de descarga Cd, que suele ser del orden de 0,60; y a partir del cual
se puede obtener el caudal que sale por el orificio:
Q = C d A 0 2gH
b) Deposito de nivel variable:
En un determinado instante, si el nivel es h(t), la ecuación de Bernoulli

entre un punto de la superficie libre (1) y un punto del chorro de salida (2) 1 dh
es:
2 h(t) 2
1 ⎡ A ⎛ − dh ⎞⎤
2
1 ⎛ − dh ⎞
p atm + ρ⎜ ⎟ + ρgh = p atm + ρ⎢ 1 ⎜ ⎟⎥
2 ⎝ dt ⎠ 2 ⎣ C d A o ⎝ dt ⎠⎦
en donde: la velocidad del punto (1) es –dh/dt; el signo menos es debido a que dh es <0
A1 − dh
la velocidad del punto (2), por continuidad, es v2 = A1v1/A2 =
c d A o dt
2⎡ 2 ⎤
⎛ − dh ⎞ ⎛ A1 ⎞
la ecuación anterior queda como: ⎜ ⎟ ⎢⎜⎜ ⎟⎟ − 1⎥ = 2gh
⎝ dt ⎠ ⎢⎣⎝ C d A o ⎠ ⎥⎦
dh 2g
=± ⋅ dt cuya integración es:
( )
separando variables, se tiene la ecuación diferencial:
A1 2
h
Cd A o
−1
2
⎛ A1 ⎞
⎜⎜ ⎟⎟ − 1
∫ ∫ ( )
h t
t=2 ⎝ d o⎠
dh 2g C A
= ± 2
⋅ dt ⇒ h0 − h
h0 h 0 ⎛ A1 ⎞ 2g
⎜⎜ ⎟⎟ − 1
⎝ Cd A o ⎠
2
⎛ A1 ⎞
⎜⎜ ⎟ −1
⎝ C d A o ⎟⎠
el tiempo de descarga del deposito, se tiene con h=0: t desc arg a =2 ⋅ h0
2g
πD 2 πD 02
A= ; D = 2m; A0 = ; D 0 = 20mm
numéricamente, con los datos de 4 4 ,
Cd = 0,60 h 0 = 3m
_________________________________________________________________________________________________________________
el tiempo de descarga es:
2 2
⎛ D2 ⎞ ⎛ 22 ⎞
⎜ ⎟ ⎜ ⎟ −1
⎜ C D2 ⎟ −1 ⎜ 0,60 ⋅ 0,020 2 ⎟
=2 ⎝
d o ⎠ ⎝ ⎠
t desc arg a h0 = 2 3 = 13041 s
2g 2 ⋅ 9,8
es decir el deposito tarda en vaciarse 3h 37 m 21s
c) Deposito de nivel y sección variable, con velocidad de disminución de nivel constante:
dh En un determinado instante, si el nivel es h(t), la ecuación de Bernoulli entre

r un punto de la superficie libre (1) y un punto del chorro de salida (2) es:
h(t)
2
1 ⎡ A ⎛ − dh ⎞⎤
2
1 ⎛ − dh ⎞
p atm + ρ⎜ ⎟ + ρgh = p atm + ρ⎢ 1 ⎜ ⎟⎥
2 ⎝ dt ⎠ 2 ⎣ A 2 ⎝ dt ⎠⎦
en donde: la velocidad del punto (1) es –dh/dt; el signo menos es debido a que dh es <0
A1 − dh
la velocidad del punto (2), por continuidad, es v2 = A1v1/A2 =
A 2 dt
Considerando que el recipiente es de revolución, es decir en cada sección, el área viene determinada por πr2, la
ecuación anterior queda como:
2⎡ 2 ⎤
⎛ − dh ⎞ ⎛ 2 ⎞
⎜ ⎟ ⎢⎜ r ⎟ − 1⎥ = 2gh
⎝ dt ⎠ ⎢⎜⎝ r22 ⎟
⎠ ⎥
⎣ ⎦
Se quiere construir el recipiente, de tal forma, que la disminución de nivel sea constante con el tiempo; por
ejemplo, se da el dato de que la velocidad de disminución del nivel es de (2mm/min), es decir cada minuto el
nivel del deposito baja 2mm, con lo que al cabo de 6 horas el nivel bajara 720 mm. El que la velocidad de bajada
del nivel sea constante, permite establecer una escala vertical de longitudes, que mida el nivel, proporcional a la
escala de tiempos: se tiene un reloj de agua.
La velocidad “v” de bajada del nivel es constante e igual a –dh/dt; con lo que de la última ecuación se tiene la
forma que debe tener el recipiente que se utilice como reloj de agua:
v2 ⎡ 1 4 ⎤
h= ⎢ ⋅ r − 1⎥ es decir un paraboloide de revolución de cuarto orden: h ∝ r
4
2g ⎣ r24 ⎦
Numéricamente, el radio del recipiente con h=720mm es:
2gh 2 ⋅ 9,8 ⋅ 0,720

rh =720mm = ro 4 + 1 = 2⋅ + 1 = 671, 4 mm
( 0, 002 / 60 )
4 2
v2
La cantidad de agua que debemos de poner en la clepsidra, para tener un reloj de 6 horas, y que cada minuto el
nivel descienda 2mm, será el volumen del paraboloide:
2πr22 2π ⋅ 0,0022
V(h=720mm a h=0) = ...= 2gh o = 2 ⋅ 9,8 ⋅ 0, 720 = 2,832 m 3
v 0,002 / 60
_________________________________________________________________________________________________________________
4.4. Aplicación de la ecuación de BERNOULLI: flujo no viscoso entre discos horizontales.

El flujo radial entre dos discos horizontales y concéntricos, da lugar a un campo de presiones, con
valores menores que la presión atmosférica, lo que hace que la fuerza de presión, pueda ser suficiente para
sujetar el disco inferior. Se supone que el caudal es aportado por una tubería vertical de diámetro D1, y que fluye
entre los discos hasta salir a la atmósfera. Considere el flujo ideal y que la velocidad, desde R1 a R2 es
exclusivamente radial.
DETERMINE: 1. La velocidad de salida.

2. La presión en una posición radial genérica comprendida entre R1 y R2.
3. La fuerza de presión sobre la corona del disco inferior (D1 a D2).
DATOS: Fluido: densidad constante = ρ; caudal = Q

Geometría: diámetro de la tubería = D1; discos: diámetro = D2; huelgo = H.
Numéricos: ρ=900kg/m3; Q = 0,4 litros/segundo ; D1=30mm; D2=400mm; H=10mm
RESOLUCIÓN:
(1) Velocidad de salida: se conoce el caudal de salida por los discos (igual al que se aporta por la tubería
vertical) y se conoce la sección de salida (área lateral del estrecho cilindro de salida = 2πR2H) con lo que
evidentemente:
Q
v2 =
2πR 2 H
numéricamente, con los datos de Q=4 litros/segundo; R2=200mm y H=10mm, se tiene que la velocidad de salida
es:
−3
4 ⋅ 10
v2 = = 0,318 m / s
2π ⋅ 0,200 ⋅ 0,010
Q
r
H A 2
R1
R2
A2=2πR2 H
En cada sección (r=cte), el perfil de velocidades es uniforme, por ser un fluido no viscoso (precisamente por eso,
la condición de no deslizamiento no se cumple, es decir el fluido que toca las paredes no se pone a la velocidad
de ellas, lo que si se cumple es que no puede atravesarlas).
Si el fluido fuese viscoso, se cumple la condición de no deslizamiento y las partículas del flujo que tocan las
paredes de los discos, estarían a velocidad nula, al estar los discos quietos, se tiene máxima velocidad en la línea
central de los discos, y esa velocidad máxima va disminuyendo conforme el fluido avanza (ver problema 3.8)
_________________________________________________________________________________________________________________
(2) Presión en una posición radial genérica (entre R1 y R2)
Consideraremos que en las posiciones radiales entre R1 y R2 se tiene flujo irrotacional, y como la densidad es
constante, se cumple la ecuación de BERNOULLI8. La establecemos entre un punto genérico (A) y un punto en la
salida del flujo a la atmósfera (2):
2 2
p + 1 ρ v + ρgz = p atm + 1 ρ v 2 + ρgz
2 2
con lo que se tiene que la presión es una posición radial genérica es:
⎡ 2 2⎤ ⎛ 1 1 ⎞⎟
( )
ρ ⎛ Q ⎞
2
2 ρ 2 ⎢ ⎜ ⎟ ⎛ Q ⎞ ⎥ ρQ ⎜
p = p atm + v 2 − v = p atm + −⎜ ⎟ = p atm + 2 2 −
2 2 ⎢⎜⎝ 2 πR 2 H ⎟⎠ ⎝ 2πrH ⎠ ⎥ 8π H ⎜ R2 r2 ⎟
⎣ ⎦ ⎝ 2 ⎠
ρQ
2 ⎛ 1 1 ⎞⎟
p( r ) = p atm + ⎜ − [1]
2 2
8π H ⎜ R2 r2 ⎟
⎝ 2 ⎠
como en todas las posiciones radiales r<R2, el termino del paréntesis es siempre negativo, es decir que la presión
en cualquier posición radial es menor que la atmosférica(e igual a la atmosférica en la salida).
Hemos considerado como sección de inicio del flujo irrotacional, al correspondiente al radio de la tubería
vertical de aporte de caudal. La presión en la citada sección es:
ρQ
2 ⎛ 1 1 ⎞⎟ 900 ⋅ 0,004
2 ⎛ 1 1 ⎞
p( R 1 ) = p atm + ⎜ − = p atm + ⎜⎜ − ⎟ = p atm − 8060,1 Pa
2 2 ⎜ R2 R2 ⎟ 2 2 2 2⎟
8π H ⎝ 2 1 ⎠ 8π 0,010 ⎝ 0,200 0,015 ⎠
es decir que la presión vacuométrica en la sección inicial es de 8060,1 Pa
8
En el interior de los discos, en posiciones radiales próximas al eje de la tubería el flujo es rotacional, no
pudiendo establecer la ecuación de Bernoulli, conforme el flujo avanza por el interior de los discos, por ser el
huelgo muy pequeño, el fluido se va haciendo irrotacional, y al ser incompresible, se puede aplicar la ecuación
de Bernoulli. En este problema consideraremos como hipótesis simplificativa, que se tiene flujo irrotacional a
partir de radios mayores o iguales a R1
_________________________________________________________________________________________________________________
(3) Fuerza de presión sobre la corona del disco inferior (D1 a D2):
Sobre el disco inferior están actuando: por la parte exterior la presión atmosférica y por la parte interior la
presión del líquido, que es menor o igual a la atmosférica; por lo tanto la resultante de las fuerzas de presión
sobre el citado disco tiene sentido positivo de fuerza de sustentación. Sobre una determinada corona elemental,
las fuerzas de presión serán:
dFp = (p atm − p( r ) ) ⋅ dA
y de la ecuación [1] que da la presión en el interior de los discos, se tiene que la fuerza de elemental de presión
sobre el área elemental considerada (dA=2πr dr) es:
− ρQ 2
dFp = ⋅ 2πr dr
8π 2 H 2
dr
p(r) p(r)
r
patm patm dA=2πr dr
con lo que la resultante de las fuerzas de presión, sobre la corona considerada (D1 a D2), es:
⎛ D2 − D2 ⎞
R2 R2 − ρQ 2 ⎛⎜ 1 1 ⎞ − ρQ 2 ⎜ 2 1 − ln D 2 ⎟
Fp = ∫R 1
dFp = ∫R1
− ⎟ 2π r dr =... =
8π 2 H 2 ⎜⎝ R 22 r 2 ⎟⎠ 4πH 2
⎜
⎜ 2D 2 D1 ⎟
⎟
⎝ 2 ⎠
Numéricamente, a partir de los datos, la resultante de la fuerza de presión sobre el disco inferior es:
− ρQ
2 ⎛ D 22 − D12 D 2 ⎞ 900 ⋅ 0,004 2 ⎛ 0,400 2 − 0,030 2 0, 400 ⎞
Fp = ⎜ − ln ⎟= ⎜ − ln ⎟ = 23,985 N
4πH
2 ⎜ 2D 2 D1 ⎟ 4 π ⋅ 0,010 2 ⎜ 2 ⋅ 0,400 2 ⎟
⎝ 2 ⎠ ⎝ 0 , 030 ⎠
El signo positivo, evidencia que la resultante de las fuerzas de presión es una fuerza de sustentación (vertical
hacía arriba).
_________________________________________________________________________________________________________________
4.5. Aplicación de la ecuación de Bernoulli en flujo no estacionario: oscilaciones en un

tubo en U.
En flujo no viscoso, las oscilaciones de un líquido en un tubo en U, son senoidales puras; y la frecuencia
de la oscilación depende exclusivamente de la longitud del tubo mojado. Si el flujo es viscoso, las oscilaciones
se van amortiguando hasta alcanzar el equilibrio con los meniscos en las misma cota.
DETERMINE: 1. La ecuación de la cota de los meniscos en función del tiempo: z=z(t).

2. El periodo de las oscilaciones.
DATOS: fluido: incompresible y no viscoso

Longitud de la columna del tubo mojada por el líquido = L
Desequilibrio inicial: velocidad nula y diferencia de cotas entre los meniscos = H
Numéricos: L = 400mm: H = 100 mm
RESOLUCIÓN:
Ecuación de Bernoulli en régimen no estacionario es: z
∫1
2 ∂v
∂t
(
⋅ ds + 12 v 22 − v12 + ) ∫ 2 dp
1 ρ
+ g ( z 2 − z1 ) = 0
L
En el caso de las oscilaciones de una columna de líquido en un tubo en U, la
densidad es constante, y en cada instante de tiempo todos los puntos de la
línea de corriente central, deben tener la misma velocidad: las partículas del
líquido se mueven mantenido siempre la misma distancia relativa, es decir la longitud total (L) de la columna de
líquido es constante; con lo que en la ecuación anterior, aplicada a la línea de corriente central, y entre los puntos
correspondientes a cada uno de los meniscos, el término de la integral del elemento de distancia de la línea de
corriente como:
2 ∂v 2
∫ ∫
dv dv
⋅ ds = ⋅ ds = ⋅L
1 ∂t dt 1 dt
Con respecto a la línea horizontal de equilibrio, la distancia de los meniscos es la misma, es decir: lo que
asciende un menisco es lo que desciende el otro, por lo tanto se tiene que los módulos de las velocidades de los
meniscos son los mismos, con lo que no hay variación de la energía cinética:
v1 =
dz
dt
∧ v2 =
dz
dt
⇒ v12 = v 22 ⇒ 1
2
(v 2
2 )
− v12 = 0
La densidad es constante, y los dos meniscos están a la presión atmosférica, con lo que el término de trabajo de
flujo es nulo:
2 dp p 2 − p 2 p atm − p atm
∫1 ρ
=
ρ
=
ρ
=0
En cada instante de tiempo, si la cota del menisco izquierdo (punto 2) es “z” (respecto al nivel de equilibrio), la
cota del menisco derecho (punto 1) es “-z”; con lo que la diferencia entre las cotas es:
z2-z1 = z – (-z) = 2z
dv
Con todo lo anterior, la ecuación de Bernoulli queda: ⋅ L + 2g z = 0
dt
_________________________________________________________________________________________________________________
Expresando la velocidad de un punto de la línea de corriente, como: v = dz/dt; queda la ecuación diferencial:
d2z 2g ⎛ 2g ⎞ ⎛ 2g ⎞
=− ⋅z ; cuya solución es: z = A ⋅ cos⎜ ⋅ t ⎟ + B ⋅ sen⎜ ⋅ t⎟ [1 ]
dt 2 L ⎜ L ⎟ ⎜ L ⎟
⎝ ⎠ ⎝ ⎠
en donde las constantes de integración se obtienen por las condiciones de contorno: en el instante inicial el
desequilibrio entre los dos meniscos es “H” y los dos meniscos estas quietos. A partir de ese instante, los
meniscos empiezan a moverse, de tal forma que si inicialmente el m
t = 0; z = H/2 y dz/dt =0
⎛ 2g ⎞ ⎛ 2g ⎞
en la ecuación [1], si z = H/2 y t=0, queda: H / 2 = A ⋅ cos⎜ ⋅ 0 ⎟ + B ⋅ sen⎜ ⋅ 0⎟ = A
⎜ L ⎟ ⎜ L ⎟
⎝ ⎠ ⎝ ⎠
es decir la constante de integración A = H/2
La derivada temporal de la ecuación [1] da la velocidad de movimiento de los meniscos:
dz 2g ⎛ 2g ⎞ 2g ⎛ 2g ⎞
v= = −A ⋅ sen⎜ ⋅ t⎟ + B ⋅ cos⎜ ⋅ t⎟
dt L ⎜ L ⎟ L ⎜ L ⎟
⎝ ⎠ ⎝ ⎠
con la condición inicial que v = 0 en t = 0; se tiene:
2g ⎛ 2g ⎞ 2g ⎛ 2g ⎞ 2g
0 = −A ⋅ sen⎜ ⋅ 0⎟ + B ⋅ cos⎜ ⋅ 0⎟ = B ⇒ B=0
L ⎜ ⎟ ⎜ ⎟
⎝ L ⎠ L ⎝ L ⎠ L
Con todo, la ecuación que da la posición de los meniscos en función del tiempo es:
H ⎛ 2g ⎞
z= ⋅ cos⎜ ⋅t⎟
2 ⎜ L ⎟
⎝ ⎠
Es decir una función senoidal, que se repite periódicamente, con valores del termino del coseno igual a 2kπ; es
decir que el periodo es:
2g 2L
⋅ T = 2π ⇒ T=π
L g
numéricamente, si la longitud de la columna de líquido es L=400mm; el periodo de las oscilaciones periódicas

puras es:
2L 2 ⋅ 0,400
T=π =π = 0,898 s
g 9,8 z
y la amplitud máxima es H= 100 mm

t
H
Las oscilaciones de los meniscos, son de periodo constante y de amplitud constante. Esto es considerando que el
fluido es no viscoso. Si el fluido es viscoso, las oscilaciones son amortiguadas y desaparecen a lo largo del
tiempo.
_________________________________________________________________________________________________________________
4.6. Aplicación de la ecuación de Bernoulli con aceleración de arrastre: bomba rotativa.

Un tubo en forma de L, se sumerge parcialmente en un líquido el tramo vertical, dejando el tramo
horizontal a una cierta altura de la superficie libre. En el tramo vertical se hace girar con una velocidad
constante. Inicialmente el tubo esta lleno de líquido y tiene una válvula en el extremo no sumergido; en un
determinado instante abre la válvula y el líquido empieza a salir por el tramo horizontal.
DETERMINE: 1. Ecuación de Bernoulli considerando la aceleración de arrastre.

2. Caudal bombeado.
3. Velocidad máxima de giro para que el punto de mínima presión no cavite.
4. Velocidad mínima de giro para que el líquido pueda salir.
DATOS: fluido: incompresible y no viscoso

Tubo en forma de L: longitud del tramo horizontal = Lh;
longitud del tramo vertical no mojado = Lv
Presión atmosférica: Pa; presión de vapor: Pv
Numéricos: L h=250mm; Lv=400mm; Pa=1013 mbar; Pv=12 mbar
RESOLUCIÓN:
ω
3
2
Q
R v
H
patm
(1) ECUACIÓN DE BERNOULLI CON ACELERACIÓN DE ARRASTRE:

r
r dv
En la ecuación de movimiento de EULER: ρg − ∇p = ρ , tanto el gradiente de presión, como la
dt
aceleración de las partículas, están referidos a un determinado sistema de referencia. Si el sistema de referencia
es inercial (fijo o se mueve a velocidad constante), se mueve sin aceleración, con lo que las fuerzas másicas sólo
vienen afectadas por el vector aceleración genérica de campos de fuerzas externos.
En cambio, si el sistema de referencia, en el que se tiene los campos de presión y velocidad, es no

inercial (tiene una aceleración de arrastre, respecto a un sistema inercial), para poder mantener en la ecuación
anterior, el gradiente de presión y la aceleración de las partículas (respecto al sistema no inercial), hay que
r r
afectar a las fuerzas másicas, por la aceleración relativa, es decir por g − aSRNI .
Con lo que la ecuación de EULER para un sistema de referencia no inercial es:

r
r r dv
ρ ( g - aSRNI ) − ∇p = ρ
dt
_________________________________________________________________________________________________________________
La ecuación de EULER, para puntos de una línea de corriente, es la ecuación de BERNOULLI, que para
un sistema de referencia no inercial queda como:
∂v
⋅ ds + d ( v 2 ) + + g ⋅ dz − ( a arrastre ) ⋅ ds = 0
1 dp r r
∂t 2 ρ
En el problema, el termino de la aceleración de arrastre, se simplifica notablemente, debido a que el

vector velocidad de giro del sistema de referencia es constante y perpendicular al vector de posición.
Las hipótesis restrictivas son:

∂v
flujo estacionario: =0
∂t
flujo incompresible: ρ = cte.
r r r
aceleración de arrastre, exclusivamente centrípeta:
r
( )r
a arrastre = ωk × ωk × ru axial = ω2 r ⋅ u axial
r r
línea de corriente en el tramo horizontal, dirección axial: ds = dr ⋅ u axial
con lo que el sumando correspondiente a la aceleración de arrastre, en la dirección de la línea de corriente es:
r r
( a arrastre ) ⋅ ds = ω2 r ⋅ dr
d ( v 2 ) + + g ⋅ dz −ω2 r ⋅ dr = 0
1 dp
con todo la ecuación de BERNOULLI queda:
2 ρ
v2 p ω2 r 2
cuya integración da: + + gz − = cte. ≡ p + 12 ρ v 2 − 12 ρ ω 2 r 2 + ρgz = cte. [1]
2 ρ 2
(2). CAUDAL BOMBEADO.
La ecuación de Bernoulli [1] entre el punto (2) y el punto (1) (a tal profundidad que su velocidad es nula), es:
p 2 + 12 ρv 22 − 12 ρω 2 ⋅ 0 2 + ρgz 2 = p1 + 12 ρ ⋅ 0 2 − 12 ρω 2 ⋅ 0 2 + ρgz1
A su vez, como en el punto (1) la velocidad es nula, su presión es la atmosférica más la debida a la columna de
agua que tiene encima: p1=patm+ρg(zSL-z1), en donde zSL es la cota de la superficie libre; con lo que la queda:
p 2 + 12 ρv 2 + ρgz 2 = p atm + ρg(z SL − z1 ) + ρgz1
siendo H = z2 – zSL;
la presión en el punto (2) es: p 2 = p atm − 12 ρv 2 − ρgH [2]
La ecuación de Bernoulli [1] entre el punto (2) y el punto (3) (chorro libre a la presión atmosférica), es:
p 2 + 12 ρv 22 − 12 ρω 2 ⋅ 0 2 + ρgz 2 = p atm + 12 ρ ⋅ v 22 − 12 ρω 2 ⋅ R 2 + ρgz 2
La velocidad en cualquier sección del interior del tubo es la misma, e igual al caudal por la sección:
4Q
v=
πD 2
Con lo que la presión en el punto (2) es: p 2 = p atm − 12 ρω 2 R 2 [3]
_________________________________________________________________________________________________________________
De la igualdad de las dos expresiones que dan la presión en el punto (2) se obtiene la velocidad de la corriente
uniforme:
[2] = [3] p atm − 12 ρv 2 − ρgH = p atm − 12 ρω 2 R 2 v = ω 2 R 2 − 2gH [4]
πD 2
con lo que el caudal bombeado es: Q= ω 2 R 2 − 2gH [5]
4
Numéricamente, con los datos: D=30mm; R=200mm; H=250mm; ω=200(2π/60)=20,944s-1:
v = ω 2 R 2 − 2gH = 20,944 2 ⋅ 0,200 2 − 2 ⋅ 9,8 ⋅ 0,250 = 3,556 m / s
πD 2 π ⋅ 0,030 2
Q= ω 2 R 2 − 2gH = ⋅ 3,556 = 2,514 10 −3 m 3 / s = 2,514 litros / segundo
4 4
(3). VELOCIDAD MÍNIMA DE GIRO. En la ecuación que da el caudal en función de la velocidad de giro [5],
si el radicando en nulo el caudal es nulo, y si es negativo quiere decir que no hay suficiente “fuerza centrifuga”
para establecer el flujo; con lo que para que haya caudal la velocidad de giro debe ser:
2gH
ω>
R
numéricamente:
2gH 2 ⋅ 9,8 ⋅ 0,250 60

ω> = = 11,068 s −1 ≡ N > 11,068 = 105,69 rpm
R 0,200 2π
(4) VELOCIDAD MÁXIMA DE GIRO: en todo el tubo, la presión mínima se tiene en el punto (2), y su valor
viene determinado por la ecuación [3]:
ρω 2 2
p 2 = p atm − R
2
esa presión debe ser siempre mayor que la presión de vapor a la temperatura del flujo, para que el líquido no
empiece a hervir (cavitación); lo que lleva a poder obtener la velocidad máxima a la que puede girar el sistema:
ρω 2 2 2(p atm − p vapor ) / ρ

p 2 = p atm − R > pvapor ⇒ ω≤
2 R
numéricamente, con los datos: patm = 1020 mbar; pvapor = 23 mbar; ρ=1000 kg/m3; R=0,2m
2(p atm − p vapor ) / ρ 2(102000 − 2300) / 1000

⇒ ω≤ = = 70.605 s −1
R 0,2
N < 70,605(60/2π) = 674,2 rpm
Es decir que la velocidad de giro deber ser menor de 674,2 rpm; si fuese mayor el punto (2) tendría una presión
igual a la de vapor, y el líquido “herviría” (cavitación).
_________________________________________________________________________________________________________________
4.7. Aplicación de las ecuaciones de Navier-Stokes: flujo de Couette-Poiseuille

Una fina película de líquido de espesor constante, fluye estacionariamente por un plano inclinado. La atmósfera
ejerce una presión constante en todas las secciones del flujo y es despreciable el esfuerzo viscoso en la superficie
libre.
DETERMINE: 1. El perfil de velocidades: u=u(y).

2. La ecuación del caudal volumétrico.
3. A partir de los datos numéricos, el ángulo del plano inclinado.
DATOS: Flujo: estacionario, incompresible y laminar (v=w=0)

Plano inclinado: ángulo de inclinación = θ
Película: espesor = H; anchura = a
Numéricos: viscosidad = 2,3 St; H=10mm; a=1m; Q= 8 lps
RESOLUCIÓN:
y
H
u
Q
α
x
(1)PERFIL DE VELOCIDADES:
En la dirección de la corriente (x), la ecuación de Navier-Stokes, para flujo incompresible y estacionario es:
∂p ⎛ ∂ 2u ∂ 2u ∂ 2u ⎞ ⎛ ∂u ∂u ∂u ⎞
ρ⋅gx − + µ⎜⎜ 2 + 2 + 2 ⎟⎟ = ρ⎜⎜ u +v + w ⎟⎟ [1]
∂x ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂x ∂y ∂z ⎠
particularizando para el problema:

gx
componente gravitacional: g α
g x = g ⋅ senα
gradiente de presión: todos los puntos de la superficie libre están a la presión atmosférica, y considerando que
el espesor de la capa de fluido es pequeño, se puede considerar que todos los puntos del flujo están a la presión
atmosférica, es decir no hay gradiente de presión: ∇p = 0 ; con lo que se tiene que en cualquier dirección la
presión no varía, y en particular en la dirección del flujo:
∂p
=0
∂x
∂p ∂p
Estrictamente: =0 ∧ = −ρg cos α ⇒ p = p atm + ρg cos α(H − y )
∂x ∂y
Pero como el huelgo es muy pequeño, ρgcosα(H-y) << patm, y con ello p = patm = cte.
_________________________________________________________________________________________________________________
vector velocidad: el vector velocidad sólo tiene componente en la dirección “x”; con lo que v = w = 0; por lo
que por continuidad:
∂u ∂v ∂w
+ + =0 ∂u
∂x ∂y ∂z ⇒ =0 ⇒ u = u ( y)
∂x
v=w=0
quedando de la ecuación de Navier-Stokes [1], la ecuación diferencial que permite obtener el perfil de
velocidades:
d2u − ρg ⋅ senα 2
ρg ⋅ senα + µ ⋅ =0 cuya solución es: u= ⋅y +A⋅y+B
dy 2
2µ
las constantes de integración, se obtienen con las condiciones de contorno: en la superficie inclinada (que esta
fija) la velocidad es cero, y en la superficie libre la tensión de rozamiento con el aire ambiente es nula (se
desprecia el efecto de rozamiento del aire9) y con ello la velocidad es máxima:
y=0 ⇒ u=0
ρg ⋅ senα
du B=0; A= H
y=H ⇒ τ=0 ≡ =0 µ
dy
con lo que el perfil de velocidades es:
u=
−ρg ⋅ senα 2 ρg ⋅ senα
2µ
⋅y +
µ
H⋅y u=
ρg ⋅ senα
2µ
(
⋅ 2Hy − y 2 ) [2]
(2) CAUDAL VOLUMÉTRICO DE LA CAPA DE ACEITE.

r r r r
Al ser: v = u i ; dA = (a ⋅ dy ) i ; en donde “a” es la anchura de la capa de aceite, el caudal se obtiene por:
Q=
∫∫A
r r y=H
∫
v ⋅ dA = u ⋅ ady =
y =0 0
H ρg ⋅ senα
2 µ ∫ (
2Hy − y 2 ⋅ ady = ... = )
g ⋅ senα 3
3ν
aH Q=
g ⋅ senα
3ν
aH 3 [3]
(3) ÁNGULO DEL PLANO INCLINADO, PARA QUE EL CAUDAL SEA DE 8 litros/segundo:
De la ecuación del caudal [3], el ángulo del plano inclinado, para que el caudal de líquido sea de 8 lps es:
⎛ 3Qν ⎞ ⎛ −3 −4 ⎞
α = arcsen⎜ ⎟ = arcsen⎜ 3 ⋅ 8 ⋅ 10 ⋅ 2,3 ⋅ 10 ⎟ = 34,28º
⎜ g a H3 ⎟ ⎜ 9,8 ⋅ 1 ⋅ 0,010 3 ⎟
⎝ ⎠ ⎝ ⎠
la viscosidad cinemática es de 2,3 St (STOKES) = 2,3 cm2/s = 2,3 10-4 m2/s
9
Estrictamente, en la superficie libre (interfase aceite-aire) se cumple que las tensiones son iguales: en nuestro
⎛ ∂u 1 ⎞ ⎛ ∂u 2 ⎞
caso: µ1 ⎜ ⎟ = µ2 ⎜ ⎟ . Se recuerda que para fluidos inmiscibles, en la interfase, aunque las
⎝ ∂y ⎠ y =H ⎝ ∂y ⎠ y = H
tensiones sean iguales, las velocidades tangenciales son distintas, por la diferencia de viscosidad.
_________________________________________________________________________________________________________________
4.8. Aplicación de las ecuaciones de NAVIER-STOKES: flujo viscoso laminar entre discos
El flujo radial entre dos discos horizontales y concéntricos, da lugar a un campo de presiones, que en
función del Re, puede tener un gradiente de presión positivo o negativo. En el caso de flujo laminar, se obtiene
que el gradiente de presión es negativo, con lo que la presión va disminuyendo conforme aumenta el radio. Con
lo cual, si se descarga a la atmósfera, esa será la sección de mayor radio y de menor presión, con lo que en
secciones interiores, la presión será mayor que la atmosférica, que es justamente lo contrario a lo que ocurría en
el problema 4.
Se supone que el caudal es aportado por una tubería vertical de diámetro D1, y que fluye entre los
discos hasta salir a la atmósfera. Considere como hipótesis, que entre las posiciones radiales R1 a R2 el flujo es
estacionario, incompresible, viscoso y laminar.
DETERMINE: 1. El perfil de velocidades en cada posición radial: v=v(r,z)..

2. La presión en cada posición radial: p=p(r)
3. La fuerza de presión sobre la corona del disco inferior (D1 a D2).
DATOS: Fluido: densidad y viscosidad constantes; caudal constante = Q

Geometría: diámetro de la tubería = D1; discos: diámetro = D2; huelgo = H.
Numéricos: µ=0,8kg/ms; Q = 0,4 litros/segundo ; D1=30mm; D2=400mm; H=10mm
RESOLUCIÓN:
Q vr=vr(r,z)
2
r z´s
H A
r´s
R1
R2
z=H/2
z=0
z=-H/2
A2=2πR2 H
Si el flujo es laminar, las partículas se mueven en capas horizontales (vz=0) sin vorticidad (vθ=0), y la
velocidad solo tiene componente radial; además las fuerzas de inercia se pueden despreciar frene a las viscosas:
r r r r
v = v r ⋅ ur + 0 ⋅ u + 0 ⋅ uz
⇒ θ
FLUJO LAMINAR Re ↓↓ ⇒ fuerzas viscosas > fuerzas inercia
La ecuación de NAVIER-STOKES para la componente radial, para flujo incompresible y estacionario,

con velocidades tangencial y axial nulas, es:
∂p ⎡ ∂ ⎡1 ∂ (r v r ) ⎤ ∂ 2 v r ⎤ ∂v r
ρg r − + µ⎢ ⎢ ⎥ + ⎥ = ρ v [1]
∂r ⎣ ∂r ⎣ r ∂r ⎦ ∂z ⎦ ∂r
2 r
(1) (2) (3) (4) (5)
_________________________________________________________________________________________________________________
Respecto a cada uno de los 5 términos que aparecen en la ecuación anterior, se tiene:
- El término (1) es nulo, por ser horizontal la dirección radial.
- El término (2) se puede expresar por dp/dr, si consideramos que la presión prácticamente no varia con
las “z´s”“al ser el huelgo muy pequeño.
- El término (5) de fuerzas de inercia (por unidad de volumen) es despreciable frente al termino (3)(4) de
fuerzas viscosas (por unidad de volumen), por ser el flujo laminar.
- El término (3) es nulo por continuidad, pues al ser incompresible, el flujo es adivergente:
r 1 ∂ (r v r ) 1 ∂v θ ∂v z
∇⋅v = 0 = + + ∂ (r v r )
r ∂r r ∂θ ∂z ⇒ =0 [2]
∂r
vθ = v z = 0
2
dp ∂ vr
Con todo la ecuación de Navier-Stokes, queda: − +µ =0 [3]
2
dr ∂z
(1) CAMPO DE VELOCIDADES:
De la ecuación de continuidad [2], se puede deducir, que el producto “r vr” no depende del radio, y por tanto
dependerá exclusivamente de la altura “z”, es decir:
∂ (rv r ) r v r ≠ f (r )
=0 ⇒ ⇒ r v r = f (z)
∂r v r = v r ( z, r )
1
con lo que la velocidad radial es una función del tipo: vr = f (z) [4]
r
2⎛1 ⎞
2 ∂ ⎜ f (z) ⎟ 2 2
∂ vr ⎝r ⎠ = 1 ∂ f (z) = 1 d f (z)
cuya segunda derivada, respecto a “z” es: =
2 2 2 2
∂z ∂z r ∂z r dr
dp 1 d 2 f (z)
que llevado a la ecuación [3], queda: − +µ =0
dr r dz 2
r ⎛ dp ⎞ z 2
cuya integración da la función “f(z)”: f (z) = ⎜ ⎟ + Az + B
µ ⎝ dr ⎠ 2
con lo que de la ecuación [4], se puede obtener la expresión de la velocidad radial:
1 ⎛ dp ⎞ 2 A B
vr = ⎜ ⎟z + z + [5]
2µ ⎝ dr ⎠ r r
_________________________________________________________________________________________________________________
Las constantes de integración se obtienen con las condiciones de contorno: por la condición de no deslizamiento,
en las paredes, la velocidad es nula: vr=0, con z=-H/2 y z=H/2; que llevadas a la ecuación [5], resulta:
− 1 ⎛ dp ⎞ H 2
A=0 y B= ⎜ ⎟
8µ ⎝ dr ⎠ r
1 ⎛ −dp ⎞ 2
vr = ⎟ ( H − 4z )
2
Con lo que se tiene como expresión de la velocidad radial: ⎜ [6]
8µ ⎝ dr ⎠
Como las velocidades son siempre positivas, y 4z2 siempre es menor que H2; el gradiente radial de presiones
debe ser negativo, es decir la presión va disminuyendo con el radio. Como en la salida se tiene la presión
atmosférica, la presión en el interior siempre es mayor que la atmosférica (¡compare y comente con el
resultado del mismo problema con flujo no viscoso¡).
Una vez obtenido el perfil de velocidades, en cada sección, el gradiente de presión (dp/dr) se puede determinar, a
partir de la ecuación del caudal volumétrico:
r r
r r v = vr ⋅ u r z = H 1 ⎛ dp ⎞
r ... = ∫ z =− H / 2 v r 2πr dz = ∫ z = 0 ⎜ ⎟ ( 4z − H ) 2πr dz = ...
z=H / 2
Q = ∫∫ A v ⋅ dA = ... r 2 2
dA = 2πr dz ⋅ u r 8µ ⎝ dr ⎠
πH 3 ⎛ −dp ⎞ dp 6µQ 1
Q= ⎜ ⎟ r ⇒ =− [7]
6µ ⎝ dr ⎠ dr πH 3 r
Con la expresión del gradiente de presión, y la ecuación [6] se obtiene el campo de velocidades:
r 3Q H 2 − 4z 2 r
v= ⋅ ur [8]
4πH 3 r
(2) CAMPO DE PRESIONES:
La integración del gradiente de presiones, permite obtener el campo de presiones:
−6µQ 1 p = patm r =R2 −6µQ dr −6µQ R

dp = dr ⇒ ∫p = p dp = patm − p(r) = ∫r = r πH3 r = πH3 ln r
2
πH 3 r
6µQ R
p(r) = patm + ln 2 [9]
πH 3
r
En la ecuación anterior, como r<R2, el logaritmo es positivo, y con ello la presión en cualquier posición radial es
mayor que la presión atmosférica; evidentemente con r=R2, se tiene la presión atmosférica.
_________________________________________________________________________________________________________________
(3) Fuerza de presión sobre la corona del disco inferior (D1 a D2):
Sobre el disco inferior están actuando: por la parte exterior la presión atmosférica y por la parte interior la
presión del líquido, que es mayor o igual a la atmosférica; por lo tanto la resultante de las fuerzas de presión
sobre el citado disco tiene sentido hacia abajo. Sobre una determinada corona elemental, las fuerzas de presión
serán:
dFp = (p atm − p( r ) ) ⋅ dA
y de la ecuación [9] que da la presión en el interior de los discos, se tiene que la fuerza de elemental de presión
sobre el área elemental considerada (dA=2πr dr) es:
−6µQ R 2
dFp = ln ⋅ 2πr dr
πH 3 r
p(r) p(r)
dr
r
patm patm dA=2πr dr
con lo que la resultante de las fuerzas de presión, sobre la corona considerada (D1 a D2), es:
R R −6µQ R 2 −12µQ ⎡ 1 ⎛ 2 2⎛ R 2 ⎞ ⎞⎤
Fp = ∫ R 2 dFp = ∫ R 2 2π r dr =... = ⎢ ⎜ R 2 − R1 ⎜ 1 + 2 ln ⎟⎥
R1 ⎟⎠ ⎠ ⎦
ln
1 1
πH 3
r 3
H ⎣4 ⎝ ⎝
−3µQ ⎡ 2 2⎛ R 2 ⎞⎤
Fp = 3 ⎢ R 2 − R1 ⎜1 + 2 ln ⎟⎥
H ⎣⎢ ⎝ R1 ⎠ ⎦⎥
Numéricamente, a partir de los datos, la resultante de la fuerza de presión sobre el disco inferior es:
−3µQ ⎡ 2 ⎛ R ⎞ ⎤ −3 ⋅ 0,8 ⋅ 0, 410−3 ⎡ 2⎛ 0, 200 ⎞ ⎤

Fp = ⎢ R − R 2 ⎜1 + 2 ln 2 ⎟ ⎥ = 2
⎢0, 200 − 0, 015 ⎜1 + 2 ln ⎟ ⎥ = −9,3 N
H ⎢⎣
3 2 1 ⎜ R ⎥⎟ 3 ⎣ ⎝ 0, 015 ⎠ ⎦
⎝ 1 ⎠⎦ 0, 010
El signo negativo, evidencia que la resultante de las fuerzas de presión es una fuerza vertical hacía abajo.
_________________________________________________________________________________________________________________
4.9. Aplicación de la ecuación de ENERGÍA: distribución de temperaturas en flujo de

Poiseuille.
Un líquido fluye en régimen laminar entre dos placas horizontales fijas; la placa inferior es una pared
adiabática y la placa superior es una pared isoterma. El flujo esta provocado por una gradiente horizontal de
presión.
DETERMINE: 1. El perfil de velocidades : u=u(y)..
2. El gradiente de presión en función del caudal.
3. El perfil de temperaturas: T=T(y)
4. Potencia calorífica que atraviesa la pared isoterma.
DATOS: Fluido: densidad, viscosidad y conductividad térmica, constantes = ρ, µ , κ

Flujo: caudal constate = Q; laminar (v=w=0)
Geometría: huelgo: H; ancho: A; longitud: L ; tisoterma=60ºC
Numéricos: ρ=900 kg/m3; µ=0,2kg/ms; κ=0,3W/mK; Q = 1 lps; H=10mm; L=A=1m.
RESOLUCIÓN:
pared isoterma
y´s y=H/2
T(y)
Q u(y) Q x´s
y=-H/2
pared adiabática
(1) PERFIL DE VELOCIDADES:
CONTINUIDAD:
∂ρ ∂ (ρu ) ∂ (ρv ) ∂ (ρw )

+ + + =0
∂t ∂x ∂y ∂z
∂u
incompresible ⇒ ρ = cte ⇒ =0 ⇒ u ≠ u(x) ⇒ u = u ( y)
∂x
la min ar ⇒ v = w = 0
NAVIER-STOKES, en la dirección horizontal:
⎛ ∂ 2 u ∂ 2 u ∂ 2 u ⎞ ⎛ ∂u
ρg x −
∂p
+ µ⎜ + + ⎟ = ρ⎜ + u ∂u + v ∂u + w ∂u ⎞⎟
∂x ⎜ 2 2 2 ⎟ ⎝ ∂t ∂x ∂y ∂z ⎠
⎝ ∂x ∂y ∂z ⎠
d 2 u 1 ∂p
con gx=0; v=w=0; y u=u(y); la ecuación anterior queda como: = [1]
dy 2 µ ∂x
_________________________________________________________________________________________________________________
La integración de la ecuación diferencial anterior, nos permite obtener el campo de velocidades:
1 ∂p 2
u= y + C1 ⋅ y + C 2
2µ ∂x
y las constantes de integración, se obtienen con las condiciones de contorno, de que en las paredes fijas, la
velocidad de las partículas es nula: y=H/2; u=0; e y=-H/2; u=0; quedando:
1 ⎛ −∂p ⎞ 2
u= ⎟ ( H − 4y ) [2]
2
⎜
8µ ⎝ ∂x ⎠
(2) GRADIENTE DE PRESIÓN EN FUNCIÓN DEL CAUDAL:
r r y=H / 2 y=H
1 ⎛ −∂p ⎞ 2 BH 3 ⎛ −∂p ⎞
Q=
∫∫ v ⋅ dA =
∫ u ⋅ Bdy =
∫ ⎟ ( H − 4y ) ⋅ Bdy = ... =
2
⎜ ⎜ ⎟
y =− H / 2 y=0 8µ ⎝ ∂x ⎠ 12µ ⎝ ∂x ⎠
Evidentemente, para que se tenga caudal, el gradiente de presión debe ser negativo. De la expresión anterior, se
obtiene el citado gradiente:
∂p 12µQ
=− [3]
∂x 3
BH
Numerícamente con los datos: µ=0,2 kg/ms; Q= 1 litro/segundo, B=1m; y H=0,010 m; se tiene que el
gradiente horizontal de presión es:
12µQ 12 ⋅ 0, 2 ⋅ 10−3
=− =− = −2400 Pa/m
BH 3 1 ⋅ 0,0103
(3) PERFIL DE TEMPERATURAS: T=T(y):
r
= −p(∇ ⋅ v ) + κ ⋅ ∇ 2 T + Φ
dû
Ec. ENERGÍA: ρ
dt
En este caso, no hay variación local de la energía interna, por ser el flujo estacionario; y la variación convectiva
de energía interna, es nula, porque sólo hay componente “x” de la velocidad, y la energía interna sólo varía con
la distancia vertical “y”:
∂û
estacionario ⇒ =0
∂t
dû ∂û ∂û ∂û ∂û dû
= +u +v +w = ... v = w = 0 ⇒ =0
dt ∂t ∂x ∂y ∂z dt
û = û ( y)
Como el flujo es adivergente, de la ecuación de energía se tiene que la laplaciana de la temperatura absoluta es:
2 −Φ
∇ T= [4]
κ
_________________________________________________________________________________________________________________
2 2 2 2
2 ∂ T ∂ T ∂ T d T
en el problema, como T=T(y), la laplaciana es: ∇ T = + + = ... =
2 2 2 2
∂x ∂y ∂z dy
y la función de disipación viscosa, con u=u(y), v=w=0, queda:
⎡ ⎛ ∂u ⎞ 2 ⎛ ∂v ⎞
2
⎛ ∂w ⎞ ⎛ ∂u ∂v ⎞ ⎛ ∂u ∂w ⎞ ⎛ ∂v ∂w ⎞ ⎤
2 2 2 2
⎛ du ⎞
2
Φ = µ ⎢2⎜ ⎟ + 2⎜⎜ ⎟⎟ + 2⎜ ⎟ + ⎜⎜ + ⎟⎟ + ⎜ + ⎟ + ⎜⎜ + ⎟⎟ ⎥ = ... = µ⎜⎜ ⎟⎟
⎢⎣ ⎝ ∂x ⎠ ⎝ ∂y ⎠ ⎝ ∂z ⎠ ⎝ ∂y ∂x ⎠ ⎝ ∂z ∂x ⎠ ⎝ ∂z ∂y ⎠ ⎥⎦ ⎝ dy ⎠
du 1 ⎛ ∂p ⎞
el gradiente de velocidad en la dirección vertical es: = ⎜ ⎟ ⋅ y ; con lo que la función de disipación
dy µ ⎝ ∂x ⎠
2 2
⎛ du ⎞ ⎛ 1 ∂p ⎞ 1 ⎛ ∂p ⎞
2
viscosa de RAYLEIGH, queda: Φ = µ⎜ ⎟ = µ⎜ ⋅ y ⎟ = ⎜ ⎟ y2 [5]

⎝ dy ⎠ ⎝ µ ∂x ⎠ µ ⎝ ∂x ⎠
con lo que la ecuación diferencial del campo de temperaturas [4], queda como:
1 ⎛ ∂p ⎞ 2
2
d 2T
=− ⎜ ⎟ y ;
dy 2 µκ ⎝ ∂x ⎠
2
−1 ⎛ ∂p ⎞ 4
cuya integración da: T= ⎜ ⎟ y + C3 ⋅ y + C4
12µκ ⎝ ∂x ⎠
las constantes de integración, se obtienen con las condiciones de contorno en las paredes:
&
Q ⎛ ∂T ⎞ ⎛ ∂T ⎞
La pared inferior es adiabática y por tanto = κ⎜ ⎟ =0 ⇒ ⎜ ∂y ⎟ =0
A ⎝ ∂y ⎠ y=− H /2 ⎝ ⎠ y=− H /2
La pared inferior es isoterma y por tanto: T(y=H/2) = Tisoterma
Operando, se obtiene el perfil de temperaturas, en cada posición horizontal:
1 ⎛ ∂p ⎞
2
isoterma 192µκ ⎜⎝ ∂x ⎟⎠ (
T=T + 5H 4 − 8H 3 ⋅ y − 16y 4 ) [6]
2
1 ⎛ ∂p ⎞
= T(y = − H / 2) = ... = T + H4
isoterma 24µκ ⎜⎝ ∂x ⎟⎠
La temperatura en la pared adiabática, es: T
adiabática
∂p
Numerícamente con los datos: Tisoterma =333,15K; µ=0,2 kg/ms; κ=0,6W/mK; H=0,010 m; =-2400Pa/m:
∂x
1 2 4
T = 333,15 + ( −2400) 0,010 = 333,17 K = 60,02 º C
adiabática 24 ⋅ 0, 2 ⋅ 0,6
_________________________________________________________________________________________________________________
(4) POTENCIA CALORIFÍCA QUE ATRAVIESA LA PARED ISOTERMA:
A partir de la ley de FOURIER de transmisión por conducción, la potencia calorifica que atraviesa una
determina sección, viene determinada por el coeficiente de conductividad térmica (κ) y por el gradiente de
r
temperaturas: Q& = − κ∇T ⋅ A . En el problema, la potencia calorifica elemental, que atraviesa un elemnto de
área (Bdx) de la pared isoterma, viene dada por:
⎛ ∂T r ∂T r ∂T r ⎞ r
k ⎟ ⋅ (− Bdx ⋅ j ) = κ
dT
isoterma = − κ⎜
&
dQ i+ j+ B ⋅ dx
⎝ ∂x ∂y ∂z ⎠ dy
−1 ⎛ ∂p ⎞
2
⎜ ⎟ ( H + 8y )
dT
en donde el gradiente vertical de temperatura es: = 3 3
dy 24µκ ⎝ ∂x ⎠
2
⎛ dT ⎞ − 1 ⎛ ∂p ⎞ 3
y en la pared isoterma (y=H/2), se tiene queese gradiente es: ⎜ ⎟ = ⎜ ⎟ H
⎝ dy ⎠ y =H 12µκ ⎝ ∂x ⎠
con lo que la potencia calorifica que atraviesa la pared isoterma queda como:
2
⎛ dT ⎞ − 1 ⎛ ∂p ⎞ 3
Q isoterma = κ⎜
& ⎟ BL = ⎜ ⎟ H BL [7]
⎝ dy ⎠ y=H 12µ ⎝ ∂x ⎠
∂p
Numerícamente con los datos: µ=0,2 kg/ms; H=0,010 m, B=L=1m; =-2400Pa/m:
∂x
2
− 1 ⎛ ∂p ⎞ −1
&
Q isoterma =
3
⎜ ⎟ H BL = (− 2400 )2 ⋅ 0,010 3 ⋅ 1 ⋅ 1 = −2,4W
12µ ⎝ ∂x ⎠ 12 ⋅ 0, 2
es decir, a traves de la pared isoterma se “evacua” una potencia calorífica de 2,4 W. El signo negativo da el
sentido de la transfrencia de calor: hacía el exterior.
Se puede comprobar que la potencia calorifica que se “evacua” ha sido generada por las irreversibilidades
debidas a los efectos viscosos:
irreversibilidades = ∫∫∫ Φ ⋅ dV
&
Q
en nuestro caso, la función de disipación viscosa de RAYLEIGH, viene dada por la ecuación [5]:
1 ⎛ ∂p ⎞ 2
2
Φ= ⎜ ⎟ y
µ ⎝ ∂x ⎠
Con lo que la potencia calorífica generada por los efectos viscosos es:
( H / 2) − ( − H / 2)3
3
1 ⎛ ∂p ⎞ 2 LB ⎛ ∂p ⎞
2 2
∫∫∫ ∫∫∫ ∫ ∫ ∫
L B H/2
&
Q irreversibilidades = Φ ⋅ dV = Φ ⋅ dx ⋅ dy ⋅ dz = dx dz ⎜ ⎟ y dy = ⎜ ⎟
0 0 − H / 2 µ ⎝ ∂x ⎠ µ ⎝ ∂x ⎠ 3
2
& 1 ⎛ ∂p ⎞ 3
Q = ⎜ ⎟ H LB
irreversibilidades 12µ ⎝ ∂x ⎠
que coincide (en modulo) con la ecuación [7] que da la potencia calorifica que atraviesa la pared isoterma.
_________________________________________________________________________________________________________________
ANEXO 1:
ECUACIONES DE ANÁLISIS DIFERENCIAL DE FLUIDOS.
E3.1. Cinemática.
E3.2. Operador nabla.
E3.2.1. Operador nabla en coordenadas cartesianas.
E3.2.2. Operador nabla en coordenadas cilíndricas.
E3.2.3. Operador nabla en coordenadas esféricas.
E3.3. Ecuación de conservación de masa: Ec. de continuidad.
E3.4. Ecuaciones de conservación de cantidad de movimiento.
E3.4.1. Ecuación de CAUCHY.
E3.4.2. Ecuación de NAVIER-STOKES.
E3.4.3. Ecuación de EULER.
E3.4.4. Ecuación de BERNOULLI.
E3.5. Ecuaciones de conservación de energía.
E3.5.1. Ecuación de la energía.
E3.5.2. Ecuación de la entalpía.
E3.5.3. Ecuación de la entropía.
E3.5.4. Función de disipación viscosa de RAYLEIGH.
E3.6. Ecuaciones de conservación para fluidos newtonianos: Ecuaciones de NAVIER-STOKES.
E3.6.1. Ecuaciones de Navier-Stokes en coordenadas cartesianas.
E3.6.2. Ecuaciones de Navier-Stokes en coordenadas cilíndricas.
E3.6.3. Ecuaciones de Navier-Stokes en coordenadas esféricas.
z
r
E3.1. CINEMÁTICA. v
COORDENADAS CARTESIANAS:
r dz
r r r r ds
Vector desplazamiento: ds = dx ⋅ i + dy ⋅ j + dz ⋅ k
r r dx
r ds dx r dy r dz r
Vector velocidad: v = = ⋅i+ ⋅ j+ ⋅k ds dy
dt dt dt dt r
ds
r r r r r y
v = u⋅ i + v⋅ j+ w⋅k ds
x
r r r r r r
dv ∂v r r ∂v ⎛ ∂v ∂v ∂v ⎞
Vector aceleración: = + v ⋅ ∇v = +⎜u + v + w ⎟ =
dt ∂t ∂t ⎝ ∂x ∂y ∂z ⎠
⎡ ∂u ⎛ ∂u ∂u ∂u ⎞ ⎤ r ⎡ ∂v ⎛ ∂v ∂v ∂v ⎞ ⎤ r ⎡ ∂w ⎛ ∂w ∂w ∂w ⎞ ⎤ r
= ⎢ + ⎜ u + v + w ⎟⎥ i + ⎢ + ⎜ u + v + w ⎟⎥ j + ⎢ +⎜u +v +w ⎟⎥ k [1]
⎣ ∂t ⎝ ∂x ∂y ∂z ⎠ ⎦ ⎣ ∂t ⎝ ∂x ∂y ∂z ⎠ ⎦ ⎣ ∂t ⎝ ∂x ∂y ∂z ⎠ ⎦
⎛ ∂u ∂v ∂w ⎞
⎜ ∂x ∂x ∂x ⎟
⎜ ⎟
r r ⎜ ∂u ∂v ∂w ⎟
Aceleración Convectiva: v ⋅ ∇v = ( u v w ) ⋅ ⎜ =
∂y ∂y ∂y ⎟
⎜ ⎟
⎜ ∂u ∂v ∂w ⎟
⎜ ⎟
⎝ ∂z ∂z ∂z ⎠
⎛ ∂u ∂u ∂u ⎞ r ⎛ ∂v ∂v ∂v ⎞ r ⎛ ∂w ∂w ∂w ⎞ r
= ⎜u + v + w ⎟⋅ i + ⎜u + v + w ⎟⋅ j + ⎜u +v +w ⎟⋅k =
⎝ ∂x ∂y ∂z ⎠ ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂x ∂y ∂z ⎠
r r r r r r r r r
=u
(
∂ u i + vj + wk
+v
) (
∂ u i + vj + wk
+w
)
∂ u i + vj + wk( r
∂v
=u +v +w
) r
∂v
r
∂v
∂x ∂y ∂z ∂x ∂y ∂z
_________________________________________________________________________________________________________________
z
rdθ r
COORDENADAS CILINDRÍCAS: v
dz dr
r r r r
Vector desplazamiento: ds = dr ⋅ u r + rdθ ⋅ u θ + dz ⋅ u z r
ds r
r
r ds dr r dθ r dz r θ uz
Vector velocidad: v= = ⋅ ur + r ⋅ uθ + ⋅ uz r
dt dt dt dt ur
r
r r r r uθ
v = v r ⋅ ur + vθ ⋅ uθ + vz ⋅ uz
dθ
velocidad de giro: ω =
dt
dθ
velocidad tangencial : v θ = r = r ⋅ω
dt
r r
dv ∂v r r ⎡ ∂v r ⎛ ∂v r v θ ∂v r v θ2 ∂v ⎞ ⎤ r
= + v ⋅ ∇v = ⎢ + ⎜ vr + − + v z r ⎟⎥ ⋅ u r +
dt ∂t ⎣ ∂t ⎝ ∂r r ∂θ r ∂z ⎠ ⎦
⎡ ∂v ⎛ ∂v v ∂v vv ∂v ⎞ ⎤ r
Vector aceleración: + ⎢ θ + ⎜ v r θ + θ θ + r θ + v z θ ⎟⎥ ⋅ u θ + [2]
⎣ ∂t ⎝ ∂r r ∂θ r ∂z ⎠ ⎦
⎡ ∂v ⎛ ∂v v ∂v ∂v ⎞ ⎤ r
+ ⎢ z + ⎜ v r z + θ z + v z z ⎟⎥ ⋅ u z
⎣ ∂t ⎝ ∂r r ∂θ ∂z ⎠ ⎦
r
u
r
COORDENADAS ESFÉRICAS:
r
v
r
u
θ
Vector desplazamiento:
dϕ r
ϕ u
r r r r ϕ
ds = dr ⋅ u r + rdθ ⋅ u θ + rsenϕdϕ ⋅ u ϕ r
Vector velocidad: θ
r.senϕ
r r r r dθ
v = v r ⋅ ur + v θ ⋅ uθ + vϕ ⋅ uϕ
Vector aceleración:
r r
dv ∂v r r ⎡ ∂v r ⎛ ∂v r v θ ∂v r v ∂v r v θ + v ϕ ⎞⎤ r
2 2
= + v ⋅ ∇v = ⎢ + ⎜ vr + + θ − ⎟⎟ ⎥ ⋅ u r +
dt ∂t ⎢⎣ ∂t ⎜⎝ ∂r r ∂θ rsenθ ∂ϕ r ⎠ ⎥⎦
⎡ ∂v ⎛ ∂v v ∂v v ϕ ∂v θ v r v θ v ϕ2 ⎞ ⎤ r
+ ⎢ θ + ⎜ vr θ + θ θ +
⎜
+ − ⎟⎥ ⋅ u θ + [3]
⎣⎢ ∂t ⎝ ∂r r ∂θ rsenθ ∂ϕ r rtgθ ⎟⎠ ⎦⎥
⎡ ∂v ϕ ⎛ ∂v ϕ v θ ∂v ϕ v ϕ ∂v ϕ v r v ϕ v θ v ϕ ⎞⎤ r
+⎢ + ⎜ vr + + + + ⎟⎥ ⋅ uϕ
⎣ ∂t ⎝ ∂r r ∂θ rsenθ ∂ϕ r rtgθ ⎠⎦
_________________________________________________________________________________________________________________
E3.2. OPERADOR NABLA “ ∇ ”
OPERADOR NOTACIÓN SIMBÓLICA TIPO
GRADIENTE DE UN CAMPO ESCALAR ∇T vector

r
GRADIENTE DE UN CAMPO VECTORIAL ∇v tensor
r
DIVERGENCIA DE UN CAMPO VECTORIAL ∇⋅v escalar
DIVERGENCIA DE UN CAMPO TENSORIAL ∇ ⋅τ vector
r
ROTACIONAL DE UN CAMPO VECTORIAL ∇×v vector
LAPLACIANA DE UN CAMPO ESCALAR ∇ ⋅ (∇ T ) = ∆T = ∇ 2 T escalar

r r
LAPLACIANA DE UN CAMPO VECTORIAL ∇ ⋅ ∇v = ∇ 2 v vector
r r r r
ACELERACIÓN CONVECTIVA v ⋅ ∇v = (v ⋅ ∇ ) ⋅ v vector
E3.2.1. OPERADOR NABLA EN COORDENADAS CARTESIANAS.
GRADIENTE DE PRESIÓN: ∇p = ⋅i+ ⋅ j+ ⋅k [4]
∂x ∂y ∂z
⎛ ∂u ∂v ∂w ⎞
⎜ ∂x ∂x ∂x ⎟⎟
⎜
r ⎜ ∂u ∂v ∂w ⎟
GRADIENTE DE VELOCIDAD: ∇v = ⎜ [5]
∂y ∂y ∂y ⎟
⎜ ⎟
∂
⎜ u ∂v ∂w ⎟
⎜ ⎟
⎝ ∂z ∂z ∂z ⎠
r ∂u ∂v ∂w
DIVERGENCIA DE VELOCIDAD: ∇⋅v = + + [6]
∂x ∂y ∂z
DIVERGENCIA DEL TENSOR DE TENSIONES VISCOSAS:
⎛ τ xx τ xy τ xz ⎞
⎛ ∂ ∂ ∂ ⎞ ⎜ ⎟
∇⋅τ = ⎜ ⎟ ⋅ τ yx τ yy τ yz ⎟ =
⎝ ∂x ∂y ∂z ⎠ ⎜⎜
⎝ τ zx τzy τ zz ⎟⎠ [7]
⎛ ∂τ ∂τ yx ∂τzx ⎞ r ⎛ ∂τ xy ∂τ yy ∂τ zy ⎞ r ⎛ ∂τ xz ∂τ yz ∂τzz ⎞ r
= ⎜ xx + + ⎟⋅ i + ⎜ + + ⎟⋅ j + ⎜ + + ⎟⋅k
⎝ ∂x ∂y ∂z ⎠ ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂x ∂y ∂z ⎠
r r r r
ROTACIONAL DE VELOCIDAD: ∇ × v = Ω = VORTICIDAD = 2ω ω =velocidad de giro
r ⎛ ∂w ∂v ⎞ r ⎛ ∂u ∂w ⎞ r ⎛ ∂v ∂u ⎞ r
∇×v = ⎜ − ⎟⋅ i + ⎜ − ⎟⋅ j + ⎜ − ⎟⋅k [8]
⎝ ∂y ∂z ⎠ ⎝ ∂z ∂z ⎠ ⎝ ∂x ∂y ⎠
∂2T ∂2T ∂2T

LAPLACIANA DE TEMPERATURA: ∇ ⋅ (∇ T ) = + + [9]
∂x 2 ∂y 2 ∂z 2
_________________________________________________________________________________________________________________
LAPLACIANA DE VELOCIDAD:
r ⎛ ∂2 u ∂2 u ∂2 u ⎞ r ⎛ ∂2 v ∂2 v ∂2 v ⎞ r ⎛ ∂2 w ∂2 w ∂2 w ⎞ r
∇ ⋅ ∇v = ⎜ 2 + 2 + 2 ⎟ ⋅ i + ⎜ 2 + 2 + 2 ⎟ ⋅ j + ⎜ 2 + 2 + 2 ⎟ ⋅ k [10]
⎝ ∂x ∂y ∂z ⎠ ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂x ∂y ∂z ⎠
TENSOR DE VELOCIDAD DE DEFORMACIÓN:

⎛ 1 ⎛ ∂u ∂v ⎞ 1 ⎛ ∂u ∂w ⎞ ⎞⎟
⎜ ∂u
⎜ ⎜ + ⎟ ⎜ + ⎟⎟
⎜
∂x 2 ⎝ ∂y ∂x ⎠ 2 ⎝ ∂z ∂x ⎠ ⎟
= ⎜ 1 ⎛ ∂v ∂u ⎞ ∂v 1 ⎛ ∂v ∂w ⎞ ⎟
ε& = ⎜ ⎜ + ⎟ ⎜ + ⎟ ⎟⎟ [11]
⎜ 2 ∂x ∂y
⎜ ⎝ ⎠ ∂y 2 ⎝ ∂z ∂y ⎠
⎟
⎜ 1 ⎛ ∂w ∂u ⎞ 1 ⎛ ∂w ∂v ⎞ ∂w ⎟
⎜ ⎜ + ⎟ ⎜ + ⎟ ⎟
⎜ 2 ⎝ ∂x ∂z ⎠ 2 ⎝ ∂y ∂z ⎠ ∂z ⎟
⎝ ⎠
E3.2.2. OPERADOR NABLA EN COORDENADAS CILÍNDRICAS.
∂p r 1 ∂p r ∂p r
GRADIENTE DE PRESIÓN: ∇p = ⋅ ur + ⋅ uθ + ⋅ uz [12]
∂r r ∂θ ∂z
⎛ ∂v r ∂v θ ∂v z ⎞
⎜ ∂r ∂r ∂r ⎟
⎜ ⎟
r ⎜ 1 ⎛ ∂v r ⎞ 1 ⎛ ∂v θ ⎞ 1 ∂v z ⎟
GRADIENTE DE VELOCIDAD: ∇v = ⎜ ⎜ − vθ ⎟ ⎜ + vr ⎟ ⎟ [13]
r ∂θ ⎠ r ⎝ ∂θ ⎠ r ∂θ ⎟
⎜ ⎝
⎜ ∂v r ∂v θ ∂v z ⎟
⎜ ⎟
⎝ ∂z ∂z ∂z ⎠
r 1 ∂ ( r v r ) 1 ∂v θ ∂v z
DIVERGENCIA DE VELOCIDAD: ∇⋅v = + + [14]
r ∂r r ∂θ ∂z
⎛ 1 ∂ ( rτ rr ) 1 ∂τ θr ∂τzr τ θθ ⎞ r
∇⋅τ = ⎜ + + − ⎟ ⋅ ur +
⎝ r ∂r r ∂θ ∂z r ⎠
⎛ 1 ∂ ( rτ rθ ) 1 ∂τ θθ ∂τ zθ τ rθ ⎞ r
+⎜ + + + ⎟ ⋅ uθ + [15]
⎝ r ∂r r ∂θ ∂z r ⎠
⎛ 1 ∂ ( rτ rz ) 1 ∂τ θz ∂τ zz ⎞ r
+⎜ + + ⎟ ⋅ uz
⎝ r ∂r r ∂θ ∂z ⎠
ROTACIONAL DE VELOCIDAD:
r ⎛ 1 ∂v z ∂v θ ⎞ r ⎛ 1 ∂ ( rv θ ) 1 ∂v r ⎞ r
⎞ r ⎛ ∂v r ∂v z
∇×v = ⎜ − ⎟ ⋅ u r + ⎜ ∂z − ∂z
⎟ ⋅ u θ + ⎜ r ∂r − r ∂θ ⎟ ⋅ u z [16]
⎝ r ∂θ ∂z ⎠
⎠ ⎝
⎝ ⎠
∂ 2
T 1 ∂T 1 ∂ 2 T ∂ 2 T
LAPLACIANA DE TEMPERATURA: ∇ ⋅ ( ∇T ) = ∇ 2 T = 2 + + + [17]
∂r r ∂r r 2 ∂θ2 ∂z 2
_________________________________________________________________________________________________________________
LAPLACIANA DE VELOCIDAD:
r r ⎡ ∂ 2 v 1 ∂v r 1 ∂ 2 v r ∂ 2 v r 1 ⎛ ∂v ⎞ ⎤ r
∇ ⋅ ∇v = ∇ 2 v = ⎢ 2 r + + 2 + − 2 ⎜ v r + 2 θ ⎟⎥ ⋅ u r +
⎣ ∂r r ∂r r ∂θ 2
∂z 2
r ⎝ ∂θ ⎠ ⎦
⎡ ∂2 v 1 ∂v θ 1 ∂ 2 v θ ∂ 2 v θ 1 ⎛ ∂v ⎞ ⎤ r
+ ⎢ 2θ + + 2 + − 2 ⎜ v θ − 2 r ⎟⎥ ⋅ uθ + [18]
⎣ ∂r r ∂r r ∂θ 2
∂z 2
r ⎝ ∂θ ⎠ ⎦
⎡ ∂2 v 1 ∂v z 1 ∂ 2 v z ∂ 2 v z ⎤ r
+ ⎢ 2z + + + ⎥ ⋅ uz
⎣ ∂r r ∂r r 2 ∂θ2 ∂z 2 ⎦
⎛ ∂v 1 ⎛ ∂ ( v θ r ) 1 ∂v r ⎞ 1 ⎛ ∂v z ∂v r ⎞ ⎟
⎞
⎜ r
⎜ ⎜r + ⎟ 2 ⎝⎜ ∂r
+ ⎟
∂z ⎠ ⎟
⎜
∂r 2⎝ ∂r r ∂θ ⎠ ⎟
= ⎜ 1 ⎛ ∂ ( v θ r ) 1 ∂v r ⎞ ⎛1 ∂v
θ + r⎟
v ⎞ 1 ⎛ ∂v θ 1 ∂v z ⎞ ⎟
ε& = ⎜ ⎜ r + ⎟ ⎜⎜ ⎟ +
2 ⎜⎝ ∂z r ∂θ ⎟⎠ ⎟
⎟ [19]
⎜ 2⎝ ∂r r ∂θ ⎠ ⎝ r ∂θ r ⎠
⎜ ∂v ⎟
⎜ 1 ⎛ ∂v z ∂v r ⎞ 1 ⎛ ∂v θ 1 ∂v z ⎞ z ⎟
+ +
⎜ 2 ⎜⎝ ∂r ∂z ⎠
⎟ 2 ⎜⎝ ∂z r ∂θ ⎟⎠ ∂z ⎟
⎜ ⎟
⎝ ⎠
E3.2.3. OPERADOR NABLA EN COORDENADAS ESFÉRICAS.
∂p r 1 ∂p r 1 ∂p r
GRADIENTE DE PRESIÓN: ∇p = ⋅ ur + ⋅ uθ + ⋅ uϕ [20]
∂r r ∂θ rsenθ ∂ϕ
r 1 ∂ (r vr ) 1 ∂ ( senθ v θ ) 1 ∂v ϕ
2
DIVERGENCIA DE VELOCIDAD: ∇⋅v = 2 + + [21]

r ∂r r senθ ∂θ r senθ ∂ϕ
⎛ 1 ⎡ ∂ ( r 2 senθ τ rr ) ∂ ( rsenθ τ ) ∂ ( rτϕr ) ⎤ τ + τ ⎞

⎜
∇⋅τ = 2 ⎢ + θr
+ ⎥−
θθ ϕϕ
⎟ ⋅ ur r +
⎜ r senθ ⎢ ∂r ∂θ ∂ϕ ⎥ r ⎟
⎝ ⎣ ⎦ ⎠
⎛ 1 ⎡ ∂ ( r 2 senθ τ rθ ) ∂ ( rsenθ τ ) ∂ ( rτϕθ ) ⎤ τ − τ / tgθ ⎞ r
+⎜ 2 ⎟ ⋅ uθ +
rθ ϕϕ
⎢ + θθ
+ ⎥+ [22]
⎜ r senθ ⎢ ∂r ∂θ ∂ϕ ⎥ r ⎟
⎝ ⎣ ⎦ ⎠
⎛ 1 ⎡ ∂ ( r 2 senθ τ rϕ ) ∂ ( rsenθ τθϕ ) ∂ ( rτϕϕ ) ⎤ τ + τ / tgθ ⎞ r
+⎜ 2 ⎟ ⋅ uϕ +
rϕ θϕ
⎢ + + ⎥+
⎜ r senθ ⎢ ∂r ∂θ ∂ϕ ⎥ r ⎟
⎝ ⎣ ⎦ ⎠
ROTACIONAL DE VELOCIDAD: [23]
r 1 ⎛ ∂ ( rsenθ v ϕ ) ∂ ( rv θ ) ⎞ r 1 ⎛ ∂v r ∂ ( rsenθ v ϕ ) ⎞ r 1 ⎛ ∂ ( rv θ ) ∂v r ⎞ r
∇×v = ⎜ − ⎟ ⋅ ur + ⎜ − ⎟ ⋅ uθ + ⎜ − ⎟ ⋅ uϕ
r senθ ⎜
2
∂θ ∂ϕ ⎟ r senθ ⎜ ∂ϕ ∂r ⎟ r ⎝ ∂r ∂θ ⎠
⎝ ⎠ ⎝ ⎠
LAPLACIANA DE TEMPERATURA:
1 ∂ ⎛ 2 ∂T ⎞ 1 ∂⎛ ∂T ⎞ 1 ∂2T
∇2 T = ⎜r ⎟+ 2 ⎜ senθ ⎟+ 2 2 [24]
r ∂r ⎝ ∂r ⎠ r senθ ∂r ⎝
2
∂θ ⎠ r sen θ ∂ϕ2
_________________________________________________________________________________________________________________
r ⎡ 2⎛ ∂v v 1 ∂v ϕ ⎞⎤ r
∇ 2 v = ⎢∇ 2 v r − 2 ⎜ v r + θ + θ + ⎟⎥ ⋅ u r +
⎣ r ⎝ ∂θ tgθ senθ ∂ϕ ⎠⎦
⎡ 2 ∂v 1 ⎛ ∂v ϕ ⎞ ⎤ r
LAPLACIANA DE VELOCIDAD: + ⎢∇ 2 v θ + 2 r − 2 2 ⎜ v θ + 2 cos θ ⎟⎥ ⋅ uθ + [25]
⎣ r ∂θ r sen θ ⎝ ∂ϕ ⎠ ⎦
⎡ 1 ⎛ ∂v ∂v ⎞ ⎤ r
+ ⎢∇ 2 v ϕ − 2 2 ⎜ v ϕ − 2senθ r − 2 cos θ θ ⎟ ⎥ ⋅ u ϕ
⎣ r sen θ ⎝ ∂ϕ ∂ϕ ⎠ ⎦
⎛ ⎞
⎜ ∂v r 1 ⎛ ∂ ( v θ / r ) 1 ∂v r ⎞ 1
⎛ 1 ∂v ∂ ⎛ v ϕ ⎞⎞ ⎟
⎜ ⎜r + ⎟ ⎜ r
+r ⎜ ⎟⎟ ⎟⎟ ⎟
⎜ ∂r 2⎝ ∂r r ∂θ ⎠
2 ⎜ rsenθ ∂ϕ ∂r ⎜ r ⎟
⎜ ⎝ ⎝ ⎠⎠ ⎟
= ⎜⎜ 1 ⎛ ∂ ( v θ / r ) 1 ∂v r ⎞ 1 ∂v r v ⎛
1 senθ ∂ ⎛ ϕ ⎞
v 1 ∂v θ ⎞ ⎟
⎟
ε& = ⎜ ⎜ r ⎜
⎟ 2 ⎜ r ∂θ ⎜ ⎜ senθ ⎟⎟ rsenθ ∂ϕ ⎟⎟ ⎟⎟
+ + +
r
⎜ 2⎝ ∂r r ∂θ ⎠ r ∂θ
⎜
r
⎝ ⎝ ⎠ ⎠⎟
⎜ ⎛
∂ ⎛ vϕ ⎞⎞ ⎛
1 senθ ∂ ⎛ ϕ ⎞
⎞ ∂ ⎟
⎜1 1 ∂v r v 1 ∂v θ 1 v v u
⎟
⎜ 2⎜ +r ⎜ ⎟⎟ ⎜ ⎜ ⎟+ ⎟
ϕ θ
+ +
r
∂r ⎜ r ⎟ ⎟ r ∂θ ⎜ senθ ⎟ rsenθ ∂ϕ
⎟
⎜⎜ ⎜ rsenθ ∂ϕ
2⎜ ⎟ rsenθ ∂ϕ r rtgθ
⎝ ⎝ ⎝ ⎠⎠ ⎝ ⎝ ⎠ ⎠ ⎟⎟
⎠
[26]
E3.3. Conservación de masa: Ecuación de CONTINUIDAD.
∂ρ r
+ ∇ ⋅ (ρ v ) = 0 [27]
∂t
dρ r
+ρ ∇⋅v = 0 [28]
dt
COORDENADAS CARTESIANAS:
∂ρ ⎡ ∂ ( ρ u ) ∂ (ρ v ) ∂ (ρ w ) ⎤
+⎢ + + ⎥=0 [29]
∂t ⎣ ∂x ∂y ∂z ⎦
COORDENADAS CILINDRÍCAS:
∂ρ ⎡ 1 ∂ ( ρ r v r ) 1 ∂ ( ρ v θ ) ∂ ( ρ v z ) ⎤
+⎢ + + ⎥=0 [30]
∂t ⎣ r ∂r r ∂θ ∂z ⎦
COORDENADAS ESFÉRICAS:
∂ρ ⎡ 1 ∂ ( ρ r v r ) 1 ∂ ( ρ senθ v θ ) 1 ∂ (ρ v ϕ ) ⎤
2
+⎢ 2 + + ⎥=0 [31]
∂t ⎢ r ∂r r senθ ∂θ r senθ ∂ϕ ⎥
⎣ ⎦
_________________________________________________________________________________________________________________
E3.4. Conservación de cantidad de movimiento: Ecuaciones de MOVIMIENTO.

E3.4.1. ECUACIÓN DE MOVIMIENTO DE CAUCHY
r
r dv
ρ g − ∇p + ∇ ⋅ τ = ρ [32]
dt
E3.4.2.ECUACIÓN DE MOVIMIENTO PARA FLUIDOS NEWTONIANOS: Ec. DE NAVIER-STOKES

r
r r r dv
ρ g − ∇p + ( λ + µ ) ∇ ( ∇ ⋅ v ) + µ ∇ 2 v = ρ [33]
dt
2 µ
condición de STOKES: viscosidad bruta nula ⇒ µB = λ + µ = 0 ⇒ λ + µ =
3 3
r
r µ r 2r dv
ρ g − ∇p + ∇ (∇ ⋅ v ) + µ ∇ v = ρ [34]
3 dt
r
r r r dv
flujo incompresible: adivergente ⇒ ∇ ⋅ v = 0 ρ g − ∇p + µ ∇2 v = ρ [35]
dt
ECUACIONES DE NAVIER-STOKES PARA FLUJO INCOMPRESIBLE EN CARTESIANAS:

∂p ⎛ ∂2 u ∂2 u ∂2 u ⎞ ⎛ ∂u ∂u ∂u ∂u ⎞
ρ gx − + µ ⎜ 2 + 2 + 2 ⎟ = ρ⎜ + u +v +w ⎟ [36]
∂x ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂t ∂x ∂y ∂z ⎠
∂p ⎛ ∂2 v ∂2 v ∂2 v ⎞ ⎛ ∂v ∂v ∂v vu ⎞
ρ gy − + µ ⎜ 2 + 2 + 2 ⎟ = ρ⎜ + u +v +w ⎟ [37]
∂y ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂t ∂x ∂y ∂z ⎠
∂p ⎛∂ w ∂ w ∂ w⎞
2 2 2
⎛ ∂w ∂w ∂w ∂w ⎞
ρ gz − + µ ⎜ 2 + 2 + 2 ⎟ = ρ⎜ +u +v +w ⎟ [38]
∂z ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂t ∂x ∂y ∂z ⎠
ECUACIONES DE NAVIER-STOKES PARA FLUJO INCOMPRESIBLE EN CILÍNDRICAS:

∂p ⎛ ∂ ⎛ 1 ∂ ( r v r ) ⎞ 1 ∂ 2 v r 2 ∂v θ ∂ 2 v r ⎞
ρ gr − +µ ⎜ ⎜ ⎟+ 2 − 2 + 2 ⎟=
∂r ⎜ ∂r r ∂r r ∂θ ∂z ⎟⎠
⎠ r ∂θ
2
⎝ ⎝
Componente radial “r”: [39]
⎛ ∂v ∂v v ∂v v2 ∂v ⎞
= ρ ⎜ r + vr r + θ r − θ + vz r ⎟
⎝ ∂t ∂r r ∂θ r ∂z ⎠
1 ∂p ⎛ ∂ ⎛ 1 ∂ ( r v θ ) ⎞ 1 ∂ 2 v θ 2 ∂v r ∂ 2 v θ ⎞
ρ gθ − +µ ⎜⎜ ⎜ ⎟+ 2 + 2 + ⎟⎟ =
r ∂θ ⎝ ∂r ⎝ r ∂r ⎠ r ∂θ
2
r ∂θ ∂z 2 ⎠
Componente tangencial “θ”: [40]
⎛ ∂v ∂v v ∂v v v ∂v ⎞
= ρ ⎜ θ + vr θ + θ θ + θ r + vz θ ⎟
⎝ ∂t ∂r r ∂θ r ∂z ⎠
∂p ⎛ ∂ ⎛ 1 ∂ ( r vz ) ⎞ 1 ∂2 vz ∂2 vz ⎞
ρ gθ − +µ ⎜ ⎜ ⎟+ 2 + ⎟=
∂z ⎜ ∂r r ∂r ∂z 2 ⎟⎠
⎠ r ∂θ
2
Componente axial “z”: ⎝ ⎝ [41]

⎛ ∂v ∂v v ∂v ∂v ⎞
= ρ ⎜ z + vr z + θ z + vz z ⎟
⎝ ∂t ∂r r ∂θ ∂z ⎠
ECUACIONES DE NAVIER-STOKES PARA FLUJO INCOMPRESIBLE EN ESFÉRICAS:
Componente radial “r”:
⎡ ⎛ 2 ∂v r ⎞ ⎛ ⎛ ∂v ⎞⎞ ⎤
⎢ ∂⎜⎜ r ⎟⎟ ⎜ ∂⎜⎜ senθ r ⎟⎟ ⎟ ⎥
∂p 1 ⎝ ∂r ⎠ 1 ⎜ ⎝ ∂θ ⎠⎟ 1 ∂2vr ⎛ ∂v 1 ∂vϕ ⎞⎥
+ µ⎢ 2
2 v
ρg r − + 2 ⎜ ⎟ + senθ − 2 ⎜⎜ v r + θ + θ + ⎟⎟ =
∂r ⎢ ∂r r senθ ⎜ ∂θ ∂ϕ 2 ∂θ tgθ senθ ∂ϕ ⎠⎥
⎢
r ⎟ r ⎝ ⎥
⎜ ⎟
⎣⎢ ⎝ ⎠ ⎦⎥
⎛ ∂v ∂v r 1 ∂v r ∂v r v θ + v ϕ
2 2 ⎞
ρ⎜ r
+ vr + vθ +
1
vϕ − ⎟
⎜ ∂t ∂r r ∂θ rsenθ ∂ϕ r ⎟
⎝ ⎠
[42]
_________________________________________________________________________________________________________________
Componente angular “θ”:

⎡ ⎛ 2 ∂v θ ⎞ ⎛ ⎛ ∂v ⎞⎞ ⎤
⎢ ∂⎜⎜ r ⎟⎟ ⎜ ∂⎜⎜ senθ θ ⎟⎟ ⎟ ⎥
1 ∂p ⎝ ∂ ⎠ ⎜ ⎝ ∂θ ⎠⎟ 1 ∂ v θ 2 ∂v θ
2
⎛ ∂v ⎞⎥
+ µ⎢ 2
1 r 1 1
ρg θ − + 2 ⎜ +
⎟ senθ − 2 − 2 ⎜ v θ + 2 cos θ ϕ ⎟ =
r ∂θ ⎢ ∂r r senθ ⎜ ∂θ ∂ϕ 2
r ∂θ r sen θ ⎝
2 ⎜ ∂ϕ ⎟⎥
⎢
r ⎟ ⎠⎥
⎢⎣ ⎜ ⎟ ⎥⎦
⎝ ⎠
⎛ ∂v ∂v θ 1 ∂v ∂v v r v θ − v ϕ2 / tgθ ⎞⎟
ρ⎜ θ + v r
1
+ vθ θ + vϕ θ −
⎝ ⎠
[43]
Componente angular “ϕ”:
⎡ ⎛ ∂vϕ ⎞ ⎛ ⎛ ∂v ⎞ ⎞ ⎤
⎢ ∂⎜⎜ r 2 ⎟⎟ ⎜ ∂⎜ senθ ϕ ⎟ ⎟ ⎥
1 ∂p ∂r ⎠ 1 ⎜ ⎝ ⎜ ⎟
∂θ ⎠ ⎟ 1 ∂ vϕ
2
⎛ ∂vϕ ∂v r ⎞⎥
+ µ ⎢⎢ 2 ⎝
1 1
ρg ϕ − + 2 ⎜ ⎟+ − 2 2 ⎜⎜ vϕ − 2 cos θ − 2senθ ⎟ =
rsenθ ∂ϕ r ∂r r senθ ⎜ ∂θ ⎟ senθ ∂ϕ
2
r sen θ ⎝ ∂ϕ ∂ϕ ⎟⎠⎥
⎢ ⎥
⎢⎣ ⎜ ⎟ ⎥⎦
⎝ ⎠
⎛ ∂vϕ ∂vϕ 1 ∂vϕ 1 ∂vϕ v r vϕ + vϕ vϕ / tgθ ⎞⎟
ρ⎜ + vr + vθ + vϕ + [44]
⎝ ⎠
r
r dv
E3.4.3.ECUACIÓN DE EULER. ρg − ∇p = ρ
dt
ECUACIONES DE EULER PARA FLUJO INCOMPRESIBLE EN CARTESIANAS:
∂p ⎛ ∂u ∂u ∂u ∂u ⎞
ρ gx − = ρ⎜ + u +v +w ⎟ [45]
∂x ⎝ ∂t ∂x ∂y ∂z ⎠
∂p ⎛ ∂v ∂v ∂v vu ⎞
ρ gy −
= ρ⎜ + u +v +w ⎟ [46]
∂y ⎝ ∂t ∂x ∂y ∂z ⎠
∂p ⎛ ∂w ∂w ∂w ∂w ⎞
ρ gz − = ρ⎜ +u +v +w ⎟ [47]
∂z ⎝ ∂t ∂x ∂y ∂z ⎠
ECUACIONES DE EULER PARA FLUJO INCOMPRESIBLE EN CILÍNDRICAS:
∂p ⎛ ∂v ∂v v ∂v v2 ∂v ⎞
ρ gr − = ρ ⎜ r + vr r + θ r − θ + vz r ⎟ [48]
∂r ⎝ ∂t ∂r r ∂θ r ∂z ⎠
1 ∂p ⎛ ∂v ∂v v ∂v v v ∂v ⎞
ρ gθ − = ρ ⎜ θ + vr θ + θ θ + θ r + vz θ ⎟ [49]
r ∂θ ⎝ ∂t ∂r r ∂θ r ∂z ⎠
∂p ⎛ ∂v ∂v v ∂v ∂v ⎞
ρ gθ − = ρ ⎜ z + vr z + θ z + vz z ⎟ [50]
∂z ⎝ ∂t ∂r r ∂θ ∂z ⎠
ECUACIONES DE EULER PARA FLUJO INCOMPRESIBLE EN ESFÉRICAS:
∂p ⎛ ∂v ∂v ∂v ∂v v θ2 + v ϕ2 ⎞
Componente radial “r”: ρg r − = ρ⎜⎜ r + v r r + v θ r + ⎟
1 1
vϕ r −
⎟
[51]
∂r ∂t ∂r r ∂θ rsenθ ∂ϕ r
⎝ ⎠
Componente angular “θ”:
1 ∂p ⎛ ∂v ∂v θ 1 ∂v θ ∂v θ v r v θ − v ϕ / tgθ ⎞⎟
2
ρg θ − =ρ⎜ θ
+ vr + vθ +
1
vϕ − [52]
r ∂θ ⎜ ∂t ∂r r ∂θ rsenθ ∂ϕ r ⎟
⎝ ⎠
Componente angular “ϕ”:
1 ∂p ⎛ ∂vϕ ∂vϕ 1 ∂v ϕ 1 ∂vϕ v r vϕ + v ϕ vϕ / tgθ ⎞⎟
ρg ϕ − = ρ⎜ + vr + vθ + vϕ + [53]
rsenθ ∂ϕ ⎜ ∂t ∂r r ∂θ rsenθ ∂ϕ r ⎟
⎝ ⎠
∂v v 2 − v12
+ g(z 2 − z1 ) = 0
2 2 dp
E3.4.4.ECUACIÓN DE BERNOULLI. ∫1 ∂t
ds + ∫
1 ρ
+ 2
2
[54]
FLUJO ESTACIONARIO INCOMPRESIBLE p + 12 ρv 2 + ρgz = cte LC [55]

_________________________________________________________________________________________________________________
E3.5. ECUACIÓN DE CONSERVACIÓN DE ENERGÍA.

E3.5.1.ECUACIÓN DE ENERGÍA INTERNA.
duˆ r
ρ = − p ( ∇ ⋅ v ) + κ ⋅ ∇2 T + Φ [56]
dt
⎛ ∂û ⎞ dT r
c v = ⎜ ⎟ = cte ⇒ duˆ = c v dT = − p ( ∇ ⋅ v ) + κ ⋅ ∇2 T + Φ
ρc v [57]
⎝ ∂t ⎠ dt
∂T
flujos muy lentos: ecuación de conducción en fluidos en reposo: ρc v = κ∇2 T [58]
∂t
E3.5.2.ECUACIÓN DE ENTALPÍA.
dh dp
ρ = + κ ⋅ ∇2 T + Φ [59]
dt dt
⎛ ∂h ⎞ dT dp
c p = ⎜ ⎟ = cte ⇒ dh = c p dT ρc p = + κ ⋅ ∇2 T + Φ [60]
⎝ ∂t ⎠ dt dt
E3.5.3.ECUACIÓN DE ENTROPÍA.
ds
ρT = κ ⋅ ∇2 T + Φ [61]
dt
E3.5.4.ECUACIÓN DE DISIPACIÓN DE ENERGÍA DE RAYLEIGH.

r
( ) ( )
r
Φ = ∇⋅ v⋅τ − v⋅ ∇⋅τ [62]
EN FLUJO INCOMPRESIBLE:
ECUACIÓN DE DISIPACIÓN VISCOSA DE RAYLEIGH, EN COORDENADAS CARTESIANAS:
⎡ ⎛ ∂u ⎞ 2 ⎛ ∂v ⎞
2
⎛ ∂w ⎞ ⎛ ∂u ∂v ⎞ ⎛ ∂u ∂w ⎞ ⎛ ∂v ∂w ⎞ ⎤
2 2 2 2
Φ = µ ⎢2 ⎜ ⎟ + 2 ⎜ ⎟ + 2 ⎜ ⎟ +⎜ + ⎟ +⎜ + ⎟ +⎜ + ⎟ ⎥ [63]
⎢⎣ ⎝ ∂x ⎠ ⎝ ∂y ⎠ ⎝ ∂z ⎠ ⎝ ∂y ∂x ⎠ ⎝ ∂z ∂x ⎠ ⎝ ∂z ∂y ⎠ ⎥⎦
ECUACIÓN DE DISIPACIÓN VISCOSA DE RAYLEIGH, EN COORDENADAS CILINDRÍCAS: [64]
( ) ⎟⎞
2⎤
⎡ ∂v 2 2 2 2 2 ⎛
2 ⎛ ∂v θ ⎛ ∂v ⎛ ∂v 1 ∂v z ⎛ ∂v ∂v ⎞ ∂ v /r
⎛ ⎞ ⎞ ⎞ ⎞ ⎜ 1 ∂v r ⎥
Φ = µ ⎢⎢ 2 ⎜ r ⎟ + ⎜ + v ⎟ + 2⎜ z ⎟ + ⎜ θ+ ⎟ +⎜ z + r ⎟ +⎜ +r θ
⎟ ⎥
⎜ ⎟ 2⎜ r⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ∂r ⎟
⎢⎣ ⎝ ∂r ⎠ r ⎝ ∂θ ⎠ ⎝ ∂z ⎠ ⎝ ∂z r ∂θ ⎠ ⎝ ∂z ⎠ ⎜ r ∂θ ∂y ⎟ ⎥
⎝ ⎠ ⎥
⎦
ECUACIÓN DE DISIPACIÓN VISCOSA DE RAYLEIGH, EN COORDENADAS ESFÉRICAS:
⎡
2 ⎛⎜ 1 ∂v ϕ
2
⎛ ∂v 2 ⎛⎜ ∂v θ v ⎞
2
⎞
2
⎞
Φ = µ ⎢⎢ 2 ⎜ r ⎟ + ⎟
+v + 2 +v + θ ⎟ +
⎜ ⎟ ⎜ ⎟ ⎜ r tgθ ⎟⎟
⎢ ⎝ ∂r r ⎝ ∂θ r ⎜ senθ ∂ϕ
2 r
⎠ ⎠ ⎝ ⎠
⎣
⎤
( vθ / r ) + 1 ∂v r
2 2
⎛ ⎞ ⎛ ⎞
∂ ⎛⎜ v / r ⎞⎟ ⎟ ⎜ ∂ ⎛ v / senθ ⎞⎟
2
⎛ ∂ ⎞
⎜ ⎟ ⎜ 1 ∂v
r ⎝ ϕ ⎠ senθ ⎜⎝ ϕ ⎠ 1 ∂v θ ⎟ ⎥⎥
+ ⎜r ⎟ +⎜ +r ⎟ +⎜ + ⎟
⎜ ∂r r ∂θ ⎟ ⎜ rsenθ ∂ϕ ∂r ⎟ ⎜ r ∂θ rsenθ ∂ϕ ⎟ ⎥
⎝ ⎠ ⎜ ⎟ ⎜ ⎟ ⎥
⎝ ⎠ ⎝ ⎠ ⎦⎥
[65]
_________________________________________________________________________________________________________________
E3.6. ECUACIONES DE CONSERVACIÓN PARA FLUIDOS NEWTONIANOS:

ECUACIONES DE NAVIER-STOKES:
∂ρ r
CONTINUIDAD: + ∇ ⋅ ( ρv ) = 0
∂t
r
r µ r 2r ⎛ ∂v r r ⎞
NAVIER-STOKES: ρg − ∇p + ⎣⎡∇ ( ∇ ⋅ v )⎦⎤ + µ∇ v = ρ ⎜ + v ⋅ ∇v ⎟
3 ⎝ ∂t ⎠
⎛ ∂û r ˆ ⎞ r
ENERGÍA: ρ ⎜ + v ⋅ ∇u ⎟ = − p ( ∇ ⋅ v ) + κ∇ 2 T + Φ
⎝ ∂t ⎠
E3.6.1. Ecuaciones de NAVIER-STOKES en coordenadas cartesianas para flujo incompresible:
∂u ∂v ∂w
continuidad: + + =0
∂x ∂y ∂z
∂p ⎛ ∂2u ∂2u ∂2u ⎞ ⎛ ∂u ∂u ∂u ∂u ⎞

ρg x − +µ⎜ 2 + 2 + 2 ⎟ = ρ⎜ + u + v + w ⎟
∂x ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂t ∂x ∂y ∂z ⎠
∂p ⎛ ∂2v ∂2v ∂2v ⎞ ⎛ ∂v ∂v ∂v vu ⎞
movimiento: ρg y − +µ⎜ 2 + 2 + 2 ⎟ = ρ⎜ + u + v + w ⎟
∂y ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂t ∂x ∂y ∂z ⎠
∂p ⎛ ∂2w ∂2w ∂2w ⎞ ⎛ ∂w ∂w ∂w ∂w ⎞
ρg z − +µ⎜ 2 + 2 + 2 ⎟ = ρ⎜ +u +v +w ⎟
∂z ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂t ∂x ∂y ∂z ⎠
⎡ ∂uˆ ⎛ ∂uˆ ∂uˆ ∂uˆ ⎞ ⎤ ⎛ ∂ 2T ∂ 2T ∂ 2T ⎞

ρ ⎢ + ⎜ u + v + w ⎟⎥ = κ ⎜ 2 + 2 + 2 ⎟ +
⎣ ∂t ⎝ ∂x ∂y ∂z ⎠ ⎦ ⎝ ∂x ∂y ∂z ⎠
energía:
⎡ ⎛ ∂u ⎞ 2 ⎛ ∂v ⎞
2
⎛ ∂w ⎞ ⎛ ∂u ∂v ⎞ ⎛ ∂u ∂w ⎞ ⎛ ∂v ∂w ⎞ ⎤
2 2 2 2
+µ ⎢ 2 ⎜ ⎟ + 2 ⎜ ⎟ + 2 ⎜ ⎟ ⎜+ + ⎟ ⎜+ + ⎟ ⎜+ + ⎟ ⎥
⎢⎣ ⎝ ∂x ⎠ ⎝ ∂y ⎠ ⎝ ∂z ⎠ ⎝ ∂y ∂x ⎠ ⎝ ∂z ∂x ⎠ ⎝ ∂z ∂y ⎠ ⎥⎦
Ecuaciones de NAVIER-STOKES en coordenadas cartesianas:
∂ρ ∂ ( ρu ) ∂ (ρv ) ∂ ( ρw )
continuidad: + + + =0
∂t ∂x ∂y ∂z
movimiento:
∂p µ ⎛ ∂ 2 u ∂ 2 v ∂ 2 w ⎞ ⎛ ∂2u ∂2u ∂2u ⎞ ⎛ ∂u ∂u ∂u ∂u ⎞
ρg x − + ⎜ 2+ + ⎟ +µ⎜ 2 + 2 + 2 ⎟ = ρ⎜ + u + v + w ⎟
∂x 3 ⎝ ∂x ∂x∂y ∂x∂z ⎠ ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂t ∂x ∂y ∂z ⎠
∂p µ ⎛ ∂ 2 u ∂ 2 v ∂ 2 w ⎞ ⎛ ∂2v ∂2v ∂2v ⎞ ⎛ ∂v ∂v ∂v vu ⎞
ρg y − + ⎜ + 2+ ⎟ +µ⎜ 2 + 2 + 2 ⎟ = ρ⎜ + u + v + w ⎟
∂y 3 ⎝ ∂y∂x ∂y ∂y∂z ⎠ ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂t ∂x ∂y ∂z ⎠
∂p µ ⎛ ∂ 2 u ∂2v ∂2w ⎞ ⎛ ∂2w ∂2w ∂2w ⎞ ⎛ ∂w ∂w ∂w ∂w ⎞
ρg z − + ⎜ + + 2 ⎟ +µ⎜ 2 + 2 + 2 ⎟ = ρ⎜ +u +v +w ⎟
∂z 3 ⎝ ∂z∂x ∂z∂y ∂z ⎠ ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂t ∂x ∂y ∂z ⎠
⎡ ∂uˆ ⎛ ∂uˆ ∂uˆ ∂uˆ ⎞ ⎤ ⎛ ∂u ∂v ∂w ⎞ ⎛ ∂ 2T ∂ 2T ∂ 2T ⎞

energía: ρ ⎢ + ⎜ u + v + w ⎟⎥ = −p ⎜ + + ⎟+ κ⎜ 2 + 2 + 2 ⎟+Φ
⎣ ∂t ⎝ ∂x ∂y ∂z ⎠ ⎦ ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂x ∂y ∂z ⎠
_________________________________________________________________________________________________________________
E3.6.2. Ecuaciones de NAVIER-STOKES en coordenadas cilíndricas:

⎛ ∂ρ ∂ρ v θ ∂ρ ∂ρ ⎞ ⎡ 1 ∂ ( rv r ) 1 ∂v θ ∂v z ⎤
continuidad: ⎜ + vr + + vz ⎟+⎢ + + ⎥=0
⎝ ∂t ∂r r ∂θ ∂z ⎠ ⎣ r ∂r r ∂θ ∂z ⎦
movimiento:
∂p ⎛ ∂ ⎛ 1 ∂ ( r v r ) ⎞ 1 ∂ 2 v r 2 ∂v θ ∂ 2 v r ⎞ ⎛ ∂v r ∂v r v θ ∂v r v θ2 ∂v ⎞
ρ gr − +µ ⎜ ⎜ ⎟+ 2 − 2 + = ρ + + − + vz r ⎟
⎟⎟ ⎜ v
∂r ⎜ ∂r r ∂r ⎠ r ∂θ r ∂θ ∂z 2 ∂ ∂ ∂θ ∂z ⎠
2 r
⎝ ⎝ ⎠ ⎝ t r r r
1 ∂p ⎛ ∂ ⎛ 1 ∂ ( r v θ ) ⎞ 1 ∂ 2 v θ 2 ∂v r ∂ 2 v θ ⎞ ⎛ ∂v ∂v v ∂v v v ∂v ⎞
ρ gθ − +µ ⎜⎜ ⎜ ⎟+ 2 + 2 + 2 ⎟ = ρ ⎜ θ + vr θ + θ θ + θ r + vz θ ⎟
r ∂θ ⎟
⎝ ∂r ⎝ r ∂r ⎠ r ∂θ r ∂θ ∂z ⎠ ⎝ ∂t ∂r r ∂θ ∂z ⎠
2
r
∂p ⎛ ∂ ⎛ 1 ∂ ( r vz ) ⎞ 1 ∂2vz ∂2vz ⎞ ⎛ ∂v ∂v v ∂v ∂v ⎞
ρ gθ − +µ ⎟+ 2 + 2 = ρ ⎜ z + vr z + θ z + vz z ⎟
∂z ⎜⎜ ∂r ⎜ r ∂r ⎟⎟
⎠ r ∂θ ∂z ⎝ ∂ ∂ ∂θ ∂z ⎠
2
⎝ ⎝ ⎠ t r r
energía:
⎛ ∂uˆ ∂uˆ v θ ∂uˆ ∂uˆ ⎞ ⎛ 1 ∂ ( rv r ) 1 ∂v θ ∂v z ⎞ ⎛ ∂ 2T 1 ∂T 1 ∂ 2T ∂ 2T ⎞
ρ ⎜ + vr + + vz ⎟ = −p ⎜ + + ⎟+ κ⎜ 2 + + + ⎟+Φ
⎝ ∂t ∂r r ∂θ ∂z ⎠ ⎝ r ∂r r ∂θ ∂z ⎠ ⎝ ∂r r ∂r r 2 ∂θ2 ∂z 2 ⎠
E3.6.3. Ecuaciones de NAVIER-STOKES en coordenadas esféricas:

∂ρ ⎡ 1 ∂ (ρr v r ) 1 ∂ ( ρsenθ⋅ v θ ) 1 ∂ ( ρv ϕ ) ⎤
2
continuidad: +⎢ 2 + + ⎥=0
∂t ⎢ r ∂r rsenθ ∂θ rsenθ ∂ϕ ⎥
⎣ ⎦
movimiento:
⎡ ⎛ 2 ∂v r ⎞ ⎛ ⎛ ∂v ⎞⎞ ⎤
⎢ ∂⎜⎜ r ⎟⎟ ⎜ ∂⎜⎜ senθ r ⎟⎟ ⎟ ⎥
∂p 1 ⎝ ∂r ⎠ 1 ⎜ ⎝ ∂θ ⎠⎟ 1 ∂2vr ⎛ ∂v 1 ∂vϕ ⎞⎥
+ µ⎢ 2
2 v
ρg r − + 2 ⎜ ⎟ + senθ − 2 ⎜⎜ v r + θ + θ + ⎟⎟ =
∂r ⎢r ∂r r senθ ⎜ ∂θ ∂ϕ 2 ∂θ tgθ senθ ∂ϕ ⎠⎥
⎢ ⎟ r ⎝ ⎥
⎜ ⎟
⎣⎢ ⎝ ⎠ ⎦⎥
⎛ ∂v ∂v r 1 ∂v r ∂v v θ2 + v ϕ2 ⎞
= ρ⎜ r + v r ⎟
1
+ vθ + vϕ r −
⎝ ⎠
⎡ ⎛ ∂v ⎞ ⎛ ⎛ ∂v ⎞ ⎞ ⎤
⎢ ∂⎜ r 2 θ ⎟ ⎜ ∂⎜ senθ θ ⎟ ⎟ ⎥
1 ∂p 1 ⎝ ∂r ⎠+ 1 ⎜ ⎝ ∂θ ⎠⎟+ 1 ∂ 2
v 2 ∂v 1 ⎛ ∂v ⎞
+ µ⎢ 2 θ θ ϕ ⎥
ρg θ − − − ⎜ v + 2 cos θ ⎟ =
r ∂θ ⎢r ∂r r 2senθ ⎜ ∂θ ⎟ senθ ∂ϕ2 r 2 ∂θ r 2sen 2θ ⎜ θ
∂ϕ ⎟⎠⎥
⎜⎜ ⎟⎟ ⎝
⎢ ⎥
⎣ ⎝ ⎠ ⎦
⎛ ∂v ∂v 1 ∂v 1 ∂v v r vθ − vϕ2 / tgθ ⎞⎟
= ρ⎜ θ + v r θ + vθ θ + vϕ θ −
⎝ ⎠
⎡ ⎛ ∂vϕ ⎞ ⎛ ⎛ ∂v ⎞ ⎞ ⎤
⎢ ∂⎜⎜ r 2 ⎟⎟ ⎜ ∂⎜ senθ ϕ ⎟ ⎟ ⎥
1 ∂p ∂r ⎠ 1 ⎜ ⎝ ⎜ ⎟
∂θ ⎠ ⎟ 1 ∂ vϕ
2
⎛ ∂vϕ ∂v r ⎞⎥
+ µ ⎢⎢ 2 ⎝
1 1
ρg ϕ − + 2 ⎜ ⎟+ − 2 2 ⎜⎜ vϕ − 2 cos θ − 2senθ ⎟ =
rsenθ ∂ϕ r ∂r r senθ ⎜ ∂θ ⎟ senθ ∂ϕ
2
r sen θ ⎝ ∂ϕ ∂ϕ ⎟⎠⎥
⎢ ⎥
⎢⎣ ⎜ ⎟ ⎥⎦
⎝ ⎠
⎛ ∂vϕ ∂vϕ 1 ∂vϕ 1 ∂vϕ v r vϕ + vϕ vϕ / tgθ ⎞⎟
= ρ⎜ + vr + vθ + vϕ +
⎝ ⎠
energía:
⎛ ∂uˆ ∂uˆ v θ ∂uˆ v φ ∂uˆ ⎞
ρ ⎜ + vr + + ⎟=
⎝ ∂t ∂r r ∂θ rsenθ ∂φ ⎠
⎛ 1 ∂ (r2 vr ) 1 ∂ ( v θ senθ ) 1 ∂v φ ⎞ ⎡1 ∂ ∂T ⎞ 1 ∂⎛ ∂T ⎞ ∂2T ⎤
⎟ + κ ⎢ 2 ⎛⎜ r 2
1
= −p ⎜ 2 + + ⎟ + ⎜ θ ⎟ + +Φ
∂θ ⎠ r 2 sen 2 θ ∂ϕ2 ⎥⎦
sen
⎜r ∂r rsenθ ∂θ rsenθ ∂φ ⎟ ⎣ r ∂r ⎝ ∂r ⎠ r 2 senθ ∂r ⎝
⎝ ⎠
_________________________________________________________________________________________________________________
Isaac NEWTON Daniel BERNOULLI Leonhard EULER

(1642 – 1727) (1700 – 1782) (1707 – 1783)
Joseph Louis LAGRANGE Pierre Simon LAPLACE S. POISSON

(1736 – 1813) (1749 – 1827) (1781 – 1840)
Louis Marie Henri NAVIER Augustin Louis de CAUCHY George Gabriel STOKES
(1785 – 1836) (1789 – 1857) (1819 – 1903)
_________________________________________________________________________________________________________________
William FROUDE (1810-1879) Ernst MACH (1838-1919)
Osborne REYNOLDS (1842-1912) Moritz WEBER (1871-1951)
Ludwing PRANDTL (1875-1953) Theodor von KARMAN (1881-1963)
_________________________________________________________________________________________________________________

Análisis Diferencial en Mecánica de Fluidos

Cargado por

Copyright:

Formatos disponibles

Análisis Diferencial en Mecánica de Fluidos

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Análisis Diferencial en Mecánica de Fluidos

Cargado por

Copyright:

Formatos disponibles

2.

Análisis diferencial en Mecánica de Fluidos 1

Escuela Politécnica Superior de Ingeniería de Gijón

Apuntes de Mecánica de Fluidos

ANÁLISIS DIFERENCIAL EN MECÁNICA DE FLUIDOS.

1. Cinemática y dinámica de fluidos.

Julián Martínez de la Calle

ANÁLISIS DIFERENCIAL EN MECÁNICA DE FLUIDOS

1. CINEMÁTICA Y DINÁMICA DE FLUIDOS

1.1. Métodos de análisis en Mecánica de Fluidos: Lagrangiano y Euleriano. 3

2.1. Comportamiento mecánico: tensor de velocidad de deformación. 15

3.1. Ecuación diferencial de conservación de masa: ecuación de continuidad. 21

4.1. Métodos de análisis: Euleriano y Lagrangiano. 38

1. CINEMÁTICA Y DINÁMICA DE FLUIDOS

1.1. MÉTODOS EN MECÁNICA DE FLUIDOS: LAGRANGIANO Y EULERIANO.

El ámbito general de la Mecánica de Fluidos, es la interacción entre fluidos y su entorno. Además el

- La identificación de cada partícula fluida y su seguimiento en el tiempo; es decir, hay que

Al movimiento de un fluido se le denomina flujo, y en su análisis es interesante tener algún tipo de

Fig.1.1.: Trayectoria de una partícula A a lo largo del tiempo: Método Lagrangiano

El vector unitario de la velocidad, es el vector unitario de la dirección tangencial, en determinado en un

1.2. TIPOS DE ANÁLISIS EN MECÁNICA DE FLUIDOS.

1.3. CINEMÁTICA DE FLUIDOS: PROPIEDADES DEL VECTOR VELOCIDAD.

- flujo uniforme, cuando la velocidad no depende de la posición, con lo que en un determinado

1.3.1. DERIVADA SUSTANCIAL O MATERIAL.

La variación temporal del campo de presiones, viene dada por:

dx = u·dt; dy = v·dt; dz = w·dt

en notación compacta, se tiene:

el termino ∇p se denomina gradiente de presión, y su expresión depende del sistema de coordenadas:

El gradiente de velocidad, es un tensor, cuyas componentes cartesianas son:

en donde la aceleración de arrastre, tiene 4 sumandos, respectivamente:

1.3.2. EL OPERADOR NABLA.

El operador NABLA: ∇ , se utiliza para compactar las expresiones cinemáticas. Es un operador

LAPLACIANA: la divergencia de un gradiente, se denomina laplaciana; que cuando se aplica a un escalar da

r ⎛ ∂2u ∂2u ∂2u ⎞ r ⎛ ∂2 v ∂2 v ∂2 v ⎞ r ⎛ ∂2w ∂2w ∂2 w ⎞ r

1.4. DINÁMICA DE FLUIDOS: FUERZAS MACROCÓPICAS.

En el Análisis Diferencial de Fluidos, hemos considerado como volumen de control a la partícula

1.4.1. FUERZAS DE VOLUMEN.

1.4.2. FUERZAS DE SUPERFICIE.

1.4.3. TENSOR DE TENSIONES.

(2) Análogamente sobre las caras y=cte, z=cte.

⎛ ∂τzx ∂τzy ∂σzz ⎞

en donde: σ es el tensor de tensiones; “p” es la presión termodinámica y τ el tensor de tensiones viscosas.

1.5. TIPOS DE FLUJOS

- Flujo estacionario: ( ∂ #/ ∂t = 0 ∧ ∇ # ≠ 0) , en un determinado punto las propiedades del fluido no

2.1. COMPORTAMIENTO MECÁNICO: Tensor de velocidad de deformación.

que además representa la velocidad del aumento (o disminución) unitario de volumen :

Fig.1.8. Dilatación cúbica: debida a la deformación lineal.

Si se tiene un campo de velocidades genérico: u=u(x,y,z), v=v(x,y,z), w=w(x,y,z), se obtienen las

Hasta ahora hemos considerado deformaciones puramente lineales de dilatación o de contracción,

En coordenadas cartesianas, el tensor de velocidad de deformación es:

2.2. FLUIDOS STOKESIANOS: TENSOR DE TENSIONES VISCOSAS.

El tensor de tensiones viscosas, es la diferencia entre el tensor de tensiones y el tensor esférico,

2.3. FLUIDOS NEWTONIANOS.

2.3.1. ECUACIÓN DE NAVIER-POISSON.

2.3.2. TENSOR DE TENSIONES EN UN FLUIDO NEWTONIANO.

Su expresión en coordenadas cartesianas es:

- viscosidad bruta nula: κ=0 ⇒ p = p ; obteniéndose por la condición de STOKES:

2.4. COMPORTAMIENTO TÉRMICO.