Télécharger la plaquette Cast3M


fr en

Télécharger fla016.eso

Retour à la liste

Numérotation des lignes : 

fla016

  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