TheseSPerrin Vep

ECOLE NATIONALE SUPERIEURE N° d’ordre : 296CD
DES MINES DE SAINT-ETIENNE
THESE
Présentée par
Stéphane PERRIN
pour obtenir le grade de
DOCTEUR
DE L’ECOLE NATIONALE SUPERIEURE DES MINES DE SAINT-ETIENNE
(Spécialité : GENIE DES PROCEDES)
Modélisation de la cinétique de transformations non isothermes et

(ou) non isobares. Application à la déshydroxylation de la
kaolinite et à la réduction de l'octooxyde de triuranium par
l'hydrogène
Soutenue à Saint-Etienne le 19 Décembre 2002
COMPOSITION du JURY
Monsieur D. ABLITZER Président

Rapporteur
Madame F. ROUQUEROL Rapporteur
Messieurs J. M. CRIADO Examinateurs

F. NICOLAS
M. SOUSTELLE
Madame M. PIJOLAT
Mesdames C. ABLITZER-THOUROUDE Invités

M. PRIN-LAMAZE
Monsieur L. CARRARO
Thèse préparée au Centre « Sciences des Processus Industriels et Naturels »

REMERCIEMENTS
Je tiens à remercier messieurs René LALAUZE et Michel COURNIL pour m’avoir

accueilli au centre de recherche “ Sciences des Processus Industriels et Naturels ”.
Je voudrais adresser mes plus sincères remerciements Michelle PIJOLAT pour son
encadrement, sa disponibilité, ses encouragements et son apport scientifique. Je tiens
également à exprimer ma gratitude à Michel SOUSTELLE pour son encadrement scientifique
ainsi que pour son soutien. J’adresse toute ma reconnaissance à Françoise VALDIVIESO
pour l’intérêt qu’elle a porté à cette étude et pour les nombreuses discussions scientifiques
qui m’ont permis de progresser.
Je remercie également Mesdames Françoise ROUQUEROL, Carine ABLITZER-

THOUROUDE et Marie PRIN-LAMAZE, Messieurs Denis ABLITZER, François NICOLAS,
José Manuel CRIADO et Laurent CARRARO pour avoir accepté de juger du contenu
scientifique de ce travail.
Mes remerciements s’adressent également à l’équipe des mathématiciens du centre

SIMMO, Eric TOUBOUL, Anca BADEA et Céline HELBERT, avec qui la collaboration a
toujours été fructueuse.
Je n’oublierai pas non plus l’ensemble du personnel du centre SPIN dont l’aide, la
disponibilité et les compétences ont contribué à faire avancer ce travail. Je remercie plus
particulièrement Nathalie BERNE d’avoir accepté de relire cette thèse et pour son aide à la
mise en page.
Enfin je voudrais remercier les personnes qui m’ont soutenu durant ces trois années,
Tuyet, ma mère, ma famille et mes amis.
Sommaire
Introduction 1
Chapitre I : Etude bibliographique 5
I.1 Modèles de transformation et lois cinétiques 6
I.1.1 Lois cinétiques de germination 7
I.1.2 Lois cinétiques de croissance 11
I.1.3 Modèles de transformation 19
I.1.4 Influence de la pression et de la température 28
I.1.4.1 Approche utilisant la séparation des variables 29

I.1.4.1.1 Influence de la pression 30
I.1.4.1.2 Influence de la température 30
I.1.4.2 Approche sans supposer la séparation des variables 31
I.1.4.2.1 Détermination de la réactivité surfacique de croissance
en fonction de la température et des pressions partielles
des gaz 32
I.1.4.2.2 Expressions de φ(T,P) 34
I.1.5 Effet d'une distribution granulométrique 39
I.2 Modélisation en condition non isotherme 41
I.2.1 Objectifs 41
I.2.2 Classification des méthodes non isothermes 42
I.2.2.1 Programmation en température 43

I.2.2.2 Analyse thermique à vitesse contrôlée (ATVC) 48
I.3 Conclusions 50
Chapitre II : Approche non isotherme et non isobare de

la modélisation cinétique 53
II.1 Intérêts de la modélisation 54
II.2 Principe de la modélisation en conditions non isotherme

et non isobare 58
II.2.1 Cas des réactions avec germination-croissance 58
II.2.2 Cas des réactions successives 62
II.3 Choix des systèmes réactionnels 63
Chapitre III : Modélisation d’une réaction avec germination-

croissance : la déshydroxylation de la kaolinite 65
III.1 Etude bibliographique 66
III.1.1 Structure de la kaolinite 66
III.1.2 Etude thermodynamique de la réaction 67
III.1.3 Cinétique de la déshydroxylation de la kaolinite 69
III.2 Etude expérimentale 71
III.2.1 Caractérisation des échantillons 71
III.2.1.1 Diffractogrammes de la kaolinite et de la métakaolinite 71

III.2.1.2 Observations microscopiques 73
III.2.1.3 Caractérisation par adsorption gazeuse 75
III.2.1.4 Distribution granulométrique 77
III.2.2 Etude cinétique de la réaction par thermogravimétrie 79
III.2.2.1 Protocole expérimental 79

III.2.2.2 Courbes cinétiques expérimentales 80
III.2.2.2.1 Influence de la pression de vapeur d'eau 80
III.2.2.2.2 Influence de la température 81
III.2.2.3 Expériences réalisées en programmation de température 82
III.3 Modélisation de la réaction 83
III.3.1 Modèle de transformation 83
III.3.1.1 Vérification des hypothèses du modèle 83

III.3.1.1.1 Pseudo-stationnarité de la réaction 83
III.3.1.1.2 Test du φE 85
III.3.1.1.3 Test du cas limite 87
III.3.1.2 Choix du modèle 88
III.3.1.3 Expression de la vitesse 90
III.3.1.3.1 Avec une seule taille de grain 90
III.3.1.3.2 Avec une distribution granulométrique 90
III.3.1.4 Comparaison avec l'expérience 91
III.3.1.4.1 Avec une seule taille de grain 91
III.3.1.4.2 Avec une distribution granulométrique 94
III.3.2. Modélisation physico-chimique 95
III.3.2.1 Variation de φ et γ en fonction des conditions physico-chimiques 95

III.3.2.1.1 Par comparaison du modèle avec les expériences 95
III.3.2.1.2 Par la méthode des décrochements 96
III.3.2.1.3 Validation du modèle de transformation 99
III.3.2.2 Mécanisme réactionnel de croissance 100
III.3.2.3 Mécanisme réactionnel de germination 107
III.3.3 Modélisation en condition non isotherme 110
Chapitre IV : Modélisation de réactions successives :

La réduction de U3O8 par l’hydrogène 121
IV.1 Etude bibliographique 122
IV.1.1 Diagramme thermodynamique 122
IV.1.2 Structure des oxydes d'uranium étudiés 123
IV.1.2.1 L'octooxyde triuranium : U3O8 123

IV.1.2.2 Le dioxyde d'uranium : UO2 124
IV.1.3 Etude cinétique et modélisation de la réaction 126
IV.2 Etude expérimentale 134
IV.2.1 Etude de la réaction par thermogravimétrie 134
IV.2.1.1 Protocole expérimental 134

IV.2.1.2 Evolution de la pression d'hydrogène dans la thermobalance 136
IV.2.1.3 Courbes cinétiques 140
IV.2.1.3.1 Condition non isobare 140
IV.2.1.3.2 Condition non isotherme 141
IV.2.2 Caractérisation texturale des produits 142
IV.2.2.1 Caractérisation par adsorption d'azote à 77 K 142

IV.2.2.1.1 Isothermes d'adsorption et de désorption 143
IV.2.2.1.2 Mesure de la surface spécifique 143
IV.2.2.2 Caractérisation par observations microscopiques 144
IV.2.2.2.1 Microscopie à force atomique (AFM) 144
IV.2.2.2.2 Microscopie électronique à transmission (MET) 146
IV.2.2.3 Distribution granulométrique 147
IV.2.3 Caractérisation structurale des phases au cours de la réaction 148
IV.2.3.1 Diffractogrammes de U3O8 et de UO2 148

IV.2.3.2 Etude du début de la réaction 150
IV.2.3.3 Etude de la fin de la réaction 152
IV.2.3.3.1 Echantillons réduits à 550°C 152
IV.2.3.3.2 Echantillon réduit à 430°C 154
IV.3 Modélisation en conditions non isotherme ou non isobare 155
IV.3.1 Présentation et hypothèses du modèle 155
IV.3.2 Modélisation de chaque transformation 162
IV.3.2.1 Ecart à la stœchiométrie de U3O8 163

IV.3.2.2 Réduction de U3O8-x en U4O9 168
IV.3.2.2.1 Test du φE 168
IV.3.2.2.2 Modèle de transformation 170
IV.3.2.2.3 Modèle physico-chimique 172
IV.3.2.3 Réduction de U4O9 en UO2 179
IV.3.2.3.1 Modèle de transformation 179
IV.3.2.3.2 Modèle physico-chimique 181
IV.3.3 Expression de la vitesse théorique 184
IV.3.3.1 Avec une seule taille de grains 185

IV.3.3.2 Avec une distribution granulométrique 185
IV.3.4 Confrontation modèle-expérience 191
IV.3.4.1 Changement de la stœchiométrie de U3O8 192

IV.3.4.2 Condition non isobare 195
IV.3.4.3 Condition non isotherme 199
Conclusion 203
Annexes 207
Annexe A : Expression de la vitesse de croissance en considérant

les hypothèses de pseudo-stationnarité et de l'étape limitante 207
Annexe B : Pseudo-stationnarité d'un système réactionnel 213
Annexe C : Test du φΕ 217
Annexe D : Comparaison modèle-expérience. Utilisation du paramètre de modèle 219
Annexe E : Modèle de germination-croissance anisotrope à développement interne

avec une étape limitante de croissance de diffusion et
des grains cylindriques 223
Annexe F : Expression de la vitesse pour la variation de stœchiométrie de

l'octooxyde de triuranium 231
Annexe G : Modèle de germination instantanée et croissance anisotrope à

développement interne avec une réaction limitant la croissance
située à l'interface externe pour des grains sphériques (de rayon r0) 239
Références bibliographiques 247

Introduction
Le traitement thermique de solides divisés joue fréquemment, dans les procédés

industriels, un rôle crucial sur la qualité d'usage des produits.
Ces réactions se déroulent à des températures élevées, dans des fours de grandes
dimensions, et mettent en général en jeu un ou plusieurs gaz réactifs ou produits. Les
caractéristiques (composition, surface spécifique ...) des produits formés au cours de ce
traitement dépendent notamment de la manière dont la réaction s'est déroulée, et en particulier
de l'évolution de la température et des variations des pressions partielles des gaz. Afin de
mieux comprendre et également d'améliorer le fonctionnement de ces réacteurs industriels, il
est nécessaire de développer des outils permettant de simuler les phénomènes qui se déroulent
en leur sein.
Ainsi, un programme de recherches s'est fixé pour mission de mettre en place ces
outils de modélisation des fours industriels. Ce projet est réalisé grâce à la participation de
partenaires universitaires (le centre SPIN et le centre SIMMO de l'Ecole Nationale Supérieure
des Mines de Saint Etienne et le centre LSG2M de l'Ecole Nationale Supérieure des Mines de
Nancy) soutenus par plusieurs partenaires industriels (Comurhex, CEA, CRV Pechiney).
Afin d'expliquer le fonctionnement de ces fours, un certain nombre de bilans sont à

prendre en compte. Nous pouvons distinguer deux types de termes entrant dans ces bilans :
• les termes d'échanges relatifs à l'ensemble des processus non réactionnels au sein du
réacteur : écoulement du solide et des gaz, transfert des gaz au sein du lit de poudre,
échanges thermiques avec les parois, …. ;
• les termes sources ou les vitesses de réactions qui expriment les quantités de matière et de
chaleur créées ou consommées par les réactions.
Ainsi, la réalisation de ce projet peut être séparée en deux parties : une première étude
doit s'intéresser à la modélisation du réacteur, et une deuxième étude à la modélisation de la
réaction. Ces deux parties devront par la suite être reliées puisque, comme le montre la
figure 1, chaque modèle interagit sur l'autre.
-1-
Variations de température et des
pressions engendrées par la réaction
Modèle réacteur Modèle réaction
Conditions de température et de
pressions de la réaction au sein du four
Figure 1 : Schéma de l'interaction entre les modèles développés
Notre étude porte sur la détermination des termes sources, donc sur la modélisation de
la réaction lorsque les conditions physico-chimiques varient au cours du temps, comme c'est
le cas au sein d'un réacteur.
Ces réactions procèdent par la superposition de deux phénomènes : la germination qui

crée de petits îlots de la nouvelle phase, suivie par la croissance de ces petites étendues. Selon
les cas, un seul des deux processus est à prendre en considération (l'autre étant instantané), ou
bien il existe une compétition entre la germination et la croissance.
La vitesse à un instant donné dépend des conditions physico-chimiques mais peut,
dans le deuxième cas, également dépendre du passé thermique de la réaction c'est-à-dire des
conditions dans lesquelles la réaction s'est déroulée. En effet, considérons l'exemple idéalisé
de la figure 2, illustrant deux grains en cours de transformation. Chaque grain est supposé
avoir subi des conditions de transformations différentes. Ainsi, le premier grain est soumis à
des conditions physico-chimiques telles qu'il n'a pu se former qu'un seul germe. Pour le
deuxième, les conditions ont engendré la formation de deux germes. La phase A est la phase
initiale, et la phase B la phase finale. Nous supposons ici que la vitesse est proportionnelle à
l'aire de l'interface A/B.
A B
A B
B A
Figure 2 : Grains en cours de réaction, influence du passé thermique
-2-
Les deux grains sont au même degré d’avancement 0,5. Mais l’interface entre les deux
phases est double dans le deuxième cas. Ainsi, la vitesse est également deux fois plus grande
dans cet exemple.
Le degré d'avancement ne permet donc pas toujours de déterminer la vitesse et les
conditions du « passé » ont une influence sur la vitesse, ce qui rend difficile la résolution des
équations bilans.
Par ailleurs, la modélisation des réacteurs industriels est généralement basée sur
l'hypothèse d'une réaction unique à partir de laquelle est exprimée la vitesse.
Aussi, notre travail a porté sur la modélisation de deux types de transformations :
- une réaction avec compétition entre la germination et la croissance (comme le cas illustré
schématiquement par la figure 2),
- une transformation mettant en jeu trois étapes réactionnelles successives (mais pas de
compétition entre germination et croissance),
dans un champ de contraintes physico-chimiques variant avec le temps, et en tenant compte
de l'influence de ces variations sur la transformation.
Afin d'atteindre cet objectif, cette étude s'articule autour de quatre chapitres :
• Le premier chapitre est consacré à la description des études de modélisation réalisées en

condition non isotherme en cinétique du solide. Pour cela, nous détaillons préalablement
les principaux modèles développés pour décrire une réaction de type solide-gaz. Cette
première partie aborde également les carences de la modélisation développée en
programmation de température.
• Le deuxième chapitre met en avant la méthodologie utilisée afin de réaliser une approche
cinétique en conditions non isotherme et non isobare. Nous montrerons notamment que
celle-ci s'appuie sur le développement de deux modèles : un modèle de transformation et
un modèle physico-chimique. Cette méthodologie est ensuite employée dans l'étude des
réactions décrites dans les chapitres III et IV.
• Ainsi, le troisième chapitre présente un modèle pour une décomposition thermique, la

déshydroxylation de la kaolinite, et l'application de ce modèle lorsque la température
-3-
varie au cours de cette réaction. Ce modèle correspond à un cas où les processus de
germination et de croissance interviennent dans la loi de vitesse.
• Enfin, le quatrième et dernier chapitre traite de la modélisation cinétique en conditions

non isotherme et non isobare de la réduction de l'octooxyde de triuranium en dioxyde
d'uranium par l'hydrogène. Comme nous l'aborderons dans cette partie, la réaction étudiée
se divise en plusieurs transformations successives. Nous verrons ainsi comment prendre
en compte ces différentes transformations dans la modélisation non isotherme et non
isobare.
-4-
Chapitre I Etude bibliographique
C
Chhaappiittrree II
ETUDE BIBLIOGRAPHIQUE
Les réactions de transformation des solides interviennent dans de nombreux procédés

industriels.
L'étude de la cinétique de ces réactions présente d'abord un intérêt fondamental pour
comprendre les mécanismes intervenant dans ces réactions. Ainsi de nombreuses études ont
pour objectifs d'expliquer comment se déroulent ces réactions afin d'obtenir les lois de vitesse
de celles-ci en fonction du temps et également parfois en fonction des contraintes physico-
chimiques de la réaction (température, pressions des gaz).
Ensuite ces études cinétiques possèdent également un intérêt industriel puisqu'elles
doivent permettre de bien maîtriser, voire d'améliorer les procédés et les réacteurs qui mettent
en jeu ces réactions.
Ces études cinétiques ont d'abord été réalisées en conservant tout au long de ces
transformations les conditions expérimentales (température, pressions des gaz ...) constantes.
Elles ont permis le développement de modèle expliquant l'évolution de la vitesse de la
réaction au cours de celle-ci.
Dans un deuxième temps, en s'appuyant sur les lois cinétiques mises en place en
condition isotherme, de nombreux travaux ont consisté à développer des méthodes de travail
sur la cinétique des transformations de solides en faisant varier la température en cours de
réaction. L'avantage visé de ces études est en particulier de pouvoir étudier une réaction sur
un grand domaine de température en réalisant peu d'expériences.
Afin de rendre compte de ces travaux, nous présentons dans une première partie de
cette étude bibliographique les lois cinétiques et les modèles de transformation utilisés dans la
littérature ainsi que l'influence de la température et de la pression sur la vitesse de la réaction.
Nous nous intéressons essentiellement aux solides se présentant sous forme de poudres et
-5-
nous montrons également, à travers certaines études, l'influence d'une distribution

granulométrique sur les courbes cinétiques.
Dans une deuxième partie, nous développons les travaux menés sur la modélisation de
ces réactions en condition non isotherme.
Nous terminons en mettant en avant les problèmes liés à ces études.
Ces réactions hétérogènes, c'est-à-dire où les produits réagissant et formés par la

réaction appartiennent à des phases différentes, possèdent une grande variété selon le nombre
de phases qui interviennent. Dans le cadre de cette étude, nous étudions principalement deux
types de réactions avec des solides :
• les décompositions thermiques où un solide S1 va se transformer en solide S2 avec la
formation d'un gaz G : S1→ S2 + G ;
• les réactions où un solide S1 va réagir avec un gaz G1 pour former un solide S2 et un gaz
G2 : S1 + G1→ S2 + G2.
I.1 Modèles de transformation et lois cinétiques
Au cours de la transformation d'un solide A en solide B, deux processus

interviennent : la germination et la croissance.
La germination est la fabrication de petites étendues de la phase B, appelées germes.

Ces germes peuvent naître en surface ou dans le volume du solide A. Ce processus correspond
à la création de l'interface réactionnelle.
La croissance est la progression de l'interface qui conduit au grossissement des

germes formés.
Ces deux processus contribuent à l'évolution de l'interface réactionnelle qui caractérise

la vitesse de la réaction. Ainsi, les modèles de transformation décrivant ces réactions
dépendent de la manière dont ces deux phénomènes vont se dérouler.
-6-
I.1.1 Lois cinétiques de germination
D’après Galwey et al. [1], la germination peut être décrite par une succession de deux
étapes différentes :
• la transformation chimique de certains constituants de la phase initiale en constituants du
produit de réaction,
• l’association de ces constituants du produit en réseau cristallin afin de former la nouvelle
phase.
Plusieurs lois ont été développées afin de décrire la vitesse de formation des germes.
La détermination expérimentale de cette vitesse n'est pas aisée et repose en général sur des
observations microscopiques en début de réaction.
Nous développons dans ce paragraphe les principales lois décrivant la germination et
leur interprétation. Celles-ci ont été répertoriées dans divers ouvrages de cinétique hétérogène
[1,2,3,4,5,6,7].
Nous notons N le nombre de germes formés à l'instant t et N0 le nombre total de
germes qui vont se former. Pour chaque cas, les expressions du nombre de germes et de la
vitesse de formation de ces germes seront données.
1 Germination instantanée
Dans ce cas, tous les germes sont formés au même moment, au début de la réaction. La
loi de germination s'écrit donc sous la forme :
A t<0 N=0
A t ≥0 N=N0 (I.1)
2 Germination à vitesse constante ou linéaire
La vitesse de germination est indépendante du temps et quel que soit l'avancement de

la réaction, la probabilité pour qu'un germe se forme est la même.
L'expression de la loi de germination s'écrit sous la forme :
-7-
dN =k (I.2)
dt g
N =kgt (I.3)
kg représente la probabilité pour qu'un germe apparaisse pendant l'unité de temps.
Cette loi est également appelée germination d'ordre 0.
3 Germination d'ordre 1 ou en une seule étape
Cette loi correspond au cas où on suppose l'existence de sites potentiels de

germination. Ceux-ci ont la même probabilité d'être activés et donc de former un germe. La
vitesse de germination et le nombre de germes s’expriment par :
dN =k (N − N) (I.4)
dt g 0
N = N0(1−exp(−kgt)) (I.5)
où kg est la probabilité qu'un site potentiel de germination soit activé par unité de temps.
Le nombre de sites potentiels et par conséquent la vitesse de germination vont
décroître avec le temps.
4 Germination selon une loi puissance ou en plusieurs étapes
La germination est dans ce cas due à une succession d’évènements ou se déroule par
l’association de plusieurs « individus » chimiques. Les expressions de la vitesse de
germination et du nombre de germes sont données par :
dN
= Kt n −1 (I.6)
dt
N = Kt n (I.7)
où n est le nombre d’évènements nécessaires à la formation d’un germe et K représente le

produit des probabilités que chaque évènement se produise si on considère une succession
d'évènements. Dans le cas d'une association de plusieurs entités chimiques, n-1 est le nombre
d’entités chimiques qui en s’associant vont former un germe et K est la probabilité totale pour
que chaque espèce s’associe.
-8-
5 Germination en chaînes branchées
On suppose que la germination (et éventuellement la croissance du germe formé)

provoque la naissance d’autres germes. Ce phénomène est appelé branchement ou
ramification des germes. Cette théorie a été développée pour expliquer des accélérations
brutales de la vitesse en début de certaines réactions. Les équations relatives à la germination
sont de la forme :
dN
= k g N 0 exp( k b t ) (I.8)
dt
k g N0
N= (exp( k b t ) − 1) (I.9)
kb
où kg représente la probabilité qu’un germe se forme par unité de temps selon un processus
spontané et kb est la probabilité qu’un germe se forme par unité de temps due au phénomène
de « branchement ».
Delmon [2] reprend les lois précédentes en distinguant celles de probabilité uniforme,
c'est-à-dire où un germe peut se former en n'importe quel point de la surface ou du volume du
solide, aux lois de germination où les germes se forment sur des sites particuliers localisés en
des points biens définis, appelés sites potentiels. Il explique que sauf dans les transformations
polymorphiques, l'interface réactionnelle ne peut naître que sur la surface du réactif
solide. Afin de relier le nombre de germes à la surface du solide, il définit une vitesse
spécifique de germination dγ/dt comme le nombre de germes formés par unité de temps sur
l'unité de surface.
Toutes ces lois de cinétique de germination ne peuvent en aucun cas être utilisées à un
instant quelconque quand le processus de croissance est entamé. Nous verrons qu'il s'agit en
fait de réactivité de germination (stationnaire ou non) et qu'il faudra leur adjoindre, pour
atteindre la vitesse de germination, l'étendue de l'espace réactionnel concerné c'est-à-dire la
surface libre de germination à chaque instant.
Cependant trop peu d'expressions de lois se raccordent aux mécanismes de
germination proprement dits et de ce fait, ne permettent pas de prendre en compte l'influence
des variables intensives en particulier les pressions partielles de gaz sur la germination.
-9-
Certains [1,8] proposent de faire un bilan sur les germes potentiels. Pour cela, ils
écrivent que le nombre de germes susceptibles de se former à un instant t, N1(t), est égal au
nombre de germes potentiels total, N0, auquel on soustrait les germes ayant déjà réagi, N(t), et
le nombre de germes, notés N2(t), ayant été éliminés par recouvrement au cours de la
croissance des germes précédents :
N1(t)= N0 − N(t)− N2(t) (I.10)
En différenciant cette expression, il est possible d'obtenir :
−dN1=dN +dN2 (I.11)
Or la variation du nombre de germes qui réagissent à un instant t dépend du nombre de
germes potentiels restant (germination d'ordre un) :
dN =kg N1dt (I.12)
La variation des germes potentiels « ingérés », dN2, est exprimée par :
dN2 = 1 dα

N
(I.13)
 1−α 
Ainsi à partir des équations (I.10) à (I.13), il est possible de déterminer le nombre de
germes qui se forment en fonction du temps. Cependant le degré d'avancement, α, présent
dans l'équation (I.13), est caractéristique de l'évolution du volume. Cette approche de la
germination se limite donc aux germes apparaissant au sein du volume de solide. Ainsi cette
loi de germination traitant du recouvrement des germes ne permet pas de considérer une
germination en surface.
Afin de tenir compte de la modification du domaine de la réaction due à l'apparition de

germes et à la croissance de ces derniers, Soustelle et Pijolat [9] supposent une germination
en surface et définissent une surface libre de germination, notée SL, qui correspond à la
surface non transformée des grains de la poudre. Ils définissent une fréquence surfacique de
germination, notée γ, comme le nombre de germes formés par unité de surface et par unité de
temps. Celle-ci correspond à la vitesse spécifique de germination définie par Delmon [2].
Cette grandeur dépend des conditions physico-chimiques. La vitesse de germination (nombre
de germes formés par unité de temps) s'exprime alors sous la forme :
v g = dN =γS L (I.14)
dt
- 10 -
L'utilisation d'une surface libre de germination pour décrire la germination est

intéressante puisqu'elle prend en compte l'effet de la croissance des germes sur la vitesse de
germination. C'est cette approche que nous utiliserons dans notre étude.
La formation d'un germe est suivie par sa croissance. Comme pour la germination, il
est nécessaire d'expliquer ce phénomène à l'aide de lois qui permettent de décrire ce deuxième
processus.
I.1.2 Lois cinétiques de croissance
1 Description de la croissance
Comme le montre la figure I.1, la croissance peut avoir la même vitesse quelle que soit
la direction (croissance isotrope), ou alors des directions privilégiées (croissance
anisotrope). Dans ce dernier cas, la composante radiale de la vitesse est seule prise en
compte, la vitesse étant supposée infinie dans les deux directions tangentielles à la surface.
Croissance isotrope Croissance anisotrope

Figure I.1 : Mode de croissance du germe
Dans le cas d'une croissance anisotrope, celle-ci est souvent considérée comme se
déroulant selon une, deux ou trois dimensions de l'espace. Ce nombre de dimensions de
l'expansion de la nouvelle phase dépend en général de la géométrie des grains. Pour des grains
en plaquettes minces (en supposant uniquement l'attaque des champs), la croissance se
déroule selon une dimension de l'espace. Pour des grains cylindriques (cylindres
« allongés »), elle a lieu selon deux dimensions de l'espace (en négligeant la croissance au
niveau des faces des cylindres). Pour des grains de géométrie sphérique, la nouvelle phase se
développe selon les trois dimensions de l'espace. Enfin pour des parallélépipèdes, le
développement de la nouvelle phase se déroule également selon trois dimensions mais avec
des constantes différentes.
- 11 -
Un autre élément qui joue au niveau de la croissance est le sens du développement de

la nouvelle phase. En effet, si on se réfère aux travaux réalisés sur l'oxydation des métaux ou
des alliages par un gaz [10], les études, développées en général en géométrie plane et pour
une germination instantanée, présentent des mécanismes où la croissance de la couche
d'oxyde peut se faire soit vers l'intérieur du métal, soit vers l'extérieur. Ce sens dépend
notamment du type de défauts présents dans l'oxyde (lacunes, atomes en positions
interstitielles ...). On peut citer comme exemple de développement externe, le cas de
l'oxydation du magnésium par l'oxygène [11]. Ainsi comme le présente la figure I.2, la phase
solide produite peut se développer soit vers l'extérieur du grain (développement externe) soit
vers le cœur du grain (développement interne).
développement externe développement interne

Figure I.2 : Schéma du développement de la nouvelle phase
Enfin, la croissance d'un germe se déroule en plusieurs étapes élémentaires qui ont
lieu dans des zones réactionnelles différentes. Ces étapes sont des réactions interfaciales et
des étapes de diffusion [12]. Les réactions interfaciales peuvent se dérouler soit à l'interface
interne (interface solide initial/solide final), soit à l'interface externe (interface solide
final/gaz). Pour établir l'expression de la vitesse de croissance, les auteurs considèrent qu'une
étape est limitante et qu'elle fixe la vitesse de croissance.
1 Expression de la vitesse de croissance
La plupart des études cinétiques considère que le développement de la nouvelle phase

se fait aux dépens de la phase initiale.
Pour de nombreux auteurs, en supposant que l'étape limitante se déroule à l'interface

interne, la vitesse de l'avancement de l'interface s'exprime alors a priori sous la forme [1,4,6] :
dr
= kc (I.15)
dt
- 12 -
où r est une direction de la croissance du germe et kc est la « constante » cinétique pour la

croissance.
Pour Delmon [2], la vitesse de croissance, définie comme la quantité de réactif

transformée par unité de temps, est proportionnelle à l'aire de l'interface réactionnelle entre
les deux phases solides notée Si et s'écrit sous la forme :
Vc = k i S i (I.16)
où ki est la vitesse interfaciale spécifique c'est-à-dire la vitesse de croissance ramenée à
l'unité de surface de l'interface réactionnelle. Il suppose donc que l'étape limitante se situe à
l'interface entre les deux phases solides.
Dans la majorité des études, les auteurs ne s'intéressent pas à la signification de la

constante ki (ou kc) par rapport à l'étape limitant la croissance.
Remarque 1 : les équations (I.15) et (I.16) expriment la vitesse de croissance de deux

manières différentes. Ainsi les constantes ki et kc peuvent être reliées par la relation :
kc =VmAki (I.17)
où VmA est le volume molaire de la phase initiale.
Remarque 2 : de manière générale, on peut noter chez la plupart des auteurs (excepté les
études d’oxydation de métaux), l’absence de la proposition de mécanismes en étapes élémentaires
(avec les intermédiaires réactionnels des deux phases) comme en cinétique homogène.
En se basant sur des mécanismes de croissance en étapes élémentaires, Soustelle et

Pijolat [9] proposent d’exprimer la vitesse sous la forme :
Vc =n0φE (I.18)
où n0 est la quantité initiale du solide de départ ;

φ est la réactivité surfacique de croissance de l'étape limitante et s'exprime en
mol.m-2.s-1 ;
E est appelée fonction d’espace et dépend de l'évolution au cours de la réaction de la
zone réactionnelle où se produit l'étape limitant la croissance et s'exprime en m2.mol-1.
- 13 -
Comme on le verra dans la suite, cette écriture offre de nombreux avantages, dont en
particulier :
• la décomposition de la vitesse en deux fonctions séparées ;
• pour chaque étape élémentaire d’un mécanisme, on peut calculer les expressions des deux
fonctions φ et E.
Quelle que soit l’approche considérée, il apparaît que l’hypothèse d’une étape
limitante est à la base des modèles de croissance. Le paragraphe suivant montre que cette
hypothèse n’est valable que dans certaines conditions, et qu’il est possible de la vérifier
expérimentalement.
1 Hypothèse de l'étape limitante de croissance
Pour établir une loi de vitesse de croissance, les auteurs supposent donc que la vitesse
est imposée par une étape particulière du mécanisme, se produisant dans une zone
réactionnelle donnée. L’équation (I.18) est ainsi une généralisation de ce principe
communément (et parfois implicitement) admis.
Un mécanisme faisant intervenir généralement plusieurs étapes élémentaires et

plusieurs intermédiaires réactionnels, il est en général impossible de résoudre le système
d’équations différentielles déduit des bilans de matière à tout instant. Pour simplifier, il est
possible d’avoir recours à des hypothèses dont celles de la stationnarité et de l’étape limitante.
Soustelle [12] précise que pour les transformations de solides, il est nécessaire de
considérer des systèmes « pseudo-stationnaires » pour lesquels les concentrations en espèces
intermédiaires ([Xi]) ne varient pas au cours du temps et la variation du logarithme du volume
(Vi) ou de la surface des zones réactionnelles (Si) évolue peu dans le temps. Ces deux
conditions s'écrivent sous la forme :
∂[ X i ] ∂ ln( S i ) ∂ ln( Vi )
= 0 et ≈ 0 ou ≈0 (I.19)
∂t ∂t ∂t
Ces deux conditions peuvent également se traduire par la faible variation des quantités
δni/δ
de chaque intermédiaire réactionnel (δ δt≈
≈0 où ni est le nombre de mole d'un intermédiaire
- 14 -
i). Il est possible de démontrer que les vitesses des étapes élémentaires constituant le
mécanisme réactionnel d'un système pseudo-stationnaire sont liées entre elles [12].
Dans l’annexe A, nous montrons que la vitesse de croissance peut s’exprimer par
l’équation (I.18), lorsque les hypothèses de pseudo-stationnarité et de l’étape limitante sont
satisfaites.
La vérification de l’hypothèse de pseudo-stationnarité a été faite pour la première fois

