20eme Congres Franc¸ais de Mecanique Besanc¸on aout au septembre

Sophia Antipolis

Description

Niveau: Supérieur
20eme Congres Franc¸ais de Mecanique Besanc¸on, 29 aout au 2 septembre 2011 Simulations numeriques de phenomenes anisotropes pour la Fusion par Confinement Magnetique 1 A. Bonnementa,b, H. Guillarda,b, B. Nkongab,a, R. Pasquettib,a a. INRIA Sophia-Antipolis, 2004 Route des Lucioles, 06902 Sophia Antipolis b. Laboratoire J.A. Dieudonne, UMR 6621, Universite de Nice-Sophia Antipolis, 06108 Nice Resume : Afin d'etudier le comportement du plasma dans le cadre de la FCM, nous considerons un systeme fluide, obtenu a partir du modele cinetique Vlasov-Maxwell avec fermetures de Braginskii, d'ou de fortes anisotropies suivant les lignes de champ magnetique. Des etudes numeriques de phenomenes presentant cette particularite ont ete menees, tout d'abord sur une equation de diffusion, avec un terme de rayonnement, et ensuite sur des systemes de type Euler et Navier-Stokes. Les resultats obtenus sont presentes et analyses. Abstract : In order to study plasma for Magnetic Confinement Fusion, we consider a fluid system, given by the Vlasov-Maxwell kinetic model with the Braginskii closure. This system presents high anisotropy along the magnetic field lines. Some numerical studies have been done in this case, first on a diffusion equation with a radiative term and then on the Euler and Navier-Stokes like systems. Our results are presented and analyzed.

• plasma


• troisieme tenseur de la viscosite de braginskii donnee


• collisions entre particules de differentes especes


• coeur du plasma


• profils de mach et de densite


• lieu dans le plasma de bord des tokamaks


• flux de chaleur


• methode des volumes finis

Sujets

• Antipolis
• Plasma
• Flux thermique
• Méthode des volumes finis

Informations

• Publié par: profil-nechor-2012
• Publié le: 01 septembre 2011
• Nombre de lectures: 29
• Langue: Français
• Poids de l'ouvrage: 1 Mo

Extrait

`eme 20Congr`esFran¸caisdeM´ecanique

Besanc¸on,29aouˆtau2septembre2011

Simulationsnume´riquesdeph´enom`enesanisotropes 1 pourlaFusionparConfinementMagn´etique

a,b a,bb,a b,a A. Bonnement, H. Guillard, B. Nkonga, R. Pasquetti a.INRIA Sophia-Antipolis, 2004 Route des Lucioles, 06902 Sophia Antipolis b.Lrobaiota.Jer6612U,inevsrtie´A.Dieudonn´e,UMR0,silopieciN8016-SceNidentaAhiop

R´esum´e: Afind’´etudierlecomportementduplasmadanslecadredelaFCM,nousconsid´eronsunsyst`eme fluide,obtenua`partirdumod`elecine´tiqueVlasov-MaxwellavecfermeturesdeBraginskii,d’o`ude fortesanisotropiessuivantleslignesdechampmagne´tique.Des´etudesnume´riquesdephe´nom`enes pr´esentantcetteparticularite´ont´ete´men´ees,toutd’abordsurunee´quationdediffusion,avecunterme derayonnement,etensuitesurdessyste`mesdetypeEuleretNavier-Stokes.Lesre´sultatsobtenussont pr´esent´esetanalyse´s. Abstract : In order to study plasma for Magnetic Confinement Fusion, we consider a fluid system, given by the Vlasov-Maxwell kinetic model with the Braginskii closure. This system presents high anisotropy along the magnetic field lines. Some numerical studies have been done in this case, first on a diffusion equation with a radiative term and then on the Euler and Navier-Stokes like systems. Our results are presented and analyzed. Mots clefs : Simulation des plasmas, fermetures de Braginskii, diffusion anisotrope.

1 Introduction LaFusionparConfinementMagn´etique(FCM)s’appuiesurleprincipedeproductiond’´energiedes ´etoiles,lorsdelafusiondedeux´el´ements(parexemple,deuteriumettritium).Pourcela,cesdeux 8 ´ele´mentsdoiventˆetre`adestemp´eratures´eleve´es,del’ordrede10K,soit10keV,afindefavoriserla fusion.Acettetemp´erature,unnouvel´etatdelamati`ereapparaıˆt:leplasma,danslequellesatomes sontionis´es.Lafusionparconfinementmagne´tiqueconsiste`autiliserlechampmagn´etiqueafinde confinerleplasmadansdeschambresappel´eestokamaks.Desphe´nome`nesphysiquescomplexesse d´eroulentdansleplasma,pouvantjouerunroˆled´ecisifsurlesperformancesdestokamaks.Leprincipe deFCMestparexempled´eveloppe´en[1,4]. Danslecadredecetravail,ons’int´eresseplusparticulie`rementauxphe´nome`nesquiontlieudans leplasmadebord,commeparexempleen[9].Danscebut,unmod`elefluideae´t´e´etablia`partir dumode`lecine´tiqueVlasov-MaxwellaveclesrelationsdefermeturedeBraginskii[2].Cemode`leest pre´sente´dansleparagraphe2. Nosr´esultatsnume´riquessontpr´esent´esdanslesparagraphes3et4,cecipourdiffe´rentesge´ome´tries detokamaks(ITER,Tore-Supra).L’ensembledesr´esultatsa´et´eobtenugraˆcea`unem´ethodevol- umes/e´le´mentsfinis(VF/EF),impl´emente´eauseindulogicielFluidBox[11].Cetteapprochea´et´ede pluscompar´ee`aunem´ethoded’e´lementsspectraux(codeSEMd´eveloppe´parR.Pasquetti). 1.Cetravaila´et´eeffectu´edanslecadred’unprojetdelaFR-FCM,del’ANRESPOIRetdel’Actiond’envergure FUSIONdel’INRIA.A.BonnementremercieleConseilRe´gionalPACAquicofinancesathe`se.

