Télécharger fullerx.eso

Retour à la liste

Numérotation des lignes :

fullerx
  1. C FULLERX SOURCE CB215821 20/11/25 13:29:38 10792
  2. SUBROUTINE FULLERX
  3. C
  4. C--------------------------------------------------------------------
  5. C Evaluation du coefficient de diffusion de la vapeur dans un melange
  6. C de gaz pour une pression et une temperature donnee
  7. C--------------------------------------------------------------------
  8. C
  9. C---------------------------
  10. C Phrase d'appel (GIBIANE) :
  11. C---------------------------
  12. C
  13. C CHP1 = FULLERNC TAB1 ;
  14. C
  15. C------------------------
  16. C Operandes et resultat :
  17. C------------------------
  18. C
  19. C TAB1 : TABLE SET contenant l'ensemble des data 0D
  20. C CHP1 : CHPO de support CELL contenant le coefficient de diffusion
  21. C de la vapeur dans le melange de gaz (en m2/s)
  22. C
  23. C--------------------------------------------------------------------
  24. C
  25. IMPLICIT INTEGER(I-N)
  26. IMPLICIT REAL*8 (A-H,O-Z)
  27. CHARACTER*4 NOMTOT(1)
  28. CHARACTER*8 TYPE,TYPS,MTYPI,MTYPR,MOT1,MOT2,BETA(7)
  29. CHARACTER*72 CHARR
  30. LOGICAL LOGI,LOGR
  31. C
  32.  
  33. -INC PPARAM
  34. -INC CCOPTIO
  35. -INC SMCHPOI
  36. -INC SMELEME
  37. -INC SMTABLE
  38. SEGMENT IWORK
  39. REAL*8 V(NBCONS,2),W(NBCELL,2)
  40. CHARACTER*8 BETA(NBCONS)
  41. ENDSEGMENT
  42. C
  43. C- Lecture et controles des donnees d'entree,
  44. C
  45. C MTAB1 : Table SET
  46. C MTAB2 : Table INCO
  47. C MTAB3 : Table COMPONENT
  48. C MTAB4 : Table GEOINF
  49. C
  50. CALL LIRTAB('SET',MTAB1,1,IRETOU)
  51. IF (IERR.NE.0) RETURN
  52. MTYPI = 'MOT '
  53. MTYPR = 'TABLE '
  54. CALL ACCTAB(MTAB1,MTYPI,IVALI,XVALI,'INCO',LOGI,IRETI,
  55. & MTYPR,IVALR,XVALR,CHARR,LOGR,MTAB2)
  56. IF (IERR.NE.0) RETURN
  57. CALL ACCTAB(MTAB1,MTYPI,IVALI,XVALI,'COMPONENT',LOGI,IRETI,
  58. & MTYPR,IVALR,XVALR,CHARR,LOGR,MTAB3)
  59. IF (IERR.NE.0) RETURN
  60. CALL ACCTAB(MTAB1,MTYPI,IVALI,XVALI,'GEOINF',LOGI,IRETI,
  61. & MTYPR,IVALR,XVALR,CHARR,LOGR,MTAB4)
  62. IF (IERR.NE.0) RETURN
  63. C
  64. C IPT1 : Maillage POI1 des CELLs
  65. C
  66. TYPE = 'MAILLAGE'
  67. CALL ACMO(MTAB4,'CELL',TYPE,IPT1)
  68. IF (IERR.NE.0) RETURN
  69. C
  70. C MCHPO1 : Pression
  71. C MCHPO2 : Temperature
  72. C
  73. TYPE = 'CHPOINT '
  74. TYPS = ' '
  75. CALL ACMO(MTAB2,'PT',TYPE,MCHPO1)
  76. CALL LRCHT(MCHPO1,MPOVA1,TYPS,IGEO1)
  77. IF (IERR.NE.0) RETURN
  78. IF (IGEO1.NE.IPT1) THEN
  79. INDIC = 1
  80. NBCOMP = -1
  81. CALL QUEPOI(MCHPO1,IPT1,INDIC,NBCOMP,NOMTOT)
  82. IF (IERR.NE.0) RETURN
  83. CALL LRCHT(MCHPO1,MPOVA1,TYPS,IGEO1)
  84. ENDIF
  85. C
  86. TYPE = 'CHPOINT '
  87. CALL ACMO(MTAB2,'TGAS',TYPE,MCHPO2)
  88. IF (IERR.NE.0) RETURN
  89. CALL LRCHT(MCHPO2,MPOVA2,TYPS,IGEO1)
  90. IF (IGEO1.NE.IPT1) THEN
  91. INDIC = 1
  92. NBCOMP = -1
  93. CALL QUEPOI(MCHPO2,IPT1,INDIC,NBCOMP,NOMTOT)
  94. IF (IERR.NE.0) RETURN
  95. CALL LRCHT(MCHPO2,MPOVA2,TYPS,IGEO1)
  96. ENDIF
  97. C
  98. C- Creation du chpo contenant le coefficient de diffusion
  99. C- de la vapeur dans les incondensables
  100. C
  101. SEGINI, MCHPOI=MCHPO1
  102. NC = 1
  103. SEGINI MSOUPO
  104. SEGINI, MPOVAL=MPOVA1
  105. IPCHP(1) = MSOUPO
  106. NOCOMP(1) = 'DVAP'
  107. IGEOC = IPT1
  108. IPOVAL = MPOVAL
  109. NOHARM(1) = NIFOUR
  110. C
  111. C Creation du segment IWORK
  112. C Ordre des constituants : H2O, H2, He, N2, O2, CO2, CO
  113. C
  114. SEGACT IPT1
  115. NBCELL = IPT1.NUM(/2)
  116. SEGDES IPT1
  117. NBCONS = 7
  118. SEGINI IWORK
  119. C
  120. C- Liste des constituants autorises
  121. BETA(1) = 'H2O '
  122. BETA(2) = 'H2 '
  123. BETA(3) = 'HE '
  124. BETA(4) = 'N2 '
  125. BETA(5) = 'O2 '
  126. BETA(6) = 'CO2 '
  127. BETA(7) = 'CO '
  128. C
  129. C- Masse molaire (en kg/mol)
  130. V(1,1) = 0.018D0
  131. V(2,1) = 0.002D0
  132. V(3,1) = 0.004D0
  133. V(4,1) = 0.028D0
  134. V(5,1) = 0.032D0
  135. V(6,1) = 0.044D0
  136. V(7,1) = 0.028D0
  137. C
  138. C- Coefficient intervenant dans l'expression de la loi de diffusion
  139. C- de la vapeur dans le constituant de rang i lorsqu'il est seul
  140. V(1,2) = 13.1
  141. V(2,2) = 6.12
  142. V(3,2) = 2.67
  143. V(4,2) = 18.5
  144. V(5,2) = 16.3
  145. V(6,2) = 26.9
  146. V(7,2) = 18.0
  147. C
  148. SEGACT MTAB3
  149. IDIM1 = MTAB3.MLOTAB
  150. DO 40 I=1,IDIM1
  151. CALL ACEM(MTAB3,I,MOT1)
  152. MOT2 = 'R' // MOT1(1:7)
  153. TYPE = 'CHPOINT '
  154. CALL ACMO(MTAB2,MOT2,TYPE,MCHPO3)
  155. IF (IERR.NE.0) RETURN
  156. II = 0
  157. DO 10 K=1,NBCONS
  158. IF (BETA(K).EQ.MOT1) THEN
  159. II = K
  160. GOTO 20
  161. ENDIF
  162. 10 CONTINUE
  163. C Le mot %m1:4 n'est pas un nom de composante reconnu
  164. MOTERR(1:4) = MOT1(1:4)
  165. CALL ERREUR(197)
  166. RETURN
  167. 20 CONTINUE
  168. CALL LRCHT(MCHPO3,MPOVA3,TYPS,IGEO1)
  169. IF (IGEO1.NE.IPT1) THEN
  170. INDIC = 1
  171. NBCOMP = -1
  172. CALL QUEPOI(MCHPO3,IPT1,INDIC,NBCOMP,NOMTOT)
  173. IF (IERR.NE.0) RETURN
  174. CALL LRCHT(MCHPO3,MPOVA3,TYPS,IGEO1)
  175. ENDIF
  176. A1 = V(II,1)
  177. A2 = 1.D0 / V(II,1)
  178. A3 = 1.D0 / V(1,1)
  179. WAB = 2.D3 / (A2 + A3)
  180. A4 = (V(II,2)**(1.D0/3.D0) + V(1,2)**(1.D0/3.D0))**2.D0
  181. DO 30 J=1,NBCELL
  182. GIIJ = MPOVA3.VPOCHA(J,1) / A1
  183. W(J,1) = W(J,1) + GIIJ
  184. IF (II.NE.1) THEN
  185. DVIIJ = (1.D-4 * 0.00143D0) * (MPOVA2.VPOCHA(J,1) ** 1.75)
  186. & / ((MPOVA1.VPOCHA(J,1) * 1.D-5)*(WAB ** 0.5D0)*A4)
  187. W(J,2) = W(J,2) + GIIJ/DVIIJ
  188. ELSE
  189. VPOCHA(J,1) = GIIJ
  190. ENDIF
  191. 30 CONTINUE
  192. SEGDES MPOVA3
  193. 40 CONTINUE
  194. SEGDES MTAB3,MPOVA1,MPOVA2
  195. C
  196. C- Remplissage du chpo contenant le coefficient de diffusion
  197. C- de la vapeur dans les incondensables
  198. C
  199. DO 50 I=1,NBCELL
  200. VPOCHA(I,1) = (1.D0 - VPOCHA(I,1)/W(I,1)) * W(I,1) / W(I,2)
  201. 50 CONTINUE
  202. SEGDES MCHPOI,MSOUPO,MPOVAL
  203. C
  204. SEGSUP IWORK
  205. CALL ECROBJ('CHPOINT ',MCHPOI)
  206. RETURN
  207. END
  208.  
  209.  
  210.  
  211.  
  212.  
  213.  
  214.  

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