BIGFLOW - Binôme E
EXEMPLES D'UTILISATION |
1. Introduction
2. TEST1 en régime stationnaire
3. TEST1 en régime transitoire
4. TEST6
Description du test
Il s'agit d'un écoulement dans un milieu poreux saturé et isotrope.
Pour effectuer ce test, on considère un cube. L'écoulement traverse le cube par l'une des faces et ressort par la face opposée (axe des x). Toutes les autres faces sont considérées comme imperméables.
On impose donc sur la première face une condition de Dirichlet (H=1) ainsi que sur la face opposée (H=0). Sur toutes les autres faces , on impose une condition de Neumann ( flux constant) : le flux est nul.
Comment remplir INPUT1
Pour effectuer ce test simple, seul le fichier INPUT1 est nécessaire.
Nous avons effectué ce test pour Nx=101 et Ny=Nz=5.
Le fichier INPUT1 est le suivant :
Résultats du test
On obtient l'évolution suivante de H sur
la ligne centrale du cube :
On obtient bien une évolution linéaire de
la hauteur piézométrique H.
Description du test
Il s'agit du même test que celui décrit précédemment sauf que le régime n'est plus permanent mais transitoire.
Comment remplir INPUT1
Le fichier INPUT 1 est le suivant :
Résultats du test
L'évolution de H sur la ligne centrale au cours du
temps est la suivante :
On retrouve bien une solution transitoire de forme exponentielle
qui tend vers la solution stationnaire quand le temps augmente.
Retour au Sommaire
Description du test
Il s'agit d'un geyser : l'écoulement provient de deux directions opposées (axe des X), et s'écoule suivant l'axe vertical par un trou à la surface du cube.
Sur les deux faces opposées de l'axe des x, on impose une condition de Neumann : sur l'une le flux est fixé à 1, sur l'autre à -1. Sur la face supérieure, on a une condition mixte :Neumann à l'extérieur du trou (flux nul) et Dirichlet à l'intérieur du trou (H imposée). Sur les autres faces, le flux est nul.
Comment remplir INPUT1, 2 et 3
Le fichier INPUT1 pour ce test est le suivant :
Fichier INPUT2
Il faut entrer :
Fichier INPUT3
Il faut entrer les valeurs suivantes :
Résultats du test
On visualise deux types de graphe :
Dans les plans (x,z), (y,z) et (x,y), on obtient les iso-H
suivantes :
Dans les plans (x,z), (y,z) et (x,y) on obtient les évolutions
de H suivantes :
On retrouve la symétrie de l'écoulement par
rapport aux axes x et y. On observe bien une chute de H vers le trou sur
la face supérieure de la boîte. Plus on s'approche du trou,
plus les iso-H sont resserrées .