Equipo 6.
Cálculo de fuga de vapor a través de un orificio
Fórmulas utilizadas para el proyecto de fuga de vapor a través de una válvula.
1. Cálculo de densidad del fluido inicial (Redlich/Kwong)
T P
Tr= Pr=
Tc Pc
Ψα ( Tr) Pr
q= β=Ω
ΩTr Tr
(Z− β)
z=1+ β−qβ
Z (Z + β )
Teniendo el factor de compresibilidad se procede a calcular la densidad del fluido.
PM
ρ=
ZRT
ρ = densidad del fluido (g/L)
PM = peso molecular del compuesto
Z = el factor de compresibilidad
R = constante de gases ideales (L*bar/K*mol)
T = temperatura (K)
2. Cálculo del área de apertura de la válvula
Primero se procede a calcular el Cp y Cv para cierta temperatura (Tabla A-2, Cengel)
2 3
Cp=a+ bT + c T + d T
T = temperatura (K)
Cp = capacidad calorífica (KJ/kmol*K)
Se calcula el Cv con el Cp obtenido usando una R de 8.3144 J/mol*K
Cp−Cv=R
Conociendo nuestras capacidades caloríficas, calculamos k
Cp
k=
Cv
Para conocer el caudal máximo en una tubería se usa:
k+1 1/ 2
( pm )∗k 2
G=P[
RT ( )
1+ k
k−1
]
Donde:
P = Es la presión inicial (Pa)
Pm = Peso molecular (kg/mol)
R = constante de gases ideales (J/mol* K)
T = temperatura del fluido (K)
G = caudal máximo en el orificio (kg/m2*s)
Teniendo el caudal máximo en le orificio se procede a calcular el área máxima del orificio
relacionandola con el flujo másico inicial (el área real deberá ser menor a la máxima).
Amax =ḿ/G
Conociendo el área máxima se insertan valores de áreas de orificios para conocer la fuga de
vapor en t/h.
Fugav=( Aorificio∗G )∗3.6
Teniendo nuestra k y las condiciones de entrada se pueden conocer las
condiciones de salida del fluido con las siguientes ecuaciones
3. Cálculo de conducción y convección después de la válvula
Una vez obtenido el flujo de la fuga se procede a calcular la conducción y convección en la
tubería. Para la conducción se usan 3 fórmulas que dependerán del tipo de flujo que se
maneje (laminar, de transición o turbulento).
Para flujo laminar se usa la correlación de Sieder-Tate:
.33
[ ( )] ( )
.14
hD D μ
=1.86 ℜ∗Pr
k L μw
Para un flujo de transición se utiliza la correlación de Hilpert:
)[ ( ) ] ( )
.66
(
.66 .14
h ℜ −125 D μ
=.116 1+ Pr −.66
cρv ℜ L μw
Para un flujo turbulento se utiliza la correlación de Dittus-Boettler:
.14
hD μ
k
=.023 ℜ.8∗Pr .33
μw ( )
Para la convección se utiliza la ecuación de Sutherland para calcular la visocidad, se
obtiene el número de Grashof y posteriormente el número de Nusselt para asi calcular el
coeficiente de transferencia de calor externa.
Ecuación de Sutherland:
C T 3/ 2
μ=
T +S
Número de Grashof:
3 2
L ρ gβ (T O−T ∞ )
Gr= 2
μ
Número de Nusselt:
hD 1 /4
Nu= =.53(Gr D∗Pr)
k
hD
Nu= =.114 (Gr D∗Pr)1/ 3
k
Dependiendo del valor de la multiplicación del número de Prandalt y Grashof se usarán las
fórmulas
Una vez teniendo nuestros coeficientes de transferencia de calor y nuestro coeficiente
global de transferencia de calor, que se calculó a partir de las resistencias, calculamos la
pérdida de calor a través de la pared.
Q́=∆T ∗U Donde Q́ esta en W/m2
Con nuestra temperatura en la superficie exterior de la tubería se puede saber la cantidad de
vapor que se esta fugando, una vez teniendo el flujo de vapor, usamos la entalpia para
calcular asi la energia que se esta perdiendo debido a la fuga.
Energía=∆ H∗Fugav( flujo másico)
Referencias
Cao, Eduardo (2010), Heat Transfer in Process Engineering, McGraw-Hill
Van Ness (2017), Introducción a la termodinámica en ingenieria quimica
Levenspiel (1993), Flujo de fluidos e intercambio de calor
Cengel (2012), Termodinámica
Shankar, Heat transfer in Flow Through Conduits
Ríos (2011), Determinación del calor de vaporización del agua
FLUENT, Viscosity as a Function of Temperature
Alan [Link] (1990). Transmisión del calor