(à notre connaissance) par Lalauze et al. [13] pour l'oxydation du niobium par l'oxygène. Le
principe, démontré en annexe B, est simple : les courbes de vitesse d'un système pseudo-
stationnaire obtenues par deux méthodes de mesure différentes, doivent être affines
entre elles. Cette méthode a notamment été développée en réalisant des couplages
thermogravimétrie-calorimétrie. Elle a été appliquée à l’étude de la décarbonatation du
carbonate de calcium [14], de l'oxydation du magnésium par l'oxygène [11] et de l'oxydation
d'alliages de zirconium par l'oxygène et par la vapeur d'eau [15].
Il est également possible de vérifier l’hypothèse de l’étape limitante de la croissance.

En effet, en se basant sur l’expression (I.18) de la vitesse, Soustelle et Pijolat [9] proposent
une méthode expérimentale basée sur des « décrochements », c'est-à-dire la variation brutale
de l'un des paramètres physico-chimiques de l'expérience. Cette méthode, appelée test du φE,
est décrite dans l'annexe C. Pour chaque expérience, la vitesse est mesurée avant le
décrochement (notée (dα/dt)av) et après le décrochement (notée (dα/dt)ap). Comme le montre
l'annexe C, si la vitesse est décomposable sous la forme φE alors quel que soit l'instant
du décrochement le rapport (dα dt)ap (dα dt)av est le même.
Cette méthode a notamment été utilisée par Viricelle et al. [16] pour la transformation
de l'hydroxycarbonate de cérium en cérine et par Bouineau [14] pour la décomposition du
carbonate de calcium.
Remarque : le test du φE peut également être vérifié pour un cas mixte (c'est-à-dire quand il y
a deux étapes limitantes) où les deux étapes limitantes se produisent dans des zones réactionnelles de
même dimension. En effet dans ce cas de figure, d'après le théorème des lenteurs [12], il est possible
d'écrire la vitesse sous la forme :
- 15 -
dα 1
= (I.20)
dt  1 1 
 + 
 φ 1 E1 φ 2 E 2 
où φ1, φ2 sont les réactivités surfaciques de croissance de chacune des étapes considérées limitantes et
E1, E2 sont les fonctions d'espace de ces deux étapes. Sachant que ces deux étapes se produisent dans
des zones réactionnelles de même dimension, E1 est égale à E2. La vitesse est alors égale à :
dα  φ 1φ 2 
= E1 (I.21)
dt  φ 1 + φ 2 
La vitesse peut une nouvelle fois être considérée comme le produit de deux fonctions φ.E et le
test sera également vérifié dans ce cas de figure.
Ainsi, nous pouvons conclure, à l'aide du tableau I.1, que le test du φE sera vérifié
uniquement dans certains cas (sauf cas très particulier de diffusion sous champ électrique
[15]).
Cas pur, pseudo-stationnaire dα =φE

dt
Zones réactionnelles de dα =φE
mêmes dimensions dt
Cas mixte, pseudo-stationnaire
Zones réactionnelles de dα ≠φE
dimensions différentes dt
Autre cas : E n'est pas définissable dα ≠φE
dt
Tableau I.1 : Résultats du test du φE
En conclusion, il est possible de vérifier expérimentalement l'hypothèse de l'étape

limitante.
1 Choix d'une étape limitant la croissance
L'expression (I.18) pour la vitesse de croissance est celle que nous utiliserons dans la
suite car elle permet d'imaginer toutes les étapes limitantes de la croissance : diffusion ou
réaction d'interface. En effet, le tableau I.2 [9] montre quelles seraient les expressions de φ et
E dans les deux cas précédents.
- 16 -
φ E
Réaction élémentaire d'interface k0f(T,P) S/n0
Diffusion D(T)∆C/l0 GDSI /n0
Tableau I.2 : Expressions de φ et E selon l'étape limitant la croissance
Expliquons la signification des termes de ce tableau.
Si l'étape limitante est une réaction d'interface interne, c'est-à-dire à l'interface entre
les deux phases solides, on retrouve l'expression de Delmon pour la croissance (I.16) avec :
φ = ki (I.22)
Si
E= (I.23)
n0
On peut également imaginer l'étape limitant la croissance située à l'interface externe,
c'est-à-dire à la surface transformée des grains. La fonction d'espace est alors :
Se
E= (I.24)
n0
où Se représente l'aire de l'interface externe.
Si l'étape limitante est une réaction d'interface (interne ou externe), alors l'expression
de φ est donnée dans le tableau I.2 où k0 est le facteur de fréquence de l'étape élémentaire qui
limite la réaction et f(T,P) est une fonction qui dépend notamment du mécanisme réactionnel
de croissance.
Dans le cas où une étape de diffusion est l'étape limitant la croissance, on montre que
l'expression de la réactivité surfacique de croissance, en régime pseudo-stationnaire de
diffusion, est alors :
φ = D∆C (I.25)
l0
où D est le coefficient de diffusion de l'espèce diffusante au sein de la nouvelle phase formée,
∆C est la différence de concentrations de l'espèce diffusante entre l'interface interne et
l'interface externe et l0 est une longueur arbitraire. Ainsi, la réactivité surfacique de croissance
s'exprime également en mol.m-2.s-1 dans ce cas de figure.
L'expression de la fonction d'espace (en m2.mol-1) pour une diffusion limitant la
croissance est :
- 17 -
GD Si
E= (I.26)
n0
où GD est un terme sans dimension qui vient de l'expression du flux pour une géométrie
particulière. L'expression de celui-ci est donnée dans le tableau I.3 pour différentes
géométries de grains :
Géométrie GD
l0
Plaquette
x
l0
Cylindre
ri ln(re ri)
l0re
Sphère
ri(ri −re)
Tableau I.3 : Expression de GD selon la géométrie des grains
Dans le tableau I.3, x est la demi-épaisseur de la couche de produit formé pour une
plaquette, ri et re sont respectivement le rayon interne et le rayon externe de la couche de
solide produit pour une sphère ou un cylindre.
On pourra être amené à prendre pour l0 une valeur liée à la géométrie et à la taille du
grain comme par exemple, dans le cas de plaquettes, la demi-épaisseur initiale x0 ou dans le
cas de la sphère et du cylindre, le rayon initial r0.
Remarque : si la diffusion est l'étape limitante de la croissance alors celle-ci se déroule a

priori de manière anisotrope. En effet en supposant une croissance isotrope et une étape limitante de
diffusion, alors la vitesse à l'interface « triple » solide réactif, solide produit, phase gazeuse est infinie
puisque l'épaisseur de diffusion à cet endroit est nulle. Ainsi, à chaque fois que les trois phases sont en
contacts directs, on a une vitesse infinie. Ceci a pour conséquence la formation d'une couche de la
nouvelle phase autour du grain. La croissance est alors anisotrope.
En conclusion, quelle que soit l'étape limitant la croissance, diffusion ou réaction

d'interface, il est possible d'écrire la vitesse de croissance à l'aide de l'équation (I.18).
La loi décrivant la croissance de la nouvelle phase est déterminée avec les différents
éléments que nous venons de détailler (mode de croissance du germe, géométrie des grains,
développement de la nouvelle phase, étape limitante) dans le calcul de la fonction d'espace.
- 18 -
La croissance d'un germe va entraîner une diminution de la surface où de nouveaux

germes sont susceptibles de se former. Afin de décrire une réaction en cinétique hétérogène, il
paraît donc impossible de considérer séparément la germination et la croissance. Ainsi, les
modèles de transformation relatifs aux réactions, mettant en jeu des solides et décrivant la
fonction d'espace, seront obtenus en considérant les deux processus en même temps.
I.1.3 Modèles de transformation
L'évolution du degré d'avancement d'une réaction est la conséquence de l'évolution de

l'interface réactionnelle due à la germination et à la croissance. Deux échelles de temps sont
nécessaires pour décrire cette évolution : une échelle pour la germination (τ) et une échelle
pour la croissance (t). La loi de germination permet de connaître le nombre de germes qui
apparaissent par unité de temps, dN/dt. La vitesse de croissance caractérise l'évolution de
l'interface réactionnelle d'un germe né à l'instant τ. Elle est notée vc(t,τ). L'expression du
degré d'avancement α(t), proposée par plusieurs auteurs [1,5], est obtenue en considérant
l'évolution de l'ensemble des germes nés à différents instants τ et qui vont croître entre τ et t :
t
( t
)( )
α(t)=∫0σ ∫τ vc(t,x)dx dN
dt t =τ
dτ (I.27)
où σ est un facteur de forme du germe et (dN/dt)t=τ est une loi de germination du type (I.1) à
(I.13).
Certains auteurs [14,17,18] proposent d'utiliser la loi de germination (I.14). Ils

prennent ainsi en considération l’apparition de germes à une fréquence γ donnée sur la surface
« libre » de la poudre, c’est-à-dire la surface correspondant au solide initial. En prenant
l'équation (I.18) pour exprimer la vitesse de croissance, ils établissent alors une expression de
la vitesse sous la forme :
dα =φ(t) t[γ(τ)S (τ)G (t,τ)s (t,τ)]dτ
dt n0 ∫0 L d P (I.28)
où SL (τ) est la surface libre de germination à l'instant τ, sp(t,τ) est la surface à l'instant t de la
zone où se déroule l’étape limitant la croissance d'un germe né à l'instant τ et Gd est le terme
défini pour une étape limitant la croissance de diffusion. Dans le cas contraire (étape limitante
d'interface), celui-ci est égal à 1.
- 19 -
On retrouve ici une expression du type (I.18), et la vitesse peut s'écrire sous la forme
d'un produit de deux fonctions, dα/dt=φE, avec la fonction d'espace égale à :
E(t)= 1 ∫0[γ(τ)SL(τ)Gd(t,τ)sP(t,τ)]dτ
t
(I.29)
n0
Remarque : dans le cas où la vitesse de germination est très lente par rapport à la vitesse de
croissance, chaque fois qu'un germe apparaît sur un grain il recouvre complètement celui-ci. La vitesse
de la réaction s'exprime alors sous la forme :
dα n g γS L
= (I.30)
dt n0
où ng est la quantité de matière dans un grain. Il est donc possible d'écrire de nouveau la vitesse sous la
forme d’un produit de fonctions : dα/dt=ngγE avec E=SL/n0.
Selon les lois de germination et de croissance choisies, différents modèles ont été
développés par de nombreux auteurs depuis environ 70 ans. En 1966, Sharp et al. [19] ont
publié un article dans lequel ils répertorient neuf modèles de transformation présentés dans le
tableau I.4. Ils proposent d'écrire les lois de vitesse en utilisant deux fonctions, f(α) et g(α)
définies par :
dα
= kf ( α ) (I.31)
dt
α dα
g( α ) = ∫ (I.32)
0 f (α )
La définition de la constante cinétique k est différente selon les lois utilisées pour la
germination, pour la croissance et dépend des paramètres géométriques des grains. Dans la
dernière colonne de ce tableau, nous avons ajouté l'expression de la fonction d'espace selon
l'équation (I.18) pour les modèles cités.
Remarque : comme nous le verrons, il n'est pas toujours possible d'écrire l'équation (I.28) qui
exprime la vitesse de la réaction comme une fonction du degré d'avancement (équation I.31).
- 20 -
Symbole α)
f(α α)
g(α α)
E(α
A2 2(1−α)[−ln(1- α)] 1/2
[−ln(1- α)]1/2
A3 3(1−α)[−ln(1- α)]1/3 [−ln(1- α)]1/3
F1 (1- α) -ln(1- α) (s0/n0)(1- α)
R2 (1- α)1/2 2[1-(1- α)1/2] (s0/n0)(1- α)1/2
R3 (1- α)2/3 3[1-(1- α)1/3] (s0/n0)(1- α)2/3
D1 1/2α α2 (s0/n0)/2α
D2 1/−ln(1- α) (1- α)ln( 1- α)+α (2s0/n0)/−ln(1- α)
2/3
D3 (3(1−α) )/(2[1−(1- α) ])1/3
[1-(1- α)1/3]2
D4 3/(2[(1- α)-1/3-1]) 1−2α/3−(1- α)2/3 (s0/n0)/( [(1- α)-1/3-1])
(s0 étant la surface initiale des grains)
Tableau I.4 : Lois cinétiques répertoriées par Sharp et al. [19]
La signification de ces lois est décrite ci dessous.
1 Germination dans le volume et croissance isotrope (An)
Ce modèle a été développé par Avrami [8] et par Jonhson et Mehl [20]. Il a ensuite
été repris par Erofeev [21] avec une loi de germination différente. Ces auteurs considèrent une
germination dans le volume du solide suivie par une croissance isotrope des germes. Ce
modèle tient compte du recouvrement de sites potentiels à la germination et du recouvrement
des germes en cours de croissance. L'expression de la loi cinétique est de la forme :
− ln( 1 − α ) 1 n = kt (I.33)
Cette loi est notée An où n dépend de la forme des germes et de la loi de germination.
2 Germination aléatoire et croissance instantanée (F1)
Pour cette loi couramment appelée d'ordre 1 (et notée F1), Sharp et al. ne donnent pas
de signification physique. Or, si nous supposons que la germination aléatoire est suivie d'une
croissance instantanée sur chaque grain, il est aisé de démontrer que c'est précisément cette loi
qui en résulte :
ln(1 − α ) = −kt (I. 34)
Remarque : on retrouve cette loi en considérant le modèle d’Avrami pour des plaquettes et
une germination supposée instantanée c’est à dire l’équation (I.33) avec n=1.
- 21 -
3 Germination instantanée et croissance limitée par une étape d'interface interne

avec un développement interne de la nouvelle phase (Rn)
La germination a lieu sur tous les grains en même temps et chaque germe recouvre la
surface d'un grain. La croissance tangentielle est considérée comme très rapide par rapport à
la croissance radiale (croissance anisotrope). La réaction, qui est considérée limitante pour la
croissance, est située au niveau de l'interface réactionnelle entre la phase initiale et la nouvelle
phase formée. La loi cinétique s'écrit sous la forme :
1 − ( 1 − α )1 n = kt (I.35)
où n dépend de la géométrie de l'interface réactionnelle. Pour une plaquette, n=1, la loi est
appelée R1. Pour une interface de forme cylindrique, n=2, la loi est notée R2 alors que pour
une interface de forme sphérique, n=3, la loi est notée R3 et appelée couramment cœur
rétrécissant. On trouve la même expression dans le cas du cube rétrécissant.
4 Germination instantanée et croissance limitée par une étape de diffusion avec

un développement interne de la nouvelle phase (Dn)
Comme pour le cas précédent, la germination est instantanée sur tous les grains et
chaque germe occupe toute la surface d'un grain. La diffusion d'intermédiaires réactionnels au
sein de la couche de produit formé est l'étape limitante. Il est supposé que le rapport des
volumes molaires entre la phase produite et la phase initiale, noté Z (coefficient de Pilling-
Bedworth), est égal à 1.
La loi cinétique dépend de la géométrie du système. Pour une géométrie où la
diffusion a lieu selon une seule dimension (grains en plaquettes), la loi, notée D1 et appelée loi
parabolique, s'écrit sous la forme :
α 2 = kt (I.36)
Pour une diffusion se déroulant dans des grains cylindriques, l'équation cinétique notée
D2 a été obtenue par Holt et al. [22] :
(1 − α ) ln(1 − α ) + α = kt (I.37)
Cette expression pourrait également être calculée dans le cas où le coefficient de
Pilling-Bedworth, Z, est différent de 1.
- 22 -
Enfin pour une diffusion à l'intérieur de grains sphériques, la loi cinétique déterminée
par Ginstling et Brounshtein [23] est notée D4 et s'exprime sous la forme :
1 − 2α / 3 − (1 − α ) = kt
23
(I.38)
L'expression de la loi de diffusion dans des grains sphériques a également été calculée
par Valensi [24] et Carter [25] dans le cas où le volume molaire de la phase finale est
différent du volume molaire de la phase initiale (Z≠1). Ceci conduit à une variation du
volume des grains. L'équation cinétique s'écrit :
Z − [1 + ( Z − 1 )α ] − ( Z − 1 )( 1 − α ) 2 3
23
= kt (I.39)
2( Z − 1 )
Une dernière loi de diffusion est fréquemment utilisée dans la littérature. Il s'agit du
modèle de Jander [26], noté D3 et qui s'exprime par l'équation :
[1 − (1 − α ) ] 13 2
= kt (I.40)
Cependant comme le soulignent plusieurs auteurs [27,28], cette loi n'a aucune
signification physique. Elle est l'application directe de la loi parabolique, uniquement valable
en symétrie plane, à une géométrie sphérique. Elle ne peut être envisagée que pour des
couches de faible épaisseur car en début de réaction, les expressions des flux de diffusion en
coordonnées planes et sphériques sont sensiblement les mêmes.
Excepté les lois d’Avrami (An), les lois du tableau I.4 correspondent à des modèles
dans lesquels, soit la germination (Rn, Dn), soit la croissance (F1) est supposée instantanée.
Les lois d’Avrami étant établies sur la base de l’apparition aléatoire de germes en volume,
elles ne sont en principe pas utilisables dans le cas des réactions entre un solide et un gaz ou
des décompositions de solides inorganiques pour lesquelles les germes se forment à la surface
du solide. Cependant, sans doute parce que les lois d’Avrami sont simples et donnent des
courbes d’allure sigmoïde, force est de constater qu’elles ont été fréquemment choisies pour
rendre compte des courbes expérimentales quel que soit le type de transformation (et donc de
germination).
Il existe cependant dans la littérature un modèle, développé par Mampel [29] en 1940
et par ailleurs par Johnson et Mehl [20] en 1939, qui est bien adapté à la germination en
surface suivie d'une croissance isotrope. Peu utilisé dans le passé, probablement à cause de
l’importance des calculs mis en jeu par rapport aux moyens de l’époque, le modèle de
- 23 -
Mampel a été repris par Delmon (pour des grains sphériques et des plaquettes) [2]. L’intérêt
de ce modèle est qu’il permet de remonter à deux grandeurs caractéristiques de la germination
et de la croissance, notées respectivement k'g0 et ki par Delmon : k'g0 correspond à la
fréquence surfacique de germination γ (nombre de germes.s-1.m-2) et ki est la vitesse
spécifique de croissance et sera remplacé avec, l’expression de vitesse donnée par l’équation
(I.18), par φ.
Avec le modèle de Mampel, en faisant varier un seul paramètre de modèle (noté A),
qui est proportionnel au rapport γ/φ (quelle que soit la forme géométrique des grains), il est
possible de rendre compte des courbes expérimentales [2]. Des exemples de courbes de
vitesses réduites ω (vitesse divisée par vitesse maximum) en fonction du degré d'avancement
sont représentés pour différentes valeurs du paramètre de modèle sur la figure I.3.
Figure I.3 : Courbes de vitesse en fonction du degré d'avancement pour le modèle de Mampel et pour
différentes valeurs du paramètre de modèle A [2]
Grâce aux moyens de calculs actuels, il est assez facile de calculer la vitesse selon le
modèle de Mampel. Viricelle [18] et Bouineau [14] ont pu interpréter leurs courbes de vitesse
grâce à ce modèle pour la transformation de l’hydroxycarbonate de cérium en cérine et la
décomposition du carbonate de calcium, respectivement.
De manière analogue à la mise en équation du modèle de Mampel à l’aide d’un

paramètre A, et en considérant que la croissance des germes est cette fois-ci anisotrope (cf.
figure I.1), divers modèles de germination-croissance anisotrope ont été établis sur la base de
l’équation (I.28). En effet, en choisissant un sens de développement de la nouvelle phase, une
- 24 -
étape limitante pour la croissance et une forme de grains, les expressions de GD, s0 et SL sont
déterminées. Ainsi, Bouineau [14] a interprété les courbes de vitesse obtenues pour la
déshydratation du sulfate de lithium monohydraté, et Brun [17] celles de la réduction de U3O8
par l’hydrogène. Nous reviendrons plus tard sur cette réaction.
Comme pour le calcul de la vitesse d’après le modèle de Mampel, avec les modèles de
germination-croissance anisotrope, il est possible de paramétrer les courbes calculées à l’aide
d’un seul paramètre de modèle, également noté A, défini de la même manière. La
comparaison des courbes expérimentales aux lois de vitesse à l'aide du paramètre de modèle
est présentée dans l'annexe D.
Remarque : des modèles de germination-croissance anisotrope [30] ont été proposés dans le
passé à partir de l’équation (I.27). Nous ne les présentons pas ici car le choix d’une loi de germination
en surface indépendante de la quantité de surface restant disponible pour la formation de nouveaux
germes ne nous a pas semblé très réaliste.
1 Modèles de transformation
Finalement on s'aperçoit que l'ensemble des modèles pouvant être développés afin de
décrire une transformation prend en compte :
• le mode de germination : instantanée ou non ;

• le mode de croissance du germe : instantanée, anisotrope ou isotrope ;
• le sens du développement de la nouvelle phase : interne ou externe ;
• l'étape limitant la croissance : réaction interfaciale (interne ou externe) ou diffusion ;
• la géométrie des grains : sphère, cylindre, plaquette.
Le tableau I.5 rassemble à notre connaissance la plupart des modèles ainsi déterminés.
Les modèles correspondant aux lois répertoriées par Sharp sont précisés dans le
tableau I.5. De plus, dans ce tableau figurent également des références d'exemples de
réactions correspondant à d'autres modèles. On peut remarquer que c'est à partir de l'équation
(I.18) de la vitesse que certains modèles ont été construits sans qu'il n'existe (à notre
connaissance) de cas concret correspondant.
- 25 -
Grains sphériques Grains cylindriques Grains plaquettes

Croissance instantanée (F1) (F1) (F1)
Interface interne (R3) Interface interne (R2) Interface interne
Développement
Interface externe Interface externe Interface externe
interne
Germination Diffusion (D4) Diffusion (D2) Diffusion (D1)
instantanée Interface interne Interface interne Interface interne
Développement
Interface externe Interface externe [11] Interface externe
externe
Diffusion Diffusion Diffusion
Interface interne [31] Interface interne Interface interne
Développement
Interface externe [32] Interface externe Interface externe
interne
Germination- Diffusion Diffusion Diffusion
croissance
anisotrope Interface interne Interface interne Interface interne
Développement
Interface externe Interface externe Interface externe
externe
Diffusion Diffusion Diffusion
Développement Interface interne [33,34] Interface interne Interface interne

Germination- interne Interface externe Interface externe Interface externe
croissance
isotrope Développement Interface interne Interface interne Interface interne
externe Interface externe [35] Interface externe [35] Interface externe
Tableau I.5 : Modèles de transformation
Remarque : dans le cas d'une étape limitante d'interface externe ou de diffusion, il est
également nécessaire de tenir compte dans le calcul de la vitesse du coefficient de Pilling-Bedworth
(ou coefficient d'expansion volumique) Z.
Ce tableau met en évidence le grand nombre de modèles possibles pour décrire une
réaction en cinétique hétérogène. Chaque modèle entraîne un calcul spécifique afin de
déterminer l'expression de la fonction d'espace afin d’obtenir la vitesse de la réaction.
Les modèles décrits précédemment peuvent être divisés en deux groupes :

• les modèles limites où un des deux phénomènes, germination ou croissance, est
considéré comme se déroulant très rapidement par rapport à l'autre ;
• les modèles complexes de germination-croissance.
- 26 -
Dans le premier cas, le degré d'avancement caractérise complètement la vitesse de la

réaction pour des conditions physico-chimiques données (la fonction d'espace est une fonction
de α).
Plusieurs études ont montré qu'il est possible de vérifier expérimentalement dans quel
groupe de modèles se situe le système étudié [9,14]. Comme le présente la figure I.4, dans le
cas de la décomposition du carbonate de calcium [14], deux expériences sont réalisées avec au
début de la réaction des conditions physico-chimiques différentes. Pour les expériences de la
figure I.4, la pression de dioxyde de carbone est différente. Puis à un degré d'avancement
donné, on place ces deux expériences dans les mêmes conditions (ici PCO2=800Pa). Si à partir
de ce même degré d'avancement, les expériences ont la même vitesse alors la réaction peut
être décrite par un modèle limite. En effet dans ce cas, le passé n'influence pas la vitesse de la
réaction. Dans le cas contraire, comme pour la décarbonatation du carbonate de calcium, la
réaction suit un modèle complexe de germination-croissance.
Figure I.4 : test du cas limite appliqué à la décomposition du carbonate de calcium [14]
Soustelle et Pijolat [9] ont appelé ce test le test du cas limite. Flynn [36] propose une
méthode similaire afin de vérifier l'influence de l'histoire thermique sur la vitesse de réaction.
Remarque : il est également possible de savoir si le modèle correspond à un cas limite selon
la forme des courbes expérimentales. En effet, la présence d'un maximum sur les courbes de vitesse est
en général caractéristique d'une « compétition » entre la germination et la croissance et donc d'un
modèle où les deux phénomènes jouent un rôle sur la vitesse.
- 27 -
Nous avons détaillé, au sein de ce paragraphe, les modèles de transformation et leur

signification. Nous avons pu mettre en avant deux approches concernant les lois de vitesses de
ces modèles. La première se limite aux lois de Sharp alors que la deuxième prend en compte
l'ensemble des modèles de germination-croissance.
Les équations de vitesse liées à ces deux approches (I.31) et (I.28) dépendent de
grandeurs cinétiques, k ou φ selon les notations utilisées. Il paraît nécessaire de connaître
l'influence des conditions expérimentales sur ces grandeurs afin de comprendre comment la
vitesse de la réaction évolue lorsqu'on modifie celles-ci.
I.1.4 Influence de la pression et de la température
La vitesse d'une réaction peut dépendre de plusieurs paramètres : température,

pressions partielles des différents gaz ...
Afin de comprendre comment ces paramètres jouent un rôle sur la vitesse d'une
réaction, deux approches différentes ont été également développées.
La première approche, utilisée par un grand nombre d'auteurs, suppose de pouvoir

écrire la vitesse sous la forme :
dα
= f 1 ( T ) f 2 ( P1 ).... f ( α ) (I.41)
dt
Cette équation est à variables séparées T,P1,...α. Comme nous l’avons déjà évoqué,
elle n'est notamment pas utilisable pour décrire les cas complexes de germination croissance.
La deuxième approche ne suppose pas d'hypothèse sur l'expression de la vitesse en

fonction des paramètres intensifs. Elle revient à utiliser l'équation (I.18),
α/dt=φ(
dα φ(T,P
φ( 1...)E(t), et à déterminer l'influence de ces paramètres sur la réactivité surfacique
de croissance.
Remarque : pour les travaux menés en considérant la séparation des variables, la constante
cinétique est toujours notée k et dépend seulement de la température.
- 28 -
Pour les études prenant en compte la non-séparation des variables, les notations diffèrent selon
les auteurs. Afin de rendre plus compréhensible cette partie de l'étude bibliographique, nous utiliserons
l'équation (I.18). La grandeur cinétique sera à chaque fois la réactivité surfacique de croissance φ. Ces
études supposent une succession d'étapes élémentaires pour décrire les mécanismes réactionnels.
Ainsi, nous noterons ki la constante de vitesse de l'étape élémentaire i. Les constantes ki dépendent de
la température selon la loi d'Arrhenius [12].
I.1.4.1 Approche utilisant la séparation des variables
L'expression de la vitesse dans ce cas est sous la forme de l'équation (I.41). Cette
équation dépend du type de réaction.
Pour les décompositions de solides, la vitesse s'écrit en général :

dα
= k( T ) f ( α ) (I.42)
dt
Pour les réactions solide-gaz, la vitesse de croissance peut dépendre ici de la pression
du gaz réactif, P1 et de la pression du gaz produit, P2. Comme nous le verrons dans le
chapitre IV avec l'exemple de la réduction d'oxyde d'uranium par l'hydrogène, de nombreux
auteurs se contentent de chercher une loi à ordre vis-à-vis de la pression P1. Ainsi, ils écrivent
la vitesse de réaction sous la forme :
dα
= k( T )P1n f ( α ) (I.43)
dt
Etant donné que les variables sont séparées, il apparaît possible d'étudier
indépendamment l'influence de la température et de la pression en comparant les vitesses en
gardant deux variables constantes (par exemple en étudiant l'influence de la température à
pression et degré d'avancement constants).
- 29 -
I.1.4.1.1 Influence de la pression
Pour le gaz réactif, les auteurs ne justifient pas le choix de rendre compte de
l'influence de la pression à l'aide d'une loi à ordre. Ils n'expliquent pas la signification
physique de cet ordre n.
Pour les deux équations (I.42) et (I.43), la pression du gaz produit n'est pas prise en
compte. Or ceci ne peut être vrai que si la réaction se déroule loin des conditions d'équilibre
thermodynamique. Mais comme le souligne Bertrand et al. [37], si la pression est proche de
la pression d'équilibre, il est nécessaire de prendre en compte dans l'équation cinétique l'écart
à l'équilibre, c'est-à-dire l'écart entre la pression de la réaction et la pression d'équilibre. Cette
dernière dépend de la température. Dans ce cas de figure, les variables température et pression
ne sont plus séparées et cette approche paraît donc ne plus convenir.
I.1.4.1.2 Influence de la température
La plupart des études décrivent la variation de la constante cinétique en fonction de la

température selon la loi d'Arrhenius sous la forme :
k( T ) = A exp( − Ea RT ) (I.44)
où A est appelé facteur pré-exponentiel ou facteur de fréquence et Ea est l'énergie d'activation.
Celle-ci est souvent définie comme la barrière d'énergie (sans préciser ce qui doit franchir
cette barrière) qu'il faut surmonter pour que la réaction se déroule. Le but de nombreuses
études en cinétique hétérogène est la détermination du « triplet cinétique » A, Ea, f(α).
Dans les études cinétiques, la valeur de l'énergie d'activation varie avec les conditions
expérimentales. Un exemple est la décomposition du carbonate de calcium où selon les
auteurs, l'énergie d'activation varie de 146 kJ.mol-1 à 960 kJ.mol-1 [38]. Il paraît alors
nécessaire de remettre en cause l'équation (I.44) pour exprimer l'influence de la température
sur la constante cinétique puisqu'elle conduit à des résultats si disperses. De plus certaines
études conduisent également à une variation de l'énergie d'activation en fonction du degré
d'avancement de la réaction. Comme le souligne Garn [39], si la valeur de l'énergie
d'activation varie avec les conditions expérimentales, la loi d'Arrhenius n'est pas applicable
directement. Dans ce cas de figure, l'auteur suppose que la loi d'Arrhenius n'a aucun sens ou
- 30 -
que certains paramètres n'ont pas été pris en compte dans la vitesse alors qu'ils jouent un rôle
sur celle-ci.
Bien que selon les études réalisées sur une même réaction et selon les conditions
expérimentales, les valeurs de l’énergie d’activation et du facteur pré-exponentiel ne soient
pas les mêmes, certains auteurs observent une relation entre les valeurs de l'énergie
d'activation et les valeurs du facteur pré-exponentiel sous la forme :
ln A = a + bEa (I.45)
Cette relation entre A et Ea est appelée effet de compensation.
Afin de résumer cette approche sur la séparation des variables, prenons un exemple.
Dans le cas d'une réaction limitée par une étape interfaciale, et où on suppose l'adsorption
gazeuse (décrite par la théorie de Langmuir) à l'équilibre, la vitesse peut suivre une loi
homographique en fonction de la pression et s'écrire :
a( T )P
k( T , P ) = (I.46)
1 + b( T )P
Il apparaît clairement que, dans ce cas, les variables température et pression ne sont
pas séparées.
De plus, la séparation de la variable α dans l'équation (I.41) restreint cette approche
aux modèles limites.
Afin de pouvoir palier à ces problèmes, une deuxième approche a été proposée sans
supposer la séparation des variables.
I.1.4.2 Approche sans supposer la séparation des variables
Dans cette approche aucune hypothèse n'est réalisée sur la loi liant la grandeur
cinétique (φ) aux paramètres intensifs de la réaction. Les études réalisées ont essayé de
montrer comment ces différents paramètres pouvaient jouer sur la réactivité surfacique de
croissance.
- 31 -
I.1.4.2.1 Détermination de la réactivité surfacique de croissance en

fonction de la température et des pressions partielles des gaz
Delmon [2] puis Barret [3] ont développé une méthode afin d'accéder directement à
l'évolution de la réactivité surfacique de croissance en fonction de la température ou de la
pression partielle des gaz de la réaction. Cette méthode est connue sous le nom de méthode de
l'isolement.
Elle consiste à commencer une série d'expériences avec les mêmes contraintes
physico-chimiques (T1, P1,...) jusqu'à un même instant td. A cet instant, une des contraintes
(celle dont on veut évaluer l'influence) est brusquement modifiée et amenée rapidement à une
nouvelle valeur. Comme le montre la figure I.5, dans le cas de la réduction de l'oxyde de
nickel par l'hydrogène, étude réalisée par Delmon [2] avec des variations de température, la
valeur de la contrainte modifiée (ici la température) est différente après le décrochement pour
chaque expérience. Il réalise des décrochements en température entre 206,5°C et différentes
températures (Ti).
Figure I.5 : Application de la méthode de l'isolement à la réduction de l'oxyde de nickel par

l'hydrogène, influence de la température [2]
A chaque expérience, sont mesurées les vitesses avant ((dα/dt)av) et après le

décrochement ((dα/dt)ap). Le rapport de ces vitesses, dans le cas d’une variation de
température, est égal à :
( dα dt )ap ( t d ) φ ( Ti , P ,...)E( t d )
= (I.47)
( dα dt )av ( t d ) φ ( T1 , P ,...)E( t d )
- 32 -
Cette variation doit être rapide afin que la fonction d'espace reste constante durant le
temps du décrochement. Le rapport des vitesses se simplifie pour obtenir le rapport des
réactivités de croissance :
( dα dt )ap ( t d ) φ ( Ti , P ,...)
= (I.48)
( dα dt )av ( t d ) φ ( T1 , P ,...)
Ainsi, il est possible de déterminer l'évolution de φ(T) directement à partir
d'expériences. En fait ce n'est pas exactement les valeurs de φ(T) qui sont obtenues mais les
rapports φ(T)/ φ(T1). De même, en réalisant des décrochements en pression, il est possible
d'obtenir l'évolution de la réactivité surfacique de croissance en fonction de la pression
partielle d'un gaz à une constante près.
Cette méthode a été reprise par Viricelle et al. [16] pour la transformation de
l'hydroxycarbonate de cérium en cérine. Ainsi, la figure I.6 représente d'une part les courbes
du degré d’avancement (noté ici λ) en fonction du temps pour des décrochements en pression
de dioxyde de carbone et d'autre part l'évolution de la vitesse mesurée après le décrochement,
qui est proportionnelle à la réactivité surfacique de croissance, en fonction de cette pression.
Ils ont pu ainsi mettre en évidence l'influence de la pression de dioxyde de carbone sur la
réactivité surfacique de croissance de cette réaction.
Figure I.6 : Détermination de φ(PCO2) par la méthode des décrochements pour la transformation de
l'hydroxycarbonate de cérium en cérine [16]
Ils ont réalisé un travail similaire pour étudier l'influence de la pression d'oxygène sur
cette réaction.
- 33 -
Bouineau [14] a également utilisé cette méthode pour déterminer l'influence de la

pression de dioxyde de carbone sur la décarbonatation du carbonate de calcium.
Cette méthode, appelée aussi méthode des décrochements, est très intéressante car elle
peut être réalisée sans connaître le modèle de transformation de la réaction et permet donc
d'accéder directement par l'expérience à la connaissance de φ(T) et φ(P) à une constante près.
Remarque 1 : comme l'a montré Surla et al. [11] pour l'oxydation du magnésium, un autre
intérêt de cette méthode est la possibilité d'obtenir l'expression de la réactivité surfacique de croissance
en fonction de la température et de la pression même avec des expériences non reproductibles dues à la
non-reproductibilité de la fonction d'espace.
Remarque 2 : il n’est pas possible de connaître l’évolution de la fréquence surfacique de

germination, γ, en fonction des contraintes intensives de la réaction en utilisant cette méthode.
Remarque 3 : si le modèle géométrique et l'équation de la vitesse sont connus, il est

également possible de connaître l'évolution de la réactivité surfacique de croissance et de la fréquence
surfacique de germination en fonction de la température et de la pression en réalisant une série
d'expériences en conditions isotherme et isobare. L'ajustement des courbes calculées par rapport aux
courbes expérimentales permet de déterminer les valeurs de la fréquence surfacique de germination et
de la réactivité surfacique de croissance dans les conditions expérimentales où ont été réalisées les
expériences. On obtient ainsi les courbes φ(T), φ(P), γ(T) et γ(P), si un nombre suffisamment grand
d'expériences a été effectué avec différentes températures et différentes pressions partielles du gaz
étudié. Comme le montre l'annexe D, la détermination de ces grandeurs est facilitée par le passage en
vitesse réduite (vitesse divisée par vitesse à demi réaction) et la connaissance du paramètre de modèle
A (cf. I.3).
I.1.4.2.2 Expressions de φ(T,P)
La croissance est la conséquence d'une succession d'étapes élémentaires qui se

déroulent dans des zones réactionnelles différentes. Il est possible d'imaginer un mécanisme
réactionnel qui représente cet ensemble d'étapes.
- 34 -
En supposant qu'une étape est très lente par rapport aux autres, celle-ci impose sa
vitesse à la réaction. Les autres étapes élémentaires sont alors considérées comme étant à
l'équilibre. La recherche de l'étape limitante va permettre d'interpréter les courbes φ(T,P).
Dans la suite de ce chapitre nous développons quelques études qui illustrent la manière
de déterminer φ(T,P).
Pour les décompositions thermiques, Bouineau [14] classe les courbes φ(P) en deux
catégories :
• celles qui sont monotones où φ diminue avec la pression et qui sont représentées sur la
figure I.7 ;
• celles qui présentent un ou plusieurs extremum représentées sur la figure I.8.
a b
Figure I.7 : Influence de la pression, courbes décroissantes
d
c
Figure I.8 : Influence de la pression, courbes à extremum
Dans le cas où la courbe φ(P) admet un ou plusieurs extremum, elles sont

représentatives d'un effet Smith-Toppley. Ces auteurs sont les premiers à avoir mis en
évidence ce phénomène pour la déshydratation de l'oxalate de manganèse hydraté [40]. Cet
effet est souvent observé pour les réactions de déshydratation.
- 35 -
Afin d'interpréter les courbes constamment décroissantes, Barret [3] propose un

mécanisme de croissance en trois étapes pour la décomposition thermique d'oxydes
métalliques :
• 1) passage d'un ion oxygène du réseau en position d'oxygène adsorbé ;
• 2) désorption de l'oxygène ;
• 3) réorganisation du réseau pour former la nouvelle phase.
Il suppose le système quasi-stationnaire (ce qui revient en fait à la pseudo-stationnarité
que nous avons définie). L'adsorption est par hypothèse conforme à l'isotherme de Langmuir.
Il suppose également que la dernière étape n'est pas limitante.
Comme le présente le tableau I.6, selon l'étape limitante choisie pour le mécanisme, il
obtient des lois de vitesse qui décroissent linéairement avec la pression (figure I.7.a) ou des
lois de vitesse homographiques en pression (figure I.7.b).
Etape limitante Vitesse de réaction

Réaction interfaciale Peq − P
(passage de l'ion oxygène du réseau dans la couche de φ1 = k1 (figure I.7.b)
gaz adsorbé) Peq ( 1 + K a P )
Peq − P
Désorption φ2 = k2 (figure I.7.a)
1 + K a Peq
Tableau I.6 : Expression de la vitesse de réaction selon l'étape limitante pour la décomposition
thermique d'oxydes métalliques [3]
Dans le tableau I.6, Ka représente la constante d'équilibre d'adsorption, Peq est la

pression partielle d'oxygène à l'équilibre et P la pression de l'oxygène de la réaction. Les
grandeurs Ka et Peq dépendant de la température, les variables température et pression ne sont
pas séparées.
Dans son étude sur la décarbonatation du carbonate de calcium et la déshydratation du

sulfate de lithium monohydraté, Bouineau [14] a établi le même mécanisme de croissance en
trois étapes pour chacune des réactions :
• 1) passage d'une molécule du solide à l'état adsorbé ;
• 2) désorption du gaz produit ;
• 3) changement de phase et avancement de l'interface.
- 36 -
Selon l'étape limitante choisie, il obtient plusieurs lois exprimant φ(T,P). En

