Ampliación de Energía y Máquinas Térmicas. Formulario Examen AEMT - 1º Parte. 2016-2017.

ECUACIONES DE CONSERVACIÓN ESTACIONARIAS

Masa:
𝛼𝜌𝑣 𝑣𝑣
=Γ
𝑧
(1 − 𝛼)𝜌𝑙 𝑣𝑙
= −Γ
𝑧
Momento:
𝐴𝑧 𝐹𝑤𝑣 𝐴𝑧 𝜏𝑠
𝛼𝜌𝑣 𝑣𝑣2= Γ𝑣𝑠 𝐴𝑧 − − 𝛼𝑝 − − 𝛼𝜌𝑣 𝑔𝐴𝑧 𝑐𝑜𝑠𝜃
𝑧 𝑧 𝑧 𝐴𝑠
𝐴
2 𝑧
𝐹𝑤𝑙 𝐴𝑧 𝜏𝑠
(1 − 𝛼)𝜌𝑙 𝑣𝑙 = −Γ𝑣𝑠 𝐴𝑧 − − (1 − 𝛼)𝑝 + − (1 − 𝛼)𝜌𝑙 𝑔𝐴𝑧 𝑐𝑜𝑠𝜃
𝑧 𝑧 𝑧 𝐴𝑠
Energía:
𝐴𝑧 𝐴𝑧 𝐴𝑠 𝐴𝑠
𝛼𝜌𝑣 𝑣𝑣 ℎ𝑣0 = −[𝑞𝑤𝑣′′ 𝛼]
𝑧 − 𝛼𝜌𝑣 𝑔𝑣𝑣 𝐴𝑧 𝑐𝑜𝑠𝜃 + (Γℎ𝑣0 )𝑠 𝐴𝑧 − 𝑞𝑣′′ − 𝑝𝑣𝑠
𝑧 𝑧 𝑧 𝑧
𝐴𝑧 𝐴𝑧 𝐴
′′ 𝑠
𝐴𝑠
(1 − 𝛼)𝜌𝑙 𝑣𝑙 ℎ𝑙0 ′′ (1
= −[𝑞𝑤𝑙 − 𝛼)]𝑧 (1
− − 𝛼)𝜌𝑙 𝑔𝑣𝑙 𝐴𝑧 𝑐𝑜𝑠𝜃 − (Γℎ 0)
𝑣 𝑠 𝐴𝑧 + 𝑞𝑙 + 𝑝𝑣𝑠
𝑧 𝑧 𝑧 𝑧
Velocidad interfacial:
𝑣⃗𝑠 = 𝜂𝑣⃗𝑣 + (1 − 𝜂)𝑣⃗𝑙 , 𝜂 = 0 𝑠𝑖 Γ < 0
𝜂 = 1 𝑠𝑖 Γ > 0
Caída de presión por rozamiento:
𝐹𝑤𝑔 𝐹𝑤𝑙 𝑣2
∆𝑝 = − − − 𝑧𝜌𝑚 𝑔𝑐𝑜𝑠𝜃 − 𝑘𝜌𝑚
𝐴𝑧 𝐴𝑧 2

REGÍMENES DE CAUDALES EN TUBERÍAS VERTICALES (tubería entre 46º y 90º)

Velocidad de ascensión de las burbujas pequeñas:
1⁄4 ⁄
𝑔(𝜌𝑙 − 𝜌𝑣 )𝜎 𝐷∗ 𝜎 1 2
𝑣𝑠𝑏 = 1.53 [ 2 ] = 1.53 [ ]
𝜌𝑙 𝐷𝜌𝑙
Velocidad de ascensión de la burbuja Taylor:
⁄
𝑔𝐷(𝜌𝑙 − 𝜌𝑣 ) 1 2 𝜎 1⁄2
𝑣𝑇𝑏 = 0.35 [ ] = 0.35𝐷∗ [ ]
𝜌𝑙 𝐷𝜌𝑙
∗
Diámetro adimensional tubo que permite la presencia del flujo bubbly is 𝐷 ≥ 22.22,
1⁄2
𝑔(𝜌𝑙 − 𝜌𝑣 )
𝐷∗ = 𝐷 [ ]
𝜎
Criterio transición estratificación vertical.
1
𝑣𝑚 < 𝑣𝑇𝑏 , Estratificado
2
𝑣𝑚 > 𝑣𝑇𝑏 , No Estratificado
Criterio transición Bubbly a Slug.
- Para flujos másicos de 𝐺𝑚 < 2000 𝑘𝑔/𝑚2 𝑠: 𝛼𝐵𝑆 = 𝛼𝐿 .
- Para flujos másicos de 2000 < 𝐺𝑚 < 3000 𝑘𝑔/𝑚2 𝑠:
𝛼𝐵𝑆 = 𝛼𝐿 + 0.001(0.5 − 𝛼𝐿 )(𝐺𝑚 − 2000).
- Para flujos másicos de 𝐺𝑚 ≥ 3000 𝑘𝑔/𝑚2 𝑠: 𝛼𝐵𝑆 = 0.5.
𝐷∗ 8
Con: 𝛼𝐿 = 0.25 𝑚𝑖𝑛 [1.0, ( ) ]
22.22
Transición Slug a Slug-Annular.
𝑓 𝑒
Entre 𝛼𝐷𝐸 y 𝛼𝑆𝐴 : 𝛼𝐷𝐸 = max(𝛼𝐵𝑆 , 𝛼𝑆𝐴 − 0.05) y 𝛼𝑆𝐴 = max{0.5, min(𝛼𝑐𝑟𝑖𝑡 , 𝛼𝑐𝑟𝑖𝑡 , 0.9)}.
𝑓 1 𝑔𝐷(𝜌𝑙 −𝜌𝑣 ) 1⁄2
𝛼𝑐𝑟𝑖𝑡 = [ ] , caudal hacia arriba.
𝑣𝑣 𝜌𝑣
𝑓
𝛼𝑐𝑟𝑖𝑡 = 0.75, caudal hacia abajo o en contracorriente.
1⁄4
𝑒
3.2 𝑔𝜎(𝜌𝑙 − 𝜌𝑣 )
𝛼𝑐𝑟𝑖𝑡 = [ ]
𝑣𝑣 𝜌𝑣2
Transición Slug-Annular a Annular-Mist.
𝑓𝑒
𝛼𝑆𝐴 = max{0.5, min(𝛼𝑐𝑟𝑖𝑡 , 𝛼𝑐𝑟𝑖𝑡 , 0.9)}
Transición Annular Mist y pre-CHF Mist (MPR).
𝛼𝐴𝑀 = 0.9999
Post-CHF.
Si el volumen hidrodinámico tiene flujo calorífico al vapor debido a una pared circundante y la temperatura del vapor está
más de 1K sobrecalentado, entonces el régimen de caudal es post-CHF. Para transferencia de calor post-CHF, las
formulaciones anteriores también se usan para definir los regímenes equivalentes, menos para el siguiente régimen:
Transición Slug invertido y Annular invertido.
Entre 𝛼𝐵𝑆 y 𝛼𝐶𝐷 , con 𝛼𝐶𝐷 = 𝛼𝐵𝑆 + 0.2

