XIV.

- CONDENSACIÓN Y VAPORIZACIÓN
[Link]

XIV.1.- TRANSFERENCIA DEL CALOR POR CONDENSACIÓN

La condensación se produce cuando un vapor saturado que se pone en contacto con una superficie a
menor temperatura, se enfría hasta que la temperatura se hace inferior a su temperatura de satura-
ción. Si el vapor es puro, la temperatura de saturación corresponde a la presión total; si se trata de una
mezcla de vapor y gas no condensable, la temperatura de saturación corresponderá a la presión parcial
del vapor.

CONDENSACIÓN EN FORMA DE PELÍCULA.- Si el condensado se forma sobre una pared ver-

tical y la humedece, en condiciones normales se forma sobre la superficie un flujo continuo de fluido flu-
yendo el condensado hacia abajo por la acción de la gravedad, aumentando el espesor de la película de
modo continuo de arriba hacia abajo. Salvo que la velocidad del vapor sea muy alta o que la película de
líquido sea muy gruesa, el movimiento del condensado es laminar, y se transfiere la entalpía de conden-
sación desde la intercara (vapor/líquido) hacia la superficie simplemente por conducción.

La velocidad del flujo de calor depende principalmente del espesor de la película del condensado, la
cual a su vez depende de la velocidad a la que condensa el vapor y de la velocidad conque se elimina el
condensado. Cuando la placa es inclinada respecto a la posición vertical, disminuye la velocidad del con-
densado y la película de líquido se hace más gruesa, lo que origina una disminución de la velocidad de
transferencia de calor.

Los coeficientes de transmisión de calor por convección para la condensación de vapores puros en
forma de película, sobre tubos y placas, fueron obtenidos por primera vez por Nusselt en 1916.

Proceso dinámico.- Para placa vertical, a una distancia x de la parte superior de la placa de an-
chura unidad, el espesor de la película es δ; si el flujo de líquido es laminar y está producido sólo por la ac-
ción de la gravedad, se puede estimar la velocidad del líquido mediante un balance de fuerzas sobre el ele-
mento de volumen dx, Fig XIV.1.

La fuerza Fl que actúa sobre el líquido contenido en el volumen de espesor dx, a una distancia de la su-
perficie entre y y δ es:

Fl = ( δ - y ) dx ρ l g
[Link] Condensación y vaporizació[Link].-273
 Fig XIV.1.- Condensación en forma de película sobre una superficie vertical

Suponiendo que el vapor que está fuera de la capa del condensado se encuentra en equilibrio hidros-
tático, se tiene que:

dp dp dp
dx = = ; = ρv g
γv ρv g dx

Como consecuencia de este gradiente de presiones, se puede interpretar que el elemento de condensado
(δ - y) se encuentra en una atmósfera de vapor saturado, por lo que aparecerá sobre el mismo una fuer-
za sustentadora Fv (prácticamente despreciable), igual al peso del volumen de vapor desalojado de la
forma:

Fv = ( δ - y ) dx ρ v g

Las fuerzas viscosas frenan el movimiento hacia abajo (rozamiento); por lo que respecta al vapor, a
menos que fluya a una velocidad muy elevada, las tensiones de cortadura en la superficie libre son muy
pequeñas, por lo que el rozamiento debido al vapor se puede despreciar; la fuerza de rozamiento genera-
da Fr será debida únicamente a la viscosidad del líquido, y se opone al deslizamiento del condensado de la
forma:

Fr = ηl du dx
dy

En condiciones estacionarias las fuerzas hacia arriba y hacia abajo son iguales, por lo que:

( δ − y ) ( ρ l - ρ v ) g = η l du
dy

⎧δ el espesor de la capa de condensado

siendo: ⎨ ρ l la densidad del líquido; ρ v la densidad del vapor; η l la viscosidad del líquido
⎩u = u( y) la velocidad a la distancia x

La distribución de velocidades en la película de condensado se obtiene integrando la anterior:

( ρ l - ρv ) g y2
u( y ) = (δ y - )
ηl 2

La velocidad alcanza un valor máximo sobre la superficie de la película (y = δ), en la forma:

( ρ l - ρv ) g δ 2
umáx =
2 ηl

[Link] Condensación y vaporizació[Link].-274

 El flujo másico de condensado por unidad de anchura de la placa es:

δ δ ρ l ( ρl − ρv ) g y2 ρ ( ρ − ρv ) δ 3 g ( ρl − ρv ) δ 3 g
G= ∫0 ρl u dy = ∫0 (δ y - ) dy = l l =
ηl 2 3 ηl 3 νl

Proceso térmico.- Como la velocidad del flujo másico del condensado G (kg/seg) sobre la placa de-
pende de la velocidad a la que condensa el vapor, si se supone que el flujo de película es de tipo laminar y
paralelo a la superficie, y que el gradiente de temperatura es lineal de la forma:

dT = Ts - TpF
dy δ(x)

y dado que el calor se transfiere a través de la capa de condensado, únicamente por conducción, se tiene:

Ts - TpF Ts - TpF
dq = kl dx dT = kl dx = dG rl-v ⇒ dG = kl dx 1
dy δ(x) δ(x) rl-v

en la que dG es la cantidad de vapor condensado en el elemento, (δ - y) dx:

ESPESOR DE LA CAPA DE CONDENSADO.- Para hallar el espesor de la capa de condensado se

igualan los valores de dG :

dx

dG = kl Ts - TpF =
g ρl ( ρ l - ρv ) d( δ 3 )
=
g ρ l ( ρ l - ρv ) δ 2 d δ
dx rl-v δ ( x ) 3 ηl dx ηl dx

k1 Ts − TpF 3 k1 Ts − TpF x δ
ηl
rl-v g ρl ( ρ l − ρv )
dx = δ dδ ; ηl
rl-v g ρ l ( ρ l − ρv )
∫0 dx = ∫ 0 δ 3dδ

k1 Ts - T pF 4 4 kl ηl(Ts - T pF ) x
ηl x=δ ⇒ δ= 4
€ rl-v g ρl (ρl − ρv) 4 rl-v g ρ l (ρ l - ρv)

que dice que el espesor de la película de condensado aumenta en función de la raíz cuarta de la distancia
recorrida x a lo largo de la superficie.

COEFICIENTE DE TRANSFERENCIA DE CALOR.- El calor evacuado en la condensación por

conducción, en el elemento de volumen de espesor dx a la distancia x, es el mismo que el evacuado por
convección de la forma:

Ts - TpF
dq = hC x dx (Ts − T pF ) = kl dx dT = kl dx
dy d

por lo que el coeficiente de convección local por unidad de anchura hC x, es:

Ts - TpF 1 k rl-v g ρ l(ρl - ρ v) kl3

hC x = k1 = 1 = 4
δ Ts - TpF δ 4 η l(Ts - TpF ) x

hC x x rl-v g ρ l ( ρ l - ρ v ) x 3
Nu x = = 4
k1 4 η l k1 ( Ts - T pF )

[Link] Condensación y vaporizació[Link].-275

 El aumento del espesor de la película de condensado es semejante al crecimiento de la capa límite so-
bre una placa plana en la convección; un aumento de la diferencia de temperaturas (Ts - TpF) produce
una disminución de la conductancia superficial; ésto está originado por el aumento del espesor de la pelí-
cula, como resultado del incremento de la velocidad de condensación.

