Procedures | Cast3M

   1 : * FPAL      PROCEDUR  MAGN      98/09/26    21:15:03     3310           
   2 : *X FPAL (Procedure)
   3 : *
   4 : *  procedure  FPAL
   5 : *  ---------------
   6 : *
   7 : *    objet : calcule la vitesse de depot AK de particules microniques
   8 : *            en ecoulement laminaire
   9 : *
  10 : *            AK = (flux depose sur la paroi)/(concentration a l'infini)
  11 : *
  12 : *
  13 : *  SYNTAXE :    AK = FPAL  NU ROF UET NOR ROG RAP $PAROI ;
  14 : *  -------
  15 : *
  16 : *   AK    CHPOINT SCAL CENTRE PAROI  vitesse de depot (m/s)
  17 : *   NU    FLOTTANT                   viscosite cinematique du gaz (m2/s)
  18 : *   ROF   FLOTTANT                   masse volumique du gaz (kg/m3)
  19 : *   UET   CHPOINT SCAL CENTRE PAROI  vitesse de frottement (m/s)
  20 : *   NOR   CHPOINT VECT FACE          champ des normales aux faces (m)
  21 : *   ROG   POINT                      masse vol. part.x gravite (kg/m2s2)
  22 : *   RAP   FLOTTANT                   rayon des particules (m)
  23 : *   $PAROI MMODEL TYPE 'NAVIER_STOKES'   maillage paroi
  24 : *
  25 : *-----------------------------------------------------------------------
  26 : 
  27 : DEBPROC FPAL ;
  28 : ARGU NU*FLOTTANT ROF*FLOTTANT UET*CHPOINT NOR*CHPOINT ROG*POINT
  29 :      RAP*FLOTTANT $PAROI*MMODEL ;
  30 : 
  31 : * coefficient de correction de Cunningham pour les petites particules
  32 : * LL = libre parcors moyen des molecules de gaz
  33 : * KK = nombre de KNUDSEN
  34 : 
  35 : LL  = 6.53E-8 ;
  36 : KK  = LL/RAP ;
  37 : CUN = 1. + (1.257*KK) + (0.4*KK/(EXP(1.1/KK))) ;
  38 : 
  39 : *--------------------------------------- SEDIMENTATION -----------------
  40 : * vitesse de sedimentation incomplete
  41 : 
  42 : CRN  = CUN/(ROF*NU)      ;
  43 : VSI  = 2.*CRN*RAP*RAP/9. ;
  44 : 
  45 : * produit scalaire verticale.normale
  46 : 
  47 : NORP = KCHT $PAROI VECT CENTRE NOR ;
  48 : VSI2 = KCHT $PAROI VECT CENTRE ROG ;
  49 : 
  50 : AKI = PSCA  NORP VSI2 (MOTS UX UY) (MOTS UX UY) ;
  51 : 
  52 : * vitesse de depot par sedimentation
  53 : 
  54 : AKS = VSI*AKI ;
  55 : 
  56 : *--------------------------------------- DIFFUSION ---------------------
  57 : * coefficient de diffusion brownienne
  58 : 
  59 : DIF = 2.145E-22*CRN/RAP ;
  60 : 
  61 : * nombre de SCHMIDT
  62 : 
  63 : SC = NU/DIF ;
  64 : 
  65 : * epaisseur de la couche limite massique adimensionnelle (5