Procedures | Cast3M

   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 :