Procedures | Cast3M

   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 :