Esope source files

   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 :