```
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
* 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
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.
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. )
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
```