INSTITUTO TECNOLÓGICO DE MÉXICO

EN CELAYA

PORTAFOLIO DE EVIDENCIAS

MECANISMOS DE
TRANSFERENCIA
RODRÍGUEZ MARÍN ITZEL
No. CONTROL 18031869

CÁRDENAS RIVERA JOSÉ CARLOS

INGENIERÍA QUÍMICA
23 DE JUNIO DE 2021
UNIDAD
UNO
EJERCICIOS DE CLASE
1.1 Calcular la viscosidad del nitrógeno a 50°C y 854 atm, donde la presión crítica del nitrógeno es
33.5 atm y la temperatura critica de 126°K.

𝑁2 : 𝑀 = 28
𝑇𝑐 = 126.2°𝐾
𝑃𝑐 = 33.5 𝑎𝑡𝑚
𝑇 = 50°𝐶 + 273
𝑃 = 854 𝑎𝑡𝑚
𝑇 323
𝑇𝑟 = = = 2.56
𝑇𝑐 126.2
𝑃 854
𝑃𝑟 = = = 25.8
𝑃𝑐 33.5
𝑀𝑐 = 7.7𝑀1⁄2 𝑃𝑐 2⁄3 𝑇 −1⁄6

𝑀𝑐 = 7.7(28)1⁄2 (33.5)2⁄3 (126.2)−1⁄6

𝑀𝑐 = 189.0554𝑚𝑖𝑐𝑟𝑜𝑝𝑜𝑖𝑠𝑒𝑠 ≈ 189𝑥10−6 𝑝𝑜𝑖𝑠𝑒𝑠
𝜇
𝜇 ∗ (𝑡𝑎𝑏𝑙𝑎 1.3 − 1) = 2.6 𝜇∗ =
𝜇𝑐
𝜇 = 2.6(189𝑥10−6 )
𝜇 = 4.9𝑥10−4 𝑝𝑜𝑖𝑠𝑒𝑠

1.2 La viscosidad del CO2 a 45.3 atm y 40.3°C es 1800x10-7 poises. Estimar el valor de la viscosidad
a 114.6 atm y 40.3°C utilizando una Pc=72.9 atm y una Tc=304.2°K.
Estado 1 Estado 2
CO2 CO2
P = 45.3 atm P = 114.6 atm
T= 40.3°C + 273 T= 40.3°C + 273
µ = 1800x10-7 poises µ = ¿?
𝑇𝑐 = 304.3 °𝐾
𝑃𝑐 = 72.9 𝑎𝑡𝑚
Estado 1:
313.3
𝑇𝑟 = = 1.03
304.2
45.3
𝑃𝑐 = = 0.622
72.9
 𝜇 # = 1.2 (𝑇𝑎𝑏𝑙𝑎 1.3 − 2)

#
𝜇 𝜇 1800𝑥10−7
𝜇 = ∴ 𝜇𝑜 = # = = 1.5𝑥10−4 𝑝𝑜𝑖𝑠𝑒.
𝜇𝑜 𝜇 1.2
Estado 2:
313.3
𝑇𝑟 = = 1.03
304.2
114.6
𝑃𝑐 = = 1.57
72.9
𝜇
𝜇 # = 4.8 ∴ 𝜇 # =
𝜇𝑜
𝜇 = 𝜇 # × 𝜇𝑜 = 4.8(1.5𝑥10−4 ) 𝑝𝑜𝑖𝑠𝑒.
𝜇 = 7.2𝑥10−4 𝑝𝑜𝑖𝑠𝑒𝑠
1.3 Predecir la viscosidad del oxígeno, nitrógeno y metano moleculares a la presión atmosférica y
20°C. Expresar todos los resultados en centipoises.

√𝑀𝑇
𝜇 = 2.6693𝑥10−5
𝜎 2 𝑅𝜇
𝑂2 : 𝑇 = 20°𝐶 + 273 = 293°𝐾
𝑃 = 1 𝑎𝑡𝑚
𝑀 = 32
𝜎: (𝐴𝑝𝑒𝑛𝑑𝑖𝑐𝑒 𝐾, 𝑝𝑎𝑔 912 − 913)
= 3.435 Å
∈
= 113°𝐾
𝐾
𝐾 293
(𝑇) = = 2.59
∈ 113
𝐾 R𝜇
(𝑇)
∈
2.50 1.093
2.51 1.0822
2.60 1.081

= (3.933)2 = 11.98𝐾 2
√32(293)
𝜇 = 2.6693𝑥10−5
11.78 ∗ 1.0822
𝑔𝑟
𝜇 = 2.027𝑥10−4
𝑐𝑚 − 𝑠
𝜇 = 2.027𝑥10−4 𝑝𝑜𝑖𝑠𝑒
𝜇 = 0.02027 𝑐𝑝(𝑐𝑒𝑛𝑡𝑖𝑝𝑜𝑖𝑠𝑒𝑠)

𝑁2 : 𝑇 = 20°𝐶 + 273 = 293°𝐾

 𝑃 = 1 𝑎𝑡𝑚
𝑀 = 28
𝜎 = 3.681 Å 𝜎 2 = 13.5497 Å
∈
= 91.5 °𝐾
𝐾
𝐾 293
(𝑇) = = 3.2021
∈ 91.5
𝐾 R𝜇
(𝑇)
∈
3.2 1.022
3.2021 
3.30 1.014

√28(293)
𝜇 = 2.6693𝑥10−5
13.54 ∗ 1.022
𝑔𝑟
𝜇 = 1.7459𝑥10−4
𝑐𝑚 − 𝑠
𝜇 = 1.7459𝑥10−4 𝑝𝑜𝑖𝑠𝑒
𝜇 = 0.017459 𝑐𝑝(𝑐𝑒𝑛𝑡𝑖𝑝𝑜𝑖𝑠𝑒𝑠)

𝐶𝐻4 : 𝑇 = 20°𝐶 + 273 = 293°𝐾

𝑃 = 1 𝑎𝑡𝑚
𝑀 = 16.04
𝜎 = 3.822 Å 𝜎 2 = 13.5497 Å
∈
= 136.5 °𝐾
𝐾
𝐾 293
(𝑇) = = 2.1465
∈ 136.5
𝐾 R𝜇
(𝑇)
∈
2.10 1.156
2.1465 1.1476
2.20 1.138

√16(293)
𝜇 = 2.6693𝑥10−5
14.6076 ∗ 1.147
𝑔𝑟
𝜇 = 1.0908𝑥10−4
𝑐𝑚 − 𝑠
𝜇 = 1.0908𝑥10−4 𝑝𝑜𝑖𝑠𝑒
𝜇 = 0.010908 𝑐𝑝(𝑐𝑒𝑛𝑡𝑖𝑝𝑜𝑖𝑠𝑒𝑠)
1.4 Predecir la viscosidad de la siguiente gaseosa a 1 atm y 293°K, a partir de los datos de los
componentes puros.
Sustancia Fracción Molar (x) Peso Molecular (M) Viscosidad (poise)µ i j
CO2 0.133 44.01 1462 x10-7 1 1,2,3
O2 0.039 32.00 2031 x10-7 2 1,2,3
N2 0.828 20.016 1754 x10-7 3 1,2,3

𝑛
𝑥𝑖𝜇𝑖
𝑀𝑚𝑒𝑧𝑐𝑙𝑎 = ∑
∑𝑛𝑗−1 𝑥𝑗∮𝑖𝑗
𝑖−1

1 1 1 2
1 𝑀𝑖 −2 𝑀𝑖 2 𝑀𝑗 4
∮𝑖𝑗 = (1 + ) [1 + ( ) ( ) ]
√8 𝑀𝑗 𝑀𝑗 𝑀𝑖

1 1 1 2
1 44 −2 1462x10−7 2 44 4
𝐼11 = (1 + ) [1 + ( ) ( ) ] =1
√8 44 1462x10−7 44

1 1 1 2
−
1 44 2 1462x10−7 2
32 4
𝐼12 = (1 + ) [1 + ( −7
) ( ) ] = 0.7296
√8 32 2031x10 44

1 1 1 2
−
1 44 2 1462x10−7 2
28 4
𝐼13 = (1 + ) [1 + ( −7
) ( ) ] = 0.7267
√8 28 1754x10 44

1 1 1 2
1 32 −2 2031x10−7 2 44 4
𝐼21 = (1 + ) [1 + ( ) ( ) ] = 1.396
√8 44 1462x10−7 32

1 1 1 2
−
1 32 2 2031x10−7 2 32 4
𝐼22 = (1 + ) [1 + ( −7
) ( ) ] =1
√8 31 2031x10 32

1 1 1 2
1 32 −2 2031x10−7 2 28 4
𝐼23 = (1 + ) [1 + ( ) ( ) ] = 1.0073
√8 28 1754x10−7 32

1 1 1 2
−
1 28 2 1754x10−7 2
44 4
𝐼31 = (1 + ) [1 + ( −7
) ( ) ] = 1.3699
√8 44 1462x10 28
 1 1 1 2
−
1 28 2 1754x10−7 2
32 4
𝐼32 = (1 + ) [1 + ( −7
) ( ) ] = 0.9927
√8 32 2031x10 28

1 1 1 2
−
1 28 2 1754x10−7 2 28 4
𝐼31 = (1 + ) [1 + ( −7
) ( ) ] =1
√8 28 1754x10 28

𝑥1 𝜇 1 𝑥2 𝜇 2 𝑥3 𝜇 3
𝑀𝑚𝑒𝑧𝑐𝑙𝑎 = 𝑥 +𝑥 +𝑥
1 𝐼11 +𝑥2 𝐼12 +𝑥3 𝐼13 1 𝐼21 +𝑥2 𝐼22 +𝑥3 𝐼23 1 𝐼31 +𝑥2 𝐼32 +𝑥3 𝐼33

0.133(1462x10−7 )
𝑀𝑚𝑒𝑧𝑐𝑙𝑎 = +
0.133(1) + 0.039(0.7288) + 0.828(0.7266)
0.039(2031x10−7 )
+
0.133(1.396)+0.039(1)+0.828(1.0078)

0.828(1754x10−7 )
0.133(1.3699) + 0.039(0.9927) + 0.828(1)
𝑀𝑚𝑒𝑧𝑐𝑙𝑎 = 1.7144𝑥10−4 𝑝𝑜𝑖𝑠𝑒𝑠
Viscosidad en líquidos.
1.5 Determinar la viscosidad del Benceno líquido a 20°C.
𝑇𝑏 = 80.1 °𝐶
𝐵𝑒𝑛𝑐𝑒𝑛𝑜 { 𝑐𝑚3 }
𝜎 = 89 (20°𝐶)
𝑚𝑜𝑙
6.023x1023 (6.624x10−27 ) 3.8(353.1)
𝜇= 𝑒 293
89
𝑔𝑟
𝜇 = 4.3689x10−3
𝑐𝑚 − 𝑠
Distribución de la velocidad.
Balances de Cantidad de movimiento.
1.6 Coordenadas rectangulares - Condiciones de operaciones.
 Estado estacionario
 Flujo luminar
 Una sola dirección
 µ y ρ constantes
 Unidireccional
 𝑣𝑧
𝐹𝑚 = 𝑚𝑎 = 𝑚 ( )
𝑡
𝐹𝑚 = ṁvz
ṁ
𝜌= ∴ ṁ = ρv
v
Fm = ρv(flujo)vz

v̇ (flujo) = 𝑣𝑧𝐴(𝑎𝑟𝑒𝑎)

𝐹𝑚 = 𝜌𝑣𝑧𝑣𝑧∆𝑥∆𝑦
𝐹𝓉
𝓉𝑦𝑧 = ∴ 𝐹𝓉 = 𝓉𝑦𝑧 𝐴
𝐴
𝐹𝑢𝑒𝑟𝑧𝑎 𝓉 = 𝓉𝑦𝑧 ∆𝑥∆𝑦

𝐹𝑔 = 𝑚𝑔
𝑚
𝜌= ∴ 𝑚 = 𝜌𝑣
𝑣
𝐹𝑔 = 𝑔𝜌∆𝑥∆𝑦 ← 𝑒𝑐. ℎ𝑜𝑟𝑖𝑧𝑜𝑛𝑡𝑎𝑙
𝐹𝑔 = 𝑔𝜌∆𝑥∆𝑦∆𝑧𝑐𝑜𝑠𝜃
Ecuación general.
𝐸𝑛𝑡𝑟𝑎𝑑𝑎 𝑑𝑒 𝑐𝑎𝑛𝑡𝑖𝑑𝑎𝑑 𝑑𝑒 𝑚𝑜𝑣𝑖𝑚𝑖𝑒𝑛𝑡𝑜 𝑆𝑎𝑙𝑖𝑑𝑎 𝑑𝑒 𝑐𝑎𝑛𝑡𝑖𝑑𝑎𝑑 𝑑𝑒 𝑚𝑜𝑣𝑖𝑚𝑖𝑒𝑛𝑡𝑜
−
𝐹𝑚 𝐹𝓉 𝐹𝑃 𝐹𝑚 𝐹𝓉 𝐹𝑃
∑ 𝐹 𝑞𝑢𝑒 𝑎𝑐𝑡𝑢𝑎𝑛 𝑒𝑛 𝑒𝑙 𝑠𝑖𝑠𝑡𝑒𝑚𝑎
+ =0
𝐹𝑔

𝜌𝑣𝓏𝑣𝓏∆𝑥∆𝑦|𝑧 − 𝜌𝑣𝓏𝑣𝓏∆𝑥∆𝑦|𝑧+∆𝑧 + 𝓉𝑦𝑧 ∆𝑥∆𝑧|𝑦 − 𝓉𝑦𝑧 ∆𝑥∆𝑧|𝑦+∆𝑦 + 𝜌𝑔∆𝑥∆𝑦∆𝑧𝑐𝑜𝑠𝜃 = 0

[𝓉𝑦𝑧 ∆𝑥∆𝑧|𝑦 − 𝓉𝑦𝑧 ∆𝑥∆𝑧|𝑦+∆𝑦 + 𝜌𝑔∆𝑥∆𝑦∆𝑧𝑐𝑜𝑠𝜃]

=0
∆𝑥∆𝑦∆𝑧
𝓉𝑦𝑧 |𝑦 − 𝓉𝑦𝑧 |𝑦+∆𝑦
+ 𝜌𝑔𝑐𝑜𝑠𝜃 = 0
∆𝑦
−𝓉𝑦𝑧 |𝑦 + 𝓉𝑦𝑧 |𝑦+∆𝑦
= 𝜌𝑔𝑐𝑜𝑠𝜃
∆𝑦
−𝓉𝑦𝑧 |𝑦 + 𝓉𝑦𝑧 |𝑦+∆𝑦
lim = 𝜌𝑔𝑐𝑜𝑠𝜃
∆𝑦→0 ∆𝑦
𝑑𝓉𝑦𝑧
= 𝜌𝑔𝑐𝑜𝑠𝜃 ∴ ∫ 𝑑𝓉𝑦𝑧 = 𝜌𝑔𝑐𝑜𝑠𝜃 ∫ 𝑑𝑦
𝑑𝑦
𝓉𝑦𝑧 = 𝜌𝑔𝑐𝑜𝑠𝜃(𝑦) + 𝑐1
 𝐶. 𝐹. 1 𝓉𝑦𝑧 = 0 𝑦=0
0 = 0 + 𝑐1 ∴ 𝑐1 = 0
⟹ 1 . 𝓉𝑦𝑧 = 𝜌𝑔𝑐𝑜𝑠𝜃

Ley de Newton
𝑑𝑣𝓏
𝓉𝑦𝑧 = −𝜇 → (−)𝑝𝑜𝑟𝑞𝑢𝑒 𝑒𝑠𝑡𝑎 𝑑𝑖𝑠𝑚𝑖𝑛𝑢𝑦𝑒𝑛𝑑𝑜
𝑑𝑦
𝑑𝑣𝓏
−𝜇 = 𝜌𝑔𝑐𝑜𝑠𝜃
𝑑𝑦
𝜌𝑔𝑐𝑜𝑠𝜃
∫ 𝑑𝑣𝓏 = − ∫ 𝑦𝑑𝑦
𝜇
𝜌𝑔𝑐𝑜𝑠𝜃𝑦 2
𝑣𝓏 = − + 𝑐2
2𝜇
𝐶. 𝐹. 2 𝑣𝓏 = 0 𝑦=𝛿
2
𝑐2 = 𝜌𝑔𝑐𝑜𝑠𝜃(𝛿 )/2𝜇
𝜌𝑔𝑐𝑜𝑠𝜃𝑦 2 𝜌𝑔𝑐𝑜𝑠𝜃(𝛿 2 )
𝑣𝓏 = − +
2𝜇 2𝜇
𝜌𝑔𝑐𝑜𝑠𝜃 2
⟹ 2 . 𝑣𝓏 = [𝛿 − 𝑦 2 ]
2𝜇

𝐶. 𝐹. 3 𝑣𝓏 = 𝑣𝑚𝑎𝑥𝑖𝑚𝑎 𝑦=0
𝜌𝑔𝑐𝑜𝑠𝜃(𝛿 2 )
⟹ 3 . 𝑣𝓏𝑚𝑎𝑥𝑖𝑚𝑎 =
2𝜇

𝑄 = 𝑣̇ = 𝑣𝓏𝐴
𝑑𝑄 = 𝑣𝓏𝑑𝑦𝑑𝑥
𝑎 𝑤 𝛿
∫ 𝑑𝜃 = ∫ ∫ 𝑣𝓏𝑑𝑦𝑑𝑥
0 0 0
𝑎 𝑤 𝛿
𝜌𝑔𝑐𝑜𝑠𝜃 2
∫ 𝑑𝜃 = ∫ ∫ [𝛿 − 𝑦 2 ]𝑑𝑦𝑑𝑥
0 0 0 2𝜇
𝑤 𝛿
𝜌𝑔𝑐𝑜𝑠𝜃 2
𝑄=∫ ∫ [𝛿 − 𝑦 2 ]𝑑𝑦𝑑𝑥
0 0 2𝜇
 𝛿
𝜌𝑔𝑐𝑜𝑠𝜃 2
→ ∫ [𝛿 − 𝑦 2 ]𝑑𝑦
0 2𝜇

𝜌𝑔𝑐𝑜𝑠𝜃 𝛿 2
∫ [𝛿 − 𝑦 2 ]𝑑𝑦
2𝜇 0
𝛿 𝛿
∫ 𝛿 2 𝑑𝑦 − ∫ 𝑦 2 𝑑𝑦
0 0
𝛿 𝛿
𝛿 ∫ 𝑑𝑦 − ∫ 𝑦 2 𝑑𝑦
2
0 0

2 [𝑦]𝛿
𝑦3 𝛿
𝛿 −[ ]
0 3 0

𝛿3
𝛿 2 (𝛿) −
3

3
𝛿3
𝛿 −
3
3𝛿 3 − 𝛿 3 2𝛿 3
=
3 3
𝑤
2𝛿 3
→∫ 𝑑𝑥
0 3

2𝛿 3 𝑤
∫ 1 𝑑𝑥
3 0

2𝛿 3 𝑤
[𝑥]
3 0
2𝛿 3 𝑤
3
𝑤 𝛿
𝜌𝑔𝑐𝑜𝑠𝜃 2
𝑄=∫ ∫ [𝛿 − 𝑦 2 ]𝑑𝑦𝑑𝑥
0 0 2𝜇
2𝛿 3 𝑤
𝜌𝑔𝑐𝑜𝑠𝜃
𝑄= 3
2𝜇
2𝜌𝑔𝑐𝑜𝑠𝜃𝛿 3 𝑤
𝑄=
6𝜇
𝜌𝑔𝑐𝑜𝑠𝜃𝛿 3 𝑤
⟹ 4. 𝑄 =
3𝜇
 𝑄
<𝑣 >= 𝑤 𝛿
∫0 ∫0 𝑑𝑦𝑑𝑥

𝜌𝑔𝑐𝑜𝑠𝜃𝛿 3 𝑤
3𝜇
<𝑣 >=
𝑤𝛿
𝜌𝑔𝑐𝑜𝑠𝜃𝛿 2
⟹5. <𝑣 >=
3𝜇
𝐿 𝑤
𝐹𝑅 = ∫ ∫ 𝓉𝑦𝑧 |𝑧=𝛿 𝑑𝑧𝑑𝑥
0 0
𝐿 𝑤
𝑑𝑣
𝐹𝑧 = ∫ ∫ (−𝜇 ) | 𝑑𝑧𝑑𝑥
0 0 𝑑𝑦 𝑧=𝛿
𝐿 𝑤
= ∫ ∫ (𝑦𝜌𝑔𝑐𝑜𝑠𝜃)|𝑧=𝛿 𝑑𝑧𝑑𝑥
0 0

⟹ 6 . 𝐹𝑧 = 𝐿𝑤(𝛿𝜌𝑔𝑐𝑜𝑠𝜃)
1.7
Condiciones de operación
 Unidireccional
 Flujo laminar
 Estado estacionario
 ρ y µ constantes
 Isotérmico

𝜌𝑣𝑧𝑣𝑧∆𝑥∆𝑦|𝑧 − 𝜌𝑣𝑧𝑣𝑧∆𝑥∆𝑦|𝑧+∆𝑧 + 𝓉𝑦𝑧 ∆𝑥∆𝑧|𝑦 − 𝓉𝑦𝑧 ∆𝑥∆𝑧|𝑦+∆𝑦 + 𝑃𝑜∆𝑥∆𝑦|𝑧 − 𝑃𝐿∆𝑥∆𝑦|𝑧+∆𝑧 = 0

𝓉𝑦𝑧 ∆𝑥∆𝑧|𝑦 − 𝓉𝑦𝑧 ∆𝑥∆𝑧|𝑦+∆𝑦 + 𝑃𝑜∆𝑥∆𝑦|𝑧 − 𝑃𝐿∆𝑥∆𝑦|𝑧+∆𝑧
=0
∆𝑥∆𝑦∆𝑧
𝓉𝑦𝑧 |𝑦 − 𝓉𝑦𝑧 |𝑦+∆𝑦 𝑃𝑜|𝑧 − 𝑃𝐿|𝑧+∆𝑧
+ =0
∆𝑦 ∆𝑧
𝐿𝑖𝑚 𝓉𝑦𝑧 |𝑦 − 𝓉𝑦𝑧 |𝑦+∆𝑦 𝑃𝑜|𝑧 − 𝑃𝐿|𝑧+∆𝑧
+
∆𝑦 → 0 ∆𝑦 ∆𝑧
 𝑑𝓉𝑦𝑧 ∆𝑃
=− ∆𝑧 = 𝐿
𝑑𝑦 ∆𝑧
𝑑𝓉𝑦𝑧 −𝑃𝑜 − 𝑃𝐿
=
𝑑𝑦 𝐿
𝑑𝓉𝑦𝑧 𝑃𝑜 − 𝑃𝐿
=
𝑑𝑦 𝐿
𝑃𝑜 − 𝑃𝐿
∫ 𝑑𝓉𝑦𝑧 = ∫ 𝑑𝑦
𝐿
𝑃𝑜 − 𝑃𝐿
𝓉𝑦𝑧 = ( )𝑦 + 𝑐
𝐿
𝐶. 𝐹. 1 𝓉𝑦𝑧 = 0 𝑦 = 0 𝑐 = 0
𝑃𝑜 − 𝑃𝐿
⟹ 1 . 𝓉𝑦𝑧 = ( )𝑦
𝐿
𝑑𝑣𝑧
𝓉𝑦𝑧 = ±𝜇
𝑑𝑦
𝑑𝑣𝑧
𝓉𝑦𝑧 = −𝜇
𝑑𝑦
𝑃𝑜 − 𝑃𝐿 𝑑𝑣𝑧
( ) 𝑦 = −𝜇
𝐿 𝑑𝑦
𝑃𝑜 − 𝑃𝐿
∫ 𝑑𝑣𝑧 = − ∫ 𝑦𝑑𝑦
𝜇𝐿
𝑃𝑜 − 𝑃𝐿 𝑦 2
𝑣𝑧 = − +𝑐
𝜇𝐿 2
𝐶. 𝐹. 2 𝑣𝑧 = 0 𝑦=ℎ
𝑃𝑜 − 𝑃𝐿 ℎ2
0=− +𝑐
𝜇𝐿 2
(𝑃𝑜 − 𝑃𝐿)ℎ2
𝑐=
2𝜇𝐿
𝑃𝑜 − 𝑃𝐿 𝑦 2 (𝑃𝑜 − 𝑃𝐿)ℎ2
𝑣𝑧 = − +
𝜇𝐿 2 2𝜇𝐿
𝑃𝑜 − 𝑃𝐿 2
⟹ 2 . 𝑣𝑧 = [ℎ − 𝑦 2 ]
2𝜇𝐿

𝐶. 𝐹. 2 𝑣𝑧 = 𝑣𝑥𝑚𝑎𝑥𝑖𝑚𝑎 𝑦=0
 𝑃𝑜 − 𝑃𝐿 2
⟹ 3. 𝑣𝑧𝑚á𝑥 = [ℎ ]
2𝜇𝐿

𝑑𝑄 = 𝑣𝑥𝑑𝐴
𝑄 𝑤 𝐿
∫ 𝑑𝑄 = ∫ ∫ 𝑣𝑧𝑑𝑦𝑑𝑥
0 0 0
𝑤 𝐿
𝑃𝑜 − 𝑃𝐿 2
𝑄=∫ ∫ [ℎ − 𝑦 2 ]𝑑𝑦𝑑𝑥
0 0 2𝜇𝐿
𝐿
𝑃𝑜 − 𝑃𝐿 2
→ ∫ [ℎ − 𝑦 2 ]𝑑𝑦
0 2𝜇𝐿

𝑃𝑜 − 𝑃𝐿 𝐿 2
∫ [ℎ − 𝑦 2 ]𝑑𝑦
2𝜇𝐿 0

𝑃𝑜 − 𝑃𝐿 𝑤 2 𝑦3 𝐿
𝑄= ∫ [ℎ 𝑦 − ]
2𝜇𝐿 0 3 0

𝑃𝑜 − 𝑃𝐿 𝑤 3 ℎ3
𝑄= ∫ [ℎ − ] 𝑑𝑥
2𝜇𝐿 0 3
𝑃𝑜 − 𝑃𝐿 𝑤 2 3
𝑄= ∫ ℎ 𝑑𝑥
2𝜇𝐿 0 3

