La solution numérique est obtenue pour un nombre de Reynolds magnétique Rm=1000, Le temps est mesuré en unités du temps de retournement L/U; un tour complet sur la ligne d'écoulement la plus rapide prend un temps t/pi=0.532.
La solution, obtenue ici pour omega=1 and eta=0.1, correspond au trait jaune. Les deux courbes vertes correspondent à l'enveloppe exponentielle multipliant la partir purement oscillatoire Les traits verticaux bleu indiquent l'épaisseur de peau électromagnétique, tel que donnée par l'équation (7.43), mesurée à partir de y=0.
Solution numérique semblable à celle présentée plus haut (Rm=1000), sauf qu'ici on utilise cette fois des conditions frontières périodiques en x et y. La solution est maintenant suivie sur 5 temps de retournement.
La solution suivante est identique à la précédente, sauf que cette fois la condition est telle que la valeur moyennée du potentiel vecteur sur chaque ligne d'écoulement diffère d'une ligne d'écoulement à l'autre. Il s'agit encore une fois d'une solution à Rm=1000, suivie maintenant sur 20 temps d'advection L/u. La dissipation du champ magnétique se produit selon les trois étapes décrites à la section 7.3.4 des notes. La phase "rapide" se produit dans l'intervalle temporel t/pi=0.3-3, et la phase "lente" dans l'intervalle t/pi=3-300. Notez comment, après t/pi=3, les lignes de champs magnétiques se retrouvent alignées avec les lignes d'écoulement (dont la forme est tracées via les "bouchons" verts).
La première animation correspond à une inclinaison de 60 degrés (comme sur la Figure 7.12 des notes), et la seconde à un dipole perpendiculaire (inclinaison 90 degrés). Notexz comment, dans ce dernier cas, le champ magnétique est entièrement détruit. Dans les deux cas
-Last Revised 20 octobre 2008 by paulchar@astro.umontreal.ca.