* fichier : pressu.dgibi *--------------------------------------------------------- * Pressurisation d'une enceinte type Phébus * * Le maillage correspond à une enceinte cylindrique * d'environ 10 m3 avec un mur en contact avec la * paroi verticale (10 cm d'acier) * Tout le volume est initialement a 1bar et 40oC * et la température du mur est mise à 60oC * au début du calcul. On calcule la pressurisation * de cette enceinte sur 50 secondes en injectant un débit * de 50g/s de vapeur à 150oC. Ce test (un peu long) * verifie le demarrage de la condensation a 20" * verifie l'evolution moyenne de la temprature gaz * verifie la pression max a 50" * verifie la vitesse max a 50" (Convection naturelle) * verifie la masse d'eau a 50" * Auteurs : E. Studer, J.P. Magnaud Novembre 1999 * revisite Mars 2001 *-------------------------------------------------------- COMPLET = FAUX ; *COMPLET = VRAI ; *GRAPH = VRAI ; GRAPH = FAUX ; *'OPTI' 'TRAC' 'X' ; 'SI' COMPLET ; nbit=100; DT0 = 1. ; n1 = 1 ; n2 = 4 ; n3 = 4 ; n4 = 8 ; nn = 2 ; 'SINON' ; nbit= 5 ; DT0 = 10. ; n1 = 1 ; n2 = 2 ; n3 = 4 ; n4 = 4 ; nn = 1 ; 'FINSI' ; *-------------------------------------------------------- * Definition du maillage de l'enceinte cylindrique * ri = 1.052 ; sp = 0.10 ; hc = 4.163 ; p0 = 0.000 0.000 0.000 ; px = -1000.000 0.000 0.000 ; py = 0.000 -1000.000 0.000 ; pz = 0.000 0.000 1000.000 ; cd = 0.000 0.000 -20.000 ; ph0 = 0.000 0.000 hc ; phx = ri 0.000 hc ; phy = 0.000 ri hc ; fg1 = 0.25 ; fg2 = fg1 * (2.0 ** 0.5) / 2. ; p1 = (ri*fg1) 0.000 0.000 ; p2 = (ri*fg2) (ri*fg2) 0.000 ; p3 = 0.000 (ri*fg1) 0.000 ; p4 = ri 0.000 0.000 ; p5 = 0.000 ri 0.000 ; p6 = (ri+sp) 0.000 0.000 ; p7 = 0.000 (ri+sp) 0.000 ; * Hauteur de l'enceinte h1 = 4.163 ; * Vecteur de translation v1 = 0. 0. h1 ; * Creation du volume dx = ri / 2. ; nz1 = ('ENTIER' ( h1 / dx ))*nn ; mt = bas ; wall = mbas ; * Localisation d'une brèche éventuelle au bas de l'enceinte *-------------------------------------------------------------------- * Fin de la définition du maillage *-------------------------------------------------------------------- *-------------------------------------------------------------------- * Début de l'initialisation de la procédure ENCEINTE : table RXT *-------------------------------------------------------------------- rxt = 'TABLE' ; rxt.'VERSION'= 'V0' ; *-- Nom du volume fluide rxt.'vtf' = mt ; *-- Definition des murs de l'enceinte : ici un seul mur *-- en ACIER dont on traitera la thermique dans l'épaisseur *-- et que l'on initialise a 40oC *-- On definit d'abort la matériau ACIER avec sa conductivite *-- thermique LAMBDA (W/m/K) et le produit ro*Cp (J/m3/K) rxt.'THERMP' = VRAI ; rxt.'vtp' = wall ; rxt.'LAMBDA' = 15. ; rxt.'ROCP' = 3.9E6 ; rxt.'Tp0' = 60. ; rxt.'ECHAN' = 10. ; *-- Conditions initiales dans l'enceinte de test rxt.'TF0' = 40.0 ; rxt.'PT0' = 1.00000e5 ; rxt.'Yvap0' = 0.0023 ; *-- On definit un point interne au maillage pour imposer la valeur de *-- la pression rxt.'pi' = (0.0 0.0 0.5) ; *-- On indique que le calcul comporte de la vapeur d'eau rxt.'VAPEUR' = VRAI ; *-- On active le recalcul automatique du préconditionnement *-- toutes les 5 itérations rxt.'FRPREC' = 5 ; rxt.'DETMAT' = VRAI ; *-- Definition du scenario thermohydraulique Si FAUX; * Ancienne méthode rxt.'breche' = jg ; rxt.'diru1' = 0 0 1 ; rxt.'scenario' = table ; Sinon ; * Nouvelle méthode rxt.'Breches' = table ; rxt.'Breches'.'A' = table ; rxt.'Breches'.'A'.'scenario' = table ; rxt.'Breches'.'A'.'Maillage' = jg ; rxt.'Breches'.'A'.'diru' = (0. 0. 1.) ; Finsi ; *-- On impose le pas de temps (s) rxt.'DT0' = DT0 ; *-- On impose la viscosite turbulente (m2/s) rxt.'MODTURB'='NUTURB' ; rxt.'NUT' = 1.e-2 ; *-- On lance le calcul sur 20 itérations d'une seconde rxt.'GRAPH'=GRAPH ; ltfm='PROG' 40.000 65.434 73.712 81.158 86.814 90.738 ; lPT ='PROG' 1.00000E+05 1.06009E+05 1.21199E+05 1.29592E+05 1.36462E+05 1.43998E+05 ; Lqc ='PROG' 0.0000 0.0000 0.0000 3.81471E-04 2.57962E-03 4.78997E-03 ; Lmaxu='PROG' 0.0000 0.81320 2.0813 2.6901 2.2350 2.5277 ; tic=rxt.'TIC' ; Si (ERtf '>' 1.e-4) ; erreur 5 ; Finsi ; Si (ERPT '>' 1.e-3) ; erreur 5 ; Finsi ; Si (ERQc '>' 1.e-4) ; erreur 5 ; Finsi ; Si (ERum '>' 1.e-2) ; erreur 5 ; Finsi ; Si GRAPH ; $vtf=rxt.'GEO'.'$vtf' ; plan=plan1 et plan4 ; paroif = rxt.'GEO'.'paroif'; rho=rxt.'TIC'.'RHO'; rvp=rxt.'TIC'.'RVP'; tf=rxt.'TIC'.'TF'; un=rxt.'TIC'.'UN' ; trace ung plan ; trace rho plan 'TITRE'' Rho' ; trace rvp plan 'TITRE'' Rvp' ; trace tf plan 'TITRE'' Tf ' ; trace rho paroif 'TITRE'' Rho' ; trace rvp paroif 'TITRE'' Rvp' ; trace Tf paroif 'TITRE'' Tf ' ; Fcond = tic.'Fcondw'; dess evr 'TITRE' ' Rho axe '; dess evvp 'TITRE' ' rvp axe '; 'FINSI' ; 'FIN' ;
© Cast3M 2003 - Tous droits réservés.
Mentions légales