Télécharger equa1.eso

Retour à la liste

Numérotation des lignes :

  1. C EQUA1 SOURCE CHAT 05/01/12 23:44:15 5004
  2. C***********************************************************************
  3. C* *
  4. C* PROJET : Opérateur CLMI *
  5. C* NOM : equa1, equa2,equa3 *
  6. C* DESCRIPTION : Dans ce fichier, il s'agit de trois fonctions qui *
  7. C* calculent les valeurs des équations du système des *
  8. C* relations de fermeture de Cousteix *
  9. C* LANGAGE : Esope *
  10. C* AUTEUR : Guillaume VENCO - DRN/DMT/SEMT/LTMF *
  11. C* *
  12. C***********************************************************************
  13. C* *
  14. C* APPELES : Aucun *
  15. C* *
  16. C***********************************************************************
  17. C* *
  18. C* ENTREES : H : -Facteur de forme *
  19. C* G : -Paramètre de Clauser *
  20. C* HH : -Facteur H*,c orrespond au rapport Y1N/D2N *
  21. C* limite *
  22. C* ga : -Facteur de frottement gamma *
  23. C* RD1 : -Nombre de Reynolds calculé sur l'épaisseur *
  24. C* d'entrainement *
  25. C* *
  26. C* SORTIES : equa1 : -Valeur de la relation 1 *
  27. C* equa2 : -Valeur de la relation 2 *
  28. C* equa3 : -Valeur de la relation 3 *
  29. C* *
  30. C***********************************************************************
  31. C* *
  32. C* VERSION : 20/04/2000 *
  33. C* CREATION : *
  34. C* *
  35. C***********************************************************************
  36.  
  37. C* *********************************
  38. C* ********** EQUATION 1 ***********
  39. C* *********************************
  40.  
  41. DOUBLE PRECISION FUNCTION equa1(H,G,HH)
  42. IMPLICIT INTEGER(I-N)
  43. IMPLICIT REAL*8 (A-H,O-Z)
  44. REAL*8 H,G,HH,f,beta
  45. f=0.613*G-(3.6+76.86*(((1./G)-0.154)**2))/G
  46. beta=(G/f)-1
  47. equa1= HH-(beta*(H**2)+H)/(H-1)
  48. RETURN
  49. END
  50.  
  51.  
  52.  

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