0% ont trouvé ce document utile (0 vote)
41 vues90 pages

Study

Transféré par

Yassine Alami Tahiri
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd
0% ont trouvé ce document utile (0 vote)
41 vues90 pages

Study

Transféré par

Yassine Alami Tahiri
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd

Électrodialyse de la liqueur noire avec membrane bipolaire-

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]

Document en libre accès dans PolyPublie


Open Access document in PolyPublie

URL de PolyPublie:
[Link]
PolyPublie URL:
Directeurs de
recherche: Louis Fradette
Advisors:
Programme:
Génie chimique
Program:

Ce fichier a été téléchargé à partir de PolyPublie, le dépôt institutionnel de Polytechnique Montréal


This file has been downloaded from PolyPublie, the institutional repository of Polytechnique Montréal
[Link]
UNIVERSITÉ DE MONTRÉAL

ÉLECTRODIALYSE DE LA LIQUEUR NOIRE AVEC MEMBRANE BIPOLAIRE –


CARACTÉRISATION DE L’HYDRODYNAMIQUE D’UNE CELLULE ET

RELATION À LA PERFORMANCE

LAURA BOUCHRA ORABI

DÉPARTEMENT DE GÉNIE CHIMIQUE

ÉCOLE POLYTECHNIQUE DE MONTRÉAL

MÉMOIRE PRÉSENTÉ EN VUE DE L’OBTENTION

DU DIPLÔME DE MAÎTRISE ÈS SCIENCES APPLIQUÉES

(GÉNIE CHIMIQUE)

JANVIER 2018

© Laura Bouchra Orabi, 2018.


ii

UNIVERSITÉ DE MONTRÉAL

ÉCOLE POLYTECHNIQUE DE MONTRÉAL

Ce mémoire intitulé :

ÉLECTRODIALYSE DE LA LIQUEUR NOIRE AVEC MEMBRANE BIPOLAIRE –


CARACTÉRISATION DE L’HYDRODYNAMIQUE D’UNE CELLULE ET RELATION À LA

PERFORMANCE

Présenté par : ORABI Laura Bouchra

en vue de l’obtention du diplôme de : Maîtrise ès Sciences Appliquées

a été dûment accepté par le jury d’examen constitué de :

M. PERRIER Michel, Ph. D., président

M. FRADETTE Louis, Ph. D., membre et directeur de recherche

M. LABRECQUE Raynald, Ph. D., membre


iii

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.

L’objectif général de ce projet de recherche est d’investiguer si l’hydrodynamique dans la cellule


pourrait avoir une incidence sur l’encrassement de la membrane et conséquemment sur la durée
d’opération de la cellule. Plusieurs étapes étaient nécessaires afin de pouvoir atteindre cet objectif.
vi

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.

La principale recommandation pour des travaux futurs consiste à pousser l’analyse


hydrodynamique de l’électrodialyse de la liqueur noire en régime turbulent dans le but de trouver
des paramètres d’opération optimaux qui permettent de retarder voire d’empêcher l’encrassement
de la membrane dans la cellule.
vii

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.

An innovative technology of bipolar membrane electrodialysis is beginning to be studied more and


more around the world. This method has significant advantages over the different extraction
methods available. However, several limitations prevent the prolonged operation of electrodialysis
and reduce the efficiency of the cell. The major limitation is the membrane fouling. It is a
phenomenon of deposition of precipitated lignin particles on the surface of the membrane thus
causing an increase in the electrical resistance in the cell and the eventual stoppage of the operation.
It is therefore essential to be able to find a solution to remedy the membrane fouling in order to go
forward with the development of this technology.

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

TABLE DES MATIÈRES

DÉDICACE ................................................................................................................................... III

REMERCIEMENTS ..................................................................................................................... IV

RÉSUMÉ........................................................................................................................................ V

ABSTRACT .................................................................................................................................VII

TABLE DES MATIÈRES ............................................................................................................ IX

LISTE DES TABLEAUX ............................................................................................................XII

LISTE DES FIGURES ............................................................................................................... XIII

LISTE DES SIGLES ET ABRÉVIATIONS ............................................................................. XIV

LISTE DES ANNEXES ............................................................................................................XVII

CHAPITRE 1 INTRODUCTION ............................................................................................... 1

1.1 Contexte actuel de l’industrie des pâtes et papiers ........................................................... 1

1.1.1 Le procédé kraft ........................................................................................................... 2

1.1.2 Extraction de la lignine à partir de la liqueur noire ...................................................... 4

1.2 Problématique................................................................................................................... 5

1.3 Objectif principal du projet de recherche ......................................................................... 5

1.4 Plan du mémoire............................................................................................................... 6

CHAPITRE 2 REVUE DE LA LITTÉRATURE ....................................................................... 7

2.1 Précipitation de la lignine ................................................................................................. 7

2.1.1 La lignine...................................................................................................................... 7

2.1.2 La liqueur noire ............................................................................................................ 9

2.1.3 Les technologies d’extraction de la lignine à partir de la liqueur noire ..................... 10

2.1.4 Procédés commerciaux ............................................................................................... 13


x

2.2 L’électrodialyse de la liqueur noire ................................................................................ 14

2.2.1 Notions sur l’électrodialyse ........................................................................................ 14

2.2.2 Paramètres d’opération affectant l’opération de l’électrodialyse ............................... 19

2.2.3 Difficultés d’opération ............................................................................................... 20

2.3 L’hydrodynamique de la cellule ..................................................................................... 24

2.3.1 Perte de charge dans la cellule ................................................................................... 25

2.3.2 Régime d’écoulement dans la cellule ......................................................................... 26

2.3.3 Comportement rhéologique ........................................................................................ 26

2.3.4 Courbe de puissance et coefficient de Metzner et Otto .............................................. 27

2.4 Problématique et objectifs spécifiques ........................................................................... 28

CHAPITRE 3 RESSOURCES MATÉRIELLES ET MÉTHODOLOGIE .............................. 29

3.1 Montage expérimental et composantes .......................................................................... 29

3.1.1 Conception de la cellule électrochimique .................................................................. 31

3.1.2 Tuyauterie et système de pompage ............................................................................ 34

3.1.3 Système d’acquisition des données ............................................................................ 35

3.1.4 Agencement des étages .............................................................................................. 35

3.2 Caractérisation hydrodynamique de la cellule ............................................................... 36

3.2.1 Caractéristiques des solutions références ................................................................... 36

3.2.2 Perte de charge ........................................................................................................... 36

3.2.3 Fluide non newtonien ................................................................................................. 37

3.2.4 Détermination de la constante de cisaillement Ks ..................................................... 38

3.3 Procédure expérimentale de l’électrodialyse de la liqueur noire ................................... 39

3.3.1 Caractérisation de la liqueur noire ............................................................................. 39

3.3.2 Préparation des solutions pour l’électrodialyse .......................................................... 42


xi

3.3.3 Le protocole suivi pour l’électrodialyse ..................................................................... 43

CHAPITRE 4 RÉSULTATS .................................................................................................... 45

4.1 Caractérisation hydrodynamique de la cellule ............................................................... 45

4.1.1 Géométrie I : sans grille ............................................................................................. 45

4.1.2 Détermination de Ks................................................................................................... 47

