- le programme
(pour obtenir les simulations correspondant aux exemples suivants, penser à modifier le temps maximum de simulation, la donnée initiale, les conditions aux limites, et le champ de vitesse dans le programme)
- exemple 1 : convection dans un canal avec champ de vitesse
parabolique
- exemple 2 : convection dans un carré avec champ de
vitesse en "rotation rigide"
- le programme
(pour obtenir les simulations correspondant aux exemples suivants, penser à modifier le temps maximum de simulation, la donnée initiale, les conditions aux limites, et le champ de vitesse dans le programme)
- exemple 1 : convection entre deux cylindres avec champ de vitesse
en "rotation rigide" (comme exemple 2, partie 1)
- exemple 2 : convection entre deux cylindres avec écoulement de
Couette cylindrique
- exemple 3 : convection autour d'un cylindre, champ de vitesse
à potentiel complexe
- exemple 4 : même exemple, calcul sur la moitié du domaine
seulement (avec conditions de symétrie)
- exemple 5 : comme l'exemple 3, avec donnée initiale nulle, et
entrée de traceur sur une partie du bord gauche
- maillages générés avec : EMC2 (logiciel libre, avec documentation, code source et versions pré-compilées). Pour chaque exemple, le fichier décrivant la géométrie est donné (xxx.emc2_bd). On peut le charger dans emc2 pour modifier les paramètres de raffinement du maillage
- pour vous aider à générer un maillage, voici une marche à suivre pour l'exemple 3
- le programme
- exemple 1 : convection entre deux cylindres avec champ de vitesse
en "rotation rigide" (comme exemple 2, partie 1)
- exemple 2 : convection autour d'un cylindre, champ de vitesse
à potentiel complexe
- les données
- le fichier de maillage
- resultat de simulation : maillage, champ de pression
stationnaire (où l'on distingue une pression plus importante
dans le réservoir de droite), évolution des champs de
pression et de vitesse dans le tube de jonction (zoom). Note : la
différence de pression générée par le dispositif simulé ici est
très faible : moins de 0.2%
- le programme
- exemple 1 : diffraction d'un choc plan sur un cylindre
(ci-dessous isochores à gauche (lignes d'égales densité) et
superposition lignes de même nombre de Mach et champ de
vitesse). Le nombre de Mach est calculé dans visit en sélectionnant
dans le menu : Commandes/Expressions
- le programme
- exemple 1 : écoulement autour du nez d'une navette spatiale
(ci-dessous isochores, isobares, isothermes, isomach).
- le programme
- exemple 1 : ondes générées par une "goutte" tombant dans un
bassin, fond plat.
- exemple 2 : lac au repos, fond avec butte, le schéma ne
préserve pas cet état et génére des
vagues artificielles. Surface libre (en bleu) et fond (en vert).
- exemple 3 : rupture d'un barrage dans un canal, fond plat. Surface libre (en bleu) et fond (en vert).
- exemple 4 : rupture d'un barrage dans un canal, fond avec
butte. Là encore, les vagues au dessus de la butte avant l'arrivée
de l'eau du barrage sont artificielles, et dues à la non
préservation de l'état de lac au repos. Surface libre (en bleu) et fond (en vert).
- exemple 5 : rupture d'un barrage dans un canal, avec "île" au
milieu du canal. L'arrivée de l'eau du barrage sur l'île provoque
une vague qui remonte le courant. Surface libre (en bleu) et fond (en vert).
- le programme (en matlab !!)
- exemple 1 : instabilités de Kelvin-Helmoltz (tiré du livre
"Introduction au calcul scientifique par la pratique", de Danaila et
al., Dunod), marqueur passif et vorticité.
