TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 3.4, Date 04/09/01
CPVNAM=VDI;SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
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
RADIATING PARTICLES IN A FLAME SYSTEM
Wall-fired furnace: 24 burners
Dry coal, no water in solid phase .
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.
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 1. Run Title
TEXT(Wall-fired model furnace: P-1 radiation model)
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,5,5,20)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = F
* Non-default variable names
NAME(105) =CRAD ; NAME(106) =MCCO
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
* Solved variables list
SOLVE(P1 ,U1 ,U2 ,V1 ,V2 ,W1 ,W2 ,R1 )
SOLVE(R2 ,RS ,H1 ,H2 ,CRAD,YH2O,YCO2,YCO )
SOLVE(YO2 ,YCH4,COL2,CHA2,YHCN,YNOX)
* Stored variables list
STORE(EPKE,ENUT,EL1 ,CMDO,VRAT,CFIP,TMP1,TMP2)
STORE(NUSS,VREL,CD ,REYN,APRJ,DEN1,SIZE,KD )
STORE(KC ,HCOF,VCOK,ASH2,YN2 ,VCOM,CCOM,WMIX)
STORE(BURN,VCOD,CCOD,CCOK,FK1 ,MOX ,RNOX,BRN1)
STORE(BRN2,XXX ,MCCO)
* Additional solver options
SOLUTN(P1 ,Y,Y,Y,N,N,Y)
SOLUTN(R1 ,Y,Y,N,Y,N,Y)
SOLUTN(R2 ,Y,Y,N,Y,N,Y)
TURMOD(KEMODL)
************************************************************
Group 8. Terms & Devices
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)
NEWRH1 = T
************************************************************
Group 9. Properties
RHO1 = GRND5
RHO1A = 0.000000E+00 ;RHO1B = 3.496503E-03
RHO1C = 7.142857E-01
PRESS0 = 1.000000E+05
TMP1 = GRND
TMP1A = 5.000000E-01 ;TMP1B = 1.300000E-01
TMP1C = 0.000000E+00
CP1 = 1.050000E+03
RHO2 = 1.000000E+03
RHO2A = 3.000000E-01 ;RHO2B = 0.000000E+00
RHO2C = 0.000000E+00
TMP2 = GRND
TMP2A = 5.000000E-01 ;TMP2B = 5.000000E-01
TMP2C = 0.000000E+00
CP2 = 1.800000E+03
ENUL = 1.000000E-05
DRH1DP = GRND5
PRNDTL(R1 ) = 1.000000E+09 ;PRNDTL(R2 ) = 1.000000E+09
PRNDTL(CRAD) = GRND
PRT (R1 ) = 1.000000E+09 ;PRT (R2 ) = 1.000000E+09
PRT (EP ) = 1.314000E+00 ;PRT (CRAD) = 1.000000E+09
EMISS = 0.000000E+00 ;SCATT = 0.000000E+00
************************************************************
Group 10.Inter-Phase Transfer Processes
CFIPS = GRND7
RLOLIM = 1.000000E-09 ;CMDOT = GRND
CFIPA = 1.000000E-09 ;CFIPB = 2.500000E-05
CFIPC = 1.000000E+02 ;CFIPD = 0.000000E+00
CINT (H1 ) = 1.000000E-20 ;CINT (H2 ) = 1.000000E-20
************************************************************
Echo PLANT settings
PLANTBEGIN
TMP1A is gas absorption coefficient, 1/m
TMP1B is gas scattering coefficient, 1/m
TMP2A is particle emmisivity, dimensionless
TMP2B is particle reflectivity, 1-TMP2A
NAMSAT=MOSG
PRNDTL(CRAD)=GRND
LAMPR(CRAD)=1e-5*DEN1*3*(TMP1A+TMP1B+APRJ*TMP2A*(2-TMP2A))
* 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=TMP1A*(CRAD-4.*5.68e-8*TMP1**4)+$
APRJ*TMP2A*(CRAD-4.*5.68e-8*TMP2**4)
COVAL(RADSOR,H1,FIXFLU,GRND)
VAL=TMP1A*4*APRJ*(CRAD/4.-5.68e-8*TMP2**4)
COVAL(RADSOR,H2,FIXFLU,GRND)
* Radiosity sources
* from gas
PATCH(RDSTGAS,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP)
VAL=4.*5.68e-8*TMP1**4
COVAL(RDSTGAS,CRAD,TMP1A,GRND)
* from fines
PATCH(RDSTFINE,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP)
CO =APRJ*TMP2A
VAL=4.*5.68e-8*TMP2**4
COVAL(RDSTFINE,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/ $
EXP(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 11.Initialise Var/Porosity Fields
FIINIT(R1 ) = 5.000000E-01 ;FIINIT(R2 ) = 5.000000E-01
FIINIT(RS ) = 5.000000E-01
No PATCHes used for this Group
INIADD = F
************************************************************
Group 12. Convection and diffusion adjustments
No PATCHes used for this Group
************************************************************
Group 13. Boundary & Special Sources
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 1250
RESFAC = 1.000000E-03
************************************************************
Group 16. Terminate Iterations
LITER (R1 ) = 1 ;LITER (R2 ) = 1
LITER (RS ) = 1
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX, 5.000000E-01)
RELAX(U1 ,FALSDT, 1.000000E-02)
RELAX(U2 ,FALSDT, 1.000000E-02)
RELAX(V1 ,FALSDT, 1.000000E-02)
RELAX(V2 ,FALSDT, 1.000000E-02)
RELAX(W1 ,FALSDT, 1.000000E-02)
RELAX(W2 ,FALSDT, 1.000000E-02)
RELAX(R1 ,LINRLX, 5.000000E-01)
RELAX(R2 ,LINRLX, 5.000000E-01)
RELAX(KE ,LINRLX, 2.000000E-01)
RELAX(EP ,LINRLX, 2.000000E-01)
RELAX(H1 ,FALSDT, 1.000000E-02)
RELAX(H2 ,FALSDT, 1.000000E-02)
RELAX(CRAD,FALSDT, 1.000000E-02)
RELAX(YH2O,FALSDT, 1.000000E-02)
RELAX(YCO2,FALSDT, 1.000000E-02)
RELAX(YCO ,FALSDT, 1.000000E-02)
RELAX(YO2 ,FALSDT, 1.000000E-02)
RELAX(YCH4,FALSDT, 1.000000E-02)
RELAX(DEN1,LINRLX, 1.000000E-01)
RELAX(COL2,FALSDT, 1.000000E-02)
RELAX(CHA2,FALSDT, 1.000000E-02)
RELAX(CMDO,LINRLX, 3.000000E-01)
RELAX(YHCN,FALSDT, 1.000000E-02)
RELAX(YNOX,FALSDT, 1.000000E-02)
KELIN = 3
************************************************************
Group 18. Limits
VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06
VARMAX(U2 ) = 1.000000E+06 ;VARMIN(U2 ) =-1.000000E+06
VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06
VARMAX(V2 ) = 1.000000E+06 ;VARMIN(V2 ) =-1.000000E+06
VARMAX(W1 ) = 1.000000E+06 ;VARMIN(W1 ) =-1.000000E+06
VARMAX(W2 ) = 1.000000E+06 ;VARMIN(W2 ) =-1.000000E+06
VARMAX(R1 ) = 1.000000E+00 ;VARMIN(R1 ) = 1.000000E-06
VARMAX(R2 ) = 1.000000E+00 ;VARMIN(R2 ) = 1.000000E-09
VARMAX(H1 ) = 1.000000E+10 ;VARMIN(H1 ) = 0.000000E+00
VARMAX(H2 ) = 1.000000E+10 ;VARMIN(H2 ) = 0.000000E+00
VARMAX(CCOM) = 1.000000E+10 ;VARMIN(CCOM) = 0.000000E+00
VARMAX(VCOM) = 1.000000E+10 ;VARMIN(VCOM) = 0.000000E+00
VARMAX(YN2 ) = 1.000000E+00 ;VARMIN(YN2 ) = 0.000000E+00
VARMAX(ASH2) = 1.000000E+00 ;VARMIN(ASH2) = 0.000000E+00
VARMAX(YH2O) = 1.000000E+00 ;VARMIN(YH2O) = 0.000000E+00
VARMAX(YCO2) = 1.000000E+00 ;VARMIN(YCO2) = 0.000000E+00
VARMAX(YCO ) = 1.000000E+00 ;VARMIN(YCO ) = 0.000000E+00
VARMAX(YO2 ) = 1.000000E+00 ;VARMIN(YO2 ) = 0.000000E+00
VARMAX(YCH4) = 1.000000E+00 ;VARMIN(YCH4) = 0.000000E+00
VARMAX(COL2) = 1.000000E+00 ;VARMIN(COL2) = 0.000000E+00
VARMAX(CHA2) = 1.000000E+00 ;VARMIN(CHA2) = 0.000000E+00
VARMAX(TMP2) = 1.000000E+10 ;VARMIN(TMP2) = 3.000000E+02
VARMAX(TMP1) = 1.000000E+10 ;VARMIN(TMP1) = 3.000000E+02
VARMAX(CMDO) = 1.000000E+10 ;VARMIN(CMDO) = 0.000000E+00
VARMAX(YHCN) = 1.000000E+00 ;VARMIN(YHCN) = 0.000000E+00
VARMAX(YNOX) = 1.000000E+00 ;VARMIN(YNOX) = 0.000000E+00
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD = T ;USEGRX = T
NAMSAT =MOSG
GENK = T
ASAP = T
************************************************************
Group 20. Preliminary Printout
ECHO = T
************************************************************
Group 21. Print-out of Variables
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)
************************************************************
Group 22. Monitor Print-Out
IXMON = 3 ;IYMON = 4 ;IZMON = 9
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
NZPRIN = 1
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
NOWIPE = T
GVIEW(P,-8.292519E-01,3.366233E-01,4.461233E-01)
GVIEW(UP,4.125013E-01,-1.699202E-01,8.949692E-01)
> DOM, SIZE, 1.000000E+00, 1.000000E+00, 1.000000E+01
> DOM, MONIT, 5.000000E-01, 7.000000E-01, 4.250000E+00
> DOM, SCALE, 3.000000E+00, 3.000000E+00, 1.000000E+00
> DOM, SNAPSIZE, 1.000000E-02
> DOM, RELAX, 3.000000E-01
> OBJ, NAME, OUTLET
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 1.000000E+01
> OBJ, SIZE, 1.000000E+00, 1.000000E+00, 0.000000E+00
> OBJ, CLIPART, cube12t
> OBJ, ROTATION24, 1
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, 0.000000E+00
> OBJ, TEMPERATURE, -1.026000E+04
> OBJ, COEFFICIENT, 1.000000E+00
> OBJ, TEMPERATURE-2,-1.026000E+04
> OBJ, COEFFICIENT-2, 1.000000E+04
> OBJ, TURBULENCE, -1.026000E+04,-1.026000E+04
> OBJ, NAME, B4
> OBJ, POSITION, 6.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ, SIZE, 2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ, CLIPART, cube3t
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, INLET
> OBJ, DENSITY, 6.930000E-01
> OBJ, DENSITY-2, 7.000000E-02
> OBJ, VELOCITY, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, VELOCITY-2, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, TEMPERATURE, 5.250000E+05
> OBJ, TEMPERATURE-2, 9.000000E+05
> OBJ, INLET_YO2, 3.500000E-01
> OBJ, INLET_COL2, 8.000000E-01
> OBJ, TURB-INTENS, 5.000000E+00
> OBJ, NAME, B6
> OBJ, POSITION, 2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ, SIZE, 2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ, CLIPART, cube3t
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, INLET
> OBJ, DENSITY, 6.930000E-01
> OBJ, DENSITY-2, 7.000000E-02
> OBJ, VELOCITY, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, VELOCITY-2, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, TEMPERATURE, 5.250000E+05
> OBJ, TEMPERATURE-2, 9.000000E+05
> OBJ, INLET_YO2, 3.500000E-01
> OBJ, INLET_COL2, 8.000000E-01
> OBJ, TURB-INTENS, 5.000000E+00
> OBJ, NAME, B7
> OBJ, POSITION, 2.000000E-01, 0.000000E+00, 1.500000E+00
> OBJ, SIZE, 2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ, CLIPART, cube3t
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, INLET
> OBJ, DENSITY, 6.930000E-01
> OBJ, DENSITY-2, 7.000000E-02
> OBJ, VELOCITY, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, VELOCITY-2, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, TEMPERATURE, 5.250000E+05
> OBJ, TEMPERATURE-2, 9.000000E+05
> OBJ, INLET_YO2, 3.500000E-01
> OBJ, INLET_COL2, 8.000000E-01
> OBJ, TURB-INTENS, 5.000000E+00
> OBJ, NAME, B9
> OBJ, POSITION, 6.000000E-01, 0.000000E+00, 1.500000E+00
> OBJ, SIZE, 2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ, CLIPART, cube3t
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, INLET
> OBJ, DENSITY, 6.930000E-01
> OBJ, DENSITY-2, 7.000000E-02
> OBJ, VELOCITY, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, VELOCITY-2, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, TEMPERATURE, 5.250000E+05
> OBJ, TEMPERATURE-2, 9.000000E+05
> OBJ, INLET_YO2, 3.500000E-01
> OBJ, INLET_COL2, 8.000000E-01
> OBJ, TURB-INTENS, 5.000000E+00
> OBJ, NAME, B8
> OBJ, POSITION, 6.000000E-01, 0.000000E+00, 2.500000E+00
> OBJ, SIZE, 2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ, CLIPART, cube3t
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, INLET
> OBJ, DENSITY, 6.930000E-01
> OBJ, DENSITY-2, 7.000000E-02
> OBJ, VELOCITY, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, VELOCITY-2, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, TEMPERATURE, 5.250000E+05
> OBJ, TEMPERATURE-2, 9.000000E+05
> OBJ, INLET_YO2, 3.500000E-01
> OBJ, INLET_COL2, 8.000000E-01
> OBJ, TURB-INTENS, 5.000000E+00
> OBJ, NAME, B10
> OBJ, POSITION, 2.000000E-01, 0.000000E+00, 2.500000E+00
> OBJ, SIZE, 2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ, CLIPART, cube3t
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, INLET
> OBJ, DENSITY, 6.930000E-01
> OBJ, DENSITY-2, 7.000000E-02
> OBJ, VELOCITY, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, VELOCITY-2, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, TEMPERATURE, 5.250000E+05
> OBJ, TEMPERATURE-2, 9.000000E+05
> OBJ, INLET_YO2, 3.500000E-01
> OBJ, INLET_COL2, 8.000000E-01
> OBJ, TURB-INTENS, 5.000000E+00
> OBJ, NAME, B11
> OBJ, POSITION, 6.000000E-01, 0.000000E+00, 3.500000E+00
> OBJ, SIZE, 2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ, CLIPART, cube3t
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, INLET
> OBJ, DENSITY, 6.930000E-01
> OBJ, DENSITY-2, 7.000000E-02
> OBJ, VELOCITY, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, VELOCITY-2, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, TEMPERATURE, 5.250000E+05
> OBJ, TEMPERATURE-2, 9.000000E+05
> OBJ, INLET_YO2, 3.500000E-01
> OBJ, INLET_COL2, 8.000000E-01
> OBJ, TURB-INTENS, 5.000000E+00
> OBJ, NAME, B12
> OBJ, POSITION, 2.000000E-01, 0.000000E+00, 3.500000E+00
> OBJ, SIZE, 2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ, CLIPART, cube3t
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, INLET
> OBJ, DENSITY, 6.930000E-01
> OBJ, DENSITY-2, 7.000000E-02
> OBJ, VELOCITY, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, VELOCITY-2, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, TEMPERATURE, 5.250000E+05
> OBJ, TEMPERATURE-2, 9.000000E+05
> OBJ, INLET_YO2, 3.500000E-01
> OBJ, INLET_COL2, 8.000000E-01
> OBJ, TURB-INTENS, 5.000000E+00
> OBJ, NAME, B13
> OBJ, POSITION, 6.000000E-01, 0.000000E+00, 4.500000E+00
> OBJ, SIZE, 2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ, CLIPART, cube3t
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, INLET
> OBJ, DENSITY, 6.930000E-01
> OBJ, DENSITY-2, 7.000000E-02
> OBJ, VELOCITY, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, VELOCITY-2, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, TEMPERATURE, 5.250000E+05
> OBJ, TEMPERATURE-2, 9.000000E+05
> OBJ, INLET_YO2, 3.500000E-01
> OBJ, INLET_COL2, 8.000000E-01
> OBJ, TURB-INTENS, 5.000000E+00
> OBJ, NAME, B14
> OBJ, POSITION, 2.000000E-01, 0.000000E+00, 4.500000E+00
> OBJ, SIZE, 2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ, CLIPART, cube3t
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, INLET
> OBJ, DENSITY, 6.930000E-01
> OBJ, DENSITY-2, 7.000000E-02
> OBJ, VELOCITY, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, VELOCITY-2, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, TEMPERATURE, 5.250000E+05
> OBJ, TEMPERATURE-2, 9.000000E+05
> OBJ, INLET_YO2, 3.500000E-01
> OBJ, INLET_COL2, 8.000000E-01
> OBJ, TURB-INTENS, 5.000000E+00
> OBJ, NAME, B15
> OBJ, POSITION, 6.000000E-01, 0.000000E+00, 5.500000E+00
> OBJ, SIZE, 2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ, CLIPART, cube3t
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, INLET
> OBJ, DENSITY, 6.930000E-01
> OBJ, DENSITY-2, 7.000000E-02
> OBJ, VELOCITY, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, VELOCITY-2, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, TEMPERATURE, 5.250000E+05
> OBJ, TEMPERATURE-2, 9.000000E+05
> OBJ, INLET_YO2, 3.500000E-01
> OBJ, INLET_COL2, 8.000000E-01
> OBJ, TURB-INTENS, 5.000000E+00
> OBJ, NAME, B16
> OBJ, POSITION, 2.000000E-01, 0.000000E+00, 5.500000E+00
> OBJ, SIZE, 2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ, CLIPART, cube3t
> OBJ, ROTATION24, 1
> OBJ, GRID, 2
> OBJ, TYPE, INLET
> OBJ, DENSITY, 6.930000E-01
> OBJ, DENSITY-2, 7.000000E-02
> OBJ, VELOCITY, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, VELOCITY-2, 0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ, TEMPERATURE, 5.250000E+05
> OBJ, TEMPERATURE-2, 9.000000E+05
> OBJ, INLET_YO2, 3.500000E-01
> OBJ, INLET_COL2, 8.000000E-01
> OBJ, TURB-INTENS, 5.000000E+00
STOP