Télécharger la plaquette Cast3M


fr en

Télécharger turb5.eso

Retour à la liste

Numérotation des lignes : 

turb5

  1. C TURB5     SOURCE    CHAT      05/01/13    03:56:59     5004
  2. C***********************************************************************
  3. C*                                                                     *
  4. C* PROJET      : Opérateur LIMI                                        *
  5. C* NOM         : TURB5                                                 *
  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 Cousteix.    *
  9. C* LANGAGE     : Esope                                                 *
  10. C* AUTEUR      : Guillaume VENCO - DRN/DMT/SEMT/LTMF                   *
  11. C*                                                                     *
  12. C***********************************************************************
  13. C*                                                                     *
  14. C* APPELES     : lEBAl2.eso, sous-programme                            *
  15. C*               RACFGA.eso, sous-programme                            *
  16. C*               JACOB.eso,  sous-programme                            *
  17. C*               equa1 , fonction                                      *
  18. C*               equa2 , fonction                                      *
  19. C*               equa3 , fonction                                      *
  20. C*                                                                     *
  21. C***********************************************************************
  22. C*                                                                     *
  23. C* ENTREES     : D2N : -Epaisseur de quantité de mouvement,            *
  24. C*                      correspond à l'inconnue de l'équation intégrale*
  25. C*                      de quantité de mouvement                       *
  26. C*               Y1N : -Différence entre delta et delta1,              *
  27. C*                      correspond à l'inconnue de l'équation intégrale*
  28. C*                      d'entrainement                                 *
  29. C*               UEN : -Champ de vitesse à l'extérieur de la couche    *
  30. C*                      limite                                         *
  31. C*               N   : -Nombre de noeuds du maillage                   *
  32. C*                                                                     *
  33. C* SORTIES     : HH  : -Facteur H*, correspond au rapport Y1N/D2N      *
  34. C*               HN  : -Facteur de forme                               *
  35. C*               CFN : -Coefficient de frottement                      *
  36. C*               CEN : -Coefficient d'entraînement                     *
  37. C*               D1N : -Epaisseur de déplacement                       *
  38. C*               DeN : -Epaisseur de couche limite                     *
  39. C*                                                                     *
  40. C***********************************************************************
  41. C*                                                                     *
  42. C* VERSION     : 20/04/2000                                            *
  43. C* CREATION    :                                                       *
  44. C*                                                                     *
  45. C***********************************************************************
  46.  
  47.       SUBROUTINE TURB5(D2N,Y1N,UEN,HH,HN,CFN,CEN,D1N,DeN,N)
  48.  
  49. C*       *********************************
  50. C*       *** Déclaration des variables ***
  51. C*       *********************************
  52.  
  53.       IMPLICIT INTEGER(I-N)
  54.       IMPLICIT REAL*8 (A-H,O-Z)
  55.       REAL*8 D2N(N),Y1N(N),UEN(N)
  56.       REAL*8 HH(N),HN(N),CFN(N),CEN(N),D1N(N),DeN(N)
  57.       REAL*8 F,P,NU,GAN,GN,A0,B0,be
  58.       REAL*8 WX(3,3),QX(3,3)
  59.       REAL*8 RD1,X1,X2,X3,XP1,XP2,XP3,f1P,f2P,f3P,r1,r2,r3,R
  60.       PARAMETER (eps=1.d-5)
  61.  
  62.  
  63. C*       *********************************
  64. C*       ******** Sous-programme *********
  65. C*       *********************************
  66.  
  67. C---- Viscosité du fluide
  68.       NU = 1e-6
  69.  
  70.       DO 10 I=1,N
  71. c         write(6,*) 'Ue(',I,')=',UEN(I)
  72. C---- Facteur H*
  73.          HH(I)=Y1N(I)*(D2N(I)**(-1))
  74. c         write(6,*) 'HH(',I,')=',HH(I)
  75.  
  76. C---- Calcul des coefficients H,gamma et G par itérations
  77. C---- Méthode de Newton
  78.  
  79. C---- Valeurs initiales calculées à partir de la
  80. C---- formule approchée
  81.          A0=1
  82.          B0=2
  83.          be=0.613
  84.          CALL LEBAL2(HH(I),be,HN(I),A0,B0)
  85.          D1N(I)=D2N(I)*HN(I)
  86.          A0=1.D-3
  87.          B0= 10
  88.          CALL RACFGA(UEN(I),D1N(I),HN(I),GAN,A0,B0)
  89.          GN=(HN(I)-1)/(HN(I)*GAN)
  90.          RD1=UEN(I)*D1N(I)/NU
  91. c         write(6,*) 'D1N(',I,')=',D1N(I)
  92. c         write(6,*) 'H0(',I,')=',HN(I)
  93. c         write(6,*) 'GN0(',I,')=',GN
  94. c         write(6,*) 'Gamma0(',I,')=',GAN
  95.  
  96. C---- Calcul des coefficients réels par itérations
  97.          X1=HN(I)
  98.          X2=GN
  99.          X3=GAN
  100.          K=0
  101.  40      CALL JACOB(X1,X2,X3,WX,QX)
  102.          K=K+1
  103. c         write(6,*) 'it=',K
  104.          f1P= equa1(X1,X2,HH(I))
  105.          f2P= equa2(X2,X3,RD1)
  106.          f3P= equa3(X1,X2,X3)
  107.  
  108.          XP1=X1-(QX(1,1)*f1P+QX(1,2)*f2P+QX(1,3)*f3P)
  109.          XP2=X2-(QX(2,1)*f1P+QX(2,2)*f2P+QX(2,3)*f3P)
  110.          XP3=X3-(QX(3,1)*f1P+QX(3,2)*f2P+QX(3,3)*f3P)
  111. c         write(6,*) 'X1(P+1)=',XP1
  112. c         write(6,*) 'X2(P+1)=',XP2
  113. c         write(6,*) 'X3(P+1)=',XP3
  114.          r1=abs(XP1-X1)
  115.          r2=abs(XP2-X2)
  116.          r3=abs(XP3-X3)
  117.  
  118.          R=max(r1,r2,r3)
  119.          IF (R.GE.eps) GOTO 20
  120.          GOTO 30
  121.  
  122.  20      X1=XP1
  123.          X2=XP2
  124.          X3=XP3
  125.          GOTO 40
  126.  
  127.  30      HN(I)=XP1
  128.          GN=XP2
  129.          GAN=XP3
  130. c         write(6,*) 'H',HN(I)
  131. c         write(6,*) 'gamma=',GAN
  132. c         write(6,*) 'G=',GN
  133.  
  134. C---- Valeurs des formules de similitude F(G) et P(G)
  135.          P=0.074*GN-1.0957*(GN**(-1))
  136.          F=0.613*GN-((3.6+76.86*((GN**(-1))-0.154)**2)*(GN**(-1)))
  137.  
  138. C---- Coefficient de frottement
  139.          CFN(I)=2*(GAN**2)
  140. c         write(6,*) 'Cf=',CFN(I)
  141.  
  142. C---- Coefficient d'entrainement Ce
  143.          CEN(I) = P*GAN
  144. c         write(6,*) 'CE=',CEN(I)
  145.  
  146. C---- Epaisseur de couche limite
  147.          DeN(I)=D1N(I)*((F*GAN)**(-1))
  148. c         write(6,*) 'delta=',DeN(I)
  149.  
  150.  10   CONTINUE
  151.  
  152.          RETURN
  153.          END
  154. C*       *********************************
  155. C*       ****  Fin Sous-Programme  *******
  156. C*       *********************************
  157.  
  158.  
  159.  
  160.  
  161.  
  162.  
  163.  
  164.  

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