comparant aux valeurs expérimentales de φ, il trouve pour les deux réactions la même étape
limitante, l'étape 1, et obtient l'expression suivante :
K2  P
φ = k1  1 −  (I.49)
(K 2 + P )  Pe 
où K2 est la constante d'équilibre de l'étape 2, P est la pression du gaz étudié et Pe est la
pression d'équilibre de ce même gaz à la température de la réaction.
Pour les réactions solide-gaz, la réactivité surfacique de croissance peut dépendre de

la pression du gaz réactif, P1 et de la pression du gaz produit, P2. L'étude de l'influence des
paramètres intensifs doit permettre de déterminer φ(T,P1,P2).
De nouveau Barret [3] s'appuie sur l'écriture d'un mécanisme de croissance relatif à la
réduction d'oxyde métallique par l'hydrogène pour rendre compte de l'influence des différents
paramètres sur cette réaction. Il considère une nouvelle fois le système quasi-stationnaire (ou
plutôt pseudo-stationnaire) et suppose un mécanisme en quatre étapes :
• 1) adsorption de l'hydrogène ;
• 2) formation d'une molécule d'eau adsorbée à partir d'un anion superficiel ;
• 3) désorption de l'eau ;
• 4) réorganisation du réseau du métal et renouvellement de l'interface.
Cette dernière étape est considérée comme étant à l'équilibre et n'est donc pas l'étape
limitante.
Remarque : dans les deux mécanismes écrits par Barret, décomposition d'oxydes métalliques
et réduction d'oxydes métalliques par l'hydrogène, celui-ci suppose l'absence de couche protectrice
c'est-à-dire que le solide produit ne constitue pas une couche imperméable au gaz. Ainsi, il ne fait pas
intervenir, dans ses mécanismes réactionnels, d'étapes de diffusion au sein de la nouvelle phase
formée.
L'influence de la pression et de la température dépend de l'étape limitante. Les

différentes lois trouvées, selon l'étape choisie comme limitante, sont rassemblées dans le
tableau I.7 avec Ki représentant la constante d'équilibre de l'étape élémentaire i et dépendant
de la loi de Van't Hoff, K0 égale au produit de ces constantes.
- 37 -
Etape limitante Vitesse de réaction

1 K 0 P1 − P2
Adsorption d'hydrogène φ1 = k 1
K 1 P2 ( 1 + K 2 ) + K 2 K 3
K 0 P1 − P2
Formation d'une molécule d'eau adsorbée φ2 = k2
K 2 P2 + P1 K 0 + K 2 K 3
1 K 0 P1 − P2
Désorption de la vapeur d'eau φ3 = k3
K 3 P1 K 1 K 4 ( 1 + K 2 ) + 1
Tableau I.7 : Expression de la vitesse de réaction selon l'étape limitante pour la réduction d'oxyde
métallique par l'hydrogène [3]
Brun [17] a également développé un mécanisme de croissance pour la réduction de

l'octooxyde de triuranium par l'hydrogène. Elle imagine quatre étapes élémentaires :
adsorption de l'hydrogène, désorption de l'eau, diffusion des défauts au sein de la nouvelle
phase formée et réorganisation pour former la nouvelle phase.
En considérant que l'étape limitante est la désorption de l'eau, elle obtient l'expression
de la réactivité surfacique de croissance suivante :
1
k 2 K 1 K 4 2 PH 2  
φ=  1 − PH 2O / PH 2  (I.50)
2  1 
 1 1   K1K 2 K4 2 
 1 + K 1 2 K 4 8 PH 2 
 
Ce mécanisme sera présenté plus en détail dans le quatrième chapitre.
Comme on le voit à travers ces exemples, l'influence de la pression et de la

température dépend du mécanisme réactionnel et de l'étape limitant la croissance. Les
variables température et pression ne sont pas séparées puisque les constantes cinétiques (ki) et
thermodynamiques (Ki) dépendent de la température.
Remarque : afin de rendre compte de l'évolution de la grandeur cinétique caractéristique de la

germination (la fréquence surfacique de germination, γ), Brun [17] et Bouineau [14] ont également
développé des mécanismes réactionnels en plusieurs étapes élémentaires pour la germination. Ils
obtiennent ainsi, comme pour la croissance, les fonctions γ(T,P) qui dépendent de l'étape limitante
choisie pour la germination.
Dans les études que nous allons développer dans ce manuscrit, nous décrirons la
croissance à l'aide d'une succession d'étapes élémentaires. Pour prendre en compte l'influence
de la pression et de la température, il sera nécessaire de déterminer l'étape limitante du
- 38 -
mécanisme réactionnel. On remarquera que la détermination de la zone où se déroule cette

étape est également nécessaire pour la détermination d'un modèle de transformation. Nous
chercherons également à déterminer un mécanisme réactionnel et une étape limitante pour la
germination.
L'écriture d'un mécanisme de croissance et d'un mécanisme de germination et la

détermination d'une étape limitante pour chacun de ces mécanismes sera appelée
modélisation physico-chimique.
I.1.5 Effet d'une distribution granulométrique
Les différents modèles de transformation décrits dans le paragraphe I.1.3 sont

développés pour des poudres constituées de grains de même forme et de même taille. Afin de
se rapprocher des solides réels, certains auteurs ont travaillé sur l'effet engendré par une
distribution de tailles de grains sur les courbes cinétiques. Les fonctions mathématiques
représentant les distributions granulométriques sont en général des lois normales ou log-
normales en fonction du rayon des grains :
1 ( r − µ )2
g( r ) = exp( − ) (loi normale) (I.51)
σ 2π 2σ 2
1 (ln r − ln µ )2
g( r ) = exp( − ) (loi log-normale) (I.52)
σ 2π r 2σ 2
où µ est le rayon moyen, σ est la déviation standard qui caractérise l'étalement de la
granulométrie et r est le rayon des grains de la poudre ; g(r) représente alors la proportion de
grain dont le rayon est compris entre r et r+dr.
La prise en compte de cette distribution granulométrique dans la modélisation se fait
en calculant le degré d'avancement des grains de chaque taille en fonction du temps, α(t,r), et
en intégrant sur l'ensemble des tailles des grains pour obtenir le degré d'avancement global de
la réaction, α(t) :
∞
α( t ) = ∫0 α( t , r )g( r )d ( r ) (I.53)
Cette étude bibliographique montre l'influence de la granulométrie pour des modèles

particuliers.
- 39 -
Delmon [2] a étudié l'influence de la granulométrie de la poudre dans le cas d'un cube
rétrécissant. Il montre que les variations de µ et σ ont une influence sur la forme des courbes
cinétiques.
McIlvried et Massoth [28] ont étudié l'influence d'une distribution granulométrique

suivant les lois normales et log-normales pour les modèles R3 et D4. Ils déduisent de la forme
des courbes obtenues que seulement dans le cas où la loi granulométrique est log-normale et
le modèle diffusionnel (D4), la distribution granulométrique a un effet réel sur les courbes
cinétiques. Ils montrent également que la forme de la courbe cinétique obtenue par l'équation
de Jander (D3), qui n'a pas de signification physique, peut être obtenue à partir de la loi D4 et
d'une distribution granulométrique. Ils expliquent ainsi pourquoi le mécanisme D3 est souvent
trouvé comme loi décrivant le mieux certaines réactions.
Miyokawa et Masuda [41] ont étudié l'effet d'une distribution granulométrique sur
trois modèles : cube contractant, modèle de Jander et modèle d'Avrami. Ils regardent
l'influence de la distribution granulométrique sur les valeurs de l'énergie d'activation et du
facteur pré-exponentiel et constatent que ces grandeurs varient avec la distribution
granulométrique seulement pour le modèle d'Avrami. Lahiri [42] réalise une étude similaire
pour un modèle de cœur rétrécissant (R3) en condition non isotherme. Il observe que Ea et A
varient avec la distribution granulométrique et propose d'interpréter l'effet de compensation
(lnA=a+bEa) comme provenant des effets de la distribution granulométrique.
Koga et Criado [43,44] ont réalisé l'étude du rôle d'une distribution de taille de grains
sur les modèle D1, D2, D3, D4, R1, R2 et R3. Ils constatent que la détermination de l'énergie
d'activation ne dépend pas de la distribution granulométrique par contre celle-ci a une
influence sur les valeurs du facteur pré-exponentiel et sur la forme des courbes cinétiques.
Ainsi, la figure I.9 où sont tracées, pour les modèles R3 et D4, les courbes cinétiques de α en
fonction du temps réduit (temps divisé par le temps à demi réaction) pour différentes valeurs
de σ, montre que la distribution granulométrique a une influence sur la fin de courbes. Ils
observent également que plus l'étalement de la granulométrie est grand (σ élevée), plus il est
difficile de faire le choix entre les modèles de transformation. Johnson et Kessler [45]
obtiennent le même type de courbe (figure I.9) pour le modèle D4 au cours de leur étude sur
l'influence d'une distribution granulométrique sur la déshydroxylation de la kaolinite.
- 40 -
Figure I.9 : Courbes cinétiques R3 et D4 avec différentes distributions granulométriques [43]
Hutchinson et al. [46] ont étudié la décomposition de PbN6 avec des poudres de
granulométries différentes. En choisissant un modèle de germination-croissance où la
croissance est du type cœur rétrécissant (R3), ils arrivent à bien interpréter les courbes
cinétiques expérimentales en tenant compte des différentes distributions granulométriques des
poudres étudiées.
A travers cette étude bibliographique, il apparaît que la distribution granulométrique

d'une poudre a souvent un effet non négligeable sur la courbe cinétique d'une réaction. Par
conséquent, nous essaierons d'estimer le rôle d'une distribution de taille de grains sur la
cinétique des réactions que nous allons étudier par la suite.
Après ce rappel sur les modèles relatifs aux réactions solide-gaz, voyons à présent
comment ceux-ci sont utilisés en condition non isotherme.
I.2 Modélisation en condition non isotherme
I.2.1 Objectifs
La modélisation en condition non isotherme a pour but de rendre compte des modèles
de transformation décrits précédemment lorsque la température varie au cours de la réaction.
- 41 -
Les premières études réalisées en programmation de température datent des années

1960. Les avantages de faire varier la température au cours de la réaction sont nombreux
comme le souligne Carrasco [47].
Tout d'abord ces expériences permettent d'étudier et de modéliser la réaction sur un
large domaine de température. Une expérience suffit pour connaître l'évolution des grandeurs
cinétiques en fonction de la température alors qu'en condition isotherme, il fallait réaliser
toute une série d'expériences à différentes températures pour balayer le même domaine. Ainsi,
la modélisation et l'étude en condition non isotherme apparaissent plus rapides.
De plus, étant donné que moins d'expériences sont nécessaires, les quantités de produit
nécessaires sont également moins importantes. Cela présente surtout un avantage dans le cas
de réactifs coûteux ou difficiles à obtenir.
Ensuite dans le cas d'échantillons présentant de grandes hétérogénéités, il est évident
que l'utilisation d'un seul échantillon pour étudier tout un domaine de température simplifiera
l'étude.
Enfin un dernier avantage provient du problème qui existe en condition isotherme pour
amener l'échantillon à la température de réaction. En fait le début de l'expérience ne se
déroule pas à température constante, or la modélisation considère que durant toute la réaction,
la température ne varie pas. Dans le cas de la modélisation en condition non isotherme, ce
problème ne se présente pas.
Ces raisons ont rendu les expériences en programmation de température très utilisées.
Cependant il faut observer que l'étude directe en programmation de température ne

permet pas de prendre en compte un changement éventuel de régime cinétique avec la
température, ni d'observer la formation de nouveaux composés lorsque la température varie.
I.2.2 Classification des méthodes non isothermes
Nous distinguons deux types d'études réalisées à température variable :

• les études en programmation en température qui consistent en général à étudier la
réaction au cours d'une montée linéaire en température ;
• les études effectuées en analyse thermique à vitesse contrôlée.
- 42 -
I.2.2.1 Programmation en température
Galwey et Brown [1] ont fait une synthèse de ces méthodes. Celles-ci utilisent en
général des programmations linéaires en température en fonction du temps et appliquent le
système d'équations suivant :
∂α =kf(α) (I.31)
∂t
 E 
k = Aexp − a  (I.44)
 RT 
∂α = 1 ∂α (I.54)
∂T β ∂t
où β est la pente de la montée programmée en température.
Cependant comme nous l'avons vu précédemment, ces équations conduisent à deux

hypothèses :
• l'utilisation de la loi d'Arrhenius (I.44) de cette manière suppose la séparation des
variables température et pression, or ce n'est pas toujours le cas ;
• l'équation (I.31) est très restrictive car cette équation est seulement vraie pour des régimes
cinétiques correspondant à des modèles limites (germination instantanée ou croissance
instantanée).
De plus, au départ de l'utilisation des méthodes non isothermes, certains auteurs

écrivaient la fonction f(α) sous la forme :
f (α ) = (1 − α )
n
(I.55)
où n est appelé improprement l’ordre de la réaction. Le but est alors de déterminer la valeur
de n. Cette écriture provient de la cinétique homogène. De nombreux auteurs [48,49,50,51]
insistent sur l'impossibilité d'utiliser cette expression en cinétique hétérogène. Ils constatent
l'absence de signification physique de l'ordre de la réaction écrit sous la forme de l'équation
(I.55). En réalité n n'a de signification que si sa valeur est égale à 1, 2/3, 1/2 et 0 qui
correspond respectivement aux modèles F1, R3, R2, R1, et dans ce cas n n’est pas un ordre.
Ces méthodes d'étude en programmation de température peuvent être classées selon la

méthode du calcul développée pour obtenir la loi ou les paramètres cinétiques de la réaction.
- 43 -
On distinguera les méthodes dérivées des méthodes intégrales. Pour chaque cas, nous
donnerons quelques exemples d'études réalisées ainsi.
1) Méthodes dérivées
Elles consistent à reprendre les équations (I.31), (I.44) et (I.54) sous la forme logarithmique :
 dα 
ln dT  = ln A  − E a
  β  RT
(I.56)
 f(α )
 
Friedman [52] propose ainsi de représenter ln(dα/dt)i en fonction de 1/Ti, mesuré pour
un même degré d'avancement αi, pour différentes montées en température et utilise l'équation
(I.56) sous la forme :
(
ln d α
dt
)=ln Afβ(α ) − RT
E
a (I.57)
Les lignes parallèles, obtenues pour différents αi, ont une pente –Ea/R et une ordonnée
à l'origine ln(Af(αi)). Il est possible de déterminer A par extrapolation lorsque αi tend vers 0.
Cette méthode permet la détermination des paramètres cinétiques sans connaître f(α). Cela
suppose que f(α) existe, donc qu'il s'agisse d'un cas limite.
Il applique cette méthode à l'étude de la décomposition de polymère phénolique.
Kissinger [53] a étudié la décomposition de solides en programmation linéaire de

température en utilisant l’analyse thermique différentielle. Il détermine ainsi la température
correspondant à la vitesse maximale de la réaction (Tm). A partir des équations (I.31), (I.44) et
en supposant que la réaction obéisse à un ordre suivant l'équation (I.55), il obtient :
dα = Aexp − Ea  (1−α )n (I.58)

dt  RT 
Il détermine alors les équations donnant l’expression de α, dα/dT, d2α/dT2 et établit
une relation entre l’énergie d’activation, la vitesse de montée en température (β) et Tm :
 β 
d ln  2 
 Tm  = − E a (I.59)

d 1  R

 Tm 
En réalisant différentes expériences, à différentes montées en température, il en déduit
l’énergie d’activation. Il obtient la valeur de « l’ordre » de la réaction en utilisant les courbes
d’analyse thermique différentielle. Il applique ces résultats aux décompositions de plusieurs
- 44 -
solides (magnésite, calcite, brucite ...) et trouve des valeurs pour les différents paramètres
cinétiques (E,A,n) proches de celles trouvées en utilisant des méthodes isothermes.
Freeman et Caroll [54], à partir de l'équation (I.56), obtiennent la relation :
∆ln( dα )=n∆ln(1−α) − a ∆( 1 T )
E
(I.60)
dT R
Ainsi ils déterminent l'énergie d'activation en traçant [∆ ln(dα dT) / ∆ ln(1 − α)] en
fonction de [∆(1 T ) / ∆ ln(1 − α)] ou en traçant [∆ ln(dα dT ) / ∆ ln(1 T )] en fonction de
[∆ ln(1 − α) / ∆(1 T )] .
2) Méthodes intégrales
Ces méthodes sont basées sur l'intégration de l'équation (I.31) sous la forme :
α
g ( α )= ∫ d α = A ∫ exp  − E a dT
 (I.61)
0 f (α ) β  RT 
En effectuant un changement de variable, on obtient :
g ( α )= AE a p ( x )
βR
(I.62)
E ∞ e− x
avec x = a et p ( x )= ∫ dx
RT x x2
Afin de déterminer g(α), il faut évaluer la fonction p(x). Trois méthodes sont utilisées :
• utilisation de valeurs numériques de p(x) ;
• approximation de p(x) à l'aide d'une somme de termes d'une série ;
• approximation de p(x) à l'aide d'une fonction.
Ozawa [55] s’appuie au départ de son raisonnement sur une loi à "ordre" pour
exprimer f(α) et il fait une approximation de l'intégrale p(x) à l'aide d'une fonction
mathématique. Il établit une relation permettant de déterminer l’énergie d’activation. L’auteur
démontre qu’à un avancement donné, il existe une relation liant la pente (βi) de la montée en
température et la température (Ti) :
−log β 1 −0 .4657 Ea = −log β 2 −0 .4657 Ea = ... (I.63)

RT1 RT2
La valeur de l’énergie d’activation trouvée ainsi et la forme de la courbe expérimentale
permettent de déterminer l’ordre de la réaction et le facteur pré-exponentiel.
- 45 -
La méthode est appliquée à deux transformations : la décarbonatation de l'oxalate de

calcium en carbonate de calcium et la dégradation du nylon 6.
Coats et Redfern [56] utilisent eux aussi une équation avec un "ordre" pour décrire
une réaction du type aA(s) → bB(s)+cC(g). Les auteurs établissent, en faisant une approximation
sur p(x) et en le remplaçant par une somme de termes d'une série, les équations suivantes :
Ea
Y =− avec
2.3RT
 1−(1−α )n 
Y =log 2  pour n≠1 (I.64)
 T (1−n) 
 −log10(1−α ) 
Y =log  pour n=1

 T2 
En traçant les courbes Y=f(1/T), ils choisissent la valeur de n correspondant à une
droite et la pente de celle-ci leur permet d’accéder à l’énergie d’activation.
Ils utilisent cette méthode afin de déterminer les paramètres cinétiques de trois
réactions : la déshydratation de CaC2O4,H2O puis la décarbonatation de CaC2O4 et enfin la
décarbonatation de CaCO3.
Pan et al. [57] combinent les méthodes de Coats et Redfern [55] et d'Ozawa [56] pour
étudier en condition non isotherme la déshydratation de l'oxalate de manganèse dihydraté. Ils
aboutissent à une loi R2, f(α)=(1-α)1/2 et à une énergie d'activation égale à 162kJ/mol.
Carrasco [47] s’appuie sur une équation du type :

dα
− = k ( T ) f ( α )q( α ,T ) (I.65)
dt
D'après l'auteur, le terme q(α,T) doit prendre en compte les interactions entre la
température et le degré d'avancement. Ce terme dépend de plusieurs facteurs comme la
modification des propriétés physiques de l'échantillon durant la réaction. Il admet qu’étant
difficile à déterminer, celui-ci est égal à 1. Il suppose que f(α) s'écrit comme une loi à "ordre".
Il obtient, comme Kissinger [53], des équations donnant l’expression de α, dα/dT, d2α/dT2 et
établit la relation :
Ea =nRTm2
(ddTα )
m
(I.66)
(1−α m )
- 46 -
 dα 
où Tm, αm, et   sont la température, le degré d’avancement et la vitesse correspondant à
 dT  m
la vitesse maximum de la réaction.

Il applique cette méthode à la décomposition du polystyrène, compare les valeurs
trouvées pour l'énergie d'activation à d'autres méthodes et constate des écarts importants sur
ces valeurs selon la méthode utilisée.
Tomashevitch et al. [58] ont repris la méthode par intégration mais en opérant avec
une programmation en température selon :
1 1
= + γt (I.67)
T T0
Ils obtiennent une expression de f(α) sans approximation :
f ( α )= ∂ α 1 (I.68)
∂t  exp  − E aγ t 
A exp  − E a  
 RT 0   R 
Ils calculent alors f(α) grâce à l'expression précédente pour différentes valeurs de E
avec différentes pentes « réciproques » γ. Quelle que soit la pente, f(α) est le même pour une
valeur donnée de α. Ils tracent f(α) et choisissent la valeur de l’énergie d’activation qui
correspond le mieux à cette hypothèse, c'est-à-dire pour que, quel que soit Ea, toutes les
courbes f(α) soient confondues. Enfin ils déterminent le meilleur modèle qui correspond à la
courbe f(α) trouvée.
Ils appliquent cette technique à la réaction Fe2O3+Li2CO3=2LiFeO2+CO2.
Li et Tang [59] essaient aussi d'utiliser une méthode intégrale sans approximation.
Pour cela, ils passent par le logarithme de la vitesse avant d'intégrer :
( )
g ( α )= ∫0 ln d α d α = − a ∫0 d α + G ( α )
α
dt
E α
R T
(I.69)
α
avec G ( α )=α ln A + ∫0 ln ( f ( α ) ) d α
α  dα  α dα
Le tracé de ∫0 ln  d α en fonction de ∫0 , pour différentes pentes et à un
 dt  T
avancement donné, permet d'obtenir la valeur de l’énergie d’activation Ea.

Ils appliquent leur méthode à la réaction SrCO3=SrO+CO2.
- 47 -
Afin de vérifier le bon accord entre toutes ces méthodes, des études ont consisté à
vérifier les paramètres cinétiques obtenus en utilisant différentes méthodes. Ainsi, Koga et
Tanaka [60] ont étudié la déshydratation du sulfate de lithium monohydraté en utilisant à la
fois des méthodes isothermes et non isothermes. Des paramètres d'Arrhenius différents sont
obtenus selon les méthodes utilisées. Par certaines méthodes, ils obtiennent une énergie
d'activation qui évolue avec α. Cette variation de l'énergie d'activation en fonction de
l'avancement de la réaction confirme que les méthodes et les modèles de transformation
utilisés ne sont pas toujours adaptés aux réactions étudiées.
Achar et al. [61] recommandent l'utilisation simultanée d'expériences isothermes et

non isothermes. Ils se servent d'expériences isotherme-isobare dans la détermination de f(α)
puis des montées en température pour obtenir les paramètres d'Arrhenius. Ils appliquent cette
méthode à la déshydroxylation de la kaolinite.
Ce paragraphe donne quelques exemples du traitement des données cinétiques en

programmation de température. Les défauts de ces méthodes proviennent essentiellement de
leur point de départ : l'application uniquement à des modèles limites et l'utilisation de la
loi d'Arrhenius pour décrire l'évolution de la constante cinétique, k, en fonction de la
température, ce qui suppose notamment la séparation des variables, température et pression.
I.2.2.2 Analyse thermique à vitesse contrôlée (ATVC)
Cette méthode, appelé également ATCE (Analyse thermique contrôlée par

l’échantillon) et développée par Rouquérol [38], consiste à asservir le chauffage de
l'échantillon à une grandeur expérimentale (pression d'un gaz, débit gazeux, variation de
masse, puissance thermique ...) afin que celle-ci reste constante tout au long de la réaction.
Cette technique permet de séparer des réactions successives, d'obtenir des solides de
compositions précises et surtout d'éviter les problèmes de transferts de chaleur, de masse et de
gaz au sein du réacteur d'étude à condition que la vitesse de réaction soit suffisamment lente.
Criado et al. [62] ont calculé l'allure des courbes d'analyse thermique à vitesse
contrôlée, α(T), pour les lois répertoriées par Sharp et al. (tableau I.4). La forme de ces
courbes permet de différencier les différents modèles de transformation (An, Dn, Rn, F1).
- 48 -
Cependant, comme le montre la figure I.10, qui représente les courbes obtenues en ATVC
pour les modèles R2 et R3, il n'y a par contre pas beaucoup de différence entre les lois d'une
même famille.
Figure I.10 : Courbe α(T) obtenue en ATVC pour les modèles R2 et R3 [62]
Ces études, menées en analyse thermique à vitesse contrôlée, se limitent pour l'instant
aux modèles limites et supposent la séparation des variables au niveau de la constante
cinétique.
La détermination de l'énergie d'activation peut être réalisée à l'aide de la méthode des

créneaux de vitesse développée par Rouquérol et Rouquérol [63] qui est similaire à la
méthode des décrochements. L'obtention de cette énergie d'activation a lieu sans avoir besoin
de connaître le modèle (à condition que ce soit un cas limite). Cette méthode consiste à faire
varier la vitesse de la réaction entre deux valeurs C1 et C2 au cours de l'expérience et à
mesurer à chaque fois les températures avant (T1) et après (T2) le saut de vitesse. La valeur de
l'énergie d'activation est alors obtenue par la formule :
E a = R T1T2 ln C 2  (I.70)

T2 −T1  C1 
Remarque : Flynn et al. [64] utilisent également cette méthode de décrochements en analyse
thermique classique en arguant l'avantage par rapport aux méthodes isothermes où le passé réactionnel
peut influencer la vitesse de réaction. Avec cette méthode, tous les échantillons ont la même histoire
thermique. En utilisant les équations (I.31) et (I.44), ils obtiennent l'expression de l'énergie d'activation
en réalisant des sauts en température :
E a = R T1T2 ln v2  (I.71)

T2 −T1  v1 
- 49 -
T2 et T1 étant les températures avant et après le décrochement, v2 et v1 étant les vitesses obtenues avant
et après le décrochement. La valeur de la vitesse après le décrochement est extrapolée à l'instant du
décrochement.
L'analyse thermique à vitesse contrôlée est de plus en plus utilisée pour étudier la
cinétique d'une réaction. Par exemple, Tiernan et al. [51] ont étudié par ATVC la réduction
de CuO et NiO par l'hydrogène en gardant constante la vitesse de production de vapeur d'eau.
Par la méthode des créneaux, ils obtiennent une variation de l'énergie d'activation avec le
degré d'avancement de la réaction. Ils en déduisent par un raisonnement qualitatif la présence
et l'effet simultané de la germination et de la croissance sur les courbes cinétiques. Ils mettent
ainsi en évidence la nécessité de prendre en compte les deux phénomènes. Ils suggèrent aussi
qu'il faut mesurer une énergie d'activation pour la germination et une pour la croissance.
Par conséquent, il serait très intéressant de prendre en compte les modèles complexes
de germination-croissance en analyse thermique à vitesse contrôlée.
I.3 Conclusions
Dans cette revue bibliographique, nous avons d’abord tenté de comprendre l’origine des
lois de vitesse utilisées pour l’étude cinétique des réactions de transformation des solides. En
recherchant les hypothèses utilisées lors de l’établissement de ces lois, il a été possible de
recenser les différentes sortes de modèles en se basant sur les équations (I.14) et (I.18) et sur
les hypothèses suivantes :
• la pseudo-stationnarité de la croissance,
• la germination aléatoire en surface,
• l’étape limitante de la croissance,
• le mode de croissance des germes (isotrope ou anisotrope),
• le sens de développement de la nouvelle phase,
• la forme géométrique du solide initial.
- 50 -
Remarque : un logiciel (« CIN ») regroupe l’ensemble des modèles du tableau I.5, permettant
ainsi d’avoir rapidement accès à l’allure théorique des courbes (α(t) et dα/dt(α)) et à la détermination
de φ et/ou γ.
Les critiques que nous avons pu émettre lors de cette analyse s’ajoutent à celles déjà
exprimées par de nombreux auteurs. En effet, Brown et Galwey [65] faisaient remarquer en
2000 le manque d’une théorie générale et la nécessité de la prise en considération de la
complexité des phénomènes. Baram [66] reproche notamment l'utilisation de modèle en
cinétique se limitant en général à des cas où la germination est instantanée.
De plus comme nous l'avons montré, de nombreuses études supposent la séparation