4.1.3 Géométrie II : avec grille #1 ...................................................................................... 48

4.1.4 Géométrie avec grille #2 ............................................................................................ 50

4.1.5 Analyse générale pour les trois géométries ................................................................ 53

4.2 Résultats de la caractérisation de la liqueur noire .......................................................... 56

4.2.1 Total des solides dissous ............................................................................................ 56

4.2.2 Masse volumique et teneur en lignine ........................................................................ 56

4.2.3 Comportement rhéologique de la liqueur noire.......................................................... 56

4.3 Résultats de l’électrodialyse de la liqueur noire............................................................. 57

CHAPITRE 5 CONCLUSION ET RECOMMANDATIONS ................................................. 61

BIBLIOGRAPHIE ........................................................................................................................ 63
xii

LISTE DES TABLEAUX

Tableau 1-1 : Composition typique des solutés de la liqueur blanche ............................................. 2

Tableau 2-1 - Différentes applications potentielles de la lignine ..................................................... 8

Tableau 2-2 - Pourcentage massique des matières organiques dans la liqueur noire ...................... 9

Tableau 2-3 - Avantages et désavantages des méthodes d'extraction de la lignine ....................... 12

Tableau 2-4 - Les différents mécanismes d’encrassement, leur fondement phénoménologique et


leur effet sur le transfert de masse .......................................................................................... 23

Tableau 3-1 - Constituants de la cellule électrochimique .............................................................. 31

Tableau 3-2 - Caractéristiques des solutions utilisées pour la caractérisation ............................... 36

Tableau 3-3 - Propriétés du CMC 1% ............................................................................................ 38

Tableau 3-4 - pH des échantillons préparés ................................................................................... 42

Tableau 4-1 - Relation Ne-Re pour géométrie I ............................................................................. 47

Tableau 4-2 - Relation Ne-Re pour géométrie II ........................................................................... 49

Tableau 4-3 - Relation Ne-Re pour géométrie III .......................................................................... 52

Tableau 4-4 - Sommaire des relations Ne-Re pour les trois géométries ........................................ 54

Tableau 4-5 - Ks pour les trois géométries .................................................................................... 55

Tableau 4-6 - Valeurs de Ks pour certains mélangeurs industriels ................................................ 55

Tableau 4-7 - Masse volumique et teneur en lignine de la liqueur noire utilisée .......................... 56

Tableau 4-8 - Paramètres d'opération ............................................................................................. 57

Tableau 4-9 - Débits d'opération pour l'électrodialyse de la liqueur noire .................................... 58


xiii

LISTE DES FIGURES

Figure 1-1 : Schéma bloc du procédé kraft ...................................................................................... 4

Figure 2-1 - Le processus d'électrodialyse ..................................................................................... 15

Figure 2-2 - Membrane bipolaire reproduit de (Patil, 2016) .......................................................... 17

Figure 2-3 - Schématisation du processus de l'électrodialyse de la liqueur noire à membrane


bipolaire .......................................................................................................................................... 18

Figure 2-4 - Profils de vitesse et de concentration dans la cellule d'ED ........................................ 19

Figure 2-5 - Mécanismes d'encrassement ...................................................................................... 22

Figure 3-1- Schéma du système complet ...................................................................................... 30

Figure 3-2 - Séparateur vu en haut ................................................................................................. 32

Figure 3-3 - Séparateur vu à côté ................................................................................................... 32

Figure 3-4 - Agencement des étages dans la cellule ...................................................................... 35

Figure 4-1 - Perte de charge en fonction du débit pour la géométrie I .......................................... 45

Figure 4-2 - Ne en fonction de Re pour la géométrie I .................................................................. 46

Figure 4-3 - Ks en fonction de Re pour la géométrie I .................................................................. 47

Figure 4-4 - Perte de charge en fonction du débit pour la géométrie II ......................................... 48

Figure 4-5 - Ne en fonction de Re pour la géométrie II ................................................................. 49

Figure 4-6 - Ks en fonction de Re pour la géométrie II ................................................................. 50

Figure 4-7 - Perte de charge en fonction du débit pour la géométrie III ........................................ 51

Figure 4-8 - Ne en fonction de Re pour la géométrie III................................................................ 52

Figure 4-9 - Ks en fonction de Re pour la géométrie III ................................................................ 53

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

LISTE DES SIGLES ET ABRÉVIATIONS

SYMBOLE DESCRIPTION UNITÉ

A Aire transversale de l’écoulement m2

Dh Diamètre hydraulique m

f Facteur de friction ---

k Indice de consistance [Link]

KP Constant de géométrie ---

Ks Constante de géométrie (coefficient de Metzner et Otto) ---

L Longueur caractéristique m

m Masse Kg

n Indice de comportement ---

P Perte de charge Pa

Q Débit m3/s

v Vitesse d’écoulement m/s

V Volume m3
xv

NOMBRES ADIMENSIONNELS :

SYMBOLE DESCRIPTION

Ne Nombre de Newton

Re Nombre de Reynolds

SYMBOLES GRECS :

SYMBOLE DESCRIPTION UNITÉ

∆P Perte de charge Pa

γ̇ Taux de cisaillement 1/s

γ̇ eff Taux de cisaillement effectif 1/s

μ Viscosité Pa.s

μeff Viscosité effective (viscosité du procédé) Pa.s

ρ Masse volumique kg/m3

τ Contrainte de cisaillement Pa
xvi

ABRÉVIATIONS

ABRÉVIATION DESCRIPTION

AEM Membrane anionique

BPM Membrane bipolaire

CEM Membrane cationique

ED Électrodialyse

EDMBP Électrodialyse à membrane bipolaire

LN Liqueur noire

P&P Pâte et papier

SDT Solides dissous totaux


xvii

LISTE DES ANNEXES

ANNEXE A – FICHE TECHNIQUE DES POMPES UTILISÉES ............................................. 68

ANNEXE B – DÉTAILS DU CALCUL DE L’ERREUR ........................................................... 71


1

CHAPITRE 1 INTRODUCTION

1.1 Contexte actuel de l’industrie des pâtes et papiers


Avec l’essor de la technologie des communications électroniques et la décroissance de
l’utilisation du papier, l’industrie des pâtes et papier (P&P) connait actuellement la plus importante
crise de son histoire. Les produits forestiers contribuaient à plus que 3% du produit intérieur brut
du Canada il y a deux décennies, alors qu’en 2014, ce chiffre était inférieur à 1%. Durant la même
période, plus que 40% des employés dans ce secteur ont perdu leur travail au Canada. (Kennan,
Parkinson, & Jang, 2014)

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.

Le bioraffinage correspond à une série de procédés industriels qui permettent de transformer la


matière ligneuse en une variété de produits chimiques commerciaux. Dans une bioraffinerie, les
principaux constituants de base du bois (cellulose, hémicellulose, lignine et matières extractibles)
servent à la fabrication de produits à valeur ajoutée tels que des biocombustibles, des produits
chimiques spéciaux, des matières premières de base pour d’autres industries et des précurseurs de
matières plastiques ainsi que d’autres polymères. Les activités de bioraffinage complètent et
diversifient donc la production de pâtes et papiers. (Ressources naturelles, 2013)

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.

1.1.1 Le procédé kraft