El valor medio del coeficiente de convección hC correspondiente a una placa vertical de anchura uni-
dad y altura L es:

1 L rl-v g ρ l ( ρ l - ρ v ) k13 4
hC =
L
∫0 hCx dx = 0,943 4
ηl (Ts - TpF ) L
= h
3 CL

rl-v g ρ l ( ρ l - ρv ) L3
Nu = 0,943 4
€ ηl ( Ts - TpF ) kl

Aunque el análisis anterior está hecho específicamente para una placa vertical, el desarrollo es váli-
do para las superficies interior y exterior de tubos verticales, si éstos tienen diámetros grandes en com-
paración con el espesor de la película; estos resultados se pueden extender también a placas inclinadas
un ángulo θ, respecto al plano horizontal, sustituyendo g por g sen θ.

Tabla XIV.1.- Valores del coeficiente de transmisión de calor por convección en la condensación de vapores puros
Fluido Materiales Ts - TpF W/m2 ºC
Agua Tubos horizontales de diámetro exterior (25-75) mm 3-20 11.000-23.000
Agua Superficie vertical de altura 3 m 3-20 5.700-11.000
Etanol Superficie vertical de altura 0,15 m 10-55 1.100-1.900
Benceno Tubo horizontal de diámetro exterior 25 mm 15-45 1400-2.000
Etanol Tubo horizontal de diámetro exterior 50 mm 5-20 1.700-2.600
Amoniaco Anillo horizontal de 50 a 75 mm 1-4 1.400-2.600

XIV.2.- CONDENSACIÓN EN PELÍCULA LAMINAR SOBRE PLACAS Y TUBOS VERTICA-

LES

Determinación del tipo de régimen del flujo en la película.- Se pueden considerar 3 regímenes
de flujo en película, laminar, ondulatorio y turbulento.

Para números de Reynolds bajos el flujo es laminar y la superficie de la película presenta un aspecto
liso; a medida que aumenta el nº de Reynolds se forman ondas en la superficie de la película; al seguir au-
mentando aún más el nº de Reynolds estas ondas toman una forma compleja ondulatoria en tres dimen-
siones.

Las ondas hacen que el líquido se mezcle ligeramente, pero el flujo en la base sigue siendo laminar,
hasta que a velocidades relativamente altas, el flujo se vuelve turbulento por toda la película debido a la
inestabilidad originada por los esfuerzos cortantes.

El número de Reynolds de una película descendente se puede definir en función de la velocidad uF y

del diámetro hidráulico dh de la película:

uF ρ l 4 A hC A ( Ts - TpF )
u ρ d p Q
Re = F l h = = G = = = u F A ρl = 4G
ηl ηl rl-v rl-v p ηl

[Link] Condensación y vaporizació[Link].-276

siendo G el gasto másico.

Para tubos verticales: A = π d L ; p = π d ; Re = 4 G = 4 G

p ηl π d ηl

Para tubos horizontales: A = π d L ; p = L ; Re = 4 G = 4 G

p ηl L ηl

⎧W anchura
Para placas inclinadas: A = W L ; p = W ; Re = 4 G = 4 G , con ⎨ L longitud del tubo
p ηl W ηl ⎩

En la condensación, las propiedades del líquido condensado se evalúan a la temperatura media entre
la temperatura de saturación Ts y la de la pared TpF, mientras que las propiedades del vapor y del calor
latente de condensación rl-v se evalúan a la temperatura de saturación Ts.

a.1) Condensación sobre placa vertical y tubos verticales, en régimen laminar.- La expre-
sión para placa vertical se puede extender tanto a las superficies interiores como a las exteriores de tu-
bos verticales, siempre que su diámetro sea muy grande en comparación con el espesor de la película.
El tubo se puede considerar como una placa vertical de área (π d L) y perímetro (π d).

rl-v g ρ l (ρl - ρv) kl3 rl-v g ρ l2 kl3

La solución analítica es: hC(vertical)= 0,943 4 ≅ 0,943 4
η l (Ts - TpF ) L η l (Ts - T pF ) L

Resultados experimentales han demostrado que esta ecuación es conservadora, de forma que los re-
sultados obtenidos con ella están un 20% por debajo de los valores medidos, por lo que se propone la si-
guiente ecuación experimental:

rl-v g ρl ( ρ l - ρ v ) kl3 rl-v g ρ l2 kl3

hC(vertical)= 1,13 4 ≅ 1,13 4
η l ( Ts - T pF ) L η l (Ts - TpF ) L

El número de Reynolds del condensado en la parte inferior del tubo vertical es:

4 hC L ( Ts - T pF ) g ρ l ( ρ l - ρv ) kl3 ( Ts - TpF )3 L3
Re = = 4,52 4 < 1800
η l rl-v ηl5 rl-v
3

a.2) Otras expresiones para la condensación sobre placa vertical y tubos verticales, en ré-
gimen laminar

Para el caso particular del agua sobre tubos verticales, el flujo se vuelve:
- Laminar ondulatorio para ReL ≈ 30
- La transición a flujo turbulento en la región exterior de la película se da para ReL ≈ 1000
- La turbulencia se completa para ReL ≈ 1800

⎧ 4 4 L Ja l 3 g
⎪ Re L = 3 4 ( Prl ) v 2
Para ReL < 30, la película es laminar: ⎨ 1
c pl ( Tsat - TpF )
⎪ Jal = Nº de Jakob =
⎩ rl-v

El número de Nu local para el agua es:

3 3 Re
Nulaminar = ; 0 < Re < 30
4
[Link] Condensación y vaporizació[Link].-277
 El número de Nu medio es:

Prl 1 v12
Nu = 4 4 3
3 4 Jal L g

a.3) Condensación sobre placa vertical o tubos verticales, en película laminar ondulato-
ria.- En película laminar ondulatoria: 30 < ReL < Retrans, siendo el número de Nu local para el caso parti-
cular del agua:

Nulaminar ondul. = 0 ,822 Re -0 ,22 ; 30 < Re < Retrans.

en la que la transición a la turbulencia sucede, para el agua, cuando:

Retrans. = 5800 Prl-1,06