𝑃𝑜 − 𝑃𝐿 2 3
𝑄= [ ℎ 𝑤]
2𝜇𝐿 3
(𝑃𝑜 − 𝑃𝐿)𝑤ℎ3
⟹ 4. 𝑄 =
3𝜇𝐿
𝑄
<𝑣 >= 𝑤 ℎ
∫0 ∫0 𝑑𝑦𝑑𝑥

(𝑃𝑜 − 𝑃𝐿)𝑤ℎ3
3𝜇𝐿
<𝑣 >=
𝑤ℎ
(𝑃𝑜 − 𝑃𝐿)ℎ2
⟹ 5. < 𝑣 > =
3𝜇𝐿
1.8
Condiciones
 Unidireccional
 Flujo laminar
 Estado estable
 ρ y µ constantes
 Isotérmico

𝜌𝑣𝑧𝑣𝑧∆𝑥∆𝑦|𝑧 − 𝜌𝑣𝑧𝑣𝑧∆𝑥∆𝑦|𝑧+∆𝑧 + 𝓉𝑦𝑧 ∆𝑥∆𝑧|𝑦 − 𝓉𝑦𝑧 ∆𝑥∆𝑧|𝑦+∆𝑦 = 0

𝓉𝑦𝑧 ∆𝑥∆𝑧|𝑦 − 𝓉𝑦𝑧 ∆𝑥∆𝑧|𝑦+∆𝑦
=0
∆𝑥∆𝑦∆𝑧
𝐿𝑖𝑚 𝓉𝑦𝑧 |𝑦 − 𝓉𝑦𝑧 |𝑦+∆𝑦
=0
∆𝑦 → 0 ∆𝑦
𝑑𝓉𝑦𝑧 𝑑𝑣𝑧
=0 → 𝓉𝑦𝑧 = ±𝜇
𝑑𝑦 𝑑𝑦
𝑑𝑣𝑧
𝓉𝑦𝑧 = −𝜇
𝑑𝑦
𝑑 𝑑𝑣𝑧
(−𝜇 )=0
𝑑𝑦 𝑑𝑦
𝑑 𝑑𝑣𝑧
( )=0
𝑑𝑦 𝑑𝑦
𝑑𝑣𝑧
∫𝑑( )=0
𝑑𝑦
𝑑𝑣𝑧
= 𝑐1
𝑑𝑦

∫ 𝑑𝑣𝑧 = 𝑐1 ∫ 𝑑𝑦

𝑣𝑧 = 𝑐1𝑦 + 𝑐2
𝐶. 𝐹. 1 𝑣𝑧 = 𝑣 𝑦 = 0
𝐶. 𝐹. 2 𝑣𝑧 = 0 𝑦 = ℎ
𝑣 = 𝑐2 0 = 𝑐1ℎ + 𝑐2
 𝑐2 𝑣
𝑐1 = − =−
ℎ ℎ
𝑣
𝑣𝑧 = − 𝑦 + 𝑣
ℎ
𝑦
⟹ 1. 𝑣𝑧 = 𝑣 [1 − ]
ℎ
𝑑𝑣𝑧
𝓉𝑦𝑧 = −𝜇
𝑑𝑦
𝑑 𝑦
𝓉𝑦𝑧 = −𝜇 (𝑣 (1 − ))
𝑑𝑦 ℎ
1
𝓉𝑦𝑧 = −𝜇𝑣 −
ℎ
𝜇𝑣
⟹ 2. 𝓉𝑦𝑧 =
ℎ
𝐶. 𝐹. 3 𝑣𝑧 = 𝑣𝑚á𝑥 𝑦=0
⟹ 3. 𝑣𝑚á𝑥 = 𝑣

𝑑𝑄 = 𝑣𝑧𝑑𝑦𝑑𝑥
𝑄 𝑤 ℎ
∫ 𝑑𝑄 = ∫ ∫ 𝑣𝑧𝑑𝑦𝑑𝑥
0 0 0
𝑤 ℎ
𝑦
𝑄 = ∫ ∫ 𝑣 [1 − ] 𝑑𝑦𝑑𝑥
0 0 ℎ
𝑤 ℎ
1 ℎ
𝑄 = ∫ [𝑣 ∫ 𝑑𝑦 − ∫ 𝑦𝑑𝑦] 𝑑𝑥
0 0 ℎ 0
𝑤
ℎ
𝑄 = ∫ 𝑣 [ℎ − ] 𝑑𝑥
0 2
𝑤
𝑣ℎ 𝑣ℎ
𝑄=∫ 𝑑𝑥 ∴ 𝑄 = 𝑤
0 2 2
𝑤𝑣ℎ
⟹ 4. 𝑄 =
2
𝑄 𝑤𝑣ℎ/2 𝑣
<𝑣 >= 𝑤 ℎ ∴ =
∫0 ∫0 𝑑𝑦𝑑𝑥 𝑤ℎ 2

𝑣
⟹ 5. < 𝑣 > =
2
1.9
Condiciones de frontera
 Unidireccional
 Flujo laminar
 ρ y µ constantes
 Estado estacionario
 Isotérmico

𝜌𝑣𝑧𝑣𝑧2𝜋𝑟∆𝑟|𝑧 − 𝜌𝑣𝑧𝑣𝑧2𝜋𝑟∆𝑟|𝑧+∆𝑧 + 𝓉𝑟𝑧 2𝜋𝑟𝐿|𝑟 − 𝓉𝑟𝑧 2𝜋𝑟𝐿|𝑟+∆𝑟 + 𝑃𝑜2𝜋𝑟∆𝑟|𝑧 − 𝑃𝐿2𝜋𝑟∆𝑟|𝑧−∆𝑧

+ 𝜌𝑔(2𝜋𝑟∆𝑟𝐿) = 0
𝓉𝑟𝑧 2𝜋𝑟𝐿|𝑟 − 𝓉𝑟𝑧 2𝜋𝑟𝐿|𝑟+∆𝑟 + 𝑃𝑜2𝜋𝑟∆𝑟|𝑧 − 𝑃𝐿2𝜋𝑟∆𝑟|𝑧−∆𝑧 + 𝜌𝑔(2𝜋𝑟∆𝑟𝐿) = 0
2𝜋∆𝑟𝐿
𝓉𝑟𝑧 𝑟|𝑟 − 𝓉𝑟𝑧 𝑟|𝑟+∆𝑟 𝑃𝑜𝑟|𝑧 − 𝑃𝐿𝑟|𝑧−∆𝑧
= + 𝜌𝑔𝑟
∆𝑟 𝐿
𝑑𝓉𝑟𝑧 𝑟 (𝑃𝑜 − 𝑃𝐿)𝑟 𝐿
= + 𝜌𝑔𝑟 ( )
𝑑𝑟 𝐿 𝐿
𝑑𝓉𝑟𝑧 𝑟 𝑟(𝑃𝑜 − 𝑃𝐿 + 𝜌𝑔𝐿)
=
𝑑𝑟 𝐿
𝒫𝑜 − 𝒫𝐿 = 𝑃𝑜 − 𝑃𝐿 + 𝜌𝑔𝐿
𝑑𝓉𝑟𝑧 𝑟 𝑟(𝒫𝑜 − 𝒫𝐿)
=
𝑑𝑟 𝐿
(𝒫𝑜 − 𝒫𝐿)
∫ 𝑑𝓉𝑟𝑧 𝑟 = ∫ 𝑟𝑑𝑟
𝐿
(𝒫𝑜 − 𝒫𝐿) 𝑟 2
𝓉𝑟𝑧 𝑟 = + 𝑐1
𝐿 2
𝐶. 𝐹. 1 𝓉𝑟𝑧 = 0 𝑟 = 0 ∴ 𝑐1 = 0
(𝒫𝑜 − 𝒫𝐿) 𝑟 2
𝓉𝑟𝑧 =
𝐿 2𝑟
 (𝒫𝑜 − 𝒫𝐿)𝑟
⟹ 1. 𝓉𝑟𝑧 =
2𝐿
(𝒫𝑜 − 𝒫𝐿)𝑟 𝑑𝑣𝑧
= −𝜇
2𝐿 𝑑𝑟
(𝒫𝑜 − 𝒫𝐿)
∫ 𝑑𝑣𝑧 = − ∫ 𝑟𝑑𝑟
2𝐿𝜇
(𝒫𝑜 − 𝒫𝐿)𝑟 2
𝑣𝑧 = − + 𝑐2
4𝐿𝜇
𝐶. 𝐹. 2 𝑣𝑧 = 0 𝑟 = 𝑅
(𝒫𝑜 − 𝒫𝐿)𝑅 2
𝑐2 =
2𝐿𝜇
(𝒫𝑜 − 𝒫𝐿)𝑟 2 (𝒫𝑜 − 𝒫𝐿)𝑅 2
𝑣𝑧 = − +
4𝐿𝜇 4𝐿𝜇
𝒫𝑜 − 𝒫𝐿 2
⟹ 2. 𝑣𝑧 = (𝑅 − 𝑟 2 )
4𝐿𝜇

𝐶. 𝐹. 3 𝑣𝑧 = 𝑣𝑧𝑚𝑎𝑥 𝑟 = 0
𝒫𝑜 − 𝒫𝐿 2
⟹ 3. 𝑣𝑧𝑚𝑎𝑥 = 𝑅
4𝐿𝜇

𝑄 = 𝑣𝑧𝐴
𝑑𝑄 = 𝑣𝑧𝑑𝐴
𝑑𝑄 = 𝑣𝑧𝑟𝑑𝑟𝑑𝜃
𝑄 2𝜋 𝑅
𝒫𝑜 − 𝒫𝐿 2
∫ 𝑑𝑄 = ∫ ∫ [𝑅 − 𝑟 2 ]𝑟𝑑𝑟𝑑𝜃
0 0 0 4𝐿𝜇

𝒫𝑜 − 𝒫𝐿 2𝜋 𝑅 2
𝑄= ∫ ∫ [𝑅 − 𝑟 2 ]𝑟𝑑𝑟𝑑𝜃
4𝐿𝜇 0 0

𝒫𝑜 − 𝒫𝐿 2𝜋 𝑅 4
𝑄= ∫ 𝑑𝜃
4𝐿𝜇 0 4

𝒫𝑜 − 𝒫𝐿 𝑅 4
𝑄= ( ) 2𝜋
4𝐿𝜇 4

𝜋(𝒫𝑜 − 𝒫𝐿)𝑅 4
⟹ 4. 𝑄 =
8𝐿𝜇

𝑄 =< 𝑣𝑧 > 𝐴
 𝑑𝑄 =< 𝑣𝑧 > 𝑑𝐴
𝑑𝑄 =< 𝑣𝑧 > 𝑟𝑑𝑟𝑑𝜃
𝑄
∫0 𝑑𝑄
𝑑𝑄 = 2𝜋 𝑅
∫0 ∫0 𝑟𝑑𝑟𝑑𝜃
𝑄
< 𝑣𝑧 >= 2𝜋 𝑅
∫0 ∫0 𝑟𝑑𝑟𝑑𝜃

𝜋(𝒫𝑜 − 𝒫𝐿)𝑅 4
8𝐿𝜇
< 𝑣𝑧 > =
𝜋𝑅 2
(𝒫𝑜 − 𝒫𝐿)𝑅 4
⟹ 5. < 𝑣𝑧 > =
8𝐿𝜇
1.10 Condiciones de operación.
 Unidireccional
 Flujo laminar
 ρ y µ constantes
 Estado estacionario
 Isotérmico

𝜌𝑣𝑧𝑣𝑧∆𝑥∆𝑦|𝑧 − 𝜌𝑣𝑧𝑣𝑧2𝜋∆𝑥∆𝑦|𝑧+∆𝑧 + 𝓉𝑟𝑧 2𝜋𝑟𝐿|𝑟 − 𝓉𝑟𝑧 2𝜋𝑟𝐿|𝑟+∆𝑟 = 0

𝓉𝑟𝑧 2𝜋𝑟𝐿|𝑟+∆𝑟 − 𝓉𝑟𝑧 2𝜋𝑟𝐿|𝑟
=0
2𝜋𝐿∆𝑟
𝓉𝑟𝑧 𝑟|𝑟+∆𝑟 − 𝓉𝑟𝑧 𝑟|𝑟
=0
∆𝑟
𝑙𝑖𝑚 𝓉𝑟𝑧 𝑟|𝑟+∆𝑟 − 𝓉𝑟𝑧 𝑟|𝑟
=0
∆𝑟 → 0 ∆𝑟
𝑑𝑣𝑧
𝓉𝑟𝑧 = ±𝜇
𝑑𝑟
𝑑 𝑑𝑣𝑧
(−𝜇 )=0
𝑑𝑟 𝑑𝑟
 𝑑 𝑑𝑣𝑧
(−𝑟 )=0
𝑑𝑟 𝑑𝑟
𝑑𝑣𝑧
𝑟 = 𝑐1
𝑑𝑟
𝑑𝑣𝑧 𝑐1
=
𝑑𝑟 𝑟
𝑐1
∫ 𝑑𝑣𝑧 = ∫ 𝑑𝑟
𝑟
𝑣𝑧 = 𝑐𝑧𝑙𝑛𝑟 + 𝑐2
𝐶. 𝐹. 1 𝑣𝑧 = 𝑣 𝑟 = 𝐾𝐹
𝐶. 𝐹. 2 𝑣𝑧 = 0 𝑟=𝑅
𝑣 = 𝑐1𝑙𝑛𝐾𝑅 + 𝑐2 ∴ 𝑣 = 𝑐1𝑙𝑛𝐾𝑅 − 𝑐1𝑙𝑛𝑅
0 = 𝑐1𝑙𝑛𝑅 + 𝑐2 ∴ 𝑐2 = −𝑐1𝑙𝑛𝑅
𝑣 = 𝑐1[𝑙𝑛𝐾𝑅 − 𝑙𝑛𝑅]
𝑣 𝑣
𝑐1 = ∴ 𝑐1 =
𝑘𝑅 𝑙𝑛𝑘
𝑙𝑛 𝑅
𝑣
𝑐2 = − 𝑙𝑛𝑅
𝑙𝑛𝑘
𝑣 𝑣
𝑣𝑧 = 𝑙𝑛𝑟 − 𝑙𝑛𝑅
𝑙𝑛𝑘 𝑙𝑛𝑘
𝑣
𝑣𝑧 = (𝑙𝑛𝑟 − 𝑙𝑛𝑅)
𝑙𝑛𝑘
𝑣 𝑟
⟹ 1. 𝑣𝑧 = (𝑙𝑛 )
𝑙𝑛𝑘 𝑅
1.11 Condiciones de operación
 Unidireccional
 Flujo laminar
 ρ y µ constantes
 Isotérmico
 Estado estacionario
 𝜌𝑣𝑧𝑣𝑧∆𝑥∆𝑦|𝑧 − 𝜌𝑣𝑧𝑣𝑧∆𝑥∆𝑦|𝑧+∆𝑧 + 𝓉𝑥𝑧 ∆𝑦∆𝑧|𝑥 − 𝓉𝑥𝑧 ∆𝑦∆𝑧|𝑥+∆𝑥 + 𝑃𝑜∆𝑥∆𝑦|𝑧
− 𝑃𝐿∆𝑥∆𝑦|𝑧+∆𝑧 + 𝜌𝑔∆𝑥∆𝑦∆𝑧 = 0
𝓉𝑥𝑧 ∆𝑦∆𝑧|𝑥 − 𝓉𝑥𝑧 ∆𝑦∆𝑧|𝑥+∆𝑥 + 𝑃𝑜∆𝑥∆𝑦|𝑧 − 𝑃𝐿∆𝑥∆𝑦|𝑧+∆𝑧 + 𝜌𝑔∆𝑥∆𝑦∆𝑧
=0
∆𝑥∆𝑦∆𝑧
𝓉𝑥𝑧 |𝑥 − 𝓉𝑥𝑧 |𝑥+∆𝑥 𝑃𝑜|𝑧 − 𝑃𝐿|𝑧+∆𝑧
+ + 𝜌𝑔 = 0
∆𝑥 ∆𝑧
𝑃𝑜 − 𝑃𝐿 −𝓉𝑥𝑧 |𝑥 + 𝓉𝑥𝑧 |𝑥+∆𝑥
∆𝑧 = 𝐿 ∴ + 𝜌𝑔 =
𝐿 ∆𝑥
𝑃𝑜 − 𝑃𝐿 𝑙𝑖𝑚 −𝓉𝑥𝑧 |𝑥 + 𝓉𝑥𝑧 |𝑥+∆𝑥
+ 𝜌𝑔 =
𝐿 ∆𝑥 → 0 ∆𝑥
𝑃𝑜 − 𝑃𝐿 𝑑𝓉𝑥𝑧
+ 𝜌𝑔 =
𝐿 𝑑𝑥
→ 𝒫𝑜 − 𝒫𝐿 = 𝑃𝑜 − 𝑃𝐿 + 𝜌𝑔
𝒫𝑜 − 𝒫𝐿 𝑑𝓉𝑥𝑧
=
𝐿 𝑑𝑥
𝒫𝑜 − 𝒫𝐿
∫ 𝑑𝓉𝑥𝑧 = ∫ 𝑑𝑥
𝐿
𝒫𝑜 − 𝒫𝐿
𝓉𝑥𝑧 = 𝑥 + 𝑐1
𝐿
𝐶. 𝐹. 1 𝓉𝑥𝑧 = 0 𝑥=0 𝑐1 = 0
𝒫𝑜 − 𝒫𝐿
⟹ 1. 𝓉𝑥𝑧 = 𝑥
𝐿
𝑑𝑣𝑧
𝓉𝑥𝑧 = ±𝜇
𝑑𝑥
𝑑𝑣𝑧 𝒫𝑜 − 𝒫𝐿
−𝜇 = 𝑥
𝑑𝑥 𝐿
𝒫𝑜 − 𝒫𝐿
∫ 𝑑𝑣𝑧 = − ∫ 𝑥𝑑𝑥
𝜇𝐿
𝒫𝑜 − 𝒫𝐿 2
𝑣𝑧 = − 𝑥 + 𝑐2
2𝜇𝐿
𝐶. 𝐹. 2 𝑣𝑧 = 0 𝑥=𝐵
𝒫𝑜 − 𝒫𝐿 2
𝑐2 = 𝐵
2𝜇𝐿
𝒫𝑜 − 𝒫𝐿 2 𝒫𝑜 − 𝒫𝐿 2
𝑆𝑢𝑠𝑡. 𝑣𝑧 = − 𝑥 + 𝐵
2𝜇𝐿 2𝜇𝐿
 𝒫𝑜 − 𝒫𝐿 2
⟹ 2. 𝑣𝑧 = (𝐵 − 𝑥 2 )
2𝜇𝐿

𝐶. 𝐹. 3 𝑣𝑧 = 𝑣𝑧𝑚𝑎𝑥 𝑥=0
𝒫𝑜 − 𝒫𝐿 2
⟹ 3. 𝑣𝑧𝑚𝑎𝑥 = (𝐵 )
2𝜇𝐿

𝑄 = 𝑣𝑧𝐴
𝑑𝑄 = 𝑣𝑧𝑑𝐴
𝑑𝑄 = 𝑣𝑧𝑑𝑥𝑑𝑦
𝑄 𝑤 𝐵
𝒫𝑜 − 𝒫𝐿 2
∫ 𝑑𝑄 = ∫ ∫ (𝐵 − 𝑥 2 )𝑑𝑥𝑑𝑦
0 0 0 2𝜇𝐿

𝑤 = 𝑎𝑛𝑐ℎ𝑜 𝑑𝑒 𝑡𝑢𝑏𝑒𝑟𝑖𝑎
𝒫𝑜 − 𝒫𝐿 𝑤 3 𝐵 3
𝑄= ∫ (𝐵 − ) 𝑑𝑦
2𝜇𝐿 0 3

𝒫𝑜 − 𝒫𝐿 𝑤 2𝐵 3
𝑄= ∫ 𝑑𝑦
2𝜇𝐿 0 3

𝐵 3 𝒫𝑜 − 𝒫𝐿 𝑤
𝑄= ∫ 𝑑𝑦
3𝜇𝐿 0

2𝐵 3 𝑤(𝒫𝑜 − 𝒫𝐿)
⟹ 4. 𝑄 =
3𝜇𝐿

𝑄 =< 𝑣 > 𝐴
𝑑𝑄 =< 𝑣 > 𝑑𝐴
𝑄
<𝑣 >= 𝑤 𝐵
∫0 ∫0 𝑑𝑥𝑑𝑦

2𝐵3 𝑤(𝒫𝑜 − 𝒫𝐿)

3𝜇𝐿
< 𝑣𝑧 > =
𝑤𝐵
2𝐵 2 (𝒫𝑜 − 𝒫𝐿)
⟹ 5. < 𝑣𝑧 > =
3𝜇𝐿

𝐵2 (𝒫𝑜 − 𝒫𝐿)
𝑣𝑧𝑚𝑎𝑥 =
2𝜇𝐿
1
< 𝑣𝑧 > = 𝑣𝑚𝑎𝑥
3
1.12 Condiciones de operación
 Unidireccional
 Flujo laminar
 Isotérmico
 ρ y µ constantes
 Estado estacionario

𝜌𝑣𝑧𝑣𝑧2𝜋𝑟∆𝑟|𝑧 − 𝜌𝑣𝑧𝑣𝑧2𝜋𝑟∆𝑟|𝑧+∆𝑧 + 𝓉𝑟𝑧 2𝜋𝑟∆𝑧|𝑟 − 𝓉𝑟𝑧 2𝜋𝑟∆𝑧|𝑟+∆𝑟 + 𝜌𝑔2𝜋𝑟∆𝑟∆𝑧 = 0

𝓉𝑟𝑧 2𝜋𝑟∆𝑧|𝑟 − 𝓉𝑟𝑧 2𝜋𝑟∆𝑧|𝑟+∆𝑟 + 𝜌𝑔2𝜋𝑟∆𝑟∆𝑧
=0
2𝜋∆𝑟∆𝑧
𝓉𝑟𝑧 𝑟|𝑟 − 𝓉𝑟𝑧 𝑟|𝑟+∆𝑟
+ 𝜌𝑔𝑟 = 0
∆𝑟
𝑙𝑖𝑚 −𝓉𝑟𝑧 𝑟|𝑟 + 𝓉𝑟𝑧 𝑟|𝑟+∆𝑟
= 𝜌𝑔𝑟
∆𝑟 → 0 ∆𝑟
𝑑𝓉𝑟𝑧 𝑟
= 𝜌𝑔𝑟
𝑑𝑟

∫ 𝑑𝓉𝑟𝑧 𝑟 = 𝜌𝑔 ∫ 𝑟𝑑𝑟

𝑟2
𝓉𝑟𝑧 𝑟 = 𝜌𝑔 + 𝑐1
2
𝐶. 𝐹. 1 𝓉𝑟𝑧 = 0 𝑟 = 𝑘𝑅
𝜌𝑔𝑘 2 𝑅 2
0= + 𝑐1
2
𝜌𝑔𝑘 2 𝑅 2
𝑐1 = −
2
𝑟 2 𝜌𝑔𝑘 2 𝑅 2
𝓉𝑟𝑧 𝑟 = 𝜌𝑔 −
2 2
𝜌𝑔𝑟 𝜌𝑔𝑘 2 𝑅 2
𝓉𝑟𝑧 = −
2 2𝑟
𝜌𝑔 𝑘 2𝑅2
⟹ 1. 𝓉𝑟𝑧 = [𝑟 − ]
2 𝑟
 𝑑𝑣𝑧
𝓉𝑟𝑧 = ±𝜇
𝑑𝑟
𝑑𝑣𝑧
𝓉𝑟𝑧 = −𝜇
𝑑𝑟
𝑑𝑣𝑧 𝜌𝑔 𝑘 2𝑅2
−𝜇 = [𝑟 − ]
𝑑𝑟 2 𝑟

𝜌𝑔 𝑘 2𝑅2
𝑑𝑣𝑧 = − [𝑟 − ] 𝑑𝑟
2𝜇 𝑟

𝜌𝑔 𝑘 2𝑅2
∫ 𝑑𝑣𝑧 = − ∫ [𝑟 − ] 𝑑𝑟
2𝜇 𝑟

𝜌𝑔 𝑟 2
𝑣𝑧 = − [ − 𝑘 2 𝑅 2 𝑙𝑛𝑟] + 𝑐2
2𝜇 2

𝐶. 𝐹. 2 𝑣𝑧 = 0 𝑟=𝑅
𝜌𝑔 𝑅 2
0=− [ − 𝑘 2 𝑅 2 𝑙𝑛𝑅] + 𝑐2
2𝜇 2

𝜌𝑔 𝑅 2
𝑐2 = [ − 𝑘 2 𝑅 2 𝑙𝑛𝑅]
2𝜇 2

𝜌𝑔 𝑟 2 𝜌𝑔 𝑅 2
𝑣𝑧 = − [ − 𝑘 2 𝑅 2 𝑙𝑛𝑟] + [ − 𝑘 2 𝑅 2 𝑙𝑛𝑅]
2𝜇 2 2𝜇 2

𝜌𝑔 𝑟 2 𝑅2
𝑣𝑧 = [− + 𝑘 2 𝑅 2 𝑙𝑛𝑟 + − 𝑘 2 𝑅 2 𝑙𝑛𝑅]
2𝜇 2 2

𝜌𝑔 𝑅 2 𝑟 2 2 2 2 𝑟
𝑣𝑧 = [ − + 𝑘 𝑅 𝑙𝑛 ]
2𝜇 2 2 2 𝑅
𝜌𝑔 2 𝑟
⟹ 2. 𝑣𝑧 = [𝑅 − 𝑟 2 + 𝑘 2 𝑅 2 𝑙𝑛 ]
4𝜇 𝑅
1.13 Condiciones de operación
 Unidireccional
 Flujo laminar
 Isotérmico
 ρ y µ contantes
 Estado estable
𝜌𝑣𝑧𝑣𝑧∆𝑦∆𝑥|𝑧 − 𝜌𝑣𝑧𝑣𝑧∆𝑦∆𝑥|𝑧+∆𝑧 + 𝓉𝑦𝑧 ∆𝑥∆𝑧|𝑦+∆𝑦 − 𝓉𝑦𝑧 ∆𝑥∆𝑧|𝑦 + 𝜌𝑔∆𝑥∆𝑦∆𝑧𝑐𝑜𝑠𝜃 = 0
𝓉𝑦𝑧 ∆𝑥∆𝑧|𝑦+∆𝑦 − 𝓉𝑦𝑧 ∆𝑥∆𝑧|𝑦 + 𝜌𝑔∆𝑥∆𝑦∆𝑧𝑐𝑜𝑠𝜃
=0
∆𝑥∆𝑦∆𝑧
𝓉𝑦𝑧 ∆𝑥∆𝑧|𝑦+∆𝑦 − 𝓉𝑦𝑧 ∆𝑥∆𝑧|𝑦
+ 𝜌𝑔𝑐𝑜𝑠𝜃 = 0
∆𝑦

