Télécharger fla016.eso

Retour à la liste

Numérotation des lignes :

  1. C FLA016 SOURCE TTMF3 12/07/05 21:15:38 7425
  2. FUNCTION FLA016(X,ROE,A,B,C,DD,EE,FF)
  3. C------------------------------------------------------------------------
  4. C Calcul de f(x) dans le cas où le débit n'est pas limité par les
  5. C conditions d'entrée : recherche de la racine par dichotomie
  6. C (voir fla015.eso)
  7. C------------------------------------------------------------------------
  8. C
  9. C---------------------------
  10. C Parametres Entree/Sortie :
  11. C---------------------------
  12. C
  13. C E/ X : flottant : Debit au travers du PAR (kg/s)
  14. C E/ ROE : flottant : Densité du mélange à l'entrée du PAR (kg/m3)
  15. C E/ A : flottant : Coefficient devant X^2
  16. C E/ B : flottant : Coefficient devant X
  17. C E/ DD : flottant : Coefficient
  18. C E/ EE : flottant : Coefficient
  19. C E/ FF : flottant : Coefficient
  20. C /S FLA016 : flottant : Valeur de la fonction f(X)
  21. C
  22. C------------------------------------------------------------------------
  23. C
  24. C Langage : FORTRAN 77
  25. C
  26. C Mise en oeuvre : H. Paillère (1997, TTMF)
  27. C
  28. C------------------------------------------------------------------------
  29. IMPLICIT INTEGER(I-N)
  30. IMPLICIT REAL*8 (A-H,O-Z)
  31. C
  32. C Polynome de degré 6 tirée de Bernouilli dont la plus grande racine au
  33. C dela des poles (EE-DD) et FF, si elle existe, est le débit cherché
  34. P61 = A*X*X + B*X + C*ROE*ROE
  35. P62 = (X-EE+DD)*(X-EE+DD)*(X-FF)*(X-FF)
  36. P63 = C*ROE*ROE*(X-EE)*(X-EE)*X*X
  37. C
  38. FLA016 = P61*P62 - P63
  39. C
  40. RETURN
  41. END
  42.  
  43.  
  44.  
  45.  
  46.  
  47.  
  48.  

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