3
v12
4 Jal L ) 0 ,82 ⇒ Nu = { Prl g 0 ,18
Para el caso general: Re L = ( }
Prl v12 4 Jal L
3
g

EFECTO DEL SUBENFRIAMIENTO DEL LÍQUIDO.- En el caso de líquidos con calores latentes
de cambio de estado relativamente bajos, como muchos refrigerantes, es necesario hacer una corrección
para tener en cuenta el efecto del subenfriamiento; la correlación propuesta por Rohsenow tiene en
cuenta esta variación de la entalpía media del vapor al condensar y enfriarse posteriormente a la tem-
peratura media de la capa de condensado, para lo cual sustituye el calor latente rl-v a Ts por otro de la
forma:

⎪⎧ Pr > 0 ,5
* = r +
rl−v l-v 0 ,68 c pl ( Ts - T pF ) , para: ⎨ Ja = c Ts - T pF < 1 ( Rohsenov )
⎪ l pl *
⎩ rl−v

* = r + ( 0 ,683 - 0 ,228 0 ,228

rl−v l-v ) c pl ( Ts - TpF ) = rl-v {1 + ( 0 ,683 - ) Jal }
Prl Prl

El coeficiente de convección es:

* g ρ ( ρ - ρ ) k3
rl-v l l v l
hC = 0,943 4
η l ( Ts - TpF ) L

EFECTO DEL SOBRECALENTAMIENTO DEL VAPOR EN LA CONDENSACIÓN.- Para ha-

llar la influencia del sobrecalentamiento del vapor en la condensación sobre el coeficiente de transferen-
cia de calor hC , se puede utilizar la misma expresión obtenida para un tubo horizontal que es de la for-
ma:

* g ρ ( ρ - ρ ) k3
rl-v l l v l * = r +c
hC = 0,725 4 ; rl-v l-v pl ( Tsobrec - Ts ) + 0 ,35 c pl (Ts - TpF )
η l ( Ts - TpF ) d

El efecto del sobrecalentamiento del vapor aumenta el coeficiente de convección y, por lo tanto, la
transferencia de calor, pero disminuye la cantidad del condensado.

[Link] Condensación y vaporizació[Link].-278

XIV.3.- CONDENSACIÓN EN PELÍCULA LAMINAR SOBRE PLACAS Y TUBOS INCLINADOS

En la condensación sobre placas inclinadas un ángulo θ respecto a la horizontal se utilizan expresio-

nes para hC y Re similares a la anterior de la forma:

rl-v g ρ l ( ρ l - ρ v ) kl3 sen θ

hC sup. inclinada = 1,13 4
η l ( Ts - T pF ) L

g ρ l ( ρ l - ρ v ) kl3 ( Ts - TpF ) 3 L3 sen θ

Re = 4,52 4 < 1800
η l5 rl-v
3

XIV.4.- CONDENSACIÓN EN PELÍCULA LAMINAR SOBRE UN TUBO HORIZONTAL

a) El coeficiente de convección medio de un vapor puro saturado que está condensando so-
bre el exterior de un tubo horizontal, de forma que el espesor de película sea nulo en la parte supe-
rior del tubo, es:

4 hC dx ( Ts - TpF ) g ρ l ( ρ l - ρ v ) kl3 ( Ts - TpF ) 3 d 3

Re = = 9,11 4 < 3600
η l rl-v ηl5 rl-v
3

rl-v g ρ l ( ρ l - ρv ) kl3 rl-v g ρ l2 kl3

hC horiz = 0,725 4 = 0,725 4
ηl ( Ts - TpF ) d η l ( Ts - T pF ) d

1 = 1 ln ( r2/r1 ) 1
+ +
UA 2 π r1 L hCF1 2 πkL 2 π r2 L hC horiz

b) Condensación sobre una batería de N tubos horizontales, en régimen laminar.- Si la con-

densación se produce sobre N tubos horizontales dispuestos de tal modo que el condensado de un tubo
cae directamente sobre el tubo que tiene debajo, que es lo que sucede en la mayoría de los condensado-
res, se puede estimar la conductancia superficial para el conjunto mediante la expresión anterior, susti-
tuyendo el diámetro d por, Nd; este método proporciona resultados conservativos, porque es inevitable
una cierta turbulencia; para régimen laminar, Re < 3600:

rl-v g ρ l2 kl3 h
hC( N tubos ) = 0,725 4 = C( 11/4
tubo )
ηl (Ts - TpF ) d N N

Una expresión que tiene en cuenta que todo el calor del subenfriamiento se utiliza para una conden-
sación adicional de la capa límite entre N tubos horizontales, colocados unos encima de otros, es de la
forma:
rl-v g ρ l2 kl3
hC( N tubos ) = 0,728 {1 + 0 ,2 Ja1 ( N - 1 )} 4
ηl ( Ts - TpF ) d N

que concuerda bastante bien con los resultados experimentales, siempre que:

(N - 1) Ja1 < 2

Relación entre los coeficientes de convección hC horiz y hC vert.- Para una diferencia de tempe-
raturas determinada, el coeficiente de convección es mayor cuando se coloca el tubo en posición horizon-
[Link] Condensación y vaporizació[Link].-279
tal, que cuando se coloca vertical, ya que el camino recorrido por el condensado es menor, resultando la
película más delgada, por lo que:

hChoriz 0,725
= 4 L = 0,77 4 L
hCvert 0,943 d d

que implica el que para una diferencia de temperaturas dada, Ts - TpF, resulte:
hC horiz
Para: L = 2,86 d ⇒ = 1
hCvert
Para: L = 100 d ⇒ hC horiz = 2 ,43 hCvert

por lo que la disposición de tubos condensadores horizontales se prefiere a la vertical.

g ρl ( ρ l - ρv ) kl3 (Ts - TpF )3 N 3 d 3

Re = 4 G = 9,11 4 < 3600
ηl L N η 5l rl-v
3

XIV.5.- CONDENSACIÓN EN RÉGIMEN TURBULENTO

a) El régimen se considera turbulento cuando Re > 1800. Es difícil que se llegue a alcanzar el flujo
turbulento en un tubo horizontal durante la condensación, pero sí es posible se pueda conseguir en la parte
más baja de un tubo vertical.

Cuando ésto sucede, el coeficiente de transferencia térmica hC crece al aumentar la longitud L de la

superficie condensada, debido a que también aumenta el grado de turbulencia.

El número de Reynolds para superficies verticales isotérmicas en régimen turbulento es:

g ρ l ( ρ l - ρv ) kl3 (Ts - TpF )3 L3

Re = 0,00296 { } 5/9 > 1800
η 5l rl-v
3

Cuando por el exterior de la pared el vapor está en reposo, o cuando por el interior del tubo la velocidad
del vapor es muy pequeña, el coeficiente de convección local correspondiente al flujo turbulento del conden-
sado en placa vertical se determina mediante la relación:

g ρl2 kl3
hC x = 0,056 ( 4 u )0,2 3 Prl
ηl η l2

siendo u la velocidad en la base, del orden de: Rebase = 4 u

ηl
Si el vapor no está en reposo, y el movimiento del mismo no perturba la formación de la película, se
mantiene la formulación anterior.

Ts + TpF
Las propiedades del condensado se evalúan a la temperatura: T =
2
b) El número de Nu local para el caso particular del agua es:

0 ,65
Nuturb = 3 ,8.10-3 Re 0 ,4 Prl ; Retransv < Re

Caso general:

[Link] Condensación y vaporizació[Link].-280

 L = xtrans + Prl0,35
( Re0L ,6 - Retrans
0 ,6
-3 ) Ja
)
2
vl v 2 2 ,4 ( 3 ,8.10 l
3 3 l
g g

Número de Nu medio:

3
vl2
g Prl 9,12.10-3 Ja l ( L - xtrans ) 0,6
Nu = { + Retrans }10/6
L 4 Jal v 2
3 l Prl0,35
g

vl2 Prl c pl ( Tsat - T pared )

siendo: xtrans = 3 Re 1,22
trans ; Jal =
g 4 Jal rl-v

Efecto de la velocidad del vapor en placas y tubos verticales

a) Si la velocidad del vapor no condensado es considerable frente a la del líquido condensado, se puede
utilizar la siguiente ecuación:

hC λ ρl
Prl = 0,046
c pl G m ρv

siendo λ el coeficiente de rozamiento correspondiente a la velocidad media del vapor, viniendo dado el
gasto másico del vapor en la forma:

2 +G
Ginf 2
inf Gsup + G sup
Gm =
3

