Télécharger la plaquette Cast3M


fr en

Télécharger flama1.eso

Retour à la liste

Numérotation des lignes : 

flama1

  1. C FLAMA1    SOURCE    CHAT      05/01/13    00:03:14     5004
  2.       SUBROUTINE FLAMA1(NC,TSEUIL,H0H2,H0O2,H0H2O,TEMP,
  3.      &                  RYH2, RYO2, RYH2O,
  4.      &                  DE, NRYH2, NRYO2, NRYH2O,
  5.      &                  LOGNEG, MESERR, VALER)
  6. C************************************************************************
  7. C
  8. C PROJET            :  CASTEM 2000
  9. C
  10. C NOM               :  FLAMA1
  11. C
  12. C DESCRIPTION       :  voir  FLAMA
  13. C
  14. C                      Calculs
  15. C
  16. C LANGAGE           :  FORTRAN 77 (avec extensions F90)
  17. C
  18. C AUTEUR            :  A. BECCANTINI, DRN/DMT/SEMT/LTMF
  19. C
  20. C************************************************************************
  21. C
  22. C
  23. C APPELES (E/S)     :  /
  24. C
  25. C APPELES (Calcul)  :  /
  26. C
  27. C
  28. C************************************************************************
  29. C
  30. C ENTREES  :
  31. C
  32. C            NC = nombre de points centre
  33. C
  34. C            TSEUIL(NC) = champ de la temperature de seuil
  35. C
  36. C            H0H2, H0O2, HOH2O = enthalpie de formation a 0K (J/kg)
  37. C
  38. C            TEMP(NC)  = champ de temperature
  39. C
  40. C            RYH2(NC)  = masse volumique d'hydrogene (initial)
  41. C
  42. C            RYO2(NC)  = masse volumique d'oxygene
  43. C
  44. C            RYH2O(NC) = masse volumique d'eau
  45. C
  46. C SORTIES  :
  47. C
  48. C            DE(NC) = variation d'energie thermique a suite de la reaction
  49. C
  50. C            NRYH2(NC)  = masse volumique d'hydrogene (final)
  51. C
  52. C            NRYO2(NC)  = masse volumique d'oxygene
  53. C
  54. C            NRYH2O(NC) = masse volumique d'eau
  55. C
  56. C************************************************************************
  57. C
  58. C HISTORIQUE (Anomalies et modifications éventuelles)
  59. C
  60. C HISTORIQUE :  Créée le 1.3.99
  61. C
  62. C************************************************************************
  63. C
  64. C
  65. C**** Les variables
  66. C
  67.       IMPLICIT INTEGER(I-N)
  68.       INTEGER I1, NC
  69.       REAL*8 MH2,MO2,MH2O, H0H2, H0O2, H0H2O, VALER
  70.      &       ,TSEUIL(NC),RYH2(NC),RYO2(NC),RYH2O(NC), DE(NC)
  71.      &       ,NRYH2(NC),NRYO2(NC),NRYH2O(NC), TEMP(NC)
  72.      &       ,TEMP0, TSEU0, RYH2I, RYO2I, RYH2OI
  73.      &       ,XH2, XO2, XH2O, XMAX, CELL, DRYH2, DRYO2, DRYH2O, DMOL
  74.       CHARACTER*(8) MESERR
  75.       LOGICAL LOGNEG
  76.       PARAMETER(MH2=2.016D-3,MO2=31.998D-3,MH2O=18.015D-3)
  77. C
  78.       DO I1 = 1, NC
  79.          TEMP0 = TEMP(I1)
  80.          TSEU0 = TSEUIL(I1)
  81.          RYH2I = RYH2(I1)
  82.          RYO2I = RYO2(I1)
  83.          RYH2OI = RYH2O(I1)
  84.          IF(RYH2I .LT. 0.0D0)THEN
  85.             LOGNEG = .TRUE.
  86.             VALER = RYH2I
  87.             MESERR = 'RYH2    '
  88.             GOTO 9999
  89.          ENDIF
  90.          IF(RYO2I .LT. 0.0D0)THEN
  91.             LOGNEG = .TRUE.
  92.             VALER = RYO2I
  93.             MESERR = 'RYO2    '
  94.             GOTO 9999
  95.          ENDIF
  96.          IF(RYH2OI .LT. 0.0D0)THEN
  97.             LOGNEG = .TRUE.
  98.             VALER = RYH2OI
  99.             MESERR = 'RYH2O   '
  100.             GOTO 9999
  101.          ENDIF
  102.          IF(TEMP0 .GE. TSEU0)THEN
  103. C
  104. C********** A la recherche de l'espece majortitaire
  105. C
  106. C           XH2 = mole / m^3 d'H2
  107. C
  108.             XH2 = RYH2I / MH2
  109.             XO2 = RYO2I / MO2
  110.             XH2O = RYH2OI / MH2O
  111.             XMAX = MAX(XH2, XO2, XH2O)
  112.             IF(XH2 .GE. (2.0D0 * XO2))THEN
  113. C
  114. C           H2 majoritaire -> O2 brule completement
  115. C
  116.                CELL = MIN(XO2,1.0D-8*XMAX) * MO2
  117.                DRYO2 = CELL - RYO2I
  118.                DMOL = DRYO2 / MO2
  119.                DRYH2 = 2.0D0 * DMOL * MH2
  120.                DRYH2O = -2.0D0 * DMOL * MH2O
  121.                NRYO2(I1) = CELL
  122.                NRYH2(I1) = RYH2I + DRYH2
  123.                NRYH2O(I1) = RYH2OI + DRYH2O
  124.             ELSE
  125. C
  126. C           O2 majoritaire -> H2 brule completement
  127. C
  128.                CELL = MIN(XH2,1.0D-8*XMAX) * MH2
  129.                DRYH2 = CELL - RYH2I
  130.                DMOL = DRYH2 / MH2 / 2.0D0
  131.                DRYO2 = DMOL * MO2
  132.                DRYH2O = -2.0D0 * DMOL * MH2O
  133.                NRYH2(I1) = CELL
  134.                NRYO2(I1) = RYO2I + DRYO2
  135.                NRYH2O(I1) = RYH2OI + DRYH2O
  136.             ENDIF
  137.             DE(I1) = (DRYH2 * H0H2) + (DRYO2 * H0O2) + (DRYH2O * H0H2O)
  138.             DE(I1) = -1.0D0 * DE(I1)
  139.          ENDIF
  140.       ENDDO
  141. C
  142.  9999 CONTINUE
  143. C
  144.       RETURN
  145.       END
  146.  
  147.  
  148.  
  149.  
  150.  
  151.  
  152.  
  153.  
  154.  

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