Comportement de corrosion dans centrales REP
Comportement de corrosion dans centrales REP
THESE
Présentée par
Bruna MARTIN CABAÑAS
1
REMERCIEMENTS
Je remercie tout d’abord messieurs Christophe Varé, Hubert Catalette et madame Ellen
Mary Pavageau son successeur pour m’avoir accueilli dans le groupe T29 au sein du
département MMC (Matériaux et Mécanique des Composants) sur le site d’EDF - Les
Renardières. Ainsi que messieurs Michel Cassir et son prédécesseur Daniel Lincot pour
leur accueil au LECIME (Laboratoire d’Electrochimie, Chimie des Interfaces et
Modélisation pour l’Energie) de l’ENSCP.
Je remercie Philippe Barboux, qui a accepté d’être mon directeur de thèse ainsi que
Gregory Lefèvre et Michel Fédoroff qui m’ont encadré durant la thèse. Merci à Stéphanie
Leclercq qui, par sa disponibilité, ses remarques avisées et ses nombreuses relectures
m’a permis d’avancer sans m’égarer.
Je tiens à remercier les membres du jury pour avoir accepté de juger ce travail.
J’exprime toute ma gratitude à Patricia Vigne pour son aide, son soutien et ses précieux
conseils tout au long de ces trois années. Merci pour ces petits riens qui font un grand
tout !
J’adresse un merci chaleureux à Fabien Houlé et Jia Bai pour leurs travaux de stage de fin
d’études et leur souhaite bonne chance pour la suite. Merci à Sophie Delaunay, ma
voisine de bureau, pour nos échanges divers et variés au court de nos thèses, je lui
souhaite bonne continuation.
J’adresse également ma reconnaissance à toutes les personnes qui m’ont permis, de près
ou de loin, de mener à bien ce travail, en particulier Johannes Lützenkirchen, Mathieu
Guingo, Philippe Le Bec, Michel Mahé, Carine Mansour, Philippe Vermaut, Sandrine
Zanna, aux personnes du laboratoire de chimie… ainsi qu’à toutes celles qui ont fait de
certains moments, des instants de partage : Anissa, Perrine, Mathilde, Laurence, Antoine,
Aurélie, Nina, Arnaud, Jean-Baptiste, Alexandre, Sacha, Bandiougou, Xavier, Manu,
Christophe, Isabelle, Véronique, Chaveli…
J’exprime toute mon affection à ma famille et à mes amis qui ont partagé avec moi ce
moment de vie.
Enfin, à mon mari pour m’avoir soutenue, écoutée, épaulée et supportée sans relâche
durant ces trois années… merci !!!
A Louis…
2
I. INTRODUCTION............................................................................................................ 7
II. LES PARTICULES COLLOÏDALES DANS LE CIRCUIT PRIMAIRE DES
CENTRALES A EAU SOUS PRESSION............................................................................ 11
II.1. PRINCIPE DE FONCTIONNEMENT DES CENTRALES REP .............................................. 11
II.2. CONDITIONNEMENT CHIMIQUE DU CIRCUIT PRIMAIRE ............................................... 13
II.2.1. Contrôle de la réaction nucléaire .................................................................... 13
II.2.2. Contrôle du pH................................................................................................. 13
II.2.3. Contrôle du redox............................................................................................. 15
II.3. CONTAMINATION DU CIRCUIT PRIMAIRE .................................................................... 16
II.3.1. Corrosion des matériaux métalliques .............................................................. 16
II.3.2. Produits de corrosion dans le circuit ............................................................... 17
II.3.2.1. Mode de formation ................................................................................... 17
II.3.2.2. Nature des produits de corrosion.............................................................. 19
II.3.3. Activation des produits de corrosion ............................................................... 20
II.3.4. Purification du fluide primaire ........................................................................ 22
II.3.5. Mécanismes régissant le comportement des produits de corrosion (activés ou
non) dans le fluide primaire ............................................................................................. 23
II.3.5.1. Solubilisation/précipitation ...................................................................... 23
II.3.5.2. Adhésion/dépôt......................................................................................... 23
II.3.5.3. Sorption .................................................................................................... 25
II.3.5.4. Récapitulatif ............................................................................................. 26
II.3.6. Comportement dans le circuit lors des phases de fonctionnement .................. 27
II.3.6.1. Lors du fonctionnement ........................................................................... 27
II.3.6.2. Lors de la mise à l’arrêt............................................................................ 27
II.4. REDUCTION DU TERME SOURCE ................................................................................. 28
II.4.1. Choix des particules à étudier.......................................................................... 29
II.4.2. Choix des ions à étudier ................................................................................... 29
II.5. OBJECTIF DE LA THESE .............................................................................................. 29
III. REACTIVITE AUX INTERFACES SOLIDE-LIQUIDE. BIBLIOGRAPHIE .. 31
III.1. LES SYSTEMES COLLOÏDAUX ................................................................................. 31
III.2. L’INTERFACE SOLIDE-LIQUIDE............................................................................... 31
III.3. CHARGE DE SURFACE ............................................................................................ 32
III.3.1. Propriétés acide-base des oxyhydroxydes et charge de surface ...................... 32
III.3.2. Modèles d’acido-basicité et de charge de surface ........................................... 33
III.3.2.1. Description du modèle 2-pK monosite..................................................... 33
III.3.2.2. Description du modèle 1-pK multisites.................................................... 34
III.3.2.3. Description du modèle 1-pK monosite..................................................... 35
III.3.2.4. Choix du modèle ...................................................................................... 36
III.4. STRUCTURE DE L’INTERFACE SOLIDE-LIQUIDE ...................................................... 36
III.4.1. Modèle de Helmholtz........................................................................................ 36
III.4.2. Modèle de Gouy et Chapman ........................................................................... 37
III.4.3. Modèle de Stern................................................................................................ 38
3
III.4.4. Modèle de Grahame ......................................................................................... 39
III.4.5. Potentiel zêta .................................................................................................... 40
III.5. PARAMETRES CARACTERISANT L’ACIDO-BASICITE DE SURFACE ............................ 42
III.6. SORPTION .............................................................................................................. 43
III.6.1. Sorption par échange d’ions ............................................................................ 44
III.6.2. Sorption par complexation de surface ............................................................. 44
III.6.3. Sorption par précipitation/dissolution ............................................................. 45
III.6.4. Sorption par précipitation de surface .............................................................. 45
III.7. ETUDE BIBLIOGRAPHIQUE SUR LA SPECIATION DU NICKEL ET DU COBALT EN
SOLUTION ET DE LA SORPTION DE LEURS IONS SUR LES OXYDES ............................................ 46
III.7.1. Spéciation du nickel et du cobalt en solution ................................................... 46
III.7.2. Sorption d’ions sur les oxydes.......................................................................... 48
III.7.2.1. Sorption du cobalt .................................................................................... 49
III.7.2.2. Sorption du nickel .................................................................................... 53
III.7.3. Conclusion générale sur la sorption du cobalt et du nickel............................. 55
III.8. INFLUENCE DE LA TEMPERATURE SUR LES PROPRIETES ACIDE-BASE ET LA CHARGE
DE SURFACE ........................................................................................................................... 56
III.8.1. Influence de la température sur la charge de surface...................................... 56
III.8.2. Influence de la température sur la sorption ..................................................... 59
III.9. CONCLUSION ......................................................................................................... 63
IV. MATERIAUX ET METHODES .............................................................................. 65
IV.1. MATERIAUX .......................................................................................................... 65
IV.2. CARACTERISATION DES POUDRES .......................................................................... 65
IV.2.1. Diffraction des rayons X .................................................................................. 65
IV.2.2. Détermination de la surface spécifique par la méthode Brunauer, Emmet et
Teller (BET)...................................................................................................................... 66
IV.2.3. Taille et morphologie des particules ................................................................ 66
IV.2.3.1. Microscope électronique à balayage (MEB) ............................................ 66
IV.2.3.2. Diffraction des rayons X (DRX) .............................................................. 68
IV.2.3.3. Microscope électronique en transmission (MET) .................................... 69
IV.2.3.4. Diffusion de la lumière (DLS) ................................................................. 70
IV.2.3.5. Récapitulatif des tailles obtenues par les différentes méthodes ............... 71
IV.2.4. Impuretés présentes dans les poudres .............................................................. 71
IV.2.4.1. ICP-AES et ICP-MS................................................................................. 72
IV.2.4.2. Lavage des poudres .................................................................................. 74
IV.3. METHODES ............................................................................................................ 74
IV.3.1. Détermination des propriétés acide-base ........................................................ 74
IV.3.1.1. Titrage de masse à 25°C........................................................................... 74
IV.3.1.2. Zétamétrie................................................................................................. 74
IV.3.1.3. Titrage potentiométrique .......................................................................... 77
IV.[Link]. Principe ............................................................................................... 77
IV.[Link]. Hydratation des poudres...................................................................... 78
IV.[Link]. Conditions expérimentales .................................................................. 78
4
IV.3.1.4. Modélisation des courbes de titrage ......................................................... 78
IV.3.2. Sorption ............................................................................................................ 80
IV.3.2.1. Conditions expérimentales ....................................................................... 80
IV.3.2.2. Modélisation des courbes de sorption ...................................................... 81
IV.3.3. XPS ................................................................................................................... 82
IV.3.3.1. Principe de fonctionnement et appareillage ............................................. 82
IV.3.3.2. Echantillons d’oxydes témoins................................................................. 83
IV.4. POTENTIEL D’ECOULEMENT................................................................................... 84
IV.4.1. Appareillage et principe de fonctionnement .................................................... 84
IV.4.2. Conditions expérimentales ............................................................................... 86
V. ETUDE DES MECANISMES DE SORPTION DANS LE CIRCUIT PRIMAIRE
DES CENTRALES REP........................................................................................................ 87
V.1. SYSTEME NI2+ - FERRITE DE COBALT ......................................................................... 87
V.1.1. Cinétique de sorption ....................................................................................... 87
V.1.2. Influence du pH sur la quantité de nickel sorbée............................................. 87
V.1.3. Détermination de la nature des processus de sorption .................................... 89
V.1.3.1. XPS........................................................................................................... 89
V.1.3.2. Influence du nickel sur le PIE de CoFe2O4 .............................................. 91
V.1.4. Modélisation des courbes de sorption.............................................................. 92
V.1.4.1. Choix des espèces..................................................................................... 92
V.1.4.2. Calcul des constantes ............................................................................... 94
V.1.5. Conclusion........................................................................................................ 97
V.2. SYSTEME CO2+ - FERRITE DE NICKEL ......................................................................... 98
V.2.1. Cinétique de sorption ....................................................................................... 98
V.2.2. Influence du pH sur la quantité de cobalt sorbée ............................................ 98
V.2.3. Détermination de la nature des processus de sorption .................................. 100
V.2.3.1. XPS......................................................................................................... 100
V.2.3.2. Influence du cobalt sur le PIE de NiFe2O4 ............................................. 100
V.2.4. Modélisation des courbes de sorption............................................................ 102
V.2.5. Conclusion...................................................................................................... 104
V.3. SYSTEMES CO2+ ET NI2+ - MAGNETITE ..................................................................... 105
V.3.1. Cinétiques de sorption.................................................................................... 105
V.3.2. Influence du pH sur les quantités de cobalt et de nickel sorbées................... 106
V.3.3. Influence du cobalt et du nickel sur le PIE de la magnétite........................... 109
V.3.4. Modélisation des courbes de sorption............................................................ 111
V.4. MESURES A 50°C .................................................................................................... 115
V.4.1. Ferrite de cobalt............................................................................................. 115
V.4.2. Ferrite de nickel ............................................................................................. 116
V.4.3. Magnétite........................................................................................................ 118
V.5. MESURES A PLUS HAUTE TEMPERATURE.................................................................. 120
V.6. COMPARAISON DES TROIS SYSTEMES ....................................................................... 121
V.7. APPLICATION AU CIRCUIT PRIMAIRE ........................................................................ 122
5
VI. COMPORTEMENT DES PARTICULES COLLOÏDALES DANS LE CIRCUIT
PRIMAIRE DES CENTRALES REP ................................................................................ 125
VI.1. THEORIE DLVO ET IMPORTANCE DE L’HYDRAULIQUE........................................ 125
VI.2. PRESENTATION DES CALCULS AVEC LE MODULE LAGRANGIEN ........................... 128
VI.3. DETERMINATION DES POTENTIELS DE SURFACE DES MATERIAUX MASSIFS .......... 130
VI.3.1. Etude bibliographique.................................................................................... 130
VI.3.2. Résultats expérimentaux................................................................................. 133
VI.3.2.1. Mesure du potentiel d’écoulement ......................................................... 134
VI.3.2.2. Adhésion de particules de latex chargées............................................... 135
VI.4. CALCULS EFFECTUES A PARTIR DU MODULE LAGRANGIEN .................................. 137
VI.4.1. Récapitulatif des différents PIE...................................................................... 138
VI.4.2. Résultats des calculs....................................................................................... 138
VI.5. APPLICATION AU CIRCUIT PRIMAIRE .................................................................... 140
VII. CONCLUSION GENERALE ET PERSPECTIVES............................................ 147
ANNEXES............................................................................................................................. 151
REFERENCES ..................................................................................................................... 165
LISTE DES FIGURES......................................................................................................... 175
LISTE DES TABLEAUX .................................................................................................... 183
6
I. Introduction
La dégradation des matériaux constitutifs du circuit primaire des centrales nucléaires à eau
sous pression (REP), bien que faible, engendre la formation de produits de corrosion (PC)
sous forme d’ions et de particules colloïdales. Deux problématiques majeures se posent alors,
la modification du flux neutronique (AOA : Axial Offset Anomaly) (Hawkes-2004) dans le
cas de leur dépôt sur les gaines du combustible, et l’augmentation de la dosimétrie du circuit
suite à leur activation lors de leur passage dans le cœur, puis leur dépôt dans certaines parties
du circuit.
En effet, celles-ci ont la particularité de n’être ni retenues par les cartouches filtrantes, car leur
diamètre est trop faible, ni par les résines échangeuses d’ions, puisque, malgré leur charge de
surface, leur taille, leur encombrement stérique et la vitesse de passage sur les lits de résine
sont des conditions défavorables à cet échange. Les principaux processus, dont une meilleure
connaissance permettrait de réduire la contamination des centrales nucléaires à eau sous
pression, sont les phénomènes de corrosion, la formation des particules colloïdales, le dépôt
de ces particules sur les différentes parties du circuit primaire, et les interactions entre les
espèces ioniques et les oxydes (particules et dépôts). La charge de surface, des particules et
des parois du circuit, et les interactions électrostatiques sont des facteurs importants
contrôlant le dépôt de particules. C’est pourquoi, au cours d’un travail précédent (Barale-
2006), l’évolution de la charge de surface en fonction du pH a été mesurée entre 5°C et 350°C
pour trois types de particules, prises comme modèles représentatifs des produits de
corrosion présents dans le circuit primaire : le ferrite de cobalt, le ferrite de nickel et la
magnétite.
Dans le cadre de cette thèse, il a été décidé de continuer à travailler sur des particules
modèles, dans la mesure où des prélèvements in-situ de produits de corrosion ont rarement été
effectués. Ainsi les solides choisis sont des oxydes synthétiques de haute pureté et de
granulométrie faible.
7
L’étude des interactions particules/ions et particules/parois du circuit primaire constitue le
cœur de ce travail. Nous nous sommes d’abord intéressés aux interactions possibles entre les
oxydes sous forme particulaire ou déposés et les ions présents dans le fluide primaire, à partir
de l’étude des mécanismes régissant la sorption (évolution de la sorption en fonction du pH,
zétamétrie). La sorption d’ions a deux conséquences potentielles : la contamination des
oxydes de surface par des ions radioactifs et la modification de la charge de surface des
particules, ce qui peut faire évoluer leur taux de dépôt. Le processus de sorption sur ces
oxydes ayant été peu étudié, il est nécessaire de déterminer son importance et d’extraire les
constantes le caractérisant. Les ions étudiés ont été choisis du fait de leur impact sur la
dosimétrie : les isotopes radioactifs 58Co et 60Co, responsables respectivement d’environ 30 %
et 50 % de la radioactivité totale (Guinard-2003), proviennent respectivement de l’activation
du 58Ni et 59Co. Pour ces raisons, ce sont les cations du nickel et du cobalt qui seront étudiés
dans ce travail.
Etant donné que le dépôt des produits de corrosion peut avoir lieu sur les parois du circuit
primaire, les paramètres de surface de matériaux métalliques représentatifs du circuit (Inconel
690, Zircaloy 4 et acier inoxydable 304L) doivent être déterminés. De cette façon les
paramètres électrostatiques, à la fois des particules et des parois sur lesquelles elles vont se
déposer, seront connus. Pour cela deux méthodes sont utilisées : la zétamétrie (potentiel
d’écoulement) et la mesure de la vitesse d’adhésion de particules de latex chargées sur la
surface.
Dans la première partie de cet exposé seront donnés le contexte de l’étude, un rappel sur la
réactivité à l’interface solide-liquide ainsi qu’une étude bibliographique sur les données
connues à ce jour. Ensuite, la caractérisation des oxydes retenus (ferrite de cobalt, ferrite de
nickel, magnétite), qui a été complétée, afin d’avoir la meilleure connaissance possible de
leurs propriétés physico-chimiques, sera présentée. La sorption des cations d’intérêt nucléaire
(Co2+, Ni2+) en fonction du pH et de leur concentration, réalisée à différentes températures
(25°C et 50°C) sera présentée par la suite. Dans cette même partie, sera détaillée la
modélisation des données de sorption acquises à 25°C, qui a permis l’obtention des constantes
de complexation et de solubilité des espèces formées lors du processus de sorption. Une
dernière partie sera consacrée à une tentative de description du comportement global des
particules dans le circuit primaire. Pour cela les résultats obtenus sur les mesures de charge de
surface des matériaux massifs constitutifs des éléments du circuit (Inconel 690, Zircaloy 4,
8
acier inoxydable 304L), seront donnés. Ensuite, les calculs et les premières tendances prenant
en compte l’hydraulique, pour prévoir le dépôt de particules sur les parois, seront donnés.
Enfin, dans la conclusion, une synthèse des résultats sera effectuée ainsi qu’un ensemble de
propositions d’études pour améliorer la connaissance du comportement des produits de
corrosion dans le circuit primaire.
9
10
II. Les particules colloïdales dans le circuit primaire des
centrales à eau sous pression
11
Le cœur du réacteur est composé de 157 assemblages combustible (pour une centrale REP
900 MWe). Ces assemblages sont constitués de 264 crayons composés chacun de 272
pastilles d’UO2 fritté (ou de MOX, mélange UO2/PuO2) (figure II - 3). L’uranium est enrichi
en 235U à environ 3,5 % et l’eau est utilisée comme fluide caloporteur et modérateur.
La réaction de fission nucléaire a lieu dans le circuit primaire. Celle-ci est contrôlée à la fois
par la chimie du fluide primaire (bore) et par les barres de commandes. L’eau du circuit
primaire est maintenue sous forme liquide à l’aide de pressuriseurs (figure II - 4). La chaleur
produite au niveau du cœur est transmise au circuit secondaire par échange thermique au
niveau des générateurs de vapeur (GV – figure II - 4). L’eau du circuit secondaire passe alors
sous forme vapeur actionnant ainsi la turbine puis l’alternateur créant alors du courant
électrique. Cette vapeur est ensuite recondensée par échange thermique avec le circuit de
refroidissement.
12
Au cours de son fonctionnement, différents cycles sont effectués par le réacteur nucléaire. Un
cycle est composé de différentes phases, le fonctionnement normal, la mise à l’arrêt et le
redémarrage. La physico-chimie du fluide primaire évolue lors des transitoires d’arrêt et de
redémarrage (Taunier-2008). Avant de préciser en quoi consistent ces évolutions, les
spécificités chimiques du milieu primaire vont être décrites.
Le circuit primaire d’une centrale REP est conditionné par trois constituants : l’acide borique
(H3BO3), la lithine (LiOH) et l’hydrogène (H2), chacun répondant à un besoin de
fonctionnement ou de conditionnement particulier.
10
5 B+ 01n→37 Li+ 24 He + γ Equation II - 1
La teneur en bore dans le circuit primaire est imposée par la neutronique. Celle-ci varie donc
afin que la réactivité du milieu reste constante. En effet, lorsque le combustible est neuf, la
concentration en acide borique introduite dans le fluide est maximale (2500 ppm environ). Au
cours du fonctionnement du réacteur le combustible s’appauvrit en 235U impliquant ainsi la
diminution de la concentration en acide borique du fluide primaire. La corrosion des divers
matériaux du circuit est alors infime. Toutefois en cas de perte d’étanchéité aux joints des
pompes primaire, une forte concentration en acide borique peut entraîner la corrosion d’aciers
non alliés du circuit primaire (Gras-1987).
II.2.2. Contrôle du pH
1
Notons que le pH de neutralité de l’eau à 300°C est de 5,7.
13
lithine a été choisie. Préalablement à son utilisation, le lithium est enrichi à 99,9 % en isotope
7
Li afin de limiter la formation de tritium par réaction de l’isotope 6Li avec les neutrons selon
l’équation II – 2.
6
3 Li+ 01n→24 He+13H Equation II - 2
Lors des phases de fonctionnement du réacteur, la chimie du fluide doit être comprise dans la
zone rose de la figure II – 5. Hors de celle-ci, il est possible, en fonction des domaines de
concentration, de rester en fonctionnement pendant une durée maximale bien définie selon les
spécifications techniques d’exploitation.
14
Lors des différentes phases de fonctionnement la chimie du fluide est soumise à de grandes
variations au niveau du redox et de la température, et à de plus faibles pour le pH (tableau II –
1).
Un arrêt à froid peut également être envisagé mais différentes stratégies sont proposées en ce
qui concerne le lithium et l’oxygène, selon l’objectif visé (contamination minimum,
dissolution des produits de corrosion) (Bachet-2008).
Le dernier paramètre à contrôler concerne les conditions redox du milieu. En effet, sous
l’effet des rayonnements, la radiolyse de l’eau (Brun-1994) conduit à la formation
d’hydrogène et de radicaux libres oxydants par la réaction suivante :
1
H 2 O ⇄ H 2 + O2 Equation II - 3
2
L’oxygène ainsi produit est particulièrement néfaste pour la corrosion des parties métalliques
du circuit. Afin de limiter la corrosion des matériaux du circuit par l’oxygène formé par
recombinaison des radicaux libres, un ajout d’hydrogène de 25 à 35 cm3/kg est alors effectué.
Cet ajout permet de réduire la teneur en oxygène dissous à une valeur inférieure à 0,1 ppm.
15
II.3. Contamination du circuit primaire
Dans une centrale à eau pressurisée, le gainage du combustible est la partie soumise aux plus
fortes contraintes (température, pression, rayonnement). Le Zircaloy 4, du fait de ses
propriétés particulièrement bien adaptées à ces contraintes (faible absorption des neutrons,
bonnes propriétés mécaniques et bonne tenue à la corrosion) (Saint Paul-1993, Ly-2009) a été
retenu. Les alliages de zirconium sont également utilisés pour les tubes guides ou les grilles
d’assemblage. Ces alliages subissent deux types de corrosion. Le premier type est la corrosion
sous contrainte du fait du changement de densité du combustible et de l’apparition de gaz de
fission tels que l’iode, pouvant mener à l’endommagement de la gaine. Le second type est la
corrosion généralisée des surfaces en contact avec le fluide primaire.
16
Les tubes des générateurs de vapeur (GV), conçus en alliages à base nickel (Inconel 600 et
690), sont eux aussi soumis à la corrosion sous contrainte influencée notamment par la nature
du matériau, la température et les sollicitations mécaniques (Boudot-1992). L’alliage 600
étant fortement sujet à la corrosion sous contrainte et présentant une teneur en cobalt élevée,
l’alliage 690 lui a été préféré par la suite.
Les aciers inoxydables présents dans les tuyauteries et utilisés pour le revêtement de la cuve,
ainsi que les stellites (vannes et pompes) sont également soumis à la corrosion (Lister-1975a).
Globalement, les matériaux du circuit primaire ont été choisis car ils sont peu sensibles à la
corrosion. Ils se passivent en général spontanément et les couches d’oxyde formées sur les
Inconels ou les aciers inoxydables sont de l’ordre de la centaine de nanomètres et excèdent
rarement le micron. Néanmoins, compte tenu de la surface de la tuyauterie déployée qui est
très élevée (tableau II – 3), la quantité de produits de corrosion générée, bien que faible est
non négligeable. D’après les caractéristiques du tableau II – 3, il apparaît clairement que ce
sont les surfaces en Inconel qui présentent la plus grande surface (entre 15000 et 28000 m2),
et qui relâcheront le plus d’espèces, contribuant majoritairement à la contamination, selon un
mécanisme décrit plus loin. Ainsi la corrosion des tubes GV en Inconel conduit à la formation
de plusieurs dizaines de kilogrammes de produits de corrosion chaque année pour chaque
tranche.
17
Figure II - 6 : Vue schématique de la coupe de la couche d’oxyde formée à la surface
d’un Inconel 690 au contact du milieu primaire d’un REP (Laghoutaris-2009)
Des études menées par EDF R&D ont montré que deux régimes de relâchement se succèdent
(figure II – 7) lorsqu’un tube en Inconel est placé dans les conditions du primaire (325°C,
155 bar, 1000 ppm de bore, 2 ppm de lithium). Tout d’abord, un régime transitoire est observé
au cours duquel il y a un fort relâchement de nickel, puis dans un second temps, un régime
stationnaire est atteint, où le relâchement est plus faible, du fait de la formation d’une couche
passivante riche en chrome et de la diffusion lente du nickel au travers de celle-ci.
18
II.3.2.2. Nature des produits de corrosion
La composition des différents oxydes dépend du type d’alliage du circuit mais également de
la chimie du fluide (pH, concentration en H2). A partir de calculs d’équilibre chimique entre
le métal et les oxydes formés en surface, des estimations sur la quantité et la forme des oxydes
présents ont été réalisées (Bachet-2008). Les résultats sont donnés dans le tableau II - 4.
L’élément majoritaire selon ces calculs est le nickel et il peut se présenter sous différentes
formes, Ni métallique, NiO ou encore ferrites ou chromites qui ont une structure de type
spinelle (Annexe 1). Une étude menée par Kanbe et al. à partir d’une boucle d’essai, mène à
la même conclusion quant à la présence d’espèces à base de chrome, de fer et de nickel en
majorité (Kanbe-1983).
Expérimentalement, peu de prélèvements d’oxydes ont été effectués in situ. Une synthèse des
résultats disponibles, pour le parc nucléaire français, a été réalisée (Cattant-2008), et indique
qu’à la surface des alliages 600 la composition de l’oxyde est d’environ 50 % d’oxyde de
chrome, 30 % d’oxyde de fer (en majorité sous forme de magnétite) et 20 % d’oxyde de
nickel. Les éléments Co et Zn sont présents en très faible quantité (environ 2 %). La
composition chimique des différents oxydes n’est pas donnée précisément car les analyses ont
été réalisées par ICP-AES, ce qui permet uniquement de connaître les quantités de chaque
élément. Une étude menée aux Etats Unis sur la composition des couches d’oxyde présentes
sur les tubes de générateur de vapeur (Parsi-2009) indique que celles-ci sont formées de
Ni0,9Cr1,0 Fe1,1O4 . Ce qui est en accord avec la figure II – 6.
Concernant les espèces en solution, les tableaux II – 5 et II – 6 résument les résultats obtenus
en terme de répartition entre les produits de corrosion solubles et insolubles définis comme
ayant un diamètre supérieur à 10 nm, selon une étude menée en 1996 (Strasser-1996).
19
Tableau II - 5 : Concentrations des produits de corrosion insolubles en fonction de leur
diamètre (Strasser-1996)
Diamètre/Espèce (ppb) Cr Mn Fe Ni Co
d > 3 µm 0,15 0,012 0,8 0,19 0,003
0,45 µm < d < 3 µm 0 0 ,008 0,4 0,03 0,0011
0,01 µm < d < 0,45 µm 0,003 0,004 0,7 0,02 0,0002
Somme 0,153 0,024 1,9 0,24 0,0043
Les principaux éléments présents en solution sont le fer, le nickel et le cobalt. L’érosion ou
l’abrasion de certains éléments sont à l’origine de la présence d’autres contaminants tels que
110m
Ag et 124Sb.
Les produits de corrosion relâchés dans le milieu primaire ne contribuent pas directement à la
contamination du circuit, mais ils engendrent les espèces radioactives à l’origine de cette
contamination (De-Regge-1988).
20
Figure II - 9 : Mécanismes de transferts des produits de corrosion sous et hors flux
A
X + n→ A+1X + γ Equation II - 4
A
X + n →A Y + p Equation II - 5
21
Tableau II - 7 : Principaux radionucléides présents dans le circuit primaire
58
Isotope d’origine Ni 59Co 58
Fe 54
Fe 50
Cr 109Ag 121Sb 123
Sb
58
Isotope radioactif Co 60Co 59
Fe 54Mn 51Cr 110mAg 122Sb 124
Sb
Type de désintégration β+, γ β-, γ β-, γ β+, γ β+, γ β-, β+, β-, β,γ
-
γ β+, γ
58 60 59 54 51 110
Isotope issu de la Fe Ni Co Cr V Cd, 122Te, 124
Te
110
désintégration de l’isotope Pd 122Sn
radioactif
Temps de demi-vie 71 j 5,3 a 45 j 313 j 28 j 249 j 2,7 j 60 j
Le tableau II - 7 montre clairement que le 58Ni relâché par les tubes GV va s’activer lors de
son passage sous flux en 58Co qui est l’un des radionucléides les plus contaminants avec le
60
Co, responsables respectivement de 30% et 50% des doses reçues par les travailleurs
(Berger-2002). La nocivité du 58Co est due à sa forte concentration dans le circuit primaire, du
fait de la quantité importante de 58Ni relâchée par les tubes GV. Quant à celle du 60Co, elle est
due à sa période de demi-vie d’environ 5,3 ans. Différents prélèvements ont été réalisés sur
sites et ont permis de déterminer une tendance sur la répartition soluble/particulaire de
l’activité. Celle en 60Co est majoritairement transportée sous forme soluble, tandis que celle
en 58Co l’est d’avantage sous forme particulaire (Guinard-1995). Toutefois des écarts
importants sont trouvés selon les lieux de prélèvement.
L’importance de la taille des particules dans le transfert d’activité a également été montrée
(Guinard-1995). Selon cette étude, la radioactivité est principalement transportée par les
petites particules. En effet, 50% de l’activité serait transportée par des particules comprises
entre 0,2 et 0,45 µm de diamètre représentant pourtant seulement 6 % en masse des produits
de corrosion observés. Il est alors important de s’intéresser aux particules dans cette gamme
de diamètre afin de bien comprendre leur comportement et de pouvoir réduire leur impact.
Les produits de corrosion étant présents en grande quantité dans le fluide, une purification est
nécessaire. Celle-ci est réalisée en continu lors du fonctionnement de la centrale ainsi que lors
des arrêts de tranche. Cela est effectué par la mise en dérivation d’une certaine quantité de
fluide (environ 3,5 %) dirigée vers le circuit de contrôle volumétrique et chimique (RCV).
Après diminution de la température entre 25°C et 45°C, l’eau passe successivement sur des
cartouches filtrantes permettant la rétention de 99,8 % des particules en suspension ayant un
diamètre supérieur au micron, puis à travers des résines échangeuses d’ions (Gressier-2008).
L’étude des particules ayant un diamètre inférieur au micron, couramment appelées particules
colloïdales, est importante car celles-ci ne sont ni retenues par les filtres, ni par les résines
échangeuses d’ions et peuvent alors circuler librement dans le circuit. A la suite de ces
différentes étapes, le fluide est re-conditionné puis réinjecté dans le circuit.
22
II.3.5. Mécanismes régissant le comportement des produits de corrosion
(activés ou non) dans le fluide primaire
Dans le fluide, les produits de corrosion peuvent exister sous trois formes : espèces ioniques,
particules ou agglomérats de particules. Différents mécanismes de formation ont lieu selon les
phases de fonctionnement. Tous ces produits de corrosion peuvent réagir entre eux ou avec
les parois du circuit primaire selon plusieurs mécanismes en fonction de leurs caractéristiques
physico-chimiques ou de la chimie du fluide.
II.3.5.1. Solubilisation/précipitation
Des calculs ont été réalisés (Bachet-2008) afin de déterminer d’une part la composition à
l’équilibre des oxydes formés sur les tubes GV, et d’autre part les phases solides constituant le
dépôt sur les gaines d’assemblage combustible en fonction de la température et de la teneur en
hydrogène dissous. Puis des calculs de solubilité des différentes phases ont mis en évidence
d’importantes variations avec la température, la concentration en hydrogène ou en oxygène
dissous et le pH.
Les parties du circuit riches en chrome sont peu affectées par la mise en arrêt ou le démarrage
tandis que celles riches en nickel, fer ou cobalt le sont plus fortement.
II.3.5.2. Adhésion/dépôt
L’adhésion des produits de corrosion sur les parois du circuit est reliée à leur charge
respective, comme représenté schématiquement sur la figure II – 10.
23
Répulsion Attraction Répulsion
PIE 1
PIE 1 PIE 1
La formation de bulles à la surface des matériaux est un autre phénomène pouvant expliquer
le dépôt sur les parois des produits de corrosion contenus dans le fluide (Jones-2008). Si
l’hypothèse que les particules en suspension sont hydrophobes est faite, elles auront donc
tendance à se disposer à l’interface fluide primaire/bulle. Lorsque le diamètre de la bulle
augmente, celle-ci se détache progressivement de la surface. Les particules alors en contact
avec le matériau, et dans le cas où leurs charges sont opposées, auront une préférence pour la
surface ce qui formera ainsi un dépôt. Des expérimentations menées avec de l’hématite et de
l’aluminium montrent qu’entre pH 7 et 9 les conditions sont favorables à l’attraction.
24
Figure II - 11 : Formation des bulles induisant le dépôt de produits de corrosion (Jones-
2008)
Les données connues sur les mécanismes d’adhésion et de dépôt des produits de corrosion
modèles seront complétées lors de cette étude, afin d’approfondir les connaissances de ces
phénomènes.
Il faut noter que l’hydraulique joue un rôle très important sur ce mécanisme. Son influence sur
l’adhésion et le dépôt sera donc étudiée dans cette thèse.
II.3.5.3. Sorption
Le troisième mécanisme qui régit le comportement des espèces dans le fluide est la sorption.
Ce phénomène comprend la complexation et la précipitation de surface. Dans le circuit
primaire, peu de données sont disponibles sur la sorption d’ions tels que le cobalt, le fer ou le
nickel sur les produits de corrosion. Cependant plusieurs études ont été menées sur celle de
l’acide borique sur la magnétite et le ferrite de nickel, résumées par Guzonas (Guzonas-2002).
Mansour a regardé la sorption des sulfates sur la magnétite dans le circuit secondaire
(Mansour-2007). Dans le cadre du stockage des déchets radioactifs, la sorption du strontium
et de l’europium (Catalette-2004) ainsi que du lanthane et de l’ytterbium (Marmier-1994) sur
plusieurs oxydes de fer a été étudiée.
Ce phénomène dépend fortement du pH. Sur la figure II – 12, sont représentées la sorption du
cobalt et du nickel sur la magnétite. Les fronts de sorption de chacune des espèces ne se font
pas pour un même pH. Il est alors possible, que lors des modifications du pH du fluide, la
sorption ou la désorption d’une espèce, radioactive ou non, ait lieu.
25
Ni/Fe3O4 Co/Fe3O4
II.3.5.4. Récapitulatif
26
II.3.6. Comportement dans le circuit lors des phases de fonctionnement
Cette dernière étape doit être contrôlée car elle peut conduire à une augmentation de la
contamination du fluide primaire et elle est réalisée à différents moments de la mise à l’arrêt
en fonction de l’objectif visé :
- Si on se fixe comme objectif de dissoudre un maximum de ferrite de nickel pour éviter son
activation sous flux par la suite, il convient de maintenir une concentration en hydrogène
supérieure à 10 [Link]-1 le plus longtemps possible. Par contre, il y a un risque de
re-contamination à l’oxygénation (ou éventuellement lors de la déshydrogénation chimique) si
une quantité importante de fer (II) dissous est présente et incorpore du cobalt en précipitant.
Cette re-contamination sera plus sensible là où la surface de contact entre l’eau et les parois
est importante, c’est-à-dire au niveau des tubes de générateurs de vapeur et des échangeurs.
- Si l’objectif est de minimiser le risque de re-contamination, il faut réduire au plus tôt
l’hydrogène à un niveau bas (de l’ordre de 3 [Link]-1 à 200°C), et si possible maintenir une
27
concentration en lithium élevée. Dans ce cas, l’activité sous flux libérée dans le circuit
primaire devrait être minimisée.
- Une solution intermédiaire consiste à réduire au plus tôt l’hydrogène pour éviter de
dissoudre une quantité importante de fer (II), mais sans maintenir une teneur en lithium
élevée.
Lors d’une mise à l’arrêt avec oxygénation, une forte activité volumique en 58Co est observée
(figure II – 14) lors des premiers cycles effectués après le remplacement des tubes de
générateurs de vapeur.
12
10
Co (GBq/m2)
6
58
0
1 2 3 4 5 6 7 8 9 10 11 12 13
Au cours des cycles 1 et 2, une augmentation de la teneur en cobalt est visible, puis celle-ci
diminue pour atteindre un état stationnaire à des teneurs faibles à partir du septième cycle (à
l’exception d’une gamme de tube GV).
Des observations réalisées aux Etats Unis ont indiqué un comportement similaire de la teneur
en 58Co. Après le démarrage d’une tranche, le suivi de la concentration en 58Co lors des
différents cycles montre une forte augmentation puis une stagnation.
28
d’ions et de particules colloïdales pouvant être activés par le flux neutronique lors de leur
passage dans le cœur.
Malgré une étape de purification certains produits de corrosion, les particules colloïdales
notamment, ne sont pas récupérés et restent donc dans le fluide.
La problématique liée à la dosimétrie est de plus en plus contraignante malgré une forte baisse
ces dernières années. Celle-ci est passée de 2,44 [Link]/tranche en 1991 à 0,63 [Link]/tranche en
2007 (Ranchoux-2008). Pour abaisser la dosimétrie du parc nucléaire des améliorations sont
continuellement nécessaires. Diverses voies d’études ont été mises en place dans ce but :
- améliorer les outils de prédiction et d’optimisation de la contamination ;
- utiliser des matériaux avec des relâchements les plus faibles possibles ;
- augmenter le rendement de purification.
Après avoir travaillé sur ces divers points, il s’est avéré que les possibilités de gains
s’amenuisent et qu’il est nécessaire de s’intéresser à l’origine même de la contamination,
c’est-à-dire la création d’espèces radioactives ou susceptibles de le devenir ainsi qu’à leur
dépôt sur les parois hors flux du circuit.
Les produits de corrosions activés sont présents dans les dépôts sous forme de mélanges
Ni / NiO et d’oxydes mixtes de type spinelle, principalement des ferrites de nickel (Dacquait-
2002). La solubilité de ces dépôts évolue de plusieurs ordres de grandeur durant la mise à
l’arrêt, ce qui explique l’intérêt qui leur est porté.
Au vue des divers éléments résumés dans le paragraphe II.3.2.2 et d’études réalisées
précédemment (Barale-2006, Mansour-2007) le choix des particules s’est porté sur le ferrite
de cobalt, le ferrite de nickel et la magnétite.
Il a été vu que lors de la mise à l’arrêt du réacteur un pic de concentration en 58Co, provenant
de l’activation du 58Ni, apparaît pour les premiers cycles suivant le remplacement des
générateurs de vapeur. Pour cela, ce sont donc le nickel et le cobalt qui seront étudiés dans ce
travail.
29
nickel, ferrite de cobalt, magnétite). C’est en effet un mécanisme qui a été très peu étudié et
qu’il est important d’investiguer. Les phénomènes d’adhésion et de dépôt des produits de
corrosion au sein du circuit primaire ont été également étudiés afin d’élargir les connaissances
disponibles à ce jour, en prenant en compte en particulier les phénomènes hydrauliques.
L’acquisition de données physico-chimiques (jusqu’à 50°C) pour les produits de corrosion
classiquement rencontrés et pour les parois du circuit (Zircaloy, Inconel, acier inoxydable),
dans les conditions chimiques du milieu primaire a été effectuée.
Cette étude doit permettre, in-fine, par une meilleure compréhension des mécanismes
régissant le comportement des espèces dans le fluide primaire, de réduire le terme source.
30
III. Réactivité aux interfaces solide-liquide. Bibliographie
Divers processus traduisent la réactivité d’une surface au contact d’un liquide. Les termes
dissolution et désorption sont utilisés lorsqu’il s’agit d’un transfert d’espèces du solide vers le
liquide, et le terme sorption lorsque le transfert se fait du liquide vers le solide.
Toutefois il faut noter que la notion de sorption englobe différents mécanismes tels que la
précipitation ou la complexation de surface. Ceux-ci seront décrits par la suite.
L’étude des phénomènes de sorption aux interfaces solide-liquide a une place importante dans
le domaine de l’environnement. En effet, ces phénomènes jouent un rôle fondamental dans la
distribution des espèces chimiques entre les phases solides et liquides dans les milieux
naturels (lacs, mers, sols…). Ces phénomènes affectent particulièrement la réactivité des
surfaces et les propriétés électrostatiques des particules en suspension (Stumm-1970, Stumm-
1991). Dans le cadre du stockage de déchets nucléaires, la connaissance des mécanismes de
sorption permet notamment une meilleure compréhension des phénomènes de rétention et de
migration des éléments radioactifs en milieu naturel (Cremel-2007, Song-2009, Tertre-2005).
Ces phénomènes entrent également en jeu dans les mécanismes d’adhésion des particules sur
la surface, comme par exemple la formation de dépôts dans les circuits des centrales
nucléaires à réacteur à eau sous pression (Barale-2008, Sawicki-2008). La connaissance des
phénomènes de sorption à la surface d’un matériau contribue également à mieux comprendre
les processus de corrosion en milieux aqueux. En effet, la présence de certains éléments en
solution dans les centrales REP peut favoriser la corrosion de certains matériaux du circuit
(Coriou-1962, Kritzer-1999).
Le terme colloïde désigne une dispersion stable de particules dans une phase de nature
différente, tel qu’une dimension au moins des particules est comprise entre 1 nm et 1 µm.
L’une des caractéristiques des suspensions colloïdales est que les forces de gravité deviennent
négligeables, ce qui conduit à une vitesse de décantation très faible, voire nulle, des
suspensions. Une autre caractéristique est leur grand rapport surface/volume, ce qui favorise
les interactions aux interfaces solide-liquide, telles que les interactions électrostatiques et les
sorptions d’espèces présentes dans la solution aqueuse.
Les oxydes mis en solution acquièrent une charge de surface du fait de leur ionisation au
contact de l’eau.
D’un point de vue macroscopique une différence de potentiel existe entre le solide et la
solution. Toutefois, la faible mobilité des ions présents dans le solide ou la constante de
solubilité de celui-ci, ne permettent pas forcément la migration des ions du solide vers le
31
liquide ou la dissolution de l’oxyde. Cependant, la différence de potentiel chimique de
l’oxygène peut être abaissée par adsorption de l’eau à la surface de l’oxyde ou dissociation
des molécules adsorbées.
D’un point de vue structural, l’organisation atomique de la surface d’un matériau, par
exemple une particule colloïdale dans une solution aqueuse, est différente de celle de son
cœur. Les atomes présents à la surface sont entourés de moins d’atomes que ceux du cœur et
peuvent alors réagir avec le milieu environnant et s’hydrater.
Des groupements fonctionnels peuvent ainsi se former à la surface d’un solide lors de sa mise
en solution (Conway-1977, Jolivet-1994). La figure III - 1 représente un modèle d’hydratation
de la surface de la silice.
32
Une charge globale finale positive favorisera la sorption d’anions tandis qu’une négative celle
de cations. Le pH pour lequel la charge globale est neutre est appelé point de charge nulle
(PCN ou PZC). Celui-ci sera décrit plus longuement par la suite.
L’ajout d’anions ou de cations en solution, peut modifier la charge de surface par l’apparition
de phénomènes de sorption. Cette charge de surface perturbe la distribution des ions en
solution en contact avec le solide, créant ainsi l’interface solide liquide (Hunter-2001). Le
terme sorption prend en compte plusieurs phénomènes comme la complexation de surface et
la précipitation de surface (figure III – 3).
La charge de surface des oxydes peut être décrite par un modèle de complexation de surface.
Ces derniers possèdent différents degrés de complexité. Les plus simples sont les modèles
2-pK monosite et 1-pK monosite tandis que le plus élaboré est le 1-pK multisites (MUSIC –
MUltiSItes Complexation). Ces différents modèles sont décrits dans les paragraphes suivants.
Ce modèle, introduit par Stumm et al. (Stumm-1970) suppose l’existence d’un seul type de
site à la surface du solide. Il suppose aussi que ce site a des propriétés amphotères et peut se
protoner ou se déprotoner. Il décrit la charge de surface comme provenant de la protonation et
de la déprotonation des sites hydroxyles (Hunter-1981, Hunter-2001, Kosmulski-2001) selon
les réactions suivantes :
33
où K+ et K- représentent les constantes d’équilibre acide-base et X un atome métallique de la
surface. Les réactions précédentes sont supposées réversibles.
Donc :
[≡ XOH 2+ ]
K+ = Equation III - 3
[≡ XOH ][ H + ]
[≡ XO − ][ H + ]
K− = Equation III - 4
[≡ XOH ]
Lorsque la surface se charge, une attraction ou une répulsion électrostatique s’exerce sur les
protons, selon l’équation de Boltzmann :
−eΨ0
[H ] = [H ]
+ +
solution
exp
kT
Equation III - 5
A la surface la charge est égale au bilan des concentrations des espèces chargées, soit dans ce
cas :
Au point de charge nulle (PZC) la somme des concentrations en sites chargés positivement est
égale à celle en sites chargés négativement et les charges électrostatiques s’annulent (Hunter-
2001). Par arrangement des équations III – 3, III – 4, III – 5 et III – 6 l’équation III – 7 est
obtenue.
(pK − + pK + )
1
PZC = Equation III - 7
2
34
domaines de pH usuels). Il prend en compte les hétérogénéités de la surface en considérant
différents types de groupements (Kosmulski-2001) ayant des constantes de protonation et
déprotonation différentes (Hiemstra-1989a, Hiemstra-1989b, Hiemstra-1991, Pivovarov-
1997).
Ce modèle suppose qu’un atome d’oxygène puisse être lié à l’atome métallique de la surface
par une, deux ou trois liaisons. L’affinité du proton avec ces trois types de sites est alors bien
différente. La charge de l’ion métallique dans le solide est partagée avec les ligands qui
l’entourent ce qui implique les réactions suivantes (Hiemstra-1991) :
ou
ayant pour constantes d’équilibre Kn,1 et Kn,2, et où n est égal au nombre d’ions métalliques
coordinnés à un même type de sites et υ la valence de la liaison, définie comme la charge de
l’atome (z) divisée par son nombre de coordination (N). Dans ce modèle, la charge d’un site
de surface peut être fractionnaire.
Le modèle 1-pK monosite est un cas particulier du modèle 1-pK multisite décrit
précédemment. L’hypothèse prise en compte est que le solide présente une phase
cristallographique unique et donc un seul type de groupement de surface (Bolt-1982).
Par exemple, pour un oxyde ayant un groupement de surface –OH mono-coordiné, l’équation
acido-basique est la suivante :
≡ XOH −1 2
+ H + ⇄ ≡ XOH 21 2 Equation III - 10
35
III.3.2.4. Choix du modèle
Le modèle MUSIC est le plus réaliste des trois présentés précédemment pour décrire la
répartition des charges à la surface d’un solide. Cependant du fait du grand nombre de
données nécessaires, et difficiles à obtenir expérimentalement, pour déterminer la nature des
faces cristallographiques présentes, il reste assez peu utilisé. Le modèle 2-pK monosite étant
plus largement utilisé que le 1-pK monosite, il a été choisi pour interpréter nos données.
Dans ce modèle, la double couche est considérée comme un condensateur plan (Hunter-1981,
Hunter-2001). La charge et le potentiel de surface sont reliés par l’équation :
ε 0ε r
C= (F.m-2) Equation III - 13
d
36
solide. La capacitance C est ajustée par calcul à partir de données expérimentales. En règle
générale, la mesure d’une capacitance est rarement effectuée car les solides sont peu
conducteurs.
Dans le modèle de Helmholtz, il est supposé que les ions restent figés sur un plan. Gouy et
Chapman ont modifié ce dernier en tenant compte de l’agitation thermique dispersant les ions
qui compensent la charge de surface. Ces ions sont alors localisés dans une couche appelée
couche diffuse. A l’intérieur de la couche diffuse, lorsque la distance à la surface diminue, la
concentration en contre-ions augmente et celle en co-ions baisse.
Dans ce modèle, tout comme dans le précédent l’électroneutralité de l’interface doit être
respectée, d’où la relation σ0 + σd = 0 où σd est la charge de la couche diffuse. σ0 est calculé
par la relation suivante :
ZF Ψ
σ 0 = ( 8RT εε 0 I ) sinh
12 0
Equation III - 14
2 RT
37
L’inverse de l’épaisseur de la double couche, aussi appelé paramètre de Debye et Hückel, κ
est défini par l’équation suivante :
2Z 2 F 2 I 1 2
κ = Equation III - 15
εε 0 RT
Ce modèle est également appelé modèle de la couche diffuse ou Diffuse Layer Model (DLM).
Il est souvent utilisé pour calculer le terme électrostatique à partir de la valeur de la force
ionique sans ajustement des données expérimentales. Toutefois, il ne tient pas compte de la
taille des ions qui limite leur accumulation à la surface en n’empêchant pas leur approche à
une distance plus faible que celle imposée par leur rayon.
Afin de créer ce modèle, Stern a combiné ceux de Helmholtz et Gouy et Chapman avec
comme supposition que les ions ont une taille finie. Cela implique donc que le déplacement
des ions est limité. Cette distance limite est appelée plan de Stern et notée β. Les ions ayant
une affinité avec la surface se positionnent au plan de Stern et compensent une fraction de la
charge de surface. L’autre fraction est compensée par les ions présents dans la couche diffuse.
En notant σβ la charge au plan de Stern, la relation imposant l’électroneutralité de l’interface
est donc la suivante : σ0 + (σd + σβ) = 0.
Ce modèle est également appelé Basic Stern Model (BSM) (Stern-1924, Westall-1980).
38
Figure III - 6 : Schéma de l’interface solide-liquide selon le modèle de Stern
(Ψ : potentiel, x : distance à la surface)
Dans ce modèle (Westall-1980, Grahame-1947) la couche diffuse est décalée par rapport au
plan de Stern. Il y a donc présence de trois couches. Les ions ayant une affinité avec la surface
se placent au plan interne d’Helmholtz (PIH ou plan de Stern), situé à une distance β de la
surface. La couche diffuse débute alors au plan externe d’Helmholtz (PEH), situé à une
distance d de la surface. La distance entre la surface et le plan externe d’Helmholtz est de
quelques nanomètres. Ainsi, les ions faiblement hydratés traversent l’interface pour se placer
au PIH plus facilement que les ions fortement hydratés, qui eux se placent au PEH.
Ce modèle est également appelé modèle de triple couche ou Triple Layer Model (TLM).
39
III.4.5. Potentiel zêta
Lorsqu’on applique un champ électrique à une suspension de particules dont la surface est
chargée, celles-ci vont être mises en mouvement dans une direction qui dépend du signe de la
charge. Au cours du déplacement d’une particule dans un liquide, une partie des contre-ions
est entraînée par la particule tandis qu’une autre est immobile dans le liquide. La limite entre
ces deux parties est appelée plan de cisaillement (Conway-1977).
Les deux zones délimitées par le plan de cisaillement ne respectent pas l’électroneutralité.
Toutefois elle l’est sur l’intégralité de l’interface solide-liquide puisque les contre-ions
compensent la charge de surface de la particule. La différence de potentiel entre le plan de
cisaillement et la solution est appelée potentiel zêta (potentiel ζ). Le potentiel ζ est, en valeur
absolue, toujours inférieur au potentiel de surface σ0.
40
Les phénomènes électrocinétiques sont des processus irréversibles, dans lesquels une force
électrique induit un flux mécanique (ou inversement), tangentiellement à la surface chargée.
Le potentiel zêta est calculé à partir de coefficients en relation entre ces forces et le flux qui
en résulte (Kosmulski-2009). Différentes méthodes peuvent être utilisées pour déterminer le
potentiel zêta : l’électrophorèse, l’électro-osmose, le potentiel d’écoulement, le potentiel de
sédimentation et les méthodes électro-acoustiques (Kosmulski-2009).
- Par électrophorèse, c’est la vitesse des particules dans un champ électrique qui est
déterminée. Le potentiel ζ n’est pas mesuré directement mais calculé à partir de l’équation de
Smoluchowski :
- Lorsque la phase solide est fixe et qu’un flux électrique induit la mise en mouvement du
fluide, on parle d’électro-osmose. Cette méthode peut être utilisée pour des particules
sédimentant rapidement ou étant dans une solution colorée. Toutefois les limites sont la force
ionique (devant être inférieure à 0,1 M) et la forte quantité de matière nécessaire pour
l’expérience.
- Pour la mesure par potentiel d’écoulement, la phase solide est fixe et c’est le liquide en
mouvement en contact avec le matériau qui crée un champ électrique. La différence de
potentiel, entre les points d’entrée et de sortie du fluide, au niveau du solide est mesurée.
Cette méthode sera décrite plus en détail par la suite.
- Les méthodes précédentes considèrent que les phénomènes étudiés sont indépendant du
temps. Les méthodes électro-acoustiques prennent en compte cette variable en considérant
qu’un champ électrique alternatif induit des vibrations sur les colloïdes et sur les ions de
l’électrolyte support, variables et périodiques dans le temps. Cette contribution du temps est
surtout significative à fortes forces ioniques et pour de faibles concentrations en solide.
41
Dans cette étude, la surface des particules sera étudiée par électrophorèse et celle des solides
représentatifs du circuit par potentiel d’écoulement.
Les propriétés acido-basiques d’un oxy-hydroxyde en suspension dans une solution aqueuse
sont le plus souvent déterminées par un titrage acide-base. Par la mesure du pH au cours
d’ajouts successifs d’un acide ou d’une base à la suspension, le nombre de protons fixés ou
relargués par le solide en fonction du pH sont alors déterminés. La charge de surface σ0, est
prise comme égale à ce nombre, rapporté à la surface du solide. Elle s’exprime en [Link]-2,
mol.L-1 ou en C.m-2. Dans le modèle 2-pK, la charge de surface correspond à la différence :
[≡ XOH 2+ ] − [≡ XO − ] .
Lorsque la charge de surface, σ0, est nulle, le pH correspondant est alors appelé Point de
Charge Nulle (PCN) ou Point of Zero Charge (PZC) (Hunter-2001, Parks-1965). Il s’agit en
général d’une charge de surface moyenne nulle, car, sauf cas particulier, au PZC, il y a égalité
de la concentration en sites positifs et négatifs. Toutefois, en raison d’une charge initiale non
nulle au départ du titrage ou à cause de l’adsorption d’ions autres que H+, le pH où le bilan
des protons fixés ou relargués est nul, appelé PZNPC (Point of Zero Net Proton Charge), peut
ne pas correspondre à une charge réellement nulle.
Une méthode pour déterminer le PZC est d’effectuer des titrages à plusieurs forces ioniques
différentes. La charge de surface dépend de la force ionique de la solution, sauf quand elle est
nulle. On observe en général que les courbes de titrage à forces ioniques différentes se
croisent en un point appelé PZSE (Point of Zero Salt Effect). Ce point peut être pris comme le
véritable PZC.
Différents paramètres acido-basiques doivent être connus pour l’utilisation des modèles
présentés précédemment. Il a été vu que, dans le modèle 2-pK monosite, les deux paramètres
caractérisant les propriétés acide-base, sont les deux constantes pK+ et pK-. Dans ce modèle le
PZC est égal à ½(pK+ + pK-). Dans le modèle 1-pK monosite, PZC = pK. Un autre paramètre
important est la densité de site Ns, c'est-à-dire la charge maximale que peut prendre la
surface, quand les sites actifs acide-base sont saturés. Les autres paramètres dépendent du
modèle électrostatique de l’interface.
42
- Le modèle CCM nécessite la connaissance de la valeur de la capacitance C pour une force
ionique donnée. Cependant, ce paramètre est difficile à déterminer expérimentalement et
plusieurs valeurs peuvent être trouvées dans la littérature (Hiemstra-1991, Westall-1980). De
plus ce modèle ne permet pas de prendre en compte les effets de la force ionique. Il est
utilisable pour de fortes forces ioniques et/ou de faibles valeurs de potentiel de surface.
- Le modèle DLM n’implique pas de paramètre électrostatique inconnu, mais il n’est valable
que pour des faibles forces ioniques, de faibles charges de surface, et ne tient pas compte de la
taille des ions.
L’ensemble des paramètres acide-base est généralement obtenu par modélisation des courbes
de titrage acide-base. La modélisation consiste à faire une simulation numérique de la courbe
de titrage en optimisant les valeurs des paramètres de l’oxyde et de la double couche, de façon
à avoir l’écart minimal par rapport à la courbe expérimentale.
Le tableau III - 1 indique les paramètres nécessaires à l’utilisation du modèle 2-pK monosite
et les différents modèles électrostatiques.
Ces différents modèles ont couramment été utilisés pour décrire l’adsorption d’ions sur des
surfaces d’oxydes (Dzombak-1990, Stumm-1987, Marmier-1994).
III.6. Sorption
Le terme « sorption » recouvre tous les processus de transfert d’une espèce dissoute dans une
solution sur la surface d’un solide en contact avec cette solution : lorsqu’il s’agit de la fixation
d’une particule sur une surface, on parle généralement « d’adhésion ». Plusieurs processus de
sorption sont possibles :
43
- L’échange ionique
- La précipitation/dissolution
- La complexation de surface
- La précipitation de surface
L’échange ionique consiste en un échange entre cations (ou anions) en solution et cations (ou
anions) fixés à la surface d’un solide (comme par exemple sur une résine).
Dans ce modèle, ce sont les substitutions d’ions dans le réseau cristallin qui impliquent une
charge de surface. Il n’y a pas de transfert de charge entre la solution et l’espèce solide.
Plusieurs modèles pour décrire les phénomènes présents dans les milieux naturels ont été
proposés (Trémillon-1965). L’échange ionique (cas d’un cation) peut être représenté par
l’équation III – 17 et l’échange de ligand par la III – 18.
Dans les milieux naturels, les phases solides rencontrées possèdent des groupes fonctionnels
ionisables (-OH, -COOH, -NH2 …). En milieu aqueux, les groupements fonctionnels se
comportent comme des ligands classiques en solution (Schindler-1987). Comme nous l’avons
vu, ils peuvent subir des réactions acido-basiques conduisant à la formation d’une charge de
surface par protonation/déprotonation, jouer le rôle de ligands vis-à-vis des cations
métalliques ou des anions. Ces réactions entraînent une modification de la charge de surface.
Les groupements de surface sont généralement notés ≡SOH ou ≡XOH, ≡S et ≡X représentant
un atome de la surface du solide. Ils peuvent représenter différents groupements comme par
exemple des sites hydroxylés tels que ≡AlOH d’oxydes d’aluminium, des groupements
≡(CO2)OH de carbonates métalliques ou encore des groupements carboxyliques ≡(CO)OH
d’une surface organique.
La sorption d’ions sur ces sites –OH peut amener la formation de complexes de sphère
externe (Rahnemaie-2006) ou de sphère interne. Dans le premier cas les interactions
électrostatiques entrent en jeu ; l’ion se place au niveau du plan de Stern (PIH) après avoir
perdu sa sphère d’hydratation. Dans le deuxième cas, les groupements –OH de surface
forment une liaison à caractère covalent avec l’ion. (figure III - 9).
La formation d’une liaison entre la surface et l’ion conduit à un complexe monodentate, la
formation de deux liaisons donne un complexe bidentate (figure III - 9). La complexation de
surface a souvent été utilisée pour décrire les phénomènes présents à la surface d’oxydes
(Dzombak-1990, Stumm-1987).
44
Figure III - 9 : Schéma de l’interface solide-liquide dans le cas d’un hydroxyde. Le plan
0 correspond à la surface du solide, le 1 au PIH et le 2 au PEH. Les grandes sphères
blanches représentent les atomes d’oxygène et les petites les atomes d’hydrogène
A cet équilibre correspond une constante de complexation Kc. On voit dans cet exemple,
comment la sorption par complexation de surface modifie la charge de surface, qui passe de la
neutralité à négative.
La précipitation de surface est un mécanisme par lequel une nouvelle phase solide se forme à
la surface du solide. (Farley-1985). Ce mécanisme peut être détecté en mesurant les
45
isothermes de sorption, qui représentent la quantité d’élément sorbée en fonction de la
concentration en ce même élément en solution. Lorsqu’il y a une précipitation de surface, les
isothermes ont une allure caractérisée par l’absence de limite au niveau de la quantité sorbée
(Dumont-Ravat-1998) :
Dans un premier temps la spéciation du nickel (Ni2+) et du cobalt (Co2+) en solution aqueuse
a été calculée à 25°C (figure III – 11) à l’aide du code JChess 2.0 (Van der Lee-2001) et de la
base de données EQ3/6 (Wolery-1992).
L’hydrolyse de Ni2+ ou Co2+ et la formation des divers complexes hydroxo peuvent être
généralisées par la réaction :
Pour chacun de ces deux éléments, six complexes solubles ont été pris en compte : MOH + ,
M (OH ) 2 aq , M (OH ) −3 , M (OH ) 24 − , M 2 (OH ) 3 + et M 4 (OH ) 44 + (Plyasunova-1998a, Plyasunova-
1998b), ainsi qu’un précipité : M (OH ) 2 .
46
la cinétique. L’hydroxyde étant plus stable en solution que l’oxyde jusqu’à une température
de 250°C (Lister-1975b, Pourbaix-1963), la spéciation du cobalt et du nickel a été réalisée en
prenant en compte la formation des hydroxydes au lieu de celle des oxydes.
La présence de Ni(OH ) 2 aq et Co(OH ) 2 aq est prédite à partir de pH 8 et celle des
hydroxydes entre pH 8,5 et 13,5 à partir d’une concentration en cobalt ou en nickel de 10-5 M.
A très haut pH (supérieur à 13,5) le modèle prédit la dissolution du précipité pour former
successivement M (OH ) −3 et M (OH ) 24 − .
47
Figure III - 11 : Spéciation du nickel et du cobalt en solution à différentes concentrations
en Ni2+ ou Co2+ (a-b : 10-6 M, c-d : 10-5 M, e-f : 10-4 M) à 25°C
La sorption d’ions sur les oxydes entre en jeu dans divers phénomènes comme par exemple la
migration et le transport d’éléments dans les milieux naturels (Lemarchand-2005, Peak-2003,
Ridley-1999, Sanyal-2000) ou dans les circuits des centrales nucléaires (Mansour-2007,
Mansour-2009, Missana-2003, Venkataramani-1985). Ces processus et leurs mécanismes ont
donc été étudiés pour un grand nombre d’éléments et de matériaux.
48
Le bore, qui est un nutriment nécessaire aux plantes en faible quantité (Peak-2003), se sorbe
facilement sur les acides humiques (Lemarchand-2005), la calcite (Sanyal-2000) ou encore
certains oxydes de fer et d’aluminium (Goldberg-1985, Goldberg-1993). Goldberg et al.
(Golberg-1985) ont étudié la sorption du bore sur l’hématite et la goethite et ont proposé la
formation du complexe ≡ XH 2 BO3 par la réaction représentée par l’équation III – 21, avec
pour constante logK = 5,51 ± 0,31.
L’étude de la sorption du plomb sur l’alumine gamma (Hohl-1976) a amené les auteurs à
proposer la formation de deux types de complexes, un monodentate positif et un bidentate
neutre (0,1 M en NaClO4) :
Dans le domaine du nucléaire, l’activation des produits de corrosion, explique l’intérêt porté à
ce phénomène. Missana et al. (Missana-2003), ont regardé la sorption de l’uranium VI sur la
magnétite. Une courbe en « cloche » est observée avec un maximum pour un pH compris
entre 5 et 8. Venkataramani et al. (Venkataramani-1985) se sont eux intéressés à l’influence
de l’irradiation γ sur la sorption de Li+, Cs+ et Sr2+ sur la magnétite et l’oxyde de titane.
Dans le domaine du nucléaire, et plus précisément dans le cadre de ce travail, les éléments
d’intérêt sont le cobalt et le nickel. Les différentes études menées sur ces éléments sont donc
plus largement décrites par la suite.
Tamura et al. ont étudié la sorption des ions cobalt sur la magnétite entre pH 5 et 8 (Tamura-
1983, Tamura-1997a). Deux types de complexes de surface formés par réactions des ions
cobalt avec les sites ≡FeOH ont été proposés :
49
Le front de sorption s’étire de pH 5,7 à pH 7,5 pour une concentration en cobalt de
1,7.10-6 M, de pH 5,7 à pH 8,0 pour une concentration de 1,7.10-5 M et de pH 6,2 à pH 7,8
pour celle de 1,7.10-4 M, la concentration en sel de fond (NaNO3) étant de 0,1 M (figure III -
13a). Un léger décalage du front de sorption vers les forts pH est visible lors de
l’augmentation de la concentration en sel de fond (figure III – 13b) et un décalage vers les
faibles pH lors de l’augmentation de la température (figure III - 13c).
50
a b
Figure III - 13 : a). Pourcentage de Co2+ sorbé sur la magnétite (m/V = 10 g/L) en
fonction du pH pour différentes concentrations en cobalt à 25°C (Tamura-1983). b).
Influence de la force ionique sur le pourcentage de Co2+ sorbé à 25°C. (concentration en
cobalt 1,7.10-7 M) (Tamura-1983) c). Pourcentage de Co2+ sorbé en fonction du pH
(25°C) pour une concentration en cobalt de 1,7.10-7 M et en sel de fond de 0,1 M à
différentes températures (Tamura-1983)
Tewari et al. (Tewari-1975) ont étudié la sorption du cobalt sur les oxydes tels que Al2O3 ,
MnO2 et TiO2 avec une concentration en cobalt de 4.10-4 M. La sorption du cobalt augmente
significativement entre pH 6 et pH 8. En revanche, lorsque le pH dépasse 8, celle-ci diminue
et est remplacée par la précipitation de Co(OH ) 2 . L’augmentation de la quantité sorbée avec
le pH et la température apparaît comme étant due à l’augmentation de l’hydrolyse de Co(II)
suivie de la précipitation de surface et de la formation d’hydroxyde de cobalt insoluble. A
plus haute température, Co(OH ) 2 réagit avec Al2O3 pour former CoAl2O4 . De façon
similaire, les spinelles de nickel et de cobalt peuvent être formés par adsorption du cobalt et
du nickel sur la magnétite à haute température.
51
Une étude de sorption du cobalt a été réalisée sur des nanoparticules de magnétite (Uheida-
2006). Celles-ci ont une surface spécifique dix fois plus grande et une densité de sites cinq
fois supérieure à celles de la magnétite particulaire. De pH 4 à pH 12, seul le phénomène de
sorption est observé, autrement dit, il n’y a pas de précipitation pour ces conditions
expérimentales (figure III - 14). La réaction générale est la suivante :
Il a également été montré que l’adsorption de Co 2+ sur différents oxydes de fer varie
lorsqu’une forte concentration (0,1 M) en ions Cl − est ajoutée à la solution (Misak-1996). Le
front de sorption peut être décalé d’environ une unité de pH (vers les pH acides) avec le seul
ajout de 0,1 M de NaCl dans la solution. De plus, la sorption de Co 2+ et Co(OH ) + sur la
magnétite peut être accompagnée de celle de l’hydroxyde de cobalt dissous (Motl-2003).
Des études menées par Towle et al. (Towle-1995) et Persson et al. (Persson-1995) par EXAFS
sur la sorption du cobalt à l’interface eau-alumine alpha (à pH 8,1), eau- ZnO et eau- ZnS ,
montrent la formation de complexes monodentate et bidentate ainsi que celle d’un précipité
d’hydroxyde de cobalt. Une conclusion similaire pour les complexes formés a été obtenue par
l’étude de la sorption du cobalt sur MnO2 et Fe2O3 (Tamura-1997b).
L’utilisation de NiO pour décontaminer des solutions aqueuses du cobalt est en cours d’étude
(Naeem-2009). Il a été démontré la sorption préférentielle de Co(OH ) + par rapport à Co 2+
avec pour mécanisme :
52
De plus, la sorption augmente avec la concentration en cobalt, la température ainsi que le pH.
Dans une étude, il a également été montré que l’espèce hydrolysée, Co(OH ) + , devient
significative en solution au-delà de pH 9 (Kragten-1978), ce qui est en désaccord avec les
prédictions JChess ou cette espèce n’est présente de pH 8 à pH 10 qu’en faible quantité.
Un front de sorption allant de pH 5,5 à pH 8 a été observé lors d’une étude de sorption des
ions nickel sur la magnétite (figure III - 15) (Marmier-1999). Différents modèles ont été
utilisés pour déterminer les constantes de réaction (tableau III - 2).
53
Réaction de complexation de surface :
Tableau III - 2 : Résumé des constantes de sorption de Ni2+ sur la magnétite calculées
avec différents modèles (Marmier-1999)
CCM DLM NEM
pK+ 4,3 ± 0,2 4,6 ± 0,2 3,5 ± 0,3
pK- 8,3 ± 0,2 8,2 ± 0,2 9,6 ± 0,3
+
pKa1 ( ≡ Fe(OH ) 2 ) 4,3 ± 0,2 4,6 ± 0,2 3,5 ± 0,3
−
pKa2 ( ≡ FeO ) 8,3 ± 0,2 8,2 ± 0,2 9,6 ± 0,3
−
pK ( ≡ FeONi (OH ) 2 ) 14,8 ± 0,1 14,7 ± 0,1 15,9 ± 0,1
(CCM : Constant Capacitance model. DLM : Double Layer model. NEM : NonElectrostatic
Model, Surface spécifique : 1,8 m.g-1 )
Tableau III - 3 : Paramètres utilisés pour la modélisation des courbes de sorption avec le
modèle TLM (Kosmulski-1997)
LogK≡SiO- -5,71
LogK≡SiONa -7,24
C1 2,2 F.m-2
C2 0,2 F.m-2
Densité de sites de surface (Ns) 6 [Link]-2
54
Figure III - 16 : Effet de la force ionique sur l’adsorption du nickel sur les nanotubes
d’oxyde de fer en fonction du pH ([Ni2+] = 6 mg/L, m/V = 0,75 g/L, T = 25 ± 2°C (Chen-
2009)
Dans le cas du nickel comme dans celui du cobalt, des complexes monodentate et bidentate
sont proposés lors de la sorption à la surface d’oxydes (Tamura-1997b : MnO2 et Fe2O3 ),
ainsi que la formation d’un précipité d’hydroxyde de nickel (Baumgarten-1997, Beukes-2000,
Chen-2009) à partir de pH 9.
Un grand nombre d’études sur la sorption du cobalt et du nickel à la surface d’oxydes sont
disponibles dans la littérature. Comme rapporté dans les paragraphes précédents, ce
phénomène est souvent accompagné d’un mécanisme de précipitation.
Les espèces proposées dans la majorité des études menées sont les hydroxydes en ce qui
concerne la précipitation et ≡ XOM + , ≡ XOMOH ou (≡ XO)2 M pour la complexation de
surface. Il peut être noté que les comportements en solution, en présence (ou non) d’oxydes,
de ces deux ions sont similaires.
55
III.8. Influence de la température sur les propriétés acide-base et la charge
de surface
De façon générale, les mesures en température ne sont pas évidentes à mettre en œuvre
expérimentalement, notamment à cause de la mesure du pH (Lvov-2003). Cependant, des
mesures de PZC en température ont été réalisées pendant la deuxième moitié du XXième siècle.
La méthode la plus couramment utilisée pour réaliser ces expériences, est le titrage acide-
base.
Les diverses études réalisées, donnent les PZC en température pour les oxydes tels que la
magnétite (Tewari-1971, Blesa-1984), la goethite (Fokkink-1989), les hydroxydes de nickel et
de cobalt (Tewari-1976), l’oxyde de nickel (Tewari-1976) ou l’oxyde de chrome (Tewari-
1976), jusqu’à des températures de 90°C. La tendance générale observée pour ces diverses
espèces est la diminution du PZC avec l’augmentation de la température. Ces PZC ont été
extrapolés (Schoonen-1994) jusqu’à 350°C mais avec un grand nombre d’incertitudes du fait
des données expérimentales indisponibles au-delà de 150°C.
Les PZC de la magnétite, obtenus entre 25°C et 90°C puis extrapolés à plus haute température
par divers auteurs (Marmier-1999, Parks-1965, Regazzoni-1983, Blesa-1984) ont été recensés
par Wesolowski (Wesolowski-2000). Plus récemment des mesures de PZC ont été réalisées
jusqu’à 290°C (Wesolowski-2000) et 320°C (Barale-2008). La figure III - 17, montre
l’évolution du PZC de la magnétite en fonction de la température pour différentes études
(Wesolowski-2000). La courbe représentant ½ pKe, où pKe est la constante d’ionisation de
l’eau est également donnée. Sur ces courbes, on observe la diminution du PZC avec
l’augmentation de la température (jusqu’à 250°C). Une autre étude a montré que pour une
température supérieure à 250°C, le PZC augmentait à nouveau (Barale-2006).
56
Une modélisation de la variation du PZC en fonction de la température a été proposée par
Schoonen (1994), d’après le modèle 2-pK avec comme équations :
H 2 O ⇄ OH − + H + Ke Equation III - 33
Par association de ces trois équations, l’équilibre chimique suivant est obtenu :
donc :
T
∆GT = ∆G298 − (T − 298)∆S 298 + ∆C p (T − 298 − T ln ) Equation III - 36
298
où ∆GT, ∆G298, ∆S298 et ∆Cp sont respectivement les variations d’enthalpie libre à la
température T, d’enthalpie libre et d’entropie à 298 K et de capacité calorifique. Les
paramètres de cette équation sont considérés constants entre 25°C et 350°C. D’après celle-ci,
il est prédit que le PZC diminue jusqu’à 200°C puis réaugmente entre 200°C et 350°C.
Le modèle 1-pK pour décrire l’évolution du PZC d’oxydes (magnétite, rutile, zircone) en
température a également été utilisé (Machesky-1998, Wesolowski-2000, Bérubé-1968,
Barale-2008). La relation obtenue est alors :
0 298 ∆C p −1 0 ∆C p (1 + ln298) ∆C p
PZC = − ∆H 298 − T + ∆S298 − + log T Equation III - 37
2,303R 2,303R R
Dans le cas du rutile, la différence entre le 1/2pKe et le PZC est constante quelque soit la
température, par contre pour la magnétite celle-ci augmente avec la température. D’après les
études faites jusqu’à 320°C (Barale-2008), le PZC devient plus élevé que ½ pKe au-dessus
d’environ 180°C pour le ferrite de cobalt et de 250°C pour la magnétite.
57
Des études à 235°C ont été effectuées (Jayaweera-1993, Jayaweera-1994) sur des oxydes tels
que TiO2 , Fe3O4 , ZrO2 ou Cr2O3 . Par exemple, pour la magnétite le PZC obtenu en
température est de 6,1 ce qui est donc plus faible que celui à température ambiante (6,5).
Toutefois l’effet de la température est plus ou moins marqué selon les oxydes. Dans la figure
III - 18, l’évolution du potentiel zêta de PdO , Fe2O3 , Fe3O4 et TiO2 en fonction du pH à
235°C est donnée.
Figure III - 18 : Evolution du potentiel zêta de PdO, Fe2O3, Fe3O4 et TiO2 à 235°C
(Jayaweera-1993, Jayaweera-1994)
A température ambiante un grand nombre de données sur les PIE de matériaux divers est
disponible. Celles-ci ont été répertoriées par Kosmulski (Kosmulski-2009). Les valeurs à
25°C et 235°C de certains oxydes sont reportées dans le tableau III - 4 :
58
III.8.2. Influence de la température sur la sorption
Ils concluent que la sorption du calcium a lieu pour des pH de plus en plus acides avec
l’augmentation de la température, mais que le processus de sorption de Ca 2+ ne varie pas en
fonction de celle-ci.
Tewari et Lee (Tewari-1975) ont réalisé l’étude de la sorption de Co 2+ sur NiFe2O4 , TiO2 et
ZrO2 . Ils ont montré que la saturation est obtenue pour de plus faibles concentrations en
cation à haute température. Le pH en température a été extrapolé à partir de celui à
température ambiante.
59
Figure III - 20 : Adsorption de Co(II) en fonction de sa concentration en solution à
l’équilibre (pH de référence à 25°C = 6,5) sur NiFe2O4 (Tewari-1975)
Bauer et al. (Bauer-2005) se sont intéressés à la sorption de l’europium (III) sur la smectite
entre 25°C et 80°C. Les trois courbes obtenues sont similaires mais compte tenu de la
formation d’une seconde phase à partir de 60°C, il n’est pas exclu que différentes réactions
entrent en jeu. De plus le pH indiqué en abscisse des courbes est celui mesuré à la température
de l’essai.
Tertre et al. ont étudié la sorption de l’Eu(III) sur la kaolinite (Tertre -2006) et celle de Ni(II),
Ln(III) et Eu(III) sur la montmorillonite (Tertre-2005, Tertre-2006) entre 25°C et 150°C. La
tendance observée dans ces études est une augmentation de la sorption avec la température
pour un même pH, avec un décalage du front de sorption vers les pH acides.
60
Figure III - 22 : Log Kd de l’europium sur la kaolinite à différentes températures
(Tertre-2006)
61
Figure III - 24 : Concentration en Ni(II) sorbé sur la montmorillonite à différentes
températures (pH = 5,8, I = 0,01 M, : 20°C, : 37°C, O : 67°C) (Song-2009)
L’influence de la température sur la sorption des anions a été moins étudiée que sur celle des
cations. Toutefois Mansour (Mansour-2007, Mansour-2010) a regardé l’influence de la
température pour le système sulfates - magnétite. A 25°C il se forme majoritairement des
complexes de sphère externe. Entre 50°C et 275°C (sans présence d’hydrogène) d’autres
formes de complexes (sphère interne monodentate et bidentate) sont présents.
62
III.9. Conclusion
Les différentes études de sorption de cations en température indiquent une même tendance, à
savoir la diminution du PZC avec l’augmentation de la température jusqu’à 250°C puis
l’inversement du phénomène pour des températures supérieures.
La sorption de cations à la surface d’un solide augmente avec la hausse du PZC en fonction
de la température. Les résultats obtenus par Mansour (Mansour-2007) montrant une
diminution de la sorption des anions avec l’augmentation de la température est cohérente avec
l’évolution du PZC.
63
64
IV. Matériaux et méthodes
Les matériaux utilisés au cours des expériences ont été choisis à partir des différentes
particules supposées présentes dans le fluide primaire (Taunier-2008). Il est en effet
impossible de travailler directement sur des particules issues du circuit primaire du fait de leur
forte radioactivité et de la difficulté de les extraire dans les conditions de fonctionnement du
réacteur. Il faut noter qu’afin de s’affranchir d’une éventuelle oxydation de la surface des
solides étudiés, toutes les expériences ont été réalisées en boite à gants sous atmosphère
inerte.
IV.1. Matériaux
Dans le but de bien comprendre les mécanismes entrant en jeu dans les phénomènes de
sorption, le choix des matériaux s’est porté sur les plus purs disponibles dans le commerce. Il
s’agit de cette façon, de s’affranchir au maximum de l’influence que peuvent avoir les
impuretés sur les phénomènes étudiés. Des études menées précédemment (Barale-2006,
Mansour-2007, Catalette-1998) sur les mêmes poudres (mêmes lots de fabrication) ont permis
de connaître les caractéristiques principales des produits utilisés. Les poudres choisies sont :
la magnétite ( Fe3O4 ), le ferrite de cobalt ( CoFe2O4 ) et le ferrite de nickel ( NiFe2O4 ) et
proviennent de chez Nanostructured and Amorphous Materials et Alfa Aesar (annexe 2).
Des échantillons massifs représentatifs de certains éléments du circuit primaire ont également
été utilisés : Inconel 690 (tubes de générateurs de vapeur, tuyauteries), Zircaloy-4 (gaines
combustibles) et acier inoxydable 304L (tuyauteries).
Dans ce chapitre, les méthodes d’analyses utilisées uniquement dans les études précédentes
(Barale-2006, Mansour-2007) seront décrites succinctement. Celles ayant été également
utilisée lors de ce travail seront plus longuement détaillées.
Les ferrites de cobalt et de nickel ainsi que la magnétite ont été analysés par diffraction des
rayons X (DRX) à l’aide d’un diffractomètre Rigaku Geigerflex ayant une anticathode en
cobalt. La radiation utilisée correspond à la radiation Kα1 du cobalt.
Cette méthode non destructive permet d’identifier les phases cristallines d’un matériau par
comparaison du diffractogramme obtenu avec des fichiers de référence : JCPDS (Joint
Committe on Powder Diffraction Standards). Les diagrammes de diffraction obtenus
65
correspondent bien aux phases attendues. Des traces de Co3O4 dans le ferrite de cobalt et
d’hématite dans la magnétite ont été détectées.
La méthode BET (Brunauer, Emmett et Teller) est utilisée pour déterminer la surface
spécifique de matériaux. Par mesure de la quantité de gaz (dans cette étude, l’azote) adsorbé à
leur surface et en faisant l’hypothèse que le gaz est physisorbé sur une monocouche, il est
possible de calculer l’aire spécifique selon la méthode mise au point par Brunauer, Emmet et
Teller (Brunauer-1938).
Avant analyse, les poudres ont été dégazées 24 h à 110°C.
A partir de la surface spécifique, on peut estimer la taille des particules, en faisant une
hypothèse sur leur forme. Pour des particules supposées sphériques, le diamètre d se calcule
par :
6
d= Equation IV - 1
Sρ
où S est la surface spécifique et r la masse spécifique. Les résultats sont reportés dans le
tableau IV - 2.
Il est possible de déterminer la taille de particules par plusieurs méthodes. Les résultats
obtenus peuvent toutefois être très différents en fonction de celle utilisée.
Les poudres ont été lavées préalablement aux analyses. Le protocole de lavage est décrit en
annexe 3.
La première méthode utilisée pour évaluer la taille et la forme des particules est la
microscopie électronique à balayage (MEB). Le principe du MEB consiste en un faisceau
d’électrons balayant la surface de l’échantillon à analyser qui, en réponse, réémet des
électrons secondaires permettant de déterminer la topographie de la surface.
66
Figure IV - 1 : Observations au MEB du ferrite de cobalt (a : Barale-2006)
67
Pour le ferrite de nickel, des agglomérats pouvant atteindre 20 µm et constitués de particules
de 10 à 50 nm sont observés.
La magnétite est constituée de particules polyédriques, de 150 à 500 nm (en accord avec
Philippini-2006), soudées les unes aux autres, formant ainsi des agglomérats de 10 à 50 µm.
0,9 λ
Dp = Equation IV - 2
L cosθ
68
IV.2.3.3. Microscope électronique en transmission (MET)
Un échantillon suffisamment mince est mis sous un faisceau d’électrons focalisés. En mode
image, les électrons qui traversent l’échantillon passent ensuite à travers un système de
lentilles magnétiques permettant de projeter l’image de l’échantillon sur un écran fluorescent
ou sur un capteur CCD, transformant alors l’image électronique en image optique (figure IV -
4).
Les électrons sont plus ou moins absorbés selon l’épaisseur, la densité ou la nature chimique
de l’échantillon.
Pour obtenir un cliché de diffraction, c’est le comportement ondulatoire des électrons qui est
utilisé. Lorsque ceux-ci rencontrent de la matière sous forme cristallisée, ils sont diffractés
dans différentes directions dépendantes de l’organisation atomique. L’image de diffraction
permet de déterminer la structure cristalline de la zone observée (figure IV - 4).
69
20 nm 50 nm
Les particules observées ont un diamètre qui varie entre 30 et 150 nm et appartiennent à des
agglomérats de 1 à 2 µm.
Les particules en suspension dans un liquide ne sont jamais immobiles. Elles sont soumises au
mouvement Brownien et entrent en collision les unes avec les autres ainsi qu’avec le milieu.
Le principe de la méthode DLS (Dynamic Light scattering) est le suivant : une source
lumineuse éclaire les particules (dans ce cas un faisceau laser) qui diffusent la lumière dans
toutes les directions ; une constante de diffusion est alors obtenue. Une caractéristique
importante à prendre en compte est que plus les particules sont petites, plus elles ont un
déplacement rapide. En utilisant l’équation d’Einstein-Stokes (équation IV - 3) qui relie la
constante de diffusion au rayon d’une particule, il est alors possible de mesurer sa taille.
L’hypothèse sous-jacente est cependant que les particules sont sphériques.
kB T
D= Equation IV - 3
6πηr
70
où D est la constante de diffusion, kB la constante de Boltzmann, T la température, η la
viscosité et r le rayon de la particule en mouvement.
Cette méthode donne une répartition des tailles avec un intervalle très large. L’inconvénient
principal est que le temps de passage aux ultrasons de la solution lors de sa préparation joue
un grand rôle sur le résultat obtenu.
Les tailles de particules obtenues par les méthodes décrites dans ce paragraphe sont bien du
même ordre de grandeur. Pour ces trois poudres toutes les analyses ont montré qu’il y a
formation d’agglomérats. Les résultats sont synthétisés dans le tableau IV - 2.
La disparité des résultats n’est pas étonnante car chaque technique ne rend pas compte de la
même manière de la distribution de taille et de l’état d’agrégation. Le MEB et la DLS ne
visualisent pas avec précision les objets les plus petits quand ils sont agrégés. Du fait de
l’agrégation, toute la surface des particules n’est pas accessible pour l’adsorption et
l’extrapolation du diamètre des particules à partir du BET donne une valeur largement
surdimensionnée. Au contraire la DRX donne une valeur des domaines de diffraction,
cristallites de base, mais ne rend pas compte de leur agglomération. D’autre part cette
moyenne (qui donne en fait une épaisseur moyenne dans une direction qui varie sur une
sphère), minore celle observée directement en MET.
En fait, même si l’état de dispersion n’est pas véritablement homogène, ce qui reste important
est la surface spécifique qui interviendra dans le calcul du nombre d’atomes adsorbés par
unité de surface.
Les matériaux, bien qu’acquis sous la forme la plus pure, contiennent quelques impuretés.
Celles-ci, en plus ou moins grande quantité, ont un rôle non négligeable sur les propriétés de
surface, notamment pour la détermination du point de charge nulle et du point isoélectrique.
La présence d’oxyde de cobalt dans le ferrite et d’hématite dans la magnétite a été exposée
lors de la caractérisation des oxydes par diffraction des rayons X. Une autre méthode peut
71
s’avérer utile pour la détermination des impuretés en présence, la spectrométrie d’émission
atomique.
Lors de la synthèse, des éléments tels que le sodium ou le calcium peuvent se retrouver piégés
dans le solide. Ceux-ci proviennent de la synthèse des poudres. L’analyse par DRX ne détecte
que des phases dont la concentration est supérieure à quelques %. La spectrométrie
d’émission atomique par plasma à couplage inductif (ICP-AES) permet, après dissolution du
matériau, d’obtenir une information sur les éléments présents.
La magnétite et les deux ferrites ont été dissous par une attaque acide (2/3 HCl, 1/3 HNO3
concentrés) avant d’être analysés par ICP-AES.
La technique d’ionisation des éléments par couplage inductif associé à un plasma (Inductively
Coupled Plasma) consiste à faire passer une solution aqueuse contenant le ou les éléments à
analyser à travers une torche plasma. Un plasma est un gaz ionisé, caractérisé par la présence
d’électrons libres ayant une température comprise entre 6000 et 10000K. Le plasma à
couplage inductif est formé à partir d’une décharge électrique qui excite un gaz plasmagène
(généralement de l’argon) circulant à travers une série de tubes de quartz (torche) entourés
d’une spire (bobine d’induction). Celle-ci est reliée à un générateur de radiofréquence. Le
passage de courant électrique dans la bobine d’induction permet la formation d’un champ
électromagnétique. Une étincelle se forme donc dans la torche qui est alimentée en argon, est
permet l’ionisation de celui-ci selon la réaction suivante :
Ar ⇄ Ar + + e − Equation IV - 4
Les électrons, accélérés par le champ magnétique, amplifient le phénomène jusqu’à ionisation
totale du gaz. Le liquide à analyser est alors ionisé à son tour dès son arrivée dans la torche.
Une analyse élémentaire est donc possible.
L’échantillon est prélevé à l’aide d’une pompe péristaltique. Il passe ensuite à travers un
nébuliseur lui permettant de passer sous forme d’aérosol constitué de gouttelettes de liquide
d’environ une dizaine de micromètres de diamètre. Ces gouttelettes entrent ensuite dans le
plasma et sont alors ionisées.
72
L’analyse des éléments est rendue quantitative par un étalonnage réalisé avant chaque série de
mesure avec des solutions étalons commerciales, ainsi que par des points de contrôle effectués
à partir de solutions commerciales différentes des précédentes.
Lorsque l’échantillon à analyser est présent en très faible quantité (< 10 ppb), l’ICP-MS
(Inductively Coupled Plasma – Mass Spectroscopy) est utilisé. Le système de torche est le
même que pour l’ICP-AES mais le système d’analyse est dans ce cas un spectromètre de
masse.
A la sortie de la torche, une partie du plasma est échantillonée par un premier orifice (1 mm
de diamètre) au sommet du premier cône (sampler), puis se détend sous l’effet du vide
(1-2 mbar). Ceci permet de passer de la pression atmosphérique au vide secondaire du
spectromètre de masse. L’échantillon voit également sa température réduite par ce système.
Le passage à travers un deuxième cône échantillonneur (skimmer) permet d’écrêter à nouveau
le signal puis celui-ci est focalisé par un faisceau de lentilles ioniques et dirigé vers un
spectromètre de masse quadripolaire.
Celui-ci joue le rôle d’un filtre en détectant les éléments en fonction de leur rapport
masse/charge. Seuls les éléments désirés atteignent alors le détecteur. Un multiplicateur
d’électrons permet d’augmenter l’intensité du signal et donc d’avoir une très grande
sensibilité.
Tout comme pour l’ICP-AES, l’analyse des éléments est rendue quantitative par un
étalonnage et des points de contrôle. De plus la présence d’un étalon interne (89Y) permet de
s’affranchir d’une dérive éventuelle de l’appareil.
L’ICP-MS utilisé est le 810/820-MS de chez Varian Incorporation.
Les limites de détection et de quantification de ces deux méthodes d’analyse sont données en
annexe 4.
Les principales impuretés présentes dans les trois oxydes utilisés sont indiquées dans le
tableau IV-3.
Des impuretés dues aux méthodes de fabrications des poudres sont également présentes mais
en plus faibles quantités.
73
IV.2.4.2. Lavage des poudres
Le pourcentage massique de certains éléments étant assez élevé, le lavage des poudres est
alors nécessaire avant toute expérience.
Celui-ci consiste en la mise en contact du solide avec des solutions acides et/ou basiques afin
de désorber les éléments en surface de la poudre puis à des rinçages à l’eau ultra pure. Les
protocoles optimisés (Barale-2006) sont donnés en annexe 3. Ceux-ci sont jugés suffisants
lorsque le PZC des oxydes reste constant d’un lavage à l’autre. Les variations de point de
charge nulle des solides étudiés sont données dans le tableau IV - 4.
IV.3. Méthodes
Le point de charge nulle est mesuré grâce à un titrage de masse. Cette méthode repose sur le
fait qu’une augmentation de la masse de solide dans une suspension, contenant au départ un
électrolyte à pH fixé, provoque une évolution de son pH avec la variation de la concentration
en poudre. Cette évolution, lorsque la concentration augmente, tend vers une valeur stable
correspondant au PZC (Tewari-1976), et ce pour les trois milieux de départ (acide, basique et
neutre). Rappelons que le PZC correspond au pH où la charge de surface est nulle, c’est à dire
que le nombre de groupements fonctionnels positifs est égal au nombre de groupements
fonctionnels négatifs. Les courbes de titrage de masse sont données en annexe 5.
IV.3.1.2. Zétamétrie
La mesure de la mobilité électrophorétique des particules, c’est à dire leur vitesse lorsqu’elles
sont soumises à un champ électrique, permet de déterminer leur potentiel zêta. L’équation de
base, pour une particule sphérique, est de la forme :
ν = µE E Equation IV - 5
74
où ν est la vitesse des particules, µE leur mobilité électrophorétique et E le champ électrique
appliqué.
En présence de particules poreuses ayant une double couche fine, il a été démontré que le
champ appliqué est modifié. Overbeek (Overbeek-1952) a démontré que l’équation IV – 5 est
toujours valable mais que le champ E n’est alors que local. La relation entre la vitesse de la
particule et le champ électrique appliqué devient donc :
ε 0ε rζ
ν = E Equation IV - 6
η
Soit :
ε 0ε rζ
µE = Equation IV - 7
η
Dans le cas d’une double couche épaisse (Hückel-1924), le champ électrique est peu affecté
par les particules. L’équation est donc :
2ε 0ε rζ
µE = Equation IV - 8
3η
La contradiction entre les deux formules a été résolue par Henry (Henry-1931) qui a pris en
compte la façon dont les particules influent sur les lignes de champ à leur voisinage (figure IV
– 6).
Afin d’étudier l’effet de l’épaisseur de la double couche sur la mobilité, Henry a calculé µE
pour des particules sphériques ayant une épaisseur de la double couche arbitraire et en
supposant que la densité de charge n’est pas affectée par le champ appliqué. Cette hypothèse
est valide à condition que le potentiel ζ soit suffisamment faible. Pour déterminer la mobilité
électrophorétique d’une particule avec un faible potentiel ζ, il est alors nécessaire de résoudre
l’équation de Stokes (équation IV – 9) avec comme contraintes le fait que la vitesse tend vers
zéro à forte distance de la particule et que la force appliquée sur celle-ci est nulle.
η∇ 2ν − ∇p = ρe ∇Ψ Equation IV - 9
75
où ρe est la densité volumique de charge, η la viscosité du fluide, Ψ le potentiel électrique
local, ν la vitesse du fluide et p la pression appliquée. Avec ∇ ν = 0 puisque le fluide est
incompressible.
La mobilité électrophorétique est donc obtenue avec la relation :
2ε 0ε rζ
µE = f1 (κa) Equation IV - 10
3η
1
f1 (κa) = 1 + Equation IV - 11
2[1 + (δ /κa)]
3
avec :
2,5
δ= Equation IV - 12
[1 + [Link](−κa)]
Lorsque l’épaisseur de la double couche est plus faible que le rayon de la particule (κa >> 1),
les lignes de champ électriques sont considérées parallèles à la surface de la particule (figure
IV - 6). f1 (κa) vaut alors 2/3 et la mobilité électrophorétique est exprimée par l’équation de
Schmoluchowski (Schmoluchowski-1921) :
ε 0ε rζ
µE = Equation IV - 13
η
Lorsque κa << 1, c’est à dire lorsque l’épaisseur de la double couche est supérieure au rayon
de la particule, les lignes de champ électrique sont très peu perturbées (figure IV - 6).
L’équation IV – 14 est donc obtenue.
2 ε 0ε rζ
µE = Equation IV - 14
3 η
76
a b
Figure IV - 6 : Effet de la taille de la particule sur les lignes de champ électrique. (a) :
κa << 1, (b) : κa >> 1. La ligne en pointillée est à une distance 1/κ de la surface de la
particule (Hunter-2001)
La mesure de µE permet donc de calculer le potentiel ζ qui ne peut être mesuré directement.
IV.[Link]. Principe
Dans un premier temps, un titrage sans solide, appelé blanc, est effectué. Il permet de mesurer
la perte en protons due à l’électrolyte et à l’appareillage. On réalise ensuite un titrage en
présence de la poudre étudiée en conservant les mêmes conditions expérimentales. La
soustraction de la courbe du blanc à celle du titrage de la poudre permet de déterminer la
quantité de protons réellement adsorbés ou désorbés par la surface du solide, et ainsi de
remonter à la charge de surface. Celle-ci est calculée, en utilisant le modèle 2-pK (chapitre
III) à partir de l’équation d’électroneutralité IV – 15.
où :
- [≡ XOH 2+ ] et [≡ XO − ] sont respectivement les concentrations en sites positifs et négatifs à
la surface
- [ H + ] et [OH − ] sont les concentrations mesurées à l’aide de l’électrode pH
- [ B + ] et [ A− ] sont respectivement les concentrations en base et en acide ajoutés lors du
titrage
77
La courbe représentant la charge de surface en fonction du pH permet d’obtenir la densité de
sites et le PZC (voir chapitre III). La modélisation de celle-ci permet la détermination des
constantes acide-base du solide dans le modèle 2-pK.
Les titrages acide-base des trois poudres étudiées ont été réalisés à 25°C dans deux études
précédentes (Barale-2006, Mansour-2007).
Les concentrations en solide dans le réacteur sont de 10 g/L pour la magnétite et de 20 g/L
pour les ferrites et celle en KCl de 10-2 M. La température est stabilisée à 25°C ± 0,1°C à
l’aide d’un bain thermostaté. Après 24 h d’hydratation, le pH est fixé à 4 par ajout de HCl
0,1 M pour les ferrites et le dosage se fait par ajouts successifs de 10 µL de KOH 0,1 M à
l’aide d’une chaîne de titrage automatique (Dosimat 765 de chez Metrohm et logiciel
développé au LECIME). Pour la magnétite, deux titrages sont effectués en sens inverse. Le
pH après hydratation est de 5,9. Le titrage est alors réalisé en deux parties afin d’éviter la
dissolution de la poudre à pH très acide, respectivement avec KOH 0,1 M et HCl 0,1 M.
78
Les paramètres ajustables dans le modèle 2-pK/BSM sont les deux constantes acide-base pK+
et pK-, la densité de site, la capacitance et les deux constantes de sorption des ions du sel de
fond pour la surface, pKK et pKCl.
Ceci constitue un nombre trop important de paramètres pour qu’ils puissent être ajustés à
partir de la courbe de titrage.
Ainsi, la densité de site de surface a été choisie à 10 [Link]-2 pour chacun des solides. A partir
d’images obtenues au MEB (Barale-2008), il n’a pas été possible de déterminer la nature des
faces exposées à la solution, une valeur du même ordre que celle de la densité moyenne des
atomes d’oxygène pouvant être adsorbés à la surface des oxydes a donc été adoptée. Par
exemple, la densité de sites de la face (113) de l’hématite est de 4,1 [Link]-2 pour les
groupements hydroxyles monocoordinés et doublement coordinés (Barron-1996).
où :
K−
KK = Equation IV - 20
K'K
K+
K Cl = Equation IV - 21
K ' Cl
D’après Christl (Christl-1999), des constantes égales pour K+ et Cl- et pour les trois solides
ont été prises en compte :
pK’K = pK’Cl = 2,5
79
Les deux constantes d’acidité sont reliées par :
log K + + log K −
= PZC Equation IV - 22
2
Finalement, les paramètres ajustés sont la différence entre pK, ∆pK et la capacitance.
L’optimisation de ces paramètres s’effectue de façon à minimiser le carré des écarts entre les
valeurs expérimentales et les valeurs calculées à travers le WSOS/DF (weighted sum of
squares divided by degree of freedom).
IV.3.2. Sorption
Les expériences de sorption ont été réalisées à 25°C et 50°C. La méthode utilisée est celle des
batchs. Une quantité fixe de poudre (20 mg) est introduite dans un flacon en PEHD
(polyéthylène haute densité) puis une solution contenant le sel de fond à la concentration
souhaitée est ajoutée. Comme il a été vu précédemment, un temps de contact minimum est
nécessaire entre le solide et la solution afin que l’hydratation de la poudre soit optimum.
L’échantillon est donc mis sous agitation durant le temps nécessaire (24 h) avant que
l’élément d’intérêt ne soit ajouté à la solution. Pour les différentes séries de mesure, la
concentration en solide choisie est de 0,4 g/L2 et celle en sel de fond (KCl3) de 10-4 M. Le
2
Une concentration de 400 ppm a été choisie pour deux raisons. Premièrement, la zétamétrie ne
permet de travailler que pour une concentration en solide comprise entre 100 et 1000 ppm.
Deuxièmement, une quantité suffisante en poudre doit être utilisée pour que les différents phénomènes
80
cobalt ou le nickel sont ensuite ajoutés à la solution pour obtenir la concentration souhaitée
puis celle-ci est placée 5 min sous ultrasons et laissée sous agitation le temps nécessaire.
Celui-ci est déterminé préalablement par une étude de cinétique. Le pH final des échantillons
est mesuré puis ils sont filtrés avec un filtre en acétate de cellulose (0,2 µm) et acidifiés (1%
en acide nitrique) pour stabiliser les solutions de cations avant l’analyse.
Pour les échantillons réalisés à 50°C une étuve munie d’une table d’agitation a été utilisée.
Les analyses des solutions de sorption ont été réalisées par ICP-AES et ICP-MS. Les
principes de fonctionnement de ces deux appareils ainsi que les limites de détection et de
quantification déterminées pour les éléments étudiés sont donnés en annexe 4.
Les paramètres obtenus grâce à la modélisation des courbes de titrage acide-base (tableau IV-
5) ont été utilisés pour la modélisation des courbes de sorption par le logiciel ECOSAT
(Keizer-2009). Le logiciel FITEQL n’a pu être utilisé pour la modélisation des courbes de
sorption car il ne prend pas en compte le phénomène de précipitation.
Le modèle choisi pour réaliser la modélisation est également le 2-pK/BSM. Différentes
réactions de complexation et de précipitation de surface sont implémentées dans le logiciel et
testées pour chacun des systèmes.
soient observés, mais elle ne doit pas être trop forte pour qu’en cas de présence de plusieurs
mécanismes l’un d’eux soit pas prépondérant par rapport à l’autre, comme par exemple la précipitation
vis-à-vis de la complexation.
3
KCl étant considéré comme indifférent vis-à-vis des oxydes, c’est donc ce sel de fond qui a été choisi
pour réaliser les différentes expériences (Hunter-2001).
81
M 2+ + OH − ⇄ MO + H + Equation IV - 27
Les constantes de complexation et les produits de solubilité de ces espèces sont ajustés
manuellement. Le choix de la courbe la plus représentative du système se fait par le calcul des
écarts entre les points expérimentaux et ceux résultants de la modélisation.
Le choix des complexes et des réactions de précipitation sera discuté dans le chapitre V. Les
constantes des réactions seront données dans le même chapitre.
Il faut noter que la force ionique varie avec la concentration en chlorure de cobalt ou de nickel
en solution. Cette variation est prise en compte lors de la modélisation.
IV.3.3. XPS
L’analyse de la surface d’un matériau par spectroscopie de photoélectrons induits par rayons
X permet de connaître :
- les éléments présents en surface (à l’exception de H et He) ;
- les états chimiques de ces éléments ;
- la proportion de ces états chimiques pour chaque élément.
82
où :
- hν est l’énergie des photons incidents
- E c et E l (1s) sont respectivement les énergies cinétique et de liaison de l’orbitale 1s du
photoélectron éjecté
- Φspectro est le travail de sortie du spectromètre
La valeur de Φspectro doit être déterminée en utilisant des solides de référence. Pour des
échantillons isolants, un terme supplémentaire relié à l’effet de charge doit être rajouté.
L’appareil utilisé est un spectromètre XPS Escalab 250 de chez Thermo Fisher Scientific
(figure IV - 8). La source de photons utilisée est une raie monochromatisée Kα émise par des
cibles d’aluminium (énergie de 1486,6 eV). La tension utilisée est de 15 kV, pour une
puissance de 150 Watts. L’angle d’émergence est de 90°. Les analyses peuvent être couplées,
le cas échéant, à des séquences d’abrasion ionique par bombardement d’ions argon de
manière à obtenir un profil des espèces en surface. Cette méthode est intéressante dans le cas
où l’épaisseur de la couche à analyser est supérieure à la profondeur de résolution de l’XPS
(environ 10 nm).
Analyseur hémisphérique
Chambre d’analyse
Sas d’introduction
Une fois les spectres XPS obtenus, une décomposition est nécessaire. Les énergies des pics
observés sont comparées avec les valeurs disponibles dans des tables de références. Afin
d’avoir une meilleure comparaison, des échantillons témoins, de divers oxydes pouvant être
obtenus à la surface des matériaux, ont été analysés. Les énergies des pics utilisées par la suite
pour la décomposition des spectres, sont données dans le tableau IV – 6.
83
Tableau IV - 6 : Energies caractéristiques des espèces étudiées - (a) : Grosvenor-2006,
(b) : Galtayries-1999, (c) : Kim-1998, (d) : McIntyre-1975
Espèces Ni (OH ) 2 NiO Co(OH ) 2 CoO
Raie Ni 2p3 Ni 2p3 Co 2p3 Co 2p3
Energie
855,8 eV 854 eV 781,2 eV 781 eV
Ech. Témoins
Energie 855,5 eV (a) 853,9 eV (a), (b) 780,7 eV (c) 779,7 eV (c)
Bibliographie 855,6 eV (d) 854 eV (d) 781 eV (d) 780 eV (d)
Les spectres de NiO et Ni (OH ) 2 indiquent la présence de multiplet (Ni 2p3) pour ces deux
espèces. Dans le tableau IV – 6, les énergies d’un seul pic de ce multiplet sont comparées. Les
énergies de l’oxyde et de l’hydroxyde d’un même élément, étant distinctes d’environ 1 eV,
leur différenciation est donc possible.
La connaissance des propriétés de surface des parois du circuit primaire est nécessaire si l’on
veut pouvoir les comparer aux charges de surface des produits de corrosion afin de mieux
comprendre leur comportement dans le fluide. Dans ce but, leur étude par la méthode du
potentiel d’écoulement a été réalisée.
Ces expériences ont été réalisées au KIT de Karlsruhe à l’aide d’un Zétamètre SurPASS de
chez Anton Paar (figure IV - 9).
84
Deux échantillons identiques sont placés l’un en face de l’autre dans l’appareil. Lorsqu’une
pression est appliquée sur le liquide placé entre ces deux plaques, les charges de la partie
mobile de la double couche sont mises en mouvement par le fluide (figure IV – 10). Il se
forme alors un courant d’écoulement Is. La différence de potentiel électrique résultant de ce
déséquilibre de charges génère un courant de conduction, Ic, opposé au flux de la solution et
donc à Is. Un état stationnaire est rapidement atteint, se traduisant par l’égalité de ces deux
courants :
I = Ic + Is = 0 Equation IV - 29
∆U η
ζ= × × κB Equation IV - 30
∆P εε 0
85
Figure IV - 11 : Schéma de la mise en place des échantillons dans le zétamètre
L’Inconel 690, le Zircaloy 4 et l’acier inoxydable 304L ont été étudiés. Ceux-ci ont été polis
(granulométrie 600 µm). Pour chaque matériau deux échantillons ont été vieillis 80 jours en
autoclave avec un milieu primaire (1000 ppm de bore, 2 ppm de lithine, 30 [Link]-1
d’hydrogène, 360°C4). L’électrolyte utilisé est NaNO3 (10-2 et 10-3 M). La distance entre les
deux plaques d’échantillon est de 100 µm.
4
La température utilisée est plus haute que celle du circuit primaire afin d’accélérer le phénomène
d’oxydation.
86
V. Etude des mécanismes de sorption dans le circuit
primaire des centrales REP
Préalablement aux mesures proprement dites de la quantité sorbée en fonction du pH, il est
nécessaire de connaître la cinétique de la réaction afin de déterminer le temps de contact
nécessaire pour être à l’équilibre. Pour cela, des essais en batch ont été réalisés à des temps
variables, tous les paramètres étant égaux par ailleurs. La concentration initiale en nickel a été
choisie à 10-6 M et le pH à 6 dans le but de limiter le risque de précipitation. La figure V – 1
représente l’évolution de la quantité de nickel sorbée sur le ferrite de cobalt, exprimée en
atomes de nickel par nm2 de surface, en fonction du temps de contact. Celui-ci varie entre 1 h
et 24 h.
0.05
0.04
([Link] )
-2
0.03
sorbé
0.02
Ni
0.01
0
0 5 10 15 20 25
Un palier de sorption est atteint à partir de 7 heures. Pour la suite des expériences, la solution
de chlorure de nickel a donc été laissée en contact avec le ferrite de cobalt durant 8 heures
avant d’être filtrée.
Le temps de contact optimum ayant été déterminé, les mesures de sorption en fonction du pH
ont été réalisées pour quatre concentrations initiales en chlorure de nickel, à savoir 10-6 M,
10-5 M, 5.10-5 M et 10-4 M. Les figures V – 2 et V – 3 regroupent les résultats.
87
Les courbes observées ont une allure caractéristique de la sorption de cations sur les oxydes
(Christl-1999, Ridley-1999) : il y a augmentation de la quantité sorbée en fonction du pH –
avec une variation en forme de sigmoïde. Les ions nickel sont donc bien sorbés sur le ferrite
de cobalt dans ces conditions.
Un front de sorption s’étalant de pH 5,5 à pH 8 est observé pour les quatre concentrations. La
concentration sorbée exprimée en atomes par unité de surface (figure V – 2) atteint 8,5 [Link]-2
pour la concentration la plus élevée en solution. Pour les deux plus faibles concentrations
(10-6 M et 10-5 M) les concentrations sorbées sont aussi exprimées en moles par litre de
solution (unité généralement utilisée pour les modélisations). Ces courbes (figure V – 3)
montrent qu’à ces concentrations, le nickel est sorbé pratiquement à 100 % à partir de pH 7,5
et pH 9, respectivement.
A cause de la faible solubilité du nickel au-dessus de pH 8, comme le prévoient les calculs de
spéciation effectués avec le logiciel JChess (III.7), la complexation de surface peut ne pas être
le seul phénomène expliquant la sorption et il faut envisager d’autres mécanismes tels que la
précipitation.
10
-6
10 M
10-5 M
8
5.10-5 M
-4
10 M
([Link] )
-2
6
sorbé
4
Ni
0
4 5 6 7 8 9 10
pH
88
-6 -5
1 x 10 1.2 x 10
-5
1 x 10
8 x 10-7
-6
8 x 10
(M)
-7
6 x 10
(M)
sorbé
-6
sorbé
6 x 10
Ni
Ni
-7
4 x 10
-6
4 x 10
2 x 10-7 -6
2 x 10
0 0
4 5 6 7 8 9 10 4 5 6 7 8 9 10
pH pH
Dans la suite du texte et sur les figures, la concentration en cation sur la phase solide (obtenue
par différence entre la concentration initiale et la concentration finale en solution) rapportée à
la surface du solide, quelque soit le mécanisme de cette sorption, sera notée « concentration
sorbée ». Pour comprendre ces mécanismes, il est nécessaire de caractériser le ferrite de
cobalt après sorption par des techniques complémentaires (XPS, zétamétrie) en comparant les
résultats à ceux obtenus avant sorption.
V.1.3.1. XPS
Dans le but de déterminer l’espèce formée à la surface de celui-ci, des analyses XPS ont été
réalisées avant et après sorption dans les conditions suivantes :
- 25°C,
- [ KCl ] = 10-4 M,
- [ NiCl2 ] = 10-3 M, 10-5 M, 10-6 M et 0 M,
- [ CoFe2O4 ] = 0,4 g.L-1,
- Temps de contact : 8 h,
- pH 10,5.
Les échantillons ont été filtrés (0,2 µm) puis séchés avant d’être introduits dans l’appareil. La
figure V – 4 regroupe les spectres XPS du nickel pour les quatre essais.
89
-3
10 M
-5
10 M
-6
10 M
sans NiCl
2
Figure V - 4 : Spectres XPS du nickel avant et après sorption ([NiCl2] = 10-3 M, 10-5 M,
10-6 M et 0 M) sur le ferrite de cobalt à pH 10,55 (B.E. représente l’énergie de liaison des
électrons)6
Les spectres des échantillons sans nickel ou avec une faible concentration (10-6 M)
n’indiquent pas de bande correspondant à cet élément. Pour les concentrations supérieures ou
égales à 10-5 M, un multiplet, correspondant à Ni2p3 de l’hydroxyde de nickel (Grosvenor-
2006), apparaît sur les spectres entre 856 eV et 863 eV. Son intensité augmente avec la
concentration en nickel. Cette analyse montre l’apparition de nickel à la surface du ferrite de
cobalt pour des concentrations en chlorure de nickel supérieures ou égales à 10-5 M. Celui-ci
peut se trouver sous forme d’hydroxyde de nickel.
Les spectres du fer et du cobalt ont également été réalisés sur ces échantillons. Sur ceux du
fer, le pic à 711,1 eV correspond à la raie 2p3 de l’élément (Grosvenor-2004). Celui du cobalt
indique une structure en trois pics avec des énergies comprises entre 775 eV et 792 eV. Le pic
à 780 eV correspond à Co2p3 tandis que ceux à 783 eV et 787 eV sont des pics Auger du fer
(2p3) qui interfèrent avec le spectre du cobalt (Tewari-1975). La sorption du nickel ne
modifie pas ces spectres.
Les observations XPS indiquent donc la possibilité d’un processus de précipitation à partir
d’une concentration en nickel en solution de 10-5 M. En deçà de cette concentration, le nickel
sorbé n’est pas détecté. Les résultats sont en accord avec les prédictions JChess (paragraphe
5
Le pH a été ajusté à 10,5 lors des expériences de sorption afin d’être placé dans des conditions
favorables à la précipitation.
6
L’ordonnée des spectres XPS représente des coups/s, cependant aucune valeur n’est indiquée car les
courbes ont été superposées à différents niveaux pour une meilleure lisibilité.
90
III.7.1) quant à la possible précipitation de l’hydroxyde de nickel Ni (OH ) 2 et non de l’oxyde
NiO .
De façon à ce que la comparaison avec les courbes de sorption soit possible, les mêmes
conditions expérimentales ont été choisies :
- 25°C,
- Temps d’hydratation 24 h,
- [ KCl ] = 10-4 M,
- [ NiCl2 ] = 10-3 M, 5.10-5 M, 10-5 M, 10-6 M et 0 M,
- [ CoFe2O4 ] = 0,4 g.L-1.
40
Potentiel zéta (mV)
20
sans NiCl
2
-20 10 M
-6
-5
10 M
-5
5.10 M
-40 -4
10 M
5 6 7 8 9 10 11 12
pH
Figure V - 5 : Variation du potentiel zêta du ferrite de cobalt en fonction du pH pour
différentes concentrations initiales en chlorure de nickel (25°C, [KCl] = 10-4 M,
[NiCl2] = 10-4 M, 5.10-5 M, 10-5 M, 10-6 M et 0 M, [CoFe2O4] = 0,4 g.L-1)
91
En l’absence ou pour une faible concentration (10-6 M) en nickel, le point isoélectrique obtenu
est de 6,9 ± 0,1, ce qui est concordant avec les mesures antérieures pour le ferrite de cobalt
(Barale-2006). A une concentration de 10-5 M, les courbes montrent une légère augmentation
du potentiel zêta entre pH 8 et pH 10, sans modification significative du PIE. Pour les deux
concentrations les plus hautes (5.10-5 M et 10-4 M), on observe un décalage positif entre pH 5
et pH 7, puis une forte augmentation du potentiel zêta avec un maximum aux alentours de pH
9. Pour ces deux concentrations le PIE est de 11.
Les analyses XPS et la zétamétrie ont permis de montrer la présence d’un précipité
d’hydroxyde de nickel Ni (OH ) 2 en surface du ferrite de cobalt, pour les expériences de
sorption dont la concentration initiale en chlorure de nickel est supérieure ou égale à 10-5 M et
pour des pH supérieurs à 7. A des pH et des concentrations plus faibles, un mécanisme de
sorption faisant appel à la formation de complexes de surface peut être envisagé. Les
méthodes précédentes n’ont pas permis d’identifier la nature de ces complexes de surface. La
zétamétrie a simplement montré que, à pH inférieur à 7, la sorption conduit à une charge de
surface plus positive.
Les analyses de surface réalisées par XPS et l’évolution du point isoélectrique imposent la
présence d’hydroxyde de nickel : Ni (OH ) 2 , pour les conditions de pH et de concentration où
la précipitation se produit.
92
Par ailleurs, les complexes de surface suivants peuvent être envisagés :
avec une préférence pour les complexes positifs ou neutres, afin de respecter l’évolution du
potentiel zêta à pH inférieur à 7.
La modélisation a été réalisée comme décrit dans le paragraphe IV.3.2. Les paramètres acide-
base, sont ceux pris en compte par le modèle 2-pK/BSM (Basic Stern Model) (tableau IV –
5), d’après les titrages réalisés par Barale (Barale-2006).
Pour chaque courbe de sorption, l’ajustement par rapport aux points expérimentaux a été
réalisé de façon à obtenir la plus faible erreur possible, en faisant varier les constantes de
formation des espèces envisagées.
Le produit de solubilité de Ni (OH ) 2 a été considéré comme une variable ajustable bien que
des valeurs soient disponibles dans la littérature. En effet, celle-ci varie fortement avec les
conditions de synthèse du précipité, ou son vieillissement (Marthur-1998, Osinska-2004,
Plyasunova-1998b).
Parmi les quatre complexes de surface proposés ci-dessus et testés pour la modélisation des
courbes de sorption, seuls deux, ≡ XONi + et ≡ XONiOH , ont permis un ajustement correct
par rapport aux points expérimentaux à 10-6 M et un seul, ≡ XONi + , à 10-4 M (figure V – 6).
Pour réaliser cette modélisation, les paramètres choisis sont ceux donnés dans le tableau IV –
5. Le modèle pris en compte est le 2-pK/BSM.
-4
1 x 10
expérimental
+
XONi
-5 XONiOH
8 x 10
(XO) Ni
2
-
XONi(OH)
(M)
2
-5
6 x 10
sorbé
Ni
-5
4 x 10
-5
2 x 10
0
3 4 5 6 7 8 9 10
pH
93
Afin de faire un choix entre ≡ XONi + et ≡ XONiOH , pour une concentration initiale en
chlorure de nickel de 10-6 M, la charge de surface induite par chacun d’eux a été comparée à
celle observée par zétamétrie pour un pH de 7. Les équations prises en compte pour ce calcul
sont les suivantes :
-8 +
2 x 10 XONi
XONiOH
4
2
-8
-2 x 10
-8
-4 x 10
-8
-6 x 10
-8
-8 x 10
-1 x 10-7
6.4 6.6 6.8 7 7.2 7.4 7.6
pH
La charge de surface induite par le complexe ≡ XONi + est plus positive que celle imposée par
≡ XONiOH et correspond mieux à l’évolution observée par zétamétrie. C’est donc le
complexe monodentate positif qui est choisi pour modéliser les données de sorption.
Les courbes obtenues pour les quatre concentrations indiquent un front de sorption s’étalant
de pH 5,5 à pH 8. Pour les deux concentrations les plus faibles (10-5 M et 10-6 M), les courbes
94
se modélisent avec le seul phénomène de complexation de surface, ce qui est en accord avec
le produit de solubilité de l’hydroxyde. Pour les deux concentrations supérieures (5.10-5 M et
10-4 M) l’ajustement des courbes est le meilleur lorsque la précipitation et la complexation
sont prises en compte simultanément. Il est également visible que la quantité de nickel
précipité augmente avec la concentration initiale en ions. Cependant la courbe réalisée à
10-5 M indique une moins bonne corrélation entre la modélisation et les points expérimentaux
que celles aux trois autres concentrations. Dans les quatre cas, la modélisation est effectuée
avec le complexe de surface ≡ XONi + et le précipité Ni(OH ) 2 .
95
-6 -5
1 x 10 1.2 x 10
a b
8 x 10-7
1 x 10-5
-6
8 x 10
6 x 10-7
(M)
(M)
sorbé
sorbé
-6
6 x 10
Ni
Ni
4 x 10-7
-6
4 x 10
+
2 x 10-7 modélisation (>XONi )
+
expérimental modélisation (>XONi )
2 x 10-6 expérimental
0
4 5 6 7 8 9 10 0
4 5 6 7 8 9 10
pH
pH
-5 -4
5 x 10 1 x 10
c Ni(OH) d
2 précipité Ni(OH) 2 précipité
+
>XONi
4 x 10
-5
8 x 10-5 >XONi+
modélisation
modélisation
expérimental
expérimental
-5
(M)
3 x 10 6 x 10-5
(M)
sorbé
sorbé
Ni
-5
Ni
2 x 10 4 x 10-5
-5
1 x 10 -5
2 x 10
0
0
4 5 6 7 8 9 10
4 5 6 7 8 9 10
pH pH
Figure V - 8 : Evolution de la quantité de nickel sorbée en moles par litre de solution (a,
b, c, d : 10-6 M, 10-5 M, 5.10-5 M, 10-4 M respectivement) sur le ferrite de cobalt. (valeurs
expérimentales ( ) et modélisation de la quantité sorbée totale () par complexation
(- -) et par précipitation (- -)) (25°C, [KCl] = 10-4 M, [CoFe2O4] = 0,4 g.L-1, 8 h)
Pour les pH les plus élevés, le complexe ≡ XONi + ne devrait plus être présent dans la solution
du fait de la grande concentration en ions OH − . Il est probable que ces deux espèces
réagissent l’une avec l’autre (équation V – 3) pour former ≡ XONiOH .
96
Tableau V - 1 : Logarithme des constantes de complexation et de solubilité, déterminées
à l’aide du logiciel ECOSAT (2-pK/BSM) pour le système nickel / ferrite de cobalt
Espèce
≡ XONi + Ni (OH ) 2
Concentration initiale
10-6 M -2,3 /
-5
10 M -2,3 /
-5
5.10 M -2,3 -18
-4
10 M -2,2 -18
Dans la littérature, des valeurs pour le logKNi(OH)2 sont comprises dans un intervalle allant de
-10,5 à -19 (Baes-1986, Baumgarten-1997, Chen-2009, Plyasunova-1998b, Zhu-2002). La
constante obtenue par modélisation des courbes de sorption est cohérente avec les données
bibliographiques7.
V.1.5. Conclusion
Les diverses expériences et analyses réalisées sur le système Ni 2 + / CoFe2O4 indiquent que la
complexation de surface est le seul phénomène entrant en jeu pour des concentrations
inférieures à 10-5 M. Ce phénomène est accompagné de la précipitation de Ni(OH ) 2 à partir
de pH 7 lorsque la teneur initiale en nickel en solution est égale ou supérieure à 10-5 M. Le
complexe ≡ XONi + se forme à partir de pH 5,5 puis réagit probablement avec les OH − de la
solution pour former ≡ XONiOH pour des pH plus élevés. Les constantes de ces espèces
déterminées à l’aide du logiciel ECOSAT sont logK≡XONi+ = -2,25 ± 0,05 et logKNi(OH)2 = -18.
7
Rappelons que ce grand intervalle pour la valeur de logKNi(OH)2 dépend du mode et du temps de
formation de l’espèce ainsi que du temps depuis lequel elle a été formée (Marthur-1998, Osinska-
2004, Plyasunova-1998b).
97
V.2. Système Co2+ - ferrite de nickel
Les expérimentations sur le système Co 2+ / NiFe2O4 , ont été réalisées dans les mêmes
conditions expérimentales que pour le système Ni 2+ / CoFe2O4 , de façon à ce qu’une
comparaison puisse être effectuée.
A 25°C, le temps de contact minimum pour la sorption du cobalt sur le ferrite de nickel est de
5 heures (figure V – 9). Tout comme pour les expériences de sorption du nickel sur le ferrite
de cobalt, 8 heures de contact seront retenues pour les expériences de sorption.
0.02
0.015
([Link] )
-2
0.01
sorbé
Co
0.005
0
0 5 10 15 20 25
La sorption du cobalt à différentes concentrations initiales est représentée sur la figure V - 10.
Pour ce système, comme pour le précédent, la présence simultanée de deux phénomènes
(complexation de surface et précipitation) peut être supposée, du fait que le palier est atteint
pour des pH élevés, de la valeur de la constante de solubilité de l’hydroxyde de cobalt, mais
également des prédictions JChess.
Le front de sorption est visible à partir de pH 6 et un palier de sorption est atteint à partir de
pH 9. En comparaison au système précédent, la sorption débute pour des pH légèrement plus
élevés. Les concentrations sorbées en [Link]-2 sont plus faibles pour le cobalt sur le ferrite de
nickel que pour le nickel sur le ferrite de cobalt. Ceci peut être expliqué par la différence entre
les surfaces spécifiques des deux solides (60 et 17 m2.g-1, respectivement). Le rapport des
concentrations sorbées pour les deux éléments (pour des concentrations en solution et des pH
comparables) est approximativement égal au rapport inverse des surfaces spécifiques (environ
0,3).
98
3
-6
10 M
2.5
10-5 M
-5
5.10 M
2
([Link] )
-4
10 M
-2
1.5
sorbé
Co
0.5
0
4 5 6 7 8 9 10 11 12
pH
Figure V - 10 : Variation de la concentration en cobalt sorbé sur le ferrite de nickel en
fonction du pH (25°C, [KCl] = 10-4 M, [CoCl2] = 10-4 M, 5.10-5 M, 10-5 M et 10-6 M,
[NiFe2O4] = 0,4 g.L-1, 8 h)
-5
1 x 10
-6
1 x 10
8 x 10-6
-7
8 x 10
-6
(M)
6 x 10
(M)
-7
6 x 10
sorbé
sorbé
Co
Co
-6
4 x 10
-7
4 x 10
2 x 10-6
2 x 10-7
0
0
5 6 7 8 9 10 11
5 6 7 8 9 10 11
pH
pH
99
V.2.3. Détermination de la nature des processus de sorption
V.2.3.1. XPS
Afin de vérifier l’hypothèse faite quant à la présence d’un précipité à la surface de l’oxyde
aux fortes concentrations en chlorure de cobalt, l’analyse de la poudre avant et après sorption
(conditions identiques aux expériences de sorption décrites précédemment) du chlorure de
cobalt à 10-3 M, a été réalisée par XPS. Les spectres sont représentés sur la figure V – 12.
sans CoCl
2
-4
10 M
B. E. (eV)
La bande comprise entre 780 et 790 eV sur le spectre avant sorption correspond à un pic
Auger du nickel. Afin de décomposer le spectre après sorption du cobalt et de pouvoir
comparer les deux résultats, les paramètres de cette bande (énergies de début et de fin de la
bande) ont été maintenus constants. Sur le spectre après sorption du chlorure de cobalt à
10-3 M, un multiplet correspondant au cobalt 2p3 peut être observé entre 778 eV et 801 eV.
La décomposition du spectre indique la présence possible d’hydroxyde de cobalt Co(OH ) 2 à
781,5 eV et de Co2p1 à 797,5 eV. Cette observation est en accord avec les prédictions JChess
et la constante de solubilité de l’hydroxyde de cobalt.
Afin de supprimer les pics Augers du nickel, un changement de source pourrait être envisagé
pour des expériences ultérieures.
L’influence du cobalt sur le point isoélectrique du ferrite de nickel est représentée sur la
figure V – 13. Les mesures ont été réalisées pour les conditions expérimentales suivantes :
100
- 25°C,
- Temps d’hydratation 24 h,
- [ KCl ] = 10-4 M,
- [ CoCl2 ] = 10-3 M, 10-5 M, 10-6 M et 0 M,
- [ NiFe2O4 ] = 0,4 g.L-1.
40
30
20
Potentiel zéta (mV)
10
0 sans CoCl
2
-6
10 M
-10 -5
10 M
-5
5.10 M
-20 -4
10 M
-30
5 6 7 8 9 10 11 12
pH
Figure V - 13 : Variation du potentiel zêta du ferrite de nickel en fonction du pH pour
différentes concentrations initiales en chlorure de cobalt (25°C, [KCl] = 10-4 M,
[CoCl2] = 10-4 M, 5.10-5 M, 10-5 M, 10-6 M et 0 M, [NiFe2O4] = 0,4 g.L-1)
Le PIE de la poudre en l’absence de chlorure de cobalt est de 9,1. Les mesures antérieures
(Barale-2006) ont donné une valeur de 9,4, ce qui est proche, sachant que le ferrite de nickel
présente une certaine instabilité. Un décalage des courbes et un déplacement du PIE à 9,6 est
observé pour les deux concentrations les plus faibles (10-6 M et 10-5 M). Ceci peut être
interprété par la formation d’un complexe de surface qui rend la charge de surface plus
positive. Ce décalage est plus important (PIE à 10,6) pour les concentrations supérieures, mais
les courbes ne présentent pas de maximum comme cela a été observé dans le cas du nickel
sorbé sur le ferrite de cobalt. Un point d’inflexion est simplement observé, ce qui ne permet
pas d’affirmer aussi nettement que dans le cas précédent, la formation d’une nouvelle phase.
Cependant, comme la valeur du PIE correspond à celle donnée dans la littérature pour
Co(OH ) 2 (Jurkiewicz-1986), à cause de la présence d’un point d’inflexion et compte tenu
des résultats de l’analyse XPS, l’hypothèse de la formation de cet hydroxyde, à partir d’une
concentration en cobalt de 5.10-5 M et d’un pH compris entre 7 et 9, a été faite.
101
Une expérience menée par Tewari et al. (Tewari-1975) sur l’influence des ions cobalt sur le
potentiel zêta d’un ferrite de nickel (figure V - 14) donne un résultat similaire à celui de cette
étude, avec cependant une influence beaucoup plus marquée par la présence d’un maximum.
Ce résultat confirme la formation d’hydroxyde de cobalt à la surface du ferrite de nickel.
Figure V - 14 : Influence des ions Co2+ sur la mobilité électrophorétique et le PIE d’un
ferrite de nickel (Tewari-1975)
Comme dans le cas du premier système étudié, compte tenu des résultats obtenus par
zétamétrie et par les analyses de surface réalisées par XPS, pour les concentrations et les pH
les plus élevés, la présence d’un hydroxyde, Co(OH ) 2 , comme précipité de surface, a été
imposée.
La constante de solubilité de Co(OH ) 2 a été considérée comme variable bien que sa valeur
soit disponible dans la littérature. En effet, celle-ci tout comme celle de l’hydroxyde de
nickel, varie fortement avec les conditions de synthèse du précipité, ou son vieillissement
(Marthur-1998, Plyasunova-1998a).
Dans ce cas également, parmi les quatre complexes proposés ci-dessus et testés pour la
modélisation des courbes de sorption, seuls deux, ≡ XOCo + et ≡ XOCoOH , ont permis un
ajustement correct par rapport aux points expérimentaux. Le calcul des charges de surface ne
permettant pas de privilégier l’un plutôt que l’autre, il a été choisi de prendre ≡ XOCo + pour
la modélisation comme dans le cas du système nickel/ferrite de cobalt.
102
-5
1 x 10
a 1 x 10
-6
modélisation (>XOCo )
+
b
expérimental modélisation (>XOCo +)
8 x 10-6 expérimental
-7
8 x 10
-6
(M)
6 x 10
(M)
-7
6 x 10
sorbé
sorbé
Co
Co
-6
4 x 10
-7
4 x 10
-7 2 x 10-6
2 x 10
0 0
4 5 6 7 8 9 10 4 5 6 7 8 9 10
pH pH
-5 -4
5 x 10 1.2 x 10
c Co(OH)
2 précipité
d Co(OH)
2 précipité
-4 +
-5 >XOCo+ 1 x 10 >XOCo
4 x 10
modélisation modélisation
expérimental expérimental
-5
8 x 10
(M)
-5
3 x 10
(M)
sorbé
sorbé
-5
6 x 10
Co
Co
-5
2 x 10
-5
4 x 10
1 x 10-5
-5
2 x 10
0 0
4 5 6 7 8 9 10 4 5 6 7 8 9 10
pH pH
Le front de sorption des courbes modélisées est compris entre pH 6 et pH 8,5, ce qui est
légèrement décalé vers les hauts pH par rapport au premier système. Pour les mesures
réalisées à faibles concentrations (10-6 M et 10-5 M), les courbes sont modélisées uniquement
avec le phénomène de complexation de surface. Pour les concentrations supérieures, la
précipitation de Co(OH ) 2 apparaît dès pH 8. Pour ce système également, pour une
concentration en chlorure de cobalt supérieure à 10-5 M, les deux phénomènes (complexation
et précipitation de surface) sont présents simultanément.
Pour les pH élevés, la présence de ≡ XOCoOH , formé par réaction du complexe ≡ XOCo +
avec les ions hydroxydes de la solution (équation V – 4), est probable dans ce cas aussi, mais
elle ne peut pas être modélisée.
103
Les constantes déterminées à l’aide du logiciel ECOSAT sont données dans le tableau V – 2.
Dans la littérature, logKCo(OH)2 est compris dans l’intervalle -9,8 à -19,25 (Baes-1986,
Baumgarten-1997, Feitknecht-1963, Plyasunova-1998a, Zhu-2002), avec dans ce cas
également, l’influence du mode de formation de l’hydroxyde.
La comparaison de l’allure des courbes de ce système avec celles du précédent indique que le
phénomène de sorption a lieu pour des pH un peu plus élevés (6 au lieu de 5,5) ainsi que
l’apparition du précipité (8 au lieu de 7). Cette observation est en accord avec les données
calculées avec JChess (figure III – 11), où le passage de la forme dissoute de l’hydroxyde à sa
forme précipité est observé à des concentrations plus faibles en solution pour le nickel que
pour le cobalt.
V.2.5. Conclusion
104
V.3. Systèmes Co2+ et Ni2+ - magnétite
La magnétite est, comme les ferrites, un produit de corrosion présent dans le circuit primaire
des centrales à eau sous pression. Le comportement des ions nickel et cobalt en présence de
cet oxyde a également été étudié.
Les cinétiques de sorption des ions cobalt et nickel sur la magnétite ont été réalisées dans les
mêmes conditions que sur les ferrites. Le temps de contact minimum afin d’atteindre le
plateau est de 8 heures environ (à pH 6) (figure V –16 et figure V - 17). Pour les deux cations,
les nombres de sites occupés correspondant au plateau de sorption sont identiques
(0,8 [Link]-2, à pH 6 pour une concentration initiale de 10-6 M pour les deux cations). Le
comportement des ions cobalt et nickel en présence de magnétite semble donc être identique.
Pour la suite des expériences, les poudres seront laissées en contact avec la solution sous
agitation pendant 8 heures, tout comme pour les ferrites.
0,8
([Link] )
-2
0,6
sorbé
0,4
Co
0,2
0
0 5 10 15 20 25 30
Temps de contact (h)
105
1
0,8
([Link] )
-2 0,6
sorbé
0,4
Ni
0,2
0
0 5 10 15 20 25
Temps de contact (h)
La sorption des ions cobalt et nickel sur la magnétite a été réalisée pour des conditions
identiques à celles choisies pour les ferrites :
- 25°C,
- [ KCl ] = 10-4 M,
- [ CoCl2 ] ou [ NiCl2 ] = 10-4 M, 10-5 M et 10-6 M,
- [ Fe3O4 ] = 0,4 g.L-1,
- Temps de contact : 8 h,
- Hydratation préalable : 24 h.
106
2
1 x 10
-6
10 M
-5
8 x 101 10 M
-4
10 M
([Link] )
-2
6 x 101
sorbé
4 x 101
Co
1
2 x 10
0
4 5 6 7 8 9 10 11 12
pH
Figure V - 18 : Variation de la concentration en cobalt sorbé sur la magnétite en
fonction du pH (25°C, [KCl] = 10-4 M, [CoCl2] = 10-4 M, 10-5 M et 10-6 M,
[Fe3O4] = 0,4 g.L-1, 8 h)
-6
1.2 x 10
1 x 10-6
8 x 10-7
(M)
sorbé
6 x 10-7
Co
-7
4 x 10
2 x 10-7
0
5 6 7 8 9 10 11
pH
107
2
1 x 10
-6
10 M
-5
8 x 101 10 M
-4
10 M
([Link] )
-2
1
6 x 10
sorbé
4 x 101
Ni
1
2 x 10
0
4 5 6 7 8 9 10 11 12
pH
Figure V - 20 : Variation de la concentration en nickel sorbé sur la magnétite en fonction
du pH (25°C, [KCl] = 10-4 M, [CoCl2] = 10-4 M, 10-5 M et 10-6 M, [Fe3O4] = 0,4 g.L-1, 8 h)
-6
1.2 x 10
1 x 10-6
8 x 10-7
(M)
sorbé
-7
6 x 10
Ni
-7
4 x 10
-7
2 x 10
0
5 6 7 8 9 10 11
pH
Pour les deux cations, la courbe à 10-4 M indique un maximum de sorption pour 90 [Link]-2. Il
semble donc qu’il y ait un processus de précipitation important car la densité de site de
108
surface de 10 [Link]-2 prise en compte pour la modélisation des propriétés acide-base de la
magnétite (tableau IV – 5), est largement dépassée. La courbe à 10-5 M atteint un palier pour
9 [Link]-2, ce qui est proche du nombre de sites de surface, mais, compte tenu des pH élevés et
des observations faites sur les systèmes nickel/ferrite de cobalt et cobalt/ferrite de nickel, à
cette concentration également la présence d’un précipité de surface est probable.
Le front de sorption pour le cobalt est compris entre pH 8 et pH 10,5, celui du nickel entre pH
8 et pH 9. Ce décalage est dû à la solubilité plus faible du nickel par rapport à celle du cobalt
(Bachet-2008, Plyasunova-1998a, Plyasunova-1998b), qui implique que la précipitation de la
première espèce ait lieu pour des pH moins basiques que la seconde. Notons également que,
dans les deux cas, le front de sorption débute à des pH plus élevés que dans le cas des ferrites,
où celui-ci apparaissait dès pH 5,5 – 6.
Pour ces systèmes, les analyses XPS n’ont pas été effectuées par manque de disponibilité de
l’appareil. Les conditions expérimentales étant identiques à celles utilisées pour l’étude des
systèmes ions/ferrites, il est probable que les précipités formés en surface soient Ni (OH ) 2 et
Co(OH ) 2 . La zétamétrie de la magnétite en présence de ces cations a été réalisée afin de
connaître l’évolution de la charge de surface avec la concentration en ions.
En l’absence de cations en solution, le point isoélectrique de la magnétite est de 6,4, ce qui est
en accord avec la littérature (Mansour-2007).
L’évolution du potentiel zêta en présence de nickel (figure V – 23) présente des similarités
avec celle en présence de cobalt, avec, cependant, un effet généralement plus important de la
concentration. En présence de chlorure de nickel à 10-6 M, aucun décalage significatif de la
variation du potentiel zêta n’est observé. De même, pour une concentration de 10-5 M jusqu’à
pH 6,5. Au-delà de ce pH, l’effet est plus marqué qu’en présence de cobalt pour les trois
concentrations supérieures (10-5 M, 5.10-5 M et 10-4 M), quant à la valeur du maximum
109
(excepté pour 5.10-5 M) et au décalage du PIE. Celui-ci varie de 9,8 à 10,5 lorsque la
concentration de nickel passe de 10-5 M à 10-4 M et se rapproche ainsi du PIE de Ni (OH ) 2
(11 – Parks-1965). L’effet plus marqué du nickel par rapport au cobalt peut être attribué à la
différence de solubilité des deux hydroxydes (Plyasunova-1998a, Plyasunova-1998b).
60
40
Potentiel zéta (mV)
20
0
sans CoCl
2
-6
10 M
-5
-20 10 M
-5
5.10 M
-4
10 M
-40
3 4 5 6 7 8 9 10 11
pH
Figure V - 22 : Variation du potentiel zêta de la magnétite en fonction du pH pour
différentes concentrations initiales en chlorure de cobalt (25°C, [KCl] = 10-4 M,
[CoCl2] = 10-4 M, 5.10-5 M, 10-5 M, 10-6 M et 0 M, [Fe3O4] = 0,4 g.L-1)
110
60
40
Potentiel zéta (mV)
20
0
sans NiCl
2
-6
10 M
-5
-20 10 M
5.10-5 M
-4
10 M
-40
4 5 6 7 8 9 10 11 12
pH
Figure V - 23 : Variation du potentiel zêta de la magnétite en fonction du pH pour
différentes concentrations initiales en chlorure de nickel (25°C, [KCl] = 10-4 M,
[NiCl2] = 10-4 M, 5.10-5 M, 10-5 M, 10-6 M et 0 M, [Fe3O4] = 0,4 g.L-1)
Les courbes de sorption ont été modélisées à l’aide du logiciel ECOSAT et sont représentées
sur les figures V - 24 pour le cobalt et V - 25 pour le nickel.
Plusieurs complexes ont été testés ( ≡ XOM + , (≡ XO) 2 M , ≡ XOMOH et ≡ XOM (OH ) −2 ) afin
de trouver le meilleur ajustement des courbes. Par analogie avec les systèmes ferrites, et
compte tenu des résultats de zétamétrie, les précipités choisis pour effectuer la modélisation
sont Ni(OH ) 2 et Co(OH ) 2 . Ceci est en accord avec Tewari et al. (Tewari-1972) qui
observent la formation d’une couche d’hydroxyde de cobalt à la surface de la magnétite à
25°C, pH 8 et pour une concentration initiale en cobalt de 10-6 M, par des expériences menées
en batch. A 200°C, tout autres paramètres égaux par ailleurs, l’apparition de l’hydroxyde de
cobalt est visible pour une concentration inférieure à 10-7 M.
La prise en compte d’un seul complexe pour réaliser les calculs n’a pas permis d’obtenir un
ajustement correct des courbes. Deux complexes ont donc été testés simultanément : ≡ XOM +
et (≡ XO) 2 M 8.
8
M représente soit l’ion cobalt soit l’ion nickel.
111
Pour une concentration de 10-6 M (figure V – 24), seul le phénomène de complexation est
présent. L’apparition d’un précipité de Co(OH ) 2 est visible dès 10-5 M, à partir de pH 9, et sa
quantité augmente avec la hausse de la concentration en chlorure de cobalt.
-6 -5
1.2 x 10 1.2 x 10
a +
b Co(OH)
-6
>XOCo -5
2
1 x 10 (>XO) Co 1 x 10
2 >XOCo +
modélisation (>XO) Co
2
expérimental
8 x 10-7
8 x 10-6 modélisation
(M)
expérimental
(M)
sorbé
sorbé
-7 -6
6 x 10 6 x 10
Co
Co
4 x 10-7 4 x 10-6
-7 -6
2 x 10 2 x 10
0 0
5 6 7 8 9 10 11 5 6 7 8 9 10 11
pH pH
-4
1.2 x 10
c
Co(OH) 2
-4
1 x 10 +
>XOCo
(>XO) 2Co
8 x 10
-5 modélisation
(M)
expérimental
sorbée
-5
6 x 10
Co
4 x 10-5
2 x 10-5
0
5 6 7 8 9 10 11
pH
Pour la sorption du nickel (figure V – 25) sur la magnétite le comportement est similaire à
celui en présence de cobalt (figure V – 24). Seul le phénomène de complexation est présent à
10-6 M puis il y a apparition d’un précipité d’hydroxyde de nickel, à partir de pH 8,5, pour les
concentrations supérieures.
112
-6 -5
1.2 x 10 1,2 x 10
a b
XONi+ Ni(OH)
2
(XO)2Ni -5
1 x 10-6 1 x 10 >XONi
+
modélisation
expérimental (>XO) Ni
2
-6 modélisation
8 x 10-7 8 x 10 expérimental
(M)
(M)
sorbé
sorbé
-6
6 x 10-7 6 x 10
Ni
Ni
-7 -6
4 x 10 4 x 10
-7 -6
2 x 10 2 x 10
0 0
4 5 6 7 8 9 10 5 6 7 8 9 10 11
pH pH
-4
1.2 x 10
Ni(OH)
c 2
1 x 10-4 >XONi
+
(>XO) Ni
2
-5 modélisation
8 x 10 expérimental
(M)
sorbé
6 x 10-5
Ni
-5
4 x 10
-5
2 x 10
0
4 5 6 7 8 9 10
pH
113
Tableau V - 3 : Logarithme des constantes de complexation et de solubilité, déterminées
à l’aide du logiciel ECOSAT (2-pK/BSM) pour les différents systèmes
Espèces
+
Concentrations ≡ XONi (≡ XO) 2 Ni Ni(OH ) 2 ≡ XOCo + (≡ XO) 2 Co Co(OH ) 2
Initiales (M)
10-6 -5 -10,5 / -5 -11 /
-5
10 -5,9 -10,5 -15,7 -5 -10,5 -15
-4
10 -5,7 -11,5 -15,2 -5,5 -10 -15
-5,45 ± -15,45 ± -5,25 ±
LogKmoyen -11 ± 0,5 -10,5 ± 0,5 -15
0,45 0,25 0,25
La sorption du cobalt sur la magnétite a été étudiée par Tamura et al. (Tamura-1983, Tamura-
1997a). Les complexes pris en compte dans leurs études sont identiques à ceux utilisés dans
celle-ci, cependant les constantes obtenues ont des valeurs différentes :
Il faut noter que les fronts de sorption n’apparaissent pas non plus aux mêmes valeurs de pH
(pH 6 pour Tamura et al., pH 7 pour cette étude). Toutefois, les conditions expérimentales des
deux études ne sont pas identiques. En effet, Tamura et al. ont fait varier la concentration en
cobalt de 1,7.10-8 M à 1,7.10-4 M et la force ionique de 10-4 M à 10-1 M en NaNO3 .
Dans le cas de la magnétite, trois espèces (un précipité et deux complexes) sont nécessaires à
la modélisation des courbes de sorption. En comparaison avec les ferrites, un complexe
bidentate neutre (≡ XO) 2 M (M = Ni ou Co), a été ajouté. La réactivité de surface des ferrites
semble donc différente de celle de la magnétite, comme il a été observé dans les courbes de
sorption et dans les mesures de zétamétrie.
On note que les constantes de précipitation des hydroxydes trouvées dans le cas de la
magnétite (environ -15 et -15,5) sont différentes des valeurs trouvées dans le cas des ferrites
(-18 et -17,5). Ceci pourrait être dû à des conditions différentes de formation des précipités,
ou bien, cela résulte simplement de la modélisation, qui, comme il a été vu, ne semble pas
prendre en compte tous les processus.
114
V.4. Mesures à 50°C
La figure V – 26 représente la cinétique de sorption des ions nickel sur le ferrite de cobalt à
50°C et pH 6 pour une concentration en chlorure de nickel de 10-6 M. Tout comme à 25°C, un
palier de sorption est atteint au bout d’environ 8 heures. Ce temps de contact sera également
pris en compte pour la réalisation des courbes de sorption.
Le palier correspond à une quantité de nickel sorbée d’environ 0,05 [Link]-2 à 50°C pour
seulement 0,04 [Link]-2 à 25°C (figure V – 1). Cette évolution peut être due à la formation
d’un précipité dès cette température. En effet, la solubilité du nickel diminue avec la hausse
de la température (Bachet-2008).
0.06
0.05
0.04
([Link] )
-2
0.03
sorbé
Ni
0.02
0.01
0
0 5 10 15 20 25
Temps de contact (h)
9
Des essais à 120°C et 320°C ont été réalisés, mais n’ont pas abouti (voir paragraphe V.5.).
115
Aucun décalage du front de sorption vers les pH acides par rapport à la courbe à 25°C n’est
visible. Une différence très faible peut cependant être présente mais peu visible du fait de
l’incertitude des mesures et de la faible différence de température. La présence d’un décalage
peut cependant être supposée pour une température plus élevée, du fait de l’évolution du PZC
en fonction de la température (Tableau V – 4).
-6
1 x 10
25°C
50°C
8 x 10-7
-7
(M)
6 x 10
sorbé
Ni
4 x 10-7
-7
2 x 10
0
3 4 5 6 7 8 9 10
pH
T
Le logiciel ECOSAT ne pouvant modéliser que des expériences à 25°C, il n’a pas été possible
de modéliser la courbe de sorption obtenue à 50°C.
Après 5 heures de contact entre le ferrite de nickel et les ions cobalt, l’équilibre est atteint.
Tout comme pour les cas précédents, le temps d’expérience sera de 8 heures. Contrairement à
la sorption du nickel sur le ferrite de cobalt à 50°C, aucune différence significative de la
quantité sorbée n’est visible sur la figure V – 28 (≈ 0,018 [Link]-2) par rapport aux mesures à
116
25°C. Ce qui est concordant avec le fait que la solubilité du cobalt est légèrement plus élevée
que celle du nickel (Baes-1986) et de l’influence de la température sur la solubilité.
0.02
([Link] )
-2
0.015
0.01
sorbé
Co
0.005
0
0 5 10 15 20 25
Temps de contact (h)
Les courbes de sorption en fonction du pH réalisées aux deux températures (25°C et 50°C)
(figure V – 29) sont très légèrement (0,2 unités de pH) décalées l’une par rapport à l’autre.
Ceci peut être dû à la diminution du PZC avec la hausse de la température (tableau V - 5).
Toutefois ce décalage vers les pH acides avec la hausse de la température étant très faible, il
n’est pas exclu que ce soit l’incertitude sur la mesure du pH qui en soit la cause.
-6
1.2 x 10
25°C
-6
1 x 10 50°C
8 x 10-7
(M)
sorbé
6 x 10-7
Co
4 x 10-7
-7
2 x 10
0
3 4 5 6 7 8 9 10
pH
T
117
Tableau V - 5 : PZC du ferrite de nickel en fonction de la température (milieu KCl
10-2 M) (Barale-2008)
T (°C) 0 10 25 40 50 70 90 125
PZC 8,7 8,4 8,2 8,1 7,7 7,4 7,5 7,2
V.4.3. Magnétite
Dans le but de connaître l’évolution du phénomène de sorption des ions cobalt et nickel sur la
magnétite en fonction de la température, les expériences ont été effectuées à 50°C avec une
concentration initiale en cobalt et en nickel de 10-6 M, afin d’éviter, ou du moins de limiter, le
phénomène de précipitation.
Les cinétiques de sorption des deux espèces sur la magnétite à 50°C sont représentées sur les
figures V – 30 et V – 31. Le palier de sorption est atteint en quelques heures, le même temps
de contact qu’à 25°C (8 heures) sera donc utilisé.
0.8
([Link] )
-2
0.6
sorbé
0.4
Co
0.2
0
0 1 2 3 4 5 6 7
Temps de contact (h)
Figure V - 30 : Evolution de la quantité de cobalt sorbée sur la magnétite en fonction du
temps (50°C, [KCl] = 10-4 M, [CoCl2] = 10-6 M, [Fe3O4] = 0,4 g.L-1, pH 6)
118
1
0,8
([Link] )
-2 0,6
sorbé
0,4
Ni
0,2
0
0 1 2 3 4 5 6 7
Temps de contact (h)
Les courbes de sorption en fonction du pH ont par la suite été effectuées, et sont représentées
sur les figures V - 32 et V - 33. Sur la figure V - 32, un décalage du front de sorption vers les
pH acides est visible. Ce résultat est en accord avec les données de la littérature (Tamura-
1983, figure III – 13c). De plus, la sorption des cations est favorisée pour pH > PZC. Le point
de charge nulle étant plus faible à 50°C qu’à 25°C (tableau V – 6), le décalage observé est
cohérent, cependant il semble trop grand pour être uniquement dû à l’évolution du PZC. Dans
les deux cas (cobalt et nickel) les incertitudes de mesures doivent entrer en jeu.
-6
1.2 x 10
50°C
-6 25°C
1 x 10
-7
8 x 10
(M)
sorbé
-7
6 x 10
Co
-7
4 x 10
-7
2 x 10
0
4 5 6 7 8 9 10 11 12
pH
T
: 25°C,
Figure V - 32 : Variation de la concentration en cobalt sorbé sur la magnétite (
: 50°C, [KCl] = 10 M, [CoCl2] = 10 M, [Fe3O4] = 0,4 g.L , 8 h)
-4 -6 -1
119
-6
1.2 x 10
50°C
1 x 10-6 25°C
-7
8 x 10
(M)
sorbé
6 x 10-7
Ni
4 x 10-7
-7
2 x 10
0
4 5 6 7 8 9 10 11 12
pH
T
: 25°C,
Figure V - 33 : Variation de la concentration en nickel sorbé sur la magnétite (
: 50°C, [KCl] = 10 M, [CoCl2] = 10 M, [Fe3O4] = 0,4 g.L , 8 h)
-4 -6 -1
Des expériences à plus hautes températures (120°C et 320°C) ont été testées dans le but d’être
le plus proche possible des conditions du circuit primaire. Des autoclaves permettant les
études de sorption ont été utilisés. Cependant des difficultés ont été rencontrées et n’ont pas
permis la réalisation de ces essais.
L’autoclave utilisé pour une température de 120°C étant réalisé en Téflon, il permet de rester
fortement inerte vis-à-vis de la solution et des particules étudiées. Toutefois, une sorption des
cations étudiés a été observée sur le système de filtration à chaud en acier inoxydable, et a
empêché d’effectuer les mesures de sorption sur les poudres.
D’autres essais ont été réalisés à 320°C, sur un autoclave en titane. Cependant, ce métal étant
fortement réactif (Schindler-1972) avec les cations, les mesures de sorption souhaitées n’ont
pas pu être effectuées.
120
V.6. Comparaison des trois systèmes
Les deux systèmes ferrites ont un comportement similaire à 25°C. En effet, pour une
concentration faible (10-6 M) en cation la complexation est le seul phénomène observé et
lorsque la concentration augmente (supérieure à 10-5 M), un précipité de surface (hydroxyde
de cobalt ou hydroxyde de nickel) se forme. Pour le système nickel sur ferrite de cobalt, la
présence en très faible quantité d’un précipité dès 10-5 M est possible mais difficile à
déterminer avec exactitude.
A 50 °C, le front de sorption de la courbe du système nickel sur ferrite de cobalt débute pour
des pH plus acides que celui à 25°C, ce qui peut être dû à l’évolution du PZC avec la
température. Cependant les comportements des deux systèmes sont similaires.
Il a été observé précédemment que les ions cobalt et nickel interagissent de façon similaire
avec la magnétite. Dans ce cas aussi l’analogie pour les expériences en température semble
faisable.
Lorsque l’on compare les ferrites avec la magnétite, des différences apparaissent. Les fronts
de sorption débutent à des pH plus acides avec les ferrites qu’avec la magnétite. Au cours de
la modélisation, pour le premier système, le meilleur ajustement des données expérimentales
est fait avec un seul un complexe de surface ( ≡ XOM + ) tandis que ≡ XOM + et (≡ XO) 2 M
121
sont nécessaires pour ajuster les données de la magnétite. Les constantes moyennes de
complexation et de solubilité sont données dans le tableau V – 7.
Il est observé que les réactivités de surface des ferrites, d’une part, et de la magnétite, d’autre
part, sont différentes, bien que les solides aient la même structure cristalline. Cette différence
est probablement due à la présence d’atomes de cobalt ou de nickel dans le réseau cristallin
des ferrites.
Afin de voir si une plus grande quantité d’oxyde tend à favoriser la complexation de surface,
des expériences avec une concentration en nickel ou en cobalt de 10-5 M avec 4 g.L-1, en
complément aux expériences avec 0,4 g.L-1, ont été effectuées (figures V – 34 et V – 35).
122
1
-1 -6
0,4 g.L / 10 M
0.8 -1 -5
4 g.L / 10 M
([Link] )
-2
0.6
sorbé
0.4
Ni
0.2
0
4 5 6 7 8 9 10 11 12
pH
Figure V - 34 : Evolution de la quantité de nickel sorbée sur la magnétite en fonction du
pH (25°C, [KCl] = 10-4 M, [NiCl2] = 10-6 M / [Fe3O4] = 0,4 et [NiCl2] = 10-5 M /
[Fe3O4] = 4 g.L-1, 8 h)
0.6
sorbé
0.4
Co
0.2
0
4 5 6 7 8 9 10 11 12
pH
Figure V - 35 : Evolution de la quantité de cobalt sorbée sur la magnétite en fonction du
pH (25°C, [KCl] = 10-4 M, [CoCl2] = 10-6 M / [Fe3O4] = 0,4 et [CoCl2] = 10-5 M /
[Fe3O4] = 4 g.L-1, 8 h)
123
Sur les figures V – 34 et V – 35, les courbes pour lesquelles la concentration en magnétite est
de 4 g.L-1 et celle en chlorure de nickel ou de cobalt de 10-6 M avec 0,4 g.L-1, atteignent une
quantité d’atomes sorbés équivalente (≈ 0,9 [Link]-2 pour les deux cations). La valeur du palier
correspond au fait que 100 % des espèces en solution sont sorbées, comme il est visible sur
les figures V – 19, V – 21, V – 24 et V – 25. Un décalage de la courbe réalisée à 10-5 M /
4 g.L-1 vers les pH acides par rapport à celle à 10-6 M / 0,4 g.L-1 est visible. Ceci montre
qu’une plus grande quantité de solide favorise la sorption à des pH moins élevés, donc moins
favorables à la précipitation.
10
Notons que le pH augmente avec la diminution de la température. Pour un conditionnement de 2000
ppm en bore et 1,5 ppm en lithine le pH à 300°C est de 5,7, ce qui correspond à un pH de 6,5 à 25°C.
124
VI. Comportement des particules colloïdales dans le circuit
primaire des centrales REP
Dans les chapitres précédents, les interactions entre les ions cobalt et nickel et des particules
modèles de certains produits de corrosion (ferrite de cobalt, ferrite de nickel, magnétite) ont
été étudiées.
Dans cette partie, il s’agit de s’intéresser aux interactions entre ces mêmes produits de
corrosion et les parois du circuit primaire. En effet, les interactions entre les produits de
corrosion en suspension dans le fluide primaire et les matériaux constitutifs de celui-ci jouent
un rôle essentiel dans la contamination surfacique du circuit. Pour cela, un travail en commun
avec le département MFEE de EDF R&D a été initié afin de prendre en compte l’hydraulique
et donc la cinétique des particules. Le département MFEE a en effet développé dans son code
d’écoulement hydraulique (Code_Saturne) (Saturne-2007), un module permettant de prendre
en compte le dépôt de particules (module Lagrangien). Ce module nécessite la connaissance
du diamètre des particules présentes dans le fluide, de leur charge de surface ainsi que de celle
des matériaux métalliques du circuit.
Dans ce chapitre, le module Lagrangien sera présenté, ainsi que les résultats obtenus sur la
mesure des charges de surface des matériaux représentatifs du circuit primaire des centrales
REP.
Lors du transport, les particules sont en mouvement dans le fluide et sont ainsi plus ou moins
proches des parois métalliques. La probabilité qu’une particule se dépose à la surface de
celles-ci est d’autant plus forte que la distance particule/parois est faible. Cependant, les
interactions électrostatiques interviennent également dans le dépôt des particules. La théorie
DLVO (Derjaguin, Landau, Verwey and Overbeek) permet de prendre en compte ce
phénomène.
125
Figure VI - 1 : Schéma de l’évolution du potentiel en fonction de la distance
paroi/particule (partie haute du schéma : charges de même signe, partie basse du
schéma : charges de signes opposés) (Cordier-2007) (ordonnée : énergie d’interaction)
Selon la théorie DLVO, l’énergie potentielle d’interaction Ep, entre une particule et une paroi
séparées par une distance d, est constituée de la somme de l’énergie d’interaction de Van der
Waals, toujours attractive, et d’une contribution électrostatique, qui peut être soit attractive
soit répulsive suivant le signe des potentiels de surface des deux partenaires. La modélisation
de la théorie DLVO permet donc de déterminer en fonction de la distance à la paroi si, pour
un pH donné, les conditions entre particules et paroi sont plutôt attractives ou plutôt
répulsives. Le potentiel total d’interaction Ep est donné par les formules suivantes :
Ep A1w 2
=− + 4πεε 0ϕ1ϕ 2 e− d λ Equation VI - 1
Rp 6d
−1 2
2000F 2 I
λ = Equation VI - 2
εε 0 RT
En se basant uniquement sur cette théorie, donc sur les charges de surface, il pourrait être
possible de prédire le dépôt des produits de corrosion sur les matériaux métalliques
constitutifs du circuit primaire.
126
Par exemple, l’encrassement d’échangeurs de chaleur étant souvent un souci dans l’industrie,
le dépôt de particules de TiO2 sur des plaques ondulées en acier inoxydable, a été étudié par
des expériences en boucle d’essais (Grandgeorge-1998). Il a été observé que lorsque le
phénomène d’encrassement est dès le départ limité par un mécanisme d’adhésion, il le restera.
Toutefois, s’il est limité par le transfert de masse, donc par la vitesse du fluide, une transition
est visible lorsque les particules accumulées sur l’acier induisent un coefficient d’adhésion
significatif.
Des travaux menés sur l’adhésion de particules d’hématite sur de l’acier inoxydable 316L (les
particules d’hématite traversent un « lit » de particules d’acier) en fonction du pH et à 25°C
ont montré un taux de dépôt élevé de pH 6 à pH 8 (Kuo-1980). Par comparaison de la fraction
de particules ayant adhéré à la surface de l’acier avec la mobilité électrophorétique des deux
matériaux, il apparaît que l’adhésion est plus forte lorsque les mobilités électrophorétiques
sont de signes opposés, c’est-à-dire lorsque les charges des deux matériaux sont de signes
différents (figure VI – 2).
Ryde et al. (Ryde-2000), ont eux aussi montré l’influence de la vitesse du fluide dans le cas
du dépôt de particules d’hématite sur du verre.
Une étude menée sur le dépôt de particules d’hématite sur des parois en polypropylène en
conditions dynamiques a montré que la simple considération du signe de la charge des
particules et des parois n’expliquait pas l’adhésion dans tous les cas (Cerovic-2009). Les
127
expériences en conditions dynamiques, menées dans une boucle d’essai, ont été réalisées dans
les conditions suivantes :
- 25°C,
- particules d’hématite de 320 nm de diamètre,
- gamme de pH : 4-11,
- nombre de Reynolds : 3300-1770011.
Lors de ces essais il a été montré que l’adhésion a principalement lieu lorsque les conditions
sont attractives. Cependant, dans des conditions répulsives, l’adhésion a également été
observée en de plus faibles proportions. De plus, plus la vitesse du fluide était grande, plus
l’adhésion en conditions répulsives était importante.
Ces études démontrent bien que l’hydraulique doit être pris en compte pour évaluer le
comportement des produits de corrosion dans le fluide primaire.
Afin de prendre en compte l’hydraulique lors des prédictions d’adhésion des produits de
corrosion sur les parties métalliques du circuit primaire, le module Lagrangien a été utilisé
pour suivre le cheminement des particules dans l’écoulement (Guingo-2008).
Les conditions prises en compte pour les calculs sont les suivantes12 :
- Vitesse du fluide : 5,73 m.s-1,
- Diamètre de la conduite : 1,69.10-2 m,
- Distance entre les sections 1 et 2 : 10 cm,
- Nombre de particules injectées à l’entrée de la conduite : 100 000,
- Pas de temps : 3.10-6 s,
- Température du fluide : 315°C,
- Température de la conduite : 306,4°C.
Pour la réalisation des calculs, il a été fixé que lorsqu’une particule entre en contact avec la
paroi, elle est définitivement déposée. Le phénomène de remise en suspension n’a pas été pris
11
Notons qu’un nombre de Reynolds compris dans cette gamme correspond à un écoulement turbulent
12
Les conditions sont choisies les plus proches possible de celle de l’écoulement dans les tubes GV.
128
en compte. Les données d’entrée nécessaires pour réaliser les calculs sont les charges de
surface des particules et des parois, ainsi que le diamètre des particules.
Dans un premier temps, il convient de déterminer le temps minimum que mettent les
particules avant d’être en équilibre dynamique avec le fluide. Ce paramètre, appelé temps de
relaxation, est noté τ p+ et est une valeur adimensionnée (équation VI – 3).
d p2 uτ2 ρp
τ = +
× Equation VI - 3
p
18ν 2
f ρf
Dans un deuxième temps, il est nécessaire de déterminer le taux de dépôt des particules, qui
est représenté par k p (équation VI – 4).
mp U m dc N
kp = = × ln 1 Equation VI - 4
C 4 (x 2 − x1 ) N2
où :
- m p est la masse de particules déposées par seconde et par unité de surface,
- C la masse de particules dans le fluide,
- U m la vitesse de l’écoulement,
- dc le diamètre de la tuyauterie,
- x1 et x 2 les positions des sections entre lesquelles est calculé le taux de dépôt,
- N1 et N 2 le nombre de particules encore en suspension aux sections 1 et 2
De façon à avoir un taux de dépôt normalisé, k p est divisé par la vitesse de frottement pour
obtenir k p+ (équation VI – 5).
kp
k p+ = Equation VI - 5
uτ
Une compilation des données expérimentales disponibles dans la littérature est donnée dans la
figure VI – 3. Celle-ci représente l’évolution du taux de dépôt en fonction du temps de
relaxation. La figure VI – 3 montre une augmentation significative du taux de dépôt avec
l’inertie des particules.
129
Figure VI - 3 : Compilation des données expérimentales disponibles dans la littérature
des taux de dépôt des particules en fonction du temps de relaxation (Papavergos-1984)
Pour réaliser les calculs, il a donc été nécessaire de déterminer les charges de surface des
matériaux massifs. Pour cela deux méthodes ont été utilisées : la zétamétrie (par mesure du
potentiel d’écoulement) et la mesure de la vitesse d’adhésion de particules de latex chargées à
la surface des matériaux.
Tewari et al. (Tewari-1972) se sont intéressés au potentiel zêta du Zircaloy. A 25°C (pH 10,
électrolyte : LiOH ) celui-ci est de 90 mV tandis qu’il est de 650 mV à 200°C. Sa surface est
donc chargée très positivement en température ce qui favorise la possibilité d’interaction avec
des espèces chargées négativement présentes dans le fluide primaire.
L’effet de la composition des aciers inoxydables ainsi que de leur état de surface sur le point
isoélectrique, a été déterminé (Boulangé-Petermann-1995). Le PIE obtenu dans cette étude est
de 3,6 pour les aciers inoxydables « bruts13 » et de 4,2 pour les surfaces polies. Toutefois le
PIE des échantillons massifs n’est pas identique à celui des oxydes dispersés en solution
(Boulangé-Petermann-1995). La figure VI – 4 représente la variation du potentiel zêta de
l’acier inoxydable 304L en fonction du pH pour trois échantillons ayant subi des traitements
différents (Tableau VI – 1).
13
Aciers inoxydables non polis, utilisés dans leur état de réception
130
Figure VI - 4 : Potentiel zêta et PIE de l’acier inoxydable 304L ayant subi différents
traitement de surface (0,01 M NaCl, : Echantillon 1 brut, : Echantillon 2 brut, :
Echantillon 2 poli) (Boulangé-Petermann-1995)
131
Figure VI - 5 : Influence du pH sur la concentration en particules de latex dans une
suspension mise en contact avec de l’acier inoxydable 304L ( ), du Zircaloy 4
), 316L (
( ) - (Lefèvre-2009)
) et de l’Inconel 690 (
Les valeurs de point isoélectrique obtenues pour les quatre matériaux étudiés (aciers
inoxydables 304L, 316L, Zircaloy 4 et Inconel 690) sont comprises entre 2,4 et 3,4. Il est
intéressant de noter que les PIE des alliages recouverts d’une légère couche de passivation,
sont plus faibles que ceux mesurés directement sur les oxydes sous forme de poudre
(respectivement 2,4 – 3,4 contre 6 – 11). Le phénomène de charge des surfaces de matériaux
métalliques est proche de celui des surfaces hydrophobes (or, Téflon…), où l’accumulation
d’ions hydroxydes mènerait à un potentiel de surface négatif sur une grande gamme de pH (3
– 14) (Lefèvre-2009).
Kallay et al. ont, avec ce même type de particules de latex, déterminé le PIE du cuivre (9,8)
(Kallay-1987) à partir d’une concentration en particules de latex de 8,8 mg.L-1 dans une
solution à 2.10-3 M de NaOH . Celles-ci ont également été utilisées ainsi qu’un second type
de particules (polystyrène Amidine ayant des groupements NH2 en surface, la rendant
positive) afin de définir les PIE de plusieurs métaux tels que l’étain ou le platine (Kallay-
1991). Le tableau VI – 2 résume ces valeurs de point isoélectriques. Par la suite, les deux
types de latex seront notés Amidine et Carboxyle.
132
Tableau VI - 2 : Points isoélectriques de métaux obtenus par mesure de la vitesse
d’adhésion de particules de latex à 20°C (Concentration en particules de latex :
8,8 mg.L-1, 10-3 M en NaOH, rAmidine = 38 nm, rCarboxyle = 148 nm) (Kallay-1991)
Métal Pureté (%) pH initial Charge des particules de latex PIE
(- : Carboxyle, + : Amidine)
Mo 99,95 11 - 2,8 ± 0,15
Pt 99,98 11 - 3,3 ± 0,15
Te 99,5 11 - 3,3 ± 0,15
Si 99,5 11 - 3,7 ± 1
Sn 99,8 11 - 4,3 ± 0,15
3 + 4,3 ± 0,15
Cr 99,95 11 - 6,7 ± 0,15
3 + 6,6 ± 0,5
Co 99,8 11 _ 5±2
3 + 8,9 ± 0,15
Ni 99,99 11 - 11 ± 0,15
3 + 7 ± 0,15
Dans cette étude, la différence de conditionnement des métaux influe sur le PIE. Pour le
nickel par exemple, après un conditionnement à pH 11, le point isoélectrique est de 11 ce qui
correspond à celui de Ni (OH ) 2 , tandis que pour un pH initial de 3 la valeur est similaire à
celle de NiO . La surface est donc bien sensible au milieu dans lequel se trouve l’échantillon.
Kleimann et al. (Kleimann-2006) arrivent à la conclusion que la nature du substrat a une
influence sur le dépôt des particules chargées.
Deux méthodes ont été utilisées pour déterminer la charge de surface des matériaux choisis
(Inconel 690, Zircaloy 4, acier inoxydable 304L). La première est la zétamétrie par mesure du
potentiel d’écoulement (décrite dans le chapitre IV), la seconde est l’adhésion de particules de
latex chargées sur la surface des matériaux métalliques.
133
VI.3.2.1. Mesure du potentiel d’écoulement
Comme décrit dans le chapitre IV.4, deux échantillons de chaque matériau (Zircaloy 4,
Inconel 690 et acier inoxydable 304L) ont été vieillis dans le but d’évaluer l’influence du
vieillissement sur le potentiel de surface des matériaux.
Rappelons les conditions de vieillissement en autoclave, prises le plus proche possible des
conditions de fonctionnement du réacteur :
- 80 jours,
- 360°C (température plus élevée que celle du réacteur afin d’accélérer la formation des films
d’oxydes en surface des matériaux),
- 1000 ppm de bore,
- 2 ppm de lithine,
- 30 [Link]-1 d’hydrogène.
La mesure par potentiel d’écoulement est très sensible à la rugosité du matériau analysé, en
particulier à ses extrémités. En effet, l’électrolyte passant entre deux plaques du matériau
distantes de seulement 100 µm, la moindre aspérité lors de l’entrée en contact du fluide avec
le matériau, a un effet sur le flux de liquide. Pour s’affranchir de ce phénomène, chaque
mesure a été effectuée de façon à ce que le fluide passe dans un sens puis dans l’autre entre
les deux plaques.
La figure VI - 6 donne un exemple des courbes obtenues par cette méthode, pour l’Inconel
690 avant et après vieillissement en autoclave. Le tableau VI - 3 résume les PIE des différents
matériaux obtenus par cette méthode, et les courbes correspondantes pour les autres matériaux
testés sont données en annexe 6.
40
vieilli
non vieilli
20
Potentiel zéta (mV)
-20
-40
-60
3 3.5 4 4.5 5 5.5 6 6.5 7
pH
134
Tableau VI - 3 : PIE de l’Inconel 690, du Zircaloy 4 et de l’acier inoxydable 304L avant
et après vieillissement en autoclave, obtenus par zétamétrie
Inconel 690 Zircaloy 4 Acier inoxydable 304L
Matériau non vieilli 4,2 ± 0,1 / 4,5 ± 0,1
Matériau vieilli 4,3 ± 0,1 3,8 ± 0,1 /
Le tableau VI – 3 résume les valeurs des points isoélectriques des trois matériaux métalliques
étudiés. Pour l’Inconel 690, le PIE est le même pour le matériau oxydé et le non oxydé. Pour
les deux autres matériaux les PIE déterminés sont proches de ceux de l’Inconel 690 (≈ 4). Il
semble donc qu’il n’y ait pas d’influence de la présence d’un film d’oxyde, dans le cas de
l’Inconel 690, à la surface du matériau sur les valeurs de point isoélectrique.
Cette méthode étant très sensible à la rugosité des échantillons et n’ayant pas pu être réalisée
pour tous les matériaux, une seconde méthode a été utilisée pour déterminer les PIE
manquants. L’adhésion de particules de latex sur les solides a été choisie car c’est une
méthode plus simple à mettre en œuvre mais également moins coûteuse que la zétamétrie,
bien que moins précise.
La seconde méthode qui a été mise en œuvre afin de confirmer les résultats obtenus, est la
mesure de la vitesse d’adhésion de particules chargées (dans ce cas des particules de latex :
Amidine et Carboxyle) à la surface des trois matériaux métalliques (vieillis et non vieillis).
Par adhésion sur la surface, de particules de latex également chargées mais de façon uniforme
en fonction du pH, il est possible de déterminer le potentiel de surface du solide étudié. Pour
cela, deux types de latex sous forme de particules ont été choisis, l’Amidine dont la charge est
positive jusqu’à pH ≈ 11 et le Carboxyle dont la charge est négative dès pH ≈ 2. Le protocole
suivi pour ces expériences est le suivant :
- Les échantillons massifs (identiques à ceux utilisés pour la zétamétrie par potentiel
d’écoulement) sont mis dans 10 mL de solution (10-3 M en NaNO3 ),
- Le pH est ajusté à la valeur souhaitée, dans ce cas 5,5 et 3,3,
- Les particules de latex sont ajoutées à la solution (200 mg.L-1)14,
- Le mélange est ensuite mis sous agitation à température ambiante durant 96 heures,
- Les échantillons peuvent ainsi être rincés avec une solution de même pH, ne contenant pas
de particules de latex, dans le but de n’avoir à la surface plus que celles ayant réagi avec la
surface du matériau.
Une fois séchés, les échantillons ont été observés au MEB. Les figures VI – 7 et VI – 8
donnent les résultats obtenus pour l’Inconel 690 avec les particules de Carboxyle. Les
observations MEB des deux autres matériaux sont données en annexe 7.
14
Cette concentration assez élevée en comparaison aux données disponibles dans la littérature, permet
d’avoir un temps de contact moins long et une agitation plus faible.
135
Sur les images MEB obtenues pour les échantillons étudiés à pH 3,3 (figure VI – 7 : a-b,
figure VI – 8 : a-b), les particules de latex sont visibles, tandis qu’elles ne le sont pas sur les
échantillons à pH 5,5 (figure VI – 7 : c, figure VI – 8 : c). Pour chaque matériau (Inconel 690,
acier inoxydable 304L et Zircaloy 4) l’adhésion des particules de latex n’a lieu qu’à un des
deux pH (à pH 3,3 pour le Carboxyle et à pH 5,5 pour l’Amidine). Sur les différents clichés
MEB où les particules de latex sont visibles, aucune adhésion préférentielle sur les cristaux,
du fait de leur taille, n’est observée.
a b
1 µm 1 µm
1 µm
Figure VI - 7 : Observations MEB de l’Inconel 690 non vieilli après mise en contact avec
des particules de Carboxyle (a-b : pH 3,3, c : pH 5,5)
136
a b
1 µm
1 µm
1 µm
Figure VI - 8 : Observations MEB de l’Inconel 690 vieilli après mise en contact avec des
particules Carboxyle (a-b : pH 3,3, c : pH 5,5)
Les résultats obtenus par cette méthode sont en accord avec les points isoélectriques
déterminés à l’aide de la zétamétrie. Il semble donc bien que la présence d’un film d’oxyde en
surface du matériau influe peu sur la valeur de son PIE.
A partir de ces résultats, il est alors possible de tenter une prédiction du comportement des
produits de corrosion dans le circuit primaire par les calculs réalisés à l’aide du module
Lagrangien de Code_Saturne.
137
VI.4.1. Récapitulatif des différents PIE
Les tableaux VI – 4 et VI – 5 résument les différents PIE déterminés dans cette thèse, et qui
seront utilisés dans les calculs.
Tableau VI - 4 : Récapitulatif des PIE des matériaux métalliques (Inconel 690, Zircaloy
4, acier inoxydable 304L) déterminés dans la thèse
Inconel 690 Zircaloy 4 Acier inoxydable 304L
Matériau non vieilli 4,2 ± 0,1 / 4,5 ± 0,1
Matériau vieilli 4,2 ± 0,1 3,8 ± 0,1 /
0 6,9 ± 0,1
6,9 ± 0,1
-6
CoFe2O4 10
10-5 6,9 ± 0,1 ///
5.10-5 11
10-4 11
0 9,1 ± 0,1
10-6 9,6 ± 0,1
NiFe2O4 10-5 /// 9,6 ± 0,1
5.10-5 10,5 ± 0,1
10-4 10,5 ± 0,1
0 6,4 ± 0,1 6,4 ± 0,1
10-6 6,4 ± 0,1 6,4 ± 0,1
Fe3O4 10-5 10,1 ± 0,1 6,4 ± 0,1
5.10-5 10,5 ± 0,1 9,8 ± 0,1
10-4 10,5 ± 0,1 10,2 ± 0,1
Les diamètres des produits de corrosion présents dans le circuit primaire n’étant pas connus
précisément, plusieurs valeurs ont été retenues pour les calculs (50, 200 et 400 nm).
Le tableau VI – 6 résume, pour les différents calculs réalisés avec le module Lagrangien de
Code_Saturne, le diamètre des particules, les potentiels de surface de la paroi et des
particules15 et le taux de dépôt normalisé k p+ calculé entre les sections 1 et 2 distantes de
10 cm.
15
ϕ1 et ϕ 2 sont attribués indifféremment à la paroi ou aux particules
138
Tableau VI - 6 : Paramètres choisis pour réaliser les calculs avec le module Lagrangien
et taux de dépôt normalisés en résultant (tous les résultats pour les charges de même
signe sont donnés pour des potentiels (ϕi) négatifs, cependant pour les valeurs positives,
les taux de dépôt sont identiques)
Diamètre (nm) ϕ1 (mV) ϕ 2 (mV) k p+
50 -1 +1 1,69.10-2
-1 -1 1,69.10-2
-10 -10 1,69.10-2
-20 -20 2,07.10-3
-26 -26 1,39.10-5
-27 -27 0
200 -10 +10 8,64.10-3
-30 +30 8,64.10-3
-1 -1 8,64.10-3
-10 -10 8,64.10-3
-17 -17 7,66.10-5
-18 -18 0
-30 -30 0
400 -10 -10 7,69.10-3
-12 -12 7,69.10-3
-20 -20 0
- l’adhésion peut avoir lieu même lorsque les signes des charges de surface sont identiques.
- pour un même diamètre de particule et des charges de surface de signes identiques, le taux
de dépôt diminue lorsque la valeur absolue des potentiels augmente.
- pour un même diamètre de particule et des charges de surface de signes opposés faibles en
valeur absolue (inférieur à un certain seuil), le taux de dépôt est égal à celui obtenu pour des
potentiels de même signe. Lorsque les potentiels sont au-dessus de ce seuil, les valeurs de k p+
sont différentes.
- pour des diamètres de particule différents mais des charges égales, le taux de dépôt diminue
avec l’augmentation du diamètre.
Il semble donc que les petites particules se déposent plus rapidement que les grosses. Cette
observation est cohérente avec le fait que lorsque les particules ont une taille inférieure au
micron, ce qui est le cas pour ces calculs, le phénomène permettant leur dépôt est le
mouvement Brownien, et ce dernier est d’autant plus intense que les particules sont petites.
139
Figure VI - 9 : Compilation des données expérimentales disponibles dans la littérature
(.) des taux de dépôt des particules en fonction du temps de relaxation (Papavergos-
1984), : Calculs réalisés pour différents paramètres avec le modèle utilisé dans la thèse
(Guingo-2008)
Pour les paramètres choisis dans ces calculs, le temps de relaxation est environ égal à 1
(figure VI – 9). Le modèle défini par Guingo et al. (Guingo-2008) donne des résultats en
bonne corrélation avec les données de la littérature. Etant donné que le modèle utilisé dans la
thèse est celui-ci, les taux de dépôt déterminés semblent concordants avec l’allure des courbes
de la figure VI – 9.
En effet, pour des particules de faible inertie, donc de diamètre faible (inférieur au
micromètre), le taux de dépôt k p+ augmente lorsque le temps de relaxation τ p+ diminue. Selon
l’équation VI – 3, τ p+ est proportionnel au diamètre de particules. Donc k p+ augmente avec la
diminution du diamètre des particules injectées dans l’écoulement.
Selon une étude précédente (Barale-2006), en tenant compte uniquement des PIE des oxydes
(ferrites de cobalt et de nickel et magnétite) et des surfaces du circuit primaire (Inconel 690,
Zircone), il a été prédit que le dépôt des oxydes sur les tubes des générateurs de vapeur est
favorisé contrairement à celui sur les gaines du combustible.
Les points isoélectriques des principaux matériaux constitutifs du circuit primaire ont été
déterminés à l’aide de deux méthodes distinctes. Les valeurs des PIE sont, dans le cas des
deux techniques, proches de 4 et subissent très peu l’influence de l’épaisseur du film d’oxyde
en surface.
140
A partir des différents résultats obtenus au cours de la thèse, une première approche du
comportement des produits de corrosion dans le circuit primaire est possible.
Par zétamétrie, il a été possible d’obtenir les points isoélectriques des trois oxydes étudiés : le
ferrite de nickel, le ferrite de cobalt et la magnétite, en fonction de la quantité d’ions nickel ou
cobalt en solution (chapitre V). Les PIE de trois matériaux métalliques constitutifs du circuit
primaires des réacteurs REP ont également été déterminés dans ce chapitre. Par comparaison
de ces différentes données expérimentales avec les calculs effectués avec le module
Lagrangien, il est possible de tenter une prédiction qualitative, sur le comportement de ces
produits de corrosion dans le fluide primaire. Cependant, ces résultats doivent être analysés
avec précaution pour deux raisons majeures. Premièrement, les expériences ont été effectuées
à 25°C, ce qui est loin de la température du fluide primaire en fonctionnement.
Deuxièmement, le modèle utilisé par le module Lagrangien n’est pas entièrement validé.
sans NiCl
2
50 -4
10 M
I690 vieilli
I690 non vieilli
25
Potentiel zéta (mV)
-25
-50
2 4 6 8 10 12
pH
Figure VI - 10 : Evolution du potentiel zêta et du PIE de l’Inconel 690 avant et après
vieillissement ([NaNO3] = 10-3 M) en autoclave et du ferrite de cobalt ([KCl] = 10-4 M) en
présence de chlorure de nickel en fonction du pH ([CoFe2O4] = 0,4 g.L-1, [NiCl2] = 0 M et
10-4 M, 25°C)
Les figures VI – 10, VI – 11 et VI – 12 représentent les évolutions des potentiels zêta des
produits de corrosion et de l’Inconel 690 (noté I690) avant et après vieillissement en fonction
du pH. Les courbes correspondant aux deux autres matériaux sont données en annexe 8. Sur
141
les trois figures les lignes interrompues horizontales vertes et bleues correspondent aux
potentiels limites au-delà desquels le dépôt n’a plus lieu si les charges de surface des deux
éléments sont de même signe. Les lignes bleues correspondent à un diamètre de particule de
200 nm, les vertes à un diamètre de 50 nm. Pour les trois produits de corrosion, seules les
concentrations de 0 M et 10-4 M en chlorure de nickel ou de cobalt ont été représentées, et ce
par souci de lisibilité des courbes.
En se basant sur les résultats rapportés dans le tableau VI – 7, il semble que le dépôt du ferrite
de cobalt sur l’Inconel 690 soit très largement favorisé quelque soit le pH. Seuls des pH très
acides (< 3,7) dans le cas du matériau non vieilli pourraient empêcher ce dépôt, ou très
basique (> 9,2) pour le ferrite n’ayant pas subi la sorption du nickel.
142
40
20
-20
sans CoCl2
-40
-4
10 M
I690 vieilli
I690 non vieilli
-60
2 4 6 8 10 12
pH
25
Potentiel zéta (mV)
-25
-50
2 4 6 8 10 12
pH
143
Tableau VI - 8 : Détermination des gammes de pH pour lesquelles le dépôt de particules
de ferrite de nickel, avant ou après sorption de chlorure de cobalt (0 M et 10-4 M), sur
l’Inconel 690 vieilli (V) ou non (NV) en autoclave, a lieu (diamètre de la particule :
50 nm). O signifie que le dépôt est favorisé, N que non (ND (non défini) : dépend de la
continuité des courbes)
pH < 4,1 4,1 - 6 6 – 9,2 9,2 - 11 > 11 / /
I690 V
NiFe2O4 O O O O ND /
0M
9,2 –
pH < 4,1 4,1 - 6 6 – 7,5 7,5 – 9,2 > 10,5 /
10,5
I690 V
NiFe2O4 O O O O O O /
-4
10 M
pH < 3,7 3,7 – 4,2 4,2 – 4,9 4,9 - 6 6 – 9,2 9,2 - 11 > 11
I690 NV
NiFe2O4 N O O O O O N
0M
7,5 –
pH < 3,7 3,7 – 4,2 4,2 – 4,9 4,9 – 7,5 > 10,5 /
10,5
I690 NV
NiFe2O4 N O O O O O /
-4
10 M
Dans le cas du ferrite de nickel en présence de l’Inconel 690 (tableau VI – 8), la tendance est
la même que pour le ferrite de cobalt. Le dépôt n’a pas lieu, pour le matériau non vieilli, pour
des pH très acides (< 3,7), voire très basiques (> 11) lorsque le ferrite n’a pas subi la sorption
de chlorure de cobalt. Ce comportement est identique à celui observé auparavant pour le
ferrite de cobalt.
Lors de la mise en contact de particules de magnétite avec l’Inconel 690 (tableau VI – 9), le
dépôt a lieu sur une très large gamme de pH, de façon similaire, ce qui a été observé pour les
ferrites de nickel et de cobalt. Lorsque l’Inconel a été vieilli dans les conditions du circuit
primaire, la prédiction indique que le dépôt n’a lieu que pour la magnétite seule et à pH
compris entre 8 et 9,2. Pour le matériau non vieilli, comme pour les ferrites, le dépôt a lieu
pour des pH très acides (< 3,7).
144
Tableau VI - 9 : Détermination des gammes de pH pour lesquelles le dépôt de particules
de magnétite, avant ou après sorption de chlorure de cobalt ou de nickel (0 M et 10-4 M),
sur l’Inconel 690 vieilli (V) ou non (NV) en autoclave, a lieu (diamètre de la particule :
50 nm). O signifie que le dépôt est favorisé, N que non (ND (non défini) : dépend de la
continuité des courbes)
pH < 4,1 4,1 – 5,6 5,6 - 6 6 – 6,4 6,4 - 8 8 – 9,2 > 9,2
I690 V Fe3O4
O O O O O N ND
0M
pH < 4,1 4,1 - 6 6 – 9,2 9,2 – 9,5 9,5 - 11 > 11 /
I690 V Fe3O4
O O O O O ND /
10-4 M NiCl2
pH < 4,1 4,1 – 5,8 5,8 - 6 6 – 9,2 9,2 – 10,5 > 10,5 /
I690 V Fe3O4
O O O O O O /
10-4 M CoCl2
pH < 3,7 3,7 – 4,2 4,2 – 4,9 4,9 – 5,6 5,6 – 6,4 6,4 - 8 >8
I690 NV
Fe3O4 N O O O O O ND
0M
pH < 3,7 3,7 – 4,2 4,2 – 4,9 4,9 – 9,5 9,5 - 11 > 11 /
I690 NV
Fe3O4 N O O O O O /
-4
10 M NiCl2
pH < 3,7 3,7 – 4,2 4,2 – 4,9 4,9 – 5,8 5,8 – 10,5 > 10,5 /
I690 NV
Fe3O4 N O O O O O /
-4
10 M CoCl2
Dans le cas de l’Inconel 690, il semble que le dépôt soit défavorisé sur les surfaces non
vieillies et à pH inférieur à 3,7, ou dans le cas du nickel uniquement, sur les surfaces vieillies
pour un pH élevé (> 10). Cependant, les mesures de potentiel d’écoulement étant sensibles à
la force ionique du milieu, la détermination des PIE à différentes concentrations en sel de
fond, permettrait d’obtenir une meilleure prédiction du comportement des produits de
corrosion dans le fluide primaire. En comparaison avec l’étude menée précédemment et ne
tenant pas compte de l’hydraulique (Barale-2006), le dépôt des oxydes sur les parois du
circuit primaire semble favorisé pour les trois matériaux étudiés (Inconel 690, Zircaloy 4 et
acier inoxydable 304L), ce qui est donc en accord avec les dépôts observés en centrale.
La chimie du fluide primaire est très peu modifiable. En fonctionnement, le pH à 300°C est de
7,2. Il semble donc qu’à ce pH, le dépôt des trois produits de corrosion étudiés, ayant subi la
sorption du cobalt ou du nickel, soit possible, et qu’il est donc très difficile de modifier la
composition chimique pour limiter ce phénomène en phase de fonctionnement.
145
146
VII. Conclusion générale et perspectives
Dans un premier temps, les mécanismes de sorption des cations (cobalt et nickel) sur les trois
oxydes représentatifs ont été étudiés. La variation de la quantité de cations sorbée sur les
particules colloïdales en fonction du pH, modélisée à l’aide du logiciel ECOSAT, les analyses
XPS, ainsi que le suivi de l’évolution du point isoélectrique des solides avec l’augmentation
de la concentration en cation en solution (zétamétrie), ont permis de proposer les différents
mécanismes entrant en jeu (expériences réalisées à 25°C). Dans le cas des ferrites, les courbes
de sorption peuvent être modélisées avec un seul complexe de surface ≡ XOM + (M
représentant Co ou Ni). La sorption est le seul mécanisme intervenant entre les espèces
solubles et particulaires, pour des concentrations en ions cobalt ou nickel en solution,
inférieures à 10-5 M. Pour des concentrations plus élevées et à partir de pH 8, un mécanisme
de précipitation de surface est observé. Les précipités formés sont des hydroxydes : Co(OH ) 2
et Ni (OH ) 2 . Dans le cas de la magnétite, la complexation de ≡ XOM + et la précipitation de
M (OH )2 sont accompagnées de la formation d’un second complexe de surface : (≡ XO)2 M .
La modélisation des courbes de sorption a permis de déterminer les constantes de ces
différentes espèces à 25°C. Celles-ci sont résumées dans le tableau VII – 1.
Les expériences réalisées pour une température plus élevée (50°C) ont indiqué que la sorption
des cations a lieu pour des pH plus acides qu’à température ambiante, ce qui est en accord
avec l’évolution des PZC des oxydes avec la hausse de la température.
En comparant les résultats obtenus sur les systèmes ferrites, il apparaît que ceux-ci ont un
comportement similaire à la fois à 25°C et à 50°C. Leur comparaison avec les systèmes de la
magnétite indique une différence majeure : deux complexes de surface ( ≡ XOM + et
(≡ XO)2 M ) semblent apparaître lors de la mise en contact de cet oxyde avec les cations,
tandis qu’un seul ( ≡ XOM + ) se forme dans le cas des ferrites.
147
Tableau VII - 1 : Logarithmes des constantes moyennes de complexation et de solubilité,
déterminées à l’aide du logiciel ECOSAT (2-pK/BSM) pour la sorption du cobalt et du
nickel sur les trois produits de corrosion étudiés (25°C)
Oxydes
Concentrations CoFe2O4 NiFe2O4 Fe3O4
Initiales (M)
Ni : -5,45 ± 0,45
logK≡XOM+ -2,25 ± 0,05 -2,5 ± 0,5
Co : -5,25 ± 0,25
Ni : -11 ± 0,5
logK(≡XO)2M / /
Co : -10,5 ± 0,5
Ni : -15,45 ± 0,25
logKM(OH)2 -18 -17,5
Co : -15
Les mécanismes d’interaction entre les produits de corrosion et les cations en solution étant à
présent connus (au moins à basse température), l’étude du dépôt des particules sur les parois
du circuit primaire a été effectuée. La première partie de ce travail a été la détermination des
points isoélectriques des parois suivant deux méthodes distinctes : la zétamétrie et la mesure
de la vitesse d’adhésion de particules de latex chargées. Il a été démontré que les conditions
du milieu primaire (360°C, 1000 ppm de bore, 2 ppm de lithine, 30 [Link]-1 d’hydrogène,
80 jours) ont une influence très faible sur la valeur du PIE de l’Inconel 690, du Zircaloy 4 et
de l’acier inoxydable 304L. Dans les cas des trois oxydes, les valeurs sont d’environ 4.
A partir des divers résultats obtenus dans ce travail, il est alors possible de tenter une
prédiction du dépôt des produits de corrosion, présentant une sorption de surface ou non, sur
les parois du circuit. L’utilisation d’un code de calcul EDF, Code_Saturne, incluant un
module de dépôt de particules, a permis de prendre en compte à la fois les paramètres
physico-chimiques de la surface et des produits de corrosion et l’hydraulique du fluide
primaire. Bien que ce modèle ne soit pas entièrement validé (en ce qui concerne le module
Lagrangien), les calculs réalisés permettent une première compréhension du phénomène de
dépôt.
Plusieurs grandes tendances ont pu être obtenues à partir du couplage de ces données à l’aide
du module Lagrangien :
- lorsque les signes des potentiels de la surface et des particules sont opposés, le dépôt a
toujours lieu.
- lorsque ceux-ci sont identiques, le dépôt a lieu jusqu’à une valeur de potentiel
électrostatique qui est d’autant plus grande (en valeur absolue) que le diamètre des particules
est petit. Le dépôt en conditions électrostatiques répulsives est dû au mouvement Brownien.
Il est donc apparu que l’hydraulique joue un rôle important dans le dépôt, et que la simple
prédiction à partir de la comparaison des charges de surface n’est pas satisfaisante.
148
Pour les trois matériaux représentatifs du circuit primaire (Zircaloy 4, Inconel 690 et acier
inoxydable 304L), il semble que l’adhésion des produits de corrosion est très largement
présente entre pH 3,7 et pH ≈10, quelque soit le PC et la sorption de cobalt ou de nickel qu’il
a subi précédemment. La chimie du fluide primaire étant très peu modifiable (en
fonctionnement le pH à 300°C est de 7,2), il semble donc qu’à ce pH, le dépôt des trois
produits de corrosion étudiés, présentant éventuellement une sorption de cations à leur
surface, soit possible, et qu’il soit donc très difficile de jouer sur la chimie pour limiter ce
phénomène en phase de fonctionnement.
Les expériences ayant été faites à 25°C (et 50°C en ce qui concerne la sorption), il semble par
la suite important de les réaliser pour des températures plus élevées, et ce jusqu’à 320°C, dans
le cas de la sorption du cobalt et du nickel. Toutefois, les comportements des systèmes
Co 2+ / NiFe2O4 et Ni 2+ / CoFe2O4 étant similaires jusqu’à 50°C, le suivi en température de
l’un des deux systèmes devrait être suffisant. Pour les deux systèmes mettant en jeu la
magnétite, cette remarque est également valable. La réalisation de ces expériences en
température est nécessaire dans le but de se rapprocher le plus possible des conditions réelles
de fonctionnement. Il semble également indispensable par la suite, de déterminer la
compétition pouvant avoir lieu entre les deux cations vis-à-vis d’une même surface, puisque
dans les faits ceux-ci sont présents dans le fluide primaire simultanément.
Les logiciels permettant la modélisation des données de sorption ne permettent pas, à notre
connaissance, de prendre en compte la température. L’amélioration de ces codes est donc un
élément essentiel pour la détermination des constantes de complexation et de solubilité des
différentes espèces formées.
Le module Lagrangien de Code_Saturne doit également être amélioré et validé pour que les
résultats obtenus pour les calculs soient le plus prédictifs possible. En particulier la surface du
cylindre utilisée pour les calculs, qui est considérée plane. Hors dans les centrales, des
aspérités sont présentes sur les parties métalliques du circuit. Le modèle nécessite donc
l’intégration de paramètres de rugosité au niveau de la surface prise en compte. Notons que ce
point est actuellement en cours d’étude à EDF.
Le maillage utilisé correspond aux caractéristiques d’un tube de générateur de vapeur
uniquement. De ce fait, l’augmentation des dimensions du cylindre ainsi que le nombre de
Reynolds extrêmement élevé dans le cas d’une tuyauterie primaire (autre qu’un tube GV) ne
permettent pas de réaliser les calculs souhaités. Des améliorations dans ces cas également
doivent être faites.
149
Concernant l’apport de ce travail à la maîtrise de la contamination, nous avons constaté que la
complexation de surface serait très largement favorisée par rapport à la précipitation de
surface dans les conditions de fonctionnement, du fait des rapports surface/solution en jeu.
Par ailleurs, du fait des fortes contraintes sur la chimie du fluide primaire, de sa plage de
variation réduite, ainsi que des phénomènes observés qui sont stables dans la plage concernée,
il semble difficile de limiter la sorption ou l’adhésion par le levier chimique.
Le comportement des produits de corrosion dans le fluide primaire, tant entre eux que vis-à-
vis des parois métalliques du circuit résulte de la combinaison de plusieurs phénomènes
complexes. Les travaux réalisés dans cette thèse ont permis de mieux comprendre les
différents mécanismes entrant en jeu. Ceux-ci étant maintenant définis de manière isolée il
semble indispensable de les étudier simultanément et à plus haute température, afin de se
rapprocher autant que possible des conditions physico-chimiques réelles des centrales REP.
Les données acquises pourront ensuite alimenter les codes dédiés à la modélisation de la
contamination du circuit primaire comme le code Oscar développé par le CEA, ainsi que la
base de données EQ3/6 du code de calcul JChess pour la prise en compte de la complexation
de surface.
150
ANNEXES
La structure spinelle tire son nom du minéral MgAl2O4 qui cristallise dans le système
cubique. Les anions oxygène forment un réseau cubique à faces centrées. Dans ce type
d’empilement, deux types de sites interstitiels existent, tétraédriques et octaédriques entourés
respectivement par 4 et 6 atomes d’oxygène. Dans une maille élémentaire, il existe 64 sites
tétraédriques (T) et 32 octaédriques (O) dont seulement 8 et 16 sont respectivement occupés
comme présenté dans la figure 1-1.
Si les 8 ions divalents se trouvent dans les sites tétraédriques et les 16 ions trivalents dans les
sites octaédriques, le spinelle est dit normal. C’est le cas du minéral MgAl2O4.
Cependant, il existe des spinelles dans lesquels les 8 ions divalents occupent 8 des 16 sites
octaédriques et les 16 ions trivalents se répartissent entre sites octaédriques et sites
tétraédriques. Le spinelle est dit inverse.
151
Annexe 2 : Données de pureté des différents produits chimiques utilisés
Poudres :
CoFe2O4 : Pureté = 98 %
Lot : 1510FY
Nanostructured and Amorphous materials INC
NiFe2O4 : Pureté = 98 %
Lot : 5058
Nanostructured and Amorphous materials INC
Fe3O4 : Pureté = 99,997 %
Lot : 23868
Nanostructured and Amorphous materials INC
Autres produits :
152
Annexe 3 : Lavage du ferrite de cobalt (CoFe2O4 ), du ferrite de nickel
(NiFe2O4) et de la magnétite (Fe3O4).
Le protocole de lavage des poudres a été mis au point et optimisé lors de deux études (Barale-
2006, Mansour-2008).
- environ 15 g de poudre sont mélangés dans 500 mL de liquide, de pH donné, dans un flacon
en polyéthylène haute densité (PEHD) ;
- l’ensemble est agité pendant 10 minutes à température ambiante ;
- la suspension est décantée magnétiquement pendant 10 minutes : le ferrite de nickel, le
ferrite de cobalt et la magnétite se magnétisent en présence d’un champ magnétique ;
- le surnageant est aspiré puis son pH est mesuré et il est analysé par ICP-AES ;
Ces opérations sont répétées tant que le pH du surnageant évolue. Une fois que le lavage est
jugé suffisamment efficace, la poudre est filtrée sur filtre Millipore® (acétate de cellulose -
0,45 µm). Elle est ensuite séchée à l’étuve entre 50°C et 100°C puis stockée dans un
dessicateur contenant du Silica gel.
Le ferrite de cobalt est lavé deux fois avec une solution de KOH à 0,1 M, puis une fois avec
du HNO3 0,1 M et enfin abondamment à l’eau Millipore®. Les analyses ICP-AES montrent
que du cobalt et du fer passent en solution laissant penser à une dissolution du ferrite de
cobalt. Aucune différence significative n’a été détectée, avant et après lavage, sur le rapport
Co / Fe . La concentration en solution de chaque élément reste constante ou diminue au cours
des trois lavages. Les principales impuretés sont le sodium, le calcium et le magnésium.
Le ferrite de nickel est lavé avec une solution contenant 10-4 M d’acide nitrique (HNO3), puis
abondamment avec de l’eau Millipore®. Entre le premier et le dernier lavage avec une
solution contenant 10-4 M de HNO3, les concentrations en nickel, calcium et magnésium en
solution augmentent. Celles des autres éléments diminuent, c’est le cas du fer, où elle reste
constante. Ces analyses laissent penser qu’une petite partie du nickel du ferrite se dissout.
Comme pour le ferrite de cobalt celles-ci ne montrent pas de différence significative du
rapport Ni / Fe . Les autres impuretés n’ont pas fait l’objet d’une quantification.
La magnétite est lavée deux fois avec une solution de KOH à 10-2 M puis abondamment avec
de l’eau Millipore®. Les impuretés relâchées en solution sont le sodium, le calcium,
l’aluminium et le barium.
153
Annexe 4 : Détermination des limites de détection et de quantification lors
des analyses ICP-AES et ICP-MS.
Les limites de détection et de quantification des éléments (cobalt et nickel) ont été mesurées
et sont données dans le tableau 4 - 1.
154
Annexe 5 : Courbes de détermination du point de charge nulle par titrage
de masse.
La détermination du point de charge nulle des oxydes par titrage de masse a été réalisée lors
d’une étude précédente (Barale-2006). Les courbes sont données dans cette annexe (figure 5 –
1, 5 – 2 et 5 – 3). Pour réaliser ces expériences, KCl à 10-2 M a été utilisé.
Les trois courbes se rejoignent pour un pH d’environ 6,4 à partir d’une concentration en
oxyde de 50 g.L-1. Le point de charge nulle du ferrite de cobalt est donc 6,4 ± 0,1.
155
Dans le cas du ferrite de nickel le PZC déterminé est 8,1 ± 0,1. La stabilisation des courbes
est atteinte pour une concentration en oxyde de 50 g.L-1.
Pour la magnétite, la valeur du point de charge nulle est de 6,5 ± 0,1 et est atteinte pour une
concentration de 75 g.L-1.
156
Annexe 6 : Potentiel zêta du Zircalloy 4 vieilli et de l’acier inoxydable 304L
non vieilli déterminé par zétamétrie.
La détermination du point isoélectrique de l’Inconel 690 (vieilli et non vieilli), a été détaillée
dans le chapitre VI. Dans cette annexe sont données les courbes, obtenues par la même
méthode et dans les mêmes conditions expérimentales, du Zircaloy 4 vieilli et de l’acier
inoxydable 304L non vieilli.
40
20
Potentiel zéta (mV)
-20
-40
-60
-80
2 3 4 5 6 7 8
pH
Figure 6 - 1 : Evolution du potentiel zêta du Zircaloy 4 vieilli en autoclave en fonction du
pH (10-3 M en NaNO3)
10
0
Potnetiel zéta (mV)
-5
-10
-15
-20
-25
3 4 5 6 7 8 9 10 11
pH
157
Annexe 7 : Observations MEB des échantillons massifs après adhésion des
particules de latex.
a b
1 µm
1 µm
c d
1 µm
1 µm
Figure 7 - 1 : Observations MEB de l’Inconel 690 après mise en contact avec les
particules de latex (Amidine) ( a-b : vieilli, c-d : non vieilli, a-c : pH 3,3, b-d : pH 5,5)
158
a b
1 µm 500 nm
c d
500 nm 1 µm
Figure 7 - 2 : Observations MEB du Zircaloy 4 vieilli après mise en contact avec les
particules de latex ( a-b : pH 3,3, c-d : pH 5,5, a-c : Carboxyle, b-d : Amidine)
159
a b
1 µm 1 µm
c d
1 µm 1 µm
Figure 7 - 3 : Observations MEB de l’acier inoxydable 304L vieilli après mise en contact
avec les particules de latex ( a-b : pH 5,5, c-d : pH 3,3, a-c : Amidine, b-d : Carboxyle)
160
a b
1 µm 1 µm
c d
1 µm 1 µm
Figure 7 - 4 : Observations MEB de l’acier inoxydable 304L non vieilli après mise en
contact avec les particules de latex ( a-b : pH 5,5, c-d : pH 3,3, a-c : Amidine, b-d :
Carboxyle)
161
Annexe 8 : Superposition des courbes de zétamétrie du Zircaloy 4 et de
l’acier inoxydable 304 L avec celles des produits de corrosion.
32.5
Potentiel zéta (mV)
-32.5
sans NiCl
2
-4
-65 10 M
Zircaloy 4 vieilli
2 4 6 8 10 12
pH
30
Potentiel zéta (mV)
-30
sans CoCl2
-60
10-4 M
Zircaloy 4 vieilli
2 4 6 8 10 12
pH
Figure 8 - 2 : Evolution du potentiel zêta et PIE du Zircaloy 4 après vieillissement
([NaNO3] = 10-3 M) en autoclave et du ferrite de nickel ([KCl] = 10-4 M) en présence de
chlorure de cobalt en fonction du pH ([NiFe2O4] = 0,4 g.L-1, [CoCl2] = 0 M et 10-4 M,
25°C)
162
32.5
-32.5
-4
Ni 10 M
-65 -4
Co 10 M
Zircaloy 4 vieilli
2 4 6 8 10 12
pH
Figure 8 - 3 : Evolution du potentiel zêta et PIE du Zircaloy 4 après vieillissement
([NaNO3] = 10-3 M) en autoclave et de la magnétite ([KCl] = 10-4 M) en présence de
chlorure de nickel ou de cobalt en fonction du pH ([Fe3O4] = 0,4 g.L-1, [NiCl2] et
[CoCl2] = 0 M et 10-4 M, 25°C)
50
25
Potentiel zéta (mV)
163
sans CoCl
2
35 -4
10 M
acier inoxydable
304L non vieilli
-17.5
2 4 6 8 10 12
pH
50
25
Potentiel zéta (mV)
164
REFERENCES
Anthoni-1994
S. Anthoni, P. Ridoux, C. Caramel, International Conference of Chemistry in Water
Reactors, Nice, 1994
Bachet-2008
M. Bachet, Document EDF H-T29-2007-02237-FR, 2008
Baes-1986
C.F. Baes, R.E. Mesmer, « The Hydrolysis of Cations », Krieger Publishing
Compagny, Malabar, 1986
Barale-2006
M. Barale, Thèse de doctorat, Université Paris VI, 2006
Barale-2008
M. Barale, C. Mansour, F. Carette, E.M. Pavageau, H. Catalette, G. Lefèvre, M.
Fedoroff, G. Cote, J. Nucl. Materials, 381, 2008, 302
Barron-1996
V. Barron, J. Torrent, J. Colloid Interface Sci., 177, 1996, 407
Bauer-2005
A. Bauer, T. rabung, F. Claret, T. Schäfer, G. Buckau, T. Fanghänel, Applied Clay
Sci., 30, 2005, 1
Baumgarten-1997
E. Baumgarten, U. Kirchhausen-Düsing, J. Colloid Interface Sci., 194, 1997, 1
Berger-2002
M. Berger, Document EDF ENFIN/0200210A, 2002
Bérubé-1968
Y.G. Bérubé, P.L. De Bruyn, J. Colloid Interface Sci., 27, 1968, 305
Beukes-2000
J.P. Beukes, E.W. Giesekke, W. Elliott, Minerals Engineering, 13, 2000, 1573
Blakey-2003
B.C. Blakey, D.F. James, Colloids and Surfaces A : Physicochem. Eng. Aspects, 231,
2003, 19
Blesa-1984
M.A. Blesa, N. M. Figliolia, A.J.G. Maroto, A.E. Regazonni, J. Colloid Interface, Sci.,
101, 1984, 410
Bolt-1982
G.H Bolt, W.H. Van Riemsdijk, « Soil chemistry B. Physical-Chemical models 2nd
edition », G.H Bolt (Ed), Elsevier Amsterdam, 1982, 459-505
Boudot-1992
R. Boudot, Doc. EDF HT-40/NEQ1273-B, 1992
165
Boulangé-Petermann-1995
L. Boulangé-Petermann, A. Doren, B. Baroux, M.N. Bellon-Fontaine, J. Colloid
Interface Sci., 171, 1995, 179
Bourg-1983
A. Bourg, Documents du BRGM, n°63, 1983
Brun-1994
C. Brun, A. Long, M.C. Thiry, N. Lacoudre, International Conference of Chemistry in
Water Reactors, Nice, 1994
Brunauer-1938
S. Brunauer, P.H. Emmett, E. Teller, J. Am. Chem. Soc., 60, 1938, 309
Carette-2008
F. Carette, L. Guinard, B. Pieraggi, International Conference on Water Chemistry of
Nuclear Reactor Systems, VGB, Berlin, 2008
Catalette-1998
H. Catalette, J. Dumonceau, P. Ollar, J. Contaminant Hydrology, 35, 1998, 151
Catalette-2004
H. Catalette, Thèse de doctorat, Université Paris XI, 2004
Cattant-2008
F. Cattant, Document EDF, H-T20-2007-03096-FR, 2008
Cerovic-2009
L. Cerovic, G. Lefèvre, A. Jaubertie, M. Fédoroff, S. Milonjic, J. Colloids Interface
Sci., 330, 2009, 284
Chen-2009
C. Chen, J.H. Dadong Shao, J. Li, X. Wang, J. Hazardous Materials, 164, 2009, 923
Christl-1999
I. Christl, R. Kretzschmar, Geochim. Cosmochim. Acta, 63, 1999, 2929
Conway-1977
B.E. Conway, Adv. Colloid Interface Sci., 8, 1977, 91
Cordier-2007
H. Cordier, A. Douce, B. Gaudron, Document EDF, H-I81-2006-03163-FR, 2007
Coriou-1962
H. Coriou, L. Grall, J. Meunier, M. Pelras, H. Willermoz, J. Nucl. Materials, 7, 1962,
320
Crawford-1996
R.J. Crawford, I.H. Harding, D.E. Mainwaring, J. Colloid Interface Sci., 181, 1996,
561
Cremel-2007
S. Cremel, Thèse de doctorat, Université Nancy-I, 2007
Dacquait-2002
F. Dacquait, C. Andrieu, M. Berger, J.L. Bretelle, A. Rocher, International Conference
– Water Chemistry in Nuclear Reactors Systems, Avignon, 2002
166
Demmitt-1972
T.F. Demmitt, A.B. Johnson, International Corrosion Conference, Mars 1972, Saint
Louis
De-Regge-1988
P. De Regge, K. Dinov, K. De Ranter, Japan Atomic Industrial Forum, Avril 1988,
Tokyo
Dumont-Ravat-1998
C. Dumont-Ravat, Thèse de doctorat, Université Reims Champagne-Ardenne, 1998
Dyrssen-1980
D. Dyrssen, M. Wedborg, « Chemistry and biochemistry of Estuaries », E. Olausson
and I. Cato (EDs), Wiley ans sons, 1980, 71-119
Dzombak-1990
D.A. Dzombak, F.M.M. Morel, « Surface complexation modeling :Hydrous ferric
oxide », John Wiley and sons , New York, 1990
Farley-1985
K.J. Farley, D.A. Dzombak, F.M.M. Morel, J. Colloid Interface Sci., 106, 1985, 226-
246
Feitknecht-1963
W. Feitknecht, P. Schindler, « Pure and applied chemistry », Butterworths, London,
1963
Fokkink-1989
L.G.J. Fokkink, A. De Keizer, J. Lyklema, J. Colloid Interface Sci., 127, 1989, 116
Galtayries-1999
A. Galtayries, J. Grimblot, J. Elec. Spec. Rel. Phen., 98-99, 1999, 267
Gimenez-2007
J. Gimenez, M. Martinez, J. de Pablo, M. Rovira, L. Duro, J. Hazardous Materials,
141, 2007, 575
Goldberg-1985
S. Goldberg, R.A. Glaubig, Soil Sci. Soc. Am. J., 49, 1985, 1374
Goldberg-1993
S. Goldberg, H.S. Forster, E.L. Heick, Soil Sci., 156, 1993, 316
Grahame-1947
D.C. Grahame, Chem. Rev., 41, 1947, 441
Grandgeorge-1998
S. Grandgeorge, C. Jallut, B. Thonon, Chem. Eng. Sci., 53, 1998, 3051
Gras-1987
J. M. Gras, Materials Science and Engineering, 88, 1987, 47-52
Gressier-2008
F. Gressier, Thèse de doctorat, ENSMP, 2008
Grosvenor-2004
A.P. Grosvenor, B.A. Kobe, M.C. Biesinger, [Link].C. Smart, N.S. McIntyre, Surf.
Interface Anal., 36, 2004, 1564
167
Grosvenor-2006
A.P. Grosvenor, M.C. Biesinger, [Link].C. Smart, N.S. McIntyre, Surf. Sci., 600, 2006,
1771
Guillodo-2004
M. Guillodo, P. Combrade, B. dos Santos, T. Muller, G. Berthollon, N. Engler, C.
Brun, G. Turluer, International Conference of Water Chemistry of Nuclear Reactor
Systems, San Fransisco, 2004
Guinard-1995
L. Guinard, Document EDF HT-45/95/024/A, 1995
Guinard-2003
L. Guinard, Document EDF, ENTERP030024A, 2003
Guinard-2005
L. Guinard, Document EDF, ENTERP040302A, 2005
Guingo-2008
M. Guingo, J.P. Minier, Phys. Fluids, 20, 2008, 053303
Gunnarsson-2000
M. Gunnarsson, A.M. Jacobsson, S. Ekberg, Y. Albinsson, E. Ahlberg, J. Colloid
Interface Sci., 231, 2000, 326
Guzonas-2002
D. Guzonas, D. Webb, Rapport EPRI 1003384, 2002
Hawkes-2004
J.M. Hawkes, Thèse de doctorat, Georgia Institute of Technology, 2004
Henry-1931
D.C. Henry, Proc. Roy. Soc., London, A133, 1931
Herbelin-1999
A.L. Herbelin, J.C. Westall, FITEQL 4.0, Department of Chemistry, Oregon State,
University report 99-01, 1999
Hermansson-1977
H.P. Hermansson, Chem. Scripta, 12, 1977, 102
Hiemstra-1989a
T. Hiemstra, W.H. Van Riemsdijk, G.H. Bolt, J. Colloid Interface Sci., 133, 1989, 91
Hiemstra-1989b
T. Hiemstra, J.C.M. De Wit, W.H. Van Riemsdijk, J. Colloid Interface Sci., 133, 1989,
105
Hiemstra-1991
T. Hiemstra, W.H. Van Riemsdijk, Colloids and Surfaces, 59, 1991, 7
Hohl-1976
H. Hohl, W. Stumm, J. Colloid Interface Sci., 55, 1976, 281
Hückel-1924
E. Hückel, Physik Z., 25, 1924, 204
Hunter-1981
R.J. Hunter, « Zeta potential in colloid science – Principles and applications »,
Academic Press, London, 1981
168
Hunter-2001
R.J. Hunter, « Foundations of colloid science », Oxford University Press, New York,
2001
Jayaweera-1993
P. Jayaweera, S. Hettiarachchi, Rev. Sci. Instrum., 64, 1993, 524
Jayaweera-1994
P. Jayaweera, S. Hettiarachchi, H. Ocken, Colloids Surfaces A : Physycochem. Eng.
Aspects, 85, 1994, 19
Jolivet-1994
J.P. Jolivet, « De la solution à l’oxyde », CNRS Edition : EDP Science, Paris, 1994
Jones-2008
B.G. Jones, H. Bindra, Q. Rao, Rapport EPRI 1016184, 2008
Jurkiewicz-1986
K. Jurkiewicz, Int. J. Mineral Processing, 17, 1986, 67
Kallay-1987
N. Kallay, Z. Torbic, E. Barouch, J. Jednacak-Biscan, J. Colloid Interface Sci., 118,
1987, 431
Kallay-1991
N. Kallay, Z. Torbic, M. Golic, E. Matijevic, J. Physical Chem., 95, 1991, 7028
Kanbe-1983
H. Kanbe, T. Inoue, T. Tomizawa, H. Koyama, Nucl. Technol., 60, 1983, 367
Keizer-2009
M.G. Keizer, W.H. Van Riemsdijk, ECOSAT 4.9, department of soil quality,
Wageningen University, Netherlands,
Kim-1998
K-R. Kim, U-C. Kim, J. Ind. Eng. Chem., 4, 1998, 270
Kleimann-2006
J. Kleimann, G. Lecoutre, G. Papastavrou, S. Jeanneret, P. Galletto, G.J.M. Koper, M.
Borkovec, J. Colloid Interface Sci., 303, 2006, 460
Kosmulski-1997
M. Kosmulski, J. Colloid Interface Sci., 190, 1997, 212
Kosmulski-2001
M. Kosmulski, « Chemical properties of material surfaces », Marcel Dekker Inc., New
York, 2001
Kosmulski-2009
M. Kosmulski, « Surface charging and points of zero charge », CRC Press, Boca
Raton, 2009
Kratgen-1978
J. Kratgen, « Atlas of Metal-Ligand equilibria in aqueous solution », John Wiley &
Sons, New York, 1978
Kritzer-1999
P. Kritzer, N. Boukis, E. Dinjus, J. Supercritical Fluids, 15, 1999, 205
169
Kuo-1980
R.J. Kuo, E. Matijevic, J. Colloid Interface Sci., 78, 1980, 407
Laghoutaris-2009
P. Laghoutaris, Thèse de doctorat, ENSMP, 2009
Lefèvre-2006
G. Lefèvre, A. Hamza, M. Fédoroff, F. Carrette, H. Cordier, Colloids and Surfaces A :
Physicochem. Eng. Aspects, 280, 2006, 32
Lefèvre-2009
G. Lefèvre, L. Cerovic, S. Milonjic, M. Fédoroff, J. Finne, A. Jaubertie, J. Colloid
Interface Sci., 337, 2009, 449
Lemarchand-2005
E. Lemarchand, J. Schott, J. Gaillardet, Geochim. Cosmochim. Acta, 69, 2005, 3533
Lister-1975a
D.H. Lister, Nucl. Sci. Eng., 59, 1976, 406
Lister-1975b
D.H. Lister, Nucl. Sci. Eng., 58, 1975, 239
Litton-1996
G.M. Litton, T.M. Olson, Coll. Surf. A : Physicochem. Eng. Aspects, 107, 1996, 273
Lvov-2003
S.N. Lvov, X.Y. Zhou, G.C. Ulmer, H.L. Barnes, D.D. ManDonald, S.M. Ulyanov,
L.G. Benning, D.E. Granstaff, M. Manna, E. Vicenzi, Chem. Geol., 198, 2003, 141
Ly-2009
A. Ly, Thèse de doctorat, INPG, 2009
Machesky-1998
M.L. Machesky, D.J. Wesolowski, D.A. Palmer, K. Ichiro-Hayashi, J. Colloid
Interface Sci., 200, 1998, 298
Mansour-2007
C. Mansour, Thèse de doctorat, Université Paris VI, 2007
Mansour-2009
C. Mansour, G. Lefèvre, E.M. Pavageau, H. Catalette, M. Fedoroff, S. Zanna, J.
Colloid Interface Sci., 331, 2009, 77
Mansour-2010
C. Mansour, M. Fedoroff, G. Lefèvre, A. Pages, E.M. Pavageau, H. Catalette, , S.
Zanna, J. Colloid Interface Sci., in press, 2010
Marmier-1994
N. Marmier, Thèse de doctorat, Université de Reims Champagne-Ardenne, 1994
Marmier-1999
N. Marmier, A. Delisée, F. Fromage, J. Colloid Interface Sci., 211, 1999, 54
Marthur-1998
B.S. Marthur, B. Venkataramani, Colloids Surfaces A, 140, 1998, 403
McIntyre-1975
N.S. McIntyre, M.G. Cook, Anal. Chem., 47, 1975, 2208
170
Misak-1996
N.Z. Misak, H.F. Ghoneimy, T.N. Morcos, J. Colloid Interface Sci., 184, 1996, 31
Missana-2003
T. Missana, M. Garcia-Gutierrez, V. Fernandez, Geochim. Cosmochim. Acta, 67,
2003, 2543
Motl-2003
A. Motl, F. Sebesta, Czechoslovak Journal of Physics, 53, Suppl. A, 2003, A515
Naeem-2009
A. Naeem, M.T. Saddique, S. Mustafa, S. Tasleem, K.H. Shah, W. Waseem, J.
Hazardous Materials, 172, 2009, 124
Nordstrom-1979
D.K. Nordstrom, L.M. Plummer, T.L.M. Wigley, T.J. Wolery, J.W. Ball, E.A. Jenne,
R.L. Bassett, D.A. Crerar, T.M. Florence, B. Fritz, M. Haffman, G.R. Hochen, G.M.
Lafon, S.V. Mattigod, R.E. McDuff, F. Morel, M.M. Reddy, G. Sposito, J. Thrailkall,
« Chemical modeling in aqueous systems », E.A. Jenne (Ed), ASC Symposium series
n°93, ACS, Washington, 1979, 857-892
Ohshima-1994
H. Ohshima, J. Colloid Interface Sci., 168, 1994, 269
Osinska-2004
M. Osinska, T. Stefanowicz, D. Paukzsta, J. Hazardous Materials A, 112, 2004, 177
Overbeek-1952
[Link].G. Overbeek, « Colloid science », H.R. Kruyt Ed., Vol. 1, Elsevier, Amsterdam,
1952
Papavergos-1984
P.G. Papavergos, A.B. Hedley, Chem. Eng. Res. Des., 62, 1984, 275
Parks-1965
G.A. Parks, Chem. Rev., 65, 1965, 177
Parsi-2009
A. Parsi, A. Byers, R.L. DeVito, Rapport EPRI 1018720, 2009
Peak-2003
D. Peak, G.W. Luther, D.L. Sparks, Geochim. Cosmochim. Acta, 67, 2003, 2551
Persson-1995
P. Persson, G.A. Parks, G.E. Brown Jr., Physica B, 208 & 209, 1995, 453
Philippini-2006
V. Philippini, A. Naveau, H. Catalette, S. Leclercq, J. Nucl. Materials, 348, 2006, 60
Pivovarov-1997
S. Pivovarov, J. Colloid Interface Sci., 196, 1997, 321
Plyasunova-1998a
N.V. Plyasunova, Y. Zhang, M. Muhammed, Hydrometallurgy, 48, 1998, 43
Plyasunova-1998b
N.V. Plyasunova, Y. Zhang, M. Muhammed, Hydrometallurgy, 48, 1998, 153
171
Pourbaix-1963
M. Pourbaix, « Atlas d’équilibres électrochimiques », Gauthier-Villars & Cie Ed.,
Paris, 1963
Rahnemaie-2006
R. Rahnemaie, T. Hiemstra, W.H. Van Riemsdijk, J. Colloid Interface Sci., 293, 2006,
312
Ranchoux-2008
G. Ranchoux, L. Guinard, Document EDF ENTERP080167, 2008
Regazzoni-1983
A.E. Regazzoni, M.A. Blesa, A.J.G. Maroto, J. Colloid Interface Sci., 91, 1983, 560
Ridley-1999
M.K. Ridley, M.L. Machesky, D.J. Wesolowski, D.A. Palmer, Geochim. Cosmochim.
Acta, 63, 1999, 3087
Ryde-2000
N.P. Ryde, E. Matijevic, Colloid Surf. A : Physicochem. Eng. Aspects, 165, 2000, 59
Saint Paul-1993
P. Saint Paul, Journal de Physique IV, Colloque 9, supplément au Journal de Physique
III, 3, 1993
Sanyal-2000
A. Sanyal, M. Nugent, R.J. Reeder, J. Bijma, Geochim. Cosmochim. Acta, 64, 2000,
1551
Saturne-2007
[Link]
Sawicki-2008
J.A. Sawicki, J. Nucl. Materials, 374, 2008, 248
Schindler-1972
P.W. Schindler, H. Gamsjäger, Kolloid Zeitschrift u. Zeitschrift Polymere, 250, 1972,
759
Schindler-1987
P.W. Schindler, W. Stumm, « Aquatic surface chemistry », Stumm W. (Ed), John
Wiley and sons, New York, 1987, 83-110
Schmoluchowski-1921
M. von Schmoluchowski, « Handbuch der Electrizität und des Magnetismus »,
(Graetz), Barth, Leipzig, vol. II, 1921
Schoonen-1994
M.A.A. Schoonen, Geochim. Cosmochim. Acta, 58, 1994, 2845
Song-2009
X. Song, S. Wang, L. Chen, M. Zhang, Y. Dong, Applied Radiation Isotopes, 67,
2009, 1007
Stern-1924
O. Stern, Z. Elektrochem., 30, 1924, 508
Strasser-1996
A.A. Strasser, H.J. Schroeder, K.D. Sheppard, Rapport EPRI TR-107255, 1996
172
Stumm-1970
W. Stumm, C.P. Huang, S.R. Jenkins, Croat. Chem. Acta, 42, 1970, 223-245
Stumm-1987
W. Stumm, « Aquatic surface chemistry », John Wiley and sons, New York, 1997
Stumm-1991
W. Stumm, « Chemistry of the solid-water interface – Processes at the minéral-water
and particles-water interface in natural systems », Wiley – Interscience publication,
John Wiley and sons, New York, 1991
Stutzmann-1997
A. Stutzmann, Document EDF D.5001/NRN/R84442-ind.2, 1997
Tamura-1983
H. Tamura, E. Matijevic, L. Meites, J. Colloid Interface Sci., 92, 1983, 303
Tamura-1997a
H. Tamura, N. Katayama, R. Furuichi, J. Colloid Interface Sci., 195, 1997, 192
Tamura-1997b
H. Tamura, R. Furuichi, J. Cooloid Interface Sci., 195, 1997, 241
Taunier-2008
S. Taunier, A. Rocher, Document EDF EDEECH080138-A, 2008
Taunier-2009
S. Taunier, Document EDF EDEECH090242-A
Tertre-2005
E. Tertre, G. Berger, S. Castet, M. Loubet, E. Giffaut, Geochim. Cosmochim. Acta,
21, 2005, 4937
Tertre-2006
E. Tertre, G. Berger, E. Simoni, S. Castet, E. Giffaut, M. Loubet, H. Catalette,
Geochim. Cosmochim. Acta, 70, 2006, 4563
Tewari-1971
P.H. Tewari, A.W. Mclean, J. Colloid Interface Sci., 40, 1971, 267
Tewari-1972
P.H. Tewari, A.B. Campbell, Symposium on oxide-electrolyte interfaces, Miami
Beach, 1972
Tewari-1975
P.H. Tewari, W. Lee, J. Colloid Interface sci., 52, 1975, 77
Tewari-1976
P.H. Tewari, A.B. Campbell, J. Colloid Interface Sci., 55, 1976, 531
Towle-1995
S.N. Towle, J.A. Bargar, P. Persson, G.E. Brown Jr., G.A. Parks, Physica B, 208 &
209, 1995, 439
Tremaine-1980
P.R. Tremaine, J.C. LeBlanc, J. Solution Chem., 9, 1980, 415
Trémillon-1965
B. Trémillon, « Les séparations par les résines échangeuses d’ions », Gauthier-Villar,
1965
173
Uheida-2006
A. Uheida, G. Salazar-Alvarez, E. Björkman, Z. Yu, M. Muhammed, J. Colloid
Interface Sci., 298, 2006, 501
Van der Lee-2001
J. Van der Lee, JChess 2.0, Ecole des mines de Paris, Centre d’informatique
géologique, [Link] 2001
Venkataramani-1985
B. Venkataramani, A.R. Gupta, R.M. Iyer, J. Radioanal. Nucl. Chem., 96, 1985, 129
Westall-1980
J. Westall, H. Hohl, Adv. Colloid Interface Sci., 12, 1980, 265-294
Wesolowski-2000
D.J. WesolowskiM.L. Machesky, D.A. Palmer, L.M. Anovitz, Chem. Geol., 167,
2000, 193
Wolery-1992
T. Wolery, EQ3/6 (V.8-R.6) database, 1992
Zhu-2002
C. Zhu, Chem. Geol., 188, 2002, 23
174
LISTE DES FIGURES
175
FIGURE III - 7 : SCHEMA DE L’INTERFACE SOLIDE-LIQUIDE SELON LE MODELE DE GRAHAME
(Ψ : POTENTIEL, X : DISTANCE A LA SURFACE).................................................................. 39
FIGURE III - 8 : SCHEMA DU PLAN DE CISAILLEMENT A L’INTERFACE SOLIDE-LIQUIDE D’UNE
PARTICULE EN MOUVEMENT DANS UN FLUIDE ................................................................... 40
FIGURE III - 9 : SCHEMA DE L’INTERFACE SOLIDE-LIQUIDE DANS LE CAS D’UN HYDROXYDE. LE
PLAN 0 CORRESPOND A LA SURFACE DU SOLIDE, LE 1 AU PIH ET LE 2 AU PEH. LES
GRANDES SPHERES BLANCHES REPRESENTENT LES ATOMES D’OXYGENE ET LES PETITES LES
ATOMES D’HYDROGENE ..................................................................................................... 45
FIGURE III - 10 : ISOTHERME DE SORPTION CARACTERISTIQUE D’UNE PRECIPITATION DE
SURFACE (A : PHENOMENE D’ADSORPTION SEUL, B : PHENOMENE DE PRECIPITATION DE
SURFACE, [X] EST LA CONCENTRATION EN ELEMENT ........................................................ 46
FIGURE III - 11 : SPECIATION DU NICKEL ET DU COBALT EN SOLUTION A DIFFERENTES
2+ 2+ -6 -5 -4
CONCENTRATIONS EN NI OU CO (A-B : 10 M, C-D : 10 M, E-F : 10 M) A 25°C ..... 48
FIGURE III - 12 : A). MODELES DE COMPLEXES DE SURFACE (1 : 1) ET (1 : 2) DE CO2+ SORBES
SUR LA MAGNETITE (TAMURA-1997A, B). CONCENTRATION CALCULEE DES COMPLEXES DE
SURFACE (1 :1) ET (1 :2) EN FONCTION DU PH POUR UNE CONCENTRATION EN COBALT DE
10-7 M, EN SEL DE FOND (NANO3) DE 0,1 M A 25°C (TAMURA-1997A) ........................... 50
FIGURE III - 13 : A). POURCENTAGE DE CO2+ SORBE SUR LA MAGNETITE (M/V = 10 G/L) EN
FONCTION DU PH POUR DIFFERENTES CONCENTRATIONS EN COBALT A 25°C (TAMURA-
1983). B). INFLUENCE DE LA FORCE IONIQUE SUR LE POURCENTAGE DE CO2+ SORBE A 25°C.
(CONCENTRATION EN COBALT 1,7.10-7 M) (TAMURA-1983) C). POURCENTAGE DE CO2+
-7
SORBE EN FONCTION DU PH (25°C) POUR UNE CONCENTRATION EN COBALT DE 1,7.10 M
ET EN SEL DE FOND DE 0,1 M A DIFFERENTES TEMPERATURES (TAMURA-1983) ............... 51
FIGURE III - 14 : POURCENTAGE DE CO2+ SORBE SUR FE3O4 OU γ-FE2O3 EN FONCTION DU PH
(UHEIDA-2006) ([CO2+] = 8,4.10-5 M, M/V = 2,5 G/L, [NANO3] = 0,1 M, TEMPS DE
CONTACT = 15 MIN, TEMPERATURE = 22 ± 1°C) (UHEIDA-2006) ..................................... 52
FIGURE III - 15 : POURCENTAGE DE NI2+ SORBE SUR LA MAGNETITE EN FONCTION DU PH
(MARMIER-1999) ([NI2+] = 8,4.10–5 M, M/V = 8 G/L, [NANO3] = 0,1 M, TEMPS DE
CONTACT = 24 H, TEMPERATURE AMBIANTE).................................................................... 53
FIGURE III - 16 : EFFET DE LA FORCE IONIQUE SUR L’ADSORPTION DU NICKEL SUR LES
2+
NANOTUBES D’OXYDE DE FER EN FONCTION DU PH ([NI ] = 6 MG/L, M/V = 0,75 G/L, T =
25 ± 2°C (CHEN-2009) ..................................................................................................... 55
FIGURE III - 17 : ESTIMATION DU PZC DE LA MAGNETITE EN TEMPERATURE A PARTIR DE
DIVERSES ETUDES (WESOLOWSKI-2000 (CORRESPOND AUX COURBES NOTEES « THIS
STUDY »), BLESA-1984, SCHOONEN-1994)....................................................................... 56
FIGURE III - 18 : EVOLUTION DU POTENTIEL ZETA DE PDO, FE2O3, FE3O4 ET TIO2 A 235°C
(JAYAWEERA-1993, JAYAWEERA-1994) ........................................................................... 58
FIGURE III - 19 : FRACTION DE CA(II) ADSORBE SUR LE RUTILE EN FONCTION DE PHZNPC – PH A
DIFFERENTES TEMPERATURES – 25°C (), 50°C (), 100°C (), 150°C (), 200°C (),
250°C (). LES COURBES PLEINES CORRESPONDENT A LA SORPTION ET CELLES EN
POINTILLES A LA DESORPTION (RIDLEY-1999) .................................................................. 59
FIGURE III - 20 : ADSORPTION DE CO(II) EN FONCTION DE SA CONCENTRATION EN SOLUTION A
L’EQUILIBRE (PH DE REFERENCE A 25°C = 6,5) SUR NIFE2O4 (TEWARI-1975) ................. 60
176
FIGURE III - 21 : POURCENTAGE D’EU(III) SORBE SUR LA SMECTITE EN FONCTION DU PH A
25°C, 60°C ET 80°C (BAUER-2005) ................................................................................. 60
FIGURE III - 22 : LOG KD DE L’EUROPIUM SUR LA KAOLINITE A DIFFERENTES TEMPERATURES
(TERTRE-2006) ................................................................................................................. 61
FIGURE III - 23 : LOG KD DU NICKEL SUR LA MONTMORILLONITE A DIFFERENTES
TEMPERATURES (TERTRE-2005)........................................................................................ 61
FIGURE III - 24 : CONCENTRATION EN NI(II) SORBE SUR LA MONTMORILLONITE A DIFFERENTES
TEMPERATURES (PH = 5,8, I = 0,01 M, : 20°C, : 37°C, O : 67°C) (SONG-2009)..... 62
FIGURE III - 25 : SORPTION DES SULFATES SUR LA MAGNETITE A DIFFERENTES TEMPERATURES
(MANSOUR-2007) ............................................................................................................. 62
FIGURE IV - 1 : OBSERVATIONS AU MEB DU FERRITE DE COBALT (A : BARALE-2006) ............ 67
FIGURE IV - 2 : OBSERVATIONS AU MEB DU FERRITE DE NICKEL (A : BARALE-2006) ............. 67
FIGURE IV - 3 : OBSERVATIONS AU MEB DE LA MAGNETITE.................................................... 68
FIGURE IV - 4 : PRINCIPE DE FONCTIONNEMENT DU MICROSCOPE ELECTRONIQUE EN
TRANSMISSION................................................................................................................... 69
FIGURE IV - 5 : OBSERVATIONS AU MET DU FERRITE DE COBALT............................................ 70
FIGURE IV - 6 : EFFET DE LA TAILLE DE LA PARTICULE SUR LES LIGNES DE CHAMP ELECTRIQUE.
(A) : κA << 1, (B) : κA >> 1. LA LIGNE EN POINTILLEE EST A UNE DISTANCE 1/κ DE LA
SURFACE DE LA PARTICULE (HUNTER-2001) ..................................................................... 77
FIGURE IV - 7 : PRINCIPE DE LA PHOTOEMISSION ...................................................................... 82
FIGURE IV - 8 : APPAREILLAGE XPS UTILISE LORS DE L’ETUDE ............................................... 83
FIGURE IV - 9 : ZETAMETRE UTILISE POUR EFFECTUER LES MESURES DE POTENTIEL
D’ECOULEMENT ................................................................................................................. 84
FIGURE IV - 10 : ORIGINE DE LA DIFFERENCE DE POTENTIEL .................................................... 85
FIGURE IV - 11 : SCHEMA DE LA MISE EN PLACE DES ECHANTILLONS DANS LE ZETAMETRE ..... 86
FIGURE V - 1 : EVOLUTION DE LA QUANTITE DE NICKEL SORBEE SUR LE FERRITE DE COBALT EN
FONCTION DU TEMPS (25°C, [KCL] = 10 M, [NICL2] = 10 M, [COFE2O4] = 0,4 G.L , PH
-4 -6 -1
= 6) ................................................................................................................................... 87
FIGURE V - 2 : VARIATION DE LA CONCENTRATION EN NICKEL SORBE SUR LE FERRITE DE
COBALT EN FONCTION DU PH (25°C, [KCL] = 10 M, [NICL2] = 10 M, 5.10 M, 10 M
-4 -4 -5 -5
[COFE2O4] = 0,4 G.L-1, 8 H). LA CONCENTRATION SORBEE EST EXPRIMEE EN MOLES PAR
LITRE DE SOLUTION ........................................................................................................... 89
FIGURE V - 4 : SPECTRES XPS DU NICKEL AVANT ET APRES SORPTION ([NICL2] = 10-3 M, 10-5
M, 10-6 M ET 0 M) SUR LE FERRITE DE COBALT A PH 10,5 (B.E. REPRESENTE L’ENERGIE DE
LIAISON DES ELECTRONS) .................................................................................................. 90
FIGURE V - 5 : VARIATION DU POTENTIEL ZETA DU FERRITE DE COBALT EN FONCTION DU PH
POUR DIFFERENTES CONCENTRATIONS INITIALES EN CHLORURE DE NICKEL (25°C, [KCL] =
10-4 M, [NICL2] = 10-4 M, 5.10-5 M, 10-5 M, 10-6 M ET 0 M, [COFE2O4] = 0,4 G.L-1) ...... 91
177
FIGURE V - 6 : VARIATION DE LA CONCENTRATION EN NICKEL SORBE SUR LE FERRITE DE
COBALT EN FONCTION DU PH POUR DIFFERENTS COMPLEXES DE SURFACE (25°C, [KCL] =
10-4 M, [NICL2] = 10-4 M, [COFE2O4] = 0,4 G.L-1, 8 H) .................................................... 93
FIGURE V - 7 : EVOLUTION DE LA CHARGE DE SURFACE DU FERRITE DE COBALT EN PRESENCE DE
-6
CHLORURE DE NICKEL A 10 M EN FONCTION DU PH ........................................................ 94
FIGURE V - 8 : EVOLUTION DE LA QUANTITE DE NICKEL SORBEE EN MOLES PAR LITRE DE
-6 -5 -5 -4
SOLUTION (A, B, C, D : 10 M, 10 M, 5.10 M, 10 M RESPECTIVEMENT) SUR LE FERRITE
DE COBALT. (VALEURS EXPERIMENTALES () ET MODELISATION DE LA QUANTITE SORBEE
TOTALE () PAR COMPLEXATION (- -) ET PAR PRECIPITATION (- -)) (25°C, [KCL] = 10
-4
= 6) ................................................................................................................................... 98
FIGURE V - 10 : VARIATION DE LA CONCENTRATION EN COBALT SORBE SUR LE FERRITE DE
NICKEL EN FONCTION DU PH (25°C, [KCL] = 10 M, [COCL2] = 10 M, 5.10 M, 10 M
-4 -4 -5 -5
........................................................................................................................................ 105
FIGURE V - 17 : EVOLUTION DE LA QUANTITE DE NICKEL SORBEE SUR LA MAGNETITE EN
FONCTION DU TEMPS (25°C, [KCL] = 10 M, [NICL2] = 10 M, [FE3O4] = 0,4 G.L , PH 6)
-4 -6 -1
........................................................................................................................................ 106
FIGURE V - 18 : VARIATION DE LA CONCENTRATION EN COBALT SORBE SUR LA MAGNETITE EN
FONCTION DU PH (25°C, [KCL] = 10 M, [COCL2] = 10 M, 10 M ET 10 M,
-4 -4 -5 -6
178
FIGURE V - 19 : VARIATION DE LA CONCENTRATION EN COBALT SORBE EN MOLES PAR LITRE DE
SOLUTION SUR LA MAGNETITE EN FONCTION DU PH (25°C, [KCL] = 10 M, [COCL2] = 10
-4 -6
M, [COCL2] = 10-4 M, 5.10-5 M, 10-5 M, 10-6 M ET 0 M, [FE3O4] = 0,4 G.L-1) ............... 110
FIGURE V - 23 : VARIATION DU POTENTIEL ZETA DE LA MAGNETITE EN FONCTION DU PH POUR
DIFFERENTES CONCENTRATIONS INITIALES EN CHLORURE DE NICKEL (25°C, [KCL] = 10
-4
M, [NICL2] = 10-4 M, 5.10-5 M, 10-5 M, 10-6 M ET 0 M, [FE3O4] = 0,4 G.L-1) ................ 111
FIGURE V - 24 : VARIATION DE LA CONCENTRATION EN COBALT (A, B, C : 10-6 M, 10-5 M, 10-4
M RESPECTIVEMENT) SORBE SUR LA MAGNETITE. (DONNEES EXPERIMENTALES () DE LA
MODELISATION DE LA QUANTITE SORBEE TOTALE () PAR COMPLEXATION (- - ET - -) ET
PAR PRECIPITATION (- -)) (25°C, [KCL] = 10 M, [FE3O4] = 0,4 G.L , 8 H) .................. 112
-4 -1
= 6) ................................................................................................................................. 115
FIGURE V - 27 : VARIATION DE LA CONCENTRATION EN NICKEL SORBE EN MOLES PAR LITRE DE
SOLUTION SUR LE FERRITE DE COBALT EN FONCTION DU PH ( : 25°C, : 50°C, [KCL] =
10-4 M, [NICL2] = 10-6 M, [COFE2O4] = 0,4 G.L-1, 8 H) .................................................. 116
FIGURE V - 28 : EVOLUTION DE LA QUANTITE DE COBALT SORBEE EN MOLES PAR LITRE DE
SOLUTION SUR LE FERRITE DE NICKEL EN FONCTION DU TEMPS (50°C, [KCL] = 10 M,
-4
........................................................................................................................................ 118
FIGURE V - 31 : EVOLUTION DE LA QUANTITE DE NICKEL SORBEE SUR LA MAGNETITE EN
FONCTION DU TEMPS (50°C, [KCL] = 10 M, [NICL2] = 10 M, [FE3O4] = 0,4 G.L , PH 6)
-4 -6 -1
........................................................................................................................................ 119
FIGURE V - 32 : VARIATION DE LA CONCENTRATION EN COBALT SORBE SUR LA MAGNETITE ( :
25°C, : 50°C, [KCL] = 10-4 M, [COCL2] = 10-6 M, [FE3O4] = 0,4 G.L-1, 8 H) ............. 119
179
FIGURE V - 33 : VARIATION DE LA CONCENTRATION EN NICKEL SORBE SUR LA MAGNETITE ( :
25°C, : 50°C, [KCL] = 10-4 M, [COCL2] = 10-6 M, [FE3O4] = 0,4 G.L-1, 8 H) ............. 120
FIGURE V - 34 : EVOLUTION DE LA QUANTITE DE NICKEL SORBEE SUR LA MAGNETITE EN
FONCTION DU PH (25°C, [KCL] = 10 M, [NICL2] = 10 M / [FE3O4] = 0,4 ET [NICL2] =
-4 -6
180
FIGURE VI - 12 : EVOLUTION DU POTENTIEL ZETA ET PIE DE L’INCONEL 690 AVANT ET APRES
VIEILLISSEMENT ([NANO3] = 10 M) EN AUTOCLAVE ET DE LA MAGNETITE ([KCL] = 10
-3 -4
181
10-4 M) EN PRESENCE DE CHLORURE DE COBALT EN FONCTION DU PH ([NIFE2O4] = 0,4 G.L-
1
, [COCL2] = 0 M ET 10-4 M, 25°C)................................................................................. 164
FIGURE 8 - 6 : EVOLUTION DU POTENTIEL ZETA ET PIE DE L’ACIER INOXYDABLE 304L AVANT
VIEILLISSEMENT ([NANO3] = 10 M) EN AUTOCLAVE ET DE LA MAGNETITE ([KCL] = 10
-3 -4
182
LISTE DES TABLEAUX
183
TABLEAU V - 3 : LOGARITHME DES CONSTANTES DE COMPLEXATION ET DE SOLUBILITE,
DETERMINEES A L’AIDE DU LOGICIEL ECOSAT (2-PK/BSM) POUR LES DIFFERENTS
SYSTEMES ........................................................................................................................ 114
TABLEAU V - 4 : PZC DU FERRITE DE COBALT EN FONCTION DE LA TEMPERATURE (MILIEU KCL
10-2 M) (BARALE-2008).................................................................................................. 116
TABLEAU V - 5 : PZC DU FERRITE DE NICKEL EN FONCTION DE LA TEMPERATURE (MILIEU KCL
10-2 M) (BARALE-2008).................................................................................................. 118
TABLEAU V - 6 : VALEURS DE PZC DE LA MAGNETITE A DIFFERENTES TEMPERATURES
(MANSOUR-2007) ........................................................................................................... 120
TABLEAU V - 7 : LOGARITHME MOYEN DES CONSTANTES DE COMPLEXATION ET DE SOLUBILITE,
DETERMINEES A L’AIDE DU LOGICIEL ECOSAT (2-PK/BSM) POUR LES DIFFERENTS
SYSTEMES ........................................................................................................................ 122
TABLEAU VI - 1 : COMPOSITION DES ACIERS INOXYDABLES ETUDIES PAR BOULANGE-
PETERMANN ET AL. (POURCENTAGE MASSIQUE) ............................................................. 131
TABLEAU VI - 2 : POINTS ISOELECTRIQUES DE METAUX OBTENUS PAR MESURE DE LA VITESSE
D’ADHESION DE PARTICULES DE LATEX A 20°C (CONCENTRATION EN PARTICULES DE
-1 -3
LATEX : 8,8 MG.L , 10 M EN NAOH, RAMIDINE = 38 NM, RCARBOXYLE = 148 NM) (KALLAY-
1991)............................................................................................................................... 133
TABLEAU VI - 3 : PIE DE L’INCONEL 690, DU ZIRCALOY 4 ET DE L’ACIER INOXYDABLE 304L
AVANT ET APRES VIEILLISSEMENT EN AUTOCLAVE, OBTENUS PAR ZETAMETRIE .............. 135
TABLEAU VI - 4 : RECAPITULATIF DES PIE DES MATERIAUX METALLIQUES (INCONEL 690,
ZIRCALOY 4, ACIER INOXYDABLE 304L) DETERMINES DANS LA THESE ........................... 138
TABLEAU VI - 5 : RECAPITULATIF DES PIE DES PRODUITS DE CORROSION DETERMINES DANS LA
THESE EN FONCTION DE LA CONCENTRATION EN IONS NICKEL OU COBALT EN SOLUTION 138
TABLEAU VI - 6 : PARAMETRES CHOISIS POUR REALISER LES CALCULS AVEC LE MODULE
LAGRANGIEN ET TAUX DE DEPOT NORMALISES EN RESULTANT (TOUS LES RESULTATS POUR
LES CHARGES DE MEME SIGNE SONT DONNES POUR DES POTENTIELS (ϕI) NEGATIFS,
CEPENDANT POUR LES VALEURS POSITIVES, LES TAUX DE DEPOT SONT IDENTIQUES) ...... 139
TABLEAU VI - 7 : DETERMINATION DES GAMMES DE PH POUR LESQUELLES LE DEPOT DE
PARTICULES DE FERRITE DE COBALT, AVANT OU APRES SORPTION DE CHLORURE DE NICKEL
(0 M ET 10-4 M), SUR L’INCONEL 690 VIEILLI (V) OU NON (NV) EN AUTOCLAVE, A LIEU
(DIAMETRE DE LA PARTICULE : 50 NM). O SIGNIFIE QUE LE DEPOT EST FAVORISE, N QUE
NON (ND (NON DEFINI) : DEPEND DE LA CONTINUITE DES COURBES) .............................. 142
TABLEAU VI - 8 : DETERMINATION DES GAMMES DE PH POUR LESQUELLES LE DEPOT DE
PARTICULES DE FERRITE DE NICKEL, AVANT OU APRES SORPTION DE CHLORURE DE COBALT
(0 M ET 10-4 M), SUR L’INCONEL 690 VIEILLI (V) OU NON (NV) EN AUTOCLAVE, A LIEU
(DIAMETRE DE LA PARTICULE : 50 NM). O SIGNIFIE QUE LE DEPOT EST FAVORISE, N QUE
NON (ND (NON DEFINI) : DEPEND DE LA CONTINUITE DES COURBES) .............................. 144
TABLEAU VI - 9 : DETERMINATION DES GAMMES DE PH POUR LESQUELLES LE DEPOT DE
PARTICULES DE MAGNETITE, AVANT OU APRES SORPTION DE CHLORURE DE COBALT OU DE
-4
NICKEL (0 M ET 10 M), SUR L’INCONEL 690 VIEILLI (V) OU NON (NV) EN AUTOCLAVE, A
LIEU (DIAMETRE DE LA PARTICULE : 50 NM). O SIGNIFIE QUE LE DEPOT EST FAVORISE, N
QUE NON (ND (NON DEFINI) : DEPEND DE LA CONTINUITE DES COURBES) ....................... 145
184
TABLEAU VII - 1 : LOGARITHMES DES CONSTANTES MOYENNES DE COMPLEXATION ET DE
SOLUBILITE, DETERMINEES A L’AIDE DU LOGICIEL ECOSAT (2-PK/BSM) POUR LA
SORPTION DU COBALT ET DU NICKEL SUR LES TROIS PRODUITS DE CORROSION ETUDIES
(25°C) ............................................................................................................................. 148
TABLEAU 4 - 1 : LD ET LQ DU COBALT ET DU NICKEL LORS DES ANALYSES PAR ICP-AES ET
ICP-MS........................................................................................................................... 154
185