NTP 457: Discos de ruptura (II): dimensionado

Disques de rupture (II): dimensionnement
Bursting discs (II): sizing Discos de ruptura

Redactor:

Emilio Turmo Sierra
Ingeniero Industrial

CENTRO NACIONAL DE CONDICIONES DE TRABAJO

Introducción

Para calcular el área necesaria del dispositivo de protección de un recipiente por alivio de presión mediante la descarga de fluido, se debe partir del tipo de fenómeno a proteger.

1. Para la protección contra explosiones por combustión de polvos, se pueden emplear los nomogramas presentados en la NTP-428, teniendo en cuenta todos los parámetros necesarios y las excepciones indicadas.

2. Para las explosiones por combustión de gases y vapores inflamables (deflagraciones) existen unos nomogramas similares, incluidos en el código 68 del National Fire Protection Association (NFPA), correspondientes a los gases: metano, propano, gas de coque e hidrógeno. En este caso se debe conocer el valor de la constante K_G del gas considerado, mediante el ensayo correspondiente y cuyo valor se determina según la expresión de la ¨ley cúbica¨ (dp/dt)_máx · V^1/3 = K_G en la que (dp/dt)_máx es el gradiente de presión máximo y V el volumen del recipiente normalizado de ensayo. Los valores de K_G , para que sean válidos y comparativos, se deben calcular a partir del ensayo realizado en un recipiente de igual volumen y con un foco de ignición de igual energía. Utilizar una de las gráficas de los cuatro gases patrones indicados cuya K_G se aproxime más por exceso o realizar una interpolación. En caso de no disponer de datos de ensayo se puede utilizar el nomograma del hidrógeno, ya que el área adicional de venteo suele ser pequeña. Se debe partir del volumen del recipiente, de la presión reducida P_red según la resistencia del recipiente y de la presión estática P_stat del dispositivo de venteo. Debe tenerse en cuenta que estos nomogramas sirven para unas condiciones iniciales de no turbulencia en el momento de la ignición, no turbulencia por obstáculos internos, un foco de ignición £ 10 J, situación a presión atmosférica y con una razón longitud L / diámetro D £ 5.

   K_G metano = 55 bar.m.s^-1

   K_G propano = 75 bar.m.s^-1

   K_G gas de coque = 140 bar.m.s^-1

   K_G hidrógeno = 550 bar.m.s^-1

3. Para la protección contra estallidos por aumentos de presión, se efectúa un dimensionado a partir de unas expresiones deducidas de la mecánica de fluidos y adoptadas por los códigos de diseño más importantes. Entre estos hay que destacar la Norma ISO 6718: 1991, el Código ASME, la Norma API RP 520, la Norma BS 2915: 1990, la Norma AD-Merkblatt A 1, los Códigos NFPA (30, 58), etc.

Se hace la distinción entre líquidos por un lado y gases y vapores por otro. El caudal de alivio requerido puede venir impuesto por reglamentos específicos o se puede obtener examinando las posibles causas de sobrepresión, tales como fallo de la energía eléctrica, fallo de la refrigeración, incendio, etc, y haciendo el cálculo para cada una de las contingencias. Se puede recurrir a las Instrucciones Técnicas Complementarias del Reglamento sobre Almacenamiento de Productos Químicos, Reglamento de aparatos a presión y a las Normas y Códigos citados anteriormente, para disponer de recomendaciones o fórmulas de cálculo válidas para válvulas de seguridad y discos de ruptura. Por ejemplo, en el caso de incendio se deduce a partir del calor recibido y del calor o entalpía de vaporización del fluido (ver NTP-346). El caudal de alivio generalmente deberá ser el mayor de todos ellos. En bibliografía especializadas también se citan programas de cálculo para diversas situaciones. Un ejemplo sería el caso de flujo en doble fase, el cual entre otros documentos se trata en la Norma ISO 4126.-1:1991.

En los discos de ruptura, la presión máxima de ruptura no debe exceder a la presión de diseño del recipiente y la presión mínima de ruptura debe exceder a la presión máxima de trabajo o servicio en un cierto margen, que depende del tipo de disco seleccionado.

Dimensionado de discos de ruptura para líquidos

Según la norma ISO 6718: 1991, el área de la sección del disco de ruptura debería ser:

A₀ = 0,6211 q_m / [f_u · a · (Dp · r)^1/2 ]

Siendo:

A₀ = área del disco de ruptura (mm²).

