〰️ Récurrence de Fermi-Pasta-Ulam
Simulez le modèle β de FPU : une chaîne de 32 oscillateurs non linéaires. L'énergie commence dans le mode 1, se répand vers les modes supérieurs, puis revient — le célèbre paradoxe de récurrence FPU. Observez le spectre d'énergie modale évoluer en temps réel.
Ce que ça démontre
En 1955, Enrico Fermi, John Pasta, Stanislaw Ulam (et la programmeuse Mary Tsingou) ont réalisé l'une des toutes premières simulations informatiques en physique. Ils s'attendaient à ce qu'une chaîne d'oscillateurs faiblement non linéaires atteigne progressivement l'équilibre thermique — l'énergie se répartissant également entre tous les modes normaux. Au lieu de cela, après avoir visité quelques modes supérieurs, l'énergie est revenue presque entièrement au premier mode : la récurrence FPU. Cela contredisait l'hypothèse alors dominante selon laquelle une faible non-linéarité suffit à elle seule à la thermalisation.
Comment utiliser
Réglez β (non-linéarité) à zéro pour une chaîne parfaitement linéaire — le mode 1 conserve toute son énergie pour toujours. Augmentez β pour observer l'énergie se déverser dans les modes 2, 3, 4 … puis revenir. Augmentez l'énergie initiale E₀ pour une non-linéarité plus forte et des récurrences plus rapides (mais moins exactes). Suivez davantage de modes pour voir le spectre complet. La ligne pointillée jaune marque le niveau d'équipartition — la hauteur que chaque barre atteindrait en cas de véritable thermalisation. Cliquez sur Mode 2 pour démarrer avec le mode 2 excité à la place.
Le saviez-vous ?
Le paradoxe FPU a été résolu dans les années 1960 par la découverte des solitons : des ondes non linéaires localisées qui se traversent l'une l'autre sans se briser. Dans la chaîne FPU, l'énergie se concentre en structures proches des solitons qui la ramènent périodiquement au mode initial. À une énergie suffisamment élevée, le système finit par se thermaliser — mais sur des échelles de temps bien plus longues que prévu initialement. Le problème a stimulé tout le domaine de la dynamique non linéaire et des systèmes intégrables.
À propos de la récurrence de Fermi-Pasta-Ulam
Cette simulation intègre le modèle β de FPU : une chaîne de N=32 particules reliées par des ressorts ayant des forces de rappel à la fois linéaires et cubiques (β). Les équations du mouvement sont ẍⅰ = (xᵢ₋₁ − xᵢ)(1 + β(xᵢ₋₁ − xᵢ)²) + (xᵢ₋₁ − xᵢ)(1 + β(xᵢ₋₁ − xᵢ)²), intégrées avec un schéma de Runge-Kutta d'ordre 4 pour plus de précision. L'énergie est initialement chargée dans un seul mode normal en déplaçant les particules selon son motif de fonction propre sinusoïdale. Les énergies modales sont calculées via une transformée en sinus des déplacements des particules.
Trois visualisations se mettent à jour simultanément : la chaîne de particules montre les déplacements réels colorés selon la vitesse ; le spectre en barres montre comment l'énergie est répartie parmi les N premiers modes (ligne pointillée jaune = niveau d'équipartition) ; et les courbes d'historique temporel montrent l'énergie de chaque mode sur les 300 derniers échantillons de temps. Quand β = 0, seul le mode 1 possède de l'énergie et elle y reste. Avec β > 0, l'énergie fuit vers les modes 2, 3, 4 … puis revient — la récurrence. La récurrence est détectée quand le mode 1 porte >70 % de l'énergie totale.
Questions fréquentes
Qu'est-ce que la récurrence de Fermi-Pasta-Ulam ?
C'est le phénomène surprenant selon lequel l'énergie injectée dans le mode normal le plus bas d'une chaîne d'oscillateurs non linéaires se répand vers des modes supérieurs mais revient ensuite presque entièrement au premier mode, au lieu de se répartir uniformément (se thermaliser) comme le prédit la mécanique statistique.
Que contrôle le paramètre β ?
β fixe la force de la non-linéarité cubique dans la force du ressort. Pour β=0, la chaîne est parfaitement linéaire et l'énergie ne quitte jamais le mode 1. À mesure que β augmente, le couplage entre les modes s'accroît, l'énergie se répand plus vite, les récurrences deviennent moins exactes, et à très haut β le système peut finir par se thermaliser sur de longues échelles de temps.
Qu'est-ce que l'énergie modale et comment est-elle calculée ?
Un mode normal est un motif d'oscillation collectif de toute la chaîne. Son énergie est calculée en projetant les déplacements et vitesses des particules sur la fonction propre sinusoïdale du mode via une transformée en sinus discrète, puis en évaluant l'énergie cinétique plus potentielle dans cette coordonnée modale.
Pourquoi le résultat FPU a-t-il surpris les physiciens ?
La mécanique statistique classique et le théorème d'équipartition prédisent que les systèmes faiblement couplés devraient atteindre un équilibre thermique, avec l'énergie également partagée entre tous les degrés de liberté. Fermi s'attendait à ce que cela se produise rapidement. Au lieu de cela, l'énergie revenait sans cesse au premier mode, suggérant qu'une faible non-linéarité ne garantit pas la thermalisation aux échelles de temps accessibles.
Que représente la ligne pointillée jaune dans le spectre ?
Elle marque le niveau d'énergie d'équipartition — la hauteur que chaque barre atteindrait si l'énergie totale était également répartie entre tous les modes suivis. Une véritable thermalisation amènerait toutes les barres à ce niveau. La récurrence FPU maintient la barre 1 périodiquement bien plus haute que le niveau d'équipartition.
Comment les solitons ont-ils été liés au paradoxe FPU ?
En 1965, Zabusky et Kruskal ont montré que la limite continue de la chaîne FPU est l'équation de Korteweg-de Vries (KdV), qui admet des solutions solitoniques. Ces ondes localisées voyagent sans dispersion et se traversent l'une l'autre de façon élastique. Dans la chaîne discrète, des structures proches des solitons transportent l'énergie sans la partager, ce qui explique la récurrence.
Que se passe-t-il à très haute énergie ou grand β ?
Avec une forte non-linéarité, le système finit par échapper au régime quasi intégrable des solitons. Les récurrences deviennent incomplètes et irrégulières, et sur de très longues durées de simulation le système s'approche bien de l'équipartition — la véritable thermalisation. L'énergie de ce passage est appelée le seuil FPU ou seuil de stochasticité.
La simulation est-elle numériquement précise ?
RK4 (Runge-Kutta d'ordre quatre) est utilisé avec un pas de temps petit, conservant l'énergie hamiltonienne totale avec une grande précision. Vous pouvez le vérifier en observant le relevé d'énergie totale — il devrait rester quasiment constant tout au long de la simulation.
Pourquoi N=32 particules ?
Le calcul FPU original de 1955 utilisait N=32 ou N=64. Avec 32 particules intérieures, la période de récurrence est assez courte pour être visible en temps réel, et la structure modale est facile à afficher. Moins de particules réduit le nombre de modes et rend le spectre plus grossier ; davantage de particules ralentit le calcul.
Qu'est-ce que le modèle β de FPU par rapport au modèle α ?
Le modèle α ajoute un terme non linéaire quadratique (αΔx²) à la force du ressort ; le modèle β utilisé ici ajoute un terme cubique (βΔx³). Le modèle β est plus simple à analyser car il préserve la symétrie de réflexion x → −x de la chaîne. Les deux modèles montrent le phénomène de récurrence, mais avec des périodes de récurrence et des motifs de couplage de modes différents.
Quel est le lien avec la physique et les mathématiques modernes ?
Le problème FPU a lancé l'étude des systèmes intégrables, du chaos hamiltonien et de la physique computationnelle en tant que discipline. Il a influencé le théorème KAM (Kolmogorov-Arnold-Moser) sur la persistance des orbites quasi périodiques, et s'est relié aux études sur les respirations discrètes, les vagues scélérates et la dynamique des réseaux non linéaires encore actives aujourd'hui.
Une chaîne d'oscillateurs non linéaires montrant l'énergie revenant au mode 1 — le paradoxe de récurrence FPU.
2D · HTML5 Canvas 2D · 60 FPS cible · fonctionne entièrement côté client, sans installation