DATE     99/12/09
                                                                                
    Procedure ENCEINTE
    __________________


    SYNTAXE  :
    __________


    ENCEINTE NDT RXT ;

 et il faudra faire attention lors de la poursuite par un
 éventuel calcul de combustion.

- Conditions aux limites pour l'injection et les parois
  pilotées thermiquement

 rxt.'scenario'.'t'      = LISTREEL ; Liste des temps en secondes.
 rxt.'scenario'.'qeau'   = LISTREEL ; Liste des débits de
                                      vapeur injectés en kg/s.

 rxt.'scenario'.MOT   = LISTREEL ; Liste des débits d'incondensables
               injectés en kg/s. MOT peut être égal à 'qair', 'qhe',
               'qh2', 'qo2', 'qn2', 'qco', 'qco2'.

 rxt.'scenario'.'tinj'   = LISTREEL ; Liste des températures
                          d'injection en fonction du temps.

 rxt.'scenario'.MOT   = LISTREEL ; Liste de la température
 paroi (CLIMT) en fonction du temps. La valeur de MOT doit
 être cohérente avec le nom donné dans la table TABMUR à
 l'indice TPINI.

Cas avec ASPERSION

 En cas d'aspersion le logique rxt.'VAPEUR' doit etre = VRAI
 En cas d'aspersion (rxt.'ASPERS' = VRAI) rajouter:

 rxt.'aspinj' = MAILLAGE ; Maillage de la surface d'injection
                           de la phase dispersee
 rxt.'toitf'  = MAILLAGE ; Maillage de la surface superieur (par
                           rapport a la verticale) du volume
                           fluide
 rxt.'rod'    = MAILLAGE ; Densite de la phase dispersee
 rxt.'Cpd'    = MAILLAGE ; Chaleur specifique de la phase dispersee
 rxt.'scenario'.'vzinj'  = LISTREEL ; Liste des vitesses (verticale)
                                      d'injection de la phase disp.
 rxt.'scenario'.'xdinj'  = LISTREEL ; Liste de la fraction volumique
                                      de la phase disp. a l'inject.
 rxt.'scenario'.'tdinj'  = LISTREEL ; Liste de la temperature de la
                                      phase disp. a l'injection.
 rxt.'scenario'.'tdinj'  = LISTREEL ; Liste du diametre de la phase
                                      disp. a l'injection.

* Définitions concernant l'algorithme de résolution

 rxt.'discr'  = MOT ; Choix de la discrétisation, les choix
possibles sont 'MACRO', 'LINE'et 'QUAF'. Par défaut,
l'option MACRO est utilisée.
 rxt.'ALGO'   = MOT ; Choix de l'algorithme de résolution du
système Pression/Vitesse. MOT peut être 'IMPL'dans le cas
d'un algorithme implicite ou 'SPLT'dans le cas d'un
algorithme de projection. Par défaut, on choisit
l'algorithme de projection.
 rxt.'KPRESS' = MOT ; Choix pour la discrétisation de la
pression. MOT doit être égal à ``CENTRE'pour l'algorithme de
projection et ``CENTREP1'pour l'algorithme implicite.
 rxt.'STAB' = FLOTTANT ; Ce réel correspond au facteur de
stabilisation des éléments MACRO d'ancienne génération.
 rxt.'nimpr'= ENTIER ; Indice gérant les impressions
intermédiaires en cours de résolution du système. Si
nimpr=0, il n'y a pas d'impression (Option par Défaut). Si
nimpr=1, on imprime le nombre d'itérations et les maximums
des résidus.

* Coefficients de contrôle du modèle

 rxt.'epsi'   =  FLOTTANT ; précision des maillages lors de
 l'utilisation de l'opérateur ELIM.
 rxt.'DT0'    =  FLOTTANT ; valeur du pas de temps, si on ne
 spécifie pas cette donnée, le pas de temps est calculé par
 la procédure à partir du maximum de la vitesse et d'une
 dimension caractéristique égale à la racine cubique du volume.
 rxt.'FROG'= FLOTTANT ; facteur correctif  sur le
 terme de flottabilité dans l'équation de quantité de mouvement.
 rxt.'FROC'= FLOTTANT ; facteur correctif sur le débit de
 condensation, ce facteur englobe les manques de la
 modélisation actuelle (Vitesse nulle en paroi) et les
 simplications du modèle (thermique paroi ne comportant pas
 suffisamment de mailles).
 rxt.'alfg'= FLOTTANT ; coefficient de relaxation sur le
 terme source de flottabilité dans l'équation de la quantité
 de mouvement.
 rxt.'alfc'= FLOTTANT ; coefficient de relaxation sur le
 débit de condensation dans l'équation de transport de la
 densité de vapeur d'eau.
 rxt.'FRPREC'= ENTIER ; Fréquence de recalcul de
 préconditionnement. La valeur est donnée en terme de nombre
 de pas de temps pour l'instant.

* Propriétés physiques du mélange gazeux

 rxt.'muair'  =  1.800e-5 ;
 rxt.'alf'    =  1.800e-5 ;
 rxt.'lamair' =  2.580e-2 ;
 rxt.'db'     =  1.000e-5 ;
 rxt.'KH'     =  1.000e+1 ;
 rxt.'MUT'    =  1.000e-2 ;
 rxt.'Cpair' =  1000.0 ;
 rxt.'Cph2o' =  1700.0 ;
 rxt.'Cphe'  = 14550.0 ;
 rxt.'Lv'    =     2.3E6 ;
 rxt.'Rgair' = FLOTTANT ;
 rxt.'Rgh2o' = FLOTTANT ;
 rxt.'Rghe'  = FLOTTANT ;

* Instruction pour lancer la procédure ENCEINTE

 ENCEINTE NPAS RXT ; NPAS nombre de pas de temps à réaliser
 et RXT table contenant les informations nécessaires à la
 résolution du problème (cf. ci-dessus).


 LIMITATIONS DU MODELE ACTUEL

* modèle gaz parfait
* modèle de turbulence grossier (adherence a la paroi , Nut = 500. * Nu)
* Pas de possibilite de plan de symetrie

                                                                                
                                                                                
     Description :
     _____________

 The ENCEINTE procedure computes an injection of steam into a closed
 containment :

     Contents :
     __________

    NDT   ENTIER : Number of time steps
    RXT   TABLE  : contening the following informations

    for more details see the french version