Télécharger rayo_abs-axi-1.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : rayo_abs-axi-1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. ***********************************************************************
  6. * *
  7. * Rayonnement thermique en milieu absorbant dans une cavité sphérique*
  8. * (pas de couplage avec d'autres modes de transfet d'énergie) *
  9. * *
  10. * Comparaison à un calcul analytique *
  11. * Ref: Siegel&Howell Ed.3 p609-615 *
  12. * *
  13. * On evalue la puissance perdue par un gaz absorbant de température *
  14. * uniforme 2273K et de coefficient d'absorption 100./m contenu dans *
  15. * cavité de rayon 0.01m à la température uniforme de 1273K. *
  16. * *
  17. * Calcul 2D axisymétrique *
  18. * Remarque: le calcul des facteurs de forme n'utilise pas l'option *
  19. * convexe 'CONV' *
  20. ***********************************************************************
  21. *
  22. *
  23. option echo 1 dime 2 elem qua4 mode axis ;
  24. graph = faux ;
  25. *
  26. *-------------------------------
  27. * Maillage
  28. *-------------------------------
  29. *
  30. R1 = 1.E-2 ;
  31. O = 0. 0. ;
  32. p1 = 0. ( -1. * R1 ) ;
  33. p2 = R1 0. ;
  34. p3 = 0. R1 ;
  35.  
  36. d = 1.E-3 ;
  37.  
  38. p1p2 = cerc p1 O p2 DINI d DFIN d;
  39. p2p3 = cerc p2 O p3 DINI d DFIN d;
  40. sphe_ext = p1p2 et p2p3 ;
  41.  
  42. si graph ;
  43. trac sphe_ext ;
  44. finsi ;
  45.  
  46. cavite = sphe_ext ;
  47. tout = cavite;
  48.  
  49. *-------------------------------
  50. * Propriétés physiques
  51. *-------------------------------
  52.  
  53. e_wall = 0.5 ;
  54. abso0 = -100. ;
  55. T_wall = 1273. ;
  56. T_gas = 2273. ;
  57.  
  58. *-------------------------------
  59. * Modèle de rayonnement
  60. *-------------------------------
  61.  
  62. mrt = mode sphe_ext thermique rayonnement 'CAVITE' ;
  63. e = mate mrt 'EMIS' e_wall 'CABS' abso0 'TABS' T_gas;
  64.  
  65. * opti 'IMPI' 1 ;
  66.  
  67. *-------------------------------
  68. * Facteurs de forme et matrice de rayonnement
  69. *-------------------------------
  70.  
  71. fft = ffor mrt e;
  72. chamr = raye mrt fft e ;
  73.  
  74. *-------------------------------
  75. * gaz absorbant : calcul du terme R*Tg4
  76. *-------------------------------
  77.  
  78. tg = manu chpo tout 1 'T' T_gas natu 'DIFFUS';
  79. tg_cavi = redu tg cavite ;
  80. tge_cavi = chan 'CHAM' tg_cavi mrt 'GRAVITE' ;
  81. crg= rayn mrt chamr tge_cavi ;
  82. fg = crg*tg_cavi ;
  83.  
  84. *-------------------------------
  85. * paroi : calcul du terme R*Tw4
  86. *-------------------------------
  87.  
  88. tp = manu chpo tout 1 'T' T_wall natu 'DIFFUS';
  89. t_cavi = redu tp cavite ;
  90. te_cavi = chan 'CHAM' t_cavi mrt 'GRAVITE' ;
  91. cr(<3pan style="color: #000000; font-weight: bold;">= rayn mrt chamr te_cavi ;
  92. fw = cr *t_cavi ;
  93.  
  94. *-------------------------------
  95. * puissance perdue par les frontieres:
  96. *-------------------------------
  97.  
  98. puis_n = fw - fg ;
  99. t1 = manu chpo tout 1 'T' 1.0 natu 'DIFFUS' ;
  100. lm1 = MOTS 'Q' ;
  101. lm2 = MOTS 'T' ;
  102. puis_t = xty puis_n t1 lm1 lm2 ;
  103.  
  104. *-------------------------------
  105. * puissance théorique perdue par les frontieres:
  106. *-------------------------------
  107.  
  108. abso0 = -1 * abso0 ;
  109. b0 = 5.670E-8 ;
  110. EMG1 = b0 * ( T_gas ** 4. ) ;
  111. EMS1 = b0 * ( T_wall ** 4. ) ;
  112. AIRE1 = 4. * 3.1416 * R1 * R1 ;
  113. tau1 = 2. / ( ( 2. * abso0 * R1 ) ** 2. ) ;
  114. tau0 = ( 2. * abso0 * R1 ) + 1. ;
  115. tau0 = tau0 * ( 'EXP' ( -2. * abso0 * R1 ) ) ;
  116. tau1 = ( 1. - tau0 ) * tau1 ;
  117. abso1 = 1. - tau1 ;
  118. denom0 = ( 1. / e_wall ) + ( 1. / abso1 ) - 1. ;
  119. res0 = -1. * AIRE1 * ( EMG1 - EMS1 ) / denom0 ;
  120.  
  121.  
  122. *-------------------------------
  123. * Erreur methode 1
  124. *-------------------------------
  125.  
  126. si ( ( 'ABS' res0 ) '>' 1.E-5 ) ;
  127. err1 = ( 'ABS' ( res0 - puis_t ) ) / ( 'ABS' res0 ) ;
  128. err1 = err1 * 100. ;
  129. sinon ;
  130. err1 = 0. ;
  131. finsi ;
  132. 'MESS' ' Puissance absorbée calculée methode 1 = ' puis_t ;
  133. 'MESS' ' Puissance absorbee théorique = ' res0 ;
  134. 'MESS' ' Erreur methode 1 en pourcentage = ' err1 ' % ' ;
  135.  
  136.  
  137. *-------------------------------
  138. * calcul du flux rayonné au moyen de l'évaluation de la
  139. * température de rayonnement (option 2 de RAYE)
  140. *-------------------------------
  141.  
  142. trad = raye mrt fft e te_cavi 1e-7 ;
  143. * mess 'trad ' (mini trad) (maxi trad);
  144.  
  145. * calcul du coefficient d'echange
  146.  
  147. hrad = HRCAV mrt e te_cavi trad ;
  148. * mess 'hrad ' (mini hrad) (maxi hrad);
  149.  
  150. trad_n1= chan 'CHPO' mrt trad ;
  151. trad_n = nomc trad_n1 'T' 'NATU' 'DIFFUS';
  152.  
  153. * pour la condition de convection
  154. cr = cond mrt hrad;
  155. f = conv mrt hrad trad_n;
  156.  
  157. * flux rayonne
  158.  
  159. fray = (cr * tp)- f;
  160. fray_tot = maxi (resu fray);
  161.  
  162. mess ' flux par methode 2 ' fray_tot ;
  163.  
  164. *-------------------------------
  165. * Erreur methode 2
  166. *-------------------------------
  167.  
  168. si ( ( 'ABS' res0 ) '>' 1.E-5 ) ;
  169. err2 = ( 'ABS' ( res0 - fray_tot) ) / ( 'ABS' res0 ) ;
  170. err2 = err2 * 100. ;
  171. sinon ;
  172. err2 = 0. ;
  173. finsi ;
  174. 'MESS' ' Puissance absorbée calculée methode 2 = ' fray_tot ;
  175. 'MESS' ' Puissance absorbee théorique = ' res0 ;
  176. 'MESS' ' Erreur methode 2 en pourcentage = ' err2 ' % ' ;
  177.  
  178.  
  179. si (( err1 '<' .2 ) et ( err2 '<' .2 ));
  180. 'ERRE' 0 ;
  181. sinon ;
  182. 'ERRE' 5 ;
  183. finsi ;
  184.  
  185. fin ;
  186.  
  187.  
  188.  
  189.  
  190.  
  191.  

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