TALK=T;RUN(1,1)
  DISPLAY
  This is an application of the two-fluid model to
  1D turbulent combustion in reciprocating engines.

  Fluid 1 is defined as unburned combustible, and
  fluid 2 as fully burned combustion products.

  In a cylinder closed by a reciprocating piston, hot
  gas spreads from low z to high as a consequence entirely
  of the relative motion induced between hot and cold gas
  fragments. Cylinder walls and piston surface are impervious
  to heat and mass. Initially, the piston is at bottom dead
  centre, both fluids at rest, fluid 2 occupies first 10% of
  z-length, fluid 1 occupies the remainder.

  Inter-fluid friction is proportional to the square of the
  difference of velocities of the two fluids at each point.
  Inter-fluid mass transfer takes place, from fluid 1 to
  fluid 2, at a rate proportional to: r1*r2*velocity difference.
  Each fluid has density proportional to p**(1/gamma), where
  proportionality constant for fluid 2=0.3*that for fluid 1.
  ENDDIS
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(1D PISTON-IN-CYLINDER; 2-FLUID    : 978 )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 14
 ************************************************************
  Group 2. Time dependence
 STEADY = F
    * Set overall time and no. of steps
 TFIRST =0. ;TLAST =0.261799
 FSTEP = 1 ;LSTEP = 24
 TFRAC(1)=0.041667 ;TFRAC(4)=0.166667
 TFRAC(7)=0.291667 ;TFRAC(10)=0.416667
 TFRAC(13)=0.541667 ;TFRAC(16)=0.666667
 TFRAC(19)=0.791667 ;TFRAC(22)=0.916667
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 1
 XULAST =1.
 XFRAC(1)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 1
 YVLAST =1.
 YFRAC(1)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 20
 ZWLAST =0.10625
 ZFRAC(1)=0.05 ;ZFRAC(4)=0.2
 ZFRAC(7)=0.35 ;ZFRAC(10)=0.5
 ZFRAC(13)=0.65 ;ZFRAC(16)=0.8
 ZFRAC(19)=0.95
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = F
 EQUVEL = F
 NAME(1)=P1 ;NAME(7)=1W
 NAME(8)=2W ;NAME(9)=R1
 NAME(10)=R2 ;NAME(148)=DEN2
 NAME(149)=MDOT ;NAME(150)=DEN1
    * 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,N)
 SOLUTN(1W,Y,Y,N,Y,N,N)
 SOLUTN(2W,Y,Y,N,Y,N,N)
 SOLUTN(R2,Y,Y,N,N,N,Y)
 SOLUTN(DEN2,Y,N,N,N,N,Y)
 SOLUTN(MDOT,Y,N,N,N,N,Y)
 SOLUTN(DEN1,Y,N,N,N,N,Y)
 DEN1 = 150
 DEN2 = 148
 INTMDT = 149
 ************************************************************
  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(1W,Y,Y,Y,Y,Y,Y)
 TERMS(2W,Y,Y,Y,Y,N,Y)
 TERMS(R2,Y,Y,N,Y,N,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWRH1 = T
 NEWRH2 = T
 U1AD =0. ;V1AD =0.
 W1AD = ZMOVE
 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 = GRND3 ;TMP1 =0.
 EL1 =0.
 TSURR =0. ;TEMP0 =0.
 PRESS0 =1.0E+05
 DVO1DT =0. ;DRH1DP = GRND3
 RHO1A =2.682695E-04 ;RHO1B =0.714286
 RHO1C =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 RHO2 = GRND3 ;TMP2 =0.
 EL2 =0.
 DVO2DT =0. ;DRH2DP = GRND3
 RHO2A =8.048087E-05 ;RHO2B =0.714286
 RHO2C =0.
 ENUL =1.0E-05 ;ENUT =0.
 DRH2DP = GRND3
 PHINT(1W)=-2.022E+04 ;PHINT(2W)=-2.022E+04
 PHINT(R2)=-2.022E+04
 PRNDTL(1W)=1. ;PRNDTL(2W)=1.
 PRT(1W)=1. ;PRT(2W)=1.
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 CFIPS = GRND2
 RLOLIM =0. ;CMDOT =-1.
 CFIPA =0. ;CFIPB =0.
 CFIPC =100. ;CFIPD =0.
 CINT(P1)=1. ;CINT(1W)=1.
 CINT(2W)=1. ;CINT(R2)=1.
 PHINT(P1)=-2.022E+04 ;PHINT(1W)=-2.022E+04
 PHINT(2W)=-2.022E+04 ;PHINT(R2)=-2.022E+04
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.0E-10 ;FIINIT(1W)=1.0E-10
 FIINIT(2W)=1.0E-10 ;FIINIT(R1)=1.
 FIINIT(R2)=0. ;FIINIT(DEN2)=1.0E-10
 FIINIT(MDOT)=1.0E-10 ;FIINIT(DEN1)=1.0E-10
 
 PATCH(FLAME ,INIVAL, 1, 1, 1, 1, 1, 2, 1, 1)
 INIT(FLAME ,R1 ,0. ,-1. )
 INIT(FLAME ,R2 ,0. ,1. )
 INIADD = T
 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(SHSO ,CELL , 1, 1, 1, 1, 1, 1, 1, 1)
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 10 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 ISWR1 = 1 ;ISWR2 = 10000
 SELREF = T
 RESFAC =1.0E-04
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=-10 ;LITER(1W)=10
 LITER(2W)=10 ;LITER(R2)=1
 ENDIT(P1)=1.0E-03 ;ENDIT(1W)=1.0E-03
 ENDIT(2W)=1.0E-03 ;ENDIT(R2)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,0.5)
 RELAX(1W,FALSDT,2.083333E-05)
 RELAX(2W,FALSDT,2.083333E-05)
 RELAX(R1,LINRLX,1.)
 RELAX(R2,LINRLX,0.3)
 RELAX(DEN2,LINRLX,1.)
 RELAX(MDOT,LINRLX,1.)
 RELAX(DEN1,LINRLX,1.)
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
 VARMAX(1W)=1.0E+06 ;VARMIN(1W)=-1.0E+06
 VARMAX(2W)=1.0E+06 ;VARMIN(2W)=-1.0E+06
 VARMAX(R1)=1. ;VARMIN(R1)=1.0E-06
 VARMAX(R2)=0.999 ;VARMIN(R2)=1.0E-06
 VARMAX(DEN2)=1.0E+10 ;VARMIN(DEN2)=-1.0E+10
 VARMAX(MDOT)=1.0E+10 ;VARMIN(MDOT)=-1.0E+10
 VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 PARSOL = F
 IZW1 = 20
 AZW1 =24. ;BZW1 =0.05
 CZW1 =2. ;DZW1 =0.
 ISG62 = 1
 SPEDAT(SET,GXMONI,TRANSIENT,L,F)
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  Group 21. Print-out of Variables
 INIFLD = F ;SUBWGR = T
    * 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(1W,Y,N,Y,Y,Y,Y)
 OUTPUT(2W,Y,N,Y,Y,Y,Y)
 OUTPUT(R1,N,N,N,N,N,N)
 OUTPUT(R2,Y,N,Y,Y,Y,Y)
 OUTPUT(DEN2,N,N,N,N,N,N)
 OUTPUT(MDOT,Y,N,Y,N,N,N)
 OUTPUT(DEN1,N,N,N,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 1 ;IYMON = 1 ;IZMON = 10
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = 2
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 10 ;NUMCLS = 5
 NTPRIN = 3 ;ISTPRF = 1 ;ISTPRL = 100000
 NZPRIN = 3 ;IZPRF = 1 ;IZPRL = 10000
 XZPR = F ;YZPR = F
 IPLTF = 2 ;IPLTL = -1 ;NPLT = 1
 ISWPRF = 1 ;ISWPRL = 100000
 ITABL = 3 ;IPROF = 1
 ABSIZ =0.5 ;ORSIZ =0.2
 NTZPRF = 1 ;NCOLPF = 50
 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
 
 PATCH(NZD1 ,PROFIL, 1, 1, 1, 1, 10, 10, 1, 24)
 PLOT(NZD1 ,P1 ,0. ,1.0E+07 )
 
 PATCH(NZD2 ,PROFIL, 1, 1, 1, 1, 10, 10, 1, 24)
 PLOT(NZD2 ,1W ,-1. ,-1. )
 PLOT(NZD2 ,2W ,-1. ,-1. )
 PLOT(NZD2 ,R2 ,0. ,0. )
 PLOT(NZD2 ,MDOT,0. ,0. )
 
 PATCH(LONGPLOT,PROFIL, 1, 1, 1, 1, 1, 20, 1, 24)
 PLOT(LONGPLOT,1W ,-1. ,-1. )
 PLOT(LONGPLOT,2W ,-1. ,-1. )
 PLOT(LONGPLOT,MDOT,0. ,0. )
 
 PATCH(LONPL ,PROFIL, 1, 1, 1, 1, 1, 20, 1, 24)
 PLOT(LONPL ,R2 ,0. ,1. )
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP