Télécharger trkg.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : trkg.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. * Enceinte 2D Axisymetrique Cas de Torrance et Rockett
  6. *
  7. 'OPTI' 'DIME' 2 'ELEM' 'CU20' ;
  8. 'OPTI' 'MODE' 'AXIS' ;
  9.  
  10. opti trac 'PSC' ;
  11. COMPLET= VRAI ;
  12. COMPLET= FAUX ;
  13. GRAPH = VRAI ;
  14. GRAPH = FAUX ;
  15.  
  16. echl=2.e-1;
  17. *echl=5.e-1;
  18. *echl=1.e-1;
  19. 'SI' COMPLET ;
  20. nbit=100 ;
  21. nbit=20 ;
  22. DT0 = 1. ;
  23. ds1= 0.02*echl;
  24. ds2= 0.2*echl ;
  25. n1=5 ;
  26. n2=10;
  27. n3=20;
  28. 'SINON' ;
  29. nbit=100 ;
  30. DT0 = 40.;
  31. ds1= 0.02*echl;
  32. ds2= 0.2*echl ;
  33. n1=5 ;
  34. n2=10;
  35. n3=20;
  36. 'FINSI' ;
  37.  
  38. r0=0.1*echl ;
  39. r1=1.*echl ;
  40. rm=0.5*echl ;
  41. hm=0.5*echl ;
  42. h1=1.*echl ;
  43.  
  44. eps=1.e-2*echl ;
  45. eps=1.e-4*echl ;
  46. epsi=1.e-4*echl ;
  47.  
  48. p1=eps 0 ;
  49. p11=r0 0 ;
  50. p2=rm 0. ;
  51. p3=r1 0 ;
  52. p4=r1 hm;
  53. p5=r1 h1 ;
  54.  
  55.  
  56. breche = p1 'DROI' n1 p11;
  57. bas= p11 'DROI' dini ds1 dfin ds2 p2 d dini ds2 dfin ds1 p3 ;
  58. paroid=p3 'DROI' dini ds1 dfin ds2 p4 'DROI' dini ds2 dfin ds1 p5 ;
  59.  
  60. plaf=inve ( (breche et bas ) plus (0 h1));
  61. axe =inve ( paroid moins ((r1-eps) 0)) ;
  62. elim (breche et bas et paroid et plaf et axe) epsi ;
  63. mt= 'DALLER' (breche 'ET' bas) paroid plaf axe ;
  64. Si GRAPH ; trace mt ; Finsi ;
  65.  
  66. rxt = 'TABLE' ;
  67. rxt.'vtf'=mt ;
  68. rxt.'axe' = axe ;
  69. rxt.'epsi'= epsi ;
  70. rxt.'pi' = 0.5 0.5 ;
  71. rxt.'DT0' = DT0 ;
  72. rxt.'DISCR'= 'QUAF';
  73. rxt.'KPRE' = 'MSOMMET';
  74. rxt.'TYPINV'=1 ;
  75. rxt.'GRAPH' = GRAPH ;
  76. rxt.'MODTURB'='NUTURB' ;
  77. rxt.'NUT' =1.e-7 ;
  78. rxt.'CORTEMP'=FAUX ;
  79.  
  80. rxt.'TF0' = 20.0 ;
  81. rxt.'PT0' = 1.0e5 ;
  82.  
  83. rxt.'TIMP1' = table ;
  84. rxt.'TIMP1'.'MAILLAGE' = breche ;
  85. rxt.'TIMP1'.'t' = prog 0.0 900000.0 ;
  86. rxt.'TIMP1'.'TIMP' = prog 21. 21. ;
  87. rxt.'TIMP1'.'ECHAN' = 1.e2 ;
  88.  
  89. rxt.'TIMP2' = table ;
  90. rxt.'TIMP2'.'MAILLAGE' = (bas et paroid et plaf) ;
  91. rxt.'TIMP2'.'t' = prog 0.0 900000.0 ;
  92. rxt.'TIMP2'.'TIMP' = prog 20.0 20.0 ;
  93. rxt.'TIMP2'.'ECHAN' = 1.e2 ;
  94. execrxt nbit rxt ;
  95.  
  96.  
  97.  
  98. $vtf=rxt.'GEO'.'$vtf' ;
  99. vtf=doma $vtf maillage;
  100. rho=rxt.'TIC'.'RHO' ;
  101. tf =rxt.'TIC'.'TF' ;
  102. un =rxt.'TIC'.'UN';
  103.  
  104. $axe=rxt.'GEO'.'$axe';
  105. axe=doma $axe maillage ;
  106. ixe=inve axe ;
  107. evauz = evol chpo (exco un 'UY') ixe ;
  108. evarh = evol chpo rho ixe ;
  109. evatf = evol chpo tf ixe ;
  110.  
  111.  
  112. luz= (extr evauz 'ORDO' );
  113. lrh= (extr evarh 'ORDO' );
  114. ltf= (extr evatf 'ORDO' );
  115. list luz ; list lrh ; list ltf ;
  116.  
  117. luzr= prog
  118. 3.07237E-42 -2.34123E-03 2.56887E-03 6.54417E-03 1.62324E-02
  119. 2.40888E-02 2.50937E-02 2.65036E-02 2.59309E-02 2.53052E-02
  120. 2.60378E-02 2.60663E-02 2.63601E-02 2.57449E-02 2.55996E-02
  121. 2.49561E-02 2.43348E-02 2.32284E-02 2.10423E-02 1.78512E-02
  122. 1.39740E-02 8.69525E-03 3.46774E-03 7.10408E-04 -5.08355E-43;
  123.  
  124. lrhr=prog
  125. 1.1842 1.1845 1.1848 1.1852 1.1855
  126. 1.1859 1.1862 1.1865 1.1868 1.1870
  127. 1.1872 1.1874 1.1875 1.1876 1.1877
  128. 1.1877 1.1878 1.1878 1.1878 1.1878
  129. 1.1879 1.1879 1.1879 1.1880 1.1882;
  130.  
  131. ltfr=prog
  132. 20.994 20.905 20.829 20.740 20.651
  133. 20.565 20.482 20.404 20.337 20.275
  134. 20.223 20.177 20.142 20.115 20.095
  135. 20.082 20.072 20.065 20.058 20.054
  136. 20.050 20.045 20.036 20.020 20.001;
  137.  
  138. ERUZ=(SOMM (abs (luz - luzr)))*0.04 ;
  139. MESS 'ERUZ=' ERUZ ;
  140. Si (ERUZ '>' 2.e-3) ; erreur 5 ; Finsi ;
  141.  
  142. ERRH=(SOMM (abs (lrh - lrhr)))*0.04 ;
  143. MESS 'ERRH=' ERRH ;
  144. Si (ERRH '>' 2.e-3) ; erreur 5 ; Finsi ;
  145.  
  146. ERTF=(SOMM (abs (ltf - ltfr)))*0.04/20. ;
  147. MESS 'ERTF=' ERTF ;
  148. Si (ERTF '>' 2.e-3) ; erreur 5 ; Finsi ;
  149.  
  150.  
  151.  
  152.  
  153.  
  154. Si GRAPH ;
  155.  
  156. dess evauz 'TITRE' 'Vitesse sur l axe';
  157. dess evarh 'TITRE' 'Densite melange sur l axe';
  158. dess evatf 'TITRE' 'Temperature sur l axe';
  159.  
  160. ung= vect un 1. ux uy jaune;
  161. trace ung vtf;
  162. opti isov suli ;
  163. trace rho vtf 'TITRE' ' Rho' ;
  164. trace tf vtf (prog 20. pas 0.001 20.02);
  165. trace tf vtf (prog 20. 20.0001 20.001 20.005 20.01 20.1) ;
  166. *trace rh2 vtf 'TITRE' ' Rh2' ;
  167. *trace rco2 vtf 'TITRE' ' Rco2' ;
  168. *trace tf vtf 'TITRE' ' Tf ' ;
  169. Finsi ;
  170.  
  171. FIN ;
  172.  
  173.  
  174.  
  175.  
  176.  
  177.  
  178.  
  179.  
  180.  
  181.  
  182.  
  183.  
  184.  

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