Descargar como doc, pdf o txt
Descargar como doc, pdf o txt
Está en la página 1de 54

Problemario de Fenómenos de Transporte

INSTITUTO POLITECNICO NACIONAL

UNIDAD PROFESIONAL INTERDISCIPLINARIA DE


BIOTECNOLOGIA

PROBLEMARIO DE LA ASIGNATURA

FENÓMENOS DE TRANSPORTE

ELABORADO POR:

M. EN C. MARÍA GUADALUPE ORDORICA MORALES

2008

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

1. Estimación de viscosidad de un gas denso


Estimar la viscosidad del nitrógeno a 68 °F y 1000 psig
N2
Pc = 33.5 atm T=68F
Tc = 126.2 K P = 1000 psig
μc = 180 x10-6 g/cms 68 32
1000 psi 1 atm = 68 atm C 20 C
14.6061 psi 1.8
20 C 20 273.15 293.15 K
T T 293.15 K 2.3229
Tr c
126.2 K
Con los valores obtenidos de Tr y Pr, se obtiene el valor de μr (Gráfico Uyehara)
PP 68 atm
Pr c 2.029
r 1.12 33.5 atm
rc

r c

6 4
g g
1.12180 10 2.016 10
cm s cm s
2.016 x10 -4 g 1 lbm 1 kg 1 cm = 1.355 x10-5 lbm
cm s 0.453593 kg 1000 g 0.0328 ft ft s

2. Estimación de viscosidad de fluoruro de metilo (CH 3F) a 370 °C y 120


atm. CH3F
M = 34 g/mol ρc = 0.300 g/cm3
Pc = 58.0 atm T = 370 °C = 643.15 K
Tc = 4.55 °C =277.7 K P = 120 atm
T 643.15 K
c 7.70M 1 2 Pc 2 3Tc 1 6 Tr 2.3159
Tc 277.7 K
c 7.70 34 1 2 58 2 3 277.7 1 6 c
Pr 120 atm
263.38 10 6 poise P 2.206
Pc 58atm
Con los valores obtenidos de Tr y Pr, se obtiene el valor de μr (Gráfico Uyehara)
r 1.145
r
rc
c
6 4
1.145 263.38 10 g 3.015 10 g
cm s cm s
3.015 x10-4 g 1 lbm 1 kg 1 cm = 2.0269 x10-5 lbm
cm s 0.453593 kg 1000 g 0.0328 ft ft s

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

3. Viscosidad de gases de baja densidad


Predecir la viscosidad de oxígeno molecular, nitrógeno y metano a 20 °C y a presión atmosférica,
expresar los resultados en mPa·s.

T σ /K KT
Compuesto M
(K) (Å) (K) (y) (poise) (cpoise)
(x)
Oxígeno 32.00 293.15 3.433 113 2.5924 1.0818 2.0277 x10-4 0.0202
Nitrógeno 28.02 293.15 3.681 91.5 3.2038 1.0217 1.7475 x10-4 0.0174
Metano 16.04 293.15 3.822 137 2.1397 1.1488 1.0907 x10-4 0.0109

O2 N2 CH4
KT 293.15 2.5924 KT 293.15 3.2536 KT 293.15 2.1397
113 90.1 137
2
y 2y y y1 m x x1 y y y1
m 1
y 1

x 2x x 2x
x x1

1 1
MT
2.6693 10 5
2

O2
1.081 1.093
2.5934 2.50 1.093 1.0818
2.6 2.5
2.6693 10 5
32 293.15 2.0277 10 4
poise
2
3.433 1.0818
N2

1.014 1.022
3.2038 3.20 1.022 1.0217
3.3 3.2

2.6693 10 5
28.02 293.15 1.7475 10 4
poise
3.681 2 1.0217
CH4

1.138 1.156
2.1397 2.1 1.156 1.1488
2.2 2.1
16.04 293.15
2.6693 10 5 1.0907 10 4 poise
2
3.822 1.1488

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

Utilizando nomogramas para viscosidad de gases

O2 2000 10 7 poise 0.02 centipoise


7
N2 1760 10 poise 0.0176 centipoise CH4

1000 10 7 poise 0.01centipoise

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

z zmax

max

x 0 vz 0

4. Flujo de una película descendente


Deducir el perfil de velocidad y la velocidad media, situando en el origen de coordenadas de forma
que x se mida a partir de la pared (es decir x = 0 corresponde a la pared y x = σ a la superficie libre
de la película). Demostrar que la distribución de velocidad viene dada por
2 2
g x 1 x
vz cos
2
Demostrar como se puede llegar a la distribución de velocidad de la ecuación anterior a partir de ecuación:
2

v 2 1 x
g cos 0 x x
z
2 x x
L
z0 x
0 vz 0

Entradas Salidas
x x x
Transporte viscoso xz LW xz LW
z0 z L
Transporte cinético vz 2 W x vz 2
W x
Volumen W x L g cos

2 2
LW LW v W x v W x W x Lg cos0
xz x xz x x z z 0 z z L

1 xz LW
x
xz LW
x x
W x L g cos
xz x x x
LW xz LW W x Lg cos
LW x LW x
xz x x x
xz g cos d xz g cos
x dx
d xzg cosdx
xzg cos x C1
Condiciones de frontera : xz 0x
C1 g cos
xzg cos xg cos
Ecuación de Newton
dv
z
xz dx

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

igualando ecuaciones
dv z
dx g cos x g cos dvz g cos x dx

dv g cos xdx
z

