Lecture du maillage

Actuellement, PRISSMA est capable de lire deux types de maillage : le format dit neutral du logiciel Gambit (le mailleur associé à Fluent) et le format de maillage AVBP ^3.1. Le code sait lire des maillages hybrides, mais les maillages non coïncidents tels que des maillages hexaèdres/tétraèdres ne sont pas supportés.

Pour lire les données du maillage, on utilise la routine predatas. Le script commence par vérifier la présence du fichier de configuration predatas.choices dans le répertoire de travail. Au cas où il ne le serait pas, le script en crée un modèle dont il faudra modifier le contenu et relancer predatas. Ce fichier est au format texte et se présente sous la forme suivante :

'./meshes/cyl_coarse.neu'  ! Mesh file name
'./QUADRA'                 ! Quadrature definition directory
2                          ! Mesh type (1=AVBP, 2=GAMBIT)
'./INFILES/'               ! Output path
SNDOM                      ! Quadrature type [SNDOM|FVM|TNDOM|PERSO]
4                          ! Quadrature order
2                          ! Number of directions if PERSO quadrature
1. 0. 0.                   ! 1st direction for PERSO quadrature
-1. 0. 0.                  ! 2nd direction ...

Le fichier permet de définir :

ligne 1

le chemin vers le fichier de maillage ;

ligne 2

le chemin vers le répertoire contenant les données sur les modèles de quadrature ;

ligne 3

le type de maillage, on utilise 1 pour AVBP et 2 pour Gambit ;

ligne 4

le chemin vers le répertoire des fichiers d'entrée ;

ligne 5

le modèle de quadrature (à compléter) ;

ligne 6

l'ordre de la méthode SN pour fixer le nombre de directions, les valeurs possibles sont 2 pour S2, 4 pour S4, 6 pour S6, 8 pour S8, 9 pour LC11, 12 pour S12 ;

ligne 7 et suivantes

si l'on choisit l'option PERSO pour le modèle de quadrature à la ligne 5, on définit à cette ligne le nombre de directions et leur vecteur (x,y,z) aux lignes suivantes.

À partir des paramètres lus dans ce fichier, le maillage est lu et un ensemble de fichiers d'entrée sont créés dans le répertoire INFILE. Ces fichiers contiennent la géométrie du maillage, les directions et leur chemin associé de parcours des cellules du maillage.

L'étape suivante consistera à lire la solution initiale issue d'un calcul de CFD.