1

e`me 20Congr`esFranc¸aisdeM´ecaniqueBesan¸con,29aoˆutau2septembre2011 2Mod`elefluide Lebutesticid’e´tablirlemode`lefluideutilis´epoure´tudierlesphe´nom`enesphysiquesayantlieudans leplasmadeborddestokamaks.Consid´eronsleplasmaconstitue´d’unme´langed’´electronsetd’ions. Danscecas,l’´etatduplasmanousestdonne´graˆcea`lafonctiondedistributionf s (t,x,v), l’indices pouvantsere´fe´rerauxe´lectronsouauxions.Onobtientdoncl’´equationdeBoltzmann: ∂F s f s +v.r x f s +.r v f s =C s ,(1) ∂t m s ou`F s samedeluesstefolacrdeLeroneztxerec´eesuruneparticm s ee´eisalcol,nx, avec une vitessevet au tempst, etC s eunirenortbP.uooisnllsiescoantleerlmiesmolde´tdescription fluidedu plasma, on neconsid`erepaslesfonctionsdedistributionsmaisleursmomentssurl’espacedesvitesses.Ainsi,en prenantlestroispremiersmomentsdel’e´quationdeBoltzmann(1),onobtientdenouvellese´quations pourlesdensit´esdeparticules,n s ,´eesyennomsessetivsel,u s etle,se,reig´snee s :

∂ n s +r ∙(n s u s ) =0 (2.a) ∂t ∂ m s n s u s +r ∙(m s n s <u s ⊗u s >+P s ) =Z s en s (E+u s ∧B) +R s (2.b) ∂t ∂ n e+r ∙(n eu+P .u) =−rE.u+R.u+Q(2.c) s ss ss ss ∙q s +Z s en s ss ss ∂t

Dans ces expressions,R s seluetedcsloilisonsentreparticulseedid´ffreneetesstpl`eescoemse,mtqenr´ui Q s tnescaleelah´gru´eeneer´rlpasp’ee`ecrepr´esapsertuaselucitr,`caedcsuaessionolliclessaveP s est le tenseur de pression etqefsdmeerlareontisee´Braperutnnod,leelruL.seufldxceahi[2]nskiragi s permettentd’obtenirlesformesdecesdiffe´rentstermes. Nosde´veloppementsalgorithmiquespourr´esoudrelesyst`eme(2)sontencourspourlage´ome´triedu tokamakITER.Ci-apre`s,nousnousint´eressons`adesprobl`emesmode`les. 3 Diffusionanisotrope avec rayonnement 3.1Mod`ele Danslestokamaks,l’e´nergies’´evacuantduplasmadecoeuresttr`esimportanteetpeutdoncaffecter lescomposantsdelastructuredutokamak.Deplus,lapre´senced’unrayonnementdoiteˆtrepriseen compte.Nous´etudionslesconse´quencesdetelsprocessusradiatifssurlesprofilsdetempe´rature.Dans lecasuni-dimensionnel,ceprobl`emeae´te´aborde´en[5].Enmode´lisantlerayonnementparunterme depuitsdechaleurnonlin´eaire,onconside`rel’e´quationdetempe´rature: 3∂ 2 2 T+r ∙q=−CT exp(−((T−T)/ΔT(3)) ) 2∂t quiestobtenuea`partirdel’´equationdel’e´nergie(2.c)eto`uTdeuqyareceruitirp´emateresattlon- nement et ΔTep.un´cera-tyt Usuellement,ladensit´edefluxdechaleurqapeelalrdtse´nnoreieonnfdeoiurFoud-ccanodnletcoi tivit´ethermiquedumate´riau,K:q=−KrTduasecslan.Dplasma,noussommeespn´rseneec’dnu champmagne´tique,etladiffusionestalorsanisotrope: q=−K k r k T−K ⊥ r ⊥ T.(4) Les notationsket⊥e`fe´rpsertnermpireauchaneiducalelteepprrapa`ellsapoesntuatnmocxitceemev magn´etique.Onparledediffusionfortementanisotropecarletermeparall`eleestnettementsupe´rieur

2