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.  
  48. RXT.'vtf'=mt ;
  49. RXT.'axe' = axe ;
  50. RXT.'breche'=breche ;
  51. RXT.'diru1'=0. 1. ;
  52. RXT.'epsi' = 1.E-5 ;
  53. RXT.'pi' = 0.5 0.5 ;
  54.  
  55. RXT.'DT0' = DT0 ;
  56. RXT.'DETMAT'= VRAI ;
  57. RXT.'RENU' = 'RIEN' ;
  58. RXT.'DISCR'= 'LINE';
  59. RXT.'KPRE' = 'MSOMMET';
  60. RXT.'ALGO' = 'EFM1';
  61. RXT.'GRAPH' = GRAPH ;
  62. RXT.'IMPR' = 0 ;
  63.  
  64. RXT.'MODTURB' = 'KEPSILON' ;
  65. RXT.'FPAROI' = VRAI ;
  66. RXT.'YP' = 0.05 ;
  67.  
  68. RXT.'VAPEUR' = VRAI ;
  69. RXT.'HE' = VRAI ;
  70.  
  71. RXT.'TF0' = 20.0 ;
  72. RXT.'PT0' = 1.0e5 ;
  73. RXT.'Yvap0' = 0.1 ;
  74. RXT.'Yhe0' = 0.0 ;
  75.  
  76. RXT.'scenario' = table ;
  77. RXT.'scenario'.'t' = prog 0.0 100.0 ;
  78. RXT.'scenario'.'qair' = prog 0.005 0.005 ;
  79. RXT.'scenario'.'qeau' = prog 0.0011 0.0011 ;
  80. RXT.'scenario'.'qhe' = prog 0.0011 0.0011 ;
  81. RXT.'scenario'.'tinj' = prog 20.0 20.0 ;
  82.  
  83. EXECRXT nbit RXT ;
  84.  
  85. un=RXT.'TIC'.'UN';
  86. mess ' Max U ' (maxi un) ' Min U ' (mini un) ;
  87. test = abs ((maxi un) - 1.02 );
  88. mess ' Test =' test;
  89. Si (test '>' 0.1 ) ; erreur 5 ; Finsi ;
  90.  
  91. tf =RXT.'TIC'.'TF' ;
  92. mess ' Max Tf' (maxi tf) ' Min Tf' (mini tf) ;
  93. test = abs ((maxi tf) - 20.1 );
  94. mess ' Test =' test;
  95. Si (test '>' 1. ) ; erreur 5 ; Finsi ;
  96.  
  97.  
  98.  
  99. $vtf=RXT.'GEO'.'$vtf' ;
  100. vtf=doma $vtf maillage;
  101. kn= RXT.'TIC'.'KN';
  102. en= RXT.'TIC'.'EN';
  103. rho=RXT.'TIC'.'RHO' ;
  104. rhe=RXT.'TIC'.'RHE';
  105. un =RXT.'TIC'.'UN';
  106.  
  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.  
  113.  
  114. Si GRAPH ;
  115.  
  116. dess evauz 'TITRE' 'Vitesse sur l axe';
  117. dess evarh 'TITRE' 'Densite melange sur l axe';
  118. dess evatf 'TITRE' 'Temperature sur l axe';
  119. dess evahe 'TITRE' 'Densite Helium sur l axe';
  120.  
  121. ung= vect un 0.1 ux uy jaune;
  122. trace ung vtf;
  123. opti isov suli ;
  124. trace rho vtf 'TITRE' ' Rho' ;
  125. trace rhe vtf 'TITRE' ' Rhe' ;
  126. trace tf vtf 'TITRE' ' Tf ' ;
  127. trace kn vtf 'TITRE' ' KN ' ;
  128. trace en vtf 'TITRE' ' EN ' ;
  129. Finsi ;
  130.  
  131. FIN ;
  132.  
  133.  
  134.  
  135.  
  136.  

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