Le procédé kraft est un processus de transformation du bois en pâte puis éventuellement en


papier. Aujourd'hui, le procédé kraft est le procédé de fabrication de la pâte chimique le plus
répandu dans le monde en raison des meilleures propriétés de résistance de la pâte obtenue par
rapport aux autres procédés, de son application à toutes les espèces de bois, ainsi qu'aux systèmes
efficaces de récupération des produits chimiques qui ont été développés et mis en œuvre. (INERIS,
2017)

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)

Tableau 1-1 : Composition typique des solutés de la liqueur blanche

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)

- Besoins énergétiques moins élevés que les digesteurs batch

- Plus faible volume des courants entrants et sortants

- Opération en régime permanent

- Efficacité énergétique plus élevée

- Plus de flexibilité dans la configuration du procédé.

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)

En fin de combustion, des cendres composées essentiellement de sulfure de sodium et de carbonate


de sodium sont dissoutes dans l’eau pour former la liqueur verte. Cette liqueur est régénérée par
caustification afin de récupérer la liqueur blanche nécessaire pour la cuisson. (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é

Liqueur Cuisson Liqueur Chaudière de


Évaporateurs

Liqueur noire faible


(Digesteur) noire récupération
blanche
forte

Liqueur verte vers


caustification

Unité de Unité de Séchage et Produit


lavage Pâte blanchissement Pâte couchage fini
lavée blanchie

Figure 1-1 : Schéma bloc du procédé kraft

1.1.2 Extraction de la lignine à partir de la liqueur noire

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)

L’extraction de la lignine à partir de la liqueur noire présente plusieurs avantages :

- Augmenter la capacité de l’usine par le désengorgement de la chaudière de récupération :

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.

L’électrodialyse à membrane bipolaire (EDMBP) est un procédé membranaire dont la force


motrice est une différence de potentiel électrique. Cette technologie permet de baisser le pH à l’aide
de la dissociation de l’eau. Toutefois, ce processus est encore à l’échelle laboratoire et plusieurs
défis sont encore à l’étude afin de l’optimiser. Parmi les défis figure l’encrassement de la membrane
ionique qui oblige l’arrêt de l’opération après un certain temps pour un nettoyage.

1.3 Objectif principal du projet de recherche


Sachant que l’extraction de la lignine de la liqueur noire permet une opération à plus haute
charge des usines de pâte kraft, les procédés de séparation de la lignine sont d’un grand intérêt.
Considérant que les procédés à base d’électricité sont d’un intérêt certain pour le Québec et, dans
une moindre mesure, pour le Canada, l’électrodialyse à membrane bipolaire présente un potentiel
très intéressant économiquement pour une application industrielle. Enfin, ce type de procédé est
connu pour être limité dans son opération par l’encrassement de la membrane bipolaire, un
problème qui limite l’applicabilité de la technologie. Certains essais ont permis de penser que
l’hydrodynamique de l’écoulement dans la cellule pourrait avoir une influence sur la performance
et la durée de l’opération car elle aurait un effet sur la vitesse de déposition des particules sur les
membranes et donc sur l’encrassement responsable de la performance inadéquate des cellules de
ce type. Dans ce contexte, l’objectif de ce projet est de déterminer si les paramètres
hydrodynamiques de l’opération d’une cellule d’électrodialyse à membrane bipolaire ont un impact
sur la durée de l’opération possible de l’électrodialyse de la liqueur noire.
6

1.4 Plan du mémoire


Le présent mémoire est organisé en 5 chapitres. De plus amples détails sur les concepts
présentés à l’introduction sont fournis dans la revue de la littérature (chapitre 2)

Le chapitre 3 du mémoire décrit la méthodologie adaptée afin d’atteindre l’objectif principal du


projet. Ce chapitre est divisé en trois parties, chaque partie étant reliée à un objectif spécifique. La
première partie traite la conception de la cellule d’électrodialyse et de la construction du montage
expérimental. La deuxième partie traite de la caractérisation hydrodynamique de la cellulealors que
la dernière partie se concentre sur la méthodologie de la réalisation des expériences.

Le Chapitre 4 du mémoire présente la synthèse des principaux résultats du projet.

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

CHAPITRE 2 REVUE DE LA LITTÉRATURE

2.1 Précipitation de la lignine

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)

L’acidification est considérée un processus efficace et économique afin d’extraire la lignine de la


liqueur noire. Durant ce processus, le pH de la liqueur noire est réduit jusqu’à 9 environ, ce qui
permet à la lignine de se précipiter sous forme de particules qui peuvent être séparées par la suite
par la filtration après une étape de coagulation. (Arkall, Olsson, & Wallberg, 2014)

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)

Le mécanisme de la précipitation de la lignine est principalement gouverné par la protonation des


groupements phénoliques ionisés sur la molécule de la lignine. (Zhu W. , 2013)

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,

1- Gaz de synthèse (par exemple, par gazéification)

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)

Une représentation non-exhaustive de la grande variété d'applications potentielles de la lignine est


illustrée dans le Tableau 2-1 ci-dessous (Gosselink, 2011)

Tableau 2-1 - Différentes applications potentielles de la lignine

Lignine

Produits de gaz de
Hydrocarbures Phénols Produits oxydés Macromolécules
synthèse

• Méthanol • Benzène • Substituants du • Vanilline • Fibres de carbone


phénol
• Éthanol • Toluène • Acides • Substituants de
• Crésols aromatiques polymère
• Liquide Fisher- • Xylène
Tropsch • Eugénols • Aldéhydes • Charbon activé
• Cyclohexane
• Syringols • Quinones • Adhésives
• Styrène
• Guaiacols • Cyclohexanol • Résines
9

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)

- Les propriétés thermiques, telles que le pouvoir calorifique, la température de transition


vitreuse, la polydispersité, la masse molaire, etc.

- Le degré de pureté et la teneur en cendres

- La structure et les groupes fonctionnels

2.1.2 La liqueur noire

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

Matériel organique % massique

Lignine 30-40

Acides hydroxyliques 25-35

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.

2.1.3 Les technologies d’extraction de la lignine à partir de la liqueur noire

Il existe plusieurs méthodes qui permettent l’extraction de la lignine de la liqueur noire,


soient, la précipitation acide, la précipitation par CO2, l’ultrafiltration, la méthode électrochimique
et la méthode biologique.

Les technologies biologiques impliquent l’utilisation de champignons et des bactéries. Cependant,


l’efficacité des microorganismes est limitée par la présence de composés toxiques dans la liqueur
noire, par sa très forte alcalinité et force ionique. (Da Re & Papinutti, 2011) Le traitement
biologique reste ainsi un cas d’étude et ne sera pas présenté en détails.

[Link] La précipitation acide avec l’acide sulfurique

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

[Link] La précipitation acide avec le CO2

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'ultrafiltration est un type de filtration membranaire dans lequel la pression hydrostatique


force un liquide contre une membrane semi-perméable. Une membrane semi-perméable est une
fine couche de matériau capable de séparer des substances lorsqu'une force motrice est appliquée
à travers cette dernière. (APEC, 2017) L’ultrafiltration de la liqueur noire a fait l’objet de plusieurs
recherches en Suède (Holmqvist, Wallberg, & Jonsson, 2005).

