TPAspen Complet PDF

Modlisation des
Grands Systmes Chimiques

Professeur Georges HEYEN
Universit de Lige, Anne acadmique 2012-2013
Notes de Travaux Pratiques
Laboratoire dAnalyse et Synthse des Systmes Chimiques
Notes de TPs rdiges par Carine GERKENS, Alexandre LERUTH, Grgoire LEONARD
Septembre 2012
Travaux Pratiques SYST0004
Introduction
Logiciels de simulation utiliss
Au cours de ces Travaux Pratiques, nous utiliserons plusieurs logiciels permettant de
construire et dutiliser des modles de procds chimiques : Aspen Plus, Aspen Hysys et
ValiEngineering (Belsim). Les logiciels ASPEN Advanced System for Process ENgineering
sont conus par AspenTech, entreprise base au Massachusetts, USA. Fonde en 1981,
AspenTech est ne dun projet de recherche commun entre linstitut de technologie du
Massachusetts (MIT) et le dpartement de lEnergie des Etats-Unis dAmrique (DOE).
Diffrents logiciels ont t conus par AspenTech pour lindustrie des procds, laidant
amliorer ses performances en matire dobjectifs oprationnels. Des modles de simulation
prcis permettent en effet aux entreprises daccrotre lefficacit et la rentabilit de leurs
procds. La gamme propose par AspenTech comprend des logiciels de modlisation de tous
types de procds, de conception et dvaluation dchangeurs, danalyse nergtique et
systmique, dvaluation conomique, de dveloppement de procds et de supports
oprationnels ainsi quune documentation dtaille pour chacun dentre eux.
Aspen Plus et Aspen Hysys sont les principaux logiciels dAspenTech orients gnie des
procds. Aspen Plus est un outil de modlisation ayant trait la conception, loptimisation et
le suivi des performances de la plupart des procds industriels. Aspen Hysys est un outil de
modlisation ayant notamment trait la conception, loptimisation et au suivi des
performances des procds ptrochimiques (production, raffinage), gaziers et de sparation
dair.
Nous utiliserons galement les logiciels Vali de la socit Belsim. Belsim est une entreprise
belge, spin-off de lUniversit de Lige issue du Laboratoire dAnalyse et de Synthse des
Systmes Chimiques (LASSC). Elle est spcialise dans la validation de donnes. La
validation est une mthode de traitement et danalyse des mesures ralises dans un procd.
Elle vrifie la cohrence des mesures par rapport au modle. Elle permet galement de
dtecter et de corriger des erreurs, et dvaluer les paramtres de performance du procd.
Vali dsigne la batterie de logiciels de validation conus par Belsim.
Etude dun procd avec le logiciel Aspen Plus

Nous aborderons tout dabord par quelques exemples les fonctionnalits de base disponibles
dans le logiciel Aspen Plus. En exploitant ensuite l'exemple d'une boucle de production
d'ammoniac, nous illustrerons les concepts principaux discuts dans le cours thorique en
faisant appel au logiciel d'analyse et de simulation de procds chimiques Aspen Plus. Nous
donnerons l'ordinateur une description du procd, le simulerons et tudierons sa rponse
des variations de conditions opratoires.
Aprs avoir considr le racteur de faon isole, nous l'intgrerons dans une squence
d'appareils (pompes, changeurs, ), dont les conditions de fonctionnement seront varies au
moyen dune tude paramtrique. Les recyclages de matire et d'nergie seront traits par la
mthode des coupures et le calcul squentiel modulaire. La gestion des spcifications et les
problmes numriques de convergence seront abords.
Septembre 2012
Ensuite, la simulation de colonnes de distillation illustrera les techniques de rsolution globale

de grands systmes d'quations. Les logiciels Aspen Plus et Hysys seront tous deux utiliss
pour cet exemple.
Pour terminer, nous examinerons le problme de la validation de donnes exprimentales qui
nous permettra de vrifier la cohrence des mesures et d'identifier les paramtres d'une
installation existante. Nous nous servirons pour cela du logiciel de validation de mesures
BelsimVali.
Droulement des travaux pratiques

Les travaux pratiques seffectuent de manire indpendante en dehors des sances de TP. Les
logiciels Aspen Plus et Hysys peuvent tre install sur vos ordinateurs personnels (CD
dinstallation demander lassistant). Lobjectif est dassurer une bonne comprhension de
la matire et des logiciels employs, tout en vous permettant de travailler votre rythme.
Toutes les trois semaines, vous recevrez les noncs de trois exercices concernant les logiciels
Aspen Plus et Hysys. La dernire partie concernant la validation de donne sera organise de
faon plus classique et stalera sur les deux dernires sances de TPs. Les travaux pratiques
seront structurs comme suit :
-
1re Partie : Introduction au logiciel Aspen Plus et simulation de squences simples

de procds
2me Partie : Boucles dans un procd et spcifications de design
3me Partie : Optimisation de procds industriels, simulation de colonnes distiller

et introduction au logiciel Hysys.
4me Partie : Introduction la validation de donnes et au logiciel Vali de Belsim.
Par groupes de deux, vous devrez raliser une application personnelle conscutive chaque
partie. Le rapport relatif chaque partie devra tre rendu en version papier au plus tard trois
semaines aprs le TP. Lassistant sera votre disposition pour rpondre aux questions si
ncessaire. Lobjectif est de vous permettre de vous familiariser aux logiciels de simulation
couramment employ dans le monde industriel tout en vous incitant en explorer les
nombreuses possibilits par vous-mme.
Ces rapports devront tre rdigs en respectant les rgles de lart : introduction et conclusion,
illustration graphique pertinente des phnomnes dcrits, mise en page, orthographe, Il
sagira de rpondre aux questions poses en tant synthtique, sans ajouter de rsultats
inutiles et en gardant lesprit que les rapports seront corrigs par une personne qui connat
les noncs.
Septembre 2012
Exercice 1.1 : Echangeur de chaleur avec flash1

Notre premier procd consistera en le refroidissement puis la sparation dun flux
dhydrocarbures. Aprs refroidissement 180K dans un changeur de chaleur, le flux subit
une sparation de ses phases liquide et vapeur dans un sparateur flash comme reprsent cidessous.
3
HTX
FLASH
Figure 1 : Schma (flowsheet) de lexercice 1
Simulation
o
Lancez linterface du logiciel Aspen Plus V7.3.2 :

o Tous les programmes/AspenTech/Process Modeling V7.3.2 /Aspen Plus /Aspen
Plus V7.3.2. => Vous arrivez sur la page daccueil du logiciel.
o Sur la page daccueil, ouvrez une nouvelle simulation (New). Choisissez un
modle (template) de type User, General with metric units (Units SI avec
degrs Celsius et bar).
Vous arrivez sur linterface du logiciel, compose de deux parties principales: Properties
et Simulation.
o Choisissez linterface Simulation en bas gauche de votre cran. Nous
reviendrons aux proprits par la suite.
Le premier travail effectuer est celui de la conception du Flowsheet ou Process Flow

