1 : C TOCOND SOURCE CHAT 05/01/13 03:42:42 5004 2 : SUBROUTINE TOCOND 3 : C 4 : C-------------------------------------------------------------------- 5 : C Gestion de la condensation en paroi : évaluation du coef. d'échange 6 : C thermique en paroi, du débit de condensation, du titre vapeur 7 : C-------------------------------------------------------------------- 8 : C 9 : C--------------------------- 10 : C Phrase d'appel (GIBIANE) : 11 : C--------------------------- 12 : C 13 : C TSAT1 XV3 ALAT0 ALAT3 ECHAN1 = VARI TOCOND TAB1 MAIL1 NOMMODE; 14 : C ou 15 : C TSAT1 XV3 = VARI TOCOND TAB1 ; 16 : C 17 : C------------------------ 18 : C Opérandes et résultat : 19 : C------------------------ 20 : C 21 : C TAB1 : TABLE SET contenant l'ensemble des data 0D 22 : C MAIL1 : Maillage contenant l'ensemble des points faces associés 23 : C au modele 24 : C NOMMODE: Mot définissant le modele de condensation utilisé 25 : C ALAT0 : Enthalpie liquide 26 : C ALAT3 : Débit de condensation / Flux thermique 27 : C ECHAN1 : Coefficient d'échange global en paroi 28 : C TSAT1 : Température de saturation 29 : C XV3 : Titre vapeur 30 : C 31 : C-------------------------------------------------------------------- 32 : C 33 : IMPLICIT INTEGER(I-N) 34 : IMPLICIT REAL*8 (A-H,O-Z) 35 : CHARACTER*4 NOMTOT(1) 36 : CHARACTER*8 TYPE,TYPS,MTYPI,MTYPR,MODCON(3) 37 : CHARACTER*72 CHARR,MOT1 38 : LOGICAL LOGI,LOGR 39 : DATA MODCON/'TAGAMI ','UCHIDA ','CHILTON '/ 40 : C 41 : -INC CCOPTIO 42 : -INC SMCHPOI 43 : -INC SMELEME 44 : -INC SMTABLE 45 : POINTEUR MPOVA7.MPOVAL,MPOVA8.MPOVAL,MPOVA9.MPOVAL,MPOVAA.MPOVAL, 46 : & MPOVAB.MPOVAL,MPOVAC.MPOVAL,MPOVAD.MPOVAL,MPOVAE.MPOVAL, 47 : & MPOVAF.MPOVAL, 48 : & MPOMIA.MPOVAL,MPOMIB.MPOVAL,MPOMIC.MPOVAL,MPOMID.MPOVAL, 49 : & MPOMIE.MPOVAL,MPODV1.MPOVAL, 50 : & MPOUT1.MPOVAL,MPOUT2.MPOVAL,MPOUT3.MPOVAL 51 : C 52 : TYPS = ' ' 53 : TMAX = 500.D0 54 : HMAX = 5000.D0 55 : BETA = 0.92D0 56 : GAMA = 0.5D0 57 : TTAG = 0.D0 58 : TAGA = 0.D0 59 : TIME2 = 0.D0 60 : C 61 : C- Lecture et controles des données d'entrée, 62 : C 63 : C MTAB1 : Table SET 64 : C Lecture de la table SET 64 : CALL LIRTAB('SET',MTAB1,1,IRETOU) 65 : IF (IERR.NE.0) RETURN 66 : C Lecture du maillage de points envoyé 66 : CALL LIROBJ('MAILLAGE',IPT2,0,IRETOU) 67 : IF (IERR.NE.0) RETURN 68 : C Lecture du modele de condensation utilise 69 : ICOR = 0 69 : CALL LIRMOT(MODCON,3,ICOR,0) 70 : IF (IERR.NE.0) RETURN 71 : 72 : C MTAB2 : Table INCO 73 : C MTAB3 : Table COMPONENT 74 : C MTAB4 : Table GEOINF 75 : C MTAB5 : Table DATA 76 : C MTAB6 : Table PHYSICAL 77 : C MTAB7 : Table CELLAY 78 : C 79 : MTYPI = 'MOT ' 80 : MTYPR = 'TABLE ' 80 : CALL ACCTAB(MTAB1,MTYPI,IVALI,XVALI,'INCO',LOGI,IRETI, 81 : & MTYPR,IVALR,XVALR,CHARR,LOGR,MTAB2) 82 : IF (IERR.NE.0) RETURN 82 : CALL ACCTAB(MTAB1,MTYPI,IVALI,XVALI,'COMPONENT',LOGI,IRETI, 83 : & MTYPR,IVALR,XVALR,CHARR,LOGR,MTAB3) 84 : IF (IERR.NE.0) RETURN 84 : CALL ACCTAB(MTAB1,MTYPI,IVALI,XVALI,'GEOINF',LOGI,IRETI, 85 : & MTYPR,IVALR,XVALR,CHARR,LOGR,MTAB4) 86 : IF (IERR.NE.0) RETURN 86 : CALL ACCTAB(MTAB1,MTYPI,IVALI,XVALI,'DATA',LOGI,IRETI, 87 : & MTYPR,IVALR,XVALR,CHARR,LOGR,MTAB5) 88 : IF (IERR.NE.0) RETURN 88 : CALL ACCTAB(MTAB1,MTYPI,IVALI,XVALI,'PHYSICAL',LOGI,IRETI, 89 : & MTYPR,IVALR,XVALR,CHARR,LOGR,MTAB6) 90 : IF (IERR.NE.0) RETURN 91 : 92 : C --------------------------------------------------------------------C 93 : C Données utilisées quelque soit le cas traité 94 : C condensation ou non ! 95 : 96 : C 97 : C IPT1 : Maillage POI1 des CELL 98 : TYPE = 'MAILLAGE' 98 : CALL ACMO(MTAB4,'CELL',TYPE,IPT1) 99 : IF (IERR.NE.0) RETURN 100 : C 101 : C MPOVA1 : Pression 102 : C MPOVA2 : Température 103 : C MPOVA3 : Densité 104 : C MPOVA4 : Densité d'H2O 105 : C MPOVA5 : Titre vapeur 106 : C 107 : TYPE = 'CHPOINT ' 107 : CALL ACMO(MTAB2,'PT',TYPE,MCHPO1) 108 : IF (IERR.NE.0) RETURN 108 : CALL LRCHT(MCHPO1,MPOVA1,TYPS,IGEO1) 109 : IF (IGEO1.NE.IPT1) THEN 110 : INDIC = 1 111 : NBCOMP = -1 111 : CALL QUEPOI(MCHPO1,IPT1,INDIC,NBCOMP,NOMTOT) 112 : IF (IERR.NE.0) RETURN 112 : CALL LRCHT(MCHPO1,MPOVA1,TYPS,IGEO1) 113 : ENDIF 114 : C 115 : TYPE = 'CHPOINT ' 115 : CALL ACMO(MTAB2,'TGAS',TYPE,MCHPO2) 116 : IF (IERR.NE.0) RETURN 116 : CALL LRCHT(MCHPO2,MPOVA2,TYPS,IGEO1) 117 : IF (IGEO1.NE.IPT1) THEN 118 : INDIC = 1 119 : NBCOMP = -1 119 : CALL QUEPOI(MCHPO2,IPT1,INDIC,NBCOMP,NOMTOT) 120 : IF (IERR.NE.0) RETURN 120 : CALL LRCHT(MCHPO2,MPOVA2,TYPS,IGEO1) 121 : ENDIF 122 : C 123 : TYPE = 'CHPOINT ' 123 : CALL ACMO(MTAB2,'RGAS',TYPE,MCHPO3) 124 : IF (IERR.NE.0) RETURN 124 : CALL LRCHT(MCHPO3,MPOVA3,TYPS,IGEO1) 125 : IF (IGEO1.NE.IPT1) THEN 126 : INDIC = 1 127 : NBCOMP = -1 127 : CALL QUEPOI(MCHPO3,IPT1,INDIC,NBCOMP,NOMTOT) 128 : IF (IERR.NE.0) RETURN 128 : CALL LRCHT(MCHPO3,MPOVA3,TYPS,IGEO1) 129 : ENDIF 130 : C 131 : TYPE = 'CHPOINT ' 131 : CALL ACMO(MTAB2,'RH2O',TYPE,MCHPO4) 132 : IF (IERR.NE.0) RETURN 132 : CALL LRCHT(MCHPO4,MPOVA4,TYPS,IGEO1) 133 : IF (IGEO1.NE.IPT1) THEN 134 : INDIC = 1 135 : NBCOMP = -1 135 : CALL QUEPOI(MCHPO4,IPT1,INDIC,NBCOMP,NOMTOT) 136 : IF (IERR.NE.0) RETURN 136 : CALL LRCHT(MCHPO4,MPOVA4,TYPS,IGEO1) 137 : ENDIF 138 : C 139 : TYPE = 'CHPOINT ' 139 : CALL ACMO(MTAB2,'XV',TYPE,MCHPO5) 140 : IF (IERR.NE.0) RETURN 140 : CALL LRCHT(MCHPO5,MPOVA5,TYPS,IGEO1) 141 : IF (IGEO1.NE.IPT1) THEN 142 : INDIC = 1 143 : NBCOMP = -1 143 : CALL QUEPOI(MCHPO5,IPT1,INDIC,NBCOMP,NOMTOT) 144 : IF (IERR.NE.0) RETURN 144 : CALL LRCHT(MCHPO5,MPOVA5,TYPS,IGEO1) 145 : ENDIF 146 : C 147 : C RELINT : Facteur de relaxation pour les paramètres physiques internes 148 : C MOT1 : Type de correlation 149 : C TIME2 : Temps physique en fin de pas de temps 150 : C TTAG : Temps physique caractéristique pour TAGAMI-UCHIDA 151 : C TAGA : Coefficient constant intervenant dans TAGAMI-UCHIDA 152 : C 153 : MTYPI = 'MOT ' 154 : MTYPR = 'FLOTTANT' 154 : CALL ACCTAB(MTAB1,MTYPI,IVALI,XVALI,'RELINT',LOGI,IRETI, 155 : & MTYPR,IVALR,RELINT,CHARR,LOGR,IRETR) 156 : IF (IERR.NE.0) RETURN 157 : 158 : C 159 : C- Répartition vapeur/eau liquide au sein de chaque compartiment et 160 : C- données thermodynamiques associées à la vapeur 161 : C 162 : C MPOVA6 : Titre vapeur actualisé 163 : C MPOVA7 : Densité vapeur 164 : C MPOVA8 : Pression partielle de vapeur 165 : C MPOVA9 : Température de saturation 166 : C 167 : C Le calcul de TSAT et XV se fait par défaut 167 : CALL TOCON1(TMAX,RELINT,MPOVA2,MPOVA4,MPOVA5, 168 : & MPOVA6,MPOVA7,MPOVA8,MPOVA9) 169 : C --------------------------------------------------------------------C 170 : 171 : 172 : C Si on veut calculer les grandeurs liées à la condensation 173 : C passons par ici ! 174 : C --------------------------------------------------------------------C 175 : C DEBUT CALCUL DE LA CONDENSATION 176 : C --------------------------------------------------------------------C 177 : IF (IPT2.NE.0.AND.ICOR.NE.0) THEN 178 : 179 : C IPT3 : Maillage POI1 des WALL 180 : C IPT4 : Maillage SEG2 de CELLAYW (connectivites CELL-LAYERW) 181 : C 182 : TYPE = 'MAILLAGE' 182 : CALL ACMO(MTAB4,'WALL',TYPE,IPT3) 183 : IF (IERR.NE.0) RETURN 183 : CALL ACMO(MTAB4,'CELLAYW',TYPE,IPT4) 184 : IF (IERR.NE.0) RETURN 185 : 186 : C si on utilise le modele de TAGAMI 187 : IF (ICOR.EQ.1) THEN 188 : MTYPI = 'MOT ' 189 : MTYPR = 'FLOTTANT' 189 : CALL ACCTAB(MTAB2,MTYPI,IVALI,XVALI,'TIME2',LOGI,IRETI, 190 : & MTYPR,IVALR,TIME2,CHARR,LOGR,IRETR) 191 : IF (IERR.NE.0) RETURN 191 : CALL ACCTAB(MTAB6,MTYPI,IVALI,XVALI,'TTAG',LOGI,IRETI, 192 : & MTYPR,IVALR,TTAG,CHARR,LOGR,IRETR) 193 : IF (IERR.NE.0) RETURN 193 : CALL ACCTAB(MTAB6,MTYPI,IVALI,XVALI,'TAGA',LOGI,IRETI, 194 : & MTYPR,IVALR,TAGA,CHARR,LOGR,IRETR) 195 : IF (IERR.NE.0) RETURN 196 : ENDIF 197 : 198 : C 199 : C MPOVAA : Coeff d'echange par conduction 200 : C MPOVAB : Température des murs 201 : C MPOVAC : Aire de la surface de contact CELL/WALL 202 : C 203 : TYPE = 'CHPOINT ' 203 : CALL ACMO(MTAB5,'ECHAC',TYPE,MCHPOA1) 204 : IF (IERR.NE.0) RETURN 204 : CALL REDUIR(MCHPOA1,IPT2,MCHPOA) 205 : IF (IERR.NE.0) RETURN 205 : CALL LRCHT(MCHPOA,MPOVAA,TYPS,IGEO1) 206 : IF (IERR.NE.0) RETURN 207 : C 208 : TYPE = 'CHPOINT ' 208 : CALL ACMO(MTAB2,'TWAL',TYPE,MCHPOB) 209 : IF (IERR.NE.0) RETURN 209 : CALL LRCHT(MCHPOB,MPOVAB,TYPS,IGEO1) 210 : IF (IGEO1.NE.IPT3) THEN 211 : INDIC = 1 212 : NBCOMP = -1 212 : CALL QUEPOI(MCHPOB,IPT3,INDIC,NBCOMP,NOMTOT) 213 : IF (IERR.NE.0) RETURN 213 : CALL LRCHT(MCHPOB,MPOVAB,TYPS,IGEO1) 214 : ENDIF 215 : C 216 : TYPE = 'CHPOINT ' 216 : CALL ACMO(MTAB5,'AREAC',TYPE,MCHPOC1) 217 : IF (IERR.NE.0) RETURN 217 : CALL REDUIR(MCHPOC1,IPT2,MCHPOC) 218 : IF (IERR.NE.0) RETURN 218 : CALL LRCHT(MCHPOC,MPOVAC,TYPS,IGEO1) 219 : IF (IERR.NE.0) RETURN 220 : C 221 : C- Evaluation du coefficient d'échange en paroi 222 : C- et du débit de condensation 223 : C 224 : IF (ICOR.EQ.1 .OR. ICOR.EQ.2) THEN 225 : C 226 : C MPOVAD : Facteur correctif du coef d'échange pour relaxaion béton 227 : C 228 : TYPE = 'CHPOINT ' 228 : CALL ACMO(MTAB6,'UCHI',TYPE,MCHPOD1) 229 : IF (IERR.NE.0) RETURN 229 : CALL REDUIR(MCHPOD1,IPT2,MCHPOD) 230 : IF (IERR.NE.0) RETURN 230 : CALL LRCHT(MCHPOD,MPOVAD,TYPS,IGEO1) 231 : IF (IERR.NE.0) RETURN 232 : C 233 : C MPOUT1 : Coefficient d'échange global en paroi 234 : C MPOUT2 : Débit de condensation / flux thermique 235 : C MPOUT3 : Enthalpie liquide 236 : C 236 : CALL UCHIDA(ICOR,IPT1,IPT2,IPT3,IPT4, 237 : & MPOVA7,MPOVA3,MPOVA4,MPOVAA,MPOVA1, 238 : & MPOVA8,MPOVA9,MPOVA2,MPOVAB,MPOVAC, 239 : & MPOVAD,HMAX,BETA,GAMA,TTAG,TIME2,TAGA, 240 : & MPOUT1,MPOUT2,MPOUT3) 241 : C 242 : ELSE 243 : C 244 : C MPOVAD : Ancien flux thermique global 245 : C 246 : TYPE = 'CHPOINT ' 246 : CALL ACMO(MTAB2,'QW',TYPE,MCHPOD1) 247 : IF (IERR.NE.0) RETURN 247 : CALL REDUIR(MCHPOD1,IPT2,MCHPOD) 248 : IF (IERR.NE.0) RETURN 248 : CALL LRCHT(MCHPOD,MPOVAD,TYPS,IGEO1) 249 : IF (IERR.NE.0) RETURN 250 : C 251 : MTYPI = 'MOT ' 252 : MTYPR = 'FLOTTANT' 252 : CALL ACCTAB(MTAB6,MTYPI,IVALI,XVALI,'GRAVITY',LOGI,IRETI, 253 : & MTYPR,IVALR,G1,CHARR,LOGR,IRETR) 254 : IF (IERR.NE.0) RETURN 255 : C 256 : C MPODV1 : Coeff de diffusion de la vapeur dans le mélange 257 : C 257 : CALL ECROBJ('TABLE ',MTAB1) 257 : CALL TODIFF 257 : CALL LIROBJ('CHPOINT ',MCHDV1,1,IRET) 258 : IF (IERR.NE.0) RETURN 258 : CALL LRCHT(MCHDV1,MPODV1,TYPS,IGEO1) 259 : IF (IGEO1.NE.IPT1) THEN 260 : INDIC = 1 261 : NBCOMP = -1 261 : CALL QUEPOI(MCHDV1,IPT1,INDIC,NBCOMP,NOMTOT) 262 : IF (IERR.NE.0) RETURN 262 : CALL LRCHT(MCHDV1,MPODV1,TYPS,IGEO1) 263 : ENDIF 264 : C 265 : C MPOMIA : Lambda du mélange 266 : C MPOMIB : Mu du mélange 267 : C MPOMIC : Nu du mélange 268 : C MPOMID : Constante des gaz parfait des incondensables 269 : C MPOMIE : Densité des incondensables 270 : C 270 : CALL TOMIXT(MTAB2,MTAB3,MTAB6,IPT1,MPOVA3, 271 : & MPOMIA,MPOMIB,MPOMIC,MPOMID,MPOMIE) 272 : C 272 : CALL CHILTO(IPT1,IPT2,IPT3,IPT4, 273 : & MPOVA1,MPOVA2,MPOVA3,MPOVA7,MPOVA8,MPOVA9, 274 : & MPOVAA,MPOVAB,MPOVAC,MPOVAD, 275 : & MPOMIA,MPOMIC,MPOMID,MPODV1, 276 : & G1,HMAX, 277 : & MPOUT1,MPOUT2,MPOUT3) 278 : C 279 : ENDIF 280 : 281 : C 282 : C- Relaxation du coeff d'échange global et du débit de condensation 283 : C 284 : C MPOVAE : Ancien Coefficient d'échange global en paroi 285 : C MPOVAF : Ancien Débit de condensation / flux thermique 286 : C 287 : TYPE = 'CHPOINT ' 287 : CALL ACMO(MTAB2,'ECHA',TYPE,MCHPOE1) 288 : IF (IERR.NE.0) RETURN 288 : CALL REDUIR(MCHPOE1,IPT2,MCHPOE) 289 : IF (IERR.NE.0) RETURN 289 : CALL LRCHT(MCHPOE,MPOVAE,TYPS,IGEO1) 290 : IF (IERR.NE.0) RETURN 291 : C 292 : TYPE = 'CHPOINT ' 292 : CALL ACMO(MTAB2,'ALAT',TYPE,MCHPOF1) 293 : IF (IERR.NE.0) RETURN 293 : CALL REDUIR(MCHPOF1,IPT2,MCHPOF) 294 : IF (IERR.NE.0) RETURN 294 : CALL LRCHT(MCHPOF,MPOVAF,TYPS,IGEO1) 295 : IF (IERR.NE.0) RETURN 296 : C 297 : SEGACT IPT2 298 : NBLAY = IPT2.NUM(/2) 299 : SEGDES IPT2 300 : DO 10 I=1,NBLAY 301 : MPOUT1.VPOCHA(I,1) = RELINT * MPOUT1.VPOCHA(I,1) 302 : & + (1.D0 - RELINT) * MPOVAE.VPOCHA(I,1) 303 : MPOUT2.VPOCHA(I,1) = RELINT * MPOUT2.VPOCHA(I,1) 304 : & + (1.D0 - RELINT) * MPOVAF.VPOCHA(I,1) 305 : 10 CONTINUE 306 : C 307 : C- Ibid chpo contenant le coefficient d'échange global 308 : C 309 : SEGINI, MCHPOI=MCHPO1 310 : NC = 1 311 : SEGINI MSOUPO 312 : IPCHP(1) = MSOUPO 313 : NOCOMP(1) = 'SCAL' 314 : IGEOC = IPT2 315 : IPOVAL = MPOUT1 316 : NOHARM(1) = NIFOUR 317 : SEGDES MCHPOI,MSOUPO,MPOUT1 317 : CALL ECROBJ('CHPOINT',MCHPOI) 318 : C 319 : C- Ibid chpo contenant le débit de condensation 320 : C 321 : SEGINI, MCHPOI=MCHPO1 322 : NC = 1 323 : SEGINI MSOUPO 324 : IPCHP(1) = MSOUPO 325 : NOCOMP(1) = 'SCAL' 326 : IGEOC = IPT2 327 : IPOVAL = MPOUT2 328 : NOHARM(1) = NIFOUR 329 : SEGDES MCHPOI,MSOUPO,MPOUT2 329 : CALL ECROBJ('CHPOINT',MCHPOI) 330 : C 331 : C- Ibid chpo contenant l'enthalpie liquide 332 : C 333 : SEGINI, MCHPOI=MCHPO1 334 : NC = 1 335 : SEGINI MSOUPO 336 : IPCHP(1) = MSOUPO 337 : NOCOMP(1) = 'SCAL' 338 : IGEOC = IPT2 339 : IPOVAL = MPOUT3 340 : NOHARM(1) = NIFOUR 341 : SEGDES MCHPOI,MSOUPO,MPOUT3 341 : CALL ECROBJ('CHPOINT',MCHPOI) 342 : C 343 : C Menage des chpo utilisé pour le calcul de la condensation 344 : C 345 : SEGDES MPOVAA,MPOVAB,MPOVAC,MPOVAD,MPOVAE 346 : 347 : ENDIF 348 : C --------------------------------------------------------------------C 349 : C FIN CALCUL DE LA CONDENSATION 350 : C --------------------------------------------------------------------C 351 : 352 : C 353 : C- Création/remplissage du chpo contenant le titre vapeur 354 : C 355 : SEGINI, MCHPOI=MCHPO1 356 : NC = 1 357 : SEGINI MSOUPO 358 : IPCHP(1) = MSOUPO 359 : NOCOMP(1) = 'SCAL' 360 : IGEOC = IPT1 361 : IPOVAL = MPOVA6 362 : NOHARM(1) = NIFOUR 363 : SEGDES MCHPOI,MSOUPO,MPOVA6 363 : CALL ECROBJ('CHPOINT',MCHPOI) 364 : C 365 : C- Ibid chpo contenant Tsat 366 : C 367 : SEGINI, MCHPOI=MCHPO1 368 : NC = 1 369 : SEGINI MSOUPO 370 : IPCHP(1) = MSOUPO 371 : NOCOMP(1) = 'SCAL' 372 : IGEOC = IPT1 373 : IPOVAL = MPOVA9 374 : NOHARM(1) = NIFOUR 375 : SEGDES MCHPOI,MSOUPO,MPOVA9 375 : CALL ECROBJ('CHPOINT',MCHPOI) 376 : C 377 : C- Ménage 378 : C 379 : SEGDES MPOVA1,MPOVA2,MPOVA3,MPOVA4,MPOVA5 380 : SEGSUP MPOVA7,MPOVA8 381 : C 382 : RETURN 383 : END 384 : 385 : 386 : 387 :
© Cast3M 2003 - All rights reserved.
Disclaimer