Philippe Dessante HDR

N DORDRE 1515
Mmoire

Optimisation, conception et modlisation des systmes
lectrotechniques

Prsent pour obtenir

LHABILITATION A DIRIGER DES RECHERCHES
DE LUNIVERSITE PARIS-SUD 11

par

Philippe Dessante

Soutenue le 22 juin 2012 devant la Commission dexamen :
M. EICHWALD Olivier, Professeur, Laboratoire Laplace
M HECQUET Michel, Professeur, Ecole Centrale de Lille
M MULTON Bernard, Professeur, ENS de Cachan
M. VANNIER Jean-Claude, Professeur, Dpartement nergie, Suplec
M. BERGER Grard, directeur de recherche, LPGP, Universit Paris Sud
M. ROGER Daniel, Professeur, LSEE - Universit d'Artois - Bthune

Rapporteurs :
M. EICHWALD Olivier, Professeur, Laboratoire Laplace
M HECQUET Michel, Professeur, Ecole Centrale de Lille
M MULTON Bernard, Professeur, ENS de Cachan

Travaux effectus au dpartement nergie de Suplec, Gif-sur-Yvette
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Optimisation, conception et
modlisation des systmes
lectrotechniques
Philippe Dessante
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
2
FicUvr 1 Amazing Stories aot 1926
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Premire partie
Fiche synthtique
3
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Chapitre 1
CV scientifque
1.1 tat civil
Philippe Dessante
N le 28 dcembre 1971, 40 ans.
Clibataire, vie maritale, un enfant.
Professeur adjoint au dpartement nergie,
quipe daccueil E3S EA4454, Suplec.
Ml : [email protected]
Adresse : Dpartement nergie, Suplec
3 rue Joliot Curie
Plateau de Moulon
91192 Gif-sur-Yvette
1.2 Titres universitaires
2000 : Thse de doctorat de lUniversit de Versailles - Saint-Quentin-en-
Yvelines ;
1996 : DEA de Modlisation mathmatique, simulation et applications la
physique, Universit de Versailles - Saint-Quentin-en-Yvelines ;
1995 : Matrise dingnierie mathmatique, Universit de Montpellier II ;
1995 : Licence de mathmatiques, Universit de Montpellier II ;
1993 : DEUG A, Universit de Montpellier II.
5
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
6 Cu~vi1vr 1 : CV scientifque
1.3 Parcours professionnel
1996-2000 : thse de Doctorat au laboratoire de physique des gaz et des
plasmas ;
depuis 2001 : professeur assistant puis professeur adjoint au dpartement
nergie de Suplec (EA 4454).
1.4 Activits denseignement
Responsabilits denseignement :
charg du cours mthodes numriques et optimisation en deuxime anne
de Suplec (33h eq. TD) ;
charg de cours en troisime anne de Suplec en optimisation et me-
thodes numriques (28h eq. TD) ;
TD dlectrotechnique en premire anne de Suplec (6h eq. TD) ;
TD de mthodes numriques et optimisation en deuxime anne de Sup-
lec (6h eq. TD) ;
TP dlectrotechnique en premire anne de Suplec (50h eq. TD) ;
encadrement de projets dtudes dlves ingnieurs Suplec (50h eq. TD) ;
cours en formation continue danalyse numrique, doptimisation,
dlectrotechnique, et de physique des dcharges (15 eq. TD) ;
co-responsable du mastre spcialis Ingnieur danaire pour les nou-
veaux marchs de llectricit (IANME).
1.5 Activits lies ladministration
Membre du comit de pilotage du projet nergie, Suplec.
1.6 Activits lies la recherche
1.6.1 Prix
"PLASMA PHYSICS INNOVATION PRIZE" dcern par la Socit Euro-
penne de Physique, "for breakthrough developments and applications of
basic plasma physics tools to address environmental concerns" en colla-
boration avec Emmanuel Marode, Christophe Laux, Djamal Djermoune,
Christian Deniset, Pierre Sgur, Franois Bastien, et Anne Bourdon.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
1.6 Activits lies la recherche 7
1.6.2 Organisations de congrs scientifques
Organisation de congrs scientihques :
lectrotechnique du futur (2003) : EF2003, Suplec 2003 ;
Journes du club EEA, Suplec 2006 ;
Journe dconomie exprimentale, Suplec 2007 ;
sixime Confrence de la Socit Franaise dlectrostatique, SFE 2008.
1.6.3 Comits de lecture
Participation des comits de lecture :
revues darticles pour la revue IET Generation, Transmission and Distribu-
tion ;
revues darticles pour le congrs Power System Computation Conference :
PSCC 2005 ;
revues darticles pour la confrence de la socit franaise dlectrosta-
tique : SFE 2008 ;
revues darticles pour la revue IEEE TDEI Electrostatics ;
revues darticles pour la revue Journal of Electrostatics.
1.6.4 Participation des jurys de thse
Examinateur de la thse dAssia Djellalli : Optimisation technico - cono-
mique dun rseau dnergie lectrique dans un environnement drgul.
Jury : J. Robert, F. Bernot, A. Jaafari, M. Meunier, Ph. Dessante.
Examinateur de la thse de Marcelo Saguan : Analyse conomique des
architectures de march lectrique. Application au Market Design du temps
rel .
Jury : R. Belmans, C. von Hirschhausen, S. Saussier, Ph. Dessante, M. Mas-
soni, J.-M. Glachant.
Examinateur de la thse de Vincent Rious : Le dveloppement du rseau de
transport dans un systme lectrique libralis, un problme de coordination
avec la production.
Jury : R. Green, J. Doucet, S. Saussier, Ph. Dessante, J.-M. Coulondre, . Le
Boulanger, J.-M. Glachant.
Examinateur de la thse de Florent Maupas : Analyse de limpact cono-
mique de lala olien sur la gestion de lquilibre dun systme lectrique.
Jury : J.-M. Glachant, J. Usaola, Ch. Weber, E. Neau, U. Stridbaeck, Ph.
Dessante.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Examinateur de la thse de Julien Capeillere : Modlisations numriques
bidimensionnelles du transport des particules et photons dans des gaz ioniss.
Application au xnon sous excitation laser et la propagation dune dcharge
monoflamentaire barrires dilectriques dans lazote.
Jury : O. Eichwald, J-F. Loiseau, J-M. Pouvesle, Ph. Dessante, M.C. Bordage,
N. Sewraj.
Examinateur de la thse dHerman Bayem : Apport des mthodes probabi-
listes pour linsertion de la production olienne dans les systmes lectriques
.
Jury : J. Usaola, V. Miranda, G. Kariniotakis, Franoise Lamnabhi, Ph. Des-
sante, J.-C. Vannier.
Examinateur de la thse de Michelle Nsoumbi : tude des mcanismes
dinfammation dun matriau isolant en prsence dun point chaud dorigine
lectrique.
Jury : N. Ben Jemaa, S. Berthelot, S. Bourbigot, Ph. Dessante, D. Diallo, Z.
Khatir, P. Lagonotte, Ph. Test.
1.6.5 Partenariats industriels et projets de ples de comptiti-
vit
Projet Comette, (2010-) Ph. Dessante, P. Vidal. Modlisations et calculs de
performances pour la propulsion lectrique automobile.
Projet CSDL Complex Systems Design Lab, (2010-) Ph. Dessante, P. Vidal,
G. Sandou. Rduction de modles complexes et optimisation.
Projet Preface, (2009-) Ph. Dessante, P. Vidal, B. Lorcet, M. Hennebel. Mo-
dlisation de type circuit du passage de la foudre dans les composants dun
avion.
Leroy Somer (2009-) Ph. Dessante, P. Vidal, Jean-Claude Vannier. Optimi-
sation de systmes de motorisation.
Projet Sefora, (2008-2010) Ph. Dessante, P. Vidal. tude doptimisation
dun actionneur pour application freinage. Optimisation dune chane de
motorisation (actionneur, rducteur, liaison hlicodale).
Projet O2M, (2008-2010), Ph. Dessante, P. Lefranc, S. MRad. tude sur le
dveloppement dun modle dinusif dun bras donduleur dun alterno-
dmarreur.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
1.6 Activits lies la recherche 9
PSA (2007-2010) Ph. Dessante, E. Odic, Ph. Test, R. Meyer, M. Nsoumbi.
tude des mcanismes dinfammation dun matriau isolant en prsence
dun point chaud dorigine lectrique. Mesures exprimentales dchauf-
fement sur des pistes lectriques, modlisation et simulation des phno-
mnes.
Virax (2007) Ph. Dessante, P. Vidal. Optimisation du systme de motorisa-
tion pour un outil portatif. Modlisation analytique du moteur, du rduc-
teur picyclodal et de la batterie. Optimisation de la masse du systme
complet sous les contraintes de fonctionnement/performance.
RTE (2006-2008) : Ph. Dessante, F. Maupas. tude du mcanisme dajuste-
ment en prsence de lala olien. Dveloppement dun modle du mca-
nisme dajustement. Optimisation de lajustement en fonction du besoin
et des moyens mis disposition. Application de ce modle en prsence de
lala olien.
Sodielec (2006) : Ph. Dessante, M. Petit. tude dun algorithme destima-
tion dtat sur un rseau de distribution. Mise au point thorique de lal-
gorithme. Calcul de load fow. Application de lestimateur dtat au cas du
load fow sur un rseau de distribution lectrique.
RTE (2006-2007) : Ph. Dessante, V. Rious. tude sur les investissements
sur le rseau de transport dans un environnement de march de llec-
tricit. Coordination production-transport. Modlisation des acteurs, op-
timisations des dcisions dinvestissement.
Renault (2006-2009) : Ph. Dessante, M. Petit. tude dun rseau de bord
lectrique automobile, tude et conception dun modle dalternateur. ta-
blissement des paramtres du modle par mesures exprimentales. Mod-
lisation du rseau en rgime transitoire.
CRE (2005-2007) : Ph. Dessante, M. Saguan. tude des architectures de
march de llectricit. Modlisation analytique des dinrents designs de
march. Optimisation du choix des acteurs. Analyse conomique des r-
sultats.
EDF (2005-2009) : Ph. Dessante, S. Plumel, M. Petit, H. Bayem. tude sur
lopportunit dutiliser des mthodes probabilistes pour les tudes dint-
gration de lolien dans les rseaux lectriques, et en particulier sur les
tudes de raccordement. tude bibliographique. Modles probabilistes de
vent et de production olienne. Consquences de linsertion de lolien
dans les rseaux.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Basta (2005-2006) : Ph. Dessante, P. Vidal, M. Bensetti. tude des perfor-
mances et optimisation dune gnratrice embarque. Simulations par l-
ments hnis en trois dimensions. Calculs de couples. Vrihcations et tudes
des paramtres par des mesures exprimentales.
Sodern (2004-2005) : Ph. Dessante, J-C. Vannier, B. Bonafos. Dhnition des
caractristiques dun blindage magntique plusieurs tages. Analyse du
comportement dun blindage contre le champ magntique terrestre. Mise
au point dune optimisation dun systme trois niveaux dcran en re-
cherchant la masse minimale pour une application sur satellite. Simulation
de linfuence mutuelle de deux blindages voisins.
PSA (2003-2004) : Ph. Dessante, B. Lorcet. tude de cosimulation dun r-
seau de bord lectrique automobile, simulation sous Simulink et Saber.
Modlisation dun rseau de bord, tude de cas critiques.
1.7 Encadrements
1.7.1 Thses
Marcelo Saguan : 50% Analyse conomique des architectures de march
lectrique. Application au Market Design du temps rel , thse de lUniver-
sit Paris-Sud 11, soutenue le 28 avril 2007, partenariat Facult Jean Mon-
net/CRE, directeur de thse : Jean-Michel Glachant (50%). Actuellement en
poste de consultant chez Microeconomix.
Jury : Ronnie Belmans, Christian von Hirschhausen, Stphane Saussier,
Philippe Dessante, Michel Massoni, Jean-Michel Glachant.
[50][47][46][45].
Vincent Rious : 50% Le dveloppement du rseau de transport dans un sys-
tme lectrique libralis, un problme de coordination avec la production,
partenariat Facult Jean Monnet/RTE, thse de lUniversit Paris-Sud, sou-
tenue le 30 octobre 2007, directeur de thse : Jean-Michel Glachant (50%).
Actuellement en poste de consultant chez Microeconomix.
Jury : Richard Green, Joseph Doucet, Stphane Saussier, Philippe Des-
sante, Jean-Marc Coulondre, ric Le Boulanger, Jean-Michel Glachant.
[38][11][30][9] [31][29][22][32] [23][34][28][33][19] [24]
Florent Maupas : 50% Analyse de limpact conomique de lala olien sur
la gestion de lquilibre dun systme lectrique, partenariat Facult Jean-
Monnet/RTE, thse de lUniversit Paris-Sud, soutenue le 2 juillet 2008,
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
1.7 Encadrements 11
directeur de thse : Jean-Michel Glachant (50%). Actuellement ingnieur
chez RTE.
Jury : Jean-Michel Glachant, Julio Usaola, Christoph Weber, Emmanuel
Neau, Ulrik Stridbaeck, Philippe Dessante. [40]
Herman Bayem : 40% Apport des mthodes probabilistes pour linsertion
de la production olienne dans les systmes lectriques , partenariat EDF,
thse de lUniversit Paris-Sud soutenue le 23 novembre 2009, directeur
de thse : Jean-Claude Vannier (10%), co-encadrant : Marc Petit (50%). Ac-
tuellement ingnieur chez EDF.
Jury : Julio Usaola, Vladimiro Miranda, Georges Kariniotakis, Franoise
Lamnabhi, Philippe Dessante, Jean-Claude Vannier. [44] [37] [36]
Michelle Nsoumbi : 33% tude des mcanismes dinfammation dun ma-
triau isolant en prsence dun point chaud dorigine lectrique, partenariat
LGEP/PSA, thse de lUniversit Paris-Sud, soutenue le 28 septembre 2010,
directeur de thse : Philippe Test (33%), co-encadrant : Emmanuel Odic
(33%). Actuellement ingnieur chez PSA.
Jury : Noureddine Ben Jemaa, Sylvie Berthelot, Serge Bourbigot, Philippe
Dessante, Demba Diallo, Zoubir Khatir, Patrick Lagonotte, Philippe Test.
Henri Borsenberger : 50% Modlisation et optimisation robuste de sys-
tmes complexes : application aux rseaux nergtiques, thse de lUniver-
sit Paris-Sud/Suplec, hnance par la Fondation Suplec en cours, direc-
teur de thse : Philippe Dessante, co-encadrant : Guillaume Sandou (50%).
[21] [59][8]
Romaric Lanfried : 33% Contribution ltude de la transition dcharge lu-
minescente - arc lectrique dans lair et largon au voisinage de la pression at-
mosphrique, thse de lUniversit Paris-Sud/Suplec soutenue en octobre
2011, hnancement Carnot, directeur de thse : Emmanuel Odic (33%), co-
encadrant : Philippe Test (33%).[16]. Actuellement en poste dATER.
Jury : Francisco Alves, Laurent Chemartin, Jean-Marie Cormier, Alain
Gleizes, Michael J. Kirkpatrick, Christophe Laux, Olivier Lesaint, Emma-
nuel Odic, Philippe Test.
Benjamin Dagus : 50% Optimisation de systmes dentranement, parte-
nariat Leroy-Somer, thse de lUniversit Paris-Sud/Suplec , directeur de
thse : Jean-Claude Vannier (50%), soutenance prvue en octobre 2012.
[18]
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Florent Sainct : 30% tude de la ractivit de dcharges lectriques nanose-
conde la pression atmosphrique en milieu riche en hydrocarbures. valua-
tion de lefcacit nergtique de production dhydrogne, partenariat cole
Centrale de Paris, laboratoire EM2C, directeur de thse : Christophe Laux
(40%), co-encadrant Emmanuel Odic (30%), soutenance prvue en 2013.
Khaled Almaksour : 33% tude de lmission de courant noir par une surface
soumise un champ lectrique intense et de ses mcanismes de transition
larc, partenariat laboratoire LGEP, thse de lUniversit Paris-Sud/Suplec
, directeur de thse : Philippe Test (33%), co-encadrant Emmanuel Odic
(33%), soutenance prvue en 2013.
Dany Prieto : 30% Modlisation et optimisation des machines synchro-
rluctantes aimants permanents et de leur lectronique Application aux
vhicules lectriques et aux compresseurs, partenariat Leroy-Somer, thse de
lUniversit Paris-Sud/Suplec , directeur de thse : Jean-Claude Vannier
(40%), co-encadrant Pierre Vidal (30%), soutenance prvue en 2014.
Ivan Kravtzon : 40% Design et contrle dune nouvelle confguration de
groupe lectrogne associe un dispositif de stockage dnergie permettant
de minimiser la consommation et lempreinte environnementale, partena-
riat Leroy-Somer, thse de lUniversit Paris-Sud/Suplec , directeur de
thse : Daniel Sandarnac (20%), co-encadrant Pierre Lefranc (40%), soute-
nance prvue en 2014.
Participation lencadrement des thses de Miguel Lopez, Dario Morales,
Xavier Jannot et Christophe Gutfrind.[12][41][27] [17]
1.7.2 Stagiaires post-doc
Co-encadrement avec Claude Marchand et Jean-Claude Vannier de Moha-
med Bensetti, stagiaire de post-doctorat pour une tude sur lanalyse et
loptimisation des performances dune gnratrice embarque. [48]
Co-encadrement avec Pierre Lefranc de Sabrine MRad, stagiaire de post-
doctorat pour une tude sur le dveloppement dun modle dinusif dun
bras donduleur dun alterno-dmarreur. [26]
1.7.3 Stagiaires Master 2
Jesus Periscal, Mise en place du logiciel de base pour la simulation dun
march lectrique dans le cadre de lconomie exprimentale, stage de hn
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
1.8 Synthse scientifque 13
dtude de master, 2005.
Fabien Petit, conomie exprimentale, outil pour ltude de lefcacit du
market coupling, stage de hn dtude dingnieur Suplec, 2006. [39][33]
Fabien Petit, Efcacit des mthodes dallocation implicite des capacits de
transport du rseau lectrique dans un rseau fortement maill, stage de
master EGIR, 2007.
Pablo de la Cruz Martinez, Estimation dtat sur un rseau de distribution
lectrique, stage de master, 2007.
Ridga Ben Abdallah, Simulation du dplacement des charges lectriques au
sein dune dcharge lectrique laide dune mthode ELLAM, stage de mas-
ter, 2010.
1.8 Synthse scientifque
Aprs ma thse soutenue en octobre 2000 enectue sur la modlisation dune
dcharge lectrique de type streamer, je suis entr au dpartement nergie de
Suplec. Jy ai rejoint une quipe travaillant sur les machines lectriques sous la
direction de Monsieur Jean-Claude Vannier. Mes travaux portent sur la modli-
sation et loptimisation des systmes lectrotechniques.
Jai tout dabord mis proht mes connaissances sur la modlisation du champ
lectrique et de la propagation des dcharges pour ltude de la tenue en tension
en rgime nominal et en coup de foudre de bagues pour lalimentation dun pod.
Cette tude a ensuite porte sur des calculs de tempratures enectus laide de
circuits lectriques quivalents et par lments hnis.
Par la suite, lobjectif de mes recherches sest ax sur deux points : la mod-
lisation numrique et loptimisation. Jai donc dvelopp des modles physiques
et numriques par lments hnis appliqus la modlisation des machines lec-
triques ou phnomnes lectrotechniques. Ces modles ont aliment, par le biais
de simplihcations, des algorithmes doptimisation pour la conception de ma-
chines. Le dveloppement de ces algorithmes a ensuite t appliqu loptimi-
sation technico-conomique des grands rseaux de transport dlectricit. Dans
les paragraphes suivants, je dvelopperai une prsentation succincte des travaux
raliss.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Avec la collaboration de Bernard Bonnafos, jai tudi les systmes de blin-
dage lectromagntiques en vue dune utilisation dans une horloge atomique
place dans un satellite. Le fonctionnement correct de lhorloge ncessite davoir
un champ magntique le plus faible et le plus homogne possible sur laxe de ces
blindages. Le systme tudi comportait trois blindages imbriqus. La gomtrie
du problme ne permettant pas lutilisation de symtrie, le calcul du champ ma-
gntique rsultant sur laxe de lhorloge a d tre men grce lutilisation dun
code de calcul par lments hnis en trois dimensions. Lattnuation du champ
magntique calcul ntant pas assez importante, il a fallu introduire une com-
pensation active[6][54]. Cette dernire tait compose dun certain nombre de
bobinages cocylindriques avec le blindage passif. En fonction du champ rsi-
duel sans compensation active (mesure par un capteur), une optimisation sur
le nombre de bobines et sur les densits de courant appliquer dans celles-ci a
t introduite. Ceci ma permis destimer la meilleure distribution des courants
dans les bobines ahn de minimiser le champ magntique extrieur sur laxe de
lenceinte de lhorloge.
Jai par la suite appliqu ces comptences acquises lors de cette tude lop-
timisation des systmes de motorisation. Dans de nombreuses applications in-
dustrielles (spatiales, aronautiques, automobiles, outillage portatif), il est n-
cessaire de minimiser la masse du systme ou les pertes, tout en gardant les
performances constantes. La premire tape a consist modliser analytique-
ment les lments (moteur, gnrateur, batterie, convertisseur, rducteur de vi-
tesse et transformateur de mouvement ventuel), partir de simplihcations des
modles numriques lments hnis, pour relier les performances et les dimen-
sions. La seconde tape a consist mettre au point un algorithme doptimisation
relativement ouvert ahn de prendre en compte dinrentes contraintes et din-
rents objectifs en intgrant les modles prcdemment mis au point et valids.
Ainsi, diverses tudes ont port sur la minimisation de la masse ou le volume du
systme[52][51][49][48][43]. Une premire optimisation multiobjectif (minimi-
sation du volume et des pertes) a t mene sur une pompe carburant pour un
satellite[43].
Pour valider cette approche, une tude enectue dans le contexte dun pro-
jet de ple de comptitivit a port sur loptimisation globale et multiobjectif
dun systme de motorisation et dentranement pour une application aronau-
tique. Le systme complexe compos dlments mcaniques et lectriques (r-
ducteur, transformateur de mouvement et moteur) a t modlis par des fonc-
tions simples permettant de dvelopper des algorithmes doptimisation perfor-
mants. La double optimisation (simple et par front de Pareto) des pertes et de la
longueur du systme a bien montr limportance de lagrgation des composants
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
dans loptimisation.
Ces travaux se poursuivent actuellement au travers de la thse de Benja-
min Dagus portant sur loptimisation globale dun systme de motorisation.
Les recherches sorientent dune part vers la paralllisation des codes de calcul
doptimisation et dautre part vers lincorporation de modles par lments hnis.
Ces comptences en optimisation mont permis de participer aux activits du
groupe rseau du dpartement nergie de Suplec, dans le cadre de loptimisa-
tion technico-conomique des grands rseaux dnergie. Jai tout dabord la-
bor conjointement avec Sophie Plumel des algorithmes doptimisation du plan
de production pour la gestion des congestions au sein des marchs de llec-
tricit. Une premire optimisation est ralise sous la forme dun quilibre de
Nash-Cournot pour dterminer les prix et les quantits changes sur le mar-
ch de llectricit. Une deuxime optimisation est alors conduite pour sassurer
des contraintes physiques du rseau lectrique. Ces contraintes rendent nces-
saire la modlisation complte et non linaire du rseau laide dun calcul de
load-fow[14][12][57].
Jai ensuite travaill llaboration des modles de designs des marchs de
llectricit en partenariat avec Jean Michel Glachant du groupe Rseau Jean-
Monnet de lUniversit de Paris-Sud. Lors de cette tude, il a t dvelopp un
modle analytique probabiliste du march du temps rel. Ce modle a permis
ltude plus approfondie de linfuence de la dynamique de la succession des
marchs en environnement drgul. Une analyse conomique a ensuite t me-
ne avec Marcelo Saguan dans le cadre de la thse quil a soutenue en avril
2007[50][47][46][45].
Le partenariat avec le groupe rseau Jean-Monnet ma conduit coencadrer
une thse sur les problmes dinvestissements sur le rseau lorsque sa gestion
est dconnecte de la production. Il sagissait de faire une suite doptimisations
qui permettait dobtenir une tude sur les possibilits dinvestissement pour le
gestionnaire du rseau de transport (GRT). En enet, le temps de fabrication dune
ligne de transport dlectricit est en gnral suprieur celui de la construction
dune centrale quelle doit relier. Le but de ce travail conduit dans le cadre de la
thse de Vincent Rious a t dtudier entre autres les signaux de coordination
possibles qui peuvent aider le GRT dcider de ses investissements en collabo-
ration avec les producteurs dlectricit[38].
Dans le mme cadre de coopration, jai co-encadr la thse de Florent Mau-
pas sur limpact de lala olien sur le mcanisme dajustement du rseau lec-
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
trique. La premire tape t la ralisation dun modle du mcanisme dajus-
tement. Il sagit dune srie doptimisations de la gestion dynamique du parc
de production pour conserver lquilibre onre-demande tout au long de la jour-
ne. Cette gestion se fait en fonction des groupes de production disponibles et
de leur contraintes physiques et en fonction de lala de consommation. Aprs
avoir valid ce modle de gestion du mcanisme dajustement, nous avons tu-
di limpact de lala olien sur celui-ci. Lanalyse des rsultats nous a permis
de montrer les avantages et les inconvnients de lincorporation de lolien dans
les marchs et lquilibrage de llectricit.[40].
En parallle de ce travail, la thse de Herman Bayem a port sur lopportunit
dutiliser des mthodes probabilistes pour les tudes dintgration de lolien et
du photovoltaque dans les rseaux lectriques, et en particulier sur les tudes
de raccordement. Il sagit, partir de modles probabilistes ou des relevs de
vents et densoleillement rels, de remonter une production lectrique proba-
biliste. On obtient alors une probabilit de production pour une olienne ou un
panneau que lon peut agrger en tenant compte des corrlations entre chaque
olienne pour obtenir celle dune ferme ou dune rgion. Avec ces donnes,
nous pouvons dterminer limpact (probabilit de congestions par exemple) de
la production alatoire sur le rseau lectrique[44]. Cette tude a fait lobjet de
deux applications sur des rseaux insulaires.[36]
Dans la continuit de ces travaux sur les alas et les incertitudes, il est apparu
ncessaire de dvelopper des algorithmes doptimisation robustes tenant compte
de ces incertitudes sur les fonctions objectifs ou les contraintes. Un travail explo-
ratoire a donc t men lors de la thse dHenri Borsenberger en collaboration
avec Guillaume Sandou du dpartement Automatique de Suplec. La gestion
dun rseau de production, distribution ou stockage dnergie est devenue un
enjeu tant conomique que technique, devant respecter des contraintes environ-
nementales et lgislatives. Il est donc apparu ncessaire doptimiser cette gestion
ahn de maximiser le proht de fonctionnement du rseau tout en respectant les
contraintes imposes. Cette tche sest complexihe du fait de la diversihcation
des sources dnergie qui peuvent coexister dans un mme rseau (lectricit,
gaz, houl, vapeur...). Par ailleurs, la modlisation de ce type de systmes doit in-
tgrer les incertitudes lies la mconnaissance ou la simplihcation du modle
(en vue de son optimisation) ou des incertitudes provenant de la nature prvi-
sionnelle de la planihcation du fonctionnement du systme (demande relle des
consommateurs, fuctuations conomiques). La robustesse de cette gestion opti-
male est ds lors imprative, de sorte que les performances du systme puissent
tre garanties en dpit des dinrentes sources dincertitudes. Lors de ce travail
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
de thse, trois types dincertitudes ont t pris en compte : sur la demande des
consommateurs, sur la capacit de production maximale, et sur les cots de pro-
duction. Dinrents algorithmes doptimisation ont t labors, en sappuyant
sur les mthodes existantes dans le cas des mthodes dterministes (program-
mation dynamique et relaxation lagrangienne).
En parallle de ces travaux sur loptimisation des systmes lectrotechniques,
jai continu, en partenariat avec Emmanuel Odic, apporter ma connais-
sances sur la physique des plasmas et des modles associs aux activits du
dpartement[55]. Nous avons ainsi dvelopp des modles de calcul et enectu
des simulations de tempratures au sein de la colonne de plasma qui se forme
aprs le passage de la dcharge de type streamer, en vue de modliser lactivit
chimique[13][53].
En partenariat avec Philippe Test du laboratoire LGEP nous avons men une
tude des dfaillances svres susceptibles de conduire un amorage de feu
sur circuit imprim lors de la thse de Michle Nsoumbi. Ce travail bas sur
des mesures exprimentales et des simulations lectrothermiques ma permis de
mettre proht mon exprience en modlisation par lments hnis acquise lors
de travaux sur la modlisation des machines lectriques. Les expriences ont t
conduites jusqu rupture dune piste de circuit imprim aboutissant ou non un
dpart de feu. Le courant de fuite circulant dans le substrat sest rvl tre un
paramtre pertinent pour valuer la dgradation du circuit imprim. Un modle
numrique lectrothermique en trois dimensions dvelopp sur cette structure
a permis de dcrire et de valider la dominance dun processus thermique jus-
qu la rupture de la piste. Un enregistrement dtaill des paramtres lectriques
synchroniss une cinmatographie rapide au moment de la rupture a permis
didentiher un mcanisme impliquant un dcollement de piste au niveau du d-
faut li un emballement de la temprature et du courant de fuite localis. Ce
mcanisme a t associ la prsence dun canal chaud au travers du substrat
conduisant un chaunement important et localis du matriau par enet Joule
susceptible de conduire linfammation, ainsi quun entretien de chaunage du
substrat aprs rupture de la piste contribuant entretenir et propager le feu,
notamment par mission de composs infammables.
Cette collaboration, nous a amen mettre en commun nos comptences,
tant au niveau des domaine de recherche que des mthodes utilises. En enet,
le Laboratoire de Gnie Electrique de Paris a dvelopp une comptence dans
le domaine des arcs lectriques (exprience et modlisation) et le dpartement
Energie de Suplec a dvelopp une comptence dans ltude exprimentale et
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
par simulation de dcharges lectriques pr-disruptives dans plusieurs milieux
gazeux. Les dcharges tudies par les deux groupes se situent dans les domaines
dinrents, mais complmentaires, des plasmas thermiques et des plasmas froids
(hors quilibre thermodynamique). Ltude de la transition dune dcharge lumi-
nescente (forte tension, faible intensit du courant) en un arc lectrique (faible
tension, forte intensit du courant) prsente un intrt acadmique majeur tant
du point de la physique des dcharges (incluant larc lectrique) que du point
de vue de linteraction plasma / surface. Cest pourquoi nous avons dbut une
thse sur ce sujet de la transition dcharge glowarc. Les phnomnes physiques
qui rgissent la physique de la dcharge lectrique sont sensiblement dinrents
de ceux qui gouvernent la physique de larc lectrique et jusqu prsent les
mcanismes conduisant la transition ne sont pas bien connus. Ltude de la
transition dun mode lautre prsente donc un intrt scientihque particulier.
Cette thse, toujours en cours, nous a permis de dcrire les phnomnes prsents
lors de la transition glow-arc grce des mesures lectriques et visuelles laide
dune camra rapide. Cette description va par la suite nous permettre de rali-
ser un modle physique des phnomnes qui devrait conduire une simulation
numrique de validation et dinterprtation.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
FicUvr 1.1 Amazing Stories aot 1926
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Deuxime partie
Synthse des travaux de recherche
21
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Chapitre 2
Introduction
Llectrotechnique est une science relativement ancienne, mais toujours
dactualit. Que ce soit le transport dlectricit, les machines lectriques ou
llectrostatique, ce sont des domaines de recherche tudis depuis longtemps,
pourtant, on observe une augmentation de lintrt et de la demande socitale
pour la recherche en lectrotechnique. La consommation dlectricit mondiale
et en France crot de plusieurs pour-cent par an et les rcents dveloppement
de la drgulation du march de llectricit imposent une gestion des cots et
une optimisation des moyens de productions. Ces problmatiques sont compl-
tes par lintroduction de plus en plus massive des nergies renouvelables qui
apportent de lincertitude supplmentaire dans la rsolution des optimisations
technico-conomiques.
Cette entre des nergies renouvelables est due une forte demande soci-
tale de rduction des missions de gaz enet de serre. De la mme manire,
le monde des transports assiste une lectrihcation accrue des organes allant
jusquau remplacement des moteurs thermiques par des moteurs lectriques.
Des contraintes dencombrement et de poids entranent une recherche impor-
tante dans loptimisation des systmes de motorisation lectriques. Grce aux
connaissances acquises dans loptimisation des systmes dnergies nous avons
dvelopp des algorithmes et des mthodes doptimisation adapts aux systmes
de motorisation.
Ce document sarticule autour de trois parties principales. Dans la premire,
nous nous intresserons la problmatique de loptimisation et de la modlisa-
tion technico-conomique des rseaux lectriques. Nous y prsenterons les col-
laborations avec le monde conomique par le biais de co-encadrements de thse
avec le groupe rseau Jean-Monet. Larrive massive des nergies renouvelables
nous a par la suite pousse nous intresser lincorporation des incertitudes
dans les modles et les optimisations.
Dans une deuxime partie, nous traiterons de loptimisation des systmes de
23
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
24 Cu~vi1vr 2 : Introduction
motorisation. Aprs une introduction de la problmatique et la prsentation de
la mthodologie, nous prsenterons deux exemples plus dvelopps : une optimi-
sation industrielle dune gamme de moteurs lectriques, o nous prsenterons
de nouvelles techniques pour rduire les cots de fabrication et dutilisation.
Loptimisation de la chane de motorisation ne peut se rduire uniquement au
traitement du moteur, cest pourquoi notre approche doptimisation systme se
base sur la prise en compte complte du systme de motorisation : moteur, r-
ducteur de vitesse, convertisseurs de mouvement et de puissance, source dner-
gie. . .Pour cela, il faut sassurer que les modles pris pour chaque composant
puissent tre traits par loptimisation. Nous prsenterons donc une modlisa-
tion et une premire optimisation dun convertisseur de puissance en vue de son
incorporation dans loptimisation globale de la chane de motorisation.
Pour faire de loptimisation, il faut pouvoir disposer de modles de com-
portement, que ce soit pour les objectifs de la minimisation comme pour les
contraintes. En enet optimiser cest bien souvent se positionner en limite de ha-
bilit pour le systme, il faut alors bien comprendre les mcanismes physiques et
assurer sa robustesse. Ce point fera lobjet de la dernire partie. Dans un premier
exemple nous traiterons de linfammation des cartes lectroniques soumises
un stress important o nous prsenterons des rsultats exprimentaux coupls
une modlisation lectrothermique des mcanismes conduisant linfamma-
tion des circuits imprims. Enhn et toujours en rapport avec la habilit nous
exposerons dans la dernire partie de ce chapitre, les rcents avancements des
modles en physique des dcharges par le biais dun travail sur la comprhen-
sion des mcanismes de transition entre la dcharge glow luminescente et larc
lectrique.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
25
FicUvr 2.1 Amazing Stories septembre 1926
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
26 Cu~vi1vr 2 : Introduction
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Chapitre 3
Optimisation et modlisation
technico-conomique des rseaux
lectriques
3.1 Introduction
Depuis quelques dizaines dannes on assiste une rforme importante de
lindustrie lectrique. Lintroduction de la concurrence sinscrit dans la ligne
des mouvements de rformes des secteurs conomiques relevant traditionnelle-
ment de monopoles. Le dpartement nergie une tradition de recherche im-
portante dans le secteur des rseaux lectriques. Il nous a paru important de
suivre ces rformes, de les comprendre et mme de les devancer. Lintroduction
de lconomie, domaine tranger aux comptences axes techniques du dpar-
tement nous a conduit nous associer au laboratoire GRJM de la facult Jean-
Monnet de Sceaux. Ce partenariat a donn lieu de nombreuses collaborations
o nous apportions nos comptences techniques du monde de lnergie compl-
tes par le savoir-faire conomique du groupe de recherche Jean-Monnet. Plu-
sieurs thses ont t soutenues, toutes au sein du projet pluridisciplinaire ner-
gie de Suplec et co-encadres par Jean-Michel Glachant professeur en science
conomique dans le groupe rseau Jean-Monnet de luniversit Paris 11. Jai moi-
mme particip lencadrement de trois dentre elles, celles de Marcelo Saguan,
Vincent Rious et Florent Maupas. Tout au cours de cette collaboration, il nous
est apparu ncessaire de ne pas uniquement se focaliser sur les aspects tech-
niques des rseaux lectriques, mais de marier ceux-ci la nouvelle approche
conomique qui modihe en profondeur les problmes prexistants.
27
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
28
Cu~vi1vr 3 : Optimisation et modlisation technico-conomique des rseaux
lectriques
3.2 Designs de marchs
Nos travaux se sont tout dabord orients vers lanalyse, la description et
la modlisation des formes que pouvait prendre le de march de llectricit au
travers de la thse de Marcelo Saguan. Cette recherche sinscrit dans lanalyse
conomique des architectures de march lectrique. La construction dune ar-
chitecture de march est une condition ncessaire pour la cration dun march
de gros dlectricit. Le manque de matrise rationnelle dans le choix des op-
tions optimales de design pousse ainsi une recherche plus approfondie sur les
dinrents designs, prenant en compte les spcihcits du bien lectricit. Cette
recherche sest appuye sur un cadre danalyse modulaire permettant de sparer
les problmes du market design en autant de modules distincts. Lanalyse modu-
laire a montr le rle cl du module du temps rel qui constitue le noyau central
de toute architecture de march lectrique et o tous les changes physiques
entre acteurs de march se ralisent. Grce des simulations numriques, nous
avons montr que le design du module du temps rel nest pas neutre vis--vis
de la squence des marchs dnergie et de la dynamique concurrentielle. Les
designs qui scartent du type march et utilisent des systmes de prix p-
nalisant les transactions du temps rel devraient tre vits, dans la mesure du
possible, car ils provoquent des distorsions, des inemcacits et peuvent crer des
barrires lentre. Paralllement, la fermeture des marchs forward, dtermine
par la position temporelle de la gate closure, devrait tre le plus proche possible
du temps rel ahn de diminuer lampleur de ces distorsions. Ce cadre forma-
lis a permis ensuite dvaluer les consquences de lintgration de modules du
temps rel entre deux zones de contrle. Les rsultats des simulations num-
riques ont montr que lintgration des modules du temps rel est fondamentale
pour la cration de marchs rgionaux. De surcrot, lintgration de ces modules
doit tre suivie dune harmonisation adquate des designs, ahn que les enets de
lintgration proviennent des caractristiques conomiques des systmes lec-
triques intgrs et non des rgles du design.
3.3 Linvestissement dans les rseaux
Aprs avoir tudi les designs de march, il nous a paru intressant de nous
intresser aux problmes soulevs par cette ouverture concurrentielle. Ce travail
a t ralis lors de la thse de Vincent Rious co-encadre avec Jean-Michel Gla-
chant du groupe rseau Jean-Monnet. Pour assurer un accs libre et non discri-
minatoire aux infrastructures essentielles (rseau lectrique) la rforme a gn-
ralement ncessit une sparation verticale entre la production en concurrence
et le rseau lectrique en monopole naturel. La gestion de ce rseau a t dl-
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
3.4 Gestion de lalatoire et des incertitudes dans les rseaux dnergie 29
gue au gestionnaire de rseau de transport (GRT). La sparation des activits
de production et de transport a fait apparatre des problmes de coordination,
en particulier entre les investissements en production et en rseau. En enet les
temps caractristiques de construction dune centrale lectriques sont beaucoup
plus courts que ceux de cration de lignes de raccordement sur le rseau de
transport. Cette problmatique a fait lobjet de la thse de Vincent Rious co-
encadre avec Jean-Michel Glachant du Groupe Rseau Jean Monnet.
Ce travail cest donc articul autour dune analyse des mcanismes lis la
sparation des activits et des problmatiques lies. Un des objectifs de cette
thse a t de montrer que linvestissement du rseau est le cur de la coordi-
nation entre la production et le transport. La place centrale de linvestissement
du rseau dans cette coordination sexplique alors par les dinrences de rythme
temporel dinvestissement entre le transport et la production. tant donne cette
dinrence de dynamique temporelle, il faut faire connatre aux producteurs
lavance les prvisions de dveloppement du systme lectrique et les ventuels
signaux de localisation. Enhn il nous est apparu ncessaire de montrer comment
sorganise la coordination long terme entre la production et le transport dans
un systme lectrique libralis et du rle proactif que doit y tenir le GRT pour
coordonner de faon satisfaisante ces deux activits.
Pour raliser ces objectifs, il a t ncessaire de construire un simulateur
ainsi que des procdures doptimisation permettant de dhnir les meilleures so-
lutions pour une stratgie coordonne dinvestissement. La construction de ces
outils a demand un large panel de connaissances tant aussi bien sur le plan lec-
trique des rseaux (calcul de la rpartition de puissance par un load-fow), cono-
mique (signaux de localisation, investissement, incitations) que mathmatiques
(rsolution de systmes dquations non linaires, optimisation sous contraintes
mono et multi-objectif).
3.4 Gestion de lalatoire et des incertitudes dans les
rseaux dnergie
3.4.1 Introduction
Ce travail sur la coordination entre linvestissement en rseau et en produc-
tion nous a amen nous pencher sur lintgration de plus en plus importante et
encourage des moyens de production alatoires que sont les fermes oliennes et
photovoltaques, ces tudes ont fait lobjet de trois thses. La premire, celle de
Florent Maupas, co-encadre avec Jean-Michel Glachant du groupe rseau Jean-
Monnet, portant sur limpact de lolien sur les marchs de llectricit. Cette
thse mixte en conomie et gestion des rseaux eu pour objet de rechercher
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
30
lectriques
les impacts de lintgration de lolien, notamment sur lquilibre production-
consommation travers le mcanisme dajustement en temps rel.
La seconde, ralise par Herman Bayem et co-encadre avec Marc Petit et
Jean-Claude Vannier du dpartement nergie de Suplec sest intresse la
gestion locale des fermes oliennes et photovoltaques travers les conditions
dintgrations de celles-ci sur les rseaux de distribution. Nous avons recher-
ch de nouveaux moyens de dtermination des conditions de raccordement des
moyens de production alatoires. En enet actuellement, les tudes de raccorde-
ment se basent sur des cas dterministes o ne sont regards que les cas critiques
pour le rseau. En introduisant une dtermination probabiliste des moyens de
production et de la consommation, nous avons tudi comment modiher ces
tudes de raccordement.
Cette approche probabiliste et la gestion de lolien nous a amen nous
questionner sur les techniques doptimisation en prsence dincertitudes. Ainsi
une troisime thse sur ce sujet est en cours de ralisation. Henri Borsenberger
sest intress aux techniques doptimisations robustes en prsence dincerti-
tudes. Sa thse co-encadre avec Guillaume Sandou du dpartement Automa-
tique de Suplec, porte sur les modihcations apporter aux algorithmes dopti-
misation dterministes pour les rendre performants dans les cas stochastiques.
Cest thse bien que traitant dun sujet mathmatique est applique la probl-
matique de la gestion dun rseau dnergie. Ce dernier pouvant tre lectrique
ou de chaleur. Nous avons introduit les problmatiques dincertitudes sur les
cots de production, la quantit de production et enhn la demande.
3.4.2 Impact de la flire olienne sur les marchs de llectricit
Dans un premier temps, nous avons regard limpact de lintgration o-
lienne grande chelle sur lquilibre de production-consommation et les mar-
chs de llectricit au cours de la thse de Florent Maupas, co-encadre avec
Jean-Michel Glachant du groupe rseau Jean-Monet. Depuis les annes 1990, la
production dlectricit partir doliennes a t encourage dans un certain
nombre de pays par la mise en place de mcanismes de soutien. Ainsi, la hlire
olienne sest massivement dveloppe et reprsente dans certains pays une part
signihcative de la capacit totale de production lectrique. Cette hlire de pro-
duction transforme lnergie mcanique du vent en nergie lectrique qui est
ensuite injecte dans le rseau lectrique. Les oliennes exploitent ainsi une res-
source locale abondante sans mettre de gaz enet de serre. La hlire olienne
apporte donc une rponse des enjeux socitaux de premier plan tels que la
rduction des missions de gaz enet de serre ou la scurisation de lapprovi-
sionnement en lectricit. Cependant, la production lectrique des oliennes h-
rite de certaines caractristiques de sa ressource primaire : la puissance produite
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
par les oliennes est variable, peu prvisible et non contrlable. Ces caractris-
tiques distinguent la hlire olienne des hlires de production dites classiques
pour lesquelles la production est trs largement contrlable et prvisible.
cause de ces caractristiques, lintgration de la hlire olienne dans lquilibre
production-consommation soulve des enjeux conomiques et techniques im-
portants.
Ces enjeux de lintgration de la production olienne sont lis une spcih-
cit de la gestion du rseau lectrique. chaque instant, la production injecte
dans le rseau doit tre gale la consommation soutire du rseau. La pro-
duction olienne, par son caractre imprvisible, constitue ainsi un ala pour
les systmes lectriques. La gestion de cet ala revt non seulement un en-
jeu technique - car la sret du rseau en dpend - mais galement un enjeu
conomique - car la gestion dun ala a un cot. Cependant, le dh technique
pos par la production olienne nest pas dune nature entirement nouvelle.
La problmatique de gestion de lquilibre du systme nest pas apparue lors
de lmergence de la hlire olienne. Dune part, sur un rseau, la consomma-
tion dlectricit change sans pravis et constitue donc un ala pour lquilibre
production-consommation. Dautre part, les centrales des producteurs peuvent
subir des indisponibilits fortuites qui remettent galement en cause cet qui-
libre. Depuis la libralisation du secteur lectrique, la gestion de lquilibre du
rseau est organise autour de deux tapes : les marchs de llectricit et le
module du temps rel.
Dans le cadre de cette intgration de plus en plus importante de lolien
dans le mix nergtique de production, il est ncessaire de connatre limpact
de cet ala olien sur les marchs de llectricit. Ceux-ci permettent datteindre
un quilibre prvisionnel entre lonre et la demande en lectricit. Cependant
cet quilibre prvisionnel entre lonre et la demande ne garantit pas, en temps
rel, lquilibre physique entre la production et la consommation. Cest le mo-
dule du temps rel qui garantit en temps rel lquilibre physique des injections
et des soutirages. Le dh technique impos par le caractre imprvisible de la
production olienne est donc relev par deux modules de ralisation de lqui-
libre : dune part les marchs de llectricit, en particulier ceux organiss du
jour pour le lendemain et en infrajournalier et, dautre part, le module du temps
rel durant lequel un oprateur du systme coordonne les injections et les sou-
tirages. Lenjeu conomique de la gestion des systmes lectriques consiste
minimiser les cots de fourniture de la consommation, tout en respectant les
contraintes techniques et de sret du systme. Dans une situation hctive din-
formation parfaite sur les conditions du temps rel, nous pouvons considrer
que la programmation optimale des centrales lectriques est idale. Dans la si-
tuation relle dexploitation du systme lectrique, les centrales sont program-
mes dans une situation dinformation imparfaite sur les conditions du temps
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
32
lectriques
rel. La prsence de lala olien renforce encore le caractre incertain du temps
rel. Ds lors, le cot optimal de fourniture de la consommation sloigne de
son idal. Lenjeu conomique de la gestion de lala olien se situe prcisment
au niveau de la dinrence entre les cots de fourniture optimaux rel et idal.
Lenjeu conomique de la gestion de lala olien est donc li lemcacit de la
programmation des centrales. Cette programmation rsulte de deux types de d-
cisions : celles prises sur les marchs de llectricit dont la dimension physique
est importante et celles prises sur le module du temps rel. Par consquent, les
tapes dterminantes pour lenjeu conomique de la gestion de lala olien sont
le march du jour pour le lendemain, les marchs infrajournaliers et le module
du temps rel. Cest pourquoi lobjectif de ce travail est de dhnir une mthodo-
logie danalyse conomique de lenjeu conomique de la gestion de lala olien
et de lappliquer au cas particulier du systme lectrique franais. Lenjeu co-
nomique de la gestion de lala olien est fortement li la programmation des
centrales lectriques. Il faudra donc analyser comment la programmation des
centrales est infuence par ce nouvel ala. Une approche par simulation semble
alors simposer.
Nous avons dvelopp notre propre simulateur de la gestion de lala o-
lien. Celui-ci se dmarque des simulateurs existants sur plusieurs plans car il
modlise lensemble des contraintes techniques des centrales lectriques, la dua-
lit des modules en infrajournalier, les aspects dynamiques de la programma-
tion ainsi que la situation dincertitude au moment de programmer les centrales.
Les dimcults daccs aux donnes sont partiellement rsolues en exploitant les
onres dposes en infrajournalier par les producteurs lectriques franais qui
participent au module du temps rel. Ces onres donnent des informations im-
portantes sur les conditions hnancires et techniques dutilisation infrajourna-
lire des centrales lectriques. Cet outil de simulation modlise les deux mo-
dules de gestion de lquilibre en infrajournalier que sont les marchs infra-
journaliers et le module du temps rel en prenant en compte lorganisation de
ces deux modules, de la situation dinformation imparfaite pour la programma-
tion des centrales lectriques en infrajournalier ainsi que des contraintes tech-
niques dexploitation du systme. Lutilisation du simulateur a mis en vidence
une augmentation du cot de fourniture de la consommation pour grer en in-
frajournalier un ala sur lquilibre production-consommation. Ce surcot est
dhni comme le cot de gestion de lala. Lanalyse de ce cot a permis den
identiher le principal dterminant qui est la fexibilit limite des centrales lec-
triques. Cest cause du manque de fexibilit des centrales que lincertitude et la
dualit des modules infrajournaliers provoquent des pertes demcacit. Lappli-
cation du simulateur au cas du systme lectrique franais actuel a montr len-
jeu conomique de lintgration dune capacit olienne de 3 12GW. Le cot
de gestion de lala olien pour son responsable dquilibre varie largement en
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
fonction des hypothses retenues pour la liquidit des marchs infrajournaliers.
Les travaux raliss ont montr que les choix mthodologiques qui caractrisent
le simulateur dvelopp taient justihs. Le choix dune reprsentation pousse
de la fexibilit limite des centrales a des consquences importantes sur la com-
plexit du simulateur et requiert lutilisation de donnes prives des producteurs.
Les travaux ont cependant montr que cette reprsentation est trs importante
puisque le cot de gestion de lala olien est trs fortement infuenc par la fexi-
bilit limite des centrales. De plus, les rsultats montrent quil est ncessaire de
rendre compte de la fexibilit limite des centrales pour simuler les enets de la
dualit des modules en infrajournalier et les enets des incertitudes en infrajour-
nalier. Ce choix mthodologique est cependant moins fort si le simulateur est
appliqu un systme o les centrales trs fexibles abondent. Le choix mtho-
dologique retenu permet de matriser les temps de calcul. Nous avons montr
que sous certaines conditions, la reprsentation des centrales lectriques par des
agrgats de centrales permettait des gains de temps de calculs importants et
naltrait pas signihcativement la qualit des estimations numriques du cot de
gestion dun ala.
3.4.3 Intgration des moyens de productions incertains
Un autre aspect de laugmentation de la pntration des nergies renouve-
lable est la problmatique des conditions dinsertion des fermes oliennes. Cette
tude a t ralise lors de la thse dHerman Bayem. Le contexte politique,
conomique et nergtique europen est actuellement favorable une insertion
importante dnergie renouvelable (EnR) sur les rseaux lectriques. Laccrois-
sement de la production dnergie renouvelable et la drglementation des mar-
chs de llectricit font natre, dans le domaine de la gestion et de lexploitation
des rseaux, des problmes scientihques et techniques nouveaux. Ces problmes
sont induits par linsertion de nouvelles sources dnergie dans les systmes lec-
triques, non conus priori pour les accueillir. Lune des consquences de cette
arrive massive de nouvelles nergies est la modihcation de la structure des r-
seaux qui passera dune structure hirarchique avec des moyens de production
conventionnels de grande taille et centraliss une structure horizontale avec
de la production dcentralise (notamment renouvelable) dans les rseaux de
distribution.
Le nouveau systme lectrique est constitu de moyens de production
conventionnels (thermique, nuclaire, hydraulique), de moyens de production
non conventionnels (micro turbines, cognration, pile combustible, gother-
mie, biomasse, petite hydraulique, olien et photovoltaque), des rseaux de
transport et de distribution. En 2004 les nergies fossiles (ptrole, gaz) repr-
sentaient 65,8% de la production mondiale dlectricit et le nuclaire 15,8%. Les
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
34
lectriques
moyens de production conventionnels (centrales thermiques et nuclaires) sont
adapts au fonctionnement des systmes lectrique puisque leurs production est
contrlable mais ils utilisent les nergies qui prsentent plusieurs inconvnients
lis aux considrations environnementales : rserves limites, missions de gaz
enet de serre, traitement des dchets (notamment nuclaires). Ces considra-
tions environnementales ont en partie dict lvolution des systmes lectriques
vers une intgration massives des moyens de production non conventionnels
et particulirement de lolien et du photovoltaque qui feront lobjet de ce tra-
vail. La part dnergie renouvelable hors hydraulique dans la production mon-
diale est passe de 1,1% en 2000 (Papaefthymiou, 2006) plus de 2% en 2004. La
caractristique fondamentale de la production renouvelable est sa dpendance
aux conditions climatiques (le vent pour les fermes oliennes et lensoleillement
pour les fermes photovoltaques). Lexploitant du systme lectrique a donc un
pouvoir de contrle limit sur la quantit dlectricit en sortie des units de pro-
duction dnergie renouvelable. Plusieurs tudes montrent quun comportement
fuctuant peut tre caractris par des variables alatoires traduisant les varia-
tions sur une priode donne. Il en est de mme de la puissance produite par
les fermes oliennes et photovoltaques. Plusieurs paramtres du systme lec-
trique ont galement un comportement stochastique (disponibilit des moyens
de production conventionnelle, des ouvrages de rseau, la demande. . .). Ainsi la
puissance produite par les moyens conventionnels peut tre modlise par une
variable alatoire qui tient compte de la disponibilit des groupes de production.
Il est donc lgitime de penser quune modlisation probabiliste du systme lec-
trique serait adquate pour caractriser son fonctionnement. Et, par consquent,
que lanalyse des systmes lectriques par les mthodes probabilistes serait par-
ticulirement approprie.
Linsertion dune installation olienne ou photovoltaque sur un rseau lec-
trique (de transport ou de distribution) peut entraner des contraintes lies
dinrents aspects tels que : courants en rgime permanent et congestion de r-
seaux, plan de tension, courants de court-circuit, plan de protection, comporte-
ment dynamique et contribution aux services de rglage de tension et frquence,
stabilit des fermes oliennes, photovoltaques et du rseau lors de dfauts, qua-
lit de la tension, etc. La ralisation dtudes dimpact des nergies renouvelables
sur les rseaux est donc ncessaire pour analyser ces contraintes, anticiper les
problmes lis au dveloppement futur de ces nergies et rechercher des solu-
tions appropries. Ces tudes reposent en particulier sur la modlisation des uni-
ts de production dlectricit dorigine renouvelable. Les premires tudes din-
sertion dolien ou de photovoltaque dans les systmes lectriques ont t rali-
ses laide de mthodes dterministes ; ceci principalement cause de labsence
de modlisation probabiliste approprie et de limportance initialement relative-
ment rduite de ces units de production dans le parc de production. Ces ana-
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
lyses dterministes sappuient sur lexamen dun nombre restreint de situations
considres priori comme problmatiques ( les pires cas ) pour lesquelles on
vrihe la tenue du systme lectrique. On fait lhypothse implicite que les autres
situations pouvant se produire sont moins contraignantes. Lapproche probabi-
liste est une autre faon daborder le problme. Dans son principe, elle revient
considrer tous les cas possibles avec leur probabilit doccurrence ahn desti-
mer le risque de ne pas respecter une contrainte du systme. Les consquences
du non-respect de la contrainte seront bien sr mettre en regard de leur gra-
vit ou svrit pour le systme. Le gestionnaire de rseau devra au pralable
tablir une politique de risque. Lapproche probabiliste devrait ainsi permettre
de balayer lensemble des conhgurations (ou cas) possibles en tenant compte
des alas lis la production renouvelable, la disponibilit des units conven-
tionnelles et des lignes, la demande et donc de cerner plus hnement les risques
encourus avec le niveau de svrit et la probabilit doccurrence des situations
contraignantes. Lobjectif est alors de rechercher des solutions nouvelles tech-
niquement et conomiquement viables tout en garantissant la scurit des per-
sonnes et des biens. Bien sr, en plus dune modlisation probabiliste du systme
lectrique, lutilisation de ce type dapproche requiert le dveloppement dune
nouvelle mthodologie pour les tudes dinsertion et lutilisation doutils ap-
propris. Dautre part, limplmentation des solutions nouvelles auxquelles elles
pourraient conduire devra probablement sappuyer sur des moyens volus de
conduite et de contrle des fermes et des rseaux et pourrait ncessiter une vo-
lution du cadre rglementaire et contractuel.
Lanalyse des systmes lectriques dans ce contexte requiert la prise en
compte des dinrents alas induits par les nouveaux moyens de production et
par les autres lments du systme. Lobjectif principal de ce travail a t dva-
luer lapport de lapplication des mthodes probabilistes aux tudes dimpact de
lolien et du photovoltaque sur les systmes lectriques par rapport aux m-
thodes dterministes. La ralisation de ce travail est passe par trois tapes :
modlisation probabiliste du systme lectrique : il sagit ici de caractri-
ser la variation de plusieurs paramtres du systme par des variables ala-
toires, de dvelopper des mthodes de calcul des distributions de proba-
bilit de ces variables alatoires et danalyser les corrlations ventuelles
entre elles.
dveloppement de mthodes probabilistes pour les tudes de raccorde-
ment et dintgration des nergies renouvelables dans les systmes lec-
triques : les tudes de raccordement concernent lanalyse de limpact
dune ferme olienne ou photovoltaque raccorde au rseau de distribu-
tion ; les tudes dintgration concernent lanalyse de limpact dun parc
dolien et de photovoltaque sur un systme maill.
tudes de cas et comparaisons des mthodes probabilistes et dterministes.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
36
lectriques
Les nergies renouvelables (en particulier lolien et le photovoltaque) ont
des impacts sur le systme qui peuvent tre classs en deux catgories :
les impacts locaux qui concernent la qualit de tension, le plan de protec-
tion et la capacit daccueil du rseau,
les impacts globaux qui concernent la gestion de la production tous les
horizons de temps et le comportement dynamique du systme.
Les tudes dimpact des EnR ont pour principal objectif de vriher que toutes
les contraintes lies au fonctionnement normal des systmes sont respectes. Du
fait de leurs particularits, certaines des contraintes ont t amnages pour les
units de production renouvelable. Elles sont par exemple exemptes de rglage
primaire et secondaire de frquence dans les rseaux franais. Si le comporte-
ment stochastique de la puissance olienne ou photovoltaque est naturellement
mis en vidence, elle nest pas le seul paramtre du rseau avoir un comporte-
ment alatoire. Ainsi les paramtres suivants ont t modliss comme des va-
riables alatoires et leur distribution de probabilit a t calcule : la puissance
en sortie dune centrale constitue de plusieurs units de production conven-
tionnelle. Le principal facteur induisant lala est la disponibilit des units, la
disponibilit des lignes, la consommation dont la variation sur une anne est
modlise par une gaussienne, les productions olienne et photovoltaque qui
sont modlises par des distributions de probabilit non paramtres issues de
la transformation des distributions de probabilit des sources primaires. Ces
dinrentes distributions de probabilit et la structure de corrlation qui les lie
forment le modle probabiliste du systme lectrique. La gestion et lexploita-
tion dun systme lectrique passe par la ralisation danalyses dadquation, de
scurit statique et dynamique. Toutes ces analyses sont bases sur le calcul de la
rpartition de puissance. Ce calcul peut tre rsolu pour un point de fonctionne-
ment par la mthode dite dterministe. Vu les dinrents alas qui caractrisent le
fonctionnement rel dun systme, il est trs vite devenu ncessaire dintroduire
les notions de probabilit dans les calculs de rpartition et par consquent dans
les analyses des systmes lectriques. Il a t propos une mthode probabiliste
pour tudier limpact des EnR dans un systme lectrique. La principale carac-
tristique de cette mthode est quelle tient compte des alas lis aux variations
des dinrents paramtres du rseau sur la priode dtude. Ces alas sont pris
en compte par le biais des distributions de probabilit pralablement calcules.
Cette mthode est applicable aux trois types danalyse (adquation, scurit sta-
tique, scurit dynamique). Lapplication de la mthode probabiliste aux tudes
dimpact des EnR sur la scurit statique montre quelle permet de mieux simu-
ler le fonctionnement du systme sur une longue priode. Ltude probabiliste
permet de relativiser limpact des EnR sur un systme lectrique donn. En enet
limpact est dhni en termes de risque de dfaillance sur une longue priode et
non en termes de dfaillance sur des points de fonctionnement dits pire cas
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
comme cest le cas pour les tudes dterministe. De plus ltude dun grand
nombre de points de fonctionnement permet denectuer une analyse plus d-
taille de limpact des EnR. Les points de fonctionnement dfaillants ont ainsi
pu tre analyss et caractriss, les distributions de probabilit des tensions aux
nuds, des transits, de la rserve du systme ont t dtermines et analyses.
La mthode probabiliste a t applique ltude de raccordement dune ferme
olienne un rseau de distribution. Lanalyse du transit dans la ligne en amont
du point de raccordement montre que celui-ci est invers et que les congestions
sont rares. Compare celle obtenue par la mthode dterministe, la capacit de
la ferme raccordable obtenue par la mthode probabiliste est plus grande moyen-
nant un risque dcrtement de la production. Pour un risque dcrtement rai-
sonnable, cette possible augmentation de la capacit raccordable de la ferme
compare la capacit raccordable obtenue en dterministe, permet daccroitre
le productible annuel. Pour le cas dtude dun systme insulaire, lintgration
de 70 MW dEnR (soit un taux de pntration de 22,2%) ne modihe presque pas
la rserve et augmente globalement les tensions aux nuds de raccordement.
Les transits sont inverss dans les lignes voisines des nuds de raccordement.
De mme le taux de pntration maximal en puissance obtenu par la mthode
probabiliste est plus lev que celui obtenu par la mthode dterministe.
3.4.4 Modlisation et optimisation robuste de systmes com-
plexes : application aux rseaux nergtiques
Lintgration olienne et photovoltaque apporte des problmatiques de plus
en plus complexes, tant du point de vue physique et industriel, que du point
de vue mathmatique o, comme nous lavons vu, les probabilits et les incer-
titudes prennent une part importante dans les techniques doptimisation et de
rsolution. Cest pourquoi il nous est apparu ncessaire de continuer dans cette
voie de recherche en dveloppement des techniques doptimisation robuste des
systmes incertains. Cette tude a fait lobjet de la thse dHenri Borsenberger
co-encadre avec Guillaume Sandou du dpartement Automatique de Suplec.
La gestion optimale court terme (de quelques heures quelques jours)
des rseaux dnergie apparat comme un enjeu industriel et socital majeur.
En enet, dun point de vue purement technique, lmergence de nouvelles tech-
nologies (installations aux rendements levs, moyens de prdiction accrus par
exemple), combine la croissance de la puissance de calcul permet damlio-
rer la gestion des installations de production. Dun point de vue conomique,
louverture des marchs de lnergie induit un contexte fortement concurren-
tiel qui oblige les producteurs optimiser la gestion de leur parc de production.
Enhn, dun point de vue environnemental, il sagit de rduire les missions de
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
38
lectriques
polluants. Une telle rduction peut tre facilite par une gestion optimale et rai-
sonne de la production.
La gestion dun rseau de production, distribution ou stockage dnergie
est devenue un enjeu tant conomique que technique, devant respecter des
contraintes environnementales et lgislatives. Il apparat donc la ncessit dop-
timiser cette gestion ahn de maximiser le proht de fonctionnement du rseau
tout en respectant les contraintes imposes. Cette tche sest complexihe du
fait de la diversihcation des sources dnergie, qui peuvent coexister dans un
mme rseau (lectricit, gaz, houl, vapeur. . .). Par ailleurs, la modlisation de
ce type de systme doit intgrer les incertitudes lies la mconnaissance ou
la simplihcation du modle (en vue de son optimisation) ou des incertitudes
provenant de la nature prvisionnelle de la planihcation du fonctionnement du
systme (demande relle des consommateurs, fuctuations conomiques. . .). La
robustesse de cette gestion optimale est ds lors imprative, de sorte que les per-
formances du systme puissent tre garanties en dpit des dinrentes sources
dincertitude.
Il existe lheure actuelle relativement peu de travaux concernant la mod-
lisation et optimisation de rseaux multi-nergtiques, du fait de lmergence de
ce type de problmatique. Ainsi, ce sujet suscite un grand intrt de la part des
industriels, confronts la complexit et interconnexion, de plus en plus accrues,
de leurs systmes. Des travaux sur ce sujet ont t mens par le dpartement Au-
tomatique lors de la thse de Guillaume Sandou portant sur la modlisation et
optimisation dun rseau dnergie thermique. Les rsultats scientihques obte-
nus lissue de cette thse sont encourageants et ont permis de discerner les
verrous existants et donc de mieux cibler les travaux de recherche. En enet, il en
est ressorti la ncessit de prendre en compte les incertitudes prsentes dans le
systme, notamment dans la dtermination de la commande optimale du rseau.
De ce fait, un enjeu scientihque consiste en la dtermination dune commande
optimale robuste vis--vis des incertitudes.
Dans ce cadre, cette tude vis la mise au point de stratgies de gestion
optimale dun rseau dnergie. La dtermination de cette commande optimale
doit, dans la mesure du possible, utiliser des procdures systmatiques et gn-
riques, facilitant ainsi la transposition de la mthodologie dveloppe divers
types de rseau. Cette commande devra galement tre robuste vis--vis des in-
certitudes. En outre, un enjeu majeur est de mettre au point des algorithmes
numriques de rsolution qui soient rapides, hables et robustes.
Poursuivant les travaux initis dans le cadre du projet nergie de Suplec sur
loptimisation des rseaux dans un contexte dterministe, le but de cette thse a
t de rechercher des mthodologies assurant un pilotage robuste du rseau et
ce malgr plusieurs types dincertitudes.
Lapproche dterministe des problmes doptimisation sur des systmes in-
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
certains donne des rsultats imparfaits. Les optimums obtenus sounrent de deux
dfauts :
des rsultats prsentant des risques. Les algorithmes doptimisation
poussent gnralement le systme en limite de fonctionnement. Les pa-
ramtres du systme faisant varier le point de fonctionnement du systme
et/ou ses contraintes, il est possible que la solution initialement faisable
dans le cas dterministe se retrouve tre infaisable. On corrige gnrale-
ment ce dfaut en imposant une marge arbitraire au systme qui se traduit
gnralement par une moins bonne performance de la solution. La solu-
tion ainsi obtenue devient donc sous-optimale.
des rsultats sous-optimaux. Lorsque les paramtres varient, de manire
gnrale le lieu de loptimum varie aussi. Les paramtres incertains nont
pratiquement aucune chance dtre gaux leurs estimations. Par cons-
quent le lieu de loptimum obtenu a priori par la rsolution du problme
dterministe nest pas celui de loptimum rel que lon ne peut obtenir
qua posteriori aprs observation des paramtres.
Loptimisation de systmes incertains est communment appel optimisa-
tion stochastique ou optimisation robuste. Il est tonnant de constater quel
point les deux termes optimisation et robuste peuvent tre antagonistes. En ef-
fet, loptimisation du fonctionnement dun systme conduit gnralement se
placer en limite de contraintes, alors que la robustesse cherche au contraire
sen loigner en prenant des marges de scurit. Les premiers dveloppements
en matire doptimisation stochastique datent des annes 60. Les algorithmes et
la thorie ont connu un essor rcent ces 20 dernires annes particulirement
dans les problmes en chimie et en conomie ([AP98], [BN02], [BN99], [Sah04]).
Les mthodes dveloppes dans ce domaine font souvent appel des proprits
inhrentes leurs applications ce qui les rend parfois dlicates appliquer aux
problmes dnergie.
La premire dimcult rencontre dans des problmes doptimisation stochas-
tique par rapport cas dterministes est la complexihcation de problmes dj
complexes. En enet, dans la plupart des problmes poss chaque incertitude
ajoute une dimension de complexit au systme, ce qui a des rpercussions in-
vitables sur le temps de calcul.
La deuxime dimcult est un problme thorique. En gnral, un optimum
stochastique est loptimum qui minimise lesprance du critre, ou son cart
type. Dans les problmes de planihcation pour les rseaux dnergie en parti-
culier, on sintresse gnralement lesprance des cots de production. Ce-
pendant le calcul de cette esprance nest pas ais. Les systmes ne sont pas
linaires et le temps de calcul pour valuer la fonction critre est une barrire
pour le calcul de lesprance sur le domaine dincertitude. Par ailleurs, lobjectif
tant doptimiser une suite de dcisions chelonnes dans le temps, ces dci-
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
40
lectriques
sions peuvent tre amenes tre modihes selon les rajustements des prvi-
sions. Ces rajustements utilisent les nouvelles donnes acquises au cours du
temps : nouvelle donnes mto, cart de la consommation par rapport aux
prvisions. . . Lintgration du rajustement dans des problmes doptimisation
robuste fait appel la notion de recours, qui est un thme central en optimisa-
tion stochastique. Les travaux dans ce domaine font appels des mthodes de
programmation dynamique stochastique ([PCW00][GR]) et des mthodes multi-
stage stochastic programming [DR97].
Nous avons tudi les incertitudes suivantes, motives par des problmes
industriels rels :
incertitude sur la demande des consommateurs. Ce type dincertitude est
prsent sur la quasi-totalit des rseaux dnergie, car la demande de
consommation ne peut tre que prdite. La contrainte dquilibre consom-
mation / production est une contrainte forte des rseaux lectriques pour
en assurer le bon fonctionnement et la stabilit.
incertitude sur la capacit de production maximale. Ce type dincertitude
est prsent lors de lutilisation dnergies renouvelables telles que lnergie
photovoltaque ou lolien dans les rseaux lectriques, ou encore dans les
rseaux de chaleur avec lutilisation dincinrateurs de dchets.
incertitude sur les cots de production. Ce type dincertitude est prsent
sur presque toutes les installations, car les caractristiques de production
sont souvent identihes partir de points de mesure. Un cas particulier
qui pourra avoir son importance dans lavenir est celui des usines de co-
gnrations. En enet, si le prix de revente devient index sur le prix du
march de llectricit, alors de telles installations rentreront exactement
dans ce cadre.
Les mthodes dveloppes pour prendre en compte ces incertitudes appa-
raissent comme des extensions de la mthode choisie pour la rsolution du pro-
blme dterministe, savoir lutilisation de la programmation dynamique et de
la relaxation lagrangienne. Ceci est un point fort de nos travaux car les schmas
classiques doptimisation utiliss dans le cas dterministe restent utilisables.
3.5 Conclusion
Dans ce chapitre, nous avons trait de la modlisation et de loptimisation
technico conomique des rseaux lectriques de transport ou de distribution.
Nous avons dans un premier temps tudi les enets de la drgulation et de la s-
paration des entreprises de rseau de de production. Puis, au vu de linsertion de
plus en plus importante des nergie renouvelables photovoltaques et oliennes,
nous avons cherch comprendre les incidences induites par cette augmenta-
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
3.5 Conclusion 41
tion. Nous nous sommes intresss aux tudes de raccordement des moyens de
productions alatoire, puis nous avons cherch dvelopper de nouvelles tech-
niques doptimisation qui prennent en compte ces incertitudes. En parallle ces
travaux, nous avons appliqu les algorithmes doptimisation dvelopps au cas
des machines lectriques et plus particulirement des systmes de motorisation.
Ce dveloppement fait lobjet du chapitre suivant.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
42
lectriques
FicUvr 3.1 Amazing Stories juillet 1926
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Chapitre 4
Optimisation et modlisation des
machines lectriques
4.1 Principes de loptimisation des systmes lectro-
techniques
Dans beaucoup dapplications industrielles la question de la masse dun sys-
tme de motorisation est cruciale. On peut citer par exemple la conception des
satellites ou des aronefs, mais cest une ralit dans lensemble des industries
aronautique et spatiale. Dans de tels systmes une diminution de la masse ou
du volume induit une baisse de la consommation de carburant cest dire une
diminution du cot dexploitation et une hausse des prix de vente des appareils.
Lindustrie automobile est confront la mme problmatique, avec en compl-
ment la question du volume ncessaire lincorporation dans le vhicule et du
cot.
La recherche dune minimisation des pertes ou du cot nergtique est
bien sur une autre approche de loptimisation dun systme de motorisation. Cet
optimum est en gnral en contradiction avec le prcdent, ce qui nous amne
nous intresser des optimisations multiobjectifs en introduisant des fronts de
Pareto.
Nous avons choisi de nous intresser loptimisation et la modlisation des
systmes au dpartement nergie. Un systme est un ensemble fonctionnel de
composants runis au sein dune mme application. Par exemple pour un sys-
tme de motorisation, lensemble de la chane de traction est reprsente sur la
hgure 4.1 : moteur lectrique, rducteur de vitesse, convertisseur de mouvement,
source dnergie, convertisseurs lectriques. . .
La prise en compte de lensemble des composants de la chane de traction
est primordiale comme le montre la hgure 4.2 o une optimisation systme est
43
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
44 Cu~vi1vr 4 : Optimisation et modlisation des machines lectriques
LOAD
bat
U
bat
i
M
i
M
u
r

