
Numerical simulation is a computational method used to model and analyze physical systems mathematically.La simulation numérique remplace les équations différentielles continues décrivant la physique par des schémas aux différences, où le domaine de calcul est divisé en une grille discrète ou un nuage de particules. À chaque pas de temps, les valeurs des champs (densité, vitesse, température) sont mises à jour selon des formules approchées respectant les lois de conservation. Cette approche permet d'étudier des systèmes non linéaires comme la turbulence ou la formation des galaxies, pour lesquels il n'existe pas de solutions analytiques.
Les origines remontent aux calculs numériques manuels d'orbites effectués par les astronomes du XIXe siècle. L'impulsion décisive est venue des recherches sur la turbulence et la physique nucléaire au milieu du XXe siècle (travaux de von Neumann et Ulam sur la méthode de Monte-Carlo). En astrophysique, les travaux pionniers de Hénon et Aarseth sur la modélisation des amas stellaires dans les années 1960 ont été novateurs. Henri Poincaré a jeté les bases de la théorie des systèmes dynamiques, cruciale pour comprendre le chaos dans les calculs numériques. Avec l'essor des supercalculateurs dans les années 1980-2000, des simulations cosmologiques à grande échelle (Millennium, EAGLE, Illustris) sont apparues.
Le principe repose sur la discrétisation des équations aux dérivées partielles. Par exemple, en hydrodynamique, on utilise la loi de conservation de la masse, de la quantité de mouvement et de l'énergie, écrite sous la forme ∂U/∂t + ∇·F = S, où U est le vecteur des variables conservatives, F les flux, S les sources. Le cœur de l'algorithme est la condition de Courant-Friedrichs-Lewy (CFL), qui assure la stabilité : le nombre de Courant C = c·Δt/Δx ≤ 1, où c est la vitesse de propagation des perturbations, Δt et Δx les pas de temps et d'espace. Les codes modernes utilisent un raffinement adaptatif du maillage et une parallélisation sur des milliers de cœurs. Pour les problèmes de gravitation, on utilise la transformée de Fourier rapide ou des méthodes hiérarchiques en arbre.
Création de prototypes virtuels en aviation, prévision climatique, découverte des ondes gravitationnelles par comparaison de signaux réels avec des modèles simulés.
Dans le projet IllustrisTNG, l'évolution d'un cube de 300 millions d'années-lumière de côté, contenant des millions de galaxies, a été reproduite, et les scientifiques peuvent « pénétrer » à l'intérieur de n'importe laquelle d'entre elles.






























