12.01.2015 Views

Chauffage Compressionnel de l'Environnement des Disques ...

Chauffage Compressionnel de l'Environnement des Disques ...

Chauffage Compressionnel de l'Environnement des Disques ...

SHOW MORE
SHOW LESS

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.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!