r
r

v
x

em
C
r
C
DC/DC Converter
Lead Screw Device
Speed Reducer
DC Motor
FicUvr 4.1 Reprsentation schmatique dune chane de traction lectrique
ralise. Nous avons en fonction de contraintes spcihques une application,
enectu une optimisation de la masse dun ensemble rducteur de vitesse et
moteur. Sur la hgure sont tracs les masses totales et de chaque composant en
fonction dune contrainte variable sur le rapport de transformation maximal du
rducteur de vitesse. Cest dire que chaque point de ces courbes est une op-
timisation qui arrive en butte sur cette contrainte variable. En permettant un
rapport de rduction maximal plus important, loptimisation rduit le volume de
la machine, mais augmente celui du rducteur. Il sen suit une rduction globale
de la masse. La masse totale du systme (en vert) est donne par la somme des
masses du rducteur et du moteur. Si lon augmente (relche) la contrainte sur le
rapport de transformation, on diminue la masse totale du systme : la masse du
moteur (en rouge) diminue tandis que celle du rducteur (en bleu) augmente. Op-
timiser seulement sur un point de fonctionnement la masse du moteur nest pas
recommand dans ce cas, il faut prendre en compte lensemble des composants
du systme. noter dans cet exemple le point optimal minimum de la masse
totale du systme : partir dune contrainte sur le rapport de rduction au-del
de 9,5.10
4
, laugmentation de la masse du rducteur nest plus compense par la
diminution de la masse du moteur.
Composer des optimisations systmes implique de raliser des modles mul-
tiphysiques prenant en compte tous les aspects des composants : lectrique, m-
canique, thermique. . .Ces modles se doivent dtre relativement simples, ou tout
du moins dtre rapidement calculables pour tre implments dans les codes
doptimisation.
Dans la suite de cet expos, nous allons prsenter deux exemples doptimi-
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
4.2 Optimisation des cots de fabrication et du rendement dune gamme de
machine 45
5 5.5 6 6.5 7 7.5 8 8.5 9 9.5 10
4
x 10
0.04
0.06
0.08
0.1
0.12
0.14
0.16
0.18
0.2
r
max
M
m
,

