C TURB6     SOURCE    CHAT      05/01/13    03:57:03     5004
C***********************************************************************
C*                                                                     *
C* PROJET      : Opérateur LIMI                                        *
C* NOM         : TURB6                                                 *
C* DESCRIPTION : Ce sous-programme calcule les relations de fermeture  *
C*               des méthodes intégrales pour les couches limites      *
C*               turbulentes. Utilisation des formules de Head et      *
C*               Ludwieg&Tillmann.                                     *
C* LANGAGE     : Esope                                                 *
C* AUTEUR      : Guillaume VENCO - DRN/DMT/SEMT/LTMF                   *
C*                                                                     *
C***********************************************************************
C*                                                                     *
C* APPELES     : Aucun                                                 *
C*                                                                     *
C***********************************************************************
C*                                                                     *
C* ENTREES     : D2N : -Epaisseur de quantité de mouvement,            *
C*                      correspond à l'inconnue de l'équation intégrale*
C*                      de quantité de mouvement                       *
C*               Y1N : -Différence entre delta et delta1,              *
C*                      correspond à l'inconnue de l'équation intégrale*
C*                      d'entrainement                                 *
C*               UEN : -Champ de vitesse à l'extérieur de la couche    *
C*                      limite                                         *
C*               N   : -Nombre de noeuds du maillage                   *
C*                                                                     *
C* SORTIES     : HH  : -Facteur H*, correspond au rapport Y1N/D2N      *
C*               HN  : -Facteur de forme                               *
C*               CFN : -Coefficient de frottement                      *
C*               CEN : -Coefficient d'entraînement                     *
C*               D1N : -Epaisseur de déplacement                       *
C*                                                                     *
C***********************************************************************
C*                                                                     *
C* VERSION     : 20/04/2000                                            *
C* CREATION    :                                                       *
C*                                                                     *
C***********************************************************************

      SUBROUTINE TURB6(D2N,Y1N,UEN,HH,HN,CFN,CEN,D1N,N)

C*       *********************************
C*       *** Déclaration des variables ***
C*       *********************************

      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
      REAL*8 D2N(N),Y1N(N),UEN(N)
      REAL*8 HH(N),HN(N),CFN(N),CEN(N),D1N(N),NU

C*       *********************************
C*       ******** Sous-Programme *********
C*       *********************************

C---- Viscosité du fluide
      NU = 1e-6


      DO 10 I=1,N
c         write(6,*) 'Ue(',I,')=',UEN(I)
C---- facteur H*
         HH(I)=Y1N(I)*(D2N(I)**(-1))
c        write(6,*) 'HH(',I,')=',HH(I)
C---- Coefficient d'entrainement
         CEN(I)=0.0306*((HH(I)-3.)**(-0.653))
c         write(6,*) 'CEN(',I,')=',CEN(I)
C---- Détermination du facteur de forme H
         HN(I)=0.7+(((1./1.535)*(HH(I)-3.3))**(-1./2.715))
c         write(6,*) 'H(',I,')=',HN(I)
C---- Coefficient de frottement
         CFN(I)=0.246*(10**(-0.678*HN(I)))
     &         *((UEN(I)*D2N(I)/NU)**(-0.268))
c         write(6,*) 'CFN(',I,')=',CFN(I)
C---- Calcul de D1N
         D1N(I) = HN(I)*D2N(I)
 10   CONTINUE


      RETURN
      END

C*       *********************************
C*       ****  Fin Sous-Programme  *******
C*       *********************************