des variables température et pression sans la justifier. Ainsi, certaines études montrent que
l'énergie d'activation varie selon la méthode cinétique non isotherme utilisée et selon le degré
d'avancement de la réaction.
Certains auteurs remettent en cause les études menées en programmation de

température. Ainsi Boldyreva [67] critique celles-ci pour leur manque de signification
physique. Il explique que le seul avantage de ces méthodes est leur rapidité pour
éventuellement caractériser un procédé mais elles ne donnent pas d'informations sur l'aspect
physico-chimique de la réaction. Ainsi, seules les études isothermes peuvent permettre la
modélisation cinétique d'une réaction.
De plus, Galwey et Brown [1] expliquent que l'utilisation systématique de ces
méthodes ne permet pas de détecter des cas cinétiques différents. Ils critiquent l'utilisation de
ces méthodes comme des « boites noires » sans compréhension des phénomènes. D'ailleurs,
ces méthodes d'études directes en programmation de température ne permettent pas d'observer
un changement de régime cinétique ou l'apparition de nouveau produit formé lorsque la
température varie. Les modèles développés supposent également qu’une seule réaction se
déroule. Cependant il n’est parfois pas possible de séparer « thermiquement » deux réactions
successives. Comment modéliser alors ce type de réaction ?
Galwey et Brown [65] dénoncent également le trop grand nombre d'études dont
l'intérêt est la résolution mathématique des lois cinétiques appliquées en condition non
isotherme plutôt que l'étude physico-chimique des phénomènes mis en jeu.
Enfin il ne faut pas sous-estimer l'importance de la caractérisation des échantillons. Au
niveau expérimental, plusieurs auteurs [65,68] insistent sur la nécessité de réaliser des
- 51 -
expériences de caractérisation et en particulier des observations microscopiques des solides en

début et en cours de transformation. Ces observations peuvent permettre de choisir un modèle
de transformation selon la forme des solides, la formation des germes, le développement de la
nouvelle phase, les gaz produits par la réaction, les espèces adsorbées en surface ...
Dans le chapitre suivant, nous proposons une nouvelle méthodologie, afin de rendre
compte des réactions se déroulant en conditions non isotherme mais aussi non isobare, qui
permet de résoudre une partie des problèmes soulignés précédemment. Ainsi, nous essaierons
de développer une modélisation en conditions non isotherme et non isobare valable pour
l'ensemble des modèles (modèles limites et modèles complexes de germination croissance) et
également applicable aux systèmes mettant en jeu des réactions successives.
- 52 -
Chapitre II Approche non isotherme et non isobare de la modélisation cinétique
C
Chhaappiittrree IIII
APPROCHE NON ISOTHERME ET NON

ISOBARE DE LA MODELISATION CINETIQUE
Le but de cette étude est de développer des modèles permettant de décrire les réactions
solide-gaz en conditions non isotherme et non isobare. Afin de réaliser ce travail, nous allons
d'abord détailler une méthodologie pour mettre en place ces modèles. Celle-ci sera ensuite
appliquée pour les deux réactions étudiées : la déshydroxylation de la kaolinite et la réduction
de l'octooxyde de triuranium par l'hydrogène.
Pour rendre compte de l'évolution d'une réaction, nous allons développer d'une part un
modèle de transformation permettant de décrire l'évolution des zones réactionnelles au
cours du temps et de déterminer la zone réactionnelle où se déroule l'étape limitant la
croissance et d'autre part un modèle physico-chimique qui expliquera à l'aide de mécanismes
réactionnels comment évoluent les vitesses surfaciques de croissance et de germination en
fonction des paramètres intensifs de la réaction.
L'obtention de ces deux modèles nous permettra ensuite de connaître comment la
vitesse de réaction varie lorsque les paramètres physico-chimiques sont modifiés au cours du
temps.
Afin de présenter cette approche de l'étude cinétique de réaction en conditions non

isotherme et non isobare, nous montrerons d'abord l'intérêt de cette étude. Nous détaillerons
ensuite la méthodologie qui sera utilisée afin de modéliser ces réactions. Nous terminerons en
présentant les systèmes réactionnels choisis pour ce travail.
- 53 -
II.1 Intérêts de la modélisation
Les objectifs de ce travail sont de pouvoir décrire, à l'aide de modèles, des réactions
solide-gaz lorsque la température ou la pression des gaz varie au cours du temps. Cette étude
entre dans le cadre d'un programme de recherche dont le but est la modélisation de fours
industriels où les conditions de température et de pression fluctuent. Elle sera donc utilisée par
la suite dans ce programme de recherche afin d’introduire les termes sources issus des
réactions.
Par ailleurs l'étude bibliographique du premier chapitre a montré que les modèles
développés en condition non isotherme utilisent des lois « simples » pour exprimer la
variation des constantes cinétiques en fonction de la température et de la pression :
• la température influe toujours sur la constante cinétique selon une loi d'Arrhenius ;
• pour les décompositions solides, ces études ne prennent pas généralement en compte
l'influence de la pression ; pour les réactions solide-gaz, la vitesse dépend de la pression
en gaz réactant selon une loi à « ordre ».
Comme nous l'avons montré, il est possible de déterminer l'évolution des grandeurs
cinétiques (réactivité surfacique de croissance et fréquence surfacique de germination) en
fonction de la température et de la pression. Nous prendrons en compte ces lois, dont les
expressions seront différentes pour chaque réaction, dans l'utilisation de la modélisation non
isotherme et non isobare.
De plus l'étude bibliographique a également permis de conclure que les méthodes

développées actuellement en condition non isotherme ne sont pas satisfaisantes pour atteindre
précisément nos objectifs. En effet, d’une part, elles font appel aux lois dans lesquelles la
vitesse est une fonction du degré d’avancement, et donc ne permettent pas de mettre en
équation le cas des réactions de germination-croissance (isotrope ou anisotrope), pour
lesquelles la vitesse est également fonction du « passé » du solide, c’est-à-dire des variations
de température et de pression. D’autre part, ces méthodes supposent des transformations à
réaction unique et ne peuvent pas être utilisées simplement dans le cas de transformations
successives.
- 54 -
Un autre aspect important de cette étude et qui sera développé dans les prochains
chapitres est la difficulté de réaliser des expériences réellement en conditions isotherme-
isobare. En effet, afin de démarrer une réaction, un des paramètres est modifié : montée rapide
en température, envoi du gaz réactant dans le réacteur.
Au départ, il existe donc toujours une durée pendant laquelle une des contraintes
intensives varie. Or, souvent le modèle, qui interprète ces réactions, considère un système
isotherme-isobare et ne prend pas en compte cette variation. Comme nous le verrons, il sera
possible à partir du modèle développé dans un premier temps en conditions isotherme et
isobare de tenir compte cette fois, du changement initial de température ou de pression.
Ainsi, notre travail consiste à développer des modèles permettant de décrire les
systèmes réactionnels en conditions non isotherme et non isobare aussi bien pour les modèles
limites que pour les modèles complexes de germination-croissance ou encore pour des
réactions successives.
1 Cas des transformations avec germination-croissance
Nous allons illustrer, à l’aide d’un exemple simple, l’influence du passé du solide sur
la vitesse de réaction à un degré d’avancement donné.
Prenons le modèle de germination-croissance anisotrope avec un développement
interne et une réaction limitant la croissance située à l'interface interne. Ce modèle a
notamment permis à Bouineau [14] de décrire la déshydratation du sulfate de lithium
monohydraté. Les grains sont supposés sphériques avec un rayon égal à 0,2 µm. Nous
supposons une réaction pouvant être décrite à l'aide de ce modèle. Deux expériences sont
réalisées. La première se déroule en conditions isotherme et isobare avec une température T1
et une pression P1 correspondant à des valeurs de la fréquence surfacique de germination de
γ1=109 germes.m-2.s-1 et de la réactivité surfacique de croissance de φ1=5.10-7 mol.m-2.s-1. La
deuxième expérience débute dans des conditions différentes avec une température T2 et une
pression P2 correspondant à γ2= γ1/1,5 et φ2= φ1/3. Pour cette seconde expérience, on réalise
un rapide décrochement à un temps égal à 2000 secondes afin de revenir aux conditions
expérimentales de la première expérience (T1,P1). Ainsi après le décrochement, les deux
expériences se déroulent dans les mêmes conditions. Or, comme le montre la figure II.1 qui
représente pour ces deux expériences la vitesse en fonction du degré d'avancement, la vitesse
- 55 -
de la réaction n'est pas la même après le décrochement. Ainsi, le passé réactionnel influe sur
la vitesse de réaction.
0.0003 T ,P
1 1
0.00025
T ,P
d α /dt (s )
-1
1 1
0.0002
0.00015
0.0001
T ,P
2 2
-5
5 10 Expérie nce iso therme-iso bare
Expérie nce avec un décrochement
0
0 0.2 0.4 0.6 0.8 1
α
Figure II.1 : Influence du passé de la réaction, modèle complexe de germination-croissance
Si on prend le même modèle pour la croissance des grains avec les même paramètres
pour la réactivité surfacique de croissance, mais cette fois-ci avec une germination
instantanée, nous pouvons également calculer la vitesse de réaction pour une expérience
isotherme et une expérience avec un décrochement. Pour ce modèle limite, comme le montre
la figure II.2, le passé de la réaction n'a pas d'influence sur la vitesse.
0.0007 Expérience isotherme-isobare

T ,P
1 1 Expérience avec un décrochement
0.0006
0.0005
d α /dt (s )
-1
T ,P
1 1
0.0004
0.0003
0.0002
T ,P
2 2
0.0001
0
0 0.2 0.4 0.6 0.8 1
α
Figure II.2 : Influence du passé de la réaction, modèle limite de germination instantanée
- 56 -
Cet exemple montre bien que le passé « thermique » de la réaction a une influence sur
la vitesse dans les cas où il existe une compétition entre la germination et la croissance. Il
n'est pas alors possible d'écrire la vitesse de la réaction sous la forme dα/dt =kf(α).
Remarque : les deux figures II.2 et II.3 correspondent en fait aux expériences réalisées pour le
test du cas limite (cf. I.3).
Ainsi dans le chapitre III, nous essaierons de développer un modèle de germination-

croissance en conditions non isotherme et non isobare pour ce type de réaction où le passé
influe sur la vitesse
1 Cas des transformations avec réactions successives
Comme nous l'avons vu dans le premier chapitre, les méthodes utilisées en

programmation linéaire de température sont parfois trop systématiques. Elles ne permettent
pas d'interpréter la vitesse d'une réaction qui se déroule en plusieurs transformations
successives qui ne peuvent pas être séparées thermiquement.
Ainsi, il est intéressant de voir comment modéliser ce type de réactions, qui comprend
des phases solides intermédiaires, en conditions non isotherme et non isobare.
Dans ce cas, plusieurs lois de vitesse sont à prendre en considération correspondant à
chaque transformation. Ces différentes réactions successives pourront alors être modélisées,
soit comme des cas limites, soit comme des cas de germination-croissance.
Remarque : l'étude par analyse thermique à vitesse contrôlée peut permettre parfois de
séparer ces réactions successives. Il est alors possible dans ce cas de modéliser chaque réaction
indépendamment.
Ainsi, nous montrerons au chapitre IV comment modéliser ce type de réaction en

conditions non isotherme et non isobare.
Nous venons de montrer les raisons qui motivent cette étude. La première étape est de
mettre en place une méthodologie pour développer la modélisation de réactions en conditions
non isotherme et non isobare.
- 57 -
II.2 Principe de la modélisation en conditions non isotherme et

non isobare
Comme nous l'avons évoqué dans le premier chapitre, la modélisation cinétique d'une
transformation mettant en jeu des solides doit expliquer l'évolution de la vitesse d'une réaction
en fonction du temps et des paramètres physico-chimiques. Le modèle prendra alors en
compte à la fois les variations géométriques du système réactionnel et les étapes élémentaires
mises en jeu au cours de la transformation.
Voyons d'abord la méthodologie utilisée pour les réactions correspondant à des
modèles de germination-croissance.
II.2.1 Cas des réactions avec germination-croissance
La modélisation en condition non isotherme ou non isobare découle d'une part, du

modèle de transformation qui permet d'obtenir l'évolution de la fonction d'espace en fonction
du temps, et d'autre part, de l'obtention des paramètres φ et γ en fonction de la température et
de la pression.
Afin d'obtenir un modèle cinétique décrivant ces réactions réalisées à température et
pressions variables, nous allons procéder en plusieurs étapes :
1 Etude expérimentale
Ce travail a deux objectifs : caractériser les échantillons étudiés d'un point de vue
géométrique et structural et déterminer les courbes cinétiques expérimentales en conditions
isotherme-isobare. Celles-ci seront obtenues par thermogravimétrie.
2 Vérifications expérimentales des hypothèses de la modélisation
Nous supposons que la vitesse peut se mettre sous la forme de l'équation I.18,
dα/dt=φE, c'est-à-dire que le mécanisme de croissance peut s’écrire comme une succession
d’étapes élémentaires et qu’une de ces étapes peut être considérée comme limitante. Ceci
suppose que le système réactionnel est pseudo-stationnaire. Cette hypothèse sera donc
- 58 -
vérifiée à l'aide de la méthode présentée dans l'annexe B en réalisant le couplage

thermogravimétrie-calorimétrie. Le montage expérimental utilisé pour ce couplage est
également présenté dans cette annexe.
De plus, comme nous l'avons vu au premier chapitre le test du φE, décrit dans
l'annexe C, permet de conclure si la vitesse peut se décomposer sous la forme de l’équation
I.18. Ce test sera également vérifié.
Enfin un choix pour le modèle de transformation va être de savoir si la réaction suit un
modèle limite ou un modèle complexe de germination-croissance. Le test du cas limite nous
donnera un élément de réponse.
3 Détermination du modèle de transformation
Les hypothèses précédentes étant vérifiées, nous allons chercher un modèle qui, à
partir d'hypothèses sur la germination et la croissance, permet de décrire comment la nouvelle
phase va se développer sur les grains de la phase initiale. Afin de choisir quel modèle pourra
décrire une réaction, nous nous appuierons à la fois sur la forme des courbes cinétiques
expérimentales et sur les expériences de caractérisations des produits réalisées avant et
pendant la réaction.
Nous calculerons la vitesse et la fonction d'espace relatives à ce modèle dans le cas de
grains de même taille mais également pour une poudre présentant une distribution
granulométrique. Pour cela nous calculerons le degré d'avancement pour chaque taille de
grain α(r,t) à l'aide du modèle et nous intégrerons cette expression sur l'ensemble des tailles
de grain.
4 Comparaison avec l'expérience
Les courbes de vitesses obtenues dépendent de deux paramètres, la fréquence

surfacique de germination et la réactivité surfacique de croissance. Comme nous l'avons
présenté dans le premier chapitre, nous ajusterons ces paramètres pour que les courbes
théoriques correspondent aux courbes cinétiques expérimentales (cf. annexe D). Nous
vérifierons le bon accord entre le modèle et l'expérience en comparant les courbes dérivées
α/dt en fonction de α qui permettent de mieux apprécier les écarts
représentant la vitesse dα
entre le modèle et l'expérience plutôt que les courbes α(t).
- 59 -
5 Détermination des variations de φ et γ en fonction des conditions physico-

chimiques
La comparaison des courbes de vitesse théoriques et expérimentales permet d'obtenir

les valeurs de la fréquence surfacique de germination et de la réactivité surfacique de
croissance dans les différentes conditions où les expériences ont été réalisées (cf. annexe D).
Il est également possible de déterminer l'évolution de la réactivité surfacique de croissance en
fonction d'une grandeur physico-chimique, directement à partir d'expériences sans connaître
le modèle géométrique, à l'aide de la méthode des décrochements, décrite dans le premier
chapitre.
Remarque : si la variation de la réactivité surfacique de croissance en fonction d’une des

contraintes intensives est connue à la fois à l’aide de la comparaison avec le modèle et par la méthode
des décrochements, nous pouvons comparer le résultat obtenu à partir de ces deux méthodes. Si celui-
ci est le même, alors le modèle de transformation choisi semble bon puisqu’il est en accord avec les
valeurs obtenues par la méthode des décrochements qui est uniquement expérimentale.
6 Détermination du modèle physico-chimique
Cette partie de la modélisation consiste à décrire les mécanismes réactionnels de

germination et de croissance à l'aide de plusieurs étapes élémentaires. Afin d'écrire ces étapes,
nous utiliserons les notations de Kröger [69] qui rendent compte de la structure et des défauts
des solides. En considérant une étape limitante, il est possible de calculer les lois d'évolution
de ces grandeurs par rapport aux conditions physico-chimiques. La détermination de l'étape
limitante de ces mécanismes se fera en comparant les lois obtenues φ(T,P) et γ(T,P) avec les
valeurs de φ et γ déterminées précédemment.
7 Réalisation d'expériences en conditions non isotherme ou non isobare
Ces expériences seront également réalisées par thermogravimétrie en programmation

de température ou de pression. Comme nous le verrons, il sera parfois nécessaire de reprendre
les expériences réalisées en conditions « isotherme et isobare » puisqu'au départ de ces
expériences, il y a toujours un des paramètres physico-chimiques qui varie.
- 60 -
8 Calcul de la vitesse théorique obtenue en conditions de température ou de

pression variables et comparaison aux courbes expérimentales
Les étapes 1 à 2 concernent le développement d'un modèle en conditions isotherme

et isobare. Elles sont nécessaires pour l'obtention de la vitesse dans des conditions différentes.
L'étape 3 permet de connaître la fonction d'espace (E(t)). Les variations des grandeurs
cinétiques en fonction des contraintes intensives (φ(T,P et γ(T,P)) de la réaction sont obtenues
par les étapes 4 et 2.
La connaissance de ces variations va permettre de calculer l'évolution de la fonction
d'espace lorsque les paramètres physico-chimiques varient. Il sera bien sûr nécessaire de
savoir comment la température ou la pression partielle d'un gaz évolue au cours de la réaction.
L'expression de la vitesse calculée peut alors s'écrire sous la forme :
dα =φ(T,P)E(t,φ(t),γ(t)) (II.1)
dt
Le degré d'avancement de la réaction est obtenu par intégration de la vitesse. Ainsi, si

l'étude en conditions isotherme-isobare a été bien effectuée, tous les paramètres de l'équation
de la vitesse sont connus. Il est alors aisé d'obtenir l'expression de la vitesse de la réaction en
conditions non isotherme et non isobare. Il reste ensuite à comparer cette vitesse calculée à la
vitesse expérimentale obtenue pour des variations de température ou de pression des gaz au
cours des réactions étudiées (étape 2).
Remarque : si les expressions des fonctions φ(T,P) et γ(T,P) sont connues mais pas les
paramètres relatifs à celles-ci (énergies d'activation apparentes …), il est possible de déterminer ces
paramètres en comparant directement le modèle aux expériences à température ou pression variable.
La figure II.3 présente un schéma qui résume la méthode pour déterminer un modèle
en conditions non isotherme et non isobare.
- 61 -
2 Vérification de la
pseudo-stationnarité du 7 Expériences non
1 Caractérisations "φE" et test du cas limite isotherme et/ou non-
isobare
3 Modèle de transformation
1 Expériences - grains de même taille : E(t)
"isotherme-isobare" - distribution granulométrique : Edg(t)
5
8 Modélisation en
4 condition non
γ(T,P) et φ(T,P) isotherme et non
isobare
6 4 6
Mécanisme de Méthode des Mécanisme de

germination décrochements croissance
Figure II.3 : Schéma de la méthodologie pour obtenir un modèle décrivant des réactions en
conditions non isotherme et non isobare
II.2.2 Cas des réactions successives
Nous pouvons également proposer une méthodologie pour déterminer les modèles
cinétiques relatifs à des réactions successives en conditions non isotherme et non isobare.
Celle-ci est proche de la méthodologie développée pour un cas de germination-
croissance.
La première partie concerne toujours l'obtention des courbes cinétiques expérimentales
et la caractérisation des échantillons, mais une étude plus complète sera réalisée sur
l'évolution de la structure du solide au cours de la réaction. Nous essaierons, ainsi, de mettre
en évidence les phases solides intermédiaires qui se forment.
Ensuite nous appliquerons la démarche précédente à savoir la détermination d'un
modèle de transformation et d'un modèle physico-chimique pour chaque réaction.
Une étape supplémentaire consistera à relier le degré d'avancement global mesuré
expérimentalement au degré d'avancement réel de chaque transformation.
Enfin nous terminerons en comparant les expériences réalisées en conditions non
isotherme et non isobare au modèle cinétique dans ces conditions.
- 62 -
En conclusion, l'étude cinétique de réactions en conditions non isotherme et non

isobare passe par plusieurs étapes. Il est nécessaire de vérifier si cette méthodologie est
applicable à l'aide de cas concrets.
II.3 Choix des systèmes réactionnels
Afin d'appliquer l'approche cinétique développée dans ce chapitre, deux

transformations différentes ont été étudiées : l'une consiste en un cas de germination-
croissance (cf. II.2.1) et l'autre en plusieurs réactions successives (cf. II.2.2).
Le chapitre III fera l'objet de la décomposition thermique d'une argile : la

déshydroxylation de la kaolinite qui est un cas de germination-croissance. L'équation bilan de
cette réaction s'écrit :
Si2Al2O5(OH)4 = Si2Al2O7 +2H2O
Nous avons choisi de travailler sur cette transformation car actuellement aucune étude
de modélisation cinétique n'a permis d'expliquer cette réaction sur tout le domaine
réactionnel.
Le chapitre IV traitera de la réduction de l'octooxyde de triuranium par l'hydrogène :

U3O8 + 2H2 = 3UO2 + H2O
Cette réaction comprend en fait plusieurs réactions successives. Elle intervient au sein
d'un réacteur industriel dans la fabrication du combustible nucléaire. Un des partenaires du
programme de recherche étant la société Comurhex qui réalise cette transformation, l'étude de
cette réduction portait donc un intérêt supplémentaire. De plus, celle-ci avait déjà été étudiée
en conditions isotherme et isobare par Brun [17].
Pour ces deux réactions, l'étude développée comportera trois étapes :

• une étude bibliographique ;
• une étude expérimentale ;
• le développement d'un modèle cinétique applicable en conditions non isotherme et non
isobare.
- 63 -
- 64 -
Chapitre III La déshydroxylation de la kaolinite
C
Chhaappiittrree IIIIII
MODELISATION D'UNE REACTION AVEC

GERMINATION-CROISSANCE :
LA DESHYDROXYLATION DE LA KAOLINITE
L'étude de la déshydroxylation de la kaolinite a été réalisée en collaboration avec Kais

Nahdi de l'université de Bizerte (Tunisie) dans le cadre de sa thèse.
Cette réaction est une décomposition thermique qui intervient dans la fabrication de
nombreux produits céramiques : porcelaine, briques, filtres céramiques ...
L'équation bilan de cette transformation s'écrit sous la forme :
La nouvelle phase formée est appelée métakaolinite.
Après une présentation bibliographique de cette réaction, nous exposerons dans un

premier temps les expériences effectuées dans le cadre de ce travail, puis dans un deuxième
temps, nous développerons un modèle afin d'interpréter les courbes cinétiques expérimentales.
Ce modèle sera d'abord présenté en conditions isotherme-isobare, puis nous verrons son
utilisation lorsque la température varie au cours de la réaction.
Mais pour débuter ce travail, présentons ce qu'est la kaolinite.
- 65 -
III.1 Etude bibliographique
III.1.1 Structure de la kaolinite
La kaolinite est un minéral argileux de formule chimique Si4Al4O10(OH)8 constitué

par l'empilement de feuillets. Comme le montre la figure III.1, un feuillet de kaolinite
comprend une couche de tétraèdes à cœur de silicium, SiO4, reliée à une couche d'octaèdre à
cœur d'aluminium. Deux tiers des sites octaédriques sont occupés par des atomes
d'aluminium.
Les groupements hydroxydes occupent certains sommets des octaèdres :
• trois hydroxydes sur quatre occupent des positions au bord des feuillets et sont appelés
hydroxydes externes ;
• un hydroxyde sur quatre est situé entre le plan des atomes de silicium et le plan des
atomes d'aluminium et est appelé hydroxyde interne.
Figure III.1 : Structure de la kaolinite
La kaolinite présente un système cristallographique triclinique. Les paramètres

cristallographiques de la maille de la kaolinite sont [70] :
a=5,16 Å b=8,95 Å c=7,41 Å
α=91,7° β=104,9° γ=89,8°
- 66 -
Les cristaux de kaolinite se présentent généralement sous la forme de plaquettes

hexagonales allongées (parfois réduites à des losanges). Ces plaquettes sont délimitées par des
faces basales et des faces latérales. Les dimensions de ces grains sont comprises entre 0,05 et
2 µm en ce qui concerne le diamètre des faces basales et entre 300 et 2000 Å pour l'épaisseur
(faces latérales).
La kaolinite est généralement associée à l'état naturel à d'autres minéraux tels que le
mica ou le quartz.
L'étude par analyse thermique différentielle de la kaolinite montre plusieurs

transformations [70] :
• une déshydratation avec départ de l'eau physisorbée autour de 100°C (pic

endothermique) ;
• une déshydroxylation vers 500°C (pic endothermique) qui correspond à la

transformation de la kaolinite en métakaolinite selon l'équation bilan :
Le produit ainsi formé appelé métakaolinite n'a pas de structure cristalline ordonnée ;
• une recristallisation avec deux pics exothermiques vers 980°C et 1200°C, correspondant
dans un premier temps à la formation d'une phase transitoire, puis à la formation d'une
nouvelle phase, la mullite.
Nous nous intéressons à la deuxième réaction : la déshydroxylation de la kaolinite.
III.1.2 Etude thermodynamique de la réaction
L'étude thermodynamique de la réaction a été réalisée afin de connaître les conditions

de températures et de pressions favorables à la déshydroxylation de la kaolinite.
Le diagramme d'équilibre thermodynamique, obtenu à l'aide des données de la
littérature (« Thermochemical properties of inorganic substances ») [71], est présenté sur la
- 67 -
figure III.2 avec le tracé de la pression de vapeur d'eau à l'équilibre en fonction de la

température. Les domaines de stabilité de chaque phase sont indiqués sur ce diagramme.
(hPa) 1000 Kaolinite
10
(H20)eq
0.1
0.001
P
-5
Métakaolinite
10
-7
10
200 300 400 500 600 700
T(°C)
Figure III.2 : Diagramme d'équilibre thermodynamique de la déshydroxylation de la kaolinite
1 Choix des conditions expérimentales
Au vu de ce diagramme, nous avons décidé de travailler dans un domaine de

températures de réaction comprises entre 450°C et 500°C et de pressions de vapeur d'eau
comprises entre 10-3 hPa et 10 hPa (domaine hachuré sur la figure III.2).
1 Des données thermodynamiques fiables ?
Les calculs thermodynamiques de la pression de vapeur d'eau à l'équilibre ont

également été réalisés à l'aide du logiciel Thermodata développé à l'Institut National
Polytechnique de Grenoble (programme Coach). Cependant, les données de ce logiciel sur
l'enthalpie standard de formation de la métakaolinite diffèrent d'environ 5% par rapport aux
données précédentes. Comme le montre la figure III.3, cette différence modifie
considérablement les courbes de la pression d'équilibre en fonction de la température.
A partir de cette figure, il apparaît qu'avec les données du logiciel Thermodata, dans
les conditions de températures et de pressions de vapeur d'eau choisies pour les expériences,
une grande partie de celles-ci ne devraient pas se dérouler.
- 68 -
4
10
100 Kaolinite
(hPa) 1
(H2O ) eq
0.01
0.0001
Métakaolinite
P
-6
10
D onnées "therm ochemic al properties
-8 of inorganic substances"
10
D onnées "therm odata" [71]
-10
10
300 350 400 450 500 550 600 650
T(°C)
Figure III.3 : Diagramme thermodynamique ; comparaison entre des données de différentes origines
Ce résultat met en avant une certaine variabilité des données thermodynamiques

accessibles. Nous verrons par la suite quelles données nous semblent correctes.
III.1.3 Cinétique de la déshydroxylation de la kaolinite
Les mécanismes proposés pour décrire cette réaction font en général intervenir
plusieurs étapes dont une correspond à la dissociation d'ions hydroxydes, et une autre à la
formation d'eau. Selon certains auteurs, les groupements hydroxydes sont supposés adjacents
et la déshydroxylation est dite « homogène » [72-74], alors que d'autres pensent que les
hydroxydes qui réagissent ne sont pas proches et la déshydroxylation est dite « hétérogène »
[75-78]. Dans le deuxième cas, il est nécessaire de supposer une diffusion et une production
d'eau à la surface des grains.
Les études cinétiques réalisées sur cette réaction, soit en condition isotherme, soit en
programmation linéaire de température, utilisent les expressions mathématiques de la vitesse
qui correspondent seulement à des cas limites de germination ou croissance instantanée (en
particulier les lois de Sharp et al [19]). Les résultats de ces travaux diffèrent selon la nature
de l'atmosphère gazeuse. Les études menées sous pression atmosphérique [72,77,78,79,80]
proposent un modèle de type F1, c'est-à-dire un modèle de germination aléatoire et croissance
instantanée. Les travaux réalisés sous vide ou sous faible pression de vapeur d'eau sont plus
- 69 -
en accord avec un modèle de type D3 [81-84] ou D2 [61]. Comme nous l'avons évoqué dans le
premier chapitre, la loi D3, bien que souvent utilisée dans la littérature, n'a pas de signification
physique. La loi D2 correspond à un modèle de germination instantanée avec un mécanisme
de croissance limité par la diffusion selon deux directions (géométrie cylindrique).
Dans tous les cas, quelle que soit la loi utilisée pour décrire la réaction, il ne semble
pas possible de rendre complètement compte des expériences réalisées. Le modèle ne permet
de décrire l'expérience que jusqu'à un degré d'avancement environ égal à 0,7. Certains auteurs
[85,86] proposent un changement de mécanisme réactionnel afin de justifier cette non
correspondance du modèle en fin de réaction.
Cette réaction a également été étudiée par analyse thermique à vitesse controlée
[70,87,88]. Un modèle a ainsi été proposé en imaginant une compétition entre les lois F1 et D3
[70]. La forme des courbes obtenues pour des kaolinites différentes a permis de montrer que
la présence de défauts avait une influence sur le modèle à prendre en compte pour la
déshydroxylation. Il semblerait que pour une kaolinite bien cristallisée, la loi F1 soit
prépondérante, alors que dans le cas contraire la loi cinétique correspondrait plus à la loi D3.
L'influence de la pression de vapeur d'eau a également été étudiée par ATVC [88] pour
des pressions de 10-3 hPa et 5 hPa. L'énergie d'activation apparente déterminée à l'aide de la
méthode des créneaux de vitesse est différente pour ces deux pressions. Afin d'interpréter ce
résultat, il est proposé un changement de mécanisme réactionnel avec la pression de vapeur
d'eau : en dessous de 10-3 hPa, la diffusion est l'étape limitante, alors qu'au-dessus de 5 hPa, la
désorption de vapeur d'eau devient l'étape limitante.
Par diffraction des rayons X [70,77,78], il a été observé que la nouvelle phase obtenue,
la métakaolinite est très désordonnée voire amorphe. En fin de transformation, il reste au sein
de la métakaolinite environ 11% des groupements hydroxydes initiaux [89]. De plus, l'étude
par RMN de l'aluminium 27 [70,90] montre l'apparition d'ions aluminium de coordinances IV
et V au cours de la transformation. Cependant, il est important de noter que malgré cette
désorganisation de la structure du solide, la structure en feuillets persiste dans la
métakaolinite [89].
Il ressort de cette étude bibliographique que plusieurs modèles cinétiques de

déshydroxylation ont été proposés selon les conditions expérimentales ou encore selon la
nature de la kaolinite étudiée. Ces modèles sont des cas limites de germination ou de
- 70 -
croissance instantanée, et de plus, il apparaît qu'un modèle unique ne permet pas de décrire la
totalité de la transformation. Nous allons donc essayer de développer un modèle cinétique en
partant d'une étude expérimentale isotherme-isobare avec une kaolinite de Géorgie qui
présente l'avantage d'être une kaolinite pure (à plus de 98%). Puis, nous essaierons d'étendre
ce modèle à des expériences non isothermes.
III.2 Etude expérimentale
Les échantillons utilisés proviennent d'un gisement de Géorgie aux Etats Unis et sont
commercialisés par la firme Sigma.
Les expériences réalisées ont deux objectifs :

• la caractérisation des échantillons ;
• l'obtention des courbes cinétiques expérimentales.
III.2.1 Caractérisation des échantillons
Les techniques expérimentales utilisées afin de décrire la structure et la géométrie des

échantillons sont les suivantes :
• diffraction des rayons X ;
• observations microscopiques ;
• étude de la porosité et de la surface spécifique par adsorption gazeuse ;
• granulométrie.
III.2.1.1 Diffractogrammes de la kaolinite et de la métakaolinite
Les diffractogrammes de la poudre initiale et de la poudre en fin de transformation ont

été obtenus à l'aide d'un diffractomètre Siemens D5000, dont le générateur est un tube
anticathode de cuivre dont la longueur d'onde est celle de la raie Kα du cuivre (λ=1,5418 Å).
Ces diffractogrammes sont présentés sur les figures III.4 et III.5. Les pics théoriques de la
- 71 -
kaolinite (fiche JCPDS 29-1488) sont également représentés sur le diffractogramme de la

figure III.4.
150
140
130
120
110
100
Lin (Counts)
90
80
70
60
50
40
30
20
10
11 20 30 40 50 60
2-Theta - Scale
Figure III.4 : Diffractogramme des rayons X du solide initial et pics théoriques de la kaolinite
27
26
25
24
23
22
21
20
19
18
17
Lin (Counts)
16
15
14
13
12
11
10
9
8
7
6
5
4
3
2
1
0
10 20 30 40 50 60
2-Theta - Scale
Figure III.5 : Diffractogramme des rayons X du solide en fin de transformation
Les pics du diffractogramme du produit initial (figure III.4) montrent que le solide est
bien cristallisé, et correspondent aux pics théoriques de la kaolinite. De plus, la poudre ne
présente pas d’impuretés visibles par diffraction des rayons X.
En fin de transformation, le solide obtenu (figure III.5) correspond à une phase
amorphe ce qui est en accord avec les données bibliographiques sur la métakaolinite.
- 72 -
III.2.1.2 Observations microscopiques
Des observations microscopiques ont été effectuées afin de connaître la forme

géométrique et la taille des grains de la poudre étudiée.
1 Microscopie électronique à transmission
Ces observations ont été réalisées à l'aide d'un appareil Philips CM200 en dispersant la
poudre préalablement dans de l'éthanol. Ainsi, un grain de kaolinite est présenté sur la figure
III.6.
100nm
Figure III.6 : Cliché d'un grain de kaolinite observé par microscopie électronique à transmission
La forme de la face observée sur ce cliché est un hexagone allongé et, d'après la
bibliographie (cf. III.1.1), correspond sûrement à la face basale du grain. On peut donc
supposer que les grains se présentent bien sous la forme de plaquettes hexagonales de fine
épaisseur. Ceci va être vérifié par les observations de la poudre réalisées par microscopie
électronique à balayage.
1 Microscopie électronique à balayage
Des clichés de la poudre ont été obtenus par microscopie électronique à balayage à
effet de champs de type JEOL JSM6500F. Ces micrographies ont montré à la fois la forme et
la taille des faces basales représentées par les clichés de la figure III.7, mais également
l'épaisseur des plaquettes en observant des empilements de plaquettes comme le montre la
figure III.8.
- 73 -
Figure III.7 : Observation des faces basales des grains par microscopie électronique à balayage
Figure III.8 : Observation des faces latérales (épaisseurs) des grains par microscopie électronique à
balayage
Ces micrographies confirment que les cristaux de kaolinite sont des plaquettes
hexagonales. L'épaisseur des plaquettes semble être proche pour tous les grains, environ égale
à 0,05 µm (figure III.8). Par contre, l'aire de la surface basale de ces plaquettes varie
beaucoup selon les grains observés. Ainsi, dans un premier temps, les grains seront considérés
de même taille puis, dans un deuxième temps, une distribution de la taille des surfaces basales
sera prise en compte au niveau du modèle afin de se rapprocher de la réalité.
De plus, nous avons également observé des grains de métakaolinite obtenus après la
déshydroxylation de la kaolinite. Ceux-ci sont représentés sur la figure III.9.
Comme on peut le noter sur cette figure, la taille et la forme des grains ne semblent
pas être modifiées après la déshydroxylation.
- 74 -
Figure III.9 : Observation des grains de métakaolinite par microscopie électronique à balayage
Dans la modélisation, il sera plus simple de représenter les grains par des cylindres de
rayon r0 de faible épaisseur e, comme le montre la figure III.10.
r0
e
Figure III. 10 : Géométrie des grains de kaolinite utilisés pour le modèle
Dans la première partie de la modélisation, nous considèrerons les grains de même

dimension. Le rayon moyen de ces grains doit donc être estimé.
III.2.1.3 Caractérisation par adsorption gazeuse
Des isothermes d'adsorption-désorption d'azote à 77 K ont été réalisées sur la poudre

de départ, dégazée préalablement à 150°C, à l'aide d'un appareil Micromeritics ASAP 2000.
1 Isotherme d'adsorption-désorption d'azote
L'isotherme d'adsorption-désorption d'azote correspondant à la kaolinite est

représentée sur la figure III.11.
- 75 -
80
Volum e adsorbée (cm .g )

-1
70
adsorption
3
60 désorption
50
40
30
20
10
0
0 0.2 0.4 0.6 0.8 1
P/P
0
Figure III.11 : Isotherme d'adsorption-désorption d'azote sur la kaolinite de Géorgie
L'isotherme d'adsorption est de type II et est caractéristique d'une surface non poreuse
ou présentant des macropores.
Ainsi, les grains de kaolinite seront considérés comme non poreux.
1 Mesure de la surface spécifique
L’isotherme d’adsorption étant de type II, l’équation BET est applicable. La surface
spécifique de la kaolinite a donc été calculée à l'aide de la méthode BET et est égale à
18 (±1) m2.g-1.
En considérant les grains cylindriques, tous de même taille et d'épaisseur e égale à
50 nm, il est possible d'estimer un rayon moyen rm de ces cylindres sachant que :
2πrm e + 2πrm2
S sp = (III.1)
ρπrm2 e
où Ssp est la surface spécifique et ρ la masse volumique de la kaolinite (ρ=2,65.106 g.m-3).

En tenant compte de l'erreur sur la mesure de la surface spécifique, nous obtenons, à
partir de l'équation (III.1), une valeur pour le rayon moyen des grains égale à 0,3 (±0,1) µm,
valeur que nous utiliserons lorsque nous considérerons la poudre comme constituée des grains
de même taille.
Remarque : une mesure de la surface spécifique de la métakaolinite a également été réalisée.

Celle-ci est de 19,6 (±1) m2.g-1. Ainsi, la surface spécifique ne varie pas au cours de la réaction.
- 76 -
III.2.1.4 Distribution granulométrique
La détermination de la distribution granulométrique dans le domaine

submicrométrique est difficile. Nous avons tenté de l'obtenir en utilisant un granulomètre
Coulter LS 130 en voie humide dans de l’éthanol après dispersion de la poudre à l’aide
d’ultrasons.
Remarque sur les mesures de distribution granulométrique : deux théories permettent de

relier les mesures optiques à la taille des particules : la théorie de Fraunhofer pour les particules
supérieures à 1 µm et la théorie de Mie qui s’applique aux petites particules. Ces théories optiques ont
été développées en considérant des particules sphériques. L’utilisation la théorie de Mie nécessite la
connaissance de l’indice de réfraction du matériau et du milieu. Nous prendrons la valeur 1,564 pour
celui de la kaolinite [91].
Les courbes granulométriques obtenues avec les deux théories optiques de Fraunhofer
et de Mie sont représentées sur la figure III.12.
5 Mie
4
Fraunhofer
% Vo lum e
0
0.1 1 10
Diam ètre (µ m )
Figure III.12 : Courbes granulométriques expérimentales
La théorie de Fraunhofer ne peut être utilisée pour des grains de tailles inférieures à un
micron et, vu les observations microscopiques, ne s'applique donc pas à la mesure de la taille
des grains de la kaolinite. La théorie de Mie est par contre applicable à la mesure de particules
submicrométriques. D'après la figure III.12, deux populations de grains sont présentes : une
- 77 -
de diamètre supérieur à 1 µm et qui correspond très probablement à des agrégats, et une de

diamètre inférieur à 1 µm attribuable aux grains individuels de la poudre.
Cependant, les théories optiques de détermination de la distribution granulométrique
d'une poudre sont basées sur des calculs où les particules sont sphériques. Or, les grains de
kaolinite ont une géométrie de plaquettes de faible épaisseur à face hexagonale. On peut alors
se poser la question de la validité d'utiliser la théorie de Mie pour une géométrie de grains très
éloignée de la sphère.
Devant cette incertitude, nous avons préféré estimer la distribution granulométrique à
l'aide du rayon moyen obtenu par la méthode BET et des observations microscopiques.
Celles-ci ont montré que l'épaisseur des plaquettes ne varie presque pas selon les grains, ce
qui n'est pas le cas de la surface des faces basales.
Afin d'estimer une distribution du rayon de ces faces, nous avons utilisé une loi
normale en fonction du rayon (équation (I.51)), le rayon moyen choisi est celui trouvé par la
méthode BET. Le rayon des plus gros grains est déterminé à l’aide des observations
microscopiques. Cette distribution des rayons des grains est représentée sur la figure III.13 et
comparée à la distribution obtenue à l'aide de la théorie de Mie (uniquement la population de
grains dont la taille est inférieure à un micron).
Cette première partie de l'étude expérimentale nous a permis de connaître la géométrie