f_u = factor de corrección por viscosidad. Si la viscosidad es igual o menor que la del agua a 20° C, f_u = 1. Este factor depende del número de Reynolds Re y se puede obtener del gráfico de la figura 1. El número de Reynolds se obtiene de la expresión Re = 0,3134 · q_m / [m (A₀ )^1/2 ] siendo m = viscosidad dinámica en Pa . s, lo cual nos indica que se ha de considerar un valor de partida de q_m o A₀ y hacer tanteos iterativos.

Figura 1. Factor de corrección f_u de la viscosidad dinámica del líquido según el Número de Reynolds Re.

Dp = diferencia de presión entre el sistema a presión y la presión en el extremo del sistema de descarga o la atmósfera en bar.

q_m = flujo o caudal másico a descargar en kg/h.

a = coeficiente de descarga de la boquilla de salida (0,62 o según reglamentos o normas determinadas).

r = densidad en kg/m³.

Ejemplo

Dimensionar el disco de ruptura que se requiere para proteger un recipiente que contiene un líquido con las condiciones siguientes: caudal de alivio requerido 2000 L/min; densidad r = 1300 kg/m³ ; presión máxima de trabajo permitida = 15 bar manométrico; coeficiente de descarga a = 0,62; f_u = 1.

Solución:

La normativa española del Reglamento de Aparatos a Presión y del Reglamento sobre Almacenamiento de Productos Químicos, igual que la norma ISO 6718:1991 y otras, establece que el caudal de descarga debe ser tal que, durante la actuación del dispositivo de protección, la presión máxima de alivio alcanzada no sobrepase el 10% de la presión máxima de trabajo permitida, con lo que: presión máxima de alivio permitida = 16,5 bar y flujo máximo de alivio requerido = 156.000 kg/h

El área de disco de ruptura debería ser como mínimo:

A₀ = 0,6211[156.000 / [1 · 0,62 (16,5 · 1300)^1/2 ]] = 1067 mm²

equivalente a un diámetro de 36,9 mm. Se adoptaría el diámetro nominal normalizado superior DN 40 (1 1/2”)

Dimensionado de discos de ruptura para gases y vapores

Considerando p la presión absoluta en el interior del recipiente y p_b la presión absoluta exterior o contrapresión, cuando p_b/p = 1 no hay flujo. Al incrementar la presión p el flujo va aumentando y decrece la razón p_b/p hasta llegar a un valor crítico que se llama razón crítica r_crit = [2 / (k + 1)]^k/(k-1) siendo k = c_p/c_v la razón de los calores específicos a presión y volumen constante. En ese punto el flujo alcanza la velocidad del sonido. La consecuencia es que se deben considerar dos casos de flujo: el sónico (condiciones críticas) cuando p_b/p£ [2 / (k + 1)]^k/(k-1) y el subsónico (condiciones subcríticas) cuando p_b/p > [2 / (k + 1)]^k/(k-1) .

a) Fluidos compresibles (gases y vapores). Capacidad de descarga con flujo subsónico (subcrítico).

El área de la sección del disco de ruptura según fórmulas derivadas de la mecánica de fluidos sería:

A₀ = 3,469 [[q_m/(C k_b a )] [v / p]^1/2 ]

o su equivalente:

A₀ = [q_m / (C k_b a p )] [T Z / M]^1/2 (ISO 6718:1991)

siendo:

A₀ = área del disco de ruptura (mm²)

