Libro de Mecánica de Fluidos

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

Universidad Autónoma de Nuevo León

Facultad de Ingenierı́a Mecánica y Eléctrica

Notas de Clase
Transferencia de Calor

Profesor
Dr. Fausto Alejandro Sánchez Cruz

Nuevo León, México

Abril de 2020
Índice general

1. Introducción 4
1.1. Concepto de fluido . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.2. El fluido como un medio continuo . . . . . . . . . . . . . . . . . . . . . 5
1.3. Propiedades del campo de velocidades . . . . . . . . . . . . . . . . . . . 7
1.3.1. Descripciones Lagrangiana y Euleriana del movimiento . . . . . 7
1.3.2. Campo de Velocidades . . . . . . . . . . . . . . . . . . . . . . . 10
1.4. Propiedades termodinámicas de un fluido . . . . . . . . . . . . . . . . . 11
1.5. Viscosidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
1.5.1. Variación de la viscosidad con la Temperatura . . . . . . . . . . . 14
1.6. Fluidos Newtonianos . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.7. Tensión superficial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
1.8. Patrones de flujo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19

2. Distribución de presión en un fluido 25


2.1. Equlibrio de un elemento de fluido . . . . . . . . . . . . . . . . . . . . . 25
2.2. Presión Absoluta y Relativa . . . . . . . . . . . . . . . . . . . . . . . . . 27
2.3. Distribución de la presión hidrostática . . . . . . . . . . . . . . . . . . . 28
2.3.1. Efecto de la variación de la gravedad . . . . . . . . . . . . . . . 29
2.3.2. Presión hidrostática en lı́quidos . . . . . . . . . . . . . . . . . . 30
2.3.3. Barómetro de Mercurio . . . . . . . . . . . . . . . . . . . . . . . 30
2.3.4. Presión hidrostática en gases . . . . . . . . . . . . . . . . . . . . 31
2.4. Manometrı́a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
2.5. Fuerzas hidrostáticas . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
2.5.1. Fuerzas hidrostáticas sobre superficies curvadas . . . . . . . . . . 36
2.6. Flotación y estabilidad . . . . . . . . . . . . . . . . . . . . . . . . . . . 42

1
2

2.6.1. Estabilidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
2.7. Distribución de la presión en un fluido en movimiento como cuerpo rı́gido 43
2.7.1. Aceleración uniforme . . . . . . . . . . . . . . . . . . . . . . . . 44

3. Relaciones integrales para un volumen de control 46


3.1. Layes fı́sicas fundamentales de la mecánica de fluidos . . . . . . . . . . . 46
3.1.1. Flujo volumétrico y flujo másico . . . . . . . . . . . . . . . . . . 47
3.2. Teorema de transporte de Reynolds . . . . . . . . . . . . . . . . . . . . . 48
3.2.1. Aproximaciones unidimensionales al término de flujo . . . . . . . 49
3.3. Conservación de masa. Ecuación de Continuidad . . . . . . . . . . . . . 50
3.4. Conservación de la cantidad de movimiento . . . . . . . . . . . . . . . . 52
3.4.1. Aproximaciones unidimensionales al término de flujo de cantidad
de movimiento . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
3.5. Ecuación de la Energı́a . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
3.5.1. Aproximación unidimensional del término de flujo de cantidad de
energı́a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
3.5.2. Ecuación de Bernoulli en estado permanente . . . . . . . . . . . 58

4. Flujo en Conductos 61

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


3

BIBLIOGRAFIA
White, F. M. Fluid Mechanics. Mc. Graw-Hill
Streeter, V. L. Mecánica de Fluidos. Mc. Graw-Hill.
Shames I. H. Mecánica de Fluidos. Mc. Graw-Hill.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


Capı́tulo 1

Introducción

1.1. Concepto de fluido


¿Qué es un fluido?

Un fluido es un material que no puede soportar esfuerzos cortantes aplicados exter-


namente, y como respuesta a ellos, se deforma de manera continua.

El flujo del fluido es el efecto de la deformación continua del fluido que es sometido
a la acción de los esfuerzos.

Los fluidos no tienen forma definida.

Debido a su capacidad para deformarse con el más mı́nimo esfuerzo, los fluidos
muestran un comportamiento complejo cuando fluyen.

Podemos observar flujo laminar o turbulento

¿Cómo estudiamos la mecánica de fluidos?

Métodos analı́ticos

Métodos experimentales

Métodos computacionales

Métodos combinados

4
5

Figura 1.1: Respuesta de sólidos y fluidos a los esfuerzos. Tomada de White, F. M. Fluid Mechanics. Mc.
Graw-Hill.

1.2. El fluido como un medio continuo


Métodos para derivar las ecuaciones:

Estadı́stico (Leyes de la mecánica; Teorı́a de la probabilidad) ESTO NO

Suposición de medio continuo ESTO SÍ

Hipótesis de medio continuo: Espacio completamente lleno (no se consideran molécu-


las) ρ , ν , T , etc. varı́an como funciones continuas
¿Dificultades? Cálculo de propiedades del fluido
Ejemplo: La densidad

m
ρ= (1.1)
V
¿De qué tamaño es conveniente tomar V?
Muy Pequeño

m (t)
ρ (t) = ; No es constante por el
V
movimiento molecular

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


6

Figura 1.2: Volumen de control muy pequeño.

Muy Grande

Figura 1.3: Volumen de control muy grande

ρ = ρ promedio
ρ = ρ (V)

Entonces, para que no influya el movimiento molecular, δV ∗ =?


Experimentalmente se encuentra que δV ≈ 10−9 mm3 .
Para que la suposición de medio continuo se cumpla debe de ocurrir que:

λ�a�D (1.2)

donde

a : Longitud del volumen


D : Escala macroscópica
λ : Trayectoria libre media

En una medición de densidad del aire a condiciones normales:


Para el caso del aire en condiciones normales 10−9 mm3 =⇒ 3 × 107 moléculas.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


7

Figura 1.4: El tamaño del volumen de control y la densidad. Tomada de White, F. M. Fluid Mechanics. Mc.
Graw-Hill.

¿Cuándo no podemos suponer un medio continuo?



Ondas de choque, D ∼ λ
No se cumple en
Gases enrarecidos, λ ∼ D

1.3. Propiedades del campo de velocidades


1.3.1. Descripciones Lagrangiana y Euleriana del movimiento
Descripción Lagrangiana




 Se sigue una porción de fluido



 y se observan sus cambios
Lagrangiano 



 Las variables independiantes son:

 x0 , y0 , z0 , t
Para este caso caso:

r = r(t)
v = v(t)
ρ = ρ(t)
T = T (t)

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


8

Figura 1.5: Descripción de Lagrange

donde,

dr dv
v= ; a = ; etc. (1.3)
dt dt
¡¡La descripción Lagrangiana casi no se usa en el estudio de la mecánica de fluidos!!

Descripción Euleriana

Figura 1.6: Descripción de Euler


ρ = ρ (x, y, z, t)
Campos escalares (1.4)
T = T (x, y, z, t)

v = v (x, y, z, t)
Campos Vectoriales (1.5)
a = a (x, y, z, t)

¡¡Sı́ se usa en el estudio de la mecánica de fluidos!!

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


9

La importancia en la relación entre ambos enfoques se deba a que:

Leyes básicas ⇒ Sistema Lagrangiano


Mecánica de fluidos ⇒ Sistema Euleriano
¿Cómo ligar los dos enfoques?

Sea α una variable arbitraria del campo (ρ, T, u, etc.), donde

α = α (x, y, z, t) (1.6)

Entonces,

∂α ∂α ∂α ∂α
δα = δt + δx + δy + δz
∂t ∂x ∂y ∂z
(δα lagrangiana)

o bien,

δα ∂α ∂α δx ∂α δy ∂α δz
= + + +
δt ∂t ∂x δt ∂y δt ∂z δt
en el lı́mite ∂t → 0
δx δy δz
=u =v =w
δt δt δt

Dα ∂α ∂α ∂α ∂α
= +u +v +w (1.7)
Dt ∂t
���� ∂x ∂y ∂z
����������������������������������������
cambio local cambio convectivo


: Derivada material, total o sustancial de α
Dt
Entonces,

Izquierda (lagrangiano) = Derecha (euleriano)

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


10

1.3.2. Campo de Velocidades


