Télécharger copain.eso

Retour à la liste

Numérotation des lignes :

  1. C COPAIN SOURCE PV 20/03/30 21:16:38 10567
  2. SUBROUTINE COPAIN(IPT1,IPT2,IPT3,IPT4,
  3. & MPOVA1,MPOVA2,MPOVA3,MPOVA4,MPOVA5,MPOVA6,
  4. & MPOVA7,MPOVA8,MPOVA9,MPOVAA,
  5. & MPOVAB,MPOVAC,MPOVAD,MPOVAE,
  6. & G1,BETA,HMAX,MPOVACP,MPOVAF,MPOVAI,MPOVAJ,
  7. & MPOUT1,MPOUT2,MPOUT3,MPOUT4)
  8. C--------------------------------------------------------------------
  9. C Evaluation du coef. d'echange thermique en paroi et du debit de
  10. C condensation lorsque la correlation retenue est celle de COPAIN
  11. C--------------------------------------------------------------------
  12. C
  13. C E/ IPT1 : Maillage des points CELL
  14. C E/ IPT2 : Maillage des points LAYERW
  15. C E/ IPT3 : Maillage des points WALL
  16. C E/ IPT4 : Maillage de SEG2 CELL/LAYERW
  17. C E/ MPOVA1 : Pression totale
  18. C E/ MPOVA2 : Temperature du melange gazeux
  19. C E/ MPOVA3 : Densitee du melange - (en kg/m3)
  20. C E/ MPOVA4 : Densitee partielle de vapeur (en kg/m3)
  21. C E/ MPOVA5 : Pression partielle de vapeur
  22. C E/ MPOVA6 : Temperature de saturation
  23. C E/ MPOVA7 : Coefficient d'echange par conduction
  24. C E/ MPOVA8 : Temperature des murs
  25. C E/ MPOVA9 : Surface des parois
  26. C E/ MPOVAA : Flux de chaleur en paroi
  27. C E/ MPOVAB : Lambda du melange
  28. C E/ MPOVAC : Nu du melange
  29. C E/ MPOVAD : Constante des gaz parfait des incondensables
  30. C E/ MPOVAE : Coeff de diffusion de la vapeur
  31. C E/ MPOVAF : Constante des gaz parfait du melange a Tgas
  32. C E/ MPOVAI : Masse molaire des incondensables
  33. C E/ MPOVAJ : Masse molaire de la vapeur d'eau
  34. C E/ G1 : Gravite (aeleration de la pesanteur)
  35. C E/ BETA : Fraction de l'energie associee a la condensation
  36. C transferee au mur (defaut 92%)chilton
  37. C E/ HMAX : Borne sup pour le coefficient de transfert thermique
  38. C E/ MPOVACP: Coefficient calorifique de l'eau
  39. C /S MPOUT1 : Coefficient d'echange global en paroi
  40. C /S MPOUT2 : Debit de condensation / flux thermique
  41. C /S MPOUT3 : Enthalpie liquide
  42. C /S MPOUT4 : Coefficient d'échange convectif en paroi
  43. C
  44. C--------------------------------------------------------------------
  45. C Les segments MPOUT... ne sont pas desactives
  46. C--------------------------------------------------------------------
  47. C 03/07/02 : Le debit de condensation est mis a 0. si il est positif
  48. C FD On se sert de QW2 update afin d'en lisser les variations
  49. C
  50. C
  51. C--------------------------------------------------------------------
  52. C
  53. C Langage : ESOPE + FORTRAN77
  54. C Auteurs :
  55. C
  56. C-----------------------------------------------------------------------
  57. IMPLICIT INTEGER(I-N)
  58. IMPLICIT REAL*8 (A-H,O-Z)
  59. POINTEUR MPOVA7.MPOVAL,MPOVA8.MPOVAL,MPOVA9.MPOVAL,MPOVAA.MPOVAL
  60. POINTEUR MPOVAB.MPOVAL,MPOVAC.MPOVAL,MPOVAD.MPOVAL,MPOVAE.MPOVAL
  61. POINTEUR MPOVAF.MPOVAL,MPOVAI.MPOVAL,MPOVAJ.MPOVAL,MPOVACP.MPOVAL
  62. POINTEUR MPOUT1.MPOVAL,MPOUT2.MPOVAL,MPOUT3.MPOVAL,MPOUT4.MPOVAL
  63. C
  64.  
  65. -INC PPARAM
  66. -INC CCOPTIO
  67. -INC SMCHPOI
  68. -INC SMCOORD
  69. -INC SMELEME
  70. -INC SMTABLE
  71. SEGMENT REDIR
  72. INTEGER IND1(NNGOT),IND2(NNGOT),IND3(NNGOT)
  73. ENDSEGMENT
  74. C
  75. C- Tableau d'adressage indirect
  76. C
  77. SEGACT IPT1,IPT2,IPT3,IPT4
  78. SEGACT MPOVA1,MPOVA2,MPOVA3,MPOVA4,MPOVA5,MPOVA6,MPOVA7,MPOVA8,
  79. & MPOVA9,MPOVAA,MPOVAB,MPOVAC,MPOVAD,MPOVAE,MPOVACP,MPOVAF,
  80. & MPOVAI,MPOVAJ
  81. NNGOT = nbpts
  82. SEGINI REDIR
  83. NBLAY = IPT2.NUM(/2)
  84. DO 10 I=1,NBLAY
  85. K = IPT2.NUM(1,I)
  86. IND2(K) = I
  87. 10 CONTINUE
  88. NBWAL = IPT3.NUM(/2)
  89. DO 20 I=1,NBWAL
  90. K = IPT3.NUM(1,I)
  91. IND3(K) = I
  92. 20 CONTINUE
  93. NBCEL = IPT1.NUM(/2)
  94. DO 30 I=1,NBCEL
  95. K = IPT1.NUM(1,I)
  96. IND1(K) = I
  97. 30 CONTINUE
  98. UST = 1.D0/3.D0
  99. DST = 2.D0/3.D0
  100. C
  101. C- Pour chaque couple CELL/LAYERW on calcule les grandeurs liees a la
  102. C condensation
  103. C
  104. SEGINI, MPOUT1=MPOVA7
  105. N = MPOVA7.VPOCHA(/1)
  106. NC = 1
  107. SEGINI MPOUT2,MPOUT3,MPOUT4
  108. N1 = IPT4.NUM(/2)
  109. DO 70 I=1,N1
  110. C
  111. C Coeff d'echange : correction pour les murs en beton
  112. C Evaluation du debit de condensation
  113. C
  114. IPOS = IND2(IPT4.NUM(2,I))
  115. IF (IPOS.NE.0) THEN
  116. II1 = IND1(IPT4.NUM(1,I))
  117. JPOS = IND3(IPT4.NUM(2,I))
  118. PTOT1 = MPOVA1.VPOCHA(II1,1)
  119. XTGAS1 = MPOVA2.VPOCHA(II1,1)
  120. RTOT1 = MPOVA3.VPOCHA(II1,1)
  121. RVAP3 = MPOVA4.VPOCHA(II1,1)
  122. PVP1 = MPOVA5.VPOCHA(II1,1)
  123. XTSAT1 = MPOVA6.VPOCHA(II1,1)
  124. TMUR1 = MPOVA8.VPOCHA(JPOS,1)
  125. AREA1 = MPOVA9.VPOCHA(IPOS,1)
  126. QW2 = MPOVAA.VPOCHA(IPOS,1)
  127. XLBD1 = MPOVAB.VPOCHA(II1,1)
  128. XNU1 = MPOVAC.VPOCHA(II1,1)
  129. RINC1 = MPOVAD.VPOCHA(II1,1)
  130. DV2 = MPOVAE.VPOCHA(II1,1)
  131. RMEL = MPOVAF.VPOCHA(II1,1)
  132. WNC = MPOVAI.VPOCHA(II1,1)
  133. WVAP = MPOVAJ.VPOCHA(II1,1)
  134. XCP = MPOVACP.VPOCHA(II1,1)
  135. C
  136. XAL = XLBD1/(RTOT1*XCP)
  137. XMU = RTOT1*XNU1
  138. PSATP1 = PSATT0(TMUR1)
  139. HL = HLS0(PTOT1,TMUR1)
  140. HV0 = HVS0(PVP1,XTGAS1)
  141. HV = HVS0(PSATP1,TMUR1)
  142. HLAT = HV - HL
  143. HLAT0 = HV0 - HL
  144. DT1 = XTSAT1 - TMUR1
  145. DT2 = TMUR1 - XTGAS1
  146. ROITP1 =(PTOT1 - PSATP1) / (RINC1 * TMUR1)
  147. XVAPI = PSATP1/PTOT1
  148. XINCI = (PTOT1-PSATP1)/PTOT1
  149. WVAP = 0.01801528D0
  150. WMEL = XVAPI*WVAP+XINCI*WNC
  151. RMEL1 = 8.31446210D0 / WMEL
  152. rtot22 = rovap0(psatp1,tmur1) + roitp1
  153. RTOT2 = PTOT1 / (RMEL1 * TMUR1)
  154. ROITP2 = RTOT1 - RVAP3
  155. XNC = ROITP2 / RTOT1
  156. XNCIT = ROITP1 / RTOT22
  157. C
  158. XUN = (XNCIT - XNC) / 0.01D0
  159. XUN = MAX(XUN,0.D0)
  160. XUN = MIN(XUN,1.D0)
  161. XDEU = (XTSAT1 - TMUR1) / 0.2D0
  162. XDEU = MAX(XDEU,0.D0)
  163. XDEU = MIN(XDEU,1.D0)
  164. FXUN = (XUN**2.D0) * (3.D0 - 2.D0*XUN)
  165. FXDEU = (XDEU**2.D0) * (3.D0 - 2.D0*XDEU)
  166. C
  167. XNCI = FXUN*XNCIT + (1.D0 - FXUN) * XNC
  168. RAPP = (1.D0 - XNC / XNCIT) * FXDEU
  169. C
  170. C Condition implicite de positivité de Hcond (plus d'incondensable en paroi)
  171. RAPP = MAX(0.D0,RAPP)
  172. THET = 0.8254D0 + 0.616D0*((PVP1-PSATP1)/(PTOT1-PSATP1))
  173. FACT1 = (-1.D0*DT2/XTGAS1)
  174. FACT2 = ((XNCI-XNC)/(((WNC)/(WNC-WVAP))-XNC))
  175. XBETA = ABS(FACT1 + FACT2)
  176. HCONV = THET*0.13D0*(((G1*XBETA)/(XNU1*XAL)) ** UST) * XLBD1
  177. HCONV = MIN(HCONV,HMAX)
  178. C
  179. IF (DT2 .LT. -0.2D0) THEN
  180. HCOND = THET*0.13D0*(((G1*XBETA)/(XNU1)) ** UST)
  181. & *((RTOT1*(DV2**DST)*RAPP*HLAT0)/((-1.D0)*DT2))
  182. HCOND = MAX(HCOND,0.D0)
  183. DEBIC1 = -1.D0*(AREA1*DT2)*HCOND/HLAT0
  184. DEBIC1 = MAX(DEBIC1,0.D0)
  185. MPOUT1.VPOCHA(IPOS,1) = HCONV + HCOND
  186. QW2 = MPOUT1.VPOCHA(IPOS,1) * DT2
  187. MPOUT2.VPOCHA(IPOS,1) = DEBIC1 / QW2
  188. ELSE
  189. MPOUT1 . VPOCHA(IPOS,1) = HCONV
  190. MPOUT2 . VPOCHA(IPOS,1) = 0.D0
  191. ENDIF
  192. HTOT1 = MPOUT1.VPOCHA(IPOS,1)
  193. IF (HTOT1 . LT . MPOVA7.VPOCHA(IPOS,1)) THEN
  194. MPOUT1 . VPOCHA(IPOS,1) = MPOVA7.VPOCHA(IPOS,1)
  195. HCONV = HCONV + MPOVA7.VPOCHA(IPOS,1) - HTOT1
  196. ENDIF
  197. MPOUT4 . VPOCHA(IPOS,1) = HCONV
  198. MPOUT3 . VPOCHA(IPOS,1) = HL
  199. ENDIF
  200. 70 CONTINUE
  201. C
  202. C- Menage
  203. C
  204. SEGSUP REDIR
  205. SEGDES IPT1,IPT2,IPT3,IPT4,MPOVA1,MPOVA2,MPOVA3,MPOVA4,MPOVA5,
  206. & MPOVA6,MPOVA7,MPOVA8,MPOVA9,MPOVAA,MPOVAB,MPOVAC,MPOVAD,
  207. & MPOVAE,MPOVAF,MPOVAI,MPOVAJ,MPOVACP
  208. C
  209. RETURN
  210. END
  211.  
  212.  
  213.  

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