Comment Déterminer Ki Pour Les Gaz D'alimentation Supercritiques ? Maîtriser Les Calculs De Vaporisation Des Réacteurs

Pour un gaz d'alimentation supercritique, le facteur de vaporisation (Kᵢ) n'est pas trouvé dans les tables de pression de vapeur—il est déterminé en utilisant un cadre thermodynamique non symétrique. La relation gouvernante est Kᵢ = (Hᵢ · γᵢ) / (φᵢ · P)*, où Hᵢ est la constante de Henry du gaz dans le solvant, γᵢ* adopte la convention non symétrique (γᵢ* → 1 quand xᵢ → 0), φᵢ est le coefficient de fugacité en phase vapeur, et P est la pression totale. Cette expression remplace la valeur K conventionnelle de la loi de Raoult car la pression de vapeur du gaz pur n'est pas définie au-dessus de la température critique.

Déterminer Kᵢ pour des alimentations supercritiques comme l'hydrogène ou l'oxygène repose entièrement sur la loi de Henry avec un coefficient d'activité non symétrique. La valeur Kᵢ = yᵢ/xᵢ devient (Hᵢ · γᵢ*) / (φᵢ · P), liant la solubilité, les interactions moléculaires et la non-idéalité de la phase gazeuse—une conséquence directe du traitement du gaz dissous comme un soluté qui ne se condense jamais en un liquide pur.

Les Fondements Thermodynamiques pour les Gaz Supercritiques

Pourquoi l'Équilibre Vapeur-Liquide Traditionnel Échoue

