Pfe GM 0498

UNNERSIT CHEIKH ANTA DIOP DE DAKAR

1 GM 1

COLE SUPERIEURE POLYTECHNIQUE

,
Centre de This

Dpartement Gnie Mcanique

1 PROJET DE FIN D'ETUDES i

en vue de l'obtention du diplme d'ingnieur de conception

SUJET:

! Etude de faisabilit d'une unit de 1

1 production d'hydrogne par lectrolyse de 1

1 l'eau fonctionnant grce une centrale 1
1 lectro-solaire thermodynamique. i

1

Auteurs: M. Mbaye SECKet M. Moussa DIOP

1 Directeur: M. Banda NDOYE 1

Juillet 2007

L-'/I/I/I/I/171/1/I/I/I/I/I/AIPI/I/I/I/I/I/I/.... /I/...... /I/I/.......w/1/I/I/...../I/Al'YI/I/#'I7.... /...../...../I/...../...../I/I/I/I/I/.... /.... / ...../.... /I'/I/..... /..... /..... /.... /..... /I/I/I/I'/...../...../..,../I/...../...../ .... /I/I7...../I/...../I/...... /I/..... ...../I/lj
P.F.E.Prsent par M. MbayeSECKet M. MoussaDIOP
REMERCIEMENTS
E.S.P. This (2006- 2007)
Apres avoir rendu grce Dieu, pour nous avoir donn le courage et la grce d'accomplir ce
travail qui est certes le fruit de notre travail mais aussi de l'appui et des conseils de quelques
personnes, nous remercions :
> Nos parents, pour leurs sacrifices et inlassables efforts qu'ils ne cessent de dployer
afin que nous russissions et que nous nous panouissons;
> Monsieur Banda NDOYE, Directeur du projet (enseignant l'ESP), pour son
encadrement et ses conseils;
> Monsieur Cheikh THIAM, ingnieur AIR LIQUIDE, pour sa disponibilit et ses
conseils;
> Monsieur Mamadou SARR, professeur l'ESP, pour sa disponibilit et sa.
collaboration ;
> Monsieur Papa Mody NDIAYE, pour ses conseils;
> Tous ceux qui, de prs ou de loin ont contribu l'accomplissement de ce travail,
particulirement tout le personnel de l'ESP (centre de This) et les camarades de
promotion.
SOMMAIRE
Le but de ce travail est de faire l'tude de faisabilit technique et financire d'une
centrale lectro-solaire thermodynamique devant produire de l'hydrogne trs pur par
lectrolyse de l'eau. Cette tude est effectue pour assurer la demande au Sngal et dans la
sous rgion.
Dans un premier temps, nous avons procd une tude de march au Sngal, auprs
de AIR LIQUIDE, et dans la sous rgion, travers Internet. D'aprs les rsultats de cette
tude, nous nous sommes rendus compte qu'il existe des privs mconnus dans la
commercialisation du produit. La non rglementation du commerce de l'hydrogne nous
pousse estimer la demande future 2 000 m
3
par an.
En second lieu, sur la base d'une tude hibliographique, nous avons numr toutes
les utilisations actuelles de l'hydrogne, prsent d'autres mthodes de production et fait le
choix de celle adopter parmi dautres pour sun stockage qui est nergivore.
Dans la partie technique, il tait question d'tudier l'unit afin de choisir tous les
appareils entrant dans la chane de production ct de modliser le systme. La simulation pour
le dimensionnement de la chane, initialement prvue avec fORTRAN, a t annule par
manque de temps. Nanmoins, nous avons ralis un dimensionnernent bas sur le choix de
rendements thoriques des lments du systme.
Quant l'tude financire, elle constitue la phase dans laquelle a t tudie la
rentabilit du projet travers la dtermination d'indices tels que la valeur actuelle nette
(VAN) et le dlai de rcupration des capitaux investis (DRe). Ceux-ci ont rvl que le
projet n'tait financirement pas rentable.
A la suite de ce travail, les rsultats obtenus montrent qu'un projet de cette nature
n'est pas rentable dans les pays en dveloppement o la demande en hydrogne est trop
limite. Cependant, les avantages lis au caractre non polluant de l'hydrogne restent
confirms.
Mots-cls: faisabilit - hydrogne - lectre-solaire - thermodynamique - modlisation -
rentabilit - VAN- DRC - dimensionnemcnt - pollution-
..
If
P.F.E. Prsent par M. MbayeSECKet M. MoussaDIOP
LISTE DES TABLEAUX
E.S.P. This (2006- 2007)
Tableau 1: les utilisations de l'hydrogne 6
Tableau 2 : Caractristiques actuelles des dispositifs concentration .34
Tableau 3 : caractristiques du HITEC et son stockage 38
Tableau 4 : puissance des appareils 71
iii
P.F.E. Prsentpar M. Mbaye SEeK et M. MoussaDIOP
LISTES DES FIGITRES
E.S.P. This(2006- 2007)
..
Figure 1 : Prsentation du systme 15
Figure 2 : Electrolyseur lectrodes monopolaire 19
Figure 3 : Electrolyseur lectrodes bipolaires 20
Figure 4 : Synoptique de l'alimentation stabilise 21
Figure 5 : Systmes rflecteur cylindro-parabolique 28
Figure 6 : Systme rflecteur parabolique 29
Figure 7 : Systme de centrale tour. 30
Figure 8 : Schma simplifi du systme primaire .39
Figure 9 : Calcul de concentration .41
Figure 10 : Pertes par rflexion .44
Figure Il : Schma d'un absorbeur. .45
Figure 12 : Schma simplifi du systme secondaire 57
Figure 13 : Ecoulement en prsence de chicane 62
iv
P.F.E. Prsentpar M. Mbaye SECKet M. Moussa DIOP
TABLE DES MATIERES
E.S.P. This (2006- 2007)
REMERCIEMENTS i
SOMMAIRE ~ .ii
LISTE DES TABLEAUX .iii
LISTE DES FIGURES '" .iv
TABLE DES MATIERES v
INTRODUCTION 1
CHAPITRE 1 : ETUDE BIBLIOGRAPHIQUE ~ 2
1-1 : GENERALITES 2
1-2 : LES PROPRIETES DE L'HYDROGENE 2
1-3 : LES DIFFERENTES UTILISATIONS DE L'HYDROGENE 3
1-4: LES AUTRES METHODES DE PRODUCTION D'HyDROGENE 6
1-4-1 : La production partir de combustibles fossiles 6
1-4-1-1 : Par vaporeformage 6
1-4-1-2 : Par Oxydation partielle 7
1-4-1-3 : Par Reformage auto-thermique 8
1-4-2 : La transformation thermochimique de la biomasse 8
1-4-3 : Par photolyse de l'eau 9
1-5 : LE STOCKAGE DE L'HYDROGENE 10
1-5-1 : Sous forme liquide basse temprature (liqufi) 11
1-5-2: Sous forme gazeuse sous pression (comprim) 11
1-5-3 : Sous forme solide (hydrure mtallique) 12
1-6 : CONCLUSION PARTIELLE 12
CHAPITRE II : PRODUCTION D'HYDROGENE PAR ELECTROLYSE DE L'EAU A
L'AIDE D'UNE CENTRALE ELECTRO-SOLAlRE THERMODYNAMIQUE 13
11-1 : PRESENTATION DU SYSTEME 13
11-2: L'ETUDE DU SYSTEME 16
11-2-1 : Systme tertiaire 16
11-2-1-1 : Electrolyse de l'eau 16
v
P.F.E. Prsent par M. Mbaye SECK et M. MoussaDIOP E.S.P. This(2006- 2007)
11-2-1-2: La cellule d'lectrolyse 16
a: L'lectrolyte 16
b : Les lectrodes 17
c : Le diaphragme 18
11-2-1-3 : L'lectrolyseur 18
a : Conception des appareils industriels 18
II-2-1-4 : L'alimentation lectrique 20
a : Alimentation stabilise 21
b : Le transformateur triphas 21
11-2-2 : Systme secondaire 22
11-2-2-1 : L'alternateur 22
11-2-2-2 : La turbine 25
a : Dfinition 25
b : Le fonctionnement 26
c : Les type de turbines 27
d : gnrateur de vapeur 28
II-2-3 : Le systme primaire 28
11-2-3-1 : Le concentrateur 28
11-2-3-2 : Le rflecteur cylindro-parabolique 28
11-2-3-3 : Le rflecteur parabolique 29
11-2-3-4 : Les hliostats 29
11-2-4 : Les aspects techniques 30
11-2':4-1 : Concentration 30
11-2-4-2 : La Filire parabolique 31
11-2-4-3 : La filire centrale tour (hliostats) , 32
11-2-4-4 : Filire cylindro-parabolique 33
a : Rception du rayonnement. 35
b : L'absorbeur 35
c : Le rflecteur 36
d : La concentration effective 36
e : Le fluide caloporteur 36
d : Le stockage du fluide caloporteur 37
vi
P.F.E. Prsent par M. Mbaye SECK et M. Moussa DIOP E.S.P. This (2006- 2007)
CHAPITRE III: MODELISATIONET DIMENSIONNEMENTDU SySTEME 39
III-I : MODELISATION DU SYSTEME PRIMAIRE 39
III-1-1 : Introduction 39
III-1-2 : Modle de l'ensoleillement 39
III-I-2-1 : Rappel des donnes astronomiques 40
III-1-2-2 : Flux incident du rayonnement solaire 40
III-1-3 : Modle du miroir cylindro-parabolique 41
III-1-3-1 : Calcul de la concentrationproduite par le miroir cylindro-parabolique 41
III-I-3-2 : Calcul des pertes nergtiques dues la gomtrie de la parabole 43
III-1-4 : Modle de l'absorbeur 44
a : Les pertes par rflexion 44
111-1-4-2 : 'Deuxime modle 45
a: Transfert de chaleur entre l'absorbeur et le fluide caloporteur 47
b : Transfert de chaleur entre l'absorbeur et l'enveloppe de verre 48
c : Transfert de chaleur entre l'enveloppe de verre et l'environnement.. 48
d : Diffrente de temperature entre le fluide et l'absorbeur 49
III-1-5 : Modle du dbit du HITEC 51
III-1-6: Modle de stockage 52
III-I-6-1 : Modle du stockage chaud 52
III-I-6-2 : Modle du stockage froid 54
a : Pertes du stockage 54
III-2) DU SYSTEME SECONDAIRE 57
III-2-1 : Introduction 57
III-2-2 : Modle des Echangeurs 58
1II-2-2-1 : Description du modle 58
1II-2-2-2 : Equations exprimant les changes 59
a: Coefficient d'change global.. 61
b : Coefficient d'change KH travers le film intrieur 61
c : Coefficient d'change KE travers le film externe 62
III-2-3 : Modle du gnrateur de vapeur 63
111-2-3-1 : L'conomiseur EQ 65
111-2-3-2: L'vaporateur EV 66
111-2-3-3 : Le surchauffeur ZU 66
vii
P.F.E. Prsentpar M. Mbaye SECK et M. Moussa DIOP E.S.P. This (2006 - 2007)
111-2-4 : Modle de la turbine 67
111-2-4-1 : Introduction : 67
111-2-4-2 : Equation de la chute isentropique 67
lII-2-5 : Le stockage de l'hydrogne 68
111-3 : DIMENSIONNEMENT 69
III-3-1 : Dtermination de la charge 69
111-3-1-1 : Puissance de l'lectrolyseur 69
111-3-1-2: Quantit d'eau ncessaire 70
111-3-1-3 : La puissance des moteurs pour la compression 70
111-3-Z: La charge totale 70
111-3-3 : La puissance des autres appareils 71
CHAPITRE IV: ETUDE FINANCIERE e . 73
IV-1 : ETUDE DE COUT 73
IV-1-1 : Cot des quipements de production 73
/v-1-2: Cot du matriel et mobilier de bureau 73
/v-1-3: Cot des matires premires 74
IV-1-4 : Cot du personnel 74
IV-1-5.: Cot des. quipements divers et de maintenance 74
IV-1-6 : Cot de construction Gnie civil 75
IV-2 : ETUDE FINANCIERE 75
IV-2-1 : Les investissements 75
IV-2-2 : Les projections financires 75
IV-2-3 : L 'amortissement 76
IV-2-4: Les recettes d'exploitation 76
IV-2-5 : Dtrmination du prix de vente du produit 76
IV-2-5-1 : Dtermination des cots variables 77
IV-2-5-2 : Dtermination des cots fixes 77
IV-2-6 : La rentabilit du projet 78
IV-2-6-1 : Le Dlai de Rcupration du Capital Investi (DRC) 78
IV-2-6-2 : La Valeur Actuelle Nette (VAN) 78
IV-2-7: Etudefinancire en tenant compte de l'oxygne produite par l'lectrolyse 79
IV-3 : CONCLUSION PARTIELLE 79
viii
P.F.E. Prsent par M. MbayeSECKet M. MoussaDIOP E.S.P. This (2006- 2007)
CONCLUSION GENERALE 80
BIBLIOGRAPHIES 82
1 .
ANNEXES. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. . . . . . . . . . . . . . . . . . . . . . .. 84
ix
INTRODUCTION
E.S.P. This (2006- 2007)
Depuis la rvolution industrielle, la production nergtique internationale ne cesse
d'augmenter afin de satisfaire les besoins immenses dans ce domaine. En effet, gnralement
cette production nergtique quantitative se fait avec des dgagements de gaz effet de serre
qui constitue des menaces srieuses pour la sante publique par destruction de la couche
d'ozone. Face cette situation, les gens s'entendent sur le fait que l'utilisation des sources
d'nergie renouvelable devrait remplacer les combustibles actuels qui seront bientt puiss.
L'hydrogne semble tre la source et le vecteur d'nergie idal. Il est disponible en trs
grande quantit dans l'eau. Les produits de sa combustion sont non polluants car on rgnre
l'eau. Il conduit une meilleure utilisation des ressources ce qui signifie qu'on en retire une
plus grande quantit d'nergie que la plupart des combustibles actuels. Il peut servir de
moyen de transport et de stockage d'nergie.
Pour des raisons conomiques autant qu'cologiques, l'heure de l'hydrogne en tant que
vecteur nergtique semble venue. Pour qu'il puisse s'imposer, il faudra apporter des
rponses aux questions cruciales qui se posent encore, la premire tant : comment produire
cet hydrogne en n'mettant que peu (ou pas) de gaz effet de serre?
La seconde question porte sur l'infrastructure mettre en place, jusqu'au stockage. C'est dans
ce contexte que plusieurs axes de recherches sont lancs.
Les mthodes permettant de produire de l'hydrogne sont nombreuses, mais pour rpondre.
nos questions, nous avons port notre choix sur la production par lectrolyse de l'eau. Pour
l'alimentation de notre unit en courant lectrique, nous allons tudier plus particulirement
l'nergie solaire (nergie renouvelable et inpuisable) thermodynamique qui vise
transformer l'nergie rayonne par le soleil en chaleur temprature leve, puis celle-ci en
nergie mcanique (et ,lectrique) travers un cycle thermodynamique.
Aprs une tude bibliographique, nous passerons l'tude et au choix de tous les lments de
notre unit de production qui sera subdivise en trois sous systmes, ensuite nous passerons
la modlisation et ferons un dimensionnement approximatif de chaque sous systme, afin
nous terminerons par l'tude financire.
1
P.F.E. Prsent par M. Mbaye SECKet M. MoussaDIOP
CHAPITRE 1: ETUDE BIBLIOGRAPHIQUE
E.S.P. This (2006- 2007)
1-1 : Gnralits [1]
Le dihydrogne, communment appel hydrogne, a t dcouvert ere-(9n)par Boyle.
Cavendish tudia ses proprits et Lavoisier qui au cours d'une expriencebservant que la
combustion de l'hydrogne produit de l'eau, lui donna son nom hydrogne; du grec hydro,
eau et gne producteur.
Bien qu'on le trouve en abondance dans l'univers (gaz interstellaire), il n'existe pas l'tat
naturel dans l'environnement terrestre : on le trouve plutt l'tat d'atome li d'autres
individus chimiques souvent par une liaison dite hydrogne (C, 0) et est prsent, par exemple
dans l'eau, le sucre, les protines ou les hydrocarbures.
Pourquoi 1'hydrogne?
Certaines rponses cette question relvent autant d'un optimisme excessif que d'un
pessimisme exagr. Comme souvent dans de tels dbats, la vrit se situe entre ces deux
extrmes. Elle rsulte de tendances lourdes : explosion de la demande nergtique,
tarissement des ressources fossiles et rchauffement de la plante imputable un effet de serre
excessif.
1-2 : Les proprits de l'hydrogne [1]
En phase gazeuse
Masse volumique du gaz (1,013 bar et 15 OC) : 0.085 kg/rn", gaz le plus lger
Masse volumique du gaz (1,013 bar au point d'bullition) : 1.312 kg/m"
Facteur de compressibilit (Z) (1,013 bar et 15 OC) : 1.001
Densit (air = 1) (1,013 bar et 21C) : 0.0696
Volume spcifique (1,013 bar et 21 OC) : 11.986 m
3/kg
Combustion haut pouvoir calorifique lev (essence, 45 Mj/kg) : 120 Mj/kg
En phase liquide.
Masse volumique phase liquide (1,013) bar au point d'bullition) : 70.973 kgf m
3
Equivalent gaz/liquide (1,013 bar et 15 OC) : 844 vol/vol
Point d'bullition (1,013 bar) : -252.8 "C
En phase solide
Point de fusion / -259 -c
2
Chaleur latente de fusion (1,013 bar, au point triple) : 58.158 kJ/kg
tonnamment, l'hydrogne est un mtal: lorsqu'il est sous trs haute pressions et trs basse
temprature, il se cristallise. Dans le tableau priodique des lments, il est d'ailleurs dans la
colonne des mtaux alcalins.
N'tant pas prsent l'tat solide sur terre, il n'est toutefois pas considr comme un mtal en
chimie.
Autres donnes
Gaz incolore, inodore, et insipide
Solubilit dans l'eau (1,013 bar et OC) : 0.0214 vol/vol
Concentration dans l'air: 0.00005 % vol
Temprature d'auto inflammation: 560 "C
1-3 : Les diffrentes utilisations de l'hydrogne [1] [2].
Les voitures hydrogne peuvent tre une faon importante de rsoudre notre dpendance au
ptrole et de stopper le rchauffement de la plante. Dans le futur, les nergies renouvelables
(nergies olienne, solaire, hydraulique) seront produites puis stockes dans de l'hydrogne et
utilises dans des voitures pile combustible. En outre l'utilisation de rservoirs
d'hydrogne et de voitures piles combustible fournit aux constructeurs automobiles de
nouvelles possibilits dans le design des voitures et leur production.
Une voiture hydrogne utilise de l'hydrogne la place de l'essence. L'hydrogne est stock
dans un rservoir soit sous pression, sous forme liquide ou dans des hydrures mtalliques.
Depuis le rservoir, l'hydrogne est conduit la pile combustible o il est converti en
lectricit grce un procd lectrochimique. L'lectricit est ensuite transmise aux moteurs
lectriques dans les roues o il nourrit le moteur.
Les voitures l'hydrogne ont une efficience suprieure aux machines combustion
traditionnelles. Habituellement, l'efficience des machines combustion est autour de 20-25%,
alors qu'une pile combustible PEM ( membrane changeuse de protons) aura une efficience
jusqu' 60% suprieure. Une pile combustible produit de l'eau qui est la seule mission, ce
qui en fait une technologie nergtique propre et efficiente, spcialement si l'hydrogne a t
produit partir d'nergies renouvelables.
Le challenge le plus important des voitures hydrogne est le stockage de l'hydrogne. Un
3
P.F.E. Prsent par M. Mbaye SECK et M. Moussa DIOP E.S.P. This (2006 - 2007)
autre challenge pour les voitures pile combustible est le fait qu'il n'y a que quelques
stations service. C'est le problme de la "poule et de l'uf'. Les compagnies nergtiques ne
construiront pas de station service hydrogne tant qu'il n'y aura pas de voiture fournir en
hydrogne, et les consommateurs n'achteront pas de voitures hydrogne tant qu'ils ne
pourront pas "faire le plein" partout o ils conduisent. Pour dvelopper l'utilisation de
l'hydrogne dans le transport routier, il y a eu beaucoup d'initiatives comme l'autoroute de
l'Hydrogne. En Californie le but est d'tablir 50-100 stations service hydrogne avant 2010.
Avec le Lien Hydrogne les pays du nord (le Danemark, la Sude, la Norvge et le Nord de
l'Allemagne) seront relis ensemble par un rseau tablissant un couloir de station service
hydrogne.
L'hydrogne est utilis principalement:
~ pour ses proprits d'agent rducteur dans l'industrie ptrolire pour la dsulfuration
de l'essence et du gazole (37% de la production totale) ;
~ dans l'industrie chimique pour la fabrication d'ammoniaque suivant la raction
chimique exothermique en phase gazeuse (50%) :
de mthanol (8%), de polymres comme les polyurthanes et les polycarbonates ;
~ dans l'industrie du verre ou des semi-conducteurs et dans la mtallurgie pour la
transformation de l'acier.
Ces informations sont ainsi regroupes dans le tableau 1.
industrie applications
L'utilisation la plus importante de l'hydrogne se trouve
dans la synthse de l'ammoniac (50 %) suivant la
chimie
raction chimique exothermique en phase gazeuse
N2 + 3H2 ~ 2NH
3
Alimentaireet boissons
L'hydrogne pur est utilis pour la production de
plastiques, polyester et nylon. H2 est aussi utilis dans
4
P.F.E. Prsent par M. Mbaye SECKet M. MoussaDIOP E.S.P. This(2006- 2007)
Verre, chaux et ciment
Industrie des mtaux
Laboratoires et analyses
l'hydrognation des ammes et acides gras (huiles
alimentaires)
L'hydrogne est principalement utilis pour la cration
en combinaison avec l'azote d'une atmosphre
rductrice au-dessus du bain d'tain dans le procde
FLOAT. On utilise l'hydrogne pour le traitement
thermique (flamme oxy -hydrogne) du verre creux et
des performes en fibre optique.
Atmosphre rductrice pour des procds de traitement
thermique. 1
L'hydrogne est utilis comme gaz porteur en
chromatographie en phase gazeuse et dans de
nombreuses techniques analytiques. Les plus communes
sont l'utilisation dans les' flammes des dtecteurs
ionisation de flamme (FID) ou des dtecteurs
photomtrie de flamme (FPD). On utilise aussi des
mlanges contenant de l'hydrogne dans les appareils
d'tincelage et dans la mesure des hydrocarbures
totaux.
Soudage, coupage et projection Traitement thermique de divers mtaux
thermique
Ptrole et divers
lectronique
Espace et aronautique
Dsulfuration du fuel et de l'essence
L'hydrogne est gnralement utilis comme gaz de
balayage lors des tapes de dpt de silicium ou de
croissance cristalline, aussi bien que dans le brassage
sous atmosphre et le recuit de film de cuivre. L'emploi
d'hydrogne avec des gaz inertes protecteurs (H2 dilu
dans de l'azote) permet thoriquement d'liminer la
totalit de l'oxygne ainsi que ses inconvnients dans
un milieu essentiel pour les procdes haute ou
moyenne tempratures.
L'hydrogne est utilis sous forme liquide dans les
ergols pour la propulsion des tages cryogniques des
5
P.F.E.Prsent par M. MbayeSECKet M. MoussaDIOP E.S.P. This (2006- 2007)
fuses Ariane.
Avec l'hydrogne on fabrique de l'eau oxygne H
202
divers utilise comme dcolorant ou dsinfectant ou pour le
blanchiment de la pte papier.
Tableau 1 : Les utilisations de l'hydrogne
1-4 : Les autres mthodes de production d'hydrogne [2]
Honnis l'lectrolyse de l'eau, il existe d'autres mthodes de production d'hydrogne que
nous avons notes :
1-4-1: La production partir de combustibles fossiles
Dans le pass, la production industrielle de l'hydrogne partir des hydrocarbures s'est
oriente vers les options les plus performantes d'un point de vue conomique. Sans, surprise,
la matire premire utilise de faon dominante est, de loin, le gaz naturel, constitu
principalement de mthane. D'autres nergies fossiles sont parfois utilises puisqu'elles
contiennent aussi de l'hydrogne qui peut tre produit par la raction d'oxydation partielle,
toutes ces solutions qui partent des nergies fossiles passent par la production de gaz de
synthse, mlange d'oxyde de carbone (CO) et d'hydrogne (H2).
Les technologies dominantes aujourd'hui pour la production de gaz de synthse partir des
hydrocarbures sont le vaporefonnage (reformage la vapeur), l'oxydation partielle et le
reformage auto thermique, combinaison des deux prcdentes, utilis surtout pour la
production de carburant liquide de synthse.
/-4-1-1 : Par vaporeformage
Le reformage la vapeur' consiste transformer les charges lgres d'hydrocarbure en gaz de
synthse (mlange H2, CO, C02, CI-4, H20 et autres impurets) par raction avec la vapeur
d'eau sur un catalyseur au nickel. Cette transformation a lieu haute temprature (840
950C) et pression modre (de l'ordre de 20 30 bars). Elle peut tre suivie par diffrentes
oprations qui conduisent la production d'hydrogne mais donc aussi l'obtention de
6
P.f.E. Prsent par M. Mbaye SECKet M. MoussaDIOP E.S.P. This (2006- 2007)
carburant de synthse. Dans tous les cas, la charge d'une unit d vaporeformage peut tre du
gaz naturel (charge de rfrence) ou du mthane.
~ Le vaporeformage du gaz naturel
Le gaz naturel contient essentiellement du mthane. Toutefois, il doit gnralement tre
dsulfur avant d'tre dirig vers l'unit de vaporeformage. Pour maximiser la production de
l'hydrogne, les deux principales ractions chimiques mettre en uvre sont la production de
gaz de synthse et la conversation du CO.
Reaction 1:"CI-Lt + H
2
0 CO + 3H2
Reaction 2: CO+ H20 CO
2
+H
2
Bilan des deux ractions: CH4 + 2H
20
C02 + 4H2
La premire raction correspond au vaporeformage proprement dit. Elle est endothermique et
se caractrise par un rapport H
2
/CO de l'ordre de 3. La deuxime raction correspond la
conversion du CO. Elle est lgrement exothermique et plus ou moins complte, selon qu'elle
est effectue en 10u 2 tapes. Globalement, le bilan des deux ractions est endothermique.
Les deux dernires tapes consistent sparer le CO
2
et l'hydrogne puis liminer les
dernires traces d'impurets.
/-4-1-2 : Par Oxydation partielle
L'oxydation partielle peut tre effectue sur un ventail de produits beaucoup plus large allant
du gaz naturel aux rsidus lourds, et mme au charbon. D'un point de vue conomique, le
traitement des charges lourdes se justifie quand le surinvestissement consenti est compens
par un cot rduit de la matire premire, le coke de ptrole par exemple. En revanche, les
investissements et les conditions opratoires sont plus contraignants.
A haute temprature (1200 1500 OC) et pression leve (20 90 bars ou plus), en prsence
d'un oxydant (l'air ou plus couramment l'oxygne) et d'un modrateur de temprature (la
vapeur d'eau), l'oxydation partielle des hydrocarbures conduit, l'instar du vaporeformage,
la production de gaz de synthse. En revanche la raction est exothermique et se droule (en
gnral) sans catalyseur.
~ Oxydation partielle de rsidu sous vide
Les deux principales ractions sont la production de gaz de synthse (raction globale
simplifie) et la conversion du monoxyde de carbone
7
Un rsidu sous vide type contient plus d'une cinquantaine d'atomes de carbone. Il est
symbolis ici par la chane hydrocarbone CnH
m
qui correspond une raction HIC gal
rn/n.
Raction 1 : CnH
m
+ ( ~ ) 0 2 nCO + (m)H
2
2 2
Raction 2 : nCO + nH
20
nC02 + nH2
Bilan des deux ractions: CnH
m
+ nH
20
+ ( ~ ) 0 2 - nC02 + (n + m)H
2
2 2
La premire raction (temprature comprise entre 1300 et 1400 OC) correspond l'oxydation
partielle proprement dit. Elle est exothermique et se caractrise par un rapport HlCO de
l'ordre de 0.75 (pour un rapport HIC de l'ordre de 1.5). Elle reprsente le rsultat global d'un
grand nombre de ractions qui conduisent la formation de CO et d'H
2
mais aussi de : H
20,
C et CH.
1-4-1-3: Par Reformage auto-thermique
Le procd auto-therme de reformage auto-thermique est probablement la plus intressante
des prochaines volutions des procds de production d'hydrogne partir d'hydrocarbure.
Son principal atout est de compenser (dans un mme racteur lit fixe) les ractions
endothermiques du vaporeformage par les ractions exothermiques de l'oxydation partielle.
Ce procd suscite beaucoup d'intrt pour les grosses capacits et notamment la filire GtL
(filire qui permet de convertir le gaz naturel en produit liquide ptrolier et notamment en
carburant diesel) car le rapport COIHpeut tre ajust avant entre dans le racteur. Comme
pour l'oxydation partielle, le cot de l'oxygne est ici un paramtre important. Grce des
membranes cramiques qui restent valider en production industrielle, il serait possible
d'viter la dispendieuse tape du fractionnement d'air. La filire GtL est riche de potentiel car
elle offre une source alternative de trs bonne qualit aux carburants conventionnels issus du
ptrole et est aussi une voie alternative la filire gaz naturel liqufi pour la mise en
exploitation.
1-4-2: La transformation thermochimique de la biomasse.
La biomasse est constitue de tous les vgtaux qui se dveloppent la surface de la terre.
Elle effectue le captage et le stockage de l'nergie solaire qu'il est ensuite possible de restituer
sous forme de combustion ou de carburant, sans impact sur l'effet de serre, par des
transformations plus ou moins efficaces aux plans nergtique et conomique.
La transformation thermochimique, qui conduit la gazification des
8
substances organiques et vgtales, convient particulirement bien la valorisation des
produits lignocelh;llogiques comme le bois et les pailles.
Procd de gazification de la biomasse
Il comporte une succession d'oprations qui ncessitent simultanment le transfert de grandes
quantits de chaleur et le contrle de la proportion et du temps de contact des ractifs en
prsence.
Apres l'opration de schage trs endothermique, la thermolyse consiste en la dgradation
thermique des produits qui aux alentours de 600C ont perdu 70 % de leur masse devenue
gazeuse alors que 30 % de la masse d'origine demeure solide et reprsente le charbon de bois,
essentiellement constitu de carbone (C). La gazification du carbone se poursuit jusque vers
900 et 1000 OC l'aide d'un ractif comme l'air, l'oxygne (Oz) ou de la vapeur d'eau (HzO),
qui ralise une oxydation plus ou moins complte du carbone en monoxyde de carbone (CO)
ou en gaz carbonique (COz). L'utilisation d'air, ce stade, est peu coteuse mais introduit
dans le produit final de l'azote indsirable. L'oxygne provenant de la distillation de l'air, est
cher est dangereux. L'utilisation d'oxygne ou d'air favorise la formation d'oxyde de
. .
carbone. Avec la vapeur d'eau en tant que ractif, une quantit supplmentaire d'hydrogne
(Hz) est rcupre. Toutefois, un moyen de chauffage annexe doit tre employ, par exemple
par combustion de carbone ou de gaz de purge rsiduels ou encore un apport d'nergie
externe. L'obtention d'un gaz de bonne qualit et non corrosif ncessite l'limination des
acides et des goudrons rsultant de la thermolyse. Il faut pour cela porter ces produits aux
environs de 1200 ou 1300 "C ( tage haute temprature) afin de raliser leur dgradation
thermique (craquage) ou effectuer vers 800 900 "C une raction catalyse qui permet
d'viter la fusion des cendres et la cration de mchefers constates entre 900 et 1000 "C.
L'utilisation d'un tage haute temprature permet aussi de raliser simultanment le
reformage du mthane (CH4) rsiduel en CO et Hz. Le gaz obtenu peut tre raffin pour en
extraire l'hydrogne qui deviendra par cette filire le combustible propre et renouvelable.
1-4-3: Par photolyse de l'eau
Photolyse (dcomposition chimique d'un corps par la lumire)
Des cellules photo lectrochimiques (PEC) illumines par la lumire solaire peuvent
dcomposer l'eau en hydrogne et en oxygne. De telles cellules emploient des lectrodes
photo actives immerges dans un lectrolyte aqueux ou dans l'eau. Sous leur forme la plus
9
P.F.E. Prsent.par M. MbayeSECKet M. MoussaDIOP E.S.P. This (2006- 2007)
simple, elles peuvent tre dcrites comme des dispositifs photovoltaque/lectrolyte intgrs
ou monolithiques. Une large varit de processus PEC et de photo catalyse susceptibles de
dissocier l'eau sont tudis de par le monde.
Bien qu'intrinsquement simples, les aspects scientifiques et de dveloppement des matriaux
de cellules PEC pour la production d'hydrogne en sont encore au stade prliminaire mais
fascinant de la recherche fondamentale.
Quatre tapes principales sont mises en jeu dans le processus photo lectrochimique de
dissociation de l'eau. La premire est la gnration d'une charge lectronique la surface de
la photoanode soumise au rayonnement solaire, produisant des paires lectron-trou. La
deuxime est l'oxydation de l'eau la photoanode par les trous, produisant des molcules
d'oxygne (Oz). La troisime tape est celle du transport des ions (H+) et des lectrons de la
photoanode la cathode, respectivement via l'lectrolyte et la connexion lectrique (la
cathode pouvant tre une photocathode induite par la lumire). Enfin, la rduction des ft la
cathode forme des molcules d'hydrogne (Hz) l'aide des lectrons.
Le choix des matriaux de type n de la photoanode est dterminant pour les cellules PEC,
avec deux critres : l'existence d'une bande interdite parfaite et la possibilit d'viter la
recombinaison des porteurs de charge. L'nergie minimale thorique du photon devant tre
absorb par une photoanode parfaite est de 1.23 lectronvolt (eV). C'est la force
lectromotrice minimale thorique ncessaire aux cellules PEC pour dissocier la molcule
d'eau.
15 : Lestockage de l'hydrogne [3]
L'hydrogne est un produit trs volatil et inflammable, ce qui explique les difficults
rencontres pour son stockage. En effet, il doit offrir un haut degr de scurit et d'autre part
des facilits d'usage en terme de densit d'nergie et de dynamique de stockage/dstockage
pour permettre aux diffrentes applications de fonctionner dans des conditions techniques
acceptables.
Aujourd'hui, il existe principalement trois modes de stockage. Ces diffrentes modes
prsentent chacun leurs avantages et leurs inconvnients selon des critres conomiques,
nergtiques, de capacit massique et volumique, de scurit de stockage/dstockage.
10
P.F.E. Prsent par M. Mbaye SECKet M. MoussaDIOP E.S.P. This (2006- 2007)
1-5-1: Sous forme liquide basse temprature (liqufi)
Conditionner l'hydrogne sous forme liquide est une solution a priori attrayante. C'est
d'ailleurs sous cette forme qu'il est utilis dans le domaine spatial. Mais l'hydrogne est,
aprs l'hlium, le gaz le plus difficile liqufier. Cette voie est sduisante en termes de
capacit mais pose des problmes nergtiques importants et des cots levs qui rendent son
application plus difficile pour le grand public. En particulier, l'nergie ncessaire pour la
liqufaction pnalise gravement le rendement global de ce mode de stockage.
1-5-2: Sous forme gazeuse sous pression (comprim)
Le conditionnement de l'hydrogne sous forme gazeuse est une option prometteuse. Les
contraintes sont toutefois nombreuses.
Leger et volumineux, l'hydrogne doit tre comprim au maximum pour rduire
l'encombreinent des rservoirs. Toutefois, pour atteindre des capacits comptitives, les
rservoirs doivent rsister de trs fortes tensions jusqu' 700 bars. Mais cette compression a
un cot.
Le risque de fuite d'hydrogne doit tre galement pris en considration compte tenu du
caractre inflammable et explosif de ce gaz dans certaines conditions. Or en raison de la petite
taille de sa molcule. L'hydrogne est capable de traverser de nombreux matriaux, y compris
certains mtaux. Il en fragilise, de plus, certains en les rendant cassants. L'tude du stockage
haute pression consiste donc, pour l'essentiel, prouver la rsistance des matriaux
l'hydrogne sous pression. Ces matriaux doivent tre rsistants mais relativement lgers
(mobilit oblige).
Les rservoirs mtalliques, utiliss actuellement, se rvlent encore coteux et lourds au
regard de la quantit de gaz qu'ils peuvent emporter. Des rservoirs non mtalliques mais en
matriaux polymres sont en cours d'laboration pour rpondre ces contraintes.
En consquence, les procds de liqufaction comprennent obligatoirement les phases
suivantes:
purification et compression de l'hydrogne,
pr refroidissement (cascade classique),
conversion ortho-para l'aide de catalyseurs (oxyde de chrome),
liqufaction,
stockage en rservoir.
11
P.F.E. Prsent par M. Mbaye SECKet M. MoussaDIOP E.S.P. This (2006- 2007)
1-5-3: Sous forme solide (hydrure mtallique)
Cette mthode repose sur la mise au point de matriaux haute capacit de stockage :hydrures
mtalliques rversibles dans des conditions de temprature et de pression voisines des
conditions ambiantes et matriaux carbons (nanostructures, charbon activs etc.... ) pour le
stockage de l'hydrogne par adsorption/compression pression rduite. D'autres matriaux
tels que les composs molculaires large cavit, prsentent galement des performances
intressantes, Les principales difficults concernent la mise au point de ces matriaux et leurs
systmes associs de stockage satisfaisant aux critres principaux dfinis prcdemment. De
nombreuses recherches sont menes pour dvelopper la conception des cots raisonnables.
Une meilleure comprhension des phnomnes physico-chimiques d'interaction entre
hydrogne et matriaux s'avre ncessaire. D'autre part, la matrise de la mise en uvre de
ces matriaux dans les, rservoirs adapts doit tre dvelopp.
De nombreux lments peuvent fixer l'hydrogne de faon rversible en quantit
intressantes.
Les critres de slection pour un bon hydrure sont:
absorption et dsorption de H2 temprature et pression modres,
dure de stockage peu prs illimite,
manipulation facile,
capacit de stockage (d'absorption) leve, mme aprs un nombre de cycles,
sous pression modre, les hydrures absorbent l'hydrogne; en heure de pointe, on rduit la
pression, augmente la temprature et l'on dgage ainsi l'hydrogne.
1-6 : Conclusion partielle
L'hydrogne est trs difficile stocker, d'normes moyens sont utiliss pour ne stocker
que de faibles quantits. C'est pour quoi les constructeurs d'automobiles ont d'normes
problmes pour le suffisamment et surement.
Son utilisation dans le transport n'est plus considre comme impossible puisque en Irlande
30% des bus fonctionnent actuellement avec de l'hydrogne
C'est un produit trs utilis dans le monde surtout pour la production de l'ammoniac. En effet
50% de la production mondiale d'hydrogne est utilise pour produire de l'ammoniac.
12
CHAPITRE II : PRODUCTION D'HYDROGENE PAR
ELECTROLYSE DE L'EAU AL'AIDE D'UNE CENTRALE ELECTRO-
SOLAIRE THERMODYNAMIQUE
11-1: Prsentation du systme
Le systme peut tre divis en trois sous systmes selon l'tat de l'nergie, il se prsente
comme suit:
Le systme primaire
Dans ce systme l' riergie est essentiellement thermique. Les rayonnements solaires sont
concentrs dans un rcepteur qui transmet cette chaleur un fluide caloporteur. Ce fluide
transmet son nergie thermique au systme secondaire par l'intermdiaire d'un changeur de
chaleur.
Le systme secondaire
L'nergie cde par le systme primaire est recueillie et transforme ce stade en nergie
thermodynamique. Celui-ci nous permet d'obtenir une nergie mcanique via une turbine. Un
alternateur transforme cette nergie en lectricit et la transmet au systme tertiaire.
Systme tertiaire
L'nergie lectrique reue est transforme en nergie chimique par un lectrolyseur. Nous
obtenons ainsi de l'hydrogne et de l'oxygne la sortie.
Le schma suivant nous rsume les diffrentes phases de notre systme.
13
P.F.E. Prsent par M. Mbaye SECK et M. Moussa DIOP E.S.P. This (2006 - 2007)
Il
Il
Il
Il
.---...... 1
Il IL _
1
1
Condenseur
1
1
1
1
1
1
Pompe
G.S i----"------
r - - - - - - - - - - - - - - - - - ~ - - - - - I I
Il
Il
l'
Il
l'
1
Pompe
Stock chaud froid
Pompe
Stock
1
1
1
L------__<ll-----------II
1
1
1
_______________ ~ ~ ~ - - - - - - - - - - - - - - - - - - - 1
r - - - - - - - - - - ~ - - - - - - - - - - - ~ - - - - - - - - - - I
1
1
1
1
1 Systme primaire 1
1 Systme secondaire 1
1 Systme tertiaire
Figure 1 : Prsentation du systme
15
P.F.E. Prsent par M. MbayeSECKet M. MoussaDIOP , E.S.P. This (2006-2007)
1
112 : L'tude du systme
Le choix des lments du systme dpend de la quantit d'hydrogne produire et de leurs
cots. Pour plus de cohrence, nous commenons l'tude du systme par le systme tertiaire,
puis le systme secondaire et terminer avec le systme tertiaire.
11-2-1: Systme tertiaire
1 I ~ 2 - 1 - 1 : Electrolyse de l'eau
L'lectrolyse de l'eau consiste faire passer un courant lectrique dans une saumure
d'hydroxyde de potassium. C'est un procd aujourd'hui coteux de production d'hydrogne.
Il peut tre mis en uvre industriellement, dans les installations de plusieurs mgawatts,
lorsque l'lectricit est disponible faible cot. Elle pourrait se dvelopper aussi bien dans de
petites installations utilisant des nergies renouvelables [4].
11-2-1-2 ': La cellule d'lectrolyse
Une cellule d'lectrolyse est constitue de deux lectrodes (anode et cathode) conducteurs
lectroniques relies un gnrateur de courant continu et d'un lectrolyte (milieu conducteur
ionique).
a : L'lectrolyte
C'est une solution aqueuse conductrice de courant lectrique laquelle on peut faire subir une
dcomposition.
L'lectrolyte de l'eau est en gnral une solution aqueuse d'hydroxyde de potassium dont la
concentration variable en fonction de la temprature de fonctionnement est:
~ de 25 % en masse pour une temprature de 80 90 OC
~ de 30 35 % 120 C
~ de 40 % 160 -c.
Ces concentrations correspondent la conductivit lectrique maximale de la solution la
temprature considre. Le choix se porte de prfrence sur des solutions de potasse au lieu
de solution de soude, car :
./ Pour une temprature donne leur conductivit est suprieure celle des solutions de
soude;
./ Elles ont en gnrale une teneur plus faible en impuret chlorure lie la fabrication
d'hydroxyde;
16
P.F.E. Prsentpar M. MbayeSECKet M. MoussaDIOP E.S.P. This (2006-2007)
./ L e u ~ pression' de vapeur saturante est inferieure pour une mme temprature, ce qui
permet de rduire le dbit d'eau de rfrigration destin condenser la vapeur
entraine par l'hydrogne [4].
~ Qualit de l'eau d'alimentation
L'eau introduite dans l'lectrolyseur doit tre la plus pure possible car les impurets
demeurent dans l'appareil et s'accumulent au fil de l'lectrolyse. Elles finissent par perturber
les ractions lectrolytiques:
./ Formation de boues;
./ Action des chlorures qui dtruisent la protection par nickelage ou le nickel massif des
anodes.
Aprs une premire puration, l'eau traverse un filtre charbon actif, puis subit une
dminralisation totale par passage sur un filtre couche mixte changeur d'ions.
~ Qualit de l'hydroxyde de potassium
Les concentrations en chlorure et en sulfate doivent tre inferieures 0.01 % en masse. Il est
parfois ajout directement du silicate de sodium ou de potassium en concentration
correspondant l'quilibre de la dissolution de l'amiante la temprature considre. Il est
souhaitable d'viter les renouvellements trop frquents de la solution d'hydroxyde [4].
. b : Les lectrodes
une bonne lectrode doit tre peu coteuse, avoir une bonne rigidit mcanique, tre rsistante
la corrosion, avoir une bonne stabilit long terme mme lorsque l'on inverse la polarit du
systme ou qu'on laisse le systme en circuit ouvert, avoir une bonne conductivit lectrique
et thermique, tre sensible aux poisons, montrer une grande surface relle, avoir une bonne
activit catalytique et tre relativement facile fabriquer.
~ L'anode
C'est l'lectrode par o entre le courant lectrique, elle est relie la borne positive du
gnrateur. L'oxygne est recueilli du cot de l'anode.
~ La cathode
C'est l'lectrode par o sort le courant lectrique; elle est relie la borne ngative du
gnrateur. C'est elle qui recueille l'hydrogne.
17
E.S.P. This(2006-2007)
En milieu acide, on admet gnralement que la dcomposition de l'eau s'effectue suivant les
ractions:
la cathode :
2H++ 2e- - H2
l'anode:
2H20 - 20H + 2H++ Ze'
Puis:
1
2 0 H - H 2 0 + ~ 2
En solution alcaline, par contre, les ions OH" prdominent et le phnomne est habituellement
exprim par :
la cathode
2H20 + 2e- - H2 +20H-
l'anode
20H- - .!:02 + 2e- + H
20
. 2 .
Dans les deux cas, la raction globale s'crit:
1
H20-H2+-02
2
Avec ~ H = 285 kJ/mol H20 enthalpie de dissociation de l'eau [4].
e : Le diaphragme
Membrane poreuse dispose entre l'anode et la cathode, le diaphragme spare les bulles
d'hydrogne des bulles d'oxygne en permettant le passage des ions. Le diamtre moyen des
pores demeure de prfrence inferieure 1 um. Il est ncessairement constitu d'un matriau
hydrophile [4].
11-2-1-3: L'lectrolyseur
a : Conception des appareils industriels
Les premiers appareils d'lectrolyse disposaient d'lectrodes monopolaire (figure 2) :
chaque anode est connecte au ple plus, chaque cathode au ple moins. Cette disposition
conduit utiliser des transforedresseurs fonctionnant sous basse tension et forte intensit, de
cot onreux et engendrant des chutes ohmiques de contact importantes. Une premire
amlioration consiste mettre en srie des cuves d'lectrolyse, mais la liaison de cuve cuve
ncessite toujours des conducteurs de forte section.
18
P.f.E. Prsent par M. Mbaye SECK et M. Moussa DIOP E.S.P. This (2006 -2007)
Le perfectionnement ultime est obtenu grce aux lectrodes bipolaires (figure 3) : une cellule
d'lectrolyse est dfinie par deux lectrodes se faisant face et constituant les fonds. Chaque
lectrode est cathode sur une face et anode sur l'autre face. Les liaisons externes de cuve
cuve sont supprimes, la conduction s'opre l'intrieur de l'lectrode au travers de son
paisseur: la chute ohmique devient minimale. Cette conception se retrouve aujourd'hui sur
tous les appareils industriels. Cette conception introduit cependant une difficult : l'lectrode
prsente une face en milieu oxydant (l'anode) et une face en milieu rducteur (la cathode). En
milieu basique, le choix du matriau d'lectrode devient alors fort limit. Jusqu' 90 "C,
l'acier doux muni d'un nickelage soign du ct anodique convient parfaitement; au-del, les
nickelages ne prsentent pas une dure de vie suffisante, il devient ncessaire de faire appel
au nickel massif [4].
e ----.....-----.....-----..--------.
G-.....---I....--..-- --..----.----....--....lr--f-I
=- ,=-- -- ==- . ~ - ,
- - - -"
- -
-
.......
-_.
Cathodes
~
Anodes
.............np...........
':Iec::trolyte
---
- ~ - ~
--
-
-
-
-
-
--
-
--
~ ) montage.n ........
Figure 2 : Electrolyseur lectrodes monopolaire
19
P.F.E. Prsent par M. MbayeSECKet M. MoussaDIOP E.S.P. This (2006-2007)
lectrolyte
r..
les flOhea V$rticalea svmbotisem la circulation de 1'"
Figure 3 : Electrolyseur lectrodes bipolaires
b : Contrle des tempratures et des concentrations
Un lectrolyseur industriel fonctionne en gnral sous une tension comprise entre 1,8 et 2,2 V,
situe au-dessus de la tension thermoneutre. La consommation nergtique est comprise entre
4,5 et 5 kWhlNm3 d'hydrogne. Le rendement nergtique est donc toujours infrieur 1 et il
convient d'liminer en permanence la chaleur dgage lie aux irrversibilits. De plus, les
ractions qui s'oprent aux lectrodes engendrent des dsquilibres de concentration qu'il
convient de compenser. La circulation de l'lectrolyte dans tous les compartiments anodiques
et cathodiques assure le refroidissement des cellules et le maintien des concentrations
ioniques. L'cart de temprature entre l'entre et la sortie des cellules d'lectrolyse est
gnralement compris entre 5 et 15 "C. Dans ces conditions, l'cart de concentrations
ioniques entre l'anolyte et le catholyte demeure limit.
Pour un souci d'quilibre au niveau des intensits, nous allons utiliser trois cellules
d'lectrolyte en parallle et chacune est relie une phase de l'alternateur [4].
11-2-1-4 : L'alimentation lectrique
Il sera compos de deux parties :
L'alimentation stabilise: elle nous permettra d'alimenter notre lectrolyseur
qui consomme une tension continue de l'ordre de 1,8 2,2 volts,
Le' transformateur triphas_: comme son nom l'indique, il permet de
transformer la haute tension en tension d'utilisation 220/380.
20
..
a : Alimentation stabilise
Une alimentation est un montage transformant la tension alternative du secteur en une tension
continue basse tension.
Une alimentation secteur est compose d'un transformateur, d'un redresseur, d'un filtre et
d'une stabilisation/rgulation selon les besoins.
C'est un dispositif lectronique destin raliser la transformation alternatif-continu. Cette
transformation trouve son importance en la possibilit qu'elle offre d'alimenter, partir du
mme rseau de distribution lectrique, la fois des machines courant continu et des
machines courant alternatif. Il n'est plus ncessaire d'avoir deux sources de production
lectrique, ce qui est particulirement intressant pour des questions de cot mais aussi de
poids et d'encombrement. Le dveloppement des semi-conducteurs de puissance a largement
contribu l'expansion de ces dispositifs qui autorisent la transition de faon sure, avec des
rendements excellents tout en restant peu encombrants et peu onreux.
Le principe de fonctionnement consiste une modification priodique du circuit lectrique
entre les connections d'entre (rseau) et de sortie (rcepteur) du dispositif redresseur, de
faon recueillir en sortie des tensions et des courants d'ondulations suffisamment faibles
pour tre ngliges.
Abaisseur Redresseur Filtrer Rguler
Figure 4 : Synoptique de l'alimentation stabilise
b : Le transformateur triphas
Tout comme sur les lignes monophases, on utilise des transformateurs pour lever ou
abaisser la tension des lignes triphases.
Pour changer la tension d'un systme triphas, on peut fort bien utiliser trois transformateurs
monophass; chaque primaire tant li une phase mais il est plus avantageux d'utiliser un
seul transfo triphas. On gagnerait la fois sur l'encombrement et sur le poids de fer utilis.
Il est constitu d'un circuit magntique trois (3) noyaux et chaque noyau supporte un
enroulement primaire et un enroulement secondaire.
21
P.F.E. Prsent par M. Mbaye SECK et M. Moussa DIOP
E.S.P. This (2006 -2007)
Les trois enroulements primaires du transfo triphas peuvent tre couples en toile ou en
triangle ; tandis que les trois enroulements secondaires peuvent tre couples en toile, en
triangle ou en zigzag.
Il transforme la tension la borne de l'alternateur qui est de 24 kY en 380 Y, avec un
rendement de 95 %.
11-2-2:Systme secondaire
11-2-2-1 : L'alternateur
Un alternateur est une machine lectromagntique rotative fournissant des tensions de
frquences proportionnelles leur vitesse de rotation qui convertit l'nergie mcanique
fournie par un moteur ( turbine pour notre cas) en nergie lectrique courant alternatif. Ces
machines sont moins coteuses et ont un meilleur rendement que les machines courant
continu (dynamo) qui dlivrent des tensions continues (95 % au lieu de 85 %).
L'alternateur est aussi appel Gnrateur Synchrone. Pour cette raison, les lettres
conventionnelles de reprages sont G.S.
Pour alimenter notre charge (lectrolyseur et moteurs) nous aurons besoin d'un alternateur
triphas. Il est constitu de trois enroulements monophass, dcals d'un tiers de l'intervalle,
compris entre deux ples conscutifs de mme nom.
Cette machine est constitue essentiellement de deux parties: une partie fixe appele stator
qui a pour rle de produire une f..m. induite (induit) et une partie mobile appele rotor
(inducteur) destin produire le flux magntique.
Principe
Les alternateurs commerciaux sont construits avec un inducteur fixe ou un inducteur rotatif.
L'inducteur est compos de deux ou de plusieurs ples produisant un champ magntique
constant.
Un alternateur inducteur fixe : les ples saillants produisent le champ magntique qui est
coup par les conducteurs situs sur l'induit. L'induit porte un enroulement triphas dont les
bornes sont connectes trois bagues montes sur l'arbre. Un groupe de balais recueille la
tension triphase qui est applique la charge.
La valeur de la tension triphase dpend de la vitesse de rotation et de l'intensit du champ
magntique. La frquence dpend de la vitesse et du nombre de ples de l'inducteur. Les
alternateurs inducteur fixe sont utiliss pour des puissances infrieures 5 KYA. Pour des
22
P.F.E. Prsent par M. Mbaye SECK et M. Moussa DIOP E.S.P. This (2006 -2007)
puissances plus importantes, il est plus conomique, plus scuritaire et plus pratique
d'employer.un inducteur tournant.
Un alternateur inducteur tournant: cette construction est plus avantageuse car elle
permet d'alimenter directement le circuit d'alimentation sans passer par les bagues de fortes
dimensions qui seraient requises avec un induit tournant. De plus, l'isolement des bobinages
du stator est grandement simplifi du fait qu'ils ne sont soumis aucune force centrifuge.
Une gnrateur courant continu, appele excitatrice, habituellement monte sur le mme
arbre que l'alternateur, fournit le courant d'excitation aux lectro-aimants inducteurs.
Nombre de ples: Le nombre de ples d'un alternateur est impos par la vitesse du rotor et
par la frquence du courant produire. Ainsi, considrons un conducteur de l'induit devant
. .
lequel se dplacent les ples nord et les ples sud du rotor. Si la tension induite dans ce
conducteur prend une srie de valeurs positives pendant le passage d'un ple nord, elle
prendra une srie de valeurs gales, mais ngatives, pendant le passage d'un ple sud. Chaque
fois qu'une paire de ples se dplace devant un conducteur, la tension induite dcrit donc un
cycle complet. On en dduit que la frquence est donne par l'quation:
pn
f = 120
O:
f =frquence de la tension induite [Hz]
P = nombre de ples du rotor
n = vitesse du rotor [tr/min]
. .
Stator: Le stator est l'induit. Il est constitu d'enroulements qui vont tre le sige de courant
lectrique alternatif induit par le champ magntique en mouvement de l'inducteur.
Il se compose d'un noyau feuillet ayant la forme d'un cylindre vide et comportant des
encoches dans lesquelles sont logs les conducteurs d'un enroulement triphas.
L'enroulement est toujours raccord en toile et le neutre est accessible pour permettre sa
mise la terre. On prfre la connexion en toile celle en triangle pour les raisons
suivantes:
~ La tension par phase tant seulement iI.../3, ou 58 % de celle entre les lignes, on
peut rduire l'paisseur de l'isolant dans les encoches. Cela permet de grossir
la section des conducteurs et, par consquent, la puissance de la machine.
23
~ Lorsque l'alternateur est en charge, la tension induite par phase se dforme de
sorte que la forme d'onde n'est plus tout fait sinusodale. Cette distorsion est
due principalement la prsence des tensions de troisime harmonique qui se
superposent la tension fondamentale. Avec une connexion en toile, les
troisimes harmoniques n'apparaissent pas entre les fils de ligne, car elles
s'annulent. Par consquent, si l'on utilise une connexion en triangle, ces
tensions s'additionnent et provoquent la circulation d'un courant dans le
triangle et, par consquent, occasionnent des pertes joules supplmentaires.
Rotor: Le rotor est l'inducteur :
Il peut tre un aimant permanent, la rgulation de la tension de sortie se fera en rgulant la
vitesse de rotation de l'alternateur, (la frquence du courant variera galement). C'est le
principe de la dynamo de vlo, qui est .en fait un petit alternateur.
Plus couramment un lectroaimant assure l'induction. Ce bobinage est aliment en courant
continu, soit l'aide d'un collecteur bague rotatif (un double bague avec balais) amenant
une source extrieure,soit par un excitateur diodes tournantes et sans balais. Un systme de
rgulation permet l'ajustement de la tension et de la phase du courant produit.
Si l'on fait tourner le rotor, les lignes de flux produites par les ples inducteurs balaient les
trois enroulements du stator et induisent dans ceux-ci des tensions triphases. Les rotors sont
ples saillants ou ples lisses selon qu'ils sont entrans basse vitesse par des turbines
hydrauliques ou haute vitesse par des turbines fi vapeur.
Rotor ples saillants: Comme les alternateurs sont raccords directement aux turbines et
puisqu'une frquence de 50 Hz (ou 60 Hz) est impose, on doit placer un grand nombre de
ples sur le rotor. Les alternateurs basse vitesse ont donc toujours un grand diamtre, de
faon donner l'espace ncessaire pour loger le nombre de ples requis.
Rotor ples lisses: Les turbines vapeur tournent toujours haute vitesse: 3600 r/min ou
1800 r/min sur les rseaux de 60 Hz, et 3000 ou 1500 r/min sur ceux 50 Hz. Par consquent,
les alternateurs qu'elles entrainent comportent 2 ou 4 ples seulement, on les dsigne sous le
nom de turbo-alternateurs. Leur rotor une forme cylindrique, car les ples sont forms en
plaant des bobines dans les encoches tailles dans la masse du rotor. Les forces centrifuges
intenses dues la haute vitesse de rotation imposent une limite au diamtre du rotor; comme,
d'autre part, les grandes puissances ncessitent un gros rotor, on est oblig de lui donner une
forme allonge.
24
Excitatrice: L'excitatrice d'un alternateur constitue un de ses lments les plus importants.
En effet, le champ doit non seulement induire une tension approprie, mais il doit aussi
pouvoir varier rapidement lorsque la charge varie brusquement.
Excitation sans balais: cause de l'usure des balais et de la poussire conductrice qu'ils
dgagent, il faut effectuer une maintenance constante des bagues et du collecteur, sinon on
risque des courts circuits. Pour viter ce problme, on utilise de nos jours les systmes
d'excitation sans balais dans lesquels un alternateur-excitatrice et un groupe de redresseurs
fournissent le courant continu l'alternateur principal.
Si on compare avec le systme d'excitation avec balais, on constate qu'ils sont identiques,
sauf que le redresseur remplace le collecteur, les bagues et les balais. Le courant de
commande Ic provenant de l'excitatrice pilote rgularise lx, comme dans le cas d'une
excitatrice courant conventionnelle.
L'alternateur-excitatrice et les redresseurs sont monts en bout d'arbre et tournent ensemble
avec l' alternateur principal.
Choix de l'alternateur: Nous allons choisir un alternateur avec un rendement de 95 %. Il
sera inducteur t o u r n ~ t qui tourne une vitesse de 3.000 tr/min sur une frquence de 50 Hz,
ce qui impose le nombre de ple p = 1. On prfre la connexion en toile d'aprs les raisons
cites plus haut avec un rotor ple lisse qui est toujours utilis pour les turbines vapeur ou
gaz.
11-2-2-2 : La turbine
a : Dfinition
Une turbine vapeur est constitue d'un grand nombre de roues (une centaine pour un modle
de puissance) portant des ailettes. La vapeur sous pression traverse d'abord les roues de petit
diamtre avant d'atteindre les roues de plus grand diamtre. La turbine tourne alors en
entranant l'alternateur qui lui est accoupl.
Une turbine est essentiellement une roue destine transformer le mouvement d'un corps
fluide (eau ou gaz) en un mouvement de rotation. Les centrales lectriques utilisent
principalement deux types de turbines couples un alternateur : les turbines eau et
vapeur.
2S
Dans notre cas nous utiliserons la turbine vapeur.
E.S.P. This (2006-2007)
b : Lefonctionnement
Bien que les turbines vapeur soient construites selon deux principes diffrents ( action ou
raction: voir ci-aprs), leurs lments essentiels sont similaires. Elles se composent de
tuyres ou de jets, et d'ailettes. La vapeur s'coule dans les tuyres, dans lesquelles elle se
dilate. Ainsi, sa temprature diminue et son nergie cintique augmente. La vapeur en
mouvement exerce une pression contre les ailettes, entranant leur rotation. La disposition des
jets et des ailettes, fixes ou stationnaires, dpend du type de turbine. la sortie du dernier
condenseur (changeur thermique), l'eau peut tre de nouveau vaporise et surchauffe. L'eau
ou la vapeur en sortie est alors ramene vers la chaudire et la pompe alimentaire , qui
compresse de l'eau l'tat liquide. Il s'agit d'une turbine auxiliaire intgre au cycle
thermodynamique de la turbine principale utilisant de la vapeur soutire dans celle-ci.
Les turbines vapeur sont des machines relativement simples dont la seule partie mobile
importante est le rotor. Elles possdent toutefois un quipement annexe, ncessaire leur
fonctionnement. Parmi celui-ci, un palier de tourillon supporte l'arbre et un palier de bute le
positionne de manire axiale. Un systme d'huile assure le graissage des paliers; des joints
rduisent les pertes de vapeur tout au long de son trajet. Enfin, un systme d'tanchit
empche la vapeur de s'chapper l'extrieur de la turbine et l'air d'y entrer. La vitesse de
rotation est commande par des soupapes situes aux entres d'admission de la machine et
pilotes par des systmes de rgulation lectroniques ou mcaniques. Les turbines raction
dveloppent une pousse axiale considrable, du fait de la chute de pression sur les ailettes
mobiles. Cette pousse est gnralement compense par l'utilisation d'un piston d'quilibrage
La turbine vapeur utilise des principes thermodynamiques. Lorsque la vapeur se dilate, sa
temprature et donc son nergie interne diminuent. Cette rduction de l'nergie interne
s'accompagne d'une augmentation de l'nergie cintique sous la :forme d'une acclration des
particules de vapeur. Cette transformation rend une grande partie de l'nergie disponible.
Ainsi, une rduction de 100 kJ de l'nergie interne, du fait de la dilatation, peut provoquer un
accroissement de la vitesse des particules de vapeur de l'ordre de 2800 km/h, de telles
vitesses, l'nergie' disponible est importante. Lorsque la pression de la vapeur d'eau en sortie
de la turbine est gale la pression atmosphrique, la turbine est dite condensation.
Aujourd'hui, les turbines vapeur sont gnralement limites une temprature maximale de
580 "C dans le premier tage, et une pression maximale d'admission de 170 180 bars.
26
P.F.E. Prsentpar M. Mbaye SECKet M. MoussaDIOP E.S.P. This(2006-2007)
C : Les type de turbines
~ Turbines action
La forme la plus simple de turbine vapeur est la turbine action, dans laquelle les jets sont
fixs sur la partie intrieure de l'enveloppe de la turbine, et les ailettes places sur le bord des
roues tournantes montes sur un arbre central. La vapeur se dplaant dans une tuyre fixe
passe sur les ailettes incurves, qui absorbent une partie de l'nergie cintique de la vapeur
dilate, faisant ainsi tourner la roue et l'arbre sur lesquels elles sont montes. Cette turbine est
conue de manire ce que la vapeur entrant par une extrmit de la turbine se dilate travers
une succession de tuyres jusqu' ce qu'elle ait perdu la majeure partie de son nergie interne.
~ Turbines raction
Dans la turbine .raction, une partie de l'nergie mcanique est obtenue par l'impact de la
vapeur sur les ailettes. La partie la plus importante est obtenue par l'acclration de la vapeur
lors de son passage dans la roue de la turbine, o elle se dilate. Une turbine de ce type se
compose de deux jeux d'ailettes, l'un fixe, l'autre mobile. Ces ailettes sont disposes de telle
faon que chaque paire joue le rle de tuyre, travers laquelle la vapeur se dilate lors de son
passage. Dans chaque tage, une faible quantit d'nergie thermique est convertie en nergie
cintique. La vapeur se dtend dans les aubes fixes, puis entrane, les aubes mobiles disposes
sur la roue ou le tambour de la turbine. Les ailettes d'une turbine raction sont en gnral
montes sur un tambour, qui fait alors office d'arbre.
Les turbines raction ncessitent en gnral davantage d'tages que les turbines action. Il a
pu tre dmontr que, pour le mme diamtre et la mme gamme nergtique, une turbine
raction a besoin de deux fois plus d'tages pour obtenir un rendement maximal. Les grosses
turbines, qui sont gnralement action, utilisent une certaine raction la base du trajet de
vapeur pour assurer un dbit efficace travers les auges. Nombre de turbines, qui sont
normalement raction, disposent d'un premier tage de commande d'impulsion, qui permet
d'envisager la rduction du nombre total d'tages ncessaires. Les arbres des turbines de
chaque tage sont relis entre eux au moyen d'accouplements.
Choix: D'aprs l'tude de toutes ces informations, nous avons port notre choix sur une
turbine axiale action. Ainsi on n'aura pas prvoir l'utilisation d'un piston d'quilibrage. La
turbine action cote aussi moins chre.
27
d : gnrateur de vapeur
C'est lui qui fait la transmission de la chaleur du liquide caloporteur vers l'eau distille
l'entre pour transformer celle-ci en vapeur la sortie. C'est un gros changeur divis en sous
changeurs en sries qui sont: l'conomiseur, l'vaporateur et le surchauffeur.
11-2-3 : Lesystme primaire
Ce systme est compos essentiellement du concentrateur, de l'absorbeur ou rcepteur et du
fluide caloporteur et le stockage du fluide.
11-2-3-1 : Le' concentrateur
Un capteur concentration est un capteur solaire comportant un systme optique (rflecteur,
lentilles, ...) destin concentrer sur l'absorbeur le rayonnement reu (d'o le nom de capteur
concentration). Naturellement ces concentrateurs doivent suivre le mouvement du soleil. En
utilisant des rflecteurs afin de concentrer les rayons du soleil sur l'absorbeur d'un capteur
solaire, cela permet de diminuer grandement la taille de l'absorbeur, ce qui rduit les pertes de
chaleur et a u g m e ~ t e son efficacit hautes tempratures. Un autre avantage de ce systme est
que les rflecteurs sont sensiblement moins coteux, par unit de surface, que les capteurs.
Ce qui nous amne porter notre choix sur les rflecteurs.
Parmi les rflecteurs on peut citer trois types qui sont les plus utiliss:
- Rflecteur cylindro-parabolique
- Rflecteur parabolique
- hliostats
Les rflecteurs permettent d'obtenir une temprature leve au foyer, cependant ils n'utilisent
que les rayons directs du soleil et il faut que les capteurs soient orients en permanence vers le
soleil, on utilise en gnral pour cela des systmes asservis [A] [B] [c].
11-2-3-2 : Le rflecteur cylindro-parabolique
Un capteur cylindro-parabolique est un capteur concentration foyer linaire utilisant un
rflecteur cylindrique de section parabolique. Dans un concentrateur cylindro-parabolique, le
fluide caloporteur (eau, huile thermique ou gaz) peut tre port environ 400 oc.
28
F.S.r (2006-200'7)
;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.; / /////// x..xcccccccccccccoxcccccc :.:.:.:.:.:.:.;.:.. ..xccx ccc ;.:.:.:.:.y ..:.:,o:.:,o:,o:..

Figure 5 : Systemes il rflecteur cylmdro-perebolique
: Le rfle...:ftln' panmHqaw
Un c,qkur pawbn1ique est un capteur concentration utilisanr un rtkekur :;;,1 forme de
oarabolc ,.cl" rvolution -': oui ".on"('j1'l"O j,x:;; ravons solaires dans un fover nonctuel 1")""'11<: ;._'1,..- ...... _i.... :,,;....... _...... .1 ............ _ -1 ... .... .... .... ..... .....,,'.. 1.... .-.,"" ... ...oh.-..>_ '--.,1: .. -... _...... ... ...... , oJ ...... oJ
concenrraieurs paraboliques, on peut obtenir des temperatures leves (jusqu' IS00 cC).
Un capteur parabolique est un capteur concentration utiU:imlt un rflecteur en forme de
centrale tour possde un ensemble de grands miroirs qui suit le mouvement du soleil. Les
miroirs concentrent les rayons du soleil sur le rcepteur en haut de la grande tourv Un
ordinateur ganb les miroirs aligns aD que les rayons du soleil. qui sont rflchis, visent
29
thermodvuamiouc
...
solaire doit remplir nimes
'/ La concentration du rayonnement sur !' entre du recepteur.
.../ son absorption sur les parois du rcepteur et 1a transformation de son nergie
if sa dlivrance un cycle thermodynamique associ un allerTllteUr pour la
p.),:,. transformer Ie solaire i'11 (;!':li('l""l'. 1'1 suffit ")1 nrincinc de 1" "P'-
" v.: e..- .:.;.s , v', ,.. (,,;n . . CH _ ..oC s.,_. _ .<-) .L,. (.( .: __ ..._ io V t .._ ""'.i
thermodynamique performant, il faut diminuer fortement la surface de rccption pour
devant k rcepteur une optique qui concentre sur celui-ci le rayonnement capt sur une
surface bien suprieu.n:. (in caractrise la performance du systme par le chittre
" concentration " qui le rapport de la surface de collecte sur la surface du capteur,
:JO
E.S.P. This (2006-2007)
Dans la pratique cette optique est systmatiquement constitue de miroirs. C'est la nature
gomtrique des surfaces rflchissantes mises en uvre qui va dfinir la " filire" utilise.
11-2-4-2: La Filire parabolique
~ Les principes
Cette filire est un peu la rfrence en ce sens qu'elle met en uvre la meilleure surface
rflchissante possible : une parabole de rvolution. Avec ce miroir, tout rayon incident
parallle l'axe optique passe, aprs rflexion, par un mme point. On appelle ce point le
" foyer". Pour fonctionner correctement, un tel miroir doit viser en permanence le soleil (dont
l'image se forme alors au foyer). On y parvient en animant son axe d'une double rotation.
. .
Mais ceci implique, sans doute, une limitation de taille, la surface de collecte ne pouvant
gure dpasser ici quelques centaines de mtres carrs dans des conditions conomiques
raisonnables. C'est l que rside la spcificit de la filire qui n'a-d'intrt que si l'on parvient
utiliser le faible dbit de chaleur rcolt dans une installation compacte, mobile avec le
miroir (et si possible sans liaison avec l'extrieur autre qu'lectrique) et installe dans une
position voisine du foyer.
L'optique peut ici, thoriquement, dlivrer des concentrations suprieures 10000, ce qui est
norme. Dans la pratique, on ralise ventuellement des miroirs moins performants (en
s'autorisant des marges d'erreur sur la qualit gomtrique des surfaces mises en uvre). Ils
peuvent alors tre conomiques tout en dlivrant des concentrations d'un excellent niveau.
Plusieurs ralisations rcentes de qualit dlivrent ainsi des concentrations de l'ordre de 4000.
~ La technologie
La ralisation de surfaces rflchissantes paraboliques de rvolution pose des problmes
particuliers auxquels il est possible d'apporter des rponses varies. Presque toutes sont
aujourd'hui bases sur la mise en uvre de verre argent en face arrire comme surface
rflchissante. C'est de loin celle qui prsente le meilleur rapport qualit-prix: excellent
coefficient de rflexion, bonne tenue aux intempries, prix modr d'un produit industriel
d'utilisation massive et ancienne.
La mise en uvre de ce matriau en vue de raliser une parabole de rvolution peut tre
effectue selon deux voies principales:
31
E.S.P. This (2006-2007)
..
Juxtaposition de nombreux trapzes plans pour approcher au mieux la surface parabolique
thorique,
Utilisation de verre m ~ n c e cintr en double courbure qui permet, thoriquement, d'obtenir une
parabole parfaite.
La premire solution n'est valable que si elle est plus conomique car elle est obligatoirement
moins performante.
Dans les deux cas, un problme central sera le collage des miroirs lmentaires sur la surface
support qui leur imposera la gomtrie et le positionnement adquat. En effet, il y a de
difficiles problmes de compatibilit chimique entre les colles et les vernis qui protgent
l'argenture des miroirs.
11-2-4-3 : La filire centrale tour (hliostats)
~ Les principes .
Ici, on ambitionne d'emble de dployer des surfaces de miroir trs importantes. On renonce
donc utiliser une vraie parabole (qui devrait tre mobile ce qui est inimaginable en grande
taille) et on la remplace par une pseudo-parabole, foyer fixe, d'une taille pouvant atteindre
plusieurs centaines de milliers de mtres carrs, et constitue de centaines ou de milliers de
miroirs plans orientables. On appelle ces miroirs " hliostats " (en grec : qui fixe le soleil).
Leur ensemble constitue le " champ d'hliostats ". On utilise ici l'image de la pseudo-parabole
parce que, chaque instant, les rayons solaires atteignant le centre de chaque hliostat sont
rflchis par le dispositif en direction d'un point unique, fixe dans le temps, et qui joue donc
le rle du foyer de la l?arabole. La surface rflchissante de cette " parabole" tant fatalement
dploye au sol, son" foyer" se trouve en altitude. Pour y disposer le rcepteur, on est amen
construire une tour de grande hauteur qui donne son nom la filire.
Pour faire fonctionner ce dispositif, il faut raliser un " tir crois " des hliostats et pour cela
animer chaque miroir d'un mouvement de suivi du soleil particulier. Comme dans le cas
prcdent, ce mouvement peut tre obtenu par la combinaison de deux rotations.
Il faut encore avoir optimis la position des miroirs (ou plutt la distribution de leur densit
sur le terrain), la hauteur du foyer (qui conditionne la hauteur de tour et donc son cot) et la
gomtrie de l'ouverture du rcepteur. A travers un calcul difficile, on minimise le rapport
32
P.f.E. Prsent par M. Mbaye SECK et M. Moussa DIOP E.S.P. This (2006 -2007)
cot d'investissement/rendement thermo-optique de l'ensemble champ d'hliostats-rcepteur.
Les optima trouvs, en particulier en ce qui concerne l'altitude du foyer, se traduisent par des
courbes trs plates qui laissent pas mal de libert au concepteur (autrement dit, il peut faire
varier dans une large fourchette la hauteur de tour sans dgrader gravement le rendement).
Si ce travail a t bien men, et si l'ouverture du rcepteur est situe dans le plan focal de
l'optique, on peut obtenir par ce moyen des concentrations de l'ordre de 700 (cas de Thmis
par exemple 42,5 d latitude).
~ La technologie
On va donc mettre en uvre ici des hliostats qui sont des miroirs plans anims d'un
mouvement continu obtenu partir d'une double rotation.
La premire question qui se pose lors de leur conception est celle de leur taille. Sur le plan
optique, on peut dmontrer que les performances d'une pseudo-parabole s'amliorent avec la
finesse de son dcoupage en miroirs lmentaires. On a donc de ce point de vue intrt
construire des hliostats de petite taille. Du point de vue structurel, les petits hliostats sont
galement avantageux dans la mesure o les efforts qu'ils supportent (et qui sont
essentiellement ceux dus au vent) sont proportionnels au cube des dimensions linaires de
l'appareil: cela signifie qu'ils augmentent proportionnellement plus vite que la surface
dveloppe, et donc que la rigidit de la structure (pour une tolrance donne sur sa
dformation relative) doit augmenter avec la taille. Autrement dit, les petits hliostats peuvent
tre construits plus" souples" que les grands et mettre en uvre moins de matriau au mtre
carr. De ce point de vue, ils devraient tre plus conomiques.
11-2-4-4 : Filire cylindro-parabolique
~ Les principes
On va donc mettre en uvre ici des miroirs cylindro-paraboliques d'orientation est-ouest
(dans ce cas le mouvement de suivi du soleil se limite une rotation si lente qu'elle peut tre
assure sans automatisme), ou nord-sud, ce qui suppose toujours une rotation unique, mais
plus grande vitesse et qui doit donc tre automatise.
33
Les avantages recherches portent surtout sur Ia simplification de la motorisation el de la
c:onuniln,.:; du mouvement En substance, le de deux axes un seul,
ralisations de Luz les plus rcentes). En second lieu, Ll surface des miroirs ,;;:.<1 plus incline
par rapport au vecteur soleil quand celui-ci ,;;:<1. bas sur lhorizon, ct ce qu'on appelle le
" rendement cosinus " du champ de miroirs en est dgrad.
Le nnrou cylindre-parabolique est constitu en gena! de p<'lrrneaux en verre sans fer
Celle: Cl est supnorte par des D":' lnes en charpente mtallique fondes au sol 1/\1 rBI rel.
J .... ... t . _. .. t }
Iacilir de j'orlf'ntat1.un parmi d'autres avantages, nous permet de choisir les cylindre-
1
P.F.E. Prsent par M. Mbaye SECKet M. MoussaDIOP E.S.P. This (2006-2007)
paraboliques comme rflecteurs pour la centrale. Il faut noter que le paramtre cot
d'investissement a beaucoup influ sur ce choix et il est conseill d'utiliser ces rflecteurs
pour les centrales de petites tailles.
a : Rception du rayonnement.
Il faut maintenant absorber le rayonnement concentr sur une surface noire (ou dans un
volume semi-transparent, noir dans son paisseur), et de ce seul fait absorbante pour le spectre
visible (et donc pour le spectre solaire qui en est nergtiquement trs proche), et transfrer la
chaleur ainsi gnre un fluide caloporteur.
Le rcepteur est ici compltement diffrent des autres beaucoup de points de vue et d'abord
parce qu'il est linaire au lieu de surfacique. Ensuite parce que, dispos au foyer d'un systme
optique peu performant, on le fait travailler, comme nous l'avons dj signal, temprature
modre (infrieure 400C), ce qui permet d'utiliser comme fluide caloporteur une huile de
synthse. Enfin parce que, pour la mme raison, il est protg de pertes thermiques par deux
dispositions originales' et efficaces:
./ le placement de l'absorbeur dans un tube de verre sans fer vid de son air, ce
qui supprime radicalement les pertes convectives et limite les pertes radiatives.
./ le revtement du tube chaudire d'une couche de cermet (il s'agit d'un
composite cramique/mtal projet par plasma) qui joue le rle d'un absorbeur
slectif (96% d'absorptivit pour 20% d'missivit la temprature de travail)
bien plus performant que la peinture utilise dans les cas prcdents ; cette
disposition est rendue possible ici par la mise sous vide de l'environnement du
tube chaudire (le cermet ne rsisterait pas la prsence d'oxygne).
La mise au point de l'enveloppe vide est particulirement remarquable en ce qu'elle suppose
la mise en uvre de quantits impressionnantes de joints d'tanchit et de joints de dilatation
(soufflets inox) qui tous assurent leurs fonctions dans un environnement mtorologique non
protg.
b : L'absorbeur Afin
d'augmenter encore la temprature du liquide caloporteur, on utilise des capteurs isols avec
vitrage. Pour cela, l'absorbeur doit tre enferm dans un tube de verre vid de son air (tout en
35
E.S.P. This (2006 -2007)
laissant le rayon solaire entrer, la vitre empche les rayons infrarouges diffuss de sortir: on
ralise ainsi l'effet de serre ). Le verre, fragile, peut tre remplac par certains matriaux
plastiques ayant la proprit de raliser l'effet de serre (c'est le cas de matriaux comme le
polycarbonate, le mthacrylate et le tedlar). L'inconvnient principal des matriaux plastiques
est la dgradation 'de certaines de leurs proprits, avec le temps pour certains, ou le prix lev
pour d'autres.
Le foyer est un tube de 2 cm de diamtre, de mme longueur que le rflecteur. Il faudrait
utiliser un tube de cuivre (bon coefficient de conduction)
L'nergie est ensuite cde au fluide caloporteur et transfre vers un rservoir de stockage
d'nergie [5].
c : Le rflecteur
Le capteur doit permettre de transmettre le plus de chaleur possible au fluide. Pour cela, il faut
que la plaque soit'mtallique car seuls les mtaux ont de bons coefficient de conduction de la
chaleur (cuivre: 300, aluminium: 200, acier: 60). En effet, plus le coefficient de conduction
est faible et plus la plaque doit tre paisse, ce qui augmente le prix.
Le plus efficace serait une plaque d'inox avec une face polie miroir et protge par un film
plastique [5].
d : La concentration effective
La thorie dveloppe suppose que le miroir est parfait et qu'il n'y a pas d'obstacle entre le
soleil et le miroir. En ralit, de nombreux facteurs diminuent l'nergie disponible sur le plan
focal. Les pertes sont dues la rflexion, la convection, l'absorption, aux irrgularits
gomtriques de la surface du rflecteur, aux imperfections de rglage, etc. tous ces petits
dtails prennent une importance considrable si l'on veut une concentration leve.
Le rapport entre l'nergie effectivement reue et celle qui est calcule partir du facteur de
concentration est not F. la concentration effective est CE = F C. dans les meilleurs conditions,
on ne dpasse pas F ='0,6 [3].
e : Le fluide caloporteur
C'est un fluide en mouvement qui reoit de la chaleur en un point de son circuit et en cde en
un autre point. Le rle jou par les fluides caloporteurs est fondamental.
36
P.F.E. Prsentpar M. Mbaye SECKet M. MoussaDIOP E.S.P. This (2006-2007)
Le choix de ces fluides doit tre effectu en tenant compte de certains critres:
Soit un fluide froid qui se rchauffe de Tl T2 par change de chaleur avec un fluide chaud
la temprature T2 constante. Si l'changeur tait de trs grande dimension, on aurait T'2 ~ T2
et en supposant qu'il n'ya pas de pertes de chaleur, le flux de chaleur chang pour valeur :
Q = mC (T2 - Tl)
m: dbit de fluide caloporteur
C : chaleur massique a pression constante de ce fluide.
Pour un dbit donn m et une diffrence de temprature impose T2 - Tl ; on voit que Q
=f(C)
Donc un fluide caloporteur doit avoir une chaleur massique leve.
L'quation de dbit relative au fluide caloporteur est: m = pYS
. .
S : section de passage du fluide,
V : vitesse de circulation du fluide,
p : masse volumique du fluide pour les conditions temprature etde pression considres.
Onadonc S =mNp.
Donc pour un dbit et une vitesse de circulation impose on a S = f (l/p).
Donc pour diminuer les dimensions des installations, on a intrt prendre un liquide comme
fluide caloporteurde prfrence un gaz.
En rsum un fluide caloporteur doit avoir une chaleur massique et une masse volumique
forte. Cependant, d'autres lments de choix interviennent:
le cot du fluide,
les phnomnes de corrosion,
les risques de conglation et d'bullition [6].
Choix: Pour toutes ces raisons on peut choisir un sel nomm HITEC comme fluide
caloporteur ; ses proprits sont trs adaptes et il est trs utilis dans beaucoup de centrales
solaires travers le monde. C'est un mlange de sels (40% NaNO, 7% NaN03 et 53% KN03)
dont le point de fusion est voisin de 150 "C.
d : Le stockage du fluide caloporteur
Il nous permet de disposer d'une nergie suffisante pour le dmarrage de la turbine chaque
fois qu'il est ncessaire et une forte autonomie l'absence du soleil (rservoir chaud). Ces
37
P.F.E. Prsent par M. Mbaye SECK et M. Moussa DIOP E.S.P. This (2006-2007)
rservoirs sont surtout des tampons pour le fluide caloporteur; on peut contrler la
temprature du fluide a l'entre de l'absorbeur via le rservoir du fluide froid.
Les rservoirs doivent pouvoir garder les tempratures des fluides respectifs constantes
pendant un temps acceptable.
Pour ces raisons leurs parois seront constitues de l'intrieur vers l'extrieur d'un Acier
nuanc A, d'un calorifuge (laine de roche C) et d'un feuillard galvanis G [7].
!
Chaleur spcifique Conductivit thermique Masse
1 Libelle Cp (j/kg, OC) (W/m. OK) volumique
p (Kg/m
3
)
2082-0,72835T
RITEC 1561,2 0,5
(TOC)
nuance acier A 475 45 7850
Calorifuge : Laine
837 0,0372
de roche
Feuillard galvanis 61 43
Tableau 3 : Caractristiques du HITEC et son stockage [A] [B]
113 : Conclusion partielle
Le fonctionnement de la turbine risque de poser des problmes si le constructeur de la
turbine ne prend pas en compte le fait que la turbine s'arrte chaque soir aprs disparition du
soleil pour redmarrer le lendemain. C'est un fonctionnement pas trs habituel pour les
turbines, cette information va tre donne lors de la commande.
Un problme d'obtention de miroir peut se poser si on veut fabriquer les miroirs sur place
puisqu'il faut d'abord une main d'uvre qualifie difficile trouver car ces miroirs sont
presque inconnus du grand public. En plus la mise d'une couche d'argent sur la face
rflchissante du miroir ncessite une trs grande prcision.
38
P.F.E. Prsentpar M. Mbaye SECK et M. Moussa DIOP E.S.P. This (2006 -2007)
CHAPITRE III: MODELISATION ET DIMENSIONNEMENT DU
SYSTEME
1111 : Modlisation du systme primaire
111-1-1 : Introduction
Modliser le systme primaire consiste tablir successivement les modles de
l'ensoleillement solaire, du miroir cylindro-parabolique, du rcepteur solaire et du stockage.
Stockage frd Stockage chd
La, figure 8 : Schma simplifi du systme primaire.
111-1-2: Modle de l'ensoleillement
Nous prsentons ici le modle du flux solaire qui arrive effectivement la surface du miroir
par un ciel sans nuage. Ce flux varie en fonction de la direction du soleil et de la couche
d'atmosphre traverse. Les paramtres sont donc d'une part, la latitude du lieu, la date. et
l'heure, d'autre p ~ r t , un effet global d'attnuation atmosphrique. C'est ce dernier qui est le
plus difficile cerner. Aussi proposons-nous une modlisation dans les conditions moyennes
sngalaise d'ensoleillement. Notre modle permet de calculer le rayonnement solaire direct
diffrentes dates et heures sur le site.
39
P.F.E. Prsent par M. Mbaye SECK et M. Moussa DIOP E.S.P. This(2006 -2007)
1 .
Ill-1-2-1: Rappel des donnes astronomiques
Dclinaison solaire 0 : C'est l'angle que fait la direction du soleil avec le plan quatorial.
Elle varie au cours de l'anne, passant par un maximum de 2327' au solstice d't, un
minimum de -2327' au solstice d'hiver et devenant nulle aux quinoxes.
Au jour J de l'anne, compt depuis le premier janvier, cette dclinaison peut tre calcule
l'aide de la formule de type sinusodale suivante :
8 = 23,45 x sin : : ~ (J +284) (3.1.2.1.1a)
O est exprim en degrs.
Dans la pratique, on utilise la formule
stn = 0,4 x sin 360 U- 80) (3.1.2.1.1b)
365
Angle horaire du soleil ID : Il est dlimit par le mridien cleste plan et le cercle horaire de
l'astre plan: Les conventions habituellement retenues lui donnent la valeur 0 lorsque le soleil
passe au plan mridien d'un lieu: il est alors midi vrai au soleil. L'intervalle de temps entre
deux passages successifs du soleil au mridien tant de 24 heures, l'angle horaire co au temps
solaire vrai t est donc :
co = 15 x (t -12) en degrs.
(3.1.2.1.2)
(3.1.2.1.3)
. sinw
sma = cos x -.
cash
Hauteur angulaire du soleil (h) et Azimut du soleil (a) : Si on connait la latitude 0 du lieu
d'observation, les relations de trigonomtrie sphrique permettent de relier les coordonnes
horizontales et horaires par les formules suivantes:
sin h = sin0 x sin +cos0 x cos x cosw
Position du soleil: Soit Sle vecteur reprsentant la position du soleil au jour J l'heure t et
la latitude 0. Ses coordonnes sont les suivantes:
Sx =sin0 x cos x cosw - cos0 x sin =cosh x cosa
Sy = -cos x sinw =cos h x sina
Sz = -(cos0 x cos x cos+sin0 x sin ) =sinh
(3.1.2.1.4)
(3.1.2.1.5)
(3.1.2.1.6)
Ill1-2-2 : Flux incident du rayonnement solaire
L'intensit du rayonnement solaire direct au sol pour la position Sdu soleil est dtermine
par la formule suivante [3]:
40
P.F.E. Prsent par M. Mbaye SECKet M. MoussaDIOP - E.S.P. This (2006-2007)
E = s, X EXP( -.1 h )
Axsm
en w/m
2
(3.1.2.2.1)
Eo : constante solaire
A: constante traduisant l'influence de l'paisseur de l'atmosphre et dont les valeurs
courantes se situent entre 2 et 7.
Compte tenu des conditions climatiques moyennes au Sngal, de l'influence de divers
facteurs, d'un coefficient de trouble moyen gale 0.1 et d'une paisseur d'eau condensable
de 1 4 cm; on peut crire les quations suivantes :
Emoyen = 1230 x EXP ( .-; ): dans les conditions normales. (3.1.2.2.2)
4xsm h+2
Emax = 1210 x EXP ( .-; ): par ciel trs pur. (3.1.2.2.3)
xstn h+1
Emin = 1260 x EXP ( ~ 1 ( ): en zone industrielle. (3.1.2.2.4)
2,3xsm h+3
O E est exprim en w/m
2
et h en degrs.
Compte tenu de l'altitude et des conditions atmosphriques au Sngal, on peut retenir la
formule suivante:
Emax = 1210 X EXP ( '-;h )
xsm +1
111-1-3: Modle du miroir cylindro-parabolique
(3.1.2.2.5)
111-1-3-1: Calcul de la concentration produite par le miroir cylinro-
parabolique
ai)
P'
y
Figure 9 : calcul de concentration [3]
41
Considerons la parabole (P), de foyer F, de sommet S et soit un pinceau plat de rayons
lumineux issus du soleil, d'ouverture e convergeant au point I. Apres reflexion, il coupe la
trace (n) du plan horizontal passant par F suivant un sommet AB. Un pinceau plat, de mme
ouverture tombant en S coupe (x) suivant le segment A'B'. Si l'on choisit F comme origine
des coordonnes et des axes Fx et Fy comme indiqu sur la figure 9, on trouve, en appelant p
et <1> les coordonnes polaires du point 1et en admettant avec une erreur du second ordre que F
est le milieu de A ~ , nous aurons:
..
Soit l' et F' les projections orthogonales de B respectivement sur AI et FI. On peut crire:
I
BI' ~ ~ p
BF' ~ -p
2
Nous avons:
tt F'B F'B
Sm(- - <1 = -'=2x (-)
2 FB AB
Sin ( ~ - <1 = Cos <1>
2
- E
F'B .: -p
Cos <p = 2x (-) = 2x (_2)
AB AB
Puis:
AB=2L
Cos li>
(3.1.3.1.1a)
(3.1.3.1.1b)
(3.1.3.1.2)
f=A'S=B'S; A'B' <f e
La fraction de l'energie reue sur le segment focal AB est:
A'B' fCos li>
~ $ = - = - -
AB E
(3.1.3.1.3)
en tenant compte de l' equation de la polaire:
il vient:
2xf
P =1+Cosli>'
42
;: _ Coslj>(1+Cos lj
':l<lJ-
2
--_._- -_ .. _------
E.S.P. This (2006 -2007)
(3.1.3.1.4)
Si 1 est la puissance reue du soleil par unit d'air terreste et H la longueur du miroir,
l'energie correspondant au diametre d'ouverture D du miroir est:
= 21f,e H 2 F X Coslj> (1+Coslj d<t>
o 1+COS<j> 2
f,e
=21 fH 0 Cos<t>d<t>
ID = 21 HfSine
Cette nergie est reue sur un rectangle d'aire Hf e.
La puissance par unit' d'aire est donc:
ID = 2 Hf Sine 1
H CE: H CE:.
Ce qui entrane :
ID _ 2 Sine t
--
H ft E
Et le facteur de concentration est:
C= 2 xstne.
E
(3.1.3.1.5a)
(3.1.3.1.5b)
(3.1.3.1.5c)
(3.1.3.1.4)
(3.1.3.1.5)
111-1-3-2: Calcul des pertes nergtiques dues la gomtrie de la parabole
On peut appeler efficacit du miroir le rapport 1) entre l'nergie dans l'image de gauss et
l'nergie incidente reue par la surface d'ouverture normale aux rayons lumineux:
(3.1.3.2.1)
Si e = 90, 1) = 1/2 = 0,5 alors l'image de Gauss ne recueille que la moiti de l'nergie
incidente
43
111-1-4: Modle de l'absorbeur
E.S.P. This (2006-2007)
On va prsenter deux modles, le premier utilise l'optique gomtrique et l'autre des
quations thermodynamiques.
111-1-4-1 : Premier modle
Dans ce modle on utilise l'optique gomtrique, pour considrer les pertes par rflexion et
dduire ensuite le rayonnement transmis On suppose que le verre n'est pas absorbant et que
le corps noir (mtal absorbant la chaleur) est parfait (il ne rflchit pas la lumire).
a : Les pertes par rflexion
Figure 10 : pertes par rflexion
R nverre -nair
nverre+nair
{
n verre = 1.5: indice de refraction du verre
Avec: nair = 1 : ndice de refraction du verre
R : coefficient de convection verre - air
15-1
R = (-'-) =0.04=4%
1.5+1
On suppose que le verre n'est pas absorbant et que le corps noir est parfait (il ne rflchit pas
la lumire). Si on dsigne par p le coefficient global delccaveetiod, alors on peut crire:
ri.{'lt:I'- \ o ~
pl = RI + R (l-R/ 1+R
3
(1_R)2 1+ R
S
(1-R)21+.......... (3.l.4.la)
p = R + R (1_R)2 + R
3
(l-R)2 + R
S
(l-R)2 + .
En factorisant par R ou voit apparaitre le dveloppement limit de (l-R)/ (1+R)
44
(3.l.4.lb)
(1 ,;r.(:::.-::.:i:.;:j
q! T,<ux".
cc .: : .; =.; t
"', .... 'M';'';'; l ". J 1DI
==-:: ..;t i
::: _-:: ': .; .; =:+...
1a tbenncdynarnique OH a :
0,1.4.2])
4S
(3.1.4.2.3)
Pv Cv et Tl sont la densit, la chaleur spcifique et la temprature du fluide caloporteur
respectivement o Pv Cvdpendent de la temprature Tl'
Par ailleurs, si li est le dbit volumique du fluide caloporteur, on peut crire:
Q1(Z,t)
QI (z, r) = t =Pl Cl VT
I
(z, t)
Aprs l'insertion de l'quation
trouve:
(3.1.4.2.2) et (3.1.4.2.3) dans l'quation (3.1.4.2.1) on
(3.1.4.2.4)
Aprs des simplifications l'quation (3.1.4.2.4) devient:
La condition la limite de l'quation (3.1.4.2.5) est:
TI(O,t)=T e
O: Te est la temprature d'entre du fluide caloporteur.
La condition initiale s'ecrit :
T
amb
est la temprature ambiante.
(3.1.4.2.5)
(3.1.4.2.6)
(3.1.4.2.7)
Par l'analogie avec l'quation de la temprature du tube absorbeur T
z
est obtenue:
a(t:.Qz(z,t)) ) ) ()
at =(qabsorbe(t - qinterieur(Z, t - qgagne Z, t ) ~ z
Avec:
qinterieur le transfert thermique entre le tube en cuivre et l'enveloppe en verre
qabsorbe l'nergie solaire absorbe
De la thermodynamique on a :
Qz(Z, t) =pzczAzzTz(z, t)
46
(3.1.4.2.8)
(3.1.4.2.9)
P.F.E. Prsentpar M. Mbay SECK et M. Moussa DIOP
E.S.P. This (2006 -2007)
Pz , Cz et T
z
sont la densit, la chaleur spcifique et la temprature du tube absorbeur
respectivement. La substitution de l'quation (3.1.4.2.9) dans l'quation (3.1.4.2.8) entraine:
A JT
z
(z,t) ( ) () ()
pzcz Z Jt = qabsorbe t - qinterieur Z, t - qgagne Z, t
La condition initiale de l'quation (3.1.4.2.10) est:
Tz(z,O) = T
amb
Et de mme pour la temprature du tube en verre T
3
, on peut crire:
A JT3(z,t) ( ) ( )
P3
C3
3 Jt = qinterieur Z, t - qexterieur Z, t
(3.1.4.2.10)
(3.1.4.2.11)
(3.1.4.2.12)
P3 ,C3 et T
3
sont la densit, la chaleur spcifique et la temprature d'enveloppe en verre
respectivement.
qexterieur est le transfert de chaleur entre le verre et l'environnement extrieur.
La condition initiale de l'quation (3.1.4.2.12) est:
T
3
(z, 0) = T
amb
(3.1.4.2.13)
, a: Transfert de chaleur entre l'absorbeur et le fluide caloporteur
En considrant la convection force dans le tube du fluide caloporteur. qgagne est calcule
par l'quation de Dittus-Boelter pour l'coulement totalement dvelopp dans un tube
circulaire lisse. Par consquent le nombre de Nusselt local est donne par :
(3.1.4.2.14)
O : n=O.4 pour: T
z
> Tl
n = 0.3 pour
PTlIRel est le nombrede Prandtl et le nombre de Reynolds respectivement.
Donc: qgagne est donne par :
(3.1.4.2.15)
47
E.S.P. This(2006-2007)
Avec: h
v
le coefficient d'change convectif entre le fluide caloporteur et l'absorbeur.
b : Transfert de chaleur entre l'absorbeur et l'enveloppe de verre
ce transfert qinterieure s'effectue par convection et le rayonnement, il est donc obtenue par la
relation:
qinterieure =qinterieure convection + qinterieure rayonnement
L'change par convection est calcul par les relations tablies pour une convection naturelle
entre deux cylindres horizontaux et concentriques.
2Trg
qinterieure convection =In(D3) (T2 - T3 )
.: D2
(3.1.4.2.16)
g est le cofficient de conductivit thermique du vide dans le tube en verre. On peut
assimiler cette convection une conduction puisqu'il n'y a pas d'air ce niveau.
qinterieure rayonnement est donne par la formule suivante:
qinterieure rayonnement
O:
(3.1.4.2.17)
) Eabs est l'missivit de l'absorbeur, elle dpend de la temprature de l'absorbeur
) E
g
est l'missivit de l'enveloppe du verre
) 6 est la constante de Stefan Boltzmann [8].
c : Transfert de chaleur entre l'enveloppe de verre et l'environnement
qextrieur =qextrieur convection + qextrieur rayonnement
Dans le cas des cylindres fins dans l'air, Collis et Williams donnent le nombre de nusselt
Nud [6] :
48
N
ud
= 0.3+ 0.5Re 0.6
E.S.P. This (2006-2007)
Ainsi, nous pouvons calculer le coefficient de convection he sachant que :
Donc:
Avec:
A: coefficient de conductivit thermique du verre = 1,35 w /m.K [8] ;
D : le diamtre du cylindre.
d : Diffrence de temperature entre le fluide et l'absorbeur
Maintenant.que nous connaissons la puissance reue par la parabole, nous admettons que celle
ci est intgralement transmise au tube se trouvant au foyer.
Nous allons maintenant dterminer le gain de temprature pour le fluide caloporteur circulant
dans la conduite.
Pour cela considrons l'quation de la chaleur en rgime transitoire en deux dimensions en
admettant que les propriets physiques C et . sont independants de T [6] [7]:
Avec:
);> a =- est.la diffusivit thermique en m

2/s;
. pc
);> q une source de chaleur en W/m
3
En rgime permanent on peut approximer l'quation par :
49
(3.1.4.2.18)
(3.1.4.2.19)
Puisque z =u. t
E.S.P. This (2006-2007)
Avec:
~ u = uo la vitesse (constante) du fluide.
~ p sa masse volumique.
~ C, sa chaleur massique.
~ A sa conductivit thermique.
P=q/pC
Considrant d'ordre 1 tous les termes dans la direction de l'coulement, et d'ordre e tous les
termes dans la direction normale l'coulement, cette quation devient :
er (02T)
peu oz = oz 2 +p
a
2r
+Q_ Kar =0
az 2 az
Avec:
~ Q =PIA
~ K = peul A
Pour rsoudre l'quation diffrentielle posons :: =u ,
On obtient:
u= ~ +A. exp(Kz)
Q A A
!1T = T(z) - T(O) = -z +-. exp(Kz) +-
K K K
avec:
(3.1.4.2.20)
(3.1.4.2.21)
(3.1.4.2.21)
~ !1T la diffrence de temprature entre l'absorbeur et le fluide caloporteur. En z=L, la
variation de temprature est nulle :
A = Q
K.exp (KL)
50
~ Les pertes de charges
Elles s'crivent :
H = f L
2
D2g
Avec:
~ H: les pertes de charges ; m
~ f: coefficient des pertes de charges ;
~ L: la longueur de la conduite ; m
~ D: le diamtre hydraulique ; mm
~ u: la vitesse du fluide. mis
E.S.P. This (2006-2007)
Pour faire augmenter ou diminuer cet cart de temprature, nous pouvons jouer sur plusieurs
paramtres :
Le dbit: plus il sera grand, plus l'cart de temprature sera faible (et les pertes de
charges augmenteront).
Le diamtre de notre tuyau: plus il sera faible plus l'cart de temprature sera grand
(mais les pertes de charges aussi).
On peut galement tenter de construire un tuyau en serpentin afin que le fluide
caloporteur passe plusieurs fois dans la ligne focale [A] [B] [C].
111-1-5: Modle du dbit du HITEC
Dans un premier temps, le dbit du HITEC l'entre du rcepteur ou absorbeur solaire est
donn par la relation:
Puissance totale incidente-pertes totales
Qa = (3.1.5.1)
CPX(TH5N-THEN)
Avec:
~ CP: capacit calorifique;
~ THSN: temprature nominale de sortie du HITEC ;
~ THEN: temprature nominale d'entre du HITEC.
51
(3.1.5.3)
(3.1.5.2)
Nous allons ensuite procder par itration pour trouver le dbit, les pertes et la temprature
la sortie du rcepteur solaire. Le test de convergence est fait sur la temprature de sortie du
HITEC ; il doit vrifier la condition suivante
(TH5N - nominale de sortie) 2 < 1C
Finalement, le dbit Qa s'crit:
Q
QtxCPx(TH-temperature d'arrive du HITEC)
a =
CPx(TH5N-THEN)
Apres tre sortie du tube absorbeur une grande partie du HITEC est achemine vers le
gnrateur de vapeur; l'autre est soutir pour maintenir constante la temprature du HITEC
froid .
Soit: Q1 et Q2 les dbits achemins respectivement vers le gnrateur de vapeur (1l<1) et
vers le ballon froid , on peut crire:
Q1 = Il x Qa
Q2 = Qa - Q1
111-1-6: Modle de stockage
);> Hypothses
La temprature du HITEC est homogne dans le rservoir de stockage.
La variation d'nergie interne du HITEC est gale aux pertes thermiques totales du stockage.
Les pertes thermiques du rservoir sont celles d'une cuve cylindrique isotherme.
);> Description du modle
Les deux rservoirs de stockage ont un volume identique de x m
3
chacun. Le stockage
chaud et le stockage froid contiennent du sel dont la temprature est fixer.
Nous allons prsenter un modle nergtique de ce systme de stockage, obtenu en crivant
l'quation diffrentielle du bilan nergtique et la loi de conservation de la masse. Il nous
permettra d'tudier la gestion de la centrale pour un stockage thermique et une prvision
mtorologique donns. Le modle des pertes thermiques permet d'tudier le refroidissement
du stockage aprs arrt de la centrale.
111-1-6-1 : Modle du stockage chaud
Le stockage chaud est caractris par la masse MH et la temprature TH de stockage du
liquide HITEC.
52
P.F.E. Prsent par M. MbayeSECKet M. MoussaDIOP E.S.P. This(2006-2007)
Le systme considr est le rservoir de stockage de liquide chaud.
Selon la loi de conservation de la masse, nous pouvons crire l'quation diffrentielle
suivante:
:t(MH)= QHE - QHS
Soit:
MH =MHO + (QHE - QHS) x t
(3.1.6.1.1)
(3.1.6.1.2)
(3.1.6.1.3)
(3.1.6.1.10)
Avec:
);> MHO la masse du HITEC t = 0, t est le temps coul;
);> QHE est le dbit massique entrant dans le rservoir et QHS celui sortant.
Le bilan nergtique instantan dans le rservoir s'crit sous la forme :
d '.
dt (MH X CP x TH) = QHE x CP x THE - QHS x CP x TH - Wp
Avec:
);> THE et TH sont respectivement les tempratures l'entre et la sortie ou
l'intrieur du systme.
);> Wp reprsente les pertes nergtiques.
);> CP est la capacit calorifique du liquide et est suppose constante durant tout le
processus. '
Par suite on obtient:
.!!:..(MH x TH) = QHE x THE - QHS x TH - Wp (3.1.6.1.4)
dt CP
d d Wp
MH x dtTH +TH x.dtMH = QHE x THE - QHS x TH - CP
En remplaant .!!:.. MH par son expression dans l'quation (3.1.6.1.1), nous obtenons
dt
MH x .!!:.. TH = QHE x THE - QHS x TH - Wp - TH x (QHE - QHS)
~
MH x .!!:..TH =QHE x THE - QHE x TH _ Wp
dt CP
MH x .!!:..TH = -QHE x (TH - (THE _ Wp '))
dt QHEXCP
!:..-TH= _ QHE>< (TH _ (THE _ Wp )}
dt MH QHEXCP
Soit
TH=A +(THO-A) X e(BXt)
Avec:
53
A =THE- W"
QHBXCP
B = _QHE
MH
THO =temprature du stockage "chaud" a t =O.
E.S.P. This (2006 -2007)
111-1-6-2 : Modle du stockage froid
Les calculs sont analogues ceux appliqus au stockage chaud les hypothses restent
valables aussi au niveau du systme rservoir liquide froid.
La conservation de la masse donne :
d
dt MH =QHE +QHB - QH5
Soit:
MH =MHS +(QHES+QHBP - QHSS). t
QHB est le dbit massique soutir du dbit sortant de l'absorbeur.
Et la conservation de l'nergie s'exprime par :
(3.1.6.2.1)
(3.1.6.2.2)
d . .
dt (MH X cex TH) = QH X CP X THE +QHBX CP X THB - QHSX CP X THS - Wp
Soit :
TH = D +(THO - D) X e(CXt)
Avec:
D = (QHE X THE + QHB x TRB - WPJ!(QHE +QHBJ
. CP
c =_QHE+QHB
MH
TH =temperature du stockage "froid" t = 0
(3.1.6.2.2)
a: Pertes du stockage
Les pertes thermiques totales du stockage correspondent une diminution de l'nergie interne
du HITEC. Le bilan nergtique du rservoir au temps t s'crit:
ct
Wp = -.MH X CP x -TH (3.1.6.2.3)
dt"
Le flux de chaleur Wp entre le HITEC la temprature TH et l'air ambiant TAM est
conservatif et s'crit:
54
1
P.F.E. Prsentpar M. Mbaye SECK et M. Moussa DIOP
Wp =L(K X S) X (TH - TAM)
Avec:
----------
E.S.P. This (2006 -2007)
(3.1.6.2.4)
(3.1.6.2.6)
Icx XS) =Kt XSl +K2 XS2
O:
KI, K2 dsignent respectivement les coefficients d'change global du cylindre et des deux
disques.
SI, S2 d dsignent respectivement les surfaces d'changes du cylindre et des deux disques.
51 = 2-rrRL
~ Coefficient d'change global
Pour le cylindre, le coefficient d'change global rapport la surface interne SIs' crit:
2:... = ~ +.!!...ln RA +.!!...ln RC +.!!...ln RG + _R_
Kl KH A R C RA G RC RG.KO (3.1.6.2.5)
Pour les deux disques fermant le cylindre
2:... = ~ + EA + EC + EG +2:...
K2 KH A C G KO
Avec
Ki: le coefficient d'change global rapport la surface interne SI
w
KR et KO sont respectivement les coefficients de convection interne et externe en ~ 2
m .oc
K2 : le coefficient d'change global rapport la surface interneS2 en ~ o
m. C
R , RA, RG et RC rsistances thermiques respectivement global, de l'acier, laine
de roche et du calorifuge en m
2
. KJW
IlA ,IlC et IlG les coefficients de conductivit thermique de l'acier, du
calorifuge et du laine de roche en ~ .
m.OC
55
m
~ Coefficient d'change du HITEC
On se place dans le cas du liquide circulant faible vitesse dans un cylindre horizontal. La
relation de Colburn et Hougen donne alors:
Nu == 0.64x (GrGr.Pr )0.25
\D == 0.64 x (B.d
3
p Z:.CP.t:.t)0.25
Avec:
K : coefficient d'change exprim en
G: 9.81m/s
2
, intensit du champ terrestre
D : diamtre intrieur du tube
il. : conductivit du fluide
p : masse volumique du fluide
p: Coefficient de dilatation en volume
Il : viscosit dynamique du fluide
w/m OC
Kg/m
3
lIoC
Kg/m/s
tstt : Diffrence de temprature entre la paroi et le fluide OC
Nu : Nombre de Nusselt
Gr : Nombre de Grashof
Pr : Nombre de Prandtl
L'application numrique au HITEC donne :
KH == 4.4546 x (2082 - 0.7288 X TH) 0.5 X (TH - TA) 0.25 X Exp( -561.75) (3.1.6.2.6)
TH+273.15
~ Coefficient d'change externe
Le coefficient d'change externe comprend le coefficient de convection KF et le coefficient
de rayonnement KR [7].
Si l'on chauffe l'air la temprature ambiante l'aide du cylindre horizontal de diamtre D,
la relation de Colburn et Hougen donne [7]:
KF =6.16 X ( ~ ) 0.25 (3.1.6.2.7)
O D est exprim en mm.
Selon la loi de STEFAN-BOLTZMANN, le coefficient de transfert par rayonnement externe
entre le rservoir la temprature TG et l'air ambiant la temprature TAM peut s'crire
KR
=6.E(TG+273.15) 4-(TAM+273.15) 4 (3 1 6 )
.. .2.7
TG-TAM
56
Avec:
E.S.P. This (2006 -2007)
~ 6: constante de STEFAN-BOLZMANN =5.67.10 -8 w/m
2
K 4
~ E : Emissivit de la surface =0.26 pour l'acier galvanis.
Nous obtenons donc le coefficient d'change externe:
KO =KF+KR
En rsolvant l'quation obtenue en galisant (3.1.6.2.3) et (3.1.6.2.4)
-MH x CP x !!:..TH =L(K x S) x (TH - TAM) ,
dt
il vient:
TH =TAM + (THO- TAM.EXP(L(KXS) t)
MHxCPH
(3.1.6.2.8)
(3.1.6.2.9)
Nous obtenons ainsi la temprature du HITEC aprs refroidissement l'instant 1. Nous
pouvons alors calculer successivement les pertes thermiques Wp pour chaque pas de temps.
111-2) Modlisation du systme secondaire
Gnrateur de vapeur
vaporateur
Figure 12 : schma du systme secondaire
111-2-1 : Introduction
Ce systme se compose de deux circuits principaux: le circuit HITEC H et le circuit
EAU-VAPEUR. E-V . Le modle du premier fait intervenir la masse et la temprature
tandis que le modle du second fait, de plus intervenir la pression.
57
P.F.E. Prsent par M. Mbaye SECKet M. MoussaDIOP E.S.P. This(2006-2007)
Il s'agit d'un modle tabli en rgime permanent en crivant des quations algbriques qui
traduisent les bilans globaux aux bornes des lments.
Les bilans de masses s'expriment par des galits entre les dbits d'entre et de sortie. De la
mme faon, les bilans d'nergie sont effectus globalement entre les entres et les sorties des
lments.
La conservation de la quantit de mouvement conduit des quations algbriques traduisant
soit les effets de pertes de charge, soit les phnomnes de dtente dans les turbines.
Le modle dcrit dans ce chapitre permet de reprsenter, d'une manire simple, les diffrents
rgimes de fonctionnement du systme secondaire lorsque sa charge varie.
111-2-2 : Modle desEchangeurs
Un changeur se caractrise par les dbits et les tempratures des deux fluides l'entre et
la sortie. Nous n'tudions ici que l'aspect thermique des changes en nous limitant au rgime
permanant de fonctionnement.
Ill-2-2-1 : Description du modle
Ce modle permet de dterminer les nouvelles conditions de fonctionnement d'un changeur
en connaissant un tat de rfrence (rgime nominal de fonctionnement) qui est dfini par les
dbits des deux fluides et les quatre tempratures. Pour calculer le nouvel tat, il faut
connaitre les deux dbits et deux des quatre tempratures.
Normalement, le coefficient d'change global K est un paramtre local, qui peut varier le
long de l'changeur. Pour simplifier le modle, nous admettons que ce coefficient reste
constant tout le long de l'changeur. D'ailleurs, nous posons que les capacits calorifiques
massiques des fluides ne varient pas d'un pas d'un tat de fonctionnement l'autre. La
nomenclature suivante est relative aux quations qui traduisent le modle dcrit.
CPE
CPH
de
di
NOMENCLATURE
capacit calorifique du fluide froid (eau ou vapeur)
capacit calorifique du fluide chaud (HITEC)
diamtre extrieur du tube
.diamtre intrieur du tube
58
Jlkg.OC
Jlkg.oC
mm
mm
dh
Di
diamtre hydraulique
diamtre intrieur de la calandre
mm
mm
Fe, fi facteur d'encrassement extrieur et intrieur
K, KG coefficient d'change global
KE, KH coefficient .des films externe et interne
KEI coefficient du film externe longitudinal
espacement des chicanes
pas pas des tubes
PM primtre mouill
QE, QH dbit massique du fluide froid et du fluide chaud
re, Tl rayons extrieur et intrieur du tube
Pr nombre de Prandtl
Nu Nombre de Nusselt
Re Nombre de Reynolds
Sp .section de.passage
S surface totale d'change
TEE, TES temprature du fluide froid l'entre et la sortie
THE, THS temprature du fluide l'entre et la sortie
TM cart logarithmique des tempratures
V vitesse du fluide
W puissance change
p masse volumique du fluide
Conductivit du fluide
A conductivit du matriau du tube
!l viscosit dynamique du fluide la temprature du fluide
. .
ut viscosit dynamique du fluide la temprature de la paroi.
w
w
mm
mm
mm
kg/s
mm
mis
watt
kg/rn"
W/m.oC
W/m.oC
kg/mis
kg/mis
111-2-2-2: Equations exprimant les changes
Nous allons utiliser un changeur contre courant a cause de son rendement plus lev que
celui d'un changeur courants parallles
La chaleur change We travers la surface S peut s'exprimer de trois manires [7] :
59
P.F.E. Prsent.par M. MbayeSECKet M. MoussaDIOP E.S.P. This(2006-2007)
We = K X S X TM
We = QH X CPH X (THE - THS) Chaleur cde par1e HITEC
We = QE X CPE X (TES - TEE) Chaleur reue par la vapeur
Avec:
TM - (THE-TES)-(THS-TEE) _ (THE-THS)-(TES-TEE)
- THE-TES - THE-TES
Log(THS-TEE) Log(THS-TEE)
Ou encore:
(3.2.2.2.1)
(3.2.2.2.2)
(3.2.2.2.3)
(3.2.2.2.4)
. 1 1
TM
- we(QHxCPH QEXCPE)
- THE-TES
Log(THS-TEE)
En introduisant la valeur de TM, nous obtenons donc :
L (THE-TES) KXS QExCPE 1
Dg THS-TEE = QEXCPE X (QHXCPH - )
Posons:
R =QEXCPE
qHXCPH
KXS
U = Exp( X (R - 1)
QFxCPE
Il en rsulte:
TH - TEE = (THE - TES)/U
En galisant (3.2.2.2.2) et (3.2.2.2.3), il vient la relation:
THE - THS =R X (TES - TEE)
(3.2.2.2.5)
(3.2.2.2.6)
(3.2.2.2.7)
Connaissant les deux tempratures l'entre THE et TEE , nous pouvons dterminer les deux
tempratures la sortie suivantes.
D'aprs (3.2.2.2.6) et (3.2.2.2.7), nous pouvons crire:
TES = TEE + (THE - TEE)
En introduisant l'quation prcdente dans (3.2.2.2.6), il vient:
THS =TEE + ~ . x (THE - TEE)
UxR-l
60
P.F.E. Prsentpar M. Mbaye SECKet M. MoussaDIOP E.S.P. This (2006-2007)
a: Coefficient d'change global
Le coefficient d'change global reprsente la conductance thermique interpose entre les
deux fluides. Il comprend principalement trois termes [7]:
- coefficient de convection fluide chaud-paroi, KR,
- coefficient de conduction travers la paroi, en fonction de la conductivit du
matriau et de l'paisseur de la paroi,
- coefficient de convection paroi-fluide froid, KF.
Pour le cylindre, le coefficient d'change rapport la surface externe S s'crit:
1 re 1 re . re (re) 1
KG =ri X KH + ri X tl + a X Log ri +te + KE (3.2.2.2.8)
KG est le coefficient d'change global avec encrassement. Lorsque l'changeur est neuf, les
rsistances d'encrassement intrieure fi et extrieure fe sont nulles et on dfmit un coefficient
d'change global p r o p ~ e K en ngligeant le terme de conduction:
2. = re X ~ +.2...
K ri KH KE
b: Coefficient d'change KH travers le film intrieur
La formule gnrale de Colbum s'crit [7]:
Nu = 0.023 x Re
O
.
8
x Pr 0.4
Avec:
Re =p.V.di = 4.QH
Il n.di1l
P
IlCP
r=T
Nous avons donc :
KH = ~ Nu = 0 023 x ~ x Re 0.8 x Pr 0.4
di' di
Ou encore:
KH = 0.023 x ~ x (_4._) 0.8 X Pr 0.4 x QH 0.8
dl n.dl.1l
En rgime' nominal :
KHN =0.023 x ~ x (_4._) 0.8 X Pr 0.4 x QH 0.8
dl n.dl.1l
Et pour tout tat en faisant le rapport des deux quations prcdentes on a :
61
(3.2.2.2.9)
(3.2.2.2.10)
KH =KHN (QH) 0.8
QHN
E.S.P. This(2006-2007)
(3.2.2.2.11)
..
c: Coefficient d'change KE travers le film externe
~ Ecoulement sans chicane
Considrons que la formule gnrale de Colburn reste valable dans le cas d'un coulement
longitudinal condition d'utiliser le diamtre hydraulique dh, au lieu de di. Il vient donc [7]:
"KEI =0.023 x ~ x Re 0.8 x Pr 1/3
dh
O:
dh =4xsection de passage =4.sp
perimetre mouill PM
Re = pxvxdh = 4XQE = QExdh
Il Il
XPM
Il
XSp
Pr =Il
XCP
~ Ecoulement en prsence de chicanes

7 7 7 7
Di!
Figure 13 : Ecoulement en prsence de chicane
(3.2.2.2.12)
(3.2.2.2.13)
(3.2.2.2.14)
pour: 2000 < Re < 40000
Nous utilisons les formules donnes par Grgorig concernant les courants croiss; celle
relative au nombre de Nusselt est:
Nu =K1X,K2 X Re m x Pr 0.31
0.908 '
Avec:
Re =p.V.de =QExde
Il IlXSp
Pr = Il.
cP
62
KE=-xNu
de
s - lx(pas-de)xDi
P - Pas
O:
E.S.P. This (2006-2007)
(3.2.2.2.15)
1: espacement des chicanes
Pas : pas des tubes
De : diamtre intrieur de la calandre
Di : Diamtre intrieur de la calandre
KI : facteur d l'inclinaison du courant moyen par rapport l'horizontale.
K2 et m sont des constantes.
111-2-3 : Modle du gnrateur de vapeur
NOMENCLATURE relative ce modele
NOMENCLATURE
EV dsigne vaporateur, GV dsigne gnrateur de vapeur, EQ conomiseur et ZU
surchauffeur.
QHGV est le dbit massique du fluide HITEC dans le gnrateur de vapeur
QEEGV est le dbit massique d'entre de l'eau dans le gnrateur de vapeur
QEEEQ est le dbit massique d'entre de l'eau dans l'conomiseur.
QEEEV est le dbit massique d'entre de l'eau dans l'vaporateur.
QVEZU est le dbit massique d'entre de la vapeur dans le surchauffeur.
QVSGV est le dbit massique de sortie de la vapeur dans le gnrateur de vapeur.
THE et THS dsignent les tempratures d'entre et de sortie du HITEC de mme TEE
et TES dsignent les tempratures d'entre et de sortie de l'eau. TV dsigne
temprature de vapeur.
Ces notations de tempratures peuvent tre suivies d'une notation d'appareil
comme EQ qui donne par exemple THEEQ qui dsigne la temprature d'entre du
fluide HITEC dans l'conomiseur. On peut ajouter la lettre N pour signifier
nominal donc:
TEEEV est la temprature d'eau entrant dans l'vaporateur.
63
P.F.E.Prsentpar M. MbayeSECKet M. MoussaDIOP E.S.P. This (2006-2007)
TEEEQ est la temprature d'eau entrant dans l'conomiseur.
TVSEV est la temprature de vapeur sortante de l'vaporateur.
TSV est la temprature de vapeur sortante de du gnrateur de vapeur.
THEGV est la temprature du HITEC entrant dans le gnrateur de vapeur.
THSEV est la temprature du HITEC sortant de l'vaporateur
THSGV est la temprature du HITEC sortant du gnrateur de vapeur
THSEQ est la temprature du HITEC sortant de l'
TESEQ est la temprature d'eau sortante de l'conomiseur
THEEQ est la temprature du HITEC entrant dans l'conomiseur
TEMEV est la temprature moyenne de HAUSBRAND de l'eau dans l'vaporateur.
TMEQ est la temprature moyenne de HAUSBRAND dans l'conomiseur.
TMZU est la temprature moyenne de HAUSBRAND dans le surchauffeur.
THEEV est la temprature du HITEC entrant dans l'vaporateur.
KEQ et SEQ sont le coefficient global d'change et la surface d'change dans
l'conomiseur
KEV et SEV sont le coefficient global d'change et le surface d'change dans
l'vaporateur
KZU le coefficient global d'change dans le surchauffeur
WZU est le travail chang travers la surface d'change SZU
WEQ est le travail chang travers la surface d'change SEQ
WEV est le travail chang travers la surface d'change SEV
CPVZU est la capacit calorifique de la vapeur dans le surchauffeur.
~ Hypothses
Le gnrateur de vapeur comprend l'conomiseur, l'vaporateur, et le surchauffeur
l'intrieur d'une seule calandre.
Le HITEe liquide circule dans les tubes et cde ses calories l'eau sous-sature pour la
transformer en vapeur surchauffe.
Pour simplifier le modle, nous avons pos les hypothses suivantes:
64
P.F.E.Prsentpar M. MbayeSECKet M. MoussaDIOP E.S.P. This (2006-2007)
Les rgimes hydrauliques des deux fluides sont considrs comme quasiment
permanents. A chaque instant, le gnrateur de vapeur est travers par deux dbits
QHGV et QEEGV uniformes tout au long de l'changeur. Nous avons donc:
QEEGV=QEEEQ=QEEEV=QVEZU=QVSGV
La temprature du HITEC THEGV l'entre du gnrateur de vapeur est constante et
gale 400C.
L'vaporateur EV ne contient que de l'eau en bullition pression constante d'o:
TEEEV=TVSEV=TSV
Le coefficient d'change global K varie comme la puissance 0.8 du dbit de vapeur.
D'aprs (3.2.2.2.11):
K =KN x ( QV ) 0.8
QVN
La chaleur massique de la vapeur CPV est constante et gale au rapport de l'enthalpie
la temprature en rgime nominal:
CPV = HVSN-HVEN
TVSN-TVEN
~ Description du modle
Le modle des changeurs tabli au paragraphe prcdent permet de dterminer les diffrents
, ,
rgimes de fonctionnement du gnrateur de vapeur.
Nous avons 4 variables d'entre: les deux tempratures d'entre THEGV, TEEGV, la
pression PGV et le dbit d'eau QEEGV qui correspond au besoinde la turbine.
Ce modle permet de calculer les tempratures des deux fluides l'entre et la sortie de
l'conomiseur, de l'vaporateur et du surchauffeur [7].
111-2-3-1: L'conomiseur EQ
D'aprs les quations d'change, nous pouvons crire en adaptant les notations:
WEQ =KEQ x SEQ x TMEQ =KEQN x SEQ X TMEQ X (QEEGV) 0.8
QEEGVN
=QBEEQ X CPEEQ X (TESEQ - TEEEQ)
=QHGV X CPH X (THEEQ - THSEQ)
Ou encore:
WEQ = QEEEQ X (HSE - HEEEQ) = QHGV X CPH X (THSEV - THSGV)
L'cart de temprature ~ T H E Q au niveau du HITEC s'crit:
65
(3.2.3.1)
(3.2.3.2)
(3.2.3.3)
(3.2.3.4)
P.F.E. Prsent par M. Mbaye SECK et M. MoussaDIOP
~ T H E Q =THEEQ - THSGV =QEEGVX(HSE-HEEEQ)
. .' QHGVxCPH
E.S.P. This(2006-2007)
(3.2.3.4)
(3.2.3.5)
m-z-s-z . L'vaporateur EV
La temprature de l'mulsion eau-vapeur reste constante et gale la temprature saturante
TSV, nous avons donc:
TVSEV=1EEEV=TSV
En adaptant les notations, les quations peuvent s'crire:
WEV = KEV x SEV x TEMEV = QEEEV x LEV = QHGV x CPH x (THEEV - THSEV) et
TMEV - THEEV-TH5EV
- THEEV-TSV
Log(THSEV-TSV)
En raisonnant comme pour la modlisation des changeurs. Nous pouvons crire:
L
THEEV-TSV KEVxSEV
Dg -
THSEV-TSV QHGVxCPH
Posons:
Ul = EXpC-KEVXSEV)
QHGVxCPH
!1THEV=THEEV-THSEV
Il s'ensuit :
THSEV = TSV +!J.THEV x ...E2....
l-Ul
TEEV =THSEV +DTHEV
Calculons le coefficient d'change global de l'vaporateur en rgime nominal:
KEVN =QEEEVNxLEVN
SEVxTMEVN
D'aprs (3.2.2.2.11), nous avons:
KEV = KEVN X CQEEGV ) 0.8
QEEGVN
(3.2.3.6)
(3.2.3.7)
(3.2.3.8)
(3.2.3.9)
(3.2.3.10)
Ill-Z-3-3 : te surchauffeur ZU
En effectuant les calculs utiliss pour modliser les changeurs, nous obtenons:
WZU =KZU x SZU x TMZU
= QEEGV x CPVZU x (TVSGV - TSV) =QEEGV x (HVSGV- HSV)
= QHGVx CPH X (THEGV - THEEV) (3.2.3.11)
66
TMZU =(THEGV-THEEV)-(TVSGV-TSV)
L (THEGV-TVSGV)
og THEEV-TSV
Posons:
U2 =EXP( KZUxSZU X (QEEGVXCPVZU - 1))
. QEEGVXCPVZU QHGVxCPH
Il en rsulte:
TVSGV = THEGV - (THEEV - TSV) x U2
Nous avons:
KZUN =QEEGVNX(HVSGVN-HSVN)
SZUxTMZUN
Et d'aprs (3.2.2.2.12) :
KZU =KZUNx ( QEEGV ) 0.8
QEEGVN
Ou encore:
KZU X SZU = KZUN X SZUN X (QEEGV) 0.8
QEEGVN
. .
KSZU =KSZUN X ( QEEGV ) 0.8
QEEGVN
E.S.P. This (2006-2007)
(3.2.3.12)
(3.2.3.13)
(3.2.3.14)
(3.2.3.15)
Et
CPVZU = HVSZUN-HSVN
TVSZUN-TSVN
111-2-4: Modle de la turbine
(3.2.3.16)
111-2-4-1 : Introduction
C'est une turbine action de type axial 1 tage tournant 3 000 tr/mn.
~ Hypothses
Nous envisageons un modle simplifi ne faisant intervenir que des valeurs globales des
diffrents paramtres' de l'coulement: le dbit de vapeur QV, la pression PV et la
temprature TV.
Nous considrons que le rendement isentropique et le rendement interne de la turbine sont
constants et gaux leur valeur en rgime nominal.
111-2-4-2 : Equation de la chute isentropique
Une dtente isentropique correspond un coulement thorique rversible ou l'change
calorifique avec l'extrieur est nul (Q=O).
67
En transformant l'expression du premier principe de la thermodynamique afin de n'y faire
apparaitre comme travail que le travail externe sur l'arbre WI, nous trouvons dans le cas d'un
coulement adiabatique permanent dans un canalmobile :
WI + Q= f dh + II v
2
. 2
Ou encore:
WI =f dh + II v
2
2
(2.4.2.2.1)
(2.4.2.2.2)
En pratique, l'nergie cintique du fluide aux deux extrmits de la machine est petite et plus
. v2
souvent encore, de mme importance. De la sorte ll- O.
2
Le terme WI s'identifie alors la variation d'enthalpie statique f dh. il vient donc:
WI =f dh = HE - HSS (2.4.2.2.3)
Dans une transformation adiabatique rversible ne possdant pas ncessairement le caractre
de performance, il vient :
dh =dU +d(pV) =dU + pdV + VdP
dh = dW + dQ ,:. pdV + Vdp = VdP
Nous avons donc:
(2.4.2.2.4)
(2.4.2.2.5)
WI =f VdP (2.4.2.2.6)
Avec PV
N
=PE. =PS. VS
N
=constante, nous trouvons la chute isentropique WI :
WI =f:
S
E
VE. dP
n-
N
WI = HE - HSS = - x PE x VE x (1 - - )
N-l PE
(2.4.2.2.1)
111-2-5: Le stockage de l'hydrogne
De faon gnrale, le stockage de l'hydrogne est un vritable problme car pour la plupart
des systmes consommant de l'hydrogne, la production de ce gaz ne se fait pas sur place.
Des mthodes scuritaires de stockage d'hydrogne en grande quantit s'imposent. Dans le
cas de notre installation, le stockage en gaz comprim suffit, car la quantit stocke n'est pas
trs importante et le cot est moindre compar aux autres types de stockage. La compression
de l'hydrogne gazeux est donc un moyen d'augmenter la quantit stocke dans un volume
fixe de rservoir. Pour cela, on utilise de rservoir de forme cylindrique en matriau
composite renforcs avec de l'aluminium.
Nous remarquons que la pression de stockage ncessaire pour l'alimentation correspond la
pression de remplissage du rservoir.
68
Donnes du rservoir :
Pression maximale supportable par le rservoir Pmax = 250 bars
Volume du rservoir V = 50 m
3
Dure de vie du rservoir : 15 ans
Dure de fonctionnement du rservoir = 24h/24h
Prix du rservoir en moyenne = 1 189 000 Fcfa
E.S.P. This (2006-2007)
Donnes sur les conduites :
Nous aurons une conduite liant l'lectrolyseur et le rservoir de stockage. La conduite peut
tre approxime 10 m de longueur et 2 cm de diamtre.
Prix par mtre de ce type de conduite avec protection extrieure l'aluminium gale Il 429
Fcfa/m.
Dure de fonctionnement de la conduite: 24h / 24h
Dure de vie de la conduite : 15 ans.
111-3 : Dimensionnement
La mthodologie consiste estimer le besoin en hydrogne, puis calculer la puissance
nominale de l'lectrolyseur, puis la charge globale en incluant la puissance absorbe par les
moteurs et en fin de proche en proche, on dterminera les puissances des appareils de
production d'nergie.
111-3-1 : Dtermination de la charge
Au Sngal 79.2m
3
d'hydrogne sont vendus en moyenne par AIR LIQUIDE et cette
quantit quivaut celle qui est vendue dans les pays de la sous rgion comme le MALI, la
MAURITANIE, le BURKINA etc.
Notre installation devra donc prendre en compte ce march potentiel plus celui des pays
limitrophes qui n'ont pas t cits encore non cits.
Nous estimons la charge globale 2000m
3
111-3-1-1: Puissance de l'lectrolyseur

Les lectrolyseurs ont en gnral un rendement de 0.8. Des donnes concernant des
lectrolyseurs industriels donnent les indications suivantes: on a besoin environ de 5 kWh
pour produire 1 m
3
d'hydrogne.
Donc pour ces 2000m
3
on aura besoin de: Ee = 10.000 kWh par an
69
Pe =5 x 2000
52x5 x 8
Pe=4.81 KW
puisque l'unit ne fonctionne que 8b/jour. Par suite,
La puissance absorbe par l' lectrolyseur est de :
Pe
P
ab
=-
1)
P
ab
=4.81/0.8 =6 KW
111-3-1-2: Quantit d'eau ncessaire
pour produire un m
3
d'hydrogne, il nous faut 0.81 d'eau distille.
On aura besoin donc de ~ a u = 0.8 x 2000
Veau = 1600 l
111-3-1-3: La puissance des moteurs pour la compression
L'nergie requise pour la compression de l'hydrogne de la pression ambiante une pression
de 20 MPa est de l'ordre de 7 % de l'nergie contenue dans l'hydrogne qui est de 3.55 kWh /
m
3
La puissance des moteurs pour la compression est donc :

3.55
Pcom = 2000 x 52 x 5 x 8 x 0.07
P
com
=0.24 Kw
111-3-2 : La charge totale
La charge totale est constitue de l'lectrolyseur, de l'clairage et des moteurs asynchrones
(rendement. = 0,76) qui alimentent les compresseurs (2) et les pompes (4). Nous en
dnombrons une entre le stockage froid est le miroir, une entre la source d'eau et le
compresseur, une entre le compresseur et le gnrateur de vapeur. La dernire pompe
alimente l'lectrolyseur.
Puissance lectrolyseur Pab = 6 KW
Puissance des moteurs P. =6 X 0.24 =232 KW
asy 0,62 '
Divers = 15% (6+2,32) = 1,25 KW
Puissance correspondanr la charge totale est :
PT =Pab + Pasy + Pd= 9, 57KW
70
P.F.E. Prsentpar M. Mbaye SECKet M. MoussaDIOP
C'est la puissance requise la sortie du transformateur.
E.S.P. This (2006-2007)
..
111-3-3 : Lapuissance des autres appareils
Pour la suite, du processus de dimensionnement des autres lments de notre unit, nous
avons recens les rendements moyens thoriques des diffrentes parties de l'unit afin
d'estimer la puissance de chaque appareil. Ces informations sont regroupes dans le tableau
suivant:
libell Rendement Puissance (KW)
Charge totale 9.57
transformateur 0.95 10.07
alternateur 0.95 10.6
Cycle 0.4 26.51
thermodynamique
absorbeur 0.4 66.3
Tableau 4 : Puissance des appareils
En dfinitif, nous aurons le rendement global thorique de toute l'installation:
_ Puissanceentre _ 9.57
1'Jglobal
Puissancesortie 66.3
1'Jglobal =15 %
Puisque, nous connaissons la puissance par unit d'aire de la parabole, nous allons dterminer
maintenant la longueur minimale (Hmin) qui nous permettra de fournir notre turbine la
puissance optimale pour vaincre la charge.
Pour une unit d'aire A = 1 m
2
, P
-_.. 2 xsin9 xl
on a une puissance
E
Donc, pour la
puissance souhaite P = 64.3 KW, nous aurons une surface A = 1t x A'B' xH.
71
P.F.E. Prsent par M. Mbaye SECKet M. Moussa DIOP
Donnes:
A'B' : diamtre du tube = 2 cm
E : ouverture =0.0093
1* : la luminance = 1000
e :angle =90
Ceci nous permettra de calculer (Hmin) :
H . 66.3xE
mm 2xAf8fxsinOx7fx[*
. H - 66.3Xl0
3XO.0093
A.N. min - ----.----
2xO.02xsmOx7fxl000
H
mi n
= 4.906 m ~ 5 m
III-4 : Conclusion partielle
E.S.P. This(2006-2007)
Cette modlisation n'est pas exhaustive, elle n'a concern que certaines parties du
systme a cause du manque de temps. Une simulation aussi n'a pas t faite pour les mmes
raisons.
On note la prsence de beaucoup d'inconnus cause d'un manque de donnes trs
significatives sur plusieurs appareils du systme.
72
P.F.E.Prsent par M. Mbaye SECKet M. MoussaDIOP E.S.P. This (2006-2007)
CHAPITRE IV : ETUDE FINANCIERE
IV1 : Etude de cot
Avant de faire une analyse de rentabilit financire, il faut d'abord avoir une ide sur les cots
et les recettes. Les cots des divers lments ont t obtenus auprs des fournisseurs aprs des
recherches li travers l'internet.
=
=
230*10*660 = 1 518000 Fcfa
= 26400 Fcfa
30*660 = 19800 Fcfa
33000 Fcfa
= 164825740 Fcfa
= 297000 Fcfa
400000 Fcfa
= 15 552 000 Fcfa
10070 W * 3600 = 32860800 Fcfa
= 1 189000 Fcfa
171 435 Fcfa
150480 Fcfa
230340 Fcfa
5600*660 = 3 696 000 Fcfa
= 8000 Fcfa
=220 977 995 Fcfa
= 22 097 796 Fcfa
=243 075 791 Fcfa
15m * 11429 =
228*660 =
10 070 W * 40*0.62*660
3*150*660
IV-1-1: Cot des quipements de production
Miroir cylindro-parabolique :
Tube sous vide: 40*660
Alimentation stabilise: EURO-CB EPS-1012
Alternateur: A13N137 50*660
Turbine:
Pompe (3):
Stockage (2): 2* 200 000
Gnrateur de vapeur :=
Electrolyseur :
Rservoir de stockage pour l'hydrogne:
Conduites liants l'lectrolyseur au rservoir:
Distillateur d'eau: AQUATRON rf: 75338
Suiveur solaire avec rgulateur 349*660
Transformateur 24KV/358V
Matriels lectriques
Total
Cot d'installation (10%)
Cout total
IV-1-2 : Cot du matriel et mobilier de bureau
Voiture de service
2 ordinateurs 2*400000
Une photocopieuse OKI FAX 5680
6 tables de bureau 6*30 000
= 7000000 Fcfa
8000000 Fcfa
= 500000 Fcfa
= 180000 Fcfa
73
15 chaises 15*5000
2 armoires de rangement 2*50000
=
=
75000 Fcfa
100000 Fcfa
Total = 15 855 000 Fcfa
lV-1-3: Cot des matires premires
Les donnes suivantes sont les cots unitaires de matires premires. Elles permettront de
dterminer les cots de matires premires pour chaque anne. Il s'agit des prix.
Electrodes en nickel lourd (2/2 jours) : 52*5*2000Fcfa
Carburant (vhicule) :(30*0,1*5*52) litres*550Fcfa
Hydroxyde de potassium (0,8*2000 *20%) litres*2 OOOFcfa =
Divers
Total =
520000 Fcfa
= 429 000 Fcfa
640000Fcfa
50000 Fcfa
1639000 Fcfa
IV-1-4: Cot du personnel
Nous estimons les salaires mensuels bruts du personnel comme suit:
Directeur (ingnieur) 500000 Fcfa
Deux ouvriers 2*80000 160000 Fcfa
Secrtaire 125 000 Fcfa
Technicien de surface 70000 Fcfa
Gardien 75000 Fcfa
Total 930000 Fcfa/mois
Charge sociale (10%) 93000 Fcfa
Total 1 023000 Fcfa
Donc le montant total annuel des salaires bruts sera: 12276000 Fcfa/an
IV-1-5: Cot des quipements divers et de maintenance
Nous estimons ces cots 2 000000 de Fcfa par an. Ces quipements concernent les outils
(cl, marteau, pince... ), les bouteilles d'huile ...
74
P.F.E. Prsent parM. MbayeSECKet M. MoussaDIOP E.S.P. This (2006-2007)
IV-1-6: Cot de construction Gnie civil
A partir des donnes fournies par l'Institut Sngalais de Normalisation values en mtre
carr bti, nous allons estimer le cot de la construction 120000 Fcfa pour un btiment
industriel et 80 000 Fcfa pour un btiment usage domestiques et 17 500 Fcfa pour la voirie
(bitume).
IV-2 : Etude financire
L'tude financire a pour but de veiller la rentabilit des capitaux investis et l'quilibre
financier des services chargs de grer le projet.
Sur la base des informations disponibles, l'analyse financire cherchera savoir si l'entreprise
peut recouvrer les fonds investis partir des recettes gnres par le projet.
20000000 Fcfa
243 075 791 Fcfa
15 855 000 Fcfa
2 000 000 Fcfa
280930791 Fcfa
=
=
=
=
Total
IV-2-1 : Les investissements
Un investissement est un engagement des ressources d'une organisation ralis dans le
prsent dans l'espoir de raliser des bnfices futurs.
Avant tout investissement dans un projet, il faut s'assurer de sa rentabilit et de sa faisabilit.
Pour notre unit, nous allons considrer la rentabilit conomique.
L'investissement de notre projet se prsente comme suit:
Construction Gnie civil
Equipement de production
Matriel et mobilier de bureau
Matriel divers
IV-2-2 : Les projections financires
Le taux d'imposition actuel appliqu aux entreprises sngalaises est de 25 % et nous allons
considrer dans la suite de nos calculs un taux d'actualisation de 12%. En effet, le taux
d'actualisation est le taux de rendement employ pour convertir une somme, payable ou
recevable ultrieurement, sa valeur actuelle.
7S
P.F.E. Prsent par M. Mbaye SECKet M. MoussaDIOP E.S.P. This (2006 -2007)
IV-2-3 : L'amortissement (annexe 2)
L'amortissement est une constatation comptable de la perte subie sur la valeur d'actif des
immobilisations qui se dprcient par l'effet du temps, de l'usage et du changement
technologique.
C'est une dpense qui n'entraine aucune sortie d'argent; elle est donc dductible d'impts;
elle influence les flux montaires en diminuant l'impt pay.
L'amortissement dgressif peut tre utilis pour les biens d'quipements neufs et autres sauf
les immeubles d'habitation, les chantiers et les locaux servant l'exploitation de la profession
acquis ou fabriqus depuis le 1erjanvier 1987 par les entreprises industrielles.
Nous allons l'utiliser dans le cadre de ce projet, car il varie avec le temps permettant de payer
moins d'impts lors des premires annes o les recettes sont faibles et plus d'impts pendant
les dernires annes o les recettes sont plus substantielles.
Dans l'amortissement dgressif le taux d'amortissement est obtenu partir du taux linaire
affect d'un coefficient gal :
2 si la dure normale d'utilisation est infrieure ou gale 5 ans;
2.5 Si la dure normale d'utilisation est suprieure ou gale 5 ans.
Lorsque l'amortissement calcul partir du taux devient inferieur au quotient de la valeur
rsiduelle sur le nombre d'annes d'utilisation normale restant courir, il convient de
terminer le tableau d'amortissement en prenant constantes gales au quotient de la valeur
rsiduelle sur le nombre d'annes d'utilisation restantes.
IV-2-4: Les recettes d'exploitation (annexes 3 et 5)
Elles reprsentent le montant des ventes qui s'est concrtement traduit en cash. Pour pouvoir
les dterminer, il faut connaitre le prix de vente unitaire des produits.
C'est ainsi que nous allons dterminer le prix de vente unitaire en nous basant sur la notion du
seuil de rentabilit. En effet, le seuil de rentabilit nous permettra d'avoir une approximation
du prix de revient afin d'en dduire un prix de vente. Il faut noter que cette mthode n'est pas
exacte car elle ne tient pas compte de la dprciation de la valeur dans le temps.
IV-2-5 : Dtermination du prix de vente du produit
Nous supposons que la vente annuelle est constante et est gale 2 000 m
3
par annes.
76
P.F.E. Prsent-par M. Mbaye SECK et M. Moussa DIOP E.S.P. This (2006 -2007)
D'aprs les informations obtenues travers AIR LIQUIDE, le prix de vente du mtre cube
(m') d'hydrogne (hors taxe) sur le march sngalais est de 29 500 Fcfa. Soit:
CV : cot variable;
CF : cot fixe ;
Q : quantit de produit.
IV-2-S-1 : Dtermination des cots variables
Ce sont les cots qui sont lis directement au volume de production. Ils se composent
essentiellement des cots de l'nergie, de maintenance et de la matire premire. Nous aurons
donc:
Equipement divers et maintenance :
Cot annuel matires premires/Q :
Main d'uvre 355 500/2000
Total = CV
2000000/5/2000 = 200 Fcfa
1 639000/2 000 = 820 Fcfa
= 167,75 Fcfa
= 1187,75 Fcfa
IV-2-S-2 : Dtermination des cots fixes
Ce sont tous les cots qui sont indpendants de la quantit produite, ils sont constitus de la
dprciation dans le temps et des salaires du personnel.
Amortissement
Gnie civil 20 000 000 / 20
quipement de production 243 075 791 / 10
Cot de fonct +mat divers 17 855 000 / 5
Total
Personnel de la direction :
Total = CF
1 000 000 Fcfa
24 307 591 Fcfa
3 571 000 Fcfa
28 878 591 Fcfa
625000 Fcfa
29 503 591 Fcfa
1
Nous fixons le prix de vente unitaire de notre produit 29500 Fcfa, car par rapport la
concurrence, le prix du march est de 29 500 Fcfa. On calculera donc toutes recettes partir
de ce prix..
77
IV-2-6: La rentabilit du projet
Nous utilisons deux critres de rentabilit des projets d'investissements industriels le DRC et
la VAN.
IV-2-6-1 : Le Dlai de Rcupration du Capital Investi (DRe) (annexe 1)
C'est le dlai ncessaire pour que les recettes du projet quilibrent le montant des dpenses
d'investissement.
debourss requis
DRC =-------------
total des encaissements annuels nets
Pour ce projet, nous avons suppos que les entres de fond sont constantes d'une anne une
autre. Ce critre teint compte de l'impact que peut avoir un projet sur la liquidit de
l'entreprise. C'est aussi une faon de tenir compte du risque, si on accepte l'hypothse que le
risque augmente avec le temps.
Pour notre tude, le capital investi est de 280 930 791 Fcfa,
A la fin de la dix-huitime anne, le capital rcupr est de 279902029 Fcfa. Il reste donc
rcuprer 1 028162 Fcfa.
Le temps requis pour rcuprer cette somme est de : 1 028762/716277.6
Avec 6 629 835,39 le flux net actualis de la quatorzime anne.
Le dlai de rcupration du capital est donc de 18, 14 ans.
IV-2-6-2 : La Valeur Actuelle Nette (VAN) (annexes 3)
Elle s'value par la valeur actualise des rentres de fonds moins la valeur actualise des
sorties de fonds pour la dure du projet. Elle traduit aussi la richesse supplmentaire apporte
par le projet.
VAN =-investissement
+Lrecettes nettes x (; , i%,n) - Lcots nets x (; , i%,n)
n n
1% =12% = taux d'actualisation. C'est le taux d'intrt dont on se sert pour dterminer la
VAN d'un montant futur en l'actualisant.
n = dure du projet;
P = valeur actuelle;
Fn = valeur.future en l'anne n ;
78
(p/Fn, i%, n) = valeur actuelle P du montant futur Fn l'anne n au taux de i%.
(p/Fn, i%, n) = li (l +it
Si VAN > 0 le projet est rentable
Si VAN < 0 le projet n'est pas rentable
Si les projets sont exclusifs, choisir celui dont la VAN est la plus leve mme si les valeurs
sont ngatives.
Il existe deux mthodes pour calculer la VAN :
La mthode fiscale ;
La mthode des cash-flows que nous avons utilise, car plus aise pour notre cas.
VAN =Valeur actuelle des cash-flows - investissement.
En appliquant cette relation, nous avons trouv VAN = - 76 197 373, 5 Fcfa.
IV-2-' : Etude financire en tenant compte de l'oxygne produite par
l'lectrolyse (annexes 4 et 5)
Dans cette partie nous allons considrer les 1 000 m
3
d'oxygne que nous avons produit la
suite de l'lectrolyse. Le prix concurrentiel que nous utiliserons dans nos calcul est de
1.2*660 = 792 Fcfa hors taxe.
Donc d'aprs l'annexe x nous avons:
DRC = 17,95
VAN = - 75 483 039, 5 Fcfa
IV-3 : Conclusion partielle
Au vu des rsultats de cette analyse de rentabilit financire, il en rsulte que le projet n'est
pas du tout rentable car tous les indicateurs utiliss le confirment.
79
P.f.E. Prsent par M. Mbaye SECK et M. Moussa DIOP
CONCLUSION GENERALE
E.S.P. This (2006 -2007)
Le niveau atteint par la communaut scientifique dans .le domaine de la chimie,
l'lectrochimie et l'lectrolyse en particulier reste encore un niveau bas surtout par rapport
au besoin pressant de remplacement des nergies fossiles par les nergies propres.
Cette tude nous a permis de comprendre combien il est encore difficile de produire de
l'hydrogne en grande quantit. La production par lectrolyse demande des quipements trs
complexes, et une connaissance scientifique pluridisciplinaire. C'est ainsi qu'on a eu
l'occasion de mettre en application beaucoup de nos connaissances acquises lors de notre
formation. Elle nous a.rvl toutes les difficults qu'il y a mener une tude de bout en bout,
surtout quant il s'agit de faire une modlisation sans vraiment disposer toutes les donnes. Ce
manque de donnes nous a empch d'aller jusqu'au bout de l'tude de certains lments de
la chane comme la turbine, l'alternateur, le gnrateur de vapeur et l'lectrolyseur.les
caractristiques de ces lments sont obtenus en gnral par commande accompagne des
donnes de fonctionnement.
Nous dplorons vraiment l'absence dans ce rapport de simulation informatique qui
aurait donn une vision plus claire du comportement et de la performance du systme ou
certaines grandeurs au fonctionnement. Le manque de temps en est la principale cause, car il
.
faut noter que les cours de la cinquime anne sont trs importants, nous les de prs et
faisons beaucoup de dplacement pour chercher des donnes obtenir.

Il faut noter que c'est trs coteux de produire de l'hydrogne par lectrolyse,
l'investissement est trs lourd mme avec une centrale lectro-solaire thermodynamique. Si
l'objectif d'utiliser la centrale est d'amoindrir. le cot du mtre cube, il faut noter que
l'installation des centrales d'nergies renouvelables cote chre et la demande en hydrogne
est encore trs faible.
"-
Faut-il donc produire pour que la demande suive? Ou faut-il attendre que la demande
soit suffisante pour la satisfaire? Dans tous les cas la comptitivit de l'hydrogne devra
passer par la diminution du cot de l'lectricit et de l'amlioration du matriel de stockage
toujours trs cher.
Mais quelqu'un disait que la qualit a un prix mais la non-qualit cote plus chre .
80
On devrait peut tre y rflchir et voir si c'est plus chre de s'offrir une nergie propre ou de
subir la pollution des nergies fossiles et faire un gain d'argent.
Au demeurant, la valorisation de l'oxygne produite lors de l'lectrolyse contribuerait
relever de faon significative le niveau de rentabilit de ce procd.
De ce fait, comme toute uvre scientifique de recherche, ce travail comporte
certainement des insuffisances. Il conviendrait alors de le parfaire en modlisant toutes parties
du systme, en comparant les cots de conception et de ralisation de certains appareils
comme l'lectrolyseur afin de diminuer le cot de l'investissement.
81
BIBLI06RAPHIES
[1] Documentotlon AIR LIQUIDE
Par Cheikh THIAM, ingnieur polytechnicien.
E.S.P. This (2006-2007)
[2] L'hydrogne: les nouvelles technologies de l'nergie.
CLEFS CEA (n 50/51 hiver 2004-2005)
[3] Energie solaire et stockage d'nergie (2
e
dition)
Roger DUMON : MASSON
[4] Technique de L'ingnieur
[5] Capteur solaire
K. EDWARDS / traduit par L. MARLOT
[6] Cours de transfert de chaleur
Professeur: M. Salif GAYE
[7] Kamol L1AMSIRIWATANA, modlisation du fonctionnement d'ensemble de
la centrale lectre-solaire Thmis, thse de docteur-ingnieur, universit de
Perpignan.
[8] Wassila CHEKIROU, Analyse thermique pour l'absorbeur d'un
concentrateur solaire cylindro- parabolique, thse, Dpartement de
physique" Facult.des Sciences, Universit Mentouri Constantine
[7] Cours de gestion de la production et de gestion des projets.
Professeur: M. Ngor SARR
[8] Cours d'nergie renouvelable.
82
P.F.E. Prsentpar M. MbayeSECKet M. MoussaDIOP
Professeur: M. Papa M NDIAYE
E.S.P. This (2006-2007)
[9] Lerayonnement solaire: conversion thermique et application
R. BERNARD G. MENGUY M. SCHWARTZ
(2
e
dition augmente)
[10] LECTRONIQUE" Theodore WILDI"
Avec la collaboration de Gilbert SYBILLE
SITES \'ISITES
[A] hHp:1/perso.wanadoo.frIjean.boubour1ccptfr.htm
[B] http://www.energiesolaire.host.sk/
[Cl http://www.mueller-solartechnik.de/
[D]http://www.bysun.fr
[E]http://www.priceminister.com
[F] http://www.thermoptim.org
[G]http://wwwcgLebay.
[H] http://www.hitachi.fr/products/business.jsp
83
ANNEXES
E.S.P. This (2006 -2007)
ANNEXE 1 : tableaux des flux financiers
ANNEXE 2 : Tableaux des amortissements dgressifs
ANNEXE 3 : Compte d'exploitation prvisionnel
ANNEXE 4 : Tableau des flux financiers en tenant compte de l'oxygne
ANNEXE 5 : Compte d'exploitation prvisionnel en tenant compte de
L'oxygne
84
Tableaux des flux financiers (Annexe 1)
Libell 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017
Investissement
Gnie civil 20000000
Matriel de production 243075791
Matriel et mobilier de bureau 15855000
Imprvus 2000000
Total investissement 280930791
Matires premires 1639000 1639000 1639000 1639000 1639000 1639000 1639000 1639000 1639000 1639000
Personnel 11160000 11160000 11160000 11160000 11160000 11160000 11160000 11160000 11160000 11160000
Mat divers et maintenance 2000000 2000000 2000000 2000000 2000000 2000000 2000000 2000000 2000000 2000000
Total charges d'exploitation 14799000 14799000 14799000 14799000 14799000 14799000 14799000 14799000 14799000 14799000
Recettes 59000000 59000000 59000000 59000000 59000000 59000000 5900000o 59000000 59000000 59000000
Bnfice brut -26209948 -7848411 5533284 15346623 22582268 27942610 31930227 34907655 37138263 38814563
Bnfice net -26209948 -7848411 4149963,29 11509967,6 16936700,9 20956957,2 23947670,6 26180741 27853697,3 29110922,4
Fluxfinancier net 44 201000 44 201000 42817679 40364344 38555433 37215348 36218443 35474086 34916434 34497359
Fluxfinancier net actualis 39465178,57 35236766,6 30476778,1 25652270,4 21877388,1 18854453,3 16383384,3 14327388,5 12591216,2 11107226,4
Cumul des flux fin net actualis 39465178,57 74701945,2 105178723 130830994 152708382 171562835 187946219 202273608 214864824 225972051
a
Annexe 1 : suite
Libell 2018 2019 2020 2021 2022 2023 2024 2025 2026 2027
Investissement 0 0 0 0 0 0 0 0 0 0
Gniecivil 0 0 0 0 0 0 0 0 0 0
Matrielde production O' 0 0 0 0 0 0 0 0 0
Matriel et mobilier de bureau 0 0 0 0 0 0 0 0 0 0
Imprvus 0 0 0 0 0 0 0 0 0 0
Total investissement 0 0 0 0 0 0 0 0 0 0
Personnel 11160000 11160000 11160000 11160000 11160000 11160000 11160000 11160000 11160000 11160000
Total charges d'exploitation 14799000 14799000 14799000 14799000 14799000 14799000 14799000 14799000 14799000 14799000
Recettes 59000000 59000000 59000000 59000000 59000000 59000000 59000000 5 9 0 o o 0 59000000 59000000
Bnfice brut 40078013,08 41033026,4 41756971,7 42307367,9 42727100,1 43048223,9 43294754,8 43484722,4 43631689,2 43745877,2
Bnfice net 30058509,81 30774769,8 31317728,8 31730525,9 32045325,1 32286168 32471066,1 32613541,8 32723766,9 32809407,9
Flux financier net 34181496,73 33942743,4 33761757,1 33624158 33519225 33438944 33377311,3 33329819,4 33293077,7 33264530,7
Flux financier net actualis 9826363,512 8712256,81 7737323,33 6880168,91 6123837,08 5454616,12 4861216,5 4334196,05 3865551,95 3448426,29
Cumul des fluxfin net actualis 235798414,1 244510671 252247994 259128163 265252000 270706616 275567833 279902029 283767581 287216007
b
Tableaux des amortissements dgressifs (Annexe 2)
libell montant dure (ans) taux constant taux dgressif
gniecivil 20000000 20 0,05 0,125
equipde prod 243075791 10 0,1 0,25
mobilier de bureau' 15855000 5 0,2 0,4
mat divers 2000000 5 0,2 0,4
libell 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017
Gniecivil
Annuits 2500000 2187500 1914062,5 1674804,69 1465454,1 1282272,34 1121988,3 981739,759 859022,29 751644,503
Cumul 2500000 4687500 6601562,5 8276367,19 9741821,29 11024093,6 12146081,9 13127821,7 13986844 14738488,5
Valeur rsiduelle 17500000 15312500 13 398438 11 723 633 10258179 8975906 7853918 6872 178 6013156 5261512
Annuits 60768947,75 45576710,8 34182533,11 25636899,8 19227674,9 14420756,2 10815567,1 8111675,34 6083756,5 4562817,38
Cumul 60768947,75 106345659 140528191,7 166165092 185392766 199813523 210629090 218740765 224824521 229387339
Valeur rsiduelle 182306843,3 136730132 102547599,3 76910699,5 57683024,6 43262268,5 32446701,3 24335026 18251269,5 13688452,1
Cotde fonctionnement
Annuits 6342000 3805200 2283120 1369872 821923,2 493153,92 295892,352 177535,411 106521,247 63912,748
Cumul 6342000 10147200 12430320 13800192 14622115,2 15115269,1 15411161,5 15588696,9 15695218,1 15759130,9
Valeur rsiduelle 9513 000 5707800 3424680 2054808 1232885 739731 443839 266303 159782 95869
Matriels divers
Annuits 800000 480000 288000 172800 103680 62208 37324,8 22394,88 13436,928 8062,1568
Cumul 800000 1280000 1568000 1740800 1844480 1906688 1944012,8 1966407,68 1979844,61 1987906,76
Valeur rsiduelle 1200000 720000 432000 259200 155520 93312 55987,2 33592,32 20155,392 12093,2352
Total amortissement 70410947,75 52049410,8 38667715,61 28854376,5 21618732,2 16258390,4 12270772,6 9293345,39 7062736,97 5386436,79
c
Tableau des amortissements
dgressifs (Suite annexes 2)
libell 2018 2019 2020 2021 2022 2023 2024 2025 2026 2027
Gnie civil
Annuits 657688,9404 575477,823 503543,095 440600,208 385525,182 337334,534 295167,718 258271,753 225987,784 197739,311
Cumul 15396177,42 15971655,2 16475198,3 16915798,5 17301323,7 17638658,3 17933826 18192097,7 18418085,5 18615824,8 ..
Valeur rsiduelle 4603822,583 4028344,76 3524801,67 3084201,46 2698676,27 2361341,74 2066174,02 1807902,27 1581914,49 1384175,18
Annuits 3422113,033 2566584,77 1924938,58 1443703,94 1082777,95 812083,464 609062,598 456796,948 342597,711 256948,283
Cumul 232809451,9 235376037 237300975 238744679 239827457 240639541 241248603 241705400 242047998 242304946
Cot de fonctionnement
Annuits 38347,64882 23008,5893 13805,1536 8283,09214 4969,85529 2981,91317 1789,1479 1073,48874 644,093245 386,455947
Cumul 15797478,53 15820487,1 15834292,3 15842575,4 15847545,2 15850527,1 15852316,3 15853389,8 15854033,9 15854420,3
Matriels divers
Annuits 4837,29408 2902,37645 1741,42587 1044,85552 626,913313 376,147988 225,688793 135,413276 81,2479653 48,7487792
Cumul 1992744,059 1995646,44 1997387,86 1998432,72 1999059,63 1999435,78 1999661,47 1999796,88 1999878,13 1999926,88
Total amortissement 4122986,916 3167973,56 2444028,26 1893632,09 1473899,9 1152776,06 906245,152 716277,603 569310,836 455122,799
d
Compte d'exploitation
prvisionnel (Annexe 3)
Dsignation 2008 2009 2010 2011 2012 2103 2014 2015 2016
Matires premires 1639000 1639000. 1639000 1639000 1639000 1 6 ~ 9 0 0 0 1639000 1639000 1639000
Personnel 11160000 11160000 11160000 11160000 11160000 11160000 11160000 11160000 11160000
Divers 2000000 2000000 2000000 2000000 2000000 2000000 2000000 2000000 2000000
Amortissement 70410947,75 52049410,8 38667715,6 28854376,5 21618732,2 16258390,4 12270772,6 9293345,39 7062736,97
Total charge 85209947,75 66848410,8 53466715,6 43653376,5 36417732,2 31057390,4 27069772,6 24092345,4 21861737
Chiffre d'affaire 59000000 59000000 59000000 59000000 59000000 59000000 59000000 59000000 59000000
Valeur ajoute 15049947,75 3311589,19 16693284,4 26506623,5 33742267,8 39102609,6 43090227,4 46067654,6 48298263
Rsultat brut 26209947,75 7848410,81 5533284,39 15346623,5 22582267,8 27942609,6 31930227,4 34907654,6 37138263
Impts 0 0 1383321,1 3836655,87 5645566,96 6985652,4 7982556,86 8726913,65 9284565,76
Rsultat net 26209947,75 7848410,81 4149963,29 11509967,6 16936700,9 20956957,2 23947670,6 26180741 27853697,3
Cash-flow 44201000 44201000 42817678,9 40364344,1 38555433 37215347,6 36218443,1 35474086,3 34916434,2
Valeur actuelle des cash-
flows 39465178,57 35236766,6 30476778,1 25652270,4 21877388,1 18854453,3 16383384,3 14327388,5 12591216,2
Cumul 39465178,57 74701945,2 105178723 130830994 152708382 171562835 187946219 202273608 214864824
valeur actuelle nette pour 9 ans
valeur actuelle des investissements =248 394 533 + 17 855 000/(1,12)"5 =291062 197,5cfa
valeur actuelle nette des cash-flow = 214864824 Fcfa
VAN (9 ans) = - 76 197373,5 Fcfa
e
Tableau des flux financiers en tenant compte de l'oxygne (annexes 4)
Libell 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017
Investissement
Gnie civil 20000000
Matriel de production 243075791
Matriel et mobilier de bureau 15855000
Imprvus 2000000
Total investissement 280930791
Personnel 11160000 11160000 11160000 11160000 11160000 11160000 11160000 11160000 11160000 11160000 1
Total charges d'exploitation 14799000 14799000 14799000 14799000 14799000 14799000 14799000 14799000 14799000 14799000 1
Recettes 59780000 59780000 59780000 59780000 59780000 59780000 59780000 59780000 59780000 59780000
Bnfice brut -25429948 -7068411 6313 284 16126623 23362268 28722610 32710227 35687655 37918263 39594563 4
Bnfice net -25429948 -7068411 4734963,29 12094967,6 17521700,9 21541957,2 24532670,6 26765741 28438697,3 29695922,4 3C
Fluxfinancier net 44981000 44981000 43402679 40949344 39140433 37800 348 36803443 36059086 35501434 35082359 3
Fluxfinancier net actualis 40161607,14 35858577,8 30893169,5 26024048,5 22209332,9 19150832,5 16648008,6 14563660,2 12802173,1 11295580,7 9 ~
Cumul des flux fin net actualis 40161607,14 76020184,9 106913354 132937403 155146736 174297568 190945577 205509237 218311410 229606991 2:
f
Tableau des flux financiers en tenant compte de l'oxygne (Suite annexe 4)
Libell 2018 2019 2020 2021 2022 2023 2024 2025 2026 2027
Investissement
Gnie civil
Matriel de production
Matriel et mobilier de bureau
Imprvus
Total investissement
Personnel 11160000 11160000 11160000 11160000 11160000 11160000 11160000 11160000 11160000 11160000
Total chargesd'exploitation 14799000 14799000 14799000 14799000 14799000 14799000 14799000 14799000 14799000 14799000
Recettes 59780000 59780000 59780000 59780000 59780000 59780000 59780000 59780000 59780000 59780000
Bnfice brut 40858013,08 41813026,4 42536971,7 43087367,9 43507100,1 43828223,9 44074754,8 44264722,4 44411689,2 44525877,2
Bnfice net 30643509,81 31359769,8 31902728,8 32315525,9 32630325,1 32871168 33056066,1 33198541,8 33308766,9 33394407,9
Fluxfinancier net 34766496,73 34527743,4 34346757,1 34209158 34104225 34023944 33962311,3 33914819,4 33878077,7 33849530,7
Fluxfinancier net actualis 9994537,033 8862411,74 7871390,23 6999871,5 6230714,4 5550042,29 4946418,44 4410269,21 3933474,41 3509071,35
Cumul des flux fin net actualis 239601528 248463940 256335330 263335201 269565916 275115958 280062377 284472646 288406120 291915192
g
compte d'exploitation prvisionnel en tenant compte de l'oxygne (Annexe 5)
Dsignation 2008 2009 2010 2011 2012 2103 2014 2015 2016
Matires premires 1639000 1639000 1639000 1639000 1639000 1639000 1639000 1639000 1639000
Personnel 11160000 11160000 11160000 11160000 11160000 11160000 11160000 11160000 11160000
Divers 2000000 2000000 2000000 2000000 2000000 2000000 2000000 2000000 2000000
Amortissement 40161607,14 35858577,8 30893169,5 26024048,5 22209332,9 19150832,5 16648008,6 14563660,2 12802173,1
Total charge 54960607,14 50657577,8 45692i69,5 40823048,5 37008332,9 33949832,5 31447008,6 29362660,2 27601173,1
Chiffre d'affaire 59780000 59780000 59780000 59780000 59780000 59780000 59780000 59780000 59780000
Valeur ajoute 15979392,86 20282422,2 25247830,5 30116951,5 33931667,1 36990167,5 39492991,4 41577339,8 43338826,9
Rsultat brut 4819392,857 9122422,19 14087830,5 18956951,5 22771667,1 25830167,5 28332991,4 30417339,8 32178826,9
Impts 1204848,214 2280605,55 3521957,61 4739237,88 5692916,79 6457541,88 7083247,84 7604334,95 8044706,73
Rsultat net 3614544,643 6841816,65 10565872,8 14217713,6 17078750,4 19372625,6 21249743,5 22813004,9 24134120,2
Cash-flow 43776151,79 42700394,5 41459042,4 40241762,1 39288083,2 38523458,1 37897752,2 37376665 36936293,3
Valeur actuelle des cash-
flows 39085849,81 34040493 29509727,5 25574367,3 22293113,5 19517182,8 17143018,5 15095808,1 13319597,6
Cumul 39085849,81 73126342,8 102636070 128210438 150503551 170020734 187163752 202259561 215579158
valeur actuelle nette pour 9 ans
valeur actuelle des investissements =280930791 + 17 855000/(1,12)/\5 =291062 197,5 cfa
valeur actuelle nette des cash-flows = 215 579 158 Fcfa
VAN (9 ans) = - 75 483 039,5 Fcfa
h

Pfe GM 0498

Transféré par

Informations du documentcliquez pour développer les informations du document

Droits d'auteur :

Formats disponibles

Pfe GM 0498

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

Pfe GM 0498

Transféré par

Droits d'auteur :

Formats disponibles

UNNERSIT CHEIKH ANTA DIOP DE DAKAR

;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.;.; / /////// x..xcccccccccccccoxcccccc :.:.:.:.:.:.:.;.:.. ..xccx ccc ;.:.:.:.:.y ..:.:,o:.:,o:,o:..

);> a =- est.la diffusivit thermique en m

~ Ecoulement en prsence de chicanes

111-3-1-1: Puissance de l'lectrolyseur

La puissance des moteurs pour la compression est donc :

Vous aimerez peut-être aussi