x2
g cos x
2
vz C2
Condiciones de frontera : vz 0 x 0
x2
g cos x
2
C2 0 vz
2 2
g cos x2
v 2
x
2
z 2
g cos 2
x x2
vz 2
2
2 2
g x 1 x
v cos
z 2

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

5. Flujo laminar en una rendija estrecha


Un fluido viscoso circula con flujo laminar por una rendija formada por dos paredes planas
separadas una distancia 2B. Efectuar un balance diferencial de cantidad de movimiento y obtener
las expresiones para las distribuciones de densidad de flujo de cantidad de movimiento y de
velocidad.
P0 PL
xz x
L
v P0 PL B2 1 x 2

z
2 L B
en las que P pg h pg z

Entradas Salidas
x x x
Transporte viscoso xz LW xz LW
2 z 0 z L
Transporte cinético vz W x vz 2 W x
Presión P0 x W PL x W
Volumen W x Lg

xz x xz x x 2 0
LW LW v W x v2 W x p x W pL x W W x L g 0
z z
z 0 z L

xz LW x xz LW p x W pL xW W xL g
x x 0
W xL W xL
xz
xz p p
1 g
x x x 0 L

x L
xz x x xz x
g p0 p L
d xz p 0p L
g
x L dx L
p p p p
0 L L
g 0 L L

xz xz 1
d dx x g x C

Condiciones de frontera xz 0 x 0 C 0
1

xz P0 PL x
L

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

Ecuación de Newton
dvz
xz
dx
dvz P0 PL x
dx L
dvz P0 PL x dx
L
dv P P x dx vz P P x 2 C2
0 L 0 L

z L L 2
Condiciones de frontera : v 0 x 2B
z

P P 4 B2
C2 0 L

L 2
P P x2 P P 4 B2
0 L 0 L

vz L 2 L 2
PP z
v 0 L 4 B2 x2 2 L
2 22
B P0 PL B2
v
2 4 B x 2
B B
z

2 L
2 2
P0 PL x
vz B 4
2 L B

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

6. Flujo laminar en una película que desciende por el exterior de un tubo circular
En una experiencia de absorción de gases, un fluido viscoso asciende por el interior de un pequeño
tubo circular, para descender después por la parte exterior del mismo. Aplicar un balance de
cantidad de movimiento a una envoltura de película de espesor r, tal como se indica en la figura.
Obsérvese que las flechas de entrada de cantidad de movimiento se toman siempre en la dirección
r positiva al efectuar el balance, aun cuando en este caso ocurre que la cantidad de movimiento
fluye en la dirección r negativa. Demostrar que la distribución de velocidad en la película
descendiente (despreciando los efectos finales) es:
v g R2 1 r
2
4 R

Entradas Salidas
r r r
Transporte viscoso rz L2 r rz L2 r
2 z 0 2 z L
Transporte cinético vz 2 r r vz 2 r r
Fuerza de gravedad 2 r r Lg

L2 r L2 r v2 2r r
rz r rz r r z z0

rz L2 r r
rz L2 r r r
2 rr L g
2r L 2r L
rz r rz r rd
1 r g
r r r rz

r dr
vz 2 2 r r z L 2 r rL g 0

r g

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

r
r d rz g r dr r rz g 22
C1
rzg
r2 C1

2r r

Condiciones de frontera
2
rzg rg aR

22 r
Ecuación de Newton
dv
rz rz

dr g aR 2
dvrz gr
dr 2 2r
dvrz gr g aR 2
dr 2 2 r
2
dvg aR r dr
rz2r

g 2
r2
vrz aR ln r C2
2 2
g R2
2

C
Condiciones de frontera vz 0 r R 2 aR ln R
2 2
g r 2g
v aR 2 ln r aR 2 ln R R2
rz
22 2 2
g r2 R2
v aR 2 ln r aR 2 ln R
rz
222
2
g r
R2 r2
vrz aR ln
2 R 2 2
2
g R2 r r2
v
rz 2a ln 1
4 R R

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

7. Flujo en tubos concéntricos con movimiento axial del cilindro exterior


Considerar el sistema representado en la figura, en el que la varilla cilíndrica se mueve con
velocidad V. La varilla y el cilindro son coaxiales. Hallar la distribución de velocidad en estado
estacionario y la velocidad volumétrica de flujo. Este tipo de problemas se presentan en el
recubrimiento de alambres con barniz

vz vz v vz vz p 1 vz 1 2vz 2
vz
vr vr r 2 2 2 gz
t r r z z r
r r r z
vr 0 v 0 vz ?
dv
0 1 v z
1d z d d z
0
r r 0
v
r
r r r dr dr r
d z dr dr
v
dr 0 dr
dr
r dvz C
dr 1

dv C1 dr d C dr
r
v r

z z 1

vz C1 ln r C21 v V
z max

Condiciones d frontera CL1 r KR


vz 0
e r R
CL2
C1 ln KR C2 Vmax 2 C1 ln R C2 0 3
Se despeja C2 de 3
en 2

C2 C1 ln R se sustituye Vmax C1 ln KR
C1 ln R
V
KR
max C1 ln Vmax C1 ln K
R
C Vmax C Vmax ln R
2
1
ln ln K
K
V V V
max Vmax ln R max ln r ln R max r
vz ln r vz ln
ln ln ln ln R
K K K K
r
ln
vz R
V ln K
max
M. en C. María Guadalupe Ordorica Morales
Problemario de Fenómenos de Transporte

Para obtener rz

Ecuación de Newton
dvz
rz
dr
r
Vmax ln
r d R rz Vmax d r
ln
dr ln K ln K dr R

V d rz Vmax 1
z max ln r ln R ln K
dr ln K r

