Télécharger enc2dFP.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : enc2dFP.dgibi
  2. *
  3. * Enceinte 2D Axisymetrique
  4. * Relachement d'un melange d'hydrogene et d'helium
  5. * comme avec enc2d.dgibi
  6. * Ici on teste le bon fonctionnement des fonctions de parois
  7. * et de la vitesse normale imposée.
  8. * Testé avec LINE-MSOMMET et MACRO-MSOMMET
  9. * Les couples QUAF-MSOMMET et QUAF-CENTREP1 ne marchent pas
  10.  
  11. 'OPTI' 'DIME' 2 'ELEM' 'CU20' ;
  12. 'OPTI' 'MODE' 'AXIS' ;
  13.  
  14. COMPLET= VRAI ;
  15. COMPLET= FAUX ;
  16. GRAPH = VRAI ;
  17. GRAPH = FAUX ;
  18.  
  19. 'SI' COMPLET ;
  20. nbit=100 ;
  21. nbit=20 ;
  22. DT0 = 1. ;
  23. n1=5 ;
  24. n2=10;
  25. n3=20;
  26. 'SINON' ;
  27. nbit=20 ;
  28. DT0 = 1. ;
  29. n1=5 ;
  30. n2=10;
  31. n3=20;
  32. 'FINSI' ;
  33.  
  34.  
  35. p1=0 0 ;
  36. p11=0.4 0 ;
  37. p2=0.1 0. ;
  38. p3=1 0 ;
  39. p4=1 2;
  40. p5=0 2 ;
  41.  
  42.  
  43. breche = p1 'DROI' n1 p2 ;
  44. bas= p2 'DROI' n2 p3 ;
  45. paroid=p3 'DROI' n3 p4 ;
  46. plaf=p4 'DROI' (n1 + n2) p5 ;
  47. axe = p5 'DROI' n3 p1 ;
  48. vtp = bas 'ET' paroid 'ET' plaf ;
  49. mt= 'DALLER' (breche 'ET' bas) paroid plaf axe ;
  50.  
  51. rxt = 'TABLE' ;
  52. rxt.'VERSION'= 'V0' ;
  53. rxt.'vtf'=mt ;
  54. rxt.'axe' = axe ;
  55. rxt.'breche'=breche ;
  56. rxt.'diru1'=0. 1. ;
  57. rxt.'epsi'= 1.e-4 ;
  58. rxt.'pi' = 0.5 0.5 ;
  59. rxt.'DT0' = DT0 ;
  60. rxt.'DETMAT'= VRAI ;
  61. rxt.'RENU' = 'RIEN' ;
  62. rxt.'DISCR'= 'LINE';
  63. rxt.'DISCR'= 'MACRO';
  64. rxt.'KPRE' = 'MSOMMET';
  65. *rxt.'DISCR'= 'QUAF';
  66. *rxt.'KPRE' = 'CENTREP1';
  67. rxt.'MODTURB' = 'LMEL' ;
  68. rxt.'LMEL' = 0.05 ;
  69. rxt.'FPAROI' = VRAI ;
  70. rxt.'YP'=5.e-3 ;
  71. rxt.'GRAPH' = GRAPH ;
  72. rxt.'HE' = VRAI ;
  73. rxt.'H2' = VRAI ;
  74. rxt.'CO' = VRAI ;
  75.  
  76. rxt.'TF0' = 20.0 ;
  77. rxt.'PT0' = 1.0e5 ;
  78. rxt.'Yhe0' = 0.0 ;
  79. rxt.'Yh20' = 0.0 ;
  80. rxt.'Yco0' = 0.0 ;
  81.  
  82. rxt.'scenario' = table ;
  83. rxt.'scenario'.'t' = prog 0.0 100.0 ;
  84. rxt.'scenario'.'qair' = prog 0.0 0.0 ;
  85. rxt.'scenario'.'qhe' = prog 0.0011 0.0011 ;
  86. rxt.'scenario'.'qh2' = prog 0.0011 0.0011 ;
  87. rxt.'scenario'.'qco' = prog 0.0011 0.0011 ;
  88. rxt.'scenario'.'tinj' = prog 20.0 20.0 ;
  89.  
  90. EXECRXT nbit rxt ;
  91.  
  92. un=rxt.'TIC'.'UN';
  93. une=noel un rxt.'GEO'.'$vtf';
  94. un =elno une rxt.'GEO'.'$vtf';
  95. menvf=rxt.'GEO'.'menvf';
  96. uenv=redu un menvf;
  97. evunv1=evol chpo uenv 'UX' menvf;
  98. evunv2=evol chpo uenv 'UY' menvf;
  99.  
  100. Si GRAPH; dess (evunv1 et evunv2); Finsi;
  101.  
  102. Lunv1= extr (lect 1 pas 10 141) (extr evunv1 'ORDO');
  103. Lunv2= extr (lect 1 pas 10 141) (extr evunv2 'ORDO');
  104.  
  105. list lunv1;
  106. list lunv2;
  107.  
  108.  
  109. LRunv1=Prog
  110. -4.77452E-02 2.58242E-04 3.73809E-03 4.63039E-04 0.18522
  111. 0.86811 0.98377 0.22216 -1.96206E-02 -8.11362E-03
  112. 1.92531E-04 -1.41814E-02 -8.63366E-02 -0.46183 -4.77452E-02;
  113.  
  114. LRunv2=prog
  115. 0.97878 3.3448 4.1101 4.8215 2.4430
  116. 0.14722 1.73137E-02 -0.12874 -0.45017 -0.12690
  117. -5.95941E-02 -3.44642E-02 1.34808E-03 3.06570E-02 0.97878;
  118.  
  119.  
  120. ERPS1=SOMM( abs (LRunv1 - Lunv1) ) ;
  121. ERPS2=SOMM( abs (LRunv2 - Lunv2) ) ;
  122.  
  123. Mess 'ERPS1=' ERPS1 'ERPS2=' ERPS2 ;
  124.  
  125. Si (ERPS1 '>' 5.e-4) ; erreur 5 ; Finsi ;
  126. Si (ERPS2 '>' 5.e-4) ; erreur 5 ; Finsi ;
  127.  
  128.  
  129.  
  130. mess ' Max U ' (maxi un) ' Min U ' (mini un) ;
  131. test = abs ((maxi un) - 1.9 );
  132. mess ' Test =' test;
  133. *Si (test '>' 0.1 ) ; erreur 5 ; Finsi ;
  134.  
  135. tf =rxt.'TIC'.'TF' ;
  136. mess ' Max Tf' (maxi tf) ' Min Tf' (mini tf) ;
  137. test = abs ((maxi tf) - 29.456 );
  138. mess ' Test =' test;
  139. *Si (test '>' 0.5 ) ; erreur 5 ; Finsi ;
  140.  
  141.  
  142.  
  143. $vtf=rxt.'GEO'.'$vtf' ;
  144. vtf=doma $vtf maillage;
  145. rho=rxt.'TIC'.'RHO' ;
  146. rhe=rxt.'TIC'.'RHE';
  147. rh2=rxt.'TIC'.'RH2';
  148. rco=rxt.'TIC'.'RCO';
  149. un =rxt.'TIC'.'UN';
  150.  
  151. axe = chan line axe;
  152. ixe=inve axe ;
  153. evauz = evol chpo (exco un 'UY') ixe ;
  154. evarh = evol chpo rho ixe ;
  155. evatf = evol chpo tf ixe ;
  156. evahe = evol chpo rhe ixe ;
  157. evah2 = evol chpo rh2 ixe ;
  158.  
  159.  
  160. Si GRAPH ;
  161.  
  162. dess evauz 'TITRE' 'Vitesse sur l axe';
  163. dess evarh 'TITRE' 'Densite melange sur l axe';
  164. dess evatf 'TITRE' 'Temperature sur l axe';
  165. dess evahe 'TITRE' 'Densite Helium sur l axe';
  166. dess evah2 'TITRE' 'Densite Hydrogene l axe';
  167.  
  168. ung= vect un 0.1 ux uy jaune;
  169. trace ung vtf;
  170. opti isov suli ;
  171. trace rho vtf 'TITRE' ' Rho' ;
  172. trace rhe vtf 'TITRE' ' Rhe' ;
  173. *trace rh2 vtf 'TITRE' ' Rh2' ;
  174. *trace rco vtf 'TITRE' ' Rco' ;
  175. *trace tf vtf 'TITRE' ' Tf ' ;
  176. Finsi ;
  177.  
  178. FIN ;
  179.  
  180.  
  181.  
  182.  
  183.  
  184.  
  185.  
  186.  
  187.  
  188.  
  189.  

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