TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2018, Date 08/01/18
CPVNAM=VDI; SPPNAM=Core
************************************************************
Echo DISPLAY / USE settings
DISPLAY
Library case 644:
Tutorial:
Combustion model - CHEMKIN
ENDDIS
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 1. Run Title
TEXT(Combustion: CHEMKIN )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,1,30,50)
* Cylindrical-polar grid
CARTES=F
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(16)=CH4 ;NAME(17)=O2
NAME(18)=CO2 ;NAME(19)=H2O
NAME(20)=N2 ;NAME(145)=PRPS
NAME(146)=TEM1 ;NAME(147)=ENUL
NAME(148)=DEN1 ;NAME(149)=SPH1
NAME(150) =KOND
* Solved variables list
SOLVE(P1,V1,W1,CH4,O2,CO2,H2O,TEM1)
* Stored variables list
STORE(KOND,SPH1,DEN1,ENUL,PRPS,N2)
* Additional solver options
SOLUTN(P1,Y,Y,Y,N,N,Y)
************************************************************
Group 8. Terms & Devices
NEWRH1 = T
NEWENL = T
************************************************************
Group 9. Properties
RHO1 = GRND9
ENUL = GRND9
CP1 = GRND9
ENUT =0.
DVO1DT =3.41E-03
PRNDTL(CH4)=-GRND9 ;PRNDTL(O2)=-GRND9
PRNDTL(CO2)=-GRND9 ;PRNDTL(H2O)=-GRND9
PRNDTL(TEM1)=-GRND9
PRLH1A =0. ;PRLH1B =0.
PRLH1C =0.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(P1)=0. ;FIINIT(W1)=50.
FIINIT(O2)=0.24 ;FIINIT(N2)=0.76
FIINIT(PRPS)=-1. ;FIINIT(TEM1)=500.
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
PATCH(CHEMK1, VOLUME, 0, 0, 0, 0, 0, 0, 1, 1)
COVAL(CHEMK1, CH4, GRND9, GRND9)
COVAL(CHEMK1, O2, GRND9, GRND9)
COVAL(CHEMK1, CO2, GRND9, GRND9)
COVAL(CHEMK1, H2O, GRND9, GRND9)
COVAL(CHEMK1, TEM1, GRND9, GRND9)
This patch is attached to object OXIN
INLET(NOCPCK1, LOW, -1, 0, 0, 0, 0, 0, 1, 1)
VALUE(NOCPCK1, P1, GRND9)
VALUE(NOCPCK1, V1, 0.)
VALUE(NOCPCK1, W1, 300.)
VALUE(NOCPCK1, CH4, 0.)
VALUE(NOCPCK1, O2, 0.26)
VALUE(NOCPCK1, CO2, 0.)
VALUE(NOCPCK1, H2O, 0.)
VALUE(NOCPCK1, TEM1, GRND9)
This patch is attached to object FUIN1
INLET(NOCPCK2, LOW, -1, 0, 0, 0, 0, 0, 1, 1)
VALUE(NOCPCK2, P1, GRND9)
VALUE(NOCPCK2, V1, 0.)
VALUE(NOCPCK2, W1, 100.)
VALUE(NOCPCK2, CH4, 1.)
VALUE(NOCPCK2, O2, 0.)
VALUE(NOCPCK2, CO2, 0.)
VALUE(NOCPCK2, H2O, 0.)
VALUE(NOCPCK2, TEM1, GRND9)
This patch is attached to object FUIN2
INLET(NOCPCK3, LOW, -1, 0, 0, 0, 0, 0, 1, 1)
VALUE(NOCPCK3, P1, GRND9)
VALUE(NOCPCK3, V1, 0.)
VALUE(NOCPCK3, W1, 50.)
VALUE(NOCPCK3, CH4, 1.)
VALUE(NOCPCK3, O2, 0.)
VALUE(NOCPCK3, CO2, 0.)
VALUE(NOCPCK3, H2O, 0.)
VALUE(NOCPCK3, TEM1, GRND9)
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 4000
RESFAC =1.0E-04
************************************************************
Group 16. Terminate Iterations
LITER(TEM1)=50
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,1. )
RELAX(V1 ,FALSDT,1.0E-03 )
RELAX(W1 ,FALSDT,1.0E-03 )
RELAX(CH4 ,FALSDT,1.0E-03 )
RELAX(O2 ,FALSDT,1.0E-03 )
RELAX(CO2 ,FALSDT,1.0E-03 )
RELAX(H2O ,FALSDT,1.0E-03 )
RELAX(TEM1,FALSDT,5.0E-05 )
************************************************************
Group 18. Limits
VARMAX(CH4)=1. ;VARMIN(CH4)=0.
VARMAX(O2)=1. ;VARMIN(O2)=0.
VARMAX(CO2)=1. ;VARMIN(CO2)=0.
VARMAX(H2O)=1. ;VARMIN(H2O)=0.
VARMAX(N2)=1. ;VARMIN(N2)=0.
VARMAX(TEM1)=3000. ;VARMIN(TEM1)=273.
************************************************************
Group 19. EARTH Calls To GROUND Station
PARSOL = F
CHMKIN = T
ISG62 = 1
CHSOB =16.
CHSOC =1.
CSG4 ='ch4'
SPEDAT(SET,OUTPUT,NOFIELD,L,T)
SPEDAT(SET,NOCPCK2,TINLET,R,350.)
SPEDAT(SET,NOCPCK3,TINLET,R,350.)
SPEDAT(SET,NOCPCK1,TINLET,R,350.)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,MAXINC,TEM1,R,10.)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=0.7057 ;EX(V1)=5.717
EX(W1)=213.399994 ;EX(CH4)=0.4881
EX(O2)=0.01056 ;EX(CO2)=0.07457
EX(H2O)=0.06178 ;EX(N2)=0.349
EX(PRPS)=0.984 ;EX(TEM1)=1068.
EX(ENUL)=1.645 ;EX(DEN1)=2.936E-04
EX(SPH1)=1.297E+07 ;EX(KOND)=1.123E+04
************************************************************
Group 21. Print-out of Variables
OUTPUT(N2 ,Y,N,Y,N,Y,Y)
OUTPUT(DEN1,Y,N,Y,N,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 5 ;IZMON = 27
NPRMON = 100000
NPRMNT = 1
TSTSWP = -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000
YZPR = T
ISWPRF = 1 ;ISWPRL = 100000
No PATCHes used for this Group
************************************************************
Group 24. Dumps For Restarts
GVIEW(P,0.7682,0.145669,0.623418)
GVIEW(UP,-0.062667,0.986203,-0.153218)
GVIEW(VDIS,7.906215)
GVIEW(CENTRE,-2.33753,0.307689,14.231001)
> DOM, SIZE, 1.000000E-01, 3.000000E+00, 2.400000E+01
> DOM, MONIT, 5.000000E-02, 4.500000E-01, 1.019857E+01
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, INCREMENT, 1.000000E-02, 1.000000E-02, 1.000000E-02
> GRID, AUTO, T F F
> GRID, RSET_X_1, 1, 1.000000E+00
> GRID, RSET_Y_1, 10, 1.000000E+00
> GRID, RSET_Y_2, 4, 1.000000E+00
> GRID, RSET_Y_3, 10, 1.000000E+00
> GRID, RSET_Y_4, 6, 1.000000E+00
> GRID, RSET_Z_1, 6, 1.000000E+00
> GRID, RSET_Z_2, 44, 1.200000E+00
> OBJ, NAME, BLOC1
> OBJ, POSITION, 0.000000E+00, 1.000000E+00, 0.000000E+00
> OBJ, SIZE, TO_END, 4.000000E-01, 1.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polcu8
> OBJ, TYPE, BLOCKAGE
> OBJ, MATERIAL, 198,Solid with smooth-wall friction
> OBJ, NAME, OXIN
> OBJ, POSITION, 0.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ, SIZE, TO_END, 4.000000E-01, 0.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, poldef
> OBJ, TYPE, USER_DEFINED
> OBJ, PATCHES, NOCPCK1
> OBJ, NAME, FUIN1
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ, SIZE, TO_END, 1.000000E+00, 0.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, poldef
> OBJ, TYPE, USER_DEFINED
> OBJ, PATCHES, NOCPCK2
> OBJ, NAME, FUIN2
> OBJ, POSITION, 0.000000E+00, 1.400000E+00, 0.000000E+00
> OBJ, SIZE, TO_END, 1.000000E+00, 0.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, poldef
> OBJ, TYPE, USER_DEFINED
> OBJ, PATCHES, NOCPCK3
> OBJ, NAME, OUTLET
> OBJ, POSITION, 0.000000E+00, 0.000000E+00, AT_END
> OBJ, SIZE, TO_END, TO_END, 0.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polcubet
> OBJ, TYPE, OUTLET
> OBJ, PRESSURE, P_AMBIENT
> OBJ, TEMPERATURE, T_AMBIENT
> OBJ, COEFFICIENT, 1000.
> OBJ, NAME, NWALL
> OBJ, POSITION, 0.000000E+00, AT_END, 0.000000E+00
> OBJ, SIZE, TO_END, 0.000000E+00, TO_END
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polwire
> OBJ, TYPE, NULL
> OBJ, NAME, LWALL
> OBJ, POSITION, 0.000000E+00, 2.400000E+00, 0.000000E+00
> OBJ, SIZE, TO_END, TO_END, 0.000000E+00
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, polwire
> OBJ, TYPE, NULL
STOP