z
M. en C. María Guadalupe Ordorica Morales
Problemario de Fenómenos de Transporte

8. Una cañería de agua consiste en un conducto de presión hecho de concreto de 18 in de


diámetro. Calcule la caída de presión en un tramo de 1 milla de longitud debido a la fricción en la
pared del conducto, si éste transporta 15.0 ft 3/s de agua a 50 °F.
A 50 y 60 °F el peso específico de H2O es 62.4 lbf/ft3.

Se analiza la ecuación general de energía


P
1
z V
1
2
h h h P z
2
V2
2

A L 2
1
2g r
2g
P P
1 2 hL

P1 P 2 hL
Datos del problema
42 102 ft
concreto

H2O 62.4lb f ft 3
5280 ft
z 1 milla 5280 ft
1 milla
1 ft
D 18 in 1.5 ft
12 in
1 ft
r 9 in 0.75 ft
12 in
3
15 ft s
T 50 F
Cálculo de viscosidad, # de Reynolds y velocidad de flujo.
g 2.20 10 3 l m
30.48 cm l m

b b
1.27 10 2 8.53 10 4
cm s 1g 1 ft ft s
2
Ar 0.75 ft 2 1.7671 ft 2
V A 3
s ft
15 ft 8.4882
1.7671 ft 2 s
Re V DP

3
ft 1.5 ft 62.4 lb f
8.4882

Re s ft 930513 flujo turbulento


8.53 10 4 lbm
ft s

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

Del diagrama de Moody

f 0.06 D 1.5 ft 3.57


42 102 f
t

h f L
L V2
D 2g
h 0.06 5280 ft 8.4882 ft s 2 236.2866 ft
L
1.5 ft 2 32.2 ft s2
P hL 2

lb l 1 ft
b
62.4
f
3
14744.2891 f
2

12i
P 236.2866 ft

ft ft
n
lbf
P 105.34
2
in
M. en C. María Guadalupe Ordorica Morales
Problemario de Fenómenos de Transporte

9. En la figura se observa una parte del sistema de protección contra incendios en el cual una
bomba saca agua a 60 °F de un recipiente y la transporta al punto B, con una rapidez de flujo de
1500 gal/min. Calcule la altura h, requerida para el nivel del agua en el tanque, con el fin de
mantener 5.0 lb/pulg2 relativa presión en el punto A.
A 50 y 60 °F el peso específico
del agua es de 63.4 lbf/ft3
Para tubo de acero calibre 40 de
8 in el diámetro interno es
0.6651 ft
Para tubo de acero calibre 40 de
10 in el diámetro interno es
0.8350 ft.

Pz V2 h hh P z V2
1 1 2 2

1 A r L 2
V 1500 gal min 2g 2g
P
P2 5.0lb in2 z 2

1
D1 0.8350 ft
3
D2 0..6651 ft 62.4 lbf 1 ft lbf
0.03611

62.4lbf ft
3
ft3 12in in3
z2 25 ft 5.0 lb
in2 1 ft
z1 ? z1 138.4657 in 11.53 ft

0.03611 lbf 12in


in3

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

10. Predicción de conductividades caloríficas de gases a baja densidad.


a) Calcular la conductividad calorífica del argón a 100 °C y 1 atm de presión, utilizando la
teoría de Chapman-Enskog y las constantes de Lennard-Jones deducidas de los datos
de viscosidad.

M k kT k
(Å) (°K)
Ar 39.944 3.418 124 3.00 1.039

T 100 C 373.15 K
k 1.9891 10 4 TM
2
k

k 1.9891 10 4
373.15
39.
944 3.418 2

1.039
5
k 5.008 10 cal
cm s K

b) Calcular las conductividades caloríficas de óxido nítrico (NO) y del metano (CH 4) a 300 °K
y presión atmosférica, utilizando los siguientes datos para las mismas condiciones

M 10 7 Cp
(g/mol) (g/cm s) (g/mol °K)
NO 30.01 1929 7.15
CH4 16.04 1116 8.55
5 k Cp
R

4 M
cal R
1.987 mol K
5 1929 10 7
kNO 7.15 1.987 5 1116 10 7
4 30.01 k
CH4 8.55 1.987
4 16.04
kNO 619.34 10
7
cal
cm s K k 767.68 10 7 cal
CH4
cm s K
M. en C. María Guadalupe Ordorica Morales
Problemario de Fenómenos de Transporte

11. Predicción de la conductividad calorífica de un gas denso.


Predecir la conductividad calorífica del metano (CH 4) a 110.4 atm y 52.8 °C por los dos métodos
siguientes:
a) Utilizando el diagrama de Owens, tomando las propiedades críticas que sean necesarias
T 52.8 C 325.95 K
P 110.4 atm
Tc Pc kc x10-6
(°K) (atm) (cal/s cm °K)
CH4 190.7 45.8 158

T P
Tr Pr
Tc Pc
325.95
110.4
Tr 190.7 1.709
P 2.41
r 45.8
Diagrama de Owens

k 0.77
r

kr k
kc
k kr kc
6
cal
k 0.77 158 10
s cm K
cal 3600 s 100 cm 1 Kcal
k 1.2166 10 4
s cm K 1h 1m 1000 cal
k 0.04379 Kcal
hm K

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

12. Conducción de calor desde una esfera a un fluido estancado


Una esfera caliente de radio R esta suspendida en una gran masa de fluido en reposo. Se desea
estudiar la conducción de calor en el fluido que rodea la esfera. Se supone que los efectos de la
convección libre son despreciables.
a) Plantear la ecuación diferencial que describe la temperatura T del fluido circundante en
función de r, la distancia desde el centro de la esfera. La conductividad calorífica k del
fluido es constante.
b) Integrar la ecuación diferencial y utilizar las siguientes condiciones límite, para determinar
las constantes de integración.
CL1: para r R T TR
CL2: para r T T

