Télécharger enc2dke.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : enc2dke.dgibi
  2. *
  3. * Enceinte 2D Axisymetrique
  4. * Relachement d'un melange d'helium et de vapeur
  5. * modèle K-Epsilon formulation Semi explicite (EFM1)
  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=4000 ;
  16. DT0 = 8.e-3 ;
  17. n1=5 ;
  18. n2=10;
  19. n3=20;
  20. GRAPH = VRAI ;
  21. 'SINON' ;
  22. nbit=100 ;
  23. DT0 = 8.e-3 ;
  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. RXT = 'TABLE' ;
  47. rxt.'VERSION'= 'V0' ;
  48.  
  49. RXT.'vtf'=mt ;
  50. RXT.'axe' = axe ;
  51. RXT.'breche'=breche ;
  52. RXT.'diru1'=0. 1. ;
  53. RXT.'epsi' = 1.E-5 ;
  54. RXT.'pi' = 0.5 0.5 ;
  55.  
  56. RXT.'DT0' = DT0 ;
  57. RXT.'DETMAT'= VRAI ;
  58. RXT.'RENU' = 'RIEN' ;
  59. RXT.'DISCR'= 'LINE';
  60. RXT.'KPRE' = 'MSOMMET';
  61. RXT.'ALGO' = 'EFM1';
  62. RXT.'GRAPH' = GRAPH ;
  63. RXT.'IMPR' = 0 ;
  64.  
  65. RXT.'MODTURB' = 'KEPSILON' ;
  66. RXT.'FPAROI' = VRAI ;
  67. RXT.'YP' = 0.05 ;
  68.  
  69. RXT.'VAPEUR' = VRAI ;
  70. RXT.'HE' = VRAI ;
  71.  
  72. RXT.'TF0' = 20.0 ;
  73. RXT.'PT0' = 1.0e5 ;
  74. RXT.'Yvap0' = 0.1 ;
  75. RXT.'Yhe0' = 0.0 ;
  76.  
  77. RXT.'scenario' = table ;
  78. RXT.'scenario'.'t' = prog 0.0 100.0 ;
  79. RXT.'scenario'.'qair' = prog 0.005 0.005 ;
  80. RXT.'scenario'.'qeau' = prog 0.0011 0.0011 ;
  81. RXT.'scenario'.'qhe' = prog 0.0011 0.0011 ;
  82. RXT.'scenario'.'tinj' = prog 20.0 20.0 ;
  83.  
  84. EXECRXT nbit RXT ;
  85.  
  86. un=RXT.'TIC'.'UN';
  87. mess ' Max U ' (maxi un) ' Min U ' (mini un) ;
  88. test = abs ((maxi un) - 1.02 );
  89. mess ' Test =' test;
  90. Si (test '>' 0.1 ) ; erreur 5 ; Finsi ;
  91.  
  92. tf =RXT.'TIC'.'TF' ;
  93. mess ' Max Tf' (maxi tf) ' Min Tf' (mini tf) ;
  94. test = abs ((maxi tf) - 20.1 );
  95. mess ' Test =' test;
  96. Si (test '>' 1. ) ; erreur 5 ; Finsi ;
  97.  
  98.  
  99.  
  100. $vtf=RXT.'GEO'.'$vtf' ;
  101. vtf=doma $vtf maillage;
  102. kn= RXT.'TIC'.'KN';
  103. en= RXT.'TIC'.'EN';
  104. rho=RXT.'TIC'.'RHO' ;
  105. rhe=RXT.'TIC'.'RHE';
  106. un =RXT.'TIC'.'UN';
  107.  
  108. ixe=inve axe ;
  109. evauz = evol chpo (exco un 'UY') ixe ;
  110. evarh = evol chpo rho ixe ;
  111. evatf = evol chpo tf ixe ;
  112. evahe = evol chpo rhe 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.  
  122. ung= vect un 0.1 ux uy jaune;
  123. trace ung vtf;
  124. opti isov suli ;
  125. trace rho vtf 'TITRE' ' Rho' ;
  126. trace rhe vtf 'TITRE' ' Rhe' ;
  127. trace tf vtf 'TITRE' ' Tf ' ;
  128. trace kn vtf 'TITRE' ' KN ' ;
  129. trace en vtf 'TITRE' ' EN ' ;
  130. Finsi ;
  131.  
  132. FIN ;
  133.  
  134.  
  135.  
  136.  
  137.  
  138.  
  139.  

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