Diagram (PFD) , le schma du procd. Pour ce faire, il faut dessiner les diffrentes
units (ou blocs, correspondant aux appareils et aux oprations raliser) et les relier par
des flux (correspondants aux changes de matire et d'nergie entre les units et vers
l'environnement).
o Utilisez la palette de modles en bas de votre cran pour insrer les blocs Heater
(Heat exchangers) et Flash2 (Separators).
o Positionnez les deux blocs sur votre schma. Une fois les blocs placs, la touche
Escape vous permet de retrouver un curseur normal.
Dfinissez les flux de matire ncessaires au procd. Pour cela, cliquez dans la palette
des modles sur loption Material streams. A ce moment apparaissent en rouge sur le
schma les flux de matire qui doivent tre ncessairement relis chacun des blocs, i.e.
chacune des oprations unitaires. En bleu, les flux de matire possibles mais non
1
Les exercices 1.1 et 1.2 proposs dans ces notes sont librement inspirs du 3me TP de Technique des
Procds donn lUniversit Technique de Munich par le Pr. J. Stichlmair, anne 2007-2008.
Septembre 2012
ncessaires. En reliant des flux dun bloc un autre, ils seront automatiquement dfinis
comme flux dentre et/ou de sortie.
o Reliez les blocs comme indiqu la figure 1. Attention bien slectionner les
flches rouges du schma pour que le lien entre le flux et le bloc soit bien tabli.
o
Une fois le schma termin, revenez linterface Properties du procd en bas gauche
de votre cran. Le panneau de navigation sur la gauche de votre cran vous aide connatre les
lments ncessaires complter avant de dmarrer la simulation. Les ronds rouges et blancs
indiquent que des donnes doivent encore tre fournies. En navigant au moyen de
lexplorateur sur votre gauche, encodez les composants du procd.

o Tapez le nom que vous souhaiter donner aux composs dans la case Component
ID . La plupart du temps, entrer le nom anglais du composant suffit ce
quAspen Plus le reconnaisse et affiche son nom et sa formule. Si ce nest pas le
cas, utilisez loption de recherche Find pour trouver toute substance sur base de
son nom commun, de son numro CAS ou de sa formule brute.
o Les 4 composs encoder sont les suivants: mthane, thane, propane, n-butane.
o
Dans longlet Methods du navigateur sur la gauche de lcran, vous pouvez dfinir la
mthode thermodynamique de base utiliser pour le calcul des proprits. Choisissez le
modle de contribution de groupes Unifac comme modle thermodynamique. Nous ne
rentrerons pas dans le dtail des autres possibilits de cette page pour le moment.
Maintenant que la partie Properties est entirement complte (il ny a en principe plus
de ronds rouge et blanc dans le navigateur), nous revenons la partie Simulation o nous
voyons quil reste spcifier les flux (Streams) et les units (Blocks).
o Spcifiez les caractristiques du flux dentre du procd: 20C, 2bar.
Composition : 40 kmol/h de mthane, 30 kmol/h dthane, 20 kmol/h de propane
et 10 kmol/h de n-butane
o Spcifiez les caractristiques de lchangeur de chaleur : il a pour mission
damener le flux une temprature de 180K. Sa perte de charge est de 0.5bar.
Elle peut tre spcifie par lencodage dun chiffre ngatif (-0.5bar) ou en
encodant la pression absolue de sortie : 1.5bar. Notons que ce type de bloc
(Heater) ne considre quun seul flux de matire dans lchangeur de chaleur.
o Spcifiez les caractristiques du ballon flash : perte de charge de 0.5 bar,
fonctionnement adiabatique (change de chaleur nul avec lenvironnement, pas
dapport ni perte de chaleur). Il sagit donc dun flash dtermin par sa pression et
son change de chaleur nul avec lextrieur (Heat Duty=0).
Ouvrez maintenant le panneau de contrle (touche F7 ou bouton Control Panel dans la

barre doutils Run en haut de votre cran, onglet Home). Lancez la simulation laide
du bouton Play dans cette mme barre doutils (ou touche F5).
Cliquez sur le bouton Check Results dans cette mme barre doutils et parcourez les
rsultats en utilisant le navigateur. Les rsultats sont visibles dans les onglets Results
pour les flux et dans les onglets Results et Stream Results pour les blocs.
o Pour afficher les rsultats de la simulation sur le schma, cliquez dans longlet Home sur
Stream Summary. Cliquez sur Stream Table pour intgrer le rsum de vos rsultats au
flowsheet. Plus simplement, vous pouvez ajouter sur le schma les valeurs de
temprature, pression et fraction vapeur dans longlet Flowsheet/Modify.
Septembre 2012
Application
1.
Au vu de ce premier exercice, quelles sont les tapes principales de la

construction dun modle de procd dans Aspen Plus ?
2.
Quelle est la quantit de chaleur (en GJ/hr) qui est retire au flux dentre dans cet
changeur ? A quoi sert de refroidir le flux dentre dans le procd ?
3. Un ballon flash possde deux degrs de libert. Si au lieu dimposer un flash

adiabatique on spcifie une temprature de flash de 190K (en conservant la perte
de charge), comment volue la puret en mthane du flux gazeux? Discutez
brivement.
Exercice 1.2 : Analyse des proprits des composs

A ct dAspen Plus, il existe dautres programmes tels quAspen Properties, dont le rle est
lanalyse des proprits des composs. Cette analyse des proprits de corps purs et de
mlanges est possible dans Aspen Plus indpendamment des calculs dun procd. Dmarrer
une nouvelle simulation sur un modle General with metric units .
Dans la partie Properties de linterface, encodez leau comme composant, et choisissez

Ideal comme mthode pour le calcul des proprits.
Cliquez sur le bouton Pure dans la barre doutils Analysis, onglet Home.
o Construisez la courbe de tension de vapeur de leau pure sous la forme dun
diagramme P-T. La proprit choisir est PL, la pression de vapeur dun compos
pur.
o Slectionnez leau comme composant et construisez la courbe de tension de vapeur
de leau pour des tempratures variant entre 0 et 450C. 51 points suffisent pour
une bonne approximation.
o Cliquez sur le bouton Run analysis. Les rsultats apparaissent sous la forme dun
graphique. Vous pouvez les retrouver dans un tableau : dans le navigateur sur la
gauche de votre cran, cliquez sur longlet Results de lanalyse.
Application
4. Aux alentours des usines chimiques, il existe souvent diffrents rseaux de vapeur,
avec par exemple de la vapeur 1barg, 5barg et 20 barg2. Quelle temprature
doit avoir la vapeur au minimum dans chacun de ces rseaux pour ne pas
condenser ? Adaptez votre diagramme P-T au domaine de temprature adquat et
dduisez-en les tempratures de vapeur.
barg = bar gauge = pression la jauge ou au manomtre. Cest une pression relative qui
sajoute la valeur de la pression atmosphrique, i.e. : 5barg = 6bara en pression absolue
2
Septembre 2012
5. Construisez le diagramme T-xy du mlange Ethanol/Eau pression

atmosphrique (rajoutez les composs ncessaires dans le modle). Pour ce faire,
cliquez sur le bouton Binary dans la barre doutils Analysis, onglet Home et
compltez ce qui est ncessaire. Le rsultat obtenu est-il correct ? Discutez.
Exercice 2 : Raction de synthse dammoniac

Nous allons maintenant examiner un procd de synthse dammoniac. Ce procd sera
construit par tapes successives au fur et mesure des TPs. En premier lieu, nous tudions le
racteur chimique de synthse en reprenant depuis le dbut certaines notions du logiciel
Aspen Plus que nous avions laisses de ct dans un premier temps.
REOUTG
REACTEUR
REIN
REOUTL
Figure 2 : Racteur de synthse dammoniac
Simulation
Lancez le logiciel Aspen Plus, et dmarrez un modle de type General with metric units.
Une vrification de la licence seffectue (vous devez tre reli au rseau de lUniversit
de Lige). Enregistrez une premire fois le document. Celui-ci senregistre avec
lextension .apv. Il sera parfois utile dutiliser le fichier de back-up .bkp qui permet de
faire passer le mme modle dune version une autre dAspen Plus.
Nous passons dans la partie Simulation de linterface. En bas de la page, dans la palette
des modles, choisissez le modle de racteur REquil et placez cette unit au centre de
lcran. REquil modlise un racteur dans lequel lquilibre est fix au moyen une
approche stchiomtrique. Il calcule simultanment les quilibres chimiques et les
quilibres entre phases (pas plus de deux). Dans notre cas, cela nous permettra de calculer
la raction qui a lieu en phase gazeuse en supposant que lquilibre est atteint.
o Vous pouvez renommer le bloc REACTEUR en effectuant un clic droit sur
celui-ci ou en utilisant la touche de raccourci Ctrl+ M.
Insrez les flux dentre et de sortie cette unit en slectionnant Material Stream dans
lextrmit infrieure gauche de lcran. Plusieurs types de flux existent :
- Flux de matire
- Flux thermiques (puissance d'un rebouilleur, par exemple)
- Flux mcaniques (puissance d'un compresseur, par exemple)
Septembre 2012
Connectez trois flux de matire au racteur, un lentre et deux en sortie

selon le schma ci-dessus. Vous pouvez les renommer respectivement
REIN , REOUTL et REOUTG .
De retour dans la partie Properties, dans la rubrique Components du navigateur, entrez

les composants intervenant dans le procd de synthse de lammoniac : argon, mthane,
azote, hydrogne et ammoniac.
Il faut maintenant dfinir le modle qui va tre utilis pour le calcul des proprits des
composants purs ainsi que des mlanges. Un modle idal ne conviendrait pas dans ce
procd o la pression montera jusque 280 bars. Nous utiliserons donc l'quation d'tat de
Redlich-Kwong-Soave pour estimer les coefficients d'quilibre liquide-vapeur et l'enthalpie des
fluides.
o Dans la rubrique Methods du navigateur, choisissez la mthode RK-Soave (elle

ne fait pas partie des mthodes employes pour dcrire les procds communs,
vous devez donc spcifier All dans le filtre de mthode).
o Les paramtres dinteraction binaire sont disponibles mais Aspen Plus vous
demande de les approuver avant de pouvoir les utiliser.
o Pour amliorer les rsultats, on peut introduire les valeurs ajustes des
coefficients dinteraction binaire, lues dans une table de donnes, ou
dtermines au moyen dun logiciel dajustement paramtrique sur base
dinformations exprimentales. On encodera donc les valeurs suivantes dans
Properties/ Parameters/ Binary Interactions/ RKSKBV-1 :
Substance i Substance j Valeur de KAij
N2
NH3
0.224
H2
NH3
-0.89
AR
NH3
0.346
CH4
NH3
0.197
Les autres paramtres seront fixs 0.
Rem. : Normalement, des coefficients d'interaction sont recommands pour les systmes nonidaux pour chacun des binaires faisant intervenir les substances majeures. Dans notre cas,
comme l'ammoniac sera le constituant principal de la phase liquide, nous introduisons les 4
paramtres ci-dessus relatifs chacune des autres substances en paire avec l'ammoniac.
Revenez dans linterface Simulation pour la dfinition des flux et des units. Dans la
rubrique Streams, le flux dentre dans le racteur (REIN) est une temprature de
380C, une pression de 280 bar, et son dbit vaut 0,5 kmol/sec. La composition du feed
est la suivante (en fractions molaires, le total doit donc valoir 1) :
Composs
AR
CH4
N2
H2
NH3
Septembre 2012
Fractions molaires :
0.01
0.01
0.245
0.735
0
Dans la rubrique Blocks, le racteur fonctionne de faon adiabatique et la pression de

280 bar.
Pour signifier au programme les ractions qui devront tre prises en compte dans le
racteur, il faut dfinir celles-ci dans longlet Reactions du bloc REACTEUR3. En
loccurrence, la seule raction encoder est celle de synthse de lammoniac :
N2 + 3H2 2NH3
Rem. : Le racteur REquil permet de calculer des ractions qui natteignent pas ltat
dquilibre. Pour cela, on peut imposer un degr davancement dtermin (molar extent) ou
calculer la raction une temprature dquilibre diffrente de la temprature du bloc. En
spcifiant une approche de temprature de +5C, la constante dquilibre Kp serait donc
calcule une temprature plus leve que la temprature dquilibre du racteur. Mais nous
supposerons ici que lquilibre est bien atteint la temprature du racteur.
Afin de choisir le systme dunits dans lequel vous voulez voir apparatre les variables,
slectionnez Setup/ Specifications/ Global. Donnez un nom votre simulation et
choisissez Si-CBAR comme set dunits. Les flux entrant et sortant du racteur
contiennent plusieurs substances, nous allons donc demander au programme de nous
indiquer par dfaut les fractions molaires des diffrents flux. Dans la rubrique
Setup/Report Option/Stream, cochez la case Mole dans la catgorie Fraction Basis.
Vous pouvez au choix galement dcider dobtenir les fractions et les flux massiques.
Maintenant que toutes les donnes ont t entres, lancez le calcul et observez les
rsultats. Vous vrifierez aisment les rsultats de calcul pour les diffrents flux dans
longlet Streams du navigateur. Vous constaterez que la fraction molaire en NH3 (12.4%)
est relativement faible dans le flux de sortie REOUTG. Ceci fera lobjet des exercices
suivants, o nous introduirons le principe du recyclage de flux.
Etude paramtrique
Nous allons examiner linfluence de la temprature du flux dentre sur la fraction molaire en
NH3 la sortie du racteur. Pour ce faire, nous allons mener une tude de sensibilit.
Dans la rubrique Model Analysis Tool accessible partir du navigateur, crez une
nouvelle tude de sensibilit que vous nommez TEMP.
Dans longlet Vary, dfinissez une nouvelle variable TEMP. Il sagit dune variable de
flux : la temprature du flux dentre dans le racteur, faire varier entre 250 et 500C
(25 points suffisent).
Dans longlet Define, dfinissez une variable XNH3 et spcifiez-la comme tant la
fraction molaire en NH3 du flux de sortie en phase gazeuse.
3
Il est galement possible de dfinir des ractions dans la rubrique Reactions du navigateur. Cela a
lavantage de dfinir un set de ractions soit dlectrolytes (Chemistry) soit gnrales (Reactions) qui pourra
tre utilis dans diffrentes units du procd sans devoir tre redfini chaque unit. Dans le cas du modle
REquil, il est cependant ncessaire de spcifier les ractions lors de la spcification du racteur.
Septembre 2012
Dans longlet Tabulate, demandez au programme de lister dans la colonne 1 les valeurs
de XNH3 obtenues lorsque la temprature de REIN varie. Sauvegarder vos rsultats.
Application
6. Selon ltude de sensibilit, la raction est-elle endo- ou exothermique ? En restant
dans les bornes poses pour cette tude de sensibilit, quelle temprature de flux
dentre choisiriez-vous pour maximiser la synthse dammoniac ?
7. En pratique, le procd industriel de synthse dammoniac (bas sur le procd
Haber-Bosch) fonctionne des tempratures de lordre de 450-600C. Est-ce en
accord avec les rsultats de votre simulation ? Quelle pourrait tre la raison
cela ?
8. De la mme manire que cela a t fait pour la temprature du flux dentre,
ralisez une tude paramtrique pour examiner linfluence de la pression du
racteur sur la synthse dammoniac. Utilisez pour cela un racteur adiabatique.
Quelle pression prconisez-vous ?
9. Discutez de lintrt et du risque du point de vue industriel dune simulation
utilisant le logiciel Aspen Plus en regard de vos rponses aux deux points
prcdents.
Septembre 2012
10
Exercice 3 : Calcul dune squence dunits

Introduction
La synthse de l'ammoniac implique une raction dont nous avons tudi les conditions
opratoires. Pour raliser un procd complet, le racteur doit tre accompagn en amont
d'une section de prparation et en aval dune section de purification.
Le mlange ractionnel entrant dans le procd dans des conditions de temprature et de
pression (310 K et 30 bars) fort diffrentes des conditions ractionnelles (850 K et 280 bars),
il conviendra de prvoir en amont un dispositif de compression-change de chaleur que nous
dterminerons par la simulation. Dautre part, par l'analyse des rsultats de la simulation du
racteur, nous constatons que la conversion est relativement faible, ce qui implique que
l'ammoniac produit se retrouve en faible proportion la sortie du racteur. A la pression de
travail qui est la ntre, il serait trs simple et efficace de condenser par refroidissement un
liquide riche en ammoniac, et de poursuivre par une simple sparation liquide-vapeur. Ces
modifications sont schmatiquement reprsentes la figure 3 :
Figure 3 : schma du procd amlior de synthse dammoniac
Simulation
Nous allons commencer par tudier la condensation dun flux liquide riche en ammoniac la
sortie du racteur. La premire tape de ce travail consistera dterminer la temprature
adquate pour la liqufaction. Une telle recherche peut se raliser trs facilement avec le
logiciel de simulation. Reprenez la simulation ralise lors de lexercice prcdent et
dsactivez les tudes paramtriques (clic droit sur ltude de sensibilit => hide).
Dans linterface Simulation, affichez le schma du procd (Onglet View/ Flowsheet).

Rajoutez un bloc Flash2 de type Separators. Placez-le droite de votre REACTEUR et
renommez ce ballon FLASH. Insrez deux flux de matire sortants ce ballon.
Renommez ces flux FLASHLIQ (flux en pied) et FLASHVAP (flux en tte). Il reste
alors raccorder le flux de sortie du racteur REOUTG lentre du ballon. Vous pouvez
faire cette connexion en ralisant un clic droit sur le nom du flux => Reconnect
Destination.
Septembre 2012
11
Vous pouvez accder la spcification du ballon FLASH en double-cliquant dessus. La

pression est de 278,76 bars et le taux de vaporisation est de 1, car nous dsirons pour le
moment calculer le point de rose.
Lancez le calcul et observez les rsultats. Le point de rose prdit par le programme est
de 46C et prs de 7900kW doivent tre enlevs au ballon flash pour permettre la
condensation4.
Application
10. La courbe de refroidissement peut tre examine en effectuant une tude de
sensibilit sur le taux de vaporisation ou sur la temprature du ballon Flash. Faites
varier un de ces deux paramtres au choix et observez les rsultats obtenus pour le
deuxime. Reprsentez la courbe de refroidissement (T en fonction de Vfrac, ou
vice-versa).
Rem. : Attention bien dfinir la bonne variable pour le calcul de sensibilit : Temp si cest
la temprature spcifie dans le flash (calcul de la fraction vapeur), TCalc si cest la
temprature calcule (spcification de la fraction vapeur). Idem pour Vfrac et Vcalc.
11. Reprsentez une deuxime variable dans votre tude de sensibilit en fonction de
la temprature de condensation: la fraction molaire de NH3 dans le flux condens.
Quelle temprature de condensation choisiriez-vous ? Quel est le compromis
trouver ?
Nous allons maintenant rajouter la prparation des ractifs avant leur entre dans le racteur.
La compression amont implique un rapport de pression de l'ordre de 9, ce qui sort largement
du cadre des performances habituelles des compresseurs, ce qui nous amne dfinir deux
tages avec refroidissement intermdiaire en vue d'accrotre les performances de l'ensemble.
Le schma est reprsent la figure 4.
Figure 4 : Schma de la synthse dammoniac avec train de compression
Les rsultats peuvent varier lgrement selon les conditions initiales du calcul. Pour remettre zro
toutes les conditions initiales et rinitialiser la simulation, cliquez sur le bouton reinitialize sur la droite du
bouton Play dans le haut de votre cran.
Septembre 2012
12
De retour sur la fentre du schma, rajoutez un changeur de chaleur HTX3 (modle

HEATER) entre le racteur et le sparateur flash (voir figure 5). Le fluide sort de cet
changeur la temprature de 10C et sans perte de pression (0bar). Lchangeur
fonctionne en mode biphasique liquide-vapeur.
FLASHVAP
FLASH
HTX3
REOUTG
FLASHIN
REACTEUR
REIN
FLASHLIQ
REOUTL
Figure 5 : ajout dun changeur de chaleur avant le sparateur flash
En vous conformant au schma ci-dessous (figure 6), ajoutez deux compresseurs (modle
COMPR) COMP1 et COMP2. Ce type de compresseur est le modle de base pour les
compresseurs isentropiques et polytropiques (ce modle peut galement servir pour les
turbines). Aprs chacun des compresseurs, ajoutez un changeur de chaleur (HEATER),
et nommez-le HTX1 resp. HTX2. Ajoutez un flux dentre de matire que vous
nommerez F1, et reliez les nouvelles units par des flux nomms de F2 F4.
FLASHVAP
FLASH
HTX3
F1
REOUTG
COMP1
COMP2
FLASHIN
REACTEUR
FLASHLIQ
F3
REIN
F2
HTX1
F4
HTX2
REOUTL
Figure 6 : Schma du rsultat final de lexercice 3
Les caractristiques du compresseur

COMP1 sont les suivantes : compresseur
isentropique, augmentation de pression de 60 bars, efficacit isentropique de 0,75 et
efficacit mcanique de 0,95. Elles valent galement pour COMP2 si ce nest que la
sortie du 2me compresseur se fait 280 bars. Aucun des deux changeurs de chaleur ne
provoque de perte de charge. La temprature la sortie de HTX1 vaut 40C, celle la
sortie de HTX2 vaut 380C.
Septembre 2012
13
Spcifiez les caractristiques du flux dentre F1 : mme composition que REIN (1mol% dargon ; 1mol-% de mthane ; 24,5mol-% dazote et 73,5mol-% dhydrogne),
pression de 30 bars, temprature de 20C et dbit de 0,5kmol/sec. Etant donn que REIN
nest plus le flux dentre du procd, vous pouvez supprimer les valeurs encodes car
elles seront de toute faon recalcules par Aspen Plus. Si vous ne supprimez pas les
valeurs donnes REIN, le problme sera sur-spcifi, mais cela nempchera pas le
calculateur de trouver une solution.
Vous pouvez galement ajouter des flux de chaleur sortants aux changeurs de chaleur en
cliquant sur Heat Streams en bas de votre cran (menu droulant ct de material
streams). Si vous ne le faites pas, il sera quand mme possible de connatre la quantit de
chaleur change en consultant simplement les rsultats du bloc. Cette option ne change
rien au calcul, mais permet la visualisation des flux.
Adaptez les spcifications du racteur et du flash comme suit : la perte de charge dans le
racteur adiabatique vaut maintenant 1,24 bar (indiquez une pression ngative pour une
perte de charge) et le flash a lieu dans un ballon qui fonctionne de faon adiabatique et
sans pertes de charge. Vous pouvez alors relancer le calcul (dsactivez toutes les tudes
de sensibilit encore apparentes). Vrifiez les rsultats. La fraction molaire de NH3 dans
FLASHLIQ est de 0,979.
Application
12. Examinez le fichier *.his gnr et regardez linformation que vous pouvez en
tirer. Faites un bref rsum des principaux points disponibles.
13. Le modle que vous avez construit est le cas de base partir duquel on peut
explorer les conditions opratoires de la synthse d'ammoniac. Pour chacune des
configurations suivantes, tudiez l'volution par rapport au cas de base du dbit
massique de liquide condens dans le flash, de la puret molaire de l'ammoniac
dans le flux condens, et des besoins en nergie (chaleur, froid, lectricit) des
diffrents quipements :
- pression au racteur 30 et 500 bars (racteur adiabatique)
- temprature de raction isotherme du racteur 50 et 750C (perte de charge de
1.24 bar dans le racteur)
- temprature de condensation du ballon flash 200 et 100C (perte de charge
ngligeable).
Si ncessaire, adaptez les autres quipements pour former un procd cohrent
avec la configuration donne.
14. a. Dans quel cas minimise-t-on la consommation nergtique du procd ?
Attention, toutes les nergies nont pas la mme valeur, 1 kW de travail mcanique
ou dlectricit a plus de valeur que 1 kW de chaleur : 1kWl = 1kWmc = 3kWth !
b. Dans quel cas maximise-t-on la production en NH3 liquide?
c. Dans quel cas minimise-t-on le ratio kW/kmol dNH3 liquide?
Septembre 2012
14
Exercice 4 : Recyclage de matire

Lorsqu'on examine les rsultats obtenus lors des simulations prcdentes, on remarque que le
rendement de la raction est assez faible, ce qui nous impose de recycler les ractifs nayant
pas ragi pour ne pas gaspiller le gaz de synthse. Cette boucle ncessite lajout dun
troisime tage de compression, qui doit compenser les pertes de charge dans la boucle et
faire circuler le gaz dans le racteur.
L'effluent gazeux du sparateur est partiellement recycl l'admission du dernier tage de
compression. Une partie du gaz est purge pour viter l'accumulation des substances qui ne
sont pas limines par raction et condensation, c'est--dire l'argon et le mthane. En effet, la
solubilit de ces substances dans l'ammoniac est insuffisante pour qu'elles s'liminent
facilement par mise en solution dans l'ammoniac liquide. L'absence de purge conduirait
l'accumulation des inertes jusqu' un taux qui diminuerait considrablement la productivit du
racteur.
Ce type de problme est beaucoup plus complexe
modliser car il est ncessaire d'introduire une
procdure itrative. Elle aura pour but de trouver
les valeurs des variables coupes satisfaisant la
solution. Pour ce faire, nous avons dans le logiciel
un certain nombre d'outils mathmatiques notre
disposition : diffrentes mthodes de rsolution y
correspondent.
Dans notre cas, nous pouvons faire une simple
coupure sur les dbits partiels sur le flux
FLASHIN (le flux d'entre du sparateur liquide
vapeur). On pourrait galement appliquer
lalgorithme de Mottard pour confirmer ce choix.
Figure 8 : schma du recyclage de matire
Reprenez le modle obtenu la fin de lexercice 3 (pour rappel, le racteur est spcifi
adiabatique et occasionnant une perte de charge de 1.24 bar, et le flash est adiabatique et sans
perte de charge). A laide de la souris, slectionnez un large cadre autour de la chane de
compression et placez-la en haut gauche de la partie racteur et flash du procd en
maintenant le clic de la souris enfonc. Cela nous permettra dintroduire entre lchangeur de
chaleur HTX2 et le flux REIN divers lments ainsi que reprsent la figure 9 :
-
un mlangeur de flux (mixer) MIX (que vous ferez tourner de 90, utilisez pour
cela le clic droit -> rotate right). La perte de charge y est nulle.
un compresseur COMP3 isentropique, avec une efficacit isentropique de 0,75 et

une efficacit mcanique de 0,95. La pression de sortie est de 280 bars. Ramenez
la pression de sortie de COMP2 270 bars. Ce 3me tage de compression aura
pour but de compenser les pertes de charge du flux provenant du flash et qui sera
recycl par lintermdiaire du mlangeur de flux.
Septembre 2012
15
un changeur de chaleur (Heater) que nous nommerons HTX3 aprs avoir

renomm le condenseur avant le flash en HTX4. La temprature de sortie de ce
nouvel changeur de chaleur est de 380C, et sa perte de charge est de 1,07bar.
Ces units sont relies entre elles par les flux F5 F7. Reconnectez le flux REIN la sortie de
HTX3. Si ce nest dj fait, dsactiver ltude de sensibilit ralise lexercice 3 avant de
lancer le calcul. Pour ce faire, dans Model Analysis Tool, Sensitivity, effectuez un clic droit
sur ltude S-1, et cliquez sur Hide. Vous pouvez lancer le calcul pour vous assurer quaucune
erreur na t commise dans lencodage. La quantit dammoniac liquide en sortie du flash
na presque pas t influence par les modifications puisque le recyclage na pas encore eu
lieu.
SPLIT
PURGE
FLASHVAP
RECYC
F5
HTX2
MIX
FLASH
HTX4
F6
FLASHIN
COMP3
FLASH2G
REIN
F7
HTX3
FLASHLIQ
REACTEUR
FLASH2
FLASH2L
REOUTL
Figure 9 : Boucle de recyclage dans Aspen Plus

Dans le but damliorer la puret du produit final, ajoutez un second sparateur FLASH2 en
sortie (liquide) du premier sparateur FLASH, et nommez les nouveaux flux sortants
FLASH2L et FLASH2G. Ce sparateur est adiabatique et fonctionne la pression de 10 bars.
Enfin, actualisez la perte de charge dans le REACTEUR 4.29 bars en remplacement de la
perte de charge prcdente. De mme, la perte de charge dans lchangeur HTX4, nglige
dans un premier temps, vaut dsormais 1,07bars.
Remarque : la perte de charge augmente dans le racteur car le dbit y transitant a
augment. Il serait possible dautomatiser cette modification en spcifiant un coefficient de
perte de charge au lieu de la perte de charge elle-mme.
Nous allons maintenant crer le circuit de recyclage. Ajoutez un sparateur de flux SPLIT
(modle FSPLIT) au dessus du ballon FLASH et reliez FLASHVAP SPLIT. Crez ensuite
un flux de recyclage RECYC entre le sparateur SPLIT et le mlangeur MIX (ne connectez
pas ce flux au mlangeur maintenant !) ainsi quun flux PURGE en seconde sortie du
sparateur (flche bleue car la seconde sortie est facultative). La fraction de purge est de 5%.
Vous pouvez lancer le calcul.
Le recyclage que nous allons maintenant crer introduira dans le programme de rsolution une
procdure itrative puisque le flux gazeux RECYC (obtenu aprs raction et condensation)
Septembre 2012
16
sera ncessaire pour calculer la composition du gaz dentre dans le racteur. Afin de veiller
la bonne convergence de cette procdure itrative, nous allons faire une coupure du flux
FLASHIN (le flux d'entre du sparateur liquide vapeur). Nous pourrions appliquer
lalgorithme de Mottard pour confirmer ce choix (cfr. cours thorique).
Nous fixons une valeur initiale pour les caractristiques (T, P, composition) du flux coup. Le
simulateur calculera, partir de cette valeur initiale, les appareils suivants : FLASH, SPLIT,
MIX, COMP3, HTX3, REACTEUR, HTX4. A la fin de ce calcul, de nouvelles valeurs auront
t trouves pour le flux initial FLASHIN. Elles serviront de nouvelles valeurs initiales pour
ce flux et le processus se reproduira en boucle jusqu ce que la diffrence entre la valeur
calcule et la valeur initiale dune itration diminue et passe sous un seuil appel critre de
convergence.
Pour initialiser FLASHIN, nous partons des caractristiques de ce flux calcules avant de
fermer la boucle. Alternativement, vous pouvez repartir des rsultats suivants pour
linitialisation du flux de coupe :
Maintenant que le flux de coupe est initialis, vous pouvez connecter le flux RECYC au
mlangeur MIX comme 2me entre. Dans Convergence, Tear, choisissez comme flux de
coupe le flux FLASHIN afin que le calculateur sache sur quel critre se baser pour valuer si
le problme converge. Lancer le calcul. Si des erreurs persistent, il peut tre utile de
rinitialiser la simulation avant de lancer le calcul (bouton reinitialize dans la barre doutils
Home, Run), cela permet deffacer les valeurs initiales que vous navez pas spcifies mais
quAspen Plus enregistre automatiquement chaque run.
En analysant les rsultats on peut remarquer que l'unit FLASH2 permet dsormais
d'atteindre une puret de plus de 99,9 %. La quantit dammoniac produite au final augmente
galement.
Pour terminer, afin daffiner terme lvaluation du cot du procd, nous allons modliser
les changeurs de chaleur HTX1 et HTX2 par des changeurs de chaleur double flux. Pour
ce faire, supprimer ces deux changeurs et remplacez-les par des changeurs de type HeatX.
Renommez ces deux nouveaux changeurs comme prcdemment HTX1 et HTX2, reliez-les
avec les diffrents flux existants qui constituent les flux chauds de lchangeur (Attention
bien relier le flux chaud lentre du flux chaud dans lchangeur, et non celle du flux
froid !).
Il faut maintenant connecter un utilitaire froid ces changeurs. Crez quatre nouveaux flux
(1, 2, 3 et 4) qui assureront le passage de leau de refroidissement au sein de lchangeur,
comme reprsent la figure 10.
Septembre 2012
17
SPLIT
PURGE
F1
COMP1
COMP2
1
FLASHVAP
RECYC
F3
HTX4
F2
HTX1
F4
FLASH
F5
HTX2
FLASHIN
MIX
REOUTG
FLASH2G
FLASHLIQ
COMP3
F6
REIN
REACTEUR
FLASH2
F7
HTX3
FLASH2L
REOUTL
Figure 10 : Schma du procd avec boucle de recyclage

Dans Components, ajoutez la substance WATER. Les paramtres binaires proposs par
Aspen pour les mlanges avec de leau peuvent tre accepts tels quels. Fixer la fraction
molaire en eau 0 dans le flux F1. Ensuite, spcifiez les flux dentres 1 et 3 des changeurs
de chaleur en considrant un dbit de 1kmol/sec deau pure 20C, 2bar.
Ensuite, dans HTX1 et HTX2, SETUP, spcifiez les tempratures de sortie des flux chauds
40C pour chaque changeur et la diffrence de temprature minimale 10C (pour savoir
quoi correspond cette temprature, consultez laide dAspen Plus en cliquant lendroit o
vous encodez cette temprature et en appuyant ensuite sur la touche F1 de votre clavier).
Dans longlet Block Options de ces changeurs, indiquez que le modle thermodynamique
choisi pour leau de refroidissement est idal. Cela permet de ne pas devoir utiliser lquation
de RK-Soave pour le calcul du ct froid des changeurs de chaleur. Lancez le calcul.
Enregistrez le procd une fois converg, il servira de base aux prochains exercices.
Application
1. Comment calcule-t-on un rendement de raction chimique ? Dans le cas prsent, quel
est le rendement obtenu dans le racteur la fin de lexercice 3 ? Que devient-il
aprs implmentation des amliorations dcrites dans lexercice 4 et comment cela
sexplique-t-il ?
2. Si on calculait de la mme faon que pour le racteur un rendement de procd, quelle
en serait lvolution entre lexercice 3 et lexercice 4 ? Discutez brivement.
3. La convergence et la vitesse de convergence (nombre ditrations) dpendent-elles fort
des valeurs initiales de dbit encodes dans FLASHIN ? Observez-vous une
diffrence selon que vous approchez la vraie valeur du dbit par une valeur suprieure
ou infrieure, proche ou loigne ?
4. Alors que le rapport H2/N2 est de 3 l'alimentation, pourquoi diffre-t-il lgrement de
cette valeur dans la boucle alors que la raction est stchiomtrique?
5. Que se passe-t-il si on change lgrement le rapport H2/N2 l'alimentation dbit total
et taux de purge fixs ? Reprsentez graphiquement linfluence de ce rapport sur la
production de NH3 et sur le travail du compresseur COMP 3 et discutez-en.
Septembre 2012
18
Exercice 5 : Spcifications de design

La rsolution squentielle modulaire implique que les quipements sont calculs en
considrant que leurs variables dentre sont connues, ce qui implique que les flux entrant
dans le procd sont des constantes, de mme que les paramtres des units. Or, nous avons
vu que lorsquil faut introduire une coupure dans une boucle, une unit de promotion de
convergence est dfinie. Son but est de raliser des simulations successives de la squence de
calcul jusqu convergence, cela en modifiant les variables dentre de la boucle ouverte.
Dans le mme ordre dide, il est possible de provoquer des changements de valeurs de
certaines variables dentre dun procd en utilisant des units mathmatiques dun autre
type : les units de contrle. De telles units sont capables de faire varier une ou plusieurs
variables dentre de linstallation en vue de satisfaire une consigne (cas des Design Spec) ou
doptimiser une fonction objectif complexe (cas de lOptimization), ces fonctions pouvant
dpendre de variables de sortie ou de variables internes du procd.
350000.0
Puissance fournie COMP3 en fonction du taux de purge

WCOMP
0.05
0.06
0.07
0.08
0.09
0.1
0.11
0.12
300000.0
250000.0
Puissance (WATT)
DEB
0.04
0.05
0.06
0.1
PURGE FRAC
0.07
0.08
0.09
0.1
0.11
PURGE FRAC
Concentration en inertes en fonction du taux de purge
0.08
CH4+AR
dinertes
0.06
0.04
Dbit de REIN en fonction du taux de purge
Taux
Dbit (KMOL/HR)
6500.0 7000.0 7500.0 8000.0 8500.0
Lalgorithme utilis pour vrifier les spcifications de design doit effectuer la rsolution dune
ou plusieurs quations algbriques. Cette procdure ncessite une rsolution itrative du
problme pos lunit de promotion de convergence. Nous allons dans la suite illustrer une
telle application en rsolvant une srie de petits problmes de complexit croissante. Mais
examinons dabord par une tude paramtrique linfluence du taux de purge (c.--d. par le
taux de division fix en paramtre de lunit SPLIT) sur certaines variables du procd :
0.04
0.05
0.06
0.07 0.08 0.09

PURGE FRAC
0.1
0.11
0.12
Figure 11 : Evolution de diffrents paramtres en fonction du taux de purge

Septembre 2012
19
0.12
Une augmentation du taux de purge diminue le dbit recycl, ce qui diminue la puissance
fournir au compresseur. Si le dbit purg diminue, les inertes ont tendance saccumuler dans
la boucle, ce que montre bien le troisime diagramme. Lvolution des variables est monotone
et progressive. Elle nest pas linaire, mais on se rend compte quun algorithme bas sur les
linarisations successives, comme la mthode de la corde, a de bonnes chances de trouver
rapidement la solution.
Comme premier exemple, nous supposerons que la puissance du compresseur de recyclage est
limite, ce qui impose implicitement une valeur maximale au dbit de recyclage. Cependant,
celui-ci est fix indirectement par le taux de purge. Nous voudrons donc fixer la valeur dune
variable de sortie du procd (la puissance du compresseur COMP3) en fonction dun
paramtre influenant la solution du problme doptimisation (le taux de purge).
Ensuite, nous fixerons la teneur en mthane lentre du racteur en jouant sur le taux de
purge (en librant la contrainte sur la puissance du compresseur : en effet, on ne peut pas
vrifier simultanment deux spcifications en ne jouant que sur un seul degr de libert). Ceci
consiste fixer la valeur dune variable interne du procd en jouant sur un paramtre dunit.
Nous fixerons ensuite la teneur totale en inertes (mthane + argon) en jouant toujours sur le
taux de purge.
Pour terminer, nous rsoudrons un problme impliquant deux spcifications et deux degrs de
libert : nous ajusterons le taux de purge et le dbit dalimentation en gaz de synthse pour
obtenir un taux dinertes fix lentre du racteur et imposer la production dammoniac. Ces
problmes seront rsolus en partant de ltat du procd obtenu la fin de lExercice 4.
Lors du calcul itratif, le programme recherchera la convergence en suivant ces tapes :
- Donner une valeur initiale au taux de purge, par exemple 5% (la premire valeur
choisie par le programme est celle fixe par lutilisateur);
- Effectuer le calcul du procd ;
- Relever la valeur de la variable que lon veut spcifier ;
- Si la solution correspond la spcification ( une tolrance prs), accepter la solution ;
- Sinon, modifier le taux de purge et recommencer.
Problme 1 : ajuster le taux de purge pour fixer la puissance du compresseur

La puissance du compresseur de recyclage est calcule par lunit COMP3. La solution
converge dpend des paramtres de toutes les units de la boucle, et en particulier du taux de
division du sparateur SPLIT. Supposons que nous souhaitons fixer la puissance du
compresseur de recyclage 270kW.
Pour imposer une spcification quant au design des quipements (ou tout autre paramtre de
la simulation) il est possible de crer une design spec dans flowsheeting options. Dfinissez
dans longlet define une variable que vous nommez W3 comme tant la puissance (Net Work)
du compresseur COMP3. Dans longlet spec, nous souhaitons fixer la variable W3
270000W avec une tolrance de 100W. Pour atteindre cette spcification, nous ferons donc
varier le taux de purge. Dans longlet vary, spcifiez le taux de purge comme variable
manipule, ainsi que ses bornes (de 4 12%) et si vous le souhaitez, la taille pour le pas de
convergence (0,5% par exemple).
Septembre 2012
20
Dans convergence, conv options, default, sequencing, prcisez que la spcification sera
calcule lextrieur de la boucle (la spcification est une variable de sortie du procd, c.-d. que sa valeur na aucune influence sur le calcul de la boucle). Rinitialisez la simulation et
lancez le calcul. Dans les rsultats du second solveur (convergence, convergence, solver2),
vous pouvez consulter lhistorique de la variable manipule au cours des itrations. Vrifiez
dans Design Spec, Results, que la puissance du compresseur de recyclage est bien dans les
limites de tolrance de la valeur demande.
Problme 2 : ajuster le taux de purge pour fixer la teneur en mthane lentre

du racteur
Nous allons prsent fixer une autre spcification : la teneur en mthane lentre du
racteur. Pour ce faire, dans la design spec, nous dfinirons la fraction molaire en mthane C1
dans le flux REIN et nous la fixerons 3% (Tolrance : 0,1%). Faites varier le taux de purge
dans les mmes limites que pour le problme 1 et avec le mme pas. La spcification reste
calcule lextrieur de la boucle. Noubliez pas de cacher ou deffacer la design spec du
Problme 1 avant de lancer votre calcul.
Problme 3 : ajuster le taux de purge pour fixer la teneur en inertes lentre du

racteur
Nous allons prsent fixer le taux dinertes total (CH4 + Ar). Pour ce faire, dfinissez une
nouvelle variable Ar dans la design spec : la fraction molaire dargon du flux REIN. Nous
fixerons le taux dinertes 9%. Cest donc la somme des deux variables Ar + C1 qui doit
valoir 0,09. La tolrance est de nouveau de 0,1%, le nesting reste le mme.
Problme 4 : ajuster le taux de purge et le dbit dalimentation pour fixer la

teneur en inertes lentre du racteur et la production dammoniac
Nous allons ici devoir dfinir une seconde spcification en plus de la spcification prcdente.
Elle sera calcule en mme temps. Dans longlet define, dfinissez le dbit molaire
dammoniaque NH3PRO du flux FLASH2L. Dans longlet spec, fixez la valeur de la
spcification 750 kmol/h avec une tolrance de 1. Dans longlet vary, spcifiez le dbit
molaire dalimentation (F1) comme tant la variable manipule, ainsi que ses bornes (1500
10000) mais laissez Aspen libre de choisir le pas le plus adapt. Dans convergence, conv
options, default, sequencing, choisissez loption Outside simultaneous pour le nesting des
spcifications de design maintenant que vous en avez deux. Elles seront ainsi calcules en
mme temps, tout en restant en dehors de la boucle.
Rsolution des erreurs

En cas de problmes de convergence, il sagit avant tout de bien observer les messages
derreurs qui se retrouvent sur le control panel (accessible depuis le bouton de raccourci en
haut de lcran). La plupart des erreurs proviennent de distractions lors de lencodage des
donnes du problme. Cependant, il arrive que dautres erreurs empchent le calcul de
converger. Ci-dessous, quelques astuces vous permettront de dbuguer un problme
classique :
Septembre 2012
21
Nhsitez pas aller observer les rsultats obtenus pour les flux et les blocs, mme avec
erreurs, pour voir vers quel type de solution le problme converge. Si un flux a un dbit
beaucoup plus grand que ce quoi vous vous attendiez, il peut par exemple y avoir un
problme dans la purge, ou bien dans une temprature de condensation,
Nhsitez pas non plus aller voir les rsultats et historiques des solveurs dans
Convergence/convergence. Quelles variables nont pas converg, pour quelles variables
est-ce que les erreurs sont les plus grandes, ?
Converge block solverX not converged in X iterations : Cela arrive que le nombre
ditrations utilis par dfaut ne suffise pas. La plupart du temps, on le voit en observant
les rsultats de convergence (Tear ou Design Spec). Si lerreur se rapproche de 0, mais
sans y arriver cause dun nombre insuffisant ditrations, il est possible daugmenter le
nombre limite ditrations du solveur concern dans Convergence/Conv options/Methods.
Si lerreur ne se rapproche pas de 0, elle est probablement due un problme mal
conditionn. Cependant, rien ne sert daugmenter de trop le nombre ditrations ! Un
problme ne convergeant pas en 30 itrations convergera rarement en 40 itrations.
Either solution outside bounds or spec function is not monotonic : Vrifiez le choix du
nesting pour la spcification de design. La spcification est-elle calcule dans la boucle
ou en dehors ? Vrifiez les valeurs de la design spec au cours du calcul. Si elle se
rapproche systmatiquement dune limite que vous avez fixe au domaine de variation, la
solution se trouve peut-tre au-del de cette limite. Adaptez le domaine de variation.
Loose tear stream tolerance : vous pouvez diminuer la tolrance du flux de coupe,
mais cela namliorera le problme que si vos donnes initiales sont trs prcises !
De manire gnrale, Aspen Plus est trs sensible aux variations des conditions initiales
du problme. Une fois le calcul converg, mme avec des erreurs (pour autant que les
rsultats paraissent crdibles), vous pouvez copier les rsultats obtenus pour le flux de
coupe (T, P, Dbit, Composition) et les coller dans la partie Input de ce flux pour quils
servent destimation lors du prochain calcul. Cela facilitera la convergence. Vous pouvez
par exemple faire un copi-coll des rsultats obtenus pour le flux Flashin la fin de
lExercice 4 (aprs convergence, clic droit sur le flux, Reconcile).
Il est parfois intressant de lancer le calcul sans rinitialiser la simulation. Cela permet
dutiliser des valeurs initialises lors dun run prcdent. Cependant, les rsultats
varieront sensiblement selon lhistorique des calculs prcdents.
Une simulation peut trs bien comporter des erreurs ou des mises en garde, et cependant
livrer des rsultats exploitables et cohrents ! Cela dpend de la localisation des erreurs,
ainsi que de leur impact sur les rsultats, et cest l quintervient votre sens critique pour
savoir quels rsultats sont acceptables, et lesquels ne le sont pas.
Septembre 2012
22
Application
6. Ralisez un tableau rcapitulatif des rponses aux Problmes 1 4 reprenant chaque
fois la valeur de la fraction de purge, ainsi que la valeur obtenue pour la spcification
de design. Indiquez galement les valeurs obtenues dans le cas de base avant
implmentation des spcifications.
7. Discutez des rsultats prsents la question 1 : sont-ils en accord avec les graphiques
de la figure 11 prsents dans lintroduction (vous navez pas besoin de refaire ltude
paramtrique)?
8. Donnez larchitecture de la squence de calcul (voir le Control Panel) propose par

Aspen pour chaque problme et dans le cas initial. Pour chaque cas, prcisez le rle de
chaque solveur.
9. Le logiciel converge-t-il vers une solution lorsque vous choisissez loption Inside
pour le nesting des spcifications de design dans le Problme 2 ? Pourquoi ? Illustrer
votre rponse laide de larchitecture de rsolution.
Septembre 2012
23
Exercice 6 : Recyclages multiples

En examinant les rsultats avant ltude des spcifications de design (simulation obtenue la
fin de lExercice 4), on peut constater que la chaleur fournir dans l'changeur HTX3 est
infrieure celle qui doit tre limine dans l'changeur HTX4. De plus, on constatera
aisment que les niveaux de temprature des flux en prsence sont totalement compatibles en
vue de raliser l'change en interne au procd (c.--d. sans faire appel un utilitaire). Pour
cela, il suffit de modifier lgrement le flowsheet de manire prchauffer l'alimentation du
racteur partir de l'excdent thermique de sa sortie. Cela est reprsent la figure 12. Le
reste du refroidissement pourra tre accompli en faisant appel des utilitaires froids : les
niveaux de temprature sont tels qu'on pourra mme se permettre de produire de la vapeur.
En procdant de la sorte, on ajoute un recyclage de chaleur dans le procd, et on introduit
une boucle supplmentaire. En effet, l'ensemble {RECUP, REACTEUR, VAPO} (voir figure
12) constitue un bloc irrductible : aucun des appareils ne peut tre calcul de prime abord,
sans connatre les sorties d'un autre appareil. On doit donc tenir compte d'une coupure
supplmentaire dans le procd, rechercher par la mthode de Mottard.
SPLIT
PURGE
FLASHVAP
RECYC
FLASH
HTX4
FLASHIN
F5
MIX
FLASHLIQ
F9
F6
5
COMP3
F7
REOUTG
VAPO
REIN
REACTEUR
RECUP
REOUTL
F8
Figure 12 : Recyclages de matire et de chaleur

Cette coupure peut tre fixe l'entre du racteur. En effet, cet endroit, le dbit et la
composition du gaz seront connus au moment d'entamer les calculs du bloc irrductible. Seul
l'tat thermique (enthalpie ou temprature) du fluide sera inconnu, et la coupure ne portera
donc que sur une seule variable (coupure thermique). Nous aurons donc deux boucles de
calcul imbriques lune dans lautre. La premire boucle rsoudra le systme compos des
units REACTEUR, VAPO, RECUP. La deuxime boucle aura pour objectif de faire
converger le calcul centr sur le flux de coupe FLASHIN.
Nous repartons du procd tel quil tait la fin de lExercice 4. Dans le cas o cela na pas
t fait, commencez par rconcilier le flux de coupe FLASHIN (clic droit sur le flux,
reconcile). Cela facilitera le calcul de lchangeur RECUP, car les valeurs initiales de dbit et
de composition donnes FLASHIN dans lExercice 4 sont trop loignes de la valeur finale
pour obtenir une bonne convergence une fois la deuxime boucle implmente. Commencez
par supprimer lchangeur HTX3 pour le remplacer par un changeur de type HeatX que vous
nommerez RECUP. Reliez les flux F7 et REIN cet changeur. Faites attention, ceux-ci
Septembre 2012
24
doivent constituer la partie froide de cet changeur ! Vous pouvez faire pivoter licne de
lchangeur dun quart de tour sur la gauche pour un dessin plus simple.
Nous allons maintenant modliser lunit VAPO destine produire de la vapeur deau 40
bars. En effet, la temprature en sortie du racteur est tellement leve quelle nous permet de
produire de la vapeur en vue, par exemple, dalimenter un cycle de Rankine (cycle
thermodynamique permettant de produire du travail grce un fluide caloporteur changeant
de phase) ou mme d'exporter cette chaleur vers un autre endroit du site. Pour modliser cette
unit, slectionnez une nouvelle fois le modle dchangeur de chaleur HeatX que vous
placerez entre le REACTEUR et RECUP. Nommez cet changeur VAPO. Dconnectez
REOUTG de HTX4 (clic droit, reconnect destination) pour le reliez VAPO (flux chaud).
Crez un flux de matire F8 entre VAPO et RECUP. Reliez RECUP et HTX4 par un flux de
matire F9. Finalement, insrez un utilitaire froid (crez deux flux, 5 et 6, de part et dautre de
VAPO). Spcifiez le flux 5 comme suit : Eau pure ; 20C ; 40bar ; 400kmol/h.
Spcifiez ensuite RECUP. La temprature de sortie du flux froid (flux entrant dans le
racteur) est fixe 380C, la diffrence de temprature minimale 10C et la perte de charge
1,07 bar de chaque ct. Enfin, pour lchangeur VAPO, la temprature de sortie du flux
chaud F8 est fixe 450C, la diffrence de temprature minimale 10C et la perte de
charge 1.07 bar pour le flux chaud uniquement. Comme pour les changeurs de chaleur
HTX1 et HTX2, dans Block Options, indiquez IDEAL comme mthode thermodynamique de
lutilitaire froid pour VAPO.
En thorie, la coupure thermique pourrait tre nglige pour ce calcul vu que la temprature
du flux REIN est impose par la consigne donne RECUP. Cependant, le problme ne
converge pas facilement car aucune valeur initiale de dbit na t donne ni REIN ni
REOUTG ou F8. Nous allons donc donner une valeur initiale REIN en nous basant sur le
flux F7 calcul dans lExercice 4. La pression vaudra 278,93bars et la temprature sera de
380C. Le dbit vaut 8000 kmol/h et les fractions molaires sont :
Ar
CH4
N2
0,042
0,040
0,221
H2
NH3
0,657
0,040
Nous allons maintenant dfinir la mthode de convergence employe pour le calcul des deux
boucles. Dans Convergence, Convergence, dfinissez une nouvelle boucle de calcul. Cliquez
sur New et nommez-la BOUCLE2. Choisissez Broyden comme type de mthode (cfr. cours
thorique). Indiquez dans la section Tear Streams, le flux coup, savoir FLASHIN. De
mme, dfinissez une seconde boucle de calcul BOUCLE1 et choisissez Broyden comme
type de mthode. La coupure se fait sur le flux REIN et il faut ici prciser dans State
Variables le type enthalpy (coupure thermique). Gardez les valeurs de tolrance par dfaut.
Remarque : Il aurait t possible de dfinir galement cette coupure dans Convergence, Tear,
Specifications comme cela avait t fait auparavant. Nanmoins, dfinir chaque boucle
sparment permet dventuellement utiliser diffrentes mthodes de rsolution (Wegstein,
Broyden, Newton, ) pour les boucles ainsi que de dfinir lordre dans lequel nous
souhaitons quAspen Plus effectue le calcul.
Dans notre cas, la boucle interne du procd (recyclage thermique) devra correspondre la
boucle interne du calcul afin de sassurer dune bonne convergence. Pour indiquer au
Septembre 2012
25
programme dans quel ordre de boucle il doit effectuer ses itrations, dans Convergence, Conv
Order, slectionnez BOUCLE1 puis BOUCLE2. En effet, lajout seffectue de la boucle la
plus encapsule la boucle la moins encapsule.
Maintenant que nous avons dfinis ces diffrentes units, nous lanons un premier calcul. En
observant les rsultats relatifs lunit RECUP, nous constatons que cet changeur de chaleur
ncessite une surface dchange de 456 m (Thermal results/Exchanger details). Nous
pouvons maintenant changer le mode de lunit RECUP pour passer en mode simulation (au
lieu du mode design). Nous fixons la surface dchange 456 m. Nous faisons de mme pour
VAPO o la surface est fixe 26m. Lancez le calcul et sauvegardez vos rsultats.
NB :
- Si vous avez des problmes de convergence de type bilan de masse (bien que cette
erreur nimpacte le plus souvent pas significativement les rsultats), vous pouvez
spcifier un critre de tolrance plus exigeant (par exemple 10-5) dans la BOUCLE1
pour viter quun bilan imprcis dans cette boucle nait des rpercussions nfastes sur
la convergence de la BOUCLE2. Vous pouvez faire de mme avec la BOUCLE2.
- Si vos changeurs de chaleur ne convergent pas, vous pouvez lancer la convergence
en mode pas pas et observer les rsultats aprs le calcul de lchangeur. Les
tempratures dentre et de sortie sont-elles raisonnables ? Si non, pourquoi ? Une
erreur frquente est un dbit dun des deux flux trop faible ou trop important cause
dune erreur dunit lors de lencodage des donnes. Ce procd de convergence
pas pas peut galement tre utile pour dautres types dunits.
Application
10. Quelle est la temprature minimale que peut atteindre F8 si on veut prchauffer
lalimentation du racteur 400C avec les changeurs de chaleur donns ? Quelle est
la quantit de vapeur quon produit ainsi dans lunit VAPO ?
11. Si lon souhaite garder lentre dans le racteur 400C, est-il possible de produire
plus de vapeur si celle-ci est saturation (non surchauffe) ? Si oui, quelle quantit ?
Conseil : utiliser une spcification de design pour fixer la temprature de REIN et
revenez en mode design pour VAPO.
12. Quel travail peut-on rcuprer en dtendant cette vapeur dans une turbine, la pression
de sortie tant fixe 50 mbar ? Produit-on plus de travail en dtendant la quantit de

vapeur surchauffe de la question 10 ou la quantit maximale de vapeur
sature calcule la question 11 ?
Remarque : Pour insrer une turbine, insrez un compresseur et choisissez le mode turbine.
Afin dautoriser une phase liquide dans la turbine en plus de la phase vapeur, dans longlet
convergence de la turbine, vous pouvez choisir Vapor-liquid.
13. En rsum, quelles coupures a-t-on dans le procd et de quel type sont-elles ?
Justifiez leur nombre par la mthode de Mottard. Pourrait-on mettre la coupure
matire sur le flux REIN, et serait-ce plus appropri que la configuration actuelle ?
14. Comment varie le besoin nergtique de compression (somme des trois compresseurs)
par tonne de NH3 liquide produit en fonction du taux de purge (entre 6 et 12%) ?
Reprsentez graphiquement votre rponse en GJ/Tonne NH3.
Septembre 2012
26
Exercice 7 : Optimisation
Introduction
La mise au point dun procd de fabrication implique de prendre des dcisions de deux
types :
- Le choix de la structure du procd, cest--dire du type doprations qui vont
senchaner et lordre de ces oprations ;
- Le choix des paramtres opratoires (temprature, pression, dbits).
Lobjectif est bien sr de concevoir un ensemble dappareils permettant de fabriquer le
produit dsir en respectant la demande et les spcifications du march, dans des conditions
de scurit convenables, et un cot minimum. Les premires dcisions font intervenir des
variables discrtes, que lon peut reprsenter par des nombres entiers (faut-il un, deux ou trois
lits de catalyseur ? Combien dtages de compression ? Lagencement A ou lagencement
B ?). Le second type de problme fait intervenir uniquement des variables continues. Les
techniques doptimisation interviennent dans les dcisions prendre. Nous les appliquerons
ici loptimisation de variables continues uniquement, les problmes variables discrtes
tant dapproche nettement plus complexe et dpassant le cadre dun cours dintroduction.
Dans la boucle de synthse dammoniac tudie jusqu prsent, un poste de cot important
est la compression du gaz de synthse. Cest pourquoi nous allons essayer doptimiser les
paramtres de fonctionnement du train de compression dcrit par le schma suivant :
Figure 13 : schma du train de compression du gaz de synthse

Lorsquon se pose un problme doptimisation, il faut respecter la dmarche suivante :
-
Dfinir un modle de procd ;

Dfinir la fonction objectif optimiser ;
Identifier les variables de commande, en fixer les valeurs initiales et les bornes de
variation ;
Rechercher les contraintes dgalit et dingalit satisfaire ;
Septembre 2012
27
Paramtrer loptimiseur ;
Dfinir une squence de calcul valuant le procd en fonction des variables
indpendantes ;
Rsoudre.
Modle de procd
Nous nous basons sur le modle de procd tel que ralis au cours des exercices prcdents
et partiellement reprsent la figure 13. Cependant, comme nous ferons varier la surface des
changeurs de chaleur lors de loptimisation, nous devons utiliser cette variable dans notre
modle. Nous fixerons les paramtres des compresseurs et la valeur initiale des surfaces
dchange en mode simulation comme dcrit ci-dessous :
-
COMP1 : calcul pression de sortie impose (90 bar);

HTX1 : calcul en mode simulation aire et diffrence de temprature minimale fixes
(28 m et 10C) ;
COMP2 : calcul pression de sortie impose (270 bar);
HTX2 : calcul en mode simulation aire et diffrence de temprature minimale
fixes (52 m et 10C) ;
COMP3 : calcul pression de sortie impose (280 bar).
Enfin, pour les flux deau 1 et 3, nous changeons les valeurs initiales du dbit 2 kmol/s.
Fonction objectif
Dans notre cas, la fonction objectif sera le cot annuel dexploitation de ce train de
compression. Il sagira de le minimiser. Ce cot dexploitation comporte deux termes :
-
Un terme proportionnel correspondant aux consommations deau et dnergie ;

Un terme fixe correspondant lamortissement et lentretien du matriel, que nous
supposerons gal une fraction du cot dinstallation des quipements.
Concernant les cots variables (OPEX, operational expenditure), il nous faut compter le cot
li leau de refroidissement ainsi que celui li la demande nergtique du procd. Nous
supposerons que le cot de leau de refroidissement est de 0.02 /m. Ceci correspond au
pompage et la purification deau de rivire ou lutilisation dune boucle ferme et dune
tour de rfrigration atmosphrique (appoint deau frache pour compenser lvaporation,
traitement anticorrosion et anti-algues, pompage). Il faudra tenir compte du fait que leau est
disponible 20C et ne peut tre retourne une temprature suprieure 45C. Quant
lnergie, nous supposerons que son cot est de 0.07 /kWh. Linstallation est suppose
fonctionner 8000 h/an.
Pour valuer les cots fixes (CAPEX, capital expenditure), il faudra disposer dune
estimation du cot de linstallation. Les mthodes habituelles reposent sur la technique de
factorisation. On dispose dabaques qui donnent le cot de base dun appareil standard en
fonction dun paramtre de taille (surface pour un changeur de chaleur par exemple). Des
facteurs de correction multiplicatifs tiennent compte du choix des matriaux, de linfluence de
la pression ou de tempratures leves, ou de dtails de conception (par exemple, pour les
changeurs de chaleur : plaques tubulaires fixes, ou ttes flottantes, ou faisceau de tubes en
U). Le prix ainsi calcul correspond la fourniture de lappareil. Il doit tenir compte des frais
dinstallation (fondations, charpente, montage, tuyauterie, instrumentation, isolation
thermique, peinture) ce qui se fait galement au moyen de facteurs proportionnels.
Septembre 2012
28
Pour les compresseurs COMP1, COMP2 et COMP3 la variable de base de la corrlation de

cot est la puissance de la machine W, exprime en kW. Le prix des fournitures des
quipements C1 et le prix dachat du moteur lectrique dentranement C2 seront donns
respectivement par5 :
C1(k) = 11.0 W0.6 et C2(k) = 0.18 W
On obtiendra le cot de lquipement install en multipliant ce prix de base par le facteur
dinstallation de 2.40 :
C3(k) = 2.40 (C1 + C2)
Pour les changeurs de chaleur HTX1 et HTX2, le paramtre de taille est la surface dchange
A (m) obtenue lors du calcul en mode design. Elle est la base de la corrlation de cot. Le
prix des fournitures des quipements sera donn par :
C4(k) = 1.10 A0.75
On obtiendra le cot de lquipement install en multipliant ce prix de base par le facteur
dinstallation de 2.85 :
C5(k) = 2.85 C4
Enfin, il faudrait prendre en compte lindexation et les conditions conomiques la date de
construction pour convertir le cot de rfrence du matriel install en cot damortissement
annuel. Plus de dtails seront donns au cours de lavant-projet dusine en dernire anne.
Dans le cadre de cet exercice, nous supposerons que lamortissement et les frais annuels
dentretien reprsentent un septime du cot du matriel install.
C6(k) = (C3+C5)/7
Reprenez le flowsheet obtenu la fin de lexercice 6 (avant les applications). Avant dcrire
la routine permettant de calculer la valeur de notre fonction objectif dans Aspen, nous allons
ajouter un mlangeur MIXW. Ce mlangeur aura comme flux entrants trois flux mcaniques
W1, W2 et W3 ayant pour sources les trois compresseurs COMP1, COMP2 et COMP3 et un
flux sortant WTOT. Le rajout de flux mcaniques sur le schma du procd se fait de la mme
manire que pour les flux de matire, mais en prcisant Work streams au lieu de Material
streams en bas gauche de votre cran.
Le calcul de la fonction objectif se fera en ajoutant un bloc de calcul. Dans le menu, choisir
flowsheeting options, calculator, et crer un nouveau calculateur CALC1. Dans cette feuille,
nous allons dfinir les variables du procd intervenant dans le calcul de la fonction objectif.
Les cinq premires variables ci-dessous seront des variables importes, cest--dire des
variables dont la valeur obtenue dans le procd sera importe dans le calculateur. Pour
dfinir une nouvelle variable, cliquez sur new, donnez le nom de la variable, et dfinissez-la.
Noubliez pas dindiquer quil sagit ici de variables importes.
5
Ces corrlations tant approximatives, il est dconseill de les utiliser pour des applications ultrieures
lorsquil sagira dtablir des estimations prcises.
Septembre 2012
29
Le travail (Work-Power) du flux WTOT que nous appellerons galement WTOT ;

Les dbits massiques deau (Mass-flow) des flux 1 et 3 appels EAU1 et EAU2 ;
Les surfaces calcules (Area calc) des changeurs HTX1 et HTX2 que nous
appellerons respectivement A1 et A2 ;
La sixime variable, la valeur de la fonction objectif COST, est une variable exporte, cest-dire que sa valeur doit pouvoir tre exporte vers un flux ou une unit du flowsheet pour
pouvoir tre utilise lors de loptimisation. Or cette variable nappartient aucun flux ou
unit. Pour palier ce problme, nous ajoutons dans le flowsheet, sans le relier aucune autre
unit, un nouveau mlangeur FOBJCALC. Ce mlangeur a un flux de travail entrant FOBJ et
un flux de travail sortant FOBJBIS. Notre variable COST est le travail du flux FOBJ. Le
flux FOBJ devant tre spcifi, donnez-lui la valeur que vous voulez, elle sera de toute faon
modifie aprs le calcul de la variable COST. Vous pouvez maintenant ajouter la variable
COST dans CALC1, en spcifiant bien que sa valeur devra tre exporte vers le flowsheet :
-
La variable COST, travail du flux FOBJ contenant la valeur de la fonction objectif.
Dans la feuille calculate du calculateur, nous crivons la routine de calcul de notre fonction
COST. Celle-ci doit tre crite en langage fortran. Les noms de variable en fortran doivent
rpondre une srie de critres :
-
Pas plus de six caractres ;

Commencer par une lettre, les caractres suivants tant alphanumriques ;
Ne pas commencer par IZ ou ZZ ;
Pas de diffrence entre les majuscules et les minuscules ;
Par dfaut, les variables commenant par A-H ou O-Z sont des variables double
prcision et les variables commenant par I-N sont des variables entires.
Les six premires colonnes sont rserves en fortran et ne peuvent pas accueillir le code
proprement dit. Ncrivez donc rien sur les six premires colonnes de chaque ligne. Dans
longlet calculate du calculateur, choisissez fortran et encodez les relations dcrites
prcdemment (calcul du cot total en k):
ccomp=2.4*(11.0*(WTOT/1000)**0.6+0.18*WTOT/1000)
coutA1=2.85*1.10*A1**0.75
coutA2=2.85*1.10*A2**0.75
cannu=(coutA1+coutA2+ccomp)/7
cener=(0.07/1000)*8000*WTOT/1000
Dm3h=(EAU1+EAU2)/1000
ceau=(0.02/1000)*8000*Dm3h
COST=cannu+cener+ceau
Dans longlet squence du calculateur, choisissez use import/export variables. Le calculateur
effectuera alors son calcul en mme temps que la squence de calcul du procd, selon quil
rencontrera des valeurs importer ou exporter. Vous pouvez lancer le calcul et lire la valeur
du cot dexploitation annuel avant optimisation (un peu moins de 4000k).
Septembre 2012
30
Variables de commande
On peut identifier six variables de commandes qui seront dfinies plus tard dans Aspen lors
du paramtrage de loptimiseur. Ce seront les variables qui seront varies pour minimiser la
fonction objectif :
-
P1, la pression de sortie du compresseur COMP1 ;

P2, la pression de sortie du compresseur COMP2 ;
A1, la surface de lchangeur HTX1 ;
A2, la surface de lchangeur HTX2 ;
D1, le dbit molaire du flux 1 ;
D2, le dbit molaire du flux 3.
Au moment voulu, on fixera les bornes de ces variables de commande.

Contraintes du procd
Idalement, le modle doit permettre dvaluer le fonctionnement du procd pour toute
combinaison des variables de commande. Il peut toutefois exister des combinaisons de
variables qui conduisent un fonctionnement aberrant ou irraliste du procd (par exemple,
le rfrigrant devient plus chaud que le flux de procd refroidir). Cest pour viter cela que
nous avons la possibilit de fixer des contraintes.
Contraintes dgalit
Aucune contrainte dgalit nest requise.
Contraintes dingalit
On peut par exemple supposer que la puissance du troisime tage de compresseur COMP3
soit limite. En effet, le but principal de ce compresseur tant de vaincre les pertes de charges
lies au recyclage, lappareil dont on dispose est dun niveau infrieur aux deux premiers
compresseurs. Sa puissance est donc limite 300 kW. Nous imposerons galement une
borne suprieure pour la temprature de rejet de leau de refroidissement (flux 2 et 4), soit
45C.
Pour ajouter une contrainte, choisissez model analysis tools, constraints, new et dfinissez les
contraintes
- C1 contrainte sur la puissance de COMP3. La variable utilise est la puissance NETWork du flux W3 qui devra tre infrieure 350.000W avec une tolrance de 1000W.
- C2 et C3 contraintes sur la temprature de sortie de leau de refroidissement (flux 2 et
4), infrieure 45C avec une tolrance de 0,1C.
Paramtrage de loptimiseur
Nous allons prsent paramtrer loptimisation des variables de commande. Pour ce faire,
choisissez model analysis tools, optimisation, new et nommez cette optimisation OPTIM1.
Dans longlet define de loptimisation, dfinissez la variable COST (valeur de notre fonction
objectif, qui correspond au travail du flux FOBJ) ainsi que les six variables de commandes.
Les variables P1, P2, A1 et A2 sont des variables de bloc, D1 et D2 des variables de flux.
Choisissez les variables spcifies (variables de commande) et non les variables calcules.
Dans longlet objective and constraints de loptimisation, indiquez la variable optimiser.
Dans notre cas, nous voulons minimiser la variable COST. Slectionnez les trois contraintes
dingalit.
Septembre 2012
31
Dans longlet vary de loptimisation, nous devons lister nouveau les six variables de
commande avec leurs bornes et la taille du pas. Dans le menu droulant de variable number,
choisissez new et indiquez les caractristiques de chaque variable. Rfrez-vous au tableau
rcapitulatif page suivante. Vous pouvez si vous le souhaitez indiquer un pas de variation,
mais ce nest pas exig par Aspen.
Variable
P1 (Pres)
P2 (Pres)
A1 (Area)
A2 (Area)
D1 (mole-flow)
D2 (mole-flow)
Bloc/flux
Comp1
Comp2
HTX1
HTX2
1
3
Unit
bar
bar
m
m
kmol/hr
kmol/hr
Valeur initiale
90
270
28
52
7200 (2kmol/sec)
7200 (2 kmol/sec)
Minimum
40
100
20
20
3000
3000
Maximum
200
350
60
60
15000
15000
Squence de calcul
Afin de raliser loptimisation, nous allons devoir ajouter un optimiseur. Pour ce faire,
slectionnez convergence, convergence, new. Nommez-le OPTIM et choisissez la mthode
complex (voir mthode du simplex cours thorique). Dans longlet optimization de
loptimiseur, slectionnez loptimisation choisie, soit OPTIM1. Dans longlet parameters de
loptimiseur, fixer la tolrance relative 0,01. Dans convergence, conv order, placez cette
unit de convergence en dernire position de la squence, de faon englober les boucles 1 et
2.
Rsolution
Rinitialisez la simulation et lancez le calcul. Vrifiez que vos contraintes sont respectes, et
que la valeur de votre fonction objectif a bien diminu.
Application
1. Quelle est la rpartition des cots annuels, nergtique et lis la consommation deau
(cannu, cener, ceau, calculer en dehors dAspen), avant et aprs optimisation ?
Quelles conclusions pouvez-vous en tirer ?
2. Que se passe-t-il si on dispose pour HTX1 dun changeur dont la surface maximale
est de 40m ?
3. Que se passe-t-il si on augmente 400 kW la puissance maximale de COMP3 ?
4. Dans un tableau, rsumez les valeurs des variables de commande dans les quatre
configurations tudies: cas de base, cas de base aprs optimisation, et les deux cas de
lapplication ci-dessus. La fonction COST volue-t-elle de faon conforme vos
attentes dans ces diffrents cas ? Discutez.
Remarque :
- Vous rencontrerez parfois des erreurs de convergence en changeant les valeurs des
contraintes. En effet, loptimiseur ne fonctionnera pas si la configuration initiale ne
satisfait pas aux contraintes. Dans ce cas, vous de modifier vos valeurs initiales
pour que les contraintes soient respectes.
- Des erreurs naffectant pas la convergence peuvent survenir. Si elles se rglent
delles-mmes lors du calcul, pas besoin de vouloir les supprimer tout prix !
Septembre 2012
32
Exercice 8 : Distillation
La distillation est un des procds de sparation les plus frquemment utiliss dans lindustrie
chimique, que ce soit pour liminer les impurets indsirables dans les alimentations des
procds, ou pour obtenir des produits dune puret correspondant aux besoins du march.
Cette opration ncessite des investissements importants en matriel et une grande
consommation dnergie thermique. Cest pourquoi la conception et la conduite des colonnes
de distillation fait appel des modles de calcul labors.
Une unit de distillation se compose de colonnes assurant un contact intime entre une phase
vapeur et un reflux liquide. Le contact entre les phases est favoris par lutilisation de
plateaux (perfors, calottes ou clapets) ou de garnissages (anneaux de Raschig par
exemple). Le reflux est obtenu en condensant les vapeurs en tte de colonne et en en recyclant
une partie vers la colonne. En pied de colonne, un rebouilleur vaporise une partie du liquide
prlev, ce qui produit la vapeur circulant dans lappareil.
La modlisation classique des colonnes de distillation repose sur le concept dtage
dquilibre ou plateau thorique. On admet que le liquide et la vapeur quittant un plateau
thorique sont en quilibre (mmes pression et temprature, fugacits gales pour toutes les
substances). On tient compte ventuellement dun correctif cette approche idalise en
introduisant un facteur defficacit (coefficient de Murphree par exemple). Un plateau rel
dont lefficacit serait de 70% se comporte comme un tage dquilibre o seulement 70% de
la vapeur parcourant la colonne entre en contact prolong avec 70% de la phase liquide,
permettant ltablissement dun quilibre. Les 30% restant nentrent pas en contact efficace.
Le modle dun tage dquilibre thorique (figure 14) pour un mlange de N constituants
comportera donc 2N+3 variables dtat :
-
1 dbit total de liquide descendant vers le

plateau infrieur ;
N fractions molaires en phase liquide ;
1 temprature de plateau ;
1 dbit total de vapeur montant vers le plateau
suprieur ;
N fractions molaires en phase vapeur.
Figure 14 : modle dtage thorique
Les quations du modle dtage dquilibre seront galement au nombre de 2N+3 :

-
N quations de bilans de matire du plateau i (sortie = entre pour chaque substance) ;

1 quation de bilan thermique (enthalpie entrante = enthalpie sortante) ;
N quations dquilibre liquide-vapeur (fugacit liquide = fugacit vapeur pour chaque
substance) ;
Somme des fractions molaires du liquide = 1 ;
Somme des fractions molaires de la vapeur = 1.
Des variables additionnelles peuvent tre introduites pour dcrire des soutirages (saigne
dune partie du liquide si sur un plateau, ou prlvement dune partie de la vapeur Si) et pour
tenir compte dappoint ou de retrait de chaleur Q. Elles sont compenses par un nombre gal
de spcifications.
Septembre 2012
33
Si nous remarquons quune colonne de distillation peut se reprsenter par une batterie
dtages dquilibre entre lesquels circulent des flux liquide et vapeur, il apparat toutefois
que le calcul de ces ensembles ne peut se traiter aisment par une mthode squentielle
modulaire, combine une technique de promotion de convergence sur des flux coups. En
effet, lapplication dune mthode de dcomposition comme celle de Mottard conduirait
couper un flux entre chaque paire de plateaux adjacents, ce qui ne peut se faire faute dune
bonne connaissance de valeurs initiales et au vu de la complexit de la tche. Il est donc
ncessaire denvisager une rsolution globale du systme dquations. Cependant, si on
considre quune colonne moyenne peut traiter un mlange dune dizaine de substances et
comporter une trentaine de plateaux, on est confront la rsolution de systmes de grande
taille : 30 x (2 x 10 + 3) = 690 quations pour lexemple choisi. Des techniques de rsolution
tenant compte de la structure de la matrice Jacobienne du problme sont donc ncessaires.
Elles reposent sur la constatation que cette matrice est creuse : la plupart des lments sont
nuls, car les variables dun plateau ninterviennent que dans les quations de ce plateau et
dans celles des plateaux adjacents.
A la figure 15 sont illustrs des schmas de procds importants o plusieurs colonnes de
distillation sont fortement couples. Pour faciliter le calcul de ces procds, il est utile de
disposer de programmes de simulation capables de rsoudre les modles de plusieurs
colonnes.
Figure 15 : schmas de procds impliquant des colonnes de distillation

Pour nous familiariser avec la simulation des colonnes de distillation, nous utiliserons
plusieurs modles disponibles dans Aspen. Tout dabord, nous ferons appel au modle
DSTWU, modle simplifi bas sur les corrlations de Winn (calcul du nombre dtages
minimum), Underwood (calcul du reflux minimum) et Gilliland (calcul du reflux requis pour
un nombre dtage fix ou du nombre dtage requis pour un reflux spcifi). Par la suite,
nous utiliserons le modle RADFRAC, qui calcule les quilibres plateau par plateau.
Septembre 2012
34
Problme 1 : train de sparation des hydrocarbures

Notre premire application sera dvaluer un train de sparation dhydrocarbures. Le mlange
sparer sera constitu des cinq hydrocarbures suivants: propane, isobutane, normal-butane,
isopentane et normal-pentane. Nous souhaitons sparer les cinq composants de manire
obtenir cinq produits relativement purs. Pour ce faire, quatre colonnes de distillation seront
ncessaires. Ces colonnes peuvent tre disposes de plusieurs manires et raliser les
sparations suivant des squences diffrentes, comme le montre la figure suivante :
Figure 16 : Possibilits pour un train de sparation des hydrocarbures

La premire squence limine les substances en commenant par la plus lgre. La seconde
commence par la plus lourde. La troisime squence spare dabord le propane, puis spare
les C4 des C5, et enfin partage les C4 et les C5 dans les deux dernires colonnes. Dautres
squences sont possibles et nous en examinerons une par la suite. Le choix de la squence
optimale dpend des volatilits relatives et de la composition de lalimentation. Suivant les
spcifications choisies, la taille des colonnes (nombre dtages dquilibre) et la
consommation dnergie (bouilleur et condenseur) variera fortement. Pour une sparation
donne, une augmentation du taux de reflux permettra de diminuer le nombre de plateaux et
donc la hauteur de la colonne, mais cette augmentation du reflux ncessitera un plus gros
condenseur, un plus gros rebouilleur, un plus grand diamtre de colonne et des
consommations accrues dutilits (vapeur et eau de refroidissement). Un diminution du taux
de reflux impliquera par contre laugmentation du nombre de plateaux et donc de la taille de
la colonne. Tout design est donc un compromis entre les investissements lis la taille des
appareils, et les frais dexploitation lis aux consommations dutilits.
Pour tudier ce procd, dmarrez une nouvelle simulation gnrale avec units mtriques.
Dans Set-up/Specifications, donnez-lui un titre et spcifiez le systme dunits SI-Cbar.
Septembre 2012
35
Insrez les 5 composs suivants dans la liste des composs :

- propane (C3)
- isobutane (IC4)
- n-butane (NC4)
- isopentane (IC5)
- n-pentane (NC5)
Pour ce type de mlange, lquation de RK-Soave dj utilise prcdemment est adquate
pour dcrire le modle thermodynamique. Attention, elle ne figure pas parmi les mthodes les
plus communes dAspen, il faut slectionner toutes les mthodes pour la trouver. Vous
pouvez accepter les paramtres proposs par Aspen tels quels.
Passez maintenant linterface simulation. Dans Set-up/Report options/Streams, cochez la
case adquate pour obtenir les fractions molaires des composs dans les flux. Vous pouvez
dessiner le schma du procder en vous inspirant de la figure 17. Pour ce faire, insrez quatre
colonnes de type DSTWU, nommez-les et connectez-les entre elles par des flux de matire.
Comme nous lavons vu dans lintroduction, ce type de colonne se base sur les corrlations de
Winn pour le calcul du nombre dtages minimum, de Underwood pour le calcul du reflux
minimum et de Gilliland pour le calcul du reflux requis pour un nombre dtage fix ou du
nombre dtage requis pour un reflux spcifi.
Figure 17 : schma du train de sparation des hydrocarbures

Ajoutez les spcifications du flux dalimentation FEED. Le flux dalimentation que nous
souhaitons sparer prsente les caractristiques suivantes :
- Dbit : 100 kmol/h
- Le flux se trouve temprature de bulle
- Pression : 10 bar
Septembre 2012
36
Composition (fractions molaires)

o propane : 0.05
o isobutane : 0.15
o n-butane : 0.25
o isopentane : 0.20
o n-pentane : 0.35
Il ne reste plus qu spcifier le mode de calcul des quatre colonnes. Pour chacune des quatre
colonnes, il faut tout dabord spcifier si on choisit de travailler nombre de plateau ou taux
de reflux fix. Nous travaillerons ici taux de reflux fix de sorte que le programme calculera
le nombre de plateaux ncessaires.
Une valeur positive dans la case reflux ratio signifie que lon fixe un taux de reflux, tandis
quune valeur ngative signifie que lon fixe un rapport entre le taux de reflux souhait et le
taux de reflux minimum (R/Rmin). Nous opterons pour cette seconde option. Les rapports de
taux de reflux sont indiqus dans la table ci-dessous, ainsi que les composs cls spcifier
pour le calcul. Ces constituants cl sont les deux substances qui forment la charnire entre les
fractions lourdes et les lgres. On a ainsi
- Le composant cl lger (constituant le plus lourd que lon retrouve en quantit
apprciable dans le distillat) ainsi que la fraction du dbit dalimentation de cette
substance qui se retrouve dans le distillat ;
- Le composant cl lourd (constituant le plus lger que lon retrouve en quantit
apprciable dans le rsidu) ainsi que la fraction du dbit dalimentation de cette
substance qui se retrouve dans le distillat.
Spcifications des colonnes du train de sparation
Colonne Cl lger
Fraction Cl lourd Fraction R/Rmin
COL1
n-butane
0.99
isopentane 0.01
-1.25
COL2
propane
0.99
isobutane 0.01
-1.5
COL3
isobutane 0.99
n-butane
0.01
-1.3
COL4
isopentane 0.99
n-pentane 0.01
-1.3
Les fractions de composs lger et lourd vont dterminer la puret obtenue dans la colonne,
tandis que la taille de la colonne sera lie au taux de reflux. Enfin, nous fixerons la pression
au condenseur et au rebouilleur 8 bar pour chaque colonne. Pour chacune des colonnes, nous
choisirons loption condenseur total. Vous pouvez lancer le calcul et vrifier que vous
parvenez bien sparer les diffrents hydrocarbures.
Application :
5. Quelle est linfluence de la pression de fonctionnement de la premire colonne ?
Conseilleriez-vous une pression plus basse (6 bar par exemple) ou plus leve (10 bar
par exemple) ? Justifiez.
6. Quelle est linfluence du taux de reflux sur la gomtrie de la colonne ainsi que sur la
consommation nergtique de celle-ci ? Illustrez votre propos en calculant la colonne
COL2 avec 3 taux de reflux diffrents. Quel taux de reflux conseilleriez-vous ?
Pourquoi ?
Septembre 2012
37
Problme 2 : calcul dtaill de la premire colonne de sparation

Pour suivre, nous allons valuer plus en dtail le fonctionnement de la premire colonne du
train de sparation. Le modle employ pour le calcul rigoureux dune suite dtages
dquilibre constitue un lment essentiel de tout programme de simulation de procd. Ce
modle se nomme RADFRAC dans Aspen Plus. Comme dcrit dans lintroduction de cet
exercice, le calcul est effectu selon le principe dtages thoriques de sparation, i.e., on se
base sur une colonne plateaux pour laquelle on suppose que sur chacun des plateaux, ltat
dquilibre sera totalement atteint. Les flux quittant chaque tage de sparation sont donc en
quilibre thermique et chimique entre eux. Un systme dquations non linaires peut tre
obtenu partir des quations de bilan et dquilibre sur chaque tage. Dans une approche
squentielle-modulaire, ce systme va tre rsolu comme un seul module, cest--dire
simultanment. Pour cela, une mthode de calcul telle que la mthode de Newton-Raphson
par exemple sera employe.
La gomtrie des colonnes (le nombre de plateaux et les tages de Feed) doit tre spcifie
avant la simulation. Ce type de problme nest en gnral pas trivial et sera souvent rsolu par
une approche itrative dans laquelle diffrentes valeurs seront testes, souvent de manire
assez fastidieuse. Les tages sont numrots de haut en bas, le condenseur et le rebouilleur
correspondant respectivement au premier et au dernier tage de la colonne.
Une colonne de type RADFRAC comprend au moins trois flux de matire : un flux
dalimentation, le distillat (vapeur ou liquide condens) et le liquide de fond. En plus de ces
trois flux, un certain nombre de flux additionnels peuvent tre ajout la colonne, tels que :
- des flux dalimentation supplmentaires ;
- des soutirages de matire sortants dun tage (les soutirages liquide ou vapeur
permettent de dvier une partie de la matire vers lextrieur ou vers une autre
colonne);
- des flux de chaleur, entrant ou sortant.
Repartez de la simulation prcdente pour laquelle vous effacez le flowsheet (sauf le flux
FEED). Dessinez une colonne RADFRAC et connectez-la au flux dalimentation FEED.
Ajoutez-lui un distillat liquide condens DISTIL et un flux de fond BOTTOM ainsi que
reprsent la figure 18.
RADFRAC
DISTIL
FEED
BOTTOM
Figure 18 : colonne RADFRAC
Septembre 2012
38
Conservez les valeurs utilises pour le flux dalimentation FEED du train de sparation. La
colonne fonctionne en mode de calcul quilibre. Pour le nombre dtages, nous pouvons nous
rfrer au nombre dtage calcul prcdemment pour la premire colonne du train de
sparation, savoir 27 (ce nombre dtage comprend le condenseur et le rebouilleur).
Nous choisirons un condenseur total (toute la vapeur est condense, une partie du liquide est
reflu, lautre constitue le distillat) et un rebouilleur de type Kettle. Pour les spcifications
opratoires, nous pouvons galement nous rfrer aux valeurs obtenues prcdemment de
sorte que nous fixerons le dbit de distillat 45 kmol/h et le taux de reflux molaire 2,25.
Dans longlet Streams de la colonne, nous devons spcifier que lalimentation se fait au
sommet du 14me plateau. Le distillat et le flux de fond sortent tous deux sous forme liquide,
le condenseur tant un condenseur total. Dans longlet Pressure, menu droulant View, nous
spcifions le profil de pression de la colonne (il est galement possible de spcifier les pertes
de charge au condenseur et dans la colonne, ou bien les pertes de charge par section de
garnissage). La pression en tte de colonne (tage 1) est fixe 7,95bar et la pression en pied
de colonne (tage 27) 8,05bar.
Il nous est encore possible de donner des spcifications au condenseur : nous pouvons
ventuellement choisir la temprature de condensation (subcooled temperature) mais dans
notre cas nous fixons le nombre de degrs de sous-refroidissement 5C, cest--dire que la
temprature de condensation sera infrieure la temprature de bulle du distillat de 5C.
Enfin, nous fixons les efficacits des tages de la colonne en nous rendant dans le menu
efficiencies de la colonne. Nous choisirons les efficacits de Murphree, pour lesquelles nous
spcifierons les efficacits par tage. Dans longlet Vapor-Liquid, vous pouvez fixer
lefficacit 0,95 sur le premier et le dernier tage. Les efficacits des autres tages seront
interpoles par Aspen. Lancez le calcul.
Application
7. Une fois le calcul termin, observez les rsultats obtenus pour la colonne dans les
menus Results et Profiles du block RADFRAC. Reprsentez graphiquement les
profils suivants et discutez brivement de chacun de ces graphes.
a. le profil de temprature ;
b. le profil de fraction molaire pour le propane et le n-pentane dans la phase
vapeur ;
c. le profil de fraction molaire pour lisopentane et lisobutane dans la phase
vapeur.
8. Il est galement possible de simuler cette colonne en utilisant le mode cintique ratebased pour lequel la gomtrie et le garnissage de la colonne ainsi que les limitations
de transfert de matire et de chaleur sont prises en compte (ce qui nest pas le cas pour
le mode dquilibre o lon suppose ltablissement dun quilibre parfait sur chaque
plateau thorique). Cette modlisation est plus complexe, mais aussi plus prcise.
Testez le mode rate-based (Set up/Configuration/calculation type) en considrant un
nombre dtages inchang. Dans longlet Pack rating, crez une section de garnissage
structur de type Mellapak 250Y entre le 2me et le 26me tage. La colonne fait 1m de
section droite et est remplie avec 10m de garnissage. Cochez la case permettant le
calcul rate-based dans longlet Pack-rating/rate-based et acceptez les valeurs par
dfaut. Comment voluent la consommation nergtique de la colonne et la qualit de
la sparation des composs par rapport au modle RADFRAC ? Discutez.
Septembre 2012
39
Problme 3 : tude de la faisabilit dune colonne de sparation C4-C5

Dans cette tude, lobjectif est de vrifier les performances dune colonne. La colonne de
sparation principale C4-C5 dune unit doit tre arrte pour son entretien. On envisage de
dtourner temporairement son alimentation vers une petite colonne de rserve que vous
pouvez modliser par un bloc RADFRAC. Cette colonne est quipe de 14 plateaux : 1
condenseur total, 12 plateaux calotte et un rebouilleur Kettle. Lalimentation se situe au
sommet du plateau 8. Lexprience indique que lefficacit de Murphree des plateaux est de
90%.
Comme prcdemment, nous utilisons lquation de Redlich-Kwong-Soave pour estimer les
proprits thermodynamiques. Le dbit dalimentation est de 100 kmol/h. Le mlange est
disponible sa temprature de bulle, une pression de 8,2 bar. La composition molaire est la
mme que pour les deux premires applications. Dans ces conditions, la pression au
condenseur est de 8,1 bar, et la perte de charge de la colonne est de 0,15bar. Le condenseur
peut traiter un dbit de 175 kmol/h.
Comme notre modle comporte deux apports dnergie thermique (condenseur et rebouilleur)
et un soutirage (le distillat), nous devrons choisir trois spcifications. En se basant sur les
donnes, nous pouvons proposer :
- Le condenseur est total et ne sous-refroidit pas le condensat (il est donc obtenu sa
temprature de bulle) ;
- La colonne principale (qui doit tre arrte) travaille avec un rapport
distillat/alimentation D/F=0.489. Nous travaillerons avec le mme rapport.
- Le dbit de vapeur en tte de colonne entrant dans le condenseur est de 175 kmol/h.
Ce paramtre ne peut pas tre fix directement dans la colonne, mais nous pouvons
spcifier le dbit de reflux qui est la diffrence entre le dbit circulant dans le
condenseur et le dbit du distillat (48.9kmol/h), soit 126,1 kmol/h ;
Application
9. Dans la colonne qui doit tre arrte, le distillat contient moins de 7mol-%
disopentane et le rsidu moins de 3mol-% de n-butane. Peut-on tenir ces
spcifications avec la petite colonne en travaillant au mme rapport D/F? Si lon
veut amliorer la sparation et avoir moins de 3mol-% disopentane dans le
distillat, quel est le dbit maximal de distillat que lon pourra obtenir ?
Remarque : Vous pouvez crer une design spec sur le ratio D/F pour rsoudre ce problme.
10. Comment voluent les profils de temprature et de concentration si linjection du
mlange ne se fait plus sur le 8me plateau, mais sur le 5me ou le 11me ? Quelle
position vous semble la plus intressante ? Pourquoi ?
Septembre 2012
40
Exercice 9 : Introduction au logiciel HYSYS

Plusieurs logiciels de simulation commerciaux sont disponibles sur le march : Aspen Plus,
Hysys (rachet par Aspen), PRO II, Prosim, Belsim, etc. Ils ont de nombreux points
communs, et quelques particularits. Comme le logiciel Hysys est galement disponible sur
notre systme informatique, une introduction ses possibilits vous est prsente ci-aprs.
Hysys est principalement employ pour la modlisation de procds ptrochimiques. Comme
pour Aspen Plus, l'interface de dfinition de problme est graphique (ce nest pas toujours le
cas : certains logiciels se rsument des blocs dquations). Par contre, au contraire dAspen
Plus, les calculs d'appareils se droulent de manire interactive : ds qu'un appareil est
spcifi et calculable, il est calcul. Notons quil est possible de raliser des simulations du
comportement dynamique des procds, ce qui nest possible dans Aspen Plus quen
convertissant la simulation statique en mode dynamique et en utilisant ensuite le logiciel
Aspen Dynamics.
Le problme que nous allons traiter est identique au problme prsent dans lexercice
prcdent : il s'agit de la simulation d'une colonne de remplacement pour la sparation C4-C5.
Dmarrage dune nouvelle simulation
La premire chose faire quand on entre dans Hysys est de choisir son ensemble d'units.
Pour cela, il vous faut slectionner Preferences dans le menu Tools. Choisissez le set d'units
dsir dans longlet Variables. Hysys ne vous permettra pas de changer aucun des trois
ensembles dunits prsents par dfaut. Dans notre cas, nous allons choisir le set EuroSI.
Vous pouvez fermer la fentre des prfrences, celles-ci ont t enregistres sans que vous
nayez le confirmer.
Remarque : il est possible de modifier un set de base dunits. Pour cela, il suffit dappuyer
sur le bouton Clone Unit Set aprs avoir positionn le curseur sur le set de base voulu. Un
nouvel ensemble dunits va apparatre automatiquement avec comme nom par dfaut
NewUser que lon peut changer. Dans ce nouvel ensemble dunits, il vous est possible de
changer dunit pour chaque variable au cas par cas.
Pour crer une nouvelle simulation, il suffit de slectionner New Case dans le menu File. Les
tapes ncessaires la ralisation dun modle restent similaire celles dcouvertes dans
Aspen Plus :
- entrer les composants et choisir les mthodes de calcul des proprits
thermodynamiques (dans la fentre Simulation Basis Manager)
- raliser le shma du procd, paramtrer les flux dentre et les blocs du systme (dans
la fentre Simulation Environment).
Ltape suivante est donc l'introduction des composants du mlange. Slectionner la page
Components en cliquant sur son signet dans la fentre Simulation Basis Manager. Pour
slectionner les composants dsirs pour la simulation, vous pouvez faire une recherche par
nom et double-cliquer sur le composant dans la liste des composants disponibles dans la base
de donnes. Slectionnez ainsi les composants suivants dans la base de donnes Hysys :
propane (C3), i-butane (i-C4), n-butane (n-C4), i-pentane (i-C5), n-pentane (n-C5). Il est
Septembre 2012
41
possible de visualiser les proprits des composants par la touche View Component. Fermez
la fentre de choix des composants.
Vous pouvez maintenant crer un Fluid Package contenant les composants et les mthodes
pour le calcul des proprits des fluides (quations d'tat). Pour ce faire, dans longlet Fluid
Package, pressez la touche Add et choisissez l'quation d'tat voulue. Dans notre cas, ce sera
Peng - Robinson. Vous pouvez constater la prsence du nouveau Fluid package que lon
vient de crer, Basis-1 (nom par dfaut) avec le nombre de composants (NC) et lquation
dtat utilise. Maintenant que la base est dfinie, il est possible dinstaller les flux et les
oprations dans lenvironnement de simulation.
Environnement de simulation
Pour entrer dans cet environnement, pressez le bouton Enter Simulation Environment. Vous
pourrez revenir la dfinition du Fluid Package tout moment en cliquant sur Enter Basis
Environment dans longlet Simulation. Dans lenvironnement de simulation, il y a trois vues
initiales: PFD (Process Flow Diagrams, cest le flowsheet); Workbook (accessible depuis
Tools/Workbook) et Summary (accessible galement via longlet Tools).
Nous allons crer un flux par l'intermdiaire du Workbook. Pour ce faire, une fois le
Workbook ouvert, encodez alim , le nom du nouveau flux, dans la case New de la page
Material Streams du Workbook et pressez sur <Enter>. Hysys va crer automatiquement le
nouveau flux. Il vous reste ensuite remplir les cases avec les proprits que vous connaissez
(voir Problme 3 de lexercice 8).
Vapour fraction:
0
Pressure (bar):
8.2
Molar flow (kmol/h):
100
Pour entrer les compositions, double-cliquez sur la cellule molar flow ou bien slectionnez le
signet Compositions. Les compositions en fractions molaires sont les suivantes :
C3 : 5 %
iC4 : 15 % nC4 : 25 % iC5 : 20 % nC5 : 35 %
Le flux est maintenant compltement dfini, ce qui permet Hysys de calculer les proprits
restantes. Nous allons ajouter la colonne de distillation. Pour cela, entrez dans la page Unit
Ops du Workbook et cliquer sur Add UnitOp. Il faut ensuite ajouter lunit dsire : dans
notre cas cette unit est Distillation, une colonne de distillation.
La premire page de l'Input Expert apparat pour nous aider spcifier la colonne. Lorsque
vous installez une colonne, Hysys fournit certaines valeurs par dfaut qu'il est possible de
modifier. C'est le cas pour le nombre d'tages (valeur par dfaut = 10). Dans notre cas, nous
souhaitons une colonne 12 tages. Dans la cellule Inlet Streams, indiquer le flux
d'alimentation que vous avez dfini pralablement (alim) : Hysys va par dfaut placer
l'alimentation au milieu de la colonne. Vous pouvez la placer sur le 6me tage (indiqu par
6_Main TS).
Prcisez encore que le condenseur est total. Il ne reste alors quun flux de distillat liquide en
tte de colonne. Il vous reste nommer la colonne et les flux :
Septembre 2012
42
Condenser Energy Stream:

Ovhd Liquid Outlet:
Reboiler Energy Stream:
Bottom Liquid Outlet:
QC
DIST
QB
RES
Pressez le bouton Next pour avancer la page suivante o vous conservez la configuration par
dfaut du rebouilleur. Dfinissez ensuite sur la page suivante le profil de pression. Dans notre
cas, le profil de pression est suppos linaire. Entrez les valeurs suivantes:
Condenser Pressure : 7.95 bar
Condenser Pressure Drop : 0.05bar
Reboiler Pressure Drop : 0.05bar
Reboiler Pressure :
8.05 bar
Presser le bouton Next pour avancer la page Optional Estimate. Bien que Hysys ne require
pas d'estimation pour faire converger la colonne, une bonne estimation permettra d'atteindre la
solution plus rapidement. De la mme manire, une estimation des tempratures dans le
modle RADFRAC dAspen Plus peut galement aider la convergence, mme si en principe
vous nen avez pas eu besoin lors du problme prcdent. La temprature estime au
condenseur est de 60C, celle au bouilleur de 90C. Enfin, la dernire page de lInput
Expert, vous pouvez spcifier que le flux liquide de distillat vaut 45 kmol/h, et que le taux de
reflux est gal 2,25 en premire approximation. Pressez ensuite le bouton Done pour
terminer la spcification des connections de la colonne.
Dans longlet Monitor sur la gauche de la fentre de la colonne, dcochez loption Active du
taux de reflux afin que le taux de reflux spcifi soit uniquement une approximation et pas
une spcification. Le nombre de degrs de libert augmente de 1. Il faut donc rajouter une
spcification. Pour cela, ajoutez une spcification dans longlet Specs sur votre gauche. On
prcisera que le flux de vapeur Column Vapour Flow passant du premier tage au condenseur
(spcifiez pour ltage 1_Main TS) soit de 175 kmol/h. La nouvelle spcification apparat
dans la liste. Vous pouvez ensuite retourner la page Monitor et lancer le calcul en appuyant
sur le bouton Run. Le taux de reflux calcul vaut 2,89.
Rsolution
La colonne converge trs rapidement comme lindique le Control Panel dans le bas de votre
fentre. Nhsitez pas lagrandir pour y lire dventuels messages derreur. Le profil de
temprature est fournit dans le coin suprieur droit de la page de la colonne, onglet Monitor.
Vous pouvez avoir accs des informations plus dtailles en tudiant le Flowsheet
Summary, accessible depuis le menu Flowsheet. Vous pouvez galement visualiser la
colonne et les rsultats depuis la fentre PFD, accessible dans le menu Tools. Comme dans
Aspen Plus, double-cliquer sur un lment permet den voir les dtails. Enfin, les rsultats
sont galement accessibles via le Workbook.
Septembre 2012
43
Figure 19 : Flowsheet (PFD) dune colonne de distillation dans Hysys
Ajout dunits supplmentaires

Nous allons prsent utiliser la chaleur contenue dans le distillat pour prchauffer
lalimentation de la colonne dans un changeur de chaleur. Comme le distillat ne contient pas
assez de chaleur pour prchauffer suffisamment lalimentation de la colonne, nous rajouterons
un deuxime changeur ainsi quune pompe pour vaincre les pertes de charge lies ces
changeurs. Le diagramme rsultant est reprsent la figure 20 :
Figure 20 : Schma de la colonne de distillation avec rcupration de chaleur
A partir du PFD, slectionnez un changeur de chaleur (Heat Exchanger) dans la palette

(accessible laide de la touche <F4>) et positionnez-le dans votre diagramme. Doublecliquez sur lchangeur pour en spcifier les caractristiques. Vous le nommerez HTX1. Du
ct des tubes, le flux dentre est Alim1, et vous pouvez nommer le flux de sortie Alim2. Du
ct Shell, le flux dentre est le flux Dist, celui de sortie sappellera Dist2. Dans longlet
Parameters de lchangeur de chaleur, vous pouvez choisir un design dchangeur de type
Weighted (il sagit dun modle de calcul de lchangeur par sections lquilibre), et
spcifier une perte de charge de 0,1 bar de chaque ct de lchangeur (Tube et Shell). Le
coefficient de transfert global de chaleur UA vaut 10.000kJ/C-h.
Il convient maintenant de dfinir le nouveau flux dalimentation Alim1 (laissez le flux initial
Alim connect tant que la boucle ne sera pas boucle). Le flux Alim1 se trouve la
Septembre 2012
44
temprature de 20C et une pression de 8,2bar. Son dbit et sa composition sont les mmes
que prcdemment. Pour spcifier cela, double-cliquez sur le flux Alim1 dans le PFD.
Vous constaterez que la temprature du flux Alim2 atteinte la sortie de lchangeur (32C)
est insuffisante pour quil soit inject directement dans la colonne en remplacement du flux
Alim. Nous allons donc rajouter un deuxime changeur de chaleur, ainsi quune pompe entre
les deux changeurs pour compenser les pertes de charge de ces changeurs, sans quoi la
pression du flux serait infrieure celle de la colonne. Effectuez ces ajouts sur le PFD de la
mme manire que prcdemment. Attention, pour lchangeur, il sagit cette fois-ci dun
bloc de type Heater.
Vous nommerez la pompe Pump et spcifierez quelle est alimente par le flux Alim2 et le
flux de sortie se nomme Alim3. Il faut galement lui attacher un flux de travail W. Dans la
rubrique Parameters, la pompe provoque une augmentation de pression de 0,2bar et son
efficacit adiabatique est fixe 75%.
Les caractristiques du second changeur, HTX2, sont les suivantes : le flux entrant est
Alim3, le flux sortant Alim4. Le flux dnergie se nomme Q. La perte de charge vaut 0.1bar.
On demande un delta T de 58,14C dans lchangeur car nous souhaitons que la temprature
du flux Alim4 soit de 90C.
En prchauffant lalimentation par le distillat, nous avons introduit un recyclage. Il est donc
ncessaire dintroduire une unit de promotion de convergence qui, tout comme dans Aspen
Plus, va aider le problme converger en utilisant le flux de coupe que vous allez lui
spcifier. Dans Tools/Workbooks/Case, cliquez sur Add Unit Op. De la mme faon que nous
avons ajout une colonne, nous ajoutons un bloc de type Logical, Recycle (on peut aussi le
rajouter via la palette, cest un carr sur pointe avec un R vert lintrieur).
Contrairement Aspen, lunit de promotion de convergence ncessite deux flux sur le
schma. Par un processus itratif, elle va tenter de faire diminuer les carts entre ces deux
flux. Dans notre cas, ces flux sont Alim4 en entre et Alim en sortie.
Vous pouvez maintenant lancer le calcul (relancez le calcul de la colonne si ncessaire) et
examiner les rsultats.
Application
11. Pourquoi na-t-on pas rajout la pompe aprs le deuxime changeur, mais bien
entre les deux changeurs de chaleur? Que se passe-t-il si on le fait ?
12. Serait-ce une bonne ide de rcuprer la chaleur du flux Res dans lchangeur
HTX1 au lieu de celle du flux Dist ? Quy gagne-t-on concrtement ?
13. En revenant la configuration o la chaleur du flux Dist est rcupre dans HTX1,
et en examinant votre rponse la question prcdente concernant le flux Res,
serait-il possible de rcuprer de lnergie des deux flux ? Sans pour autant faire la
simulation, quelle configuration proposeriez-vous pour cela ?
Septembre 2012
45

TPAspen Complet PDF

Transféré par

Droits d'auteur :

Formats disponibles

TPAspen Complet PDF

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

TPAspen Complet PDF

Transféré par

Droits d'auteur :

Formats disponibles

Modlisation des

Grands Systmes Chimiques

Notes de Travaux Pratiques

Laboratoire dAnalyse et Synthse des Systmes Chimiques

Travaux Pratiques SYST0004

Etude dun procd avec le logiciel Aspen Plus

Travaux Pratiques SYST0004

Ensuite, la simulation de colonnes de distillation illustrera les techniques de rsolution globale

Droulement des travaux pratiques

1re Partie : Introduction au logiciel Aspen Plus et simulation de squences simples

2me Partie : Boucles dans un procd et spcifications de design

3me Partie : Optimisation de procds industriels, simulation de colonnes distiller

4me Partie : Introduction la validation de donnes et au logiciel Vali de Belsim.

Travaux Pratiques SYST0004

Exercice 1.1 : Echangeur de chaleur avec flash1

Figure 1 : Schma (flowsheet) de lexercice 1

Lancez linterface du logiciel Aspen Plus V7.3.2 :

Le premier travail effectuer est celui de la conception du Flowsheet ou Process Flow

Travaux Pratiques SYST0004

lexplorateur sur votre gauche, encodez les composants du procd.

Ouvrez maintenant le panneau de contrle (touche F7 ou bouton Control Panel dans la

Travaux Pratiques SYST0004

Au vu de ce premier exercice, quelles sont les tapes principales de la

3. Un ballon flash possde deux degrs de libert. Si au lieu dimposer un flash

Exercice 1.2 : Analyse des proprits des composs

Dans la partie Properties de linterface, encodez leau comme composant, et choisissez

Travaux Pratiques SYST0004

5. Construisez le diagramme T-xy du mlange Ethanol/Eau pression

Exercice 2 : Raction de synthse dammoniac

Figure 2 : Racteur de synthse dammoniac

Travaux Pratiques SYST0004

Connectez trois flux de matire au racteur, un lentre et deux en sortie

De retour dans la partie Properties, dans la rubrique Components du navigateur, entrez

o Dans la rubrique Methods du navigateur, choisissez la mthode RK-Soave (elle

Travaux Pratiques SYST0004

Dans la rubrique Blocks, le racteur fonctionne de faon adiabatique et la pression de

Travaux Pratiques SYST0004

Travaux Pratiques SYST0004

Exercice 3 : Calcul dune squence dunits

Figure 3 : schma du procd amlior de synthse dammoniac

Dans linterface Simulation, affichez le schma du procd (Onglet View/ Flowsheet).

Travaux Pratiques SYST0004

Vous pouvez accder la spcification du ballon FLASH en double-cliquant dessus. La

Figure 4 : Schma de la synthse dammoniac avec train de compression

Travaux Pratiques SYST0004

De retour sur la fentre du schma, rajoutez un changeur de chaleur HTX3 (modle

Figure 5 : ajout dun changeur de chaleur avant le sparateur flash

Figure 6 : Schma du rsultat final de lexercice 3

Les caractristiques du compresseur

Travaux Pratiques SYST0004

Travaux Pratiques SYST0004

Exercice 4 : Recyclage de matire

un compresseur COMP3 isentropique, avec une efficacit isentropique de 0,75 et

Travaux Pratiques SYST0004

un changeur de chaleur (Heater) que nous nommerons HTX3 aprs avoir

Figure 9 : Boucle de recyclage dans Aspen Plus

Travaux Pratiques SYST0004

Travaux Pratiques SYST0004

Figure 10 : Schma du procd avec boucle de recyclage