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.'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. rxt.'RENU' = 'RIEN' ;
  61. rxt.'DISCR'= 'LINE';
  62. rxt.'DISCR'= 'MACRO';
  63. rxt.'KPRE' = 'MSOMMET';
  64. *rxt.'DISCR'= 'QUAF';
  65. *rxt.'KPRE' = 'CENTREP1';
  66. rxt.'MODTURB' = 'LMEL' ;
  67. rxt.'LMEL' = 0.05 ;
  68. rxt.'FPAROI' = VRAI ;
  69. rxt.'YP'=5.e-3 ;
  70. rxt.'GRAPH' = GRAPH ;
  71. rxt.'HE' = VRAI ;
  72. rxt.'H2' = VRAI ;
  73. rxt.'CO' = VRAI ;
  74.  
  75. rxt.'TF0' = 20.0 ;
  76. rxt.'PT0' = 1.0e5 ;
  77. rxt.'Yhe0' = 0.0 ;
  78. rxt.'Yh20' = 0.0 ;
  79. rxt.'Yco0' = 0.0 ;
  80.  
  81. rxt.'scenario' = table ;
  82. rxt.'scenario'.'t' = prog 0.0 100.0 ;
  83. rxt.'scenario'.'qair' = prog 0.0 0.0 ;
  84. rxt.'scenario'.'qhe' = prog 0.0011 0.0011 ;
  85. rxt.'scenario'.'qh2' = prog 0.0011 0.0011 ;
  86. rxt.'scenario'.'qco' = prog 0.0011 0.0011 ;
  87. rxt.'scenario'.'tinj' = prog 20.0 20.0 ;
  88.  
  89. EXECRXT nbit rxt ;
  90.  
  91. un=rxt.'TIC'.'UN';
  92. une=noel un rxt.'GEO'.'$vtf';
  93. un =elno une rxt.'GEO'.'$vtf';
  94. menvf=rxt.'GEO'.'menvf';
  95. uenv=redu un menvf;
  96. evunv1=evol chpo uenv 'UX' menvf;
  97. evunv2=evol chpo uenv 'UY' menvf;
  98.  
  99. Si GRAPH; dess (evunv1 et evunv2); Finsi;
  100.  
  101. Lunv1= extr (lect 1 pas 10 141) (extr evunv1 'ORDO');
  102. Lunv2= extr (lect 1 pas 10 141) (extr evunv2 'ORDO');
  103.  
  104. list lunv1;
  105. list lunv2;
  106.  
  107.  
  108. LRunv1=Prog
  109. -4.77452E-02 2.58242E-04 3.73809E-03 4.63039E-04 0.18522
  110. 0.86811 0.98377 0.22216 -1.96206E-02 -8.11362E-03
  111. 1.92531E-04 -1.41814E-02 -8.63366E-02 -0.46183 -4.77452E-02;
  112.  
  113. LRunv2=prog
  114. 0.97878 3.3448 4.1101 4.8215 2.4430
  115. 0.14722 1.73137E-02 -0.12874 -0.45017 -0.12690
  116. -5.95941E-02 -3.44642E-02 1.34808E-03 3.06570E-02 0.97878;
  117.  
  118.  
  119. ERPS1=SOMM( abs (LRunv1 - Lunv1) ) ;
  120. ERPS2=SOMM( abs (LRunv2 - Lunv2) ) ;
  121.  
  122. Mess 'ERPS1=' ERPS1 'ERPS2=' ERPS2 ;
  123.  
  124. Si (ERPS1 '>' 5.e-4) ; erreur 5 ; Finsi ;
  125. Si (ERPS2 '>' 5.e-4) ; erreur 5 ; Finsi ;
  126.  
  127.  
  128.  
  129. mess ' Max U ' (maxi un) ' Min U ' (mini un) ;
  130. test = abs ((maxi un) - 1.9 );
  131. mess ' Test =' test;
  132. *Si (test '>' 0.1 ) ; erreur 5 ; Finsi ;
  133.  
  134. tf =rxt.'TIC'.'TF' ;
  135. mess ' Max Tf' (maxi tf) ' Min Tf' (mini tf) ;
  136. test = abs ((maxi tf) - 29.456 );
  137. mess ' Test =' test;
  138. *Si (test '>' 0.5 ) ; erreur 5 ; Finsi ;
  139.  
  140.  
  141.  
  142. $vtf=rxt.'GEO'.'$vtf' ;
  143. vtf=doma $vtf maillage;
  144. rho=rxt.'TIC'.'RHO' ;
  145. rhe=rxt.'TIC'.'RHE';
  146. rh2=rxt.'TIC'.'RH2';
  147. rco=rxt.'TIC'.'RCO';
  148. un =rxt.'TIC'.'UN';
  149.  
  150. axe = chan line axe;
  151. ixe=inve axe ;
  152. evauz = evol chpo (exco un 'UY') ixe ;
  153. evarh = evol chpo rho ixe ;
  154. evatf = evol chpo tf ixe ;
  155. evahe = evol chpo rhe ixe ;
  156. evah2 = evol chpo rh2 ixe ;
  157.  
  158.  
  159. Si GRAPH ;
  160.  
  161. dess evauz 'TITRE' 'Vitesse sur l axe';
  162. dess evarh 'TITRE' 'Densite melange sur l axe';
  163. dess evatf 'TITRE' 'Temperature sur l axe';
  164. dess evahe 'TITRE' 'Densite Helium sur l axe';
  165. dess evah2 'TITRE' 'Densite Hydrogene l axe';
  166.  
  167. ung= vect un 0.1 ux uy jaune;
  168. trace ung vtf;
  169. opti isov suli ;
  170. trace rho vtf 'TITRE' ' Rho' ;
  171. trace rhe vtf 'TITRE' ' Rhe' ;
  172. *trace rh2 vtf 'TITRE' ' Rh2' ;
  173. *trace rco vtf 'TITRE' ' Rco' ;
  174. *trace tf vtf 'TITRE' ' Tf ' ;
  175. Finsi ;
  176.  
  177. FIN ;
  178.  
  179.  
  180.  
  181.  
  182.  
  183.  
  184.  
  185.  
  186.  

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