MIZZI Benoit
MIZZI Benoit
Ecole doctorale :
Mécanique, Energétique, Génie civil, Procédés (MEGeP)
Unité de recherche :
Laboratoire de Génie Chimique (L.G.C.)
Directeur(s) de Thèse :
M. MICHEL MEYER
M. LAURENT PRAT
Rapporteurs :
M. LUC MARCHAL, UNIVERSITE DE NANTES
Mme VIOLAINE ATHES, AGROPARISTECH
Membre(s) du jury :
1 M. CHRISTOPHE GOURDON, INP TOULOUSE, Président
2 M. LAURENT PRAT, INP TOULOUSE, Membre
2 M. MICHEL MEYER, INP TOULOUSE, Membre
A toute ma famille de sang ou de cœur que je vois plus ou moins souvent.
A mes parents, pour tout ce qu’ils donnent sans rien attendre en retour, MERCI.
J’ai des idées simples, je suis un homme simple ; donc tu comprends ce que je raconte.
(Joey Starr, Pensées répliques & autres tac-tac, 2016)
i
Remerciement
Pour commencer ces remerciements, je tiens à citer mes deux directeurs de thèse
Michel Meyer et Laurent Prat. Merci à tous les deux pour votre encadrement durant ces
trois ans pour m’avoir fait grandir scientifiquement même si nos réunions communes
me prouvent qu’il reste encore du chemin à parcourir. Merci Michel de t’être soucié
aussi bien de ma santé mentale et physique, merci Laurent pour ta bonne humeur
perpétuelle et tes conseils avisés dans les moments les plus difficiles.
Je tiens à remercier les deux directeurs du Laboratoire de Génie Chimique que j’ai
connus, Béatrice Biscans et Pierre Aimar. Au travers vous je tiens aussi à remercier les
deux directions du LGC, merci de faire du LGC un laboratoire performant mais surtout
un endroit convivial ou on se sent bien. Merci plus particulièrement à Pierre Aimar qui
lors d’une discussion alors que je n’étais encore que votre élève, m’a sûrement
convaincu de faire une thèse.
iii
d’échanger avec vous tous que ce soit dans la sphère professionnelle ou extra
professionnelle. Merci pour votre accueil, votre aide, vos conseils, vos soutiens, vos
coups de mains, votre écoute, merci pour nos discussions nos cafés nos soirées et autres.
Parmi tous les membres du personnel, je tiens plus particulièrement à remercier Alain
pour ces moments de « récréations ». Merci plus particulièrement aussi à Frank pour sa
disponibilité quotidienne sa bonne humeur, nos discussions rugbystiques nos cafés
échangés même si je pense que je t’en dois encore pas mal. Merci à vous et à tous les
autres. Je tiens aussi à remercier le personnel de l’IFPEN à Lyon, le second site auquel je
suis rattaché. Merci pour votre accueil et votre écoute en toutes circonstances et je tiens
plus particulièrement à remercier Martine Claire et Camille de s’être si bien occupé de
moi mais aussi Pierre pour sa collaboration sur un sujet et des expérimentations
épineuses.
Je tiens à remercier aussi les personnes qui m’ont aidé pour la rédaction de ce
manuscrit surtout d’un point de vue orthographe. Merci à Leslie et Martine, et je tiens à
m’excuser encore une fois auprès de Damien et Michel de ne pas avoir ouvert la boite à
« s » mais promis je vais m’y mettre.
Travailler sur deux sites m’a permis de rencontrer deux groupes de thésards que
je tiens à remercier. Merci à tous les thésards de L’IFPEN pour leur accueil lors de mes
séjours sur Lyon ou Paris, merci pour ces bons moments partagés ces soirées parfois
inoubliables et pour ce voyage formidable au Canada. Adeline, merci pour ce stage
effectué à l’IFPEN sur mon sujet de thèse, merci pour cette collaboration et pour le
travail réalisé. Je tiens à remercier maintenant les thésards du LGC, les anciens, Max
Mazu Cécile Séda Jennifer Martial Mylène Sara Claire, les nouveaux, Claire Omar Fatine
Manu Guillaume Mathieu Lucille, et ceux de ma génération. Merci aussi pour ces bons
moments, ces parties de tarot un peu folles parfois, ces soirées aux thèmes multiples et
variés. Merci aussi pour votre soutien votre aide, merci de former cette petite famille
fort sympathique du LGC. Je terminerai les remerciements des thésards en citant ceux de
ma génération, Alexandre Léo Doriane Kévin Ségolène Pierre Jérémy, en vous félicitant
de votre travail accompli pendant ces trois mêmes années que moi et merci pour
certains d’entre vous d’avoir fait partie de l’association des doctorants du LGC,
l’ALAMBIC, pour cette aventure parfois éprouvante mais ô combien réjouissante.
iv
Le manque de place m’a amené à aménager dans un bureau du département PSI
du LGC, merci tout d’abord de m’avoir accueilli chez vous, même si toutefois je pense
que je me suis un peu imposé. Merci à Guillaume B. Guillaume W. Lucille Alexandre et
autres stagiaires et doctorants voisins qui sont passés par ce bureau 2-r3-1 pour créer
une ambiance conviviale et parfois un peu trop agité.
Merci à mes amis du rugby et autres, de soirée d’un jour d’une semaine… Durant ces
trois ans on croise et on sympathise avec beaucoup de personnes et merci à eux de
m’avoir permis de penser, de m’investir, dans autre chose que ma thèse et d’avoir fait
semblant de s’intéresser à mon travail. Merci à tous ceux avec qui j’ai joué à Muret et à la
Vallée du Girou, merci aux entraineurs et notamment Thierry et Didier d’avoir compris
ma situation dans les moments difficiles ou dans la dernière ligne droite.
Et je ne pourrais pas terminer de remercier mes amis sans vous citer vous,
William Paul Etienne et Florian pour votre soutien quasi quotidien et vos
divertissements tout aussi plaisant, les virées en bateau, celles à Biarritz et Barcelone,
les soirées au Club à l’Ubu et tous les autres établissements de cette belle ville de
Toulouse qui est chère à notre cœur. Encore une fois merci pour tous ces moments
partagés ensemble et souhaitons qu’il y en ait encore beaucoup. Merci aussi et surtout à
Pierre Benji Yann Louis-Jo Tutu pour tous les bon moments passés, les voyages
organisés dans différentes ville du sud-ouest pour des évènements qu’on appelle Féria,
les virées au ski ou à Madrid, les sorties explosives sur Marmande au dolce et au TNT.
Malgré l’éloignement et ces rares moments partagés ensemble, me permettent de me
rappeler où sont mes racines et le bonheur qu’ils m’apportent lors de nos retrouvailles
pour ces différents évènements intenses et festifs.
Une thèse c’est des partenaires des encadrants des amis, ses lots de bons et de
moins bon moments, ses soutiens ses aides ses conseils. Mais pour mener à bien un tel
projet il faut quelque chose d’essentiel et de primordial. J’ai jusqu’à maintenant parlé de
ma famille de cœur mais la plus importante est la famille de sang. Je terminerai donc par
remercier toutes ma famille, mes grands-parents, mes oncles mes tantes mes cousins et
surtout mes parents et ma sœur. Merci à vous de constituer mes racines et de me
permettre de grandir et de m’élever dans tous les projets que j’entreprends. Merci à ma
petite sœur Béatrice pour sa joie de vivre et ses discussions toujours positives, merci à
mes parents pour tout et d’avoir fait en sorte que je sois l’homme que je suis, je ne
v
pourrais jamais vous rendre tous ce que vous m’avez donné mais la fierté dans les yeux
de ma mère et sûrement dans celle de mon père trop pudique pour la montrer m’ont
réconforté dans l’idée que ce que je fais est bien à vos yeux et c’est ce qui m’importe le
plus.
vi
Table des matiè res
RÉSUMÉ/ABSTRACT .................................................................................................................. xi
INTRODUCTION GÉNÉRALE ..................................................................................................... 1
CHAPITRE I ÉTAT DE L’ART ...................................................................................................... 7
I.1 Cadre scientifique de l’étude.................................................................................... 7
I.1.1 Intensifications des procédés .............................................................................. 7
I.1.1.1 Les équipements pour l’intensification des procédés ....................... 8
I.1.1.2 Les méthodes pour l’intensification des procédés ............................. 9
I.1.2 Production par fermentation ............................................................................ 11
I.1.3 Les procédés de séparations réactives ......................................................... 13
I.2 Les méthodes de conception en génie des procédés ................................... 15
I.2.1 Conception globale d’un procédé chimique ................................................ 15
I.2.2 Méthode de conception d’une unité d’un procédé chimique. .............. 17
I.2.3 Méthode de conception d’un procédé de distillation réactive ............. 19
I.3 Production de l’acide succinique......................................................................... 21
I.3.1 État des lieux de la production et du marché actuel de l’acide
succinique .................................................................................................................................... 21
I.3.2 Production de l’acide succinique par voie fermentaire .......................... 23
I.3.2.1 Synthèse de l'acide succinique ................................................................ 24
I.3.2.2 Procédé de récupération et de purification de l’acide succinique
biosourcé .............................................................................................................................. 27
I.4 Extraction liquide-liquide ...................................................................................... 34
I.4.1 Définition et champs d’application ................................................................. 34
I.4.2 Outils thermodynamiques et dimensionnement ...................................... 36
I.4.3 Technologies ........................................................................................................... 42
I.4.4 Extraction liquide-liquide réactive ................................................................. 48
I.4.4.1 Méthodes conceptuelle de conception d’extraction liquide-
liquide réactive ............................................................................................................................ 49
I.4.4.2 Méthodes géométriques de conception d’extraction liquide-
liquide réactive (Olán-Acosta et al., 2014) ........................................................................ 53
I.4.4.3 Analyse de faisabilité................................................................................... 55
I.5 Conclusion état de l’art............................................................................................ 57
CHAPITRE II UNE METHODOLOGIE POUR LA CONCEPTION D’UN PROCEDE D’EXTRACTION
LIQUIDE-LIQUIDE REACTIF ................................................................................................................... 61
vii
II.2.1 Equilibre chimique d’une réaction équilibrée .......................................... 64
II.2.2 Equilibre physique d’un système polyphasique ...................................... 66
II.2.2.1 Caractérisation de l’état d’équilibre de phases ................................. 66
II.2.2.2 Application aux systèmes liquide-liquide ........................................... 67
II.2.3 Equilibre physico-chimique d’un système réactif polyphasique :
introduction au diagramme liquide-liquide réactif ..................................................... 68
II.2.3.1 Formulation mathématique des équilibres physico-chimiques . 68
II.2.3.2 Variance du système.................................................................................... 71
II.2.3.3 La notion de compositions réactives .................................................... 72
II.2.3.4 Le diagramme liquide-liquide réactif ................................................... 74
a. Exemple : réactions de type R1 + R2 ↔ P1 + P2 ............................ 75
b. Exemple de deux réactions successives............................................... 76
II.3 Description détaillée de la méthode de conception ..................................... 78
II.3.1 Analyse de faisabilité .......................................................................................... 79
II.3.1.1 Introduction des courbes de raffinats réactifs : les courbes de
résidus réactifs en distillation réactive .............................................................................. 80
II.3.1.2 Génération et analyse des courbes de raffinat réactif .................... 81
II.3.1.3 Détermination des taux de solvant limites et analyses des
régimes permanents .................................................................................................................. 85
II.3.2 Synthèse et dimensionnement : détermination du nombre d’étages
théoriques.................................................................................................................................... 91
II.3.2.1 Saisie des données de conception .......................................................... 92
II.3.2.2 Génération des profils de composition ................................................ 92
II.3.2.3 Génération des profils de composition : cas d’une colonne à
reflux .............................................................................................................................. 94
II.3.2.4 Le critère d’arrêt et l’analyse des résultats ........................................ 95
II.3.3 Conception et validation.................................................................................... 97
II.3.3.1 Limitation cinétique des équilibres chimiques entre phases ...... 97
II.3.3.2 Simulation d’une colonne d’extraction liquide-liquide réactive 99
II.4Résumé de la méthode ..........................................................................................101
CHAPITRE III EXTRACTION LIQUIDE–LIQUIDE REACTIVE APPLIQUEE A L’ACIDE
SUCCINIQUE BIOSOURCE. ...................................................................................................................105
III.1 Stratégie A : l’alcool est l’agent extractant.....................................................107
III.1.1 Analyse de faisabilité .......................................................................................109
III.1.2 Synthèse et conception ...................................................................................112
III.1.3 Simulation ............................................................................................................113
viii
III.1.3.1 Influence du nombre d’étages théoriques ........................................115
III.1.3.2 Influence du taux de solvant ..................................................................117
III.1.3.3 Influence de la pureté du solvant .........................................................118
III.1.3.4 Influence du volume de rétention ........................................................119
III.1.4 Conclusion............................................................................................................121
III.2 Stratégie B : le solvant est inerte par rapport à la réaction ....................122
III.2.1 Analyse de faisabilité .......................................................................................122
III.2.2 Synthèse et conception ...................................................................................127
III.2.3 Simulation ............................................................................................................130
III.2.3.1 Influence du nombre d’étages théoriques ........................................131
III.2.3.2 Influence du taux de solvant ..................................................................132
III.2.3.3 Influence de la composition du solvant .............................................133
III.2.3.4 Influence du volume de rétention ........................................................134
III.2.4 Conclusion............................................................................................................135
III.3 Comparaison de procédé de récupération de l’acide succinique
biosourcé .................................................................................................................................135
III.3.1 Distillation du mélange binaire eau-acide succinique ........................137
III.3.2 Extraction liquide-liquide simple avec régénération du solvant ...139
III.3.3 Extraction liquide-liquide réactive avec régénération du solvant .141
III.3.3.1 Procédé A, extraction liquide réactive avec régénération du
solvant, le solvant est le réactif de la réaction ...............................................................141
III.3.3.2 Procédé B, extraction liquide réactive avec régénération du
solvant, le solvant est inerte .................................................................................................142
III.4 Conclusion : comparaison des différents procédés de récupération de
l’acide succinique. ................................................................................................................145
CONCLUSION GÉNÉRALE...............................................................................................149
PERSPECTIVES .................................................................................................................152
ANNEXES ............................................................................................................................155
NOMENCLATURE .............................................................................................................176
LISTES DES FIGURES ......................................................................................................178
LISTE DES TABLEAUX ....................................................................................................184
RÉFÉRENCES .....................................................................................................................188
ix
RES SUMES /ABSTRACT
Le couplage fonctionnel des opérations de séparation et de réaction ainsi que les
bioprocédés sont deux axes de recherche largement explorés. Cependant, l’industrie du
génie des procédés a du mal à se tourner vers des technologies de ce type car il demeure
un réel manque de connaissances et d’outils de conception pour ce genre de procédés.
Une méthodologie de conception générale pour l'extraction liquide-liquide
réactive est introduite dans cette étude. Elle est composée de trois étapes différentes:
l'analyse de faisabilité, la synthèse ou dimensionnement du procédé et la validation par
simulation. Cette méthodologie conduit à des paramètres structuraux et opératoires de
la colonne étudiée à partir seulement des informations concernant le comportement
physico-chimique du système étudié, en exploitant les équations d’équilibre chimique et
entre phase ainsi que les bilans matières. Les résultats de cette méthode sont un bon
point de départ pour une étude d'optimisation ou d'un processus de calcul
d'investissement. Cette méthodologie a été appliquée à différentes études de cas:
regroupant deux stratégies différentes d'extraction avec plusieurs solvants pour
récupérer l'acide succinique dans un milieu de fermentation.
xi
INTRODUCTION GES NES RALE
La notion de développement durable définie dans le rapport Brundtland en 1987
prend une ampleur de plus en plus significative dans l’industrie chimique. En particulier,
la recherche dans le domaine de l’ingénierie des bioprocédés ou du génie biochimique
connait une forte expansion depuis une trentaine d’années. Les bioprocédés peuvent
s’inscrire dans une grande diversité de contextes industriels et s’inscrivent dans le
développement durable. En effet, l’avantage des bioprocédés est de remplacer les
matières premières fossiles, tout en cherchant à assurer une production en molécules
d’intérêt économiquement et techniquement viable.
L’une des voies pour répondre à ce frein est d’utiliser les concepts, méthodes et
outils développés en Génie des Procédés, notamment ceux de l’intensification des
procédés. Comparés aux procédés déjà existants, les procédés durables doivent
répondre aux mêmes exigences de productivité tout en étant plus sûrs, moins couteux et
plus respectueux de l’environnement. Les procédés de séparation réactive sont tout à
fait caractéristiques de ce domaine de recherche. Les procédés de séparation réactive
combinent dans un même appareil réaction et séparation et présentent différents
avantages économiques comme techniques.
1
fortement dilués comme le sont les milieux fermentaires. Pour l’extraction liquide-
liquide réactive, les données bibliographiques sont moins abondantes, cependant le
couplage fonctionnel réaction séparation est très intéressant et plus adapté dans ce cas,
la forme transformée de l’acide succinique pouvant faciliter la séparation liquide-liquide
mais aussi le produit de la réaction pouvant être directement valorisé.
Dans le premier chapitre le cadre scientifique général de cette étude est détaillé.
La notion d’intensification des procédés est introduite. Les axes de recherche qui
s’inscrivent dans cette notion sont brièvement présentés au travers de la production par
fermentation en général et des équipements de séparation couramment rencontrés. Les
méthodes de conception globale en génie des procédés sont introduites afin de
connaître les moyens de concevoir un procédé sur lequel notre méthode pourra
s’appuyer. Une partie du chapitre est destinée aux bioprocédés de production de l’acide
succinique existants.
Enfin dans le chapitre III, la méthode est appliquée à deux exemples différents
mettant en jeu deux stratégies d’extraction différentes qui présentent un intérêt
particulier de par le nombre de constituants et de par la nature de certains constituants.
Les différents aspects de la méthode sont ainsi illustrés. Une partie du chapitre est
dédiée à la comparaison de différents procédés de récupération de l’acide succinique
2
pour montrer si l’extraction liquide-liquide réactive est une alternative viable et qui
répond aux enjeux de développement durable.
3
Chapitre I
État de l’art
Ce chapitre introductif se propose de définir le cadre dans lequel s’inscrivent ces
travaux de recherche.
Tout d’abord dans la section I.1 est présenté le cadre scientifique de l’étude
notamment la tendance qui émerge en génie des procédés, l’intensification des procédés.
Les concepts qui nous intéressent et qui sont développés dans ce paragraphe sont la
production par voie fermentaire et les procédés de séparations réactives.
Dans la section I.2, nous nous attacherons à présenter les différentes méthodes de
conception de procédés en génie des procédés. Il s’agit d’une étape complexe et essentielle
avant la réalisation physique du procédé. Cette partie permettra de définir et de présenter
les différentes étapes de ces méthodologies et de limiter le cadre de notre étude.
5
Chapitre I État de l’art
7
L’intensification des procédés peut se présenter sous forme de méthodes ou de
technologies mais se résument principalement à diminuer les volumes utiles en
augmentant la productivité et la sélectivité. Ces mesures d’action sont regroupées dans
divers domaines d’investigation regroupés sur la Figure I-2. Les paragraphes suivants
présentent une liste non exhaustive des nouvelles technologies et méthodes qui
permettent d’illustrer l’étendu du champ d’investigation relatif à l’intensification des
procédés.
8
transfert de matière garantissant une meilleure pureté du produit. Ils sont surtout
utilisés pour des réactions très exothermiques.
• Les garnissages catalytiques (Cybulski, 1997) comme les garnissages
KATAPAK sont introduits dans des unités comme la distillation réactive pour assurer
une bonne redistribution radiale et un échange de matières efficace.
• Les microréacteurs (Mills et al. 2007) sont des réacteurs chimiques de
tailles relativement faibles, le diamètre des canaux qui les composent et où a lieu la
réaction étant de l’ordre de 10 à 100µm. L’arrangement de ces microcanaux permet de
servir différentes fonctions comme l’échange de chaleur, l’échange de matière, la
séparation ou encore le mélange. L’avantage important de ce type de système est de
pouvoir intégrer toutes ces opérations dans une même micro unité avec une rétention
globale très faible et donc intrinsèquement plus sûre.
9
que d’autres comme l’osmose inverse et la distillation membranaire consomme
beaucoup d’énergie.
10
I.1.2 Production par fermentation
11
végétales en produits normalement synthétisés par l’industrie pétrolière. Les matières
premières végétales peuvent varier, des plantes entières (céréales, canne à sucre,
betteraves, …), des produits de l’industrie agroalimentaire, le bois ou encore des algues.
Cette matière première est source d’hydrates de carbone (plus spécialement les sucres)
que les cellules mises en œuvre dans le bioprocédé vont transformer dans l’étape de
synthèse en produit d’intérêt.
Dans les études à l’échelle laboratoire, du glucose pur est utilisé mais les
procédés industriels nécessitent une matière première moins coûteuse. La biomasse,
matière première des bioprocédés, peut être catégorisée en quatre principales
catégories : la cellulose, l’hémicellulose, la lignine et les huiles. La cellulose et
l’hémicellulose sont dégradées par voies enzymatiques en sources carbonées alors que
la lignine n’est pas dégradable, la lignine étant donc peu utilisée dans les procédés
industriels. Les huiles sont au stade d’étude comme matière première pour les procédés
de fermentation. La Figure I-3 regroupe les différentes catégories de biomasse avec
quelques exemples de matières premières.
Biomasse
Sucrose
Sorbitol
Dans l’étape de purification, ces produits sont récupérés. Les produits ainsi
obtenus sont destinés aux domaines de l’énergie, de la chimie, de l’agroalimentaire et
des matériaux. Les marchés visés sont aussi divers, les produits de masse à faible valeur
ajoutée (carburants solvants) et les molécules produites à faible volume et à haute
valeur ajoutée (pharmaceutique, cosmétiques, …)
12
Le Tableau I-2 fait le parallèle entre les voies de production biologique et celles à
partir d’énergie fossile. La différence remarquable concerne sur les conditions
opératoires de chaque voie. En effet toutes les réactions citées nécessitent de hautes
températures et pressions, et les voies biologiques faisant intervenir des micro-
organismes vivants, les conditions sont plus douces et donc moins coûteuses en énergie.
L’un des leviers pour intensifier les étapes de réactions consiste à travailler avec
des appareils multifonctionnels. Un réacteur est rendu multifonctionnel dans l’optique
d’améliorer ses performances (rendement, productivité) en couplant la réaction et la
séparation. Cette synergie permet une amélioration des performances du procédé
(Seader and Henley, 2011, Treybal, 1951).
13
Deux fonctionnements peuvent être distingués dans le couplage de réaction et de
séparation :
14
nouveaux procédés et les méthodes pour concevoir le procédé susceptible d’intégrer ces
unités.
Phase Gaz
Adsorption Absorption
Réactive Réactive
Distillation
Réactive
Phase Liquide
Phase Solide Cristallisation réactive Aqueuse
Extraction
Réactive Extraction
Adsorption
Réactive
Réactive
15
Critères Critères
techniques économiques
Conception de procédés
durables
Critères
Critères sociaux
environnementaux
Figure I-6 : Les différentes étapes dans la conception des procédés durables (Azzaro-Pantel, 2010)
16
le procédé et vise à la synthèse du procédé et de toutes les unités qui le composent. Cette
étape permet aussi d’analyser et de choisir entre les différentes alternatives possibles
avec la plupart du temps des critères économiques. L’étape suivante est la réalisation du
procédé, la réalisation du flow-sheet et l’étude du procédé avec des outils de simulation.
Enfin la dernière étape concerne la construction physique du procédé, l’opérabilité et la
maintenance de celui-ci.
Les deux étapes concernées par cette étude sont les étapes de conception et de
développement. Bien que les dépenses majeures soient réalisées dans l’étape suivante,
ces deux étapes conditionnent les choix et donc, ces dépenses. En résumé, lorsque les
étapes de conception et de développement sont bien réalisées, les décisions prises ont
généralement un impact positif sur le résultat économique global de l'ensemble du
projet. Malheureusement les méthodes de conception et de développement utiles à ces
étapes sont la plupart du temps, difficilement exploitables, trop nombreuses, diversifiées
et non synthétiques. Du fait de l’importance de ces étapes et de ce manque de méthodes
hiérarchisées et synthétiques, la phase de conception de procédés chimiques est
extrêmement importante dans le processus de développement du procédé et a un
intérêt particulier pour nous.
17
Zone de l’unité
Zone opératoire (technologie)
réaction
(thermodynamique)
Zone opératoire
Zone opératoire
séparation d’une séparation
(thermodynamique) réactive
•Expérimentations
Validation •Simulation modèle
MESH et NEQ
18
L’étape de synthèse permet de déterminer le design de procédé pour atteindre
les objectifs fixés à partir des informations disponibles grâce à l’étape précédente
d’analyse (Hamid et al., 2010).
La méthodologie développée par Théry (2002) est très intéressante du fait de son
caractère séquentiel et hiérarchisé.
19
qu’avec une quantité restreinte de données thermodynamiques relatives au système
étudié. Cette approche se réfère au trois premiers points d’une procédure en cinq
étapes.
20
Données thermodynamiques du système
L'acide succinique est reconnu comme l'une des molécules d'avenir majeures.
Présent dans la nature chez les animaux et les végétaux, elle est largement utilisée
comme intermédiaire chimique pour des produits de la vie quotidienne dans des
applications alimentaires, pharmaceutiques, cosmétiques ou plastiques.
21
Figure I-10: Principales molécules dérivées de l'acide succinique (Zeikus, 1980)
22
Les paragraphes suivants présentent les différents micro-organismes, identifiés
et utilisés pour synthétiser l'acide succinique. Des recherches sont réalisées sur les
conditions opératoires de la fermentation pour optimiser les productivités et les
rendements des micro-organismes. Les milieux qui découlent de cette fermentation sont
très complexes et comportent, en plus de l'acide succinique, de nombreux sous-produits
comme des monoacides, des alcools, de la biomasse, les micro-organismes et leurs
milieux de culture. L’un des verrous essentiel à l’optimisation de ces procédés
industriels par voies biochimiques concerne la récupération et la séparation du produit
final, l'acide succinique. Il est résumé aussi dans le paragraphe I.3.2.2 les nombreuses
techniques et opérations unitaires envisageables pour récupérer l’acide succinique des
milieux fermentaires. Citons notamment la précipitation, l’électrolyse ou l’extraction
liquide-liquide avec une amine comme agent de séparation.
23
I.3.2.1 Synthèse de l'acide succinique
La synthèse de l’acide succinique par voie fermentaire met en jeu des micro-
organismes. Cette partie vise à présenter les micro-organismes et leurs performances en
termes de fermentation, de rendement et de productivité. Quatre micro-organismes font
l’objet de recherches approfondies et présentent de bonnes aptitudes pour la synthèse
de l’acide succinique, il s’agit d’Actinobacillus Succinogenes, Mannheimia
Succiniciproducens, Anaerobiospirillum Succiniproducens et Escherichia coli (Song and
Lee, 2006). Les deux premiers micro-organismes ont été isolés à partir de panse de
bovin. Les deux suivants quant à eux ont été isolés à partir de système digestif et de
matière fécale de mammifères. Ces micro-organismes ont été retenus pour leurs hauts
rendements où leurs fortes productivités. Ils ne présentent pas les deux caractéristiques
simultanément mais des travaux sont en cours sur leurs génomes afin de remplacer ou
de modifier certains gènes et d’améliorer l’une des performances sans diminuer l’autre
(Lee et al, 2005, Beauprez et al, 2010).
a. Actinobacillus succinogenes
24
b. Mannheimia succiniciproducens
c. Anaerobiospirillum succiniproducens
d. Escherichia coli
Escherichia coli est un micro-organisme encore peu utilisé car il ne produit pas
majoritairement de l’acide succinique. Les expérimentations effectuées sur ce micro-
organisme, ont déterminé que celui-ci utilise six schémas réactionnels pour transformer
25
la source carbonée en différentes espèces chimiques comme l’éthanol, l’acide formique
l’acide acétique et lactique. Des modifications du génome ont été entreprises afin de
favoriser certains chemins réactionnels, comme celui produisant l’acide succinique.
Grace à ces travaux, la fermentation du glucose permet d’obtenir un meilleur rendement
de l’ordre de 1,0 mole d’acide par mole de glucose et une productivité dans la normale
de l’ordre de 1,0 g.L-1.h-1 (Vemuri et al, 2002, Lin et al, 2005, Lin et al., 2010). Ce nombre
élevé de chemins réactionnels conduit à des jus de fermentation très complexes avec
une grande mixité de produits, ce qui rend le traitement de ces jus et la récupération de
l’acide succinique très difficile. Malgré cela et bien que les autres micro-organismes
produisent en majorité de l’acide succinique, E. coli est majoritairement choisie dans les
procédés industriels car elle nécessite des nutriments relativement simples et donc
moins coûteux ce qui réduit à la fois le coût et la complexité du procédé.
26
Concentration Rendement Productivité
Micro-organisme Substrat
(g/l) (-) (g/l/h)
A succinogenes
glucose 33,9 0,86 0,88
(Souche naturelle)
130Z
glucose 66,4 0,67 0,79
(Souche naturelle)
E coli HL27659K
(Souche génétiquement glucose 58,3 0,55 0,98
modifiée)
Tableau I-4 : Récapitulatif production de l'acide succinique par les principaux micro-organismes (Song and Lee,
2006)
a. Précipitation
27
succinate sont ensuite récupérés par filtration. L’acide succinique est obtenu par un
traitement avec un acide fort et subit finalement une dernière étape de cristallisation.
28
Figure I-12 – Procédé de précipitation Berglund & al
Figure I-11 – Procédé de précipitation Datta & al (1992)
(2001)
b. Electrolyse
29
Figure I-13 présente le principe de l'électrolyseur. Le sel de succinate récupéré après le
passage dans l'électrolyseur, sera converti dans un "water-splitting" électrolyse, Figure
I-14). Une membrane bipolaire formera à partir de l'eau des ions hydrogène pour
reformer l'acide succinique et des ions hydroxyde pour éliminer le contre ion Sodium. La
Figure I-15 donne le schéma simplifié de ce système. Une dernière étape de filtration et
de purification est mise en place pour récupérer de l'acide succinique avec une pureté de
80% et un taux de récupération compris entre 75% et 80%.
Figure I-13 - Schéma d'une électrolyse Figure I-14 - Schéma d'une “water splitting” électrolyse
30
Figure I-15 – Procédé d’électrolyse Glassner & al (1995)
c. Couplage estérification-séparation
L’acide succinique peut, par une réaction d'estérification, réagir avec un alcool et
donner un ester. Cette transformation chimique modifie les propriétés physico-
chimiques de l'espèce. Une séparation par extraction liquide-liquide peut alors être mise
en place. Deux cas de figures sont envisageables : opérer l'étape de réaction et de
séparation en deux étapes distinctes, ou effectuer la réaction et la séparation
simultanément.
Dans le premier procédé, développé par S.S. Mims (1978), l'acide succinique
réagit avec un alcool court, dans une phase aqueuse ; cette phase aqueuse est ensuite
mise en contact avec un solvant organique afin d'extraire l'ester. Dans l'exemple datant
de 1978 pris pour ce paragraphe, dans un premier temps la réaction d'estérification se
fait en milieu homogène, l'alcool étant court il est soluble dans l'eau. Dans un second
temps la solution aqueuse de toutes ces espèces est mise en contact avec le solvant, du
Benzène. Cette méthode demande trois étapes de réaction séparation, ce qui fait huit
opérations unitaires comme on peut le voir sur la Figure I-16. Par ailleurs, le benzène est
aujourd'hui classé CMR et doit être remplacé.
31
Dans un second procédé de V.P. Kuceski (1995), l’utilisation d'alcools immiscibles
avec l'eau est préférée. La solution aqueuse d'acide carboxylique est mise en contact
avec la phase organique composée d'alcool à chaine longue, le procédé de Kuceski utilise
de l’isoctyl alcool (C8H17OH). L’acide est transféré de la phase aqueuse à la phase
organique et réagit dans la phase organique avec l’alcool pour donner l’ester. Après un
passage dans ce procédé, environ 30% de l’acide succinique est récupéré. Les
performances de séparation ne sont pas aussi bonnes que les autres procédés présentés
dans ce paragraphe. L’acide succinique non récupéré dans la phase aqueuse peut l’être
grâce à des recyclages ou à d’autres procédés de séparation comme la cristallisation. Des
expériences ont été réalisées avec d’autres alcools avec des chaines carbonées plus
longues. Il en découle que plus l’ester formé a une masse molaire importante plus la
séparation est complète (V.P. Kuceski (1995)). Ceci justifie aussi l’utilisation des alcools
lourds, d’une part car ils ne sont pas miscibles avec l’eau et d’autre part ils augmentent
la masse molaire de l’ester formé.
Figure I-16 - Procédé de récupération des acides carboxyliques S.S. Mims (1968)
Extraction réactive
Une autre méthode pour récupérer l'acide succinique est l'extraction liquide-
liquide réactive. Le principe de cette méthode consiste à faire réagir l'acide avec un
réactif présent dans le solvant. Les propriétés physico-chimiques de l'acide seront ainsi
modifiées et le transfert de matière entre les deux phases favorisé.
32
Les solvants couramment utilisés sont des amines, la Tri-octylamine (TOA),
tripropylamine (TPA), des diluants comme des alcools à chaine longue sont utilisés pour
améliorer la solvatation des complexes formés entre l'acide et l'amine dans le solvant
(Krzyżaniak et al, 2013).
Le Tableau I-5 regroupe les performances des unités de séparations citées ci-
dessus. Il résume les taux de récupération des différentes unités de séparation et donne
quelques remarques sur l’unité.
33
I.4 Extraction liquide-liquide
L’extraction liquide-liquide est une opération unitaire classique dans l’industrie
du génie des procédés. C’est un procédé qui permet de séparer un ou plusieurs
constituants d’un mélange en faisant intervenir leurs différences de distribution entre
deux liquides non miscibles. La force de séparation dépend du transfert de matière d’un
constituant dans une phase vers une autre phase.
Le procédé d’extraction liquide-liquide est définie comme une unité où est mis en
contact deux phases qui ne peuvent pas ou peuvent mais seulement dans une proportion
très limitée se mélanger. Chaque phase est multi constituant et on distingue aussi les
courants liquide qui entrent et sortent de l’extracteur (Figure I-17).
Alimentation F Extrait E
A+C B+C
Phase C Phase
raffinat extrait
Raffinat R Solvant S
A B
34
Pour décrire le principe de l’extraction liquide-liquide simplement, prenons le cas
le plus simple traité avec l’extraction liquide. Considérons un mélange ternaire (A, B, C),
respectivement le diluant, le solvant et le soluté. L’alimentation est le courant qui
contient le constituant à séparer (C), et le constituant majoritaire du courant le diluant
A. Le solvant est le courant alimenté dans l’extracteur, il est composé du constituants B
qui peut être pur ou contenir quelques traces de soluté et de diluant. Le solvant est
recyclé à partir d’une unité de régénération qui souvent ne permet pas d’avoir une
séparation tranchée ni de récupérer du solvant pur. Ces deux phases en contact, raffinat
et extrait, permettent le transfert du soluté entre les deux phases et de le récupérer ainsi
dans la phase extrait. Cette phase extrait est récupérée en sortie de l’extracteur dans un
courant où le constituant B est majoritaire et où on retrouve le soluté (C). Ce courant
peut être traité par une unité de régénération pour recycler le constituant B vers
l’extracteur et obtenir un courant de production riche en soluté. Le dernier courant est
le raffinat qui est l’alimentation débarrassée d’une grande partie du soluté qui la
compose.
Dans l’industrie l’extraction liquide-liquide est une alternative aux autres unités
de séparation comme la distillation, la cristallisation ou l’adsorption. Les motivations
pour employer cette technique plutôt qu’une autre sont souvent d’ordre économique,
technologique ou physico-chimique. L’extraction est le plus souvent utilisée lorsque les
constituants à traiter sont sensibles thermiquement (antibiotique) ou que les points
d’ébullition sont proches (séparations d’aromatiques). La liste suivante présente de
façon non exhaustive les différentes utilisations de l’extraction liquide-liquide dans
l’industrie :
35
• Séparation de composés thermosensibles notamment dans les
bioprocédés.
• Séparation de mélange azéotropique, exemple de l’extraction de l’acide
acétique ou formique avec du MTBE comme solvant.
• Extraction de composés organiques dans une solution de sel,
(caprolactame)
• Extraction de sels d’une solution polymère, (polyols)
• Extraction de métaux dans les minerais et les eaux usées, (extraction du
cuivre)
• Récupération des combustibles nucléaires (procédé Purex CEA, 2008)
L’extraction liquide-liquide joue un rôle très important dans divers secteurs
industriels comme le montre la liste ci-dessus, cependant son intégration dans le
paysage industriel est inégale. Par exemple en pétrochimie le procédé est optimisé
tandis que dans le nucléaire ou la métallurgie nous sommes encore à l’étape de
recherche du procédé ou de solvant plus efficace (Cusack, Roger W; Glatz, Donald J.,
1996).
36
Figure I-18: Différents types d’équilibre liquide-liquide avec leurs fréquences
37
Figure I-19: Influence de la température sur l’isotherme de solubilité d’un mélange ternaire
Figure I-20: Évolution de l’isotherme de solubilité en Figure I-21 Évolution de l’isotherme de solubilité en
fonction de la température pour un système sans fonction de la température pour un système avec une
température critique de solubilité température critique de solubilité du binaire A-S
38
Figure I-22: Évolution de la température qui transforme un diagramme de type II en diagramme de type I
Dans la littérature (Teixeira et al., 2009, Wibowo and Ng, 2002), on trouve des
méthodes permettant de représenter des systèmes quaternaire et quinaire. Cette
représentation est appliquée dans le domaine de la parfumerie est n’est pas adaptée à la
détermination d’un design de colonne. Un outil est présenté dans le paragraphe
74II.2.3.4, pour le dimensionnement d’une extraction liquide-liquide mettant en jeu un
système avec plus de trois constituants.
Maintenant que nous avons introduit l’outil indispensable pour tous calculs de
dimensionnement de colonne, nous allons présenter la méthode de dimensionnement de
la colonne d’extraction liquide-liquide pour traiter un mélange ternaire quelconque (A,
B, C) (Minotti et al., 1996, Wuithier and Giraud, 1965, Wuithier and Giraud, 1972a,
Wuithier and Giraud, 1972b).
39
• Bilan matière global
; + 7 = <= + %>
; + <? = <= + %=
Ces bilans peuvent être transformés en écrivant la différence entre deux courants
à un niveau donné de la colonne :
Les bilans matières partiels et globaux peuvent être traduits en une relation de
colinéarité permettant de déterminer le nombre d’étage théorique, la règle du bras de
levier.
; + 7 = <= + %> = E
40
Chaque courant est matérialisé par son débit et par ses compostions de chaque
constituants. Les courants peuvent donc être placés sur le diagramme liquide-liquide au
moyen des compositions des constituants.
Figure I-23: Succession de droites d’équilibres et de droites opératoires pour la détermination du nombre d’étages
théoriques
41
I.4.3 Technologies
42
Figure I-24: Organigramme de présélection d’un extracteur liquide-liquide (Technique de l’ingénieur J2766)
Les appareils utilisés dans les procédés d’extraction par solvant sont des
contacteurs liquide-liquide dans lesquels le mélange intime des deux phases crée une
aire interfaciale à travers laquelle le soluté peut se transférer. On distingue deux grandes
classes d’appareil :
Dans les premiers, une phase est dispersée dans l’autre et amenée à l’équilibre
dans le mélangeur, puis les phases sont séparées par décantation et envoyées à contre-
43
courant dans les étages adjacents. On distingue notamment les mélangeurs décanteurs
et les extracteurs centrifuges.
44
On distingue dans cette classe les colonnes gravitaires qui sont des colonnes
simples à concevoir avec un prix de revient faible. Cependant, l’absence de moyen
d’agitation ne permet pas de maintenir une aire interfaciale d’échange en garantissant
l’écoulement à contre-courant. Les colonnes rencontrées industriellement ont des étages
théoriques de 1,5 à 10 m et les débits spécifiques totaux sont compris entre 10 et 100
m3h-1m-2 (Technique de l’ingénieur J2766). On retrouve dans les colonnes différents
garnissages permettant d’accroître l’aire interfaciale et le mélange axial. Le Tableau I-6
montre les modèles de colonnes rencontrés industriellement ainsi que leur débit
surfacique et quelques remarques sur l’unité (Stichlmair, 1980).
Débits
Colonnes Remarques
surfaciques
45
Débits
Colonnes Remarques
surfaciques
46
Colonnes à disques rotatifs (RDC)
Colonnes Rushton
Colonnes pulsées
Colonnes Kühni
47
I.4.4 Extraction liquide-liquide réactive
48
Deux méthodes de design sont examinées dans cette partie, les méthodes
conceptuelles proposées par Minotti et al (1998) et Samant et Ng, (1998b), puis la
méthode géométrique d’Olán-Acosta et al, (2014).
Les méthodes conceptuelles reposent sur l’écriture des bilans le long de la colonne. À
partir des compositions et des débits des alimentations en mélange réactionnel et en
agent de séparation, il s’agit de calculer le design de la colonne qui permettra d’atteindre
les objectifs de séparation en termes de compositions ou de taux de récupération dans
l’extrait et le raffinat.
a. Minoti et al (1998)
49
liquide-liquide sont calculés en fonction de la composition du mélange de la température
et de la pression.
>
* = N *D 1.6
DQ=
Où +DT représente la vitesse de réaction. Cette réaction est valable pour tous les
composés participants à la réaction.
Pour chaque étage théorique n, on compte Nc bilans matière partiels (équation 1.2),
Nc relations d’équilibre liquide-liquide (1.3), Nc équations sur la production molaire et
l’avancement de la réaction (1.5) et une équation sur la rétention liquide (1.6). A cela
s’ajoute les 2(N+1) équations de sommation. Ce système à pression et température
fixées pour la colonne entière est constitué de :
• 3Nc*N+N+2(N+1) équations
• 2Nc(N+1)+Nc*N+2(N+1)+N+Nr*N+1 : N, S, xS, F, xF, En, yn, Rn, xn, Hn, Pn et ζnj avec
n=1, 2, …, N et j=1, 2, …, Nr.
50
Il y a donc 2Nc+NrN+1 degrés de liberté. Les alimentations en mélange réactionnel et
en solvant sont connues, S, F, xS et xF sont des données du problème. Il reste donc à
spécifier NrN+1 paramètres pour calculer complètement le système. Minotti et al (1998)
ont choisi de se fixer l’avancement de chaque étage théorique et le nombre N d’étage
théorique. Pour les calculs une méthode pas à pas est employée, à partir des valeurs aux
bornes du système S F xS xF.
l’équation 1.5.
− Il faut ensuite déterminer les profils de composition de chaque constituant le long
de la colonne ainsi que les débits. Pour cela les bilans matières, les équilibres entre
phases et les équations de sommation sont résolus. Pour la colonne à contre-courant une
méthode itérative est utilisée. A partir d’une estimation de la composition F>I (ou F=II ), les
profils dans la colonne sont calculés. Le critère d’arrêt est défini pour qu’il s’active
lorsque la composition F=II (ou F>I ) calculée vérifie le bilan matière global.
− Enfin avec l’équation 1.7, la rétention liquide est calculée.
Cette méthode de design développée par Minotti et al (1998) n’est à notre sens pas
exactement une véritable méthode de design. Le fait est que pour calculer le profil de
compositions dans la colonne, le nombre d’étage théorique N et l’avancement dans
chaque étage UDT qui sont les paramètres à déterminer dans une méthode de conception,
b. Samant et Ng (1998b)
51
sont le raffinat et l’extrait représentés respectivement par les débits R et E, et les
compositions F>I et F=II .
La méthode proposée par Samant et Ng (1998) fait appel aux compositions réactives
décrites dans la partie II.2.3.3.
Maintenant le terme relatif à la réaction est éliminé des équations 1.8 et 1.9 en
utilisant l’équation 1.10 pour obtenir :
7_ + ;^ = <^ + %^ 1.17
7_`AH + ;^ `AG = <
a= `A,=
II
+ %a> `A,>
I
1.18
De la même manière, le bilan matière réactif autour d’un seul étage s’écrit :
b
< b I a II a I
D@= `A,D@= + %>B= `A,>B= = <D `A,D + %D `A,D
II
1.19
52
Les courants sortant de chaque étage i sont à l’équilibre entre phases :
`A,D
II
= 2Zc (`A,D
I
) 1.20
>P >P@>V
N FA,=
I
=1 N `A,D
I
=1 1.21
AQ= AQ=
>P >P@>V
N FA,>
II
=1 N `A,D
II
=1 1.22
AQ= AQ=
53
méthode est en fait une adaptation de la méthode graphique pour les systèmes non
réactifs, dans le mode de représentation qu’il utilise et qui est illustrée au II. 3. b. ii.
Les bilans matière sur la colonne d’extraction réaction sont écrits via les
différences entre les courants entrants et sortants :
Le point P est constant pour tous les étages N et fait office de pôle dans la
construction graphique de la colonne. Ses coordonnées dans le diagramme liquide-
liquide réactif sont déterminées à partir des bilans matière partiels.
Une fois ce pôle placé sur le diagramme liquide-liquide réactif, le nombre d’étages
théoriques est déterminé de la même manière que celle établie pour les systèmes non-
réactifs. A partir du point représentant l’extrait E1 dans le diagramme, on alterne
conodales et droites opératoires jusqu’à ce qu’une conodale permette d’obtenir un
raffinat de composition en soluté inférieure à celle de Rn. La procédure de cet
enchainement est décrite dans les 5 étapes ci-dessous
54
Figure I-27 : Construction graphique pour une colonne d'extraction liquide-liquide réactive (Olán-Acosta et al., 2014)
La valeur de ce taux de solvant est bornée donc par des valeurs minimales et
maximales. Ces deux taux sont désignés par un taux inférieur et maximum. Cet intervalle
comprend toutes les valeurs possibles du taux de solvant et est régi par l’équilibre
liquide-liquide réactif. Toutes les valeurs de cet intervalle ne garantissent pas pour
autant la faisabilité du procédé d’extraction liquide-liquide réactive. Le taux de solvant
minimum est introduit et correspond à la valeur clé de cet intervalle. Si la valeur du taux
de solvant opératoire est comprise entre le taux de solvant inférieur et le taux de solvant
minimum, l’objectif de séparation sera atteint en un nombre infini d’étage théorique. Si
la valeur du taux de solvant opératoire est comprise entre le taux de solvant minimum et
le taux de solvant maximum, l’objectif de séparation sera atteint en un nombre fini N
d’étage théorique.
55
Les études de faisabilité faites par les différents auteurs, consistent à déterminer ce
taux de solvant minimum.
a. Samant et Ng (1998b)
Au point fixe donc, les compositions `^ et k^ sont à l’équilibre entre phases et les
bilans matières s’écrivent :
Où <^l 10 %^l sont les débits molaires réactifs de chaque phase au point fixe. Ces
trois équations sont réécrites pour faire apparaitre le rapport S/F, puis résolues pour
déterminer le taux de solvant minimum.
56
7 ` G − ` eP
m n = eP
; AD' ` − `H
1.30
7 `G − `p
m n = p
; oAD ` − `H
1.31
Figure I-28 : Limite de faisabilité d'une colonne réactive Figure I-29 : Détermination graphique du taux de solvant
à contre-courant minimum
7 ` G − `Jqrs
m n =
; oAD `Jqrs − `H
1.32
57
appelés à être de plus en plus présents dans le paysage industriel chimique. Le procédé
de distillation réactive est largement étudié dans ce contexte de développement durable,
et une méthodologie de conception a été développée dans les travaux de Théry et al
(2002). Cependant la distillation réactive peut être contre-indiquée pour traiter des
milieux fermentaires qui on le rappelle sont des milieux très dilués.
Le procédé qui nous intéresse dans cette étude est l’extraction liquide-liquide
réactive. Couplées au sein d’une même unité, l’extraction liquide-liquide et la réaction
d’estérification représentent le procédé intensifié retenu pour la récupération d’acides
carboxyliques produits par voie fermentaire. L’acide succinique qui est un acide
dicarboxyliques a été choisi comme molécule modèle.
58
Chapitre II
Une méthodologie pour la
conception d’un procédé d’extraction
liquide-liquide réactif
59
Chapitre II Une méthodologie pour la conception d’un procédé
d’extraction liquide-liquide réactif
ANALYSE DE
FAISABILITÉ
BILAN MATIÈRES SUR LA COLONNE
D’EXTRACTION LIQUIDE-LIQUIDE RÉACTIVE
SIMULATION ET OPTIMISATION
VALIDATION
ESSAIS EXPERIMENTAUX
61
− les données sur les réactions chimiques, coefficients stœchiométriques et
cinétiques chimiques pour chaque réaction,
− les objectifs de séparation en termes de taux de conversion ou taux de
récupération et composition du produit final.
Le dimensionnement de colonne est donc réalisé par le biais des trois étapes
hiérarchisées. L’analyse de faisabilité, l’étape de synthèse et l’étape validation
composent la méthode de conception d’une colonne d’extraction liquide-liquide réactive.
Il est introduit au fur et à mesure de la démarche de conception des hypothèses sur les
phénomènes physiques et chimiques de plus en plus complexes.
62
La méthode de synthèse pour l’extraction liquide-liquide simple est ici adaptée au
système réactif et la méthode détaillée pour les systèmes réactifs est exploitée avec les
données de l’analyse de faisabilité.
Simulation de
la colonne à
Données Résultats l’échelle
Outils de Essais
synthèse, Outils de expérimentaux sur
méthode simulation colonne à échelle
Short Cut laboratoire
Amélioration
Transfert de du modèle
données (NET,
profils de
composition)
Résultats
Figure II-2 : Processus de développement d'un outil de simulation et relation entre essais expérimentaux et méthode
Short Cut.
63
II.2 Thermodynamique des systèmes réactifs équilibrés entre phases
>P
>V
64
Où UT représente l’avancement de la réaction chimique et dans le cas particulier
d’une réaction équilibrée, désigne plus précisément un état de l’équilibre chimique. En
intégrant cette équation 2.2, on en déduit l’expression du nombre de moles 4A de
l’espèce i au cours de la réaction en fonction du nombre de moles initial de l’espèce i
dans le système et de l’avancement de la réaction chimique :
>P
Pour une valeur particulière de l‘avancement U = UZc , la situation est celle d’un
état d’équilibre. Ce dernier est unique et ne dépend pas de la façon dont on y est
parvenu. La condition d’équilibre - X,e w = 0 implique que la relation 2.4 est satisfaite
pour la valeur particulière de l’avancement U = UZc .
∆V wX,e \UZc ] = 0 (2.5)
>P
AQ= AQ=
>P
∆V w : (z) = N SA yA:
AQ=
65
Donc d’après la relation ∆V wX,e (U) = 0 :
>P
Ce qui nous donne l’expression du quotient réactionnel qui dans le cas particulier
de l’équilibre chimique prend le nom de constante d’équilibre.
>P
∆V w : (z)
€(z) = 1F• ‚− ƒ = } 5A,Zc
~r
%z
(2.8)
AQ=
>P
>l >P
66
De cette formulation de l’enthalpie libre, on peut déduire que dans un système à
l’équilibre physique, les potentiels chimiques de chaque constituant sont égaux dans les
différentes phases où ils sont présents.
2A = MA FA 2A
‹ ‹ ‹ Œr
(2.15)
67
chimique, l’équation 2.20 de l’équilibre physique pourra être utilisée dans ces systèmes
d’équations non linéaires dans le cadre de calculs d’équilibre chimique entre phases.
• Les fractions molaires de chaque constituant dans chacune des 2 phases : FAI et
FAII
Cela représente (2 ∙ …P ) inconnues.
68
• Les avancements de chaque réaction : UV
Cela représente (…V ) inconnues.
FAI
Fraction molaire phase I
∀• ∈ )‘, ’“ 6
FAII
Fraction molaire phase II
∀• ∈ )‘, ’“ 6
UV
Avancement de réaction
•V =
…X:
adimensionnel
∀” ∈ )‘, ’” 6
•I
Taux de partage
molaire de phase I
Tableau II-1 : Variables d’état du système liquide-liquide réactif
Avec FA: la fraction molaire du constituant / dans le mélange initial, et •I et •II les
taux de partages molaires de la phase I et de la phase II. On a donc : •I = 1 − •II . …A
lV–&
>™
69
On a donc le bilan de matière total qui donne, puisque l’on a ∑AQ= FAI = ∑AQ= FAII =
› › > >
1 :
>™ >›
>™
FA: + N\˜A,V ∙ •+ ]
VQ=
>™ >›
i=1,2,…,Nc (2.20)
= ¡•I ∙ FAI + (1 − •I ) ∙ FAII ¢ ∙ œ1 + N ••V ∙ N ˜T,V Ÿ
VQ= TQ=
Avec MAI et MAII les coefficients d’activité du constituant / dans les phases I et II, FAI
et FAII les fractions molaires du constituant / dans les phases I et II.
70
Type
Formulation des équations
d’équations
>™
∙ œ1 + N •UV ∙ N ˜T,V Ÿ = 0
constituant
∀• ∈ )‘, ’“ 6
VQ= TQ=
Eq. chimique Pour des raisons de stabilité numérique (éviter les divisions par 0), on
∀” ∈ )‘, ’” 6 reformule :
>› >›
La résolution numérique ainsi que le tracé des courbes sont réalisés dans un
environnement Matlab. Les calculs thermodynamiques (évaluation des coefficients
d’activités MA,T ) sont réalisés en faisant appel à Simulis (en utilisant la toolbox Simulis
pour Matlab). Pour résoudre le système d’équations algébriques en chaque point
(modèle de flash LL réactif), on utilise la fonction FSOLVE de Matlab, basée sur
l’algorithme de « trust-region dogleg » (Powell, 1970).
71
La variance V d’un système correspond au nombre minimal de variables
intensives indépendantes qu’il faut fixer arbitrairement pour que l’état du système soit
entièrement déterminé. Considérons toujours le même système comprenant Nc
constituants Nr réactions et Np phases en équilibre. Pour décrire ce système, les
fractions molaires de chaque constituant dans chaque phase doivent être précisées ainsi
que sa température et sa pression, soit …† ∙ …• + 2 variables intensives. Ces variables ne
sont pas totalement indépendantes. La somme des fractions dans chaque phase est égale
à 1, ce qui représente Np relations. Le système étant à l’équilibre de phases et chimique,
l’égalité des activités donne …†(…• − 1) + …+ relations supplémentaires. La variance
s’exprime donc de la manière suivante :
¦ = …† ∙ …• + 2 − …• − …†(…• − 1) − …+
¦ = …† − …+ − …• + 2
(2.23)
¦ = …† − …+ (2.24)
La variance d’un système réactif est inférieure à celle d’un système non réactif.
Cette caractéristique permet de représenter graphiquement plus facilement les
systèmes réactifs avec plus de 3 constituants. Pour cela, il convient d’utiliser les
compositions réactives détaillées dans le paragraphe suivant.
72
de phases en absence de réaction et sont par conséquent définies comme les
coordonnées naturelles exprimant les compositions dans un système en équilibre
thermodynamique et en présence de réactions chimiques.
S(>P@>VB=),= ⋯ S(>P@>VB=),>V
[=| ⋮ SA,>V ⋮ •
S>P,= ⋯ S>P,>V
(2.27)
FA − SAX [ @= FVZ'
`A =
1 − SXYX
X
[ @= FVZ'
i=1,2,…,Nc (2.29)
Notons que les composés de références ne doivent pas être choisis de manière
aléatoire, mais ils doivent être choisis de manière à ce que la matrice [ soit inversible.
Le choix des composés de référence doit permettre d’avoir un déterminant de la matrice
[ non-nul. Cela implique les contraintes suivantes :
73
• aucune des lignes ou des colonnes de la matrice ne peut être égale à zéro. C’est-à-dire
que chaque réaction chimique doit mettre en jeu au moins un constituant de
référence. Un composé inerte ne peut pas être choisi comme composé de référence,
• deux composés de référence ne peuvent pas avoir les mêmes coefficients
stœchiométriques dans chaque réaction,
• afin d’éviter que le dénominateur de l’expression 2.29 ne soit nul, les règles suivantes
devront être suivies :
o Si SXYX
X
>0, choisir un constituant de référence parmi les réactifs
o Si SXYX
X
<0, choisir un constituant de référence parmi les produits
o Si SXYX
X
=0, choisir indifféremment un constituant de référence parmi les
réactifs ou les produits.
L’intérêt de ce nouveau formalisme réside dans l’invariance par rapport à la
réaction chimique : les compositions réactives expriment la conservation de la masse à
chaque valeur d’avancement de la réaction. Elles expriment ainsi les compositions du
système comme si aucune réaction n’avait eu lieu. A une composition réactive
correspond plusieurs compositions réelles à différentes valeurs de l’avancement, il en
résulte que :
L’autre intérêt est que la somme des compositions réactives X et Y est égale à 1
comme pour les compositions réelles x et y.
>P@>V
N `A = 1 (2.31)
AQ=
74
pour déterminer le domaine de composition réactif, celle-ci étant représentée ci-
dessous.
75
− `±« = F±« + Fe«
− `±¬ = F±¬ + Fe«
− `e¬ = Fe¬ − Fe¬
4. Positionnement des sommets correspondant aux espèces pures :
Calculons les compositions réactives X=(`±« , `±¬ ) pour chaque espèce pure :
− %? pure : X= (1, 0)
− *= pure : X= (1, 1)
− *? pure : X= (0, 0)
Le domaine de composition réactive peut alors être représenté sur la Figure II-3 :
t= + t? ↔ t· + t¸
t¸ + t? ↔ t· + t¹
76
En appliquant comme précédemment la procédure de détermination du domaine
de composition réactive, on obtient le domaine présenté dans la Figure II-4 mais avec
une particularité : l’espèce A4 n’est pas située sur un sommet. L’espèce A4 est à la fois
produite dans la première réaction et réactive dans la seconde. Elle ne peut donc exister
à l’état pur dans la solution puisque dès son apparition elle donne lieu à une autre
réaction.
Figure II-4 : Domaine de composition réactive pour l'exemple de deux réactions successives.
77
On applique la procédure permettant de tracer le diagramme liquide-liquide
réactif. Le constituant de référence choisis ici est l’acétate de butyle.
Il en résulte le diagramme d’équilibre liquide-liquide réactif présenté ci-dessous.
Tous les points de ce diagramme sont à l’équilibre chimique. Si le mélange considéré est
dans la zone monophasique, le milieu réactionnel reste homogène. Si le mélange
considéré est dans la zone diphasique, il démixte en deux phases représentées par les
points rouges (raffinat) et bleus (extrait) situés aux extrémités des conodales (droites
bleues).
L’équilibre liquide-liquide qui s’exprime selon les compositions et les coefficients
d’activité, FAI MAI − FAII MAII = 0, induit l’équilibre chimique dans les deux phases.
>P >P
Figure II-5: Représentation en coordonnées réactives des équilibres chimiques entre phases pour le système Acide
Acétique/Butanol/Eau/Acétate de Butyle.
78
enfin une étape de simulation validation. Dans chaque partie seront développés les
outils et les équations décrivant le système sous différentes hypothèses simplificatrices
et permettant d’avancer dans la méthode et dans la détermination du dimensionnement
de la colonne d’extraction liquide-liquide réactive.
Toutes ces données calculées lors de l’analyse de faisabilité sont des paramètres
essentiels pour l’étape de dimensionnement de la colonne.
79
II.3.1.1 Introduction des courbes de raffinats réactifs : les courbes de résidus
réactifs en distillation réactive
Afin de mieux comprendre l’intérêt des courbes de raffinats réactifs, on
s’intéresse dans un premier temps aux courbes de résidus réactifs qui sont analogues
aux courbes de raffinats, mais appliquées à la distillation réactive.
Les courbes de résidus réactifs ont été développées par Barbosa et Doherty
(1988) pour un mélange soumis à une réaction équilibrée traité par distillation simple
(type distillation de Rayleigh). Plusieurs auteurs (Ung et Doherty 1995, Brehelin 2006,
Niang et Mikitenko 1998, Pham and Doherty, 1999) reprennent ces travaux pour
développer le modèle décrivant la distillation réactive simple d’un mélange soumis à
plusieurs réactions chimiques équilibrées et à la distillation réactive liquide-liquide
vapeur. Ils calculent les courbes de résidus réactifs ainsi que les caractéristiques de
l’équilibre chimique entre phases, des limites de séparations d’une façon relativement
simple.
Une courbe de résidu réactif est le lieu des points de compositions liquides en
équilibre avec la phase vapeur lors d’une distillation simple mettant en jeu un système
soumis à une ou plusieurs réactions équilibrées.
80
Le réseau de courbes de résidu réactif permet entre autres de déterminer les
points singuliers et les frontières de distillation du système.
-`A
= `A − kA
-¿
i=1, 2, …, Nc-Nr (2.32)
N OA = 1 (2.35)
AQ=
Tout comme les courbes de résidu réactif reposent sur la distillation simple, les
courbes de raffinat réactif reposent sur une extraction simple comme présentée sur la
Figure II-7. C’est un procédé d’extraction à un étage théorique à rétention constante. On
retrouve une déverse matérialisée par le courant E et une alimentation en solvant frais.
81
Le liquide dans le décanteur est à l’équilibre liquide-liquide et réactif et chimique.
L’alimentation de solvant est considérée comme pure en agent extractant et un courant
de phase légère (courant E) est récupéré. Les courbes de raffinat réactif peuvent alors
être calculées pour les deux phases et représentent les profils liquides de chaque phase
dans une extraction liquide à courants croisés à taux de solvant infini.
Figure II-7 : Procédé d’extraction liquide-liquide réactive simple, décanteur réactif pour le calcul des courbes
de raffinat réactif.
82
Les Nc bilans matières sur chaque constituant sont exprimés de la façon
suivante :
>V
-ÀA -UV
W − ;ÂA + <FAII + N WSA,V =0
-0 -0
i=1,…,Nc (2.36)
AQ=
L’algorithme développé sous Matlab utilise une fonction “ODE” (Shampine and
Reichelt, 2006) qui résout les systèmes d’équations différentielles de la forme y′=f(t, y).
Il est nécessaire de formuler le système d’équations de manière à n’avoir qu’un seul
terme différentiel par équation.
>V
ÀA − ÁFAI
FAII =
1−Á
i=1,…,Nc (2.38)
;0 ;-0
¿= -¿ =
W W
(2.39)
-3A
− ÂA + <Ã; FAII = 0
-¿
i=1,…,Nc (2.40)
>P
83
L’équation de sommation sur la variable F II :
>P
1 − N FAII = 0 (2.43)
AQ=
N 3A − W |1 − N N SA,V UV • = 0 (2.44)
AQ= AQ= VQ=
Durant l’étude des courbes de raffinat réactive, les équations ont été écrites en
fonction des compositions réactives (annexe B). Le programme pour résoudre le
système d’équations en fonction des compositions réactives n’a pas été réalisé. Le
principal problème rencontré lors de la création de ces programmes (en fonction des
compositions normales et des compositions réactives) est la possibilité de se retrouver
avec un état monophasique. Afin de pallier à ce problème, une méthode a été développée
dans le cadre de la distillation réactive et pourrait être transposée à notre cas
d’extraction liquide-liquide réactive, il s’agit de la méthode des taux (Brehelin, 2006). Le
but de cette méthode est de gérer la présence ou l’absence d’une phase, elle se base sur
la propriété qu’ont les phases non existantes physiquement de voir la somme de leurs
fractions molaires des constituants être inférieure à l’unité. Ainsi en ajoutant une
variable de phase notée pour chaque phase prenant la valeur 0 dans le cas d’une phase
existante et la valeur ¿ = 1 − ∑>P
AQ= FA lorsque la phase n’est pas présente. Cette nouvelle
84
Les résultats de ces calculs de raffinat réactif sont représentés sur le diagramme
liquide-liquide réactif. A partir d’une composition global w de la rétention, l’intégration
est réalisée dans le sens direct et indirect (Figure II-8). Les courbes de raffinat réactif
comme les courbes de résidu réactif permettent de déterminer les zones d’extraction
respectivement de distillation et l’identification de frontières de séparation. Même si
avec les courbes de raffinat réactif on n’identifie pas de zone d’extraction, elles nous
indiquent tout de même que toutes les compositions atteignables sont sur l’enveloppe
de phase. Les courbes de raffinat réactif ne donnent pas plus d’information que le
diagramme liquide-liquide réactif dans le cas où la représentation du système est en
deux dimensions. Nous verrons plus tard dans l’application de cette méthode de
dimensionnement que les courbes de raffinat réactif sont vraiment utiles pour des
systèmes avec une variance supérieure à 3.
85
solvant peuvent être déterminées avec la règle du bras de levier connaissant
parfaitement le comportement de l’équilibre chimique entre phases du système avec le
diagramme liquide-liquide réactif.
Figure II-9 : Quantités minimales et maximales de solvant qui garantissent la présence de deux phases.
86
(a) (b)
Figure II-10: Schéma d’une colonne d’extraction liquide-liquide réactive à contre-courant multi-étage (a)
simple et (b) à reflux.
87
La détermination du taux de solvant minimum est réalisée de manière analogue à
celle effectuée pour l’extraction liquide-liquide simple. La conodale passant par le point
du diagramme liquide-liquide réactif représentant l’alimentation correspond à l’extrait
le plus riche que l’on peut obtenir avec une colonne à contre-courant (<Ä ). Sur la Figure
II-11, le croisement entre les droites alimentation-solvant (droite F-D) et la droite entre
la phase aqueuse et la phase organique (droite RN-El) définit le point ML qui permet de
calculer le taux de solvant minimum avec le système d’équations suivant :
; + 7 = EÄ
;•A + 7`A = EÄ `JÅ ,A
;•A + 7`A = (; + 7)`JÅ,A
;(•A − `JÅ ,A ) = 7(`JÅ,A − `A )
7 (•A − `JÅ,A )
m n =
(2.45)
; oAD (`JÅ ,A − `A )
7 7 7
m n = 1.6 ∗ m n ↔ 7 = ; ∗ 1.6 ∗ m n
; ; oAD ; oAD
B C
El
F
ML
RN
A S
Figure II-11 : Représentation graphique de la détermination du taux de solvant minimum.
88
II
débit d’alimentation (z, F) et la composition du solvant (F>B= ). Une autre donnée à notre
disposition correspond à la composition de la phase lourde (F>I ) fixée avec l’analyse des
courbes de raffinat réactives.
FA,=
I
∙ MÈ/ − FA,=
II
∙ MÈÈ
/ =0
i=1,…,Nc (2.49)
>P
N FÈÈ/,1 − 1 = 0 (2.50)
AQ=
Pour la colonne à reflux, on retrouve les mêmes équations, bilan matière global,
bilan matière partiel et équations d’équilibre chimique entre phases, plus les bilans
matière globaux et partiels sur chaque unité autre que la colonne, à savoir le séparateur,
le reflux et le mélangeur de solvant en pied de colonne. On a ainsi un système de 4Nc+9
équations avec autant de variables inconnues.
89
• Colonne
7 − %> − <= + ; + %: = 0 (2.51)
FA,=
I
∙ MÈ/ − FA,=
II
∙ MÈÈ
/ =0 i=1,…,Nc (2.54)
>P
N FÈÈ/,…+1 − 1 = 0 (2.55)
AQ=
>P
N FÈÈ/,1 − 1 = 0 (2.56)
AQ=
• Séparateur
<= − %: − 71 − *1 = 0 (2.57)
N F/,71 − 1 = 0 (2.59)
AQ=
• Reflux
Pour le reflux on définit le taux de reflux ¿V comme le rapport du débit %: par le
dé bit *1 et est une donnée du système.
¿V *1 − %: = 0 (2.60)
FA,eZ − FA,ŒÉ = 0 (2.61)
• Mélangeur
7`ÈÈ…+1 − 71`/,71 − 75`/,75 = 0 i=1,…,Nc (2.62)
90
II.3.2 Synthèse et dimensionnement : détermination du nombre
d’étages théoriques
Dans cette étape de notre méthode de conception, l’objectif est de déterminer le
nombre d’étages théoriques suivant l’organigramme présenté sur la Figure II-12.
Donné es du systè me
• Modèle Thermodynamique
• Composition et débits des débits
entrants et sortants du procédé (F,
S, E1, RN)
j=2
données : FTI
• j=j+1
•
Oui
NET=j
91
II.3.2.1 Saisie des données de conception
Les paramètres suivants servent de données de base à l’étape de synthèse du
procédé :
92
Figure II-13 : Calcul du bilan de matière pour déterminer le nombre d'étages théoriques.
FA,T
I
∙ MÈ/ − FA,T
II
∙ MÈÈ
/ =0 i=1,…,Nc (2.63)
N FÈ/,Í − 1 = 0 (2.66)
AQÊ
P
N FÈÈ/,Í − 1 = 0 (2.67)
AQÏ
Sur chaque étage où les compositions sont calculées, un test est réalisé avec un
critère d’arrêt afin de déterminer le nombre d’étages théoriques. Ce test est détaillé dans
le paragraphe II.3.2.4 suivant.
93
II.3.2.3 Génération des profils de composition : cas d’une colonne à reflux
Considérons donc maintenant le système d’extraction liquide-liquide réactif à
reflux. A partir d’une alimentation connue et des valeurs aux frontières, et supposant
que l’équilibre liquide-liquide est atteint sur chaque étage, les jeux d’équation définit ci-
dessous permettent respectivement de modéliser les zones d’enrichissement et
d’appauvrissement.
<T ∙ `A,T
II
− %T@= ∙ `A,T@=
I
− <= ∙ `A,=
II
+ %: ∙ `A,ŒÉ = 0
i=[1,Nc-Nr] j=[1,f]
>P
N FA,T
I
−1=0
AQÊ
P
N FA,T
II
−1=0
AQÏ
<T ∙ `A,TB=
II
− %T@= ∙ `A,T@=
I
− %> ∙ `A,>
I
+ 7 ∙ `A>B= = 0
i=[1,Nc-Nr] j=[f,N]
>P
N FA,T
I
−1=0
AQÊ
P
N FA,T
II
−1=0
AQÏ
94
II.3.2.4 Le critère d’arrêt et l’analyse des résultats
La procédure de détermination du nombre d’étages théoriques est complète
lorsque la composition de la phase II (F=II ) fixée lors de l’analyse de faisabilité est
atteinte. La composition du constituant désiré est alors choisie comme critère d’arrêt.
Lorsque la composition du profil de composition de la phase deux atteint cette valeur, la
procédure est complète et on peut alors déterminer le nombre d’étages théoriques.
0,04
0,035
0,03
Composition Phase II
du constituant A
0,025
0,02
0,015
Courbe d'équilibre
0,01 Courbe opératoire
Construction type McCabe et Thiele
0,005
Critère d'arret
0
0 0,05 0,1 0,15 0,2
Composition phase I du constituant A
Figure II-15 : Construction type McCabe et Thiele pour illustrer le critère d’arrêt de la détermination du nombre
d’étages théoriques.
95
0,10
0,09
Composition phase II du constituant A
0,08
0,07
0,06
0,05
0,04
enrichissement distribution
0,03
appauvrissement distribution
0,02 enrichissement opératoire
Alimentation appauvrissement opératoire
0,01
0,00
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1
Composition phase I du constituant A
0,10
0,09
Composition phase II du constituant A
0,08
0,07
0,06
0,05
0,04
0,03
Critère d'arret
0,02 Distribution
Opératoire
0,01
0,00
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1
Composition phase I du constituant A
Figure II-16 : Construction type McCabe et Thiele pour illustrer le critère d’arrêt de la détermination du nombre
d’étage théorique et la position de l’alimentation, colonne à reflux.
96
II.3.3 Conception et validation
>› >›
Õ >› >› Ø
ÑI
•V = Ô (MA ∙ FA ) r,™ − Ð@V ∙ } (MAI ∙ FAI )@£r,™ ×
: ∙ ¿ ∙ ÐV ∙ }
I I £
Ò0é X Ô ×
r=1, 2, …, Nr (2.69)
AQ= AQ=
Ó £r,™ ¥: £r,™ ¤: Ö
• Si réaction en phase I et II :
Õ >› >› Ø
ÑX
•V = ∙ ¿ ∙ ÔÐV ∙ } (MAI ∙ FAI )£r,™ − Ð@V ∙ } (MAI ∙ FAI )@£r,™ ×
Ò0é:X Ô ×
r=1, 2, …, Nr (2.70)
AQ= AQ=
Ó £r,™ ¥: £r,™ ¤: Ö
•I ∗ Ò0é:X ∗ E
¦I =
Ù
(2.71)
Ò0é:X ∗ E
¦X =
Ù
(2.72)
On obtient finalement :
Õ >›
>› Ø
∙ } (MAI ∙ FAI )@£r,™ ×
AQ=
×
£r,™ ¤: Ö
Õ >› >› Ø
E
•V = ∙ ¿ ∙ ÔÐV ∙ } (MA ∙ FA )£r,™ − Ð@V ∙ } (MAI ∙ FAI )@£r,™
I I ×
Ù Ô ×
r=1, 2, …, Nr (2.74)
AQ= AQ=
Ó £r,™ ¥: £r,™ ¤: Ö
98
Type
Expressions des résidus
d’équations
>™ >™ >›
Õ >› >› Ø
E
•V − ∙ ¿ ∙ ÔÐV ∙ } \MA ∙ FA ] ×=0
Cinétique
£r,™ @£
I I
− Ð@V ∙ } \MAI ∙ FAI ] r,™
Ù Ô ×
chimique
AQ= AQ=
∀” ∈ )‘, ’” 6 Ó £r,™ ¥: £r,™ ¤: Ö
Tableau II-3: Système d’équations de l’équilibre liquid-liquide réactif hors équilibre chimique
Pour chaque étage j, les variables sont décrites dans le tableau 1 et les données du
système sont regroupées dans le tableau 2.
<Í
Débit molaire de la phase II sortant de mol.s-
l’étage j 1
%Í
Débit molaire de la phase I sortant de mol.s-
l’étage j 1
99
FÈÍ Composition molaire courant %Í Nc
FÈÈ
Í Composition molaire courant <Í - Nc
%A,T@= ∙ FA,T@=
I
+ <A,TB= ∙ FAB=
II
− %T ∙ FA,T
I
− <T ∙ FA,T
II
+ àA,T = 0
i=1, 2, …, Nc
(2.75)
j=1, 2, …, NET
MAI ∙ FA,T
I
− MAII ∙ FA,T
II
=0
i=1, 2, …, Nc
(2.76)
j=1, 2, …, NET
1 − N FA,T
I
=0 (2.77)
1 − N FA,T
II
=0 (2.78)
>V Õ >›
ÑX
àA,T = Ô I £r,™
: ∙ ¿ ∙ N ÐV ∙ } \MA,T ∙ FA,T ] − Ð@V
I
Ò0é X Ô
VQ= AQ=
Ó £r,™ ¥:
Ø
r=1, 2, …, Nr (2.79)
>›
∙ } \MA,T
I
∙ FA,T
I
]
@£r,™
×
AQ=
×
£r,™ ¤: Ö
100
II.4 Résumé de la méthode
Dans ce chapitre, nous avons présenté en détail la méthodologie de conception
d’une colonne d’extraction liquide-liquide réactive développée durant cette thèse. Cette
méthode ne se limite pas aux systèmes avec un nombre de constituants et de réactions
limites, ni sur des interprétations graphiques. En effet, comme nous l’avons montré avec
la notion de variance, celle-ci doit être inférieure ou égale à trois pour avoir une
représentation en deux dimensions. Or dans notre exemple choisi (les systèmes
diphasiques réactifs avec de l’acide succinique, et suivant la stratégie de séparation, avec
ou sans inerte), la représentation graphique peut être plus complexe à interpréter.
La méthode de conception s’appuie sur trois étapes majeures où des outils ont été
développés afin d’aboutir aux paramètres structuraux et opératoires d’une colonne.
Dans ce chapitre nous avons adapté des outils au système réactif, et introduit un nouvel
outil comme les courbes de raffinat réactif issu de l’analogie avec les courbes de résidu
réactif en distillation réactive.
101
DIAGRAMME LIQUIDE-LIQUIDE RÉACTIF Taux de solvant infini
et
COURBES D’EXTRACTION RÉACTIVE Colonne entièrement réactive
Introduction de la cinétique
SIMULATION ET OPTIMISATION Et
Étude de sensibilité
Figure II-17: Organigramme détaillé de la méthode de conception d’une extraction liquide-liquide réactive
102
Chapitre III
Application à la séparation d’un
acide carboxylique issu d’un milieu
fermentaire : l’acide succinique
La seconde stratégie consiste à prendre un alcool léger, ici le méthanol, et d’utiliser un tiers
corps comme agent extractant. Les composés sélectionnés ici comme agent extractant sont
des hydrocarbures lourds comme le para-xylène ou le n-octane.
103
Chapitre III Extraction liquide–liquide réactive appliquée à l’acide
succinique biosourcé.
Dans cette étude de cas, la méthodologie de conception est appliquée à des systèmes
réactifs complexes. Les différents systèmes étudiés sont composés de l’acide succinique
d’eau et de l’alcool choisi, des produits de la réaction d’estérification entre l’acide
succinique et l’alcool et d’un tiers corps suivant la stratégie d’extraction envisagée. Le
Tableau III-1 suivant regroupe toutes les espèces chimiques mises en jeu dans cette
étude.
º¸ Wá »¸
Acide Succinique
118.09
(SA)
º¹ Wâ »¸
Mono-Méthyl
132.11
Succinate (MMS)
ºá W=: »¸
Di-Méthyl Succinate
146.14
(DMS)
ºá W=: »¸
Mono-Éthyl Succinate
146.14
(MES)
ºâ W=¸ »¸
Di-Éthyl Succinate
174.19
(DES)
ºâ W=¸ »¸
Mono-Butyle
174.19
Succinate (MBS)
º=? W?? »¸
Di-Butyle Succinate
230.30
(DBS)
º=? W?? »¸
Mono-Octyl Succinate
230.30
(MOS)
º?: W·â »¸
Di-Octyl Succinate
342.51
(DOS)
105
L’estérification de l’acide succinique avec un alcool est une succession de deux
réactions d’estérification : dans un premier temps, la formation du mono-ester et dans
un second temps la formation du di-ester. Aucune autre réaction secondaire n’est
considérée quels que soient les systèmes étudiés. Les seules réactions considérées sont
donc:
t7 + %»W ↔ W? » + E,4,<901+
Les cinétiques chimiques sont issues des travaux d’Orjuela & al, (2012) sur
l’estérification de l’acide succinique par éthanol. Les paramètres des lois cinétiques sont
considérés comme indépendants de l’alcool choisi. Une telle hypothèse est conséquente,
mais aucune étude sur l’estérification de l’acide succinique par différents alcools n’a été
trouvée dans la littérature. Pour compléter notre méthodologie de conception lors de la
dernière étape de validation et d’expérimentation, en plus d’un pilotage sur colonne
pilote, une étude de l’estérification de l’acide succinique avec différents alcools pourra
être menée.
Le Tableau III-2 donne les paramètres du modèle pseudo homogène exprimé par
rapport aux activités des constituants décrit par les équations suivantes.
106
<5= (FJf MJf )\F´¬ Y M´¬ Y ]
+= = Ð= 1F• m n ä(FHÊ MHÊ )(FŒ@Y´ MŒ@Y´ ) − æ
%z K ®å,=
Réaction 1 Réaction 2
Paramètres
Estérification Acide Succinique Estérification Mono-ester
èé (J.mol-1)
k (mol.s-1) 4,03E+07 5,46E+08
46900 61400
Keq 48.9 10.14
Tableau III-2: Constantes d’équilibre chimiques pour les deux réactions consécutives.
Pour faire ce choix, nous allons interpréter les diagrammes liquide-liquide réactifs
pour juger de la démixtion du système en fonction de la réactivité de l’alcool. Le
composant principal de la phase aqueuse est l'eau alors que pour la phase organique ce
sont les molécules non polaires ou légèrement polaires. Plus la masse molaire de l’alcool
sera importante plus le mono-ester et le di-ester seront non polaires. En conséquence, la
composition de la phase organique va également augmenter et la zone diphasique du
diagramme s’élargir. Lorsque le caractère polaire de l’alcool est assez important, comme
par exemple l’octanol, celui-ci forme aussi l’essentiel de la phase organique.
107
butanol et l’1-octanol. Comme indiqué précédemment, pour un alcool plus lourd, le 1-
octanol (C8), la zone diphasique sera plus grande que pour les alcools plus légers: 1-
butanol (C4) et éthanol (C2).
Figure III-1 : Diagramme liquide-liquide réactif pour différents alcools, éthanol, 1-butanol et 1-octanol.
108
III.1.1 Analyse de faisabilité
Figure III-2 : Diagramme liquide-liquide réactif pour le système Nc=5 Nr = 2. Acide Succinique, 1-Octanol, Eau,
mono-ester, di-ester représenté en compositions réactives.
109
Etant donné que ce système a une variance de 2, l'enveloppe de phase est une
courbe, et les courbes de raffinat réactif sont confondues avec les courbes d'équilibre du
diagramme liquide–liquide réactif. Par conséquent, les courbes de raffinat réactif ne
nous donnent pas plus d’informations que le diagramme liquide–liquide réactif pour le
système avec une variance inférieure ou égale à 2.
On peut observer tout de même que les courbes de raffinat réactif commencent
au point critique et se terminent au niveau des mélanges binaires eau-alcool
correspondant à chaque phase. Le point critique peut être considéré comme le nœud
instable et les mélanges binaires eau-alcool peuvent être considérés comme des nœuds
stables.
Pour notre étude de cas, l'objectif de séparation est fixé sur la composition du
raffinat récupéré en pied de colonne. On cherche à ce que celle-ci soit maximale voire
quasi pure en eau. Plus la teneur en eau dans le courant raffinat est élevée, plus le taux
de récupération en molécules C4 est élevé (C4 représente toutes les molécules à 4 atomes
de carbone, l'acide succinique et les deux esters). On détermine donc cette composition
directement avec le diagramme liquide-liquide réactif comme le point de composition le
plus proche du sommet eau. Le Tableau III-3 donne la composition déterminée.
Le choix que nous avons fait est d’avoir le meilleur rendement d’extraction en
molécules valorisables (molécules C4). C’est un choix pertinent si la matière première
(mélange fermentaire traité par extraction liquide-liquide réactive) représente le
principal coût pour le produit fini.
110
Enfin, il faut bien préciser qu’une étude détaillée du procédé nécessiterait
d’utiliser des outils d’analyses technico-économiques permettant de déterminer le cas le
plus favorable économiquementparlant.
L’analyse des régimes permanents permet de déterminer la valeur des limites grâce
à la spécification de séparation. La méthode est la suivante :
111
III.1.2 Synthèse et conception
112
Dans cette partie nous avons déterminé le design d’une colonne d’extraction
liquide-liquide réactive pour récupérer l’acide succinique d’un milieu aqueux. Le
système diphasique réactif dans cet exemple comporte 5 constituants et deux réactions.
III.1.3 Simulation
113
1
0,9
0,8
0,7
0,6
0,5
0,4
0,3
0,2
0,1
0
Taux de
Taux de
récupération en Composition en C4
conversion
tête E1
équilibre 0,923 1,000 0,658
cinétique 0,078 0,450 0,184
Figure III-5: Influence de la cinétique chimique sur les performances de la colonne d’extraction
Les cinétiques de réaction ont un fort impact sur les performances de la colonne, le
taux de conversion, le taux de récupération et la composition en tête de colonne sont
réduit respectivement de 92%, 55% et 72%. L’influence de la cinétique est aussi
identifiable sur l’analyse des profils de composition de la colonne,
1 1
0,9 0,9
0,8 0,8
0,7 0,7
Composition
Composition
0,6 0,6
0,5 0,5
0,4 0,4
0,3 0,3
0,2 0,2
0,1 0,1
0 0
0 1 2 3 4 5 0 2 4 6
Étage Étage
Profil de composition simulation à l’équilibre, phase Profil de composition simulation à l’équilibre, phase
aqueuse organique
114
1 1
0,9 0,9
0,8 0,8
0,7 0,7
Composition
Composition
0,6 0,6
0,5 0,5
0,4 0,4
0,3 0,3
0,2 0,2
0,1 0,1
0 0
0 1 2 3 4 5 0 2 4 6
Étage Étage
Profil de composition simulation cinétique, phase Profil de composition simulation cinétique, phase
aqueuse organique
Figure III-6: Influence de la cinétique sur les profils de composition dans la colonne d’extraction liquide-liquide
réactive
On remarque que les profils de composition sont plus plat ce qui traduit une
conversion de l’acide succinique plus lente et les performances de la colonne sont donc
limitées.
115
Le taux de solvant est donc de 0,15, le débit de solvant est pur en 1-octanol et la
rétention liquide est de 0,035 m3. Les figures III-7 et III-8 illustrent les effets du nombre
d’étages théoriques sur le taux de conversion de l’acide succinique, sur le taux de
récupération dans le courant extrait 1 en tête de colonne et sur la composition dans ce
même courant.
Composition
Composition molaire
Débit molaire
alimentation (%) Taux de Rétention
alimentation solvant (%)
solvant (-) liquide (m3)
Acide (mol.s-1)
Eau Octanol
succinique
20% 80% 2 0,15 100% 0,035
Tableau III-7 : Paramètres de simulation pour l’étude de l’influence du nombre d’étages théoriques
1 1
0,8 0,8
Critères (-)
Critères (-)
0,6 0,6
0,4
0,4
0,2
0,2
0
0 0 2 4 6
0 2 4 6
NET NET
Taux de conversion taux de recup en tete E1 Taux de conversion taux de recup en tete E1
compo C4 compo C4
Figure III-7 : Influence du nombre d’étages théoriques Figure III-8 : Influence du nombre d’étages théoriques
sur les critères choisis, simulation à l’équilibre sur les critères choisis, simulation hors équilibre
Ces résultats peuvent aussi être confrontés aux résultats de l’étape de conception.
En considérant l’équilibre thermodynamique sur chaque étage, quatre plateaux
théoriques sont nécessaires pour atteindre un taux de conversion de 92,30%, un taux de
récupération quasi-total de 99,99% et une composition de 65,77%. Or les simulations
réalisées montrent que seulement deux plateaux sont nécessaires. Les raisons de ce
phénomène peuvent être expliquées avec l’observation de l’annexe F sur l’orientation et
116
la sensibilité des conodales. L’orientation des conodales surtout dans la partie basse du
diagramme liquide-liquide réactif, peut conduire à des écarts de composition importants
qui ont pour conséquence de tronquer le nombre d’étages théoriques.
Composition molaire
Composition
alimentation (%) Rétention
NET (-) molaire
Acide liquide (m3)
Eau solvant (%)
succinique
20% 80% 4 100% 0,035
Tableau III-8 : Paramètres de simulation pour l’étude de l’influence du taux de solvant
1 1
0,8 0,8
Critères (-)
0,6 0,6
Critères
0,4 0,4
0,2 0,2
0 0
0 1 2 0 0,5 1 1,5 2
taux de solvant taux de solvant
Taux de conversion Taux de conversion
taux de recup en tete E1 taux de recup en tete E1
compo C4 compo C4
Figure III-9 : Influence du taux de solvant sur les Figure III-10 : Influence taux de solvant sur les critères
critères choisis, simulation à l’équilibre choisis, simulation hors équilibre
L’octanol étant à la fois agent extractant et réactif, plus le taux de solvant sera
important et meilleure seront la conversion et la récupération. En revanche, plus la
quantité d’octanol introduite dans la colonne sera importante, plus la composition sera
117
faible. La quantité d’acide succinique extraite sera la même mais celle-ci sera dans un
mélange avec une grande quantité d’alcool ceux qui peut avoir une conséquence sur la
suite du procédé, les quantités d’octanol à régénérer seront plus importante.
Composition molaire
alimentation (%) Taux de Rétention
NET (-)
Acide solvant (-) liquide (m3)
Eau
succinique
20% 80% 0,15 4 0,035
Tableau III-9 : Paramètres de simulation pour l’étude de l’influence de la pureté du solvant
118
1 1
0,8 0,8
Critères (-)
Critères (-)
0,6 0,6
0,4 0,4
0,2 0,2
0 0
0 0,2 0,4 0,6 0,8 1 0 0,2 0,4 0,6 0,8 1
Composition du solvant Composition du solvant
Taux de conversion Taux de conversion
taux de recup en tete E1 taux de recup en tete E1
compo C4 compo C4
Figure III-11 : Influence de la composition du solvant Figure III-12 : Influence de la composition du solvant
sur les critères choisis, simulation à l’équilibre sur les critères choisis, simulation hors équilibre
L’ajout d’eau, qui est un produit des réactions d’estérification, dans le débit de
solvant déplace l’équilibre vers la formation d’acide succinique et d’alcool. Ce qui
impacte fortement le taux de conversion car en ajoutant de l’eau on favorisera la
formation d’acide succinique. Le taux de conversion global de l’acide succinique (i.e. la
part d’acide succinique convertie en mono-ester ou di-ester) diminue. Si l’on étudie
l’influence de la composition sur le taux de conversion de l’acide succinique en ester, on
s’aperçoit que plus la composition est riche en eau, plus le taux de conversion de l’acide
succinique en di-ester diminue (baisse de 24% à 5%), alors que le taux de conversion en
mono-ester est relativement proche d’une valeur moyenne de 65%. Cette observation a
une importance dans la conception globale du procédé en prenant en compte la
régénération du solvant. Cette analyse permet de décider jusqu’où pousser la
séparation, la régénération du solvant suivant les objectifs de récupération de l’acide
fixés.
119
1
0,9
0,8
0,7
0,6
Critères (-)
0,5
0,4
0,3
0,2
0,1
0
0,0000001 0,00001 0,001 0,1 10 1000 100000
Rétention (m3)
120
HEPT (m) Volume étage (m3)
Min Max Min Max
colonne garnie 0,4 0,67 0,012 0,021
colonne pulsée 0,17 0,29 0,005 0,009
colonne plateaux perforées 0,83 1,25 0,026 0,039
colonne à plateaux pulsée 0,83 1,25 0,026 0,039
colonne Scheibel 0,11 0,2 0,003 0,006
Colonne RDC 0,29 0,4 0,009 0,012
colonne Khuni 0,13 0,2 0,004 0,006
colonne Karr 0,14 0,29 0,004 0,009
Contacteur RTL 0,08 0,17 0,002 0,005
Tableau III-10: Ordre de grandeur des rétentions liquides par rapport aux différentes technologies d’extraction
III.1.4 Conclusion
L’étape de vérification nous permet aussi d’affiner cette configuration comme pour
l’exemple choisi où la limitation cinétique est très contraignante. Ici les opérations
unitaires d’extraction liquide avec une forte capacité de rétention seront privilégiées.
L’option d’une cascade de mélangeurs décanteurs est aussi à prévoir lors de l’étape de
construction du procédé.
Nous avons vu aussi que suivant l’alcool choisi dans cette stratégie, la forme et la
taille du diagramme liquide-liquide réactif sont impactées. Le dimensionnement peut
l’être aussi et au lieu de prendre des alcools de moins en moins polaires, plus lourds et
donc qui réagit moins bien avec l’acide succinique, une autre stratégie d’extraction
s’offre à nous avec l’emploie d’un alcool léger pour la réaction d’estérification et d’un
solvant inerte pour la réaction pour extraire les esters. L’exploitation de la méthode de
dimensionnement sur un exemple à 6 constituants et 2 réactions est présentée dans la
partie 3.2 de ce chapitre.
121
III.2 Stratégie B : le solvant est inerte par rapport à la réaction
Le premier exemple présenté avait pour vocation d’illustrer la méthode proposée
avec une application sur un système avec cinq constituants et deux réactions.
L’utilisation des compositions réactives permet d’avoir une représentation des résultats
en deux dimensions et permet de comprendre la démarche de la méthode et d’identifier
les résultats que l’on peut attendre. Dans cette partie, la méthode est appliquée à un
système avec six constituants et deux réactions. Le schéma réactionnel est similaire,
l’estérification de l’acide succinique et du mono-ester avec un alcool léger, le méthanol
et un constituant inerte par rapport aux réactions qui jouent le rôle d’agent extractant.
Ici le solvant choisi est le para-xylène, un hydrocarbure aromatique immiscible avec
l’eau. Le choix de ce constituant est arbitraire mais une étude de l’influence de différents
hydrocarbures sur le diagramme liquide-liquide réactif permettrait de définir le solvant
qui permet d’obtenir la séparation la plus efficace.
Le système choisi pour cette partie possède six constituants et deux réactions. Le
méthanol, l'alcool choisi est miscible en toute proportion dans l'eau et les esters ne sont
probablement pas assez apolaires pour générer une phase apolaire en grande
proportion. L’alcool ne peut donc pas être dans cet exemple à la fois le réactif et l’agent
extractant. Il est nécessaire d'utiliser un constituant inerte par rapport aux réactions
comme solvant. La conception d'une colonne d'extraction réactive apparaît comme le
véritable défi pour ce système. Notre méthodologie de conception est appliquée de la
même manière, mais la représentation est plus complexe en raison de la forme de
pentaèdre du domaine de composition réactive.
122
Figure III-14 : Diagramme liquide-liquide du système acide succinique /méthanol/eau/mono-ester/di-
ester/para-xylène.
Plane III
Plane II
Plane IV
Plane I
123
Cette représentation de la succession de plans réactifs permet d'identifier une
courbe critique, formée avec tous les points critiques de chaque plan d’équilibre liquide-
liquide réactif. Cette ligne remarquable commence au point critique du système para-
xylène-Eau-Méthanol et se termine au point critique du système ternaire Acide
Succinique-Mono-ester-para-xylène (figures III-15-b et III-15-c)
Figure III-16 : Analyse de faisabilité, courbe de raffinat réactif du système Acide succinique
/méthanol/eau/mono-ester/di-ester/para-xylène.
Pour la phase aqueuse, les courbes vertes, toutes les Courbes De Raffinats
Réactifs commencent à un nœud instable représenté par les points critiques situés sur la
ligne rouge, et se terminent au niveau du nœud stable symbolisé par la phase aqueuse
du système alcool-eau-para-xylène réactif symbolisé par le point vert.
124
Le nœud principal est le point commun à toutes les courbes d'extraction et peut
être repéré près du sommet eau. Cette composition est thermodynamiquement possible
et parce qu'il est près du sommet de l'eau, il correspond à la composition la plus riche en
eau et donc la plus pauvre en molécules organiques. Du fait de cette caractéristique, ce
point pourrait être choisi comme une spécification de séparation étant donné que le flux
aqueux du procédé est faible en molécules valorisables. En effet, le taux de récupération
obtenu ici serait le meilleur. D’autres options sont possibles, comme choisir la
composition la plus riche en molécules C4 dans la phase organique, ce choix
correspondrait au système permettant de minimiser la consommation énergétique du
procédé de récupération.
125
Composition molaire
Débit
Courant Acide Mon-octyl Di-octyl Para-
(mol.s-1) Méthanol Eau
Succinique succinate succintae Xylène
F 2 0.156 0.220 0.624 0 0 0
S 4 0 0 0 0 0 1
RN 1.67 0.00564 0.00604 0.988 0.000467 0 0
E1 4.33 0.00210 0.00111 0.00199 0.0138 0.0133 0.968
NTS=5
Tableau III-12 : Résultats de l’analyse de régime permanent de l’analyse de faisabilité sur l’exemple de colonne
étudiée avec Nc=6 et Nr=2.
Figure III-17-b : Plan I de représentation de l’analyse Figure III-17-c : Plan II de représentation de l’analyse
de faisabilité de faisabilité
126
coplanaires, tandis que les points représentant les compositions de l'extrait et du
raffinat sont de part et d'autre du plan. Pour la figure III-17-c, c’est le contraire, les
points représentatifs des compositions de l'extrait et du raffinat sont coplanaires et sur
l'enveloppe de phase tandis que les points S et F sont de part et d'autre du plan vert.
Sur la Figure III-18, la succession des lignes d'équilibre et des droites opératoires
sont représentées sur la représentation en trois dimensions. Avec cette représentation,
il est difficile de déterminer le nombre d'étages théoriques. On utilise alors
préférentiellement les représentations des figures III-19 et III-20.
127
Figure III-18 : Détermination du nombre d’étages théoriques pour le système avec Nc=6 et Nr=2
0,1
3 2 1
0,01
4
0,001
SA SA-E1
Figure III-19 : Profil de composition réactive en acide succinique, détermination du NET, NET=5
128
0,04
0,025
succinique
3
0,02
0,015
4
0,01
5
0,005
F
S 0
0R 0,05 0,1 0,15 0,2
N
Composition réactive aqueuse en acide succinique
courbe opératoire Courbe d'équilibre Détermination of NET
Figure III-20 : Diagramme de distribution pour la détermination du NET d’une façon analogue à McCabe
Thièle
129
III.2.3 Simulation
• le taux de conversion,
• le taux de récupération de molécules valorisables dans le courant de tête de
colonne
• et la composition en constituants valorisables.
Les paramètres étudiés sont :
Compositions molaire
Taux de Rétention
alimentation du solvant (%) NET (-)
solvant (-) liquide (m3)
Para-Xylène
100% 2 5 0,035
Tableau III-13 : Paramètres de simulation
Compositions molaires
Débit
AS MeOH Eau m-ester d-ester P-X
(mol.s-1)
F 0,156 0,22 0,624 0 0 0 2
Tableau III-14 : Composition et débit de la charge traitée par extraction liquide réactive
130
sur les performances de la colonne d’extraction liquide-liquide réactive en termes de
taux de conversion, de taux de récupération et de composition. Les résultats et la
comparaison de ces deux simulations sont présentés sur la Figure III-21.
1
0,9
0,8
0,7
0,6
0,5
0,4
0,3
0,2
0,1
0
Taux de récupération
Taux de conversion Composition en C4
en tête E1
équilibre 0,883 0,802 0,058
cinétique 0,016 0,026 0,002
Figure III-21: Influence de la cinétique chimique sur les performances de la colonne d’extraction
On remarque que les performances sont toutes réduites de plus de 95%. Les
cinétiques de réactions ont donc un fort impact sur ce système avec un alcool léger
comme réactif des réactions d’estérifications et un solvant non polaire.
131
lorsque l’on a un faible nombre d’étages théoriques. Cette forte influence du NET est en
fait liée au temps de séjour dans la colonne : chaque étage théorique a un volume fixé ;
en augmentant le NET, on augmente donc le volume de la colonne, donc le temps de
séjour, donc la conversion lorsque l’on est limité cinétiquement.
On constate de plus que les valeurs limites du taux de conversion et du taux de
récupération sont beaucoup plus faibles ; íÎî = 1,7% et ¿Îî« = 2,6%. Ceci s’explique par
f
1 0,03
0,8 0,025
Critères (-)
0,02
Critères (-)
0,6
0,015
0,4
0,01
0,2 0,005
0 0
1 2 3 4 5 6 7 8 9 0 5 10
Nombre d'étage Nombre d'étage
Taux de conversion Taux de conversion
taux de récupération en tête E1 taux de récupération en tête E1
composition C4 composition C4
Figure III-22 : Influence du nombre d’étages théoriques Figure III-23 : Influence du nombre d’étages théoriques
sur les critères choisis, simulation à l’équilibre sur les critères choisis, simulation hors équilibre
132
d’estérification, ce qui aura tendance à déplacer l’équilibre chimique dans le sens
indirect et donc aura un impact négatif sur le taux de conversion. L’augmentation de la
quantité de solvant dans la colonne a aussi pour effet de diluer le courant d’extrait donc
la composition en molécules C4 en tête de colonne diminue pour atteindre une
composition faible de `ÎIIî ,= = 5% pour des taux de solvant de l’ordre de S/F=4.
1 0,035
0,03
0,8
0,025
Critères (-)
Critères (-)
0,6 0,02
0,015
0,4
0,01
0,2 0,005
0 0
0 2 4 6 0 2 4 6
Taux de solvant Taux de solvant
Taux de conversion Taux de conversion
taux de recup en tete E1 taux de recup en tete E1
compo C4 compo C4
Figure III-24 : Influence de la composition du solvant Figure III-25 : Influence de la composition du solvant
sur les critères choisis, simulation à l’équilibre sur les critères choisis, simulation hors équilibre
En prenant en compte les effets des phénomènes de cinétique chimique (Figure III-25),
les conclusions sur le taux de récupération et la composition en tête de colonne sont les
mêmes. En revanche le taux de conversion baisse lorsque le taux de solvant croît : le
taux de conversion atteint pour un taux de solvant élevé (S/F=4) est de τõ«î = 2,87%
ö
alors que pour le taux de solvant S/F=2, le taux de conversion est de τõ«î = 2,91%. La
ö
variation est infime mais ceci traduit encore une limitation cinétique. En effet, une
augmentation du débit de solvant entraine une baisse du temps de séjour dans la
colonne, et donc une baisse de conversion lorsque l’on est fortement limité en transfert.
On précise aussi que la réaction est considérée dans la phase aqueuse et que donc le
débit de la phase organique aura moins d’influence sur le taux de conversion que si la
réaction a lieu en phase organique.
133
donc être étudiée sur les trois critères d’efficacité de la colonne. Les figures III-26 et III-
27 présentent l’influence de la composition du solvant sur le taux de conversion, le taux
de conversion en tête de colonne et la composition.
1 0,03
0,8 0,025
Critères (-)
Critères (-)
0,02
0,6
0,015
0,4
0,01
0,2
0,005
0 0
0,5 0,6 0,7 0,8 0,9 1 0,5 0,6 0,7 0,8 0,9 1
Composition du solvant Composition du solvant
Taux de conversion Taux de conversion
taux de recup en tete E1 taux de recup en tete E1
compo C4 compo C4
Figure III-26 : Influence de la composition du solvant Figure III-27 : Influence de la composition du solvant
sur les critères choisis, simulation à l’équilibre sur les critères choisis, simulation hors équilibre
134
1
0,8
Critères (-)
0,6
0,4
0,2
0
0,0000001 0,00001 0,001 0,1 10 1000 100000
Rétention (m3)
III.2.4 Conclusion
L’étape de vérification nous permet aussi d’affiner cette configuration comme pour
l’exemple choisi où la limitation cinétique est très contraignante. Ici les opérations
unitaires d’extraction liquide-liquide avec une forte capacité de rétention seront
privilégiées. L’option d’une cascade de mélangeurs décanteurs est aussi à prévoir lors de
l’étape de construction du procédé.
135
− La distillation simple, pour séparer l’eau de l’acide succinique grâce à l’apport
d’énergie thermique et de la différence de volatilité entre les deux espèces.
− Une extraction liquide-liquide simple avec régénération du solvant. La propriété
de séparation mise en jeu ici est la différence de solubilité. Le soluté du système sera
plus soluble dans le solvant que dans le diluant et un transfert de matières s’opérera
entre les phases aqueuses et organiques pour séparer le soluté et le récupérer après
régénération du solvant.
− Le dernier procédé est l’extraction liquide-liquide réactive avec régénération du
solvant. Le soluté va réagir dans une réaction d’estérification pour modifier ses
propriétés de solubilité et améliorer sa séparation. La récupération se fera après
régénération du solvant. Plusieurs options sont possibles comme nous l’avons vu, le
solvant peut être aussi le réactif de la réaction mise en jeu, ou le solvant peut être inerte
et le réactif est un constituant supplémentaire ajouté.
137
Ces résultats sont les paramètres d’entrée d’une colonne de distillation à condenseur
total que nous allons simuler maintenant. Pour le calcul, le logiciel fait une série
d’hypothèses. La colonne est alors considérée comme un ensemble d’étages théoriques où on
suppose :
o un régime permanent
o il n’y a pas de réaction
o un système bi-phasique (liquide/vapeur).
o Le système est à l’équilibre thermodynamique
Alimentation Alimentation
Distillat Résidu
25°C Liq-bouillant
DEBITS
SUCCINIC ACID 170,047 170,047 1,700 168,346
MASSIQUES
PARTIELS
WATER 103,768 103,768 103,663 0,105
(KG/HR)
FRACTIONS SUCCINIC ACID 0,2 0,2 0,002 0,996
MOLAIRES WATER 0,8 0,8 0,998 0,004
FRACTIONS SUCCINIC ACID 0,621 0,621 0,016 0,999
MASSIQUES WATER 0,379 0,379 0,984 0,001
DEBIT MOLAIRE TOTAL
7,2 7,2 5,76856 1,43144
(KMOL/HR)
DEBIT MASSIQUE TOTAL
273,815 273,815 105,363 168,452
(KG/HR)
TEMPERATURE (K) 273.15 500 373,237 571,242
PRESSION (ATM) 1 1 1
Tableau III-18 : Bilan matières du procédé de distillation pour la récupération de l’acide succinique
138
Ces résultats seront donc repris à la fin de ce paragraphe pour la comparaison des
différents procédés étudiés.
5
Décanteur
4 8
Colonne à distiller 7
1
Eau
Charge
10
6
9
Sortie Production Molécules C4
2 3 Mélangeur
Raffinat
Figure III-29 : Schéma de simulation du procédé d’extraction pour la récupération d’acide succinique
La simulation de ce procédé est réalisée sur ProSimPlus 3®. Une première étape
de détermination du pré-design des colonnes d’extraction et de distillation a été menée.
Le taux de solvant et le nombre d’étages théoriques ont été déterminés respectivement
et la simulation shortcut a été réalisée. Les résultats de cette étape se retrouvent dans
les tableaux suivants et on retrouve en annexe H le tableau des courants pour ce
procédé. Après optimisation du procédé, les différents résultats de simulation sont
regroupés dans les tableaux suivants.
Colonne d’extraction
Plateau théorique 10
Taux de solvant 1,04
Colonne de distillation
Taux de reflux R 0,016
16 11 plateaux d’enrichissement
(condenseur compris)
Plateau théorique
(bouilleur et condenseur compris)
5 plateaux d’appauvrissement
(bouilleur compris)
Chaleur récupérable au
172 kW
condenseur
Chaleur à fournir au bouilleur 245 kW
Décanteur
Taux de partage 0.842
Tableau III-20 : Paramètres des différentes opérations unitaires du procédé d’extraction liquide
140
III.3.3 Extraction liquide-liquide réactive avec régénération du solvant
Comme il a été vu dans les parties 3.1 et 3.2 de ce chapitre, deux stratégies
d’extraction réactive sont possibles. La première est de considérer le solvant comme
réactif de la réaction et l’autre stratégie est d’employer un solvant inerte par rapport à la
réaction et qui fait seulement office d’agent extractant. Les procédés pour ces deux
stratégies ont été étudiés, le procédé A avec l’octanol comme solvant et le procédé B
avec un hydrocarbure lourd comme solvant, le para-xylène, et le méthanol comme
réactifs de la réaction.
141
Taux Composition Débit Débit Consommation Energétique Du
Récupération Molaire Total Partiel Procédé
- - kg/h kg/h kW kWk/kg de produit
99% 0,996 358,283 357,538 83,35 0,233
Tableau III-23 : Résultats de simulation du procédé d’extraction liquide réactive pour la récupération de
l’acide succinique
Ces résultats seront repris dans la conclusion de cette partie pour la comparaison
des différents procédés. Une nuance doit néanmoins être apportée à ces résultats, le
produit prix en compte ici n’est pas le même, c’est une mélange d’acide succinique de
mono-ester et de di-ester. La comparaison est faite sans prendre en compte l’énergie
nécessaire à hydrolyser les ester pour retrouver l’acide succinique mais celle–ci devra
être prise en compte pour une étude plus poussée.
142
Charge 1 2 Alimentation Alcool
7 8
Appoint Para-Xylène
Mélangeur 3 6
Mélangeur
9
Colonne à distiller
Colonne d'extraction réactive
Production
Raffinat
Figure III-30 : Schéma de simulation du procédé d’extraction réactive avec 6 constituants pour la
récupération d’acide succinique
La simulation de procédé est réalisée sur ProSimPlus 3®. Dans un premier temps
une simulation shortcut a permis de déterminer les paramètres structuraux de la
colonne de distillation. Dans un second temps, l’optimisation du procédé a été réalisée
pour atteindre l’objectif de séparation qui est un courant de production avec une
composition de 99,6%. Comme le procédé de séparation est réactif, la composition prise
en compte ici pour l’objectif de production est la composition réactive de l’acide
succinique. Celle-ci se calcule comme la somme de la composition de l’acide succinique
et des deux esters. Les tableaux ci-dessous regroupent les configurations des colonnes
du procédé ainsi que les résultats de simulation et on retrouve ne annexe H le tableau
des courants du procédé.
143
Colonne d’extraction
Plateau théorique 5
Taux de solvant 2
Colonne de distillation
Taux de reflux R 0,663
26 20 plateaux d’enrichissement
(condenseur compris)
Plateau théorique
(bouilleur et condenseur compris)
6 plateaux d’appauvrissement
(bouilleur compris)
Chaleur récupérable au
328,02 kW
condenseur
Chaleur à fournir au bouilleur 384,93 kW
Tableau III-24 : Paramètres des différentes opérations unitaires du procédé d’extraction liquide réactive
On remarque ici aussi que le produit pris en compte n’est pas de l’acide
succinique pur mais un mélange d’acide succinique et d’esters. Quelques remarques
supplémentaires sur cette stratégie de récupération de l’acide succinique avec le para-
xylène comme solvant. L’objectif de séparation fixé qui est d’avoir un courant de
production avec une composition en molécules C4 de 99,6 % nécessite un grand apport
d’énergie. Ceci est dû à la température d’ébullition du solvant qui est inférieure aux
températures d’ébullition de l’acide succinique et des esters comme le montre le tableau
ci-dessous.
Température d’ébullition en °C
Méthanol 64,7
Eau 100
Para-Xylène 137,85
di-méthyl succinate 196,2
Mono-méthyl succinate 223,06
Acide succinique 235
Tableau III-26 : Température d’ébullition des constituants du système Nc=6
144
Ces résultats seront repris dans la conclusion de cette partie pour la comparaison
des différents procédés.
− Le débit partiel de molécules valorisables. Pour les procédés non réactifs est pris
en compte le débit partiel d’acide succinique et pour les procédés réactifs la
somme des débits partiels d’acide succinique et des esters.
− La consommation énergétique au travers de deux valeurs, la quantité de chaleur
à apporter au bouilleur des colonnes du procédé, et cette même quantité de
chaleur ramenée à la production d’un kilogramme heure de produit.
Les débits de production exprimés ici en kmol de molécule C4 par heure sont tous
du même ordre sauf pour le dernier procédé. La cause est le taux de récupération de
molécules C4 de la colonne d’extraction réactive qui n’est que de 80%.
1,600
1,400
1,200
1,000
0,800
0,600
0,400
0,200
0,000
Extraction Extraction
Extraction liquide liquide
Distillation
liquide réactive réactive
Nc=5 Nc=6
débit partiel kmol/h 1,426 1,440 1,438 1,130
Figure III-31 : Débit partiel du courant de production de molécule C4
145
réactive à 6 constituants, nécessite le plus d’énergie. La raison est le choix du solvant,
celui-ci ayant une température d’ébullition non adaptée au système. Comme nous avons
pu le démontrer auparavant, le choix du solvant est donc primordial non seulement
pour l’unité d’extraction réactive mais aussi pour l’unité de régénération du solvant. Des
tests ont été réalisés avec d’autres solvants de type n-alcanes long (n-hexadecane), qui
ont la particularité d’avoir un point d’ébullition élevé (286°C). Nous pouvons donc peut
penser qu’ils sont davantage adaptés pour l’unité de régénération, mais ceux-ci forment
des mélanges azéotropiques qui complexifient le schéma de séparation.
400,00
350,00
300,00
250,00
200,00
150,00
100,00
50,00
0,00
Extraction Extraction
Extraction
Distillation liquide liquide
liquide
réactive Nc=5 réactive Nc=6
kW 100,31 245,00 83,35 384,93
Figure III-32 : Consommation énergétique globale des procédés étudiés
2,500
2,000
1,500
1,000
0,500
0,000
Extraction Extraction
Extraction liquide liquide
Distillation
liquide réactive réactive
Nc=5 Nc=6
kWh/kg de produit 0,595 1,44 0,23 2,39
Figure III-33 : Consommation énergétique par kg/h de molécule produit des procédés étudiés
146
La quantité d’énergie la plus faible nécessaire est pour le procédé d’extraction
réactive avec 5 constituants. Cette première analyse montre que l’extraction réactive
peut être un procédé viable pour la récupération de l’acide succinique. Cependant
comme nous l’avons précisé au début de ce paragraphe, cette analyse économique
nécessite d’être plus poussée pour affirmer de façon certaine que l’extraction liquide
réactive est réellement un procédé viable pour la récupération de l’acide succinique.
147
CONCLUSION GES NES RALE
Les bioprocédés et les couplages fonctionnels sont deux axes de recherche
prédominants en génie des procédés et dans la mouvance actuelle axée sur le
développement durable. Les bioprocédés ont l’avantage certain de pouvoir substituer
aux matières premières fossiles des biomasses en mettant en œuvre des micro-
organismes pour produire des molécules d’intérêt. Cependant, avec la fermentation qui
est l’étape de transformation des bioprocédés, les milieux fermentaires à traiter pour
récupérer les molécules d’intérêt sont complexes et nécessitent la plupart du temps des
procédés de séparation très énergivores et très coûteux. Le recours au couplage
fonctionnel des opérations de réaction et de séparation est une solution intéressante
puisqu’un tel agencement permet une intégration thermique et une capacité des
équipements.
La première étape, l’analyse de faisabilité repose sur trois points essentiels pour
le bon déroulement de la méthode de conception et pour une bonne représentation des
résultats obtenus. Tout d’abord, la construction du diagramme liquide-liquide réactif qui
est l’outil de base pour la compréhension et la représentation des résultats. Il présente
149
les zones monophasiques et diphasiques du système supposé à l’équilibre chimique
délimité par la courbe de démixtion et est représenté avec le système de compositions
réactives. Le jeu de compositions réactives permet de réduire la dimension de
représentation du système. Un système avec Nc constituants représenté en Nc-1
dimensions, et Nr réactions verra sa dimension de représentation réduite à Nc-Nr-1 en
choisissant Nr constituants de référence. Ensuite les courbes de raffinats réactifs sont un
outil développé sur le même principe des courbes de résidus réactifs en distillation
réactive. Cet outil est une aide au concepteur pour déterminer un objectif de séparation.
Nous avons vu aussi que cet outil n’avait un rôle intéressant que pour les systèmes où
Nc-Nr>3, en effet les courbes de raffinats réactifs sont la succession d’états d’équilibres
liquide-liquide réactifs et se trouvent donc sur la courbe ou la surface de démixtion.
Dans le cas où Nc-Nr≤3, les courbes de raffinat réactif redéfinissent la courbe de
démixtion en 2D, dans le cas où Nc-Nr>3, les courbes de raffinat réactif définissent des
chemins de composition sur la surface de démixtion en trois dimensions. Enfin le
dernier point de cette étape est l’analyse des régimes permanents qui permet de
sélectionner les proportions de réactif et de solvant nécessaires à la réalisation de
l’objectif de séparation.
150
de la cinétique a été étudiée et nos résultats n’ont pu être validés de façon
expérimentale.
Enfin, au niveau de la validation, au cours de notre étude nous avons été souvent
confrontés à un manque de données expérimentales qui nous a empêché de valider les
résultats obtenus à partir de la méthode de conception. Il est sûr que de telles données
doivent être acquises et seront très utiles dans cette étape de validation notamment sur
les cinétiques chimiques en milieu diphasique et sur les équilibres liquide-liquide
réactifs. Des essais sur une colonne pilote à l’échelle laboratoire sont aussi envisageables
afin de valider nos modèles et démontrer que ceux-ci sont conformes à la réalité.
151
PERSPECTIVES
A l’issue de ce travail, de nombreuses perspectives sont envisageables tant d’un
point de vue modélisation, simulation qu’expérimentation :
152
D’un point de vue expérimental, les perspectives et les voies de recherche sont
aussi nombreuses.
Une réflexion avec une approche technologique pourra être menée aussi sur la
mise en place d’une colonne d’extraction liquide-liquide réactive. Cette réflexion portera
153
essentiellement sur la catalyse de la réaction d’estérification. Comment sera amené le
catalyseur, sous sa forme homogène ou hétérogène, fixé à un interne de colonne ou libre
dans la colonne. Toutes ces questions devront trouver une réponse pour pouvoir mettre
en place physiquement le procédé d’extraction liquide-liquide réactif à l’échelle pilote
dans l’optique de valider expérimentalement les résultats de notre méthode de
conception ou à l’échelle industrielle pour réaliser une production d’acide succinique.
* *
*
154
ANNEXES
Annexe A : Données thermodynamiques du système acide acétique,
butanol, acétate de butyle et eau.
Le système est à température ambiante et pression ambiante, 293.15K et 1atm. Le
butanol est à la fois le réactif de la réaction d’estérification et l’agent de séparation. Les données
thermodynamiques de ce système sont reprises des travaux d’Olán-Acosta et al (2014) et
décrites en détail dans le tableau ci-dessous.
155
Annexe B : Établissement des équations des courbes de raffinat réactif
en fonction des compositions réactives
Soit l’expression du système liquide-liquide réactif de rétention H, alimenté par un courant F et
avec un soutirage E (Figure II-7). Ce système est le siège d’une ou plusieurs réactions.
L’expression du bilan matière global et du bilan matière partiel s’écrit :
-W
= ; − < + Sˆ U
-0
Equation 1
-(W ∙ ÀA )
= ; ∙ ÂA − < ∙ OA + SA U
-0
Equation 2
-W
= ; − < + Sˆ U
-0
Equation 3
-(W ∙ ÀA )
= ; ∙ ÂA − < ∙ OA + SA U
-0
Equation 4
-(W ∙ Ài )
= ; ∙ Âi − < ∙ Oi + •i U
-0
Equation 5
-(W ∙ Ài )
U = •i@= m − ; ∙ Âi + < ∙ Oi n
-0
Equation 6
On peut substituer cette expression de l’avancement dans les équations 3 et 4 et obtenir ainsi :
-W -Ài -W
= ; − < + Sˆ •i@= mW + Ài − ; ∙ Âi + < ∙ Oi n
-0 -0 -0
Equation 7
-W -(Ài )
\1 − Sˆ •i@= Ài ] = ;\1 − Sˆ •i@= Âi ] − <(1 − Sˆ •i@= Oi ) + Sˆ •i@= W
-0 -0
Equation 8
et
-ÀA -W -Ài -W
W + ÀA = ; ∙ ÂA − < ∙ OA + SA •i@= mW + Ài − ; ∙ Âi + < ∙ Oi n
-0 -0 -0 -0
Equation 9
-(ÀA − SA •i@= Ài ) -W
W + \ÀA − SA •i@= Ài ] = ;\ÂA − SA •i@= Âi ] − (< ∙ OA − SA •i@= Oi )
-0 -0
Equation 10
On trouve alors :
156
-\ÀA − SA •i@= Ài ]
W
-0
\ÀA − SA •i@= Ài ]
+ ‚;\1 − Sˆ •i@= Âi ] − <\1 − Sˆ •i@= Oi ]
(1 − Sˆ •i@= Ài )
Equation 13
-(Ài )
+ Sˆ •i@= W ƒ = ;\ÂA − SA •i@= Âi ] − <\OA − SA •i@= Oi ]
-0
ou,
\ÀA − SA •i@= Ài ]
+ \;^ − <^ ]
(1 − Sˆ •i@= Ài )
Equation 14
ou,
et,
\ÀA − SA •i@= Ài ]
-| \1 − Sˆ •i@= Ài ]•
û (1 − Sˆ •i@= Ài ) \ÀA − SA •i@= Ài ] -\1 − Sˆ •i@= Ài ]þ
Wú − ý
-0 (1 − Sˆ •i@= Ài ) -0
ù ü
\ÂA − SA •i@= Âi ] \ÀA − SA •i@= Ài ]
Equation 16
= ;^ ‚ − ƒ
(1 − Sˆ •i@= Âi ) (1 − Sˆ •i@= Ài )
\OA − SA •i@= Oi ] \ÀA − SA •i@= Ài ]
− <^ ‚ − ƒ
(1 − Sˆ •i@= Oi ) (1 − Sˆ •i@= Ài )
- A
W\1 − Sˆ •i@= Ài ] = ;^ (•A − A) − <^ (kA − A)
-0
Equation 17
soit,
157
- A
= ;^ (•A − A) − <^ (kA − A)
-¿
-0
-¿ =
Equation 18
W(1 − Sˆ •i@= Ài )
-W -(Ài )
\1 − Sˆ •i@= Ài ] = ;\1 − Sˆ •i@= Âi ] − <(1 − Sˆ •i@= Oi ) + Sˆ •i@= W
-0 -0
Equation 19
-W -(Ài )
\1 − Sˆ •i@= Ài ] − Sˆ •i@= W = ;\1 − Sˆ •i@= Âi ] − <(1 − Sˆ •i@= Oi )
-0 -0
Equation 20
-W -(1 − Sˆ •i@= Ài )
\1 − Sˆ •i@= Ài ] − W = ;\1 − Sˆ •i@= Âi ] − <(1 − Sˆ •i@= Oi )
-0 -0
Equation 21
-W(1 − Sˆ •i@= Ài )
= ;\1 − Sˆ •i@= Âi ] − <(1 − Sˆ •i@= Oi )
-0
Equation 22
On trouve :
-W
= ;^ − <^
-0
Equation 26
Or comme on travaille à rétention liquide constante la relation entre ;^ et <^ est établie comme
pour le cas non réactif. Ainsi l’équation 18 se réécrit :
-W
= 0 ↔ ;^ = <^
-0
Equation 27
- A
= •A − kA
-¿
;^ -0
-¿ =
Equation 28
158
Annexe C : Calcul des équilibres chimiques entres phases, tracé des
diagrammes liquide-liquide réactif.
Le diagramme fonctionnel du programme permettant de calculer et de tracer le
diagramme liquide-liquide réactif est présenté et brièvement détaillé ici.
Données thermodynamiques
Maillage du domaine de
composition
Choix de Nc composition
Pour utiliser ce programme l’utilisateur doit fournir les données thermodynamiques sur
Simulis qui est appelé à partir du programme Matlab principal.
Les calculs des équilibres liquide-liquide réactif sont réalisés pour toutes les
compositions déterminées au maillage du domaine de composition.
159
Annexe D : Calcul des courbes de raffinats réactifs
Le diagramme fonctionnel du programme de calcul des courbes de raffinats réactifs est
présenté et brièvement détaillé ici.
Données thermodynamiques
Pour utiliser ce programme, l’utilisateur doit fournir les données thermodynamique sur
Simulis. L’utilisateur choisit ensuite une composition ou une série de composition de mélange
pour tracer respectivement sur le diagramme liquide-liquide réactif une courbe de raffinats
réactifs ou un réseau de courbes de raffinats réactifs.
160
Annexe E : Analyses des régimes permanents
Les données brutes issues du diagramme liquide-liquide réactif et du réseau de courbes
de raffinats réactifs ainsi que les données relatives à l’alimentation du procédé d’extraction
réactive sont les suivantes :
Données thermodynamiques
161
Annexe F : Dimensionnement à partir des valeurs aux frontières.
Les données issues de l’analyse de faisabilité et exploitées dans ce programme sont :
− Les compositions et les débits de tous les courants entrant et sortant du procédé
d’extraction réactive.
− Les compositions à l’équilibre liquide-liquide réactif avec l’extrait et le raffinat
Données thermodynamiques
Profil
satisfaisant
Oui
Design de la colonne
Les calculs réalisés avec ce programme permettent de visualiser et d’analyser les profils
de compositions dans la colonne. Les paramètres de la colonne comme les compositions du
raffinat et de l’extrait ou le taux de solvant peuvent être modifié jusqu’à obtenir des profils
satisfaisant et un design de colonne.
162
Annexe G : Vérification des outils de simulation
Il s’agit ici de présenter l’exactitude des données calculées avec l’outil à l’équilibre en les
comparants aux données issues de la méthode de design. On vérifiera aussi l’outil de
simulation avec cinétique en imposant un temps de séjour très élevé et ainsi retrouver
les données calculées avec l’outil à l’équilibre
Pour vérifier les outils de simulation, nous allons nous appuyer sur les résultats de
l’étape de synthèse et de conception. Les figures 5 et 6 présentent les profils de
compositions au sien de la colonne, l’étage 0 représente l’étage d’alimentation de la
charge à traiter et l’étage 5 l’étage d’alimentation en solvant.
1 1
0,8 0,8
Compositions (-)
Compositions (-)
0,6 0,6
0,4 0,4
0,2 0,2
0 0
0 2 4 6 0 2 4 6
Étage Étage
Acide Succinque Octanol Eau Acide Succinque Octanol Eau
M-ester Di-ester M-ester Di-ester
Figure 1 : Profils de composition du courant aqueux Figure 2 : Profils de composition du courant organique
dans la colonne d’extration réative avec 5 constituants dans la colonne d’extration réative avec 5 constituants
calculés à l’étape de prédesign calculés à l’étape de prédesign
163
1 1
0,8 0,8
Composition (-)
Compositions (-)
0,6 0,6
0,4 0,4
0,2 0,2
0
0
0 2 4 6
0 2 4 6
Étage
Étage
Acide Succinque Octanol Eau Acide Succinque Octanol Eau
M-ester Di-ester M-ester Di-ester
Figure 3 : Profils de composition du courant aqueux Figure 4 : Profils de composition du courant organique
dans la colonne d’extration réative avec 5 constituants dans la colonne d’extration réative avec 5 constituants
calculés par simulation calculés par simulation
A première vue, les profils de compositions ne sont pas identiques mais cela ne traduit
pas forcement une erreur de calculs mais plutôt une erreur de sensibilité sur le calcul
des équilibres liquide-liquide réactifs. L’orientation des conodales fait qu’une petite
variation de composition de la phase aqueuse entraîne une forte variation de
composition de la phase organique.
1 1
Compositions (-)
0,8 0,8
0,6 0,6
Compositions (-)
0,4 0,4
0,2 0,2
0 0
0 1 2 3 0 1 2 3 4
Étage Étage
Enfin pour valider notre code de calcul cinétique, nous allons comparer les résultats de
deux simulations, l’une lancée avec l’outil à l’équilibre et une seconde avec l’outil
164
cinétique auquel on impose un temps de séjour infini, ce qui permet de retirer toute
sensibilité à la cinétique. Ces profils sont présentés sur les figures 11 à 14. La
comparaison des profils de composition montre une parfaite analogie entre les profils
calculés avec les deux modèles différent. Les modèles à l’équilibre et hors équilibre
chimique sont cohérent entre eux pour le cas limite d’un temps de séjour infinie et
justifie sont utilisations.
1 1
Compositions (-)
Compositions (-)
0,8 0,8
0,6 0,6
0,4 0,4
0,2 0,2
0 0
0 2 4 6 0 2 4 6
Étage Étage
SA Oc-OH W SA Oc-OH W
m-ester d-ester m-ester d-ester
Figure 7 ; Profils de composition du courant aqueux Figure 8: Profils de composition du courant organique
dans la colonne d’extration réative avec 5 constituants dans la colonne d’extration réative avec 5 constituants
calculés par simulation à l’équilibre calculés par simulation à l’équilibre
1 1
Compositions (-)
Compositions (-)
0,8 0,8
0,6 0,6
0,4 0,4
0,2 0,2
0 0
0 2 4 6 0 2 4 6
Étage Étage
SA Oc-OH W SA Oc-OH W
m-ester d-ester m-ester d-ester
Figure 9: Profils de composition du courant aqueux Figure 10 : Profils de composition du courant
dans la colonne d’extration réative avec 5 constituants organique dans la colonne d’extration réative avec 5
calculés par simulation hors équilibre constituants calculés par simulation hors équilibre
165
Annexe G : tableaux des courants des procédés de récupération de l’acide succinique
166
Procédé d’extraction liquide-liquide réactive avec régénération du solvant, procédé A.
Solvant Solvant
Charge Raffinat Solvant Extrait Distillat Production Eau
appoint recyclé
1 2 3 4 5 6 7 9 10
SUCCINIC ACID 170,047 0 0 13,252 0,000 13,252 0 0 0
DEBITS 1-OCTANOL 0 0,21 236,878 19,174 18,430 0,744 0,036 218,484 18,394
MASSIQUES WATER 103,768 122,99 0,954 10,847 10,847 0,000 9,893 0 0,954
PARTIELS MONO-OOCTYL
(KG/HR) 0 0 0,424 225,917 0,425 225,492 0 0 0,424
SUCC
DI-OCTYL SUCCINA 0 0 0 118,794 0 118,794 0 0 0
SUCCINIC ACID 0,621 0 0 0,034 0 0,037 0 0 0
1-OCTANOL 0 0,0017 0,994 0,049 0,621 0,002 0,004 1 0,930
FRACTIONS WATER 0,379 0,9983 0,004 0,028 0,365 0 0,996 0 0,048
MASSIQUES MONO-OOCTYL 0,014293
0 0 0,002 0,582 0,62937 0 0 0,021
SUCC 2
DI-OCTYL SUCCINA 0 0 0 0,306 0 0,331565 0 0 0
DEBIT MASSIQUE TOTAL (KG/HR) 273,815 123,195 238,256 387,985 29,702 358,283 9,930 218,484 19,772
TEMPERATURE (K) 298,15 298,15 298,15 298,15 367,56 585,13 367,56 298,15 367,56
PRESSION (ATM) 1 1 1 1 1 1 1 1 1
167
Procédé d’extraction liquide-liquide réactive avec régénération du solvant, procédé B.
Appoint
Alimentation Alimentation
Charge alcool colonne Raffinat Solvant Extrait Distillat Para- Production
Xylène
1 2 3 4 5 6 7 8 9
SUCCINIC ACID 132,627 0 132,627 26,080 0,000 9,242 0 0 41,744
DEBITS METHANOL 0 50,751 50,751 0,116 2,746 1,323 2,746 0 0
MASSIQUES WATER 80,960 0 80,960 107,064 1,450 1,331 1,450 0 0
PARTIELS MONO-METHYL SUCC 0 0 0 0,268 1,838 67,764 1,838 0 116,305
(KG/HR) DIMETHYL SUCCINA 0 0 0 0 0,000 72,391 0,000 0 2,149
PARA-XYLENE 0 0 0 0,042 1528,776 3821,920 1528,190 0,586 0,527
SUCCINIC ACID 0,621 0,000 0,502 0,195 0,000 0,002 0,000 0,000 0,260
METHANOL 0,000 1,000 0,192 0,001 0,002 0,000 0,002 0,000 0,000
FRACTIONS WATER 0,379 0,000 0,306 0,802 0,001 0,000 0,001 0,000 0,000
MASSIQUES MONO-METHYL SUCC 0,000 0,000 0,000 0,002 0,001 0,017 0,001 0,000 0,724
DIMETHYL SUCCINA 0,000 0,000 0,000 0,000 0,000 0,018 0,000 0,000 0,013
PARA-XYLENE 0,000 0,000 0,000 0,000 0,996 0,962 0,996 1,000 0,003
DEBIT MASSIQUE TOTAL (KG/HR) 213,588 50,751 264,339 133,570 1534,810 3973,970 1534,224 0,586 160,724
TEMPERATURE (K) 298,15 298,15 298,15 298,15 298,15 298,15 350,393 298,15 479,376
PRESSION (ATM) 1 1 1 1 1 1 1 1 1
168
Annexe H : Equilibre thermodynamique diphasique des binaires et
ternaires des différents systèmes étudiés.
1 1
0,9 0,9
0,5 0,5
0,4 0,4
0,3 0,3
0,2 0,2
0,1 0,1
0 0
0 0,5 1 0 0,5 1
xAS xAS
0,8 0,8
0,7 0,7
0,6 0,6
0,5 0,5
0,4 0,4
0,3 0,3
0,2 0,2
0,1 0,1
0 0
0 0,5 1 0 0,5 1
xAS xAS
169
1 1
0,9 0,9
Fraction liquide bulle Oct
0,8 0,8
0,7 0,7
0,6 0,6
0,5 0,5
0,4 0,4
0,3 0,3
0,2 0,2
0,1 0,1
0 0
0 0,5 1 0 0,5 1
xOCT xEAU
0,8 0,8
0,7 0,7
0,6 0,6
0,5 0,5
0,4 0,4
0,3 0,3
0,2 0,2
0,1 0,1
0 0
0 0,5 1 0 0,5 1
xEAU xMOE
170
Equilibres liquide-liquide des mélanges ternaires du système acide
succinique/méthanol/eau/mono-octyl succinate/di-octyl succinate.
171
Equilibre liquide-vapeur des quinze mélanges binaires des constituants du
système acide succinique/méthanol/eau/mono-méthyl succinate/di-méthyl
succinate/para-xylène.
1 1
0,9 0,9
0,8 0,8
0,7 0,7
0,6 0,6
yAS
yAS
0,5 0,5
0,4 0,4
0,3 0,3
0,2 0,2
0,1 0,1
0 0
0 0,5 1 0 0,5 1
xAS xAS
yAS
0,5 0,5
0,4 0,4
0,3 0,3
0,2 0,2
0,1 0,1
0 0
0 0,5 1 0 0,5 1
xAS xAS
172
1 1
0,9 0,9
0,8 0,8
0,7 0,7
0,6 0,6
yMeOH
yAS
0,5 0,5
0,4 0,4
0,3 0,3
0,2 0,2
0,1 0,1
0 0
0 0,5 1 0 0,5 1
xAS xMeOH
yMeOH
0,5 0,5
0,4 0,4
0,3 0,3
0,2 0,2
0,1 0,1
0 0
0 0,5 1 0 0,5 1
xMeOH xMeOH
yEAU
0,5 0,5
0,4 0,4
0,3 0,3
0,2 0,2
0,1 0,1
0 0
0 0,5 1 0 0,5 1
xMeOH xEAU
173
1 1
0,9 0,9
0,8 0,8
0,7 0,7
0,6 0,6
yEAU
yEAU
0,5 0,5
0,4 0,4
0,3 0,3
0,2 0,2
0,1 0,1
0 0
0 0,5 1 0 0,5 1
xEAU xEAU
yMMS
0,5 0,5
0,4 0,4
0,3 0,3
0,2 0,2
0,1 0,1
0 0
0 0,5 1 0 0,5 1
xMMS xMMS
0,5
0,4
0,3
0,2
0,1
0
0 0,5 1
xDMS
174
Equilibres liquide-liquide des mélanges ternaires du système acide succinique/méthanol/eau/mono-méthyl succinate/di-
méthyl succinate/para-xylène.
175
NOMENCLATURE
E Débit molaire de la phase organique mol.s-1
F Débit molaire d’alimentation de la charge mol.s-1
€Zc
H Rétention liquide du procédé d’extraction liquide réactive simple mol
Constante d’équilibre chimique
M Mélange des débits F et S ou des débits RN et E1 mol.s-1
…X:
m Variable définie pour les courbes de raffinats réactifs
…†
Quantité totale initiale de mole mol
…+
Nombre total de constituants
à
Nombre total de réaction
Conversion cinétique mol.s-1
R Débit molaire de la phase aqueuse mol.s-1
S Débit molaire d’alimentation du solvant mol.s-1
(S/F) Taux de solvant
(S/F) min Taux de solvant minimum
¦
t Temps s
Variance
`
w Composition molaire de la rétention H
F
Composition molaire réactive
Composition molaire
z Composition molaire de la charge F
%&'
Z Composition molaire réactive de la charge F
Rendement de fermentation
T Température K
P Pôle opératoire
n Nombre de mole d’un constituant mol
y Composition molaire de la phase vapeur
Y Composition molaire réactive de la phase vapeur
r Cinétique de réaction mol.s-1
Zr Terme cinétique de réaction pour les bilans matières
U Rétention liquide totale d’un étage de la colonne m-3
k Constante cinétique chimique mol.s-1
M
Ea Énergie d’activation J.mol-1.K-1
Á
Coefficient d’activité
U
Taux de partage, fraction de la phase I dans le mélange
¿
Avancement de réaction
d
Temps adimensionnel
Coefficient stœchiométrique
•
Matrice des coefficients stœchiométriques des constituants de
¿
référence
í
Taux de récupération
Taux de conversion
Exposants
I Phase aqueuse
II Phase organique
176
Indices
i Constituants
j étage
L Limite
N Dernier étage de la colonne
r réaction
ref référence
eq Equilibre
Acronymes
DMS Di-Methyl Succinate
DOS Di-octyl Succinate
MeOH Méthanol
MMS Mono-Méthyl Succinate
MOS Mono-Octyl Succinate
NET Nombre d’étages théorique
OcOH Octanol
P-X Para-Xylène
SMILES Simplified Molecular Input Line Entry Specification
AS Acide succinique
R-OH Alcool
ME Mono-ester
DE Di-ester
C4 Molécules avec 4 carbone (acide succinique, mono-ester, di-ester)
177
LISTES DES FIGURES
Chapitre I
178
Figure I-24: Organigramme de présélection d’un extracteur liquide-liquide (Technique de
l’ingénieur J2766) .......................................................................................................................................... 43
Figure I-25: Exemples de mélangeurs-décanteurs ................................................................................ 44
Figure I-26: Schéma d’une colonne d’extraction liquide-liquide réactive ............................. 48
Figure I-27 : Construction graphique pour une colonne d'extraction liquide-liquide réactive
(Olán-Acosta et al., 2014)............................................................................................................................ 55
Figure I-29 : Limite de faisabilité d'une colonne réactive à contre-courant .................................. 57
Figure I-30 : Détermination graphique du taux de solvant minimum ............................................. 57
Chapitre II
179
Figure II-15 : Construction type McCabe et Thiele pour illustrer le critère d’arrêt de la
détermination du nombre d’étages théoriques. ...................................................................................... 95
Figure II-16 : Construction type McCabe et Thiele pour illustrer le critère d’arrêt de la
détermination du nombre d’étage théorique et la position de l’alimentation, colonne à reflux.
............................................................................................................................................................................. 96
Figure II-17: Organigramme détaillé de la méthode de conception d’une extraction
liquide-liquide réactive ............................................................................................................................102
Chapitre III
180
Figure III-13 : Influence de la rétention liquide sur les critères choisis ...............................120
Figure III-14 : Diagramme liquide-liquide du système acide succinique
/méthanol/eau/mono-ester/di-ester/para-xylène......................................................................123
Figure III-15 – Décomposition de l’enveloppe de phase .............................................................123
Figure III-16 : Analyse de faisabilité, courbe de raffinat réactif du système Acide
succinique /méthanol/eau/mono-ester/di-ester/para-xylène. ..............................................124
Figure III-17-a : Représentation 3d de l’analyse de faisabilité .................................................126
Figure III-18 : Détermination du nombre d’étages théoriques pour le système avec Nc=6
et Nr=2 ............................................................................................................................................................128
Figure III-19 : Profil de composition réactive en acide succinique, détermination du NET,
NET=5 .............................................................................................................................................................128
Figure III-20 : Diagramme de distribution pour la détermination du NET d’une façon
analogue à McCabe Thièle .......................................................................................................................129
Figure III-21: Influence de la cinétique chimique sur les performances de la colonne
d’extraction ...................................................................................................................................................131
Figure III-22 : Influence du nombre d’étages théoriques sur les critères choisis,
simulation à l’équilibre.............................................................................................................................132
Figure III-23 : Influence du nombre d’étages théoriques sur les critères choisis,
simulation hors équilibre ........................................................................................................................132
Figure III-24 : Influence de la composition du solvant sur les critères choisis, simulation
à l’équilibre ...................................................................................................................................................133
Figure III-25 : Influence de la composition du solvant sur les critères choisis, simulation
hors équilibre ...............................................................................................................................................133
Figure III-26 : Influence de la composition du solvant sur les critères choisis, simulation
à l’équilibre ...................................................................................................................................................134
Figure III-27 : Influence de la composition du solvant sur les critères choisis, simulation
hors équilibre ...............................................................................................................................................134
Figure III-28 : Influence de la rétention liquide sur les critères choisis ...............................135
Figure III-29 : Schéma de simulation du procédé d’extraction pour la récupération
d’acide succinique ......................................................................................................................................139
Figure III-30 : Schéma de simulation du procédé d’extraction réactive avec 6
constituants pour la récupération d’acide succinique .................................................................143
Figure III-32 : Débit partiel du courant de production de molécule C4 .................................145
181
Figure III-33 : Consommation énergétique globale des procédés étudiés...........................146
Figure III-34 : Consommation énergétique par kg/h de molécule produit des procédés
étudiés ............................................................................................................................................................146
182
LISTE DES TABLEAUX
Chapitre I
Chapitre II
Chapitre III
185
Tableau III-23 : Résultats de simulation du procédé d’extraction liquide réactive pour la
récupération de l’acide succinique ......................................................................................................142
Tableau III-24 : Paramètres des différentes opérations unitaires du procédé d’extraction
liquide réactive ............................................................................................................................................144
Tableau III-25 : Résultats de simulation du procédé d’extraction liquide réactive pour la
récupération de l’acide succinique ......................................................................................................144
Tableau III-26 : Température d’ébullition des constituants du système Nc=6 ..................144
186
RES FES RENCES
A
A Tuchlenski, A.B., 2001. Reactive Distillation Industrial Applications: Process Design and Scale-
Up. Chemical Engineering Science 56, 387–394. doi:10.1016/S0009-2509(00)00240-2
Agar, D.W., 1999. Multifunctional reactors: Old preconceptions and new dimensions. Chemical
Engineering Science, 1st International symposium on multifunctional reactors 54, 1299–1305.
doi:10.1016/S0009-2509(99)00040-8
B
Barbosa, D., Doherty, M.F., 1988. The simple distillation of homogeneous reactive mixtures.
Chemical Engineering Science 43, 541–550. doi:10.1016/0009-2509(88)87015-5
Barbosa, D., Doherty, M.F., 1987. A New Set of Composition Variables for the Representation of
Reactive-Phase Diagrams. Proceedings of the Royal Society of London A: Mathematical,
Physical and Engineering Sciences 413, 459–464. doi:10.1098/rspa.1987.0126
Beauprez, J.J., De Mey, M., Soetaert, W.K., 2010. Microbial succinic acid production: natural
versus metabolic engineered producers. Process Biochemistry 45, 1103–1114.
Berglund, K.A., Yedur, S., Dunuwila, D.D., 2008. Succinic acid production and purification.
CA2301177 C.
Bevia, F.R., Rico, D.P., 1983. Quaternary liquid-liquid equilibria: Experimental determination and
correlation of equilibrium data. Part III. New methods of representation and correlation of
liquid-liquid equilibria in quaternary systems. Fluid Phase Equilibria 10, 115–124.
doi:10.1016/0378-3812(83)80008-9
188
C
CEA, 2008. Le traitement-recyclage du combustible nucléaire usé [WWW Document]. CEA/CEA.
URL http://www.cea.fr/Pages%2fdomaines-recherche%2fenergies%2fenergie-
nucleaire%2ftraitement-recyclage-combustible-nucleaire-use.aspx (accessed 8.25.16).
Cusack, Roger W; Glatz, Donald J., 1996. Apply Liquid-Liquid Extraction to Today’s Problems -
ProQuest. Chemical Engineering Journal 103, 94–103.
D
D. Barbosa, M.F.D., n.d. Design and minimum reflux calculations for multicomponent reactive
distillation columns.
Datta, R., Glassner, D.A., Jain, M.K., Roy, J.R.V., 1992. Fermentation and purification process for
succinic acid.
Dudáš, J., Kotora, M., Bradáč, M., Markoš, J., 2014. Design consideration of dimethyl succinate
production process. Chem. Pap. 68, 1667–1677. doi:10.2478/s11696-014-0580-4
E
Engasser, J.M., Biocatyle et bioprocédés : Poux, M., Cognet, P., Gourdon, C., Génie des procédés
durables du concept à la concrétisation industrielle, DUNOD, Paris, 2010, p. 365-389. ISBN
978-2-10-051605-6
F
F. Ruiz Beviá, D.P.R., 1984. Liquid-liquid extraction: A graphical method for equilibrium stage
calculations for quaternary systems. Fluid Phase Equilibria 15, 257–265. doi:10.1016/0378-
3812(84)87011-9
189
G
Garzón, C.S.L., Straathof, A.J., 2014. Recovery of carboxylic acids produced by fermentation.
Biotechnology advances 32, 873–904.
Ghosh, T.K., Lin, H.D., Hines, A.L., 1993. Hybrid adsorption-distillation process for separating
propane and propylene. Ind. Eng. Chem. Res. 32, 2390–2399. doi:10.1021/ie00022a024
Glassner, D.A., Elankovan, P., Beacom, D.R., Berglund, K.A., 1995. Purification process for
succinic acid produced by fermentation. Appl Biochem Biotechnol 51–52, 73–82.
doi:10.1007/BF02933412
Gundepsen, T., Naess, L., 1988. The synthesis of cost optimal heat exchanger networks: An
industrial review of the state of the art. Computers & Chemical Engineering 12, 503–530.
doi:10.1016/0098-1354(88)87002-9
Guettler, M.V., Rumler, D., Jain, M.K., 1999. Actinobacillus succinogenes sp. nov., a novel
succinic-acid-producing strain from the bovine rumen. International journal of systematic
bacteriology 49, 207–216.
H
Hamid, M.K.A., Sin, G., Gani, R., 2010. Integration of process design and controller design for
chemical processes using model-based methodology. Computers & Chemical Engineering,
Selected Paper of Symposium ESCAPE 19, June 14-17, 2009, Krakow, Poland 34, 683–699.
doi:10.1016/j.compchemeng.2010.01.016
Hauan, S., Ciric, A.R., Westerberg, A.W., Lien, K.M., 2000. Difference points in extractive and
reactive cascades. I – Basic properties and analysis. Chemical Engineering Science 55, 3145–
3159. doi:10.1016/S0009-2509(99)00591-6
Hoffmaster, W.R., Hauan, S., 2004. Difference points in reactive and extractive cascades. III—
Properties of column section profiles with arbitrary reaction distribution. Chemical
Engineering Science 59, 3671–3693. doi:10.1016/j.ces.2004.05.013
Hoffmaster, W.R., Hauan, S., 2005. Difference points in reactive and extractive cascades: IV—
Feasible regions for multisection columns with kinetic reactions and side streams. Chemical
Engineering Science 60, 7075–7090. doi:10.1016/j.ces.2005.06.011
190
Huss, R.S., Chen, F., Malone, M.F., Doherty, M.F., 1999. Computer-aided tools for the design of
reactive distillation systems. Computers & Chemical Engineering 23, S955–S962.
doi:10.1016/S0098-1354(99)80230-0
K
Klaassen, R., Feron, P.H.M., Jansen, A.E., 2005. Process Developments in Membrane Contactors
and ReactorsMembrane Contactors in Industrial Applications. Chemical Engineering
Research and Design 83, 234–246. doi:10.1205/cherd.04196
Krzyżaniak, A., Leeman, M., Vossebeld, F., Visser, T.J., Schuur, B., de Haan, A.B., 2013. Novel
extractants for the recovery of fermentation derived lactic acid. Separation and purification
technology 111, 82–89.
Kuceski, V.P., 1995. Purification of alpha, omega alkanedioic acids by partial esterification.
US2824122 A.
L
Lawson, K.W., Lloyd, D.R., 1997. Membrane distillation. Journal of Membrane Science 124, 1–
25. doi:10.1016/S0376-7388(96)00236-0
Lee, J.W., Hauan, S., Lien, K.M., Westerberg, A.W., 2000. Difference points in extractive and
reactive cascades. II — Generating design alternatives by the lever rule for reactive systems.
Chemical Engineering Science 55, 3161–3174. doi:10.1016/S0009-2509(99)00590-4
Lee, P.C., Lee, S.Y., Chang, H.N., 2008. Succinic acid production by Anaerobiospirillum
succiniciproducens ATCC 29305 growing on galactose, galactose/glucose, and
galactose/lactose. J. Microbiol. Biotechnol. 18, 1792–1796.
Lee, P.C., Lee, S.Y., Hong, S.H., Chang, H.N., 2003. Batch and continuous cultures of
Mannheimia succiniciproducens MBEL55E for the production of succinic acid from whey
and corn steep liquor. Bioprocess Biosyst Eng 26, 63–67. doi:10.1007/s00449-003-0341-1
Lee, P.C., Lee, S.Y., Hong, S.H., Chang, H.N., 2002. Isolation and characterization of a new
succinic acid-producing bacterium, Mannheimia succiniciproducensMBEL55E, from bovine
rumen. Appl. Microbiol. Biotechnol. 58, 663–668. doi:10.1007/s00253-002-0935-6
Lee, P.C., Lee, W.G., Kwon, S., Lee, S.Y., Chang, H.N., 1999. Succinic acid production by
Anaerobiospirillum succiniciproducens: effects of the H 2/CO 2 supply and glucose
concentration. Enzyme and Microbial Technology 24, 549–554.
191
Lee, S.J., Lee, D.-Y., Kim, T.Y., Kim, B.H., Lee, J., Lee, S.Y., 2005. Metabolic engineering of
Escherichia coli for enhanced production of succinic acid, based on genome comparison and
in silico gene knockout simulation. Applied and environmental microbiology 71, 7880–7887.
Lee, S.Y., Kim, J.M., Song, H., Lee, J.W., Kim, T.Y., Jang, Y.-S., 2008. From genome sequence to
integrated bioprocess for succinic acid production by Mannheimia succiniciproducens.
Applied microbiology and biotechnology 79, 11–22.
Leybros., J., Extraction liquide-liquide – Choix, calcul et conception des appareils, Techniques de
l’ingénieur référence J2766, décembre 2004.
Li, Q., Wang, D., Wu, Y., Li, W., Zhang, Y., Xing, J., Su, Z., 2010. One step recovery of succinic
acid from fermentation broths by crystallization. Separation and Purification Technology 72,
294–300. doi:10.1016/j.seppur.2010.02.021
Lin, H., Bennett, G.N., San, K.-Y., 2005. Metabolic engineering of aerobic succinate production
systems in Escherichia coli to improve process productivity and achieve the maximum
theoretical succinate yield. Metab. Eng. 7, 116–127. doi:10.1016/j.ymben.2004.10.003
Lin, S.K.C., Du, C., Blaga, A.C., Camarut, M., Webb, C., Stevens, C.V., Soetaert, W., 2010. Novel
resin-based vacuum distillation-crystallisation method for recovery of succinic acid crystals
from fermentation broths. Green Chemistry 12, 666. doi:10.1039/b913021g
Linnhoff, B., Hindmarsh, E., 1983. The pinch design method for heat exchanger networks.
Chemical Engineering Science 38, 745–763. doi:10.1016/0009-2509(83)80185-7
Luque, R., Lin, C.S.K., Du, C., Macquarrie, D.J., Koutinas, A., Wang, R., Webb, C., Clark, J.H.,
2009. Chemical transformations of succinic acid recovered from fermentation broths by a
novel direct vacuum distillation-crystallisation method. Green Chem. 11, 193–200.
doi:10.1039/B813409J
M
Mills, P.L., Quiram, D.J., Ryley, J.F., 2007. Microreactor technology and process miniaturization
for catalytic reactions—A perspective on recent developments and emerging technologies.
Chemical Engineering Science, 8th International Conference on Gas-Liquid and Gas-Liquid-
Solid Reactor Engineering 62, 6992–7010. doi:10.1016/j.ces.2007.09.021
Mims, S.S., 1978. Process for treatment of adipic acid mother liquor. US4076948 A.
192
Minotti, M., Doherty, M.F., Malone, M.F., 1998. Design for Simultaneous Reaction and
Liquid−Liquid Extraction. Ind. Eng. Chem. Res. 37, 4748–4755. doi:10.1021/ie970786i
Minotti, M., Doherty, M.F., Malone, M.F., 1996. A Geometric Method for the Design of Liquid
Extractors. Ind. Eng. Chem. Res. 35, 2672–2681. doi:10.1021/ie950687+
N
Niang, M., Mikitenko, P., 1998. Analyse de faisabilité d’une distillation réactive par la méthode
des courbes de résidus. Oil & Gas Science and Technology 53, 439–462.
doi:10.2516/ogst:1998037
O
Olán-Acosta, M.A., Rico-Ramírez, V., González-Alatorre, G., Castrejón-González, E.O., 2014. A
graphical design method for reaction–extraction processes in quaternary systems. Chemical
Engineering Research and Design 92, 2027–2040.
Orjuela, A., 2010. Separation of succinic acid from fermentation broths and esterification by a
reactive distillation method. Michigan State University.
Orjuela, A., Kolah, A., Hong, X., Lira, C.T., Miller, D.J., 2012a. Diethyl succinate synthesis by
reactive distillation. Separation and Purification Technology 88, 151–162.
doi:10.1016/j.seppur.2011.11.033
Orjuela, A., Orjuela, A., Lira, C.T., Miller, D.J., 2013. A novel process for recovery of
fermentation-derived succinic acid: Process design and economic analysis. Bioresource
Technology 139, 235–241. doi:10.1016/j.biortech.2013.03.174
Orjuela, A., Yanez, A.J., Evans, J., Hassan, A.M., Miller, D.J., Lira, C.T., 2011a. Phase equilibria
in binary mixtures with monoethyl succinate. Fluid Phase Equilibria 309, 121–127.
Orjuela, A., Yanez, A.J., Peereboom, L., Lira, C.T., Miller, D.J., 2011b. A novel process for
recovery of fermentation-derived succinic acid. Separation and Purification Technology 83,
31–37. doi:10.1016/j.seppur.2011.08.010
Orjuela, A., Yanez, A.J., Santhanakrishnan, A., Lira, C.T., Miller, D.J., 2012b. Kinetics of mixed
succinic acid/acetic acid esterification with Amberlyst 70 ion exchange resin as catalyst.
Chemical Engineering Journal 188, 98–107. doi:10.1016/j.cej.2012.01.103
Oxley, P., Brechtelsbauer, C., Ricard, F., Lewis, N., Ramshaw, C., 2000. Evaluation of Spinning
Disk Reactor Technology for the Manufacture of Pharmaceuticals. Ind. Eng. Chem. Res. 39,
2175–2182. doi:10.1021/ie990869u
193
P
Pappu, V.K.S., Kanyi, V., Santhanakrishnan, A., Lira, C.T., Miller, D.J., 2013. Butyric acid
esterification kinetics over Amberlyst solid acid catalysts: The effect of alcohol carbon chain
length. Bioresource Technology 130, 793–797. doi:10.1016/j.biortech.2012.12.087
Pham, H.N., Doherty, M.F., 1990a. Design and synthesis of heterogeneous azeotropic
distillations—I. Heterogeneous phase diagrams. Chemical Engineering Science 45, 1823–
1836. doi:10.1016/0009-2509(90)87058-Z
Pham, H.N., Doherty, M.F., 1990b. Design and synthesis of heterogeneous azeotropic
distillations—II. Residue curve maps. Chemical Engineering Science 45, 1837–1843.
doi:10.1016/0009-2509(90)87059-2
Pham, H.N., Doherty, M.F., 1990c. Design and synthesis of heterogeneous azeotropic
distillations—III. Column sequences. Chemical Engineering Science 45, 1845–1854.
doi:10.1016/0009-2509(90)87060-6
Powell, M. J. D., "A Fortran Subroutine for Solving Systems of Nonlinear Algebraic
Equations," Numerical Methods for Nonlinear Algebraic Equations, P. Rabinowitz, ed., Ch.7,
1970.
R
RICE, R.G., 1981. Adsorptive Distillation. Chemical Engineering Communications 10, 111–126.
doi:10.1080/00986448108910929
S
Samant, K.D., Ng, K.M., 1998a. Synthesis of extractive reaction processes. AIChE journal 44,
1363–1381.
Samant, K.D., Ng, K.M., 1998b. Design of multistage extractive reaction processes. AIChE journal
44, 2689–2702.
194
Shampine, L.F., Reichelt, M.W., 2006. The MATLAB ODE Suite. SIAM Journal on Scientific
Computing. doi:10.1137/S1064827594276424
Sie, S.T., Krishna, R., 1998. Process development and scale up: 1. Process development strategy
and methodology. Reviews in Chemical Engineering 14, 47–88.
Smith, R.M., 1999. Supercritical fluids in separation science – the dreams, the reality and the
future. Journal of Chromatography A 856, 83–115. doi:10.1016/S0021-9673(99)00617-2
Smith, W.R., 1980. The Computation of Chemical Equilibria in Complex Systems. Ind. Eng.
Chem. Fund. 19, 1–10. doi:10.1021/i160073a001
Solomons, T.W. Graham; Fryhle, Craig B. (2003), Organic Chemistry (8th ed.), Wiley, ISBN 0-
471-41799-8
Song, H., Lee, S.Y., 2006. Production of succinic acid by bacterial fermentation. Enzyme and
Microbial Technology 39, 352–361.
Srikhirin, P., Aphornratana, S., Chungpaibulpatana, S., 2001. A review of absorption refrigeration
technologies. Renewable and Sustainable Energy Reviews 5, 343–372. doi:10.1016/S1364-
0321(01)00003-X
Stankiewicz, A.I., Moulijn, J.A., others, 2000. Process intensification: transforming chemical
engineering. Chemical Engineering Progress 96, 22–34.
Stichlmair, J., 1980. Leistungs- und Kostenvergleich verschiedener Apparatebauarten für die
Flüssig/Flüssig-Extraktion. Chemie Ingenieur Technik 52, 253–255.
doi:10.1002/cite.330520313
Subramaniam, B., McHugh, M.A., 1986. Reactions in supercritical fluids - a review. Ind. Eng.
Chem. Proc. Des. Dev. 25, 1–12. doi:10.1021/i200032a001
T
Taylor, R., Krishna, R., 2000. Modelling reactive distillation. Chemical Engineering Science 55,
5183–5229. doi:10.1016/S0009-2509(00)00120-2
Teixeira, M.A., Rodríguez, O., Mata, V.G., Rodrigues, A.E., 2009. Perfumery quaternary diagrams
for engineering perfumes. AIChE J. 55, 2171–2185. doi:10.1002/aic.11811
Thery, R., 2002. Analyse de faisabilité, synthèse et conception de procédés de distillation réactive.
Toulouse, INPT.
195
Thery, R., Meyer, X.-M., Joulia, X., 2005. Analyse de faisabilité, synthèse et conception des
procédés de distillation réactive : état de l’art et analyse critique. Can. J. Chem. Eng. 83, 242–
266. doi:10.1002/cjce.5450830211
Thery, R., Meyer, X.M., Joulia, X., Meyer, M., 2005. Preliminary Design of Reactive Distillation
Columns. Chemical Engineering Research and Design 83, 379–400. doi:10.1205/cherd.04112
Treybal, R.E., 1951. Liquid Extraction. McGraw Hill Book Company Inc.
Turnberg, S.L., Hauan, S., 2009. Difference points in reactive and extractive cascades, V:
Determining reaction distribution from separatrix manifolds. Chemical Engineering Science
64, 4846–4855. doi:10.1016/j.ces.2009.07.019
U
Ung, S., Doherty, M.F., 1995a. Vapor-liquid phase equilibrium in systems with multiple chemical
reactions. Chemical Engineering Science 50, 23–48.
Ung, S., Doherty, M.F., 1995b. Calculation of residue curve maps for mixtures with multiple
equilibrium chemical reactions. Ind. Eng. Chem. Res. 34, 3195–3202.
doi:10.1021/ie00037a003
V
Vemuri, G.N., Eiteman, M.A., Altman, E., 2002. Succinate production in dual-phase Escherichia
coli fermentations depends on the time of transition from aerobic to anaerobic conditions. J.
Ind. Microbiol. Biotechnol. 28, 325–332. doi:10.1038/sj/jim/7000250
Victor H. Agreda, Lee R. Partin, 1984. Reactive distillation process for the production of methyl
acetate. US4435595 A.
W
Wang, Y., Gong, C., Sun, J., Gao, H., Zheng, S., Xu, S., 2010. Separation of ethanol/water
azeotrope using compound starch-based adsorbents. Bioresource Technology 101, 6170–
6176. doi:10.1016/j.biortech.2010.02.102
Wibowo, C., Ng, K.M., 2002. Visualization of high-dimensional phase diagrams of molecular and
ionic mixtures. AIChE J. 48, 991–1000. doi:10.1002/aic.690480509
Wu, R.C., Xu, Y.Z., Song, Y.Q., Luo, J.A., Liu, D., 2011. A novel strategy for salts recovery from
1,3-propanediol fermentation broth by bipolar membrane electrodialysis. Separation and
Purification Technology 83, 9–14. doi:10.1016/j.seppur.2011.06.028
196
Wuithier, P., Giraud, A., 1972a. Le Pétrole: raffinage et génie chimique. tome I. Editions Technip,
Paris.
Wuithier, P., Giraud, A., 1972b. Le Pétrole: raffinage et génie chimique. tome II. Editions Technip,
Paris.
Wuithier, P., Giraud, A., 1965. Le pétrole: raffinage et génie chimique. Technip.
Z
Zeikus, J.G., 1980. Chemical and Fuel Production by Anaerobic Bacteria. Annual Review of
Microbiology 34, 423–464. doi:10.1146/annurev.mi.34.100180.002231
Zeikus, J.G., Jain, M.K., Elankovan, P., 1999. Biotechnology of succinic acid production and
markets for derived industrial products. Appl Microbiol Biotechnol 51, 545–552.
doi:10.1007/s002530051431
197
Les lauriers c’est pour les autres. Les lauriers je m’en bats les c… (Joey Starr, Pensées
répliques & autres tac-tac, 2016)
198
RES SUMES /ABSTRACT
Le couplage fonctionnel des opérations de séparation et de réaction ainsi que les
bioprocédés sont deux axes de recherche largement explorés. Cependant, l’industrie du
génie des procédés a du mal à se tourner vers des technologies de ce type car il demeure
un réel manque de connaissances et d’outils de conception pour ce genre de procédés.
Une méthodologie de conception générale pour l'extraction liquide-liquide
réactive est introduite dans cette étude. Elle est composée de trois étapes différentes:
l'analyse de faisabilité, la synthèse ou dimensionnement du procédé et la validation par
simulation. Cette méthodologie conduit à des paramètres structuraux et opératoires de
la colonne étudiée à partir seulement des informations concernant le comportement
physico-chimique du système étudié, en exploitant les équations d’équilibre chimique et
entre phase ainsi que les bilans matières. Les résultats de cette méthode sont un bon
point de départ pour une étude d'optimisation ou d'un processus de calcul
d'investissement. Cette méthodologie a été appliquée à différentes études de cas:
regroupant deux stratégies différentes d'extraction avec plusieurs solvants pour
récupérer l'acide succinique dans un milieu de fermentation.