Télécharger la plaquette Cast3M


fr en

Télécharger kuet.procedur

Retour à la liste

Numérotation des lignes : 

  1. * KUET      PROCEDUR  MAGN      98/09/26    21:16:48     3310           
  2. *X KUET (Procedure)
  3. *
  4. *    procedure  KUET
  5. *    ---------------
  6. *
  7. *    calcule la vitesse de frottement a la paroi U* pour un ecoulement
  8. *  laminaire
  9. *             U* = SQRT (TAU/RO) = SQRT (NU*dV/dY)
  10. *
  11. *      TAU  cisaillement
  12. *      RO   masse volumique
  13. *      NU   viscosite cinematique
  14. *      V    vitesse parallele a la paroi
  15. *      Y    distance a la paroi
  16. *
  17. *
  18. *    SYNTAXE :  UET = KUET  NU UN NOR $DOMTOT $PAROI ;
  19. *    --------
  20. *       UET   CHPOINT SCAL CENTRE PAROI  vit. de frottement
  21. *       NU    FLOTTANT                   viscosite cinematique
  22. *       UN    CHPOINT VECT SOMMET        champ de vitesse
  23. *       NOR   CHPOINT VECT FACE          champ des normales aux faces
  24. *    $DOMTOT  MAILLAGE                   table domaine total
  25. *    $PAROI   MAILLAGE              table domaine paroi support de UET
  26. *
  27. *-----------------------------------------------------------------------
  28.  
  29. DEBPROC KUET ;
  30. ARGU NU*FLOTTANT UN*CHPOINT NOR*CHPOINT $DOMTOT*MMODEL $PAROI*MMODEL ;
  31.  
  32. VXM = MAXI  VX  'ABS' ;
  33. VYM = MAXI  VY  'ABS' ;
  34.  
  35. * cas particulier où V // (OY) : VX = 0
  36. SI ( VXM EGA 0. 1.E-6 ) ;
  37.  
  38.   GVY = KOPS  VY 'GRAD' $DOMTOT ;
  39.   GXV = EXCO  UX GVY  ;
  40.   GVY = KCHT  $DOMTOT SCAL CENTRE GXV ;
  41.   GVS = ELNO  $DOMTOT GVY ;
  42.   GVP = KCHT  $PAROI  SCAL SOMMET GVS ;
  43.  
  44. SINON ;
  45.   SI ( VYM EGA 0. 1.E-6 ) ;
  46. * cas particulier où V // (OX) : VY = 0
  47.  
  48.     GVX = KOPS  VX 'GRAD' $DOMTOT ;
  49.     GYV = EXCO  UY GVX  ;
  50.     GVY = KCHT  $DOMTOT SCAL CENTRE GYV ;
  51.     GVS = ELNO  $DOMTOT GVY ;
  52.     GVP = KCHT  $PAROI  SCAL SOMMET GVS ;
  53.  
  54.   SINON ;
  55. * cas general où  NON (VX=0 OU VY=0 )
  56.  
  57.     NV = PSCA  VN VN (MOTS UX UY) (MOTS UX UY) ;
  58.     NV = NV**0.5 ;
  59.     NV = KCHT  $DOMTOT SCAL SOMMET NV ;
  60.     GVC = KOPS NV 'GRAD' $DOMTOT ;
  61.     GVS = ELNO  $DOMTOT GVC ;
  62.     GVVP = KCHT  $PAROI VECT SOMMET GVS ;
  63.     NORP = KCHT  $PAROI VECT CENTRE NOR ;
  64.     NORS = ELNO  $PAROI NORP ;
  65.     GVP  = PSCA  NORS GVVP (MOTS UX UY) (MOTS UX UY) ;
  66.  
  67.   FINSI ;
  68. FINSI ;
  69.  
  70. DVDY = ABS  GVP ;
  71. UET0 = (NU*DVDY)**0.5 ;
  72. UET1 = KCHT $PAROI SCAL SOMMET UET0 ;
  73. UET  = NOEL $PAROI UET1 ;
  74.  
  75. FINPROC UET ;
  76.  
  77.  
  78.  
  79.  

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