Télécharger la plaquette Cast3M


fr en

Télécharger fla008.eso

Retour à la liste

Numérotation des lignes : 

fla008

  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