Algoritmo de cálculo

1. Realizar los balances de energía tanto para condensación en gota como para condensación en
película, teniendo en cuenta que el proceso está en estado estable, por lo tanto:

𝑄̇𝑎𝑔𝑢𝑎 ≅ 𝑄̇𝑣𝑎𝑝𝑜𝑟

𝑇𝑠𝑎𝑙
1
𝑚̇𝑎𝑔𝑢𝑎 ∙ ∙∫ (𝑎 + 𝑏𝑇 + 𝑐𝑇 2 + 𝑑𝑇 3 )𝑑𝑇 = 𝑚̇𝑣𝑎𝑝𝑜𝑟 ∙ 𝜆𝑣
𝑃𝑚 𝐻 𝑇𝑒𝑛𝑡
2𝑂

Dónde:

𝜆𝑣 = 2275 kJ/kg a P = 585 mmHg (Presión atmosferica Manizales)

𝑎 = 1.82964 ∙ 101
𝑏 = 4.72118 ∙ 10−1
𝑐 = −1.33876 ∙ 10−3
𝑑 = 1.31424 ∙ 10−6
𝑃𝑚 𝐻 𝑂 = 18.01528 𝑔/𝑚𝑜𝑙
2

Y las temperaturas utilizadas en la integral son:

 Condensación en película

T1 = Temperatura de entrada del agua (𝐾)

T3 = Temperatura de salida del agua (𝐾)

 Condensación en gota

T4 = Temperatura de entrada del agua (𝐾)

T6 = Temperatura de salida del agua (𝐾)

Se calcula el error para asegurar que el proceso se encuentra en estado estable, se admite un
error del 5 %.

𝑄̇𝑣𝑎𝑝𝑜𝑟 − 𝑄̇𝑎𝑔𝑢𝑎
%𝑒𝑟𝑟𝑜𝑟 = ∙ 100
𝑄̇𝑣𝑎𝑝𝑜𝑟

2. Calcular los coeficientes peliculares de condensación

2.1. Para la condensación en gota

Para la condensación en gota se determina los coeficientes tanto experimental como teórico:

 Para la parte experimental se utiliza:

 𝑄̇𝑣𝑎𝑝𝑜𝑟
ℎ𝑜 =
𝐴0 ∙ ∆𝑇𝑣𝑎𝑝𝑜𝑟𝑔

Dónde:

𝐴0 = 𝜋 ∙ 𝐷𝑜 ∙ 𝐿
∆𝑇𝑣𝑎𝑝𝑜𝑟𝑔 = 𝑇𝑠𝑎𝑡 − 𝑇𝑝𝑝
𝑇7 + 𝑇8
𝑇𝑝𝑝 =
2
T𝑠𝑎𝑡 = Temperatura de saturación del agua a P = 585 mmHg (𝐾)
T𝑝𝑝 = Temperatura promedio de pared (𝐾)
T7 = Temperatura del equipo a la izquierda del tubo de condensación en gota (𝐾)
T8 = Temperatura del equipo a la derecha del tubo de condensación en pelicula (𝐾)

 Para la parte teórica se utiliza:

ℎ𝑑𝑐 = 51.104 + 2044𝑇𝑠𝑎𝑡 (°C); 22°C < 𝑇𝑠𝑎𝑡 < 100°C

ℎ𝑑𝑐 = 255.510; 𝑇𝑠𝑎𝑡 > 100°C

2.2. Para la condensación en película

Para la condensación en película se calcula el coeficiente experimental, suponiendo que la

transferencia de calor se da por convección a través del condensado y que la capa de líquido es
muy delgada:

𝑄̇𝑣𝑎𝑝𝑜𝑟
ℎ𝑜 =
𝐴0 ∙ ∆𝑇𝑣𝑎𝑝𝑜𝑟𝑝

Dónde:

𝐴𝑡𝑟𝑎𝑛𝑠 ≈ 𝐴𝑒𝑥𝑡𝑒𝑟𝑖𝑜𝑟𝑡𝑢𝑏𝑜
𝐴0 = 𝜋 ∙ 𝐷𝑜 ∙ 𝐿
∆𝑇𝑣𝑎𝑝𝑜𝑟𝑝 = 𝑇𝑠𝑎𝑡 − 𝑇2
T2 = Temperatura de pared de la tuberia para la condensación en pelicula (𝐾)
3. Calcular el coeficiente pelicular convectivo, para el agua de enfriamiento

 Para la condensación en gota

- Flujo en el interior de los tubos
 Régimen laminar: Ecuación de Sieder y Tate
1
𝐷𝑖 3 𝜇 0.14
𝑁𝑢 = 1.86 (𝑅𝑒𝑃𝑟 ) ( )
𝐿 𝜇𝑠
𝑁𝑢 ∙ 𝑘
ℎ𝑖 =
𝐷𝑖
 Régimen turbulento: Ecuación de Sieder y Tate
