Dans une première partie, chaque étudiant
s’initie à un langage de programmation bien adapté au calcul scientifique. On propose python ou le C++. Voici un
didacticiel pour C++ et voici
un didacticiel pour python. On a choisit ces deux langages, car de nombreuses bibliothèques scientifiques sont disponibles. (pour python il y a scipy, matplotlib,..., pour C++ il y a gsl, root, armadillo... ).
Dans une deuxième partie, en travaillant par binomes (ou seuls), les étudiants choisissent un projet d’expérimentation numérique qu’ils développent. Une liste de projets possibles est donnée ci-dessous. On peut distinguer trois parties dans ce travail:
- Phase de modélisation: partant du problème posé, il faut réfléchir à la façon de modéliser le problème, c’est à dire le mettre en équation, quel algorithme mettre en oeuvre, définir les objectifs que l’on peut obtenir de façon réaliste,...
Très souvent on connait les résultats au problème posé pour des cas particuliers. Il faut bien étudier au préalable ces cas particuliers, car ils permettront de tester le programme.
- Phase de programmation: on écrit le programme dans le langage choisit. (On teste le programme au fur et à mesure en utilisant les cas particuliers)
- Phase d’exploitation: on étudie et discute les résultats de la simulation. On rédige un “document scientifique” (en utilisant Lyx).
Outre la programmation, il pourra être utile de faire du
calcul formel au cours du projet. Pour cela on suggère l’utilisation du logiciel libre
xcas.
1 Systèmes dynamiques avec peu de degrés de liberté
1.1 L’équation du télégraphe en dimension 1.
1.2 Modèles de dynamiques prédateurs, proix et généralisation.
1.3 Application logistique, fractale de Julia et de Mandelbrot.
- Voir Cours.
- Zoomdans la fractale de Mandelbrot.
1.5 Dynamique sur une surface hyperbolique
1.6 Billard polygonal
Est équivalent à
masses libres sur un cercle qui se collisionnent.
Question ouvert: existence d’une orbite périodique.
1.7 Application de Gauss et fractions continues en théorie des nombres.
1.8 Moulin de Lorenz et flot de Lorenz
1.9 Billard de sinai avec obstacles ronds ou obstacles rectangles
Observer la diffusion abnormal, cf results de pascal Hubert
1.11 « Dripping faucet » , « le robinet qui goutte »
Faire expérience, et étudier le modèle. Référence: the Dripping faucet as a model chaotic system, shaw.
Modèle simple: un oscillateur vertical
, raideur 1, masse qui augmente
et si
on impose
et
(ou autre d’après conservation impulsion et énergie).
Réaliser une expérience avec rendu sonore (la goutte tombe sur une plaque de métal).
1.14 Théorie Astronomique du climat
A partir de données issues de carottes de glace de l’antartique, on reconstruit la teneur en CO2 de l’atmosphère sur les 400000 ans précédents, ainsi que la température moyenne de la Terre. A l’aide d’un programme, en fait une décomposition en fréquences de ces données, afin d’identifier des cycles dominants, leur amplitude et période. Par ailleurs des calculs astronomiques nous donnent les paramètres de l’orbite terrestre sur cette même durée. On déduira l’ensoleillement moyen de la Terre, et on essayera de reconnaitre les même cycles.
1.15 (4) Relativité générale:
Objectif:
étude du mouvement d’une planète ou d’un rayon lumineux autour d’une étoile ou d’un trou noir. Dans le cadre de la relativité générale, l’étoile déforme l’espace temps autour d’elle. Une planète ou un rayon, lumineux ne subit pas directement une force (gravitation) de la part de l’étoile, mais avance ”le plus droit possible” dans cet espace courbe. On dit que sa trajectoire est une géodésique. Dans un premier temps: étude des géodésiques sur une surface (2D) courbe; Puis étude de la trajectoire d’une planète. Observation du mouvement elliptique (Newton), mouvement du périhélie, inflexion des rayons lumineux, décalage gravitationel de la lumière vers le rouge, effet du trou noir. Pour les indications physiques, voirle cours de
mécanique analytique et surtout
l'exercice 2 et sa solution du TD6 de mécanique analytique (L3). Voici des
notes manuscriptes qui indiquent comment programmer et tester ce projet.
1.16 (3) Simulateur de Vol d’un planeur
Intégration numérique du vol d’un planeur en temps réel dans le plan vertical; On agit sur la position du manche à balet en temps reel.
Enonce et
solutions. Visualisation du paysage en
trois dimensions et
solutions.
1.17 Simulateur de Vol d’une montgolfiere
Idem pour une montgolfière. Le pilote agit sur les gaz. Un vent d’est derive l’engin vers des montagnes.
1.18 (2) Pendule chaotique
1.19 Dynamique du pendule pulsé
L’objet suivant (les 3 disques bleus) tourne librement autour de son axe. Mais un aimant extérieur (rouge) lui impose une force périodique. Il en résulte un comportement chaotique.
On fera des images stroboscopiques des trajectoires, afin d’observer la transition de l’ordre vers le chaos. Observation du phénomènes des résonnances, et accrochage de fréquences.
1.20 Boule du jeu Billard
étude assez complete du mouvement de la boule dans jeu de billard, soumise au chocs, aux frottements,...
1.21 Etude dynamique d’un billard parfait
On étudie les trajectoires d’une particules dans un billard au bord quelconque, sans frottement et se réfléchissant parfaitement sur les bords. Observation du chaos pour certaines formes du billard.
1.22 Modèle des anneaux de Saturne
Modèle dynamique de poussières tournant autour de Saturne, faisant apparaitre les gaps par effets de résonances.
1.23 Optique géométrique
Comportement d’un rayon lumineux à travers un système optique. Comportement à travers une interface, sur un miroir, aberration d’une lentille. Effets non linéaires et chaos.
1.24 Oscillation de deux masses couplées
Phénomène de battement. Section de Poincaré et chaos.
1.25 Bouteille magnétique
- Calcul trajectoire particule dans champ magnétique donné - champ uniforme, puis dipôle (incliné ou non)
- Effet "bouteille magnetique" (?), vitesse de dérive - éventuellement, calcul du champ electromagnétique rayonné par les particules.
1.26 Astronomie
- Voyage d’une sonde spatiale à travers le système solaire. Observer le phénomène de ”billard planétaire”
- Mouvement du soleil et des planètes dans le ciel, éphémérides.
- (*) Mouvement relatif de la Terre et du Soleil puis Construction d’un cadran solaire. Equation du temps. Documents
- Couplages entre les effets de marée et rotation d’une planète.
- Prévision des éclipses de Lune et de Soleil.
3 EDP et problèmes d’ondes
3.1 (5) Chaos Quantique
(Avec le modèle du Harper Pulsé, décrit
ici).
La mécanique quantique stipule que les propriétés spatiales d’un objet sont décrites par une onde, appelées onde quantique. L’équation de Schrödinger gouverne l’évolution des cette onde au cours du temps.
Il y a des cas où cette onde peut être concentrée spatialement, et former un paquet d’onde. On observe alors que ce paquet d’onde se déplace au cours du temps et reste localisé pendant un certain temps. Ce paquet d’onde peut avoir une taille microscopique, et à grande échelle, il apparait comme un point, une particule. Il est donc de considérer ce paquet d’onde comme une particule dite “classique”, qui correspond au centre du paquet d’onde. Le mouvement de cette particule classique est décrite par une équation de mouvement classique de type Newton (m a=F) ou Hamilton (voir ci-dessous). C’est ainsi que l’on a une correspondance entre la mécanique quantique et la mécanique classique.
Le but de ce projet est de simuler dans un modèle simple le mouvement de ces paquets d’ondes quantiques, et simultanément d’observer le mouvement de la particule classique, afin de comprendre cette correspondance.
Il apparait que certaines équations de mouvement classiques sont simples en apparence, et déterministe, mais engendre un mouvement très complexe de la particule, imprévisible, dit chaotique. On se pose alors la question de savoir ce que devient le paquet d’onde dans ce cas? C’est le problème du chaos quantique. D’un point de vue plus mathématique, cette situation est paradoxale à priori car l’équation de Schrodinger est linéaire et ne peut donc engendrer du chaos en principe. Le modèle qui suit permettra d’observer cette situation.
La première étape sera d’observer le comportement classique de la particule, et ensuite le comportement quantique.
3.2 (*2 puis 3) Formation de nuages sur une montagne par condensation de l’humidité
Ecoulement d’une atmosphère humide sur un massif montagneux, et formation de nuages par condensation de l’humidité.
Enoncé
3.3 Son: réalisation d’un accordeur pour intrument musical
On détecte une note. On mesure sa fréquence et on indique s’il est écartée de l’échelle du tempérament égal où le La (A) est à
.
Enoncé.
3.4 Son: détection d’une note musicale
- détection d’une note musicale Enoncé.
- détection du rythme
- écriture automatique de la partition avec le logiciel lilypond.
3.5 Son: détection d’une polyphonie
C’est à dire détecter la présence de plusieurs notes simultanées et identifier leur timbre.
3.6 Evolution d’une onde dans une cavité bidimensionelle (Billard).
Pour cela on calcule d’abord les modes propres de la cavité grace à la méthode de superposition d’ondes planes, et de la "décomposition QR "d’une matrice. Ref: "Betcke-Trefethen Reviving the method of particular solutions 2005". Ensuite, par surperposition d’ondes stationnaires on construit un paquet d’onde que l’on fait évoluer dans la cavité. On pourra comparer l’évolution du paquet d’onde à la dynamique d’un rayon dans la même cavité.
3.7 (4) Etats stationnaires et effet tunnel en mécanique quantique
Etude d’un potentiel confinant (polynomial). Résolution de l’équation stationnaire dans la base de l’oscillateur
Harmonique.
Calcul et dessin des etats stationnaires. Evolution d’un paquet d’onde atomique dans une molécule; comparaison
avec la trajectoire classique; étalement du paquet, effet tunnel...
3.8 Solitons dans une chaine d’atomes
Que se passe t-il dans un système où les atomes (sur une ligne)
interagissent avec leur voisins avec des forces non harmoniques?
Force harmonique=-k.r
Force non harmonique=a*(exp(-b.r)-1)
où r est le déplacement de la position d’équilibre de la
distance entre l’atome et son voisin. On représentera la propagation
d’une déformation. Dans le cas non harmonique apparaissent des ondes
d’un caractère spécial qui se déplacent sans dispersion: des
solitons. Mots clefs: solitons, solitons de Toda, (theory of
non-linear lattices springer Verlag, 1989).
3.9 Ondes solitaires à la surface de l’eau
Modele de KdV . Resolution par Runge-Kutta.
3.10 Transformée par ondelettes d’un signal sonore
Representation du signal dans le plan temps-frequence, avec les ondelettes de Morlet.
Applications à la musique: reconnaissance des notes, du timbre
d’instruments, d’accords...
3.11 Microscope en ondelettes d’une image (méthode de Daubuchie)
on veut représenter une
image en s’affranchissant plus ou moins des détails de celle-ci. On peut
arriver à quantifier le plus ou moins grâce à une double
transformée rapide en ondelettes. Applications à la compression
d’images.
3.12 (*3) Vibration d’une membrane de tambour de forme quelconque
3.13 Evolution d’un paquet d’onde en milieu dispersif
On résoud l’équation dépendant du temps par différences
finies.
On observe le déplacement du paquet d’onde son étalement
(dispersion),
et les phénomènes de réflexion et transmission sur une barrière de potentiel.
3.14 Localisation des ondes dans un milieu désordonné
- On prend par l’exemple d’ondes quantiques d’électrons dans un milieu désordonné (métal contenant des impuretés) se traduisant par un potentiel V(x) aléatoire. On étudie la forme des ondes stationnaires, et la propagation des paquets d’ondes. On observe le résultatremarquable (et curieux) qu’à une dimension ces ondes sont toujours localisées dans une région de l’espace : ”localisation de Anderson”.
3.15 Equation de la chaleur
- jouer avec conditions initiales et aux limites, permettant d’introduire d’autres termes dans l’équation (écoulement)
3.16 Détection et reconnaissance d’une vibration
Supposons un micro posé sur la table. On tape légèrement avec le doigt à un endroit A de la table. Cela crée une vibration sonore dans la table qui est détectée par le micro. L’ordinateur enregistre ce signal noté
. Si on tape à un autre endroit B, l’ordinateur enregistre un autre signal
. Dans la suite si on tape en A ou B l’ordinateur enregistre un signal
et en faisant un produit scalaire
et
l’ordinateur pourra savoir quel pooint on a touché.
A partir de cette idée, on peut faire imaginer des lettres A,B,C.. sur la table à des endroits précis, dans un premier temps l’ordinateur enregistre le signal obtenu si on tape sur chacun de ces endroits (dans l’ordre) et dans un deuxième temps, on tape un texte avec ces lettres imaginaire et le programme reconnait les signal donc les lettres et écrit le texte sur l’écran.
Etapes du projet
- Partir de l’exemple donné sur l’utilisation du micro. Observer qu’il affiche un signal bruité en permanence et que si on tape sur la table, ce signal devient plus intense. Modifier le programme pour qu’il extrait un signal sur une durée
si le signal observé dépasse un certain seuil.
4 Physique statistique et thermodynamique
4.1 (*4) Mouvement de particules dans une enceinte
Simulation du mouvement et des chocs de particules dans une enceinte;
Cela permet l’observation de propriétés statistiques: apparition du désordre,
(entropie, loi de Bolztmann), équilibre thermodynamique, loi des gaz
parfait PV=nRT, transformation adiabatiques,...
4.2 (3) Mécanismes de l’aimantation et ferromagnétisme
(Voir Cours de M1: système dynamiques, chapitre sur les dynamiques de Markov ).
On considère un réseau
périodique dont les sites sont notés
et dont les variables sont
et modélisent des spins (up/down). Si deux sites sont voisins on note
. Une configuration des spins est un champ donné:
. Son énergie ferromagnétique est:
- Combien y a t-il de sites? Quelles configurations donnent l’énergie minimale? et maximale?
- L’algorithme d’évolution suivant est appelé algorithme de montecarlo.
est un paramètre fixé qui est l’inverse de la température:
.
- A l’instant
, on part d’une configuration
choisie au hasard.
- On choisit un site
au hasard, et on note
la configuration identique à
sauf au site
où le spin est opposé:
(spin opposé).
- A l’instant
,
- Si
on choisit la configuration
.
- Si
on choisit la configuration
avec la probabilité
(et on reste donc avec
avec la probabilité complémentaire).
- On revient en (b) pour poursuivre l’évolution du champ de spins.
Cet algorithme correspond à une matrice stochastique réversible pour la mesure de Boltzmann:

- Ecrire un programme qui fait évoluer et représente la configuration du champ de spins
.
- L’aimantation de la configuration
est
. Représenter l’aimantation au cours du temps. Quand l’équilibre statistique est atteind, calculer la moyenne
et la variance
. Représenter
et
en fonction de
.