Pour les épreuves d'admissibilité, l'usage de calculatrices électroniques de poche à alimentation autonome, non imprimantes et sans document d'accompagnement, est autorisé, une seule à la fois étant admise sur la table ou le poste de travail, et aucun n'échange n'est autorisé entre les candidats.

Dans ce problème, on se propose d'étudier la modélisation d'un propulseur plasma susceptible d'être embarqué sur des satellites comme moteur d'appoint de contrôle d'attitude. Le principe de ce propulseur est d'injecter un gaz neutre à la rencontre d'un faisceau d'électrons ionisant. Les ions produits sont fortement accélérés et leur vitesse de sortie du propulseur fournit l'impulsion et donc la poussée nécessaire au satellite.
Pour que l'ionisation soit efficace, il est nécessaire de confiner les électrons par un champ magnétique produit par un dispositif qui ne sera pas étudié dans ce problème.

Les différentes parties du problème sont relativement indépendantes, à condition d'admettre les expressions données dans l'énoncé. Dans la partie I. on étudie le rôle essentiel du champ magnétique appliqué dans le processsus d'ionisation des atomes neutres de Xénon injectés dans le moteur. Une première modélisation simple du fonctionnement en régime stationnaire est abordée dans la partie II. et conduit aux courbes de la figure 4. Les insuffisances de cette modélisation conduisent alors à une amélioration du modèle qui fait l'objet de la partie III. Compte tenu des résultats obtenus par simulation figure 5 , on développe dans la partie IV. une approche asymptotique simplifiée afin d'obtenir l'expression analytique des paramètres caractéristiques importants. Pour garantir un bon fonctionnement du propulseur, il est important que le flux d'ions soit stable et ne subisse pas, notamment, d'oscillations temporelles, fréquentes dans les plasmas. On étudiera donc une de ces oscillations en fin de problème dans la partie .

Le schéma d'ensemble du propulseur est représenté figure 1.

Sur la figure 2, est représenté un schéma de principe du même dispositif, en coupe (2.a.) et en représentation quasi uni dimensionnelle simplifiée (2.b.) dont le principe sera précisé dans la partie .

On distingue essentiellement trois zones, en allant de la droite vers la gauche :

région d'ionisation des neutres par les électrons

Figure 2

Un des paramètres cruciaux caractérisant l'efficacité d'un propulseur spatial est la vitesse de sortie du fluide propulsif par rapport au moteur.

Un dispositif à plasma permet d'atteindre typiquement u de l'ordre de alors que les propulseurs classiques à combustion chimique ne permettent que d'atteindre des vitesses d'éjection dix fois plus faibles.
5. Sachant que dans les moteurs à plasma étudiés ici les débits de masse sont de l'ordre de la dizaine de milligrammes par seconde, expliquer pourquoi ces moteurs sont utilisés pour les contrôles d'attitude des satellites, mais pas pour le lancement et la mise sur orbite initiale.

Les données numériques nécessaires aux différentes applications demandées sont regroupées cidessous.

On considère ici le réacteur, de longueur d, modélisé de façon quasi-bidimensionnelle comme suit : la propagation des neutres et des ions produits se fait suivant l'axe , le champ électrique étant aussi dirigé suivant :

Le champ magnétique radial est supposé ici uniforme et dirigé suivant , la direction orthoradiale étant ici modélisée comme deuxième direction des coordonnées cartésiennes (Oy).
I.1. À quelle condition sur les dimensions du réacteur cette approximation est-elle justifiée?