en la que Ginf es el gasto másico del vapor en la parte inferior del tubo y Gsup es el gasto másico del va-
por en la parte superior del tubo.

Las propiedades del líquido se calculan a la temperatura: T = 0 ,25 Ts + 0 ,7 5 TpF

b) Una expresión que se aplica a un vapor saturado seco a la presión de 1 atm, que desciende por un
tubo vertical a una velocidad comprendida entre 40 y 80 m/seg, es:

3 1,21
hC = 3400 + 100 uFv
L

siendo L la longitud del tubo en metros y uFv la velocidad del vapor en m/seg.

c) Otra expresión en función de los números de Reynolds Rel y Rev es:

⎧ 4 G1
hC d ρ l 0,8 ⎪ Rel = π d η ; Rel > 5.000
3 l
= 0,026 Prl ( Rel + Rev ) ; ⎨ 4 Gv
kl ρv
⎪ Rev = π d η ; Rev > 20.000
⎩ v

d) Otra expresión en función del número de Jakob, que tiene en cuenta el efecto del arrastre del vapor a
una distancia x del borde superior del tubo vertical, a la velocidad uv es:

[Link] Condensación y vaporizació[Link].-281

 kl2 u v 16 Prl g x Prl *
ρ l ν l rl-v
hC = (1 + 1+ ) con: =
8 νl x Ja l uv2 Jal kl ( Tsat − T pF )

En ausencia de arrastre de vapor:

* g ρ k3
rl−v l l
hC = 4
4 x (Tsat - T pF ) ν l

* = r + ( 0 ,683 - 0 ,228 0 ,228

rl−v l-v ) c pl ( Ts - TpF ) = rl-v + {1 + ( 0 ,683 - ) Jal }
Prl Prl

observándose que el efecto de arrastre del vapor consiste en aumentar el coeficiente de transferencia de
calor.

XIV.6.- CONDENSACIÓN EN RÉGIMEN TURBULENTO EN EL INTERIOR DE TUBOS HO-

RIZONTALES

En aplicaciones prácticas de condensadores en refrigeración y en sistemas de aire acondicionado, se

observa que el vapor condensa en el interior de los tubos con una cierta velocidad; el fenómeno de la con-
densación y su posible tratamiento analítico, en forma simple, es complicado.

El valor del coeficiente de transmisión de calor para la condensación de refrigerantes en el interior de

tubos horizontales, si se considera que el vapor tiene una velocidad despreciable, es:

* g ρ 2 k3
rl-v
hC = 0,555 4 l l * = r + 3 c
con: rl-v ( T - T pF )
ηl ( Ts - TpF ) d l-v 8 pl s

ρ v uv d
Revapor = < 35.000
ηv

XIV.7.- CONDENSACIÓN EN FORMA DE GOTAS

Cuando una superficie sobre la que va a condensar un vapor está contaminada con una sustancia
que impide que el condensado moje la superficie, el vapor condensará en forma de gotas, en lugar de ha-
cerlo como una película continua, fenómeno que se conoce como condensación en forma de gotas.

En estas condiciones una gran parte de la superficie no se ve cubierta con una película aislante y los
coeficientes de transferencia de calor pueden ser de cuatro a ocho veces más elevados que en la conden-
sación en forma de película.

Hasta ahora, la condensación en forma de gotas sólo se ha obtenido de modo fiable, con vapor de
agua. Para calcular el coeficiente de transmisión de calor por convección, se recomienda suponer una
condensación en forma de película porque, incluso en el caso del vapor de agua, sólo se puede esperar que
se obtenga la condensación en forma de gotas bajo condiciones cuidadosamente controladas, que no pue-
den mantenerse siempre en la práctica.

XIV.8.- TRANSFERENCIA DE CALOR POR EVAPORACIÓN DE LÍQUIDOS EN REPOSO

EVAPORACIÓN EN PELÍCULA DESCENDENTE SOBRE UNA PARED VERTICAL.- La

[Link] Condensación y vaporizació[Link].-282

 transferencia de calor a través de una evaporación en película descendente es casi idéntica a la conden-
sación en película descendente, salvo que se lleva en dirección contraria, es decir, la película siempre tie-
ne un espesor finito al comienzo disminuyendo el número de Re al descender por la pared a medida que
va progresando la evaporación; en determinadas situaciones la velocidad de evaporación es pequeña
comparada con la velocidad del flujo de la película y la variación del número de Re podría ser pequeña.

Como la temperatura TpF > Ts el líquido que se encuentra más cerca de la pared estará sobrecalen-
tado, y comienzan a formarse burbujas por nucleación sobre la pared iniciándose la ebullición, o tam-
bién, las burbujas contenidas en el líquido pueden crecer de forma explosiva (cavitación). Sea cual fuere
el mecanismo, lo cierto es que el proceso de transmisión de calor se complica.

Una película evaporativa puede ser, en principio, laminar ondulatoria o turbulenta. Si el arrastre de
vapor es mínimo, se pueden utilizar los siguientes números de Nu locales:

Nulaminar = 3 4 ; 0 < Re < 30

3 Re

Nulaminar ond.= 0 ,822 Re −0 ,22 ; 30 < Re < Retrans

Nuturbulento= 3 ,8.10-3 Re 0 ,4 Prl0 ,65 ; Re transv< Re, siendo: Retransv = 5800 Prl-1 ,06

El número de Reynolds se obtiene a partir de:

L Prl Re L d Re Prl ReL d Re Prl 1,22 1,22

= -
4 Ja l
∫ Re Nu
= -
4 Jal
∫ Re -0,22
= -
4 Ja l
(Re L - Re0 )
3
vl2 0 0 0,822 Re
g

4 Jal L g
Re 1,22 1,22
L = Re0 - 3
Prl vl2
€
en la que el número de Re inicial de la película es:
4 G0
Re0 =
π d ρl ν l
siendo G0 la cantidad de líquido inicial disponible para la evaporación en película descendente.

Prl 3 vl2/g
El nº de Nu medio es: Nu = ( Re 0 - Re L )
4 Jal L

EVAPORACIÓN NUCLEADA EN RECIPIENTES CON UN LÍQUIDO EN REPOSO.- Estudiar

el proceso de la evaporación nucleada en un líquido en reposo es extremadamente difícil, ya que en la ac-
tualidad no se pueden explicar muy bien ni el mecanismo exacto de la formación de burbujas, ni el movi-
miento de las mismas, totalmente aleatorio; por lo tanto no existe un análisis satisfactorio de la transfe-
rencia de calor en este tipo de evaporación; el mecanismo principal de la transferencia de calor es la con-
vección natural en régimen turbulento, causada por la agitación producida por las burbujas.

Para su estudio se puede dividir el fenómeno en las siguientes partes:

Convección libre (Régimen I, Fig XIV.2): La ecuación general es:

Q
= C k (Gr Pr )n ( Ts - T pF ), en la que n y C son constantes que se toman de las Tablas XIV.3.
A L
[Link] Condensación y vaporizació[Link].-283
 Fig XIV.2.- Vaporización nucleada en recipientes con líquido en reposo