Entradas Salidas
2 r 2 r r
4 r qr 4 r qr
r
r r
4 r 2 qr 4 r 2 qr 0

4 r 2 qr r
4 r 2 qr r r 0

4r
2 r 2
r q r q
r r r r 0 1
r

r 2 qr r2 q r
0 d r 2 qr 0
r r r

r dr
Se sustitye qr por la ley de Fourier
qr k dT
dr
2
d k dT 0
r
dr dr
2
k dT
dr 0 dr
dr
r 2 k dT C
dr 1

dT C1 d
kr2 r
dT C1 r 2 dr
k
T C
1 C2
kr

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

Se aplican las condiciones iniciales


CL1: para r R T TR
CL2: para r T T
TR C1 C2

k R
T C
1 C T C
2 2
k
T C1 T
R
k R
C1 T TR k R
TT k R
T R T k
r
T T T T R
R
r
Se sustituye en la ley de Fourier para obtener una ecuación para qr
q dT
k dr
r
d R
k T T T R
q dr r
r
kT T R
R
q
r2
r

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

13. Calentamiento viscoso en el flujo a través de una rendija


Deducir una expresión para la distribución de temperatura T(x) en un fluido viscoso que circule
con flujo laminar por el espacio comprendido entre dos grandes láminas paralelas tal como se
indica en la figura. Ambas láminas se mantienen a temperatura constante T0. Téngase en
cuenta el calor generado por disipación viscosa. Desprecie la variación de k y μ con la
temperatura.

De las ecuaciones de variación (coordenadas rectangulares)


T
2 T
2 2 T vx 2 vy 2 vz 2
T T T T
pC v x v y v z k
v

t x y x2 y2 z 2 2 xy z
z
2 2 2
vx vy vx vz vy vz
y x z x z y
quedando
2 z 2
T v
k
2 0
x x 2
v x V T Vb x C1
vz Vb
z
b b x b x k b
2 2
T V V 2

k b
T
b

x C1 x
x2 b k
2 b
T Vb 2 2
T Vb C 1 x C2
x x k b x
2
2k b
T Vb
x
x 0 T T0
x k b
CL1
Se aplican condiciones límite CL2
x b T T1

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

CL1
To C 2
CL2
2
T V b2 Cb T
b

1
2k b 1
2
o

Vb
T1 To
C1 2k
b
V 2 T T V 2

C b bo
C b

1 1
2k b b 2k b

Se sustituyen valores de C1 y C2

2
V V2
T 2 b
x 2 b
x T
o
b 2k b
1 2
x x
k T To V
b 1
2k
b b

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

14. Temperatura máxima de un lubricante


Un aceite actúa como lubricante de dos superficies cilíndricas como las de la figura. La velocidad angular del
cilindro exterior es de 7908 rpm. El radio del cilindro exterior es de 5.06 cm y la luz entre los 2 cilindros,
0.027 cm ¿Cuál es la máxima temperatura en el aceite si se sabe que la temperatura de ambas paredes es de
70 °C? Las propiedades físicas del aceite son:
Viscosidad 92.3 cp
Densidad 1.22 g/cm3
Conductividad calorífica 0.0055 cal/s cm °C

Del problema anterior, aplicando las ecuaciones de variación se obtiene

2 2 Se aplican condiciones iniciales


T vz
k 0 CL1 x o T To
2
x CL2 x b T Tb
x
V T o C2
v xV vz b
2
z
B b x b T V b2 Cb T
b
2 2
T V b
2k b 1 o
k b
2
Vb
x2 b Tb To
T Vb 2
C1 2k
k b b 2
x x T T Vb
2 V2
T Vb C b b o
C 2k b
1 1
x k b
x
2k b b
2 Se sustituyen valores de C1 y C2
T V x C 2
b V V 2
2
T b
x b
x T
x k b 1
o
2
2k b 2k b
V
b
T To 1 2
x x
T b x C1 x Vb 1
k 2k
2
b b
V
T x2 C x C
b

2
2k b 1

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

15. Transmisión de calor en el acoplamiento de una barra de combustible nuclear.


Considere una barra larga de combustible nuclear que esta rodeada por una plancha anular de un
revestimiento de aluminio, tal como se indica en la figura. Debido al proceso de fisión, se produce
calor en el interior de la barra de combustible; el desarrollo de calor depende de la posición,
variando la intensidad del manantial calorífico de acuerdo con la expresión aproximada:
r 2
S S 1 b
n n0 R
F
Siendo Sn0 el calor producido por unidad de volumen y unidad de tiempo para r=0, y r la distancia
desde el eje de la barra de combustible, si la superficie externa de la vaina de aluminio esta en
contacto con un liquido regrigerante cuya temperatura es Tr y el coeficiente de transmisión de
calor en la interfase vaina-refrigerante es hL. Las conductividades caloríficas de la barra y la vaina
son kF y kC.

Balance 1 (Barra de combustible)


Entradas Salidas
r r r
2 r L qr 2 r L qr
Volumen 2 r L r Sn
2 r L qr r 2 r L qr r r 2 r L r Sn 0

r r r
2 r L qr 2 r L qr 2 r L r Sn
2 L r 2L r
rq rq
r r r r r
r Sn 1
r

r 2 qr r2 q r
r Sn rd r
r Sn
r r r
q
r dr
r dqr Sn r dr
r q Sn r2 C q Sn r C
1
r r
1
2 2 r

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

