TALK=T;RUN(1,1)
DISPLAY
F Ricou (London University PhD thesis 1959) injected gases into
a porous-walled cylindrical chamber with a circular exit opposite
the injection orifice. He varied the inflow rate of air through
the porous wall until the pressure in the chamber equalled that
outside, and took this to be the rate of entrainment of air into
jet.
When the flow was turbulent, he found that the rate of entrain-
ment per unit length equalled 0.32 times the mass-injection rate
divided by the nozzle diameter, if the injected fluid had the
same density as the surrounding fluid.
Please press RETURN to continue
When the densities differed he found that the 0.32 factor should
be replaced by:
0.32 * (surrounding-fluid density/injected-fluid density)**0.5
This Q1 file causes PHOENICS to simulate the conditions of the
experiment. The Reynolds number and density ratio can be varied.
The Prandtl mixing-length model of turbulence is used, with the
low-Reynolds-Number modification activated by IENUTA=6 .
PHOTON commands are supplied.
ENDDIS
PHOTON USE
p
phi
up z
gr ou x 1
msg contours of axial velocity
con w1 x 1 fi;0.001;pause;con off; red
msg contours of concentration of injected material
con a x 1 fi;0.001;pause;con off;red
msg enut contours
con enut x 1 fi;0.001;pause;con off;red
msg velocity vectors
vec x 1 sh;pause;con off;vec off;red
msg streamlines
stre 2d x 1 sh;int 20
enduse
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(Ricou's Jet In Cylindrical Chamber )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 281
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = F
NX = 1
XULAST =1.0E-02
XFRAC(1)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 40
YVLAST =0.1143
YFRAC(1)=0.025 ;YFRAC(5)=0.125
YFRAC(9)=0.225 ;YFRAC(13)=0.325
YFRAC(17)=0.425 ;YFRAC(21)=0.525
YFRAC(25)=0.625 ;YFRAC(29)=0.725
YFRAC(33)=0.825 ;YFRAC(37)=0.925
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 40
ZWLAST =0.32385
ZFRAC(1)=0.025 ;ZFRAC(5)=0.125
ZFRAC(9)=0.225 ;ZFRAC(13)=0.325
ZFRAC(17)=0.425 ;ZFRAC(21)=0.525
ZFRAC(25)=0.625 ;ZFRAC(29)=0.725
ZFRAC(33)=0.825 ;ZFRAC(37)=0.925
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(5)=V1
NAME(7)=W1 ;NAME(16)=A
NAME(149)=LEN1 ;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,N,N,N,Y)
SOLUTN(V1,Y,Y,N,N,N,Y)
SOLUTN(W1,Y,Y,N,N,N,Y)
SOLUTN(A,Y,Y,N,N,N,Y)
SOLUTN(LEN1,Y,N,N,N,N,Y)
SOLUTN(ENUT,Y,N,N,N,N,Y)
VIST = 150
LEN1 = 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(V1,Y,Y,Y,Y,Y,Y)
TERMS(W1,Y,Y,Y,Y,Y,Y)
TERMS(A,N,Y,Y,Y,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 = GRND6
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =0. ;DRH1DP =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
EL1A =0.05 ;EL1B =0.075 ;EL1C =1.0E-02
EL1D =0. ;EL1E =50.
ENUL =1.0E-05 ;ENUT = GRND2
ENUTA =0.2 ;ENUTB =4. ;ENUTC =0.
IENUTA = 6
PRNDTL(V1)=1. ;PRNDTL(W1)=1.
PRNDTL(A)=0.7
PRT(V1)=1. ;PRT(W1)=1.
PRT(A)=0.7
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1.0E-10 ;FIINIT(V1)=1.0E-10
FIINIT(W1)=1.0E-10 ;FIINIT(A)=1.0E-10
FIINIT(LEN1)=1.0E-10 ;FIINIT(ENUT)=0.018288
PATCH(INITIAL ,INIVAL, 1, 1, 1, 8, 1, 40, 1, 1)
INIT(INITIAL ,W1 ,0. ,50. )
INIT(INITIAL ,A ,0. ,1. )
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(IN ,LOW , 1, 1, 1, 8, 1, 1, 1, 1)
COVAL(IN ,P1 , FIXFLU ,50. )
COVAL(IN ,V1 ,0. ,0. )
COVAL(IN ,W1 ,0. ,50. )
COVAL(IN ,A ,0. ,1. )
PATCH(WALL ,NORTH , 1, 1, 40, 40, 1, 40, 1, 1)
COVAL(WALL ,P1 , FIXFLU ,0.4 )
COVAL(WALL ,W1 ,0. ,0. )
COVAL(WALL ,A ,0. ,0. )
PATCH(OUTLET ,HIGH , 1, 1, 1, 32, 40, 40, 1, 1)
COVAL(OUTLET ,P1 ,1.0E+05 ,0. )
COVAL(OUTLET ,V1 ,0. ,0. )
COVAL(OUTLET ,W1 ,0. ,0. )
COVAL(OUTLET ,A ,0. , SAME )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 100 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-02
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(V1)=10
LITER(W1)=10 ;LITER(A)=20
ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03
ENDIT(W1)=1.0E-03 ;ENDIT(A)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(V1,FALSDT,0.1)
RELAX(W1,FALSDT,1.0E-03)
RELAX(A,FALSDT,1.0E+09)
RELAX(LEN1,LINRLX,1.)
RELAX(ENUT,LINRLX,0.5)
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(A)=1.0E+10 ;VARMIN(A)=-1.0E+10
VARMAX(LEN1)=1.0E+10 ;VARMIN(LEN1)=-1.0E+10
VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
DWDY = T
PARSOL = F
IENUTA = 6
ISG62 = 1
SPEDAT(SET,GXMONI,PLOTALL,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(A,Y,N,Y,Y,Y,Y)
OUTPUT(LEN1,Y,N,Y,N,N,N)
OUTPUT(ENUT,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 19 ;IZMON = 10
NPRMON = 20 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NYPRIN = 4 ;IYPRF = 1 ;IYPRL = 10000
NZPRIN = 4 ;IZPRF = 1 ;IZPRL = 10000
XZPR = F ;YZPR = F
IPLTF = 1 ;IPLTL = -1 ;NPLT = 1
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 1 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
PATCH(EXIT ,PROFIL, 1, 1, 1, 40, 39, 39, 1, 1)
PLOT(EXIT ,W1 ,0. ,0. )
PLOT(EXIT ,A ,0. ,0. )
PATCH(PIPE ,CONTUR, 1, 1, 1, 40, 1, 40, 1, 1)
PLOT(PIPE ,W1 ,0. ,10. )
PLOT(PIPE ,A ,0. ,10. )
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP