Télécharger la plaquette Cast3M


fr en

Télécharger fpal.procedur

Retour à la liste

Numérotation des lignes : 

  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<Cste<20)
  66.  
  67. YPLUS =  5./(SC**0.333) ;
  68.  
  69. * epaisseur de la couche limite massique dimensionnee
  70.  
  71. UET1 = UET**(-1.)    ;
  72. YC = (YPLUS*NU)*UET1 ;
  73.  
  74. * vitesse de depot par diffusion
  75.  
  76. YC1 = YC**(-1.) ;
  77. AKD = DIF*YC1 ;
  78.  
  79. *-------------------------- VITESSE DE DEPOT TOTALE --------------------
  80.  
  81. AK0 = AKS + AKD ;
  82. AK = KCHT $PAROI SCAL CENTRE  AK0 ;
  83.  
  84. FINPROC AK ;
  85.  
  86.  
  87.  
  88.  

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