Libro Mecanica
Libro Mecanica
Libro Mecanica
1.1 Introducción
Incluso si se confirmara la existencia de la ”materia oscura”, el estado de
agregación más abundante en el universo serı́a el fluido, puesto que la materia
en forma de plasma estelar generalización del estado fluido es la más frecuente
del cosmos.
De los 6.1 x 1024 kg de masa que se atribuyen a nuestro planeta, casi una ter-
cera parte se encuentra en forma fluida en su núcleo denso de hierro y nı́quel fun-
didos. En nuestro entorno terrestre mas próximo se comprueba que los océanos
ocupan tres cuartas partes de la superficie de la tierra, totalizando la masa de
agua en el planeta 1.4 x 1021 kg. La atmósfera se extiende centenares de kms
sobre nuestras cabezas, aunque la práctica totalidad de sus 5.1 x 1018 kg de
masa se encuentra en los primeros 17 km.
Descendiendo a las dimensiones humanas, encontramos que en la mayorı́a de
los procesos geológicos, biológicos y en la práctica totalidad de las realizaciones
tecnológicas los flujos de fluidos intervienen en alguna medida. En el capı́tulo
I.2 de esta Memoria se da cuenta de abundantes ejemplos.
Dentro del campo especı́fico de la técnica se constata además la abundancia
de dispositivos mecánicos que intercambian energı́a con fluidos en movimiento:
bombas, compresores, turbinas hidráulicas o de gas, bombas o motores de de-
splazamiento positivo, accionamientos, hidroneumáticos, etc...
Nuestra sociedad, eminentemente tecnológica, se caracteriza por una con-
tinua realimentación entre las disciplinas técnicas y cientı́ficas y una mutua
potenciación de las diversas técnicas entre si. Las materias que nos ocupan no
son una excepción:
El rápido desarrollo de las técnicas Informáticas y de Materiales han impul-
sado decididamente en los últimos años el avance del conocimiento de los flujos
y el diseño de Turbomaquinaria, y viceversa.
Los procesos de moldeo de plásticos, aleaciones o extrusión de nuevos mate-
riales precisa de la descripción detallada de los flujos involucrados.
Por consiguiente, no cabe ninguna duda del interés de explicar y predecir el
comportamiento dinámico de los fluidos, tanto en sus aspectos fı́sicos fundamen-
tales, contenidos en la Mecánica de los Fluidos, como en sus aplicaciones en
la ingenierı́a, caso de la Ingenierı́a Fluidodinámica.
1
1.2 Concepto de fluido
Aproximación intuitiva
2
Observaciones
3
• Existe un volumen limite δϑ* bajo el cuál las variaciones moleculares son
importantes y sobre el cual puede haber variaciones importantes en la
agregación de las moléculas
δm
ρ= lim (1)
δϑ→δϑ∗ δϑ
δϑ*≈10−9 mm3 para lı́quidos y gases a presión atmosférica. e.g. aire condi-
ciones estándar ⇒ 3x107 moléculas son suficientes para definir una densidad
casi constante de acuerdo a (1).
4
1.4 Dimensiones y unidades
Dimensión: Medida por la cual una variable fı́sica se expresa cuantitativa-
mente. Ej. Longitud expresa distancia, desplazamiento, anchura ...
Sistemas de Unidades
F L t T M
Europeo Kilopondio (kp) Metro (m) Segundo (s) Kelvin (◦ K) U.T.M.
Británico (BG) Libra fuerza (lbf) Pie (ft) Segundo (s) Rankine (◦ R) slug
Algunas conversiones son:
1 ft = 0.3048 m 1 lbf = 4.4482 N
1s=1s 1 slug = 14.5939 kg
1 o K = 1.8 o R
Ecuación de Dimensiones
5
Homogeneidad dimensional de las ecuaciones
6
1.5 Propiedades termodinámicas de un fluido
• Presión
• Temperatura
• Densidad
• Energı́a Interna
Trabajo requerido para llevar el sistema de masa m desde el origen a una cierta
posición ~r en contra del campo gravitatorio (−m~g · ~r).
• Energı́a cinética
1
Energı́a por unidad de masa, e = û + 2 V2 - (~g · ~r).
Gases a altas temperaturas y bajas presiones (lejos del punto crı́tico) verifican
la ley de los gases perfectos
7
p=ρRT R = cp - cv = constante de los gases perfectos
R = Λ/Mg Λ= 8314 m2 / (s2 ·o K) = 49700 ft2 / (s2 ·o R)
donde Λ es una constante universal, Mg es la masa molecular del gas y cp , cv
son los calores especı́ficos a presión y volumen constante, respectivamente.
du
τ =µ
dy
du
es la velocidad de deformaciòn, τ el esfuerzo cortante y µ es el coeficiente
dy
de viscosidad dinámica. Observe que esta es una ecuación lineal con que pasa
por el origen y cuya pendiente es µ.
8
8.8x10−6 kg/(m·s) → Hidrógeno
µ=
1.5 kg/m·s) → Glicerina
Gases: µ ↑ con T
µ T0 T2
ln ≈ a + b + c 02
µ0 T T
9
• Conductividad Térmica
k = k (p , T) [k] = julios/(s m o K)
• Tensión Superficial
10
Entrefases más comunes: agua - aire, Hg - aire
Si θ < π/2 se dice que el lı́quido Si θ > π/2 se dice que el lı́quido
moja el sólido no moja el sólido
(Agua-vidrio θ ≈ 0o ) (mercurio-vidrio θ ≈ 130o )
• Presión de vapor
En un flujo acotado por una superficie sólida, las interacciones moleculares cau-
san que el fluido en contacto con la superficie busque encontrar un equilibrio
de momento y energı́a. El equilibrio se da en lı́quidos y gases en condiciones
normales y en el punto de contacto
Vf = Vwall y Tf = Twall
11
• Velocidad del sonido
agasideal = (k R T)1/2
12
Problemas Resueltos
Problema 1
Solución
(a)
g 1kg 100cm kg kg
µ = 0.01 = 0.001 = 10−3
cm · s 1000g 1m m·s m·s
(b)
kg 1slug 0.3048m slug
µ = 10−3 = 2.09 · 10−5
m · s 14.59kg 1f t ft · s
Problema 2
Solución
13
Problema 3
Solución
−1 1.49 1.49
1. L t = L2/3
[1] ⇒ = L1/3 t−1 ⇒ [1.49] = L1/3 t−1
n n
es decir, la fórmula no es dimensionalmente consistente pues 1.49 deberı́a ser un
número puro
Solución
14
Figura Problema 4
Problema 5
Figura Problema 5
Solución
15
ferencia.
h: altura de capilaridad
Demuetre:
1. La velocidad disminuirá exponencialmente con el tiempo.
2. La particula se detendrá después de haber viajado una distancia
x = mV0 /K.
Solución
16
mV0
De este modo, V decae exponencialmente con t, y cuando t → ∞, x =
K
Problema 7
Solución
p = (32 lbf/in2 ) (144 in2 /f t2 ) + 2116 lbf/f t2 = 4608 +2116 ≈ 6724 lbf/f t2
p 6724f t2
ρaire = = = 0.00732 slug/f t3
RT (1717f t · lbf /(slug ·◦ R))(535◦ R)
Por lo que el peso total del aire en la rueda es:
Solución
17
r s
8σ 8(0.073N/m)
Cálcule dmax = = ≈ 0.00156 m ≈ 1.6 mm
πρg π(998kg/m3 )(9.81m/s2
Respuesta (b)
18
Problemas Propuestos
Problema 1
Problema 2
Problema 3
Problema 4
¿Qué diámetro de tubo de vidrio limpio se requiere para que el ascenso capilar
de agua a 20◦ C en el tubo sea menor que 1.0 mm?
Problema 5
Un avión de propulsión a chorro vuela a 550 mph a una altitud de 35000 ft,
donde la temperatura es de -66 ◦ F. Determinar la razón de la velocidad del
avión, V , a la velocidad del sonido, a, a la altura especificada (número de
Mach). Suponer k = 1.40.
19
Capı́tulo 2
Fluidostática
−∇p + ρ~g = 0
Esta ecuación es una ecuación diferencial en derivadas parciales (EDP) lineal
para la presión que posee un carácter vectorial, es decir, equivale a tres ecua-
ciones escalares.
20
Figure 2: Sistema de referencia situado en la superficie libre donde la gravedad
está dirigida hacia abajo
Consecuencias:
21
Z2 Z2
dp = − ρ g dz → p2 − p1 = −ρg(z2 − z1 )
1 1
∆p
∆p = −ρg∆z → ∆z = −
ρg
Consecuencias:
B∼
= 0.0065 K/m z ∈ [0, 1.1 × 103 ]; T0 = 15◦ C
g/RB
Bz g
p = pa 1 − = 5.26
To RB
22
tomando z = 0 en el nivel del mar y p(z = 0) = patm .
2.3 Manómetros
∆p
Ya que en un lı́quido ∆z = − , un cambio de altura de un lı́quido equivale a
ρg
un cambio de presión, una columna estática de uno o más lı́quidos puede usarse
para medir diferencias de presión entre dos puntos
Manómetros
Para resolver una columna estática de fluidos incompresibles se deben aplicar
adecuadamente las tres reglas siguientes:
Ejemplo 2.1:
Solución:
23
pA − pB = (pA − pC ) + (pC − pD ) + (pD − pE ) + (pE − pF ) + (pF − pB )
d1 = 10” , d2 = 3” , d3 = 4” , d4 = 5” , d5 = 8”
Superficies planas
24
centro de presiones CP.
Ejemplo 2.2:
Solución:
25
Figure 6: Compuerta del Ejemplo 2.2
26
P
c) Bz = Fz = F cos θ = 30700 lbf ( Fz = 0)
X
Fx = 0 = Bx + Fx − P ⇒ Bx = P − F sin θ = 6300 lbf ( Fx = 0)
Superficies Curvas
Principio de Arquı́medes
R R
∇p = ρ~g ⇒ V
∇p dV = V
ρ~g dV = ρ~g V
⇒ F~e = −ρ~g V
pd~s = −F~e
R R
∇p dV =
V S
Cuerpos Flotantes: Solo una porción del cuerpo esta sumergida (figura 8).
En este caso
Fe = ρgVdisp = ρb gVb
27
Figure 8: Cuerpo flotante. Observe que el peso actúa en CG y la flotación en
CF (o B).
y la fuerza de empuje equilibra el peso del cuerpo. Además, ambas fuerzas son
colineales pues la condición de equilibrio implica que no hay momento de giro
neto. En el caso en que Fe = ρb gVb se habla de cuerpos flotantes neutros los
cuales permanecerán en equilibrio en cualquier punto en el que se encuentren
depositados en el fluido (figura 9).
En esta posición, el cuerpo se inclina un ángulo pequeño ∆θ, con lo cual cambia
28
el centro de flotación a B’. La intersección de la lı́nea vertical trazada desde
B’ con el eje de simetrı́a del cuerpo nos proporciona un punto M denominado
metacentro, el cual es independiente de ∆θ para ángulos pequeños.
Figure 10: Ilustración de las fuerzas que actúan sobre un cuerpo flotante cuando
se perturba su posición de equilibrio.
29
Cálculo de la estabilidad de flotación
Este cálculo hace uso del concepto ”waterline area” o área subtendida por la
lı́nea de flotación. Ésta es el área sombreada en azul en la figura 11.
Ejemplo 2.3:
Solución:
30
Figure 13: Esquema de la barcaza del ejemplo
IO b(2L)3 /12 H L2 H
MG = − GB = − = −
Vsub 2LbH 2 3H 2
luego la barcaza será estable si M G > 0 → L2 > 3H 2 /2 o L > 1.225H.
31
Problemas Resueltos
Problema 1
Solución
Aire → p ≈ pA
AS = 3 m × 1 m = 3 m2
32
pB − pA 8000
1. pB = pA = γoil h → γoil = = = 8000 N/m3
h 1
2. pC = pD + γagua · hagua pD = pB + γoil · hoil
pC = pB + γoil · hoil + γagua · hagua = 58000 + 8000 · 1 + 9806 · 1 = 75806 P a
3. FH = 0 por simetrı́a
4. FV = pC · AS + γagua · vol
V ol =
1
V1 = V3 = 0.5 × 0.7 × 1 m3 , V2 = 0.7 × 2 × 1 m3
2
FV = 75806 × 3 + 9806 [(3 × 1.5 × 1) − (2 × 0.5 × 0.5 × 0.7 × 1) − (0.7 × 2 × 1)]
FV = 254389.5 N
33
Problema 2
Solución
1. R = γagua · ȳ · πr2
Ixx πr4 /4 r2
2. yCP = − sin φ 2
= − sin φ 2
= − sin φ
ȳ · πr ȳ · πr 4ȳ
r2 πr4
3. M = |R| · |yCP | sin 90◦ = γagua ȳπr 2
· sin φ = γagua
4ȳ 4
Independiente de ȳ
34
r2
D = ȳ + |yCP | sin φ = ȳ + sin2 φ
4ȳ
35
Problema 3
2. La presión en el punto C.
3. La presión en el punto B.
4. La altura h00 .
Nota: γagua = 9800 N/m3 , γaceite = 8000 N/m3 , h = 100 cm, h0 = 50 cm,
d = 20 cm, M = 1 kg.
Solución
Mg
1. p0 = π = 312 P a
d2
4
2. pC − patm = +γH2 O · h = 9800 P a
3. pB − pC = −γH2 O · h0 = −4900 P a ⇒ pB = pC − 4900 = 4900 P a
pB − p0 4900 − 312
4. pB − p0 = γaceite · h00 ⇒ h00 = = ≈ 0.57 m
γaceite 8000
36
Problema 4
Solución
b = L = 30 cm
37
1. FP = γH2 O h · L · b → FP ≈ 869 N
1 3 √
π 12 bL 2 1 L2
2. yCP Agua = − sin =− → yCP ≈ 5.4 mm
4 hb L 2 12 h
3. FAire = pAire L · b → FAire ≈ 900 N
4. yCP Aire = 0 Pues la presión del aire es uniforme y FAire actúa en el CG
5. MAire B = MH2 O B ⇒ FAire · `Aire = FP · `CP Ambas fuerzas actúan
normal a la válvula
√ !
L L L 2L
`Aire = ; `CP = + |yCP Agua | = 1+
2 2 2 12 h
2 3
√ √ !
L
b 1 2 h
L b 2 2
pAire = γH2 O h
L b +γH2 O
⇒ pAire = γH2 O h + L
2
2
h 2 12 12
√
pAire 2
h= −L ≈ 0.985 m
γH2 O 12
38
Problema 5
Solución
pA = pair + γoil z
pair = pA − γoil z
pA y pair deben convertirse a psf desde psi multiplicando por 144, para
tener unidades BG consistentes
pair = 2 × 144 − 54 × 3.5 = 99 psf → pair = 0.6175 psi
pA + γoil (2 f t) − SG1 γH2 O h = 0
pA + γoil (2 f t)
SG1 = → SG1 = 4.23
γH2 O h
39
Problema 6
La figura muestra una compuerta rectangular AB que separa dos fluidos distin-
tos. La compuerta tiene una anchura b y puede girar en una bisagra en el punto
B. La columna izquierda está limitada por una cavidad con aire presurizado
pair , mientras que la derecha se encuentra abierta a la atmósfera. En estas
condiciones determine:
Solución
40
Definiendo hCG i = hA + 16 cos 30
1 FA
γA = − pair γA = 311.04 lbf /f t3
hCG i A
b) 60◦ forma la compuerta con la horizontal
Ixx 1 3
yCP i = −γA sin 60 Ixx = bAB = 4 × 122 = 576 f t4
FA 12
yCP i = −0.31 f t
1 FB
FB = pCG d A ; hB = − 16 cos 30◦ = 55.64 f t
γH2 O SGB A
Ixx
d) yCP d = −SGB γH2 O sin 60 = −0.19 f t
FB
e) FT = FA − FB ya que van en sentidos opuestos
41
Problema 7
1. La presión en el aire pA .
2. La presión en la entrefase agua-mercurio. (Punto N).
Solución
1. pN = pM
pN = pA + SGoil γH2 O 3 m + γH2 O (4.6 − 3)
pN = pA + 0.8 × 9806 × 3 + 9806 × 1.6 = PA + 39224
pN = γHg × 0.3 = 13.6 × 9806 × 0.3 = 40008.48 P a (manométrica).
pN = pM ⇒ pA +39224 = 40008.48 → pA = 784.48 P a (manométrica)
2. pN = 40008.48 P a manométrica
42
Problema 8
Solución
1. pM = pN
8.5 7 3
pM = patm + 62.4 − 62.4 − 0.9 (62.4) = patm − 6.24
12 12 12
7 10
pN = patm + SGA 62.4 − SGA 62.4 = patm − SGA × 15.6
12 12
6.24
Luego patm − 6.24 = patm − SGA × 15.6 → SGA = = 0.4
15.6
43
Problema 9
Solución
44
Problema 10
Tome el peso especı́fico del agua como 62.4 lbf /f t3 . La anchura de la compuerta
es 8 f t.
Solución
L = 10 f t; b = 8 f t; A = 80 f t2
45
1. pCG = 100 × 144 − 62.4 × (35 − 20) = 13464 psf .
1
b L3
2. yCP = − sin 90◦12 γH2 O = −0.03862 f t. Posición CP
F
MA
Fext = = 1.068800 × 106 lbf
5
46
Problema 11
Tome el peso especı́fico del agua como 1000 kg/m3 . La anchura de la compuerta
es 2.5 m y su longitud 3 m.
Solución
1 1
b L3 2.5 × 33
2. yCP = −γH2 O sin 60 12 = −9806 sin 60 12 = −0.00961 m =
F 4.969 × 106
−9.61 mm
`A = (L/2) − |yCP | = 1.5 − 9.61 × 10−3 ≈ 1.49 m desde A
3. MA = F · `A sin 90 = 4.969 × 106 × 1.49 = 7.404 × 106 N · m
47
L = 3 m; b = 2.5 m; A = 7.5 m2
MA 7.404 × 106
Fext = = ≈ 4.936 × 106 N
L/2 1.5
48
Problema 12
Tome el peso especı́fico del agua como 1000 kg/m3 . La anchura de la compuerta
es 3 m.
Ignore la fricción y el peso de la compuerta.
Solución
h h
L = 9 m; b = 3 m; Lsub = ; Asub = b Lsub ; hCG sub =
sin 60 2
h h
1. F = γH2 O hCG sub Asub = 9806 3 = 16985 h2 N
2 sin 60
3
1 1 h
b L3 3
12 sin 60
2. yCP = − sin 60γH2 O 12 = − sin 60×9806 = −0.192h m
F 16985 h2
49
Lsub 1
3. `A = − |yCP | = h − 0.192 = 0.385h
2 2 sin 60
Equilibrio de torques:
h3 = 49.55 → h = 3.67 m
4. F ≈ 228769 N
5. yCP = −0.7 m
50
Problema 13
1. Halle el valor de h.
2. Si la presión absoluta en el tanque se duplica, ¿Cuál será la nueva difer-
encia de alturas entre el nivel del mercurio en el manómetro?
Solución
1. p = 150 kP a ; pB = 100 kP a
p − pB
p − pB = 50 kP a = ρHg g h → h=
ρHg g
50000 P a
h= = 0.375 m
13600 × 9.8
2. Si pA = 300 kP a
300000 − 100000
h0 = = 1.5 m
13600 × 9.8
51
Problema 14
El aire que se encuentra en la parte superior del lı́quido está sometido a una
presión de 40 kP a manométrica. La gravedad especı́fica del lı́quido del tanque
es de 0.8. A la derecha de la compuerta AB, la cual puede girar en torno a una
bisagra situada en B, se encuentra un tanque con un gas presurizado a 60 kP a
manométrica. Si la compuerta es rectangular y mide 2 m de ancho, con y1 = 3 m
y y2 = 1 m, halle las siguientes variables:
1. La fuerza de presión ejercida por el gas a la derecha sobre la compuerta.
2. La distancia del centro de presiones a la derecha respecto de la bisagra B.
3. La fuerza de presión ejercida por el lı́quido a la izquierda sobre la com-
puerta.
4. La distancia del centro de presiones a la izquierda respecto de la bisagra
B.
5. La fuerza P, ejercida en la parte inferior de la compuerta, necesaria para
mantener la compuerta en su lugar.
Solución
b = 2m L = 1m γ` = SG` γH2 O
1. Fd = pCGd · A = 60000 × 2 × 1 = 120000 N
L
2. yCP d = 0 pues la presión es uniforme `d = = 0.5 m
2
52
1
b L3 1 2 × 13
4. yCP i = −γ` sin 90 12
◦
= −0.8 × 9800 = −9.69 × 10−3 m
Fi 12 134880
P
5. MB = 0 → Fi `i − Fd `d − P L = 0
L 1
Con `i = + |yCP i | = + 9.69 × 10−3 = 0.51 m
2 2
Fi `i − Fd `d 134880 × 0.51 − 120000 × 0.5
P = =
L 1
P = 8788.8 N
53
Problema 15
Solución
pE = patm pC = 90 kP a Manométrica
54
(pC − pE ) − γH2 O (0.9)
= SGH
γH2 O (0.6)
90588
pA = 90588 P a = 90.6 kP a = P si = 13.14 psi
6894.8
1 psi = 6894.8 P a
55
Problema 16
Solución
b = 2m Asubd = d1 b Asubi = d2 b
3.2
1. Fd = pCGd Asubd = ( + γ` hCGd )d1 b = 8800 ×
patm
× 3.2 × 2 = 90112 N
2
56
1
b d31 1 2(3.2)3
◦ 12
2. yCP d = γ` sin 90 = −8800 = −0.53 m
Fd 12 90112
d1
`d = (h − d1 ) + + (yCP d ) = (4 − 3.2) + 1.6 + 0.53 = 2.93 m
2
2
patm
3. Fi = pCGi Asubi = ( +γH2 O ·hCGi )d2 b = 9800× ×(2×2) = 39200 N
2
1
b d32 1 2(2)3
4. yCP i = −γH2 O sin 90 12
◦
= −9800 = −0.33 m
Fi 12 39200
d2
`i = (h − d2 ) + + (yCP i ) = 2 + 1 + 0.33 = 3.33 m
2
P
5. M0 = 0 W L + Fi `i − Fd `d = 0 W = γW V ol
57
Problema 17
Solución
PD = Patm PB = 120 kP a
58
2. pA − pB = SGoil γH2 O · 1.3 → pA = pB + SGoil γH2 O · 1.3
pA = 120000 + 0.9 × 9800 × 1.3 = 131466 P a
131466
pA = 131.47 kP a = = 19.1 psi
6894.8
1 psi = 6894.8 P a
59
Problema 18
Solución
` = longitud de ` = 2m
la compuerta 1
L = 2 m + 2 = 2.5 m
sumergida 4
b = 1m
60
Trabajamos con presiones manométricas
MF O 13066.67
T= = = 5226.67 N
L 2.5
T + E = W → T + γH2 O V ol = γ · V ol
(γ − γH2 O ) V ol = T
T 5226.67
V ol = = = 0.38 m3
γ − γH2 O 23600 − 9800
61
Problema 19
El manómetro de la figura consiste en tres fluidos: agua a la izquierda, un aceite
(con SGoil ) en el centro y un lı́quido desconocido (con SGl ) a la derecha. Halle
el valor de la presión absoluta en el punto A
Solución
62
Problema 20
Solución
L
h0 = sin (90 − θ)
2
L
h0 = cos θ
2
F = pCG A
L L
pCG = pA + γ` h1 − γ` h2 − γ` cos θ = pA + γ` h1 − h2 − cos θ
2 2
L
F = pA + γ` h1 − h2 − cos θ b L
2
63
1
b L3 bL3
yCP = −γ` sin (90 − θ) 12 = −γ` cos θ
F 12F
bL3
L L
M =F − |yCP | sin 90 = F − γ` cos θ
2 2 12
64
Problema 21
Solución:
65
Problema 22
Nota: Suponga que se impide el movimiento del lı́quido B hacia la parte supe-
rior del tanque
Solución
66
• pA − pC = −γA H
pA = pC − γA H
• pB − pA = (pB − pD ) + (pD − pA ) = γB h + γA (H − h)
pB − pC +
γA
H H − γA h
= γB h + γA
pB − pC = h(γB − γA ) →
pC − pB
→ h=
γA − γB
67
Problema 23
Solución
h
LS =
sin θ
bh
A = b · LS =
sin θ
Fuerza de presión
h bh bh2
Fp = pCG A = γ · =γ
2 sin θ 2 sin θ
Punto de aplicación
1 3 1 h3
bLS b
12 sin3 θ
yCp = −γ sin θ 12 = −γ
sin
θ 2
Fp b
h
γ
2sin θ
1 h
yCp = −
6 sin θ
LS h 1 h h 1 1 h
`= − |yCp | = − = 1− =
2 2 sin θ 6 sin θ 2 sin θ 3 3 sin θ
P
Condición de equilibrio MA = 0
Fp · ` 1 bh2 1 h
Fp · ` = F L → F = = ·γ ·
L L 2 sin θ 3 sin θ
68
1 bh3
F = γ 2
6L sin θ
69
Problema 24
Solución
70
P
Condición de equilibrio F =0 F + W = Fe
F + γS a3 = γa2 d →
F + γS a3
→ d=
γa2
Condición de estabilidad
1 3
I0 12 a a 1 1 2 1
MG = − GB = − (a − d) = a d − (a − d)
VSub 2
a d 2 12 2
Si MG > 0 Estable
MG < 0 Inestable
Para mejorar la estabilidad se podrı́a aumentar el área de la lı́nea de flotación
conservando el volumen
71
Problema 25
Solución
72
L
pCG,oil = pentre f ase − γoil = 13401.6 − 0.8 · 62.4 · 3 = 13251.84 psf
4
Fw = pCG,w × Aw = 13588.8 · 36 = 489197 lbf
FT · `T = Fw · `w + Foil · `oil ⇒
Fw · `w + Foil · `oil 489197 · 9.014 + 477066 · 3.011
`T = =
FT 966263
`T = 6.05 f t
73
Problema 26
Solución
74
Problema 27
Solución
75
Puntos de aplicación respecto a B.
Aceite: 3
1 L 1
b · 1 · 23
12 2 12 1
yCP,oil = −γoil sin 90 = −0.8 · 9800 =− m
Foil 15680 3
L L 4 1 1 8
`oil = + − |yCP,oil | = + 1 − =3− = m
2 4 2 3 3 3
Agua:
1
· 1 · 23
yCP,w = −9800 · 12 = −0.128 m
50960
L
`w = − |yCP,w | = 1 − 0.128 = 0.872 m
4
Fuerza horizontal que actua en CG de AB=L0
L L L
F = γoil + γw (b L0 ) = γw (1 + 0.8)b L0
2 2 2
F = 9800 · 2 · 1.8 · 1 · L0 = 35280 L0 [N ]
L0
`= → Brazo de F respecto a B
2
P
Equilibrio ⇒ TB = 0
L0
Fw · `w + Foil · `oil = F ` = 35280 L0
2
Fw · `w + Foil · `oil 50960 · 0.872 + 15680 · 8/3
L02 = =
17640 17640
L02 = 4.89 m2 → L0 = 2.21 m
76
Problema 28
Solución
ρs gV
s = ρw gV
w + ρl gV
l
ρ
0.8 w ρ
b (3L)(6L) = w ρ
b L (6L) + 1.2 wb ẏ(6L)
1.4 L
2.4 L = L + 1.2 ẏ → 1.4 L = 1.2 ẏ → ẏ = = 1.17L
1.2
Para que el bloque esté a ras del agua, se debe hallar SGs , y en este caso ẏ = 2 L
ρs Vs = ρw Vw + ρl Vl ρs = SGs ρw
SGs
ρ
wb (3L)
(6L)
= ρ
w b L
(6L)
+ 1.2 ρ
w b (2L)
(6L)
3.4
SGs = = 1.13
3
77
Problema 29
La compuerta AB, la cual puede girar en torno a una bisagra situada en B, cierra
un depósito el cual se encuentra abierto a la atmósfera. El depósito contiene
agua en su parte inferior y aceite en su parte superior como muestra la figura.
Dicha compuerta posee una longitud L, una anchura b y un espesor t. En estas
condiciones determine la densidad del material de la compuerta para que ésta
comience a abrirse.
Datos: L = 4 m, b = 1 m, t = 0.1 m, SGoil = 0.8, γagua = 9800 N/m3 , θ = 45◦ .
Solución
P
La compuerta comenzará a abrirse cuando TB = 0
L L
Fpo `o + Fpw `w = W cos θ `W =
2 2
`o = CP0 B ; `w = CPw B
Aceite:
L L
hCG,o = sin θ pCG,o = γo sin θ = 5543.7 P a
4 4
L L 4 4
Fpo = γo sin θ b = 0.8 · 9800 · sin 45◦ · 1 ·
4 2 4 2
78
L
Fpo = 11087.4 N `o = + |yCP,o |
4
3
1 L 1
b √ 1 · 23
12 2 2 12
yCP,o = −γ0 sin 45 = −0.8 · 9800 · ≈ −0.33 m
Fpo 2 11087.4
4
`o = + 0.33 = 1.33 m
4
Agua:
L L
pCG,w = pCG,o + γo sin θ + γw sin θ
4 4
√
4 2
pCG,w = 5543.7 + 5543.7 + 9800 = 18017 P a
4 2
L 4
Fpw = pCG,w b = 18017 · 1 · = 36034 N
2 2
L L 3L
`w = + + |yCP,w | = + |yCP,w |
2 4 4
3
1 L 1
b √ 1·8
12 2 2 12
yCP,w = −γw sin 45 = −9800 = −0.128 m
Fpw 2 36034
4
`w = 3 + 0.128 = 3.128 m
4
P
Equilibrio TB = 0
79
Problema 30
Solución
Resolviendo el manómetro
1 16
pair = −γ` 16 = −62.4 · = 83.2 lbf /f t2
12 12
80
1 3
Io L L 1
MG = − GB = 12 − (L − H)
Vsub 3 2
0.7
L
L 1
MG = − (L − 0.7 L) = 0.119 L − 0.15 L < 0
12 · 0.7 2
Del resultado anterior se puede concluir que la flotación es inestable
81
Problema 31
La compuerta angular ABC, la cual puede girar en torno a una bisagra situada
en B, cierra un depósito el cual se encuentra abierto a la atmósfera. Los lados
AB y BC de la compuerta tienen longitud L. El depósito se encuentra lleno de
agua como muestra la figura. En estas condiciones encuentre la fuerza P que
debe aplicarse en el punto A para que la compuerta comience a abrirse.
Solución
P
La compuerta comenzará a abrirse cuando TB = 0
L L
F2 = P L + F1 `1 pues `2 =
2 2
Lado inclinado:
L
pCG,1 = γw hCG,1 = γw cos θ
2
L L
F1 = pCG,1 b L = γw cos θ b L ; `1 = + |yCP,1 |
2 2
82
1 3
b L 1
yCP,1 = − ◦((
( − θ) 12
(( = − 6 L
(
γsin(90
w( (
γ
(w(L/2
( (cos
( θ b L
◦
pues sin(90 − θ) = cos θ
L L L
`1 = − =
2 6 3
Lado horizontal:
pCG,2 = γw L cos θ
L
F2 = γw L cos θ b L `2 = (es horizontal)
2
P
Equilibrio TB = 0
L
F2 = F1 `1 + P L
2
L L2 L
γw b L2 cos θ = γw b cos θ + P
L
2 2 3
1 1 1
P = γw b L2 cos θ − = γw b L2 cos θ
2 6 3
1
P = · 62.4 · 1 · 22 · cos 30◦
3
P = 72.05 lbf
83
Problema 32
Una esfera sólida de aluminio (SG = 2.8) y otra de cobre (SG = 8.9) pesan lo
mismo cuando están sumergidas en un lı́quido. Si el diámetro de la esfera de
aluminio es de 3 in ¿cuál debe ser el diámetro de la esfera de cobre si el lı́quido
es agua? Si la esfera de cobre tiene un diámetro de 1.5 in, ¿cuál será la densidad
del lı́quido?
Solución
π 3 π 3 π 3 π 3
SGAl γw D − SG` γw DAl = SGCu γw DCu − SG` γw DCu
6 Al 6 6 6
1/3
3 3 SGAl − SG`
DAl (SGAl − SG` ) = DCu (SGCu − SG` ) → DCu = DAl
SGCu − SG`
3 3 3 3
SGAl DAl − SGCu DCu = SG` DAl − DCu
3 3
SGAl DAl − SGCu DCu 2.8 · 33 − 8.9 · 1.53
SG` = 3 − D3 =
DAl Cu 33 − 1.53
75.6 − 30.0375
SG` = = 1.93
23.625
ρ` = 1.93ρw = 1.93 · 1.94 = 3.74 slug/f t3
84
Problema 33
Solución
P
Para que la compuerta comience a abrirse TB = 0
85
Gas:
L 4
Fpg = pgas b = 30 · 144 · · 1 = 8640 lbf
2 2
L
`g = = 1 f t
4
Lı́quido:
L L 4 4
Fpl = γl + pgas b = 2 · 62.4 · + 30 · 144 ·1
4 2 4 2
L L 3L
Fpl = 8889.6 lbf `l = + + |yCP,l | = + |yCP,l |
2 4 4
3
1 L 1
b 1 · 23
◦ 12 2 12
yCP,l = −γl sin 90 = −2.0 · 62.4 · ≈ −0.010 f t
Fpl 8889.6
4
`l = 3 + 0.010 = 3.010 f t
4
Agua:
L L
Fpw = γw hCG,w L b = γw h− L b; `w = + |yCP,w |
2 2
1
b L3 1 3
L2
12 b L 1
yCP,w = −γw sin 90 = −
γ
w =−
Fpw 12 w (h − L/2)
γ Lb 12 (h − L/2)
L 1 L2
`w = +
2 12 (h − L/2)
Equilibrio Fpg · `g + Fpl · `l = Fpw `w
L2
L L 1
8640 · 1 + 8889.6 · 3.01 = γw h − Lb +
2 2 12 (h − L/2)
L2
1 L L 1
L
35397.7 · = h− + h−
−
L γw b 2 2 2 12 (h L/2)
35397.7 · 1 4 4 1 4
= h− + · 16 = (h − 2)2 +
4 · 62.4 · 1 2 2 12 3
4 1
141.8 − + 2 = h → h = 72.2 f t
3 2
86
Problema 34
En un estuario donde el agua dulce se mezcla con agua de mar ocurre una
condición de salinidad estratificada donde el agua dulce está en la parte supe-
rior y la salada en la inferior. La entrefase entre ambas se llama una haloclina.
Una idealización de ella es la presentada en la figura donde en cada lado de ella
hay una densidad constante. Una esfera de diámetro 35 cm y que pesa 50 lbf
flotarı́a alrededor de dicha haloclina. Halle el porcentaje del volumen de la es-
fera sumergido en agua salada.
Solución
Hallamos γs de la esfera
W 50
W = γs V → γs = = π
V D3
6
1 in 1 f t
D = 35 cm · = 1.148 f t
2.54 cm 12 in
50
γs = π = 63.12 lbf /f t3
1.1483
6
Posición de flotación W = Fe,dulce + Fe,salada
γs V = γdulce (1 − d) V + γsalada d V
87
γs = 1.01
γ
w = 1.0 w (1 − d) + 1.025
γ γ
wd
0.01
1.01 − 1.0 = d(1.025 − 1.0) → d = = 0.4
0.025
El 40% del volumen de la esfera está sumergido en el agua salada
88
Problema 35
Figura Problema 35
Solución
3. FB = W = m g
FB = 790 · 9.8 · 0.1 · 0.1 · 0.085 + 1000 · 9.8 · 0.1 · 0.1 · 0.015 = 8.05 N
89
FB 8.05
m= = = 0.821 kg
g 9.8
m 0.821
m = ρV → ρ= = = 821 kg/m3
V 0.13
90
Problema 36
Solución
91
Problema 37
Solución
h h
Beta = 90 − α, sin α = , Ls =
Ls sinα
Resultante de fuerzas de presión sobre compuerta es el empuje Fe
Compuerta en equilibrio → La suma de torques sobre B es cero.
92
P
TB = 0
FL ( (( + F Ls (
sin(Beta (( = W L (
sin(Beta sin(Beta
((
e
2 2
h
Con Fe = γbt Ls , W = γs btL Fe = γbt
sinα
h 1 h L2
Entonces F L + γbt × × = γs bt
sinα 2 sinα 2
2 2
1 h L
γbt = γs bt − FL
2 sin2 α 2
L2 1
sin2 α (γs bt − F L) = γbt h2
2 2
v v
u 1 2
u 1
u γbth u 64 × × 4 × 0.5 × 36
2 2
u u
sin α = u =u
t L2 t 100
γs bt − FL 500 × × 4 × 0.5 − 4500 × 10
2 2
r
2304
sin α = = 0.6788 → α = 42.75◦
50000 − 45000
93
Problemas Propuestos
Problema 1
Problema 2
Un bloque cúbico de madera puede girar sobre uno de sus lados tal y como
muestra la figura. El bloque se encuentra en equilibrio cuando se sumerge en
agua a la profundidad que se muestra. Despreciando la fricción en el pivote,
encuentre la gravedad especı́fica de la madera.
94
Problema 3
Problema 4
95
Problema 5
Datos: Tome el peso especı́fico del agua como 62.4 lbf /f t3 , SGHg = 13.6, θ
= 30◦ , H = 20 ft, h1 = 1 ft, h2 = 2 ft, h3 = 10 ft, la longitud de la compuerta
es 4 ft y su anchura 2 ft.
Explique razonadamente cada una de las etapas que usted aplicó para la re-
solución del problema.
Problema 6
96
Problema 7
Problema 8
97
Problema 9
Una estructura está sujeta al fondo del océano como se muestra en la figura. En
una pared inclinada hay una escotilla de diámetro D articulada en uno de los
bordes. Determinar la presión mı́nima del aire, p1 , que debe haber dentro del
contenedor para abrir la escotilla. Ignorar su peso y la fricción en la articulación.
πR4
Datos: ρ = 1025 kg/m3 , h = 10 m, D = 2 m, Ixx = ( R = D/2), θ =
4
30◦ .
Problema 10
98
Capı́tulo 3
Cinemática
El campo de velocidades
a) Descripción Lagrangiana:
99
b) Descripción Euleriana:
i) Lı́nea fluida: Lugar geométrico definido por las posiciones de una familia de
partı́culas que ocupan inicialmente una lı́nea continua
x` = x x0 (α), t x0 ∈ LF (t0 )
Una lı́nea fluida cerrada permanecerá cerrada siempre.
ii) Superficie fluida: Lugar geométrico definido por las posiciones de las partı́culas
fluidas que inicialmente ocupan una superficie
100
xs = x x0 (α, β), t
iii) Volumen fluido: Lugar geométrico definido por la posición de un conjunto
de partı́culas fluidas que inicialmente ocupan un volumen
xV = x x0 (α, β, γ), t
b) Movimiento virtual
c) Trayectoria (pathline)
101
d~x = (dx1 , dx2 , dx3 ) ~υ = (υ1 , υ2 , υ3 )
dx1 dx2 dx3
= = ⇐ Ecuación que define la lı̀nea de corriente
υ1 υ2 υ3
iii) Tubo de corriente: Volumen encerrado por una superficie de corriente tal
que la lı́nea cerrada sobre la que se apoya es cerrada
~υ · ~n = 0
102
e) Traza (streakline)
La traza es la lı́nea que se forma cuando se inyecta tinta en un punto del flujo.
La figura muestra un flujo cuya dirección varı́a con t. En este caso las lı́neas de
corriente son las rectas de color verde, la trayectorı́a es el camino rojo recorrido
por la cometa y la traza es la curva gris que representa el humo emitido por la
chimenea. Es evidente que en este caso las tres lı́neas son diferentes.
103
Flujo estacionario ⇒ Traza = trayectoria = lı́nea de corriente que pasan por un
punto del espacio
Ejemplo 3.1:
Solución:
g) Vorticidad y circulación
~ (~x, t) = ∇ × ~υ (~x, t)
ω
104
3.2 Derivada material. Aceleración de una partı́cula fluida
∂φ
En la descripción Euleriana del campo φ(~x, t), indica la evolución
∂t x=cte
temporal de las propiedades de diversas partı́culas que pasan por un punto en
instantes consecutivos, es decir, no suministra información sobre la evolución
de la partı́cula fluida que estaba en ~x en el instante t.
En el caso de la velocidad
∂υ υ(~x, t + ∆t) − υ(~x, t)
(~x, t) = lim
∂t ∆t→0 ∆t
no coincide con la definición de aceleración de una partı́cula fluida dada en el
contexto lagrangiano:
∂~x(~x0 , t)
~υ = ~x0 f ijo
∂t
∂~υ (~x0 , t)
~a = ~x0 f ijo
∂t
En el cálculo de la variación de φ siguiendo
la evolución de la partı́cula fluida
~x = ~x(~x0 , t) ⇒ φ(~x, t) = φ ~x(~x0 , t), t
Solución:
Du ∂u ∂u ∂u ∂u
ax = = +u +υ +w = 3 + 3t · 0 + xz · 0 + ty 2 · 0 = 3
Dt ∂t ∂x ∂y ∂z
105
Dυ ∂υ ∂υ ∂υ ∂υ
ay = = +u +υ +w = 0 + 3t · z + xz · 0 + ty 2 · x = 3tz + xty 2
Dt ∂t ∂x ∂y ∂z
Dw ∂w ∂w ∂w ∂w
az = = +u +υ +w = y 2 +3t·0+xz ·2ty +ty 2 ·0 = y 2 +2xyzt
Dt ∂t ∂x ∂y ∂z
Finalmente
D~υ
~a = = (3, 3tz + xty 2 , y 2 + 2xyzt)
Dt
106
Problemas Resueltos
Problema 1
Solución
y
u = u0 sin ω t − ; v = v0 ; w = 0
v0
~
DV ~
∂V
1. = ~ · ∇)V
+ (V ~
Dt ∂t
∂u ∂u ∂u ∂u
ax = +u +v +
w =
∂t ∂x ∂y ∂z
y u0 ω y
= u0 ω cos ω t − +0 −
v cos ω t − =0
v0 v
0 v0
∂v ∂v ∂v ∂v
ay = + u + v + w = 0
∂t ∂x ∂y ∂z
∂w + u ∂w
+ v ∂w ∂w
az = +
w =0
∂t ∂x ∂y ∂z
107
î ĵ k̂
~ = ∂
∂ ∂ ∂u
~ =∇×V
2. ω = î 0 + ĵ 0 + k̂
∂x
∂y ∂z ∂y
u v0 0
u0 ω y
ω
~ = 0, 0, cos ω t −
v0 v0
~
5. V ~ =0
incompresible ⇔ ∇ · V
~ = ∂u
∇·V + ∂v
∂w
+
= 0 → ~
V Incompresible
∂x ∂y ∂z
108
Problema 2
Solución
∂vi X ∂vi
u = yzt, v = xzt, w = xyt ai = + vj
∂t j
∂xj
~a = yz + xz 2 t2 + xy 2 t2 , xz + yz 2 t2 + yx2 t2 , xy + zy 2 t2 + zx2 t2
î ĵ k̂
= î ∂w − ∂v − ĵ ∂w − ∂u +
∂ ∂ ∂
~ = ∇ × ~v =
2. ω
∂x ∂y ∂z
∂y ∂z ∂x ∂z
u v w
∂v ∂u
k̂ −
∂x ∂y
109
Problema 3
Solución
~v = e−t (1, z, y)
∂u ∂u ∂u ∂u
1. ax = +u +v +w = −e−t + 0 = −e−t
∂t ∂x ∂y ∂z
∂v ∂v ∂v ∂v
ay = +u +v +w = −ze−t + e−t (1 · 0 + z · 0 + y · e−t )
∂t ∂x ∂y ∂z
= e−t (ye−t − z)
∂w ∂w ∂w ∂w
az = +u +v +w = −ye−t + e−t (1 · 0 + z · e−t + y · 0)
∂t ∂x ∂y ∂z
= e−t (ze−t − y)
110
Problema 4
Solución
~ = u0 y 2 z î + v0 z 2 ĵ + e−kt2 k̂
V
~
DV ~
∂V
1. = ~ · ∇)V
+ (V ~
Dt ∂t
∂u ∂u ∂u ∂u
ax = +u +v +w =
∂t ∂x ∂y ∂z
2
= 0 + 0 + v0 z 2 2u0 yz + e−kt u0 y 2
∂v ∂v ∂v ∂v
ay = +u +v +w =
∂t ∂x ∂y ∂z
2
= 0 + 0 + 0 + e−kt 2zv0
∂w ∂w ∂w ∂w
az = +u +v +w =
∂t ∂x ∂y ∂z
2
= −2kte−kt + 0 + 0 + 0
î ĵ k̂
∂ ∂ ∂
~
~ = ∇ × V =
2. ω = −î 2v0 z + ĵ u0 y 2 − k̂ 2u0 yz
∂x ∂y ∂z
u v w
−2v0 z , u0 y 2 , −2u0 yz
ω
~ =
111
~
5. V ~ =0
incompresible ⇔ ∇ · V
~ = ∂u
∇·V + ∂v
∂w
+
= 0 → ~
V Incompresible
∂x ∂y ∂z
112
Problema 5
Solución
~
DV ~
∂V
1. = ~ · ∇)V
+ (V ~
Dt ∂t
∂u ∂u ∂u ∂u
ax = +u +v +w =
∂t ∂x ∂y ∂z
= u0 ω cos ωt + 0 + 0 + 0
∂v ∂v ∂v ∂v
ay = +u +v +w =
∂t ∂x ∂y ∂z
= −v0 ω sin ωt + 0 + 0 + 0
∂w ∂w ∂w ∂w
az = +u +v +w =
∂t ∂x ∂y ∂z
−k(x2 +y 2 ) 2 2
= 0 + u0 sin ωt(−2xk)e + v0 cos ωt(−2yk)e−k(x +y ) + 0
î ĵ k̂
~ = ∂ ∂ ∂ = î (−2yk)e−k(x2 +y2 ) − ĵ (−2xk)e−k(x2 +y2 ) +
~ = ∇×V
2. ω ∂x
∂y ∂z
u v w
k̂ 0
2
+y 2 ) 2
+y 2 )
ω
~ = (−2yk)e−k(x , (2xk)e−k(x , 0
113
4. No es homogéneo porque depende del espacio
~
5. V ~ =0
incompresible ⇔ ∇ · V
~ = ∂u
∇·V + ∂v
∂w
+
= 0 → ~
V Incompresible
∂x ∂y ∂z
114
Problema 6
~ = xe−kt2 î + xz 2 ĵ + ln yt k̂
V
Solución
~ = xe−kt2 î + xz 2 ĵ + ln yt k̂
V
~
DV ~
∂V
1. = ~ · ∇)V
+ (V ~
Dt ∂t
∂u ∂u ∂u ∂u
ax = +u +v +w =
∂t ∂x ∂y ∂z
2 2 2
= −2ktxe−kt + xe−kt e−kt + 0 + 0
∂v ∂v ∂v ∂v
ay = +u +v +w =
∂t ∂x ∂y ∂z
2
= 0 + xz 2 e−kt + 0 + 2xz ln yt
∂w ∂w ∂w ∂w
az = +u +v +w =
∂t ∂x ∂y ∂z
y t 1 xz 2
= + 0 + xz 2 + 0 = +
yt
y t t y
î ĵ k̂
~ = ∂
∂ ∂ t
~ =∇×V
2. ω = − î − 2xz + ĵ 0 − k̂ z 2 − 0
∂x
∂y ∂z y t
u v w
1 2
ω
~ = − 2xz , 0 , z
y
115
4. No es homogéneo porque depende del espacio
~
5. V ~ 6= 0
No es incompresible ⇔ ∇ · V
~ = ∂u + ∂v
∇·V
∂w
= e−kt2 + 0 + 0 6= 0
+
∂x ∂y ∂z
116
Problemas Propuestos
Problema 1
~ = −y 2 x2
V î + ĵ + 0 k̂
x2 + y 2 x2 + y 2
Responda las siguientes preguntas:
2. Halle su vorticidad.
3. Dicho campo de velocidades, ¿es estacionario? ¿por qué?
4. Dicho campo de velocidades, ¿es homogéneo? ¿por qué?
Problema 2
2. Halle su vorticidad.
3. Dicho campo de velocidades, ¿es estacionario? ¿por qué?
4. Dicho campo de velocidades, ¿es homogéneo? ¿por qué?
117
Problema 3
~ = y p
V î + ln x2 + y 2 ĵ + 0 k̂
x2 +y 2
Problema 4
~ = y − x î +
V
x
ĵ + 0 k̂
2
x +y 2 x + y2
2
2. Halle su vorticidad.
3. Dicho campo de velocidades, ¿es estacionario? ¿por qué?
4. Dicho campo de velocidades, ¿es homogéneo? ¿por qué?
118
Capı́tulo 4
Dinámica
Las leyes de la mecánica describen que sucede cuando hay una interacción entre
el sistema y su entorno.
Leyes de la mecánica
119
4. Si se añade calor al sistema dQ o el sistema realiza trabajo dW , su energı́a
dE cambia de acuerdo a la 1a Ley de la termodinámica
dE dQ dW
= −
dt dt dt
dQ
5. Finalmente, la Segunda Ley de la termodinámica dS ≥
T
Todas estas leyes se complementan con las relaciones de estado
p = p(ρ, T ); e = e(ρ, T )
Las leyes básicas se reformulan para ser aplicadas a esta región local llamada
V.C. Todo lo que se necesita saber es el flujo en esta región y frecuentemente
simples asunciones serán suficientemente aproximadas siendo las condiciones
fuera del VC irrelevantes.
Si VC(t) la cuestión no es tan fácil pues las fronteras del VC dependen del
tiempo. El caso general se trata en el Teorema de Transporte de Reynolds
(TTR).
120
1a forma Tma transporte de Reynolds
d R R ∂Ψ
ΨdV = + ∇ · (Ψ~υ ) dV
dt V C V C(t) ∂t
d R R ∂Ψ H
ΨdV = dV + Ψ~υ · ~n da
dt V C V C(t) ∂t SC(t)
d R d R H
Ψ(~x, t) dV = Ψ(~x, t)dV + Ψ(~υ − ~υc ) · ~n da
dt V F (t) dt V C(t) SC(t)
121
donde:
~υ : Velocidad del fluido en SC
~υc : Velocidad del volumen de control en el movimiento virtual
Ejemplo 4.1:
Solución:
122
como en el globo ~υ −~υc = 0, pues el V C se expande con el globo, y en la entrada
~1
~υ − ~υc = V Z Z
dm d ~1 · ~n da
= ρ dV + ρ1 V
dt sist. dt
V C(t) A1
El flujo másico definido como ṁ = ρV1 A1 [kg/s] sólo tiene lugar en la entrada
por lo que podemos aproximar la masa en el globo por ρb Vb (t) con Vb (t) el
volumen del globo en el instante t. Se debe tener en cuenta que en la entrada
~1 y ~n1 son antiparalelos, por lo que V
V ~1 · ~n1 = −V1
dm d 4 3
= ρb πR (t) − ρ1 A1 V1
dt sist. dt 3
dm
Pero por conservación de la masa = 0 luego
dt sist
d 3
(ρb R) = ρ1 A1 V1
dt 4π
4.2 Conservación de la masa. Ecuación de continuidad
R dF (t)
Si F (t) = Ψ(~x, t)dV es una cantidad conservada entonces = 0 y
V F (t) dt
según el TTR-1:
Z
∂Ψ ∂Ψ
∀V C + ∇ · (Ψ~υ ) dV = 0 ⇔ + ∇ · (Ψ~υ ) = 0
∂t ∂t
V C(t)
R
Si Ψ = ρ, m(t) = ρ(~x, t)dV . Por hipótesis la masa del sistema se conserva
V F (t)
(primera ley de la Mecánica) luego el campo de densidades verificará
123
∂ρ
+ ∇ · (ρ~υ ) = 0
∂t
Ejemplo 4.2:
Solución:
d
R
Flujo estacionario ⇒ dt ρ dV = 0
V C(t)
I
ρ(~υ − ~υc )~n da = 0
SC
124
Aquı́ el V C está fijo y por definición de tubo de corriente no hay flujo lateral,
luego
ρ2 V2 A2 − ρ1 V1 A1 = 0 ⇒ ρ2 V2 A2 = ρ1 V1 A1 = ṁ = cte
i.e. el flujo másico es constante en cada sección. Si ρ = cte
Q = V1 A1 = V2 A2
y el caudal es constante en cada sección
Generalización del resultado anterior
utilizando el TTR-1
Z Z
[(ρφ),t +(ρφυi ),i ]dV = ρ[φ,t + φ,i υi ] + φ[ρ,t + (ρυi ),i ] dV
VF VF
Dφ
ρ
Dt
125
Ejemplo 4.3:
El tanque de la figura se llena con agua por dos entradas 1D y la parte superior
se encuentra llena de aire. La altura del agua es h.
dh
1. Hallar una expresión para el cambio en la altura del agua
dt
dh
2. Calcular si D1 = 1 in, D2 = 3 in, V1 = 3 f t/s, V2 = 2 f t/s y At =
dt
2 f t . Suponer H2 O a 20◦
2
Solución:
126
dh
(b) Q1 = 0.016 f t3 · s−1 ; Q2 = 0.093 f t3 · s−1 ⇒ = 0.057 f t · s−1
dt
Reto: Intente repetir el ejercicio pero tomando un VC que contenga solamente
agua y cuya superficie superior se mueva conforme el nivel de agua en el depósito
aumenta.
Por tanto, todo el volumen fluido que entra al volumen de control debe salir o
incrementar el volumen VC.
Nota: Si hay dos o más fluidos incompresibles con diferentes densidades, esta
ecuación es diferente de la ecuación de conservación de la masa.
Ejemplo 4.4:
Solución:
Con las asunciones habituales de flujo 1D, como los fluidos son incompresibles:
−ρ1 Q1 − ρ2 Q2 + ρm Qm = 0
127
Con las dos incógnitas de la densidad y el caudal de la mezcla. La ecuación
adicional nos la proporciona la conservación del volumen:
Z Z
d
dV + ~v − ~vc · ~n dS = 0
dt V C(t) S1 +S2 +Sa +
S
p
Es decir, −Q1 − Q2 + Qm = 0
128
que constituye la 2a Ley de Newton para el medio continuo.
T~ = τ · ~n
donde ~n es la normal exterior a SF y τ es el tensor de esfuerzos.
D~υ
ρ = ∇ · τ + f~
Dt
Ejemplo 4.5:
129
Solución:
π
= (p − po )A1 = (p − po ) D12 = 177 lbf.
4
Ejemplo 4.6:
Solución:
130
Por ser tubo de corriente las únicas contribuciones no nulas estén en la entrada
y salida X
F~ = −ρ1 A1 V1 ~υ1 + ρ2 A2 V2 ~υ2 = −ṁ1 ~υ1 + ṁ2 ~υ2
Pero por la ecuación de continuidad ṁ1 = ṁ2 = ṁ
X
F~ = ṁ(V
~2 − V
~1 )
Ejemplo 4.7:
Solución:
(a) Suponemos también entrada y salida 1D, luego estamos en un tubo de cor-
F~ = ṁ2 V
~2 − ṁ1 V
~1 ṁ1 = ṁ2 = ṁ = ρAV V1 = V2 = V
riente por lo que
Fx = ṁ(V cos θ − V ) Fy = ṁV sin θ
r !
q θ θ 1 − cos θ
(b) F = Fx2 + Fy2 = 2ṁ V sen sin =
2 2 2
Fy
φ = π − arctan
Fx
131
Ejemplo 4.8:
Solución:
Nos montamos en el SR del carrito. Se deben satisfacer la conservación de la
masa y el momento lineal.
Si las entradas y salidas las consideramos 1D ṁout = ṁin , i.e.,
ρ1 A1 V1 + ρ2 A2 V2 = ρj Aj (Vj − Vc )
1
ρ1 = ρ2 = ρj ; A1 = A2 = A ; ⇒ V1 + V2 = 2(Vj − Vc )
2
Si además resulta que la deflexión es simétrica y se desprecia el peso del agua
⇒ V1 = V2 i.e., V1 = V2 = (Vj − Vc ) = 5 m/s
P
Fy = Ry = ṁ1 υ1 + ṁ2 υ2 − ṁj υj = 0
132
!
P~ ~
dV
a
Entonces la 2 Ley de Newton queda F = m~ai = m + ~arel
dt
X ~
dV
F~ − m ~arel = m
dt
es decir, aparece una fuerza adicional −m~arel debido a la aceleración relativa
entre los SR’s. La ecuación integral de cantidad de movimiento se expresa
Z Z I
X
~ d
F− ~arel dm = ρ ~υ dV + ρ ~υ (~υ − ~υc ) · ~n da
dt
CV (t) CV (t) SC(t)
Ejemplo 4.9:
Solución:
133
El flujo en el V C es no estacionario pues Vj − Vc varı́a con el tiempo. La
conservación de la masa se escribe
Z
d
ρ dV + ρ1 A1 V1 + ρ2 A2 V2 − ρj Aj (Vj − Vc ) = 0
dt
CV
Para evaluar las integrales, se supone que la masa y el momento del fluido son
pequeñas comparadas con las del carro en el V C, luego
Z Z
d dVc
ρ u dV ≈ 0 ; ax rel dm ≈ mc
dt dt
CV CV
Por tanto
dVc dVc ρ Aj
0 − mc = 0 − ρ Aj (Vj − Vc )2 ⇒ = (Vj − Vc )2 (1)
dt dt mc
Integrando (1) se puede hallar el tiempo necesario para que el carro alcance el
95% de la velocidad del chorro
ZVc Zt
−2 ρAj ρAj 1 1
(Vj − Vc ) dVc = dt = t= −
mc mc Vj − Vc Vj
0 0
ρ Aj
Vj t
Vc mc Vc 19mc
→ = si = 0.95 ⇒ t = = 9.5 s
Vj ρ Aj Vj ρ Aj Vj
1 + Vj t
mc
134
4.4 Conservación de la energı́a
υ2
Z Z I I Z
d
ρ +e dV = ρ f~M · ~υ dV + (τ · ~n)· ~υ da+ Qc da+ Qrad dV
dt 2
VF VF SF SF VF
υ2
2 D +e
υ2
Z Z
d υ 2
ρ +e dV = ρ dV ; E = +e
dt 2 Dt 2
VF VF
I Z
(τ · ~n) · ~υ da = ∇ · τ · ~υ dV ;
SF VF
I I Z
Qc da = − ~q · ~n da = − ∇ · ~q dV
SF SF VF
υ2
D
= ∇ · τ · ~υ + ρ f~M · ~υ − ∇ · ~q + Qrad
ρ +e
Dt 2
D υ2
h i
∗ ∗
ρ = ∇ · (−pI + τ ) · ~υ + p ∇ · ~υ − T r(τ · ∇~υ ) + ρ f~ · ~υ
Dt 2
∗
T r(τ · ∇~υ ) = Φ > 0 se denomina función de disipación de Rayleigh (siempre
sustrae energı́a del flujo).
135
Ecuación diferencial de la energı́a interna
d R H dQ dW
ρ E dV + ρ E(~υ − ~υc ) · ~n da = −
dt V C(t) SC(t) dt dt
donde Q > 0 implica que se añade calor al sistema y W > 0 denota trabajo
realizado por el sistema.
υ2
Aquı́ E = e + + g z y el término de trabajo se descompone en tres partes
2
dW
≡ Ẇ = ẆS + Ẇp + Ẇν
dt
p υ2
Z I
d
ρEdV + ρ e+ + + g z (~υ − ~υc ) · ~n da = Q̇ − ẆS − Ẇν
dt ρ 2
V C(t) SC(t)
υ2
I
ρ ĥ + + g z (~υ − ~υc ) · ~n da = Q̇ − ẆS
2
SC(t)
136
En el caso en que el V C tenga una serie de entradas y salidas 1D el término de
ingestión se evalua como:
υ2
I
ρ ĥ + + g z (~υ − ~υc ) · ~n da =
2
SC(t)
X υ2 υ2
X
ĥ + +gz ṁout − ĥ + +gz ṁin
out
2 out in
2 in
υ2
ĥ, y g z son promedios sobre cada sección transversal. En el caso de tener
2
un gas perfecto la entalpı́a es ĥ = cp T .
Ejemplo 4.10:
Una máquina que trabaja con flujo estacionario toma aire en la sección 1 y lo
descarga en las secciones 2 y 3. Las propiedades en las secciones son las siguien-
tes:
Solución:
137
Para hallar p3 utilizamos la conservación de la masa
donde, w representa el trabajo por unidad de masa realizado por las máquinas.
138
p
Sustituyendo ĥ = e + y dividiendo por g,
e
p1 e1 V2 p2 e2 V2
+ + 1 + z1 = + + 2 + z2 − hq + hw
ρg g 2g ρg g 2g
Con w = hw g y q = hq g son el trabajo y el calor por unidad de masa. Reorde-
nando:
p1 V2 p2 V2 e2 − e1
+ 1 + z1 = + 2 + z2 + −hq + hw
ρg 2g ρg 2g g
| {z } | {z }
altura disponible altura de pérdidas por fricción hf
p1 V2 p2 V2
+ 1 + z1 = + 2 + z2 + hw + hf
ρg 2g ρg 2g
Finalmente, Ẇ = ρgQhw .
Ejemplo 4.11:
Solución:
139
Q2 Q2
V2 = = = 61.1 f t/s
A2 (π/4)D22
Aplicando la ecuación de la energı́a para flujo incompresible
V12 V22
p1 p2
+ + z1 = + + z2 + hw + hf
ρw g 2g ρw g 2g
p1 V2 p2 V2
+ 1 + z1 = + 2 + z2 = cte Ecuación de Bernoulli
ρg 2g ρg 2g
que es una particularización de la ecuación de la energı́a en flujo estacionario y
sin fricción.
140
5. Ausencia de máquinas (i.e., sin trabajo de impulsión o extracción)
6. Ausencia de transferencia de calor
p V2
Def : Energı́a mecánica total h0 = z + + 2g
ρg
p V2
Energı́a hidráulica (potencial) z+ = h0 −
ρg 2g
Ejemplo 4.12:
Considere el tanque con orificio de la figura. Hallar una relación entre la veloci-
dad de descarga V2 y la altura de la superficie libre h en el tanque. Asumir flujo
sin fricción.
Solución:
141
Ejemplo 4.13:
Solución:
p1 V2 p2 V2 1
+ 1 +
z
1 = + 2 + 2 2
2 ⇒ p1 = p2 + ρ(V2 − V1 )
z
ρg 2g ρg 2g 2
142
Def : Momento angular total de un medio continuo
Z
ρ ~r × ~υ dV
V F (t)
Z Z
d X
ρ (~r × ~v ) dV + (~r × ~v ) ρ (~v − ~vc ) · ~n da = (~r × F~ )O
dt
CV (t) SC(t)
Observación:
Ejemplo 4.14:
143
1. Encuentre una expresión para el torque T que debe ser resistido por el
soporte en A, en términos de las propiedades del flujo en las secciones 1 y
2 y las distancias h1 y h2 .
2. Calcule el valor del torque si D1 = D2 = 3in, p1 = 100 lbf/in2 manométrica,
p2 = 80 lbf/in2 manométrica, V1 = 40 ft/s, h1 = 2 in, h2 = 10 in, y
ρ = 1.94 slugs/ft3
Solución:
Z X X
~ ρ V
~r × V ~ · ~n dA = ~
~r × V ṁout − ~
~r × V ṁin
out in
CS
X
~ A = T~A + ~r1 × (−p1 A1~n1 ) + ~r2 × (−p2 A2~n2 )
M
= ~r2 × V~2 (+ṁout ) + ~r1 × V ~1 (−ṁin )
144
TA + p1 A1 h1 − p2 A2 h2 = ṁ (h2 V2 − h1 V1 )
⇓
π 2
2. A1 = A2 = (3) = 7.07 in2 = 0.0491 ft2
4
10
TA = f t [80(7.07) lbf + 3.81(40) lbf ] −
12
2
f t [100(7.07) lbf + 3.81(40) lbf ]
12
145
Problemas Resueltos
Problema 1
Solución
π 2
1. Conociendo Qi y Ai = d Qi = Q4 = 5 f t3 /s Q2 = Q3 = 2, 5 f t3 /s
4 i
d1 = 6/12 f t ; d2 = 3/12 f t ; d3 = 4/12 f t ; d4 = 5/12 f t
146
2.
Pi V2
3. Hay flujo sin fricción. Es posible aplicar Bernoulli entre: + i =
ρg 2g
Pj Vj2
+
ρg 2g
Entre 1 y 2 Entre 1 y 3
ρ ρ
P1 = P2 + (V22 − V12 ) P1 = P3 + (V32 − V12 )
2 2
ρ 2 ρ 2
(V− V12 ) = P3 +
Igualando obtengo. P2 + (V− V12 )
2 2 2 3
2 (P3 − P2 )
ρ= ≈ 1.95 slug/f t3
V22 − V32
P1 V2 P4 V2
4. Aplicando Bernoulli entre 4 y 1 + 1 = + 4
ρg 2g ρg 2g
1
P4 = P1 + ρ(V12 − V42 ) ≈ 2946 psf
2
Fx +P1 A1 −P2 A2 sin 30◦ −P3 A3 sin 50◦ = ρQ2 (V2 sin 30+V3 sin 50)−ρQ1 V1
147
Problema 2
Nota: Asumir flujo sin fricción y despreciar la diferencia de alturas entre los
puntos 1 y 5.
Solución
π
Q = V1 A1 = V1 D12
4
148
Q π 2
1. Q = cte en todo el sistema, luego Vi = ; Ai = D i =
Ai 4 i
0
1, 2, 3, 4 , 4, 5
2. Tomando el codo , el cual es un tubo de corriente. Aplicamos el Teorema
del momentum.
F̂ = ṁ(V̂out − V̂in ) Es la Ecuación vectorial. Sistema abierto P = Patm
en él.
x → Fx = ṁ(−V5 cos θ − V40 )
y → Fy = ṁ · V5 sin θ
2
~
F = Fx2 + Fy2 = ṁ2 [(V5 cos θ + V40 )2 + (V5 sin θ)2 ] ← Todo es conocido
~
F
ṁ =
[(V5 cos θ + V40 )2 + (V5 sin θ)2 ]1/2
3. Fx y Fy se hallan sustituyendo ṁ
ṁ
4. ṁ = ρQ ⇒ ρ = reemplazando, obtenemos la densidad del fluido.
Q
5. Aplicando Bernoulli entre 4 y i = 1, 2, 3 P4 = Patm P40 = P5 = Patm
P4 V2 Pi V2 1
+ 4 = + i ⇒ Pi = P4 + ρ(V42 − Vi2 )
ρg 2g ρg 2g 2
149
Problema 3
3. El caudal Q3 .
4. La presión en el punto 3, p3 .
5. La densidad de la mezcla, ρ3 .
6. La velocidad en 2, V2 .
150
Solución:
π π
⇒ Q3 = V3 D32 = Q4 = V4 D42
4 4
ṁ3
e) Por definición ṁ3 = ρ3 Q3 → ρ3 = ρ4 = ρ3
Q3
d) Conservación de energı́a (Bernoulli) V C2
P3 V2 P4 V2 1
+ 3 +z
3 = + 4 +z4 → P3 = P4 + ρ(V42 − V32 )
ρ3 g 2g ρ3 g 2g 2
f) Conservación de volumen en V C1 Q2 = Q3 − Q1
π Q3 − Q1
Luego como Q1 = V1 D12 V2 = π 2
4 D
4 2
ṁ1 ṁ2
f) ρ1 = ; ρ2 =
Q1 Q2
151
Problema 4
152
Solución
p1 V2 p2 V2
+ 1 + z1 = + 2 + z2
ρg 2g ρg 2g
2 2 2
V2 − V1 = −2g(z2 − z1 ) + p1
ρ
1/2
2
d2
4 !
2 p1 − 2gH
ρ
V22 1− = p1 − 2gH → V2 =
4
d1 ρ d2
1−
d1
1/2
2
43512 − 2 × 9.8 × 3
V2 = 1000 = 5.35 m/s
4
2.5
1−
7.5
1
V1 = 5.35 × = 0.59 m/s
9
153
Problema 5
Solución
∆p = p1 − p2 = h (γoil − γCO2 )
∆p = 0.1 × 9.8(0.8 × 1000 − 3)
∆p = 781 P a
154
2
π π d2
2. Continuidad V1 d21 = V2 d22 → V1 = V2
4 4 d1
p1 V2 p2 V2 2
+ 1 = + 2 → V22 − V12 = ∆p
ρg 2g ρg 2g ρCO2
1/2
4 !
d2 2 ∆p 2∆p
V22 1− = ⇒ V2 = 4 !
d1 ρCO2
d2
ρCO2 1−
d1
1/2
2 × 781
V2 = 4
! = 27.6 m/s
6
3 1−
8
3. 2
6
V1 = 27.6 × = 15.5 m/s
8
π π
4. ṁ = ρCO2 V1 d21 = 3 × 15.5 (0.08)2 = 0.23 kg/s
4 4
5. Momento lı́neal
π
~2 − V~1
−F + p1 d21 = ṁ V
4
~1 = π (p2 + ∆p)d2 − p2 d2 = π d2 ∆p + π p2 d2 − d2
~2 − V
F + ṁ V 1 2
4 4 1 4 1 2
~1 − π d2 ∆p
F + ṁ V ~2 − V
p2 = 4 1 → p1 = p2 + ∆p
π
p2 (d21 − d22 )
4
π
10 + 0.23(27.6 − 15.5) − × 0.082 × 781
p2 = 4 = 4024 P a (gage)
π
(0.08 − 0.062 )
2
4
155
Problema 6
Solución
1. Manómetro p1 + γH2 O h0 =
P
2 + γHg h
156
p1 = γHg h − γH2 O h0
p1 = 9800(13.6 × 0.25 − 0.05) =
p1 = 32830 P a (gage)
2
π π d2
2. Continuidad V1 d21 = V2 d22 → V1 = V2
4 4 d1
p1 V2 p2 V2 2
+ 1 + z1 = + 2 + z2 → V22 − V12 = −2g(z2 − z1 ) + p1
ρg 2g ρg
2g ρ
Pero: (z2 − z1 ) = H
1/2
2
p1 − 2gH
4 !
2 d2 2 p1 ρ
V2 1 − = − 2gH → V2 =
4
d1 ρ d2
1−
d1
1/2
2
32830 − 2 × 9.8 × 3
V2 = 100 = 2.7 m/s
4
1
1−
2
3. 2
1
V1 = 2.7 = 0.675 m/s
2
π π
4. ṁ = ρV1 d21 = 1000 × 0.675 (0.1)2 = 5.3 kg/s
4 4
5. Momento lı́neal placa (eje y)
−W = ṁ ~out
V ~in = −ṁVin = −ṁV2
− V
157
Problema 7
Tome los datos del problema de acuerdo con el último dı́gito de su cédula:
Solución:
V22 V12
Q̇ − Ẇ = ρQ Cp T2 + − Cp T1 +
2 2
1
Q̇ − Ẇ + ρQ(V12 − V22 ) = Cp (T2 − T1 )ρQ
2
1 1 1
T2 = T1 + Q̇ − Ẇ + ρQ(V12 − V22 ) ×
Cp 2 ρQ
Q̇ − Ẇ
Si desprecio energı́a cinética T20 = T1 +
ρ Q Cp
158
El error es insignificante, del orden del 0.4 %
2. Q2 = 0.556 m3 /s ; V1 = 110.6 m/s ; V1 = 144.5 m/s
◦ ◦
T2 = 18.69 C ; T2 = 19.720 C
El error es insignificante, del orden del 5.5 %
159
Problema 8
Solución
C.M. ṁ = cte
p1 80×144×250×π 18 2
ṁ1 = ρ1 V1 A1 = RT1 V1 A1 = 1715×847×4 12 = 3, 5 slug/s
ṁRT 3,5×1715×2460
V = pπ 2 → V2 = 80×144× π 2 = 725, 3 ft/s
4D 4 ×1,5
160
3,5×1715×2013 3,5×1715×1114
V3 = 40×144× π 2 = 1187 ft/s ; V4 = 14,7×144× π 2 = 4022 ft/s
4 ×1,5 4 ×1
C.E. entre 1 y 2
V2 2 V1 2
Q̇ = ṁ cp T2 + − cp T1 −
2 2
2
2502
725, 3 ft · lbf
Q̇ = 3, 5 6003 × 2460 + − 6003 × 847 − = 34, 7 × 106
2 2 s
1
Q̇ = 34, 7 × 106 778,2 = 44592Btu/s → Recibe el aire
C.E. entre 2 y 3
V3 2 V2 2
−Ẇs = ṁ cp T3 + − cp T2 −
2 2
11822 725, 32
ft · lbf
→ −Ẇs = 3, 5 6003 × 2013 + − 6003 × 2460 − ≈ −7, 85×106
2 2 s
1
Ẇs = 7, 85 × 106 550 = 14266, 5hp → Realizado por el aire sobre la turbina
C.M.L entre 1 y 4 presión manométrica
π π
Fx + p1 D1 2 − p4 D4 2 = ṁ(V4 − V1 )
4 4
π 2 π
Fx = ṁ(V4 − V1 ) − p1 D1 = 3, 5(1022 − 250) − 80 − 14, 7) × 182
4 4
Fx = −3414, 8lbf → F = −Fx = 3414, 8lbf → Gas sobre la turbina
161
Problema 9
Solución
γ 62.4
ρ= = = 1.94 slug/f t3
g 32.2
VC1
162
Por conservación del momento lineal
P~
F = ṁ V¯2 − V¯1 ~1 = 0
V V̂2 = −V2 ĵ
π π
T − W1 = ρ|V2 | D2 (−V2 ) → W1 = T + ρV22 D2
4 4
2 π 2
W1 = 500 + 1.94 · 25.4 · · 0.1 ≈ 510 lbf
4
VC2
163
Problema 10
Solución
γ 62.4
ρ= = = 1.94 slug/f t3
g 32.2
VC1
164
p0 V02 p1 V12
+ + z0 = + + z
1
ρ g 2g
ρg
2g
√ √
V1 = 2gL = 2 · 32.2 · 8 = 22.7 f t/s
VC2
Eje y: FA − w1 − ww = ṁ (−(−V2y ))
FA = 642 lbf
165
Problema 11
Solución
p1 = patm
166
VC
pa V 2 p2 V22
+ 1 + z1 = + + z
2
ρ g 2g ρg 2g
2 pa 1
V12 = 2 g H + → pa = ρ V22 − ρ g H
ρ 2
P~ ¯ − V¯in
Por conservación del momento lineal F = ṁ Vout
~in = 0
V V̂out = −V2 ĵ
π
T − w = ṁ (−V2 ) = −ρ V22 d2
4
r
2 (w − T )4 4(w − T )
V2 = 2
→ V2 =
ρπd ρ π d2
r
4(200 − 10)
V2 = ≈ 41 f t/s
1000 · π · 0.0122
1
pa = · 1000 · 412 − 1000 · 9.8 · 0.425 = 836335 P a
2
167
Problema 12
Solución
V2 = V
Resolviendo el VC
P~ ¯ − V¯in
Por conservación del momento lineal F = ṁ Vout
168
~out = 0
V V̂in = V2 î
π
−F = ṁ (−V2 ) = −ρ V22 d2
4
π 2 2
F = ρ V2 d
4
Aplicando conservación de la energı́a entre 1 y 2 (Bernoulli)
ρ = 780 kg/m3
169
Problema 13
Solución
CM
2 2
π π d1 1
V1 d21 = V2 d22 → V2 = V1 = 25 = 6.25 f t/s
4 4 d2 2
π π
Q = V1 d21 = 25 52 = 490.9 f t3 /s
4 4
170
1. Aplicando conservación de la energı́a entre 1 y 2 (Bernoulli)
p1 V2 p2 V2
+ 1 +
z
1 = + 2 +
z
2 + hf
ρg 2g ρg 2g
p1 − p2 1 226 1
hf = + V 2 − V22 = − + 252 − 6.252 = 5.48 f t
ρg 2g 1 62.4 2 · 32.2
P~
F = ṁ V¯2 − V¯1
Por conservación del momento lineal de 1 a 2
π π
F + p1 d21 − p2 d22 = ρ Q (V2 − V1 )
4 4
γ π π
F = Q (V2 − V1 ) − p1 d21 + p2 d22 =
g 4 4
62.4 π π
F = · 490.9(6.25 − 25) − 720 52 + 946 102 =
32.2 4 4
171
Problema 14
Solución
172
p1 V2 p2 V2
+ 1 +z
1 = + 2 + z
2 + hf
ρg 2g ρg 2g
V2 V2 650000 25.462 9.952
p1
p2 = γ −hf + + 1 − 2 = 9800 −10 + + − =
γ 2g 2g 9800 2 · 9.8 2 · 9.8
p2 = 826924 P a = 826.9 kP a
P~
F = ṁ V¯2 − V¯1
2. Por conservación del momento lineal de 1 a 2
π π
Fx + p2 d22 + p1 d21 cos 60 = ṁ (−V2 − V1 cos 60)
4 4
π 2
Fy + p1 d1 sin 60 = ṁ (−V1 sin 60)
4
π π
Fx = −ṁ (V2 + V1 cos 60) − p2 d22 − p1 d21 cos 60
4 4
π π
Fx = −5000(9.95 + 25.46 cos 60) − 826924 · 0.82 − 650000 · 0.52 cos 60
4 4
Fx = −592871 N ≈ −592.9 kN
π
Fy = −ṁ (V1 sin 60) − p1 d21 sin 60
4
π
Fy = −5000(25.46 sin 60) − 650000 0.52 sin 60
4
Fy = −220773.4 N ≈ −220.8 kN
173
Problema 15
Solución
p3 = p1 = 100 kP a(abs)
p2 = 2340 P a(abs)
p2 V2 p3 V2
+ 2 +
z
2 = + 3 +
z
3 + hf
ρg 2g ρg 2g
174
Y la ecuación de conservación de energı́a queda
p2 2.562 2 p3 V2 V2
+ V3 = + 3 + 0.9 3 →
ρg 2g ρg 2g 2g
p2 − p3 V2
= 1 + 0.9 − 2.562 3 →
ρ g 2g
r r
2 p2 − p3 2 2340 − 100000
V3 = = = 6.48 m/s
ρ 1.9 − 2.562 1000 1.9 − 2.562
π π
Q = d23 V3 = · 0.42 · 6.48 = 0.81 m3 /s
4 4
p1 V2 p2 V2
+ 1 + z1 = + 2 + z2 donde hf ≈ 0 en 1 hasta 2
ρg 2g ρg 2g
H = 4.07 m
175
Problema 16
Solución
176
Se asume flujo sin fricción y z1 ≈ z2
Aplicando CM en VC
Q3 = Q4 + Q5 = 2 f t3 /s
Q2 2
Q3 = Q2 en la caida → V2 = π = 2 = 91.67 f t/s
D 2 π 2
4 2
4 12
Aplicando conservación de la energı́a en la caida del chorro de 2 a 3 (Sin fricción)
Aplicando CM en la Boquilla
2 2
π π D2 1
V1 D12 = V2 D22 → V1 = V2 = 91.67 = 22.92 f t/s
4 4 D1 2
p1 V2 p2 V22
+ 1 + z1 = + + z2
ρg 2g ρg 2g
" 4 #
1 1 D2
→ p1 = ρ V22 − V12 = ρV22 1 −
2 2 D1
4 !
1 1
→ p1 = · 2 · 91.672 1 − = 7878 psf = 54.7 psi
2 2
F = 6834 lbf
177
Problema 17
Solución
178
Aplicando CM en la Boquilla
2
π π D1
V1 D12 = V2 D22 → V2 = V1 = 10 × 22 = 40 f t/s
4 4 D2
p1 V2 p2 V22
+ 1 + z
1 = + + z
2
ρg 2g ρg
2g
1 1
→ p1 = ρ V22 − V12 → p1 = 2 402 − 102 = 1500 psf = 10.42 psi
2 2
P~ P ~out − P ṁin V
~in
Por conservación del momento lineal F = ṁout V
out in
Q3 /2 Q3
V4 = = = V3 = V5
1/2 A3 A3
~4 = V4 (− cos 45◦ , sin 45◦ ) = V3
V ~5 = V5 (cos 30◦ , sin 30◦ ) = V3
V
Q3
CML Fx = ρ[(−V3 cos 45) + (V3 sin 30)] −
ρQ (−V
3 3x )
2
1 1
→ Fx = ρ Q3 V3 − cos 45 + cos 30 = 0.08ρ Q3 V3
2 2
Q3
Fy = ρ [V3 sin 45 + V3 sin 30] − ρ Q3 (−V3 )
2
1 1
→ Fx = ρ Q3 V3 sin 45 + sin 30 + 1 = 1.60ρ Q3 V3
2 2
Fx = 0.08 · 2 · 0.87 · 59.43 = 8.3 lbf
179
Problema 18
Datos: cp = 6000 f t lbf /(slug ◦ R), R = 1715 f t lbf /(slug ◦ R), V1 = 200 f t/s,
Ẇ = 0. Las variables en las secciones 1 y 2 están dadas por la tabla de más
arriba.
Solución
pi × 144
ρi =
R(Ti + 460)
2
π π 1
Q1 = V1 D12 = 200 = 39.27 f t3 /s
4 4 2
π
ρ1 Q1 = ρ2 Q2 ⇒ ρ2 V2 D22 →
4
4 ρ1 Q1
V2 = ṁ1 = 0.177 slug/s = ṁ2
π ρ2 D22
4 · 0.177
V2 = = 737.5 f t/s
π 2.75 × 103 (1/3)2
V2 V2
Q̇ = ṁ2 cp T2 + 2 − ṁ1 cp T1 + 1
2 2
180
737.52 2002
Q̇ = 0.177 6000 · 610 + − 0.177 6000 · 560 +
2 2
f t · lbf
Q̇ = 97693.7 ≈ 125.5 Btu/s
s
p1r = p1 abs −pat = 30−14.7 = 15.3 psi p2r = p2 abs −pat = 20−14.7 = 5.3 psi
π π
Fx + p1r D12 cos 30 − p2r D22 cos 45 = ṁV2 cos 45 − ṁV1 cos 30
4 4
2 2
π 1 π 1
→ Fx + 15.3 × 144 × cos 30 − 5.3 × 144 × cos 45 =
4 2 4 3
= 0.177 · 737.5 cos 45 − 0.177 · 200 cos 30
π π
Fy − p1r D12 sin 30 − p2r D22 sin 45 = ṁV2 sin 45 − ṁV1 sin 30(−1)
4 4
2 2
π 1 π 1
→ Fy − 15.3 × sin 30 − 5.3 × 144 × sin 45
4 2 4 3
= 0.177 · 737.5 sin 45 + 0.177 · 200 sin 30
181
Problema 19
Solución
p1 V2 p2 V2
q
+ 1 +
z
1 = + 2 + z2 ⇒ V2 = V12 − 2gh
ρg 2g ρg 2g
P~ ~out − V
~in )
Conservación momento lineal en el VC, F = ṁ(V
182
En el eje y Vouty = 0 y Viny = V2 . Además la única fuerza que actua sobre el
fluido es el peso de la placa W ~ = −M ĝ. Finalmente, ṁ = ρV1 A
π
q
−M g = ṁ(0 − V2 ) = −ṁV2 = −ρV1 d2 V12 − 2gh
4
2
Mg Mg
q
2 2
= V1 − 2gh ⇒ V1 − 2gh =
ρV1 π4 d2 ρV1 π4 d2
" 2 #
1 2 Mg
h= V −
2g 1 ρV1 π4 d2
" 2 #
1 2 2 × 9.8
h= 10 −
2 × 9.8 1000 × 10 × π4 × 0.0252
h = 4.29 m
183
Problema 20
Solución
p1 V2 p2 V2
q
+ 1 + z1 = + 2 + z
2 ⇒ V 2 = V12 + 2gh = 8.07 m/s
ρg
2g ρg
2g
P~ ~out − V
~in )
Conservación momento lineal en el VC, F = ṁ(V
En el eje y Vouty = 0 y Viny = −V2 . Además la única fuerza que actua sobre el
fluido es la fuerza F~ = F̂. Finalmente, ṁ = ρV1 A
184
π π
q
F = ṁ(0 − (−V2 )) = +ρV1 d2 V2 = ρV1 d2 V12 + 2gh
4 4
π p
F = 1000 × 2.5 × × (0.015)2 2.52 + 2 × 9.8 × 3
4
F = 3.56 N
185
Problema 21
Solución
Flujo estacionario, incompresible y sin fricción (se desprecian los efectos vis-
cosos).
π 2
Conservación de la Masa entre 1 y 2: Q = V1 A1 = V2 A2 con A = 4d
L 10−3 m3 /L 1 3
Q = 1000 × = m /s
min 60 s/min 60
Q 1/60 Q 1/60
V1 = 2
= = 2.12 m/s; V2 = = = 8.49 m/s
πD /4 π0.12 /4 2
πd /4 π0.052 /4
Conservación de la Energı́a entre 1 y 2:
p1 V2 p2 V2
+ 1 +
z
1 = + 2 +
z
2
ρg 2g ρg 2g
p2 p1 1 1
= + (V 2 − V22 ) ⇒ p2 = p1 + ρ(V12 − V22 )
ρg ρg 2g 1 2
186
1
p2 = 200000 + × 1000 × (2.122 − 8.492 ) = 166207 P a
2
P~ ~out − V
~in ) y ṁ = ρV1 π D2
Conservación momento lineal entre 1 y 2, F = ṁ(V 4
F = −1138.4 N
La fuerza que ejerce el agua sobre la contracción es la misma pero en sentido
contrario:
F = 1138.4 N
187
Problema 22
Dos depósitos grandes que contienen agua poseen cada uno un pequeño orificio
de igual área A, tal y como se muestra en la figura. Un chorro de agua emerge
del depósito izquierdo impactando en una placa plana vertical, la cual cubre el
orificio del depósito derecho (ver figura). Asumiendo flujo sin fricción determine
el mı́nimo valor de la altura h en el depósito izquierdo para que la placa no se
mueva horizontalmente.
Solución
p0 V02 p1 V2 p
+ + z0 = + 1 + 1 ⇒ V1 =
z 2gh (1)
ρg 2g ρg 2g
Aplicando la distribución hidrostática de presión se halla p2 = γH (2)
P~ ~out − V
~in )
Conservación momento lineal en VC de la placa, F = ṁ(V
188
En el eje x, Voutx = 0 y Vinx = V1 .La fuerza que actúa sobre la placa es p2 A en
el sentido de derecha a izquierda.
H 30
h= = = 15 f t
2 2
189
Problema 23
Solución
p1 V2 p2 V22 2
+ 1 +
z
1 = + + z2 ⇒ V22 = V12 + p1 − 2gz2
ρg 2g ρg
2g ρ
190
r
2
V2 = (4 × 144) + 32 − 2 × 32.2 × 8 = 8.35 f t/s
2
Conservación de la masa
2
π π 1
Q2 = V2 D22 = 8.35 × × = 0.0455 f t3 /s
4 4 12
2
π π 3
Q1 = V1 D12 = 3 × × = 0.1473 f t3 /s
4 4 12
Q3 0.1018
Q3 = Q1 − Q2 = 0.1018 f t3 /s ⇒ V3 = = = 4.62 f t/s
πD32 /4 π(2/12)2 /4
Conservación de la energı́a entre 1 y 3
p1 V2 p3 V2 1
+ 1 +
z
1 = + 3 + z3 ⇒ p3 = p1 + (V12 − V32 ) − ρgz3
ρg 2g ρg 2g 2
2
p3 = 4 × 144 − × (32 − 4.672 ) − 2 × 32.2 × (−1) = 627.6 psf = 4.36 psi
2
Conservación de momento lineal F~ + F~p1 + F~p2 + F~p3 = ṁ2 V~2 + ṁ3 V~3 − ṁ1 V~1
V~1 = (V1 , 0); V~2 = (V2 cos 30, V2 sin 30); V~3 = (−V3 , 0)
Eje x
π π 2 π
Fx + p1 D12 − p2D2 cos 30 + p3 D32 = ρQ2 V2 cos 30 − ρQ3 V3 − ρQ1 V1
4 4 4
π 2 π
Fx = ρQ2 V2 cos 30 − ρQ3 V3 − ρQ1 V1 − D − p3 D32
4 1 4
2 2
π 3 π 2
Fx = 2 (0.0455·8.35 cos 30−0.1018·4.67−0.1473·3)−4·144 −627.6
4 12 4 12
Fx = −43.1 lbf
Eje y
π 2
Fy − p2D sin 30 = ρQ2 V2 sin 30
4 2
Fy = 2 × 0.0455 × 8.35 × sin 30 = 0.38 lbf
191
Problema 24
Solución:
p1 V2 p2 V2 2∆ p
+ 1 +
z
1 = + 2 +
z 2 2
2 ⇒ V2 − V1 = ; ∆p = p1 − p2
ρg 2g ρg 2g ρ
de la conservación de la masa encontramos
D22
A1 V 1 = A2 V 2 ⇒ V 1 = V2
D12
r
2∆ p
Resolviendo V2 = , luego
ρ (1 − (D2 /D1 )4 )
s
2ρ ∆p
ṁ = ρA2 V2 = A2
(1 − (D2 /D1 )4 )
192
Problema 25
Solución:
193
Q3 4
Q3 = V3 S3 = V1 S1 + V2 S2 = 3 × 1 + 2 × 0.5 = 4 f t3 /s → V3 == = 4 m/s
S3 1
P~ P ~out − P ṁin V
~in
Finalmente, conservación momento lineal: F = out ṁout V in
Eje x
194
Problema 26
Solución:
195
b) Pb1 = F 0 V1 = 313.5 Kw
p1 V2 p2 V2 V22 − V12
+ 1 +
z
1 + hb2 = + 2 +2 → hb2 =
z = 1.06 m
ρg 2g ρg 2g 2g
Pues p1 ≈ p2 y z1 ≈ z2 .
π
Ahora Pb2 = g ṁ hb2 = ρg V1 D12 hb2 = 319.7 Kw
4
e) Pb1 6= Pb2 ya que Pb1 está basado en el movimiento del bote como un todo a
velocidad V1 . Sin embargo, Pb2 está basado en el movimiento del agua a través
V1 + V2
del bote, cuyo valor medio es > V1 por tanto Pb2 > Pb1
2
No obstante Pb1 es una aproximación más sencilla a la potencia de la bomba.
196
Problema 27
Solución:
197
1. aplicando Bernoulli en 1 - 2
P1 v12 P2 v2 1
+ + z1 = + 2 + z2 → P2 = ρgz1 − v22
ρg
2g
ρg 2g 2
Q 0.5
v2 = π = π = 15.9 m/s
D 2 · 0.22
4 4
1000
P2 = 1000 × 9.8 × 45 − 15.92 = 314595 P a = 314.6 kP a
2
2. Aplicando Bernoulli en 2 - C z2 = zc
2
v22
r
PC vC P2 2
+ = + → vC = v22 + (P2 − PC )
ρg 2g ρg 2g ρ
r
2
vC = 15.92 + (314595 − 200000) = 21.95 m/s
1000
π 2 π
QC = vC DC = 21.95 × × 0.152 = 0.39 m3 /s
4 4
3. QA = Q − QC = 0.5 − 0.39 = 0.11 m3 /s
QA 0.11
4. vA = π = π = 14 m/s
2
DA · 0.12
4 4
Aplicando Bernoulli en 2 - A zA = z2
PA v2 P2 v2 1
+ A = + 2 → PA = P2 + ρ(v22 − vA
2
)
ρg 2g ρg 2g 2
198
1
PA = 314595 + 1000 · (15.92 − 142 ) = 343000 P a = 343 kP a
2
5. Momento lineal
π 2 π 2
En x Fx +PA DA −PC DC = ρQA VAx +ρQC VCx = ρQC VC −ρQA VA
4 4
π 2 2
Fx = ρ(QC VC − QA VA ) + (PC DC − PA D A )=
4
π
Fx = 1000(0.39 × 21.95 − 0.11 × 14) + (200000 × 0.152 − 343000 × 0.12 )
4
Fx = 7806.9 N = 7.8 kN → T sobre el agua
Fx0 = −Fx = −7.8 kN ← agua sobre T (fuerza soportada por T)
π π
En y Fy − P2 D2 = −ρQV2y = ρQ V2
4 4
π π 2 π
Fy = ρQ V2 + P2 D = 1000 × 0.5 × 15.9 + 314595 × × 0.22
4 4 4
π
Fy = 1000(0.39 × 21.95 − 0.11 × 14) + (200000 × 0.15 − 343000 × 0.12 )
2
4
Fy = 17833 N = 17.8 kN ↑ T sobre el agua
Fy0 = −Fy = −17.8 kN ↓ agua sobre T (fuerza soportada por T)
199
Problema 28
A través de una tuberı́a bien aislada (los efectos de transferencia de calor con
el entorno son despreciables), se hace circular aire mediante una bomba (ver
figura). Se desea que 50 f t3 /s de aire fluyan a través de la sección A. La presión
absoluta de entrada es pA = 10 psi y la temperatura es 60 ◦ F . ¿Cuál es la
potencia requerida por la bomba en hp? La temperatura del aire a la salida
es de 90◦ F . Tome cp = 6003 f t lbf /slug ◦ r, R = 1716 f t lbf /slug ◦ R, 1 hp =
550 f t lbf /s.
Solución
V2 = 183.3 f t/s
200
Conservación de la energı́a entre 1 y 2 con Q̇ = 0:
V2 V2
Q̇ − ẆS = ṁ ĥ2 + 2 − ĥ1 − 1
2 2
2
p1 10 × 144 π 10
ṁ = V1 A1 = × 91.7 · = 0.0807 slug/s
RT1 1716 · (60 + 460) 4 12
Como ĥ = cp T
V2 V2
ẆS = ṁ cp T1 + 1 − cp T2 − 2
2 2
91.72 183.32
ẆS = 0.0807 · 6003 × 520 + − 6003 × 550 −
2 2
f t · lbf
ẆS = −16764 < 0 Pues la bomba proporciona energı́a al fluido
s
16764
ẆS = = 30.5 hp
550
201
Problema 29
Solución
En el V C que se desplaza con usted, las gotas de agua lo alcanzan con una
velocidad V~ = (−Vx , −Vy ). Planteando la ecuación de conservación de la masa
Z I
d ~ · ~nda = 0
ρd dV + ρd V
dt V C SC
R
Pero md = V C ρd dV y como la lluvia entra solo por el área superior As y
frontal Af
Z Z
dmd dmd
− ρd Vy da − ρd Vx da = 0 ⇒ = ρd Vy As + ρd Vx Af
dt As Af dt
202
ya que D = Vx t. Por consiguiente, dado que ṁs y D son constantes, md (t) será
menor cuanto menor sea t. ¡A correr se ha dicho!
203
Problema 30
Solución
204
" 2 #
1 2 h2 2 h1
F = ρgbh1 1 − − ρh1 bV1 −1
2 h1 h2
205
Problema 31
Solución
206
ṁt
V (t) = −Ve ln 1 − − gt
M0
El primer término es positivo y, si la masa de combustible quemada es una
fracción importante de la masa inicial, la velocidad final del cohete puede ser
mayor que Ve .
207
Problemas Propuestos
Problema 1
Problema 2
Problema 3
208
Problema 4
209
Problema 5
Problema 6
Una bomba extrae agua de un gran tanque, tal y como se muestra en la figura.
¿Cuál es la potencia que debe entregar la bomba al flujo en hp? Suponga flujo
sin fricción en la tuberı́a.
Problema 7
210
Problema 8
Un flujo de alta velocidad en un canal (V1 .h1 ) puede saltar a una condición de
menor energı́a (V2 .h2 ) como se muestra en la figura. La distribución de presión
en las secciones 1 y 2 es aproximadamente hidrostática y la fricción con la pared
es despreciable. Aplique las relaciones de conservación de la masa y momento
para encontrar h2 y V2 en función de (V1 .h1 ).
Problema 9
211
Problema 10
Problema 11
212
Capı́tulo 5
Reactores medioambientales
213
Aplicando el TTR-1 a la parte izquierda del Principio de Conservación de Mása
de la especie A y la expresión anterior se obtiene:
Z Z Z
∂ρA
+ ∇ · (ρA ~υ ) dV = − ∇ · ~jA dV + ω̇A dV
∂t
V F (t) V F (t) V F (t)
Z I Z Z
d
ρA dV + ρA (~υ −~υc )·~nda = − ∇·~jA dV + ω̇A dV = T DA +ṘA
dt
V C(t) SC(t) V F (t) V F (t)
Ejemplo 5.1:
214
Solución:
Entonces
Z I
d dρb
ρb dV + ρb ~υ · ~n da = V + ρb,out Q −
ρb,in
Q = 0
dt dt
VC SC
y reordenando
dρb Q
= −ρb
dt V
Luego, integrando e imponiendo las condiciones iniciales
Q
ρb (t) = ρb0 exp − t
V
Se encuentra que la contaminación decrece exponencialmente. Al comienzo el
suelo se libera del benceno muy rápido pero el proceso es más lento al disminuir
la concentración de benceno.
215
especie se encuentran las ecuaciones empleadas habitualmente en la bibliografı́a.
Es necesario resaltar que los libros normalmente presentan las ecuaciones de los
reactores mediambientales a partir de un balance de masa especı́fico para cada
tipo. En la presentación actual éstas se derivan de la expresión general im-
poniendo las condiciones de operación de cada tipo de reactor.
Z I I Z
d
C dV + C (~v − ~vc ) · ~n da = Jd da + r dV
dt
V C(t) SC(t) SC(t) V C(t)
216
Existen tres tipos de reactores ideales
Z Z
d d dC
C dV = r dV → V C = rV → = r(C)
dt dt dt
V C(t) V C(t)
Este reactor se concibe esencialmente como un tubo donde los reactivos avanzan
a velocidad constante y su transformación en productos se produce a medida que
217
aquellos progresan en el reactor. En esta configuración se asume mezcla com-
pleta en dirección radial pero no en dirección axial. Este reactor opera siempre
en estado estacionario. En estas condiciones la ecuación integral consiste en la
igualdad de los términos de ingestión y reacción
Z I Z
r dV = C ~v · ~n da = ∇ · (C ~v ) dV
VC SC VC
Z Z
1D dC
∇ · (C~v ) dV = r dV → ∇ · (C~v ) = r −→ U =r
dx
UC VC
218
Z I Z
d
C dV + C v̄ · n̄ da = r dV
dt
VC SC VC
Los tres tipos de reactores operan con reacciones de diferentes de diferente orden.
Habitualmente se consideran tres tipos de orden cero, uno y dos atendiendo al
exponente de C en la reacción:
Obsérvese que como las dimensiones de la reacción, [r] = M L−3 t−1 , son fijas,
las dimensiones de la constante de la reacción varı́an dependiendo del orden
• Reactor Batch
i) C(t) = C(0) − k0 tC(t) = 0 t ≤ C(0)/k0
C(t) = 0 t > C(0) k0
dC
= r(C) →
dt
ii) C(t) = C(0) e−k1 t
C(0)
iii) C(t) = (Ejercicio)
1 + 2k2 t C(0)
C(0)
Se define el tiempo caracterı́stico Θ del reactor como: Θ=
r (C(t = 0))
• Reactor PFR
219
k0 Cin U
1. C(x) = Cin − x x≤
U k0
Cin U
C(x) = 0 x>
k0
k1 x
−
2. C(x) = Cin e U
Cin
3. C(x) = (Ejercicio)
1 + 2k2 Cin x/U
L A V V
× = Θ= : Tiempo medio de residencia
U A Q Q
Θ es el tiempo caracterı́stico del reactor y su parámetro de operación más im-
portante. Las soluciones para los tres tipos de reacciones se escriben:
Cin
1. C(L) = Cin − k0 Θ Θ≤
k0
Cin
C(L) = 0 Θ>
k0
2. C() = Cin e−k1 Θ
Cin
3. C(L) =
1 + 2k2 Cin Θ
• Reactor CMFR
P
Obsérvese que en este reactor Qout = Q el caudal que circula por el sistema,
out
luego:
dC 1 1 X V
= r(C) − C + Qin Cin ; Θ=
dt Θ V in Q
Como el CMFR puede operar en modo transitorio o estacionario, escribamos
las soluciones en el último caso donde dC
dt = 0
1 P
1. C = Qin Cin − k0 Θ
Q in
220
P
Qin Cin
1 in
2. C =
Q (1 + k1 Θ)
s
P
−1 + 1 + 8k2 Θ Qin Cin /Q
in
3. C = (Ejercicio)
4k2 Θ
1 1 P
• Sin reacción A = ; B= Qin Cin
Θ V in
1 1 P
• Orden 0: A = ; B = −k0 + Qin Cin
Θ V in
1 1 P
• Orden 1: A = k1 + ; B= Qin Cin
Θ V in
• Orden 2: La ecuación es no lineal y se conoce como ecuación de Riccati
(ver siguiente sección)
Ecuación diferencial:
" #
dC 1 X X
= −2k C 2 + Qin Cin − C Qout
dt V in out
221
Esta es una ecuación diferencial cuadrática de primer orden, la cual se conoce
como ecuación de Riccati. Se expresa en la forma general:
dC
= A (t) C 2 + B (t) C + D (t)
dt
En nuestro caso los coeficientes A, B y D son constantes (no dependen de t).
Identificado términos en el caso del reactor CMFR:
1 P 1 1 P
A = −2k < 0 B=− Qout = − < 0 D= Qin Cin
V out Θ V in
los cuales se pueden evaluar conociendo los caudales, las concentraciones a la
entrada y la constante de la reacción.
dZ
+ (B + 2ACp ) Z = −A
dt
Realizando la identificación E = B + 2ACp y F = −A, se obtiene la forma
canónica de la ecuación de primer orden cuya solución es conocida:
F
Z(t) = Z0 e−Et + 1 − e−Et
E
222
Sustituyendo, la expresión para la función Z(t) en términos de los coeficientes
A y B es la siguiente:
1 A h i
Z(t) = e−(B+2ACp )t − 1 − e−(B+2ACp )t
C0 − Cp B + 2ACp
Esto finaliza la presentación de las soluciones de los reactores ideales con reac-
ciones hasta orden dos.
223
Problemas Resueltos
Problema 1
Solución
224
Problema 2
Dos tuberı́as, cada una de 10 cm2 de área transversal, introducen agua en una
cámara de mezcla, mientras que una tercera la saca de ella. La primera trans-
porta agua saturada en oxı́geno (C1 = 9 mg/L), y la segunda transporta agua
desoxigenada (C2 = 0 mg/L). Un agitador en la cámara mezcla rápidamente las
dos corrientes de tal modo que la concentración en el tanque es espacialmente
uniforme. Asumiendo que el sistema se encuentra en estado estacionario halle
el valor de la concentración en la salida C3 .
Datos: V1 = 20 cm/s, V2 = 10 cm/s, A3 = 30 cm2 .
Solución
Q1 + Q2 = Q3 → Continuidad fluido
Conservación de masa Q2
C1 Q1 + C2 Q2 = C3 Q3 = C3 (Q1 + Q2 )
C1 Q1 + C2 Q2
C3 =
Q1 + Q2
Q1 = V1 A1 = 10 × 20 = 200 cm3 /s
Q2 = V2 A2 = 10 × 10 = 100 cm3 /s
225
Problema 3
Solución
Qout = Qin + Q1 + Q2 = 50 + 3 + 10
Qout = 63 MGD
Conservación E. Coli
226
Problema 4
Nota: El caudal de trabajo del CMFR es el mismo que el del PFR y la reacción
quı́mica de segundo orden en ambos reactores posee la misma constante k2 .
Nota: Las modificaciones preferidas al diseño son aquellas que no implican
cambio de las dimensiones de los reactores.
Indicación: Tenga mucho cuidado con las unidades
Solución:
227
8k2 Θ0 Cin + 1 − 1
r
V 1000
3. C = CMFR Θ0 = = = 2000 s
4k2 Θ0 Q 0.5
r
8 × 2 × 2000 × 0.02 + 1 − 1
C= = 1.52×10−3 kg/m3 = 1.52×10−3 mg/L
4 × 2 × 2000
228
Problema 5
Usted está diseñando un reactor que utiliza cloro en un reactor de flujo uniforme
PFR para desinfectar patógenos en agua potable. Se requiere un tiempo de
residencia mı́nimo de 30 minutos para reducir la concentración de patógenos
desde 100 patógenos/L a un valor de 2 patógenos/L, a través de una reacción
de segundo orden (r(C) = −2k2 C 2 ). Usted planea utilizar un flujo volumétrico
de 1000 galones/min.
Solución:
Cin
Cin Cin −1
1. C(L) = → 2k2 Cin Θ = −1 → k2 = C
1 + 2k2 Cin Θ C 2Cin Θ
100
−1 L
Con Θ = 30 min k2 = 2 ≈ 8.2 × 10−3
p p · min
2 × 100 × 30
L
V
Θ= → V = ΘQ = 30 × 1000 = 30000 gal
Q
2. CMFR
r
8k2 Θ Cin + 1 − 1
C= → (4k2 Θ C + 1)2 = 8k2 Θ Cin + 1
4k2 Θ0
229
16k22 Θ2 C 2 + 8k2 Θ C + 1 = 8k2 Θ Cin + 1
Cin − C 100 − 2
2k2 Θ C 2 + C = Cin → Θ = = =
2C 2 k2 2 × 22 × 8.2 × 10−3
1494 min
230
Problema 6
Tome los datos del problema según el último dı́gito de su código de estudiante:
Dı́gito Cin [mg/L] Cmax [mg/L] V [m3 ] Qmin [m3 /s] k [m3 /(kg s)]
1-2-3 100 5 10000 0.2 25 × 10−3
4-5-6 50 2 5000 0.1 10 × 10−3
7-8-9-0 200 10 15000 0.3 15 × 10−3
Solución:
V
Para Qmin → Θmin =
Qmin
CMFR PFR
L 103 L/m3
Debo pasar k2 a [k2 ] = multiplicándola por 6
mg · s 10 mg/kg
L
1. k2 = 25 × 10−6 = 2.5 × 10−5 ; Θmin = 5 × 104 s = 13.89 h
mg · s
CCM F R = 6.13 mg/L CP F R = 0.40 mg/L
Elijo el PFR, pues CCM F R > Cm x
L
2. k2 = 10 × 10−6 = 10−5 ; Θmin = 50000 s = 13.89 h
mg · s
CCM F R = 6.59 mg/L CP F R = 0.98 mg/L
Elijo el PFR, pues CCM F R > Cm x
3. k2 = 15 × 10−6 = 1.5 × 10−5 mg·s
L
; Θmin = 50000 s = 13.89 h
CCM F R = 11.22 mg/L CP F R = 0.66 mg/L
Elijo el PFR, pues CCM F R > Cm x
231
Problema 7
Solución
C0 = 1 g/L C = 0, 2 g/L
Q
− +kt
dC Q V
=− + k0 ·C → C(t) = C0 · e → Solución
dt V
Resuelvo para t
C(t) Q 1 C
ln =− +k ·t → t=− ln
C0 V Q C0
+k
V
1
t=− ln 0.2 = 193562 s ≈ 2.24 day
3 0.2
+
5 × 105 86400
232
Problema 8
Solución
233
Problema 9
Solución
dC Q dC Q
Cin − C) − k C
= ( → =− +k C
dt V dt V
234
Problema 10
Solución
Q1 + Q2 100 + 150
A = k1 + = 0.02 + = 0.045 S −1
V 100000 L
Q1 C1 + Q2 C2 100 · 5 + 150 · 10
B= = = 0.2 mg/L s
V 10000
235
B
C(t) = C0 e−At + 1 − e−At
Sustituyendo en la solución
A
0.2
C(t) = 12 e−0.045t + 1 − e−0.045t = 7.56 e−0.045t + 4.44
0.045
Si Cf = 8 mg/L
8 − 4.44
8 = 7.56 e−0.045t + 4.44 → = e−0.045t →
7.56
ln (0.47)
ln (0.47) = −0.045t → t=− = 16.8 s
0.045
236
Problema 11
Solución
Q1 + Q2 60 + 120
A = k1 + = 0.5 + = 0.68 min−1
V 1000 L
Q1 C1 + Q2 C2 60 · 1 + 120 · C2 g
B= = = 0.06 + 0.12 C2
V 1000 L L min
237
B
C(t) = C0 e−At + 1 − e−At
Sustituyendo en la solución
A
0.06 + 0.12 C2
1 − e−0.68 t
C(t) =
0.68
Pues C0 = 0
Si Cf = 3 g/L en 5 min
0.06 + 0.12 C2
1 − e−0.68·5 → 3 = (0.06 + 0.12 C2 ) 1.42
3=
0.68
1 3
C2 = − 0.06 = 17.1 g/L
0.12 1.42
238
Problema 12
Solución
dC Q
Ecuación diferencial = −2 k1 C 2 + (
Cin − C)
dt V
donde Q = 0.1 f t3 /s × 86400 ≈ 8640 f t3 /day
Q 1 8640
A = −2k = −0.16 f t3 /mg day B=− =− =− = −8.64 day −1
V Θ 1000
Solución partı́cular:
Solución: t en dı́as.
−1
1 A
e−(B+ 1 − e−(B+
2Acp )t 2A
cp )t
C(t) = cp + −
C0 − cp B + 2Acp
−1
1 8.64 t −0.16
C(t) = e − 1 − e8.64 t =
100 −8.64
−1
C(t) = 0.01 e8.64 t − 0.0185 1 − e8.64 t =
−1
C(t) = −0.0185 + 0.0285 e8.64 t
239
Hallando el tiempo t para que C = C0 /2 = 50 mg/L
1 0.0385
t= ln ≈ 0.0348 day ≈ 0.835 horas ≈ 50 min
8.64 0.0285
240
Problema 13
Solución
dC Q
Ecuación diferencial = −2 k C 2 + (Cin − C)
dt V
donde Q = 250 f t3 /h × 24 = 6000 f t3 /day
Q 1 6000
A = −2k = −0.36 f t3 /mg day B=− =− =− = −0.12 day −1
V Θ 5 × 104
V Q mg
Θ= = 8.33 day D = Cin = 0.12 × 80 = 9.6
Q V L day
Solución partı́cular:
√ √
1 + 8k ΘCin + 1 1 + 8 · 0.18 · 8.33 · 80 + 1
Cp = − =− = 5 mg/L
4kΘ 4 · 0.18 · 8.33
−1
C(t) = 5 + −0.097 + 0.119 e3.72 t
241
Hallando el tiempo t para que C = C0 /2 = 25 mg/L
242
Problema 14
Solución
0.15
C(t) = 12e−0.045t + (1 − e−0.045t ) → 6 = 3.33 + 8.67e−0.045t
0.045
| {z }
3.33
Resolviendo para t
6 − 3.33 −1 6 − 3.33
ln = −0.045t ⇒ t = ln = 26.17 s
8.67 0.045 8.67
243
Problema 15
Dos tuberı́as, cada una de 20 cm2 de área transversal, introducen agua en una
cámara de mezcla, de volumen V , mientras que una tercera la saca de ella.
La primera transporta cierto quı́mico venenoso con una cierta concentración, y
la segunda transporta agua limpia (C2 = 0 mg/L). Un agitador en la cámara
mezcla rápidamente las dos corrientes de tal modo que la concentración en el
tanque es homogéna. Adicionalmente, el quı́mico decae siguiendo una reacción
de segundo orden de constante k2 = 10−6 L/(mgs). La concentración a la salida
es C3 . Asumiendo que el sistema se encuentra en estado estacionario halle el
valor de la concentración en la primera entrada C1 .
Solución
Q1 = V1 A1 = 20 × 10 = 200 cm3 /s
Q2 = V2 A2 = 20 × 15 = 300 cm3 /s
Q = Q1 + Q2 = 500 cm3 /s
Se trata de un reactor CMFR estacionario, luego (C = C3 ):
Q Q1 C1 +
Q2
C2
0 = −2k2 C 2 − C+
V V
Q QC + 2k2 C 2 V
2k2 C 2 V + C = Q1 C1 ⇒ C1 =
V Q1
Cambio de unidades: V = 5 m3 = 5000 L, Q = 500 cm3 /s = 0.5 L/s,
Q1 = 200 cm3 /s = 0.2 L/s
244
Problema 16
Solución
p P p
−1 + 1 + 8k2 Θ( in Qin Cin )/Q −1 + 1 + 8 · 2.5 × 10−6 · 104 · 2/0.5
Cp = =
4 k2 Θ 4 · 2.5 × 10−6 · 104
Cp = 3.416 mg/L
Solución: t en s. B + 2ACp = −10−4 + 2(−5 × 10−6 ) × 3.416 = −1.34 × 10−4
245
−1
1 −(B+2A Cp )t A −(B+2A Cp )t
C(t) = Cp + e − 1−e
C0 − Cp B + 2A Cp
1 −4 −5 × 10−6 −1
−(−1.34×10−4 ) t
C(t) = 3.416+ e−(−1.34×10 ) t − 1 − e
2 − 3.416 −1.34 × 10−4
h −4
−4
i−1
C(t) = 3.416 + −0.706 e1.34×10 t − 0.037 1 − e1.34×10 t
h −4
i−1
C(t) = 3.416 + −0.037 − 0.669 e1.34×10 t
Como 1 h = 3600 s
246
Problema 17
Solución
p P p
−1 + 1 + 8k2 Θ( in Qin Cin )/Q −1 + 1 + 8 · 0.02 · 40 · 1500/250
Cp = =
4 k2 Θ 4 · 0.02 · 40
Cp = 1.65 mg/L
Solución: t en s. B + 2ACp = −0.025 + 2(−0.04) × 1.65 = −0.157
−1
1 A
e−(B+2A Cp )t − 1 − e−(B+2A Cp )t
C(t) = Cp +
C0 − Cp B + 2A Cp
247
−1
1 −0.04
C(t) = 1.65 + e0.157 t − 1 − e0.157 t
1 − 1.65 −0.157
0.157 t
−1
C(t) = 1.65 + −1.54 e − 0.25 1 − e0.157 t
−1
C(t) = 1.65 + −0.25 − 1.29 e0.157 t
248
Problemas Propuestos
Problema 1
Problema 2
Con el objeto de desinfectar un cierto volumen de agua se inyecta cloro con una
concentración de 1 g/L. El cloro reacciona con los patógenos y decae siguiendo
una cinética de primer orden de constante k = 0.18 day −1 . Si el agua circula
con un flujo volumétrico de 300 m3 /day y el tratamiento de desinfección no
debe durar más de 6 horas, encuentre el volumen de agua tratada si la mı́nima
concentración permisible de cloro a la salida es de 0.3 g/L.
249
Problema 3
¿Cuál es el flujo volumétrico máximo en L/s que se puede tratar si la salida con-
tiene un 10% de los metales pesados contenidos en el flujo de entrada? ¿Cómo
se modifica dicho flujo volumétrico si la constante de la reacción se dobla?
Problema 4
Problema 5
250
Datos: k = 0.06 s−1 , Q = 100 L/s.
Problema 6
Problema 7
Problema 8
Dos tuberı́as, cada una de 10 cm2 de área transversal, introducen agua en una
cámara de mezcla, de volumen V , mientras que una tercera la saca de ella.
La primera transporta cierto quı́mico venenoso (C1 = 9 mg/L), y la segunda
transporta agua limpia (C2 = 0 mg/L). Un agitador en la cámara mezcla
rápidamente las dos corrientes de tal modo que la concentración en el tanque
es espacialmente uniforme. Adicionalmente, el quı́mico decae siguiendo una
reacción de segundo orden de constante k2 = 10−6 L/(mg s). Asumiendo que el
sistema se encuentra en estado estacionario halle el valor de la concentración en
la salida C3 .
Problema 9
251
concentración C del contaminante venenoso en la cámara de tratamiento si ésta
se puede considerar como un reactor completamente mezclado. Con ayuda de
las fórmulas vistas en clase, escriba la solución para C en función del tiempo
si la concentración inicial en la cámara es C(t = 0) = C0 . Grafique aproxi-
madamente la curva C(t) y encuentre el tiempo que ha transcurrido para que
la concentración en la cámara sea de 10 mg/L.
Problema 10
Suponga que la PTAR de la UAO se puede modelar como un reactor tipo CMFR
de volumen V por el que fluye un cierto caudal Q. El agua se trata mediante
una reacción de segundo orden de constante k con el fin de reducir la concen-
tración de contaminantes desde un valor inicial C0 hasta un valor final C(t) en
un cierto tiempo. Con los datos proporcionados más abajo escriba la función
que expresa la evolución temporal de la concentración de contaminantes en la
PTAR C(t). ¿En qué unidades debe expresar el tiempo?
Problema 11
Problema 12
252
0 mg/L). Un agitador en la cámara mezcla rápidamente las dos corrientes de
tal modo que la concentración en el tanque es homogénea. En el reactor el
patógeno es sometido a un proceso quı́mico que emplea una reacción de se-
gundo orden de constante k2 . Los caudales en las tuberı́as de entrada son Q1
y Q2 , respectivamente. Si la concentración inicial del patógeno es C0 , halle el
valor de la concentración en el reactor después de 15 s. Asuma el reactor como
un CMFR.
Problema 13
Problema 14
Problema 15
253
(en L/s) que debe utilizarse para reducir la salinidad del agua del estanque a
500 mg/L en un año. Asuma el estanque como un CMFR.
254
Capı́tulo 6
Razones:
Insuficiente caracterización
Ecuaciones constitutivas o propiedades fı́sicas desconocidas (flujos turbulentos,
fluidos no newtonianos...)
Solución inabordable
Complejidad de las ecuaciones básicas (matemática y/o numérica). Necesidades
de computo desbordantes (flujos turbulentos)
Costo excesivo
Incluso si la solución matemática es posible, quizás, demande una elevada in-
versión de recursos y tiempo que lo conviertan en inabordable para una industria
en particular
Necesidad de fiabilidad
Incluso si el problema está caracterizado convenientemente, el cómputo es accesi-
ble y rentable su solución es inevitable confrontar los resultados con las medidas
experimentales (verificación de predicciones, acotar o eliminar incertidumbres o
errores en el cálculo → industria aeronáutica y nuclear, por ejemplo)
255
• Tabulación dispositivos y componentes catalogados: renuncia a una predicción
racional y económica de un fenómeno ⇒ ensayos sistemáticos
• Ensayo de prototipos: Desarrollo y/o caracterización productos novedoso
únicos ⇒ caracterización experimental del diseño
1. Economı́a
• Ejecución en escala conveniente
• Parámetros significativos a medir reducidos
• Número de medidas mı́nimo pero suficiente
• Equipos y personal necesario a coste razonable
2. Fiabilidad
• Caracterización de los fenómenos esenciales
• Errores de medida acotados
• Reproducibilidad
3. Generalidad
• Situación experimental caracteriza diversas situaciones
• Resultados extrapolables a sistemas semejantes
• Resultados conducentes a formular leyes de comportamiento
Por ejemplo, la fuerza F que actúa sobre un cuerpo inmerso en un fluido depende
de la longitud del cuerpo L, la velocidad del flujo V , la densidad del fluido ρ y
su viscosidad µ , es decir,
F = f (L, V, ρ , µ) (1)
Si se necesitan ∼ 10 puntos experimentales para determinar una curva y F
depende de 4 variables ⇒ un ensayo sistemático requiere del orden de 104
256
medidas experimentales para caracterizar suficientemente el sistema.
Sin embargo, con el análisis dimensional, se puede reducir (1) a la forma equiva-
lente.
F ρV L
= g CF = g(Re)
ρ V 2 L2 µ
ρV L
donde CF es el Coeficiente de fuerza y Re el número de Reynolds Re =
µ
los cuales son números adimensionales
257
6.2. Principio de homogeneidad dimensional
Enunciado
ϕ(π1 , . . . , πN −M ) = 0
Desarrollo de la adimensionalización
258
efecto del flujo mediante un conjunto de ensayos de laboratorio.
M L t
ρ 1 -3 0
R 0 1 0
U∞ 0 1 -1
FD 1 1 -2
µ 1 -1 -1
259
El rango de la matriz es 3 ya que la submatriz cuadrada con determinante
distinto de cero es 3x3
1 −3 0
0 1 0 = −1 6= 0
0 1 −1
Recomendaciones:
• Se excluirá aquella variable que se desee obtener como dependiente del
resto, ya que sino dicho parámetro aparecerá en todos los grupos adimen-
sionales, generando una expresión implı́cita inadecuada para el propósito
inicial.
• Se excluirán variables sospechosas de no intervenir siempre para evitar
la influencia en los grupos adimensionales de un parámetro tal vez no
significativo.
• Aplicados ambos criterios, “como primeras variables” se eligen aquellas
cuyas dimensiones sean lo mas sencillas posible.
• Finalmente, la costumbre o experiencia aconsejará un tipo de elección
frente a otra.
Sin embargo, sea cual sea la alternativa elegida para las “primeras variables”
esta conducirá a resultados equivalentes.
260
“Primeras variables” = {R, U∞ , ρ}
Etapa 4: Cálculo números adimensionales
M 1 0 0 α 1
L −3 1 1 · β = 1
t 0 0 −1 γ −2
M α=1 → α=1
FD
L 1 = −3α + β + γ → β=2 ⇒ πF D =
ρ R2 U∞
2
t −2 = −γ → γ=2
Para el 2◦ parámetro:
0
M 1 0 0 α 1
L −3 1 1 · β 0 = −1
t 0 0 −1 γ0 −1
M α0 = 1 → α0 = 1
µ
L −1 = −3α + β 0 + γ 0 → β0 = 1 ⇒ πµ =
ρ R U∞
t −γ 0 = −1 → γ0 = 1
ρ R U∞
como Re = Número de Reynolds
µ
despejando πFD , se puede escribir
2 FD 2 ρ R U∞
= f
π ρ R2 U∞
2 π µ
donde AT = πR2 es el área transversal al flujo. Entonces
FD 2 FD 2
1 2
= f (Re) ⇒ cD (Re) = 1 2
= f (Re)
2 ρ AT U∞ π 2 ρ AT U∞ π
261
1 2
FD = cD (Re)ρ AT U∞
2
Ejemplo 6.1:
Solución:
S S0 V0 g t
L 1 1 1 1 0 rango R = 2
t 0 0 −1 −2 1
262
# magnitudes fundamentales P=2
# parámetros fundamentales N=5
# números adimensionales M=N-R=3 ⇒ φ(π, π, π3 )
S g t
πS = πg = 0 πt = 00
V0β V0β V0β
0 00
S0α S0α S0α
S
L 1=α+β → α=1 ⇒ πS =
S0
t 0 = −β → β0 = 0
g S0
L 1 = α0 + β 0 → α0 = −1 ⇒ πg =
V02
t −2 = −β 0 → β0 = 2
t V0
L 0 = α00 + β 00 → α00 = 1 ⇒ πt =
S0
t 1 = −β 00 → β 00 = −1
Sg
L α+β =1 → −2β + β = 1 → β = −1 ⇒ πS =
V02
t α + 2β = 0 → α = +2
S0 g
Análogamente πS0 =
V02
tg
L α0 + β 0 = 0 → α0 = −β 0 → α0 = 1 ⇒ πt =
V0
t −α0 − 2β = 1 → β 0 = −1
263
III) Primeras Variables {S0 , g}
r
S V0 g
πS = πV 0 = √ πt = t
S0 g S0 S0
6.4. Modelos y semejanza
Consecuencias de la adimensionalizacion
II) Universalidad
La evaluación de la relación adimensional, por ejemplo cD (Re) se realiza experi-
mentalmente. Tal representación es universal, en el sentido de que para cualquier
combinación de parámetros que den el mismo Re le corresponde idéntico valor
cD ⇒ un único conjunto de experimentos, sirve para caracterizar todos los
tipos de esferas a cualquier velocidad y en cualquier fluido.
III) Economı́a
El numero neto de ensayos se reduce dramáticamente. En nuestro ejemplo, 10
en vez de 104 . Como cD (Re) es independiente del mecanismo que se use para
variar Re se puede utilizar en el ensayo el procedimiento más económico, en
nuestro caso variar la velocidad del fluido.
264
y la semejanza dinámica implicará Rem = Rep ⇒ cD m = cD p y dividiendo
2
cD p FD p (ρ U∞ AT ) p
=1 ⇒ = 2
cD m FD m (ρ U∞ AT )m
que permite evaluar directamenteFD p en función de FD m .
Ejemplo 6.2:
Un avión tiene una longitud caracterı́stica de 125 f t y está diseñado para volar a
300 mi/h a 10000f t de altura estándar. El coeficiente de resistencia aerodinámica
se mide en un modelo y se encuentra que es igual a 0.01. Si se satisfacen las
condiciones de semejanza dinámica ¿cuál será la potencia, en hp, desarrollada
por el motor?
Solución:
π0 = π0 (π1 , · · · , πM )
265
Si π1p = π1m , · · · , πM p = πM m ⇒ π0p = π0m ⇒ semejanza dinámica
| {z }
M ecuaciones con N parámetros adimensionales
Puede ocurrir también que las condiciones experimentales deducidas sean irrea-
lizables (valor desproporcionado, coste excesivo, desbordamiento rango opera-
tivo de la instalación, inexistencia del fluido con propiedades fı́sicas determi-
nadas, ...).
Rem = Rep y F rm = F rp
Como la escala y propiedades del barco vienen dadas en el ensayo, resulta im-
posible verificar a la vez ambas relaciones. La velocidad del ensayo U∞ es el
único parámetro disponible.
Como solo se puede conservar uno de los números adimensionales la ley de escala
vendrá dada por
2 1/2
(ρ U∞ AT )p cD p
FD p = FD m 1/2
(ρ U∞ AT )m |cD
2
{zm}
6=1
donde cD p 6= cD m pues F rp 6= F rm .
266
Por tanto en semejanza parcial no es indiferente la elección de perı́metros a con-
servar. En el ejemplo cD (Re) → 0.44 si Re >> 1, luego los efectos de escala
se minimizarán si se conserva F r y se renuncia a Re.
Como resultado, los valores obtenidos para el prototipo a partir de las medi-
das en el modelo están sujetas a incertidumbres que es preciso acotar y controlar.
Tipos de semejanzas:
Ejemplo 6.3:
Solución:
267
Pm Fm V m 1
= =
Pp Fp Vp 254000
6.5. Adimensionalizacion de las ecuaciones básicas. Números
adimensionales más comunes en Mecánica de Fluidos
En esta sección se aplican las técnicas vistas a las ecuaciones básicas, por ejem-
plo, a las ecuaciones de flujo incompresible newtoniano e isotermo:
∇ · ~υ = 0 Ec. Continuidad
D~υ
ρ = ρ ~g − ∇p + µ∇2 ~υ Ec. Momento lineal
Dt
Definimos:
~υ x y z tU p + ρ gz
~υ ∗ = ; x∗ = ; y∗ = ; z∗ = ; t∗ = ; p∗ =
U L L L L ρU 2
y como ρ, U, L son constantes, las derivadas espaciales se pueden escribir:
∂u ∂(U u∗ ) U ∂u∗
= ∗
= donde ~υ = (u, υ, w)
∂x ∂(Lx ) L ∂x∗
Entonces la forma adimensional de las ecuaciones queda:
∇∗ · ~υ ∗ = 0
D~υ ∗ µ
ρ ∗ = −∇∗ p∗ + ∇∗ 2 ~υ ∗
Dt ρ UL
∗ ∂ ∂ ∂
con ∇ = , , .
∂x∗ ∂y ∗ ∂z ∗
Se observa que el único grupo adimensional no defnidos antes que aparece es el
ρU L
número de Reynolds Re = .
µ
Los números adimensionales más importantes que aparecen en Mecánica de Flu-
idos se recogen en la Tabla 1 adjunta.
268
El número de Reynolds y regı́menes de flujo
269
Parámetro Definición Significado Importancia
ρU L Inercia
Número de Reynolds Re = Siempre
µ V iscosidad
U V elocidad f lujo
Número de Mach Ma = Flujo compresible
a V elocidad sonido
U2 Inercia
Número de Froude Fr = Superficie libre
gL Gravedad
ρU 2 L Inercia
Número de Weber We = Superficie libre
σ T ensión superf icial
µcp Disipación
Número de Prandtl Pr = Convección de calor
k Conducción
cp Entalpia
Razón calores especı́ficos k= Flujo compresible
cv Energia interna
ωL Oscilación
Número de Strouhal St = Flujo oscilante
U V elocidad media
Rugosidad pared
Rugosidad relativa Turbulencia, pared rugosa
L Longitud cuerpo
p − pv P resión
Número de cavitación Ca = Cavitación
ρU 2 Inercia
p − p∞ P resión
Coeficiente de presión Cp = Aerodinámica
ρU 2 /2 Energia cinética
FD Resistencia
Coeficiente resistencia CD = Aerodinámica
ρU 2 A/2 F uerza dinámica
hf F ricción
Factor de fricción f= Flujo interno
(V 2 /2g)(L/d) Altura velocidad
τw F ricción
Coef. fricción superficial Cf = Flujo capa lı́mite
ρV 2 /2 Energia cinética
270
Problemas Resueltos
Problema 1
271
Solución:
ρLV V 2
R
Primeras variables {ρ, L, V } =f ,
ρV 2 L2 µ gL
!
R µ V2
{ρ, L, g} =f ,
gL3 ρ
p
gL3 ρ gL
ρLV V2 R
Donde = Re ; = Fr ; = CD
µ gL ρV 2 L2
2. CD , Re , Fr
2
V2
R 0 V R
3. 2 2
=f ó = f0
ρV L gL gL3 ρ gL
Vp2 Vm2 Lm
4. Fr m = Fr p = → Vm2 = Vp2
gp Lp gm Lm Lp
Vp
Vm = √ = 0.179 m/s
20
Rp Rm Vp2 L2p
5. CD m = CD p = → R p = Rm = 20(20)2 Rm
ρp Vp2 L2p ρm Vm2 L2m Vm2 L2m
Rp = 17.6 KN
R R
Igual sale si =
gL3 ρp gL3 ρ m
3
Lp
Rp = Rm = Rm (20)3 = 17.6 KN
Lm
272
Problema 2
Solución:
FD = f (d, D, V, ρ, µ)
FD FD
2. ó = Coeficiente de fuerza
ρV 2 D2 ρV 2 d2
ρV D ρV d d D
ó = Reynolds ; ó = Diámetro relativo
µ µ D d
273
FD d FD D
3. = g10 ó = g20
ρV 2 D2 D ρV 2 d2 d
d1 d2 d1 0.2 d1 D1 1
4. = → d2 = D2 = · 2 = 0.8 f t = =
D1 D2 D1 0.5 d2 D2 48
Semejanza dinámica
FD 2 ρ V 2 D2 FD 2 ρ V 2 d2
= 22 22 ó = 22 22
FD 1 ρ
V1 D1 FD 1 ρ
V1 d1
2
FD 2 6
= (48)2 = 20736 ⇒ FD 2 = 20736 · FD 1 = 311 lbf
FD 1 2
274
Problema 3
Solución
τ (d, g, ρp , ρ, µ)
Se aconseja sacar τ , µ (no es relevante) y ρp (lo piden en c))
ρp ρg 1/2 d3/2
τ
1. =f ,
ρgd ρ µ
τ
Ó = f (SG , Re)
ρgd
τ ρp
2. = f0
ρgd ρ
ρp ρp ρp 2650
3. = → ρp m = ρm = 800 = 2068 kg/m3
ρ m ρ p ρ p 1025
275
τ τ
4. = g es igual en el modelo y el prototipo.
ρgd
m ρgd
p
ρm 800
dp = dm = 150 µm = 117 µ m
ρp 1025
276
Problema 4
d(D, U, ρ, µ, σ)
Se aconseja sacar d, µ (no es relevante) y σ (dimensiones)
ρ U 2D ρ U 2D
2 σ
3. W em = W ep → = → Um = W ep
σ m σ p ρD m
s r
σm ρp D p 0.073 830
Um = Up = 20 5 = 58.8 m/s
σp ρm Dm 0.035 1000
277
d d Dp
4. = → dp = dm = 5 dm = 50 µ m
D m D p Dm
278
Problema 5
Solución:
V σ
2. µ no es relevante √ =f
gD ρgD2
s
V V VT (gD)T
3. √ = √ ⇒ =1= → (gD)T = (gD)L
gD T gD L VL (gD)L
DL gT
= =6
DT gL
279
Problema 6
Solución:
∆pρD2
∆s
= f α,
µ2 D
∆pρD2
2. ∆s no es relevante = f (α)
µ2
2 2
Lm 1 Dm ∆pρD ∆pρD
3. = = → =
Lp 5 Dp µ2 m µ2 p
2 2
Dm 1
∆pp = ∆pm = 1000 = 40 P a
Dp 5
280
Problema 7
Solución:
C k Cin V D
1. Grupos adimensionales πC = ; πk = ; πD = V 1/3
Cin Q Q
C k Cin V D 1/3
=g , V
Cin Q Q
C k Cin V
2. = g0
Cin Q
281
kp Cin p 633 3
k Cin V k Cin V L
m k L
m Cin mm
= → =
Q p Q m 18
Q m Q
m
Cin p 18
= 3 ≈ 0.083
Cin m 6
282
Problema 8
Solución
T = T (D, ρ, ω, µ) [T ] = M L2 t−2
Elección primeras variables más fácil {ρ, D, ω}
T µ
= f (3)
D5 ρω 2 D2 ρω
Prototipo: Petróleo ; Modelo: agua
Otras opciones
T µ
{ω, D, µ} 3
= f1
µDsω D2 ρω
3
ρ ω µ
{ρ, µ, ω} T = f2
µ5 D2 ρω
Tρ µ
{ρ, D, µ} = f3
µ2 D D2 ρω
Semejanza dinámica
2
µ µ µm Dp ρp
= → ωm = ωp = 0.48rad/s
D2 ρω p D2 ρω m µp Dm ρm
283
Tp ≈ 171 N · m
Potencia
Pp = Tp · ωp = 171 × 0.2 = 34.1 W
284
Problema 9
Solución
1.
[∆H] = L R=2 {D, ω}
[D] = L N =5 ∆H Q ν
=f ,
[ω] = t−1 M =N −R=3 D ωD3 ωD2
[Q] = L3 t−1
[ν] = L2 t−1
∆H 0 Q
2. =f
D ωD3
ωp Dp3
Q Q 1250
3. = → Qp = Qm = 50 = 35.7 m3 /s
ωD3 m ωD3 p ωmDm
3 1750
Pues Dp = Dm
285
Problema 10
Solución
1
ω = f (V, ρ, g, `, ρA ) `m = `p
10
1. V 2 ρA
ω`
{V, ρ, `} =f ,
V g` ρ`
V 2 ρA
r
`
{ρ, g, `} ω =f ,
g g` ρ`
V2 V
Nota: se puede expresar tambien como √
g` g`
ρa ρa `m 1
1. = → ρA m = ρA p = 6×10−3 ≈ 6 × 10−4 slug/f t2
ρ
` p ρ
` m `p 10
r
V V `m 1
2. √ = √ → Vm = Vp = 30 ≈ 9.49 f t/s
g ` m g ` p `p 10
s
√
r r
` ` `m
3. ω = ω → ωp = ωm = 6 10 ≈ 1.9 Hz
g m g p `p
286
Problema 11
Con el objeto de mantener un vuelo uniforme, los pájaros pequeños deben batir
sus alas más rápido que los grandes. Se ha sugerido que la relación entre la
frecuencia de batido, ω (batidos por segundo), y la envergadura del pájaro, L,
viene determinada por una relación funcional, ω ∼ Ln . Sabiendo que la fre-
cuencia de batido es función de la envergadura, el peso especı́fico del pájaro,
la aceleración de la gravedad y la densidad del aire, utilice el análisis dimen-
sional para determinar el exponente n. Sabiendo que un colibrı́ es 50 veces más
pequeño que un albatros y que éste bate sus alas una vez por segundo, halle
la frecuencia de batido del colibrı́ (asuma que ambos pájaros poseen el mismo
peso especı́fico promedio).
Solución
1
1. ω = f (`, ρ, g, γ) `col = `alb
50
287
Problema 12
Asuma que el caudal Q de una columna de gas que sale por una chimenea es
función de la densidad del aire ambiente, ρa , la densidad del gas en la columna,
ρg , la gravedad, g, y el diámetro de la columna, d, respectivamente. Se ha
sugerido que dicho caudal depende del diámetro de la columna siguiendo una
relación funcional Q ∼ dn . Utilice el análisis funcional para determinar el ex-
ponente n.
Solución
Q = Q (ρa , ρg , g, d) Q ∼ dn
Lp = 5 Lm para satisfacer la semejanza dinámica todos los π’s deben ser iguales.
Luego ρg m = ρg p ⇒ se debe usar CO2 .
Qm = 4 m3 /h
Semejanza Dinámica
! !
Q Q
=
g 1/2
d5/2
m d5/2g 1/2
p
5/2
dp
Qp = Qm = 4(5)5/2 = 223.6 m3 /h
dm
288
Problema 13
En ciertas circunstancias, el viento que sopla sobre una señal de tráfico rectan-
gular puede hacerla oscilar con una frecuencia ω. Asuma que ω es función de
la anchura y altura de la señal, b y h, la velocidad del viento, V , la densidad
del aire ρ y la constante elástica k del poste soporte. Encuentre una expresión
adimensional para este problema.
Solución
ω = ω (b, h, V, ρ, k)
k k
Semejanza dinámica =
b3 V 2 ρ
m b3 V 2 ρ
p
2 3 3/2 1/2
Vm bp km bp km
= → Vm = Vp
Vp bm kp bm kp
1/2
√
3/2 1
Vm = 30(10) = 30 10 = 94.87 km/h
100
ωb ωb
Por otro lado =
V p V m
289
Vp bm √ −1 1
ωp = ωm = 40 10
Vm bp 10
ωp = 1.26 Hz
290
Problema 14
Solución
Análisis dimensional:
r
ρ d u2
b µ ρd b µ
=g , √ ó u=g , √
σ d ρdσ σ d ρdσ
2
dp 1.2 × 10−3 998 0.0728
= 2 = 5.8
dm (1 × 10−3 ) 789 0.0228
2 2
ρ d u2 2 ρp 5.8
d
m up
ρd
u up σp ρm 1
= = ⇒ =
σ p σ m σp um σm ρp 5.8
291
r
up 1 998 0.0228
= = 0.26
um 5.8 789 0.0728
292
Problema 15
Se construye un modelo a escala 1:5 (el modelo es cinco veces más pequeño que
el prototipo) de la válvula que va a funcionar con agua. Para una apertura
particular cuando la velocidad de la parafina es de 3.0 m/s, ¿Cuál debe ser
la velocidad del agua en el modelo para trabajar en condiciones de semejanza
dinámica? ¿Cuál es el cociente de caudales entre prototipo y modelo? Encuentre
la caı́da de presión en el prototipo si en el modelo se miden 60 kP a.
La densidad y viscosidad de la parafina son 800 kg/m3 y 0.002 kg/(m s) respec-
tivamente; la viscosidad cinemática del agua es 10−6 m2 /s.
Solución
1. ∆p = ∆p (h, d, V, ρ, µ)
∆p h ρV d
=g ,
ρV 2 d µ
2. Lp = 5 Lm
10−6 · 800
ρV d ρV d Vm µm ρp dp
= → = = 5=2
µ p µ m Vp µp ρm dm 0.002
µm
= 10−6 : es la viscosidad cinemática del agua.
ρm
Vm = Vp · 2 = 2 · 3 = 6 m/s
2
Vp π4 d2p
Qp Vp dp 1
= π 2 = = × 52 = 12.5
Qm Vm 4 dm Vm dm 2
2
∆p
∆p
ρV p
= ⇒ ∆pp = ∆pm
ρV 2 p ρV 2 m (ρ V 2 )m
293
2
800 1
∆pp = 60 · = 12 kP a
1000 2
294
Problema 16
Solución
295
Problemas Propuestos
Problema 1
Problema 2
Sobre una presa fluye una lámina de lı́quido, cuyo caudal por unidad de longitud
q a lo largo de la presa es función de la altura de la lámina de lı́quido sobre la
presa H, su anchura b, la aceleración de la gravedad g, la densidad del lı́quido
ρ y su viscosidad µ.
296
3. En las condiciones del punto anterior, se realiza un experimento en un
modelo a escala 1/20. Asumiendo semejanza dinámica cuál será el caudal
por unidad de longitud en el prototipo si el caudal por unidad de longitud
en el modelo fue de 1 m2 /s?
Suponga que el lı́quido es agua tanto en el modelo como en el prototipo
Problema 3
Problema 4
297
Se desea realizar un experimento basado en un modelo, con una escala de lon-
gitud Lm /Lp = 1/4 y una escala de densidades ρm /ρp = 1, para determinar
el espesor de la capa de fluido que fluye sobre el obstáculo. Determine las es-
calas de viscosidad y tensión superficial que se deben satisfacer entre modelo y
prototipo para que se verifique semejanza dinámica.
Problema 5
Se desea determinar la máxima anchura de una lámina de oro que puede flotar en
agua mediante un modelo a escala que considera una lámina de plomo flotando
en aceite. Bajo condiciones de semejanza dinámica determine:
Datos: σagua = 0.073 N/m, σaceite = 0.033 N/m, ρoro = 19300 kg/m3 ,
ρplomo = 11340 kg/m3
Problema 6
298
la expresión adimensional que gobierna este problema. Un golfista profesional
golpea la bola, de diámetro D = 4.27 cm, con una velocidad V = 75 m/s y
ω = 8100 rpm (vueltas por minuto). El modelo se desea probar en un túnel de
viento con una velocidad de 25 m/s; bajo condiciones de semejanza dinámica
determine el diámetro de la bola en el modelo y su velocidad de rotación.
Nota: Los problemas resueltos en clase le dan una indicación para la elección
de las primeras variables. Las dimensiones de la velocidad de rotación son t−1 .
Problema 7
Nota: Los problemas resueltos en clase le dan una indicación para la elección
de las primeras variables.
Problema 8
299
Problema 9
Datos: ρagua = 1000 kg/m3 , µagua = 10−3 kg/(m s), σagua = 0.073 kg/s2 ,
ρoil = 800 kg/m3 , µoil = 0.03 kg/(m s), σoil = 0.033 kg/s2
Problema 10
Datos: ρagua = 1000 kg/m3 , µagua = 10−3 kg/(m s), ρoil = 800 kg/m3 , µoil =
0.03 kg/(m s), σoil = 0.033 kg/s2 , Dagua = 5×10−5 m2 /s, Doil = 1.5×10−5 m2 /s,
[h] = L t−1
300
Capı́tulo 7
301
1 2 1 1
p∞ + ρV∞ = ps + ρVs2 ⇒ ps = p∞ + ρ(V∞ 2
− Vs2 )
2 2 2
1 2 1 2
ps = p∞ + ρV∞ (1 − 4 sin2 θ) ⇒ ∆p = ps − p∞ = ρV∞ (1 − 4 sin2 θ)
2 2
Como la distribución de presiones es simétrica respecto del eje vertical, la fuerza
de presión neta en la dirección horizontal es nula, luego la fuerza de resistencia
se anula. En este caso la distribución de presiones también es simétrica respecto
del eje horizontal por lo que la fuerza de sustantación también es cero.
302
posee un valor apreciable, pues aunque la viscosidad µ sea muy pequeña el gradi-
ente de velocidad es muy alto. Como consecuencia la resistencia τ ×S es grande.
Fue Ludwig Pradtl quien desarrolló la teorı́a de la Capa Lı́mite a principios del
siglo XX y se le considera el padre de la Mecánica de Fluidos moderna.
La Capa Lı́mite:
303
• Fuera de la capa lı́mite el fluido se comporta como ideal y los efectos
viscosos son despreciables.
• La presión en el interior de la capa lı́mite se mantiene muy aproximada-
mente constante.
Por convenio se define el espesor de la capa lı́mite δ como el lugar de los pun-
tos donde la velocidad u paralela a la superficie alcanza el 99% del valor de la
velocidad externa U .
En el caso de una placa plana dicho espesor depende del régimen de flujo, lami-
nar o turbulento:
δ 5.0
Laminar: ≈ 1/2
x Rex
δ 0.16
Turbulento: ≈ 1/7
para Rex > 106
x Rex
ρU x
donde x es la distancia desde el borde de ataque de la placa y Rex = es el
µ
número de Reynolds local.
Ejemplo 7.1:
Solución:
δ 5.0 δ 2 U/ν
≈ ⇒ x = = 511 f t
x 1/2
Rex 5.02
Ux
con Rex = = 9.4 × 108 > 106 lo cual es incorrecto.
ν
Debemos usar la expresión para flujo turbulento
7/6
δ(U/ν)1/7
δ 0.16
≈ 1/7
⇒x= = 5.17 f t
x Rex 0.16
Ux
y el test Rex = = 9.5 × 108 que es correcto.
ν
304
Capı́tulo 8
Ejemplo 8.1:
Solución:
Casi todas las fórmulas para tuberı́as están basadas en la velocidad promedio
V = Q/A . La transición tiene lugar a ρ V d/µ ≈ 2300
que son velocidades muy bajas, luego la mayorı́a de flujos con aire y agua
en la ingenierı́a son turbulentos, los flujos laminares tienen lugar con fluidos
más viscosos como aceites de lubricación o glicerina.
305
Descripción del flujo en una tuberı́a circular
Un flujo interno esta constreñido por paredes. Los efectos viscosos en las pare-
des crecerán y se encontraran permeando el flujo entero.
306
• Zona de flujo totalmente desarrollado
Le
Flujo laminar: ≈ 0.06Red
d
para Red crit ≈ 2300 ⇒ Le = 138d es la máxima longitud de desarrollo posible
Le 1/6
Flujo turbulento: ≈ 4.4Red
d
la capa limite crece más rápido y es relativamente corta, por ejemplo:
Red = 106 → Le ≈ 44d
Ejemplo 8.2:
Solución:
0.0023 f t3 s−1
Q = 5 gal/min = 5 × = 0.0111 f t3 s−1
1 gal/min
Q 0.0111 f t3 s−1
V = = 2 = 8.17 f t/s
A π 1
0.5 f t
4 12
307
a 20◦ C νw = 1.01 × 10−6 m2 s−1 = 1.09 × 10−5 f t2 /s. El Red , es
308
2. Flujo imcompresible.
3. No hay máquinas.
4. No hay transferencia de calor
Aplicamos entonces el análisis integral del VC al flujo de la figura entre las sec-
ciones 1 y 2:
Ecuación de continuidad
Q1 = Q2 = cte
Q1 Q2
V1 = = V2 =
A1 A2
Ecuación de la energı́a (flujo estacionario)
p1 1 p2 1
+ V12 + g z1 = + V22 + g z2 + hf
ρ 2 ρ 2
Como V1 = V2 la altura de pérdidas por fricción se escribe:
p1 p2 p ∆p
hf = z1 + − z2 + = ∆ z1 + = ∆z +
ρg ρg ρg ρg
Ecuación de momento
Las fuerzas que actúan sobre el fluido son la fuerza de presión, el peso y la
fricción con las paredes (debido al esfuerzo cortante en la pared τw ). Observe
que ∆z = ∆L sin φ.
∆p 2τw ∆L
∆z + = hf = (1)
ρg ρg R
Ahora podemos intentar relacionar τw con las condiciones de flujo (laminar o
turbulento). Para ello aplicamos seguidamente el Análisis Dimensional.
τw = F (ρ, V, µ, d, ε)
309
LV2
hf = f Ecuación de Darcy-Weisbach
d 2g
Ecuación válida para flujos en conductos de cualquier sección y para flujo lami-
nar o turbulento.
ε
Finalmente nos resta encontrar la forma de la función F Red , en función
d
de los regı́menes de flujo.
310
64µ L V 2 32µLV
hf lam = = y hf lam ∝ V
ρV d d 2g ρ g d2
Ejemplo 8.3:
Aceite (ρ = 900 kg m−3 , ν = 2 × 10−4 m2 s−1 ) fluye hacia arriba en una tuberı́a
circular inclinada 40◦ respecto de la horizontal. Su diámetro es 6 cm, su longi-
tud 10 m y la diferencia de presiones entre la entrada y la salida es de 105 P a.
Asumiendo flujo laminar hallar (a) hf entre 1 y 2 (b) Q (c) V (d) Red . ¿Es
el flujo realmente laminar?
Solución:
µ = ρν = 0.18 kg m−1 s−1
z2 = ∆L sin 40◦ = 6.43m
(a) Aplicamos la ecuación de la energı́a asumiendo flujo totalmente desarrollado
p1 V12 p2 V22
z1 +
ρg + = z2 + + + hf
2g ρg 2g
311
Perfiles de velocidad instantáneo y promedio en un flujo turbulento
ε
Ya vimos que f = F Red , donde ε era la rugosidad absoluta de la tuberı́a.
d
En flujo turbulento el factor de fricción se estima mediiante experimentos cuyos
resultados se ajustan a correlaciones.
ε/d
Tuberı́as muy rugosas f −1/2 = −2.0 log (independiente del Red )
3.7
Régimen de transición liso-rugoso
−1/2 ε/d 2.51
f = −2.0 log + Colebrook (1939)
3.7 Red f 1/2
Todas estas correlaciones incluyendo la formula laminar, se recogen en el dia-
grama de Moody (1944). Para 2000 < Red < 4000 (transición flujo laminar-
turbulento) no hay datos fiables para f y por eso la zona correspondiente está
sombreada en el diagrama. Notar que para los Red más elevados las curvas son
casi horizontales.
312
Ejemplo 8.4:
Solución:
∆p
= hf ⇒ ∆p = ρ g hf = 62.4 lbf /f t3 × 4.5 f t = 280 lbf /f t2
ρg
313
8.5 Cálculos tipo. Determinación pérdida de carga, caudal y diámetro
El diagrama de Moody puede usarse para resolver casi todos los problemas que
implican pérdidas por fricción en tuberı́as largas. Pero muchos problemas impli-
can iteración y cálculos repetidos ya que el diagrama de Moody es un diagrama
de altura de pérdidas.
Problemas tipo:
1. Dados d, L, V o Q, ρ, µ, g hallar hf (problema altura de pérdidas).
⇒ ya visto
2. Dados d, L, hf , ρ, µ, g hallar V ó Q (problema de caudal).
3. Dado Q, L, hf , ρ, µ, g hallar d (problema de tamaño).
Ejemplo 8.5:
Solución:
1. Solución directa
314
1.1 × 105
Red
V =ν = (274800) ≈ 6.05 f t/s
d 0.5
2. Solución iterativa
Solución:
315
ε
f ≈ 0.03 → d ≈ 0.542 f t → Red ≈ 252000 → ≈ 7.38 × 10−4
d ε
fnew ≈ 0.0196 → dnew ≈ 0.498 f t → Red ≈ 274000 → ≈ 8.03 × 10−4
d
fbetter ≈ 0.0190 → d ≈ 0.499 f t
πd2 /4 d R
Rh = = = ⇒ d = 4Rh
πd 4 2
En el caso de otras secciones el radio hidráulico se calcula de la misma manera.
316
Por consiguiente, en un conducto no circular la Ecuación de Darcy-Weisbach se
expresa en función de Rh como
L V2
hf = f
4Rh 2g
la cual es aplicable a canales y ductos de sección transversal cualquiera. Del
mismo modo el diagrama de Moody puede aplicarse en este caso sustituyendo
d por 4Rh .
Las pérdidas locales o secundarias son pérdidas de forma que tienen lugar en
los cambios de sección y dirección de la corriente y en todos los accesorios de
tuberı́as. Estos elementos producen una perturbación de la corriente que orig-
ina remolinos y desprendimientos, que intensifican las pérdidas. Sin embargo,
estas pérdidas secundarias pueden ser más importantes que las primarias si la
conducción es relativamente corta. Se admite que si Ltotal /d > 1000 el error
cometido al despreciar las pérdidas secundarias es menor que el cometido al
calcular el valor del factor de fricción f usando el diagrama de Moody.
Un ejemplo: las pérdidas en una válvula pueden ser muy pequeñas cuando está
totalmente abierta, y muy importantes cuando está parcialmente abierta.
Le V 2
hf = f
d 2g
1er método:
V2
hs = ζ pérdida de carga secundaria
2g
317
Tabla de coeficientes de pérdidas por accesorios
318
Coeficiente total de pérdidas
Ejemplo 8.7:
Solución:
p1 = p2 = pa V1 = V2 ≈ 0
0 = z2 − z1 + hf t ⇒ z2 − z1 = H
319
V2 LV2
H = (ζentrada + ζvalvula + 2ζcodo + ζsalida ) + +f
2g d 2g
Q 0.04
donde V = 2
= = 5.09 m s−1
π d /4 π (0.1)2 /4
Vd 5.09 × 0.1
Red = = = 5.09 × 105
νH2 O 10−6
ε 0.046 × 10−3 M oody
= = 4.6 × 10−4 −−−−→ f = 0.0173
d 0.1
Aplicando los coeficientes de pérdidas de la tabla (10 cm = 3.94 in ≈ 4 in)
5.092 50 5.092
H = (0.5+5.7+2×0.64+1.0) +0.0173 = 11.2+11.4 = 22.6 m
2 × 9.8 0.1 2 × 9.8
◦
2 método (longitud equivalente)
Ejemplo 8.8:
Solución:
Tenemos
V2 LV2 L
hs = ζ =f ⇒ ζ=f
2g d 2g d
Hallamos primero f
Vd 4 × 0.1 ε 0.17 × 10−3
Red = = = 400000 = = 0.17 × 10−2
ν 10−6 d 0.1
Diagrama de Moody ⇒ f ≈ 0.0229
8
ζ = 0.0229 = 1.832
0.1
320
Existen tres configuraciones básicas que permiten construir las redes de dis-
tribución de fluidos.
321
Q1 = Q2 = Q3 = constante ⇒ V1 d21 = V2 d22 = V3 d23
N1 N2 N3
V12 L1 X V 2
L 2
X V 2
L 3
X
∆hA→B = f1 + ζ1j + 2 f2 + ζ2j + 3 f3 + ζ3j
2g d1 j=1
2g d2 j=1
2g d3 j=1
V12
∆hA→B = (α0 + α1 f1 + α2 f2 + α3 f3 )
2g
donde los coeficientes αi dependen solo de la geometrı́a y los accesorios.
322
1. Si ∆h es conocido, en cada tramo se halla Qi iterando desde la estimación
para fi de tuberı́as totalmente rugosas (como en el caso anterior).
P
2. Si Q = Qi es conocido, se debe determinar ∆h y como se distribuye el
caudal en cada tramo. De nuevo el cálculo es iterativo. En cada tramo
L V2 Q2 π 2 gd5
hf = f =f ; C=
4d 2g C 8L
La pérdida de carga se relaciona con el caudal del siguiente modo:
Q2 π 2 gd5i
hf = P p 2 ; Ci =
8Li
Ci /fi
Li Vi2
∆hi = fi = zi − HJ , i = 1, 2, 3
di 2g
El problema consiste en encontrar el caudal total que circula por el sistema y
como se distribuye. Este caso es más complejo que los anteriores y normalmente
323
los valores de fi se toman como los de flujo totalmente rugoso. En la iteración
se supone un valor de HJ y de la expresión anterior se calculan Vi , es decir, Qi .
Se verifica la suma de los caudales frente a cero. Si el valor de HJ es muy alto,
P
Qi será negativa, por lo que habrá que disminuirlo, y viceversa.
Red general
324
Problemas Resueltos
Problema 1
Nota: Para el aire R = 1716 f t lbf /(slug ◦ R) , µ = 4.7 × 10−7 lbf s/f t2 , y 0
◦
F = 460 ◦ R. El tubo de escape descarga a la atmósfera
Solución:
p2 (14.7 P si) × 144
1. p2 = ρRT → ρ= = = 1.74 × 10−3 slug/f t3
RT 1716 · (460 + 250)
LV2 1.5 Q 0, 2
2. hf p = f d= = 1/8 f t V = π = 2 = 16.3 f t/s
d 2g 12 d 2 π 1
4
4 8
ε 0.00085 ρV d 1.74 × 10−3 · 16.3 · 1/8
= = 0.0068 ; Red = = = 7543
d 1/8 µ 4.7 × 10−7
14 (16.3)2
Moody f ≈ 0.042 → hf p = 0.042 ≈ 19, 41 f t
1/8 2 · 32.2
V2
P
hs = ζi Elementos: 6 codos roscados de 1.5” ζcodos = 1, 225 (interpolado)
3. 2g i
Silenciador ζsil = 8.5 ; Salida ζ = 1.0
(16.3)2
hs = (8.5 + 1.0 + 6 × 1.225) = 69.52 f t
2 · 32.2
325
p1 V2 p2 V2
4. Ecuación de energı́a + 1 + z1 = + 2 + z2 + (hf p + hs )
ρg 2g ρg 2g
V1 = V2 ; z1 = z2
p1 − p2 = ρg(hf p + hs ) = 1.74 × 10−3 × 32.2(19.41 + 69.52)
326
Problema 2
Se descarga agua a 15◦ C (ρ = 999 kg/m3 , µ = 1.14 × 10−3 kg/(m · s)) desde
un depósito a una tasa de 18 litros/s utilizando dos tuberı́as de hierro forjado
(ε = 0.26 mm) conectadas por una bomba. La primera tuberı́a tiene una lon-
gitud de 20 m y un diámetro de 6 cm y la segunda una longitud de 35 m un
diámetro de 4 cm. El nivel del agua está situado a 30 m sobre el centro de la
tuberı́a, la cual tiene una entrada con aristas vivas. Las pérdidas asociadas a la
conexión de la bomba son despreciables. Determine la altura de la bomba y la
potencia mı́nima de la bomba para mantener el flujo volumétrico indicado.
Solución
µ
Datos: H2 O → ρ1 µ ⇒ ν = , Q Material Hierro forjado ε (mm)
ρ
Entrada tuberı́a arista viva ⇒ ζe = 0.5
p1 V12 p2 V22
+ + z1 + hb = + + z2 + hf hf : perdidas por fricción y accesorios
ρg 2g ρg 2g
Como p1 = p2 = patm y el depósito es grande V1 ≈ 0
V22 V22
Luego z1 + hb = + z2 + hf ⇒ hb = −(z1 − z2 ) + + hf
2g 2g
327
Q
V2 = π Donde d2 está en cm y se debe pasa a m
d22
4
Las pérdidas hf son las correspondientes a los tramos 1 y 2
hf = hf 1 + hf 2
2
L1 V1 Q
hf 1 = f1 + ζe Con V1 = π
d1 2g d2
4 1
ε ρV1 d1 ε
f1 Red1 , Con Red1 = , con ε en mm debe pasarse a cm
d1 µ1 d1
El factor de fricción se obtiene del diagrama de Moody o de
" 1.11 #
1 6.9 ε/d
√ ≈ −1.8 log +
f Red 3.7
L2 V22
En el tramo 2 hf 2 = f2 pues no hay perdidas locales
d2 2g
ε ρV2 d2 ε
f2 Red2 , Red2 = , >> f2
d2 µ2 d2
V22
hb = −z1 + + (hf 1 + hf 2 ) ⇒ Ẇb min = ρgQhb
2g
Es la potencia mı́nima porque considera eficiencia del 100%
328
Problema 3
• Halle el flujo volumétrico de agua (en litros/s) que debe utilizarse para
reducir la salinidad del agua del estanque a 500 mg/L en un año.
Para sacar el agua salada del estanque se utiliza una cierta bomba a través de
una conducción de hierro forjado (ε = 0.046 mm). El primer tramo, desde el
estanque a la bomba tiene un diámetro d1 y longitud L1 mientras que el segundo,
desde la bomba al océano, tiene un diámetro d2 y longitud L2 . El agua salada
se toma del estanque mediante una válvula de pie (ζ = 1.8) y se conduce hasta
la bomba utilizando 4 codos estándar de 90◦ y una válvula de globo (totalmente
abierta). Desde la bomba hasta el océano se encuentra una válvula de retención
de disco oscilante (abierta). Todos los accesorios se encuentran roscados.
Solución:
329
20000 0.5
Q=− ln ≈ 214.4 m3 /dia = 2.48 L/s
365 dias 25
zA = 20 m; zB = 0 m; ε = 0.046 mm; pA = pB = pa t
VB2 V2
zA + hb = + hT → hb = −zA + B + hT (2)
2g 2g
hT = h1 + h2
4Q V1 d1
V1 = = 0.56 m/s ; Re1 = = 42000 ; ε/d1 = 6.1×10−4 → f1 = 0.0233
π d21 ν
4Q V2 d2
V2 = = 1.26 m/s ; Re2 = = 63000 ; ε/d2 = 9.2×10−4 → f2 = 0.0228
π d21 ν
(Haaland)
!
V12 0.562
L1 X (1) 30
h1 = f1 + ζi = 0.0233 · + (1.8 + 4 × 0.795 + 6.3) = 0.33 m
d1 i
2g 0.075 2 × 9.8
!
V22 1.262
L2 X (2) 1200
h2 = f2 + ζi = 0.0228 · + 2.1 = 44.49 m
d2 i
2g 0.05 2 × 9.8
1.262
De (2) hb = −20 + + 44.49 + 0.33 = 24.9 m
2 × 9.8
330
Problema 4
Solución
1. Conservación [NH3 ] C = [N H3 ]
Cd Qd 3 × 50
Cd Qd +
C
in Qin = CT QT → QT = = = 300 m3 /day
CT 0.5
1
QT = 300 m3 /day ≈ 3.5 × 10−3 m3 /s = 3.5 L/s
86400
2. Esquema sistema de bombeo
3. Aplicando ecuación de energı́a en 1 - 2 z2 = 0 m
p1 V12 p2 V22 V22
+ +
z
1 = + +2 + hT + hw ⇒ hw = z1 − hT −
z
ρg 2g ρg 2g 2g
331
V22
LT P
hT = hf + hs = f + ζi
2g d i
Q 4Q 4 × 3.5 × 10−3
V2 = π = 2
= = 0.45 m/s
d2 πd π 0.12
4
Coeficientes de pérdidas de los accesorios
Tramo Accesorio # Vr. Unit ζ
Entrada admisión 1 0.8 0.8
entrante
Válvula globo 1 5.7 5.7
abierta
Codo ext. larga 4 0.23 0.92
Válvula compuerta 1 0.11 0.11
abierta P
ζi = 7.53
i
1
Haaland f= " 1.11 #!2 ≈ 0.0216
6.9 ε/d
−1.8 log10 +
Red 3.7
Vd 0.45 × 0.1
Red = = ≈ 45000
ν 10−6
0.015
ε/d = = 1.5 × 10−4
100
0.452
50
hT = 0.0216 × + 7.53 = 0.19 m
2 × 9.8 0.1
0.452
hw = 5 − 0.19 − = 4.8 m
2 × 9.8
332
Problema 5
Solución
kCout V 5 × 2 × 104
Q= = = 0.037 m3 /s = 37 L/s
Cin − Cout 10 − 5
2. Esquema sistema de bombeo
3. Aplicando ecuación de energı́a en 1 - 2 25 psi = 1.7 × 105 P a
p1 V12 p2 V 2
+ + z1 + hb = + 2 + z2 + hT ⇒
ρg 2g ρg 2g
p2 1.7 × 105
hb = z2 + + hT = 50 + + hT ≈ 67.35 + hT
ρg 1000 × 9.8
333
V2
LT P
hT = 2 f + ζi
2g d i
Q 4Q 4 × 37 × 10−3
V = π = 2
= = 18.8 m/s
d2 πd π 0.052
4
Coeficientes de pérdidas de los accesorios
Tramo Accesorio # Vr. Unit ζ
Entrada borde 1 0.5 0.5
cuadrado
Válvula globo 1/2 1 17 17.0
abierta
Codo 90◦ estándar 4 0.95 3.8
Codo 45◦ 2 0.3 0.6
Válvula disco os- 1 2.1 2.1
cilante
salida 1 P1 1.0
ζi = 25.0
i
1
Haaland f= " 1.11 #!2 ≈ 0.024
6.9 ε/d
−1.8 log10 +
Red 3.7
Vd 18.8 × 0.05 0.1
Red = = ≈ 940000, ε/d = = 2 × 10−3
ν 10−6 50
18.82
40
hT = 0.024 × + 25.0 = 797 m
2 × 9.8 0.05
334
Problema 6
Datos: ρagua = 1.94 slug/f t3 , ν = 1.217 × 10−5 f t2 /s, g = 32.2 f t/s2 . Cono
reductor de 8” a 2” ζ = 8
No considere la salida a la atmósfera como accesorio. Tome 3 cifras significativas
en el cálculo de los coeficientes de fricción.
Solución
2 2 2 !2 2
d200
A1 1 1
ζensanch. = 1− = 1 − 24 = 1− = 1− = 0.56
A2 d800 2 4
L400 = 260 f t ; L800 = 30 f t ; z1 − z2 = H = 100 f t.
p1 V12 p2 V2
C.E. entre 1 y 2 + + z1 = + 2 + z2 + hT
ρg 2g ρg 2g
V22
z1 − z2 = H = + hT
2g
V : velocidad tuberı́a 4”
2
4
C.M. V (d400 )2 = V2 (d200 )2 → V2 = V = 4V
2
2
4
V800 (d800 )2 = V (d400 )2 → V800 = V = 0.25V
8
335
V2 V8200
L X L X
hT = f + ζ + f + ζ
d 400 2g d 800 2g
V2 V2
hT = (f400 × 780 + 1.82) + (f800 × 45 + 8) · 0.252
2g 2g
V22
Ecuación de la energı́a H= + hT
2g
V2 V2 V2
100 = 16 + (f400 × 780 + 1.82) + (f800 × 45 + 8) · 0.252
2g 2g 2g
2g · 100 = V 2 16 + f400 × 780 + 1.82 + 0.252 × f800 × 45 + 8 × 0.252
336
(1) (1)
Haaland → f400 = 0.0175 f800 = 0.0171
2
π 2 π 4
Q=V d400 = 14.18 · ≈ 1.24 f t3 /s
4 4 12
337
Problema 7
Solución:
338
2. ε = 0.015 mm d1 = 10 cm = 4 in L1 = 200 m Accesorios roscados,
d2 = 7.5 cm = 3 in L2 = 30 m
Esquema instalación
f1 ≈ 0.0151 f2 ≈ 0.0152
339
Coeficientes de pérdidas de los accesorios
Tramo 1 Accesorio # Vr. Unit ζ
Válvula 1 1.3 1.3
Codo estándar 20 0.64 12.8
T Ramificada 1 1.1 1.1
T Lineal 3 0.9 2.7
Válvula globo 1/2 1 14 14
abierta
(1)
ζT = 31.9
42 7.22
200 30
hT = 0.0151 × + 31.9 + 0.0152 × + 4.315
0.1 2 × 9.8 0.075 2 × 9.8
hw = z2 + hT = 20 + 78.2 = 98.2 m
Ẇideal = γQhw = 9800 × 31.7 × 10−3 × 98.2 = 30507 W ≈ 30.5 kW
Ẇideal 30.5
Ẇreal = = ≈ 40.7 kW
η 0.75
340
Problema 8
Solución:
341
V22
Ecuación de energı́a entre 1) y 2) z1 + hb = + hT
2g
(z2 = 0)
!
V22 25.462
L X 620
hT = hf +hs = f + ζi = 0.0156 + 0.5 + 0.28
D i
2g 1/2 2 × 32.2
hT = 202.56 f t
342
Problema 9
• Halle el caudal necesario para que el agua a la salida del tanque alcance la
concentración de 5 mg/L en una semana. Exprese su resultado en L/day.
El agua se extrae del tanque utilizando una cierta bomba a través de una con-
ducción de fundición corriente nueva (ε = 0.25 mm), llevándose a una laguna
cercana. El primer tramo, desde el tanque a la bomba tiene un diámetro d1
y longitud L1 mientras que el segundo, desde la bomba a la laguna, tiene un
diámetro d2 y longitud L2 . El agua del tanque se conduce hasta la bomba
utilizando 6 codos estándar de 90◦ , 1 Te ramificada y una válvula de ángulo
(abierta). Desde la bomba hasta la laguna se encuentra una válvula de com-
puerta (abierta), una curva de retorno y dos codos de 45◦ . La tuberı́a descarga
a la presión atmosférica, unos pocos centı́metros sobre el nivel de la laguna.
Todos los accesorios se encuentran roscados.
Solución:
343
Esquema instalación
f1 ≈ 0.0271 f2 ≈ 0.0305
344
L1 P (1) V12 L2 P (2) V22
hT = f1 + i ζi + f2 + i ζi
d1 2g d2 2g
11.22 25.212
25 300
hT = 0.0271 × + 8.02 + 0.0305 × + 1.71
0.075 2 × 9.8 0.05 2 × 9.8
25.212
hw = z2 + hT = 5 + + 6098, 5 ≈ 6136 m
2 × 9.8
Ẇideal = γQhw = 9800 × 49.5 × 10−3 × 6136 = 2976574 W ≈ 2977 kW
Ẇideal 2977
Ẇreal = = ≈ 3721 kW
η 0.8
345
Problema 10
• Halle el caudal necesario para que el agua a la salida del reactor cumpla
la norma. Exprese su resultado en L/s.
El agua tratada efluente del reactor se utiliza para recuperar energı́a a través de
una turbina hidráulica y se conduce hasta un rı́o cercano. El agua se conduce a
través de una conducción de fundición asfaltada (ε = 0.1 mm). El primer tramo,
desde el reactor a la turbina tiene un diámetro d1 y longitud L1 mientras que
el segundo, desde la turbina al rı́o, tiene un diámetro d2 y longitud L2 . El agua
del reactor se conduce hasta la turbina utilizando 4 codos de 45◦ , una curva
de retorno y una válvula de globo (abierta un cuarto). Desde la turbina hasta
el rio se encuentra una válvula de compuerta (abierta), dos codos de extensión
larga y 8 codos estándar de 90◦ . La tuberı́a descarga a la presión atmosférica,
unos pocos centı́metros sobre el nivel del rio. Todos los accesorios se encuentran
roscados.
• Teniendo en cuenta que el reactor está situado 80 m por encima del rio,
encuentre la potencia extraı́da por la turbina en kW del caudal que usted
calculó en el punto a), asumiendo que la máquina tiene una eficiencia del
90%.
Solución:
346
2. ε = 0.1 mm d1 = 10 cm = 4 in L1 = 15 m Accesorios roscados
d2 = 7.5 cm = 3 in L2 = 150 m
Esquema instalación
f1 ≈ 0.0208 f2 ≈ 0.0218
347
Coeficientes de pérdida de los accesorios
Tramo 1 Accesorio # Vr. Unit ζ
Entrada 1 0.5 0.5
Codo 45◦ 4 0.29 1.16
Curva retorno 1 0.64 0.64
Válvula globo 1/4 1 40 40
abierta
(1)
ζT = 42.30
2.32 42
15 150
hT = 0.0208 × + 42.3 + 0.0218 × + 7.135
0.1 2 × 9.8 0.075 2 × 9.8
V22 42
hw = z1 − − hT = 80 − − 53.68 ≈ 25.5 m
2g 2 × 9.8
Ẇideal = γQhw = 9800 × 17.75 × 10−3 × 25.5 = 4436 W ≈ 4.4 kW
348
Problema 11
Solución
Cin
1. PFR segundo orden reacción C(L) =
1 + 2k2 Cin V /Q
349
m3 L
Pues k2 = 2 × 10−3 = 2 × 10−6 V = 180 m3 = 180000 L
kg · s mg · s
2. Esquema instalación
f1 ≈ 0.0260 f2 ≈ 0.0247
350
Coeficientes de pérdida de los accesorios
Tramo 1 2”
Accesorio # Vr. Unit ζ
Entrada 1 0.5 0.5
Codo estándar 4 0.95 3.8
Curva retorno 1 0.95 0.95
Válvula globo 1/2 1 17 17
abierta
(1)
ζT = 22.25
Tramo 2 3”
Accesorio # Vr. Unit ζ
2+1
Válvula ángulo 1 1.5
2
0.9 + 0.9
T Lı́nea 1 0.9
2
0.41 + 0.23
Codo ext. larga 3 0.96
2
Salida 1 1 1
(2)
ζT = 4.36
1.962
(1) 10
hT = 0.026 × + 22.25 × = 5.38 m
0.05 2 × 9.8
0.862
(2) 20
hT = 0.0247 × + 4.36 × = 0.41 m
0.075 2 × 9.8
Ẇideal 586.5
Ẇreal = = ≈ 782 W = 0.782 kW
η 0.75
351
Problema 12
• Teniendo en cuenta que el reactor está situado 8 m por encima del es-
tanque, encuentre la potencia real proporcionada por la turbina en kW
asumiendo que la máquina tiene una eficiencia del 91%.
Solución
Q 1 C
1. C(t) = C0 e−( V +k1 )t −→ Q = V − ln − k1
t C0
352
1 5 1
Q = 20 ln − 0.5 = 229.6 m3 /day× = 2.66 × 10−3 m3 /s
6/24 100 86400
2. Esquema instalación
p1 V12 p2 V 2 V22
+ + z1 = + 2 + z2 + hw + hT ⇒ hw = (z1 − z2 ) − − hT
ρg 2g ρg 2g 2g
0.1 0.1
ε/d1 = = 1.3 × 10−3 ε/d2 = = 10−3
75 100
1
Haaland f= " 1.11 #!2
6.9 ε/d
−1.8 log10 +
Red 3.7
f1 ≈ 0.02495 f2 ≈ 0.0251
353
Coeficientes de pérdida de los accesorios
Tramo 1 3”
Accesorio # Vr. Unit ζ
Entrada 1 0.5 0.5
0.95 + 0.64
Curva retorno 2 1.59
2
T lı́nea 3 0.9 2.7
2+1
Válvula ángulo 1 1.5
2
(1)
ζT = 6.29
Tramo 2 4”
Accesorio # Vr. Unit ζ
Válvula globo(ab) 1 5.7 5.7
Codo estándar 2 0.64 1.28
Codo ext. larga 3 0.23 0.69
(2)
ζT = 7.67
0.62
(1) 100
hT = 0.02495 × + 6.29 × = 0.73 m
0.075 2 × 9.8
0.342
(2) 10
hT = 0.0251 × + 7.67 × = 0.06 m
0.1 2 × 9.8
0.342
hw = 8 − − 0.73 − 0.06 = 7.2 m >0 OK
2 · 9.8
354
Problema 13
Suponga que la PTAR de la UAO se puede modelar como un reactor tipo CMFR
de volumen V por el que fluye un cierto caudal Q. El agua se trata mediante una
reacción de primer orden con el fin de reducir la concentración de contaminantes
desde un valor inicial C0 hasta un valor final Cf .
Solución
Q 1 C
1. C(t) = C0 e−( V +k1 )t −→ Q = V − ln − k1
t C0
1 40 1
Q = 120 − ln − 0.2 = 555.4 m3 /day× = 6.43 × 10−3 m3 /s
8/24 200 86400
355
2. z1 − z2 = 10 m accesorios bridados ε = 0.12 mm
L1 = 50 m; d1 = 15 cm = 600 L2 = 10 m; d2 = 10 cm = 400
Esquema instalación
p1 V12 p2 V 2 V22
+ + z1 = + 2 + z2 + hw + hT ⇒ hw = (z1 − z2 ) − − hT
ρg 2g ρg 2g 2g
1
Haaland f= " 1.11 #!2
6.9 ε/d
−1.8 log10 +
Red 3.7
f1 ≈ 0.0229 f2 ≈ 0.0230
356
Coeficientes de pérdida de los accesorios
Tramo 1 6”
Accesorio # Vr. Unit ζ
Admisión 1 0.8 0.8
0.64 + 0.58
T ramificada 3 1.83
2
0.3 + 0.26
Codo estándar 4 1.12
2
0.16 + 0.07
Válvula compuerta 1 0.115
2
(1)
ζT = 3.865
Tramo 2 4”
Accesorio # Vr. Unit ζ
Válvula ángulo 1 2 2
Codo estándar 3 0.3 0.9
Codo ext. larga 1 0.19 0.16
(2)
ζT = 3.09
0.362
(1) 50
hT = 0.0229 × + 3.865 × = 0.076 m
0.15 2 × 9.8
0.822
(2) 10
hT = 0.023 × + 3.09 × = 0.185 m
0.1 2 × 9.8
0.822
hw = 10 − − 0.076 − 0.085 = 9.7 m >0 OK
2 · 9.8
357
Problema 14
El agua sin tratar se extrae de un tanque utilizando una cierta bomba a través
de una conducción de hierro galvanizado (ε = 0.15 mm), llevándose a la cámara
de contacto de cloro. El primer tramo, desde el tanque a la bomba tiene un
diámetro d1 y longitud L1 mientras que el segundo, desde la bomba a la cámara,
tiene un diámetro d2 y longitud L2 . El agua del tanque se conduce hasta
la bomba utilizando 4 codos de extensión larga, una curva de retorno y una
válvula de globo (abierta a la mitad). Desde la bomba hasta la cámara de
contacto se encuentra una válvula de compuerta (abierta), una T ramificada
y 3 codos estándar. La tuberı́a descarga a la presión atmosférica, unos pocos
centı́metros sobre el nivel de la cámara de contacto. Todos los accesorios se
encuentran bridados.
Nota: Suponga que la cámara se puede modelar como un CMFR con Cin =
10 mg/L, V = 300 m3 .
Nota: La concentración de cloro en la cámara evoluciona como:
C(t) = Cin [1 − exp (−Qt/V )].
358
Solución:
300 × 103
2
Q=− ln 1 − ≈ 74.4 L/s = 0.0744 m3 /s
15 × 60 10
2. ε = 0.15 mm d1 = 10 cm = 4 in L1 = 100 m
d2 = 15 cm = 6 in L2 = 20 m
Esquema instalación
f1 ≈ 0.0220 f2 ≈ 0.0201
359
Coeficientes de pérdida de los accesorios
Tramo 1
1) Accesorio # Vr. Unit ζ
Entrada 1 0.5 0.5
Codo ext Larga 4 0.19 0.76
Curva retorno 1 0.3 0.3
Válvula globo 1/2 1 15 15
abierta
(1)
ζT = 16.56
Tramo 2
2) Accesorio # Vr. Unit ζ
Válvula compuerta 1 0.115 0.115
abierta
T ramificada 1 0.61 0.61
Codo estándar 3 0.28 0.84
(2)
ζT = 1.565
4.212
hb = + 15 + 180.27 ≈ 196.2 m
2 × 9.8
Ẇideal 143
Ẇreal = = ≈ 204.3 kW
η 0.7
360
Problema 15
• Halle el caudal necesario para que el agua a la salida del reactor haya
alcanzado una concentración de contaminante de 20 mg/L en 45 minutos.
Exprese su resultado en L/s.
El agua efluente del reactor se utiliza para recuperar energı́a a través de una
turbina hidráulica y se conduce hasta una laguna de tratamiento. El agua se
conduce a través de una conducción de acero laminado nuevo (ε = 0.05 mm).
El primer tramo, desde el reactor a la turbina tiene un diámetro d1 y longitud
L1 mientras que el segundo, desde la turbina a la laguna, tiene un diámetro
d2 y longitud L2 . El agua del reactor se conduce hasta la turbina utilizando 6
codos estándar, dos codos de extensión larga y una válvula de ángulo (abierta).
Desde la turbina hasta la laguna se encuentra una válvula de retención de disco
oscilante (abierta), dos curvas de retorno y 3 T ramificadas. La tuberı́a descarga
a la presión atmosférica, unos pocos centı́metros sobre el nivel de la laguna.
Todos los accesorios se encuentran bridados.
Nota: Suponga que la cámara se puede modelar como un CMFR con Cin =
100 mg/L, V = 500 m3 .
Nota: La concentración de cloro en la cámara evoluciona como
C(t) = Cin [1 − exp (−Qt/V )].
361
Solución:
500 × 103
20
Q=− ln 1 − ≈ 41.3 L/s = 0.0413 m3 /s
45 × 60 100
2. ε = 0.05 mm d1 = 15 cm = 6 in L1 = 50 m
d2 = 10 cm = 4 in L2 = 20 m
Esquema instalación
f1 ≈ 0.0169 f2 ≈ 0.0175
362
Coeficientes de pérdida de los accesorios
Tramo 1
1) Accesorio # Vr. Unit ζ
Entrada 1 0.5 0.5
Codo estándar 6 0.28 1.68
Codo ext. Larga 2 0.17 0.34
Válvula ángulo 1 2.0 2.0
abierta
(1)
ζT = 4.52
Tramo 2
2) Accesorio # Vr. Unit ζ
Válvula retención 1 2 2.0
Curva retorno 2 0.3 0.6
T Ramificada 3 0.64 1.92
(2)
ζT = 4.52
5.262
hS = 40 − − 14.2 ≈ 24.39 m
2 × 9.8
363
Problema 16
• A partir de los datos proporcionados más abajo estime el caudal que circula
por el tanque en f t3 /day si el tiempo de tratamiento fue de 8 horas.
Solución
dC Q Q Q Q
CMFR = Cin − C − kC = −C +k → C = C0 e−( V +k)t
dt V V V
364
Cf Q 1 Cf Q 1 Cf
ln =− + k t → − ln = +k → Q = −V ln +k
C0 V t C0 V t C0
1 5
Q = −1000 ln − 0.16 = 8827.2 f t3 /day
1/3 100
1
Q = 8827.2 · = 0.102 f t3 /s
86400
Esquema instalación
p1 V12 p2 V 2 p2
+ + z
1 = + 2 + z2 + hw + hT ⇒ hw = − − z2 − hT
ρg
2g
ρg 2g
ρg
f1 ≈ 0.0221 f2 ≈ 0.0213
365
Coeficientes de pérdida de los accesorios
Tramo 1 3”
Accesorio # Vr. Unit ζ
Admisión entrante 1 0.8 0.8
0.95 + 0.64
Codo estándar 5 3.975
2
1.4 + 1.0
T ramificada 2 2.5
2
0.16 + 0.11
Válvula compuerta 1 0.135
2
abierta
(1)
ζT = 7.41
Tramo 2 2”
Accesorio # Vr. Unit ζ
Válvula globo 1 6.9 6.9
abierta
Codo ext. larga 2 0.41 0.82
Codo 45◦ 4 0.3 1.2
Salida 1 1 1.0
(2)
ζT = 9.92
2.082
(1) 30
hT = 0.0221 × + 7.41 × = 0.68 m
3/12 2 × 32.2
4.672
(2) 300
hT = 0.0213 × + 9.92 × = 16.34 m
2/12 2 × 32.2
(1) (2)
hT = hT + hT = 0.68 + 16.34 = 17.02 f t
10 × 144
hw = − − 25 − 17.02 = −65.10 f t (OK Bomba)
62.4
f t lbf
Ẇideal = γ Q hw = 62.4 × 0.102 × 65.10 ≈ 414.35
s
366
Problema 17
Solución
367
Cout V −k V
ln = −k → Q=
Cin Q Cout
ln
Cin
4 1
Q = −0.36 · 50000/ ln = 6008.5 f t3 /day × = 250.36 f t3 /h
80 24
1
Q = 250.36 f t3 /h × = 0.07 f t3 /s
3600
Esquema instalación
p1 V12 p2 V2 V2
+ + z1 = + 2 + z
2 + hw + hT ⇒ hw = z1 − 2 − hT
ρg 2g ρg 2g 2g
368
f1 ≈ 0.0656 f2 ≈ 0.0786
Tramo 2 2”
Accesorio # Vr. Unit ζ
Válvula ángulo 1 2.4 2.4
abierta
Codo estándar 3 0.39 1.17
T lı́nea 2 0.19 0.38
(2)
ζT = 3.95
1.432
(1) 100
hT = 0.0656 × + 8.565 × = 1.10 m
3/12 2 × 32.2
3.212
(2) 200
hT = 0.0786 × + 3.95 × = 15.72 m
2/12 2 × 32.2
(1) (2)
hT = hT + hT = 1.10 + 15.72 = 16.82 f t
3.212
hw = 100 − − 16.82 ≈ 83 f t (T urbina OK)
2 · 32.2
f t lbf
Ẇideal = γ Q hw = 62.4 × 0.07 × 83 ≈ 362.5
s
f t lbf 1
Ẇreal = Ẇideal η = 362.5 · 0.91 ≈ 330 · ≈ 0.6 hp
s 550
369
Problemas Propuestos
Problema 1
Datos: ρ = 1.94 slug/f t3 , g = 32.2 f t2 /s, µ = 2.34 × 10−5 slug f t−1 s−1
Problema 2
370
• Teniendo en cuenta que el reactor está situado 10 m por encima de la
laguna de oxidación, encuentre la potencia de la bomba en kW necesaria
para impulsar el caudal que usted calculó anteriormente, asumiendo que
la máquina tiene una eficiencia del 85%.
Problema 3
371
Datos: Cin = 100 mg/L, Cout = 3.5 mg/L, V = 20 m3 , k2 = 2×10−3 L/(mg s),
ρagua = 1010 kg/m3 , ν = 10−6 m2 /s, d1 = 10 cm, d2 = 7.5 cm, L1 = 200 m,
L2 = 80 m.
Nota: Considere ambos tanques están abiertos a la atmósfera. La entrada al
primer tramo es una admisión entrante.
Problema 4
• A partir de los datos proporcionados más abajo estime el caudal que circula
por el tanque en f t3 /day.
372
Nota: Considere que tanto el tanque como el depósito están abiertos a la
atmósfera. La entrada al primer tramo es una admisión entrante.
Problema 5
• A partir de los datos proporcionados más abajo estime el caudal que circula
por el tanque en m3 /day.
373
Problema 6
Como se muestra en la figura, se bombea agua entre dos depósitos con un caudal
de 0.2 f t2 /s a través de una conducción de 400 f t de longitud y 2 in de diámetro
incluyendo accesorios. Si la rugosidad relativa ε/d = 0.001 halle la potencia en
hp que debe sumistrar la bomba.
Problema 7
374
Datos: Cin = 100 mg/L, Cout = 5 mg/L, V = 10 m3 , k2 = 10−3 L/(mg s),
ρagua = 1000 kg/m3 , ν = 10−6 m2 /s, d1 = 7.5 cm, d2 = 5 cm, L1 = 10 m,
L2 = 100 m.
Nota: Considere que tanto la laguna como el reactor están abiertos a la
atmósfera. La entrada al primer tramo es una admisión de borde cuadrado
(arista viva)
Problema 8
Suponga que la PTAR de la UAO se puede modelar como un reactor tipo CMFR
de volumen V por el que fluye un cierto caudal Q. El agua se trata mediante una
reacción de segundo orden de constante k con el fin de reducir la concentración
de contaminantes desde su valor a la entrada Cin hasta un valor Cf = Cin /10.
• Con los datos proporcionados más abajo ¿Cuál es el caudal que circula
por la PTAR en f t3 /s?
375
Problema 9
Para sacar el agua contaminada de la laguna se utiliza una cierta bomba a través
de una conducción de hierro fundido asfaltado (ε = 3 × 10−5 f t) de diámetro
d1 y longitud L1 . La tuberı́a descarga a un contenedor grande presurizado a
25 psi (gage) y contiene los siguientes elementos: entrada con admisión de borde
cuadrado, cuatro codos estándar de 90◦ , válvula de globo (abierta a la mitad),
bomba, una válvula de retención de disco oscilante (abierta) y dos codos de 45◦ .
Todos los accesorios se encuentran roscados.
Problema 10
376
Problema 11
377
en kW necesaria para impulsar el caudal que usted calculó en el primer
punto, asumiendo que la máquina tiene una eficiencia del 78%.
Datos: Cin = 100 mg/L, Cout = 3.5 mg/L, V = 20 m3 , k2 = 2×10−3 L/(mg s),
ρagua = 1010 kg/m3 , ν = 10−6 m2 /s, d1 = 10 cm, d2 = 7.5 cm, L1 = 200 m,
L2 = 80 m.
Nota: Considere que ambos tanques están abiertos a la atmósfera. La entrada
al primer tramo es una admisión entrante.
378