[Link] La méthode électrochimique

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)

L’électrodialyse conventionnelle de la liqueur noire a été étudiée par Cloutier et al.


(Cloutier, Champagne, & Labrecque, 2013) 4 paires de membranes anionique/cationique ont été
placées alternativement entre les deux électrodes. L’étude a conclu que la liqueur noire ne peut pas
être complétement acidifiée par l’électrodialyse et que cette étape doit être suivie par une étape
d’acidification afin de réduire le pH au point de précipitation de la lignine. L’électrodialyse à
12

membrane bipolaire se diffère de l’électrodialyse conventionnelle par le type des membranes


utilisées dans la cellule. Cette technologie sera étudiée en profondeur dans le cadre de ce projet.

Le Tableau 2-3 ci-dessous sommaire les avantages et les désavantages de chaque méthode décrite
ci-haut.

Tableau 2-3 - Avantages et désavantages des méthodes d'extraction de la lignine

Méthode Avantages Désavantages

Précipitation acide • Rendement très élevé • Ajout d’un produit chimique


(H2SO4)
• Génération d’un résidu acide à traiter

• Filtrabilité difficile

Précipitation au • Coûts réduits • Limitation du pH qu’on peut


CO2 atteindre donc le rendement est limité
• Meilleure filtrabilité de la
lignine

Ultrafiltration • Ne nécessite pas l’ajout de • Très faible rendement


produits chimique
• Technologie coûteuse en capital et en
opération

Méthode • Ne nécessite pas l’ajout de • Stabilité des membranes ioniques


électrochimique produit chimique
• Encrassement de la membrane
• Consommation
énergétique modérée

• Récupération du NaOH qui


pourrait être réinjecté dans
le procédé Kraft
13

2.1.4 Procédés commerciaux

[Link] LignoBoost

Le procédé LignoBoost est un procédé d’acidification de la liqueur noire développé par


Innventia en collaboration avec l’université de Chalmers en suède.

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é.

Les avantages majeurs du procédé LignoBoost sont : (Tomani, 2010)

1) Meilleur rendement en lignine produite;

2) La lignine obtenue a un teneur inférieur en cendres et en carbohydrates;

3) Investissement moins couteux à cause de la réduction de la surface du filtre et le volume du


lavage acide.

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

Le procédé LignoForce a été développé par FPInnovations, ce procédé propose une


amélioration de la filtrabilité de la lignine précipitée par l’oxydation de la liqueur noire.

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)

1) Réduire les besoins en acide pendant les étapes de la précipitation et du lavage

2) Minimisation des émissions sulfuriques pendant le traitement de la liqueur noire

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)

2.2 L’électrodialyse de la liqueur noire


L’électrodialyse d’une façon générale est un procédé membranaire basé sur le transport des
ions à travers une membrane semi-perméable avec comme force motrice un potentiel électrique.
Les membranes utilisées sont appelées des membranes échangeuses d’ions ou membranes
ioniques. (Brun, 1989) (Lamarche & Boulet, 1992)

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)

2.2.1 Notions sur l’électrodialyse

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.

En général, l’électrodialyse est effectuée dans une cellule constituée de 3 compartiments


principaux ; un compartiment pour l’électrolyte, dont la fonction est d’éliminer les produits se
formant au niveau des électrodes tout en maintenant un milieu conducteur, un deuxième pour le
concentrat et un troisième pour le diluat tel que présenté dans la figure ci-dessous.
15

Concentrat (sortie)

CEM : Membrane cationique


AEM : Membrane anionique
Concentrat (entrée)

Figure 2-1 - Le processus d'électrodialyse

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)

[Link] Les membranes échangeuses d’ions

Les membranes ioniques sont constituées d’un matériau macromoléculaire (squelette)


portant des groupes ionisables. La membrane contient des ions fixes, solidaires du squelette,
électriquement neutralisés par des charges mobiles de signe opposé appelés les contre-ions.
(Bazinet et al. , 1999) Ces derniers sont positifs dans le cas d’une membrane cationique et négatifs
16

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)

Les membranes doivent être performantes et résistantes. Leurs caractéristiques sont


fonction de quatre propriétés principales : (Quoc, 2004)

• Une imperméabilité par rapport aux co-ions ;


• Une conductivité élevée ;
• Une bonne résistance mécanique et une stabilité dimensionnelle ;
• Bonne stabilité chimique et physique.

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.

Il existe trois types différents de membranes échangeuses d'ions. (Xu, 2005)

1) Les membranes échangeuses de cation (membranes cationiques)


Les membranes échangeuses de cations contiennent des groupes chargés négativement fixés sur
une matrice polymère hydrophobe. Elles sont perméables aux cations uniquement.

2) Les membranes échangeuses d’anions (membranes anioniques)


Les membranes échangeuses d'anions contiennent des groupes chargés positivement fixés à une
matrice polymère hydrophobe. Ce type de membrane est uniquement perméable aux anions.

3) Les membranes bipolaires

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

La combinaison d’une membrane bipolaire et d’une membrane cationique ou anionique constitue


un compartiment pour l’acidification ou l’alcalinisation des solutions traitées.

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
+
+ - - -
+
+
+ - -
+
+ -
-
+
+ -
-
+
+ - -

Figure 2-2 - Membrane bipolaire reproduit de (Patil, 2016)

[Link] Les séparateurs

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

[Link] Électrodialyse à membrane bipolaire

L’électrodialyse à membrane bipolaire EDMBP est une technologie émergente pour le


traitement des solutions afin de générer des acides et des bases. Cette technologie a récemment
commencé à être examinée pour l’extraction de la lignine à partir de la liqueur noire.

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.

Liqueur noire NaOH


acidifié concentré

CEM BPM CEM

H+ OH- Anode (+)


Cathode (-)
=

Liqueur noire NaOH


Figure 2-3 - Schématisation du processus de l'électrodialyse de la liqueur noire à membrane
bipolaire
19

2.2.2 Paramètres d’opération affectant l’opération de l’électrodialyse

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)

[Link] L’intensité du courant et le voltage

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.

[Link] La vitesse d’écoulement

La vitesse d’écoulement a un impact significatif sur la performance de l’électrodialyse.

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)

Figure 2-4 - Profils de vitesse et de concentration dans la cellule d'ED


20

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.

[Link] La concentration des solutions alimentées

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.

2.2.3 Difficultés d’opération

Plusieurs difficultés d’opération sont rencontrées lors de l’électrodialyse de la liqueur noire.

[Link] Imperfection de la permsélectivité des membranes

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)

[Link] Effet Joule

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

la solution et dans la membrane. La polarisation cause une augmentation importante de la densité


du courant et autres conséquences aussi nuisibles tel que : (Granier & Bibenet, 1976)

• L’augmentation de la consommation d’énergie due à l’augmentation de la résistance


électrique des couches limites ;
• Un risque de précipitation provoquant le colmatage des membranes ;
• Une différence importante du pH de part et d’autre des membranes, ce qui cause une
détérioration de la membrane ;
• Une diminution majeure du rendement de transport et de la permsélectivité

[Link] L’encrassement

Ce phénomène pourrait être perçu comme étant une polarisation secondaire. Les causes de
dépôt sont multiples :

