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
 ************************************************************
  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