Télécharger dpressu.dgibi

Retour à la liste

Numérotation des lignes :

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

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