C = 3,948 [k [2 / (k + 1)^{(k+1) / (k-1)} ]^1/2 (ver tabla 3)

k_b = [2k / (k-1) [(p_b/ p)^2/k-(p_b/p)^(k+1)/k] / [k[2/(k+1)]^(k+1)/(k-1)]]^1/2 es un factor de corrección que reduce la capacidad de descarga debido a incrementos de la contra presión. Para flujo crítico k_b = 1. La tabla 1 da valores calculados de k_b.

TABLA 1. Factor de corrección k_b de la capacidad de descarga según la contrapresión.

k = exponente isoentrópico = c_p/c_v

M = masa molecular en kg/kmol

p = presión de alivio en bar absolutos

p_b = contrapresión inmediatamente aguas abajo del área de la sección transversal mínima en bar absolutos

q_m = flujo máximo a descargar o la capacidad requerida de descarga del disco de ruptura en kg/h

T = temperatura absoluta de alivio en K

v = volumen específico a la presión y temperatura de alivio en m³/kg

Z = factor de compresibilidad (ver fig. 2).

Fig. 2: Factor de compresibilidad Z en función de la presión y temperatura reducidas

a = coeficiente de descarga según tabla 2 y debiéndose cumplir 0,5 A₁ £ A₀ £ A₁ siendo A₁ el área de la sección transversal del conducto de entrada en mm² . Para configuraciones distintas de las de la tabla 2 se requieren ensayos. El coeficiente establecido experimentalmente debe multiplicarse por 0,9 antes de aplicar las fórmulas de cálculo de A₀.

Tabla 2. Coeficientes de descarga a

Nº	Tipo de boquilla		Coeficiente de descarga a (fluidos compresibles)
1		Boquilla saliente	0,68
2		Boquilla encastrada o rebajada y también una brida maciza (formando bloque) con diseño de configuración no hidrodinámica	0,73
3		Brida maciza (formando bloque) de configuración hidrodinámica, p.e. con bordes de entrada redondeados o achaflanados y también con orificio rebordeado hacia fuera.	0,80

Ejemplo

Dimensionar el disco de ruptura que se requiere para proteger un recipiente a presión que contiene gas nitrógeno en las condiciones de proceso siguientes: presión máxima de trabajo permitida: 5 bar manométrico; caudal de alivio requerido: 300.000 L/min en las condiciones de salida; temperatura del flujo de salida: 200 °C; contrapresión del sistema durante el venteo: 4,5 bar manométrico; coeficiente efectivo de descarga a = 0,62; coeficiente isoentrópico k = c_p/c_v = 1,404; masa molecular del N₂ : M = 28; temperatura crítica T_c = -146, 95 °C = 126,2 K; presión crítica P_c = 33,94 bar

Solución

Presión máxima de alivio: p = 6,513 bar abs.; p_b = 5,513 bar abs.

p_b/p = 0,846; r_crit= [2 / (1,404 + 1)]^{1,404/(1,404-1)} = 0,528

Al ser 0,846 > 0,528 se deduce que el flujo es subsónico.

Se deben calcular los parámetros que intervienen en la fórmula que da A₀. El valor q_m en kg/h se calcula a partir de los 300.000 L/min dados en las condiciones de salida a 200 °C y 6,513 bar, aplicando las leyes de los gases perfectos:

p V / T = p₀ V₀ /T₀

(6,513 / 1,013) ·18.000.000 / (273,15 + 200) = 1·V₀/ 273,15

V₀ = 66.796.853 L/h = 66 796 m³/h (a 0 °C y 1 atm).

En estas condiciones 1 mol ocupa 22,4 L:

(66 796 853 L/h) / (22,4 L/mol) = 2.982.002,4 mol/h

Y al ser 28 el peso molecular del nitrógeno, tendremos:

q_m = 83 496,066 kg/h

El volumen específico v en las condiciones de salida (200 °C y 6,513 bar) será:

v = (18.000 m³/h) / (83.496,066 kg/h) = 0,215 579 m³/kg

La constante C = 2,70 (Tabla 3 para k = 1,404). De la tabla 1 para k = 1,4 y p_b / p = 0,84 se obtiene el valor aproximado k_b = 0,752. Aplicando la fórmula correspondiente daría:

k_b = 0,740

La sección mínima del disco de ruptura sería:

A₀ = 3,469 [[83.496,066/(2,70 · 0,740 · 0,62 )] [0,215579/ 6,513]^1/2 ] = 42.539,8 mm²

que corresponde a un diámetro de 232,7 mm. Se tomaría un diámetro nominal DN 250 (10”).

b) Fluidos compresibles (gases y vapores). Capacidad de descarga con flujo sónico (crítico)

En este caso la razón de las presiones p_b/p £ r_crit y la fórmula a emplear es según la norma ISO 6718:1991:

A₀ = 3,469 [[q_m / (C a )] [v / p]^1/2 ]

o su equivalente A₀ = [q_m / (C a p )] [T Z / M]^1/2 siendo los parámetros los mismos que en el caso anterior, excepto el k_b que por ser igual a la unidad no aparece.

Ejemplo

Dimensionar el disco de ruptura que se requiere para proteger un recipiente a presión que contiene gas helio en las condiciones siguientes: presión máxima de trabajo permitida: 10 bar manométrico; caudal de alivio requerido: 1.500.000 L/min medidos a 1 atm y 15,5° C; contrapresión del sistema durante el venteo: 0 bar manométrico (presión atmosférica); temperatura del flujo: 40° C; coeficiente de descarga: a = 0,62; coeficiente isoentrópico: k = c_p / c_v = 1,66; masa molecular del He: M = 4; temperatura crítica T_c = -267,95° C; presión crítica P_c = 2,275 bar

Solución

p_b = 1,013 bar abs.; p = 12,013 bar abs.