1 𝜇 0.14
𝑁𝑢 = 0.027𝑅𝑒 0.8 𝑃𝑟 3 ( )
𝜇𝑠
Dónde:
𝐺𝐷𝑖 𝑚̇𝑎𝑔𝑢𝑎
𝑅𝑒 = → 𝐺= → 𝑁𝑡 = 1
𝜇 𝐿 ∙ 𝑁𝑡
𝐶𝑝 ∙ 𝜇
𝑃𝑟 =
𝑘
𝐷𝑖 = Diámetro interno de los tubos (m)
𝐿 = Longitud de los tubos(m)
𝜇, 𝐶𝑝, 𝑘 = Propiedades del agua a la temperatura media global (Tmg )
𝑇4 + 𝑇6
𝑇𝑚𝑔 =
2
𝜇𝑠 = Viscosidad a T5
T5 = Temperatura de pared de la tuberia para la condensación en gota (𝐾)

4. Calcular los coeficientes con las correlaciones matemáticas encontradas en la bibliografía.

 Para la condensación en película

- Por la teoría de Nusselt:
1
3 3
𝑘𝑝𝑒𝑙 ∙ 𝜌𝑝𝑒𝑙 ∙ 𝜆𝑣 ∙ 𝑔 4
ℎ̅ = 0.725 ∙ ( )
𝜇𝑝𝑒𝑙 ∙ 𝐷𝑜 ∙ ∆𝑇𝑓
Dónde:
𝜇, 𝜌, 𝑘 = Propiedades del agua a la temperatura de película(T𝑝𝑒𝑙 ), u. inglesas
𝑇𝑠𝑎𝑡 + 𝑇𝑝𝑝
𝑇𝑝𝑒𝑙 =
2
 𝑇7 + 𝑇8
𝑇𝑝𝑝 =
2
𝑇𝑠𝑎𝑡 − 𝑇𝑝𝑝
∆𝑇𝑓 =
2
𝐷𝑜 = Diámetro externo de los tubos (m)
Para el régimen laminar se encuentra la correlación de McAdams para tubos
horizontales.
3100
ℎ̅ = 1 1 → 1800 < 𝑅𝑒 < 2100 → ∆𝑇𝑓 , 𝑟𝑎𝑛𝑔𝑜 𝑑𝑒 10 𝑎 150 °𝐹
𝐷𝑜 4 ∙ ∆𝑇𝑓 3

5. Calcular los coeficientes totales de transmisión de calor: UC , UD y el factor R d .

5.1. Para la condensación en gota

 UD :
𝑄̇𝑎𝑔𝑢𝑎
𝑈𝐷 =
𝐴𝑖 ∙ 𝑀𝐿𝐷𝑇𝐺
Dónde:
𝐴𝑖 = 𝜋 ∙ 𝐷𝑖 ∙ 𝐿
(𝑇𝑠𝑎𝑡 − 𝑇4 ) − (𝑇𝑠𝑎𝑡 − 𝑇6 )
𝑀𝐿𝐷𝑇𝐺 =
𝑇 −𝑇
ln (𝑇𝑠𝑎𝑡 − 𝑇4 )
𝑠𝑎𝑡 6

 UC :
ℎ𝑖𝑜 ∙ ℎ𝑜
UC =
ℎ𝑖𝑜 + ℎ𝑜
Dónde:
𝐷𝑖
ℎ𝑖𝑜 = ℎ𝑖 ∙
𝐷𝑜
𝑘 1 𝜇 0.14
ℎ𝑖 = 0.027 ∙ ∙ 𝑅𝑒 0.8 𝑃𝑟 3 ( )
𝐷𝑖 𝜇𝑠
𝜇, 𝑘 = Propiedades del agua a la temperatura media global (Tmg )
𝜇𝑠 = Viscosidad a T5
  Rd:
U C − UD
Rd =
UC ∙ UD

5.2. Para la condensación en película

 UD :
𝑄̇𝑎𝑔𝑢𝑎
𝑈𝐷 =
𝐴𝑖 ∙ 𝑀𝐿𝐷𝑇𝑃
Dónde:
𝐴𝑖 = 𝜋 ∙ 𝐷𝑖 ∙ 𝐿
(𝑇𝑠𝑎𝑡 − 𝑇1 ) − (𝑇𝑠𝑎𝑡 − 𝑇3 )
𝑀𝐿𝐷𝑇𝑃 =
𝑇 −𝑇
ln (𝑇𝑠𝑎𝑡 − 𝑇1 )
𝑠𝑎𝑡 3

 UC :
