Télécharger rayo-axi-4.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : rayo-axi-4.dgibi
  2.  
  3. *******************************************
  4. *******************************************
  5.  
  6. ***************************************************
  7. * *
  8. * Ce jeu de données permet la vérification du *
  9. * calcul des facteurs de forme dans le cas *
  10. * axisymétrique. On calcule le flux dû au *
  11. * rayonnement entre 2 sphères concentriques, la *
  12. * température de chacune des sphères étant *
  13. * homogène. Le résultat est comparé à la solution *
  14. * analytique, voir *
  15. * *
  16. * Jean Crabol, Transfert de Chaleur, tome 2, *
  17. * Masson, 1990, pp.169-175 *
  18. * *
  19. ***************************************************
  20.  
  21. **** Options ****
  22.  
  23. option dime 2 mode axis elem qua4 ;
  24.  
  25. graph = faux ;
  26.  
  27. ****************************************************
  28. **** DEFINITION DE LA GEOMETRIE DU PROBLEME ****
  29. ****************************************************
  30.  
  31. **** Dimensions de la sphere interieure ****
  32.  
  33. RIint = 97.2e-3;
  34. MRIint = -97.2e-3;
  35. RIext = 99.0e-3;
  36. MRIext = -99.0e-3;
  37.  
  38. **** Dimensions de la sphere exterieure ****
  39.  
  40. REint = 130.0e-3;
  41. MREint = -130.0e-3;
  42.  
  43. REext = 134.5e-3;
  44. MREext = -134.5e-3;
  45.  
  46. *** Nombres de divisions ****
  47.  
  48. n1 = 10 ; n2 = 2 ; n3 = 4 ;
  49.  
  50. *** Le centre ***
  51.  
  52. O = 0. 0. ;
  53.  
  54. **** Points sur l'enceinte interieure ****
  55.  
  56. P1 = 0. RIint;
  57. P2 = RIint 0.;
  58. P3 = 0. MRIint;
  59. P4 = 0. MRIext;
  60. P5 = RIext 0.;
  61. P6 = 0. RIext;
  62.  
  63. **** Points sur l'enceinte exterieure ****
  64.  
  65. Q1 = 0. REint;
  66. Q2 = REint 0.;
  67. Q3 = 0. MREint;
  68. Q4 = 0. MREext;
  69. Q5 = REext 0.;
  70. Q6 = 0. REext;
  71.  
  72. *** Les lignes (droites et arcs) ****
  73.  
  74. ARCP1 = cerc n1 P1 O P2 ;
  75. ARCP2 = cerc n1 P2 O P3 ;
  76. ARCII = ARCP1 et ARCP2 ;
  77.  
  78. ARCP3 = cerc n1 P4 O P5 ;
  79. ARCP4 = cerc n1 P5 O P6 ;
  80. ARCIE = ARCP3 et ARCP4 ;
  81.  
  82.  
  83. ARCQ1 = cerc n1 Q1 O Q2 ;
  84. ARCQ2 = cerc n1 Q2 O Q3 ;
  85. ARCEI = ARCQ1 et ARCQ2 ;
  86.  
  87. ARCQ3 = cerc n1 Q4 O Q5 ;
  88. ARCQ4 = cerc n1 Q5 O Q6 ;
  89. ARCEE = ARCQ3 et ARCQ4 ;
  90.  
  91. P6P1 = droi n2 P6 P1 ;
  92. P3P4 = droi n2 P3 P4 ;
  93.  
  94.  
  95. Q6Q1 = droi n2 Q6 Q1 ;
  96. Q3Q4 = droi n2 Q3 Q4 ;
  97.  
  98.  
  99. Q1P6 = droi n3 Q1 P6 ;
  100. P4Q3 = droi n3 P4 Q3 ;
  101.  
  102. *** Les surfaces ****
  103.  
  104. SPHERINT = P6P1 ARCII P3P4 ARCIE DALLE ;
  105.  
  106. SPHEREXT = Q6Q1 ARCEI Q3Q4 ARCEE DALLE ;
  107.  
  108. SPHERAIR = Q1P6 ARCIE P4Q3 ARCEI DALLE ;
  109.  
  110. iARCIE=inve ARCIE ;
  111. iARCEI=inve ARCEI ;
  112. cavite = iARCIE et iARCEI ;
  113.  
  114. TOUTacie = SPHERINT et SPHEREXT ;
  115. TOUT = TOUTacie ;
  116.  
  117. si(graph) ;
  118. titr 'Le maillage du modele' ;
  119. trac TOUT ;
  120. finsi ;
  121.  
  122. **** Rayonnement ****
  123.  
  124. MRI = MODELI iARCIE thermique RAYONNEMENT 'CAVITE' CONS 'CAV1';
  125. MRE = MODELI iARCEI thermique RAYONNEMENT 'CAVITE' CONS 'CAV1';
  126. mrt = mri et mre ;
  127.  
  128. vei = 0.9 ;
  129. vee = 0.3 ;
  130. ei = mate mri 'EMIS' vei ;
  131. ee = mate mre 'EMIS' vee ;
  132. chemis = ei et ee ;
  133.  
  134. tref = 0. ;
  135.  
  136. *** On testera sur le champ de température suivant ***
  137.  
  138. T_inter = 700.0 ;
  139. T_exter = 1000.0 ;
  140.  
  141. T_inter = T_inter + 273.0 ;
  142. T_exter = T_exter + 273.0 ;
  143.  
  144. chptint = manu chpo spherint 1 T T_inter nature discret ;
  145. chptext = manu chpo spherext 1 T T_exter nature discret ;
  146.  
  147. chptmpt = chptint et chptext ;
  148.  
  149. *** Les facteurs de forme ****************************
  150.  
  151. ff = ffor mrt chemis;
  152.  
  153. *** Matrice de rayonnement ***************************
  154.  
  155. mr = raye mrt ff chemis ;
  156.  
  157. *** Conductivité due au rayonnement ******************
  158.  
  159. cte_sb = 5.67e-8 ;
  160.  
  161. cr = rayn mrt mr (chan CHAM chptmpt mrt GRAVITE) cte_sb ;
  162.  
  163. *** Flux résultant ***********************************
  164.  
  165. flux12 = cr * chptmpt ;
  166.  
  167. flux1 = redu flux12 iARCIE ;
  168. flux2 = redu flux12 iARCEI ;
  169.  
  170. rsfl1 = maxi (resu flux1) ;
  171. rsfl2 = maxi (resu flux2) ;
  172.  
  173.  
  174. *** Solution analytique *******************************
  175.  
  176. solflux = cte_sb * 4 * pi * RIext * RIext * ((T_exter**4)-(T_inter**4));
  177. denom1 = 1.0 / vei ;
  178. denom2 = (RIext * RIext * ((1.0/vee) - 1.0)) / (REint*REint) ;
  179. denom = denom1 + denom2 ;
  180. solflux = solflux / denom ;
  181.  
  182. diffrel = 100 * ((abs rsfl2) - solflux) / solflux ;
  183.  
  184. opti echo 0 ;
  185. mess 'Flux arrivant sur la sphère intérieure' rsfl1 ;
  186. mess 'Flux partant de la sphère extérieure ' rsfl2 ;
  187. mess 'Solution analytique pour le flux = ' solflux ;
  188. mess 'Erreur obtenue est de ' diffrel '%';
  189. opti echo 1 ;
  190.  
  191. *** Test si c'est OK *********************************
  192.  
  193. si((abs diffrel) > 1.) ;
  194. erre 5 ;
  195. finsi;
  196.  
  197. *** Bye **********************************************
  198.  
  199. fin ;
  200.  
  201.  
  202.  
  203.  
  204.  
  205.  
  206.  
  207.  
  208.  

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