```
TALK=T;RUN(1,1)
DISPLAY
A piston with a flat crown having an open cavity with a
central cusp moves inside a flat topped cylinder. The solution
domain is represented by a 2-part moving grid, the first
part of which expands and contracts in the Z-direction.
The second part of the grid moves with the piston, but does
not change its size. This part of the grid includes the
piston cavity, which is represented by porosities. Since the
model simulates only the compression and expansion parts of the
four-stroke cycle, there are no inflow or outflow boundaries.
The flow is represented in two-dimensional axi-symmetric
polar coordinates in the Y-Z plane but with swirl velocity
in the X-direction. An initial field of U-velocity is set up
to represent 'solid body' rotation of the fluid. The density
varies isentropically with pressure. This is based on the
PHOENICS-81 PDR/CHAM UK/6.
ENDDIS
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************

TEXT(Case 1 AXI-Symmetric Piston Bowl        )

************************************************************
************************************************************

IRUNN = 1 ;LIBREF = 14
************************************************************
Group 2. Time dependence
* Set overall time and no. of steps
TFIRST =1.157407E-04 ;TLAST =0.033681
FSTEP = 1 ;LSTEP = 10
Method of pairs used for grid setting.
TFRAC(1)=-7. ;TFRAC(2)=0.034483
TFRAC(3)=18. ;TFRAC(4)=0.017241
TFRAC(5)=40. ;TFRAC(6)=3.448276E-03
TFRAC(7)=4. ;TFRAC(8)=0.017241
TFRAC(9)=7. ;TFRAC(10)=0.034483
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = F
NX = 1
XULAST =0.1
XFRAC(1)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 11
YVLAST =0.0689
YFRAC(1)=0.073 ;YFRAC(2)=0.145
YFRAC(3)=0.218 ;YFRAC(4)=0.29
YFRAC(5)=0.363 ;YFRAC(6)=0.421
YFRAC(7)=0.458 ;YFRAC(8)=0.508
YFRAC(9)=0.653 ;YFRAC(10)=0.753
YFRAC(11)=0.88632
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 11
ZWLAST =0.1707
ZFRAC(1)=0.227 ;ZFRAC(2)=0.518
ZFRAC(3)=0.673 ;ZFRAC(4)=0.776
ZFRAC(5)=0.825 ;ZFRAC(6)=0.853
ZFRAC(7)=0.882 ;ZFRAC(8)=0.926
ZFRAC(9)=0.953 ;ZFRAC(10)=0.982
ZFRAC(11)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(3)=U1
NAME(5)=V1 ;NAME(7)=W1
NAME(149)=PRPS
* 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(U1,Y,Y,N,N,N,Y)
SOLUTN(V1,Y,Y,N,N,N,Y)
SOLUTN(W1,Y,Y,N,N,N,Y)
SOLUTN(PRPS,Y,N,N,N,N,N)
PRPS = 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(U1,Y,Y,Y,Y,Y,Y)
TERMS(V1,Y,Y,Y,Y,Y,Y)
TERMS(W1,Y,Y,Y,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWRH1 = T
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.03E+05
DVO1DT =0. ;DRH1DP = GRND3
RHO1A =1.4839E-04 ;RHO1B =0.765
RHO1C =0.
EMISS =0. ;SCATT =0.
ENUL =1.0E-05 ;ENUT =0.1
PRNDTL(U1)=1. ;PRNDTL(V1)=1.
PRNDTL(W1)=1.
PRT(U1)=1. ;PRT(V1)=1.
PRT(W1)=1.
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1.0E-10 ;FIINIT(U1)=407.130005
FIINIT(V1)=1.0E-10 ;FIINIT(W1)=1.0E-10
FIINIT(PRPS)=-1.

PATCH(CMP0 ,INIVAL, 1, 1, 8, 11, 6, 11, 1, 1)
INIT(CMP0 ,PRPS,0. ,198. )
FSWEEP = 1
IURINI = -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(HEAD ,LWALL , 1, 1, 1, 11, 1, 1, 1, 10)
COVAL(HEAD ,U1 , GRND2 ,0. )
COVAL(HEAD ,V1 , GRND2 ,0. )

PATCH(WALL ,NWALL , 1, 1, 11, 11, 1, 5, 1, 10)
COVAL(WALL ,U1 , GRND2 ,0. )
COVAL(WALL ,W1 , GRND2 ,0. )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 20 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-04
************************************************************
Group 16. Terminate Iterations
LITER(P1)=10 ;LITER(U1)=10
LITER(V1)=10 ;LITER(W1)=10
ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
ENDIT(V1)=1.0E-03 ;ENDIT(W1)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(U1,FALSDT,1.)
RELAX(V1,FALSDT,0.1)
RELAX(W1,FALSDT,0.1)
RELAX(PRPS,LINRLX,1.)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=5.0E+04
VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
PARSOL = F
IZW1 = 5
AZW1 =150.796005 ;BZW1 =0.06985
CZW1 =3.7265 ;DZW1 =0.
ISG62 = 1
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
SPEDAT(SET,MATERIAL,198,L,T)
************************************************************
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(U1,Y,N,Y,Y,Y,Y)
OUTPUT(V1,Y,N,Y,Y,Y,Y)
OUTPUT(W1,Y,N,Y,Y,Y,Y)
OUTPUT(PRPS,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 1 ;IZMON = 5
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = 1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 100000
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 = 2 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20

PATCH(PROF ,PROFIL, 1, 1, 1, 11, 3, 3, 10, 10)
PLOT(PROF ,U1 ,0. ,0. )
PLOT(PROF ,V1 ,0. ,0. )
PLOT(PROF ,W1 ,0. ,0. )

PATCH(CONT ,CONTUR, 1, 1, 1, 11, 1, 11, 10, 10)
PLOT(CONT ,U1 ,0. ,20. )
PLOT(CONT ,W1 ,0. ,20. )
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP
```