• Entartrage par formation d’hydroxyde métallique due à l’élévation du pH et la


précipitation de sels ;
• Encrassement par des macro-ions de masse moléculaire élevée ;
• Empoisonnement par des ions organiques hydrophobes et peu mobiles.

L’encrassement des membranes cause un accroissement de la résistance électrique dans la cellule


et une diminution du rendement ce qui altère réversiblement ou irréversiblement les membranes.
(Granier & Bibenet, 1976)

[Link].1 Les différentes formes de l’encrassement

L’encrassement de la membrane pourrait prendre plusieurs formes : (Field, 2010)

1) L’adsorption : apparait quand une interaction spécifique entre la membrane et la solution


ou les particules a lieu.
2) Blocage des pores : que ça soit partiel ou total, il entraine une réduction du flux.
3) La déposition : un dépôt de particule peut avoir lieu causant ainsi une résistance hydraulique
additionnelle. Ce phénomène est appelé «résistance de gâteau»
4) La formation du gel : pour certaines macromolécules, le niveau de la polarisation de la
concentration peut causer une formation d’un gel à la surface de la membrane. Ce
phénomène est observé dans le cas de molécules biologiques.
22

[Link].2 Influents clés de l’encrassement

L’encrassement entraine une augmentation de la résistance électrique et hydraulique ce qui


implique un courant plus faible pour un potentiel électrique constant ou bien une augmentation du
potentiel si le courant est constant.

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)

Plusieurs paramètres influencent l’encrassement de la membrane tel que :

- La nature et la concentration du soluté et du solvant


- Le type de membrane
- La distribution de la taille des pores
- Les caractéristiques et le matériau de la membrane
- Les hydrodynamiques à la surface de la membrane

[Link].3 Les mécanismes de l’encrassement

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)

Figure 2-5 - Mécanismes d'encrassement

Le Tableau 2-4 ci-dessous résume ces mécanismes et leur effet sur le transfert de masse. (Field,
2010)
23

Tableau 2-4 - Les différents mécanismes d’encrassement, leur fondement phénoménologique et


leur effet sur le transfert de masse

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

Formation de gâteau Formation d’un gâteau de • La résistance devient une


(figure 2-4.d) particules à la surface de la résistance du gâteau plus celle
membrane de la membrane

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

• Composante complètement irréversible : correspond à un colmatage en profondeur qui ne


peut être éliminé par simple lavage.

Cette distinction permet de définir la notion du débit critique.

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)

2.3 L’hydrodynamique de la cellule


L’hydrodynamique est une branche de la physique qui a pour objet l’étude des fluides en
mouvement. (Larousse, 2017) Les fluides en mouvement possèdent une énergie cinétique qui
1
augmente avec l’augmentation de la vitesse du fluide selon 𝐸𝑘 = 2 𝑚. 𝑣 2 . C’est à partir de cette

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

2.3.1 Perte de charge dans la cellule

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.

Il faut d’abord commencer par définir l’expression du nombre de Reynolds et du nombre de


puissance (ou nombre de Newton) pour un fluide newtonien. Le nombre de Reynolds est défini
comme étant le rapport entre les forces d’inertie et les forces visqueuses.

𝜌𝑣𝐷ℎ Équation 2-1


𝑅𝑒 =
𝜇

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)

Le nombre de Newton pourrait être exprimé par l’équation suivante :

𝐷 1 Équation 2-2
𝑁𝑒 = 2𝑓 = ∆𝑃 ( ) ∗ ( )
𝐿 𝜌𝑣 2

f étant le facteur de friction de Fanning.

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

permet conséquemment le calcul de la viscosité. Ce qui permettra de calculer un nombre de


Reynolds et une perte de charge.

2.3.2 Régime d’écoulement dans la cellule

Le régime d’écoulement dépend de plusieurs paramètres tels que la vitesse du fluide, sa


viscosité et sa masse volumique ainsi que le diamètre de la conduite.

Il est nécessaire de déterminer le nombre de Reynolds critique (laminaire/turbulent) pour identifier


le régime d’écoulement dans une géométrie car ce régime aura un impact sur le comportement du
fluide. Le régime d’écoulement est défini par la suite par rapport au nombre de Re.

Si Re < ReC → le régime d’écoulement est laminaire

Si Re > ReC → le régime d’écoulement est turbulent

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.

2.3.3 Comportement rhéologique

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

𝜏 étant la contrainte de cisaillement, 𝜇 la viscosité du fluide et 𝛾̇ le taux de cisaillement

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.

Les fluides non newtoniens peuvent être distincts selon 2 groupes :

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

2.3.4 Courbe de puissance et coefficient de Metzner et Otto

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:

𝛾̇ 𝑒𝑓𝑓 = 𝐾𝑠 ∗ 𝑁 Équation 2-4

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)

La valeur de Ks est déterminée selon les étapes suivantes : (Fradette, 1999)

1. Établir des courbes de puissances Ne vs. Re pour un fluide newtonien

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)

3. Calculer un nombre de Reynolds effective à l’aide de la relation obtenue par la courbe de


puissance en (1) selon la relation

4. Calculer une viscosité effective à partir de Reeff selon l’équation :

𝐷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.

5. Utiliser la courbe d’écoulement du fluide non newtonien pour obtenir la vitesse de


cisaillement effective correspondant à cette viscosité.

6. Déterminer Ks à l’aide de l’Équation 2-4.

2.4 Problématique et objectifs spécifiques


L’encrassement des membranes se présente comme un défi majeur pour le développement de la
technologie d’électrodialyse car l’encrassement oblige l’arrêt fréquent de l’opération pour nettoyer
les membranes. L’hydrodynamique dans la cellule semble avoir un impact sur la couche limite qui
se crée à la surface des membranes et donc sur le dynamisme du dépôt des particules de la lignine
précipitée. La promotion de la turbulence, à notre avis, va améliorer la performance de l’opération
et conséquemment prolonger la durée de l’opération.

À 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 :

1- Concevoir et Construire une cellule d’électrodialyse à membrane bipolaire et du montage


associé

2- Caractériser le comportement hydrodynamique de la cellule conçue

3- Déterminer l’impact de l’hydrodynamique dans la cellule sur la vitesse d’encrassement de


la membrane

La réalisation de ces trois objectifs nous permettra d’atteindre l’objectif principal du projet.
29

CHAPITRE 3 RESSOURCES MATÉRIELLES ET MÉTHODOLOGIE

Afin d’atteindre les objectifs du présent projet de recherche, le travail a été divisé en trois phases :

- Phase I : Conception et construction du montage expérimental

- Phase II : Caractérisation hydrodynamique de la cellule construite

- Phase III : Traitement de la liqueur noire

3.1 Montage expérimental et composantes


Le montage expérimental a été conçu et construit dans le cadre de cette maîtrise. Ce montage est
constitué de trois parties :

• La cellule d’électrodialyse : qui représente le cœur du processus électrochimique. Elle


représente le réacteur dans lequel se tiennent toutes les réactions

• Le système de pompage : incluant les pompes permettant la circulation du liquide dans le


système, la tuyauterie ainsi que toute l’instrumentation nécessaire pour monitorer
l’opération

• Le système d’acquisition des données : qui correspond à l’interface entre les capteurs et
l’ordinateur

