Examples simples de calculs Saint Venant 1D

PYL
à jour 08/03/14
en liaison avec le cours de M1 "Mécanique des Fluides en Milieu Naturel"
cours du Master sciences et technologies UPMC Mention Sciences de l'Ingénieur (SDI),



Code explicite simple qui ne marche pas très bien
FACULTATIF
pour référence la méthode différences finies centrées qui marche à peu près
fichier différences finies simples en C, onde1.c, onde2.c et onde3.c.

compilation

cc -std=C99 -lm onde2.c

lancer :
./a.out | gnuplot

l'affichage se fait par gnuplot, voir le cours...


Rupture de Barrage
OBLIGATOIRE

Programme en C simple fichier (en cours on a regardé un fichier plus compliqué). Cas de rupture de barrage flux Rusanov (essayer HLC) compilation

on le compile et on crée l’exécutable db
cc -O3 -ffast-math -std=C99 -lm svdb.c -o db

pour lancer le programme:
 ./db

un fichier appelé solxhQt.OUT avec x h Q et t est créé, pour le tracer, on lance gnuplot
 set ylabel "t"; set xlabel "x"; set hidden3d;
 sp’solxhQt.OUT’ u 1:4:2 not w l linec 3
On obtient la courbe suivante hauteur d’eau en fonction de x et de t.








le dessin animé de droite est obtenu en enlevant le commentaire /* */ à la fin
./db | gnuplot

 pour la solution exacte voir le cours...




Ressaut hydrodynamique fixe
FACULTATIF

Programme en C simple fichier.
Exemple du cours, modifier les condiftions pour le faire avancer...

cc -lm sv_ressaut.c ;./a.out | gnuplot







Arrivée d'un Tsunami sur la côte
FACULTATIF
Une perturbation initiale se casse en deux vagues dont une déferle sur la côte
Programme en C simple fichier.

cc -lm sv_plage.c ;./a.out | gnuplot









Ecoulement sur un obstacle, influence du Froude
FACULTATIF

Programme en C simple fichier.
Ecoulement transcritique sur une bosse; le flot est accéléré, passe Fr=1 au sommet, un choc se forme derrière la bosse.

 cc -lm sv_transc.c ;./a.out | gnuplot







FACULTATIF
Programme en C avec reconstruction hydrostatique fichier zippé. compilation
make
exécution de plusiseurs cas les uns après les autres
./run_sv.sh

lire le fichier de la version en cours de la partie Résolution numérique de Saint Venant dans les fluides et les granulaires: fichier pdf.