Télécharger enc2dQ.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : enc2dQ.dgibi (similar to enc2d except quadratic elements QUAF)
  2. *
  3. * Adiabatic and 2D axisymetric containment
  4. * H2, He and CO2 injections
  5. *
  6. 'OPTI' 'DIME' 2 'ELEM' 'CU20' 'MODE' 'AXIS' 'TRAC' 'PSC' 'DENS' 1. ;
  7. *
  8. COMPLET = FAUX ;
  9. GRAPH = FAUX ;
  10. 'SI' COMPLET ;
  11. nbit = 100 ;
  12. DT0 = 1. ;
  13. n1 = 5 ;
  14. n2 = 10 ;
  15. n3 = 20 ;
  16. 'SINO' ;
  17. nbit = 40 ;
  18. DT0 = 0.1 ;
  19. n1 = 3 ;
  20. n2 = 5 ;
  21. n3 = 10 ;
  22. 'FINSI' ;
  23. *
  24. * Mesh
  25. *
  26. p1 = 5.D-4 0. ;
  27. p2 = 0.1 0. ;
  28. p3 = 1. 0. ;
  29. p4 = 1. 2. ;
  30. p5 = 5.D-4 2. ;
  31. *
  32. breche = p1 'DROI' n1 p2 ;
  33. bas = p2 'DROI' n2 p3 ;
  34. paroid = p3 'DROI' n3 p4 ;
  35. plaf = p4 'DROI' (n1 + n2) p5 ;
  36. axe = p5 'DROI' n3 p1 ;
  37. vtp = bas 'ET' paroid 'ET' plaf ;
  38. mt = 'DALLER' (breche 'ET' bas) paroid plaf axe ;
  39. *
  40. *- Data for execrxt.procedur
  41. *
  42. rxt = 'TABLE' ;
  43. rxt . 'VERSION' = 'V0' ;
  44. rxt . 'vtf' = mt ;
  45. rxt . 'axe' = axe ;
  46. rxt . 'epsi' = 1.e-4 ;
  47. *
  48. rxt . 'pi' = 0.5 0.5 ;
  49. *
  50. rxt . 'DISCR' = 'QUAF' ;
  51. rxt . 'KPRE' = 'MSOMMET' ;
  52. rxt . 'DT0' = DT0 ;
  53. rxt . 'CORTEMP' = faux ;
  54. *
  55. rxt . 'MODTURB' = 'LMEL' ;
  56. rxt . 'LMEL' = 0.05 ;
  57. *
  58. rxt . 'HE' = VRAI ;
  59. rxt . 'H2' = VRAI ;
  60. rxt . 'CO2' = VRAI ;
  61. rxt . 'TF0' = 20.0 ;
  62. rxt . 'PT0' = 1.0e5 ;
  63. rxt . 'Yhe0' = 0.0 ;
  64. rxt . 'Yh20' = 0.0 ;
  65. rxt . 'Yco20' = 0.0 ;
  66. *
  67. rxt . 'Breches' = 'TABLE' ;
  68. rxt . 'Breches' . 'A' = 'TABLE' ;
  69. rxt . 'Breches' . 'A' . 'Maillage' = breche ;
  70. rxt . 'Breches' . 'A' . 'diru' = (0. 1.) ;
  71. rxt . 'Breches' . 'A' . 'scenario' = 'TABLE' ;
  72. rxt . 'Breches' . 'A' . 'scenario' . 't' = 'PROG' 0.0 100.0 ;
  73. rxt . 'Breches' . 'A' . 'scenario' . 'qair' = 'PROG' 0.0 0.0 ;
  74. rxt . 'Breches' . 'A' . 'scenario' . 'qhe' = 'PROG' 0.0011 0.0011 ;
  75. rxt . 'Breches' . 'A' . 'scenario' . 'qh2' = 'PROG' 0.0011 0.0011 ;
  76. rxt . 'Breches' . 'A' . 'scenario' . 'qco2' = 'PROG' 0.0011 0.0011 ;
  77. rxt . 'Breches' . 'A' . 'scenario' . 'tinj' = 'PROG' 20.0 20.0 ;
  78. *
  79. rxt . 'GRAPH' = GRAPH ;
  80. rxt . 'DETMAT' = VRAI ;
  81. rxt . 'RENU' = 'RIEN' ;
  82. *
  83. *- Transient
  84. *
  85. EXECRXT nbit rxt ;
  86. *
  87. *- Tests
  88. *
  89. * probleme en 32 bits avec gfortran 14.1
  90. 'SI' ('NON' COMPLET) ;
  91. ERR1 = 0 ;
  92. 'SI' (2**31 > 0) ;
  93. un = rxt . 'TIC' . 'UN';
  94. 'MESS' ' Max U ' ('MAXI' un) ' Min U ' ('MINI' un) ;
  95. test = 'ABS' (('MAXI' un) - 2.1261) ;
  96. 'MESS' ' Test =' test ;
  97. 'SI' (test '>' 0.6 ) ; ERR1 = ERR1 '+' 1 ; 'FINS' ;
  98. 'FINSI';
  99. *
  100. tf = rxt . 'TIC' . 'TF' ;
  101. 'MESS' ' Max Tf' (maxi tf) ' Min Tf' (mini tf) ;
  102. test = abs ((maxi tf) - 20.972);
  103. 'MESS' ' Test =' test;
  104. 'SI' (test '>' 1.0) ; ERR1 = ERR1 '+' 1 ; 'FINS' ;
  105. list test;
  106. list (test '>' 1.0);
  107. list err1;
  108. *
  109. 'SI' ('NEG' ERR1 0) ;
  110. 'ERRE' 5 ;
  111. 'FINS' ;
  112. 'FINS' ;
  113. *
  114. *- Plots devoted to enc2d... cases
  115. *
  116. 'SI' rxt . 'GRAPH' ;
  117. tbt = rxt . 'TBT' ;
  118. tic = rxt . 'TIC' ;
  119. *
  120. $vtf = rxt . 'GEO' . '$vtf' ;
  121. vtf = 'DOMA' $vtf 'MAILLAGE' ;
  122. 'SI' ('EXIS' tic 'TP') ;
  123. $vtp = rxt . 'GEO' . '$vtp' ;
  124. vtp = 'DOMA' $vtp 'MAILLAGE' ;
  125. 'FINS' ;
  126. *
  127. un = tic . 'UN';
  128. ung = 'VECT' un 0.1 'UX' 'UY' 'JAUN' ;
  129. tf = tic . 'TF' ;
  130. rho = tic . 'RHO' ;
  131. rair = tic . 'RAIR' ;
  132. 'SI' tbt . 'THE' ; rhe = tic . 'RHE' ; 'FINS' ;
  133. 'SI' tbt . 'TH2' ; rh2 = tic . 'RH2' ; 'FINS' ;
  134. 'SI' tbt . 'TCO' ; rco = tic . 'RCO' ; 'FINS' ;
  135. 'SI' tbt . 'TCO2' ; rco2 = tic . 'RCO2' ; 'FINS' ;
  136. 'SI' tbt . 'VAPEUR' ; rvap = tic . 'RVAP' ; 'FINS' ;
  137. *
  138. axe = 'CHAN' rxt . 'DISCR' axe ;
  139. ixe = 'INVE' axe ;
  140. *
  141. evauz = 'EVOL' 'CHPO' ('EXCO' un 'UY') ixe ;
  142. 'DESS' evauz
  143. 'TITR' 'Velocity with the z axis' 'MIMA'
  144. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' m/s' ;
  145. *
  146. evatf = 'EVOL' 'CHPO' tf ixe ;
  147. 'DESS' evatf
  148. 'TITR' 'Gas temperature with the z axis' 'MIMA'
  149. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' C' ;
  150. 'TRAC' tf vtf vtp 'TITR' ' Temperature' ;
  151. evarh = 'EVOL' 'CHPO' rho ixe ;
  152. 'DESS' evarh
  153. 'TITR' 'Gas density with the z axis' 'MIMA'
  154. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' kg/m3' ;
  155. 'TRAC' rho vtf ('CHAN' 'LINE' vtp) ung 'TITR' ' Density & velocity' ;
  156. *
  157. 'SI' tbt . 'THE' ;
  158. evahe = 'EVOL' 'CHPO' rhe ixe ;
  159. 'DESS' evahe
  160. 'TITRE' 'Helium density with the z axis' 'MIMA'
  161. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' kg/m3' ;
  162. 'TRAC' rhe vtf vtp 'TITR' ' He density' ;
  163. 'FINS' ;
  164. 'SI' tbt . 'TH2' ;
  165. evah2 = 'EVOL' 'CHPO' rh2 ixe ;
  166. 'DESS' evah2
  167. 'TITR' 'Hydrogen density with the z axis' 'MIMA'
  168. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' kg/m3' ;
  169. 'TRAC' rh2 vtf vtp 'TITR' ' H2 density' ;
  170. 'FINS' ;
  171. 'SI' tbt . 'TCO' ;
  172. evaco = 'EVOL' 'CHPO' rco ixe ;
  173. 'DESS' evaco
  174. 'TITR' 'CO density with the z axis' 'MIMA'
  175. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' kg/m3' ;
  176. 'TRAC' rco vtf vtp 'TITR' ' CO density' ;
  177. 'FINS' ;
  178. 'SI' tbt . 'TCO2' ;
  179. evaco2 = 'EVOL' 'CHPO' rco2 ixe ;
  180. 'DESS' evaco2
  181. 'TITR' 'CO2 density with the z axis' 'MIMA'
  182. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' kg/m3' ;
  183. 'TRAC' rco2 vtf vtp 'TITR' ' CO2 density' ;
  184. 'FINS' ;
  185. 'SI' tbt . 'VAPEUR' ;
  186. evavap = 'EVOL' 'CHPO' rvap ixe ;
  187. 'DESS' evavap
  188. 'TITR' 'Steam density with the z axis' 'MIMA'
  189. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' kg/m3' ;
  190. 'TRAC' rvap vtf vtp 'TITR' ' Steam density' ;
  191. 'FINS' ;
  192. *
  193. * Wall temperature or air density (in cases enc2d with wall temperature
  194. * only air is involved ; hence rair=rho and rho is plotted above))
  195. 'SI' ('EXIS' tic 'TP') ;
  196. 'TRAC' tic . 'TP' vtp 'TITR' ' Wall temperature' ;
  197. 'SINO' ;
  198. evarair = 'EVOL' 'CHPO' rair ixe ;
  199. 'DESS' evarair
  200. 'TITR' 'Air density with the z axis' 'MIMA'
  201. 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' kg/m3' ;
  202. 'TRAC' rair vtf vtp 'TITR' ' Air density' ;
  203. 'FINS' ;
  204. *
  205. 'SI' ('EXIS' tic 'KN') ;
  206. 'TRAC' tic . 'KN' vtf vtp 'TITR' ' KN' ;
  207. 'TRAC' tic . 'EN' vtf vtp 'TITR' ' EN' ;
  208. 'FINS' ;
  209. 'FINS' ;
  210. *
  211. FIN ;
  212.  
  213.  
  214.  
  215.  
  216.  
  217.  
  218.  
  219.  
  220.  
  221.  
  222.  
  223.  
  224.  
  225.  
  226.  
  227.  
  228.  

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