Esope source files

   1 : C LAM2      SOURCE    CHAT      05/01/13    01:12:18     5004           
   2 : C***********************************************************************
   3 : C*                                                                     *
   4 : C* PROJET      : Opérateur LIMI                                        *
   5 : C* NOM         : lAM2                                                  *
   6 : C* DESCRIPTION : Ce sous programme calcule les relations de fermeture  *
   7 : C*               des méthodes intégrales pour les couches limites      *
   8 : C*               laminaires dans le cas de la méthode à 2 équations.   *
   9 : C* LANGAGE     : Esope                                                 *
  10 : C* AUTEUR      : Guillaume VENCO - DRN/DMT/SEMT/LTMF                   *
  11 : C*                                                                     *
  12 : C***********************************************************************
  13 : C*                                                                     *
  14 : C* APPELES     : NEWT2.eso, sous-programme                             *
  15 : C*                                                                     *
  16 : C***********************************************************************
  17 : C*                                                                     *
  18 : C* ENTREES     : D2N : -Epaisseur de quantité de mouvement,            *
  19 : C*                      correspond à l'inconnue de l'équation intégrale*
  20 : C*                      de quantité de mouvement                       *
  21 : C*               D3N : -Epaisseur d'énergie cinétique.                 *
  22 : C*                      correspond à l'inconnue de l'équation intégrale*
  23 : C*                      d'énergie cinétique.                           *
  24 : C*               N   : -Nombre de noeuds du maillage                   *
  25 : C*                                                                     * 
  26 : C* SORTIES     : H32 : -Facteur H32, correspond au rapport D3N/D2N     *
  27 : C*               H   : -Facteur de forme                               *
  28 : C*               BN  : -Coefficient de la relation permettant de       *
  29 : C*                      calculer le coefficient de frottement          *
  30 : C*               DN  : -Coefficient de la relation permetteant de      *
  31 : C*                      calculer le coefficient d'énergie cinétique    *
  32 : C*                                                                     *
  33 : C***********************************************************************
  34 : C*                                                                     *
  35 : C* VERSION     : 21/04/2000                                            *
  36 : C* CREATION    :                                                       *
  37 : C*                                                                     *
  38 : C***********************************************************************
  39 : 
  40 :       SUBROUTINE LAM2(D2N,D3N,H32,H,BN,DN,D1N,N)
  41 : 
  42 : C*       *********************************
  43 : C*       *** Déclaration des variables ***
  44 : C*       *********************************
  45 : 
  46 :       IMPLICIT INTEGER(I-N)
  47 :       IMPLICIT REAL*8(A-H,O-Z)
  48 :       REAL*8 D2N(N),D3N(N)
  49 :       REAL*8 H32(N),H(N),BN(N),DN(N),D1N(N)
  50 :       REAL*8 a0,a1,a2
  51 :       REAL*8 A00,B00
  52 : 
  53 : C*       *********************************
  54 : C*       ******** Sous-programme *********
  55 : C*       *********************************
  56 : 
  57 :       DO I=1,N
  58 : 
  59 : C---- Calcul du facteur H32
  60 :          H32(I) = D3N(I)*(D2N(I)**(-1))
  61 : c         write(6,*) 'H32=',H32(I)
  62 : C---- Test pour la poursuite du calcul
  63 :          IF(H32(I).LE.1.5215d0) THEN
  64 :             write(6,*) 'ERREUR: On ne peut pas déterminer H'
  65 :          ENDIF 
  66 :      
  67 : C---- Calcul des coefficients de la relation donnant H
  68 :          a0 = -H32(I)-(50.84951/H32(I))+46.8818
  69 :          a1 = -23.78186
  70 :          a2 = (4.02923**2)
  71 :          
  72 : C---- Calcul de H en fonction de H32 par la méthode de dichotomie
  73 : C---- Recherche de la racine comprise entre les valeurs A00 et B00         
  74 :          A00=1.512
  75 :          B00=4.              
  75 :          CALL newt2(a0,a1,a2,A00,B00,H(I))
  76 : 
  77 : C---- Calcul du coefficient b(H)
  78 :          BN(I)= 2.99259*((((1./H(I))-(1./8.05846))**1.7)
  79 :      &          -((1./8.05846)**1.7))
  80 : 
  81 : C---- Calcul du coefficient d(H)
  82 :          DN(I)=H32(I)*(BN(I)-(H(I)-1)*(-0.06815+4.336355
  83 :      &          *(((1./H(I))-(1./4.02923))**2.095065)))*0.5
  84 : 
  85 : C---- Epaisseur de déplacement
  86 :          D1N(I) = H(I)*D2N(I)
  87 :       ENDDO
  88 : 
  89 :       RETURN
  90 :       END
  91 : 
  92 : C*       *********************************
  93 : C*       ****  Fin Sous-Programme  *******
  94 : C*       ********************************
  95 :  
  96 :