Télécharger la plaquette Cast3M


fr en

Télécharger enc2dpp.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  enc2d.dgibi
  2. *
  3. *  Enceinte 2D Axisymetrique
  4. *  Relachement d'un melange d'hydrogene et d'helium
  5. *
  6. 'OPTI' 'DIME' 2 'ELEM' 'CU20' ;
  7.  'OPTI' 'MODE' 'AXIS' ;
  8.  'OPTI' 'TRAC' 'PSC';
  9.  
  10.  COMPLET= VRAI                                          ;
  11.  COMPLET= FAUX                                          ;
  12.  GRAPH  = VRAI                                          ;
  13.  GRAPH  = FAUX                                          ;
  14.  
  15.  'SI' COMPLET ;
  16.   nbit=100    ;
  17.   nbit=20     ;
  18.   DT0 = 1.    ;
  19.   n1=80 ;
  20.   n2=160;
  21.   n3=320;
  22.  'SINON'      ;
  23.   nbit=20     ;
  24.   DT0 = 1.    ;
  25.   n1=5 ;
  26.   n2=10;
  27.   n3=20;
  28.  'FINSI'      ;
  29.  
  30. p1=0 0 ;
  31. p11=0.4 0 ;
  32. p2=0.1 0. ;
  33. p3=1 0 ;
  34. p4=1 2;                                                                                
  35. p5=0 2 ;
  36.  
  37.  
  38. breche = p1  'DROI' n1 p2 ;
  39. bas= p2 'DROI' n2 p3 ;
  40. paroid=p3 'DROI' n3 p4 ;
  41. plaf=p4 'DROI' (n1 + n2) p5 ;
  42. axe = p5 'DROI' n3 p1 ;
  43. vtp = bas 'ET' paroid 'ET' plaf ;
  44. mt= 'DALLER' (breche 'ET' bas) paroid plaf axe ;
  45.  
  46. Si GRAPH ;
  47. trace mt Titr (chai ' NBEL= ' (nbel mt));
  48. Finsi ;
  49.  
  50. rxt = 'TABLE' ;
  51. rxt.'VERSION'= 'V0'                                           ;
  52. rxt.'vtf'=mt ;
  53. rxt.'axe'   = axe ;
  54. rxt.'breche'=breche ;
  55. rxt.'diru1'=0. 1. ;
  56. rxt.'epsi'= 1.e-4 ;
  57. rxt.'pi'  = 0.5 0.5 ;
  58. rxt.'DT0'  = DT0     ;
  59. rxt.'DETMAT'= VRAI ;
  60.  
  61. rxt.'TKPR'  = VRAI   ;
  62.  
  63. rxt.'RENU'  = 'RIEN' ;
  64. rxt.'DISCR'= 'LINE';
  65. rxt.'KPRE' = 'MSOMMET';
  66. rxt.'MODTURB'  = 'LMEL'  ;
  67. rxt.'LMEL'  = 0.05  ;
  68. rxt.'GRAPH'  = GRAPH ;
  69. rxt.'HE'     = VRAI ;
  70. rxt.'H2'     = VRAI ;
  71. rxt.'CO'     = VRAI ;
  72.  
  73.  rxt.'TF0'    =       20.0                         ;
  74.  rxt.'PT0'    =        1.0e5                       ;
  75.  rxt.'Yhe0'   =        0.0                         ;
  76.  rxt.'Yh20'   =        0.0                         ;
  77.  rxt.'Yco0'   =        0.0                         ;
  78.  
  79.  rxt.'scenario'          = table                              ;
  80.  rxt.'scenario'.'t'      = prog   0.0 100.0                   ;
  81.  rxt.'scenario'.'qair'   = prog   0.0   0.0                   ;
  82.  rxt.'scenario'.'qhe'    = prog   0.0011  0.0011              ;
  83.  rxt.'scenario'.'qh2'    = prog   0.0011  0.0011              ;
  84.  rxt.'scenario'.'qco'    = prog   0.0011  0.0011              ;
  85.  rxt.'scenario'.'tinj'   = prog  20.0  20.0                   ;
  86.  
  87. EXECRXT  nbit rxt ;
  88.  
  89. un=rxt.'TIC'.'UN';
  90. mess ' Max U ' (maxi un) ' Min U ' (mini un) ;
  91. test = abs ((maxi un) - 1.9    );
  92. mess ' Test =' test;
  93.  Si (test '>' 0.1  ) ; erreur 5 ; Finsi ;
  94.  
  95. tf =rxt.'TIC'.'TF' ;
  96. mess ' Max Tf' (maxi tf) ' Min Tf' (mini tf) ;
  97. test = abs ((maxi tf) - 27.    );
  98. mess ' Test =' test;
  99.  Si (test '>' 0.5  ) ; erreur 5 ; Finsi ;
  100.  
  101.  
  102.  
  103. $vtf=rxt.'GEO'.'$vtf' ;
  104. vtf=doma $vtf maillage;
  105. rho=rxt.'TIC'.'RHO' ;
  106. rhe=rxt.'TIC'.'RHE';
  107. rh2=rxt.'TIC'.'RH2';
  108. rco=rxt.'TIC'.'RCO';
  109. rair=rxt.'TIC'.'RAIR';
  110. un =rxt.'TIC'.'UN';
  111.  
  112. DISCR=rxt.'DISCR' ;
  113. axe = chan DISCR axe ;
  114. ixe=inve axe ;
  115. evauz = evol chpo (exco un 'UY') ixe ;
  116. evarh = evol chpo rho ixe ;
  117. evatf = evol chpo tf  ixe ;
  118. evahe = evol chpo rhe ixe ;
  119. evah2 = evol chpo rh2 ixe ;
  120.  
  121.  
  122. Si GRAPH ;
  123.  
  124.  dess evauz 'TITRE' 'Vitesse sur l axe';
  125.  dess evarh 'TITRE' 'Densite melange sur l axe';
  126.  dess evatf 'TITRE' 'Temperature sur l axe';
  127.  dess evahe 'TITRE' 'Densite Helium sur l axe';
  128.  dess evah2 'TITRE' 'Densite Hydrogene l axe';
  129.  
  130. ung= vect un 0.1 ux uy jaune;
  131. trace ung vtf;
  132. opti isov suli ;
  133. trace tf  vtf 'TITRE' ' Tf ' ;
  134. trace rho vtf 'TITRE' ' Rho' ;
  135. trace rhe vtf 'TITRE' ' Rhe' ;
  136. trace rh2 vtf 'TITRE' ' Rh2' ;
  137. trace rco vtf 'TITRE' ' Rco' ;
  138. trace rair vtf 'TITRE' ' Rair' ;
  139. Finsi ;
  140.  
  141. FIN ;
  142.  
  143.  
  144.  
  145.  
  146.  
  147.  
  148.  

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