Télécharger pressuw.dgibi

Retour à la liste

Numérotation des lignes :

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

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