l.

Martínez: TERMODINAMICA BASICA Y APLICADA

Transporte de masa

Se van a estudiar aquí los procesos de transporte de masa asociados a un gradiente de

concentración. Aunque se dice transporte de masa, se trata del transporte de especies
químicas por difusión molecular (y por convección másica en el capítulo siguiente), y no de
estudiar el movimiento medio de la masa global (típico de la Mecánica de fluidos).

La Termodinámica enseña que los sistemas tienden a evolucionar hacia el estado de

equilibrio, en el cual desaparecen las velocidades relativas, y la temperatura, y (si no hay
campos de fuerza exteriores) la presión y el potencial químico, son uniformes. La
Termodinámica de la evolución enseña que lejos del equilibrio pueden desarrollarse procesos
que tiendan a disminuir la entropía (p.e. difusión de especies desde bajas concentraciones a
altas), pero sólo si algún otro proceso acoplado genera más entropía de la que pierde el
primero.

La difusión es muy efectiva en las mezclas gaseosas y mucho menos en las líquidas, debido a
la gran diferencia de camino libre medio. En estado sólido la difusión es despreciable,
aunque se logran difusiones superficiales a alta temperatura (de utilidad en tratamientos
superficiales de materiales).

Recordemos que para medir la composición de una mezcla se puede utilizar las fracciones
molares xi, las fracciones másicas yi, las densidades pi, o las concentraciones ci, definidas
como:

El balance másico de una especie genérica i, considerada como un fluido con densidad pi y
velocidad local (macroscópica) Vi, en un volumen de control áxdydz, donde se supone que
puede haber una fuente o sumidero de especies wipor reacciones químicas, será:
Cap. 1 1 : DlFUSlON TERMICA Y MASICA 259

Sumando para todas las especies e introduciendo la definición de velocidad másica media v',
se obtiene la bien conocida ecuación del balance másico total (5.22) o (5.26):

*+cv-(~~c~)=o,
dt
ycon V I -
P'Cpici + - +aPatv . ( ~ v ) = o + - = -DP~~tv . v +
(1 1.21)

Nótese que análogamente se podría haber introducido la velocidad molar media

v" = cici / z c i , la cual tendría en general otro valor.
Ley de Fick

La relación más sencilla entre el flujo de masa de un componente y el gradiente de

concentración que lo ocasiona Vpi, fue propuesta por el fisiólogo alemán Fick, en 1855:
-. 1 - c ) p . c .
j.1 =p.(c.
1
r dr
o bien xiVdi= -DiVxi
O bien yiVdi = -DiVyi
o bien ciGi = -DiVci

donde se supone que sólo hay gradiente de la especie i, y no hay gradiente de otras especies,
ni de temperatura, ni de presión, ni campos de fuerza externos. El coeficiente de difusión Di
depende de la especie i, de la temperatura, de la presión (poco) y de las demás especies j. A
Vdi, definida por (1 1.22), se le llama velocidad de difusión. En primera aproximación, las
cuatro formas de (1 1.22) son equivalentes.

Ecuación de la difusión

A partir de (1 1.20), con ayuda de (1 1.22), se llega a:

= wi +V.(DiVpi)-piV.v

que también se puede escribir como:

En todo el desarrollo posterior se va a introducir la simplificación de considerar que la

variación de densidad global p debida al transporte de especies es despreciable (implica
V .v' = O), lo cual será válido para pequeñas concentraciones o cuando la difusión de una
especie se compense másicamente con la difusión de las otras. Además, se considerará que la
difusividad Di de una especie en una determinada mezcla es constante, con lo que se obtiene
finalmente:
 l. Martínez: TERMODINAMICA BASICA Y APLICADA

Nótese la gran similitud entre los procesos de difusión de especies y de difusión de energía
térmica (conducción de calor), lo que ha motivado el presentarlos en este mismo capítulo. El
a
equivalente a la ley de Fick (1 1.22) es la ley de Fourier, = -kVT y el equivalente a la
ecuación de la difusión (1 1.24) es la ecuación del calor, $/&=a V2T+@l(pc),deducida en
(1 1.3).

Como la variación de entropía por variación de composición es dS=Oli/T)dni la producción

de entropía
*
será Sgen= :v(pi I T ) , siendo el flujo molar de sustancia (mo1.m-2.s-l), y
como ji = ciGdi = -DiVci y para mezclas ideales V(//,/~IT)=RVX~/X~,
se llega a que S,,=
ciDi[V(/ilr)]2/R.

Difusión unidimensional plana

Este tema se va a desarrollar con un ejemplo: la difusión de vapor de agua desde una
superficie líquida hacia una atmósfera no saturada. Se supone que no hay convección (p.e.
sería aplicable a la evaporación en un tubo de ensayo y no en una cubeta con aire circulante),
que se ha alcanzado el régimen estacionario, que el proceso es isotermo y las
simplificaciones usuales (difusividad constante, gases ideales, aire no soluble en agua y sin
reacciones químicas).

Si se supone conocida la fracción másica de vapor y1 a una altura dada z l , (Fig. 11.5), y
sabiendo que a nivel del líquido z, la fracción másica será la de saturación, las ecuaciones de
balance másico total y de la especie vapor (i.e. vapor de agua) serán:

donde se ha eliminado la notación vectorial por ser el problema unidimensional, se ha

tomado p=cte y además pv=piVi es decir, se ha despreciado la velocidad propia del otro
componente (aire en calma). El objetivo es determinar la ley de evaporación, es decir,
calcular voy y(z) y dz,ldt (Fig. 11.5) en función de las propiedades de los fluidos y las
condiciones de contorno impuestas:
Cap. 1 1 : DlFUSlON TERMICA Y MASICA 261

(donde Mw=MvlMa=0,62es el cociente de masas molares del vapor y el aire), lo cual se hace
fácilmente con (1 1.26), obteniéndose:

y un
que en el caso normal de y<<l conduce a v ~ = D ~ ( ~ ~ - ~ ~ ) / (perfil
z ~ - de
z ~concentraciones
)
prácticamente lineal.

Fig. 11.5. Perfil de concenuación de vapor de agua en las

inmediaciones de agua líquida ('y' es la fracción másica) en
la difusión plana estacionaria en aire en calma.

Pero la evaporación del líquido requiere un aporte de energía qo=pvoAhlvo(siendo Ahlvola

entalpía de vaporización a la temperatura de la interfase), que bien puede llegar por radiación
desde un cuerpo lejano, o por conducción a través del líquido (o enfriándose éste), o por
conducción a través del gas. En este último caso, el balance energético aplicado a la mezcla
estacionaria (5.25), será:

donde se ha tenido en cuenta el estado estacionario, la ausencia de efectos viscosos, que el

aire está en calma ( p ~ ~ ~ p ~que ~se ha ~ llamado
= p v,=v ~ como
~ ~ antes,
~ ~y que ~ se~han~ ~ ~
introducido las condiciones de contorno térmicas en z, para la integración primera. Nótese
que el cp que aparece es sólo el del vapor, mientras que la k es la de la mezcla (prácticamente
la del aire). Integrando otra vez se llega a:

habiendo supuesto que el término entre corchetes es próximo a la unidad (en este caso el
perfil de temperaturas será lineal). Como se han supuesto conocidas las condiciones de
contorno en concentraciones, se ha determinado v, con (1 1.28), y con (1 1.30) se determinaría
el salto de temperaturas. Si se suponen conocidas las condiciones de contorno térmicas, de
(11.30) se obtendría v, y de (1 1.28) la diferencia de concentración. En ambos casos hay que
tener presente el equilibrio en la interfüse, dado por la primera relación en (1 1.27).
 l. Martínez: TERMODINAMICA BASlCA Y APLICADA

Difusión unidimensional esférica

Este problema tiene más aplicaciones prácticas que el anterior, pues sirve para modelizar la
evaporación de combustibles en carburadores (motores de gasolina), inyectores en cámara
cerrada (motores de gasolina y de gasoil), inyectores en cámara abierta (turbinas de gas,
calderas y hornos), así como la evaporación de gotitas de agua (de gran interés en
Meteorología) y de otras sustancias de interés en procesos químicos (extinción de fuegos,
secado de leche en polvo, etc).

Con las mismas hipótesis del apartado anterior, y la geometría de la Fig. 11.6, las ecuaciones
de conservación de la masa y las especies dan:

que, con ayuda de las condiciones de contorno impuestas:

da lugar a los siguientes resultados:

ln--- - yo - varo - -- P
dro -
- cte
v0 = - (1 1.34)
l-ym Di dt plíquido o'

De esta última expresión se deduce que el área de la gota disminuye linealmente con el
tiempo. En el límite y<<l (que es válido si las temperaturas no son muy altas), vo=
Di(yo-yw)r0, el perfil de concentración decae hiperbólicamente, y el tiempo que tarda en
evaporarse toda la gota (Langmuir-19184) es:

Fig. 11.6. Perfil de concentración de vapor de líquido en la evaporación estacionaria

de una gota en aire en calma ('y' es la fracción másica).

4. Langmuir, I., Phys. Rev. 12, p 368, 1918.

Cap. 1 1 : DlFUSlON TERMICA Y MASICA

Para el balance energético se procede como en (1 1.29), obteniéndose:

donde se ha tenido en cuenta que vr2=cte. Integrando otra vez, se obtiene:

habiendo hecho la misma linealización que en (11.30), aunque en este caso el perfil de
temperaturas a que da lugar ya no varía linealmente sino hiperbólicamente con el radio.
Como se ha supuesto conocidas las condiciones de contorno en concentraciones, se ha
determinado vo con (1 1.34), y con (1 1.37) se determinaría el salto de temperaturas. Si se
suponen conocidas las condiciones de contorno térmicas, de (1 1.37) se obtendría vo y de
(11.34) la diferencia de concentración. En ambos casos hay que tener presente el equilibrio
en la interfase, dado por la primera relación en (1 1.27). El tiempo que tarda la gota en
evaporarse totalmente viene dado por (1 1.35), aunque puede ponerse en función del salto de
temperaturas en lugar del de concentraciones. Si no es apropiada la linealización asumida en
(1 1.35), P.e. porque las temperaturas son altas, puede usarse:

En la Fig. 11.7 se representa la temperatura que alcanza el líquido en una gota de agua en
presencia de una atmósfera de aire, y el tiempo que tarda en evaporarse hasta la mitad del
radio, en función de la humedad relativa ambiente.

Fig. 11.7. Evaporación por difusión de una gota de agua pura de 0,l mm de radio en aire en calma a 20 "C. La
temperatura del líquido en el equilibrio de enfriamiento por evaporación T ~ í q , ~esqindependiente
, del
radio, y el tiempo que tarda en disminuir el tamaño a la mitad t l e s~ proporcional al cuadrado del
radio @.e. para un radio de 1 mm, con @=O, el tiempo pasa a ser de 5000 S).
 l. Martínez: TERMODlNAMlCA BASlCA Y APLICADA

RECAPITULACION

1. Se presentan los tres mecanismos clásicos de transmisión de calor: conducción,

convección y radiación y las leyes básicas que los rigen: ley de Fourier, ley de Newton y
ley de Stefan-Boltzmann.

2. Con ayuda del balance energético y la ley de Fourier se deduce la ecuación del calor,
Zldt = aV2T, analizando en detalle sus diferentes condiciones de contorno.

3. Se estudian los problemas elementales de conducción unidimensional y estacionaria en

pared plana, cilíndrica y esférica, y en particular el fenómeno del radio crítico de
aislamiento.

4. Se presentan algunos problemas elementales de conducción unidimensional no

estacionaria (excitación periódica y en escalón), de tanta importancia práctica (y para la
metrología térmica).

5. Se da una descripción general de los métodos numéricos que en la realidad se usan para
resolver todos los problemas prácticos de conducción.

6. Se analiza el problema de la difusión de especies en una mezcla con gradientes de

concentración, y se ve que es matemáticamente semejante al de la conducción de calor.

7. Se resuelven dos problemas sencillos de transferencia de masa: la difusión

unidimensional plana y la esférica.

PROBLEMAS

11.1. Determinar el flujo de calor y el perfil de temperatura en una pared plana compuesta
de 10 cm de ladrillo de k=0,5 W.m-l.K-l y 1 cm de chapa de hierro de k=50
W.m-1.K-1, sabiendo que las temperaturas extremas son de 1200 K en el lado del
ladrillo y de 300 K en el otro lado.
Sol.: q=45 kW/m2, T=309 K.

11.2. Determinar el flujo de calor y el perfil de temperatura en una pared plana compuesta
que representa el suelo de una bodega de un buque frigorífico de 5 mm de hierro de
espesor del casco, 10 mm de cámara de aire, 20 mm de espesor de madera, 30 mm de
espesor de corcho y 1 mm de pared de cinc, sabiendo que el agua está a 10 "C y la
bodega a -20 'C.
Sol.: q=26,5 kW/m2, las temperaturas en las interfases son 10; 10; -6; -93 y -20 'C.

11.3. Una varilla cilíndrica de 1 cm de diámetro está rodeada de u n baño térmico a 130 "C.
El material de la varilla, de k=0,5 W.m-1.K-1,está sufriendo un proceso reactivo que
libera 10 W/cm3. Determinar el perfil de temperatura.
Sol.: q=50 kW/m2, perfil parabólico con el máximo T0=255 'C en el centro.
Cap. 1 1 : DlFUSlON TERMICA Y MASICA 265

11.4. Para mover un eje a 1000 rpm en un cojinete de fricción, de bronce, de 60 mm de

agujero, 100 mm de diámetro exterior y 500 mm de ancho, se necesita un par de 10
N.m para vencer el rozamiento. Calcular la potencia consumida y el perfil de
temperatura en función de la temperatura exterior.
Sol.: q=1,77 kWlm2, A T 4 , 9 'C.

11.5. Una aleta rectangular de aluminio tiene 5 mm de espesor, 50 mm de longitud y 100

mm de anchura. La temperatura en la raiz es de 130 OC y el coeficiente de
convección con el aire atmosférico es h=50 W.m-2.K-1. Determinar el perfil de
temperatura, el flujo de calor y el rendimiento de aleta.
Sol.: &=56 W, q=0,92.

11.6. Se trata de calcular el perfil de temperatura en una varilla de aluminio de 30 cm de

longitud y 3 mm de diámetro, uno de cuyos extremos está en contacto con aceite
hirviendo a 200 OC, en una atmósfera a 20 OC con la cual se estima que el coeficiente
de convección térmica es de 15 W.m-2.K-1.Se pide:
a) Establecer el balance energético de un trozo de varilla.
b) Determinar analíticamente el perfil de temperatura y el flujo total de calor.
c) Calcular la temperatura en el otro extremo y el flujo total de calor.
d) Generación de entropía
Sol.: dZTldK2-(hplkA)(T-T,)=O;b) (T-Tm)l(To-T,)=coshmx-tanhmL.senhmr, &=kA(T,-T,)
.m.tanhmL, con m2=(hp)l(kA);c) TL=38OC, ~ = 2 , W. 5

