Chauffage Compressionnel de l'Environnement des Disques ...
Chauffage Compressionnel de l'Environnement des Disques ...
Chauffage Compressionnel de l'Environnement des Disques ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Annexe B<br />
Sillage MHD d’un cylindre<br />
tel-00011431, version 1 - 20 Jan 2006<br />
Cette annexe s’attache à fournir un support plus rigoureux et plus général aux conclusions<br />
<strong>de</strong> la section 3.2.3. Nous montrons comment, grâce à un développement incompressible, on<br />
peut déterminer l’écoulement MHD autour d’un cylindre infini, quel que soit le régime <strong>de</strong><br />
viscosité.<br />
B.1 Principe général et développement incompressible<br />
Nous reprenons pour cela, la même situation que celle décrite en section 3.2.3 : un cylindre<br />
infini d’axe parallèle au champ magnétique moyen. Un flot arrivant avec une vitesse v c dans<br />
le repère du cylindre, loin en amont du cylindre.<br />
L’écoulement qui nous intéresse est la solution du système d’équations <strong>de</strong> la MHD en<br />
présence d’une viscosité <strong>de</strong> compression : l’équation d’Euler, l’équation <strong>de</strong> conservation <strong>de</strong> la<br />
masse, l’équation d’induction et une équation <strong>de</strong> fermeture type polytrope.<br />
L’invariance le long <strong>de</strong> l’axe vertical simplifie l’expression <strong>de</strong> la force visqueuse <strong>de</strong> compression<br />
et <strong>de</strong> la force <strong>de</strong> Lorentz. La force visqueuse s’écrit en effet simplement :<br />
⃗F η = 1 3 ⃗ ∇ ⊥<br />
(<br />
η 0<br />
⃗ ∇⊥ .⃗v ⊥<br />
)<br />
(B.1)<br />
Dans cette géométrie, il ne peut pas y avoir <strong>de</strong> champ magnétique perturbé perpendiculaire<br />
au champ <strong>de</strong> référence. Les lignes <strong>de</strong> champ restent verticales. Elles se contentent <strong>de</strong> se<br />
compresser et d’éviter le cylindre. La force <strong>de</strong> tension <strong>de</strong>s lignes <strong>de</strong> champ est donc nulle et<br />
la force <strong>de</strong> Lorentz s’exprime simplement comme le gradient <strong>de</strong> la pression magnétique :<br />
⃗F M = ⃗j × ⃗ B = 1 2 ⃗ ∇B 2<br />
(B.2)<br />
Par ailleurs, la composante verticale <strong>de</strong> l’équation d’induction donne une relation simple<br />
entre le champ et la <strong>de</strong>nsité :<br />
(⃗v. ∇) ⃗ ∇(B/ρ) ⃗ = 0<br />
(B.3)<br />
Si on impose une <strong>de</strong>nsité et un champ constants loin en aval <strong>de</strong> l’obstacle, alors naturellement<br />
le champ magnétique est partout proportionnel à la <strong>de</strong>nsité et les forces <strong>de</strong> pression et <strong>de</strong><br />
Lorentz peuvent simplement se regrouper pour donner le terme :<br />
⃗∇P tot = v 2 m ⃗ ∇ρ<br />
(B.4)<br />
où vm 2 = c 2 s + vA 2 est la vitesse <strong>de</strong> propagation <strong>de</strong>s on<strong>de</strong>s magnétosonores rapi<strong>de</strong>s qui, dans<br />
cette géométrie, ne peuvent se propager que perpendiculairement au champ magnétique.