Télécharger copain.eso

Retour à la liste

Numérotation des lignes :

  1. C COPAIN SOURCE SERRE 16/12/13 21:15:04 9254
  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. -INC CCOPTIO
  65. -INC SMCHPOI
  66. -INC SMCOORD
  67. -INC SMELEME
  68. -INC SMTABLE
  69. SEGMENT REDIR
  70. INTEGER IND1(NNGOT),IND2(NNGOT),IND3(NNGOT)
  71. ENDSEGMENT
  72. C
  73. C- Tableau d'adressage indirect
  74. C
  75. SEGACT IPT1,IPT2,IPT3,IPT4
  76. SEGACT MPOVA1,MPOVA2,MPOVA3,MPOVA4,MPOVA5,MPOVA6,MPOVA7,MPOVA8,
  77. & MPOVA9,MPOVAA,MPOVAB,MPOVAC,MPOVAD,MPOVAE,MPOVACP,MPOVAF,
  78. & MPOVAI,MPOVAJ
  79. NNGOT = XCOOR(/1)/(IDIM+1)
  80. SEGINI REDIR
  81. NBLAY = IPT2.NUM(/2)
  82. DO 10 I=1,NBLAY
  83. K = IPT2.NUM(1,I)
  84. IND2(K) = I
  85. 10 CONTINUE
  86. NBWAL = IPT3.NUM(/2)
  87. DO 20 I=1,NBWAL
  88. K = IPT3.NUM(1,I)
  89. IND3(K) = I
  90. 20 CONTINUE
  91. NBCEL = IPT1.NUM(/2)
  92. DO 30 I=1,NBCEL
  93. K = IPT1.NUM(1,I)
  94. IND1(K) = I
  95. 30 CONTINUE
  96. UST = 1.D0/3.D0
  97. DST = 2.D0/3.D0
  98. C
  99. C- Pour chaque couple CELL/LAYERW on calcule les grandeurs liees a la
  100. C condensation
  101. C
  102. SEGINI, MPOUT1=MPOVA7
  103. N = MPOVA7.VPOCHA(/1)
  104. NC = 1
  105. SEGINI MPOUT2,MPOUT3,MPOUT4
  106. N1 = IPT4.NUM(/2)
  107. DO 70 I=1,N1
  108. C
  109. C Coeff d'echange : correction pour les murs en beton
  110. C Evaluation du debit de condensation
  111. C
  112. IPOS = IND2(IPT4.NUM(2,I))
  113. IF (IPOS.NE.0) THEN
  114. II1 = IND1(IPT4.NUM(1,I))
  115. JPOS = IND3(IPT4.NUM(2,I))
  116. PTOT1 = MPOVA1.VPOCHA(II1,1)
  117. XTGAS1 = MPOVA2.VPOCHA(II1,1)
  118. RTOT1 = MPOVA3.VPOCHA(II1,1)
  119. RVAP3 = MPOVA4.VPOCHA(II1,1)
  120. PVP1 = MPOVA5.VPOCHA(II1,1)
  121. XTSAT1 = MPOVA6.VPOCHA(II1,1)
  122. TMUR1 = MPOVA8.VPOCHA(JPOS,1)
  123. AREA1 = MPOVA9.VPOCHA(IPOS,1)
  124. QW2 = MPOVAA.VPOCHA(IPOS,1)
  125. XLBD1 = MPOVAB.VPOCHA(II1,1)
  126. XNU1 = MPOVAC.VPOCHA(II1,1)
  127. RINC1 = MPOVAD.VPOCHA(II1,1)
  128. DV2 = MPOVAE.VPOCHA(II1,1)
  129. RMEL = MPOVAF.VPOCHA(II1,1)
  130. WNC = MPOVAI.VPOCHA(II1,1)
  131. WVAP = MPOVAJ.VPOCHA(II1,1)
  132. XCP = MPOVACP.VPOCHA(II1,1)
  133. C
  134. XAL = XLBD1/(RTOT1*XCP)
  135. XMU = RTOT1*XNU1
  136. PSATP1 = PSATT0(TMUR1)
  137. HL = HLS0(PTOT1,TMUR1)
  138. HV0 = HVS0(PVP1,XTGAS1)
  139. HV = HVS0(PSATP1,TMUR1)
  140. HLAT = HV - HL
  141. HLAT0 = HV0 - HL
  142. DT1 = XTSAT1 - TMUR1
  143. DT2 = TMUR1 - XTGAS1
  144. ROITP1 =(PTOT1 - PSATP1) / (RINC1 * TMUR1)
  145. XVAPI = PSATP1/PTOT1
  146. XINCI = (PTOT1-PSATP1)/PTOT1
  147. WVAP = 0.01801528D0
  148. WMEL = XVAPI*WVAP+XINCI*WNC
  149. RMEL1 = 8.31446210D0 / WMEL
  150. rtot22 = rovap0(psatp1,tmur1) + roitp1
  151. RTOT2 = PTOT1 / (RMEL1 * TMUR1)
  152. ROITP2 = RTOT1 - RVAP3
  153. XNC = ROITP2 / RTOT1
  154. XNCIT = ROITP1 / RTOT22
  155. C
  156. XUN = (XNCIT - XNC) / 0.01D0
  157. XUN = MAX(XUN,0.D0)
  158. XUN = MIN(XUN,1.D0)
  159. XDEU = (XTSAT1 - TMUR1) / 0.2D0
  160. XDEU = MAX(XDEU,0.D0)
  161. XDEU = MIN(XDEU,1.D0)
  162. FXUN = (XUN**2.D0) * (3.D0 - 2.D0*XUN)
  163. FXDEU = (XDEU**2.D0) * (3.D0 - 2.D0*XDEU)
  164. C
  165. XNCI = FXUN*XNCIT + (1.D0 - FXUN) * XNC
  166. RAPP = (1.D0 - XNC / XNCIT) * FXDEU
  167. C
  168. C Condition implicite de positivité de Hcond (plus d'incondensable en paroi)
  169. RAPP = MAX(0.D0,RAPP)
  170. THET = 0.8254D0 + 0.616D0*((PVP1-PSATP1)/(PTOT1-PSATP1))
  171. FACT1 = (-1.D0*DT2/XTGAS1)
  172. FACT2 = ((XNCI-XNC)/(((WNC)/(WNC-WVAP))-XNC))
  173. XBETA = ABS(FACT1 + FACT2)
  174. HCONV = THET*0.13D0*(((G1*XBETA)/(XNU1*XAL)) ** UST) * XLBD1
  175. HCONV = MIN(HCONV,HMAX)
  176. C
  177. IF (DT2 .LT. -0.2D0) THEN
  178. HCOND = THET*0.13D0*(((G1*XBETA)/(XNU1)) ** UST)
  179. & *((RTOT1*(DV2**DST)*RAPP*HLAT0)/((-1.D0)*DT2))
  180. HCOND = MAX(HCOND,0.D0)
  181. DEBIC1 = -1.D0*(AREA1*DT2)*HCOND/HLAT0
  182. DEBIC1 = MAX(DEBIC1,0.D0)
  183. MPOUT1.VPOCHA(IPOS,1) = HCONV + HCOND
  184. QW2 = MPOUT1.VPOCHA(IPOS,1) * DT2
  185. MPOUT2.VPOCHA(IPOS,1) = DEBIC1 / QW2
  186. ELSE
  187. MPOUT1 . VPOCHA(IPOS,1) = HCONV
  188. MPOUT2 . VPOCHA(IPOS,1) = 0.D0
  189. ENDIF
  190. HTOT1 = MPOUT1.VPOCHA(IPOS,1)
  191. IF (HTOT1 . LT . MPOVA7.VPOCHA(IPOS,1)) THEN
  192. MPOUT1 . VPOCHA(IPOS,1) = MPOVA7.VPOCHA(IPOS,1)
  193. HCONV = HCONV + MPOVA7.VPOCHA(IPOS,1) - HTOT1
  194. ENDIF
  195. MPOUT4 . VPOCHA(IPOS,1) = HCONV
  196. MPOUT3 . VPOCHA(IPOS,1) = HL
  197. ENDIF
  198. 70 CONTINUE
  199. C
  200. C- Menage
  201. C
  202. SEGSUP REDIR
  203. SEGDES IPT1,IPT2,IPT3,IPT4,MPOVA1,MPOVA2,MPOVA3,MPOVA4,MPOVA5,
  204. & MPOVA6,MPOVA7,MPOVA8,MPOVA9,MPOVAA,MPOVAB,MPOVAC,MPOVAD,
  205. & MPOVAE,MPOVAF,MPOVAI,MPOVAJ,MPOVACP
  206. C
  207. RETURN
  208. END
  209.  
  210.  

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