11.7. Para estudiar la influencia de la vestimenta sobre la transmisión de calor del cuerpo
humano se propone el siguiente ejercicio. Determinar la masa de aislante necesario
para mantener un cilindro vertical de 1,7 m de altura y 30 cm de diámetro a 20 OC por
encima de la temperatura ambiente, sabiendo que debe disipar 100 W, para los casos
de lana (k=0,05 W.m-l.K-l, p=100 kg.m-3, como modelo del pelo, y cuero (k=0,15
W.m-1.K-1,p=1000 kg.m-3) como modelo de la piel. Supóngase que la transmisión de
calor en el aire puede variar desde 5 W.m-2.K-1hasta 100 W.m-2.K-1si hay viento.
Sol.: Para h=5, ecuerO=18mm, mcuer,=29 kg, e1,,=6 mm, mr,,=l kg; para h=100,
eCuerO=48 rnrn,mcuero=77kg, e1,,,=16 mm, m1,,,=2,5 kg.

11.8. La pared de un horno de tratamientos térmicos tiene una conductividad térmica k=0,7
W.m-2.K-1, una densidad p=1200 kg.m-3, una capacidad térmica c=1130 J.kgl.K-1
y 20 cm de espesor. Para un cierto proceso se necesita que la temperatura interior
varíe cíclicamente entre 100 OC y 750 OC cada 3 horas. Se pide:
a) Modelización del problema.
b) Amplitud de la variación de temperatura en la cara externa.
c) Retraso de la onda térmica.
Sol.: a) se supone proceso periódico (cada 6 horas) y variación sinusoidal, y pared
térmicamente gruesa, despreciando la influencia externa;