Dans toute la suite du problème on se place dans un cadre non-relativiste et on considère le référentiel lié au satellite comme galiléen pendant la durée de transit des électrons dans le dispositif.
I.2. On considère d'abord le mouvement des électrons en l'absence de champ magnétique. En plus de la force électrique, ils sont soumis à une «force de frottement visqueux» d'expression où m est la masse des électrons, ū leur vitesse et .
a. Rappeler la signification d'un tel teme dans l'équation du mouvement. À quoi correspond la durée ? En régime permanent on peut alors écrire où est la mobilité. Exprimer en fonction de et .
b. On suppose que les électrons se meuvent dans un champ électrique uniforme crée par une différence de potentiel de 300 V imposée sur la distance . Les électrons sont injectés dans le réacteur en (par la sortie) avec une vitesse négligeable. Déterminer la valeur numérique de et de la vitesse atteinte en régime permanent. Si on suppose au contraire que le terme de frottement est négligeable, quelle est la vitesse maximale atteinte par les électrons, et la durée de leur traversée du réacteur (de jusqu'à ) ? Application numérique. Que conclure des résultats précédents quant à l'efficacité du processus d'ionisation dans une telle situation?
I.3. On considère à présent le mouvement d'ions de charge te ou d'électrons dans un champ uniforme et permanent, en l'absence de champ électrique, les particules arrivant dans la zone de champ perpendiculairement à ce champ, avec une vitesse de module .
a. Rappeler sans démonstration la nature des trajectoires.
b. Exprimer le rayon R de ces trajectoires en fonction de (ici m peut être suivant les cas la masse des électrons ou des ions).
c. En supposant que est la vitesse atteinte par les particules si elles sont accélérées par le champ électrique de la question précédente sur la distance d. en déduire que oú K est une constante dont on précisera l'expression et m la masse des particules. Application numérique: , calculer pour les ions de Xénon. et pour les électrons.
d. Conclure quant à l'influence du champ magnétique sur les ions dans les conditions de l'expérience. Quelle est la nouvelle contrainte qui doit être vérifiée dans le cadre de la question I.1. ? Est-ce le cas ?
I.4. On s'intéresse à présent au mouvement des électrons en présence des champs et croisés (figure 3) et en tenant compte du terme de relaxation collisionnel (désormais et jusqu'à la fin du problème, m désigne la masse de l'électron). Le champ des vitesses électroniques est décrit par en régime stationnaire. On néglige en première approximation l'accélération convective des électrons dans les équations du mouvement. On pose .
a. Écrire l'équation du mouvement et la projeter sur les trois axes de coordonnées.
b. En deduire , une relation entre les composantes et , et montrer que l'on peut en déduire: avec
c. Application numérique : calculer et comparer et . Calculer l'ordre de grandeur de et de .
d. Conclure quant à l'efficacité du processus d'ionisation dans cette nouvelle configuration.
e. En analysant les équations du mouvement en l'absence de collisions, expliquer le rôle déterminant du couplage champ électrique-champ magnétique pour piéger les électrons tout en accélérant les ions.
f. On cherche à vérifier à présent la validité de l'approximation consistant à négliger l'accélération dans les équations du mouvement des électrons. Exprimer l'accélération en fonction du champ des vitesses et de ses dérivées spatiales. Évaluer un ordre de grandeur des termes de l'accélération suivant et ; conclure sur la validité de l'approximation.
g. Le terme collisionnel traduit l'influence des collisions des électrons sur les neutres et doit donc faire intervenir la vitesse relative ( ) au lieu de . Montrer, toujours par une étude d'ordre de grandeur, que ceci ne modifie pas de façon notable l'étude ci-dessus.

On cherche ici à retrouver directement par une étude microscopique simplifiée le mouvement de dérive des électrons dans les champs croisés et dans le cas particulier où . Les champs et sont ici toujours supposés uniformes et perpendiculaires et on ne s'intéresse qu'au mouvement d'un électron entre deux collisions, projeté perpendiculairement à . On note sa vitesse perpendiculairement à sa valeur juste après une collision avec un neurre. On note l'intervalle de temps moyen entre deux collisions. On décompose en où ne dépend pas du temps et satisfait à une équation où intervient plus le champ électrique.
a. Exprimer alors en fonction de et .
b. On introduit la variable complexe avec: . Montrer que et en déduire l'expression de puis celle de en prenant pour origine la position juste après la demière collision.
c. En déduire movennant quelques hypothèses que l'on précisera, dans le cas , que le déplacement moyen entre deux collisions peut s écrire
d. Retrouver alors les expressions de et déterminées en I.4.b. (dans le cadre de l'approximation ).

Les équations qui régissent le comportement des électrons, des ions et des neutres ainsi que l'expression du champ électrique sont couplées et relativement complexes. On s'intéresse ici à une mise en équation simplifiée conduisant par une résolution numérique à une première solution décrite par les courbes de la figure 4 . On note la densité particulaire des ions Xénon monochargés et celle des électrons.
Les hypothèses de ce premier modèle sont les suivantes :
(i) quasi-neutralité du plasma ( ) vérifiée localement
(ii) toutes les grandeurs ne dépendent que de l'abcisse x
(iii) la distribution de champ magnétique est imposée de l'extérieur (figure 4.f)
(iv) les neutres, émis à l'anode en , ont une vitesse uniforme et constante parallèle à l'axe du réacteur Ox
(v) les ions positifs sont crées avec une vitesse initiale par impact des électrons sur les neutres. Leur vitesse reste alors parallèle à l'axe Ox et ils ne subissent aucune collision après leur création. Leur masse est notée M (ainsi que celle des neutres)
(vi) le mouvement des électrons est donné par les équations déterminées en I.4.b., en particulier selon l'axe
(vii) le nombre de collisions efficaces (conduisant à ionisation) par unité de volume et de temps est décrit par où est une fonction connue (non précisée ici) de la température des électrons décrivant la collision au niveau microscopique.

II.1. On note et les flux surfaciques de chaque espèce (ou densités volumiques de courant particulaire) comptés algébriquement suivant l'axe orienté Ox .
a. Exprimer en fonction de et .
b. En effectuant un bilan sur les neutres, déterminer une équation différentielle reliant et .
c. Relier les variations spatiales de et , et en déduire l'expression de et à l'aide de et , et des conditions aux limites et (on précisera la valeur choisie pour au vu des hypothèses du modèle).
La connaissance de nécessite alors celle de et de .
II.2. Pour déterminer , on utilise l'hypothèse (v).
a. On s'interesse aux ions produits par unité de temps dans la tranche de section unité et comprise entre les abcisses et . Exprimer la vitesse de ces ions lorsqu'ils parviennent à l'abcisse en fonction de la différence de potentiel , de , e, et M.
b. Relier alors par un bilan que l'on précisera le nombre de ces ions (produits entre et et observés en x ) par unité de volume à , dx et .
c. En déduire que peuts'exprimer par:

II.3. L'énergie des électrons est essentiellement sous forme thermique et vaut en moyenne pour un électron situé au voisinage de l'abcisse x .
a. Justifier qualitativement cette affirmation au vu des questions I.4.c. et I.4.d..
b. On considère que l'évolution de l'énergie est donnée par l'équation : où est une constante phénoménologique (longueur de relaxation). Justifier la façon dont intervient le champ dans cette équation. Que traduit le terme ?

Cette équation permet alors de déterminer à partir de , lui-même obtenu à partir de (vi) et de II.2.c. .
II.4. L'ensemble de la démarche ci-dessus permet alors la résolution du problème stationnaire par une méthode numérique itérative par relaxation dont le détail ne sera pas abordé ici. Les résultats sont représentés par les courbes de la figure 4 (on fera attention aux unités utilisées et à la double échelle, potentiel à gauche, champ électrique à droite, de la figure 4 (a)).
a. L'ionisation est-elle efficace ?
b. Pourquoi passe-t-il par un maximum ?
c. Le mouvement réel des ions n'étant pas tout à fait unidimensionnel, pourquoi choisit-on de localiser la zone de champ magnétique plutôt vers la sortie du réacteur?
d. Commenter la forme de la partie droite de la courbe 4 (b) compte tenu de lallure de : que retrouve-t-on ici comme information concernant les ions?
II.5. On s'intéresse dans cette question à l'étude, à partir des courbes de la figure 4, de la validité de certaines approximations effectuees dans la modelisation.
a. On note . En utilisant une équation de Maxwell, relier aux variations spatiales du champ électrique. En déduire un ordre de grandeur de pour . Compte tenu de la courbe 4 (b) peut-on considérer que l'hypothèse (i) est correctement vérifiée ici?
b. On s'intéresse à l'hypothèse (vi) et on suppose que les électrons sont assimilables à un gaz parfait de température et de densité . Quelle force n'a pas été prise en compte dans le modèle pour l'étude du mouvement des électrons ? Donner son équivalent volumique en fonction de , et de leurs variations spatiales, puis en utilisant les courbes, montrer que l'on ne peut négliger cette force dans la zone .

On est donc conduit à reprendre les équations ci-dessus afin d'améliorer la modélisation et si possible. en même temps, de préciser la forme du bilan énergétique. C'est l'objet de la troisième partie.

On considère à présent que les neutres, les ions et les électrons satisfont aux équations suivantes en régime stationnaire :
(III.1.)

(III.2.)

(III.3.)

(III.4.)

avec

a. A quoi correspondent les équations (III.1.) ? Quelle hypothèse est implicitement conservée ici?
b. Expliquer la forme du premier membre de l'équation (III.3.). À quoi correspond physiquement le dernier terme de cette équation?
c. En utilisant les résultats des parties et II, commenter l'équation (III.4.) : on identifiera l'origine de chaque terme, on précisera quelles sont les approximations qui ont été effectuées (termes manquants) et pourquoi. Pourquoi appelle-t-on fréquence effective de diffusion axiale?
III.2. Afin de pouvoir réaliser un bilan énergétique correct concernant les électrons en utilisant le premier principe de la thermodynamique, on s'intéresse d'abord à l'érude simplifiée du bilan associé à une ionisation. Lors d'une ionisation, un électron rencontre un neutre et provoque la création d'un ion et d'un nouvel électron.
a. Si l'électron incident a une énergie de l'ordre de la dizaine d'électronvolts, et l'atome neutre de Xénon une vitesse de l'ordre de , comparer leurs quantités de mouvement.
b. Justifier alors que l'on a considéré en II que l'ion était créé avec une vitesse égale à celle des neutres.
c. En déduire que lors de l'ionisation, l'énergie nécessaire à la création de la nouvelle paire (électron, ion) est essentiellement fournie par l'électron.
d. En comparant alors l'énergie disponible après le choc pour les deux électrons et l'énergie de l'électron incident, montrer que tout se passe comme si l'électron incident perdait et l'électron créé l'était avec une énergie nuile.
e. En déduire que dans le bilan d'énergie, on peut oublier le terme associé aux électrons créés en incluant seulement un terme associé aux pertes d'énergie de la forme - par unité de volume et de temps.
f. En fait. on écrit ce terme de perte où est un coefficient sans dimension de l'ordre de 2 à 3 . En considérant que toutes les collisions n'aboutissent pas à une ionisation, proposer une interpretation de .
III.3. On considère à présent le bilan énergétique issu du premier principe, en négligeant l'énergie cinétique macroscopique des électrons devant leur énergie interne.
a. Quelle est l'expression de la puissance volumique cédée par le champ électrique aux électrons?
b. On considère le système ouvert constitué des électrons présents entre x et évoluant dans l'intervalle de temps ( ). Construire un système fermé adapté que l'on précisera en tenant compte des conclusions de III.2.e. et III.2.f.. En lui appliquant le premier principe, établir :

III.4. Une manipulation fastidieuse mais sans difficulté des équations précédentes permet d'obtenir les dérivées spatiales de chaque variable en fonction de l'ensemble des variables du système. On obtient ainsi par exemple: avec
a. Pourquoi dit-on que correspond à une transition subsonique-supersonique pour les ions?
b. À quelle condition cette transition peut-elle avoir lieu de façon régulière (sans divergence des dérivées spatiales) ?
III.5. La résolution numérique des équations du modèle amélioré fournit alors les courbes de la figure 5 (On notera que la figure 5 (b) modifie les dimensions du réacteur, différentes dans cette simulation de celles utilisées dans la partie II, E désignant le point de sortie, A l'anode de potentiel et le point de passage sonique). Sur la figure .
On remarque que l'ionisation est extrêmement localisée ici, en présence d'un champ qui n'est pourtant négligeable nulle part. Proposer une explication à partir des courbes.
En s'aidant des courbes et , expliquer pourquoi la zone située entre 0 et 15 mm est appelée région de diffusion. Commenter les différences apparaissant entre les deux modèles dans cette région.

Afin de préciser les paramètres importants du modèle, on cherche, en exploitant les courbes, à trouver des expressions simplifiées des variables dans différents domaines, c'est le cadre de l'étude asymptotique menée dans la partie IV.

Compte-tenu des résultats de la figure , l'ionisation est limitée à une couche de faible épaisseur en amont du point où l'écoulement des ions devient supersonique. On peut considérer en outre qu'en amont de cette couche d'ionisation, dans la zone B-D (figure 2 (b)), la fréquence effective de diffusion axiale des électrons est à peu près constante et égale à . Les ions présents dans cette zone ont été produits en aval et donc .
IV.1.
a. Que deviennent alors les équations (III.1.), (III.3.) et (III.5.) compte-tenu de ces approximations? En proposer des intégrales premières.
b. Les densités volumiques de courants particulaires et ionique et électronique sont toutes deux négatives dans la zone B - D considérée ; on note la valeur de en B , la valeur de en et la valeur de la vitesse ionique associée en .
Une étude plus approfondie et non abordée ici montrerait que et que à la transition entre la zone de quasi-neutralité et la gaine très mince chargée existant au voisinage de l'anode. On considère que . Déduire des équations précédentes la relation :

c. En utilisant également l'équation (III.4.) simplifiée, établir alors que :

d. En déduire par élimination de entre ces deux équations que :

précisera en fonction des données en .
En utilisant la relation entre et et l'expression de données en IV.1.b., et en exploitant les courbes de la figure 5 , donner une estimation numérique de en prenant .
e. En intégrant (IV.3.), en déduire, en supposant que les ions sont à vitesse subsonique entre
D et B, que :
(IV.4.)

IV.2. L'estimation numérique effectuée en IV.1.d. montre que . On s'intéresse donc ici aux expressions approchées des différents paramètres dans la zone de diffusion caractérisée par .
a. Quel est le comportement asymptotique de dans cette zone ?
b. En déduire que et tendent vers des valeurs constantes et dans cette zone que l'on exprimera en fonction de et .
c. Comment est alors le champ électrique dans cette zone ?
d. Déduire alors de l'équation (III.4.) simplifiée l'évolution de la pression dans cette zone et comparer à la courbe de la figure . On utilisera dans la suite que, compte tenu de cette figure, .
e. À quoi est dû le mouvement des électrons dans cette zone ? Pourquoi est-elle appelée région de diffusion?

Les mesures expérimentales confirment bien l'existence de cette zone ainsi que ses caractéristiques.
IV.3. Afin de déterminer la position du début de la couche d'ionisation, on considère à présent que 'est plus tout à fait nul, mais varie avec suivant :

tant que reste inférieur à quelques .
a. En utilisant les courbes, peut-on considérer que reste constant ?
b. En utilisant alors l'expression obtenue en IV.2.d. pour la pression, en déduire une relation entre et x où interviennent et m .
c. En reprenant le bilan d'énergie (III.5.) et les équations de conservation (III.1.) dans ce cas, établir alors :

d. Déduire de (IV.5.) que :

e. variant beaucoup plus rapidement avec que dans le début de la zone d'ionisation, déduire alors de (IV.6.) et (IV.7.) que l'on a à peu près: où est la valeur de dans la région de diffusion et une expression que l'on précisera.
f. Compte tenu de la courbe caractérisant l'évolution de , en déduire en précisant le raisonnement la position de l'entrée dans la couche d'ionisation en fonction de A et de . On prend . A l'aide des valeurs numériques données en introduction et de la valeur de tirée de la courbe , calculer et et comparer les valeurs obtenues à celles de la figure 5 .

Dans cette partie, on souhaite étudier un exemple d'oscillations qui peuvent naître spontanément dans le plasma et entraver le bon fonctionnement du moteur.

Cet exemple est celui d'oscillations de la région d'ionisation, région d'extension faible en et située à l'abcisse près de la sortie du réacteur (figures 2 (b) et 5 (i)). On suppose toujours valable I'hypothèse de neutralité locale dans le réacteur. On note 1 l'extension typique en x de la zone d'ionisation, .

À quel phénomène est relié ?
2. On veut décrire les perturbations de et de par rapport à leurs valeurs stationnaires d'équilibre et . On écrit alors avec . Déterminer l'expression de à partir des équations (V.1.).
3. En déduire qu'au premier ordre on a les deux équations :