M
r

M
m
M
r
[ M
m
+ M
r
] / 2
FicUvr 4.2 Variation des masses optimales du systme en fonction de la
contrainte sur le rapport de transformation
sation. Le premier ralis lors de la thse de Xavier Jannot, prsente une forme
originale de problme doptimisation : la rduction des cots de fabrication dune
gamme dalternateurs. Il sagit encore une fois doptimisation systme, bien que
celui-ci nest pas cette fois ci compos dlments physiques, mais plutt de
types dalternateurs dans une mme gamme.
Par la suite, nous prsenterons une premire ralisation doptimisation dun
convertisseur de puissance. Cette optimisation multi-objective servira de base
pour une conception future complte du systme de motorisation : batterie,
convertisseur dnergie, rducteur de vitesse, convertisseur de mouvement, mo-
teur.
4.2 Optimisation des cots de fabrication et du ren-
dement dune gamme de machine
Avec la rcente augmentation du prix des matires premires les construc-
teurs de machines lectriques cherchent amliorer leurs dimensionnement ahn
de rduire la quantit de matriaux actifs utiliss ainsi que les cots de fabrica-
tion. Cette diminution du poids et de la quantit de matire utilise infuence
directement le rendement nergtique de la machine. Cette augmentation des
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
pertes est en contradiction avec les nouvelles rglementations europennes qui
recommandent au contraire une augmentation des rendements. Les construc-
teurs sont donc confronts un choix qui peut tre reprsent par une optimi-
sation sous contraintes : rduire le cot en matire des machines en respectant
des contraintes de rendement.
Une autre problmatique introduite lors de cette tude vient du fait que les
constructeurs tendent harmoniser et mutualiser les conceptions au sein dune
mme gamme. Pour un mme type de machine (ici un alternateur), une partie
de la conception est gard le mme pour toute une gamme de puissance. Dans
lapplication propose par la suite, certaines caractristiques comme le rayon de
la machine, le nombre de paires de ples sont gardes constants pour une mme
gamme de puissance, la dinrence se faisant sur la longueur de tle utile.
Plusieurs stratgies doptimisation ont t proposes, la premire qui pourra
servir de rfrence est doptimiser chaque machine de la gamme sparment. La
seconde dite classique est celle utilise habituellement dans le monde industriel.
Elle consiste optimiser la machine de plus forte puissance puis la dcliner en
passant des modles de plus basse puissance. Cette dclinaison se faisant en
optimisant les paramtres pouvant varier pour chaque machine de la gamme.
Une nouvelle approche que nous avons propose et qui apporte des gains
supplmentaires en terme de rduction des cots de fabrication est doptimiser
globalement la gamme en tenant compte de la rpartition des productions sui-
vant les puissances.
Pour raliser ces optimisations, il est ncessaire davoir un modle phy-
sique des composants, cette modlisation sera dveloppe dans le prochain para-
graphe. Nous prsenterons par la suite les rsultats des optimisations de gamme
en mono-objectif puis en multi-objectif en prenant en compte les rendements.
4.2.1 Modlisation de la machine
Le type de machine tudie ici est une gnratrice synchrone ple saillant.
Ahn de raliser loptimisation, nous avons besoin dun modle global et phno-
mnologique de la machine. La gomtrie type est prsente sur la hgure 4.3.
Les dents du stator sont de section constante (excepte dans les pieds et les fond
dencoches). Le bobinage du rotor est ordonn et le nombre de paires de ples
est hx 2.
Un modle lectromagntique tenant compte de la saturation a t dve-
lopp, et permet de dterminer partir de donnes gomtriques, les densits
de fux magntiques, la force lectromotrice et les courants dans les bobinages.
Les rsultats de ce modle servent ensuite de point dentre une modlisation
des pertes. Les pertes joules, les pertes fers ainsi que les pertes mcaniques sont
values partir dun premier point de temprature. Les valeurs de ces pertes
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
machine 47
FicUvr 4.3 Gomtrie de la machine tudie
dpendant grandement de la temprature, celle-ci est ensuite recalcule puis par
itrations successives (mthode du point hxe), les pertes et la temprature sont
rajustes suivant lorganigramme de la hgure 4.4.
Un modle complet de la machine et des interactions entre composants est
donn dans le chapitre douvrage [1].
4.2.2 Optimisation dune gamme dalternateurs
Lapproche de conception optimale nous permet de trouver le meilleur sys-
tme tout en respectant les spcihcations aussi appeles contraintes. Loptimisa-
tion dune gamme de machines lectriques bases sur le modle prsent dans le
paragraphe prcdent est complexe, il y a en enet des variables mixtes (entires
et relles), la fonction objectif est non linaire et surtout non drivable analyti-
quement. Cest pourquoi le choix de lalgorithme doptimisation sest port sur
une mthode dalgorithme gntique et dvolution dinrentielle. Ce choix est
renforc par la ncessit par la suite de raliser des optimisations multi-objectifs
traites laide de lalgorithme NGSA2 [KAST02].
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
FicUvr 4.4 Organigramme de calcul du couplage lectriquethermique
Optimisation mono-objectif
Plusieurs optimisations ont t ralises ahn de minimiser dinrentes fonc-
tions objectifs. Une gamme de trois machines a t considre couvrant des puis-
sances variant de 125kVA 180kVA. La fonction objectif dans chacune des opti-
misations ralises dans ce paragraphe est reprsentative du cot des matires
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
machine 49
utilises lors de la fabrication :
WRSG
cost
=
weight.material
cost
(4.1)
Les contraintes peuvent tre spares en deux types. Les premires portent
directement sur les variables dentres du problme, aussi appel paramtres de
loptimisation. Ce sont les bornes de recherche et sont hgures dans le tableau
4.1.
WRSGs parameters Initial design Range
Outer stator diameter [mm] 390 [390,420]
Inner stator diameter [mm] 270.5 [255,290]
Slot diameter [mm] 317.5 [300,360]
Tooth width [mm] 10.2 [5,24]
Rotor pole width [mm] 86.5 [70,100]
Rotor pole opening factor 0.7045 [0.6,0.8]
Machine length [mm] 410 [210,420]
Conductor number 6 [5,12]
Relative cost 1
T~nir 4.1 Paramtre de la machine optimiser et leurs bornes
Le second type de contraintes porte sur des grandeurs calcules par la mo-
dlisation. Les tempratures du stator et du rotor ne doivent pas augmenter au-
dessus de 125K. Les rendements doivent tre suprieurs une certaine valeur
dpendant de la puissance de la machine considre. Ces contraintes sont indi-
ques dans le tableau 4.2.
Quantity Constraint
Stator temperature rise < 125K
Rotor temperature rise < 125K
125kVA WRSG emciency > 91.6%
T~nir 4.2 Contraintes de loptimisation sur les valeurs de sortie
Une premire approche est doptimiser le cot de chaque machine prise s-
parment. Cela conduit la ralisation de trois calculs dont le rsultat nous don-
nera le cot minimum de fabrication dune machine si celle-ci tait la seule tre
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
produite. Dun point de vue industriel cette approche nest pas ralisable car elle
conduit trois sections de rotor dinrentes. Cela oblige le constructeur dve-
lopper des chanes de production dinrentes pour chaque machine et ne permet
aucune mutualisation des processus des fabrication. Les rsultats des optimisa-
tions pour ces trois machines de puissance dinrente sont prsents dans le ta-
bleau 4.3. Daprs ces rsultats, les machines 165kVA et 180kVA optimises sont
en limite de contrainte pour les tempratures de rotor et de stator ainsi que sur
le diamtre extrieur du stator. Pour la machine 125kVA, la bute en contrainte
se fait cette fois-ci sur la temprature du rotor et sur le rendement. Ces rsultats
sont en accord avec lexprience de conception de machines : les machines les
plus volumineuses sont en gnral limites par la hausse de temprature, alors
que les machines plus petites sont connues pour avoir des rendements moins
performants.
variables 125 kVA 165 kVA 180 kVA
Outer stator diameter [mm] 393.7 420 420
Inner stator diameter [mm] 255.2 276.9 278.9
Slot diameter [mm] 311 326.1 331.5
Tooth width [mm] 11.1 12.9 12.4
Rotor pole width [mm] 89 100 96.9
Rotor pole opening factor 0.748 0.7415 0.7415
Machine length [mm] 285.2 331.7 349.9
Conductor number 8 6 6
Stator temperature rise [K] 83.0 125 125
Rotor temperature rise [K] 124.8 125 125
Emciency [%] 91.60 91.79 92.34
Relative cost 0.7129 0.8693 0.9589
T~nir 4.3 Paramtres des machines suite aux optimisations individuelles
Approche classique de loptimisation dune gamme
Les rsultats prcdents, bien que garantissant pour chaque puissance de ma-
chine la meilleure solution, ne peuvent tre utiliss dans les tapes de production
industrielle. Les solutions proposes imposeraient de modiher la chaine de mon-
tage et de production lors du changement de type de machine ce qui induirait des
cots de production bien plus importants que le gain fait sur chaque machine.
Pour la rduction de ces cots une contrainte forte est dimposer une section
identique pour les trois machines de puissance dinrente. Cest dire que les
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
machine 51
trois machines ne seront dinrentes que par le nombre de conducteurs dans les
encoches et leur longueur. Lapproche classique, prsente ici pour des raisons
de comparaison, se base sur loptimisation de la machine la plus contraignante
(la plus puissante en gnrale) puis le nombre de conducteurs et la longueur
des autres machines sont optimiss en imposant la mme section. En se basant
sur les rsultats prcdents pour la machine 180kVA et en optimisant le nombre
de conducteur et la longueur des machines 165kVA et 125kVA, on obtient les
rsultats du tableau 4.4.
variables 125 kVA 165 kVA 180 kVA
Outer stator diameter [mm] 420 420 420
Inner stator diameter [mm] 278.9 278.9 278.9
Slot diameter [mm] 331.5 331.5 331.5
Tooth width [mm] 12.4 12.4 12.4
Rotor pole width [mm] 96.9 96.9 96.9
Rotor pole opening factor 0.7415 0.7415 0.7415
Machine length [mm] 261.8 349 349.9
Stator temperature rise [K] 87.1 97.8 125
Rotor temperature rise [K] 90.3 106.9 125
Emciency [%] 91.61 92.53 92.34
Relative cost 0.7526 0.9568 0.9589
T~nir 4.4 Paramtres des machines aprs optimisation classique de la gamme :
optimisation de la machine de plus forte puissance, puis dclinaison pour les
machines de plus faible puissance
Cette optimisation nous donne des rsultats moins optimaux que celles
mono-objectif ralises pour chaque machine prise sparment, mais permet
une mutualisation des processus de fabrication entrainant une diminution des
cots. Les cots relatifs sont comparer la conception initiale des machines
dans le tableau 4.1.
Une autre possibilit serait de prendre comme optimisation initiale celle
dune machine de plus faible puissance, puis augmenter la longueur jusqu trou-
ver celle correspondant aux puissances 165kVA et 180kVA. Lalgorithme dopti-
misation nous montre que cela nest pas faisable, car dans ce cas aucune solution
optimale nest trouve. En enet laugmentation de la longueur de la machine est
plus grande que la limite haute hxe dans les contraintes. Cela conhrme bien
lapproche faite par les constructeurs qui permet en plus de garder la gamme
de machine dans un volume donn. Mais hxer la section partir de la machine
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
la plus puissante nest pas la meilleure optimisation lorsque lon considre la
gamme entire comme nous allons le voir dans la section suivante.
Optimisation simultane dune gamme de machine
Une nouvelle approche que nous avons propose est doptimiser lensemble
des paramtres des trois types de machines considres. Il sagit donc de ra-
liser une seule optimisation prenant en compte lensemble des paramtres. La
fonction objectif globale devient :
range
cost
=
3
i=1
i
.WRSG
cost,i
(4.2)
La pondration
i
correspond la rpartition des ventes lintrieur de la
gamme.
125 kVA 165 kVA 180 kVA
Cas I 8/10 1/10 1/10
Cas II 1/10 8/10 1/10
Cas III 1/10 1/10 8/10
T~nir 4.5 Pondrations des ventes pour les trois scnarios envisags
Le nombre de variables doptimisation est augment, puisquil comprend
maintenant trois longueurs de machine et trois nombres de conducteurs par
encoche (une pour chaque puissance). Les autres paramtres restent inchangs
(voir le tableau 4.1). Cela nous conduit raliser une optimisation sur douze
paramtres sous les contraintes du tableau 4.2. Nous avons tudi trois cas cor-
respondants aux dinrentes pondrations du tableau 4.5 sur les volumes de vente
des machines au sein la gamme. Les rsultats de ces trois optimisation sont re-
ports dans les tableaux 4.6,4.7 et 4.8.
Dans les trois cas prsents, le gain relatif en terme de cot de production
est report sur la dernire ligne du tableau. Pour les trois rpartitions de vente,
nous obtenons grce cette nouvelle approche, une baisse des cots de fabrica-
tion. Le gain tant le plus important lorsque les volumes de vente sont en faveur
des machines de petite puissance. La troisime conhguration apporte un gain
relativement faible (mais non ngligeable pour le constructeur), car elle favorise
la machine de forte puissance, ce qui est dj le cas pour la technique dopti-
misation classique. Nous pouvons aussi observer une rpartition des contraintes
atteignant les bornes hxes suivant les distributions, les puissances plus faibles
que celle favorise tant en bute sur le rendement minimum, alors que les puis-
sances plus grandes le sont sur la temprature.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
machine 53
125 kVA 165 kVA 180 kVA
Outer stator diameter [mm] 412
Inner stator diameter [mm] 267.6
Slot diameter [mm] 324.1
Tooth width [mm] 12.1
Rotor pole width [mm] 94.7
Rotor pole opening factor 0.7472
Machine length [mm] 266.2 360.7 394.8
Rotor temperature rise [K] 110.4 124.9 124.9
Emciency [%] 91.60 92.53 92.40
Relative cost 0.726 0.9372 1.0135
Relative average cost per machine 0.7759
Saving compared to classical approach - 2.24 %
T~nir 4.6 Paramtres optimaux des trois types de machine dans le cas du
scnario I
125 kVA 165 kVA 180 kVA
Machine length [mm] 266.1 306.6 355
Stator temperature rise [K] 78.6 125 111.7
Rotor temperature rise [K] 89.5 125 125
Emciency [%] 91.89 92.00 92.59
Relative cost 0.7753 0.8712 0.9857
scnario II
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
125 kVA 165 kVA 180 kVA
Machine length [mm] 261.3 348 352.1
Rotor temperature rise [K] 92.5 110.1 125
Emciency [%] 91.60 92.51 92.35
Relative cost 0.7478 0.95 0.9596
scnario III
4.2.3 Optimisation multi-objectif du cot et du rendement
Nous avons par la suite travaill sur une optimisation multi-objectif base sur
le mme modle et les mmes contraintes. Nous optimisons conjointement, le
rendement et le cot de fabrication. La premire partie de cette section se plaait
du cot constructeur de machine, ou seul le cot de fabrication tait important ;
ici on se place cot consommateur ou utilisateur de la machine. Celui-ci doit
faire un choix entre le cot dachat (reli au cot de fabrication) et le rendement
long terme de la machine. Un exemple de calcul doptimisation multi-objectif
est reprsent sur la hgure 4.5, o le front de Pareto pour la machine 125kVA
est trac en trait plein. Cette optimisation est ralise laide dun algorithme
volutionnaire coupl un mode de rsolution de type NGSA2 [KAST02] et est
base sur le scnario 1 de la section prcdente. Nous avons aussi report les
caractristiques rendementcot pour les autres machines de la gamme pour
chaque point du front de Pareto.
Le premier objectif est de minimiser le cot de fabrication de la gamme de
machine dans le cadre dune rpartition donne. Le second objectif est le ren-
dement de la machine 125kVA (et uniquement celui-ci). On peut tracer un front
de pareto de cette double optimisation : en trait plein sur la hgure 4.5. Dautre
part, pour chaque point de ce front de pareto, on obtient un jeu de paramtre
pour lensemble des machines de la gamme, on peut donc en dduire un cot et
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
4.3 Prototypage virtuel en lectronique de puissance 55
FicUvr 4.5 Front de Pareto, cot en fonction du rendement pour la machine de
125kVA. Les rendements des autres machines de la gamme sont reports sur la
hgure.
un rendement. On trouve alors les nuages de points : indiqu par des o pour les
machines 180kVA et des + pour les machines 165kVA. Ces points ne forment
pas une courbe, mais sont rpartis en plusieurs paquets. On observe une coupure
vers 1.0 en cot et .93 en rendement.
La hgure 4.5 peut tre une aide prcieuse pour aider au choix des carac-
tristiques dune gamme de machine. En enet, on saperoit quil vaut mieux
pnaliser un peu le rendement (au-dessous de 0.93) de la machine 125kVA ou
son cot (au-dessous de 1) ; les rendements des autres machine de la gamme sen
trouvant largement amliors. Dautres recherches sur ce sujet sont en cours et
ont donn lieu une publication actuellement en cours dexamen.
4.3 Prototypage virtuel en lectronique de puissance
4.3.1 Introduction
Aprs avoir tudi et appliqu ces techniques doptimisation au cas des
chanes de traction, il nous est apparu ncessaire de nous intresser un autre
composant des systmes lectrotechniques, les convertisseurs de puissance. En
enet avant leur intgration complte dans le processus doptimisation systme, il
est ncessaire davoir des modles simples et des applications testes. Cest pour-
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
quoi nous avons travaill avec Pierre Lefranc et Xavier Jannot sur loptimisation
et le prototypage virtuel dun convertisseur de puissance dc-dc de type hacheur
abaisseur. Cette tude a fait lobjet dun article accept dans le journal IET Power
Electronics. Ce type de convertisseur est largement utilis, que ce soit dans les
applications domestiques, automobiles ou aronautiques. Les contraintes lec-
triques, thermiques, mcaniques, dencombrement, et de cot imposes par les
constructeurs requirent une approche de dimensionnement par optimisation
multiphysique qui doit prendre en compte les divers impacts de lenvironnement
sur le convertisseur, mais aussi ceux induits par le convertisseur lui-mme. Di-
verses approches ont dj t proposes dans la littrature, ainsi des techniques
doptimisation ont t dveloppes pour le prototypage dun convertisseur ha-
cheur multi-phase [HAD05], mais aussi pour un abaisseur (buck) DC-DC pour
une application automobile [TD03] et pour des convertisseurs DC-DC faible ten-
sion [XTF97]. Dautres auteurs se sont focaliss sur des techniques de prototy-
page virtuel en plusieurs tapes : ceci permet aux concepteurs de commencer
donner des premiers lments de rponse sur la dhnition du convertisseur.
On parle alors de pre-dimensionnement. Dans un second temps, le concepteur
amne ses choix, les modlisations, les optimisations ahn de dhnir de plus en
plus prcisment le systme [KCPC09]. Le pr-dimensionnement est trs impor-
tant dans le cas dapplications industrielles avec une forte contrainte sur le cot
(automobiles par exemple). Le but de ce travail a t de proposer une technique
de prototypage virtuel laide dune optimisation bi-objectif dans le contexte
dune application industrielle. Pour cet exemple, les donnes des composants du
convertisseur (MOSFETS, diodes et changeurs de chaleur. . .) sont issus directe-
ment des donnes techniques des constructeurs. Nous nous sommes intresss
en guise dapplication et dexemple au prototypage dun hacheur abaisseur (24V-
12V) dune puissance dun kW sans hltre dentre. Les variations de la tension
dentre et de la puissance de sortie ont t prises en compte dans la formulation
des contraintes : la tension dentre v
i
varie entre 20V et 30V, et la puissance de
sortie entre 100W et 1kW.
La mthodologie de prototypage virtuel applique ici au convertisseur abais-
seur peut tre applique aux autres convertisseurs continu-continu en adaptant
les contraintes. Comme prcdemment, un modle simple (permettant loptimi-
sation) a du tre dvelopp, et il sera prsent dans la prochaine section o il
servira de point dentre des mthodes doptimisations dont les rsultats seront
exposs dans la section 4.3.3.
4.3.2 Modlisation du hacheur abaisseur
Les convertisseurs de type hacheur sont trs souvent utiliss dans des appli-
cations industrielles et grand public (alimentation de carte CPU par exemple).
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Nom de variable Dhnition
K MOSFET
D diode
L inductor
C capacitor
R load resistor
v
i
input voltage
v
o
output voltage
P output power, P = v
2
o
/R
d duty cycle
T switching period
F switching frequency
i
L
inductor current ripple
i
Lm
minimum inductor current value
i
LM
maximum inductor current value
v
o
output voltage ripple
T~nir 4.9 Notations utilises pour le convertisseur abaisseur
On trouvera un schma de principe du hacheur tudi sur la hgure 4.6. Le sys-
tme est constitu de composants classiques : un MOSFET K, une diode D, une
inductance L et un condensateur C. Les formes dondes (sans pertes, cest dire
en utilisant des composants parfaits) sont prsentes sur la hgure 4.7, o i
Lm
et
i
LM
sont dhnis respectivement par les valeurs minimale et maximale du cou-
rant dinductance. La hgure 4.8 montre les ondulations de la tension de sortie
v
o
. On utilisera par la suite les notations du tableau 4.9.
FicUvr 4.6 Schma lectrique de principe du hacheur
Plusieurs hypothses ont t prises en compte pour llaboration du mo-
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
FicUvr 4.7 Forme donde du courant dans linductance et le MOSFET
dle du convertisseur. Elles nous permettent, sans nous loigner de la ralit, de
simpliher les modles utiliss et de les incorporer ainsi dans une approche dop-
timisation multi-critre et multi-physique. On considrera donc les hypothses
suivantes :
la tension dentre v
i
est constante,
les ondulations de tension v
o
sont ngligeables par rapport la valeur
moyenne de v
0
,
la valeur de la rsistance de la bobine est ngligeable,
le courant i
L
volue linairement,
les temps de monte et de descente dans le MOSFET sont ngligeables face
au temps de conduction,
il ny a pas de saturation dans le matriau magntique de linductance L.
partir de ces hypothses, nous avons labor un modle fonctionnel du
convertisseur dont les quations de base sont :
d =
v
o
v
i
(4.3)
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
FicUvr 4.8 Ondulations de la tension et du courant de sortie
i
LM
=
P
dv
i
+
v
i
(1 d) d
2LF
(4.4)
i
lm
=
P
dv
i
v
i
(1 d) d
2LF
(4.5)
v
o
=
v
i
d (1 d)
8LCF
2
(4.6)
Ces quations servent de point de dpart la modlisation des pertes dans
les composants du convertisseur. Nous avons considr les pertes par conduction
et par commutation dans linterrupteur MOSFET et par conduction seulement
dans la diode. Pour le MOSFET, lvaluation des pertes par conduction se fait en
fonction des tensions dentre et de sortie, de la puissance de sortie, de la priode
de fonctionnement, de linductance, mais aussi de lvaluation de la temprature
de jonction entre le composant et le boitier T
jK
. Finalement une relation pour
les pertes par conduction dans le MOSFET est donn par une quation du type :
P
Kcond
= f
1
(v
i
, v
o
, P, T, L, T
jK
) . (4.7)
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
De la mme manire, les pertes par commutation sont fonction des tensions,
de la puissance, de linductance et de la priode de fonctionnement :
P
Ksw
= f
2
(v
i
, v
o
, P, L, T) . (4.8)
Les pertes totales dans le MOSFET sont donnes par la somme des pertes par
conduction et commutation :
P
K
= P
Kcond
+P
Ksw
= f
1
(v
i
, v
o
, P, T, L, T
jK
) +f
2
(v
i
, v
o
, P, L, T) . (4.9)
Les pertes dans la diode sont des pertes par conduction, et sont values
partir dune modlisation des donnes constructeurs en fonction de la tension
ses bornes, du courant, de sa temprature de jonction (T
jd1
), de la priode de
dcoupage et du rapport cyclique d.
P
d1cond
= P
d2cond
= f
3
(T, d, i
d
, v
d
, T
jd1
) . (4.10)
Linductance est le troisime composant tre modlis. Il sagit, partir du
modle fonctionnel et des pertes (dans le fer P
iron,vol
vol
core
et le cuivre P
copper
),
de caractriser son volume et llvation de temprature au sein de linductance.
On obtient ainsi des quations analytiques permettant de calculer les dinrents
paramtres de cette dernire.
Enhn le dissipateur de puissance est un lment passif qui doit tre dimen-
sionn. Il est reprsent par une rsistance thermique value en fonction de
sa gomtrie. Un modle analytique a t dvelopp pour tenir compte de la
longueur du dissipateur dans la valeur de la rsistance thermique globale. Pour
cela, les donnes constructeur sont utilises pour raliser une approximation du
type :
R
th
(L
hs
) = a + b.e
c.L
hs
(4.11)
Une modle lectrothermique a t mis au point en tenant compte des tem-
pratures de jonction et des rsistances thermiques entre les composants et le
botier puis entre le botier et le dissipateur. Ce modle est reprsent sur la
hgure 4.9.
4.3.3 Dfnition de loptimisation et prototypage virtuel du ha-
cheur
Pour raliser une optimisation, il faut dhnir les objectifs ainsi que les
contraintes. Comme dans de nombreuses applications lectrotechniques, la ten-
dance actuelle est la rduction des volumes et des pertes. Nous dhnirons deux
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
FicUvr 4.9 Modle lectro-thermique
objectifs, qui nous le verrons dans les rsultats sont contradictoires ou en op-
position. Il sagit pour le premier de proposer un dimensionnement le plus petit
possible (optimisation du volume du convertisseur) et pour le second doptimiser
le rendement (diminuer les pertes).
Le rendement est dhni classiquement par le rapport entre la puissance de
sortie et la somme de celle-ci et des pertes :
=
P
P + P
K
+ P
dcond
+ P
iron,vol
vol
core
+ P
copper
(4.12)
Il dpend fortement du point de fonctionnement comme constat sur la h-
gure 4.10 qui montre la variation de celui-ci en fonction de la puissance et de
la tension dentre. Les spcihcations techniques imposent en gnral un rende-
ment minimum qui est dhni par :
F
1
= min
v
i
,P
(4.13)
Nous reprenons ici, cette dhnition pour le premier objectif de notre optimi-
sation. Le volume total du convertisseur est notre deuxime objectif, il est dhni
par le volume du plus petit paralllpipde entourant le dissipateur (hgure 4.11).
En ngligeant le volume du circuit imprim, le volume total est dhni par :
F
2
= max (L
hs
, L
L
) . max (E
hs
.E
L
) . (H
hs
+ H
L
) (4.14)
Les huit variables doptimisation sont donnes dans le tableau 4.10, ainsi que
leur domaine de variation respectifs.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
P [W]
V
e