𝑙𝑖𝑚 𝓉𝑦𝑧 ∆𝑥∆𝑧|𝑦+∆𝑦 − 𝓉𝑦𝑧 ∆𝑥∆𝑧|𝑦

= −𝜌𝑔𝑐𝑜𝑠𝜃
∆𝑟 → 0 ∆𝑦
𝑑𝓉𝑦𝑧
= −𝜌𝑔𝑐𝑜𝑠𝜃
𝑑𝑦

∫ 𝑑 𝓉𝑦𝑧 = −𝜌𝑔𝑐𝑜𝑠𝜃 ∫ 𝑑 𝑦

𝓉𝑦𝑧 = −𝜌𝑔𝑐𝑜𝑠𝜃𝑦 + 𝑐1

𝐶. 𝐹. 1 𝓉𝑦𝑧 = 0 𝑦=𝛿

𝑐1 = 𝜌𝑔𝑐𝑜𝑠𝜃𝑦𝛿
⟹ 1. 𝓉𝑦𝑧 = 𝜌𝑔𝑐𝑜𝑠𝜃(𝛿 − 𝑦)

𝑑𝑣𝑧
𝓉𝑦𝑧 = ±𝜇
𝑑𝑦
𝑑𝑣𝑧
𝓉𝑦𝑧 = 𝜇
𝑑𝑦
𝑑𝑣𝑧
𝜇 = 𝜌𝑔𝑐𝑜𝑠𝜃(𝛿 − 𝑦)
𝑑𝑦
𝜌𝑔𝑐𝑜𝑠𝜃
∫ 𝑑𝑣𝑧 = ∫(𝛿 − 𝑦)𝑑𝑦
𝜇
𝜌𝑔𝑐𝑜𝑠𝜃 𝑦2
𝑣𝑧 = [𝛿𝑦 − ] + 𝑐2
𝜇 2

𝐶. 𝐹. 2 𝑣𝑧 = 0 𝑦 = 0 ∴ 𝑐2 = 0
𝜌𝑔𝑐𝑜𝑠𝜃 𝑦2
⟹ 2. 𝑣𝑧 = [𝛿𝑦 − ]
𝜇 2

𝐶. 𝐹. 3 𝑣𝑧 = 𝑣𝑧𝑚𝑎𝑥 𝑦=𝛿
𝜌𝑔𝑐𝑜𝑠𝜃 2 𝛿 2
𝑣𝑧𝑚𝑎𝑥 = [𝛿 − ]
𝜇 2
𝜌𝑔𝑐𝑜𝑠𝜃 2
⟹ 3. 𝑣𝑧𝑚𝑎𝑥 = 𝛿
2𝜇
 𝑑𝑄 = 𝑣𝑧𝑑𝐴
𝑤 = 𝑎𝑛𝑐ℎ𝑜
𝑄 𝑤 𝛿
𝜌𝑔𝑐𝑜𝑠𝜃 𝑦2
∫ 𝑑𝑄 = ∫ ∫ [𝛿𝑦 − ] 𝑑𝑦𝑑𝑥
0 0 0 𝜇 2

𝜌𝑔𝑐𝑜𝑠𝜃 𝑤 𝛿 3 𝛿 3
𝑄= ∫ [ − ] 𝑑𝑥
𝜇 0 2 6
𝜌𝑔𝑐𝑜𝑠𝜃 3
⟹ 4. 𝑄 = 𝛿 𝑤
3𝜇
𝑄
<𝑣 >= 𝑤 𝛿
∫0 ∫0 𝑑𝑦𝑑𝑥
𝜌𝑔𝑐𝑜𝑠𝜃 3
3𝜇 𝛿 𝑤
<𝑣 >=
𝑤𝛿
𝜌𝑔𝑐𝑜𝑠𝜃𝛿 2
⟹ 5. < 𝑣 > =
3𝜇
UNIDAD
DOS
EJERCICIOS DE CLASE
Ecuaciones de variación para sistemas isotérmicos.
2.1 Condiciones de operación:
Estado estable
Unidireccional
ρ y µ constantes
Isotérmico
Flujo laminar

Tabla 3.4-1
𝜕𝜌 𝜕 𝜕 𝜕
+ (𝜌𝑣𝑥) + (𝜌𝑣𝑦) + (𝜌𝑣𝑧) = 0
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝑣𝑧 = 𝑓(𝑦, 𝑧)
𝑣𝑥 = 0
𝑣𝑦 = 0
𝑑𝑣𝑧
=0
𝑑𝑧
Tabla 3.4-2
Componente x: 0 = 𝜌𝑔𝑥
Componente y: 0 = 𝜌𝑔𝑦
𝜕𝜏
Componente z: 0 = − 𝜕𝑦 + 𝜌𝑔𝑧 𝑐𝑜𝑠𝜃
𝑑𝜏𝑦𝑧
= 𝜌𝑔𝑧 𝑐𝑜𝑠𝜃
𝑑𝑦
∫ 𝑑𝜏𝑦𝑧 = 𝜌𝑔𝑧 𝑐𝑜𝑠𝜃 ∫ 𝑑𝑦
𝜏𝑦𝑧 = 𝜌𝑔𝑧 𝑐𝑜𝑠𝜃𝑦 + 𝑐1
𝐶. 𝐹. 1 𝜏𝑦𝑧 = 0 𝑦 = 0 𝑐1 = 0

 1. 𝜏𝑦𝑧 = 𝜌𝑔𝑧 𝑐𝑜𝑠𝜃𝑦

𝑑𝑣𝑧
𝜏𝑦𝑧 = ± 𝜇
𝑑𝑦
𝑑𝑣𝑧
𝜏𝑦𝑧 = − 𝜇
𝑑𝑦
𝑑𝑣𝑧
−𝜇 = 𝜌𝑔𝑧 𝑐𝑜𝑠𝜃𝑦
𝑑𝑦
𝜌𝑔𝑧 𝑐𝑜𝑠𝜃
∫ 𝑑𝑣𝑧 = − ∫ 𝑦 𝑑𝑦
𝜇
𝜌𝑔𝑧 𝑐𝑜𝑠𝜃 2
𝑣𝑧 = − 𝑦 + 𝑐2
2𝜇
𝐶. 𝐹. 2 𝑣𝑧 = 0 𝑦 = 𝛿

[27]
 𝜌𝑔𝑧 𝑐𝑜𝑠𝜃 2
𝑐2 = 𝛿
2𝜇
𝜌𝑔𝑧 𝑐𝑜𝑠𝜃 2 𝜌𝑔𝑧 𝑐𝑜𝑠𝜃 2
𝑣𝑧 = − 𝑦 + 𝛿
2𝜇 2𝜇
𝜌𝑔𝑧 𝑐𝑜𝑠𝜃 2
 2. 𝑣𝑧 = [𝛿 − 𝑦 2 ]
2𝜇
𝐶. 𝐹. 3 𝑣𝑧 = 𝑣𝑧𝑚𝑎𝑥 𝑦 = 0
𝜌𝑔𝑧 𝑐𝑜𝑠𝜃 2
 3. 𝑣𝑧𝑚𝑎𝑥 = [𝛿 ]
2𝜇
𝑄 𝑊 𝛿
∫ 𝑑𝑄 = ∫ ∫ 𝑣𝑧𝑑𝑦𝑑𝑥
0 0 0
𝑊 𝛿
𝜌𝑔𝑧 𝑐𝑜𝑠𝜃 2
𝑄=∫ ∫ [𝛿 − 𝑦 2 ]𝑑𝑦𝑑𝑥
0 0 2𝜇
𝜌𝑔𝑧 𝑐𝑜𝑠𝜃
 4. 𝑄 = [𝑊𝛿 3 ]
3𝜇
𝑄
<𝑣 >= 𝑊 𝛿
∫0 ∫0 𝑑𝑦𝑑𝑥
𝜌𝑔𝑧 𝑐𝑜𝑠𝜃
[𝑊𝛿 3 ]
3𝜇
<𝑣 >= 𝑊 𝛿
∫0 ∫0 𝑑𝑦𝑑𝑥
𝜌𝑔𝑧 𝑐𝑜𝑠𝜃𝛿 2
 5. < 𝑣 > =
3𝜇

2.2 Condiciones de operación.

Estado estacionario
Unidireccional
Isotérmico
ρ y µ constantes
Flujo laminar

Tabla 3.4-1
𝑣𝑧 = 𝑓(𝑦, 𝑧)
𝑣𝑥 = 0
𝑣𝑦 = 0
𝑑𝑣𝑧
=0
𝑑𝑧
Tabla 3.4-2
Componente x: 0 = 𝜌𝑔𝑥
Componente y: 0 = 𝜌𝑔𝑦
𝑑2 𝑣𝑧
Componente z: 0=𝜇 𝑑𝑦 2

[28]
 𝑑2 𝑣𝑧
=0
𝑑𝑦 2
𝑑 𝑑𝑣𝑧
∫ ( )=0
𝑑𝑦 𝑑𝑦
𝑑𝑣𝑧
= 𝑐1
𝑑𝑦
𝑣𝑧 = 𝑐1𝑦 + 𝑐2 → 𝑒𝑐. 𝑔𝑒𝑛𝑒𝑟𝑎𝑙
𝐶. 𝐹. 1 𝑣𝑧 = 𝑣 𝑦 = 0
𝐶. 𝐹. 2 𝑣𝑧 = 0 𝑦 = ℎ
𝑣 = 𝑐1(0) + 𝑐2
𝑣
𝑐2 = −
ℎ
𝑣
𝑣𝑧 = − 𝑦 + 𝑣
ℎ
𝑦
 1. 𝑣𝑧 = 𝑣 (1 − )
ℎ
𝑑𝑣𝑧
𝜏𝑦𝑧 = ± 𝜇
𝑑𝑦
𝑑 𝑦
𝜏𝑦𝑧 = − 𝜇 [𝑣 (1 − )]
𝑑𝑦 ℎ
𝜇𝑣
 2. 𝜏𝑦𝑧 =
ℎ
𝐶. 𝐹. 3 𝑣𝑧 = 𝑣𝑧𝑚𝑎𝑥 𝑦 = 0
 3. 𝑣𝑧𝑚𝑎𝑥 = 𝑣
𝑄 𝑊 𝛿
∫ 𝑑𝑄 = ∫ ∫ 𝑣𝑧𝑑𝑦𝑑𝑥
0 0 0
𝑊 𝛿
𝑦
𝑄 = ∫ ∫ [𝑣 (1 − )] 𝑑𝑦𝑑𝑥
0 0 ℎ
𝑣𝑤ℎ
 4. 𝑄 =
2
𝑄
<𝑣 >= 𝑊 𝛿
∫0 ∫0 𝑑𝑦𝑑𝑥
𝑣𝑤ℎ
<𝑣 >= 2
𝑤ℎ
𝑣
 5. < 𝑣 > =
2

[29]
2.3 Condiciones de operación.
Estado estacionario
Unidireccional
Isotérmico
ρ y µ constantes
Flujo laminar

Tabla 3.4-1
𝑣𝑧 = 𝑓(𝑟, 𝑧)
𝑣𝑟 = 0
𝑣𝜃 = 0
𝑑𝑣𝑧
=0
𝑑𝑧
Tabla 3.4-3
Componente r: 0 = 𝜌𝑔𝑟
Componente Ѳ: 0 = 𝜌𝑔𝜃
1 𝑑 𝑑𝑣𝑧
Componente z: 0 = 𝜇 [𝑟 𝑑𝑟 (𝑟 )]
𝑑𝑟

𝑑 𝑑𝑣𝑧
(𝑟 )=0
𝑑𝑟 𝑑𝑟
𝑑𝑣𝑧
𝑟 = 𝑐1
𝑑𝑟
𝑐1
𝑑𝑣𝑧 = 𝑑𝑟
𝑟
𝑣𝑧 = 𝑐1𝑙𝑛𝑟 + 𝑐2
𝐶. 𝐹. 1 𝑣𝑧 = 𝑣 𝑟 = 𝑘𝑅
𝐶. 𝐹. 2 𝑣𝑧 = 0 𝑟=𝑅
𝑣 = 𝑐1𝑙𝑛𝑟 + 𝑐2 … 𝑠𝑢𝑠𝑡𝑖𝑡𝑢𝑖𝑟
𝑣 = 𝑐1𝑙𝑛𝑘𝑅 + 𝑐2
0 = 𝑐1𝑙𝑛𝑅 + 𝑐2
𝑐2 = −𝑐1𝑙𝑛𝑅
𝑣 = 𝑐1𝑙𝑛𝑘𝑅 − 𝑐1𝑙𝑛𝑅
𝑣 = 𝑐1(𝑙𝑛𝑘𝑅 − 𝑙𝑛𝑅)
𝑣
𝑐1 =
𝑙𝑛𝐾
𝑣
𝑐2 = − 𝑙𝑛𝑅
𝑙𝑛𝐾

[30]
 𝑣 𝑣
𝑣𝑧 = 𝑙𝑛𝑟 − 𝑙𝑛𝑅
𝑙𝑛𝐾 𝑙𝑛𝐾
𝑣
𝑣𝑧 = [𝑙𝑛𝑟 − 𝑙𝑛𝑅]
𝑙𝑛𝐾
𝑣 𝑟
 1. 𝑣𝑧 = 𝑙𝑛
𝑙𝑛𝐾 𝑅
2.4 Condiciones de operación.
Estado estacionario
Unidireccional
Isotérmico
ρ y µ constantes
Flujo laminar

Tabla 3.4-1
𝑣𝑧 = 𝑓(𝑟, 𝑧)
𝑣𝑟 = 0
𝑣𝜃 = 0
𝑑𝑣𝑧
=0
𝑑𝑧
Tabla 3.4-3
Componente r: 0 = 𝜌𝑔𝑟
Componente Ѳ: 0 = 𝜌𝑔𝜃
𝑑𝑝 1 𝑑
Componente z: 0 = − 𝑑𝑧 − [𝑟 𝑑𝑟 (𝑟𝜏𝑟𝑧 )] + 𝜌𝑔𝑧
∆𝑃 1 𝑑
0=− − (𝑟𝜏𝑟𝑧 ) + 𝜌𝑔𝑧 ∆𝑧 − 𝐿
∆𝑧 𝑟 𝑑𝑟
𝑃𝑙 − 𝑃𝑜 1 𝑑
0=− − (𝑟𝜏𝑟𝑧 ) + 𝜌𝑔𝑧
𝐿 𝑟 𝑑𝑟
1𝑑 𝑃𝑜 − 𝑃𝑙
(𝑟𝜏𝑟𝑧 ) = + 𝜌𝑔𝑧
𝑟 𝑑𝑟 𝐿
1𝑑 𝑃𝑜 − 𝑃𝑙 + 𝜌𝑔𝑧 𝐿
(𝑟𝜏𝑟𝑧 ) =
𝑟 𝑑𝑟 𝐿
𝒫𝑜 − 𝒫𝑙 = 𝑃𝑜 − 𝑃𝑙 + 𝜌𝑔𝑧 𝐿
1𝑑 𝒫𝑜 − 𝒫𝑙
(𝑟𝜏𝑟𝑧 ) =
𝑟 𝑑𝑟 𝐿
𝑑 𝒫𝑜 − 𝒫𝑙
(𝑟𝜏𝑟𝑧 ) = ( )𝑟
𝑑𝑟 𝐿
Integrando…
𝒫𝑜 − 𝒫𝑙 2
𝑟𝜏𝑟𝑧 = ( ) 𝑟 + 𝑐1
2𝐿
𝐶. 𝐹. 1 𝜏𝑟𝑧 = 0 𝑟=0 𝑐1 = 0
[31]
 𝒫𝑜 − 𝒫𝑙 2
𝑟𝜏𝑟𝑧 = ( )𝑟
2𝐿
𝒫𝑜 − 𝒫𝑙
 2. 𝜏𝑟𝑧 = ( )𝑟
2𝐿
𝑑𝑣𝑧
𝜏𝑟𝑧 = −𝜇
𝑑𝑟
𝑑𝑣𝑧 𝒫𝑜 − 𝒫𝑙
−𝜇 = 𝑟
𝑑𝑟 2𝐿
𝒫𝑜 − 𝒫𝑙 2
𝑣𝑧 = − 𝑟 + 𝑐2
4𝜇𝐿
𝐶. 𝐹. 2 𝑣𝑧 = 0 𝑟 = 𝑅
𝒫𝑜 − 𝒫𝑙 2
𝑐2 = − 𝑅
4𝐿𝜇
𝒫𝑜 − 𝒫𝑙 2 𝒫𝑜 − 𝒫𝑙 2
𝑣𝑧 = − 𝑟 + 𝑅
4𝜇𝐿 2𝐿
𝒫𝑜 − 𝒫𝑙 2
 3. 𝑣𝑧 = [𝑅 − 𝑟 2 ]
4𝜇𝐿
𝐶. 𝐹. 3 𝑣𝑧 = 𝑣𝑧𝑚𝑎𝑥 𝑟 = 0
𝒫𝑜 − 𝒫𝑙 2
 4. 𝑣𝑧𝑚𝑎𝑥 = [𝑅 ]
4𝜇𝐿
𝑄 2𝜋 𝑅
∫ 𝑑𝑄 = ∫ ∫ 𝑣𝑧𝑟𝑑𝑟𝑑𝜃
0 0 0
2𝜋 𝑅
𝒫𝑜 − 𝒫𝑙 2
𝑄=∫ ∫ [𝑅 − 𝑟 2 ]𝑑𝑟𝑑𝜃
0 0 4𝜇𝐿
𝜋(𝒫𝑜 − 𝒫𝑙)𝑅 4
 4. 𝑄 =
8𝜇𝐿
𝑄
< 𝑣 > = 2𝜋 𝑅
∫0 ∫0 𝑟𝑑𝑟𝑑𝜃
𝜋(𝒫𝑜 − 𝒫𝑙)𝑅 4
8𝜇𝐿
<𝑣 >=
𝜋𝑅 2
(𝒫𝑜 − 𝒫𝑙)𝑅 2
 5. < 𝑣 > =
8𝜇𝐿
2.5 Condiciones de operación.
Estado estacionario
Unidireccional
Isotérmico
ρ y µ constantes
Flujo laminar
Tabla 3.4-1
𝑣𝑧 = 𝑓(𝑦, 𝑧)
𝑣𝑥 = 0
𝑣𝑦 = 0
[32]
 𝑑𝑣𝑧
=0
𝑑𝑧
Tabla 3.4-2 (D,E,F)
Componente x: 0 = 𝜌𝑔𝑥
Componente y: 0 = 𝜌𝑔𝑦
𝑑𝑃 𝑑2 𝑣𝑧
Componente z: 0 = − 𝑑𝑧 + 𝜇 − 𝜌𝑔𝑧 𝑐𝑜𝑠𝜃
𝑑𝑥 2
𝑑𝑃 𝑑 2 𝑣𝑧
+ 𝜌𝑔𝑧 𝑐𝑜𝑠𝜃 = 𝜇
𝑑𝑧 𝑑𝑥 2
2
∆𝑃 𝑑 𝑣𝑧
+ 𝜌𝑔𝑧 𝑐𝑜𝑠𝜃 = 𝜇 ∆𝑧 − 𝐿
∆𝑧 𝑑𝑥 2
𝑃𝑙 − 𝑃𝑜 𝑑 2 𝑣𝑧
+ 𝜌𝑔𝑧 𝑐𝑜𝑠𝜃 = 𝜇
𝐿 𝑑𝑥 2
𝑃𝑜 − 𝑃𝑙 𝑑 2 𝑣𝑧
− + 𝜌𝑔𝑧 𝑐𝑜𝑠𝜃 = 𝜇
𝐿 𝑑𝑥 2
𝑃𝑜 − 𝑃𝑙 𝑑 2 𝑣𝑧
−( − 𝜌𝑔𝑧 𝑐𝑜𝑠𝜃) = 𝜇
𝐿 𝑑𝑥 2
𝑃𝑜 − 𝑃𝑙 − 𝜌𝑔𝑧 𝑐𝑜𝑠𝜃𝐿 𝑑 2 𝑣𝑧
−( )=𝜇
𝐿 𝑑𝑥 2
𝒫𝑜 − 𝒫𝑙 = 𝑃𝑜 − 𝑃𝑙 − 𝜌𝑔𝑧 𝐿
𝒫𝑜 − 𝒫𝑙 𝑑 2 𝑣𝑧
−( )=𝜇
𝐿 𝑑𝑥 2
𝒫𝑜 − 𝒫𝑙 𝑑 2 𝑣𝑧
( ) = −𝜇
𝐿 𝑑𝑥 2
𝒫𝑜 − 𝒫𝑙 𝑑 𝑑 2 𝑣𝑧
( )= (−𝜇 )
𝐿 𝑑𝑥 𝑑𝑥 2
𝑑𝑣𝑧
𝜏𝑥𝑧 = −𝜇
𝑑𝑥
𝒫𝑜 − 𝒫𝑙 𝑑
= 𝜏
𝐿 𝑑𝑥 𝑥𝑧
𝒫𝑜 − 𝒫𝑙
𝜏𝑥𝑧 = 𝑥 + 𝑐1
𝐿
𝐶. 𝐹. 1 𝜏𝑥𝑧 = 0 𝑣 = 0 𝑐1 = 0
𝒫𝑜 − 𝒫𝑙
 1. 𝜏𝑥𝑧 = 𝑥
𝐿
𝑑𝑣𝑧 𝒫𝑜 − 𝒫𝑙
−𝜇 = 𝑥
𝑑𝑥 𝐿
𝒫𝑜 − 𝒫𝑙
𝑑𝑣𝑧 = − 𝑥𝑑𝑥
𝜇𝐿
𝒫𝑜 − 𝒫𝑙 𝑥 2
𝑣𝑧 = − + 𝑐2
𝜇𝐿 2
𝐶. 𝐹. 2 𝑣𝑧 = 0 𝑥 = 𝐵
𝒫𝑜 − 𝒫𝑙 𝐵 2
𝑐2 =
𝜇𝐿 2
[33]
 𝒫𝑜 − 𝒫𝑙 𝑥2 𝒫𝑜 − 𝒫𝑙 𝐵 2
𝑣𝑧 = − +
𝜇𝐿 2 𝜇𝐿 2

𝒫𝑜 − 𝒫𝑙 2
 2. 𝑣𝑧 = [𝐵 − 𝑥 2 ]
2𝜇𝐿
𝐶. 𝐹. 3. 𝑣𝑧 = 𝑣𝑧𝑚𝑎𝑥 𝑥 = 0
𝒫𝑜 − 𝒫𝑙 2
 3. 𝑣𝑧𝑚𝑎𝑥 = [𝐵 ]
2𝜇𝐿

𝑑𝑄 = 𝑣𝑧𝑑𝐴
𝑄 𝑊 𝐵
∫ 𝑑𝑄 = ∫ ∫ 𝑣𝑧𝑑𝑥𝑑𝑦
0 0 0
𝑊 𝐵
𝒫𝑜 − 𝒫𝑙 2
𝑄=∫ ∫ [𝐵 − 𝑥 2 ]𝑑𝑦𝑑𝑥
0 0 2𝜇𝐿
2(𝒫𝑜 − 𝒫𝑙)𝐵3 𝑤
 4. 𝑄 =
3𝜇𝐿
𝑄
<𝑣 >= 𝑊 𝐵
∫0 ∫0 𝑑𝑥𝑑𝑦
2(𝒫𝑜 − 𝒫𝑙)𝐵3 𝑤
3𝜇𝐿
<𝑣 >=
𝑤𝐵
(𝒫𝑜 − 𝒫𝑙)𝐵 2
 5. < 𝑣 > =
3𝜇𝐿

2.6 Condiciones de operación.

Estado estacionario
Unidireccional
Isotérmico
ρ y µ constantes
Flujo laminar

Tabla 3.4-1
𝑣𝜃 = 𝑓(𝑟, 𝜃)
𝑣𝑟 = 0
𝑣𝑧 = 0
𝑑𝑣𝜃
=0
𝑑𝜃
Tabla 3.4-3
𝑑𝑣 2 𝜃 𝑑𝑃
Componente r: −𝜌 = − 𝑑𝑟
𝑟
𝑑
Componente Ѳ: 0 = 𝜇𝑑
𝑑𝑝
Componente z: 0 = − 𝑑𝑧 + 𝜌𝑔𝑧

[34]
 𝑑 1𝑑
( 𝑟𝑣𝜃) = 0
𝑑𝑟 𝑟 𝑑𝑟
1 𝑑
∫𝑑( 𝑟𝑣𝜃) = 0
𝑟 𝑑𝑟
1𝑑
𝑟𝑣𝜃 = 𝑐1
𝑟 𝑑𝑟
𝑑
𝑟𝑑𝜃 = 𝑐1𝑟
𝑑𝑟
𝑑(𝑟𝑣𝜃) = 𝑐1𝑟𝑑𝑟
𝑐1𝑟 2
𝑟𝑣𝜃 = + 𝑐2
2
𝑐1𝑟 𝑐2
→ 𝑣𝜃 = +
2 𝑟
𝐶. 𝐹. 1 𝑣𝜃 = 0 𝑟 = 𝑘𝑅
𝐶. 𝐹. 2 𝑣𝜃 = Ω𝑜𝑅 𝑟 = 𝑅
𝑐1𝑘𝑅 𝑐2
(0 = + )−𝑘
2 𝑘𝑅
𝑐1𝑅 𝑐2
(Ω𝑜𝑅 = + )
2 𝑅
2
𝑐1𝑘 𝑅 𝑐2
0=− +
2 𝑘𝑅
𝑐1𝑅 𝑐2
Ω𝑜𝑅 = +
2 𝑅
___________________________________________
𝑐1𝑅 𝑐1𝑘 2 𝑅
Ω𝑜𝑅 = −
2 2
2Ω𝑜 = 𝑐𝑧(1 − 𝑘 2 )
2Ω𝑜
𝑐1 =
1 − 𝑘2
𝑐1𝑘𝑅 𝑐2
0= +
2 𝑘𝑅
2Ω𝑜 2 2
𝑐1𝑘 2 𝑅 2 1 − 𝑘 2 𝑘 𝑅
𝑐2 = − =
2 2
Ω𝑜𝑘 2 𝑅 2
𝑐2 = −
1 − 𝑘2
𝑐1𝑟 𝑐2
−−→ 𝑣𝜃 = +
2 𝑟
2Ω𝑜 Ω𝑜𝑘 2 𝑅 2
2𝑟 − 2
𝑣𝜃 = 1 − 𝑘 + 1 − 𝑘
2 𝑟
Ω𝑜 Ω𝑜𝑘 2 𝑅 2
𝑣𝜃 = 𝑟−
1 − 𝑘2 (1 − 𝑘 2 )𝑟
Ω𝑜 𝑘 2𝑅2
 1. 𝑣𝜃 = [𝑟 − ]