El campo de velocidades es la propiedad del flujo más importante en la mecánica de
fluidos

v = v (x, y, z, t)

Determinar el campo de velocidades es prácticamente resolver el problema de estu-


dio.

Otras propiedades del flujo se obtienen directamente del campo de velocidades, e.g.
el campo de presiones.

En estudio de la transferencia de calor se determina la temperatura y el flujo de calor


una vez que se conoce el campo de velocidades

El campo de velocidades es una función vectorial de variable vectorial:

v (x, y, z, t) = u (x, y, z, t) i + v (x, y, z, t) j + w (x, y, z, t) k (1.8)

donde

u : componente de la velocidad en x
v : componente de la velocidad en y
w : componente de la velocidad en z

Otras propiedades cinemáticas del flujo están relacionadas con el campo de velocida-
des:
1. Vector desplazamiento �
r= vdt (1.9)

2. Aceleración
dv
a= (1.10)
dt
3. Flujo volumétrico �
Q= (v · n) dA (1.11)

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


11

4. Razón de expansión volumétrica


1 dV
=∇·v (1.12)
V dt
5. Velocidad angular local
1
ω= ∇×v (1.13)
2

1.4. Propiedades termodinámicas de un fluido


El flujo de fluidos está ligado a las propiedades de éstos. Las propiedades termodinámi-
cas de un fluido son:

Presión, P (Pa, N/m2 )

Densidad, ρ (kg/m3 )

Temperatura, T (K)

Cuando se realiza el análisis con intercambio de trabajo y flujo de calor es necesario


considerar:

Energı́a interna, e (J)

Entalpı́a, h = �
u + P/ρ (J/kg K)

Entropı́a, s (J/kg K)

Calor especı́fico, c p ó cv (J/kg K)

Además, los efectos de fricción y la transferencia de calor se modelan como fenómenos


de transporte y para ello es necesario definir:

Viscosidad dinámica, µ (kg/m s)

Conductividad térmica, k (W/m K)

Todas son propiedades termodinámicas y en general pueden variar en un sistema. Va-


mos a analizar algunas de ellas con un poco más de detalle.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


12

Presión (P)

Fuerza normal local por unidad de área. De manera formal:


dF
. P= (1.14)
dA
Generalmente necesitamos conocer la fuerza que se ejerce en una área de interés, para
ello:

dF = �
PdA,
∴ F = PdA. (1.15)
A

Densidad (ρ)

Es la masa por unidad de volumen. Es altamente variable en gases y casi constante en


lı́quidos. De manera formal:
dm
ρ = , (1.16)
dV
m : masa
V : volumen.

Temperatura (T)

Es una medida de la cantidad de energı́a interna que tiene un cuerpo

La temperatura es proporcional a la energı́a interna molecular

Es común el uso de la escala Celsius (◦ C) y la escala Fahrenheit (◦ F), cuya relación


está dada por:
◦ ◦
F = C · (9/5) + 32,

C = (◦ F − 32) · 5/9.

En el análisis de procesos termodinámicos es necesario utilizar escalas absolutas de


temperatura, Kelvin (K) en SI y Rankine (R) en el sistema británico:

K = C + 273,15,
K = (◦ F + 459,67) · 5/9,

R = F + 459,67.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


13

Peso especı́fico (γ)

Es el peso por unidad de volumen. Ya que el peso es W = mg, entonces


� �
γ = ρg N/m3 . (1.17)

Energı́a Interna (e)

Es la energı́a total acumulada en un cuerpo por unidad de masa:


1
u + V 2 + gz,
e = � (1.18)
2
u : energı́a interna molecular (función de P y T en una sustancia pura),

V : módulo de la velocidad,
z : coordenada vertical.

1.5. Viscosidad
Viscosidad (µ)

Es una propiedad que caracteriza el comportamiento mecánico de un fluido en par-


ticular.

“Cuando se ejerce un esfuerzo cortante en el fluido, éste se deforma con una rapidez
proporcional al esfuerzo e inversamente proporcional al coeficiente de viscosidad
µ”, es decir que:

δθ
τ∝ . (1.19)
δt
De nuestra figura anterior:
δuδt
tan δθ = ,
δy
en el caso de que se trate de un cambio infinitesimal: tan δθ = dθ, por lo que:
du dt
dθ ≈ ,
dy
dθ du
≈ .
dt dy

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


14

Figura 1.7: Rapidez de deformación. Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

Es decir que
du
τ∝ . (1.20)
dy
Los fluidos que llamamos Newtonianos cumplen que la constante de proporcionalidad
es µ, es decir

du
τ = µ , (1.21)
dy
Ley de Newton (1687).

Figura 1.8: Propiedades de fluidos comunes. Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

1.5.1. Variación de la viscosidad con la Temperatura


La temperatura tiene un efecto importante sobre la viscosidad
En los gases la viscosidad aumenta con la temperatura

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


15

Dos aproximaciones comunes para el comportamiento de la viscosidad en los gases


son:
� �n
µ T
= ; Ley de potencias (1.22)
µ0 T0
µ (T/T 0 )3/2 (T 0 + S )
= ; Ley de Sutherland (1.23)
µ0 T0 + S
donde µ0 es la viscosidad a una temperatura de referencia T 0 , usualmente T 0 = 273 K.
Para aire:

n ≈ 0,7 (1.24)
S = 110 K

En los lı́quidos la viscosidad decrece con la temperatura


En modelo simple es
µ �T � � T �2
0 0
ln ≈a+b +c (1.25)
µ0 T T
donde para el agua con T 0 = 273,16 K,

µ0 = 0,001792 kg/m · s (1.26)


a = −1,94
b = −4,8
c = 6,74

1.6. Fluidos Newtonianos


Los fluidos que no siguen la relación lineal
du
τ=µ , (1.27)
dy
se llaman fluidos no newtonianos y son estudiados con detalle por la Reologı́a
Dependiendo de su comportamiento, existen algunas clasificaciones útiles:
Fluido dilatante: Aumenta su resistencia al flujo con el aumento del esfuerzo cortante
aplicado
Fluido Pseudoplástico: Disminuye su resistencia al flujo con el aumento del esfuerzo
cortante aplicado

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


16

Fluido Plástico: se comporta como el fluido Pseudoplástico, pero mucho más notable.
Plástco Ideal de Bringham: Resiste el esfuerzo y empieza a fluir una vez que se
supera un valor crı́tico, después su comportamiento es lineal

Figura 1.9: Rapidez de deformación y esfuerzo cortante. Tomada de White, F. M. Fluid Mechanics. Mc.
Graw-Hill.

Algunos otros fluidos tienen comportamientos más complicados, e.g. pueden mostrar
dependencia con el estado transitorio:
Fluido Reopéctico: Requiere que se aumente paulatinamente el esfuerzo cortante para
mantener constante la rapidez de deformación
Fluido Tixitrópico: Requiere que se disminuya paulatinamente el esfuerzo cortante
para mantener constante la rapidez de deformación

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


17

Figura 1.10: Realción entre el esfuerzo cortante y el tiempo bajo rapidez de deformación constante

1.7. Tensión superficial


La tensión superficial se manifiesta en la interface de un lı́quido con un segundo lı́quido
o un gas.
Si se hiciera una cortada de tamaño dL en la superficie interfacial, se experimentarı́an
fuerzas iguales y opuestas de magnitud σdL normales a la dirección del corte y paralelas
a la superficie interfacial.
σ es el coeficiente de tensión superficial
Los coeficientes de dos interfaces comunes son:

σ = 0,073 N m para aire-agua


σ = 0,48 N m para aire-mercurio

Si la interface está curvada, entonces se presenta una diferencia de presión a través de


la interface, donde la presión es mayor del lado cóncavo.
La diferencia de presiones a través de la interface es igual a las fuerzas de tensión
superficial.
en el caso de la figura a, el balance es entonces:

2RL ΔP = 2σL

y entonces
σ
ΔP = (1.28)
R
Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL
18

Figura 1.11: Diferencia de presión en la interface de un lı́quido. Tomada de White, F. M. Fluid Mechanics.
Mc. Graw-Hill.

Para la figura b, el incremento de la presión en el interior de la gota está balanceado


con la fuerza de tensión superficial en el anillo, es decir

πR2 ΔP = 2πRσ

y entonces

