Study
Study
Titre:
caractérisation de l'hydrodynamique d'une cellule et relation à la
Title:
performance
Auteur:
Laura Bouchra Orabi
Author:
Date: 2018
Type: Mémoire ou thèse / Dissertation or Thesis
Orabi, L. B. (2018). Électrodialyse de la liqueur noire avec membrane bipolaire-
Référence: caractérisation de l'hydrodynamique d'une cellule et relation à la performance
Citation: [Mémoire de maîtrise, École Polytechnique de Montréal]. PolyPublie.
[Link]
URL de PolyPublie:
[Link]
PolyPublie URL:
Directeurs de
recherche: Louis Fradette
Advisors:
Programme:
Génie chimique
Program:
RELATION À LA PERFORMANCE
(GÉNIE CHIMIQUE)
JANVIER 2018
UNIVERSITÉ DE MONTRÉAL
Ce mémoire intitulé :
PERFORMANCE
DÉDICACE
À Walid
Le compagnon de mes aventures
et Karim
Je suis fière de l’homme que tu deviens
iv
REMERCIEMENTS
Je remercie tout d’abord mon directeur de recherche, M. Louis Fradette de m’avoir prodigué
des conseils et des commentaires scientifiques toujours judicieux, pour la confiance qu’il m’a
accordée, de m’avoir supporté, encouragé et orienté tout au long de ce projet. Sans vous, ce
mémoire n’aurait jamais pu être mené à terme.
Merci aux membres de jury de prendre le temps d’évaluer mon projet de recherche pour l’obtention
du diplôme de maîtrise.
J'adresse également mes sincères remerciements à toute l'équipe technique du département de génie
chimique. À Daniel Pilon, Martine Lamarche, Gino Robin et Jean Huard.
Je remercie spécialement Ahmad Al-Haiek pour son soutien et sa présence. Sans qui cette maîtrise
n'aurait surement pas été la même.
J'adresse également une pensée particulière à mes amis et partenaires de bureau: Fatma, Adrian,
Christine, Olivier, Rémi, David, Igor, Jean-Michel, Faezeh, Bing, Anna. Vous êtes tous une réelle
source d’inspiration.
J’adresse mes plus affectueuses pensées à Aida et Ghassan pour leur soutien et leur présence. Pour
m’avoir encouragé tout au long de cette maîtrise. Je suis chanceuse de vous avoir dans ma vie.
J’aimerais surtout remercier du fond de mon cœur mon conjoint Walid d’être toujours présent pour
m’épauler à travers toutes ces belles années. Merci également à mon fils, Karim, ton amour et ton
soutien font toute la différence dans ma vie. Sans vous deux, je ne serai jamais où je suis en ce
moment. Je vous aime.
v
RÉSUMÉ
L’industrie de pâte et papier canadienne existe depuis les années 1800 et a connu des évolutions
révolutionnaires au fil du temps. En 1990, le Canada est devenu le premier producteur mondial de
pâte et papier dans le monde. Cette industrie est demeurée une composante fondamentale de
l’économie canadienne jusqu’à l’émergence des technologies électroniques au début des années
2000. Pour s’adapter avec la nouvelle réalité, l’industrie des P&P devra pouvoir rénover son
modèle d’affaire afin de diversifier son portfolio de produits et de varier ses sources de revenus.
L’extraction de la lignine à partir de la liqueur usée de cuisson dans une usine de P&P s’avère une
voie potentiellement rentable pour convertir un résidu du procédé en bénéfice. La lignine est
produite en quantité excédant 50 Mtonnes dans le monde chaque année. De nombreuses études
examinent la possibilité de transformer ce polymère naturel en divers produits à valeur ajoutée tels
que, les fibres de carbone, le charbon actif, les adhésifs, les phénols dans le but de trouver des
substituants aux produits dérivés du pétrole tout en créant une nouvelle richesse pour les usines.
Plusieurs méthodes d’extraction de la lignine ont été développées au cours des dernières années.
En effet, la dissolution de la lignine dans la liqueur noire est liée au pH. La lignine se retrouve
dissoute avec un pH fortement alcalin et commence à précipiter lorsque le pH diminue en bas de
10.9. C’est pourquoi les méthodes d’extraction reposent sur l’acidification de la liqueur noire.
Une méthode novatrice d’électrodialyse à membrane bipolaire commence à être étudiée de plus en
plus à travers le monde en tant que technique d’acidification de liqueur noire avec récupération de
NaOH. Cette méthode présente des avantages importants par rapport aux différentes méthodes
d’extraction disponibles. Cependant, plusieurs limitations empêchent l’opération prolongée de
l’électrodialyse et réduisent l’efficacité de la cellule. La limitation majeure est l’encrassement de
la membrane. Il s’agit d’un phénomène de déposition des particules de lignine précipitée sur la
surface de la membrane causant ainsi une augmentation dans la résistance électrique dans la cellule
et causant l’arrêt éventuel de l’opération. Il est donc primordial de pouvoir trouver une solution
pour remédier à l’encrassement de la membrane dans l’objectif d’aller vers l’avant avec le
développement de cette technologie.
Tout d’abord, une cellule électrochimique à membrane bipolaire a été conçue et construite à partir
de zéro. Par la suite, une caractérisation hydrodynamique détaillée a été réalisée à l’aide de
certaines solutions références dans le but de mieux comprendre le comportement du fluide à
l’intérieur de la cellule. Finalement, l’électrodialyse de la liqueur noire a été effectuée pour
différents régimes d’écoulement.
Les résultats obtenus ont montré que la promotion de la turbulence dans la cellule a permis de
prolonger la durée d’opération d’environ 25%. Ce résultat pourrait être expliqué à l’aide de la
notion des couches limites de diffusion. En effet, durant les opérations membranaires, il y a une
couche limite de diffusion qui se forme à la surface de la membrane. L’épaisseur de cette couche
est liée au degré de turbulence près de la membrane. Plus la turbulence est importante, moins le
serait l’épaisseur de la couche et vice-versa.
ABSTRACT
The Canadian Pulp and Paper industry has begun in the 1800s, and has undergone revolutionary
evolution over the years. In 1990, Canada became the world's largest producer of pulp and paper
in the world. This industry remained a fundamental component of the Canadian economy until the
emergence of electronic technologies in the 2000s.
To adapt to the new reality, the pulp and paper industry will have to renovate its business model to
diversify its product portfolio and to vary its sources of revenue.
The extraction of lignin from spent cooking liquor in a P&P plant is a potentially cost-effective
way to convert a process residue into a profit. Lignin is produced in quantities exceeding 50
Mtonnes worldwide each year. Many studies are investigating the possibility of converting this
natural polymer into various value-added products such as carbon fiber, activated carbon,
adhesives, phenols in order to find substituents for petroleum products while creating a new wealth
for the factory.
Several methods of lignin extraction have been developed in recent years. In fact, the dissolution
of lignin in the black liquor is related to the pH. Lignin is dissolved in an alkaline pH and begins
to precipitate as the pH decreases below 10.9. Therefore, the extraction methods are based on the
acidification of the black liquor.
The general objective of this research project is to investigate whether the hydrodynamics in the
cell could have an impact on the fouling of the membrane and consequently on the duration of the
operation. Several steps were needed to achieve this goal. First, a bipolar membrane
electrochemical cell was designed and constructed from scratch. Subsequently, a detailed
viii
hydrodynamic characterization was carried out using some reference solutions to better understand
the behavior of the fluid inside the cell. Finally, electrodialysis of the black liquor was performed
for laminar and turbulent flow regimes.
The obtained results showed that the promotion of turbulence in the cell allowed the duration of
operation to be extended by approximately 25%. This result could be explained by the concept of
diffusion boundary layer. In fact, there is a diffusion boundary layer that forms on the surface of
the membrane during membrane operations. The thickness of this layer is related to the degree of
turbulence near the membrane. The greater the turbulence, the less is the thickness of the layer and
vice versa.
The main recommendation for future work is to push the hydrodynamic analysis of black liquor
electrodialysis in a turbulent regime in order to find optimal operating parameters that can delay or
even prevent the membrane fouling in the cell.
ix
REMERCIEMENTS ..................................................................................................................... IV
RÉSUMÉ........................................................................................................................................ V
ABSTRACT .................................................................................................................................VII
1.2 Problématique................................................................................................................... 5
2.1.1 La lignine...................................................................................................................... 7
BIBLIOGRAPHIE ........................................................................................................................ 63
xii
Tableau 2-2 - Pourcentage massique des matières organiques dans la liqueur noire ...................... 9
Tableau 4-4 - Sommaire des relations Ne-Re pour les trois géométries ........................................ 54
Tableau 4-7 - Masse volumique et teneur en lignine de la liqueur noire utilisée .......................... 56
Figure 4-7 - Perte de charge en fonction du débit pour la géométrie III ........................................ 51
Figure 4-10- La viscosité en fonction du taux de cisaillement pour différents pH de la liqueur noire
........................................................................................................................................................ 57
Figure 4-11 - Évolution de la résistance électrique pour les différentes géométries avec régimes
laminaire et turbulent ..................................................................................................................... 59
xiv
Dh Diamètre hydraulique m
L Longueur caractéristique m
m Masse Kg
P Perte de charge Pa
Q Débit m3/s
V Volume m3
xv
NOMBRES ADIMENSIONNELS :
SYMBOLE DESCRIPTION
Ne Nombre de Newton
Re Nombre de Reynolds
SYMBOLES GRECS :
∆P Perte de charge Pa
μ Viscosité Pa.s
τ Contrainte de cisaillement Pa
xvi
ABRÉVIATIONS
ABRÉVIATION DESCRIPTION
ED Électrodialyse
LN Liqueur noire
CHAPITRE 1 INTRODUCTION
Afin de garantir son avenir, l’industrie des P&P devra pouvoir s’adapter avec la nouvelle
réalité en réinventant ses modèles d’affaires afin de trouver de nouvelles sources de revenus. De
plus, les préoccupations internationales par rapport au changement climatique incitent les
industries à revalider leurs procédés dans le but de réduire leur empreinte carbone.
L’avenir de ces usines passe notamment par la diversification de leur portfolio afin de produire des
biocarburants de seconde génération, des bioproduits à forte valeur ajoutée et des biomatériaux. Ce
concept est appelé le bioraffinage.
Une des voies potentiellement rentables dans une bioraffinerie est l’exploitation de la
lignine comme matière première pour fabriquer des produits à valeur ajoutée. La lignine est l’un
des polymères naturels aromatiques les plus abondants. Elle constitue entre 18-35% du poids du
bois. (Saito, et al., 2012) La lignine est un sous-produit de l’industrie papetière et se retrouve
dissoute dans la liqueur noire après la cuisson du bois. Elle est produite en quantité excédant les 50
Mtonnes dans le monde chaque année. (Smolarski, 2012) Cependant, la majorité de ce sous-produit
2
est actuellement utilisé comme combustible de faible valeur afin de répondre aux besoins
énergétiques de l’usine. (Saito, et al., 2012) Par conséquent, une attention importante a été accordée
aux possibilités d’extraction de ce sous-produit afin de l’utiliser à grande échelle dans diverses
applications.
Le procédé kraft commence par la cuisson des copeaux de bois dans un digesteur à haute
température (150 °C) et haute pression (7 à 7.5 bars) dans une solution appelée « La liqueur
blanche ». La composition typique des solutés contenus dans la liqueur blanche est présentée dans
le Tableau 1-1 : Composition typique des solutés de la liqueur blancheci-dessous. (Grace &
Malcolm, 1989)
Composition % Total
NaOH 53
Na2S 21
Na2SO4 5
Na2CO3 15
Na2S2O3 3
Na2SO3 3
3
L’objectif de cette cuisson est de dissoudre la lignine qui lie les fibres de cellulose ensemble. (EPA,
1983)
Il existe deux types de système de cuisson, batch et continu. La plupart des usines kraft employaient
le procédé batch. Le premier digesteur en continu à l’échelle industrielle a été mis en service la
première fois en Suède en 1948. A l’arrivée des années 2000, 65% des usines de pâtes et papiers
se sont converties pour ce type de digesteur. (Savelin & Gavle, 2014) Les digesteurs en continu
offrent plusieurs avantages : (NRCAN, 2008) (Savelin & Gavle, 2014)
Quand la cuisson est complétée, le contenu du digesteur est transféré dans un réservoir de collecte
appelé « blow tank » qui fonctionne à pression atmosphérique.
Par la suite, la pâte est séparée de la liqueur de cuisson pendant des étapes de lavage en série. La
pâte subit par la suite un processus de blanchissement avant d’être pressée et séchée en produit fini.
La liqueur récupérée maintenant appelée « Liqueur noire » est concentrée dans des multiples
évaporateurs afin d’atteindre un pourcentage de 65% en solide. Cette liqueur est brûlée dans une
chaudière de récupération afin de produire la vapeur nécessaire pour l’usine. (EPA, 1983)
Une schématisation simplifiée en blocs du procédé kraft est présentée dans la Figure 1-1 ci-
dessous.
4
Copeaux de bois
Vapeur pour le
procédé
La lignine est le deuxième polymère aromatique le plus abondant dans la nature après la
cellulose. La quantité totale de la lignine disponible dans la biosphère dépasse 300 milliards de
tonnes avec une capacité de régénération annuelle d'environ 20 milliards de tonnes. (Smolarski,
2012)
La lignine possède un fort pouvoir calorifique (26 MJ/Kg). D’un autre côté, la chaudière de
récupération est souvent limitée par sa charge thermique qui est directement liée à la quantité de
lignine dans la liqueur noire. Donc, l’extraction d’une partie de lignine permettra de réduire la
charge thermique dans la chaudière et par conséquent, augmenter la capacité de l’usine.
- Augmenter les revenus de l’usine par la transformation de la lignine en produits à valeur ajoutée.
La lignine est une matière brute très versatile qui pourrait être transformée en divers produits à
valeur ajoutée tels que les fibres de carbone, le charbon activé, des phénols, etc.
5
1.2 Problématique
Plusieurs méthodes d’extraction de la lignine à partir de la liqueur noire sont actuellement
à l’étude. Deux méthodes ont été examinées à Polytechnique Montréal soient, la précipitation acide
(Kannangara M. , 2015) et l’électrodialyse avec une membrane bipolaire (Haddad, 2016). En effet,
la solubilisation de la lignine dans la liqueur noire est directement reliée au pH de la solution. C’est
pourquoi les méthodes d’extraction reposent sur l’acidification de la liqueur noire dans le but de
précipiter la lignine avant de la filtrer.
Le Chapitre 5 présente la synthèse des travaux de recherche effectués dans le cadre de la maîtrise
et amorce la discussion sur la portée des découvertes scientifiques qui en découle. Cette dernière
section ouvre finalement la discussion sur de nouvelles avenues de recherche découlant du projet.
7
La solubilité de la lignine dans la liqueur noire est principalement gouvernée par le pH. En
effet, le pH initial de la liqueur noire est de l’ordre de 13, la lignine commence à précipiter lorsque
le pH est inférieur à 10.9. Donc, les processus qui ont été mis au point afin d’extraire la lignine
reposent sur l’acidification de la liqueur noire. (Nagy et al. , 2010)
2.1.1 La lignine
La structure de la lignine est complexe. Toutefois, ses éléments structurels sont connus.
D’un point de vue chimique, la lignine est un polymère formé de trois unités phénoliques, soient
l’alcool coniférylique, l’alcool sinapylique et l’alcool p-coumarylique. (Brodin, 2009) Ces
groupements sont ionisés à pH alcalin.
Les unités aromatiques respectives de ces alcools sont le p-hydroxyphényl (H), le guaiacyl (G) et
le syringyle (S). Ces unités sont liées ensemble par des liens éther et des liens carbone-carbone
pour former le réseau tridimensionnel de la lignine. (Adler, 1977)
Les propriétés de la lignine sont gouvernées d’une part par ses différentes fonctions chimiques, à
savoir l’hydroxyle, le méthoxyle, le carbonyle et le carboxyle et d’autre part par l’agencement de
ses unités (guaiacyle, syringyle et p-hydroxyphényle) (Sarkanen & Ludwig, 1971)
La lignine extraite de la liqueur noire est appelée la lignine technique. Cette lignine se diffère par
sa structure et ses propriétés physicochimiques qui sont fortement liées à la source à partir de
laquelle la lignine est extraite ainsi qu’à la méthode d’extraction. (Lora & Glasser, 2002) Par
rapport à la lignine native, la lignine kraft contient généralement beaucoup plus de groupements
8
hydroxyles phénoliques et de groupements carboxyles. La lignine semble être une matière première
polyvalente pour de nombreuses applications. Les utilisations potentielles de la lignine ont été
classées dans différents groupes listé ci-dessous,
2- Macromolécules
3- Aromatiques
Dans le premier groupe, la lignine est utilisée comme source de carbone pour la production de gaz
de synthèse produit par la gazéification par exemple. Le deuxième groupe utilise la nature
macromoléculaire de la lignine où cette dernière sera transformée en adhésifs, fibres de carbone,
ou des polymères comme les mousses de polyuréthane. Le troisième groupe utilise des technologies
pour cliver la structure de la lignine en monomères sans sacrifier les cycles aromatiques afin de
produire des monomères aromatiques tels que le benzène, le toluène, le xylène (BTX), le phénol et
la vanilline. (Gosselink, 2011)
Lignine
Produits de gaz de
Hydrocarbures Phénols Produits oxydés Macromolécules
synthèse
Les propriétés de la lignine extraite ont une incidence sur l’application envisagée après l’extraction.
D’une façon générale, les propriétés de la lignine qui sont souvent considérées sont : (Sameni,
Krigstin, Rosa, Leao, & Sain, 2014)
Dans un procédé kraft, la liqueur épuisée de la cuisson des copeaux de bois est séparée de
la pâte et traitée dans un circuit de récupération de produits chimiques. Cette liqueur est connue
sous le nom de la liqueur noire.
La liqueur noire est un résidu liquide hautement visqueux dont le pH frôle le 13, elle est constituée
de plusieurs composés inorganiques qui sont ajoutés lors de la fabrication de la pâte, tels que le
NaOH, le Na2S, le Na2CO3, etc. (Theliander, 2007), et des composés organiques. Les composés
organiques sont les substances ligneuses du bois qui lui attribuent sa couleur noire en milieu
fortement alcalin. (Demirbas, 2008)
La liqueur noire contient entre 15% et 40% massique de solides, parmi lesquelles, la lignine
représente entre 30% à 40%. (Cardoso et al. , 2009)
Les pourcentages massiques typiques des composants organiques de la liqueur noire sont présentés
au Tableau 2-2 ci-dessous. (Theliander, 2007)
Tableau 2-2 - Pourcentage massique des matières organiques dans la liqueur noire
Lignine 30-40
Extractibles 3-5
Acide formique ~5
10
Acide acétique ~3
Méthanol ~1
Après l’évaporation d’une partie de l’eau de la liqueur noire généralement effectuée dans un
évaporateur à effets multiples, les liqueurs de cuisson sont régénérés et les substances organiques
sont brûlées dans une chaudière dite de récupération afin de produire de l’énergie nécessaire pour
l’usine. (Theliander, 2007)
Souvent, cette chaudière correspond à un point d’engorgement dans la production d’une usine de
kraft. Pour augmenter la capacité de l’usine sans impliquer des investissements majeurs pour
augmenter la capacité de la chaudière, une extraction d’une quantité partielle de la lignine à partir
de la liqueur noire est envisageable. Cette lignine sera utilisée par la suite comme matière première
pour fabriquer des produits à valeur ajoutée tel que discuté précédemment.
Le principe de précipitation de la lignine par l’acide sulfurique repose sur le procédé inverse
du procédé kraft qui dissout la lignine dans une solution de NaOH. En effet, le pH de la liqueur
noire est abaissé, ce qui entraine une protonation de la lignine et sa précipitation en conséquence.
Le rendement de la précipitation acide est élevé puisque le pH peut être réduit à des valeurs très
basses. Toutefois, cette méthode génère un filtrat acide qui engendre des problèmes d’opération
pour son traitement subséquent. (Ohman & Theliander, 2007)
11
Le procédé consiste à utiliser le dioxyde de carbone, un gaz faiblement acide, qui sera
transformé en acide carbonique avant de protoner la lignine. Cette étape est suivie par des étapes
de coagulation, de filtration et de lavage. (Diop, 2014)
La précipitation par CO2 s’avère plus avantageuse que la précipitation acide puisqu’elle offre des
coûts d’opération beaucoup moins élevés et génère un précipité plus facilement filtrable par rapport
à l’acidification avec l’acide sulfurique. (Kannangara, Marinova, Fradette, & Paris, 2015) La
limitation de son rendement correspond à son désavantage majeur puisque le pH ne peut pas être
abaissé à moins que 9. (Loutfi et al. , 1991)
Au cours du processus d’acidification, que ce soit avec l’acide sulfurique ou le CO2, de nombreux
chercheurs ont signalé une résistance élevée à la filtration, ce qui affecte la pureté de la lignine
ainsi que le coût et la productivité du processus.
[Link] L’ultrafiltration
L’électrodialyse pourrait être définie comme étant un procédé membranaire, avec lequel les
ions sont transportés par une membrane semi-perméable sous l’influence d’un potentiel électrique.
(Lenntech, 2017)
Le Tableau 2-3 ci-dessous sommaire les avantages et les désavantages de chaque méthode décrite
ci-haut.
• Filtrabilité difficile
[Link] LignoBoost
Dans le procédé LignoBoost, un courant de liqueur noire est extrait du courant principal à l’étape
de l’évaporation. Ce courant est acidifié par l’injection du CO2 à une température de 60-80 °C
jusqu’à l’atteinte d’un pH de 10-10.5. L’acidification amène à la précipitation de la lignine. Cette
lignine est par la suite filtrée en utilisant un filtre presse. Le procédé LignoBoost diffère du procédé
d’acidification conventionnel à une étape. En effet, durant le processus d’acidification
conventionnel, la lignine est directement lavée après la filtration, alors que durant le processus
LignoBoost, la lignine filtrée est redispersée dans l’eau de lavage acidifiée au lieu d’être
directement lavée. La boue résultante est de nouveau déshydratée et lavée avec de l'eau de lavage
acidifiée pour produire des gâteaux de lignine pratiquement pur. Cette dispersion supplémentaire
permet le contrôle du pH jusqu’à un pH semblable à celui de la solution du lavage. Ce qui aboutit
à une concentration minimale en ions hydrogène lors du lavage. (Ohman & Theliander, 2007). Les
valeurs exactes/optimales du procédé restent inconnues pour des raisons de confidentialité.
Domtar a démarré avec succès une usine de LignoBoost à l'échelle commerciale dans son usine de
Plymouth en Caroline du Nord, aux États-Unis. C'est la première installation commerciale d'une
usine LignoBoost dans le monde. (Metso, 2013) Une autre usine de LignoBoost a été intégrée par
Stora Enso en Finlande en 2015 dans son usine Sunila. Cette usine possède une capacité annuelle
de production de 370 000 tonnes en pâtes softwood. (Valmet, 2015).
14
[Link] LignoForce
En effet, la liqueur noire est oxydée dans des conditions contrôlées avant l'étape de l’acidification.
L’étape d’oxydation a conduit à une augmentation du taux de filtration de la lignine sans
compromettre la pureté ou la structure de la lignine obtenue. Les avantages du procédé LignoForce
sont : (Lamfeddal, P, Maki, & Paleologou, 2012)
Ce procédé a été intégré dans l’usine de pâtes et papiers de West Fraser en Hinton, Alberta en 2015.
(Tuck & Hunter, 2015)
Au départ, l’électrodialyse a été employée à grande échelle pour le dessalement d’ eaux de mer
pour les rendre aptes à la consommation. Parmi les premières utilisations industrielles figure aussi
la production d’acides et de bases à partir de leur sel. (Gavach, 1999)
Le principe de l’ED est basé sur la sélectivité des membranes pour le transfert de certains
ions sous l’effet d’un champ électrique. Ces transports ioniques vont avoir une incidence directe
sur la nature et la concentration des électrolytes qui composent les milieux traités.
Concentrat (sortie)
Les espèces ioniques qui sont dans la cellule vont migrer en direction de l’électrode de signe opposé
à leur charge sous l’effet du champ électrique. La cellule électrochimique est composée d’une
succession de membranes perméables aux cations (membrane cationique) et perméables aux anions
(membrane anionique) disposées alternativement. Dans le compartiment du diluat, les cations vont
dans le sens du courant électrique traversent la membrane cationique, mais sont bloqués par la
membrane anionique. Ce qui amène à une augmentation de la concentration en cations dans ce
compartiment. De même, les anions vont quitter le diluat en traversant la membrane anionique en
direction de l’anode pour arriver au concentré où ils seront bloqués par une membrane cationique.
(Quoc, 2004)
pour une anionique. Les contre-ions s’échangent avec les électrolytes fixés sur la membrane en la
traversanttandis que les ions de même signe que les groupes ioniques appelés co-ions sont
repoussés par la force électrostatique. C’est ce mécanisme qui détermine la sélectivité des
membranes. (Davis et al. , 1997)
Ces propriétés sont liées à deux paramètres, soient la nature et structure du polymère de base ainsi
que le type et la concentration des charges fixes.
Les membranes bipolaires correspondent à un type spécial de membrane échangeuse d'ions. Elles
sont composées de deux couches de polymère portant des charges fixes, une couche anionique et
une couche cationique. Contrairement aux membranes utilisées pour des fins de séparation, rien ne
devrait être transporté d'un côté à l'autre de la membrane. Le seul transfert désiré est celui des
molécules d’eau en vue d’une réaction dans la jonction bipolaire de la membrane dans laquelle
l’eau est séparée en ions hydroxyles et en protons. Sous l’effet du champ électrique, les protons et
les ions OH- migrent à travers la face cationique et anionique respectivement. (Wilhelm, 2001)
17
Outre les caractéristiques exigées pour les membranes monopolaires, les membranes bipolaires
doivent avoir une perméabilité élevée à l’eau et une vitesse élevée de dissociation de l’eau.
(Hurwitz, 1999)
Un schéma d’une membrane bipolaire est illustré dans la Figure 2-2 ci-dessous.
𝑶 → 𝑯+layer
𝑯𝟐Interphase + 𝑶𝑯−
Couche
Anion selective Couche
Cation selective
anionique
layer cationique
layer
+ -
+ -
-
+
+ -
-
- +
+ -
+ -
OH- + H+
+ -
-
H 2O + H 2O
-
+ + -
-
Anode + - Cathode
+
+ - - -
+
+
+ - -
+
+ -
-
+
+ -
-
+
+ - -
Un séparateur devrait toujours être placé entre deux membranes adjacentes. Les deux
fonctions principales d’un séparateur sont : (Kemperman, 2000)
• Fournir un support structurel aux membranes fragiles et les garder séparées les unes des
autres ;
• Créer des chemins d’écoulement pour les solutions et assurer un certain degré de turbulence
dans le volume entre les séparateurs
L’épaisseur du séparateur définit l’épaisseur des canaux dans lesquelles l’électrolyte va circuler.
18
D’une façon générale, une membrane bipolaire est utilisée avec une membrane cationique,
anionique ou les deux dépendamment des besoins du procédé. Dans le cas de la liqueur noire, une
membrane bipolaire est utilisée avec une membrane cationique. 4 compartiments sont impliqués
soient; le compartiment de la liqueur noire, le compartiment de la soude caustique et deux
compartiments d’électrolyte pour le rinçage des électrodes. Lorsque le courant électrique passe
dans le système, une réaction de dissociation de l’eau prend place dans la jonction de la membrane
bipolaire pour créer des protons H+ et des ions hydroxydes OH-. Les protons migrent vers la cathode
dans le compartiment de la liqueur noire causant ainsi la protonation des groupements phénoliques
ionisés sur la molécule de la lignine tandis que les ions OH- migrent vers l’anode au compartiment
du NaOH pour augmenter la concentration de cette dernière. Le schéma de l’électrodialyse de la
liqueur noire est présenté dans la Figure 2-3 ci-dessous.
La performance de l’électrodialyse est affectée par plusieurs facteurs, tels que la nature
chimique des membranes et leur disposition, les concentrations et la conductivité des solutions
d'alimentation, le courant, la tension (voltage) et la vitesse d'écoulement. (Cloutier et al. , 2013)
La densité de courant est la principale force motrice dans l’ED parce qu'elle détermine la
quantité d'équivalents de sel transporté à travers les membranes. Pour un procédé de séparation
donné, l'utilisation d'une densité de courant élevée nécessite une surface de membrane moins élevée
et un temps de traitement plus court, rendant ainsi le procédé plus attrayant d’un point de vue
économique. Cependant, une densité de courant plus élevée implique une augmentation de la
consommation d'énergie spécifique en raison d'une augmentation de l'irréversibilité
thermodynamique du processus. (Walker, 2010) Il y a donc un compromis entre le coût de l'énergie
électrique et le coût de la membrane. La plupart des équipements d'électrodialyse sont opérés à
voltage constant ou à courant constant ou une combinaison de ces deux modes. En mode de voltage
constant, une tension constante est maintenue à travers la cellule alors que le courant varie avec le
temps, tandis qu'en mode courant constant, une densité de courant est appliquée et le voltage varie
pour maintenir la densité de courant fixe. Typiquement, le voltage augmentera au fil du temps pour
compenser l’augmentation de la résistance afin de maintenir l’ampérage constant.
W.S Walker a publié dans son doctorat une schématisation approximative (Figure 2-4) des profils
de vitesse et de concentration dans un pair d’une cellule d’électrodialyse. (Walker, 2010)
Il est supposé que les séparateurs provoquent un écoulement bien mélangé et que la condition limite
sans glissement à la surface de la membrane provoque des couches limites (symétriques) de
diffusion, d'épaisseur δ, dans laquelle un mélange très faible ou négligeable se tient. Ce concept est
d’une grande importance car ça va permettre d’expliquer les résultats obtenus plus tard dans ce
mémoire.
Les concentrations des solutions alimentées et leur conductivité doivent être maintenues
dans des limites prédéterminées. Toute concentration en dehors de ces limites peut diminuer
l'efficacité du processus ou demander un ajustement des conditions pour maintenir la performance.
Le taux de transport des co-ions dans les membranes idéales serait nul. Mais évidemment,
les membranes ne sont pas rigoureusement imperméables aux co-ions.
Le transport de ces derniers peut-être moins de 1% et il peut aussi atteindre 10% dépendamment
de la concentration en ions dans la solution. (Brun, 1989), (Davis et al. , 1997)
Une partie importante de l’énergie mise en jeu pour l’ED est dissipée sous forme de chaleur
provoquant un échauffement important du produit et du matériel. (Mafart & Beliard, 1992) Cet
échauffement est nuisible, non seulement parce qu’il représente une perte d’énergie, mais aussi
parce qu’il limite la tension applicable jusqu’à une détérioration éventuelle de la membrane. (Brun,
1989)
[Link] La polarisation
La polarisation est un phénomène électrochimique qui limite la migration des ions à travers
les membranes. (Atkins, 1978) Ce phénomène est causé par l’existence d’une couche limite en
écoulement laminaire au voisinage des membranes, et par la différence de mobilité des ions dans
21
[Link] L’encrassement
Ce phénomène pourrait être perçu comme étant une polarisation secondaire. Les causes de
dépôt sont multiples :
L’encrassement peut être vu comme étant la réduction de la surface active de la membrane, ce qui
entraine une diminution du flux au deçà de la capacité théorique de la membrane pour un potentiel
électrique donné. (Field, 2010)
La surface active d’une membrane correspond à ses pores. Donc, les mécanismes
d’encrassement sont principalement liés à celles-ci et les procédures qui entrainent la réduction du
nombre des pores actives. Généralement, quatre mécanismes d’encrassement sont observés comme
illustrés dans la Figure 2-5 ci-dessous. (Field, 2010)
Le Tableau 2-4 ci-dessous résume ces mécanismes et leur effet sur le transfert de masse. (Field,
2010)
23
Mécanisme
Fond phénoménologique Effet sur le transfert de masse
d’encrassement
Blocage complet des pores Les particules sont plus • Réduction de la surface active
(figure 2-4.a) grandes que les pores de membrane
• Dépend de la vitesse de
l’alimentation
Blocage interne des pores Les particules sont de taille • Augmentation de la résistance
(figure 2-4.b) inférieure par rapport aux de la membrane à cause de la
pores et sont adsorbées ou réduction de la taille des pores
déposées sur les parois des
• Indépendant de la vitesse de
pores
l’alimentation
Blocage des pores par les Les particules à la surface • Effet similaire au blocage
particules (figure 2-4.c) pourront bloquer complet, mais moins sévère
partiellement ou combler un
pore
La résistance élevée due à l’encrassement qui survient durant l’opération pourrait être divisée en
trois parties :
• Composante réversible : a lieu lors de l’opération, mais n’est pas présente après la coupure
du courant.
• Composante réversible par lavage : correspond au dépôt de matériel qui est seulement
éliminée par le nettoyage de la membrane.
24
Le débit critique
C’est le débit limite à partir duquel l’encrassement est observé pour une concentration et un
potentiel donné. (Field, 2010)
En plus du débit critique, il est important de considérer s’il y a un seuil de flux à partir duquel
l’encrassement devient élevé. Ceci permet d’examiner si les conditions d’opération pourront être
optimisées pour obtenir un encrassement modéré facilement éliminé par un simple nettoyage
périodique. (Field, 2010)
notion que l’idée d’étudier l’impact de l’hydrodynamique dans la cellule a été proposée. En effet,
nous croyons que la promotion de la turbulence dans la cellule va permettre d’améliorer la
performance de l’opération en retardant l’encrassement de la membrane. Une récente étude nous
a permis de valider cette conviction. Enciso et al. ont récemment mené une étude sur l’impact de
l’hydrodynamique dans une cellule d’électrodialyse à l’aide des simulations numériques de
dynamique des fluides (CFD). Bien que la cellule a été simulée en régime laminaire, ils ont trouvé
que la performance de la cellule a été améliorée lorsque la vitesse du fluide augmente. (Enciso,
Delgadillo, Domínguez, & Rodríguez-Torres, 2017)
Cette étude est d’une grande importance puisque les expériences qui sont réalisées dans le cadre
de ce projet permettront la validation des résultats obtenus en simulation numérique.
25
La perte de charge désigne les pertes irréversibles d’énergie de pression que subit un fluide
lors de son passage dans une géométrie quelconque.
Il n’existe pas une formule explicite qui permette de calculer la perte de charge dans une cellule
électrochimique. Toutefois, il est évident qu’un écoulement dans une cellule va produire une perte
de charge supplémentaire par rapport à un écoulement dans une conduite vide, Il est donc
primordial de pouvoir exprimer la perte de charge dans la cellule.
La perte de charge dans n’importe quelle géométrie est fonction du débit, de la masse volumique
du fluide, de se viscosité, et des dimensions physiques du système exprimés par le biais du rayon
hydraulique et de sa longueur caractéristique.
Où Dh est le diamètre hydraulique qui correspond au rapport entre la surface et le périmètre d’une
conduite et 𝑣 est la vitesse moyenne. La notion du rayon hydraulique est utilisée dans ce cas-là car
elle permet de calculer le nombre de Reynolds même si la géométrie de l’écoulement n’est pas
cylindrique. Ce diamètre correspond au diamètre équivalent d'un tuyau circulaire qui donne le
même bilan de forces que pour un conduit non circulaire. (Bird, Stewart, & Lightfoot, 2006)
𝐷 1 Équation 2-2
𝑁𝑒 = 2𝑓 = ∆𝑃 ( ) ∗ ( )
𝐿 𝜌𝑣 2
Cette équation est d’un grand intérêt étant donné qu’elle fait apparaître une perte de charge
mesurable dans le cas du montage employé. (Arzat, 1999)
Plus loin, une constante représentative de la géométrie de la cellule sera déterminée. Cette constante
permet le calcul d’une vitesse de cisaillement caractéristique pour n’importe quelle condition et
26
Il faut noter que les résultats expérimentaux seront associés aux deux régimes laminaire et
turbulent. Évidemment, il y a un régime transitoire entre les deux mais ce dernier ne sera pas pris
en considération.
Les propriétés physiques et chimiques d’un fluide ont un impact sur son comportement lors
d’un écoulement. Il existe deux types de comportement visqueux, newtonien et non newtonien
dépendamment de la relation entre la viscosité du fluide et le taux de cisaillement. Cette relation
est connue sous le nom de la loi de Newton et pourrait être écrite comme suit :
𝜏 = 𝜇 ∗ 𝛾̇ Équation 2-3
Un fluide newtonien est défini comme étant un fluide dont la viscosité reste constante quelle que
soit l’intensité du taux de cisaillement (𝛾̇ ). Si la viscosité d’un fluide dépend du cisaillement, nous
sommes dans le cas de fluide non newtonien.
1- Les fluides rhéofluidifiants : sont les fluides dont la viscosité diminue lorsque le taux de
cisaillement augmente
27
2- Les fluides rhéoépaississants : sont les fluides dont la viscosité augmente lorsque le taux
de cisaillement augmente
Metzner et Otto ont développé en 1957 une méthode en agitation mécanique qui consiste à
calculer une constante représentative de la géométrie et indépendante des conditions d’écoulement
(Ks) dans le cas d’une cuve agitée avec un liquide rhéofluidifiant en régime laminaire. En effet, Ils
ont mesuré le nombre de puissance pour une variété de mélangeurs en régime laminaire pour des
fluides newtoniens, puis ils ont répété les mêmes mesures avec des fluides rhéofluidifiants. Ils ont
trouvé que le nombre de puissance n'a pas été affecté par le comportement non-newtonien du fluide
et que les viscosités newtonienne et apparente étaient égales pour le même nombre de puissance et
le même nombre de Reynolds. (Paul, Atiemo-Obeng, & Kresta, 2003)
Selon Metzner et Otto, le mouvement du fluide dans une cuve peut être représenté par une vitesse
de cisaillement moyenne ou effective. Cette vitesse est reliée à la vitesse de rotation de l’agitateur
(N) par la relation suivante:
Depuis, la notion a été étendue et plusieurs études suggéraient que la notion de Ks peut être utilisé
pour d’autre type de mélangeurs tels que les mélangeurs statiques. (Fradette, 1999), (Arzat, 1999),
(Paul, Atiemo-Obeng, & Kresta, 2003)
2. Calculer Ne pour un fluide non newtonien par l’entremise des mesures de couple obtenus
(la perte de charge en fonction du débit)
𝐷2 ∗ 𝑁 ∗ 𝜌
𝜇𝑒𝑓𝑓 = Équation 2-5
𝑅𝑒𝑒𝑓𝑓
28
Cette viscosité correspond à celle d’un fluide newtonien qui exigerait la même puissance
d’agitation aux mêmes conditions.
À nos connaissances, il n’existe aucun travail expérimental relié à l’étude de l’impact des
paramètres hydrodynamiques sur l’extraction de la lignine à partir de la liqueur noire dans une
cellule d’électrodialyse.
À l’issue de cette analyse, les objectifs spécifiques du projet sont définis comme suit :
La réalisation de ces trois objectifs nous permettra d’atteindre l’objectif principal du projet.
29
Afin d’atteindre les objectifs du présent projet de recherche, le travail a été divisé en trois phases :
• Le système d’acquisition des données : qui correspond à l’interface entre les capteurs et
l’ordinateur
Une cellule électrochimique est composée d’une succession d’étages disposés d’une façon
particulière afin de laisser circuler les différentes solutions dans la cellule sans qu’elles soient
mélangées.
Dimensions en mm
Élément (longueur X largeur X Matériau
épaisseur)
Suivant les dimensions fournies au tableau précédent, il est possible de calculer la surface active.
Cette surface correspond à la surface des membranes où la réaction va avoir lieu :
𝑆𝑢𝑟𝑓𝑎𝑐𝑒 𝑎𝑐𝑡𝑖𝑣𝑒 = 0.0635 𝑚 ∗ 0.0635 𝑚 = 0.00403 𝑚2
La conception du séparateur a été modifiée plusieurs fois afin d’éliminer toute fuite à l’intérieur de
la cellule. La conception optimale est illustrée dans la Figure 3-2 ci-dessous.
Afin de faire varier le degré de turbulence à l’intérieur de la cellule, trois géométries sont
proposées :
La seule différence entre les deux grille #1 et #2 est l’ouverture de leurs mailles qui sont de l’ordre
de 1.013 mm et 1.778 mm respectivement. Ces grilles sont présentées dans les figures ci-dessous.
Il faut noter que bien que l’ouverture des mailles de la grille #1 est plus petite que celle de la grille
#2, la grille #1 présente moins de turbulence par rapport à la grille #2 pour deux raisons :
- La grille #2 est plus épaisse (2.38 mm) que la grille #1 (0.9 mm) occupant ainsi tout
l’espace disponible dans la cavité du séparateur contrairement à la grille #1
Grâce à la conception proposée, les trois solutions circulent à l’intérieur de la cellule sans être
mélangées. Lorsque le courant électrique passe dans le système, la réaction de la dissociation d’eau
prend lieu et les ions migrent tel que discuté dans la section [Link]. La migration des ions permet
l’acidification du diluat et l’alcalinisation du concentré.
À noter que toutes les composantes de la cellule ont été commandées en feuilles qui ont été
découpées par la suite à l’aide d’une découpeuse laser Trotec Speedy 300.
34
pompes et la cellule)
pression différentiel
• Plage de mesure 0-500
kPa
• Précision de 2.5%
L’acquisition des données permet le suivi de mesures de variables opératoires tel que la tension, le
courant, la température, la pression, le débit, la conductivité et le pH avec un ordinateur. Un
système d’acquisition des données comprend des capteurs, du matériel de mesure et un ordinateur
avec un logiciel programmable. Le système de l’acquisition des données ainsi que la
programmation derrière, a été conçu par M. Daniel Pilon de l’équipe technique du département de
génie chimique à l’école Polytechnique.
Les valeurs du voltage (V), courant (Amp), débits (l/min), perte de charge (Pa), conductivité (s),
pH et température sont enregistrées dans un fichier Excel à toutes les secondes.
Un étage correspond à un élément dans la cellule électrochimique. Cet élément pourrait être un
séparateur, une membrane ou un joint.
Afin de faire la caractérisation de la cellule conçue, trois solutions newtoniennes ont été utilisés.
Les détails de ces solutions sont présentés dans le tableau ci-dessous.
Le glucose 30% a été choisi car il possède la même viscosité que la liqueur noire avant le début de
l’électroacidification (~0.003 Pa.s). À noter que la viscosité de toutes les solutions incluant la
liqueur noire a été testée avec un rhéomètre Physica MCR 501. L’eau et le glucose 42% ont été
choisis comme des valeurs autour du glucose 30%.
Afin de bâtir les courbes de puissance pour les différentes géométries, les étapes suivantes ont été
suivies :
a) Mesurer la perte de charge dans la cellule pour des débits entre 0.9 l/min et 2.8 l/min avec
un pas de 0.1 l/min pour chaque solution et pour trois géométries différentes tel que proposé
auparavant :
𝑄
𝑣 = Équation 3-1
𝐴𝑖𝑟𝑒 𝑡𝑟𝑎𝑛𝑠𝑣𝑒𝑟𝑎𝑙𝑒
𝑆𝑆𝐸
𝑅2 = 1 − Équation 3-3
𝑆𝑆𝑇
Avec
newtonien utilisé dans le cadre de cette étude est le carboxyméthylcellulose (CMC) dilué dans l’eau
avec une concentration de 1% . Il s’agit d’un gel rhéofluidifiant qui suit une loi de puissance. Ce
fluide a été choisi car il a le même comportement rhéologique que la liqueur noire qu’on désire
avoir à la fin de l’expérience (10.2 < pH < 10.3)
Les propriétés rhéologiques du CMC 1% ont été mesurés à l’aide d’un rhéomètre Physica MCR
501 Le CMC 1% possède les propriétés suivantes.
Masse volumique k n
1. Établir les courbes de puissances (Ne vs. Re) à l’aide des fluides newtoniens présentés dans
le Tableau 3-2 - Caractéristiques des solutions utilisées pour la caractérisation
2. Effectuer des mesures de couple débit/perte de charge pour le fluide non newtonien (CMC
1%)
𝜌 ∗ 𝐷ℎ ∗ 𝑣
𝜇𝑒𝑓𝑓 = Équation 3-7
𝑅𝑒𝑒𝑓𝑓
1
𝐷ℎ 𝜇𝑒𝑓𝑓 𝑛−1 Équation 3-8
𝐾𝑠 = ( ) ∗ ( )
𝑣 𝑘
À partir de la constante Ks calculée pour chaque géométrie, la viscosité effective pour n’importe
quel fluide non newtonien est obtenue par l’entremise de la courbe d’écoulement de ce fluide selon
l’équation.
𝑄
𝛾̇ 𝑒𝑓𝑓 = 𝐾𝑠 ∗ ( ) Équation 3-9
𝑣
Il est important de caractériser la liqueur noire avant de procéder avec les expériences. Chaque
élément caractérisé a un impact sur le comportement de cette solution.
Le total des solides dissous a été déterminé selon la méthode TAPPI T650 om-09. Un
échantillon de liqueur noire a été séché à une température de 105°C pour 24 h. Par la suite,
l’échantillon a été pesé et le total des solides dissous a été calculé selon l’équation :
𝑚𝑠𝑒𝑐 − 𝑚𝑎𝑙𝑢𝑚
%𝑇𝑆𝐷 = ∗ 100 Équation 3-10
𝑚𝐿𝑁 − 𝑚𝑎𝑙𝑢𝑚
En effet, le % de solide dans la liqueur noire a un incident direct sur la viscosité et la conductivité
de la liqueur noire. L’objectif est d’avoir un fluide newtonien au départ avec une conductivité
optimale.
40
La masse volumique de la liqueur noire a été déterminée à l’aide d’une fiole jaugée. La
fiole jaugée a été rempli avec un volume précis de liqueur noire. Par la suite, ce volume a été pesé
à l’aide d’une balance. La masse volumique a été calculée par la suite selon :
𝑝𝑜𝑖𝑑𝑠
𝜌𝐿𝑁 = Équation 3-11
𝑣𝑜𝑙𝑢𝑚𝑒
La lignine est l'un des principaux constituants organiques de la liqueur noire. Il contient des
groupes phénoliques qui absorbent la lumière ultraviolette à une longueur d'onde spécifique de
sorte qu'un spectrophotomètre UV peut être utilisé pour sa détection.
Un échantillon de liqueur noire doit être dilué à 15% TSD et ajusté à un pH de 13. Puis 1 g de
l'échantillon est transféré dans une fiole jaugée de 50 ml et rempli jusqu'à la marque avec 0.1 M de
NaOH pour faire une solution originale. 1 ml et 2 ml de cette solution d'origine sont transférés dans
deux flacons volumétriques de 50 ml et remplis jusqu'à la marque avec de l'eau déminéralisée pour
faire deux solutions de concentration différente (solution 1 et 2). Pour la détermination par UV,
des cuves de verre de quartz spéciales doivent être utilisées sinon la détection n'est pas précise en
raison de l’interférence avec des matériaux de cuvette tels que le plastique.
Les deux échantillons sont analysés et la concentration peut être déterminée en utilisant la loi de
Lambert-Beer.
𝐼0
𝐴 = 𝑙𝑜𝑔10 ( ) = 𝜀 ∙ 𝑐 ∙ 𝑑 Équation 3-12
𝐼
La concentration en lignine peut être calculée en se basant sur cette formule avec un coefficient
d'absorption de 23,7 g / (l ∙ cm) et un diamètre de la cuvette de 1 cm. À partir de là, l’équation ci-
dessous est utilisée.
41
𝑔 𝐴
𝑐[ ] = ∙𝐷 Équation 3-13
𝑙 23.7
𝑉𝐷
𝐷= Équation 3-14
𝑉0
Cependant, plusieurs études ont découvert que la teneur en lignine peut être surestimée en raison
des impuretés présentes telles que les produits de dégradation des hydrates de carbone. En se basant
sur des expériences avec la lignine pure et la détermination de la concentration des impuretés, des
facteurs de correction ont été définis de sorte qu'un résultat plus précis peut être calculé par les
équations suivantes.
𝐴1
𝑋1 (𝑐𝑜𝑟𝑟. ) = 10.54 Équation 3-15
𝑆
𝐴2
𝑋2 (𝑐𝑜𝑟𝑟. ) = 5.27 Équation 3-16
𝑆
Où A1 et A2 sont les valeurs d'absorbance des solutions 1 et 2 à 280 nm, S est le % du solide dans
la solution originale et X (corr.) est le pourcentage de lignine dans la liqueur noire.
Une caractérisation rhéologique de la liqueur noire s’avère une étape indispensable pour pouvoir
comprendre le comportement hydrodynamique de cette solution et l’évolution de ce comportement
avec la diminution du pH.
Afin de faire la caractérisation, 14 échantillons de solution avec un total de solide dissout de 20%
ont été préparé. Chacun de ces échantillons a été acidifié jusqu’à un certain pH en utilisant une
solution de l’acide sulfurique 1N. Puis la viscosité de chaque échantillon a été mesurée à l’aide
d’un rhéomètre Physica MCR 501 avec un taux de cisaillement allant de 1 à 1500 1/s avec une
rampe logarithmique.
Les pH des différents échantillons préparés sont présentés dans le tableau ci-dessous.
42
E1 E2 E3 E4 E5 E6 E7
Trois solutions sont utilisées pour faire l’électrodialyse de la liqueur noire soient, la liqueur
noire, une solution de NaOH et la solution de rinçage des électrodes.
La première étape pour préparer la liqueur noire est la dilution. La liqueur noire lourde après
les étapes d’évaporation a un pourcentage de solide dissout entre 50% et 60%. Cette liqueur est
diluée jusqu’à un pourcentage de solide dissout de 20% en ajoutant l’eau déionisée selon
l’équation :
𝐶1 ∗ 𝑉1 = 𝐶2 ∗ 𝑉2 Équation 3-18
C’est V1 (volume de liqueur noire brut) qui doit être déterminée dans cette équation, C1 étant la
concentration de la liqueur noire au départ, C2 étant 20% et V2 est le volume final désiré.
Après la dilution, la liqueur noire diluée est filtrée sous vide avec un filtre numéro 113, diamètre
150mm afin d’enlever les solides non-dissouts, le savon, etc.
Une fois la filtration terminée, la liqueur noire est prête pour être utilisée dans l’électrodialyse.
La solution basique de NaOH agit comme le concentrat dans l’électrodialyse. Cette solution
est préparée par la dissolution des granules de la soude caustique dans de l’eau distillée afin
d’obtenir une concentration désirée de 5%.
43
Afin de dissoudre les granules de Na2SO4 dans de l’eau distillée, un agitateur TBD est utilisé
tournant à une vitesse de rotation N = 22.5 rpm
Une fois les trois solutions sont préparées, il faut suivre le protocole suivant pour effectuer
l’électrodialyse de la liqueur noire :
1- Monter les étages de la cellule d’électrodialyse et connecter les conduits entre la cellule et
le montage
2- Démarrer labview
6- Faire circuler de l’eau dans le système pour purger les bulles d’air qui peuvent fausser les
résultats et pour s’assurer qu’il n’y a pas de fuite
9- Mettre en marche le courant électrique qui restera constant tout au long de l’opération
10- Monitorer et enregistrer les résultats obtenus pour les différents paramètres
Deux débits d’opération ont été choisis pour chaque géométrie. Un débit qui correspond au régime
laminaire et un autre pour le régime turbulent.
Il faut noter que la résistance électrique a été choisie comme paramètre observé pour déterminer
l’impact de l’hydrodynamique dans la cellule sur l’encrassement de la membrane qui va se
construire au fur et à mesure que la réaction avance. La résistance électrique pourrait être exprimée
selon l’équation suivante.
44
𝑉
𝑅= Équation 3-19
𝐼
Puisque la cellule est opérée à courant constant, lorsque la lignine commence à précipiter, des
dépôts de particules vont commencer à s’accumuler sur la membrane. Ceci devrait causer une
augmentation dans le voltage afin de maintenir un courant constant. Cette augmentation de voltage
cause une augmentation de résistance selon la réaction précédente.
À la fin de l’opération, les cuves sont vidées dans des barils de collecte des résidus et de l’eau
propre est circulée plusieurs fois afin de bien nettoyer le système.
45
CHAPITRE 4 RÉSULTATS
Ce chapitra présente les résultats au même ordre de la méthodologie présentée au troisième chapitre
À noter que les calculs de l’erreur sont expliqués en détails dans l’Annexe B.
Tel que discuté dans la section 3.2.1, trois fluides newtoniens de références ont été choisis pour
calculer les paramètres rhéologiques soient, l’eau, glucose 30% et glucose 42%
Les pertes de charges en fonction du débit obtenues pour ces fluides sont présentées à la figure 4-
1 ci-dessous
100000
Perte de charge (Pa)
10000
ΔP - Eau
ΔP - glucose 30
1000 ΔP - Glucose 42
100
0.8
Débit (l/min)
100
Ne
10
1
25 250 2500
Re
Géométrie Relation R2
4.1.2 Détermination de Ks
1000
Ks
100
10
10 100 1000
Reeff
En calculant la moyenne des valeurs de Ks obtenues, il a été trouvé que la constante pour la
géométrie I : Ks = 275
48
La perte de charge en fonction du débit pour la géométrie II est illustrée dans la Figure 4-4 ci-
dessous.
100000
Perte de charge (Pa)
10000 ΔP - Eau
ΔP - glucose 30
ΔP - Glucose 42
1000
0.8
Débit (l/min)
100
Ne
10
1
25 250 2500
Re
Géométrie Relation R2
[Link] Détermination de Ks
La figure ci-dessous présente les valeurs de Ks en fonction de Reeff pour la géométrie II.
50
1000
Ks
100
10
10 100 1000
Reeff
En calculant la moyenne des valeurs de Ks obtenues, il a été déterminé que la constante pour la
géométrie II : Ks = 308
La perte de charge en fonction du débit pour la géométrie II est illustrée dans la Figure 4-7 ci-
dessous.
51
100000
Perte de charge (Pa)
10000 ΔP - Eau
ΔP - glucose 30
ΔP - Glucose 42
1000
0.8
Débit (l/min)
100
Ne
10
1
25 250 2500
Re
Géométrie Relation R2
[Link] Détermination de Ks
La figure ci-dessous présente les valeurs de Ks en fonction de Reeff pour la géométrie III.
53
10000
1000
Ks
100
10
10 100 1000
Reeff
En calculant la moyenne des valeurs de Ks obtenues, il a été déterminé que la constante pour la
géométrie III : Ks = 551
Depuis les résultats obtenus dans les sections précédentes, il est possible de retirer les conclusions
suivantes :
Conclusion 1 : La perte de charge dans la cellule pour géométrie III est plus élevée que pour
géométrie II qui est de son tour plus élevée que pour géométrie I.
Ce résultat est tout à fait sensé, puisque pour chaque géométrie on ajoute un élément qui perturbe
l’écoulement de façon plus importante. Cet élément va générer plus de friction pour la géométrie
III que pour la géométrie II. La géométrie I n’a pas de grille, c’est pourquoi elle présente la perte
de charge la moins élevée.
Conclusion 2 : Le nombre de Reynolds critique pour géométrie III est plus grand que pour
géométrie II. Et ce dernier est plus grand pour la géométrie I.
54
Cette augmentation dans le nombre de Reynolds critique pourrait être induite par le fait qu’on
n’impose pas un facteur correctif relié à la taille des mailles pour chaque grille pour les géométries
II et III, ce qui cause un nombre de Reynolds différent pour chaque géométrie.
Conclusion 3 : Pour le même nombre de Reynolds, le nombre de Newton est plus élevé pour la
géométrie qui présente plus de turbulence. Cette conclusion est directement reliée à la relation entre
les deux nombres adimensionnels. Le tableau ci-dessous sommaire les relations Ne vs. Re obtenues
pour chaque géométrie.
Tableau 4-4 - Sommaire des relations Ne-Re pour les trois géométries
Komar Kawatra mentionne dans son livre « Advances in Comminution » que le nombre de Newton
augmente lorsque la turbulence est plus élevée pour le même nombre de Reynolds. (Kawatra, 2006)
Cette déclaration concorde avec les résultats que nous avons obtenus
Nous remarquons depuis le tableau précédent aussi que la constante Kp augmente avec
l’augmentation de la turbulence dans le système. Ce résultat est sensé puisque Kp dicte la relation
entre Ne et Re et l’augmentation de Kp aboutit à l’augmentation de Ne lorsque Reynolds augmente.
D’un autre côté, KNE ne semble pas avoir une tendance particulière telles que les autres constantes.
Conclusion 4 : Ks suit la même tendance que Reynolds critique. Autrement dit, Ks pour la
géométrie III est plus élevée que pour la géométrie II. Et Ks pour la géométrie II est plus élevée
que pour la géométrie I. Le tableau ci-dessous sommaire les valeurs de Ks pour les trois géométries.
55
Géométrie Ks
Grille #1 308
Grille #1 551
Dans le but d’illustrer et de fournir une représentation de valeurs typiques de Ks, le tableau suivant
présente quelques valeurs de référence pour des agitateurs mécaniques souvent retrouvés dans
l’industrie.
Mélangeur Ks Référence
En comparant les valeurs de Ks obtenues avec les valeurs présentées dans le tableau précédent, il
est clair que les valeurs de Ks obtenues pour les trois géométries sont beaucoup plus élevées par
rapport aux mélangeurs qu’on retrouve dans l’industrie d’un facteur d’au moins cinq par rapport à
la valeur la plus élevée qui provient d’un mélangeur statique, la seule géométrie jamais étudiée qui
ressemble un peu aux grilles utilisées ici.
Par opposition aux mélangeurs statiques conçus pour favoriser le mélange sans trop augmenter la
perte de charge, les grilles engendrent une très forte perturbation de l’écoulement, bien supérieure
56
à celles du SMX car elles ont été choisies pour favoriser le nettoyage de la surface des membranes
et non pour engendrer du mélange au coût le plus bas possible.
Un lot de liqueur noire faible a été reçu d’une usine canadienne de pâte et papier. En suivant le
protocole décrit dans la section [Link], il a été trouvé que le total des solides dissous dans la liqueur
noire est de l’ordre de 20%. Aucune dilution n’est nécessaire dans ce cas-là puisqu’il a été décidé
d’utiliser un pourcentage de solide de 20% pour effectuer les expériences d’électrodialyse.
La masse volumique et le teneur en lignine ont été déterminés selon les sections [Link] et [Link]
Selon les résultats de mesure de viscosité obtenus pour les 14 échantillons il a été trouvé que :
Ce résultat est plausible puisque la lignine commence à précipiter à partir d’un pH de 10.9. À
un pH de 10.5, la lignine précipitée va commencer à avoir une incidence sur le comportement
rhéologique de la liqueur noire
- Bien que le comportement de la liqueur noire reste newtonien entre les pH 13.14 et 10.5, la
viscosité augmente de 0.003 jusqu’à 0.008 entre ces deux pH.
57
1
10 100 1000 10000
La viscosité (Pa.s)
0.1 pH = 10.67
pH = 10.29
pH = 9.85
pH = 9.52
0.01
pH = 9.16
0.001
Taux de cisaillement (1/s)
Figure 4-10- La viscosité en fonction du taux de cisaillement pour différents pH de la liqueur noire
Courant 1.2 A
Les débits d’opération pour chaque géométrie ont été choisis comme suit :
Il faut noter que pour la géométrie #3 il n’a pas été possible d’opérer le système à un débit turbulent
(3.8 l/min) à cause que les pompes sont saturées à cause de la perte de charge induite à ce débit-là.
Il est important de noter que le débit a été calculé en prenant en compte une viscosité initiale de la
liqueur noire égale à 0.003 Pa.s qui correspond à la viscosité de la liqueur noire avant le début de
l’acidification.
Les résultats de l’électrodialyse pour les différentes situations sont illustrés dans la figure ci-
dessous sous la forme de la résistance électrique en fonction du temps.
59
Résistance (Ω) 5
4
Géométrie I - Laminaire
3 Géométrie II - Laminaire
Géométrie III - Laminaire
2
Géométrie I - Turbulent
1 Géométrie II - Turbulent
0
0 20 40 60 80
Temps (min)
Figure 4-11 - Évolution de la résistance électrique pour les différentes géométries avec
régimes laminaire et turbulent
Ces résultats montrent que l’impact de l’hydrodynamique est très faible voir inaperçu en régime
laminaire puisque la résistance suit sensiblement la même tendance pour les trois géométries.
Toutefois, la cellule semble présenter une meilleure performance en régime turbulent en
comparaison avec le régime laminaire puisque la résistance augmente d’une façon moins
importante qu’en régime laminaire. La performance de la cellule s’est grandement détériorée à
partir d’une résistance approchant le 3.5 Ohm. En comparant les deux régimes laminaire et
turbulent, il est notable qu’en régime laminaire cette résistance a été atteinte après 60 minutes
environs alors qu’en régime turbulent, la même valeur a été atteinte au bout de 75 minutes environ.
Donc, la durée d’opération est 25% plus longue dans le cas du régime turbulent par rapport au
régime laminaire. Ce temps plus long implique moins d’arrêt pour effectuer le nettoyage, ce qui
permet de réduire les coûts d’opération.
Il est notable aussi depuis la figure que l’impact des différentes géométries semble être plus visible
en régime turbulent qu’en régime laminaire. Ceci peut être expliqué par la notion des couches
limites de diffusion dans les procédés membranaires. En effet, pendant l’électrodialyse, une couche
limite se forme près de la membrane avec un mélange très faible ou négligeable. L’épaisseur de
60
cette couche est fonction de la turbulence près de la membrane. Plus la turbulence est importante,
moins le serait l’épaisseur de la couche limite et vice-versa. Dans le cas de la présente étude, la
grille #2 qui génère plus de turbulence présente une meilleure performance par rapport à la grille
#1 et la géométrie sans grille surtout avec l’avancement de la réaction.
Le travail sur l’électrodialyse n’a pas été élaboré d’avantage à cause de la contrainte de temps.
L’objectif de cette partie était strictement de voir si l’hydrodynamique a ou non une incidence sur
l’encrassement.
61
Une nouvelle technologie qui permet d’extraire la lignine à partir de la liqueur noire a été
examinée dans le cadre de ce projet de recherche. Dans une première étape, une cellule
électrochimique à membrane bipolaire a été conçue et construite. La cellule a été opérée avec 4
compartiments; un compartiment acide, un compartiment basique et deux compartiments
d’électrolyte qui servent à rincer les électrodes. La surface active de la cellule est 0.004 m2. Aussi,
trois configurations de la cavité ont été testées afin de faire varier le niveau de turbulence dans la
cellule. Des éléments de mélange ont été ajouté dans deux cas alors que la cavité vide a été utilisée
pour établir la performance de référence généralement utilisée.
Il a été trouvé aussi qu’en régime laminaire, l’impact de la cavité dans la cellule est moins
important que dans le cas du régime turbulent.
La principale recommandation pour les futurs travaux serait d’optimiser les différents paramètres
hydrodynamiques de l’opération afin d’améliorer la performance de l’électrodialyse de la liqueur
noire. Ce projet de recherche, lorsque discuté en détails, laisse également la place à d’autres
recommandations tels que :
62
- Déterminer l’impact des géométries différentes aux grilles utilisées ici, en régimes
laminaire et turbulent, sur la réduction de l’encrassement
- Ajouter un système de pulsation de débit qui permettra de faire des pulsations de débit sur
certaines plages de temps et examiner l’impact de cette pulsation sur l’encrassement de la
membrane; combiner avec l’utilisation des régimes d’écoulement.
- La cellule électrochimique est comme une boîte noire, ça serait donc intéressant de faire
une étude CFD de ce qu’il se passe à l’intérieur de la cellule afin de mieux comprendre la
circulation des fluides et l’impact des paramètres hydrodynamiques
- Procéder avec les étapes suivant l’électroacidification telles que la maturation et la filtration
de la lignine précipitée pour en évaluer les qualités finales ; comparer avec la qualité de la
lignine déjà publiée et obtenue par acidification directe ou avec du CO2.
- Réaliser des recherches poussées sur les procédés de transformation de la lignine extraite
en produits à valeur ajoutée
- Réaliser une analyse de cycle de vie (ACV) afin de quantifier les impacts
environnementaux de ce procédé ainsi que de valider son efficacité d’un point de vue
environnemental.
63
BIBLIOGRAPHIE
• Adler, E. (1977). Lignin chemistry-past, present and future. 11, 169-218. Wood Science
Technology.
• APEC. (2017). What is ultrafiltration (UF) in terms of membrane filter technology? APEC
WATER.
• Arkall, A., Olsson, J., & Wallberg, O. (2014). Process performance in lignin separation from
softwood black liquor by membrane filtration. Chem Eng Res Des.
• Arzat, A. (1999). Détermination de la viscosité en ligne des sauces de couchage de papier. École
Polytechnique de Montréal.
• Atkins, P. (1978). Physical chemistry. San Fransisco: Freeman, W.H and Co.
• Bird, R. B., Stewart, W., & Lightfoot, E. (2006). Transport Phenomena. Wiley.
• Brodin, I. (2009). Chemical properties and thermal behaviour of kraft lignins. Stockholm,
Sweden: KTH Royal Institute of Technology.
• Cloutier, J., Champagne, P., & Labrecque, R. (2013). Electrolytic Treatment of Black Liquor and
Caustic Recovery. United States of America: PEERS conference.
• Da Re, L., & Papinutti, L. (2011). Black liquor decolorization by selected white-rot fungi.
Applied Biochemistry and Biotechnology, 165, 406.
• Demirbas, A. (2008). Reconvery of oily products from organic fraction of black liquor via
pyrolysis. Energy Source, 30, 1849.
• Enciso, R., Delgadillo, J., Domínguez, O., & Rodríguez-Torres, A. (2017). Analysis and
validation of the hydrodynamics of an electrodialysis cell using computational fluid dynamics.
ELSEVIER.
• EPA. (1983, September). Review of New Source Performance Standards for Kraft Pulp Mills.
EPA - 450/3-83-017, U.S Environmental Protection Agency.
• Field, R. (2010). Fundamentals of Fouling. Dans K. Peinemann, S. K. Pereira Nunes, & S. Pereira
Nunes (Éd.), Membrane Technology: Membranes for Water Treatments (Vol. 4). Weinheim,
Germany: Wiley-VCH Verlag GmbH & Co. KGaA.
• Fradette, L. (1999). Étude de dispersion dans un mélangeur statique Sulzer SMX. QC., Canada:
École Polytechnique de Montréal.
• Gosselink, R. (2011). Lignin as a renewable aromatic : resource for the chemical industry.
Wageningen University.
• Grace, T., & Malcolm, E. (1989). Pulp and Paper manufacture : Alkaline pulping 3rd ed, . TAPPI
Press.
• Granier, P., & Bibenet, J. (1976). Les procédés de séparation par membrane: électrodialyse,
osmose inverse, ulrafiltration. Série synthèse bibliographiques, Centre de documentation
internationale des industries utilisatrices de produits agricoles.
• Holmqvist, A., Wallberg, O., & Jonsson, A.-S. (2005). Ultrafiltration of Kraft Black Liquor from
Two Swedish Pulp Mills. ELSEVIER Volume 83, Issue 8.
• Hurwitz, H. (1999, Mars 23-24). Principes de l'lectrodialyse avec membranes bipolaires. Saint-
Hyacinthe: Symposium Applications alimentaires de l'électrodialyse.
• INERIS. (2017). Processus de fabrication de la pâte kraft (pâte au sulfate). Récupéré sur Ineris:
[Link]
65
• Kannangara, M., Marinova, M., Fradette, L., & Paris, J. (2015). Effect of Mixing Hydrodynamics
on The Particle And Filtration Properties of Precipated Lignin . Canada: Elsevier.
• Kennan, G., Parkinson, D., & Jang, B. (2014). Paper trail: The decline of Canada’s forestry
industry.
• Kuhad, R., & Singh, A. (2007). Lignocellulose biotechnology: future prospects. I.K.
International Pvt Ltd. .
• Lamarche, F., & Boulet, M. (1992). Électrodialyse et alimentation: le courant passe! 5, 14-15.
Alimentech.
• Lamfeddal, K., P, H.-H., Maki, K., & Paleologou, M. (2012). The LignoForce system : A new
process for the production of high-quality lignin from black liquor. Journal of Science &
Technology for Forest Products and Processes VOL.2, NO.4.
• Lora, J., & Glasser, W. (2002). Recent industrial applications of lignin: a sustainable alternative
to nonrenewable materials. 10(39-48). J. Polym. Environ.
• Mafart, P., & Beliard, E. (1992). L'électrodialyse. Dans Génie industriel alimentaire Tome 2.
Techniques et Documentation (pp. 185-208). Paris: Lavoisier.
• Metso. (2013). Metso-supplied world's first commercial LignoBoost plant successfully starts up
at Domtar in the USA. Metso.
66
• Minge, Z., Luhong, Z., Bin, J., & Yuguo, Y. (2008). Calculation of Metzner Constant for Double
Helical Ribbon Impeller by Computational Fluid Dynamics Method. Chinese Journal of
Chemical Engineering.
• NRCAN. (2008). Benchmarking energy use in Canadian pulp and paper mills. Camada: Natural
Resources Canada’s Office of Energy Efficiency.
• Ohman, F., & Theliander, H. (2007). Filtration properties of lignin precipitated from black liquor
. TAPPI Journal, 6, 3.
• Öhman, F., Wallmo, H., & Theliander, H. (2007). An improved method for washing lignin
precipitated from kraft black liquor - The key to a new Biofuel. Filtration, 7, 309-315.
• Patil, R. A. (2016). Production of Acetic Acid in Kraft Pulp Mill Biorefinery Using. Orono, USA:
University of Maine.
• Paul, E., Atiemo-Obeng, V., & Kresta, S. M. (2003). Handbook of Industrial Mixing. Wiley-
Interscience.
• Ressources naturelles, Q. (2013). La filière du bioraffinage : une chimie verte. Ministère du Forêt,
Faune et Parcs au Québec (mffp).
• Saito, T., Brown, R. H., Hunt, M. A., Pickel, D. L., Pickel, J. M., Messman, J. M., . . . Naskar,
A. K. (2012). Turning renewable resources into value-added polymer: development of. Green
Chemistry, 3239-3454.
• Sameni, J., Krigstin, S., Rosa, D., Leao, A., & Sain, M. (2014). Thermal Characteristics of Lignin
Residue from Industrial processes. Bioressources 9(1).
• Sarkanen, K., & Ludwig, C. (1971). Lignins: Occurence, Formation, Structure and Reactions.
New York: Wiley.
• Smolarski, N. (2012). High-Value Opportunities for Lignin: Unlocking its potential. Frost &
Sullivan.
67
• Theliander, H. (2007). Recovery of cooking chemicals; the treatment and burning of black liquor.
In Ljungberg textbook. Pulp and Paper Chemistry and Technology. Book 2. Pulping Chemistry
and Technology Ed(s). Ek M, Gellerstedt G and Henriksson G. Stockholm: Fiber and Polymer
Technology, KTH.
• Tomani, P. (2010). The LignoBoost process. Cellulose Chemistry and Technology, 44, 53-58.
• Tuck, D., & Hunter, C. (2015). Extracting Lignin for New Forest Products. NRCAN.
• Uhl, V., & Gray, J. (1996). Mixing theory and practice. USA: Academic Press.
• Valmet. (2015). Valmet-supplied LignoBoost plant now handed over to Stora Enso's Sunila mill
in Finland. Suède: Valmet.
• Xu, T. (2005). Ion exchange membranes: State of their development and perspective. Journal of
Membrane Science, 1-29.
• Zhu, W. (2013). Equilibrium of Lignin Precipitation - The effects of pH, Temperature, Ion
Strength and Wood origin. Suède: Chalmers Univeristy.
• Zhu, W. (2015). Precipitation of Kraft Lignin - Yield and Equilibrium. Gotherburg, Sweden:
Chalmers University of Technology.
68
Afin de calculer l’erreur sur ∆P, le système a été roulé avec une conduite de longueur connue. La
perte de charge a été calculée et mesurée. Puis l’erreur a été calculé en se basant sur l’écart entre
les deux valeurs tout en prenant en compte l’erreur du capteur (2.5%)
L’erreur sur ∆P
∆𝑃𝑐𝑎𝑙𝑐𝑢𝑙é𝑒 − ∆𝑃𝑚𝑒𝑠𝑢𝑟é𝑒
𝜕(∆𝑃) = ∆𝑃 ∗ [( ) ∗ 100 + 2.5%]
∆𝑃𝑐𝑎𝑙𝑐𝑢𝑙é
L’erreur sur Ne
∆𝑃 𝐷ℎ
𝑁𝑒 = ∗
𝜌. 𝑣 2 2𝐿
𝜕𝑁𝑒 𝐷ℎ 𝜕∆𝑃 2 𝜕𝑣 2
=( √
)∗ ( ) + 2( )
𝑁𝑒 𝜌. 2𝐿 ∆𝑃 𝑣
Calcul de R2
𝑆𝑆𝐸
𝑅2 = 1 −
𝑆𝑆𝑇