Le programme de simulation.



Le problème numérique.

PHOENICS.

Le fichier Q1.

Le fichier GROUND.

PHOTON.

Script d'automatisation.

Traitement sur XMGR.



Sommaire du compte rendu.





LE PROBLEME NUMERIQUE.

Les équations à résoudre sont celles de conservation de la masse, de la quantité de mouvement et de l'énergie en régime laminaire.
Pour exprimer rho de manière simple, on l'exprime en fonction de T, grâce à l'
approximation de Boussinesq. Ceci revient à négliger les variations de rho, sauf dans le terme de flottabilité.

Sommaire.


PHOENICS.

Parabolic Hyperbolic Or Elliptic Numerical Integration Code Series.

PHOENICS est un code de calcul pour la simulattion d'écoulement de fluide et de transfert de chaleur.

Il résoud en particulier des problèmes:

STRUCTURE DE PHOENICS:

Sommaire.


Le fichier Q1.

C'est le fichier de données de PHOENICS.

Les conditions initiales et aux limites:

Le maillage:

Nous avons testé 2 sortes de maillages afin d'en voir succintement les effets (l'influence du maillage a été étudiée plus en détails par le binôme Naud-Souleres):

Pour le créer, nous avons défini deux zones en X de largeurs égales et avec le même maillage. En fait, chaque zone en X et la zone en Y possèdent le même genre de maillage (maillage plus rapprochés aux bords).

Sommaire.


Le fichier GROUND.

Nous avons rajouté dans le fichier ground.f:


Sommaire.


PHOTON.

C'est le post processeur graphique de PHOENICS, qui nous permet de visualiser les résultats

maillage

vecteurs vitesses

évolution de T, V, ...

Sommaire.


Script d'automatisation.

La recherche du Rayleigh critique nous obligeant a relancer plusieurs fois le même calcul en changeant juste la valeur de HBAS, nous avons décidés d'automatiser cette tâche . Nous avons donc écrit un script en shell avec le binôme Béon - Galateau. Ce script est donc écrit en Korn-Shell et fait appel à des subroutines écrites avec les utilitaires awk et sed. Il a les fonctionnalités suivantes :


Modification de la température dans Q1;
Lancement de runsat et earexe avec sauvegarde du nombre de    Rayleigh dans un fichier;
Extraction des paramétres de calcul du nombre de Nusselt et de    Ux_Max dans les fichiers RESULTet PHI;
Calcul du nombre de Nusselt et création d'un fichier comportant les nombres de Nusselt et de Rayleigh.

Toutes ces opérations sont incluses dans une boucle permettant l'incrémentation de la température.
Il ne reste ensuite plus qu'à exploiter les résultats sous XMGR.

Sommaire.


Traitement sur XMGR.

La recherche du Rayleigh critique implique le traitement de quelques donnés telles que tracé le nombre de Nusselt en fonction du nombre de Rayleigh. Nous avons tracé ces courbes grâce au logiciel XMGR

Sommaire.