Télécharger vapdif.procedur

Retour à la liste

Numérotation des lignes :

  1. * VAPDIF PROCEDUR MAGN 00/10/12 21:17:54 3973
  2. 'DEBPROC' VAPDIF PM*FLOTTANT TM*CHPOINT YV*CHPOINT
  3. YH2*CHPOINT YHE*CHPOINT YO2*CHPOINT
  4. YN2*CHPOINT YCO2*CHPOINT YCO*CHPOINT;
  5. *___________________________________________________________________
  6. * Procedure donnant le coefficient de diffusion de la vapeur dans le
  7. * melange
  8. *___________________________________________________________________
  9. * E/ PM : Pression totale (Pa)
  10. * TM : Temperature du melange (K)
  11. * YV : Fraction massique de vapeur
  12. * YH2 : Fraction massique de hydrogene
  13. * YHE : Fraction massique de helium
  14. * YO2 : Fraction massique de oxygene
  15. * YN2 : Fraction massique de azote
  16. * YCO2 : Fraction massique de CO2
  17. * YCO : Fraction massique de CO
  18. * /S DV : coefficient de diffusion en m2/s
  19. *___________________________________________________________________
  20. * Constante absolue des gaz parfaits (produit du nombre d'Avogadro par
  21. * la constante de Boltzmann) en J.mol-1.K-1
  22. RABS = 6.022*1.38066;
  23. * Masse molaire en kg.mol-1
  24. MVAP = 0.018 ;
  25. MH2 = 0.002 ;
  26. MHE = 0.004 ;
  27. MO2 = 0.032 ;
  28. MN2 = 0.028 ;
  29. MCO2 = 0.044 ;
  30. MCO = 0.028 ;
  31.  
  32. * Constantes relatives des différents gaz en J.kg-1.K-1
  33. RVAP = RABS/MVAP;
  34. RH2 = RABS/MH2;
  35. RHE = RABS/MHE;
  36. RO2 = RABS/MO2;
  37. RN2 = RABS/MN2;
  38. RCO2 = RABS/MCO2;
  39. RCO = RABS/MCO;
  40.  
  41. * Constante relative du mélange
  42. Rmel = (YV*RVAP) + (YH2*RH2) + (YHE*RHE) + (YO2*RO2) + (YN2*RN2)
  43. + (YCO2*RCO2) + (YCO*RCO) ;
  44.  
  45. * Masse molaire du melange
  46. Mmel = RABS*(Rmel**(-1.));
  47.  
  48. * Fractions molaires
  49. XV = YV * Mmel / MVAP ;
  50. XH2 = YH2 * Mmel / MH2 ;
  51. XHE = YHE * Mmel / MHE ;
  52. XO2 = YO2 * Mmel / MO2 ;
  53. XN2 = YN2 * Mmel / MN2 ;
  54. XCO2 = YCO2 * Mmel / MCO2 ;
  55. XCO = YCO * Mmel / MCO ;
  56.  
  57. * Calcul du coefficient de diffusion DV
  58. DV = (YV * 0.0) ;
  59. TR = TM / 273.0 ;
  60. N = 1.75 ;
  61. A = 7.8 ;
  62. DVH2 = ( TR ** N ) * A / PM ;
  63. DV = DV + (XH2 * (INVE DVH2)) ;
  64. N = 1.75 ;
  65. A = 7.3 ;
  66. DVHE = ( TR ** N ) * A / PM ;
  67. DV = DV + (XHE * (INVE DVHE)) ;
  68. N = 1.71 ;
  69. A = 2.4 ;
  70. DVO2 = ( TR ** N ) * A / PM ;
  71. DV = DV + (XO2 * (INVE DVO2)) ;
  72. N = 1.75 ;
  73. A = 2.27 ;
  74. DVN2 = ( TR ** N ) * A / PM ;
  75. DV = DV + (XN2 * (INVE DVN2)) ;
  76. N = 1.75 ;
  77. A = 1.55 ;
  78. DVCO2 = ( TR ** N ) * A / PM ;
  79. DV = DV + (XCO2 * (INVE DVCO2)) ;
  80. N = 1.75 ;
  81. A = 2.2 ;
  82. DVCO = ( TR ** N ) * A / PM ;
  83. DV = DV + (XCO * (INVE DVCO)) ;
  84. N = 0.0 ;
  85. A = 27.7 ;
  86. DVV = ( TR ** N ) * A / PM ;
  87. DV = DV + (XV * (INVE DVV)) ;
  88.  
  89. DV = INVE DV ;
  90. 'RESPROC' DV ;
  91. 'FINPROC' ;
  92.  
  93.  

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