ΔPgota =
R
y como una burbuja tiene dos interfaces el aumento de presión en ellas es ΔPburbu ja =
2ΔPgota = 4σ/R.
Para el caso de la figura c, una superficie cuyos radios de curvatura principales son R1
y R2 , entonces:
� �
ΔP = σ R−1 −1
1 + R2 (1.29)

Otro parámetro importante es el águlo de contacto θ que aparece en las interfaces que
está en contacto con superficies sólidas como en la figura siguiente.

Figura 1.12: Ángulo de mojado entre un lı́quido y un sólido. Tomada de White, F. M. Fluid Mechanics. Mc.
Graw-Hill.

Si el ángulo es menor a 90◦ se dice que el fluido moja al sólido

Si el ángulo es menor a 90◦ se dice que el fluido no moja.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


19

El ángulo de mojado es altamente dependiente de las condiciones fı́sico-quı́micas


de la interface sólido-lı́quido.

Ejemplo Derive una expresión para determinar el cambio de altura h en un tubo circular
que contiene un fluido con tensión superficial σ y águlo de contacto θ. Determine h si se
trata de agua y θ = 0◦ ; compara con el caso de mercurio y θ =130◦ .

Figura 1.13: Altura de equilibrio. Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

1.8. Patrones de flujo


Los patrones de flujo se pueden visualizar por medio de varias técnicas.
Cuatro tipos básicos de lı́nea pueden ser empleadas para visualizar flujos:
1. Lı́nea de corriente (streamline). Es una lı́nea que es tangente al vector de velocidad
en un cierto instante
2. Senda (pathline). Es la trayectoria que verdaderamente sigue una partı́cula dada de
fluido
3. Lı́nea de emisión (streakline). Es un lugar geométrico que está compuesto por
partı́culas que con anterioridad han pasado por un punto en particular.
4. Lı́nea fluida (timeline). Es un grupo de partı́culas de fluido que forman una lı́nea
en un cierto instante

La lı́nea de corriente se calcula convenientemente de forma matemática, las demás


se pueden generar experimentalmente

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


20

Una lı́nea fluida se genera por ejemplo por medio de emisión de burbujas desde un
alambre

La traza se puede obtener fotografiando partı́culas suspendidas que son arrastradas


por el flujo

En estado permanente la lı́nea de corriente, la traza, lı́nea de emisión y la lı́nea de


tiempo son idénticas entre sı́.

Figura 1.14: Lı́neas de corriente. Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

Considerando un vector arbitrario como el de la siguiente figura, un arco de longitud dr


que se pueda considerar lı́nea de corriente debe ser paralelo a V, entonces sus respectivas
componentes (arco-vector) deberán estár en la misma proporción:
dx dy dz dr
= = = (1.30)
u v w V

Figura 1.15: Vector de velocidad y componentes de un fragmento de la longitud de arco de la lı́nea de


corriente. Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


21

Si las componentes de la velocidad (u, v, w) se conocen como funciones de la posición


y el tiempo, esas ecuaciones se pueden integrar para encontar las lı́neas de corriente que
pasan inicialmente por un punto (x0 , y0 , z0 , t0 ).
La traza es el despalzamiento de una partı́cula y se determina de la integración de las
componentes de la velocidad:
� � �
x= u dt; y = v dt; z = w dt (1.31)

Si las componentes de la velocidad (u, v, w) se conocen como funciones de la posición


y el tiempo, la integración comienza en un punto inicial (x0 , y0 , z0 , t0 ).
Ejemplo Se tiene una distribución de velocidad bidimensional:

u = K x, v = −Ky, w = 0

donde la K es una constante positiva. Determine y grafique las lı́neas de corriente inclu-
yendo las direcciones.

Figura 1.16: Lı́neas de corriente. Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

Las lı́neas de emisión se obtienen experimentalmente con humo, tinta o burbujas.


Métodos para la visualización de flujos incluyen los siguientes:
1. Inyección de tinta, humo o burbujas
2. Viruta o polvo arrastrados en la superficie libre de lı́quidos
3. Partı́culas trazadoras de densidad neutral
4. Técnicas ópticas que detectan cambios de densidad: shadowgraph, Schlieren e in-
terferometrı́a
5. Hilos o lanas sujetos a superficies que limitan el flujo

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


22

Figura 1.17: Descripción del flujo. Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

6. Sustancias que se evaporan desde superficies sólidas


7. Fluidos liminicentes o aditivos, etc.
Problema Un bloque con peso W desliza hacia abajo en un plano inclinado lubricado
con una pelı́cula delgada de aceite (ver figura). El espesor de la pelı́cula es h y área de
contacto con es A. Suponiendo una distribución de velocidad lineal en el interior de la
pelı́cula de fluido, determine una expresión para determinar la velocidad terminal (acele-
ración cero) V del bloque.

Figura 1.18: Bloque que desliza en un plano inclinado. Tomada de White, F. M. Fluid Mechanics. Mc.
Graw-Hill.

Problema Una placa delgada está separada de dos placas planas fijas por lı́quidos muy
viscosos µ1 y µ2 , respectivamente. La separaciones h1 y h2 no son iguales. El área de
contacto entre la placa central y cada uno de los fluidos es A. Suponiendo un perfil de
velocidad lineal en cada fluido, determinar la fuerza F necesaria para jalar la placa a una
velocidad V. ¿Es necesario que exista alguna relación entre las dos viscosidades, µ1 y µ2 ?
Problema Un disco de radio R rota a una velocidad angular Ω dentro de un contenedor
con forma de disco que está lleno de un aceite de viscosidad µ. Suponiendo un perfil
de velocidad lineal y despreciando el esfuerzo cortante sobre los bordes, determine una

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


23

Figura 1.19: Placa horizontal que separa dos fluidos. Tomada de White, F. M. Fluid Mechanics. Mc. Graw-
Hill.

ecuación para la potencia que requiere el disco para girar con velocidad constante.

Figura 1.20: Disco giratorio inmerso en un aceite. Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

Problema Determine una expresión para la altura capilar h para un fluido con tensión
superficial σ y ángulo de contacto θ entre dos placas paralelas separadas una distancia W.
¿Cuánto vale h para el caso de agua a 20 ◦ C con W = 0.5 mm

Figura 1.21: Altura de equilibrio entre placas paralelas. Tomada de White, F. M. Fluid Mechanics. Mc.
Graw-Hill.

Problema Obtener una expresión que permita determinar el ascenso capilar h(x) que
habrá entre dos placas que forman un ángulo 2θ, la tensión superficial es σ y la altura de
las placas sobre el nivel del depósito es h0 .

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


24

Figura 1.22: Alturas de equilibrio entre placas oblicuas. Tomada de White, F. M. Fluid Mechanics. Mc.
Graw-Hill.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


Capı́tulo 2

Distribución de presión en un fluido

2.1. Equlibrio de un elemento de fluido


Cuando el fluido está en reposo, llamada condición hidrostática, la variación de la
presión en el fluido se debe únicamente al efecto del peso

Cuando el fluido tiene un movimiento de cuerpo rı́gido, la presión se puede calcular


fácilmente debido a que el fluido se encuentra libre de esfuerzos cortantes.

En un fluido en reposo el esfuerzo normal sobre cualquier plano tiene un valor único
llamado presión del fluido

La presión o cualquier otro esfuerzo no causan efecto en el fluido a menos que varı́e
espacialmente.

Digamos que tenemos la presión como

p = p (x, y, z, t) (2.1)

De acuerdo a la figura, la fuerza neta en la dirección horizontal es


� �
∂p ∂p
dF x = p dy dz − p + dx dy dz = − dx dy dz
∂x ∂x
Aplicando esto a las direcciones y y z:
∂p
dFy = − dx dy dz
∂y
∂p
dFz = − dx dy dz
∂z

25
26

Figura 2.1: Fuerzas debidas a la presión que se ejercen sobre un elemento de fluido. Tomada de White, F.
M. Fluid Mechanics. Mc. Graw-Hill.

es decir que vectorialmente tenemos:


� �

− ∂p→− ∂p→ − ∂p→ −
d F presión = − i + j + k dx dy dz
∂x ∂y ∂z


o bien, si f es la fuerza por unidad de volumen:


f presión = −∇p (2.2)

Si consideramos el efecto de la gravedad, →


−g :



