TEXT(CO Diffusion Flame
TITLE
  DISPLAY
  The problem considered is the turbulent diffusion flame
  of a round jet of carbon monoxide in a co-flowing stream of
  air. The jet diameter is 7.62mm and the jet to air stream
  velocity ratio is 10. For convenience, the calculation is
  carried out with the elliptic solver in a domain 15
  diameters downstream.
  ENDDIS
    GROUP 1. Run title and other preliminaries
REAL(WINF,WINO,KEINO,EPINO,KEINF,EPINF,DIAM,PRADO,PRADI,CD)
REAL(TEMFU,TEMOX)
DIAM=0.00762
PRADI=0.5*DIAM;PRADO=5.*PRADI;CD=0.1643;WINF=151.0;WINO=15.1
KEINO=(0.05*WINO)**2; EPINO=CD*KEINO**1.5/(0.1*(PRADO-PRADI))
KEINF=(0.05*WINF)**2; EPINF=CD*KEINF**1.5/(0.1*PRADI)
TEMFU=15.+273.;TEMOX=15.+273.
REAL(YFU1IN,YOX1IN,YN21IN,YH21I,YCO1I,YH2O1I,YCO21I)
REAL(YFU2IN,YOX2IN,YN22IN,YH22I,YCO2I,YH2O2I,YCO22I)
  ** fuel stream
YCO1I=1.0
  ** oxidiser stream
YOX2IN=0.232;YN22IN=1.-.232
    GROUP 3. X-direction grid specification
CARTES=F;XULAST=0.1
    GROUP 4. Y-direction grid specification
INTEGER(NYF,NYO,NYG);NYF=4;NYO=8;NYG=NYF+NYO
NREGY=2;NY=12
IREGY=1;GRDPWR(Y,NYF,PRADI,1.0)
IREGY=2;GRDPWR(Y,NYO,(PRADO-PRADI),1.5)
    GROUP 5. Z-direction grid specification
NZ=12;GRDPWR(Z,NZ,15.*DIAM,1.5)
    GROUP 7. Variables stored, solved & named
SOLVE(P1,V1,W1);STORE(VIST,DEN1,TMP1,H1,YSUM)
SOLUTN(P1,P,P,Y,P,P,P);SOLUTN(V1,P,P,P,P,P,N)
SOLUTN(W1,P,P,P,P,P,N);TURMOD(KEMODL)
    GROUP 9. Properties of the medium (or media)
ENUL=3.E-5
  *** START OF EXTENDED SCRS MODEL SETTINGS
PRESS0=1.0000E+05
INTEGER(NSPEC,NELEM);NSPEC=7;NELEM=4
INTEGER(NCSTEP,NCREAC);NCSTEP=-1;NCREAC=1
SCRS(SYSTEM,NCSTEP,NCREAC,NELEM,FASTC)
SCRS(SPECIES,CH4,O2,CO,H2,CO2,H2O,N2)
   ** Define fuel & oxidiser composition & temperatures
SCRS(FUIN,0.0,0.0,YCO1I,0.0,0.0,0.0,YN21IN,TEMFU)
SCRS(OXIN,0.0,YOX2IN,0.0,0.0,0.0,0.0,YN22IN,TEMOX)
SCRS(PROP,CHEMKIN,SCCO)
MESG(Infinite-rate model   2CO  +  O2 > 2CO2
  *** END OF EXTENDED SCRS MODEL SETTINGS
    GROUP 11. Initialization of variable or porosity fields
INIADD=F; FIINIT(W1)=WINO
PATCH(INIT,INIVAL,1,NX,1,NYF,1,NZ,1,LSTEP)
INIT(INIT,W1,0.0,WINF); FIINIT(KE)=KEINF; FIINIT(EP)=EPINF
    GROUP 13. Boundary conditions and special sources
  ** Fuel Stream Inlet Conditions
REAL(YIN);INTEGER(KK1)
INLET(SCRSF,LOW,1,NX,1,NYF,1,1,1,LSTEP)
VALUE(SCRSF,P1,GRND1); VALUE(SCRSF,W1,WINF)
VALUE(SCRSF,F,1.0); VALUE(SCRSF,KE,KEINF)
VALUE(SCRSF,EP,EPINF)
  ** Oxidiser Stream Inlet Conditions
INLET(SCRSO,LOW,1,NX,NYF+1,NYG,1,1,1,LSTEP)
VALUE(SCRSO,P1,GRND1); VALUE(SCRSO,W1,WINO)
VALUE(SCRSO,F,0.0); VALUE(SCRSO,KE,KEINO)
VALUE(SCRSO,EP,EPINO)
  ** Free Boundary Conditions
PATCH(FREEB,NORTH,1,NX,NYG,NYG,1,NZ,1,LSTEP)
COVAL(FREEB,W1,ONLYMS,WINO);COVAL(FREEB,F,ONLYMS,0.0)
COVAL(FREEB,P1,1.E3,0.0)
  ** Exit Boundary Conditions
OUTLET(OUT,HIGH,1,NX,1,NYG,NZ,NZ,1,LSTEP)
COVAL(OUT,P1,1.E3,0.0)
VALUE(OUT,V1,0.0); VALUE(OUT,W1,0.0); VALUE(OUT,F,0.0)
    GROUP 15. Termination of sweeps
LSWEEP=250
    GROUP 16. Termination of iterations
    GROUP 17. Under-relaxation devices
KELIN=1.0
REAL(RLXFAC); RLXFAC=5.*ZWLAST/WINO/NZ
RELAX(V1,FALSDT,RLXFAC); RELAX(W1,FALSDT,RLXFAC)
RELAX(KE,LINRLX,0.4); RELAX(EP,LINRLX,0.4)
RELAX(DEN1,LINRLX,0.3); RELAX(F,LINRLX,0.5)
    GROUP 18. Limits on variables or increments to them
VARMIN(TMP1)=1.E-10; VARMIN(DEN1)=1.E-10
OUTPUT(TMP1,P,P,P,P,Y,Y); OUTPUT(DEN1,P,P,P,P,Y,Y)
    GROUP 20. Preliminary print-out
ECHO=T
    GROUP 21. Print-out of variables
NYPRIN=1;NZPRIN=1
    GROUP 22. Spot-value print-out
TSTSWP=-1;IYMON=NYF+2;IZMON=NZ-1;ITABL=3;NPLT=10
    GROUP 23. Field print-out and plot control
    GROUP 24. Dumps for restarts