Afficher cette notice en
1 : $$$$ EXECRXT NOTICE FD144363 24/03/27 21:15:03 11877 2 : DATE 24/03/27 3 : 4 : Procedure EXECRXT Voir aussi : 5 : ---------------- 6 : EXECRXT ENTI1 TAB1 ; 7 : 8 : 9 : 10 : Commentaire : 11 : _____________ 12 : 13 : ENTI1 : nombre de pas de temps 14 : 15 : TAB1 : table RXT contenant la description du probleme fluide 16 : 17 : Objet : 18 : _______ 19 : 20 : La procedure EXECRXT calcule, a partir d'un etat initial, l'evolution 21 : au cours du temps d'un melange gazeux dans une enceinte fermee. 22 : 23 :
SOMMAIRE DE LA NOTICE
---------------------
1. Generalites
2. Objets geometriques
3. Modelisation Physique
3.1 Loi d'etat pour la vapeur
3.2 Thermique parois
3.3 Temperature de paroi imposee
3.4 Conditions initiales
3.5 Conditions aux limites
3.6 Turbulence
3.7 Aspersion
3.8 Condensation en masse (en test)
3.9 PAR (recombineur catalytique passif)
3.10 Condensation en parois
4. Parametres de fonctionnement
4.1 Table RXT
4.2 Limitations du modele actuel
4.3 Quelques recettes
5. Resultats du calcul
5.1 Table rxt.'TIC'
5.2 Table rxt.'GEO'
1. Generalites
==============
24 : 25 : L'etat initial est uniforme en espace lorsqu'il est donne par 26 : l'utilisateur, issu d'un calcul precedent sinon. 27 : 28 : L'air est toujours present dans le melange gazeux. Il peut aussi contenir un 29 : ou plusieurs des constituants suivants: vapeur d'eau, H2, N2, He, O2, CO et 30 : CO2. Les gaz incondensables sont modelises par la loi des gaz parfaits. C'est 31 : aussi le choix par defaut de la vapeur (version V0), un modele gaz reel etant 32 : en test (version V1). 33 : 34 : En presence de vapeur la condensation en paroi peut apparaitre si les 35 : conditions locales sont reunies (Pvap > Psat). Le modele de condensation est 36 : de type Chilton-Colburn et associe a une correlation d'echange de type 37 : convection naturelle le long d'une plaque plane verticale. 38 : 39 : On distingue 4 types de conditions aux limites suivant la nature de la 40 : frontiere du domaine fluide : zones d'injection (breches), ventilation forcee 41 : et clapet de decharge (sorties), parois thermiques et parois inertes : 42 : 43 : - sur les zones d'injection ou breches, on impose des conditions aux limites 44 : de type valeur imposee pour la vitesse, la temperature du melange, la densite 45 : du melange et de ses constituants. Il faut donc y preciser le debit massique 46 : de chaque constituant du melange (kg/s) et la temperature d'entree (oC) (voir 47 : entree 'scenario' des Breches) ; 48 : 49 : - sur les sorties, on impose un debit (ventilation) et on precise les 50 : parametres de la perte de charge (clapet de decharge) (voir entree 'scenario' 51 : des Sorties) ; 52 : 53 : - sur les parois thermiques, la vitesse est nulle (suf si fonction de paroi 54 : (entree FPAROI)) et la temperature evolue au cours du temps (entree TIMP), est 55 : a priori constante (entree ECHANP), le resultat d'un calcul de thermique paroi 56 : (entree THERMP). Pour ce dernier cas, il est possible de coupler la resolution 57 : des equations de l'enrergie paroi et fluide de façon implicite (entree 58 : THERCO). En presence de vapeur les parois thermiques sont susceptibles de 59 : condenser. Elles sont par contre impermeables pour tous les incondensables. 60 : 61 : - sur les parois inertes, la vitesse est nulle, et elles sont impermeables 62 : pour toutes les autres inconnues (temperature, vapeur et gaz incondensables). 63 : Les parois inertes correspondent au maillage obtenu par difference entre 64 : l'enveloppe du volume fluide et les parois, breches et sorties. 65 : 66 : Par suite, les conditions aux limites sont correctement definies. 67 : 68 : La turbulence des mouvements de gaz est modelisee soit par une viscosite 69 : tourbillonnaire constante, soit par un modele de longueur de melange soit par 70 : un modele K-epsilon (entree MODTURB). En absence de l'entree MODTURB, 71 : l'ecoulement est laminaire. 72 : 73 : Un modele d'aspersion est disponible. 74 : 75 : Un modele de condensation en masse est en test. 76 : 77 : Au moyen du fichier d'extension dgibi, l'utilisateur transmet a CAST3M les 78 : donnees du scenario etudie. Regroupees dans la table notee RXT a differents 79 : indices qui sont precises dans cette notice, le transitoire est alors calcule 80 : par la procedure EXECRXT avec la table RXT et le nombre de pas de temps ndt en 81 : donnees d'entree : 82 : 83 : EXECRXT ndt rxt ; 84 : 85 : La table RXT est completee au moment de l'execution par trois tables : 86 : - la sous table rxt.'GEO' contient les modeles et objets geometriques crees a 87 : partir des donnees fournies. 88 : - la sous table rxt.'TBT' est la table de travail proprement dite et contient 89 : les autres objets crees necessaires au calcul hormis les inconnues. 90 : - la sous table rxt.'TIC' contient les inconnues au dernier temps connu 91 : (temps calcule ou condition initiale) ainsi que les champs variables en temps. 92 : 93 : Les entrees de RXT fournis par l'utilisateur ne sont donc pas modifiees. 94 : 95 : Les indices de RXT sont les presentes ci-dessous. 96 : 97 :
2. Objets geometriques
======================
98 : 99 : rxt . 'vtf' = GEO1 ; maillage fluide (OBLIGATOIRE) 100 : rxt . 'pi' = POI1 ; point interieur du domaine fluide ou sera imposee 101 : la pression (OBLIGATOIRE). 102 : rxt . 'axe' = GEO2 ; axe de revolution si 2D AXI 103 : 104 : On peut definir un nombre quelconque de Breches en indiquant le nom de la 105 : breche, son maillage et la direction du champ de vitesse a la breche autant 106 : de fois que necessaire : 107 : 108 : rxt . 'Breches' = 'TABLE' ; data des Breches 109 : rxt . 'Breches' . MOT1 = 'TABLE' ; data de la breche MOT1 110 : rxt . 'Breches' . MOT1 . 'Maillage' = GEO1 ; maillage surfacique de MOT1 111 : rxt . 'Breches' . MOT1 . 'diru' = POI1 ; direction de la vitesse 112 : 113 : avec MOT1 le nom de la breche choisie par l'utilisateur. 114 : 115 : On peut definir un nombre quelconque de Sorties en indiquant le nom de la 116 : sortie, son maillage et la direction du champ de vitesse a la sortie autant 117 : de fois que necessaire : 118 : 119 : rxt . 'Sorties' = 'TABLE' ; data des Sorties 120 : rxt . 'Sorties' . MOT1 = 'TABLE' ; data de la sortie MOT1 121 : rxt . 'Sorties' . MOT1 . 'Maillage' = GEO1 ; maillage surfacique de MOT1 122 : rxt . 'Sorties' . MOT1 . 'diru' = POI1 ; direction de la vitesse 123 : 124 : avec MOT1 le nom de la sortie choisi par l'utilisateur. 125 : 126 : On peut definir un nombre quelconque de PAROIS en indiquant le nom MOT1 de 127 : la paroi et son maillage autant de fois que necessaire : 128 : 129 : rxt . 'PAROIS' = 'TABLE' ; data des PAROIS 130 : rxt . 'PAROIS' . MOT1 = 'TABLE' ; data de la paroi MOT1 131 : rxt . 'PAROIS' . MOT1 . 'vtp' = GEO1 ; maillage volumique de MOT1 132 : 133 : avec MOT1 le nom de la paroi choisi par l'utilisateur. 134 : 135 : On peut definir un nombre quelconque de surface a temperature imposee en 136 : indiquant le nom MOT1 de la surface et son maillage autant de fois que 137 : necessaire : 138 : 139 : rxt . 'TIMP' = 'TABLE' ; data des surfaces a TIMP 140 : rxt . 'TIMP' . MOT1 = 'TABLE' ; data de la surface MOT1 141 : rxt . 'TIMP' . MOT1 . 'MAILLAGE' = GEO1 ; maillage surfacique de MOT1 142 : 143 : avec MOT1 le nom de la paroi a temperature imposee choisi par l'utilisateur. 144 : 145 : Exemple : On definit les breches A et B de la facon suivante : 146 : 147 : rxt . 'Breches' = 'TABL' ; 148 : rxt . 'Breches' . 'A' = 'TABL' ; 149 : rxt . 'Breches' . 'A' . 'Maillage' = surfacea ; 150 : rxt . 'Breches' . 'A' . 'diru' = (0. 0. 1.) ; 151 : rxt . 'Breches' . 'B' = 'TABL' ; 152 : rxt . 'Breches' . 'B' . 'Maillage' = surfaceb ; 153 : rxt . 'Breches' . 'B' . 'diru' = (0. -1. 0.) ; 154 : 155 :
3. Modelisation Physique
========================
156 :
3.1 Loi d'etat pour la vapeur
-----------------------------
157 : 158 : rxt . 'VERSION' = (MOT1) ; 159 : avec MOT1=('V0') si gaz parfait (par defaut), 'V1' si gaz reel (en test) 160 : 161 :
3.2 Thermique parois
--------------------
162 : 163 : A. Activation 164 : ------------- 165 : 166 : Certaines parois sont maillees en volume et on y calcule l'evolution de la 167 : temperature due aux echanges avec le fluide et le milieu ambiant. Les 168 : informations globales associees a la thermique paroi sont les suivantes : 169 : 170 : rxt . 'THERMP' = (LOG1) ; VRAI si thermique paroi (FAUX par defaut) 171 : rxt . 'THERCO' = (LOG2) ; VRAI si couplage paroi/fluide implicite (FAUX) 172 : rxt . 'ECHAN' = FLOT1 ; coefficient d'echange fluide / paroi (W/m2/oC) 173 : 174 : En cas de condensation, la valeur du coefficient d'echange convectif FLOT1 175 : est ecrasee par HTC, celle issue de la correlation associee au modele de 176 : condensation, sauf si FLOT1 lui est superieure : HT = max(FLOT1,HTC). La 177 : valeur fournie est donc un seuil. 178 : 179 : 180 : B. Proprietes des PAROIS 181 : ------------------------ 182 : 183 : Pour chaque paroi MOT1, on precise les proprietes de la paroi par : 184 : 185 : rxt . 'PAROIS' . MOT1 . 'vtp' = GEO1 ; maillage de MOT1 186 : rxt . 'PAROIS' . MOT1 . 'LAMBDA' = FLOT1 ; condu. thermique de MOT1 (W/m/oC) 187 : rxt . 'PAROIS' . MOT1 . 'ROCP' = FLOT2 ; rho*Cp de MOT1 (J/m3/oC) 188 : rxt . 'PAROIS' . MOT1 . 'Tp0' = FLOT3 ; temperature initiale de MOT1 (oC) 189 : 190 : Par exemple, on definit ainsi les parois 'Peinture' et 'Acier' : 191 : 192 : rxt . 'PAROIS' = 'TABL' ; 193 : rxt . 'PAROIS' . 'Peinture' = 'TABL' ; 194 : rxt . 'PAROIS' . 'Peinture' .'vtp' = mpeinture ; 195 : rxt . 'PAROIS' . 'Peinture' .'LAMBDA' = 0.2 ; 196 : rxt . 'PAROIS' . 'Peinture' .'ROCP' = 1.E6 ; 197 : rxt . 'PAROIS' . 'Peinture' .'Tp0' = 60. ; 198 : rxt . 'PAROIS' . 'Acier' = 'TABL' ; 199 : rxt . 'PAROIS' . 'Acier' . 'vtp' = macier ; 200 : rxt . 'PAROIS' . 'Acier' . 'LAMBDA' = 15. ; 201 : rxt . 'PAROIS' . 'Acier' . 'ROCP' = 3.9E6 ; 202 : rxt . 'PAROIS' . 'Acier' . 'Tp0' = 60. ; 203 : 204 :
3.3 Temperature de paroi imposee
--------------------------------
205 : 206 : A. Surfaces internes 207 : -------------------- 208 : 209 : Les parois internes de l'enceinte peuvent echanger via un coefficient 210 : d'echange avec le fluide, l'evolution de la temperature de la paroi etant 211 : connue en fonction du temps. 212 : 213 : rxt . 'TIMP' = 'TABL' ; data des surfaces a T imposee 214 : rxt . 'TIMP' . MOT1 = 'TABL' ; data de la surface MOT1 215 : rxt . 'TIMP' . MOT1 . 'MAILLAGE' = GEO1 ; maillage de la surface MOT1 216 : rxt . 'TIMP' . MOT1 . 'ECHAN' = FLOT1 ; coef. d'echange (W/m2/oC) 217 : rxt . 'TIMP' . MOT1 . 't' = LREE1 ; temps (s) 218 : rxt . 'TIMP' . MOT1 . 'TIMP' = LREE2 ; temperature imposee (oC) 219 : 220 : avec MOT1 le nom de la paroi a temperature imposee choisi par l'utilisateur. 221 : Cela permet en particulier de modeliser des condenseurs. 222 : 223 : En cas de condensation, la valeur du coefficient d'echange convectif FLOT1 224 : est ecrasee par HTC, celle issue de la correlation associee au modele de 225 : condensation, sauf si FLOT1 lui est superieure : HT = max(FLOT1,HTC). La 226 : valeur fournie est donc un seuil. 227 : 228 : L'ancienne structure de donnee ECHANP permettant d'imposer une surface a 229 : temperature constante au cours du temps sera supprimee a la prochaine version 230 : de CAST3M. 231 : 232 : rxt . 'ECHANP' = 'TABL' ; 233 : rxt . 'ECHANP' . 'MAILLAGE' = GEO1 ; 234 : rxt . 'ECHANP' . 'ECHAN' = FLOT1 ; 235 : rxt . 'ECHANP' . 'TMUR' = FLOT2 ; 236 : 237 : 238 : B. Surface externe 239 : ------------------ 240 : 241 : Les parois externes de l'enceinte peuvent echanger via un coefficient 242 : d'echange avec le milieu ambiant (exterieur) a condition d'activer ECHEXT. 243 : 244 : rxt . 'ECHEXT' = (LOG1) ; VRAI si echange externe paroi/milieu ambiant (FAUX) 245 : 246 : rxt . 'parext' = GEO1 ; maillage de la surface externe 247 : rxt . 'HEXT' = FLOT1 ; coefficient d'echange avec l'exterieur (W/m2/oC) 248 : rxt . 'TPEXT' = FLOT2 ; temperature exterieure (oC) 249 : 250 : Exemple : Si l'enceinte est semi immergee on peut definir deux zones (une 251 : pour la partie immergee et une pour la partie emmergee) en utilisant la 252 : procedure PROCHEXT. Les deux zones sont identifiees par la hauteur d'eau qui 253 : est une donnee du modele. Cette procedure est a appeler a partir d'une 254 : procedure PERSO (voir entree PRCPERSO). On renvoie a la notice de la procedure 255 : PROCHEXT pour renseigner correctement l'appel dans la procedure PERSO. Les 256 : donnees fournies peuvent dependre du temps (voir les jeux de donnees pressuhx1 257 : et pressuhx2.dgibi). 258 : 259 :
3.4 Conditions initiales
------------------------
260 : 261 : A. Constituants du melange 262 : -------------------------- 263 : 264 : L'air est toujours present dans le melange. 265 : 266 : rxt . 'VAPEUR' = (LOG1) ; presence de vapeur (FAUX) 267 : rxt . 'H2' = (LOG2) ; presence d'hydrogene (FAUX) 268 : rxt . 'HE' = (LOG3) ; presence d'helium (FAUX) 269 : rxt . 'N2' = (LOG4) ; presence d'azote (FAUX) 270 : rxt . 'O2' = (LOG5) ; presence d'oxygene (FAUX) 271 : rxt . 'CO' = (LOG6) ; presence d'oxyde de carbone (FAUX) 272 : rxt . 'CO2' = (LOG7) ; presence de gaz carbonique (FAUX) 273 : 274 : 275 : B. Composition du melange 276 : ------------------------- 277 : 278 : Le mélange initial est constant dans le domaine fluide. 279 : 280 : rxt . 'TF0' = FLOT1 ; temperature initiale du melange (oC) 281 : rxt . 'PT0' = FLOT2 ; pression totale initiale (Pa) 282 : rxt . MOT1 = FLOT3 ; fraction massique d'un constituant present 283 : avec MOT1='Yvap0','Yhe0','Yh20','Yo20','Yn20','Yco0','Yco20'. 284 : 285 : Remarque : L'air assure le complement a 1 des fractions massiques fournies. 286 : 287 : 288 : C. Proprietes du melange 289 : ------------------------ 290 : 291 : * Proprietes physiques du melange gazeux prises par defaut 292 : * Constantes des gaz (J/kg/K) 293 : Rgh2 = 4156.5 294 : Rghe = 2078.25 295 : Rgo2 = 259.8 296 : Rgn2 = 296.9 297 : Rgco2 = 188.9 298 : Rgco = 296.9 299 : Rgvap = 461.513 300 : Rgair = 287.1 301 : 302 : muair = 1.800e-5 : viscosite dynamique (air) (kg/m/s) 303 : alf = 1.800e-5 : diffusivite thermique (m2/s) 304 : lamair= 2.580e-2 : conductivite thermique (W/m/oC) 305 : db = 1.000e-5 : diffusion Browniene (m2/s) 306 : Cpvap = 1700.0 : chaleur specifique a pression constante 307 : pour la vapeur d'eau (J/oC/kg) 308 : Lv = 2.3E6 : chaleur latente (J/kg) 309 : 310 : Les chaleurs specifiques des autres gaz sont donnees par la 311 : procedure CALCP. 312 : 313 :
3.5 Conditions aux limites
--------------------------
314 : 315 : Les frontieres entre fluide et paroi etant traites en B/ et C/, nous 316 : indiquons les conditions aux entrees (breches) et aux sorties fluide. 317 : 318 : A. Breches 319 : ---------- 320 : 321 : Pour chaque breche MOT1, le debit et la temperature associes a differents 322 : temps sont regroupes dans la sous table scenario de la breche concernee : 323 : 324 : rxt.'Breches'.MOT1.'scenario' = 'TABL' ; scenario de la breche MOT1 325 : rxt.'Breches'.MOT1.'scenario'.'t' = LREE1 : temps du scenario (s) 326 : rxt.'Breches'.MOT1.'scenario'.'qair' = LREE2 : debit d'air (kg/s) 327 : rxt.'Breches'.MOT1.'scenario'. MOT2 = LREE3 : debit d'une espece (kg/s) 328 : rxt.'Breches'.MOT1.'scenario'.'tinj' = LREE4 : temperature d'injection (oC) 329 : 330 : avec MOT2='qeau','qhe','qh2','qo2','qn2','qco' ou 'qco2' suivant la 331 : composition du melange. Comme l'air est present par defaut, la donnee 'qair' 332 : est toujours indiquee. 333 : 334 : Par exemple, on definit une breche A de vapeur d'eau a 150oC et 50 g/s par : 335 : 336 : rxt.'Breches'.'A'.'scenario' = 'TABL' ; 337 : rxt.'Breches'.'A'.'scenario'.'t' = 'PROG' 0.0 1000.0 ; 338 : rxt.'Breches'.'A'.'scenario'.'qair' = 'PROG' 0.000 0.000 ; 339 : rxt.'Breches'.'A'.'scenario'.'qeau' = 'PROG' 0.050 0.050 ; 340 : rxt.'Breches'.'A'.'scenario'.'tinj' = 'PROG' 150.0 150.0 ; 341 : 342 : Pour la vapeur, la separation vapeur/liquide a la breche est realisee a la 343 : temperature d'ebullition. Par suite, la qualite de la vapeur a la breche est : 344 : xsi = (h_{vap,inj}-h_{liq,ebu}) / (h_{vap,sat}-h_{liq,ebu}), 345 : avec h_{vap,inj} = h(P_{vap,inj},T_{inj}), 346 : h_{vap,ebu} = hvap(P,Tsat(P)), 347 : h_{liq,ebu} = hliq(P,Tsat(P)), 348 : P, la pression totale dans l'enceinte. 349 : 350 : La partie liquide va directement dans le puisard dont la masse est 351 : initialisee avec rxt . 'Mliq0' (en kg, 0. par defaut). 352 : 353 : 354 : B. Sorties 355 : ---------- 356 : 357 : Pour chaque sortie MOT1, deux modeles sont disponibles : on modelise une 358 : ventilation forcee (entree 'Qext') ou un clapet de decharge (entree 'Pext'). 359 : Dans le premier cas, le debit de ventilation forcee est donne en fonction du 360 : temps ; dans le deuxieme cas, il est calcule en fonction de la difference 361 : entre la pression de l'enceinte et l'exterieur. Il faut alors renseigner les 362 : donnees exterieures relatives a la perte de charge : 363 : P - Pext = Kext (U ** Bext). 364 : Si les conditions sont reunies, le col est sonique (debit critique atteint). 365 : 366 : rxt . 'Sorties' . MOT1 . 'scenario' = 'TABL' ; scenario de MOT1 367 : rxt . 'Sorties' . MOT1 . 'scenario' . 't' = LREE1 ; temps (s) 368 : rxt . 'Sorties' . MOT1 . 'scenario' . 'Section' = FLOT1 ; section (m2) 369 : 370 : Cas a : 'Qext' : ventilation forcee 371 : 372 : rxt . 'Sorties' . MOT1 . 'scenario' . 'Qext' = LREE2 ; debit (kg/s) 373 : 374 : Cas b : 'Pext' : clapet de decharge 375 : 376 : rxt . 'Sorties' . MOT1 . 'scenario' . 'Pext' = LREE2 ; P exterieur (Pa) 377 : rxt . 'Sorties' . MOT1 . 'scenario' . 'Kext' = LREE3 ; Perte de charge 378 : rxt . 'Sorties' . MOT1 . 'scenario' . 'Bext' = LREE4 ; Perte de charge 379 : 380 : Remarque : Si la condition de sortie n'est pas active pendant tout le 381 : scenario, il suffit d'annuler par une procedure PERSO sa section (<1.e-10). 382 : 383 : Exemple d'une sortie appelee S : 384 : 385 : rxt . 'Sorties' = 'TABL' ; 386 : 387 : rxt . 'Sorties' . 'S' = 'TABL' ; 388 : rxt . 'Sorties' . 'S' . 'Maillage' = ms ; 389 : rxt . 'Sorties' . 'S' . 'diru' = (0. 0. 1.) ; 390 : rxt . 'Sorties' . 'S' . 'scenario' = 'TABL' ; 391 : rxt . 'Sorties' . 'S' . 'scenario' . 'Section' = 4.68e-4 ; 392 : rxt . 'Sorties' . 'S' . 'scenario' . 't' = 'PROG' 0.0 1000.0 ; 393 : 394 : Cas a : 'Qext' 395 : rxt . 'Sorties' . 'S' . 'scenario' . 'Qext' = 'PROG' 1. 1. ; 396 : 397 : Cas b : 'Pext' 398 : rxt . 'Sorties' . 'S' . 'scenario' . 'Pext' = 'PROG' 1.e5 1.e5 ; 399 : rxt . 'Sorties' . 'S' . 'scenario' . 'Kext' = 'PROG' 1.000 1.000 ; 400 : rxt . 'Sorties' . 'S' . 'scenario' . 'Bext' = 'PROG' 0.75 0.75 ; 401 : 402 : 403 : C. Flux de masse diffuse impose 404 : ------------------------------- 405 : 406 : On peut imposer les debits de masse issus d'une surface MOT1 sous forme de 407 : flux diffusif lorsque le debit est tres faible. 408 : 409 : rxt . 'XFIMP' = 'TABL' ; data des flux diffusifs 410 : rxt . 'XFIMP' . MOT1 = 'TABL' ; data de la surface MOT1 411 : rxt . 'XFIMP' . MOT1 . 'MAILLAGE' = GEO1 ; maillage 412 : rxt . 'XFIMP' . MOT1 . 't' = LREE1 ; temps (s) 413 : rxt . 'XFIMP' . MOT1 . MOT2 = LREE2 ; debit (kg/s) 414 : 415 : avec MOT1 le nom de la paroi a flux diffusif imposee choisi par l'utilisateur, 416 : MOT2='qeau','qair','qhe','qh2','qo2','qn2','qco','qco2' suivant l'espece 417 : concernee. 418 : 419 :
3.6 Turbulence
--------------
420 : 421 : L'absence de cette entree signifie que l'on est en laminaire. 422 : 423 : < rxt . 'MODTURB' = MOT1 ; > modele de turbulence choisi 424 : avec MOT1='NUTURB','LMEL' ou 'KEPSILON'. 425 : 426 : a/ rxt . 'MODTURB' = 'NUTURB' ; Nu turbulent constant 427 : rxt . 'NUT' = FLOT1 ; valeur de Nut (m2/s) 428 : 429 : b/ rxt . 'MODTURB' = 'LMEL' ; longueur de melange 430 : rxt . 'LMEL' = FLOT1 ; valeur de lm (m) 431 : 432 : c/ rxt . 'MODTURB' = 'KEPSILON' ; k-epsilon 433 : < rxt . 'FPAROI' = LOGIQUE ; > fonction de paroi (FAUX) 434 : < rxt . 'YP' = FLOT1 ; > valeur de y+ 435 : 436 : Le modele KEPSILON n'est disponible qu'avec l'algorithme semi-explicite EFM1 437 : (voir entree ALGO). 438 : 439 :
3.7 Aspersion
-------------
440 : 441 : rxt . 'ASPER' = (LOG1) ; VRAI si presence d'aspersion (FAUX) 442 : 443 : En cas d'aspersion le logique rxt.'VAPEUR' doit etre = VRAI 444 : 445 : rxt.'aspinj' = GEO1 ; Maillage de la surface d'injection 446 : de la phase dispersee 447 : rxt.'toitf' = GEO2 ; Maillage de la surface superieure (par 448 : rapport a la verticale) du volume 449 : fluide 450 : rxt.'rod' = GEO3 ; Densite de la phase dispersee 451 : rxt.'Cpd' = GEO4 ; Chaleur specifique de la phase dispersee 452 : rxt.'scenario'.'vzinj' = LREE1 ; Liste des vitesses (verticales) 453 : d'injection de la phase disp. 454 : rxt.'scenario'.'xdinj' = LREE2 ; Liste de la fraction volumique 455 : de la phase disp. a l'inject. 456 : rxt.'scenario'.'tdinj' = LREE3 ; Liste de la temperature de la 457 : phase disp. a l'injection. 458 : rxt.'scenario'.'ddinj' = LREE4 ; Liste du diametre de la phase 459 : disp. a l'injection. 460 : 461 :
3.8 Condensation en masse (en test)
-----------------------------------
462 : 463 : rxt . 'CONDMAS' = (LOG1) ; Activation de la condensation en masse (FAUX). 464 : 465 : 466 :
3.9 PAR (recombineur catalytique passif)
----------------------------------------
467 : 468 : Les recombineurs ne peuvent etre defini qu'en presence a minima des 469 : constituants suivants : H2, N2, O2, VAPEUR. 470 : 471 : On definit autant de PAR que necessaire (rang=1,n). 472 : 473 : rxt . 'RECOMB' = 'TABLE' ; data des PARs 474 : rxt . 'RECOMB' . rang = 'TABL' ; data du rangieme PAR 475 : rxt . 'RECOMB' . rang . 'PAREXT' = MAIL1 ; paroi du PAR 476 : rxt . 'RECOMB' . rang . 'ENTREE' = MAIL2 ; entree du PAR 477 : rxt . 'RECOMB' . rang . 'SORTIE' = MAIL3 ; sortie du PAR 478 : rxt . 'RECOMB' . rang . 'direntr' = POI1 ; normale a l'entree 479 : rxt . 'RECOMB' . rang . 'dirsort' = POI2 ; normale a la sortie 480 : rxt . 'RECOMB' . rang . 'NSECTION' = (ENTI1) ; nombre de PAR 481 : rxt . 'RECOMB' . rang . 'A' = (FLOT1) ; coef. taux de combustion 482 : rxt . 'RECOMB' . rang . 'B' = (FLOT2) ; coef. taux de combustion 483 : rxt . 'RECOMB' . rang . 'XH2ON' = (FLOT3) ; condition allumage 484 : rxt . 'RECOMB' . rang . 'XH2OFF' = (FLOT4) ; condition extinction 485 : rxt . 'RECOMB' . rang . 'BETA' = (FLOT5) ; parametrage 486 : 487 : Par defaut, le nombre d'unite NSECTION dans le bati est 1. 488 : 489 : Le taux de consommation d'H2 est calcule par la loi constructeur SIEMENS. 490 : Les valeurs par defaut correspondent au modele SIEMENS FR90/1-150 : 491 : - A et B : Coefficient de la loi lineaire en pression 492 : (0.48D-8 kg/s/Pa et 0.58D-3 kg/s par defaut) 493 : - XH2ON : Fraction molaire d'H2 au dela de laquelle le 494 : recombineur entre en fonctionnement (0.005 par defaut) 495 : - XH2OFF : Fraction molaire d'H2 en deça de laquelle le 496 : recombineur passe en arret (0.005 par defaut) 497 : - BETA : Zone de transition affectee par la correction 498 : pour les faibles concentrations en H2 (0.2 par defaut) 499 : 500 :
3.10 Condensation en parois
---------------------------
501 : 502 : rxt . 'MODCOND' = (MOT1) ; Modele de condensation en paroi 503 : avec MOT1=('CHIL0') ou 'CHIL1'. 504 : 505 : a/ rxt . 'MODCOND' = 'CHIL0' : Jv = Jstand = kc rho (Yv - Yvsat) 506 : Modele de type Chilton-Colburn valable pour les tres faibles concentration de 507 : vapeur. La correlation utilisee est une correlation de convection naturelle. 508 : 509 : b/ rxt . 'MODCOND' = 'CHIL1' : Jv = Jstand / (1-Yvsat) 510 : Modele deduit de la loi de Fick en masse volumique avec diffusivite de la 511 : vapeur dans le melange calculé comme dans RALOC. 512 : 513 :
4. Parametres de fonctionnement
===============================
514 :
4.1 Table RXT
-------------
515 : 516 : Il y a la possibilite d'executer une procedure personnelle en tete de la 517 : boucle en temps de l'algorithme. 518 : 519 : rxt . 'PERSO' = (LOG1) ; activation d'une procedure PERSO (FAUX) 520 : rxt . 'PRCPERSO' = MOT1 ; nom de la procedure 521 : rxt . 'TABPERSO' = TAB1 ; data associees 522 : 523 : rxt . 'DT0' = FLOT1 ; valeur du pas de temps (s) 524 : rxt . 'IMPR' = (ENTI1) ; (0) ou 1 si impressions supplementaires 525 : rxt . 'epsi' = (FLOT1) ; critere de precision des maillages (1.E-4) 526 : 527 : rxt . 'FRPREC' = (ENTI1) ; frequence de calcul du preconditionnement (1) 528 : rxt . 'FCPRECT' = (ENTI2) ; frequence de recalcul du preconditionneur base 529 : sur le rang des pas de temps (1) 530 : rxt . 'FCPRECI' = (ENTI3) ; frequence de recalcul du preconditionneur base 531 : sur le rang des iterations internes (1) 532 : 533 : Attention : si FPREC est donne, les valeurs par defaut pour FCPRECT et 534 : FCPRECI sont egales a 5 au lieu de 1. 535 : 536 : rxt . 'ALGO' = (MOT1) ; Type d'algorithme ('IMPL') 537 : avec MOT1='IMPL' si semi-implicite et 'EFM1' si semi explicite. 538 : 539 : rxt . 'DISCR' = (MOT1) ; discretisation de la vitesse et des champs 540 : scalaires (sauf la pression) ('MACRO'). 541 : rxt . 'KPRE' = (MOT2) ; discretisation de la pression ('CENTRE') 542 : 543 : Attention : tous les couples MOT1/MOT2 ne sont pas licites (voir les 544 : notices de NAVI et de MODEL). 545 : 546 : rxt . 'STAB' = (FLOT1) ; facteur de stabilisation des elements MACRO 547 : d'ancienne generation (1.). 548 : 549 : rxt . 'MAXELPRE' = (ENTI1) ; nombre maximum d'elements pour lequel on utilise 550 : une methode directe (Cholesky) pour inverser la 551 : matrice de pression (150000). 552 : rxt . 'TYPINV' = (ENTI2) ; methode de resolution iterative (voir KRES) (3) 553 : rxt . 'RESIDU' = (FLOT1) ; residu des methodes iteratives (1.E-20). 554 : rxt . 'NTIT' = (ENTI3) ; Nombre max d'iterations pour les methodes 555 : iteratives (800) 556 : 557 : rxt . 'GRAPH' = (LOG1) ; trace de controle (FAUX) 558 : rxt . 'DETMAT' = (LOG1) ; Destruction des matrices en fin du calcul (VRAI) 559 : rxt . 'REINIT' = (LOG1) ; Re-initialisation du scenario (FAUX) 560 : rxt . 'CORTEMP' = (LOG1) ; Controle du bilan d'energie (VRAI) 561 : rxt . 'CORLIM' = (LOG1) ; Controle du bilan de masse incluant CLIM (VRAI) 562 : 563 :
4.2 Limitations du modele actuel
--------------------------------
564 : 565 : Pas de modelisation gaz reel de la vapeur d'eau 566 : Pas de modelisation de la turbulence en implicite 567 : Pas de modelisation de la condensation en masse 568 : Pas de plan de symetrie possible 569 : Modele numerique adapte aux grands pas de temps 570 : Parallelisme embryonnaire 571 : 572 :
4.3 Quelques recettes
---------------------
573 : 574 : A. Modification des conditions aux limites 575 : ------------------------------------------- 576 : 577 : Lorsque par exemple une breche cesse d'emettre, la nature des conditions aux 578 : limites change a la breche. Pour prendre en compte cela, le plus simple est 579 : de sauvegarder le premier calcul, de decrire le 'nouveau' probleme sans la 580 : breche et de surcharger la table des inconnues par celle de l'ancien calcul. 581 : 582 : 1/ On recupere la table rxt.'TIC' du premier calcul 583 : 'OPTI' 'REST' 'MONFIC.sauv' ; 'REST' ; 584 : tic = rxt . 'TIC' ; 585 : 2/ On recupere eventuellemnt le maillage 586 : opti rest 'MONMAIL.sauv'; 587 : rest ; 588 : 3/ On decrit entierement le nouveau probleme 589 : rxt = 'TABL' ; 590 : rxt . 'vtf' = mon maillage ; 591 : rxt . ... etc 592 : 4/ On initialise avec la table tic precedente 593 : rxt . 'TIC' = tic ; 594 : 595 : 596 : B. Taille du fichier de sauvegarde 597 : ---------------------------------- 598 : 599 : Il est possible de reduire la taille du fichier de sauvegarde obtenu 600 : par l'opérateur 'SAUV' en detruisant les matrices (voir rxt . 'DETMAT'). 601 : 602 :
5. Resultats du calcul
======================
603 :
5.1 Table rxt.'TIC'
-------------------
604 : 605 : La table TIC contient en particulier : 606 : - des listes temporelles (LT) 0D (champs moyens) qui permettent de tracer des 607 : evolutions lorsqu'elles sont associees a la liste des temps de calcul stockee 608 : a l'indice 'LTPS' de TIC. 609 : - des champs multi-D (CHPO) calcules au dernier temps calcule. 610 : 611 : Pour chaque indice de la table, nous indiquons le type de l'indice, son nom 612 : ainsi que le type du contenu de l'indice. 613 : 614 : 615 : A. Suivi temporel 0D 616 : -------------------- 617 : 618 : Indice Objet 619 : Type & Valeur Type Commentaires 620 : 621 : MOT Tps FLOTTANT : temps physique 622 : MOT NUPADT ENTIER : rang du pas de temps Tps 623 : MOT DT FLOTTANT : pas de temps courant 624 : 625 : MOT LTPS LISTREEL : LT, liste des temps calcules 626 : 627 : MOT LMAXU LISTREEL : LT, module max de la vitesse (m/s) 628 : MOT Rhomn LISTREEL : LT, densite du melange (kg/m3) 629 : MOT Rhomv LISTREEL : LT, densite de vapeur dans le melange (kg/m3) 630 : MOT Rhom? LISTREEL : LT, densite de ? dans le melange (kg/m3) 631 : avec '?'='air','he','h2','o2','n2','co','co2' 632 : 633 : MOT PT LISTREEL : LT, pression thermodynamique 634 : MOT Tfm LISTREEL : LT, temperature moyenne du melange (oC) 635 : MOT Remn LISTREEL : LT, energie interne specifique melange (J/m3) 636 : 637 : MOT Econv LISTREEL : LT, energie extraite par convection (J/m3) 638 : MOT Econd LISTREEL : LT, energie extraite par condensation (J/m3) 639 : MOT Hcond LISTREEL : LT, enthalpie extraite par condensation (J/m3) 640 : 641 : MOT Easpe LISTREEL : LT, energie extraite par l'aspersion (J/m3) 642 : MOT Haspe LISTREEL : LT, enthalpie extraite par l'aspersion (J/m3) 643 : MOT Qaspe LISTREEL : LT, debit aspersion (kg/s) 644 : 645 : MOT Rgpm LISTREEL : LT, constante gaz parfaits melange (J/kg/K) 646 : MOT Cvm LISTREEL : LT, capacite calorifique a volume cst (J/kg/K) 647 : MOT Cpm LISTREEL : LT, capacite calorifique a p cst (J/kg/K) 648 : MOT Gamm LISTREEL : LT, rapport Cpm/Cvm 649 : 650 : MOT Minj LISTREEL : LT masse injectee (kg) 651 : MOT Mrest LISTREEL : LT masse restante (kg) 652 : 653 : MOT Mcond LISTREEL : LT masse de vapeur condensee (total) (kg) 654 : MOT Mcondw LISTREEL : LT masse de vapeur condensee (THERMP) (kg) 655 : MOT Mcond0 LISTREEL : LT masse de vapeur condensee (ECHANP) (kg) 656 : MOT Mcond1i LISTREEL : LT masse de vapeur condensee (TIMP) (kg) 657 : 658 : MOT Qc LISTREEL : LT, Debit de vapeur condensee (total) (kg/s) 659 : MOT Qcw LISTREEL : LT, Debit de vapeur condensee (THERMP) (kg/s) 660 : MOT Qc0 LISTREEL : LT, Debit de vapeur condensee (ECHANP) (kg/s) 661 : MOT Qc1i LISTREEL : LT, Debit de vapeur condensee (TIMP) (kg/s) 662 : 663 : MOT mdTf LISTREEL : LT, Correction locale sur TF (min) (oC) 664 : MOT MdTf LISTREEL : LT, Correction locale sur TF (max) (oC) 665 : MOT Mdr? LISTREEL : LT, Correction de masse pour l'espece ? (kg) 666 : avec '?'='vap','air','he','h2','o2','n2','co','co2' 667 : 668 : 669 : *--- Parois a temperature imposee (TIMP) ---* 670 : 671 : MOT TIMP TABLE : informations relatives a TIMP 672 : 673 : TIMP . nomi TABLE, data de la paroi nomi 674 : MOT Ltbpi LISTREEL : LT, temperature de la paroi nomi (oC) 675 : 676 : *--- Breches (Breches) ---* 677 : 678 : MOT Breches TABLE : informations relatives aux breches 679 : 680 : Breches . nomi TABLE, data de la breche nomi 681 : MOT guj LISTREEL : LT, vitesse a l'injection (m/s) 682 : MOT Qj LISTREEL : LT, debit d'injection (kg/s) 683 : MOT Tinj LISTREEL : LT, temperature d'injection (kg/s) 684 : MOT Hj LISTREEL : LT, enthalpie a l'injection (J/kg) 685 : MOT Ej LISTREEL : LT, energie a l'injection (J/kg) 686 : MOT Ksi LISTREEL : LT, partition liquide/vapeur (kg/s) 687 : MOT Qlj LISTREEL : LT, debit liquide a l'injection (kg/s) 688 : MOT Hlj LISTREEL : LT, enthalpie liquide a l'injection (J/kg) 689 : 690 : Usage : rxt . 'Breches' . nomi . 'Tinj' contient la LT de la 691 : temperature imposee a la breche nomi. Les temps associes sont 692 : ceux de rxt . 'TIC' . 'LTPS' 693 : 694 : *--- RECOMBINEUR ---* 695 : 696 : MOT RECOMB TABLE : informations issues des PARs 697 : 698 : RECOMB . i TABLE, data du ieme PAR 699 : MOT DEB LISTREEL : LT, debit (kg/s) 700 : MOT Uin LISTREEL : LT, vitesse a l'entree (m/s) 701 : MOT Uout LISTREEL : LT, vitesse a la sortie (m/s) 702 : MOT Tin LISTREEL : LT, temperature a l'entree (oC) 703 : MOT Tout LISTREEL : LT, temperature a la sortie (oC) 704 : MOT Tpla LISTREEL : LT, temperature des plaques (oC) 705 : MOT EFF LISTREEL : LT, efficacite du recombineur 706 : MOT Hin LISTREEL : LT, enthalpie a l'entree (J/kg) 707 : MOT Hout LISTREEL : LT, enthalpie a la sortie (J/kg) 708 : MOT Y?E LISTREEL : LT, fraction massique de l'espece ? a l'entree 709 : MOT Y?S LISTREEL : LT, fraction massique de l'espece ? a la sortie 710 : avec '?'='H2O','HE','H2','O2','N2','CO','CO2'. 711 : 712 : Usage : rxt . 'RECOMB' . 2 . 'DEB' contient la LT des debits traversant le 713 : deuxieme recombineur. Les temps associes sont ceux de rxt . 'TIC' . 'LTPS' 714 : 715 : 716 : B. Grandeurs multi-D 717 : -------------------- 718 : 719 : Indice Objet 720 : Type & Valeur Type Commentaires 721 : 722 : MOT UN CHPOINT : vitesse au pas de temps courant (m/s) 723 : MOT UNM CHPOINT : vitesse au pas de temps precedent (m/s) 724 : MOT PRES CHPOINT : pression (Pa) 725 : MOT TF CHPOINT : temperature au pas courant (oC) 726 : MOT TFNM CHPOINT : temperature au pas precedent (oC) 727 : MOT TP CHPOINT : temperature paroi au pas courant (oC) 728 : MOT TPNM CHPOINT : temperature paroi au pas precedent (oC) 729 : MOT RHO CHPOINT : densite au pas courant (kg/m3) 730 : MOT RHONM CHPOINT : densite au pas precedent (kg/m3) 731 : MOT RVP CHPOINT : densite de vapeur (kg) 732 : MOT YVAP CHPOINT : Fraction massique de vapeur 733 : 734 : MOT R? CHPOINT : densite de ? (kg) 735 : MOT Y? CHPOINT : FLOTTANT si l'espece est absente 736 : avec ?='AIR','HE','H2','O2','N2','CO','CO2' 737 : 738 : MOT Mu CHPOINT : viscosite dynamique du melange (kg/m/s) 739 : MOT NU CHPOINT : viscosite cinematique du melange (m2/s) 740 : MOT NUEFF CHPOINT : viscosite cinematique effective (m2/s) 741 : 742 : MOT KHW CHPOINT : coefficient d'echange (convectif) (THERMP) 743 : MOT KH0 CHPOINT : coefficient d'echange (convectif) (ECHANP) 744 : MOT KH1i CHPOINT : coefficient d'echange (convectif) (TIMP) 745 : 746 : 747 : *--- Aspersion ---* 748 : 749 : MOT VN CHPOINT : vitesse de la phase dispersee (m/s) 750 : MOT TD CHPOINT : temperature de la phase dispersee (oC) 751 : MOT XD CHPOINT : Fraction volumique de la phase dispersee 752 : MOT DD CHPOINT : Diametre de la phase dispersee (m) 753 : 754 :
5.2 Table rxt.'GEO'
-------------------
755 : 756 : La table GEO contient les MMODEL NAVIER-STOKES des modelisations 757 : choisies ainsi que les maillages associes. 758 : 759 : Indice Objet 760 : Type & Valeur Type Commentaires 761 : 762 : MOT epsi FLOTTANT : tolerance sur les maillages 763 : MOT $vtf MMODEL : modele du domaine fluide 764 : MOT $menvf MMODEL : modele pour la frontiere fluide 765 : MOT $axe MMODEL : axe ou plan de symetrie 766 : MOT $vtp MMODEL : paroi thermique 767 : MOT $mtp1i MMODEL : paroi temperature imposee (TIMP) 768 : MOT $mtp0 MMODEL : paroi temperature constante (ECHANP) 769 : MOT Pimp MAILLAGE : POI1 ou la pression est imposee. 770 : 771 :
© Cast3M 2003 - Tous droits réservés.
Mentions légales