et la taille des grains de la poudre. La deuxième partie présente comment les courbes
cinétiques expérimentales ont été obtenues.
16 Estim ée
14
M ie
12
% Volum e
10
8
6
4
2
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
r ( µ m)
Figure III.13 : Courbe granulométrique estimée des rayons des grains et comparée à celle obtenue à
l'aide du granulomètre
- 78 -
III.2.2 Etude cinétique de la réaction par thermogravimétrie
III.2.2.1 Protocole expérimental
L'étude cinétique expérimentale de la déshydroxylation a été réalisée à l'aide d'une

thermobalance symétrique SETARAM MTB10-8 en atmosphère statique. Celle-ci est décrite
sur la figure III.14.
Jauges Acquisition
pression informatique
de ∆m
Introduction de H2O Groupe

de vide
H2O
Bain
thermorégulé
Echantillon
Four symétrique
mobile
Thermocouple
Creuset de référence
Enceinte
thermostatée
Figure III.14 : Schéma de la thermobalance utilisée pour l'étude de la déshydroxylation de la
kaolinite
Les échantillons de 20 mg de kaolinite sont introduits dans la thermobalance.

L'échantillon est chauffé à 150°C sous vide dynamique afin d'éliminer l'eau adsorbée.
La vapeur d'eau est ensuite introduite dans la thermobalance à la valeur choisie en
fixant la température d'un ballon d'eau thermorégulé. Différentes pressions de vapeur d'eau
ont été étudiées : 2,5 hPa, 3,8 hPa, 5 hPa, 10 hPa. Certaines expériences ont également été
réalisées sous vide dynamique ; on considère alors la pression de vapeur d'eau résiduelle
environ égale à 10-3 hPa. Cette valeur correspond à la pression limite détectable dans notre
système. La vapeur d'eau étant un gaz difficile à éliminer, nous supposons que la quantité de
gaz restant au sein de la thermobalance est alors de la vapeur d'eau. La pression est contrôlée
à l'aide de deux jauges de pression (de type Membranovac MV11052 et Varian WV100-2).
Parallèlement, le four est chauffé en position basse jusqu'à la température d'étude de la
réaction comprise entre 450°C et 500°C. Une fois que la pression de vapeur d'eau est fixée et
- 79 -
que le four a atteint sa température de consigne, celui-ci est remonté rapidement afin de
démarrer la réaction. Un thermocouple placé sous le creuset indique qu'il faut à peu près
300 secondes pour que l'échantillon atteigne la température du four. Le début de la réaction
n'est donc pas en condition isotherme. Cependant, dans la première partie de la
modélisation, nous supposerons que la réaction se déroule en condition isotherme tout au long
de l'expérience. Nous reviendrons par la suite sur la non isothermicité du début de la réaction.
A partir des expériences de thermogravimétrie décrites précédemment, nous obtenons

l'évolution de la masse de l'échantillon au cours du temps.
En laissant l'échantillon pendant une longue période à 450°C, il a été possible de
déterminer la perte de masse totale de la réaction. Celle-ci est de 13,24%. Elle est inférieure à
la perte de masse théorique de cette réaction, 13,98%. Ceci peut d'abord être la conséquence
de la présence d'impuretés (quartz, mica), non détectables par diffraction des rayons X, dans
les échantillons de kaolinite. De plus, comme nous l'avons vu dans l'étude bibliographique de
la réaction, il reste des groupements hydroxydes au sein de la métakaolinite.
Le degré d'avancement de cette réaction est obtenu par la formule :
∆m( t )
α( t ) = (III.2)
∆mt
où ∆m(t) est la perte de masse à un instant t, et ∆mt est la perte de masse totale (calculée sur la
base de 13,24% de la perte de masse théorique).
III.2.2.2 Courbes cinétiques expérimentales
A partir de ce mode opératoire, les courbes cinétiques expérimentales ont été obtenues
pour différentes pressions de vapeur d'eau et différentes températures. Pour toutes ces
expériences, la partie non isotherme en début de réaction ne sera pas représentée.
III.2.2.2.1 Influence de la pression de vapeur d'eau
Pour cette première série d'expériences, la température est fixée à 450°C. Les
expériences réalisées à différentes pressions de vapeur d'eau sont représentées sur la figure
III.15 avec les courbes α(t).
- 80 -
1 2,5hPa
-3
10 hPa 5hPa
0.8
3,8hPa
0.6
10hPa
α
0.4
0.2
0
4 4 4 4 4
0 1 10 2 10 3 10 4 10 5 10
temps (s)
Figure III.15 : Courbes cinétiques expérimentales à 450°C pour différentes pressions de vapeur d'eau
Il ressort de cet ensemble d'expériences que l'augmentation de la pression de vapeur

d'eau ralentit la réaction. Ceci était prévisible vu que la vapeur d'eau est produite par la
réaction et que les conditions expérimentales ne sont pas très éloignées des conditions
d'équilibre thermodynamique (à 450°C la pression d'équilibre est de l'ordre de 25 hPa [71]).
III.2.2.2.2 Influence de la température
De même l'effet de la température a été étudié à l'aide d'une deuxième série

d'expériences. Celles-ci sont effectuées à une pression de vapeur d'eau de 2,5 hPa et sont
représentées par les courbes du degré d'avancement en fonction du temps sur la figure III.16.
1 500°C
490°C
480°C
0.8 470°C 460°C
450°C
0.6
α
0.4
0.2
0
0 1000 2000 3000 4000 5000 6000 7000 8000
temps (s)
Figure III.16 : Courbes cinétiques expérimentales à PH2O=2,5 hPa pour différentes températures
- 81 -
Remarque : comme nous l’avons vu au III.1.2, la deuxième série de données

thermodynamiques (logiciel Thermodata) ne prévoit pas que la réaction se produise dans les
conditions de températures et de pressions choisies, or c’est le cas. Donc, nous pouvons rejeter ces
données. Cependant, ceci ne nous permet pas de conclure de la validité des données de la première
série.
III.2.2.3 Expériences réalisées en programmation de température
Certaines expériences ont également été réalisées en programmation de température.

Le protocole expérimental est similaire à celui décrit précédemment mais au lieu de remonter
rapidement le four pour démarrer la réaction, nous programmons une montée linéaire de la
température.
Un exemple de ce type d'expérience est présenté sur la figure III.17 pour une vitesse
de programmation de température de 3°C.min-1 jusqu'à 500°C suivi d'un palier isotherme. La
pression partielle de vapeur d'eau est fixée à 2,5 hPa dès la température ambiante.
1 520
500
0.8
480
T(°C)
0.6
T(°C)
460
α
0.4 α 440
420
0.2
400
0
4
0 2000 4000 6000 8000 1 10
temps (s)
Figure III. 17 : Courbe cinétique expérimentale pour une montée en température de 3°C.min-1 jusqu'à
500°C (PH2O=2,5 hPa)
Afin de rendre compte des courbes expérimentales obtenues dans ces différentes
conditions, nous allons mettre en place un modèle en utilisant des hypothèses simplificatrices
et en considérant d'abord la forme des courbes obtenues en condition « isotherme ».
- 82 -
III.3 Modélisation de la réaction
La méthodologie utilisée pour le développement d'un modèle cinétique applicable en

conditions non isotherme et non isobare est celle décrite dans le chapitre II.
La première étape est la détermination d'un modèle de transformation.
III.3.1 Modèle de transformation
La construction du modèle repose toujours sur certaines hypothèses qu'il est nécessaire
de vérifier préalablement. Nous verrons successivement :
• la pseudo-stationnarité,
• le test du φE,
• le test du cas limite.
III.3.1.1 Vérification des hypothèses du modèle
III.3.1.1.1 Pseudo-stationnarité de la réaction
L'étude expérimentale de la pseudo-stationnarité du système réactionnel a été effectuée

par couplage thermogravimétrie-calorimétrie avec un appareil TG-DSC111 décrit en annexe
B.
Les expériences de thermogravimétrie couplée à la calorimétrie ont été réalisées pour
des échantillons de 15 mg de kaolinite sous flux d'hélium à 450°C. La vitesse de montée de la
température est de 15°C.min-1 de la température ambiante jusqu'à 400°C, puis de 5°C.min-1
jusqu'à 450°C.
Les courbes calorimétrique et thermogravimétrique en fonction du degré d'avancement
(défini par rapport à la perte de masse) sont tracées sur la figure III.18 en ajustant les deux
échelles (mg.s-1 et mW) de manière à superposer au mieux les deux courbes.
Si le système est pseudo-stationnaire alors la relation entre la vitesse de perte de masse
(dm/dt) et la puissance thermique (dQ/dt) s'exprime par (cf. équation (B.8)) :
- 83 -
dQ ∆H 0 dm
= (III.3)
dt 2M H 2O dt
où ∆H0 est l'enthalpie standard de la réaction et MH2O est la masse molaire de l'eau.
Remarque : le changement de régime de température, c'est-à-dire le passage d'une montée de

5°C.min-1 à un palier à 450°C, entraîne une perturbation sur le signal du calorimètre. Ainsi, la
comparaison entre la vitesse de perte de masse et la puissance thermique n'est possible qu'à partir d'un
degré d'avancement égal à 0,2.
0.0001
vitesse de perte de m asse 0
0
puiss ance therm ique

dm/dt (mg.s )
-4
-1
-1 10 -1
dQ /dt (mW )
-0.0002
-2
-0.0003
-0.0004 -3
-0.0005
-4
-0.0006
0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

α (m/ ∆ m)
Figure III. 18 : Comparaison entre la vitesse de perte de masse et la puissance thermique pour la
déshydroxylation de la kaolinite à 450°C
D'après la figure III.18, il existe bien une affinité orthogonale entre la vitesse de perte
de masse et la puissance thermique pour un degré d'avancement supérieur à 0,2. Donc, le
système est pseudo-stationnaire au moins à partir de ce degré d'avancement.
De plus, à partir de l'équation (III.3) et du rapport entre la vitesse de perte de masse et
la puissance thermique obtenue à partir de la figure III.18, nous avons pu déterminer la valeur
de l'enthalpie standard à 450°C :
∆H0=230 (±15) kJ.mol-1
Cette valeur est inférieure à celle obtenue à l'aide les données thermodynamiques [71]
(∆H0=270 kJ.mol-1).
Afin de vérifier la pseudo-stationnarité sur l'ensemble de la réaction, nous avons

réalisé le couplage calorimétrie-thermogravimétrie pour des expériences réalisées en
- 84 -
programmation en température tout au long de la réaction. Ainsi, le signal relatif à la

puissance thermique n'est pas perturbé en début de réaction.
Afin de vérifier le test de pseudo-stationnarité en programmation de température, nous
supposons que l'enthalpie standard ne varie pas avec la température dans les conditions de
l'expérience. Ainsi, d'après l'équation (III.3), si le système est pseudo-stationnaire, les courbes
de vitesse de perte de masse et de puissance thermique doivent de nouveau présenter une
affinité orthogonale.
La figure III.19 représente ainsi les deux vitesses obtenues par thermogravimétrie et
par calorimétrie en fonction du degré d'avancement pour une programmation en température
de 8°C.min-1 jusqu'à 600°C sous hélium.
T=300°C T=600°C
0.0005 0
0
dQ/dt
-5
dm/dt (mg.s )
-1
-0.0005
dQ/dt (m W )
-0.001 dm/dt
-10
-0.0015
-0.002
-15
-0.0025
-0.003 -20
-0.0035
0 0.2 0.4 0.6 0.8 1
α (= ∆ m/m)
Figure III.19 : Comparaison de la vitesse de perte de masse et de la puissance thermique pour une
programmation de température de 8°C.min-1
Les courbes, représentatives de la vitesse de perte de masse et de la puissance

thermique, présentent une affinité orthogonale entre elles ; par conséquent, nous pouvons
considérer que le système est pseudo-stationnaire sur tout le domaine de degré
d'avancement et pour des températures comprises entre 300°C et 600°C.
III.3.1.1.2 Test du φE
La vérification expérimentale du test du φE est réalisée par thermogravimétrie à l'aide

du montage décrit sur la figure III.14. Nous utilisons la méthode développée dans l'annexe C.
- 85 -
Ainsi, une série d'expériences est effectuée avec au départ les mêmes conditions de
température, 450°C, et de pression de vapeur d'eau, 2,5 hPa. A un certain instant de la
réaction, la pression est modifiée rapidement à l'aide du bain thermostaté pour obtenir une
pression de vapeur d'eau égale à 10 hPa. Quelques exemples de ces décrochements sont
présentés sur la figure III.20.
Remarque : la durée nécessaire à un changement de pression est de l'ordre de 200 secondes

pour passer de 2,5 à 10 hPa, ce qui est faible devant la durée de l'expérience (>6000 s).
10hPa
0.8
10hPa
10hPa
0.6
α
2,5hPa
0.4 10hPa
0.2
0
0 1000 2000 3000 4000 5000 6000 7000
temps (s)
Figure III.20 : Courbes du degré d'avancement en fonction du temps pour des décrochements en
pression de vapeur d'eau de 2,5 à 10 hPa
Pour chaque expérience, les vitesses, avant (vg) et après (vd) le décrochement, sont
mesurées, et les rapports entre les deux vitesses sont calculés. Ces résultats sont présentés
dans le tableau III.1.
α 0,25 0,5 0,57 0,75 0,84 0,9

-1
vd (s ) 1,67.10-4 4,4510-5 4,26.10-5 1,82.10-5 8,87.10-6 2,07.10-6
vg (s-1) 6,42.10-4 2,11.10-4 1,81.10-4 7,09.10-5 4,43.10-5 1,01.10-5
vd/vg (±0,04) 0,26 0,21 0,23 0,25 0,20 0,20
Tableau III.1 : Rapport des vitesses à partir de la méthode des décrochements en pression de 2,5 hPa
à 10 hPa
A partir de ce tableau et vu l'erreur estimée sur les mesures, il est possible de conclure
que le rapport des vitesses ne varie pas quel que soit l'instant du décrochement. Ainsi, la
α/dt=φ
vitesse peut s'écrire sous la forme dα φE.
- 86 -
III.3.1.1.3 Test du cas limite
Ce test, décrit dans le premier chapitre, permet de savoir si la réaction correspond ou

non à un cas complexe de germination-croissance en observant l'influence du passé sur la
réaction. Ainsi, deux expériences ont été réalisées par thermogravimétrie avec la
thermobalance MTB10-8. La première est effectuée selon le protocole expérimental décrit
pour obtenir les courbes cinétiques expérimentales en condition « isotherme » pour une
pression de vapeur d'eau 2,5 hPa et une température de 460°C. Pour la deuxième expérience,
de nouveau, le même mode opératoire est appliqué au départ, avec une pression de vapeur
d'eau de 2,5 hPa, mais cette fois-ci une température de 450°C. Au bout d'un certain temps,
cette température est modifiée pour obtenir une température de 460°C.
Les vitesses de ces deux expériences sont comparées sur la figure III.21.
0.0007
T=460°C
0.0006
0.0005
d α /dt (s )
-1
T=450°C
0.0004
T=460°C
0.0003
0.0002
0.0001
0
0 0.2 0.4 0.6 0.8 1
α
Figure III.21 : Test du cas limite (PH2O=2,5 hPa)
La vitesse est différente après le décrochement pour les deux expériences (qui sont
alors dans les mêmes conditions expérimentales). La vitesse dépend donc du passé de la
réaction et nous en déduisons que le modèle de transformation est un cas complexe de
germination-croissance (et non un modèle limite comme on l'a vu dans l'étude
bibliographique).
Une fois ces trois tests réalisés, il est possible de choisir un modèle de transformation
en tenant compte de la forme des courbes cinétiques expérimentales obtenues dans les
conditions "isotherme" et isobare.
- 87 -
III.3.1.2 Choix du modèle
La figure III.22 décrit la courbe de la vitesse en fonction du degré d'avancement, pour

une expérience isotherme à 480°C et avec une pression de vapeur d'eau fixée à 2,5 hPa. Une
nouvelle fois, la partie non isotherme en début de réaction n'est pas représentée.
0.0012
0.001
0.0008
d α /dt (s )
-1
0.0006
0.0004
0.0002
0
0 0.2 0.4 0.6 0.8 1
α
Figure III.22 : Courbe cinétique dα/dt en fonction de α pour T=480°C et PH2O=2,5 hPa
Deux éléments de cette courbe vont être utiles pour construire le modèle :
• la présence d'un maximum en accord avec une compétition entre la germination et la
croissance (cf. test du cas limite au III.3.1.1.3),
• la forme de la courbe caractéristique d'une étape limitante de diffusion. Ce résultat est en
accord avec la bibliographie pour les modèles de type D2 (cf. III.1.3).
Compte tenu de la nature même de la transformation (déshydroxylation), le sens du
développement de la nouvelle phase, la kaolinite, ne peut être qu'interne.
Nous ne connaissons pas le volume molaire de la métakaolinite. Donc, il n'a pas été
possible de calculer le coefficient d'expansion volumique Z pour cette réaction. Cependant,
étant donné que les observations microscopiques montrent que la taille des grains semble être
la même au début et en fin de réaction et que la surface spécifique ne varie pas non plus au
cours de la déshydroxylation, nous supposons que la taille des grains reste constante durant la
réaction, et donc que le coefficient d'expansion volumique est égal à 1.
Pour interpréter les courbes cinétiques expérimentales, nous proposons un modèle de
germination-croissance anisotrope. Ainsi que nous l'avions remarqué au I.1.3 (tableau I.2),
seule la croissance anisotrope est compatible avec une étape de diffusion.
- 88 -
L'étape limitant la croissance est une étape de diffusion au sein de la couche de

métakaolinite formée. La diffusion se déroule vraisemblablement dans la direction des
feuillets, parallèlement aux faces basales, de manière prépondérante.
En tenant compte des observations microscopiques, les grains sont assimilés à des
cylindres non poreux de rayon r0=0,3 µm et d'épaisseur e=0,05 µm. Le schéma d'un grain en
cours de transformation est présenté sur la figure III.23.
Métakaolinite
Kaolinite
Figure III. 23 : Schéma d'un grain en cours de transformation
Les figures III.24.a et III.24.b décrivent la poudre au cours de la déshydroxylation de

la kaolinite. La germination a lieu à la surface des grains avec une fréquence surfacique de
germination γ constante. Le germe apparaissant sur un grain entraîne alors quasi
instantanément la déshydroxylation de toute la surface du grain. En effet, la vitesse de
croissance tangentielle de ces germes est supposée très grande devant la vitesse de croissance
radiale (croissance anisotrope). Les germes n'apparaissent pas tous au même moment (figure
III.24.a). Comme tous les germes ne se forment pas au même instant, en cours de
transformation, il est à la fois possible d'avoir des grains de kaolinite non transformés, des
grains en cours de réaction et des grains qui ont complètement réagi (figure III.24.b).
(a) Kaolinite (b)
Kaolinite
Métakaolinite
Métakaolinite (germe)
Figure III. 24 : Poudre de kaolinite en cours de déshydroxylation
A partir des hypothèses de ce modèle, il est possible de calculer les courbes théoriques
du degré d'avancement et de la vitesse en fonction du temps.
- 89 -
III.3.1.3 Expression de la vitesse
III.3.1.3.1 Avec une seule taille de grain
Le calcul de la vitesse pour un modèle de germination-croissance anisotrope à

développement interne avec une étape limitante de diffusion pour la croissance, et des grains
cylindriques de même taille, est développé dans l’annexe E (cf. équation E.21). Il conduit à
l’expression de la vitesse en conditions isotherme et isobare suivante :
dα =VmKφγs0 exp(−γs τ)
t
−4
dt r0 ∫ 0
ln(1− β(t,τ))
dτ (III.4)
0
où VmK est le volume molaire de la kaolinite. β(t,τ) est le degré d’avancement d’un grain né à
τ et est solution de l'équation :
4VmKφ
β +(1−β)ln(1−β)= (t −τ) (III.5)
r0
En résolvant le système d'équations (III.4) et (III.5), par l'utilisation d'un calcul

numérique réalisé sous MATLAB 5.3 et développé en collaboration avec Eric Touboul du
centre SIMMO de l'Ecole des Mines de Saint Etienne, il est possible d'obtenir la vitesse dα
dt
en fonction du temps et, par intégration de celle-ci, l’expression du degré d’avancement.
Remarque : l'équation (III.4) correspond à l'écriture dα/dt=φE avec une fonction d'espace
égale à :
Vm γs t −4
E(t)= K 0 ∫exp(−γs0τ) dτ (III.6)
r0 0 ln(1− β(t,τ))
III.3.1.3.2 Avec une distribution granulométrique
L’introduction d’une distribution granulométrique dans le calcul de la vitesse est

simple à réaliser. A partir du calcul précédent, nous pouvons déterminer l’expression du degré
d’avancement pour une taille de grain r. Notons celui ci α(t,r). Il est possible de réaliser ce
- 90 -
calcul pour n’importe quel rayon et d’en déduire l’expression du degré d’avancement pour
une poudre de grains de tailles différentes (cf. équation I.66).
La vitesse peut s’écrire sous la forme :
dα( t ) rmax dα( t , r )
=∫ g( r )d ( r ) (III.7)
dt rmin dt
rmin est le rayon des plus petits grains et rmax le rayon des plus gros.
La distribution granulométrique des rayons des particules (g(r)) a été présentée sur la
figure III.13.
De nouveau, l’ensemble des calculs a été réalisé à l’aide du logiciel Matlab 5.3.
III.3.1.4 Comparaison avec l'expérience
III.3.1.4.1 Avec une seule taille de grain
Les courbes théoriques obtenues à partir des équations (III.4) et (III.5) dépendent de
deux paramètres inconnus : γ et φ. La détermination de ces paramètres est réalisée à l’aide de
la méthode présentée dans l'annexe D en utilisant les vitesses réduites et le paramètre de
modèle A. Les équations de la vitesse réduite sont présentées dans l’annexe E. De plus, dans
cette annexe, figure également un exemple de comparaison en vitesse réduite d'une courbe
expérimentale et d'une courbe théorique, avec la recherche du paramètre A, ainsi que la
détermination de l'erreur sur A, φ et γ. Comme le montre le tableau III.2, le paramètre de
modèle (proportionnel au rapport γ/φ) augmente avec la pression de vapeur d’eau.
PH2O (hPa) 10-3 2,5 3,8 5 10

A 15±2 25±5 45±10 55±15 600±? *
Tableau III. 2 : Valeur de A en fonction de la pression de vapeur d’eau
Ainsi, nous avons pu déterminer les valeurs de φ et γ par ajustement entre les courbes
calculées et les courbes expérimentales. La confrontation entre ces courbes de vitesse en
fonction du degré d'avancement pour différentes pressions de vapeur d’eau, et pour une
température de 450°C, est présentée sur la figure III.25.
- 91 -
Remarque 1 : seule la partie isotherme des courbes expérimentales est représentée sur les
figures qui montrent la comparaison entre le modèle et l’expérience.
Remarque 2 : l'erreur sur le paramètre de modèle de l'expérience à 10 hPa n'a pu être

