```
TALK=T;RUN(1,1)
DISPLAY
The problem considered is the calculation of turbulent flow and
heat transfer in a sudden pipe expansion. The geometry consists of
an axisymmetric sudden pipe expansion of diameter ratio Do/Di=
1.945. The Reynolds number based upon the bulk velocity and
diameter of the larger pipe is 2.0E5. The thermal field is
calculated for a constant heat flux through the side wall and an
adiabatic back step. The laminar Prandtl number is taken as 1.0
and the turbulent Prandtl number as 0.86.

The calculation may be performed with one of the following low-Re
turbulence models: the Lam-Bremhorst k-e model; the Chen-Kim k-e
model; the 2-Layer k-e model; or the Wilcox k-omega model. The
numerical integration is taken right down to the wall, although
the mesh employed is still relatively coarse and so the solution
cannot be considered as grid independent.
ENDDIS
AUTOPLOT USE
file
phi 5

da 1 w1 y m
shift z -.972 1
divide z .243 1
col9 1
level y 0;level z 0
msg Velocity (W1) profile
msg Press e to END
ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************

TEXT(CHEN-KIM_2D ABRUPT PIPE EXPANS  :T213   )

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

IRUNN = 1 ;LIBREF = 14
************************************************************
Group 2. Time dependence
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = F
NX = 1
XULAST =0.1
XFRAC(1)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 20
YVLAST =0.5
YFRAC(1)=0.097686 ;YFRAC(2)=0.18509
YFRAC(3)=0.262211 ;YFRAC(4)=0.329049
YFRAC(5)=0.385604 ;YFRAC(6)=0.431877
YFRAC(7)=0.467866 ;YFRAC(8)=0.493573
YFRAC(9)=0.508997 ;YFRAC(10)=0.514139
YFRAC(11)=0.529887 ;YFRAC(12)=0.565305
YFRAC(13)=0.616078 ;YFRAC(14)=0.680379
YFRAC(15)=0.757069 ;YFRAC(16)=0.83376
YFRAC(17)=0.898061 ;YFRAC(18)=0.948834
YFRAC(19)=0.984252 ;YFRAC(20)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 25
ZWLAST =4.858612
ZFRAC(1)=0.053061 ;ZFRAC(2)=0.097959
ZFRAC(3)=0.134694 ;ZFRAC(4)=0.163265
ZFRAC(5)=0.183673 ;ZFRAC(6)=0.195918
ZFRAC(7)=0.2 ;ZFRAC(8)=0.20275
ZFRAC(9)=0.211 ;ZFRAC(10)=0.22475
ZFRAC(11)=0.244 ;ZFRAC(12)=0.26875
ZFRAC(13)=0.299 ;ZFRAC(14)=0.33475
ZFRAC(15)=0.376 ;ZFRAC(16)=0.42275
ZFRAC(17)=0.475 ;ZFRAC(18)=0.518296
ZFRAC(19)=0.574469 ;ZFRAC(20)=0.636807
ZFRAC(21)=0.703519 ;ZFRAC(22)=0.773686
ZFRAC(23)=0.846734 ;ZFRAC(24)=0.922269
ZFRAC(25)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(5)=V1
NAME(7)=W1 ;NAME(12)=KE
NAME(13)=EP ;NAME(14)=H1
NAME(145)=PRPS ;NAME(147)=REYN
NAME(148)=LTLS ;NAME(149)=WDIS
NAME(150)=ENUT
* 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(V1,Y,Y,N,N,N,N)
SOLUTN(W1,Y,Y,N,N,N,N)
SOLUTN(KE,Y,Y,N,N,N,N)
SOLUTN(EP,Y,Y,N,N,N,N)
SOLUTN(H1,Y,Y,Y,N,N,N)
SOLUTN(PRPS,Y,N,N,N,N,N)
SOLUTN(REYN,Y,N,N,N,N,Y)
SOLUTN(LTLS,Y,Y,Y,N,N,Y)
SOLUTN(WDIS,Y,N,N,N,N,N)
SOLUTN(ENUT,Y,N,N,N,N,Y)
VIST = 150
PRPS = 145
************************************************************
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(V1,Y,Y,Y,Y,Y,Y)
TERMS(W1,Y,Y,Y,Y,Y,Y)
TERMS(KE,N,Y,Y,Y,Y,N)
TERMS(EP,N,Y,Y,Y,Y,N)
TERMS(H1,N,Y,Y,Y,Y,Y)
TERMS(LTLS,N,N,Y,N,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWENT = 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 =1. ;TMP1 =0.
EL1 = GRND4
TSURR =0. ;TEMP0 =0.
PRESS0 =0.
DVO1DT =0. ;DRH1DP =0.
EMISS =0. ;SCATT =0.
EL1A =0. ;EL1B =0.
EL1C =0.
ENUL =1.321694E-06 ;ENUT = GRND3
ENUTA =0. ;ENUTB =0.
ENUTC =0.
IENUTA = 4
PRNDTL(V1)=1. ;PRNDTL(W1)=1.
PRNDTL(KE)=1. ;PRNDTL(EP)=1.
PRNDTL(H1)=1. ;PRNDTL(LTLS)=1.
PRT(V1)=1. ;PRT(W1)=1.
PRT(KE)=0.75 ;PRT(EP)=1.15
PRT(H1)=0.86 ;PRT(LTLS)=1.
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1.3E-04 ;FIINIT(V1)=0.
FIINIT(W1)=1. ;FIINIT(KE)=4.5E-03
FIINIT(EP)=1.378912E-03 ;FIINIT(H1)=1.0E-10
FIINIT(PRPS)=-1. ;FIINIT(REYN)=1.0E-10
FIINIT(LTLS)=1.0E-10 ;FIINIT(WDIS)=1.0E-02
FIINIT(ENUT)=1.0E-10

PATCH(CMP2 ,INIVAL, 1, 1, 11, 20, 1, 7, 1, 1)
INIT(CMP2 ,PRPS,0. ,198. )
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(KESOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1)
COVAL(KESOURCE,KE , GRND4 , GRND4 )
COVAL(KESOURCE,EP , GRND4 , GRND4 )

PATCH(KECHEN ,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1)
COVAL(KECHEN ,EP , FIXFLU , GRND4 )

PATCH(INLET ,LOW , 1, 1, 1, 10, 1, 1, 1, 1)
COVAL(INLET ,P1 , FIXFLU ,1. )
COVAL(INLET ,V1 ,0. ,0. )
COVAL(INLET ,W1 ,0. ,1. )
COVAL(INLET ,KE ,0. ,4.5E-03 )
COVAL(INLET ,EP ,0. ,2.268462E-03 )
COVAL(INLET ,H1 ,0. ,0. )

PATCH(OUTLET ,HIGH , 1, 1, 1, 20, 25, 25, 1, 1)
COVAL(OUTLET ,P1 ,1.0E+05 ,0. )
COVAL(OUTLET ,H1 ,0. , SAME )

PATCH(WFUNNORT,NWALL , 1, 1, 20, 20, 8, 25, 1, 1)
COVAL(WFUNNORT,W1 , GRND2 ,0. )
COVAL(WFUNNORT,KE ,1. ,0. )

PATCH(HEATIN ,NORTH , 1, 1, 20, 20, 8, 25, 1, 1)
COVAL(HEATIN ,H1 , FIXFLU ,0.017002 )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 250 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-04
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(V1)=10
LITER(W1)=10 ;LITER(KE)=20
LITER(EP)=20 ;LITER(H1)=20
LITER(LTLS)=20
ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03
ENDIT(W1)=1.0E-03 ;ENDIT(KE)=1.0E-03
ENDIT(EP)=1.0E-03 ;ENDIT(H1)=1.0E-03
ENDIT(LTLS)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(V1,FALSDT,0.194344)
RELAX(W1,FALSDT,0.194344)
RELAX(KE,FALSDT,0.194344)
RELAX(EP,FALSDT,0.194344)
RELAX(H1,FALSDT,1.)
RELAX(PRPS,LINRLX,1.)
RELAX(REYN,LINRLX,1.)
RELAX(LTLS,FALSDT,1.0E+09)
RELAX(WDIS,LINRLX,1.)
RELAX(ENUT,LINRLX,1.)
KELIN = 1
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
VARMAX(H1)=1.0E+09 ;VARMIN(H1)=0.
VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
VARMAX(REYN)=1.0E+10 ;VARMIN(REYN)=-1.0E+10
VARMAX(LTLS)=1.0E+10 ;VARMIN(LTLS)=-1.0E+10
VARMAX(WDIS)=1.0E+10 ;VARMIN(WDIS)=-1.0E+10
VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=1.0E-10
************************************************************
Group 19. Data transmitted to GROUND
GENK = T
PARSOL = F
IENUTA = 4
ISG62 = 1
SPEDAT(SET,MATERIAL,198,L,T)
************************************************************
Group 20. Preliminary Printout
************************************************************
Group 21. Print-out of Variables
INIFLD = F ;SUBWGR = F
* 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(V1,Y,N,Y,Y,Y,Y)
OUTPUT(W1,Y,N,Y,Y,Y,Y)
OUTPUT(KE,Y,N,Y,Y,Y,Y)
OUTPUT(EP,Y,N,Y,Y,Y,Y)
OUTPUT(H1,Y,N,Y,Y,Y,Y)
OUTPUT(PRPS,Y,N,Y,N,N,N)
OUTPUT(REYN,Y,N,Y,N,N,N)
OUTPUT(LTLS,Y,N,Y,Y,Y,Y)
OUTPUT(WDIS,Y,N,N,N,N,N)
OUTPUT(ENUT,Y,N,Y,N,N,N)
WALPRN = T
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 8 ;IZMON = 9
NPRMON = 10000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
XZPR = F ;YZPR = F
IPLTF = 1 ;IPLTL = 2000 ;NPLT = 10
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
No PATCHes yet used for this Group
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP
```