* fichier : dp3.dgibi ************************************************************************ * Section : Nautilus ************************************************************************ *--------------------------------------------------------- * Cas test de performace EXECRXT // RESOU et Réolutions * * Depressurisation 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.86bar et 90oC * et la température du mur est mise à 40oC * au début du calcul. On calcule la depressurisation * de cette enceinte sur 50 secondes en n'injectant pas de * vapeur. La temperature d injection est de 90oC. * Ce test (un peu long) * verifie le demarrage de la condensation " * 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 * Ce test peut utiliser la procedure enceinte * 1/ il faut l'avoir (ce n'est pas donne a tout le monde) * 2/ faire enceinte = mot 'ENCEINTE' ; *-------------------------------------------------------- COMPLET = FAUX ; GRAPH = FAUX ; *COMPLET = VRAI ; *GRAPH = VRAI ; TCPT = FAUX; TKPR = VRAI; TRESOU = VRAI; IMPARA = VRAI; 'SI' COMPLET ; nbit=100; DT0 = 1. ; n1 = 1 ; n2 = 4 ; n3 = 4 ; n4 = 8 ; nn = 2 ; n1 = 2 ; n2 = 8 ; n3 = 8 ; n4 = 16 ; nn = 4 ; n1 = 3 ; n2 = 12 ; n3 = 12 ; n4 = 24 ; nn = 6 ; n1 = 4 ; n2 = 16 ; n3 = 16 ; n4 = 32 ; nn = 8 ; 'SINON' ; nbit= 5 ; DT0 = 10. ; n1 = 1 ; n2 = 2 ; n3 = 4 ; n4 = 4 ; nn = 1 ; 'FINSI' ; *-------------------------------------------------------- * Definition du maillage de l'enceinte cylindrique * 'OPTI' 'DIME' 3 'ELEM' 'CU20' 'ECHO' 1 ; ri = 1.052 ; sp = 0.10 ; hc = 4.163 ; epsi = 1.000e-2 ; ; epsi = 1.000e-5 ; ; 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 ; l1 = 'DROI' p0 p1 n1 ; l2 = 'DROI' p1 p2 n1 ; l3 = 'DROI' p2 p3 n1 ; l4 = 'DROI' p3 p0 n1 ; l5 = 'CERC' p4 p0 p5 (2*n1) ; l6 = 'CERC' p6 p0 p7 (2*n1) ; basf0= 'DALL' l1 l2 l3 l4 'PLAN' ; basf1=('REGL' (l2 'ET' l3) l5 n2) ; l44= cote 2 basf1; ax4= (inve l4) et l44 ; l11= cote 4 basf1; ax1= l11 et (inve l1) ; basf = basf0 'ET' ('REGL' (l2 'ET' l3) l5 n2) ; 'ELIM' basf epsi ; basf = basf 'ET' ('SYME' basf 'DROI' p0 p3) ; ax11 = ('SYME' ax1 'DROI' p0 p3) 'ET' (inve ax1) ; 'ELIM' basf epsi ; basf = basf 'ET' ('SYME' basf 'DROI' p0 p1) ; ax44 = (inve ax4) 'ET' ('SYME' ax4 'DROI' p0 p4) ; 'ELIM' basf epsi ; basm = 'REGL' l5 l6 n3 ; basm = basm 'ET' ('SYME' basm 'DROI' p0 p3) ; 'ELIM' basm epsi ; basm = basm 'ET' ('SYME' basm 'DROI' p0 p1) ; 'ELIM' basm epsi ; * Creation du volume dx = ri / 2. ; nz1 = ('ENTIER' ( h1 / dx ))*nn ; bas = basf 'VOLU' nz1 'TRAN' v1 ; mbas = basm 'VOLU' nz1 'TRAN' v1 ; mbas = mbas 'COUL' 'ROUG' ; plan1 = ax11 'TRAN' nz1 v1 ; plan4 = ax44 'TRAN' nz1 v1 ; 'ELIM' (bas et plan1 et plan4) epsi ; mt = bas ; wall = mbas ; elim (mt et wall) epsi ; * Localisation d'une brèche éventuelle au bas de l'enceinte pjg = 'POIN' basf 'PROC' (0.000 0.000 0.000) ; jg = ('ELEM' basf 'APPUIE' 'LARGEMENT' pjg) 'COUL' 'VERT' ; *-------------------------------------------------------------------- * Fin de la définition du maillage *-------------------------------------------------------------------- *-------------------------------------------------------------------- * Début de l'initialisation de la procédure ENCEINTE : table RXT *-------------------------------------------------------------------- rxt = 'TABLE' ; rxt.'VERSION'= 'V0' ; rxt.'TCPT' = TCPT ; rxt.'TKPR' = TKPR ; rxt.'TRESOU' = TRESOU ; rxt.'IMPARA' = IMPARA ; list (nbno mt); *opti donn 5; *-- Nom du volume fluide rxt.'vtf' = mt ; rxt.'epsi' = epsi ; *-- 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' = 40. ; rxt.'ECHAN' = 10. ; *-- Conditions initiales dans l'enceinte de test rxt.'TF0' = 90.0 ; rxt.'PT0' = 1.86076e5 ; rxt.'Yvap0' = 0.2728 ; *-- 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) * rxt.'scenario'.'t' = prog 0.0 1000.0 ; * rxt.'scenario'.'qeau' = prog 0.000 0.000 ; * rxt.'scenario'.'qair' = prog 0.000 0.000 ; * rxt.'scenario'.'tinj' = prog 90.0 90.0 ; *-- 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 ; TEMPS; EXECRXT 2 rxt ; TEMPS; EXECRXT (nbit - 2 ) rxt ; TEMPS; list rxt.TIC.'Tfm' ; list rxt.TIC.'PT' ; list rxt.TIC.'Qc' ; list rxt.TIC.'LMAXU'; 'SI' (NON COMPLET) ; ltfm=Prog *90.000 82.759 65.078 60.601 56.047 51.777 ; *90.000 82.759 68.655 63.051 58.082 53.682 ; 90.000 82.759 68.586 62.980 57.998 53.590 ; lPT =Prog *1.86076E+05 1.83833E+05 1.73114E+05 1.55674E+05 1.49937E+05 1.47211E+05 ; *1.86076E+05 1.83833E+05 1.74918E+05 1.61022E+05 1.54803E+05 1.51174E+05 ; 1.86076E+05 1.83833E+05 1.74882E+05 1.60905E+05 1.54659E+05 1.51012E+05 ; Lqc =Prog *0.0000 7.74979E-02 4.99805E-02 3.61752E-02 3.16373E-02 2.90635E-02 ; *0.0000 5.58061E-02 4.29392E-02 3.41083E-02 3.01362E-02 2.76068E-02 ; 0.0000 5.62319E-02 4.32212E-02 3.43710E-02 3.03769E-02 2.78135E-02 ; Lmaxu=Prog *0.0000 0.0000 0.30666 0.45360 0.55440 0.30930 ; *0.0000 0.0000 0.30666 0.45360 0.49687 0.29288 ; 0.0000 0.0000 0.30666 0.45360 0.49797 0.29336 ; tic=rxt.'TIC' ; ERtf=SOMM( abs (ltfm - tic.'Tfm') )/ 80. ; ERPT=SOMM( abs (lPT - tic.'PT' ) ) /1.e5 ; ERQc=SOMM( abs (lqc - tic.'Qc' ) ) ; ERum=SOMM( abs (Lmaxu - tic.'LMAXU' ) )/2. ; Mess 'ERtf=' ERtf 'ERPT=' ERPT 'ERQc=' ERQc 'ERum=' ERum ; ierr = 0 ; Si (ERtf '>' 1.e-4) ; ierr = ierr + 1 ; Finsi ; Si (ERPT '>' 1.e-4) ; ierr = ierr + 1 ; Finsi ; Si (ERQc '>' 1.e-4) ; ierr = ierr + 1 ; Finsi ; Si (ERum '>' 1.e-3) ; ierr = ierr + 1 ; Finsi ; 'FINSI' ; Si GRAPH ; $vtf=rxt.'GEO'.'$vtf' ; mt =doma $vtf maillage ; Mpl1=chan 'QUAF' plan1 ; Mpl4=chan 'QUAF' plan4 ; ELIM (mt et Mpl1 et Mpl4) epsi ; $mpl1= mode Mpl1 'NAVIER_STOKES' MACRO ; $mpl4= mode Mpl4 'NAVIER_STOKES' MACRO ; plan1= doma $mpl1 maillage ; plan4= doma $mpl4 maillage ; plan=plan1 et plan4 ; paroif = rxt.'GEO'.'paroif'; rho=rxt.'TIC'.'RHO'; rvp=rxt.'TIC'.'RVP'; tf=rxt.'TIC'.'TF'; un=rxt.'TIC'.'UN' ; un1=redu un plan ; ung= vect un1 1. ux uy uz jaune ; trace ung plan ; opti isov suli ; 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'; trace Fcond paroif 'TITRE' ' Flux de condensation Kg / m**2 ' ; axe = p0 d nz1 (p0 plus v1) ; axe = chan axe 'QUAF' ; elim (axe et mt) epsi ; evr= evol 'CHPO' rho axe ; dess evr 'TITRE' ' Rho axe '; evt = evol 'CHPO' Tf axe ; dess evt 'TITRE' ' Température axe '; evvp= evol 'CHPO' rvp axe ; dess evvp 'TITRE' ' rvp axe '; 'FINSI' ; 'SI' (IERR '>' 0) ; 'MESS' 'Il y a des problemes !!!' ; 'ERRE' 5 ; 'SINO'; 'MESS' 'Tout s est bien passe!' ; 'FINS' ; 'FIN' ;