d F grav = −ρ→
−g dx dy dz

o bien


f grav = ρ→
−g (2.3)
Las fuerzas por efectos viscosos en el fluido son:
� 2 �

− ∂ → − ∂2 →
− ∂2 →
− →

f visc = µ 2
V + 2 V + 2 V = µ∇2 V (2.4)
∂x ∂y ∂y
Aplicando la segunda ley de Newton:
�→ − → − →
− →

ρ→
−a = f = f presión + f grav + f visc (2.5)

o bien:
ρ→ −g + µ∇2→
−a = −∇p + ρ→ −
V

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


27

y entonces
�− → � →

∇p = ρ →
g − −a + µ∇2 V (2.6)

Examinando nuestra última ecuación tenemos al menos cuatro casos especiales:


1. Fluido en reposo o con velocidad constante. La aceleración y los efectos visco-
sos desaparecen. La p sólo depende de la densidad y la gravedad. Esta se conoce como
condición hidrostática.
2. Traslación y rotación como cuerpo rı́gido. El efecto viscoso desaparece y la pre-
�− → �
sión depende del término ρ → g − −a .
� �


3. Movimiento irrotacional ∇ × V = 0 . El efecto viscoso desaparece y en este caso
existe una integral exacta llamada ecuación de Bernoulli que permite determinar la distri-
bución de presión.
4. Movimiento viscoso arbitrario. No existe simplificación al problema y no existen
reglas generales.

2.2. Presión Absoluta y Relativa


Es común medir la presión de dos formas.
a. La presión absoluta es la que se mide con respecto a la presión nula
b. La presión relativa está referida a la presión atmosférica local

Si la presión es superior a la presión atmosférica se denomina presión manométrica.

Si la presión es inferior a la presión atmosférica se denomina presión de vacı́o.

La figura ilustra la relación entre las presiones absoluta y relativa.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


28

Figura 2.2: Presión absoluta y relativa. Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

2.3. Distribución de la presión hidrostática


Cuando el fluido está en reposo o se mueve a velocidad constante, entonces →
−a = 0 y


∇2 V = 0, por lo tanto,
∇p = ρ→
−g (2.7)

Esta corresponde a la distribución de la presión hidrostática y es adecuada para


fluidos en reposo independientemente de su viscosidad.

∇p representa la magnitud y dirección de la máxima razón de incremento espacial


de la propiedad p.

∇p es perpendicular a cualquier superficie con p constante.

En la condición hidrostática, el máximo aumento de presión ocurre en la dirección


de la gravedad.

Si definimos un sistema cartesiano con el eje z apuntando verticalemente hacia arriba,


entonces
−g = −g→
→ −
k

y en ese caso tenemos que


∂p→− ∂p→ − ∂p→ − →

∇p = i + j + k = −ρg k (2.8)
∂x ∂y ∂z

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


29

es decir que
∂p ∂p ∂p
= 0; = 0; = −ρg
∂x ∂y ∂z
Entonces es posible usar derivada total:
dp
= −ρg (2.9)
dz
e integrando: � 2
p2 − p 1 = − ρg dz (2.10)
1
Observaciones:

En un fluido estático, la distribución de presión varı́a solamente con la distancia


vertical y es independiente de la forma del contenedor

Entonces, la presión es la misma en todos los puntos alojados en un mismo plano


horizontal

La presión aumenta con la profundidad en el fluido.

Figura 2.3: Presión hidrostática en vasos comunicados. Tomada de White, F. M. Fluid Mechanics. Mc.
Graw-Hill.

2.3.1. Efecto de la variación de la gravedad


Para un planeta esférico con densidad uniforme, la variación de la gravedad es:
� r �2
0
g = g0 (2.11)
r

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


30

donde r0 es el radio del planeta y g0 es la gravedad en la superficie. Para la tierra, , y para


casos tı́picos de ingenierı́a, la máxima diferencia existirı́a entre lo más profundidad del
mar, cerca de 11 km y la altura de transportes supersónicos, cerca de 20 km. Esto da una
cambio en la gravedad de alrededor de 0.6 %.
Es entonces común despreciar las variaciones de la gravedad.

2.3.2. Presión hidrostática en lı́quidos


Los lı́quidos se comportan prácticamente como fluidos incompresibles
Generalmente se desprecia la variación de la densidad
Entonces, al integrar nuestra ecuación para la presión hidrostática considerando gra-
vedad y densidad constantes:

P2 − P1 = −ρg (z2 − z1 ) = −γ (z2 − z1 ) (2.12)

donde γ = ρg es el peso especı́fico del fluido.

2.3.3. Barómetro de Mercurio


Una aplicación práctica de la ecuación de la hidrostática es el berómetro que se emplea
en la medición de la presión atmosférica.
Se llena un tubo con mercurio y se invierte en un depósito, como en la figura

Figura 2.4: Barómetro de Mercurio. Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

Al ser invertido el tubo se produce un vacı́o en la parte superior

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


31

La presión atmosférica fuerza al mercurio a ascender por la columna una distancia h,


La altura h se puede determinar mediante una balance de fuerzas en la dirección verti-
cal:
Pa − 0 = −γHg (0 − h)

o bien
Pa
h= (2.13)
γHg
A nivel del mar Pa = 101, 350 Pa, γHg = 133100 N/m3 , entonces h = 0,761 m ó 761
mm.

2.3.4. Presión hidrostática en gases


Los gases son compresibles con la densidad cercanamente proporcional a la presión
La densidad se debe considerar variable si la integración de la ecuación de la hidrostáti-
ca se hace bajo cambios importantes en la presión.
La ecuación de gas ideal establece que:
P
Pv = RT =
ρ
P
ρ = (2.14)
RT
donde R es la constante del gas y T es la temperatura absoluta.
La ecuación de la hidrostáica para este caso es:

dP P
= −ρg = − g
dz RT
o bien � �
g (z2 − z1 )
P2 = P1 exp − (2.15)
RT 0
Nuestra última ecuación es una buena aproximación, sin embargo, en la atmósfera la
temperatura varı́a linealmente con la altura:

T ≈ T 0 − Bz (2.16)

donde T 0 = 288,16 K, B = 0,0065 K/m.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


32

Entonces para determinar la variación de la presión en la atmósfera tenemos que:


dP P P
= −ρg = − g = − g
dz RT R (T 0 − Bz)
o bien � �
dP g dz
=−
P R (T 0 − Bz)
y entonces:
� �g/(RB)
Bz
P = Pa 1 − (2.17)
T0
donde para el aire g/ (RB) = 5,26, R = 287 m2 /s2 K.

2.4. Manometrı́a
El cambio en elevación, z2 − z1 , en un lı́quido es equivalente al cambio en la presión
(P2 − P1 ) / (ρg)
Entonces una o más columnas estáticas de fluido se pueden utilizar para medir las
diferencias de presión entre dos puntos
Cuando se tienen varios fluidos, se debe utilizar la densidad de cada uno de ellos con
la longitud que corresponde al fluido en particular
Consideremos la diferencia entre la presión de los puntos 1 y 5 en el dibujo:

Figura 2.5: Estratificación de fluidos debido a la densidad. Tomada de White, F. M. Fluid Mechanics. Mc.
Graw-Hill.

� �
P5 − P1 = g −ρ0 (z2 − z1 ) − ρw (z3 − z2 ) − ρG (z4 − z3 ) − ρ M (z5 − z4 )

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


33

O escrito de otra manera:

Paba jo = Parriba + ρg |Δz| (2.18)

Ejemplo La aplicación clásica de un manómetro ocurre cuando dos tubos en Ü”de


igual longitud se emplean para medir la diferencia en presión entre dos puntos ubicados
ene le mismo plano horizontal. Generalmente se mide la diferencia de presión del fluido
antes y después de un dispositivo. Determine la diferencia entre las presiones en los puntos
a y b empleado la manometrı́a.

Figura 2.6: Caı́da de presión en un dispositivo. Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

Ejemplo Se emplea un manómetro en B para determinar la presión en A en un flujo de


agua. Si la presión en B es de 87 kPa, estime la presión en A, suponga que ambos fluidos
están a 20◦C.

Figura 2.7: Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

Ejemplo En la figura, en el manómetro se leen 1.5 kPa. Los fluidos están a 20◦ C.
Determine las elevaciones en z de los niveles de los lı́quidos en los tubos B y C.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