déterminée. En effet, pour cette expérience nous n'avons pas accès à la position du maximum de
vitesse (voir figure III.25). Cette expérience est proche du cas limite de germination instantanée (γ
grand). Les différences entre les courbes de vitesse pour les valeurs de fréquence surfacique de
germination élevées interviennent au début de celles-ci. Or, nous ne pouvons pas considérer ce début
puisqu'il se déroule en condition non isotherme. Donc, il n'est pas possible de déterminer une valeur
précise de A.
La valeur très élevée de γ à 10 hPa de vapeur d'eau n'est pas un artefact, même s'il existe une
incertitude importante sur sa détermination. En effet, nous avons tenté de réaliser des expériences en
établissant une pression supérieure à la pression d'équilibre, l'échantillon étant déjà à la température
voulue pour la réaction. Dans ces conditions, après diminution de la pression de vapeur d'eau à la
valeur de l'expérience inférieure à la pression d'équilibre, nous avons observé que la forme des courbes
était différente, en particulier le maximum de vitesse avait disparu. Cette modification de la forme des
courbes traduit en fait le passage à un modèle limite de transformation, dans lequel la germination est
instantanée, en accord avec une très forte augmentation de γ quand la pression de vapeur d'eau
augmente.
0.0006 -3
10 hPa
Expériences
0.0005
Modèle
d α /dt (s )
2,5 hPa
-1
0.0004
3,8 hPa
0.0003
5 hPa
0.0002
0.0001 10 hPa
0
0 0.2 0.4 0.6 0.8 1
α
Figure III.25 : Comparaison des courbes cinétiques expérimentales et théoriques pour différentes
pressions de vapeur d’eau (T=450°C)
Pour les pressions de vapeur d’eau comprises entre 10-3 et 5 hPa, le modèle représente
assez bien la vitesse expérimentale sauf pour la fin des courbes. On peut penser que ces écarts
- 92 -
proviennent de la granularité de la poudre qui n’a pas encore été prise en compte, et peut être
aussi du début non isotherme des expériences.
Pour l’expérience réalisée à 10 hPa, le modèle n’interprète pas bien le début de
l’expérience. Cette expérience est proche du cas limite de germination instantanée.
Comme le montre la figure III.26, les courbes théoriques ont également été comparées
aux expériences réalisées à différentes températures pour une même pression de vapeur d’eau
(2,5 hPa).
0.0024
Expériences
0.002 500°C Modèle
0.0016
d α /dt (s )
-1
490°C
0.0012 480°C
470°C
0.0008
460°C
0.0004
450°C
0
0 0.2 0.4 0.6 0.8 1
α
Figure III. 26 : Comparaison des courbes cinétiques expérimentales et théoriques pour différentes
températures (PH2O=2,5 hPa)
L’ajustement a été plus délicat. En effet, plus les expériences se déroulent à haute
température, plus une partie importante de la réaction se déroule en condition non isotherme.
Par exemple à 500°C, la partie isotherme de l’expérience débute seulement à partir d’un degré
d’avancement égal à 0,2. Ceci explique pourquoi le modèle ne décrit pas bien l’expérience à
490°C et 500°C : le modèle considère que la réaction se déroule en condition isotherme, or
avec ces températures, on s’éloigne beaucoup de ces conditions.
Pour les courbes obtenues à une température inférieure, nous pouvons faire les mêmes
remarques que précédemment à savoir que le modèle décrit assez bien l’expérience, mais qu'il
existe des écarts en fin d’expérience.
Afin d’essayer de comprendre cette différence entre les courbes théoriques et
expérimentales en fin de réaction, nous avons considéré une distribution granulométrique
dans le calcul de la vitesse.
- 93 -
III.3.1.4.2 Avec une distribution granulométrique
Les valeurs de φ et γ, nécessaires au calcul, sont celles obtenues précédemment pour

une population de grains de même taille.
Les figures III.27 et III.28 présentent la confrontation entre les courbes calculées et
les courbes expérimentales pour les différentes pressions de vapeur d’eau à 450°C (figure
III.27) et pour les différentes températures avec une pression de 2,5 hPa (figure III.28).
0.0006 -3
10 hPa
Expériences
0.0005
Modèle
2,5hPa
0.0004
d α /dt (s )
-1
3,8hPa
0.0003
5hPa
0.0002
0.0001
10hPa
0
0 0.2 0.4 0.6 0.8 1
α
pressions de vapeur d’eau (T=450°C) en considérant une distribution granulométrique
0.0025
Expériences
Modèle
0.002 500°C
d α /dt (s )
-1
0.0015 490°C
480°C
0.001
470°C
460°C
0.0005
450°C
0
0 0.2 0.4 0.6 0.8 1
α
températures (PH2O=2,5 hPa) en considérant une distribution granulométrique
- 94 -
Comme le montrent ces deux figures, la prise en compte des différentes tailles des
grains dans le calcul de la vitesse, donne un meilleur accord entre les courbes calculées et les
courbes expérimentales, en particulier en ce qui concerne la fin de la réaction.
Cependant, pour des températures élevées où la réaction se déroule rapidement, une
nouvelle fois le modèle n’interprète pas bien la courbe cinétique expérimentale. Ceci est la
conséquence du traitement des courbes en condition isotherme alors qu’une partie a lieu à
température non constante. Il est donc nécessaire de prendre en compte dans le modèle, la
variation initiale de température.
Mais pour cela, il faut connaître l’évolution de la réactivité surfacique de croissance et
celle de la fréquence surfacique de germination en fonction de la température. Afin d'obtenir
les lois correspondantes, nous présentons les résultats de la modélisation physico-chimique.
III.3.2. Modélisation physico-chimique
III.3.2.1 Variation de φ et γ en fonction des conditions physico-

chimiques
III.3.2.1.1 Par comparaison du modèle avec les expériences
L’ajustement des courbes théoriques sur les courbes expérimentales nous a permis
d’obtenir les valeurs de la fréquence surfacique de germination, et de la réactivité surfacique
de croissance en fonction de la pression de vapeur d’eau, et de la température. Ces valeurs
sont présentées dans les figures III.29 et III.30.
La réactivité surfacique de croissance diminue avec la pression de vapeur d’eau, alors

que la fréquence surfacique de germination est constante pour les basses pressions de vapeur
d’eau (PH2O<5 hPa), puis augmente pour la pression la plus forte (PH2O=10 hPa). Pour cette
pression, nous rappelons que le modèle ne permet pas de bien interpréter le début de
l’expérience.
- 95 -
-7 11
1 10 2.5 10
-8 11 γ
γ (germ e.m .s )
8 10 φ 2 10
-1
φ (mole.m .s )
-1
-2
-8
-2
11
6 10 1.5 10
-8 11
4 10 1 10
-8 10
2 10 5 10
0 0
0 2 4 6 8 10
0 2 4 6 8 10
P (hPa)
P (hPa) H2O
H2O
Figure III.29 : Valeurs de φ et γ en en fonction de la pression de vapeur d’eau (T=450°C)
11
-7 2 10
7 10
-7
6 10 11
γ (germes.m .s )
φ -1
1.6 10
φ (mole.m .s )
-1
5 10
-7 γ
-2
-2
11
-7 1.2 10
4 10
-7
3 10 10
8 10
-7
2 10
-7 10
1 10 4 10
0
450 460 470 480 490 500 450 460 470 480 490 500
T(°C) T(°C)
Figure III.30 : Valeurs de φ et γ en en fonction de la température
La réactivité surfacique de croissance et la fréquence surfacique de germination

augmentent avec la température. Cependant, les valeurs trouvées pour les températures les
plus élevés (490°C et 500°C) ne sont pas très précises puisque, pour ces températures,
considérer la réaction comme se déroulant à température constante ne semble plus du tout
valable. D’ailleurs, le modèle développé ne permet pas d’interpréter correctement ces courbes
expérimentales.
III.3.2.1.2 Par la méthode des décrochements
La méthode des décrochements (ou de l'isolement, cf. I.1.4.2.1) a également été

utilisée dans la recherche des variations de la réactivité surfacique de croissance en fonction
de la pression de vapeur d'eau (φ(PH2O)) et de la température (φ(T)).
- 96 -
1 φ(PH2O)
Une série d'expériences a été réalisée par thermogravimétrie, en débutant celles-ci à

450°C sous 2,5 hPa de vapeur d'eau. Puis, au même degré d'avancement (α=0,3), la pression
de vapeur d'eau est modifiée en changeant la température du bain thermorégulé, et fixée à une
nouvelle valeur. Ces expériences sont présentées sur la figure III.31.
0.6 3,8hPa
5hPa
0.5
8hPa 10hPa
0.4
α
0.3
0.2
2,5hPa
0.1
0
0 1000 2000 3000 4000 5000 6000
t (s)
Figure III. 31 : Décrochements en pression de vapeur d'eau (T=450°C)
Le rapport des vitesses, avant et après le décrochement, permet d'accéder aux valeurs
de φ(PH2O)/φ(2,5 hPa) à 450°C présentées sur la figure III.32.
2.5
2
)/ φ (2.5hPa)
1.5
H2O
1
φ( P
0.5
0
0 2 4 6 8 10
P (hPa)
H2O
Figure III. 32 : Evolution de φ en fonction de la pression de vapeur d'eau par la méthode des
décrochements (T=450°C)
- 97 -
1 φ(T)
De nouveau, des expériences sont commencées dans les mêmes conditions de

température, 450°C, et de pression de vapeur d'eau, 2,5 hPa. Cette fois ci, comme le montre la
figure III.33, la température est modifiée rapidement pour être fixée à une nouvelle valeur.
1 500°C
490°C
0.8 480°C
470°C
0.6
460°C
α
0.4
0.2
450°C
0
0 1000 2000 3000 4000 5000
temps (s)
Figure III.33 : Décrochements en température (PH2O=2,5 hPa)
Le rapport des vitesses conduit à la courbe φ(T)/φ(450°C) en fonction de la

température pour une pression de vapeur d'eau de 2,5 hPa tracée sur la figure III.34.
14
12
φ (T)/ φ (450°C)
10
0
450 460 470 480 490 500
T(°C)
Figure III. 34 : Evolution de φ en fonction de la température par la méthode des décrochements
(PH2O=2,5 hPa)
- 98 -
Les variations de la réactivité surfacique de croissance en fonction des conditions

expérimentales ont été obtenues de deux façons différentes : d'une part par comparaison entre
le modèle de transformation et l'expérience, et d'autre part directement par l'expérience à
l'aide de la méthode des décrochements. Voyons à présent si ces deux méthodes donnent les
mêmes résultats, ce qui validera le modèle de germination-croissance utilisé.
III.3.2.1.3 Validation du modèle de transformation
La comparaison des valeurs de φ, déterminées par les deux méthodes, permet de

prouver la validité du modèle développé puisque la méthode des décrochements est purement
expérimentale. Ainsi, la figure III.35 représente, pour les différentes pressions de vapeurs
d’eau, le tracé du rapport φ(PH2O)/ φ(2,5 hPa) obtenu par la méthode des décrochements en
fonction de φ(PH2O) obtenue à l’aide du modèle de transformation. Il existe une relation
linéaire dont la pente est l’inverse de la réactivité surfacique de croissance à 2,5 hPa.
2.5
-3
10 hPa
)/ φ (2,5hPa)
1.5
1 2,5hPa
H2O
φ (P
3,8hPa
0.5
10hPa 5hPa
0
-8 -8 -8 -8 -7
0 2 10 4 10 6 10-2 -1 8 10 1 10
φ (P ) (mole.m .s )
H2O
(modèle de transformation)
Figure III.35 : Comparaison des valeurs de la réactivité surfacique de croissance obtenues par le
modèle de transformation et par la méthode des décrochements pour différentes pressions de vapeur
d’eau
De même les valeurs de la réactivité surfacique de croissance ont été comparées pour
différentes températures, sur la figure III.36, en traçant φ(T)/φ(450°C) obtenu par la méthode
des décrochements en fonction de φ(T) obtenue par confrontation entre les expériences et le
modèle de transformation. Comme nous l’avons évoqué, les valeurs de la réactivité surfacique
- 99 -
de croissance à 490°C et 500°C sont très approximatives puisque ces expériences doivent être
modélisées en condition non isotherme. Ainsi, nous comparerons uniquement les valeurs
obtenues pour des températures inférieures.
5 4 80°C
4
φ (T)/ φ (450°C)
470°C
3
4 60°C
2
450°C
0
-8 -7 -7 -7 -7 -7
0 5 10 1 10 1.5 10 2 10 2.5 10 3 10
-2 -1
φ (T) (mole.m .s )
(m odèle de transformation)
Figure III.36 : Comparaison des valeurs de la réactivité surfacique de croissance obtenues par le
modèle de transformation et par la méthode des décrochements pour différentes températures
De nouveau, il existe une relation linéaire entre les valeurs de la réactivité surfacique
de croissance obtenues par les deux méthodes pour différentes températures.
Ce résultat valide le modèle de germination-croissance anisotrope et les valeurs de la
réactivité surfacique de croissance obtenues à l'aide de ce modèle.
Essayons de trouver à présent un mécanisme de croissance qui permette d’interpréter
l'évolution de cette grandeur en fonction des paramètres physico-chimiques.
III.3.2.2 Mécanisme réactionnel de croissance
Afin d'écrire un mécanisme réactionnel relatif à la croissance de la couche de

métakaolinite, les notations de Kröger [69] ont été utilisées. De plus les indices K ou MK
indiquent la phase dans laquelle se trouve l'espèce considérée, kaolinite ou métakaolinite
respectivement, et int ou ext l'interface où se situe cette espèce, interne ou externe
respectivement. Les sites d'adsorption sont notés s et les molécules d'eau adsorbées H2O-s.
Le mécanisme de croissance proposé comprend six étapes élémentaires :
- 100 -
1− OHOHK → OH'i MKint +VOH

•
int ← K int
2− OHOHKint → •
← OOH Kint +H iMKint
'
3− OH'i MKint → OH'i MKext

H•i → H•i
MK int MK ext
4− OH'i MKext +H•i +s →

← H2O−s
MK ext
5− H2O−s →
← H2O(g) +s
•
6− O'OHK +VOH → O×
K
int int ← OMKint
Les deux premières étapes se déroulent à l'interface interne (entre la kaolinite et la
métakaolinite), et correspondent à la création de défauts : ions hydroxydes en positions
interstitielles et protons dans la métakaolinite, lacunes d'hydroxydes et ions oxygène en
positions d'hydroxydes dans la kaolinite. La troisième étape est la diffusion des défauts au
sein de la métakaolinite, de l'interface interne vers l'interface externe. L'étape (4) a lieu à
l'interface externe et correspond à la formation d’une molécule d'eau qui est désorbée pendant
l'étape (5). L'étape (6) traduit l'avancement de l'interface interne avec la création de la phase
métakaolinite.
D'après le modèle de transformation, l'étape limitant la croissance est une étape de
diffusion à travers la couche de métakaolinite formée. Il s'agit de l'étape (3). Deux défauts
diffusent : les ions hydroxydes et les protons. Celui qui possède le plus faible coefficient de
diffusion est vraisemblablement l'ion hydroxyde en position interstitielle. En supposant
l'électroneutralité de la métakaolinite, les protons vont diffuser à la même vitesse que les ions
hydroxydes.
Les étapes (1), (2), (4), (5) et (6) sont supposées être à l'équilibre et on note Ki la
constante d'équilibre de l'étape i. Ainsi, en considérant l'étape de diffusion des ions
hydroxydes comme étape limitante, l'expression de la réactivité surfacique de croissance peut
s'écrire selon l'équation (III.8) :
φ = D ([OHi'MK int ] −[OHi'MKext ]) (III.8)