[
V
]
Efficiency Map for the range of working conditions

100 200 300 400 500 600 700 800 900 1000
20
21
22
23
24
25
26
27
28
29
30
0.95
0.955
0.96
0.965
0.97
0.975
FicUvr 4.10 Valeurs du rendement dans le plan puissancetension de sortie
FicUvr 4.11 Dhnition de la gomtrie du disspipateur
Ces variables dhnissent une optimisation valeurs mixtes entires et
relles. Les variables discrtes ne peuvent prendre quun nombre limit de va-
leurs dans le domaine autoris, ainsi le nombre de tours de la bobine de lin-
ductance est un entier entre 1 et 200, tandis que les valeurs de la capacit du
condensateur ne peuvent prendre que des valeurs normalises ({1, 1.5, 2.2, 3.3,
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Optimisation variable min value max value
F [kHz] : switching frequency 10 20
C [F] : capacitor value 10
6
10
1
e [mm] : air gap 0 3,3
A [cm2 ] : winding square area 0 100
S [m2 ] : square cross section of the feritte core 10
6
36 10
4
s [m2 ] : cross section of the elementary wire 10
7
16 10
6
n : wiring turn number of the inductor 1 200
L
hs
[mm] : heatsink length 50 200
T~nir 4.10 Dhnitions et domaine de variation des paramtres de loptimisa-
tion
4.7, 6.8} multiplies par des puissances de dix.
Imposer des contraintes, permet de sassurer que les solutions du problme
sont ralistes et ralisables. Dans cette application, nous avons tenu compte des
limites physiques, des conditions dutilisation et des contraintes gomtriques.
Ces contraintes sont exposes dans la liste suivante.
CCM Le hacheur tant destin tre utilis sous condition de conduction
continue (CCM), la valeur minimale du courant dans linductance doit tre posi-
tive :
i
Lm
> 0 (4.15)
Ondulation de la tension de sortie Les spcihcations des applications im-
posent en gnral une valeur minimale des ondulations de la tension de sortie.
Ici nous imposons le taux dondulation k
vo
= 0, 05 :
v
o
< k
vo
v
o
(4.16)
Densit de courant Ahn dviter des tempratures trop importantes dans le
cuivre, une densit de courant maximale de 10 A.mm
2
est impose.
Temprature de jonction des semi-conducteurs Les tempratures de jonction
des diodes et du MOSFET ne doivent pas excder 130
C ahn dassurer un fonc-

tionnement normal.
Induction magntique dans le noyau de linductance La densit de fux ma-
gntique B
M
dans le noyau de fer de linductance est borne 0,3T pour sas-
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Quantit Contrainte
CCM condition i
Lm
> 0
Output voltage ripple v
o
< k
vo
v
o
Current density J < 10 A.mm
2
Junction temperatures (T
jK
, T
jD
< 130
C
Induction magntique dans la ferrite B
M
< 0, 3T
Geometric constraint n.s < k
f
.A
Inductor temperature T
core
< 80
C
T~nir 4.11 Contraintes de loptimisation du convertisseur abaisseur
surer que linductance ne soit jamais sature pour un matriau ferrite de type
3F3.
Gomtrie de linductance La gomtrie de linductance est caractrise par
cinq variables (e lentrefer entre la bobine et le noyau de ferrite, A laire carre
de la section de bobinage, S la section du noyau, s la section dun hl de la bobine,
et n le nombre de tours dans la bobine). Ces variables sont indpendantes, mais
on doit sassurer de la faisabilit de linductance en imposant que le bobinage
sinsre dans la fentre A du noyau :
n.s < k
f
.A (4.17)
o k
f
est pris ici gal 0,75.
Temprature de linductance Pour assurer une stabilisation de la temprature
au sein de linductance, celle-ci est limite 80
C. En enet en dessous de cette

valeur critique les pertes diminuent alors quelles augmentent pour des temp-
ratures suprieures (caractristique du matriaux 3F3 du noyau magntique).
Ces contraintes sont rsumes dans les tableau 4.11.
4.3.4 Rsultats de loptimisation
Choix de la technique doptimisation Loptimisation dun hacheur abaisseur
est une tche complexe, elle implique des variables entires et relles. De plus le
cas prsent ici est celui dune optimisation multi-objectif, cest pourquoi nous
nous sommes orient vers un algorithme volutionnaire coupl une proc-
dure du type Multiobjective Optimisation Evolutionnary Approches (MOEA) et
plus particulirement de lalgorithme SGA2 [EKL] [KAST02].
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Rsultats La hgure 4.12 montre le rsultat du front de Pareto reliant les ob-
jectifs de volume en abscisse et de rendement en ordonne. Le rendement va-
rie entre 92,6% et 94,6% (convertisseur avec le moins de perte). Le volume est
lui compris entre une valeur minimale de 1, 4 10
6
mm
3
et un encombrement
maximal de 7, 6 10
6
mm
3
. Rappelons que les ralisations faisables sont toutes
situes en dessous de cette courbe. Pour dmontrer la robustesse de loptimisa-
tion, lalgorithme a t appliqu plusieurs fois. Nous obtenons les courbes de la
hgure 4.13 qui se superposent.
0 1 2 3 4 5 6 7 8
x 10
3
92
92.5
93
93.5
94
94.5
95
Buck converter external volume [m
3
]
E
f
f
i
c
i
e
n
c
y

[
%
]
FicUvr 4.12 Front de Pareto pour loptimisation multi-objectif volume
rendement
Ce front de Pareto nous permet dans un premier temps de qualiher les conh-
gurations possibles et daider au choix pour le constructeur. Sa forme trs incur-
ve nous indique quun choix logique serait de situer les paramtres autour du
point dabscisse 2 10
6
mm
3
. Mais le front de Pareto nest pas le seul outil dis-
ponible, il est aussi intressant dobserver les variations des autres paramtres
ou contraintes de loptimisation le long de cette courbe. Ces observations per-
mettent aussi de prendre en compte des contraintes non exprimes ici car dim-
ciles quantiher, telle que le cot (automobile) ou la habilit (aronautique).
Infuence du volume extrieur et du rendement Ces deux rsultats se
lisent directement sur le front de Pareto. Si le volume est le paramtre le plus
important le choix sera fait sur la partie gauche du front ; si cest le rendement
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
0 0.002 0.004 0.006 0.008 0.01
92
92.5
93
93.5
94
94.5
95
3
]
E
f
f
i
c
i
e
n
c
y

[
%
]

front 1
front 2
front 3
front 4
FicUvr 4.13 Calculs successifs du front de Pareto
on choisira plutt sur la droite. Bien sr, il faut tenir compte de la forte courbure
du front et ne pas trop pnaliser le deuxime paramtre.
Infuence des tempratures de jonction Sur la hgure 4.14, les tempra-
tures de jonction de la diode et du MOSFET sont confrontes au front de Pareto.
Dun point de vue industriel, on sait que des tempratures suprieures 100
C
pnalisent de manire importante la dure de vie des composants [EF97]. La lec-
ture de cette hgure permet donc de complter le choix sans avoir modiher les
contraintes de loptimisation en y incluant une temprature de jonction maxi-
male plus basse.
Infuence de la frquence de commutation La frquence de commutation
et le volume extrieur du systme sont troitement lis comme le montre la
hgure 4.15. Cette liaison se fait par le volume de linductance. Cette interaction
est principalement due au volume de linductance.
De la mme manire, on peut reprsenter les variations de longueur du dis-
sipateur, des paramtres de linductance et du condensateur. Comme vu pour
les autres variables et rsultats de loptimisation ces variations et tendances
peuvent elles aussi aider la dcision pour le concepteur en faisant apparatre
des contraintes non spcihes dans un premier niveau de cahier des charges.
Cet exemple dmontre que le front de Pareto seul ne peut pas donner un
rsultat optimal. Bien souvent toutes les contraintes ne peuvent tre modli-
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
0 0.002 0.004 0.006 0.008 0.01
118
120
122
124
126
128
130
T
e
m
p
e
r
a
t
u
r
e

[
C
]
3
]

