TALK=T;RUN(1,1)
PHOTON USE
p
text
1
1
Particles with different
0.18841E+04 0.23725E+04 CR
text
1
1
size and init. velocity
0.18841E+04 0.22422E+04 CR
text
1
1
are introduced into a radial
0.18841E+04 0.21119E+04 CR
text
1
1
impellar in a rotating BFC
0.18841E+04 0.19817E+04 CR
text
1
1
coordinate system.
0.18841E+04 0.18514E+04 CR
upause 4
gr out z 1
gr out z m
gr out y 1
gr out y m
upause 2
GEOMETRY READ
T00001
GEOMETRY READ
T00002
GEOMETRY READ
T00003
GEOMETRY READ
T00004
GEOMETRY READ
T00005
GEOMETRY READ
T00006
GEOMETRY READ
T00007
GEOMETRY READ
T00008
GEOMETRY READ
T00009
GEOMETRY READ
T00010
GEOMETRY READ
T00011
upause 5
vec y 3
upause 2
text
3
Press e to end
0.18841E+04 0.14265E+04 CR
enduse
DISPLAY
The continuous-phase set-up is that of PHOENICS Library Case
524.The geometry is that of the B-type impeller tested by Mizuki
et al. (1974). The solution domain represents the
passage between two blades and a 'vaneless space' beyond the outer
radius of the impeller, which is represented by applying cyclic
boundary conditions on the East and West surfaces at IZ=11 and 12.
The grid given in file GRID1 was generated by linear
interpolation from boundary values using a conical coordinate
system, and is non-orthogonal.
Particles of several sizes are introduced at the inlet.
Since the coordinate system is a rotating one, Coriolis and
centrifugal forces acting on the particles are automatically
introduced by GENTRA-Earth.
ENDDIS
30
*__XP __YP__ZP__UP__VP__WP__DIAM _DENSTY_FRATE
0.057 0.015 0.0 0.0 -1.0 10.0 7.0E-05 1000.0 1.E-06
0.057 0.015 0.0 0.0 -2.0 10.0 6.5E-05 1000.0 1.E-06
0.057 0.015 0.0 0.0 -3.0 10.0 6.0E-05 1000.0 1.E-06
0.057 0.015 0.0 0.0 -4.0 10.0 5.5E-05 1000.0 1.E-06
0.057 0.015 0.0 0.0 -5.0 10.0 5.0E-05 1000.0 1.E-06
0.057 0.015 0.0 0.0 -6.0 10.0 4.5E-05 1000.0 1.E-06
0.057 0.015 0.0 0.0 -7.0 5.0 4.0E-05 1000.0 1.E-06
0.057 0.015 0.0 0.0 -8.0 5. 3.5E-05 1000.0 1.E-06
0.057 0.015 0.0 0.0 -9.0 5.0 3.0E-05 1000.0 1.E-06
0.057 0.015 0.0 0.0 -10.0 5.0 2.5E-05 1000.0 1.E-06
0.057 0.015 0.0 0.0 -11.0 5.0 2.0E-05 1000.0 1.E-06
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(Particles in radial impeller: G204 )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 1
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 5
XULAST =1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 5
YVLAST =1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 12
ZWLAST =1.
************************************************************
Group 6. Body-Fitted Coordinates
BFC = T ;NONORT = T
NCRT = 1
RSTGEO = F ;SAVGEO = F
UUP = F ;VUP = F ;WUP = F
NGEOM =CHAM
NAMXYZ =CHAM
ANGMIN =20.
DOMAIN( 1, 6, 1, 6, 1, 13)
* Set fixed sub-domain
FIXDOM( 1, 0, 0, 0, 0, 0, 0)
FIXDOM( 2, 0, 0, 0, 0, 0, 0)
FIXDOM( 3, 0, 0, 0, 0, 0, 0)
FIXDOM( 4, 0, 0, 0, 0, 0, 0)
FIXDOM( 5, 0, 0, 0, 0, 0, 0)
FIXDOM( 6, 0, 0, 0, 0, 0, 0)
FIXDOM( 7, 0, 0, 0, 0, 0, 0)
FIXDOM( 8, 0, 0, 0, 0, 0, 0)
FIXDOM( 9, 0, 0, 0, 0, 0, 0)
FIXDOM( 10, 0, 0, 0, 0, 0, 0)
LIJ = F ;LJK = T ;LIK = F
* X-cyclic boundaries switched
XCYIZ( 11, T ) ; XCYIZ( 12, T )
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(3)=U1
NAME(5)=V1 ;NAME(7)=W1
NAME(145)=MOMZ ;NAME(146)=MOMY
NAME(147)=MOMX ;NAME(148)=WCRT
NAME(149)=VCRT ;NAME(150)=UCRT
* 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(MOMZ,Y,N,N,N,N,Y)
SOLUTN(MOMY,Y,N,N,N,N,Y)
SOLUTN(MOMX,Y,N,N,N,N,Y)
SOLUTN(WCRT,Y,N,N,N,N,N)
SOLUTN(VCRT,Y,N,N,N,N,N)
SOLUTN(UCRT,Y,N,N,N,N,N)
************************************************************
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
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.24 ;TMP1 =0. ;EL1 =0.
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =0. ;DRH1DP =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
ENUL =1.451613E-05 ;ENUT =1.032258E-03
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)=1.0E-10
FIINIT(V1)=1.0E-10 ;FIINIT(W1)=33.929176
FIINIT(MOMZ)=1.0E-10 ;FIINIT(MOMY)=1.0E-10
FIINIT(MOMX)=1.0E-10 ;FIINIT(WCRT)=1.0E-10
FIINIT(VCRT)=1.0E-10 ;FIINIT(UCRT)=1.0E-10
No PATCHes yet used for this Group
INIADD = F
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(BFC1 ,LOW , 1, 5, 5, 5, 1, 1, 1, 1)
COVAL(BFC1 ,P1 , FIXFLU ,42.072178 )
COVAL(BFC1 ,U1 ,0. , GRND1 )
COVAL(BFC1 ,W1 ,0. , GRND1 )
COVAL(BFC1 ,WCRT,0. ,33.929176 )
COVAL(BFC1 ,VCRT,0. ,-27.83449 )
COVAL(BFC1 ,UCRT,0. ,0. )
PATCH(BFC2 ,LOW , 1, 5, 4, 4, 1, 1, 1, 1)
COVAL(BFC2 ,P1 , FIXFLU ,42.072178 )
COVAL(BFC2 ,U1 ,0. , GRND1 )
COVAL(BFC2 ,W1 ,0. , GRND1 )
COVAL(BFC2 ,WCRT,0. ,33.929176 )
COVAL(BFC2 ,VCRT,0. ,-32.609707 )
COVAL(BFC2 ,UCRT,0. ,0. )
PATCH(BFC3 ,LOW , 1, 5, 3, 3, 1, 1, 1, 1)
COVAL(BFC3 ,P1 , FIXFLU ,42.072178 )
COVAL(BFC3 ,U1 ,0. , GRND1 )
COVAL(BFC3 ,W1 ,0. , GRND1 )
COVAL(BFC3 ,WCRT,0. ,33.929176 )
COVAL(BFC3 ,VCRT,0. ,-37.384926 )
COVAL(BFC3 ,UCRT,0. ,0. )
PATCH(BFC4 ,LOW , 1, 5, 2, 2, 1, 1, 1, 1)
COVAL(BFC4 ,P1 , FIXFLU ,42.072178 )
COVAL(BFC4 ,U1 ,0. , GRND1 )
COVAL(BFC4 ,W1 ,0. , GRND1 )
COVAL(BFC4 ,WCRT,0. ,33.929176 )
COVAL(BFC4 ,VCRT,0. ,-42.160145 )
COVAL(BFC4 ,UCRT,0. ,0. )
PATCH(BFC5 ,LOW , 1, 5, 1, 1, 1, 1, 1, 1)
COVAL(BFC5 ,P1 , FIXFLU ,42.072178 )
COVAL(BFC5 ,U1 ,0. , GRND1 )
COVAL(BFC5 ,W1 ,0. , GRND1 )
COVAL(BFC5 ,WCRT,0. ,33.929176 )
COVAL(BFC5 ,VCRT,0. ,-46.935356 )
COVAL(BFC5 ,UCRT,0. ,0. )
PATCH(GXIN ,LOW , 1, 5, 1, 5, 1, 1, 1, 1)
PATCH(GXOUT ,HIGH , 1, 5, 1, 5, 12, 12, 1, 1)
COVAL(GXOUT ,P1 ,100. ,0. )
COVAL(GXOUT ,U1 ,0. ,0. )
COVAL(GXOUT ,V1 ,0. ,0. )
COVAL(GXOUT ,W1 ,0. ,0. )
PATCH(ROTA ,PHASEM, 1, 5, 1, 5, 1, 12, 1, 1)
COVAL(ROTA ,U1 , FIXFLU , GRND1 )
COVAL(ROTA ,V1 , FIXFLU , GRND1 )
COVAL(ROTA ,W1 , FIXFLU , GRND1 )
PATCH(WFUN1 ,SWALL , 1, 5, 1, 1, 1, 12, 1, 1)
COVAL(WFUN1 ,U1 , GRND2 ,-67.858353 )
COVAL(WFUN1 ,W1 , GRND2 ,0. )
PATCH(WFUN2 ,NWALL , 1, 5, 5, 5, 1, 12, 1, 1)
COVAL(WFUN2 ,U1 , GRND2 ,0. )
COVAL(WFUN2 ,W1 , GRND2 ,0. )
PATCH(WFUN3 ,WWALL , 1, 1, 1, 5, 1, 10, 1, 1)
COVAL(WFUN3 ,V1 , GRND2 ,0. )
COVAL(WFUN3 ,W1 , GRND2 ,0. )
PATCH(WFUN4 ,EWALL , 5, 5, 1, 5, 1, 10, 1, 1)
COVAL(WFUN4 ,V1 , GRND2 ,0. )
COVAL(WFUN4 ,W1 , GRND2 ,0. )
PATCH(GENPAT ,CELL , 1, 5, 1, 5, 1, 12, 1, 1)
COVAL(GENPAT ,U1 , FIXFLU , GRND )
COVAL(GENPAT ,V1 , FIXFLU , GRND )
COVAL(GENPAT ,W1 , FIXFLU , GRND )
XCYCLE = T
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 15 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
LITER(P1)=10 ;LITER(U1)=1
LITER(V1)=1 ;LITER(W1)=1
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,0.4)
RELAX(U1,FALSDT,5.0E-04)
RELAX(V1,FALSDT,5.0E-04)
RELAX(W1,FALSDT,5.0E-04)
RELAX(MOMZ,LINRLX,0.7)
RELAX(MOMY,LINRLX,0.7)
RELAX(MOMX,LINRLX,0.7)
RELAX(WCRT,LINRLX,1.)
RELAX(VCRT,LINRLX,1.)
RELAX(UCRT,LINRLX,1.)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
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(MOMZ)=1.0E+10 ;VARMIN(MOMZ)=-1.0E+10
VARMAX(MOMY)=1.0E+10 ;VARMIN(MOMY)=-1.0E+10
VARMAX(MOMX)=1.0E+10 ;VARMIN(MOMX)=-1.0E+10
VARMAX(WCRT)=1.0E+10 ;VARMIN(WCRT)=-1.0E+10
VARMAX(VCRT)=1.0E+10 ;VARMIN(VCRT)=-1.0E+10
VARMAX(UCRT)=1.0E+10 ;VARMIN(UCRT)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
PARSOL = F
ISG62 = 1
ANGVEL =628.318054
ROTAZB =-1.
SPEDAT(SET,GENTRA,RESTART,C,NONE)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
IG( 4) = 2
IG( 9) = 10
IG( 10) = 1
IG( 11) = 20
IG( 12) = 1
IG( 14) = 1
IG( 15) = 1000
IG( 16) = 5
IG( 17) = 1
IG( 18) = 15
IG( 19) = 3
IG( 20) = 30
RG( 51) =0.0731
RG( 52) =942.
RG( 53) =8.65E-04
RG( 54) =1.76
RG( 55) =-0.206
RG( 56) =0.9
RG( 57) =1.1E+17
RG( 58) =0.495
RG( 71) = GRND1
RG( 72) =4131.799805
RG( 74) = GRND1
RG( 75) = GRND1
RG( 76) = GRND1
RG( 77) =1.
RG( 80) =1.
RG( 82) =4131.799805
RG( 83) =1.0E+04
RG( 84) =1.
RG( 86) =18.
RG( 87) =28.9
RG( 88) = GRND1
RG( 89) =1007.
RG( 90) = GRND1
RG( 91) =273.149994
RG( 92) =273.149994
RG( 93) =1.
RG( 94) =3.335E+05
RG( 95) = GRND1
RG( 96) =0.0263
RG( 97) = GRND1
CG( 6) ='GHIS'
CG( 8) ='Q1'
CG( 9) ='NONE'
CG( 10) ='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,Y,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(MOMZ,Y,N,Y,N,N,N)
OUTPUT(MOMY,Y,N,Y,N,N,N)
OUTPUT(MOMX,Y,N,Y,N,N,N)
OUTPUT(WCRT,Y,N,Y,N,N,N)
OUTPUT(VCRT,Y,N,Y,N,N,N)
OUTPUT(UCRT,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 2 ;IYMON = 2 ;IZMON = 6
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = 1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 15 ;NUMCLS = 5
NXPRIN = 2 ;IXPRF = 1 ;IXPRL = 10000
NYPRIN = 2 ;IYPRF = 1 ;IYPRL = 10000
NZPRIN = 2 ;IZPRF = 1 ;IZPRL = 10000
XZPR = F ;YZPR = F
IPLTF = 1 ;IPLTL = -1 ;NPLT = 1
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
PATCH(SUCTION ,CONTUR, 5, 5, 1, 5, 1, 12, 1, 1)
PLOT(SUCTION ,P1 ,0. ,20. )
PATCH(PRESSURE,CONTUR, 1, 1, 1, 5, 1, 12, 1, 1)
PLOT(PRESSURE,P1 ,0. ,20. )
PATCH(MIDWAY ,CONTUR, 1, 5, 2, 2, 1, 12, 1, 1)
PLOT(MIDWAY ,P1 ,0. ,20. )
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP