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

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