Evaporación nucleada (Régimen II y III, Fig XIV.2).- A medida que la temperatura TpF de la pa-
red bañada por el fluido va aumentando, aparecen pequeñas burbujas de vapor que ascienden, de forma
que unas condensan en el propio líquido y otras llegan a la superficie libre, en donde revientan y sueltan
el vapor; la tensión superficial del líquido ofrece una gran resistencia a la formación de burbujas dentro
del propio líquido, por lo que inicialmente éstas se forman a partir de puntos de nucleación sobre la su-
perficie de la pared, en la que existen pequeñas imperfecciones en su acabado, o en bolsas de gas, en las
que los efectos de tensión superficial se hallan minimizados.

Tabla XIV.3.- Valores de la tensión superficial

líquido-vapor para varios líquidos
Líquido Temp. saturación Tensión superficial
ºC ( 1000 σ *) N/m
Agua 0 75,6
Agua 15,56 73,2
Tabla XIV.2 Agua 20,00 72,6
Valores del coeficiente C para distintas superficies-líquidos Agua 37,78 69,7
Alcohol n-butílico-cobre 0,0030 Agua 40,00 69,4
Alcohol isopropílico-cobre 0,0025 Agua 60,00 66,0
Alcohol etílico-cromo 0,0027 Agua 80,00 62,5
Agua-acero inoxidable 0,0133 Agua 93,34 60,1
Agua-acero inoxidable, pulido mecánicamente 0,0132 Agua 100,00 58,8
Agua-acero inoxidable, esmerilado y pulido 0,0080 Agua 150,00 48,2
Agua-acero inoxidable, picado y recubierto con teflón 0,0058 Agua 160,00 46,1
Agua-acero inoxidable esmerilado y pulido 0,0080 Agua 200,00 37,6
Agua-cobre 0,0130 Agua 226,70 31,9
Agua-cobre pulido con esmeril y tratado con parafina 0,0147 Agua 250,00 26,4
Agua-cobre pulido esmeril 0,0128 Agua 293,30 16,2
Agua-cobre grabado 0,0068 Agua 300,00 14,7
Agua-platino 0,0130 Agua 350,00 3,7
Benceno-cromo 0,0100 Agua 360,00 1,5
n-Pentano-níquel pulido con esmeril 0,0127 Agua 374,11 0,0
n-Pentano-cobre amolado 0,0049 Sodio 881,10 11,2
n-Pentano-cobre esmerilado 0,0074 Potasio 760,00 62,7
n-Pentano-cobre pulido esmeril 0,0154 Mercurio 357,20 39,4
Tetracloruro de carbono-cobre pulido esmeril 0,0070 Benceno 80,00 27,7
Tetracloruro de carbono-cobre 0,0130 Freón 11 44,40 8,5

A medida que el metal de la pared se calienta aún más, las burbujas se forman libremente y la eva-
poración es intensa, con gran turbulencia y altos valores del coeficiente hC de transferencia de calor. La
[Link] Condensación y vaporizació[Link].-284
formulación que se propone se puede aplicar a la evaporación de líquidos de un solo componente, sobre
superficies limpias.

Estudiando el significado de algunos parámetros relacionados con la convección forzada, Rohsenow

propuso la siguiente relación empírica, mediante la cual se puede estudiar el flujo de calor en el régimen
de evaporación nucleada:

Q g ( ρ l - ρ v ) c pl ( TpF - Ts ) 3
= ηl rl-v { }
A σ* rl-v Prl1,7 C

en la que:
cpl es el calor específico del líquido a la temperatura de saturación, y viene dado en [Link]/kg°C
Q/A es el flujo de calor en W/m2
rl-v es el calor latente de vaporización en [Link]/kg
Prl es el número de Prandtl correspondiente a la fase líquida, a la temperatura de saturación
C es una constante que se determina de acuerdo con datos experimentales, Tabla XIV.2
ΔT= TpF - Ts es la diferencia de temperaturas entre las correspondientes a la pared y al fluido, en con-
diciones de saturación
ηl es la viscosidad dinámica del líquido, a la temperatura de saturación, en kg/[Link]
ρ l y ρv son las densidades correspondientes a la fase líquida y de vapor saturado seco
σ* es la tensión superficial correspondiente a la superficie de separación líquido-vapor en N/m

Se ha encontrado, para superficies sucias o contaminadas, que el exponente del número de Prandtl
varía entre 0,8 y 2 en vez de ser 1,7 como aparece en la citada ecuación.

El coeficiente C es el único parámetro de la ecuación que permite su ajuste teniendo en cuenta la in-
fluencia que ejercen sobre el fenómeno, tanto el líquido a vaporizar, como la superficie calefactora; sus
valores para algunas combinaciones (pared-fluido) vienen resumidos en la Tabla XIV.2.

Pico de calor.- Al final de la etapa de evaporación nucleada, en el fondo del recipiente se ha alcan-
zado el momento de la máxima transferencia de calor por cuanto hasta aquí el líquido mojaba la superfi-
cie del recipiente, y lo refrigeraba, pero a partir de este instante, la refrigeración comienza a hacerse a
través de la película de vapor y esta nueva situación implica un cambio brusco en el coeficiente de con-
vección, por cuanto el del vapor es mucho menor que el del líquido y, por lo tanto, esta variación se tra-
duce en un aumento brusco de la temperatura de la pared calentada, por cuanto se sigue aplicando ca-
lor; este aumento de temperatura de la pared se conoce como pico de calor, cuya importancia y riesgo
radican en que existe la posibilidad de que el flujo de calor aplicado sobrepase un cierto valor, (al pasar
del régimen nucleado al de evaporación en película inestable, comienza a formarse la película de vapor,
(inestable), que separa el líquido de la pared, aparece el pico de calor, y como la superficie calefactora está
mal refrigerada por el vapor, se puede llegar a quemar).

Para estudiar este fenómeno se han desarrollado ecuaciones, como la de Zuber, que determina el
pico de calor Qmáx en la evaporación nucleada, mediante la expresión:

Qmáx σ * g ( ρl - ρ v ) ρl + ρv
= π ρv rl-v 4
A 24 ρ v2 ρl
[Link] Condensación y vaporizació[Link].-285
según la cual, si el pico de flujo de calor fuese alto, sería de desear que los valores de rl-v, ρv, g y σ* fuesen
también altos; el valor de rl-v, para el caso particular del agua es elevado y, en consecuencia, el pico de
calor que se obtiene en la evaporación del agua también es elevado.

Esta ecuación indica también la influencia que ejerce la gravedad en la vaporización y en el pico de
flujo térmico, por lo que un campo gravitacional reducido tiende a disminuir el pico de flujo térmico.

Evaporación en película (Régimen IV, Fig XIV.2).- Cuando se alcanza el pico de flujo de calor, ter-
mina la región de evaporación nucleada, y empieza la región de evaporación en película inestable.

No existe ningún análisis que permita calcular el flujo de calor en función de la diferencia de tempera-
turas (TpF - Ts) en la región inestable, ya que la pared está refrigerada alternativamente por vapor y por
líquido, sin ningún orden, hasta que se alcanza el punto mínimo de la curva de evaporación, en donde co-
mienza la región de evaporación en película estable, y la refrigeración de la pared, se sabe, se realiza
únicamente por el vapor.

En las regiones IV y V que se muestran en la Fig XIV.2, de evaporación en película inestable y en

