Les autres caractéristique du problème (FIPREP)

Une fois le maillage défini, il s'agit de donner les caractéristiques du ou des fluides utilisés et de choisir la méthode de résolution. Cela se fait soit via l'interface graphique, soit par ajouts à l'intérieur du fichier essai.FDREAD. Les principales commandes sont :

• PROBLEM(...) pour définir si le problème est LINEAR ou NONLINEAR, STEADY ou TRANSIANT, ISOTHERMAL, ENERGY ou si la température est à la source du mouvement BUOYANCY...
• TIMEINTEGRATION(NSTEP=nn, DT=xx,BACKWARD/FORWARD/TRAPEZ) pour un problème transitoire (TRANSIANT) avec nn le nombre de pas de temps, xx la taille du pas de temps et ensuite le schéma d'integration en temps choisi.
• EXECUTION(NEWJOB/RESTART) pour indiquer si on part de zéro ou si on récupère la solution du calcul précédent.
• PRESSURE(MIXTE/PENALTY,CONT/DISCONTINUOUS) pour choisir le type d'approximation de la pression.
• VISCOSITY(CONST=r) pour indiquer la viscosité du fluide.
• DENSITY(CONST=r) pour indiquer la densité du fluide.
• SOLUTION(S.S.=n/...) pour choisir la méthode de résolution

Il est toujours bon d'adimensionner son problème, ne serait-ce que pour avoir une idée sur le type d'écoulement étudié. Voici ce que donne l'écoulement autour d'un cylindre pour un Reynolds de 40 en reprenant le maillage décrit dans FIMESH (on a aussi le même problème présenté en annexe en maillant avec FI-GEN) :

FIPREP
PROBLEM(NONLINEAR)
EXECUTION(NEWJOB)
DENSITY(CONSTANT=1.)
VISCOSITY(CONSTANT=.025)
PRESSURE(PENALTY=1.E-6)
ICNODE(STOKES)
SOLUTION(Q.N.=6)
STRATEGY(S.S.=3)
DATAPRINT(NORMAL)
PRINTOUT(ALL)
ENTITY(FLUID,NAME="fluide")
RENUMBER
END

par défaut le problème est 2-D, le fait qu'il soit NONLINEAR implique que le terme de convection est présent.

la PENALTY permet d'enlever la variable <<pression>> en la liant à la vitesse par avec ici
0.0cm on initialise les algorithmes itératifs de résolution par la solution du problème de Stokes (qui est donc calculé en premier)
0cm on combine 3 itérations de Picard puis 6 itérations de quasi Newton pour résoudre finalement le problème (toutes les itérations ne seront pas faites si le calcul converge avant).
0cm RENUMBER optimise la numérotation des n uds du maillage pour réduire la largeur de bande de la matrice