34

Figura 2.8: Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

2.5. Fuerzas hidrostáticas


Para determinar las fuerzas debidas a la presión hidrostática que actúan en una super-
ficie sumergida considere la figura.

Figura 2.9: Fuerzas hidrostáticas en una superficie plana sumergida. Tomada de White, F. M. Fluid Mecha-
nics. Mc. Graw-Hill.

Se trata de una placa plana inclinada un ángulo θ.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


35

La presión hidrostática en el interior de fluido es:

P = Pa + ρgh (2.19)

Entonces la fuerza hidrostática en una de las placas es:


� � �
F= PdA = (Pa + ρgh) dA = Pa A + ρg hdA

En este caso, de la figura:


h = ξ sin θ

Ahora la definición del centroide es



1
ξCG = ξdA
A
Para la placa plana θ es constante, entonces:
� �
F = Pa A + ρg hdA = Pa A + ρg sin θ ξdA
F = Pa A + ρg sin θξCG A

pero como sin θξCG = hCG , entonces

F = Pa A + ρghCG A = A (Pa + ρghCG )


F = PCG A (2.20)

La fuerza sobre la cara de cualquier superficie sumergida en un fluido es igual a


la presión que hay en el centro de gravedad de dicha cara multiplicada por su área,
independientemente de la forma de la placa o de su ángulo de inclinación θ.
La fuerza sobre la cara actúa en un punto llamado centro de presiones CP cuyas coor-
denadas son (xCP , yCP ).
Para una superficie sumergida:

F = ρghCG A (2.21)
I xy sin θ
xCP = − (2.22)
hCG A
I xx sin θ
yCP = − (2.23)
hCG A

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


36

donde I xx es el momento de inercia del área de la placa respecto a su eje central x calculado
en el plano de la placa, y I xy es el momento de inercia de la placa calculado en el plano de
la placa respecto a ejes que pasan por el centro de gravedad.
La figura muestra los momentos de inercia respecto al centroide para varias formas
planas.

Figura 2.10: Momentos de inercia. Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

Ejemplo Una compuerta de 5 ft de ancho está articulada en la lı́nea L y descansa sobre


una pared lisa en la lı́nea A. Calcular (a) la fuerza sobre la compuerta debida a la presión
del agua; (b) la fuerza horizontal P que se ejerce sobre la pared en A; (c) las reacciones en
B.

2.5.1. Fuerzas hidrostáticas sobre superficies curvadas


Las fuerzas sobre superficies curvadas se calculan separando las componentes hori-
zontal y vertical.

La componente horizontal de la fuerza en una superficie curvada es igual a la fuerza


sobre el área plana formada por la proyección en la dirección vertical de la superficie

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


37

Figura 2.11: Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

Figura 2.12: Balance de fuerzas. Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

curvada.

La componente vertical de la presión sobre una superficie curvada es igual al peso


de la columna de fluido, aire y lı́quido sobre la superficie curvada.

Ejemplo Una presa tiene forma parabólica z/z0 = (x/x0 )2 como se muestra en la figura,
con x0 = 10 ft y z0 = 24 ft. El fluido es agua ρg = 62,4 lb/ft3 y la presión atmosférica se
puede despreciar. Calcula las componentes horizontal y vertical sobre la presa y la posición
del centro de presiones. La anchura de la presa es de 50 ft.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


38

Figura 2.13: Componentes de la fuerza hidrostática en una superficie curvada. Tomada de White, F. M.
Fluid Mechanics. Mc. Graw-Hill.

Para calcular la componente horizontal determinamos la fuerza que se ejerce sobre la


proyección vertical de la superficie curva. Esta es igual a la fuerza que se ejerce sobre el
punto ubicado en el centroide del rectángulo que resulta de proyectar la superficie curva
en la dirección vertical:
z0
hCG = = 12ft
2 � �
FH = gρhCG A proy = 62,4lbf/ft3 (12ft) (24ft) (50ft)
F H = 899 × 103 lbf

La lı́nea de acción de F H está ubicada por debajo del CG del rectángulo proyectado por
la superficie curva en la dirección vertical. Para un rectángulo el CP está ubicado según
Iyy sin θ bL3
zCP = − ; Iyy =
hCG A 12
y para este caso:

wz30 50 (24)3
Iyy = = = 57600
12 12
57600 sin (π/2)
zCP = − = −4
12 (24) (50)
Como el centro de masa está ubicado en z = 12, el centro de masa está en z = 12−4 = 8
respecto al origen ubicado en el fondo de la presa.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


39

La componente vertical la calculamos como el peso del agua:

FV = gρV = gρwAba jo curva


� z0
Aba jo curva = x (z) dz
0

donde x (z) se lee ”x como función de z”. En este caso


� �2 � �1/2
z x z
= ⇒ x (z) = x0
z0 x0 z0
por lo tanto
� z0 � �1/2 � �z
z x0 2 3/2 0
Aba jo curva = x0 dz = 1/2 z
0 z0 z0 3 0
2
Aba jo curva = x0 z0
3
y entonces
� � � �
2 2
FV = gρw x0 z0 = 62,4 (50) (10) (24)
3 3
FV = 4. 992 × 105 lbf

La magnitud de la fuerza que ejerce el fluido sobre la pared curvada es entonces:


� � �� � � � �1/2
2 2 1/2 5 2 3 2
F = FV + F H = 4. 992 × 10 + 899 × 10
F = 1. 028 3 × 106 lbf

Sabemos que la fuerza horizontal actúa a lo largo de z = 8 ft pero no sabemos a lo


largo de que valor de x actúa la fuerza vertical. Para determinarlo hay que calcular la
ubicación del centro de gravedad del fluido contenido en la zona que está por arriba de la
pared parabólica.
Observaremos a nuestra parábola como la función
� �1/2
z
x (z) = x0
z0
y viéndola de esa manera, lo que nos interesa es determinar la posición en x por la que
pasa la fuerza FV , que desde la nueva perspectiva corresponderı́a a la posición vertical. La
forma general de determinar la posición del centro de gravedad es:

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


40

�b
x 0
x f (x)dx
xhorizontal = �b
0
f (x)dx

1 b� �
2 0
f (x) 2 dx
xvertical = �b
0
f (x)dx

Para determinar la posición en x (del centro de gravedad) por la que pasa la fuerza FV ,
hay que determinar la xvertical :
� � � �1/2 �2

1 z0 � � 1 z0 z
f (z) 2 dz 2 0
x0 z0
dz
2 0
xvertical = � z0 = 2
0
f (z) dz xz
3 0 0
� � �z0 � �
1 2 z0 1 2 1 2 1 2 1 2
x
2z0 0 0
zdz x z
2z0 0 2 x z
2z0 0 2 0
0
xvertical = 2
= 2
= 2
xz
3 0 0
xz
3 0 0
xz
3 0 0
3
xvertical = x0
8
Entonces la fuerza vertical actúa a lo largo de la dirección x = 38 x0 = 38 (10) = 3. 75 ft.
Ası́, la fuerza tiene las componentes

F H = 899 × 103
FV = 4. 992 × 105

y pasa por el punto


(x1 , z1 ) = (3. 75, 8)

La fuerza actúa sobre la pared en un punto (centro de presiones) que resulta de prolon-
gar la dirección de la resultante fuerza hasta que se alcanza la pared. Este punto se puede
determinar determinando el cruce de la parábola que define a la pared con la lı́nea recta en
la que actúa la fuerza.
La ecuación de la lı́nea recta en la que actúa la fuerza la determinamos empleando la
pendiente del vector fuerza y el punto por el que pasa:

FV 4. 992 × 105
pendiente = m = − =− = −0,555 28
FH 899 × 103

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


41

y la ecuación de la recta es:

z − z1 = m (x − x1 )
z − 8 = −0,555 28 (x − 3. 75)
z = −0,555 28x + 10,0823

Entonces es necesario determinar el punto en donde nuestra recta intersecta a la parábo-


la � �2
x � x �2
z = z0 = 24 = 0,24x2
x0 10
Resolviendo el sistema de ecuaciones simultáneas se tiene que el punto de intersección
es (xCP , zCP ) = (5,42, 7,08)

Figura 2.14: Ubicación del centroide del área. Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

Figura 2.15: Ubicación del centro de presión. Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


42

2.6. Flotación y estabilidad