b) AT = ATOexp -x
1 E11 - =22 OC;
 l. Martínez: TERMODINAMICA BASICA Y APLICADA

Por un hilo de acero de resistividad eléctrica Q.m y 1 mm de diámetro se hace

pasar una corriente eléctrica de 2 A. El hilo está forrado de un aislante de 1 mm de
espesor y k=l ~.m-1.K-l, y la convección con el aire atmosférico es de h=5
W.m-2.K-l. Se pide:
a) Generación de calor por unidad de volumen.
b) Flujo radial de calor.
c> Perfil de temperaturas.
Sol.: a) 9 =4,6. lo6 W.m-3; b) ~ l L = 3 , 6Wlm; C)Tia=365 K, Texl=364,4K.

Se trata de dimensionar el aislamiento de un tanque de almacenamiento de gas

natural licuado (GNL, pCNL=500kgIm3) de 20 000 toneladas de capacidad, a presión
atmosférica, para que la vaporización por transmisión de calor desde el ambiente no
supere el 0,04% por día. Para minimizar estas pérdidas se va a adoptar una forma
cilíndrica de área mínima (de radio R=(VI2n)ln, siendo V el volumen, que ha de ser
un 20% mayor que el ocupado por el GNL).
El aislante de la base es principalmente (aparte del soporte estructural) lana de vidrio
de k=0,04 W.m-l.K-l y 70 cm de espesor, y el de las paredes y el techo es perlita de
k=0,03 W.m-1.K-l. Se pide:
Diferencia de temperatura entre el exterior y el interior.
Máximo flujo de calor admisible.
Flujo de calor por el suelo.
Espesor requerido de aislante en paredes y techo.
Repetir c) añadiendo el efecto de la vasija donde está realmente contenido el GNL,
que es un vaso abierto de 50 cm de espesor .de hormigón.
f) Repetir d) con el efecto del hormigón.
Sol.: a) TGNL=l12 K, Td=288 K; b) ~ = 5 1 , 4kW; c) ~,,,~,=12,3 kW; d) e=0,82 m; e)
~ ~ ~ ~ ~kW; ~ = f) e=0,82
1 2 , m.
1

