talk=t;run(1,1)
PHOTON USE
p
*1 3 0.5
up z
con temp x 1 y 1 m z 7 14 fil;.001
con temp x 1 y 1 m z 15 m sh;int 10
vec x 1 y 1 m z 15 m sh
vec x 1 y 1 m z 7 14
set vec ref 4.
vec x 1 y 1 m z 1 6 sh
gr x 1 y 1 m z 7 14
gr ou z 1
gr ou y 1 z 1 25
gr ou y m z 1 25
gr ou z 26 y 16 m
gr ou z 26 y 1 5
gr ou z m y 6 15
gr ou y 16 z 26 m
gr ou y 6 z 26 m
msg Velocity vectors and
msg temperature contours
msg
msg
msg Press
msg to continue
pause;vec off;con off;red
con temp x 1 fil;.001
gr x 1 y 1 m z 7 14
gr ou z 1
gr ou y 1 z 1 25
gr ou y m z 1 25
gr ou z 26 y 16 m
gr ou z 26 y 1 5
gr ou z m y 6 15
gr ou y 16 z 26 m
gr ou y 6 z 26 m
msg Temperature
msg
msg
msg Press
msg to continue
pause;con off;red
con YO2 x 1 fil;.001
gr x 1 y 1 m z 7 14
gr ou z 1
gr ou y 1 z 1 25
gr ou y m z 1 25
gr ou z 26 y 16 m
gr ou z 26 y 1 5
gr ou z m y 6 15
gr ou y 16 z 26 m
gr ou y 6 z 26 m
msg Oxygen
msg
msg
msg Press
msg to continue
pause;con off;red
con YN2 x 1 fil;.001
gr x 1 y 1 m z 7 14
gr ou z 1
gr ou y 1 z 1 25
gr ou y m z 1 25
gr ou z 26 y 16 m
gr ou z 26 y 1 5
gr ou z m y 6 15
gr ou y 16 z 26 m
gr ou y 6 z 26 m
msg Nitrogen
msg
msg
msg Press
msg to continue
pause;con off;red
con YH2 x 1 fil;.001
gr x 1 y 1 m z 7 14
gr ou z 1
gr ou y 1 z 1 25
gr ou y m z 1 25
gr ou z 26 y 16 m
gr ou z 26 y 1 5
gr ou z m y 6 15
gr ou y 16 z 26 m
gr ou y 6 z 26 m
msg Hydrogen
msg
msg
msg Press
msg to continue
pause;con off;red
con YH2O x 1 fil;.001
gr x 1 y 1 m z 7 14
gr ou z 1
gr ou y 1 z 1 25
gr ou y m z 1 25
gr ou z 26 y 16 m
gr ou z 26 y 1 5
gr ou z m y 6 15
gr ou y 16 z 26 m
gr ou y 6 z 26 m
msg Water vapour
msg
msg
msg Press
msg to continue
pause;con off;red
con YCO2 x 1 fil;.001
gr x 1 y 1 m z 7 14
gr ou z 1
gr ou y 1 z 1 25
gr ou y m z 1 25
gr ou z 26 y 16 m
gr ou z 26 y 1 5
gr ou z m y 6 15
gr ou y 16 z 26 m
gr ou y 6 z 26 m
msg Carbon dioxide
msg
msg
msg Press
msg to continue
pause;con off;red
con YCO x 1 fil;.001
gr x 1 y 1 m z 7 14
gr ou z 1
gr ou y 1 z 1 25
gr ou y m z 1 25
gr ou z 26 y 16 m
gr ou z 26 y 1 5
gr ou z m y 6 15
gr ou y 16 z 26 m
gr ou y 6 z 26 m
msg Carbon monoxide
msg
msg
msg Press
msg to continue
enduse
GROUP 1. Run title and other preliminaries
TEXT(System-level coke-combustion simulation
DISPLAY
System level lump coke packed bad combustion.
One-phase, 2D, one space, lump coke packed bed
combustion, air supply underneath and above bed.
ENDDIS
REAL(HINCL,CINCL,NINCL,AINCL,GASCON)
REAL(AIRO2,AIRN2)
REAL(MN2,MC,MO2,MH2,MCO,MCO2,MH2O)
GASCON=8.3143e3
HINCL=0.05
CINCL=0.95
NINCL=0.0
AINCL=1.-CINCL-HINCL-NINCL
AIRO2=0.232
AIRN2=0.768
MN2=28.; MC=12.; MO2=32.; MH2=2.; MCO=28.; MCO2=44.; MH2O=18.
REAL(HCCO2,HCCO,HHH2O,HCHX,HCOCO2)
HCCO2 =3.279E7
HCCO = 9.208E6
HHH2O = 1.209E6
REAL(FS,BURNRATE, WIN)
** FS is the mass of fuel per unit mass of air/fuel mixture
to convert all carbon and oxygen to carbon monoxide.
FS=0.232/(0.232 + CINCL*16.0/12.0)
** The heat of combustion per unit mass of co is hcco2 minus hcco
the mass of c per unit mass of co, ie 12/28
HCOCO2=(12.0/28.0)*(HCCO2-HCCO)
** The heat of coal combustion per unit mass of carbon
HCHX=(CINCL*HCCO2+HINCL*HHH2O)*(HINCL*MH2+CINCL*MC+NINCL*MN2)/MC
** The rate of burning
BURNRATE=9.
** The inlet gas velocity
WIN=2.
GROUP 4. Y-direction grid specification
NY=20;GRDPWR(Y,NY,1.,1.0)
GROUP 5. Z-direction grid specification
NZ=30;GRDPWR(Z,NZ,1.5,1.0)
GROUP 7. Variables stored, solved & named
** FCL - carbon element mass fraction;
** The specific heat at constant pressure, CP =1.100E3
H = CP*T + HCHX*YCHX + HCOCO2*YCO * HHH2*YH2
SOLVE(P1,V1,W1,H1,FCL)
REAL(CP,TFUEL,HGIN,TGIN)
TGIN = 350.0
CP= 1100.
HGIN = CP*TGIN
STORE(RMIX,HSUB,TEMP,YN2,YH2,YO2,YCO,YCO2,YH2O)
STORE(DEN1,FLIM,FRAC,GO,GC,GH,GOFU,GOPA)
SOLUTN(P1 ,Y,Y,Y,P,P,P)
SOLUTN(FCL,Y,Y,Y,P,P,P)
SOLUTN(H1 ,Y,Y,Y,P,P,P)
GROUP 8. Terms (in differential equations) & devices
TERMS(FCL,N,Y,N,P,P,P)
TERMS(H1 ,N,Y,N,P,P,P)
GROUP 9. Properties of the medium (or media)
REAL(RHOIN1,WAIR)
PRESS0=1.e5
WAIR=32.
RHOIN1=PRESS0*WAIR/(8314.*TGIN)
RHO1=GRND
DEN1=PRESS0/(RMIX*TEMP+tiny)
DEN1=AMIN1(VARMAX(140),AMAX1(0.0,DEN1,VARMIN(140)))
ENUL=0.0
GROUP 11. Initialization of variable or porosity fields
FIINIT(W1)=1.
STORE(VPOR,NPOR,HPOR)
FIINIT(FCL)=FS
FIINIT(P1)=1.3E-4
INIADD=F
real(porbed);porbed=0.75
PATCH(POROS,INIVAL,1,NX,1,NY,7,14,1,1)
COVAL(POROS,VPOR,0.0,PORBED)
COVAL(POROS,NPOR,0.0,PORBED)
COVAL(POROS,HPOR,0.0,PORBED)
PATCH(SCHIMNEY,INIVAL,1,NX,1,5,26,NZ,1,1)
COVAL(SCHIMNEY,VPOR,0.0,0.)
PATCH(NCHIMNEY,INIVAL,1,NX,16,NY,26,NZ,1,1)
COVAL(NCHIMNEY,VPOR,0.0,0.0)
PATCH(POR1,INIVAL,1,NX,1,1,6,6,1,1)
COVAL(POR1,hpor,0.0,0.0)
PATCH(POR2,INIVAL,1,NX,4,5,6,6,1,1)
COVAL(POR2,hpor,0.0,0.0)
PATCH(POR3,INIVAL,1,NX,10,10,6,6,1,1)
COVAL(POR3,hpor,0.0,0.0)
PATCH(POR4,INIVAL,1,NX,11,11,6,6,1,1)
COVAL(POR4,hpor,0.0,0.0)
PATCH(POR5,INIVAL,1,NX,14,15,6,6,1,1)
COVAL(POR5,hpor,0.0,0.0)
PATCH(POR6,INIVAL,1,NX,20,20,6,6,1,1)
COVAL(POR6,hpor,0.0,0.0)
GROUP 13. Boundary conditions and special sources
** Air to the fireplace bottom
INLET(INLET,EAST,1,1,1,NY,1,6,1,1)
VALUE(INLET,P1 , RHOIN1*WIN)
VALUE(INLET,W1 , WIN)
VALUE(INLET,FCL , 0.0)
VALUE(INLET,H1 , HGIN)
** Frictional momentum transfer
PATCH(BEDFRIC,VOLUME,1,NX,1,NY,7,14,1,1)
COVAL(BEDFRIC,W1,100.,0.0)
COVAL(BEDFRIC,v1,100.,0.0)
** Chimney outlet
PATCH(OUTLET,HIGH,1,NX,1,NY,NZ,NZ,1,1)
COVAL(OUTLET,P1,fixp,0.0)
** Air to the open flame area
PATCH(EASTAIR,EAST,1,NX,1,NY,15,25,1,1)
COVAL(EASTAIR,P1 , FIXFLU , 2.*RHOIN1*WIN)
COVAL(EASTAIR,FCL , ONLYMS, 0.0)
COVAL(EASTAIR,H1 , ONLYMS, HGIN)
**
Carbon mass transfer related sources:
------------------------------------
PATCH(CARGAS,VOLUME,1,NX,1,NY,1,NZ,1,1)
(1) Transfer of mass leading to increase of gas flow rate:
- VPOR is volume fraction of lump coal
VAL=:BURNRATE:*(1.-VPOR)*(:FS:-FCL)
COVAL(CARGAS,P1,FIXFLU,GRND)
(2) Transfer of carbon leading to increase of mixture
fraction at the same rate:
- CO=1. signifies that mass tarnsfer brings in
material which is 100% carbon
COVAL(CARGAS,FCL,ONLYMS,1.)
(3) Transfer of enthalpy and heat leading to increase of
gas enthalpy at the same rate:
- Interphase gas temperature is assumed as TEMP.
- HSUB = HCOCO2*YCO * HH2*YH2
VAL=:CP:*TEMP+:HCHX:+HSUB
COVAL(CARGAS,H1,ONLYMS,GRND)
GROUP 15. Termination of sweeps
LSWEEP=100
RESFAC=0.0001
GROUP 16. Termination of iterations
LITHYD=10
VARMAX(FCL)=1.0;VARMIN(FCL)=0.0
VARMIN(TEMP)=TGIN;VARMAX(TEMP)=3000.
VARMIN(DEN1)=0.001;VARMAX(DEN1)=3.
GROUP 17. Under-relaxation devices
RELAX(P1,LINRLX,0.175)
RELAX(W1,FALSDT,0.1)
RELAX(V1,FALSDT,0.1)
RELAX(FCL,FALSDT,0.1)
RELAX(DEN1,LINRLX,0.1)
GROUP 22. Monitor print-out
IZMON=NZ/2;IYMON=NY/2;UWATCH=T
GROUP 23. Field print-out and plot control
NPLT=1;NYPRIN=1;NZPRIN=1
NYPRIN=1;IYPRF=1;IYPRL=30
TSTSWP=-1
namsat=mosg
** Cell-wise composition parameters
--------------------------------
FLIM=:AIRO2:/(:AIRO2:+:CINCL:*:MO2:/:MC:+$
:HINCL:*:MO2:/(2*:MH2:))
GO=:AIRO2:*(1-FCL)
GC=:CINCL:*FCL
GH=:HINCL:*FCL
GOPA=GC*:MO2:/(2*:MC:)/(1-GO+GC*:MO2:/(2*:MC:)+TINY)
GOFU=(GH*:MO2:/(2*:MH2:)+GC*:MO2:/:MC:)/$
(1.-GO+GH*:MO2:/(2*:MH2:)+GC*:MO2:/:MC:+TINY)
FRAC=(GO-GOPA)/(GOFU-GOPA+TINY)
** For all regions
---------------
YN2=:NINCL:*FCL+:AIRN2:*(1.-FCL)
** Region 1
--------
YH2O=:HINCL:*FCL*:MH2O:/:MH2:
IF(FCL.LE.FLIM)
YCO2=:CINCL:*FCL*:MCO2:/:MC:
IF(FCL.LE.FLIM)
YO2 =:AIRO2:*(1-FCL)-:CINCL:*FCL*:MO2:/:MC:-$
:HINCL:*FCL*:MO2:/(2.*:MH2:)
IF(FCL.LE.FLIM)
YCO=0.0
IF(FCL.LE.FLIM)
YH2=0.0
IF(FCL.LE.FLIM)
HSUB=0.0
IF(FCL.LE.FLIM)
RMIX=:GASCON:*(YO2/:MO2:+YH2O/:MH2O:+YCO2/:MCO2:+$
YN2/:MN2:)
IF(FCL.LE.FLIM)
** Region 2
--------
YH2O=:HINCL:*FCL*:MH2O:/:MH2:*FRAC*(1-GOFU)/(1-GO+TINY)
IF(FCL.GT.FLIM.AND.FRAC.GE.0.)
YCO2=:CINCL:*FCL*:MCO2:/:MC:*FRAC*(1-GOFU)/(1-GO+TINY)
IF(FCL.GT.FLIM.AND.FRAC.GE.0.)
YO2=0.0
IF(FCL.GT.FLIM.AND.FRAC.GE.0.)
YCO=:CINCL:*FCL*:MCO:/:MC:*(1-FRAC)*$
(1-GOPA)/(1-GO+TINY)
IF(FCL.GT.FLIM.AND.FRAC.GE.0.)
YH2=:HINCL:*FCL*(1-FRAC)*(1-GOPA)/(1-GO+TINY)
IF(FCL.GT.FLIM.AND.FRAC.GE.0.)
YH2O=AMAX1(0.,YH2O)
YCO=AMAX1(0.,YCO)
YCO2=AMAX1(0.,YCO2)
HSUB=YCO*:HCOCO2:+YH2*:HHH2O:
IF(FCL.GT.FLIM.AND.FRAC.GE.0.)
RMIX=:GASCON:*(YH2O/:MH2O:+YCO/:MCO:+YCO2/:MCO2:+$
YH2/:MH2:+YN2/:MN2:)
IF(FCL.GT.FLIM.AND.FRAC.GE.0.)
** Region 3
--------
YH2O=0.0
IF(FCL.GT.FLIM.AND.FRAC.LT.0.)
YCO2=0.0
IF(FCL.GT.FLIM.AND.FRAC.LT.0.)
YO2=0.0
IF(FCL.GT.FLIM.AND.FRAC.LT.0.)
YCO=:AIRO2:*(1-FCL)*2*:MCO:/:MO2:
IF(FCL.GT.FLIM.AND.FRAC.LT.0.)
YH2=:HINCL:*FCL
IF(FCL.GT.FLIM.AND.FRAC.LT.0.)
HSUB=YCO*:HCOCO2:+YH2*:HHH2O:
IF(FCL.GT.FLIM.AND.FRAC.LT.0.)
RMIX=:GASCON:*(YCO/:MCO:+YH2/:MH2:+YN2/:MN2:)
IF(FCL.GT.FLIM.AND.FRAC.LT.0.)
** Calculation of absolute gas temperature
--------------------------------------
TEMP=(H1-HSUB)/:CP:
TEMP=AMIN1(VARMAX(147),AMAX1(100.,TEMP,VARMIN(147)))
store(ysum)
YSUM=YN2+YO2+YCO+YCO2+YH2O+YH2
STOP