1. Un cuerpo sumergido en un fluido experimenta una fuerza de flotación igual al peso
del lı́quido que desplaza
2. Un cuerpo que flota desplaza su propio peso

F = ρgVdesplazado (2.24)

Figura 2.16: Cuerpo que flota en un fluido. Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

2.6.1. Estabilidad
Un cuerpo puede estar flotando en una posición que resulta inestable y en ese caso se
acomodará en la primera oportunidad.

La menor perturbación provocará que el cuerpo tome otra posición

Los ingenieros deben diseñar evitando inestabilidad en la flotación

Para determinar si un cuerpo flota de forma estable se debe perturbar y observar si


se presenta un momento que restaure el cuerpo a su posición original.

Para un cuerpo simétrico una forma de probar la estabilidad es:

1. Se calcula la posición inicial de flotación. Se calculan el centro de masa G y el centro


de empuje B.

2. Se inclina el cuerpo un pequeño ángulo Δθ y aparece ası́ una nueva lı́na de flotación.
Se calcula el nuevpo centro de flotación B� . La vertical trazada desde B� corta la
lı́nea de simetrı́a en el punto M, llamado metacentro, que es independiente de Δθ si
éste es pequeño.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


43

3. Si el punto M está por encima de G, esto es, si la altura metacéntrica es positiva,


aparecerá un momento restaurador y entonces decimos que la posición original es
estable. Si M está por debajo de G, el cuepo es inestable y volcará a la mı́nima per-
turbación. Cuanto mayor sea la distancia MG más estable será la posición original.

La altura metacéntrica es una propiedad de la sección transversal para un peso dado


y su valor es un indicador de la estabilidad del cuerpo.

2.7. Distribución de la presión en un fluido en movimien-


to como cuerpo rı́gido
Como ya se ha mencionado antes, bajo ciertas circunstancias un fluido puede moverse
como un cuerpo rı́gido en el que el movimiento relativo entre las partı́culas es nulo.
Cuando no existe movimiento relativo entre las partı́culas de fluido tampoco existen
fricciones internas y el efecto viscoso desaparece.
En esas condiciones, la conservación de la cantidad de movimiento establece un ba-
lance entre la presión, la gravedad y la aceleración de las partı́culas:

∇P = ρ (g − a) (2.25)

En el caso general de movimiento en el que existe traslación y rotación


Si el centro de rotación es O y la velocidad de traslación de ese punto es V0 , la veloci-
dad de un punto arbitrario del cuerpo es:

V = V 0 + Ω × r0 (2.26)

donde Ω es la velocidad angular y r0 es la posición del punto P.


Derivando se obtiene la aceleración del punto P:
d
a = V
dt
dV 0 � � dΩ
a = + Ω × Ω × r0 + × r0 (2.27)
dt
���� ���������������������� dt
����������
Aceleración Aceleración Aceleración lineal
lineal centrı́peta debida a cambios en Ω
Ejemplo Un tanque de 1m de profundidad se encuentra en caı́da libre. Despreciando
la fricción con el aire, calcule la presión en el fondo del tanque si pa = 101 kPa.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


44

2.7.1. Aceleración uniforme


Cuando la aceleración no involucra movimiento de rotación, solo existe aceleración
lineal y entonces de acuerdo a
∇P = ρ (g − a) (2.28)

solo es necesario llevar a cabo la suma vectorial de g y −a para determinar la dirección del
gradiente de la presión
El gradiente de la presión es la dirección en la que los cambios de la presión son más
intensos
Las superficies de presión constante son perpendiculares a la dirección del gradiente

Figura 2.17: Superficie libre en un fluido que se mueve con aceleración constante. Tomada de White, F. M.
Fluid Mechanics. Mc. Graw-Hill.

De acuerdo a la figura, el ángulo que forman la horizontal y las superficies de presión


constante es: � �
−1 ax
θ = tan (2.29)
g + az
donde a x y az son las componentes horizontal y vertical de a, respectivamente.
Evidentemente, una de estas lı́neas inclinadas corresponde a la lı́nea de la superficie
libre. Para determinarla hay que satisfacer la condición de que el lı́quido bajo la superficie
libre se conserva aun después de que la superficie libre se inclina, siempre que el lı́quido
no se derrame.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


45

La variación máxima en la presión la encontraremos si nos movemos en la dirección


del gradiente y está dada por:
dP
= |∇P| = |ρ (g − a)| (2.30)
ds
donde |∇P| es el módulo del gradiente, y para el caso mostrado en la figura:
� �1/2
|∇P| = ρ a2x + (gz + az )2 (2.31)

El cambio en la presión lo podemos determinar integrando la variación de la presión


en la dirección del gradiente:
� s2
P= |∇P| ds = |∇P| Δs (2.32)
s1

Ejemplo Una taza de café está sujeta a una aceleración de 7 m/s2 . La taza tiene una
profundidad de 10 cm, un diámetro de 6 cm, y la superficie libre del café en reposo se
encuentra a 7 cm. (a) Suponiendo un movimiento de cuerpo rı́gido del café, determine si
el café se derrama. (b) calcule la presión en el punto A cuando la densidad del café es
1010 kg/m3 .

Figura 2.18: Superficie libre en una taza que se mueve con aceleración constante. Tomada de White, F. M.
Fluid Mechanics. Mc. Graw-Hill.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


Capı́tulo 3

Relaciones integrales para un volumen


de control

3.1. Layes fı́sicas fundamentales de la mecánica de fluidos


Para el análisis de los sistemas que involucran el flujo de fluidos, en general se hace
uso de alguna de las siguientes técnicas:

1. Volumen de control o análisis integral

2. Análisis diferencial

3. Análisis experimental

¡Los tres son iguales en importancia!.


El análisis diferencial se desarrolló inicialmente por Euler y Lagrange en el siglo
XVIII. Conduce a la solución de las ecuaciones de Navier-Stokes, ecs. diferenciales en
derivadas parciales (x, y, z, t), no lineales. La complejidad matemática y el alcance de las
computadoras limitan este análisis.
El análisis experimental es muy costoso es todos los sentidos, tiempo, dinero y esfuer-
zo. Los resultados son limitados a condiciones particulares, pero resulta completamente
necesario cuando se plantea una nueva solución o teorı́a.
El análisis con volúmenes de control fue planteado originalmente por Euler pero fue
hasta los años 40 que se desarrolló de manera formal. El análisis con volúmenes de control
es el método más útil para fines de uso práctico. éste da respuestas ingenieriles, algunas
veces “burdas” pero siempre útiles.

46
47

3.1.1. Flujo volumétrico y flujo másico


Para realizar los análisis se requiere determinar el flujo de masa o de volumen que
atraviesa una superficie imaginaria en el flujo.
¿Cómo determinamos el volumen de fluido que atraviesa la superficie S ?
En general la velocidad del flujo V varı́a localmente, entonces debemos integrar V
sobre el área elemental dA.
Adicionalmente, el flujo atraviesa el área con un ángulo θ respecto a la dirección nor-
mal a la superficie.
Si n es el vector normal a dA, la proyección del vector de velocidad V en la dirección
del vector n es
� �
V · n = ��V �� |n| cos θ (3.1)

entonces, el volumen de fluido que atravieza dA en un tiempo dt es


� �
dV = V · n dA dt
dV � �
= dQ = V · n dA (3.2)
dt
dQ : flujo volumétrico o gasto a través de dA

Entonces, el flujo volumétrico que atraviesa la superficie S es:


� � � �
Q= V · n dA = Vn dA (3.3)
S S

donde Vn es la componente de V en la dirección de n.


Es importante mencionar que el uso del producto punto permite distinguir entre flujos
de entrada y salida.
Como n es el vector unitario que apunta hacia afuera del volumen de control, cuando
� � � �
V · n es positivo se trata de una salida y cuando V · n es negativo se trata de una entrada.
Para determinar el flujo másico se multiplica el flujo volumétrico por la densidad:
� � � �
·
m= ρ V · n dA = ρVn dA (3.4)
S S

Si la densidad es constante:
·
m = ρQ (3.5)

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


48

3.2. Teorema de transporte de Reynolds


Si B una propiedad cualquiera del fluido, por ejemplo energı́a, cantidad de movimiento,
etc., y sea
dB
β= , (3.6)
dm
el valor intensivo de B (cantidad de B por unidad de masa m). Entonces,

