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