Télécharger chilto.eso

Retour à la liste

Numérotation des lignes :

  1. C CHILTO SOURCE SERRE 16/12/13 21:15:02 9254
  2. SUBROUTINE CHILTO(IPT1,IPT2,IPT3,IPT4,
  3. & MPOVA1,MPOVA2,MPOVA3,MPOVA4,MPOVA5,
  4. & MPOVA6,MPOVA7,MPOVA8,MPOVA9,MPOVAA,
  5. & MPOVAB,MPOVAC,MPOVAD,MPOVAE,G1,HMAX,
  6. & MPOUT1,MPOUT2,MPOUT3,MPOUT4)
  7. C--------------------------------------------------------------------
  8. C Evaluation du coef. d'échange thermique en paroi et du débit de
  9. C condensation lorsque la corrélation retenue est celle de CHILTON
  10. C--------------------------------------------------------------------
  11. C
  12. C E/ IPT1 : Maillage des points CELL
  13. C E/ IPT2 : Maillage des points LAYERW
  14. C E/ IPT3 : Maillage des points WALL
  15. C E/ IPT4 : Maillage de SEG2 CELL/LAYERW
  16. C E/ MPOVA1 : Pression totale
  17. C E/ MPOVA2 : Température du mélange gazeux
  18. C E/ MPOVA3 : Densité du mélange - (en kg/m3)
  19. C E/ MPOVA4 : Densité partielle de vapeur (en kg/m3)
  20. C E/ MPOVA5 : Pression partielle de vapeur
  21. C E/ MPOVA6 : Température de saturation
  22. C E/ MPOVA7 : Coefficient d'échange par conduction
  23. C E/ MPOVA8 : Température des murs
  24. C E/ MPOVA9 : Surface des parois
  25. C E/ MPOVAA : Flux de chaleur en paroi
  26. C E/ MPOVAB : Lambda du mélange
  27. C E/ MPOVAC : Nu du mélange
  28. C E/ MPOVAD : Constante des gaz parfait des incondensables
  29. C E/ MPOVAE : Coeff de diffusion de la vapeur
  30. C E/ G1 : Gravité (accélération de la pesanteur)
  31. C E/ HMAX : Borne sup pour le coefficient de transfert thermique
  32. C /S MPOUT1 : Coefficient d'échange global en paroi
  33. C /S MPOUT2 : Débit de condensation / flux thermique
  34. C /S MPOUT3 : Enthalpie liquide
  35. C /S MPOUT4 : Coefficient d'échange par convection
  36. C
  37. C--------------------------------------------------------------------
  38. C Les segments MPOUT... ne sont pas désactivés
  39. C--------------------------------------------------------------------
  40. C 03/07/02 : Le débit de condensation est mis à 0. si il est positif
  41. C FD On se sert de QW2 updaté afin d'en lisser les variations
  42. C
  43. C
  44. C--------------------------------------------------------------------
  45. C
  46. C Langage : ESOPE + FORTRAN77
  47. C Auteurs : F.DABBENE 12/01
  48. C
  49. C--------------------------------------------------------------------
  50. IMPLICIT INTEGER(I-N)
  51. IMPLICIT REAL*8 (A-H,O-Z)
  52. POINTEUR MPOVA7.MPOVAL,MPOVA8.MPOVAL,MPOVA9.MPOVAL,MPOVAA.MPOVAL
  53. POINTEUR MPOVAB.MPOVAL,MPOVAC.MPOVAL,MPOVAD.MPOVAL,MPOVAE.MPOVAL
  54. POINTEUR MPOUT1.MPOVAL,MPOUT2.MPOVAL,MPOUT3.MPOVAL,MPOUT4.MPOVAL
  55. C
  56. -INC CCOPTIO
  57. -INC SMCHPOI
  58. -INC SMCOORD
  59. -INC SMELEME
  60. -INC SMTABLE
  61. SEGMENT REDIR
  62. INTEGER IND1(NNGOT),IND2(NNGOT),IND3(NNGOT)
  63. ENDSEGMENT
  64. C
  65. C- Tableau d'adressage indirect
  66. C
  67. SEGACT IPT1,IPT2,IPT3,IPT4
  68. SEGACT MPOVA1,MPOVA2,MPOVA3,MPOVA4,MPOVA5,MPOVA6,MPOVA7,MPOVA8,
  69. & MPOVA9,MPOVAA,MPOVAB,MPOVAC,MPOVAD,MPOVAE
  70. NNGOT = XCOOR(/1)/(IDIM+1)
  71. SEGINI REDIR
  72. NBLAY = IPT2.NUM(/2)
  73. DO 10 I=1,NBLAY
  74. K = IPT2.NUM(1,I)
  75. IND2(K) = I
  76. 10 CONTINUE
  77. NBWAL = IPT3.NUM(/2)
  78. DO 20 I=1,NBWAL
  79. K = IPT3.NUM(1,I)
  80. IND3(K) = I
  81. 20 CONTINUE
  82. NBCEL = IPT1.NUM(/2)
  83. DO 30 I=1,NBCEL
  84. K = IPT1.NUM(1,I)
  85. IND1(K) = I
  86. 30 CONTINUE
  87. UST = 1.D0/3.D0
  88. DST = 2.D0/3.D0
  89. C
  90. C- Pour chaque couple CELL/LAYERW on calcule les grandeurs liees a la
  91. C condensation
  92. C
  93. SEGINI, MPOUT1=MPOVA7
  94. N = MPOVA7.VPOCHA(/1)
  95. NC = 1
  96. SEGINI MPOUT2,MPOUT3,MPOUT4
  97. N1 = IPT4.NUM(/2)
  98. DO 70 I=1,N1
  99. C
  100. C Coeff d'échange : correction pour les murs en béton
  101. C Evaluation du débit de condensation
  102. C
  103. IPOS = IND2(IPT4.NUM(2,I))
  104. IF (IPOS.NE.0) THEN
  105. II1 = IND1(IPT4.NUM(1,I))
  106. JPOS = IND3(IPT4.NUM(2,I))
  107. XLBD1 = MPOVAB.VPOCHA(II1,1)
  108. XNU1 = MPOVAC.VPOCHA(II1,1)
  109. RINC1 = MPOVAD.VPOCHA(II1,1)
  110. DV2 = MPOVAE.VPOCHA(II1,1)
  111. XTGAS1 = MPOVA2.VPOCHA(II1,1)
  112. RTOT1 = MPOVA3.VPOCHA(II1,1)
  113. RVAP3 = MPOVA4.VPOCHA(II1,1)
  114. PVP1 = MPOVA5.VPOCHA(II1,1)
  115. XTSAT1 = MPOVA6.VPOCHA(II1,1)
  116. PTOT1 = MPOVA1.VPOCHA(II1,1)
  117. TMUR1 = MPOVA8.VPOCHA(JPOS,1)
  118. AREA1 = MPOVA9.VPOCHA(IPOS,1)
  119. QW2 = MPOVAA.VPOCHA(IPOS,1)
  120. IF (TMUR1 .LT. XTSAT1) THEN
  121. PSATP1 = PSATT0(TMUR1)
  122. RVATP1 = ROVAP0(PSATP1,TMUR1)
  123. ROITP1 = (PTOT1 - PSATP1) / (RINC1 * TMUR1)
  124. XMASS1 = ABS((RVATP1+ROITP1-RTOT1)/RTOT1)
  125. HCONV = 0.13D0*XLBD1*(((G1*XMASS1)/XNU1/XNU1)**UST)
  126. XKC1 = HCONV*(DV2**DST)/(XLBD1*(XNU1**(-1.D0*UST)))
  127. DEBIC1 = XKC1*AREA1*(RVAP3-RVATP1)/(PSATP1/PTOT1-1.D0)
  128. ELSE
  129. RVATP1 = ROVAP0(PVP1,TMUR1)
  130. ROITP1 = (PTOT1 - PVP1) / (RINC1 * TMUR1)
  131. XMASS1 = ABS((RVATP1+ROITP1-RTOT1)/RTOT1)
  132. HCONV = 0.13D0 *XLBD1 * (((G1*XMASS1)/XNU1/XNU1) ** UST)
  133. DEBIC1 = 0.D0
  134. ENDIF
  135. DEBIC1 = MIN(DEBIC1,0.D0)
  136. C
  137. C---> C'est ici que l'on fixe les échanges d'énergie comp/mur :
  138. C---> Etat thermodynamique du comp = PVAP TGAS
  139. C---> Etat thermodynamique du puisard = PTOT TMUR
  140. C---> La chaleur Latente découle de la difference comp/puisard
  141. C---> pour les états thermo ci-dessus CHOISIS A PRIORI.
  142. C
  143. HL = HLS0(PTOT1,TMUR1)
  144. HV = HVS0(PVP1,XTGAS1)
  145. HLAT = HV - HL
  146. C--------------------------------------------------------------------
  147. DT1 = XTSAT1 - TMUR1
  148. DT2 = TMUR1 - XTGAS1
  149. HCONV = MIN(HCONV,HMAX)
  150. IF (DT1.GT.0.1D0) THEN
  151. IF (DT2.LT.-0.1D0) THEN
  152. MPOUT1.VPOCHA(IPOS,1) = (DEBIC1*HLAT)/(AREA1*DT2) + HCONV
  153. QW2 = MPOUT1.VPOCHA(IPOS,1) * DT2
  154. MPOUT2.VPOCHA(IPOS,1) = -1.D0 * DEBIC1 / QW2
  155. ELSE
  156. MPOUT1.VPOCHA(IPOS,1) = HCONV
  157. MPOUT2.VPOCHA(IPOS,1) = 0.D0
  158. ENDIF
  159. ELSE
  160. MPOUT1.VPOCHA(IPOS,1) = HCONV
  161. MPOUT2.VPOCHA(IPOS,1) = 0.D0
  162. ENDIF
  163. C
  164. HTOT1 = MPOUT1.VPOCHA(IPOS,1)
  165. IF (HTOT1 . LT . MPOVA7.VPOCHA(IPOS,1)) THEN
  166. MPOUT1 . VPOCHA(IPOS,1) = MPOVA7.VPOCHA(IPOS,1)
  167. HCONV = HCONV + MPOVA7.VPOCHA(IPOS,1) - HTOT1
  168. ELSE
  169. MPOUT1 . VPOCHA(IPOS,1) = HTOT1
  170. ENDIF
  171. MPOUT3 . VPOCHA(IPOS,1) = HL
  172. MPOUT4 . VPOCHA(IPOS,1) = HCONV
  173. ENDIF
  174. 70 CONTINUE
  175. C
  176. C- Ménage
  177. C
  178. SEGSUP REDIR
  179. SEGDES IPT1,IPT2,IPT3,IPT4
  180. C
  181. RETURN
  182. END
  183.  
  184.  
  185.  
  186.  
  187.  
  188.  
  189.  

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