Télécharger enc2D-therco.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : enc2D-therco.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. * Enceinte 2D Axisymetrique
  6. * Refroidissement d'un gaz chaud (AIR)
  7. * Interaction GAz/paroi/temperature exterieure
  8. * Teste thermique paroi et coefficient d'echange
  9.  
  10. 'OPTI' 'DIME' 2 'ELEM' 'CU20' ;
  11. 'OPTI' 'MODE' 'AXIS' ;
  12.  
  13. COMPLET= VRAI ;
  14. COMPLET= FAUX ;
  15. GRAPH = VRAI ;
  16. GRAPH = FAUX ;
  17. opti trac 'PSC';
  18.  
  19. 'SI' COMPLET ;
  20. nbit=200 ;
  21. DT0 = 2. ;
  22. n1=3 ;
  23. n2=5;
  24. n3=10;
  25. ds1= 0.02 ;
  26. ds2= 0.2 ;
  27. 'SINON' ;
  28. nbit=20 ;
  29. DT0 = 2. ;
  30. n1=3 ;
  31. n2=5;
  32. n3=10;
  33. ds1= 0.02 ;
  34. ds2= 0.2 ;
  35. 'FINSI' ;
  36.  
  37.  
  38. p1=5.e-4 0 ;
  39. p2=0.5 0. ;
  40. p3=1 0 ;
  41. p4=1 2;
  42. p41=0.5 2. ;
  43. p5=5.e-4 2 ;
  44.  
  45.  
  46. bas= p1 'DROI' DINI ds1 DFIN ds2 p2 d DINI ds2 DFIN ds1 p3 ;
  47. paroid=p3 'DROI' n3 p4 ;
  48. plaf=inve (bas plus (0. 2.) );
  49. axe = p5 'DROI' n3 p1 ;
  50. elim (bas et plaf et paroid et axe) 1.e-5 ;
  51. vtp = (paroid trans 2 (0.1 0.)) coul bleu;
  52. parext = (cote 3 vtp) coul rouge ;
  53.  
  54.  
  55. mt= 'DALLER' bas paroid plaf axe ;
  56. Si GRAPH;
  57. trace (mt et vtp);
  58. Finsi ;
  59. rxt = 'TABLE' ;
  60. rxt.'VERSION'= 'V0' ;
  61. rxt.'vtf'=mt ;
  62. rxt.'axe' = axe ;
  63.  
  64. rxt.'TF0' = 100. ;
  65. rxt.'PT0' = 1.0e5 ;
  66.  
  67. rxt.'THERMP'=VRAI ;
  68. rxt.'THERCO'=VRAI ;
  69. rxt.'vtp'=vtp;
  70. rxt.'ROCP'=1.e4 ;
  71. rxt.'LAMBDA'=40. ;
  72. rxt.'Tp0' = 30. ;
  73. rxt.'ECHAN'=50. ;
  74.  
  75. rxt.'ECHEXT' = VRAI ;
  76. rxt.'HEXT'=50. ;
  77. rxt.'TPEXT'=20. ;
  78. rxt.'parext'=parext ;
  79.  
  80. rxt.'epsi'= 1.e-4 ;
  81. rxt.'pi' = 0.5 0.5 ;
  82. rxt.'DT0' = DT0 ;
  83. rxt.'RENU' = 'RIEN' ;
  84. rxt.'DISCR'= 'QUAF';
  85. rxt.'KPRE' = 'MSOMMET';
  86. rxt.'MODTURB' = 'NUTURB' ;
  87. rxt.'NUT' = 0.005 ;
  88. rxt.'GRAPH' = GRAPH ;
  89.  
  90. EXECRXT 0 rxt ;
  91. $vtp=rxt.'GEO'.'$vtp' ;
  92. $vtf=rxt.'GEO'.'$vtf' ;
  93. *rxt.'LAMBDA'=kcht $vtf scal sommet 40. ;
  94.  
  95. EXECRXT nbit rxt ;
  96.  
  97. un=rxt.'TIC'.'UN';
  98. tf =rxt.'TIC'.'TF' ;
  99.  
  100. $vtf=rxt.'GEO'.'$vtf' ;
  101. $vtp=rxt.'GEO'.'$vtp' ;
  102. vtf=doma $vtf maillage;
  103. vtp=doma $vtp maillage;
  104.  
  105. rho=rxt.'TIC'.'RHO' ;
  106. un =rxt.'TIC'.'UN';
  107. tp =rxt.'TIC'.'TP';
  108.  
  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. Si (NON COMPLET) ;
  116. list rxt.TIC.'Tfm' ;
  117. list rxt.TIC.'PT' ;
  118. list rxt.TIC.'LMAXU';
  119.  
  120.  
  121. ltfm=Prog
  122. 100.00 79.731 77.454 71.786 67.833
  123. 64.355 61.029 57.925 55.136 52.636
  124. 50.398 48.391 46.573 44.903 43.359
  125. 41.927 40.602 39.376 38.240 37.187
  126. 36.207;
  127.  
  128. lPT =Prog
  129. 1.00000E+05 98068. 92332. 91994. 91665.
  130. 90512. 89494. 88730. 88003. 87312.
  131. 86707. 86175. 85692. 85248. 84837.
  132. 84455. 84099. 83770. 83464. 83182.
  133. 82918.;
  134.  
  135. Lmaxu=Prog
  136. 0.0000 0.0000 0.29154 0.46506 0.52465
  137. 0.48430 0.42198 0.35131 0.28534 0.22824
  138. 0.19253 0.17409 0.16273 0.15506 0.14841
  139. 0.14168 0.13494 0.12844 0.12257 0.11747
  140. 0.11309;
  141.  
  142. tic=rxt.'TIC' ;
  143. ERtf=SOMM( abs (ltfm - tic.'Tfm') )/ 80. ;
  144. ERPT=SOMM( abs (lPT - tic.'PT' ) ) /1.e5 ;
  145. ERum=SOMM( abs (Lmaxu - tic.'LMAXU' ) ) ;
  146.  
  147. Mess ' ERtf=' ERtf ' ERPT=' ERPT ' ERum=' ERum ;
  148.  
  149. *Si (ERtf '>' 2.e-4) ; erreur 5 ; Finsi ;
  150. *Si (ERPT '>' 2.e-4) ; erreur 5 ; Finsi ;
  151. *Si (ERum '>' 2.e-3) ; erreur 5 ; Finsi ;
  152.  
  153. evuz= extr evauz ordo ;
  154. list evuz ;
  155. lrr= prog
  156. 4.88090E-38 1.63432E-02 5.27145E-02 8.21404E-02 0.10017
  157. 0.11203 0.11314 0.10972 9.63041E-02 8.71612E-02
  158. 6.94869E-02 6.24003E-02 4.82692E-02 4.89579E-02 3.89615E-02
  159. 4.48514E-02 3.37560E-02 3.37986E-02 2.10628E-02 2.61349E-03
  160. -1.42428E-37;
  161.  
  162. ER=SOMM( (evuz - lrr)*(evuz - lrr) ) *0.05;
  163. mess ' Ecart ' er ;
  164. *Si ( er > 2.e-4) ; erreur 5 ; finsi ;
  165.  
  166. Finsi ;
  167.  
  168.  
  169. Si GRAPH ;
  170.  
  171. dess evauz 'TITRE' 'Vitesse sur l axe';
  172. dess evarh 'TITRE' 'Densite melange sur l axe';
  173. dess evatf 'TITRE' 'Temperature sur l axe';
  174.  
  175. ung= vect un 0.1 ux uy jaune;
  176. trace ung vtf;
  177. trace tp vtp 'TITRE' ' Tp ' ;
  178. opti isov suli ;
  179. trace tf vtf 'TITRE' ' Tf ' ;
  180. trace rho vtf 'TITRE' ' Rho' ;
  181. Finsi ;
  182.  
  183. FIN ;
  184.  
  185.  
  186.  
  187.  
  188.  
  189.  
  190.  
  191.  
  192.  
  193.  
  194.  
  195.  
  196.  
  197.  
  198.  

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