Saltar al contenido principal

k = exponente isoentrópico

La importancia de  k  para válvula de seguridad

Editado por Alejandro Ruzza 

El dimensionamiento de válvulas de seguridad diseñadas para la descarga de gases o vapores, según la Colección “E” de la lspesl, requiere el conocimiento del exponente isoentrópico k en condiciones de descarga.

La aplicación descuidada del capítulo “E.1” de la Colección “E” de lspesl, relativo al dimensionamiento de las válvulas de seguridad, puede conducir a una sobreestimación de la capacidad de descarga de las válvulas y los discos de ruptura.

Este artículo da algunas pautas para estimar el valor de k para gases reales y
destaca el error al considerar k igual a la relación de calores específicos Cp/Cv

Un primer y grave error a evitar es utilizar la fórmula de la Colección 'E', válida para gases o vapores, en situaciones en las que un descarga de dos fases de líquido y gas/vapor tiene lugar. En tales casos, de hecho, los diámetros calculados serán sin duda subdimensionados en comparación con la necesidad real.
Un segundo error, que en muchos casos puede conducir a la subdimensionar el sistema de seguridad, es dar al exponente isoentrópico k el valor de la relación Cp/Cv. Si bien el primer punto será objeto de una serie de artículos posteriores, aquí nos gustaría dar algunos consejos útiles para calcular el exponente isoentrópico y mostrar, en casos concretos, la magnitud del error que se puede cometer.

Salida isoentrópica a través de una tobera

 

La fórmula [ 1 ] que se usa en la colección "E", así como en otros italianos [ 2 ] y extranjero [ 3 ] standards, para el cálculo de las válvulas de seguridad que deben descargar gases o vapores, es el de la salida isoentrópica por una tobera en condiciones críticas de salto, que para un gas ideal es:

Fórmula lspesl Colección “E”

donde la experienciaansisobre el coeficiente C viene dado por:

expansisobre el coeficiente C

"Ser" k el exponente de la isoentrópica expansien ecuación: pxv^k=costo

FluidoP1 (bar)T1 (ºC)q' (kg/h)q (kg/hora)(q'/q) x 100
Metano125014721466100.4
Metano2320023142267102.1
Propano1210022612181103.7
hexano1217830992740113.1
hexano2322065195111127.5
Heptano1221532322821114.4

q'= caudal calculado con k = Cp/Cv (20 °C, 1 atm)
q = caudal calculado con k = (Cp/Cv) • (Z/Zp)

Introduciendo el coeficiente experimental k de salida de la válvula de seguridad, que considera globalmente el rendimiento real de salida de la válvula, un coeficiente de seguridad de 0.9 y el factor de compresibilidad Z1 para el fluido real, llegamos a la formulación de la colección “E”:

(1) [ 1 ]

El exponente isoentrópico k se puede expresar como:

[ 2 ] [ 2 ]

Para una gas idealpara lo cual P x V / R x T = 1 , se demuestra que k es igual a la relación Cp/Cv entre los calores específicos a presión y volumen constantes.

Para una gasolina de verdad, k se puede expresar (ver Apéndice B) por:

[ 3 ] [ 3 ]

donde Z es el factor de compresibilidad definido por Z=P x V / R x T y Zp es el “factor de compresibilidad derivado”. Al aplicar la fórmula [ 3 ], según la colección “E”, los valores de Cp/Cv, Z y Zp deben evaluarse en las condiciones de descarga P1 y T1.

El factor de compresibilidad derivado Zp se define en la fórmula [ 4 ] como:

[ 3.1 ]

El factor de compresibilidad Z se puede expresar como:

[ 4 ][ 4 ]

y de manera similar, se puede expresar como:

[ 5 ][ 5 ]

donde los valores de Z^0, Z^1, Zp^0, Zp^1 están tabulados en el Apéndice A en función de Pr y Tr.

In [ 4 ] y [ 5 ], Ω es el factor acéntrico de Pitzer definido por:

[ 10 ] [ 10 ]

Donde Pr^SAT es la presión de vapor reducida correspondiente a un valor de temperatura reducida Tr=T/Tc=0,7. El Apéndice A muestra los valores de Ω de algunos fluidos. Z e Zp también se puede derivar directamente de una ecuación de estado analítica.

Un ejemplo numérico

 

Pasando a un ejemplo numérico, supongamos que necesitamos calcular la capacidad de descarga de una válvula de seguridad bajo las siguientes condiciones:

Fluidon-butano
Estado fisicovapor sobrecalentado
Masa molecularM58,119
Configurar presiónP19,78 bares
Presión demasiada10%
Temperatura del fluidoT400 K
Coeficiente de salida0,9
Diámetro del orificioDo100 mm

la presión de descarga viene dada por:

siendo para n-Butano: Tc=425,18 K y Pc=37,96 bar, tenemos:

y usando las tablas del Apéndice A, tenemos:

Conociendo el volumen específico del vapor en las condiciones de descarga (P1, T1) igual a 0,01634 m^3/kg (0,0009498 m^3/g-mol), también podríamos haber calculado Z a partir de:

Dada la relación de los calores específicos a presión y volumen constantes, en condiciones de descarga (P1, T1), igual a 1,36, de la fórmula [ 3 ] tenemos:

147060

Aplicando la fórmula [1], con cálculo del caudal

Aplicando fórmula [ 1 ], que se resolvió para el cálculo del caudal, tenemos un valor de caudal de descarga de 147.060 kg / h.

174848

Aplicando la fórmula [1], usando el valor de Cp/Cv a 1 atm y 20 °C

Si en su lugar hubiéramos utilizado el valor de Cp/Cv a 1 atm y 20 °C, habríamos tenido k = 1,19 y de fórmula [ 1 ] un caudal de descarga de 174.848 kg / h.

Esto nos hubiera llevado a sobrestimar la descarga capacidad de la válvula de seguridad en alrededor 19%

ADVERTENCIA:

El error que se puede cometer al asignar el valor Cp/Cv a k puede ser mucho mayor que en este ejemplo.

MÁS DEL 20%

Para dar una idea, en la siguiente tabla se muestran los caudales de un orificio de 18 mm para otros hidrocarburos saturados, calculados en los dos casos. Los cálculos se realizaron con especial developed software.

FluidoP1 (bar)T1 (ºC)q' (kg/h)q (kg/hora)(q'/q) x 100
Metano125014721466100.4
Metano2320023142267102.1
Propano1210022612181103.7
hexano1217830992740113.1
hexano2322065195111127.5
Heptano1221532322821114.4

El software no utiliza fórmulas. [ 4 ] [ 5 ] pero, a partir de la modificada Ecuación de estado de Redlich y Kwong, calcula el valor del exponente isoentrópico utilizando correlaciones termodinámicas.

Apéndice A y B
derivación de fórmulas

BESA estará presente en el IVS - IVS Industrial Valve Summit 2024