Passer au contenu principal

k = exposant isoentropique

L'importance de  k  pour soupape de sécurité

édité par Alessandro Ruzza 

Le dimensionnement des soupapes de sûreté destinées à évacuer des gaz ou des vapeurs, selon lspesl Collection « E », nécessite la connaissance de l'exposant isoentropique k aux conditions d'évacuation.

Une application négligente du chapitre « E.1 » de la Collection « E » de lspesl, concernant le dimensionnement des soupapes de sûreté, peut conduire à une surestimation de la capacité de décharge des soupapes et des disques de rupture.

Cet article donne quelques lignes directrices pour estimer la valeur de k pour les gaz réels et
met en évidence l'erreur en considérant k égal au rapport des chaleurs spécifiques Cp/Cv

Une première et grossière erreur à éviter est d'utiliser la formule de la Collection 'E', valable pour les gaz ou les vapeurs, dans les situations où un décharge biphasée de liquide et de gaz/vapeur a lieu. Dans de tels cas, en effet, les diamètres calculés seront sans doute sous-dimensionnés par rapport au besoin réel.
Une deuxième erreur, qui dans de nombreux cas peut conduire à la sous-dimensionner le système de sécurité, est de donner à l'exposant isoentropique k la valeur du rapport Cp/Cv. Alors que le premier point fera l'objet d'une série d'articles ultérieurs, nous souhaitons ici donner quelques indications utiles pour calculer l'exposant isoentropique et montrer, dans des cas concrets, l'ampleur de l'erreur qui peut être commise.

Écoulement isoentropique à travers une buse

 

La formule qui est utilisé dans la collection "E", ainsi que dans d'autres italiens et étranger standards, pour le calcul des soupapes de sécurité qui doivent évacuer des gaz ou des vapeurs, est celui du débit isoentropique à travers une tuyère dans des conditions de saut critique, qui pour un gaz parfait est :

Formule lspesl Collection « E »

où l'expansisur le coefficient C est donnée par :

expansisur le coefficient C

va k l'exposant de l'isoentropique expansisur l'équation : pxv^k=coût

FluideP1 (bar)T1 (°C)q' (kg/h)q (kg/heure)(q'/q) x 100
Méthane125014721466100.4
Méthane2320023142267102.1
Propane1210022612181103.7
hexane1217830992740113.1
hexane2322065195111127.5
heptane1221532322821114.4

q'= débit calculé avec k = Cp/Cv (20 °C, 1 atm)
q = débit calculé avec k = (Cp/Cv) • (Z/Zp)

En introduisant le coefficient expérimental k de débit de la soupape de sécurité, qui prend globalement en compte les performances de débit réelles de la soupape, un coefficient de sécurité de 0.9 et le facteur de compressibilité Z1 pour le fluide réel, on arrive à la formulation de la collection « E » :

(1)

L'exposant isoentropique k peut s'exprimer comme suit:

Pour une gaz parfait, Pour qui P x V / R x T =1 , il est démontré que k est égal au rapport Cp/Cv entre les chaleurs spécifiques à pression et volume constants.

Pour une vrai gaz, k peut être exprimée (voir annexe B) ​​par :

où Z est le facteur de compressibilité défini par Z=P x V / R x T et Zp est le "facteur de compressibilité dérivé". Lors de l'application de la formule , selon la collection "E", les valeurs de Cp/Cv, Z et Zp doivent être évaluées aux conditions de décharge P1 et T1.

Le facteur de compressibilité dérivé Zp est défini dans la formule comme:

Le facteur de compressibilité Z peut être exprimé par :

et de même, peut être exprimé comme:

où les valeurs de Z^0, Z^1, Zp^0, Zp^1 sont tabulées à l'annexe A en fonction de Pr et Tr.

In et , Ω est le facteur acentrique de Pitzer défini par :

Où Pr^SAT est la pression de vapeur réduite correspondant à une valeur de température réduite Tr=T/Tc=0,7. L'annexe A montre les valeurs Ω de certains fluides. Z e Zp peut aussi être déduit directement d'une équation d'état analytique.

Un exemple chiffré

 

En prenant un exemple numérique, supposons que nous devions calculer la capacité de décharge d'une soupape de sécurité dans les conditions suivantes :

Fluiden-Butano
État physiquevapeur surchauffée
Masse moléculaireM58,119
Pression de réglageP19,78 bar
Surpression10%
Température du fluideT400 K
Coefficient d'efflux0,9
Diamètre de l'orificeDo100 mm

la pression de refoulement est donnée par :

étant pour le n-Butane : Tc=425,18 K et Pc=37,96 bar, on a:

et en utilisant les tableaux de l'annexe A, nous avons :

Connaissant le volume spécifique de la vapeur aux conditions de décharge (P1, T1) égal à 0,01634 m^3/kg (0,0009498 m^3/g-mole), on aurait également pu calculer Z à partir de :

Étant donné le rapport des chaleurs spécifiques à pression et volume constants, aux conditions de décharge (P1T1), égal à 1,36, de la formule nous avons:

147060

Application de la formule [1], avec calcul du débit

Appliquer la formule , qui a été résolu pour le calcul du débit, nous avons une valeur de débit de décharge de 147.060 kg / h.

174848

En appliquant la formule [1], en utilisant la valeur de Cp/Cv à 1 atm et 20 °C

Si nous avions plutôt utilisé la valeur de Cp/Cv à 1 atm et 20 °C, nous aurions eu k = 1,19 et de la formule un débit de décharge de 174.848 kg / h.

Cela nous aurait conduit à surestimer la décharge capacité de la soupape de sécurité d'environ 19%

AVERTISSEMENT:

L'erreur que l'on peut commettre en attribuant la valeur Cp/Cv à k peut être beaucoup plus élevée que dans cet exemple.

PLUS DE 20%

A titre indicatif, le tableau suivant donne les débits d'un orifice de 18 mm pour les autres hydrocarbures saturés, calculés dans les deux cas. Les calculs ont été effectués avec spécialement developed logiciel.

FluideP1 (bar)T1 (°C)q' (kg/h)q (kg/heure)(q'/q) x 100
Méthane125014721466100.4
Méthane2320023142267102.1
Propane1210022612181103.7
hexane1217830992740113.1
hexane2322065195111127.5
heptane1221532322821114.4

Le logiciel n'utilise pas de formules  mais, à partir de la version modifiée Équation d'état de Redlich et Kwong, calcule la valeur de l'exposant isoentropique à l'aide de corrélations thermodynamiques.

Annexes A et B
dérivation de formules

BESA sera présent au IVS - IVS Industrial Valve Summit 2024