Download fla014.eso
Back to the list
1 : C FLA014 SOURCE TTMF3 12/07/05 21:15:36 7425
2 : FUNCTION FLA014(T,TC,INDIC,TRAV)
3 : C---------------------------------------------------------------------
4 : C Calcul du débit au travers du PAR pour une température de plaque de
5 : C TC et un controle de mh2 par un déficit en O2/H2 à l'entrée du PAR
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 E/ INDIC : entier : Déficit en O2 (2) ou en H2 (3)
15 : C /S FLA014 : flottant : Débit total traversant le PAR (kg/s)
16 : C
17 : C------------------------------
18 : C Variables de TRAV utilisées :
19 : C------------------------------
20 : C
21 : C E/ PRESSION : Pression à l'entrée du PAR (Pa)
22 : C E/ TEMPENT : Température à l'entrée du PAR (K)
23 : C E/ XiMOY : Fraction molaire à l'entrée du PAR
24 : C E/ M(7) : Masse molaire des constituants du mélange (kg/mol)
25 : C (dans l'ordre N2, O2, H2 et H2O)
26 : C E/ U : Vitesse minimale dans le PAR (= 0.01 m/s)
27 : C E/ S : Surface des plaques (m2)
28 : C E/ SP : Section de passage entrée/sortie fluide du PAR (m2)
29 : C E/ L : Hauteur des plaques (m)
30 : C E/ LCH : Hauteur de la cheminée (m)
31 : C E/ DH : Diamètre hydraulique entre les plaques (m)
32 : C E/ CK : Demi somme des coefficients de perte de charge à
33 : C l'entrée et à la sortie du PAR
34 : C E/ G : Gravité (=9.81 m/s2)
35 : C E/ AL : Cste fonction de CP(i)
36 : C E/ EPS_CON : Seuil de convergence pour Newton (|f(x)|CCOPTIO
63 : segment trav
64 : integer iKALP,iMODEL
65 : real*8 e,L,Lch,Dh,S,sp,Ck
66 : real*8 mc,Cpc
67 : real*8 g,R,deltah
68 : real*8 M(nbesp),cstmod(ncst)
69 : real*8 Cpi(nbesp),al
70 : real*8 eps_mh2,eps_dt,eps_con,u
71 : real*8 XH2MOY,XO2MOY,XN2MOY,XH2OMOY,PRESSION,TEMPENT
72 : real*8 XHEMOY,XCO2MOY,XCOMOY
73 : endsegment
74 : C
75 : MAXNEWT = 10000
76 : C
77 : PE = PRESSION
78 : TE = TEMPENT
79 : XE(1) = XN2MOY
80 : XE(2) = XO2MOY
81 : XE(3) = XH2MOY
82 : XE(4) = XH2OMOY
83 : XE(5) = XHEMOY
84 : XE(6) = XCO2MOY
85 : XE(7) = XCOMOY
85 : CALL FLA012(XE,OME,TRAV)
86 : ROE = FLA001(PE,TE,XE,TRAV)
87 : MUE = FLA004(TE,XE,TRAV)
88 : CPE = FLA003(OME,TRAV)
89 : C
90 : ibou = cpi(/1)
91 : sum = 0.D0
92 : do 10 i=1,ibou
93 : sum = sum + xe(i)*m(i)
94 : 10 continue
95 : ME = sum
96 : C
97 : C Coefficient de la partie linéaire en x
98 : C = -0.5D0*sp*sp*g*(L/2.D0+Lch)
99 : A = -Ck/C*1.D-6
100 : B = -48.D0*mue*sp*L/(Dh*Dh)/C*1.D-3
101 : C
102 : C Coefficient de la partie non-linéaire en x
103 : C Initialisation de DD et EE avec mh2 controlé par le déficit en H2 à
104 : C l'entrée du PAR (indic=3) (resp. par le déficit en O2, indic=2)
105 : IF (INDIC .EQ. 3) THEN
106 : DD = 1.D3*S*fla009(pe,Te,Tc,Te,xe,trav)*
107 : & (Tc-Te)/(Te*(Cpe-al*ome(3)))
108 : EE = 1.D0 - Me*ome(3)/(2*M(3))
109 : ELSE
110 : DD = 1.D3*S*fla009(pe,Te,Tc,Te,xe,trav)*(Tc-Te)/
111 : & (Te*(Cpe-2*al*M(3)*ome(2)/M(2)))
112 : EE = 1.D0 - Me*ome(2)/M(2)
113 : ENDIF
114 : C
115 : C Méthode de Newton : x_(k+1) = x_k - f(x_k)/f'(x_k)
116 : C (l'initialisation X=DD correspond d'après l'équation de bilan
117 : C d'énergie du gaz au débit qu'on aurait pour une température
118 : C en sortie du PAR égale à 2*TE)
119 : RESU = 1.D0
120 : X = DD
121 : I0 = 0
122 : 200 CONTINUE
123 : IF (I0 .GE. MAXNEWT) GOTO 201
124 : IF ((ABS(RESU)) .LT. (EPS_CON*100.D0)) GOTO 201
125 : I0 = I0 + 1
125 : CALL FLA017(x,y,yp,roe,A,B,DD,EE)
126 : RESU = Y
127 : X = X - Y/YP
128 : GOTO 200
129 : 201 CONTINUE
130 : C
131 : C On maintient quoiqu'il arrive un débit résiduel dans le PAR
132 : C de vitesse débitante u
133 : IF (I0 .GE. MAXNEWT) THEN
134 : WRITE(IOIMP,*) 'fla014 : non convergence : ',(X*1.D-3),TC
135 : INTERR(1) = MAXNEWT
135 : CALL ERREUR(151)
136 : X = 0.D0
137 : ENDIF
138 : BOBO = X * 1.D-3
139 : IF (U*ROE*SP .GT. BOBO) THEN
140 : BOBO = U*ROE*SP
141 : ENDIF
142 : FLA014 = BOBO
143 : C
144 : RETURN
145 : END
146 :
147 :
148 :
149 :
150 :
151 :
152 :
153 :
154 :
155 :