* fichier : pressutq.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 1.3bar et 123oC * et la température du mur est mise à 100oC * au début du calcul. On calcule la pressurisation * de cette enceinte sur 50 secondes en injectant un débit * de 10g/s de vapeur à 130oC. 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) * Auteurs : E. Studer, J.P. Magnaud Novembre 1999 * revisite Avril 2002 *-------------------------------------------------------- COMPLET = FAUX ; *COMPLET = VRAI ; GRAPH = VRAI ; GRAPH = FAUX ; '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 ; * 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' ; *-- Nom du volume fluide rxt.'vtf' = mt ; * Definition du bidon : rxt.'TIMP1' = table ; rxt.'TIMP1'.'MAILLAGE' = mur ; rxt.'TIMP1'.'ECHAN' = 10. ; *-- Conditions initiales dans l'enceinte de test rxt.'TF0' = 123.0 ; rxt.'PT0' = 1.30000e5 ; rxt.'Yvap0' = 1.0E-8 ; *-- On positionne une brèche rxt.'breche' = jg ; rxt.'diru1' = 0. 0. 1. ; *-- 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 rxt.'scenario' = table ; *-- Conditions a la breche (Obligatoire pour l'instant) *-- 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 ; 'SI' (NON COMPLET) ; 115.21 ; 1.32798E+05 1.33293E+05 ; 0.0000 0.0000; 1.0398 1.9786 ; tic=rxt.'TIC' ; Si (ERtf '>' 1.e-3) ; erreur 5 ; Finsi ; Si (ERPT '>' 1.e-4) ; erreur 5 ; Finsi ; Si (ERQc '>' 1.e-4) ; erreur 5 ; Finsi ; Si (ERum '>' 1.e-2) ; erreur 5 ; Finsi ; 'FINSI' ; Si GRAPH ; $vtf=rxt.'GEO'.'$vtf' ; plan=plan1 et plan4 ; paroif = rxt.'GEO'.'mtp1'; 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 = rxt.'TIC'.'Fcond1'; dess evr 'TITRE' ' Rho axe '; dess evvp 'TITRE' ' rvp axe '; 'FINSI' ; 'FIN' ;
© Cast3M 2003 - Tous droits réservés.
Mentions légales