Télécharger pressu.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : pressu.dgibi
  2. *
  3. * Containment pressurization (Phebus size)
  4. *
  5. * 3D mesh of a 14.5 m3 cylindrical containment with a 10cm in depth vertical wall.
  6. * Initial pressure and temperature are 1bar and 40oC with the wall at 60oC.
  7. * Steam is injected at a constant 50g/s mass flow rate and a temperature of 150oC.
  8. *
  9. * Transient is computed for 50s. We check the pressure, the averaged temperature,
  10. * the velocity and the condensation mass flow rate (that starts at 20s) evolutions.
  11. *
  12. * This test case is similar as pressu2.dgibi except for the mass conservation
  13. * of the scalar fields : we use the new paradigm (see 'CORLIM' and 'CORTEMP')
  14. *
  15. * Auteurs : E. Studer, J.P. Magnaud Novembre 1999
  16. *
  17. 'OPTI' 'DIME' 3 'ELEM' 'CU20' 'TRAC' 'PSC' ;
  18. 'DENS' 1. ;
  19. *
  20. COMPLET = FAUX ;
  21. GRAPH = FAUX ;
  22. 'SI' COMPLET ;
  23. nbit = 100 ;
  24. DT0 = 1. ;
  25. n1 = 1 ;
  26. n2 = 4 ;
  27. n3 = 4 ;
  28. nn = 2 ;
  29. 'SINO' ;
  30. nbit = 5 ;
  31. DT0 = 10. ;
  32. n1 = 1 ;
  33. n2 = 2 ;
  34. n3 = 4 ;
  35. nn = 1 ;
  36. 'FINS' ;
  37. epsi = 1.D-2 ;
  38. *
  39. *
  40. * Mesh
  41. *
  42. *
  43. * ri=cavity radius, h1=cavity height, sp=wall depth,
  44. * fg1=break to cavity radius ratio
  45. ri = 1.052 ;
  46. h1 = 4.163 ;
  47. sp = 0.10 ;
  48. fg1 = 0.25 ;
  49. fg2 = fg1 * (2.0 ** 0.5) / 2. ;
  50. *
  51. p0 = 0.000 0.000 0.000 ;
  52. p1 = (ri*fg1) 0.000 0.000 ;
  53. p2 = (ri*fg2) (ri*fg2) 0.000 ;
  54. p3 = 0.000 (ri*fg1) 0.000 ;
  55. p4 = ri 0.000 0.000 ;
  56. p5 = 0.000 ri 0.000 ;
  57. p6 = (ri+sp) 0.000 0.000 ;
  58. p7 = 0.000 (ri+sp) 0.000 ;
  59. *
  60. * Fluid and structure basement surfaces
  61. * (built by symetry)
  62. l1 = 'DROI' p0 p1 n1 ;
  63. l2 = 'DROI' p1 p2 n1 ;
  64. l3 = 'DROI' p2 p3 n1 ;
  65. l4 = 'DROI' p3 p0 n1 ;
  66. l5 = 'CERC' p4 p0 p5 (2*n1) ;
  67. l6 = 'CERC' p6 p0 p7 (2*n1) ;
  68. basf0 = 'DALL' l1 l2 l3 l4 'PLAN' ;
  69. basf1 = ('REGL' (l2 'ET' l3) l5 n2) ;
  70. *
  71. l44 = 'COTE' 2 basf1;
  72. ax4 = ('INVE' l4) 'ET' l44 ;
  73. l11 = 'COTE' 4 basf1;
  74. ax1 = l11 'ET' ('INVE' l1) ;
  75. basf = basf0 'ET' ('REGL' (l2 'ET' l3) l5 n2) ;
  76. 'ELIM' basf epsi ;
  77. basf = basf 'ET' ('SYME' basf 'DROI' p0 p3) ;
  78. ax11 = ('SYME' ax1 'DROI' p0 p3) 'ET' ('INVE' ax1) ;
  79. 'ELIM' basf epsi ;
  80. basf = basf 'ET' ('SYME' basf 'DROI' p0 p1) ;
  81. ax44 = ('INVE' ax4) 'ET' ('SYME' ax4 'DROI' p0 p4) ;
  82. 'ELIM' basf epsi ;
  83. basm = 'REGL' l5 l6 n3 ;
  84. basm = basm 'ET' ('SYME' basm 'DROI' p0 p3) ;
  85. 'ELIM' basm epsi ;
  86. basm = basm 'ET' ('SYME' basm 'DROI' p0 p1) ;
  87. 'ELIM' basm epsi ;
  88. *
  89. * Fluid and structure volumes
  90. * (built by translation)
  91. nz1 = ('ENTI' (h1 '/' (ri '/' 2.))) '*' nn ;
  92. v1 = 0. 0. h1 ;
  93. mt = basf 'VOLU' nz1 'TRAN' v1 ;
  94. wall = (basm 'VOLU' nz1 'TRAN' v1) 'COUL' 'ROUG' ;
  95. plan1 = ax11 'TRAN' nz1 v1 ;
  96. plan4 = ax44 'TRAN' nz1 v1 ;
  97. 'ELIM' (mt et wall 'ET' plan1 'ET' plan4) epsi ;
  98. *
  99. * Break at the basement if any
  100. pjg = 'POIN' basf 'PROC' (0. 0. 0.) ;
  101. breche = ('ELEM' basf 'APPUIE' 'LARGEMENT' pjg) 'COUL' 'VERT' ;
  102. *
  103. *
  104. * Data for execrxt.procedur
  105. *
  106. *
  107. rxt = 'TABLE' ;
  108. rxt . 'VERSION' = 'V0' ;
  109. rxt . 'vtf' = mt ;
  110. rxt . 'epsi' = epsi ;
  111. rxt . 'pi' = 0. 0. 0.5 ;
  112. *
  113. rxt . 'DISCR' = 'MACRO';
  114. rxt . 'KPRE' = 'CENTRE';
  115. rxt . 'DT0' = DT0 ;
  116. *
  117. * New paradigm (V2022)
  118. rxt . 'CORLIM' = faux ;
  119. rxt . 'CORTEMP' = faux ;
  120. *
  121. rxt . 'MODTURB' = 'NUTURB' ;
  122. rxt . 'NUT' = 1.D-2 ;
  123. *
  124. rxt . 'THERMP' = VRAI ;
  125. rxt . 'vtp' = wall ;
  126. rxt . 'LAMBDA' = 15. ;
  127. rxt . 'ROCP' = 3.9E6 ;
  128. rxt . 'Tp0' = 60. ;
  129. rxt . 'ECHAN' = 10. ;
  130. *
  131. rxt . 'VAPEUR' = VRAI ;
  132. rxt . 'TF0' = 40.0 ;
  133. rxt . 'PT0' = 1.0e5 ;
  134. rxt . 'Yvap0' = 0.0023 ;
  135. *
  136. rxt . 'Breches' = 'TABLE' ;
  137. rxt . 'Breches' . 'A' = 'TABLE' ;
  138. rxt . 'Breches' . 'A' . 'scenario' = 'TABLE' ;
  139. rxt . 'Breches' . 'A' . 'Maillage' = breche ;
  140. rxt . 'Breches' . 'A' . 'diru' = (0. 0. 1.) ;
  141. rxt . 'Breches' . 'A' . 'scenario' . 't' = 'PROG' 0.0 1000.0 ;
  142. rxt . 'Breches' . 'A' . 'scenario' . 'qeau' = 'PROG' 0.050 0.050 ;
  143. rxt . 'Breches' . 'A' . 'scenario' . 'qair' = 'PROG' 0.000 0.000 ;
  144. rxt . 'Breches' . 'A' . 'scenario' . 'tinj' = 'PROG' 150.0 150.0 ;
  145. *
  146. * Ancienne syntaxe
  147. * rxt . 'breche' = breche ;
  148. * rxt . 'diru1' = 0. 0. 1. ;
  149. * rxt . 'scenario' = 'TABLE' ;
  150. * rxt . 'scenario' . 't' = 'PROG' 0.0 1000.0 ;
  151. * rxt . 'scenario' . 'qeau' = 'PROG' 0.050 0.050 ;
  152. * rxt . 'scenario' . 'qair' = 'PROG' 0.000 0.000 ;
  153. * rxt . 'scenario' . 'tinj' = 'PROG' 150.0 150.0 ;
  154. *
  155. rxt . 'GRAPH' = GRAPH ;
  156. rxt . 'DETMAT' = VRAI ;
  157. rxt . 'FRPREC' = 5 ;
  158. rxt . 'RENU' = 'RIEN' ;
  159. *
  160. *
  161. * Transient (with restart after 2 time steps)
  162. *
  163. *
  164. EXECRXT 2 rxt ;
  165. EXECRXT (nbit - 2) rxt ;
  166. *
  167. *
  168. * Tests
  169. *
  170. *
  171. 'SI' ('NON' COMPLET) ;
  172. ERR1 = 0 ;
  173. 'LIST' rxt.TIC.'Tfm' ;
  174. 'LIST' rxt.TIC.'PT' ;
  175. 'LIST' rxt.TIC.'Qc' ;
  176. 'LIST' rxt.TIC.'LMAXU';
  177.  
  178. ltfm = 'PROG'
  179. 40.000 65.434 77.251 86.108 92.289 95.977 ;
  180. lPT = 'PROG'
  181. 1.00000E+05 1.06009E+05 1.16438E+05 1.26781E+05 1.36502E+05 1.45161E+05;
  182. Lqc = 'PROG'
  183. 0.0000 0.0000 0.0000 3.12085E-04 2.66573E-03 5.14421E-03;
  184. Lmaxu = 'PROG'
  185. 0.0000 0.81320 2.0813 2.9715 2.2448 2.5356 ;
  186.  
  187. tic = rxt . 'TIC' ;
  188. ERtf = 'SOMM' ('ABS'(ltfm - tic . 'Tfm')) '/' 80. ;
  189. ERPT = 'SOMM' ('ABS'(lPT - tic . 'PT' )) '/' 1.e5 ;
  190. ERQc = 'SOMM' ('ABS'(lqc - tic . 'Qc' )) ;
  191. ERum = 'SOMM' ('ABS'(Lmaxu - tic . 'LMAXU')) '/' 2. ;
  192. 'MESS' 'ERtf=' ERtf ' ' 'ERPT=' ERPT ' '
  193. 'ERQc=' ERQc ' ' 'ERum=' ERum ;
  194. 'SI' (ERtf '>' 1.e-4) ; err1 = err1 '+' 1 ; 'FINS' ;
  195. 'SI' (ERPT '>' 1.e-3) ; err1 = err1 '+' 1 ; 'FINS' ;
  196. 'SI' (ERQc '>' 1.e-4) ; err1 = err1 '+' 1 ; 'FINS' ;
  197. 'SI' (ERum '>' 1.e-2) ; err1 = err1 '+' 1 ; 'FINS' ;
  198. 'SI' ('NEG' ERR1 0) ;
  199. 'ERRE' 5 ;
  200. 'FINS' ;
  201. 'FINS' ;
  202. *
  203. *
  204. * Plots devoted to pressu... cases
  205. *
  206. *
  207. 'SI' GRAPH ;
  208. tbt = rxt . 'TBT' ;
  209. tic = rxt . 'TIC' ;
  210. *
  211. $vtf = rxt . 'GEO' . '$vtf' ;
  212. vtf = 'DOMA' $vtf 'MAILLAGE' ;
  213. *
  214. Mpl1 = 'CHAN' 'QUAF' plan1 ;
  215. Mpl4 = 'CHAN' 'QUAF' plan4 ;
  216. 'ELIM' (vtf 'ET' Mpl1 'ET' Mpl4) epsi ;
  217. $mpl1 = 'MODE' Mpl1 'NAVIER_STOKES' 'MACRO' ;
  218. $mpl4 = 'MODE' Mpl4 'NAVIER_STOKES' 'MACRO' ;
  219. plan1 = 'DOMA' $mpl1 'MAILLAGE' ;
  220. plan4 = 'DOMA' $mpl4 'MAILLAGE' ;
  221. plan = plan1 'ET' plan4 ;
  222. cplan = 'CONT' plan ;
  223. *
  224. 'SI' ('EXIS' tic 'TP') ;
  225. $vtp = rxt . 'GEO' . '$vtp' ;
  226. vtp = 'DOMA' $vtp 'MAILLAGE' ;
  227. 'FINS' ;
  228. paroif = rxt . 'GEO' . 'paroif';
  229. cparoif = 'CONT' paroif ;
  230. *
  231. axe = p0 d nz1 (p0 plus v1) ;
  232. axe = 'CHAN' axe 'QUAF' ;
  233. 'ELIM' (axe 'ET' mt) epsi ;
  234. *
  235. un = tic . 'UN';
  236. unp = 'REDU' un plan ;
  237. ung = 'VECT' un 0.5 'UX' 'UY' 'UZ' 'JAUN' ;
  238. ungp = 'VECT' unp 0.5 'UX' 'UY' 'UZ' 'JAUN' ;
  239. tf = tic . 'TF' ;
  240. rho = tic . 'RHO' ;
  241. rair = tic . 'RAIR' ;
  242. 'SI' tbt . 'THE' ; rhe = tic . 'RHE' ; 'FINS' ;
  243. 'SI' tbt . 'TH2' ; rh2 = tic . 'RH2' ; 'FINS' ;
  244. 'SI' tbt . 'TCO' ; rco = tic . 'RCO' ; 'FINS' ;
  245. 'SI' tbt . 'TCO2' ; rco2 = tic . 'RCO2' ; 'FINS' ;
  246. 'SI' tbt . 'VAPEUR' ; rvp = tic . 'RVP' ; 'FINS' ;
  247. *
  248. evauz = 'EVOL' 'CHPO' ('EXCO' un 'UZ') axe ;
  249. 'DESS' evauz
  250. 'TITR' 'Velocity with the z axis' 'MIMA'
  251. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' m/s' ;
  252. 'TRAC' ung plan ('CONT' plan) 'TITR' ' Velocity' ;
  253. *
  254. evatf = 'EVOL' 'CHPO' tf axe ;
  255. 'DESS' evatf
  256. 'TITR' 'Gas temperature with the z axis' 'MIMA'
  257. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' C' ;
  258. 'TRAC' tf plan cplan 'TITR' ' Temperature' ;
  259. 'TRAC' tf paroif cparoif 'TITR' ' Temperature' ;
  260. *
  261. evarh = 'EVOL' 'CHPO' rho axe ;
  262. 'DESS' evarh
  263. 'TITR' 'Gas density with the z axis' 'MIMA'
  264. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' kg/m3' ;
  265. 'TRAC' rho plan ('CONT' plan) ungp 'TITR' ' Density & velocity' ;
  266. 'TRAC' rho paroif cparoif 'TITR' ' Density' ;
  267. *
  268. evavap = 'EVOL' 'CHPO' rvp axe ;
  269. 'DESS' evavap
  270. 'TITR' 'Steam density with the z axis' 'MIMA'
  271. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' kg/m3' ;
  272. 'TRAC' rvp plan cplan 'TITR' ' Steam density' ;
  273. *
  274. Fcond = rxt . 'TIC' . 'Fcondw';
  275. 'TRAC' fcond paroif cparoif 'TITR' ' Fcond kg/m**2' ;
  276. *
  277. * Wall temperature
  278. 'SI' ('EXIS' tic 'TP') ;
  279. 'TRAC' tic . 'TP' vtp 'TITR' ' Wall temperature' ;
  280. 'FINS' ;
  281. 'FINSI' ;
  282.  
  283. 'FIN' ;
  284.  
  285.  
  286.  
  287.  
  288.  
  289.  
  290.  
  291.  
  292.  
  293.  
  294.  
  295.  
  296.  
  297.  
  298.  
  299.  
  300.  
  301.  
  302.  
  303.  
  304.  

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