p_b / p = 1,013 / 12,013 = 0,084

r_crit = [2 / (1,66 + 1)]^{1,66/(1,66-1)} = 0,488

p_b / p = 0,084 < r_crit = 0,488. El flujo es sónico

Para aplicar la expresión A₀ = 3,469 [[q_m / (C a )] [v / p]^1/2 ] es necesario calcular los parámetros en los pasos siguientes:

El flujo másico q_m en kg/h se puede conocer a partir del número de moles que corresponden al caudal de alivio requerido. Aplicando la ley de los gases perfectos p V / T = p₀ V₀ /T₀.

El volumen V₀ a 0° C y 1 atm sería:

V₀ = p V T₀ / (p₀ T) =1 · 1.500.000 · 273,15 / [1 · (273,15 + 15,5)] = 1 419 452,6 L

y en flujo volumétrico 1.419.452,6 L/min. Sabiendo que 1 mol de helio = 4g ocupa 22,4 L a 0° C y 1 atm, el flujo másico de helio sería:

q_m = 15.208, 421 kg / h

El valor de C se puede calcular con la fórmula dada en el apartado a) Fluidos compresibles (gases y vapores). Capacidad de descarga con flujo subsónico (subcrítico) o con la tabla 3 para k = 1,66: C = 2,86; p = 12,013 bar absolutos (valor a sustituir en A₀); p = 12,013 bar a / (1,013 bar/atm) = 11,856 atm absolutas; T =40° C = 273,15 + 40 = 313,15 K

El caudal de alivio requerido de 1.500.000 L/min. = 90.000 m³/h (a 1 atm y 15,5°C). Aplicando la ley de los gases perfectos p V / T = p₀ V₀ /T₀:

11,856 · V / (273,15 + 40) = 1 · 90.000 / (273,15 + 15,5); V = 8 235,30 m³/h

El volumen específico v sería:

v = (8.235,30 m³/h) / (15.208,421 kg/h) = 0,541 m³/kg

La sección mínima del disco de ruptura debería ser:

A₀ = 3,469 [[15.208,421 / (2,86 · 0,62)] [0,541 / 12,013]^1/2 ] = 6.316,81 mm²

Equivalente a un diámetro de 89,68 mm. Se tomaría el diámetro nominal inmediato superior DN 100 (4”).

Bibliografía

(1) LEES, F.P.
Loss Prevention in the Process Industries
Oxford, ButterworthHeinemann, 1996, 3 Vols.

(2) PARRY, C.F.
Relief Systems Handbook
Rugby, England, Institution of Chemical Engineers, 1992

(3) TURMO, E.
NTP 427 y 428. Paramentos débiles para el venteo de alivio de explosiones (I y II)
Instituto Nacional de Seguridad e Higiene en el Trabajo, 1996

(4) CEJALVO, A.
NTP 346. Válvulas de seguridad (II). Capacidad de alivio y dimensionado
Instituto Nacional de Seguridad e Higiene en el Trabajo, 1994

(5) Reglamento sobre Instalaciones de almacenamiento de gases licuados del petróleo (GLP) en depósitos fijos
Madrid. Mº de Industria, Comercio y Turismo

(6) Reglamento sobre Almacenamiento de Productos Químicos y sus Instrucciones Técnicas Complementarias
Madrid, Mº de Industria, Comercio y Turismo

(7) Reglamento de Aparatos a Presión y sus Instrucciones Técnicas Complementarias
Madrid, Mº de Industria, Comercio y Turismo

(8) Código NFPA 68.
Guide for Venting of Deflagrations

(9) Código NFPA 30.
Flammable and Combustible Liquids Code

(10) Código NFPA 58.
Standard for the Storage and Handling of Liquefied Petroleum Gases
Quincy, Massachusetts, National Fire Protection Association, 1989.

(11) Norma ISO 6718: 1991.
Bursting discs and bursting disc devices

(12) Norma ISO 41261:1991.
Safety valves. Part 1. General requirements

(13) Código ASME. Section VIII. Division 1.
Pressure Relief Devices

(14) Norma API RF 520.
Sizing, Selection and Instalation of Pressure Relieving Devices in Refineries. Part I. Sizing and Selection

(15) Norma API 521.
Guide for pressure relieving and depressuring systems

(16) Norma API 2000.
Venting atmospheric and low pressure storage tanks

(17) Norma BS 2915: 1990.
Specification for Bursting Discs and Bursting Disc Devices

(18) Norma ADMerkblatt A1.
Safety devices against excess pressure. Bursting Safety Devices. Berlin