```
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(TEMFU,TEMOX)
DIAM=0.00762
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
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
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
```