dB = β dm,

además, sabemos que la masa está relacionada con la densidad como sigue:
dm
= ρ ⇒ dm = ρ dV,
dV
es decir que:
dB = ρβ dV, (3.7)

por lo tanto, la cantidad total de B en un volumen de control con volumen V es:


� �
BVC = dB = ρβ dV. (3.8)
VC VC

Ahora, las variaciones de cualquier propiedad B de un sistema abierto que es encerrado


por un volumen de control VC está dado por:
d d �� � �� �
(Bsis ) = (BVC ) + flujos de salida de B − flujos de entrada de B .
dt dt
������������
Variación temporal
de B en el VC
(3.9)
Este es el teorema de transporte de Reynolds.
En el caso unidimensional se tendrı́a que:
�� � � �
d d
(Bsis ) = ρβ dV + (ρβAV) sal − (ρβAV)ent .
dt dt VC

En el caso general multidimensional cada elemento diferencial dA tendrá una veloci-


dad diferente V que formará un ángulo θ con el vector unitario normal a dA que llamare-
mos n.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


49

Recuérdese que:
V · n = |V| |n| cos θ. (3.10)
Si definimos a n de tal manera que apunte hacia el exterior del VC:

V · n = Vn en flujo saliente,
V · n = −Vn en flujo entrante.

Ası́, los términos de flujo se pueden representar como integrales simples que involu-
cran a V · n, como sigue:
�� � � �

flujos de B = ρβ V · n dA.
SC

Además, si consideramos la posibilidad de que el volumen de control se mueva con


una velocidad constante V s , el efecto dinámico resultante es que los flujos no son propor-
cionales a V sino a Vr que es la velocidad relativa entre el flujo y el volumen de control.
Entonces:

Vr = V − V s , (3.11)
V : Velocidad del flujo,
Vs : Velocidad del VC,
Vr : Velocidad relativa.

La forma general del teorema de transporte de Reynolds es:


�� � � � �
d d
(Bsis ) = ρβ dV + ρβ V r · n dA. (3.12)
dt dt VC SC

3.2.1. Aproximaciones unidimensionales al término de flujo


Cuando los flujos son aproximadamente unidimensionales (las propiedades del flujo
son casi constantes en la secciones transversal de entrada y salida):
� � � �� ��
� �
ρβ V · n dA = ρi βi Ai Vr,i sal − ρi βi Ai Vr,i ent . (3.13)
SC

Por ejemplo, observando la figura ?? vemos que hay entradas de flujo en las secciones
1 y 4, y salidas en las secciones 2, 3 y 5. En este caso particular:
� � �
ρβ V r · n dA = ρ2 β2 A2 Vr,2 + ρ3 β3 A3 Vr,3 + ρ5 β5 A5 Vr,5 − ρ1 β1 A1 Vr,1 − ρ4 β4 A4 Vr,4 .
SC

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


50

Figura 3.1: Flujos unidimensionales a través de un volumen de control. Tomada de White, F. M. Fluid
Mechanics. Mc. Graw-Hill.

3.3. Conservación de masa. Ecuación de Continuidad


La conservación de masa en un sistema establece que en el sistema la masa permanece
constante:

m sis = cte,
dm sis
⇒ = 0. (3.14)
dt
El teorema de transporte de Reynolds está definido para una cierta propiedad B. Pode-
mos aplicar ese teorema a las ecuaciones de conservación para obtener ecuaciones inte-
grales aplicables a fluidos.
En el caso de la ecuación de conservación de masa:

B = m,
m : masa.

De la definición de β, (3.6), tenemos:


dB dm
β= = = 1.
dm dm
Entonces, según (3.12) y (3.14), la variación de la masa en un volumen de control está
dada como: �� � �
d d � �
(m sis ) = 0 = ρ dV + ρ V r · n dA. (3.15)
dt dt VC SC

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


51

Cuando el volumen de control está fijo:


� � � �

ρ dV + ρ V · n dA = 0. (3.16)
∂t VC SC

Cuando el VC tiene sólo un cierto número de entradas y salidas:


� �� ��
∂ � �
ρdV + ρi Ai Vr,i sal − ρi Ai Vr,i ent = 0. (3.17)
∂t VC
Cuando el flujo es permanente ∂/∂t = 0, entonces:
� � �
ρ V · n dA = 0. (3.18)
SC

La ec. (3.18) indica que en estado permanente los flujos de entradas y salidas se com-
pensan. Además, si las salidas son unidimensionales:
�� � �� �
ρi Ai Vr,i sal = ρi Ai Vr,i ent .

En nuestro ejemplo anterior, en caso de estado permanente la ecuación de continuidad


es:
ρ2 A2 V2 + ρ3 A3 V3 + ρ5 A5 V5 = ρ1 A1 V1 + ρ4 A4 V4 .

Notar que:
·
ρAV = m,
ρAV : flujo másico a través de una cierta sección.

Por otra parte, si el fluido es incompresible, ρ ≈ cte, entonces de la (3.16)


� � �
ρ V · n dA = 0, (3.19)
�S C � �
∴ V · n dA = 0, (3.20)
SC

y si además las entradas y salidas son unidimensionales:


�� � �� �
Ai Vr,i sal = Ai Vr,i ent , o bien,
� �
Q sal = Qent , (3.21)
Q : Flujo volumétrico o caudal.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


52

Ejemplo A través de un tubo con diámetro de 10 cm circula agua a 50 lt/min. Si la


tuberı́a tiene una reducción para conducir el agua en un tubo con diámetro 2 cm, ¿a qué
velocidad circula en esta tuberı́a?
Ejemplo El tanque de la figura se encuentra abierto al ambiente. La altura del agua
dentro del tanque es h. (a) Encuentre una expresión para calcular el cambio instantáneo
de h. (b) Suponga que el área transversal del tanque es At = 0.3 m2 . Si en el tiempo t = 0
se midió que h = 1 m, si Q1 = 0,001 m3 /s ¿cuánto deben valer Q2 para que después de
10 min te tenga que h = 70 cm. (c) Si el diámetro de la entrada 1 es 2 cm y en 2 es 4 cm,
¿cuánto valen las velocidades de entrada y salida?

Figura 3.2: Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

3.4. Conservación de la cantidad de movimiento


Aplicaremos el teorema de transporte de Reynolds con

B = mV : cantidad de movimiento,

y en este caso de la ec. (3.6)

dB d � �
β= = mV = V,
dm dm
y de la ec. (3.12) es:
�� � � � �
d d
(Bsis ) = ρβ dV + ρβ V r · n dA,
dt dt VC SC

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


53

entonces: �� � �
d � � d � �
mV = ρV dV + ρV V r · n dA. (3.22)
dt dt VC SC
La ecuación de conservación de cantidad de movimiento es:
� d � �
F = mV ,
dt �� � �
� d � �
∴ F = ρV dV + ρV V r · n dA. (3.23)
dt VC SC

Observaciones importantes:

1. La velocidad está referida a un sistema inercial (sin aceleración).

2. La resultante de las velocidades debe incluir todas las fuerzas que actúan sobre la
masa dentro del volumen de control, incluyendo peso.

3. La ecuación es vectorial, por lo tanto son debemos tener una ecuación por cada
coordenada considerada en el análisis.

3.4.1. Aproximaciones unidimensionales al término de flujo de canti-


dad de movimiento
Cuando los flujos son aproximadamente unidimensionales (las propiedades del flujo
son casi constantes en la secciones transversal de entrada y salida):
� � � �� ��
� �
ρV V · n dA = ρi Vn,i Ai Vr,i sal − ρi Vn,i Ai Vr,i ent , (3.24)
SC
Vn,i : Velocidad normal a la superficie de salida o entrada, (3.25)

por lo tanto:
·
mi ≈ ρi Ai Vn,i . (3.26)
Es decir que la ec. de la cantidad de movimiento con entradas y salidas unidimensio-
nales queda como
� �� � �� · � �� · �
d
F= ρV dV + mi Vr,i − mi Vr,i .
dt VC
sal ent

Ejemplo Considere el conducto de la figura. Suponiendo flujos de entrada y salida


uniformes y estado permanente, encuentre una expresión para determinar la fuerza total
sobre el volumen de control.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


54