1 − 𝑘2 𝑟

[35]
 Tabla 3.4-6
𝑑 𝑣𝜃 𝑑 𝑣𝜃
𝜏𝑟𝜃 = −𝜇 [𝑟 ] = −𝜇𝑟
𝑑𝑟 𝑟 𝑑𝑟 𝑟
Ω𝑜 𝑘 2𝑅2
𝑑 1 − 𝑘 2 [𝑟 − 𝑟 ]
𝜏𝑟𝜃 = −𝑟 [ ]
𝑑𝑟 𝑟

𝑘 2𝑅2
𝜇𝑟Ω𝑜 𝑑 𝑟 ]
𝜏𝑟𝜃 =− [𝑟 −
1 − 𝑘 2 𝑑𝑟 𝑟

𝜇𝑟Ω𝑜 𝑑 𝑘 2𝑅2
𝜏𝑟𝜃 = − [1 − ]
1 − 𝑘 2 𝑑𝑟 𝑟2
𝜇𝑟Ω𝑜 2𝑘 2 𝑅 2
𝜏𝑟𝜃 = − [ ]
1 − 𝑘2 𝑟3
𝜇Ω𝑜 2𝑘 2 𝑅 2
 2. 𝜏𝑟𝜃 = − [ ]
1 − 𝑘2 𝑟2
𝜉 = 𝜏𝑟𝜃 × 𝐴 × 𝐵𝑟𝑎𝑧𝑜 𝑑𝑒 𝑝𝑎𝑙𝑎𝑛𝑐𝑎
𝜇𝑟Ω𝑜 2𝑘 2 𝑅 2
𝜉= [ ] × 2𝜋𝑅𝐿 × 𝑅
1 − 𝑘2 𝑟3

2.7 Condiciones de operación.

Estado estacionario
Unidireccional
Isotérmico
ρ y µ constantes
Flujo laminar

Tabla 3.4-1
𝑣𝜃 = 𝑓(𝑟, 𝜃)
𝑣𝑟 = 0
𝑣𝑧 = 0
𝑑𝑣𝜃
=0
𝑑𝜃
Tabla 3.4-3
𝑑𝑣 2 𝜃 𝑑𝑃
Componente r: −𝜌 =
𝑟 𝑑𝑟
𝑑 1 𝑑
Componente Ѳ: 0 = 𝜇 𝑑𝑟 𝑟 𝑑𝑟 (𝑟𝑣𝜃)
𝑑𝑝
Componente z: 0 = − 𝑑𝑧 + 𝜌𝑔𝑧
𝑑 1𝑑
( (𝑟𝑣𝜃)) = 0
𝑑𝑟 𝑟 𝑑𝑟
𝑑
(𝑟𝑣𝜃) = 𝑐1𝑟
𝑑𝑟
[36]
 𝑐1𝑟 2
𝑟𝑣𝜃 = + 𝑐2
2
𝑐1𝑟 𝑐2
−−→ 𝑣𝜃 = +
2 𝑟
𝐶. 𝐹. 1 𝑣𝜃 = Ω𝑜𝑘𝑅 𝑟 = 𝑘𝑅
𝐶. 𝐹. 1 𝑣𝜃 = 0 𝑟=𝑅
𝑐1𝑘𝑅 𝑐2
(Ω𝑜𝑘𝑅 = + )−𝑘
2 𝑘𝑅
𝑐1𝑅 𝑐2
0= +
2 𝑅
2
𝑐1𝑘 𝑅 𝑐2
−Ω𝑜𝑘 2 𝑅 = − +
2 𝑅
𝑐1𝑅 𝑐2
0= +
2 𝑅
________________________________________________
𝑐1𝑘 2 𝑅 𝑐1𝑅
−Ω𝑜𝑘 2 𝑅 = − +
2 2
2 2
−2Ω𝑜𝑘 𝑅 = −𝑐1𝑘 𝑅 + 𝑐1𝑅
2Ω𝑜𝑘 2
→ 𝑐1 = −
1 − 𝑘2
𝑐1𝑅 𝑐2
0= +
2 𝑅
𝑐1𝑅 2
𝑐2 = −
2
2Ω𝑜𝑘 2 𝑅 2
𝑐2 =
1 − 𝑘2 2
Ω𝑜𝑘 2 𝑅 2
𝑐2 =
1 − 𝑘2
𝑐1𝑟 𝑐2
𝑣𝜃 = +
2 𝑟
2Ω𝑜𝑘 2 Ω𝑜𝑘 2 𝑅 2
− 2𝑟 2
𝑣𝜃 = 1 − 𝑘 + 1 − 𝑘
2 𝑟
Ω𝑜𝑘 2 Ω𝑜𝑘 2 𝑅 2
𝑣𝜃 = − 𝑟+
1 − 𝑘2 𝑟(1 − 𝑘 2 )
2
Ω𝑜𝑘 𝑘 2𝑅2
𝑣𝜃 = [−𝑟 + ]
1 − 𝑘2 𝑟
𝑑 𝑣𝜃
𝜏𝑟𝜃 = 𝜇𝑟 ( )
𝑑𝑟 𝑟
Ω𝑜𝑘 2 𝑘 2𝑅2
[−𝑟 +
𝑑 1−𝑘 2 𝑟 ]
𝜏𝑟𝜃 = 𝜇𝑟 [ ]
𝑑𝑟 𝑟

[37]
 𝑘 2𝑅2
μrΩ𝑜𝑘 2 𝑑 −𝑟 + 𝑟
𝜏𝑟𝜃 = [ ]
1 − 𝑘 2 𝑑𝑟 𝑟

μrΩ𝑜𝑘 2 2𝑘 2 𝑅 2
𝜏𝑟𝜃 = (− )
1 − 𝑘2 𝑘3
2μrΩ𝑜𝑘 4 𝑅 2
 1. 𝜏𝑟𝜃 = 2
𝑟 (1 − 𝑘 2 )
𝜉 = 𝜏𝑟𝜃 × 𝐴 × 𝐵𝑟𝑎𝑧𝑜 𝑑𝑒 𝑝𝑎𝑙𝑎𝑛𝑐𝑎
 𝜉 = 𝜏𝑟𝜃 [2𝜋𝑘𝑅𝐿] × 𝑘𝑅

2.8 El No. Reynolds es una función de ρ, µ, v, L, característica, establecer el número de Reynolds

mediante análisis dimensional.
∅(𝜌, 𝜇, 𝑣, 𝐿) = 0
𝑁𝑢𝑚𝑒𝑟𝑜 𝑑𝑒 𝑔𝑟𝑢𝑝𝑜𝑠 𝜋 = 4 − 3 = 1
𝜋1 = 𝜌𝑎 𝜇 𝑏 𝑣 𝑐 𝐿 𝜋 = 𝑀0 𝐿0 𝑇 0
𝑀0 𝐿0 𝑇 0 = (𝑀𝐿−3 )𝑎 (𝑀𝐿−1 𝑇 −1 )𝑏 (𝐿𝑇 −1 )𝑐 (𝐿)
𝑀: 0=𝑎+𝑏 → 𝑎=1
𝐿: 0 = −3𝑎 − 𝑏 + 𝑐 + 1 → 𝑏 = −1
𝑇: 0 = −𝑏 − 𝑐 →𝑐=1
𝜋1 = 𝜌1 𝜇 −1 𝑣 1 𝐿
𝜌𝑣𝐿
 𝜋 = 𝑅𝑒 =
𝜇
2.9 Suponiendo que la fuerza de arrastre ejercido sobre un cuerpo sumergido en una corriente
fluida es función de la densidad, la viscosidad y la velocidad del fluido, y una longitud
característica del cuerpo, desarrollar la ecuación general.
∅(𝐹, 𝜌, 𝜇, 𝑣, 𝐿) = 0
𝑁𝑢𝑚𝑒𝑟𝑜 𝑔𝑟𝑢𝑝𝑜𝑠 𝜋 = 5 − 3 = 2
𝜋1 = 𝜌𝑎 𝜇 𝑏 𝐿𝑐 𝐹
𝜋2 = 𝜌𝑎 𝜇 𝑏 𝐿𝑐 𝜇
𝑀0 𝐿0 𝑇 0 = (𝑀𝐿−3 )𝑎 (𝐿𝑇 −1 )𝑏 (𝐿)𝑐 (𝑀𝐿𝑇 −2 )
𝑀: 0 = 𝑎 + 1 → 𝑎 = −1
𝐿: 0 = −3𝑎 + 𝑏 + 𝑐 + 1 → 𝑏 = −2
𝑇: 0 = −𝑏 − 2 → 𝑐 = −2
[38]
 𝐹
 𝜋1 =
𝜌𝑣 2 𝐿2

𝑀0 𝐿0 𝑇 0 = (𝑀𝐿−3 )𝑎1 (𝐿𝑇 −1 )𝑏1 (𝐿)𝑐1 (𝑀𝐿−1 𝑇 −1 )

𝑀: 0 = 𝑎1 + 1 → 𝑎1 = −1
𝐿: 0 = −3𝑎1 + 𝑏1 + 𝑐1 − 1 → 𝑐1 = −1
𝑇: 0 = −𝑏1 − 1 → 𝑏1 = −1
𝐹
 𝜋2 = 𝑅𝑒 =
𝜌𝑣𝐿
𝐹
𝑘𝑅𝑒 =
𝜌𝑣 2 𝐿2
2
𝐹 = 𝜌𝑣 2 𝐿2 (𝑘𝑅𝑒)
2
𝜌𝑣 2 2 𝑣 2 𝐿2
𝐹 = 2𝑘𝑅𝑒 ( ) 𝐿 = 2𝑘𝑅𝑒 ( )𝜌
2 2

𝐿2 = 𝑅 𝐶𝑝 = 2𝑘𝑅𝑒
𝑅𝑣 2
𝐹 = 2𝑘𝑅𝑒 𝜌
2
𝑅𝑣2
 𝐹 = 𝐶𝑝 𝜌
2

2.10 Desarrollar una expresión que de las pérdidas de carga en una tubería horizontal, para un
flujo turbulento incompresible la perdida de carga viene dada por lucida de presión y es una
medida de la resistencia presentada al flujo a través de la tubería. La resistencia es una función
del diámetro de la tubería, viscosidad y la densidad del fluido, la longitud de la tubería, la
velocidad y la rugosidad (K) del fluido.

∅(∆𝑃, 𝐷, 𝜇, 𝜌, 𝐿, 𝑣, 𝑘) = 0
𝑁𝑜. 𝑔𝑟𝑢𝑝𝑜𝑠 𝜋 = 7 − 3 = 4
𝜋1 = 𝐷𝑎 𝜌𝑏 𝑣 𝑐 ∆𝑃
𝜋2 = 𝐷𝑎1 𝜌𝑏1 𝑣 𝑐1 𝜇
𝜋3 = 𝐷𝑎2 𝜌𝑏2 𝑣 𝑐2 𝐿
𝜋4 = 𝐷𝑎3 𝜌𝑏3 𝑣 𝑐3 𝑘

[39]
 𝜋1: 𝑀0 𝐿0 𝑇 0 = (𝐿)𝑎 (𝑀𝐿−3 )𝑏 (𝐿𝑇 −1 )𝑐 (𝑀𝐿−1 𝑇 −2 )
𝑀: 0 = 𝑏 + 1 → 𝑏 = −1
𝐿: 0 = 𝑎 − 3𝑏 + 𝑐 − 1 →𝑎=0
𝑇: 0 = −𝑐 − 2 → 𝑐 = −2
∆𝑃
 𝜋1 :
𝜌𝑣 2
𝜋2: 𝑀0 𝐿0 𝑇 0 = (𝐿)𝑎1 (𝑀𝐿−3 )𝑏1 (𝐿𝑇 −1 )𝑐1 (𝑀𝐿−1 𝑇 −1 )
𝑀: 0 = 𝑏1 + 1 → 𝑏1 = −1
𝐿: 0 = 𝑎1 − 3𝑏1 + 𝑐1 − 1 → 𝑎1 = −1
𝑇: 0 = −𝑐1 − 1 → 𝑐1 = −1
1
 𝜋2 = = 𝑅𝑒
𝐷𝜌𝑣
𝜋3: 𝑀0 𝐿0 𝑇 0 = (𝐿)𝑎2 (𝑀𝐿−3 )𝑏2 (𝐿𝑇 −1 )𝑐2 (𝐿)
𝑀: 0 = 𝑏2 → 𝑏2 = 0
𝐿: 0 = 𝑎2 − 3𝑏2 + 𝑐2 + 1 → 𝑎2 = −1
𝑇: 0 = −𝑐2 + 0 → 𝑐2 = 0
𝐿
 𝜋3 =
𝐷
𝜋4: 𝑀0 𝐿0 𝑇 0 = (𝐿)𝑎3 (𝑀𝐿−3 )𝑏3 (𝐿𝑇 −1 )𝑐3 (𝐿)
𝑀: 0 = 𝑏3 → 𝑏3 = 0
𝐿: 0 = 𝑎3 − 3𝑏3 + 𝑐3 + 1 → 𝑎3 = −1
𝑇: 0 = −𝑐3 + 0 → 𝑐3 = 0
𝑘
 𝜋4 =
𝐷
∆𝑃 𝐿 𝑘
 2
= 𝛽(𝑅𝑒) ( ) ( )
𝜌𝑣 𝐷 𝐷

[40]
2.11 La velocidad del flujo U, a través de una placa de orificio tiene que ser determinada como
una función del diámetro del tubo que ser determinada como una función del diámetro del
tubo, densidad, diámetro del orificio, viscosidad y caída de presión. Determine el número de
grupos adimensionales.
∅(𝑈, 𝜌, 𝐷𝑜, 𝜇, ∆𝑃, 𝐷𝑡) = 0

𝑁𝑜. 𝑔𝑟𝑢𝑝𝑜𝑠 𝜋 = 6 − 3 = 3
𝜋1 = 𝜌𝑎 ∆𝑃𝑏 𝐷𝑡 𝑐 𝑈
𝜋2 = 𝜌𝑎1 ∆𝑃𝑏1 𝐷𝑡 𝑐1 𝐷𝑜
𝜋3 = 𝜌𝑎2 ∆𝑃𝑏2 𝐷𝑡 𝑐2 𝜇
𝜋1 : 𝑀0 𝐿0 𝑇 0 = (𝑀𝐿−3 )𝑎 (𝑀𝐿−1 𝑇 −2 )𝑏 (𝐿)𝑐 (𝐿3 𝑇 −1 )
𝑀: 0 = 𝑎 + 𝑏 → 𝑎 = 1/2
𝐿: 0 = −3𝑎 − 𝑏 − 𝑐 + 3 → 𝑐 = −2
1
𝑇: 0 = −2𝑏 − 1 → 𝑏 = −
2
1
𝜌2 𝑈
𝜋1 = 1
∆𝑃2 𝐷𝑡 2
𝜋2 : 𝑀0 𝐿0 𝑇 0 = (𝑀𝐿−3 )𝑎1 (𝑀𝐿−1 𝑇 −2 )𝑏1(𝐿)𝑐1 (𝐿)
𝑀: 0 = 𝑎1 + 𝑏1 → 𝑎1 = 0
𝐿: 0 = −3𝑎1 − 𝑏1 + 𝑐1 + 1 → 𝑐1 = −1
𝑇: 0 = −2𝑏1 → 𝑏1 = 0
𝐷𝑜
𝜋1 =
𝐷𝑡
𝜋3 : 𝑀0 𝐿0 𝑇 0 = (𝑀𝐿−3 )𝑎2 (𝑀𝐿−1 𝑇 −2 )𝑏2 (𝐿)𝑐2 (𝑀𝐿−1 𝑇 −1 )
𝑀: 0 = 𝑎2 + 𝑏2 + 1 → 𝑎2 = −1/2
𝐿: 0 = −3𝑎2 − 𝑏2 + 𝑐2 − 1 → 𝑐2 = −1
𝑇: 0 = −2𝑏2 − 1 → 𝑏2 = −1/2
𝜇
𝜋3 =
𝜌1/2 ∆𝑃1/2 𝐷𝑡

[41]
2.12 El caudal a través de una placa de orificio es determinado como una función del Dt, ρ, Do, µ,
ΔP. Determinar una ecuación con análisis dimensional.

𝑄 = 𝑓(𝐷, 𝜌, 𝐷𝑜, 𝜇, ∆𝑃) = 0

𝐿3 𝑇 −1 = 𝑘(𝐷)𝑎 (𝜌)𝑏 (𝐷𝑜)𝑐 (𝜇)𝑑 (ΔP)𝑒
𝐿3 𝑇 −1 = 𝑘(𝐿)𝑎 (𝑀𝐿−3 )𝑏 (𝐿)𝑐 (𝑀𝐿−1 𝑇 −1 )𝑑 (M𝐿−1 𝑇 −2 )𝑒

𝑀: 0 = 𝑏 + 𝑑 + 𝑒
𝐿: 3 = 𝑎 − 3𝑏 + 𝑐 − 𝑑 − 𝑒
𝑇: − 1 = −𝑑 − 2𝑒
Variables importantes: Do (c) µ (d)

1−𝑑 1−𝑑 −𝑑 − 1 1−𝑑

𝑒= 0=𝑏+𝑑+ 3=𝑎−3 +𝑐−𝑑−
2 2 2 2
−2𝑑 − 1 + 𝑑 3 3 1 𝑑
𝑏= 3=𝑎− 𝑑+ +𝑐−𝑑− −
2 2 2 2 2
−𝑑 − 1 2=𝑎+𝑑+𝑐
𝑏=
2
𝑎 =2−𝑑−𝑐

−𝑑−1 1−𝑑
𝑄 = 𝑘(𝐷)2−𝑑−𝑐 (𝜌) 2 (𝐷𝑜)𝑐 (𝜇)𝑑 (ΔP) 2
𝑑
𝑐
𝜇 𝐷𝑜 1 1
𝑄 = 𝑘𝐷2 ( 1 1 ) ( ) (∆𝑃)2 (𝜌)2
𝐷
𝐷𝜌2 ∆𝑃2
𝑄𝜌1/2 𝜇 𝑑 𝐷𝑜 𝑐
= 𝑘 ( 1/2 1/2 ) ( )
𝐷∆𝑃1/2 𝐷𝜌 ∆𝑃 𝐷

2.13 Determinar la Presión dinámica ejercida sobre un cuerpo totalmente sumergido en la

corriente de un fluido incompresible al suponer que la presión es función de la densidad y la
velocidad.
𝑓(𝑃, 𝜌, 𝑣) = 0
𝑁𝑜. 𝐺𝑟𝑢𝑝𝑜𝑠 𝑑𝑒 𝜋 = 3 − 3 = 0
𝑃 = 𝑘(𝜌)𝑎 (𝑣)𝑏
𝑀𝐿−1 𝑇 −2 = 𝑘(𝑀𝐿−3 )𝑎 (𝐿𝑇 −1 )𝑏
𝑀: 1 = 𝑎 →𝑎=1
𝐿: − 1 = −3𝑎 + 𝑏
𝑇: − 2 = −𝑏 →𝑏=2
𝑃 = 𝑘(𝜌)1 (𝑣)2
 𝑃 = 𝑘𝜌𝑣 2

[42]
2.14 Para el caso de un líquido ideal, expresar el caudal (Q) a través de un orificio en función de
la densidad del líquido, el diámetro del orificio y la densidad de las presiones.
Método pi 𝑓(𝑄, 𝜌, 𝐷𝑜, ∆𝑃) = 0
𝑁𝑜. 𝐺𝑟𝑢𝑝𝑜𝑠 𝑑𝑒 𝜋 = 4 − 3 = 1
𝜋1: 𝑀0 𝐿0 𝑇 0 = (𝑀𝐿−3 )𝑎 (𝐿)𝑏 (𝑀𝐿−1 𝑇 −2 )𝑐 (𝐿3 𝑇 −1 )
𝑀: 0 = 𝑎 + 𝑐 → 𝑎 = 1/2
𝐿: 0 = −3𝑎 + 𝑏 − 𝑐 + 3 → 𝑏 = −2
𝑇: 0 = −2𝑐 − 1 → 𝑐 = −1/2
𝜌1/2 𝑄
 𝜋4 =
𝐷𝑜2 ∆𝑃1/2

∆𝑃
 𝑄 = 𝑘𝐷𝑜2 √
𝜌

Método Rayligh 𝑓(𝜌, 𝐷𝑜, ∆𝑃) = 𝜌𝑎 𝐷𝑜𝑏 ∆𝑃𝑐

𝐿3 𝑇 −1 = (𝑀𝐿3 )𝑎 (𝐿)𝑏 (𝑀𝐿−1 𝑇 −2 )𝑐
𝑀: 0 = 𝑎 + 𝑐 → 𝑎 = −1/2
𝐿: 3 = −3𝑎 + 𝑏 − 𝑐 → 𝑏 = 2
𝑇: − 1 = −2𝑐 → 𝑐 = 1/2
 𝑄 = 𝑘𝜌−/2 𝐷𝑜2 ∆𝑃1/2

∆𝑃
 𝑄 = 𝑘𝐷𝑜2 √
𝜌

2.15 Suponiendo que la potencia comunicada a una bomba es función del peso específico del
fluido del Q y la altura comunicada a la corriente, establecer una ecuación por análisis
dimensional.
Método pi
𝑓(𝑃, ℘, 𝑄, 𝐻) = 0
𝑁𝑢𝑚𝑒𝑟𝑜 𝑑𝑒 𝑔𝑟𝑢𝑝𝑜𝑠 𝜋 = 4 − 3 = 1

𝜋1 : ℘𝑎 𝑄 𝑏 𝐻 𝑐 𝑃
𝑀0 𝐿0 𝑇 0 = (𝑀𝐿−2 𝑇 −2 )𝑎 (𝐿3 𝑇 −1 )𝑏 (𝐿)𝑐 (𝑀𝐿2 𝑇 −3 )
𝑀: 0 = 𝑎 + 1 → 𝑎 = −1
𝐿: 0 = −2𝑎 + 3𝑏 + 𝑐 + 2 → 𝑐 = −1
𝑇: 0 = −2𝑎 − 𝑏 − 3 → 𝑏 = −1
𝑃
𝜋1 :
℘𝑄𝐻
 𝑃 = 𝑘℘𝑄𝐻

[43]
 Metodo Rayleight
𝑃 = 𝑓(℘𝑄𝐻)

(𝑀𝐿2 𝑇 −3 ): ℘𝑎 𝑄 𝑏 𝐻 𝑐
(𝑀𝐿2 𝑇 −3 ) = (𝑀𝐿−2 𝑇 −2 )𝑎 (𝐿3 𝑇 −1 )𝑏 (𝐿)𝑐

𝑀: 1 = 𝑎 → 𝑎 = 1
𝐿: 2 = −2𝑎 + 3𝑏 + 𝑐 → 𝑐 = 1
𝑇: − 3: − 2𝑎 − 𝑏 → 𝑏 = 1

 P = k℘𝑄𝐻

2.16 Si se realiza el experimento de Reynolds con un tubo Pl de 38 mm ¿Qué velocidad de flujo

ocurrirá en la transición?
𝐷𝜌𝑣
𝐷 = 38𝑚𝑚 𝑅𝑒 = = 2300
𝜇
𝑘𝑔
𝜌 = 996.1
𝑇 = 300°𝐾 { 𝑚3 }
−6 −6
𝑘𝑔
𝜇 = 875.75𝑥10 𝑃𝑎 − 𝑠 = 875.75𝑥10
𝑚
𝜌𝐷𝑣
𝑅𝑒 =
𝜇
𝑘𝑔
(996.1 ) (0.038𝑚)𝑣
𝑚3
2300 =
𝑘𝑔
875.75𝑥10−6 𝑚
𝑚
𝑣 = 0.053
𝑠
𝑇 = 27°
𝜈 = 0.8788𝑥10−6
𝐷𝑣
𝑅𝑒 =
𝜈
(0.038𝑚)𝑣
2300 =
0.878810−6
0.053 𝑚/𝑠

[44]
2.17 Estudie el flujo de aire a 30 m/seg. A lo largo de una placa plana. ¿A qué distancia del borde
de ataque ocurrirá la transición?
𝑚
𝑣 = 30
𝑠
𝜌𝑣𝐿
𝑅𝑒𝑥 = 𝜌𝑣𝑥 𝑅𝑒 =
𝜇
𝑎𝑖𝑟𝑒 𝑇 = 27°𝐶 = 300°𝐾
𝜈 = 1.5689𝑥10−5
𝑣𝑥
𝑅𝑒𝑥 =
𝜈
30𝑚
5
2𝑥10 = 𝑠 𝑥
1.5689𝑥10−5
 𝑥 = 0.1045𝑚

2.18 Los modernos aviones se han perfeccionado tanto que el 75% del arrastre parasito (porción
del arrastre total del avión que no esté directamente asociado con la producción de empuje de
sustentación) puede atribuirse a la friccion a lo largo de las superficies externas. Para un avión
subsónico de propulsión, el coeficiente de arrastre parasitario basado en el área del ala es de
0.011. Determine el arrastre de la fricción en un avión como este.

𝐶𝐷 = 0.75(0.011)
𝑚𝑖𝑙𝑙𝑎𝑠 𝑓𝑡
𝑎. 𝑣∞ = 500 = 733.33
ℎ 𝑠
ℎ = 35000 𝑓𝑡
𝐴 = 2400 𝑓𝑡 2
𝑠𝑙𝑢𝑔 𝑙𝑚
𝜌 = 0.000737 = 0.0237 𝑓𝑡 3 1𝑠𝑙𝑢𝑔 = 32.174 𝑙𝑏𝑚
𝑓𝑡 3

𝜌𝑣∞ 2
𝐹 = 𝐴𝐶𝐷
2
0.0237𝑙𝑏𝑚 733.33𝑓𝑡 2
( )( )
𝑓𝑡 3 𝑠
𝐹 = (2400𝑓𝑡 2 )(0.75 ∗ 0.011)
2
𝑓𝑡
𝐹 = 126177.6529 𝑙𝑏𝑚 − 2
𝑠
𝑙𝑏𝑚 ∗ 𝑓𝑡
𝑝𝑎𝑟𝑎 𝑐𝑜𝑛𝑣𝑒𝑟𝑡𝑖𝑟 𝑎 𝑢𝑛𝑖𝑑𝑎𝑑𝑒𝑠 𝑓𝑢𝑒𝑟𝑧𝑎 𝑠𝑒 𝑑𝑖𝑣𝑖𝑑𝑒 𝑒𝑛𝑡𝑟𝑒 𝑔𝑐 = 32.174
𝑙𝑏𝑓 − 𝑠 2

 𝐹 = 3921.72 𝑙𝑏𝑓

[45]
 𝑚𝑖𝑙𝑙𝑎𝑠 𝑓𝑡
𝑏. 𝑣∞ = 200 = 293.33
ℎ 𝑠
ℎ = 0 𝑓𝑡
𝐴 = (0.75 ∗ 0.011) 𝑓𝑡 2
𝑠𝑙𝑢𝑔 𝑙𝑚
𝜌 = 0.002378 = 0.0765 𝑓𝑡 3
𝑓𝑡 3

𝜌𝑣∞ 2
𝐹 = 𝐴𝐶𝐷
2
0.07565𝑙𝑏𝑚 293.33𝑓𝑡 2
( )( )
2 𝑓𝑡 3 𝑠
𝐹 = (2400𝑓𝑡 )(0.75 ∗ 0.011)
2
𝑓𝑡
65164.27 𝑙𝑏𝑚 − 2
𝐹= 𝑠
𝑙𝑏𝑚 ∗ 𝑓𝑡
32.174
𝑙𝑏𝑓 − 𝑠 2

 𝐹 = 2025.37 𝑙𝑏𝑓

2.19 Calcule la fuerza de arrastre de una antena cuya longitud es de 3ft y diámetro promedio de
0.2 plg a una velocidad de 60 mph.
𝐿 = 3 𝑓𝑡
𝐷 = 0.2 𝑝𝑙𝑔
𝑚𝑖𝑙𝑙𝑎𝑠 𝑓𝑡
𝑣∞ = 60 = 88
ℎ 𝑠
𝜌𝑣∞ 2 𝜌𝐷𝑣 𝐷𝑣
𝐹 = 𝐴𝐶𝐷 𝑅𝑒 = 𝑅𝑒 =
2 𝜇 𝜈

𝜌 = 0.0764 𝑙𝑏𝑚/𝑓𝑡 3
𝑇 = 60 °𝐹 { }
𝜐 = 0.159𝑥10−3 𝑓𝑡 2 /𝑠
88𝑓𝑡
𝐷𝑣 ( 𝑠 ) (0.2/12𝑓𝑡)
𝑅𝑒 = = = 9224.31
𝜈 0.159𝑥10−3 𝑓𝑡 2 /𝑠

0.07564𝑙𝑏𝑚 88𝑓𝑡 2
( )( )
2 𝑓𝑡 3 𝑠
𝐹 = (0.2/12𝑓𝑡 )(3 ∗ 1.1)
2
𝑓𝑡
16.27 𝑙𝑏𝑚 − 2
𝐹= 𝑠
𝑙𝑏𝑚 ∗ 𝑓𝑡
32.174
𝑙𝑏𝑓 − 𝑠 2

 𝐹 = 0.5057 𝑙𝑏𝑓

[46]
2.20 La frecuencia del desprendimiento de los vórtices para un cilindro puede predecirse por
medio de la ecuación:

𝑣𝑥 19.7
𝑓 = 0.198 ( ) (1 − )
𝐷 𝑅𝑒

Donde f es la fuerza de los vórtices desprendidos desde un lado del cilindro y Re es el número
de Reynolds de corriente libre. Determínese la frecuencia de desprendimiento de vórtices de
un alambre de 0.25 plg de diámetro en un viento de 20 mph. ¿Para qué velocidad del viento
será cero la frecuencia?
𝐷 = 0.25 𝑝𝑙𝑔 = 0.0208𝑓𝑡
𝑚𝑖𝑙𝑙𝑎𝑠 𝑓𝑡
𝑉𝑥 = 20 = 29.33
ℎ 𝑠

0.159𝑥10−3 𝑓𝑡2
𝑇 = 60°𝐹 − 𝜐 =
𝑠
𝑓𝑡
𝑉𝑥𝐷 29.33 𝑠 (0.0208𝑓𝑡)
𝑅𝑒 = =
𝜐 0.159𝑥10−3 𝑓𝑡 2
𝑠
𝑅𝑒 = 3836.88
29.33𝑓𝑡
𝑠 19.7
𝑓 = 0.198 ( ) (1 − )
0.0208𝑓𝑡 3836.88

277.7555
𝑓=
𝑠
𝑆𝑢𝑝𝑜𝑛𝑖𝑒𝑛𝑑𝑜 𝑓 = 0
𝑣𝑥 19.7
0 = 0.198 ( ) (1 − )
𝐷 𝑅𝑒
19.7
0= 1−
𝑅𝑒
𝑅𝑒 = 19.7
𝐷𝑣𝑥
𝑅𝑒 =
𝜈
𝑉𝑥(0.0208𝑓𝑡)
19.7 =
0.159𝑥10−3 𝑓𝑡 2
𝑠
 𝑉𝑥 = 0.1505 𝑓𝑡/𝑠

[47]
UNIDAD
TRES

[48]
EJERCICIOS DE CLASE
Ejercicio 3.1

Ta
T1
T2
T3
T4 Tb

Ley de conservación de la energía

𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑠𝑎𝑙𝑖𝑑𝑎 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑝𝑟𝑜𝑑𝑢𝑐𝑐𝑖ó𝑛

{ 𝑑𝑒 𝑒𝑛𝑒𝑟𝑔𝑖𝑎 } − { 𝑑𝑒 𝑒𝑛𝑒𝑟𝑔𝑖𝑎 } + { 𝑑𝑒 𝑒𝑛𝑒𝑟𝑔𝑖𝑎 }=0
𝑐𝑎𝑙𝑜𝑟𝑖𝑓𝑖𝑐𝑎 𝑐𝑎𝑙𝑜𝑟𝑖𝑓𝑖𝑐𝑎 𝑐𝑎𝑙𝑜𝑟𝑖𝑓𝑖𝑐𝑎

𝑄 − 𝑄𝐸 = 0
𝑄
𝑞= ∴𝑄 =𝑞∙𝐴
𝐴
𝑥2 − 𝑥1 = 𝑒𝑠𝑝𝑒𝑠𝑜𝑟
𝑥3 − 𝑥2 = 𝑟
𝑥4 − 𝑥3 = 𝑟
Región 12 -

l
Q
Q

l 𝑄𝐸 −𝑄𝑆 = 0
𝐴 = 𝑊𝐿
1 12
[𝑞12 𝑤𝐿|𝑥 − 𝑞𝑥+∆𝑥 𝑤𝐿 = 0]
𝑊𝐿 𝑥

[49]
 1 12
[𝑞 | − 𝑞𝑥12 |𝑥 = 0]
∆𝑥 𝑥 𝑥+∆𝑥
𝑞𝑥12 |𝑥+∆𝑥 − 𝑞𝑥12 |𝑥
lim =0
∆𝑥→0 ∆𝑥
𝑑𝑞𝑥12
= 0 ∴ 𝑞𝑥12 = 𝑐1 − 𝑐1 = 𝑞0
𝑑𝑥
⟹ 1. 𝑞𝑥12 = 𝑞0 (1)
1 23
[𝑞 | − 𝑞𝑥23 |𝑥 = 0]
∆𝑥 𝑥 𝑥+∆𝑥
𝑞𝑥23 |𝑥+∆𝑥 − 𝑞𝑥23 |𝑥
lim =0
∆𝑥→0 ∆𝑥
𝑑𝑞𝑥23
= 0 ∴ 𝑞𝑥23 = 𝑐1 − 𝑐1 = 𝑞0
𝑑𝑥
⟹ 2. 𝑞𝑥23 = 𝑞0 (1)
1 34
[𝑞 | − 𝑞𝑥34 |𝑥 = 0]
∆𝑥 𝑥 𝑥+∆𝑥
𝑞𝑥34 |𝑥+∆𝑥 − 𝑞𝑥34 |𝑥
lim =0
∆𝑥→0 ∆𝑥
𝑑𝑞𝑥34
= 0 ∴ 𝑞𝑥34 = 𝑐1 − 𝑐1 = 𝑞0
𝑑𝑥
⟹ 3. 𝑞𝑥34 = 𝑞0 (1)

____________________________________
𝑞𝑥12 = 𝑞0
𝑞𝑥23 = 𝑞0
𝑞𝑥34 = 𝑞0
𝑞𝑥12 = 𝑞𝑥23 = 𝑞𝑥34 = 𝑞0
____________________________________
𝑞𝑥12 = 𝑞0
𝑑𝑇
𝑞 = ±𝑘 𝐿𝑒𝑦 𝑑𝑒 𝐹𝑜𝑢𝑟𝑖𝑒𝑟
𝑑𝑥
𝑑𝑇
𝑞0 = −𝑘12
𝑑𝑥

[50]
 𝑇2
𝑞0
∫ 𝑑𝑇 = − ∫ 𝑑𝑥
𝑇1 𝑘12
𝑞0
𝑇2 − 𝑇1 = − (𝑥 − 𝑥1 )
𝑘12 2
𝑞0
⇒ 4. 𝑇1 − 𝑇2 = 12 (𝑥2 − 𝑥1 )
𝑘
𝑞𝑥23 = 𝑞0
𝑑𝑇
𝑞 = ±𝑘 𝐿𝑒𝑦 𝑑𝑒 𝐹𝑜𝑢𝑟𝑖𝑒𝑟
𝑑𝑥
𝑑𝑇
𝑞0 = −𝑘 23
𝑑𝑥
𝑇3
𝑞0
∫ 𝑑𝑇 = − ∫ 𝑑𝑥
𝑇2 𝑘 23
𝑞0
𝑇3 − 𝑇2 = − (𝑥 − 𝑥2 )
𝑘 23 3
𝑞0
⇒ 5. 𝑇2 − 𝑇3 = 23 (𝑥3 − 𝑥2 )
𝑘
𝑞𝑥34 = 𝑞0
𝑑𝑇
𝑞 = ±𝑘 𝐿𝑒𝑦 𝑑𝑒 𝐹𝑜𝑢𝑟𝑖𝑒𝑟
𝑑𝑥
𝑑𝑇
𝑞0 = −𝑘 34
𝑑𝑥
𝑇4
𝑞0
∫ 𝑑𝑇 = − ∫ 𝑑𝑥
𝑇3 𝑘 34
𝑞0
𝑇4 − 𝑇3 = − (𝑥 − 𝑥3 )
𝑘 34 4
𝑞0
⇒ 6. 𝑇3 − 𝑇4 = 34 (𝑥4 − 𝑥3 )
𝑘
𝑞 = ℎ∆𝑇
𝑞𝑥 = 𝑞0
⇒ 7. 𝑞0 = ℎ𝑖 (𝑇2 − 𝑇1 )
⇒ 8. 𝑞0 = ℎ𝑒 (𝑇𝑎 − 𝑇𝑏 )
𝑞0
𝑇𝑎 − 𝑇1 =
ℎ𝑖

[51]
 𝑞0
𝑇4 − 𝑇𝑏 =
ℎ𝑒

𝑞0
𝑇2 − 𝑇1 = (𝑥 − 𝑥1 )
𝑘12 2
𝑞0
𝑇3 − 𝑇2 = 23 (𝑥3 − 𝑥2 )
𝑘
𝑞0
𝑇4 − 𝑇3 = 34 (𝑥4 − 𝑥3 )
𝑘
𝑞0
𝑇2 − 𝑇1 =
ℎ𝑖
𝑞0
𝑇𝑎 − 𝑇𝑏 =
ℎ𝑒
___________________________________
𝑞0 𝑞 𝑞 𝑞 𝑞
𝑇𝑎 − 𝑇𝑏 = 0
+ 𝑘 12 0
(𝑥2 − 𝑥1 ) + 23 0
(𝑥3 − 𝑥2 ) + 34 (𝑥4 − 𝑥3 ) + ℎ0 6
ℎ𝑖 𝑘 𝑘 𝑒

1
𝑞0 = [ ] (𝑇𝑎 − 𝑇𝑏)
1 𝑥2 − 𝑥1 𝑥3 − 𝑥2 𝑥4 − 𝑥3 1
+ + + +
ℎ𝑖 𝑘12 𝑘 23 𝑘 34 ℎ𝑒
𝑞0 = 𝑈(𝑇𝑎 − 𝑇𝑏 )
𝑄 𝑄
𝑞0 = = 𝑈(𝑇𝑎 − 𝑇𝑏 ) ∴ 𝑄 = 𝐴𝑈(𝑇𝑎 − 𝑇𝑏 )
𝐴 𝐴
Ejercicio 3.2

Ta T0
T1
T2
T3

Tb

[52]
 𝐴 = 2𝜋𝑟𝐿
𝑄
𝑄𝐸 − 𝑄𝑆 = 0 𝑞= 𝑄 = 𝑞𝐴
𝐴
𝑅𝑒𝑔𝑖𝑜𝑛 01
[𝑞𝑟01 2𝜋𝐿|𝑡 − 𝑞𝑟01 2𝜋𝑟𝐿|𝑟+∆𝑟 = 0]1
2𝜋𝐿∆𝑟
𝑞𝑟01 𝑟|𝑟+∆𝑟 − 𝑞𝑟01 𝑟|𝑟 /∆𝑟
𝑑𝑞𝑟01 𝑟
=0
𝑑𝑟
𝑞𝑟01 𝑟 = 𝑐1 𝑐1 = 𝑟0 𝑞0

𝑑𝑞𝑟12 𝑟
=0
𝑑𝑟
𝑞𝑟12 𝑟 = 𝑐2 𝑐2 = 𝑟0 𝑞0

𝑑𝑞𝑟23 𝑟
=0
𝑑𝑟
𝑞𝑟23 𝑟 = 𝑐3 𝑐3 = 𝑟0 𝑞0
⇒ 1. 𝑞𝑟01 𝑟 = 𝑟0 𝑞0
⇒ 2. 𝑞𝑟12 𝑟 = 𝑟0 𝑞0
⇒ 3. 𝑞𝑟23 𝑟 = 𝑟0 𝑞0
𝑞𝑟01 𝑟 = 𝑞𝑟12 𝑟 = 𝑞𝑟23 𝑟 = 𝑟0 𝑞0
𝑑𝑇
𝑞𝑟01 𝑟 = 𝑟0 𝑞0 𝑞 = ±𝑘
𝑑𝑟
𝑟0 𝑞0
𝑞𝑟01 =
𝑟
𝑑𝑇 𝑟0 𝑞0
−𝑘 01 =
𝑑𝑟 𝑟
𝑇1
𝑟0 𝑞0 𝑟1 𝑑𝑟
∫ 𝑑𝑇 = − 01 ∫
𝑇0 𝑘 𝑟0 𝑟
𝑟0 𝑞0 𝑟1
𝑇1 − 𝑇0 = − 𝑙𝑛
𝑘 01 𝑟0
𝑟0 𝑞0 𝑟1
⇒ 4. 𝑇0 − 𝑇1 = 𝑙𝑛
𝑘 01 𝑟0

[53]
 𝑟0 𝑞0
𝑞𝑟12 =
𝑟
𝑑𝑇 𝑟0 𝑞0
−𝑘12 =
𝑑𝑟 𝑟
𝑇2
𝑟0 𝑞0 𝑟2 𝑑𝑟
∫ 𝑑𝑇 = − ∫
𝑇1 𝑘12 𝑟1 𝑟
𝑟0 𝑞0 𝑟2
𝑇2 − 𝑇1 = − 𝑙𝑛
𝑘12 𝑟1
𝑟0 𝑞0 𝑟2
⇒ 5. 𝑇1 − 𝑇2 = 𝑙𝑛
𝑘12 𝑟1
𝑟0 𝑞0
𝑞𝑟23 =
𝑟
𝑑𝑇 𝑟0 𝑞0
−𝑘 23 =
𝑑𝑟 𝑟
𝑇3
𝑟0 𝑞0 𝑟3 𝑑𝑟
∫ 𝑑𝑇 = − 23 ∫
𝑇2 𝑘 𝑟2 𝑟
𝑟0 𝑞0 𝑟3
𝑇3 − 𝑇2 = − 𝑙𝑛
𝑘 23 𝑟2
𝑟0 𝑞0 𝑟3
⇒ 6. 𝑇2 − 𝑇3 = 𝑙𝑛
𝑘 23 𝑟2

𝑞𝑟 𝑟 = 𝑟0 𝑞0
𝑞 = ℎ∆𝑇
𝑞𝑟 𝑟0 = 𝑟0 𝑞0
𝑟0 𝑞0
𝑞𝑟 =
𝑟0
𝑟0 𝑞0
= ℎ𝑖 (𝑇𝑎 − 𝑇𝑖 )
𝑟0
𝑟0 𝑞0
⇒ 7. 𝑇𝑎 − 𝑇0 =
𝑟0 ℎ𝑖

𝑞𝑟 𝑟 = 𝑟0 𝑞0
𝑞𝑟 𝑟3 = 𝑟0 𝑞0
𝑟0 𝑞0
𝑞𝑟 =
𝑟3
𝑟0 𝑞0
= ℎ𝑒 (𝑇3 − 𝑇𝑏 )
𝑟3
[54]
 𝑟0 𝑞0
⇒ 8. 𝑇3 − 𝑇𝑏 =
𝑟3 ℎ𝑒
𝑟0 𝑞0 𝑟1
𝑇0 − 𝑇1 = 𝑙𝑛
𝑘 01 𝑟0
𝑟0 𝑞0 𝑟2
𝑇1 − 𝑇2 = 𝑙𝑛
𝑘12 𝑟1
𝑟0 𝑞0 𝑟3
𝑇2 − 𝑇3 = 𝑙𝑛
𝑘 23 𝑟2
𝑟0 𝑞0
𝑇𝑎 − 𝑇0 =
𝑟0 ℎ𝑖
𝑟0 𝑞0
𝑇3 − 𝑇𝑏 =
𝑟3 ℎ𝑒
_____________________________
𝑟0 𝑞0 𝑟1 𝑟0 𝑞0 𝑟2 𝑟0 𝑞0 𝑟3 𝑟0 𝑞0 𝑟0 𝑞0
𝑇𝑎 − 𝑇𝑏 = 𝑙𝑛 + 𝑙𝑛 + 𝑙𝑛 + +
𝑘 01 𝑟0 𝑘12 𝑟1 𝑘 23 𝑟2 𝑟0 ℎ𝑖 𝑟3 ℎ𝑒
𝑟 𝑟 𝑟 −1
1 𝑙𝑛 𝑟1 𝑙𝑛 𝑟2 𝑙𝑛 𝑟3 1
𝑟0 𝑞0 [ + 010 + 121 + 232 + ] (𝑇𝑎 − 𝑇𝑏 )
𝑟0 ℎ𝑖 𝑘 𝑘 𝑘 𝑟3 ℎ𝑒

2𝜋𝐿(𝑇𝑎 − 𝑇𝑏 )
𝑞0 2𝜋𝑟0 𝐿 = 𝑟 𝑟 𝑟
𝑙𝑛 𝑟1 𝑙𝑛 𝑟2 𝑙𝑛 𝑟3
1 0 1
+ + 12 + 232 +
1
𝑟0 ℎ𝑖 𝑘 01 𝑘 𝑘 𝑟3 ℎ𝑒

𝑟0 2𝜋𝐿(𝑇𝑎 − 𝑇𝑏 )
𝑞 2𝜋𝑟0 𝐿 = 𝑟1 𝑟 𝑟
𝑟0 0 𝑙𝑛 𝑟 𝑙𝑛 𝑟2 𝑙𝑛 𝑟3
1 0 1
+ + 121 + 232 +
[ 𝑟0 ℎ𝑖 𝑘 01 𝑘 𝑘 𝑟3 ℎ𝑒 ]

2𝜋𝑟0𝐿(𝑇𝑎 − 𝑇𝑏 )
𝑞0 𝐴 = 𝑟 𝑟 𝑟
𝑙𝑛 𝑟1 𝑙𝑛 𝑟2 𝑙𝑛 𝑟3
1 1
𝑟0 [ + 010 + 121 + 232 + ]
𝑟0 ℎ𝑖 𝑘 𝑘 𝑘 𝑟3 ℎ𝑒

𝑄
𝑞= 𝑞𝐴 = 𝑄
𝐴
1
𝑄= 𝑟1 𝑟2 𝑟3 ∆(𝑇𝑎 − 𝑇𝑏 )
𝑙𝑛 𝑙𝑛 𝑙𝑛
1 𝑟 𝑟 𝑟 1
𝑟0 [ + 010 + 121 + 232 + ]
𝑟0 ℎ𝑖 𝑘 𝑘 𝑘 𝑟3 ℎ𝑒

𝑄 = 𝑈𝐴(𝑇𝑎 − 𝑇𝑏 )
[55]
Conducción de estado permanente
Ejercicio 3.3

Fluye vapor saturado a 40 psia, a 5 ft/s a través de una tubería de 1 ½ in de acero de

cedula 40. El coeficiente de transferencia de calor convectivo para condensar el
vapor en la superficie interior es de 1500 Btu/ h ft2 °F. El aire circundante se encuentra
a 80°F y el coeficiente de la superficie exterior es de 3 Btu/h ft2 °F.
Determine lo siguiente:

a. La pérdida de calor por caída de 10 ft de tubo sin aislamiento

b. La pérdida de calor por cada 10 ft de tubo aislado con 2 in de 85% de aislante
de magnesia.
c. La masa del vapor condensado en 10 ft de tubo sin aislamiento.
Tb
Acero

Cedula 40I1/2

hi=1500
V=5 ft/s he=3
Ta ρ=40 lb/pgls

a) Encontrar el diámetro externo e interno.

