```
TEXT(Kent-Bilger H2 Diffusion Flame
TITLE
DISPLAY
The problem considered is the turbulent diffusion flame
of a round jet of hydrogen in a co-flowing stream of
air as studied experimentally by Kent and Bilger [14th Int.
Symposium on Combustion, p1353, 1975]. The jet diameter is
7.62mm and the jet to air stream velocity ratio is 10. For
convenience. The calculation may carried out with the elliptic
solver in a domain 15 diameters downstream, or with the
parabolic solver for 160 diameters downstream.
ENDDIS
GROUP 1. Run title and other preliminaries
REAL(TEMFU,TEMOX);CHAR(CPDF)
DIAM=0.00762
PARAB=T
IF(PARAB) THEN
ELSE
ENDIF
CD=0.1643;WINF=151.0;WINO=15.1
KEINO=(0.05*WINO)**2; EPINO=0.09*KEINO**2/3.E-5
TEMFU=15.+273.;TEMOX=15.+273.
REAL(YFU1IN,YOX1IN,YN21IN,YH21I,YCO1I,YH2O1I,YCO21I)
REAL(YFU2IN,YOX2IN,YN22IN,YH22I,YCO2I,YH2O2I,YCO22I)
** fuel stream
YH21I=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)
IF(PARAB) THEN
+ NYF=25;NYO=15;NYG=NYF+NYO
+ NREGY=2;NY=40
ELSE
+ NYF=6;NYO=10;NYG=NYF+NYO
+ NREGY=2;NY=16
ENDIF
GROUP 5. Z-direction grid specification
IF(PARAB) THEN
+ NZ=550;AZDZ=GRND2;DZW1=0.05
ELSE
+ NZ=16;GRDPWR(Z,NZ,15.*DIAM,1.5)
ENDIF
GROUP 7. Variables stored, solved & named
SOLVE(P1,V1,W1);STORE(VIST,DEN1,TMP1,H1,MH2,MO2,MH2O)
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 8. Terms (in differential equations) & devices
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,H2,CO,H2O,CO2,N2)
** Define fuel & oxidiser composition & temperatures
SCRS(FUIN,0.0,0.0,YH21I,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,SCH2)
MESG( Enter required combustion model: default mixed-is-burnt
MESG( The options are:
MESG(  BURN  - Mixed-is-burnt infinite-rate model
MESG(  DDEL  - Infinite-rate model with Double-Delta PDF
CASE :CPDF: OF
WHEN BURN,4
+ MESG(Infinite-rate model
WHEN DDEL,4
+ MESG(Infinite-rate model with Double-Delta PDF
+ SCRS(PDF,DDELTA)
ENDCASE
MESG(2H2  +  O2 > 2H2O
*** 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
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,KE,ONLYMS,KEINO);COVAL(FREEB,EP,ONLYMS,EPINO)
COVAL(FREEB,P1,1.E3,0.0)
** Exit Boundary Conditions
IF(PARAB) THEN
+ IPARAB=1
ELSE
+ 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)
ENDIF
GROUP 15. Termination of sweeps
IF(PARAB) THEN
+ LITHYD=40
ELSE
+ LSWEEP=200
ENDIF
GROUP 16. Termination of iterations
GROUP 17. Under-relaxation devices
KELIN=1
REAL(RLXFAC); RLXFAC=8.*ZWLAST/WINF/NZ
IF(PARAB) THEN
+ DWDY=T; RELAX(F,LINRLX,0.4); RELAX(DEN1,LINRLX,0.2)
+ RELAX(V1,FALSDT,1.E-2); RELAX(W1,FALSDT,1.E-2)
+ RELAX(KE,LINRLX,0.4); RELAX(EP,LINRLX,0.4)
ELSE
+ 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)
ENDIF
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
OUTPUT(CH4,N,N,N,N,N,N); OUTPUT(CO,N,N,N,N,N,N)
OUTPUT(CO2,N,N,N,N,N,N)
GROUP 22. Spot-value print-out
TSTSWP=-1
IF(PARAB) THEN
+ IYMON=NYF;IZMON=1;ITABL=2;NZPRIN=100;NPLT=5;IPLTL=LITHYD
+ IDISPA=10
ELSE
+ IYMON=NYF+2;IZMON=NZ-1;NPLT=10;ITABL=3;NZPRIN=1
ENDIF
GROUP 23. Field print-out and plot control
GROUP 24. Dumps for restarts
```