1
Ampliación de Energía y Máquinas Térmicas. Formulario Examen AEMT - 1º Parte. 2016-2017.

REGÍMENES DE CAUDALES EN TUBERÍAS HORIZONTALES (tubería entre 0º y 45º)

Velocidad crítica:
1 (𝜌𝑙 − 𝜌𝑣 ) ∙ 𝑔 ∙ 𝛼𝑣 ∙ 𝐴
𝑣𝑐𝑟𝑖𝑡 = ∙√ ∙ (1 − 𝑐𝑜𝑠(𝜃) )
2 𝜌𝑣 ∙ 𝐷ℎ ∙ 𝑠𝑒𝑛(𝜃)
𝜽 es el ángulo entre la vertical y el nivel del líquido estratificado:
𝛼𝑣 ∙ 𝜋 = 𝜃 − 𝑠𝑖𝑛𝜃𝑐𝑜𝑠𝜃
Criterio transición estratificación horizontal.
- Si |𝑣𝑣 -𝑣𝑙 | > 𝑣𝑐𝑟𝑖𝑡 y 𝐺𝑚 > 3000 𝑘𝑔/𝑚2 𝑠  el régimen es no estratificado.
- Si |𝑣𝑣 -𝑣𝑙 | < 𝑣𝑐𝑟𝑖𝑡 y 𝐺𝑚 < 3000 𝑘𝑔/𝑚2  el régimen es de transición a estratificado.
- Si |𝑣𝑣 -𝑣𝑙 | < 0.5𝑣𝑐𝑟𝑖𝑡 y 𝐺𝑚 < 2500 𝑘𝑔/𝑚2 𝑠 el régimen es completamente estratificado.
Criterio transición Bubbly a Slug.
𝑘𝑔
0.25 𝐺𝑚 ≤ 2000 2
𝑚 𝑠
𝑘𝑔
𝛼𝐵𝑆 = 0.25 + 0.00025 ∙ (𝐺𝑚 − 2000) 2000 ≤ 𝐺𝑚 ≤ 3000 2
𝑚 𝑠
𝑘𝑔
{ 0.5 𝐺𝑚 > 3000 2
𝑚 𝑠
Transición Slug a Slug-Annular
𝛼𝐷𝐸 = 0.75
Transición Slug-Annular a Annular-Mist.
𝛼𝑆𝐴 = 0.8
Transición Annular Mist y pre-CHF Mist (MPR).
𝛼𝐴𝑀 = 0.9999

ECUACIONES DE CIERRE EN FUNCIÓN DEL RÉGIMEN DE CAUDAL PARA TUBERÍA VERTICAL Y HORIZONTAL EN
MODELOS DE 6 ECUACIONES
FRICCIÓN PARED (𝑭𝒘𝒍 y 𝑭𝒘𝒗 ).
No se diferencia según regímenes:
2𝜌𝑘
𝐹𝑤𝑘 = 𝑓𝑤𝑘 · |𝑣𝑘 | · 𝑣𝑘 · 𝑉
𝐷ℎ
El factor de fricción 𝑓𝑤𝑘 :
64
𝑝𝑎𝑟𝑎 𝑅𝑒𝑘 ≤ 2000
𝑓𝑤𝑘 = { 𝑅𝑒𝑘
0.0055 + 0.55𝑅𝑒𝑘 −1/3 𝑝𝑎𝑟𝑎 𝑅𝑒𝑘 > 2000
NOTA:
- Para condiciones de régimen de caudal pre-CHF (Bubbly, Slug, Slug/Annular), se asume que el vapor no entra
en contacto con la pared, y por tanto siendo 𝑭𝒘𝒗 = 𝟎. El régimen Mist pre-CHF (MPR) se considerará como
post-CHF a estos efectos.
- Para condiciones de régimen de caudal post-CHF, se asume que el líquido no entra en contacto con la pared,
y por tanto siendo 𝑭𝒘𝒍 = 𝟎.

Para régimen Annular/Mist:

𝛼 − 𝛼𝑆𝐴 2𝜌𝑙
𝐹𝑤𝑙 = (1 − ) · 𝑓𝑤𝑙 · |𝑣𝑙 | · 𝑣𝑙 · 𝑉
𝛼𝐴𝑀 − 𝛼𝑆𝐴 𝐷ℎ
𝛼 − 𝛼𝑆𝐴 2𝜌𝑣
𝐹𝑤𝑣 = ( ) · 𝑓𝑤𝑣 · |𝑣𝑣 | · 𝑣𝑣 · 𝑉
𝛼𝐴𝑀 − 𝛼𝑆𝐴 𝐷ℎ
siendo 𝛼𝑆𝐴 y 𝛼𝐴𝑀 el criterio de fracción de huecos para cambio de Slug a Annular, y de Annular/Mist a Mist pre-CHF (MPR),
respectivamente.

′′
TRANSFERENCIA CALOR PARED (𝒒′′ 𝒘𝒗 y 𝒒𝒘𝒍 ).
La transferencia de calor en la pared se produce por convección (en 𝑘𝐽/𝑚2 𝑠), así:
′′ = ℎ (𝑇 − 𝑇 )
𝑞𝑤𝑣 𝑤𝑣 𝑣 𝑤
′′
𝑞𝑤𝑙 = ℎ𝑤𝑙 (𝑇𝑙 − 𝑇𝑤 )
Los coeficientes de convección de la pared (en 𝑘𝐽/𝑚2 𝑠𝐾) con el vapor y con el líquido se tienen que calcular según el
régimen en el que nos encontremos, acorde a las ecuaciones siguientes.
Líquido monofásico
𝑘𝑙
ℎ𝑤𝑙 = 7.86
𝐷ℎ
𝑘𝑙 , coeficiente de conductividad del líquido (tablas). 𝐷ℎ , diámetro hidráulico.

2
Ampliación de Energía y Máquinas Térmicas. Formulario Examen AEMT - 1º Parte. 2016-2017.

Bubbly (BBY), Slug (SLG), Slug/Annular o Churn (SLG/ANM).

- Líquido (correlación de Thom y Chen):
ℎ𝑤𝑙 = ℎ𝑓𝑐 + ℎ𝑛𝑏 · (1 − 𝐹𝑔𝑎𝑚 )
Con
𝑘𝑙
ℎ𝑓𝑐 = 0.023𝐹𝑐ℎ𝑒𝑛 𝑅𝑒 0.8 𝑃𝑟𝑙0.4
𝐷ℎ 𝑙
−1
1.0 𝑠𝑖 𝜒𝑇𝑇 < 0.1
𝐹𝑐ℎ𝑒𝑛 = { −1 −1
2.34(𝜒𝑇𝑇 + 0.213)0.736 𝑠𝑖 𝜒𝑇𝑇 > 0.1
−1
𝑥 0.9 𝜌𝑙 0.5 𝜇𝑣 0.1
𝜒 𝑇𝑇 =( ) ( ) ( )
1−𝑥 𝜌𝑣 𝜇𝑙
𝑥, calidad de caudal.
𝑃𝑟𝑙 , numero de Prandtl del líquido. 𝑅𝑒𝑙 , número de Reynolds del líquido.
−7
𝑒 2.4143·10 ·𝑝
ℎ𝑛𝑏 = 5.6782 · 10−3 { } (𝑇𝑤 − 𝑇𝑠𝑎𝑡 )
722
𝑝 es la presión del sistema.

- Vapor:
ℎ𝑤𝑣 = ℎ𝑛𝑏 · 𝐹𝑔𝑎𝑚
Siendo 𝐹𝑔𝑎𝑚 la fracción de calor generada en el vapor
𝑘
[−7.86 𝑙 · (𝑇𝑤 − 𝑇𝑙 )]
𝐷ℎ
𝐹𝑔𝑎𝑚 = 1 −
ℎ𝑛𝑏 · (𝑇𝑤 − 𝑇𝑠𝑎𝑡 )
Annular/mist (ANM).
- Líquido:
𝑘𝑙 ℎ𝑤𝑏
ℎ𝑤𝑙 = 7.86 +
𝐷ℎ max(1, 𝑇𝑤 − 𝑇𝑙 )
Con:
2 ′′
𝑇𝑤 − 𝑇𝑚𝑖𝑛 𝑞𝑐ℎ𝑓
ℎ𝑤𝑏 = ( )
𝑇𝑐ℎ𝑓 − 𝑇𝑚𝑖𝑛 𝑇𝑤 − 𝑇𝑙
′′
El flujo de calor igual al Critical Heat Flux, 𝑞𝑐ℎ𝑓 , se calcula interpolando con las tablas Groeneveld. Para la temperatura a
la que se alcanza el CHF, 𝑇𝑐ℎ𝑓 , se utiliza la correlación de Thom y Chen para Ebullición Nucleada. La temperatura mínima
de ebullición en película, 𝑇𝑚𝑖𝑛 , se obtiene con la correlación de Groeneveld.

- Vapor:
𝑘𝑣 0.6774 0.333
ℎ𝑤𝑣 = 0.07907 𝑅𝑒 𝑃𝑟𝑣
𝐷ℎ

Mist pre-CHF (MPR).

Se considera vapor monofásico.
Inverted Annular (IAN).
Se utiliza la correlación de Bromley + radiación (Stefan Boltzmann).
Transición Inverted Annular/Inverted Slug (IAN/ISL).
Se utiliza la correlación de Bromley + radiación (Stefan Boltzmann).
Inverted Slug (ISL).
Se utiliza la correlación de Bromley + radiación (Stefan Boltzmann).
Mist (MST) y Mist post-CHF (MPO).
Se consideran vapor monofásico.
Vapor monofásico.
𝑘𝑣
ℎ𝑤𝑣 = 0.07907 𝑅𝑒𝑣0.6774 𝑃𝑟𝑣0.333
𝐷ℎ

ÁREA INTERFACIAL (𝑨𝒔 ).

Bubbly (BBY).
𝑉
𝐴𝑠 = 3 · 𝛼 ·
𝑟𝑏𝑢𝑏𝑏𝑙𝑦
𝑟𝑏𝑢𝑏𝑏𝑙𝑦 es el radio medio de la burbuja bubbly, y se asume que es dependiente del criterio del número de Weber, con 𝑊𝑒𝑏
igual a 10:
0.5 · 𝑊𝑒𝑏 · 𝜎
𝑟𝑏𝑢𝑏𝑏𝑙𝑦 =
𝜌𝑙 · (𝑣𝑣 − 𝑣𝑙 )2

3
Ampliación de Energía y Máquinas Térmicas. Formulario Examen AEMT - 1º Parte. 2016-2017.

Slug (SLG).
𝑉
𝐴𝑠 = 3 · 𝛼 ·
𝑟𝑠𝑙𝑢𝑔
Con
𝑟𝑠𝑙𝑢𝑔 = 𝐹𝑏 · 𝑟𝑏𝑢𝑏𝑏𝑙𝑦 + (1 − 𝐹𝑏 ) · 𝑟𝑇𝑐
1
( )
3
1−𝛼 𝛼𝐵𝑆 ·(1−𝛼)
𝑟𝑇𝑐 = max {0.5 · 𝐷ℎ , [(𝛼 − 𝛼𝐵𝑆 · ( )) · 𝑉 · 0.2387] }, 𝐹𝑏 =
1−𝛼𝐵𝑆 𝛼·(1−𝛼𝐵𝑆 )

𝑟𝑇𝑐 es el radio medio de las Taylor Caps, y 𝛼𝐵𝑆 la fracción de huecos para el criterio de cambio de régimen de caudal Bubbly
a Slug.
Slug/Annular o Churn (SLG/ANM).
Igual que Slug.
Annular/Mist (ANM).
𝑇𝑚𝑖𝑛 − 𝑇𝑤 1−𝛼 𝑉
𝐴𝑠 = | |·( ) . 4 · √𝛼 ·
𝑇𝑚𝑖𝑛 𝛼 𝐷ℎ
Mist Pre-CHF (MPR).
Como Mist.
Inverted Annular (IAN).
𝑉
𝐴𝑠 = 4 · √(1 − 𝛼) ·
𝐷ℎ
Transición Inverted Annular/Inverted Slug (IAN/ISL).
Interpolación entre las correlaciones Inverted Annular e Inverted Slug.
Inverted Slug (ISL).
Interpolación entre la correlación Slug y la correlación Mist, pero cambiando (1 − 𝛼) por 𝛼 para el cálculo de 𝐴𝑠 .
Mist (MST).
𝑉
𝐴𝑠 = 3 · (1 − 𝛼) ·
𝑟𝑑𝑟𝑜𝑝𝑝𝑙𝑒𝑡
0.5 · 𝑊𝑒𝑏 · 𝜎
𝑟𝑑𝑟𝑜𝑝𝑝𝑙𝑒𝑡 =
𝜌𝑙 · (𝑣𝑣 − 𝑣𝑙 )2
Se asume que es dependiente del criterio del número de Weber, con 𝑊𝑒𝑏 igual a 4.
Mist post-CHF (MPO).
Como Mist.

NOTA: para una fase pura, vapor o líquido monofásico, 𝑨𝒔 = 𝟎.

FRICCIÓN INTERFACIAL (𝝉𝒔 ).

Fórmula general de la fricción interfacial para regímenes pre-CHF:
1
𝜏𝑠 = · 𝐶𝐷 · 𝜌𝑙 · |𝑣𝑣 − 𝑣𝑙 | · (𝑣𝑣 − 𝑣𝑙 ) · 𝐴𝑠
8
Para cada régimen se calcula un 𝐶𝐷 distinto.
Bubbly (BBY).
24
𝐶𝐷 = (1.0 + 0.1𝑅𝑒𝑏0.75 )
𝑅𝑒𝑏
donde
2 · 𝑟𝑏𝑢𝑏𝑏𝑙𝑦 · 𝜌𝑙 · |𝑣𝑣 − 𝑣𝑙 |
𝑅𝑒𝑏 =
𝜇𝑚𝑏
(𝜇 +0.4𝜇𝑙 )
−2.5 𝑣
𝜇𝑚𝑏 = 𝜇𝑙 (1 − 𝛼) (𝜇𝑣 +𝜇𝑙 )

Slug (SLG).
𝐶𝐷 = 0.45(1 − 𝛼)2
Slug/Annular o Churn (SLG/ANM).
8
𝐶𝐷 = (1 − 𝛼)2
3
Annular/Mist (ANM).
En este caso cambia la fórmula para el cálculo de la fricción interfacial:
1
𝜏𝑠 = · 𝑓𝑖𝑤 · 𝜌𝑣 · |𝑣𝑣 − 𝑣𝑙 | · (𝑣𝑣 − 𝑣𝑙 ) · 𝐴𝑠
2
Con:
𝑓𝑖𝑤 = 0.005(1 + 75(1 − 𝛼))
Mist Pre-CHF (MPR).
Como Mist.

4
Ampliación de Energía y Máquinas Térmicas. Formulario Examen AEMT - 1º Parte. 2016-2017.

Inverted Annular (IAN).

En este caso cambia la fórmula, deriva de la utilizada en régimen ANM.
1
𝜏𝑠 = · 𝑓𝑖𝑤 · 𝜌𝑙 · |𝑣𝑣 − 𝑣𝑙 | · (𝑣𝑣 − 𝑣𝑙 ) · 𝐴𝑠
2
Con 𝑓𝑖𝑤 = 0.005(1 + 75(1 − 𝛼)).
Transición Inverted Annular/Inverted Slug (IAN/ISL).
Interpolación entre las correlaciones Inverted Annular e Inverted Slug.
Inverted Slug (ISL).
Interpolación entre la correlación Slug y la correlación Mist, pero cambiando (1 − 𝛼) por 𝛼 para el cálculo de 𝜇𝑚 .
Mist (MST).
24
𝐶𝐷 = (1.0 + 0.15𝑅𝑒𝑑0.687 )
𝑅𝑒𝑑
(𝜇𝑙 +0.4𝜇𝑣 )
2·𝑟𝑑𝑟𝑜𝑝𝑝𝑙𝑒𝑡 ·𝜌𝑣 ·|𝑣𝑣 −𝑣𝑙 | 0.5·4·𝜎 −2.5
Con 𝑅𝑒𝑑 = , 𝑟𝑑𝑟𝑜𝑝𝑝𝑙𝑒𝑡 = y 𝜇𝑚 = 𝜇𝑣 𝛼 (𝜇𝑙 +𝜇𝑣 )
.
𝜇𝑚 𝜌𝑙 ·(𝑣𝑣 −𝑣𝑙 )2
Mist post-CHF (MPO).
Como Mist.

NOTA: para una fase pura, vapor o líquido monofásico, 𝝉𝒔 = 𝟎.

′′
TRANSFERENCIA CALOR INTERFACIAL (𝒒′′
𝒗 y 𝒒𝒍 ).
De forma general (en 𝑘𝐽/𝑚2 𝑠):
𝑞𝑣′′ = ℎ𝑣 · (𝑇𝑠𝑎𝑡 − 𝑇𝑣𝑠 )
𝑞𝑙′′ = ℎ𝑙 · (𝑇𝑠𝑎𝑡 − 𝑇𝑙𝑠 )
Donde ℎ𝑣 y ℎ𝑙 son los coeficientes de transferencia de calor interfacial (en 𝑘𝐽/𝑚2 𝑠𝐾), que se calcularán según el régimen.
𝑇𝑣𝑠 y 𝑇𝑙𝑠 se calcularán en función de la temperatura de pared 𝑇𝑤 con la siguiente aproximación:
𝑇𝑣 + 𝑇𝑤
𝑆𝑖 𝑇𝑤 > 𝑇𝑠𝑎𝑡 , 𝑇𝑣𝑠 = 𝑚𝑖𝑛 ( , 𝑇𝑣 + 2.8), 𝑇𝑙𝑠 = 𝑇𝑙 = 𝑇𝑠𝑎𝑡 𝐸𝑣𝑎𝑝𝑜𝑟𝑎𝑐𝑖ó𝑛
2
𝑇𝑙 + 𝑇𝑤
𝑆𝑖 𝑇𝑤 < 𝑇𝑠𝑎𝑡 , 𝑇𝑙𝑠 = 𝑚𝑎𝑥 ( , 𝑇𝑙 − 2.8), 𝑇𝑣𝑠 = 𝑇𝑣 = 𝑇𝑠𝑎𝑡 𝐶𝑜𝑛𝑑𝑒𝑛𝑠𝑎𝑐𝑖ó𝑛
2
Nota: lógicamente, si estamos en evaporación 𝒒′′ 𝒍 = 𝟎 y si estamos en condensación 𝒒′′ = 𝟎.
𝒗

Bubbly (BBY).
- Líquido: convección (+ conducción, conducción se desprecia).
𝑘𝑙
ℎ𝑙 = 1.1284√ (𝑣𝑣 − 𝑣𝑙 )𝜌𝑙 𝐶𝑝
𝐷𝑏
𝐷𝑏 , diámetro medio de la burbuja: 𝐷𝑏 = 2𝑟𝑏𝑢𝑏𝑏𝑙𝑦 .
- Vapor.
ℎ𝑣 = 1.0
Slug (SLG).
- Líquido
𝑘𝑙
ℎ𝑙 = 1.1284√ (𝑣𝑣 − 𝑣𝑙 )𝜌𝑙 𝐶𝑝
𝐷𝑏
𝐷𝑏 , diámetro medio de la burbuja: 𝐷𝑏 = 2𝑟𝑠𝑙𝑢𝑔 .
- Vapor.
𝑘𝑣
ℎ𝑣 = (2.0 + 0.74𝑅𝑒𝑣0.5 𝑃𝑟𝑣0.33 )
𝐷ℎ
Slug/Annular o Churn (SLG/ANM).
Interpolación entre Slug y Annular Mist (ANM).
- Líquido
ℎ𝑙,𝑠𝑙𝑔/𝑎𝑛𝑚 = 𝐹𝑏 · ℎ𝑙,𝑎𝑛𝑚 + (1 − 𝐹𝑏 ) · ℎ𝑙,𝑠𝑙𝑢𝑔
- Vapor
ℎ𝑣,𝑠𝑙𝑢𝑔/𝑎𝑛𝑚 = 𝐹𝑏 · ℎ𝑣,𝑎𝑛𝑚 + (1 − 𝐹𝑏 ) · ℎ𝑣,𝑠𝑙𝑢𝑔
𝛼−𝛼𝐵𝑆
Con 𝐹𝑏 = , siendo 𝛼𝐵𝑆 y 𝛼𝑆𝐴 el criterio de fracción de huecos para cambio de Bubbly a Slug y de Slug a Annular,
𝛼𝑆𝐴 −𝛼𝐵𝑆
respectivamente.

Annular/Mist (ANM).
- Líquido.
−2/3
ℎ𝑙 = 𝑓𝐻𝑀 𝜌𝑙 𝐶𝑝 |𝑣𝑣 − 𝑣𝑙 |𝑃𝑟𝑙 𝑓𝑎𝑖

5
Ampliación de Energía y Máquinas Térmicas. Formulario Examen AEMT - 1º Parte. 2016-2017.

𝑓𝐻𝑀 = 3.850𝑅𝑒 −2/3 , 𝑠𝑖 𝑅𝑒𝑙 < 1000

𝑓𝐻𝑀 = 0.5402𝑅𝑒 −0.38 , 𝑠𝑖 𝑅𝑒𝑙 ≥ 1000
𝑓𝑎𝑖 = 𝑚𝑎𝑥[0.0, min(1.0,111.1(0.999 − 𝛼))]
- Vapor
−2/3
ℎ𝑣 = 𝑓𝑖𝑤 𝜌𝑣 𝐶𝑝,𝑓𝑖𝑙𝑚 (𝑣𝑣 − 𝑣𝑙 )𝑃𝑟𝑣
𝑓𝑖𝑤 = 0.005(1 + 75(1 − 𝛼))
𝐶𝑝,𝑓𝑖𝑙𝑚 = 𝐶𝑝 (𝑝, 𝑇𝑓𝑖𝑙𝑚 ), es decir, se evalúa a la presión del sistema y a la temperatura de la película, 𝑇𝑓𝑖𝑙𝑚 =
𝑇 +𝑇
𝑚𝑖𝑛 ( 𝑣 𝑤 , 𝑇𝑣 + 2.8).
2
Mist Pre-CHF (MPR).
Como Mist.
Inverted Annular (IAN).
- Líquido
𝑘𝑙
ℎ𝑙 = 17.77
𝐷ℎ
- Vapor
𝑘𝑣
(2.0 + 0.55𝑅𝑒𝑣0.5 𝑃𝑟𝑣0.33 )
𝐷ℎ
ℎ𝑣 =
𝐹𝑦𝑐
0
ℎ𝑣0 −ℎ𝑣𝑠𝑎𝑡
𝐷ℎ , diámetro hidráulico de la tubería y 𝐹𝑦𝑐 = 1.0 + 0.5 .
ℎ𝑙0 −ℎ𝑣𝑠𝑎𝑡
0

𝑇𝑣 +𝑇𝑤
Cabe notar que 𝑘𝑣 y ℎ𝑣0 se evalúan a la presión del sistema y a la temperatura de la película, 𝑇𝑓𝑖𝑙𝑚 = 𝑚𝑖𝑛 ( , 𝑇𝑣 + 2.8).
2

Slug/Annular o Churn (SLG/ANM).

Se utiliza la correlación de Ley-Ryley para líquido y Yuen-Chen para vapor.
Inverted Slug (ISL).
Se utiliza la correlación de Ley-Ryley para líquido y Yuen-Chen para vapor.
Mist.
- Líquido
𝑘𝑙
ℎ𝑙 = 8.88
𝑟𝑑𝑟𝑜𝑝𝑝𝑙𝑒𝑡
𝑟𝑑𝑟𝑜𝑝𝑝𝑙𝑒𝑡 , radio de las gotas.
- Vapor
𝑘𝑣
(2.0 + 0.55𝑅𝑒𝑣0.5 𝑃𝑟𝑣0.33 )
𝐷ℎ
ℎ𝑣 =
𝐹𝑦𝑐
0
ℎ𝑣0 −ℎ𝑣𝑠𝑎𝑡
𝐷ℎ , diámetro hidráulico de la tubería y 𝐹𝑦𝑐 = 1.0 + 0.5 .
ℎ𝑙0 −ℎ𝑣𝑠𝑎𝑡
0

𝑇𝑣 +𝑇𝑤
Cabe notar que 𝑘𝑣 y ℎ𝑣0 se evalúan a la presión del sistema y a la temperatura de la película, 𝑇𝑣𝑠 = 𝑚𝑖𝑛 ( , 𝑇𝑣 + 2.8).
2

Mist post-CHF (MPO).

Como Mist.

′′
NOTA: para una fase pura, vapor o líquido monofásico, 𝒒′′
𝒗 y 𝒒𝒍 son nulos.

CONDENSACIÓN/EVAPORACIÓN (𝚪).
No se diferencia según regímenes:
(𝑞⃗𝑣′′ − 𝑞⃗𝑙′′ )𝑠 · 𝑛⃗⃗𝑣𝑠 𝐴𝑠
(ℎ𝑣0 − ℎ𝑙0 )𝑠 𝑘𝑔
Γ= ( 3 )
𝑉 𝑚 𝑠
0
Para condensación (ℎ𝑣0 )𝑠 = 𝑥 · ℎ𝑣,𝑠𝑎𝑡 y (ℎ𝑙0 )𝑠 = ℎ𝑙,𝑠𝑎𝑡
0
.
Para evaporación (ℎ𝑙0 ) = (1 − 𝑥) · ℎ𝑙,𝑠𝑎𝑡
𝑠
0 0
y (ℎ𝑣0 )𝑠 = ℎ𝑣,𝑠𝑎𝑡 .
𝑥, calidad de caudal.

NOMENCLATURA.
𝛼𝐵𝑆 , la fracción de huecos para el criterio de cambio de régimen de caudal bubbly a slug.
𝛼𝑆𝐴 , criterio de fracción de huecos para cambio régimen Slug a Annular.
𝜎, tensión superficial del fluido (tablas).
𝜇𝑣 y 𝜇𝑙 , viscosidades de vapor y líquido (tablas).
𝐴𝑤 , área de las paredes del volumen de control.

6
Ampliación de Energía y Máquinas Térmicas. Formulario Examen AEMT - 1º Parte. 2016-2017.

𝐶𝑝 , calor específico a presión constante (tablas).

𝐷𝑏 , diámetro medio de la burbuja: 𝐷𝑏 = 2𝑟𝑏𝑢𝑏𝑏𝑙𝑦 o 𝐷𝑏 = 2𝑟𝑠𝑙𝑢𝑔 , dependiendo del caso.
𝐷𝑑 , diámetro de las gotas: 𝐷𝑑 = 2𝑟𝑑𝑟𝑜𝑝𝑝𝑙𝑒𝑡 .
𝐷ℎ , diámetro hidráulico.
𝑔, aceleración de la gravedad.
𝑘𝑙 , coeficiente de conductividad del líquido (tablas).
𝑘𝑣 , coeficiente de conductividad del vapor (tablas).
𝐶𝑝 𝜇𝑘
𝑃𝑟𝑘 , número de Prandtl de la fase 𝑘, 𝑃𝑟𝑘 = .
𝑘𝑘
′′
𝑞𝑐ℎ𝑓 , flujo de calor igual al Critical Heat Flux (tablas Groeneveld)
𝑟𝑏𝑢𝑏𝑏𝑙𝑦 , es el radio medio de la burbuja bubbly.
𝑟𝑑𝑟𝑜𝑝𝑝𝑙𝑒𝑡 , radio de las gotas en régimen Mist.
𝑟𝑠𝑙𝑢𝑔 , es el radio medio de la burbuja slug.
𝑟𝑇𝑐 , es el radio medio de las Taylor Caps.
𝜌𝑘 |𝑣𝑘 |𝐷ℎ
𝑅𝑒𝑘 , número de Reynolds de la fase 𝑘, 𝑅𝑒𝑘 = .
𝜇𝑘
𝑇𝑐ℎ𝑓 , temperatura a la que se produce el Critical Heat Flux (tablas Groeneveld).
𝑇𝑚𝑖𝑛 , temperatura mínima de ebullición en película.
𝑇𝑠𝑎𝑡 , temperatura de saturación (tablas).
𝑇𝑤 , temperatura de la pared.
𝑉, volumen.
𝑊𝑒𝑏 , número de Weber.
𝑥, calidad de caudal.

ECUACIONES DE CIERRE EN CONDENSACIÓN EN PELÍCULA

TRANSFERENCIA CALOR PARED (𝒒′′

𝒘 ).

La transferencia de calor en la pared se produce por convección, así:

′′ = ℎ · (𝑇
𝑞𝑤 𝑤 𝑠𝑎𝑡 − 𝑇𝑤 ), 𝑘𝐽⁄𝑚2 𝑠

Los coeficientes de convección de la pared con el líquido se tienen que calcular según el régimen en el que nos
encontremos, acorde a las ecuaciones siguientes.
∗
ℎ𝑓𝑔 = ℎ𝑓𝑔 + 0.68𝐶𝑝𝑙 (𝑇𝑠𝑎𝑡 − 𝑇𝑤 ) + 𝐶𝑝𝑣 (𝑇𝑣 − 𝑇𝑠𝑎𝑡 ), (𝑘𝐽/𝑘𝑔)

Flujo laminar (calmado) sobre placas verticales.

∗ 1⁄4
𝑔𝜌𝑙 (𝜌𝑙 − 𝜌𝑣 )𝑘𝑙3 ℎ𝑓𝑔
ℎ𝑣𝑒𝑟𝑡 = 0.943 [ ] , (𝑘 𝐽⁄𝑚2 𝑠𝐾), 0 < 𝑅𝑒 < 30, 𝜌𝑙 ≫ 𝜌𝑣
𝜇𝑙 (𝑇𝑠𝑎𝑡 − 𝑇𝑠 )𝐿

Flujo laminar ondulado sobre placas verticales.

1⁄3 0.820
3.70𝐿𝑘𝑙 (𝑇𝑠𝑎𝑡 − 𝑇𝑤 ) 𝑔
𝑅𝑒𝑣𝑒𝑟𝑡,𝑤𝑎𝑣𝑦 = [4.81 + ∗ ( 2) ] , 𝜌𝑙 ≫ 𝜌𝑣
𝜇𝑙 ℎ𝑓𝑔 𝜈𝑙

1⁄3
𝑅𝑒𝑘𝑙 𝑔
ℎ𝑣𝑒𝑟𝑡,𝑤𝑎𝑣𝑦 ≅ ( ) , (𝑘𝐽⁄𝑚2 𝑠𝐾), 30 < Re < 1800, 𝜌𝑙 ≫ 𝜌𝑣
1.08𝑅𝑒1.22 − 5.2 𝜈𝑙2

Flujo turbulento sobre placas verticales.

1⁄3 4⁄3
0.0690𝐿𝑘𝑙 𝑃𝑟 0.5 (𝑇𝑠𝑎𝑡 − 𝑇𝑤 ) 𝑔
𝑅𝑒𝑣𝑒𝑟𝑡,𝑡𝑢𝑟𝑏𝑢𝑙𝑒𝑛𝑡𝑜 =[ ∗ ( 2) − 151𝑃𝑟 0.5 + 253] , 𝜌𝑙 ≫ 𝜌𝑣
𝜇𝑙 ℎ𝑓𝑔 𝜈𝑙

1⁄3
𝑅𝑒𝑘𝑙 𝑔
ℎ𝑣𝑒𝑟𝑡,𝑡𝑢𝑟𝑏𝑢𝑙𝑒𝑛𝑡𝑜 = ( ) , (𝑘𝐽⁄𝑚2 𝑠𝐾), 𝜌𝑙 ≫ 𝜌𝑣
8750 + 58𝑃𝑟 −0.5 (𝑅𝑒 0.75 − 253) 𝜈𝑙2

7
Ampliación de Energía y Máquinas Térmicas. Formulario Examen AEMT - 1º Parte. 2016-2017.

Placas inclinadas.

ℎ𝑖𝑛𝑐𝑙𝑖𝑛𝑎𝑑𝑜 = ℎ𝑣𝑒𝑟𝑡 (𝑐𝑜𝑠𝜃)1⁄4 (𝑘𝐽⁄𝑚2 𝑠𝐾), (𝑙𝑎𝑚𝑖𝑛𝑎𝑟)

Condensación en tubos horizontales.

∗ 1⁄4
𝑔𝜌𝑙 (𝜌𝑙 − 𝜌𝑣 )𝑘𝑙3 ℎ𝑓𝑔
ℎℎ𝑜𝑟𝑖𝑧,𝑁 𝑡𝑢𝑏𝑜𝑠 = 0.729 [ ] , (𝑘𝐽⁄𝑚2 𝑠𝐾)
𝜇𝑙 (𝑇𝑠𝑎𝑡 − 𝑇𝑤 )𝑁𝐷ℎ

TASA DE CONDENSACIÓN (𝚪).

′′ · 𝐴
𝑞𝑤 𝑤
∗
𝑚̇ ℎ𝑓𝑔
Γ=− =− (𝑘𝑔/𝑚3 𝑠)
𝑉 𝑉
𝑻𝒔𝒂𝒕 +𝑻𝒘
NOTA: Todas las propiedades del líquido se deben evaluar a la temperatura de película 𝑻𝒇𝒊𝒍𝒎 = ( ). El 𝒉𝒇𝒈 y la
𝟐
𝝆𝒗 se deben evaluar a la temperatura de saturación 𝑻𝒔𝒂𝒕 .

NOMENCLATURA.

𝜃, ángulo respecto de la vertical.

𝜇𝑣 y 𝜇𝑙 , viscosidades de vapor y líquido (tablas).
𝐴𝑤 , área superficial de transferencia de calor (de la lámina o del tubo/s).
𝐶𝑝𝑙 , calor específico del líquido (a 𝑇𝑓𝑖𝑙𝑚 ) a presión constante (tablas).
𝐶𝑝𝑣 , calor específico del vapor a presión constante (tablas).
𝐷ℎ , diámetro hidráulico.
𝑔, aceleración de la gravedad.
ℎ𝑓𝑔 , calor latente de vaporización.
∗
ℎ𝑓𝑔 , calor latente de vaporización modificado.
𝑘𝑙 , coeficiente de conductividad del líquido (tablas).
𝐿, altura vertical de la placa.
𝑚̇, caudal másico de concentrado.
𝑁, número de tubos.
𝐶𝑝 𝜇
𝑃𝑟, número de Prandtl del líquido, 𝑃𝑟 = .
𝑘
𝜌|𝑣|𝐷ℎ
𝑅𝑒, número de Reynolds del líquido, 𝑅𝑒 = .
𝜇
𝜈𝑙 , viscosidad dinámica del líquido.
𝑇𝑓𝑖𝑙𝑚 , temperatura promedio de la película.
𝑇𝑠𝑎𝑡 , temperatura de saturación (tablas).
𝑇𝑣 , temperatura del vapor.
𝑇𝑤 , temperatura de la pared (𝑇𝑠 , en apuntes de teoría).
𝑉, volumen del volumen de control.