CL1 r 0 qr 0
Se aplican condiciones límite CL2 r Rf T TF
C1 0
qr Sn r
2
qr . k dT
dr
dT r
k dr Sn 2
Sn
k dT 2 r dr T
Sn 2
4 k r C2
Sn 2
TF4 k R F C2
Sn
C2 TF 4 k RF 2
Sn Sn
T 4 k r 2 TF 4 k R F 2

Sn 2 2
T TF 4 k RF r

Balance 2 (Revestimiento de Aluminio)

Entradas Salidas
r r r
2 r L qr 2 r L qr
2 r L qr r 2 r L qr 0
r r

r r
2 r L qr 2 r L qr
r 0
2 L r
rq r q
r r r r r
0 1
r

r 2 qr r2 q r
0 r d r 0
r r r
q
r dr
r dqr 0 dr

r qr C3 qr C3
r

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

qr k
dT
dr

dT C 3
k dr r
dr
k dT C3 r T
C 3
k ln r C4
CL3 r R f T TF
Se aplican condiciones límite CL4 r Rc T Tc
T C3 ln R C
f 4
F
k
C
T 3 ln R C
4
c
k c

Despenjando C4 de ambas ecuaciones


C4 TF C3 ln R f
k
C
C T 3 ln R
4 c
k c

T C3 ln R T C3 ln R
f
F
k c
k c
C C C
TF Tc 3 ln Rc 3 ln R f 3 ln Rc ln R f
k k k
C3 k TF Tc
R
ln c R
f

k TF Tc
R
ln c
R
C4 TF f
ln R f
k
C4 TF TF Tc ln R f
R
ln c
Rf Balance 3 (Refrigerante)
T TF Tc ln r T TF Tc ln
R f
R F R
ln c ln c
Rf R
f

T T Rf
T T F c
ln
F
R r
ln c
R
f

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

16. Conducción de calor un anillo circular


El calor fluye a través de una pared anular cuyo radio interno es r0 y el externo r1. La conductividad
calorífica varía linealmente con la temperatura desde k0 a la temperatura T0 hasta k1 a la temperatura
T1. Deducir una expresión para el flujo de calor a través de la pared situada en r = r0.

Entradas Salidas
r r r
2 r L qr 2 r L qr
2 r L qr r 2 r L qr 0
r r

r r
2 r L qr r 2 r L qr 0
2 L r
rq r q
r r r r r
0 1
r

r 2 qr r2 q r
0 r d r
0
r r r
q
r dr
r dqr 0 dr
r
rq C qr C1
1
r
qr k dT

dr
k dT C
1

dr r
k dT C1 dr
r
C
T 1 ln r C2
k
CL1 r ro T To k ko
Se aplican condiciones límite CL2 r r1 T T1 k k1

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

C1
T
k
ln r C
o o 2
o
C1
T
k
ln r C
1 1 2
1

Despenjando C2 de la segunda ecuación


C T C3 ln r
2 k1 1 1

Se sustituye en la primer ecuación


T C1 ln r T C3 ln r
o k o 1 k 1
o 1

T T C1 ln r ln r
o 1
k1 1 o

C k 1 To T1
1 r
ln 1r
0

k To T1
r
ln 1r
C 2 T
0
ln r
1
k1 1

C T To T1 ln r
2 1 r 1
ln 1r
o

T k 1 T0 T1 ln r T To T1 ln r
r 1 r 1
k ln 1r
o
ln 1r
o

To T1 k
1
T ln r1 ln r T1
r1 r k
ln
0
q k dT
r
dr
kT T T T
d
q k 1 0 1
ln r T o 1
ln r1
r
dr r 1 r
k ln 1r ln 1r
o o

To T1 k
qr
r1 r r
ln
o

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

QqA
A 2 rL
Q To T1 k2 r L
r1 r r
ln
o
1
Q k To T1 2 L ln r1 ln ro
kk
k 1 2 2
2 1
k k
Q 1
T T2 L ln r ln r
2 o 1 1 o

r r
1 o 0

ro
r
1 1
ro
r1 ro

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

17. Estimar DAB para el sistema argón-oxígeno a 293.2 K y 1 atm de presión total. Utilizar

Tc Pc M /K
(K) (atm) (Å) (K)
A- Argón 151 48 39.94 3.418 124
B - Oxigeno 154.4 49.7 32 3.433 113
1.823
a) La ecuación de Slattery

P DAB T
2.745 10 4
1 T T
PcA PcB 1 TcA TcB 5 1 1 2 cA cB
12

M
3

MA B
1.823 1 5
T 1 1 12
2.745 10 4 PcA PcB 3
TcA TcB 12
T T MA M B
cA cB
DAB
P
4 1.823 1 5
293.2 1 1
1
2
D
AB 2.745 10 48 49.7 3
151 154.4 12

1 39.94 32
5
1
1
5
4
.
4
D cm2
AB 0.1888
s
b) La ecuación de Chapman - Enskog
3
1 1
T 3.418 3.433
M AB
3.4255
D 0.0018583 MA B
2
AB 2 AB
P AB AB
Donde : 124 114 118.3722
K
AB
AB

2 KT 293.2 2.4769
AB
A B
AB
118.3722
K K K
Se interpola con los valores siguientes

KT AB AB
0.9996 1.012
AB
2.4769 2.40 0.9996 1.002464
2.40 1.012 2.50 2.40
2.50 0
.
9
9
9
6

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

3
1 1
293.2
DAB 0.0018583 39.94 32
2
1 3.4255 1.0024644
D cm2
AB 0.1881
s

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

