PHOTON USE
  p
 
 
 
  gr ou z 1
  ve z 1 sh
  msg                    Velocity distribution:
  MSG Press  to continue
  pause
  clear
  msg                    Fuel distribution :
  con fuel z 1 fil;0.01
  msg Press e to END
  ENDUSE
 
   *  GROUP 1. Run title and other preliminaries.
 
  DISPLAY
      This case describes how the SCRS option in conjuction
      with PLANT can be used to predict the turbulent diffusion
      of three gases with different physical properties and
      eddy break-up model to represent the chemical mass source.
  ENDDIS
 
TEXT(Three component mixing of different gases:133
REAL(M1IN,M2IN)
REAL(HIN,WMU1,WMU2,WMU3,CPF,CPS,CPT)
REAL(UIN1,UIN2,UIN3)
UIN1=15.;UIN2=5.;UIN3=15.
HIN=7.5e05
M1IN=0.06;M2IN=0.02
    AIR      NH3      NO2
WMU1=29.;WMU2=17.;WMU3=46.
CPF=1500.;CPS=1500.;CPT=1500.
CP1=GRND10
CP1A=CPF; CP1B=CPS;CP1C=CPT
   *  GROUP 3. X-direction grid specification.
CARTES=F
XULAST= 1.570E+00
NX=20
GRDPWR(X,NX,XULAST, 1.000E+00)
INTEGER(NXF01,NXL01); NXF01=1; NXL01=NX
   *  GROUP 4. Y-direction grid specification.
RINNER= 2.500E+00
NY=15
INTEGER(NYF01,NYL01); NYF01=1; NYL01=5
YFRAC(1)= 5.000E-02;YFRAC(2)= 1.000E-01
YFRAC(3)= 1.500E-01;YFRAC(4)= 2.000E-01
YFRAC(5)= 2.500E-01
   *  SUBGRD(Y,6,10, 2.500E-01, 1.000E+00)
INTEGER(NYF02,NYL02); NYF02=6; NYL02=10
YFRAC(6)= 3.000E-01;YFRAC(7)= 3.500E-01
YFRAC(8)= 4.000E-01;YFRAC(9)= 4.500E-01
YFRAC(10)= 5.000E-01
   *  SUBGRD(Y,11,15, 2.500E-01, 1.000E+00)
INTEGER(NYF03,NYL03); NYF03=11; NYL03=15
YFRAC(11)= 5.500E-01;YFRAC(12)= 6.000E-01
YFRAC(13)= 6.500E-01;YFRAC(14)= 7.000E-01
YFRAC(15)= 7.500E-01
   *  GROUP 7. Variables stored, solved & named.
        OXID =  M1
        FUEL =  M2
        PROD =  M3
        MIXF =  M2+M3,
        i.e.
        OXID =  1-MIXF
        PROD =  1-FUEL-OXID
SOLVE(MIXF,FUEL)
SOLVE(P1); SOLUTN(P1,Y,Y,Y,N,N,N)
SOLVE(U1,V1)
STORE(OXID,PROD,H1,RHO1,TMP1)
SOLVE(KE,EP)
ENUT=PRKOLM;EL1=KE15DEP;KELIN=0
PATCH(KESOURCE,PHASEM,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(KESOURCE,KE,GRND4,GRND4)
COVAL(KESOURCE,EP,GRND4,GRND4)
GENK=T
   *  GROUP 9. Properties of the medium (or media).
REAL(RHOR);RHOR=0.75
RHO1= 3GASES;RHO1B=WMU1;RHO1C=WMU3;RHO1A=WMU2
TMP1=SCRSNONEQ;TMP1B=CPT;TMP1C=CPF;TMP1A=CPS
TMP2A=1.-1.E-6;TMP2B=0.0
CP1A=CPS; CP1B=CPT; CP1C=CPF
PRESS0=1.e05
ENUL= 1.000E-05
PRNDTL(MIXF)=1.
   *  GROUP 11. Initialization of variable or porosity fields.
FIINIT(U1)= -(UIN1+UIN2+UIN3)/3.
FIINIT(MIXF)= 1.-M1IN
FIINIT(FUEL)= M2IN
FIINIT(H1)=HIN
FIINIT(RHO1)=RHOR
FIINIT(OXID)=M1IN
FIINIT(PROD)=1.-M1IN-M2IN
FIINIT(KE)= 1.000E+01
FIINIT(EP)= 2.000E+01
   *  GROUP 13. Boundary conditions and special sources.
PATCH(INL1,EAST,NXL01,NXL01,NYF01,NYL01,1,1,1,1)
COVAL(INL1,P1,FIXFLU,RHOR* UIN1)
COVAL(INL1,U1,ONLYMS,   -UIN1)
COVAL(INL1,V1,ONLYMS,   0.000E+00)
COVAL(INL1,MIXF,ONLYMS, 1.e-5)
COVAL(INL1,FUEL,ONLYMS, 1.e-5)
COVAL(INL1,KE,ONLYMS, 1.000E+01)
COVAL(INL1,EP,ONLYMS, 2.000E+01)
COVAL(INL1,H1,ONLYMS, 1.000E+05)
 
PATCH(INL2,EAST,NXL01,NXL01,NYF02,NYL02,1,1,1,1)
COVAL(INL2,P1,FIXFLU,RHOR* UIN2)
COVAL(INL2,U1,ONLYMS,   -UIN2)
COVAL(INL2,V1,ONLYMS,   0.000E+00)
COVAL(INL2,MIXF,ONLYMS, 0.05)
COVAL(INL2,FUEL,ONLYMS, 0.025)
COVAL(INL2,KE,ONLYMS, 2.000E+01)
COVAL(INL2,EP,ONLYMS, 4.000E+01)
 
PATCH(INL3,EAST,NXL01,NXL01,NYF03,NYL03,1,1,1,1)
COVAL(INL3,P1,FIXFLU,RHOR* UIN3)
COVAL(INL3,U1,ONLYMS,   -UIN3)
COVAL(INL3,V1,ONLYMS,   0.000E+00)
COVAL(INL3,MIXF,ONLYMS, 1.e-5)
COVAL(INL3,FUEL,ONLYMS, 1.e-5)
COVAL(INL3,KE,ONLYMS, 1.000E+01)
COVAL(INL3,EP,ONLYMS, 2.000E+01)
 
PATCH(OUT,WEST,NXF01,NXF01,NYF01,NYL03,1,1,1,1)
COVAL(OUT,P1,FIXP, 0.000E+00)
COVAL(OUT,MIXF,ONLYMS,SAME)
COVAL(OUT,FUEL,ONLYMS,SAME)
COVAL(OUT,KE,ONLYMS,SAME)
COVAL(OUT,EP,ONLYMS,SAME)
 
PATCH(WALLS,SWALL,NXF01,NXL01,NYF01,NYF01,1,1,1,1)
COVAL(WALLS,U1,LOGLAW, 0.000E+00)
COVAL(WALLS,KE,LOGLAW,LOGLAW)
COVAL(WALLS,EP,LOGLAW,LOGLAW)
 
PATCH(WALLN,NWALL,NXF01,NXL01,NYL03,NYL03,1,1,1,1)
COVAL(WALLN,U1,LOGLAW, 0.000E+00)
COVAL(WALLN,KE,LOGLAW,LOGLAW)
COVAL(WALLN,EP,LOGLAW,LOGLAW)
    
     PLANTBEGIN
   *  Chemical reaction sources
PATCH(REAC,PHASEM,1,NX,1,NY,1,NZ,1,LSTEP)
    COF=4.0*EP/(KE+1.e-10)
COVAL(REAC,FUEL,GRND,0.0)
     PLANTEND

   *  GROUP 15. Termination of sweeps.
LSWEEP=100
RESREF(P1)= 1.000E-06
RESREF(U1)= 1.000E-06
RESREF(V1)= 1.000E-06
RESREF(MIXF)= 1.000E-06
RESREF(FUEL)= 1.000E-06
RESREF(KE)= 1.000E-06
RESREF(EP)= 1.000E-06
   *  GROUP 16. Termination of iterations.
LITER(U1)=1; ENDIT(U1)= 1.000E-06
LITER(V1)=1; ENDIT(V1)= 1.000E-06
LITER(MIXF)=1; ENDIT(MIXF)= 1.000E-06
LITER(FUEL)=1; ENDIT(FUEL)= 1.000E-06
LITER(KE)=1; ENDIT(KE)= 1.000E-06
LITER(EP)=1; ENDIT(EP)= 1.000E-06
   *  GROUP 17. Under-relaxation devices.
RELAX(P1,LINRLX,0.5)
real(fac);fac=10.
REAL(SCALEL,SCALEU);SCALEL= 1.000E-01;SCALEU= 30.000E+00
RELAX(U1,FALSDT, fac*SCALEL/SCALEU)
RELAX(V1,FALSDT, fac*SCALEL/SCALEU)
RELAX(KE,FALSDT, fac*SCALEL/SCALEU)
RELAX(EP,FALSDT, fac*SCALEL/SCALEU)
RELAX(FUEL,FALSDT,0.01)
   *  GROUP 19. Data communicated by satellite to GROUND.
NAMSAT=MOSG
   *  GROUP 20. Preliminary print-out.
   *  GROUP 21. Print-out of variables.
OUTPUT(P1,Y,Y,N,Y,Y,Y)
OUTPUT(U1,Y,Y,N,Y,Y,Y)
OUTPUT(V1,Y,Y,N,Y,Y,Y)
OUTPUT(MIXF,Y,Y,N,Y,Y,Y)
OUTPUT(PROD,Y,Y,N,Y,Y,Y)
OUTPUT(OXID,Y,Y,N,Y,Y,Y)
OUTPUT(FUEL,Y,Y,N,Y,Y,Y)
OUTPUT(H1,Y,Y,N,Y,Y,Y)
OUTPUT(KE,Y,Y,N,Y,Y,Y)
OUTPUT(EP,Y,Y,N,Y,Y,Y)
   *  GROUP 22. Spot-value print-out.
IXMON=16
IYMON=7
   *  GROUP 23. Field print-out and plot control.
ITABL=3
tstswp=-1
dmpstk=t
DISTIL=T
EX(P1)=1.159E+01; EX(U1)=1.252E+01; EX(V1)=2.344E-01
EX(KE)=1.666E+01; EX(EP)=5.617E+01; EX(H1)=7.500E+05
EX(TMP1)=5.000E+02; EX(RHO1)=6.991E-01; EX(PROD)=6.618E-03
EX(OXID)=9.927E-01; EX(FUEL)=6.621E-04; EX(MIXF)=7.280E-03
 LIBREF=133
STOP