Una pared de hormigón de 0,5 m de espesor está en una atmósfera a 300 K y a partir
de un cierto instante empieza a recibir una irradiación de 500 Wlm2 que absorve en
su totalidad. Se pide:
Perfil de temperatura estacionario suponiendo que las pérdidas son por convección
con h=10 W.m-2.K-1.
Perfil de temperatura estacionario suponiendo que las pérdidas son por convección
con h=10 W.m-2.K-1 y por radiación como cuerpo negro.
c) Evolución del perfil de temperatura en régimen transitorio.
Sol.: a) T1=306,4 K, T2=343,2 K; b) T1=305,4K, T2=339,0 K.

Considérese una fuente volumétrica de energia de intensidad = 4, e x p [ p ( -

~ T,)]
por unidad de volumen. Determinar bajo qué condiciones se alcanzana un régimen
estacionario cuando el contorno está a T, en los casos siguientes:
Pared plana de espesor L y conductividad k.
Cilindro de radio R y conductividad k.
Cap. 11: DlFUSlON TERMICA Y MASICA

c) Esfera de radior R y conductividad k.

Sol.: a) #&'L2/k=0,88; b) @JR2/k=2; c) @JR2/k=3,32.

11.13. El núcleo de un reactor nuclear consta de varillas cilíndricas de óxido de uranio de

conductividad térmica k=2 W.m-l.K-l y 1 cm de diámetro, que liberan lo8 W.m3 de
energía nuclear, encapsuladas en una vaina cerámica de 2 mm de espesor y k=20
W.m-l.K-l. El coeficiente de convección forzada con el agua es de 5 kW.m-2.K-1.
Determinar el perfil radial de temperatura en la varilla.
Sol.: ATol=312 K, AT12=21K, AT12=36K.

11.14. Se trata del estudio del perfil de temperatura en régimen periódico en un estanque de
agua en reposo, en contacto con una atmósfera gaseosa cuya temperatura varía según
la ley T,=Trn+ATsen(2nt/~), con Tm=5 "C y un coeficiente de convección h=10
W.m-2.K-1, en función de AT y T. En particular, se desea conocer el AT mínimo para
la formación de hielo en los casos -1 hora, -1 día y -1 mes.
Sol.: AThra=370C,ATdía=ll 0C,ATms=60C.

