Opérateur ENCEINTE DATE 99/12/09 Procedure ENCEINTE __________________ SYNTAXE : __________ ENCEINTE NDT RXT ; OBJET : ----- La procedure ENCEINTE Calcule une injection de vapeur dans une enceinte fermee et/ou l'aspersion par gouttes d'eau: Commentaires ____________ NDT ENTIER : Nombre de pas de temps RXT TABLE : contenant les informations suivantes La constitution d'un jeu de données revient à remplir la table RXT et ensuite à appeler la procédure ENCEINTE avec cette table RXT. Les indices et leur signification sont les suivants : rxt = TABLE ; Définition de la table de travail (OBLIGATOIRE) * Définition des objets géométriques rxt.'vtf' = MAILLAGE ; Définition du maillage fluide (OBLIGATOIRE) rxt.'breche' = MAILLAGE ; Définition de la localisation de la brèche (OBLIGATOIRE pour le moment). rxt.'pi' = POIN ; Définition d'un point intérieur au domaine fluide ce qui sera utilisé pour imposer la pression en ce point (OBLIGATOIRE). * Définition des murs pour la thermique paroi rxt.'MATMUR' = TABLE ; Définition d'une table pour les matériaux dans le cas de murs traités en thermique paroi rxt.'MATMUR'.'NMAT' = ENTIER ; Nombre de matériaux différents rxt.'MATMUR'.'LAMBDA' = TABLE 'VECTEUR' ; Table des conductivités thermiques (W/m/K)pour chaque matériaux 'NMAT' valeurs rxt.'MATMUR'.'LAMBDA' . 1 = FLOTTANT ; rxt.'MATMUR'.'LAMBDA' . 2 = FLOTTANT ;... rxt.'MATMUR'.'ROCP' = TABLE 'VECTEUR' ; Table des $\rho*Cp$ de chacun des matériaux (J/m3/K) rxt.'MATMUR'.'ROCP' . 1 = FLOTTANT ; rxt.'MATMUR'.'ROCP' . 2 = FLOTTANT ; rxt.'TABMUR' = TABLE ; Définition de la table pour les murs du cas test considéré. Les murs peuvent être de deux catégories différentes : les murs dont on voudra connaître l'évolution thermique au cours du temps (type THERMP) et les murs qui ne servent que de conditions aux limites pour l'équation en température du fluide Tf (type CLIMT). rxt.'TABMUR'.'NMUR' = ENTIER ; Nombre de murs rxt.'TABMUR'.'TYPE' = TABLE 'VECTEUR' ; Table définissant le type des murs rxt.'TABMUR'.'TYPE' . 1 = 'THERM' ; Mur dont on voudra traiter l'évolution thermique interne en fonction du temps, il faudra dans ce cas spécifier le matériau et les conditions initiales. rxt.'TABMUR'.'TYPE' . 2 = 'CLIMT' ; Mur qui servira de conditions aux limites pour l'équation d'évolution de la température du fluide. Ce mur sera traité comme une condition de flux imposé (opérateur FIMP). rxt.'TABMUR'.'MAILLAGE' = TABLE 'VECTEUR' ; Table contenant les maillages associés aux murs rxt.'TABMUR'.'MAILLAGE' . 1 = MAILLAGE ; Définition des maillages rxt.'TABMUR'.'MAILLAGE' . 2 = MAILLAGE ;... rxt.'TABMUR'.'MATER' = TABLE 'VECTEUR' ; Table definissant les matériaux utilisés pour les murs de type THERMP rxt.'TABMUR'.'MATER' . 1 = ENTIER ; Numéro du matériau défini dans la table MATMUR rxt.'TABMUR'.'TPINI' = TABLE 'VECTEUR' ; Table définissant les conditions initiales en température pour les murs. Les murs de type THERMP ont une température initiale constante (FLOTTANT). Pour les murs de type CLIMT, on associe le nom d'une table de la table scenario. rxt.'TABMUR'.'TPINI' . 1 = FLOTTANT ; rxt.'TABMUR'.'TPINI' . 2 = MOT ; * Définition des logiques pour spécifier les options rxt.'ASPERS' = LOGIQUE ; Variable logique définissant la présence d'aspersion dans le cas test. rxt.'VAPEUR' = LOGIQUE ; Variable logique définissant la présence de vapeur dans le cas test. rxt.'GRAPH' = LOGIQUE ; Variable logique permettant l'obtention de graphes 2D sur les variables OD et les variables moyennes. rxt.'PERMANENT'= LOGIQUE ; Variable logique permettant la reconnaissance d'un cas test évoluant vers un régime permanent. Dans ce cas, on peut redéfinir le critère de convergence rxt.CRIT, par défaut il est imposé à $10^{-4}$. rxt.'FRPREC'= ENTIER ; fréquence de recalcul du préconditionnement imposée par l'utilisateur exprimée en terme de nombre de pas de temps. Si on veut que la procédure calcule elle-même la fréquence de repréconditionnement, il ne faut pas spécifier cette valeur et définir l'indice suivant. rxt.'FACPREC'= ENTIER ; écart minimum en terme d'itérations nécessitant le recalcul du préconditionnement. rxt.'CORVIT'= LOGIQUE ; Variable logique permettant de ne pas localiser le correction 0D dans des zones du maillage ou la vitesse est faible. rxt.'CORLOC'= LOGIQUE ; Variable logique permettant de localiser la correction 0D dans des zones au prorata de la présence de la variable à corriger. rxt.'MSCONS'= LOGIQUE ; Variable logique permettant de traiter les sources de masse et d'énergie en variables conservatives. rxt.'RESTART'= LOGIQUE ; Variable logique permettant de forcer le recalcul du préconditionnement lors d'une reprise afin d'éviter les problèmes de renumérotation. * Définition du scénario de calcul - Conditions initiales dans le volume fluide rxt.'scenario' = TABLE ; Table définissant le scénario rxt.'scenario'.'TF0' = FLOTTANT ; Température initiale de fluide en °C. rxt.'scenario'.'PT0' = FLOTTANT ; Pression totale initiale en Pascal. rxt.'scenario'.'Yvap0' = FLOTTANT ; Fraction massique initiale de vapeur d'eau dans le volume fluide rxt.'scenario'.MOT = FLOTTANT ; Fraction massique des incondensables initialement présents dans le volume fluide. MOT peut être égal à 'Yhe0', 'Yh20', 'Yo20', 'Yn20', 'Yco0', 'Yco20'. Il est à remarquer que l'air est considéré dans ce modèle comme un composant à part entière 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