TALK=T;RUN(1,1) PHOTON USE p phi;;;;; set con fi dep 3 con yo2 x 1 fi;0.05 gr ou x 1 gr ou x 1 y 25 32 z 1 5 te 1 oxygen mass-fraction contours 0.13265E+04 0.19215E+04 CR msg Pressto continue pause cl con fue x 1 fi;0.05 gr ou x 1 gr ou x 1 y 25 32 z 1 5 te 1 contours of volume fraction of solid 0.13265E+04 0.21215E+04 CR msg Press to continue pause cl con yco2 x 1 fi;0.05 gr ou x 1 gr ou x 1 y 25 32 z 1 5 te 1 carbon dioxide contours 0.13265E+04 0.21215E+04 CR msg Press to continue pause cl con tmp1 x 1 fi;0.05 gr ou x 1 gr ou x 1 y 25 32 z 1 5 te 1 gas temperature contours 0.13265E+04 0.23215E+04 CR msg Press to continue pause cl con tmp2 x 1 fi;0.05 gr ou x 1 gr ou x 1 y 25 32 z 1 5 te 1 coal temperature contours 0.13265E+04 0.23215E+04 CR msg Press to continue pause cl con mixf x 1 fi;0.05 gr ou x 1 gr ou x 1 y 25 32 z 1 5 te 1 mixture fraction contours 0.13265E+04 0.23215E+04 CR msg Press to continue pause enduse DISPLAY The case considered is the 2d turbulent flow of combusting pulverised coal particles in the near-field region of a power-station burner. The burner comprises 3 air streams, tertiary air enters through a central air stream, swirling primary air carrying coal particles enters through an inner annulus located above the tertiary-air stream, and swirling secondary air enters through an outer annulus located above the primary air stream. The main combustion process is simulated using an IPSA-based 'equilibrium' two-phase combustion model. A 2nd post-processing run can be performed which computes thermal NOx based on CHO equilibrium as determined by the CREK module. =================================================== Coal-combustion model; solid phase - carbon only reactions: C (s) + 0.5 O2 > CO (exothermic ) CO + 0.5 O2 > CO2 (exothermic ) C(s) + CO2 > 2CO (endothermic) C(s) + H2O > CO + H2 (endothermic) H2 + 0.5 O2 > H2O (exothermic ) =================================================== ENDDIS noxdbegin ns 7 nr 3 noxdend ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(2-D Coal Combustion With NOX ) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 14 ************************************************************ Group 2. Time dependence STEADY = T ************************************************************ Group 3. X-Direction Grid Spacing CARTES = F NX = 1 XULAST =1.0E-02 XFRAC(1)=1. ************************************************************ Group 4. Y-Direction Grid Spacing NY = 32 YVLAST =1. YFRAC(1)=0.03993 ;YFRAC(2)=0.078285 YFRAC(3)=0.114752 ;YFRAC(4)=0.148828 YFRAC(5)=0.179533 ;YFRAC(6)=0.2032 YFRAC(7)=0.219075 ;YFRAC(8)=0.23495 YFRAC(9)=0.250825 ;YFRAC(10)=0.2667 YFRAC(11)=0.282575 ;YFRAC(12)=0.29845 YFRAC(13)=0.314325 ;YFRAC(14)=0.3302 YFRAC(15)=0.351836 ;YFRAC(16)=0.379905 YFRAC(17)=0.411056 ;YFRAC(18)=0.444393 YFRAC(19)=0.479456 ;YFRAC(20)=0.515959 YFRAC(21)=0.553704 ;YFRAC(22)=0.592547 YFRAC(23)=0.632375 ;YFRAC(24)=0.6731 YFRAC(25)=0.700059 ;YFRAC(26)=0.735036 YFRAC(27)=0.773852 ;YFRAC(28)=0.815391 YFRAC(29)=0.859082 ;YFRAC(30)=0.904567 YFRAC(31)=0.9516 ;YFRAC(32)=1. ************************************************************ Group 5. Z-Direction Grid Spacing PARAB = F NZ = 25 ZWLAST =5.28735 ZFRAC(1)=0.010869 ;ZFRAC(2)=0.021739 ZFRAC(3)=0.032608 ;ZFRAC(4)=0.043477 ZFRAC(5)=0.054347 ;ZFRAC(6)=0.060154 ZFRAC(7)=0.073215 ;ZFRAC(8)=0.091938 ZFRAC(9)=0.11565 ;ZFRAC(10)=0.143931 ZFRAC(11)=0.176481 ;ZFRAC(12)=0.213073 ZFRAC(13)=0.253521 ;ZFRAC(14)=0.297675 ZFRAC(15)=0.345406 ;ZFRAC(16)=0.396601 ZFRAC(17)=0.451163 ;ZFRAC(18)=0.509004 ZFRAC(19)=0.570048 ;ZFRAC(20)=0.634224 ZFRAC(21)=0.701467 ;ZFRAC(22)=0.771718 ZFRAC(23)=0.844924 ;ZFRAC(24)=0.921033 ZFRAC(25)=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)=GAS NAME(10)=FUE ;NAME(11)=SHAD NAME(12)=KE ;NAME(13)=EP NAME(14)=H1 ;NAME(15)=H2 NAME(16)=MIXF ;NAME(17)=C2 NAME(133)=PRPS ;NAME(134)=NPOR NAME(136)=EPKE ;NAME(137)=RADY NAME(138)=RADZ ;NAME(139)=ENUT NAME(140)=RHO1 ;NAME(141)=TMP2 NAME(142)=TMP1 ;NAME(143)=YH2O NAME(144)=YH2 ;NAME(145)=YN2 NAME(146)=YCO2 ;NAME(147)=YCO NAME(148)=YO2 ;NAME(149)=CFIP NAME(150)=MDOT * 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(GAS,Y,Y,N,N,N,Y) SOLUTN(FUE,Y,Y,N,N,N,Y) SOLUTN(SHAD,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,Y,N,N,Y) SOLUTN(H2,Y,Y,Y,N,N,Y) SOLUTN(MIXF,Y,Y,N,N,N,Y) SOLUTN(C2,Y,Y,N,N,N,Y) SOLUTN(PRPS,Y,N,N,N,N,N) SOLUTN(NPOR,Y,N,N,N,N,N) SOLUTN(EPKE,Y,N,N,N,N,N) SOLUTN(RADY,Y,Y,N,N,N,N) SOLUTN(RADZ,Y,Y,Y,N,N,N) SOLUTN(ENUT,Y,N,N,N,N,Y) SOLUTN(RHO1,Y,N,N,N,N,Y) SOLUTN(TMP2,Y,N,N,N,N,Y) SOLUTN(TMP1,Y,N,N,N,N,Y) SOLUTN(YH2O,Y,N,N,N,N,Y) SOLUTN(YH2,Y,N,N,N,N,Y) SOLUTN(YN2,Y,N,N,N,N,Y) SOLUTN(YCO2,Y,N,N,N,N,Y) SOLUTN(YCO,Y,N,N,N,N,Y) SOLUTN(YO2,Y,N,N,N,N,Y) SOLUTN(CFIP,Y,N,N,N,N,Y) SOLUTN(MDOT,Y,N,N,N,N,Y) DEN1 = 140 VIST = 139 INTMDT = 150 EPOR = 0 ;HPOR = 0 ;NPOR = 134 ;VPOR = 0 TEMP1 = 142 TEMP2 = 141 PRPS = 133 ************************************************************ 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(GAS,Y,Y,Y,Y,Y,Y) TERMS(FUE,Y,Y,Y,Y,N,Y) TERMS(SHAD,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(MIXF,N,Y,Y,Y,Y,Y) TERMS(C2,N,Y,Y,Y,N,Y) TERMS(RADY,N,N,Y,N,Y,N) TERMS(RADZ,N,N,Y,N,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 = GRND7 ;TMP1 =0. ;EL1 = GRND4 TSURR =0. ;TEMP0 =0. ;PRESS0 =1.0E+05 DVO1DT =0. ;DRH1DP =0. RHO1A =0.95 ;RHO1B =0.05 ;RHO1C =0. EMISS =0. ;SCATT =0. RADIA =0.5 ;RADIB =0.02 EL1A =0. ;EL1B =0. ;EL1C =0. RHO2 =1350. ;TMP2 =0. ;EL2 =0. DVO2DT =0. ;DRH2DP =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(GAS)=-2.022E+04 ;PHINT(FUE)=-2.022E+04 PHINT(SHAD)=-2.022E+04 ;PHINT(H1)=GRND7 PHINT(H2)=GRND7 ;PHINT(MIXF)=1. PHINT(C2)=-2.022E+04 PHNH1A =0. ;PHNH1B =0. ;PHNH1C =0. PHNH2A =0. ;PHNH2B =0. ;PHNH2C =0. PRNDTL(U1)=1. ;PRNDTL(U2)=1. PRNDTL(V1)=1. ;PRNDTL(V2)=1. PRNDTL(W1)=1. ;PRNDTL(W2)=1. PRNDTL(GAS)=1. ;PRNDTL(FUE)=1. PRNDTL(SHAD)=1. ;PRNDTL(KE)=1. PRNDTL(EP)=1. ;PRNDTL(H1)=1. PRNDTL(H2)=1. ;PRNDTL(MIXF)=1. PRNDTL(C2)=1. ;PRNDTL(RADY)=-GRND1 PRNDTL(RADZ)=-GRND1 PRT(U1)=1. ;PRT(U2)=1. PRT(V1)=1. ;PRT(V2)=1. PRT(W1)=1. ;PRT(W2)=1. PRT(GAS)=1. ;PRT(FUE)=1. PRT(SHAD)=1. ;PRT(KE)=1. PRT(EP)=1.314 ;PRT(H1)=1. PRT(H2)=1. ;PRT(MIXF)=1. PRT(C2)=1. ;PRT(RADY)=1.0E+10 PRT(RADZ)=1.0E+10 CP1 =1100. ;CP2 =1100. ************************************************************ Group 10.Inter-Phase Transfer Processes CFIPS = GRND1 RLOLIM =1.0E-06 ;CMDOT = GRND3 CFIPA =1.0E-03 ;CFIPB =0. CFIPC =1.0E+05 ;CFIPD =0. CMDTA =1.0E+04 ;CMDTB =0. CMDTC =0.154804 ;CMDTD =0. CINT(P1)=1. ;CINT(U1)=1. CINT(U2)=1. ;CINT(V1)=1. CINT(V2)=1. ;CINT(W1)=1. CINT(W2)=1. ;CINT(GAS)=1. CINT(FUE)=1. ;CINT(SHAD)=1. CINT(H1)=1. ;CINT(H2)=1. CINT(MIXF)=0. ;CINT(C2)=0. 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(GAS)=-2.022E+04 PHINT(FUE)=-2.022E+04 ;PHINT(SHAD)=-2.022E+04 PHINT(H1)=GRND7 ;PHINT(H2)=GRND7 PHINT(MIXF)=1. ;PHINT(C2)=-2.022E+04 ************************************************************ Group 11.Initial field variables (PHIs) FIINIT(P1)=1.0E-10 ;FIINIT(U1)=10. FIINIT(U2)=10. ;FIINIT(V1)=0. FIINIT(V2)=0. ;FIINIT(W1)=28.07303 FIINIT(W2)=21.354288 ;FIINIT(GAS)=0.99999 FIINIT(FUE)=1.0E-05 ;FIINIT(SHAD)=1.0E-05 FIINIT(KE)=1.970238 ;FIINIT(EP)=13.250978 FIINIT(H1)=5.863E+05 ;FIINIT(H2)=3.756981E+07 FIINIT(MIXF)=0. ;FIINIT(C2)=1. FIINIT(PRPS)=-1. ;FIINIT(NPOR)=1. FIINIT(EPKE)=1.0E-10 ;FIINIT(RADY)=1.0E-10 FIINIT(RADZ)=1.0E-10 ;FIINIT(ENUT)=1.0E-10 FIINIT(RHO1)=0.652786 ;FIINIT(TMP2)=338.555542 FIINIT(TMP1)=1.0E-10 ;FIINIT(YH2O)=1.0E-10 FIINIT(YH2)=1.0E-10 ;FIINIT(YN2)=1.0E-10 FIINIT(YCO2)=1.0E-10 ;FIINIT(YCO)=1.0E-10 FIINIT(YO2)=1.0E-10 ;FIINIT(CFIP)=1.0E-10 FIINIT(MDOT)=4.309906E-05 PATCH(INIT1 ,INIVAL, 1, 1, 1, 6, 1, 25, 1, 1) INIT(INIT1 ,U1 ,0. ,0. ) INIT(INIT1 ,U2 ,0. ,0. ) INIT(INIT1 ,W1 ,0. ,41.934875 ) INIT(INIT1 ,KE ,0. ,4.396334 ) INIT(INIT1 ,EP ,0. ,74.533218 ) INIT(INIT1 ,H1 ,0. ,5.863E+05 ) INIT(INIT1 ,RHO1,0. ,0.652786 ) PATCH(INIT2 ,INIVAL, 1, 1, 7, 14, 1, 25, 1, 1) INIT(INIT2 ,U1 ,0. ,10. ) INIT(INIT2 ,U2 ,0. ,10. ) INIT(INIT2 ,W1 ,0. ,21.355698 ) INIT(INIT2 ,KE ,0. ,1.140164 ) INIT(INIT2 ,EP ,0. ,15.750166 ) INIT(INIT2 ,H1 ,0. ,3.724111E+05 ) INIT(INIT2 ,RHO1,0. ,1.027704 ) PATCH(BLOCK1 ,INIVAL, 1, 1, 25, 32, 1, 5, 1, 1) INIT(BLOCK1 ,PRPS,0. ,199. ) PATCH(WALL1 ,INIVAL, 1, 1, 24, 24, 1, 5, 1, 1) INIT(WALL1 ,V1 ,0. ,0. ) INIT(WALL1 ,V2 ,0. ,0. ) INIT(WALL1 ,NPOR,0. ,0. ) INIADD = F FSWEEP = 1 IURINI = -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(RADISO ,FREEVL, 0, 0, 0, 0, 0, 0, 1, 1) COVAL(RADISO ,H1 , GRND1 , GRND1 ) COVAL(RADISO ,RADY, GRND1 , GRND1 ) COVAL(RADISO ,RADZ, GRND1 , GRND1 ) PATCH(KESOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1) COVAL(KESOURCE,KE , GRND4 , GRND4 ) COVAL(KESOURCE,EP , GRND4 , GRND4 ) PATCH(WALL1-NW,NWALL , 1, 1, 24, 24, 1, 5, 1, 1) COVAL(WALL1-NW,U1 , GRND2 ,0. ) COVAL(WALL1-NW,W1 , GRND2 ,0. ) COVAL(WALL1-NW,KE , GRND2 , GRND2 ) COVAL(WALL1-NW,EP , GRND2 , GRND2 ) PATCH(WALL1-SW,SWALL , 1, 1, 25, 25, 1, 5, 1, 1) COVAL(WALL1-SW,U1 , GRND2 ,0. ) COVAL(WALL1-SW,W1 , GRND2 ,0. ) COVAL(WALL1-SW,KE , GRND2 , GRND2 ) COVAL(WALL1-SW,EP , GRND2 , GRND2 ) PATCH(INPRMA ,LOW , 1, 1, 7, 14, 1, 1, 1, 1) COVAL(INPRMA ,U1 ,0. ,10. ) COVAL(INPRMA ,V1 ,0. ,0. ) COVAL(INPRMA ,W1 ,0. ,21.355698 ) COVAL(INPRMA ,GAS , FIXFLU ,22.33717 ) COVAL(INPRMA ,KE ,0. ,1.140164 ) COVAL(INPRMA ,EP ,0. ,15.750166 ) COVAL(INPRMA ,H1 ,0. ,3.724111E+05 ) COVAL(INPRMA ,MIXF,0. ,0. ) PATCH(INPRMC ,LOW , 1, 1, 7, 14, 1, 1, 1, 1) COVAL(INPRMC ,U2 ,0. ,10. ) COVAL(INPRMC ,V2 ,0. ,0. ) COVAL(INPRMC ,W2 ,0. ,21.355698 ) COVAL(INPRMC ,FUE , FIXFLU ,12.724513 ) COVAL(INPRMC ,H2 ,0. ,3.756981E+07 ) COVAL(INPRMC ,C2 ,0. ,1. ) PATCH(INSEC ,LOW , 1, 1, 15, 24, 1, 1, 1, 1) COVAL(INSEC ,U1 ,0. ,10. ) COVAL(INSEC ,V1 ,0. ,0. ) COVAL(INSEC ,W1 ,0. ,28.07303 ) COVAL(INSEC ,GAS , FIXFLU ,18.663469 ) COVAL(INSEC ,KE ,0. ,1.970238 ) COVAL(INSEC ,EP ,0. ,13.250978 ) COVAL(INSEC ,H1 ,0. ,5.863E+05 ) COVAL(INSEC ,MIXF,0. ,0. ) PATCH(INTER ,LOW , 1, 1, 1, 6, 1, 1, 1, 1) COVAL(INTER ,U1 ,0. ,0. ) COVAL(INTER ,V1 ,0. ,0. ) COVAL(INTER ,W1 ,0. ,41.934875 ) COVAL(INTER ,GAS , FIXFLU ,27.869837 ) COVAL(INTER ,KE ,0. ,4.396334 ) COVAL(INTER ,EP ,0. ,74.533218 ) COVAL(INTER ,H1 ,0. ,5.863E+05 ) COVAL(INTER ,MIXF,0. ,0. ) PATCH(OUTLET ,HIGH , 1, 1, 1, 32, 25, 25, 1, 1) COVAL(OUTLET ,U1 ,0. ,0. ) COVAL(OUTLET ,U2 ,0. ,0. ) COVAL(OUTLET ,V1 ,0. ,0. ) COVAL(OUTLET ,V2 ,0. ,0. ) COVAL(OUTLET ,W1 ,0. ,0. ) COVAL(OUTLET ,W2 ,0. ,0. ) COVAL(OUTLET ,GAS ,100. ,0. ) COVAL(OUTLET ,FUE ,1.35E+06 ,0. ) COVAL(OUTLET ,KE ,0. , SAME ) COVAL(OUTLET ,EP ,0. , SAME ) COVAL(OUTLET ,H1 ,0. , SAME ) COVAL(OUTLET ,H2 ,0. , SAME ) COVAL(OUTLET ,MIXF,0. , SAME ) COVAL(OUTLET ,C2 ,0. ,0. ) COVAL(OUTLET ,RADY,0. , SAME ) COVAL(OUTLET ,RADZ,0. , SAME ) PATCH(WALNR ,NORTH , 1, 1, 24, 24, 1, 5, 1, 1) COVAL(WALNR ,RADY,1. ,1.21226E+04 ) PATCH(WALLR ,LOW , 1, 1, 25, 32, 6, 6, 1, 1) COVAL(WALLR ,RADZ,1. ,1.21226E+04 ) IURVAL = -1 XCYCLE = F EGWF = T WALLCO = GRND2 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 350 ;ISWC1 = 1 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1 ISWR1 = 1 ;ISWR2 = 10000 SELREF = T RESFAC =1.0E-05 ************************************************************ 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(GAS)=1 LITER(FUE)=1 ;LITER(SHAD)=1 LITER(KE)=20 ;LITER(EP)=20 LITER(H1)=20 ;LITER(H2)=20 LITER(MIXF)=20 ;LITER(C2)=20 LITER(RADY)=2 ;LITER(RADZ)=2 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(GAS)=1.0E-03 ENDIT(FUE)=1.0E-03 ;ENDIT(SHAD)=1.0E-03 ENDIT(KE)=1.0E-03 ;ENDIT(EP)=1.0E-03 ENDIT(H1)=1.0E-03 ;ENDIT(H2)=1.0E-03 ENDIT(MIXF)=1.0E-03 ;ENDIT(C2)=1.0E-03 ENDIT(RADY)=1.0E-03 ;ENDIT(RADZ)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(P1,LINRLX,1.) RELAX(U1,FALSDT,1.) RELAX(U2,FALSDT,1.) RELAX(V1,FALSDT,1.) RELAX(V2,FALSDT,1.) RELAX(W1,FALSDT,1.) RELAX(W2,FALSDT,1.) RELAX(GAS,LINRLX,0.4) RELAX(FUE,LINRLX,0.4) RELAX(SHAD,LINRLX,0.4) RELAX(KE,LINRLX,0.5) RELAX(EP,LINRLX,0.5) RELAX(H1,FALSDT,1.0E+09) RELAX(H2,FALSDT,1.0E+09) RELAX(MIXF,FALSDT,1.0E+09) RELAX(C2,FALSDT,1.0E+09) RELAX(PRPS,LINRLX,1.) RELAX(EPKE,LINRLX,1.) RELAX(RADY,FALSDT,1.0E+09) RELAX(RADZ,FALSDT,1.0E+09) RELAX(ENUT,LINRLX,1.) RELAX(RHO1,LINRLX,1.) RELAX(TMP2,LINRLX,1.) RELAX(TMP1,LINRLX,1.) RELAX(YH2O,LINRLX,1.) RELAX(YH2,LINRLX,1.) RELAX(YN2,LINRLX,1.) RELAX(YCO2,LINRLX,1.) RELAX(YCO,LINRLX,1.) RELAX(YO2,LINRLX,1.) RELAX(CFIP,LINRLX,1.) RELAX(MDOT,LINRLX,1.) 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(GAS)=1. ;VARMIN(GAS)=1.0E-06 VARMAX(FUE)=1. ;VARMIN(FUE)=1.0E-12 VARMAX(SHAD)=1. ;VARMIN(SHAD)=1.0E-12 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)=-1.0E+10 VARMAX(H2)=1.0E+10 ;VARMIN(H2)=-1.0E+10 VARMAX(MIXF)=0.154804 ;VARMIN(MIXF)=0. VARMAX(C2)=1. ;VARMIN(C2)=0. VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10 VARMAX(NPOR)=1.0E+10 ;VARMIN(NPOR)=-1.0E+10 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10 VARMAX(RADY)=1.0E+10 ;VARMIN(RADY)=-1.0E+10 VARMAX(RADZ)=1.0E+10 ;VARMIN(RADZ)=-1.0E+10 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10 VARMAX(RHO1)=1.0E+10 ;VARMIN(RHO1)=0.069686 VARMAX(TMP2)=1.0E+10 ;VARMIN(TMP2)=-1.0E+10 VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=-1.0E+10 VARMAX(YH2O)=1.0E+10 ;VARMIN(YH2O)=-1.0E+10 VARMAX(YH2)=1.0E+10 ;VARMIN(YH2)=-1.0E+10 VARMAX(YN2)=1.0E+10 ;VARMIN(YN2)=-1.0E+10 VARMAX(YCO2)=1.0E+10 ;VARMIN(YCO2)=-1.0E+10 VARMAX(YCO)=1.0E+10 ;VARMIN(YCO)=-1.0E+10 VARMAX(YO2)=1.0E+10 ;VARMIN(YO2)=-1.0E+10 VARMAX(CFIP)=1.0E+10 ;VARMIN(CFIP)=-1.0E+10 VARMAX(MDOT)=1.0E+10 ;VARMIN(MDOT)=-1.0E+10 ************************************************************ Group 19. Data transmitted to GROUND GENK = T PARSOL = F CONWIZ = T ISG62 = 1 SPEDAT(SET,LONGNAME,epke,C,NGNAME,epke,c,micro-mixing_rate) SPEDAT(SET,LONGNAME,rady,C,NGNAME,rady,c,y-direction_flux_sum) SPEDAT(SET,LONGNAME,radz,C,NGNAME,radz,c,z-direction_flux_sum) SPEDAT(SET,GXMONI,PLOTALL,L,T) SPEDAT(SET,MATERIAL,199,L,T) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=57.52 ;EX(U1)=2.407 EX(U2)=2.496 ;EX(V1)=1.866 EX(V2)=1.714 ;EX(W1)=35.029999 EX(W2)=29.549999 ;EX(GAS)=0.95 EX(FUE)=0. ;EX(SHAD)=0. EX(KE)=10.85 ;EX(EP)=231.100006 EX(H1)=1.147E+06 ;EX(H2)=3.604E+07 EX(MIXF)=0.02047 ;EX(C2)=0. EX(PRPS)=0. ;EX(NPOR)=0.9913 EX(EPKE)=0. ;EX(RADY)=1.106E+05 EX(RADZ)=2.517E+05 ;EX(ENUT)=0.05748 EX(RHO1)=0.3856 ;EX(TMP2)=635.400024 EX(TMP1)=1042. ;EX(YH2O)=0. EX(YH2)=0. ;EX(YN2)=0.7139 EX(YCO2)=0.0713 ;EX(YCO)=0. EX(YO2)=0.1556 ;EX(CFIP)=0. EX(MDOT)=0. ************************************************************ 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,N,Y,Y,Y,Y) OUTPUT(U1,Y,N,Y,Y,Y,Y) OUTPUT(U2,Y,N,Y,Y,Y,Y) OUTPUT(V1,Y,N,Y,Y,Y,Y) OUTPUT(V2,Y,N,Y,Y,Y,Y) OUTPUT(W1,Y,N,Y,Y,Y,Y) OUTPUT(W2,Y,N,Y,Y,Y,Y) OUTPUT(GAS,Y,N,Y,Y,Y,Y) OUTPUT(FUE,Y,N,Y,Y,Y,Y) OUTPUT(SHAD,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(MIXF,Y,N,Y,Y,Y,Y) OUTPUT(C2,N,N,N,N,N,N) OUTPUT(PRPS,Y,N,Y,N,N,N) OUTPUT(NPOR,Y,N,Y,N,N,N) OUTPUT(EPKE,Y,N,Y,N,N,N) OUTPUT(RADY,Y,N,Y,Y,Y,Y) OUTPUT(RADZ,Y,N,Y,Y,Y,Y) OUTPUT(ENUT,Y,N,Y,N,N,N) OUTPUT(RHO1,Y,Y,Y,Y,Y,Y) OUTPUT(TMP2,Y,Y,Y,Y,Y,Y) OUTPUT(TMP1,Y,Y,Y,Y,Y,Y) OUTPUT(YH2O,Y,N,Y,N,N,N) OUTPUT(YH2,Y,N,Y,N,N,N) OUTPUT(YN2,Y,N,Y,N,N,N) OUTPUT(YCO2,Y,N,Y,N,N,N) OUTPUT(YCO,Y,N,Y,N,N,N) OUTPUT(YO2,Y,N,Y,N,N,N) OUTPUT(CFIP,Y,N,Y,N,N,N) OUTPUT(MDOT,Y,N,Y,N,N,N) ************************************************************ Group 22. Monitor Print-Out IXMON = 1 ;IYMON = 16 ;IZMON = 18 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1 UWATCH = T ;USTEER = T HIGHLO = F ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 ;NUMCLS = 5 NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000 NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000 XZPR = F ;YZPR = F IPLTF = 1 ;IPLTL = -1 ;NPLT = 10 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