11.15. El modelo de transmisión de calor para el calentamiento de agua en una olla de

aluminio sobre una placa de hierro (cocina) expuesta a la llama de un hogar es el
siguiente. Se supone flujo unidimensional (vertical); la parte inferior del hierro se
supone a 600 "C, el espesor de hierro es de 5 mm y su conductividad k=50
W.m-1.K-1; debido al mal contacto físico entre la placa y la olla, se supone que existe
una capa de 50 pm de aire (k=0,02 W.m-I.K-l) entre ellas; el espesor de aluminio de
la olla es de 2 nim y su conductividad k=200 W.m-1.K-1; el coeficiente de
convección térmica con el agua, que está en ebullición con una temperatura lejos del
fondo de 100 "C, se supone que es de 4000 W.m-2.K-l. Se pide:
a> Esquema del perfil de temperatura.
b) Cálculo del flujo unitario de calor.
c) Cálculo de las temperaturas intermedias.
d) Corrección correspondiente a tomar una dependencia de la forma klT1I2=cte para la
conductividad del aire.
Sol.: b) 9 ~ 1 7 5kW/m2; C) T1=582,5 "C, T2=145,5 "C, T3=143,7 OC; d) kaire(20 OC)=0,02
W.m-l.K-l, kaire(T1)=0,034W.m-1.K-1, kaire(T2)=0,024W.m-1.K-1, por lo que se
puede tomar un valor medio de 0,029 W.m-1.K-1 y volver a resolver, obteniéndose
q=234 kW/m2, T1=576 "C, T2=162 "C,y T3=160 "C, o bien resolver el problema con
k variable:

11.16. Se trata de analizar el régimen estacionario de una instalación de agua caliente

sanitaria que consta de u n crilentador, donde la corriente de agua recibe 20 kW, y de
una tubería de 1 cm de diámetro y 10 m de longitud cuya transmitancia térmica
global con el cambiente se supone que es U = 200 W.m-2.K-1 referido al área del
tubo. Se pide:
268 l. Martínez: TERMODINAMICA BASlCA Y APLICADA

a) Temperatura de salida del calentador, T1 (en función del gasto másico de agua).
b) Perfil de temperatura a lo largo del conducto.
c) Temperatura de salida del conducto, T2.
d) Gasto másico que daría máxima temperatura de salida del conducto y valor de ésta.
e) Comentar la viabilidad y conveniencia del punto de trabajo anterior.
Sol.: a) Ti = To + Q 1(mcp) = 288 + 20000/(h .4200), con en kg/s y T en K [ver Fig.
P-11.161; b) T(x) = To + [~I(mc,)]exp[-hnDx/(hcp)]; c) T2 = To +
[ ~ l ( m c ~ ) ] e x ~ [ - h [ver ) ] P-11.161; d) mlT2,mh = hnDL/cp =
~ ( ~Fig.
200.x0,01.10/4200 = 0,015 k g l ~ , T2,rnsr = To + Q l ( h n ~ ~ c=) 288 +
20000/(200.7r.O,Ol.lO.e) = 405 K; e) Como la presión de suministro doméstico es del
orden de 0,3 MPa, la TI habrá de ser menor de T,, = 410 K y la T2 correspondiente
será menor de 370 K, aunque realmente no se necesitará que sea mayor de unos
330 K y sí se deseará más caudal, así que el punto de trabajo puede ser el T2 = 330 K
y m = 0,09 kg/s (5,5 litros por minuto de agua a 57 OC).

m [kg/sl
Fig. P-11.16

11.17. Estimar el tiempo que tardaría en evaporarse 1 cm de alcohol etílico en el fondo de

un tubo de ensayo de 1 cm de diámetro y 10 cm de altura. Datos: Di=1P5 m2/s,
pvapor=10 kPa, plí, = 790 kg/m3.
Sol.: t = 36 horas.

11.18. En un cilindro vertical cerrado por ambos extremos y con una válvula intermedia
(célula de Loschmidt), se ha dispuesto un gas ligero en la mitad superior y uno más
pesado en la mitad inferior. Determinar la evolución del perfil de concentraciones al
abrir la válvula que separa ambos gases.
Sol.: y,=sen(mlL), siendo x=O el punto medio y L la longitud total del tubo.