Prise en compte du temps local dans la modélisation des ceintures de radiation terrestres
Auteur / Autrice : | Damien Herrera |
Direction : | Vincent Maget, Sébastien Bourdarie |
Type : | Thèse de doctorat |
Discipline(s) : | Astrophysique, sciences de l'espace, planétologie |
Date : | Soutenance le 09/10/2017 |
Etablissement(s) : | Toulouse, ISAE |
Ecole(s) doctorale(s) : | École doctorale Sciences de l’univers, de l’environnement et de l’espace (Toulouse) |
Partenaire(s) de recherche : | Laboratoire : Office national d'études et recherches aérospatiales. Département Physique, Instrumentation, Environnement, Espace (DPHY) (Toulouse, Haute-Garonne) |
Equipe de recherche : Équipe d'accueil doctoral Physique spatiale et instrumentation (Toulouse) | |
Jury : | Président / Présidente : Fabrice Mottez |
Examinateurs / Examinatrices : Bénédicte Escudier | |
Rapporteur / Rapporteuse : Dominique Delcourt, Viviane Pierrard |
Mots clés
Mots clés contrôlés
Résumé
Depuis le début de l’ère spatiale avec le lancement du satellite Spoutnik 1 en 1957, les ceintures de radiationterrestres n’ont cessé de faire l’objet d’études du fait de leur dangerosité pour les satellites mais aussi pour l’êtrehumain. En effet, lors d’une forte activité solaire, l’injection de particules dans cet environnement radiatif peut induiredes flux jusqu’à 1000 fois plus élevés que par temps calme. Par conséquent, il est important d’en comprendrela physique ainsi que la dynamique au cours de ce que l’on appelle un orage géomagnétique. Dans ce but, le Département Physique Instrumentation Environnement et Espace (DPhIEE) de l’ONERA développe depuis maintenantplus de 20 ans la famille de modèles Salammbô reproduisant de façon robuste et en trois dimensions la dynamiquedes particules piégées dans ces ceintures. Néanmoins, bien que précis au-delà d’environ 100 keV, la physique et leshypothèses prises en compte dans ce modèle restent insuffisantes en deçà. En effet, aux basses énergies, les ceintures de radiation ne peuvent plus être considérées comme homogènes autour de la Terre. L’objectif de cette thèse a donc été de prendre en compte une quatrième dimension, le temps magnétique local (MLT), afin de mieux reproduire l’évolution des structures fines lors d’un orage géomagnétique. La première partie s’est portée sur l’optimisation du schéma numérique. L’ajout d’une quatrième dimension induit, via l’apparition d’un terme d’advection, une forte diffusion numérique qu’il convient de limiter, tout en tenant compte du temps de calcul. L’équation statistique implémentée a alors été discrétisée selon un schéma de type Beam-Warming du second ordre couplé à un limiteur Superbee, garantissant une propagation satisfaisante de la distribution initiale. Une fois les problèmes numériquesmaitrisés, les différents mécanismes physiques pilotant la dynamique des particules piégées ont été implémentésdans le code, avec une attention toute particulière sur la dépendance en MLT de l’interaction onde-particule. Laprise en compte des champs électriques magnétosphériques fut également nécessaire. En effet, ils constituent l’undes moteurs principaux du mouvement des particules de basses énergies. Le modèle Salammbô 4D a ensuite étévalidé par comparaison avec le modèle 3D déjà existant sur une simulation de l’orage magnétique de Mars 2015.Les résultats ont montré une bonne restitution de la dynamique des ceintures de radiation, avec en plus l’accès à laphase principale de l’orage. Cet évènement a ensuite été modélisé à plus basse énergie pour constater la dynamiqueasymétrique des électrons piégés avec le rôle prépondérant du champ électrique de convection. La comparaison avecles données du satellite THEMIS a montré une bonne modélisation des différents processus physiques, notammentcelui de « dropout » par traversée de la magnétopause. Enfin, la mise en place d’une condition limite dynamiquemodulée par les paramètres du vent solaire et dépendante du MLT ouvre de nombreuses perspectives.