1 900 0.95
𝐷𝑒 = 1 900 𝑖𝑛 ∴ 𝑟1 = = 0.95 𝑖𝑛 𝑎 𝑓𝑡 ∴ = 0.07916 𝑓𝑡
2 12
1 610 0.805
𝐷𝑖 = 1 610 𝑖𝑛 ∴ 𝑟0 = = 0.805 𝑖𝑛 𝑎 𝑓𝑡 ∴ = 0.06708 𝑓𝑡
2 12
68
𝑇𝑎 = 267.26°𝐹 {212}
512
𝐵𝑡𝑢
𝐾𝑎𝑐𝑒𝑟𝑜 = 𝑘01 = 24.8
ℎ − 𝑓𝑡 − °𝐹
𝑇𝑎 − 𝑇𝑏
𝑞0 𝑟1 =
𝑟1 𝑟2 𝑟3
1 𝑙𝑛 𝑟𝑜 𝑙𝑛 𝑤 𝑟2 1
𝑟1
+ 01 + 12 + 23 +
𝑟0 ℎ1 𝑘 𝑘 𝑘 𝑟3 ℎ2
2𝜋𝐿(𝑇𝑎 − 𝑇𝑏 )
𝑄=
𝑟1
1 𝑙𝑛 𝑟0 1
+ 01 +
𝑟0 ℎ1 𝑘 𝑟3 ℎ3

[56]
 2𝜋(10)(267.26 − 80)
𝑄=
0.07916
1 𝑙𝑛 0.06708 1
+ +
0.06708(1500) 24.8 0.07916(3)
𝐵𝑡𝑢
𝑄 = 2783.1824
ℎ
b) B
2𝜋𝐿(𝑇𝑎 − 𝑇𝑏 )
𝑄=
𝑟1 𝑟2
1 𝑙𝑛 𝑟0 𝑙𝑛
+ 01 + 12 𝑟1 + 1
𝑟0 ℎ1 𝑘 𝑘 𝑟2 ℎ𝑒
𝑟1 = 0.07916 𝐹𝑡
𝑟0 = 0.06708 𝐹𝑡
2
𝑟2 = 0.07916 + = 0.2458 𝐹𝑡
12
𝑘 01 = 24.8
𝑘12 = 0.041

2𝜋(10)(267.26 − 80)
𝑄=
0.07916 0.2458
1 𝑙𝑛 0.06708 𝑙𝑛 0.07916 1
+ + 0.041 +
0.06708(1500) 24.8 0.2458(3)
𝐵𝑡𝑢
𝑄 = 405.6081
ℎ
c) C
𝑄̂ − 𝑊
̂ = ∆𝐻
̂ + ∆𝐸𝑐 + ∆𝐸𝑝

𝑄̂ = ∆𝐻
̂

̂
𝑄 = 𝑚̇∆𝐻
Vapor
𝑄 = 𝑚̇(𝐻𝑟 − 𝐻𝑒) = 𝑚̇𝜆
2783.1824
𝑚̇ =
993.6
𝑘𝑔
𝑚̇ = 2.9811
ℎ

[57]
Conducción de estado permanente.
Ejercicio 3.4

La pared de un horno, que consta de 0.25 m de ladrillo refractario, 0.20 m de caolín y

una capa exterior de tabique del usado en la construcción de 0.10 m, expuesta al gas
del horno, a una temperatura de 1 370 K, con aire de 300 K adyacente a la pared
exterior. Los coeficientes de transferencia de calor convectivo interior y exterior son
115 y 23 W/m2 K, respectivamente. Determine la perdida de calor por pie cuadrado
de pared y la temperatura de la superficie exterior de la pared en estas condiciones.

0.25 0.20 0.15

T0 R C L hi=115 w/m2-°K

Hi he=23 w/m2-°K

370 °K
He

300°K

𝑇𝑎 − 𝑇𝑏
𝑞=
1 𝑥1 − 𝑥0 𝑥2 − 𝑥1 𝑥3 − 𝑥2 1
+ + + +
ℎ𝑖 𝑘 01 𝑘12 𝑘 23 ℎ𝑒
𝑊
𝑘 01 = 𝑘𝑅 = 1.13
𝑚2 °𝐾
𝑊
𝑘12 = 𝑘𝐶 = 0.09
𝑚2 °𝐾
𝑊
𝑘 23 = 𝑘𝐿 = 0.66
𝑚2 °𝐾
1370 − 300
𝑞=
1 0.25 0.23 0.1 1
+ + + +
115 1.13 0.09 0.66 23
𝑤
𝑎. 𝑞 = 404.208 2
𝑚