l0
où D est le coefficient de diffusion des hydroxydes au sein de la métakaolinite, et comme
nous l'avons vu au premier chapitre (équation I.25), l0 est une longueur arbitraire nécessaire
pour avoir la réactivité surfacique de croissance en mol.m-2.s-1, nous avons choisi ici l0 égal au
rayon moyen des grains (l0=r0=0,3 µm).
- 101 -
Les étapes (1), (2), (4), (5) et (6) étant à l'équilibre, on peut écrire les équations :
•
K1 =[VOH K int
][OH 'iMK int ] (III.9)
K2 =[O'OH K int ][H i•MK int ] (III.10)
[H 2O−s]
K4 = •
(III.11)
[H iMKext ][OH 'iMKext ][s]
[s]PH 2O
K5 = (III.12)
[H 2O−s]
K6 = 1 (III.13)
•
[O'OH K int ][VOH K int
]
L'électroneutralité dans chaque phase et à chaque interface permet d'écrire :

•
[VOH K int
] = [O'OH K int ] (III.14)
[OH 'iMK int ] = [H i•MK int ] (III.15)
[OH 'iMKext ] = [H i•MKext ] (III.16)
Les équations (III.9) à (III.16) permettent de déterminer la concentration des différents

intermédiaires aux deux interfaces, et en particulier celles en ions hydroxydes. Ainsi,
l'équation (III.8) peut s'écrire :
 1/2
PH 2O 
φ= D 
K1K 2 K6 1−  (III.17)
r0  K1K 2 K 4 K 5 K 6 
 
où PH2O est la pression de vapeur d'eau divisée par la pression de référence (1 bar).
Le produit des constantes d’équilibre K1K2K4K5K6 est égal à la pression d’équilibre,
d'après le mécanisme écrit précédemment dont la réaction bilan met en jeu une molécule
d'eau :
×
2OHOHKint =OOMK +H2O(g)
int
Nous noterons Peq la pression de vapeur d’eau à l’équilibre. L’équation (III.17) peut
donc également s’écrire :
 P 1/2 
φ=D K 1 K 2 K 6  1 − H1/2
2O  (III.18)
r0  Peq 
 
- 102 -
La pression d’équilibre peut être déterminée à l’aide du diagramme thermodynamique

de la figure III.2. Cependant, comme nous l’avons vu, ces données thermodynamiques ne sont
pas connues de manière très précise. Donc, nous allons déterminer la pression d’équilibre de
la réaction en confrontant la loi (III.18) aux valeurs de la réactivité surfacique de croissance
obtenues pour différentes pressions de vapeur d'eau.
1 Comparaison avec les valeurs expérimentales de la réactivité surfacique de

croissance
• φ(PH2O) :
L'équation (III.18) peut être confrontée aux valeurs expérimentales obtenues, soit par
comparaison avec l’expérience, soit par la méthode des décrochements.
Ainsi, la figure III.37 présente l’ajustement de la loi précédente par rapport aux valeurs
de φ en fonction de la pression de vapeur d’eau pour une température de 450°C.
-7
1 10 T=450°C
-8
8 10
φ (mole.m .s )
-1
1/2
P 2O
φ = D K 1K 2 K 6 (1− H1/2 )
-2
-8 r0 P
6 10 eq
-8
4 10
-8
2 10
0
0 2 4 6 8 10
P (hPa)
H2O
Figure III.37 : Comparaison des valeurs expérimentales de φ à différentes pressions de vapeur d’eau
avec la loi déterminée à l’aide du mécanisme réactionnel de croissance
La loi, obtenue à l’aide du mécanisme réactionnel, interprète bien l’évolution de la

réactivité surfacique de croissance en fonction de la pression de vapeur d’eau. L’ajustement
de cette loi est réalisé à l'aide de deux paramètres, et donc à 450°C :
D K K K =9,2.10−8 (±3.10−8) mole.m−2.s−1
r0 1 2 6
- 103 -
Peq =10,5 (±0,7) hPa

Cette valeur de la pression d'équilibre est légèrement inférieure à celle des données
thermodynamiques (Peq=25 hPa) [71], mais elle est très supérieure à celle de Thermodata
(Peq=5.10-4 hPa, cf. figure III.3).
Remarque : à chaque fois que nous ajustons les paramètres d'une loi par rapport à des valeurs
expérimentales, nous déterminons une erreur sur ces paramètres. Celle-ci est calculée à l'aide d'un
programme d'ajustement du logiciel Kaleidagraph basé sur une méthode des moindres carrées . Cette
erreur ne prend pas en compte les barres d'erreur obtenues sur les valeurs expérimentales.
• φ(T) :
Dans la modélisation en condition non isotherme, nous avons besoin de connaître

l’expression de la réactivité surfacique de croissance en fonction de la température.
Dans l’équation (III.18), le coefficient de diffusion dépend de la température et peut
s’écrire :
D=D0 exp(−Ea RT)
où D0 est le facteur pré-exponentiel et Ea est l’énergie d’activation.
De même, les constantes d’équilibres dépendent de la température selon la loi de Van’t
Hoff :
Ki =Ki0 exp(−∆Hi0 RT)
où ∆Hi0 est l’enthalpie standard de l’étape i, et Ki0 est le terme entropique :
Ki0 =exp(−∆Si0 R) .
Ainsi, l’équation (III.18) peut être écrite en fonction de la température :
 E a + ∆ H 10 2 + ∆ H 20 2 + ∆ H 60 2   PH 2 O 
1/2
D0 K 10 K 20 K 60
φ= exp  -  1− 1/2 (III.19)
r0  RT   Peq ( T ) 
Dans l'équation (III.19), la pression d'équilibre intervient. Il est donc nécessaire de

connaître l'expression de celle-ci en fonction de la température. La pression d’équilibre
dépend également de la température selon :
Peq =(exp( ∆S 0 R)exp(−∆H 0 RT) )1 2 (III.20)
où ∆H0 et ∆S0 sont respectivement l’enthalpie standard et l’entropie standard de la réaction.
- 104 -
La valeur de ∆H0 a été déterminée expérimentalement (cf. III.3.1.1.1) et est égale à

230 kJ.mol-1. Il reste donc à déterminée ∆S0. D'après l'équation (III.20), l’entropie standard de
la réaction est égale à :
∆S 0 =2Rln(Peq(T))+∆H 0 / T (III.21)
Or, nous connaissons la valeur de la pression d'équilibre à 450°C. Celle-ci est égale à
10,4 hPa. D'après l'équation (III.21), nous pouvons en déduire l'entropie standard :
∆S0=242 (± 35) J.K-1.mole-1
Ainsi, à partir de ces éléments, nous connaissons la variation de la pression à

l’équilibre en fonction de la température (Peq(T)) et nous l'introduirons dans les calculs.
Il ressort, de l'expression de φ (cf. équation III.17), que la loi d'Arrhenius n'est pas
suivie pour la réactivité surfacique de croissance à cause du terme exprimant l'écart à
l'équilibre.
Afin de déterminer les paramètres cinétiques relatifs à la température, nous ajustons,

sur la figure III.38, la loi (III.19) sur les valeurs expérimentales obtenues à une pression de
vapeur d’eau de 2,5 hPa.
-7
3.5 10
 E + ∆H10 2+ ∆H 02 2+ ∆H 06 2   PH 2O 
1/2
D0 K10K 02K 06
exp - × 1− 1/2
-7
3 10 φ=   
r0  RT   Peq (T) 
φ (mole.m .s )
-1
-7
2.5 10
-2
-7
2 10
-7
1.5 10
-7
1 10
-8
5 10
0
720 725 730 735 740 745 750 755
T(K)
Figure III.38 : Comparaison des valeurs expérimentales de φ à différentes températures avec la loi
déterminée à l’aide du modèle physico-chimique
La loi (III.19) est en bon accord avec les valeurs expérimentales de φ en fonction de la
température. L’ajustement se fait au moyen de deux paramètres, et en conséquence :
- 105 -
D0 K10K02K06
3,2.107 ≤ ≤1,1.109 mole.m−2.s−1
r0
Ea +∆H10 2+∆H02 2+∆H06 2=215 (±15) kJ.mole−1
Remarque : l'erreur sur le paramètre D0 K10K 02K 06 r0 est très grande. Ceci peut paraître
surprenant. En fait, ce résultat est dû au fait que les erreurs sur les deux paramètres de la loi (III.19)
sont liées. Prenons un exemple pour justifier ce propos. Nous notons Ac le paramètre D0 K10K 02K 06 r0
et Eac le paramètre Ea +∆H10 2+∆H02 2+∆H06 2 . Prenons deux couples de paramètres (Ac1,Eac1) et
(Ac2,Eac2) qui permettent d'obtenir un ajustement valable de la loi (III.19) par rapport aux valeurs
expérimentales. Alors, il est possible d'écrire :
 PH 2 O 
1/2
 PH 2 O 
1/2
 E ac 1    E ac 2  
Ac 1 exp  -  1− 1/2 ≈ A exp  -  1− 1/2
 RT   P ( T )  c2  RT   Peq ( T ) 
 eq  
A partir de cette équation, nous pouvons écrire :
 E −E 
Ac 1 ≈ Ac2 exp  - ac 2 ac 1 
 RT 
Supposons alors qu'il existe un erreur sur l'énergie d'activation apparente (Eac) de 15 kJ.mole-1
et prenons la température à 480°C.
Si l'erreur est de +15 kJ.mole-1 alors Ac1≈0,09Ac2
Si l'erreur est de –15 kJ.mole-1 alors Ac1≈11Ac2
Ceci montre bien qu'une faible erreur sur Eac entraîne une forte erreur sur Ac.
Nous connaissons à présent l'expression de la réactivité surfacique de croissance en

fonction de la température, et nous l'utiliserons dans la modélisation non isotherme de la
cinétique de la réaction.
De même, il nous est nécessaire de déterminer les variations de la fréquence

surfacique de germination en fonction de la température.
- 106 -
III.3.2.3 Mécanisme réactionnel de germination
Nous proposons le mécanisme réactionnel suivant pour décrire le processus de

germination. Il comprend quatre étapes qui se déroulent à la surface des grains (indice
« ext ») :
1− OHOHext →
← OOHext +H i
' •
2− OHOHext +H•i +s →
•
← H2O−s+VOHext
3− H2O−s →
← H2O+s
4− nO'OHext +nVOHext →
• ×
← nOOMK
Dans l'étape 4, n représente le nombre d'unités de construction du réseau cristallin de

la kaolinite nécessaire à la formation d'un germe (n≥2).
L’étape 1 correspond à la formation de défauts de structure à partir d’un hydroxyde en

surface qui se dissocie pour donner un ion oxygène et un proton en position interstitielle. La
deuxième étape correspond à la réaction entre ce proton, un autre hydroxyde et un site
d'adsorption libre pour former une molécule d’eau adsorbée, ainsi qu’une lacune en position
hydroxyde. La troisième étape est la désorption d’une molécule d'eau. Enfin, la dernière étape
est la formation du germe par réaction entre les défauts formés dans les étapes 1 et 2.
Remarque : ce mécanisme s'inspire de ceux proposés par Bouineau [14] et Brun [17] pour
les réactions de déshydratation du sulfate de lithium monohydraté, de décomposition du carbonate de
calcium et de réduction de l'octooxyde de triuranium par l'hydrogène.
Nous supposons une étape limitante pour la germination. Il peut s’agir a priori de
n’importe laquelle des quatre étapes. Pour chaque cas, nous avons calculé la vitesse de
germination en considérant les autres étapes à l’équilibre. Ces expressions sont rassemblées
dans le tableau III.3. Nous notons kgi la constante cinétique de l’étape i, Kgi la constante
d'équilibre de l’étape i. Le produit des constantes d’équilibre des quatre étapes est égal à la
pression d’équilibre de la germination notée Pgeq ( Pgeq =Kg1Kg2Kg3K1/n
g4 ).
- 107 -
Etape limitante Expression de la fréquence surfacique de germination

 PH 2 O 
Etape 1 limitante γ 1 = k g 1 1−  (III.22)
 K g1 K g2 K g3 K g41/n 
 
 K K K 1/n   PH 2 O 
Etape 2 limitante γ 2 = k g2 1  g1 g3 g4   1−  (III.23)
1/n
K g4 1/n
(K g1 K g4 +1 )  K g3 + PH O   K g1 K g2 K g3 K g4 1/n 
 2  
 K g1 K g2 K g4 1/n   PH 2 O 

Etape 3 limitante γ 3 = k g 3 
 1 − (III.24)
K K
 g1 g2 g4 K 1/n
+ 1   K K K K 1/n 
g1 g2 g3 g4

 K g1 K g2     
n n
   PH 2 O 
γ 4 = k g1 1−
   K g1 K g2 K g3 K g4 1/n  
Etape 4 limitante (III.25)
 P
 H 2O     
Tableau III.3 : Expression de la fréquence surfacique de germination selon l’étape limitante
Toutes les lois obtenues à partir du mécanisme de germination mettent en avant une
diminution de la fréquence surfacique de germination en fonction de la pression de vapeur
d’eau. Or, expérimentalement (cf. figure III.29), on observe une fréquence surfacique de
germination constante aux faibles pressions de vapeur d’eau, puis qui augmente aux plus
fortes pressions. Aucune de ces lois ne peut rendre compte de cette augmentation.
Toutefois comme le montre la figure III.39, les lois (III.22) et (III.24) peuvent
expliquer que la fréquence surfacique de germination ne varie pas en fonction de la pression
de vapeur d'eau, si on considère que la germination se déroule loin de l’équilibre
thermodynamique de ce processus, et si on se limite aux basses pressions.
11
2.5 10
11
2 10
γ (germes.m .s )
-1
-2
11
1.5 10
11
1 10
10
5 10
0
0 2 4 6 8 10 12
P (hPa)
H2O
Figure III.39 : Confrontation de la loi cinétique de germination en fonction de la pression de vapeur

d’eau avec les valeurs déterminées à l’aide du modèle de transformation (T=450°C)
- 108 -
Pour interpréter les variations de la fréquence surfacique de germination aux basses

pressions de vapeurs d’eau, nous allons donc considérer que l’étape 1 ou l'étape 3 du
mécanisme réactionnel est l’étape limitante de la germination. Il en résulte, en considérant
l'étape 1 limitante, que la fréquence surfacique de germination dépend de la température selon
une loi de type Arrhenius (à condition de négliger l'écart à l'équilibre) :
E ag
γ = Ag exp( − ) (III.26)
RT
Si l'étape 3 est limitante, on retrouve la même loi, mais il faut supposer que
g4 〉〉 1 ou K g1 K g2 K g4 〈〈 1 .
K g1 K g2 K 1/n 1/n
Comme le montre la figure III.40, cette loi est en accord avec l’évolution de la
fréquence surfacique de germination en fonction de la température.
11
1 10
10
9 10
γ (germes.m .s )
-1
10
8 10 Eag
γ=Ag exp(−
-2
)
10 RT
7 10
10
6 10
10
5 10
10
4 10
10
3 10
720 725 730 735 740 745 750 755
T (K)
Figure III.40 : Confrontation de la loi cinétique de germination en fonction de la température avec
les valeurs déterminées à l’aide du modèle de transformation (PH2O=2,5 hPa)
L’ajustement de la loi par rapport aux valeurs de la fréquence surfacique de

germination permet la détermination des valeurs des paramètres cinétiques de cette loi :
8.1014 ≤Ag ≤2,5.1016 germe.m−2.s−1
Eag =74,8 (±11) kJ.mol−1
Nous connaissons donc à présent l'expression de γ en fonction de la température.

Utilisons maintenant les lois exprimant la réactivité surfacique de croissance et la
fréquence surfacique de germination en fonction de la température dans la modélisation en
condition non isotherme.
- 109 -
III.3.3 Modélisation en condition non isotherme
La connaissance du modèle de transformation et des mécanismes de germination et de

croissance va nous permettre d'obtenir la vitesse lorsque la température varie au cours de la
réaction.
L’expression de la vitesse en conditions non isotherme et non isobare, dont le calcul
est présenté dans l’annexe E, est donnée par :
dα =φ(t) γ(τ)s exp(−s τ γ(x)dx ) −4VmK

t
dt ∫ 0 0 ∫0
r0 ln(1−
1
β(t,τ))
dτ (III.27)
0
β +(1−β)ln(1−β)= 4VmK ∫τφ(t)dt

t
avec :
r0
Nous avons déterminé les expressions de la réactivité surfacique de croissance et de la

fréquence surfacique de germination en fonction de la température : φ(T) et γ(T) (équations
(III.19) et (III.26)). De plus, pour chaque expérience, nous connaissons (grâce au
thermocouple de mesure) la température au niveau de l’échantillon au cours du temps : T(t).
Ainsi, il est possible d’en déduire les variations des deux grandeurs cinétiques en fonction du
temps : φ(t) et γ(t).
Connaissant ces fonctions, il est possible de calculer, à l'aide du système d'équations

(III.27), la vitesse de la réaction lorsque la température varie.
Nous prenons également en compte, dans ce calcul, la distribution de tailles de grains
de la figure III.13 et calculons, à chaque fois, la vitesse pour un ensemble de grains de tailles
différentes.
1 Expériences où la température varie rapidement au début de la réaction
Pour ces expériences, la modélisation a d'abord été réalisée en supposant la

température constante tout au long de la réaction. Cependant, comme nous l'avons évoqué
précédemment, le fait de considérer que la température est constante dès le début de
l'expérience, ne permet pas de bien décrire certaines courbes cinétiques expérimentales (en
particulier pour des températures de réaction élevées). Ainsi, nous avons reconsidéré ces
courbes avec le modèle appliqué en condition non isotherme.
- 110 -
Un exemple est présenté sur la figure III.41 où nous comparons la courbe de vitesse
d'une expérience réalisée à 490°C sous 2,5 hPa et la courbe théorique obtenue avec le modèle
non isotherme. La variation de la température enregistrée depuis le début de l'expérience est
également représentée sur cette figure.
0.0015 Température 500
480
d α /dt (s )
-1
0.001
460
T (°C)
440
0.0005
Expérience
Modèle 420
0 400
0 0.2 0.4 0.6 0.8 1
α
Figure III.41 : Comparaison des courbes cinétiques expérimentale et théorique, obtenues à l'aide du
modèle non isotherme, à 490°C et pour PH2O=2,5 hPa
La courbe est obtenue directement sans ajustement par rapport à la courbe

expérimentale. Nous pouvons conclure à un assez bon accord entre la vitesse théorique et la
vitesse expérimentale sur tout le domaine réactionnel, ce qui n'était pas le cas avec le modèle
isotherme. Il est donc nécessaire de prendre en compte la variation initiale de température
pour décrire la vitesse de cette expérience. Cependant, il faut quand même noter un écart entre
les deux courbes pendant la montée en température.
Nous pouvons appliquer ce modèle non isotherme à l’ensemble des expériences
réalisées à différentes températures. Ainsi, la figure III.42 présente les courbes calculées et
expérimentales de la vitesse en fonction du degré d’avancement, pour différentes
températures et avec la même pression de vapeur d'eau, 2,5 hPa.
Ainsi, quelle que soit la température de la déshydroxylation, le modèle, qui prend en

compte à présent les variations initiales de température, permet d'assez bien décrire l'évolution
de la vitesse de la réaction tout au long de la réaction (l'accord avec l'expérience est d'autant
meilleur que la température est plus faible).
- 111 -
0.0015 490°C
Expériences
M odèle
480°C
d α /dt (s )
-1
0.001
470°C
0.0005 460°C
0
0 0.2 0.4 0.6 0.8 1
α
Figure III.42 : Comparaison des courbes cinétiques expérimentales et théoriques, obtenues à l'aide
du modèle non isotherme, pour différentes températures (PH2O=2,5 hPa)
Voyons, pour terminer, l'utilisation du modèle pour interpréter les courbes cinétiques
expérimentales obtenues dans le cas d'une programmation linéaire de la température en
fonction du temps.
1 Expériences en programmation de température
Le système d'équations (III.27) nous permet de connaître l'expression de la vitesse

pour une variation d'une des contraintes physico-chimiques au cours de la réaction et donc, en
particulier, pour les expériences effectuées avec des programmations en température. Ainsi,
les figures III.43, III.44 et III.45 représentent les courbes cinétiques calculées et
expérimentales pour différentes programmations en température :
• une programmation de 1°C.min-1 suivie d'un palier à 450°C sous 2,5 hPa de vapeur d'eau
(figure III.43) ;
(figure III.44) ;
(figure III.45).
Comme on le voit sur ces trois figures, il existe à chaque fois des écarts entre les
courbes cinétiques théoriques et les courbes cinétiques expérimentales.
- 112 -
0.00028 460
0.00024 450
0.0002 T (°C)
440
d α /dt (s )
-1
T(°C)
0.00016
430
0.00012
-5 420
8 10
-5 Expérience 410
4 10
Modèle
0 400
0 0.2 0.4 0.6 0.8 1
α
Figure III.43 : Comparaison des courbes cinétiques expérimentale et théorique pour une montée
linéaire en température de 1°C.min-1 suivie d'un palier à 450°C avec PH2O=2,5 hPa
0.00035 500
T(°C)
0.0003
480
0.00025
d α /dt (s )
-1
460
T(°C)
0.0002
0.00015 440
0.0001 Expérience 420

-5 Modèle
5 10
400
0
0 0.2 0.4 0.6 0.8 1
α
Figure III. 44 : Comparaison des courbes cinétiques expérimentale et théorique pour une montée
0.0008 T(°C) 520
0.0007
500
0.0006
d α /dt (s )
480
-1
0.0005
T (°C)
0.0004 460
0.0003
Expérience 440
0.0002 Modèle
420
0.0001
0 400
0 0.2 0.4 0.6 0.8 1
α
Figure III.45 : Comparaison des courbes cinétiques expérimentale et théorique pour une montée
- 113 -
Les différences entre les vitesses calculée et expérimentale peuvent être expliquées en
considérant l'imprécision sur les valeurs de la réactivité surfacique de croissance et surtout sur
celles de la fréquence surfacique de germination. Plusieurs causes sont explicitées dans ce qui
suit : l'incertitude sur l'estimation de γ et φ, le choix de la fonction γ(T) et la méthode de
détermination de γ.
• Incertitude estimée sur γ et φ :
Comme l'ont montré les figures III.29 et III.30, une erreur a été déterminée sur les
valeurs expérimentales de la fréquence surfacique de germination et de la réactivité surfacique
de croissance. On peut alors se demander comment celle-ci peut intervenir sur les lois φ(T) et
γ(T).
Pour cela, nous avons cherché à exprimer, de nouveau, les variations de ces deux
grandeurs en fonction de la température, toujours avec les lois (III.19) et (III.26), mais en
prenant cette fois-ci pour les températures la plus basse et la plus haute, les valeurs situées aux
extrémités des barres d'erreur, comme le montre la figure III.46. Nous obtenons ainsi deux
nouvelles courbes de variation de φ et γ en fonction de la température. Ces courbes
correspondent aux variations maximale et minimale de φ et γ dans le domaine de température
considéré. Ces variations maximale et minimale, ainsi obtenues, sont notées respectivement 1
et 2 pour la fréquence surfacique de germination et a et b pour la réactivité surfacique de
croissance.
11
1 10 -7
10
3.5 10 a
9 10 -7
3 10
γ (germes.m .s )
1
-1
φ (mole.m .s )
-1
10
8 10 2.5 10
-7
b
-2
-2
10 -7
7 10 2 10
10 -7
6 10 1.5 10
10 -7
5 10 2 1 10
-8
4 10
10 5 10
10 0
3 10
720 725 730 735 740 745 750 755 720 725 730 735 740 745 750 755
T (K)
T(K)
Figure III.46 : Estimation des lois de γ(T) et φ(T) en tenant compte de l'erreur expérimentale
Nous avons alors associé dans la loi de vitesse (équations III.27), les lois de
germination et de croissance correspondant, d'une part aux variations de φ et γ les plus
- 114 -
grandes en fonction de la température (1 et a), et d'autre part à celles correspondant aux

variations de φ et γ les plus faibles (2 et b). Nous appellerons alors les vitesses obtenues 1-a et
2-b. Le tableau III.4 rassemble les valeurs des paramètres correspondant aux lois de
germination et de croissance.
1 2
-2 -1
Ag (germe.m .s ) 2.1021 5,55.1010
γ Eag (kJ.mole-1) 150 0
a b
D0 K10K 02K 06 r0 (mole.m-2.s-1) 5,8.109 6,5.106
φ
Ea +∆H10 2+∆H02 2+∆H06 2 (kJ.mole-1) 231 190
Tableau III.4 : Valeurs des paramètres des lois cinétiques exprimant γ et φ en tenant compte de la
barre d'erreur expérimentale
Les vitesses théoriques, ainsi calculées, sont comparées, d'une part aux courbes
cinétiques expérimentales, et d'autre part aux courbes calculées précédemment sur les figures
III.47, III.48, III.49.
Remarque : l'erreur estimée sur la réactivité surfacique de croissance est plus faible que celle
estimée sur la fréquence surfacique de germination. Les différentes courbes φ(T) sont assez proches.
On peut donc penser que l'erreur estimée sur γ va davantage jouer un rôle sur la forme des courbes
cinétiques.
0.0003
Expérience
0.00025 Modèle
Modèle "2-b"
d α /dt (s )
0.0002 Modèle "1-a"

-1
0.00015
0.0001
-5
5 10
0
0 0.2 0.4 0.6 0.8 1
α
Figure III.47 : Comparaison des courbes cinétiques expérimentale et théoriques en tenant compte de
l'erreur sur les valeurs de φ et γ pour une montée linéaire en température de 1°C.min-1 suivie d'un
palier à 450°C avec PH2O=2,5 hPa
- 115 -
0.00032
0.00028
0.00024
d α /dt (s )
-1
0.0002
0.00016
Expérience
0.00012
Modèle
-5
8 10 Modèle "2-b"
4 10
-5 Modèle "1-a"
0
0 0.2 0.4 0.6 0.8 1
α
0.0008
0.0007
0.0006
d α /dt (s )
-1
0.0005
0.0004
Expérience
0.0003
Modèle
0.0002 Modèle "2-b"
0.0001 Modèle "1-a"
0
0 0.2 0.4 0.6 0.8 1
α
La figure III.47, correspondant à une programmation en température de 1°C.min-1

jusqu'à 450°C, montre un bon accord entre la courbe expérimentale et la courbe de vitesse
théorique calculée avec des variations de φ et γ les plus grandes (courbe « 1-a ») en fonction
de la température.
Avec ces mêmes lois, la figure III.48 montre, pour une programmation en température
de 1°C.min-1 jusqu'à 480°C, que la vitesse théorique est en meilleur accord avec la courbe
cinétique expérimentale, bien que la vitesse maximale expérimentale soit un peu supérieure.
- 116 -
Pour l'expérience réalisée avec une programmation de température jusqu'à 500°C, bien
qu'il subsiste entre les deux courbes des écarts, il semble que le "modèle 1-a " donne encore
le meilleur résultat.
De manière générale, la différence entre les courbes théoriques issues des deux
« modèles extrêmes » montre bien que les erreurs estimées sur les valeurs de φ et surtout de γ
jouent un rôle important sur la forme des courbes cinétiques de vitesse en condition non
isotherme.
D'autres éléments vont également contribuer à une approximation sur l'obtention de la

réactivité surfacique de croissance et de la fréquence surfacique de germination.
• Terme d'écart à l'équilibre pour la loi de la fréquence surfacique de

germination :
Une deuxième source d'erreur provient du « choix » de la loi décrivant les variations
de γ en fonction de la température. En effet, nous avons vu que dans les conditions des
expériences « isothermes », la germination se déroule loin des conditions d'équilibre, et nous
supposons alors que le terme d'écart à l'équilibre (1-P/Pgeq) est égal à 1. Ces expériences ont
été réalisées pour des températures supérieures à 450°C. Lorsque la température diminue, il
semble normal de supposer que la pression d'équilibre diminue également. Ainsi, il n'est pas
évident, pour des températures inférieures à 450°C, que le terme prenant en compte l'écart à
l'équilibre soit toujours négligeable.
Ainsi, la loi de germination choisie (équation III.26) n'est peut être pas correcte à des
températures inférieures à 450°C. Nous pensons qu'elle décrit alors moins bien les variations
de γ en début d'expérience non isotherme, ce qui entraîne une erreur sur la vitesse calculée
pour ces températures.
• Incertitude non estimée sur γ et φ :
Une dernière cause des différences observées entre les vitesses expérimentales et
théoriques provient de l'incertitude sur les valeurs φ et γ du fait de la méthode d'obtention de
ces grandeurs. En effet, nous déterminons celles-ci à l’aide de courbes qui sont considérées
- 117 -
isothermes alors qu’en fait ce n’est pas le cas en début de réaction. Ceci entraîne donc une
nouvelle imprécision, en particulier sur les valeurs de la fréquence surfacique de germination.
Pour la réactivité surfacique de croissance, les valeurs sont a priori correctes puisque
nous les avons comparées à celles obtenues par la méthode des décrochements et constaté un
bon accord.
Nous pouvons donc conclure que les écarts entre modèle et expérience proviennent
essentiellement de la difficulté de connaître précisément les variations de la fréquence
surfacique de germination en fonction de la température.
Une méthode pour résoudre ces problèmes serait d'obtenir les variations de ces
grandeurs directement à partir d'une expérience non isotherme. Il est toujours nécessaire de
connaître le modèle de transformation qui nous permet d'obtenir l'expression de la fonction
d'espace au cours du temps et donc de la vitesse. En condition non isotherme, celle-ci peut
être écrite sous la forme :
dα/dt=φ(T)E(t,φ(T),γ(T))
Ne connaissant pas les expressions de φ(T) et γ(T), il pourrait être possible de
remonter à ces dernières par comparaison des courbes cinétiques théoriques et expérimentales
et à l'aide d'un programme d'optimisation qui déterminerait pour chaque température les
valeurs de φ et γ.
Ainsi, les valeurs de la fréquence surfacique de germination et de la réactivité
surfacique de croissance seront obtenues directement sans approximation (excepté celles sur
les hypothèses physiques du modèle de transformation) et sur tout le domaine de température
de la programmation.
1 Conclusions
L'objectif de cette étude sur la déshydroxylation de la kaolinite était de développer un

modèle de germination-croissance en condition non isotherme, en appliquant la méthodologie
présentée au chapitre II.
Le modèle proposé est un modèle de germination-croissance anisotrope avec un
développement interne de la nouvelle phase ; l'étape limitant la croissance est une étape de
diffusion parallèlement aux feuillets de la structure ; les grains sont assimilés à des cylindres
de faible épaisseur et de rayons différents.
- 118 -
Pour appliquer ce modèle en condition non isotherme, il était nécessaire de connaître

les lois φ(T) et γ(T).
La variation de la réactivité surfacique de croissance a été obtenue à l'aide de deux
méthodes, d'une part à l'aide de la comparaison des expériences par rapport au modèle de
transformation, et d'autre part par la méthode des décrochements. La correspondance entre ces
deux méthodes tend à montrer que les valeurs déterminées sont correctes. De plus, un
mécanisme de croissance, qui suppose que l'étape limitante est la diffusion des ions
hydroxydes, a permis de bien interpréter les variations de cette grandeur en fonction de la
température et de la pression de vapeur d'eau.
Par contre, la détermination des valeurs de la fréquence surfacique de germination ne
repose que sur la comparaison entre le modèle de transformation en condition isotherme et
des expériences qui ne sont pas isothermes au début. Ceci engendre une première erreur sur γ,
et s'ajoute à l'incertitude déjà importante sur les valeurs de cette grandeur. De plus, le
mécanisme de germination proposé ne permet pas de bien interpréter les variations de la
fréquence surfacique de germination en fonction de la pression de vapeur d'eau, et suppose
notamment que la germination se produit loin de l'équilibre, ce qui n'est peut être plus le cas
pour les basses températures (T<450°C).
Cette forte incertitude dans la détermination de la fréquence surfacique de germination
entraîne des écarts entre les courbes cinétiques expérimentales et les courbes calculées en
condition non isotherme, et met en avant la nécessité de déterminer la fonction γ(T)
directement à partir d'expériences non isothermes.
Nous avons traité un premier exemple de modélisation cinétique en condition non

isotherme avec un cas de germination-croissance.
Dans le dernier chapitre de ce manuscrit, nous allons voir le développement d'un
modèle, à la fois en conditions non isotherme et non isobare, pour une réaction qui se déroule
selon plusieurs transformations successives : la réduction de U3O8 par l’hydrogène.
- 119 -
Conclusion
L'objet de ce travail était d'étudier la cinétique des réactions solide-gaz en conditions

non isotherme et non isobare, et de mettre en place des modèles susceptibles de décrire la
vitesse d'une réaction dans ces conditions. Nous avons appliqué cette étude à deux réactions :
la déshydroxylation de la kaolinite et la réduction de l'octooxyde de triuranium en dioxyde
d'uranium par l'hydrogène.
La méthodologie mise en place pour étudier ces systèmes repose sur le développement
de deux modèles : un modèle de transformation susceptible de décrire l'évolution des zones
réactionnelles au cours de la réaction et un modèle physico-chimique qui traduit, à l'aide de
mécanismes réactionnels, l'évolution de la réactivité surfacique de croissance et de la
fréquence surfacique de germination en fonction des contraintes intensives de la réaction.
L'utilisation de ces deux modèles nous a permis d'en déduire une expression de la vitesse
théorique pour n'importe quelle modification de température ou de pression partielle pendant
la transformation.
Pour chaque réaction, les modèles développés s'appuient sur une étude expérimentale.
Ainsi, pour la déshydroxylation de la kaolinite, la forme des courbes expérimentales, obtenues
par thermogravimétrie, nous a permis de développer un modèle de germination-croissance
anisotrope à développement interne et dont l'étape limitante est une étape de diffusion.
L'observation des grains et la structure de la kaolinite nous ont fait considérer ces grains
comme des cylindres de faible épaisseur. La diffusion se déroule parallèlement aux plans des
feuillets à travers la surface latérale des cylindres. Nous avons également essayé de prendre en
compte la granularité de la poudre dans la modélisation cinétique. Le mécanisme de
croissance développé a mis en évidence que l'étape limitante est la diffusion des ions
hydroxydes au sein de la métakaolinite. Nous en avons déduit une loi relative à la réactivité
surfacique de croissance. Pour la fréquence surfacique de germination, nous avons montré
qu'il était plus difficile de déterminer une bonne estimation de cette grandeur (incertitudes
importantes des valeurs de γ, détermination avec des courbes expérimentales non isothermes
considérées isothermes, processus de germination supposé se dérouler loin de l'équilibre).
- 120 -
A partir du modèle de transformation et du modèle physico-chimique, nous avons
déterminé l'expression de la vitesse en conditions non isobare et non isotherme, et comparé
celle-ci aux vitesses expérimentales. Nous avons alors constaté des écarts entre les courbes
calculées et les expériences, et mis en évidence que ceux-ci étaient liés à l'incertitude sur la
connaissance de la fréquence surfacique de germination.
Ce résultat renforce notre précédente remarque sur le problème de la fréquence

surfacique de germination. Dans le cas où la germination est rapide, seule la réactivité
surfacique de croissance est à prendre en compte et la modélisation en conditions non
isotherme et non isobare est plus aisée.
L'étude des deux systèmes réactionnels montre que le succès de la méthodologie

utilisée repose sur la bonne connaissance de φ et γ, c'est-à-dire des grandeurs cinétiques
relatives à la germination et à la croissance. Nous avions mis en évidence que l'incertitude sur
les valeurs de γ, en particulier, est un facteur crucial. Ainsi, notre étude soulève tout l'intérêt
de traiter le problème inverse, c'est-à-dire de rechercher les variations des grandeurs
cinétiques à partir d'expériences réalisées en conditions non isotherme et non isobare. Pour
cela, il est toutefois nécessaire de connaître le modèle géométrique et donc l'expression de la
vitesse de réaction. Ce travail d'optimisation des grandeurs cinétiques nécessite également un
traitement mathématique poussé.
Enfin, ne perdons pas de vue le projet auquel appartient ce travail, la modélisation de

réacteurs industriels. Si les modèles développés interprètent bien les réactions étudiées alors
l'étape suivante de ce travail consiste à coupler ces modèles de réactions (qui fournissent les
termes sources pour n'importe quel type de variation de température ou de pression) avec les
modèles décrivant le réacteur.
- 121 -
Annexe A : Expression de la vitesse de croissance en considérant
les hypothèses de pseudo-stationnarité et de l'étape limitante
1 Vitesse d'une étape élémentaire
Une réaction élémentaire ne se déroule que dans une seule zone réactionnelle. La
vitesse (en mole.s-1), V, d'une étape élémentaire est supposée être égale à :
V =n0φE (A.1)
où n0 est le nombre de mole mise en jeu (en supposant ici des coefficients stœchiométriques
égaux à 1), φ est la réactivité surfacique de l'étape considérée et E la fonction d'espace
caractérisant la zone réactionnelle où a lieu cette étape.
1 Vitesse de croissance d'une succession d'étapes élémentaires
La croissance est le résultat d'une succession d'étapes élémentaires qui se déroulent

dans différentes zones réactionnelles.
Nous allons montrer que si une étape est supposée limitante et que si le système
réactionnel est pseudo-stationnaire, alors la vitesse de la réaction peut s'écrire comme
dans le cas d'une seule étape élémentaire sous la forme : V =n0φE .
Pour cela nous allons démontrer deux théorèmes.
Théorème 1 : si une des étapes d'un système linéaire a ses constantes des deux réactions
inverses finies, toutes les constantes de vitesse des autres étapes étant infinies, alors la
concentration des produits intermédiaires est constante et chaque espèce intermédiaire est en
équilibre avec certains réactants et certains produits de la réaction.
Afin de démontrer ce théorème, nous supposerons un système linéaire simple, et que

tous les nombres stœchiométriques sont égaux à 1 dans toutes les étapes.
Etudions le système de n+1 étapes élémentaires mettant en jeu n intermédiaires notés
Xi. Chaque étape peut correspondre, soit à une réaction chimique, soit à une diffusion dans
l'approximation de Wagner.
- 122 -
….⇔ X1
….+ X1 ⇔ X2 + ….
…………………….
….+ Xm-1 ⇔ Xm +...
…………………….
….+ Xi-1 ⇔ Xi +...
…………………….
….+ Xj-1 ⇔ Xj +...
……………………
….+ Xn-1 ⇔ Xn +...
Dans une réaction élémentaire mettant en jeu des intermédiaires réactionnels, on

appelle facteur de vitesse, noté f, le produit de la constante de vitesse par les concentrations,
élevées à la puissance des ordres partiels, des réactants et produits, qui interviennent dans les
étapes élémentaires.
Nous notons alors kk la constante de vitesse, et fk le facteur de vitesse de l'étape k, et
k'k f'k ceux de l'étape inverse. Il est alors possible d'écrire la réactivité surfacique de la
première étape sous la forme :
 f' 
φ1 = f1 − f1' [X 1 ] = f1 1− 1 [X 1 ]  (A.2)
 f1 
La loi d'action de masse permet d'écrire :
f 1' 1
= (A.3)
f 1 [ X 1 ]e
où [X1]e représente la concentration de l'intermédiaire X1 si l'étape 1 est à l'équilibre. On en
déduit l'expression de φ1 :
φ1= f1 1− 1 
[X ]
(A.4)
 [ X 1 ]e 
De même pour n'importe quelle étape k, il est possible d'écrire :
 f ' [X k ] 
φk = fk [X k −1 ] 1− k  (A.5)
 fk [X k −1 ] 
En utilisant la loi d'action de masse, on obtient :

f k'  [ X k −1 ] 
=  (A.6)
f k  [ X k ]  e
- 123 -
D'où l'expression de φk :
φk = f k [X k −1 ] 1− k
[X ] [X k −1 ]e 
X k −1 ] [X k ]e 
(A.7)
 [
Remarque : si l'étape k est une diffusion alors l'expression de la réactivité est :
φk = D ([X k −1 ] −[X k ] ) (A.8)

x
où D est le coefficient de diffusion et x la longueur de diffusion. On peut donc également écrire la
réactivité spécifique sous la forme :
φk = D [X k −1 ] 1− k 
[X ]
(A.9)
x  [ X k −1 ] 
Pour la dernière étape, la vitesse est égale à :
φn+1 = f n+1[X n ] 1− k e 

[X ]
(A.10)
 [X k ] 
Supposons que les constantes kk et k'k soient toutes infinies, sauf celles de l'étape i :
ki et k'i sont donc finies. Donc de même, fk et f'k sont infinis contrairement à fi et f'i. L'étape i
correspond alors à l'étape limitante.
Regroupons l'ensemble des étapes à constantes infinies en deux sous-ensembles. Le
sous-ensemble A inclut toutes les étapes m précédant l'étape i (m<i). Le sous-ensemble B
contient toutes les étapes j qui suivent l'étape i (j>i).
Considérons le sous-ensemble A. L'étape 1 fait partie de ce sous-ensemble. Sa vitesse

ne peut être infinie, or, sa constante de vitesse est infinie. Pour que cette vitesse soit finie, la
seule solution est d'avoir une forme indéterminée pour l'équation (A.4) ∞x0. Ceci entraîne que
l'étape 1 est à l'équilibre et que :
[ X 1 ] = [ X 1 ]e (A.11)
L'étape 2 appartient également au sous-ensemble A. Par le même raisonnement, que
pour l'étape 1, appliqué à l'équation (A.7) (avec k=2), il est possible d'obtenir :
[ X2 ]  [ X1 ] 
  =1 (A.12)
[ X 1 ]  [ X 2 ]  e
A partir des équations (A.3), (A.6), (A.11) et (A.12), on peut en déduire la

concentration en X2 :
- 124 -
[X 2 ] = f 2 f1
(A.13)
f 2' f 1'
En continuant ce raisonnement de proche en proche, on obtient pour toutes les étapes
m<i :
m
[X m ] = ∏ f k' (A.14)
k =1 fk
Et pour notamment pour l'étape m=i-1, il est possible d'écrire :
i −1
[X i −1 ] = ∏ f k' (A.15)
k =1 fk
Ainsi, tous les intermédiaires Xm, formés par une étape du sous-ensemble A, sont en
équilibre avec les réactants et les produits intervenant dans la somme des étapes qui précèdent
l'étape de sa consommation.
Considérons à présent le sous-ensemble B, et reprenons le même raisonnement mais,

cette fois-ci, en partant de la dernière étape (n+1). Pour celle-ci, à partir de l'équation (A.10),
on peut montrer comme pour l'étape 1 que :
f n' +1
[ X n ] = [ X n ]e = (A.16)
fn
Ceci conduit de proche en proche à en déduire pour toutes les étapes j+1 (j>i), la
concentration en Xj :
n +1
[X j ]= ∏ fk'
(A.17)
k = j +1 f k
Et finalement pour l'étape i+1, la concentration en Xi est déterminée par :

n +1
[X i ] = ∏ f k'
(A.18)
k =i + 1 fk
Ainsi, tous les intermédiaires Xj, formés par une étape du sous-ensemble B, sont en
équilibre avec les réactants et les produits intervenant dans la somme des étapes qui suivent
l'étape de sa formation.
Par conséquent tous les intermédiaires réactionnels ont bien une concentration
constante et chaque espèce est en équilibre avec certains réactants et certains produits
de réaction. On parle, dans ce cas, de régime cinétique de cas pur ou de processus
limitant.
- 125 -
Théorème 2 : dans l'approximation pseudo-stationnaire d'un processus limitant, la vitesse de
la réaction peut se mettre sous la forme du produit :
V =n0φi Ei (A.19)
φi est la réactivité de l'étape limitante, les autres étapes étant à l'équilibre et Ei est la
fonction d'espace liée à la zone où se déroule cette étape limitante.
D'après le théorème précédent, les réactivités de toutes les étapes ont une forme
indéterminée (∞x0), sauf celle de l'étape limitante, qui, d'après l'équation (A.5) appliquée à
k=i, s'écrit :
 f ' [ Xi ] 
φ i = f i [ X i −1 ]  1 − i  (A.20)
 f i [ X i −1 ] 
Or, d'après les expressions des concentrations des intermédiaires de réaction Xi-1 et Xi
(équations (A.15) et (A.18)), nous pouvons écrire :
 n +1
f k' 
i −1
fk 

f i'
∏ 

k = i +1 f k
φi = f i ∏ 1 − (A.21)
f k'  fi i −1
f 
k =1 
 ∏ f k' 

 k =1 k 
Ainsi, la réactivité de l'étape limitante s'écrit :
i −1
fk  n +1 ' 
f
φi = f i ∏ 1 − ∏ k 
  (A.22)
k =1 f k'  k =1 f k 
Ecrivons à présent les bilans sur les variations des quantités de chaque espèce [12] :
dn X 1
= (φ 1 E1 − φ 2 E 2 )n0 (A.23)
dt
...................
dn X i
= (φ i Ei − φ i +1 Ei +1 )n0 (A.24)
dt
................
dn X n
= (φ n E n − φ n+1 E n+1 )n0 (A.25)
dt
Tous les produits φkEk (k≠i) sont des expressions indéterminées.
La pseudo-stationnarité permet d'écrire que toutes les quantités en intermédiaires
réactionnels sont constantes :
- 126 -
dn X k
∀k , =0
dt
On en déduit des équations (A.23) à (A.25) :
φ 1 E1 = φ 2 E 2 = .... = φ i Ei = .... = φ n E n (A.26)
Ainsi, toutes les vitesses sont égales.
On peut mesurer la vitesse à l'aide d'une étape quelconque, par exemple la première :
V = n0φ 1 E1 (A.27)
Dans la plupart des cas, cette forme est indéterminée. Cependant, l'égalité des vitesses
nous permet d'écrire que cette vitesse est égale à la vitesse de l'étape i (qui est déterminée).
V = n0φ i Ei
Nous avons ainsi démontré le théorème 2.
- 127 -
Annexe B : Pseudo-stationnarité d'un système réactionnel
1 Méthode pour savoir si un système est pseudo-stationnaire
Si un système est pseudo-stationnaire, les courbes de vitesse obtenues par deux

méthodes de mesure différentes doivent être affines entre elles. C'est ce que nous allons
démontrer, en prenant par exemple le cas où on mesure en même temps la vitesse de variation
de masse par thermogravimétrie et la puissance thermique par calorimétrie.
La thermogravimétrie permet d'accéder à la vitesse de variation de masse

caractéristique des étapes d'adsorption et désorption :
d∆m =M des×vdes −M ads×vads (B.1)
dt
où Mdes et Mads sont respectivement les masses molaires des espèces désorbées et adsorbées, et
vdes et vads sont respectivement les vitesses des étapes élémentaires de désorption et
d’adsorption.
La calorimétrie mesure la puissance thermique qui est la somme des produits des
vitesses vi par la variation de l'enthalpie de l'étape i pour toutes les étapes élémentaires :
dQ
= vi∆Hi
dt ∑
(B.2)
i
Nous supposons qu'il n'y a qu'une étape élémentaire d'adsorption de gaz et qu'une
étape élémentaire de désorption gazeuse. Nous prenons alors la première étape élémentaire
comme étape d'adsorption, et la deuxième étape élémentaire comme étape de désorption :
vads=v1 (B.3)
vdes=v2 (B.4)
Si le système est pseudo-stationnaire, il est possible d'écrire que l'étape 1 est liée à
toutes les autres :
∀i, v1=λivi (B.5)
- 128 -
On obtient alors les expressions suivantes pour la vitesse de variation de masse et la
puissance thermique :
d∆m =(λ M − M ) v (B.6)
2 des ads 1
dt
dQ  
dt  ∑
=  λ i ∆ H i  v1 (B.7)
i 
Les équations (B.6) et (B.7) permettent de relier les vitesses thermogravimétrique et

calorimétrique sous la forme :
d∆m = (λ2 M des −M ads ) dQ (B.8)
dt  λ ∆H  dt
∑ i i 
i 
Ainsi, si le système est pseudo-stationnaire, les deux mesures sont

proportionnelles au cours du temps. Il existe une affinité orthogonale entre la courbe
représentant la vitesse de variation de masse et celle représentant le flux de chaleur.
Remarque : l'enthalpie de la réaction est égale à :

 
∆H = ∑λi ∆H i  (B.9)
 i 
Ainsi, le rapport d'affinité de l'équation (B.8) permet de déterminer la valeur de
l'enthalpie de la réaction qui peut être comparée aux données des tables thermodynamiques.
1 Description du couplage thermogravimétrie-calorimétrie
Ce travail a été réalisé avec un appareil Setaram DSC 111 qui réalise en parallèle les
mesures de la variation de masse et du flux de chaleur au cours du temps. Il fonctionne en
système dynamique avec une entrée pour le gaz vecteur GV, et une entrée pour le gaz
auxiliaire. Les débits des gaz sont réglés à l'aide de débitmètres massiques situés en amont de
l'appareil et permettant donc de fixer les pressions partielles des différents gaz.
Le schéma du montage est représenté sur la figure B.1.
- 129 -
Gv
dm/dt
Gaz
auxiliaire
dQ/dt
Débitmètre Réactif
massique
Fluxmètres
Thermocouple
Creuset de référence Sortie des gaz

Figure B.1 : Schéma du couplage thermogravimétrie-calorimétrie
- 130 -
- 131 -
Annexe C : Test du φΕ
Nous développons ici une méthode décrite par Soustelle et Pijolat [9] pour montrer la
possibilité d'écrire la vitesse sous la forme du produit φE. φ ne dépend que des paramètres
intensifs de la réaction : φ(T,P..). E caractérise la géométrie du système réactionnel, et plus
particulièrement la zone réactionnelle de l'étape limitante à un instant t.
Cette méthode est basée sur une variation d'une des conditions expérimentales au
cours de la réaction comme la méthode de l'isolement [2].
Elle consiste à réaliser une série d'expériences avec au départ les mêmes conditions
expérimentales, et à modifier une des variables réactionnelles pour l'amener rapidement à une
autre valeur. On peut par exemple réaliser un décrochement en température, et passer d'une
température T1 à une température T2.
Comme le montre la figure C.1, pour chaque expérience, l’instant où le décrochement
est réalisé, est différent.
α (dα/dt)ap (t1) T2
T2
(dα/dt)ap (t0)
(dα/dt)av(t1)
(dα/dt)av (t0) T1
t0 t1 Temps
Figure C.1 : Schéma de la méthode utilisée pour vérifier l'unicité de la zone réactionnelle
La vitesse est supposée pouvoir être écrite sous la forme dα/dt=φ(T,P,...).E.

On mesure alors la vitesse de la réaction avant le décrochement (dα/dt)av et la vitesse
après (dα/dt)ap. La durée du décrochement doit être suffisamment rapide afin de supposer que
la fonction d'espace ne varie pas pendant celui-ci.
Pour la première expérience, où le décrochement a eu lieu au temps t0, les expressions
des vitesses sont alors :
(dα dt ) av ( t 0 ) = φ(T1 , P,...)E( t 0 ) (C.1)
(dα dt ) ap ( t 0 ) = φ(T2 , P,...)E ( t 0 ) (C.2)
- 132 -
Le rapport des vitesses est alors égal à :
(dα dt ) ap ( t 0 ) φ(T2 , P,...)E( t 0 ) φ(T2 , P,...)
= = (C.3)
(dα dt ) av ( t 0 ) φ(T1 , P,...)E( t 0 ) φ(T1 , P,...)
En réalisant une nouvelle expérience, et en réalisant cette fois-ci le décrochement en

température au temps t1, il est de nouveau possible de connaître les expressions des vitesses :
(dα dt ) av ( t 1 ) = φ(T1 , P,...)E( t 1 ) (C.4)
(dα dt ) ap ( t 1 ) = φ(T2 , P,...)E ( t 1 ) (C.5)
Le rapport des vitesses est une nouvelle fois égal à :

(dα dt ) ap ( t 1 ) φ(T2 , P,...)E( t 1 ) φ(T2 , P,...)
= = (C.6)
(dα dt ) av ( t 1 ) φ(T1 , P,...)E( t 1 ) φ(T1 , P,...)
En conclusion, si la vitesse peut s'écrire sous la forme du produit φE et si le

système est pseudo-stationnaire, alors le rapport des vitesses doit être constant quel que
soit l'instant où le décrochement est effectué.
Ainsi, il est possible de vérifier expérimentalement que la vitesse peut s'écrire sous la
forme φE.
- 133 -
Annexe D : Comparaison modèle-expérience. Utilisation du
paramètre de modèle
1 Comparaison modèle-expérience
Une fois la vitesse calculée à l'aide d'un modèle, il est possible de comparer les courbes
cinétiques théoriques et expérimentales. Pour un cas limite, un seul paramètres, φ ou γ, est à
ajuster pour confronter ces courbes. Par contre, dans un cas complexe de germination-
croissance, la vitesse théorique de la réaction dépend de deux inconnues : la fréquence
surfacique de germination et la réactivité surfacique de croissance. Si la réaction s'est déroulée
en conditions isotherme-isobare, φ et γ sont constants. Afin de confronter le modèle aux
expériences, il faut déterminer les valeurs de ces deux grandeurs pour chaque expérience
réalisée. Pour obtenir ces valeurs, il est possible d'ajuster les courbes calculées à partir du
modèle afin de le faire correspondre le mieux possible aux courbes expérimentales.
Cependant, il n'est pas aisé d'ajuster deux paramètres à la fois.
Pour simplifier les expressions de la vitesse et du degré d’avancement, certains

[2,9,31,32,33] proposent d'effectuer un changement de variables au niveau du temps, et
d'introduire le temps adimensionné θ défini par :
Vm A φ
θ = t (D.1)
x0
où x0 représente la dimension caractéristique du grain (le rayon initial pour une sphère ou un
cylindre, la demi-épaisseur pour une plaquette), et VmA est le volume molaire de la phase
solide initiale.
Ce changement de variable fait apparaître un paramètre A qui relie les vitesses de
germination, de croissance et les caractéristiques géométriques de la poudre. A est appelé
paramètre de modèle, et est défini par :
γ x0 s0
A = (D.2)
φ Vm A
où s0 est la surface d'un grain.
- 134 -
Pour tous les modèles, il est alors possible de déterminer les expressions du degré
d'avancement en fonction du temps adimensionné, α(θ), et la vitesse réduite, dα/dθ. Celles-ci
ne dépendent que du paramètre de modèle A.
Afin de comparer le modèle à l’expérience, il est nécessaire de définir une vitesse
adimensionnée pour le modèle et l'expérience définie par :
dα
ω 0.5théorique = dθ (D.3)
 dα 
 
 dθ  α =0 ,5
dα
ω 0.5 exp érimentale = dt (D.4)
 dα 
 
 dt α =0 ,5
Si le modèle décrit bien l'expérience, sur tout le domaine réactionnel, les vitesses dα/dt
et dα/dθ sont proportionnelles puisque θ est proportionnel à t :
Vm A φ
θ= t
x0
On obtient ainsi l'égalité :
dα dα
ω 0.5théorique ( α ) = dθ = dt = ω 0.5 exp érimentale ( α ) (D.5)
 dα   dα 
   
 dθ  α =0 ,5  dt α =0 ,5
Or, la courbe théorique ne dépend que du paramètre A. Il suffit donc de faire
correspondre les courbes théoriques, ω0.5(α), par rapport aux courbes expérimentales, en
ajustant le paramètre A uniquement.
1 Détermination de la réactivité surfacique de croissance et de la fréquence

surfacique de germination
Pour un cas limite, la vitesse ne dépend que d'une seule grandeur cinétique (φ ou γ). La
comparaison du modèle aux courbes expérimentales permet d'accéder directement à cette
grandeur.
Pour un cas complexe de germination-croissance, il est plus facile en général de
travailler en vitesses réduites. Comme nous venons de le voir, la comparaison des courbes
cinétiques théoriques et expérimentales en vitesses réduites permet la détermination du
paramètre de modèle A. Une fois le paramètre A connu, la courbe α(θ) est calculée. Or,
- 135 -
l'expérience permet d'accéder à la courbe α(t). A partir de ces deux courbes, la fonction θ(t)
Vm A φ
est tracée. D'après la définition de θ (équation D.1), celle-ci est une droite de pente .
x0
Connaissant VmA et x0, la valeur de φ est déduite dans les conditions où a été réalisée
l'expérience.
Connaissant φ et A, la valeur de γ est obtenue à partir de la définition de A :
Vm A
γ = Aφ (D.6)
x0 s0
Par conséquent, l'ajustement des courbes calculées par rapport aux courbes
expérimentales permet de déterminer les valeurs de la fréquence surfacique de germination et
de la réactivité surfacique de croissance dans les conditions où ont été réalisées les
expériences.
- 136 -
- 137 -
Annexe E : Modèle de germination-croissance anisotrope à
développement interne avec une étape limitante de croissance
de diffusion et des grains cylindriques
Ce modèle a été développé pour interpréter les courbes cinétiques expérimentales de la

déshydroxylation de la kaolinite.
Nous considérerons dans un premier temps un seul grain, puis nous verrons le calcul
pour une poudre de grains de même taille.
1 Calcul du degré d’avancement et de la vitesse pour un grain
Considérons un grain de kaolinite en cours de croissance :
r0
ri
Métakaolinite
e Kaolinite
Figure E.1 : grain en cours de croissance
Notons : e l'épaisseur du cylindre

ri le rayon interne, si la surface interne : si =2πri e
r0 le rayon initial du grain, s0 la surface initiale : s0 =2πr0 e
VK le volume de la kaolinite
VK0 le volume initial de la kaolinite
VmK le volume molaire de la kaolinite
n0 la quantité initiale de kaolinite
Le degré d'avancement s'exprime en fonction du volume sous la forme :
α =VK −0VK
0
(E.1)
VK
Or, les expressions de ces volumes sont :
VK0 =πr02e (E.2)
VK =πri2e (E.3)
- 138 -
Il est donc possible, à partir des équations (E.1), (E.2) et (E.3), d'exprimer le rayon
interne et la surface interne en fonction du degré d'avancement :
ri =r0(1−α)1/ 2 (E.4)
si =s0(1−α)1/ 2 (E.5)
De plus, comme nous l'avons vu dans le chapitre I, l'expression de la vitesse de

croissance d'un grain pour une étape limitante de diffusion s'exprime par :
dα =φ(t)GD si (E.6)
dt n0
où GD est un terme sans dimension qui tient compte de la diffusion. Pour des cylindres, ce
coefficient s'écrit (cf tableau I.3) :
l0
GD = (E.7)
r0
ri ln
ri
où l0 est une longueur arbitraire. Nous choisissons, ici, celle-ci égale à r0.
En remplaçant ri à l'aide de l'équation (E.4), le coefficient GD est égal à :

r0
GD = (E.8)
r0(1−α) ln(1−α)−1/ 2
1/ 2
D'où l'expression de la vitesse, en remplaçant l'expression du coefficient GD (E.8) dans

l'équation (E.6) :
dα =φ r0 s0(1−α)1/ 2
(E.9)
dt r0(1−α)1/ 2 ln(1−α)−1/ 2 n0
Celle-ci se simplifie sous la forme :

d α = φ s0 −2 (E.10)
dt n0 ln( 1−α )
D’après la définition du volume molaire, il est possible d’écrire :
s0 VmK
=2 (E.11)
n0 r0
A partir des équations (E.10) et (E.11), la vitesse peut s'écrire sous la forme :
d α = − 4 Vm K φ 1 (E.12)
dt r0 ln( 1 −α )
- 139 -
En intégrant l’équation (E.12) entre 0 et t, on obtient la relation entre le degré
d'avancement et le temps :
α +(1−α)ln(1−α)= 4VmK ∫0φ(t)dt

t
(E.13)
r0
La transformation est terminée au temps tf pour α=1. Il est possible de calculer cette
durée tf à partir de l’équation (E.13) :
tf r0
∫0 φ(t)dt =
4VmK
(E.14)
En conditions isotherme et isobare, la réactivité surfacique de croissance est

constante et tf est égal à :
r0
tf = (E.15)
4VmKφ
1 Calcul du degré d’avancement et de la vitesse pour une poudre
Dans le cas d'une poudre, tous les grains ne commencent pas à réagir en même temps.
Entre un instant τ et un instant τ+dτ, il s'est formé un nombre de germes :
dNτ =γ(τ)SL(τ)dτ (E.16)
où SL(τ) est l'aire de la surface de l'ensemble des grains qui n'ont pas commencé à réagir au
temps τ.
• Calcul de la surface libre de germination SL(ττ)
Soit N, le nombre de grains non attaqués à un instant τ, et N0 est le nombre initial de

grains.
Dans le cas d'un modèle où la croissance est anisotrope, un seul germe recouvre
l'ensemble du grain. Ainsi, la surface libre de germination correspond à l'ensemble des grains
non transformés :
S L ( τ ) = N ( τ )s 0 (E.17)
- 140 -
Or, la variation du nombre de grains non attaqués pendant dτ s’exprime à l’aide de la
fréquence surfacique de germination par :
dN = −γs 0 N ( τ )dτ (E.18)
L’intégration de l’équation précédente entre 0 et τ en conditions isotherme-isobare

donne :
N ( τ ) = N 0 exp( −γs0τ ) (E.19)

La surface libre de germination s’écrit donc :
SL(τ)= s0 N0 exp(−γs0τ) (E.20)
En conditions non isotherme et non isobare, la fréquence surfacique de germination,

γ, n'est plus constante et va varier en fonction du temps. Ainsi, la surface libre de germination
est égale à :
τ
S L( τ )= s0 N 0 exp( − s0 ∫0 γ ( x ) dx ) (E.21)
Les germes, nés entre τ et τ+dτ,

τ, se trouvent à un degré d'avancement noté β (t,ττ) à
dβ ( t , τ )
l'instant t. Leur vitesse est égale à .
dt
La vitesse globale de la réaction à un instant t est due à l'ensemble des germes nés
entre l'instant initial de la réaction et t :
dα = 1 γ(τ)S (τ) dβ(t,τ)dτ

t
dt N0 ∫0 L (E.22)
dt
En utilisant les équations (E.21) pour exprimer la surface libre de germination, et

(E.12) pour exprimer la vitesse correspondant à un grain, l’expression de la vitesse globale
devient :
dα = γ(τ)s exp(−s τ γ(x)dx ) −4VmKφ(t)

t
dt ∫0 0 ∫0
1 dτ
0
r0 ln(1− β(t,τ)) (E.23)
Or, d’après l’équation (E.13), le degré d’avancement d’un grain peut être obtenu par :
β +(1−β)ln(1−β)= 4VmK ∫τφ(t)dt

t
(E.24)
r0
Ainsi, afin de connaître la vitesse de la réaction, il faut résoudre le système
d’équations (E. 23) et (E.24).
- 141 -
L'équation (E.23) peut être écrite sous la forme φ.E, avec la fonction d'espace égale à :
t τ −4VmK
E(t)=∫γ(τ)s0 exp(−s0 ∫0 γ(x)dx ) 1 dτ
0
r0 ln(1− β(t,τ))
Remarque pour la prise en compte de la distribution granulométrique : pour une

diffusion limitante, la réactivité surfacique de croissance est définie par :
φ = D ∆C
l0
l0 est une longueur arbitraire. Dans le calcul de la vitesse, cette valeur est égale au rayon moyen des
grains : 0,3 µm. Cependant, dans le cas où nous considérons une distribution granulométrique, la
valeur de φ n’est pas modifiée selon le rayon des grains. Donc, la fonction d’espace pour une classe de
grains de rayon r sera égale à E’(r,t).
r0
E '(r , t ) = E (r0 , t ) ×
r
En conditions isotherme-isobare, la réactivité surfacique de croissance et la

fréquence surfacique de germination sont constantes, et le système d’équations (E.23) et
(E.24) s’écrit alors :
dα =VmKφγs0 exp(−γs τ)
t
−4
r0 ∫0
dτ (E.25)
dt 0
ln(1− β(t,τ))
4VmKφ
β +(1− β)ln(1− β)= (t −τ) (E.26)
r0
• Passage en temps adimensionné
Le temps adimensionné θ est défini par :

VmKφ Vm φ
θ= t η= K τ
r0 r0
Ce passage en temps adimensionné fait apparaître un paramètre A défini par :
γ r0 s 0
A=
φ Vm K
Nous pouvons alors calculer la vitesse en temps adimensionné (en condition isotherme
et isobare) qui est égale à :
- 142 -
θ
dα =2A exp(−Aη) −2 dη
dθ ∫0 ln(1− β))
(E.27)
avec β +(1− β)ln(1− β)=4(θ −η) (E.28)
• Comparaison avec l'expérience
Une nouvelle fois, l’expression de la vitesse en temps adimensionné ne dépend que du

paramètre A.
Les courbes théoriques, ω0,5(α), sont alors ajustées par rapport aux courbes
expérimentales en faisant varier seulement le paramètre de modèle et en déterminant la valeur
qui donne le meilleur accord (cf. annexe D). Un exemple est ainsi présenté sur la figure E.2
en comparant à une expérience de déshydroxylation réalisée à 450°C sous 10-3 hPa de vapeur
d'eau avec des courbes calculées obtenues pour différentes valeurs de A.
4
A=20 Expérience (T = 450°C,P
-3
=10 hP a)
H2O
3.5
A=15 Modèle
3
2.5
0,5
A=10
ω
1.5
A=1
1
0.5
0
0 0.2 0.4 0.6 0.8 1
α
Figure E.2 : Confrontation du modèle avec l’expérience. Détermination du paramètre A
A partir de l’exemple de la figure E.2, nous voyons comment déterminer le paramètre

de modèle. A 450°C et sous 10-3 hPa de pression de vapeur d’eau, ce paramètre est égal à 15.
Comme nous l'avons montré dans l'annexe D, à partir de la droite θ(t) de pente
VmKφ
, on accède la valeur de φ :
r0
φ = 9,210-8 mole.m2.s-1
- 143 -
Connaissant A et φ, la fréquence surfacique de germination est alors égale à :
Aφ Vm K
γ= = 4,75.1010 germes.m2.s-1
r0s 0
• Estimation de l'erreur sur les grandeurs cinétiques
Afin d'estimer l'erreur sur les grandeurs cinétiques, nous cherchons d'abord les valeurs
« limites » de A. Comme le présente la figure E.3, nous encadrons la courbe expérimentale
(expérience réalisée à 450°C sous 10-3 hPa de vapeur d'eau) à l'aide de courbes calculées, pour
lesquelles nous considérons que le modèle est encore satisfaisant. Ainsi, nous déterminons
deux valeurs de A « limites » :
A=15±2
3.5 A=17
3 Expérience
2.5 Modèle
A=13
2
0.5
ω
1.5
0.5
0
0 0.2 0.4 0.6 0.8 1
α
Figure E.3 : Confrontation du modèle avec l’expérience. Détermination de l'erreur sur A
Pour chacune des courbes calculées, nous pouvons calculer une valeur de φ, et nous en
déduisons une erreur sur cette grandeur :
φ = 9,2.10-8 (±0,6.10-8) mole.m2.s-1
L'erreur sur la fréquence surfacique de germination est obtenue à partir de l'erreur sur
les valeurs de A, et de celle sur les valeurs de la réactivité surfacique de croissance. On
obtient alors :
γ = 4,75.1010 (±0,5.109) germes.m2.s-1
- 144 -
Références bibliographiques
[1] A.K.Galwey, M.E.Brown, "Thermal Decomposition of ionics solids", Elsevier,

Amsterdam, 1999
[2] B.Delmon, "Introduction à la cinétique hétérogène", Technip, Paris, 1969
[3] P.Barret, "Cinétique hétérogène", Gauthier Villars, Paris, 1973
[4] M.E.Brown, D.Dollimore, A.K.Galwey, "Reactions in the solid state", Comprehensive

Chemical Kinetics, Elsevier, Amsterdam, Vol.22, 1980
[5] P.W.M Jacobs, F.C.Tomkins, "Chemistry of the solid state", Butterworth, London, ,
Chap7, 1955,
[6] F.C.Tomkins, "Treatise on solid state chemistry", Hannay, New York, Vol 4, Chap 4,
1976
[7] J.Szekely, J.W.Evans, H.Y.Sohn, "Gas-solid reactions", Academic press, New York, 1976
[8] Avrami, J.Chem.Phys., 7 (1939) 1103, 8 (1940) 212, 9 (1941) 177
[9] M.Soustelle, M.Pijolat, Solide State Ionics, 95 (1997) 33
[10] P.Sarrazin, A.Galerie, J.Fouletier, "Les mécanismes de la corrosion sèche : Une approche
cinétique", EDP Sciences, Paris, 2000, p.68
[11] K.Surla, F.Valdivieso, M.Pijolat, M.Soustelle, M.Prin, Solid State ionics, 143 (2001) 355
[12] M.Soustelle, "Modélisation macroscopique des transformations physico-chimiques",

Masson, Paris, 1990
[13] R.Lalauze, A.Souchon, M.Soustelle, Oxydation of metals, 10(2) (1976) 105
[14] V.Bouineau, Thèse, "Expérimentation et modélisation des réactions de décomposition

isotherme-isobare des solides : application au sulfate de lithium monohydrate et au carbonate
de calcium", Saint Etienne, 1998
[15] M.Tupin, Thèse, "Approfondissement des mécanismes d'oxydation de deux alliages de

zirconium : zircaloy-4 et ZrNbO, sous oxygène et sous vapeur d'eau. Comparaison des
régimes cinétiquement limitant", Saint Etienne, 2002
[16] J.P.Viricelle, M.Pijolat, M.Soustelle, J.Chem. Soc. Faraday Trans., 91(24) (1995) 4437
[17] C.Brun, Thèse, "Modélisation de la cinétique de réduction d'oxydes d'uranium par

l'hydrogène, réactivités de germination et de croissance", Saint Etienne, 1997
- 145 -
[18] J.P.Viricelle, Thèse, "Etude de la transformation thermique de l'hydroxycarbonate de
cérium III en dioxyde de cérium IV", Saint Etienne, 1994
[19] J.H. Sharp, G.W. Brindley and B.N.N. Achar, J. Amer. Ceram. Soc., 49 (1966) 379
[20] W.A Johnson, R.F.Mehl, Trans.AIME, 135 (1939) 416
[21] B.V.Erofeev, C.R. Dokl. Akad. Sci. URSS, 52 (1946) 511
[22] J.B.Holt, J.B.Cutler, M.E.Wadsworth, J. Amer. Ceram. Soc., 45 (1962) 133
[23] A.M.Ginstling, B.I.Brounshtein, Zh. Prikl. Khim, 23 (1950) 1327
[24] G.Valensi, C.R Acad. Sci., C, 202 (1936) 309
[25] R.E.Carter, J. Chem. Phys, 34 (1961) 2010, 35 (1961) 1137
[26] W.Z.Jander, Z. Anorg. Allg. Chem., 163 (1927) 1
[27] H.Schmalzried, "Treatise on solid state chemistry", Hannay, New York, Vol 4, Chap 5,
1976
[28] H.G.McIlvried, F.E.Massoth, Ind. Eng. Chem. Fundam 12 [2] (1973) 225
[29] K.L.Mampel, Z. Phys. Chem., A 187 (1940) 43 et 235
[30] E.G.Prout et F.C.Tomkins, Trans. Faraday. Soc., 40 (1944) 488
[31] F.Valdivieso, V.Bouineau, M.Pijolat, M.Soustelle, Solid State ionics, 101 (1997) 1299
[32] C.Brun, F.Valdivieso, M.Pijolat, M.Soustelle, Phys.Chem.Chem.Phys., 1 (1999) 471
[33] V.Bouineau, M.Pijolat, M. Soustelle, Journal of the European Ceramic Society, 18

(1998) 1319
[34] J.P.Viricelle, M.Pijolat, M.Soustelle, J.Chem. Soc. Faraday Trans., 91(24) (1995) 4431
[35] J.P Dvoraczek, D.Delafosse, P.Barret, J.C.Colson, B.Sorbo, Bulletin de la société

chimique de France, 10 (1970) 3421
[36] J.H.Flynn, Journal of Thermal Analysis, 34 (1988) 367
[37] G.Bertrand, M.Lallemant, G.Watelle, Journal of Thermal Analysis, 13 (1978) 525
[38] J.Rouquérol, Journal of Thermal Analysis, 5 (1973) 203
[39] P.D.Garn, Journal of Thermal Analysis, 13 (1978) 581
[40] B.Topley, M.L.Smith, J.Chem.Soc., 1935 321
- 146 -
[41] K.Miyokawa, I.Masuda, Thermochimica Acta, 86 (1985) 113
[42] A.K.Lahiri, Thermochimica Acta, 40 (1980) 289
[43] N.Koga, J.M.Criado, Journal of Thermal Analysis, 49 (1997) 1477
[44] N.Koga, J.M.Criado, J. Amer. Ceram. Soc, 81 [11] (1998) 2901
[45 ] H.B.Johnson, F.Kessler, J. Amer. Ceram. Soc, 52 [4] (1969) 199
[46] R.W.Hutchinson, S Kleinberg, F.P.Stein, J. Phys. Chem, 77[7] (1973) 870
[47] F.Carrasco, Thermochica Acta, 213 (1993) 115
[48] W.Gomes, Nature, 192 (1961) 865
[49] W.B.Hillig, "Kinetics of high temperature processes", W.D Kingery, 1950, p.311
[50] A.W.Coats, J.P.Redfern, Analyst 88 (1963) 906
[51] M.J.Tiernan, P.A.Barnes, G.M.B.Parkes, J. Phys. Chem., 103 (1999) 338
[52] H.L.Friedman, Journal of Polymer Science, 50 (1965) 183
[53] H.E.Kissinger, Analytical Chemistry, 29 (1957) 1702
[54] E.S.Freeman, B.Caroll, J. Phys. Chem., 73 (1969) 751
[55] T.Ozawa, Bulletin of the Chemical Society Japan, 38 (1965) 1881
[56] A.W.Coats, J.P.Redfern, Nature, 201 (1964) 68
[57] Y.Pan, X.Guan, Z.Feng, Y.Wu, X.Li, Journal of Thermal Analysis and Calorimetry, 55
(1999) 877
[58] K.V.Tomashevitch et al, Thermochimica Acta, 323 (1998) 101
[59] C.R.Li, T.B.Tang, Thermochimica Acta, 325 (1999) 43
[60] Koga, H.Tanaka, Thermochimica Acta, 183 (1991) 125
[61] B.N.N.Achar, G.W.Brindley, J.H. Sharp, Proceedings of the International Clay

Conference (Jerusalem), 1 (1966) 67
[62] J.M.Criado, A.Ortega, F.Gotor, Thermochimica Acta, 157 (1990) 171
[63] F.Rouquérol, J.Rouquérol, Journal of Thermal Analysis, 1 (1971) 373
[64] J.H.Flynn, B.Dickens, Thermochimica Acta, 15 (1976) 1
- 147 -
[65] A.K.Galwey, M.E.Brown, Journal of Thermal Analysis and Calorimetry, 60 (2000) 863
[66] J.Baram, V.Erukhimovitch, Thermochimica Acta, 291 (1997) 81, 323 (1998) 43
[67] E.V.Boldyreva, Thermochimica Acta, 110 (1987) 107
[68] M. Maciejewski, Journal of Thermal Analysis, 33 (1988)1269
[69] F.A.Kröger, "The Chemistry of Imperfect Crystals", North Holland Publ.Comp,

Amsterdam, 1964, p 796
[70] P.Dion, Thèse, "Déshydroxylation de la kaolinite par analyse thermique à vitesse de

transformation sontrôlée. Etude de la métakaolinite", Orléans, 1994
[71] O.Knacke, O.Kubaschewski, H.Kesselmann, “Thermochemical Properties of Inorganic

Substances”, 2nd edition, Springer-Verlag, Berlin, 1991
[72] F.Toussaint, J.J.Fripiat, M.C.Gastuche, J.Phys.Chem., 67 (1963) 26
[73] N.H.Brett, K.J.D.Mackenzie, J.H.Sharp, Rev.Chem.Soc., 24 (1970) 185
[74] R.L.Frost, A.M.Vassallo, Clays and Clay minerals, 44 (1996) 635
[75] H.F.WTaylor, Clay Miner.Bull., 5 (1962) 44
[76] M.Gabor, J.Wajand, L.Popal, Z.G.Szabo, React.Solids, Proc. Int. Symp., 8th, 1977
[77] G.W.Brindley, M.Nakahira, J.Amer.Ceram.Soc., 40 (1957) 346
[78] P.R.Suitch, J.Amer.Ceram.Soc., 69 (1986) 61
[79] P.Murray, J.White, Trans.Brit.Ceram.Soc., 48 (1949) 187
[80] F.Vaughan, Clay Miner. Bull., 2 (1955) 265
[81] J.B.Holt, I.B.Cutler, M.E.Wadsworth, J.Amer.Ceram. Soc., 45 (1962) 133
[82] G.W.Brindley, J.H.Sharp, J.H.Paterson, B.N.N.Achar, The American Mineralogist, 52

(1967) 201
[83] I.Horvath, Thermochimica Acta, 85 (1985) 193
[84] S.A.T.Redfern, Clay Miner., 22 (1987) 447
[85] J.B.Holt, J.B.Culter, M.E.Wadsworth, J.Amer.Ceram.Soc.,45 (1962) 133
[86] C.Otero, M.Letellier, B.C.Gerstein, J.Fripiat, Dev. Sedimentol., 35 (1982) 73
[87] A.Ortega, F.Rouquérol, S.Akhouayri, Y.Laureiro, J.Rouquérol, Appl.Clay Sci., 8 (1993)

207
- 148 -
[88] K.Nahdi, P.Llewellyn, F.Rouquérol, J.Rouquérol, N.K.Ariguib, M.T.Ayedi,
Termochimica Acta 390 (2002) 123-132
[89] E.Garcia-Diaz, Thèse, "Réactivité pouzzolanique des métakaolinites corrélations avec les
caractéristiques minéralo gitologiques des kaolinites", Alès, 1995
[90] K.J.D.Mackenzie, I.W.M.Brown, R.H. Meinhold, M.E. Bowden, J.Amer.Ceram.Soc., 68

(1985) 293
[91] Handbook of Chemistry and Physic, 79th edition, 1998-1999, CRC Press, p4-141
[92] Kofstad, "Non stoichiometry, diffusion and electrical conductivity in binary metal
oxides", Wiley interscience, New York, 1972, p300-318.
[93] B.Belboech, "Nouveau traité de chimie minérale" Tome XV quatrième fascicule ChapII,
Masson et Cie, Paris, 1967, p558
[94] F.Grønvold, Nature, 162 (1948) 70
[95] A.J.Andresen, Symp.Reactor Materials, Stockholm, 1959
[96] B.O.Loopstra, Act.Cryst., 17 (1964) 651
[97] P.Pério "Nouveau traité de chimie minérale" Tome XV deuxième fascicule ChapVI,
Masson et Cie, Paris, 1961, p264
[98] L.M.Kovba, N.I.Komarevtseva, E.U.Kuz'micheva, Radiokhimiya, 21(5) (1979) 754
[99] B.T.Willis, Nature, 197 (1963) 755
[100] B.Belboech, C.Piekaski, P.Pério, Acta Cryst., 14 (1961) 837
[101] R.E.De Marco, M.G.Mendel, J.Inorg.Nucl.Chem., 14 (1953) 132
[102] S.A.Morrow, S.Graves, L.Tomlinson, Trans.Far.Soc., 57 (1960) 1400
[103] R.M.Dell, V.J.Wheeler, Trans.Far.Soc., 58 (1962) 485
[104] A.H.Le Page, A.G.Fane, J.Inorg.Nucl.Chem., 36 (1974) 87
[105] K.J.Notz, M.G.Mendel, J.Inorg.Nucl.Chem., 14 (1950) 55
[106] Y.Balaji Rao, R.B.Yadav, R.Narazana Swany, B.Gopalan, S.Syamsundar, Journal of

thermal analysis, 44 (1995) 1439
[107] T.Challier, Thèse, "Elaboration et caractérisation des formes cristallines du trioxyde

d'uranium. Réactions d'intercalation et d'insertion dans les composés d'uranium hexavalent",
Montpellier, 1991
- 149 -
[108] H.W.G.Heynen, C.G.M.M.Camp-Van Berkel, H.S. Van Der Baan, Journal of catalysis,
48 (1977) 386
[109] T.Gilardi, Thèse, "Etude par analyse thermique à vitesse de transformation contrôlée
des mécanismes d'oxydation et de réduction des oxydes d'uranium", Marseille, 1993
[110] D.Chaillot, Thèse, "Etude du mécanisme de réduction du difluorure d'uranyle par le

dihydrogène en présence de vapeur d'eau", Marseille, 2000
[111] S.Aronson, J.C.Clayton, J.Inorg.Nucl.Chem., 7 (1958) 384
[112] A.Genin, Thèse, Nancy, à paraître
[113] J.Crank, “Mathematics of Diffusion”, Oxford University Press, London, 1956
- 150 -
ECOLE NATIONALE SUPERIEURE N° d’ordre : 296CD
DES MINES DE SAINT-ETIENNE
Nom : Stéphane PERRIN
Titre : Modélisation de la cinétique de transformations non isothermes et (ou) non

isobares. Application à la déshydroxylation de la kaolinite et à la réduction de
l'octooxyde de triuranium par l'hydrogène
Spécialité : Génie des Procédés
Mots clés : Modélisation

Cinétique
Germination
Croissance
Déshydroxylation
Kaolinite
Réduction
Oxyde d’uranium
Résumé :
L’objectif de ce travail est de pouvoir décrire des transformations mettant en jeu des
solides et des gaz en conditions non isotherme et non isobare, à l’aide de modèles cinétiques.
Afin d’atteindre cet objectif, nous avons mis en place une méthodologie. Deux
processus essentiels doivent être pris en compte : la germination et la croissance. Les germes
sont supposés se former (à température et pression constantes) en surface des grains avec une
vitesse constante par unité de surface, γ, appelé fréquence surfacique de germination (nombre
de germes.m-2.s-1). La vitesse de croissance est caractérisée par une réactivité surfacique de
croissance, φ (en mol.m-2.s-1). Avec un modèle de transformation approprié, il est possible
d’obtenir les variations de γ et φ en fonction de la température et de la pression qui sont
ensuite utilisées dans le calcul de la vitesse en conditions non isotherme et non isobare.
Afin de valider la méthode développée, deux réactions ont été étudiées. Pour la
première, la déshydroxylation de la kaolinite, un modèle de germination-croissance anisotrope
où l’étape limitant la croissance est une étape de diffusion, a été développé afin d’interpréter
les courbes cinétiques expérimentales. Cependant les courbes de vitesse calculées à partir de
ce modèle ne permettent pas de décrire la réaction pour certaines variations de température.
Ce résultat met en avant la difficulté de déterminer précisément la fréquence surfacique de
germination ce qui engendre une approximation importante sur les courbes cinétiques.
La deuxième réaction est la réduction de l’octooxyde de triuranium par l’hydrogène.
Nous avons montré que celle-ci se déroulait selon trois transformations successives. Nous
avons développé un modèle cinétique pour chacune de ces réactions en considérant cette fois-
ci la germination comme instantanée. Enfin en comparant ce modèle aux courbes de vitesse
expérimentales, nous avons constaté un très bon accord aussi bien pour une variation de
température que pour un changement de la pression partielle d’hydrogène au cours de la
réaction.

TheseSPerrin Vep

Transféré par

Droits d'auteur :

Formats disponibles

TheseSPerrin Vep

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

TheseSPerrin Vep

Transféré par

Droits d'auteur :

Formats disponibles

ECOLE NATIONALE SUPERIEURE N° d’ordre : 296CD

DES MINES DE SAINT-ETIENNE

pour obtenir le grade de

(Spécialité : GENIE DES PROCEDES)

Modélisation de la cinétique de transformations non isothermes et

Soutenue à Saint-Etienne le 19 Décembre 2002

Monsieur D. ABLITZER Président

Madame F. ROUQUEROL Rapporteur

Messieurs J. M. CRIADO Examinateurs

Mesdames C. ABLITZER-THOUROUDE Invités

Thèse préparée au Centre « Sciences des Processus Industriels et Naturels »

Je tiens à remercier messieurs René LALAUZE et Michel COURNIL pour m’avoir

Je remercie également Mesdames Françoise ROUQUEROL, Carine ABLITZER-

Mes remerciements s’adressent également à l’équipe des mathématiciens du centre

Chapitre I : Etude bibliographique 5

I.1 Modèles de transformation et lois cinétiques 6

I.1.1 Lois cinétiques de germination 7

I.1.2 Lois cinétiques de croissance 11

I.1.3 Modèles de transformation 19

I.1.4 Influence de la pression et de la température 28

I.1.4.1 Approche utilisant la séparation des variables 29

I.1.5 Effet d'une distribution granulométrique 39

I.2 Modélisation en condition non isotherme 41

I.2.2 Classification des méthodes non isothermes 42

I.2.2.1 Programmation en température 43

Chapitre II : Approche non isotherme et non isobare de

II.1 Intérêts de la modélisation 54

II.2 Principe de la modélisation en conditions non isotherme

II.2.1 Cas des réactions avec germination-croissance 58

II.2.2 Cas des réactions successives 62

II.3 Choix des systèmes réactionnels 63

Chapitre III : Modélisation d’une réaction avec germination-

III.1 Etude bibliographique 66

III.1.1 Structure de la kaolinite 66

III.1.2 Etude thermodynamique de la réaction 67

III.1.3 Cinétique de la déshydroxylation de la kaolinite 69

III.2 Etude expérimentale 71

III.2.1 Caractérisation des échantillons 71

III.2.1.1 Diffractogrammes de la kaolinite et de la métakaolinite 71

III.2.2 Etude cinétique de la réaction par thermogravimétrie 79

III.2.2.1 Protocole expérimental 79

III.3 Modélisation de la réaction 83

III.3.1 Modèle de transformation 83

III.3.1.1 Vérification des hypothèses du modèle 83

III.3.2. Modélisation physico-chimique 95

III.3.2.1 Variation de φ et γ en fonction des conditions physico-chimiques 95

III.3.3 Modélisation en condition non isotherme 110

Chapitre IV : Modélisation de réactions successives :

IV.1 Etude bibliographique 122

IV.1.1 Diagramme thermodynamique 122

IV.1.2 Structure des oxydes d'uranium étudiés 123

IV.1.2.1 L'octooxyde triuranium : U3O8 123

IV.1.3 Etude cinétique et modélisation de la réaction 126

IV.2 Etude expérimentale 134

IV.2.1 Etude de la réaction par thermogravimétrie 134

IV.2.1.1 Protocole expérimental 134

IV.2.2 Caractérisation texturale des produits 142

IV.2.2.1 Caractérisation par adsorption d'azote à 77 K 142

IV.2.3.1 Diffractogrammes de U3O8 et de UO2 148