0 0.002 0.004 0.006 0.008 0.01
92
92.2
92.4
92.6
92.8
93
93.2
93.4
93.6
93.8
94
94.2
94.4
94.6
94.8
95
E
f
f
i
c
i
e
n
c
y

[
%
]
T
jK
T
jd
FicUvr 4.14 Tempratures de jonction des diode et du MOSFET le long du front

de Pareto
0 0.002 0.004 0.006 0.008 0.01
1
2
3
4
5
6
7
8
9
10
x 10
4
F
r
e
q
u
e
n
c
y

[
H
z
]
3
]

0 0.002 0.004 0.006 0.008 0.01
92
92.2
92.4
92.6
92.8
93
93.2
93.4
93.6
93.8
94
94.2
94.4
94.6
94.8
95
E
f
f
i
c
i
e
n
c
y

[
%
]
F
FicUvr 4.15 Frquence de commutation le long du front de Pareto

ses (ou leur prise en compte conduirait des calculs lourds). Il faut une ana-
lyse supplmentaire sur tous les paramtres et rsultats pour aider la dcision,
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
apportant ainsi un clairage important sur la faisabilit du systme. Ceci nous
montre que loptimisation nest pas seulement un problme mathmatique de
formulation du problme puis de rsolution par des algorithmes adapts. Par cet
exemple, nous montrons que lanalyse a posteriori des rsultats doptimisation
doit rester trs critique. Pour cela il est ncessaire de continuer produire des
modles comportementaux, fonctionnels et de comprhension. Cet aspect sera
lobjet du prochain chapitre de cet expos.
4.4 Conclusion
Le thme principal de ce chapitre t de prsenter notre recherche sur
loptimisation des machines lectriques. Aprs avoir introduit la problmatique,
o nous avons montr limportance de loptimisation systme, nous avons pr-
sent deux exemples de loptimisation des lments des chanes de motorisa-
tion lectrique. Dans le premier, nous avons introduit une nouvelle technique
doptimisation dune gamme de machines suivant leur puissance. Cette nouvelle
approche permet des gains relativement importants face aux volumes pouvant
tre concerns par les fabricants. Nous avons aussi trait au cours de cet exemple
dune double optimisation sur les cots de fabrication et les rendements. Le front
de Pareto construit, permet de choisir entre linvestissement ou les cots de fa-
brication mais aussi entre les cots dutilisation.
Enhn, sur un deuxime exemple, nous avons introduit loptimisation dun
convertisseur de puissance, et ce dans le but de prendre en compte cet aspect
dans loptimisation globale des systmes.
Dans le prochain chapitre, nous nous intresserons la modlisation et la
comprhension des phnomnes physiques, en prenant comme hl conducteur la
habilit des systmes lie aux contraintes des optimisations.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
4.4 Conclusion 69
FicUvr 4.16 Amazing Stories juin 1926
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Chapitre 5
Modlisation des phnomnes
physiques
5.1 Introduction
Optimiser cest se placer en limite de habilit. Sur les exemples prcdents,
les rsultats des optimisations ont bien montr limportance des contraintes
dans lnonc du problme. Que ce soit sur loptimisation dune gamme ou dun
convertisseur de puissance, les contraintes en temprature, en rendement ou
en fonctionnement bornent les problmes poss. Il est donc impratif de com-
prendre et de poser des contraintes ralistes et compatibles avec les applications
industrielles. De plus pour optimiser, il nous faut comprendre et dvelopper des
modles permettant de simuler les points de fonctionnement des systmes. Ces
modles bien que ncessairement simples pour pouvoir tre inclus dans la pro-
cdure doptimisation doivent sappuyer sur des modlisations complexes qui
seront rduites par la suite. Pour ces deux raisons, il nous est apparu ncessaire
de poursuivre dans llaboration et la comprhension de modles physiques des
composants.
Cest ce que nous avons ralis par le biais de collaborations avec Philippe
Test du Laboratoire de Gnie lectrique de Paris (LGEP), Emmanuel Odic et
Mike Kirkpatrick du dpartement nergie de Suplec. Cette collaboration
donn lieu plusieurs thses, dont deux que nous allons prsenter plus en d-
tail. Dans la premire, nous nous sommes intresss ltude des mcanismes
pouvant donner lieu linfammation dune carte de circuit imprim.
71
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
72 Cu~vi1vr 5 : Modlisation des phnomnes physiques
5.2 tude des mcanismes dinfammation dun ma-
triau isolant en prsence dun point chaud
dorigine lectrique
La multiplication des calculateurs au sein des vhicules terrestres ou ariens
a entran une hausse de la puissance embarque et de la densit de composants
dans les circuits, augmentant ainsi la probabilit dincidents dorigine lectrique
et lectronique ; ceux-ci pouvant occasionner des dfaillances du point de vue
de la habilit du circuit mais aussi une dgradation svre du systme lectrique
et lectronique, et conduire la mise feu de celui-ci. Petit petit, larchitec-
ture lectrique et lectronique sest complexihe au point denvisager des ar-
chitectures multiplexes dans le but de rduire la quantit de cbles lectriques.
Au cours de leur fonctionnement, les systmes lectriques embarqus sont sou-
mis de fortes contraintes environnementales : mcaniques (vibrations), clima-
tiques/chimiques (variations de tempratures, prsence deau, de poussires, de
polluants divers. . .). Aussi, ils doivent subir un certain nombre de validations en
vue de garantir leur robustesse et leur durabilit. Ces validations sappuient sur
des normes dessais qui sont rgulirement mises jour en fonction de lvolu-
tion des technologies et du retour dexprience, particulirement des construc-
teurs automobiles. Parmi les paramtres prendre en considration lors de la
validation hgurent les phnomnes dchaunement, au niveau notamment de la
carte lectronique et plus prcisment du circuit imprim. Lchaunement dun
circuit imprim peut enectivement conduire dinrentes consquences dont des
dysfonctionnements, des pannes, voire, dans le pire des cas, une infammation
de la carte. Lors de sa thse sinscrivant dans le cadre dune action de prven-
tion et de matrise des risques, mademoiselle Michle Nsoumbi sest intresse
ltude lchaunement dun circuit imprim dans le but den apprhender les
modes dapparition et den valuer les consquences, et de mettre en place,
terme, des dispositifs ou rgles de conception prventifs. Ainsi lobjectif de ce
travail est de contribuer lidentihcation et la comprhension des mcanismes
physiques lorigine dun chaunement et qui conduisent linfammation dun
circuit imprim. Si ce sujet parait, au premier abord, simple et comprhensible
par tous, il savre quil na que peu ou pas t tudi et quil met en uvre de
nombreux phnomnes physiques. La complexit rside dans lamalgame de ces
phnomnes lectriques, thermiques, chimiques et dans la dimcult de les
isoler pour comprendre leurs enchanement et interactions. Dans cette optique,
la prsente tude a considr quune carte lectronique se composait principale-
ment de deux lments savoir un substrat (matriau composite isolant) et des
pistes conductrices de cuivre confront diverses situations dchaunement.
Lapproche choisie pour ce travail repose sur lhypothse quun point chaud
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
5.2 tude des mcanismes dinfammation dun matriau isolant en prsence
dun point chaud dorigine lectrique 73
dorigine lectrique, rsultant dune dfaillance sur un composant et/ou une
piste, est susceptible damorcer une infammation. Une premire tape a consist
caractriser les matriaux tudis dun point de vue lectrique, et dans une
certaine mesure, dun point de vue chimique. Ainsi, ahn de prdire lvolution
du caractre isolant du substrat en prsence dun point chaud, des mesures de
rsistivit lectrique en volume ont t menes sur trois circuits imprims de d-
nomination FR-4. Les essais ont montr que celle-ci dcroissait fortement avec la
temprature, en prsentant des disparits de comportement selon la provenance
des chantillons. Les dinrences de comportement observes entre les trois ma-
triaux ont pu tre attribues, dune part aux proprits lectriques des consti-
tuants de la matrice poxy, et dautre part, au caractre anisotrope du substrat
notamment d au tissage des hbres de verre. Pour une grande majorit des ma-
triaux, il sest avr que plus la fraction volumique de hbres tait faible, plus
rapide tait la dcroissance de la rsistivit lectrique avec la temprature. Les
mesures de caractrisation chimique, menes haute temprature (300
C ), ont
notamment permis de dceler la prsence dhydrogne et dactylne. Ces es-
pces chimiques, extrmement infammables, seraient susceptibles de jouer un
rle dans linitiation et/ou la propagation dune infammation.
La seconde tape a consist tudier la dgradation dun circuit imprim
simple tel que reprsent schmatiquement sur la hgure 5.1 et en photo aprs
infammation sur la hgure 5.2.
FicUvr 5.1 Schma du circuit imprim tudi
FicUvr 5.2 Photo recto (a) et verso (b) du circuit imprim tudi aprs infam-
mation
Contrairement aux dfaillances de habilit, couramment observes, lies
des dfauts apparaissant sur de longues dures (plusieurs centaines ou milliers
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
dheures) et pouvant mettre en jeu des phnomnes dlectro-migration, nous
nous sommes intresss des dfauts lis des surintensits conduisant rapi-
dement des tempratures de piste et de substrat leves. La dfaillance a t
exprimentalement modlise par :
un dfaut gomtrique contrl de la piste,
une surcharge en courant, elle aussi contrle, conduisant ainsi un
chaunage progressif du circuit.
Un suivi temporel des grandeurs lectriques a t ralis. On peut voir sur la
hgure 5.3, lvolution du courant dalimentation (entre les deux extrmits de la
piste) et des courants de fuite (entre la piste et la contre-lectrode au verso, cest
dire un courant traversant le matriau du circuit) en fonction du temps et de
la temprature du substrat pour deux paisseurs (0,4mm et 1,6mm).
De plus, au cours de lexprience, des cartographies thermiques de la surface
du circuit imprim (substrat et piste) ou de lpaisseur de la plaque au niveau du
dfaut ont t ralises laide dune camra thermique. La hgure 5.4 illustre un
rsultat de cartographie thermique ainsi que les prohls de tempratures le long
de la piste.
Enhn, une camra rapide a permis dobserver la rupture de la piste et, le cas
chant, les dparts de feu ainsi que leurs propagations sur la carte (voir hgure
5.5 pour un exemple.)
Une tude prliminaire du comportement dun circuit imprim face un
chaunement caractris a permis didentiher divers paramtres et conditions
opratoires favorables linfammation du circuit imprim. Parmi eux, on peut
citer :
la dynamique daugmentation de lintensit du courant dans la piste,
lpaisseur des cartes,
la nature des cartes,
la prsence dun conducteur polaris proximit de la piste,
la gomtrie de la piste.
La suite de cette tude a port sur les phnomnes lectriques aggravant ra-
pidement la dgradation du circuit imprim. Les courants de fuite circulant dans
le substrat entre deux conducteurs polariss ont ainsi t tudis. Dans cet ob-
jectif, une carte double face a t utilise, comportant sur le recto, une piste avec
un rtrcissement en son centre, modlisant par un point chaud la dfaillance
dun composant lectronique sur une carte et sur le verso, une contre-lectrode
polarise. Parmi les dinrentes conhgurations dlectrodes testes, seule celle
nomme contre- lectrode tendue a fait lobjet dune tude approfondie en
raison de sa correspondance avec les cartes lectroniques usuelles. Dans cette
conhguration, deux types de matriaux et trois paisseurs de carte ont t utili-
ss. En alimentant la piste avec une intensit de courant progressive, pouvant at-
teindre jusqu quatre fois lintensit nominale (hgure 5.3), trois comportements
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
FicUvr 5.3 volution du courant dalimentation et du courant de fuite en fonc-
tion du temps (haut) et de la temprature du PCB (bas)
ont principalement t observs lissue des essais et ont conduit dhnir dif-
frents types de rupture de piste savoir :
la simple coupure de piste ;
la coupure de piste accompagne dune infammation visible (quelques se-
condes) ;
la coupure de piste accompagne dune infammation visible et entretenue
(dure suprieure 10 secondes).
Le courant de fuite, circulant au travers dun substrat entre les deux pistes po-
larises, a t identih comme paramtre indicateur de ltat de dtrioration
du matriau. En enet, au cours dun chaunement croissant, trois phases se dis-
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
FicUvr 5.4 Cartographie de temprature et prohls le long de la piste
tinguent principalement :
une premire phase durant laquelle lintensit du courant de fuite est
faible (< 10 A) ;
une deuxime phase, dite de conduction, au cours de laquelle on observe
une augmentation signihcative du courant de fuite (< 1 mA) ;
et une dernire phase demballement lissue de laquelle interviennent la
coupure et linfammation de la carte.
Les rsultats exprimentaux enrichis par les cartographies thermiques de surface
et les mesures de rsistivit lectrique en fonction de la temprature ont t
confronts aux rsultats dune simulation numrique par lments hnis ralise
laide du logiciel COMSOL.
La Figure 5.6 illustre la gomtrie du modle considr pour les simulations.
Cette gomtrie se compose de 3 volumes :
au premier plan, la piste est dhnie par trois sections (2mm, 1,5mm et
2mm) et par son paisseur (33m) ;
au second plan, le substrat dune paisseur de 0,4mm, dune longueur de
60mm et dune largeur de 25mm, sur lequel repose la piste ;
le plan de masse (33m dpaisseur), situ au verso et occupant la mme
surface que le substrat, qui nest pas visible sur ce schma en raison de sa
faible paisseur par rapport au substrat.
Une tude temporelle couplant les phnomnes lectriques et thermiques
a t mene, un exemple de calcul est reprsent sur la hgure 5.7. Lvolution
de la cartographie de temprature simule montre lchaunement localis au
niveau du dfaut de la piste, mais aussi la relativement bonne rpartition de la
temprature le long de la piste.
En complment on trouvera sur la hgure 5.8 une comparaison des temp-
ratures exprimentales et simules le long dun prohl transversal la piste au
niveau du dfaut qui montre une cohrence des rsultats. Cette bonne adqua-
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
FicUvr 5.5 Images enregistres par la camra rapide au cours de linfammation
dun circuit imprim de 0,8mm
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
FicUvr 5.6 Gomtrie utilise pour le modle numrique
FicUvr 5.7 Distribution de la temprature au cours du temps dans un substrat
z = 0,2 mm pour une intensit de courant de chaune comprise entre 10A et 20A
et un coemcient dchange thermique gal 8W.m
2
.K
1
tion entre modle numrique et exprimental nest par contre pas respecte sur
un temps plus long comme le montre la hgure 5.9. Ceci met en vidence lexis-
tence de phnomnes non pris en compte dans la simulation numrique et qui
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
FicUvr 5.8 Comparaison entre les prohls de temprature exprimentaux et
ceux calculs numriquement transversalement au dfaut, pour une valeur de h
gale 8W/m
2
/K
conduisent llvation de temprature au sein du matriau.
Cette confrontation a permis de valider un modle physique bas sur le ca-
ractre thermique de la dgradation. La deuxime phase a ainsi pu tre interpr-
te par une chute de la rsistivit lectrique dans le volume du matriau due
laugmentation de temprature (phnomne non pris en compte dans le modle
contrairement la dpendance de la rsistivit du cuivre en fonction de la temp-
rature). Lemballement juste avant la rupture de la piste semble li lexistence
dun dfaut trs localis permettant le passage dun courant de densit leve
susceptible de crer les conditions propices lapparition dune infammation.
De plus, lissue de lemballement et aprs rupture de la piste, la persistance du
courant de fuite entre lextrmit de la piste potentiel nul et la contre-lectrode
tendue (polarise) apparat comme un facteur aggravant permettant lentretien
et la propagation de linfammation, dune part en maintenant une zone chaude
en surface du circuit, et dautre part, en contribuant lmission de composs
organiques volatils et infammables susceptibles dentretenir et de propager le
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
FicUvr 5.9 Courbes dvolution en fonction du temps de chaune de la temp-
rature la surface dun substrat pour un courant de chaune initial gal 4 A -
Comparaison entre les rsultats exprimentaux et le modle numrique
feu.
Le modle numrique nous a aussi permis davoir accs une cartographie
de la rpartition du courant de fuite au sein du matriaux sur la hgure 5.10. On
note que la densit de courant est concentre autour du dfaut, ce qui corrobore
les hypothses faites lors de ce travail.
Le courant de fuite est donc apparu non seulement comme un indicateur de
ltat du circuit imprim, mais aussi comme un facteur aggravant dune part, en
facilitant le dpart de feu et dautre part, en contribuant son entretien et sa
propagation sur le circuit imprim.
laide des cartographies thermiques, une estimation de la surface chauf-
fe, travers laquelle le courant de fuite tait susceptible de circuler, a montr
que la contribution de celui-ci lchaunement global de la carte tait ngli-
geable au cours de cette phase. La modlisation numrique en 3-D de lchauf-
fement dun circuit imprim double face, dune paisseur de 0,4mm, a permis de
calculer la distribution de la tension et de la temprature lintrieur de celui-ci.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
FicUvr 5.10 Modle dune distribution de la densit de courant selon laxe z
un instant t dans un substrat
Une bonne corrlation a t obtenue entre les rsultats exprimentaux et ceux
issus du modle numrique. La transition menant la dernire phase, attribue
au dcollement de la piste du substrat, a galement pu tre modlise. Au cours
de cette phase, un emballement de la temprature est visualis, aussi bien la
surface de la carte que dans le volume, o des points chauds apparaissent.
ce moment-l, lintensit du courant de fuite, trs leve, est reprsentative de
ltat de dgradation du substrat. Dans ces conditions, lenregistrement du cou-
rant de fuite peut tre considr comme un paramtre indicateur de ltat de
dtrioration du matriau. Dans le cas de cartes double face, la persistance du
courant de fuite apparat comme un facteur aggravant permettant lentretien de
linfammation sur la carte.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
FicUvr 5.11 Amazing Stories avril 1926
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
5.3 Physique des dcharges lectriques 83
5.3 Physique des dcharges lectriques
5.3.1 Description et modlisation des phnomnes physiques
des dcharges lectriques dans les plasmas froids
Introduction
Nous avons tudi dans la section prcdente les infammations dues des
surintensits et des courants de fuite au sein du matriau. Dans la prsente sec-
tion, nous allons aborder un sujet plus physique, mais tout aussi important dans
llaboration des composants de systmes optimiser. En enet loptimisation qui
conduit rduire les masses ou les volumes, rduit aussi les distances entre les
lments conducteurs. La rduction de cette distance augmente ainsi les champs
lectriques et les probabilits de prsence de dcharge partielle ou pire darc lec-
trique. Ces derniers peuvent tre destructeurs pour les conducteurs et conduire
un vieillissement prmatur des composants. Cest pourquoi dans le cadre des
problmatiques de habilit des systmes tudies au dpartement nergie de
Suplec, nous nous sommes intresss aux conditions et aux mcanismes de
passage dun arc entre deux lectrodes.
Le terme de claquage, souvent utilis en franais comme synonyme darc
lectrique, exprime le passage dun courant important dans un milieu initiale-
ment non-conducteur. Larc lectrique est caractris par un fort courant et une
faible tension aux lectrodes, il est prcd dune phase de pr-ionisation de les-
pace interlectrode. Celle-ci peut prendre plusieurs formes, soit une dcharge de
type hlamentaire aussi appele streamer ou dinuse de type glow, en fonction des
conditions de pression, de temprature, de champ lectrique. . .
Dans le cadre des tudes des dcharges hlamentaire prdisruptives, nous
avons continu les travaux de modlisation numrique des streamers ; en tra-
vaillant dune part sur la partie mathmatique et algorithmique des schmas
numriques de rsolution, mais aussi sur les rsultats physiques donns par les
codes de calculs.
Physique des dcharges prdisruptives
Lair est un composant naturel de notre environnement, cest pourquoi de
nombreuses tudes ont t ralises sur les plasmas dans lair pression at-
mosphrique ainsi que sur leurs applications. Ltude gnrale et synthtique
prsente ici, porte sur la description phnomnologique et la modlisation nu-
mrique des tats possibles des dcharges dans lair ainsi que des applications
possibles. Elle servira de plus dintroduction ltude de la section suivante sur
ltude des mcanismes de transition glow-arc dans largon faible pression.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
La conduction lectrique dans un gaz est un phnomne hautement dpen-
dant de la pression de celui-ci. La raison principale est la modihcation du libre
parcourt moyen des lectrons avec la pression. 10
1
Torr le libre parcourt
moyen dun lectron est de lordre de 30 cm, soit de lordre de grandeur dune
distance interlectrode typiquement tudie en laboratoire. 10 bar, il se r-
duit 4m, soit infrieur a une distance interlectrode typique. Lapparence
visuelle et lespace occup par une dcharge lectrique entre deux lectrodes est
de ce fait largement dpendante de la pression. faible pression, une dcharge
lectrique va apparaitre comme tant tendue spatialement, on parle alors de
dcharge dinuse ou glow. A contrario, a haute pression, les dcharges sont des
phnomnes locaux conhns spatialement. Pour toutes les gammes de pression,
un autre paramtre primordial est lnergie gagne entre deux collisions qui d-
terminera le type de ces dernires et le comportement macroscopique de la d-
charge. Lnergie gagne par les lectrons est au premier ordre gal au produit
du champ lectrique qui induit le dplacement lectronique par le libre par-
court moyen. Cest pourquoi il est habituel de caractriser les dcharges par le
rapport
E
N
du champ lectrique E par la densit N des espces lourdes (
1
N
tant
proportionnel au libre parcourt moyen). Ce rapport tant appel champ rduit.
premire vue, une rgle de similarit peut sappliquer ici : les phnomnes
basse pression peuvent tre retrouvs pour toute autre pression tant que le
champ rduit et la temprature restent constants. Mais cette rgle de simila-
rit, ne prend pas en compte des enets locaux de chaunage du gaz par le cou-
rant de dcharge. Llvation de temprature qui en dcoule est de plus en plus
concentr au centre de la dcharge pour des pressions leves car la conduction
thermique est de moins en moins emcace. Cest pourquoi pour des pressions
de plusieurs atmosphres, une dcharge darc lectrique transitoire ou continue
apparait ds lapplication dune tension seuil. Larc lectrique est caractris par
une temprature de plasma leve (de lordre de 10
3
10
4
K et par un quilibre
entre les tempratures des ions et des lectrons (tat dnomm LTE ou quilibre
thermodynamique local).
pression atmosphrique les tudes prcdentes ont montres que larc
lectrique tait toujours prcd dune phase dionisation aussi appel d-
charge prdisruptive prparant le passage de celui-ci. Durant cette phase dio-
nisation qui peut prendre plusieurs formes, la temprature des lectrons est
plus importante que celle des ions ou des atomes/molcules du gaz. Le plasma
ainsi obtenu est dit dans un tat NTP (non thermal plasma). Lavantage des NTP
est la grande ractivit chimique induite par la forte temprature des lectrons
face celle des neutres ou des ions. Les collisions lectroniques avec les par-
ticules lourdes sont frquentes et produisent des espces ractives (collisions
innlastiques). De nombreuses applications environnementales sont bases sur
la grande ractivit chimique de ces plasmas NTP, on peut citer sans tre ex-
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
haustif le traitement de leau, la production dozone, la strilisation ou dconta-
mination, la dpolution et lammlioration de la combustion. Nanmoins luti-
lisation des ces dcharges impose un moyen de controler leur tendance natu-
relle passer larc lectrique, il est alors ncessaire de bien comprendre les
mcanismes physiques lis au dveloppement de cette phase dionisation. Dans
un premier temps, nous ferons une description des dinrentes structures que
peuvent prendre ces phases dionisation, puis nous prsenterons les modles
utiliss ainsi que les rsultats de simulations.
Dveloppement spatial et structure de la phase dionisation
On peut classer les phases de prionisation avant le passage larc dun es-
pace interlectrode en deux types de dveloppement spatial. Le premier est un
mode dinus avec une apparence quasi uniforme. On se rferre cette extension
spatiale comme tant une dcharge glow, dinuse ou corrona. Un exemple dune
telle dcharge est prsente sur la photographie de la hgure 5.16 o une dcharge
glow est obtenue dans largon faible pression.
Le second type est une dcharge hlamentaire, fortement non-uniforme dans
lespace et prsentant un aspect transitoire et des branchements importants. Ces
dcharges sont appelles en gnral streamers, dards ou hlamentaires et sont g-
nralement obtenues haute pression dans les gaz lectrongatifs comme lair
et particulirement avec des lectrodes fort rayon de courbure. Ce paragraphe
traitera plus particulirement de ce type de dcharge, la section 5.3.2 sintres-
sera plus particulirement la transition entre un glow et un arc lectrique dans
largon faible pression.
Au dessous de la tension seuil de claquage Dans une conhguration dlec-
trodes uniformes de type plan-plan, tant que la tension applique reste en des-
sous de la tension seuil de claquage (passage larc), la thorie classique de
Townsend base sur une suite davalanches lectroniques sumt expliquer les
phnomnes. Les avalanches lectroniques successives sont induites par des
lectrons secondaires issus de la premire avalanche (par exemple par enet dar-
rachement dlectrons la cathode par des ions). On parle dune dcharge auto-
maintenue car le courant reste continu indhniment.
Dans une conhguration non uniforme, par exemple une pointe faisant face
un plan, la rgion dionisation appelle glow est situe proche de llectrode
faible rayon de courbure ou le champ lectrique est le plus intense. Une telle
conhguration est donne sur la rgion A de la hgure 5.12 pour un systme pointe
positive-plan dans lair.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Plasma Phys. Control. Fusion 51 (2009) 124002 E Marode et al
Figure 1. Coronas: left, sketch of the glow region (A) versus streamers region (B). Top right:
the dc glow current, and the increased total current due to streamer. Bottom right: the repetitive
streamers intermittent voltage and current. Solid lines, under spark threshold; dotted lines above
spark threshold with insets for details.
as parallel plate conguration and non-electronegative gas, hoping to keep complete control of
all parameters and discharge phenomena. They considered that the Loeb [10] and Meek and
Craggs [11] approach, in air with various gap congurations, was not sufciently mastered,
too complicated, just able to lead to descriptive observations rather than scientic results.
In non-uniform gaps, indeed, strong geometrical eld gradients around curved electrodes
and negative ion formation lead to various complicated discharge shapes, surrounding such
highly stressed electrodes, following the way Loeb liked to qualify such gaps. The generic
name coronas encompassing these various discharge manifestations was then introduced.
They may be diffuse or lamentary, with uctuations in time and space, more or less statistically
reproducible.
Both approaches however, have been the roots of the actual concept of high pressure gas
discharge and NTPs. In the parallel gap conguration, as far as the gap voltage remains near
the threshold of the spark breakdown, the Townsend theory of ionization growth does well.
Each electron avalanche is followed by another avalanche through secondary effects from
the rst one, the secondary effect being, for example, ions of the rst avalanche providing the
secondary electron by hitting the cathode. The discharge current may thus endlessly continue
and the discharge is said to be self-sustained. Such a theory is the physical basis of a number
of uniform and largely spread discharges often qualied as glows. For non-uniform gaps,
the glow ionization region is conned near the stressed electrode where the geometrical eld
is large enough, as seen, for example, in the case of a positive point-to-plane discharge in
air, gure 1, region A. Region B is the ion drift region with no ion production. As explained
in [12, 13] the eld distribution is mainly space charge controlled while the current is limited
by the ions drift.
4
FicUvr 5.12 Description dune dcharge pointe-plan
Au dessus de la tension seuil de claquage Si la tension entre les lectrodes
dpasse la tension seuil de claquage, la morphologie de la dcharge change ra-
dicalement et ceci pour les deux conhgurations dlectrodes mentionnes (uni-
formes et non-uniformes). Dans le cas dune gomtrie plane, on remarque deux
phnomnes propagatifs hlamentaires partant de lespace interlectrodes et se
dirigeant chacune vers une lectrode. Ces phnomnes sont associs des ondes
dionisation appels streamer, il est largement admis celles-ci sont dues la mo-
dihcation du champ lectrique cause par la charge despace laisse par les ava-
lanches prcdentes (ions relatiment moins mobiles que les lectrons). Londe
dionisation dirige vers lanode appelle en anglais anode directed streamer
avance dans le sens de drive des lectrons dans le champ lectrique, mais celle
dirige vers la cathode (cathode directed streamer ou streamer positif) avance
une vitesse comparable celle des lectrons dans un sens oppos celle impose
par le champ lectrique. Il est donc ncessaire davoir un mcanisme crant des
lectrons en amont du streamer (entre celui-ci et la cathode) pour provoquer de
nouvelles avalanches. Des thories bases sur la dinusion arrire des lectrons
ont t proposes, mais ncessitent des tempratures lectroniques largement
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
suprieures celle rencontres dans la dcharge. Il est maintenant commun-
ment admis que cet lectron secondaire est cr par un phnomne de photoio-
nisation partir dun photon mis partir des avalanches prcdentes.
Dans le cas dun systme non-uniforme pointe posivite - plan, un seul strea-
mer se dplace depuis la pointe positive vers le plan. Sur la hgure 5.12, le hlament
positif (C) merge de la zone dite de glow (A) et parcourt lespace interlectrode.
Sur cette hgure, a t aussi reprsent lvolution du courant pour des tensions
de lordre de la tension seuil de claquage (courbes en bleu continues). Dans ce
cas, le streamer est un phnomne rptitif et en fonction des conditions ex-
primentales, on peut obtenir une rptition sporadique et spontane ou une
rptition controle pouvant aller jusqu des frquences de plusieurs dizaines
de kiloHerz. Dans le cas dun streamer rptif, on observe une volution du cou-
rant prsentant une augmentation et un pic au moment o le streamer atteint
la cathode, puis une dcroissance lente jusqu rptition de londe dionisation
suivante.
Si la tension est suprieure la tension de claquage (courbe en rouge et en
pointills), il va se former aprs le passage du streamer un arc lectrique. Dans
ce cas le premier pic est suivi dun second beaucoup plus intense. Lchelle de
temps entre les deux phnomnes est de lordre de 50ns quelques dizaines de
s. Cest durant ce temps caractristique que le plasma laiss par le passage du
streamer est dans un tat NTP.
Distribution de la densit dnergie jE
Les mcanismes aux sein des NTP sont largement caractriss par la dis-
tribution de la densit dnergie j.E, produit de la densit de courant par le
champ lectrique. Cette densit dnergie peut tre dcompose en trois parties
lies respectivement aux ions positifs, aux ions ngatifs et aux lectrons. Cette
dernire part tant prpondrante. Cette energie va se dissiper sous plusieurs
formes : chaunage du gaz et des lectrodes (principalement pour les ions), mais
surtout sous forme de collisions inlastiques qui participeront terme elles aussi
au chaunage du gaz. Cet chaunement global est primordial pour expliquer la
formation de larc lectrique. cause de cet lvation de temprature, le canal
de dcharge subit une expansion, la densit N dcroit et le champ rduit E/N
(quivalent lnergie moyenne des lectrons) augmente. De plus les processus
chimiques sont aussi infuencs au travers de la loi dArrhenius. Il est donc n-
cessaire dvaluer hnement la distribution des lectrons et des espces charges
au sein du canal conducteur. La densit des lectrons dpend fortement du rayon
de ce celui-ci, en enet plus la dcharge sera constrainte dans un espace rduit,
plus la densit des lectrons sera importante ainsi que le densit de courant.
Par consquent, plus le rayon sera faible, plus llvation de temprature et de
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
pression sera importante.
Pour rpondre a ces questions et tablir une cartographie des distributions
des espces en prsence dans la dcharge, une modlisation hne du plasma et
des mcanismes de propagation du streamer est ncessaire. Cest pourquoi nous
nous sommes intresss la simulation en deux dimensions dune dcharge
pointe positive plan. Notre travail cest concentr sur lanalyse des rsultats de
codes de calculs dj existant et sur lammlioration des algorithmes numriques
mis en place notamment par des travaux exploratoires sur des mthodes mlants
des lments hnis et mthode des caractristiques (ELLAM).
Modlisation de londe dionisation streamer dans lair
Le modle utilis dans cette tude se base sur un jeu dquations fuides hy-
drodynamiques. Pour chaque espce, une quation de convection-dinusion en
coordonnes cylindriques est crite en se basant sur le modle de lquation 5.1 :
n
t
+

x
nw
x
D
x
n
x
+
1
r
nw
r
D
r
n
r
= S (5.1)
o, pour chaque espce considre : n est la densit, w
x
et w
r
les vitesses
de conduction axiale et radiale, D
x
et D
r
les coemcients de dinusion axiale et
radiale. S est un terme source de cration et destruction. Il prend en compte
pour les lectrons et les ions positifs de lionisation (via le premier coemcient
de Townsend) ; et de la photoionisation et pour les lectrons et les ions ngatifs
de lattachement. Ce modle sert de base de nombreux travaux sur la simu-
lation dun streamer, en gnral cest sur la modlisation de la photoionisation
que ceux-ci dinrent. Ici une approche sappuyant sur une quation intgrale du
type de lquation 5.2a t choisi. Cette quation est simplihe par lhypothse
dEddington se basant sur la supposition que le rayonnement est faiblement ani-
sotropique proche de son lieu dmission (ici la tte du streamer).
S
ph
= c
ph
nu
r,
, t
dd (5.2)
On dhnit

comme la direction dun photon la position r au temps t, c est
la clrit de la lumire,
ph
est le coemcient de photoionisation dpendant de

la frquence du photon considr.
nu
r,
, t
est la distribution de la densit

de photons la frquence .
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Ce systme dquation est ferm par la connaissance du champe lectrique
E donn par la rsolution de lquation de Poisson :
E =

0
(5.3)
o est la densit nette de charge.
Ce modle complet est rsolu par lapplication dune mthode de pas frac-
tionns qui nous permet de traiter les quations dinremment suivant leur type :
dinrences hnies pour la dinusion et caractristiques pour la convection. Las-
pect volutif et propagatif du streamer impose la cration dun maillage mobile
suivant le dplacement de la dcharge dans lespace et le temps. Les travaux ac-
tuels sur ce modle portent sur lapplication de mthodes de rsolution mixtes
par lments hnis et caractristiques.
Comparaison des rsulats de simulation avec les expriences
FicUvr 5.13 Comparaison du courant de dcharge et de lmission lumineuse
la camra streak, simulation gauche et expriences droite
La hgure 5.13 montre les rsultats de la simulation en comparaison avec
ceux issus de lexprience. Sur la partie haute de la hgure, droite est reporte
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
une photographie streak dune dcharge. Chaque image reprsente un instant
donn (toute les 3 nanosecondes) lmission lumineuse de la dcharge et le rayon
total est de 600 m. Sur la mme hgure, droite en haut pour les mmes ins-
tants est report les rsultats de la simulation pour lmission photonique. Les
rsultats de la simulation sont en bon accord avec les rsultats exprimentaux,
notamment sur le temps de propagation du streamer (40 ns) pour la jonction
la cathode et donc sur la vitesse de londe dionisation. De plus, aprs ltablis-
sement du canal conducteur, une zone lumineuse la cathode aparait dans les
deux cas appel un streamer secondaire ou secondary streamer en anglais.
Sur la mme hgure, est report gauche le courant de dcharge simul, com-
par droite avec celui issu de lexprience. Si la forme gnrale semble la mme,
on peut remarquer quelques dinrences. Le premier pic de courant simul, nap-
parait pas sur le courant exprimental. En enet dans le cas expriemental, la
mesure est ralise sur un streamer stablissant dans un milieu prionis par le
passage prcdent dun autre streamer. Dans le cas de la simulation, ce premier
pic de courant est reprsentatif de ltablissement de la zone de glow autour de
la pointe et des premires avalanches lectronique dans un milieu non ionis.
Considrations sur le rayon du streamer
Nous lavons prcis prcdemment, le rayon du canal conducteur form par
le passage de londe dionisation est trs importante pour laspect dexpansion
hydrodynamique et le possible passage larc lectrique aprs la phase de prio-
nisation du gaz. Il est donc ncessaire davoir une bonne estimation de la valeur
de ce rayon. Exprimentalement les mesures, stalent suivant les techniques
utilises, sur une vaste plage entre 10m et 170m. Les rsultats de la simula-
tion donnent pour leur part un rayon de lordre de 100m, cest dire dans la
fourchette haute de lestimation exprimentale. Pour rduire ce rayon simul,
nous nous sommes intresss aux valeurs des coehcients de la photoionisation
et des mcanismes pris en compte. Lajout de lmission induite dans le calcul du
terme de photoionisation na pas permis lui seul la rduction de la valeur du
rayon de dcharge simul. Par contre, nous nous sommes aperu que le rayon
variait en fonction du coemcient dabsorption des photons mis par la tte de
streamer. Pour de faibles valeur de ce coemcient dabsorption, le rayon obtenu
par simulation est proche de 50m. Ceci montre la ncessit davoir des donnes
prcises des coehcients et des mcanismes pris en compte.
De plus, le modle bas ici sur des quations fuides une tendance intrin-
sque crer des densits dont la somme est infrieure une particule. Si cela
nest pas pnalisant lors de la phase de propagation du streamer, cela doit poser
des problmes lors de la formation de la zone de glow et de ltablissement de la
charge despace au dmarrage du streamer la pointe. En enet, ce moment trs
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
peu dlectrons sont prsent dans le gaz (un seul au dpart du streamer), et la
simulation cr des fractions dlectrons dans les zones de plus fort champ lec-
trique, cest dire sur la surface de la pointe. On assiste alors un talement de
la distribution de la densit des lectrons autour de la pointe, talement spatial
qui sera projet dans lespace interlectrode lors de la propagation et sera donc
une des causes de la valeur du rayon. De plus ce phnomne purement num-
rique et simulatoire peut aussi expliquer le premier pic de courant relativement
important dans la simulation.
Dernirement, des mesures exprimentales ont dtermin un rayon du cur
de la dcharge denviron 10m, montrant lcart avec les modlisations. Des
questions se posent naturellement sur la dtermination exprimentale et la dh-
nition du rayon de dcharge. En enet, faut il se baser sur des mesures lumineuses
par spectroscopie ou plutot sur des mesures de la distribution de la densit des
lectrons. Nanmoins un faible rayon de dcharge suporte lhypothse que la
transformation de celle-ci en un arc lectrique est principalement due une ex-
pansion hydrodynamique et une augmentation du champ rduit
E
N
.
Extensions du modle lectrique 2D
La modlisation 2D prsente ici traite du dveloppement du streamer en
coordonnes cylindriques axiales et radiales. Il ne tient pas encore compte des
phnomnes dexpansion hydrodynamique ultrieurs au passage de la dcharge
au sein de lespace interlectrode. Cette extension est prvue dans un cadre plus
vaste de ralisation dun logiciel de simulation comportant une phase de simu-
lation du passage de la dcharge dans un gaz qui se couplerait des codes de
simulation dexpansion hydrodynamique et de cintique chimique.
Nanmoins un code de calcul une dimension radial t dj mis au point
partir des rsultats du modle 2D prsent ici. Sa ralisation sappuie sur un
dcouplage de lquation de convection-dinusion en trois quations de conser-
vation (masse, quantit de mouvement et nergie). Ce modle est reli partir
des donnes de temprature et de pression une simulation de la dynamique
chimique des ractions.
Un autre aspect des dcharges prdisruptives ou NTP sur lequel nous avons
travaill est ltude des mcanismes de transision des dcharges non hlamen-
taires ou glow faible pression dans largon et fera lobjet du paragraphe suivant.
5.3.2 tude des mcanismes de transitions glow-arc dans lar-
gon faible pression
Ltude de la transition dune dcharge luminescente ou glow (forte tension,
faible intensit du courant) en un arc lectrique (faible tension, forte intensit
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
du courant) prsente un intrt acadmique majeur tant du point de la physique
des dcharges (incluant larc lectrique) que du point de vue de linteraction
plasma/surface. En enet les phnomnes physiques qui rgissent la physique de
la dcharge lectrique sont sensiblement dinrents de ceux qui gouvernent la
physique de larc lectrique et jusqu prsent les mcanismes conduisant la
transition ne sont pas bien connus. Ltude de la transition dun mode lautre
prsente donc un intrt scientihque particulier.
Ltude de la transition glow vers arc et arc vers glow est un sujet largement
tudi depuis plus de 40 ans, suivant des approches la fois exprimentales et
de simulations. Ces tudes sont encourages par de nombreuses applications in-
dustrielles comme le traitement de surface ou la prvention du passage larc
pour les dcharges dans le CO
2
(lasers CO
2
). Ltude des dcharges glow est
trs ancienne et de nombreux travaux ont t raliss dans dinrents gaz (air
pression atmosphrique, azote, oxygne, air et argon). Plus prcisment des re-
cherches importantes ont t faites sur les transitions glow-arc dans des gammes
de pression allant de quelques Torrs la pression atmosphrique dans des gaz
rares. Ltude que nous prsentons ici, et qui est lobjet de la thse de Romarick
Landfried co-encadre avec Emmanuel Odic, Mike Kirkpatrick (Suplec dpar-
tement nergie) et Philippe Test du Laboratoire de Gnie lectrique de Paris,
sintresse aux transitions dans largon faible pression (jusqu 100mbar).
Avant de raliser une tude par simulation numrique, il est toujours n-
cessaire de comprendre les mcanismes mis en jeu dans le processus physique,
ainsi que de caractriser les paramtres dune telle simulation. Cest pourquoi la
premire partie de cette tude porte sur un travail exploratoire. Aprs avoir
prsent rapidement le dispositif exprimental, nous aborderons les rsultats et
les hypothses que nous en avons dduites sur les mcanismes de transition. Ces
hypothses serviront dans un futur proche alimenter un modle numrique de
la dcharge.
Dispositif exprimental
La cellule de dcharge est constitue de deux lectrodes cylindriques co-
axiales en tungstne. Elles sont places lintrieur dune enceinte hermtique
(hgure 5.14) en acier inox emplie dargon une pression de 100mbar. La distance
interlectrode est variable et peut tre ajuste entre 5mm et 45mm. Lalimenta-
tion lectrique dont le schma est reprsent sur la hgure 5.15 mise en jeu
dans le dispositif de transition est dcompose en deux parties.
Dune part une composante continue permettant de gnrer une dcharge
glow de type luminescente telle que celle photographie en pause longue sur
la hgure 5.16. Le courant de dcharge est maintenu constant via lajustement
de la tension par une capacit de 10F qui se dcharge travers une rsistance
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
FicUvr 5.14 Photographie de lenceinte de dcharge
FicUvr 5.15 Schma lectrique de lalimentation utilise
de 180k. La deuxime partie permet de superposer au courant de dcharge
un chelon de courant dune amplitude variable dans lespace intrelectrode.
Cet chelon permet de faire transiter la dcharge luminescente glow vers un
arc photographi en pause longue sur la hgure 5.17. Lamplitude de ce pulse
de forme rectangulaire peut tre ajuste jusqu 400mA pour une dure de 2ms
grce la dcharge dune capacit de 10F travers une rsistance de 18k
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
relie lanode par un interrupteur rapide. Dans ces conditions la seconde partie
de lalimentation ( droite sur la hgure 5.15) agit comme une source de courant.
Lorsque linterrupteur est ferm, le courant rsultant est hx par la tension du
condensateur et la valeur de la rsistance (18k).
En plus de mesures lectriques (courant et tension de dcharge), une camra
rapide est couple au dispositif et permet de raliser des observations visuelles
de la transition.
FicUvr 5.16 Photographie en pause longue dune dcharge luminescente glow
Transitions irrversibles
Illustrons tout dabord le cas simple dune transition irrversible, ce qui
nous permettra de poser quelques dhnitions phnomnologiques de la dcharge
glow luminescente et de larc. Lorsque la premire partie de lalimentation est
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
FicUvr 5.17 Photographie en pause longue dune transition glow-arc
utilise, le dispositif exprimental permet ltablissement dune dcharge de type
glow caractrise par une tension leve (370V sur la partie (a) de la hgure 5.18)
et un courant faible (10ma). Visuellement laide dun appareil photo classique
et une pause longue, nous voyons une dcharge dinuse (hgure 5.16) qui remplit
lespace interlectrode. Limage (a) prise la camra rapide ne permet pas de
voir la partie dinuse cause du temps de pause court (16s), mais nous montre
une autre caractristique de la dcharge glow, son talement la cathode.
Au temps 1500s, la seconde partie de lalimentation est connecte via le
switch et un chelon de courant est appliqu aux bornes de la dcharge. Lam-
plitude de ce courant, de lordre de 400mA sur une dure de 2ms, est sumsante
pour provoquer une transition larc de la dcharge. Larc est caractris lec-
triquement par une chute de la tension (ici de 370V 49V sur la partie (b) de
la hgure 5.18), et une lvation du courant de dcharge (passage de 10mA
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
380mA). Le gaz devient largement conducteur et la puissance consomme dans
la dcharge passe de 3,7W 18,6W. Visuellement larc est une constriction de
la dcharge dinuse en un hlament visible sur la photographie de la hgure 5.17.
La pause longue utilise ici avec un appareil photographique classique, nous
montre la superposition de la dcharge glow et de larc. Sur limage (b) de la
hgure 5.18 prise laide de la camra rapide, le hlament est bien visible ainsi que
le rtrcissement du pied cathodique en un spot trs lumineux.
Aprs les 2ms de pulse la dconnexion de la seconde partie de lalimentation
permet la transition inverse de larc vers une dcharge luminescente qui induit
une hausse de la tension aux bornes des lectrodes et une diminution du courant
de dcharge.
Nous avons caractris ce phnomne dj largement observ de transition
irrversible par opposition aux phnomnes que nous allons dcrire dans le pro-
chain paragraphe.
!!

impulsion de courant denviron 400 mA et une distance inter-lectrode de 25 mm. Sur la
partie de gauche de la figure 3.14 lvolution au cours du temps de la tension entre les
lectrodes ainsi que lvolution de lintensit dans la dcharge sont prsentes. Sur la partie de
droite se trouvent trois photographies extraites dun film (camra rapide) et correspondant aux
instants nots (a), (b) et (c) qui sont aussi indiqus approximativement sur la figure de gauche.
Cathode
Anode
t =1250s

Cathode
Anode
t =2500s

Cathode
Anode
t =3750s

(a) 10 mA/370 V (b) 380 mA/49 V (c) 9 mA/285V
Figure 3.14 : Enregistrement typique de la tension et du courant de dcharge en fonction du temps pour une
impulsion de 2 ms et denviron 400 mA damplitude et pour une distance inter-lectrode de 25 mm.
P = 100 mbar. Photographie de la dcharge (16,6 s de temps dexposition) en rgime de glow (a+c) et en
rgime darc (c).

La tension passe de 370 V 49 V pendant que le courant crot de 10 mA 380 mA en
approximativement 100 s. La puissance moyenne de la dcharge passe de 3,7 W 18,6 W.
Le systme dimagerie rapide montre pour le rgime de glow (figure 3.14 (a) et (c)) un pied
cathodique diffus et une colonne positive diffuse. La colonne diffuse napparat pas sur les
images prsentes du fait de sa faible luminosit et du choix du temps dexposition impos
par la luminosit de larc lectrique. En revanche, le caractre diffus de la colonne est
parfaitement observable lil nu. Pour le rgime darc (figure 3.14(b)) un spot cathodique et
une colonne filamentaire sont systmatiquement observs. Pour les deux rgimes, un spot
anodique a t observ. Ce type de transition a dj t observ prcdemment par diffrents
auteurs
20
pression atmosphrique. Il faut galement noter que la transition dun glow vers un
arc lectrique peut se produire sans que le pied cathodique noccupe la totalit de la surface
cathodique.

20
S Watanabe, S Saito, K Takahashi and T Onzawa -Welding International, 2003, 17, 8, pp. 593597

FicUvr 5.18 Oscillogramme du courant et de la tension et photographies des
phases de la transition.
Transitions spontanes
Les transitions irrversibles prsentes dans la section prcdente ne sont
possibles qu partir dun chelon de courant, li la connexion de la seconde
partie de lalimentation, relativement lev (dpendant comme nous le verrons
par la suite des conditions gomtriques et de pression). Pour des amplitudes de
courant plus faibles (par exemple 280mA pour largon 100mbar pour un es-
pace interletrode de 25mm), il se produit une suite de transitions glow-arc et
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
arc-glow tout au long de lapplication (2ms) de lchelon de courant. Ces transi-
tions sont appeles spontanes car elles ne se produisent pas forcment au dbut
de lchelon, mais comme le montre loscillogramme de la hgure 5.19 peuvent
apparatre en grand nombre pendant lapplication de lchelon de courant.
!"

associes une modification structurelle du pied cathodique. Comme nous lavons prcis
dans le cas des transitions irrversibles, il nest pas non plus ncessaire dans le cas de
transitions spontanes que la totalit de la surface dlectrode soit occupe par la zone
cathodique pour quune transition survienne.

Figure 3.15 : volution au cours du temps de la tension de dcharge et de lintensit de dcharge lors de
lapplication dun pulse dintensit environ gale 200 mA (P = 100 mbar et d = 25 mm)

Pour des lectrodes de mme nature ayant cependant une gomtrie en pointe,
Watanabe et al.
16 20
nont pas observ pression atmosphrique dans largon lexistence de
transitions spontanes se produisant intensit constante. Sur la figure 3.17 nous avons
rappel les rsultats prsents par Watanabe et al. concernant la caractristique U/I pour une
distance inter-lectrode de 1 mm : aucune zone relative aux transitions spontanes (deux
valeurs de tension possibles pour une mme valeur de lintensit) napparat.
FicUvr 5.19 Oscillogramme du courant et de la tension pour un chelon de
280mA, apparition de transitions spontanes
La hgure 5.20 est un zoom de cet oscillogramme autour des temps 1700s et
2000s avec les photographies la camra rapide (temps de pause de 16,6s) des
phases associes. La partie (a) correspond au temps avant lchelon de courant,
la dcharge est en phase glow (faible courant et forte tension). Lapplication de
lchelon de courant modihe la morphologie de la dcharge luminescente, mais
sans transition larc, laspect visuel, photographie (b), est plus hlamentaire,
mais le pied cathodique reste tendu (lectrode du haut sur la photographie).
Au niveau de loscillogramme, on observe une lgre baisse de la tension et une
forte hausse du courant (environ 200mA dans ce cas). La transition complte
larc se ralise (partie (c) de la hgure 5.20) lorsquil y a une forte chute de
la tension et une constriction importante du pied cathodique en un spot lumi-
neux. Il faut noter que lors de cette transition, laspect gnral de la colonne
positive de dcharge reste inchang. Un autre fait important lors de ces transi-
tions est la variation de la puissance consomme dans la dcharge, de 56W en
mode luminescent moins de 15W pour larc. Enhn prcisons que ces transitions
spontanes ne sont pas anectes par le circuit extrieur dalimentation qui reste
inchang pendant les 2ms de temps dapplication de lchelon et sont caract-
riss par une chute de la tension (200V) qui semble associe aux changements
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
FicUvr 5.20 Zoom de loscillogramme du courant et de la tension pour un
chelon de 280mA, apparition de transitions spontanes et photographies ca-
mra rapide de la dcharge
morphologiques du pied cathodique.
Pour mieux comprendre les mcanismes entrant en jeu dans ces transitions,
nous avons focalis la suite de ce travail sur deux points, une tude statique
des caractristiques de chaque mode de dcharge (arc, glow, glow-arc) en fonc-
tion des paramtres gomtriques, de pression et damplitude de lchelon de
courant. La seconde partie traitera des mcanismes dynamiques de la transition
elle-mme.
tude paramtrique des conditions stables
Dans un premier temps, nous nous sommes intresss au rgime de glow
uniquement. En faisant varier le courant de dcharge, nous avons observ les
changements morphologiques de la dcharge grce aux images camra rapide,
mais aussi la valeur de la tension aux bornes des lectrodes. La distance inter-
lectrode est garde constante 25mm et la pression est de 100mbar. Ces obser-
vations sont reportes sur la hgure 5.21 avec la caractristique courant-tension,
mise en relation avec les photographies de la dcharge. Notons tout dabord qu
cause des temps de pause et de louverture de lobjectif, la dcharge luminescente
napparat pas sur le clich, mais est bien visible lil nu. Notons aussi la pr-
sence de stries dj rencontres dans la littrature pour ces situations. Ces points
ne reprsentent pas une exprience en continu, mais une srie dexpriences ra-
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
lises avec des chelons de courant dinrents. Si on observe maintenant le pied
cathodique (en haut sur limage), on saperoit que malgr laspect hlamentaire
prononc et ce mme pour les plus fortes valeurs de courant, celui-ci reste large-
ment tal sur lensemble de llectrode et nest pas modih par laugmentation
du courant. La tension diminue lgrement avec laugmentation du courant, sui-
vant en ce sens la taille du hlament qui tend rejoindre les deux lectrodes pour
un courant de 50mA. Le point important ici est la coexistence des deux modes
de dcharges glow et hlament sans passage proprement dit larc.
FicUvr 5.21 Caractristique courant-tension associe des images camra ra-
pide dune dcharge glow
Par la suite en gardant la mme distance interlectrode (25mm) et la mme
pression (100mbar), nous avons augment lamplitude de lchelon de courant
entre 5mA et 400mA. Pour chaque exprience, nous avons relev des tensions
stables dont les valeurs moyennes sur la priode de stabilit sont reportes
sur la hgure 5.22. Ainsi pour des courants entre 5mA et 50mA nous retrouvons
uniquement des dcharges de type glow luminescente. Une seconde zone (II)
est prsente entre les valeurs de courant 50mA et environ 280mA o des tran-
sitions spontanes et rversibles sont observes pour chaque exprience. Enhn
une troisime zone (III) a t dtermine pour des valeurs de courant suprieures
280mA pour lesquelles seules des transitions irrversibles vers larc sont pr-
sentes.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
FicUvr 5.22 Caractristique courant-tension pour une distance intrlectrode
de 25mm
tude dynamique de la transition
Nous nous sommes ensuite intresss ltude dynamique de la transition
en elle-mme. Tout dabord nous avons cherch caractriser lectriquement
et morphologiquement les changements ayant lieu partir de lapplication de
lchelon de courant par le circuit extrieur. Les mesures lectriques pour une
distance interlectrode de 25mm, une pression de 100mbar et un chelon de
25mA sont reportes sur loscillogramme de la hgure 5.23, avec en correspon-
dance des intervalles nots sur le haut de cette hgure, les images issues de la
camra rapide (temps de pause pour chaque image de 16,6s) sur la hgure 5.24.
Dans un premier temps, durant environ 25s, on remarque une hausse de la ten-
sion ainsi que ltablissement du courant la valeur hxe par lalimentation. Par
la suite la tension et le courant se stabilisent sur un plateau. Visuellement, le
hlament semble se propager depuis lanode vers la cathode pour former un tat
stationnaire avec jonction ou non des deux lectrodes en fonction de lamplitude
de lchelon de courant. Un temps de quelques dizaines de s est ncessaire pour
ltablissement dun hlament au sein de la colonne positive.
Nous avons ensuite enectu une tude paramtrique sur la distance inter-
lectrode et relev les temps de transition glow-arc pour trois chelons de cou-
rant (100mA, 150mA et 250mA). Comme nous pouvons le voir sur les rsultats
agrgs sur la hgure 5.25, les temps de transitions augmentent linairement en
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
FicUvr 5.23 Courant et tension de dcharge lors de la transition (extension de
la colonne hlamentaire
FicUvr 5.24 Images de la camra rapide en rapport avec les intervalles de la
hgure 5.23 lors dune transition
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
fonction de la distance, et ce quel que soit lamplitude de lchelon de courant.
En extrapolant ces courbes vers une distance hctive interlectrode nulle, les trois
courbes tendent vers une valeur commune proche de 50ns. Ce temps est peut
tre caractristique des phnomnes se produisant aux lectrodes (typiquement
un changement dmission lectronique la cathode). La dpendance linaire
entre le temps de transition et la distance interlectrode dune part et le courant
de dcharge dautre part semble suggrer un mcanisme propagatif dpendant
de la distance et du degr dionisation de la colonne. Ce phnomne propagatif a
une vitesse constante seulement fonction du courant de dcharge. Finalement un
temps de propagation de lordre de quelques centaines de nanosecondes semble
tre caractristique dans les conditions de nos expriences.
FicUvr 5.25 Dure de la transition en fonction de la distance interlectrode
Finalement se pose la question de savoir quel est lordre temporel des mca-
nismes de transition glow-arc. Est-ce que les phnomnes la cathode (constric-
tion du pied cathodique) prennent place avant les modihcations de la colonne
positive ou est-ce que le hlament doit tre totalement propag pour que la tran-
sition au pied cathodique ait lieu ? Pour rpondre partiellement cette question,
un exemple signihcatif illustr sur la hgure 5.26 a t trouv. Lamplitude de
lchelon de courant est de 80mA pour une distance interlectrode de 25mm et
une pression de 100mbar. Dans ces conditions, on peut observer quun temps
plus long est ncessaire pour tablir un tat similaire ceux des hgures 5.19
et 5.20. Les transitions spontanes lies une baisse signihcative de la tension
napparaissent qu partir de 2000s. Les images prises la camra rapide sur la
hgure 5.26 montrent que dinrentes structures existent avant ce temps : la d-
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
charge peut tre partiellement hlamentaire avec un pied cathodique sous forme
de spot sur les photographies (a) et (c), mais aussi partiellement hlamentaire avec
un pied cathodique dinus. Les modihcations du pied cathodique sont donc pos-
sibles mme si la dcharge na pas atteint un aspect compltement hlamentaire.
En enet sur les images le pied cathodique oscille entre un spot et un mode dinus
durant le temps de propagation du hlament au sein de lespace interlectrode.
En conclusion nous pouvons avancer lhypothse que le dveloppement hla-
mentaire complet nest pas une condition ncessaire pour ltablissement dun
spot cathodique, mais aussi que celui-ci nest pas une condition sumsante pour
mener une transition complte larc. Larc lectrique ncessite donc une nou-
velle dhnition structurelle : la colonne de dcharge doit tre totalement hlamen-
taire et le pied cathodique doit tre sous la forme dun spot (constriction du pied
cathodique).
FicUvr 5.26 Tension et courant de dcharge pour une srie de transitions glow-
arc et arc-glow spontanes et rversibles
Pour tre totalement complet ce modle comportemental expos ici devra
tre tendu un modle comprenant et expliquant en dtail les mcanismes
internes qui conduisent ces transitions.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
FicUvr 5.27 Amazing Stories avril 1926
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Chapitre 6
Conclusion gnrale et perspectives
scientifques
6.1 Conclusion
Ce mmoire en vue de lobtention de lhabilitation diriger les recherches
fait lobjet dun dcoupage en 3 parties. La premire sur la modlisation et lop-
timisation des grands rseaux lectriques, la deuxime sur loptimisation sys-
tmes des machines lectriques et la dernire sur la modlisation des phno-
mnes physiques. Ces trois approches de llectrotechnique, premire vue dif-
frentes, sont associes par la volont de raliser des optimisations et modlisa-
tions globales des systmes. Pour cela, il faut dvelopper des algorithmes perfor-
mants, tant du point de vue des rsultats que de la rapidit de calcul. Cela passe
par une phase de dveloppement sur des cas simples et des exemples concrets
dj rsolus. Pour bien poser les problmes doptimisation, il faut de plus avoir
une connaissance approfondie des phnomnes physiques. Cest pourquoi nous
avons notamment continu travailler sur les mcanismes de claquage dans les
gaz isolants, ceci en vu de bien comprendre les contraintes poss lors des opti-
misations.
6.2 Perspectives scientifques
Mes travaux de recherche, tels que dj prsents, ont ports sur le dvelop-
pement de modles et dalgorithmes doptimisation appliqus plusieurs do-
maines dtude. court et moyen termes, mes perspectives de recherche se
structurent autour de deux principaux axes de travail : loptimisation dans le
domaine de lnergie, en particulier des systmes dentranement, et la modli-
sation par simulations numriques applique la physique des plasmas.
105
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
106 Cu~vi1vr 6 : Conclusion gnrale et perspectives scientifques
Lobjet du prsent document est de prsenter lintrt port pour ces thma-
tiques, de dvelopper les orientations de recherche envisages dans le cadre de
thses en cours et programmes, et de mentionner les collaborations mises en
place, ainsi que les moyens qui seront dploys.
6.2.1 Optimisation des systmes dnergie
Le durcissement progressif du contexte rglementaire relatif aux missions
de gaz enet de serre et les attentes croissantes de la socit pour le dvelop-
pement dinnovations propres, conduisent les industriels sorienter vers luti-
lisation de moteurs lectriques. Au-del de la motivation environnementale, ces
technologies sont galement utilises en raison de leurs performances cono-
miques. Les moteurs lectriques investissent les dinrents compartiments des in-
dustries aronautique, spatiale, navale et automobile, et les besoins de recherche,
notamment doptimisation, saccroissent fortement.
Ainsi, les avionneurs tendent remplacer la plupart des actionneurs hydrau-
liques par des chanes de motorisations lectriques, la conception navale cherche
remplacer les modes de propulsion thermiques par des groupes hybrides (puis-
sance thermique et propulsion lectrique). Ces deux derniers concepts sont sou-
vent reconnus dans les termes avions ou bateaux tout (plus) lectriques. Les
constructeurs automobiles et leurs fournisseurs ont dans un premier temps cher-
ch remplacer ou crer des fonctions o llectrihcation des organes tient une
place centrale. Depuis quelques annes compte tenu de la ncessit de rduire la
consommation et les rejets de polluants atmosphriques, les moteurs lectriques
prennent leur place dans la propulsion du vhicule. La question de la masse ou
du volume est critique dans laronautique, alors que le cot de fabrication est
un lment clef dans lautomobile.
Compte tenu de ce contexte, divers axes de recherches sont proposs. La
prise en compte de nouveaux modles physiques comportementaux plus per-
formants et plus prcis va ncessiter une amlioration des algorithmes dopti-
misation qui devront tre eux aussi plus prcis et surtout plus rapides. En enet,
la prise en compte de la mcanique (transformateurs et rducteurs) et de llec-
tronique de puissance est encore ralise de manire trop simple. Pour la partie
mcanique, il va nous falloir faire des rapprochements avec des spcialistes de
la conception mcanique pour raliser des modles plus complexes des trans-
formateurs de mouvement. La modlisation de llectronique de puissance en
vue de son incorporation dans loptimisation devrait elle se faire rapidement en
collaboration avec Pierre Lefranc du dpartement nergie de Suplec.
La complexihcation des modles va ncessiter de travailler sur les algo-
rithmes doptimisation ahn de les rendre plus performants. De ce point de vue, la
premire amlioration apporter est de bnhcier de laugmentation du nombre
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
6.2 Perspectives scientifques 107
de curs de calcul au sein des ordinateurs modernes ainsi que du calcul sur
architectures parallles. Un tel cluster de machines est dj en place Suplec,
nous projetons donc de parallliser les codes dj existants sur cette architec-
ture. Les algorithmes stochastiques du type algorithme gntique forment un
trs bon exemple au vu de leur caractre massivement parallle. Des premiers
essais dacclration de calcul seront raliss trs prochainement lors de la thse
de Benjamin Dagus.
Une seconde piste de recherche est de sorienter vers des optimisations ou la
fonction objectif nest plus sous forme analytique (ou simple), mais reprsente
comme le rsultat de calculs par lments hnis (calcul du couple dune machine
par exemple). De telles optimisations ne peuvent tre menes laide des al-
gorithmes classiques, il est donc ncessaire dinventer de nouvelles approches.
Une piste possible est lutilisation de mthodes o le nombre dvaluations de
la fonction cot est minimis pour rendre loptimisation dans un temps de cal-
cul raisonnable. Ces mthodes font appel des techniques dapproximation de
la fonction cot, cette approximation tant ramne au fur et mesure des it-
rations uniquement aux endroits de loptimum probable. Ces dveloppements
feront lobjet dune collaboration avec le dpartement Signaux et Systmes lec-
troniques et la chaire de Systmes analogiques complexes Suplec.
Ces nouvelles mthodes doptimisation, aprs validation sur les systmes
dentranement lectriques, seront appliques llectronique de puissance, en
particulier la problmatique du placement routage de composants sur une carte
lectronique. En enet, le placement de composants a un enet direct sur les cou-
plages inductifs et sur lchaunement du dispositif. Il est donc ncessaire par le
biais de simulations et doptimisations de trouver le placement des composants
minimisant les fuites et lchaunement. Lanalyse thermique et lectromagn-
tique ne peut tre faite de manire prcise que par lments-hnis, ce qui corres-
pond bien au thme de recherche prsent. La problmatique ici sera renforce
par le grand nombre de degrs de liberts prsents sur le systme.
6.2.2 Modlisation des processus plasmas
La physique des plasmas et plus prcisment des dcharges lectriques est
relativement bien connue, mais sa modlisation hne par le biais de simulations
numriques est encore une question non totalement rsolue. En enet, le couplage
de la dynamique des dcharges avec le circuit lectrique, et de la dynamique des
gaz avec la cintique chimique des ractions restent encore des enjeux scienti-
hques majeurs. Cest pourquoi dans le cadre dun dpt de projet ANR, nous
nous proposons dtudier ces couplages au sein dune collaboration entre plu-
sieurs laboratoires. Ce projet vise la ralisation dun logiciel de simulation
comportant une phase de simulation du passage de la dcharge dans un gaz qui
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
se couplerait des codes de simulation dexpansion hydrodynamique et de ci-
ntique chimique. Toutes ces modlisations seront alimentes par des donnes
exprimentales ralises pour ce projet. Le contexte de cette tude sinscrit dans
un cadre environnemental de la rduction des composs polluants dans des gaz.
Ma participation ce projet portera sur la modlisation de la phase de dcharge.
La multiplication de la puissance des calculateurs rend possible la simulation en
trois dimensions de la dcharge, ainsi que lutilisation de codes de calculs par
lments hnis ou par des mthodes mixtes du type ELLAM (Euler Lagrange Lo-
cal Adjoint Methods). Ces algorithmes couplent des mthodes de rsolution par
lments hnis avec des mthodes des caractristiques, permettant ainsi dutili-
ser des gomtries complexes tout en gardant les avantages des mthodes des
caractristiques pour la modlisation des quations fuides de conservation. En
complment, il sera ncessaire de parallliser ces codes de calcul.
Une deuxime collaboration vient de dbuter avec le laboratoire EM2C de
lcole Centrale de Paris par le biais du co-encadrement dune thse portant sur
ltude de lemcacit des plasmas froids pour la production dhydrogne. Lutili-
sation de lhydrogne comme vecteur dnergie prend une importance croissante
en raison de la ncessit du stockage dnergie pour lapplication aux trans-
ports et lutilisation de sources dnergie renouvelables. La production emcace
dhydrogne devient donc un enjeu majeur. Plusieurs techniques et procds
de production existent aujourdhui, mais la recherche de technologies alterna-
tives plus emcaces dun point de vue nergtique demeure trs active. Ainsi les
plasmas froids pression atmosphrique pourraient constituer une technolo-
gie prometteuse. En parallle avec une approche exprimentale qui permettra
de mieux dcrire les phnomnes mis en jeux dans les racteurs plasma, un
modle thermique et physique dcrivant les mcanismes dinteractions entre la
dcharge et les hydrocarbures lgers sera dvelopp et valid face ces exp-
riences. Cette double approche de modlisation et dexprimentation permettra
ensuite doptimiser lemcacit nergtique de la production dhydrogne.
Enhn, une thse sur ltude du courant noir par une surface soumise un
champ lectrique intense et des mcanismes de transition larc vient de d-
buter, en partenariat avec Philippe Test du Laboratoire de Gnie lectrique de
Paris et Emmanuel odic du dpartement nergie de Suplec. Dans divers dispo-
sitifs mettant en uvre des champs lectriques levs et des pressions rduites
tels que les acclrateurs de particules (acclrateur de neutres pour ITER par
exemple), il apparat parfois des courants non prvus. Ces courants sont appels
courants noirs et sont susceptibles de fortement dgrader les performances du
dispositif. Divers phnomnes peuvent contribuer lmission de ces courants.
Lobjectif de ltude exprimentale sera dune part didentiher les paramtres
dominants conduisant ltablissement du courant noir dont on proposera un
modle physique, et dautre part dlaborer des solutions permettant de rduire
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
6.2 Perspectives scientifques 109
ce phnomne. Dans un second temps les conditions de passage larc seront
tudies par le biais dune modlisation qui permettra de proposer un scnario
expliquant la transition brutale de lmission lectronique larc lectrique.
6.2.3 Conclusions
Ces perspectives de recherche sinscrivent dans le cadre vaste de loptimi-
sation et la modlisation des systmes dnergie. Deux domaines sont prsen-
ts, le premier traite de loptimisation des systmes de conversion de lnergie
(entranement lectrique et lectronique de puissance) et le second de la phy-
sique des plasmas. La complexihcation de ces modles passe bien sr par une
phase danalyse physique de comprhension des phnomnes. Une caractris-
tique commune ces deux axes de recherche est la ncessaire augmentation
des cots de calculs lis la complexihcation des modles en cours de dvelop-
pement. Pour remdier cette dimcult, diverses pistes sont envisages : algo-
rithmes innovants et paralllisation des codes de calcul.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
FicUvr 6.1 Amazing Stories dcembre 1926
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Troisime partie
Bibliographie
111
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Chapitre 7
Production scientifque
7.1 Tableau rcapitulatif
Revues avec comit de lecture : 8 ;
actes de confrences dites et confrences invites : 6 ;
congrs avec actes : 44 ;
communications : 2 ;
articles accepts en cours de publication : 1 ;
articles soumis en cours de relecture : 2.
7.2 Actes de confrences dites, confrences invi-
tes et chapitres de livres
[1] X~virv J~NNo1, Puiiivvr Drss~N1r, Pirvvr Viu~i r1 Jr~NCi~Uur V~N
Nirv, Optimization of wound rotor synchronous generators based on genetic
algorithm, in Computational Methods for the Innovative Design of Elec-
trical Devices., vol. 327 de Studies in Computational Intelligence, Springer
Verlag, 2010, p. 87113.
[2] E. M~vour, D. DJrvxoUNr, P. Drss~N1r, C. DrNisr1, P. ScUv, F. B~s
1irN, A. BoUvuoN r1 C. L~Ux, Physics and applications of atmospheric non-
thermal air plasma with reference to environment, in 36th EPS Conference
on Plasma Physics, Soha, Bulgaria, 2009, 17 pages.
[3] A. Goiux~N, M. Goiux~N, E. Ouic r1 P. Drss~N1r, Partial discharges in-
ception and ageing efects in a gas insulated HV equipment, in Proc. Gaseous
dielectrics 2008, Cardin Royaume-Uni, 2008, p. 49 56.
[4] N. Briirc~vur, P. Drss~N1r, P. Viu~i r1 J.C. V~NNirv, Optimisation of a
drive system and its epicyclic gear set, in Intelligent Computer Techniques
113
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
114 BIBLIOGRAPHIE
in Applied Electromagnetics, vol. 119, Springer Berlin Heidelberg, Berlin,
Heidelberg, 2008, p. 241248.
[5] P. Drss~N1r, J. V~NNirv r1 C. Rivoii, Optimization of a linear brushless DC
motor drive, in Recent Developments of Electrical Drives, Springer Nether-
lands, Dordrecht, 2006, p. 127136.
[6] P. Drss~N1r, J. V~NNirv r1 B. BoN~ios, Active shielding optimisation, in
Computer Engineering in Applied Electromagnetism, Springer-Verlag, Ber-
lin/Heidelberg, 2005, p. 3540.
7.3 Revues internationales avec comit de lecture
[7] P. Lriv~Nc, X. J~NNo1 r1 P. Drss~N1r, Virtual prototyping and pre-sizing
methodology for buck DC-DC converters using genetic algorithms, Power
Electronics, IET, 5 (2012), p. 41 52.
[8] G. S~NuoU, H. BovsrNnrvcrv r1 P. Drss~N1r, Unit commitment with pro-
duction cost uncertainty : A recourse programming method, Journal of Energy
and Power Engineering, 5 (2011), p. 164172.
[9] V. RioUs, J.M. Gi~cu~N1, Y. Prvrz r1 P. Drss~N1r, Linsufsance des si-
gnaux de localisation pour la coordination entre la production et le transport
dlectricit dans les systmes lectriques libraliss, Revue conomique, 60
(2009), p. 819829.
[10] E. M~vour, D. DJrvxoUNr, P. Drss~N1r, C. DrNisr1, P. ScUv, F. B~s
1irN, A. BoUvuoN r1 C. L~Ux, Physics and applications of atmospheric
non-thermal air plasma with reference to environment, Plasma Physics and
Controlled Fusion (IF 2,4), 51 (2009), p. 124002, 15 pages.
[11] V. RioUs, J.M. Gi~cu~N1, Y. Prvrz r1 P. Drss~N1r, The diversity of design
of TSOs, Energy Policy (IF 2,4), 36 (2008), p. 33233332.
[12] J. MrNuoz~, R. Lovrz, D. Mov~irs, E. Lovrz, P. Drss~N1r r1 R. Mov~c~,
Minimal loss reconfguration using genetic algorithms with restricted popula-
tion and addressed operators : Real application, IEEE Transactions on power
systems (1,94), 21 (2006), p. 948954.
[13] E. Ouic, M. Du~iN~U1, A. Goiux~N, M. Goiux~N r1 P. Drss~N1r, Study in
space and time of the gas temperature variations in dielectric barrier discharge
reactors, Journal of Advanced Oxydation Technologies (IF 0,79), 8 (2005),
p. 133141.
[14] S. Dri~uvrUr, P. Drss~N1r, F. BvoU~vr r1 P. B~s1~vu, Using cluster ana-
lysis in power system planning under uncertainty, Engineering Intelligent
Systems for Electrical Engineering and Communications (IF 0,2), 11 (2003),
p. 18792.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
BIBLIOGRAPHIE 115
7.3.1 Confrences internationales avec comit de lecture
[15] P. MoiiNi, R. H~NN~, T. P~Uixirv, B. Div~ssrN, P. Drss~N1r, M. Bri
u~J, D. P~v~N r1 N. B~icoN, Polyimide and FEP charging behavior under
multienergetic Electron-Beam irradiation : Experiments and a simple model,
in 14th International Symposium on Electrets ISE14, Montpellier, France,
2011, 7 pages.
[16] R. L~Nuiviru, R. ANui~Urv, P. Drss~N1r, M. Kivcxv~1vicx, T. Lrni~Nc,
E. Ouic r1 P. Trs1, Investigation of a glow discharge structure in 100 mbar
argon, in 30th International Conference on Phenomena in Ionized Gases
ICPIG, Belfast, Northern Ireland, UK, sept. 2011, 3 pages.
[17] C. GU1iviNu, J. V~NNirv, P. Viu~i r1 P. Drss~N1r, Analytical study of an
optimized limited motion actuator used in engine combustive fow regula-
tion, in XV International Symposium on Electromagnetic Fields in Mecha-
tronics, Electrical and Electronic Engineering ISEF2011, Funchal, Portugal,
sept. 2011, 9 pages.
[18] B. D~cUs, P. Drss~N1r, P. Viu~i, J. V~NNirv, J. S~iN1Micuri r1 J. Tuo
x~s, Optimization and comparison of optimal saliency permanent magnet
synchronous machines for electric vehicle application, in XV International
Symposium on Electromagnetic Fields in Mechatronics, Electrical and Elec-
tronic Engineering ISEF2011, Funchal, Portugal, sept. 2011, 8 pages.
[19] V. RioUs, J. Gi~cu~N1 r1 P. Drss~N1r, Transmission network investment as
an anticipation problem, in DIME Workshop "The changing Governance of
Network Industries", Naples, Italie, avril 2010, 32 pages.
[20] M. Pr1i1 r1 P. Drss~N1r, Modeling and experimental set-up of an automo-
tive electrical network for transient studies, in International Conference on
Electrical Machines, ICEM2010, Rome, Italy, sept. 2010, 6 pages.
[21] H. BovsrNnrvcrv, P. Drss~N1r r1 G. S~NuoU, Unit commitment with pro-
duction cost uncertainty, a recourse programming method, in IFAC Confe-
rence on Control Methodologies and Technology for Energy Emciency, Vi-
lamoura, Portugal, mars 2010, 5 pages.
[22] V. RioUs, J. Gi~cu~N1, Y. Prvrz r1 P. Drss~N1r, The role of transmission
investment in the coordination between generation and transmission in the li-
beralized power systems, in 32nd IAEE International Conference, San Fran-
cisco, California, USA, juin 2009, 9 pages.
[23] V. RioUs r1 P. Drss~N1r, Real gains from fow-based methods for allocating
power transmission capacity in europe, in 6th International Conference on
the European Energy Market, Leuven Belgique, avril 2009, 6 pages.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
116 BIBLIOGRAPHIE
investment in the coordination between generation and transmission in the
liberalized power systems, in 13th Annual Conference of The International
Society for New Institutional Economic, Berkeley, USA, juin 2009, 9 pages.
[25] M. Pr1i1, P. Drss~N1r r1 E. L~Uv~iN, Modlisation et simulation en rgime
dynamique dun rseau de bord automobile, in lectrotechnique du Futur,
EF2009, Compigne France, sept. 2009, 8 pages.
[26] S. MR~u, P. Lriv~Nc, P. Drss~N1r, P. Cuiozzi, G. BioNuri, M. F~xrs r1
P. M~ssoN, A compact transient electrothermal model for integrated power
systems : Automotive application, in IECON09 Annual Conference of the
IEEE Inductrial Electronics Society, Porto : Portugal, nov. 2009, p. 3755
3760.
[27] X. J~NNo1, K. G~Uurivov, N. ViNri, P. Viu~i r1 P. Drss~N1r, Optimal
design of a wound rotor synchronous generator using genetic algorithm, in
ISEF 2009 - XIV International Symposium on Electromagnetic Fields in Me-
chatronics, Electrical and Electronic Engineering, Arras, France, sept. 2009,
6 pages.
[28] V. RioUs, J. Us~oi~, M. S~cU~N, J. Gi~cu~N1 r1 P. Drss~N1r, Assessing
available transfer capacity on a realistic european network : Impact of as-
sumptions on wind power generation, in Proceedings of the Conference :
"Building Networks for a Brighter Future" 1st International scientihc confe-
rence, Rotterdam Pays-Bas, 2008, 6 pages.
[29] V. RioUs, Y. Prvrz r1 P. Drss~N1r, The efciency of short run and long
run locational signals to coordinate generation location with lumpy transmis-
sion investments, in Proceedings of the 31st IAEE International Conference.
Bridging Energy Supply and Demand : Logistics, Competition and Envi-
ronment, Istanbul Turquie, juin 2008, 21 pages.
investment in the coordination between generation and transmission in the
liberalized power systems, in 7th Conference on Applied Infrastructure Re-
search, Berlin Germany, oct. 2008, 9 pages.
[31] V. RioUs, J. M. Gi~cu~N1, Y. Prvrz r1 P. Drss~N1r, Le rle de linvestisse-
ment rseaux dans la coordination entre la production et le transport dlec-
tricit dans les systmes lectriques libraliss, in 57me Congrs Annuel
de lAssociation Franaise de Sciences Economiques (AFSE), Paris, France,
sept. 2008, 10 pages.
[32] V. RioUs, J. Gi~cu~N1, Y. Prvrz r1 P. Drss~N1r, Anticipation for efcient
electricity transmission network investments, in Infrastructure Systems and
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
BIBLIOGRAPHIE 117
Services : Building Networks for a Brighter Future (INFRA), 2008 First In-
ternational Conference on, Rotterdam Pays-Bas, 2008, 6 pages.
[33] V. RioUs, Y. PuUiviN, F. Pr1i1, P. Drss~N1r r1 M. S~cU~N, A contribution in
experimental economics to classify market power behaviors in an oligopolistic
power market, in 5th International Conference on the European Electricity
Market, Lisbon, PORTUGAL, mai 2008, IEEE, p. 271276.
[34] V. RioUs, Y. Prvrz r1 P. Drss~N1r, Is combination of nodal pricing and ave-
rage participation tarif the best solution to coordinate the location of power
plants with lumpy transmission investments ?, in Proceedings of The Econo-
mics of Energy Market, Toulouse France, juin 2008, 23 pages.
[35] W. FvriiN, M. Pr1i1, J. V~NNirv, L. Brv1ur1 r1 P. Drss~N1r, Evaluation of
thermal heating in LV-cable in presence of harmonics, using FEM methods,
in International Conference on Electrical Machines, ICEM2008, Villamoura
Portugal, 2008, 6 pages.
[36] H. B~vrx, M. Pr1i1, P. Drss~N1r, F. DUioUvu r1 R. Briuoxxr, Wind inte-
gration study for an island network using a probabilistic tool, in Proceedings
of European Wind Energy European Wind Energy Conference & Exhibi-
tion, EWEC, Bruxelles Belgique, 2008, 8 pages.
[37] H. B~vrx, Y. PuUiviN, P. Drss~N1r r1 J. Brc1, Probabilistic computation of
a wind farm generation based on wind turbine dynamic modeling, in 10th
International Conference on Probabilistic Methods Applied to Power Sys-
tems, PMAPS, Rincon, PR, mai 2008, IEEE, p. 517522.
[38] V. RioUs r1 P. Drss~N1r, A limited place for competition to develop the power
transmission network, in 9th IAEE European Energy Conference : Energy
Markets and Sustainability in a Larger Europe, Florence, Italie, juin 2007,
19 pages.
[39] F. Pr1i1, Y. PuUiviN, P. Drss~N1r r1 J. Gi~cu~N1, A contribution of ex-
perimental economics toward characterization of the use of market power in
oligopolistic markets, in 9th IAEE European Energy Conference : Energy
Markets and Sustainability in a Larger Europe, Florence, Italie, juin 2007,
12 pages.
[40] F. M~Uv~s, P. Drss~N1r r1 J. Gi~cu~N1, Defnition and controllability of
wind power integration cost, in 9th IAEE European Energy Conference :
Energy Markets and Sustainability in a Larger Europe, Florence, Italie, juin
2007, 15 pages.
[41] M. Lovrz, P. Drss~N1r, D. Mov~irs, J. V~NNirv r1 D. S~u~vN~c, Optimisa-
tion of a small non controlled wind energy conversion system for Stand-Alone
applications, in International Conference on Renewable Energies and Po-
wer Quality, ICREPQ, Seville Espagne, mars 2007, 5 pages.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
118 BIBLIOGRAPHIE
[42] N. Briirc~vur, J. V~NNirv, P. Viu~i r1 P. Drss~N1r, Multi objective opti-
misation of an electromagnetic valve actuator, in Sensors and Actuators for
Advanced Powertrain Application, SENSACT2007, Gif sur Yvette, France,
oct. 2007, 6 pages.
[43] N. Briirc~vur, P. Drss~N1r, P. Viu~i r1 J. V~NNirv, Optimisation of a drive
system and its epicyclic gear set, in International Simposium on Electroma-
gnetique Field, ISEF2007, Prague, Rpublique Tchque, sept. 2007, 6 pages.
[44] H. B~vrx, M. Pr1i1, P. Drss~N1r, F. DUioUv r1 R. Briuoxxr, Probabilistic
characterisation of wind farms for grid connection studies, in European Wind
Energy Conference & Exhibition, EWEC, Milan, Italie, mai 2007, 8 pages.
[45] M. S~cU~N, J. Gi~cu~N1 r1 P. Drss~N1r, Assessing economic consequences
of the european electricity market design choices : The case of balancing me-
chanism, in 29th IAEE (International Association for Energy Economics)
International Conference, Berlin Germany, 2006, 20 pages.
[46] M. S~cU~N, N. Krsrvic, P. Drss~N1r r1 J. Gi~cu~N1, Market power in po-
wer markets : Game theory vs. agent-based approach, in IEEE/PES Transmis-
sion and Distribution Conference and Exposition Latin America, Caracas,
VENEZUELA, aot 2006, p. 10811086.
[47] M. S~cU~N, J. Gi~cu~N1 r1 P. Drss~N1r, Risk management and optimal
hedging in electricity forward markets coupled with a balancing mechanism,
in 2006 International Conference on Probabilistic Methods Applied to Po-
wer Systems, Stockholm, SWEDEN, juin 2006, IEEE, p. 917922.
[48] M. BrNsr11i, P. Drss~N1r, P. Viu~i r1 J. V~NNirv, Study and optimisation of
an electrical generator with original structure, in XVII International Confe-
rence on Electrical Machines, ICEM, Chania, Crete, Greece, 2006, 5 pages.
[49] P. Drss~N1r, J. V~NNirv r1 P. Viu~i, Optimisation of a linear brushless dc
moteur drive and the associated power supply, in Fith Symposium on the
all-electric ship, AES2005, Versailles, France, oct. 2005, 5 pages.
[50] M. S~cU~N, S. PiUxri, P. Drss~N1r, J. Gi~cu~N1 r1 P. B~s1~vu, Gene-
tic algorithm associated to game theory in congestion management, in 2004
International Conference on Probabilistic Methods Applied to Power Sys-
tems, PMAPS 2004, Ames, IA, USA, 2004, p. 41520.
[51] P. Drss~N1r, J. V~NNirv r1 C. Rivoii, Optimisation of a linear brushless DC
motor drive, in 16th International Conference on Electrical Machines, ICEM
2004, Cracow, Poland, sept. 2004, Institute of Mechatronics and Information
Systems, 5 pages.
[52] J. V~NNirv, P. Viu~i r1 P. Drss~N1r, Optimisation dun ensemble conver-
tisseur, actionneur lectrique et pompe pour un systme dalimentation en
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
7.4 Communications 119
carburant, in lectrotechnique du Futur, EF2003, Suplec, Gif-sur-Yvette,
France, 2003, 7 pages.
[53] M. Du~iN~U1, E. Ouic, P. BUiN, P. Hrvzoc, M. Goiux~N, A. Goiux~N,
P. Drss~N1r r1 C. K~vixi, From acoustic waves to the gas temperature in
a point-to-plane dielectric barrier discharge in air at atmospheric pressure, in
16th International Symposium on Plasma Chemistry, ISPC, Taormina, Italy,
2003, 6 pages.
[54] P. Drss~N1r, J. V~NNirv r1 B. BoN~ios, Active shielding optimisation, in
11th International Symposium on Electromagnetic Fields in Electrical En-
gineering, ISEF 2003, Maribor, Slovnie, sept. 2003, p. 631636.
[55] P. Drss~N1r r1 E. M~vour, Modlisation dune dcharge lectrique de type
streamer en gomtrie pointe-plan, in lectrotechnique du Futur, EF2003,
Suplec, Gif-sur-Yvette, France, 2003, 6 pages.
[56] B. BoN~ios, P. Drss~N1r, J. V~NNirv, P. LocNoNNr, J. K~vczrvsxi r1
N. S1irnic, Analysis and parameter identifcation of a PM linear actuator for
space application, in International Simposium on Applied Electromagntics
and Mechanism, ISEM 2003, Versailles, France, 2003, p. 252253.
[57] S. Dri~uvrUr, P. Drss~N1r, F. BvoU~vr r1 P. B~s1~vu, Using cluster analy-
sis in power system planning under uncertainty, in International Conference
on Probabilistic Methods Applied to Power Systems, PMAPS 2002, Naples,
Italie, sept. 2002, p. 479483.
[58] E. M~vour, P. Drss~N1r, N. Drscu~xvs r1 C. DrNisr1, Diagnostics and
modeling of high pressure streamer induced discharges, in Spring Meeting
on Diagnostics of Non-Equilibrium High Pressure Plasmas, Bad Honnef,
Germany, 2001, Univ. Essen, p. 3948.
7.4 Communications
[59] H. BovsrNnrvcrv, P. Drss~N1r r1 G. S~NuoU, Optimisation de la produc-
tion dans un rseau dnergie soumis des incertitudes sur le cot de produc-
tion, in Recherche Oprationnelle et dAide la Dcision, ROADEF 2010,
Toulouse, France, fv. 2010, 2 pages.
[60] J. V~NNirv, P. Viu~i, P. Drss~N1r r1 F. DUcU, Ensemble convertisseur ac-
tionneur lectrique et pompe pour un systme dalimentation en carburant
dun satellite, in Journes 2004 du club EEA, Cergy-Pontoise, France, mars
2004.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
120 BIBLIOGRAPHIE
FicUvr 7.1 Amazing Stories septembre 1926
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Chapitre 8
Rfrences bibliographiques
[AP98] Joaquin Acevedo and Efstratios N. Pistikopoulos. Stochastic opti-
mization based algorithms for process synthesis under uncertainty.
Computers & Chemical Engineering, 22(4-5) :647671, 1998.
[BN99] A. Ben-Tal and A. Nemirovski. Robust solutions of uncertain linear
programs. Operations Research Letters, 25(1) :113, August 1999.
[BN02] Aharon Ben-Tal and Arkadi Nemirovski. Robust optimization ? me-
thodology and applications. Mathematical Programming, 92 :453
480, May 2002.
[DR97] Darinka Dentcheva and Werner Rmisch. Optimal power generation
under uncertainty via stochastic programming. In : Stochastic Pro-
gramming Methods and Technical Applications (K. Marti And P. Kall
Eds.), Lecture Notes in Economics and Mathematical Systems, 458 :22
56, 1997.
[EF97] E.A.Amerasekera and F.N.Najm. Failure Mechanisms in Semiconduc-
tor Devices, 2nd Edition. WILEY, ISBN : 978-0-471-95482-8, 1997.
[EKL] E.Zitzler, K.Deb, and L.Thiele. Comparison of multiobjective evolu-
tionary algorithms : Empirical results. In Evolutionary Computation
8, EC2000.
[GR] Nicole Grwe-Kuska and Werner Rmisch. Stochastic unit commit-
ment in Hydro-Thermal power production planning.
[HAD05] H.N.Nagaraja, A.Patra, and D.Kastha. Design optimization of coupled
inductor multiphase synchronous buck converter. In IEEE Internatio-
nal Conference on Industrial Technology, ICIT05, Hong Kong, 14-17
December 2005. DOI : 10.1109/ICIT.2005.1600735.
[KAST02] K.Deb, A.Pratab, S.Agrawal, and T.Meyarivan. A fast elitist non-
dominated sorting genetic algorithm for multi-objective optimiza-
121
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
122 BIBLIOGRAPHIE
tion : Nsga-ii. IEEE Trans. on Evolutionary Computation, 6(2) :182197,
April 2002.
[KCPC09] K.Ejjabraoui, C.Larouci, P.Lefranc, and C.Marchand. A new pre-
sizing approach of dc-dc converters, application to a boost converter
for the automotive domain. In 35
th
Annual Conference of the IEEE
Industrial Electronics Society, IECON09, Porto, Portugal, November
2009.
[PCW00] A.B. Philpott, M. Craddock, and H. Waterer. Hydro-electric unit com-
mitment subject to uncertain demand. European Journal of Operatio-
nal Research, 125(2) :410424, September 2000.
[Sah04] Nikolaos V Sahinidis. Optimization under uncertainty : State-of-the-
art and opportunities. COMPUTERS AND CHEMICAL ENGINEERING,
28 :971983, 2004.
[TD03] T.C.Neugebauer and D.J.Perreault. Computer-aided optimization of
dc-dc converters for automotive applications. IEEE Trans. on Power
Electronics, 18(3), May 2003.
[XTF97] X.Zhou, T.G.Wang, and F.C.Lee. Optimizing design for low voltage
dc-dc converters. In IEEE Applied Power Electronics Conference and
Exposition, APEC97, Vancouver, Canada, 23-27 February 1997. DOI :
10.1109/APEC.1997.575633.
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Table des matires
I Fiche synthtique 3
1 CV scientifque 5
1.1 tat civil . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.2 Titres universitaires . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.3 Parcours professionnel . . . . . . . . . . . . . . . . . . . . . . . 6
1.4 Activits denseignement . . . . . . . . . . . . . . . . . . . . . . 6
1.5 Activits lies ladministration . . . . . . . . . . . . . . . . . . 6
1.6 Activits lies la recherche . . . . . . . . . . . . . . . . . . . . 6
1.6.1 Prix . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.6.2 Organisations de congrs scientihques . . . . . . . . . . 7
1.6.3 Comits de lecture . . . . . . . . . . . . . . . . . . . . . 7
1.6.4 Participation des jurys de thse . . . . . . . . . . . . . 7
1.6.5 Partenariats industriels et projets de ples de comptitivit 8
1.7 Encadrements . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.7.1 Thses . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.7.2 Stagiaires post-doc . . . . . . . . . . . . . . . . . . . . . 12
1.7.3 Stagiaires Master 2 . . . . . . . . . . . . . . . . . . . . . 12
1.8 Synthse scientihque . . . . . . . . . . . . . . . . . . . . . . . . 13
II Synthse des travaux de recherche 21
2 Introduction 23
3 Optimisation et modlisation technico-conomique des rseaux lec-
triques 27
3.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
3.2 Designs de marchs . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.3 Linvestissement dans les rseaux . . . . . . . . . . . . . . . . . 28
3.4.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . 29
123
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
124 TABLE DES MATIRES
3.4.2 Impact de la hlire olienne sur les marchs de llectricit 30
3.4.3 Intgration des moyens de productions incertains . . . . 33
3.4.4 Modlisation et optimisation robuste de systmes com-
plexes : application aux rseaux nergtiques . . . . . . 37
3.5 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
4 Optimisation et modlisation des machines lectriques 43
4.1 Principes de loptimisation des systmes lectrotechniques . . . 43
4.2 Optimisation des cots de fabrication et du rendement dune
gamme de machine . . . . . . . . . . . . . . . . . . . . . . . . . 45
4.2.1 Modlisation de la machine . . . . . . . . . . . . . . . . 46
4.2.2 Optimisation dune gamme dalternateurs . . . . . . . . 47
4.2.3 Optimisation multi-objectif du cot et du rendement . . 54
4.3 Prototypage virtuel en lectronique de puissance . . . . . . . . . 55
4.3.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . 55
4.3.2 Modlisation du hacheur abaisseur . . . . . . . . . . . . 56
4.3.3 Dhnition de loptimisation et prototypage virtuel du ha-
cheur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
4.3.4 Rsultats de loptimisation . . . . . . . . . . . . . . . . . 64
4.4 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
5 Modlisation des phnomnes physiques 71
5.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
5.2 tude des mcanismes dinfammation dun matriau isolant en
prsence dun point chaud dorigine lectrique . . . . . . . . . . 72
5.3 Physique des dcharges lectriques . . . . . . . . . . . . . . . . 83
5.3.1 Description et modlisation des phnomnes physiques
des dcharges lectriques dans les plasmas froids . . . . 83
5.3.2 tude des mcanismes de transitions glow-arc dans lar-
gon faible pression . . . . . . . . . . . . . . . . . . . . 91
6 Conclusion gnrale et perspectives scientifques 105
6.1 Conclusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
6.2 Perspectives scientihques . . . . . . . . . . . . . . . . . . . . . . 105
6.2.1 Optimisation des systmes dnergie . . . . . . . . . . . 106
6.2.2 Modlisation des processus plasmas . . . . . . . . . . . 107
6.2.3 Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . 109
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
TABLE DES MATIRES 125
III Bibliographie 111
7 Production scientifque 113
7.1 Tableau rcapitulatif . . . . . . . . . . . . . . . . . . . . . . . . . 113
7.2 Actes de confrences dites, confrences invites et chapitres
de livres . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 113
7.3 Revues internationales avec comit de lecture . . . . . . . . . . 114
7.3.1 Confrences internationales avec comit de lecture . . . 115
7.4 Communications . . . . . . . . . . . . . . . . . . . . . . . . . . . 119
8 Rfrences bibliographiques 121
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
126 TABLE DES MATIRES
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Liste des tableaux
4.1 Paramtre de la machine optimiser et leurs bornes . . . . . . . 49
4.2 Contraintes de loptimisation sur les valeurs de sortie . . . . . . 49
4.3 Paramtres des machines suite aux optimisations individuelles . 50
4.4 Paramtres des machines aprs optimisation classique de la
gamme : optimisation de la machine de plus forte puissance, puis
dclinaison pour les machines de plus faible puissance . . . . . 51
4.5 Pondrations des ventes pour les trois scnarios envisags . . . 52
4.6 Paramtres optimaux des trois types de machine dans le cas du
scnario I . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
scnario II . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
scnario III . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
4.9 Notations utilises pour le convertisseur abaisseur . . . . . . . . 57
4.10 Dhnitions et domaine de variation des paramtres de loptimi-
sation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
4.11 Contraintes de loptimisation du convertisseur abaisseur . . . . 64
127
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
128 LISTE DES TABLEAUX
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
Table des fgures
1 Amazing Stories aot 1926 . . . . . . . . . . . . . . . . . . . . . 2
1.1 Amazing Stories aot 1926 . . . . . . . . . . . . . . . . . . . . . 19
2.1 Amazing Stories septembre 1926 . . . . . . . . . . . . . . . . . . 25
3.1 Amazing Stories juillet 1926 . . . . . . . . . . . . . . . . . . . . 42
4.1 Reprsentation schmatique dune chane de traction lectrique 44
4.2 Variation des masses optimales du systme en fonction de la
contrainte sur le rapport de transformation . . . . . . . . . . . . 45
4.3 Gomtrie de la machine tudie . . . . . . . . . . . . . . . . . . 47
4.4 Organigramme de calcul du couplage lectriquethermique . . . 48
4.5 Front de Pareto, cot en fonction du rendement pour la machine
de 125kVA. Les rendements des autres machines de la gamme
sont reports sur la hgure. . . . . . . . . . . . . . . . . . . . . . 55
4.6 Schma lectrique de principe du hacheur . . . . . . . . . . . . 57
4.7 Forme donde du courant dans linductance et le MOSFET . . . 58
4.8 Ondulations de la tension et du courant de sortie . . . . . . . . . 59
4.9 Modle lectro-thermique . . . . . . . . . . . . . . . . . . . . . . 61
4.10 Valeurs du rendement dans le plan puissancetension de sortie . 62
4.11 Dhnition de la gomtrie du disspipateur . . . . . . . . . . . . 62
4.12 Front de Pareto pour loptimisation multi-objectif volume
rendement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
4.13 Calculs successifs du front de Pareto . . . . . . . . . . . . . . . . 66
4.14 Tempratures de jonction des diode et du MOSFET le long du
front de Pareto . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
4.15 Frquence de commutation le long du front de Pareto . . . . . . 67
4.16 Amazing Stories juin 1926 . . . . . . . . . . . . . . . . . . . . . 69
5.1 Schma du circuit imprim tudi . . . . . . . . . . . . . . . . . 73
5.2 Photo recto (a) et verso (b) du circuit imprim tudi aprs in-
fammation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
129
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
130 TABLE DES FIGURES
5.3 volution du courant dalimentation et du courant de fuite en
fonction du temps (haut) et de la temprature du PCB (bas) . . . 75
5.4 Cartographie de temprature et prohls le long de la piste . . . . 76
5.5 Images enregistres par la camra rapide au cours de linfam-
mation dun circuit imprim de 0,8mm . . . . . . . . . . . . . . 77
5.6 Gomtrie utilise pour le modle numrique . . . . . . . . . . 78
5.7 Distribution de la temprature au cours du temps dans un sub-
strat z = 0,2 mm pour une intensit de courant de chaune com-
prise entre 10A et 20A et un coemcient dchange thermique gal
8W.m
2
.K
1
. . . . . . . . . . . . . . . . . . . . . . . . . . . 78
5.8 Comparaison entre les prohls de temprature exprimentaux et
ceux calculs numriquement transversalement au dfaut, pour
une valeur de h gale 8W/m
2
/K . . . . . . . . . . . . . . . . . 79
5.9 Courbes dvolution en fonction du temps de chaune de la tem-
prature la surface dun substrat pour un courant de chaune
initial gal 4 A - Comparaison entre les rsultats exprimen-
taux et le modle numrique . . . . . . . . . . . . . . . . . . . . 80
5.10 Modle dune distribution de la densit de courant selon laxe z
un instant t dans un substrat . . . . . . . . . . . . . . . . . . . 81
5.11 Amazing Stories avril 1926 . . . . . . . . . . . . . . . . . . . . . 82
5.12 Description dune dcharge pointe-plan . . . . . . . . . . . . . 86
5.13 Comparaison du courant de dcharge et de lmission lumineuse
la camra streak, simulation gauche et expriences droite . 89
5.14 Photographie de lenceinte de dcharge . . . . . . . . . . . . . . 93
5.15 Schma lectrique de lalimentation utilise . . . . . . . . . . . . 93
5.16 Photographie en pause longue dune dcharge luminescente glow 94
5.17 Photographie en pause longue dune transition glow-arc . . . . 95
5.18 Oscillogramme du courant et de la tension et photographies des
phases de la transition. . . . . . . . . . . . . . . . . . . . . . . . 96
5.19 Oscillogramme du courant et de la tension pour un chelon de
280mA, apparition de transitions spontanes . . . . . . . . . . . 97
5.20 Zoom de loscillogramme du courant et de la tension pour un
chelon de 280mA, apparition de transitions spontanes et pho-
tographies camra rapide de la dcharge . . . . . . . . . . . . . 98
5.21 Caractristique courant-tension associe des images camra
rapide dune dcharge glow . . . . . . . . . . . . . . . . . . . . . 99
5.22 Caractristique courant-tension pour une distance intrlec-
trode de 25mm . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
5.23 Courant et tension de dcharge lors de la transition (extension
de la colonne hlamentaire . . . . . . . . . . . . . . . . . . . . . . 101
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
TABLE DES FIGURES 131
5.24 Images de la camra rapide en rapport avec les intervalles de la
hgure 5.23 lors dune transition . . . . . . . . . . . . . . . . . . 101
5.25 Dure de la transition en fonction de la distance interlectrode . 102
5.26 Tension et courant de dcharge pour une srie de transitions
glow-arc et arc-glow spontanes et rversibles . . . . . . . . . . 103
5.27 Amazing Stories avril 1926 . . . . . . . . . . . . . . . . . . . . . 104
6.1 Amazing Stories dcembre 1926 . . . . . . . . . . . . . . . . . . 110
7.1 Amazing Stories septembre 1926 . . . . . . . . . . . . . . . . . . 120
8.1 Amazing Stories aot 1926 . . . . . . . . . . . . . . . . . . . . . 132
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2
FicUvr 8.1 Amazing Stories aot 1926
t
e
l
-
0
0
7
1
8
5
5
5
,

v
e
r
s
i
o
n

1

-

1
7

J
u
l

2
0
1
2

Philippe Dessante HDR

Transféré par

Droits d'auteur :

Formats disponibles

Philippe Dessante HDR

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

Philippe Dessante HDR

Transféré par

Droits d'auteur :

Formats disponibles

N DORDRE 1515

C ahn dassurer un fonc-

C. En enet en dessous de cette

FicUvr 4.14 Tempratures de jonction des diode et du MOSFET le long du front

FicUvr 4.15 Frquence de commutation le long du front de Pareto

est le coemcient de photoionisation dpendant de

est la distribution de la densit

E donn par la rsolution de lquation de Poisson :

Vous aimerez peut-être aussi