Production d'Acide Sulfurique (Procédé de Contact)
📝 Situation du Projet et Enjeux Industriels
Vous avez récemment intégré le bureau d'études "Process Engineering" d'un vaste complexe pétrochimique stratégique, situé au cœur de la "Vallée de la Chimie". Cette installation majeure produit annuellement plus de 500 000 tonnes d'acide sulfurique (\(H_2SO_4\)) concentré à 98% (qualité technique), un intrant absolument critique pour l'ensemble de l'industrie chimique environnante, notamment pour la fabrication d'engrais phosphatés, le décapage des métaux et le traitement des minerais. L'unité de production actuelle repose sur le standard industriel éprouvé : le Procédé de Contact Double Absorption (DCDA).
Récemment, la direction technique a identifié un besoin urgent de modernisation. Face au durcissement drastique des normes environnementales européennes (IED/BREF) concernant les rejets atmosphériques de dioxyde de soufre (\(SO_2\)) et dans une optique de maximisation de l'efficacité énergétique (récupération de vapeur haute pression), une refonte de la section réactionnelle est à l'étude. Le projet se concentre spécifiquement sur le remplacement et l'optimisation du catalyseur au pentoxyde de vanadium (\(V_2O_5\)) dopé aux sulfates alcalins (K, Cs) utilisé dans le convertisseur principal à quatre lits.
Votre mission est de haute technicité : modéliser et valider le comportement thermodynamique de l'étape la plus critique du procédé, à savoir l'oxydation catalytique du dioxyde de soufre (\(SO_2\)) en trioxyde de soufre (\(SO_3\)). Cette réaction est le véritable "cœur battant" de l'usine. Contrairement à une simple combustion, il s'agit d'un équilibre chimique délicat. Vous devez déterminer les fenêtres opératoires (températures d'entrée et de sortie des lits) qui permettront de garantir un taux de conversion global supérieur à 99,7%, seuil indispensable pour respecter les limites d'émission sans surcharger le système de lavage des gaz de queue (scrubber). De plus, la gestion de l'exothermicité massive de cette réaction est vitale : une dérive thermique pourrait entraîner la désactivation irréversible du catalyseur par frittage ou la corrosion accélérée des équipements.
En votre qualité d'Ingénieur Procédés Chimie confirmé, vous avez la responsabilité de dimensionner les paramètres thermodynamiques clés du premier étage du convertisseur catalytique. Votre note de calcul servira de base à la spécification des échangeurs de chaleur intermédiaires. Vous devrez calculer avec précision l'enthalpie de réaction standard, déterminer l'évolution de la constante d'équilibre en fonction de la température via la relation de Van't Hoff intégrée, et valider le taux de conversion théorique maximal atteignable pour assurer la sécurité thermique (risque d'emballement) et la viabilité économique de l'installation.
"Attention, point de vigilance critique : la réaction d'oxydation du SO2 est extrêmement exothermique. Le profil de température doit être surveillé de près. Si la température locale au sein du lit catalytique dépasse les 620°C (893 K), nous risquons une désactivation irréversible du catalyseur (perte de surface active) et, paradoxalement, une chute drastique du rendement à l'équilibre conformément au principe de Le Chatelier. Soyez d'une rigueur absolue sur le calcul du Kp et du bilan thermique."
Afin de mener à bien l'étude thermodynamique du réacteur, nous nous appuyons sur un ensemble de données validées et de normes internationales. Les calculs seront effectués en considérant le comportement des fluides comme celui des gaz parfaits. Cette hypothèse est justifiée par les conditions opératoires (Pression proche de l'atmosphérique et températures élevées, éloignant les gaz de leur point de liquéfaction).
📚 Référentiel Normatif & Standards
Le dimensionnement s'inscrit dans le cadre strict des réglementations suivantes :
- IUPAC Green Book (3e éd.) : Convention de signes et notations thermodynamiques standard (Enthalpie de formation \(\Delta_{\text{f}} H^\circ\), Entropie \(S^\circ\)).
- Directive Seveso III (2012/18/UE) : Gestion des risques liés aux substances dangereuses. L'exothermicité de la réaction impose une analyse de sécurité thermique (emballement).
- Tables JANAF (NIST) : Source de référence pour les propriétés thermochimiques des espèces gazeuses.
Les valeurs ci-dessous correspondent à l'état standard de référence (\(P^\circ = 1\) bar, \(T = 298,15\) K). Elles sont tirées des tables JANAF et serviront de base pour l'application des lois de Kirchhoff et de Hess.
| Espèces Chimiques (État Gazeux) | ||
| Composé | Enthalpie Standard de Formation \(\Delta_{\text{f}} H^\circ\) (kJ/mol) | Entropie Standard Molaire \(S^\circ\) (J/mol·K) |
|---|---|---|
| Dioxyde de Soufre \(\text{SO}_2(g)\) | -296,8 | 248,2 |
| Dioxygène \(\text{O}_2(g)\) | 0 (État standard de l'élément) | 205,2 |
| Trioxyde de Soufre \(\text{SO}_3(g)\) | -395,7 | 256,8 |
📐 Conditions Opératoires (Entrée du Lit n°1)
Le gaz d'alimentation provient directement de l'unité de combustion du soufre liquide, après refroidissement, filtration des poussières et séchage strict à l'acide sulfurique concentré (pour éliminer toute trace d'eau corrosive). Ses caractéristiques à l'entrée du convertisseur sont les suivantes :
- Température d'entrée (\(T_e\)): 700 K (427°C) - Température d'amorçage du catalyseur.
- Pression Totale (\(P_{\text{tot}}\)): 1,0 bar (soit \(10^5\) Pa) - Légère surpression pour vaincre les pertes de charge.
- Composition du gaz (volumique):
- \(\text{SO}_2\): 10% (Limité pour contrôler la chauffe)
- \(\text{O}_2\): 11% (Excès stœchiométrique pour favoriser l'équilibre)
- \(\text{N}_2\): 79% (Inerte, gaz vecteur thermique)
- Constante des gaz parfaits (\(R\)): 8,314 J/mol·K
⚖️ Équation de Réaction Étudiée
Il s'agit d'une oxydation catalytique hétérogène en phase gazeuse :
| Variable | Symbole | Type | Unité |
|---|---|---|---|
| Avancement réactionnel | \(\xi\) (xi) | Variable d'état extensive | mol |
| Enthalpie libre de réaction standard | \(\Delta_{\text{r}} G^\circ(T)\) | Potentiel Thermodynamique | kJ/mol |
| Constante d'équilibre thermodynamique | \(K_p\) (ou \(K^\circ\)) | Grandeur adimensionnelle | - |
E. Protocole de Résolution
Pour mener à bien ce dimensionnement thermodynamique, nous suivrons une méthodologie rigoureuse, allant des propriétés intrinsèques des molécules jusqu'au comportement macroscopique du réacteur.
Bilan Enthalpique Standard
Détermination de la chaleur de réaction à 298K (Loi de Hess) pour quantifier l'exothermicité du système.
Équilibre & Loi d'Action de Masse
Calcul de l'entropie, de l'enthalpie libre de Gibbs \(\Delta_{\text{r}} G^\circ(T)\) et déduction de la constante d'équilibre \(K_p\) à la température de fonctionnement (700K).
Conversion à l'Équilibre
Établissement du tableau d'avancement et résolution de l'équation d'équilibre pour trouver le taux de conversion maximal théorique.
Bilan Thermique & Sécurité
Vérification de l'énergie libérée par rapport à la capacité thermique des flux pour éviter l'emballement thermique du réacteur.
Production d'Acide Sulfurique (Procédé de Contact)
🎯 Objectif Scientifique
L'objectif fondamental de cette première étape de dimensionnement est de quantifier avec précision l'énergie thermique libérée (ou absorbée) par la réaction d'oxydation du dioxyde de soufre. Dans un contexte industriel, la connaissance exacte de cette enthalpie de réaction (\(\Delta_{\text{r}} H^\circ\)) est le paramètre zéro du design thermique. Elle conditionne non seulement la taille des échangeurs de chaleur mais surtout la stratégie de sécurité du procédé : une réaction fortement exothermique présente des risques d'emballement thermique ("runaway") pouvant conduire à la destruction du catalyseur ou à la fusion du réacteur. Nous allons donc déterminer la signature thermique de la réaction dans les conditions standard.
📚 Référentiel Théorique
1er Principe ThermodynamiqueLoi de HessÉtats StandardPour calculer l'enthalpie d'une réaction globale, il est inutile et coûteux de procéder à des mesures calorimétriques directes pour chaque nouveau procédé. Nous exploitons la propriété de l'enthalpie \(H\) d'être une fonction d'état : sa variation ne dépend que de l'état final et de l'état initial, et non du chemin chimique parcouru. C'est le fondement de la Loi de Hess. Nous allons décomposer notre réaction en une somme algébrique de réactions de formation à partir des corps simples, pour lesquelles les valeurs sont tabulées (JANAF). La rigueur impose de bien respecter la stœchiométrie : l'équation bilan implique 2 moles de \(\text{SO}_2\), le résultat sera donc associé à cet avancement spécifique.
La loi de Hess stipule que l'enthalpie standard de réaction \(\Delta_{\text{r}} H^\circ\) est égale à la somme des enthalpies standard de formation des produits pondérée par leurs coefficients stœchiométriques, moins la somme correspondante pour les réactifs. Rappel crucial : par convention, l'enthalpie standard de formation d'un corps pur simple dans sa phase la plus stable (ici \(\text{O}_2\) gazeux) est strictement nulle.
Étape 1 : Données d'Entrée (Extrait Tables JANAF à 298K)
| Composé Chimique | Enthalpie de Formation \(\Delta_{\text{f}} H^\circ\) (kJ/mol) |
|---|---|
| Trioxyde de soufre \(\text{SO}_3(g)\) | -395,7 |
| Dioxyde de soufre \(\text{SO}_2(g)\) | -296,8 |
| Dioxygène \(\text{O}_2(g)\) | 0,0 |
Vérifiez systématiquement l'état physique des constituants dans vos tables. L'enthalpie de formation de \(\text{SO}_3\) gaz est différente de celle de \(\text{SO}_3\) liquide ou solide. Ici, à haute température dans le réacteur, tout est en phase gazeuse.
Calcul Détaillé de l'Enthalpie Standard de Réaction
1. Décomposition Algébrique des Produits :
L'équation bilan produit 2 moles de \(\text{SO}_3\). Nous appliquons donc le coefficient stœchiométrique \(\nu = 2\) à l'enthalpie de formation molaire du \(\text{SO}_3\).
Interprétation : L'état final du système possède une énergie de formation globale de -791,4 kJ.
2. Décomposition Algébrique des Réactifs :
L'équation consomme 2 moles de \(\text{SO}_2\) et 1 mole de \(\text{O}_2\). Nous sommons leurs enthalpies respectives pondérées par ces coefficients.
Interprétation : L'état initial du système possède une énergie de formation globale de -593,6 kJ.
3. Soustraction Finale (Loi de Hess) :
Nous appliquons maintenant la formule de Hess : \(H_{\text{produits}} - H_{\text{réactifs}}\). Notez comment le signe moins de la formule interagit avec le signe négatif de l'enthalpie des réactifs.
Interprétation : Le résultat final est négatif (-197,8 kJ). Cela signifie que le système perd de l'énergie thermique au profit de son environnement. La réaction est donc exothermique.
La valeur obtenue de -197,8 kJ/mol est fondamentale. Le signe négatif indique sans ambiguïté que la réaction est exothermique. Cela signifie que le réacteur va naturellement s'échauffer au fur et à mesure que la conversion progresse. Dans l'industrie, cette chaleur n'est pas perdue : elle est récupérée pour produire de la vapeur haute pression, qui alimente ensuite des turbines pour générer de l'électricité. L'unité de production d'acide sulfurique est ainsi souvent excédentaire en énergie.
La valeur de -197,8 kJ correspond à environ 99 kJ par mole de \(\text{SO}_2\) transformée. C'est une densité énergétique élevée, comparable à certaines réactions de combustion modérées. Cet ordre de grandeur est cohérent avec la forte stabilité de la molécule de \(\text{SO}_3\) par rapport au \(\text{SO}_2\). Si nous avions trouvé une valeur positive, cela aurait contredit l'expérience industrielle qui nécessite un refroidissement constant des réacteurs.
Attention à l'unité : le résultat est donné en "kJ par mole d'avancement". Comme l'équation est écrite pour 2 moles de \(\text{SO}_2\), cette énergie correspond à la transformation de 2 moles. Il ne faut pas confondre cette valeur avec l'enthalpie molaire de combustion du soufre. De plus, cette valeur est calculée à 298 K. Bien que l'enthalpie varie légèrement avec la température (loi de Kirchhoff), pour une première approximation d'ingénierie, on considère souvent \(\Delta H\) comme constant, mais une étude fine nécessiterait d'intégrer les capacités calorifiques \(C_p\).
🎯 Objectif Scientifique
La thermodynamique nous impose des limites. Même avec un catalyseur parfait et un temps infini, une réaction réversible ne consommera jamais 100% des réactifs. Elle s'arrête à un état d'équilibre dynamique. Notre objectif est de calculer la position de cet équilibre à la température de fonctionnement industriel (700 K). Nous devons déterminer la constante d'équilibre \(K_p\) pour savoir si la réaction est thermodynamiquement favorisée (K >> 1) ou non à cette température.
📚 Référentiel Théorique
Enthalpie Libre de GibbsRelation de Van't HoffApproximation d'EllinghamLa constante d'équilibre est directement liée à l'enthalpie libre standard de réaction \(\Delta_{\text{r}} G^\circ(T)\). Cependant, les tables ne donnent les valeurs qu'à 298 K. Pour projeter ces valeurs à 700 K, nous allons utiliser l'approximation d'Ellingham. Elle suppose que l'enthalpie standard (\(\Delta_{\text{r}} H^\circ\)) et l'entropie standard (\(\Delta_{\text{r}} S^\circ\)) varient peu avec la température dans la plage considérée. C'est une hypothèse robuste tant qu'il n'y a pas de changement de phase des espèces.
La relation fondamentale liant la thermodynamique à l'équilibre chimique est : \(\Delta_{\text{r}} G^\circ(T) = -R \cdot T \cdot \ln(K)\). Si \(\Delta G\) est très négatif, K est très grand (réaction quasi-totale). Si \(\Delta G\) est positif, K est petit (réaction limitée).
📋 Données d'Entrée (Q1 + Tables)
| Paramètre | Valeur |
|---|---|
| \(\Delta_{\text{r}} H^\circ\) (Q1) | \(-197\,800 \text{ J/mol}\) |
| \(S^\circ(\text{SO}_3)\) | \(256,8 \text{ J/mol}\cdot\text{K}\) |
| \(S^\circ(\text{SO}_2)\) | \(248,2 \text{ J/mol}\cdot\text{K}\) |
| \(S^\circ(\text{O}_2)\) | \(205,2 \text{ J/mol}\cdot\text{K}\) |
| Température (T) | \(700 \text{ K}\) |
Étape 1 : Calcul de l'Entropie Standard de réaction \(\Delta_{\text{r}} S^\circ_{298}\)
Nous devons d'abord déterminer comment l'ordre du système change. Pour cela, nous appliquons la loi de Hess aux entropies molaires absolues (S°). Nous multiplions l'entropie de chaque molécule par son coefficient stœchiométrique.
Interprétation : L'entropie diminue (-188 J/K). C'est logique car nous passons de 3 moles de gaz (désordre élevé) à 2 moles de gaz (ordre plus élevé). Le "coût entropique" de la réaction est défavorable.
Attention ! L'entropie est donnée en Joules (J), alors que l'enthalpie calculée à la Q1 est en Kilojoules (kJ). Avant de les combiner dans la formule de Gibbs, vous devrez impérativement convertir les kJ en J (multiplier par 1000) ou les J en kJ (diviser par 1000). Oublier cette conversion fausse le résultat d'un facteur 1000.
Étape 2 : Calcul de l'Enthalpie Libre \(\Delta_{\text{r}} G^\circ\) à 700 K
Nous utilisons maintenant la température de fonctionnement (700 K) pour calculer l'énergie libre disponible. Nous convertissons ici l'enthalpie \(\Delta H\) (-197,8 kJ) en Joules (-197 800 J) pour l'homogénéité.
Interprétation : \(\Delta G\) reste négatif à 700K (-66,2 kJ/mol). La réaction est donc toujours thermodynamiquement spontanée à cette température, bien que moins favorisée qu'à 298K.
Étape 3 : Déduction de la Constante d'Équilibre \(K_p\)
Nous inversons la relation logarithmique fondamentale pour isoler Kp.
Calcul de l'argument de l'exponentielle :
Calcul final de \(K_p\) :
Enfin, nous prenons l'exponentielle \(e^{11,375}\) pour obtenir la valeur adimensionnelle de la constante.
Interprétation : La constante est très grande (> 1). L'équilibre est fortement déplacé vers la formation de produits.
Le calcul montre qu'à la température de démarrage du catalyseur (700 K), la constante d'équilibre est encore très élevée (\(8,7 \times 10^4\)). Cela signifie que la limite thermodynamique est très repoussée : il est théoriquement possible de convertir la quasi-totalité du \(\text{SO}_2\). Le choix de cette température est donc validé sur le plan thermodynamique : elle permet un excellent rendement potentiel.
La valeur de \(K_p\) diminue exponentiellement avec la température pour les réactions exothermiques (Loi de Van't Hoff). Si nous avions effectué le calcul à 1000 K, le terme entropique \(-T\Delta S\) aurait dominé, rendant \(\Delta G\) positif et \(K_p\) inférieur à 1. Le résultat obtenu est donc parfaitement cohérent avec la théorie : "Plus il fait chaud, moins on produit", mais à 700 K, nous sommes encore dans une zone très favorable.
Le principal piège ici est de croire que parce que \(K_p\) est grand, la réaction sera complète et rapide. La thermodynamique ne dit rien sur la vitesse. À 700 K, sans catalyseur, la réaction est infiniment lente. C'est la présence du \(V_2O_5\) qui permet d'atteindre cet équilibre en un temps raisonnable. De plus, l'approximation d'Ellingham néglige les variations de \(C_p\) avec la température, ce qui introduit une petite erreur (acceptable en pré-dimensionnement) sur la valeur exacte de \(K_p\).
🎯 Objectif Scientifique
Le calcul de \(K_p\) est une étape intermédiaire abstraite. L'ingénieur doit traduire cette constante abstraite en une grandeur concrète et mesurable : le taux de conversion (\(\alpha\)). Quel pourcentage exact des molécules de \(\text{SO}_2\) entrantes ressortira sous forme de \(\text{SO}_3\) ? Nous allons relier la constante \(K_p\) aux pressions partielles des gaz via un bilan matière rigoureux (tableau d'avancement).
📚 Référentiel Théorique
Loi de DaltonLoi d'Action de MasseNous partons d'un mélange complexe (\(\text{SO}_2, \text{O}_2, \text{N}_2\)). L'azote est inerte mais il dilue le mélange, ce qui affecte les pressions partielles. Nous allons dresser un bilan molaire en fonction de l'avancement \(\alpha\), exprimer la fraction molaire de chaque gaz, multiplier par la Pression Totale pour avoir la Pression Partielle, et injecter le tout dans l'expression de K.
Étape 1 : Construction du Tableau d'Avancement (Bilan Matière)
Pour simplifier les calculs de fractions, nous raisonnons sur une base de 100 moles de gaz entrant dans le réacteur. Nous définissons \(\alpha\) comme le taux de conversion du réactif limitant (\(\text{SO}_2\)), tel que \(0 \le \alpha \le 1\).
1.1. Détermination des Variations Molaires (Stœchiométrie)
L'équation de réaction est \(2\text{SO}_2 + \text{O}_2 \rightleftharpoons 2\text{SO}_3\). Cela signifie que pour 2 moles de \(\text{SO}_2\) consommées, 1 mole d'\(\text{O}_2\) disparaît et 2 moles de \(\text{SO}_3\) apparaissent.
Interprétation : L'oxygène est consommé deux fois moins vite que le dioxyde de soufre.
1.2. Calcul du Nombre de Moles Total à l'Équilibre
La quantité de matière totale varie car la réaction modifie le nombre de molécules gazeuses. Nous sommons les variations :
Interprétation : Le nombre total de moles diminue de \(5\alpha\). C'est une contraction de volume.
1.3. Tableau Récapitulatif
| Espèce | État Initial (n0) | Variation | État Équilibre (n_éq) | Fraction Molaire (x_i) |
|---|---|---|---|---|
| \(\text{SO}_2\) | 10 | \(-10\alpha\) | \(10(1-\alpha)\) | \(\frac{10(1-\alpha)}{n_{\text{tot}}}\) |
| \(\text{O}_2\) | 11 | \(-5\alpha\) | \(11 - 5\alpha\) | \(\frac{11 - 5\alpha}{n_{\text{tot}}}\) |
| \(\text{SO}_3\) | 0 | \(+10\alpha\) | \(10\alpha\) | \(\frac{10\alpha}{n_{\text{tot}}}\) |
| \(\text{N}_2\) | 79 | 0 | 79 | \(\frac{79}{n_{\text{tot}}}\) |
| TOTAL | 100 | \(-5\alpha\) | \(n_{\text{tot}} = 100 - 5\alpha\) | 1 |
Remarquez que le nombre total de moles diminue (\(100 \rightarrow 100 - 5\alpha\)). C'est une contraction de volume (3 moles de gaz donnent 2 moles). Selon le principe de Le Chatelier, une augmentation de pression favoriserait donc la réaction. Ici, nous restons à 1 bar.
Étape 2 : Développement Analytique de Kp
Nous partons de la définition de \(K_p\) basée sur les pressions partielles. La pression partielle \(P_i\) est le produit de la fraction molaire \(x_i\) par la pression totale \(P_{\text{tot}}\).
1. Formulation selon la Loi d'Action de Masse :
2. Simplification algébrique (avec \(P_{\text{tot}} = 1\) bar) :
Les termes de pression \(P_{\text{tot}}\) au numérateur et au dénominateur se simplifient. Il reste un facteur \(1/P_{\text{tot}}\) au dénominateur, mais comme \(P_{\text{tot}} = 1\), ce facteur vaut 1. Nous substituons ensuite chaque fraction molaire \(x_i\) par son expression du tableau (moles / \(n_{\text{tot}}\)).
3. Expression Finale en fonction de \(\alpha\) :
En remplaçant \(n_{\text{tot}}\) par \(100 - 5\alpha\), on obtient l'équation finale à résoudre :
Étape 3 : Résolution Numérique
Nous avons une équation où un terme très grand (87 116) est égal à une fraction. Pour que cela soit possible, le dénominateur \((1-\alpha)^2\) doit être extrêmement petit, ce qui signifie que \(\alpha\) est très proche de 1. En utilisant une méthode par itération (en testant 0.99, 0.991, 0.992...), nous trouvons la racine :
Interprétation : Le taux de conversion atteint 99,2%. Presque tout le réactif est consommé.
Le taux de conversion théorique de 99,2% calculé est extrêmement prometteur. Il suggère qu'avec un temps de séjour suffisant et en maintenant la température à 700 K, on peut transformer la quasi-totalité du \(\text{SO}_2\) polluant en \(\text{SO}_3\) valorisable dès le premier passage. Cela valide la faisabilité économique du procédé à cette température.
Ce résultat est cohérent avec la très grande valeur de \(K_p\) trouvée précédemment. Si \(K_p\) est grand, le rapport produits/réactifs doit être grand, ce qui implique mathématiquement que \(\alpha\) doit tendre vers 1. Si nous avions trouvé un \(\alpha\) de 50% avec un \(K_p\) de 87 000, il y aurait eu une erreur de calcul manifeste.
Attention, ce calcul est effectué "à l'équilibre". Dans la réalité industrielle, l'équilibre n'est jamais parfaitement atteint car cela nécessiterait un réacteur infini. De plus, ce calcul suppose que la température reste fixe à 700 K. Or, la réaction chauffe énormément ! Si la température monte, \(K_p\) chute et \(\alpha\) avec. Ce résultat de 99,2% est donc une limite théorique idéale (isotherme) et non la conversion réelle en sortie de réacteur adiabatique.
🎯 Objectif Scientifique
Nous abordons ici la question de sécurité la plus critique. Les réacteurs industriels de grande taille (plusieurs mètres de diamètre) se comportent presque comme des enceintes adiabatiques : les pertes thermiques par les parois sont négligeables face aux flux d'énergie internes. La chaleur dégagée par la réaction n'a nulle part où aller, elle s'accumule dans le gaz. Nous devons calculer la température finale adiabatique. Si celle-ci dépasse la température de frittage du catalyseur (environ 900 K), le réacteur sera détruit.
📚 Référentiel Théorique
Calorimétrie AdiabatiqueCapacité Thermique des GazLe principe est celui de la conservation de l'enthalpie (\(\Delta H_{\text{système}} = 0\)) dans un système isolé. Nous allons décomposer la transformation en deux étapes fictives pour le calcul : 1. La réaction se produit à température constante (700K), libérant une quantité de chaleur \(Q_{\text{réac}}\). 2. Cette chaleur \(Q_{\text{réac}}\) est intégralement réinjectée dans le mélange gazeux final pour élever sa température de \(T_{\text{initiale}}\) à \(T_{\text{finale}}\).
L'énergie sensible est la quantité de chaleur nécessaire pour changer la température d'un corps sans changer sa phase. Elle est donnée par la relation \(Q = n \cdot C_p \cdot \Delta T\), où \(C_p\) est la capacité calorifique molaire.
📋 Données d'Entrée
| Paramètre | Valeur |
|---|---|
| \(\Delta_{\text{r}} H^\circ\) | \(-197,8 \text{ kJ/mol}\) |
| Avancement \(\alpha\) | 0,992 |
| \(C_{p,\text{mix}}\) (Moyen) | \(\approx 32 \text{ J/mol}\cdot\text{K}\) |
| \(T_{\text{initiale}}\) | 700 K |
Calculs Détaillés
1. Calcul de la Chaleur Libérée (Source) :
Nous calculons d'abord l'énergie totale fournie par la réaction chimique. Nous multiplions le nombre de moles d'avancement (le nombre de fois que la réaction s'est produite) par l'enthalpie de réaction unitaire.
Interprétation : La réaction libère environ 981 kJ pour 100 moles de gaz injecté.
2. Calcul de l'Inertie Thermique (Puits) :
Nous déterminons ensuite la capacité calorifique totale du mélange gazeux qui va absorber cette chaleur. C'est le produit du nombre total de moles final par la capacité calorifique molaire moyenne.
Interprétation : Il faut 3040 Joules pour augmenter la température de l'ensemble du gaz de 1 degré Kelvin.
3. Calcul de l'Élévation de Température \(\Delta T\) :
Nous appliquons la relation calorimétrique fondamentale \(Q = C \Delta T\), que nous inversons pour trouver \(\Delta T\). Nous divisons l'énergie disponible par la capacité thermique.
Interprétation : La température du gaz s'élève de plus de 322 degrés.
4. Température Finale :
Nous ajoutons cette élévation à la température initiale de 700 K.
Interprétation Critique : La température finale adiabatique atteindrait théoriquement 1022 K (750°C). Or, nous savons que le catalyseur se dégrade irréversiblement au-delà de 900 K. De plus, à 1022 K, l'équilibre serait totalement inversé (Kp deviendrait minuscule), stoppant la réaction bien avant d'atteindre ce point.
Ce calcul de bilan thermique est la pierre angulaire de la conception du réacteur. Il démontre formellement qu'il est physiquement impossible de réaliser la conversion complète en une seule passe adiabatique. L'énergie libérée est telle qu'elle porterait le système à des températures destructrices (1022 K) bien au-delà de la résistance du catalyseur (900 K) et de l'équilibre thermodynamique favorable.
Le saut de température calculé (\(\Delta T \approx 322\) K) est tout à fait cohérent avec l'enthalpie de réaction très élevée trouvée à la question 1. Une réaction très exothermique doit nécessairement entraîner une élévation de température massive si le système n'est pas refroidi. Cela valide a posteriori l'ensemble de la démarche thermodynamique.
La capacité calorifique \(C_p\) du gaz n'est pas constante mais augmente avec la température. Le calcul effectué avec un \(C_p\) moyen est une approximation. Un calcul plus fin (intégrale du \(C_p(T)\)) donnerait une température finale légèrement plus basse, mais la conclusion resterait identique : le seuil de sécurité est largement dépassé. C'est pourquoi l'ingénierie impose la fragmentation du réacteur en 4 lits successifs avec refroidissement intermédiaire, comme vu sur le schéma initial.
📄 Livrable Final (Note de Synthèse)
| Ind. | Date | Objet de la modification | Rédacteur |
|---|---|---|---|
| A | 01/10/23 | Création du document préliminaire | Ing. Jr. |
| B | 24/10/23 | Calculs thermodynamiques & Validation Refroidissement | Ing. Procédés |
- Modèle Gaz Parfaits applicable (Haute T, Basse P).
- Enthalpies de formation issues de la base JANAF.
- Catalyseur : Pentoxyde de Vanadium sur support silice.
| Débit molaire Entrée | 100 kmol/h (Base calcul) |
| Ratio SO2/O2 | 0,91 (Excès d'air) |
| Température Entrée | 700 K |
Validation des performances du premier lit catalytique.
Dépt. Procédés
Dir. Technique
Laisser un commentaire