Mathematical Modeling and Optimization of Liquid-Junction ...
Mathematical Modeling and Optimization of Liquid-Junction ... Mathematical Modeling and Optimization of Liquid-Junction ...
259 C C C C FUNCTION CALCULATES COMPLETE ELLIPTIC INTEGRALS OF THE FIRST KIND WITH COMPLEMENTARY MODULUS AS A PARAMETER. IF(EMl.NE.O.) GO TO 10 PRINT 200 200 FORMAT( * PROGRAM STOPPED IN FUNCTION EK(EM1). * 1 /* COMPLEMENTARY MODULUS EQUAL TO ZERO.*) STOP 10 EK - 1.38629436112 + EMl*(0.09666344259 + EMl*(0.03590092383 + 1 EM1*(0.03742563713 + EM1*0.01451196212») 2 - ALOG(EMl)*(0.5 + EM1*(0.12498593597 + EMl*(0.06880248576 + 3 EM1*(0.03328355346 + EM1*0.00441787012»» RETURN END C C C C SUBROUTINE PRINT(NCOUNT) SUBROUTINE FOR MANIPULATING AND PRINTING RESULTS OF PROGRAM LJCMPY DIMENSION RXNMOD(14),CR(4) COMMON/SC/ PHI(250),P(250),AN(250),DFN(250),DFP(250),FN(250), 1 FP(250) COMMON/RUN/ EX,H,H2,HD2,EXH,NJJ,INDEX(7,250),KERR,ERRSUB,IJ COMMON/IN/ ALM,DELT,BD,SO,CURRNT,QII,SPOSN,TET,ANI2 COMMON/IN2/ BE,B3,RSC,GSC,RSOLN,GSOLN,GV,GT,GC,SI,DEL1,DEL3 COMMON/SOLI CBULK(4),Z(4),CSOL(4,250),F(4),HS,DC(4,250) COMMON/FLUX/ F3(250),F4(250) COMMON/IFC/ RATE(13),EQUIL(13),S(13,13),CONC(13),RXN(13),RXS(13) 1 ,CIF(5) COMMON/IFC2/ PHISS,PHISC,PHIHP,BKISS COMMON/SMQ/ SMQ,QI,QSL,SMQL COMMON/CHG/ SQL(250) COMMON/M/ MODE,MRPRNT,MREAL COMMON/CONV/ CL,CPOT,CCHG,COND COMMON/CE/ DIST,EXCOB,NREACT,CLIMJ,CLIM4 COMMON/OPT/ MOPT,SFLUX,ACE,RESIR COMMON/STOP/ NSTOP COMMON/END/ CURAN(40),VOLT(40),VNOT(40),VWIR(40),VWCE(40) 1 ,CHGISS(40),CHGIHP(40),CHG(40) 2 ,POWER1(40),POWER2(40),POWER3(40),POWER4(40) 190 FORMAT( * J*,10X,*PHI(J)*,10X,*CSOL(I)*,10X,*CR(I)* 1 ,6X,*F(I)/HS*, 8X,*I*) 195 FORMAT( /* INTERFACE ••• */ 15X,*PHISC*, 1 10X,*PHISS*, 9X,*PHIHP*) DO 2 1-1,4 F(I) - O. 2 CONTINUE
260 NJ - NJJ + 1 C C ZPOT LIQUID-JUNCTION VOLTAGE AT EQUILIBRIUM (DIMENSIONLESS) C DIST DISTANCE BETWEEN ELECTRODES (CM) C COND BULK ELECTROLYTE CONDUCTIVITY (1./0HM-CM) C EXCOB BUTLER-VOLMER EXCHANGE-CURRENT DENSITY FOR C REACTION AT COUNTERELECTRODE WITH NO MASS-TRANSFER C LIMITATIONS (MA/CM2) C NREACT NUMBER OF ELECTRONS TRANSFERED IN REACTION AT C COUNTERELECTRODE (EQUIV./MOLE) C CLIMJ,CLIM4 MASS-TRANSFER LIMITING CURRENTS FOR SPECIES 3 C AND 4 AT COUNTERELECTRODE (MA/CM2) C RMSC THE RATIO OF ELECTRON CONCENTRATIONS IN THE METAL C CURRENT-COLLECTOR TO THAT IN THE SEMICONDUCTOR C IF(NCOUNT.EQ.1) ZPOT - 45.6052*CPOT IF(NCOUNT.GT.l) GO TO 6 PRINT 140 PRINT 141,DIST,COND PRINT 142,EXCOB,NREACT PRINT 143,CLIMJ,CLIM4 140 FORMAT( /* ELECTROLYTE AND COUNTERELECTRODE CHARACTERISTICS*) 141 FORMAT( /* DIST -*F9.4* CM*15X*COND ~*F9.4* 1/0HM-CM*) 142 FORMAT( * EXCOB -*F9.4* MA/CM2*llX*NREACT ·*I4,5X* EQUIV./MOLE*) 143 FORMAT( * CLIMJ -*F9.4* MA/CM2*11X*CLIM4 -*F9.4* MA/CM2*/) 6 CONTINUE C C C CALCULATION OF TOTAL CELL VOLTAGE INCLUDING SEMICONDUCTOR, COUNTERELECTRODE, AND SOLUTION RESISTANCE CONTRIBUTIONS. C C VW - -(PHI(NJ) - PHISC + PIlISS - PHI(l»*CPOT IF(MODE.EQ.1) CURRENT - CURRNT*CL IF(MODE.EQ.2) CURRENT - (FP(IJ)/BD-FN(IJ»/H*CL IF(SO.NE.O.) GO TO 7 IF(CURRENT.NE.O.) GO TO 7 ZPOT - -vw 7 VIR - DIST*CURRENT/COND IF(MOPT.NE.O) VIR - CURRENT*RESIR IJJ - IJ - 1 C3 - -F3(IJJ)/HS/BE*CL C4 • -F4(IJJ)/HS/BE/B3*CL IF(CURRENT.EQ.O.) C3 - O. IF(CURRENT.EQ.O.) C4 - o. IF(MOPT.NE.O) C3 - C3*ACE IF(MOPT.NE.O) C4 - C4*ACE BETA - 0.50 IF«1.+e3/CLIMJ)/(1.+e4/CLIM4).LE.1.0E-10) PRINT 300 IF«1.+e3/CLIMJ)/(1.+C4/CLIM4).LE.1.0E-10) NSTOP ~ 1 IF«1.+C3/CLIMJ)/(1.+e4/CLIM4).LE.1.0E-10) RETURN ETACE - CPOT/(FLOAT(NREACT»*ALOG«1.+C3/CLIM3)/(1.+C4/CLIM4»
- Page 219 and 220: 208 C MSOL - 1; ILLUMINATION AT SEM
- Page 221 and 222: 210 C 15 IF(JCOUNT.GE.JCMAX) GO TO
- Page 223 and 224: 212 C C C IF(CC(13,1).LT.-DLIH2) CC
- Page 225 and 226: 214 CALL PRINT(NCOUNT) IF(JCOUNT.GE
- Page 227 and 228: 216 4 CONTINUE 5 CONTINUE IF(MDIM.G
- Page 229 and 230: 218 C C C C C C C GC GT GV DL1 DL2
- Page 231 and 232: 220 C DPOT2 - (DPOT3*EPSOL/DL3 - FA
- Page 233 and 234: 222 C C C RATE(lO) - RMOD(lO)*GISS*
- Page 235 and 236: 415 FORMAT( * LOW, INTERMEDIATE, AN
- Page 237 and 238: 226 24 PRINT 240 240 FORMAT( 1* SEM
- Page 239 and 240: 228 C C C C C C C C C C C C 10 SUBR
- Page 241 and 242: 230 COMMON/IN2/ BE,B3,RSC,GSC,RSOLN
- Page 243 and 244: 232 C C C C 00(1,6) - -1./4. BB(l,1
- Page 245 and 246: 234 C C C C C C C C C C C C C C C C
- Page 247 and 248: 236 GG(16)-ALOG(EQUIL(9)*CONC(10)*E
- Page 249 and 250: C C RXN(5) CHARGE-TRANSFER WITH LOW
- Page 251 and 252: 240 BB(24,24) - -1.0 2018 CONTINUE
- Page 253 and 254: 242 C C C C C C C C C IF(ABS(RXN(2
- Page 255 and 256: 244 C GM(6) - F3(J-l) - F3(J) AM(6,
- Page 257 and 258: 246 C C C C C C C C C C C C C CONTI
- Page 259 and 260: 248 C C C C C C C C C B(l,l) - -1.
- Page 261 and 262: 250 COMMON/BA/ N,NJ COMMON/LI/ AA(2
- Page 263 and 264: 252 8 IF (J-NJ) 11,9,9 9 DO 10 I a
- Page 265 and 266: 254 C C CALCULATION OF POTENTIAL DR
- Page 267 and 268: 256 lS/SOLN, SOLN/SEMI*) 202 FORMAT
- Page 269: 258 C C C IF (INDEX.GE.20) GO TO 35
- Page 273 and 274: 262 C C C CHG(NCOUNT) - SUMQ DPHDYZ
- Page 275 and 276: 264 RXNMOD(5) - RXN(5)/BD RXNMOD(6)
- Page 277 and 278: 266 0.2 1.00 1 30 41 41 40 000 1 o
- Page 279 and 280: 268 2. PROGRAM RCALe This program w
- Page 281 and 282: 270 C C C C C C PROGRAM RCALC (INPU
- Page 283 and 284: 272 C DZN - DZI*WOVRAZ DELBT - (CT-
- Page 285 and 286: 274 51 AZWDTH - CMPLX(AZN,O.) BZWDT
- Page 287 and 288: 276 U2C1 • U2C1 + 2.*HCT1 U2C2
- Page 289 and 290: 278 1 NATMAX - NATMAX*2 WOLD - WIMA
- Page 291 and 292: 280 DO 10,N-1,MAX TNM2 - TN TNM3 -
- Page 293 and 294: 282 IF(NPRINT.EQ.l) PRINT 220, RESC
- Page 295 and 296: 284 CURAVG - SUMAB/(DZIMAG-CZIMAG)
- Page 297 and 298: o 0 0 0 0 40 1.000e-06 1.500e+Ol 3.
259<br />
C<br />
C<br />
C<br />
C<br />
FUNCTION CALCULATES COMPLETE ELLIPTIC INTEGRALS OF THE FIRST<br />
KIND WITH COMPLEMENTARY MODULUS AS A PARAMETER.<br />
IF(EMl.NE.O.) GO TO 10<br />
PRINT 200<br />
200 FORMAT( * PROGRAM STOPPED IN FUNCTION EK(EM1). *<br />
1 /* COMPLEMENTARY MODULUS EQUAL TO ZERO.*)<br />
STOP<br />
10 EK - 1.38629436112 + EMl*(0.09666344259 + EMl*(0.03590092383 +<br />
1 EM1*(0.03742563713 + EM1*0.01451196212»)<br />
2 - ALOG(EMl)*(0.5 + EM1*(0.12498593597 + EMl*(0.06880248576 +<br />
3 EM1*(0.03328355346 + EM1*0.00441787012»»<br />
RETURN<br />
END<br />
C<br />
C<br />
C<br />
C<br />
SUBROUTINE PRINT(NCOUNT)<br />
SUBROUTINE FOR MANIPULATING AND PRINTING RESULTS OF PROGRAM<br />
LJCMPY<br />
DIMENSION RXNMOD(14),CR(4)<br />
COMMON/SC/ PHI(250),P(250),AN(250),DFN(250),DFP(250),FN(250),<br />
1 FP(250)<br />
COMMON/RUN/ EX,H,H2,HD2,EXH,NJJ,INDEX(7,250),KERR,ERRSUB,IJ<br />
COMMON/IN/ ALM,DELT,BD,SO,CURRNT,QII,SPOSN,TET,ANI2<br />
COMMON/IN2/ BE,B3,RSC,GSC,RSOLN,GSOLN,GV,GT,GC,SI,DEL1,DEL3<br />
COMMON/SOLI CBULK(4),Z(4),CSOL(4,250),F(4),HS,DC(4,250)<br />
COMMON/FLUX/ F3(250),F4(250)<br />
COMMON/IFC/ RATE(13),EQUIL(13),S(13,13),CONC(13),RXN(13),RXS(13)<br />
1 ,CIF(5)<br />
COMMON/IFC2/ PHISS,PHISC,PHIHP,BKISS<br />
COMMON/SMQ/ SMQ,QI,QSL,SMQL<br />
COMMON/CHG/ SQL(250)<br />
COMMON/M/ MODE,MRPRNT,MREAL<br />
COMMON/CONV/ CL,CPOT,CCHG,COND<br />
COMMON/CE/ DIST,EXCOB,NREACT,CLIMJ,CLIM4<br />
COMMON/OPT/ MOPT,SFLUX,ACE,RESIR<br />
COMMON/STOP/ NSTOP<br />
COMMON/END/ CURAN(40),VOLT(40),VNOT(40),VWIR(40),VWCE(40)<br />
1 ,CHGISS(40),CHGIHP(40),CHG(40)<br />
2 ,POWER1(40),POWER2(40),POWER3(40),POWER4(40)<br />
190 FORMAT( * J*,10X,*PHI(J)*,10X,*CSOL(I)*,10X,*CR(I)*<br />
1 ,6X,*F(I)/HS*, 8X,*I*)<br />
195 FORMAT( /* INTERFACE ••• */ 15X,*PHISC*,<br />
1 10X,*PHISS*, 9X,*PHIHP*)<br />
DO 2 1-1,4<br />
F(I) - O.<br />
2 CONTINUE