Esope source files

   1 : C FLAM1     SOURCE    CHAT      05/01/13    00:03:07     5004           
   2 :       SUBROUTINE FLAM1(N,RHO,CV,R,YH2,YO2,YN2,YH2O,T,
   3 :      &                 ZH2,ZO2,ZN2,ZH2O,Q,Dt,NCVOPT)
   4 : C
   5 : C LAMINAR COMBUSTION OF HYDROGEN
   6 : C
   7 :       IMPLICIT INTEGER(I-N)
   8 :       IMPLICIT REAL*8 (A-H,O-Z)
   9 : C
  10 :       DIMENSION RHO(*),YH2(*),YO2(*),YN2(*),YH2O(*)
  11 :       DIMENSION CV(*),T(*),Q(*),R(*)
  12 :       DIMENSION ZH2(*),ZO2(*),ZN2(*),ZH2O(*)
  13 :       REAL*8 MH2,MO2,MN2,MH2O
  14 :       DATA Rg,MH2,MO2,MN2,MH2O/8.313D0,2.0D-3,32.0D-3,28.0D-3,18.0D-3/
  15 :       DATA Cf,Ta/3.3D5,1.0D4/
  16 :       DATA RRH2,RRO2,RRN2,RRH2O/4156.5D0,259.8D0,296.9D0,461.8D0/
  17 :       DATA ho_H2O/-1.327D7/
  18 : C
  19 :         DO I=1,N
  20 : C
  21 :         YH2(I) = ABS(YH2(I))
  22 :         YO2(I) = ABS(YO2(I))
  23 :         YH2O(I) = ABS(YH2O(I))
  24 :         YN2(I) = 1.D0 - YH2(I) - YO2(I) - YH2O(I)
  24 :         CALL XCV(T(I),CVH2,CVN2,CVO2,CVH2O,NCVOPT)
  25 :         CV(I) = YH2(I)*CVH2 + YO2(I)*CVO2 + YN2(I)*CVN2 + YH2O(I)*CVH2O
  26 :         R(I) = YH2(I)*RRH2 + YO2(I)*RRO2 + YN2(I)*RRN2 + YH2O(I)*RRH2O
  27 :         XH2 = YH2(I)/MH2
  28 :         XO2 = YO2(I)/MO2
  29 :         RH2 = YH2(I)*RHO(I)
  30 :         RO2 = YO2(I)*RHO(I)
  31 :         RH2O = YH2O(I)*RHO(I)
  32 :         IF ( XH2 .GT. 2.0D0*XO2 ) THEN
  33 :           ROPREC = RO2
  34 :           RO2 = RO2/(1.0D0 + DT*Cf*RH2/MH2*EXP(-Ta/T(I)))
  35 :           DRHO = RO2-ROPREC
  36 :           RH2 = RH2 + 2.0D0*MH2/MO2*DRHO
  37 :           RH2O = RH2O - 2.0D0*MH2O/MO2*DRHO
  38 :           Q(I) = 2.0D0*MH2O/MO2*DRHO*ho_H2O/DT
  39 :         ELSE
  40 :           ROPREC = RH2
  41 :           RH2 = RH2/(1.0D0 + 2.0D0*DT*Cf*RO2/MO2*EXP(-Ta/T(I)))
  42 :           DRHO = RH2-ROPREC
  43 :           RO2 = RO2 + MO2/(2.0D0*MH2)*DRHO
  44 :           RH2O = RH2O - MH2O/MH2*DRHO
  45 :           Q(I) = MH2O/MH2*DRHO*ho_H2O/DT
  46 :         ENDIF
  47 :           ZH2(I) = RH2/RHO(I)
  48 :           ZO2(I) = RO2/RHO(I)
  49 :           ZH2O(I) = RH2O/RHO(I)
  50 :           ZN2(I) = 1.D0 - ZH2(I) - ZO2(I) - ZH2O(I)
  51 :         ENDDO
  52 : C
  53 :        RETURN
  54 :        END
  55 :  
  56 :