Télécharger la plaquette Cast3M


fr en

Télécharger turb6.eso

Retour à la liste

Numérotation des lignes : 

turb6

  1. C TURB6     SOURCE    CHAT      05/01/13    03:57:03     5004
  2. C***********************************************************************
  3. C*                                                                     *
  4. C* PROJET      : Opérateur LIMI                                        *
  5. C* NOM         : TURB6                                                 *
  6. C* DESCRIPTION : Ce sous-programme calcule les relations de fermeture  *
  7. C*               des méthodes intégrales pour les couches limites      *
  8. C*               turbulentes. Utilisation des formules de Head et      *
  9. C*               Ludwieg&Tillmann.                                     *
  10. C* LANGAGE     : Esope                                                 *
  11. C* AUTEUR      : Guillaume VENCO - DRN/DMT/SEMT/LTMF                   *
  12. C*                                                                     *
  13. C***********************************************************************
  14. C*                                                                     *
  15. C* APPELES     : Aucun                                                 *
  16. C*                                                                     *
  17. C***********************************************************************
  18. C*                                                                     *
  19. C* ENTREES     : D2N : -Epaisseur de quantité de mouvement,            *
  20. C*                      correspond à l'inconnue de l'équation intégrale*
  21. C*                      de quantité de mouvement                       *
  22. C*               Y1N : -Différence entre delta et delta1,              *
  23. C*                      correspond à l'inconnue de l'équation intégrale*
  24. C*                      d'entrainement                                 *
  25. C*               UEN : -Champ de vitesse à l'extérieur de la couche    *
  26. C*                      limite                                         *
  27. C*               N   : -Nombre de noeuds du maillage                   *
  28. C*                                                                     *
  29. C* SORTIES     : HH  : -Facteur H*, correspond au rapport Y1N/D2N      *
  30. C*               HN  : -Facteur de forme                               *
  31. C*               CFN : -Coefficient de frottement                      *
  32. C*               CEN : -Coefficient d'entraînement                     *
  33. C*               D1N : -Epaisseur de déplacement                       *
  34. C*                                                                     *
  35. C***********************************************************************
  36. C*                                                                     *
  37. C* VERSION     : 20/04/2000                                            *
  38. C* CREATION    :                                                       *
  39. C*                                                                     *
  40. C***********************************************************************
  41.  
  42.       SUBROUTINE TURB6(D2N,Y1N,UEN,HH,HN,CFN,CEN,D1N,N)
  43.  
  44. C*       *********************************
  45. C*       *** Déclaration des variables ***
  46. C*       *********************************
  47.  
  48.       IMPLICIT INTEGER(I-N)
  49.       IMPLICIT REAL*8 (A-H,O-Z)
  50.       REAL*8 D2N(N),Y1N(N),UEN(N)
  51.       REAL*8 HH(N),HN(N),CFN(N),CEN(N),D1N(N),NU
  52.  
  53. C*       *********************************
  54. C*       ******** Sous-Programme *********
  55. C*       *********************************
  56.  
  57. C---- Viscosité du fluide
  58.       NU = 1e-6
  59.  
  60.  
  61.       DO 10 I=1,N
  62. c         write(6,*) 'Ue(',I,')=',UEN(I)
  63. C---- facteur H*
  64.          HH(I)=Y1N(I)*(D2N(I)**(-1))
  65. c        write(6,*) 'HH(',I,')=',HH(I)
  66. C---- Coefficient d'entrainement
  67.          CEN(I)=0.0306*((HH(I)-3.)**(-0.653))
  68. c         write(6,*) 'CEN(',I,')=',CEN(I)
  69. C---- Détermination du facteur de forme H
  70.          HN(I)=0.7+(((1./1.535)*(HH(I)-3.3))**(-1./2.715))
  71. c         write(6,*) 'H(',I,')=',HN(I)
  72. C---- Coefficient de frottement
  73.          CFN(I)=0.246*(10**(-0.678*HN(I)))
  74.      &         *((UEN(I)*D2N(I)/NU)**(-0.268))
  75. c         write(6,*) 'CFN(',I,')=',CFN(I)
  76. C---- Calcul de D1N
  77.          D1N(I) = HN(I)*D2N(I)
  78.  10   CONTINUE
  79.  
  80.  
  81.       RETURN
  82.       END
  83.  
  84. C*       *********************************
  85. C*       ****  Fin Sous-Programme  *******
  86. C*       *********************************
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  

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