N˚ d'ordre

N˚d’ordre : 2453
THESE
Pr´esent´ee
pour obtenir le titre de
DOCTEUR DE L’INSTITUT NATIONAL POLYTECHNIQUE DE
TOULOUSE
Ecole doctorale : Transferts, Dynamiques des Fluides, Syst`emes & Proc´ed´es
Sp´ecialit´e:Energ´etique et Transferts
Par
Idriss AHMED HARED
Optimisation d’un proc´ed´e de pyrolyse en four tournant :
application `a la production de charbons actifs
Soutenue le 14 f´evrier 2007
a` l’Ecole des Mines d’Albi-Carmaux devant le jury compos´ede:
Pr. FLAMANT Gilles PROMES, Perpignan Pr´esident, rapporteur
Dr. MARIAS Fr´ed´eric ENSGTI, Pau Rapporteur
Dr. VAN DE STEENE Laurent CIRAD, Montpellier Examinateur
Pr. LACROIX Marcel UNIVERSITE DE SHERBROOKE Examinateur
Pr. SALVADOR Sylvain EMAC, Albi Examinateur
Dr. DIRION Jean-Louis EMAC, Albi ExaminateurRemerciements
Je profite de cette occasion pour remercier toutes les personnes qui m’ont aid´e`ar´ealiser mes travaux
de th`ese.
Je tiens `a remercier M. Didier LECOMTE pour m’avoir accueilli au sein de l’ancien Centre Energ´etique
et Environnement de l’Ecole des mines d’Albi-Carmaux devenu depuis Centre RAPSODEE. Je tiens a`
remercier l’actuel directeur du Centre RAPSODEE M. Jacques FAGES.
J’exprime mes vifs remerciements `a mon directeur de th`ese M. Marcel LACROIX pour ses pr´ecieux
conseils et ses encouragements.
Je tiens `a remercier particuli`erement mon co-directeur de th`ese M. Sylvain SALVADOR pour sa tr`es
grande disponibilit´e, son aide pr´ecieuse et ses conseils. Merci de m’avoir communiquer le goutˆ pour la
recherche.
Je voudrais remercier ´egalement M. Jean-Louis DIRION, mon encadrant de th`ese pour son aide et ses
conseils et notamment en mati`ere de calcul num´erique.
Je remercie M. Gilles FLAMANT et M. Fr´ed´eric MARIAS d’avoir accept´ed’ˆetre rapporteurs de mes tra-
vaux de th`ese. Merci pour leurs critiques constructives et leurs analyses pertinentes. Je remercie ´egalement
M. Laurent VAN DE STEENE d’avoir ´et´emembredecejuryetdem’avoirapport´e des suggestions pour
am´eliorer le contenu du manuscrit.
Je voudrais associer `a ces remerciements, les techniciens du Centre RAPSODEE et en particulier Jean-
Claude POUSSIN pour son aide technique pr´ecieuse et sa sympathie. Merci ´egalement `a Ludivine MOGA,
Denis MARTY, Jean-Marie SABATHIER, Bernard AUDUC et Sylvie DELCONFETTO pour leur aide et
leur disponibilit´e.
Je n’oublierai pas de remercier tous mes coll`egues doctorants et post-doctorants et en particulier mes
coll`egues de bureau Cedric, Marcio, Hicham et Damien.
Je remercie mon ´epouse Nadira qui m’a beaucoup soutenu pendant ces travaux de th`ese. Merci `atoute
ma famille.
Merci `a tous ceux qui, de pr`es ou de loin ont contribu´e`a ce travail.A Nadira et Deksane
A ma familleTable des mati`eres
Introduction 1
1Pr´esentation du proc´ed´e5
1.1 Le four tournant : un r´eacteurpolyvalent..................... 5
1.2 G´en´eralit´essurlesfourstournants.... 6
1.2.1 Les diff´erentstypesdetuberotatif .......... 7
1.2.2 Les m´ethodes d’apport de l’´energie.......... 8
1.3 Pr´esentation du four tournant ´etudi´e............. 10
1.3.1 Caract´eristiquestechniquesdufourtournant............... 10
1.3.2 Chauffagedufour......... 12
2 Bibliographie 15
2.1 Mouvementdelachargedansunfourtournant ................. 16
2.1.1 Lestypesdemouvementenfourtournant...... 16
2.1.2 Le m´elange transversal et la s´egr´egation ....... 18
2.1.3 Letransportaxial.............................. 21
2.2 Lestransfertsdechaleurdanslesfourstournants...... 26
2.2.1 M´ecanismesdetransfertdechaleur.......... 27
2.2.2 Les ph´enom`enesdeconvectiondanslesfourstournants......... 27
2.2.3 La conduction paroi-solide ......................... 29
2.2.4 Les ´echangesradiatifs....... 31
2.3 Etatdel’artsurlapyrolysedubois... 32
2.3.1 Etudequalitativedelapyrolysedebois.................. 34
i2.3.2 Mod`eles cin´etiquesdelapyrolysedebiomase.............. 39
2.4 Etat de l’art sur l’´etudedescharbonsactifs.......... 46
2.4.1 Influence des param`etres op´eratoires de synth`ese... 48
2.4.2 Les m´ecanismesd’actiondel’acidephosphorique............. 53
2.5 Conclusion ........................... 58
3 Etude de la cin´etique de pyrolyse du bois impr´egn´e par de l’acide phospho-
rique 61
3.1 Caract´erisation du bois utilis´e........................... 61
3.2 Exp´eriences en thermogravim´etrie.... 62
3.3 Mod´elisationdelapyrolyseduBIAP.. 67
3.3.1 Choix d’un sch´ema r´eactionnel....................... 67
3.3.2 Proc´edure d’identification param´etrique........ 70
3.4 Conclusion................. 73
4 Etude du d´eveloppement de la porosit´e des charbons actifs 75
4.1 Description du r´eacteur `a lit travers´e....................... 76
4.2 Protocole exp´erimental.......... 78
4.3 D´eveloppement de la porosit´e...... 80
4.3.1 Influence de la temp´erature ........................ 81
4.3.2 Influence du temps ........ 84
4.4 Cin´etique de d´eveloppement des surfaces sp´ecifiques............... 86
4.4.1 D´etermination d’un sch´ema r´eactionnel........ 86
4.4.2 Evaluation du sch´ema r´eactionnel..................... 8
4.4.3 Limites du mod`ele........ 90
4.5 Conclusion ................ 91
5 Elaboration d’un mod`ele global de four tournant de production de charbon
actif 93
5.1 Construction du mod`eledefourtournant..................... 95
5.1.1 Ecoulementdelachargedanslefourtournant.... 95
ii5.1.2 Equations de conservation de la mati`ere ................. 96
5.1.3 Mod´elisationthermique ................ 99
5.2 R´esolutiondusyst`eme..........10
5.3 R´esultats............................13
5.4 Conclusion.................121
6 Exploitation du mod`ele de four tournant et recherche de profils optimaux
de temp´erature 123
6.1 Influences des conditions op´eratoires .......................123
6.1.1 Influence des temp´eraturesdeszones.........124
6.1.2 Influence du d´ebitd’alimentationensolide......126
6.1.3 Influence de la vitesse de rotation et de l’inclinaison . ..........128
6.2 Limites de la m´ethode de r´esolution..............131
6.3 Lin´earisation de la temp´eraturedeparoi132
6.4 Optimisation du profil de temp´eraturedelaparoi................134
6.4.1 G´en´eralit´essurl’optimisation.............134
6.4.2 Rechercheduprofildeparoioptimal...................136
6.5 Conclusion.................142
Conclusion 144
Annexe A La m´ethode SQP 160
Annexe B Analyse ATG-DSC 165
Annexe C M´ethodes de r´esolution des ´equations diff´erentielles 167
Annexe D M´ethode de mesure de la porosit´e des charbons actifs 169
iiiNomenclature
A Surface d’´echange convective gaz-paroi par m (m)cgp
A Surface d’´echange convective gaz-charge par m (m)cgs
2A Surface occup´ee par le gaz dans une section du four tournant (m )g
A Surface d’´echange radiative gaz-paroi par m (m)rgp
A Surface d’´echange radiative gaz-charge par m (m)rgs
A Surface d’´echange radiative paroi-charge par m (m)rps
2A Surface occup´ee par la charge dans une section du four tournant (m )s
Cp Chaleur sp´ecifique de la charge (J/(kg.K))s
G D´ ebit massique de l’esp`ece d’indice i dans le four tournant (kg/s)i
h Hauteur de la charge (m)
2h Coefficient de convection gaz-paroi W/(m .K)cgp
2h Coefficient de convection gaz-charge (W/(m .K))cgs
2h Coefficient d’´echange de chaleur par conduction entre la paroi et la charge (W/(m .K))cps
2h Inverse de la r´esistance de contact entre la paroi et la charge (W/(m .K))c
k Conductivit´e thermique du gaz (W/(m.K))g
m Masse de bois (kg)Bois
n Vitesse de rotation (rad/s)
iv3Q D´ ebit volumique de la charge (m /s)s
T Temp´erature du gaz (K)g
T Temp´erature de la paroi du four tournant (K)p
T Temp´erature de la charge (K)s
V Vitesse axiale de la charge dans le four tournant (m/s)S
X Taux de remplissage volumique de la charge (%)r
y Masse normalis´ee de l’esp`ece d’indice ii
BIAP Bois impr´egn´e par l’acide phosphorique
BP Bois phosphat´e
DDiam`etre du four tournant (m)
k Distance transversale parcourue par une particule en une cascade (m)
L Longueur du four tournant (m)
R Rayon du four tournant (m)
2S Surface sp´ecifique (m /g)
s Distance axiale parcourue par une particule en une cascade (m)
X Taux d’impr´egnation du bois par l’acide phosphorique (%)
z Distance depuis l’entr´ee du four tournant (m)
α Angledechargecorrespondanta`un rayon r dans la couche passive (rad)
α Coefficient d’absorption du gazg
β Angle de repos dynamique de la charge (rad)
m ˙ D´ ebit massique de l’esp`ece d’indice i transf´er´e`a l’interface charge-gaz par m (kg/(m.s)i
˙R Taux de cr´eation ou de consommation de l’esp`ece d’indice i par m (kg/(m.s)i
Emissivit´e de la paroip
v Emissivit´edelacharges
λ Conductivit´e thermique de la charge (W/(m.K)