On veut d’abord résoudre le problème d’advection scalaire par un schéma de transport semi-Lagrangien.
Ecrire ce schéma dans un domaine périodique en completant le programme :
Quelle observation peut-on faire sur le transport de colorant ?
Quel critère doit satisfaire dt pour que le schéma soit stable ?
Modifier le programme
Ajouter une visualisation de l’obstacle en utilisant la commande ``contour’’ de Matplotlib.
En utilisant la commande ``streamplot’’ de Matplotlib, ajouter une visualisation des lignes de courant de l’écoulement à la fin de la simulation.