Tous les travaux

Calcul scientifique

Orbites

Une simulation d'orbites en C, stable dans le temps grâce à un intégrateur symplectique, rejouée dans le navigateur.

RôleAlgorithmique · Visualisation
Année2024
Stack
CEuler symplectiquep5.jsJSON
v ← v + a·dt · x ← x + v·dt

Problème

Simuler des orbites planétaires stables dans le temps. Une intégration naïve (Euler explicite) dérive : l'énergie n'est pas conservée et les planètes spiralent.

Approche

Implémentation en C d'un intégrateur d'Euler symplectique : on met à jour la vitesse, puis la position. Trajectoires exportées en JSON, visualisées en p5.js.

Difficultés résolues

  • Stabilité à long terme : le schéma symplectique garde l'énergie quasi constante.
  • Choix du pas de temps : compromis précision / coût de calcul.
  • Pont C → web : sérialiser les trajectoires en JSON pour les rejouer.

Résultat

Des orbites stables sur de longues durées, explorables dans le navigateur via une visualisation p5.js.