TALK=T;RUN(1,1) PHOTON USE p 2 2 0.9 up z use patgeo vec x 3 sh pause;cl use patgeo con tmp1 z 1 fil;.001 con tmp1 x 2 fil;.001 con tmp1 x 4 fil;.001 con tmp1 z m fil;.001 con tmp1 z 10 fil;.001 vec x 3 msg(Gas temperature contours pause;cl use patgeo con yhcn z 1 fil;.001 con yhcn x 2 fil;.001 con yhcn x 4 fil;.001 con yhcn z m fil;.001 con yhcn z 10 fil;.001 vec x 3 msg(HCN mass fractions pause;cl use patgeo con ynox z 1 fil;.001 con ynox x 2 fil;.001 con ynox x 4 fil;.001 con ynox z m fil;.001 con ynox z 10 fil;.001 vec x 3 msg(NOX mass fractions pause;cl use patgeo con tmp2 z 1 fil;.001 con tmp2 x 2 fil;.001 con tmp2 x 4 fil;.001 con tmp2 z m fil;.001 con tmp2 z 10 fil;.001 vec x 3 msg(Solid temperature contours pause;cl use patgeo con VRAT z 1 fil;.001 con VRAT x 2 fil;.001 con VRAT x 4 fil;.001 con VRAT z m fil;.001 con VRAT z 10 fil;.001 vec x 3 msg(RATe of Volatalization, kg/s pause;cl use patgeo con VCOM z 1 fil;.001 con VCOM x 2 fil;.001 con VCOM x 4 fil;.001 con VCOM z m fil;.001 con VCOM z 10 fil;.001 vec x 3 msg(Rate of Volatiles, CH4, COMbustion, kg/s pause;cl use patgeo con CCOM z 1 fil;.001 con CCOM x 2 fil;.001 con CCOM x 4 fil;.001 con CCOM z m fil;.001 con CCOM z 10 fil;.001 vec x 3 msg(Rate of CO COMbustion, kg/s pause;cl use patgeo con BURN z 1 fil;.001 con BURN x 2 fil;.001 con BURN x 4 fil;.001 con BURN z m fil;.001 con BURN z 10 fil;.001 vec x 3 msg(Rate of char BURNing, kg/s pause;cl use patgeo con YCH4 z 1 fil;.001 con YCH4 x 2 fil;.001 con YCH4 x 4 fil;.001 con YCH4 z m fil;.001 con YCH4 z 10 fil;.001 vec x 3 msg(Volatiles, CH4, mass fractions pause;cl use patgeo con YCO2 z 1 fil;.001 con YCO2 x 2 fil;.001 con YCO2 x 4 fil;.001 con YCO2 z m fil;.001 con YCO2 z 10 fil;.001 vec x 3 msg(CO2 mass fractions pause;cl use patgeo con YCO z 1 fil;.001 con YCO x 2 fil;.001 con YCO x 4 fil;.001 con YCO z m fil;.001 con YCO z 10 fil;.001 vec x 3 msg(CO mass fractions pause;cl use patgeo con YH2O z 1 fil;.001 con YH2O x 2 fil;.001 con YH2O x 4 fil;.001 con YH2O z m fil;.001 con YH2O z 10 fil;.001 vec x 3 msg(H2O mass fractions pause;cl use patgeo con YO2 z 1 fil;.001 con YO2 x 2 fil;.001 con YO2 x 4 fil;.001 con YO2 z m fil;.001 con YO2 z 10 fil;.001 vec x 3 msg(O2 mass fractions pause;cl use patgeo con YN2 z 1 fil;.001 con YN2 x 2 fil;.001 con YN2 x 4 fil;.001 con YN2 z m fil;.001 con YN2 z 10 fil;.001 vec x 3 msg(N2 mass fractions pause;cl use patgeo con COL2 z 1 fil;.001 con COL2 x 2 fil;.001 con COL2 x 4 fil;.001 con COL2 z m fil;.001 con COL2 z 10 fil;.001 vec x 3 msg(COAL particle content pause;cl use patgeo con CHA2 z 1 fil;.001 con CHA2 x 2 fil;.001 con CHA2 x 4 fil;.001 con CHA2 z m fil;.001 con CHA2 z 10 fil;.001 vec x 3 msg(CHAR particle content pause;cl use patgeo con ASH2 z 1 fil;.001 con ASH2 x 2 fil;.001 con ASH2 x 4 fil;.001 con ASH2 z m fil;.001 con ASH2 z 10 fil;.001 vec x 3 msg(ASH particle content pause;cl use patgeo con R2 z 1 fil;.001 con R2 x 2 fil;.001 con R2 x 4 fil;.001 con R2 z m fil;.001 con R2 z 10 fil;.001 vec x 3 msg(Particle volume fractions pause;cl use patgeo con VREL z 1 fil;.001 con VREL x 2 fil;.001 con VREL x 4 fil;.001 con VREL z m fil;.001 con VREL z 10 fil;.001 vec x 3 msg(Particle relative velocity pause;cl use patgeo con SIZE z 1 fil;.001 con SIZE x 2 fil;.001 con SIZE x 4 fil;.001 con SIZE z m fil;.001 con SIZE z 10 fil;.001 vec x 3 msg(Particle diameter pause;cl use patgeo con REYN z 1 fil;.001 con REYN x 2 fil;.001 con REYN x 4 fil;.001 con REYN z m fil;.001 con REYN z 10 fil;.001 vec x 3 msg(Particle Reynolds number pause;cl use patgeo con CD z 1 fil;.001 con CD x 2 fil;.001 con CD x 4 fil;.001 con CD z m fil;.001 con CD z 10 fil;.001 vec x 3 msg(Particle drag coefficient pause;cl use patgeo con NUSS z 1 fil;.001 con NUSS x 2 fil;.001 con NUSS x 4 fil;.001 con NUSS z m fil;.001 con NUSS z 10 fil;.001 vec x 3 msg(Particle Nusselt number pause;cl use patgeo con HCOF z 1 fil;.001 con HCOF x 2 fil;.001 con HCOF x 4 fil;.001 con HCOF z m fil;.001 con HCOF z 10 fil;.001 vec x 3 msg(Interphase heat transfer coefficient pause;cl use patgeo con DEN1 z 1 fil;.001 con DEN1 x 2 fil;.001 con DEN1 x 4 fil;.001 con DEN1 z m fil;.001 con DEN1 z 10 fil;.001 vec x 3 msg(Gas density pause;cl use patgeo con EPKE z 1 fil;.001 con EPKE x 2 fil;.001 con EPKE x 4 fil;.001 con EPKE z m fil;.001 con EPKE z 10 fil;.001 vec x 3 msg(Reciprocal time scale. 1/s pause;cl ENDUSE DISPLAY FUEL-DUST FLAMES in a FURNACE Wall-fired furnace: 24 burners Dry coal. 1 Kinetically released volatiles; 2 Combustion of: * Volatiles, CH4, by EBU with a provison for kinetic ; * Char, C(s), by harmonic blend of kinetic and diffusion and * Carbon monoxide, CO, by harmonic blend of kinetic and EBU. 3 Thermal NOX and 4 Fuel NOX formations. 5 Radiosity model for radiating particles ENDDIS PLANTBEGIN NAMSAT=MOSG PRNDTL(CRAD)=GRNDLAMPR(CRAD)=:ENUL:*DEN1*3./4.*(0.5+5.e-1*6.*R2/SIZE) * Interphase mass transfer VRAT-volatile release; BURN-char burning INTMDT=VRAT+BURN * Phase temperatures TEMP1=H1/(1059.+0.25*(TMP1-300.)) TEMP2=H2/1800. ========================================= * Devolatilization related sources/sinks # Reaction: 1kg Coal -> Ykg Vol + (1-Y)kg Char PATCH(PLSR1,CELL,1,NX,1,NY,1,NZ,1,LSTEP) ** Coal VAL=-VRAT/RHO2A COVAL(PLSR1,COL2,FIXFLU,GRND) ** Char VAL=VRAT*(1.-RHO2A)/RHO2A COVAL(PLSR1,CHA2,FIXFLU,GRND) ====================================== CH4-volatile oxidation * YCH4 supply via volatile release PATCH(PLSR5,CELL,1,NX,1,NY,1,NZ,1,LSTEP) VAL=VRAT*1.0 COVAL(PLSR5,YCH4,FIXFLU,GRND) Reaction: 2CH4 + 3O2 -> 2CO +4H2O PATCH(PLSR6,CELL,1,NX,1,NY,1,NZ,1,LSTEP) CO=96./32*VCOM/(YO2+tiny) COVAL(PLSR6,YO2,GRND,0.0) CO=VCOM/(YCH4+tiny) COVAL(PLSR6,YCH4,GRND,0.0) VAL=72./32.*VCOM COVAL(PLSR6,YH2O,FIXFLU,GRND) VAL=56./32.*VCOM COVAL(PLSR6,YCO,FIXFLU,GRND) ========================================== Carbon-Char oxidation Reaction: C(s) + (x+1)/2*O2 -> xCO2 + (1-x)CO * Sink of char,C(s), in a 2nd phase PATCH(PLSR8,CELL,1,NX,1,NY,1,NZ,1,LSTEP) VAL=-BURN COVAL(PLSR8,CHA2,FIXFLU,GRND) * Sink/sources of 1st phase components PATCH(PLSR9,CELL,1,NX,1,NY,1,NZ,1,LSTEP) VAL=44./12.*XXX*BURN COVAL(PLSR9,YCO2,FIXFLU,GRND) VAL=28./12.*(1.-XXX)*BURN COVAL(PLSR9,YCO,FIXFLU,GRND) VAL=-32./12.*(1.+XXX)/2.*BURN COVAL(PLSR9,YO2,FIXFLU,GRND) =========================================== Carbon monoxide oxidation Reaction: 2CO + O2 -> 2CO2 PATCH(PLSR7,CELL,1,NX,1,NY,1,NZ,1,LSTEP) CO=CCOM/(YCO+tiny) COVAL(PLSR7,YCO,GRND,0.0) CO=CCOM*32./56./(YO2+tiny) COVAL(PLSR7,YO2,GRND,0.0) VAL=CCOM*88./56. COVAL(PLSR7,YCO2,FIXFLU,GRND) ========================================== * Enthalpy sources due to heats of combustion # Volatiles PATCH(PLSR3,CELL,1,NX,1,NY,1,NZ,1,LSTEP) VAL=VCOM*5.5e7+CCOM*7.45e6 COVAL(PLSR3,H1,FIXFLU,GRND) # Char PATCH(PLSR10,CELL,1,NX,1,NY,1,NZ,1,LSTEP) VAL=BURN*(XXX*3.28e7+(1.-XXX)*9.21e6) COVAL(PLSR10,H1,FIXFLU,GRND) * Enthalpy sources due to interphase heat transfer PATCH(PLSR4,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP) CO=HCOF/(1059.+0.25*(TMP1-300.)) VAL=H2*(1059.+0.25*(TMP1-300.)) /1800. COVAL(PLSR4,H1,GRND,GRND) CO=HCOF/1800. VAL=1800./(1059.+0.25*(TMP1-300.)) *H1 COVAL(PLSR4,H2,GRND,GRND) * Enthalpy sources due to radiation PATCH(RADSOR,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP) VAL=4.*0.5*(CRAD-5.68e-8*TMP1**4) COVAL(RADSOR,H1,FIXFLU,GRND) VAL=4.*5.e-1*6.*R2/SIZE*(CRAD-5.68e-8*TMP2**4) COVAL(RADSOR,H2,FIXFLU,GRND) * Radiosity source PATCH(RDIOSITY,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP) CO=4.*(5.e-1*6.*R2/SIZE+0.5) VAL=5.68e-8*(5.e-1*6.*R2/SIZE*TMP2**4+0.5*TMP1**4 )/$ (5.e-1*6.*R2/SIZE+0.5) COVAL(RDIOSITY,CRAD,GRND,GRND) ============================================ * Sources to counter the transfer to phase 1 PATCH(PLSR90,CELL,1,NX,1,NY,1,NZ,1,LSTEP) VAL=CMDO*COL2 COVAL(PLSR90,COL2,FIXFLU,GRND) VAL=CMDO*CHA2 COVAL(PLSR90,CHA2,FIXFLU,GRND) =========================================== PATCH(THERMNOX,VOLUME,1,NX,1,NY,1,NZ,1,1) VAL=RNOX COVAL(THERMNOX,YNOX,FIXFLU,GRND) PATCH(FUELHCN,CELL,1,NX,1,NY,1,NZ,1,LSTEP) VAL=27./18.*(1.e-2*VRAT+1.e-2*BURN) COVAL(FUELHCN,YHCN,FIXFLU,GRND) PATCH(FUELNOX,PHASEM,1,NX,1,NY,1,NZ,1,LSTEP) CO=3.e12*YHCN/EXP(2.5e8/(8314.*AMAX1(300.,TMP1))) VAL=0.0333*YO2 COVAL(FUELNOX,YNOX,GRND,GRND) PATCH(NOXRED,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP) CO=6.*R2/SIZE*40./EXP(1.5e8/(8314.*AMAX1(300.,TMP2))) COVAL(NOXRED,YNOX,GRND,0.0) =========================================== store(HCOF,VRAT,VCOM,VCOD,VCOK,CCOM,CCOD,CCOK,WMIX,KD,KC,KT,BURN) store(BRN1,BRN2) # Char burning rate constants: ** Mixture molecular mass WMIX=28.8 ** Diffusion KD=1.55*NUSS*1.78e-5*12./8314./AMAX1(300.,TMP2)/SIZE ** Kinetic KC=0.1309/EXP(26850./4./AMAX1(300.,TMP2)) ** Harmonic blending KT=1./(1./KD+1./KC) # Rates: # Volatile release, kg/s VRAT=2.e3*RHO2A*COL2/EXP(2829./AMAX1(300.,TMP2))*$ RHO2*R2*VOL*RHO2A # Volatile combustion, kg/s ** EBU VCOD=4.*EPKE*AMIN1(YCH4,32.*YO2/96.) ** Kinetic VCOK=1.15e9*DEN1*YO2**1.3/EXP(24444./AMAX1(300.,TMP1))/$ (YCH4**0.3+tiny) ** Harmonic blending VCOM=1./(1./VCOD+1./VCOK)*VOL*R1*DEN1 # Carbon monoxide,CO, combustion, kg/s ** EBU CCOD=4.*EPKE*AMIN1(YCO,32.*YO2/56.) ** Kinetic CCOK=5.42e9*DEN1*YCO*YO2**0.25*YH2O**0.5/ E$ XP(15152./AMAX1(300.,TMP1)) ** Harmonic blending CCOM=1./(1./CCOD+1./CCOK)*VOL*R1*DEN1 # Char,C, burning out, kg/s BRN1=6./SIZE*KT*(P1+PRESS0)*YO2*WMIX/32.*R2*VOL BRN2=4.*DEN1*R1*VOL*EPKE*AMIN1(CHA2*MASS2/MASS1,$ 12./(32.*(1.+XXX)/2.)*YO2) BURN=1./(1./BRN1+1./BRN2) # CO/CO2 of C-combustion store(MCCO) MCCO=2500.*EXP(-6249./AMAX1(300.,TMP2)) XXX=1./(1.+44./28.*MCCO) # Interphase heat transfer coefficient NUSS=2.+0.65*REYN**0.5 HCOF=6.*R2*0.0458*NUSS/SIZE**2. # Ash contents of 2nd phase ASH2=1.-CHA2-COL2 # Nitrogen contents of 1st phase YN2=1.-YCH4-YO2-YCO-YCO2-YH2O store(FK1,MOX,RNOX) *** Reaction constant FK1=1.8e8*EXP(-38370./AMAX1(TMP1,300.)) *** Mole O atom concentration MOX=3.97e5*(YO2*DEN1/32./AMAX1(TMP1,300.))**.5 *EXP($ -31090./AMAX1(TMP1,300.)) *** NOX reaction rate RNOX=30.*2.*FK1*MOX*YN2*DEN1/28. PLANTEND ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(Wall-fired model furnace ) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 2. Time dependence STEADY = T ************************************************************ Group 3. X-Direction Grid Spacing CARTES = T NX = 5 XULAST =1. XFRAC(1)=0.2 ;XFRAC(2)=0.4 XFRAC(3)=0.6 ;XFRAC(4)=0.8 XFRAC(5)=1. ************************************************************ Group 4. Y-Direction Grid Spacing NY = 5 YVLAST =1. YFRAC(1)=0.2 ;YFRAC(2)=0.4 YFRAC(3)=0.6 ;YFRAC(4)=0.8 YFRAC(5)=1. ************************************************************ Group 5. Z-Direction Grid Spacing PARAB = F NZ = 20 ZWLAST =10. ZFRAC(1)=0.05 ;ZFRAC(2)=0.1 ZFRAC(3)=0.15 ;ZFRAC(4)=0.2 ZFRAC(5)=0.25 ;ZFRAC(6)=0.3 ZFRAC(7)=0.35 ;ZFRAC(8)=0.4 ZFRAC(9)=0.45 ;ZFRAC(10)=0.5 ZFRAC(11)=0.55 ;ZFRAC(12)=0.6 ZFRAC(13)=0.65 ;ZFRAC(14)=0.7 ZFRAC(15)=0.75 ;ZFRAC(16)=0.8 ZFRAC(17)=0.85 ;ZFRAC(18)=0.9 ZFRAC(19)=0.95 ;ZFRAC(20)=1. ************************************************************ Group 6. Body-Fitted Coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = F EQUVEL = F NAME(1)=P1 ;NAME(3)=U1 NAME(4)=U2 ;NAME(5)=V1 NAME(6)=V2 ;NAME(7)=W1 NAME(8)=W2 ;NAME(9)=R1 NAME(10)=R2 ;NAME(11)=RS NAME(12)=KE ;NAME(13)=EP NAME(14)=H1 ;NAME(15)=H2 NAME(105)=CRAD ;NAME(106)=MCCO NAME(107)=KT ;NAME(108)=XXX NAME(109)=BRN2 ;NAME(110)=BRN1 NAME(111)=RNOX ;NAME(112)=MOX NAME(113)=FK1 ;NAME(114)=CCOK NAME(115)=CCOD ;NAME(116)=VCOD NAME(117)=BURN ;NAME(118)=WMIX NAME(119)=CCOM ;NAME(120)=VCOM NAME(121)=YN2 ;NAME(122)=ASH2 NAME(123)=VCOK ;NAME(124)=YH2O NAME(125)=YCO2 ;NAME(126)=YCO NAME(127)=YO2 ;NAME(128)=YCH4 NAME(129)=HCOF ;NAME(130)=KC NAME(131)=KD ;NAME(132)=SIZE NAME(133)=DEN1 ;NAME(134)=APRJ NAME(135)=COL2 ;NAME(136)=CHA2 NAME(137)=REYN ;NAME(138)=CD NAME(139)=VREL ;NAME(140)=NUSS NAME(141)=TMP2 ;NAME(142)=TMP1 NAME(143)=CFIP ;NAME(144)=VRAT NAME(145)=CMDO ;NAME(146)=EL1 NAME(147)=ENUT ;NAME(148)=EPKE NAME(149)=YHCN ;NAME(150)=YNOX * Y in SOLUTN argument list denotes: * 1-stored 2-solved 3-whole-field * 4-point-by-point 5-explicit 6-harmonic averaging SOLUTN(P1,Y,Y,Y,N,N,Y) SOLUTN(U1,Y,Y,N,N,N,Y) SOLUTN(U2,Y,Y,N,N,N,Y) SOLUTN(V1,Y,Y,N,N,N,Y) SOLUTN(V2,Y,Y,N,N,N,Y) SOLUTN(W1,Y,Y,N,N,N,Y) SOLUTN(W2,Y,Y,N,N,N,Y) SOLUTN(R1,Y,Y,N,Y,N,Y) SOLUTN(R2,Y,Y,N,Y,N,Y) SOLUTN(RS,Y,Y,N,N,N,Y) SOLUTN(KE,Y,Y,N,N,N,N) SOLUTN(EP,Y,Y,N,N,N,N) SOLUTN(H1,Y,Y,N,N,N,Y) SOLUTN(H2,Y,Y,N,N,N,Y) SOLUTN(CRAD,Y,Y,N,N,N,Y) SOLUTN(MCCO,Y,N,N,N,N,Y) SOLUTN(KT,Y,N,N,N,N,Y) SOLUTN(XXX,Y,N,N,N,N,Y) SOLUTN(BRN2,Y,N,N,N,N,Y) SOLUTN(BRN1,Y,N,N,N,N,Y) SOLUTN(RNOX,Y,N,N,N,N,Y) SOLUTN(MOX,Y,N,N,N,N,Y) SOLUTN(FK1,Y,N,N,N,N,Y) SOLUTN(CCOK,Y,N,N,N,N,Y) SOLUTN(CCOD,Y,N,N,N,N,Y) SOLUTN(VCOD,Y,N,N,N,N,Y) SOLUTN(BURN,Y,N,N,N,N,Y) SOLUTN(WMIX,Y,N,N,N,N,Y) SOLUTN(CCOM,Y,N,N,N,N,Y) SOLUTN(VCOM,Y,N,N,N,N,Y) SOLUTN(YN2,Y,N,N,N,N,Y) SOLUTN(ASH2,Y,N,N,N,N,Y) SOLUTN(VCOK,Y,N,N,N,N,Y) SOLUTN(YH2O,Y,Y,N,N,N,Y) SOLUTN(YCO2,Y,Y,N,N,N,Y) SOLUTN(YCO,Y,Y,N,N,N,Y) SOLUTN(YO2,Y,Y,N,N,N,Y) SOLUTN(YCH4,Y,Y,N,N,N,Y) SOLUTN(HCOF,Y,N,N,N,N,Y) SOLUTN(KC,Y,N,N,N,N,Y) SOLUTN(KD,Y,N,N,N,N,Y) SOLUTN(SIZE,Y,N,N,N,N,Y) SOLUTN(DEN1,Y,N,N,N,N,Y) SOLUTN(APRJ,Y,N,N,N,N,Y) SOLUTN(COL2,Y,Y,N,N,N,Y) SOLUTN(CHA2,Y,Y,N,N,N,Y) SOLUTN(REYN,Y,N,N,N,N,Y) SOLUTN(CD,Y,N,N,N,N,Y) SOLUTN(VREL,Y,N,N,N,N,Y) SOLUTN(NUSS,Y,N,N,N,N,Y) SOLUTN(TMP2,Y,N,N,N,N,Y) SOLUTN(TMP1,Y,N,N,N,N,Y) SOLUTN(CFIP,Y,N,N,N,N,Y) SOLUTN(VRAT,Y,N,N,N,N,Y) SOLUTN(CMDO,Y,N,N,N,N,Y) SOLUTN(EL1,Y,N,N,N,N,Y) SOLUTN(ENUT,Y,N,N,N,N,Y) SOLUTN(EPKE,Y,N,N,N,N,Y) SOLUTN(YHCN,Y,Y,N,N,N,Y) SOLUTN(YNOX,Y,Y,N,N,N,Y) DEN1 = 133 VIST = 147 TEMP1 = 142 TEMP2 = 141 LEN1 = 146 ************************************************************ Group 8. Terms & Devices * Y in TERMS argument list denotes: * 1-built-in source 2-convection 3-diffusion 4-transient * 5-first phase variable 6-interphase transport TERMS(P1,Y,Y,Y,N,Y,Y) TERMS(U1,Y,Y,Y,Y,Y,Y) TERMS(U2,Y,Y,Y,Y,N,Y) TERMS(V1,Y,Y,Y,Y,Y,Y) TERMS(V2,Y,Y,Y,Y,N,Y) TERMS(W1,Y,Y,Y,Y,Y,Y) TERMS(W2,Y,Y,Y,Y,N,Y) TERMS(R1,Y,Y,Y,Y,Y,Y) TERMS(R2,Y,Y,Y,Y,N,Y) TERMS(RS,Y,Y,Y,Y,N,Y) TERMS(KE,N,Y,Y,Y,Y,N) TERMS(EP,N,Y,Y,Y,Y,N) TERMS(H1,N,Y,Y,Y,Y,Y) TERMS(H2,N,Y,Y,Y,N,Y) TERMS(CRAD,N,N,Y,Y,Y,N) TERMS(YH2O,N,Y,Y,Y,Y,N) TERMS(YCO2,N,Y,Y,Y,Y,N) TERMS(YCO,N,Y,Y,Y,Y,N) TERMS(YO2,N,Y,Y,Y,Y,N) TERMS(YCH4,N,Y,Y,Y,Y,N) TERMS(COL2,N,Y,Y,Y,N,N) TERMS(CHA2,N,Y,Y,Y,N,N) TERMS(YHCN,N,Y,Y,Y,Y,N) TERMS(YNOX,N,Y,Y,Y,Y,N) DIFCUT =0.5 ;ZDIFAC =1. GALA = F ;ADDDIF = F NEWRH1 = T NEWENT = T EQDVDP = F ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1 ************************************************************ Group 9. Properties used if PRPS is not stored, and where PRPS = -1.0 if it is! RHO1 = GRND5 ;TMP1 = GRND ;EL1 = GRND4 TSURR =0. ;TEMP0 =0. ;PRESS0 =1.0E+05 DVO1DT =0. ;DRH1DP = GRND5 RHO1A =0. ;RHO1B =3.496503E-03 ;RHO1C =0.714286 TMP1A =0. ;TMP1B =0. ;TMP1C =0. EMISS =0. ;SCATT =0. RADIA =0. ;RADIB =0. EL1A =0. ;EL1B =0. ;EL1C =0. RHO2 =1000. ;TMP2 = GRND ;EL2 =0. DVO2DT =0. ;DRH2DP =0. TMP2A =0. ;TMP2B =0. ;TMP2C =0. ENUL =1.0E-05 ;ENUT = GRND3 ENUTA =0. ;ENUTB =0. ;ENUTC =0. IENUTA = 0 PHINT(U1)=-2.022E+04 ;PHINT(U2)=-2.022E+04 PHINT(V1)=-2.022E+04 ;PHINT(V2)=-2.022E+04 PHINT(W1)=-2.022E+04 ;PHINT(W2)=-2.022E+04 PHINT(R1)=-2.022E+04 ;PHINT(R2)=-2.022E+04 PHINT(RS)=-2.022E+04 ;PHINT(H1)=-2.022E+04 PHINT(H2)=-2.022E+04 PRNDTL(U1)=1. ;PRNDTL(U2)=1. PRNDTL(V1)=1. ;PRNDTL(V2)=1. PRNDTL(W1)=1. ;PRNDTL(W2)=1. PRNDTL(R1)=1.0E+09 ;PRNDTL(R2)=1.0E+09 PRNDTL(RS)=1. ;PRNDTL(KE)=1. PRNDTL(EP)=1. ;PRNDTL(H1)=1. PRNDTL(H2)=1. ;PRNDTL(CRAD)=GRND PRNDTL(YH2O)=1. ;PRNDTL(YCO2)=1. PRNDTL(YCO)=1. ;PRNDTL(YO2)=1. PRNDTL(YCH4)=1. ;PRNDTL(COL2)=1. PRNDTL(CHA2)=1. ;PRNDTL(YHCN)=1. PRNDTL(YNOX)=1. PRT(U1)=1. ;PRT(U2)=1. PRT(V1)=1. ;PRT(V2)=1. PRT(W1)=1. ;PRT(W2)=1. PRT(R1)=1.0E+09 ;PRT(R2)=1.0E+09 PRT(RS)=1. ;PRT(KE)=1. PRT(EP)=1.314 ;PRT(H1)=1. PRT(H2)=1. ;PRT(CRAD)=1.0E+09 PRT(YH2O)=1. ;PRT(YCO2)=1. PRT(YCO)=1. ;PRT(YO2)=1. PRT(YCH4)=1. ;PRT(COL2)=1. PRT(CHA2)=1. ;PRT(YHCN)=1. PRT(YNOX)=1. CP1 =1050. ;CP2 =1800. ************************************************************ Group 10.Inter-Phase Transfer Processes CFIPS = GRND7 RLOLIM =1.0E-09 ;CMDOT = GRND CFIPA =1.0E-09 ;CFIPB =2.5E-05 CFIPC =100. ;CFIPD =0. CMDTA =0. ;CMDTB =0. CMDTC =0. ;CMDTD =0. CINT(P1)=1. ;CINT(U1)=1. CINT(U2)=1. ;CINT(V1)=1. CINT(V2)=1. ;CINT(W1)=1. CINT(W2)=1. ;CINT(R1)=1. CINT(R2)=1. ;CINT(RS)=1. CINT(H1)=1.0E-20 ;CINT(H2)=1.0E-20 PHINT(P1)=-2.022E+04 ;PHINT(U1)=-2.022E+04 PHINT(U2)=-2.022E+04 ;PHINT(V1)=-2.022E+04 PHINT(V2)=-2.022E+04 ;PHINT(W1)=-2.022E+04 PHINT(W2)=-2.022E+04 ;PHINT(R1)=-2.022E+04 PHINT(R2)=-2.022E+04 ;PHINT(RS)=-2.022E+04 PHINT(H1)=-2.022E+04 ;PHINT(H2)=-2.022E+04 ************************************************************ Group 11.Initial field variables (PHIs) FIINIT(P1)=1.0E-10 ;FIINIT(U1)=1.0E-10 FIINIT(U2)=1.0E-10 ;FIINIT(V1)=1.0E-10 FIINIT(V2)=1.0E-10 ;FIINIT(W1)=1.0E-10 FIINIT(W2)=1.0E-10 ;FIINIT(R1)=0.5 FIINIT(R2)=0.5 ;FIINIT(RS)=0.5 FIINIT(KE)=0.25 ;FIINIT(EP)=1.437814 FIINIT(H1)=5.25E+05 ;FIINIT(H2)=9.0E+05 FIINIT(CRAD)=0. ;FIINIT(MCCO)=1.0E-10 FIINIT(KT)=1.0E-10 ;FIINIT(XXX)=1.0E-10 FIINIT(BRN2)=1.0E-10 ;FIINIT(BRN1)=1.0E-10 FIINIT(RNOX)=1.0E-10 ;FIINIT(MOX)=1.0E-10 FIINIT(FK1)=1.0E-10 ;FIINIT(CCOK)=1.0E-10 FIINIT(CCOD)=1.0E-10 ;FIINIT(VCOD)=1.0E-10 FIINIT(BURN)=1.0E-10 ;FIINIT(WMIX)=1.0E-10 FIINIT(CCOM)=1.0E-10 ;FIINIT(VCOM)=1.0E-10 FIINIT(YN2)=1.0E-10 ;FIINIT(ASH2)=1.0E-10 FIINIT(VCOK)=1.0E-10 ;FIINIT(YH2O)=0. FIINIT(YCO2)=0. ;FIINIT(YCO)=0. FIINIT(YO2)=0.35 ;FIINIT(YCH4)=0. FIINIT(HCOF)=1.0E-10 ;FIINIT(KC)=1.0E-10 FIINIT(KD)=1.0E-10 ;FIINIT(SIZE)=1.0E-10 FIINIT(DEN1)=1.0E-10 ;FIINIT(APRJ)=1.0E-10 FIINIT(COL2)=0.8 ;FIINIT(CHA2)=0. FIINIT(REYN)=1.0E-10 ;FIINIT(CD)=1.0E-10 FIINIT(VREL)=1.0E-10 ;FIINIT(NUSS)=1.0E-10 FIINIT(TMP2)=1.0E-10 ;FIINIT(TMP1)=1.0E-10 FIINIT(CFIP)=1.0E-10 ;FIINIT(VRAT)=1.0E-10 FIINIT(CMDO)=1.0E-10 ;FIINIT(EL1)=1.0E-10 FIINIT(ENUT)=1.0E-10 ;FIINIT(EPKE)=1.0E-10 FIINIT(YHCN)=0. ;FIINIT(YNOX)=0. No PATCHes yet used for this Group INIADD = F FSWEEP = 1 NAMFI =CHAM ************************************************************ Group 12. Patchwise adjustment of terms Patches for this group are printed with those for Group 13. Their names begin either with GP12 or & ************************************************************ Group 13. Boundary & Special Sources PATCH(IN1 ,SOUTH , 4, 4, 1, 1, 2, 2, 1, 1) COVAL(IN1 ,CRAD,1. ,142. ) PATCH(IN2 ,SOUTH , 2, 2, 1, 1, 2, 2, 1, 1) COVAL(IN2 ,CRAD,1. ,142. ) PATCH(IN3 ,SOUTH , 2, 2, 1, 1, 4, 4, 1, 1) COVAL(IN3 ,CRAD,1. ,142. ) PATCH(IN4 ,SOUTH , 4, 4, 1, 1, 4, 4, 1, 1) COVAL(IN4 ,CRAD,1. ,142. ) PATCH(IN5 ,SOUTH , 4, 4, 1, 1, 6, 6, 1, 1) COVAL(IN5 ,CRAD,1. ,142. ) PATCH(IN6 ,SOUTH , 2, 2, 1, 1, 6, 6, 1, 1) COVAL(IN6 ,CRAD,1. ,142. ) PATCH(IN7 ,SOUTH , 4, 4, 1, 1, 8, 8, 1, 1) COVAL(IN7 ,CRAD,1. ,142. ) PATCH(IN8 ,SOUTH , 2, 2, 1, 1, 8, 8, 1, 1) COVAL(IN8 ,CRAD,1. ,142. ) PATCH(IN9 ,SOUTH , 4, 4, 1, 1, 10, 10, 1, 1) COVAL(IN9 ,CRAD,1. ,142. ) PATCH(IN10 ,SOUTH , 2, 2, 1, 1, 10, 10, 1, 1) COVAL(IN10 ,CRAD,1. ,142. ) PATCH(IN11 ,SOUTH , 4, 4, 1, 1, 12, 12, 1, 1) COVAL(IN11 ,CRAD,1. ,142. ) PATCH(IN12 ,SOUTH , 2, 2, 1, 1, 12, 12, 1, 1) COVAL(IN12 ,CRAD,1. ,142. ) PATCH(KESOURCE,PHASEM, 1, 5, 1, 5, 1, 20, 1, 1) COVAL(KESOURCE,KE , GRND4 , GRND4 ) COVAL(KESOURCE,EP , GRND4 , GRND4 ) PATCH(PLSR1 ,CELL , 1, 5, 1, 5, 1, 20, 1, 1) COVAL(PLSR1 ,COL2, FIXFLU , GRND ) COVAL(PLSR1 ,CHA2, FIXFLU , GRND ) PATCH(PLSR5 ,CELL , 1, 5, 1, 5, 1, 20, 1, 1) COVAL(PLSR5 ,YCH4, FIXFLU , GRND ) PATCH(PLSR6 ,CELL , 1, 5, 1, 5, 1, 20, 1, 1) COVAL(PLSR6 ,YH2O, FIXFLU , GRND ) COVAL(PLSR6 ,YCO , FIXFLU , GRND ) COVAL(PLSR6 ,YO2 , GRND ,0. ) COVAL(PLSR6 ,YCH4, GRND ,0. ) PATCH(PLSR8 ,CELL , 1, 5, 1, 5, 1, 20, 1, 1) COVAL(PLSR8 ,CHA2, FIXFLU , GRND ) PATCH(PLSR9 ,CELL , 1, 5, 1, 5, 1, 20, 1, 1) COVAL(PLSR9 ,YCO2, FIXFLU , GRND ) COVAL(PLSR9 ,YCO , FIXFLU , GRND ) COVAL(PLSR9 ,YO2 , FIXFLU , GRND ) PATCH(PLSR7 ,CELL , 1, 5, 1, 5, 1, 20, 1, 1) COVAL(PLSR7 ,YCO2, FIXFLU , GRND ) COVAL(PLSR7 ,YCO , GRND ,0. ) COVAL(PLSR7 ,YO2 , GRND ,0. ) PATCH(PLSR3 ,CELL , 1, 5, 1, 5, 1, 20, 1, 1) COVAL(PLSR3 ,H1 , FIXFLU , GRND ) PATCH(PLSR10 ,CELL , 1, 5, 1, 5, 1, 20, 1, 1) COVAL(PLSR10 ,H1 , FIXFLU , GRND ) PATCH(PLSR4 ,VOLUME, 1, 5, 1, 5, 1, 20, 1, 1) COVAL(PLSR4 ,H1 , GRND , GRND ) COVAL(PLSR4 ,H2 , GRND , GRND ) PATCH(RADSOR ,VOLUME, 1, 5, 1, 5, 1, 20, 1, 1) COVAL(RADSOR ,H1 , FIXFLU , GRND ) COVAL(RADSOR ,H2 , FIXFLU , GRND ) PATCH(RDIOSITY,VOLUME, 1, 5, 1, 5, 1, 20, 1, 1) COVAL(RDIOSITY,CRAD, GRND , GRND ) PATCH(PLSR90 ,CELL , 1, 5, 1, 5, 1, 20, 1, 1) COVAL(PLSR90 ,COL2, FIXFLU , GRND ) COVAL(PLSR90 ,CHA2, FIXFLU , GRND ) PATCH(THERMNOX,VOLUME, 1, 5, 1, 5, 1, 20, 1, 1) COVAL(THERMNOX,YNOX, FIXFLU , GRND ) PATCH(FUELHCN ,CELL , 1, 5, 1, 5, 1, 20, 1, 1) COVAL(FUELHCN ,YHCN, FIXFLU , GRND ) PATCH(FUELNOX ,PHASEM, 1, 5, 1, 5, 1, 20, 1, 1) COVAL(FUELNOX ,YNOX, GRND , GRND ) PATCH(NOXRED ,VOLUME, 1, 5, 1, 5, 1, 20, 1, 1) COVAL(NOXRED ,YNOX, GRND ,0. ) PATCH(OB1 ,HIGH , 1, 5, 1, 5, 20, 20, 1, 1) COVAL(OB1 ,U1 ,0. ,0. ) COVAL(OB1 ,U2 ,0. ,0. ) COVAL(OB1 ,V1 ,0. ,0. ) COVAL(OB1 ,V2 ,0. ,0. ) COVAL(OB1 ,W1 ,0. ,0. ) COVAL(OB1 ,W2 ,0. ,0. ) COVAL(OB1 ,R1 ,1. ,0. ) COVAL(OB1 ,R2 ,1.0E+04 ,0. ) COVAL(OB1 ,KE ,0. , SAME ) COVAL(OB1 ,EP ,0. , SAME ) COVAL(OB1 ,H1 ,0. , SAME ) COVAL(OB1 ,H2 ,0. , SAME ) PATCH(^OB2 ,SOUTH , 2, 0, 0, 0, 0, 0, 1, 1) COVAL(^OB2 ,U1 ,0. ,0. ) COVAL(^OB2 ,U2 ,0. ,0. ) COVAL(^OB2 ,V1 ,0. ,10. ) COVAL(^OB2 ,V2 ,0. ,10. ) COVAL(^OB2 ,W1 ,0. ,0. ) COVAL(^OB2 ,W2 ,0. ,0. ) COVAL(^OB2 ,R1 , FIXFLU ,6.93 ) COVAL(^OB2 ,R2 , FIXFLU ,0.7 ) COVAL(^OB2 ,KE ,0. ,0.25 ) COVAL(^OB2 ,EP ,0. ,1.437814 ) COVAL(^OB2 ,H1 ,0. ,5.25E+05 ) COVAL(^OB2 ,H2 ,0. ,9.0E+05 ) COVAL(^OB2 ,YO2 ,0. ,0.35 ) COVAL(^OB2 ,COL2,0. ,0.8 ) PATCH(^OB3 ,SOUTH , 3, 0, 0, 0, 0, 0, 1, 1) COVAL(^OB3 ,U1 ,0. ,0. ) COVAL(^OB3 ,U2 ,0. ,0. ) COVAL(^OB3 ,V1 ,0. ,10. ) COVAL(^OB3 ,V2 ,0. ,10. ) COVAL(^OB3 ,W1 ,0. ,0. ) COVAL(^OB3 ,W2 ,0. ,0. ) COVAL(^OB3 ,R1 , FIXFLU ,6.93 ) COVAL(^OB3 ,R2 , FIXFLU ,0.7 ) COVAL(^OB3 ,KE ,0. ,0.25 ) COVAL(^OB3 ,EP ,0. ,1.437814 ) COVAL(^OB3 ,H1 ,0. ,5.25E+05 ) COVAL(^OB3 ,H2 ,0. ,9.0E+05 ) COVAL(^OB3 ,YO2 ,0. ,0.35 ) COVAL(^OB3 ,COL2,0. ,0.8 ) PATCH(^OB4 ,SOUTH , 4, 0, 0, 0, 0, 0, 1, 1) COVAL(^OB4 ,U1 ,0. ,0. ) COVAL(^OB4 ,U2 ,0. ,0. ) COVAL(^OB4 ,V1 ,0. ,10. ) COVAL(^OB4 ,V2 ,0. ,10. ) COVAL(^OB4 ,W1 ,0. ,0. ) COVAL(^OB4 ,W2 ,0. ,0. ) COVAL(^OB4 ,R1 , FIXFLU ,6.93 ) COVAL(^OB4 ,R2 , FIXFLU ,0.7 ) COVAL(^OB4 ,KE ,0. ,0.25 ) COVAL(^OB4 ,EP ,0. ,1.437814 ) COVAL(^OB4 ,H1 ,0. ,5.25E+05 ) COVAL(^OB4 ,H2 ,0. ,9.0E+05 ) COVAL(^OB4 ,YO2 ,0. ,0.35 ) COVAL(^OB4 ,COL2,0. ,0.8 ) PATCH(^OB5 ,SOUTH , 5, 0, 0, 0, 0, 0, 1, 1) COVAL(^OB5 ,U1 ,0. ,0. ) COVAL(^OB5 ,U2 ,0. ,0. ) COVAL(^OB5 ,V1 ,0. ,10. ) COVAL(^OB5 ,V2 ,0. ,10. ) COVAL(^OB5 ,W1 ,0. ,0. ) COVAL(^OB5 ,W2 ,0. ,0. ) COVAL(^OB5 ,R1 , FIXFLU ,6.93 ) COVAL(^OB5 ,R2 , FIXFLU ,0.7 ) COVAL(^OB5 ,KE ,0. ,0.25 ) COVAL(^OB5 ,EP ,0. ,1.437814 ) COVAL(^OB5 ,H1 ,0. ,5.25E+05 ) COVAL(^OB5 ,H2 ,0. ,9.0E+05 ) COVAL(^OB5 ,YO2 ,0. ,0.35 ) COVAL(^OB5 ,COL2,0. ,0.8 ) PATCH(^OB6 ,SOUTH , 6, 0, 0, 0, 0, 0, 1, 1) COVAL(^OB6 ,U1 ,0. ,0. ) COVAL(^OB6 ,U2 ,0. ,0. ) COVAL(^OB6 ,V1 ,0. ,10. ) COVAL(^OB6 ,V2 ,0. ,10. ) COVAL(^OB6 ,W1 ,0. ,0. ) COVAL(^OB6 ,W2 ,0. ,0. ) COVAL(^OB6 ,R1 , FIXFLU ,6.93 ) COVAL(^OB6 ,R2 , FIXFLU ,0.7 ) COVAL(^OB6 ,KE ,0. ,0.25 ) COVAL(^OB6 ,EP ,0. ,1.437814 ) COVAL(^OB6 ,H1 ,0. ,5.25E+05 ) COVAL(^OB6 ,H2 ,0. ,9.0E+05 ) COVAL(^OB6 ,YO2 ,0. ,0.35 ) COVAL(^OB6 ,COL2,0. ,0.8 ) PATCH(^OB7 ,SOUTH , 7, 0, 0, 0, 0, 0, 1, 1) COVAL(^OB7 ,U1 ,0. ,0. ) COVAL(^OB7 ,U2 ,0. ,0. ) COVAL(^OB7 ,V1 ,0. ,10. ) COVAL(^OB7 ,V2 ,0. ,10. ) COVAL(^OB7 ,W1 ,0. ,0. ) COVAL(^OB7 ,W2 ,0. ,0. ) COVAL(^OB7 ,R1 , FIXFLU ,6.93 ) COVAL(^OB7 ,R2 , FIXFLU ,0.7 ) COVAL(^OB7 ,KE ,0. ,0.25 ) COVAL(^OB7 ,EP ,0. ,1.437814 ) COVAL(^OB7 ,H1 ,0. ,5.25E+05 ) COVAL(^OB7 ,H2 ,0. ,9.0E+05 ) COVAL(^OB7 ,YO2 ,0. ,0.35 ) COVAL(^OB7 ,COL2,0. ,0.8 ) PATCH(^OB8 ,SOUTH , 8, 0, 0, 0, 0, 0, 1, 1) COVAL(^OB8 ,U1 ,0. ,0. ) COVAL(^OB8 ,U2 ,0. ,0. ) COVAL(^OB8 ,V1 ,0. ,10. ) COVAL(^OB8 ,V2 ,0. ,10. ) COVAL(^OB8 ,W1 ,0. ,0. ) COVAL(^OB8 ,W2 ,0. ,0. ) COVAL(^OB8 ,R1 , FIXFLU ,6.93 ) COVAL(^OB8 ,R2 , FIXFLU ,0.7 ) COVAL(^OB8 ,KE ,0. ,0.25 ) COVAL(^OB8 ,EP ,0. ,1.437814 ) COVAL(^OB8 ,H1 ,0. ,5.25E+05 ) COVAL(^OB8 ,H2 ,0. ,9.0E+05 ) COVAL(^OB8 ,YO2 ,0. ,0.35 ) COVAL(^OB8 ,COL2,0. ,0.8 ) PATCH(^OB9 ,SOUTH , 9, 0, 0, 0, 0, 0, 1, 1) COVAL(^OB9 ,U1 ,0. ,0. ) COVAL(^OB9 ,U2 ,0. ,0. ) COVAL(^OB9 ,V1 ,0. ,10. ) COVAL(^OB9 ,V2 ,0. ,10. ) COVAL(^OB9 ,W1 ,0. ,0. ) COVAL(^OB9 ,W2 ,0. ,0. ) COVAL(^OB9 ,R1 , FIXFLU ,6.93 ) COVAL(^OB9 ,R2 , FIXFLU ,0.7 ) COVAL(^OB9 ,KE ,0. ,0.25 ) COVAL(^OB9 ,EP ,0. ,1.437814 ) COVAL(^OB9 ,H1 ,0. ,5.25E+05 ) COVAL(^OB9 ,H2 ,0. ,9.0E+05 ) COVAL(^OB9 ,YO2 ,0. ,0.35 ) COVAL(^OB9 ,COL2,0. ,0.8 ) PATCH(^OBA ,SOUTH , 10, 0, 0, 0, 0, 0, 1, 1) COVAL(^OBA ,U1 ,0. ,0. ) COVAL(^OBA ,U2 ,0. ,0. ) COVAL(^OBA ,V1 ,0. ,10. ) COVAL(^OBA ,V2 ,0. ,10. ) COVAL(^OBA ,W1 ,0. ,0. ) COVAL(^OBA ,W2 ,0. ,0. ) COVAL(^OBA ,R1 , FIXFLU ,6.93 ) COVAL(^OBA ,R2 , FIXFLU ,0.7 ) COVAL(^OBA ,KE ,0. ,0.25 ) COVAL(^OBA ,EP ,0. ,1.437814 ) COVAL(^OBA ,H1 ,0. ,5.25E+05 ) COVAL(^OBA ,H2 ,0. ,9.0E+05 ) COVAL(^OBA ,YO2 ,0. ,0.35 ) COVAL(^OBA ,COL2,0. ,0.8 ) PATCH(^OBB ,SOUTH , 11, 0, 0, 0, 0, 0, 1, 1) COVAL(^OBB ,U1 ,0. ,0. ) COVAL(^OBB ,U2 ,0. ,0. ) COVAL(^OBB ,V1 ,0. ,10. ) COVAL(^OBB ,V2 ,0. ,10. ) COVAL(^OBB ,W1 ,0. ,0. ) COVAL(^OBB ,W2 ,0. ,0. ) COVAL(^OBB ,R1 , FIXFLU ,6.93 ) COVAL(^OBB ,R2 , FIXFLU ,0.7 ) COVAL(^OBB ,KE ,0. ,0.25 ) COVAL(^OBB ,EP ,0. ,1.437814 ) COVAL(^OBB ,H1 ,0. ,5.25E+05 ) COVAL(^OBB ,H2 ,0. ,9.0E+05 ) COVAL(^OBB ,YO2 ,0. ,0.35 ) COVAL(^OBB ,COL2,0. ,0.8 ) PATCH(^OBC ,SOUTH , 12, 0, 0, 0, 0, 0, 1, 1) COVAL(^OBC ,U1 ,0. ,0. ) COVAL(^OBC ,U2 ,0. ,0. ) COVAL(^OBC ,V1 ,0. ,10. ) COVAL(^OBC ,V2 ,0. ,10. ) COVAL(^OBC ,W1 ,0. ,0. ) COVAL(^OBC ,W2 ,0. ,0. ) COVAL(^OBC ,R1 , FIXFLU ,6.93 ) COVAL(^OBC ,R2 , FIXFLU ,0.7 ) COVAL(^OBC ,KE ,0. ,0.25 ) COVAL(^OBC ,EP ,0. ,1.437814 ) COVAL(^OBC ,H1 ,0. ,5.25E+05 ) COVAL(^OBC ,H2 ,0. ,9.0E+05 ) COVAL(^OBC ,YO2 ,0. ,0.35 ) COVAL(^OBC ,COL2,0. ,0.8 ) PATCH(^OBD ,SOUTH , 13, 0, 0, 0, 0, 0, 1, 1) COVAL(^OBD ,U1 ,0. ,0. ) COVAL(^OBD ,U2 ,0. ,0. ) COVAL(^OBD ,V1 ,0. ,10. ) COVAL(^OBD ,V2 ,0. ,10. ) COVAL(^OBD ,W1 ,0. ,0. ) COVAL(^OBD ,W2 ,0. ,0. ) COVAL(^OBD ,R1 , FIXFLU ,6.93 ) COVAL(^OBD ,R2 , FIXFLU ,0.7 ) COVAL(^OBD ,KE ,0. ,0.25 ) COVAL(^OBD ,EP ,0. ,1.437814 ) COVAL(^OBD ,H1 ,0. ,5.25E+05 ) COVAL(^OBD ,H2 ,0. ,9.0E+05 ) COVAL(^OBD ,YO2 ,0. ,0.35 ) COVAL(^OBD ,COL2,0. ,0.8 ) XCYCLE = F EGWF = T WALLCO = GRND2 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 1250 ;ISWC1 = 1 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1 ISWR1 = 1 ;ISWR2 = 10000 SELREF = T RESFAC =1.0E-03 ************************************************************ Group 16. Terminate Iterations LITER(P1)=20 ;LITER(U1)=10 LITER(U2)=10 ;LITER(V1)=10 LITER(V2)=10 ;LITER(W1)=10 LITER(W2)=10 ;LITER(R1)=1 LITER(R2)=1 ;LITER(RS)=1 LITER(KE)=20 ;LITER(EP)=20 LITER(H1)=20 ;LITER(H2)=20 LITER(CRAD)=20 ;LITER(YH2O)=20 LITER(YCO2)=20 ;LITER(YCO)=20 LITER(YO2)=20 ;LITER(YCH4)=20 LITER(COL2)=20 ;LITER(CHA2)=20 LITER(YHCN)=20 ;LITER(YNOX)=20 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03 ENDIT(U2)=1.0E-03 ;ENDIT(V1)=1.0E-03 ENDIT(V2)=1.0E-03 ;ENDIT(W1)=1.0E-03 ENDIT(W2)=1.0E-03 ;ENDIT(R1)=1.0E-03 ENDIT(R2)=1.0E-03 ;ENDIT(RS)=1.0E-03 ENDIT(KE)=1.0E-03 ;ENDIT(EP)=1.0E-03 ENDIT(H1)=1.0E-03 ;ENDIT(H2)=1.0E-03 ENDIT(CRAD)=1.0E-03 ;ENDIT(YH2O)=1.0E-03 ENDIT(YCO2)=1.0E-03 ;ENDIT(YCO)=1.0E-03 ENDIT(YO2)=1.0E-03 ;ENDIT(YCH4)=1.0E-03 ENDIT(COL2)=1.0E-03 ;ENDIT(CHA2)=1.0E-03 ENDIT(YHCN)=1.0E-03 ;ENDIT(YNOX)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(P1,LINRLX,0.5) RELAX(U1,FALSDT,1.0E-02) RELAX(U2,FALSDT,1.0E-02) RELAX(V1,FALSDT,1.0E-02) RELAX(V2,FALSDT,1.0E-02) RELAX(W1,FALSDT,1.0E-02) RELAX(W2,FALSDT,1.0E-02) RELAX(R1,LINRLX,0.5) RELAX(R2,LINRLX,0.5) RELAX(RS,LINRLX,0.4) RELAX(KE,LINRLX,0.2) RELAX(EP,LINRLX,0.2) RELAX(H1,FALSDT,1.0E-02) RELAX(H2,FALSDT,1.0E-02) RELAX(CRAD,FALSDT,1.0E-02) RELAX(MCCO,LINRLX,1.) RELAX(KT,LINRLX,1.) RELAX(XXX,LINRLX,1.) RELAX(BRN2,LINRLX,1.) RELAX(BRN1,LINRLX,1.) RELAX(RNOX,LINRLX,1.) RELAX(MOX,LINRLX,1.) RELAX(FK1,LINRLX,1.) RELAX(CCOK,LINRLX,1.) RELAX(CCOD,LINRLX,1.) RELAX(VCOD,LINRLX,1.) RELAX(BURN,LINRLX,1.) RELAX(WMIX,LINRLX,1.) RELAX(CCOM,LINRLX,1.) RELAX(VCOM,LINRLX,1.) RELAX(YN2,LINRLX,1.) RELAX(ASH2,LINRLX,1.) RELAX(VCOK,LINRLX,1.) RELAX(YH2O,FALSDT,1.0E-02) RELAX(YCO2,FALSDT,1.0E-02) RELAX(YCO,FALSDT,1.0E-02) RELAX(YO2,FALSDT,1.0E-02) RELAX(YCH4,FALSDT,1.0E-02) RELAX(HCOF,LINRLX,1.) RELAX(KC,LINRLX,1.) RELAX(KD,LINRLX,1.) RELAX(SIZE,LINRLX,1.) RELAX(DEN1,LINRLX,0.1) RELAX(APRJ,LINRLX,1.) RELAX(COL2,FALSDT,1.0E-02) RELAX(CHA2,FALSDT,1.0E-02) RELAX(REYN,LINRLX,1.) RELAX(CD,LINRLX,1.) RELAX(VREL,LINRLX,1.) RELAX(NUSS,LINRLX,1.) RELAX(TMP2,LINRLX,1.) RELAX(TMP1,LINRLX,1.) RELAX(CFIP,LINRLX,1.) RELAX(VRAT,LINRLX,1.) RELAX(CMDO,LINRLX,0.3) RELAX(EL1,LINRLX,1.) RELAX(ENUT,LINRLX,1.) RELAX(EPKE,LINRLX,1.) RELAX(YHCN,FALSDT,1.0E-02) RELAX(YNOX,FALSDT,1.0E-02) KELIN = 3 OVRRLX =0. EXPERT = F ;NNORSL = F ************************************************************ Group 18. Limits VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10 VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06 VARMAX(U2)=1.0E+06 ;VARMIN(U2)=-1.0E+06 VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06 VARMAX(V2)=1.0E+06 ;VARMIN(V2)=-1.0E+06 VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06 VARMAX(W2)=1.0E+06 ;VARMIN(W2)=-1.0E+06 VARMAX(R1)=1. ;VARMIN(R1)=1.0E-06 VARMAX(R2)=1. ;VARMIN(R2)=1.0E-09 VARMAX(RS)=1. ;VARMIN(RS)=1.0E-06 VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10 VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10 VARMAX(H1)=1.0E+10 ;VARMIN(H1)=0. VARMAX(H2)=1.0E+10 ;VARMIN(H2)=0. VARMAX(CRAD)=1.0E+10 ;VARMIN(CRAD)=-1.0E+10 VARMAX(MCCO)=1.0E+10 ;VARMIN(MCCO)=-1.0E+10 VARMAX(KT)=1.0E+10 ;VARMIN(KT)=-1.0E+10 VARMAX(XXX)=1.0E+10 ;VARMIN(XXX)=-1.0E+10 VARMAX(BRN2)=1.0E+10 ;VARMIN(BRN2)=-1.0E+10 VARMAX(BRN1)=1.0E+10 ;VARMIN(BRN1)=-1.0E+10 VARMAX(RNOX)=1.0E+10 ;VARMIN(RNOX)=-1.0E+10 VARMAX(MOX)=1.0E+10 ;VARMIN(MOX)=-1.0E+10 VARMAX(FK1)=1.0E+10 ;VARMIN(FK1)=-1.0E+10 VARMAX(CCOK)=1.0E+10 ;VARMIN(CCOK)=-1.0E+10 VARMAX(CCOD)=1.0E+10 ;VARMIN(CCOD)=-1.0E+10 VARMAX(VCOD)=1.0E+10 ;VARMIN(VCOD)=-1.0E+10 VARMAX(BURN)=1.0E+10 ;VARMIN(BURN)=-1.0E+10 VARMAX(WMIX)=1.0E+10 ;VARMIN(WMIX)=-1.0E+10 VARMAX(CCOM)=1.0E+10 ;VARMIN(CCOM)=0. VARMAX(VCOM)=1.0E+10 ;VARMIN(VCOM)=0. VARMAX(YN2)=1. ;VARMIN(YN2)=0. VARMAX(ASH2)=1. ;VARMIN(ASH2)=0. VARMAX(VCOK)=1.0E+10 ;VARMIN(VCOK)=-1.0E+10 VARMAX(YH2O)=1. ;VARMIN(YH2O)=0. VARMAX(YCO2)=1. ;VARMIN(YCO2)=0. VARMAX(YCO)=1. ;VARMIN(YCO)=0. VARMAX(YO2)=1. ;VARMIN(YO2)=0. VARMAX(YCH4)=1. ;VARMIN(YCH4)=0. VARMAX(HCOF)=1.0E+10 ;VARMIN(HCOF)=-1.0E+10 VARMAX(KC)=1.0E+10 ;VARMIN(KC)=-1.0E+10 VARMAX(KD)=1.0E+10 ;VARMIN(KD)=-1.0E+10 VARMAX(SIZE)=1.0E+10 ;VARMIN(SIZE)=-1.0E+10 VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=-1.0E+10 VARMAX(APRJ)=1.0E+10 ;VARMIN(APRJ)=-1.0E+10 VARMAX(COL2)=1. ;VARMIN(COL2)=0. VARMAX(CHA2)=1. ;VARMIN(CHA2)=0. VARMAX(REYN)=1.0E+10 ;VARMIN(REYN)=-1.0E+10 VARMAX(CD)=1.0E+10 ;VARMIN(CD)=-1.0E+10 VARMAX(VREL)=1.0E+10 ;VARMIN(VREL)=-1.0E+10 VARMAX(NUSS)=1.0E+10 ;VARMIN(NUSS)=-1.0E+10 VARMAX(TMP2)=1.0E+10 ;VARMIN(TMP2)=300. VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=300. VARMAX(CFIP)=1.0E+10 ;VARMIN(CFIP)=-1.0E+10 VARMAX(VRAT)=1.0E+10 ;VARMIN(VRAT)=-1.0E+10 VARMAX(CMDO)=1.0E+10 ;VARMIN(CMDO)=0. VARMAX(EL1)=1.0E+10 ;VARMIN(EL1)=-1.0E+10 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10 VARMAX(YHCN)=1. ;VARMIN(YHCN)=0. VARMAX(YNOX)=1. ;VARMIN(YNOX)=0. ************************************************************ Group 19. Data transmitted to GROUND NAMSAT =MOSG GENK = T PARSOL = F ISG62 = 1 SPEDAT(SET,GXMONI,PLOTALL,L,T) SPEDAT(SET,OBJNAM,!OB1,C,OUTLET) SPEDAT(SET,OBJTYP,!OB1,C,OUTLET) SPEDAT(SET,ARATIO,!OB1,R,1.) SPEDAT(SET,OBJNAM,^OB2,C,B4) SPEDAT(SET,OBJTYP,^OB2,C,INLET) SPEDAT(SET,OBJNAM,^OB3,C,B6) SPEDAT(SET,OBJTYP,^OB3,C,INLET) SPEDAT(SET,OBJNAM,^OB4,C,B7) SPEDAT(SET,OBJTYP,^OB4,C,INLET) SPEDAT(SET,OBJNAM,^OB5,C,B9) SPEDAT(SET,OBJTYP,^OB5,C,INLET) SPEDAT(SET,OBJNAM,^OB6,C,B8) SPEDAT(SET,OBJTYP,^OB6,C,INLET) SPEDAT(SET,OBJNAM,^OB7,C,B10) SPEDAT(SET,OBJTYP,^OB7,C,INLET) SPEDAT(SET,OBJNAM,^OB8,C,B11) SPEDAT(SET,OBJTYP,^OB8,C,INLET) SPEDAT(SET,OBJNAM,^OB9,C,B12) SPEDAT(SET,OBJTYP,^OB9,C,INLET) SPEDAT(SET,OBJNAM,^OBA,C,B13) SPEDAT(SET,OBJTYP,^OBA,C,INLET) SPEDAT(SET,OBJNAM,^OBB,C,B14) SPEDAT(SET,OBJTYP,^OBB,C,INLET) SPEDAT(SET,OBJNAM,^OBC,C,B15) SPEDAT(SET,OBJTYP,^OBC,C,INLET) SPEDAT(SET,OBJNAM,^OBD,C,B16) SPEDAT(SET,OBJTYP,^OBD,C,INLET) SPEDAT(SET,FACETDAT,NUMOBJ,I,13) ************************************************************ Group 20. Preliminary Printout ************************************************************ Group 21. Print-out of Variables INIFLD = F ;SUBWGR = F * Y in OUTPUT argument list denotes: * 1-field 2-correction-eq. monitor 3-selective dumping * 4-whole-field residual 5-spot-value table 6-residual table OUTPUT(P1,Y,Y,Y,Y,Y,Y) OUTPUT(U1,Y,Y,Y,N,Y,Y) OUTPUT(U2,Y,Y,Y,N,Y,Y) OUTPUT(V1,Y,Y,Y,N,Y,Y) OUTPUT(V2,Y,Y,Y,N,Y,Y) OUTPUT(W1,Y,Y,Y,Y,Y,Y) OUTPUT(W2,Y,Y,Y,Y,Y,Y) OUTPUT(R1,Y,Y,Y,Y,Y,Y) OUTPUT(R2,Y,Y,Y,Y,Y,Y) OUTPUT(RS,Y,N,Y,Y,Y,Y) OUTPUT(KE,Y,N,Y,Y,Y,Y) OUTPUT(EP,Y,N,Y,Y,Y,Y) OUTPUT(H1,Y,N,Y,Y,Y,Y) OUTPUT(H2,Y,N,Y,Y,Y,Y) OUTPUT(CRAD,Y,N,Y,Y,Y,Y) OUTPUT(MCCO,Y,N,Y,N,N,N) OUTPUT(KT,Y,N,Y,N,N,N) OUTPUT(XXX,Y,N,Y,N,N,N) OUTPUT(BRN2,Y,N,Y,N,N,N) OUTPUT(BRN1,Y,N,Y,N,N,N) OUTPUT(RNOX,Y,N,Y,N,N,N) OUTPUT(MOX,Y,N,Y,N,N,N) OUTPUT(FK1,Y,N,Y,N,N,N) OUTPUT(CCOK,Y,N,Y,N,N,N) OUTPUT(CCOD,Y,N,Y,N,N,N) OUTPUT(VCOD,Y,N,Y,N,N,N) OUTPUT(BURN,Y,N,Y,N,N,N) OUTPUT(WMIX,Y,N,Y,N,N,N) OUTPUT(CCOM,Y,N,Y,N,N,N) OUTPUT(VCOM,Y,N,Y,N,N,N) OUTPUT(YN2,Y,N,Y,N,N,N) OUTPUT(ASH2,Y,N,Y,N,N,N) OUTPUT(VCOK,Y,N,Y,N,N,N) OUTPUT(YH2O,Y,N,Y,Y,Y,Y) OUTPUT(YCO2,Y,N,Y,Y,Y,Y) OUTPUT(YCO,Y,N,Y,Y,Y,Y) OUTPUT(YO2,Y,N,Y,Y,Y,Y) OUTPUT(YCH4,Y,N,Y,Y,Y,Y) OUTPUT(HCOF,Y,N,Y,N,N,N) OUTPUT(KC,Y,N,Y,N,N,N) OUTPUT(KD,Y,N,Y,N,N,N) OUTPUT(SIZE,Y,N,Y,N,N,N) OUTPUT(DEN1,Y,N,Y,N,N,N) OUTPUT(APRJ,Y,N,Y,N,N,N) OUTPUT(COL2,Y,N,Y,Y,Y,Y) OUTPUT(CHA2,Y,N,Y,Y,Y,Y) OUTPUT(REYN,Y,N,Y,N,N,N) OUTPUT(CD,Y,N,Y,N,N,N) OUTPUT(VREL,Y,N,Y,N,N,N) OUTPUT(NUSS,Y,N,Y,N,N,N) OUTPUT(TMP2,Y,N,Y,N,N,N) OUTPUT(TMP1,Y,N,Y,N,N,N) OUTPUT(CFIP,Y,N,Y,N,N,N) OUTPUT(VRAT,Y,N,Y,N,N,N) OUTPUT(CMDO,Y,N,Y,N,N,N) OUTPUT(EL1,Y,N,Y,N,N,N) OUTPUT(ENUT,Y,N,Y,N,N,N) OUTPUT(EPKE,Y,N,Y,N,N,N) OUTPUT(YHCN,Y,N,Y,Y,Y,Y) OUTPUT(YNOX,Y,N,Y,Y,Y,Y) ************************************************************ Group 22. Monitor Print-Out IXMON = 3 ;IYMON = 4 ;IZMON = 9 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1 UWATCH = T ;USTEER = T HIGHLO = F ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 ;NUMCLS = 5 NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000 NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000 NZPRIN = 1 ;IZPRF = 1 ;IZPRL = 10000 XZPR = F ;YZPR = F IPLTF = 1 ;IPLTL = -1 ;NPLT = -1 ISWPRF = 1 ;ISWPRL = 100000 ITABL = 3 ;IPROF = 1 ABSIZ =0.5 ;ORSIZ =0.4 NTZPRF = 1 ;NCOLPF = 50 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20 No PATCHes yet used for this Group ************************************************************ Group 24. Dumps For Restarts SAVE = T ;NOWIPE = F NSAVE =CHAM STOP