18. Estímese DAB para una mezcla constituida por 80 moles por ciento de metano y 20 moles por
ciento de etano a 136 atm y 313 K. El valor experimental de (PD AB)° a 293 K es 0.163 atm cm2/s.

Tc Pc xi
(K) (atm)
A-Metano 190.5 45.8 0.8
B -Etano 305.4 48.2 0.2

Datos
T 313 K
P 136 atm
P DAB 293 K 0.163 atm cm2
s
b
T
P DAB T P DAB 293K
293
2
313 0.1823
0.1838 atm cm
P DAB 313 K 0.163 s
293
Pc 'xi Pci 46.28 atm
Tc 'xi Tci 213.4 K
P' P 136 2.938

r P ' 46.28
c

T' T 313 1.4667

r
Tc ' 213.4
M. en C. María Guadalupe Ordorica Morales
Problemario de Fenómenos de Transporte

PD AB
0.73
P DAB
atm cm2
0.1838 2
0.73

s cm
DAB 9.8657 10 4

136 atm s

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

19. Se conoce el valor de DAB (0.151 cm2/s) para el sistema CO2 –aire a 293 K y 1 atm. Extrapólese
DAB para 1500 K utilizando los métodos siguientes.
Tc Pc M /K
(K) (atm) (Å) (K)
A- CO2 304.2 72.9 14.01 3.996 190
B - Aire 132 36.4 28.91 3.617 97.0
1.823
a) Ecuación de Slattery.

P DAB T
2.745 10 4
1
PcA PcB 1 TcA TcB 5 1 1 2 T
cA
T
cB
3 12

MM
A B
1.823 1 5 1
T 1 1 2
3
2.745 10 4 PcA PcB TcA TcB 12
T T MA M B
D cA cB
AB
P
4 1.823 1 5 1
1500 1 1 2
D
AB 2.745 10 72.9 36.4 3 304.2 132 12
304.2 132 14.01 28.91
D cm 2
AB
2.2956
s
b) Ecuación de Chapman - Enskog

1 1
T3 3.996 3.617
M AB
3.8065
D 0.0018583 MA B
2
2 AB
AB
P AB AB
Donde : 190 97 135.7571
K
AB
AB

2 K T 293.2 11.0491
AB
A B
AB
118.3722
K K K
Se interpola con los valores siguientes

KT AB AB
0.6640 0.7424
AB
11.0491 10 0.7424 0.7341
10.0 0.7424 20 10
20.0 0
.
6
6
4
0
M. en C. María Guadalupe Ordorica Morales
Problemario de Fenómenos de Transporte

3
1 1
1500
D 0.0018583 44.01 28.91
AB
1 3.8065 2 0.7341
D cm 2
s
AB 2.4297

c) Gráfico de Slattery

Pc ' xi Pci 0.5 72.9 0.5 36.4 54.65 atm Tc ' xi Tci 0.5
304.2 0.5 132 218.1 K
P' P 1 0.1829
r
Pc ' 54.65
T' T 1500 6.8775
r
Tc ' 218.1

P DAB
No es posible obtener un valor para ya que la gráfica es para gases densos, y se está P
DAB
trabajando con un gas ideal, además de que el valor de T r obtenido no se encuentra en el gráfico.

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

20. Difusión de metano a través de helio


Un tubo contiene CH4 y He gaseoso a 101.32 kPa de presión y 298 K. En un punto, la presión
parcial del metano es PA1 = 60.79 kPa y en otro 0.02 m de distancia P A2 = 20.26 kPa. Si la presión
total es contante en todo el tubo. Calcule el flujo específico de CH 4 (metano) en estado
estacionario para contradifusión equimolar.

Z2= 0.02 m
PA1 = 20.26 kPa
Z1= 0
PA1 = 60.71 kPa

P1 = 101.32 kPa = 1 atm


T=298K
Tc Pc M /K
(K) (atm) (Å) (K)
A- CH4 190.7 45.8 16.04 3.822 137
B - He 5.26 2.26 4.003 2.576 10.2
1.823
P DAB T
2.745 10 4
1
PcA PcB 1 TcA TcB 5 1 1 2 T
cA
T
cB

M M
3 12

A B
1
T 1.823 1 5 1 1 2
4 3 12
2.745 10 PcA PcB TcA TcB
T T MA M
cA cB B
DAB
P
4 1.823 1 5 1
298 1 1 2
DAB 2.745 10 45.8 2.26 3 190.75.26 12
190.7 5.26 16.04 4.003
2
cm 1m 2 m 2
D AB 0.76334 7.6334 10 5
2
s 10000 cm s

D AB P P P
NA ln A2

RT z 2 z P P 1 A1

7.6334 10 5 101325 101325 20260 4


kg mol
N ln 1.0818 10
A
8314.3 298 0.02 101325 60790 m2 s
C
JA NA C
A
NA NB

J A NA x A NA
kg mol
1.0818 10 4 0.5 1.0818 10 4 5.409 10 5 2 A ms

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

21. Contradifusión equimolar de NH3 y N2 en estado estable.


A través de un tubo recto de vidrio de 2.0 ft (0.610 m) de longitud con diámetro interno de 0.080 ft
(24.4 mm). Se produce una contradifusión de amoniaco gaseoso (A) y nitrógeno gaseoso (B) a 298
K y 101.3 kPa. Ambos extremos del tubo están conectados a grandes cámaras de mezclado
colocadas a 101.32 kPa. La presión parcial de NH 3 en una cámara es constante e igual a 20.0 kPa y
en la otra cámara la presión es 6.666 kPa. La difusividad a 298 K y 101.32 kPa es 2.30 x10 -5 m2/s.