Le facteur K classique, yᵢ/xᵢ = Pᵢˢᵃᵗ/P, repose sur une pression de vapeur de composant pur.
Pour un composant dont la température dépasse sa température critique (par exemple, l'hydrogène à 298 K), Pᵢˢᵃᵗ n'existe pas. Le gaz ne peut pas se liquéfier par la seule pression, donc la convention symétrique standard devient dénuée de sens.

Les ingénieurs concevant des réacteurs gaz-liquide à l'échelle pilote doivent donc passer à un cadre où le composant supercritique est traité comme un soluté dissous, et non comme une vapeur condensable. La condition d'équilibre fᵢⱽ = fᵢᴸ tient toujours, mais la fugacité en phase liquide doit référencer la loi de Henry plutôt que la fugacité du liquide pur.

La Convention Non Symétrique et la Loi de Henry

La fugacité en phase liquide pour un soluté supercritique s'écrit fᵢᴸ = γᵢ · Hᵢ · xᵢ.*
Ici Hᵢ est la constante de Henry (une propriété de la paire soluté-solvant à la température du système), et γᵢ* est un coefficient d'activité normalisé de sorte que γᵢ* → 1 lorsque la fraction molaire xᵢ approche zéro. Cette convention non symétrique ancre l'état standard à la dilution infinie, exactement là où la loi de Henry gouverne.

Hᵢ capture le comportement de solubilité limite—plus l'interaction entre le gaz et le liquide est forte, plus Hᵢ est petit et plus la concentration dissoute est élevée. γᵢ* corrige ensuite les écarts par rapport au comportement de dilution idéal lorsque xᵢ n'est pas négligeable.

L'Expression Complète pour Kᵢ

Égaliser les fugacités des phases vapeur et liquide donne φᵢ · yᵢ · P = γᵢ · Hᵢ · xᵢ.*
En réarrangeant, on obtient directement le facteur de vaporisation :

[ K_i = \frac{y_i}{x_i} = \frac{H_i \cdot \gamma_i^*}{\phi_i \cdot P} ]

Chaque terme est accessible : Hᵢ provient de corrélations de la littérature ou d'un ensemble limité de mesures de solubilité en unité pilote ; γᵢ* des modèles d'énergie de Gibbs excédentaire (ex : NRTL, UNIQUAC) adaptés à la convention non symétrique ; φᵢ d'une équation d'état cubique (Peng‑Robinson, Soave‑Redlich‑Kwong) ; et P est simplement la pression de fonctionnement du réacteur. Cette expression est le remplacement logique, et nécessaire, de l'analogue de la loi de Raoult chaque fois qu'un flux d'alimentation entre au-dessus de sa température critique.

Application de l'Équation dans un Contexte d'Unité Pilote

Obtention de la Constante de Henry

Hᵢ est la pierre angulaire du calcul.
Elle est fortement dépendante de la température et, pour des pressions modérées à élevées, doit être corrigée pour la pression via l'équation de Krichevsky‑Kasarnovsky [ \ln H_i(P) = \ln H_i(P_{\text{sat}}) + \frac{\bar{v}i^\infty}{RT}(P - P{\text{sat}}) ], où (\bar{v}_i^\infty) est le volume molaire partiel du soluté à dilution infinie.

Pour de nombreux gaz d'unité pilote (H₂, O₂, CO) dans des solvants organiques courants, les valeurs fiables de Hᵢ sont tabulées ou peuvent être prédites par des méthodes de contribution de groupes. En l'absence de données, quelques expériences de solubilité à l'échelle laboratoire établissent rapidement la constante de Henry dans les conditions de température du réacteur.

Estimation du Coefficient d'Activité (γᵢ*)

Même dans des solutions relativement diluées, γᵢ peut s'écarter de l'unité si le solvant est complexe ou si plusieurs solutés interagissent.*
Lorsque la fraction molaire en phase liquide du gaz supercritique reste inférieure à environ 0,001, l'hypothèse γᵢ* = 1 introduit souvent une erreur inférieure à 5 %—acceptable pour le dimensionnement initial d'une unité pilote.

Pour des systèmes plus concentrés, des modèles comme UNIQUAC ou NRTL peuvent être re-paramétrés pour la convention non symétrique. Les paramètres d'interaction binaire nécessaires sont généralement régressés à partir de données de solubilité des gaz ou estimés par similarité moléculaire. L'arbre de décision de l'ingénieur d'unité pilote se divise donc ainsi : utiliser γᵢ* = 1 pour une évaluation rapide, ou investir dans un modèle d'activité complet lorsque la prédiction précise du transfert de matière est critique.

Calcul du Coefficient de Fugacité en Phase Vapeur (φᵢ)

Aux pressions modérées courantes dans de nombreux réacteurs pilotes (5–20 bar), φᵢ diffère rarement de 1 de plus de quelques pourcents.
Cependant, les hydrogénations ou oxydations à haute pression fonctionnent souvent à 50–200 bar, où la non-idéalité de la phase gazeuse devient significative. Une équation d'état cubique, adaptée au mélange gazeux et à la température, fournit φᵢ directement.

Dans le cadre Peng‑Robinson, par exemple, φᵢ est une fonction de la température réduite, de la pression et du facteur acentrique. Pour l'hydrogène pur, la correction φᵢ peut augmenter Kᵢ de 10 à 30 % à 100 bar par rapport à l'hypothèse du gaz parfait ; l'ignorer conduirait donc à sous-estimer la fraction molaire en phase gazeuse et à fausser les calculs de force motrice de transfert de matière.

Le Rôle de la Pression

La pression apparaît explicitement au dénominateur de l'expression de Kᵢ, rendant Kᵢ inversement proportionnel à P lorsque tous les autres termes restent constants.
Cela signifie que lorsque la pression du réacteur augmente, la fraction molaire à l'équilibre en phase vapeur pour une concentration donnée en phase liquide diminue—un reflet direct de la loi de Henry : plus de gaz se dissout, donc moins reste dans l'espace de tête ou la phase bulle.

Dans la conception d'unités pilotes, cet effet de pression doit être couplé avec la dépendance à la pression souvent négligée de Hᵢ lui-même. À hautes pressions, Hᵢ augmente, contrecarrant partiellement l'effet de P. Un modèle complet de Kᵢ couple donc la correction de Krichevsky‑Kasarnovsky avec le dénominateur, garantissant que le facteur de vaporisation reflète avec précision l'amélioration de la solubilité dans les conditions de fonctionnement réelles.

Comprendre les Compromis et les Limitations

L'Hypothèse de Solution Diluée

Traiter γᵢ* = 1 est tentant mais pas toujours sûr. Dans des mélanges d'hydrogénation très concentrés ou lorsque les fractions de gaz dissous dépassent quelques pourcents molaires, le coefficient d'activité peut s'écarter de l'unité de 10 à 20 %. L'erreur résultante sur Kᵢ se propage dans les estimations de concentration en phase liquide et finalement dans les prédictions de volume de réacteur. Pour une faisabilité préliminaire, l'erreur est souvent tolérable ; pour une mise à l'échelle finale, elle ne l'est pas.

Sensibilité à la Température et à la Pression

Les constantes de Henry pour les gaz légers doublent ou triplent généralement sur un changement de température de 50 K, et le coefficient de fugacité évolue avec la pression. Une valeur de Kᵢ valide à un point de fonctionnement peut devenir inexacte si le profil de température du réacteur change—phénomène courant lors des démarrages ou des épisodes de désactivation du catalyseur. L'équipe de l'unité pilote doit traiter Kᵢ comme une variable dynamique, et non comme un nombre statique, lors de la modélisation du réacteur sur l'ensemble de son enveloppe de fonctionnement.

Mélanges Multicomposants

Lorsque l'alimentation en gaz contient plusieurs composants supercritiques (par exemple, H₂ et CO dans le gaz de synthèse), le coefficient de fugacité φᵢ devient une propriété du mélange, et les paramètres d'interaction croisée dans l'équation d'état sont importants. Utiliser des estimations de φᵢ pour des composants purs peut introduire un biais systématique. Dans la mesure du possible, les ingénieurs d'unité pilote devraient valider les prédictions de Kᵢ contre quelques mesures de pression totale et de composition vapeur, affinant les règles de mélange de l'équation d'état si nécessaire.

Faire le Bon Choix pour Votre Objectif d'Unité Pilote

Le niveau de détail approprié pour déterminer Kᵢ dépend entièrement de l'usage que vous faites de ce nombre.

Si votre objectif principal est un criblage de faisabilité rapide : Utilisez la forme simplifiée Kᵢ ≈ Hᵢ/P avec des valeurs de Hᵢ de la littérature et supposez γᵢ* = φᵢ = 1. Cela donne des estimations de l'ordre de grandeur pour la séparation gaz-liquide et les débits minimaux de solvant en un après-midi.
Si votre objectif principal est une modélisation précise du transfert de matière : Employez l'expression complète avec un Hᵢ corrigé en pression (via la relation de Krichevsky‑Kasarnovsky), une équation d'état de pointe pour φᵢ, et au moins un modèle d'activité de premier ordre pour γᵢ*. Vos simulations de colonnes à bulles ou de réacteurs agités vous récompenseront avec des valeurs KLa qui résistent à l'examen de la mise à l'échelle.
Si votre objectif principal est le fonctionnement à haute pression (au-dessus de 30 bar) : Ne négligez jamais la correction du coefficient de fugacité. Une équation d'état cubique—ajustée aux données binaires publiées ou aux points de rosée mesurés—est l'exigence minimale pour éviter une sous-estimation systématique des concentrations en phase gazeuse.
Si votre objectif principal est le criblage de solvants ou les alimentations multicomposants : Investissez dans des essais de solubilité expérimentaux pour régresser les paramètres Hᵢ et γᵢ*. L'effort se rentabilise en réduisant l'incertitude du dimensionnement du réacteur et en évitant les surprises dans les unités de séparation en aval.

La détermination précise de Kᵢ pour les alimentations supercritiques n'est pas un exercice théorique—c'est la pièce maîtresse qui relie la solubilité, le transfert de matière et le dimensionnement du réacteur pour aboutir à une conception d'unité pilote fiable.

Tableau Récapitulatif :

Paramètre	Symbole	Rôle dans l'équation $K_i$	Méthode d'Estimation Clé
Facteur de Vaporisation	$K_i$	Mesure l'équilibre gaz-liquide ($y_i/x_i$)	Calculé comme $(H_i \cdot \gamma_i^*) / (\phi_i \cdot P)$
Constante de Henry	$H_i$	Définit le comportement de solubilité limite	Correction par l'équation de Krichevsky-Kasarnovsky
Coefficient d'Activité	$\gamma_i^*$	Corrige la non-idéalité de la phase liquide	NRTL / UNIQUAC (convention non symétrique)
Coefficient de Fugacité	$\phi_i$	Corrige la non-idéalité de la phase vapeur	Équation d'État Cubique (ex : Peng-Robinson)
Pression du Système	$P$	Représente la pression de fonctionnement du réacteur	Dénominateur explicite ; inversement proportionnel à $K_i$

L'optimisation des modèles thermodynamiques des réacteurs gaz-liquide nécessite des tests précis à l'échelle pilote. LABPARK fournit des Unités Pilotes Opérationnelles Pédagogiques et Professionnelles de premier plan en génie chimique, bioprocédés & biotechnologie, et traitement de l'eau & environnement pour les universités, instituts de recherche et entreprises. Renforcez vos recherches et votre formation avec des systèmes prêts pour la mise à l'échelle—contactez LABPARK dès aujourd'hui pour discuter de votre projet !