PHOTON USE
  p;;;;;
 
  msg grid
  gr x 1
  msg Press  to continue
  pause;  cl
  msg 1st-phase (lighter-fluid) velocity vectors
  vec x 1 sh;  gr ou x 1
  msg Press  to continue
  pause;  cl
  msg 2nd-phase (dense fluid) velocity vectors
  set vec comp; - v2 w2;  vec x 1 sh;  gr ou x 1
  msg Press  to continue
  pause;  cl
  msg contour of 2nd-phase (heavier-fluid) volume fraction
  con r2 x 1 fi;0.01;  gr ou x 1
  msg Press  to continue
  pause;  cl
  msg contour of turbulence kinetic energy
  con ke x 1 fi;0.01;  gr ou x 1
  msg Press  to continue
  pause;  cl
  msg contour of reciprocal of turbulent time scale
  con epke x 1 fi;0.01;  gr ou x 1
  msg Press e to END
  enduse
    GROUP 1. Run title and other preliminaries
TEXT(ROUND JET;K-E 2PHS TURB MODULATION: W573
TITLE
  DISPLAY
    The case considered is a submerged turbulent round jet of air
    issuing into a low-velocity free stream. The jet is laden with
    solid particles, and the nozzle-to-stream velocity ratio is
    205. This flow was studied experimentally by Modaress et al
    [AIAA J, Vol.22, No.5, p624, 1984] and numerically by Hamill
    and Malin [PHOENICS J, Vol.4, Supp.II, p212, 1991]. The
    calculation may be made with the standard k-e model, or
    alternatively with the Mostafa-Mongia k-e variant which allows
    for gas turbulence modulation due to the presence of particles.
    The pipe Reynolds number is 1.33E4, the density ratio is 2538,
    the particle-to-air mass-flow ratio is 0.85, and the particle
    diameter is 50 microns.
  ENDDIS
REAL(TKEIN,EPSIN,GMIXL,MAS1,MAS2,R1IN,R2IN,AJET)
REAL(WJET,REYNO,DIAM,TKEF,EPSF,GRAD,PI,WFREE);CHAR(CTURB)
MESG( Enter the required turbulence-modulation option;
MESG( default NONE
MESG( The options are:
MESG(  MOST  - Mostafa-Mongia k-e modulation sources
MESG(  NONE  - Standard k-e model
MESG(
READVDU(CTURB,CHAR,NONE)
REYNO=1.33E4;DIAM=0.02;GRAD=0.5*DIAM;PI=3.14159;WFREE=0.05
AJET=PI*GRAD*GRAD;MAS1=3.76E-3/AJET;MAS2=3.2E-3/AJET
RHO1=1.178;RHO2=2990
R1IN=MAS1*RHO2/(MAS2*RHO1+MAS1*RHO2);R2IN=1-R1IN
WJET=MAS1/(RHO1*R1IN);TKEIN=WJET*WJET*0.001
    Inlet dissipation rate=.1643*k**1.5/Lm
GMIXL=0.1*GRAD;EPSIN=0.1643*TKEIN**1.5/GMIXL
    GROUP 3. X-direction grid specification
CARTES=F;XULAST=0.1
    GROUP 4. Y-direction grid specification
NREGY=2
IREGY=1;GRDPWR(Y,10,GRAD,1)
IREGY=2;GRDPWR(Y,10,DIAM*6,2)
    GROUP 5. Z-direction grid specification
GRDPWR(Z,20,DIAM*25,1.7)
    GROUP 7. Variables stored, solved & named
SOLVE(P1,V1,W1);SOLUTN(P1,P,P,Y,P,P,P)
ONEPHS=F;SOLVE(V2,W2,R1,R2)
    GROUP 9. Properties of the medium (or media)
ENUL=WJET*DIAM/REYNO
  ** Te, Tp are characteristic time scales of turbulence
TURMOD(KEMODL);STORE(ENUT,LEN1,GEN1,TE,TP)
PRT(V2)=1.E10;PRT(W2)=1.E10
    GROUP 10. Inter-phase-transfer processes and properties
CFIPS=GRND7;CFIPA=WJET/10000;CFIPB=50E-6
    Activate the following line if required
STORE(CFIP,REYN,CD,APRJ,VREL)
    GROUP 11. Initialization of variable or porosity fields
FIINIT(R1)=0.9999;FIINIT(R2)=0.0001
FIINIT(KE)=TKEIN;FIINIT(EP)=EPSIN
FIINIT(W1)=WJET;FIINIT(W2)=WJET
TKEF=(0.01*WFREE)**2;EPSF=0.09*TKEF*TKEF/ENUL
    GROUP 13. Boundary conditions and special sources
    Outer Boundary-- free stream
PATCH(HIGHY,NORTH,1,1,NY,NY,1,NZ,1,1)
COVAL(HIGHY,P1,1.0E3,0.0)
COVAL(HIGHY,W1,ONLYMS,0.0);COVAL(HIGHY,V1,ONLYMS,0.0)
COVAL(HIGHY,KE,ONLYMS,TKEF);COVAL(HIGHY,EP,ONLYMS,EPSF)
 
PATCH(LOWZ,LOW,1,1,NY/2+1,NY,1,1,1,1)
COVAL(LOWZ,P1,1.0E3,0.0)
COVAL(LOWZ,W1,ONLYMS,0.0);COVAL(LOWZ,V1,ONLYMS,0.0)
COVAL(LOWZ,KE,ONLYMS,TKEF);COVAL(LOWZ,EP,ONLYMS,EPSF)
 
    Inlet Boundary-- uniform flow
INLET(UNIFORM,LOW,1,1,1,NY/2,1,1,1,1)
VALUE(UNIFORM,P1,MAS1);VALUE(UNIFORM,W1,WJET)
VALUE(UNIFORM,P2,MAS2);VALUE(UNIFORM,W2,WJET)
VALUE(UNIFORM,KE,TKEIN);VALUE(UNIFORM,EP,EPSIN)
 
OUTLET(OUT,HIGH,1,NX,1,NY,NZ,NZ,1,1)
COVAL(OUT,P1,1.E3,0);COVAL(OUT,P2,1E3*RHO2/RHO1,0)
  ** turbulence-modulation sources
CASE :CTURB: OF
WHEN NONE,4
+ MESG(No turbulence modulation sources
+ RELAX(KE,LINRLX,0.5);RELAX(EP,LINRLX,1.0)
WHEN MOST,4
+ MESG(Mostafa-Mongia turbulence-modulation sources
+ PATCH(KEDISP,CELL,1,NX,1,NY,1,NZ,1,LSTEP);STORE(TE,TP)
+ COVAL(KEDISP,KE,GRND1,ZERO);COVAL(KEDISP,EP,GRND1,ZERO)
+ RELAX(KE,LINRLX,0.3);RELAX(EP,LINRLX,0.3)
ENDCASE
    GROUP 16. Termination of iterations
LSWEEP=300
    GROUP 17. Under-relaxation devices
RELAX(P1,LINRLX,.6);RELAX(r1,LINRLX,.6);RELAX(R2,LINRLX,.6)
RELAX(V1,FALSDT,ZWLAST/WJET);RELAX(W1,FALSDT,ZWLAST/WJET)
RELAX(V2,FALSDT,ZWLAST/WJET);RELAX(W2,FALSDT,ZWLAST/WJET)
KELIN=3
    GROUP 18. Limits on variables or increments to them
VARMIN(R1)=1E-10;VARMIN(R2)=1E-10
    GROUP 22. Monitor print-out
IZMON=NZ/2;IYMON=2;NPLT=10;TSTSWP=-1;ITABL=3
    GROUP 24. Dumps for restarts