Télécharger la plaquette Cast3M


fr en

Télécharger lam2.eso

Retour à la liste

Numérotation des lignes : 

lam2

  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.
  76.          CALL newt2(a0,a1,a2,A00,B00,H(I))
  77.  
  78. C---- Calcul du coefficient b(H)
  79.          BN(I)= 2.99259*((((1./H(I))-(1./8.05846))**1.7)
  80.      &          -((1./8.05846)**1.7))
  81.  
  82. C---- Calcul du coefficient d(H)
  83.          DN(I)=H32(I)*(BN(I)-(H(I)-1)*(-0.06815+4.336355
  84.      &          *(((1./H(I))-(1./4.02923))**2.095065)))*0.5
  85.  
  86. C---- Epaisseur de déplacement
  87.          D1N(I) = H(I)*D2N(I)
  88.       ENDDO
  89.  
  90.       RETURN
  91.       END
  92.  
  93. C*       *********************************
  94. C*       ****  Fin Sous-Programme  *******
  95. C*       ********************************
  96.  
  97.  
  98.  

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