Télécharger fla008.eso

Retour à la liste

Numérotation des lignes :

  1. C FLA008 SOURCE SERRE 12/07/10 21:15:13 7431
  2. FUNCTION FLA008(P,X,A1,B1,XH2OFF,BETA)
  3. C----------------------------------------------------------------------
  4. C This subroutine evaluate the consumption hydrogen rate (in kg/s)
  5. C for a passive autocatalytic recombiner; a Siemens law type is used.
  6. C----------------------------------------------------------------------
  7. C
  8. C---------------------------
  9. C Paramètres Entrée/Sortie :
  10. C---------------------------
  11. C
  12. C E/ P : Averaged containment pressure (Pa)
  13. C E/ X(7) : Molar fraction (1=N2 ,2=O2 ,3=H2 ,4=H2O)
  14. C E/ A1 : First Siemens like law coefficient
  15. C E/ B1 : Secund Siemens-like law coefficient
  16. C E/ XH2OFF : H2 molar fraction to switch off the PAR
  17. C E/ BETA : Link from the low Xh2 to the high Xh2
  18. C /S FLA008 : Consumption hydrogen rate (in kg/s)
  19. C
  20. C---------------------------------------------------------------------
  21. C
  22. C Langage : FORTRAN 77
  23. C
  24. C Mise en oeuvre : H. Paillère (1997, TTMF)
  25. C 99 : Siemens law modification : relaxation to zero with tanh
  26. C (old one: fla008=bobo*(0.308D-8*p+0.406D-3))
  27. C 03 : Modification of the linear pressure evolution and
  28. C coefficients are now user data. Default values change to
  29. C (0.48D-8*p+0.58D-3) for a Siemens FR90/1-150 PAR.
  30. C 05 : Positivity is assumed even we test it in fla010 and fla019
  31. C
  32. C---------------------------------------------------------------------
  33. IMPLICIT INTEGER(I-N)
  34. IMPLICIT REAL*8 (A-H,O-Z)
  35. REAL*8 X(7)
  36. C
  37. C Threshold value for the Siemens Law (bobo=min(Xh2,2Xo2,0.08))
  38. BOBO = 2.D0 * X(2)
  39. IF (X(3) .LT. BOBO) THEN
  40. BOBO = X(3)
  41. ENDIF
  42. IF (0.08D0 .LT. BOBO) THEN
  43. BOBO = 0.08D0
  44. ENDIF
  45. C
  46. C Siemens Law
  47. SIE008 = (A1*P+B1) * BOBO * TANH((20.D0/BETA)*(BOBO-XH2OFF))
  48. IF (SIE008 .LT. 0.D0) THEN
  49. FLA008 = 0.D0
  50. ELSE
  51. FLA008 = SIE008
  52. ENDIF
  53. C
  54. RETURN
  55. END
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.  
  64.  
  65.  
  66.  
  67.  

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