película estable, la superficie de calefacción está separada del líquido por una capa de vapor a través de
la cual se tiene que transmitir el calor, necesitándose para ello grandes gradientes de temperatura,
puesto que la conductividad térmica del vapor es baja; por lo tanto, al no estar la pared suficientemente
refrigerada por el vapor, y continuar el aporte de energía calorífica a la misma, hace que el gradiente de
temperatura para el vapor aumente mucho, por lo que, cuando intervienen temperaturas elevadas, se
procura evitar la transferencia de calor en esta región, ya que por encima de un determinado gradiente
de temperaturas entre el líquido y la superficie calefactora, éste es incapaz de humedecerla, y es enton-
ces cuando se presenta una considerable reducción en el valor del coeficiente hC de transmisión de calor,
por cuanto sólo existe vapor en contacto con la pared.

Si en ese momento no se interrumpe el aporte de energía a la superficie metálica, la temperatura de

la pared aumentará para así acomodarse a la menor capacidad de la superficie para transferir calor,
hasta que la radiación emitida por la misma, más la transferencia de calor por evaporación pelicular,
sea igual a la energía incidente.

A partir de aquí, y por obtenerse temperaturas muy elevadas en la pared calefactora, se puede pro-
ducir la inutilización de la misma.

Para placa horizontal, y vaporización del pentano, tetracloruro de carbono, benceno, alcohol etílico,
etc., se puede considerar la siguiente ecuación (Berenson):

kv3 ρ v ( ρ l - ρ v ) g rl-v 0,4 c p(sat) ΔT

hC = 0,425 (1 + )
4 σ* rl-v
η v ΔT
g ( ρ l - ρv )

Flujo de calor mínimo (Régimen V, Fig XIV.2).- Zuber y Tribus proponen la siguiente ecuación
para el flujo de calor en la zona inestable, correspondiente al mínimo de la curva de vaporización:

Q ( ρ l - ρv ) g σ*
〉 = 0,09 ρ v rl-v 4
A mín ρ l + ρv g ( ρ l - ρv )

[Link] Condensación y vaporizació[Link].-286

 EVAPORACIÓN EN LA SUPERFICIE EXTERIOR DE UN HILO HORIZONTAL CALIENTE
O EN UN TUBO HORIZONTAL, POR CUYO INTERIOR CIRCULA UN FLUIDO CALEFACTOR
SUMERGIDO EN UN LÍQUIDO.

Cuando un líquido en reposo se calienta mediante un hilo por el que pasa una corriente eléctrica, o
mediante un tubo por cuyo interior circula un fluido calefactor, se puede llegar a un proceso de vaporiza-
ción similar al descrito anteriormente, con la diferencia de que el hilo o el tubo (pared calefactora), se en-
cuentra rodeado completamente por el líquido, estando así refrigerado por éste.

Si se llega a la temperatura de saturación del líquido, éste comienza a vaporizar en etapas idénticas
a las ya vistas en recipientes con líquidos en reposo, con la diferencia de que la evaporación en película
inestable y en película estable se realiza a través de una corona de vapor; cuando se empieza a formar
esta corona cambia bruscamente el coeficiente de convección y como ya se ha alcanzado el máximo hC
cuando refrigeraba el líquido, resulta que la temperatura del hilo aumenta bruscamente por una defi-
ciente refrigeración, a través de la película de vapor, hecho que puede detectarse porque la energía co-
municada al hilo experimenta una alteración.

El coeficiente de transmisión de calor hC para este tipo de evaporación en película estable sobre la
superficie exterior de un cilindro horizontal en ausencia de radiación, viene dado por la ecuación de
Bromley de la forma:

kv3 ρv ( ρ l - ρv ) g rl-v 0 ,4 c p( sat ) ΔT

hC = 0,62 4 (1 + )
η v de ΔT rl-v

ecuación que es válida para cuando la transferencia de calor a través de la película de vapor se realiza
por conducción pura, no interviniendo para nada la radiación.

Cuando la temperatura de la pared sea suficientemente alta, los efectos de la radiación son impor-
tantes, y el coeficiente promedio de transferencia de calor se calcula teniendo en cuenta los fenómenos
de radiación y de convección. Para ello se define un nuevo coeficiente de transmisión de calor por convec-
ción de la forma:

hC
hC* = hC 3 + hr
hC*

Fig XIV.3.- Vaporización de agua a la presión atmosférica y temperatura Ts calentada por un hilo de platino

[Link] Condensación y vaporizació[Link].-287

 que es difícil de utilizar, ya que hC * está en forma implícita, por lo que se puede hacer uso de otras ecua-
ciones, válidas para: u F < gd :

h hr
hC* = hC + hr { 3 + 1 r ( 1 )} , con: 0 < < 10 ; ± 0 ,3% de error
4 4 hC h hC
2,62 + r
hC
4 - T4 )
σ (T pF
3 hr hr s
hC* = hC + , con : < 1 ; ± 5% de error , y: hr =
4 hC 1
( + 1 - 1 ) ( TpF - Ts )
ε α
siendo ε la emisividad del tubo, α la absortividad del líquido y σ la constante de Stefan-Boltzman.
Cuando se cumpla que:

hr << hC ⇒ hC* = hC + 3 hr
4
Estas consideraciones son válidas para vaporizaciones a la presión atmosférica y una diferencia de
temperaturas, TpF - Ts < 1400°C.

Si: u F > 2 g d , se utiliza :

kv ρv uF rl-v 0,4 c p( sat ) ΔT

h C = 2,7 (1 + )
€ d e ΔT rl-v

ecuación que ha sido comprobada, en régimen laminar, para el benceno, tetracloruro de carbono, etanol,
etc., con diámetros comprendidos entre 10 y 16 cm y velocidades entre 0 y 4,5 m/seg.
€
El pico de calor viene dado por:

Qmáx σ * g ( ρ l - ρv ) ρl
〉 = 0 ,18 ρ v rl-v 4
A má x ρ v2 ρ l + ρv

En la Fig XIV.3 se han tenido en cuenta las siguientes consideraciones:

- Región I, de convección natural, en la que el calor se transmite al líquido a calentar, el cual se va ele-
vando hasta la superficie de separación líquido-vapor

- Región II, de evaporación nucleada, en la que las burbujas condensan en el líquido caliente

- Región III, también de evaporación nucleada, en la que las burbujas ascienden hasta la superficie de
separación líquido-vapor

- Región IV, en la que se tiene una evaporación nucleada parcial y evaporación en película nucleada
inestable, comienza en el máximo de la curva y acaba en el mínimo)

- Región V, de evaporación en película estable, que se corresponde con el comienzo de la zona de vapor
recalentado

- Región VI, en la que empieza a intervenir la radiación, junto con la convección, siendo ambas signifi-
cativas.

[Link] Condensación y vaporizació[Link].-288

XIV.9- EVAPORACIÓN DE LÍQUIDOS EN FLUJO FORZADO EN EL INTERIOR DE TUBOS

Fenomenología de la evaporación en tubos horizontales.- En la Fig XIV.4 se muestran los mo-

delos de flujos posibles en un flujo horizontal o inclinado con paredes isotermas.

Fig XIV.4.- Modelos de flujo horizontal de dos fases

Fig XIV.5.- Diagrama de flujo adiabático de dos fases para flujo horizonta

- El flujo en burbujas consiste en que la fase de vapor está en forma de burbujas aisladas inmersas en
la fase líquida, pudiendo ser pequeñas y esféricas, o grandes en forma de casquete esférico.