L
NB

x PA1 PA2

PB1 PB2
→ PA1 20.0 kPa
NA z
P 6.666 kPa
z1 z2 A2

L 0.61 m
a) Calcule la difusión del NH3 en lb mol/h y kg mol/s

D AB P PP
NA ln A2

RT z 2 z PP 1 A1

2.30 10 5 101325 101325 6666 7


kg mol
NA ln 2.3411 10
8314.3 298 0.61 101325 20000 m2 s
Ar 20.0122 2 4.6759 10 4 m2
7
kg mol 4 2
10 kg mol
NA A 2.3411 10 4.6759 10 m 1.0946 10

m2 s s
kg mol 2.2046 lb 3600 s lb mol
1.0946 10 10 8.688 10 7
s 1 1h h
kg
b) Calcule la difusión del N2
NA NB

c) Calcule las presiones parciales en un punto situado a 1.0 ft (0.305 m) en el tubo y


grafíquese PA, PB y P en función de la distancia z

6.666 20
P 0.305 20 13.333 kPa
0.61
P 13333 Pa 1.33 104 Pa

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

22. Difusión de A a través de B en reposo y efecto del tipo de límite sobre el flujo específico.
Se difunde amoniaco gaseoso a través de N 2 en estado estacionario, donde N 2 es el gas que no se
difunde, puesto que es insoluble en uno de los límite. La presión total es 1.013 x10 5 Pa y la
temperatura marca 298 K. La presión parcial de NH 3 e un punto es 1.333 x10 4 Pa y en el otro
punto, situado a una separación de 20 mm es 6.666 x10 3 Pa. El valor de DAB para la mezcla a 1.013
x105 Pa es 2.30 x10-5 m2/s.
a) Calcule el flujo específico de NH3 en kg mol/s m2
Datos :
m2
DAB 2.30 10 5
s
5
P 1.013 10 Pa
T 298K
PA1 1.333 104 Pa
PA2 6.666 103 Pa
D AB P P P
N
A ln A2

RT z 2z PP 1 A1

2.30 10 101325
5
ln 101325 1.333 10 4
3.4332 10 6 kg mol
NA 3
8314.3 298 0.02 101325 6.666 10 m s
2

b) Haga lo mismo que en (a) pero suponiendo que el N 2 tambien se difunde, esto es, ambos
límites son permeables a los dos gases y el flujo específico es una contradifusión
equimolar. ¿En qué caso es mayor al flujo específico?
d xA C A
N cDAB N A NB d z C
A
NB 0 NA DAB d PA
RT d z
N
A
z
2 D P
A2

NA dz
AB dP
A
z RT P
1 A1
D
NA z z P P
2
AB

1 RT A2 A1

NA
D P P
AB A2 A1

RT z2 z1
2.30 10 5 6.666 103 1.333 104
6
kg mol
NA 8314.3 298 3.0931 10
0.02 m2 s

La difusión es mayor en el caso planteado en el primer inciso.

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

23. Sublimación de pequeñas esferas de yodo en aire estático


Una esfera de yodo, 1 cm de diámetro, se encuantra en aire estático a 40 °C y a 747 mmHg de
presión. A esta temperatura la presión de vapor del yodo es de 1.03 mmHg. Se desea determinar
la difusividad del sistema yodo-aire midiendo el índice de sublimación. Estimar la difusividad para
el sistema aire-yodo a la temperatura y presión dadas anteriormente.

Datos :
T 40 C 313.15 K
P 747 mmHg 0.9828 atm
Tc M /K
(K) (Å) (K)
A-Aire 132 28.7 3.617 97.0
B - Yodo 800 253.82 4.982 550

1 1
T 3 3.617 4.982
M AB
4.2995
D 0.0018583 MA B
2
2 AB
AB
P AB AB
Donde : 97 550 230.9761
K
AB
AB

2 K T 313.15 1.3557
AB
A B
AB
230.9761
K K K
Se interpola con los siguientes valores

KT AB AB
1.233 1.253 1.3557 1.35 1.253 1.25072
AB
1.35 1.253 1.40 1.35
1.40 1.233

313.15 3 1 1
DAB 0.0018583 28.7 253.82
0.9828 4.2995 2 1.25072
cm2 s
DAB 0.08925

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

24. Deducción alternativa de la difusión a través de una película estancada.


En las ecuaciones (1) se obtuvo una expresión para calcular la velocidad de evaporación al
diferenciar el perfil de concentración encontrado anteriormente. Demostrar que los mismos
resultados pueden deducirse sin tener que encontrar el perfil de concentración. Nótese que en
estado estacionario NAz es una constante según la ecuación (2), luego la ecuación (3) puede
integrarse para obtener la ecuación (1).

N cDAB d xA cDAB d xB cD x B2

Az x
AB
ln --------- (1)
z z
1
1 xA1 dz z z B1 dz z z1 z 2 z x B1
1 1
d NAz 0 ------- (2)
dz
N cDAB d xA --------- (3)
Az
1 xA d z
Balance Entradas -
Salidas = 0
z z z
S N Az - S N Az = 0
1 C2 C1 z 1 cDAB ln 1 xA1
S N
S N Az z Az z z 0
S z C2 C1 z2 cDAB ln 1 xA2
lim N Az z z
N Az
z
0 d AZ 0 C1 z2 z1 cDAB ln 1 xA2 ln 1 xA1
N xA xB 1xB 1 xA
z 0 z dz C1 z2 z1 cDAB ln xB 2 ln xB1
dN Az 0 dz AB B2
c x
N C1 D
AZ
C ln
1
Por Ley de Fick N cDAB d xA z2 z1 x
B1
AZ
1 xA d z
cD d x cD AB x B2
C1 AB 1 A
N ln
xA d z
dx Az
z2 z1 x
C d zcD A B1
1 AB 1 x A