Une schématisation du système complet est illustrée dans la figure ci-dessous.


30

Figure 3-1- Schéma du système complet


31

3.1.1 Conception de la cellule électrochimique

Il n’existe pas un protocole prédéfini pour réaliser la conception des cellules


d’électrodialyse. Il existe toutefois une référence qui propose des lignes directrices. Le design de
la cellule a été entièrement réalisé en prenant en compte ces dernières. (Kemperman, 2000)

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.

Une cellule d’électrodialyse à membrane bipolaire pourrait comprendre 4 compartiments et plus :

• 1 compartiment pour le diluat -compartiment acide- (ou plus)

• 1 compartiment pour le concentré -compartiment basique- (ou plus)

• 2 compartiments pour le rinçage des électrodes.

Les éléments suivants ont été pris en considération lors de la conception :

• Les débits des solutions à traiter


• Les conformités chimique et mécanique des matériaux utilisés
• La perte de charge dans le système
• La fuite entre les étages

Les différents constituants de la cellule sont présentés dans le Tableau 3-1ci-dessous.

Tableau 3-1 - Constituants de la cellule électrochimique

Dimensions en mm
Élément (longueur X largeur X Matériau
épaisseur)

Plaque extérieure anode/cathode 152.4 x 101.6 x 25.4 Acrylique

Anode/cathode 63.5 x 63.5 x 2 Titane platinisé

Séparateur 152.4 x 101.6 x 2.38 Acrylique


32

Grille #1 63.5 x 63.5 x 0.9 Polypropylène

Grille #2 63.5 x 63.5 x 2.38 Polymère HDPE

Polymère inconnu (fourni


Membrane ionique 152.4 x 101.6 x 0.2
par la compagnie Fumatech)

Joint d’étanchéité 152.4 x 101.6 x 0.4 Néoprène

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.

Figure 3-2 - Séparateur vu en haut Figure 3-3 - Séparateur vu à côté


33

Afin de faire varier le degré de turbulence à l’intérieur de la cellule, trois géométries sont
proposées :

- Géométrie I : avec des séparateurs sans grillage

- Géométrie II : avec des séparateurs et grille #1

- Géométrie III : avec des séparateurs grille #2.

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.

Figure 3-4 - Grille #1 Figure 3-5 - Grille #2

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

- Le matériel composant de la grille #1 est le plastique extensible mou (polypropylène)


tandis que la grille #2 est composé de plastique HDPE rigide

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

3.1.2 Tuyauterie et système de pompage

Tel qu’illustré à la figure 3.1, le système de pompage et d’instrumentation du montage expérimental


est composé des éléments suivants :

Pièce Fabricant Modèle/type Caractéristiques

Cuves --- Plastique HDPE • Capacité de 5 litres

Tuyaux (entre les --- PVC • ½" de diamètre

cuves et les pompes)

Amorce de la REAMTOP Caoutchouc noir


pompe

Pompe centrifuge BLDC • 0 – 8 l/min


DC50F-24150A
• Tête maximale = 15 m

Tubes (entre les --- PVC • 3/8" de diamètre

pompes et la cellule)

Débitmètre SEEED RB-See-98 • Plage de débit : 0-40


L/min
• Précision de 3%

Capteurs de NXP MPX5500 • Capteur de pression type

pression différentiel
• Plage de mesure 0-500
kPa
• Précision de 2.5%

Électrodes de pH DFRobot • Plage de mesure : 0 -14


RB-Dfr-483
• Précision : ±0.1 pH (à
25°C)
35

3.1.3 Système d’acquisition des données

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.

Le logiciel Labview de National Instruments permet de commander :

- Les pompes afin d’obtenir le débit désiré

- Le courant alimenté au système ou le voltage

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.

3.1.4 Agencement des étages

Un étage correspond à un élément dans la cellule électrochimique. Cet élément pourrait être un
séparateur, une membrane ou un joint.

La figure ci-dessous, illustre l’agencement des étages dans la cellule conçue.

Anode (+) Séparateur Membrane Membrane Joint Séparateur Cathode (-)


cationique bipolaire d’étanchéité

Figure 3-4 - Agencement des étages dans la cellule


36

3.2 Caractérisation hydrodynamique de la cellule


À nos connaissances, il n’y a jamais eu une caractérisation hydrodynamique d’une cellule
électrochimique. Toutefois, dans un contexte de traitement de solution qui peuvent contenir des
impuretés, l’encrassement de la membrane s’avère un problème très commun. L’hydrodynamique
semble être une piste avec un potentiel intéressant à examiner afin d’optimiser la performance de
la cellule dans le but de retarder voire d’empêcher l’encrassement de la membrane.

3.2.1 Caractéristiques des solutions références

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.

Tableau 3-2 - Caractéristiques des solutions utilisées pour la caractérisation

Eau Glucose 30% Glucose 42%

Masse volumique (kg/m3) 1000 1116 1166

Viscosité (Pa.s) 0.00100 0.00326 0.00545

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%.

3.2.2 Perte de charge

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 :

• Géométrie I : sans grille dans le séparateur


37

• Géométrie II : avec la grille #1

• Géométrie III : avec la grille #2

b) Calculer la vitesse pour chaque débit selon :

𝑄
𝑣 = Équation 3-1
𝐴𝑖𝑟𝑒 𝑡𝑟𝑎𝑛𝑠𝑣𝑒𝑟𝑎𝑙𝑒

c) Calculer Reynolds pour chaque débit selon l’Équation 2-1

À noter que le rayon hydraulique considéré pour calculer le nombre de Reynolds


correspond au ratio entre la surface et le périmètre de l’aire transversale à l’écoulement à
l’intérieur de la cellule.

d) Calculer le nombre de Newton selon l’Équation 2-2

Déterminer la relation entre Re et Ne à l’aide de la méthode des moindres carrées pondérées en


régime laminaire sous la forme : 𝑁𝑒 = 𝐾𝑃 ∗ 𝑅𝑒 −1 + 𝐾𝑁𝐸 Équation 3-2

Kp étant une constante qui dépend de la géométrie du système

Afin de déterminer la qualité de chaque régression, le coefficient de détermination R2 a été calculé


comme indicateur selon l’équation,

𝑆𝑆𝐸
𝑅2 = 1 − Équation 3-3
𝑆𝑆𝑇

Avec

𝑆𝑆𝑇 = ∑(𝑦𝑖 − 𝑦̅)2 Équation 3-4

𝑆𝑆𝐸 = ∑(𝑦𝑖 − 𝑦̂𝑖 )2 Équation 3-5

𝑦𝑖 𝑒𝑡 𝑦̂𝑖 étant les valeurs 𝑦originales


𝑖 𝑒𝑡 𝑦
̂𝑖 etétant les valeurs
modélisées
originales et modélisées
3.2.3 Fluide non-newtonien

Le fluide non newtonien permet de déterminer les valeurs de la constante de géométrie Ks en se


basant sur les courbes de puissance développées pour les fluides newtoniens. Le fluide non
38

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)

𝜇 = 𝑘 ∗ 𝛾̇ 𝑛−1 Équation 3-6

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.

Tableau 3-3 - Propriétés du CMC 1%

Masse volumique k n

1200 kg/m3 0.226 0.666

