Les produits de la réaction sont des particules alpha (noyaux d'hélium 4) et des neutrons. La fusion nécessite une température élevée : les atomes sont alors entièrement ionisés, et pour décrire l'état de la matière on parle d'état plasma. Ce plasma est enfermé dans une «boîte » spéciale, appelée tokamak, au moyen d'un champ magnétique. Après une brève étude cinématique de la réaction (partie I), on présente une modélisation de la distribution de pression dans le tokamak (partie II). Dans la partie III, une des techniques utilisées pour chauffer le plasma est analysée. Enfin, la partie IV étudie le problème du confinement des particules chargées par un champ magnétique. Les quatre parties sont indépendantes.

Données numériques :

é é é é é é

Formulaire :

Composantes du gradient en coordonnées cylindriques (

) :

Courbe 3D :

abscisse curviligne,

vecteur unitaire tangent,

vecteur unitaire normal principal,

rayon de courbure :

I. Cinématique de la réaction

On suppose que toute l'énergie libérée par la réaction de fusion se transforme en énergie cinétique des particules créées, et on néglige les énergies cinétiques des particules incidentes.

Calculer en eV l'énergie cinétique du neutron et celle de la particule alpha.
Calculer leurs vitesses respectives.

II. Pression et densité du plasma

Dans cette partie, on suppose le plasma composé uniquement d'ions deutérium et d'électrons, confinés dans le tokamak par le champ magnétique

. Chacune des deux espèces est traitée comme un gaz.

On note la densité volumique d'électrons, leur masse volumique et leur vitesse moyenne. De même, on note la densité volumique d'ions, leur masse volumique et leur vitesse moyenne.

Écrire la densité locale de courant électrique

en fonction des paramètres.
2. On considère que le plasma est localement électriquement neutre en tout point. On note

la pression totale à l'intérieur du plasma, a priori non uniforme.
(a) Préciser les diverses densités volumiques de force agissant dans le plasma; en déduire l'équation satisfaite par le plasma en régime permanent.
(b) Montrer que les lignes de courant et les lignes de champ magnétique sont isobares.
3. Le plasma circule à l'intérieur d'un cylindre d'axe

, de rayon

, pour lequel on choisit les coordonnées cylindriques (

) adaptées (figure 1); un solénoïde crée un champ magnétique (

) uniforme dans le cylindre. Cet ensemble constitue un modèle simplifié du tokamak.

Figure 1

(a) On suppose la densité de courant

uniforme. Déterminer le champ magnétique additionnel créé par ce courant.
(b) En déduire, avec ces hypothèses, le profil de pression

à l'intérieur du cylindre, la pression devenant nulle à la paroi

.
(c) La valeur caractéristique du champ magnétique généré par le solénoїde est

, et la composante additionnelle créée par

est de l'ordre de

. Calculer la pression sur l'axe du cylindre (

).
4. On suppose que chaque espèce chargée se comporte comme un gaz parfait monoatomique.
(a) On a obtenu par chauffage la température

telle que

sur l'axe du cylindre. Calculer

en Kelvin. Évaluer la densité volumique d'ions correspondante.
(b) La condition à satisfaire pour que l'énergie produite par les réactions de fusion soit supérieure à l'énergie consommée par le tokamak est le «critère de Lawson » reliant la densité volumique

d'ions, leur température

et leur durée de confinement

; numériquement

.
Avec les valeurs numériques précédentes, évaluer la durée minimale de confinement nécessaire pour que le critère de Lawson soit vérifié.

III. Chauffage du plasma

Pour réaliser la fusion, il faut une température très élevée. Nous allons étudier une des techniques de chauffage utilisées.

Dans un plasma, une onde longitudinale de pression (du type onde « sonore ») est souvent accompagnée d'un champ électrique longitudinal ; ce type d'onde, sans champ magnétique, est appelée «onde électrostatique ». Soit une telle onde plane, donnée par son potentiel . On notera la célérité de l'onde. Le choix des origines est tel que .
(a) Dans le référentiel galiléen R du laboratoire, écrire l'équation du mouvement d'une particule chargée, de coordonnées ( ), en présence de cette onde.
(b) À l'aide d'un changement de référentiel adéquat, ramener cette équation à l'équation du mouvement dans un champ de forces indépendant du temps. En déduire une constante du mouvement, correspondant à l'énergie dans le nouveau référentiel R'.
(c) Préciser dans R' les types de mouvement possibles et décrire les différentes formes du portrait de phase auxquelles elles correspondent.
Préciser en particulier les équations des courbes qui séparent dans le plan de phase ces divers types de mouvement; on posera .
(d) Montrer au moyen du portrait de phase que les particules dont les vitesses initiales selon dans R sont comprises dans l'intervalle peuvent être «piégées» par le champ électrique. Montrer que pour une particule piégée, la valeur moyenne de la composante de sa vitesse dans R , sur un temps suffisamment long, est égale à .
On considère à présent une population de particules chargées. A l'instant initial , leur répartition spatiale est uniforme, et leurs vitesses suivant l'axe sont distribuées selon une gaussienne : la probabilité pour que soit dans l'intervalle est , où désigne la température.
(a) Tracer l'allure du graphe de la fonction .
(b) On «allume» à l'onde électrostatique . On rappelle que la température est proportionnelle au carré de la vitesse quadratique moyenne. On suppose que est petit devant un. En exploitant l'allure du graphe de au voisinage de , expliquer qualitativement pourquoi le piégeage des particules par l'onde chauffe le plasma.
(c) Pour quelle valeur de arrive-t-on à réchauffer le plus de particules?

