Télécharger flam1.eso

Retour à la liste

Numérotation des lignes :

  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)
  25. CALL XCV(T(I),CVH2,CVN2,CVO2,CVH2O,NCVOPT)
  26. CV(I) = YH2(I)*CVH2 + YO2(I)*CVO2 + YN2(I)*CVN2 + YH2O(I)*CVH2O
  27. R(I) = YH2(I)*RRH2 + YO2(I)*RRO2 + YN2(I)*RRN2 + YH2O(I)*RRH2O
  28. XH2 = YH2(I)/MH2
  29. XO2 = YO2(I)/MO2
  30. RH2 = YH2(I)*RHO(I)
  31. RO2 = YO2(I)*RHO(I)
  32. RH2O = YH2O(I)*RHO(I)
  33. IF ( XH2 .GT. 2.0D0*XO2 ) THEN
  34. ROPREC = RO2
  35. RO2 = RO2/(1.0D0 + DT*Cf*RH2/MH2*EXP(-Ta/T(I)))
  36. DRHO = RO2-ROPREC
  37. RH2 = RH2 + 2.0D0*MH2/MO2*DRHO
  38. RH2O = RH2O - 2.0D0*MH2O/MO2*DRHO
  39. Q(I) = 2.0D0*MH2O/MO2*DRHO*ho_H2O/DT
  40. ELSE
  41. ROPREC = RH2
  42. RH2 = RH2/(1.0D0 + 2.0D0*DT*Cf*RO2/MO2*EXP(-Ta/T(I)))
  43. DRHO = RH2-ROPREC
  44. RO2 = RO2 + MO2/(2.0D0*MH2)*DRHO
  45. RH2O = RH2O - MH2O/MH2*DRHO
  46. Q(I) = MH2O/MH2*DRHO*ho_H2O/DT
  47. ENDIF
  48. ZH2(I) = RH2/RHO(I)
  49. ZO2(I) = RO2/RHO(I)
  50. ZH2O(I) = RH2O/RHO(I)
  51. ZN2(I) = 1.D0 - ZH2(I) - ZO2(I) - ZH2O(I)
  52. ENDDO
  53. C
  54. RETURN
  55. END
  56.  
  57.  
  58.  

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