Télécharger lam2.eso

Retour à la liste

Numérotation des lignes :

  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