3.2.4 Détermination de la constante de cisaillement Ks

La procédure de détermination du Ks pour une cellule électrochimique suit le même


cheminement développé par Metzner et Otton pour les agitateurs mécaniques. La méthodologie
utilisée est la suivante :

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%)

3. Calculer Ne à l’aide de l’Équation 2-2

4. Calculer le nombre de Reynolds effectif en employant la relation de Ne vs. Re établie dans


l’étape 1.

5. Calculer la viscosité effective en utilisant l’équation


39

𝜌 ∗ 𝐷ℎ ∗ 𝑣
𝜇𝑒𝑓𝑓 = Équation 3-7
𝑅𝑒𝑒𝑓𝑓

Calculer Ks à l’aide de l’équation

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
𝑣

3.3 Procédure expérimentale de l’électrodialyse de la liqueur noire

3.3.1 Caractérisation de la liqueur noire

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.

[Link] Total des solides dissous

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

[Link] La masse volumique

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
𝑣𝑜𝑙𝑢𝑚𝑒

[Link] Détermination du teneur en lignine

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
𝐼

A est l'absorbance mesurée, I0 l'intensité de la lumière incidente, I l'intensité transmise, ε le


coefficient d'absorption, c la concentration de la sonde et d la longueur de la voie de la lumière.

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

Où D est la dilution, VD le volume de la LN diluée et V0 le volume de la LN d'origine.

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
𝑆

𝑋(𝑐𝑜𝑟𝑟. ) = 0.5 ∙ (𝑋1 + 𝑋2 ) Équation 3-17

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.

[Link] Caractérisation rhéologique

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

Tableau 3-4 - pH des échantillons préparés

E1 E2 E3 E4 E5 E6 E7

13.14 12.83 12.68 12.30 11.86 11.18 10.89

E8 E9 E10 E11 E12 E13 E14

10.67 10.29 10.05 9.85 9.52 9.16 9.01

3.3.2 Préparation des solutions pour l’électrodialyse

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.

[Link] Préparation de la liqueur noire

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.

[Link] Préparation de la solution NaOH

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

[Link] Préparation de la solution de rinçage

La solution de rinçage des électrodes est une solution 20% Na2SO4.

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

3.3.3 Le protocole suivi pour l’électrodialyse

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

3- Amorcer les pompes manuellement

4- Calibrer la perte de charge à zéro

5- Démarrer les pompes selon le débit désiré

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

7- Charger 3.5 litres de chaque solution dans sa cuve respective

8- Mettre en marche les pompes selon le débit désiré

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
𝐼

Où V est le voltage anode-cathode.

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

4.1 Caractérisation hydrodynamique de la cellule


Les résultats de la caractérisation hydrodynamique sont illustrés pour chaque géométrie
séparément. Les résultats seront suivis par une discussion générale.

À noter que les calculs de l’erreur sont expliqués en détails dans l’Annexe B.

4.1.1 Géométrie I : sans grille

[Link] Perte de charge en fonction du débit

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)

Figure 4-1 - Perte de charge en fonction du débit pour la géométrie I


46

[Link] Courbe de puissance pour géométrie I : Newton en fonction de Re

À partir des résultats précédents, la courbe du nombre de Ne en fonction du nombre de Re a été


construite. La limite du régime laminaire a été établie par rapport à la qualité du coefficient du
détermination R2. À partir des données expérimentales, un modèle mathématique a été établi à
l’aide de la méthode des moindres carrés pondérés. Par la suite, les deux courbes ont été comparées
en ajoutant un point à la fois et en évaluant le coefficient de détermination R2 après l’addition de
chaque point. Il a été décidé que la valeur de Reynolds pour laquelle, le coefficient R2 atteint un
maximum et il commence à régresser après sera considérée Reynolds critique. À partir de la vitesse
correspondant au nombre de Reynolds critique, le nombre de Newton reste constant quelle que soit
la vitesse.

100
Ne

10

1
25 250 2500
Re

Figure 4-2 - Ne en fonction de Re pour la géométrie I

La relation entre le nombre de Ne et le nombre de Reynolds pour un régime laminaire pour la


géométrie I est présentée au tableau 4-1 ci-dessous. Cette relation a été obtenue par la méthode des
moindres carrées pondérées.
47

Tableau 4-1 - Relation Ne-Re pour géométrie I

Géométrie Relation R2

Sans grille Ne = 1620 * Re-1 + 4.9 89%

Le nombre de Reynolds critique pour la géométrie I est 500.

4.1.2 Détermination de Ks

La constante Ks a été déterminé en suivant les étapes dans la section 3.2.4

La figure ci-dessous présente les valeurs de Ks en fonction de Reeff pour la géométrie I.

1000
Ks

100

10
10 100 1000
Reeff

Figure 4-3 - Ks en fonction de Re pour la géométrie I

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

4.1.3 Géométrie II : avec grille #1

[Link] Perte de charge en fonction du débit

La géométrie II correspond au séparateur avec la grille #1.

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)

Figure 4-4 - Perte de charge en fonction du débit pour la géométrie II

[Link] Newton en fonction de Re

À partir des résultats précédents, la courbe de puissance Ne en fonction du nombre de Re a été


construite. La limite du régime laminaire a été établie par rapport à la qualité du coefficient du
détermination R2
49

100

Ne

10

1
25 250 2500
Re

Figure 4-5 - Ne en fonction de Re pour la géométrie II

La relation entre le nombre de Ne et le nombre de Reynolds pour un régime laminaire pour la


géométrie II est présentée au ci-dessous. Cette relation a été obtenue par la méthode des moindres
carrées pondérées.

Tableau 4-2 - Relation Ne-Re pour géométrie II

Géométrie Relation R2

Grille #1 Ne = 1923 * Re-1 + 4.76 80%

Le nombre de Reynolds critique pour la géométrie I a été déterminé à 600.

[Link] Détermination de Ks

La constante Ks a été déterminé en suivant les étapes dans la section 3.2.4

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

Figure 4-6 - Ks en fonction de Re pour la géométrie II

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

4.1.4 Géométrie avec grille #2

[Link] Perte de charge en fonction du débit

La géométrie II correspond au séparateur avec la grille #2.

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)

Figure 4-7 - Perte de charge en fonction du débit pour la géométrie III

[Link] Newton en fonction de Re

À partir des résultats précédents, la courbe de puissance Ne en fonction du nombre de Re a été


construite. La limite du régime laminaire a été établie par rapport à la qualité du coefficient du
détermination R2
52

100
Ne

10

1
25 250 2500
Re

Figure 4-8 - Ne en fonction de Re pour la géométrie III

La relation entre le nombre de Ne et le nombre de Reynolds pour un régime laminaire pour la


géométrie III est présentée au ci-dessous. Cette relation a été obtenue par la méthode des moindres
carrées pondérées.

Tableau 4-3 - Relation Ne-Re pour géométrie III

Géométrie Relation R2

Grille #2 Ne = 3094 * Re-1 + 7.22 89%

Le nombre de Reynolds critique pour la géométrie I a été déterminé à 700.

[Link] Détermination de Ks

La constante Ks a été déterminé en suivant les étapes dans la section 3.2.4

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

Figure 4-9 - Ks en fonction de Re pour la géométrie III

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

