Download fla015.eso
Back to the list
1 : C FLA015 SOURCE TTMF3 12/07/05 21:15:37 7425
2 : FUNCTION FLA015(T,TC,MH2,TRAV)
3 : C---------------------------------------------------------------------
4 : C Calcul du débit au travers du PAR pour une température de plaque de
5 : C Tc et une consommation d'hydrogène mh2=mh2max
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/ MH2 : flottant : Hydrogène consommé (kg/s)
15 : C /S FLA015 : 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 E/ U : Vitesse minimale dans le PAR (= 0.01 m/s)
26 : C E/ S : Surface des plaques (m2)
27 : C E/ SP : Section de passage entrée/sortie fluide du PAR (m2)
28 : C E/ L : Hauteur des plaques (m)
29 : C E/ LCH : Hauteur de la cheminée (m)
30 : C E/ DH : Diamètre hydraulique entre les plaques (m)
31 : C E/ CK : Demi somme des coefficients de perte de charge à
32 : C l'entrée et à la sortie du PAR
33 : C E/ G : Gravité (=9.81 m/s2)
34 : C E/ AL : Cste fonction de CP(i)
35 : C E/ EPS_CON : Seuil de convergence pour Newton (|f(x)|CCOPTIO
61 : segment trav
62 : integer iKALP,iMODEL
63 : real*8 e,L,Lch,Dh,S,sp,Ck
64 : real*8 mc,Cpc
65 : real*8 g,R,deltah
66 : real*8 M(nbesp),cstmod(ncst)
67 : real*8 Cpi(nbesp),al
68 : real*8 eps_mh2,eps_dt,eps_con,u
69 : real*8 XH2MOY,XO2MOY,XN2MOY,XH2OMOY,PRESSION,TEMPENT
70 : real*8 XHEMOY,XCO2MOY,XCOMOY
71 : endsegment
72 : C
73 : MAXNEWT = 10000
74 : COFNEWT = 1.D0
75 : C
76 : PE = PRESSION
77 : TE = TEMPENT
78 : XE(1) = XN2MOY
79 : XE(2) = XO2MOY
80 : XE(3) = XH2MOY
81 : XE(4) = XH2OMOY
82 : XE(5) = XHEMOY
83 : XE(6) = XCO2MOY
84 : XE(7) = XCOMOY
84 : CALL FLA012(XE,OME,TRAV)
85 : ROE = FLA001(PE,TE,XE,TRAV)
86 : MUE = FLA004(TE,XE,TRAV)
87 : CPE = FLA003(OME,TRAV)
88 : ibou = cpi(/1)
89 : sum = 0.D0
90 : do 5 i=1,ibou
91 : sum = sum + xe(i)*m(i)
92 : 5 continue
93 : ME = sum
94 : C
95 : ROEXT = ROE
96 : C
97 : C Coefficient de l'équation de Bernouilli et de la densité en sortie
98 : A = Ck*COFNEWT*COFNEWT
99 : B = COFNEWT*48.D0*mue*sp*L/(Dh*Dh)
100 : C = -0.5D0*sp*sp*g*(L/2.D0+Lch)
101 : DD = S*fla009(pe,Te,Tc,Te,xe,trav)*(Tc-Te)/(Te*Cpe)*COFNEWT
102 : EE = al*mh2/Cpe*COFNEWT
103 : FF = mh2*Me/(2.D0*M(3))*COFNEWT
104 : C
105 : C Méthode de Dichotomie : la plus grande des racines est au-delà
106 : C des poles du denominateur (EE-DD) et FF
107 : I0 = 0
108 : XMIN = MAX(EE-DD,FF)
109 : IF (XMIN .LE. 0.D0) THEN
110 : XMIN = 0.D0
111 : ENDIF
112 : XMAX = MAX(1.D0,U*ROE*SP*500.D0)/COFNEWT
113 : RESUP = FLA016(XMAX,roext,A,B,C,DD,EE,FF)
114 : RESU = RESUP
115 : RESUM = FLA016(XMIN,roext,A,B,C,DD,EE,FF)
116 : EPS_2 = MIN(ABS(RESUM),EPS_CON)
117 : C
118 : C Recherche de la racine sur l'intervalle (XMIN,XMAX)
119 : 10 CONTINUE
120 : X = (XMAX + XMIN) / 2.D0
121 : IF (I0 .GE. MAXNEWT) GOTO 20
122 : IF ((XMAX-XMIN) .LT. EPS_CON) GOTO 30
123 : C IF (ABS(RESU) .LT. (EPS_CON*EPS_2)) GOTO 40
124 : CC IF (ABS(RESU) .LT. (RESUP*EPS_CON*EPS_2)) GOTO 40
125 : I0 = I0 + 1
126 : RESU = FLA016(x,roext,A,B,C,DD,EE,FF)
127 : IF ((RESU*RESUP) .LE. 0.D0) THEN
128 : XMIN = X
129 : ELSE
130 : XMAX = X
131 : ENDIF
132 : GOTO 10
133 : 20 CONTINUE
134 : C
135 : C I0 .GE. MAXNEWT : non convergence du Newton
136 : C Se produit à l'usage lorsque la racine et les poles sont quasi nuls.
137 : C Par suite, le débit est mis à zero.
138 : INTERR(1) = MAXNEWT
138 : CALL ERREUR(151)
139 : X = 0.D0
140 : GOTO 40
141 : 30 CONTINUE
142 : C
143 : C (XMAX-XMIN) .LT. EPS_CON : les 2 bornes de dichotomie sont égales.
144 : C Si la fonction est de même signe en XMAX et XMIN, il n'y a pas de
145 : C racine réelle au dela des poles. Par suite, le débit est nul.
146 : RESUP = FLA016(XMAX,roext,A,B,C,DD,EE,FF)
147 : RESUM = FLA016(XMIN,roext,A,B,C,DD,EE,FF)
148 : IF ((RESUM*RESUP) .GT. 0.D0) THEN
149 : X = 0.D0
150 : ENDIF
151 : 40 CONTINUE
152 : C
153 : C On maintient quoiqu'il arrive un débit résiduel dans le PAR
154 : C de vitesse débitante u
155 : MP = X * COFNEWT
156 : IF (U*ROE*SP .GT. MP) THEN
157 : MP = U*ROE*SP
158 : ENDIF
159 : FLA015 = MP
160 : C
161 : RETURN
162 : END
163 :
164 :
165 :
166 :
167 :
168 :