ℎ𝑖𝑜 ∙ ℎ𝑜
UC =
ℎ𝑖𝑜 + ℎ𝑜
Dónde:
𝐷𝑖
ℎ𝑖𝑜 = ℎ𝑖 ∙
𝐷𝑜
𝑘 1 𝜇 0.14
ℎ𝑖 = 0.027 ∙ ∙ 𝑅𝑒 0.8 𝑃𝑟 3 ( )
𝐷𝑖 𝜇𝑠
𝜇, 𝑘 = Propiedades del agua a la temperatura media global (Tmgp )
𝑇1 + 𝑇3
𝑇𝑚𝑔 =
2
𝜇𝑠 = Viscosidad a T2
1
𝑘 3 ∙ 𝜌 3 ∙ 𝜆𝑣 ∙ 𝑔 4
ℎ𝑜 = 0.725 ∙ ( )
𝜇 ∙ 𝐷𝑜 ∙ ∆𝑇𝑝
∆𝑇𝑝 = 𝑇𝑠𝑎𝑡 − 𝑇2
𝜇, 𝜌, 𝑘 = Propiedades del agua a𝑇2

 Rd:
U C − UD
Rd =
UC ∙ UD
6. Comprobar las ecuaciones de Nusselt en tubo horizontal

- Para la condensación en película

Para comprobar las ecuaciones de Nusselt, se debe cumplir ciertas condiciones como la uniformidad
de la temperatura y de régimen que se trabaje.
4 ∙ 𝐺′
𝑅𝑒 =
𝜇
𝑚̇𝑎𝑔𝑢𝑎
𝐺′ =
𝐿
 Para el cálculo del coeficiente convectivo externo, se utiliza:
1
3 3
𝑘𝑝𝑒𝑙 ∙ 𝜌𝑝𝑒𝑙 ∙ 𝜆𝑣 ∙ 𝑔 4
ℎ̅ = 0.725 ∙ ( )
𝜇𝑝𝑒𝑙 ∙ 𝐷𝑜 ∙ ∆𝑇𝑓

7. Caída de presión para la condensación en gota y en película:

 Interior de los tubos: Realizando un balance de energía y utilizando la correlación de

Bernoulli, con cambios de densidad despreciables:

2 2
𝑃1 𝑉̅1 𝑔𝑧1 𝑃2 ̅̅̅
𝑉2 𝑔𝑧2
+ + = + + + ℎ𝑓𝑠
𝜌 2𝑔𝑐 𝑔𝑐 𝜌 2𝑔𝑐 𝑔𝑐
Y asumiendo que:
𝑧1 = 𝑧2
𝑉̅1 = ̅̅̅
𝑉2 = 𝑐𝑡𝑒
1
(𝑃 − 𝑃2 ) = ℎ𝑓𝑠
𝜌 1
Dónde:
𝐿 𝑉2
ℎ𝑓𝑠 = 4 ∙ 𝑓 ∙ ∙( )∙𝜌
𝐷𝑖 2𝑔𝑐
Luego se tiene:
𝐿 𝑉2
∆𝑃 = 4 ∙ 𝑓 ∙ ∙ ( )
𝐷𝑖 2𝑔𝑐
Dónde:
𝑚̇𝑎𝑔𝑢𝑎 1
𝑉= ∙
𝜌 𝐴𝑖
ρ = densidad del agua a Tmg ó Tmgp
 El factor de fricción de fanning (f) se encuentra con correlaciones dependiendo del
régimen laminar o turbulento, o también se puede utilizar un diagrama como Moody
dependiendo de la superficie.

 Fuera de los tubos: El ∆Pvapor, es la medida registrada por el manómetro de nivel, la

caída de presión se da como:
𝑃𝑣𝑎𝑝𝑜𝑟 = 𝑃𝑎𝑡𝑚 + ∆Pvapor
∆Pvapor ≪ 𝑃𝑎𝑡𝑚
𝑃𝑣𝑎𝑝𝑜𝑟 ≅ 𝑃𝑎𝑡𝑚

8. Calcular Re para la película y graficar Nu vs Re.

ℎ𝑖 ∙ 𝐷𝑖
𝑁𝑢 =
𝑘
Los valores de k se calculan a la 𝑇𝑝𝑒𝑙 , y con los valores de h encontrados en la parte No. 2, se halla
Nusselt y se grafican contra los números de Reynolds hallados por la teoría de Nusselt, y se observa
el comportamiento de estos.

9. Implementar una ecuación de tipo experimental para el cálculo del coeficiente de condensación
pelicular ℎ𝑐 y comparar con lo reportado en la bibliografía.

