Procedures | Cast3M

   1 : * CALMU     PROCEDUR  MAGN      06/02/17    21:15:01     5314           
   2 : DEBPROC CALMU ;
   3 : 'ARGU' T/'CHPOINT'
   4 :        TF/'FLOTTANT'
   5 :        TL/'LISTREEL' ;
   6 : *-------------------------------------------------*
   7 : * calcul de la viscosite dynamique                *
   8 : *           pour les incondensables en fonction   *
   9 : * de la temperature  (J/kg/K)                     *
  10 : *                                                 *
  11 : * Entree T en K                                   *
  12 : * Sortie Mu en Kg/m/s                             *
  13 : *                                                 *
  14 : * PS : le resultat est du meme type que l'entree  *
  15 : *-------------------------------------------------*
  16 : A0 = 0.1541D-6 ;
  17 : A1 = -17.856 ;
  18 : A2 = 1.703   ;
  19 : B0 = 0.5865D-6 ;
  20 : B1 = 27.183 ;
  21 : B2 = 1.630 ;
  22 : C0 = 0.7406D-6 ;
  23 : C1 = 48.16 ;
  24 : C2 = 1.610  ;
  25 : D0 = 1.791D-6 ;
  26 : D1 = 140.812  ;
  27 : D2 = 1.469   ;
  28 : E0 = 0.136D-6 ;
  29 : E1 = -59.30   ;
  30 : E2 = 1.810     ;
  31 : F0 = 0.573D-6 ;
  32 : F1 = 46.45    ;
  33 : F2 = 1.620     ;
  34 : G0 = 2.142D-6 ;
  35 : G1 = 890.0 ;
  36 : G2 = 1.500 ;
  37 : 
  38 : 'SI' ( 'EXIST' T ) ;
  39 :     Muh2  =(((T ** A2 ) * A0 ) / ( T + A1 ))           ;
  40 :     Muhe  =(((T ** B2 ) * B0 ) / ( T + B1 ))           ;
  41 :     Muo2  =(((T ** C2 ) * C0 ) / ( T + C1 ))           ;
  42 :     Mun2  =(((T ** D2 ) * D0 ) / ( T + D1 ))           ;
  43 :     Muco2 =(((T ** E2 ) * E0 ) / ( T + E1 ))           ;
  44 :     Muco  =(((T ** F2 ) * F0 ) / ( T + F1 ))           ;
  45 :     Muvap =(((T ** G2 ) * G0 ) / ( T + G1 ))           ;
  46 :     Muair = (0.8*Mun2) + (0.2*Muo2) ;
  47 :     'RESPROC' Muh2 Muhe Muo2 Mun2 Muco2 Muco Muvap Muair ;
  48 :     'QUITTER' CALMU ;
  49 : 'FINSI' ;
  50 : 'SI' ( 'EXIST' TF) ;
  51 :     Muh2  =(((TF** A2 ) * A0 ) / ( TF+ A1 ))           ;
  52 :     Muhe  =(((TF** B2 ) * B0 ) / ( TF+ B1 ))           ;
  53 :     Muo2  =(((TF** C2 ) * C0 ) / ( TF+ C1 ))           ;
  54 :     Mun2  =(((TF** D2 ) * D0 ) / ( TF+ D1 ))           ;
  55 :     Muco2 =(((TF** E2 ) * E0 ) / ( TF+ E1 ))           ;
  56 :     Muco  =(((TF** F2 ) * F0 ) / ( TF+ F1 ))           ;
  57 :     Muvap =(((TF** G2 ) * G0 ) / ( TF+ G1 ))           ;
  58 :     Muair = (0.8*Mun2) + (0.2*Muo2) ;
  59 :     'RESPROC' Muh2 Muhe Muo2 Mun2 Muco2 Muco Muvap Muair ;
  60 :     'QUITTER' CALMU ;
  61 : 'FINSI' ;
  62 : 'SI' ( 'EXIST' TL) ;
  63 :     T00 = TL / TL ;
  64 :     A10 = A1 * T00 ;
  65 :     B10 = B1 * T00 ;
  66 :     C10 = C1 * T00 ;
  67 :     D10 = D1 * T00 ;
  68 :     E10 = E1 * T00 ;
  69 :     F10 = F1 * T00 ;
  70 :     G10 = G1 * T00 ;
  71 : 
  72 :     Muh2  =(((TL** A2) * A0) / ( TL+ A10))           ;
  73 :     Muhe  =(((TL** B2) * B0) / ( TL+ B10))           ;
  74 :     Muo2  =(((TL** C2) * C0) / ( TL+ C10))           ;
  75 :     Mun2  =(((TL** D2) * D0) / ( TL+ D10))           ;
  76 :     Muco2 =(((TL** E2) * E0) / ( TL+ E10))           ;
  77 :     Muco  =(((TL** F2) * F0) / ( TL+ F10))           ;
  78 :     Muvap =(((TL** G2) * G0) / ( TL+ G10))           ;
  79 :     Muair = (0.8*Mun2) + (0.2*Muo2) ;
  80 :     'RESPROC' Muh2 Muhe Muo2 Mun2 Muco2 Muco Muvap Muair ;
  81 :     'QUITTER' CALMU ;
  82 : 
  83 : 'FINSI' ;
  84 : 'ERRE' 5 ;
  85 : FINPROC                ;
  86 :  
  87 :