1 : C FLA019 SOURCE TTMF3 12/07/05 21:15:42 7425 2 : SUBROUTINE FLA019(T,TC,MP,TS,XS,IDEB,TRAV) 3 : C--------------------------------------------------------------------- 4 : C Calcul du débit et des conditions en sortie de recombineur 5 : C (concentrations et température) 6 : C--------------------------------------------------------------------- 7 : C 8 : C--------------------------- 9 : C Parametres Entree/Sortie : 10 : C--------------------------- 11 : C 12 : C E/ T : flottant : Temps courant (s) 13 : C E/ TC : flottant : Température des plaques (K) 14 : C /S MP : flottant : Debit total au travers du PAR (kg/s) 15 : C /S TS : flottant : Température du gaz à la sortie du PAR (K) 16 : C /S XS : flottant : Fractions molaires à la sortie du PAR 17 : C /S IDEB : entier : Flag identifiant la situation 18 : C (négatif si le débit est négatif !) 19 : C 20 : C------------------------------ 21 : C Variables de TRAV utilisées : 22 : C------------------------------ 23 : C 24 : C E/ PRESSION : Pression à l'entrée du PAR (Pa) 25 : C E/ TEMPENT : Température à l'entrée du PAR (K) 26 : C E/ XiMOY : Fraction molaire à l'entrée du PAR 27 : C E/ M(7) : Masse molaire des constituants du mélange (kg/mol) 28 : C E/ EPS_MH2 : Seuil en deça duquel on suppose que la consommation 29 : C d'H2 max est mise à 0. (kg/s) 30 : C E/ EPS_DT : Ecart entre la température des plaques et celle du 31 : C gaz (Tc-T) en deça duquel on suppose qu'il n'y a pas 32 : C d'échange convectif (K) 33 : C E/ U : Vitesse minimale dans le PAR (= 0.01 m/s) 34 : C E/ S : Surface des plaques (m2) 35 : C E/ SP : Section de passage entrée/sortie fluide du PAR (m2) 36 : C E/ MC : Masse des plaques (kg) 37 : C E/ CPC : Chaleur spécifique des plaques (J/kg/K) 38 : C 39 : C------------------- 40 : C Autres variables : 41 : C------------------- 42 : C 43 : C MH2MAX : Consommation maximale par les plaques catalytique (kg/s) 44 : C MH2 : Consommation d'H2 par les plaques pour le débit MP (kg/s) 45 : C MH2EMAX : Débit d'hydrogène disponible à l'entrée du PAR (kg/s) 46 : C MO2EMAX : Débit d'hydrogène disponible à l'entrée du PAR limité par 47 : C le débit d'O2 insuffisant (kg/s) 48 : C 49 : C--------------------------------------------------------------------- 50 : C 51 : C Langage : ESOPE + FORTRAN 77 52 : C 53 : C Mise en oeuvre : H. Paillère (1997, TTMF) 54 : C 55 : C--------------------------------------------------------------------- 56 : IMPLICIT INTEGER(I-N) 57 : IMPLICIT REAL*8 (A-H,O-Z) 58 : REAL*8 XE(7),OME(7),XS(7),OMS(7) 59 : REAL*8 MH2MAX,MH2,MP,MH2EMAX,MO2EMAX 60 : segment trav 61 : integer iKALP,iMODEL 62 : real*8 e,L,Lch,Dh,S,sp,Ck 63 : real*8 mc,Cpc 64 : real*8 g,R,deltah 65 : real*8 M(nbesp),cstmod(ncst) 66 : real*8 Cpi(nbesp),al 67 : real*8 eps_mh2,eps_dt,eps_con,u 68 : real*8 XH2MOY,XO2MOY,XN2MOY,XH2OMOY,PRESSION,TEMPENT 69 : real*8 XHEMOY,XCO2MOY,XCOMOY 70 : endsegment 71 : C 72 : C OME(7) : Fraction massique 73 : PE = PRESSION 74 : TE = TEMPENT 75 : XE(1) = XN2MOY 76 : XE(2) = XO2MOY 77 : XE(3) = XH2MOY 78 : XE(4) = XH2OMOY 79 : XE(5) = XHEMOY 80 : XE(6) = XCO2MOY 81 : XE(7) = XCOMOY 81 : CALL FLA012(XE,OME,TRAV) 82 : ROE = FLA001(PE,TE,XE,TRAV) 83 : CPE = FLA003(OME,TRAV) 84 : C 85 : C Calcul de dmh2/dt sans limitation par le débit d'H2 86 : C et/ou d'O2 à l'entrée du PAR (kg/s) 87 : IF (IMODEL .EQ. 1) THEN 88 : MH2MAX = FLA008(PE,XE,CSTMOD(1),CSTMOD(2),CSTMOD(4),CSTMOD(5)) 89 : ELSEIF (IMODEL .EQ. 2) THEN 90 : MH2MAX = FLA008(PE,XE,CSTMOD(1),CSTMOD(2),CSTMOD(4),CSTMOD(5)) 91 : ENDIF 92 : C 93 : C 1. Cas sans comsommation d'hydrogène 94 : IF (MH2MAX.LT.EPS_MH2) THEN 95 : IF ((TC-TE) .LT. EPS_DT) THEN 96 : C 1.1 et sans echange convectif 97 : MH2 = 0.D0 98 : MP = U*ROE*SP 99 : TC = TE 100 : IDEB = 11 101 : ELSE 102 : C 1.2 et avec echange convectif (bilan d'énergie des plaques) 103 : MH2 = 0.D0 104 : MP = FLA015(T,TC,MH2,TRAV) 105 : IDEB = 12 106 : ENDIF 107 : ELSE 108 : C 2. Cas avec comsommation d'hydrogène 109 : IF ((TC-TE) .LT. EPS_DT) THEN 110 : C 2.1 et sans echange convectif mais prise en compte d'un déficit 111 : C éventuel en H2/O2 et débit minimum correspondant à la vitesse U 112 : BOBO = MH2MAX 113 : MH2EMAX = U * ROE * SP * OME(3) 114 : MO2EMAX = U * ROE * SP * 2.D0*M(3)*OME(2)/M(2) 115 : IF (MH2EMAX .LT. BOBO) THEN 116 : BOBO = MH2EMAX 117 : ENDIF 118 : IF (MO2EMAX .LT. BOBO) THEN 119 : BOBO = MO2EMAX 120 : ENDIF 121 : MH2 = BOBO 122 : MP = U*ROE*SP 123 : IDEB = 21 124 : ELSE 125 : C 2.2 et avec echange convectif (bilan d'énergie des plaques) 126 : C et prise en compte d'un déficit éventuel en H2/O2 par rapport 127 : C au débit mp traversant le PAR et correspondant à MH2MAX 128 : MP = FLA015(T,TC,MH2MAX,TRAV) 129 : C 2.2.1 Déficit en O2 (indic=2) ou déficit en H2 (indic=3) 130 : BNERF = OME(3) 131 : INDIC=3 132 : IF (BNERF .GT. 2.D0*M(3)*OME(2)/M(2)) THEN 133 : INDIC = 2 134 : BNERF = 2.D0*M(3)*OME(2)/M(2) 135 : ENDIF 136 : IF (MP*BNERF .LT. MH2MAX) THEN 137 : MP = FLA014(T,TC,INDIC,TRAV) 138 : MH2 = MP*BNERF 139 : IDEB = 221 140 : ELSE 141 : C 2.2.2 Pas de déficit en O2/H2 142 : MH2 = MH2MAX 143 : IDEB = 222 144 : ENDIF 145 : ENDIF 146 : ENDIF 147 : C 148 : C Température en sortie du PAR (bilan d'énergie du gaz) 149 : IF ((TC-TE).LT.EPS_DT) THEN 150 : TS = TE 151 : ELSE 152 : TS = TE + S*FLA009(PE,TE,TC,TE,XE,TRAV)*(TC-TE) 153 : & / (MP*CPE - AL*MH2) 154 : ENDIF 155 : C 156 : C Composition du mélange en sortie 157 : OMS(1) = OME(1) 158 : OMS(2) = OME(2) - M(2)*MH2/(2.D0*M(3)*MP) 159 : OMS(3) = OME(3) - MH2/MP 160 : OMS(4) = OME(4) + M(4)*MH2/(M(3)*MP) 161 : OMS(5) = OME(5) 162 : OMS(6) = OME(6) 163 : OMS(7) = OME(7) 163 : CALL FLA013(XS,OMS,TRAV) 164 : C 165 : C Conditions de réalisabilité 166 : XS(1) = ABS(XS(1)) 167 : XS(2) = ABS(XS(2)) 168 : XS(3) = ABS(XS(3)) 169 : XS(5) = ABS(XS(5)) 170 : XS(6) = ABS(XS(6)) 171 : XS(7) = ABS(XS(7)) 172 : XS(4) = 1.D0 - XS(1) - XS(2) - XS(3) - XS(5) - XS(6) - XS(7) 173 : C 174 : RETURN 175 : END 176 : 177 : 178 : 179 : 180 : 181 : 182 : 183 : 184 :
© Cast3M 2003 - All rights reserved.
Disclaimer