0.25 0.20 0.15

L
Ta R C
T0
Hi T1
T2
370 °K Tb
T3
He

300°K
[58]
 𝑞 = ℎ∆𝑇
𝑞 = ℎ𝑒(𝑇3 − 𝑇)
404.208 = 23(𝑇3 − 300)
𝑏. 𝑇3 = 317.5740 °𝐾
𝑇2 =?
𝑞
𝑇2 − 𝑇3 = (𝑥 − 𝑥3 )
𝑘 23 2
404.200
𝑇2 − 317.5742 = (0.1)
0.6
𝑐. 𝑇2 = 378.8178 °𝐾
Conducción de estado permanente.
Ejercicio 3.5

¿Cuál sería el resultado del problema 17.30 si el fluido que circunda al alambre
aislado se mantuviera a 70 °F y su coeficiente de transferencia de calor entre el
material aislante y el fluido es de 4 Btu/h-Ft2-°F? ¿Cuál sería la temperatura superficial
del aislamiento bajo estas condiciones?
Cobre 3/16 in

T= 70 °F T= 120 °F
𝐵𝑡𝑢
𝑘 = 0.14
ℎ − 𝐹𝑡 − °𝐹
3 1
3′′ 𝑟0 = ( ) = 0.0937′′
𝐷= ∴ { 16 2 }
16
𝑟1 = 0.0937 + 4 = 4.0937′′
2𝜋𝐿(120 − 70)°𝐹
𝑄=
𝑟1
1 𝑙𝑛 𝑟0 1
+ 01 +
𝑟0 ℎ𝑖 𝑘 𝑟1 ℎ2
𝑄 2𝜋𝐿(120 − 70)°𝐹
=
𝐿 4.0937
𝑙𝑛 0.0937
𝐵𝑡𝑢
0.14
ℎ − 𝐹𝑡 − °𝐹
[59]
 𝑄 𝐵𝑡𝑢
= 11.64
𝐿 𝐹𝑡 − ℎ

𝑄 𝐵𝑡𝑢 1𝑤 𝑤
= 11.64 | | = 3.40
𝐿 𝐹𝑡 − ℎ 3.41 𝐵𝐹 𝐹𝑡
ℎ
𝑄 = 𝑖 2𝑅
𝑄
= 𝑖 2𝑅
𝐿
_______________
𝜌𝐿
𝑅=
𝐴
𝜌(𝑅𝑒𝑠𝑖𝑠𝑡𝑖𝑣𝑖𝑑𝑎𝑑) = 1.72𝑥10−6
1.72𝑥10−6 Ω − 𝑐𝑚𝐿
𝑅=
𝜋𝐷2
4
𝑅 1.72𝑥10−6 Ω − 𝑐𝑚 10 − 2𝑚 𝐹𝑡
= 2 | || |
𝐿 3 2
1𝑐𝑚 0.3098
𝜋( ) 𝐹𝑡
16(12)
4
𝑅 Ω
= 2.9429𝑥10−4
𝐿 𝐹𝑡
𝑄 = 𝑖 2𝑅
𝑄
= 𝑖 2𝑅
𝐿
_______________

𝑄
𝑄
𝑖 = √𝐿 = √
𝑅 𝑅
𝐿

3.413𝑤
𝑖=√
2.9429𝑥10−4 Ω

𝑖 = 107.6912 𝑎𝑚𝑝

[60]
Conducción de estado permanente
Ejercicio 3.6

La pérdida de calor que sufre un calentador de agua debe mantenerse en un máximo

de 900 Btu/h-ft2 de área de pared. ¿De qué grosor se requiere que sea el asbesto
(k=0.10 Btu/h-ft°F) si la superficie interior y exterior del aislamiento deben estar a 1600°F
y 500 °F respectivamente.

X0 x1

1600°F

500°F

Espesor= x1-x0
T0 T1

𝐵𝑡𝑢
𝐵𝑡𝑢 𝑄( )
𝑞 = 900 = ℎ
ℎ − 𝐹𝑡 2 𝐴
𝐵𝑡𝑢
𝑘 = 0.1
ℎ − 𝐹𝑡 − °𝐹
𝑞
𝑇0 − 𝑇1 = 01
𝑘 (𝑥1 − 𝑥0 )
𝐵𝑡𝑢
900
(1600 − 500)°𝐹 = ℎ − 𝐹𝑡 2 (𝑥 − 𝑥 )
𝐵𝑡𝑢 1 0
0.1
ℎ − 𝐹𝑡 − °𝐹
(𝑥1 − 𝑥0 ) = 0.1222 𝐹𝑡

[61]
Conducción de estado permanente
Ejercicio 3.7

Un calentador formado por alambre de nicromo enrollado hacia adelante y luego

hacia atrás, con el enrollado muy junto, está cubierto en ambos lados, con una capa
de asbesto (k=0.15 Btu/h-ft-°F, de 1/8 pulgadas de espesor y después con otra capa
de 1/8 de grosor de acero inoxidable (k=10 Btu/h-ft-°F. si la temperatura central de
esta construcción en forma de emparedado se considera constante, y el coeficiente
exterior de transferencia de calor convectivo es 3 Btu /h-ft2-°F, a 1000 °F, ¿cuánta
energía deberá suministrársele al calentador en W/ft2? ¿Cuál será la temperatura
exterior del acero inoxidable.
Tb

T1 Acero
T0 Asbesto
Ta
Ta Alambre

Asbesto
acero
𝑇𝑎 − 𝑇𝑏
𝑞=
1 𝑥0 − 𝑥𝑎 𝑥1 − 𝑥0 1
+ + +
ℎ𝑖 𝑘 0𝑎 𝑘 01 ℎ𝑒
Suponer T = 70 °F
Ta = 1000 °F
1000 − 70
𝑞=
1 1
𝐹𝑡 𝐹𝑡
8(12) 8(12) 1
𝐵𝑡𝑢 + 𝐵𝑡𝑢 + 𝐵𝑡𝑢
0.15 10 3
ℎ − 𝐹𝑡 − °𝐹 ℎ − 𝐹𝑡 − °𝐹 ℎ − 𝐹𝑡 2 − °𝐹
𝐵𝑡𝑢
2303.009
𝑞= ℎ − 𝐹𝑡 2
3.41 𝑤𝑎𝑡𝑡𝑠
𝑤
𝑞 = 675.36 2
𝐹𝑡
𝑤 𝑤
𝑞𝑡𝑜𝑡𝑎𝑙 = 2 (675.36 2 ) 𝑞𝑡𝑜𝑡𝑎𝑙 = 1350.73 2
𝐹𝑡 𝐹𝑡
𝑞 = ℎ∆𝑇 𝑞 = ℎ𝑒 (𝑇1 − 𝑇𝑏 )
𝐵𝑡𝑢 𝐵𝑡𝑢
2303.009 =3 (𝑇 − 70°𝐹)
ℎ − 𝐹𝑡 2 ℎ − 𝐹𝑡 2 − °𝐹 1

[62]
 𝑇1 = 837.66 °𝐹

Aletas de enfriamiento.
Ejercicio 3.8
Ta
Ta
r
q
q
z
Tw
z z+Δz
Ta
h Ta

𝑄𝐸 − 𝑄𝑆 ± 𝑝𝑒𝑟𝑑𝑖𝑑𝑎 𝑜 𝑝𝑟𝑜𝑑𝑢𝑐𝑐𝑖ó𝑛 𝑑𝑒 𝑐𝑎𝑙𝑜𝑟 = 0

𝑄 = ℎ𝐴 (𝑇 − 𝑇𝑎)
𝑞 = ℎ(𝑇 − 𝑇𝑎)
𝑄
𝑞= 𝐴 = 𝐴𝑞
𝐴
1
[𝑞𝑧 ∗ 2𝛽𝑤|𝑧 − 𝑞𝑧 2𝛽𝑤|𝑧+∆𝑧 − 2ℎ(𝑇 − 𝑇𝑎)∆𝑧𝑤 = 0]
2𝛽𝑤∆𝑧
𝑞𝑧 |𝑧 − 𝑞𝑧 |𝑧+∆𝑧 ℎ(𝑇 − 𝑇𝑎)
−
∆𝑧 𝛽
𝑑𝑞𝑧
lim → 0 = −ℎ (𝑇 − 𝑇𝑎)/𝛽
∆𝑧 𝑑𝑧
𝑑𝑇
𝐿𝑒𝑦 𝑑𝑒 𝐹𝑜𝑢𝑟𝑖𝑒𝑟 𝑞 = ±ℎ
𝑑𝑧
𝑑 𝑘𝑑𝑇 ℎ(𝑇 − 𝑇𝑎)
(− )=−
𝑑𝑧 𝑑𝑧 𝛽
𝑑 2 𝑇 ℎ(𝑇 − 𝑇𝑎)
=
𝑑𝑧 2 𝑘𝛽
𝐶. 𝐹. 1 𝑧=0 𝑇 = 𝑇𝑤
𝑑𝑇
𝐶. 𝐹. 2 𝑧=𝐿 =0
𝑑𝑧
𝑑2 𝑇 ℎ(𝑇 − 𝑇𝑎) 𝑇𝑤 − 𝑇𝑎 𝐿2
= ∗ ∗
𝑑𝑧 2 𝑘𝛽 𝑇𝑤 − 𝑇𝑎 𝐿2

[63]
 𝐿2 𝑑 2 𝑇 ℎ𝐿2 𝑇 − 𝑇𝑎
∗ = ( )
𝑇𝑤 − 𝑇𝑎 𝑑𝑧 2 𝑘𝛽 𝑇𝑤 − 𝑇𝑎
𝑇 − 𝑇𝑎
𝑄=
𝑇𝑤 − 𝑇𝑎
𝐿2 𝑑 2 𝑇 ℎ𝐿2
∗ 2= 𝜃
𝑇𝑤 − 𝑇𝑎 𝑑𝑧 𝑘𝛽
1 2 2
𝑇𝑤 − 𝑇𝑎 ∗ 𝑑 𝑇 = ℎ𝐿 𝜃
𝐿2 𝑑𝑧 2 𝑘𝛽
1 𝑑2𝑇 ℎ𝐿2 𝑧
∗ = 𝜃 𝜑=
𝑇𝑤 − 𝑇𝑎 𝑧 2 𝑘𝛽 𝐿
𝑑 (𝐿 )

1 𝑑2𝑇 ℎ𝐿2
∗ = 𝜃
𝑇𝑤 − 𝑇𝑎 𝑑(𝜑)2 𝑘𝛽
1 𝑑 2 (𝑇 − 𝑇𝑎 + 𝑇𝑎) ℎ𝐿2
∗ = 13 𝜃
𝑇𝑤 − 𝑇𝑎 𝑑(𝜑)2 𝑘
𝑑2 (𝑇 − 𝑇𝑎) 𝑑 2 (𝑇𝑎) ℎ𝐿2
+ = 𝜃
𝑑𝜑 2 𝑇𝑤 − 𝑇𝑎 𝑑𝜑 2 𝑇𝑤 − 𝑇𝑎 𝑘𝛽
𝑑 2 𝜃 ℎ𝐿2
= 𝜃
𝑑𝜑 2 𝑘𝛽
𝑑2𝜃
= 𝑁2𝜃
𝑑𝜑 2
𝑑2𝜃
+ 𝑁2𝜃 = 0
𝑑𝜑 2
0 = 𝑐1𝑒 𝑁𝜑 + 𝑐2𝑒 −𝑁𝜑 𝑚1 ≠ 𝑚2
𝑇 − 𝑇𝑎 𝑇𝑤 − 𝑇𝑎 𝐿
𝐶. 𝐹. 1 𝑧 = 0 𝑇 = 𝑇𝑤 𝜃 = ∴ 𝑄= = 1𝑒 = = 0
𝑇𝑤 − 𝑇𝑎 𝑇𝑤 − 𝑇𝑎 𝐿
⟹ 1. 1 = 𝑐1 + 𝑐2
𝑑 𝑇
[ = 𝑐1𝑒 𝑁𝜑 + 𝑐2𝑒 −𝑁𝜑 ]
𝑑𝑧 𝑇𝑤 − 𝑇𝑎
𝑑 𝑇 𝑞 𝑇𝑎 𝑁 𝑁
( )− ( ) = 𝑐1 ( ) 𝑒 𝑁𝜑 − 𝑐2 ( ) 𝑒 −𝑁𝜑
𝑑𝑧 𝑇𝑤 − 𝑇𝑎 𝑑𝑧 𝑇𝑤 − 𝑇𝑎 𝐿 𝐿
1 𝑑𝑇 𝑁 𝑁
= 𝑐1 ( ) 𝑒 𝑁𝜑 − 𝑐2 ( ) 𝑒 −𝑁𝜑
𝑇𝑤 − 𝑇𝑎 𝑑𝑧 𝐿 𝐿
⟹ 2. 0 = 𝑐1𝑒 𝑁 − 𝑐2𝑒 −𝑁
[64]
 𝑐1 = 1 − 𝑐2
0 = (1 − 𝑐2)𝑒 𝑁 − 𝑐2𝑒 −𝑁
0 = 𝑒 𝑁 − 𝑐2𝑒 𝑁 − 𝑐2𝑒 −𝑁
𝑐2(𝑒 𝑁 + 𝑒 −𝑁 ) = 𝑒 𝑁
𝑒𝑁
𝑐2 =
𝑒 𝑁 + 𝑒 −𝑁
𝑒𝑁
𝑐1 = 1 − 𝑁
𝑒 + 𝑒 −𝑁
𝑒 𝑁 + 𝑒 −𝑁 − 𝑒 𝑁 𝑒 −𝑁
𝑐1 = ∴ 𝑐1 =
𝑒 𝑁 + 𝑒 −𝑁 𝑒 𝑁 + 𝑒 −𝑁
𝜃 = 𝑐1𝑒 𝑁𝜑 + 𝑐2𝑒 −𝑁𝜑
𝑒 −𝑁 𝑁𝜑
𝑒 −𝑁
0= 𝑁 −𝑁
𝑒 + 𝑁 −𝑁
𝑒 −𝑁𝜑
𝑒 +𝑒 𝑒 +𝑒
𝑒 −𝑁(1−𝜑) 𝑒 𝑁(1−𝜑)
𝜃= +
𝑒 𝑁 + 𝑒 −𝑁 𝑒 𝑁 + 𝑒 −𝑁
𝑒 −𝑁(1−𝜑) + 𝑒 𝑁(1−𝜑)
𝑒 𝑁 + 𝑒 −𝑁 2
𝜃= ∗
𝑒 𝑁 + 𝑒 −𝑁 2
2
𝑐𝑜𝑠ℎ𝑁(1 − 𝜑)
𝜃=
𝑐𝑜𝑠ℎ𝑁
Aletas de enfriamiento.
Ejercicio 3.9
w

t=0.05’’

𝑄𝑎𝑙𝑒𝑡𝑎𝑠 = 𝜂ℎ𝐴𝑎𝑙𝑒𝑡𝑎𝑠 (𝑇0 − 𝑇∞ )

𝑆𝑢𝑝𝑜𝑛𝑒𝑟 𝐻 = 1 𝐹𝑡 𝑊 = 1 𝐹𝑡
𝐴 = 𝐻𝑤 − (𝑁𝑜. 𝑎𝑙𝑒𝑡𝑎𝑠)(𝑡)(𝑤)
𝐻
𝑁𝑜. 𝑎𝑙𝑒𝑡𝑎𝑠 =
𝑒𝑠𝑝𝑎𝑐𝑖𝑎𝑚𝑖𝑒𝑛𝑡𝑜 𝑒𝑛𝑡𝑟𝑒 𝑎𝑙𝑒𝑡𝑎𝑠

[65]
 1 𝐹𝑡
𝑁𝑜. 𝑎𝑙𝑒𝑡𝑎𝑠 = = 24
0.5
12 𝐹𝑡
𝐴0 = 𝐻𝑤 − 𝑁𝑜. 𝑎𝑙𝑒𝑡𝑎𝑠(𝑡)(𝑤)
0.05
𝐴0 = 1 𝐹𝑡(1𝐹𝑡) − 24 ( 𝑓𝑡) (1 𝐹𝑡)
12
𝐴0 = 0.9 𝐹𝑡 2
𝐴𝑓 = 2(𝑁𝑜. 𝑎𝑙𝑒𝑡𝑎𝑠)(𝐿)(𝑤) + 𝑁𝑜. 𝑎𝑙𝑒𝑡𝑎𝑠(𝑡)(𝑤)
1 0.05
𝐴𝑓 = 2(24) ( ) (1) + 24 ( ) (1)
12 12
𝐴𝑓 = 4.1 𝐹𝑡 2

𝐸𝑓𝑒𝑐𝑡𝑖𝑣𝑖𝑑𝑎𝑑 = 𝜂

𝑡 3 2ℎ
(𝐿 + )2 √
2 𝑘𝑡𝐿

3 𝐵𝑡𝑢
1 0.05 2 2 (2 )
( + ) √ ℎ − 𝐹𝑡 2 − °𝐹
12 2(12) 𝐵𝑡𝑢 0.05 1
24.8
ℎ − 𝐹𝑡 − °𝐹 12 12
𝐵𝑡𝑢
𝐾𝑎𝑐𝑒𝑟𝑜 = 24.8
ℎ − 𝐹𝑡 − °𝐹
𝐴𝑑𝑖𝑚𝑒𝑛𝑐𝑖𝑜𝑛𝑎𝑙 = 0.5380 𝑎𝑖𝑟𝑒
𝜂𝑎𝑖𝑟𝑒 = 0.9
𝐻2 𝑂

𝑡 3 2ℎ
(𝐿 + )2 √
2 𝑘𝑡𝐿

3 𝐵𝑡𝑢
1 0.05 2 2 (45 2 − °𝐹 )
( + ) √ ℎ − 𝐹𝑡 = 2.5521
12 2(12) 𝐵𝑡𝑢 0.05 1
24.8
ℎ − 𝐹𝑡 − °𝐹 12 12
𝑔𝑟𝑎𝑓𝑖𝑐𝑎 17.11 𝜂𝐻2𝑂 = 0.4
𝑎𝑖𝑟𝑒
𝑄 = 𝐴𝑜 ℎ(∆𝑇) + 𝐴𝑓 ℎ𝜂(∆𝑇)

𝑄 = 𝐴𝑜 ℎ(𝑇𝑤 − 𝑇∞ ) + 𝐴𝑓 ℎ𝜂(𝑇𝑤 − 𝑇∞ )
[66]
 𝑄 = 0.9(2)(𝑇𝑤 − 𝑇∞ ) + 4.1(2)(0.9)(𝑇𝑤 − 𝑇∞ )
𝑄 = 9.18(𝑇𝑤 − 𝑇∞ )
𝑄 = ℎ(𝑇𝑤 − 𝑇∞ )
𝑞 = 𝐴ℎ(𝑇𝑤 − 𝑇∞ ) ∴ 𝐴 = 1
𝑞 = ℎ(𝑇𝑤 − 𝑇∞ )
𝐻2 𝑂
𝑄 = 𝐴𝑜 ℎ(𝑇𝑤 − 𝑇∞ ) + 𝐴𝑓 ℎ𝜂(𝑇𝑤 − 𝑇∞ )

𝑄 = 0.9(45)(𝑇𝑤 − 𝑇∞ ) + 4.1(0.4)(45)(𝑇𝑤 − 𝑇∞ )
𝑄 = 114.3(𝑇𝑤 − 𝑇∞ ) ∴ 𝐴 = 1
𝑞 = 114.3(𝑇𝑤 − 𝑇∞ )
(𝑇𝑤 − 𝑇∞ )
𝑞=
1 𝑥1 − 𝑥0 1
+ +
ℎ𝑖 𝑘 01 ℎ𝑒
(𝑇𝑤 − 𝑇∞ )
𝑞=
1 𝑥1 − 𝑥0 1
+ +
ℎ𝑖 𝑘 01 ℎ𝑒
𝑞 = 1.9148(𝑇𝑤 − 𝑇∞ ) = ℎ(𝑇𝑤 − 𝑇∞ )
𝐼𝑛𝑐𝑟𝑒𝑚𝑒𝑛𝑡𝑜𝑠 − 𝐴𝑙𝑒𝑡𝑎𝑠 𝑒𝑛 𝑒𝑙 𝑎𝑖𝑟𝑒
(𝑇𝑤 − 𝑇∞ )
𝑞= = 7.6245(𝑇𝑤 − 𝑇∞ )
1 1
9.18 + 45
7.6245(𝑇𝑤 − 𝑇∞ ) − 1.9148(𝑇𝑤 − 𝑇∞ )
𝐼𝑛𝑐𝑟𝑒𝑚𝑒𝑛𝑡𝑜 𝑑𝑒 𝑐𝑎𝑙𝑜𝑟: ∗ 100 = 297.72%
1.9148(𝑇𝑤 − 𝑇∞ )
𝐴𝑙𝑒𝑡𝑎𝑠 𝑒𝑛 𝑒𝑙 𝐻2 𝑂
(𝑇𝑤 − 𝑇∞ )
𝑞= = 1.96(𝑇𝑤 − 𝑇∞ )
1 1
+
2 114.3
1.96(𝑇𝑤 − 𝑇∞ ) − 1.9148(𝑇𝑤 − 𝑇∞ )
𝐼𝑛𝑐𝑟𝑒𝑚𝑒𝑛𝑡𝑜 𝑑𝑒 𝑐𝑎𝑙𝑜𝑟: ∗ 100 = 2.36%
1.9148(𝑇𝑤 − 𝑇∞ )
𝐴𝑙𝑒𝑡𝑎𝑠 𝑒𝑛 𝑒𝑙 𝐻2 𝑂 𝑦 𝑎𝑖𝑟𝑒
(𝑇𝑤 − 𝑇∞ )
𝑞= = 1.96(𝑇𝑤 − 𝑇∞ )
1 1
+
9.18 114.3

[67]
 8.49(𝑇𝑤 − 𝑇∞ ) − 1.9148(𝑇𝑤 − 𝑇∞ )
𝐼𝑛𝑐𝑟𝑒𝑚𝑒𝑛𝑡𝑜 𝑑𝑒 𝑐𝑎𝑙𝑜𝑟: ∗ 100 = 343.38%
1.9148(𝑇𝑤 − 𝑇∞ )

Aletas de enfriamiento - aletas circulares

Ejercicio 3.10

Para disipar calor, se colocan unas aletas de aluminio de 1.5 cm de ancho y 1.0 mm
de espesor sobre un tubo de 2.5 cm de diámetro. La temperatura de superficie del
tubo es de 170°C, y la temperatura del fluido del ambiente es 25 °C. Calcular la
perdida de calor por aleta para un coeficiente de transferencia de calor de 130
W/m2-°C. Suponga que k=200 W/m-°C.

r2

2.5
𝑟1 = = 1.25 𝑐𝑚
2
𝑟2 = 𝑟1 + 𝐿 = 1.25 − 1.5 𝑟2 = 2.75 𝑐𝑚 𝐶𝑎𝑙𝑐𝑢𝑙𝑎𝑟 𝑙𝑎 𝑒𝑓𝑒𝑐𝑡𝑖𝑣𝑖𝑑𝑎𝑑
𝑎𝑙𝑒𝑡𝑎 𝑐𝑖𝑟𝑐𝑢𝑙𝑎𝑟
3
𝑡 2 2ℎ
(𝑟2 + − 𝑟0) √
2 𝑘𝑡(𝑟2 − 𝑟0)

𝑤𝑎𝑡𝑡 𝑤
ℎ = 130 𝑘 = 200
𝑚2 − °𝐶 𝑚 − °𝐶
3 𝑤𝑎𝑡𝑡
1𝑋103 2 2 (130 𝑚2 − °𝐶 )
(0.0275 + − 0.0125) √ 𝑤 = 0.568
2 200 (1𝑋10 − 3)(0.0275 − 0.0125)𝑚
𝑚 − °𝐶
𝑟2 2.75
𝑢𝑠𝑎𝑟 𝑔𝑟𝑎𝑓𝑖𝑐𝑎 17.11 = = 2.2 𝜂 = 0.8
𝑟1 1.25
𝑐𝑎𝑙𝑐𝑢𝑙𝑎𝑟 𝑝𝑒𝑟𝑑𝑖𝑑𝑎 𝑑𝑒 𝑐𝑎𝑙𝑜𝑟
𝑡 2
𝑄 = 2𝜋 [(𝑟2 + ) − 𝑟12 ] ℎ 𝜂(𝑇𝑜 − 𝑇∞)
2

[68]
 1𝑋10 − 3 2
𝑄 = 2𝜋 [(0.0275 + ) − 0.01252 ] (130)(0.8)(170 − 25)
2

𝑄 = 59.4795 𝑤𝑎𝑡𝑡𝑠
Aletas de enfriamiento
Ejercicio 3.11

Una aleta de aluminio con k=200 W/m-°C de 3mm de espesor y 7.5 cm de longitud
sobresalen una pared igual a la que se muestra en la figura. La base se mantiene a
300°C, y la temperatura del ambiente es de 50°C, con un coeficiente de transferencia
de calor de 10 W/m2-°C. Calcular la perdida de calor de la aleta por unidad de
anchura de material.

Z-ancho
7.5 cm
𝑤
𝑘 = 200
𝑚 − °𝐶
cm
3

𝑤
ℎ = 10 2
𝑚 − °𝐶
𝑇𝑤 = 300 °𝐶
𝑇𝑏 = 50 °𝐶

𝑇 = 𝑇0 𝑥=0
𝑎) { }
𝑇 = 𝑇𝐿 𝑥=𝐿
𝑇 = 𝑇0 𝑥=0
𝑏) {𝑑𝑇 }
=0 𝑥=𝐿
𝑑𝑥
𝑇 = 𝑇0 𝑥=0
𝑐) { 𝑑𝑇 }
−𝑘 = ℎ(𝑇 − 𝑇∞ ) 𝑥=𝐿
𝑑𝑥
𝑄 = 𝑘𝐴𝑚𝜃0 𝑡𝑎𝑛ℎ𝑚𝐿
𝜃0 = 𝑇0 − 𝑇∞
𝜃0 = 𝑇𝑤 − 𝑇𝑏 ∴ 𝑇𝑏 = 𝑇∞
𝐴 = 𝑧(3𝑥10−3 𝑚)
ℎ𝑃
𝑚2 =
𝑘𝐴

ℎ𝑃
𝑚=√ 𝑃(𝑝𝑒𝑟𝑖𝑚𝑒𝑡𝑟𝑜) = 2𝑧 + 𝑧𝑡
𝑘𝐴

