Télécharger Henc2d_therm1.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : Henc2d_therm1.dgibi
  2. *
  3. * Enceinte 2D Axisymetrique
  4. * Refroidissement d'un gaz chaud (AIR)
  5. * Interaction GAz/paroi/temperature exterieure
  6. * Teste thermique paroi et coefficient d'echange
  7.  
  8. 'OPTI' 'DIME' 2 'ELEM' 'CU20' ;
  9. 'OPTI' 'MODE' 'AXIS' ;
  10.  
  11. COMPLET= VRAI ;
  12. COMPLET= FAUX ;
  13. GRAPH = VRAI ;
  14. GRAPH = FAUX ;
  15.  
  16. 'SI' COMPLET ;
  17. nbit=200 ;
  18. DT0 = 2. ;
  19. n1=3 ;
  20. n2=5;
  21. n3=10;
  22. ds1= 0.02 ;
  23. ds2= 0.2 ;
  24. 'SINON' ;
  25. nbit=20 ;
  26. DT0 = 2. ;
  27. n1=3 ;
  28. n2=5;
  29. n3=10;
  30. ds1= 0.02 ;
  31. ds2= 0.2 ;
  32. 'FINSI' ;
  33.  
  34.  
  35. p1=5.e-4 0 ;
  36. p2=0.5 0. ;
  37. p3=1 0 ;
  38. p4=1 2;
  39. p41=0.5 2. ;
  40. p5=5.e-4 2 ;
  41.  
  42.  
  43. bas= p1 'DROI' DINI ds1 DFIN ds2 p2 d DINI ds2 DFIN ds1 p3 ;
  44. paroid=p3 'DROI' n3 p4 ;
  45. plaf=inve (bas plus (0. 2.) );
  46. axe = p5 'DROI' n3 p1 ;
  47. elim (bas et plaf et paroid et axe) 1.e-5 ;
  48. vtp = (paroid trans 2 (0.1 0.)) coul bleu;
  49. parext = (cote 3 vtp) coul rouge ;
  50.  
  51.  
  52. mt= 'DALLER' bas paroid plaf axe ;
  53. Si GRAPH;
  54. trace (mt et vtp);
  55. Finsi ;
  56. rxt = 'TABLE' ;
  57. rxt.'vtf'=mt ;
  58. rxt.'axe' = axe ;
  59.  
  60. rxt.'TF0' = 100. ;
  61. rxt.'PT0' = 1.0e5 ;
  62.  
  63. rxt.'THERMP'=VRAI ;
  64. rxt.'vtp'=vtp;
  65. rxt.'ROCP'=1.e4 ;
  66. rxt.'LAMBDA'=40. ;
  67. rxt.'Tp0' = 30. ;
  68. rxt.'ECHAN'=50. ;
  69.  
  70. rxt.'ECHEXT' = VRAI ;
  71. rxt.'HEXT'=50. ;
  72. rxt.'TPEXT'=20. ;
  73. rxt.'parext'=parext ;
  74.  
  75. rxt.'epsi'= 1.e-4 ;
  76. rxt.'pi' = 0.5 0.5 ;
  77. rxt.'DT0' = DT0 ;
  78. rxt.'RENU' = 'RIEN' ;
  79. rxt.'DISCR'= 'QUAF';
  80. rxt.'KPRE' = 'MSOMMET';
  81. rxt.'MODTURB' = 'NUTURB' ;
  82. rxt.'NUT' = 0.005 ;
  83. rxt.'GRAPH' = GRAPH ;
  84.  
  85. EXECRXT nbit rxt ;
  86.  
  87. un=rxt.'TIC'.'UN';
  88. tf =rxt.'TIC'.'TF' ;
  89.  
  90. $vtf=rxt.'GEO'.'$vtf' ;
  91. $vtp=rxt.'GEO'.'$vtp' ;
  92. vtf=doma $vtf maillage;
  93. vtp=doma $vtp maillage;
  94.  
  95. rho=rxt.'TIC'.'RHO' ;
  96. un =rxt.'TIC'.'UN';
  97. tp =rxt.'TIC'.'TP';
  98.  
  99. ixe=inve axe ;
  100. evauz = evol chpo (exco un 'UY') ixe ;
  101. evarh = evol chpo rho ixe ;
  102. evatf = evol chpo tf ixe ;
  103.  
  104.  
  105. Si (NON COMPLET) ;
  106. list rxt.TIC.'Tfm' ;
  107. list rxt.TIC.'PT' ;
  108. list rxt.TIC.'LMAXU';
  109.  
  110. ltfm=Prog
  111. 100.00 79.731 76.196 69.951 66.250
  112. 62.717 59.351 56.246 53.460 50.969
  113. 48.751 46.775 44.995 43.368 41.869
  114. 40.482 39.201 38.017 36.924 35.912
  115. 34.974;
  116.  
  117. lPT =Prog
  118. 1.00000E+05 98068. 92332. 91608. 91095.
  119. 90110. 89122. 88304. 87567. 86891.
  120. 86289. 85759. 85285. 84851. 84451.
  121. 84080. 83735. 83416. 83121. 82849.
  122. 82596.;
  123.  
  124. Lmaxu=Prog
  125. 0.0000 0.0000 0.29154 0.46508 0.52094
  126. 0.47754 0.41642 0.34661 0.28039 0.22271
  127. 0.18579 0.16626 0.15409 0.14588 0.13910
  128. 0.13245 0.12580 0.11938 0.11357 0.10855
  129. 0.10427;
  130.  
  131. tic=rxt.'TIC' ;
  132. ERtf=SOMM( abs (ltfm - tic.'Tfm') )/ 80. ;
  133. ERPT=SOMM( abs (lPT - tic.'PT' ) ) /1.e5 ;
  134. ERum=SOMM( abs (Lmaxu - tic.'LMAXU' ) ) ;
  135.  
  136. Mess ' ERtf=' ERtf ' ERPT=' ERPT ' ERum=' ERum ;
  137.  
  138. Si (ERtf '>' 2.e-4) ; erreur 5 ; Finsi ;
  139. Si (ERPT '>' 2.e-4) ; erreur 5 ; Finsi ;
  140. Si (ERum '>' 2.e-3) ; erreur 5 ; Finsi ;
  141.  
  142. evuz= extr evauz ordo ;
  143. list evuz ;
  144. lrr= prog
  145. 4.88090E-38 1.63432E-02 5.27145E-02 8.21404E-02 0.10017
  146. 0.11203 0.11314 0.10972 9.63041E-02 8.71612E-02
  147. 6.94869E-02 6.24003E-02 4.82692E-02 4.89579E-02 3.89615E-02
  148. 4.48514E-02 3.37560E-02 3.37986E-02 2.10628E-02 2.61349E-03
  149. -1.42428E-37;
  150.  
  151. ER=SOMM( (evuz - lrr)*(evuz - lrr) ) *0.05;
  152. mess ' Ecart ' er ;
  153. Si ( er > 2.e-4) ; erreur 5 ; finsi ;
  154.  
  155. Finsi ;
  156.  
  157.  
  158. Si GRAPH ;
  159.  
  160. dess evauz 'TITRE' 'Vitesse sur l axe';
  161. dess evarh 'TITRE' 'Densite melange sur l axe';
  162. dess evatf 'TITRE' 'Temperature sur l axe';
  163.  
  164. ung= vect un 0.1 ux uy jaune;
  165. trace ung vtf;
  166. trace tp vtp 'TITRE' ' Tp ' ;
  167. opti isov suli ;
  168. trace tf vtf 'TITRE' ' Tf ' ;
  169. trace rho vtf 'TITRE' ' Rho' ;
  170. Finsi ;
  171.  
  172. FIN ;
  173.  
  174.  
  175.  
  176.  
  177.  
  178.  
  179.  
  180.  
  181.  

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