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

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