- En el flujo taponado las burbujas tienen un diámetro apreciable y tienden a fluir por la parte supe-
rior de la tubería.

- El flujo estratificado se presenta cuando la velocidad del líquido es pequeña, de forma que éste circu-
la por la parte inferior del tubo siendo su superficie relativamente uniforme.

- Si la velocidad del vapor aumenta se forman ondas de superficie líquida, que pueden llegar a ser lo
suficientemente grandes como para formar grandes masas de líquido que humedecen la parte superior de
la tubería, dando lugar al flujo en bala o en resalte.

- Si la velocidad del vapor es aún mayor, se provoca un flujo anular o un flujo en el que el líquido se
dispersa en forma de gotitas arrastradas por el flujo de vapor.

En la Fig XIV.5 se presenta un diagrama de flujo adiabático de dos fases para flujo horizontal, (Tai-
tel y Dukler), en el que se muestran las condiciones para una serie de regímenes de flujo. Si se añade
calor, los regímenes pueden desplazarse apreciablemente, por lo que dejarían de ser fiables.

[Link] Condensación y vaporizació[Link].-289

 Fenomenología de la evaporación en tubos verticales.- Cuando la evaporación acontece en el
interior de un tubo vertical caliente, a través del cual el fluido fluye, el proceso se realiza por convección
forzada. El mecanismo, y la hidrodinámica de la evaporación, en este caso, son mucho más complejos
que los de evaporación en recipientes con líquidos en reposo, como los vistos anteriormente, por cuanto
la velocidad del fluido dentro del tubo afecta al crecimiento de las burbujas y a su separación de la super-
ficie.

Fig XIV.6.- Modelos de flujo vertical de dos fases

Todavía no existe una teoría fundamentada que permita calcular el coeficiente de transferencia de
calor para este caso; mediante técnicas fotográficas se han podido seguir las complejas trayectorias del
flujo en los diferentes regímenes de la evaporación y, mediante ellas, se ha demostrado que el comporta-
miento del flujo en la evaporación por convección forzada es apreciablemente diferente del comporta-
miento de la evaporación en recipientes, a medida que aumenta la velocidad del fluido y el título del va-
por.

En la Fig XIV.6 se muestran, cuando el tubo se calienta uniformemente, las diferentes característi-
cas del flujo y en la Fig XIV.7 su efecto sobre el coeficiente de transmisión de calor; a medida que el líqui-
do avanza se va calentando, luego ebulle y, finalmente, el vapor producido se recalienta.

En la región de entrada se transfiere al líquido subenfriado calor por convección forzada, y este régi-
men permanece así hasta que empieza la evaporación, que viene acompañada por un aumento repenti-
no del coeficiente de transferencia de calor; en la región de evaporación aparecen burbujas sobre la su-
perficie caliente (flujo en burbujas), que crecen y son arrastradas hacia la corriente principal, prevale-
ciendo entonces en un cierto tramo del tubo un régimen de flujo con burbujas en el que la fase de vapor
está en forma de burbujas aisladas inmersas en la fase líquida, pudiendo ser pequeñas y esféricas, o
grandes en forma de casquete esférico.

A medida que las burbujas individuales van ascendiendo se agrupan para formar otras más grandes,
y aparecen en el flujo unas bolsas de vapor o tapones de diámetro algo menor que el del tubo, que frenan
el ascenso del líquido, dando lugar al llamado flujo taponado; esta situación prevalece hasta valores del
título de vapor x = 0,5. Es un tipo de flujo entre el de burbujas y el anular, que se identifica como flujo
agitado y flujo con arrastre de gotas.

El intervalo varía entre:

- Burbujas grandes, con el mismo diámetro que el tubo, separadas de la pared por una delgada pelí-
cula anular y separadas entre sí por el arrastre de agua líquida que pueden contener burbujas pequeñas

- Mezclas caóticas de grandes y pequeñas burbujas

Como la fracción del volumen de vapor aumenta por encima del valor anterior, la naturaleza del flui-
[Link] Condensación y vaporizació[Link].-290
do cambia notoriamente; el vapor empieza a fluir a través del centro del tubo como una corriente conti-
nua, en tanto que el líquido se adhiere a la pared y se desplaza formando una película anular; este tipo de
flujo recibe el nombre de régimen de película anular, en el que la película de fluido es cada vez más delga-
da a lo largo del tubo, pasando el título del vapor desde un valor x = 0,5 hasta x = 0,9 pudiéndose descom-
poner en otros dos:

- Flujo semianular, en el que el flujo es altamente irregular, a veces inestable, y consiste en un nú-
cleo de burbujas de gran tamaño que se unen y separan continuamente; el flujo de líquido tiende a estar
cerca de la pared.

- Flujo anular con niebla, en el que una película de líquido relativamente gruesa, que puede contener
pequeñas burbujas, fluye a lo largo de la pared; el núcleo de vapor contiene pequeñas gotitas de líquido
que forman nubes irregulares; al final x = 0,8 las fases están casi completamente separadas en un nú-
cleo de vapor y una película líquida sobre la pared, existiendo en el núcleo de vapor algunas gotas que
provienen de las crestas de las ondas que se forman en la superficie de la película.

Fig XIV.7.- Características del flujo durante la vaporización en tubo vertical cuando el tubo se calienta uniformemente

El coeficiente de transferencia de calor sigue siendo alto siempre que la película de fluido líquido hu-
medezca la pared calefactora del tubo; después, dependiendo de las condiciones de la superficie del tubo,
de la presión, y de la velocidad másica, aparecen en la pared del mismo unas zonas secas, en las que se
produce una brusca disminución del coeficiente de transmisión de calor; esta es la región de transición
del flujo anular en neblina, al flujo en neblina.

Las zonas secas continúan creciendo hasta que el líquido restante quede en forma de finísimas gotas
de agua; el flujo en neblina persistirá hasta que el título del vapor alcance la unidad, siendo el flujo a par-
[Link] Condensación y vaporizació[Link].-291
tir de este momento totalmente de vapor, comenzando el recalentamiento. Es a partir de la transición
flujo anular en neblina-flujo en neblina, cuando se modifica el tipo de refrigeración de la pared, (ya que
mientras estaba mojada por el líquido el coeficiente de convección aumentaba hasta un máximo), pa-
sando a otra refrigeración por vapor con un coeficiente de convección mucho menor, que es el que provo-
ca una mala transferencia térmica fluido-pared del tubo, lo que implica un aumento de la temperatura
de la pared, que es lo que ya hemos definido como pico de calor.

Fig XIV.8.- Efecto del aporte de calor en condiciones de flujo térmico crítico y
el incremento de temperatura correspondiente de la pared tubular para diversas condiciones de aporte de calor

En la Fig XIV.9 se presenta un diagrama de flujo adiabático de dos fases para flujo vertical, (Taitel y
Dukler), en el que se muestran las condiciones para
una serie de regímenes de flujo. Si se añade calor,
los regímenes pueden desplazarse apreciablemen-
te, por lo que dejarían de ser fiables.

Factores que afectan al flujo térmico crítico.-

La fenomenología asociada al flujo térmico crítico
en las condiciones de flujo que se presentan en los
generadores de vapor, fósiles y nucleares, está
afectada por una gran variedad de parámetros,
siendo los principales las condiciones de operación y
la geometría del diseño.

a) Las condiciones de operación afectan al flujo tér-

