1 : C UCHIDA SOURCE CHAT 05/01/13 03:57:54 5004 2 : SUBROUTINE UCHIDA(ICOR,IPT1,IPT2,IPT3,IPT4, 3 : & MPOVA1,MPOVA2,MPOVA3,MPOVA4,MPOVA5, 4 : & MPOVA6,MPOVA7,MPOVA8,MPOVA9,MPOVAA, 5 : & MPOVAB,HMAX,BETA,GAMA,TTAG,TIME2,TAGA, 6 : & MPOVAC,MPOVAD,MPOVAE) 7 : C-------------------------------------------------------------------- 8 : C Evaluation du coef. d'échange thermique en paroi et du débit de 9 : C condensation lorsque la corrélation est UCHIDA ou TAGAMI-UCHIDA 10 : C-------------------------------------------------------------------- 11 : C 12 : C E/ ICOR : Corrélation à considérer (1=Tagami , 2=Uchida) 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 : Densité partielle de vapeur (en kg/m3) 18 : C E/ MPOVA2 : Densité du mélange - (en kg/m3) 19 : C E/ MPOVA3 : Densité H2O - vapeur et brouillard - (en kg/m3) 20 : C E/ MPOVA4 : Coefficient d'échange par conduction 21 : C E/ MPOVA5 : Pression totale 22 : C E/ MPOVA6 : Pression partielle de vapeur 23 : C E/ MPOVA7 : Température de saturation 24 : C E/ MPOVA8 : Température du mélange gazeux 25 : C E/ MPOVA9 : Température des murs 26 : C E/ MPOVAA : Surface des parois 27 : C E/ MPOVAB : Facteur correctif pour les murs en béton 28 : C E/ HMAX : Borne sup pour le coefficient de transfert thermique 29 : C E/ BETA : Fraction de l'énergie associée à la condensation 30 : C transférée au mur (défaut 92%) 31 : C E/ GAMA : Coef. de lissage au voisinage de Tsat (defaut 0.5) 32 : C E/ TTAG : Temps caractéristique pour Tagami 33 : C E/ TIME2 : Temps physique en fin de pas de temps 34 : C E/ TAGA : Partie constante du modèle de Tagami 35 : C /S MPOVAC : Coefficient d'échange global en paroi 36 : C /S MPOVAD : Débit de condensation / flux thermique 37 : C /S MPOVAE : Enthalpie liquide 38 : C 39 : C-------------------------------------------------------------------- 40 : C 41 : C Langage : ESOPE + FORTRAN77 42 : C Auteurs : F.DABBENE 02/01 43 : C 44 : C-------------------------------------------------------------------- 45 : IMPLICIT INTEGER(I-N) 46 : IMPLICIT REAL*8 (A-H,O-Z) 47 : POINTEUR MPOVA7.MPOVAL,MPOVA8.MPOVAL,MPOVA9.MPOVAL,MPOVAA.MPOVAL 48 : POINTEUR MPOVAB.MPOVAL,MPOVAC.MPOVAL,MPOVAD.MPOVAL,MPOVAE.MPOVAL 49 : C 50 : -INC CCOPTIO 51 : -INC SMCHPOI 52 : -INC SMCOORD 53 : -INC SMELEME 54 : -INC SMTABLE 55 : SEGMENT REDIR 56 : INTEGER IND1(NNGOT),IND2(NNGOT),IND3(NNGOT) 57 : ENDSEGMENT 58 : C 59 : C- Tableau d'adressage indirect 60 : C 61 : SEGACT IPT1,IPT2,IPT3 62 : NNGOT = XCOOR(/1)/(IDIM+1) 63 : SEGINI REDIR 64 : NBLAY = IPT2.NUM(/2) 65 : DO 10 I=1,NBLAY 66 : K = IPT2.NUM(1,I) 67 : IND2(K) = I 68 : 10 CONTINUE 69 : NBWAL = IPT3.NUM(/2) 70 : DO 20 I=1,NBWAL 71 : K = IPT3.NUM(1,I) 72 : IND3(K) = I 73 : 20 CONTINUE 74 : NBCEL = IPT1.NUM(/2) 75 : DO 30 I=1,NBCEL 76 : K = IPT1.NUM(1,I) 77 : IND1(K) = I 78 : 30 CONTINUE 79 : 80 : SEGINI, MPOVAC=MPOVA4 81 : N = MPOVA4.VPOCHA(/1) 82 : NC = 1 83 : SEGINI MPOVAD 84 : SEGINI MPOVAE 85 : 86 : C 87 : C- Pour chaque couple CELL/LAYERW on calcule les grandeurs liees a la 88 : C condensation 89 : C 90 : SEGACT IPT4 91 : N1 = IPT4.NUM(/2) 92 : DO 70 I=1,N1 93 : C 94 : C Evaluation du coefficient d'échange thermique lié à la condensation 95 : C 96 : IPOS = IND2(IPT4.NUM(2,I)) 97 : IF (IPOS.NE.0) THEN 98 : II1 = IND1(IPT4.NUM(1,I)) 99 : ROV = MPOVA1.VPOCHA(II1,1) 100 : RO = MPOVA2.VPOCHA(II1,1) 101 : ROH = MPOVA3.VPOCHA(II1,1) 102 : X1 = ROV / (RO - ROH) 103 : IF (ICOR.EQ.1) THEN 104 : H1 = 11.351D0 + 283.77D0*X1 105 : IF (TIME2.LE.TTAG) THEN 106 : H2 = ((TIME2/TTAG) ** 0.5) * TAGA 107 : H1 = H2 108 : ELSE 109 : H2 = TAGA 110 : E1 = EXP((TTAG-TIME2)/20.D0) 111 : H1 = (H2-H1)*E1 + H1 112 : ENDIF 113 : H1 = MIN(H1,HMAX) 114 : ELSE 115 : H1 = MIN(11.351D0 + 283.77D0*X1,HMAX) 116 : ENDIF 117 : C 118 : C Coeff d'échange : correction pour les murs en béton 119 : C Evaluation du débit de condensation 120 : C 121 : PTOT = MPOVA5.VPOCHA(II1,1) 122 : PVAP = MPOVA6.VPOCHA(II1,1) 123 : TSAT = MPOVA7.VPOCHA(II1,1) 124 : TGAS = MPOVA8.VPOCHA(II1,1) 125 : C IPOS = IND2(IPT4.NUM(2,I)) 126 : JPOS = IND3(IPT4.NUM(2,I)) 127 : TMUR = MPOVA9.VPOCHA(JPOS,1) 128 : AREA = MPOVAA.VPOCHA(IPOS,1) 129 : XCON = MPOVAB.VPOCHA(IPOS,1) 130 : HCON = MPOVA4.VPOCHA(IPOS,1) 131 : C 132 : C---> C'est ici que l'on fixe les échanges d'énergie comp/mur : 133 : C---> Etat thermodynamique du comp = PVAP TGAS 134 : C---> Etat thermodynamique du puisard = PTOT TMUR 135 : C---> La chaleur Latente découle de la difference comp/puisard 136 : C---> pour les états thermo ci-dessus CHOISIS A PRIORI. 137 : C 138 : HL = HLS0(PTOT,TMUR) 139 : HV = HVS0(PVAP,TGAS) 140 : HLAT = HV - HL 141 : C-------------------------------------------------------------------- 142 : MPOVAC.VPOCHA(IPOS,1) = MAX(XCON*H1,HCON) 143 : MPOVAE.VPOCHA(IPOS,1) = HL 144 : IF (TMUR.LT.TSAT) THEN 145 : DT = TSAT - TMUR 146 : IF (DT.LT.10.D0) THEN 147 : Z1 = 1.D0 - EXP(-1.D0*DT*DT) 148 : ELSE 149 : Z1 = 1.D0 150 : ENDIF 151 : DT2 = ABS(TGAS - TSAT) 152 : B2 = (BETA + (1.D0-BETA)*EXP(-1.D0*GAMA*DT2))*Z1 153 : MPOVAD.VPOCHA(IPOS,1) = -1.D0 * AREA * B2 / HLAT 154 : ELSE 155 : MPOVAD.VPOCHA(IPOS,1) = 0.D0 156 : ENDIF 157 : ENDIF 158 : 70 CONTINUE 159 : C 160 : C- Ménage 161 : C 162 : SEGSUP REDIR 163 : SEGDES IPT1,IPT2,IPT3,IPT4 164 : C 165 : RETURN 166 : END 167 : 168 : 169 : 170 :
© Cast3M 2003 - All rights reserved.
Disclaimer