TALK=T;RUN(1,1)
DISPLAY
The case considered is the zero-swirl gas-fired axisymmetric
combustion chamber used in the experiments of Owen et al (
16th Int. Symp. on Combustion, p105-117, 1976). The geometry
comprises a 0.122m diameter tube into which natural gas (96%CH4)
is injected at 288K axially to mix and burn with a coaxial
annular air stream at 750K. The inlet air/fuel equivalence ratio
is 0.9, and the combustor operating pressure is 3.8 atm.
The length of the solution domain is taken as 4 furnace diameters.
Calculations are made with one of 4 different combustion models:
(a) fast-chemistry model with a 1-step global irreversible
methane reaction to form CO2 and H2O; (b) as (a) but with an
assumed double-delta pdf to account for the influence of
concentration fluctuations; (c) eddy-break-up finite-rate
chemistry model with a 1-step global irreversible reaction; and
(d) as (c) but with a 2-step global reaction mechanism involving
the intermediate CO.
Radiative heat transfer is modelled via the P1 spherical-harmonics
approximation for a gray medium.
ENDDIS
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(Owen Furnace CH4 Diffusion Flame )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 302
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = F
NX = 1
XULAST =0.1
XFRAC(1)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 30
YVLAST =0.061
YFRAC(1)=0.051639 ;YFRAC(2)=0.103279
YFRAC(3)=0.154918 ;YFRAC(4)=0.206557
YFRAC(5)=0.258197 ;YFRAC(6)=0.309836
YFRAC(7)=0.361475 ;YFRAC(8)=0.413115
YFRAC(9)=0.464754 ;YFRAC(10)=0.516393
YFRAC(11)=0.542623 ;YFRAC(12)=0.568852
YFRAC(13)=0.595082 ;YFRAC(14)=0.621311
YFRAC(15)=0.647541 ;YFRAC(16)=0.67377
YFRAC(17)=0.7 ;YFRAC(18)=0.726229
YFRAC(19)=0.752459 ;YFRAC(20)=0.778688
YFRAC(21)=0.80082 ;YFRAC(22)=0.822951
YFRAC(23)=0.845082 ;YFRAC(24)=0.867213
YFRAC(25)=0.889344 ;YFRAC(26)=0.911475
YFRAC(27)=0.933607 ;YFRAC(28)=0.955738
YFRAC(29)=0.977869 ;YFRAC(30)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 40
ZWLAST =0.5
ZFRAC(1)=5.716314E-03 ;ZFRAC(2)=0.015085
ZFRAC(3)=0.026613 ;ZFRAC(4)=0.039811
ZFRAC(5)=0.054409 ;ZFRAC(6)=0.070231
ZFRAC(7)=0.087147 ;ZFRAC(8)=0.105061
ZFRAC(9)=0.123895 ;ZFRAC(10)=0.143587
ZFRAC(11)=0.164084 ;ZFRAC(12)=0.18534
ZFRAC(13)=0.207318 ;ZFRAC(14)=0.229983
ZFRAC(15)=0.253305 ;ZFRAC(16)=0.277258
ZFRAC(17)=0.301818 ;ZFRAC(18)=0.326962
ZFRAC(19)=0.352672 ;ZFRAC(20)=0.378929
ZFRAC(21)=0.405717 ;ZFRAC(22)=0.43302
ZFRAC(23)=0.460824 ;ZFRAC(24)=0.489116
ZFRAC(25)=0.517883 ;ZFRAC(26)=0.547115
ZFRAC(27)=0.5768 ;ZFRAC(28)=0.606928
ZFRAC(29)=0.63749 ;ZFRAC(30)=0.668476
ZFRAC(31)=0.699878 ;ZFRAC(32)=0.731688
ZFRAC(33)=0.763898 ;ZFRAC(34)=0.796501
ZFRAC(35)=0.82949 ;ZFRAC(36)=0.862858
ZFRAC(37)=0.896599 ;ZFRAC(38)=0.930707
ZFRAC(39)=0.965176 ;ZFRAC(40)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(5)=V1
NAME(7)=W1 ;NAME(12)=KE
NAME(13)=EP ;NAME(14)=H1
NAME(16)=CH4 ;NAME(17)=O2
NAME(18)=H2 ;NAME(19)=CO
NAME(20)=H2O ;NAME(21)=CO2
NAME(22)=N2 ;NAME(140)=EPKE
NAME(141)=SRAD ;NAME(142)=MCO2
NAME(143)=MCO ;NAME(144)=MMWT
NAME(145)=P1RS ;NAME(146)=SPH1
NAME(147)=F ;NAME(148)=TMP1
NAME(149)=DEN1 ;NAME(150)=VIST
* 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(V1,Y,Y,N,N,N,N)
SOLUTN(W1,Y,Y,N,N,N,N)
SOLUTN(KE,Y,Y,N,N,N,N)
SOLUTN(EP,Y,Y,N,N,N,N)
SOLUTN(H1,Y,N,N,N,N,Y)
SOLUTN(CH4,Y,Y,N,N,N,N)
SOLUTN(O2,Y,N,N,N,N,N)
SOLUTN(H2,Y,N,N,N,N,N)
SOLUTN(CO,Y,N,N,N,N,N)
SOLUTN(H2O,Y,N,N,N,N,N)
SOLUTN(CO2,Y,N,N,N,N,N)
SOLUTN(N2,Y,N,N,N,N,N)
SOLUTN(EPKE,Y,N,N,N,N,N)
SOLUTN(SRAD,Y,N,N,N,N,Y)
SOLUTN(MCO2,Y,N,N,N,N,Y)
SOLUTN(MCO,Y,N,N,N,N,Y)
SOLUTN(MMWT,Y,N,N,N,N,Y)
SOLUTN(P1RS,Y,N,N,N,N,Y)
SOLUTN(SPH1,Y,N,N,N,N,N)
SOLUTN(F,Y,Y,Y,N,N,Y)
SOLUTN(TMP1,Y,N,N,N,N,Y)
SOLUTN(DEN1,Y,N,N,N,N,Y)
SOLUTN(VIST,Y,N,N,N,N,Y)
DEN1 = 149
VIST = 150
TEMP1 = 148
************************************************************
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(V1,Y,Y,Y,Y,Y,Y)
TERMS(W1,Y,Y,Y,Y,Y,Y)
TERMS(KE,N,Y,Y,Y,Y,N)
TERMS(EP,N,Y,Y,Y,Y,N)
TERMS(CH4,N,Y,Y,Y,Y,Y)
TERMS(F,N,Y,Y,Y,N,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWRH1 = T
NEWENT = T
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
DENPCO = T
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 = GRND9 ;TMP1 = GRND9 ;EL1 = GRND4
TSURR =0. ;TEMP0 =0. ;PRESS0 =3.85035E+05
DVO1DT =0. ;DRH1DP =0.
RHO1A = GRND9 ;RHO1B =0. ;RHO1C =0.
TMP1A =0. ;TMP1B =0. ;TMP1C =0.
TMP2A =0. ;TMP2B =0. ;TMP2C =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
EL1A =0. ;EL1B =0. ;EL1C =0.
ENUL =4.2E-05 ;ENUT = GRND3
ENUTA =0. ;ENUTB =0. ;ENUTC =0.
IENUTA = 0
PRNDTL(V1)=1. ;PRNDTL(W1)=1.
PRNDTL(KE)=1. ;PRNDTL(EP)=1.
PRNDTL(CH4)=1. ;PRNDTL(F)=1.
PRT(V1)=1. ;PRT(W1)=1.
PRT(KE)=1. ;PRT(EP)=1.314
PRT(CH4)=1. ;PRT(F)=1.
CP1 = GRND7 ;CP2 =1.
CP1A =0. ;CP1B =0. ;CP1C =0. ;CP1D =0.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1.0E-10 ;FIINIT(V1)=1.0E-10
FIINIT(W1)=19.604807 ;FIINIT(KE)=8.64784
FIINIT(EP)=684.966309 ;FIINIT(H1)=4.724E+05
FIINIT(CH4)=1.0E-10 ;FIINIT(O2)=1.0E-10
FIINIT(H2)=1.0E-10 ;FIINIT(CO)=1.0E-10
FIINIT(H2O)=1.0E-10 ;FIINIT(CO2)=1.0E-10
FIINIT(N2)=1.0E-10 ;FIINIT(EPKE)=1.0E-10
FIINIT(SRAD)=1.0E-10 ;FIINIT(MCO2)=1.0E-10
FIINIT(MCO)=1.0E-10 ;FIINIT(MMWT)=1.0E-10
FIINIT(P1RS)=1.0E-10 ;FIINIT(SPH1)=1.0E-10
FIINIT(F)=0. ;FIINIT(TMP1)=750.
FIINIT(DEN1)=1.759791 ;FIINIT(VIST)=1.0E-10
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, 0, 0, 0, 0, 0, 0, 1, 1)
COVAL(KESOURCE,KE , GRND4 , GRND4 )
COVAL(KESOURCE,EP , GRND4 , GRND4 )
PATCH(SCRSPRSO,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1)
COVAL(SCRSPRSO,CH4 , GRND2 , GRND2 )
PATCH(SCRSF ,LOW , 1, 1, 1, 10, 1, 1, 1, 1)
COVAL(SCRSF ,P1 , FIXFLU , GRND1 )
COVAL(SCRSF ,V1 ,0. ,0. )
COVAL(SCRSF ,W1 ,0. ,0.906197 )
COVAL(SCRSF ,KE ,0. ,0.032848 )
COVAL(SCRSF ,EP ,0. ,0.310516 )
COVAL(SCRSF ,CH4 ,0. ,1. )
COVAL(SCRSF ,F ,0. ,1. )
PATCH(SCRSO ,LOW , 1, 1, 11, 20, 1, 1, 1, 1)
COVAL(SCRSO ,P1 , FIXFLU , GRND1 )
COVAL(SCRSO ,V1 ,0. ,0. )
COVAL(SCRSO ,W1 ,0. ,19.604807 )
COVAL(SCRSO ,KE ,0. ,15.373939 )
COVAL(SCRSO ,EP ,0. ,6190.065918 )
COVAL(SCRSO ,CH4 ,0. ,0. )
COVAL(SCRSO ,F ,0. ,0. )
PATCH(OUT ,HIGH , 1, 1, 1, 30, 40, 40, 1, 1)
COVAL(OUT ,P1 ,100. ,0. )
COVAL(OUT ,F ,0. , SAME )
PATCH(NWALL3 ,NWALL , 1, 1, 30, 30, 1, 40, 1, 1)
COVAL(NWALL3 ,W1 , GRND2 ,0. )
COVAL(NWALL3 ,KE , GRND2 , GRND2 )
COVAL(NWALL3 ,EP , GRND2 , GRND2 )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 600 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(V1)=10
LITER(W1)=10 ;LITER(KE)=20
LITER(EP)=20 ;LITER(CH4)=20
LITER(F)=20
ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03
ENDIT(W1)=1.0E-03 ;ENDIT(KE)=1.0E-03
ENDIT(EP)=1.0E-03 ;ENDIT(CH4)=1.0E-03
ENDIT(F)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,0.5)
RELAX(V1,FALSDT,6.375987E-04)
RELAX(W1,FALSDT,6.375987E-04)
RELAX(KE,LINRLX,0.5)
RELAX(EP,LINRLX,0.5)
RELAX(H1,LINRLX,1.)
RELAX(CH4,FALSDT,3.187994E-04)
RELAX(O2,LINRLX,1.)
RELAX(H2,LINRLX,1.)
RELAX(CO,LINRLX,1.)
RELAX(H2O,LINRLX,1.)
RELAX(CO2,LINRLX,1.)
RELAX(N2,LINRLX,1.)
RELAX(EPKE,LINRLX,1.)
RELAX(SRAD,LINRLX,1.)
RELAX(MCO2,LINRLX,1.)
RELAX(MCO,LINRLX,1.)
RELAX(MMWT,LINRLX,1.)
RELAX(P1RS,LINRLX,1.)
RELAX(SPH1,LINRLX,1.)
RELAX(F,LINRLX,0.5)
RELAX(TMP1,LINRLX,1.)
RELAX(DEN1,LINRLX,1.0E-02)
RELAX(VIST,LINRLX,1.)
KELIN = 3
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-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)=-1.0E+10
VARMAX(CH4)=1. ;VARMIN(CH4)=0.
VARMAX(O2)=1. ;VARMIN(O2)=0.
VARMAX(H2)=1. ;VARMIN(H2)=0.
VARMAX(CO)=1. ;VARMIN(CO)=0.
VARMAX(H2O)=1. ;VARMIN(H2O)=0.
VARMAX(CO2)=1. ;VARMIN(CO2)=0.
VARMAX(N2)=1. ;VARMIN(N2)=0.
VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
VARMAX(SRAD)=1.0E+10 ;VARMIN(SRAD)=-1.0E+10
VARMAX(MCO2)=1.0E+10 ;VARMIN(MCO2)=-1.0E+10
VARMAX(MCO)=1.0E+10 ;VARMIN(MCO)=-1.0E+10
VARMAX(MMWT)=1.0E+10 ;VARMIN(MMWT)=-1.0E+10
VARMAX(P1RS)=1.0E+10 ;VARMIN(P1RS)=-1.0E+10
VARMAX(SPH1)=1.0E+10 ;VARMIN(SPH1)=-1.0E+10
VARMAX(F)=1. ;VARMIN(F)=1.0E-10
VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=-1.0E+10
VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=1.0E-10
VARMAX(VIST)=1.0E+10 ;VARMIN(VIST)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
GENK = T
PARSOL = F
ISG62 = 1
CHSOB =16.
CSG4 ='sch4'
SPEDAT(SET,SCRS,IREACT,I,2)
SPEDAT(SET,SCRS,NRSYS,I,-2)
SPEDAT(SET,SCRS,NREAC,I,1)
SPEDAT(SET,SCRS,NELEM,I,4)
SPEDAT(SET,SCRS,IPDF,I,0)
SPEDAT(SET,SCRS,CEBP,R,6.)
SPEDAT(SET,SCRS,CEBS1,R,4.)
SPEDAT(SET,SCRS,CEBS2,R,4.)
SPEDAT(SET,SCRS,IFULIN,I,1)
SPEDAT(SET,SCRS,IPRMOD,I,0)
SPEDAT(SET,SCRS,ISRMOD,I,0)
SPEDAT(SET,SCRS,IPCHEM,I,0)
SPEDAT(SET,SCRS,ISCHEM,I,0)
SPEDAT(SET,SCRS,PREXP,R,1.585E+13)
SPEDAT(SET,SCRS,FUEXP,R,0.7)
SPEDAT(SET,SCRS,OXEXP,R,0.8)
SPEDAT(SET,SCRS,EACTP,R,4.8E+04)
SPEDAT(SET,SCRS,PREXS1,R,3.98E+14)
SPEDAT(SET,SCRS,FUEXS1,R,1.)
SPEDAT(SET,SCRS,OXEXS1,R,0.25)
SPEDAT(SET,SCRS,EACTS1,R,4.0E+04)
SPEDAT(SET,SCRS,PREXS2,R,0.)
SPEDAT(SET,SCRS,FUEXS2,R,1.)
SPEDAT(SET,SCRS,OXEXS2,R,1.)
SPEDAT(SET,SCRS,EACTS2,R,0.)
SPEDAT(SET,SCRS,NSPEC,I,7)
SPEDAT(SET,SCRS,CFUEL,C,CH4)
SPEDAT(SET,SCRS,COXID,C,O2)
SPEDAT(SET,SCRS,CFP1,C,CO2)
SPEDAT(SET,SCRS,CFP2,C,H2O)
SPEDAT(SET,CHEM,SIUNIT,L,T)
SPEDAT(SET,SCRS,YFU1,R,1.)
SPEDAT(SET,SCRS,YFU2,R,0.)
SPEDAT(SET,SCRS,YFU3,R,0.)
SPEDAT(SET,SCRS,YFU4,R,0.)
SPEDAT(SET,SCRS,YFU5,R,0.)
SPEDAT(SET,SCRS,YFU6,R,0.)
SPEDAT(SET,SCRS,YFU7,R,0.)
SPEDAT(SET,SCRS,TEMPFU,R,288.)
SPEDAT(SET,SCRS,YOX1,R,0.)
SPEDAT(SET,SCRS,YOX2,R,0.232)
SPEDAT(SET,SCRS,YOX3,R,0.)
SPEDAT(SET,SCRS,YOX4,R,0.)
SPEDAT(SET,SCRS,YOX5,R,0.)
SPEDAT(SET,SCRS,YOX6,R,0.)
SPEDAT(SET,SCRS,YOX7,R,0.768)
SPEDAT(SET,SCRS,TEMPOX,R,750.)
SPEDAT(SET,PRINT,NUMBER,I,10)
SPEDAT(SET,PRINT,COMMAND1,C,AVE_DEN1)
SPEDAT(SET,PRINT,COMMAND2,C,AVE_P1)
SPEDAT(SET,PRINT,COMMAND3,C,AVE_TMP1)
SPEDAT(SET,PRINT,COMMAND4,C,AVE_F)
SPEDAT(SET,PRINT,COMMAND5,C,AVE_V1)
SPEDAT(SET,PRINT,COMMAND6,C,AVE_W1)
SPEDAT(SET,PRINT,COMMAND7,C,AVE_EP)
SPEDAT(SET,PRINT,COMMAND8,C,AVE_KE)
SPEDAT(SET,PRINT,COMMAND9,C,AVE_VIST)
SPEDAT(SET,PRINT,COMMAND10,C,AVE_H1)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
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,N,Y,Y,Y,Y)
OUTPUT(V1,Y,N,Y,Y,Y,Y)
OUTPUT(W1,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,N,N,N)
OUTPUT(CH4,Y,N,Y,Y,Y,Y)
OUTPUT(O2,Y,N,Y,Y,N,N)
OUTPUT(H2,Y,N,Y,Y,N,N)
OUTPUT(CO,Y,N,Y,Y,N,N)
OUTPUT(H2O,Y,N,Y,Y,N,N)
OUTPUT(CO2,Y,N,Y,Y,N,N)
OUTPUT(N2,Y,N,Y,Y,N,N)
OUTPUT(EPKE,Y,N,Y,N,N,N)
OUTPUT(SRAD,Y,N,Y,N,N,N)
OUTPUT(MCO2,Y,N,Y,N,N,N)
OUTPUT(MCO,Y,N,Y,N,N,N)
OUTPUT(MMWT,Y,N,Y,N,N,N)
OUTPUT(P1RS,Y,N,Y,N,N,N)
OUTPUT(SPH1,Y,N,Y,N,N,N)
OUTPUT(F,Y,N,Y,Y,Y,Y)
OUTPUT(TMP1,Y,N,Y,N,Y,Y)
OUTPUT(DEN1,Y,N,Y,N,Y,Y)
OUTPUT(VIST,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 2 ;IZMON = 6
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 = -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