Télécharger pressu2.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : pressu.dgibi
  2. *---------------------------------------------------------
  3. * Pressurisation d'une enceinte type Phébus
  4. *
  5. * Le maillage correspond à une enceinte cylindrique
  6. * d'environ 10 m3 avec un mur en contact avec la
  7. * paroi verticale (10 cm d'acier)
  8. * Tout le volume est initialement a 1bar et 40oC
  9. * et la température du mur est mise à 60oC
  10. * au début du calcul. On calcule la pressurisation
  11. * de cette enceinte sur 50 secondes en injectant un débit
  12. * de 50g/s de vapeur à 150oC. Ce test (un peu long)
  13. * verifie le demarrage de la condensation a 20"
  14. * verifie l'evolution moyenne de la temprature gaz
  15. * verifie la pression max a 50"
  16. * verifie la vitesse max a 50" (Convection naturelle)
  17. * verifie la masse d'eau a 50"
  18. * Auteurs : E. Studer, J.P. Magnaud Novembre 1999
  19. * revisite Mars 2001
  20. *--------------------------------------------------------
  21. COMPLET = FAUX ;
  22. *COMPLET = VRAI ;
  23. *GRAPH = VRAI ;
  24. GRAPH = FAUX ;
  25. *'OPTI' 'TRAC' 'X' ;
  26. 'OPTI' 'TRAC' 'PSC' ;
  27.  
  28. 'SI' COMPLET ;
  29.  
  30. nbit=100;
  31. DT0 = 1. ;
  32. n1 = 1 ; n2 = 4 ; n3 = 4 ;
  33. n4 = 8 ; nn = 2 ;
  34.  
  35. 'SINON' ;
  36.  
  37. nbit= 5 ;
  38. DT0 = 10. ;
  39. n1 = 1 ; n2 = 2 ; n3 = 4 ;
  40. n4 = 4 ; nn = 1 ;
  41.  
  42. 'FINSI' ;
  43.  
  44. *--------------------------------------------------------
  45. * Definition du maillage de l'enceinte cylindrique
  46. *
  47. 'OPTI' 'DIME' 3 'ELEM' 'CU20' ;
  48.  
  49. ri = 1.052 ; sp = 0.10 ; hc = 4.163 ;
  50.  
  51.  
  52.  
  53. epsi = 1.000e-2 ; ;
  54.  
  55. p0 = 0.000 0.000 0.000 ;
  56. px = -1000.000 0.000 0.000 ;
  57. py = 0.000 -1000.000 0.000 ;
  58. pz = 0.000 0.000 1000.000 ;
  59. cd = 0.000 0.000 -20.000 ;
  60. ph0 = 0.000 0.000 hc ;
  61. phx = ri 0.000 hc ;
  62. phy = 0.000 ri hc ;
  63.  
  64. fg1 = 0.25 ;
  65. fg2 = fg1 * (2.0 ** 0.5) / 2. ;
  66.  
  67. p1 = (ri*fg1) 0.000 0.000 ;
  68. p2 = (ri*fg2) (ri*fg2) 0.000 ;
  69. p3 = 0.000 (ri*fg1) 0.000 ;
  70. p4 = ri 0.000 0.000 ;
  71. p5 = 0.000 ri 0.000 ;
  72. p6 = (ri+sp) 0.000 0.000 ;
  73. p7 = 0.000 (ri+sp) 0.000 ;
  74.  
  75. * Hauteur de l'enceinte
  76. h1 = 4.163 ;
  77. * Vecteur de translation
  78. v1 = 0. 0. h1 ;
  79.  
  80. l1 = 'DROI' p0 p1 n1 ;
  81. l2 = 'DROI' p1 p2 n1 ;
  82. l3 = 'DROI' p2 p3 n1 ;
  83. l4 = 'DROI' p3 p0 n1 ;
  84. l5 = 'CERC' p4 p0 p5 (2*n1) ;
  85. l6 = 'CERC' p6 p0 p7 (2*n1) ;
  86.  
  87. basf0= 'DALL' l1 l2 l3 l4 'PLAN' ;
  88. basf1=('REGL' (l2 'ET' l3) l5 n2) ;
  89. l44= cote 2 basf1;
  90. ax4= (inve l4) et l44 ;
  91. l11= cote 4 basf1;
  92. ax1= l11 et (inve l1) ;
  93.  
  94. basf = basf0 'ET' ('REGL' (l2 'ET' l3) l5 n2) ;
  95. 'ELIM' basf epsi ;
  96. basf = basf 'ET' ('SYME' basf 'DROI' p0 p3) ;
  97. ax11 = ('SYME' ax1 'DROI' p0 p3) 'ET' (inve ax1) ;
  98. 'ELIM' basf epsi ;
  99. basf = basf 'ET' ('SYME' basf 'DROI' p0 p1) ;
  100. ax44 = (inve ax4) 'ET' ('SYME' ax4 'DROI' p0 p4) ;
  101. 'ELIM' basf epsi ;
  102.  
  103. basm = 'REGL' l5 l6 n3 ;
  104. basm = basm 'ET' ('SYME' basm 'DROI' p0 p3) ;
  105. 'ELIM' basm epsi ;
  106. basm = basm 'ET' ('SYME' basm 'DROI' p0 p1) ;
  107. 'ELIM' basm epsi ;
  108.  
  109. * Creation du volume
  110.  
  111. dx = ri / 2. ;
  112. nz1 = ('ENTIER' ( h1 / dx ))*nn ;
  113.  
  114. bas = basf 'VOLU' nz1 'TRAN' v1 ;
  115. mbas = basm 'VOLU' nz1 'TRAN' v1 ;
  116. mbas = mbas 'COUL' 'ROUG' ;
  117. plan1 = ax11 'TRAN' nz1 v1 ;
  118. plan4 = ax44 'TRAN' nz1 v1 ;
  119. 'ELIM' (bas et plan1 et plan4) epsi ;
  120.  
  121. mt = bas ;
  122. wall = mbas ;
  123. elim (mt et wall) epsi ;
  124.  
  125. * Localisation d'une brèche éventuelle au bas de l'enceinte
  126.  
  127. pjg = 'POIN' basf 'PROC' (0.000 0.000 0.000) ;
  128. jg = ('ELEM' basf 'APPUIE' 'LARGEMENT' pjg) 'COUL' 'VERT' ;
  129.  
  130. *--------------------------------------------------------------------
  131. * Fin de la définition du maillage
  132. *--------------------------------------------------------------------
  133.  
  134. *--------------------------------------------------------------------
  135. * Début de l'initialisation de la procédure ENCEINTE : table RXT
  136. *--------------------------------------------------------------------
  137.  
  138. rxt = 'TABLE' ;
  139. rxt.'VERSION'= 'V0' ;
  140.  
  141. *-- Nom du volume fluide
  142. rxt.'vtf' = mt ;
  143. rxt.'epsi' = epsi ;
  144.  
  145. *-- Definition des murs de l'enceinte : ici un seul mur
  146. *-- en ACIER dont on traitera la thermique dans l'épaisseur
  147. *-- et que l'on initialise a 40oC
  148.  
  149. *-- On definit d'abort la matériau ACIER avec sa conductivite
  150. *-- thermique LAMBDA (W/m/K) et le produit ro*Cp (J/m3/K)
  151. rxt.'THERMP' = VRAI ;
  152. rxt.'vtp' = wall ;
  153. rxt.'LAMBDA' = 15. ;
  154. rxt.'ROCP' = 3.9E6 ;
  155. rxt.'Tp0' = 60. ;
  156. rxt.'ECHAN' = 10. ;
  157.  
  158. *-- Conditions initiales dans l'enceinte de test
  159. rxt.'TF0' = 40.0 ;
  160. rxt.'PT0' = 1.00000e5 ;
  161. rxt.'Yvap0' = 0.0023 ;
  162.  
  163.  
  164. *-- On definit un point interne au maillage pour imposer la valeur de
  165. *-- la pression
  166. rxt.'pi' = (0.0 0.0 0.5) ;
  167.  
  168. *-- On indique que le calcul comporte de la vapeur d'eau
  169. rxt.'VAPEUR' = VRAI ;
  170.  
  171. *-- On active le recalcul automatique du préconditionnement
  172. *-- toutes les 5 itérations
  173. rxt.'FRPREC' = 5 ;
  174. rxt.'DETMAT' = VRAI ;
  175.  
  176. *-- Definition du scenario thermohydraulique
  177. Si FAUX;
  178. * Ancienne méthode
  179. rxt.'breche' = jg ;
  180. rxt.'diru1' = 0 0 1 ;
  181. rxt.'scenario' = table ;
  182. rxt.'scenario'.'t' = prog 0.0 1000.0 ;
  183. rxt.'scenario'.'qeau' = prog 0.050 0.050 ;
  184. rxt.'scenario'.'qair' = prog 0.000 0.000 ;
  185. rxt.'scenario'.'tinj' = prog 150.0 150.0 ;
  186. Sinon ;
  187. * Nouvelle méthode
  188. rxt.'Breches' = table ;
  189. rxt.'Breches'.'A' = table ;
  190. rxt.'Breches'.'A'.'scenario' = table ;
  191. rxt.'Breches'.'A'.'Maillage' = jg ;
  192. rxt.'Breches'.'A'.'diru' = (0. 0. 1.) ;
  193. rxt.'Breches'.'A'.'scenario'.'t' = prog 0.0 1000.0 ;
  194. rxt.'Breches'.'A'.'scenario'.'qeau' = prog 0.050 0.050 ;
  195. rxt.'Breches'.'A'.'scenario'.'qair' = prog 0.000 0.000 ;
  196. rxt.'Breches'.'A'.'scenario'.'tinj' = prog 150.0 150.0 ;
  197. Finsi ;
  198.  
  199. *-- On impose le pas de temps (s)
  200. rxt.'DT0' = DT0 ;
  201.  
  202. *-- On impose la viscosite turbulente (m2/s)
  203. rxt.'MODTURB'='NUTURB' ;
  204. rxt.'NUT' = 1.e-2 ;
  205.  
  206. *-- On lance le calcul sur 20 itérations d'une seconde
  207. rxt.'GRAPH'=GRAPH ;
  208.  
  209. EXECRXT 2 rxt ;
  210. EXECRXT (nbit - 2) rxt ;
  211.  
  212. list rxt.TIC.'Tfm' ;
  213. list rxt.TIC.'PT' ;
  214. list rxt.TIC.'Qc' ;
  215. list rxt.TIC.'LMAXU';
  216.  
  217.  
  218. ltfm='PROG'
  219. 40.000 65.434 73.712 81.158 86.814
  220. 90.738 ;
  221.  
  222. lPT ='PROG'
  223. 1.00000E+05 1.06009E+05 1.21199E+05 1.29592E+05 1.36462E+05
  224. 1.43998E+05 ;
  225. Lqc ='PROG'
  226. 0.0000 0.0000 0.0000 3.81471E-04 2.57962E-03
  227. 4.78997E-03 ;
  228. Lmaxu='PROG'
  229. 0.0000 0.81320 2.0813 2.6901 2.2350
  230. 2.5277 ;
  231.  
  232. tic=rxt.'TIC' ;
  233. ERtf=somm( abs (ltfm - tic.'Tfm') )/ 80. ;
  234. ERPT=somm( abs (lPT - tic.'PT' ) ) /1.e5 ;
  235. ERQc=somm( abs (lqc - tic.'Qc' ) ) ;
  236. ERum=somm( abs (Lmaxu - tic.'LMAXU' ) )/2. ;
  237.  
  238. Mess 'ERtf=' ERtf 'ERPT=' ERPT 'ERQc=' ERQc 'ERum=' ERum ;
  239.  
  240. Si (ERtf '>' 1.e-4) ; erreur 5 ; Finsi ;
  241. Si (ERPT '>' 1.e-3) ; erreur 5 ; Finsi ;
  242. Si (ERQc '>' 1.e-4) ; erreur 5 ; Finsi ;
  243. Si (ERum '>' 1.e-2) ; erreur 5 ; Finsi ;
  244.  
  245. Si GRAPH ;
  246.  
  247. $vtf=rxt.'GEO'.'$vtf' ;
  248. mt =doma $vtf maillage ;
  249. Mpl1=chan 'QUAF' plan1 ;
  250. Mpl4=chan 'QUAF' plan4 ;
  251. ELIM (mt et Mpl1 et Mpl4) epsi ;
  252. $mpl1= mode Mpl1 'NAVIER_STOKES' MACRO ;
  253. $mpl4= mode Mpl4 'NAVIER_STOKES' MACRO ;
  254. plan1= doma $mpl1 maillage ;
  255. plan4= doma $mpl4 maillage ;
  256. plan=plan1 et plan4 ;
  257.  
  258. paroif = rxt.'GEO'.'paroif';
  259. rho=rxt.'TIC'.'RHO';
  260. rvp=rxt.'TIC'.'RVP';
  261. tf=rxt.'TIC'.'TF';
  262. un=rxt.'TIC'.'UN' ;
  263. un1=redu un plan ;
  264. ung= vect un1 1. ux uy uz jaune ;
  265. trace ung plan ;
  266. opti isov suli ;
  267. trace rho plan 'TITRE'' Rho' ;
  268. trace rvp plan 'TITRE'' Rvp' ;
  269. trace tf plan 'TITRE'' Tf ' ;
  270.  
  271. trace rho paroif 'TITRE'' Rho' ;
  272. trace rvp paroif 'TITRE'' Rvp' ;
  273. trace Tf paroif 'TITRE'' Tf ' ;
  274. Fcond = tic.'Fcondw';
  275. trace Fcond paroif 'TITRE' ' Flux de condensation Kg / m**2 ' ;
  276.  
  277. axe = p0 d nz1 (p0 plus v1) ;
  278. axe = chan axe 'QUAF' ;
  279. elim (axe et mt) epsi ;
  280.  
  281. evr= evol 'CHPO' rho axe ;
  282. dess evr 'TITRE' ' Rho axe ';
  283.  
  284. evt = evol 'CHPO' Tf axe ;
  285. dess evt 'TITRE' ' Température axe ';
  286.  
  287. evvp= evol 'CHPO' rvp axe ;
  288. dess evvp 'TITRE' ' rvp axe ';
  289.  
  290.  
  291.  
  292. 'FINSI' ;
  293.  
  294.  
  295.  
  296.  
  297.  
  298. 'FIN' ;
  299.  
  300.  
  301.  
  302.  
  303.  
  304.  
  305.  
  306.  
  307.  
  308.  
  309.  
  310.  
  311.  
  312.  
  313.  
  314.  
  315.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales