Télécharger la plaquette Cast3M


fr en

Télécharger Henc2d.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  Henc2d.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.  
  9.  COMPLET= VRAI                                          ;
  10.  COMPLET= FAUX                                          ;
  11.  GRAPH  = VRAI                                          ;
  12.  GRAPH  = FAUX                                          ;
  13.  
  14.  'SI' COMPLET ;
  15.   nbit=100    ;
  16.   nbit=20     ;
  17.   DT0 = 1.    ;
  18.   n1=5 ;
  19.   n2=10;
  20.   n3=20;
  21.  'SINON'      ;
  22.   nbit=20     ;
  23.   DT0 = 1.    ;
  24.   n1=5 ;
  25.   n2=10;
  26.   n3=20;
  27.  'FINSI'      ;
  28.  
  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. * Mise en place de la procedure
  47. rxt = 'TABLE' ;
  48. rxt . 'VERSION' = 'V0' ;
  49. rxt . 'vtf'     = mt ;
  50. rxt . 'axe'     = axe ;
  51. rxt . 'epsi'    = 1.e-4 ;
  52. *
  53. rxt . 'pi' = 0.5 0.5 ;
  54. *
  55. rxt . 'DISCR'   = 'LINE' ;
  56. rxt . 'KPRE'    = 'MSOMMET' ;
  57. rxt . 'MODTURB' = 'LMEL' ;
  58. rxt . 'LMEL'    = 0.05 ;
  59. *
  60. rxt . 'HE'    = VRAI ;
  61. rxt . 'H2'    = VRAI ;
  62. rxt . 'CO2'   = VRAI ;
  63. rxt . 'Yhe0'  = 0.0 ;
  64. rxt . 'Yh20'  = 0.0 ;
  65. rxt . 'Yco20' = 0.0 ;
  66. *
  67. rxt . 'TF0' = 20.0 ;
  68. rxt . 'PT0' = 1.0e5 ;
  69. *
  70. rxt . 'Breches' = 'TABLE' ;
  71. rxt . 'Breches' . 'A' = 'TABLE' ;
  72. rxt . 'Breches' . 'A' . 'Maillage' = breche ;
  73. rxt . 'Breches' . 'A' . 'diru'  = (0. 1.) ;
  74. rxt . 'Breches' . 'A' . 'scenario'          = 'TABLE' ;
  75. rxt . 'Breches' . 'A' . 'scenario' . 't'    = 'PROG' 0.0     10000.0 ;
  76. rxt . 'Breches' . 'A' . 'scenario' . 'qair' = 'PROG' 0.      0.      ;
  77. rxt . 'Breches' . 'A' . 'scenario' . 'qhe'  = 'PROG' 0.0011  0.0011  ;
  78. rxt . 'Breches' . 'A' . 'scenario' . 'qh2'  = 'PROG' 0.0011  0.0011  ;
  79. rxt . 'Breches' . 'A' . 'scenario' . 'qco2' = 'PROG' 0.0011  0.0011  ;
  80. rxt . 'Breches' . 'A' . 'scenario' . 'tinj' = 'PROG' 20.0    20.0    ;
  81. *
  82. rxt . 'DT0'    = DT0 ;
  83. rxt . 'DETMAT' = VRAI ;
  84. rxt . 'RENU'   = 'RIEN' ;
  85. rxt . 'GRAPH'  = GRAPH ;
  86. *
  87. EXECRXT nbit rxt ;
  88.  
  89.  
  90. un=rxt.'TIC'.'UN';
  91. mess ' Max U ' (maxi un) ' Min U ' (mini un) ;
  92. test = abs ((maxi un) - 1.9    );
  93. mess ' Test =' test;
  94.  Si (test '>' 0.1  ) ; erreur 5 ; Finsi ;
  95.  
  96. tf =rxt.'TIC'.'TF' ;
  97. mess ' Max Tf' (maxi tf) ' Min Tf' (mini tf) ;
  98. test = abs ((maxi tf) - 27.    );
  99. mess ' Test =' test;
  100.  Si (test '>' 0.5  ) ; erreur 5 ; Finsi ;
  101.  
  102.  
  103.  
  104. $vtf=rxt.'GEO'.'$vtf' ;
  105. vtf=doma $vtf maillage;
  106. rho=rxt.'TIC'.'RHO' ;
  107. rhe=rxt.'TIC'.'RHE';
  108. rh2=rxt.'TIC'.'RH2';
  109. rco2=rxt.'TIC'.'RCO2';
  110. un =rxt.'TIC'.'UN';
  111.  
  112. ixe=inve axe ;
  113. evauz = evol chpo (exco un 'UY') ixe ;
  114. evarh = evol chpo rho ixe ;
  115. evatf = evol chpo tf  ixe ;
  116. evahe = evol chpo rhe ixe ;
  117. evah2 = evol chpo rh2 ixe ;
  118.  
  119.  
  120. Si GRAPH ;
  121.  
  122.  dess evauz 'TITRE' 'Vitesse sur l axe';
  123.  dess evarh 'TITRE' 'Densite melange sur l axe';
  124.  dess evatf 'TITRE' 'Temperature sur l axe';
  125.  dess evahe 'TITRE' 'Densite Helium sur l axe';
  126.  dess evah2 'TITRE' 'Densite Hydrogene l axe';
  127.  
  128. ung= vect un 0.1 ux uy jaune;
  129. trace ung vtf;
  130. opti isov suli ;
  131. trace rho vtf 'TITRE' ' Rho' ;
  132. trace rhe vtf 'TITRE' ' Rhe' ;
  133. *trace rh2 vtf 'TITRE' ' Rh2' ;
  134. *trace rco2 vtf 'TITRE' ' Rco2' ;
  135. *trace tf  vtf 'TITRE' ' Tf ' ;
  136. Finsi ;
  137.  
  138.  FIN ;
  139.  
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  
  146.  
  147.  

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