IV. Confinement magnétique

On considère une particule de charge évoluant dans un champ magnétique indépendant du temps. On note les vecteurs unitaires d'un trièdre trirectangle direct de référence.
(a) Comment évolue l'énergie cinétique de la particule?
(b) On considère un champ magnétique uniforme avec . Déterminer le mouvement de la particule avec les conditions initiales et vitesse . Préciser la pulsation et le rayon de giration.
(c) Les ions et les électrons ont une énergie d'agitation thermique de l'ordre de 10 keV . Donner l'ordre de grandeur de et de pour les ions puis pour les électrons avec .
(d) Dans le plan , en assimilant le mouvement d'une charge à une spire, montrer que le moment magnétique associé s'écrit avec où , étant la vitesse dans ce plan.
(e) Évaluer l'énergie cinétique transverse en fonction de et .
On considère maintenant un champ magnétique non uniforme , avec et de l'ordre du Tesla. On suppose que le champ varie très peu, en valeur relative, sur des distances de l'ordre de . Le mouvement d'une particule de charge comprend alors un mouvement de giration « rapide» orthogonal localement à et un mouvement «lent » de vitesse . On admet que le mouvement lent est celui d'un système de vitesse , portant une charge et un moment magnétique .
(a) La force s'exerçant sur un dipôle magnétique s'écrit

Pour

, montrer que cette force s'écrit également

.
(b) On effectue l'hypothèse que

est parallèle à

, soit

. Donner l'équation différentielle que doit satisfaire

.
(c) On désigne par

l'énergie cinétique associée au mouvement lent.

Montrer que

.
On rappelle que pour toute fonction

.
(d) Soit

la vitesse du mouvement de giration. On admet que

est donné par la même expression qu'en 1.d. Effectuer un bilan global d'énergie cinétique et en déduire que

est une constante.
(e) L'hypothèse effectuée en 2.b est une excellente approximation. Cependant, l'équation du mouvement ne peut être satisfaite en général que si

comporte une composante

orthogonale à

. Compléter alors l'équation du mouvement obtenue en 2 .b pour en tenir compte.
(f) En déduire, en supposant les variations temporelles de

négligeables, que

est donné par :

où

est le rayon de courbure de la ligne de champ et

le vecteur normal unitaire. Indication : on calculera

et on utilisera l'expression de

obtenue en 1.b.
3. Pour confiner les particules chargées composant le plasma, une idée naturelle est de fermer les lignes de champ. Dans ce but on utilise un solénoïde torique. On utilisera les coordonnées cylindriques (

) de la figure 2.

Figure 2

(a) Soit

le nombre total de spires; calculer le champ magnétique créé par le solénoïde lorsqu'il est parcouru par un courant continu d'intensité

.
(b) Expliciter les deux termes dans l'expression (1) donnant

en fonction de

, contributions de

à l'énergie cinétique en notant

; expliquer pourquoi un tel champ ne peut confiner les particules indéfiniment.
(c) Écrire les équations du mouvement d'une particule chargée de coordonnées (

) dans ce champ magnétique torique.
(d) Montrer que l'équation correspondant à

conduit après intégration à une relation qui exprime la conservation d'une grandeur physique que l'on précisera.
(e) Montrer qu'après une intégration de l'équation correspondant à

, on peut se ramener dans une seconde étape à une intégrale première du mouvement ne portant que sur

. A la conservation de quelle grandeur physique correspond-elle?
(f) Dans le cas particulier où les particules circulent à

constant, retrouve-t-on la direction de

donnée en 2.f?
4. Un courant de densité

(sur la direction orthoradiale

) est généré à l'intérieur du tore. On suppose que le « rapport d'aspect »

est grand devant

étant le rayon moyen du tore et

le rayon de sa section (figure 2). On assimile localement le tore à un cylindre, et on suppose pour une approche qualitative que

est uniforme.
(a) Dans ces approximations, donner sans calculs la forme des lignes de champ à l'intérieur du solénoïde.
(b) Dans cette configuration, le confinement des particules est-il amélioré? on demande un argument qualitatif.