Télécharger dp3.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : dp3.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *---------------------------------------------------------
  5. * Cas test de performace EXECRXT // RESOU et Réolutions
  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. GRAPH = FAUX ;
  30. *COMPLET = VRAI ;
  31. *GRAPH = VRAI ;
  32.  
  33. TCPT = FAUX;
  34. TKPR = VRAI;
  35. TRESOU = VRAI;
  36. IMPARA = VRAI;
  37.  
  38. 'SI' COMPLET ;
  39.  
  40. nbit=100;
  41. DT0 = 1. ;
  42. n1 = 1 ; n2 = 4 ; n3 = 4 ;
  43. n4 = 8 ; nn = 2 ;
  44.  
  45. n1 = 2 ; n2 = 8 ; n3 = 8 ;
  46. n4 = 16 ; nn = 4 ;
  47.  
  48. n1 = 3 ; n2 = 12 ; n3 = 12 ;
  49. n4 = 24 ; nn = 6 ;
  50.  
  51. n1 = 4 ; n2 = 16 ; n3 = 16 ;
  52. n4 = 32 ; nn = 8 ;
  53.  
  54. 'SINON' ;
  55.  
  56. nbit= 5 ;
  57. DT0 = 10. ;
  58. n1 = 1 ; n2 = 2 ; n3 = 4 ;
  59. n4 = 4 ; nn = 1 ;
  60.  
  61. 'FINSI' ;
  62.  
  63. *--------------------------------------------------------
  64. * Definition du maillage de l'enceinte cylindrique
  65. *
  66. 'OPTI' 'DIME' 3 'ELEM' 'CU20' 'ECHO' 1 ;
  67.  
  68. ri = 1.052 ; sp = 0.10 ; hc = 4.163 ;
  69.  
  70.  
  71.  
  72. epsi = 1.000e-2 ; ;
  73. epsi = 1.000e-5 ; ;
  74.  
  75. p0 = 0.000 0.000 0.000 ;
  76. px = -1000.000 0.000 0.000 ;
  77. py = 0.000 -1000.000 0.000 ;
  78. pz = 0.000 0.000 1000.000 ;
  79. cd = 0.000 0.000 -20.000 ;
  80. ph0 = 0.000 0.000 hc ;
  81. phx = ri 0.000 hc ;
  82. phy = 0.000 ri hc ;
  83.  
  84. fg1 = 0.25 ;
  85. fg2 = fg1 * (2.0 ** 0.5) / 2. ;
  86.  
  87. p1 = (ri*fg1) 0.000 0.000 ;
  88. p2 = (ri*fg2) (ri*fg2) 0.000 ;
  89. p3 = 0.000 (ri*fg1) 0.000 ;
  90. p4 = ri 0.000 0.000 ;
  91. p5 = 0.000 ri 0.000 ;
  92. p6 = (ri+sp) 0.000 0.000 ;
  93. p7 = 0.000 (ri+sp) 0.000 ;
  94.  
  95. * Hauteur de l'enceinte
  96. h1 = 4.163 ;
  97. * Vecteur de translation
  98. v1 = 0. 0. h1 ;
  99.  
  100. l1 = 'DROI' p0 p1 n1 ;
  101. l2 = 'DROI' p1 p2 n1 ;
  102. l3 = 'DROI' p2 p3 n1 ;
  103. l4 = 'DROI' p3 p0 n1 ;
  104. l5 = 'CERC' p4 p0 p5 (2*n1) ;
  105. l6 = 'CERC' p6 p0 p7 (2*n1) ;
  106.  
  107. basf0= 'DALL' l1 l2 l3 l4 'PLAN' ;
  108. basf1=('REGL' (l2 'ET' l3) l5 n2) ;
  109. l44= cote 2 basf1;
  110. ax4= (inve l4) et l44 ;
  111. l11= cote 4 basf1;
  112. ax1= l11 et (inve l1) ;
  113.  
  114. basf = basf0 'ET' ('REGL' (l2 'ET' l3) l5 n2) ;
  115. 'ELIM' basf epsi ;
  116. basf = basf 'ET' ('SYME' basf 'DROI' p0 p3) ;
  117. ax11 = ('SYME' ax1 'DROI' p0 p3) 'ET' (inve ax1) ;
  118. 'ELIM' basf epsi ;
  119. basf = basf 'ET' ('SYME' basf 'DROI' p0 p1) ;
  120. ax44 = (inve ax4) 'ET' ('SYME' ax4 'DROI' p0 p4) ;
  121. 'ELIM' basf epsi ;
  122.  
  123. basm = 'REGL' l5 l6 n3 ;
  124. basm = basm 'ET' ('SYME' basm 'DROI' p0 p3) ;
  125. 'ELIM' basm epsi ;
  126. basm = basm 'ET' ('SYME' basm 'DROI' p0 p1) ;
  127. 'ELIM' basm epsi ;
  128.  
  129. * Creation du volume
  130.  
  131. dx = ri / 2. ;
  132. nz1 = ('ENTIER' ( h1 / dx ))*nn ;
  133.  
  134. bas = basf 'VOLU' nz1 'TRAN' v1 ;
  135. mbas = basm 'VOLU' nz1 'TRAN' v1 ;
  136. mbas = mbas 'COUL' 'ROUG' ;
  137. plan1 = ax11 'TRAN' nz1 v1 ;
  138. plan4 = ax44 'TRAN' nz1 v1 ;
  139. 'ELIM' (bas et plan1 et plan4) epsi ;
  140.  
  141. mt = bas ;
  142. wall = mbas ;
  143. elim (mt et wall) epsi ;
  144.  
  145. * Localisation d'une brèche éventuelle au bas de l'enceinte
  146.  
  147. pjg = 'POIN' basf 'PROC' (0.000 0.000 0.000) ;
  148. jg = ('ELEM' basf 'APPUIE' 'LARGEMENT' pjg) 'COUL' 'VERT' ;
  149.  
  150. *--------------------------------------------------------------------
  151. * Fin de la définition du maillage
  152. *--------------------------------------------------------------------
  153.  
  154. *--------------------------------------------------------------------
  155. * Début de l'initialisation de la procédure ENCEINTE : table RXT
  156. *--------------------------------------------------------------------
  157.  
  158. rxt = 'TABLE' ;
  159. rxt.'VERSION'= 'V0' ;
  160. rxt.'TCPT' = TCPT ;
  161. rxt.'TKPR' = TKPR ;
  162. rxt.'TRESOU' = TRESOU ;
  163. rxt.'IMPARA' = IMPARA ;
  164.  
  165. list (nbno mt);
  166. *opti donn 5;
  167. *-- Nom du volume fluide
  168. rxt.'vtf' = mt ;
  169. rxt.'epsi' = epsi ;
  170.  
  171. *-- Definition des murs de l'enceinte : ici un seul mur
  172. *-- en ACIER dont on traitera la thermique dans l'épaisseur
  173. *-- et que l'on initialise a 40oC
  174.  
  175. *-- On definit d'abort la matériau ACIER avec sa conductivite
  176. *-- thermique LAMBDA (W/m/K) et le produit ro*Cp (J/m3/K)
  177. rxt.'THERMP' = VRAI ;
  178. rxt.'vtp' = wall ;
  179. rxt.'LAMBDA' = 15. ;
  180. rxt.'ROCP' = 3.9E6 ;
  181. rxt.'Tp0' = 40. ;
  182. rxt.'ECHAN' = 10. ;
  183.  
  184. *-- Conditions initiales dans l'enceinte de test
  185. rxt.'TF0' = 90.0 ;
  186. rxt.'PT0' = 1.86076e5 ;
  187. rxt.'Yvap0' = 0.2728 ;
  188.  
  189. *-- On positionne une brèche
  190. * rxt.'breche' = jg ;
  191. * rxt.'diru1' = 0 0 1 ;
  192.  
  193. *-- On definit un point interne au maillage pour imposer la valeur de
  194. *-- la pression
  195. rxt.'pi' = (0.0 0.0 0.5) ;
  196.  
  197. *-- On indique que le calcul comporte de la vapeur d'eau
  198. rxt.'VAPEUR' = VRAI ;
  199.  
  200. *-- On active le recalcul automatique du préconditionnement
  201. *-- toutes les 5 itérations
  202. rxt.'FRPREC' = 5 ;
  203. rxt.'DETMAT' = VRAI ;
  204.  
  205. *-- Definition du scenario thermohydraulique
  206. * rxt.'scenario' = table ;
  207.  
  208. *-- Conditions a la breche (Obligatoire pour l'instant)
  209. * rxt.'scenario'.'t' = prog 0.0 1000.0 ;
  210. * rxt.'scenario'.'qeau' = prog 0.000 0.000 ;
  211. * rxt.'scenario'.'qair' = prog 0.000 0.000 ;
  212. * rxt.'scenario'.'tinj' = prog 90.0 90.0 ;
  213.  
  214. *-- On impose le pas de temps (s)
  215. rxt.'DT0' = DT0 ;
  216.  
  217. *-- On impose la viscosite turbulente (m2/s)
  218. rxt.'MODTURB'='NUTURB' ;
  219. rxt.'NUT' = 1.e-2 ;
  220.  
  221. *-- On lance le calcul sur 20 itérations d'une seconde
  222. rxt.'GRAPH'=GRAPH ;
  223.  
  224. TEMPS;
  225. EXECRXT 2 rxt ;
  226. TEMPS;
  227. EXECRXT (nbit - 2 ) rxt ;
  228. TEMPS;
  229.  
  230. list rxt.TIC.'Tfm' ;
  231. list rxt.TIC.'PT' ;
  232. list rxt.TIC.'Qc' ;
  233. list rxt.TIC.'LMAXU';
  234.  
  235. 'SI' (NON COMPLET) ;
  236.  
  237. ltfm=Prog
  238. *90.000 82.759 65.078 60.601 56.047 51.777 ;
  239. *90.000 82.759 68.655 63.051 58.082 53.682 ;
  240. 90.000 82.759 68.586 62.980 57.998 53.590 ;
  241.  
  242.  
  243. lPT =Prog
  244. *1.86076E+05 1.83833E+05 1.73114E+05 1.55674E+05 1.49937E+05 1.47211E+05 ;
  245. *1.86076E+05 1.83833E+05 1.74918E+05 1.61022E+05 1.54803E+05 1.51174E+05 ;
  246. 1.86076E+05 1.83833E+05 1.74882E+05 1.60905E+05 1.54659E+05 1.51012E+05 ;
  247.  
  248. Lqc =Prog
  249. *0.0000 7.74979E-02 4.99805E-02 3.61752E-02 3.16373E-02 2.90635E-02 ;
  250. *0.0000 5.58061E-02 4.29392E-02 3.41083E-02 3.01362E-02 2.76068E-02 ;
  251. 0.0000 5.62319E-02 4.32212E-02 3.43710E-02 3.03769E-02 2.78135E-02 ;
  252.  
  253. Lmaxu=Prog
  254. *0.0000 0.0000 0.30666 0.45360 0.55440 0.30930 ;
  255. *0.0000 0.0000 0.30666 0.45360 0.49687 0.29288 ;
  256. 0.0000 0.0000 0.30666 0.45360 0.49797 0.29336 ;
  257.  
  258. tic=rxt.'TIC' ;
  259. ERtf=SOMM( abs (ltfm - tic.'Tfm') )/ 80. ;
  260. ERPT=SOMM( abs (lPT - tic.'PT' ) ) /1.e5 ;
  261. ERQc=SOMM( abs (lqc - tic.'Qc' ) ) ;
  262. ERum=SOMM( abs (Lmaxu - tic.'LMAXU' ) )/2. ;
  263.  
  264. Mess 'ERtf=' ERtf 'ERPT=' ERPT 'ERQc=' ERQc 'ERum=' ERum ;
  265.  
  266. ierr = 0 ;
  267. Si (ERtf '>' 1.e-4) ; ierr = ierr + 1 ; Finsi ;
  268. Si (ERPT '>' 1.e-4) ; ierr = ierr + 1 ; Finsi ;
  269. Si (ERQc '>' 1.e-4) ; ierr = ierr + 1 ; Finsi ;
  270. Si (ERum '>' 1.e-3) ; ierr = ierr + 1 ; Finsi ;
  271.  
  272. 'FINSI' ;
  273.  
  274. Si GRAPH ;
  275.  
  276. $vtf=rxt.'GEO'.'$vtf' ;
  277. mt =doma $vtf maillage ;
  278. Mpl1=chan 'QUAF' plan1 ;
  279. Mpl4=chan 'QUAF' plan4 ;
  280. ELIM (mt et Mpl1 et Mpl4) epsi ;
  281. $mpl1= mode Mpl1 'NAVIER_STOKES' MACRO ;
  282. $mpl4= mode Mpl4 'NAVIER_STOKES' MACRO ;
  283. plan1= doma $mpl1 maillage ;
  284. plan4= doma $mpl4 maillage ;
  285. plan=plan1 et plan4 ;
  286.  
  287. paroif = rxt.'GEO'.'paroif';
  288. rho=rxt.'TIC'.'RHO';
  289. rvp=rxt.'TIC'.'RVP';
  290. tf=rxt.'TIC'.'TF';
  291. un=rxt.'TIC'.'UN' ;
  292. un1=redu un plan ;
  293. ung= vect un1 1. ux uy uz jaune ;
  294. trace ung plan ;
  295. opti isov suli ;
  296. trace rho plan 'TITRE'' Rho' ;
  297. trace rvp plan 'TITRE'' Rvp' ;
  298. trace tf plan 'TITRE'' Tf ' ;
  299.  
  300. trace rho paroif 'TITRE'' Rho' ;
  301. trace rvp paroif 'TITRE'' Rvp' ;
  302. trace Tf paroif 'TITRE'' Tf ' ;
  303. Fcond = tic.'Fcondw';
  304. trace Fcond paroif 'TITRE' ' Flux de condensation Kg / m**2 ' ;
  305.  
  306. axe = p0 d nz1 (p0 plus v1) ;
  307. axe = chan axe 'QUAF' ;
  308. elim (axe et mt) epsi ;
  309.  
  310. evr= evol 'CHPO' rho axe ;
  311. dess evr 'TITRE' ' Rho axe ';
  312.  
  313. evt = evol 'CHPO' Tf axe ;
  314. dess evt 'TITRE' ' Température axe ';
  315.  
  316. evvp= evol 'CHPO' rvp axe ;
  317. dess evvp 'TITRE' ' rvp axe ';
  318.  
  319.  
  320. 'FINSI' ;
  321.  
  322. 'SI' (IERR '>' 0) ;
  323. 'MESS' 'Il y a des problemes !!!' ;
  324. 'ERRE' 5 ;
  325. 'SINO';
  326. 'MESS' 'Tout s est bien passe!' ;
  327. 'FINS' ;
  328.  
  329. 'FIN' ;
  330.  
  331.  
  332.  
  333.  
  334.  
  335.  
  336.  
  337.  
  338.  
  339.  
  340.  
  341.  
  342.  
  343.  
  344.  
  345.  
  346.  
  347.  
  348.  
  349.  
  350.  
  351.  
  352.  

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