Figura 3.3: Cambio del vector de velocidad de un fluido en un conducto. Tomada de White, F. M. Fluid
Mechanics. Mc. Graw-Hill.

Ejemplo Suponga que un chorro fijo de agua es desviado un ángulo θ por una superficie
como lo muestra la figura. El área transversal del chorro es A y se supone que al ser
desviado el chorro no cambia la magnitud de su velocidad ni su área transversal. El chorro
fluye en estado permanente y la presión es constante en todos lados. (a) Encuentre una
expresión para las componentes F x y Fy de la fuerza que se ejerce sobre la superficie. (b)
¿Para que valor de θ se ejerce la máxima fuerza? (c) Grafique el comportamiento de F.

Figura 3.4: Cambio del vector de velocidad. Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

Ejemplo Un chorro de agua con velocidad V j choca de manera normal a una placa
plana que se mueve con una velocidad Vc , como lo muestra la figura. Encuentre la fuerza
necesaria para hacer que la placa se mueva con velocidad constante. El área del chorro es 3

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


55

Figura 3.5: Fuerza como función del ángulo que se desvı́a el flujo. Tomada de White, F. M. Fluid Mechanics.
Mc. Graw-Hill.

cm2 , V j = 20 m/s y Vc = 15 m/s. Desprecie el peso del chorro, suponga estado permanente
respecto a la placa móvil y que el chorro se divide en dos partes iguales.

Figura 3.6: Chorro que de desvı́a en la dirección vertical. Tomada de White, F. M. Fluid Mechanics. Mc.
Graw-Hill.

3.5. Ecuación de la Energı́a


En este caso aplicaremos el teorema de transporte de Reynolds a la primera ley de la
termodinámica. En este caso:

B = E,
E : Energı́a,

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


56

y la energı́a por unidad de masa es:


dE
β= = e. (3.27)
dm
La primera ley de la termodinámica es:
dQ dW dE
− = , (3.28)
dt dt
���� ���� dt
· ·
Q W

donde Q es el calor. El teorema de Reynolds, ec. (3.12):


�� � � � �
d d
(Bsis ) = ρβ dV + ρβ V r · n dA,
dt dt VC
�� S C � � � �
d d
⇒ (E sis ) = ρe dV + ρe V r · n dA,
dt dt VC SC

es decir que: �� � �
· d · � �
Q−W = ρe dV + ρe V r · n dA. (3.29)
dt VC SC

La energı́a total acumulada en un cuerpo por unidad de masa:


1
u + V 2 + gz,
e = � (3.30)
2
u : energı́a interna molecular (función de P y T en una sustancia pura),

V : módulo de la velocidad,
z : coordenada vertical.

El trabajo puede ser de varios tipos también:


· · · ·
W = W eje + W presión + W viscoso ,
·
W eje : Trabajo entregado o recibido a través de un eje mecánico,
·
W presión : Trabajo hecho por la presión en la superficie del VC,
·
W viscoso : Trabajo hecho por los esfuerzos cortantes,

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


57

considerando el trabajo hecho por unidad de tiempo (potencia), tenemos que:


· � �
dW presión = − (PdA) Vn,in = −P −V · n dA,
� � �
·
W presión = P V · n dA,
SC
·
dW viscoso = −dA� τ · V,
·
W viscoso = − τ · V dA,
SC
� � � �
· ·
∴ W = W eje + P V · n dA − τ · V dA,
SC SC

entonces la ecuación de la energı́a, ecuación (3.29), es:


� � � � �� � � � �
· · d
Q − W eje − P V · n dA + τ · V dA = ρe dV + ρe V · n dA,
SC SC dt VC
� �� � �S C � �
· · d P � �
Q − W eje + τ · V dA = ρe dV + ρ e+ V · n dA,
SC dt VC SC ρ

ahora, incluyendo los términos que incluye la energı́a, ecuación (3.30):


� �� � � � � � �
· · d 1 2 1 2 P � �
Q−W eje + τ·V dA = ρ �u + V + gz dV + ρ �u + V + gz + V · n dA.
SC dt VC 2 SC 2 ρ

Sólo estamos considerando flujos de trabajo, y en el caso de entradas y salidas el



trabajo viscoso es muy pequeño y en este análisis se puede despreciar ( S C τ · V dA ≈ 0),
por lo que la ecuación de la energı́a es:
 
 
�� � � � �   � �
· · d 1 2  P 1 2 
Q − W eje = ρ �u + V + gz dV + ρ � u + + V + gz V · n dA.
dt VC 2 SC  ρ 2 
���� 

h

Además, también podrı́amos utilizar la definición de entalpı́a (h = �


u + P/ρ) en nuestra
ecuación de energı́a.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


58

3.5.1. Aproximación unidimensional del término de flujo de cantidad


de energı́a
Cuando los flujos son aproximadamente unidimensionales (las propiedades del flujo
son casi constantes en la secciones transversal de entrada y salida):
� � �� � �� �
P 1 2 P 1 2 ·
ρ �u + + V + gz V · n dA = u + + V + gz m sal ,

SC ρ 2 ρ 2 sal
�� P 1 2

·
− u + + V + gz ment , (3.31)

ρ 2 ent

donde:
·
mi ≈ ρi Ai Vn,i . (3.32)

Es decir que la ec. de la energı́a con entradas y salidas unidimensionales queda como
�� � � � �� �
· · d 1 2 P 1 2 ·
Q − W eje = ρ �u + V + gz dV + u + + V + gz m sal

dt VC 2 ρ 2 sal
� � �
P 1 ·
− �u + + V 2 + gz ment . (3.33)
ρ 2 ent

3.5.2. Ecuación de Bernoulli en estado permanente


En realidad existe una ecuación de Bernoulli que es válida para analizar estado transi-
torio (para detalles consultar libro: White, F. M. Fluid Mechanics, McGraw Hill). En esta
clase nos enfocamos en problemas en estado permanente.
Retomando nuestra ecuación (3.33), y suponiendo estado permanente (d ( ) /dx = 0),
· ·
sin transferencia de calor (Q = 0), y sin trabajo de eje (W eje = 0):
�� P 1 2

·
�� P 1 2

·
u + + V + gz ment −
� u + + V + gz m sal = 0.
� (3.34)
ρ 2 ent ρ 2 sal

Ahora, definimos un volumen de control con una entrada y una salida, digamos condi-
ción 1 y 2, y como suponemos estado permanente, de la ecuación de continuidad:
· ·
ment = m sal ,

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


59

y si la temperatura es constante (�u = cte ⇒ �


uent = �
u sal ):
� � � �
P 1 2 P 1 2
�u + + V + gz − � u + + V + gz = 0,
ρ 2 ent ρ 2 sal
P1 1 2 P2 1 2
+ V1 + gz1 = + V2 + gz2 = cte,
ρ 2 ρ 2
Ecuación de Bernoulli.

Las limitaciones para la aplicación de la ecuación de Bernoulli son todas las condicio-
nes que hemos supuesto en su obtención, es decir:

1. Flujo permanente

2. Flujo incompresible

3. Flujo sin fricción viscosa

4. Flujo a lo largo de una lı́nea de corriente

5. Sin flujo de trabajo entre 1 y 2 (sin bombas ni turbinas entre 1 y 2).

6. Sin flujo de calor entre las condiciones 1 y 2.

Ejemplo Determine una expresión para calcular el flujo másico a través de un tubo
venturi en función de la caı́da de presión ΔP = P2 − P1 . Si el flujo es de agua, y la tuberı́a
tiene un diámetro D1 = 0,16 m, D2 = D1 /2 m y se mide que ΔP = −8 kPa, ¿cuánto vale
el gasto másico?

Figura 3.7: Tubo venturi. Tomada de White, F. M. Fluid Mechanics. Mc. Graw-Hill.

Ejemplo Determine una expresión que permita calcular la velocidad de salida del fluido
a través de una tobera ubicada una altura h por debajo del nivel del agua en el tanque
abierto a la atmósfera. Suponga que la capacidad del tanque es muy grande.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL


60

Figura 3.8: Chorro de agua debido a la presión hidrostática en un depósito. Tomada de White, F. M. Fluid
Mechanics. Mc. Graw-Hill.

Grupo de Energı́as Térmica y Renovable, Facultad de Ingenierı́a Mecánica y Eléctrica, UANL

También podría gustarte