TALK=T;RUN(1,1) PHOTON USE p PHI 1 1 2 use patgeo set con fi dep 3 con mixf x 10 fi;0.001 gr ou x 10 te 1 mixture fraction contours at ix=10 plane 0.50627E+03 0.31916E+03 CR msg Pressto continue pause text cl;con cl;redr con yco2 x 10 fi;0.001 gr ou x 10 te 1 CO2 contours at ix=10 plane 0.50627E+03 0.31916E+03 CR msg Press to continue pause text cl;con cl;redr con yco x 10 fi;0.001 gr ou x 10 te 1 CO contours at ix=10 plane 0.50627E+03 0.31916E+03 CR msg Press to continue pause text cl;con cl;redr con yo2 x 10 fi;0.001 gr ou x 10 te 1 O2 contours at ix=10 plane 0.50627E+03 0.31916E+03 CR msg Press to continue pause text cl;con cl;redr con tmp1 x 10 fi;0.001 gr ou x 10 te 1 gas phase temperature contours at ix=10 plane 0.50627E+03 0.31916E+03 CR msg Press to continue pause text cl;con cl;redr con tmp2 x 10 fi;0.001 gr ou x 10 te 1 solid phase temperature contours at ix=10 plane 0.50627E+03 0.31916E+03 CR msg Press to continue pause text cl;con cl;redr vec x 10 sh gr ou x 10 te 1 velocity vectors at ix=10 plane 0.50627E+03 0.31916E+03 CR msg Press to continue pause enduse DISPLAY The case considered is the 3d turbulent flow of combusting pulverised coal particles in a wall-fired power-station furnace. The calculation exploits two vertical planes of symmetry, so that only one quarter of the furnace is simulated. The burner arrangement comprises 3 rows of 4 horizontally-fired burners on opposing vertical walls, so that only 3 rows of 2 on a single wall are modelled in the actual solution domain. Each burner consists of 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 burner inflow boundary conditions are coded in the special ground routine furngr.f. The main combustion process is simulated using an IPSA-based 'equilibrium' 2-phase combustion model. Thermal radiation is accounted for by means of the radiosity model. A 2nd post-processing run can be performed which computes thermal NOx based on CHO equilibrium as determined by the CREK module. ENDDIS noxdbegin ns 7 nr 3 noxdend ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(3-D Coal Combustion In A Furnace ) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 14 ************************************************************ Group 2. Time dependence STEADY = T ************************************************************ Group 3. X-Direction Grid Spacing CARTES = T NX = 18 XULAST =10.509001 XFRAC(1)=0.109004 ;XFRAC(2)=0.21086 XFRAC(3)=0.302642 ;XFRAC(4)=0.373385 XFRAC(5)=0.40541 ;XFRAC(6)=0.437435 XFRAC(7)=0.469459 ;XFRAC(8)=0.501484 XFRAC(9)=0.556171 ;XFRAC(10)=0.610857 XFRAC(11)=0.665544 ;XFRAC(12)=0.72023 XFRAC(13)=0.752255 ;XFRAC(14)=0.78428 XFRAC(15)=0.816305 ;XFRAC(16)=0.84833 XFRAC(17)=0.905802 ;XFRAC(18)=1. ************************************************************ Group 4. Y-Direction Grid Spacing NY = 37 YVLAST =46.795998 YFRAC(1)=0.076177 ;YFRAC(2)=0.143336 YFRAC(3)=0.200049 ;YFRAC(4)=0.243844 YFRAC(5)=0.267796 ;YFRAC(6)=0.28218 YFRAC(7)=0.289372 ;YFRAC(8)=0.296564 YFRAC(9)=0.303756 ;YFRAC(10)=0.310948 YFRAC(11)=0.329131 ;YFRAC(12)=0.347314 YFRAC(13)=0.354506 ;YFRAC(14)=0.361698 YFRAC(15)=0.368889 ;YFRAC(16)=0.376081 YFRAC(17)=0.394265 ;YFRAC(18)=0.412448 YFRAC(19)=0.41964 ;YFRAC(20)=0.426831 YFRAC(21)=0.434023 ;YFRAC(22)=0.441215 YFRAC(23)=0.459398 ;YFRAC(24)=0.477581 YFRAC(25)=0.484773 ;YFRAC(26)=0.491965 YFRAC(27)=0.499157 ;YFRAC(28)=0.506349 YFRAC(29)=0.510292 ;YFRAC(30)=0.514234 YFRAC(31)=0.558464 ;YFRAC(32)=0.602693 YFRAC(33)=0.651723 ;YFRAC(34)=0.72342 YFRAC(35)=0.807207 ;YFRAC(36)=0.899958 YFRAC(37)=1. ************************************************************ Group 5. Z-Direction Grid Spacing PARAB = F NZ = 8 ZWLAST =6.398 ZFRAC(1)=0.04976 ;ZFRAC(2)=0.09952 ZFRAC(3)=0.250398 ;ZFRAC(4)=0.422935 ZFRAC(5)=0.604717 ;ZFRAC(6)=0.792775 ZFRAC(7)=0.896388 ;ZFRAC(8)=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(134)=HPOR ;NAME(135)=NPOR NAME(136)=EPOR ;NAME(137)=VPOR NAME(138)=PRPS ;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,N,N,N,Y) SOLUTN(H2,Y,Y,N,N,N,Y) SOLUTN(MIXF,Y,Y,N,N,N,Y) SOLUTN(C2,Y,Y,N,N,N,Y) SOLUTN(HPOR,Y,N,N,N,N,Y) SOLUTN(NPOR,Y,N,N,N,N,Y) SOLUTN(EPOR,Y,N,N,N,N,Y) SOLUTN(VPOR,Y,N,N,N,N,Y) SOLUTN(PRPS,Y,N,N,N,N,Y) 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 = 136 ;HPOR = 134 ;NPOR = 135 ;VPOR = 137 TEMP1 = 142 TEMP2 = 141 PRPS = 138 ************************************************************ 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,N,Y,Y) TERMS(H2,N,Y,Y,N,N,Y) TERMS(MIXF,N,Y,Y,Y,Y,Y) TERMS(C2,N,Y,Y,Y,N,Y) 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.93 ;RHO1B =0.07 ;RHO1C =0. EMISS =0. ;SCATT =0. RADIA =0. ;RADIB =0. 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. 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. CP1 =1000. ;CP2 =1000. ************************************************************ Group 10.Inter-Phase Transfer Processes CFIPS =1.0E+06 RLOLIM =0. ;CMDOT = GRND3 CMDTA =1000. ;CMDTB =0. CMDTC =0.157609 ;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)=0. ;FIINIT(U1)=0. FIINIT(U2)=0. ;FIINIT(V1)=2. FIINIT(V2)=2. ;FIINIT(W1)=1. FIINIT(W2)=1. ;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.33E+05 ;FIINIT(H2)=3.929812E+07 FIINIT(MIXF)=0. ;FIINIT(C2)=1. FIINIT(HPOR)=1. ;FIINIT(NPOR)=1. FIINIT(EPOR)=1. ;FIINIT(VPOR)=1. FIINIT(PRPS)=-1. ;FIINIT(ENUT)=1.0E-10 FIINIT(RHO1)=0.652786 ;FIINIT(TMP2)=1.0E-10 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)=1.0E-08 PATCH(SLOPE1 ,INIVAL, 1, 18, 1, 37, 1, 8, 1, 1) INIT(SLOPE1 ,HPOR,0. , GRND ) INIT(SLOPE1 ,NPOR,0. , GRND ) INIT(SLOPE1 ,EPOR,0. , GRND ) INIT(SLOPE1 ,VPOR,0. , GRND ) INIT(SLOPE1 ,PRPS,0. , GRND ) 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(KESOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1) COVAL(KESOURCE,KE , GRND4 , GRND4 ) COVAL(KESOURCE,EP , GRND4 , GRND4 ) PATCH(BURN1 ,CELL , 5, 8, 7, 10, 1, 1, 1, 1) COVAL(BURN1 ,U1 ,0. , GRND ) COVAL(BURN1 ,U2 ,0. , GRND ) COVAL(BURN1 ,V1 ,0. , GRND ) COVAL(BURN1 ,V2 ,0. , GRND ) COVAL(BURN1 ,W1 ,0. , GRND ) COVAL(BURN1 ,W2 ,0. , GRND ) COVAL(BURN1 ,GAS , FIXFLU , GRND ) COVAL(BURN1 ,FUE , FIXFLU , GRND ) COVAL(BURN1 ,KE ,0. , GRND ) COVAL(BURN1 ,EP ,0. , GRND ) COVAL(BURN1 ,H1 ,0. , GRND ) COVAL(BURN1 ,H2 ,0. , GRND ) COVAL(BURN1 ,MIXF,0. ,0. ) COVAL(BURN1 ,C2 ,0. ,1. ) PATCH(BURN2 ,CELL , 5, 8, 13, 16, 1, 1, 1, 1) COVAL(BURN2 ,U1 ,0. , GRND ) COVAL(BURN2 ,U2 ,0. , GRND ) COVAL(BURN2 ,V1 ,0. , GRND ) COVAL(BURN2 ,V2 ,0. , GRND ) COVAL(BURN2 ,W1 ,0. , GRND ) COVAL(BURN2 ,W2 ,0. , GRND ) COVAL(BURN2 ,GAS , FIXFLU , GRND ) COVAL(BURN2 ,FUE , FIXFLU , GRND ) COVAL(BURN2 ,KE ,0. , GRND ) COVAL(BURN2 ,EP ,0. , GRND ) COVAL(BURN2 ,H1 ,0. , GRND ) COVAL(BURN2 ,H2 ,0. , GRND ) COVAL(BURN2 ,MIXF,0. ,0. ) COVAL(BURN2 ,C2 ,0. ,1. ) PATCH(BURN3 ,CELL , 5, 8, 19, 22, 1, 1, 1, 1) COVAL(BURN3 ,U1 ,0. , GRND ) COVAL(BURN3 ,U2 ,0. , GRND ) COVAL(BURN3 ,V1 ,0. , GRND ) COVAL(BURN3 ,V2 ,0. , GRND ) COVAL(BURN3 ,W1 ,0. , GRND ) COVAL(BURN3 ,W2 ,0. , GRND ) COVAL(BURN3 ,GAS , FIXFLU , GRND ) COVAL(BURN3 ,FUE , FIXFLU , GRND ) COVAL(BURN3 ,KE ,0. , GRND ) COVAL(BURN3 ,EP ,0. , GRND ) COVAL(BURN3 ,H1 ,0. , GRND ) COVAL(BURN3 ,H2 ,0. , GRND ) COVAL(BURN3 ,MIXF,0. ,0. ) COVAL(BURN3 ,C2 ,0. ,1. ) PATCH(BURN4 ,CELL , 13, 16, 7, 10, 1, 1, 1, 1) COVAL(BURN4 ,U1 ,0. , GRND ) COVAL(BURN4 ,U2 ,0. , GRND ) COVAL(BURN4 ,V1 ,0. , GRND ) COVAL(BURN4 ,V2 ,0. , GRND ) COVAL(BURN4 ,W1 ,0. , GRND ) COVAL(BURN4 ,W2 ,0. , GRND ) COVAL(BURN4 ,GAS , FIXFLU , GRND ) COVAL(BURN4 ,FUE , FIXFLU , GRND ) COVAL(BURN4 ,KE ,0. , GRND ) COVAL(BURN4 ,EP ,0. , GRND ) COVAL(BURN4 ,H1 ,0. , GRND ) COVAL(BURN4 ,H2 ,0. , GRND ) COVAL(BURN4 ,MIXF,0. ,0. ) COVAL(BURN4 ,C2 ,0. ,1. ) PATCH(BURN5 ,CELL , 13, 16, 13, 16, 1, 1, 1, 1) COVAL(BURN5 ,U1 ,0. , GRND ) COVAL(BURN5 ,U2 ,0. , GRND ) COVAL(BURN5 ,V1 ,0. , GRND ) COVAL(BURN5 ,V2 ,0. , GRND ) COVAL(BURN5 ,W1 ,0. , GRND ) COVAL(BURN5 ,W2 ,0. , GRND ) COVAL(BURN5 ,GAS , FIXFLU , GRND ) COVAL(BURN5 ,FUE , FIXFLU , GRND ) COVAL(BURN5 ,KE ,0. , GRND ) COVAL(BURN5 ,EP ,0. , GRND ) COVAL(BURN5 ,H1 ,0. , GRND ) COVAL(BURN5 ,H2 ,0. , GRND ) COVAL(BURN5 ,MIXF,0. ,0. ) COVAL(BURN5 ,C2 ,0. ,1. ) PATCH(BURN6 ,CELL , 13, 16, 19, 22, 1, 1, 1, 1) COVAL(BURN6 ,U1 ,0. , GRND ) COVAL(BURN6 ,U2 ,0. , GRND ) COVAL(BURN6 ,V1 ,0. , GRND ) COVAL(BURN6 ,V2 ,0. , GRND ) COVAL(BURN6 ,W1 ,0. , GRND ) COVAL(BURN6 ,W2 ,0. , GRND ) COVAL(BURN6 ,GAS , FIXFLU , GRND ) COVAL(BURN6 ,FUE , FIXFLU , GRND ) COVAL(BURN6 ,KE ,0. , GRND ) COVAL(BURN6 ,EP ,0. , GRND ) COVAL(BURN6 ,H1 ,0. , GRND ) COVAL(BURN6 ,H2 ,0. , GRND ) COVAL(BURN6 ,MIXF,0. ,0. ) COVAL(BURN6 ,C2 ,0. ,1. ) PATCH(WALLW ,WWALL , 1, 1, 1, 37, 1, 8, 1, 1) COVAL(WALLW ,V1 , GRND2 ,0. ) COVAL(WALLW ,W1 , GRND2 ,0. ) COVAL(WALLW ,KE , GRND2 , GRND2 ) COVAL(WALLW ,EP , GRND2 , GRND2 ) PATCH(WALLA ,LWALL , 1, 18, 1, 6, 1, 1, 1, 1) COVAL(WALLA ,U1 , GRND2 ,0. ) COVAL(WALLA ,V1 , GRND2 ,0. ) COVAL(WALLA ,KE , GRND2 , GRND2 ) COVAL(WALLA ,EP , GRND2 , GRND2 ) PATCH(WALLB ,LWALL , 1, 4, 7, 10, 1, 1, 1, 1) COVAL(WALLB ,U1 , GRND2 ,0. ) COVAL(WALLB ,V1 , GRND2 ,0. ) COVAL(WALLB ,KE , GRND2 , GRND2 ) COVAL(WALLB ,EP , GRND2 , GRND2 ) PATCH(WALLC ,LWALL , 9, 12, 7, 10, 1, 1, 1, 1) COVAL(WALLC ,U1 , GRND2 ,0. ) COVAL(WALLC ,V1 , GRND2 ,0. ) COVAL(WALLC ,KE , GRND2 , GRND2 ) COVAL(WALLC ,EP , GRND2 , GRND2 ) PATCH(WALLD ,LWALL , 17, 18, 7, 10, 1, 1, 1, 1) COVAL(WALLD ,U1 , GRND2 ,0. ) COVAL(WALLD ,V1 , GRND2 ,0. ) COVAL(WALLD ,KE , GRND2 , GRND2 ) COVAL(WALLD ,EP , GRND2 , GRND2 ) PATCH(WALLE ,LWALL , 1, 18, 11, 12, 1, 1, 1, 1) COVAL(WALLE ,U1 , GRND2 ,0. ) COVAL(WALLE ,V1 , GRND2 ,0. ) COVAL(WALLE ,KE , GRND2 , GRND2 ) COVAL(WALLE ,EP , GRND2 , GRND2 ) PATCH(WALLF ,LWALL , 1, 4, 13, 16, 1, 1, 1, 1) COVAL(WALLF ,U1 , GRND2 ,0. ) COVAL(WALLF ,V1 , GRND2 ,0. ) COVAL(WALLF ,KE , GRND2 , GRND2 ) COVAL(WALLF ,EP , GRND2 , GRND2 ) PATCH(WALLG ,LWALL , 9, 12, 13, 16, 1, 1, 1, 1) COVAL(WALLG ,U1 , GRND2 ,0. ) COVAL(WALLG ,V1 , GRND2 ,0. ) COVAL(WALLG ,KE , GRND2 , GRND2 ) COVAL(WALLG ,EP , GRND2 , GRND2 ) PATCH(WALLH ,LWALL , 17, 18, 13, 16, 1, 1, 1, 1) COVAL(WALLH ,U1 , GRND2 ,0. ) COVAL(WALLH ,V1 , GRND2 ,0. ) COVAL(WALLH ,KE , GRND2 , GRND2 ) COVAL(WALLH ,EP , GRND2 , GRND2 ) PATCH(WALLI ,LWALL , 1, 18, 17, 18, 1, 1, 1, 1) COVAL(WALLI ,U1 , GRND2 ,0. ) COVAL(WALLI ,V1 , GRND2 ,0. ) COVAL(WALLI ,KE , GRND2 , GRND2 ) COVAL(WALLI ,EP , GRND2 , GRND2 ) PATCH(WALLJ ,LWALL , 1, 4, 19, 22, 1, 1, 1, 1) COVAL(WALLJ ,U1 , GRND2 ,0. ) COVAL(WALLJ ,V1 , GRND2 ,0. ) COVAL(WALLJ ,KE , GRND2 , GRND2 ) COVAL(WALLJ ,EP , GRND2 , GRND2 ) PATCH(WALLK ,LWALL , 9, 12, 19, 22, 1, 1, 1, 1) COVAL(WALLK ,U1 , GRND2 ,0. ) COVAL(WALLK ,V1 , GRND2 ,0. ) COVAL(WALLK ,KE , GRND2 , GRND2 ) COVAL(WALLK ,EP , GRND2 , GRND2 ) PATCH(WALLL ,LWALL , 17, 18, 19, 22, 1, 1, 1, 1) COVAL(WALLL ,U1 , GRND2 ,0. ) COVAL(WALLL ,V1 , GRND2 ,0. ) COVAL(WALLL ,KE , GRND2 , GRND2 ) COVAL(WALLL ,EP , GRND2 , GRND2 ) PATCH(WALLM ,LWALL , 1, 18, 23, 37, 1, 1, 1, 1) COVAL(WALLM ,U1 , GRND2 ,0. ) COVAL(WALLM ,V1 , GRND2 ,0. ) COVAL(WALLM ,KE , GRND2 , GRND2 ) COVAL(WALLM ,EP , GRND2 , GRND2 ) PATCH(HTLWLA ,LOW , 1, 18, 1, 6, 1, 1, 1, 1) COVAL(HTLWLA ,H1 , GRND ,6.996666E+05 ) PATCH(HTLWLB ,LOW , 1, 4, 7, 10, 1, 1, 1, 1) COVAL(HTLWLB ,H1 , GRND ,6.996666E+05 ) PATCH(HTLWLC ,LOW , 9, 12, 7, 10, 1, 1, 1, 1) COVAL(HTLWLC ,H1 , GRND ,6.996666E+05 ) PATCH(HTLWLD ,LOW , 17, 18, 7, 10, 1, 1, 1, 1) COVAL(HTLWLD ,H1 , GRND ,6.996666E+05 ) PATCH(HTLWLE ,LOW , 1, 18, 11, 12, 1, 1, 1, 1) COVAL(HTLWLE ,H1 , GRND ,6.996666E+05 ) PATCH(HTLWLF ,LOW , 1, 4, 13, 16, 1, 1, 1, 1) COVAL(HTLWLF ,H1 , GRND ,6.996666E+05 ) PATCH(HTLWLG ,LOW , 9, 12, 13, 16, 1, 1, 1, 1) COVAL(HTLWLG ,H1 , GRND ,6.996666E+05 ) PATCH(HTLWLH ,LOW , 17, 18, 13, 16, 1, 1, 1, 1) COVAL(HTLWLH ,H1 , GRND ,6.996666E+05 ) PATCH(HTLWLI ,LOW , 1, 18, 17, 18, 1, 1, 1, 1) COVAL(HTLWLI ,H1 , GRND ,6.996666E+05 ) PATCH(HTLWLJ ,LOW , 1, 4, 19, 22, 1, 1, 1, 1) COVAL(HTLWLJ ,H1 , GRND ,6.996666E+05 ) PATCH(HTLWLK ,LOW , 9, 12, 19, 22, 1, 1, 1, 1) COVAL(HTLWLK ,H1 , GRND ,6.996666E+05 ) PATCH(HTLWLL ,LOW , 17, 18, 19, 22, 1, 1, 1, 1) COVAL(HTLWLL ,H1 , GRND ,6.996666E+05 ) PATCH(HTLWLM ,LOW , 1, 18, 23, 37, 1, 1, 1, 1) COVAL(HTLWLM ,H1 , GRND ,6.996666E+05 ) PATCH(HWWA-EW ,EWALL , 10, 10, 31, 32, 1, 6, 1, 1) COVAL(HWWA-EW ,V1 , GRND2 ,0. ) COVAL(HWWA-EW ,W1 , GRND2 ,0. ) COVAL(HWWA-EW ,KE , GRND2 , GRND2 ) COVAL(HWWA-EW ,EP , GRND2 , GRND2 ) PATCH(HWWA-WW ,WWALL , 11, 11, 31, 32, 1, 6, 1, 1) COVAL(HWWA-WW ,V1 , GRND2 ,0. ) COVAL(HWWA-WW ,W1 , GRND2 ,0. ) COVAL(HWWA-WW ,KE , GRND2 , GRND2 ) COVAL(HWWA-WW ,EP , GRND2 , GRND2 ) PATCH(HWWB-EW ,EWALL , 10, 10, 33, 37, 3, 6, 1, 1) COVAL(HWWB-EW ,V1 , GRND2 ,0. ) COVAL(HWWB-EW ,W1 , GRND2 ,0. ) COVAL(HWWB-EW ,KE , GRND2 , GRND2 ) COVAL(HWWB-EW ,EP , GRND2 , GRND2 ) PATCH(HWWB-WW ,WWALL , 11, 11, 33, 37, 3, 6, 1, 1) COVAL(HWWB-WW ,V1 , GRND2 ,0. ) COVAL(HWWB-WW ,W1 , GRND2 ,0. ) COVAL(HWWB-WW ,KE , GRND2 , GRND2 ) COVAL(HWWB-WW ,EP , GRND2 , GRND2 ) PATCH(HWWC ,EWALL , 18, 18, 31, 32, 1, 8, 1, 1) COVAL(HWWC ,V1 , GRND2 ,0. ) COVAL(HWWC ,W1 , GRND2 ,0. ) COVAL(HWWC ,KE , GRND2 , GRND2 ) COVAL(HWWC ,EP , GRND2 , GRND2 ) PATCH(HWWD ,EWALL , 18, 18, 33, 37, 3, 8, 1, 1) COVAL(HWWD ,V1 , GRND2 ,0. ) COVAL(HWWD ,W1 , GRND2 ,0. ) COVAL(HWWD ,KE , GRND2 , GRND2 ) COVAL(HWWD ,EP , GRND2 , GRND2 ) PATCH(HWWUA ,EAST , 10, 10, 31, 32, 1, 6, 1, 1) COVAL(HWWUA ,U1 , GRND ,0. ) PATCH(HWWUB ,EAST , 10, 10, 33, 37, 3, 6, 1, 1) COVAL(HWWUB ,U1 , GRND ,0. ) PATCH(HTEHWA ,EAST , 10, 10, 31, 32, 1, 6, 1, 1) COVAL(HTEHWA ,H1 , GRND ,6.996666E+05 ) PATCH(HTEHWB ,EAST , 10, 10, 33, 37, 3, 6, 1, 1) COVAL(HTEHWB ,H1 , GRND ,6.996666E+05 ) PATCH(HTEHWC ,WEST , 11, 11, 31, 32, 1, 6, 1, 1) COVAL(HTEHWC ,H1 , GRND ,6.996666E+05 ) PATCH(HTEHWD ,WEST , 11, 11, 33, 37, 3, 6, 1, 1) COVAL(HTEHWD ,H1 , GRND ,6.996666E+05 ) PATCH(HTEHWE ,EAST , 18, 18, 31, 32, 1, 8, 1, 1) COVAL(HTEHWE ,H1 , GRND ,6.996666E+05 ) PATCH(HTEHWF ,EAST , 18, 18, 33, 37, 3, 8, 1, 1) COVAL(HTEHWF ,H1 , GRND ,6.996666E+05 ) PATCH(HTEWST ,WEST , 1, 1, 1, 37, 1, 8, 1, 1) COVAL(HTEWST ,H1 , GRND ,6.996666E+05 ) PATCH(OUTLET ,NORTH , 1, 18, 37, 37, 1, 8, 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 ,1.0E-04 ,0. ) COVAL(OUTLET ,FUE ,0.135 ,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. ,1. ) PATCH(GRAVTY ,PHASEM, 1, 18, 1, 37, 1, 8, 1, 1) COVAL(GRAVTY ,V1 , FIXFLU ,-9.81 ) COVAL(GRAVTY ,V2 , FIXFLU ,-9.81 ) XCYCLE = F WALLA =5.0E-03 ;WALLB =0. EGWF = T WALLCO = GRND2 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 2000 ;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)=30 ;LITER(U1)=4 LITER(U2)=4 ;LITER(V1)=4 LITER(V2)=2 ;LITER(W1)=3 LITER(W2)=3 ;LITER(GAS)=2 LITER(FUE)=2 ;LITER(SHAD)=2 LITER(KE)=20 ;LITER(EP)=20 LITER(H1)=20 ;LITER(H2)=20 LITER(MIXF)=20 ;LITER(C2)=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(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 ************************************************************ Group 17. Relaxation RELAX(P1,LINRLX,1.) RELAX(U1,FALSDT,0.037449) RELAX(U2,FALSDT,0.037449) RELAX(V1,FALSDT,0.037449) RELAX(V2,FALSDT,0.037449) RELAX(W1,FALSDT,0.037449) RELAX(W2,FALSDT,0.037449) RELAX(GAS,LINRLX,0.3) RELAX(FUE,LINRLX,0.3) RELAX(SHAD,LINRLX,0.3) RELAX(KE,FALSDT,7.489804E-03) RELAX(EP,FALSDT,7.489804E-03) RELAX(H1,FALSDT,0.074898) RELAX(H2,FALSDT,0.074898) RELAX(MIXF,FALSDT,0.018725) RELAX(C2,FALSDT,1.0E+09) RELAX(PRPS,LINRLX,1.) RELAX(ENUT,LINRLX,1.) RELAX(RHO1,LINRLX,0.3) 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,0.7) KELIN = 1 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-10 VARMAX(FUE)=1. ;VARMIN(FUE)=1.0E-10 VARMAX(SHAD)=1. ;VARMIN(SHAD)=1.0E-10 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)=3.929812E+07 VARMAX(MIXF)=0.157609 ;VARMIN(MIXF)=0. VARMAX(C2)=1. ;VARMIN(C2)=0. VARMAX(HPOR)=1.0E+10 ;VARMIN(HPOR)=-1.0E+10 VARMAX(NPOR)=1.0E+10 ;VARMIN(NPOR)=-1.0E+10 VARMAX(EPOR)=1.0E+10 ;VARMIN(EPOR)=-1.0E+10 VARMAX(VPOR)=1.0E+10 ;VARMIN(VPOR)=-1.0E+10 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10 VARMAX(RHO1)=1.0E+10 ;VARMIN(RHO1)=1.0E-10 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 NAMGRD =FURN GENK = T PARSOL = F ISG62 = 1 CSG9 ='COAL' SPEDAT(SET,FURN,XCB1,R,4.597) SPEDAT(SET,FURN,YCB1,R,13.878) SPEDAT(SET,FURN,XCB2,R,4.597) SPEDAT(SET,FURN,YCB2,R,16.926001) SPEDAT(SET,FURN,XCB3,R,4.597) SPEDAT(SET,FURN,YCB3,R,19.974001) SPEDAT(SET,FURN,XCB4,R,8.242001) SPEDAT(SET,FURN,YCB4,R,13.878) SPEDAT(SET,FURN,XCB5,R,8.242001) SPEDAT(SET,FURN,YCB5,R,16.926001) SPEDAT(SET,FURN,XCB6,R,8.242001) SPEDAT(SET,FURN,YCB6,R,19.974001) SPEDAT(SET,FURN,FGIN1,R,27.869837) SPEDAT(SET,FURN,FGIN2,R,22.33717) SPEDAT(SET,FURN,FGIN3,R,18.663469) SPEDAT(SET,FURN,FCOAL,R,10.816901) SPEDAT(SET,FURN,WGIN1,R,41.934875) SPEDAT(SET,FURN,WGIN2,R,21.355698) SPEDAT(SET,FURN,WGIN3,R,28.07303) SPEDAT(SET,FURN,WCOAL,R,21.354288) SPEDAT(SET,FURN,HGIN1,R,5.33E+05) SPEDAT(SET,FURN,HGIN2,R,3.385555E+05) SPEDAT(SET,FURN,HGIN3,R,5.33E+05) SPEDAT(SET,FURN,HCOAL,R,3.929812E+07) SPEDAT(SET,FURN,KEIN1,R,4.396334) SPEDAT(SET,FURN,KEIN2,R,1.140164) SPEDAT(SET,FURN,KEIN3,R,1.970238) SPEDAT(SET,FURN,EPIN1,R,74.533218) SPEDAT(SET,FURN,EPIN2,R,15.750166) SPEDAT(SET,FURN,EPIN3,R,13.250978) SPEDAT(SET,FURN,BX1,R,4.597) SPEDAT(SET,FURN,BX2,R,8.242001) SPEDAT(SET,FURN,BY1,R,13.878) SPEDAT(SET,FURN,BY1A,R,16.926001) SPEDAT(SET,FURN,BY1B,R,19.974001) SPEDAT(SET,FURN,DIAT,R,0.4064) SPEDAT(SET,FURN,DIAP,R,0.6604) SPEDAT(SET,FURN,DIAS,R,1.3462) SPEDAT(SET,FURN,GTHETA,R,55.) SPEDAT(SET,FURN,POINTY,R,8.238957) SPEDAT(SET,FURN,POINTZ,R,0.) SPEDAT(SET,FURN,OMEGT,R,10.) SPEDAT(SET,FURN,OMEGP,R,10.) SPEDAT(SET,FURN,OMEGS,R,10.) SPEDAT(SET,FURN,OMEGC,R,10.) SPEDAT(SET,FURN,TWALL,R,699.666626) SPEDAT(SET,FURN,OUTCOL,R,1.0E-04) SPEDAT(SET,FURN,OMGOFA,R,0.) SPEDAT(SET,FURN,XCOFA1,R,4.597) SPEDAT(SET,FURN,YCOFA1,R,23.021999) SPEDAT(SET,FURN,XCOFA2,R,8.242001) SPEDAT(SET,FURN,YCOFA2,R,23.021999) SPEDAT(SET,FURN,GROFA,R,1.417201) SPEDAT(SET,FURN,GWOFA,R,2.131711) SPEDAT(SET,FURN,GKEOFA,R,0.01136) SPEDAT(SET,FURN,GEPOFA,R,2.955652E-03) SPEDAT(SET,FURN,GHOFA,R,5.33E+05) SPEDAT(SET,FURN,DTUBE,R,0.0508) SPEDAT(SET,FURN,KWALL,R,0.03365) SPEDAT(SET,FURN,PRWALL,R,0.688) SPEDAT(SET,FURN,MUWALL,R,2.286E-05) SPEDAT(SET,FURN,G71,R,1.) SPEDAT(SET,FURN,JBCOPT,I,4) SPEDAT(SET,FURN,ANNINL,L,T) SPEDAT(SET,GXMONI,PLOTALL,L,T) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=2.634E+04 ;EX(U1)=3.091 EX(U2)=3.05 ;EX(V1)=7.427 EX(V2)=7.334 ;EX(W1)=6.17 EX(W2)=6.228 ;EX(GAS)=0.9797 EX(FUE)=0. ;EX(SHAD)=0. EX(KE)=13. ;EX(EP)=44.279999 EX(H1)=1.406E+06 ;EX(H2)=3.954E+07 EX(MIXF)=0.02949 ;EX(C2)=0. EX(HPOR)=0.9706 ;EX(NPOR)=0.9766 EX(EPOR)=0.9663 ;EX(VPOR)=0.9663 EX(PRPS)=0. ;EX(ENUT)=0. EX(RHO1)=0. ;EX(TMP2)=1401. EX(TMP1)=1406. ;EX(YH2O)=0.01858 EX(YH2)=0. ;EX(YN2)=0.7298 EX(YCO2)=0.1006 ;EX(YCO)=0. EX(YO2)=0.1308 ;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,Y,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,Y,N,N,N) OUTPUT(HPOR,Y,N,Y,N,N,N) OUTPUT(NPOR,Y,N,Y,N,N,N) OUTPUT(EPOR,Y,N,Y,N,N,N) OUTPUT(VPOR,Y,N,Y,N,N,N) OUTPUT(PRPS,Y,N,Y,N,N,N) OUTPUT(ENUT,N,N,Y,N,N,N) OUTPUT(RHO1,N,N,Y,N,N,N) OUTPUT(TMP2,Y,N,Y,N,N,N) OUTPUT(TMP1,Y,N,Y,N,N,N) 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,N,N,Y,N,N,N) OUTPUT(MDOT,N,N,Y,N,N,N) ************************************************************ Group 22. Monitor Print-Out IXMON = 4 ;IYMON = 11 ;IZMON = 7 NPRMON = 2000 ;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 = 2 ;IPLTL = -1 ;NPLT = 50 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