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(WINF,WINO,KEINO,EPINO,KEINF,EPINF,DIAM,PRADO,PRADI,CD)
REAL(TEMFU,TEMOX);CHAR(CPDF)
DIAM=0.00762
PRADI=0.5*DIAM
PARAB=T
IF(PARAB) THEN
+ PRADO=2.0*PRADI
ELSE
+ PRADO=5.*PRADI
ENDIF
CD=0.1643;WINF=151.0;WINO=15.1
KEINO=(0.05*WINO)**2; EPINO=0.09*KEINO**2/3.E-5
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
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
+ IREGY=1;GRDPWR(Y,NYF,PRADI,1.0)
+ IREGY=2;GRDPWR(Y,NYO,(PRADO-PRADI),1.3)
+ REAL(DYGDZ);DYGDZ=0.078*2.2;AZYV=0.85;ZWADD=PRADO/DYGDZ
  y/d=({z+zo}/zo)**n   with zwadd=zo & azyv=n
AZYV=0.73;ZWADD=0.035
ELSE
+ NYF=6;NYO=10;NYG=NYF+NYO
+ NREGY=2;NY=16
+ IREGY=1;GRDPWR(Y,NYF,PRADI,1.0)
+ IREGY=2;GRDPWR(Y,NYO,(PRADO-PRADI),1.5)
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
READVDU(CPDF,CHAR,BURN)
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
INIADD=F; FIINIT(W1)=WINO;PRT(F)=0.9
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