Para efectuar la realización de la ecuación experimental se deben tener en cuenta las variables que
se ven involucradas en la determinación del coeficiente pelicular como lo son la Temperatura, la
conductividad térmica, la viscosidad, el diámetro, la longitud, la velocidad de flujo y la capacidad
calorífica. Luego se aplica un método que relacione todas las variables en este caso el Teorema Pi
de Buckingham y con los datos obtenidos en el laboratorio y el método de mínimos cuadrados se
pueden encontrar los exponentes involucrados en la ecuación.

ℎ = 𝑓(𝑇, 𝑘, 𝐷, 𝐿, 𝑉, 𝐶𝑝, 𝜇)

Las dimensiones de todas las variables involucradas se presentan así:

 Temperatura (𝑇) → 𝑇
𝑀
 Viscosidad (𝜇) → 𝐿∙𝑡
 𝐸
 Conductividad térmica (𝑘) → 𝑡∙𝐿∙𝑇

 Diámetro (𝐷) → 𝐿
 Longitud (𝐿) → 𝐿
𝐿
 Velocidad de flujo (𝑉) → 𝑡
𝐸
 Capacidad calorífica 𝐶𝑝 →
𝑀∙𝑇

Dónde:
 T : temperatura
 t : tiempo
 M : masa
 E : energía
 L : longitud

Y las dimensiones del coeficiente pelicular son:

𝐸
ℎ → [𝐿2 ∙𝑡∙𝑇]

Teniendo en cuenta las dimensiones de las variables, se plantea la siguiente relación dimensional:

[ℎ] = [𝑇]𝑎 [𝜇]𝑏 [𝑘]𝑐 [𝐷]𝑑 [𝐿]𝑒 [𝑉] 𝑓 [𝐶𝑝] 𝑔

𝐸 𝑎
𝑀 𝑏 𝐸 𝑐
𝑑 𝑒
𝐿 𝑓 𝐸 𝑔
= 𝑇 ∙ ( ) ∙ ( ) 𝐿 ∙ 𝐿 ∙ ( ) ∙ ( )
𝐿2 ∙ 𝑡 ∙ 𝑇 𝐿∙𝑡 𝑡∙𝐿∙𝑇 𝑡 𝑀∙𝑇

Se supone esta ecuación como dimensionalmente homogénea, donde los exponentes del lado
derecho tienen que ser iguales a las del lado izquierdo.

 Exponentes de E : 1 = 𝑐 + 𝑔
 Exponentes de L : −2 = −𝑏 − 𝑐 + 𝑑 + 𝑒 + 𝑓
 Exponentes de t : −1 = −𝑏 − 𝑐 − 𝑓
 Exponentes de T : −1 = 𝑎 − 𝑐 − 𝑔
 Exponentes de M : 0 = 𝑏 − 𝑔

De este sistema se obtiene:

  𝑎=0
 𝑏=𝑔
 𝑐 =1−𝑔
 𝑑 = −1 − 𝑒
 𝑓=0

Por lo tanto:

ℎ = (𝑇)0 ∙ (𝜇)𝑏 ∙ (𝑘)1−𝑏 ∙ (𝐷)−1−𝑒 ∙ (𝐿)𝑒 ∙ (𝑉)0 ∙ (𝐶𝑝)𝑏

𝐶𝑝 ∙ 𝜇 𝑏 𝐿 𝑒 𝑘
ℎ=( ) ∙( ) ( )
𝑘 𝐷 𝐷
𝐿 𝑒 𝑘
ℎ = (𝑃𝑟)𝑏 ∙ ( ) ( )
𝐷 𝐷

Los valores de los exponentes b y e se encuentran correlacionando Nusselt y Prant de todos los
ensayos realizados durante la práctica.

Bibliografía
Espinoza, A. A. (s.f.). Determinación experimental de h a partir de los números de Nusselt, Reynolds
y Prandtl. Recuperado el 04 de Mayo de 2019, de
[Link]
partir-de-los-nmeros-de-nusselt-reynolds-y-prandtl

Fabio Augusto Mesa Rueda, G. A. (2006). Manuel de prácticas para el laboratorio de operaciones
unitarias II. Manizales.

Universidad Nacional de Colombia. (s.f.). Condensación en gota y en película. Recuperado el 04 de

Mayo de 2019, de [Link]
colombia/cinetica-y-catalisis/ejercicios-obligatorios/condensacion-en-gota-y-en-pelicula-
laboratorio-operaciones-unitarias/2653070/view

Universidad Nacional Experimental "Rafael María Baralt". (Septiembre de 2015). Transferencia de

calor en la condensación. Recuperado el 04 de Mayo de 2019, de
[Link]
Universidad Tecnologica Nacional. (Mayo de 2009). Fenómenos de transporte, Notas de cátedra.
Recuperado el 04 de Mayo de 2019, de
[Link]
/Unidad_tematica_6.pdf

Yaws, C. L. (s.f.). Propiedades de gases y liquidos.