Télécharger enc2dFP.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : enc2dFP.dgibi similar with enc2d but with wall functions
  2. *
  3. * Adiabatic and 2D axisymetric containment
  4. * H2, He and CO injections
  5. *
  6. * We check MACRO-MSOMMET (instead of LINE-MSOMMET in enc2d).
  7. * With wall functions, QUAF-MSOMMET and QUAF-CENTREP1 are not available.
  8. *
  9. 'OPTI' 'DIME' 2 'ELEM' 'CU20' 'MODE' 'AXIS' 'TRAC' 'PSC' ;
  10. *
  11. COMPLET = FAUX ;
  12. GRAPH = FAUX ;
  13. 'SI' COMPLET ;
  14. nbit = 100 ;
  15. DT0 = 1. ;
  16. n1 = 5 ;
  17. n2 = 10 ;
  18. n3 = 20 ;
  19. 'SINO' ;
  20. nbit = 20 ;
  21. DT0 = 1. ;
  22. n1 = 5 ;
  23. n2 = 10 ;
  24. n3 = 20 ;
  25. 'FINSI' ;
  26. *
  27. * Mesh
  28. *
  29. p1 = 0. 0. ;
  30. p2 = 0.1 0. ;
  31. p3 = 1. 0. ;
  32. p4 = 1. 2. ;
  33. p5 = 0. 2. ;
  34. *
  35. breche = p1 'DROI' n1 p2 ;
  36. bas = p2 'DROI' n2 p3 ;
  37. paroid = p3 'DROI' n3 p4 ;
  38. plaf = p4 'DROI' (n1 + n2) p5 ;
  39. axe = p5 'DROI' n3 p1 ;
  40. vtp = bas 'ET' paroid 'ET' plaf ;
  41. mt = 'DALLER' (breche 'ET' bas) paroid plaf axe ;
  42. *
  43. *- Data for execrxt.procedur
  44. *
  45. rxt = 'TABLE' ;
  46. rxt . 'VERSION' = 'V0' ;
  47. rxt . 'vtf' = mt ;
  48. rxt . 'axe' = axe ;
  49. rxt . 'epsi' = 1.e-4 ;
  50. *
  51. rxt . 'pi' = 0.5 0.5 ;
  52. *
  53. rxt . 'DISCR' = 'MACRO' ;
  54. rxt . 'KPRE' = 'MSOMMET' ;
  55. rxt . 'DT0' = DT0 ;
  56. *
  57. rxt . 'MODTURB' = 'LMEL' ;
  58. rxt . 'LMEL' = 0.05 ;
  59. rxt . 'FPAROI' = VRAI ;
  60. rxt . 'YP' = 5.e-3 ;
  61. *
  62. rxt . 'HE' = VRAI ;
  63. rxt . 'H2' = VRAI ;
  64. rxt . 'CO' = VRAI ;
  65. rxt . 'TF0' = 20.0 ;
  66. rxt . 'PT0' = 1.0e5 ;
  67. rxt . 'Yhe0' = 0.0 ;
  68. rxt . 'Yh20' = 0.0 ;
  69. rxt . 'Yco0' = 0.0 ;
  70. *
  71. rxt . 'Breches' = 'TABLE' ;
  72. rxt . 'Breches' . 'A' = 'TABLE' ;
  73. rxt . 'Breches' . 'A' . 'Maillage' = breche ;
  74. rxt . 'Breches' . 'A' . 'diru' = (0. 1.) ;
  75. rxt . 'Breches' . 'A' . 'scenario' = 'TABLE' ;
  76. rxt . 'Breches' . 'A' . 'scenario' . 't' = 'PROG' 0.0 100.0 ;
  77. rxt . 'Breches' . 'A' . 'scenario' . 'qair' = 'PROG' 0.0 0.0 ;
  78. rxt . 'Breches' . 'A' . 'scenario' . 'qhe' = 'PROG' 0.0011 0.0011 ;
  79. rxt . 'Breches' . 'A' . 'scenario' . 'qh2' = 'PROG' 0.0011 0.0011 ;
  80. rxt . 'Breches' . 'A' . 'scenario' . 'qco' = 'PROG' 0.0011 0.0011 ;
  81. rxt . 'Breches' . 'A' . 'scenario' . 'tinj' = 'PROG' 20.0 20.0 ;
  82. *
  83. rxt . 'GRAPH' = GRAPH ;
  84. rxt . 'DETMAT' = VRAI ;
  85. rxt . 'RENU' = 'RIEN' ;
  86. *
  87. *- Transient
  88. *
  89. EXECRXT nbit rxt ;
  90. *
  91. *- Tests
  92. * FD : Pourquoi le lissage ELNO(NOEL(un)) ?
  93. 'SI' ('NON' COMPLET) ;
  94. ERR1 = 0 ;
  95. un = 'ELNO' ('NOEL' rxt . 'TIC' . 'UN' rxt . 'GEO' . '$vtf')
  96. rxt . 'GEO' . '$vtf' ;
  97. 'MESS' ' Max U ' ('MAXI' un) ' Min U ' ('MINI' un) ;
  98. test = 'ABS' (('MAXI' un) '-' 5.029) ;
  99. 'MESS' ' Test =' test ;
  100. 'SI' (test '>' 0.1 ) ; ERR1 = ERR1 '+' 1 ; 'FINS' ;
  101. *
  102. tf = rxt . 'TIC' . 'TF' ;
  103. 'MESS' ' Max Tf' ('MAXI' tf) ' Min Tf' ('MINI' tf) ;
  104. test = 'ABS' (('MAXI' tf) '-' 27.293) ;
  105. 'MESS' ' Test =' test ;
  106. 'SI' (test '>' 0.5) ; ERR1 = ERR1 '+' 1 ; 'FINS' ;
  107. *
  108. *- Specific test for wall functions
  109. * Each ten points on the wall boundary we check the velocity
  110. * FD : Pourquoi la vitesse normale en paroi n'est-elle pas nulle ?
  111. evunv1 = 'EVOL' 'CHPO' un 'UX' rxt . 'GEO' . 'menvf' ;
  112. evunv2 = 'EVOL' 'CHPO' un 'UY' rxt . 'GEO' . 'menvf' ;
  113. Lunv1 = 'EXTR' ('LECT' 1 'PAS' 10 141) ('EXTR' evunv1 'ORDO') ;
  114. Lunv2 = 'EXTR' ('LECT' 1 'PAS' 10 141) ('EXTR' evunv2 'ORDO') ;
  115. 'LIST' lunv1 ;
  116. 'LIST' lunv2 ;
  117. LRunv1 = 'PROG'
  118. -4.77452E-02 2.58242E-04 3.73809E-03 4.63039E-04 0.18522
  119. 0.86811 0.98377 0.22216 -1.96206E-02 -8.11362E-03
  120. 1.92531E-04 -1.41814E-02 -8.63366E-02 -0.46183 -4.77452E-02 ;
  121. LRunv2 = 'PROG'
  122. 0.97878 3.3448 4.1101 4.8215 2.4430
  123. 0.14722 1.73137E-02 -0.12874 -0.45017 -0.12690
  124. -5.95941E-02 -3.44642E-02 1.34808E-03 3.06570E-02 0.97878 ;
  125. ERPS1 = 'SOMM' ('ABS' (LRunv1 '-' Lunv1)) ;
  126. ERPS2 = 'SOMM' ('ABS' (LRunv2 '-' Lunv2)) ;
  127. 'MESS' 'ERPS1=' ERPS1 ' ' 'ERPS2=' ERPS2 ;
  128. 'SI' (ERPS1 '>' 5.e-4) ; ERR1 = ERR1 '+' 1 ; 'FINS' ;
  129. 'SI' (ERPS2 '>' 5.e-4) ; 'ERRE' 5 ; 'FINS' ;
  130. *
  131. 'SI' ('NEG' ERR1 0) ;
  132. 'ERRE' 5 ;
  133. 'FINS' ;
  134. 'FINS' ;
  135. *
  136. *- Plots devoted to enc2d... cases
  137. *
  138. 'SI' rxt . 'GRAPH' ;
  139. tbt = rxt . 'TBT' ;
  140. tic = rxt . 'TIC' ;
  141. *
  142. $vtf = rxt . 'GEO' . '$vtf' ;
  143. vtf = 'DOMA' $vtf 'MAILLAGE' ;
  144. 'SI' ('EXIS' tic 'TP') ;
  145. $vtp = rxt . 'GEO' . '$vtp' ;
  146. vtp = 'DOMA' $vtp 'MAILLAGE' ;
  147. 'FINS' ;
  148. *
  149. un = tic . 'UN';
  150. ung = 'VECT' un 0.1 'UX' 'UY' 'JAUN' ;
  151. tf = tic . 'TF' ;
  152. rho = tic . 'RHO' ;
  153. rair = tic . 'RAIR' ;
  154. 'SI' tbt . 'THE' ; rhe = tic . 'RHE' ; 'FINS' ;
  155. 'SI' tbt . 'TH2' ; rh2 = tic . 'RH2' ; 'FINS' ;
  156. 'SI' tbt . 'TCO' ; rco = tic . 'RCO' ; 'FINS' ;
  157. 'SI' tbt . 'TCO2' ; rco2 = tic . 'RCO2' ; 'FINS' ;
  158. 'SI' tbt . 'VAPEUR' ; rvap = tic . 'RVAP' ; 'FINS' ;
  159. *
  160. axe = 'CHAN' rxt . 'DISCR' axe ;
  161. ixe = 'INVE' axe ;
  162. *
  163. evauz = 'EVOL' 'CHPO' ('EXCO' un 'UY') ixe ;
  164. 'DESS' evauz
  165. 'TITR' 'Velocity with the z axis' 'MIMA'
  166. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' m/s' ;
  167. *
  168. evatf = 'EVOL' 'CHPO' tf ixe ;
  169. 'DESS' evatf
  170. 'TITR' 'Gas temperature with the z axis' 'MIMA'
  171. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' C' ;
  172. 'TRAC' tf vtf vtp 'TITR' ' Temperature' ;
  173. evarh = 'EVOL' 'CHPO' rho ixe ;
  174. 'DESS' evarh
  175. 'TITR' 'Gas density with the z axis' 'MIMA'
  176. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' kg/m3' ;
  177. 'TRAC' rho vtf ('CHAN' 'LINE' vtp) ung 'TITR' ' Density & velocity' ;
  178. *
  179. 'SI' tbt . 'THE' ;
  180. evahe = 'EVOL' 'CHPO' rhe ixe ;
  181. 'DESS' evahe
  182. 'TITRE' 'Helium density with the z axis' 'MIMA'
  183. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' kg/m3' ;
  184. 'TRAC' rhe vtf vtp 'TITR' ' He density' ;
  185. 'FINS' ;
  186. 'SI' tbt . 'TH2' ;
  187. evah2 = 'EVOL' 'CHPO' rh2 ixe ;
  188. 'DESS' evah2
  189. 'TITR' 'Hydrogen density with the z axis' 'MIMA'
  190. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' kg/m3' ;
  191. 'TRAC' rh2 vtf vtp 'TITR' ' H2 density' ;
  192. 'FINS' ;
  193. 'SI' tbt . 'TCO' ;
  194. evaco = 'EVOL' 'CHPO' rco ixe ;
  195. 'DESS' evaco
  196. 'TITR' 'CO density with the z axis' 'MIMA'
  197. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' kg/m3' ;
  198. 'TRAC' rco vtf vtp 'TITR' ' CO density' ;
  199. 'FINS' ;
  200. 'SI' tbt . 'TCO2' ;
  201. evaco2 = 'EVOL' 'CHPO' rco2 ixe ;
  202. 'DESS' evaco2
  203. 'TITR' 'CO2 density with the z axis' 'MIMA'
  204. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' kg/m3' ;
  205. 'TRAC' rco2 vtf vtp 'TITR' ' CO2 density' ;
  206. 'FINS' ;
  207. 'SI' tbt . 'VAPEUR' ;
  208. evavap = 'EVOL' 'CHPO' rvap ixe ;
  209. 'DESS' evavap
  210. 'TITR' 'Steam density with the z axis' 'MIMA'
  211. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' kg/m3' ;
  212. 'TRAC' rvap vtf vtp 'TITR' ' Steam density' ;
  213. 'FINS' ;
  214. *
  215. * Wall temperature or air density (in cases enc2d with wall temperature
  216. * only air is involved ; hence rair=rho and rho is plotted above))
  217. 'SI' ('EXIS' tic 'TP') ;
  218. 'TRAC' tic . 'TP' vtp 'TITR' ' Wall temperature' ;
  219. 'SINO' ;
  220. evarair = 'EVOL' 'CHPO' rair ixe ;
  221. 'DESS' evarair
  222. 'TITR' 'Air density with the z axis' 'MIMA'
  223. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' kg/m3' ;
  224. 'TRAC' rair vtf vtp 'TITR' ' Air density' ;
  225. 'FINS' ;
  226. *
  227. 'SI' ('EXIS' tic 'KN') ;
  228. 'TRAC' tic . 'KN' vtf vtp 'TITR' ' KN' ;
  229. 'TRAC' tic . 'EN' vtf vtp 'TITR' ' EN' ;
  230. 'FINS' ;
  231. *
  232. 'FINS' ;
  233. *
  234. * Specific plots for wall functions
  235. 'SI' rxt . 'GRAPH' ;
  236. TAB1 = 'TABLE' ;
  237. TAB1 . 'TITRE' = TABLE ;
  238. TAB1 . 'TITRE' . 1 = 'MOT' ' UX' ;
  239. TAB1 . 'TITRE' . 2 = 'MOT' ' UY' ;
  240. TAB1 . 1 = ' REGU MARQ TRIL' ;
  241. TAB1 . 2 = ' REGU MARQ TRIU' ;
  242. 'DESS' (evunv1 'ET' evunv2)
  243. 'TITR' 'Velocity with the wall curvilinear abscissa' 'MIMA'
  244. 'GRIL' 'POIN' 'GRIS' 'TITX' 's' 'TITY' ' m/s' 'LEGE' tab1 ;
  245. 'FINS' ;
  246. *
  247. FIN ;
  248.  
  249.  
  250.  
  251.  
  252.  
  253.  
  254.  
  255.  
  256.  
  257.  
  258.  
  259.  
  260.  

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