4.1.5 Analyse générale pour les trois géométries

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

Géométrie Relation KP KNE

Sans grille Ne = 1620 * Re-1 + 4.9 1620 4.9

Grille #1 Ne = 1923 * Re-1 + 4.76 1923 4.76

Grille #1 Ne = 3094 * Re-1 + 7.22 3094 7.22

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

Tableau 4-5 - Ks pour les trois géométries

Géométrie Ks

Sans grille 275

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.

Tableau 4-6 - Valeurs de Ks pour certains mélangeurs industriels

Mélangeur Ks Référence

SMX (Sulzer) 47 (Fradette, 1999)

Ruban Hélicoïdal double 39 (Minge, Luhong, Bin, & Yuguo, 2008)

Ruban Hélicoïdal 30 (Paul, Atiemo-Obeng, & Kresta, 2003)

Rushton 12 (Paul, Atiemo-Obeng, & Kresta, 2003)

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.

4.2 Résultats de la caractérisation de la liqueur noire

4.2.1 Total des solides dissous

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.

4.2.2 Masse volumique et teneur en lignine

La masse volumique et le teneur en lignine ont été déterminés selon les sections [Link] et [Link]

Les résultats sont illustrés dans le tableau ci-dessous

Tableau 4-7 - Masse volumique et teneur en lignine de la liqueur noire utilisée

Masse volumique Teneur en lignine

1108 kg/m3 37%

4.2.3 Comportement rhéologique de la liqueur noire

Selon les résultats de mesure de viscosité obtenus pour les 14 échantillons il a été trouvé que :

- La liqueur noire se comporte comme un fluide newtonien jusqu’à un pH de 10.5 environ

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

- À partir d’un pH de 10.5 et en descendant, la liqueur noire se comporte comme un fluide


rhéofluidifiant. Ceci pourrait être expliqué par le fait que les particules de lignine précipitées
vont avoir un impact sur la viscosité de la liqueur noire.

La figure ci-dessous montre les courbes de la contrainte de cisaillement en fonction du taux de


cisaillement pour certains pH inférieurs à 10.9. Les tests de la viscosité ont été effectués à
température ambiante (23°C)

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

4.3 Résultats de l’électrodialyse de la liqueur noire


Les conditions d’opération pour les expériences d’électrodialyse sont présentées dans le tableau ci-
dessous.

Tableau 4-8 - Paramètres d'opération

Paramètre Valeur (unité)

Température des solutions au départ de l’expérience 23°C

Courant 1.2 A

Densité du courant ~300 A.m-2


58

Les débits d’opération pour chaque géométrie ont été choisis comme suit :

- Débit du régime laminaire ∝ Reynolds critique – 100

- Débit du régime turbulent ∝ Reynolds critique + 50

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à.

Les débits suivants ont été choisis pour chaque géométrie :

Tableau 4-9 - Débits d'opération pour l'électrodialyse de la liqueur noire

Géométrie I Géométrie II Géométrie III

Débit laminaire 1.75 l/min 2.3 l/min 2.65 l/min

Débit turbulent 2.6 l/min 3.2 l/min -

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

CHAPITRE 5 CONCLUSION ET RECOMMANDATIONS

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.

Après la construction de la cellule, une caractérisation hydrodynamique a été réalisée dans


le but de mieux comprendre le comportement des fluides à l’intérieur de la cellule et pour trouver
une constante représentative de chaque géométrie (Ks).

Finalement, l’électrodialyse de la liqueur noire a été effectuée dans le but d’extraire la


lignine par acidification et examiner l’impact de l’hydrodynamique sur l’opération. Et ce en
monitorant la résistance électrique dans la cellule à travers le temps puisque la résistance électrique
est liée à l’encrassement de la lignine sur la membrane. Les résultats trouvés concordent avec
l’étude CFD qui a été effectuée par Enciso et al. . En effet, l’hydrodynamique de la cellule a un
impact sur la performance de l’opération, la cellule semble être plus performante lorsque la
turbulence est plus élevée, offrant un quart de plus en temps d’opération

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

- Relier la géométrie (Ks) à la performance de réduction de l’encrassement

- Examiner l’impact de la température conjointement avec la turbulence sur la performance


de la cellule

- Réaliser l’électrodialyse en ajoutant d’autres compartiments de liqueur noire et de NaOH

- 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.

• Brun, J. (1989). Procédés de séparation par membrane. 169-186. Paris: Masson.

• Cloutier, J., Champagne, P., & Labrecque, R. (2013). Electrolytic Treatment of Black Liquor and
Caustic Recovery. United States of America: PEERS conference.

• Cloutier, J.-N., & Labrecque, R. (2013). Canada Brevet n° WO2013003941 A1.

• 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.

• Diop, A. (2014). Extraction, dépolymérisation et valorisation de la lignine Kraft de la liqueur


noire. Trois-Rivières, Québec, Canada: UQTR.
64

• 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.

• Gavach, C. (1999, Mars 23-24). Principes de l'électrodialyse. Saint-Hyacinthe, Québec:


Symposium sur les applications alimentaires de l'électrodialyse.

• 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.

• Haddad, M. (2016). Towards Sustainable Electrochemical Acidification of Kraft Black Liquor


for Lignin Extraction: Proof of Concept, Control of Membrane Fouling and Yeild Enhancement.
École Polytechnique de Montréal.

• 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. (2015). Development and integration of acid precipitation based lignin


biorefineries in Kraft pulping mills. Montreal, Quebec, Canada: École Polytechnique de
Montréal.

• Kannangara, M., Marinova, M., Fradette, L., & Paris, J. (2015). Effect of Mixing Hydrodynamics
on The Particle And Filtration Properties of Precipated Lignin . Canada: Elsevier.

• Kawatra, S. K. (2006). Advances in Comminution. SME.

• Kemperman, A. J. (2000). Handbook on bipolar membrane technology. Twente University Press.

• 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.

• Larousse. (2017). Hydrodynamique - définition. [Link].

• Lenntech. (2017). Electrodialysis. Récupéré sur Lenntech:


[Link]

• 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.

• Quoc, A. (2004). Optimisation de l'électrodialyse avec membranes bipolaires du jus de pomme


opalescent. Impact sur la qualité. Laval, Québec, Canada: Université de Laval.

• 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.

• Savelin, P., & Gavle, B. (2014). Digester developments. Sweden: Valmet .

• 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.

• Walker, W. (2010). Improving recovery in reverse osmosis desalination of inland brackish.


Austin: University of Texas.

• Wilhelm, F. (2001). Bipolar Membrane Electrodialysis. Twente University Press.

• 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

ANNEXE A – FICHE TECHNIQUE DES POMPES UTILISÉES


69
70
71

ANNEXE B – DÉTAILS DU CALCUL DE L’ERREUR


72

ANNEXE B – DÉTAILS DU CALCUL DE L’ERREUR

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 −
𝑆𝑆𝑇

𝑆𝑆𝑇 = ∑(𝑦𝑖 − 𝑦̅)2

𝑆𝑆𝐸 = ∑(𝑦𝑖 − 𝑦̂𝑖 )2

Avec 𝑦𝑖 𝑒𝑡 𝑦̂𝑖 étant les valeurs originales et modélisées

Vous aimerez peut-être aussi