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
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
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/ 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: P-1 radiation )
************************************************************
************************************************************
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.5 ;TMP1B =0.13 ;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.5 ;TMP2B =0.5 ;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(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(RDSTGAS ,VOLUME, 1, 5, 1, 5, 1, 20, 1, 1)
COVAL(RDSTGAS ,CRAD,0.5 , GRND )
PATCH(RDSTFINE,VOLUME, 1, 5, 1, 5, 1, 20, 1, 1)
COVAL(RDSTFINE,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