C1 z cDAB ln 1 xA C2
Condiciones
C1 z z1
xx
A A1

C2 z z2 xA xA2
C1 z1 cDAB ln 1 xA1 C2 4 C1 z2 cDAB
ln 1 xA2 C2 5

Se despeja C2 de ambas y se igualan

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

25. Efecto de la transferencia de masa en perfiles de concentración


a) Combine los resultados de las siguientes ecuaciones para obtener
1x A
N Az z 2
z 1

exp
1 xA1 cDAB

N
cD AB x B2

(1) ------- Az ln

z2 z1 x
B1

1 x A
1 x A2
z z 1

ln ln
1 xA1 1 xA1 z
2
z
1

1 x A
ln
1 x A2 1x
ln A1

1 x A1 z z 1

z
2 z1
De la ecuación 1
cD AB x B2
c AB
1 x A2

D
N ln ln
Az x
z2 z1 x
B1 z2 z1 1 A1

1 x A2
N Az z 2
z 1

ln
1 xA1 cDAB
Se igualan ambas
ecuaciones

b) Obtener el mismo resultado del inciso anterior integrando:


cD d x
AB
N
Az 1 xA A d
dx
z
N dz cD
Az AB 1 x
A
A Aplicando exponenciales
N Az z cDAB ln 1 xA C1 1x A N Az z z 1

exp ln exp
Condiciones 1 x
cDAB
A1

1 xA N Az z z1
z z1 xx exp
A A1
1 xA1 cDAB
C1 N Az z1 cDAB ln 1 xA1
N Az z cDAB ln 1 xA N Az z1 cDAB ln 1 xA1
A
1 x
N Az z z1 cDAB ln
1x
A1

1 x A
N z z Az 1

ln
1 xA1 cDAB
M. en C. María Guadalupe Ordorica Morales
Problemario de Fenómenos de Transporte

26. Determinación de la difusividad de un gas mediante dos bulbos (Análisis en estado


cuasiestático)
Una manera de medir las difusividades de un gas es mediante dos bulbos. El bulbo izquierdo y el
tubo desde z=-L hasta z=0 son llenados con gas A. El bulbo derecho y el tubo desde z=0 hasta z=+L
son llenados con el gas B. Al tiempo t=0, la valvula es abierta, y la difusión empieza, luego la
concentración de A en ambos bulbos cambia. Uno mide x A + en función del tiempo, y de esta
manera se deduce DAB. Se desea encontrar las ecuaciones que describan dicha difusión.

Ya que los bulbos son largos en comparación con el tubo, x A+ y xA cambian lentamente con el
tiempo. Por lo tanto la difusión en el tubo puede ser tratada como un problema de estado
cuasiestático, con las condiciones de frontera x A=xA- z=-L y xA=xA+ z=+L

- +
z=-L z=0 z=+L +
xA (t)
x =1-xA

a) Escriba un balance molar de A en un segmento Δz del tubo (de un área transversal S) y


demuestre que NAz=C1
1
S N Az
z
S N Az
z z 0
S z
N z z z AZ
lim Az N Az 0 dN 0
z 0 z dz
dN Az 0 d z
NAZ C1

b) Demuestre que la ecuación N Az cDAB


xA
xA N Az NBz se simplifica, para este z
problema en N Az cDAB d xA
d
z
0
De acuerdo con la Ley de Fick
N cD xA x N N
Az AB A Az Bz
z
N cDAB d xA
Az
dz
c) Integre la ecuación del inciso b), usando respuesta de a) obtenga C2
d x
C1 cDAB d zA

C1 d z cDAB d xA

C1 z cDAB xA C2
C2 C1 z cDAB xA

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

d) Evalue la constante con la condiciones de frontera.


xA xA z L
xA xA 1 xA z L
C2 C1 z cDAB xA
C 2 C L cD
1
AB
x
A

C2 C1 L cDAB 1 xA
C1 L cDAB xA C1 L cDAB 1 xA
C1 L C1 L cDAB 1 xA cDAB xA
2 L C1 cDAB 1 xA xA
2 L C1 cDAB 1 2 xA
C cDAB 1 2 xA
1
2L
cD
AB 1
C1 xA
L 2
N c
Az 1

N
cDAB 1
Az xA
L 2

M. en C. María Guadalupe Ordorica Morales


Problemario de Fenómenos de Transporte

Referencias

Bennett, O, Myers, J.E. Transferencia de cantidad de movimiento, calor y materia (tomos 1


y 2). Reverte, México, 2ª, edición,1998

Bird, R., Byron, W.E., Stewart, E.N., Lightfoot. Fenómenos de transporte, un estudio
sistemático de los fundamentos del transporte de materia, energía y cantidad de movimiento.
Reverte, México, 1ª. Edición, 1993.

Garcell Pyans, L. Transferencia de cantidad de movimiento, calor y masa. Ministerio


de educación Superior de Cuba – IPN. México 1998.

Geankoplis Ch. J., Procesos de transporte y principios de procesos de separación.


Compañía Editorial continental, cuarta Edición México, 2006.

Treybal, J.C. Operaciones de transferencia de masa. Mc. Graw Hill, Mèxico, 1980.

Welty, J.R. Fundamentos de transferencia de momento, calor y masa. Limusa, México 1972.

M. en C. María Guadalupe Ordorica Morales

También podría gustarte