[69]
 ℎ(2𝑧 + 2𝑡)
𝑚=√
𝑘𝑧(3𝑥10−3 )

𝑧 ≫≫> 𝑡

2ℎ𝑧
𝑚=√
𝑘𝑧(3𝑥10−3 )

𝑊
(2)10
𝑚=√ 𝑚2 − °𝐶
𝑊
200 𝑚 − °𝐶 (3𝑥10−3 )

5.77535
𝑚=
𝑚
𝑤 5.7735 5.77355
𝑄 = 200 ( ) (300 − 50)°𝐶(𝑧(3𝑥10−3 𝑚))𝑡𝑎𝑛ℎ ( (0.075𝑚))
𝑚 − °𝐶 𝑚 𝑚
𝑄 𝑤
= 353.196
𝑧 𝑚
Aletas de enfriamiento
Ejercicio 3.12

Una varilla cilíndrica de 3 cm de diámetro está insertada parcialmente en un horno y

uno de sus extremos está expuesto al aire circundante que está a 300 K. las
temperaturas en dos de sus puntos, separados 7.6 cm. Son de 399 K y 365 K. si el
coeficiente de transferencia de calor convectivo es 17 W/m2K, determine la
conductividad térmica del material del que está hecha la varilla.
Ejercicio # – Libro 17.18
X=0 X=L

T0 TL
D=3cm

Tb = T = 300 °K
7.6cm

𝑇𝐿 − 𝑇∞
= −𝑒 𝑚𝑥
𝑇0 − 𝑇∞
𝑇0 − 𝑇∞
= 𝑒 𝑚𝑥 𝑥 = 𝐿
𝑇𝐿 − 𝑇∞

[70]
 𝑇0 − 𝑇∞
= 𝑒 𝑚𝐿
𝑇𝐿 − 𝑇∞
399 − 300
= 𝑒 𝑚(0.076)
365 − 300
99
= 𝑒 𝑚(0.076)
65
99
𝑙𝑛 = 𝑚(0.076)
65
0.7207
𝑚=
0.076
5.5342
𝑚=
𝑚
ℎ𝑃
𝑚2 =
𝑘𝐴

ℎ(𝜋𝐷)
𝑚=√
𝜋(0.03)2
𝑘( )
4

𝑤
5.5342 17 2 (𝜋(0.03𝑚))
=√ 𝑚 −𝐾
𝑚 𝜋(0.03)2
𝑘( )
4
𝑤
𝐾 = 74.0078
𝑚 − °𝐾
Transferencia de calor por radiación
Ejercicio 3.13

Una caja de lámina de acero de forma cúbica de 0.70 m por lado tiene una emitancia
de 0.7. La caja contiene equipo electrónico que disipa 1,200 W de energía. Si la región
circundante se supone negra, que se encuentra a 280 K y se considera que la tapa y
las paredes de la caja radian uniformemente, ¿cuál será la temperatura de la
superficie de la caja?
0.7 m

T1
[71]
 𝐸1 − 𝐸2 𝜎(𝑇14 − 𝑇24 )
𝑄12 = =
∑𝑅 ∑𝑅
𝜎(𝑇14 − 𝑇24 ) 𝜎(𝑇14 − 𝑇24 )
𝑄12 = =
1 −∈1 1 1 −∈1 1
+ +
𝐴1 ∈1 𝐴1 𝐹12 𝐴1 ∈1 𝐴1 𝐹12
𝑄 = 𝐴1 ∈1 𝜎(𝑇14 − 𝑇24 )
𝑄
=∈1 𝜎(𝑇14 − 𝑇24 )
𝐴1
𝑄 1200𝑤 𝑊
= = 489.7959
𝐴1 (0.7)2 (5) 𝑚2
𝑊 𝑤
489.7959 2
= (0.7)(5.686𝑥10−8 2 )(𝑇14 − 2804 )
𝑚 𝑚 − 𝑘4
𝑇1 = 368.67 °𝐾
Transferencia de calor por radiación
Ejercicio 3.14

La base circular del recipiente cilíndrico de la figura se puede considerar como una
superficie rerradiante, las paredes del cilindro tienen una emitancia efectiva de 0.80 y
se mantiene a 540 °F. La parte superior del recipiente permanece abierta a la región
circundante se mantiene a 40° F. ¿Cuál es la rapidez neta de transferencia de
radiación a la región circundante?

T3

6 ft T2

T1
𝐸1 − 𝐸2 𝜎(𝑇24 − 𝑇34 )
𝑄23 = =
∑𝑅 ∑𝑅
𝜎(𝑇24 − 𝑇34 )
𝑄23 =
1 −∈2 1
𝐴2 ∈2 + 1
𝐴2 𝐹23 + 1 1
𝐴2 𝐹21 + 𝐴3 𝐹31
𝐴𝑟𝑒𝑎 𝑑𝑒 𝑢𝑛 𝑐𝑖𝑙𝑖𝑛𝑑𝑟𝑜
𝐴2 = 2𝜋𝑟𝐿
𝐴2 = 2𝜋(6)(6) = 226.19 𝑓𝑡 2

[72]
 𝐴3 = 𝜋(𝑟)2 = 𝜋(6)2 = 113.097 𝑓𝑡 2
𝐴3 = 𝐴1
12
𝐺𝑟𝑎𝑓𝑖𝑐𝑎 23.14 𝑅𝑒𝑙𝑎𝑐𝑖ó𝑛 = =2
6
𝐴3 𝐹31 = 𝐴1 𝐹13
𝐹13 = 𝐹31 = 0.38
𝐹12 + 𝐹13 = 1 𝐹12 = 1 − 0.38 = 0.62
𝐴1 𝐹12 = 𝐴2 𝐹21
113.098(0.62) = 226.19 𝐹21 𝐹21 = 0.31
𝐴2 𝐹23 = 𝐴3 𝐹32
𝐴2 𝐹21 = 𝐴1 𝐹12

_________________________________
𝐴2 𝐹23 = 𝐴2 𝐹21
226.19𝐹23 = 226.19(0.31)
𝐹23 = 0.31
𝜎(𝑇24 − 𝑇34 )
𝑄23 =
1 −∈2 1
𝐴2 ∈2 + 1
𝐴2 𝐹23 + 1 1
+
𝐴2 𝐹21 𝐴3 𝐹31
𝐵𝑡𝑢
(0.1716𝑥10−8 )(10004 − 5004 )
ℎ − 𝑓𝑡 2 − °𝑅 4
𝑄23 =
1 − 0.8 1
+ 1
226.19(0.8)
226.19(0.31) + 1 1
+
226.19(0.8) 113.097(0.38)
𝐵𝑡𝑢
𝑄23 = 140 467.9482
ℎ
Transferencia de calor por radiación
Ejercicio 3.15

Una habitación que mide 12 ft x 20 ft x 8 ft tiene una temperatura de 85°F en el piso y

65° F en el techo. Suponiendo que las paredes son rerradiantes y todas las superficies
tienen una emitancia de 0.8, determine el intercambio neto de energía entre el piso y
el techo.
T= 65°F

[73]
 8 ft

20 ft
12 ft

𝐸1 − 𝐸2 𝜎(𝑇14 − 𝑇24 )
𝑄12 = =
∑𝑅 ∑𝑅
𝜎(𝑇14 − 𝑇24 )
𝑄23 =
1 −∈1 1 1 −∈
+ 𝐴 ∈2
𝐴1 ∈1 + 1 2 2
𝐴1 𝐹12 + 1 1
𝐴1 𝐹13 + 𝐴2 𝐹23
𝐴1 = 𝐴2 𝐴1 = 20(12) = 240 𝑓𝑡 2 𝐸𝑛𝑐𝑜𝑛𝑡𝑟𝑎𝑟 𝑒𝑙 𝑓𝑎𝑐𝑡𝑜𝑟 𝐹12 = 0.4 𝐹21 = 0.4
𝐹12 + 𝐹13 = 1 𝐹13 = 1 − 0.4 = 0.6
𝐴2 𝐹23 = 𝐴3 𝐹32
𝐴2 𝐹13 = 𝐴3 𝐹31

_________________________________
𝐴2 𝐹23 = 𝐴1 𝐹13
𝐹23 = 0.6
𝐵𝑡𝑢
(0.1714𝑥10−8 )(5444 − 524.674 )
ℎ − 𝑓𝑡 2 − °𝑅 4
𝑄23 =
1 − 0.8 1
+ 1
240(0.8)
240(0.4) + 1 1
+
240(0.6) 240(0.6)
𝐵𝑡𝑢
𝑄23 = 2 609.06
ℎ
Transferencia de calor por radiación
Ejercicio 3.16

Un conducto circular de 2 ft de longitud tiene, en su centro, un termopar cuya área

es de 0.3 in2. Las paredes del conducto se encuentran a 200 °F y el termopar indica
310°F. Suponiendo que el coeficiente de transferencia de calor convectivo entre el
termopar y el gas del conducto es de 30 Btu/h-ft2-°F, calcule la temperatura real del
gas. La emitancia de las paredes del conducto es de 0.8 y la del termopar de 0.6.

[74]
 𝐸1 − 𝐸2 𝜎(𝑇𝑢𝑛𝑖𝑜𝑛4 − 𝑇24 )
𝑄𝑅 = =
∑𝑅 ∑𝑅
𝜎(𝑇𝑢𝑛𝑖𝑜𝑛4 − 𝑇24 )
𝑄𝑅 =
1 −∈𝑢𝑛𝑖𝑜𝑛 1
+
𝐴𝑢𝑛𝑖𝑜𝑛 ∈𝑢𝑛𝑖𝑜𝑛 𝐴𝑢𝑛𝑖𝑜𝑛 𝐹12
𝑄𝑅 𝜎(𝑇𝑢𝑛𝑖𝑜𝑛4 − 𝑇24 )
= = 𝑞𝑅
𝐴𝑢𝑛𝑖𝑜𝑛 1 −∈𝑢𝑛𝑖𝑜𝑛
+ 1
∈𝑢𝑛𝑖𝑜𝑛
𝜎(𝑇𝑢𝑛𝑖𝑜𝑛4 − 𝑇24 )
𝑞𝑅 =
1 −∈𝑢𝑛𝑖𝑜𝑛 +∈𝑢𝑛𝑖𝑜𝑛
∈𝑢𝑛𝑖𝑜𝑛
𝑞𝑅 =∈𝑢𝑛𝑖𝑜𝑛 𝜎(𝑇𝑢𝑛𝑖𝑜𝑛4 − 𝑇24 )
𝑞𝑐 = ℎ∆𝑇 = ℎ(𝑇𝑔𝑎𝑠 − 𝑇𝑢𝑛𝑖𝑜𝑛)
ℎ(𝑇𝑔𝑎𝑠 − 𝑇𝑢𝑛𝑖𝑜𝑛) =∈𝑢𝑛𝑖𝑜𝑛 𝜎(𝑇𝑢𝑛𝑖𝑜𝑛4 − 𝑇24 )
𝐵𝑡𝑢
30(𝑇𝑔𝑎𝑠 − 310°𝐹) = 0.6(0.1714𝑥10−8 2 4
)(769.674 − 659.674 )
ℎ − 𝑓𝑡 − °𝑅
𝑇𝑔𝑎𝑠 = 315.53 °𝐹

[75]
UNIDAD
CUATRO

[76]
EJERCICIOS DE CLASE
Ejercicio 4.1
Determine los valores de las siguientes difusividades de los gases.
CO2 / aire a 310°K y 1.5X105 Pa.
1
3 1 1 2
0.001858𝑇 2 (𝑀𝐴 𝑀𝐵 )
+
𝐷𝐴𝐵 = 2
𝑃𝜎𝐴𝐵 Ω𝐷
𝑔 𝜀𝐴
𝐴 = 𝐶𝑂2 ; 𝑀𝐴 = 44 ; = 190 𝐾 ; 𝜎 = 3.996 𝐴̇
𝑚𝑜𝑙 𝑘
𝑔 𝜀𝐴
𝐵 = 𝑎𝑖𝑟𝑒 ; 𝑀𝐵 = 29 ; = 97 𝐾 ; 𝜎 = 3.617 𝐴̇
𝑚𝑜𝑙 𝑘
𝑃 = 1.5𝑋105 𝑃𝑎 = 1.48.3 𝑎𝑡𝑚
𝐴𝑝𝑒𝑛𝑑𝑖𝑐𝑒 𝐾
𝜎𝐴 + 𝜎𝐵 3.996 + 3.617
𝜎𝐴𝐵 = = = 3.8065
2 2
𝜎𝐴𝐵 2 = 14.48
𝜀𝐴𝐵 𝜀𝐴 𝜀𝐴
= √ ( ) = √190(97) = 135.7572
𝑘 𝑘 𝑘
𝜀𝐴𝐵 1
(𝑇) = (310) = 2.2834
𝑘 135.7572
Ω𝐷 = 𝑥
2.20 − 1.041
2.2834 − 𝑥
2.30 − 1.026
Ω𝐷 = 1.02844
1
3 1 1 2
0.001858𝑇 2 ( + 29)
𝐷𝐴𝐵 = 44
(1.4803 𝑎𝑡𝑚)(14.48 𝐴)(1.02844)
𝑐𝑚2
𝐷𝐴𝐵 = 0.110038
𝑠

[77]
Ejercicio 4.2
Determinas el valor de la siguiente difusividad. Etanol/aire a 235 °K y 2.0X105 PA.
𝑇 = 325 °𝐾
𝑃 = 2𝑋105 𝑃𝑎 = 1.9738 𝑎𝑡𝑚
𝑇1 = 25°𝐶 = 298°𝐾
1𝑒𝑟 𝑚𝑒𝑡𝑜𝑑𝑜 { }
𝑃1 = 1 𝑎𝑡𝑚
1
3 1 1 2
0.001858𝑇 2 ( 𝑀𝐴 + 𝑀𝐵 )
𝐷𝐴𝐵 = 2
𝑃𝜎𝐴𝐵 Ω𝐷
𝑔 𝜀𝐴
𝐴 = 𝐸𝑡𝑎𝑛𝑜𝑙 ; 𝑀𝐴 = 46 ; = 391 𝐾 ; 𝜎 = 4.455 𝐴̇
𝑚𝑜𝑙 𝑘
𝑔 𝜀𝐴
𝐵 = 𝑎𝑖𝑟𝑒 ; 𝑀𝐵 = 29 ; = 97 𝐾 ; 𝜎 = 3.617 𝐴̇
𝑚𝑜𝑙 𝑘
𝐴𝑝𝑒𝑛𝑑𝑖𝑐𝑒 𝐾
𝜎𝐴 + 𝜎𝐵 4.435 + 3.617
𝜎𝐴𝐵 = = = 4.036
2 2
𝜎𝐴𝐵 2 = 16.2892
𝜀𝐴𝐵 𝜀𝐴 𝜀𝐴
= √ ( ) = √391(97) = 194.7485
𝑘 𝑘 𝑘
𝜀𝐴𝐵 1
(𝑇) = (298) = 1.5301
𝑘 194.7485
Ω𝐷 = 𝑥
1.5 − 1.198
1.5301 − 𝑥
1.55 − 1.182
Ω𝐷 = 1.1803
1
3 1 1 2
0.001858𝑇 2 ( + 29)
𝐷𝐴𝐵 = 46
(1.9738 𝑎𝑡𝑚)(16.2892 𝐴)(1.1803)
𝑐𝑚2
𝐷𝐴𝐵 = 0.1178
𝑠

[78]
 𝑇1 = 298°𝐾
2𝑑𝑜 𝑚𝑒𝑡𝑜𝑑𝑜 { }
𝑃1 = 1 𝑎𝑡𝑚
3
𝑃1 𝑇1 2 Ω𝐷 | 𝑇1
𝐷𝐴𝐵𝑇2𝑃2 = 𝐷𝐴𝐵𝑇1𝑃1 ( ) ( )
𝑃2 𝑇2 Ω𝐷 | 𝑇2
1
Ω𝐷 | 𝑇2 = (325) = 1.6688
194.7485
𝐴𝑝𝑒𝑛𝑑𝑖𝑐𝑒 𝐾
Ω𝐷 = 𝑥
1.65 − 1.153
1.6888 − 𝑥
1.70 − 1.140
Ω𝐷 = 1.1481
3
1 325 2 1.1883
𝐷𝐴𝐵𝑇2𝑃2 = 0.1178 ( )( )
1.9738 298 1.1481
𝑐𝑚2
𝐷𝐴𝐵𝑇2𝑃2 = 0.0703
𝑠
Ejercicio 4.3

Las bolas contra polilla se fabrican a menudo de naftalina. Cuando la naftalina se

vaporiza, se difunde en el aire atmosférico circundante. Calcule la difusividad del
sistema de naftalina y aire a 1 atm y 20°C, usando la ccuaci6n (24-33) de Wirschfelder,
Bird y Spotz. Calcule el valor a la presi6n y temperatura estipuladas, usando el valor
que aparece en el Apéndice J. Las propiedades críticas de la naftalina son:
Vc = 3.1847 ml/g
Tc = 469°C
Pc = 29 792 torr
1
3 1 1 2
0.001858𝑇 2 (𝑀𝐴 + 𝑀𝐵 )
𝐷𝐴𝐵 = 2
𝑃𝜎𝐴𝐵 Ω𝐷
𝑔 𝜀𝐴
𝐴 = 𝑁𝑎𝑓𝑡𝑎𝑙𝑒𝑛𝑜 ; 𝑀𝐴 = 28 ; = 571.34 𝐾 ; 𝜎 = 6.2385 𝐴̇
𝑚𝑜𝑙 𝑘
𝑔 𝜀𝐴
𝐵 = 𝑎𝑖𝑟𝑒 ; 𝑀𝐵 = 29 ; = 97 𝐾 ; 𝜎 = 3.617 𝐴̇
𝑚𝑜𝑙 𝑘

[79]
 1 1
𝜎 = 0.841𝑉𝑐 3 = 0.841(408.1845)3 = 6.2385
128.1702𝑔
𝑚𝑙 1𝑐𝑚3 𝑚𝑜𝑙 𝑐𝑚3
𝑉𝑐 = 3.1847 = 408.1845
𝑔 1𝑚𝑙 1 𝑚𝑜𝑙 𝑚𝑜𝑙
𝜀𝐴
= 0.77 𝑇𝑐 = 0.77(742) = 571.34
𝑘
𝑇𝑐 = 742
𝜎𝐴 + 𝜎𝐵 6.2385 + 3.617
𝜎𝐴𝐵 = = = 4.9277
2 2
𝜎𝐴𝐵 2 = 24.2827
𝜀𝐴𝐵 𝜀𝐴 𝜀𝐴
= √ ( ) = √571.34(97) = 235.41
𝑘 𝑘 𝑘
𝜀𝐴𝐵 1
(𝑇) = (293) = 1.2446
𝑘 235.41
Ω𝐷 = 𝑥
1.20 − 1.320
1.2446 − 𝑥
1.25 − 1.296
Ω𝐷 = 1.1803
1
3 1 1 2
0.001858𝑇 2 (128 29)
+
𝐷𝐴𝐵 =
(1 𝑎𝑡𝑚)(24.2827 𝐴)(1.2985)
𝑐𝑚2
𝐷𝐴𝐵 = 0.0607
𝑠

[80]
Ejercicio 4.4

Determinar la difusividad de CO2 a través de una mezcla de gases que tiene la

siguiente composición:
O2 7%
CO 10%
CO2 15%
N2 68%
1
𝐷1−𝑚𝑒𝑧𝑐𝑙𝑎 =
𝑦′ 2
𝑦′ 3
𝑦′
+𝐷 +𝐷 4 +⋯
𝐷1−2 1−3 1−4

𝑇 = 273 °𝐾 𝑃 = 1.5𝑋105 𝑃𝑎
1
𝐷𝐶𝑂2−𝑚𝑒𝑧𝑐𝑙𝑎 =
𝑦′ 2
𝑦 ′3 𝑦 ′4
+𝐷 +𝐷
𝐷𝐶𝑂2−𝑂2 𝐶𝑂2−𝐶𝑂 𝐶𝑂2−𝑁2

𝐷𝐶𝑂2−𝑂2 =?
𝐷𝐶𝑂2−𝐶𝑂 =?
𝐷𝐶𝑂2−𝑁2 =?
1
3 1 1 2
0.001858𝑇 2 ( 𝑀𝐴 + 𝑀𝐵 )
𝐷𝐴𝐵 = 2
𝑃𝜎𝐴𝐵 Ω𝐷
𝑔 𝜀𝐴
𝐴 = 𝐶𝑂2 ; 𝑀𝐴 = 44 ; = 190 𝐾 ; 𝜎 = 3.996 𝐴̇
𝑚𝑜𝑙 𝑘
𝑔 𝜀𝐴
𝐵 = 𝑂2 ; 𝑀𝐵 = 32 ; = 113 𝐾 ; 𝜎 = 3.433 𝐴̇
𝑚𝑜𝑙 𝑘
𝐴𝑝𝑒𝑛𝑑𝑖𝑐𝑒 𝐾
𝜎𝐴 + 𝜎𝐵 3.996 + 3.433
𝜎𝐴𝐵 = = = 3.7145
2 2
𝜎𝐴𝐵 2 = 13.7975
𝜀𝐴𝐵 𝜀𝐴 𝜀𝐴
= √ ( ) = √190(113) = 146.52
𝑘 𝑘 𝑘
𝜀𝐴𝐵 1
(𝑇) = (273) = 1.8631
𝑘 146.52
Ω𝐷 = 𝑥

[81]
 1.85 − 1.105
1.8631 − 𝑥
1.90 − 1.094
Ω𝐷 = 1.102
1
3 1 1 2
0.001858𝑇 2 ( + 32)
𝐷𝐴𝐵 = 44
(1 𝑎𝑡𝑚)(13.7975 𝐴)(1.102)
𝑐𝑚2
𝐷𝐴𝐵 = 0.0862
𝑠
𝑔 𝜀𝐴
𝐴 = 𝐶𝑂2 ; 𝑀𝐴 = 44 ; = 190 𝐾 ; 𝜎 = 3.996 𝐴̇
𝑚𝑜𝑙 𝑘
𝑔 𝜀𝐴
𝐵 = 𝐶𝑂 ; 𝑀𝐵 = 28 ; = 110 𝐾 ; 𝜎 = 3.590 𝐴̇
𝑚𝑜𝑙 𝑘
𝐴𝑝𝑒𝑛𝑑𝑖𝑐𝑒 𝐾
𝜎𝐴 + 𝜎𝐵 3.996 + 3.590
𝜎𝐴𝐵 = = = 3.793
2 2
𝜎𝐴𝐵 2 = 14.3868
𝜀𝐴𝐵 𝜀𝐴 𝜀𝐴
= √ ( ) = √190(110) = 144.5683
𝑘 𝑘 𝑘
𝜀𝐴𝐵 1
(𝑇) = (273) = 1.8883
𝑘 144.5683
Ω𝐷 = 𝑥
1.85 − 1.105
1.8883 − 𝑥
1.90 − 1.094
Ω𝐷 = 1.0965
1
3 1 1 2
0.001858𝑇 2 (
44 28)
+
𝐷𝐴𝐵 =
(1 𝑎𝑡𝑚)(14.3868 𝐴)(1.0965)
𝑐𝑚2
𝐷𝐴𝐵 = 0.0867
𝑠
𝑔 𝜀𝐴
𝐴 = 𝐶𝑂2 ; 𝑀𝐴 = 44 ; = 190 𝐾 ; 𝜎 = 3.996 𝐴̇
𝑚𝑜𝑙 𝑘

[82]
 𝑔 𝜀𝐴
𝐵 = 𝑁2 ; 𝑀𝐵 = 28 ; = 91.5 𝐾 ; 𝜎 = 3.681 𝐴̇
𝑚𝑜𝑙 𝑘
𝐴𝑝𝑒𝑛𝑑𝑖𝑐𝑒 𝐾
𝜎𝐴 + 𝜎𝐵 3.996 + 3.681
𝜎𝐴𝐵 = = = 3.8385
2 2
𝜎𝐴𝐵 2 = 14.7340
𝜀𝐴𝐵 𝜀𝐴 𝜀𝐴
= √ ( ) = √190(91.5) = 131.8521
𝑘 𝑘 𝑘
𝜀𝐴𝐵 1
(𝑇) = (273) = 2.0705
𝑘 131.8521
Ω𝐷 = 𝑥
2.00 − 1.075
2.0705 − 𝑥
2.10 − 1.057
Ω𝐷 = 1.06231
1
3 1 1 2
0.001858𝑇 2 (
44 28)
+
𝐷𝐴𝐵 =
(1 𝑎𝑡𝑚)(14.7340 𝐴)(1.06231)
𝑐𝑚2
𝐷𝐴𝐵 = 0.0874
𝑠
O2 7% 7% 0.0823 Y’2
CO 10% 10% 0.1173 Y’3
CO2 15%
N2 68% 68% 0.8 Y’4

1
𝐷𝐶𝑂2−𝑚𝑒𝑧𝑐𝑙𝑎 =
0.0823 0.0176 0.8
0.0862 + 0.0864 + 0.0874
𝑐𝑚2
𝐷𝐶𝑂2−𝑚𝑒𝑧𝑐𝑙𝑎 = 0.0871
𝑠

[83]
Ejercicio 4.5
Estimar la difusividad de N-Butanol / Agua a 15 °C.
1
𝐷𝐴𝐵 𝜇𝐵 7.4𝑋10−8 (Φ𝐵 𝑀𝐵 )2
=
𝑇 𝑉𝐴6
𝑇 = 15°𝐶 = 288°𝐾
Φ𝐵 = 2.6
𝑐𝑚3
𝑉𝐴 = 103.6
𝑚𝑜𝑙
𝐶4 𝐻10 𝑂 = 4(14.8) + 10(3.7) + 1(7.4) = 103.6
𝐴𝑝𝑒𝑛𝑑𝑖𝑐𝑒 𝐼 𝜇𝐵
1.794𝑋10−6 𝑃𝑎 − 𝑠 / 273 °𝐾
𝜇𝐵 / 288 °𝐾
0.993𝑋10−6 𝑃𝑎 − 𝑠 / 295 °𝐾
𝜇𝐵 = 1193.25𝑋10−6 𝑃𝑎 − 𝑠
𝜇𝐵 = 1193.25𝑋10−6 𝑃𝑜𝑖𝑠𝑒𝑠
𝑁
𝜇𝐵 = 1193.25𝑋10−6 −𝑠
𝑚2
𝐾𝑔 − 𝑚
2 −𝑠
𝜇𝐵 = 1193.25𝑋10−6 𝑠 2
𝑚
𝐾𝑔
𝜇𝐵 = 1193.25𝑋10−6
𝑚−𝑠
𝑔𝑟
𝜇𝐵 = 1193.25𝑋10−5
𝑐𝑚 − 𝑠
𝜇𝐵 = 1193.25𝑋10−2 𝑐𝑝
1
𝐷𝐴𝐵 (1.19325) 7.4𝑋10−8 (2.6(18))2
=
288 103.66
𝑐𝑚2
𝐷𝐴𝐵 = 7.5474𝑋10−6
𝑠

