Télécharger la plaquette Cast3M


fr en

Télécharger vapdif.procedur

Retour à la liste

Numérotation des lignes : 

  1. * VAPDIF    PROCEDUR  SERRE     19/10/14    21:15:02     10310          
  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 = 2.77 ;
  86. DVV = ( TR ** N ) * A / PM ;
  87. DV = DV + (XV * (INVE DVV)) ; 
  88.  
  89. DV = INVE DV ; 
  90. 'RESPROC' DV ; 
  91. 'FINPROC' ;
  92.  
  93.  
  94.  

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