DISPLAY

  Library case 643:

  Tutorial:

  Combustion model - ESCRS  

  ENDDIS
  
 TALK=T;RUN( 1, 1)
 
 ************************************************************
   Q1 created by VDI menu, Version 2018, Date 08/01/18
 CPVNAM=VDI; SPPNAM=Core
 ************************************************************
 IRUNN = 1 ;LIBREF = 0
 ************************************************************
  Group 1. Run Title
 TEXT(ESCRS - Fast Chemistry                  )
 ************************************************************
  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,40)
    * 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)=H2 ;NAME(19)=CO
 NAME(20)=H2O ;NAME(21)=CO2
 NAME(22)=N2 ;NAME(141)=TMP1
 NAME(142)=MCO2 ;NAME(143)=MCO
 NAME(144)=MMWT ;NAME(145)=SPH1
 NAME(146)=F ;NAME(147)=EPKE
 NAME(148)=DEN1 ;NAME(149)=EL1
 NAME(150) =ENUT
    * Solved variables list
 SOLVE(P1,V1,W1,F)
    * Stored variables list
 STORE(ENUT,EL1,DEN1,EPKE,SPH1,MMWT,MCO,MCO2)
 STORE(TMP1,N2,CO2,H2O,CO,H2,O2,CH4)
 STORE(H1)
    * Additional solver options
 SOLUTN(P1,Y,Y,Y,N,N,Y)
 SOLUTN(F,Y,Y,Y,N,N,Y)
 TURMOD(KECHEN)
 
 ************************************************************
  Group 8. Terms & Devices
 NEWRH1 = T
 DENPCO = T
 ************************************************************
  Group 9. Properties
 PRESS0 =3.8E+05
 RHO1 = GRND9
 RHO1A = GRND9 ;RHO1B =0.
 RHO1C =0.
 ENUL =4.2E-05
 CP1 = GRND7
 TMP1 = GRND9
 PRT(KE)=0.75 ;PRT(EP)=1.15
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(P1)=0. ;FIINIT(W1)=5.5
 FIINIT(KE)=9.458 ;FIINIT(EP)=4725.
 FIINIT(F)=0.186
   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
 
   This patch is attached to object SCRSO
 INLET(SCRSO, LOW, -1, 0, 0, 0, 0, 0, 1, 1)
 VALUE(SCRSO, P1, GRND1)
 VALUE(SCRSO, V1, 0.)
 VALUE(SCRSO, W1, 19.6)
 VALUE(SCRSO, KE, 15.4)
 VALUE(SCRSO, EP, 6190.)
 VALUE(SCRSO, CH4, 0.)
 VALUE(SCRSO, F, 0.)
 
   This patch is attached to object SCRSF
 INLET(SCRSF, LOW, -1, 0, 0, 0, 0, 0, 1, 1)
 VALUE(SCRSF, P1, GRND1)
 VALUE(SCRSF, V1, 0.)
 VALUE(SCRSF, W1, 0.91)
 VALUE(SCRSF, KE, 0.033)
 VALUE(SCRSF, EP, 0.31)
 VALUE(SCRSF, CH4, 1.)
 VALUE(SCRSF, F, 1.)
 
 EGWF = T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 400
 SARAH =1.
 RESFAC =1.0E-04
 ************************************************************
  Group 16. Terminate Iterations
 LITER(KE)=10 ;LITER(EP)=10
 ************************************************************
  Group 17. Relaxation
 RELAX(P1 ,LINRLX,1. )
 RELAX(KE ,LINRLX,0.5 )
 RELAX(EP ,LINRLX,0.5 )
 RELAX(F ,LINRLX,0.5 )
 RELAX(DEN1,LINRLX,1.0E-02 )
 KELIN = 3
 ************************************************************
  Group 18. Limits
 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.
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 GENK = T
 PARSOL = F
 IENUTA = 2
 ISG62 = 1
 CHSOB =16.
 CSG4    ='sch4'
 SPEDAT(SET,SCRS,IREACT,I,1)
 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,0.)
 SPEDAT(SET,SCRS,CEBS1,R,0.)
 SPEDAT(SET,SCRS,CEBS2,R,0.)
 SPEDAT(SET,SCRS,IFULIN,I,0)
 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,OUTPUT,NOFIELD,L,T)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=31.299999 ;EX(V1)=0.5365
 EX(W1)=12.21 ;EX(KE)=6.236
 EX(EP)=1453. ;EX(H1)=7.251E+05
 EX(CH4)=0.1186 ;EX(O2)=0.09841
 EX(H2)=0. ;EX(CO)=0.
 EX(H2O)=0.05644 ;EX(CO2)=0.06894
 EX(N2)=0.6576 ;EX(TMP1)=1410.
 EX(MCO2)=0.04011 ;EX(MCO)=0.
 EX(MMWT)=26.35 ;EX(SPH1)=1672.
 EX(F)=0.1437 ;EX(EPKE)=177.5
 EX(DEN1)=1.023 ;EX(EL1)=3.651E-03
 EX(ENUT)=4.547E-03
 ************************************************************
  Group 21. Print-out of Variables
 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,Y,Y)
 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(DEN1,Y,N,Y,N,Y,Y)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 1 ;IYMON = 10 ;IZMON = 25
 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,-1.,0.,0.)
 GVIEW(UP,0.,1.,0.)
 GVIEW(VDIS,0.236081)
 GVIEW(CENTRE,3.044919E-03,0.0305,0.244)
 
> DOM,    SIZE,        1.000000E-01, 6.100000E-02, 5.000000E-01
> DOM,    MONIT,       5.000000E-02, 2.992500E-02, 2.517498E-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,     10, 1.000000E+00
> GRID,   RSET_Y_3,     10, 1.000000E+00
> GRID,   RSET_Z_1,     40, 1.400000E+00
> DOM,    INI_AMB,    YES
> DOM,    INI_BUOY,   YES
 
> OBJ,    NAME,        SCRSO
> OBJ,    POSITION,    0.000000E+00, 3.150000E-02, 0.000000E+00
> OBJ,    SIZE,        TO_END,       1.600000E-02, 0.000000E+00
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    poldef
> OBJ,    TYPE,        USER_DEFINED
> OBJ,    PATCHES,     SCRSO
 
> OBJ,    NAME,        SCRSF
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        TO_END,       3.150000E-02, 0.000000E+00
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    poldef
> OBJ,    TYPE,        USER_DEFINED
> OBJ,    PATCHES,     SCRSF
 
> OBJ,    NAME,        OUTL3
> 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,    TURBULENCE,  SAME , SAME
 
> OBJ,    NAME,        NWAL
> OBJ,    POSITION,    0.000000E+00, AT_END,       0.000000E+00
> OBJ,    SIZE,        TO_END,       0.000000E+00, TO_END
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    polcu10
> OBJ,    TYPE,        PLATE
STOP