Physique  |  Technique

 

Raul Alonso De Libano, 2007 | Charmey, FR
Kevin Bourquenoud, 2006 | Attalens, FR

 

Notre travail de maturité porte sur la modélisation et la simulation numérique de la structure interne des naines blanches à l’aide d’une équation d’état de gaz d’électrons dégénérés. Nous avons développé une simulation basée sur la méthode d’Euler, implémentée dans Excel, afin de résoudre ces équations et d’étudier les propriétés globales de ces objets compacts, notamment leur relation masse–rayon et la limite de Chandrasekhar.

Problématique

Dans quelle mesure un modèle numérique simple, basé sur une équation d’état de gaz de Fermi dégénéré, permet-il de reproduire les propriétés fondamentales des naines blanches, telles que leur structure interne, leur relation masse–rayon et leur masse limite ?

Méthodologie

La simulation et la modélisation des naines blanches constituent le cœur de notre travail. Nous avons utilisé Excel, un outil simple et accessible.

Après une phase d’apprentissage théorique notamment grâce à Astrophysics for Physicists, nous avons implémenté la méthode d’Euler pour résoudre numériquement l’équation de Lane–Emden, permettant de déterminer la constante intervenant dans la masse de Chandrasekhar.

Nous avons ensuite développé une simulation complète de la structure interne d’une naine blanche. Le modèle repose sur la densité centrale, à partir de laquelle les équations sont intégrées couche par couche depuis le centre jusqu’à la surface (in-out), en calculant masse, rayon, densité, pression et énergie. La simulation s’arrête lorsque la pression devient nulle. En répétant ce processus pour différentes densités centrales, nous obtenons les relations globales des naines blanches.

Nous avons également étendu le modèle au cas de la rotation en utilisant les pseudo-sphéroïdes de Maclaurin. Cette approche permet de modéliser l’aplatissement de l’étoile et d’introduire l’énergie cinétique de rotation. Nous avons étudié l’équilibre global à l’aide du paramètre bêta (le rapport entre énergie de rotation et gravitationnelle) afin d’identifier des vitesses de rotation critiques.

Résultats

La simulation reproduit correctement les propriétés globales des naines blanches. La relation masse–rayon obtenue est cohérente avec la théorie.

La résolution de l’équation de Lane–Emden donne une valeur de 2.018, en excellent accord avec la théorie. La masse limite obtenue est de 1.456 masses solaires, proche de la valeur attendue (environ 1.44 masses solaire). L’analyse montre que la méthode d’Euler est d’ordre 1 avec des erreurs négligeables pour les pas utilisés.

L’introduction de la rotation met en évidence un aplatissement de l’étoile et une modification de son énergie. Lorsque la vitesse angulaire augmente, l’énergie de rotation devient significative et modifie l’équilibre. Nous mettons en évidence une vitesse critique au-delà de laquelle l’étoile devient instable. La rotation permet également une légère augmentation de la masse maximale supportée.

Discussion

Les résultats montrent qu’un modèle numérique relativement simple permet de reproduire avec une bonne précision les propriétés fondamentales des naines blanches. Malgré l’utilisation d’une méthode d’ordre faible, la précision obtenue reste satisfaisante grâce à un choix adapté du pas d’intégration.
Cependant, certaines limites subsistent. Le modèle ne prend pas en compte les effets relativistes complets, le champ magnétique ou la rotation (en dehors d’une extension partielle de notre modèle). De plus, la précision est limitée par les contraintes techniques d’Excel, notamment le nombre maximal de lignes (1’048’576), ce qui restreint la finesse des simulations.
Enfin, le temps de calcul relativement long (plusieurs heures pour certaines simulations) met en évidence les limites de cet outil pour ce type d’étude, bien qu’il reste pertinent dans une démarche pédagogique.

Conclusions

Ce travail démontre qu’un modèle basé sur des hypothèses physiques simples et une méthode numérique accessible permet de reproduire de manière fiable les propriétés essentielles des naines blanches. La cohérence entre les résultats numériques et les prédictions théoriques valide l’approche adoptée.
Ce projet met également en évidence qu’il est possible d’aborder des problématiques complexes d’astrophysique avec des outils simples, ce qui renforce son intérêt pédagogique.
Des améliorations futures pourraient inclure l’utilisation de méthodes numériques d’ordre supérieur, l’intégration d’effets relativistes complets ou l’étude plus approfondie de la rotation et des étoiles à neutrons.

 

 

Appréciation de l’expert

Prof. em. Georges Meynet

Ce travail a consisté (entre autres) à résoudre l’équation de Lane-Emden pour déterminer la structure interne de naines blanches. Le travail est original par son utilisation d’Excel, explicitant le travail d’intégration des équations. Il est remarquable par son étude des barres d’erreur, liées à la discrétisation spatiale, et par son approche des effets de la rotation. Ce travail montre une très belle maîtrise des questions numériques et physiques par ces gymnasiens. La comparaison avec les observations et des résultats obtenus dans la littérature valide leurs calculs.

Mention:

Argent

Prix spécial «Recherche sur le Jungfraujoch» décerné par l’Académie suisse des sciences naturelles & l’Institut Paul Scherrer

 

 

 

Collège du Sud, Bulle
Enseignant: Dr. Emmanuel Kohlprath