[84]
Ejercicio 4.6

Con frecuencia los suministros municipales de agua se tratan con cloro como una de
las etapas en el tratamiento de aguas residuales. Estimar el coeficiente de difusión
en líquidos del cloro en una solución infinita de agua a 289 °K.
Utilizando:
a. Ecuacion Wike-chang
1
𝐷𝐴𝐵 𝜇𝐵 7.4𝑋10−8 (Φ𝐵 𝑀𝐵 )2
=
𝑇 𝑉𝐴6
𝑐𝑚3
𝐴 = 𝐶𝑙2 𝐵 = 𝐻2 𝑂 𝑀𝐵 = 18 Φ𝐵 = 2.6 𝑉𝐴 = 48.4
𝑚𝑜𝑙

𝑇 − 𝜇𝐻2𝑂
273 − 1794𝑋10−6 𝑃𝑎 − 𝑠
289 − 1153.2𝑋10−6 𝑃𝑎 − 𝑠
293 − 993𝑋10−6 𝑃𝑎 − 𝑠
𝜇𝐵 = 1.1532 𝑐𝑝
1
𝐷𝐴𝐵 (1.1532 𝑐𝑝) 7.4𝑋10−8 ((2.6)(18))2
=
289 (48.4)6
𝑐𝑚2
𝐷𝐴𝐵 = 1.2372𝑋10−5
𝑠
b. Hayduk-Laudie

𝐷𝐴𝐵 = 13.26𝑋10−5 𝜇𝐵 −1.14 𝑉𝐴 −0.589

𝜇𝐵 =?
273 − 1794𝑋10−6 𝑃𝑎 − 𝑠
293 − 993𝑋10−6 𝑃𝑎 − 𝑠
𝜇𝐵 = 1153.2𝑋10−6 𝑃𝑎 − 𝑠
𝜇𝐵 = 1.1532 𝑐𝑝
𝑐𝑚3
𝑉𝐴 = [𝐶𝑙2] = 48.4
𝑚𝑜𝑙
𝐷𝐴𝐵 = 13.26𝑋10−5 (1.1532)−1.14 (48.4)−0.589
𝑐𝑚2
𝐷𝐴𝐵 = 1.1470𝑋10−5
𝑠

[85]
Ejercicio 4.7

El monóxido de carbono se difunde a través de una capa de 0.1 in de aire a un vaso

de ácido sulfúrico, donde es instantáneamente absorbido. La concentración de
monóxido de carbono en el borde exterior de la capa de aire es de tres moles por
ciento. Use la ecuación diferencial general de transferencia de masa para escribir la
ecuación diferencial correspondiente a este proceso de transferencia unidimensional
de masa en estado permanente. Haga una lista de las condiciones de frontera.

H2SO4

Condiciones de frontera:

 Unidireccional
 Estado estacionario
 Isotérmico
 DAB = constante
 No hay reacción química
𝜕𝑐𝐴 𝜕𝑁𝐴𝑥 𝜕𝑁𝐴𝑦 𝜕𝑁𝐴𝑧
[ + + ] = 𝑅𝐴
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝑐𝐴 𝜕𝑁𝐴𝑥 𝜕𝑁𝐴𝑦 𝜕𝑁𝐴𝑧
[ + + ] = 𝑅𝐴
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝑑𝑁𝐴𝑧
⟹ 1. =0
𝑑𝑧
𝑑𝑦𝐴
𝑁𝐴 = −𝑐𝐷𝐴𝐵 + 𝑦𝐴 (𝑁𝐴 + 𝑁𝐵 )
𝑑𝑧
𝐴 = 𝐶𝑂2 𝐵 = 𝐻2𝑆𝑂4
𝑑𝑦𝐴
𝑁𝐴 (1 − 𝑦𝐴 ) = −𝑐𝐷𝐴𝐵
𝑑𝑧
𝑐𝐷𝐴𝐵 𝑑𝑦𝐴
⟹ 2. 𝑁𝐴 = −
(1 − 𝑦𝐴 ) 𝑑𝑧
𝑑 𝑐𝐷𝐴𝐵 𝑑𝑦𝐴
2 𝑒𝑛 1 [− ]=0 𝑐 = 𝑐𝑡𝑒
𝑑𝑧 (1 − 𝑦𝐴 ) 𝑑𝑧

[86]
 𝑑 1 𝑑𝑦𝐴
[− ]=0
𝑑𝑧 (1 − 𝑦𝐴 ) 𝑑𝑧
𝐶. 𝐹. 1 𝑧 = 0 𝑦 = 0.03
𝐶. 𝐹. 2 𝑧 = 0.1′′ 𝑦 = 0
Ejercicio 4.8

El siguiente dibujo explica la difusión en la fase gaseosa en la vecindad de una

superficie catalizadora.

La componente A, se difunde a través de una película estancada a la superficie

catalizadora, en la cual se convierte instantáneamente en la especie B debido a la
reacción:
AB

Cuando la especie B se difunde en la película estancada, comienza a

descomponerse debido a una reacción de primer orden:
BA

La rapidez de formación de la componente A es igual a

Donde yB es la concentración de B expresada en forma de fracción molar.

Use la ecuación diferencial general de transferencia de masa para escribir la

ecuación diferencial que describa. Este proceso de difusión. Haga una lista de las
condiciones de frontera que podrían usarse en la solución de la ecuación
diferencial.

Condiciones:

 Estado estacionario
 Unidireccional
 Isotérmico
 DAB = constante
𝑑𝑁𝐴𝑧
= 𝑅𝐴
𝑑𝑧

[87]
 𝑑𝑁𝐴𝑧
= 𝑘1 𝑦𝐵
𝑑𝑧
𝑦𝐴 + 𝑦𝐵 = 1
𝑦𝐵 = 1 − 𝑦𝐴
𝑑𝑁𝐴
⟹ 1. = 𝑘1 (1 − 𝑦𝐴 )
𝑑𝑧
𝑑𝑦𝐴
𝑁𝐴 = −𝑐𝐷𝐴𝐵 + 𝑦𝐴 (𝑁𝐴 + 𝑁𝐵 )
𝑑𝑧
𝐴→𝐵
𝑁𝐴 1
= − = −1
𝑁𝐵 1
𝑁𝐴 = −𝑁𝐵
𝑑𝑦𝐴
⟹ 2. 𝑁𝐴 = −𝑐𝐷𝐴𝐵
𝑑𝑧
𝑑 𝑑𝑦𝐴
2 𝑒𝑛 1 [−𝑐𝐷𝐴𝐵 ] = 𝑘1 (1 − 𝑦𝐴 )
𝑑𝑧 𝑑𝑧
𝑑𝑦𝐴
−𝑐𝐷𝐴𝐵 = 𝑘1 (1 − 𝑦𝐴 )
𝑑𝑧 2
𝐶. 𝐹. 1 𝑧 = 0 𝑦𝐴 = 𝑦𝐴0
𝑧 = 𝛿 𝑦𝐴 = 0

[88]
Ejercicio 4.9

La humedad en el aire caliente y húmedo que rodea a una tubería de agua fría se
difunde continuamente hasta la superficie fría donde se condensa.

Reducir la ecuación diferencial general para la transferencia de masa y escribir la

ecuación diferencial especifica que describa este proceso de transferencia de masa
en estado estacionario.

Proporcionar dos condiciones de frontera que podrían utilizarse para resolver la

ecuación resultante diferencial.

Aire + agua  caliente
H2O

r R
z

Condiciones:

 Unidireccional
 Isotérmico
 DAB = constante
 No hay reacción química
1𝜕
(𝑟𝑁𝐴,𝑟 ) = 0
𝑟 𝜕𝑟
𝑑
⟹ 1. (𝑟𝑁𝐴,𝑟 ) = 0
𝑑𝑟
𝑑𝑦𝐴
𝑁𝐴 = −𝑐𝐷𝐴𝐵 + 𝑦𝐴 (𝑁𝐴 + 𝑁𝐵 )
𝑑𝑟
𝐴 = 𝐻2𝑂 𝐵 = 𝑎𝑖𝑟𝑒 𝑁𝐵 = 0
𝑑𝑦𝐴
𝑁𝐴 = −𝑐𝐷𝐴𝐵 + 𝑦𝐴 (𝑁𝐴 )
𝑑𝑟
𝑑𝑦𝐴
𝑁𝐴 (1 − 𝑦𝐴 ) = −𝑐𝐷𝐴𝐵
𝑑𝑟
𝑐𝐷𝐴𝐵 𝑑𝑦𝐴
⟹ 2. 𝑁𝐴 = −
1 − 𝑦𝐴 𝑑𝑟
𝑑 𝑐𝐷𝐴𝐵 𝑑𝑦𝐴
2 𝑒𝑛 1 (𝑟(− )) = 0
𝑑𝑟 1 − 𝑦𝐴 𝑑𝑟
[89]
 𝑑 𝑐𝐷𝐴𝐵 𝑑𝑦𝐴
(𝑟( )) = 0
𝑑𝑟 1 + 𝑦𝐴 𝑑𝑟
𝐶. 𝐹. 1 𝑟 = 𝑅 𝑦𝐴 = 1
𝐶. 𝐹. 2 𝑟 = ∞ 𝑦𝐴 = 0
̅̅̅
𝑃𝐴
𝑟 = ∞ 𝑦𝐴 =
𝑃
Ejercicio 4.10

En una cámara caliente de combustión se difunde oxígeno a través de aire hasta una
superficie de carbono donde reacciona para formar CO y CO2. La concentración de
oxígeno en z = 6 es de 2 1 moles por ciento. La reacción en la superficie puede
suponerse instantánea. No ocurre ninguna reacción en la película de gas. Determine
la rapidez de difusión del oxígeno por hora, a través de un área de un metro cuadrado
si:

𝑁2 = 79%
𝑦𝐴 = 0.21 { }
𝑂2 = 21%
𝑑𝑦𝐴
𝑁𝐴 = −𝑐𝐷𝐴𝐵 + 𝑦𝐴 (𝑁𝐴 + 𝑁𝐵 )
𝑑𝑧
a) se produce solamente bióxido de carbono en la superficie del carbono
2𝐶 + 𝑂2 → 2𝐶𝑂
𝑁𝐶𝑂 2
= − = −1
𝑁𝑂2 1
𝐴 = 𝑂2 𝐵 = 𝐶𝑂
𝑑𝑦𝐴
⃗⃗⃗⃗
𝑁𝐴 = −𝑐𝐷𝐴𝐵 + 𝑦𝐴 (𝑁𝐴 + 2𝑁𝐴 )
𝑑𝑧
𝑑𝑦𝐴
𝑁𝐴 (1 − 𝑦𝐴 ) = −𝑐𝐷𝐴𝐵
𝑑𝑧
𝛿 0.21
𝑑𝑦𝐴
𝑁𝐴 ∫ 𝑑𝑧 = −𝑐𝐷𝐴𝐵 ∫
0 0 1 − 𝑦𝐴
𝑐𝐷𝐴𝐵
𝑁𝐴 = − 𝐿𝑛(1.21)
𝛿

[90]
b) si se produce solamente monóxido de carbono en la superficie del carbono y
𝑑𝑦𝐴
⃗⃗⃗⃗
𝑁𝐴 = −𝑐𝐷𝐴𝐵 ⃗⃗⃗⃗𝐴 + ⃗⃗⃗⃗⃗
+ 𝑦𝐴 (𝑁 𝑁𝐵 )
𝑑𝑧
𝐶 + 𝑂2 → 𝐶𝑂2
𝑁𝐶𝑂2 1
= − = −1
𝑁𝑂2 1
𝐴 = 𝑂2 𝐵 = 𝐶𝑂2
𝑑𝑦𝐴
⃗⃗⃗⃗
𝑁𝐴 = −𝑐𝐷𝐴𝐵 + 𝑦𝐴 (𝑁𝐴 − 𝑁𝐴 )
𝑑𝑧
𝑑𝑦𝐴
𝑁𝐴 = −𝑐𝐷𝐴𝐵
𝑑𝑧
𝛿 0.21
𝑁𝐴 ∫ 𝑑𝑧 = −𝑐𝐷𝐴𝐵 ∫ 𝑑𝑦𝐴
0 0

𝑐𝐷𝐴𝐵
𝑁𝐴 = − (0.21)
𝛿
c) ocurre la siguiente reacción en la superficie del carbono;
3 C+2 O2  2 CO + CO2
𝑑𝑦𝐴
⃗⃗⃗⃗
𝑁𝐴 = −𝑐𝐷𝐴𝐵 + 𝑦𝐴 (𝑁𝐴 + 𝑁𝐵 + 𝑁𝐶 )
𝑑𝑧
4𝐶 + 3𝑂2 → 2𝐶𝑂 + 2𝐶𝑂2
𝐴 = 𝑂2 𝐵 = 𝐶𝑂2 𝐶 = 𝐶𝑂
𝑁𝐶𝑂2 2 2
=− 𝑁𝐵 = − 𝑁𝐴
𝑁𝑂2 3 3
𝑁𝐶𝑂 2 2
=− 𝑁𝐶 = − 𝑁𝐴
𝑁𝑂2 3 3
𝑑𝑦𝐴 1
⃗⃗⃗⃗
𝑁𝐴 = −𝑐𝐷𝐴𝐵 + 𝑦𝐴 (− 𝑁𝐴 )
𝑑𝑧 3
1 𝑑𝑦𝐴
𝑁𝐴 (1 + 𝑁𝐴 ) = −𝑐𝐷𝐴𝐵
3 𝑑𝑧
𝛿 0.21
𝑑𝑦𝐴
𝑁𝐴 ∫ 𝑑𝑧 = −𝑐𝐷𝐴𝐵 ∫
1
0 0 1 + 3 𝑁𝐴

𝑐𝐷𝐴𝐵 1
𝑁𝐴 = − 𝐿𝑛(1 + (0.21))
𝛿 3

[91]
Ejercicio 4.11

Calcule la rapidez de difusión del azúcar a través de una película de .S café de 0.1
cm de espesor cuando las concentraciones son de 14% y 6% respectivamente, en
cualquiera de los lados de la película. Suponga que la difusividad del azúcar en el
café en las condiciones especificadas es de 0.7X10-5 cm2 /seg. y la densidad de una
solución al 10% es de 1.013 g/cm3.

14%
0.1
6%

𝑐𝑚2
𝐷𝐴𝐵 = 0.7𝑋10−5
𝑠
𝑔𝑟
𝜌10% = 1.013 3
𝑐𝑚
𝑑𝑤𝐴
𝑛𝐴 = −𝜌𝐷𝐴𝐵 + 𝑤𝐴 (𝑛𝐴 + 𝑛𝐵 )
𝑑𝑧
𝑛𝐴 = −𝑛𝐵
𝑑𝑤𝐴
𝑛𝐴 = −𝜌𝐷𝐴𝐵
𝑑𝑧
0.1 0.6
𝑛𝐴 ∫ 𝑑𝑧 = −𝜌𝐷𝐴𝐵 ∫ 𝑑𝑤𝐴
0 0

𝜌𝐷𝐴𝐵
𝑛𝐴 = − (0.06 − 0.14)
0.1
𝑔𝑟 −5
𝑐𝑚2
𝑛𝐴 = −1.013 (0.7𝑋10 )(0.06 − 0.14)
𝑐𝑚2 𝑠
𝑔𝑟
𝑛𝐴 = 5.6728 𝑋10−7
𝑐𝑚2 − 𝑠

[92]
Ejercicio 4.12

Una bolita de naftalina de forma esférica y de 3/4 in de diámetro está suspendida en

aire estacionario. La presión y temperatura del sistema son de 1 atm y 165" F,
respectivamente. Determine la cantidad de naftalina que entra en la fase gaseosa
en un período de 24 h. La naftalina tiene un peso molecular de 128 g/mol y una presión
de vapor de 5 mm de Hg a 165" F.

R
r

A= Naftalina B=Aire
D=3/4’’ P=1 atm T= 165°F t=24h Pa=5 mm Hg
Condiciones:

 Estado estacionario
 Isotérmico
 DAB = constante
 Una dirección
 No hay reacción química
1 𝜕 2
(𝑟 𝑁𝐴,𝑟 ) = 0
𝑟 2 𝜕𝑟
𝑑 2
⇒ 1. (𝑟 𝑁𝐴,𝑟 ) = 0
𝑑𝑟
𝑑𝑦𝐴
𝑁𝐴 = −𝑐𝐷𝐴𝐵 + 𝑦𝐴 (𝑁𝐴 − 𝑁𝐵 )
𝑑𝑟
𝑁𝐵 = 0
𝑑𝑦𝐴
𝑁𝐴 = −𝑐𝐷𝐴𝐵 + 𝑦𝐴 (𝑁𝐴 )
𝑑𝑟
𝑑𝑦𝐴
𝑁𝐴 (1 − 𝑦𝐴 ) = −𝑐𝐷𝐴𝐵
𝑑𝑟
𝑐𝐷𝐴𝐵 𝑑𝑦𝐴
⇒ 2. 𝑁𝐴 = − + 𝑦𝐴
1 − 𝑦𝐴 𝑑𝑟

[93]
 𝐷𝑒𝑠𝑝𝑒𝑗𝑒 𝑑𝑒 1 ∫ 𝑑(𝑟 2 𝑁𝐴 ) = 0

𝑟 2 𝑁𝐴 = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒
2 𝑒𝑛 1
𝑟2 𝑐𝐷𝐴𝐵 𝑑𝑦𝐴
𝑁𝐴 = −
𝑟2 1 − 𝑦𝐴 𝑑𝑟
∞ 0
2
𝑑𝑟 𝑑𝑦𝐴
𝑟 𝑁𝐴 ∫ 2 = −𝑐 𝐷𝐴𝐵 ∫
𝑅 𝑟 𝑦𝐴 1 − 𝑦𝐴

1 1
𝑟 2 𝑁𝐴 ( ) = 𝑐𝐷𝐴𝐵 𝐿𝑛 ( )
𝑅 1 − 𝑦𝐴
1
[𝑟 2 𝑁𝐴 = 𝑐𝐷𝐴𝐵 𝑅𝐿𝑛 ( )] 4𝜋
1 − 𝑦𝐴
1
4𝜋𝑟 2 𝑁𝐴 = 4𝜋𝑐𝐷𝐴𝐵 𝑅𝐿𝑛 ( )
1 − 𝑦𝐴
1
𝑤𝐴 = 4𝜋𝑐𝑅𝐷𝐴𝐵 𝐿𝑛 ( )
1 − 𝑦𝐴𝑜
1 𝑎𝑡𝑚
𝑐=
𝐹𝑡 3 − 𝑎𝑡𝑚
(625 °𝑅)(0.7302 )
𝐿𝑏 − 𝑚𝑜𝑙 − °𝑅
𝐿𝑏 − 𝑚𝑜𝑙
𝑐 = 2.19𝑋10−3
𝐹𝑡 3
3
𝐷 4′′ 𝐹𝑡
𝑅= = = 0.0312 𝐹𝑡
2 2 12′′
2
𝑇1 = 298 °𝐾 𝑐𝑚
{ } 𝐷𝐴𝐵 = 0.0611
𝑃1 = 1 𝑎𝑡𝑚 𝑠
3
𝑃1 𝑇2 2 Ω| 𝑇1
𝐷𝐴𝐵𝑇2,𝑃2 = 𝐷𝐴𝐵𝑇1,𝑃1 ( ) ( )
𝑃2 𝑇1 Ω| 𝑇2
3
𝑐𝑚2 346.88 2 𝑐𝑚2 0.00107639 𝐹𝑡 2 3600 𝑠
𝐷𝐴𝐵 = 0.0611 ( ) = 0.0767 | || |
𝑠 298 𝑠 1 𝑐𝑚 1ℎ

𝐹𝑡 2
𝐷𝐴𝐵 = 0.2974
𝑠
̅̅̅̅
𝑃𝐴 5
𝑦𝐴𝑜 = = = 6.5789𝑋10−3
𝑃 760

[94]
 1
𝑤𝐴 = 4𝜋(2.19𝑋10−3 )(0.0312)(0.2974)𝐿𝑛 ( )
1 − 6.5789𝑋10−3
𝑙𝑏 − 𝑚𝑜𝑙
𝑤𝐴 = 1.6855𝑋10−6
ℎ
𝑙𝑏 − 𝑚𝑜𝑙 128 𝑙𝑏
𝑤𝐴 = 1.6855𝑋10−6 | |
ℎ 𝑙𝑏 − 𝑚𝑜𝑙
𝐿𝑏 454 𝑔𝑟
𝑤𝐴 = 2.1568𝑋10−4 | |
ℎ 1 𝑙𝑏
𝑔𝑟
𝑤𝐴 = 0.0979
ℎ
Ejercicio 4.13

Se utiliza una celda de Arnold en estado permanente en la determinación de la

difusividad del alcohol etílico en aire a 297 K y 1 atm. Si los resultados concuerdan con
el valor que aparece en el Apéndice J 1 y la celda tiene una sección transversal de
0.82 cm2 y una trayectoria de difusión de 15 cm, ¿qué cantidad de etanol deberá
suministrarse a la celda para mantener constante el nivel del líquido? A 297 K, la
presión de vapor de etanol es de 53 mm de Hg y su gravedad específica es de 0.79.

R-OH
15 cm

R-OH

P = 1 atm T= 297 °K Área = 0.82 cm2 DA = 53 mm Hg gc = 0.79

𝑑𝑦𝐴
𝑁𝐴 = −𝑐𝐷𝐴𝐵 + 𝑦𝐴 (𝑁𝐴 − 𝑁𝐵 )
𝑑𝑧
𝜕𝑐𝐴 𝜕𝑁𝐴 𝜕𝑁𝐴 𝜕𝑁𝐴
+[ + + ] = 𝑅𝐴
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝑁𝐴
=0
𝜕𝑧
𝑁𝐵 = 0
𝑑𝑦𝐴
𝑁𝐴 = −𝑐𝐷𝐴𝐵 + 𝑦𝐴 (𝑁𝐴 )
𝑑𝑧

[95]
 𝑐𝐷𝐴𝐵 𝑑𝑦𝐴
⟹ 1. 𝑁𝐴 = −
1 − 𝑦𝐴 𝑑𝑧

𝑦𝐴 + 𝑦𝐵 = 1
𝑦𝐴 = 1 − 𝑦𝐵
𝑐𝐷𝐴𝐵 𝑑𝑦𝐵
⟹ 𝑎𝑖𝑟𝑒 𝑁𝐴 = −
1 − 𝑦𝐴 𝑑𝑧
𝑑𝑦𝐵
𝑁𝐴 ∫ 𝑑𝑧 = −𝑐𝐷𝐴𝐵 ∫
1 − 𝑦𝐴
𝑦𝐵 = 1 − 𝑦𝐴
𝑧2 𝑦𝐵2
𝑑𝑦𝐵
𝑁𝐴 ∫ 𝑑𝑧 = 𝑐𝐷𝐴𝐵 ∫
𝑧1 𝑦𝐵1 𝑦𝐵
𝑦𝐵2
𝑁𝐴 (𝑧2 − 𝑧1) = 𝑐𝐷𝐴𝐵 𝐿𝑛
𝑦𝐵1
𝑃𝐵2
𝑦𝐵2 =
𝑃
𝑃𝐵1
𝑦𝐵1 =
𝑃
𝑃𝐵2
𝑐𝐷𝐴𝐵
𝑁𝐴 = 𝐿𝑛 𝑃
𝑧2 − 𝑧1 𝑃𝐵1
𝑃
𝑐𝐷𝐴𝐵 𝑃𝐵2
⟹ 𝑁𝐴 = 𝐿𝑛
𝑧2 − 𝑧1 𝑃𝐵1
𝑃
𝑐=
𝑅𝑇
1 𝑎𝑡𝑚
𝑐=
𝐿𝑡 − 𝑎𝑡𝑡𝑚
297 °𝐾 (0.0826 )
𝑚𝑜𝑙 − °𝐾
𝑚𝑜𝑙 1 𝑙𝑡
𝑐 = 0.0410 | |
𝑙𝑡 1000 𝑐𝑚3
𝑚𝑜𝑙
𝑐 = 4.1𝑋10−5
𝑐𝑚3
𝑐𝑚2
𝑇1 = 298°𝐾 ; 𝑃1 = 1 𝑎𝑡𝑚 {𝐷𝐴𝐵 = 0.132 }
𝑠

𝑐𝑚2
𝑇2 = 297°𝐾 ; 𝑃2 = 1 𝑎𝑡𝑚 {𝐷𝐴𝐵 = 0.132 }
𝑠
[96]
 𝑃𝐵2 = 760 𝑚𝑚 𝐻𝑔
𝑃𝐵1 + 𝑃𝐴 = 1 𝑎𝑡𝑚
𝑃𝐵1 = 760 − 53
𝑃𝐵1 = 707 𝑚𝑚 𝐻𝑔
𝑚𝑜𝑙 𝑐𝑚2
4.1𝑋10−5 3 (0.132 𝑠 ) 760
𝑁𝐴 = 𝑐𝑚 𝐿𝑛
15 𝑐𝑚 707
𝑚𝑜𝑙 3600𝑠
𝑁𝐴 = 2.6081𝑋10−8 | |
𝑐𝑚2 − 𝑠 ℎ
𝑚𝑜𝑙
𝑤𝐴 = 9.3893𝑋10−5 (0.82 𝑐𝑚2 )
𝑐𝑚2 − ℎ
𝑚𝑜𝑙 46 𝑔𝑟
𝑤𝐴 = 7.6992𝑋10−5 | |
ℎ 1 𝑚𝑜𝑙
𝑔𝑟
𝑤𝐴 = 3.5416𝑋10−3
ℎ
𝑤𝐴
𝑃𝐴 =
𝑉𝐴
𝑔𝑟
3.5416𝑋10−3
𝑉𝐴 = ℎ
𝑔𝑟
0.79 3
𝑐𝑚
𝑐𝑚3
𝑉𝐴 = 4.4831𝑋10−3
ℎ

[97]