1 : C CHILTO SOURCE CHAT 05/01/12 21:57:25 5004 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) 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 36 : C-------------------------------------------------------------------- 37 : C Les segments MPOUT... ne sont pas désactivés 38 : C-------------------------------------------------------------------- 39 : C 03/07/02 : Le débit de condensation est mis à 0. si il est positif 40 : C FD On se sert de QW2 updaté afin d'en lisser les variations 41 : C 42 : C 43 : C-------------------------------------------------------------------- 44 : C 45 : C Langage : ESOPE + FORTRAN77 46 : C Auteurs : F.DABBENE 12/01 47 : C 48 : C-------------------------------------------------------------------- 49 : IMPLICIT INTEGER(I-N) 50 : IMPLICIT REAL*8 (A-H,O-Z) 51 : POINTEUR MPOVA7.MPOVAL,MPOVA8.MPOVAL,MPOVA9.MPOVAL,MPOVAA.MPOVAL 52 : POINTEUR MPOVAB.MPOVAL,MPOVAC.MPOVAL,MPOVAD.MPOVAL,MPOVAE.MPOVAL 53 : POINTEUR MPOUT1.MPOVAL,MPOUT2.MPOVAL,MPOUT3.MPOVAL 54 : C 55 : -INC CCOPTIO 56 : -INC SMCHPOI 57 : -INC SMCOORD 58 : -INC SMELEME 59 : -INC SMTABLE 60 : SEGMENT REDIR 61 : INTEGER IND1(NNGOT),IND2(NNGOT),IND3(NNGOT) 62 : ENDSEGMENT 63 : C 64 : C- Tableau d'adressage indirect 65 : C 66 : SEGACT IPT1,IPT2,IPT3,IPT4 67 : SEGACT MPOVA1,MPOVA2,MPOVA3,MPOVA4,MPOVA5,MPOVA6,MPOVA7,MPOVA8, 68 : & MPOVA9,MPOVAA,MPOVAB,MPOVAC,MPOVAD,MPOVAE 69 : NNGOT = XCOOR(/1)/(IDIM+1) 70 : SEGINI REDIR 71 : NBLAY = IPT2.NUM(/2) 72 : DO 10 I=1,NBLAY 73 : K = IPT2.NUM(1,I) 74 : IND2(K) = I 75 : 10 CONTINUE 76 : NBWAL = IPT3.NUM(/2) 77 : DO 20 I=1,NBWAL 78 : K = IPT3.NUM(1,I) 79 : IND3(K) = I 80 : 20 CONTINUE 81 : NBCEL = IPT1.NUM(/2) 82 : DO 30 I=1,NBCEL 83 : K = IPT1.NUM(1,I) 84 : IND1(K) = I 85 : 30 CONTINUE 86 : UST = 1.D0/3.D0 87 : DST = 2.D0/3.D0 88 : C 89 : C- Pour chaque couple CELL/LAYERW on calcule les grandeurs liees a la 90 : C condensation 91 : C 92 : SEGINI, MPOUT1=MPOVA7 93 : N = MPOVA7.VPOCHA(/1) 94 : NC = 1 95 : SEGINI MPOUT2 96 : SEGINI MPOUT3 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 : XLBD1 = MPOVAB.VPOCHA(II1,1) 107 : XNU1 = MPOVAC.VPOCHA(II1,1) 108 : RINC1 = MPOVAD.VPOCHA(II1,1) 109 : DV2 = MPOVAE.VPOCHA(II1,1) 110 : XTGAS1 = MPOVA2.VPOCHA(II1,1) 111 : RTOT1 = MPOVA3.VPOCHA(II1,1) 112 : RVAP3 = MPOVA4.VPOCHA(II1,1) 113 : PVP1 = MPOVA5.VPOCHA(II1,1) 114 : XTSAT1 = MPOVA6.VPOCHA(II1,1) 115 : PTOT1 = MPOVA1.VPOCHA(II1,1) 116 : C IPOS = IND2(IPT4.NUM(2,I)) 117 : JPOS = IND3(IPT4.NUM(2,I)) 118 : TMUR1 = MPOVA8.VPOCHA(JPOS,1) 119 : AREA1 = MPOVA9.VPOCHA(IPOS,1) 120 : QW2 = MPOVAA.VPOCHA(IPOS,1) 121 : IF (TMUR1.LT.XTSAT1) THEN 122 : PSATP1 = PSATT0(TMUR1) 123 : RVATP1 = ROVAP0(PSATP1,TMUR1) 124 : ROITP1 = (PTOT1 - PSATP1) / (RINC1 * TMUR1) 125 : XMASS1 = ABS((RVATP1+ROITP1-RTOT1)/RTOT1) 126 : HCONV = 0.13D0*XLBD1*(((G1*XMASS1)/XNU1/XNU1)**UST) 127 : XKC1 = HCONV*(DV2**DST)/(XLBD1*(XNU1**(-1.D0*UST))) 128 : DEBIC1 = XKC1*AREA1*(RVAP3-RVATP1)/(PSATP1/PTOT1-1.D0) 129 : C********* DEBIC1 = -1.D0*XKC1*AREA1*(RVAP3-RVATP1) 130 : ELSE 131 : RVATP1 = ROVAP0(PVP1,TMUR1) 132 : ROITP1 = (PTOT1 - PVP1) / (RINC1 * TMUR1) 133 : XMASS1 = ABS((RVATP1+ROITP1-RTOT1)/RTOT1) 134 : HCONV = 0.13D0*XLBD1*(((G1*XMASS1)/XNU1/XNU1)**UST) 135 : DEBIC1 = 0.D0 136 : ENDIF 137 : DEBIC1 = MIN(DEBIC1,0.D0) 138 : C 139 : C---> C'est ici que l'on fixe les échanges d'énergie comp/mur : 140 : C---> Etat thermodynamique du comp = PVAP TGAS 141 : C---> Etat thermodynamique du puisard = PTOT TMUR 142 : C---> La chaleur Latente découle de la difference comp/puisard 143 : C---> pour les états thermo ci-dessus CHOISIS A PRIORI. 144 : C 145 : HL = HLS0(PTOT1,TMUR1) 146 : HV = HVS0(PVP1,XTGAS1) 147 : HLAT = HV - HL 148 : C-------------------------------------------------------------------- 149 : MPOUT3.VPOCHA(IPOS,1) = HL 150 : DT1 = XTSAT1 - TMUR1 151 : DT2 = TMUR1 - XTGAS1 152 : HCONV = MIN(HCONV,HMAX) 153 : CC HCONV = MAX(HCONV,MPOVA7.VPOCHA(IPOS,1)) 154 : IF (DT1.GT.0.1D0) THEN 155 : IF (DT2.LT.-0.1D0) THEN 156 : MPOUT1.VPOCHA(IPOS,1) = (DEBIC1*HLAT)/(AREA1*DT2) 157 : & + HCONV 158 : QW2 = MPOUT1.VPOCHA(IPOS,1) * DT2 159 : MPOUT2.VPOCHA(IPOS,1) = -1.D0 * DEBIC1 / QW2 160 : ELSE 161 : MPOUT1.VPOCHA(IPOS,1) = HCONV 162 : MPOUT2.VPOCHA(IPOS,1) = 0.D0 163 : ENDIF 164 : ELSE 165 : MPOUT1.VPOCHA(IPOS,1) = HCONV 166 : MPOUT2.VPOCHA(IPOS,1) = 0.D0 167 : ENDIF 168 : MPOUT1.VPOCHA(IPOS,1) = 169 : & MAX(MPOUT1.VPOCHA(IPOS,1),MPOVA7.VPOCHA(IPOS,1)) 170 : CC WRITE(6,*) IPOS,MPOUT1.VPOCHA(IPOS,1) 171 : CC WRITE(6,*) IPOS,MPOUT2.VPOCHA(IPOS,1) 172 : CC WRITE(6,*) IPOS,MPOUT3.VPOCHA(IPOS,1) 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 :
© Cast3M 2003 - All rights reserved.
Disclaimer