Télécharger la plaquette Cast3M


fr en

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