Fig XIV.9 mico crítico, son la presión, el flujo másico y el título
Diagrama de flujo vertical ascendente adiabático de dos fases
del vapor

b) Los factores geométricos de diseño incluyen:

* Las dimensiones y perfil de las secciones de paso del flujo másico
* Las obstrucciones en el recorrido del flujo másico
* El perfil del flujo térmico
* La inclinación y configuración de la superficie de paredes conductoras del flujo másico, etc.

Algunas de estas dependencias se representan en la Fig XIV.8, que muestra (para un tubo vertical
uniformemente calentado, y refrigerado interiormente por un flujo ascendente de agua) el efecto que ori-
[Link] Condensación y vaporizació[Link].-292
gina en la temperatura un incremento en el aporte de calor:

- Con flujos térmicos bajos, el flujo de agua se puede vaporizar casi por completo, antes de que se pue-
da observar una elevación de la temperatura

- Con flujos térmicos moderados y altos, la posición del flujo térmico crítico se mueve progresivamente
hacia la entrada del tubo al tiempo que se incrementa la temperatura

- Con flujos térmicos muy altos, el flujo térmico crítico se presenta con un bajo título del vapor, pu-
diendo ser la temperatura del metal lo suficientemente elevada como para que se llegue a fundir el tubo

- Con regímenes de aporte de calor extraordinariamente altos, el flujo térmico crítico se puede presen-
tar incluso en agua subenfriada

Gradiente de presiones en el flujo en 2 fases en el interior de tubos verticales.- La caída de

presión del flujo en dos fases dentro de los tubos es un proceso complicado dada la gran variedad de for-
mas de flujos existentes.

a) El gradiente de presión en el flujo de dos fases se puede calcular como sigue:

⎧ G viene dada en kg ; 2
G d ⎪ A = πd
seg 4
Re = A , en la que: ⎨ 1 x 1 - x
η ref
⎪ η ref = η v + η l
⎩

siendo (1/ηref) la inversa de la viscosidad dinámica de referencia que compensa las fluctuaciones del fac-
tor de fricción λ en el tubo, que es demasiado bajo si x < 0,7 y demasiado alto si x > 0,7

La densidad se define en la forma: ρ = x ρv + ( 1 - x) ρ l

El coeficiente de rozamiento en flujo turbulento es: λ = ( 0 ,79 ln Red − 1,64) -2

El gradiente de presión debido a la fricción en la pared es:

dp λ ( 4 G2 ) 2
( ) = - πd
dz pared 2ρd

El gradiente de presión debido a la gravedad es:

dp
( ) = - ρ g sen θ
dz grav

El gradiente de presión debido a la variación de la cantidad de movimiento es:

dp 4 G/π d 2 2 dp 4 G/π d 2 2
( )Δmov = ( ) =-( ) ρ l dx
dz ρ dz ρ dz

en la que el gradiente del título de vapor es :

dx = Q/L Q
, con dado en W
dz 4 G rl-v/π d 2 L m

El gradiente de presión total es:

[Link] Condensación y vaporizació[Link].-293

dp dp dp dp
=( ) +( ) +( )
dz dz pared dz gravedad dz Δmov

b) Otra formulación debida a Wallis, válida para el flujo con burbujas x ≤ 0,5, relaciona bastante
bien los datos experimentales de la caída de presión en dos fases, mediante la expresión:

Δ pTubo ρ G 3
= 1+ 3 l v 10-6 G
Δ plíq ρ v Gl

⎧ Δptubo es la caída de presión en el tubo, cuando coexisten las dos fases

en la que: ⎨ Δp
⎩ líquido es la caída de presión en el tubo correspondiente a la fase líquida

COEFICIENTE DE CONVECCIÓN EN LA EVAPORACIÓN EN EL INTERIOR DE TUBOS

VERTICALES PARA EL FLUJO EN 2 FASES.- Los coeficientes de transmisión de calor para líqui-
dos que se calientan, hasta que comienza la evaporación, y para el vapor recalentado se pueden deter-
minar mediante las correlaciones que definen la convección forzada dentro de tubos. Sin embargo, no
existe en la actualidad una formulación general con la exactitud necesaria, que permita determinar el
coeficiente de transmisión de calor de todos los regímenes de evaporación que pueden darse en el interior
del tubo.

- David y David han propuesto una relación empírica que permite determinar el coeficiente de
transmisión de calor promedio hC para un flujo de dos fases en el interior de un tubo, siempre que el líqui-
do humedezca la superficie interior del mismo, es decir, cuando se esté en el régimen de película anular x
≤ 0,9, en la forma:

hC di ρ d G* 0,87 0,4
= 0,06 ( l )0,28 ( i )x Prl
kl ρv η
en la que x es el título de vapor y G* = 3600 ρx ux en (kg/hm2).

Los datos relativos a las propiedades del fluido hay que tomarlos de las tablas y diagramas corres-
TpF - Ts
pondientes a la temperatura media de la capa de condensado .
2
- V. Klimenko propone una formulación válida sólo cuando las paredes del tubo estén mojadas, por
lo que lo primero que hay que realizar es comprobar esta situación, es decir, si domina la evaporación
nucleada (pared húmeda) o la evaporación en película, (pared seca) ; para ello se evalúa un parámetro Φ
que define el tipo de evaporación, función del flujo de calor por unidad de superficie de la pared interior del
tubo, de la forma:

4 G rl-v ρ ρ Q/L
Φ= { 1 + x ( l - 1 )} ( l )3 , siendo ; q =
πd q2 ρ v ρ v πd

Si, Φ < 1,6.104 (ebullición nucleada):

kpared 0 ,15 hC( 2 fases )

Nu = 7 ,4.10-3 q *0 ,6 p* 0 ,5 Prl-1/3 ( ) = LC
kl kl

siendo Lc (longitud característica), q* y p* de la forma:

σ q LC p p LC
LC = ; q*= ; p* = =
g ( ρl - ρv ) rl-v ρ v α l σ g ( ρ l - ρv ) σ
[Link] Condensación y vaporizació[Link].-294
 Si, Φ > 1,6.104 (ebullición en película):

ρ v 0 ,2 kpared 0 ,09 hC( 2 fases )

Nu = 0 ,087 Re 0 ,6 Prl1/6 ( ) ( ) = LC
ρl kl kl

v* LC ρ
Re = ; v* = G {1 + x ( l - 1)}
νl ρl ρv

evaluándose las propiedades a Tsat

Finalmente hay que calcular un coeficiente de transmisión de calor por convección de una sola fase
hC(1 fase) basado en un número de Reynolds de líquido puro:

(G/A ) d
Re =
ηl

Puede suceder que: hC(1 fase)<< hC(2 fases), por lo que no se consideraría, es decir:

hC = hC(2 fases)

Si hC(1 fase) no es pequeño, el coeficiente de convección

hC = 3 hC3 ( 1 fase )líq + hC3 ( 2 fases )

Estas correlaciones son válidas para el agua y refrigerantes con un error menor del 10%, en el campo:

- Presión: 0,61 atm ÷ 30,4 atm

- Flujo másico por unidad de superficie: 50 ÷ 2700 kg/m2seg
- Título de vapor: 0,017 ÷ 1
- Diámetro del tubo: 1,63 ÷ 41,3 mm

[Link] Condensación y vaporizació[Link].-295