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